Файл: Левшин А.Л. Поверхностные и каналовые сейсмические волны [монография].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. Спектрально-временной анализ поверхностных волн
Приведенные теоретические сейсмограммы могут служить модельными примерами исходных данных, используемых в сей смологии при решении обратной задачи — определении строения среды по наблюденному волновому полю. Вид сейсмограмм в зна чительной степени определяется параметрами источника, варьи рующими от землетрясения к землетрясению и известными ин терпретатору с малой точностью и надежностью (механизм, спектр, глубина очага). Поэтому наилучший путь интерпретации — в^ разделении влияния среды и источника, например путем опреде ления по имеющимся сейсмограммам инвариантных относитель но источника характеристик поверхностных волн. Согласно гл. 1—3, к ним относятся дисперсионные кривые 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, со) | как от дельные гребни. Это позволяет разделить интерферирующие гармоники и определить их групповые скорости. Возможность разделения сигналов лимитируется разрешающей способно
стью анализа, зависящей от ширины полосы пропускания |
фильт- |
||||||||||
1Т Іх — со \ |
от длительности |
сигнала. |
|
|
|
|
|
|
|||
ра H I ———1 и |
|
|
|
|
|
|
|||||
Разрешающая |
способность |
СВАН. |
Хорошо |
известно, |
что |
||||||
частота и время являются сопряженными понятиями и их |
невоз |
||||||||||
можно измерять |
одновременно |
сколь |
угодно |
точно. |
Если |
мы |
|||||
определим разрешающую способность |
по |
времени |
At |
и |
частоте |
||||||
Лео одним из принятых способов [88], |
то |
принцип неопределен |
|||||||||
ности можно выразить неравенством |
Асо At |
= |
qx |
q > |
0, |
где |
|||||
q зависит от способа определения |
Асо и |
At. |
Поэтому, |
когда |
мы |
говорим об амплитуде или фазе составляющей сигнала на часто
те |
со в момент времени t, фактически |
речь идет о полосе частот |
со |
+ Асо и интервале времен t + At. |
Таким образом, мы можем |
124