Файл: Левшин А.Л. Поверхностные и каналовые сейсмические волны [монография].pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 21.06.2024

Просмотров: 110

Скачиваний: 2

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

И . П О В Е Р Х Н О С Т Н Ы Е В О Л Н Ы В З Е М Л Е

Г л а в а 5

ПОВЕРХНОСТНЫЕ ВОЛНЫ В КОНТИНЕНТАЛЬНОЙ КОРЕ И ВЕРХНЕЙ МАНТИИ

ЗЕМЛИ

Континентальная кора и подстилающая ее верхняя мантия яв­ ляются одними из главных объектов применения метода поверх­ ностных волн. Наиболее часто для определения скоростного раз­ реза этих частей Земли используют дисперсию основных гармоник волн Рэлея и Лява. Для исследования коры необходимы наблюде­ ния в диапазоне периодов 10—50 сеж; определяется либо групповая скорость на трассе эпицентр — станция, либо фазовая скорость между парой станций, лежащих примерно в одном и том же азиму­ те от эпицентра. При этом землетрясения и станции подбирают так, чтобы изучаемая трасса (или участок трассы) проходила через континентальный регион примерно однотипного строения: для регионального изучения коры иногда применяют методику дробления протяженных трасс на участки с разным строением коры с определением дисперсии групповой скорости на каждом участке по большому числу наблюдений [80, 162].

Для изучения верхней мантии до глубин 300—400 км исполь­ зуют дисперсию основных гармоник в диапазоне периодов 40 — 300 сек, при этом для периодов менее 100 сек применяют те же ме­ тоды, что и при изучении коры (обязательно только наличие длиннопериодной аппаратуры); для наблюдений больших периодов пригодны только записи редких сильных землетрясений ^> 7.5). При этом трассы чаще всего оказываются смешанными, пересе­ кающими как океаны, так и разные континенты; извлечь из на­ блюденных данных сведения о дисперсии на чисто континенталь­ ном пути удается только косвенным путем, при некоторых допу­ щениях о свойствах трассы [20, 86, 123].

Для получения дисперсионных кривых применяется визуаль­ ный метод обработки записи, основанный на измерении координат экстремальных точек сигнала [77, 84, 124]; в некоторых случаях используется метод спектрального анализа зацифрованных запи­ сей на ЭВМ [84, 163]. Решение обратной задачи — определение скоростного разреза по наблюденной дисперсии — осуществля­ ется ручным [12, 35, 84], а в последнее время и машинными ме­ тодами подбора модели, для которой расчетная кривая в среднем

89

наиболее близка к наблюденной. При этом на класс моделей, как правило, накладываются жесткие ограничения: например, модель предполагается слоисто-однородной с фиксированным небольшим числом слоев [5, 120]. Часто неизвестными считаются только одиндва параметра модели (например, мощности слоев), а остальные, в том числе и существенно влияющие на дисперсию, фиксируются; поскольку фиксированные параметры фактически известны с по­ грешностями, точность и разрешающая способность метода по­ верхностных волн не получают объективной оценки.

Использование высших гармоник для изучения строения коры и, в особенности, мантии носит эпизодический характер, в основ- 4 ном в связи с недостаточной надежностью их выделения и первич­ ной обработки существующими методами. В редких случаях опре­ деляются дисперсионные кривые второй и третьей гармоник [6, 117, 118, 149—151]; большое внимание уделяется так называе­ мым каналовым волнам Sa, Lg, Rg [6, 9, 45, 59, 72, 78, 107, 155]. Эти фазы, обладающие четкими вступлениями с устойчивыми груп­ повыми скоростями, используют для определения типа коры на трассе: волны Lg, Rg не проходят даже через очень короткие участки коры океанического типа, сильно ослабевают при про­ хождении через тектонические зоны. По наблюдениям этих фаз оценивают также скорости поперечных волн в коре и мантии, хотя не вполне ясно, к каким участкам среды относятся измере­ ния, поскольку природа волн еще окончательно не определена. Ряд авторов связывает их со слоями пониженной скорости в кон­ тинентальной коре и верхней мантии [89, 107, 115, 127]; другие считают, что такие волны могут возникать в результате интерфе­ ренции высших гармоник и в моделях, лишенных внутренних^ волноводов [78, 124, 149, 151].

В этой главе исследуются возможности дальнейшего развития методики изучения строения континентальной коры и мантии с помощью поверхностных волн. Наш подход заключается в мате­ матическом моделировании (расчете) волновых полей в некоторых характерных моделях Земли и анализе возможностей использова­ ния тех или иных свойств поля для определения параметров среды. Основное внимание уделяется свойствам и перспективам исполь­ зования высших гармоник для изучения верхней мантии, новым машинным методам первичного анализа сейсмограмм, исследова­ нию неоднозначности решения обратной задачи. Наряду с модель­ ными исследованиями приводятся примеры интерпретации реаль­ ных данных.

§ 1. Спектральные характеристики

поверхностных волн

в моделях Земли Гутенберга

и Джеффриса

Модели среды и источника. Моделирование волновых полей поверхностных волн проводилось для двух классических моделей Земли — Гутенберга [33] и Джеффриса [35] (рис. 16 и табл. 6). 90



 

 

 

 

Т а б л и ц а

6

 

 

 

 

 

 

Модель

Модель

Модель

А

Модель

Модел ь

 

 

Джеффриса

Гутенберга

Буллена

Андерсона

автора

г

z, км

а,

ь,

а,

ь,

р, г см*

 

 

 

 

 

 

 

«а

 

«ь

 

 

КМ/ССК

км; сек

км, сек

к.\і/сек

 

 

0

0

6,14

3,55

6,14

3,55

2,74

 

1012,5

450

400

1—

19

6,14

3,55

6,14

3,55

2,74

 

1012,5

' 450

400

- 4 +

19

6,58

3,80

6,58

3,80

3,00

 

1012,5

450

400

2—

38

6,58

3,80

6,58

3,80

3,00

 

1012,5

450

400

2 +

38

7,75

4,37

8,20

4,65

3,32

 

135,0

60

200

3

60

7,82

4,40

8,15

4,58

3,35

 

157,5

70

180

4

80

7,88

4,44

8,06

4,48

3,36

 

202,5

90

160

5

100

7,95

4,46

8,00

4,38

3,39

 

225,0

100

100

6

140

8,07

4,50

7,86

4,35

3,41

 

247,5

110

50

7

170

8,18

4,55

7,92

4,36

3,44

 

258,75

115

60

8

200

8,26

4,60

8,04

4,40

3,47

 

270,0

120

100

9

250

8,42

4,68

8,22

4,50

3,51

 

303,75

135

150

10

300

8,58

4,76

8,55

4,62

3,55

 

337,5

150

200

11

400

8,93

4,94

8,96

4,93

3,65

 

360,0

160

250

12

500

9,66

5,32

9,60

5,32

3,92

 

382,5

170

400

13

600

10,24

5,66

10,12

5,62

4,12

 

472,5

210

600

14

800

11,01

6,13

10,88

6,14

4,48

 

1057,5

470

800

15

1000

11,43

6,36

11,44

6,37

4,68

 

1462,5

650

1000

16

1200

11,71

6,50

11,76

6,50

4,80

 

1901,25

845

1100

17

1400

11,99

6,62

12,08

6,62

4,90

 

2351,25

1045

1200

18

1600

12,26

6,73

12,28

6,75

5,00

 

2801,25

1245

1300

Обе они относительно

хороню

согласуются

с данными

объемных

сейсмических волн; наиболее важные отличия моделей — разница в скоростях упругих волн непосредственно под корой и существо­ вание или отсутствие слоя пониженной скорости в верхней мантии. Хотя в последние годы предложен ряд других моделей, несколько лучше согласующихся с теми или иными выборками наблюденных

данных, сравнение полей для моделей Г и Д представляет

интерес

с точки зрения выявления эффекта слоя пониженной

скорости

на волновую картину. Земная кора в обеих моделях представлена двумя однородными слоями одинаковой мощности. Распределе­ ние плотностей, сравнительно слабо влияющее на особенности поверхностных волн, выбрано согласно Буллену (А) [12] и оди­ наково в обеих моделях (табл. 6). Для выяснения эффекта влияния отдельных частей разреза (коры, верха мантии) на свойства волн

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

моделей

сре­

ды, частично отличающихся от опорных (например,

моделей

Г1

и Г2 на рис. 16).

 

 

91


 

z, кп

 

 

z, ни

 

 

Рис. 16. Распределение скорости поперечных волн Ъ и плотности

р по глубине

в

моделях Гутенберга (Г),

Джеффриса (Д), вспомогательных

моделях

(Г1

и

Г2)

 

 

 

 

 

 

Рис. 17. Распределение Qb

в поглощающей модели Андерсона и др. (пунктир)

и модели автора

(сплошная линия). Штриховка — область возможных частот-

нозависимых Сь

 

 

 

 

 

 

 

Для учета

поглощения

необходимо

задаться распределение!!-

QA

и QB в Земле. В расчетах

использовались две модели распреде­

ления Q — Андерсона

и др.

[103] и модель, отличающаяся от мо­

дели Андерсона большими

значениями

QB непосредственно

под

корой (рис. 17 и табл. 6). Для обеих моделей характерны сравни­

тельно высокие

Q в

земной коре, очень малые Q на глубинах

38—300 км и быстро

растущие Q на глубинах свыше 500 км. Пред­

полагалось, что

QA

и QB частотно независимы; для волн Лява

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

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

92


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

При проведении расчетов был принят ряд ограничений, обу­ словленных главным образом техническими соображениями. Диа­ пазон частот был ограничен максимальными значениями периода Т = 300 сек; волны с таким периодом при любом номере к захва­ тывают всю верхнюю мантию; при больших Т необходимо учиты-

~**Тзать эффект гравитации. Со стороны высоких частот расчетный диапазон был ограничен значением периода Т = 1 ~ 2 сек. Предполагалось, что нас интересуют такие удаления от источни­ ка, на которых более высокочастотные колебания в силу изби­ рательного поглощения сильно затухают.

Ограничения на число учитываемых гармоник необходимо вводить в связи с сильным возрастанием затрат машинного вре­ мени при расчетах полей высоких гармоник. Как правило, мак­ симальное к не превышало 10; во многих расчетах теоретических сейсмограмм, где исследовался интервал времен регистрации, со­ ответствующий групповым скоростям меньше 4,5 кмісек, число учитываемых гармоник ограничивалось тремя-четырьмя.

Максимальные эпицентралыіые расстояния г в наших задачах не превышали 10 ООО км, что вполне достаточно для изучения волн на чисто континентальных трассах; для расчета сейсмограмм

;на больших удалениях при разумных затратах машинного времени "^необходимо сузить расчетный диапазон со стороны высоких

частот.

Большая часть расчетов волн Рэлея выполнена для плоских моделей, все расчеты волн Лява проведены для сферических моде­ лей; мы в дальнейшем всюду используем индексы L и В для волн Лява и Рэлея вне зависимости от того, учитывался или нет эф­ фект сферичности.

Днсііерси і скоростей. Волны Лява. Дисперсионные кривые Vhh {Т) и Chh (Т) для первых восьми гармоник волн Лява в моде­ лях Г и Д приведены на рис. 18, 19. При оценке влияния на харак­ тер волновой картины параметров разреза помимо самих диспер­ сионных кривых удобно использовать графики частных произ­ водных фазовой скорости по этим параметрам. На рис. 20 приве-

Ч дены примеры таких графиков для модели Г, к = 1, 2; в качестве

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

Весь

расчетный диапазон групповых скоростей

CFTL (от 5 до

3 кмісек)

можно разделить на четыре участка (I — IV),

существенно

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

93