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

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

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

профилей, или кажущимися ско­

ростями,

и обозначаются

через

 

* _

 

____

 

(11.45)

 

V+ ~ sin (i — ip)

 

 

 

v* = ____ .

 

 

(11.46)

 

 

sin (i 4-ip)

 

 

 

Если

профиль

идет

вкрест

простирания

границы

раздела,

выражения для

и v*

прини­

мают вид:

 

 

 

*V

=

Г1

*V

=

V1

+

sin (i—<р) '

 

sin(i-|-<p)

 

Выясним возможные

пределы

изменения кажущихся скоростей

распространения фронтов

прелом­

ленных

волн вдоль

дневной по­

верхности.

 

 

залегаю­

При

горизонтально

щей границе раздела

угол

= О

и V.* — = —Л1— = v2.

 

 

 

sin

I

z

 

 

 

По

профилю,

идущему

в

сто­

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

уменьшается и становится равной

vi при условии, что

i

ср = 90°

или <р = 90° — i.

 

 

Если же

угол

ф > 90° — г,

сейсмический

луч,

 

преломив­

шийся на границе раздела, не сможет выйти на дневную поверх­ ность. При этом условии в направ­ лении падения границы раздела не может быть зарегистрирован при­ ход преломленных волн (рис. 42, д). Таким образом, v* может изме­

няться

ОТ V2 до

V1.

По профилю, идущему в сто­

рону

восстания

преломляющей

границы раздела, кажущаяся ско­

рость

*г>

по мере

увеличения

угла

ф

становится

больше гг и

достигает

максимального зиаче-

Рис. 42. Ход лучей преломленных волн при различном соотношении углов in <р и значения кажу­ щихся скоростей распространения фронта преломленных волн вдоль дневной поверхности по обе сто­

роны от пункта взрыва.

а — г > ф,

i -|-ф < 90°;

б — г = Ф,

г -)- <р < 90°;

в — i <ф, '{ ф- Ф < 90°;

г — i < ф,

i +ф = 90°;

д — г < ф,

i + ф > 90°.

7*

99

яия, равного -l-oc, когда угол <р равен углу i (рис. 42, б). При этом значении угла i v’+ меняет знак с плюса на минус и при даль­

нейшем увеличении угла ср увеличивается от —оо до значения, равного — 7^27 , ПРИ Ф = 90° — i (рис. 42, г).

Возможные соотношения углов i и ср и соответствующие им пути распространения преломленных волн изображены на рис. 42. На рис. 43 приведен график изменения кажущейся ско-

Рис. 43.

Кривая изменения кажущейся

скорости v*

распространения фронта

преломленных

волн вдоль

дневной

поверхности

в

зависимости от знака и ве­

личины

угла наклона

преломляющей границы (при

 

i

const).

 

рости *v от величины и знака угла ср. Здесь,

так же как и ранее,

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

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

<. = <«±^±£ = <«±

(П.47)

V

где

I = ]/ ж2 у2.

Это уравнение двух прямых, отсекающих на оси времен от­ резки, равные to. Прямые образуют с осью I угол, тангенс кото­ рого равен ±1/г?2, т. е. кажущаяся скорость ветвей линейных

100

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

Линейный непродольный годограф преломленных волн. Ура­

внение годографа преломленных волн по линейному непродоль­ ному профилю может быть получено также из уравнения поверх­ ностного годографа преломленных волн. Например, профиль па­ раллелен оси Ох и отстоит от нее на расстоянии у = d. Тогда уравнение непродольного годографа можно найти из выраже­ ния (11.20). Положив в нем у — d, имеем

tx = t0 +

4/z2 (ж2 + cZ2) — (хох

y0d)2 —

 

— (V + Уо^) ctg i],

(II. 48)

Из рассмотренных выше свойств изохрон поверхностного годо­ графа преломленных волн следует, что линейные непродольные годографы этих волн должны иметь минимумы.

Для нахождения абсциссы хмин и ординаты (ZX)M0H непродоль­ ного годографа преломленных волн напишем выражение для

значения тх в

любой точке этого годографа.

Для этого в уравне­

нии (11.23) заменим у = d. Тогда

 

 

 

 

хх = -г}- Г

 

 

+

 

_ Хо ctg Я

(П> 49)

 

 

J/ 4/г2

(а;2 —|—<Z2) — (хох-\-уй<1)г

 

J

 

Приравняв в этом

уравнении выражение, стоящее в

прямых

скобках, нулю и решая его относительно х,

найдем значение хмив-

После преобразований получим

 

 

 

 

 

 

4/1

Хо

'л +

у 4^2--- Ж()2

sec2 t

c|8z ■

<11-5о>

 

 

 

 

Подставляя это выражение в уравнение для

времени

£х, на­

ходим после

соответствующих преобразований следующее выра­

жение ДЛЯ (Zx)mhh:

 

 

 

 

 

 

 

^мин =

 

 

/4/г2 - *о( 2 + уо2)

V 4^2

— х2 sec 1 +

 

+ 2//ctgi[4/i2 —(xo2 + d?/o)]}.

 

(П.51)

Когда линейный непродольный профиль,

идущий вдоль оси Ох,

задан вкрест

простирания плоской преломляющей поверхности,

уравнение непродольного годографа примет вид (полагая в (II. 48)

уо = 0):

tx = t0 Н—— (|/ х2 cos2 гр + d2 — ж sin ср ctg г),

(II. 52)

vz

так как в данном случае

х0 Sin(? = Th-

101

Абсцисса

минимума

годографа

хМ|1н

будет

равна

(полагая

в (11.50) уо = 0)

 

 

 

 

 

 

 

Ордината

минимума

годографа

(QMnn

 

будет

равна

(полагая

в (11.51) уо = 0)

 

 

 

 

 

 

 

 

, .

d

V 1—sin2 <р

sec i

.

*2/

.

/тг ~/\

 

Ммин=— --------- ^7--------- +

7TCtgl-

(ILa4>

 

 

с 2

vvd

 

 

с 2

 

 

......... 1 -------2 --------3

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

 

Параметры разреза: tjj = 1780 м/сек, в. = 4000 м/сек,

h = 261 м, <р = 16,7°. Снос

профиля:

d = 750 м.

 

 

 

 

 

 

прямых, отраженных и преломленных

1,2,3 — непродольные годографы соответственно

 

 

 

 

 

 

 

 

 

волн.

 

 

 

 

 

В точках: Aj — xj = -4-750 м,

Ц = 0,412 сек.,

Tj — 22 •

10“6 —— ; А.„,„ — х,,„ =

 

 

 

л

*

*

*

* сек. г

мин

мин

•= 630 м,

— 0,411 сек., тМ1.„ = 0; Ао — х0 = 0,

 

(о = 0,449 сек., то = — 138-10~6 сек.;

Д2 —

х2 = —750 м,

ti = 0,628 сек.,

т = — 298 •

 

10~6 ——

 

 

*•

 

гиг

 

 

г

 

 

орт*

 

 

Выражение для тх

 

принимает вид (полагая в (11.49)

уо — 0):

1

/

X

cos2 <р

— sin ip ctg i

(11.55)

тх =---

\

Ух2 cos2 <р—|- с£2

V2

 

 

 

102

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

Если линейный непродольпый профиль параллелен оси Ох и задан вдоль простирания плоской преломляющей границы раз­ дела, то уравнение годографа преломленных волн на этом про­ филе можно получить из уравнения (11.48), если в нем положить хо = 0. Тогда имеем

t = t0

-f- — (l^x2 + <Z2 cos2 <р

dsin фctg i),

 

(II. 56)

так как

Г2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin<p = -g-,

 

 

 

 

 

Абсцисса минимума хМин годографа

будет равна пулю, если

положить в (11.50) хо = 0. Ордината

 

минимума

£МИц

годографа

будет равна (£х)о,

т. е.

 

 

 

 

 

 

 

I.....=

 

и 2 ЫН с

=А,+

Vj

 

(И.57)

 

 

 

 

 

 

 

Если в (11.49)

положить

хо — 0,

выражение

для

хх

примет

вид:

 

 

 

 

 

 

 

 

 

тх = — ,

COS2 ф

.

 

 

(П. 58)

 

v2

|/ж2-|-(/2

 

 

 

 

Для плоской и горизонтальной преломляющей поверхности уравнение непродольного годографа преломленных волн можно определить из уравнения (11.47). Для профиля, расположенного параллельно оси Ох, при замене в (11.47) у = d получим

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

нены к горизонтальной оси под

углом, тангенс которого ра­

вен ±1/г>2. Время в минимуме

гиперболы

= <о+“.

(П.вО)

г2

 

Годограф преломленных волн на дуговом профиле.

Уравне­

ние проекции дугового годографа преломленных волн

на пло-

103

с кость

Юх можно получить из уравнения поверхностного годо­

графа

преломленных волн, если

положить в нем

 

.г2 + //2 = I2 и у

12 — х2,

где I — радиус дугового профиля.

Взяв за исходное уравнение (11.20), имеем

 

 

 

 

 

 

cos i

 

 

 

(11.61)

 

 

 

 

 

 

sin i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это

уравнение

можно

 

 

 

 

 

 

изобразить

графически в

 

 

 

 

 

 

виде

замкнутой

кривой,

 

 

 

 

 

 

являющейся проекцией на

 

 

 

 

 

 

плоскость Юх сечения по­

 

 

 

 

 

 

верхностного

годографа

 

 

 

 

 

 

преломленных

волн

ци­

 

 

 

 

 

 

линдрической

поверхно­

 

 

 

 

 

 

стью с радиусом /исосью,

 

 

 

 

 

 

совпадающей с

осью

Ot.

 

 

 

 

 

 

Крайними

точками

 

этой

 

 

 

 

 

 

кривой будут точки

с

абс­

Рис. 45. Проекция

дугового

годографа

 

циссами х = ± I.

 

 

преломленных волн

на плоскость xOt в

 

Если

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

том случае,

когда

ось Ох совпадает с на­

 

что ось Ох

идет

в направ­

правлением восстания преломляющей гра­

 

 

лении

восстания границы

 

ницы раздела.

 

 

уравнении

(11.61)

уо = 0 и заменив

раздела, то, положив в

xo!2h

 

через sin ср,

полу­

чим

 

 

 

 

 

 

 

 

 

 

 

 

 

t —

+ — (X Z2 — .т2 sin2 ср — г sin ср

)

=

 

 

 

 

 

 

г2

'

 

т sin

I

 

 

 

 

 

= t0 4-

(sin i

|/ I2 — a;2sin2cp —x sin cp cos i).

(11.62)

Функция, заданная уравнением (11.62), изображена на рис. 45.

При х = -f-Z

t

будет

иметь минимальное

значение, которое

обозначим через

/мин, а при х = — I

t

будет иметь максимальное

значение,

которое обозначим через

ZMaKC.

 

 

 

 

 

 

Полагая х = -j-Z, получим следующее выражение для ZMnn — /0:

1 мин — ZM11H — t0 = (sin i cos cp — sin cp cos i). =

= _L sin (Z - cp).

(H.63)

 

104

Полагая х = —I,

получим выражение

для £макс— to-

7’макс = макс*

 

V1

(s'll I COS ф -f- sin ф COS j)

-

 

=-J-sin (i-f-ф).

 

(11.64)

При x = 0 уравнение (11.62) принимает вид:

 

7’о = to

— to =

SJ1 i =

,

(11.65)

 

 

 

1 ]

1 2

 

где to' — значение t

при

x — 0.

 

 

§ 2. ОПРЕДЕЛЕНИЕ ЭЛЕМЕНТОВ ЗАЛЕГАНИЯ ПЛОСКОЙ ПРЕЛОМЛЯЮ­ ЩЕЙ ПОВЕРХНОСТИ И СКОРОСТЕЙ РАСПРОСТРАНЕНИЯ УПРУГИХ ВОЛН В ПОКРЫВАЮЩЕЙ И ПРЕЛОМЛЯЮЩЕЙ СРЕДАХ

ПО ГОДОГРАФАМ ПРЕЛОМЛЕННЫХ ВОЛН

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

Если, например, ф < i, то изохроны поверхностного годографа представляют собой семейство эллипсов. Направление большой оси семейства эллипсов дает направление восстания или падения преломляющей поверхности. Поверхность восстания имеет напра­ вление в сторону смещения центров эллипсов относительно начала координат (пункта взрыва).

По любой изохроне с параметром tn или Тп = tn — to могут быть определены: ее большая полуось ai, малая полуось bi, полуфокусное расстояние С’г и смещение центра хг.

Из выражений для ai, bi, Ci и xi (см. главу II, § 1), следует,

что

_

С\

sin ф

_ а?!

cos i

61

аг

sin i ’

С,

cos ф '

Из совместного решения этих двух уравнений могут быть най­ дены значения углов ф и t по формулам

 

• •

1/

1 —V

 

Sill Г— У

 

suitp =

л Г

1 — dj2

J/

 

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

 

 

 

х 1/

1—6i2

 

COS <Р =

COS I - 6 1/

-71Гд-<-2-

105

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

г и

ф по

любому из выражений (11.28)

или

(11.29), определяем значение

г?г:

 

 

 

 

 

 

 

 

 

 

Г2

"L

I

1 —V

 

 

 

(П.67)

 

 

 

тп

) 1

— д1ае12

 

 

 

Зная значение угла i

и значение скорости г?2,

определяем зна­

чение скорости

п по

формуле

 

 

 

 

 

 

 

 

 

г _ .1Л 1

^i2

__у о ei2) и ^i2)

 

(11

ро\

 

1

2 У

1—Ve?

 

Тп

1— Ve?

 

'

8)

Последнюю неизвестную величину /г можно найти на основа­

нии

формулы (11.13) из

выражения

 

 

 

 

 

 

 

 

 

 

 

 

 

г-'^о

 

 

 

 

(11.69)

 

 

 

 

 

2 cos i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

которое после подстановки вместо vi и cos i

их значений, получен­

ных выше, принимает

вид:

 

 

 

 

 

 

 

 

 

 

 

/

 

Mi

1/

1 —V

 

 

 

(11.70)

 

 

 

” 2di7’n У

1—Ve,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Параметры

разреза могут

 

 

 

 

 

 

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

также

и

в

 

 

 

 

 

 

том

случае,

если

задано не

 

 

 

 

 

 

семейство изохрон, а лишь

 

 

 

 

 

 

абсциссы четырех

точек

пе­

 

 

 

 

 

 

ресечения любой

изохроны

 

 

 

 

 

 

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

годографа

с

 

 

 

 

 

 

координатными осями хОу.

 

 

 

 

 

 

 

Пусть даны

абсциссы че­

 

 

 

 

 

 

тырех точек пересечения осей

 

 

 

 

 

 

координат

 

изохроной

 

с

 

 

 

 

 

 

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

tn

Рис. 46. Изохрона

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

го­

или

Тп

-- tn

— ta

(рис.

46).

Из уравнения (II.14)-следует,

дографа преломленных волн с пара­

что

 

 

 

 

 

 

 

xs,

метром Тп.

точек пересе­

 

 

 

 

 

 

 

ylt у2 — абсциссы

 

 

 

 

 

 

 

 

чения

изохроны с

осями

координат.

tn — t0^

 

 

sin (i — ф),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

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

плоскости,

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

и

проходящей через

начало

координат

и

через точку с координатами х и у.

 

 

 

 

 

 

 

 

Для точек, лежащих на оси Ох,

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

в следующем виде:

 

 

 

 

 

 

 

 

 

 

 

tn — t0~ — sin (i — фЗ,

11

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

106

Угол грх в сторону положительного направления оси Ох счи­ таем положительным и в обратном направлении — отрицатель­ ным.

Если tn — to = Тп, то можем записать

I т

1

= ——п 1

1

2 !

sin (t -j-iE) '

Соответственно для точек, лежащих на оси Оу:

7 гЦ'1

sin (i —ip2) ’

(II. 72)

где ip2 — угол наклона преломляющей поверхности в плоскости, ей перпендикулярной и проходящей через ось Оу.

Для отношения сумм, разностей и произведений абсолютных величин отрезков xi и хо и отрезков yi и уо соответственно получим

Л =

1

^1 1

+ 1

*2_1_

_

COS 1р!

 

sin2

i—sin2 ipa

 

(11.73)

1

У1 1

+ 1

3/2 1

 

COS 1р2

sin2

i — sin2

ipi

 

 

 

 

к

 

+1 |

— | Х2 |

 

sin ipi

sin2

i — sin2

ip2

(И.74)

 

 

1

У1 1

— 1 У-2 1

 

sin ip2

sin2 i—sin2 ip!

 

 

 

 

 

 

 

Х-1%2

_

sin2 i -sin2 ip2

 

 

(П.75)

 

 

 

 

— У1У2

 

sin2 i -sin2 ipj

 

 

 

 

 

 

 

 

 

 

Следовательно, можно записать, что

 

 

 

 

 

 

 

А

 

 

cosily

В

 

 

sin ipj

 

 

 

 

 

 

С

 

cos тр2 ’

С

 

 

sin 1р2

 

 

 

Решая полученные уравнения относительно грц и ip2, находим

 

 

 

 

 

 

 

 

В л / С2 —А2

 

 

 

 

 

 

 

 

S1nip1= —

 

 

 

Т2,

 

 

 

 

 

 

 

 

. .

 

г

..

 

 

 

 

 

 

(11.7(5)

 

 

 

 

s

 

 

С2 — А2

.

 

 

 

 

 

 

 

 

 

nip2 - |/

 

 

 

 

 

 

 

По значениям углов ipx и чр2 нетрудно определить остальные

параметры разреза.

Так,

из

(11.73)

находим угол г:

 

 

.

 

 

.

 

.

Я/(С2-Л2) (В2— су

 

 

,тт 77.

 

 

S1UI

V (в2—л2)

(с2—с) •

 

 

( • )

Зная углы

ipx,

 

1р2,

i,

получим vr

 

и

г2. После

подстановки

в уравнение (11.71)

 

значений

углов t

и 1рх

имеем

следующие

107

вы ражения:

 

для г\

 

 

 

 

 

 

 

 

 

 

,,

_

Xl

1 /

-- -1)*

/ - 1

I

/ /^2 _

(J__ IL 1

Г |2 _

(Л ■ /II 7О\

tl~

тп

V

*(Н —Аг)(Сг—С)

\ С

с

С

1

<J,(1J./8)

 

ДЛЯ V2

 

х.

1

 

----------

В

 

--------- \

 

 

 

 

 

I

 

 

 

г‘ “ 77

 

(- I

 

«2 - С

-

-I2 - С) (П. 79)

и,

наконец,

для

 

 

 

 

 

 

 

 

" =

 

V (В‘-л^-С} (4 >

 

- 4 к _4^с).(11.80)

Истинный угол падения и азимут восстания а границы раздела относительно оси Ох по известным фх и фг могут быть определены на основании уравнения (II.1).

Действительно,

sin фх = sin <р cos а,

sin фг = sin <р cos (90° — а) = sin ср sin а.

Возводя оба выражения в квадрат и складывая их, получим

sin2 ср = sin2

sin2 фг

 

или

 

 

sin ср = р sin2

-f- sin2 a|?2.

(11.81)

Азимут падения может быть определен из выражения

Когда угол <р мал и в пределе стремится к пулю, семейство изохрон стремится к семейству концентрических окружностей с центром в начале координат. В этом случае даже небольшие ошибки в определении параметров a, b, С, х эллипсов или отрез­

ков xx,

Х2, yi и г/2, отсекаемых изохроной на осях координат, при­

водят,

как показывает анализ,

к очень

большим погрешностям

в определении всех параметров

разреза,

за

исключением ско­

рости

V2. Выражение же для

скорости

V2

при этом стремится

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

108

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