
книги из ГПНТБ / Применение ЦВМ и средств вычислительной техники в геологии и геофизике [сборник]
..pdfРис. 5. Сечение x=const, при расположе нии профиля вдоль простирания R.
Из рис. 5 и (24) видно, что в этом случае
е2 = f Cos2 фг — f Sin2 а.
Подставляя значение е2в (33), получим
6z = # (se a — 1). |
(34) |
Ошибка по глубине бz зависит от угла наклона а и Я. При изменении а от 0 до 90° бz возрастает от 0 до со по закону секанса.
Для выяснения влияния бокового уклонения определим величины ошибок в координатах точек отражения. Пусть пункт приема находится на выбранном профиле в точке
{х°, О, О,}.
50
Применяя формулы |
(5) — (7), получим координаты истин |
||||
ной точки отражения М |
|
|
|
||
Хі |
= |
(x<> + |
d ) f |
d, |
|
2 ( d x ° |
+ f ) |
||||
|
|
|
|||
Уі |
= |
e f |
■e. |
||
2 ( d x ° + f ) |
|||||
— |
|
dl - |
e* + V f — d*— e*. |
Из тех же формул для плоского случая будем иметь
х2 =■ ( х ° + |
d ) f |
2 ( d x ° |
+ /) |
У2 = О |
|
г2 = — |
f ) |
2 ( d x » + |
Вычисляя разности координат при учете (24), (28), (29), получим:
|
Ах = 0 |
|
(35) |
|
Ау = — £ # S in ß S in a , |
(36) |
|||
А2 = — ВН (C osa— У 1 — Sin2ß Sin2а), |
|(37) |
|||
где |
|
|
|
|
В |
2 ( х ° cos ß sin а + Н ) |
(38) |
||
х ° cos ß sin а |
2Н |
|||
|
|
|||
Величина В пропорциональна |
z x и равна І1 при х°=0. Если |
при смещении из точки |
х°=*=0 по |
профилю глубина растет, то |
В становится больше '1, |
если она |
уменьшается — В меньше 1. |
ң |
, т. е. на возможном пересечении R с ли |
||||
В= 0 при х°= |
|||||
нией профиля (рис. |
4). |
При |
В = 1 легко |
проследить |
зависи |
мость Ау и Az от а и ß. |
В целом Ау и Az растут при увеличе |
||||
нии глубины залегания |
отражающей границы и угла |
а, при- |
|||
а, град. |
|
|
Глубина (Н), |
|
|
|
1000 |
1500 |
2000 |
|
|
|
|
|
|||
5 |
|
87 |
130 |
174 |
|
10 |
|
173 |
260 |
346 |
|
15 |
|
258 |
387 |
516 |
|
4*
51
нимая максимальные значения при ß= 90°, т. е. при располо жении профиля вдоль простирания плоскости R.
Величина Ду в нашем -случае характеризует боковое укло нение. В таблице даются значения Ау в метрах при ß = 1.
ВЫВОДЫ
Аналитическое решение прямой и обратной пространст венных задач геометрической сейсмики позволяют вплотную подойти к оценке погрешности аппроксимации решения про странственных задач решениями плоских задач. Решение об ратной задачи для слоисто-однородной среды с плосконакло ненной сейсмической границей дает возможность определить ошибки, имеющие место при построении глубинного разреза по продольным линейным годографам. Анализ этого случая показывает, что
1) ошибки по глубине Az, 6z |
и боковое уклонение сейсми |
|
ческого луча Ау выражаются через величину to, а—угол |
па |
|
дения и ß — азимут падения |
плоскости отражения |
(26), |
(34)-(38). |
|
|
Следует заметить, что все эти величины на практике ча сто бывают приближенно известны, особенно при проведен ной детальной съемке;
2) ошибка в на-праРлении профиля равна нулю (35). Истинная точка отражения может быть смещена в плане от носительно предполагаемой точки в направлении, перпенди кулярном профилю на величину Ау (36);
3)боковое уклонение сейсмических лучей Ау может до стигать значительных размеров даже при небольших углах наклона границы (см. табл.). Поэтому при определении по ложения границы в плане необходимо его учитывать;
4)ошибка по глубине (34) приобретает существенную ве личину при углах а, превышающих 20—25°.
|
|
Л И Т Е Р А Т У Р А |
|
1. П у з ы р е в |
Н. Н. |
Интерпретация |
данных сейсморазведки мето |
дом отраженных |
волн. |
М., Гостоптехиздат, 195Э. |
|
2. Г у р ь я н о в В. Ж. Лучевой метод |
интерпретации годографов сей |
смических волн. Физика Земли. Изд-во АН СССР, № 9, 1965.
В. М. ГУРЬЯНОВ, О. В. КАРЕВА
ТРАНСФОРМАЦИЯ ЛИНЕЙНЫХ ГОДОГРАФОВ ОТРАЖЕННЫХ СЕЙСМИЧЕСКИХ ВОЛН
В работе [1] показана целесообразность представления линейных годографов отраженных волн в виде функции двух переменных: точки расположения источника колебаний и точ ки их приема. Получены дифференциальное уравнение пер вого порядка, которому удовлетворяет упомянутая функция, и полный интеграл этого уравнения. Нами используются ре зультаты работы [1] для трансформации годографа ОГТ в ли нию to и наоборот.
Будем считать, что имеется плоская оплошная |
двухслой |
|
ная среда, отнесенная к декартовой системе координат XOZ. |
||
Вдоль профиля Ro, |
заданного уравнением |
|
|
z = 0, |
(1) |
в точке А 1 {£, 0 } |
помещается сейсмоприемник, а |
в точке |
А 2 { т], 0 } — источник колебаний. Наблюдаемый годограф т отраженных волн в соответствии со сказанным выше пред ставляется в виде поверхности
|
* = |
т(|, л) |
|
(2) |
Эту поверхность можно |
определить по кривой на |
ней, т. е. |
||
по данным |
Коши, и полному интегралу [1,] |
принимающему |
||
в данном случае вид: |
|
|
|
|
Т(£, ті) = |
7 ( ? - Л ) 2 + |
{kl + Ъ) {кц + |
6)*, |
(3) |
где k и b — произвольные константы. Введем новые константы
53
а = |
2k |
с |
2b |
(4) |
|
уП + £2 |
+ ** |
||||
|
/1 |
|
|||
и запишем полный интеграл (3) в виде |
|
||||
'Fid, Л> а, с) = о V (I, |
л) — ( І ~ Л)2 — (al + с) X |
||||
|
X (arj + c) = 0 |
|
(б) |
Трансформация годографа ОГТ в линию to
Расположим начало координат в центре базы суммиро вания по ОГТ. В соответствии с методом ОГТ сейсмоприем ник и источник колебаний должны располагаться симметрич но относительно центра базы суммирования. Годограф ОГТ можно записать так:
£ = — X, л = * , ф(х) = т ( — X, X) |
(6) |
По полному интегралу (5) и данным Коши (6) получим функ цию т (g, л), пользуясь известной методикой [2].
Подставим данные Коши (6) в полный интеграл (5) и ре зультат продифференцируем по х.
Д' (— X, X, а, с) = ѵ2ср2 (х) — 4 х2 + а2х2— с2 = О
— = ѵ \ 2'(х) |
+ 2 (a2 — 4)х = 0 |
(7) |
|
Отсюда получим |
|
|
|
а2 = 4 _ |
ѴУ ' М |
(8) |
|
|
|
2х |
|
|
|
|
|
С2 = ѵ 2 ((р2 - ^ - ) . |
|
||
Подставляя а и с из (8) |
в |
(5), получаем однопараметриче |
|
ское семейство поверхностей |
|
|
|
Mg, л. а (*), |
с (*)] = 0 |
(9) |
|
Огибающая этого семейства и является искомой |
поверхно |
||
стью r = t(g, л): |
|
|
|
^— (а'6 + с') (ал + с) — {а\ + с) X
дх |
(10) |
Х (а'л + с ')= 0 |
Дифференцированием по х и последующим делением резуль татов можно получить
с' |
а 9 |
а' |
с |
54
Используя это соотношение, запишем (10) в виде:
2 асці + сЩ + г]) + ах\а\1 + п) + |
2 с] = 0 |
(11) |
|
Таким образом, |
(5) дает искомую функцию т(|, г)) |
в функ |
|
цій параметра х, |
связь которого с | ит) |
определяется соотно |
шением (11). Для получения линии t0 надо источник колеба
ний и сейсмоприемник совместить, |
т. е. положить в (5) и (11) |
||
| = т1= 2;. В итоге будем иметь по формуле (5) |
|
||
т(£, £) = *о(6) |
= - ( a S |
+ c), |
(12) |
|
V |
|
|
где а и с даются формулами (8), т. е. |
|
||
а» = 4 |
- ^ |
|
(8) |
|
2 х |
|
|
С2 = о2(ф2 - - р )
По равенству (11) получаем, учитывая, что а£+с=ѵт=*0,
£ = - -с* |
2 |
(13) |
Равенства (12) и (13) дают линию ^0(£) |
в функции парамет |
|
ра X и годографа ОГТ ф(х). |
ОГТ — гипербола: |
|
Пример. Пусть годограф |
Ф (х) = У Ьх2 + d
Тогда а и b постоянны
а2 = 4 — Ьѵ2, с2 = v2d
и линия to(Q — прямая:
*о(£) — — (аЪ+ с)
V
Формула (13) в этом случае представляет собой соответ ствие точек на годографе ОГТ и линии і0.
Трансформация линии U (!) в годограф ОГТ
В этом случае задана линия |
т. е. данные Коши |
È = л = Е. *(Е. о = to(Q
Решение получено в [1] и имеет вид:
t(g, л) = - (g- л )2 + m i - O t o ' + to][(4 - l ) t o ' +
V
2(g — « .(л -Е К о ' + М Е — 2 £ + л) = 0
55
Полагая | = —х, г\ = х, получаем годограф ОГТ
т ( - Х , X) = ф (X) = + t о2 (С) - (14)
МО |
(15) |
|
/о'(С> |
||
|
Равенство (15) показывает, что при трансформации линии MS) в годограф ОГТ должно соблюдаться условие
ѵ2_ t ^о(0 |
> 0 , |
|
V (О |
||
|
из которого следует, что не всякий участок линии М£) — может быть преобразован в годограф ОГТ. Более подробно этот вопрос рассмотрен ниже.
Пример. Пусть toil) — прямая.
t o i Q ^ a l + b
Тогда по (15) и (14) имеем
4x2 |
+ |
abl + |
b2)и |
а2)X2 + 6 2р , |
ф М = ( V2 |
;(16)
(17)
т. е. годограф ОГТ представляет собой гиперболу. Равенство (17) показывает, что не вся линия MS) преобразуется в годо граф ОГТ, а только те ее точки £, для которых выполнено условие
ь_ >0
а
Построение годографа ОГТ по заданной границе отражения
При рассмотрении вопроса о влиянии изменения формы границ отражения на годограф ОГТ удобно задать уравне ния границы отражения в параметрической форме с помощью некоторого параіметра ф, а затем через этот параметр выра зить все функции, описывающие линию М£) и годограф ОГТ.
В вышеизложенных выводах все функции зависели |
от пара |
|
метра £ линии t0. |
границы отражения |
являются |
Итак, пусть координаты |
||
функциями параметра ф |
|
|
* = * (ф ),£ |
= 2(ф) |
(18) |
56.
С другой стороны, имеется связь границы отражения с линией /о, а, значит, и с функцией годографа ОГТ [1]
X = } (2 £ |
»**0 (С) 'о' (С)1 |
|
(19)
Z = - о/
2
Приравнивая правые части (18) и (19), получаем неявно вы раженную связь между некоторым параметром ф границы от ражения и параметром £ линии t0, а именно
а д = £ |
ѵЧ0(0у (С) |
(20) |
4 |
а д = - ^ w j / i - ^ p
Здесь £ = £(ф).
Производную от функции іо также возьмем по параметру ф
^ОЙ(Ф)] |
£*0 ^ = К |
( 21) |
|
аф dt j |
|
Точкой обозначено дифференцирование по ф. Приведем (20) с учетом (21) к следующему виду
а д
4£
Z(ф) = і / ѵЧ* ... ?%*<> |
( 20') |
г4 4t
Найдем явную зависимость между параметром £ линии t0nф^ которая позволит связать границу отражения с годографом ОГТ.
Определим параметр £ из системы (120'). Из первого уравнения системы имеем
V2t0t0
4t
Подстановкой этого выражения во второе уравнение находим:
£ = Х ± |
(22) |
Продифференцировав £ в (22) по ф, получим
57
vH 0ta
-Zz
(2 3 )
-Z2
Подставим выражения для Z. и £ из (22) и (23) в первое урав нение из системы (20'):
Х = Х ± |
|
v 4 0t'o |
|
|
|
v 2t0i0 |
ZZ |
|
|
— |
|
|
4 I X ± |
— |
|
|
|
V - 7 |
|
После преобразований имеем следующее выражение |
|||
± v w |
-z‘ = z — |
(24) |
|
|
X |
|
|
Подставим в (20) значение радикала из (24): |
|
||
|
№) г (• |
|
(25) |
SW = * W + ^ 2jiL |
|
Ж )
Это выражение дает искомую связь между параметром грани цы отражения г|з и параметром £ линии t0.
Найдем зависимость t0 = ^ £ W ] из выражения (24)
|
*o[SWl = ^ |
| / |
1 + |
і ж ]2 |
(26) |
|
|
|
|
■X(4) • |
|
Формула |
(26) позволяет получать вид линии tQпри заданной |
||||
границе отражения (18). |
|
|
|
|
|
Чтобы установить соответствие между линией *о[£(ф)] и го |
|||||
дографом |
ОГТ, т. е. функцией ср(х), продифференцируем |
||||
А>[£,Ф)І по £ |
dta |
1 |
|
|
|
|
|
|
(27) |
||
|
|
dfy |
dijdty ’ |
||
|
|
|
|||
где -^jj- |
можно получить |
из (26), а |
-------- из |
(25). В итоге |
приведем схему вычислений по выведенным формулам. Задаются функции Х= Х(<р), Z=Z(ty) (18), описывающие
выбранную для изучения границу отражения. Параметр £ ли-
58
нии t0 и сама функция t0 выражаются через этот параметр по формулам:
т ) = в д |
+ іМ £ Ж |
(25) |
|
■*(4) |
|
u m ) } |
1 + № | 3. |
(26) |
V V |
Ljc(+) J |
|
Используя (25) и (26), вычисляется производная t0 по E
*о'І£(Ф)] |
dt0 |
1 |
(27) |
|
dty |
dl'db |
|||
|
|
И, наконец, для полученных значений £, f0, f0' подсчитываются для данного ф значения координаты х точек взрыв-прибор
*2(Ф) = Е 2(Ф )-Е (Ф ) |
(15) |
и функция ф(х) годографа ОГТ
«p W = {4- ^ + ^ 2К (Ф )]-
Отражение от границы и прием осуществляются при условиях х2^іО в (15) и 2^< —в (27), означающих, что для точек грани
цы, в которых эти условия не выполняются, граница отраже ния имеет области тени для упругой волны, т. е. отражения по методу ОГТ не существует для данного центра суммирования.
Если граница отражения задается в явном виде, т. е. Х=Х'(ф) —ф, то приведенные выше формулы приобретают сле дующий вид
|
Z = Z(X) |
(18') |
||
£(*) |
= X + Z(X)Z'(X) |
(250 |
||
Ш * ) ] = |
|
/ 1 +Z'*(X) |
(260 |
|
ікШ *)] = |
2z(x) |
(270 |
||
K l+ z'2M |
||||
|
|
|
и в формулах (14—15) берутся значения £, t0, t0' из (25'—27'). Рассмотрим конкретные формы границ и их годографы
ОГТ. Пусть граница отражения (18') имеет вид:
Z (X) = аХ + Ъ + а sin (соХ + ѵ), |
’(18") |
59