Разработка методов частотно-временного анализа поляризационных и дисперсионных свойств волновых процессов тема диссертации и автореферата по ВАК РФ 05.13.18, доктор физико-математических наук Кулеш, Михаил Александрович

• Кулеш, Михаил Александрович
• доктор физико-математических наук
• 2008

Оглавление диссертации доктор физико-математических наук Кулеш, Михаил Александрович

Обозначения и сокращения

Введение

1 Волновая динамика упругих сред

1.1 Краткая характеристика акустических волн в твердых телах

1.1.1 Объемные волны.

1.1.2 Поверхностные волны.

1.1.3 Волноводные и канализированные волны.

1.2 Плоские волны в упругой классической среде.

1.2.1 Общее решение уравнений движения для полупространства

1.2.2 Дисперсионные уравнения волн Рэлея, Стоунли, Лява и Лэмба.

1.2.3 Дисперсионные и поляризационные параметры волн

1.3 Волновые процессы в реальных средах.

1.3.1 Связь дисперсии и диссипации.

1.3.2 Геометрическая дисперсия и нелинейность

1.3.3 Анизотропия.

1.3.4 Модели микрополярного континуума.

1.4 Основные соотношения вейвлет-анализа сигналов.

1.4.1 Вейвлеты.

1.4.2 Представление вейвлет-спектра.

2 Поляризационный анализ волновых процессов

2.1 Определение поляризационных параметров во временном пространстве.

2.1.1 Метод комплексного следа.

2.1.2 Модификация метода комплексного следа.

2.1.3 Эллиптическая аппроксимация трехкомпонентного сигнала.

2.1.4 Адаптивный ковариационный анализ трехкомпонентного сигнала.

2.2 Определение поляризационных параметров в частотном пространстве

2.3 Частотно-временной поляризационный анализ.

2.3.1 Метод комплексного следа в вейвлет-пространстве

2.3.2 Эллиптическая аппроксимация трехкомпонентного сигнала в вейвлет-пространстве.

2.3.3 Частотно-временной адаптивный ковариационный метод.

2.4 Идентификация и фильтрация сигналов на основе поляризационного спектра.

3 Дисперсионный анализ волновых процессов

3.1 Моделирование распространения волн на основе преобразования Фурье.

3.2 Моделирование распространения волн с использованием асимптотического вейвлет-оператора.

3.2.1 Случай независимой дисперсии и диссипации

3.2.2 Случай связанной дисперсии и диссипации.

3.2.3 Дисперсионный оператор для поляризационных параметров

3.3 Корреляционный частотно-скоростной анализ пространственно-временной сейсмограммы.

3.4 Задача нелинейной оптимизации параметров дисперсионного оператора.

4 Программный комплекс поляризационного h дисперсионного анализа

4.1 Структура и технология реализации программного комплекса

4.1.1 Библиотечный уровень.

4.1.2 Уровень модулей.

4.1.3 Интерфейсный уровень.

4.2 Базовые модули.

5 Использование программного комплекса для обработки экспериментальных данных

5.1 Вейвлет-анализ вибродиагностических данных.

5.1.1 Вибродиагностический эксперимент на автомобильном мосту.

5.1.2 Вибродиагностический эксперимент в жилом доме

5.2 Поляризационный анализ геомагнитных данных.

5.3 Фильтрация двухкомпонентной сейсмограммы с использованием поляризационного фильтра.

5.4 Поляризационный анализ и фильтрация трехкомпонентных сейсмограмм.

5.4.1 Фильтрация сейсмограммы ударного эксперимента

5.4.2 Фильтрация сейсмограммы регионального землетрясения

5.5 Корреляционный анализ экспериментальной сейсмограммы

5.6 Определение дисперсионных кривых по экспериментальной сейсмограмме.

Введение диссертации (часть автореферата) на тему «Разработка методов частотно-временного анализа поляризационных и дисперсионных свойств волновых процессов»

Очень многие области науки и техники связаны с волновыми процессами. Понятие волнового процесса является слишком общим, охватить все его аспекты в рамках одной работы достаточно сложно, да и в этом нет необходимости, так как существует большое количество монографий и обзоров по этой теме. Поэтому в рамках данной работы за основу берется только один класс таких процессов, а именно волновая динамика упругих сред. Одним из важных направлений теории волн в упругих средах является теория сейсмических процессов [1], которые можно рассматривать как волновые процессы планетарного или регионального масштабов. Практическими приложениями теории сейсмических процессов являются геофизические методы изучения и разведки земной коры [2], проблемы анализа и предсказания землетрясений [3], сейсмические методы контроля техногенных процессов [4,5].

В СССР исследования по сейсмической разведке были начаты еще в Институте теоретической геофизики АН СССР в 1938 году. В связи с необходимостью развития минерально-сырьевой базы страны, проблемы создания новых эффективных методов сейсмической разведки рассматривались в 40-х - 60-х годах в числе важнейших проблем АН СССР. Задачи повышения точности существующих методов, а также развитие новых методов, позволяющих учитывать большее количество эффектов, не утратили своей актуальности и сегодня. По своей сути все сейсмические методы изучения свойств и состава глубоких слоев Земли базируются на интерпретации характера преломлений, отражений и изменений скорости прохождения сейсмических волн, возникающих при землетрясениях или при искусственных взрывах.

Другими интересными в инженерном плане приложениями волновой динамики упругих сред являются вибродиагностические исследования. По сравнению с анализом сейсмических волновых процессов, в основе изучения вибродиагностических данных лежит определение не свойств распространяющихся волн, а их воздействия на изучаемый объект [6,7]. Длительные вибрационные воздействия в значительной степени определяют надежность и долговечность сооружений и конструкций. В связи с этим развитие методов вибродиагностики является чрезвычайно важной и актуальной задачей для теоретического и практического изучения и прогнозирования надежности и долговечности.

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

Краткий обзор используемых сегодня моделей волновой динамики упругих сред приведен в первой главе диссертационной работы. Наиболее распространенной и изученной, по-видимому, является модель плоской волны в упругой изотропной среде, которая нашла широкое применение, как в научных исследованиях, так и в промышленных приложениях. Основные закономерности поведения таких волн изложены, например, в работах [8-10]. При распространении упругих волн в приповерхностных геологических слоях наблюдается связь диссипации энергии и дисперсии плоских волн, что учитывается в т.н. (^-моделях [1,11]. Глубоко проработана и широко применяется концепция горизонтально-слоистых (или вертикально-неоднородных) моделей строения осадочных толщ земной коры [12]. В начале 70-х годов была доказана возможность и целесообразность применения в сейсмологии и сейсморазведке моделей случайно-неоднородных сред. Развитие методов акустического каротажа скважин привело к пониманию необходимости построения тонкослоистых моделей осадочных толщ земной коры [13,14]. Развивая представления микрополярного континуума, учитывающие одновременно и трансляционные смещения, и кинематически независимые микроповороты, оказалось возможным в рамках одной модели описать процессы распространения тектонических уединенных волн, излучающих упругие сейсмические предвестники [15]. Среди последних достижений в области развития физических основ сейсмических методов необходимо отметить исследования по нелинейной сейсмике [16], в ходе которых были получены принципиально новые свидетельства о проявлении нелинейных сейсмических свойств реальных сред.

Конечной целью обработки сейсмограмм и виброграмм является измерение характеристик полезного сигнала в ситуации, когда там представлена достаточно сложная суперпозиция самых различных типов волновых движений. В дополнение к широко распространенным в волновой динамике методам идентификации волн по их частотному составу, скоростям или направлениям подхода волновых фронтов, также получили развитие методы селекции волн по характеристикам их поляризации, т.е. по чисто геометрическим характеристикам колебаний [13,17-20]. Например, отличительные особенности поляризационных свойств сейсмических волн различной природы детально рассматривались в фундаментальной монографии [21].

В задачах пассивного сейсмического мониторинга имеется потребность обнаружения вступлений Р- и 8-волн [22,23]. В связи с тем, что волны Рэлея являются диспергирующими в гетерогенной среде, их поляризационные свойства являются частотно-зависимыми [17], что может быть использовано для дисперсионного анализа. В активных сейсмических экспериментах, напротив, волны Рэлея рассматриваются в качестве шума и должны быть подавлены с целю выделения именно объемных волн [24-26]. Другим приложением поляризационного анализа является идентификация эффекта расщепления 8-волн, которое несет в себе информацию об анизотропных свойствах среды [27,28]. Применение поляризационных методов анализа позволило значительно повысить разрешающую способность методов скважинной и наземной сейсмики в сложно-построенных средах [29-32]. Во второй главе диссертации приведен обзор некоторых получивших распространение поляризационных методов [33-38].

Уже несколько десятилетий большое внимание уделяется вопросам обнаружения слабых сейсмических событий в связи с необходимостью контроля за соблюдением Договора о всеобъемлющем запрещении ядерных испытаний (ДВЗЯИ). Отработанные для этого методы с использованием частотной фильтрации, поляризационного анализа, поляризационной фильтрации, детального исследования спектрограмм, БК-анализа записей микрогрупп и пр. применяются и для обнаружения записей слабых землетрясений. Все эти процедуры направлены на увеличение отношения сигнал/шум [4,39].

К сожалению, эти методы имеют значительные ограничения, так как они не различает волны, пришедшие в один и тот же момент, но имеющие разные частоты. В связи с этим хорошие результаты дает использование спектрально-временного представления, замечательного, в частности, тем, что оно позволяет интерпретировать и изучать частотно-временные образы сигналов. Спектрально-временной подход получил мощное развитие в работах A.JI. Левшина и А.В. Ландера [13] и по праву стал одним из важнейших методов численного анализа сигналов. Сегодня разработано и активно используется большое количество различных методов спектрально-временного анализа. Среди них можно выделить оконное Фурье-преобразование, преобразование Габора, распределение Вагнера-Вилла (Wigner-Ville transform) [40], S-преобразование [41], дискретное вейвлет-преобразование [42-44], а также непрерывное вейвлет-преобра-зование [45]. Последним двум методам посвящено также большое количество обзоров, например, [46-49]. Во всех этих методах общим является то, что результат анализа некоего сигнала должен содержать в себе не только перечисление его характерных частот, но и сведения об определенных временных координатах, в которых эти частоты проявляют себя.

Вейвлет-преобразование, S-преобразование и преобразование Габора являются результатом свертки исследуемого сигнала с семейством хорошо локализованных как во времени, так и по частоте базисных функций. Поэтому частотно-временное разрешение данных методов напрямую связано с частотно-временной формой используемых функций. Распределение Вагнера-Вилла [40] базируется на вычислении корреляции сигнала с самим собой, но сдвинутым как во времени, так и по частоте, поэтому в нем нет ограничений на частотно-временное разрешение. Однако подобное распределение может содержать области интерференции, не имеющие физического смысла. Для минимизации эффекта подобной интерференции можно использовать сглаженное псевдораспределение Вагнера-Вилла [50]. К сожалению, использование операции сглаживания в данном случае приводит к частотно-временному размыванию спектра и, тем самым, снижает достоинства метода по сравнению с вейвлет-преобразованием. В общем случае, можно улучшить частотно-временное разрешение спектрально-временных методов с использованием подхода, предложенного в [51,52].

В [13] можно найти, пожалуй, первую попытку объединить поляризационный и спектрально-временной анализ. С тех пор было опубликовано большое количество работ, посвященных подобной интеграции. В [53,54] используется метод сингулярного разложения матрицы (SVD) совместно с мульти-вейвлетным анализом на базе дискретных ортогональных вейвлетов, а в [55,56] - на базе непрерывного вейвлет-преобразования. В [57,58] развивается метод, аналогичный методу Морозова [36], но с использованием S-преобразования. Комбинация комплексного ковариационного метода и непрерывного вейвлет-преобразования с использованием комплексных ортогональных вейвлетов Майера (Mayer wavelet) приведена в [59]. И, наконец, оценки некоторых поляризационных параметров с использованием дискретного вейвлет-пакетного разложения предложены в работе [22].

В качестве первого важного результата диссертационной работы, во второй главе предложены три новых метода определения спектрально-временных поляризационных параметров двухкомпонентных [60, 61] и трехкомпонентных [62-65] сигналов, а также алгоритм фильтрации сейсмограмм с применением этих методов. Данные методы разработаны на базе непрерывного вейвлет-преобразования.

Для объемных импульсных волн главными параметрами являются времена пробега, очень слабо зависящие от частоты, и амплитуды. Дисперсия для таких волн не является определяющим понятием. Поверхностные же волны, распространяющиеся в вертикально-неоднородной среде, обладают дисперсией. За счет дисперсии они несут ценную информацию о структуре и материальных параметрах среды. Например, нахождение фазовой и групповой скоростей этих волн позволяет строить скоростные [66] и структурные [67-70] модели поверхностных слоистых сред. В [71,72] развиваются методы анализа глубинной геологической структуры Земли на основе длиннопериодических волн, генерируемых при землетрясениях. Однако для интерпретации структуры и материальных параметров нужно уметь определять по исходным сейсмограммам дисперсионные и диссипативные кривые, которые в данном случае являются функциями частоты.

В литературе можно найти несколько методов определения дисперсионных кривых. В общем случае, определение фазовой скорости волны возможно уже при анализе записи с одного сейсмоприемника, но при этом нужно знать фазовый профиль волны в источнике [73]. Групповая же скорость может быть определена с использованием только одного сейсмоприемника без априорной информации об источнике [74,75]. Как фазовая, так и групповая скорости могут быть получены на основании инверсии записей с двух и более сейсмоприемников, то есть при анализе пространственно-временных сейсмограмм [1,76].

Так как дисперсия описывается не одной переменной, а функцией частоты, то также и для дисперсионного анализа интересным является спектрально-временное представление. Различные методы спектрально-временного анализа сегодня широко используются для определения дисперсионных характеристик волн [40,77-80]. В большинстве этих методов делается попытка определения дисперсионных характеристик только с использованием модулей частотно-временного спектра, который описывает частотно-временное распределение энергии, но не содержит фазовой информации о колебаниях. Это сильно ограничивает возможности анализа. Во-первых, в этом случае может быть получена только групповая скорость и, во вторых, в рамках такого анализа различные частотные моды одной волны не различимы. Однако для описания, например, волны Рэлея, возможность различать и описывать ее отдельные моды является важной.

С другой стороны, сегодня развиты методы анализа, которые используют фазовую информацию колебаний и позволяют определять волновое число как функцию частоты для всех волновых мод. Данная группа методов получила название / — к анализа и широко используется сегодня для определения фазовых скоростей распространяющихся волн [81]. Можно выделить несколько методов этой группы: / — к метод, предложенный Дж. Бургом (J. P. Burg) [82], / — к метод высокого разрешения, предложенный Дж. Капоном (J. Capon) [83] и так называемый A/I/ltiple SYgnal Classification (MUSIC) метод [84]. Все эти методы требуют большого количества исходных сигналов для стабильного определения фазовой скорости в определенном частотном интервале.

В третьей главе предложен новый дисперсионный оператор, описывающий изменения непрерывного вейвлет-спектра волны при ее распространении в среде с дисперсией и диссипацией [60,85]. Параметрами данного оператора являются частотные функции скоростей и диссипации волны, при этом рассмотрены случаи линейной и нелинейной диссипации. На основе дисперсионного оператора предложен новый корреляционный метод определения фазовой скорости волн с несколькими волновыми модами [86,87], а также новый оптимизационный метод, позволяющий определить функции скоростей и затухания для единичной волны [88]. Данные методы дисперсионного анализа сочетают в себе сильные стороны как спектрально-временного подхода, так и f — к методов. Все они базируются на анализе фазовых различий непрерывных вейвлет-спектров сигналов. Поэтому они позволяют определять с высокой точностью фазовые скорости волн (в том числе и с несколькими волновыми модами), обладая при этом хорошим частотно-временным разрешением.

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

В этом контексте можно считать, что некоторые из таких закономерностей сформулированы во второй и третьей главах. Четвертая же глава диссертации посвящена именно практическим аспектам реализации предложенных методов. К сожалению, найти программный комплекс, который сочетал бы в себе частотно-временной анализ в приложении к задачам поляризационного и дисперсионного анализа, не удалось. Исключение составляют, пожалуй, специализированные отраслевые комплексы сейсморазведки, но они, как правило, являются ноу-хау нефтяных компаний и доступ к ним затруднен. Поэтому в четвертой главе рассматривается разработанный лично автором программный комплекс GWL - Geophysical Wavelet Library. Он является свободно распространяемым программным обеспечением с открытым кодом [89], где собраны воедино все рассмотренные в предыдущих главах методы поляризационного и дисперсионного анализа.

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

Результаты работы. Таким образом, на защиту выносятся следующие наиболее существенные научные результаты, полученные лично автором: а) метод определения спектрально-временных поляризационных параметров двухкомпонентных сигналов; б) два метода расчета спектрально-временных поляризационных параметров трехкомпонентных сигналов для случаев плоско-эллиптической и объемно-эллиптической поляризации; в) дисперсионный оператор, описывающий изменение непрерывного вейвлет-спектра волны при ее распространении в среде с дисперсией и диссипацией. Параметрами данного оператора являются частотные функции скоростей и диссипации волны, при этом рассмотрены случаи линейной и нелинейной диссипации; г) корреляционный метод "частотно-скоростного" анализа для определения фазовой скорости волн с несколькими волновыми модами; д) оптимизационный метод, позволяющий определять функции скоростей и затухания для единичной волны.

Научная новизна. Все предложенные в диссертации методы являются новыми. Кроме этого, данные методы обладают значительными преимуществами по сравнению с традиционными методами поляризационного и дисперсионного анализа: а) в отличие от метода комплексного следа, предложенный метод двух-компонентного поляризационного анализа более полно характеризует поляризационные свойства и является полностью обратимым, что позволяет строить на его основе эффективные алгоритмы фильтрации; б) в отличие от стандартного ковариационного метода, в предложенном адаптивном ковариационном методе анализа трехкомпонентных сигналов отсутствует проблема выбора временного окна; в) предложенный дисперсионный оператор позволил установить новую интерпретацию функций фазовой и групповой скоростей волны при ее распространении в среде с дисперсией; г) в отличие от традиционных методов / — к анализа, предложенный в работе новый корреляционный метод может использоваться не для всего сигнала, а только для определенных частотно-временных областей, соответствующих той или иной волне; д) в отличие от традиционных методов инверсии, предложенный оптимизационный метод работает в частотно-временном пространстве, что позволяет определить с высокой точностью функции фазовой и групповой скоростей, а также функцию затухания единичной волны по имеющейся пространственно-временной сейсмограмме.

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

Теоретическая ценность. Все предложенные методы являются фрагментами единого методологического подхода частотно-временного поляризационного и дисперсионного анализа волновых процессов. Эти методы в целом составляют вклад как в математическую теорию вейвлет-анализа, так и в математическое моделирование и методы количественного анализа физических параметров волнового процесса.

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

Апробация работы. Методы, алгоритмы и программное обеспечение, предложенное в диссертационной работе, используется: а) в Центре анализа данных геомагнетизма и космического магнетизма университета г. Киото (Япония) для поляризационного анализа геомагнитных пульсаций [90]; б) в ИМСС УрО РАН как составная часть вибрационного мониторинга ответственных инженерных сооружений [91-93]; в) в Геофизическом исследовательском центре г. Потсдам (Германия) для определения поляризационных характеристик волн, представленных на сейсмограммах землетрясений [94]; г) в Геофизическом институте университета г. Потсдам при разработке методов автоматической классификации волн, представленных на сейсмограммах землетрясений [95]; д) в учебном процессе в рамках спецкурсов "Дополнительные главы теории упругости" [96] и "Современные проблемы механики" [97], проводимых на базе Научно-образовательного центра "Неравновесные переходы в сплошных средах" (г. Пермь) для студентов 5 курса и магистров 1 курса механико-математического факультета Пермского Государственного Университета.

Результаты диссертационной работы докладывались и обсуждались на российских и международных конференциях "AGU Fall Meeting" (San-Francisco, California, 2003 [98], 2004 [99], 2005 [100], 2006 [101], 2007 [102]); "EGU Joint Assembly" (Nice, France, 2003 [103], 2004 [104]); "Advanced Problems in Mechanics" (г. Санкт-Петербург, 2004 [105], 2005 [106], 2006 [107,108]); "CGU Joint Assembly" (Montreal, Canada, 2004 [109]); "SEG Annual Meeting" (Houston, Texas, 2005 [110]); "Зимняя школа по механике сплошных сред" (г. Пермь, 1999 [111], 2005 [112], 2007 [113]); "EAGE International Conference & Exhibition" (г. Санкт-Петербург, 2006 [114]); "International Conference on Scientific Computing" WORLDCOMP'07 (Las-Vegas, Nevada, 2007 [89]).

Связь исследований с научными программами. Исследования по тематике диссертации проводились при финансовой поддержке Российского фонда фундаментальных исследований (проекты 07-01-96029-р-Урал-а, 03-01-00561), Немецкого фонда академических обменов DAAD (грант А-02-14086), Американского фонда гражданских исследований и развития CRDF (грант молодым ученым №Y2-P-09-04), а также Немецкого исследовательского общества DFG (проект DFG-1114 "Mathematical methods in time series analysis and image processing").

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Кулеш, Михаил Александрович

Заключение

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

Все предложенные методы являются фрагментами единого методологического подхода частотно-временного поляризационного и дисперсионного анализа волновых процессов и обладают рядом преимуществ по сравнению с традиционными алгоритмами анализа поляризационных и дисперсионных свойств волн: а) Использование в предложенных методах вейвлет-преобразования как одного из возможных подходов к спектрально-временному анализу позволяет обрабатывать и интерпретировать сейсмические волновые поля, для которых характерно большое количество волн различного типа с различным частотным составом, параметрами поляризации и протяженными зонами их интерференции. б) Поляризационные методы (2.15) и (2.17) являются полностью обратимым, так как исследуемый сигнал связан со своими поляризационными параметрами аналитически. На основе этого можно вычислять сигнал после манипуляции с его поляризационными параметрами, что позволяет строить эффективные алгоритмы фильтрации. Фильтрация с применением адаптивного ковариационного метода (2.21) осуществляется при помощи манипулирования вейвлет-коэффициентами сигнала. в) Предложенные поляризационные методы обладают большим потенциалом по фильтрации сигналов с использованием тех или иных поляризационных параметров. Различные варианты фильтрационных алгоритмов были подробно протестированы автором ранее в работах [61,62,64]. г) Все предложенные дисперсионные операторы (3.6)-(3.13), а также базирующиеся на их основе методы частотно-скоростного анализа (3.14)-(3.18) и оптимизационный метод (3.20)-(3.24) могут использоваться не для всего вейвлет-спектра, а только для определенных частотно-временных областей, соответствующих той или иной волне. Этим свойством не обладает традиционный подход (3.2), в котором волны с одинаковой частотой неразличимы. д) Дисперсионный оператор (3.9) записан в терминах фазовой и групповой скорости. Несмотря на то, что эти параметры однозначно определяются из функции волнового числа, переход к ним позволил получить новую математическую интерпретацию фазовой и групповой скорости (рис. 3.2). е) С помощью предложенного метода частотно-скоростного анализа (3.14)-(3.18) можно определять частотно-зависимые скорости волн, состоящих из нескольких волновых мод. В отличие от традиционных методов / — к анализа, которые требуют большого количества исходных сигналов для стабильного определения фазовой скорости, новый предложенный метод стабильно работает и в случае небольшого числа исходных сигналов. ж) Оптимизационный метод (3.20)-(3.24) может применяться как для вещественных сигналов, что соответствует случаю однокомпонентной записи, так и в общем случае для комплексных, соответствующих записи двух ортогональных компонент. з) Показана связь дисперсионного оператора с параметрами поляризации (3.13), что позволяет говорить о едином методе поляризационно-дисперсионного анализа в вейвлет-пространстве.

Список литературы диссертационного исследования доктор физико-математических наук Кулеш, Михаил Александрович, 2008 год

1. Aki К., Richards P. G. Quantitative seismology: 2nd Ed. Sausalito, CA: University Science Books, 2002.

2. Соколов К. 77. Геофизические методы разведки. Д.: Недра, 1966. 464 с.

3. Моги К. Предсказание землетрясений: Пер. с англ. М.: Мир, 1988. 382 с.

4. Кедров О. К. Сейсмические методы контроля ядерных испытаний. Институт физики Земли РАН, Институт стратегической стабильности Министерства РФ по атомной энергии М.; Саранск: Тип. "Красный Октябрь", 2005. 430 с.

5. Маловичко А. А., Маловичко Д. А., Голубева И. В., Иванова Ю. В. Природная и техногенная сейсмичность Урала // Физико-технические проблемы разработки полезных ископаемых. 2005, № 1. С. 9-18.

6. Явленский К. И., Явленский А. К. Вибродиагностика и прогнозирование качества механических систем. Д.: Машиностроение, 1983. 239 с.

7. Генкин М. Д., Соколова А. Г. Виброакустическая диагностика машин и механизмов. М.: Машиностроение, 1987. 288 с.

8. Викторов И. А. Звуковые поверхностные волны в твердых телах. М.: Наука, 1981. 287 с.

9. Гринченко В. Т. Гармонические колебания и волны в упругих телах. М.: Наука, 1981. 284 с.

10. Бирюков С. В., Гуляев Ю. В., Крылов В. В., Плесский В. П. Поверхностные акустические волны в неоднородных средах. М.: Наука, 1991.416 с.

11. Kolsky Н. Stress Waves in Solids. Oxford: Clarendon, 1953.

12. Маловичко А. А. Кинематическая интерпретация данных цифровой сейсморазведки в условиях вертикально-неоднородных сред. Свердловск: УрО АН СССР, 1990. 270 с.

13. Левшин А. Л., Яновская Т. Б., Ландер А. В. Поверхностные сейсмические волны в горизонтально-неоднородной земле. М.: Наука, 1986. 278 с.

14. Берзон И. С. Сейсморазведка тонкослоистых сред. М.: Наука, 1976. 224 с.

15. Михайлов Д. Н., Николаевский В. Н. Тектонические волны ротационного типа с излучением сейсмических сигналов // Физика Земли. 2000, № 10. С. 3-10.

16. Проблемы нелинейной сейсмики. Отв. ред. А. В. Николаев, И. Н. -Галкин. М.: НаукаГГ987Г288 с.

17. Shieh C.-F., Herrmann R. В. Ground roll: Rejection using polarization filters // Geophysics. 1990. Vol. 55, №9. P. 1216-1222.

18. Алказ В. Г., Онофраш Н. И., Перелъберг А. И. Поляризационный анализ сейсмических колебаний. Кишинев: Штиинца, 1977. 110 с.

19. Reading А. М., Мао W., Gubbins D. Polarization filtering for automatic picking of seismic data and improved converted phase detection // Geophysical Journal International. 2001. Vol. 147, № 1. P. 227-234.

20. Samson J. C., Olson J. V. Data-adaptive polarization filters for multichannel geophysical data // Geophysics. 1981. Vol. 46, №10. P. 1423-1431.

21. Пасечник И. П. Характеристики сейсмических волн при ядерных взрывах и землетрясениях. М.: Наука, 1970. 192 с.

22. Любушин А. А. Вейвлет-пакетный поляризационный метод для автоматического детектирования вступлений Р и S-волн // Физика земли. 2006, № 4. С. 30-39.

23. Pisarenko V. F., Kushnir A. F., Savin I. V. Statistical adaptive algorithms for estimation of onset moments of seismic phases // Phys. Earth Planet. Interiors. 1987. Vol. 47. P. 4-10.

24. Deighan A. J., Watts D. R. Ground-roll suppression using the wavelet transform // Geophysics. 1997. Vol. 62, №6. P. 1896-1903.

25. S. H. Abdul-Jauwad M. F. Khene Two-dimensional wavelet-based ground roll filtering // 70th Annual International Meeting, SEG, Expanded Abstracts. 2000.

26. Wang Y., Singh S. C. Separation of P- and S-wavefield from wide-angle multicomponent OB С data for a basalt model // Geophysical Prospecting. 2003. Vol. 51. P. 233-245.

27. Shieh C.-F. Estimation of shear-wave splitting time using orthogonal transformation // Geophysics. 1997. Vol. 62, №2. P. 657-661.

28. Li X. L., Crampin S. Complex component analysis of shear wave splitting: theory // Geophysical Journal International. 1991. Vol. 107. P. 597-604.

29. Александров С. И. Поляризационный анализ сейсмических волн. М.: ОИФЗ РАН, 1999. 142 с.

30. Гальперин Е. И. Поляризационный метод сейсмических исследований. М.: Недра, 1977. 277 с.

31. Perelberg А. I., Hornbostel S. С. Applications of seismic polarization analysis // Geophysics. 1994. Vol. 59, № 1. P. 119-130.

32. Обрежа В. H. О поляризации сейсмических волн // Информационный бюллетень конференции "Гальперинские чтения 2001: состояние и перспективы развития метода ВСП"(29-31 октября 2001, г. Москва, ОИФЗ РАН).

33. Vidale J. Е. Complex polarization analysis of particle motion // Bulletin of the Seismological Society of America. 1986. Vol. 76, №5. P. 13931405.

34. Smith B. D., Ward S. H. On the computation of polarization ellipse parameters // Geophysics. 1974. Vol. 39, №6. P. 867-869.

35. Rene R. M., Fitter J. L., Forsyth P. M., Kim K. Y., Murray D. J., Walters J. K., Westerman J. D. Multicomponent seismic studies using complex trace analysis // Geophysics. 1986. Vol. 51, № 6. P. 1235-1251.

36. Morozov I. В., Smithson S. B. Instantaneous polarization attributes and directional filtering // Geophysics. 1996. Vol. 61, №3. P. 872-881.

37. Flinn E. A. Signal analysis using rectilinearity and direction of particle motion // Proceedings of the IEEE. 1965. Vol. 53, № 12. P. 1874-1876.

38. Kanasewich E. R. Time Sequence Analysis in Geophysics, 3rd edn. The University of Alberta Press, Edmonton, Alberta, Canada, 1981. 480 p.

39. Землетрясения и микросейсмичность в задачах современной геодинамики Восточно-Европейской платформы / Под ред. Н. В. Шарова, А. А. Маловичко, Ю. К. Щукина. Кн. 1: Землетрясения. 2007. С. 381.

40. Prosser W. Н., Seale М. D., Smith В. Т. Time-frequency analysis of the dispersion of Lamb modes // The Journal of the Acoustical Society of America. 1999. Vol. 105, №5. P. 2669-2676.

41. Stockwell R. G., Mansinha L., Lowe R. P. Localization of the complex spectrum: the S transform // IEEE Transactions on Signal Processing. 1996. Vol. 44, №4. P. 998-1001.

42. Mallat S. A wavelets tour of signal processing. Academic Press, 1999. 637 p.

43. Чуй К. Введение в вейвлеты / пер. с англ. М.: Мир, 2001. 412 с.

44. Добеши И. Десять лекций по вейвлетам. Ижевск: НИЦ "Регулярная и хаотическая динамика", 2001. 464 с.

45. Holschneider М. Wavelets: an Analysis Tool. Oxford University Press, USA, 1995. 448 p.

46. Астафьева Н. М. Вейвлет-анализ: основы теории и примеры применения // Успехи физ. наук. 1996. Т. 166, № 11. С. 1145-1170.

47. Новиков И. Я., Стечкин С. Б. Основы теории всплесков // Успехи мат. наук. 1998. Т. 53, № 6. С. 53-128.

48. Тоггепсе С., Compo G.P. A practical guide to wavelet analysis // Bulletin of the American Meteorological Society. 1998. Vol. 79, № 1. P. 61-78.

49. Дремин И M., Иванов О. В., Нечитайло В. А. Вейвлеты и их использование // Успехи физ. наук. 2001. Т. 171, № 5. С. 465-561.

50. Andria G., Savino М. Interpolated smoothed pseudo Wigner-Ville distribution for accurate spectrum analysis // IEEE Transactions on Instrumentation and Measurement. 1996. Vol. 45, №4. P. 818-823.

51. Kodera K., Gendrin R., de Villedary C. Analysis of time-varying signals with small ВТ values // Phys. Earth Planet. Interiors. 1976. Vol. 12. P. 142-150.

52. Auger F., Flandrin P. Improving the readability of time-frequency and time-scale representations by the reassignment method // IEEE Transaction on Signal Processing. 1995. Vol. 43, №5. P. 1068-1089.

53. Lilly J. M, Park J. Multiwavelet spectral and polarization analyses of seismic records // Geophysical Journal International. 1995. Vol. 122, №3. p. Ю01-1021.

54. Claassen J. P. Robust bearing estimation for three-component stations // Pure and Applied Geophysics. 2001. Vol. 158, № 1. P. 349-374.

55. Paulssen H., Levshin A. L., Lander A. V, Snieder R. Time- and frequency-dependent polarization analysis: anomalous surface waveobservations in Iberia // Geophysical Journal International. 1990. Vol. 103, №2. P. 483^96.

56. RoueffA., Chanussot J., Mars J. I. Estimation of polarization parameters using time-frequency representations and its application to waves separation // Signal Processing. 2006. Vol. 86, № 12. P. 3714-3731.

57. Pinnegar C. R. Polarization analysis and polarization filtering of three-component signals with the time-frequency S transform // Geophysical Journal International. 2006. Vol. 165, №2. P. 596-606.

58. Schimmel M., Gallart J. The inverse S-transform in filters with time-frequency localization // IEEE Transactions on Signal Processing. 2005. Vol. 53, № 11. P. 4417-4422.

59. Soma N., Niitsuma H., Baria R. Reflection technique in time-frequency domain using multicomponent acoustic emission signals and application to geothermal reservoirs 11 Geophysics. 2002. Vol. 67, №3. P. 928-938.

60. Кулеш M. А., Диалло M. С., Холыинайдер M. Вейвлет-анализ эллиптических, дисперсионных и диссипативных свойств поверхностных волн Рэлея // Акустический журнал. 2005. Т. 51, № 4. С. 500-510.

61. Diallo М. S., Kulesh М., Holschneider М., Scherbaum F., Adler F. Characterization of polarization attributes of seismic waves using continuous wavelet transforms // Geophysics. 2006. Vol. 71, №3. P. V67-V77.

62. Diallo M. S., Kulesh M., Holschneider M., Scherbaum F. Instantaneous polarization attributes in the time-frequency domain and wavefield separation // Geophysical Prospecting. 2005. Vol. 53, №5. P. 723-731.

63. Diallo M. S., Kulesh M., Holschneider M., Kurennaya K., Scherbaum F. Instantaneous polarization attributes based on an adaptive approximate covariance method // Geophysics. 2006. Vol. 71, №5. P. V99-V104.

64. Kulesh M., Diallo M. S., Holschneider M., Kurennaya K., Krüger F., Ohrnberger M, Scherbaum F. Polarization analysis in the wavelet domain based on the adaptive covariance method // Geophysical Journal International. 2007. Vol. 170, №2. P. 667-678.

65. Кулеш M. А. Мгновенные поляризационные свойства трехкомпо-нентных сигналов в базисе вейвлетов // Вестник СамГУ Естественнонаучная серия. 2008. Т. 61, № 2. С. 115-129.

66. Campillo M., Singh S. K, Shapiro N., Pacheco J., Herrman R. B. Crustal structure south of the Mexican Volcanic Belt, based on group velocity dispersion // Geofísica Internacional. 1996. Vol. 35, №4. P. 361-370.

67. Shapiro N. M., GordeevA. V., Domínguez J. Average shear-wave velocity structure of the Kamchatka Peninsula from the dispersion of surface waves // Earth Planets Space. 2000. Vol. 52, №9. P. 573-577.

68. Buttkus B. Spectral Analysis and Filter Theory in Applied Geophysics. Springer-Verlag, 2000.

69. Keilis-Borok V. I. Seismic Surface Waves in Laterally Inhomogeneous Earth. Kluwer Academic Publishers, 1989.

70. Lay T., Wallace T. C. Modern Global Seismology. Academic Press, 1995.

71. Udias A. Principles of Seismology. Cambridge University Press, 1999.

72. Dziewonski A., Block S., Landisman M. A technique for the analysis of transient seismic signals 11 Bulletin of the Seismological Society of America. 1969. Vol. 59, № 1. P. 427-^44.

73. Herrmann R. B. Some aspects of band-pass filtering of surface waves // Bulletin of the Seismological Society of America. 1973. Vol. 63, №2. P. 663-671.

74. Xia J., Miller R. D., Park C. B. Estimation of near-surface shear-wave velocity by inversion of Rayleigh waves // Geophysics. 1999. Vol. 64, №3. P. 691-700.

75. Levshin A. L., Pisarenko V. E, Pogrebinsky G. A. On a frequency-time analysis of oscillations I I Annales Geophysicae. 1972. Vol. 28. P. 211218.

76. Abbate A., Frankel J., Das P. Wavelet transform signal processing for dispersion analysis of ultrasonic signals // IEEE Ultrasonic Symposium. 1995. Vol. 1. P. 751-755.

77. Niethammer M., Jacobs L. J., Qu J., Jarzynski J. Time-frequency representations of Lamb waves // The Journal of the Acoustical Society of America. 2001. Vol. 109, №5. P. 1841-1847.

78. Pedersen H. A., Mars J. I., AmblardP.-O. Improving surface-wave group velocity measurements by energy reassignment // Geophysics. 2003. Vol. 68, №2. P. 677-684.

79. Rost S., Thomas Ch. Array seismology: methods and applications // Reviews of Geophysics. 2002. Vol. 40, №3. P. 2-1 2-27.

80. Burg J. P. Three-dimensional filtering with an array of seismometers // Geophysics. 1964. Vol. 29, №5. P. 693-713.

81. Capon J. High-resolution frequency-wavenumber spectrum analysis // Proceedings of the IEEE. 1969. Vol. 57, №8. P. 1408-1418.

82. Schmidt R., Franks R. Multiple source DF signal processing: An experimental system // IEEE Transactions on Antennas and Propagation. 1986. Vol. 34, №3. P. 281-290.

83. Kulesh M., Holschneider M., Diallo M. S., Xie Q., Scherbaum F. Modeling of wave dispersion using continuous wavelet transforms // Pure and Applied Geophysics. 2005. Vol. 162, №5. P. 843-855.

84. Kulesh M., Holschneider M., Ohrnberger M., Luck E. Modeling of wave dispersion using continuous wavelet transforms II: wavelet based frequency-velocity analysis // Pure and Applied Geophysics. 2008. Vol. 165, №2. P. 255-270.

85. Holschneider M., Diallo M. S., Kulesh M., Ohrnberger M., Luck E., Scherbaum F. Characterization of dispersive surface waves using continuous wavelet transforms // Geophysical Journal International. 2005. Vol. 163, №2. P. 463-478.

86. Kulesh M., Nose M., Holschneider M., Yumoto К Polarization analysis of a Pi2 pulsation using continuous wavelet transform // Earth Planets Space. 2007. Vol. 59, №8. P. 961-970.

87. Кулеш M. А., Матвеенко В. П., Шардаков И. Н. Использование вейвлет-анализа для обработки экспериментальных данных вибродиагностики инженерных сооружений // Проблемы машиностроения и надежности машин. 2003, № 6. С. 100-106.

88. Матвеенко В. П., Шардаков И. Н., Судаков А. И., Кулеш М. А., Цветков Р. В. Способ оценки остаточного ресурса автомобильного моста // Патент РФ №2299410. 2007.

89. Матвеенко В. П., Шардаков И. Н., Судаков А. И., Кулеш М. А., Цветков Р. В. Способ мониторинга автомобильного моста // Патент РФ №2317534. 2008.

90. Pacor F., Bindi D., Luzi L., Parolai S., Marzorati S., Monachesi G. Characteristics of strong ground motion data recorded in the Gubbiosedimentary basin (Central Italy) // Bulletin of Earthquake Engineering. 2007. Vol. 5, № 1. P. 27-43.

91. Riggelsen C., Ohrnberger M., Scherbaum F. Dynamic Bayesian networks for real-time classification of seismic signals // in: Knowledge Discovery in Databases (PKDD). 2007. P. 565-572.

92. Кулеш M. А., Шардаков И. H. Волновая динамика упругих сред: метод, материал к спецкурсу "Дополнительные главы теории упругости". Пермь: Перм. ун-т, 2007. 60 с.

93. Корепанов В. В., Кулеш М. А., Шардаков И. Н. Использование вейвлет-анализа для обработки экспериментальных вибродиагностических данных: метод, материал к спецкурсу "Современные проблемы механики". Пермь: Перм. ун-т, 2007. 64 с.

94. Diallo М. S., Holschneider М., Scherbaum F., Kulesh М. Characterization of dispersive Rayleigh waves using wavelet transform // Eos. Trans. AGU, 84(46), Fall Meet. Suppl., Abstract S22B-0442. 2003.

95. Diallo M. S., Kulesh M., Holschneider M., Scherbaum F. Instantaneous polarization attributes in the time-frequency domain: application to wave field separation // Eos Trans. AGU, 85(47), Fall Meet. Suppl., Abstract S31B-1063. 2004.

96. Kulesh M., Nose M., Holschneider M. Polarization analysis of Pi2 pulsations using continuous wavelet transform // Eos Trans. AGU, 87(52), Fall Meet. Suppl., Abstract SM43D-02. 2006.

97. Kulesh M., Holschneider M. Geophysical wavelet library: applications of the continuous wavelet transform to the polarization and dispersion analysis of signals // Eos Trans. AGU, 88(52), Fall Meet. Suppl., Abstract S43B-1300. 2007.

98. Diallo M. S., Holschneider M., Kulesh M., Scherbaum F., Adler F Characterization of the Rayleigh wave polarization attributes with continuous wavelet transform // Geophysical Research Abstracts. 2003. Vol. 5, №11237.

99. Kulesh M., Holschneider M., Diallo M. S., Kurennaya K. Elliptic properties of elastic surface waves in wavelet domain // Proceedings of the XXXIII International Summer School "Advanced Problems in Mechanics". 2005. P. 361-366.

100. Kulesh M. A., Grekova E. F., Shardakov I. N. Rayleigh waves in the isotropic and linear, reduced Cosserat continuum // Proceedings of the XXXIX International Summer School "Advanced Problems in Mechanics". 2006. P. 290-296.

101. Diallo M. S., Kulesh M., Holschneider M., Kurennaya K, Scherbaum F. Estimating polarization attributes with an adaptive covariance method in the wavelet domain // SEG Technical Program Expanded Abstracts. 2005. P. 1014-1017.

102. КулешM. А., ШардаковИ. H. Применение вейвлет-анализак вибродиагностике ответственных инженерных сооружений // Зимняя школа по механике сплошных сред (двенадцатая). Тезисы докладов. Екатеринбург: УрО РАН. 1999. С. 205.

103. Кулеш М. А., Хольшнайдер М., Диалло М., Куренная К. Эллиптические свойства поверхностных упругих волн в пространстве вейвле-тов // Зимняя школа по механике сплошных сред (четырнадцатая). Тезисы докладов. Екатеринбург: УрО РАН. 2005. С. 100.

104. Бхатнагар П. Нелинейные волны в одномерных диспергирующих системах. М.: Мир, 1983. 134 с.

105. Achenbach J. D. Wave propagation in elastic solids. Amsterdam; London: North-Holland publishing company, 1973. 443 p.

106. Лурье А. И. Теория упругости. M.: Наука, 1970. 939 с.

107. Новацкий В. Теория упругости. М.: Мир, 1975. 872 с.

108. УиземДж. Линейные и нелинейные волны. М.: Мир, 1977. 622 с.

109. Дъелесан. Э., Руайе Д. Упругие волны в твердых телах. Применение для обработки сигналов: пер. с франц. М.: Наука, 1982. 424 с.

110. Дженкинс Г., Ватте Д. Спектральный анализ и его приложения. Вып. 1. М.: Мир, 1971.316 с.

111. Рабинер Л., Гоулд Б. Теория и применение цифровой обработки сигналов. М.: Мир, 1978. 848 с.

112. Блейхуд Р. Быстрые алгоритмы цифровой обработки сигналов. М.: Мир, 1998. 448 с.

113. Корн Г., Корн Т. Справочник по математике (для научных работников и инженеров). СПб.: Лань, 2003. 832 с.

114. Гоголадзе В. Г. Отражение и преломление упругих волн. Общая теория граничных волн Рэлея // Тр. Сейсмол. ин-та АН СССР. 1947. Т. 125. С. 1-43.

115. Cheeke J., David N. Fundamentals and Applications of Ultrasonic Waves. CRC Press LLC, 2002.

116. Берзон И. С., Епинатъева А. М., Парийская Г. Н., Стародубровская С. П. Динамические характеристики сейсмических волн в реальных средах. М.: Изд. АН СССР, 1962. 511 с.

117. Wang Yanghua, Guo Лап. Modified Kolsky model for seismic attenuation and dispersion // Journal of Geophysics and Engineering. 2004, № 1. P. 187-196.

118. Jones T. D. Pore fluids and frequency-dependent wave propagation in rocks // Geophysics. 1986. Vol. 51, № 10. P. 1939-1953.

119. Lu Jian-Fei, Hanyga Andrzej. Numerical modelling method for wave propagation in a linear viscoelastic medium with singular memory // Geophysical Journal International. 2004. Vol. 159, №2. P. 688-702.

120. Винник Л. П. Сейсмические свойства мантийных плюмов // Проблемы глобальной геодинамики. Сб. статей под ред. акад. Д. В. Рундквиста. 2000, № 1. С. 104-110.

121. Nigbor R. L. Six-degree-of-freedom ground-motion measurement // Bulletin of the Seismological Society of America. 1994. Vol. 84, №5. P. 1665-1669.

122. Igel H., Schreiber U., Flaws A., Schuberth В., Velikoseltsev A., Cochard A. Rotational motions induced by the M8.1 Tokachi-oki earthquake, September 25, 2003 // Geophysical research letters. 2005. Vol. 32. P. L08309.

123. Савин Г. H., Лукашов А. А., Лыско E. M., Времеенко С. В., Ага-съев Г. Г. Распространение упругих волн в континууме Коссера со стесненным вращением // Прикладная механика. 1970. Т. 6, № 6. С. 37-41.

124. Cosserat Е., CosseratF. Theorie des corps deformables. Hermann, Paris, 1909.

125. Schwartz L. M., Johnson D. L., Feng S. Vibrational modes in granular materials // Physical Review Letters. 1984. Vol. 52, № 10. P. 831-834.

126. Eringen A. C. Microcontinuum Field Theories. I. Foundation and Solids. Springer-Verlag New York, 1998. 319 p.

127. Лялин А. Е., Пирожков В. А., Степанов Р. Д. О распространении поверхностных волн в среде Коссера // Акустический журнал. 1982. Т. 28, № 6. С. 838-840.

128. Кулеш М. А., Матвеенко В. П., Шардаков И. Н. Построение и анализ аналитического решения для поверхностной волны Рэлея в рамках континуума Коссера // Прикладная механика и техническая физика. 2005. Т. 46, № 4. С. 116-124.

129. Кулеш М. А., Матвеенко В. П., Шардаков И. Н. Дисперсия и поляризация поверхностных волн Рэлея для среды Коссера // Известия РАН, Механика твердого тела. 2007, № 4. С. 100-113.

130. Кулеш М. А., Матвеенко В. П., Шардаков И. Н. О распространении упругих поверхностных волн в среде Коссера // Доклады академии наук. 2005. Т. 405, № 2. С. 196-198.

131. Кулеш М. А., Матвеенко В. П., Шардаков И. Н. О распространении упругих поверхностных волн в среде Коссера // Акустический журнал. 2006. Т. 52, № 2. С. 227-235.

132. Кулеш М. А., Матвеенко В. П., Улитин М. В., Шардаков И. Н. Анализ волнового решения уравнений эластокинетики среды Коссера в случае плоских объемных волн // Прикладная механика и техническая физика. 2008. Т. 49, № 2. С. 196-203.

133. Кулеш М. А., Матвеенко В. П., Шардаков И. Н. Построение аналитического решения волны Лэмба в рамках континуума Коссера // Прикладная механика и техническая физика. 2007. Т. 48, № 1. С. 143-150.

134. Новиков Л. В. Основы вейвлет-анализа сигналов: учеб. пособие. СПб.: МОДУС+, 1999. 152 с.

135. Петухов А. П. Введение в теорию базисов всплесков: учеб. пособие. СПб.: Изд-во СПб ГТУ, 1999. 132 с.

136. Воробьев В. И., Грибунин В. Г. Теория и практика вейвлет-преобразования. СПб.: ВУС, 1999. 208 с.

137. Захаров В. Г. Вейвлет-анализ: теория и приложения: учеб. пособие. Ч. 1. Непрерывное вейвлет-преобразование. Пермь: Изд-во ПТУ, 2003. 100 с.

138. Taner M. Т., Koehler F., Sheriff R. E. Complex seismic trace analysis // Geophysics. 1979. Vol. 44, №6. P. 1041-1063.

139. Бендат Дж., ПирсолА. Прикладной анализ случайных данных. М.: Мир, 1989. 540 с.

140. Schimmel М., GallartJ. The use of instantaneous polarization attributes for seismic signal detection and image enhancement // Geophysical Journal International. 2003. Vol. 155, №2. P. 653-668.

141. Montalbetti J. R., Kanasewich E. R. Enhancement of teleseismic body phases with a polarization filter // Geophysical Journal of the Royal Astronomical Society. 1970. Vol. 21. P. 119-129.

142. Esmersoy C. Polarization analysis, rotation and velocity estimation in three-component VSP // in M. N. Toksoz, and R. R. Stewart, eds. Vertical seismic profiling, part B: Advanced concepts: Geophysical Press. 1984. P. 236-255.

143. Jurkevics A. Polarization analysis of three-component array data I I Bulletin of the Seismological Society of America. 1988. Vol. 78, №5. P. 1725-1743.

144. Jackson G. M, Mason I. M., Greenhalgh S. A. Principal component transforms of triaxial recordings by singular value decomposition // Geophysics. 1991. Vol. 56, №4. P. 528-533.

145. Bai C.-Y., Kennett B. L. N. Phase identification and attribute analysis of broadband seismograms at far-regional distances // Journal of Seismology. 2001. Vol. 5, №2. P. 217-231.

146. Nakamura Y. A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface // Quarterly Report of Railway Technical Research Institute. 1989. Vol. 30, № 1. P. 25-33.

147. Fah D., Kind F, GiardiniD. A theoretical investigation of average H/V ratios // Geophysical Journal International. 2001. Vol. 145, №2. P. 535549.

148. Press W. H., Teukolsky S. A., Vetterling W. Т., Flannery B. P. Numerical Recipe in С: The Art of Scientific Computing, 2nd edn. Cambridge University Press, 1992. 1020 p.

149. Lachet C., Bard P. Y. Numerical and theoretical investigations on the possibilities and limitations of Nakamura's technique // Journal of Physics of the Earth. 1994. Vol. 42. P. 377-397.

150. Scherbaum F., Hinzen K.-G., Ohrnberger M. Determination of shallow shear wave velocity profiles in the Cologne, Germany area using ambient vibrations // Geophysical Journal International. 2003. Vol. 152, №3. P. 597-612.

151. Park J., Vernon F. L., Lindberg С. R. Frequency dependent polarization analysis of high-frequency seismograms // Journal of Geophysical Research. 1987. Vol. 92, №B12. P. 12664-12674.

152. Кулеш M. А. Использование вейвлет-анализа для определения поляризационных свойств и фильтрации двухкомпонентных сигналов. Пермь: ИМСС УрО РАН Деп. в ВИНИТИ 14.03.08 №223-В2008, 2008. 29 с.

153. Herrmann R. В. Computer programs in seismology // Saint Louis University, Version 3.0. 1996.

154. Chassande-Mottin E. Méthode de réallocation dans le plan temps-fréquence pour l'analyse et le traitement des signaux non stationnaires // Ph.D. Thesis, University of Cergy-Pontoise France. 1998.

155. Borum S., Jensen K. Additive analysis/synthesis using analytically derived windows // Proc. DAFx-99, Trondheim, Norway, December 1999, Norwegian University of Science and Technology (NTNU) and COST. 1999. P. 125-128.

156. Кулеш M. А. Расчет мгновенных поляризационных свойств трех-компонентных сигналов на основе вейвлет-преобразования. Пермь: ИМСС УрО РАН Деп. в ВИНИТИ 14.03.08 №224-В2008, 2008. 25 с.

157. Cliet С., Dubesset M. Three-component recordings: Interest for land seismic source study // Geophysics. 1987. Vol. 52, №8. P. 1048-1059.

158. Benhama A., Cliet C., Dubesset M. Study and applications of spatial directional filtering in three-component recordings // Geophysical Prospecting. 1988. Vol. 36, №6. P. 591-613.

159. Gherasim M., Hoelting C., Marfurt K. Fort Worth Basin 2-D elastic synthetic depth model // SEG Technical Program Expanded Abstracts. 2004. Vol. 23, № 1. P. 1893-1896.

160. Kjartansson E. Constant Q-wave propagation and attenuation // Journal of Geophysical Research. 1979. Vol. 84, №B9. P. 4737-4748.

161. Turner G., Siggins A. F. Constant Q attenuation of subsurface radar pulses // Geophysics. 1994. Vol. 59, №8. P. 1192-1200.

162. Кулеш М. А. Моделирование распространения волн в среде с дисперсией и диссипацией с использованием вейвлет-преобразования. Пермь: ИМСС УрО РАН Деп. в ВИНИТИ 14.03.08 №225-В2008, 2008. 29 с.

163. Song S., Innanen К. Multiresolution modeling and seismic wavefield reconstruction in attenuating media // Geophysics. 2002. Vol. 67, №4. P. 1192-1201.

164. Кулеш M. А. Определение дисперсионных параметров волн с применением непрерывного вейвлет-преобразования. Пермь: ИМСС УрО РАН Деп. в ВИНИТИ 14.03.08 №226-В2008, 2008. 39 с.

165. Maraun D., Kurths J., Holschneider М. Nonstationary Gaussian processes in wavelet domain: Synthesis, estimation, and significance testing // Physical Review E (Statistical, Nonlinear, and Soft Matter Physics). 2007. Vol. 75, № 1. P. 016707.

166. Gurley K, Kijewski Т., Kareem A. First- and higher-order correlation detection using wavelet transforms // Journal of Engineering Mechanics. 2003. Vol. 129, №2. P. 188-201.

167. Kirkpatrick S., Gelatt C. D., Vecchi M. P. Optimization by simulated annealing I I Science. 1983. Vol. 220, №4598. P. 671-680.

168. Lai C. G., Rix G. J. Simultaneous inversion of Rayleigh phase velocity and attenuation for near-surface site characterization // School of Civil and Environmental Engineering, Georgia Institute of Technology. 1998, № GIT-CEE/GEO-98-2.

169. Fernández G., Periaswamy S., Sweldens Wim LIFTPACK: A software package for wavelet transforms using lifting // Wavelet Applications in Signal and Image Processing IV, Proc. SPIE 2825. 1996. P. 396-^08.

170. MisitiM., Misiti Y, Oppenheim G., PoggiJ.M. Matlab Wavelet Toolbox (Version 4.0): Tutorial and Reference Guide. The Mathworks, Natick, USA, 2007.

171. MusserD. STL Tutorial & Reference Guide: С++ Programming with the Standard Template Library. Addison-Wesley Pub. Co., Reading, Mass., USA, 1996. 400 p.

172. Frigo M., Johnson S. G. FFTW: an adaptive software architecture for the FFT // Proceedings of the 1998 IEEE International Conference on Acoustics, Speech and Signal Processing. 1998. Vol. 3. P. 1381-1384.

173. Потемкин В. MATLAB 6: Среда проектирования инженерных приложений. М.: Диалог-МИФИ, 2003. 448 с.

174. Nosé М. Automated detection of Pi2 pulsations using wavelet analysis: 2. An application for dayside Pi2 pulsation study // Earth Planets Space. 1999. Vol. 51, №1. P. 23-32.

175. Lester M., Hughes W. J., Singer H. J. Polarization patterns of Pi2 magnetic pulsations and the substorm current wedgee // Journal of Geophysical Research. 1983. Vol. 88. P. 7958-7966.

176. Lester M., Hughes W. J., Singer H. J. Longitudinal structure in Pi2 pulsations and the substorm current wedge // Journal of Geophysical Research. 1984. Vol. 89. P. 5489-5494.

177. Троицкая В. А., Гулъелъми А. В. Геомагнитные пульсации и диагностика магнитосферы // Успехи физ. наук. 1969. Т. 97, № 3. С. 453494.

178. Kosaka К., Iyemori T., Nosé M., Bitterly M., Bitterly J. Local time dependence of the dominant frequency of Pi2 pulsations at mid- and low-latitudes // Earth Planets Space. 2002. Vol. 54, №7. P. 771-781.

179. Han D., Iyemori T., Gao Y., Sano Y., Yang F., Li W., Nosé M. Local time dependence of the frequency of Pi2 waves simultaneously observed at 5 low-latitude stations // Earth Planets Space. 2003. Vol. 55, № 10. P. 601612.

180. Nosé M., Liou К., Sutcliffe P. R. Longitudinal dependence of characteristics of low-latitude Pi2 pulsations observed at Kakioka and Hermanus // Earth Planets Space. 2006. Vol. 58, №6. P. 775-783.

181. Beaty K. S., Schmitt D. R. Repeatability of multimode Rayleigh-wave dispersion studies // Geophysics. 2003. Vol. 68, №3. P. 782-790.

182. DESERT-Group Multinational geoscientific research kicks off in the Middle East I I Eos, Transactions, American Geophysical Union. 2000. Vol. 81, №50. P. 609,616-617.

183. On-line Bulletin International Seismological Centre // http://www.isc.ac.uk, Internatl. Seis. Cent., Thatcham, United Kingdom. 2001.

184. Hadiouche O., Kriiger F., Kind R. Mapping the crust in southeastern Germany using Rayleigh waves in the period range 6-16 s // Geophysical Research Letters. 1991. Vol. 18, №6. P. 1087-1090.

185. Schmidt R. O. A signal subspace approach to multiple emitter location and spectral estimation. // PhD Thesis. Stanford University, CA. 1981.

186. Schmidt R. Multiple emitter location and signal parameter estimation // IEEE Transactions on Antennas and Propagation. 1986. Vol. 34, №3. P. 276-280.

187. Wax M., Kailath T. Detection of signals by information theoretic criteria // IEEE Transactions on Acoustics, Speech, and Signal Processing. 1985. Vol. 33, №2. P. 387-392.

188. Akaike H. Information theory and an extension of the maximum likelihood principle // 2nd. International Symposium on Information Theory. 1973. P. 267-281.

189. Forbriger T. Inversion of shallow-seismic wavefields: I. Wavefield transformation// Geophysical Journal International. 2003. Vol. 153, №3. P. 719-734.

190. Forbriger T. Inversion of shallow-seismic wavefields: II. Inferring subsurface properties from wavefield transforms // Geophysical Journal International. 2003. Vol. 153, №3. P. 735-752.