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

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

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

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

Добавлен: 30.10.2024

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

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

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

Для определения лучевых параметров р2, <72, г2 луча отра­ женной волны в точке М { х, у, z } границы отражения вос­ пользуемся уравнением эйконала и законом отражения: лучи падающей и отраженной волн лежат в одной плоскости с нор­ малью в точке М к поверхности отражения и расположены по разные стороны от этой нормали, т. е.

Р22 + <7г2 + г22 = ■— ,

А1А2А3

Р\Я\Г\

( И )

РіЯ2г2

 

Р2А1 + У2А2 + /"гАз = ■— (РіАі + q\A2 + г 1А3),

 

где АI— величины, пропорциональные направляющим

коси­

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

 

Ді = — z x , Д2 = — zy, A3 = 1.

 

Разрешаем уравнение (11) относительно р2, qi, г2

 

р2 = Рі — аАі,

(12)

q2 = q\ — aA2,

r2 = r1— аA3.

 

Здесь

 

а = 2АіРі + дг?1+ V i

 

Д1* + А/ + Д32

 

Заметим, что соотношения (12) справедливы и для неод­ нородной среды в окрестности точки отражения. В этом слу­ чае значение скорости берется в этой точке.

Напишем уравнения лучей отраженной волны

 

х° = X + р— (z° — г ),

 

 

г2

 

 

 

У + — (2°

■2),

(13)

 

Н

 

 

 

*о _

 

 

Хг = U +

 

 

и присоединим к ним уравнение дневной поверхности

 

2° =

z0 (х°, у0) .

 

(14)

Подставляя (14), (9) в (13), получаем уравнения годогра­

фа отраженных волн

в функции параметров х, у,

определя­

ющих точки границы отражения:

 

 

х° = X + — [г0 (я0, у0) г],

Г2

43


У° = У + - { z о(х°, у°) — г],

(1 5 )

г2

 

г2(х, у) = ^ / (х+л:*)2+(;у—у*)2+ [z(x, у )г*]2 +

I Z0 (х °, у°) - г (х, у)

^2 ѵ ѵ

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

— (12)

 

 

 

Р2 =

X— X*

— öAi

 

 

ѴЧХ

 

 

 

<72

УУ*

аД2,

(16)

 

ѵчх

 

 

 

г2 =

ѵчг -

а А3.

 

При z0(x, J/)=Const равенства (15) представляют собой явные параметрические уравнения годографа отраженной волны.

§3. Построение границы обмена

Вслучае обменных волн будем считать, что вдоль профи­

ля z0= 0 заданы годографы поперечной волны т=т2(х, у) и продольной

т=тіі(х, у) =ах+Ьу+с. (17)

Скорости этих волн соответственно равны Ѵ2 и Vj. Тогда систему уравнений можно записать в виде [2]

 

X — Х2 + — Z = Х\ + — z,

 

 

 

Г2

 

Гх

 

 

У = У 2 + - Z

= У і + - Z ,

(18)

 

 

г2

 

гх

 

t = х2{х2,

У2) +

тут- — ах\ + Ьух + с +

—— ,

 

 

Ѵ22Г2

 

 

Ѵ^Гх

где 20 = 0, рі = а,

<7і = Ъ, гх = —

— (а2 + 62);

р2= S ’ q2= ^

’ Г2=

-

у - щ - (Рг2+^

Разрешая систему уравнений (18) относительно коорди­ нат X, у, z границы обмена в функции параметров х2, у2—ко-

44


ординат точек выхода поперечной волны на дневную поверх­ ность г0 = 0, получаем параметрические уравнения границы обмена:

х= х2 + dp2

у= у2 + dq2 z = dr2,

в которых

(х2, Уг) ах2Ьу2— с

 

d =

 

1

арг + bq2—гз

- (а? + У1)

ѴУ

Ѵ ѵ

 

§ 4. Определение величины бокового уклонения для наклонной границы отражения

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

Для упрощения вычислений выберем дневную поверхность

(рис. 2), совпадающую с плоскостью ХОУ: zQ= 0,

а источник

колебаний расположим в начале координат:

 

х# — у к. = z% — О

у0, z°} за­

Пусть годограф х% в окрестности точки В { х°,

дан уравнением

 

х02 + у02 + 2 dx° + 2 еу° + / ,

(19)

определяющим гиперболоид вращения вокруг оси, проходя­ щей через точку { х°= —d, у°— —е } плоскости ХОУ.

Решая обратную задачу по формулам (5) — (7), получим

v

+

d

 

2 (dx°+ ey°-f /)

 

-

(y° + e)f

c

У2(dx!>+ eyO+ f)

f y f - d * ~ e \

, / f _ d2 _ й 2

2 (dx°+ ey°+ /)

^

(20)

(21)

(22)

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

dx + еу z V f d2е2 + ~ f = 0,

(23)

45


Рис. 2. Поверхностный годограф т2 от плосконаклоненной границы R.

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

Cos

.(m)'=Cos<pi =

 

COS (jn)

= COS ф2 = —

(24)

46

Cos (kn) — Cos фз = ~)LL—d2__ e2,

 

У f

 

Учитывая

(24), легко получить уравнение плоскости в нор­

мальном виде, поделив (23) на — У /;

(25)

X Cos фі + у Cos ф2 + z Cos фз — Н = 0,

где

_

 

 

 

(26)

является кратчайшим расстоянием от начала координат до плоскости R. В начале координат расположен источник коле­ баний, поэтому Н — аналог величины tü, применяемой в прак­ тике интерпретации временных разрезов:

Н = ^ ,

(260

Рис. 3. Взаимосвязь углов в точке отражения М.

Углы фі, ф2, Фз можно выразить через

а — угол падения и

|3 — азимут падения плоскости R. Если

предположить, что

ось ОХ расположена в меридиональном направлении

(рис. 3),

то

ч

 

 

cos <р2

tg ß.

(27)

cos tfj

 

 

а =

Фз-

 

(28)

Добавляя сюда известное соотношение

 

CoS2 фі + CoS2 ф2 + Cos2 фз = 1,

47


п о л у ч и м

 

 

Coscpi=Sina

Cos ß,

(2 9 )

Cos q>2= Sin a

Sin ß.

 

Знак произведения зависит от величины ß (Sina^O).

Теперь рассмотрим линейный годограф, полученный от из­ вестной плоскости R. При выбранном источнике колебаний продольный профиль должен лежать на прямой, проходящей через начало координат в плоскости ХОУ.

Не нарушая общности, профиль можно расположить по оси ОХ, считая d я е в уравнении годографа тг параметрами.

При любом другом положении профиля поворотом

систем

координат ХОУі и ХОУХ вокруг вертикальных

осей Ot

и OZ можно совместить с ним ось ОХ. В новой системе коор­

динат вид уравнения годографа сохранится, изменятся толь­ ко значения коэффициентов d u e как параметров.

Годограф, наблюденный вдоль оси ОХ, получится пересе­

чением поверхности Т2с плоскостью XOZ:

 

 

т2* = I у х°2 +

2 dx° + f.

 

 

 

V

 

 

Решая

задачу

интерпретации

линейного

годографа тг*

согласно

(5) — (7),

получим границу отражения в виде пря­

мой:

 

 

 

 

 

dx z ] X j d2 -1-

f = 0.

(30)

Эта граница мыслится интерпретатором находящейся в пло­ скости XOZ. Однако истинная граница, полученная сечением плоскости отражения R плоскостью XOZ, имеет вид:

dx z V f — d2— е2 + ~ f = 0.

(31)

Обозначим истинную глубину границы через zu а полу­ ченную через Z2, и образуем разность или ошибку в определе­

нии глубин. Из (30) и (31) имеем:

=

= А ( d x + ± f ) , где Л = 7 -._ 3 |_ е, - і7Ш ^.(32)

Легко заметить,, что öz = 0 в случае е = 0

(Coscp2= ß = 0),

т. е. когда плоскость R параллельна оси ОУ и, следовательно,

профиль наблюдения расположен вкрест ее

простирания

(рис. 4). Это подтверждение известного в практике геофизи­ ческих работ факта.

48