способ построения сейсмического глубинного разреза

Формула изобретения

Способ построения сейсмического глубинного разреза, включающий развертку отображаемых точек (ОТ) разреза A(x,h) с заданным шагом по столбцам (h) и/или строкам (х), из каждой ОТ (x₀ ,h₀) производят обзор волнового поля для поиска амплитуд сигналов А сейсмических границ, возможно принадлежащих реперным и/или целевым горизонтам, на n сейсмограммах, зарегистрированных способом многократного профилирования, обеспечивающего засветку под различными углами отображаемых точек из n пунктов воздействий, предполагая наличие сейсмической границы в отображаемых точках; обзор реализуют серией решений квазипрямой задачи лучевого трассирования, находят наиболее вероятный эйконал, используя серию виртуальных годографов для получения маркерных отметок времени на сейсмотрассах; генерацию годографов производят в соответствии с предварительно заданным множеством скоростей V=V₀± V(1, 2, , /2, , ), (где V₀ - априорно заданное первоначальное значение скорости из заданного массива скоростей V(x,h), x, h - координаты глубинного разреза, V шаг приращения скорости, - число шаговых приращений скорости), а также с предварительно заданной определенной конфигурацией сейсмического объекта, который задан q ракурсами в отображаемой точке с определенным шагом поворота угла падения с центром вращения в отображаемой точке, для каждого ракурса объекта задают множество углов падения =± ₀± (1, 2, , q/2, , q), (где ₀ - априорно заданное первоначальное значение угла падения границы в отображаемой точке, - шаг приращения угла падения, q - число пошаговых приращений угла падения), задают координаты отображаемой точки, n пунктов воздействий, m сейсмоприемников и закон (функцию) расположения последовательности из m точек, прилегающих к отображаемой точке; для каждого заданного типа волны строят q mn виртуальных годографов, производят направленное псевдосинфазное накопление амплитуд многокомпонентных цугов колебаний в суммоцугах на зарегистрированных m сейсмотрассах вдоль маркерных отметок времени на годографах, причем для каждого многокомпонентного цуга колебаний считывают компонентов (обычно 20) в виде функционалов F ; после накопления в каждом из суммоцугов в виде функционалов F _m по m сигналов их группируют по признаку равенства номера ракурса q; после введения известными способами поправок на углы засветки и мьютинга, причем в соответствии с интерактивным режимом интерпретации суммоцугов отдельной сейсмограммы и коррекции заданных диапазонов V, и , q для каждой из скоростей и каждого из q суммоцугов известным способом после вычисления q функционалов статистической обработки F _m строят для каждой ОТ двумерную матрицу из элементов их энергии F^W _m( ,V), после соответствующих вычислений в зависимости от априорных углов падения и априорно заданного множества скоростей V по столбцам матрицы документируют значения функционалов в зависимости от угла падения в отображаемой точке, а по строкам - значения функционала в зависимости от скорости (или наоборот); строят отображение графика из матрицы F^W _m( ,V), на котором визуализируют значения функционалов в виде изолиний, гистограмм или иных известных способов; график представляют в виде окна обзора достоверности и обеспечивают контроль интерпретатора за процессом опознания сейсмических границ и построением сейсмического глубинного разреза; далее уровень F^W _m( ,V), в пределах окна обзора достоверности сравнивают с заданным порогом обнаружения, при превышении порога обнаружения каждый последующий уровень отображаемого функционала сравнивают с предыдущим, при превышении последующего уровня над предыдущим его документируют в блоке максимума функционала статистической обработки F^max( _max,q_max), в котором также хранят соответствующее максимальному уровню аргументы _max, и q_max, изображают отображаемую точку на глубинном разрезе с интенсивностью, обусловленной найденными параметрами F^max, _max и q_max в пределах окна обзора достоверности, переходят к документированию последующей отображаемой точки на глубинном разрезе, присваивают F^max положение F( /2,q/2) в центре матрицы окна обзора достоверности и определяют скорректированные значения углов падения и скоростей по смещению F^max относительно центра окна обзора достоверности, для чего сравнивают аргументы /2 и q/2 с _max и q_max; величину и знак скорректированных значений для V₀ и ₀ определяют по разностям:

V_0i =V_0(i-i)± V=V_0(i-1)± V( /2- _max);



где индексами i и (i-1) обозначены последующие и предыдущие априорные и текущие начальные значения диапазона изменения скорости и угла падения в отображаемой точке; последовательно в блоке F^max( _max, q_max) заменяют в автоматическом режиме предыдущие значения диапазона на скорректированные, а монотонно возрастающие с глубиной значения скорости заменяют лишь при условии совпадения знаков развертки отображаемой точки по столбцам со знаком приращения V скорости V, тем самым размещают F(q_(i-1), _(i-1)) в центр новой скорректированной сетки значений матрицы окна обзора достоверности, для последующих q× углов падения и скоростей находят скорректированные значения F(q_i, _i), среди них в пределах следующего окна обзора достоверности находят новое текущее значение F^max, документируют очередную отображаемую точку с интенсивностью, обусловленной найденными параметрами F^max, _max и q_max, размещают F^max вновь в центр очередной скорректированной сетки окна обзора достоверности, и так далее до полного построения глубинного сейсмического разреза в соответствии с заданным режимом развертки отображаемых точек по строкам и/или столбцам.

Описание изобретения к патенту

Предлагаемое изобретение относится к сейсмической разведке и может быть использовано для построения изображений глубинного разреза.

При проведении сейсморазведки сложно построенных сред с наклонными и криволинейными границами раздела последние на временном разрезе отображаются со сносом относительно их истинного расположения и не обеспечивают подобия временного разреза по отношению к глубинному. Поэтому более 50% скважин для добычи нефти и газа оказываются пустыми. Для исключения сноса из временного разреза для несложно построенных сред изредка удается получить малодостоверный глубинный разрез посредством применения процедур миграции, однако ситуации с малодостоверным прогнозом нефтегазоносности это не спасает [Бондарев В.И. Сейсморазведка (учебник). Екатеринбург. 2007, с.467-482] [1].

Известный способ построения глубинного разреза с применением метода регулируемого направленного приема (МРНП) основан на накоплениях суммоцугов посредством прямолинейного разновременного суммирования интерференционных волн. [Сейсморазведка. Справочник геофизика. - М.: Недра. 1981, с.222-223] [2].

Недостатком этого способа является то, что он неприменим для накоплений суммоцугов вдоль годографов любого заданного типа сейсмических волн.

Указанный выше аналог выбран в качестве прототипа.

Задачей изобретения является разработка способа, позволяющего обеспечить возможность построения высокодостоверного глубинного разреза с использованием единичной сейсмограммы, зарегистрированной после сейсмического воздействия из одного из пунктов на профиле наблюдений.

Технический эффект состоит в обеспечении экспресс-обработки сейсморазведочных данных в первую очередь в полевых условиях, увеличении темпов проведения геологоразведочных работ и уменьшении затрат ресурсов.

Данный технический результат достигается за счет того, что в известном способе построения сейсмического глубинного разреза, включающем развертку отображаемых точек (ОТ) разреза A(x,h) с заданным шагом по столбцам (h) и/или строкам (х), из каждой OT(x₀,h₀) производят обзор волнового поля для поиска амплитуд сигналов А сейсмических границ, возможно принадлежащих реперным и/или целевым горизонтам, на n сейсмограммах, зарегистрированных способом многократного профилирования, обеспечивающего засветку под различными углами отображаемых точек из n пунктов воздействий, предполагая наличие сейсмической границы в отображаемых точках; обзор реализуют серией решений квазипрямой задачи лучевого трассирования, находят наиболее вероятный эйконал, используя серию виртуальных годографов для получения маркерных отметок времени на сейсмотрассах; генерацию годографов производят в соответствии с предварительно заданным множеством скоростей V=V₀± V(1, 2, , )/2, , ), (где V₀ - априорно заданное первоначальное значение скорости из заданного массива скоростей V(x,h), x,h - координаты глубинного разреза, V шаг приращения скорости, - число шаговых приращений скорости), а также с предварительно заданной определенной конфигурацией сейсмического объекта, который задан q ракурсами в ОТ с определенным шагом поворота угла падения с центром вращения в отображаемой точке, для каждого ракурса объекта задают множество углов падения =± ₀± (1, 2, , q/2, , q), (где ₀ - априорно заданное первоначальное значение угла падения границы в отображаемой точке, - шаг приращения угла наклона, q - число пошаговых приращений угла падения), задают координаты отображаемой точки, n пунктов воздействий, m сейсмоприемников и закон (функцию) расположения последовательности из m точек, прилегающих к отображаемой точке; для каждого заданного типа волны строят q mn виртуальных годографов, производят направленное псевдосинфазное накопление амплитуд многокомпонентных цугов колебаний в суммоцугах на зарегистрированных m сейсмотрассах вдоль маркерных отметок времени на годографах, причем для каждого многокомпонентного цуга колебаний считывают р компонентов (обычно p 20) в виде функционалов F ; после накопления в каждом из суммоцугов в виде функционалов F _m по m сигналов их группируют по признаку равенства номера ракурса q; после введения известными способами поправок на углы засветки и мьютинга, согласно предлагаемому способу, в соответствии с интерактивным режимом интерпретации суммоцугов отдельной сейсмограммы и коррекции заданных диапазонов V, и , q для каждой из d скоростей и каждого из q суммоцугов известным способом после вычисления q функционалов статистической обработки F _m строят для каждой ОТ двумерную матрицу из элементов их энергии F^W _m( ,V), после соответствующих вычислений в зависимости от априорных углов наклона и априорно заданного множества скоростей V по столбцам матрицы документируют значения функционалов в зависимости от угла наклона в отображаемой точке, а по строкам - значения функционала в зависимости от скорости (или наоборот); строят отображение графика из матрицы F^W _m( ,V), на котором визуализируют значения функционалов в виде изолиний, гистограмм или иных известных способов; график представляют в виде окна обзора достоверности и обеспечивают контроль интерпретатора за процессом опознания сейсмических границ и построением сейсмического глубинного разреза; далее уровень F^W _m( ,V), в пределах окна обзора достоверности сравнивают с заданным порогом обнаружения, при превышении порога обнаружения каждый последующий уровень отображаемого функционала сравнивают с предыдущим, при превышении последующего уровня над предыдущим его документируют в блоке максимума функционала статистической обработки F^max( _max,q_max), в котором также хранят соответствующее максимальному уровню аргументы _max, q_max, изображают отображаемую точку на глубинном разрезе с интенсивностью, обусловленной найденными параметрами F^max _max и q_max, переходят к документированию последующей отображаемой точки на глубинном разрезе, присваивают F^max положение F( /2,q/2) в центре матрицы окна обзора достоверности и определяют скорректированные значения углов падения и скоростей по смещению F^max относительно центра окна обзора достоверности, для чего сравнивают аргументы /2, q/2 с _max, q_max; величину и знак скорректированных значений для V₀ и ₀ определяют по разностям:

V _0i=V_0(i-i)± V=V_0(i-1)± V( /2- _max);

где индексами i и (i-1) обозначены последующие и предыдущие априорные и текущие начальные значения диапазона изменения скорости и угла падения в отображаемой точке; последовательно в блоке F^max( _max,q_max) заменяют в автоматическом режиме предыдущие значения диапазона на скорректированные, а монотонно возрастающие с глубиной значения скорости заменяют лишь при условии совпадения знаков развертки ОТ по столбцам со знаком приращения V скорости V, тем самым размещают F(q_(i-1), V_(i-1)) в центр новой скорректированной сетки значений матрицы окна обзора достоверности, для последующих q× углов падения и скоростей находят скорректированные значения F(q_i, _i), среди них в пределах следующего окна обзора достоверности находят новое текущее значение F^max, документируют очередную отображаемую точку с интенсивностью, обусловленной найденными параметрами F^max _max и q_max, размещают F^max вновь в центр очередной скорректированной сетки окна обзора достоверности, и так далее до полного построения глубинного сейсмического разреза в соответствии с заданным режимом развертки отображаемых точек по строкам и/или столбцам.

Способ осуществляется следующим образом.

Построение глубинного разреза A(x,h) осуществляется разверткой отображаемых точек (ОТ) глубинного разреза по столбцам (h) и/или строкам (х). Из каждой OT(x ₀,h₀) производят обзор волнового поля на сейсмограммах, зарегистрированных способом многократного профилирования, для поиска амплитуд сигналов А сейсмических границ, возможно принадлежащих реперным и/или целевым горизонтам. Обзор реализуется серией решений квазипрямой задачи лучевого трассирования с использованием виртуальных годографов для последующего обнаружения наиболее вероятного эйконала волны заданного типа. Задают определенную конфигурацию сейсмического объекта и при каждом облучении отображаемой точки из пункта воздействия задают q ракурсов сейсмического объекта с определенным шагом поворота угла наклона с центром вращения в отображаемой точке. Для каждого ракурса объекта при =± ₀± (1, 2, , q/2, , q), (где ₀ - априорно заданное первоначальное значение угла падения границы в ОТ, - шаг приращения угла наклона, q - число пошаговых приращений угла наклона (ракурсов), и каждого из V=±V₀±Aa(1, 2, , /2, , ) значений скорости (где V₀ - априорно заданное первоначальное значение скорости из заданного массива скоростей V(x,h), x, h - координаты глубинного разреза, V шаг приращения скорости, v - число шаговых приращений скорости), а также координат ОТ, ПВ, СП и заданного закона (функции) расположения последовательности из m точек, прилегающих к отображаемой точке, рассчитывают q mn виртуальных годографов заданного типа сейсмической волны.

Определяют длины m эйконалов для каждой из n сейсмограмм, находят по m маркерных отметок времени t_µ прихода лучей на сейсмограммах, зарегистрированных от n ПВ, производят направленное синфазное накопление амплитуд, принадлежащих каждому из m цугов колебаний функционалов статистической обработки (ФСО) F , из которых на зарегистрированных m сейсмотрассах вдоль маркерных отметок времени из m годографов накапливают ФСО F _m, размещаемых в блоке суммоцугов (СЦ) в соответствии с номерами ПВ, сигналы группируют по признаку равенства номера ракурса q, для каждой из v скоростей и каждого из q супернакоплений амплитуд по фаз каждого из mn суперсигналов, специфичных для заданных позиций объекта в отображаемой точке; известным способом в необходимых случаях вычисляют q ФСО F _m, согласно предлагаемому способу строят для каждой ОТ двумерную матрицу из элементов их энергий F^W _m( ,V), после соответствующих вычислений в зависимости от априорных углов наклона и априорно заданного множества скоростей V, по столбцам матрицы документируют значения ФСО в зависимости от угла наклона в отображаемой точке, а по строкам - значения функционала в зависимости от скорости (или наоборот); строят отображение графиков из необходимого количества вышеперечисленных матриц. Интерактивный режим обнаружения волн заданного типа, интерпретации и коррекции априорно заданных диапазонов V,v и ,q осуществляют построением окон обзора достоверности (ООД) на основе визуализации до n двумерных матриц из элементов F _m, после вычисления их энергий F^W _m( ,V) в виде изолиний, гистограмм или иных способов, после чего визуализируют матрицы из сейсмотрасс от других ПВ, обеспечивают выбор наиболее информативных ООД и контроль интерпретатора за процессом опознания сейсмических границ, находят, кроме волн заданного типа, области распределения в пределах ООД иных типов волн, обеспечивают возможность осуществлять операции сложения, вычитания, умножения матриц, вводят известными способами в массивы матриц поправочные коэффициенты на углы засветки и исключают из процесса построения зашумленные эйконалы.

В пределах наиболее информативной ООД уровень каждого элемента матрицы сравнивают с заданным порогом обнаружения, при превышении порога обнаружения каждый последующий уровень сравнивают с предыдущим, при превышении последующего уровня над предыдущим его документируют в блоке F_max( _max,q_max), в котором также хранят соответствующее максимальному уровню _max и q_max. Отображают ОТ на глубинном разрезе с интенсивностью, обусловленной найденными параметрами F^max, _max и q_max. Для получения в соответствии с заданным режимом развертки ОТ по строкам и/или столбцам глубинного разреза переходят к документированию последующей ОТ на глубинном разрезе. Для этого присваивают F^max положение W( /2,q/2) в центре матрицы ООД и определяют скорректированные значения углов наклона и скоростей по смещению F^max относительно центра окна обзора достоверности, для чего сравнивают аргументы /2 и q/2 с _max и q_max; величину и знак скорректированных значений для V₀ и ₀ определяют по разностям:

V _0i=V_0(i-i)± V=V_0(i-1)± V( /2- _max);

где индексами i и (i-1) обозначены последующие и предыдущие априорные и текущие начальные значения диапазона изменения скорости и угла наклона в ОТ; последовательно в блоке F^max( _max,q_max) заменяют в автоматическом режиме предыдущие значения диапазона углов падения на скорректированные, а монотонно возрастающие с глубиной значения скорости заменяют лишь при условии совпадения знаков развертки ОТ по столбцам со знаком приращения V скорости V, тем самым размещают F(q_(i-1), _(i-1)) в центр новой скорректированной сетки значений матрицы ООД, обнуляют блок F^max( _max,q_max), снова для всех q× находят скорректированные значения F(q_i, _i), среди них в пределах следующей ООД находят новое текущее значение F^max, документируют очередную ОТ, на глубинном разрезе с интенсивностью, обусловленной найденными параметрами F^max, _max и q_max, размещают F^max вновь в центр очередной скорректированной сетки ООД, и так далее до полного построения способом обзора волнового глубинного разреза.

Пути реализации патентуемого способа приводятся в нижеследующем описании примера 1.

Пример 1. Блок-схема реализации способа построения глубинного сейсмического разреза. Описание поясняется Фиг.1, где:

Х - координатная ось вдоль профиля регистрации;

ПВ1, ПВ2, ПВ3, , ПВn - пункты возбуждения сейсмических сигналов;

OT(x₀,h₀) - отображаемая точка под n углами засветки из ПВ1, ПВ2, ПВ3, , ПВn;

1 - последовательность предполагаемых отрезков сейсмической границы под углами падения, ₀, ₀+ , ₀+2 , , ₀+q ;

t - координатная ось времени регистрации сейсмических сигналов на сейсмограмме;

2 - блок генерации mnq виртуальных годографов по амплитуд волны заданного типа по известным координатам ОТ, ПВ, СП, заданным перечнем углов наклона и скоростей V;

СП1, СП2, , СПm - сейсмоприемники на профиле регистрации сейсмограммы;

3 - отображение одного из mnqv виртуальных годографов, полученного из блока генерации 2;

t_µ - маркерные виртуальные отметки времени выемки отсчетов амплитуд фазы сейсмического сигнала на сейсмотрассе, соответствующие времени хода лучей вдоль эйконала (падающего из ПВ до ОТ и восстающего от ОТ до СП луча);

4 - блоки хранения отсчетов из каждого цуга колебаний;

5 - блок хранения текущих накоплений в суммоцуге F _m из m цугов колебаний по m амплитуд в каждом из суммоцугов (СЦ);

6 - блок СЦ отдельного хранения n суммоцугов, в которых размещают ФСО F _m из сейсмограмм в соответствии с номерами ПВ;

7 - субблоки суммоцугов, полученных в результате псевдосинфазного накопления ФСО F _m по m амплитуд вдоль каждого из mnqv виртуальных годографов 3;

8 - блок вычисления значения ФСО типа энергии супернакоплений F^W( ,V) для текущих значений углов падения отрезков q и скорости ;

9 - блок нормировки и масштабирования;

10, 11 и 12 - окна обзора достоверности;

13 - блок сравнения текущего значения ФСО с пороговым уровнем шумов;

14 - изображение сейсмического глубинного разреза;

15 - блок сравнения текущего значения ФСО с предыдущим максимальным значением и засылки большего из них в блок 15;

16 - блок хранения максимального значения ФСО F^max и соответствующего этому значению угла падения и величины скорости.

Построение сейсмического разреза в соответствии с заявляемым способом осуществляется, например, следующим образом.

На фиг.1 показана ОТ под n углами засветки из ПВ1, ПВ2, , ПВn и СП1, СП2, СП3, , СПm, расположенных на профиле регистрации и изображенных на координатной оси X. На отрезках 1 сейсмической границы под углами падения ₀, ₀+ , ₀+2 , , ₀+ ×q предположительно находятся m точек, прилегающих к данной ОТ и являющиеся источниками зарегистрированных СП вторичных волн.

Вероятность обнаружения точек отражения вдоль заданных границ обеспечивают обзором волнового поля каждой из n сейсмограмм, зарегистрированных в результате воздействий из ПВ1, ПВ2, , ПВn. Обзор реализуют серией решений квазипрямой задачи сейсморазведки по способу лучевого трассирования. Из блока 2 осуществляют генерацию q mn виртуальных годографов 3 с использованием координат ОТ, ПВ, СП, заданных величин шага и количества (m) прилегающих к ОТ точек, текущих величин углов падения в ОТ отрезков (до q) предполагаемых сейсмических границ и заданных (до ) скоростей. V

В зависимости от характера распространения падающих и восстающих лучей из блока 2 сначала с учетом величины скорости V₀ и угла падения ₀ определяют длины m эйконалов, находят по m маркерных отметок времени t_µ прихода лучей на сейсмограмме, зарегистрированной от ПВ1.

Для каждого из q текущих углов падения и каждой из v текущих значений скорости до и/или после маркерных отметок времени накапливают отсчетов амплитуд цуга колебаний в массивах F , исходящих из отрезков 1, и последовательно размещают их в m блоках цугов колебаний 4. Производят псевдосинфазные первичные накопления в массиве F _m кратностью m, который размещают в блоке суммоцуга 5, размещают его в блоке 6 в одном из массивов 7, соответствующем ПВ1. В блоке 8 определяют для ФСО F _m величину их энергии F^W _m( ,V) и при необходимости размещают в качестве одного из элементов матрицы F^W _m( ,V) системы интерактивной интерпретации энергии ФСО.

В блоке 9 производят нормировку и масштабирование ФСО ( ₀,V₀), после чего засылают, если возникает необходимость в визуальном контроле, ФСО ( ₀,V₀) в качестве первого элемента матрицы 10 окна обзора достоверности (ООД) для последующей визуализации либо в виде поля изолиний 10, либо в виде гистограмм 11 энергии суммоцуга отдельной сейсмограммы. Энергию ФСО ( ₀,V₀) сравнивают затем в блоке 12 с пороговым значением уровня шумов. Если уровень ФСО менее порогового значения, то на изображении ОТ на глубинном разрезе 13 регистрируют нулевую интенсивность. Если уровень ФСО превышает порог, то он засылается в блок 14, в котором сравнивают текущее значение ФСО с предыдущим максимальным значением F^max, хранящемся в блоке 15 хранения максимального значения ФСО F^max . Если текущее значение F^max превышает предыдущее, хранящееся в блоке 15, то содержимое блока 14 заменяют на большее значение F^max и засылают большую из них обратно в блок 15 (перед засылкой первого элемента матрицы, следовательно, блок 15 обнулен). В блок 15 засылают при этих логических условиях также найденную величину наиболее вероятного угла наклона _0i и скорости V_0i.

Затем на сейсмограмме, зарегистрированной от ПВ2, также вдоль соответствующего годографа накапливают данные из цугов 4 в m суммоцугах 6 и из каждой отдельной сейсмограммы размещают их в n субблоках 7 в виде ФСО F _m, который также может использоваться в качестве первого элемента другой матрицы (из сейсмограммы ПВ2). Подобным же образом через блоки 4-6 накапливают данные в остальных субблоках 7 из всех заданных (до n) сейсмограмм ФСО F_pm. В блоке 8 вычисляют энергию F^W _m( ₀,V₀) системы интерактивной интерпретации для последующей визуализации в блоках 10 и/или 11 известными способами и так далее до n раз.

Аналогичным образом приступают при V=V₀ к генерации следующей серии из mn виртуальных годографов, выборке данных из n сейсмограмм, засылке и обработке их в блоках 4÷15 для каждого из остальных углов падения отрезков 1: =а₀+2 и так далее до = ₀+q , при котором полностью заполняют первые строки до n матриц ООД.

Приступают к генерации виртуальных годографов при следующем значении скорости V=V₀± V и аналогично заполняют вторые строки матриц ООД и так далее до заполнения последних строк матриц ООД при V=V₀ + V и = ₀+q .

На каждое полученное таким образом (из n) изображение глубинного разреза 13 в соответствии с выбранным способом отображения засылают величины максимума ФСО из блока 15 и наиболее вероятного угла падения _0i отрезка в ОТ.

Далее при необходимости аналогичным образом осуществляют интерактивный режим интерпретации отдельных сейсмограмм, хранящихся в блоке 6 СЦ и зарегистрированных другими ПВ, коррекцию заданных диапазонов V, и , q, построением ООД на основе визуализации энергий до n матриц F_pm( ,V) в виде изолиний, гистограмм или иных способов

Документирование последующих ОТ на глубинном разрезе осуществляют в соответствии с заданным режимом развертки ОТ по строкам и/или столбцам глубинного разреза. Для каждой последующей ОТ определяют скорректированные значения углов наклона и скоростей путем задания нового центра окна обзора достоверности с координатами F( /2,q/2) и сравнения его аргументов и/2 и q/2 с аргументами предыдущего _max и q_max; величину и знак скорректированных значений для V₀ и ₀ определяют по разностям:

V _0i=V_0(i-i)± V=V_0(i-1)± V( /2- _max);

где индексами i и (1-1) обозначены последующие и предыдущие априорные и текущие начальные значения диапазона изменения скорости и угла падения в ОТ; последовательно заменяют в автоматическом режиме предыдущие значения диапазона углов падения на скорректированные, а монотонно возрастающие с глубиной значения скорости заменяют лишь при условии совпадения знаков развертки ОТ по столбцам со знаком приращения V скорости V, тем самым размещают в F^max(q _(i-1), _(i-1)) центр новой скорректированной сетки значений матрицы ООД, снова для всех q× находят скорректированные значения F(q_i, _i), среди них в пределах ООД находят текущее значение F^max, документируют очередную ОТ с интенсивностью, обусловленной найденными параметрами F^max, _max и q_max, размещают F^max вновь в центр очередной скорректированной сетки ООД. Подобным способом обзора волнового поля наносят все остальные соседние ОТ на глубинный разрез с заданным шагом приращения координат.

Использованные источники информации

1. Бондарев В.И. Сейсморазведка (учебник). - Екатеринбург. 2007, с.467-482.

2. Сейсморазведка. Справочник геофизика. - М.: Недра, 1981.