Выпуск: 4 , Том: 14 , Год издания: 2016

Предлагается новый численный алгоритм моделирования распространения сейсмических волновых полей. Конечно-разностное решение уравнений упругости происходит только в полосе, следующей за фронтом первых вступлений, которую далее будем называть "окном". Положение "окна", т. е. фронта первых вступлений, рассчитывается путем численного решения уравнения эйконала. Предлагаемый подход "оконного" моделирования позволяет существенно ускорить вычисления (так как конечно-разностные вычисления проводятся в бегущем "окне", а не во всей расчетной области) и уменьшить объем оперативной памяти, необходимый для хранения поля. На тестах показана возможность ускорения до 40 раз и экономии памяти до 50 раз. Данный подход позволяет моделировать только волновые формы первых вступлений, но этого оказывается достаточно для многих процедур сейсмической миграции. В частности, в статье рассматривается распространенная процедура миграции в обращенном времени, реализация которой требует больших объемов вычислений.