Файл: Цифровая обработка сейсмических данных..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 09.04.2024

Просмотров: 135

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

e(t)

Sg(tJ

 

U0(t)

Р и с . 26. С о о т н о ш е н и е д в у х в а р и а н т о в м о д е л и о д н о к а н а л ь н о й с е й с м о г р а м м ы .

как переходный процесс среды, связанный (при непрерывном вре­ мени) с s0 (t) соотношением [38]

Для дискретного времени, по аналогии с (2.23), будем считать

и0

(t) =

«о (0 —«о (< —АО

(2.28')

 

 

At

 

Соотношения (2.24) и (2.28') связывают два варианта модели: вариант, в котором рассматриваются элементарные волны Гюйгенса, и вариант, где рассматриваются огибающие этих волн. Одна и та же сейсмограмма может быть выражена через обе модели (рис. 26):

г - £ ( < ) — E(t — At) j =

|~ s0(t)-s0(t-At)

2 Aks0(t-6t)=SS

aku0(t-qik).

(2.29)

j _ j к

к j

 

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

Поэтому в качестве модели, непосредственно отображающей принцип Гюйгенса, будем использовать выражение, полученное из (2.15), (2.22) и (2.28) и аналогичное по своей структуре форму­ лам (2.15), (2.15') и (2.15"):

У W --= a {[z0 (t) + Щ (01 * sn 0 B (t — Атс ) +

+ « п о в (t) +

nM ( * ) } * S p e r ( 0 .

(2.30)

где

 

 

 

 

z0(t)

=

[c(t)K(t)]*u0(t);

(2.31)

n0{t)

= [cn{t)Xn(t)]

* un(t),

(2.32)

68


или

 

 

у (t) =r-a{[c(t) Xy (t)] * u0 (t) * sn o B (t — AT C ) +

 

 

 

(2.30')

V(t) = a{[c (t) ky {t)]*u(t-

Л т с ) + n (t)}.

(2.30")

М О Д Е Л Ь М Н О Г О Т Р А С С О В Ы Х

СЕЙСМОГРАММ

 

Построим эту модель для случая профильных наблюдений с взры­ вами в точках х = 0, Ах, 2Ах, ЗДж . . . и регистрацией волн на расстояниях \ = 0, ±Д£, - t2A| . . . от каждого пункта взрыва. Вместо одиночной трассы у (t) будем теперь рассматривать совокуп­ ность сейсмограмм ух (£, t). Хотя обе новые переменные х и | харак­ теризуют расстояния по профилю, они могут меняться независимо друг от друга. Поэтому ух (£, t) представляет собой трехмерную функцию, которую можно изобразить, как показано на рис. 27, а. Цифровые сейсмограммы, у которых отсчеты выбираются в моменты времени t = 0, At, 2At . . ., представляют собой трехмерные вре­ менные ряды.

Р и с .

27. О с о б е н ­

ности

м о д е л и

м н о ­

готрассовых

с е й ­

смограмм .

 

а — м н о г о т р а с с о в ы е сейсмограммы как т р е х м е р н ы е ф у н к ц и и ; о — б л о к - с х е м а о д ­ н о г о и з к а н а л о в м о ­

д е л и

м н о г о т р а с с о ­

вых

с е й с м о г р а м м .

Задержка

>>*(t,t)

Задержка

69


При переходе от одиночных трасс к сейсмограммам в рассмотрение вводится важнейшая характеристика сейсмического поля — пове­ дение волновых фронтов. Будем различать волны регулярные и не­ регулярные. У каждой из регулярных волн «приведенная» ампли-

туда —

= с (t) к

(t)

и

форма

импульса

s (t)

не зависит

от х и | ,

а

время

пробега

Qk

(?)

в толще

ниже

уровня приведения связано

с

характеристиками

среды hk

(t)

и v (х, К)

уравнением

годографа

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

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

дой конкретной трассе вместо истинных величин A,

s0

(t)

и 6 имеют

место

случайные

величины

А + ЬАХ (?),

s0

(t)

+

6s0 x (?, t)

и Qk (t)

+ &®kx (l),

где 6 4 ,

(?),

&s0x ( i , t),

8Qkx (?) — случайные

составляющие. To

же самое

можно сказать и

о всех

других пара­

метрах, которые включены в рассмотрение: Атс , s (t) и т. д.

 

Понятие нерегулярная волна также является

условным: полная

нерегулярность наблюдается только тогда, когда

расстояния

между

соседними пунктами приема превышают

радиус корреляции

волн

по

профилю — величину, которая у всех

реальных волн отлична

от

нуля.

 

 

 

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

Запишем теперь общее выражение для модели многотрассовых сейсмограмм. Очевидно, для этого могут быть использованы выраже­ ния (2.15), (2.15') и (2.15") или (2.30), (2.30') и (2.30") без каких-либо изменений в структуре, но с заменой функций одной переменной у (t), 6 (t — 0К ), s (t — Дтс ) . . . на соответствующие функции трех пере­

менных ух (g, t), б [t — QkK (?)], s [t — Axcx (?)] . . . Хотя, как мы только что отметили, практически все параметры волн содержат

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

тов вступления волн Qkx

(?) и

qkjx

(?). Поэтому

включим в

модель

только случайную компоненту 8Qkx (?) или, соответственно,

bqkjx{^,).

Уравнение годографа

Qkx (?) при ср =

0 всегда можно

записать

в виде Qk (х) + Атк (?),

где

Qk

{x) =

Qkx(l)

при ? =

0.

Вели­

чину Атк (?) будем называть кинематическим сдвигом; при обработке

этот

сдвиг

компенсируется

кинематической

поправкой.

Вели­

чина

Атк (?)

имеет важную

особенность:

значение Qk

(х) =

70


= Qkx (I) Атк

(£), вычисленное по наблюденному годографу 6f t j e (i)

в случае cp Ф 0,

с удовлетворительной степенью точности совпадает

со значением нулевого времени, которое наблюдалось бы в точке |/2 при помещении в эту точку и источника, и приемника. С учетом сказанного получаем

 

Х(|, < - А т с Л £ ) ) + Яповх(1. t)-±-nwx{l,

t)}*sper(t),

 

(2.33)

где согласно (2.5), (2.11) и (2.26),

(2.31) компонента z0x (g, f)

дается

одним из выражений

 

 

 

 

 

 

 

 

V ( 1 , 0 = р (0 2 xfc6 (П,)~! * s0(t),

 

(2.34)

или

 

 

 

 

 

 

 

 

 

Чх ( 5 , 0 = [с (<) 2

2

е (<?*,;)] * и 0 (0.

 

(2.34')

Здесь

 

 

 

 

 

 

 

 

 

=

* - [в* (ж) +

Атк

(g) + 66*Л (|)];

 

(2.35)

 

=

* - [ ? * (*) +

 

(S) +

fy**/ (£)].

 

(2.36)

По аналогии с (2.15") и (2.30") можно записать

 

 

М Е .

0 = М £ Н И * К , и ,

01 * 5 [ « - А т с Л Е ) ] + «ЛЕ,

«)}»

(2.35')

или

 

 

 

 

 

 

 

 

г/,(6,

0 = М £ Н И * ) М £ >

0 ] * » U - A ? C * ( S ) ] + » * ( S .

*)}•

(2-36')

Блок-схема одного из каналов модели многотрассовых сейсмо­ грамм изображена на рис. 27, б. Ее основной отличительной особен­ ностью по отношению к блок-схеме одноканальной модели является добавление кинематических сдвигов в каждый канал.

С Т А Т И С Т И Ч Е С К И Е ОСОБЕННОСТИ М О Д Е Л И СЕЙСМОГРАММ

Из величин и параметров, входящих в выражения (2.33)—(2.36), нас интересуют статистические свойства амплитуд и моментов всту­ пления волн на отдельной трассе, особенности формы импульсов, формы годографов и статических поправок. Построим их модели.

Модели распределения амплитуд и моментов вступления волн.

Рассмотрим введенную последовательность Ах (t) = с (t) %х (t). Мо­ делью последовательности кх (t) при фиксированных х обычно считают стационарный временной ряд с некоррелированными, нор­ мально распределенными значениями к.

Однако некоторые экспериментальные исследования показывают, что при фиксированном х величины и распределены не нормально, а по закону Джеффриса [53, 54]. Среди множества значений х, близких к нулю, встречаются отдельные довольно большие значения,

71


 

 

 

 

 

 

 

 

 

 

причем

они попадаются чаще,

чем

 

I

i

l

l

I

 

I

I

I I

э т о г о следовало

бы

ожидать при

'|1

11 Till

II ч

1^11

I V

i V ' l I

нормальном

законе.

Качественное

I N

I

I

I

I

 

'I

I

отличие

между

моделью

последо­

$

 

 

 

 

 

 

 

 

 

вательности

хх

(t)

при

распреде-

 

 

 

 

 

 

 

 

 

лениях

Гаусса

и

Джеффриса

ил-

 

I

j

 

 

 

I

I

 

 

люстрируется рис.

28.

 

 

 

' у К " Л 1

п " I '• I •1 it I' * * | 1 " "| 111111

Джеффрисовое

распределени

 

 

 

 

j

 

 

 

I

'

интуитивно

представляется

«само

п

0 0

 

 

„ „ „

 

 

собой разумеющимся» для толсто-

Р и с .

28.

П р и м е р

последовательно-

 

v

J

 

-Г,

 

 

 

*

 

с т е й в е л и ч и н

п р и

Г а у с с о в о м

(а)

и

слоистых толщ. Если пренебречь

Джеффрисовом

(б)

распределениях,

«фоном»

значений

х,

близких

 

 

 

 

 

 

 

 

 

 

к

нулю, последовательность

кх

(t)

для такой среды может рассматриваться, как стационарный случай­ ный поток [5] с независимыми значениями, т. е. случайная последо­ вательность значений %, которые появляются не в каждый данный

момент времени t — О, At, 2At . . ., а реже, и

характеризуются

определенной плотностью потока — количеством

ненулевых значе­

ний % в единицу времени или же средним значением интервала времени между двумя соседними ненулевыми х. Простейшим при­ мером потока является пуассоновский процесс, или пуассоновский поток, который представляется интуитивно весьма «подходящей» моделью для кх (t) в случае толстослоистой среды.

Чаще в качестве модели для %х (t) используется временной ряд, который можно рассматривать как частный случай потока с постоян­ ными At. Независимо от того, считается последовательность %х (t) случайным временным рядом или случайным потоком, она имеет

нулевое

математическое

ожидание

и корреляционную

функцию

Ьх (т ) =

х 2 б

(т),

равную нулю везде, кроме точки т =

0.

Двумерное

поле значений

% (х, t)

представляет

собой идеальный

временной

разрез

— совокупность

линий

нулевого времени

Qk (х) (к

= 1,

2, . . ., К),

каждой гз

которых

соответствует коэффициент

отра­

жения %k. В предположении малых углов падения будем считать линии Qk (х) стационарными случайными функциями.

Все сказанное о кх (t) относится и к %х (t) для всех типов волн, кроме многократных и поверхностных. Количество кратных волн, приходящихся на единицу времени, возрастает с t, поэтому для них К" (t) не является стационарным рядом (или стационарным потоком).

Моделью последовательности ах (£) также служит двумерный случайный процесс, некоррелированный по х и по | .

Моделью функций с (t) для разных типов волн, как уже гово­ рилось, может служить убывающая экспонента (см. рис. 22). У раз­ ных типов волн характер спадания интенсивности с временем t резко различен. Одни помехи затухают быстрее полезных волн (поверх­ ностные, иногда — поперечные и обменные волны), другие — мед­ леннее (микросейсмы и многократные волны). Последние и являются основным фактором, ограничивающим глубину исследования.

72