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

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

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

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

Добавлен: 22.10.2024

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

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

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

где ф - угол между лучом и нормалью к фронту. Согласно последнему уравнению, лучевую скорость можно рассматривать как кажущуюся ско­ рость в направлении сейсмического луча. Эта кинематическая интерпре­ тация лучевой скорости предложена Гассманом [4 J. С другой стороны, лучевая скорость является скоростью распространения упругой энергии, переносимой упругой волной, в анизотропных средах. Это энергетичес­ кая интерпретация лучевой скорости развита в монографии Ф.И. Федо­ рова [2 4 ].

Кажущаяся скорость

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

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

Рассмотрим два бесконечно близких луча, вышедших из источника и

падающих на границу раздела ОХ.

Кажущаяся скорость падающей волны

вдоль границы ОХ по определению равна

v*

НШ

( - ^ - ) =

lim

LAC

 

 

At-0

At

At-*0

At

 

Лучевая скорость определяется

равенством

v -

lim (— - ) =

lim LDC

 

 

At-.О At

At-»0 At

 

Величины

и

Lqq _ длины отрезков АС и DC на рис. 6;осталь

ные обозначения

см. на

рисунке.

Из треугольников

ACD и ABD

находим

Аг = AX sin 6 - AXcos вtgф.

Или после деления на At и перехода к пределам

lim (*—-) cos^r = lim ( — ) sin(0-i/r).

At-vO 1

At-»0 ^

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

v* = v ( в )

cos ф

(2 .1 3 )

 

sin(0- ^ )

38


Рис. 6. Лучевая схема к выво­ ду закона Бендорфа для анизо­ тропной среды

fi —нормаль к фронту; i —угол между нормалью к границе и нор­ малью к фронту (угол падения нормали); в угол падения луча; ф—угол между нормалью к фрон­ ту и лучом

п

г ш

0

А W

 

г

о

X

 

/ в

/

O

'

' Фронт

 

 

 

 

Волны

 

 

 

 

 

w

V

 

 

 

 

/

V

 

 

 

 

Формула (2 .1 3 ) представляет собой обобщение закона Бендорфа для анизотропной среды. Для случая изотропной среды, когда лучи ортого­ нальны фронту и угол ф =0, она переходит в известное уравнение.

Уравнение (2 .1 3 ) можно представить, используя формулы (2 .1 1 ), в двух видах. Во-первых, через нормальные скорости V(i)

v* =

(2 .1 4 )

 

sini

во-вторых, через лучевые скорости v (в)

*

v(fl)

,

, dlnv4

-1

(2 .1 5 )

v

=

1 -

(—JT-) ctg в

 

 

sin0

 

dа

 

 

Представление

v в виде (2 .1 4 )

ранее использовалось Гассманом

[4 ],

однако

для сейсмических задач

в большей степени интересно

уравнение (2 .1 5 ),

поскольку кажущаяся скорость в нем связана с дву­

мя другими важными кинематическими характеристиками: лучевой ско­ ростью v(0) и углом падения луча в. Уравнение (2 .1 5 ) будет ис­ пользовано в дальнейшем для решения прямых и обратных задач сейсмо­ разведки в слоистых поперечно-изотропных средах.

3. ЭФФЕКТИВНЫЕ УПРУГИЕ ПАРА ДЛИННОВОЛНОВОГО ЭКВИВАЛЕНТА ТОНКОСЛОИСТОЙ СРЕДЫ

Упругие константы длинноволнового поперечно-изотропного эквива­ лента тонкослоиотой среды - Сц»Схз, С33, С44 и Cgg —можно выра­

зить через упругие параметры тонких прослоев, их плотности и мощно­ сти [4 9 ]. В частности, для случая периодической тонкослоистой сре­ ды, прослои которой изотропны и однородны, формулы для упругих

39


констант длинноволнового поперечно-изотропного эквивалента можно представить в виде:

 

 

 

’> Л ~1

 

 

-1

n

 

 

 

v t )

 

 

 

^66 = ^ ' Mi7?!»

 

 

 

 

 

 

 

i = 1

C13

 

n

Ajt/j

n

 

V1

 

[ i = 1

 

2

h + 2^j)

 

(2 .1 6 )

 

Xi + 2nJ \ i = l

 

 

 

 

2

n

Mi^i7?! .

 

 

C11 = 2C66

13

 

 

■+ 2

I

Xi + 2Mj

 

 

 

 

 

J33

i =1

 

 

Здесь

Aj,

^

константы Ламе в

i-м слое тонкослоистой пачки;

77i = hj/ i^hi

-

удельная мощность

i-r o

прослоя в тонкослоистой пачке;

п - число таких слоев. В качестве плотности тонкослоистой среды

обычно принимают средневзвешенную величину

 

п

 

 

 

 

 

 

 

 

Р = .2

Pi7?!’

 

 

 

 

 

 

i-l

 

 

 

 

 

тонкослоистой пачки.

где р- - плотность i-ro прослоя

Константы

Си , Сгз, С3344 и

С66, обычно используемые в теории

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

Поэтому целесооб­

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

Введем эффективные упругие параметры длинноволнового эквивален­ та тонкослоистой среды в соответствии со следующими формулами:

С11 X

кр =

_ С33-

с 66

X

 

KSH = С44_

С44

1 С33

Сц/р X

С33/р

с664

С44/р

С44■х

1

_ С3з/р

VIIP

V1P

V||SH

* a:

vis

V||SH

~ViS

(2 .1 7 )

(2 .1 8 )

(2 .1 9 )

40



С

Ylp

.

 

 

( 2 . 20)

 

 

 

 

к -

C13

И

 

( 2 .2 1 )

C33

 

 

 

 

 

 

Константы

Cj:

в формулах (2 .1 7 ) - (2 .2 1 )

определяются форму­

лами (2 .1 6 ).

к р

и Kgp , представляющие собой

отношения скоростей

Параметры

продольных и поперечных волн в направлениях, параллельном и перпен­ дикулярном слоистости, называются коэффициентами анизотропии соот­ ветственно продольных и поперечных SH волн. Эти параметры рассмат­ ривались ранее в работах [ 1 , 2 7 ], где, в частности, показано, что

в то время как *р может быть как больше, так и меньше

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

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

в работе [87].

Используя формулы (2 .1 6 ) и (2 .1 9 ), можно показать,

что параметр у^

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

в про­

слоях, а следовательно, его значения заключены в пределах

0 <у|<

1

<---- - Подробное исследование параметра у в тонкослоистых средах

42

для направлений, отличных от i =0 (направление перпендикулярное слоистости), приведено в [8 8 ].

Скорость длинных упругих волн в тонкослоистой среде в направле­ нии, перпендикулярном слоистости, неоднократно рассматривалась в ря­ де работ, посвященных изучению как квазидисперсии скоростей, -так и анизотропии [1, 8 3 ]. Новым из пяти введенных эффективных упругих параметров тонкослоистой среды является параметр К. В случае од­ нородной изотропной среды он выражается формулой

где и - коэффициент Пуассона.

 

Из формулы (2 .2 2 ) следует, что

равен так называемому ко­

эффициенту распора, применяемому в ряде прикладных задач теории упругости [ 2 0 ].

Для периодической тонкослоистой среды значение К можно пред­ ставить формулой1

пvi

К2

i = 1

1 - V:

 

41