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

книги из ГПНТБ / Богданов, Алексей Иванович. Интерпретация сейсмических годографов

.pdf
Скачиваний:
6
Добавлен:
29.10.2023
Размер:
10.99 Mб
Скачать

Исключая из /2 и <р2

угол /0, находим уравнение типа

 

— /' (|ЗпЯ, рп.г, п).

Из уравнения такого

вида для различных значений п могут

быть рассчитаны и построены в тех же, что и ранее, осях координат фп Я и (3 х) семейства кривых зависимости (3Г1Я от [Зпж с различными параметрами г>орпС По-прежнему для каждого -значения п будет получено только одно семейство кривых с

параметром v0

t.

 

/2 и <р2 для п = 2 и п — 3 сопря­

Исключение параметра i0 из

жено

с большими трудностями

и не является обязательным.

Построить последнее

семейство

кривых можно так. Рассчитав

и построив по /2 семейства кривых зависимости

рпЯ от |Зпт с

параметром г0, по (р2 для тех

же значений г0 и |3ПЯ находим

значения v$nt. Эти значения

должны быть

приписаны в

осях

координат

рпЯ

и рп.т точкам пересечения

ординат рпЯ

скривыми, имеющими соответствующий параметр г0. Вычисляя по ф2 значения г?ор„/ для большого числа величин /0 и рпЯ, покрывающих достаточно равномерно координатную плоскость РпЯ и рп.т, проводим линии равных значений г’орп£ и тем самым строим второе семейство кривых с параметром v0 pn t.

Совокупности семейств кривых с параметрами г0 и v$nt в осях координат рпЯ и рпх будут представлять собой обобщен­ ные лучевые диаграммы. Для каждого значения п будет только одна такая обобщенная лучевая диаграмма.

На рис. 108 приведена начальная часть обобщенной лучевой диаграммы, построенной по изложенному приему и рассчитанной по приведенным выше формулам для п = 2. Аналогичные обоб­ щенные лучевые диаграммы, содержащие пару семейств кривых

спараметрами г0 и tfopn£, могут быть построены для п = 1, п = 3

ит. д. В качестве иллюстрации путей перехода от обобщенных

крабочим лучевым диаграммам на рис. 109 в полулогарифми­ ческом масштабе изображена зависимость средней скорости от глубины для района купола Бекшибай (Южная Эмба), построен­ ная по данным, приведенным в работе Н. Н. Пузырева [38]. Сопоставление ее с кривыми, изображенными на рис. 104,

показывает, что экспериментальной кривой,

приведенной на

рис. 109, лучше всего отвечает кривая,

соответствующая п = 2

и (Зг = 0,0014

на рис. 104. По отрезку,

отсекаемому индексом 1

логарифмической шкалы

наблюденной

кривой

на оси

абсцисс

палетки, определяем, что

l/p2v0 — 0,357.

 

 

Зная р2,

находим v0 = 2000 м/сек, что

соответствует

заданной

величине (г?ср

при

Я = 0

равна 2000

м/сек). Следовательно,

параметры

разреза

в районе купола Бекшибай равны: п = 2,

Р2 = 0,0014

,

v0 = 2000 м/сек.

 

 

 

239

Для определенных значений пг, (З2 и v0 на основании обобщен­ ной лучевой диаграммы можно построить рабочую лучевую диаграмму, которая в осях координат х и Н (в одном масштабе)

О

500

1000

<500

2000 Н

Рис.

109. Кривая зависимость 1g усР

от II для

района купола

Бекшибай

 

Параметры разреза:

(Южная Эмба).

 

 

 

п = 2,

v0 — 2000 м!сек, Ра = 0,0014.

представляет собой два семейства кривых, из которых одно с пара­ метрами г0 (лучи) и другое с параметром t (фронты волн). На

Рис. 110. Рабочая лучевая диаграмма для района купола Бекшибай.

При п = 2, 02 — 0,0014, го = 2000 м/сек.

рис. 110 изображена рабочая лучевая диаграмма для приведен­ ных выше параметров разреза в районе купола Бекшибай, полу-

240

ченная из обобщенной лучевой диаграммы, часть которой пока­ зана на рис. 108.

Перейти от обобщенной лучевой диаграммы к рабочей можно на основании элементарных расчетов. Так, например, взяв на обобщенной лучевой диаграмме, изображенной на рис. 108 и соответствующей случаю п = 2, произвольную точку А, опре­

деляем, что ее абсцисса

$2х = 0,457,

а

ее ордината

(ЗгЯ = 0,98.

Параметры i0 и v0$zt точки А равны:

i0

= 20° и

г?ор21 = 0,9.

По

известным

параметрам для

района

купола Бекшибай

0г = 0,0014—, v0

= 2000 м/сек находим следующие координаты

 

м

t точки А на рабочей лучевой диаграмме:

х w Н и параметр

0,457

 

тг

0,98

 

М’

.

0,9

„ по.

'С ~ 0Ж4 — 325 М’

Н

~ 0^014 ~ 700

1

~2$

~ 0,321 СеК-

Положение точки А на рабочей лучевой диаграмме изображено

на рис.

НО. Параметры г0 этих точек на обеих диаграммах тожде­

ственны. Взяв на обобщенной лучевой диаграмме ряд точек с рав­ ными параметрами i0 либо и пересчитав их координаты и параметры на координаты и параметры их на рабочей лучевой диаграмме, получаем положение лучей с параметрами z0 и положение фронтов волн с параметром t на рабочей лучевой диаграмме. Критерием правильности построения рабочих луче­ вых диаграмм из обобщенных является ортогональность кривых обоих семейств на рабочей лучевой диаграмме при равенстве линейных масштабов по осям х и Н1.

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

На годографе отраженных волн выбирают произвольную точку А с абсциссой ха, ординатой ta и градиентом времени та. В линейном масштабе рабочей лучевой диаграммы строят рельеф местности в интервале профиля, где был наблюден годограф (или совокупность годографов). На чертеж наносят пункт взрыва О и точку А', расположенную на дневной поверхности с абсцис­ сой ха- В точки О и А' своим началом координат накладывают две рабочие лучевые диаграммы, рассчитанные по известным для данного района параметрам п, v0, 0П, так, чтобы их оси Н были направлены вертикально вниз. Если величина ха в точке А изве­ стна, то определяют угол прихода отраженной волны в точку А по формуле

sin va = vora.

1 Рис. 108 и рис. НО были рассчитаны и построены слушателями гео­

физического

отделения Акад. нефт. промышл. А.

И. Николаевским и

В. И. Янко.

 

16 Заказ

1947.

241

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

точки Иа в точку А' и падающей волны

из

точки

О в

точку Ra

должно быть равно наблюденному времени

tn:

 

 

 

 

ta = t (^1 На) -|- t (ORa).

 

 

 

 

 

 

 

По параметрам фрон­

 

тов волн на рабочих лу­

 

чевых диаграммах срав­

 

нительно

быстро опре­

 

деляем положение

точ­

 

ки Ra, сумма

времен

 

пробега волн от которой

 

до точек О и А' удов­

 

летворяет приведенному

 

равенству.

Аналогич­

 

ным путем

может быть

 

найдено

положение

то­

 

чек, лежащих на отра­

 

жающей

 

поверхности,

 

от которых после отра­

 

жения

волны

пришли

 

к разным точкам днев­

 

ной поверхности. Поло­

 

жение отражающей

по­

 

верхности

определяется

 

как геометрическое

ме­

 

сто точек R, определен­

Рис. 111. Схема построения отражающей по­

ных по

параметрам х, t

и

т ряда

точек

про­

верхности при помощи рабочих лучевых диа­

дольного

 

годографа от­

грамм.

 

раженных

волн.

 

 

 

Знать величину т в точках продольного годографа отражен­

ных волн необязательно. Если, например,

т„ в точке А годографа

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

242

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

В США отражающие гра­ ницы раздела, покрытые непрерывными средами, стро­

ят по

параметрам tQ и

т0 го­

 

 

 

 

дографов

отраженных

волн

 

 

 

 

в пункте

взрыва, т.

е. по

 

 

 

 

центральным

лучам

[26],

 

 

 

 

так же, как и

в

случае од­

 

 

 

 

нородных

сред.

Для

этого

 

 

 

 

используют специальные ли­

 

 

 

 

нейки, например изображен­

 

 

 

 

ную

на

рис.

ИЗ.

Таким

 

 

 

 

приемом значительно

проще

 

 

 

 

и быстрее

можно

построить

 

 

 

 

отражающие границы

разде­

 

 

 

 

ла по

сравнению

с описан­

 

 

 

 

ным выше способом. Теория

 

 

 

 

этого

приема сводится к сле­

 

 

 

 

дующему.

 

 

 

 

 

 

 

 

Площадка,

после отраже­

Рис.

112.

К выводу формул для опреде­

ния

от которой

сейсмиче­

ления параметров Но

ий0 отражающей

ский луч, вышедший из пунк­

 

 

площадки.

та взрыва, снова возвращает­

 

 

 

 

ся в пункт взрыва, занимает

 

 

этому

сейсмическому

положение,

перпендикулярное

к

лучу. Угол наклона ф этой площадки к горизонту равен углу 1Н, образованному лучом с вертикалью в точке отражения (рис. 112). Если известны v0, т0 или у0*, параметр этого луча находят из уравнения

г0 =■ arc sin (v0 т0) = arc sin

Определив г0 и зная t0 = 2tH, находим угол наклона ф отра­ жающей площадки в точке отражения центрального луча из фор­

мул (IV.38),

(IV.48) и (IV.53),

которые для случаев п = 1, 2 и

3 принимают

соответственно

вид:

1g*

243

при п = 1

tg_^ = e^A£otgA;

при п = 2

 

60 Zn

 

 

Ч> = 4

sm;0 + г0;

 

при п = 3

 

 

 

cos ф = cos г0---- v° to sin2 i0.

 

По этим формулам могут быть составлены три таблицы, при

помощи которых по

известным

значениям угла i0 или

г?ото и

произведения ro0nio

находят угол наклона ф отражающей пло­

щадки.

 

 

 

Положение же точки отражения может быть определено или

при помощи ее координат х и Н,

или длины отрезка Но,

который

отсекает перпендикуляр, восстановленный к площадке в точке

отражения, на оси Н и длины этого

перпендикуляра Ro между

точкой отражения и точкой его пересечения с

 

осью Н (рис.

112).

Выражения для х, Н, Ro и Но

для п = 1,

2,

и 3 могут

быть

записаны в следующем виде.

 

 

 

 

 

 

 

 

При п = 1

 

 

 

1

 

cos г0—COS ф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X =

д—----- ,

 

 

 

 

 

 

 

 

Pl

 

 

Sin 10

 

 

 

 

 

 

 

 

 

_

 

1

sin ф — cosi0

 

 

 

 

 

 

 

~ 0!

 

sin i0

 

 

 

 

 

D

х

 

'

 

1

cosi0 — созф

,

 

 

 

 

/1 Л — .

ф

 

о

.

 

 

 

 

и

Sin

 

 

pi

Sin

l0

Sin Ф

 

 

 

 

И, = Н - ж ctg ф _ *

 

 

 

 

.

 

 

 

 

 

 

 

 

Pl

sln lo sin Ф

 

 

При n = 2

_

1 (2ф—2ip) + sin 2i0 —sin 2ф

 

 

 

X

 

 

 

02

 

 

 

1 — cos2i0

 

 

 

 

_ 1

sin2 ф—sin2 ip

1

cos2i0—ссз2ф

 

02

 

sin2

ip

 

 

 

02

 

1 —cos2i0

 

_

x

_

1

 

(2ф — 2i0)-]-sin 2i0 — sin 2ф

 

0

sinip

02

 

 

зтф(1 — cos2i0)

 

Ho - H - X cig <p = 4-

 

siny—sln<2io—p) —(2Ф—

 

0

 

T

02

 

 

31Пф(1—cos2lp)

 

При n = 3

 

 

 

 

 

 

 

 

 

 

 

 

 

_

1

9 cos ip—9 cos ф-f-cos Зф — cos3i0

 

X

03

 

 

3 sin ip—sin3i0

 

 

 

244

_ 1

3 sin <p—sin3q> — 3 sin i0-|-sin 3i0

 

Рз

 

3 sin i0—sin3i0

x

_ 1

9 cos i0— 9 cos <p-|-cos 3<p—cos3i0

sin <p

p3

sinq)(3sini0—sin3i0)

 

IIQ = II — x ctg <p = -Д- X

Рз

1 -{-cos2 cp-j-cos (3i0— <p)-|-3 [sin2 i0-|-(cos i0 — cos cp)2 — cos (<p — i0)] sin <p (3 sin i0 —sin 3i0)

По приведенным формулам для различных значений п могут быть рассчитаны и составлены таблицы, па основании которых по известным значениям i0, ср и (Зп

могут

быть определены

величи­

ны Ro

и Но.

Зная

Ro

и Но,

строим

площадку,

отвечающую

центральному

лучу,

при помощи

линейки, изображенной на рис. ИЗ.

Эта линейка

представляет

собой

 

 

 

неподвижный лимб, центр кото­

 

 

 

рого устанавливается

в точке с

Рис. 113. Линейка для

построе­

ординатой Но

под пунктом

взры­

ния отражающих

площадок (по

ва в масштабе, принятом для по­

 

 

 

строения разреза.

 

имеется алидада, которую

 

На неподвижном

лимбе

можно

повернуть относительно центра лимба на угол, равный найден­ ному значению ср. В алидаде имеется подвижная линейка, поло­ жение конца которой относительно центра лимба и алидады может изменяться при ее передвижении. Эту подвижную линейку можно оцифровать в миллиметрах, принимая за начало отсчета ее конец, и, перемещая ее в лимбе, отложить в масштабе, принятом для построения разреза, найденную величину Ro (от конца линейки до центра лимба и алидады). Линия, проведенная по короткому концу линейки, и дает искомое положение отражающей площадки, отвечающей центральному лучу с параметрами t0 и т0. Милли­ метровые деления шкалы подвижной линейки позволяют ее использовать для построения отражающих площадок в случае изменения скорости с глубиной при любом значении п.

Точность определения элементов залегания отражающей пло­ щадки по приведенным выше формулам зависит только от точности определения величины i0 или, что то же самое, от точности опре­ деления т0 или vQ,* если считать, что п, (Зп и v0< определены точно.

Поскольку величины т0 или v0* и, следовательно, i0 не всегда могут быть определены достаточно точно по наблюденному годо­ графу отраженных волн, элементы залегания ф и Hz плоской

245

отражающей площадки при известных значениях п, (Jn и г?0 можно установить по параметрам tm и х0 наблюденного годографа отра­

женных волн, которые

находят

с

большей

точностью,

чем т0,

*v

и г0 (tm и ,г0

— ордината и абсцисса

минимума

наблюденного

годографа).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Действительно, при любом значении п для

луча,

проходящего

после

отражения

к дневной

поверхности по

вертикали, имеем

(рпс.

112)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно,

при

п■

 

“ 2ф-

 

 

 

 

 

 

 

— 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

cos i0'—cos 2ср

 

 

 

 

 

 

 

 

 

■/(|~ 0!

 

sin г0'

 

 

 

 

 

 

 

 

 

• г ,

 

1

sin 2<р—sin !0'

 

 

 

 

 

 

 

 

 

0

— Pi

 

 

sin i0'

 

 

 

 

 

 

 

 

tm — t\

t2 =

о

III------------- -

. -7Г- 1 0 (1

+ f>iH0 )

=

 

 

 

 

4o Pl

‘o

'

о Pi

 

 

 

 

 

 

 

 

 

 

 

 

ё

2

 

 

 

 

 

 

 

 

 

 

 

=

1

In

sin2 ф

 

 

 

 

 

 

где

x0

и H'o — координаты

точки

отражения;

i'

— угол

выхода

из пункта взрыва луча,

падающего после отражения нормально

на дневную поверхность;

 

ti и

fa

— времена пробега

волны вдоль

падающего и отраженного лучей.

 

 

 

 

 

 

 

 

Решая совместно уравнения для х() и

tm,

находим независимо

значения углов

i'

и ср через

известные

значения- х0, tm,

и »0.

Так,

например,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ts

 

=

 

ж° Р1

 

 

 

 

 

 

 

 

 

 

ё 2

e”i) iim__j ’

 

 

 

 

 

а

”о Pi tm

Подставляя полученные значения углов г0' и ср в уравнение для //0', получим

 

е. 2 |/”

роР1/т

__ _________________

 

1

 

 

1 (®оР1)“

 

о

 

g^oP1^ш

 

 

На основании этих

формул по

известным значениям

х0, tm,

vQ и рх можно найти значения ср

и

Я',

достаточные для

опреде­

2W

ления положения плоской отражающей площадки. Решение может быть выполнено аналитически или графически при помощи семейств

кривых

(номограмм),

зависимости ср

и Н'п

от

vopi£m, для

различных значений

жфх.

 

 

 

 

 

При

п — 2

имеем

 

 

 

 

 

 

 

_ _1_

(2ф—iu')—-sin (2ф—— Zo') cos (2<р-Но')

 

 

’ 0

Р2

 

 

 

sin2 i0'

 

 

 

 

и

_ 1 sin (2ф~Но') sin (2cp—10')

 

 

 

 

0

 

fh

 

sin2j0'

 

 

tm - tj. + i2 = 2У~‘°?

 

 

 

 

1

4- -4- [О + Р2Я0') 2 ~ 1 ] =

 

 

 

Vq

p2 sin

г?ор2

> rz

U /

j

 

~ v0 p2

'

 

 

 

sin ia'

 

 

Так как выражения для х0

и tm являются трансцендентными,

то выражение для углов i0' и ф в простой аналитической форме получить нельзя. Однако по этим уравнениям можно рассчитать и построить в осях координат i0' и ф семейства кривых равных значений р2х0 и 27op2Z,n, на основании которых по известным

значениям (32 ,т0

и v0 Р2 tm находят значения углов г0' и ф. Опре­

делив углы i0'

и ф, нетрудно определить Но' аналитически или

графически по

заранее рассчитанным и построенным семействам

кривых равных значений (32W(/ в осях координат i0' и ф. При п ~ 3 найдем

_

1

cos i0' (2—[—sin2 i0') — cos2<p (2-J-sin2 2<p)

 

 

У ,

_

1

эт32ф—sin3 i0'

 

 

°

_

Рз

sin3 ip'

_ .

_

6

 

2 cos i0'+ 2 cos 2ф +sin2 2ф —sin2 i0'

-tr±t2- —.

Аналитически решить эти уравнения относительно ф и Гп

трудно. Поэтому для п =

3 задачу удобнее решать графически

тем же путем,

который

был

изложен

выше для п = 2.

По приведенным выше формулам можно рассчитать три та­ блицы для определения углов ф по известным значениям угла го и величин vofinto Для трех значений п (1, 2 и 3), кроме того,

шесть таблиц для определения $пНо и

рпхо по известным значе­

ниям углов io и ф для трех значений п (1, 2 и 3).

Результаты

расчетов

углов ф для

различных значений г0,

^орп^о и

трех

значений

п сводят в три графика зависимости ф

от i0 для

постоянных значений г?орп£о.

 

247

Результаты расчетов значении $пН0 и Р?1 //0 сводят в шесть графиков зависимости ср от i0.

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

Наличие такой зависимости ф/io от гоР« to, $пНо и $nRo при­ водит к следующим интересным и важным для практики выво­ дам, позволяющим значительно упростить количественную ин­ терпретацию наблюденных линейных продольных годографов отраженных волн для сред, в которых скорость нарастает с глу­ биной по рассматриваемому закону.

 

1.

По известному

значению величины v0 fl,, t0

сразу же находят

отношение

ф/г0 и, зная i0, определяют угол ср.

 

 

2.

По

исходному

значению

t’o|3nio непосредственно получают

величины

$пН0 п

Ro. (Для

случая

п = 1

эта возможность

была доказана ранее строго аналитически.)

 

 

 

Тем самым доказывается возможность аппроксимации фрон­

тов

волн

в рассматриваемой непрерывной среде окружностями

с координатами центра х = О,

Н = Но и

радиусами R = Ro,

являющимися функциями времени

to при постоянных значениях

Vo,

 

и

п1.

 

 

зависимости угла ф от угла i0

 

Наличие практически линейной

позволяет построить три новых графика,

дающие семейства кривых

зависимости v0 [1Л t0,

pn//n и ptl Ro

от

<p/i0

для значений n — 1,

n =

2

и п = 3 соответственно (рис.

114,

115,

116).

 

Графики зависимости vofinto,

$nHo и $nRo от ф/io можно при­

менить при интерпретации наблюденных годографов отраженных и преломленных волн в случае непрерывных сред, покрывающих границы раздела. Их также можно использовать для приближен­ ного построения обобщенных и рабочих лучевых диаграмм при п = 2 и 3. Так, например, для заданного значения vofinto по гра­ фику, соответствующему определенному значению п, находим ординату $пНо центра и радиус $nRо изохроны, а также «отно­ шение ф/io = К. Тем самым на изохроне находится точка ее пе­ ресечения с лучом, вышедшим из пункта взрыва под любым углом io, которая определяется как точка пересечения изохроны с ра­ диусом, образующим с осью $пН угол ф, равный Kio.

Эти же графики можно применить для определения угла на­ клона ф п параметров Но и Ro отражающей площадки по углу выхода io и значению io центрального луча при известных пара­ метрах vo, и п непрерывной среды.

Глубину залегания плоской горизонтальной преломляющей поверхности под пунктом взрыва, если опа покрыта непрерыв­

ной

средой и скорость изменяется по рассматриваемому закону

с известными значениями п,

и v0, можно определить по отрезку

io,

отсекаемому прямолинейной ветвью годографа преломлен-

1

Па эту возможность впервые

было указано А. А. Николаевским.

248

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