Файл: Невский, М. В. Квазианизотропия скоростей сейсмических волн.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 22.10.2024

Просмотров: 82

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Горизонтальная проекция отрезка луча в k-м слое на линию профи­ ля наблюдения определяется формулой

ЛХк 2Hktg0k

Как было показано, сейсмический луч отраженной волны для рас­ сматриваемой модели среды симметричен относительно перпендикуляра! проведенного в точке отражения к отражающей границе. Поэтому для получения полного времени пробега по лучу от источника до приемника следует просуммировать времена ЛТk и полученную сумму удвоить. То же следует сделать и для величин АХ ^ В результате получим

п211,

Т- 2 _____________

k=l vk(0k)cos0k '

п

X = 2 2Нк tg 6>к.

к = 1

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

через эффективные упругие параметры V^pk, крк, *SHk’ >Ц_к и ^ ки угол #к , что создает трудности при расчетах. Во-вторых, в указанны формулы входят значения углов преломления в каждом слое вк, связь1 между которыми, как показано, достаточно сложна. Для получения уравнений годографа в форме, удобной для вычисления, следует прежД

всего

перейти к параметрическому представлению для vk и #к по фор

мулам

(2 .1 1 ).

 

 

Подставляя в последние уравнения для X и Т параметрические выр

жения для vk

и #к через ik и замечая, что углы

ik по

формулам

(3 .5 ),

(3 .6 )

выражаются через параметр р, т.е.

iк=iк(р),

оконча­

тельно получаем параметрические уравнения годографа отраженных вd в виде:

Т(р) = 2

2

Нк |у^ 1 к (р)Ь v \ [ik(p)]l sec ik(p)+ arclg W

p) ] ,

 

к =

 

 

 

v k[ik(p))

 

n

 

 

Vk[ik(p)]

 

 

X<p)' \

s

1 Mk‘e |iltp)+arc‘e vkTi3 P T f

(3 .8 ).

 

 

 

где w

 

dVk ,

_

sin i!

 

 

 

“di- 4

’ P -

V i(ij)

 

 

82


Последние формулы дают параметрическое представление годографа отраженной волны для случая многослойной поперечно-изотропной по­

dVk

= 0, т.е. покрывающая толща является

крывающей толщи. Если

di

 

изотропной, из формул (3 .8 )

получаются известные уравнения годогра­

фа отраженной волны для изотропной покрывающей среды.

Формулы (3 .8 ) справедливы для всех типов волн в поперечно-изо­ тропных средах: Р, SV и SH Различие состоит лишь в том, что для вычисления их годографов следует выбирать соответствующие зависи­

мости V]^(i), i ^(р)

и значение параметра р Для случая поперечных

волн SH

уравнения

(3 .8 )

можно существенно упростить,

используя

формулы

(2 ,2 3 ) и (3 .5 ).

В самом деле, выражение для v

и ^§Ик

в функции параметра р по указанным формулам можно после некоторых

преобразований представить

в виде

 

 

 

 

 

 

2

2

2

г. . ..

Pv|SK1 1 + (* S H k -1) [ ] - P 2vJ.SkKSHk

tg 0SHk[ *k<P} ] " “ = -----------------------------------------------------

 

 

(k2 _ l)ir i_ n 2V2

к2 1%

 

r U D2V2

 

11

P J_sk

SHk 1JJL1

P V_LSkK SHk

 

 

P2v,

Sk(* SHk-1)2 (1-P 2vi2Sk'<SHk)1

, SHk,ik<p)1 ■ i s

k —

A

------------------------------

 

1

---------------l-)!*4

 

 

 

[ 1 - P 2v 2 sk ( . 2SHk -

Подставляя полученные выражения в уравнения (3 .8 ), окончательно получим:

Т5И(Р) - 2 I k = i

n

XSH(p) 2 X

k - 1

Н.

Vj Sk[l—p2^|HkVj^sk

PKSHkV_[SkHk

-

Г1 _ n2K2

V2 1

1/2

u p SHk vjS k J

 

где параметр p определяется формулой (3 .4a). Уравнения (3 .9 ), яв­ ляющиеся частным случаем уравнений (3 .8 ), соответствуют случаю эллиптической анизотропии и тождественны уравнениям И.И. Гурвича

[2 5 ]

и Сегонзака и Лахеррера [ 8 ], полученным иным путем.

По

(3 .5 ), (3 .6 ) и (3 .8 ) можно рассчитать годографы отраженных

волн для многослойной поперечно-изотропной среды. При этом задача вычисления точек годографа сводится к вычислению функций Т(р) и

83


Х(р) при

заданном шаге изменения параметра р или угла i ^ „Соответ­

ствующая

программа для ЭВМ М -20 составлена А.И. Дозоровым на

кафедре геофизики Геологического факультета МГУ.

2.

ПРИБЛИЖЕННЫЕ ПРЕДСТАВЛЕНИЯ УРАВНЕНИЙ ГОДОГРАФОВ

ОТРАЖЕННЫХ ВОЛН ПРИ НЕБОЛЬШИХ УДАЛЕНИЯХ ОТ ИСТОЧНИКА

Найденные уравнения годографов отраженных волн для рассматрива­ емой модели поперечно-изотропной среды оказываются достаточно слож­ ными, в особенности для Р и SV волн. Для некоторых важных в прак­ тическом отношении случаев эти уравнения можно значительно упро­ стить. Так, существенное упрощение (3 .8 ) можно получить для случая небольших углов падения 9 т.е. для небольших по сравнению с глуби­ ной залегания отражающей границы удалений от источника X Как из­ вестно [9 1 , 1 0 0 ], этот случай наиболее характерен для стандартных модификаций метода отраженных волн.

Рассмотрим вначале приближенные уравнения годографа отраженных Р волн для небольших углов падения в ^ В главе Н было показано, что индикатрисы лучевых скоростей Р волн для углов в < 2 0 -2 5 ° пред­

ставляют собой ЭЛЛИПСЫ С ПОЛУОСЯМИ

Vp р И КрУ^р.

В том случае, если максимальный

из углов падения 0^ в покрыва­

ющей толще не превосходит 2 0 -2 5 °, уравнения для нормальной скоро­ сти в k-м слое можно записать в виде

V р, (i) = V, p , [ l

+ U2 - l)sin2 i ] 2,

(3 .1 0 )

Рк

l _ F k

р р

 

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

п

 

Тр(р) « 2 2 -

 

k= 1

~2 X

 

Рк г - Р 2у|рк *Рк

 

(3 .1 1 )

п

P*Pkj_Pk Нк

Хр(р) * 2 2 -

[ 1 2у2

к= 1

 

1 Рк

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

фов волн SV

Окончательные уравнения для этих волн

имеют

вид

(3 .1 1 ), где вместо ^_р^ и кр^ следует взять Y_SVk

и “SVk

Пара­

метры /<р^ и

Kgyk выражаются через независимые' эффективные упру­

гие параметры поперечно-изотропных слоев формулами

(2 .2 8 )

и

( 2 . 2 0 ) .

 

 

 

84


Уравнения (3 ,1 1 ) существенно проще общих уравнений (3 .8 ), од­ нако они являются параметрическими. Известно, что представление го­ дографов в параметрическом виде неудобно по ряду причин 1 1 0 1 ]. Прежде всего это затрудняет аналитическое исследование их свойств, а следовательно, и выбор наиболее рациональных схем интерпретации. Поэтому желательно найти представление годографов (хотя бы прибли­ женное) в явном виде Т=Т(Х), При небольших по сравнению с глуби­ ной залегания границы удалениях от источника уравнения годографа можно записать в виде ряда Маклорена. Такое представление известно

для изотропных сред [92, 1 0 4 ].

В нашем случае следует воспользо­

ваться точными уравнениями годографа

(3 .8 ) для нахождения

коэффи­

циентов ряда

 

 

 

 

™ » - r w ip -o * 1I! Ь . ./ I

Х . Т !

т р р х р х р р т р

Р --0

■X2 +.

Ф р=0

 

[ХрР

 

 

 

 

При этом для вычисления производных Тр и Хр и производных более высоких порядков при р »0 удобно использовать приближенное представление для нормальной скорости (3 .1 0 ). Выкладки, необходи­ мые для этого, хотя и громоздки, приводят к весьма простому резуль­ тату. Так, для случая Р волн, ограничиваясь членами не выше чет­ вертой степени относительно X, имеем

Тр(Х)

То +

-2

Х "

8 х

R(X),

(3 .1 2 )

 

 

 

 

 

г«е

 

n

2Нк

 

 

 

 

 

Т0 =

2

1 Рк

 

 

 

 

 

 

к=1

 

 

 

 

 

а = (

2

2Р-2

HkV. рк )

;

 

 

 

 

к=1

 

Рк

к ± ™

 

 

 

 

 

п

 

 

 

 

п

 

 

Ь= (

2

2 ^

HkV,3ni)(

2

27(2

H.V ^ }-4

 

 

к -л

 

рк

-Lpk

k=i

Рк к± рк

Разложения вида (3.12.)

справедливы и для поперечных волн SV и

SH, В первом случае вместо

крк и

^ р к в формулах (3 .1 2 ) следует

использовать кажущиеся

коэффициенты анизотропии волн SV, т.е.

«гуь

, и скорости

V с k

,во втором -

ист1шный коэффициент анизотро­

пии

SH

волн

KS,]k

 

 

 

 

Для исследования точности представления годографов по формулам (3 .1 2 ) вычислим следующий член разложения. Поскольку ряд (3 .1 2 ) оказывается знакочередующимся, можно считать, что остаток этого

85


ряда Н(Х) по абсояютной величине не превосходит последнего члена разложения, т.е.

R(X)<| 2

 

 

-7

 

 

' I pD

 

1к=1

 

 

Продолжая это неравенство, можно получить следующую предельную

оценку погрешности представления

годографов (3 .8 )

в виде (3 .1 2 ):

2 Нь

 

 

 

 

(к)

*Pmax \

/j^Pmax^

 

X

R(X)« 0,125 - -

 

 

2кр

n 2 Hl

jLPmin \ KPmin

.LPmin ,

 

 

 

(к)

 

 

 

 

 

 

 

 

( 3 . 1 3 )

Приведем примеры оценки погрешности представления годографов

отраженных продольных волн по формуле (3 .1 3 ) для некоторых воз­

можных значений параметров для осадочных толщ.

Пусть «ртах =1>05,

кп ■

= 0 ,8 , V п

= 4 0 0 0

м/сек,

V,

. = 2 0 0 0

м/сек, 2 Н т= 3000 м ,

Fmin

_L Ртах

 

 

.LPmin

 

к

X < 1 000 м. В этом случае

Н(Х) - 0 ,0 0 3

сек. Если

Vj_ р

.^=2500 м/сек,

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

то

R(X) <; 0 ,0 0 2 4

сек. Это пока­

зывает, что полученные приближенные представления для годографов отраженных волн в поперечно-изотропной среде (3 .1 2 ) дают достаточ­ ную для практики (стандартные модификации МОВ) точность даже для сравнительно сильно дифференцированных по скоростям слоистых " поперечно-изотропных сред.

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

фов продольных волн влияние анизотропии при небольших удалениях от источника определяется параметрами ^ и К, а для волн SV - пара­ метрами у , К и кр

3.ГОДОГРАФЫ ОТРАЖЕННЫХ

ДЛИ НЕКОТОРЫХ ЧАСТНЫХ СЛУЧАЕВ СТРОЕНИЯ

ПОКРЫВАЮЩЕЙ ТОЛЩИ

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

86