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

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

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

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

Добавлен: 22.10.2024

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

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

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

откуда следует, что в данном случае К является средним между зна­ чениями К. в прослоях. Пределы изменения параметра К при усло­

вии, что 0 < <— , оказываются следующими:

0^ К < 1.

Заметим, что определение параметра К для изотропной среды мож­ но обобщить и на случай поперечно—изотропной среды путем введе­ ния двух значений коэффициентов Пуассона и i^ y , первый из

которых характеризует сжатие в плоскости изотропии при растяжении в направлении оси симметрии, второй - при растяжении в той же плос­

кости.

 

обобщение приведено в работе Ф.М. Ляховицкого и авто­

Подобное

ра [8 9 ], где формула для К имеет вид

К =

"XZ

1 -

"XZ * VXY

 

 

I /.XY

Для изотропной среды последняя переходит в формулу (2 .2 2 ). Введенные эффективные упругие параметры длинноволнового экви­

валента тонкослоистой среды позволяют представить уравнение индикат< рис нормальных скоростей ( 2 ,8 - 2 ,9 ) в более удобной форме [9 0 ,9 1 ]'

Для поперечных SH волн

(2 .2 3 )

Vsn(i) = V ^ l + U2$h 1) sin2i]^;

для квазипродольных и квазипоперечных волн SV

 

(j) V lP r, ,

, ( . 2

__ n

„ ;„ 2;

*

(2 .2 4 )

'P.SV

{1 + y'f

+ (к2

-

1)

sin^i ± F (i) ] "* ,

/ 2

P

 

 

 

 

 

 

 

 

 

 

К

F(i) = (K2 + y2 )2 sin22i + [U 2

-

2y2 +1) sin2i -

1 + y^]2

Здесь, как и прежде, знак плюс соответствует квазипродольным, а минус - квазипоперечным волнам. Уравнения (2 .2 3 ) и в особенности (2 .2 4 ) более удобны для аналитического исследования индикатрис^нор­

dV

мальных скоростей волн Р, SV и SH. Рассмотрев производные HJL. di

dv

|

. = 0 для всех трех типов

волн в попе»

можно показать, что —

 

diJ‘ i=0,n-/2

 

 

 

речно-изотропной среде.

Кроме того, для волн SV существует допол­

нительное направление i сл, ,

з

у

=0, Это

направление

где —-------|

oSV

dl

1=1 oSV

 

42


соответствует максимальному значению нормальной скорости волн SV,

где Vgy (iQsy) =Vgymax [56, 9 0 ]. Поэтому

целесообразно ввести .коэф­

фициент анизотропии квазипоперечных волн

SV в виде

^SVmax

“s v ‘ — ч Г '

Этот коэффициент не является независимым параметром поперечно­ изотропной среды, но позволяет охарактеризовать степень ее анизо­ тропности по волнам SV.,

Производная

dVp

-------также может обращаться в нуль при 0 < iQp <п /2

для некоторых типов поперечно-изотропных сред [5 6 ]. Условия относи­ тельно эффективных упругих параметров поперечно-изотропной среды,

dVp

при которых —— |. . - =0, будут найдены ниже, di i=iQp

Индикатриса лучевых скоростей поперечных волн SH после введения эффективных параметров выразится в виде

 

 

2

 

VSH {9) V1S.

1 - (-

* SH -

 

2

 

 

 

 

 

KSH

 

Лучевые скорости

квазипродольных и квазипоперечных волн SV,

как

уже указывалось,

не удается представить в явном виде v =v(6 ).

Пред­

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

 

 

4.

НОРМАЛЬНЫЕ И ЛУЧЕВЫЕ СКОРОСТИ

 

 

 

ПРИ МАЛЫХ УГЛАХ ПАДЕНИЯ

 

 

КАЖУЩИЕСЯ КОЭФФИЦИЕНТЫ АНИЗОТРОПИИ

Зависимость

v = v(0)

для случая квазипродольных и квазипопереч­

ных волн SV в

явном виде может

быть получена, в частности, при ма­

лых углах падения 6 ,

т.е. для направлений, близких к оси симметрии.

При горизонтальной или близкой к ней слоистости в осадочных толщах случай малых углов падения в представляет наибольший интерес для обычных модификаций метода отраженных' волн [9 2 ].

43


Для вывода приближенных уравнений v- %(в) рассмотрим вначале

зависимость нормальных скоростей от угла i при его

малых значениях.

Представим уравнение (2 .2 4 ) в виде

 

 

 

v P,SV(i)

= 1 — 4/2 + /3/2sin2i + q/2

1+— sin2i+ — sin4i

L VIP

j

 

q2

q2

J

q = 1 - y f ; Р =

1; « = к2 + у2;

 

 

 

а = 4а2 -

2/9q — 4q2; b = (/3 + 2q)2 -4 а 2.

 

 

 

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

v P|s V(0

~

Cn(asin2i+bsin4i)n

= 1 - q/2 + /3/2 sin2i ±

2

2q2n-l

V1P

n= 0

где Cn ~ коэффициент разложения (1 + x)

Я

 

в ряд Маклорена.

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

квазипродольных и

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

SV:

Vp(i) = V^p [1

+ /Spsin^i ] ,

(2 .2 6 )

Vgy(i) = V|g [1 + /3gysin2i]^ f

(2 .2 7 )

где K p = (K2 + 1)(K2 + 2yf - Ш - y f r 1,

 

?S V = U 2pU -vf) - K2(K2 + 2yJ)-y{]y-f (l-j^ )"1.

Если в формулах (2 .2 6 )

■V

/V

_rt

и (2 .2 7 ) обозначить /Зр = /ср — l,/8gy=/cgy—1

то они будут аналогичны уравнению (2 .2 3 ) для случая поперечных

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

для лучевых скоростей квазипродольных

и

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

получим для малых углов падения:

 

 

vp(0) и у^р

1

%

(2 .2 8 )

1 — ( ------ г---- ) sin20

 

 

 

 

 

 

Vs v (0) *>Vj_s

i - <*

. —

- Я

(2 .2 9 )

) . Л

 

г.2

 

 

 

SV

44


Таким образом, для малых углов падения индикатрисы лучевых ско­ ростей квазипродольных и квазипоперечных волн представляют собой эллипсы с полуосями соответственно Vjp, «pVjp и V^g, kSv YLS"

Величины ТГр и можно назвать кажущимися коэффициентами

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

 

2 + 1) ( К 2 + 2у 2 1)

*Р~

 

(2 .3 0 )

1

~у\

 

 

(1- у 1) - К 22 + 2у\) - у\

KSV

1 +•

(2 .3 1 )

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

казывает, что для углов в < 2 0

- 2 5 скорости,

вычисленные по точ­

ным (2 .1 1 ), (2 .2 4 )

и приближенным (2 .2 8 ) и (2 .2 9 ) формулам,

различаются не более

чем на 1

-2%. В частности,

на рис.7 приводятся

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

Задачу об аппроксимации индикатрисы лучевых скоростей в некото­ ром интервале углов 0f[o,0m] можно поставить и в другом виде.

Vp(в), м/сек

Рис. 7. Индикатрисы луче­ вых скоростей квазипродоль­ ных волн в поперечно-изо­ тропной среде

1 ,1 —расчитанные по

точным формулам

(2 .2 4 ) и

(2 . 1 1 .);

2 ,2' - ио прибли^

женной

формуле

(2 .2 8 ),

справедливой для небольших

‘ углов

20ь.25°

 

45


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

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

Vp(i) = V|p [1 + /3 sin2i ] 2,

где неизвестным является параметр /3. Для нахождения этого пара­ метра по методу наименьших квадратов следует найти минимальное значение интеграла

I - /

m [V 2p (i)~ V2 (i)]2 di,

o

r

где Vp(i) - точное выражение индикатрисы нормальных скоростей;

Vp (i) - аппроксимирующее уравнение; im - максимальное значение уг»

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

Из последнего выражения обычным способом находим уравнение дл* определения

d lm 9

/tV2 (i) - V 2 (i)]2 di ! = 0,

d/У o

r

F

 

 

откуда

 

 

 

 

*m

t V2 (i)sin2i—sin2i ]di I [

‘ m

,

/3 = 1 /

f

sin^idi ]

о

 

*

о

 

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

как было показано,

с коэффициентом /3 по формуле к .'

= /3 + 1,

к 2р — 1

 

у

(2 .3 2 )

кр = П + — ----- + ( 1 - у 2 ) <Mim)]

,

где </>(im) - функция, определяющая влияние интервала углов, в кото­ ром производится аппроксимация индикатрисы лучевых скоростей

46