Файл: Невский, М. В. Квазианизотропия скоростей сейсмических волн.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 22.10.2024
Просмотров: 82
Скачиваний: 0
Горизонтальная проекция отрезка луча в k-м слое на линию профи ля наблюдения определяется формулой
ЛХк 2Hktg0k
Как было показано, сейсмический луч отраженной волны для рас сматриваемой модели среды симметричен относительно перпендикуляра! проведенного в точке отражения к отражающей границе. Поэтому для получения полного времени пробега по лучу от источника до приемника следует просуммировать времена ЛТk и полученную сумму удвоить. То же следует сделать и для величин АХ ^ В результате получим
п211,
Т- 2 _____________
k=l vk(0k)cos0k '
п
X = 2 2Нк tg 6>к.
к = 1
Последние выражения в принципе описывают зависимость времени пробега отраженной волны от расстояния до источника. Однако для вы числения годографов они неудобны по ряду причин. Во-первых, в попе речно-изотропной среде лучевые скорости не выражаются в явном виД
через эффективные упругие параметры V^pk, крк, *SHk’ >Ц_к и ^ ки угол #к , что создает трудности при расчетах. Во-вторых, в указанны формулы входят значения углов преломления в каждом слое вк, связь1 между которыми, как показано, достаточно сложна. Для получения уравнений годографа в форме, удобной для вычисления, следует прежД
всего |
перейти к параметрическому представлению для vk и #к по фор |
|||
мулам |
(2 .1 1 ). |
|
|
|
Подставляя в последние уравнения для X и Т параметрические выр |
||||
жения для vk |
и #к через ik и замечая, что углы |
ik по |
формулам |
|
(3 .5 ), |
(3 .6 ) |
выражаются через параметр р, т.е. |
iк=iк(р), |
оконча |
тельно получаем параметрические уравнения годографа отраженных вd в виде:
Т(р) = 2 |
2 |
Нк |у^ 1 к (р)Ь v \ [ik(p)]l sec ik(p)+ arclg W |
p) ] , |
|||
|
к = |
|
|
|
v k[ik(p)) |
|
|
n |
|
|
Vk[ik(p)] |
|
|
X<p)' \ |
s |
1 Mk‘e |iltp)+arc‘e vkTi3 P T f |
(3 .8 ). |
|
||
|
|
|||||
где w |
|
dVk , |
_ |
sin i! |
|
|
|
“di- 4 |
’ P - |
V i(ij) |
|
|
82
Последние формулы дают параметрическое представление годографа отраженной волны для случая многослойной поперечно-изотропной по
dVk |
= 0, т.е. покрывающая толща является |
крывающей толщи. Если |
|
di |
|
изотропной, из формул (3 .8 ) |
получаются известные уравнения годогра |
фа отраженной волны для изотропной покрывающей среды.
Формулы (3 .8 ) справедливы для всех типов волн в поперечно-изо тропных средах: Р, SV и SH Различие состоит лишь в том, что для вычисления их годографов следует выбирать соответствующие зависи
мости V]^(i), i ^(р) |
и значение параметра р Для случая поперечных |
|||
волн SH |
уравнения |
(3 .8 ) |
можно существенно упростить, |
используя |
формулы |
(2 ,2 3 ) и (3 .5 ). |
В самом деле, выражение для v |
и ^§Ик |
в функции параметра р по указанным формулам можно после некоторых
преобразований представить |
в виде |
|
|
|||
|
|
|
|
2 |
2 |
2 |
г. . .. |
Pv|SK1 1 + (* S H k -1) [ ] - P 2vJ.SkKSHk |
|||||
tg 0SHk[ *k<P} ] " “ = ----------------------------------------------------- |
|
|
(k2 _ l)ir i_ n 2V2 |
к2 1% |
||
|
r U D2V2 |
|||||
|
■ 11 |
P J_sk |
SHk 1JJL1 |
P V_LSkK SHk |
||
|
|
P2v, |
Sk(* SHk-1)2 (1-P 2vi2Sk'<SHk)1 |
|||
, SHk,ik<p)1 ■ i s |
k — |
A |
------------------------------ |
|
1 |
---------------l-)!*4 |
|
|
|
[ 1 - P 2v 2 sk ( . 2SHk - |
Подставляя полученные выражения в уравнения (3 .8 ), окончательно получим:
Т5И(Р) - 2 I k = i
n
XSH(p) 2 X
k - 1
Н.
Vj Sk[l—p2^|HkVj^sk
PKSHkV_[SkHk
-
Г1 _ n2K2 |
V2 1 |
1/2 |
u p SHk vjS k J |
|
где параметр p определяется формулой (3 .4a). Уравнения (3 .9 ), яв ляющиеся частным случаем уравнений (3 .8 ), соответствуют случаю эллиптической анизотропии и тождественны уравнениям И.И. Гурвича
[2 5 ] |
и Сегонзака и Лахеррера [ 8 ], полученным иным путем. |
По |
(3 .5 ), (3 .6 ) и (3 .8 ) можно рассчитать годографы отраженных |
волн для многослойной поперечно-изотропной среды. При этом задача вычисления точек годографа сводится к вычислению функций Т(р) и
83
Х(р) при |
заданном шаге изменения параметра р или угла i ^ „Соответ |
ствующая |
программа для ЭВМ М -20 составлена А.И. Дозоровым на |
кафедре геофизики Геологического факультета МГУ. |
|
2. |
ПРИБЛИЖЕННЫЕ ПРЕДСТАВЛЕНИЯ УРАВНЕНИЙ ГОДОГРАФОВ |
ОТРАЖЕННЫХ ВОЛН ПРИ НЕБОЛЬШИХ УДАЛЕНИЯХ ОТ ИСТОЧНИКА
Найденные уравнения годографов отраженных волн для рассматрива емой модели поперечно-изотропной среды оказываются достаточно слож ными, в особенности для Р и SV волн. Для некоторых важных в прак тическом отношении случаев эти уравнения можно значительно упро стить. Так, существенное упрощение (3 .8 ) можно получить для случая небольших углов падения 9 т.е. для небольших по сравнению с глуби ной залегания отражающей границы удалений от источника X Как из вестно [9 1 , 1 0 0 ], этот случай наиболее характерен для стандартных модификаций метода отраженных волн.
Рассмотрим вначале приближенные уравнения годографа отраженных Р волн для небольших углов падения в ^ В главе Н было показано, что индикатрисы лучевых скоростей Р волн для углов в < 2 0 -2 5 ° пред
ставляют собой ЭЛЛИПСЫ С ПОЛУОСЯМИ |
Vp р И КрУ^р. |
В том случае, если максимальный |
из углов падения 0^ в покрыва |
ющей толще не превосходит 2 0 -2 5 °, уравнения для нормальной скоро сти в k-м слое можно записать в виде
V р, (i) = V, p , [ l |
+ U2 - l)sin2 i ] 2, |
(3 .1 0 ) |
|
Рк |
l _ F k |
р р |
|
Подставляя это уравнение в формулы (3 .8 ), после необходимых пре образований получим приближенные уравнения годографа отраженных Р волн, справедливые для небольших удалений от источника,
п |
|
Тр(р) « 2 2 - |
|
k= 1 |
~2 X |
|
Рк г - Р 2у|рк *Рк |
|
(3 .1 1 ) |
п |
P*Pkj_Pk Нк |
Хр(р) * 2 2 - |
[ 1 -р 2у2 |
к= 1 |
|
|
1 Рк |
Аналогично можно получить и приближенные уравнения для годогра
фов волн SV |
Окончательные уравнения для этих волн |
имеют |
вид |
(3 .1 1 ), где вместо ^_р^ и кр^ следует взять Y_SVk |
и “SVk |
Пара |
|
метры /<р^ и |
Kgyk выражаются через независимые' эффективные упру |
||
гие параметры поперечно-изотропных слоев формулами |
(2 .2 8 ) |
и |
|
( 2 . 2 0 ) . |
|
|
|
84
Уравнения (3 ,1 1 ) существенно проще общих уравнений (3 .8 ), од нако они являются параметрическими. Известно, что представление го дографов в параметрическом виде неудобно по ряду причин 1 1 0 1 ]. Прежде всего это затрудняет аналитическое исследование их свойств, а следовательно, и выбор наиболее рациональных схем интерпретации. Поэтому желательно найти представление годографов (хотя бы прибли женное) в явном виде Т=Т(Х), При небольших по сравнению с глуби ной залегания границы удалениях от источника уравнения годографа можно записать в виде ряда Маклорена. Такое представление известно
для изотропных сред [92, 1 0 4 ]. |
В нашем случае следует воспользо |
|||
ваться точными уравнениями годографа |
(3 .8 ) для нахождения |
коэффи |
||
циентов ряда |
|
|
|
|
™ » - r w ip -o * 1I! Ь . ./ I |
Х . Т ! |
т р р х р х р р т р |
Р --0 |
■X2 +. |
Ф р=0 |
|
[ХрР |
|
|
|
|
|
При этом для вычисления производных Тр и Хр и производных более высоких порядков при р »0 удобно использовать приближенное представление для нормальной скорости (3 .1 0 ). Выкладки, необходи мые для этого, хотя и громоздки, приводят к весьма простому резуль тату. Так, для случая Р волн, ограничиваясь членами не выше чет вертой степени относительно X, имеем
Тр(Х) |
То + |
-2 |
Х " |
8 х |
R(X), |
(3 .1 2 ) |
|||
|
|
|
|
|
|||||
г«е |
|
n |
2Нк |
|
|
|
|
|
|
Т0 = |
2 |
1 Рк |
|
|
|
|
|
||
|
к=1 |
|
|
|
|
|
|||
а = ( |
2 |
2Р-2 |
HkV. рк ) |
; |
|
|
|||
|
|
к=1 |
|
Рк |
к ± ™ |
|
|
|
|
|
|
п |
|
|
|
|
п |
|
|
Ь= ( |
2 |
2 ^ |
HkV,3ni)( |
2 |
27(2 |
H.V ^ }-4 |
|||
|
|
к -л |
|
рк |
-Lpk |
k=i |
Рк к± рк |
||
Разложения вида (3.12.) |
справедливы и для поперечных волн SV и |
||||||||
SH, В первом случае вместо |
крк и |
^ р к в формулах (3 .1 2 ) следует |
|||||||
использовать кажущиеся |
коэффициенты анизотропии волн SV, т.е. |
||||||||
«гуь |
, и скорости |
V с k |
,во втором - |
ист1шный коэффициент анизотро |
|||||
пии |
SH |
волн |
KS,]k |
|
|
|
|
Для исследования точности представления годографов по формулам (3 .1 2 ) вычислим следующий член разложения. Поскольку ряд (3 .1 2 ) оказывается знакочередующимся, можно считать, что остаток этого
85
ряда Н(Х) по абсояютной величине не превосходит последнего члена разложения, т.е.
R(X)<| 2 |
|
|
-7 |
|
|
' I pD |
|
||
1к=1 |
|
|
||
Продолжая это неравенство, можно получить следующую предельную |
||||
оценку погрешности представления |
годографов (3 .8 ) |
в виде (3 .1 2 ): |
||
2 Нь |
|
|
|
|
(к) |
*Pmax \ |
/j^Pmax^ |
|
X |
R(X)« 0,125 - - |
|
|
2кр |
n 2 Hl |
jLPmin \ KPmin |
.LPmin , |
|||
|
|
|
(к) |
|
|
|
|
|
|
|
|
|
|
( 3 . 1 3 ) |
Приведем примеры оценки погрешности представления годографов |
||||
отраженных продольных волн по формуле (3 .1 3 ) для некоторых воз |
можных значений параметров для осадочных толщ. |
Пусть «ртах =1>05, |
||||||
кп ■ |
= 0 ,8 , V п |
= 4 0 0 0 |
м/сек, |
V, |
. = 2 0 0 0 |
м/сек, 2 Н т= 3000 м , |
|
Fmin |
_L Ртах |
|
|
.LPmin |
|
к |
|
X < 1 000 м. В этом случае |
Н(Х) - 0 ,0 0 3 |
сек. Если |
Vj_ р |
.^=2500 м/сек, |
|||
а остальные параметры неизменны, |
то |
R(X) <; 0 ,0 0 2 4 |
сек. Это пока |
зывает, что полученные приближенные представления для годографов отраженных волн в поперечно-изотропной среде (3 .1 2 ) дают достаточ ную для практики (стандартные модификации МОВ) точность даже для сравнительно сильно дифференцированных по скоростям слоистых " поперечно-изотропных сред.
В заключение отметим, что влияние анизотропии скоростей на годо графы отраженных волн SV и Р при небольших удалениях от источни ка определяется кажущимися коэффициентами анизотропии кр и Это означает, согласно формулам (2 ,2 8 ) и (2 .2 9 ), что для годогра
фов продольных волн влияние анизотропии при небольших удалениях от источника определяется параметрами ^ и К, а для волн SV - пара метрами у , К и кр
3.ГОДОГРАФЫ ОТРАЖЕННЫХ
ДЛИ НЕКОТОРЫХ ЧАСТНЫХ СЛУЧАЕВ СТРОЕНИЯ
ПОКРЫВАЮЩЕЙ ТОЛЩИ
Полученные уравнения годографов отраженных волн для многослой ной поперечно-изотропной среды громоздки и сложны для аналитиче ского исследования основных особенностей годографов. Удобнее про вести такие исследования для более простых моделей сред. С этой целью мы рассмотрим случай однородного поперечно-изотропного слоя на полупространстве и случай двухслойной покрывающей среды, когда
86