Способ сейсморазведки для определения положения низкоскоростных и высокоскоростных геологических тел вдоль глубинного разреза в районах, осложненных рельефом Российский патент 2025 года по МПК G01V1/28 G01V1/30 

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

Известен способ сейсморазведки, включающий регистрацию естественного сейсмического фона по трем компонентам не менее чем двумя сейсмоприемниками до и после генерирования сейсмических колебаний с частотой 0,1-70 Гц, суждение о наличии нефтегазовых месторождений по увеличению площади под кривой спектра сейсмического фона по всем трем компонентам после генерирования сейсмических колебаний по сравнению с исходным (RU №2119677, G01V 1/00, 1998).

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

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

где К - экспериментально установленный числовой коэффициент, зависящий от слагающих пород;

λ - длина волны Рэлея, определенная по дисперсионной кривой, м. (см. патент РФ №2271554, G01V 1/00, 2005).

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

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

Технический результат достигается в способе сейсморазведки для определения положения низкоскоростных и высокоскоростных геологических тел вдоль глубинного разреза в районах, осложненных рельефом, включающем установку на исследуемом профиле X длиной L_X по меньшей мере двух вертикальных с идентичными амплитудно-частотными характеристиками сейсмометров, один из которых Z₀ установлен стационарно, а остальные - с возможностью перемещения и установкой в каждой точке измерений j с координатой x_j с шагом Δх, не превышающим минимальную глубину исследования H_min, синхронную регистрацию микросейсмических сигналов, состоящих из волн Рэлея, сейсмометрами Z₀ и Z_j в каждой точке j в течение промежутков времени регистрации T_j, продолжительность которых определяется периодом стационарности микросейсмического сигнала, вычисление спектров мощности микросейсмического сигнала сейсмометров Z_j в каждой точке измерения j, где f_k - частота в спектре, и спектров мощности S_0,j(f_k) сигналов стационарного сейсмометра Z₀ для каждого промежутка времени регистрации T_j с последующим вычислением относительного спектра мощности в каждой точке измерений j с координатой x_j и для каждой выбранной частоты f_k путем нахождения отношений полученных спектров мощности и соответствующих им по промежутку времени T_j спектров мощности по формуле вычисление среднего по всем точкам j относительного спектра мощности по формуле где n - количество точек j регистрации микросейсмического сигнала сейсмометрами Z_j, с последующим вычислением отклонения спектров мощности от среднего относительного спектра мощности по формуле определение длин волн λ_k поверхностной волны Рэлея для каждой выбранной частоты f_k в спектре по формуле где v(f_k) - средняя фазовая скорость волны Рэлея, измерение высоты уровня рельефа H_i в точках i с координатами x_i вдоль профиля X и на продолжениях профиля X от крайних точек j в обе стороны длиной L_П, которая определяется как минимальное значение между L_X и 3H_max по формуле где H_max - максимальная глубина исследования, при этом шаг Δ между точками i на профиле X составляет Δ=Δх/r₁, а на продолжениях профиля где r₁ и r₂ - числовые коэффициенты, характеризующие плотность сетки измерения высоты уровня рельефа H_i, выделение из всех точек i для каждой точки j и для каждой длины волны λ_k множества точек i, координаты которых попадают в сегмент профиля, определяемого условием

где L_k - длина сегмента усреднения уровня рельефа, определяемая из соотношения где k_L - числовой коэффициент, оцененный из численного моделирования, и равный λ_k - длина волны Рэлея, f_k - частота в спектре, вычисление в каждой точке j и для каждой длины волны λ_k средней отметки уровня рельефа как среднее взвешенное значение H_i с весом пи по всем точкам i из множества где вес mi значения H_i для крайних точек i равен а для всех остальных точек i определяется формулой построение глубинного разреза вдоль профиля X путем привязки вычисленных значений отклонения относительных спектров мощности к точкам в пространстве под каждой точкой измерения j с координатой x_j и положением по глубине H_k согласно соотношению где k_Г - экспериментально установленный числовой коэффициент, зависящий от слагающих пород, определение положения вдоль профиля х и по глубине Н низкоскоростных геологических тел на глубинном разрезе путем выделения областей с повышенным положительным значением W(x, Н) и определение высокоскоростных геологических тел по областям с низким отрицательным значением W(x, Н) на глубинном разрезе.

Измерение высоты уровня рельефа H_i в точках i с координатами x_i вдоль профиля X и на продолжениях профиля X от крайних точек j в обе стороны длиной L_П и выбор длины продолжений профиля X и шагов А между точками i на профиле X и на продолжениях профиля X позволяет получить фактические значения высот H_i в достаточном для учета рельефа множестве точек i, включая получение оценок средних отметок уровня рельефа для всех точек j измерения микросейсмического сигнала и для каждой длины волны λ_k.

Выделение из всех точек i для каждой точки j и для каждой длины волны λ_k множества точек i, координаты которых попадают в сегмент профиля позволяет отобрать для вычисления средних отметок уровней рельефа только те значения H_i, которые в данной j точке и для данной длины волны λ_k находятся в пределах сегмента усреднения пропорционального длине волны, что учитывает особенности распространения поверхностных волн.

Построение глубинного разреза вдоль профиля X путем привязки вычисленных значений отклонения относительных спектров мощности к точкам в пространстве под каждой точкой измерения j с координатой x_j и положением по глубине H_k согласно соотношению позволяет получить изображение геологической среды в координатах х и Н с учетом фактического рельефа и с учетом различия в одной и той же точке средних уровней рельефа для разных длин волн, что учитывает особенности распространения поверхностной волны Рэлея, повышая точность и достоверность результата.

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

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

Устанавливают на исследуемом профиле X длиной L_X по меньшей мере два вертикальных с идентичными амплитудно-частотными характеристиками сейсмометра, один из которых Z₀ установлен стационарно, а остальные - с возможностью перемещения и установкой в каждой точке измерений j с координатой x_j с шагом Δх, не превышающим минимальную глубину исследования H_min. Проводят синхронную регистрацию микросейсмических сигналов, состоящих из волн Рэлея, сейсмометрами Z₀ и Z_j в каждой точке j в течение промежутков времени регистрации T_j, продолжительность которых определяется периодом стационарности микросейсмического сигнала. Вычисляют спектры мощности микросейсмического сигнала сейсмометров Z_j в каждой точке измерения j, где f_k - частота в спектре, и спектры мощности сигналов стационарного сейсмометра Z₀ для каждого промежутка времени регистрации T_j. После чего вычисляют относительный спектр мощности в каждой точке измерений j с координатой x_j и для каждой выбранной частоты f_k путем нахождения отношений полученных спектров мощности и соответствующих им по промежутку времени T_j спектров мощности по формуле Затем вычисляют средний по всем точкам j относительный спектр мощности по формуле где n - количество точек j регистрации микросейсмического сигнала сейсмометрами Z_j и отклонение спектров мощности от среднего относительного спектра мощности по формуле Определяют длины волн λ_k поверхностной волны Рэлея для каждой выбранной частоты f_k в спектре по формуле где v(f_k) - средняя фазовая скорость волны Рэлея. Измеряют высоты уровня рельефа H_i в точках i с координатами x_i вдоль профиля X и на продолжениях профиля X от крайних точек x_j в обе стороны длиной L_П, которую определяют как минимальное значение между по формуле где H_max - максимальная глубина исследования, при этом шаг Δ между точками i на профиле X составляет Δ = Δх / r₁, а на продолжениях профиля где r₁ и r₂ - числовые коэффициенты, характеризующие плотность сетки измерения высоты уровня рельефа H_i. Выделяют из всех точек i для каждой точки j и для каждой длины волны λ_k множества Ω_j,k точек i, координаты которых попадают в сегмент профиля, определяемого условием где L_k - длина сегмента усреднения уровня рельефа, определяемая из соотношения где k_L - числовой коэффициент, оцененный из численного моделирования, и равный λ_k - длина волны Рэлея, f_k - частота в спектре. Вычисляют в каждой точке j и для каждой длины волны λ_k средней отметки уровня рельефа как среднее взвешенное значение H_i с весом m_i по всем точкам i из множества Ω_j,k, где вес m_i значения H_i для крайних точек i равен а для всех остальных точек i определяется формулой Проводят построение глубинного разреза вдоль профиля X путем привязки вычисленных значений отклонения относительных спектров мощности к точкам в пространстве под каждой точкой измерения j с координатой x_j и положением по глубине H_k согласно соотношению где k_Г - экспериментально установленный числовой коэффициент, зависящий от слагающих пород. Определяют положения х вдоль профиля и по глубине Н низкоскоростных геологических тел на глубинном разрезе путем выделения областей с повышенным положительным значением W(x, Н) и высокоскоростных геологических тел по областям с низким отрицательным значением W(x, Н) на глубинном разрезе.

Конкретный пример.

Сейсморазведку проводили на территории африканского континента, длина профиля X составила L_X=1750 м. Регистрация микросейсмического сигнала производилась в 60 точках, включая стационарно установленный сейсмометр Z₀. Рабочий профиль содержал n=59 точек регистрации микросейсмического сигнала, расстояние между точками измерения микросейсмического сигнала приблизительно составляло Δх=30 м, таким образом, минимальная возможная глубина исследования также составляла примерно H_min=30 м. Максимальная интересуемая глубина исследования составляла H_max=1000 м. Использовалась локальная система координат вдоль профиля так, что точка j=1 имеет координату x_j=1=0, точка j=2-x_j=2=30 м и т.д. до точки j=59 x_j=59=Lx. Для сейсморазведки применялись 32 сейсмоприемника СМ-3КВ с цифровыми регистраторами "Scout" с идентичными амплитудно-частотными характеристиками, чувствительность измерительного канала 200 вольт/(метр/сек), частотный диапазон канала 0,5-100 Гц. Регистрировалась только вертикальная компонента микросейсмического сигнала. Стационарный сейсмометр Z₀ устанавливался вблизи 30-ой точки в стороне от профиля X, остальные 31 устанавливались сначала в точках 1-31 и регистрировали микросейсмические сигналы одновременно со стационарно установленной станцией в течение промежутков времени регистрации T₁ - Т₃₁ продолжительностью 3-4 часа, затем регистрация на всех приборах кроме стационарного Z₀ была остановлена, работа 3 приборов из точек 1-3 была завершена, а 28 приборов из точек 4-31 были переустановлены в точки 32-59, после чего они синхронно со стационарным прибором Z₀ регистрировали микросейсмический сигнал в течение промежутков времени Т₃₂-Т₅₉ продолжительностью 3-4 часа. После этого приборы в точках 32-59 выводились из регистрации, последней останавливалась регистрация стационарно установленным прибором Z₀. Рассчитывали спектры мощности вертикальной компоненты микросейсмического сигнала S_j(f_k) в каждой точке регистрации j = 1, 2, …, 59, а также 59 спектров мощности вертикальной компоненты микросейсмического сигнала S_0,j(f_k), зарегистрированного стационарно установленным прибором, вычисленным для промежутков времени T₁-Т₅₉, где f_k - частоты в спектре от 1 до 15 Гц с шагом ≈ 0.00333 Гц, всего 4200 частот f_k, k = 1, 2, …, 4200. Далее вычислялись 59 относительных спектров мощности w(x_j, f_k) по одному для каждой точки измерений j с координатой x_j путем нахождения отношений полученных спектров мощности и соответствующих им по промежутку времени T_j спектров мощности по формуле j = 1, 2, …, 59. После чего вычислялся средний по всем точкам j относительный спектр мощности по формуле где n=59 - количество точек j регистрации микросейсмического сигнала сейсмометрами Z_j, не считая стационарно установленного сейсмометра Z₀, с последующим вычислением 59 отклонений спектров мощности от среднего относительного спектра мощности по формуле j = 1, 2, …, 59. Далее для каждой из 4200 частот f_k в спектре от 1 до 15 Гц определялись длины волн λ_k поверхностной волны Рэлея по формуле где v(f_k) - средняя фазовая скорость волны Рэлея, значения которой составляли от , Измерения высоты уровня рельефа H_i проводились в 175 точках i с координатами x_i вдоль линии профиля X и на его продолжениях от крайних точек в обе стороны длиной L_П=L_X=1750 м которая, согласно формуле принимает минимальное значение между при этом шаг Δ между точками i на профиле X составлял а на продолжениях профиля где коэффициенты плотности сетки r₁=r₂=2, таким образом, координатная сетка x_i точек измерений высоты уровня рельефа H_i включала: (1) - 29 точек i = 1, …, 29 с координатами от с шагом 60 м на продолжении профиля «влево» от точки x_j=1, (2) - 117 точек i = 30, …, 146 с координатами от ; с шагом 15 м на линии профиля X между точек измерения микросейсмического сигнала и (3) - 29 точек i = 147, …, 175 с координатами от x_i=147=1810 м до x_i=175=3500 м с шагом 60 м на продолжении профиля «вправо» от точки x_j=59. Далее с использованием простого компьютерного алгоритма, включающего: (1) вычисление для каждой длины волны λ_k, где k = 1, 2, …, 4200, длины сегмента L_k усреднения уровня рельефа по формуле где k_L - числовой коэффициент принятый равным k_L = 3; (2) тройной цикл - внешний по точкам j от 1 до 59, вложенный в него - по номеру длины волны к от 1 до 4200 и самый внутренний - по номерам точек i от 1 до 175, внутри которого выбирались из всех точек i для каждой точки j и для каждой длины волны λ_k множества точки i, координаты которых попадают в сегмент профиля, определяемого неравенствами (3) вычисление весовых коэффициентов m_i значений высот H_i для крайних точек i = 1 и i = 175 по формуле а для всех остальных 173-ех точек i в цикле i = 2, 3, …, 174 по формуле (4) тройной цикл по j, k, i, включающий вычисление в каждой точке j от 1 до 59 и для каждой длины волны средней отметки уровня рельефа как среднее взвешенное значение H_i с весом пи по всем i, принадлежащим множеству определенному для тех же номеров j и k. Затем строили глубинный разрез вдоль профиля X (фигура) путем привязки вычисленных значений отклонения относительных спектров мощности к точкам в пространстве под каждой из 59 точек измерения микросейсмического сигнала j с координатой x_j и положением по глубине H_k, которое определялось формулой где коэффициент глубинной привязки кг был выбран, равным 0,4, определение положения низкоскоростных геологических тел на глубинном разрезе проводилось путем выделения областей с повышенным положительным значением W(x, Н) - белая штриховая линия на фигуре 1а, определение положения высокоскоростных геологических тел выполнялось по положению областей с низким отрицательным значением W(x, Н) - темные области на фигуре 1а. Использование предложенного способа сейсморазведки для определения положения низкоскоростных и высокоскоростных геологических тел вдоль глубинного разреза в районах позволяет определять положение и глубину геологических объектов в присутствии и с учетом рельефа без использования искусственных источников сейсмического сигнала, включая определение консолидированных зон и зон пониженной прочности земной коры с целью картирования глубинных неоднородностей Земли при поиске месторождений, в частности рудных, алмазных и углеводородных, геотермальных коллекторов, для оптимального выбора мест бурения разведочных и эксплуатационных скважин, для исследования районов и площадок размещения инженерных сооружений, для определения приповерхностных и глубинных разломных структур в районах, осложненных рельефом, значительно отличающимся от плоского.

Изобретение относится к области сейсморазведки и может быть использовано для поиска месторождений, в частности рудных, алмазных и углеводородных, для определения консолидированных зон и зон пониженной прочности в районах, осложненных рельефом, отличным от плоского. Согласно заявленному решению устанавливают на исследуемом профиле X длиной L_X по меньшей мере два вертикальных с идентичными амплитудно-частотными характеристиками сейсмометра, один из которых Z₀ установлен стационарно, а остальные - с возможностью перемещения и установкой в каждой точке измерений j с координатой x_j с шагом Δх, не превышающим минимальную глубину исследования H_min. Проводят синхронную регистрацию микросейсмических сигналов, состоящих из волн Рэлея, сейсмометрами Z₀ и Z_j в каждой точке j в течение промежутков времени регистрации T_j, продолжительность которых определяется периодом стационарности микросейсмического сигнала. Вычисляют спектры мощности микросейсмического сигнала S_j(f_k) сейсмометров Z_j в каждой точке измерения j, где f_k - частота в спектре, и спектры мощности S_0,j(f_k) сигналов стационарного сейсмометра Z₀ для каждого промежутка времени регистрации T_j. После чего вычисляют относительный спектр мощности w(x_j, f_k) в каждой точке измерений j с координатой x_j и для каждой выбранной частоты f_k путем нахождения отношений полученных спектров мощности S_j(f_k) и соответствующих им по промежутку времени T_j спектров мощности S_0,j(f_k) по формуле w(x_j, f_k) = S_j(f_k)/S_0,j(f_k). Затем вычисляют средний по всем точкам j относительный спектр мощности w(x_j, f_k) по формуле w₀(f_k)=1/n∑_jw(x_j, f_k), где n - количество точек j регистрации микросейсмического сигнала сейсмометрами Z_j и отклонение спектров мощности от среднего относительного спектра мощности по формуле W(x_j, f_k)=w(x_j, f_k)-w₀(f_k). Определяют длины волн λ_k поверхностной волны Рэлея для каждой выбранной частоты f_k в спектре. Выделяют из всех точек i для каждой точки j и для каждой длины волны λ_k множества Ω_j,k точек i, координаты которых попадают в сегмент профиля, определяемого условием x_j-L_k/2≤x_i≤x_j+L_k/2, где L_k - длина сегмента усреднения уровня рельефа, определяемая из соотношения L_k=k_L×λ_k, где k_L - числовой коэффициент, оцененный из численного моделирования и равный k_L≈2-4, λ_k - длина волны Рэлея, f_k - частота в спектре. Вычисляют в каждой точке j и для каждой длины волны λ_k средней отметки уровня рельефа H₀(λ_k, _xj) как среднее взвешенное значение Hi с весом mi по всем точкам i из множества Ω_j,k, где вес m_i значения H_i для крайних точек i равен m_i=2⋅r₂⋅Δх, а для всех остальных точек i определяется формулой m_i=(x_i+1-x_i-1)/2. Проводят построение глубинного разреза W(x_j, H_k) вдоль профиля X путем привязки вычисленных значений отклонения относительных спектров мощности W(x_j, f_k) к точкам в пространстве под каждой точкой измерения j с координатой x_j и положением по глубине H_k согласно соотношению H_k=H₀(λ_k, x_j)-k_Г×λ_k, где k_Г - экспериментально установленный числовой коэффициент, зависящий от слагающих пород. Определяют положения х вдоль профиля и по глубине Н низкоскоростных геологических тел на глубинном разрезе путем выделения областей с повышенным положительным значением W(x, Н) и высокоскоростных геологических тел по областям с низким отрицательным значением W(x, Н) на глубинном разрезе. Технический результат – повышение информативности и достоверности получаемых данных. 1 ил.

Способ сейсморазведки для определения положения низкоскоростных и высокоскоростных геологических тел вдоль глубинного разреза в районах, осложненных рельефом, включающий установку на исследуемом профиле X длиной L_X по меньшей мере двух вертикальных с идентичными амплитудно-частотными характеристиками сейсмометров, один из которых Z₀ установлен стационарно, а остальные - с возможностью перемещения и установкой в каждой точке измерений j с координатой x_j с шагом Δх, не превышающим минимальную глубину исследования H_min, синхронную регистрацию микросейсмических сигналов, состоящих из волн Рэлея, сейсмометрами Z₀ и Z_j в каждой точке j в течение промежутков времени регистрации T_j, продолжительность которых определяется периодом стационарности микросейсмического сигнала, вычисление спектров мощности микросейсмического сигнала S_j(f_k) сейсмометров Z_j в каждой точке измерения j, где f_k - частота в спектре, и спектров мощности S_0j(f_k) сигналов стационарного сейсмометра Z₀ для каждого промежутка времени регистрации T_j с последующим вычислением относительного спектра мощности w(x_j, f_k) в каждой точке измерений j с координатой X_j и для каждой выбранной частоты f_k путем нахождения отношений полученных спектров мощности S_j(f_k) и соответствующих им по промежутку времени T_j спектров мощности S_0j(f_k) по формуле w(x_j, f_k)=S_j(f_k)/S_0j(f_k), вычисление среднего по всем точкам j относительного спектра мощности w(x_j, f_k) по формуле w₀(f_k)=1/nΣ_jw(x_j, f_k), где n - количество точек j регистрации микросейсмического сигнала сейсмометрами Z_j, с последующим вычислением отклонения спектров мощности от среднего относительного спектра мощности по формуле W(x_j, f_k)=w(x_j, f_k)-w₀(f_k), определение длин волн λ_k поверхностной волны Рэлея для каждой выбранной частоты f_k в спектре по формуле λ_k=v(f_k)/f_k, где v(f_k) - средняя фазовая скорость волны Рэлея, измерение высоты уровня рельефа H_i в точках i с координатами x_i вдоль профиля X и на продолжениях профиля X от крайних точек x_j в обе стороны длиной L_П, которая определяется как минимальное значение между L_X и 3H_max по формуле L_П=min{L_X, 3H_max}, где H_max - максимальная глубина исследования, при этом шаг Δ между точками i на профиле X составляет Δ=Δх/r₁, а на продолжениях профиля X-Δ=r2⋅Δх, где r₁ и r₂ - числовые коэффициенты, характеризующие плотность сетки измерения высоты уровня рельефа H_i, выделение из всех точек i для каждой точки j и для каждой длины волны λ_k множества Ω_j,k точек i, координаты которых попадают в сегмент профиля, определяемого условием

x_j-L_k/2≤x_i≤x_j+L_k/2,

где L_k - длина сегмента усреднения уровня рельефа, определяемая из соотношения L_k=k_L×λ_k, где k_L - числовой коэффициент, оцененный из численного моделирования, и равный k_L≈2-4, λ_k - длина волны Рэлея, f_k - частота в спектре, вычисление в каждой точке j и для каждой длины волны λ_k средней отметки уровня рельефа H₀(λ_k, x_j) как среднее взвешенное значение H_i с весом m_i по всем точкам i из множества Ω_j,k, где вес m_i значения H_i для крайних точек i равен m_i=2⋅r₂⋅Δх, а для всех остальных точек i определяется формулой m_i=(x_i+1-x_j-1)/2, построение глубинного разреза W(x_j, H_k) вдоль профиля X путем привязки вычисленных значений отклонения относительных спектров мощности W(x_j, f_k) к точкам в пространстве под каждой точкой измерения j с координатой x_j и положением по глубине H_k согласно соотношению H_k=H₀(λ_k, x_j)-k_г×λ_k,

где k_Г - экспериментально установленный числовой коэффициент, зависящий от слагающих пород, определение положения х вдоль профиля и по глубине Н низкоскоростных геологических тел на глубинном разрезе путем выделения областей с положительным значением W(x, Н), соответствующих области максимального положительного значения глубины Н на глубинном разрезе, определение высокоскоростных геологических тел по областям с отрицательным значением W(x, Н) на глубинном разрезе, соответствующих области минимального отрицательного значения глубины Н на глубинном разрезе.