Файл: Применение ЦВМ и средств вычислительной техники в геологии и геофизике [сборник]..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