Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Невский, М. В. Квазианизотропия скоростей сейсмических волн

.pdf
Скачиваний:
12
Добавлен:
19.10.2023
Размер:
7.28 Mб
Скачать

Эффективные скорости, определенные на базе 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

Соседние файлы в папке книги из ГПНТБ