Вестник СамГУ — Естественнонаучная серия. 2008. №2(61).
ВЫЧИСЛИТЕЛЬНАЯ МАТЕМАТИКА
УДК 550.344.5
МГНОВЕННЫЕ ПОЛЯРИЗАЦИОННЫЕ СВОЙСТВА ТРЕХКОМПОНЕНТНЫХ СИГНАЛОВ В БАЗИСЕ
ВЕЙВЛЕТОВ1
© 2008 М.А. Кулеш2
В данной работе рассматривается задача определения поляризационных параметров сейсмических сигналов, представленных на трехкомпонентной сейсмограмме. Предложен новый метод определения спектрально-временных поляризационных параметров в случае объемно-эллиптической поляризации сигнала. Использование в предложенном методе непрерывного вейвлет-пре-образования как одного из возможных подходов к спектрально-временному анализу позволяет обрабатывать и интерпретировать сейсмические волновые поля, для которых характерно большое количество волн различного типа с различным частотным составом, параметрами поляризации и протяженными зонами их интерференции. Предложенный метод может использоваться для построения эффективных алгоритмов фильтрации, что продемонстрировано на примере фильтрации экспериментальной сейсмограммы.
Введение
Конечной целью обработки сейсмограмм является измерение характеристик полезного сигнала в ситуации, когда там представлена достаточно сложная суперпозиция самых различных типов волновых движений. Например, в задачах пассивного сейсмического мониторинга имеется потребность обнаружения вступлений P-и S-волн [1]. В связи с тем, что волны Рэлея являются диспергирующими в гетерогенной среде, их поляризационные свойства являются частотно-зависимыми [2], что может быть использовано для дисперсионного анализа. В активных сейсмических экспериментах, напротив, волны Рэлея рассматриваются в качестве шума и должны быть подавлены с целью выделения именно объемных волн [3]. Другим приложением поляризационного анализа является идентификация эффекта расщепления S-волн, которое несет в себе информацию об анизотропных свойствах среды [4].
В дополнение к широко распространенным в волновой динамике методам идентификации волн по их частотному составу, скоростям или направлениям подхода волновых фронтов, также получили развитие методы селекции волн по характеристикам их поляризации, т.е. по чисто геометрическим характеристикам колебаний [5,6]. Например, применение поляризационных методов анализа позволило
1 Представлена доктором технических наук, академиком РАН В.П.Матвеенко.
2Кулеш Михаил Александрович ([email protected]), Институт механики сплошных сред УрО РАН, 614013, Россия, г.Пермь, ул. Акад. Королева, 1.
значительно повысить разрешающую способность методов скважинной и наземной сейсмики в сложнопостроенных средах [7-9].
К сожалению, используемые сегодня методы поляризационного анализа [10-13] имеют значительные ограничения, так как они не различают волны, пришедшие в один и тот же момент, но имеющие разные частоты. В связи с этим хорошие результаты дает использование спектрально-временного представления, замечательного, в частности, тем, что оно позволяет интерпретировать и изучать частотно-временные образы сигналов. Спектрально-временной подход получил мощное развитие в работах А.Л.Левшина и А.В.Ландера [6] и по праву стал одним из важнейших методов численного анализа сигналов. Сегодня разработано и активно используется большое количество различных методов спектрально-временного анализа. Среди них можно выделить оконное Фурье-преобразование, преобразование Габора, распределение Вагнера-Вилла [14], S-преобразование [15], дискретное вей-влет-преобразование [16,17], непрерывное вейвлет-преобразование [18,19]. Во всех этих методах общим является то, что результат анализа некоего сигнала должен содержать в себе не только перечисление его характерных частот, но и сведения об определенных временных координатах, в которых эти частоты проявляют себя.
В [6] можно найти, пожалуй, первую попытку интегрирования поляризационного и спектрально-временного анализа. С тех пор было опубликовано большое количество работ, посвященных подобной интеграции. В [20] используется метод сингулярного разложения матрицы ^УВ) совместно с мультивейвлетным анализом на базе дискретных ортогональных вейвлетов, а в [21] —на базе непрерывного вейвлет-преобразования. В [22] развивается метод, аналогичный методу Морозова [11], но с использованием S-преобразования. Комбинация комплексного ковариационного метода и непрерывного вейвлет-преобразования с использованием комплексных ортогональных вейвлетов Майера приведена в [23]. И, наконец, оценки некоторых поляризационных параметров с использованием дискретного вейвлет-пакетного разложения предложены в [1].
В данной работе предложен новый метод определения спектрально-временных поляризационных параметров в случае объемно-эллиптической поляризации сигнала на основе непрерывного вейвлет-преобразования. На практике предпочтение зачастую отдается дискретным вейвлетам [16,17]. Это связано с тем, что обычно используемые базисы на основе непрерывных вейвлетов не являются, строго говоря, ортонормированными. С дискретными вейвлетами таких проблем не возникает. В силу этого дискретные вейвлеты приводят обычно к более точному преобразованию и представлению сигнала, в особенности, к его обратному восстановлению. Однако предложенный в работе метод фактически базируется на свойстве избыточности непрерывного вейвлет-преобразования, связанного с непрерывным изменением временного и частотного параметров. Предложенный в работе метод может использоваться для построения эффективных алгоритмов фильтрации, что продемонстрировано на примере фильтрации экспериментальной сейсмограммы.
1. Основные соотношения вейвлет-анализа сигналов
В отличие от Фурье-преобразования, вейвлет-преобразование обеспечивает двухмерную развертку исследуемого сигнала, при этом частота и время рассматриваются как независимые переменные. В результате появляется возможность ана-
лизировать свойства сигнала одновременно в физическом (время) и масштабном (частота) пространствах.
Непрерывное вейвлет-преобразование определяется следующим образом [18,19]. Пусть s(t) — вещественный или комплексный сигнал, удовлетворяющий условию конечности энергии:
s(t) е L2(R) ^ I |s(t)|2Jt < то.
Пусть далее g(-), h(-) е L2(R) — вещественные или комплексные базисные функции специального вида, называемые вейвлетами, звездочка сверху обозначает комплексное сопряжение, а Wgs(t, a) — комплексный вейвлет-образ, зависящий от двух факторов: безразмерного масштабного a е R и размерного временного t е R. Тогда вейвлет-преобразование сигнала s(t) является его сверткой со множеством сдвинутых и скалированных вейвлетов и имеет вид:
+то +то
Wgs(t, а) = {TtDag, s) = J -ag* ) я(т) = ei(otg*(am)s(со) dm.
Операция Tt : g(t) м g(t -1) обозначает временной сдвиг материнского вейвлета, а Da : g(t) м g(x/a)/a — его масштабирование. Необходимо отметить, что вышеприведенное определение операции Da не является единственным. Например, часто используется также масштабирование вида Da : ^(т) н-> g(x/a)/ л/а с соответствующей модификацией формулы для вейвлет-преобразования. Однако в этом случае амплитуда вейвлет-коэффициентов становится пропорциональной не только амплитуде самого сигнала на том или ином масштабе, но также и масштабному параметру.
Обратное вейвлет-преобразование определяется следующим образом:
+то +то
s(t) = Mh'WgSit,^ = —!— Г f Xh(-—-\rWgs(x,a)dxda.
Cg,h J J a2 \ a I
-CO — CO
В нашем случае в качестве масштаба удобно выбрать физическую частоту, измеряемую в Гц, которая связана с безразмерным фактором: а = /о//, где /о — характерная (средняя) частота вейвлета. Тогда
+то +то
—то —то
+ TO + TO
s(t) = Mh<Wgs(t,f) = J-JJ ±h^^V-]jwgs(x,f)dxdf.
— CO —CO
(1.1)
В приведенных выше соотношениях — амплитудный коэффициент, аналогичный коэффициенту 1/2п обратного преобразования Фурье:
+ TO
CgJ% = J" — (g*(w)k(w) + g*(-w)k(-w)) d(o.
о
Достаточно подробные обзоры вейвлетообразующих функций приведены в [18,19], мы же ограничимся здесь только теми, которые будем использовать в дальнейшем (рис. 1):
Рис. 1. Характерный вид вейвлетов Морле (а) и Пауля (Ь) во временном (а), (Ь) и частотном (с) пространствах
1) Комплексный вейвлет Морле, хорошо локализованный во временном и частотном пространствах, приведен на рис. 1, а. Сплошной линией изображена действительная часть вейвлета, а пунктирной — мнимая. Фурье-спектр вей-влета Морле является действительным и изображен на рис. 1, с сплошной линией. Значения самого вейвлета и его Фурье-спектра определяются соотношениями:
g(t) = еы,е-,2'аа2), g(w) = а л/2же-{ю-2яГ'а2/2.
Ширина вейвлета в частотном пространстве определяется коэффициентом о, а средняя частота этого вейвлета f° равна единице.
2) Комплексный вейвлет Пауля приведен на рис. 1, b (сплошная линия—действительная часть, пунктирная —мнимая):
Средняя частота этого вейвлета в Фурье-пространстве также равна единице (рис. 1, с, пунктирная линия), а ширина определяется параметром p. В отличие от вейвлета Морле, Фурье-спектр вейвлета Пауля является несимметричным относительно средней частоты.
Вейвлеты Морле и Пауля являются прогрессивными, то есть спектр Фурье этих вейвлетов равен нулю для отрицательных частот. Это свойство позволяет представить вейвлет-спектр в виде суперпозиции прогрессивного W+ s(t, f) и регрессивного W-s(t, f) спектров,
- I e-(P-D«/2n
где
Wgs(t, f) - W+ s(t, f) + W- s(t, f),
W+ s(t f) - { f)' f > °' w-sit f) - i °' f > °'
Wg s(t- f) - j f < o, Wg s(t'f) - j WgS(t, f),
Для вещественных сигналов модули прогрессивного и регрессивного вейвлет-спектров являются симметричными относительно нулевой частоты. Однако, как уже отмечалось, сигнал s(t) может быть и комплексным. В этом случае прогрессивная и регрессивная составляющие вейвлет-спектра будут различаться не только по фазе, но и по модулю.
2. Адаптивный ковариационный анализ трехкомпонентного сигнала
Рассматривая сигнал ЯГ) = |5х(0> яу(Г), (0)г, соответствующий трехкомпонент-ной записи, где 5х(0, .у(0 и являются радиальной, трансверсальной и вертикальной компонентами, как все три, так и любая комбинация двух ортогональных компонент могут быть выбраны для поляризационного анализа. В каждый момент времени такой сигнал наилучшим образом (в смысле наименьших квадратов) аппроксимируется прямой (линейно-поляризованные колебания), эллипсом (эллиптически-поляризованные колебания), сферой или трехосным эллипсоидом (объемно-поляризованные колебания). Линейно-поляризованные колебания являются частным случаем эллиптически-поляризованных, которые, в свою очередь, являются подмножеством объемно-поляризованных колебаний. Для описания объемно-поляризованного сигнала могут использоваться следующие параметры:
1) большой поляризационный вектор К е К3 и большая полуось я = ||Я||,
2) средний поляризационный вектор ?в е К3 и средняя полуось г8 = || е [0, К],
3) малый поляризационный вектор Ге К3 и малая полуось г = ||Г|| е [0, г8],
а также ряд специализированных параметров:
1) эксцентриситет эллипса р = г8 /К, р е [0, 1],
2) эксцентриситет эллипсоида р1 = г/г8, р1 е [0, 1],
3) угол падения в е (-я/2, я/2] между большой полуосью и вертикальной осью г,
4) угол наклона 6 е (-я/2, п/2] между большой полуосью и горизонтальной осью х,
5) азимут у е (-я, я], который является углом между горизонтальной осью х и проекцией большой полуоси на горизонтальную плоскость.
К сожалению, математически точного априорного определения мгновенных поляризационных атрибутов трехкомпонентного сигнала не существует, поэтому любая попытка так или иначе определить поляризационные свойства всегда будет базироваться на тех или иных допущениях. В работе [11] предложен метод, в котором трехмерные движения частиц среды моделируются как плоские эллиптические движения. Однако в более общем случае подобные движения можно также моделировать трехмерной кривой, опоясывающей трехмерный эллипсоид с тремя поляризационными векторами К е К3, Г е К3 и е К3. Ковариационный метод, описанный в [6,12,13,24,25], базируется именно на такой модели сигнала.
В основе ковариационного метода лежит анализ собственных чисел и собственных векторов ковариационной матрицы [М©], вычисляемой следующим образом:
где
М© =
1+Г/2
Мхх© Мху © Их;® Мху£) Иуу © Иу;©
Их;©) Му&) И;;©
м]т© = ^ ^ о/т) - ^©)(5т(т) - цт©Мт, Л т = (х,у,г).
1-Г/2
Данная матрица вычисляется для временного окна T, где величина ^ — это центр окна. Величина цj© характеризует среднее арифметическое компоненты Sj во временном окне T.
Три собственных значения Xi© a ^2© a ^3© ковариационной матрицы [M©] и три собственных вектора Vj(^) полностью характеризуют величину и направления осей мгновенного эллипсоида, который наилучшим образом аппроксимирует движения частиц в интервале T. Это позволяет вычислить три вектора поляризации:
- большой полуоси $© = y?li©i?i©/||i?i©||,
- средней полуоси Fs© = л/^2©^2©/11^2©11,
- малой полуоси ?© = дДз©^©/^©!!.
Проблемной точкой ковариационного метода является выбор окна T. В ситуации, когда в момент времени ^ сигнал представлен суперпозицией различных волн с разной частотой, однозначных критериев выбора этого окна не существует. Естественным решением этой проблемы является использование спектрально-временного представления [6]. Однако при этом появляется проблема автоматической адаптации временного окна T в соответствии с частотным масштабом. В данной работе предложен адаптивный ковариационный метод, позволяющий эффективно решить эту проблему [26].
Будем рассматривать аналитический сигнал z(t), вычисленный покомпонентно: z(t) = SS(t) + iHsS(t), где операция H означает преобразование Гильберта:
HS(f) = { jf /Ca0' (2.1)
Тогда каждая компонента исходного сигнала sj(t) может быть аппроксимирована следующим образом:
Sj(t + т) = s+ (t + т) + sj(t + т) = (zj(t + т) + z*(t + т)) /2 ^
= (zj(t)ei^(t)x + z*(t)e~i0'«т) /2 = (2.2)
= |zj(t)| cos (Оj(t^ + argZj(t)),
где Оj(t) = д argzj(t)/dt является мгновенной частотой аналитического сигнала.
Используя аппроксимацию (2.2), можно получить два результата. Во-первых, интегралы в выражении для элементов ковариационной матрицы могут быть вычислены аналитически:
Mjm(t) = kj(0lkm(0l jsinc (^f^ Atjm(0) COS (Ajm(t)) +
+ sine (aj(t>+2a'"(t> Atjm(t))COS (A+m(0)J - \ljm\lmj, (2.3)
A±m(t) = arg Zj(t) ± argzm(t),
цjm = ft(zj(t))sinc(Atjm(t)Qj(t)/2), j, m = x,y, z,
где sinc(x) = sin(x)/x.
Во-вторых, несложно заметить, что для идеальноэллиптического сигнала, когда Ox(t) = Oy(t) = Oz(t) = O(t), аппроксимация (2.2) является гармонической с периодом 2n/0(t). Отсюда несложно показать, что при выборе окна Atjm(t) = 2nn/0(t) мы получим значения полуосей эллипса точно. В общем случае мгновенные частоты всех компонент сигнала различаются, но при этом мы можем выбирать в
качестве базовой мгновенной частоты среднее значение. Это позволяет определить адаптивное окно Atjm(t) следующим образом:
2nn 6пп
Atjm(t) =
flx(t) + fly(t) + flZ(t)
(2.4)
Atjm(t) =
2пп
iC«
4яп
О j(t) + Qm(t)'
n e N.
(2.5)
В выражениях (2.4)—(2.5) присутствует целочисленный параметр п, означающий, сколько периодов будут браться в аппроксимации (2.2) для вычисления элементов ковариационной матрицы. На рис. 2 показано, что происходит при изменении этого параметра. На рис. 2, а изображены три компоненты (Г), некоторого сигнала (сплошная линия) и аппроксимация (2.2) этих компонент (пунктирная линия). Мгновенные частоты каждой компоненты исходного сигнала различаются. На рис. 2, Ь-й показана аппроксимирующая кривая для трех разных адаптивных окон, вычисленных при помощи соотношения (2.4) для трех значений параметра п.
или
(b) n=1, At=0.04884
t, (s)
(c) n=2, At=0.10012
(d) n=3, At=0.1514
N 0
Y -5 -10 X
Y -5 -10 X
Y -5 -10 X
0
5
Рис. 2. Иллюстрация влияния параметра п в выражении (2.4) на точность аппроксимации адаптивного ковариационного метода
Можно заметить, что все эти кривые имеют соответствующее число практически полных периодов. При увеличении параметра п профиль кривых становится более объемным, повышается чувствительность метода с точки зрения определе-
ния малой оси поляризационного эллипсоида. В то же время при этом увеличивается отличие между исходным сигналом и аппроксимирующей функцией, то есть понижается точность определения остальных поляризационных осей.
Предложенный метод имеет два частных случая. Во-первых, рассмотрим двух-компонентный сигнал = [ях(Г), яу(Г)}т, у которого мгновенные частоты обеих компонент равны: Ох(Г) = Пу(0 = О(Г). В этом случае ковариационная матрица имеет вид:
M(t) =
о I Zy (О
\zAt)\2 \zx{t)\\zy{t)\ cos arg --—
1 Zx(t)
kx(0llzy(0|cos(arg^l |Zy(0P
zx(t) /
Несложно найти собственные значения этой матрицы: ^1,2(0 = \{\zx(t)\2 + \Zy(t)\2±
± Vd^WI2 - MOI2)2 + (2кх(01ку(0|со8(Аф(0))2),
где Aф(t) = arg(zX(t)zy(t)). Эти два значения в точности соответствуют поляризационным полуосям двухкомпонентного сигнала, получаемым в методе комплексного следа [10].
Во-вторых, рассмотрим трехкомпонентный сигнал также в предположении равенства мгновенных частот всех компонент. В этом случае элементы ковариационной матрицы определяются выражением:
Mjm(t) = |zj(t)||zm(t)| cos (argZj(t) - argZm(t)) = %(zj(t)zm(t)).
Задавая аналитически тестовый гармонический сигнал с полуосями R и r, собственные значения этой матрицы также несложно отыскать аналитически:
Xi(t) = R2, X2(t) = r2, X3(t) = 0.
Таким образом, предположение равенства мгновенных частот приводит к эллип-тически-поляризованному сигналу в трехмерном пространстве, и в этом случае результат, полученный новым адаптивным ковариационным методом, в точности совпадает с результатом метода [11].
Предложенный адаптивный ковариационный метод обладает рядом преимуществ:
1) Во-первых, решена проблема выбора временного окна вычисления элементов ковариационной матрицы.
2) Во-вторых, элементы ковариационной матрицы можно вычислять при помощи аналитической формулы, не содержащей суммы или интегралов. Это положительно сказывается на производительности метода.
3) В-третьих, метод является обобщением известных и хорошо зарекомендовавших себя метода комплексного следа [10] и метода Морозова и Смитсона [11].
4) И, наконец, использование адаптивного временного окна позволяет обобщить данный метод на случай его использования для вычисления поляризационных параметров в частотно-временном пространстве, что будет подробно рассмотрено далее.
3. Частотно-временной адаптивный ковариационный метод
Можно обобщить описанный выше метод (2.3) на случай частотно-временного представления. В данном случае такое обобщение допустимо благодаря возможности выбора окна интегрирования адаптивным образом, что позволяет автоматически находить оптимальное окно для каждой частоты вейвлет-спектра. Численные тесты показали, что более предпочтительным является использование соотношения (2.5), так как оно дает более сглаженную картину поляризационных свойств:
AtJm(t,f) = ^ , ^ п е N,
4пп
iW)Töm(i,/)' (3.1)
О j(t, /) = д arg W+ Sj (t, /)/dt.
Сравнивая свойство (1.2) вейвлет-преобразования с соотношениями (2.1), несложно заметить, что прогрессивный вейвлет-спектр, вычисленный при помощи прогрессивного вейвлета, соответствует вейвлет-спектру аналитического сигнала. Это свойство позволяет совершенно естественным образом записать аналог аппроксимации (2.2) для непрерывного вейвлет-спектра. Для вычисления элементов ковариационной матрицы используются только вещественные части прогрессивных вейвлет-коэффициентов, при этом, как и ранее, используем разложение вейвлет-коэффициентов в ряд Тейлора по времени в окрестности точки (t, /):
%(W+ Sj(t + т, /)) - |W+ Sj(t, /)| cos (Оj(t, /)т + arg W++ Sj(t, /)).
Используя данное разложение, соотношения для вычисления элементов ковариационной матрицы по аналогии с (2.3) примут вид [27-29]:
Mjm(t, /) = |W+ Sj (t, /)|| W++ Sm (t, /)|Х
( i Оj(t, f) j Om(t, f) \
xjsmc I —---A tjm(t,f)\ cos (A jm(t, f))+
. i О j(t, f) + Qm(t, Л , \ t+t Л
+smc M---Atjm(t, f) cos (Ajm(t, /))|- (3 2)
— ц jm (t, /)Цmj(t, /),
A±m(t, f) = arg W+ Sj(t, f) ± arg W+ Sm(t, f), j, m = x, y, z,
Ц jm(t, f) = Ä(W+ Sj(t, / ))sinc(Atjm(t, /)О j(t, f)/2),
где sinc(x) = sin(x)/x.
Три собственных значения Xi(t, f) a ^2(t, f) a ^3(t, f) ковариационной матрицы [M(t, f)], а также три собственных вектора Vj(t, f) полностью характеризуют величину и направления осей мгновенного эллипсоида в момент времени t на частоте f, который аппроксимирует движения частиц в интервале Atjm(t, f). Это позволяет вычислить три мгновенных вектора поляризации:
большой полуоси Ü(t,f) = Л1Ш7)ШЛ/\\ШЛ\\, средней полуоси ?s(t,f)= ^X2(t,f)l¥i{t,/)II, малой полуоси ?(t,f) = ф^ЛЫ', Л/Ш', Л\\-
Предложенный метод (3.2) обладает большим потенциалом по фильтрации сигналов с использованием тех или иных поляризационных параметров, вычисляемых с использованием мгновенных векторов поляризации. В основе такой фильтрации лежат два свойства. Во-первых, определяющим является возможность использовать обратное вейвлет-преобразование после манипуляций с коэффициентами вейвлет-спектра, чего нет, например, у оконного Фурье-преобразования. Во-вторых, предложенные поляризационные методы позволяют определить большое количество различных параметров, которые можно использовать для идентификации тех или иных типов волн в частотно-временной плоскости. Поэтому, объединяя эти два свойства, можно определить поляризационный фильтр следующим образом:
где N и Р — количество и множество задействованных поляризационных параметров, Рт — формальное обозначение одного какого-либо параметра, а Рт — множество допустимых значений данного параметра, соответствующее физическим свойствам волны, которую требуется отфильтровать.
4. Фильтрация сейсмограммы регионального землетрясения
С помощью адаптивного ковариационного метода (3.2) выполнена фильтрация трехкомпонентной сейсмограммы (рис. 3), записанной во время регионального землетрясения [27]. Это землетрясение (время 13:02 13.03.1989, широта 50,7198 N долгота 9,9112 Е, к = 0,7 км, М^ = 5,6 [30]) было спровоцировано разрушением солевой шахты в Тюрингене, Германия на удалении 148 км и в направлении (с обратным азимутом) уо = 128,8° относительно сейсмической станции GRA1, входящей в состав сети GRF в северной Баварии.
Необходимо отметить, что сейсмограмма не является предварительно повернутой. Поэтому восточная и северная компоненты не совпадают с радиальной и трансверсальной компонентами распространяющихся волн, то есть при анализе угловых поляризационных параметров там, где это необходимо, нужно учитывать обратный азимут у о.
Рассматриваемая сейсмограмма начинается с Р-волны примерно на 165 с. Примерно на 20 с. позже появляются поперечные волны в горизонтальной плоскости. Так как источником явилось разрушение шахты, амплитуды SH-волн не являются значительными. В интервале 180-240 с. также наблюдается рэлеевская волна с вертикальной компонентой 5г(0. Для нее на рис. 3 можно заметить дисперсию на низких частотах 0,15 Гц для первичных волн и на высоких 0,6 Гц для вторичных, что согласуется со скоростным анализом рассматриваемого события, выполненно-
Для поляризационного анализа использовался метод (3.2) с параметром временного окна п = 3. Как результат применения адаптивного ковариационного метода, на рис. 4, Ь показан частотно-временной спектр коэффициента эллиптичности рх(Г, /) с учетом знака:
/ (0 = Мк£р /),
/), если Рт е Рт, т е [1, N1,
0 в противном случае,
го в [31].
рх =
J_I_I_I_I_I_I_I_I
(«)
Рис. 3. Трехкомпонентная сейсмограмма регионального землетрясения, содержащая восточную (?), северную (?) и вертикальную (?) компоненты, а также вейвлет-спектр всех компонент. Временная ось отсчитывается от момента времени 13:00 13.03.1989
Волновые пакеты в данных частотно-временных областях имеют регрессивный характер колебаний и поэтому могут быть однозначно идентифицированы как фундаментальная мода рэлеевской волны. Применяя обратное вейвлет-преобразо-вание к спектру, отфильтрованному предварительно при помощи отрицательного коэффициента эллиптичности при условии р5(?, /) < -0.15, можно получить сейсмограмму, где остаются только волны Рэлея. Такая сейсмограмма показана на рис. 4, с, где низкочастотная компонента рэлеевской фундаментальной моды с компонентами Z - N четко различима в интервале 190-210 с.
Интересным является также применение объемного направленного фильтра с использованием азимута у(?, /). Были выделены волновые пакеты с поляризацией, параллельной направлению обратного азимута, что характерно для продольных волн: (у(?, /) - уо) е [-20°, 20°]; они приведены на рис. 5, с пунктирными линиям. Кроме этого, там же сплошными линиями показаны волновые моды с доминирующим направлением поляризации, перпендикулярным обратному азимуту: (у(?, /) -_ Уо) е [70°, 90° ] и [-70°, -90°]. В качестве примера таких волн можно наблюдать SH волну Лява в интервале 200-210 с.
Заключение
В работе предложен новый метод определения спектрально-временных поляризационных параметров трехкомпонентных сигналов. Использование в предложенном методе вейвлет-преобразования как одного из возможных подходов к спек-
Рис. 4. Регрессивный поляризационный фильтр: (а)—исходная сейсмограмма; (Ь)—частотно-временной спектр отрицательного коэффициента эллиптичности; (с)—волновые пакеты, соответствующие регрессивному характеру колебаний
трально-временному анализу позволяет обрабатывать и интерпретировать сейсмические волновые поля, для которых характерно большое количество волн различного типа с разным частотным составом, параметрами поляризации и протяженными зонами их интерференции.
Преимуществами данного адаптивного ковариационного метода являются:
1) Возможность хорошей регулировки как частотно-временного разрешения за счет подбора параметра вейвлета, так и объемной чувствительности метода за счет изменения параметра п в соотношении (3.1).
2) Допустимость его обобщения на случай анализа сигналов с любым числом компонент.
3) Отсутствие проблемы выбора временного окна вычисления элементов ковариационной матрицы, так как оно выбирается автоматически.
4) Предложенный метод может использоваться не только для трехкомпонентно-го сигнала, но и для сигналов с большим числом компонент, однако физическая интерпретация как самого сигнала, так и его поляризационных свойств представляется в этом случае не всегда однозначной.
Литература
[1] Любушин, А.А. Вейвлет-пакетный поляризационный метод для автоматического детектирования вступлений P и S-волн / А.А. Любушин // Физика земли. - 2006. - №4. - C. 30-39.
[2] Shieh, C.-F. Ground roll: Rejection using polarization filters / C.-F.Shieh, R.B.Herrmann // Geophysics. - 1990. - V. 55. - No. 9. - P. 1216-1222.
[3] Wang, Y. Separation of P- and S-wavefield from wide-angle multicomponent OBC data for a basalt model / Y.Wang, S.C.Singh // Geophysical Prospecting. -2003. - V. 51. - P. 233-245.
[4] Li, X.L. Complex component analysis of shear wave splitting: theory / X.L.Li, S.Crampin // Geophysical Journal International. - 1991. - V. 107. - P. 597-604.
[5] Алказ, В.Г. Поляризационный анализ сейсмических колебаний / В.Г. Алказ,
H.И. Онофраш, А.И. Перельберг. - Кишинев: Штиинца, 1977. - 110 c.
[6] Левшин, А.Л. Поверхностные сейсмические волны в горизонтально-неоднородной земле / А.Л.Левшин, Т.Б.Яновская, А.В.Ландер. - М.: Наука, 1986. -278 c.
[7] Александров, С.И. Поляризационный анализ сейсмических волн / С.И.Александров. - М.: ОИФЗ РАН, 1999. - 142 c.
[8] Гальперин, Е.И. Поляризационный метод сейсмических исследований / Е.И.Гальперин. - М.: Недра, 1977. - 277 c.
[9] Perelberg, A.I. Applications of seismic polarization analysis / A.I. Perelberg, S.C. Hornbostel // Geophysics. - 1994. - V. 59. - No. 1. - P. 119-130.
[10] Multicomponent seismic studies using complex trace analysis / R.M.Rene [et al.] // Geophysics. - 1986. - V. 51. - No. 6. - P. 1235-1251.
[11] Morozov, I.B. Instantaneous polarization attributes and directional filtering /
I.B. Morozov, S.B. Smithson // Geophysics. - 1996. - V. 61. - No. 3. - P. 872-881.
[12] Flinn, E.A. Signal analysis using rectilinearity and direction of particle motion / E.A.Flinn // Proceedings of the IEEE. - 1965. - V. 53. - No 12. - P. 1874-1876.
[13] Kanasewich, E.R. Time Sequence Analysis in Geophysics / E.R. Kanasewich. -University of Alberta Press, Edmonton, Alberta, 1981.
[14] Prosser, W.H. Time-frequency analysis of the dispersion of Lamb modes / W.H. Prosser, M.D. Seale, B.T. Smith // The Journal of the Acoustical Society of America. - 1999. - V. 105. - No. 5. - P. 2669-2676.
[15] Stockwell, R.G. Localization of the complex spectrum: the S transform / R.G. Stockwell, L. Mansinha, R.P. Lowe // IEEE Transactions on Signal Processing. - 1996. - V. 44. - No. 4. - P. 998-1001.
[16] Дремин, И.М. Вейвлеты и их использование / И.М.Дремин, О.В.Иванов, В.А.Нечитайло // Успехи физ. наук. - 2001. - T. 171. - No. 5. - C. 465-561.
[17] Новиков, И.Я. Основы теории всплесков / И.Я.Новиков, С.Б.Стечкин // Успехи мат. наук. - 1998. - T. 53. - No. 6. - C. 53-128.
[18] Астафьева, Н.М. Вейвлет-анализ: основы теории и примеры применения / Н.М.Астафьева // Успехи физ. наук. - 1996. - T. 166. - №11. - C. 1145-1170.
[19] Holschneider, M. Wavelets: an Analysis Tool / M. Holschneider. Oxford: Clarendon Press, 1995.
[20] Claassen, J.P. Robust bearing estimation for three-component stations / J.P. Claassen // Pure and Applied Geophysics. - 2001. - V. 158. - No. 1. -P. 349-374.
[21] Time- and frequency-dependent polarization analysis: anomalous surface wave observations in Iberia / H. Paulssen [et al.] // Geophysical Journal International. -1990. - V. 103. - No. 2. - P. 483-496.
[22] Pinnegar, C.R. Polarization analysis and polarization filtering of three-component signals with the time-frequency S transform / C. R. Pinnegar // Geophysical Journal International. - 2006. - V. 165. - No. 2. - P. 596-606.
[23] Soma N. Reflection technique in time-frequency domain using multicomponent acoustic emission signals and application to geothermal reservoirs / N. Soma, H. Niitsuma, R. Baria // Geophysics. - 2002. - V. 67. - No. 3. - P. 928-938.
[24] Jurkevics, A. Polarization analysis of three-component array data / A. Jurkevics // Bulletin of the Seismological Society of America. - 1988. -V. 78. - No. 5. - P. 1725-1743.
[25] Jackson, G.M. Principal component transforms of triaxial recordings by singular value decomposition / G.M. Jackson, I.M. Mason, S.A. Greenhalgh // Geophysics. - 1991. - V. 56. - No. 4. - P. 528-533.
[26] Instantaneous polarization attributes based on an adaptive approximate covari-ance method / M. S.Diallo [et al.] // Geophysics. - 2006. - V. 71. - No. 5. -P. V99-V104.
[27] Polarization analysis in the wavelet domain based on the adaptive covariance method / M. Kulesh [et al.] // Geophysical International Journal. - 2007. -V. 170. - No. 2. - P. 667-678.
[28] Elliptic properties of elastic surface waves in wavelet domain / M. Kulesh [et al.] // Proceedings of the XXXIII International Summer School "Advanced Problems in Mechanics". - 2005. - P. 361-366.
[29] Estimating polarization attributes with an adaptive covariance method in the wavelet domain / M. S. Diallo [et al.] // SEG Technical Program Expanded Abstracts. - 2005. - P. 1014-1017.
[30] On-line Bulletin International Seismological Centre // http://www.isc.ac.uk, In-ternatl. Seis. Cent., Thatcham, United Kingdom. 2001.
[31] Hadiouche, O. Mapping the crust in southeastern Germany using Rayleigh waves in the period range 6-16 s / O.Hadiouche, F.Kriiger, R. Kind // Geophysical Research Letters. - 1991. - V. 18. - No. 6. - P. 1087-1090.
Поступила в редакцию 13/Л7/2008; в окончательном варианте — 13/Л7/2008.
INSTANTANEOUS POLARIZATION PROPERTIES OF THREE-COMPONENT SIGNALS IN WAVELET
DOMAIN3
© 2008 M.A.Kulesh4
In the paper the problem of the determination of polarization properties of seismic signals presented on a three-component seismogram is considered. A new method that allows calculating the time-frequency polarization attributes in the case of a three-dimensional elliptical polarization is proposed. The continuous wavelet transform used by the proposed method as one of the possible approaches of the time-frequency analysis allows us to process and interpret the seismic wave fields characterized by a superposition of different waves with various spectral and polarization properties as well as with extensive interference zones. The proposed method can be used to construct the effective filtering algorithms that we demonstrate by an example of an experimental seismogram filtering.
Paper received 13/1/7/2008. Paper accepted 13/777/2008.
3Communicated by Dr. Sci. (Tech.), Acad. of RAS, Prof. V.P. Matveenko.
4Kulesh Mikhail Alexandrovich (kmaaicmm.ru), Institute of Continuum Mechanics, Ural Division, Russian Academy of Sciences, Perm, 614013, Russia.