Вычислительные проблемы моделирования природных и индустриальных процессов в Арктической зоне Российской Федерации

1. Вычислительные проблемы моделирования природных и индустриальных процессов в Арктической зоне Российской Федерации

Чл-корр. РАН, д.ф.-м.н., профессор,
заведующий кафедрой информатики и
вычислительной математики МФТИ
Игорь Борисович Петров,
[email protected]

2. План выступления

Численное решение задач освоения Арктики
Решение задач сейсморазведки, в том числе морской,
а также в условиях Арктики
Численное моделирование в геологии
Численное моделирование экспериментов по
исследованию прочностных характеристик льда

3.

Вычислительные задачи Арктики
1.Прямые и обратные задачи сейсморазведки. Задачи миграции и
инверсии.
2.Прямые и обратные задачи электроразведки;
3.Комплексные задачи сеймо- и электроразведки;
4.Расчет механических и прочностных характеристик льда, как
твердого деформируемого тела;
5.Задачи миграции (дрейфа) крупных ледовых образований (КЛО);
6.Задачи безопасности морских стационарных ледостойких
платформ (МСЛСП) и плавающих (якорных) платформ
нефтегазовой промышленности при воздействии природных
факторов (ЛО);
7.Задачи заторашивания МСЛП;
8.Моделирование процессов образования КЛО;
9.Расчет на прочность нефтегазо-проводов, наземных и поддонных;
10.Проблемы безопасного плавания судов ледового класса при
наличии КЛО;
11.Посадка самолета на плавающую льдину (ледовое поле).

4.

12.Расчет на прочность ледового покрытия при движении по нему
транспортных средств.
13.Воздействия периодических и сейсмических нагрузок на технические
сооружения в Арктике.
14.Климатические задачи Арктики.
15.Прогнозирование динамики ледовой обстановки с учетом обработки
данных наблюдений в воздушном и водном бассейнах Арктики.
16.Обнаружение “метановых бомб” на Ямале.
17.Задачи геомеханики (расчет локализации контактных поверхностей в
геологических средах).
18.Определение положения газоносных слоев вблизи скважины (задачи
безопасности скважин).
19.Расчет поля температур в КЛО с учетом фазовых переходов,
солнечной радиации и разрушения льда при их движении в
Северных морях.
20.Расчет жизненного цикла ледового острова с учетом тепловых,
радиационных и механических нагрузок.
21.Проблема связи в Арктической зоне РФ.
22.Проблема жизнеобеспечения в условиях низких температур.

5.

На пути численного моделирования данного класса
задач стоит ряд не до конца решенных проблем:
масштабный эффект: обобщение результатов лабораторных экспериментов
(измеренные механические, прочностные и др. параметры) на натурные
полномасштабные эксперименты;
неоднородность параметров льда и толщины ледяного поля: параметры
неоднородности зависят от природных условий, требуют сбора натурных
данных и статистической обработки; из-за неполноты исходных данных в
некоторых случаях используются вероятностные подходы;
сложная кристаллическая структура льда: морской лёд анизотропный
(трансверсально-изотропный) ;
морской лёд слоистый: у каждого слоя свои механические и прочностные
параметры, которые зависят от температуры, солёности, пористости;
существование скоростных режимов: на малых скоростях в зоне контакта
наблюдается перекристаллизация льда, образование микротрещин из-за
сдвиговых напряжений и таяние граничного слоя, что в совокупности называют
пластическим режимом; режим поведения на высоких скоростях называют
хрупким.

6.

одновременное наличие принципиально разных с точки зрения
численного моделирования режимов разрушения: трещины, характерные
для твердого тела, а также смесь из раскрошенного льда и более крупных
осколков; это и определяет «пилообразный» вид нагрузки на опору;
необходимость значительного измельчения расчётной сетки в области
контакта в случае тонкого ледяного поля и широкого основания опоры ведёт
к значительному увеличению требуемых вычислительных ресурсов;
ввиду того, что практический интерес представляют как пиковые, так и
средние нагрузки на опору, а максимальные скорости движения ледяных
полей порядка 0.1 м/с, то при использовании явных схем интегрирования по
времени необходимо вычислить миллионы шагов по времени;
эффективный способ распараллеливания алгоритма контактного
взаимодействия большого числа объектов на многопроцессорных
вычислительных системах;
значительные затраты машинного времени из-за того, что характерные
времена процесса существенно превышают акустические (и,
соответственно, курантовский шаг).

7. Численное решение задач освоения Арктики

8. Миграция айсбергов

9. Фотография повреждений корабля R.E. Gagnon, J. Wang Numerical simulations of a tanker collision with a bergy bit incorporating

hydrodynamics, a validated ice model and damage
to the vessel // Cold regions. Science and Technology, 2012.

10.

Столкновение ледокола с торосом

11. Воздействие килей торосов на дно и подводные трубопроводы М.А. Наумов, Д.А. Онищенко, Презентация ООО «Газпром ВНИИГАЗ»

12. Разрушение айсберга под интенсивными динамическими воздействиями.

13. Разрушение айсберга под интенсивными динамическими воздействиями.

14. Торосы и нефтедобывающая платформа

15. Торосы и нефтедобывающая платформа

16. Прорезание стойки платформы ледяным полем

17. Прорезание стойки платформы ледяным полем

18. Прорезание стойки платформы ледяным полем

19. Столкновение айсберга со стойкой стационарной нефтедобывающей платформы

19

20. Строение торосов A. Marchenko Thermodynamic consolidation and melting of sea ice ridges // Cold regions. Science and

Technology, V. 52, N. 3, 2008.

21. Ледяное тело заданной формы с заданным распределением водонасыщенных и газонасыщенных полостей и трещин

22. Решение задач сейсморазведки, в том числе морской, а также в условиях Арктики

23. Сейсморазведка - стример

• 3D
• P-волны
• Высокая
производительность

24. Сейсмика – донные станции

• 3D/4C
• Высокая
стоимость
• Высокая
информативность
данных

25. Электроразведка – донные станции

Лидер по объёмам
работ
6 компонент ЭМ поля
(важно при 3D
инверсии)
Не мельче 50 м
EMGS, http://www.emgs.no

26. Электроразведка - стримеры

PGS, http://www.pgs.com/
Высокая производительность
Не глубже 300 м
Одна осевая компонента поля: Ex
Частотная и временная область

27. Многослойная порода

28. Геология со сложными границами

29. Геология со сложными границами

30. Геология со сложными границами

31. Сейсмическая разведка Арктического шельфа

32. Волновая картина в слое льда

33. Волновая картина в слое воды

34. Волновая картина в слое морского дна

35. Волновая картина в нефтесодержащем резервуаре

36. Постановки задач

Источник во льду
Источник во льду, без резервуара
Источник на дне
Источник на дне, без резервуара

37. Волновые картины

Источник во льду, 0.135 сек.
Источник во льду, без резервуара, 0.135 сек.
Источник на дне, 0.0675 сек.
Источник в воде, без резервуара, 0.0675 сек.

38. Сейсмограммы, лед, Vy

Источник во льду
Источник во льду, без резервуара
Источник на дне
Источник на дне, без резервуара

39. Сейсмограммы, дно, Vy

Источник во льду
Источник во льду, без резервуара
Источник на дне
Источник на дне, без резервуара

40. Источник на дне

41. Источник на дне, без резервуара

42. Численное моделирование в геологии

43. Численное моделирование в геологии

Трещины и пустоты

44. Трещины и пустоты

Коридор флюидонасыщенных
вертикальных трещин

45. Коридор флюидонасыщенных вертикальных трещин

расстояние между трещинами / длина трещин
0,5
1,0
2,0
3,0
1,5
4,0

46. Коридор флюидонасыщенных вертикальных трещин

Простая флюидонасыщенная полость
Отраженная
продольная
волна
Прошедшая волна
Отраженная
волна

47. Простая флюидонасыщенная полость

Многослойная геологическая среда

48. Многослойная геологическая среда

Численное моделирование
экспериментов по исследованию
прочностных характеристик льда

49. Численное моделирование экспериментов по исследованию прочностных характеристик льда

Мотивация

50. Мотивация

Механико-математическая
модель льда
изотропная упруго-идеально-пластическая модель
критерий хрупкого разрушения по главным
напряжениям и пластического (объемного) разрушения
в случае превышения величины пластической
деформации порогового значения

51. Механико-математическая модель льда

Прочность льда на одноосное
сжатие
Сдавливание ледяного куба с постоянной скоростью,
где измеряется величина приложенной нагрузки от
осевого смещения
На внутренней картине разрушения синим цветом показано разрушение по критерию
на растяжение, красным – объёмное разрушение ячеек

52. Прочность льда на одноосное сжатие

Сечение картины разрушений
ледяного образца

53. Сечение картины разрушений ледяного образца

Картины разрушения
ледяного образца в натурных
экспериментах

54. Картины разрушения ледяного образца в натурных экспериментах

Сравнение компьютерного и
натурного эксперимента
Снаружи преобладают вертикальные трещины,
образовавшиеся в результате локальных растяжений.
Наименее подвержены разрушению конические
области, находящиеся в непосредственном контакте с
прессом.
Внутри образца преобладает объёмное разрушение, в
результате чего лёд крошится в мелкую крошку.

55. Сравнение компьютерного и натурного эксперимента

Картины разрушения
цилиндрических образцов в
натурных экспериментах
Хомуты на концах соответствуют модели закреплённой границы

56. Картины разрушения цилиндрических образцов в натурных экспериментах

Зависимости силы сжатия образца от времени
при варьировании параметров модели
прочности льда.
По оси X отложен номер временного отсчёта (интервал между отсчётами
0.00023c), по оси Y – сила сжатия пресса, Н.

57.

Изменение качественной картины разрушений
при варьировании значений параметра
максимальной пластической деформации
центральный расчёт соответствует типовому значению, слева предел увеличен,
справа - уменьшен

58.

Прочность льда на изгиб 4-х точечным
методом

59.

Зависимость силы нагрузки от времени
на индентер

60.

Картины растягивающих напряжений перед
образованием трещины нормального разрыва и
после него

61.

Внешняя и внутренняя структура
образовавшейся трещины

62.

Пример задачи о деформировании в многосвязной
области интегрирования
Морская ледостойкая стационарная платформа

63.

Моделирование столкновения танкера с причалом
Институт автоматизации проектирования РАН.
Н/рук. темы д.ф.-м.н. В.Л.Якушев
Сжиженный газ (Liquefied Natural Gas - LNG) привлекает внимание как источник
экологически чистой энергии, который может быть доставлен по морю в различные
страны, но вместе с тем представляет серьезную угрозу при катастрофах.
РИС. 1. Конечно-элементная
модель танкера.
РИС. 2. Распределение напряжений при
ударе.

64.

Штормовое воздействие на ледостойкую
платформу (ИВМ РАН)

65. Штормовое воздействие на ледостойкую платформу (ИВМ РАН)

. Схема распределения нагрузки от шасси
самолетов: a) ИЛ-76 ТД, b) C-130H Hercules.
a)
b)

66. . Схема распределения нагрузки от шасси самолетов: a) ИЛ-76 ТД, b) C-130H Hercules.

Движение одного автомобиля по ледяному покрову на поверхности
воды. Параметры эксперимента: масса автомобиля – 2.2т; скорость –
15м/с; толщина льда – 0.26м. (ИПМ РАН)

67. Движение одного автомобиля по ледяному покрову на поверхности воды. Параметры эксперимента: масса автомобиля – 2.2т; скорость –

Концентрация (компактность в %) льда по данным наблюдений
MASAM2 в Арктической зоне (a) и в выделенном Арктическом регионе
России (b). (ИВМ РАН)
a)
b)

68. Концентрация (компактность в %) льда по данным наблюдений MASAM2 в Арктической зоне (a) и в выделенном Арктическом регионе

Спасибо
за внимание!