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

книги из ГПНТБ / Невский, М. В. Квазианизотропия скоростей сейсмических волн

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

n-1

 

HkVn(i)

у

 

 

2 --------------------------------

 

 

 

 

 

k l

V k t V 2 (i)-V|sin2i 1 2

 

 

V_(i) [1+ (’

d In Vn 2

A'

dlnVn

(3 .22)

jj

) ]

cost i + arctg (

—;-----)]

n

 

 

 

ai

 

n-1

 

H^Vj^sini

d In Vn

X(i)= 2 S

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

 

 

^ H ntg[i + arctg(—

----- )].

k=1

[V2 (i)-V^ sin2i]^

 

 

Параметр i

в этих уравнениях также имеет смысл угла преломле­

ния нормали в анизотропном слое. Если скорость в направлении, пер­ пендикулярном оси симметрии в анизотропном слое, является наиболь-'

шей в покрывающей толще,

т.е. V| |n='<nV| n >Ук(1<к<п-1), значения i

изменяются от О до п/2.

Если максимальная скорость в покрываю­

щей толще наблюдается в одном из изотропных слоев, то наибольшее значение i=imax определяется уравнением

ет in

1

 

^п^тах^

Vшах "

Используя

уравнения (3 .2 1 ), исследуем в аналитическом виде

свойства годографа для данной модели среды. Рассмотрим при этом

случай

Vi <V2(i),

имеющий наибольшее практическое значение. Произ-

 

А

dX

dT

 

водные годографа

---- и — равны

 

 

 

di

di

 

— -

= 2V2(i)A(i), -p— = 2siniA(i),

 

di

*

di

 

 

где

 

 

 

 

A(i)

H iVi[V2(i)c° si -

V2(Osini]3+ H2tV2(i)+V2(i)][V2 (i)-V 2(sin2i]

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

 

 

----------------------- -•

 

[ V2(i)cosi—V2(i) sin i]2 [ V 2(i)—V2 sin2i]3i/2

 

 

 

 

(3 .2 3 )

Знаменатель производных dX и

^ так же, как и в случае одно-

 

 

 

di

di

родной поперечно-изотропной среды, не обращается в нуль и сохраня­

92

знак при it(0, п-/2)- Числитель может изменять знак, если выра-

d2v

,.

 

 

.

.

 

жение ------+V(i)

также изменяет знак при и(0,п-/ 2)

 

di2

 

.

dX

dT

 

 

 

появляются

петли

При изменении знака выражении для — и —

на годографе.

 

 

di

di

 

 

Однако в случае двухслойной среды, в отличие от одно-

 

 

 

d 2V

 

не является

достаточ-

слойной, изменение знака выражения------ + V(i)

 

 

 

di2

 

 

 

ным признаком появления петли на годографе.

Это следует из выраже­

ния (3 .2 3 ).

 

 

 

Когда

d2V

выражение для A(i) сохраняет

sign

-+V(i) =const,

 

di:2

 

 

знак и не обращается в нуль для

if (О.тг/2).

Первая производная го­

дографа при этом

равна

 

 

dT

sin i

 

 

(3 .2 4 )

1 х = "V^ID

 

 

 

иобращается в нуль в точке минимума годографа (i = 0, Х = 0 ). Вторая производная при выполнении условия (2 .1 2 ) определяется

равенством

d2 Т

V2(i) cos (i) — V2(i)sin i

dX2

2V32

(i)A(i)

 

Из результатов предыдущего раздела и (3 .2 3 ) следует, что производ-

d2T

ная------- не обращается в нуль при любых конечных значениях X, а dX2

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

Эти равенства.получаются путем предельного перехода при i-»n/2 в Формуле (3 .2 4 ).

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

93

Взаключение рассмотрим особенности лучевой картины для мэдefl|

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

всех

направлений больше скорости в изотропном слое, т.е. v

v '

---Vj

С помощью формулы (2 .1 5 )

запишем закон преломлений лучей

для рассматриваемого случая:

 

 

V

*

V1

v2

 

 

-

---■-----

 

 

(3.25)

 

 

sin

d In v2

 

 

 

 

) ctg 6>2 ]

 

 

 

 

sin6>2[ 1 - ( —

 

При этом для изотропного слоя

мы использовали равенства v^(0)s

= V^= const,

, поскольку в изотропной среде нет разницы между

направлением луча и направлением нормали к фронту и лучевой и нор мальной скоростями. Из формулы (3 .2 5 ) следует, что в тех точках,

d Inv2

 

где --------

- 0, преломление лучей происходит, как и в изотропных

dв

 

 

средах,

по закону

 

 

vJ d2>

в2>в1-

sin 02 =_£_-----sin(9j,

Обозначим угол в определяемый последним равенством, через

dIn v2

б2 . и рассмотрим случай, когда —^ — >0. Из формулы (3 .2 5 ) полу'

чаем

 

7-1

d In vr

 

sin #2 = sin в2 1 - ( -

-) ctgfl.

>sin02>sin 0j,

~!в~

 

 

V

 

 

или ^2> ^2 > ^1-

 

 

Следовательно, лучи в анизотропной среде преломляются сильнее! чем в изотропной при той же скорости.

 

din v2

 

 

 

Если — -гг—<0,

то величину утла в

можно оценить, используя

 

do

 

^

следующие преобразования:

 

 

 

 

-1

< sin#2,

 

sin #2 =sin09 I]

- ( ^ *ПV2 ctg 02

или

«

2 L

dfl

 

< ^2"

 

 

 

94

Это означает, что при уменьшении скорости с ростом угла падения лучи преломляются слабее в анизотропной среде, чем в изотропной, при прочих равных условиях. Более того, если величина

d In v2

11

de-

clg<V

настолько меньше единицы, что выполняется не­

 

равенство

 

 

 

 

у2(02 )

d In v2

 

 

 

'

<1

 

 

 

d<9

при

v2(0s) > 1 ,

из формулы (3 ,2 5 ) следует, что

 

sin#2 <sin 0-^<sin02 и ^2 <^j.

Иначе

говоря,

при сильном уменьшении скорости v2 (в) с ростом

угла падения может наблюдаться случай, когда угол преломления в анизотропном слое с большей скоростью меньше угла падения. Этот случай можно назвать аномальным преломлением лучей в анизотропной среде.

Результаты аналитического исследования лучевой картины в двух­ слойной среде иллюстрируются расчетным примером. На рис. 24, а приведена индикатриса лучевых скоростей квазипродольных волн для модели типа песчаник - глина, соответствующей третьему типу анизо­ тропии. Параметры модели указаны в табл, 3. Из рисунка следует, что

d Inv2

на индикатрисе имеются области, где—^ — как больше, так и мень­ ше нуля или равна нулю. На рис. 24,6 приведена лучевая картина,

рассчитанная для

двухслойной модели среды при скорости в верхнем

изотропном

слое,

равной 2 0 0 0 м/сек. Скорости в нижнем, поперечно­

изотропном

слое изменяются,

согласно рис. 2 4 ,а

и табл. 3, от

=

=3500 м /сек floV||p = 3 850

м/сек. На рис. 24,6

также приведены

лучи для изотропного второго слоя со скоростью,

равной v2(0)

=

=3500 м/сек. Сравнение лучей, рассчитанных для изотропного и ани­ зотропного слоев, показывает влияние анизотропии на характер пре­ ломления лучей.

Из рис. 2 4 следует, что для анизотропных сред ill типа одновре­ менно могут наблюдаться все аналитически исследованные нами случаи преломления лучей. Реальные значения параметров для тонкослоистой модели поперечно-изотропной среды, выбранные для расчета, дают воз­ можность полагать, что рассмотренные эффекты преломления лучей в анизотропных средах могут наблюдаться в реальных поперечно-изотроп­

ных средах.

 

 

Отметим, что уравнения (3 .1 4 ),

(3 .2 1 ) и

(3 .2 2 ) верны не только

для поперечно-изотропных сред, т.е.

для зависимостей V(i) вида

(2 .2 3 ), (2 .2 4 ), но и для любой функции V(i),

четной и периодической

95

Р и с . 24. Характер преломления сейсмических лучей в поперечно-изо­ тропной среде

а - лучевые скорости

волн

Р

в первом

и втором слоях. Кривая v2(0)

соответствует модели

III табл.

3; б -

преломление сейсмических лу­

чей на границе поперечно-изотропного слоя.

1 - лучи для случая попе­

речно-изотропного слоя, при v2=v2(0);

2

-

лучи для изотропного слоя,

при v2 =v'2 = V2 i= const

 

 

 

 

 

 

с периодом n/N

( N -

целое число), относительно направления перпен­

дикуляра к границам раздела.

Общие уравнения (3 .8 ) предполагают,

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

V(i)

имеют вид (2 .2 3 )

и (2 .2 4 ), поэтому справед­

ливы лишь для поперечно-изотропной среды.

 

 

 

 

4.

 

 

АНАЛИЗ ВЛИЯНИЯ АНИЗОТРОПИ

 

 

НА ФОРМУ ГОДОГРАФОВ ОТРАЖЕННЫХ ВОЛН

Для детального анализа влияния анизотропии на форму годографов отраженных волн обратимся к результатам расчетов годографов для конкретных характерных моделей поперечно-изотропной среды. Расчеты проведены на ЭВМ М -20. Рассмотрим те же два случая: однородный поперечно-изотропный слой на полупространстве и двухслойную покры­ вающую среду с одним изотропным слоем. Первый случай удобен для анализа влияния эффекта анизотропии в чистом виде, т.е. в условиях однородной среды, второй - при наличии эффекта преломления сейсми­ ческих лучей, т.е. в условиях слоистой среды.

Анализ влияния анизотропии на форму годографов отраженных волн удобно провести, сопоставляя годографы для анизотропной среды Т(Х)

с годографами для среды изотропной ТИ(Х), в которой скорости в ело» ях равны скоростям в направлении оси симметрии в анизотропной сре­ де. В силу этого для каждого анализируемого годографа для анизо­ тропной среды приводятся годографы для изотропной среды при vj(0) = = V| j =■ const.

96

Однородная поперечно-изотропная покрывающая среда

Рассмотрим последовательно годографы отраженных волн P,SV и SH для данной модели [102, 103 ].

Квазипродольные волны. На рис, 25 приведены годографы отражен­ ных квазипродольных волн для четырех типов анизотропии. Годографы соответствуют четырем тонкослоистым моделям поперечно-изотропных сред, основные параметры которых приведены в табл. 3, а индикатрисы

лучевых скоростей -

на рис.

19. Годографы приведены в обобщенных

координатах

Т(Х)

.

X

которые позволяют вьивить влияние ани-

т ---------и д = —

— ,

 

т

 

2п

 

 

 

о

 

 

 

Там же

зотропии вне зависимости от абсолютных значений Н и

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

На рисунке видно, что над пунктом взрыва годографы для изотроп­ ной и поперечно-изотропной среды совпадают, а с удалением от источ­ ника все более расходятся. Так, для первого типа анизотропии (мо­ дель I ) заметное различие между годографами для изотропной и ани­ зотропной среды, равное 0 ,0 1 Т0 , наблюдается при X >0,611 и с возраста­ нием X возрастает. Для 111 и IV типов анизотропии годографы для изотропной и анизотропной среды различаются на величину 5Т--0.01 Т0 npjj Х>0,4Н

Характер отклонения годографов для анизотропной среды от годо­ графа для изотропной среды определяется характером зависимости

vp(0)i т.е.

типом анизотропии. Это особенно наглядно проявляется

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

ги (0 -г<£)

 

 

 

'( 0 -

'( 0

 

 

 

 

 

 

 

приведенных на рис. 31

для годографов с

рис. 25. Как следует из

рис. 31, для 1 типа анизотропии величина

<r(f)

всегда положительна

и монотонно возрастает с удалением от источника. Для II типа анизо­

тропии (т=0

для малых удалений от источника, а затем возрастает.

Для типа 111

<г(<f) вначале отрицательное,

а затем становится поло­

жительным и монотонно возрастает с возрастанием £. В случае IV

типа анизотропии <т(£)

- отрицательное для всех

Для больших

Удалений от источника кривые ( т { £ ) стремятся к предельному значению

1 im f т (£ ) = к р — 1

Следовательно, относительное отличие годографов для анизотропной среды от годографа для изотропной пропорционально величине коэффи­ циента анизотропии. Поскольку годографы отраженных волн являются Неограниченно возрастающими функциями, абсолютное различие между Т(Х) И ТИ(Х) для больших удалений от источника может быть очень большим. Так, в качестве примера на рис. 26 приведены годографы отраженных волн для моделей покрывающей толщи I и Ш из табл. 3,

97

1 1257.

тг (Ю

1p(t,)= h p

Р и с . 25. Годографы отраженных волй Р для однородной поперечно-изотроп­ ной среды в обобщенных координа­ тах (т,£)

Номера кривых соответствуют но­ мерам моделей в табл. 3 и обозна­ чают тип анизотропии

Т р , с е н

Р и с . 26. Годографы от­ раженных волн Р для од­ нородной поперечно-изот­ ропной среды при различ­ ных глубинах залегания отражающей границы (спло­ шные линии)

а - для модели I ; б —

для модели III (табл.

3 );

пунктир - Vp(0) = Vi p

=

= const

 

98

рассчитанные в координатах

(Т—X), для двух значений глубин залега­

ния отражающей границы. На рис.

26 видно, что абсолютное отличие

между годографами ТИ(Х) и

Т(Х)

на удалениях X ^2Н достигает око­

ло 0,1 0 сек для модели III

и около 0 ,2 0 сек для модели I и сущест­

венно возрастает при больших значениях X,

Массовые расчеты годографов отраженных волн для различных мо­ делей однородной поперечно-изотропной покрывающей толши, проведен­ ные в нашей работе [9 7 ], показывают, что увеличение значений кр Для I—III типов анизотропии приводит к увеличению кажущейся скоро­ сти и уменьшению времен пробега для дальних ветвей годографов. Уве­

личение значений кр для I

типа анизотропии аналогичным образом

влияет на ближние ветви годографов,

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

Для III,IV типов уменьшает

значение

v* и увеличивает времена пробе­

га отраженных волн для ближних ветвей годографов.

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

значениях кажущихся коэффициентов анизотропии,

равных «р = 1,075 и

«р = 0,7, которые,

согласно расчетам главы

II, можно считать весьма

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

 

кр в реальных средах, аб­

солютное отличие годографов ТИ(Х)

и Т(Х)

 

при

Х= 0,5 Н

составляет

соответственно 0 ,0 0 4 TQ и 0 ,0 3 2 Т0, или 0

,0 0 8

и 0 ,0 6 4

сек при TQ=

=2 сек.

и квазипоперечные SV

волны. Рассмотрим годогра­

Поперечные SH

фы отраженных поперечных волн SV

и SII

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

истых моделей поперечноизотропной покрывающей среды

(табл. 6 ).

Индикатрисы нормальных и лучевых скоростей, соответствующие этим моделям, приведены на рис. 27 и 28 . Как следует из рисунков, для второй модели коэффициент анизотропии к<~у > 1,14 и на индикатрисе лучевых скоростей имеется петля; для первой модели индикатриса пред­ ставлена плавной кривой. Годографы отраженных волн SV и SH для

моделей табл. 6 приведены на рис. 29 .

Там же приведен годограф для

изотропной среды ги (£) при

= ^_S

= const-

Годографы отраженных SH волн рассматривались в целом ряде ра­

бот [2 6 , 4 4 ] при анализе кинематических особенностей сейсмических

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

повторением основных закономерностей.

 

 

Годографы Tgp (X)

и Tgy (X) совпадают над пунктом взрыва, а с

увеличением расстояния X все более расходятся. При этом первый

всегда располагается в области меньших времен пробега. Различие

между Tg|j(X) и ТИ(Х) пропорционально величине

и с

увеличе­

нием Kgpj

увеличивается. Так,

при Х=Н оно достигает 0 ,0 4

TQ для

модели

1( к5Н =*1,16)

и 0,06

Т0 для модели 2 ( k5jj= 1 ,3 8 ) и суще­

ственно возрастает с возрастанием X,

Если на индикатрисе лучевых скоростей нет петли, то годограф от­ раженной волны SV представлен плавной кривой и располагается в

99

Т абл и ц а

б

 

 

 

 

 

 

Параметры моделей тонкослоистых сред для расчета

 

годографов поперечных волн

 

 

 

 

Номер

Условное

v iP ’

V1S.

ПР

nS

m

модели

название

м /сек

м/сек

 

 

модели

 

 

 

 

 

 

 

 

 

 

 

1

Песчаник -

4 5 0 0

2760

1,50

1,67

 

 

известняк

1,0

2

Известняк -

 

2 4 8 0

 

 

 

 

глина

 

4 5 0 0

1,50

2 ,7 6

5,0

Т абл и ц а

6

(окончание)

 

 

 

 

Номер

Условное

 

i s -

 

 

 

модели

название

5

KSH

KSV

 

модели

 

м/сек

 

 

 

1

Песчаник -

 

1 990

1,16

 

 

 

известняк

1 ,1 0

1 ,1 0

0,57

2

Известняк -

 

1665

1,38

1,32

0,41

 

глина

 

1 ,1 0

области времен

Т0 <Тсу (Х)<Ти(Х)(рис.

29, а ).

При этом

различие

между Т

(X)

и Tgy (XI

с удалением от

источника вначале возрастает,

а затем

неограниченно уменьшается, т.е.

годограф Tgy (X)

асимптоти­

чески стремится к годографу Т И(Х).

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

на рис.

29 годограф

для второй модели из табл. 6 состоит из трех

ветвей:

ближней (ab),

средней (Ьс) и дальней (cd), Приближенное

уравнение годографа

(3 ,1 8 ) справедливо для ветви ab, коэффициент

наклона асимптоты по формуле

(3 .156) характеризует ветвь

cd.

Укажем, что в работе [ 61]

теоретические годографы волн

SV не

имеют петель, поскольку при расчете ошибочно использованы нормаль­ ные скорости, а не лучевые.

Следует также отметить, что при полевых сейсмических исследова­ ниях до настоящего времени не получено полного годографа волн SV с ветвью Ьс. Однако дальние и ближние ветви годографа зарегистри­ рованы при экспериментах [2 6 ].

100

V, м / с е к

Рис. 27. Индикатрисы нор­ мальных скоростей SV (сплошные линии) и SH (пунктирные) волн для мо­ делей табл. 6 ; номера кри­ вых соответствуют номе­ рам моделей

Рассмотрим соотношение годографов отраженных SV и SH волн. Вблизи от источника годограф Tgy(X) располагается в области мень­

ших времен по сравнению с годографом TQH(X).

Таким образом, вол-

 

 

 

 

оП

 

ны SV как бы обгоняют волны SH при небольших углах падения. Рас­

хождение ближних ветвей

годографов SV

и SH

волн определяется раз­

личием между к„,,

и

к

и с его возрастанием различия увеличива­

ются.

 

й '

 

 

 

Годографы Tg^X)

располагаются на меньших временах и в области

петли. При этом с

увеличением значения

к^у происходит уменьшение

времен на годографе в области петли, т.е. по ветвям Ьс и cd, а так­ же увеличение области прослеживания ветви Ьс. Для удалений X >2,5Н картина меняется: годографы SH волн располагаются в области мень­ ших времен, а годографы Tgy (X) приближаются к годографам для изотропной среды. Расхождение между годографами Тсу(Х) и Tol,(Х) достигает при Х=ЗН 0 ,0 6 5 TQ для модели 1 и 0 ,1 1 То для модбяи 2 , существенно увеличиваясь при больших значениях X.

Таким образом, анизотропия скоростей в покрывающей толще вызы­

вает расхождение годографов отраженных SV и SH

волн. Оно прояв­

ляется весьма существенно для дальних ветвей, т.е.

при Х > 2 ,5 -3 ,ОН.

При малых удалениях Х<В расхождение годографов

Tgy(X) и TSH(x)

значительно меньше, но и в этом случае оно может быть замечено при эксперименте. Кроме того, при X < Н эти годографы заметно отли­

чаются от

годографа для изотропной среды. Для удалений Н<Х<:ЗН

при наличии петли на годографах волн SV волновая картина может

 

v ( В), м/сек

Рис. 28.

Индикатрисы лу*

чевых скоростей волн SV и SH для моделей табл. 6. Обозначения те же, что и на рис. 27.

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