
книги из ГПНТБ / Невский, М. В. Квазианизотропия скоростей сейсмических волн
.pdfЭффективные скорости, определенные на базе 0 -1 0 0 0 м от источ ника способом квадратичных координат, для всех значений мощности второго слоя оказываются меньше средней скорости по вертикали и весьма близки к предельным значениям vgp. Для базы определения
10 0 0-200 0 м от источника эффективные скорости вначале превышают средние, а затем, с увеличением мощности второго слоя или с увели
чением глубины, становятся меньше средних, |
приближаясь к значениям |
|||||||
vep |
Для удалений 2 0 0 0 -3 0 0 0 , |
3 0 0 0 -4 0 0 0 , |
4 0 0 0 -6 0 0 0 и 6 0 0 0 - |
|||||
8 0 0 0 |
м эффективные скорости больше средних и изменяются весьма |
|||||||
мало с увеличением глубины, |
не более чем на 1 0 0 -1 5 0 м/сек. |
Даже |
||||||
для удалений 6 0 0 0 -8 0 0 0 |
м |
от |
источника при глубине залегания |
от |
||||
ражающей границы |
t |
= 9 0 0 |
м эффективные скорости не достига |
|||||
ют значения V2j_p во втором |
слое и тем |
более предельного макси |
||||||
мального значения |
^ м р " |
|
|
|
|
|
|
|
Следовательно, |
даже при удалениях — * 10 |
эффективные скорости |
||||||
невозможно использовать для определения |
^2|ip БО втором слое. |
Су- |
щественное отличие э ф> определенных на больших удалениях от источ ника, от значений средней и средней лучевой скорости, вычисленной в предположении прямолинейности луча, не позволяет использовать эффективные скорости для точного построения отражающих границ.
Ранее это положение отмечалось для однородной поперечно-изотропной среды; результаты расчетов, приведенные на рис. 39, подтверждают это положение для слоистой поперечно-изотропной среды.
<). в л и я н и е АНИЗОТРОПИИ с к о р о с т е й н а р е з у л ь т а т ы
ПОСТРОЕНИЯ ОТРАЖАЮЩИХ ГРАНИЦ
Для точного построения отражающей границы по годографу отражен ной волны произвольной длины, когда покрывающая толща представле на поперечно-изотропной средой, необходимо знать распределение лу чевой скорости либо в покрывающей толще v=v(6), либо вдоль по про
филю v=v(X).
В первом случае по известному закону v=v (в) можно построить лу чевые диаграммы, с помощью которых отражающая граница строится по способу полей времен или по любому другому, разработанному для случая одномерно-неоднородных сред [9 2 ].
Если же известно распределение лучевой скорости вдоль по профи лю, т.е. v~v(X), то каждая точка годографа ( Т- , X •) должна интерпре тироваться с соответствующими значением лучевой скорости V' =vj (X j). В этом случае можно использовать метод засечек, либо приближенный метод окружностей [9 2 ].
Однако в сейсморазведке построение отражающих границ обычно проводят с постоянной скоростью для всех точек годографа. Как уже отмечалось, эффективные скорости существенно отличаются от лучевых
122

и не могут быть использованы для точного построения отражающих границ. То же можно сказать и о средних скоростях по вертикали, с помощью которых на практике обычно строят отражающие границы.
Для изучения ошибок, вносимых в интерпретацию при использовании средней скорости по вертикали для построения отражающих границ, бы ли проведены следующие расчеты.
Для различных моделей однородной поперечно-изотропной покрываю щей среды и плоской горизонтальной отражающей границы, расположен ной на глубине 2000 м, рассчитывались теоретические годографы от раженных волн для удалений до 8 0 0 0 м от пункта взрыва. По годо графам проводилось построение отражающих границ различными мето дами в предположении, что средняя скорость по вертикали V известна (хотя бы по данным сейсмокаротажа) и равна V| „ Подобные расчеты проведены для годографов отраженных квазипродольных волн.
Рассмотрим результаты построения отражающих границ различными способами без учета анизотропии скоростей в покрывающей толще.
Способ полей времен
На рис. 40, а приведены результаты интерпретации годографа отра женной волны для модели покрывающей толщи I из табл. 3 ( 1 тип анизотропии). Для удалений до 5 0 0 м от пункта взрыва ( X < 0,4 Н) граница, построенная по способу полей времен с постоянной скоростью в покрывающей толще ( V=V^p= 2 5 0 0 м /сек), практически совпадает с истинной. Для больших удалений она располагается существенно вы ше истинной и имеет ^вогнутую форму. Так, для удаления 3 0 0 0 м (Х= =1,5 Н) от пункта взрыва абсолютная ошибка при построении границы составляет 4 0 0 м (20% ), а при Х = 6 8 0 0 граница выходит на поверх ность.
Следовательно, неучет изменения лучевой скорости в покрывающей толще не вызывает существенных ошибок при интерпретации достаточ но коротких годографов и приводит к очень большим ошибкам при ин терпретации дальних ветвей годографа отраженной волны.
В случае III типа анизотропии (рис. 4 0,6 ) построение отражающей границы по способу полей времен с постоянной скоростью в покрыва
ющей толще ( V = V^p= 3 5 0 0 |
м /сек) справедливо только для удаления |
До 300 м от пункта взрыва |
(Х<0,15Н). При ббльших удалениях она |
располагается примерно на 200 м ниже истинной, а при дальнейшем увеличении расстояний - существенно выше истинной. Такая сложная форма отражающей границы, построенной с постоянной скоростью в покрывающей толще, объясняется в этом случае неучетом зависимости лучевой скорости от угла, которая для данного типа анизотропии име ет минимум в области углов падения 4 0 -5 0 °.
Итак, для данного типа анизотропии неучет изменения лучевой ско рости в покрывающей толще приводит к существенным ошибкам даже при сравнительно небольшой длине годографа.
Если в покрывающей толще наблюдается обратная анизотропия (мо дель IV из табл. 3 ), то граница, построенная по способу времен с
123
в
( П В ) О |
8 0 0 |
1 6 0 0 |
2 4 0 0 |
3 2 0 0 |
4 0 0 0 |
4 8 0 0 Х , м |
1В00 |
| |
| |
| |
| |
|--------------- |
1— |
Р и с . 40 . Характер ошибок в построении отражающих границ при использовании средней скорости вместо лу чевой для однородной поперечно-изотропной среды
1 - истинное положение отражающей границы; 2 - граница, построенная по способу полей времен; 3 - по способу окружностей [9 7 ]; а - I тип анизотропии (см. табл. 3 ), б-III тип; в - IV тип; ПВ - пункт взрыва
постоянной скоростью в покрывающей толще, |
имеет выпуклую форму |
|
(см. рис. |
40, в). До удалений около 300 м |
от пункта взрыва (Х « |
<0,15Н) |
она совпадает с истинной, а при больших удалениях истинная |
|
и построенная границы все более расходятся. Так, для Х =2400 м |
||
(Х^1.2Н) |
от пункта взрыва ошибка составляет 4 0 0 м, или 20%. |
|
|
|
Способ окружностей |
Построение отражающих границ с постоянной скоростью, равной средней скорости по вертикали в покрывающей толще, проводилось так-- же и способом окружностей. При небольших удалениях от пункта взры ва границы, построенные по способу полей времен и способу окружно стей, совпадают. При больших удалениях способ окружностей дает по
124
сравнению со способом попей времен большую ошибку в построении границ (см. рис. 4 0 ,а и б). Это связано с тем, что первый основан на предположении о том, что изохрона отражения является окружно стью. Известно, что это предположение достаточно хорошо выполняется лишь для небольших углов падения и в случае изотропных сред. Влия ние анизотропии скоростей и отличие изохроны отражения от окружно сти приводят к весьма большим ошибкам при построении границ по дальним ветвям (Х> Н) годографов отраженных волн.
Исключение составляет лишь случай обратной анизотропии (см.
рис. 40, в), где способ окружностей дает результаты несколько лучшие, чем способ полей времен.
Аналогичные расчеты были проведены и для двухслойной покрываю щей среды, в которой верхний слой изотропен, а нижний представлен поперечно-изотропной средой. Индикатриса лучевых скоростей для по перечно-изотропного слоя приведена на рис. 3 9 ,о, годограф отраженной
волны для этого случая - |
на рис. 41, о. Построение отражающей грани |
цы проведено по участку |
годографа на удалениях от 5 8 0 0 до 9 0 0 0 м |
от источника при общей мощности покрывающей толщи 2 0 0 0 м. Приня тые относительные удаления X «3 —■4,5Н примерно соответствуют уда лениям, на которых регистрируются закритические отраженные волны при работах по КМПВ.
При интерпретации годографа предполагалось, что известно положе ние промежуточной границы раздела, а также пластовые скорости по вертикали в первом и втором слоях. Таким образом, были приняты предположения, обычно выполняющиеся на практике.
Построение отражающей границы выполнено по способу полей вре мен с учетом преломления на промежуточной границе. При этом изо хроны во втором слое построены с постоянной пластовой скоростью, для всех направлений равной VnjI= ^2^=3500 м/сек.
Результаты, приведенные на рис. 41, б, показывают, что неучет даже сравнительно небольшой анизотропии (кр = 1 , 1 0 ) в случае сло истой среды также приводит к очень большим ошибкам. Построенная граница закономерно поднимается с удалением от пункта взрыва и
имеет слабовогнутую форму. На удалении 5 0 0 0 |
м от источника (Х = |
= 2,5Н) ошибка построения достигает 20% (4 0 0 |
м) и существенно |
увеличивается при больших удалениях. Положение отражающей площад ки смещается по сравнению с истинными в сторону больших удалений от источника. Угол наклона построенной границы достигает 15°.
Следовательно, характер ошибок в построении отражающих границ остается в случае слоистой поперечно-изотропной толщи таким же, как й при однородной поперечно-изотропной покрывающей среде.
В -заключение отметим, что замена лучевой скорости средней или пластовой скоростью по вертикали в покрывающей толще приводит к существенным ошибкам при построении границ по закритическим вет вям годографов отраженных квазипродольных волн. В некоторых случа ях, например при обратной анизотропии при малых углах падения ( III и IV типы анизотропии), неучет изменения лучевой скорости в покры-
125
Р и с . 41. Характер ошибок в интерпретации дальней ветви годографа отраженной волны в условиях двухслойной покрывающей среды с нижним поперечно-изотропным слоем
а - теоретический годограф для двухслойной модели (показан учас ток, выбранный для интерпретации); б - результаты построения отражаю" шей границы: 1 —истинное положение; 2 - участок границы, построен ный по способу полей времен
вакяцей толще заметно сказывается и при небольших удалениях от пункта взрыва.
Неучет анизотропии скоростей приводит не только к искажению глу бин залегания отражающих границ, но и к искажению формы границы.
Аналогичные выводы будут справедливы и для годографов закритически отраженных SH волн.
Важно отметить, .что относительная ошибка в построении отражаю щей границы численно не равна коэффициенту анизотропии в покрыва ющей толще, а зависит от способа построения границы и удаления от источника. Причем с возрастанием последнего ошибка существенно возрастает и может в несколько раз превышать значение коэффициен та анизотропии в покрывающей толще.
Вьюод работы [ 44 ] о том, что ошибка в построении отражающей границы численно не превышает значения коэффициента анизотропии,
126
справедлив для небольших удалений от источника и в том случае, ког да для построения границы используется эффективная скорость вместо средней по вертикали. Кроме того, этот вывод получен для эллипти чески анизотропной, а не поперечно-изотропной среды.
7. ОН ИЗУЧЕНИИ АНИЗОТРОПИИ СК
ПРИ СОВМЕСТНОМ ИСПОЛЬЗОВАНИИ ГОДОГРАФОВ ОТРАЖЕННЫХ ВОЛН Р, SV, SH И ДАННЫХ СЕЙСМОКАРОТАЖА
Решение обратной задачи сейсморазведки для случая поперечно
изотропных сред, |
как и для изотропных, можно подразделить на два |
этапа [1 0 0 ]: 1 ) |
определение скоростных характеристик изучаемого |
разреза, 2 ) построение сейсмических границ. Первый сводится к эк спериментальному определению индикатрис скоростей или пяти незави симых упругих параметров поперечно-изотропных сред — V,
и К, полностью определяющих индикатрисы.
Как было указано, при известных индикатрисах скоростей в покры вающей толще задача построения сейсмических границ может быть Точно решена при использовании известных приемов интерпретации для изотропных сред. При построении границ не по индикатрисам, а по средним или эффективным скоростям можно допустит!, существенные ошибки, в особенности при интерпретации дальных ветвей годографов отраженных волн. В силу этого определение скоростных характеристик поперечно-изотропной среды является принципиально важной задачей, один из вариантов решения которой рассматривается в этом разделе.
Как известно [ 6 , 8, 9, 20 ], изучение индикатрис скоростей сей смических волн можно проводить при специальных скважинных иссле дованиях с использованием профиля разноудаленных от устья скважины пунктов взрыва. Нами предлагается способ определения пяти эффектив ных упругих параметров поперечно-изотропной среды по данным стан дартного сейсмокаротажа и годографов отраженных волн, зарегистри рованных на небольших удалениях от пункта взрыва.
Задача определения этих параметров с использованием годографов одного типа волн, например только продольных, является достаточно сложной. Так, при удалениях от источника, меньших глубины залегания отражающей границы, как было уже показано, форма годографа опре деляется только двумя упругими параметрами р и /Гр .Поэтому для определения остальных параметров, в частности, истинного коэффициен та анизотропии кр,требуется использовать годографы, полученные до Удалений, существенно превышающих глубину залегания отражающей границы. Однако в подавляющем большинстве практических случаев при интерпретации удаленных участков годографов отраженных волн требу ется учитывать неоднородность покрывающей толщи, например эффект преломления лучей на промежуточных границах раздела. Учет этого эффекта пока не разработан. Еще большие трудности представляет со бой подобная задача в условиях поперечно-изотропных сред.
127
В силу указанных причин для определения пяти упругих параметров поперечно-изотропной среды следует попытаться совместно использо вать годографы отраженных продольных и поперечных волн SV и SH, полученные при небольших удалениях от источника, а также данные сейсмокаротажа по Р и S волнам.
Один из вариантов решения этой задачи предложен в работе В.В. Локдика [61] . Однако, как отмечалось, уравнения годографов от раженных волн, использованные для решения обратной задачи, справед ливы лишь при слабой анизотропии одновременно по Р, SV и SH вол нам, когда различие между лучевыми и нормальными скоростями мало. Кроме того, для интерпретации годографов отраженных волн Р, SV и SH по способу В.В. Локцика [6 1 ] требуется дважды численно диффе ренцировать времена по экспериментальным годографам, что, как из вестно [8 1 ], приводит к большим вычислительным погрешностям.
Для устранения этих недостатков предлагается иной способ опреде ления пяти эффективных упругих параметров поперечно-изотропной сре дь!, принципиально применимый как при слабой, так и при сильной ани зотропии. Рассмотрим вначале случай однородной поперечно-изотропной среды с осью симметрии, перпендикулярной плоской горизонтальной от ражающей границе.
Пусть в результате трехкомпонентного сейсмического каротажа оп ределены пластовые скорости Р и S волн в вертикальном направлении (в направлении оси симметрии среды) в пределах покрывающей толщи. При этом можно найти два упругих параметра поперечно-изотропной среды:
у плР= V±P> |
|
|
v njS |
|
(3 .4 1 ) |
|
У± |
|
уплР |
|
|||
|
|
|
||||
|
|
|
|
|
|
|
По годографам отраженных волн |
Р, SV и SH, полученным вблизи |
|||||
от источника, определяются эффективные скорости. |
Из формулы (3 .2 8 ) |
|||||
находим значение |
в покрывающей толще |
|
||||
|
узф8Н ^эфБН |
|
|
(3 .4 2 ) |
||
|
\lS |
уплБ |
|
|
|
|
|
|
|
|
|
||
Зная предельные значения эффективных скоростей Р и SV волн, нахо |
||||||
дим кажущиеся коэффициенты анизотропии волн Р и |
SV согласно фор |
|||||
мулам |
(3 .2 9 ), |
(3 .3 0 ) |
|
|
|
|
к р |
veP |
УэфР , |
_ |
_ veSV |
уэфБУ |
|
= |
|
|
|
|
|
|
|
Ylp |
уплР |
^ |
YlS |
уплБ |
|
128
По значению кр при известном |
у |
по формуле (2 .3 0 ) находим па |
||
раметр К |
|
|
-L |
|
где |
|
|
|
(3 .4 3 ) |
А = |
_ i)(i_ y2 > |
|
|
|
При известных значениях К, у |
и |
согласно формуле (2 .3 1 ), |
||
определяем величину Кр |
|
|
||
кр |
В+ С 1 Л1 |
|
|
(3 .4 4 ) |
|
|
|
||
1- - |
|
|
|
|
|
|
|
|
|
где |
1 . |
|
|
|
|
|
|
|
|
b =(k2sv - m - y h y Z , |
|
|
||
С = К2 (К2 4 2у2 )+ V2 . |
|
|
||
|
1 |
1 |
|
|
Таким |
образом, |
формулы (3 .4 1 ) - (3 .4 4 ) дают принципиальную воз |
||
можность |
определения всех пяти независимых упругих параметров од |
нородной поперечно-изотропной покрывающей среды по данным сей смокаротажа-и метода отраженных волн, с использованием волн Р, SV и SH.
Найденные формулы обобщаются и на случай горизонтально-слои стой поперечно-изотропной среды. Пусть в такой многослойной среде с вертикальной осью симметрии в каждом слое по данным трехкомпо— нентного сейсмокаротажа известны пластовые скорости Р и S волн в вертикальном направлении. Эти данные дают возможность определить
параметры |
у^^ и Vp^ в каждом слое по формулам, аналогичным |
|
(3 .4 1 ). Пусть, кроме того, имеется серия годографов отраженных |
||
волн Р, SV |
и SH |
от кровли каждого поперечно-изотропного пласта. |
Для определения |
истинных коэффициентов анизотропии SH волн, а так |
же кажущихся коэффициентов анизотропии кр и Kgу следует восполь зоваться следующим уравнением, полученным из формулы (3 .3 8 ):
|
v , |
Т , |
- V 2 , Т |
кк |
ек+1 |
ок+1 |
ек ок |
|
|
(3 .4 5 ) |
|
vlk |
^ок+1 ^ок |
||
|
|||
где T0k+i и Ток - |
времена над источником по годографам, соответ |
ствующим подошве и кровле пласта, a vek+1 и vgk - предельные эф фективные скорости, определенные по этим годографам.
129
9 1257-
Формула |
(3 .4 5 ) справедлива для всех |
трех типов волн - Р, SV и |
||
S1I, Однако при использовании ее для SH |
волн получаем истинный ко |
|||
эффициент анизотропии « g p , а для волн |
Р и |
SV — кажущиеся коэф |
||
фициенты |
к р |
И Kgy в каждом толстом пласте. |
При известных значе |
|
ниях крк |
и |
ксуь определение параметров |
Кк и кр> проводится по |
|
формулам |
(3 .4 3 ) и (3 .4 4 ). |
|
|
Следовательно, предложенный способ принципиально дает возмож ность определить все пять упругих параметров каждого поперечно изотропного слоя в горизонтально-слоистой среде. Формула (3 .4 5 ) позволяет учесть и эффект преломления на границах слоев при неболь ших углах падения лучей. Отметим, что располагая только данными
по Р и SV волнам, |
можно по формулам |
(3 .4 1 ), (3 .4 3 ), (3 .4 4 ) най |
ти четыре параметра (\ р , у , кр и К), |
характеризующих анизотропию |
|
среды по волнам Р |
и SV. ^ |
|
При известных значениях упругих параметров поперечно—изотропной среды можно рассчитать индикатрисы нормальных и лучевых скоростей, которые удобно использовать для расчета лучевых диаграмм и даль нейшего построения сейсмических границ.
Гл ава ч е т в е р т а я
ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ АНИЗОТРОПИИ СКОРОСТЕЙ В РЕАЛЬНЫХ СРЕДАХ
Для экспериментального изучения анизотропии скоростей в осадоч ных толщах проведены комплексные сейсмоакустические исследования в составе экспедиции Института физики Земли (ИФЗ) АН СССР. Для этого был выбран один из районов северного борта Припятского про гиба в Белоруссии. В 1 9 6 7 -1 9 7 0 гг. Е.В. Карусом, Г.Н. Гогоненко— вым, Н.Р. Сивковым, В.В. Бандовым и другими исследователями в районах северного борта Припятского прогиба детально изучено ско ростное строение с помощью ультразвукового каротажа. Результаты этих исследований показывают, что скоростные разрезы по ряду рай онов имеют тонкослоистый характер, сопровождающийся резкой диффе
ренцированностью по скоростям [79, |
107, 1 0 8 ]. В качестве примера |
на рис. 42, а приведен такой разрез |
одной из скважин в районе север |
ного борта Припятского прогиба. Видно, что в широком интервале глу
бин (8 9 0 -1 3 0 0 |
м) он представляет собой ритмичное чередование |
|
тонких прослоев |
с |
повышенными и пониженными скоростями (от 4 7 0 0 |
до 2 5 0 0 м /сек ). |
Колебания скорости в тонких прослоях отмечаются |
на фоне сравнительно слабого изменения с глубиной среднепластовой скорости. Согласно теоретическим представлениям [1, 1 0 9 ], подобный
характер разреза должен приводить к существенной анизотропии ско ростей сейсмических волн.
Тонкая слоистость в сочетании с резкой дифференциацией среды по скоростям, а также сравнительно небольшое изменение среднепласто— вых скоростей с глубиной отличают объекты в пределах северного борта Припятского прогиба от ряда других районов, детально изучен ных с помощью УЗК. Так, в районах Краснодарского края, согласно [ 1 1 0 ], осадочные отложения также характеризуются тонкой слоисто
стью, однако дифференциация по скоростям в тонких прослоях невелика. Скорости в тонких прослоях в основном изменяются на первые сотни
метров в секунду (см. |
рис. |
4 2 ,6 ). |
|
|
|
|
|
||
В |
районах |
Саратовского и Куйбышевского Поволжья, |
в |
Баш |
|||||
кирской АССР |
скоростные |
разрезы, |
как |
показано |
на рис. |
42, в, |
|||
отличаются сложным |
строением [78, |
1 1 1 ]. На разрезах |
выделя |
||||||
ются |
толстые |
пласты с повышенными и пониженными значения |
|||||||
ми скоростей, |
сильно |
различающихся |
по |
величине. |
Тонкая |
слои |
стость в пределах пластов выражена в разной степени. В кар бонатных отложениях скорости в тонких прослоях изменяются на
1000 м/сек и более. В терригенных породах различия |
по скоро |
стям в тонких прослоях существенно меньше (см.рис. 4 2 ,в). |
Приведенные |
131