Файл: Невский, М. В. Квазианизотропия скоростей сейсмических волн.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) - К 2(К 2 + 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