Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Применение ЦВМ и средств вычислительной техники в геологии и геофизике [сборник]

..pdf
Скачиваний:
8
Добавлен:
19.10.2023
Размер:
4.59 Mб
Скачать

Рис. 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 - ^ - ) .

 

Подставляя а и с из (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) к следующему виду

а д

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

Соседние файлы в папке книги из ГПНТБ