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

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

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

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

Добавлен: 22.10.2024

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

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

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

hl

h2

•г после некото-

-заменяя тангенсы их аргументами при малых —— и —

л

Л

 

рых преобразований получим известную формулу [48]

PlV IS hl + Р2 V2S h2

a

IISH

( 2 . 6)

Pih1 + P2^2

 

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

(2 .7 )

r(Vl|SH} 24г(’ - ч)’ Щ :

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

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

грешностей по формулам

(2 .4 )

и (2 .7 )

в функции h /Л при различных

 

 

r\§

 

V jp

 

 

4 . Как следует

из

значениях п^= у----

и nD ь г-,— приведены на рис.

 

 

2S

F

V2P

 

 

 

 

 

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

при постоянном

значении

h/Л

погрешности возрастают,

особенно резко для скорости

длинных

волн

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

Однако

по­

грешности формул (2 .2 )

и (2 .6 )

не превышают и 1 % при 0,5

 

4 2 ,0 , если

Ь/Л<:

0 ,1 .

Если приближенно принять длину продольной

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

равной 6 0 - 7 0

м

(например, при V =

3 5 0 0 J-3 0 0 0

м/сек

и f = 50

гд ), то формулы

для скоростей длинных волн в периодической тонкослоистой среде бу­ дут выполняться с точностью до 1% для прослоев мощностью до 7 м; если скорости в прослоях различаются не более чем в 2 раза.

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

33

3 1257


e(v), %

Рис. 4 . Зависимость предельных относительных погрешностей от отношения п/Л

1 - по формуле (2 . 7); 2 ~ по формуле (2 .4 ). Параметр кривых

n = V1/V 2

наши результаты с результатами работ И.С. Берзон и В.И. Пасечника [8 2 , 8 3 ], где дано строгое решение задачи о скоростях распростра­ нения гармонической упругой волны в направлении, перпендикулярном

слоистости, в произвольной тонкослоистой среде.

В работе [8 3 ], в

частности, показано, что формулы для скоростей

длинных упругих волн

в тонкослоистой среде, полученные в работах Ю.В. Ризниченко

[1] и

Бейкуса [3], выполняются с высокой точностью,

если максимальная

мощность прослоя в 8 и более раз меньше длины волны, при

общей

мощности тонкослоистой пачки, превышающей длину волны. Таким об­ разом, данные работы [8 3 ] и результаты оценок, приведенных для случаев распространения волн вдоль и поперек слоистости в данном па-- раграфе, хорошо совпадают.

Полученные результаты показывают, что теория сейсмической квази^ анизотропии вполне приложима к реальным сейсмическим средам, по­ скольку, согласно многочисленным данным УЗК [7 8 - 8 0 ], мощные тол-' щи осадочных отложений в ряде районов представляют собой чередова­ ние тонких прослоев мощностью от 1 до 10 м и тоньше при сильной дифференциации скоростей в тонких прослоях.

Существенную роль в формировании квазианиэотропных свойств ре­ альных сред играют весьма тонкие прослои мощностью hj , $ 0,1 Л. Заметим, что на динамические особенности сейсмических отраженных волн в тонкослоистых средах, как показано в работе Г.Н. Гогоненкова и Ю.Г. Антипина [8 4 ], тонкие прослои даже несколько большей мощ­ ности (hj < 0 ,25 Л) практически не влияют. Поэтому можно предпо­ лагать, что для изучения квазианиэотропных свойств реальных сред требуются более детальные измерения скоростей по методике УЗК, чем для изучения динамических особенностей волновых полей в реальных тонкослоистых средах [8 5 ].


2. СКОРОСТИ УПРУГИХ ВОЛН В ПОПЕРЕЧНО-ИЗОТРОПНОЙ СРЕДЕ

Как было найдено, теория сейсмической квазианизотропии примени­ ма к реальным тонкослоистым средам. Результаты работ [1 ,3 ,2 7 ] показывают, что всякая тонкослоистая среда имеет свой длинноволновой эквивалент - поперечно-изотропную среду. Известно, что кинематические особенности упругих волн в анизотропных средах могут существенно от­ личаться от таковых в средах изотропных [2 4 , 8 6 ]. Поэтому для це­ лей сейсморазведки важно рассмотреть основные особенности кинема­ тики упругих волн в поперечно-изотропных средах и прежде всего вопрос о скоростях этих волн.

Согласно общей теории распространения упругих волн в идеально упругих анизотропных средах [ 24 ], скорости упругих волн в них зави­ сят от направления распространения волны. Кроме того, известно, что сейсмические лучи в анизотропных средах не ортогональны фронту вол­ ны, т.е. в общем случае существует различие между направлением нормали к фронту волны в данной точке и направлением сейсмического луча. Следовательно, скорости распространения упругой волны по нор­ мали к фронту и в направлении сейсмического луча будут в общем слу­ чае различны.

Скорость упругой волны по нормали к фронту мы будем называть нормальной и обозначать через V, а скорость в направлении сейсми­ ческого луча — лучевой, обозначив ее через v.

Нормальные скорости

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

Упругие свойства поперечно-изотропной среды в общем случае опи­ сываются пятью независимыми упругими параметрами С ц , С^3, С33,

С44

и

Cgg. Для

изотропной

среды, как известно, имеем два парамет­

ра;

Л

= С13, fi =

С4466 и

Л+ 2р = С33 = Сп • Скорости плоских волн

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

иугол i между нормалью к фронту и осью симметрии [24 ]: для поперечных SH волн

VSH^) = t ~ ^ 4 4 cos^ + С66 sin2i) 3^ •

(2 .8 )

35


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

д о =

- — I а + d ±

I а

[ ( a —d)^ + 4 b ^ ] ^ i

'P.SV

 

где р - плотность среды;

a^C ^siiA +:C44Cos2i; b = (С^з+ С44) cos

d = C^^sin^i + CggCOS^L

 

(2 .9 )

sin i;

Знак плюс в формуле (2 .9 ) соответствует квазипродольным, ми­ нус - квазипоперечным волнам SV.

Функции V(i) являются непрерывными дифференцируемыми четными и периодическими с периодом, равным 180 . Последние два свойства позволяют исследовать нормальные скорости упругих волн в поперечно­ изотропной среде на отрезке углов if [0,17/ 2].

Скорости упругих волн по направлению оси симметрии среды ( i = О) и перпендикулярно к ней ( i = тг/2) оказываются равными:

Vp(n-/2) Vsv ( п/ 2) = ’ VSH (7Г/2) = . ( 2. 10)

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

тям напластования тонких

слоев.

В дальнейшем мы будем обозначать

скорости в направлении этой оси

V j_, в направлении, перпендикулярной

оси симметрии - У ц ,

следуя известной работе

Ю.В. Ризниченко [1].

Следует отметить,

что

нормальные скорости

в анизотропных средах

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

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

Зависимость лучевой скорости v от угла в между лучом и осью сим­ метрии поперечно-изотропной среды можно получить из уравнения оги­ бающей семейства фронтов плоских волн, распространяющихся под все­ ми возможными углами к оси симметрии среды [8 6 ]. Для всех трех типов волн - Р, SV и SHj распространяющихся в поперечно-изотроп­ ной среде, лучевую скорость можно выразить в параметрическом виде через соответствующие нормальные скорости V(i):

v(i)

( 2 . 11)

d In V

0(i) = i + arctg

di

36


z

В

Рис. 5. Геометрическая интер­ претация нормальной V и лучевой скорости v в анизотропной среде

РВ —нормальная скорость; РА —лучевая скорость; ф угол

между нормалью к фронту и лучом

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

Исключить i из уравнений (2 .1 1 ) удается только для волн SH. В этом случае индикатриса лучевых скоростей представляет собой эллипс

С ПОЛУОСЯМИ \^gj^ H ^ | | S H '

Аргумент радиуса вектора индикатрисы лучевых скоростей выражает­ ся, согласно уравнениям ( 2 . 1 1 ), формулой

d In V

в = i + arctg (

 

_di

 

 

 

Ав

 

 

Очевидно, если производная

не меняет знака на отрезке

Ю,

тт/2], то

на кривой лучевой скорости

v(0 ) не будет особенностей ти-

па точек

 

ав

нахо-

возврата, петель и т.д. Рассмотрев производную — t

 

 

dL

 

Дим условия отсутствия особых точек на кривой лучевой скорости в виде

sign

= const

( 2 . 12)

при if [О, п/2]. Условие (2 .1 2 ) означает сохранение знака

выражения

di2

Связь между нормальной и лучевой скоростью нетрудно найти из уравнений (2 .1 1 ) и рис. 5. Соответствующая формула имеет вид

V(i)

v (0 )

cos ф ■’

37