СПОСОБ ПОЛУЧЕНИЯ АПРИОРНОГО ГОДОГРАФА ДЛЯ ВЫПОЛНЕНИЯ ЛИТОЛОГО-СТРАТИГРАФИЧЕСКОЙ ПРИВЯЗКИ Российский патент 2015 года по МПК G01V11/00 G01V1/00 G01V1/28 

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

Известен способ построения сейсмического глубинного разреза, включающий развертку отображаемых точек (ОТ) разреза A(x,h) с заданным шагом по столбцам (h) и/или строкам (х), из каждой ОТ (х₀,h₀) производят обзор волнового поля для поиска амплитуд сигналов А сейсмических границ, возможно принадлежащих реперным и/или целевым горизонтам, на n сейсмограммах, зарегистрированных способом многократного профилирования, обеспечивающего засветку под различными углами отображаемых точек из n пунктов воздействий, предполагая наличие сейсмической границы в отображаемых точках; обзор реализуют серией решений квазипрямой задачи лучевого трассирования, находят наиболее вероятный эйконал, используя серию виртуальных годографов для получения маркерных отметок времени на сейсмотрассах; генерацию годографов производят в соответствии с предварительно заданным множеством скоростей V=V₀±ΔV(1,2,…,υ/2,…,υ) (где V₀ - априорно заданное первоначальное значение скорости из заданного массива скоростей V(x,h), х, h - координаты глубинного разреза, ΔV шаг приращения скорости, υ - число шаговых приращений скорости), а также с предварительно заданной определенной конфигурацией сейсмического объекта, который задан q ракурсами в отображаемой точке с определенным шагом поворота ΔV угла падения (с центром вращения в отображаемой точке, для каждого ракурса объекта задают множество углов падения α=±α₀±Δα(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); α_0i=α_0(i-1)±δα=α_0(i-1)Δα(q/2-q_max); (1),

где индексами i и (i-1) обозначены последующие и предыдущие априорные и текущие начальные значения диапазона изменения скорости и угла падения в отображаемой точке; последовательно в блоке F^max(υ_max,q_max) заменяют в автоматическом режиме предыдущие значения диапазона (на скорректированные, а монотонно возрастающие с глубиной значения скорости заменяют лишь при условии совпадения знаков развертки отображаемой точки по столбцам со знаком приращения δV скорости V, тем самым размещают F(q_(i-1), υ_(i-1)) в центр новой скорректированной сетки значений матрицы окна обзора достоверности, для последующих q×υ углов падения и скоростей находят скорректированные значения F(q_(i-1),υ_i), среди них в пределах следующего окна обзора достоверности находят новое текущее значение F^max, документируют очередную отображаемую точку с интенсивностью, обусловленной найденными параметрами F^max, υ_max и q_max, размещают F^max вновь в центр очередной скорректированной сетки окна обзора достоверности, и так далее до полного построения глубинного сейсмического разреза в соответствии с заданным режимом развертки отображаемых точек по строкам и/или столбцам (См. патент РФ №2463626, МПК G01V 1/28, опубл. 10.10.2012).

Известен также (См. патент №2449322, МПК G01V 1/28, опубл. 27.04.12) способ построения сейсмического глубинного разреза, включающий развертку отображаемых точек разреза по столбцам (h) и/или строкам (х), из каждой отображаемой точки ОТ (х₀,h₀) производят обзор волнового поля на сейсмограммах для поиска амплитуд сигналов А сейсмических границ, возможно принадлежащих реперным и/или целевым горизонтам; обзор реализуют серией решений квазипрямой задачи лучевого трассирования с использованием виртуальных годографов, для чего задают определенную конфигурацию сейсмического объекта и при каждом облучении отображаемой точки из пункта воздействия задают q ракурсов сейсмического объекта с определенным шагом поворота Δα угла наклона α_C центром вращения в отображаемой точке, для каждого ракурса объекта α=±α₀±Δα(1,2,…,q/2,…,q), где α₀ - априорно заданное первоначальное значение угла наклона границы в отображаемой точке, рассчитывают виртуальный годограф заданного типа сейсмической волны с учетом априорно заданной скорости, координат отображаемой точки и заданного закона (функции) расположения последовательности из m точек, прилегающих к отображаемой точке, производят направленное синфазное накопление амплитуд многокомпонентных колебаний на зарегистрированных m сейсмотрассах вдоль маркерных отметок времени на годографах, причем для каждого колебания считывают ρ компонентов (ρ≤?20) в виде псевдосинфазных отсчетов величин амплитуд; после синфазного накопления вдоль q годографов по ρ_m амплитуд колебаний сигналов из каждой из n сейсмограмм сигналы группируют по признаку равенства номера шага позиционирования (ракурса) q, в каждой из q групп производят супернакопления из ρ_mn амплитуд сигналов после введения поправочных коэффициентов на угол засветки, для каждого из q супернакоплений амплитуд по ρ фаз каждого из mn сигналов, специфичных для заданных позиций объекта в отображаемой точке, вычисляют q функционалов статистической обработки, сравнивают их между собой, по максимальной величине одного из функционалов определяют наиболее вероятный номер q позиции объекта, причем после определения наиболее вероятного номера q позиции объекта для каждой отображаемой точки строят двумерную сетку значений функционалов статистической обработки типа энергий супернакоплений, W(α,V) в зависимости от априорных углов наклона α и априорно заданного множества скоростей V при V=V₀±ΔV(1,2,…,υ/2,…,υ), где V₀ - априорно заданное первоначальное значение скорости из заданного массива скоростей V(x,h), х, h - координаты глубинного разреза, ΔV - шаг приращения скорости, υ - число шаговых приращений скорости; документируют все q·υ значений энергий цугов колебаний W(α,V) в узлах прямоугольной сетки в виде двумерной матрицы, по столбцам которой документируют значения энергии в зависимости от угла наклона в отображаемой точке, а по строкам - значения энергии в зависимости от скорости (или наоборот); строят двумерное отображение графика W(α,V), график представляют в виде окна обзора достоверности, что обеспечивает визуальный контроль за процессом опознания сейсмических границ и построением сейсмического глубинного разреза; в пределах графика окна обзора достоверности определяют максимум W^max (υ_max,q_max), его сравнивают с заданным порогом обнаружения, при превышении порога обнаружения отображаемую точку наносят на глубинный сейсмический разрез с интенсивностью, пропорциональной величине максимального функционала в пределах окна обзора достоверности; а затем для построения последующей отображаемой точки на глубинном разрезе определяют скорректированные значения углов наклона и скоростей по смещению W^max относительно центра окна обзора достоверности, где находят W(υ/2,q/2) и сравнивают его аргументы υ/2 и q/2 с υ_max и q_max; величину и знак скорректированных значений для V₀ и α₀ определяют по разностям: V_0i=V_0(i-1)±δV=V_0(i-1)±V(υ/2-υ_max); α_0i=α_0(i-1)±δα=α_0(i-1)±α(υ/2-D_max), где индексами i и (i-1) обозначены последующие и предыдущие априорные и текущие начальные значения диапазона изменения скорости и угла наклона в отображаемой точке; в автоматическом режиме последовательно заменяют предыдущие значения диапазона (α,V) на скорректированные, тем самым размещают W^max (q_(i-1),υ_(i-1)) в центр новой скорректированной сетки значений матрицы окна обзора достоверности, снова для всех q·υ находят скорректированные значения W(q_i,υ_i), среди них в пределах окна обзора достоверности находят текущее значение W^max, документируют очередную отображаемую точку с интенсивностью, соответствующей величине W^max, размещают W^max вновь в центр очередной скорректированной сетки окна обзора достоверности и повторяют процесс нанесения отображаемых точек до полного построения глубинного сейсмического разреза в соответствии с заданным режимом развертки отображаемых точек по строкам и/или столбцам.

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

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

Известен также способ построения сейсмического глубинного разреза (См. патент РФ №2463628, МПК G01V 1/28, опубл. 10.10.2012), включающий развертку отображаемых точек (ОТ) разреза A(x,h) с заданным шагом по столбцам (h) и/или строкам (х), из каждой ОТ (х₀,h₀) производят обзор волнового поля для поиска амплитуд сигналов А сейсмических границ, возможно принадлежащих реперным и/или целевым горизонтам, на n сейсмограммах, зарегистрированных способом многократного профилирования, обеспечивающего засветку под различными углами отображаемых точек из n пунктов воздействий, предполагая наличие сейсмической границы в отображаемых точках; обзор реализуют серией решений квазипрямой задачи лучевого трассирования, находят наиболее вероятный эйконал, используя серию виртуальных годографов для получения маркерных отметок времени на сейсмотрассах; генерацию годографов производят в соответствии с предварительно заданным множеством скоростей V=V₀±ΔV(1,2,…, υ₂,…, υ), (где V₀ - априорно заданное первоначальное значение скорости из заданного массива скоростей V(x,h), х, 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(v_max,q_max), в котором также хранят соответствующее максимальному уровню аргументы υ_max, и q_max, изображают отображаемую точку на глубинном разрезе с интенсивностью, обусловленной найденными параметрами F^max, v_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); α_0i=α_0(i-1)±δα=α_0(i-1)±Δα(q/2-q_max); (1),

где индексами 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 вновь в центр очередной скорректированной сетки окна обзора достоверности, и так далее до полного построения глубинного сейсмического разреза в соответствии с заданным режимом развертки отображаемых точек по строкам и/или столбцам.

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

Данные недостатки обусловлены тем, что достоверность построения глубинного разреза за счет влияния сейсмического сноса полностью относится к положению объекта в плане (координатах «X» и «У»), чего нельзя сказать о вертикальной составляющей глубины «Н», т.к. используются виртуальные годографы, модельные данные скоростей, а не реальные скорости среды. Критерии оценки используют только динамическую характеристику колебательных процессов - матрицы энергии функционалов F^W ρ_m(α,V), а кинематическая характеристика в качестве критерия не используется. Кроме всего этого массив «т» и «n» данных требует немало места в оперативной памяти ЭВМ, поэтому обработка «n» сейсмограмм по «m» сейсмотрасс требует значительных затрат машинного времени, причем контроль интерпретатора за процессом опознания сейсмических границ и построением сейсмического глубинного разреза требует выезд интерпретатора в поле и больших ресурсных затрат. Следует заметить, что метод, аналогичный методу регулированного направленного приема (МРНП), связан с хорошим приемом высокочастотных сейсмических волн, не дает хорошего результата, поскольку теряет низкочастотную часть спектра, что, в итоге, сказывается на разрешающей способности и обусловливает ограничение известного способа. Данный недостаток обусловлен методикой производства работ метода регулированного направленного приема (МРНП), который является аналогом описанного способа. В МРНП “волна, у которой годограф не совпадает с линией суммирования, попадает в область подавления характеристики направленности. Это проявляется тем сильнее, чем выше ее частотный спектр. Поэтому МРНП по отношению к высокочастотным сейсмическим волнам обладает более высокой разрешающей способностью”. Бондарев В.И., Крылатков С.М. Анализ данных сейсморазведки: Учебное пособие для студентов вузов. Екатеринбург. Издательство УГГГА, 2002. - 57 с.

Известен принятый за прототип способ построения синтетической сейсмограммы для одномерного моделирования (см. Гогоненков Г.Н. «Расчет и применение синтетических сейсмограмм». М., Недра, 1972, стр.46-49, 77-85), при котором «калибруют» значения акустической диаграммы метода акустического каротажа (T_ак=f(Н_ак)) по исходному годографу сейсмокаротажа (СК) или вертикального сейсмического профилирования (ВСП) (Т_{ск_исх}=f(Н_{ск_исх})). Затем рассчитывают последовательность коэффициентов отражения (К_ОТР.) по данным метода акустического каротажа (АК). После этого производят расчет импульсной синтетической сейсмограммы, полученной посредством свертки последовательности коэффициентов отражения с исходным сигналом. Далее последовательно каждой точке глубины «Н» исходного годографа сейсмокаротажа или вертикального сейсмического профилирования находят соответствующее значение «Т» временного разреза МОГТ (ищут новую зависимость Т_{СК_ПРИВ}=f(Н_{СК_ПРИВ})). Контроль качества одномерного моделирования (привязки) осуществляют с помощью функции взаимной корреляции (ФВК) ближайшей к скважине трассы ОГТ временного разреза и рассчитанной синтетической сейсмограммы, добиваясь максимального подобия экстремумов отражений от опорных и целевых горизонтов синтетической сейсмограммы экстремумам отражений от соответствующих опорных и целевых горизонтов временного разреза МОГТ.

Недостатком известного способа является высокая трудоемкость и низкие производительность (скорость получения результата), точность и достоверность.

Данный недостаток обусловлен недостаточной достоверностью, точностью, степенью подобия соответствия экстремумов от реперных и целевых геологических горизонтов синтетической сейсмограммы экстремумам соответствующих конкретных отраженных волн временного разреза, причем в качестве критерия оценки качества используется только динамический критерий колебательных процессов - функция взаимной корреляции (ФВК). Не используются кинематические критерии - скорости, что приводит к неточности выполнения привязки, необоснованной «растяжке» или «сжатии» годографа, привязке к не той фазе временного разреза, образованию неоправданных «ступеней» на графике зависимости интервальной скорости акустического каротажа от глубины, чего не бывает на качественных исходных диаграммах акустического каротажа, особенно, если привязка выполняется только в целевых интервалах, поскольку интервальная скорость с глубиной меняется по определенному закону. Кроме того, при выполнении привязки в целевом интервале, при отсутствии контроля диапазона интервала опорными горизонтами, не всегда есть возможность правильно соотнести значения скорости привязочного годографа V_{ИНТ-ПРИВ} и тренда скорости V_CK по данным СК и ВСП в точке скважины. Следует также заметить, что в процессе получения априорного годографа для выполнения литолого-стратиграфической привязки идет последовательный подбор значений точек годографа (глубина-время), причем в годографах ВСП количество точек годографа может достигать 200 и более, последовательный подбор значений в таком случае занимает много времени, т.к. очень трудоемок.

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

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

dTi=Нi∗dT_н/Н_н, где

dTi - поправка, вводимая в каждую i-тую точку последнего годографа;

dT_н-=tн^ОГТ-tн^СК - поправка, которую нужно ввести в последний годограф для посадки на нижний опорный горизонт, т.е. совмещения во временном масштабе значения геологической отметки нижнего опорного горизонта и значения времени экстремума соответствующей отраженной волны;

tн^ОГТ - значение времени, снятое с временного разреза метода общей глубинной точки (в точке скважины) для экстремума отраженной волны, соответствующей нижнему опорному горизонту;

tн^СК - значение времени, снятое с последнего годографа сейсмокаротажа, вертикального сейсмического профилирования для геологической отметки нижнего опорного горизонта;

Hi - значение геологической отметки i-й точки последнего годографа сейсмокаротажа, вертикального сейсмического профилирования;

Н_н - значение геологической отметки нижнего опорного горизонта;

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

Между отличительными признаками и достигаемым техническим результатом существует следующая причинно-следственная связь.

В отличие от аналога и прототипа использование в предлагаемом изобретении «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» совокупности признаков в виде того, что он включает в себя следующий ряд последовательных действий, как то получение данных сейсмокаротажа, акустического каротажа, плотностного гамма-гамма каротажа, калибровки данных значений акустического каротажа, перевода во временной масштаб, расчета последовательности коэффициентов отражения, свертки их с исходным сигналом, т.е построения синтетической сейсмограммы и кривой акустической жесткости с повторением циклов, при которых независимо получают и проводят подготовку данных методов общей глубинной точки, сейсмокаротажа, вертикального сейсмического профилирования, акустического каротажа и плотностного гамма-гамма каротажа, значений геологических отметок и проводят проверку качества полученных выше данных, при этом получают глубины опорных и целевых сейсмических горизонтов, которые соответствуют времени регистрации первых вступлений проходящей волны от опорных и целевых горизонтов, и эталонные интервальные скорости по данным сейсмокаротажа и вертикального сейсмического профилирования между этими опорными и целевыми горизонтами в координатах 2Д или 3Д, затем определяют исходный годограф, после чего рассчитывают синтетическую сейсмограмму, затем проводят контроль качества, вводят постоянную временную поправку для посадки на верхний опорный горизонт литолого-стратиграфического комплекса, например ОГ «Г» и рассчитывают синтетическую сейсмограмму и приводят контроль качества, вслед за этим вводят поправку для посадки на нижний опорный горизонт литолого-стратиграфического комплекса, например, ОГ «Б» и вновь рассчитывают синтетическую сейсмограмму, причем переносят точки полученного годографа на ближайшие акустически слабые границы и повторно рассчитывают синтетическую сейсмограмму. В совокупности признаков все это значительно повышает эффективность заявленного способа, в сравнении с объектом-прототипом, поскольку дает возможность получать окончательный годограф намного быстрее, чем в объекте-прототипе, где использовали повторение циклов (подбором Н и V) для каждой точки годографа, поскольку в заявляемом способе обе предложенные поправки вводят в точки преобразуемого годографа одновременно, рассчитывается значение времени для точки забоя с использованием значений эталонных интервальных скоростей, что увеличивает скорость получения априорного годографа и повышает эффективность способа за счет использования готовых данных методов общей глубинной точки (МОГТ), сейсмокаротажа (СК), вертикального сейсмического профилирования (ВСП), акустического каротажа (АК), гамма-гамма каротажа плотностного (ГГК-П), значений геологических отметок с исключением полевого этапа, кроме того, выполняется проверка качества и подготовка методов МОГТ, СК, ВСП, АК, ГГК-П, значений геологических отметок, а также возможность использования в качестве исходного годографа данных, используемых для получения эталонных значений интервальных скоростей, т.е. годографа, полученного по данным СК и ВСП в координатах 2Д или 3Д; осредненного годографа по площади; годографа по данным АК, что значительно снижает себестоимость заявленного способа с учетом временной поправки dTi и подготовки данных - редакции годографов СК или ВСП с учетом проверки их качества; получения суммарного временного разреза МОГТ с сохранением относительных амплитуд, приведением его к нульфазовому импульсу; моделирование кривых АК и ГГК-П и возможности использования в качестве исходного годографа данных, используемых для получения эталонных значений интервальных скоростей, т.е. глубин опорных и целевых горизонтов, соответствующих времен регистрации проходящих волн, полученных по данным СК и ВСП в координатах 2Д или 3Д; осредненного годографа по площади; годографа по данным АК. Повышение точности и достоверности результатов в заявленном способе получения априорного годографа для выполнения литолого-стратиграфической привязки, а именно повышение точности соответствия значений времен экстремумов отражений от опорных и целевых горизонтов временного разреза значениям времен геологических отметок скважины во временном масштабе, определяет также временная поправка dTi, получение эталонных значений интервальных скоростей и проверка качества, а также подготовка методов МОГТ, СК, ВСП, АК, ГГК-П, значений геологических отметок, посадка на верхний опорный горизонт литолого-стратиграфического комплекса и перенос точек годографа на ближайшие глинистые пачки, и, кроме того, построение синтетической сейсмограммы в процессе получения годографа, после выполнения совокупности определенных процедур. Следует также заметить, что соответствие, степени подобия экстремумов отражений от опорных и целевых горизонтов синтетической сейсмограммы экстремумам отражений от опорных и целевых горизонтов разреза МОГТ достигается путем использования временного разреза МОГТ с сохранением относительных амплитуд, приведением его к нульфазовому импульсу; применением помимо динамического критерия оценки колебательных процессов - функции взаимной корреляции (ФВК), критерия кинематического - значения интервальной скорости. Это достигается путем: сравнения значений интервальной скорости полученного годографа после привязки V_{ИНТ-ПРИВ} с определенными ранее по данным СК и ВСП эталонными значениями интервальных скоростей, характерными для определенного интервала и геологического района, все это, как дополнительный технический результат, определяет возможность контроля качества интервальной скорости годографа. Кроме того, предложенное техническое решение «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки», в совокупности заявленных признаков, обеспечивает надежный, достоверный результат, поскольку проводится проверка качества и подготовка временного разреза МОГТ с сохранением относительных амплитуд и приведением его к нульфазовому импульсу, помимо ФВК используются другие критерии оценки, а именно данные о распределении эталонных интервальных скоростей между основными реперными и целевыми горизонтами в данном районе и данном интервале, полученные ранее по данным проведенных сейсмокаротажа и вертикального сейсмического профилирования в точке скважины в координатах 2Д или 3Д. Следует заметить, что точность и достоверность результатов повышает также то, что используют данные о диапазонах значений интервальных скоростей, полученных лабораторным путем, т.е. измерением в соответствующем интервале на образцах горных пород в скважинах нужного геологического района; использованием осредненных значений АК в определенном интервале; построение синтетической сейсмограммы в процессе получения годографа, после выполнения определенных процедур. Точность результатов достигается также за счет равномерного внесения в каждую i-тую точку годографа поправки dTi, что позволяет получить значения интервальной скорости по данным АК без «ступеней», соответственно, тренду их постепенного увеличения с глубиной. Кроме всего прочего, дополнительным техническим результатом заявляемого изобретения «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» является то, что окончательный годограф может быть использован для пересчета временного разреза в глубинный разрез, а синтетическая сейсмограмма, кривая акустического импеданса могут быть использованы для выполнения инверсии, либо при использовании нескольких таких годографов и кривых от различных скважин площади тех же самых результатов (глубинного разреза и инверсии) в координатах 3Д. Заявляемый «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» позволяет также обоснованно уточнять геологические отметки в районах сложного строения (например, в районах с аномальным строением баженовской свиты), т.к. позволяет осуществить контроль по эталонным значениям интервальных скоростей в точке скважины, полученных в координатах 2Д или 3Д и характерных для соответствующего интервала и определенного геологического района, причем контроль также можно осуществлять, например, путем сравнения с диапазоном изменения значений интервальных скоростей, полученных лабораторным путем, посредством измерений в соответствующем интервале на образцах горных пород в скважинах данного геологического района, сравнением с осредненными значениями интервальных скоростей, полученных ранее по данным акустического каротажа на данном интервале.

Проведенный заявителем анализ уровня техники, включающий поиск по патентным и научно-техническим источникам информации и выявление источников, содержащих сведения об аналогах заявленного изобретения «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки», позволил установить, что заявитель не обнаружил источник, характеризующийся признаками, тождественными всем существенным признакам заявленного технического решения. По имеющимся у заявителя сведениям, совокупность существенных признаков заявляемого изобретения «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» не известна из уровня техники, что позволяет сделать вывод о соответствии заявляемого изобретения «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» критерию “новизна”. Определение из перечня выявленных аналогов прототипа, как наиболее близкого по совокупности признаков аналога, позволил выявить совокупность существенных, по отношению к усматриваемому заявителем техническому результату, отличительных признаков в заявляемом изобретении «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки», изложенных в формуле изобретения. Следовательно, заявляемое изобретение «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» соответствует критерию "новизна".

Для проверки соответствия заявляемого изобретения «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» критерию "изобретательский уровень" заявитель провел дополнительный поиск известных решений, чтобы выявить совокупность признаков, совпадающих с отличительными от прототипа признаками заявляемого изобретения «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки». Результаты поиска показали, что заявляемое изобретение «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» не вытекает для специалиста явным образом из известного уровня техники, поскольку из уровня техники, определенного заявителем, не выявлено влияние предусматриваемых существенными признаками заявленного изобретения «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» преобразований для достижения технического результата. Следовательно, заявленное изобретение «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» соответствует критерию "изобретательский уровень".

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

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

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

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

- фиг.2 - схема получения априорного годографа для выполнения экспресс-посадки.

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

ПРИМЕР 1.

1. Первоначально, согласно схеме на фиг.1, проводили сейсмокаротаж (СК) в ближайшей скважине. Получали продольный вертикальный годограф проходящей волны Tcк=f(Н), где Тск - время регистрации первых вступлений проходящей волны. Проверяли качество годографа СК, получив зависимость интервальной скорости как функцию глубины Vинт_ск=f(Н). Т.к. качество годографа удовлетворяло, брали его как исходный.

2. Затем проводили акустический каротаж (АК) в скважине, в которой выполняли литолого-стратиграфическую привязку. Записывали dt_AK=f(Н) - интервальное время пробега продольной волны. Проводили плотностной гамма-гамма каротаж (ГГК-П), записывали диаграмму изменений рассеянного гамма излучения как функцию глубины. Проверяли качество этих данных. Подготавливали диаграмму интервального времени акустического каротажа (АК) при этом:

- последовательно сшивали интервалы по глубине;

- редактировали значения интервального времени в интервалах образования каверн.

- оцифровывали диаграмму интервального времени пробега АК с шагом дискретизации 1 м.

Если АК в скважине не проводился, моделировали из других методов ГИС, например из диаграммы метода кажущегося сопротивления (КС). Далее пересчитывали интервальное время в интервальную скорость как Vинт_AK=f(H) и переходили от показаний ГГК-П к кривой плотности пород δп=f(H).

3. После этого получали значения отметок глубин геологических границ, к которым необходимо найти и привязать соответствующие сейсмические горизонты. Проверяли геологические отбивки по диаграммам геофизических исследований скважин (ГИС) и данным сейсморазведки.

4. Получали окончательный суммарный временной разрез метода общей глубинной точки (МОГТ) с сохранением относительных амплитуд, приводили его к нульфазовому импульсу и проводили корреляцию опорных (реперных) отражающих горизонтов (ОГ) «С», «Г», «М1», «М», «Б», «А».

5. После этого через значения глубин опорных и целевых сейсмических горизонтов, соответствующие времена регистрации первых вступлений проходящей волны от опорных и целевых горизонтов получали эталонные значения интервальных скоростей по данным СК и ВСП между этими горизонтами в координатах 2Д, рассчитав, например V_{инт_СК_}(Б-A)=2dH/dT, где dH - мощность интервала (в метрах), между значениями геологических отметок, соответствующих опорным ОГ «Б» и «А» по диаграммам ГИС; dT - временная мощность интервала (в секундах) между опорными ОГ «Б» и «А», которые служили критериями точности привязки.

6. Далее в качестве исходного сигнала брали импульс U(Tj) с временного разреза МОГТ из выбранного целевого интервала временного разреза в районе скважины и оценивали частотный диапазон и форму импульса.

7. Вслед за этим проводили коррекцию (калибровку) значений V_АК=f(Н) по сейсмокаротажу. Эта процедура необходима для исключения постоянной ошибки в данных АК и приведения данных АК и СК, в одни координаты. Полученную кривую интервального времени сохраняли. После этого получали диаграмму акустической жесткости γп=Vинт_AK*δп как результат произведения диаграммы интервальной скорости Vинт_АК=f(H) и диаграммы плотности пород δп=f(H).

8. После этого осуществляли перевод исходной информации из масштаба глубин во временной масштаб V¹(Н) в V(Т).

9. Затем рассчитывали и формировали массив коэффициентов отражения при нормальном падении волны на границу раздела упругих сред P₁(T) $$K=\frac{V\left(i+1\right)\delta п\left(i+1\right)-V\left(i\right)\delta п\left(i\right)}{V\left(i+1\right)\delta п\left(i+1\right)+V\left(i\right)\delta п\left(i\right)}$$ где

индекс i соответствует параметрам слоя, из которого падает волна, индекс i+1 - параметрам слоя, находящегося под границей.

10. Вслед за этим рассчитывали синтетическую сейсмограмму как результат свертки массива коэффициентов отражения P₁(T) с исходным сигналом U(Tj).

, где

q=T/0.002 - число отсчетов исходного сигнала,

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

11. Затем осуществляли посадку на верхний опорный горизонт литолого-стратиграфического комплекса, например (ОГ) Г введением постоянной поправки dТв, поскольку чаще всего данные АК и СК были получены в разных скважинах, то по этой причине, а также из-за разных подходов к определению интервальной скорости комплексируемых методов АК, СК, и МОГТ в полученной синтетической сейсмограмме, не достигается соответствия значений экстремумов отражений от первого опорного горизонта синтетической сейсмограммы значениям экстремумов соответствующего опорного горизонта временного разреза МОГТ, поэтому разницу dT измеряли на разрезе МОГТ и вводили в исходный годограф.

12. Вслед за этим вновь рассчитывали синтетическую сейсмограмму с исправленным годографом СК. Проверяли соответствие значений экстремумов отражений от верхнего опорного горизонта литолого-стратиграфического комплекса синтетической сейсмограммы и соответствующего опорного горизонта временного разреза МОГТ, наблюдали соответствие значений экстремумов. Далее осуществляли контроль качества годографа и синтетической сейсмограммы с помощью ФВК в ближнем интервале и сравнивали с эталонными значениями интервальных скоростей.

13. Следом осуществляли посадку на нижний опорный горизонт литолого-стратиграфического комплекса, например ОГ «Б». Такое же несоответствие значений экстремумов, о котором упоминалось выше в п.11, а именно: несоответствие значений экстремумов отражений от нижнего опорного горизонта литолого-стратиграфического комплекса синтетической сейсмограммы значениям экстремумов соответствующего опорного горизонта временного разреза МОГТ, наблюдали и в интервале ОГ «Б». Измеряли разницу dT_Б=tБ^ОГТ-tБ^СК на временном разрезе МОГТ.

Рассчитывали поправку $$dTi=Hi*\frac{dT6}{H6}$$ , где

dTi - поправка, вводимая в каждую i-тую точку последнего годографа;

dT_Б-=tБ^ОГТ-tБ^СК - поправка, которую нужно ввести в последний годограф для посадки на горизонт Б, т.е. совмещения во временном масштабе значения геологической отметки ОГ Б и значения времени экстремума соответствующей отраженной волны;

tБ^ОГТ - значение времени, снятое с временного разреза метода общей глубинной точки (в точке скважины) для экстремума отраженной волны, соответствующей горизонту Б;

t_БСК - значение времени, снятое с последнего годографа СК для геологической отметки горизонта Б;

Hi - значение геологической отметки i-й точки последнего годографа СК;

Н_Б - значение геологической отметки отражающего горизонта Б по диаграммам ГИС, уточненное по данным сейсморазведки.

Поправку Тиспр=Тпосл±dTi вводили в последний полученный годограф п.11.

14. Далее вновь рассчитывали синтетическую сейсмограмму с использованием исправленного годографа. Наблюдали соответствие значений экстремумов в интервале ОГ «Б». Осуществляли контроль качества годографа и синтетической сейсмограммы с помощью ФВК в соответствующем интервале и сравнением с эталонными значениями интервальных скоростей.

15. Вслед за этим переносили точки годографа на акустически слабые границы - глины.

16. После этого рассчитывали окончательную синтетическую сейсмограмму с использованием исправленного годографа и осуществляли контроль качества годографа и синтетической сейсмограммы с помощью ФВК в соответствующем интервале и сравнением с эталонными значениями интервальных скоростей. Затем сохраняли априорный годограф, синтетическую сейсмограмму, кривую акустической жесткости.

17. Следует заметить, что для скважин с глубинами забоя в самом нижнем литолого-стратиграфическом комплексе, например ниже ОГ Б, сначала фиксировали точку годографа, соответствующую ОГ «Б» (НБ - значение геологической отметки, ТБ - значение времени экстремума отраженной волны ОГ Б, снятое с временного разреза МОГТ).

18. Затем вводили точки годографа, соответствующие глинистым пачкам.

19. После этого рассчитывали значение времени для точки забоя скважины годографа как

Тзаб=Т_Б+dT, где dT=2dH/Vинт(Б-А), где

dH - расстояние в м от значения геологической отметки ОГ «Б» до значения глубины точки забоя;

Т_Б - значение времени, снятое с временного разреза МОГТ и соответствующее времени экстремума отраженной волны от ОГ «Б»;

Vинт (Б-А) - эталонная интервальная скорость по данным СК и ВСП, полученная в п.5 для самого нижнего литолого-стратиграфического комплекса, например интервала между ОГ «Б» и ОГ «А».

20. Вслед за этим рассчитывали синтетическую сейсмограмму по окончательному годографу СК и осуществляли контроль качества годографа и синтетической сейсмограммы с помощью ФВК в соответствующем интервале и сравнением с эталонными значениями интервальных скоростей.

21. После этого сохраняли априорный годограф, синтетическую сейсмограмму и кривую акустической жесткости.

Способ получения априорного годографа для выполнения экспресс-посадки осуществляли следующим образом.

ПРИМЕР 2.

1. Первоначально проводили сейсмокаротаж (СК) в ближайшей скважине и получали исходный продольный вертикальный годограф проходящей волны Tcк=f(Н), где Тск - время регистрации первых вступлений проходящей волны. Проверяли качество годографа СК, получив зависимость интервальной скорости, как функцию глубины Vинт_CK=f (Н). Т.к. качество годографа удовлетворяло, брали его как исходный.

2. Затем получали значения отметок глубин геологических границ, которым необходимо найти соответствующие сейсмические горизонты и проверяли геологические отбивки по диаграммам геофизических исследований скважин (ГИС) и данным сейсморазведки.

3. После этого получали окончательный суммарный временной разрез метода общей глубинной точки (МОГТ) с сохранением относительных амплитуд, приводили его к нульфазовому импульсу и проводили корреляцию опорных (реперных) отражающих горизонтов (ОГ) «С», «Г», «М1», «М», «Б», «А».

4. В последующем, через значения глубин опорных и целевых сейсмических горизонтов, соответствующие времена регистрации первых вступлений проходящей волны от опорных и целевых горизонтов получали эталонные значения интервальных скоростей по данным СК и ВСП между этими горизонтами в координатах 3Д, рассчитав, например Vинт_CK_(Б-A)=2dH/dT, где dH - мощность интервала (в метрах), между значениями геологических отметок, соответствующих опорным ОГ «Б» и «А» по диаграммам ГИС; dT - временная мощность интервала (в секундах) между опорными ОГ «Б» и «А», которые служили критериями точности посадки.

5. Далее осуществляли перевод исходной информации из масштаба глубин во временной масштаб V¹(Н) в V(Т).

6. Вслед за этим осуществляли посадку на верхний опорный горизонт литолого-стратиграфического комплекса, например ОГ «Г», введением постоянной поправки dTв. Из-за разных подходов к определению интервальной скорости комплексируемых методов СК, ВСП и МОГТ при пересчете во временной масштаб не достигается соответствия значений времени геологической отметки верхнего опорного горизонта значениям времени экстремума соответствующего опорного горизонта временного разреза МОГТ. Эту разницу dT измеряли на разрезе МОГТ и вводили в исходный годограф.

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

8. После чего осуществляли посадку на нижний опорный горизонт литолого-стратиграфического комплекса, например ОГ «Б». Такое же несоответствие времен геологической отметки во временном масштабе и значения времени экстремума нижнего опорного горизонта на временном разрезе МОГТ, о котором шла речь в п.6, а именно: несоответствие значений времени геологической отметки нижнего опорного горизонта полученного годографа значениям времени экстремума соответствующего опорного горизонта временного разреза МОГТ, наблюдалось и в интервале ОГ «Б». Измеряли разницу dTБ=t_Б ^ОГТ-t_Б ^СК на временном разрезе МОГТ и

рассчитывали поправку $$dTi=Hi*\frac{dT6}{H6}$$ , где

dTi - поправка, вводимая в каждую i-тую точку последнего годографа;

dT_Б-=t_Б ^ОГТ-t_Б ^СК - поправка, которую нужно ввести в последний годограф для посадки на горизонт Б, т.е. совмещения в масштабе времен значения геологической отметки ОГ Б и значения времени экстремума соответствующей отраженной волны;

t_Б ^ОГТ - значение времени снятое с временного разреза МОГТ (в точке скважины) для экстремума отраженной волны, соответствующей горизонту Б;

t_Б ^CK - значение времени, снятое с последнего годографа СК для геологической отметки горизонта Б;

Hi - значение геологической отметки i-ой точки последнего годографа СК;

НБ - значение геологической отметки отражающего горизонта «Б» по диаграммам ГИС, уточненное по данным сейсморазведки.

Вводили поправку Тиспр=Тпосл+dTi в последний годограф п.6.

9. Затем наблюдали соответствие значения времени геологической отметки во временном масштабе и значения времени экстремума отраженной волны в интервале ОГ «Б» и вновь осуществляли контроль качества годографа сравнением с эталонными значениями интервальных скоростей.

10. Далее переносили точки годографа на акустически слабые границы - глины.

11. Вновь осуществляли контроль качества годографа сравнением с эталонными значениями интервальных скоростей и сохраняли априорный годограф.

12. Причем для скважин с глубинами забоя в самом нижнем литолого-стратиграфическом комплексе, например ниже ОГ «Б», фиксировали точку годографа, соответствующую ОГ «Б» (НБ - значение геологической отметки ОГ «Б», ТБ - значение времени экстремума отраженной волны ОГ «Б», снятое с временного разреза МОГТ).

13. После этого вводили точки годографа, соответствующие глинистым пачкам.

14. Затем рассчитывали значение времени для точки забоя скважины годографа как Тзаб=Т_Б+dT, где dT=2dH/Vинт (Б-А), где

dH - расстояние в м от значения геологической отметки ОГ Б до значения глубины точки забоя;

ТБ - значение времени, снятое с временного разреза МОГТ и соответствующее экстремуму от ОГ «Б»;

Vинт(Б-А) - интервальная скорость по данным СК и ВСП, полученная в п.4, для самого нижнего литолого-стратиграфического комплекса, например для интервала между ОГ «Б» и ОГ «А».

15. Вслед за этим осуществляли контроль качества годографа сравнением с эталонными значениями интервальных скоростей.

16. После чего сохраняли априорный годограф.

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

Изобретение относится к области геофизики и может быть использовано для получения сейсмических разрезов изображений геологической среды. Способ включает последовательные действия, при которых получают и подготавливают данные методов общей глубинной точки, сейсмического каротажа, вертикального сейсмического профилирования, акустического каротажа, плотностного гамма-гамма каротажа и проверяют качество этих данных, а также получают эталонные значения интервальных скоростей. Получают исходный годограф и рассчитывают синтетическую сейсмограмму. Затем проводят контроль качества и вводят постоянную временную поправку для посадки на верхний опорный горизонт литолого-стратиграфического комплекса. Затем вновь рассчитывают синтетическую сейсмограмму и вновь проводят контроль качества. Вслед за этим рассчитывают и вводят поправку для посадки на нижний опорный горизонт литолого-стратиграфического комплекса. После этого вновь рассчитывают синтетическую сейсмограмму и осуществляют контроль качества. Переносят точки полученного годографа на ближайшие акустически слабые границы. Повторно рассчитывают синтетическую сейсмограмму с последующим контролем качества и получают априорный годограф. Технический результат - повышение достоверности и точности соответствия горизонтов временного разреза и геологических отметок скважины. 10 з.п. ф-лы, 2 ил.

1. Способ получения априорного годографа для выполнения литолого-стратиграфической привязки путем получения и подготовки данных акустического каротажа, сейсмокаротажа, плотностного гамма-гамма каротажа, калибровки данных акустического каротажа, перевода во временной масштаб, расчета последовательности коэффициентов отражения, свертки их с исходным сигналом, т.е. построения синтетической сейсмограммы и кривой акустической жесткости с повторением циклов, отличающийся тем, что он включает в себя ряд последовательных действий, при которых независимо получают и подготавливают данные методов общей глубинной точки, сейсмокаротажа, вертикального сейсмического профилирования, акустического каротажа, плотностного гамма-гамма каротажа, геологических отметок и проверяют качество этих данных, а также получают эталонные значения интервальных скоростей, после чего получают исходный годограф и рассчитывают синтетическую сейсмограмму, затем проводят контроль качества, далее вводят временную поправку для посадки на верхний опорный горизонт литолого-стратиграфического комплекса, затем вновь рассчитывают синтетическую сейсмограмму, вновь проводят контроль качества, вслед за этим рассчитывают и вводят поправку для посадки на нижний опорный горизонт литолого-стратиграфического комплекса, по формуле
dTi=Hi∗dT_н/H_н, где
dTi - поправка, вводимая в каждую i-тую точку последнего годографа;
dT_н-=t_н ^ОГТ-tн^СК - поправка, которую нужно ввести в последний годограф для посадки на нижний опорный горизонт, т.е. совмещения во временном масштабе значения геологической отметки нижнего опорного горизонта и значения времени экстремума соответствующей отраженной волны;
tн^ОГТ - значение времени, снятое с временного разреза метода общей глубинной точки (в точке скважины) для экстремума отраженной волны, соответствующей нижнему опорному горизонту;
tн^СК - значение времени, снятое с последнего годографа сейсмокаротажа, вертикального сейсмического профилирования для геологической отметки нижнего опорного горизонта;
Hi - значение геологической отметки i-ой точки последнего годографа сейсмокаротажа, вертикального сейсмического профилирования;
Н_н - значение геологической отметки нижнего опорного горизонта;
после этого вновь рассчитывают синтетическую сейсмограмму с последующим контролем качества, и получают априорный годограф, сохраняют его, синтетическую сейсмограмму и кривую акустической жесткости.

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

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

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

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

6. Способ по п.1, отличающийся тем, что в качестве исходного годографа берут годограф, значения которого сняты со значений глубин опорных и целевых горизонтов, соответствующих времен регистрации первых вступлений проходящей волны, используемых для получения эталонных значений интервальной скорости в координатах 2Д, 3Д.

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

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

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

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