Способ построения сейсмического динамического глубинного и/или тотального временного разреза

Способ построения сейсмического динамического глубинного и/или тотального временного разреза

Иллюстрации

Показать все

Реферат

Изобретение относится к сейсмической разведке и может быть использовано для построения изображений сложно построенных сред в виде динамического глубинного K(x, h) и/или тотального (совокупного) временного разреза T(x, t) трансформацией сигналов предварительного глубинного разреза A(x, h).

Временной разрез является частным случаем поля времен, зарегистрированного способом многократного профилирования, при котором расстояние между источником и приемником равно нулю [Сейсморазведка. Справочник геофизика. - М.: Недра, 1981, с.85] /1/. Временной разрез, полученный по способу метода общей глубинной точки (МОГТ), необходим в первую очередь для получения атрибутов прямых поисков нефти и газа [Бондарев В.И. Сейсморазведка (учебник). Екатеринбург. 2007. с.467-482] /2/, обеспечения фильтрации, динамической и спектральной обработки сигналов, корреляции волн.

Однако при использовании МОГТ для построения разрезов сложно построенных сред с наклонными и криволинейными границами раздела последние на временном разрезе отображаются со сносом относительно их истинного расположения и не обеспечивают подобия временного разреза по отношению к глубинному. Поэтому более 50% скважин для добычи нефти и газа оказываются пустыми. Для исключения сноса из временного разреза изредка удается получить малодостоверный глубинный разрез посредством применения процедур миграции, однако ситуации с малодостоверным прогнозом нефтегазоносности это не спасает.

Известный способ, защищенный патентом №2415449 от 27 марта 2011 г. /3/, «Способ построения сейсмического глубинного разреза», имеющий общность с заявляемым способом по совокупности признаков, обеспечивает высокую достоверность и точность отображения сейсмических объектов и сред большой сложности в виде глубинного разреза A(x, h). В связи с большой сложностью изложения предлагаемого изобретения сущность способа далее излагается посредством применения более формализованных названий накоплений сигналов в виде функционалов и математических обозначений и состоит в том, что для получения изображения сейсмического глубинного разреза A(x, h) каждую отображаемую точку (ОТ) последовательно наносят на глубинный разрез вдоль столбцов и/или строк с заданным шагом по осям координат в результате обзора интенсивности волнового поля A_n(x, t) в этой ОТ на n сейсмограммах, зарегистрированных способом многократного профилирования. Обзор производят методом лучевого трассирования и поиском сигналов сейсмических границ, возможно пересекающих/вмещающих ОТ и принадлежащих одному из q заданных ракурсов объекта (сейсмогеологической модели); для каждого ракурса объекта при априорно заданном скоростном разрезе для каждой из n сейсмограммы строят q виртуальных годографов заданного типа волны в соответствии с трассами хода лучей вторичных волн, возможно исходящих из ОТ и прилегающей к ОТ последовательности из m точек вдоль падающей в ОТ под углом α границы, после чего из каждого цуга колебаний до и после маркерных отметок времен прихода лучей вдоль каждого годографа осуществляют псевдосинфазное считывание и накопление по ρ (обычно ρ≤20) отсчетов амплитуд сигналов в виде функционала F_ρ и из каждой из заданного количества m цугов колебаний на каждой из n зарегистрированных сейсмограмм псевдосинфазно накапливают по ρ×m сигналов в q суммоцугах (СЦ) и получают функционалы статистической обработки F_ρm. Суммарные результаты накоплений в СЦ из n сейсмограмм группируют по признаку принадлежности к каждой из q позиций и после псевдосинфазного суммирования/умножения функционалов F_ρm по ρ×m амплитуд колебаний из всех n сейсмограмм для одной и той же позиции объекта получают гиперсигнал (ГС) с ρ×m×n кратностью накоплений в виде функционала F_ρmn. Из всех заданных q позиций объекта находят наиболее вероятную путем сравнения функционалов статистической обработки F_ρmn и/или F_mn (при n=1) между собой, получают максимум вероятности величины угла наклона и обнаружения F^max. Последнюю сравнивают с заданным порогом уровня шумов, нормируют с учетом масштаба изображения и отображают точку/отрезок сейсмической границы на глубинном разрезе пропорционально величине F^max. В необходимых случаях корректируют заданные параметры обзора для достижения максимума величины вероятности обнаружения границы. Подобным способом обзора волнового поля наносят все остальные соседние ОТ на глубинный сейсмический разрез A(x, h) в соответствии с заданным шагом приращений координат.

Наиболее близким к предлагаемому является способ, применяемый при получении изображений сейсмического глубинного разреза, защищенный патентом №2415449 от 27 марта 2011 г. /3/. Данный способ избран в качестве прототипа.

Недостатком этого способа является то, что, автоматизация определения по максимуму вероятности обнаружения достигается лишь в отношении угла наклона в ОТ, в то время как для всей области обработки необходимо априорное задание скоростного разреза. Кроме того, при его применении недостаточно обеспечена необходимая точность и достоверность динамического отображения волнового поля глубинного и тотального (совокупного) временного разреза, недостаточно обеспечено получение качественных атрибутов прямых поисков нефти и газа, результатов пространственной и временной фильтрации, спектральной и динамической обработки сигналов, так как при автокорреляции, корреляции, свертке, вычислениях энергии используются до ρ отсчетов сигналов суммоцугов для каждой фазы после их усреднения, в результате чего исходное волновое поле существенно искажается.

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

Технический результат состоит в исключении сноса отображаемых на тотальном временном разрезе сейсмических границ относительно их расположения под пунктами регистрации наблюдений, комплексная автоматизация определения угла наклона сейсмических границ и монотонно возрастающей/убывающей скорости в ОТ, повышении диапазона динамических характеристик волнового поля глубинного и тотального временного разрезов.

Предлагаемый способ построения сейсмического динамического глубинного K(x, h) и/или тотального временного разреза T(x, t) и прототип имеют следующие общие признаки.

Для получения изображения предварительного сейсмического глубинного разреза A(x, h) каждую отображаемую точку (ОТ) последовательно наносят на глубинный разрез вдоль столбцов и/или строк с заданным шагом по осям координат в результате обзора интенсивности волнового поля в этой ОТ на n сейсмограммах, зарегистрированных способом многократного профилирования. Обзор производят методом лучевого трассирования и поиском сигналов сейсмических границ, возможно пересекающих/вмещающих ОТ и принадлежащих одному из q заданных ракурсов объекта (сейсмической границы, которая задана в отображаемой точке q ракурсами с определенным шагом поворота Δα угла падения α с центром вращения в отображаемой точке), для каждого ракурса объекта задают множество углов падениия α=±α₀±Δα (1, 2, …, q/2, …, q), (где α₀ - априорно заданное первоначальное значение угла падения границы в отображаемой точке, Δα - шаг приращения угла падения, q - число пошаговых приращений угла падения), а также предварительно установленной определенной конфигурации сейсмического объекта, для каждого ракурса объекта из n сейсмограмм строят q виртуальных годографов заданного типа волны в соответствии с трассами хода лучей вторичных волн, возможно исходящих из ОТ и прилегающей к ОТ последовательности из m точек вдоль падающей в ОТ под углом α границы, после чего из каждого цуга колебаний до и после маркерных отметок времен прихода лучей вдоль каждого годографа осуществляют псевдосинфазное считывание и накопление по ρ (обычно ρ≤20) отсчетов амплитуд сигналов в виде функционала F_ρ и из каждой из заданного количества m цугов колебаний на каждой из n зарегистрированных сейсмограмм псевдосинфазно накапливают по ρ×m сигналов в q суммоцугах (СЦ) и получают функционалы статистической обработки F_ρm. Суммарные результаты накоплений в СЦ из n сейсмограмм группируют по признаку принадлежности к каждой из q позиций и после псевдосинфазного суммирования/умножения функционалов F_ρm по ρ×m амплитуд колебаний из всех n сейсмограмм для одной и той же позиции объекта получают гиперсигнал (ГС) с ρ×m×n кратностью накоплений в виде функционала F_ρmn. Из всех заданных q позиций объекта находят наиболее вероятную путем сравнения функционалов статистической обработки F_ρmn и/или F_mn (при n=1) между собой, получают максимум вероятности величины угла наклона и обнаружения F^max. Последнюю сравнивают с заданным порогом уровня шумов, нормируют с учетом масштаба изображения и отображают точку/отрезок сейсмической границы на глубинном разрезе пропорционально величине F^max. В необходимых случаях корректируют заданные параметры обзора для достижения максимума величины вероятности обнаружения границы. Подобным способом обзора волнового поля наносят все остальные соседние ОТ на глубинный сейсмический разрез A(x, h) в соответствии с заданным шагом приращений координат.

Отличительные признаки предлагаемого изобретения:

Построение сейсмического динамического глубинного K(x, h) и/или тотального временного разреза T(x, t) производят на основе развертки по столбцам (оси h) и/или строкам (оси х), например энергии отображаемых точек (ОТ) заданного типа волны в пределах предварительного глубинного разреза A(x, h). Для этого в максимальной степени детализируют фазы сейсмических сигналов во временном масштабе следующим образом.

Сначала для каждой отображаемой точки OT(x₀, h} в соответствии с процедурами построения дополнительного (предварительного) глубинного разреза A(x, h) осуществляют обзор волновых полей A_n(x, t) на сейсмограммах, полученных методом многократного профилирования, посредством лучевого трассирования виртуальных годографов на n сейсмограммах, зарегистрированных m сейсмоприемниками (СП), и обеспечивающего засветку под различными углами отображаемых точек из 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 сейсмотрассах вдоль маркерных отметок времени на годографах.

Затем в соответствии с процедурами, необходимыми для построения динамических глубинного и тотального (совокупного) временного разрезов для каждого цуга колебаний считывают по единичному отсчету амплитуды фазы, последовательность единичных сигналов ЕС из m цугов колебаний накапливают в виде функционала K_m. Для каждого из углов засветки от n пунктов возбуждений массивы из m единичных отсчетов последовательно размещают в n субблоках суммофаз (СФ), в необходимых случаях вносят поправки на угол засветки, определяют область мьютинга, удаляют зашумленные эйконалы. Производят накопления СФ с кратностью mn в виде функционала K_mn.

После этого в соответствии с процедурами, необходимыми для построения дополнительного (вспомогательного) глубинного разреза, из цуга колебаний, зарегистрированного до и/или после маркерных отметок времени, считывают из каждой m-й сейсмотрассы по ρ компонентов (обычно ρ≤20) и накапливают их в виде псевдосинфазных функционалов F_ρ; после накопления сигналов, хранящихся в каждом из m функционалов F_ρ, их накаплиают по ρm сигналов в суммоцугах (СЦ) в виде функционалов F_ρm, затем их группируют по признаку равенства номера ракурса q; для каждой из q групп находят гиперсигнал ГС в виде функционала F_ρmn супернакоплениями из ρmn амплитуд после введения известными способами поправок на n углов засветки и мьютинга.

Для каждой ОТ из совокупности множеств скоростей (до υ) и углов наклона (до q) строят двумерную матрицу из элементов F_ρm и/или F_ρmn, и/или из элементов их энергий и/или , после соответствующих вычислений в зависимости от априорно заданных углов наклона α и априорно заданного множества скоростей V, по столбцам матрицы документируют значения функционалов в зависимости от угла наклона в отображаемой точке, а по строкам - значения функционала в зависимости от скорости (или наоборот); строят двумерное отображение графика для матриц F_ρm и/или F_ρmn и/или из значений их энергий и/или на котором визуализируют значения функционалов в виде изолиний, гистограмм или иных известных способов наглядного отображения; график представляют в виде окна обзора достоверности (ООД) и обеспечивают контроль интерпретатора за процессом опознания сейсмических границ и построением сейсмического глубинного разреза; уровень F_ρm и/или F_ρmn и/или значения их энергий и/или , в пределах окна обзора достоверности сравнивают с заданным порогом обнаружения, при превышении порога обнаружения каждый последующий уровень отображаемого функционала сравнивают с предыдущим, при превышении последующего уровня над предыдущим его документируют в блоке максимума функционала статистической обработки F^max(υ_max, q_max), в котором также хранят соответствующие максимальному уровню аргументы υ_max, и q_max, обеспечивающие экстраполяцию текущих значений пластовой скорости и угла падения границ между толщами горных пород и необходимые для следующих процедур построения динамических глубинного и тотального временного разрезов.

В соответствии с интерактивным режимом интерпретации единичных ЕС накоплений K_m из m фаз вдоль годографа отдельной сейсмограммы и суммофаз K_mn из годографов n сейсмограмм осуществляют (при необходимости) коррекцию заданных диапазонов V, υ и α, q следующим образом. Для каждой из υ скоростей и каждого из q ракурсов строят для каждой ОТ двумерную матрицу из элементов K_m и/или K_mn, по столбцам матрицы документируют значения функционалов в зависимости от угла падения в отображаемой точке, а по строкам - значения функционала в зависимости от скорости (или наоборот); строят отображение графика из матрицы K_m(α, V) и/или K_mn(α, V), на котором визуализируют значения функционалов в виде изолиний, гистограмм или иных известных способов; график представляют в виде окна обзора достоверности и обеспечивают контроль интерпретатора за процессом опознания заданного типа волны и построением сейсмического глубинного разреза.

После сигнала совпадения, который вырабатывается из схемы разрешения при условии смены в блоке F^max предыдущего значения максимума на последующий более высокий (при монотонном возрастании скорости с глубиной), а также при экстраполяции значения F^max для предыдущей ОТ на последующую (при постоянстве пластовой скорости), функционал K_mn или K_m заносят в буферный блок суммирования. Накопления из буферного блока суммирования единичных отсчетов амплитуд K_m или K_mn засылают в блок рассылки. Из блока рассылки накопления засылают в регистры блока масштабирования, где с учетом шкалы нормировки вырабатывают способ изображения ОТ на динамическом глубинном разрезе K(х_ОТ, h_OT) и способ изображения точки тотальной амплитуды фазы Т(х_ОТ, h_OT) на тотальном временном разрезе T(x, t), которую и засылают на суммотрассу временного разреза. Абсциссу х_ОТ точки тотальной амплитуды фазы задают в соответствии с абсциссой ОТ на глубинном разрезе, а ординату t_OT находят в соответствии с найденной наиболее вероятной величиной скорости V_OT из условия обеспечения минимального времени хода восстающего луча из отображаемой точки с коорднатами (х₀, h₀) до ближайшего, расположенного на профиле наблюдений, регистрирующего сейсмоприемника СП (х₀, 0). Обычно при определенных допущениях t_OT=h₀/V_OT. В соответствии с направлением развертки изображают отображаемую точку на дополнительном глубинном разрезе с интенсивностью, пропорциональной величине максимального функционала F^max заданного типа волны в пределах окна обзора достоверности. Затем выполняют документирование последующих ОТ дополнительного (предварительного) глубинного разреза A(x, h), присваивают F^max положение F(υ/2, q/2) в центре матрицы окна обзора достоверности и определяют скорректированные значения углов наклона и скоростей по смещению F^max относительно центра окна обзора достоверности, для чего сравнивают аргументы υ/2 и q/2 с υ_max и q_max; величину и знак скорректированных значений для V₀ и α₀ определяют по разностям:

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

Для дальнейшего увеличения степени достоверности построения динамического глубинного и/или тотального временного разреза из последовательности единичных отсчетов ЕС в блоке рассылки выделяют и засылают в отдельные регистры блока масштабирования промежуточные и полные суммы единичных отсчетов фазы сейсмических сигналов, тем самым обеспечивают возможность составления систем линейно независимых уравнений относительно неизвестной амплитуды сигнала и среднестатистической амплитуды случайной помехи, экстраполируют линейно независимые уравнения для каждой из промежуточных сумм также и на выражения для полной суммы на основе предположения о равенстве математического ожидания случайных помех для полной и частичной последовательности единичных отсчетов, решают уравнения относительно амплитуды линейно возрастающего пропорционально числу слагаемых сигнала и возрастающей в соответствии с корнем квадратным из числа слагаемых среднестатистической амплитуды случайной помехи, документируют на ОТ динамического глубинного разреза и суммотрассах тотального временного разреза найденную высокодостоверную амплитуду свободного от случайной помехи сигнала, используют найденную высокодостоверную амплитуду свободного от помехи сигнала для оценки помехозащищенности во взаимных точках и документирования тотального временного разреза.

Поставленная задача построения сейсмических динамических глубинного K(x, h) и/или тотального временного T(x, t) разрезов решена за счет того, что сначала для каждой отображаемой точки OT(x₀, h) в соответствии с процедурами построения дополнительного (предварительного) глубинного разреза A(x, h) осуществляют обзор волновых полей A_n(x, t) на сейсмограммах, полученных методом многократного профилирования. Обзор реализуют посредством лучевого трассирования виртуальных годографов на каждой из n сейсмограмм, зарегистрированных от n пунктов воздействий по способу многократных перекрытий, на которых находят сейсмотрассы, зарегистрированые m сейсмоприемниками (СП). Предполагают наличие сейсмической границы в отображаемой точке, из которой возможно исходят вторичные сейсмические волны в соответствии с предварительно заданным множеством скоростей 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 цугов колебаний накапливают в виде функционала K_m. Для каждого из углов засветки от n пунктов возбуждений массивы K_m из m единичных отсчетов последовательно размещают в n субблоках суммофаз (СФ), в необходимых случаях вносят поправки на угол засветки, определяют область мьютинга, удаляют зашумленные эйконалы. Производят накопления СФ с кратностью mn в виде функционала K_mn.

В соответствии с интерактивным режимом интерпретации единичных ЕС в виде накоплений K_m из m фаз вдоль годографа отдельной сейсмограммы и суммофаз K_mn из годографов n сейсмограмм осуществляют (при необходимости) коррекцию заданных диапазонов V, υ и α, q следующим образом. Для каждой из υ скоростей и каждого из q ракурсов строят для каждой ОТ двумерную матрицу из элементов K_m и/или K_mn, по столбцам матрицы документируют значения функционалов в зависимости от угла падения в отображаемой точке, а по строкам - значения функционала в зависимости от скорости (или наоборот); строят отображение графика из матрицы K_m(α, V) и/или K_mn(α, V), на котором визуализируют значения функционалов в виде изолиний, гистограмм или иных известных способов; график представляют в виде окна обзора достоверности и обеспечивают контроль интерпретатора за процессом опознания заданного типа волны и построением сейсмического глубинного разреза.

После этого в соответствии с процедурами, необходимыми для построения дополнительного (вспомогательного) глубинного разреза, из цуга колебаний, зарегистрированного до и/или после маркерных отметок времени, считывают из каждой m-й сейсмотрассы по ρ компонентов (обычно ρ≤20) и накапливают их в виде псевдосинфазных функционалов F_ρ; после накопления сигналов, хранящихся в каждом из m функционалов F_ρ, их накаплиают по ρm сигналов в суммоцугах (СЦ) в виде функционалов F_ρm, затем их группируют по признаку равенства номера ракурса q; для каждой из q групп находят гиперсигнал ГС в виде функционала F_ρmn супернакоплениями из ρmn амплитуд после введения известными способами поправок на n углов засветки и мьютинга.

Для каждой ОТ из совокупности множеств скоростей (до υ) и углов наклона (до q) строят двумерную матрицу из элементов F_ρm и/или F_ρmn и/или из элементов их энергий и/или , после соответствующих вычислений в зависимости от априорно заданных углов наклона α и априорно заданного множества скоростей V по столбцам матрицы документируют значения функционалов в зависимости от угла наклона в отображаемой точке, а по строкам - значения функционала в зависимости от скорости (или наоборот); строят двумерное отображение графика для матриц F_ρm и/или F_ρmn и/или из значений их энергий и/или , на котором визуализируют значения функционалов в виде изолиний, гистограмм или иных известных способов наглядного отображения; график представляют в виде окна обзора достоверности (ООД) и обеспечивают контроль интерпретатора за процессом опознания сейсмических границ и построением сейсмического глубинного разреза; уровень F_ρm и/или F_ρmn и/или значения их энергий и/или в пределах окна обзора достоверности сравнивают с заданным порогом обнаружения, при превышении порога обнаружения каждый последующий уровень отображаемого функционала сравнивают с предыдущим, при превышении последующего уровня над предыдущим его документируют в блоке максимума функционала статистической обработки F^max(υ_max, q_max), в котором также хранят соответствующие максимальному уровню аргументы υ_max, и q_max, обеспечивающие экстраполяцию текущих значений пластовой скорости и угла падения границ между толщами горных пород и необходимые для следующих процедур построения динамических глубинного и тотального временного разрезов.

После сигнала совпадения, который вырабатывается из схемы разрешения при условии смены в блоке F^max предыдущего значения максимума на последующий более высокий (при монотонном возрастании скорости с глубиной), а также при экстраполяции значения F^max для предыдущей ОТ на последующую (при постоянстве пластовой скорости), функционал K_mn или K_m заносят в буферный блок суммирования. Накопления из буферного блока суммирования единичных отсчетов амплитуд K_m или K_mn засылают в блок рассылки. Из блока рассылки накопления засылают в регистры блока масштабирования, где с учетом шкалы нормировки вырабатывают способ изображения ОТ на динамическом глубинном разрезе K(х_ОТ, h_OT) и способ изображения точки тотальной амплитуды фазы Т(х_ОТ, h_OT) на тотальном временном разрезе T(x, t), которую и засылают на суммотрассу временного разреза. Абсциссу х_ОТ точки тотальной амплитуды фазы задают в соответствии с абсциссой ОТ на глубинном разрезе, а ординату t_OT находят в соответствии с найденной наиболее вероятной величиной скорости V_OT из условия обеспечения минимального времени хода восстающего луча из отображаемой точки с коорднатами (x₀, h₀) до ближайшего, расположенного на профиле наблюдений, регистрирующего сейсмоприемника СП (х₀, 0). Обычно при определенных допущениях t_OT=h₀/V_OT. В соответствии с направлением развертки изображают отображаемую точку на дополнительном глубинном разрезе с интенсивностью, пропорциональной величине максимального функционала F^max заданного типа волны в пределах окна обзора достоверности. Затем выполняют документирование последующих ОТ дополнительного (предварительного) глубинного разреза A(x, h), присваивают F^max положение F(υ/2, q/2) в центре матрицы окна обзора достоверности и определяют скорректированные значения углов наклона и скоростей по смещению F^max относительно центра окна обзора достоверности, для чего сравнивают аргументы υ/2 и q/2 с υ_max и q_max; величину и знак скорректированных значений для V₀ и α₀ определяют по разностям:

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

Для дальнейшего увеличения степени достоверности построения динамического глубинного и/или тотального временного разреза из последовательности единичных отсчетов ЕС в блоке рассылки выделяют и засылают в отдельные регистры блока масштабирования промежуточные и полные суммы единичных отсчетов фазы сейсмических сигналов, тем самым обеспечивают возможность составления систем линейно независимых уравнений относительно неизвестной амплитуды сигнала и среднестатистической амплитуды случайной помехи, экстраполируют линейно независимые уравнения для каждой из промежуточных сумм также и на выражения для полной суммы на основе предположения о равенстве математического ожидания случайных помех для полной и частичной последовательности единичных отсчетов, решают уравнения относительно амплитуды линейно возрастающего пропорционально числу слагаемых сигнала и возрастающей в соответствии с корнем квадратным из числа слагаемых среднестатистической амплитуды случайной помехи, документируют на ОТ динамического глубинного разреза и суммотрассах тотального временного разреза найденную высокодостоверную амплитуду свободного от случайной помехи сигнала, используют найденную высокодостоверную амплитуду свободного от помехи сигнала для оценки помехозащищенности во взаимных точках и документирования тотального временного разреза.

Статистический анализ последовательностей единичных отсчетов фазы сейсмических сигналов с точки зрения теории обнаружения полезных сигналов на фоне помех допускает выделение сигнала и помехи, так как все m отсчетов амплитуды фазы сейсмического сигнала в пределах каждого из n годографов принадлежат к полезным сигналам, интенсивность которых при суммировании линейно возрастает в зависимости от числа слагаемых mn (если этот годограф зарегистрирован при проведении сейсморазведки). Сигнал случайной помехи при суммировании, как известно, возрастает пропорционально корню квадратному из числа слагаемых. [Каценбоген М.С. Характеристики обнаружения. - М.: Советское Радио. 1965], /4/:

где а энергия полезного сигнала, а е - энергия случайных шумов и помех. Таким образом, при больших mn величина второго слагаемого e√(mn) сравнительно мало меняется в зависимости от mn. Исходя из этого разделяют полную сумму последовательности единичных отсчетов Σ(A_mn) на части, например начальную Σ(A_mn/r) и конечную {Σ(A_mn)-Σ(A_mn/r)}.

Величины сигнала и случайной помехи по предложенному способу находят из предположения равной вероятности возникновения помехи как на переднем отрезке от 1 до mn/r, так и на заднем отрезке от mn/r до mn, где r - делитель mn, решая следующую систему уравнений с двумя неизвестными а и е:

Вычисления неизвестных а и е таким образом получают экстраполяцией начальной промежуточной информации в виде Σ(A_mn/r), на полный интервал Σ(A_mn). Этого часто не достаточно для достоверного определения тотальной амплитуды ТА на временном разрезе, так как на сейсмограммах, кроме заданных регулярных отраженных волн, всегда присутствуют регулярные помехи, например, в виде годографов кратных, дифрагированных, рассеянных и иных типов волн. Поэтому в предлагаемом способе используют дополнительный подход в виде обратной экстраполяции конечного интервала отсчетов от (mn-mn/r) до mn на полный (в том числе и на начальный от 1 до mn/r) и дополнительно решают систему следующих уравнений:

После вычисления неизвестных сравнивают их между собой и получают информацию об уровне помех отдельно на падающих и восстающих лучах, в первую очередь во взаимных точках, при удовлетворительном результате сравнения ТА документируют на тотальном временном разрезе, а если нет, то находят среднеарифметическую величину ТА на временном разрезе:

после чего сравнивают полученный результат со среднеарифметическим Σ(A_mn)/mn и получают информацию о наличии регулярной помехи на протяжении всего эйконала волны заданного типа. Результаты, полученные при определении величин а и также можно представить в виде матриц a(α, V) и и выявлять области распределения случайных и регулярных волн-помех сложением, вычитанием, умножением элементов этих матриц.

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

Построение изображений сложно построенных сред в виде динамического глубинного K(x, h) и/или тотального (совокупного) временного разреза T(x, t) осуществляется трансформацией результатов обзора волновых полей на сейсмограммах A_n(x, t) на основе использования в процессе построения дополнительных предварительных версий глубинного разреза A(x, h). Сущность способа состоит в том, что сначала для каждой отображаемой точки OT(x₀, h₀) предварительного глубинного разреза A(x, h) осуществляют обзор волнового поля A(x, t) на n сейсмограммах, зарегистрированных по способу многократных перекрытий в результате засветки под различными углами отображаемых точек (ОТ) из n п