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

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

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

глубинном сейсмическом зондировании при изучении анизотропии в верхней мантии Земли [15, 16, 75, 76].

А.К. Урупов [7 ] предложил использовать наблюдения по методу отраженных волн по системе радиально расположенных профилей для изучения анизотропии в горизонтальной плоскости. Зависимость эф­ фективной скорости от азимута, получаемая в результате интерпрета­ ции, характеризует совместное влияние на эффективную скорость как горизонтальной анизотропности, так и неоднородности среды.

4. ЭКСПЕРИМЕНТАЛЬНЫЕ ДАННЫЕ ПО АНИЗОТРОПИИ СКОРОСТЕЙ СЕЙСМИЧЕСКИХ ВОЛН В РЕАЛЬНЫХ СРЕДАХ

В настоящем разделе приводится сводка экспериментальных данных по анизотропии скоростей сейсмических волн, изученной в условиях естественного залегания горных пород. Сводка составлена по опубли­ кованным результатам полевых сейсмических экспериментов, выпол­ ненных непосредственно для изучения анизотропии скоростей продоль­ ных и поперечных сейсмических волн. Приводимые данные получены с применением методических приемов, описанных в предыдущем раз­ деле.

Следует отметить, что сведений по анизотропии скоростей упругих волн, изученной в сейсмическом диапазоне частот в условиях естест­ венного залегания, сравнительно немного. Тем не менее представляет­ ся интересным и важным проанализировать имеющиеся эксперименталь­ ные результаты. Рассмотрим вначале данные по анизотропии скоростей продольных сейсмических волн.

Анизотропия скоростей продольных сейсмических волн

Результаты изучения анизотропии скоростей продольных сейсмичес­ ких волн приведены в табл. 1 , где кроме значений скоростей и коэф­ фициентов анизотропии даны сведения о районе эксперимента, литоло­ гии и возрасте изученных толщ, а также о методике работ. Необходимо указать, что мы рассматриваем анизотропию скоростей для осадочных толщ земной коры. Имеющиеся в настоящее время данные по анизотро­ пии скоростей в породах верхней мантии нами не приводятся. Эти сведе­ ния содержатся в работах 115, 16, 75, 76].

Используем данные табл. 1 для обсуждения вопросов о масштабах проявления анизотропии в реальных средах, зависимости коэффициентов анизотропии от глубины и литологии пород, а также для обсуждения возможных причин анизотропии скоростей сейсмических волн в оса­ дочных толщах. Ранее подобные вопросы ставились в работе [8 ]; они касались лишь двух скважин в одном из районов Северной Сахары.

В табл. 1 видно, что около половины имеющихся данных о коэффи­ циентах анизотропии получено с помощью способа скважинного сейсмо­ графа. Следовательно, эти коэффициенты анизотропии являются средними и, как указывалось, содержат информацию как об анизотропности, так и об неоднородности исследованных разрезов. На рис. 1 представлена зави-

20

Т а б л и ц а 1

Данные по анизотропии скоростей продольных сейсмических волн в условиях естественного залегания осадочных горных пород

 

Район

Литологическая

Интервал

Средняя или

Коэффициент

 

 

характеристика

глубин, м 1

пластовая ско­

анизотропии,

 

 

пород, возраст

 

рость по верти­

к

 

 

 

 

кали, м/сек

Р

Канада

Сланцы Лорайн

Обнажение

- ■

1,40

США,

Окла­

Известняк Арб-

V

-

1,30

хома

 

рукл

0 -2 1 6 0

3 6 3 0

1,14

США, Техас

Отложения перм­

 

 

ского возраста

 

 

 

 

 

 

0 -2 4 5 0

3 4 2 0

1,13

 

 

 

0 -2 3 7 0

3420

1,19

 

 

Отложения нижне­

0 -2 3 3 0

3640

1 ,1 0

 

 

го мела

1 5 -2 0 0

 

 

США,

Колорадо

Сланцы Пьерре,

1 7 5 0 -2 1 0 0

1,18

1 0 -2 8 0

1 8 0 0 -2 2 0 0

1.14

Методика

-

-

СК с профи­ лем пунктов взрыва

СК с профилем пунктов взрыва на различных глубинах

Источ­ ник

[46]

[70]

Т а б л и ц а 1 (продолжение)

Район

Литологическая

Интервал

 

характеристика

глубин, м

 

пород, возраст

 

США, Техас

Мел Остин

4 -2 8

Средняя или пластовая скорость по верти­ кали, м/сек

2 5 9 0

Коэффициент

Методика

Источ-

анизотропии,

 

ник

кр

 

 

1.17

Сейсмопросве­

 

чивание в не­

[69]

 

 

глубоких сква­

 

 

 

жинах

 

 

Сланец Игл форд

2 0 -4 8

1830

1,33

 

 

Ю. Суматра

Сланцы и песчанис-

0 -4 0 0

1 6 9 5 -1 9 5 5

1,00

СК с профилем

[5 ]

 

тые сланцы, неоген

 

 

 

пунктов взрыва

 

 

 

 

 

 

и МПВ

 

 

 

4 0 0 -9 0 0

2 0 9 5 -2 2 6 0

1,06

 

 

 

 

9 0 0 -1 2 0 0

2 3 5 5 -2 5 5 0

1 ,1 0

 

 

 

 

1 2 0 0 -1 5 0 0

2 6 5 0 -2 8 4 5

1,13

 

 

 

 

1 5 0 0 -2 0 7 5

2 9 3 0 -3 4 1 0

1,04

СК с профилем

 

 

 

 

 

 

пунктов взрыва

 

США, Оклахома

Сланцы

0 -12 0

2 5 9 0

1 ,1 2

СК с профилем

[26]

 

 

 

 

 

пунктов взрыва

 

 

 

 

 

 

и МПВ

 

Сахара, Беррейн

I/есчано-глинистые

0 -8 0 0

2 4 2 5

1,09-1,11

СК с профилем

[ 52]

отложения,

верх­

 

 

 

пунктов взрыва

 

ний мел

 

 

 

 

 

 

Известняки,

ангид­

 

4 0 0 0

1 ,0 9 -1 ,1 1

 

 

риты, верхний мел

 

Глинистые отложе­

8 0 0 -1 2 5 0

2890

1 ,0 9 -1 ,1 4

 

 

 

ния, верхняя юра

 

»

 

 

 

Район не

Глинистые сланцы

4 9 0 -9 0 0

2070

1,075

СК и МПВ-

[ 7 3 ]

указан

 

 

9 0 0 -1 4 4 0

2 6 3 0 %

1,065

 

 

 

 

 

 

 

 

 

 

5 7 0 -1 3 0 0

2 5 4 0

1,08

 

 

Канада, Аль-

Сланцы и пес-

0 -2 7 6 0

3 500

1 ,1 1

СК с профилем

[53]

берта

чаники, мел

 

 

 

 

пунктов взрыва

 

 

 

 

0 -1 4 5 0

2 9 0 0

1 ,1 2

СК и м о в

 

 

 

 

 

1 ,0 6 -1 ,1 0

 

 

 

Эвапориты,

верх­

1 8 0 -1 0 2 0

3460

1,09

СК с профилем

[ 8 ]

 

ний мел

 

 

 

 

пунктов взрыва

 

 

Пески, континен­

1 0 2 0 -2 0 4 0

3 040

1,03

----- и Vяк

 

 

 

 

 

тальные отложения

 

 

 

 

 

 

Известняки,

юра

2 0 4 0 -2 4 5 0

3410

1 ,1 0

 

 

 

Ангидриты,

юра

2 4 5 0 -2 8 1 0

4 8 0 0

1,16

 

 

Т а б л и ц а 1 (окончание)

Район

США, НьюМехико

США, Техас

СССР, Сев. Кавказ

Литологическая

характеристика

Интервал

пород, возраст

глубин, м

Соленосные отло­

2 8 1 0 -3 4 5 0

жения, триас

 

Ангидриты и изве­

1 1 8 0 -1 5 2 0

стняки, нижняя юра

 

Соленосные отложе-

1 5 2 0 -2 0 2 0

ния, триас

 

Сланцы

2 0 2 0 -2 2 9 0

Красноцветные

0 -3 2 8

слоистые пески и

 

сланцы

 

Средняя или пластовая ско­ рость по верти­ кали, м/сек

3 9 4 0

3 9 1 0

3740

Коэффициент анизотропии,

1,02

1,17

1,05

3 3 3 0

1 ,1 2

 

1,077

Чередование пес­

0 -1 4 5 0

2 3 0 0

1,005

ков и сланцев,

 

 

 

палеоген

 

 

 

 

1 4 5 0 -2 0 5 0

3 000

1 ,10

 

2 0 5 0 -3 5 5 0

3 800

1,025

Алевролиты, ар­

■ -

2 4 6 0

1,16

гиллиты

 

 

 

Методика Источ­ ник

СК с профилем

[ 9 ]

пунктов взрыва

 

и УЗК

 

 

( 9 ]

Сейсмопросве- [20] чивание

1,00

1,10

1,Z0

1,30 n r

Р и с . 1. Зависимость среднего коэффициента анизотропии продольных волн кр от глубины

Р и с . 2. Зависимость пластового коэффициента анизотропии продоль­ ных волн кр от глубины для осадочных пород разного литологического состава

1 - карбонатные породы (известняк и ангидриты, мел); 2 - глинис­ тые отложения, сланцы; 3 - песчаники; 4 - соленосные отложения

симость средних коэффициентов анизотропии /Гр от глубины в соот­ ветствии с табл. 1. Как видим, закономерного изменения кр с глу­ биной не наблюдается. Значения достигают 1,20, а в некоторых слу­ чаях 1 ,3 -1 ,4 .

Значительно больший интерес для представления о масштабах про­ явления анизотропии в реальных средах представляют пластовые коэф­

фициенты анизотропии « р .

Они характеризуют

анизотропию скоростей

в однородных в скоростном

отношении пластах

разреза [5 , 8 , 9, 73].

При этом наиболее чистые данные получены Сегонзаком и Стоепом, применявшими ультразвуковой каротаж для изучения скоростного разреза и учета скоростных неоднородностей при интерпретации.

Зависимость пластовых коэффициентов анизотропии от глубины для осадочных пород разного литологического состава приведена на рис. 2

25

Р и с . 3. Средние и пластовые коэффициенты анизотропии продольных волн в зависимости от величины средней или пластовой скорости Vp

1 - средние коэффициенты анизотропии кр ; пластовые коэффициенты анизотропии Хр- 2 - соленосные отложения; 3 - глинистые отложения, сланцы; 4 - карбонатные толши; 5 - песчаники

по данным табл. 1 , большинство из которых получено для глинистых отложений и сланцев. Коэффициенты анизотропии для этих типов пород достигают 1 ,1 0 -1 ,1 3 . Несколько большие их значения (1 ,1 7 - 1 ,2 0 ) возможны в карбонатных породах (известняки, ангидриты). Но в це­ лом экспериментальные данные не свидетельствуют о закономерной связи между коэффициентом анизотропии и литологическим составом пород. Не наблюдается закономерного изменения пластовых коэффици­ ентов анизотропии с глубиной и для пород приближенно одинакового литологического состава.

Итак, экспериментальные данные свидетельствуют о существенном проявлении анизотропии скоростей продольных сейсмических волн в осадочных толщах. Для обсуждения вопрос^ о возможных причинах ани­ зотропии скоростей продольных волн в осадочных толщах рассмотрим зависимость коэффициентов /Гр и кр от величины средней или пласто­ вой скорости. Эта зависимость, полученная по табл. 1, приведена на рис. 3.

Прежде всего заметим, что скорости продольных волн в осадочных горных породах изменяются приблизительно от нескольких сотен мет­

ров

в секунду для слабых, сильнопористых и выветрелых пород до

6 0 0 0

м /сек в крепких разностях известняков, ангидритов, песчаников

[77].

Наибольшие значения коэффициентов анизотропии, согласно рис. 3<

приурочены к средним значениям скоростей

Vp, примерно от 2 0 0 0 до

4 0 0

0

м/сек. В высокоскоростных и низкоскоростных породах замет­

ной

анизотропии не обнаружено. Средними значениями скоростей могу*

обладать и литологически однородные толщи,

однако весьма вероятно,

26

что таким значениям соответствуют гетерогенные толщи, представлен­ ные, например, чередованием прослоев высокоскоростных и низкоско­ ростных пород. Именно в этом случае и наблюдаются достаточно боль­ шие коэффициенты анизотропии, в то время как в высокоскоростных и низкоскоростных, а следовательно, литологически однородных толщах анизотропия отсутствует. Таким образом, на основании анализа экс­ периментальных данных можно предполагать, что наиболее вероятной причиной анизотропии скоростей в осадочных толщах является гетеро­ генность. Гетерогенность проявляется, согласно данным ультразвуково­ го каротажа, в тонкой слоистости - характерном свойстве реальных осадочных сред [78 - 80 ].

Приведенные соображения, естественно, являются лишь предположи­ тельными. Ни в одной из работ по изучению анизотропии скоростей сейсмических волн не ставилась задача сопоставления скоростных раз­ резов по УЗК с соответствующими значениями пластовых коэффициен­ тов анизотропии. Поэтому прямые экспериментальные доказательства связи анизотропии скоростей с тонкой слоистостью в сейсмической практике отсутствуют.

Анизотропия скоростей поперечных сейсмических волн

К настоящему времени имеются лишь единичные опубликованные результаты по изучению анизотропии скоростей поперечных сейсмичес­ ких волн. В табл. 2 приводятся данные, полученные при сейсмораз­ ведочных работах с применением специальных способов возбуждения поперечных колебаний. Как видим, анизотропия скоростей может весь­ ма сильно проявляться для поперечных волн. Коэффициенты анизотро­

пии для

SH волн достигают значений

= 1 ,9 -2 ,1 0 . В эксперименте

Джолли

[26 ] наблюдалась существенная

анизотропия и для поперечных

волн SV.

Данные по анизотропии поперечных волн интересно сопоставить с соответствующими значениями коэффициентов анизотропии продольных

волн. Так, в экспериментах Джолли

[26]

и Уайта и Сенгбуша [69]

одновременно измерены коэффициенты анизотропии Р и SH

волн в

одних и тех же интервалах разреза.

При этом,

по данным Джолли,

для толщи глинистых сланцев

Kgjj =

1,9т 2,0,

а /<р =1,12.

Уайт и

Сенгбуш

показали,

что для

сланца

Игл форд также наблюдается пре­

вышение

кэд

над

к'р Ugpj = 1,60, «р

= 1,33),

а для мела Остин

имеет место

обратное соотношение:

кГи

= 1 ,1 1 ,

кр = 1,17. Судя

по этим данным, встречающееся в сейсмической литературе мнение о том, что коэффициенты анизотропии поперечных SH волн всегда боль­ ше коэффициентов анизотропии продольных волн [56], экспериментом не подтверждается. Оно справедливо лишь для максимальных извест­ ных значений коэффициентов анизотропии поперечных и продольных волн. Так, если максимальные значения коэффициента анизотропии Р

волн составляют

xD = 1,3 т 1,4, то для поперечных они

равны кс1|

1 ,9 *2 ,1 0 .

F

М1=

27

Т а б л и ц а 2

Данные по анизотропии скоростей поперечных сейсмических волн в условиях естественного залегания горных пород ( V, и - скорости вдоль и вкрест слоистости)

Район

Литологичес­

Интервал

 

кая характе­

глубин,м

 

ристика по­

 

 

род, возраст

 

V1SH,

V

 

VI!SV,

Методика

Источ­

IISV,

 

м/сек

м/сек

*SH

м/сек

 

ник

 

 

США, Техас

Сланец Игл

2 0 -4 8

380

-

1,60

380

Сейсмопросвечи­

 

 

форд

 

 

 

 

 

вание на Р и S

 

 

 

 

 

 

 

 

волнах

 

 

Мел Остин

4 -2 8

1070

-

1 , 1 1

1070

 

США,

Окла­

Сланцы

0 -12 0 8 0 0 -9 5 0

 

1 ,9 0 -2 ,1 0

7 8 0 -9 7 0 Трехкомпонент­

хома

 

 

 

 

 

 

 

ный сейсмокаро­

 

 

 

 

 

 

 

 

таж и МПВ

СССР,

Ле­

Глины

 

320

3 2 0

1,56

340

Трехкомпонентный

нинградская

 

 

 

 

 

 

сейсмокаротаж,

обл.

 

 

 

 

 

 

 

сейсмопросвечива­

 

 

 

 

 

 

 

 

ние и МПВ

Болгария,

К ристалличео

До 30

 

2200

1,28

 

Сейсмопросвечи­

Кырджали

кие сланцы

 

 

 

 

 

вание в горных вы­

 

 

 

 

 

 

 

 

работках на Р и S

СССР,

 

Известняк,

6 0 -2 0 0

 

 

 

 

волнах___

 

1000

1000

1 ,3 5 -1 ,4 0

1100

Трехкомпонент-

Херсонская

неоген

 

 

 

 

 

ный сейсмокаро-

обл.

 

 

 

 

 

 

 

таж и МПВ

[69]

[6 ]

[2 0 ]

[97]

Глава в т о р а я

СКОРОСТИ УПРУГИХ волн В ПЕРИОДИЧЕСКИХ ТОНКОСЛОИСТЫХ МОДЕЛЯХ ПОПЕРЕЧНО-ИЗОТРОПНЫХ СРЕД

Настоящая глава посвящена детальному анализу скоростей распро­ странения упругих волн в периодических моделях тонкослоистых сред. Несмотря на то что теоретический подход к описанию кинематики упру­ гих волн в квазианизотропных средах в основном разработан в работах Ю.В. Ризниченко, Постма и других исследователей, применение теории сейсмической квазианизотропии к конкретным сейсморазведочным за­ дачам требует дальнейших теоретических исследований и расчетов.

Нами рассмотрены вопросы о пределах применимости теории сейсми­ ческой квазианизотропии, вопросы аппроксимации индикатрис скоростей в квазианизотропной среде, введены эффективные упругие параметры тонкослоистой среды, более удобные для целей сейсморазведки, а так­ же детально проанализирован характер изменения скоростей продоль­ ных и поперечных волн в тонкослоистых средах. На основе анализа скоростей сейсмических волн в указанных моделях выделены типы ани­ зотропии скоростей в реальных тонкослоистых средах. Рассмотрены также факторы, усложняющие явление квазианизотропии скоростей, а именно: влияние анизотропности тонких прослоев и влияние отклонения структуры среды от строго периодической.

1. ПРЕДЕЛЫ ПРИМЕНИМОСТИ ТЕОРИИ СЕЙСМИЧЕСКОЙ КВАЗИАНИЗОТРОПИИ

Согласно работам Ю.В. Ризниченко, Постма и Бейкуса [1, 3 ,2 7 ], тонкослоистые среды ведут себя в отношении скоростей распростране­ ния длинных упругих волн как среды поперечно-изотропные, или квазианиэотропные, по терминологии Ю.В. Ризниченко. Формулы для скорос­ тей распространения упругих волн в тонкослоистых средах были выве­ дены в указанных работах при предположении, что мощности отдельных чередующихся прослоев h- много меньше преобладающей длины вол­ ны А, т.е. h- « А . Для выяснения вопроса о применимости теории сейсмической квазианизотропии к реальным средам представляется важным получить количественные оценки условия тонкой слоистости. Иными словами, следует прежде всего оценить точность формул для скоростей длинных упругих волн в тонкослоистой среде в зависимости от соотношения мощностей прослоев и длины волны.

С этой целью воспользуемся результатами работы С.М. Рытова [2 ], рассмотревшего задачу о распространении упругих волн в среде,

29

состоящей из двух однородных чередующихся прослоев мощностями

и h2, плотностями р ^ и р2 и упругими постоянными Л2, 'р2*

С.М. Рытовым получены дисперсионные уравнения для скоростей про­ дольных и поперечных волн в направлениях, перпендикулярном и парал­ лельном слоистости, при произвольном соотношениимощностей прослое» и длины волны Л.

Используя эти уравнения, можно оценить точность формул для ско­

ростей длинных упругих

волн в тонкослоистой среде в зависимости от

h.

самым решить поставленную задачу о предела?

соотношения -U- и тем

применимости теории сейсмической квазианизотропии.

Для случая продольных волн, распространяющихся в направлении,

перпендикулярном слоистости,

дисперсионное уравнение имеет вид

 

 

 

 

 

 

,

sin k-^h^ sin k<>h2

 

В этом

уравнении

 

 

(p + —j ------------------------ - •

 

 

 

 

 

 

 

со

 

со

 

со

 

^2Р Рп

 

 

U = Тг

 

^2 =

о

’ а = ~v

’ Р = "\7

 

V1P

 

V2P

 

V1P

 

VlPU

 

 

где VIP

 

^l+2^i "J

 

ГЛ2 + 2^2 1 y2

 

 

 

Pi

J

5 V2P = L

P2

1 - скорости продольных вол»

в прослоях, составляющих среду; Vjj>-

скорость

продольных волн в сре­

де в направлении,

перпендикулярном слоистости.

 

 

 

 

 

 

hi

 

 

 

 

 

Заменим при малых

х и

Л синусы их аргументами, а в разло­

жении косинусов, оставляя члены не выше второй степени, получим

^(hi+h

2'

 

 

 

 

 

 

 

 

1 "

 

.(^р)

*(4 )

 

IP V2PJ

( 2. 1)

UP

 

 

 

 

 

 

 

 

 

 

 

 

Из последнего уравнения после некоторых преобразований выводит­ ся известная формула для скоростей длинных продольных волн в тонкск

слоистой среде в направлении, перпендикулярном слоистости

[1 ]:

hl + h2

 

 

VJ.P

 

( 2. 2)

(pihj + p2h2H

 

 

PlV IP

P0 V2P/

 

30

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