Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Старобинец, А. Е. Состояние региональных и поисковых работ методом преломленных волн на нефть и газ за рубежом

.pdf
Скачиваний:
3
Добавлен:
19.10.2023
Размер:
3.78 Mб
Скачать

Пример конкретного использования суммирования для решения рассматриваемой задачи схематически иллюстрируется на ри с.7 и сводится к тому, что записи предварительно спрямляются в соответ­ ствии с различными значениями градиентов, а затем суммируются.

___ ГТТТЧ\. -Я ,ffi_ffl_flLJh

V V W w *

----- ““КГ

^

—/Т\ А Л /П Л \.

 

 

 

IT

 

 

-im am_________ - < rrm _

.

U U

 

. ...... J

L

Л

Л

 

 

лттм\ _

a

 

 

n

 

Рис.7, Схемы, иллюстрирующие принцип способа

суммирования

а - имитация записи двух преломленных волн с разными кажу­ щимися скоростями; б - суммотрассы при различимое градиен­

та х , принятых при спрямлении

Суммарные трассы, отвечающие градиентам S j и 9г ' П03В0ЛЯЮ* находить соответствующие водны. Основным недостатком этого спосо­ ба является необходимость выполнения требования слабого изменения (секущихся скороотей целевых волн. В противном случае будет наблю­ даться противоречие между необходимостью увеличения числа сумми­ руемых трасс для улучшения статистического эффекта и противопо­ ложным требованием уменьшения числа привлекаемых трасс из-за из­ менения градиента.

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

Второй из рассмотренных путей выполнения автоматической корреляции реализован, в частности, в объединенной системе стан­ дартных обрабатывающих и интерпретационных вычислительных прог­ рамм [65], разработанной во Франции совместными усилиями геофизи­ ческой группы SN FA (Soclete N a tio n als dee P e tro le a d'A qutane)

Исследовательского центра и Технического и Научного ведомств

CGG(Compaquie G enerals de G eopbyaique).

- 29 -

Автоматическая корреляция вступлений преломленных волн осуществляется двумя программами РЕФА (refa ) и РАПТ (rapt) , входящими в ату систему. В программа РЕФА путем выравнивания двух выбранных тр асс М и N и последующего суммирования получается контрольная трасса Р (р и с.8 ) . Корреляция других трасс Т выпол­ няется путем определения величины времени, которое максимизирует нормализованную функцию взаимной корреляции с контрольной трас­ сой Р в соответствии с выражением (2 ) . По вполне очевидным причи­ нам интегралы в зтом выражении заменяются рядами:

 

 

£

р: Г:

ф .

=

i -1

« V

 

 

(2 *)

Рио.8 . Схема, пояоняющая принцип автомати­ ческой корреляции, положенный в оонову програм­

мы "РЕФА"

Дискретная величина j , которая максимизирует P j , каждый раз запоминается, чтобы установить дискретное время корреляции. ЗапоминаЬтся также соответствующие этим времевам расстояния до источника возбуждения. Кроме того, параллельно определяются скорость и значение "Jf " и в случае перекрытий годографов вычис­ ляются начальные константы (см.раздел ’'Методы интерпретации дан­ ных МПВ"). Рассчитываются также оредняя величина и стандартное отклонение этих величин для всех пар связанных трасо. Последнее рассматривается как критерий относительного качества корреляции. В конце этой фазц обработки откоррелированная запись Выдается на магнитной ленте. В соответствии с программой РАПТ для контроля

- 30 -

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

Несколько иной процеос осуществления автоматической корре­ ляции был применен в ФРГ [3 5 J. При вычислении функции взаимной корреляции в этом случае производятся только такие взаимные пере­ мещения сейсмических трасс, которые соответствуют скорости реду­ цирования сейсмозаписей на монтаже. Таким образом, аргументом является скорость, а не величина сдвига. При практическом приме­ нении этого способа отрезок записи в определенном временном ин­ тервале коррелируется с отрезком другой записи в ином интервале. Временной интервал зависит от удаления данной трассы от источни­ ка возбуждения и скорости редуцирования. Цифровая корреляция с использованием ЭВМ осуществляется путем постепенного поиска по­ правок к скорости редуцирования. При таком подходе функция взаим­

ной корреляции выражается в следующем виде:

 

 

 

т

 

 

 

 

 

< P ( V ) ~ £

M j(V ) ,

при этом

 

 

 

J=1

J

 

 

 

 

 

 

(3)

где

п

-

число трасс;

 

 

 

/7?

-

число цифровых значений рассматриваемой трассы;

 

V-

-

принятая скорость редуцирования;

R l( ls l,...,n)- удаление трассы от

источника "возбуждения;

 

)

-цифровые значения всех трасс.

Недостатком

этого способа

является то,

что сейсмозаписи до

некоторой степени

корродируются вдоль прямой,

т .е . с постоянной

кажущейся скоростью. Этот недостаток пытаются компенсировать тем,

что все

трассы корродируются последовательно и со всеми другими,

а затем

находятся суммы корреляционных функций.

- 31 -

Введение кинематических поправок

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

с введенными таким образом кинематическими поправками группируют­ ся в непрерывные предварительные "временные разрезы” путем вычи­ тания связывающих конотент. Получение окончательных "временных разрезов" (времен 11 tg " или временных задержек на пунктах прие­ ма) осуществляется в последующем по методике, которая будет налог

жена ниже. На ри с.9, представляющем

собой информацию, получае­

мую на выходе программы КРЕЛ (CREL )

[65], в нижней части при-

.9 . Форма визуализации информации на выходе программы "КРЕЛ"

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

- 32 -

Специальные опарации выделения целевых преломленных волн на сейомозапиоях

С целью выделения целевой волны на фоне случайных или ре­ гулярных помех или других волн, например, в воне их интерферен­ ции, а также для повышения раврешенности ааписи во многих стра­ нах при обработке данных Ш1В в последние годы используются спе­ циальные приемы,основными из которых являются суммирование с вре­ менными задержками, веерная фильтрация и деконволюция.(Еще сравни­ тельно недавно для этих целей в МПВ применялась одна только по­ лосовая фильтрация .)

Суммирование с временными задержками

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

Способ для улучшения регистрации выбранных преломленных волн и подавления нежелательных волн и случайных помех, базирую­

щийся на разновременном

суммировании, был недавно запатентован

в СМ

[71]. Разработана

и соответствующая обрабатывающая аппарату­

ра. Эффект достигается

путем суммирования записей от разных

пунк­

тов

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

пос­

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

- 33 -

• Другой путь расчета временных сдвигов, вводимых при сумми­ ровании, который использует в своей основе метод Гарднера (ом. ниже), применительно к многократному перекрытию в МНВ рассматри­ вается в работе [30] .В соответствии с этим методом рассчитывают­ ся разности временных задержек на пунктах возбуждения и пунктах приема и истинное значение граничной скорости. После этого наблю­ денные времена пробега трансформируются во временные задержки на пунктах приема и перекрывающиеся записи суммируются. В результа­ те можно представить все вступления от различных преломляющих границ на одном " временном р азрезе ".

Веерная фильтрация

Каких-либо существенных отличий в применении веерной фильт­

рации к обработке

данных МНВ по сравнению с МОВ не

имеется.

Основ­

ное ее назначение

- разделение преломленных волн с

разными

кажу­

щимися скоростями

и со сравнимой интенсивностью в

зонах интерфе­

ренции. Выше уже отмечалось, что для улучшения селекции прелом­ ленных волн при использовании окоростной фильтрации в цифровой форме иногда применяются линейные группы источников возбуждения [41] . В случае, если приповерхностные отложения характеризуются латеральными литологическими изменениями и сложным отроением, скоростная фильтрация может быть успешно выполнена после пред­ варительной коррекции всех записей, предусматривающей их отне­ сение к достаточно глубокой линии приведения [3 9 ].

Преломленные волны, как правило, характеризуются значи­

тельным сдвигом

по времени

не записях, поэтому необходимо рас­

считывать очень

длинные во

времени многоканальные

скоростные

фильтры. Значительно более

короткие фильтры могут

быть

расочита­

ны с помощью трансформации

исходных записей с заданной

скоростью

редуцирования [4 1 ].

 

 

 

- 34 -

Деконволюция

Деконволюция (обратная.фильтрация) представляет особый ин­ терес при обработке данных сейсиических работ МНВ на норе, по­ скольку ее применение позволяет очистить сейсыоаапись не только от ревербераций, но также и от волн, порождаемых пульсацией га­ зовых пузырей, которые образуются при взрывах в воде. Последнее особенно важно, так как дает возможность производить взрывы на больших глубинах с целью понижения расхода ВВ. До недавнего вре­ мени этот метод цифровой фильтрации широко применялся только при работах иЮВ в пределах акваторий.

Применение деконволюции при обработке данных морской съем­

ки Ш1В,

проведенной ФРГ на Неверном

море в 1966

г . [3 8 J, показа­

ло, что

она может быть очень ценным

инструментом

при такого рода

работах.

При проведении назеиных региональных исследований Ш1В с

большими удалениями от источников

возбуждения, помещенных в во­

д е , атот

метод обработки был также

взят на вооружение[5 5 ,8 7 ].

Предполагается, что эта операция будет аффективной в отношении подавления любого прогнозируемого шума, который затрудняет выде­ ление целевых волн и при обработке данных наземных сейсмических работ МПВ _[39].

В то же время отмечается, что осуществить деконволюцию по отношению к записям 11ПВ сложнее, чем к записям ЫОВ [8 7 ]. Главная трудность при атом заключается в выделении базовых волновых им­ пульсов' иди,другими словами, в определении спектра полезного сигнала. Распространенный путь решения этой проблемы, используе­ мый в НОВ, основывается на допущениях, что сейсмический сигнал

слабо

изменяется и является

минимально-фазовым.

Известно, что в

таком

случае фазовый спектр

может

быть найден по

амплитудному.

С определением же последнего дело

обстоит проще.

При указанных

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

Для МОВ это приемлемо, а постепенное изменение базового волново­ го импульса с увеличением времени прихода отраженных волн может быть учтено. Характер записей ШШ является более сложным в том отношении, что пути пробега регистрируемых волн различаются

- 35 -

в зависимости от сложности разреза, и различные типы вступлений не группируются в простые последовательности. Поэтому такое пред­ ставление сейсмозаписи неправомочно. В работе [87], в частности, показывается, что приемлемое приближение в этом случае достигает­ ся при свертке базового волнового импульса с комбинацией единич­ ных импульсов и функционала ( —iV t )

В МОВ оценка базового импульса может быть сделана путем

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

на записи или по функции автокорреляции сейсмической записи.

В МПВ же часто никакого базового импульса не может быть выделено

на запиои и, таким образом, величина функции автокорреляции явля­

ется трудно оцениваемой. Более того,

часто оказывается, что целе­

вой сейсмический сигнал не является

минимально-фазовым. В силу

этих обстоятельств часто для того, чтобы деконволировать запись МНВ, должны быть оценены параметры физической модели, формирующей волновое поле [8 7 ]. С этой целью, например, для очищения записи от волн-помех, порождаемых пульсацией газового пузыря, вблизи места взрыва производится запись пульсации генерируемой энергии для того, чтобы зарегистрированные при атом сейсмические сигналы использовать при цифровой обработке [2 ] .

МЕТОДЫ ИНТЕРПРЕТАЦИИ ДАННЫХ МПВ

Поскольку источники получения данных о скоростях в покры­ вающих отложениях или способы их определения по данным МПВ не отличаются от принятых в СССР, основное внимание в разделе будет уделено методам построения преломляющих границ, причем тем ив них, которые наиболее широко используются за рубежом,.а у нас не применяются.

Определение скоростей в покрывающей толще

Наряду с общеизвестными стандартными методами для получе­ ния информации о скоростях, необходимой при интерпретации данных МПВ, характер распределения скоростей в покрывающей толще часто

- 36 -

определяется путей использования времен первых вступлений пре­ ломленных волн. С этой целью разработаны специальные алгоритмы и программы. Такого рода программа составлена,например, в Канаде [32] . В этой программе оценка функции скорости от глубины по первым вступлениям преломленных волн опирается не следующие допу­ щения: скорость увеличивается с глубиной (нет инверсий) и может быть описана временами первых вступлений (нет выпадающих слоев); все преломляющие границы горизонтальны и все скорости постоянны в горизонтальном направлении; анизотропия скоростей не имеет мес­ та .

Время в каждой точке на годографе рассматривается как вре­ мя вступления от отдельного слоя. Для определения мощности 1С-го

слоя

Z k ) пластовая

скорость которого

принимается

равной

кажущей­

ся скорости

V * ,

используется

общеизвестное

уравнение

 

где

V * -

кажущаяоя

скорость

между точками

к ,

к

+ 1

,

 

 

 

 

 

 

 

 

 

 

(5)

Гк

- время

прихода

волны в точку К,

-

удаление

от

ПВ.

Предусмотрено предварительное сглаживание годографов. Программа написана на языке Фортран- I I для небольшой 6 к -словной вычисли­ тельной машины.

Для подобных оценок широко используется также хорошо из­ вестный способ Вихерта-Герглотца.

Поскольку и з-за неполноты систем наблюдений необходимые для расчетов времена " te " определяются часто путем экстраполяции годографов к пункту возбуждения, для более надежного разделения вступлений, относящихся к различным границам, применяются спе­ циальные способы. Например, предложен статистический способ, основанный на допущении, что все точки годографа от одной грани­ цы должны характеризоваться некоторой вероятностью попадания на отрезок прямой [3 4 ]. Точки годографа автоматически распреде­

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

- 37 -

Вопрос о точности и достоверности подобного рода определе­ ний достаточно хорошо исследован в ряде работ [13,19 и др.].Большинство исследователей отдает себе отчет в том, что использова­ ние таких данных в качестве единственного источника оведений о скоростях может привести к серьезным ошибкам.

Аналитические и графические методы построения прелом­ ляющих границ

Во многих зарубежных странах в силу специфики решаемых там при проведении работ МПВ геологических авдач и широкого исполызования. неполных систем наблюдений наряду с методами построения преломляющих границ, применяемыми повсеместно,большое развитие получили такие методы интерпретации, которые в Советском Союзе не применяются. К числу наиболее часто употребляемых методов, ко­ торые можно очитать стандартными, относятся методы временных за­ держек (метод Гарднера),временных членов, Хейлоа и волнового фронта. Основы этих методов были заложены еще в 30-х годах, и о тех пор они рвзвивались и усовершенствовались применительно к новым геологическим задачам, которые отавились перед Ш1В. Мень­ шее распространение получили такие метода, как "А В С ", "плюоминуо" и д р ., которые по существу являются модификациями стан­ дартных методов.

Ниже рассматриваются принципы наиболее часто употребляемых методов, условия их применения и свойственные им ограничения.

Метод временных задержек

(The time-delay method)

Идея метода была выдвинута Л.Гарднером (США) в 1939 г . [24J. В настоящее время применяются многочисленные модификации метода [1 1 ,25,44].

Сущность метода временных задержек состоит в построении пре­ ломляющих границ по составляющим времен" t0 " , называемым времен-

- 38 -

Соседние файлы в папке книги из ГПНТБ