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

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

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

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

Добавлен: 21.06.2024

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

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

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

Возможности использования теоретических сейсмограмм. Рас ­ смотренные примеры демонстрируют возможности применения те­ оретических сейсмограмм для объяснения природы и изучепия свойств экспериментально наблюдаемых фаз поверхностных волн (Sa, Lg, Rg, R, G). Однако этим далеко не исчерпываются воз­ можности, открываемые численным моделированием полей по­ верхностных волн. Одним из наиболее перспективных направлений в использовании теоретических сейсмограмм является анализ магнитуд. Как известно, магнитуды являются условной характери­ стикой интенсивности землетрясений и определяются непосредст­

венно по записям

поверхностных (М) или объемных (т) волн

(см., например, [19,

83, 146]).

 

Имеющиеся калибровочные кривые магнитуд М,

приводящие

измеренные по записям амплитуды к фиксированному

расстоянию,

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

Из рассмотрения сейсмограмм на рис. 33—36 видно, что ослаб­ ление максимальной амплитуды (а следовательно, и магнитуды) по мере погружения очага под кору происходит сходным образом в моделях Г и Д. Это опровергает объяснения дефекта магнитуды при погружении очага в мантию влиянием астеносферы [92]; он связан с тем, что глубина проникания колебаний основной гар­ моники волн Лява и Рэлея с периодами меньше 30 сек, используе­ мых для измерения М, не превышает мощности коры.

Теоретические сейсмограммы могут оказаться весьма полез­ ными при физической интерпретации и развитии критерия обна­ ружения ядерных взрывов по отношению магнитуд тІМ [67, 81, 146]. Они оказываются полезными и при объяснении ряда других наблюдаемых эффектов. Так из сейсмограмм на рис. 33— 36 с очевидностью следует, что относительный вклад высших гармоник быстро возрастает по мере углубления очага: этот эффект хорошо известен экспериментально [6, 109]. Сравнение сейсмограмм для очагов в верхней и нижней частях коры свиде­ тельствует о быстром ослаблении короткопериодных частей сиг­ нала с погружением источника. Этот факт находит хорошее эк­ спериментальное подтверждение в наблюдениях землетрясений с надежно определенными глубинами очага и может быть исполь­ зован для оценок глубин очага в пределах земной коры [76, 144]. Наконец, в следующем параграфе мы продемонстрируем исполь­ зование теоретических сейсмограмм для оценки эффективности и разрешающей способности новых машинных методов обработки сейсмограмм поверхностных волн.

12Î


§ 3. Спектрально-временной анализ поверхностных волн

Приведенные теоретические сейсмограммы могут служить модельными примерами исходных данных, используемых в сей­ смологии при решении обратной задачи — определении строения среды по наблюденному волновому полю. Вид сейсмограмм в зна­ чительной степени определяется параметрами источника, варьи­ рующими от землетрясения к землетрясению и известными ин­ терпретатору с малой точностью и надежностью (механизм, спектр, глубина очага). Поэтому наилучший путь интерпретации — в^ разделении влияния среды и источника, например путем опреде­ ления по имеющимся сейсмограммам инвариантных относитель­ но источника характеристик поверхностных волн. Согласно гл. 13, к ним относятся дисперсионные кривые vkQ (Т) и CkQ (Т), функции %h (Г) и QkQ (Т) отдельных гармоник поверхностных волн. Таким образом, возникает задача разделения гармоник и определения их индивидуальных спектральных характеристик.

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

Переход

от

временной записи uq

(t)

к ее спектру

Uq (со),

как следует

из

формул (1.29), (1.70)

и

проведенных

расчетов,

не дает большого эффекта, поскольку из-за интерференции гар­ моник амплитудный спектр записи оказывается изрезанным, а

фазовый спектр — неустойчивым

в

областях малых амплитуд.

При анализе реальных записей

это

усугубляется добавлением

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

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

122


 

известна

длительное время,

в

сейсмологии

этот

подход в

цифро­

 

вой реализации стал применяться для изучения

диспергирую­

 

щих цугов только в последние

годы,

а в отечественной литерату­

 

ре

вообще

не

рассматривался.

В

связи с

этим

ниже

будет

при­

 

ведено краткое рассмотрение принципа СВАН и

реализующего

 

его

машинного

алгоритма,

а

затем

показаны результаты

обра­

 

ботки этим методом расчетных и реальных сейсмограмм

поверх­

 

ностных

волн.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Принципы спектрально-временного

анализа.

Обозначим

под­

 

лежащую анализу запись поверхностной волны в интервале

 

времени

tH

^

t

 

tk

и (t),

а

ее

преобразование

Фурье

U (х):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U(x)

 

\ e-ixtu(t)dt.

 

 

 

 

 

 

(5.1)

 

Основной

принцип

СВАН

заключается

в

замене

одномерных

 

представлений

анализируемого

сигнала

и (t)

или

U (х)

дву­

 

мерным

представлением У

(t,

со), характеризующим

спектр

сиг­

 

нала

в

окрестности

частоты

со в

моменты времени, близкие

к

L

 

У (t,

со)

определяется

так:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У (t,

со) =

±-\u(x)H

 

 

( ^ )

 

dx,

 

 

 

(5.2)

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

<

где

H

(х) — некоторая

функция,

имеющая

максимум

при х

О

 

и быстро убывающая с ростом

х; / / f ~ ~ ~ j

имеет

смысл

час­

 

тотной

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

узкополосного фильтра

с

центральной

 

частотой

со. Произведение

U (х) H

 

^ ю j

 

представляет

собой

спектр

отфильтрованного

сигнала,

а

функции

\У (t,

со) |

и

arg1 У (£,

со) j — соответственно

его

 

огибающую

и

фазу

[29].

 

 

Рассмотрим,

как

с помощью функции | У

(t, со) | можно

опреде­

лить

групповую скорость поверхностной волны.

 

 

 

 

 

 

 

Как

следует из формулы (1.29), спектр &-й

гармоники

на

расстоянии

г

от

источника

можно

представить

в

виде

 

 

 

 

 

 

 

Ukq{x)

=

\Ukq(x)\exv

 

 

i(— ѵ(х)

+%?(*)) .

 

(5.3)

где

I

Uhq

I,

vhQ,

HpQ — медленно

меняющиеся функции

часто­

ты

X.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

В акустических исследованиях речи СВАН в аналоговой форме применяют­

 

ся уже более 20 лет [135, 153]; та же аппаратура

была использована эпизо­

 

дически

сейсмологами в 1959 г. [125], в СССР получил развитие аналого­

 

вый

метод

ЧИСС [8,

37],

в котором содержатся идеи

(но

не

техники)

СВАНА.

123


Применяя преобразование (5.2) и оценивая У (t, си) на больших расстояниях г методом стационарной фазы [44], по­ лучаем

C f t Q e xp \i .со t

~ - ) + ^ Q H + - ^яs i g n dCftQ (со

 

X

vkQ M

(5.4)

d C f t Q (со)

 

 

 

dû)

 

Здесь 5 = éo (t) — корень уравнения rit— ChQ (CÛ) = 0, где Сщ

(CÙ) —

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

для

простоты положим, что это уравнение имеет только один корень.

Максимум

функции H

~ 0 3

^ на плоскости [(со,

t) достигается

там,

где аргумент

(S со)/со

обращается в нуль, т. е. на

линии

со =

5

(t).

Поскольку

функции

I UhQ

(со)|, C^Q

И dChQld(ù ме­

няются

сравнительно

медленно,

область

повышенных

значе­

ний

функций

I У

(t, со) I гребень

рельефа

этой

функции

также будет приурочена к этой

линии. Зная со (t), легко

опре­

делить

групповую

скорость. Этим и

обосновывается возмож­

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

самым

групповой

скорости

по рельефу

функции

| У

(t, со) |.

Сами

значения | У (t, со) |

вдоль гребня определяют закон ам­

плитудной модуляции, а значения arg У

(t, со) содержат

инфор­

мацию

о фазовой

скорости.

 

 

 

 

Если сигнал и

(t) содержит несколько гармоник (или состоит из

нескольких элементарных

сигналов, соответствующих

разным

участкам кривой групповой скорости одной и той же

гармоники),

они будут проявляться в рельефе функции | У (t, со) | как от­ дельные гребни. Это позволяет разделить интерферирующие гармоники и определить их групповые скорости. Возможность разделения сигналов лимитируется разрешающей способно­

стью анализа, зависящей от ширины полосы пропускания

фильт-

Іх — со \

от длительности

сигнала.

 

 

 

 

 

 

ра H I ———1 и

 

 

 

 

 

 

Разрешающая

способность

СВАН.

Хорошо

известно,

что

частота и время являются сопряженными понятиями и их

невоз­

можно измерять

одновременно

сколь

угодно

точно.

Если

мы

определим разрешающую способность

по

времени

At

и

частоте

Лео одним из принятых способов [88],

то

принцип неопределен­

ности можно выразить неравенством

Асо At

=

qx

q >

0,

где

q зависит от способа определения

Асо и

At.

Поэтому,

когда

мы

говорим об амплитуде или фазе составляющей сигнала на часто­

те

со в момент времени t, фактически

речь идет о полосе частот

со

+ Асо и интервале времен t + At.

Таким образом, мы можем

124