Файл: Применение ЦВМ и средств вычислительной техники в геологии и геофизике [сборник]..pdf

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

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

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

Добавлен: 30.10.2024

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

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

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

I(au a2,

a„

d] e - ^ ' m j — ф.(со)] X

 

X 2

a je - imTmj— (p(cö)]dco.

(24)

j = 1

 

 

 

Здесь область D представляет собой совокупность отрезков,, по которым берется интеграл: ф(Ш)— вещественная функция, к ней строится приближение: хт} = т от (£/).

Необходимые условия минимума дают

^ = 2|{ Д

ö/.os[ü)(ту—Tft)] — ф (со) cos (сот*} rfco = 0.

(25)

 

,

(k = 1 ,

2, ... , п)

 

В качестве области D выступают отрезки, заключающие весь

диапазон доминирующих частот

со, т. е. положим

£ [b, с],

Функцию ер (со)

возьмем равной 1

на отрезке [Ьі, Cj]

на котором доминирует частота волны с

годографом хт (|).

При этих условиях имеем систему уравнений для определения величин üj

in |(

 

cb

 

 

 

 

п I (zmj Tmft)

COS

 

 

c + b'

2 а'

 

 

 

 

zmj

zmk

(Tm. — zm.)

 

 

k'

2 J

J=1

c-i — b.

 

 

 

 

 

= — sin ( X

COS X

ci +

j (k

1, 2, ... , n) (26)

■mk

 

'mk

 

 

z mk

 

 

 

 

 

 

В этом случае прямо никак не учитывается выполнение вто­ рых равенств (23). Если хт/хтк = (у —k) Ат, что может быть при неравнодистантном расположении сейсмоприемников или при прямолинейности годографа, то матрицу системы (26) можно привести к диагональному виду изменением разности сЬ. Действительно, при натуральном k и соблюдении ус­ ловия

( ^ - *»*) Cj=~ = (/ - Ь) Ч г 6 At = kn

(27)

получаем равными нулю все элементы матрицы системы (26), не расположенные на главной диагонали.

Для оценки эффективности приема отдельной волны интер­ ференционной системой нередко используется так называемый коэффициент направленного действия (кнд) [2] — отношение

12


«энергии» данной волны на выходе интерференционной систе­

мы к «энергии» той же «волны при условии ее синфазного прие­ ма, т. е.

00

|5а,I(и) Hk, I(ü>) I2 da>

J

a?k, I

OO

I Sk, l (tu)

J

где

 

 

 

a*.i=

J

h(l)Ak a {i)dl.

— П

 

 

Если произведение h (g)

Ak, i

(g) знакопостоянно навеем

интервале [—а, а], получаем кнд<Д. Действительно, по опреде­ лению

\Hk,i((ü)]2= j'

j'

h(l)h(y\)Ak>i(l)Akt /,(r))exp{t'[g-(ri) —

—a—a

-ë{i)]dT\dl = 2 j'

J [/г(|)/г(г))Л й, гіа ) Л А, г(г1)соз[я(т))-

 

—a—a

-g (5 )]< V g < 2

J

)\ h { l) h { 4)Ak>l {l)Ak tl (4)}d4dl =

 

—a —a

= 2 j J

 

I (QAk, I (v})dr\dl = a2krl.

— a—a

 

 

Знак модуля в левой части предпоследнего равенства можно опустить только при знакопостоянстве произведения ft (g)

A k . i (I).

Если это произведение меняет знак на интервале [—а, а], то кнд может неограниченно возрастать, что происходит, на­ пример, в случае, когда произведение чувствительности прием­ ных элементов на интенсивность волны постоянно по модулю, меняет знак от элемента к элементу и число элементов четное. Кнд в этом случае не может служить мерой эффективности

приема волны.

Исследуем далее с позиций лучевого метода способ С. А. Нахамкина последовательных вычитаний волн, используя только основную идею способа и не накладывая ограничений, принятых в [4] и других работах упомянутого автора, на при­ нимаемые волны.

13


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

Для простоты будем считать, что в (6) погрешность еД, £) равна нулю и годографы (£) вступлений волн таковы, что их можно упорядочить по временам вступлений следующим образом:

дл я к < І'

(28)

т. е. анализируются годографы 1, 2, 3, показанные на рисунке сплошными линиями, и не рассматривается вся совокупность годографов первых вступлений (1—6). Штриховыми линиями показаны времена, соответствующие задним фронтам воли. После этих времен волна исчезает, т. е. становится тождествен­ ным нулем. Равенство нулю погрешности еД, £) означает, что рассмотрение проводится в рамках лучевого метода решения уравнений Ляме и погрешности приемной аппаратуры не учи­ тываются.

Вследствие финитности (3) функций /*, і (t) и условия (28) имеем на интервале тД£)<ДтД|) представление функ­ ции (6):

F(t

г

ДЛЯ т,‘(6 )< * < т 2(!)

 

ІО

для *<тД£).

(29)

На сейсмограммах свойство (29) функции F(t, £) выражает­

ся очень четко — после покоя следует вступление

? = тД |)

первой волны и поэтому его значение снять легко.

по ли­

Путем синфазного сложения я-трасс сейсмограммы

нии ^= ті(|) с учетом чувствительности /іД |) получаем суммо­ ленту, на которой доминирует волна

£ и , Л Т )

2 A l , l {tm)hl ( t m = 2 A lf u l

(Т) =

I

тп= 1

I

 

 

= 2 Л Д м [(-т(Е)].

(30)

 

 

I

 

Остальные волны за счет интерференции будут ослаблены. Определяем приближенно задний фронт выделенной волны

14


(штрих-линия 1 на рис.) по резкому спаду амплитуды колеба­ ний. Пусть это будет при /=Ті ( |) +Дгі. Полагаем

Fx{t, Е) =

І

ҢГ + Т1( | т ), I M

U

= 0

m= 1

 

 

 

 

для t > T i ( £ )

+ Дть

 

(31)

T. e. для t > x i (E) +Ati

в о л н ы

гасятся

лабораторным" путем.

Приближенно можно считать, что

 

 

F,{t,

Е)

« 2 Л ^ Ь1[/-т,.(Е )].

(32)

б V

ч

Годографы первых вступлений и задних фронтов волн, а — первые вступления; б— задние фронты волн.

15

Используя (32), (6) и учитывая,

что е((,

£)=0, составля­

ем разность

Аи

(S)fw{f— с*(5,)] -

F(t,lT) - a rF,(t,tr) « S

к, I

 

 

 

- a r -EAtf u #

- 1

1 (h)] =

(33)

I

 

 

 

= ^ f u t[t — Xl(lr)][Aui(lr) - M i l +

+s 4 i ( U / * . f P — t r(y ]-

ft, г fc-фі

Здесь в соответствии с (30)

 

 

 

 

 

 

 

 

A t = i A i.« (S JM U -

 

 

 

(34)

 

/7г=1

 

 

 

 

 

 

 

 

Вычитание волн с вступлением / =

Т і ( | )

будет полным толь­

ко в двух случаях:

 

 

 

 

 

 

 

 

 

1) если Ль/ ( У — аЛ і = 0

(/,

г =

 

1, 2 ...),

 

 

т. е. если интенсивности А и

(1=1,

2,

...,)

отличаются

одна

от другой лишь постоянным

множителем.

Действительно,

в

этом случае

 

 

 

 

 

 

 

 

 

Аі, I (£,-)

сИ і(5,)

 

=

а г

 

(35)

 

Лг

П

 

 

 

 

(£ш)

 

 

 

 

 

 

сі

^

 

 

 

 

 

 

ш = 1

 

 

 

 

 

 

 

при любом /;

( / = 1 , 2 , . . . ) ,

 

 

 

 

 

2) если f u l (T) =h( T)

 

 

 

 

 

т. е. если волны с вступлением / = ті(£)

имеют

одну и ту

же

форму.

 

 

 

 

 

 

 

 

В этом случае

I Ak, I (£,)

 

 

 

 

 

 

 

а.

 

 

 

 

(36)

 

 

 

 

 

 

 

Ясно, что в практических случаях для

определения

пара­

метров а г пользоваться формулами

(35),

(36)

нецелесообраз­

но.

Эти параметры подбирают из

условия, чтобы разность

F(t,

| г) — arFi(t, | г ) на интервале Ti(g)^£$6rii(£) + A tt прини­

мала в среднем наименьшее значение. Может, конечно, слу­

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

не

позволит

выделить следующую волну. Если интенсивности

А и

(I)

мало изменяются в зависимости от своего аргумен­

та,

то расстояния между нулями функции 2 A u ( l ) f u

(Т),

16


обусловленные интерференцией волн (Т), тоже мало изменяются в зависимости от | и поэтому можно сумми­ ровать трассы по кривой расположения по % нулей первой ярко выраженной фазы колебаний, а не по вступлению t = Ті(£), так как оно обычно подвержено большим аппаратур­ ным искажением. На трассах сейсмограммы можно заметить большую стабильность в расстоянии между первыми нулями функции F(t, I) от трассы к трассе.

Если процесс вычитания волн

с годографом ti( |) удался,

то на разностной сейсмограмме

четко выступит вступление

или первая доминирующая фаза волн с этим вступлением. Процесс сложения-вычитания можно повторить уже по от­

ношению к этим волнам и т. д.

В случае более сложной волновой картины, когда на одних пикетах до точки М (рис.) первыми приходят волны 1 с вступ­ лениями / = ті(|), а далее — волны 6 с вступлениями t Тб(|), разделение волн предлагаемой модификацией возможно толь­ ко, если форма первой из них очень сильно отличается от формы другой. Тогда легко найти положение точки М и про­ цесс сложения-вычитания следует осуществлять по участкам.

В ы в о д ы

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

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

 

 

 

Л И Т Е Р А Т У Р А

1.

Г о л ь ц м а н Ф. М.

Основы теории интерференционного приема ре­

гулярных волн. М.,

«Наука», 1964.

2.

Б е с п я т о е Б.

И.

Некоторые вопросы теории группирования в

сейсморазведке. Прикладная геофизика, вып. 25. Гостоптехиздат, 1960.

3.

Б е с п я т о в Б .

И.,

Ю р ч е н к о В. Г. Применение теории направлен­

ного интерференционного приема для оценки и обоснования методики груп­ пирования взрывов по способу непрерывного линейного источника («плоский

12. Заказ 1928 ,__________ 17

'Кчк

■:н;: _

О ч