Упрощенный анализ скоростей временной миграции.
Традиционный анализ скоростей миграции – весьма трудоемкий процесс.
Автором системы SPS-PC доказано, что при анализе скоростей миграции на прямых лучах не нужно получать набор мигрированных сейсмограмм для заданного перебора скоростей. Достаточно получить одну мигрированную сейсмограмму и получить правильную скорость путем анализа кривизны осей синфазности на этой сейсмограмме. Поправки вычисляются по способу средних скоростей.
Тут X0 – координата трассы-результата (точек дифракции). X1 – координата источника в пределах апертуры, X2 – координата приемника в пределах апертуры.
Пусть V0 – некоторая априорная скорость. Она всегда имеется, например, по скважине.
Кинематическая поправка в точке (X0,Z1) определяется как
(1)
(2)
(3)
(4)
Пусть истинная скорость
(5)
Если - положительная, то ось синфазности на мигрированной сейсмограмме будет переспрямлена
Если - отрицательная, то ось синфазности на мигрированной сейсмограмме будет недоспрямлена
Нетрудно разложить выражение для в ряд Маклорена по .
Получаем вот такой симпатичный ряд.
(6)
где DT(0) – уже введенные в трассы перед суммированием кинематические поправки, рассчитанные по формуле (1). Для мигрированной сейсмограммы эти поправки аппроксимируются формулой
(7)
тут: Ze – эффективная глубина; L –удаление (параметр трассы мигрированной сейсмограммы).
Ряд (6) сходится, если остаточная поправка dV не превышает V0. Ряд сходится быстро, если остаточная поправка dV не превышает V0/2.
Предложенный алгоритм в десятки раз ускоряет процедуру скоростного анализа MVD.
Процесс может быть итерационным. Если второе приближение – недостаточно точное (например, из-за малого количества членов в ряде Маклорена или больших значений dV), то процесс можно повторить, используя в качестве V0 новое полученное решение.
См. также: Вывод аппроксимационной формулы 7 .
Процедура вызывается через пункт главного меню Time Migrations->Simple Velocity Spectrum.
Появляется диалоговое окно Simple Vertical Time Migration Velocity Analysis, в котором надо указать:
Apertute (in CMP points) – база (апертура) суммирования.
Максимальное количество сейсмограмм ОГТ слева и справа от трассы-результата, участвующие в суммирование по алгоритму Кирхгофа. Чем больше это значение – тем лучше результат определения скорости.
Целое число.
Max Slope (degrees) – Угол, ограничивающий апертуру.
Его можно определить как максимальный угол наклона отражающих границ в градусах. Задается целым числом в переделах 1-89.
Целое число.
Moveout interpolation mode – способ интерполяции кинематических поправок.
Первичные поправки вычисляются с шагом Moveout Time Step.
Значения DT во всех точках трассы вычисляются путем сплайн-интерполяции одним из 4-х алгоритмов.
-
Akima's
method (continuous 1st derivatives, only);
-
Fritsch-Carlson
method (continuous 1st derivatives, only);
-
Cubic
spline (continuous 1st and 2nd derivatives)
-
Hermite
Polynomial (continuous 1st derivatives, only).
- Linear Interpolation – простая линейная интерполяция.
Лучше всего подходит интерполяция “Fritsch-Carlson method”. Это такой вариант сплайнов, который наилучшим образом сохраняет монотонность функции.
Maximum Available Moveout Gradient - максимально допустимый градиент кинематических поправок.
GRK=dTau/dt
Tau - Кинематическая поправка. t - время.
Вещественное число.
Пределы от 0.0 до 100.0. Рекомендуется 0.5
При GRK=0.0 растяжение сигнала не допускается,
При GRK>=1.0 допускается "пересечение" годографов дифрагированных волн.
Time Window – величина временного окна в мсек для расчета спектра Semblance.
Целое число.
Base Velocity – начальное приближение значений средней скорости (м/сек).
Целое число. Если 0, то начальные приближения берутся из сетки скоростей в формате GS Binary.
Percent Velocity – приращение значений средней скорости в процентах
в большую и меньшую сторону от начального приближения.
Целое положительное число.
Number of Vels – количество приращений значений средней скорости.
Половина этих приращений будет в меньшую сторону, половина – в большую.
Целое число нечетное число.
Start CMP picket – Пикет ОГТ, с которого начинается скоростной анализ.
Целое число.
Start CMP picket – Пикет ОГТ, на котором заканчивается скоростной анализ.
Целое число.
CMP picket Step – Шаг пикетов ОГТ скоростного анализа.
Целое число.
Plotting Parameters – параметры отрисовки мигрированных разверток и спектров Semblance при последующей обработке результатов анализа. Рекомендуем предварительно подобрать эти параметры анализом одной-двух точек и затем запустить анализ на весь профиль.
Plot Time Scale (mm/sec) – Временной масштаб в мм / сек.
Если ноль, то масштаб определяется размером экрана.
Целое число.
Number of Traces per cm – Масштаб по оси Х. Количество трасс в 1 см.
Если ноль, то масштаб определяется размером экрана.
Целое число.
Velocity Spectrum Width (cm) – ширина части графического окна в см
для отрисовки спектра скоростей.
Целое число.
Moveout Time Step (msec) – Шаг расчета кинематических поправок по времени в мсек.
Minimun Offset – Минимальное удаление в каждой мигрированной сейсмограмме. Неотрицательное число.
Maximum Offset – Максимальное удаление в каждой мигрированной сейсмограмме. Положительное число.
Offset Step – Шаг удалений в каждой мигрированной сейсмограмме (метры).
Каждая мигрированная сейсмограмма – это набор суммотрасс для заданного набора удалений ПП-ПВ (по абсолютной величине)
Taylor series Power – количество членов в ряде Тейлора
для аппроксимации остаточных кинематических поправок.
MVD grid сетка начальных приближений скоростей миграции во временном масштабе в формате GS Binary. Требуется, если значение Base Velocity указано равным нулю.
Внимание: если обработка проводилась от дневной поверхности, то и сетка скоростей должна быть от дневной поверхности.
Ignore Init Trace Filtering – При включенном переключателе будет
проигнорирована Предварительная нестационарная фильтрация входных трасс, что в десятки раз ускоряет процедуру.
При нажатии кнопки OK в цикле по пикетам ОГТ проводится обработка сейсмограмм в точках профиля и, при успешном завершении, появляется Графическое окно вертикального скоростного анализа.
Внимание: если обработка проводилась от дневной поверхности, то и определяемые скорости будут отсчитываться в масштабе двойного вертикального времени от дневной поверхности.