
книги из ГПНТБ / Невский, М. В. Квазианизотропия скоростей сейсмических волн
.pdf
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.