Построение изображения сейсмического разреза по модели BP2004 Benchmark

Задачи сейморазведки всегда требовали серьёзных вычислительных мощностей. В данной работе рассматривается метод миграции в обратном времени (Reverse time migration — RTM). Автором проделано тестирование программного кода [9] на модели EAGE BP2004 Benchmark. Модель представляет собой двумерный геологический разрез с заданным распределением скоростей (рис. 1). Модель создана по инициативе геофизического общества EAGE и предназначена для сравнительного тестирования программных разработок. Задача построения сейсмических изображений — ресурсоемкая и сложная задача в обработке данных. Одной из проблем миграции сейсмических данных, в том числе и методом RTM, является ограниченность и дискретность наблюдения. В работе рассмотрены проблемы, возникающие при построении изображения методом RTM, а также способы их решения. Для этого по модели EAGE BP2004 выполняется расчёт прямой задачи и полученные модельные данные применяются для тестирования RTM миграции. Отличительной особенностью программного кода является возможность расчёта с использованием неявных дифференциальных операторов (компактных схем). Произведено тестирование точности и времени исполнения алгоритма с использованием разных схем на задаче построения изображения данной модели. В работе также освещены вопросы фильтрации изображений. Рассмотрены различные методы нормализации и высокочастотная фильтрация.