Численное моделирование распространения сейсмических волн в средах с вязкоупругими включениями
Вишневский Д.М., Лисица В.В., Решетова Г.В.

Поглощение сейсмической энергии является одним из наиболее важных физических явлений, вызванным различными геологическими факторами и требующим детального изучения его проявлений в регистрируемых волновых полях. В силу ряда обстоятельств сколько-нибудь заметное поглощение сосредоточено в сравнительно небольших подобластях изучаемой геологической среды, занимающих, как правило, порядка 10–20% общего объема. В то же время, широко используемые в настоящее время подходы для численного моделирования распространения сейсмических волн в таких средах требуют значительных вычислительных ресурсов, заметно превышающих нужные для идеально упругих сред. По этой причине их использование во всей расчетной области является неэффективным. В настоящей статье изложен подход к численному моделированию упругих волн, основанный на использовании уравнений вязкоупругости только там, где это действительно необходимо, что позволяет существенно сократить потребности в вычислительных ресурсах. Для этого расчетная область представляется в виде суперпозиции подобластей, заполненных вязко- и идеально упругой средами, в каждой из которых используется своя математическая модель и своя конечно-разностная схема. Особое внимание уделяется взаимосогласованию этих двух моделей и соответствующих им конечно-разностных схем и шаблонов, в том числе и для обеспечения приемлемого уровня артефактов при переходе от одного шаблона к другому. Для организации параллельных вычислений выполняется дополнительная декомпозиция каждой из этих подобластей, обеспечивающая возможность ее размещения на отдельном узле вычислительного кластера с организацией их взаимодействия посредством библиотеки MPI. Показано, что применение такого подхода комбинированной декомпозиции расчетной области снижает время вычислений примерно в 1.7 раза по сравнению с использованием вычислительных схем, ориентированных исключительно на вязкоупругие среды.

Ключевые слова: теория упругости, среды с поглощением, конечно-разностные схемы, расщепление по пространственным подобластям, параллельные алгоритмы

Название статьи, аннотация и ключевые слова на английском языке

Вишневский Д.М., науч. сотр., e-mail: vishnevskydm@ipgg.sbras.ru;   Лисица В.В., зав. лабораторией, e-mail: lisitsavv@ipgg.sbras.ru; – Институт нефтегазовой геологии и геофизики СО РАН, просп. Коптюга, 3, 630090, Новосибирск;
Решетова Г.В., ст. науч. сотр., e-mail: kgv@sscc.nsc.ru – Институт вычислительной математики и математической геофизики СО РАН (ИВМиМГ СО РАН), просп. Лаврентьева, 6, 630090, г. Новосибирск