Конечно-разностный метод численного моделирования распространения сейсмических волн в трехмерно-неоднородных разномасштабных средах

Для численного моделирования процессов распространения сейсмических волн в трехмерно-неоднородных средах с разномасштабными неоднородностями (кавернозно-трещиноватыми резервуарами) разработан конечно-разностный метод, основанный на использовании сеток с локальным пространственно-временным измельчением. Необходимость использования таких сеток связана с огромными различиями в масштабах неоднородностей вмещающей среды (десятки и сотни метров) и микроструктуры пласта-коллектора (от долей сантиметра до первых метров). Вычисления и на грубой, и на мелкой сетках требуют использования высокопроизводительных вычислительных систем с параллельной архитектурой. Для реализации параллельных вычислений используется трехмерная декомпозиция области, когда каждый элементарный объем приписывается своему процессорному элементу. Эти процессорные элементы объединяются в две группы – для вмещающей среды (крупная сетка) и резервуара (мелкая сетка). Вычисления внутри каждой из групп производятся асинхронно за счет использования неблокирующих процедур iSend/iRecv библиотеки MPI. Обмены между группами организованы также с помощью этих процедур и осуществляются через специально выделенные мастер-процессоры. Приводятся результаты численных расчетов для реалистичных моделей карбонатных резервуаров, содержащих коридоры трещиноватости. Показано проявление ориентации этих коридоров в рассеянных волновых полях.