Локальное пространственно-временное измельчение сеток для конечно-разностного моделирования упругих волн в трёхмерно-неоднородных разномасштабных средах
В.И. Костин1, В.В. Лисица2, Г.В. Решетова3, В.А. Чеверда2
1ЗАО «Intel A/O», просп. Акад. М.А. Лаврентьева, 6/1, Новосибирск, 630090
victor.i.kostin@intel.com
2Институт нефтегазовой геологии и геофизики Сибирского отделения Российской академии наук, просп. Акад. В.А. Коптюга, 3, Новосибирск, 630090
lisitsavv@ipgg.sbras.ru
3Институт вычислительной математики и математической геофизики Сибирского отделения Российской академии наук, просп. Акад. М.А. Лаврентьева, 6, Новосибирск, 630090
kgv@nmsf.sscc.ru
Ключевые слова: сейсмические волны, конечно-разностные методы, декомпозиция области, интерполяция, группы процессорных элементов
Страницы: 45-55
Аннотация
Для численного моделирования процессов распространения сейсмических волн в трёхмерно-неоднородных средах с разномасштабными неоднородностями разработан конечно-разностный метод, использующий сетки с локальным пространственно-временным измельчением. Для описания вмещающей среды используется сравнительно грубая сетка (20÷30 точек на длину волны), в то время как мелкомасштабные неоднородности описываются на гораздо более детальной сетке (от 200 точек на длину волны). Все вычисления производятся с использованием компьютеров с параллельной архитектурой и основаны на трёхмерной декомпозиции расчётной области, когда каждый элементарный объём приписывается своему процессорному элементу. Используемые процессорные элементы формируют две группы — для вмещающей среды (крупная сетка) и скоплений мелкомасштабных неоднородностей (мелкая сетка). Взаимодействие групп организовано специальным образом, опирается на использование специальным образом выделенных мастер-процессоров и предусматривает согласование сеток. Приводятся результаты численных расчётов для реалистичных моделей карбонатных резервуаров, содержащих коридоры трещиноватости. Показано проявление ориентации этих коридоров в рассеянных волновых полях.
|
