Изменение направленности высокочастотной геоакустической эмиссии в периоды деформационных возмущений тема автореферата и диссертации по физике, 01.04.06 ВАК РФ

На правах рукописи

ЩЕРБИНА Альберт Олегович

084607160

ИЗМЕНЕНИЕ НАПРАВЛЕННОСТИ ВЫСОКОЧАСТОТНОЙ ГЕОАКУСТИЧЕСКОЙ ЭМИССИИ В ПЕРИОДЫ ДЕФОРМАЦИОННЫХ

ВОЗМУЩЕНИЙ

Специальность - 01.04.06 «Акустика»

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

2 2 ШОЛ ?П;Г

с. Паратунка, Елизовский район, Камчатский край - 2010

004607160

Работа выполнена в Институте космофизических исследований и распространения радиоволн ДВО РАН

Научный руководитель: кандидат технических наук, доцент

Марапулец Юрий Валентинович

Официальные оппоненты: доктор физико-математических наук

Дзюба Владимир Пименович

кандидат технических наук Тагильцев Александр Анатольевич

Ведущая организация: Кафедра гидроакустики и ультразвуковой тех-

ники Дальневосточного государственного технического университета

Защита диссертации состоится 02 июля 2010 года в 14 часов на заседании диссертационного совета Д 005.017.01 при Тихоокеанском океанологическом институте им. В.И. Ильичева ДВО РАН по адресу: 690041, г. Владивосток, ул. Балтийская, 43.

С диссертацией можно ознакомиться в библиотеке Тихоокеанского океанологического института им. В.И. Ильичева ДВО РАН.

Автореферат разослан « » мая 2010 года.

Ученый секретарь

диссертационного совета

доктор технических наук, профессор

Коренбаум Владимир Ильич

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность проблемы

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

В лабораторных условиях акустико-эмиссионный контроль позволяет изучать процессы деформации и разрушения материалов. В геофизике метод акустико-эмиссионного контроля широко применяется при изучении состояния горных пород. Упругие колебания могут наблюдаться в широком диапазоне длин волн (от тысяч километров сейсмических волн в земной коре до нанометров дислокационных подвижек и связанных с ними колебаний в различных средах). В геофизике принято рассматривать акустическую эмиссию в четырех частотных диапазонах: сейсмическом (от долей до 10 Гц), высокочастотном сейсмическом (10-100 Гц), сейсмоаку-стическом (100-1000 Гц) и акустическом (более 1 кГц).

На камчатском геодинамическом полигоне ИКИР ДВО РАН (Камчатский край, Елизовский район, с. Паратунка) проводятся исследования геоакустической эмиссии в частотном диапазоне от 0.1 Гц до 10 кГц с помощью гидрофонов, установленных в естественных и искусственных водоемах. В результате исследований было установлено, что в суточном временном интервале до готовящегося землетрясения, наблюдается повышение интенсивности геоакустической эмиссии. Этот эффект наиболее ярко проявляется в килогерцовом диапазоне частот и определяется деформациями пород в пунктах наблюдений на удалении первых сотен километров от эпицентров готовящихся землетрясений. Механизм возникновения описываемого объясняется следующим образом. В моменты подготовки сильных сейсмических событий под действием деформационных сил на обширной площади Земной коры с радиусом порядка сотен километров вокруг их эпицентров накапливаются напряжения. После превышения ими определенного порога наступает период деформационных возмущений с лавинообразным растрескиванием пород и проскальзыванием пластов по существующим разломам, что сопровождается характерным повышением уровня геоакустической эмиссии, проявляющимся как в увеличении амплитуды отдельных импульсов, так и в существенном увеличении их потока (числа импульсов в единицу времени).

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

В ранних работах, проводимых в ИКИР ДВО РАН, предпринимались попытки оценить направленность геоакустического излучения, однако конструктивные особенности использованных приемников не позволили в полной мере исследовать эти свойства. В дальнейшем система геоакустических наблюдений была модернизирована. Вместо пяти направленных гидрофонов была установлена комбинированная приемная система, включающая гидрофон с круговой характеристикой направленности и трехкомпонентный векторный приемник, позволяющий регистрировать в той же точке пространства три компоненты градиента звукового давления. Используя эти данные можно достаточно легко определять пеленг акустических сигналов.

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

Предмет исследования - характеристики высокочастотной геоакустической эмиссии осадочных пород.

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

Для достижения поставленной цели решены следующие задачи:

1. Разработка метода автоматизированного обнаружения геоакустических импульсов и анализа их направленности.

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

3. Проведение регулярных наблюдений, анализ полученных данных.

4. Исследование распределения интенсивности геоакустического излучения по направлениям в различные периоды сейсмической активности.

Научная новизна работы

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

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

Практическая ценность работы

Работа выполнена в соответствии с планами научных исследований ИКИР ДВО РАН, в рамках Программы фундаментальных исследований Президиума РАН № 16, проектов Президиума ДВО РАН №05-1-0-02-051, №06-1-П16-070.

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

Оригинальность программных продуктов подтверждена свидетельствами об отраслевой регистрации разработки № 9537 от 27.11.2007 и № 9766 от 10.01.2008.

Апробация работы

Результаты по теме диссертационной работы докладывались и обсуждались на:

- III международной конференции "Солнечно-земные связи и электромагнитные предвестники землетрясений» в 2004 г., с.Паратунка, Камчатский край;

- международных конференциях по мягким вычислениям и измерениям (SCM'2005, SCM'2006 и SCM'2007) в 2005 - 2007 гг., г. Санкт-Петербург;

- ежегодных научно-технических конференциях профессорско-преподавательского состава и аспирантов КамчатГТУ в 2006, 2007 гг., г.Петропавловск-Камчатский;

- ежегодной научной конференции МГУ «Ломоносовские чтения 2007» в 2007 г., г. Москва;

- международной конференции XXIV General Assembly Of The International Union of Geodesy and Geophysics (IUGG'2007), Perugia, Italia в 2007 г. (два доклада);

- IV международной конференции «Солнечно-земные связи и предвестники землетрясений» в 2007 г., с. Паратунка, Камчатский край;

- III всероссийской научной конференции «Проектирование инженерных и научных приложений в среде Matlab» в 2007 г., г. Санкт-Петербург;

- XX сессии Российского акустического общества в 2008 г., г. Москва;

- VI Всероссийском симпозиуме «Физика геосфер» в 2009 г., г.Владивосток.

Личный вклад

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

Структура работы

Диссертационная работа состоит из введения, трех глав, заключения и списка литературы. В работе содержится 97 листов машинописного текста, 49 рисунков, одна таблица. Список литературы содержит 45 наименований.

Основные положения, выносимые на защиту

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

2. Результаты статистических исследований распределения интенсивности геоакустического излучения по направлениям в сейсмически спокойные периоды и на заключительной стадии подготовки землетрясений в период 2004 -2009 гг.

3. Наличие ярко выраженной анизотропии в распределении интенсивности геоакустического излучения по направлениям в суточном интервале перед землетрясениями.

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

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

Введение

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

Сформулированы основные положения, выносимые на защиту.

Первая глава «Исследование геоакустической эмиссии векторно-фазовыми методами»

В разделе 1.1. «Особенности геоакустической эмиссии» на основе анализа публикаций, посвященных исследованию геоакустической эмиссии, показано, что увеличение интенсивности эмиссии в осадочных породах является следствием роста деформационных процессов. Поскольку прочность пород по отношению к касательным напряжениям меньше, чем к сжатию, то в области излучения геоакустической эмиссии преобладают сдвиговые источники акустического сигнала. Их ориентация обусловлена направлением максимальных касательных напряжений, ориентированных преимущественно под углом 45° к оси наибольшего сжатия. При этом максимумы продольных акустических колебаний приходятся на направления, перпендикулярные этой оси для волны сжатия и совпадающие с ней для волны разряжения (рис. 1,а). Как правило, волна сжатия по интенсивности больше волны разряжения, что связано с особенностями процессов разрушения и диссипации энергии. Диаграмма направленности поперечных колебаний повернута на 45° относительно диаграммы направленности продольных волн (рис. 1,6). Установка приемника в воде, где поперечные колебания сильно затухают, существенно облегчает задачу исследо-

а б

Рис. 1. Диаграммы направленности продольных (а) и поперечных (б) акустических колебаний сдвиговых источников. Стрелками отмечена ось наибольшего сжатия.

^ вания направленности эмиссии, так как

унъ \ \ ту, / позволяет ограничиться исследованием

только продольных волн.

Если допустить, что сдвиговые источники распределены вокруг приемника достаточно равномерно, то их ориентация при отсутствии деформационных возмущений будет произвольной. Поэтому диаграмма распределения интенсивности геоакустического излучения (потока импульсов) по направлениям в точке приема будет соответствовать примеру, приведенному на рис.2,а.

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

=*> +

©

а б

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

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

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

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

В разделе 1.2. «Выбор средств регистрации геоакустической эмиссии» обоснован выбор комбинированного приемника производства ЗАО «Геоакустика» при ФГУП ВНИИФТРИ для изучения направленности геоакустического сигнала. Данный прибор позволяет одновременно регистрировать звуковое давление и три взаимно ортогональные компоненты его градиента. Датчик рассчитан на работу в частотном диапазоне от 5 до 11000 Гц, чувствительность с учетом предварительных усилителей всенаправленного канала давления ~500 мВ-Па"1, чувствительность векторных каналов увеличивается на 6 дБ на октаву и составляет ~20 мВ Па"1 на частоте 10 кГц. Кратко рассмотрены векторно-фазовые методы определения направления прихода акустических волн.

В разделе 1.3. «Особенности применения комбинированного приемника для регистрации геоакустической эмиссии» приведена схема и условия проведения эксперимента по изучению направленных свойств геоакустического сигнала в различные периоды времени (рис. 3). Определены масштабы зоны генерации геоакустического сигнала, оценены размеры источников эмиссии. Расчеты показали, что длина регистрируемых источников колеблется в пределах 0.05-1.5 м, а их максимальное удаление - порядка 30 м. Максимальный угол г? между осью приемника и направлением на источник, при котором все еще возможна регистрация сигнала, определяется из уравнения пассивной гидролокации и равен 84° (рис. 3). В связи со

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

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

Раздел 1.4. «Анализ геоакустического сигнала» посвящен описанию геоакустического импульса. Его форма состоит из двух частей: резкое увеличение интенсивности (фронт) и достаточно продолжительная область затухания (рис. А,а). Фронт и начальный фрагмент спада геоакустического сигнала, с еще достаточно большой амплитудой, превышающей фон до 30 раз, формируют основную часть импульса, несущую важную информацию о направлении на его источник. Остальная часть периода спада в значительной степени формируется реверберацией.

а б

Рис. 4. Амплитудно-временная форма геоакустического импульса и фазовый портрет градиента давления на различных временных участках 1,2,3 (а). Параметры эллипса, описывающего область концентрации отсчетов геоакустического импульса, отображенных в полярных координатах (б). ур ур - проекции градиента давления на ось X и У, соответственно.

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

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

В разделе 1.5. «Математическое описание геоакустического импульса» представлено математическое описание геоакустических импульсов, состоящее из двух частей: описание огибающей импульсов и описание их формы в полярных координатах.

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

Значения степенных коэффициентов экспоненциальных функций получены экспериментальным путем при усреднении параметров порядка 1000 импульсов.

Представить импульс в полярных координатах можно в виде описывающего эллипса, ограничивающего пространство, заполненное отсчетами, распределенными по некоторому закону. Без учета поворота форма описывающего эллипса хорошо описывается уравнением:

Здесь Б{{) - произвольный широкополосный сигнал,

О)

- огибающая сигнала.

х21ягб +у2 т2* =1,

где R6,RV ~ соответственно большая и малая полуось описывающего эллипса.

Переписав данное уравнение из канонического в параметрический вид и добавив в него азимут источника <р получим полное уравнение описывающего эллипса:

(х = R6 cos(or) cos(^) + R „ sin(ar) sin($5) ; (2)

[>> = -R¡¡ cos(or)sin(^i) + ^l( sin(a)cos((3)

где x, y - взаимно ортогональные направления, а - вспомогательный параметр, изменяющийся в диапазоне [0, 2л].

Во второй главе «Особенности системы регистрации геоакустических сигналов» рассматриваются принципы построения аппаратной части и программного обеспечения измерительного комплекса.

В разделе 2.1 «Автоматизированный метод выявления и анализа направленности геоакустического сигнала» описывается разработанный метод автоматизированного обнаружения и анализа направленности сигнала геоакустической эмиссии, основанный на математическом описании импульса, предложенном в гл. 1.

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

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

Принципиально, разработанный метод включает два основных этапа обработки, выполняющихся по очереди в одном окне данных, размером 200 мс. После успешного выполнения этих этапов производится устранение неоднозначности направления прихода энергии сигнала и итоговая запись этого направления. Затем окно перемещается на 200 мс в случае успешного прохождения обоих этапов или 1/5 размера окна в остальных случаях.

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

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

Для этого по максимально удаленной от центра координат точке Б определяется примерный угол поворота основной оси импульса 0 и перпендикуляры к нему Д =<9 + 90° и Д =#-90° (рис. 5). Последние необходимы для определения значения радиуса я окружности, вписанной в эллипс:

Л = 1/2(г(Д) + Г(/?2)), (4)

где г(Р) - среднее удаление точек, расположенных в секторах \fi-Aip, от центра

координат, Д<р - интервал усреднения.

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

Далее происходит определение радиус-вектора центра масс М отсчетов удовлетворяющих условию 17 |> :

М = ^ т,г, / ^/и, ' / (

где г( - радиус-вектор /-го отсчета в окне данных, т, //Г его весовой коэффициент. Если угол между радиус-вектором г. и направлением в больше 90°, то он заменяется на противоположный. Истинная ось импульса совпадает с найденным радиус-вектором центра масс м.

Рис. 5. Определение направления на источник сигнала. VРх, ¥Ру - проекции

градиента давления на ось X и У, соответственно.

После определения оси импульса производится проверка его формы в полярных координатах. Для этого заново, но уже с учетом определенного выше направления оси, вычисляются значения большой re и малой R полуосей эллипса. После этого, по уточненным значениям полуосей, вычисляется эллиптичность описывающего эллипса. На основании проведенных статистических исследований, пороговое значение эллиптичности выбрано равным 0.25. Если полученное значение превышает пороговый уровень, то считается, что форма импульса существенно искажена, и он не подлежит дальнейшей обработке.

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

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

После разработки данный метод успешно прошел апробацию на синтезированном и реальном сигнале. Погрешность измерения направления составила 0.4°.

В разделе 2.2. «Общее описание комплекса» представлена структура спроектированного измерительного аппаратно-программного комплекса (рис. 6), работающего в составе системы исследования геоакустической эмиссии ИКИР.

Время, Измерительный пункт Система телеметрии ИКИР

GPS о. Микижа с. Паратунка

Рис. 6. Система исследования геоакустической эмиссии на озере Микижа ИКИР ДВО РАН. Звездочками обозначены элементы измерительного аппаратно-программного комплекса изучения направленности геоакустического сигнала; КП - комбинированный приемник.

Разработанный комплекс состоит из двух разнесенных узлов: узел регистрации, расположенный непосредственно около места установки акустических приемников,

работающий в автономном режиме, и узел обработки зарегистрированных сигналов, расположенный в лаборатории ИКИР. Узлы объединены между собой системой телеметрии, позволяющей удаленно осуществлять непрерывный контроль состояния геоакустической эмиссии на месте регистрации сигнала.

Кроме этого приведены некоторые особенности оцифровки, хранения, каталогизации и обработки данных.

Система регистрации (рис. 7), описанная в разделе 2.3. «Система регистрации», расположенная на измерительном пункте ИКИР озере Микижа, осуществляет оцифровку и каталогизацию акустического сигнала, регистрируемого вблизи дна водоема на глубине примерно 4 м. Кроме того, она производит вычисление и отображение усредненного акустического сигнала, позволяющего оперативно контролировать его интенсивность. Запись оцифрованного сигнала производится в файлы формата WAV, каждый из которых маркируется с помощью стандартных информационных полей (тегов) информацией о месте и точном времени его регистрации.

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

Для организации и управления процессом оцифровки и записи цифровых данных на диск разработано специальное программное обеспечение «Sound Processor» (свидетельство об отраслевой регистрации разработки № 9537 от 27.11.2007). Данная программа разработана на языке программирования С++ в среде разработки Microsoft Visual Studio с использованием методов объектно-ориентированного программирования, выбор которых, основан на анализе задач, решаемых программой.

Разработанный метод автоматизированного обнаружения и оценки направленности геоакустического сигнала реализован в специальном программном обеспечении «SmartPeleng» (свидетельство об отраслевой регистрации разработки № 9766 от 10.01.2008 г.), функционирующем в лабораторных условиях (рис. 8). Его описание дано в разделе 2.4. «Система обработки». Программный продукт написан на языке программирования С++, среда разработки - Microsoft Visual Studio. Вследствие отсутствия необходимости повторного использования написанного кода, в данной

ЭВМ

Рис. 7. Структурная схема программы регистрации.

программе решено отказаться от использования объектно-ориентированного программирования в пользу модульной структуры.

Архитектура программы построена в соответствии с шаблоном проектирования «Модель-представление-контроллер» (англ. Мос1е1-у1е\у-сопПх>11ег, МУС), которая отличается тем, что программа разделена на три отдельных компонента:

- модель предоставляет данные (обычно для представления), а также реагирует на запросы (обычно от контроллера), изменяя свое состояние,

- представление отвечает за отображение информации (пользовательский интерфейс),

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

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

Совокупность конфигураторов, включенных в интерфейс пользователя, и их настроек позволяет организовать три режима обработки данных:

- Ручная обработка данных. Данный режим позволяет пользователю детально изучать сигнал. Ни одна автоматизированная функция не задействуются.

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

-Пакетная обработка. В данном режиме программа обрабатывает все файлы, найденные в указанной директории. Файлы выбираются и сортируются по заданным маскам.

В третьей главе «Результаты исследования геоакустической эмиссии» производится анализ данных, полученных с помощью разработанного автоматизированного аппаратно-программного комплекса. В разделе 3.1. «Проверка работоспособности комплекса» описан эксперимент, проведенный сразу после установки комплекса на измерительном пункте, позволивший определить его работоспособность и поправку в - 8°, компенсирующую неточную ориентацию приемника по отношению к сторонам света.

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

Рис.8. Блок схема программы обработки геоакустического сигнала.

Разделы 3.3. «Анализ активности геоакустической эмиссии в сейсмически спокойные периоды» и 3.4. «Анализ активности геоакустической эмиссии в периоды возмущений» посвящены анализу распределения потока геоакустических импульсов по направлениям.

Пример сейсмически спокойного (фонового) периода, зарегистрированного 15.11.2004 в условиях отсутствия деформационных и атмосферных возмущений, приведен на рис. 9,а. Среднее значение активности геоакустической эмиссии этого периода составляет 0.19 имп./с. Из иллюстрации видно, что в целом азимутальное распределение геоакустической эмиссии достаточно равномерно. Но с направлений от -20° до 60° и от 210° до 240° наблюдается примерно двукратное повышение ее уровня.

Для сравнения в том же масштабе приведены фоновые периоды зимы 2007-2008 гг. и весны-лета 2008 года (рис. 9,б-в). Активность геоакустической эмиссии первого периода определялась с ноября 2007 г. по февраль 2008 г. и в среднем составляет 0.63 имп./с. Второй период регистрировался с мая по июнь 2008 г., среднее значение его активности равно 0.14 имп./с.

180° 1801 180" а б в

Рис. 9. Азимутальное распределение активности геоакустической эмиссии D(a) (имп./град.-с), соответствующее фоновому периоду весны 2004 г. (а), зимы 2007-2008 гг. (б) и весны-лета 2008 г (в).

Анализ приведенных графиков позволяет сделать некоторые заключения. Во-первых, в фоновые периоды источники акустической эмиссии рассредоточены в некоторой степени равномерно во всем контролируемом объеме. Во-вторых, наблюдается определенная изменчивость направленности активности геоакустической эмиссии длительного характера, что видно из сравнения ее формы в 2004 г. и в 2008 г. При этом две диаграммы 2008 года имеют и некоторые сходства, в частности, в обоих случаях в направлениях около 60°, 125° и 310° заметны локальные максимумы активности. При этом по множеству нерегулярно расположенных лепестков диаграмм направленности эмиссии можно судить о сложном характере напряжений, а по асимметрии расположения максимумов излучения можно сделать вывод о неоднородности свойств среды вокруг точки наблюдений.

Непрерывные наблюдения, проведенные в период с августа 2007 г. по декабрь 2009 г. позволили построить график сезонных изменений активности геоакустической эмиссии Q(t) (рис. 10) с усреднением в суточном временном интервале. Как видно, ее поведение хорошо аппроксимируется синусоидой, максимумы которой приходится на осень, а минимумы - на весну. Можно отметить сезонные и годовые

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

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

Рассмотрим диаграмму распределения интенсивности геоакустического излучения по направлениям данного возмущения в сравнении с фоном того же периода (рис. 11, рис. 9). Из графиков видно, что во всех случаях анизотропия направленности излучения относительно фона значительно усилилась. Причем по некоторым направлениям соотношение интенсивно-стей возмущения и фона может достигать нескольких десятков раз, тогда как по другим они практически равны.

Таблица

03.08 11.08 07.09 Дата, мм.гг

Рис. 10. Сезонные изменения геоакустической активности П(У (сплошная линия) и график синусоиды для сравнения (пунктирная линия).

Период, в течение которого наблюдалось возмущение геоакустической эмиссии, ит Информация о землетрясении (по данным КФ ГС РАН)

Дата и время, Ш" Энергетический класс Азимут, градус Расстояние, км.

10.11.2004 с 02:00 до 12:00 10.11.2004 21:47 9.1 110 14.7

14.12.2007 с 03:00 до 10:00 15.12.2007 09:00 11.6 110 175

14.05.2008 с 00:00 до 08:00 15.05.2008 05:49 11.1 100 127

14.08.2008 с 09:00 до 23:00 15.08.2008 06:03 9.1 157 89

с 18.11.2008 12:00 до 19.11.2008 12:00 19.11.2008 08:25 11.6 87 119

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

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

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

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

Рис. 11. Азимутальное распределение активности геоакустической эмиссии в периоды возмущений

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

Рассмотрим подробнее особенности формирования диаграммы распределения интенсивности геоакустического излучения по направлениям на примере случая, зарегистрированного 10.11.2004 (рис. 11,а) перед землетрясением, произошедшим в тот же день (см. таблицу). На приведенной иллюстрации хорошо заметно наличие двух основных зон излучения, расположенных по направлениям 25±10 и 225±25 градусов. Учитывая, что интенсивность волн разряжения сдвигового источника достаточно мала, считаем, что регистрируемый с указанных направлений сигнал преимущественно состоит из волн сжатия. Таким образом, ось наибольшего сжатия в породах, ориентированная перпендикулярно основным направлениям прихода акустического сигнала, наклонена под углом около 120 градусов.

Диаграммы распределения интенсивности геоакустического излучения по направлениям, приведенные на рис. 11 ,б-в, более сложные - в данные периоды наблюдаются по два хорошо отличимых направления, с которых регистрируется наибольший поток геоакустических импульсов. Наличие данных максимумов свидетельствует о том, что в эти периоды присутствовали два направления сжатия, ориентированные в первом случае под углами 115 и 160 градусов, во втором - 120 и 65 градусов. Во всех рассмотренных случаях концентрация источников геоакустического сигнала выше в северной части контролируемого пространства.

Оценим ориентации осей наибольшего сжатия как направлений, перпендикулярных положению максимума акустического излучения перед пятнадцатью землетрясениями с энергетическим классом больше 10, произошедшими за период 2004 -2009 гг. по азимутам 90 - 150 градусов на расстояниях до 250 км (рис. 12,а). Для этого, учитывая возможные искажения, восстановим каждую диаграмму в соответствии с рис.1 по одному доминирующему лепестку путем добавления остальных, расположив их по направлениям ±90 и 180 градусов от него (рис. 12,6). Полученные фигуры показывают, что ориентация диаграмм направленности геоакустической эмиссии в точке приема для всех возмущений лежит в определенном, хорошо отличимом диапазоне, несмотря на то, что доминирующие лепестки были расположены в разных квадрантах.

Анализ положений лепестков позволяет выделить два возможных направления ориентации осей наибольшего сжатия: 33 и 123 градуса. Такая неоднозначность связана с тем, что наиболее интенсивно сдвиговые источники генерируют акустические колебания в двух направлениях: в направлении оси сжатия и перпендикулярно ей. При этом большее число исследованных угловых диаграмм геоакустической эмиссии имеет максимум по оси 33 градуса, а среднее направление землетрясений, являющихся следствием накопления значительных напряжений около их очагов, равно 131 градусу. Учитывая это и то, что самое интенсивное излучение происходит в направлениях перпендикулярных оси сжатия, наиболее вероятной ориентацией оси наибольшего сжатия является направление с углом 123 градуса относительно севера.

0° о-

б

КОЛ-ВО

случаев

п

\ \ \ - 1

90° 120° 150° 180°

Рис. 12. Угловые характеристики геоакустической эмиссии в точке наблюдения, зарегистрированные

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

Полученный результат соответствует генеральному направлению действия тектонических напряжений у побережья южной Камчатки, которое в среднем равно 125 градусам. На рис. 12,в показана гистограмма распределения углов поворота осей сжатия.

Анализ изменения азимутального распределения активности геоакустической эмиссии во времени показывает, что она может иметь несколько наиболее активных направлений излучения (рис. 13).

14.08

15.08 16.08

День.месяц 11Т

17.08 о 18.11

19.11 20.11

День.месяц 11Т

б

22.11

Рис. 13. Примеры временных диаграмм азимутального распределения активности геоакустической эмиссии Да,0 перед землетрясениями: а) 15.08.2008; б) 19.11.2008. Стрелками отмечены направления от эпицентров землетрясений.

На рис. 13, а показан пример двух одновременно наблюдаемых направлений, которые в данном случае являются, по всей видимости, взаимно ортогональными лепестками искаженной диаграммы сдвигового источника (рис.1). В отдельных случаях направления могут смещаться или проявляться попеременно (рис. 13,6). Характер этих блужданий, по всей видимости, определяется источниками деформационных возмущений и неоднородностями среды вокруг точки наблюдений. Представляет интерес тот факт, что время, затрачиваемое на изменение числа и ориентации осей сжатия, сопровождающееся соответствующими изменениями диаграмм распределения интенсивности геоакустического излучения по направлениям, составляет в данном случае всего 15-30 мин.

В заключении сформулированы основные результаты работы и перечислены их возможные приложения.

ОСНОВНЫЕ ВЫВОДЫ И РЕЗУЛЬТАТЫ РАБОТЫ

Результатом работы, проведенной в рамках диссертации, можно считать следующее:

1. Разработан метод автоматизированного обнаружения геоакустических импульсов, основанный на определении структуры сигнала и реализованный в специально разработанной программе. На базе комбинированного приемника создан из-

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

2. За время исследования с 2004 по 2009 год получен представительный объем статистических данных о распределении интенсивности геоакустического излучения по направлениям. Детально рассмотрено свыше 20 случаев деформационных возмущений, предваряющих сейсмические события с энергетическим классом выше 10. Построен график многолетних сезонных изменений активности геоакустической эмиссии в сейсмически-спокойные периоды. Установлено, что минимумы активности приходятся на весну, максимумы - на осень. При этом активность геоакустической эмиссии варьирует от 0.1 до 1.2 импульсов в секунду.

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

4. По результатам анализа направленности геоакустической эмиссии перед пятнадцатью землетрясениями с энергетическим классом больше 10, произошедшими за период 2004 - 2009 гг. по азимутам 90 - 150 градусов на расстояниях до 250 км, произведена оценка ориентации оси наибольшего сжатия пород в пункте наблюдения. Установлено, что для всех случаев, диапазон ориентации оси колебался в интервале углов 105 - 140 градусов со средним значением в 123 градуса. Полученный результат совпал с генеральным направлением действия тектонических напряжений у побережья Южной Камчатки равным 125 градусам.

Список работ, опубликованных по теме диссертации

По теме диссертации опубликовано 22 работы, из них 6 в рецензируемых журналах (5 из перечня ВАК), 1 монография, 2 свидетельства об отраслевой регистрации разработок, 13 публикаций в материалах научных мероприятий.

Публикации в рецензируемых журналах и монографии

1. Б. М. Шевцов, Ю. В Марапулец., А. О. Щербина. О направленности приповерхностной высокочастотной геоакустической эмиссии в периоды деформационных возмущений // ДАН. - 2010. - Т.430, №1. - С. 119-122.

2. В. А. Гордиенко, Т. В. Гордиенко, Ю. В. Марапулец, А. О. Щербина. Изучение высокочастотной геоакустической эмиссии приемным модулем на базе векторного приемника // Вестник Московского Университета. Серия 3. Физика. Астрономия. - М. - 2009. - №4. - С. 94-98.

3. А. В. Купцов, Ю. В. Марапулец, М. А. Мищенко, О. П. Руленко, Б. М. Шевцов, А. О. Щербина. О связи высокочастотной акустической эмиссии приповерхностных пород с электрическим полем в приземном слое атмосферы // Вулканология и сейсмология. - 2007. - №5. - С.71-76.

4. И. А. Ларионов, А. О. Щербина, М. А. Мищенко. Отклик геоакустической эмиссии на процесс подготовки землетрясений в разных пунктах наблюдений // Вестник. КРАУНЦ, серия науки о Земле. - г. Петропавловск-Камчатский. - 2005. -№2.-С. 108-115.

5. А.О. Щербина, М. А. Мищенко, И. А. Ларионов. Аппаратно-программный комплекс мониторинга геоакустической эмиссии // Вестник. КРАУНЦ, серия науки о Земле. - г. Петропавловск-Камчатский. - 2005. - №2. - С.128-132.

6. Ю. В. Марапулец, А. О. Щербина. Методы исследования пространственной анизотропии геоакустической эмиссии // Электронный журнал «Техническая акустика». - 2008 - № 14. - http://ejta.org.

7. Ю. В. Марапулец, М. А. Мищенко, А. В. Шадрин, А. О. Щербина. Методы исследования высокочастотной геоакустической эмиссии: монография. - г. Петропавловск-Камчатский: КГТУ. - 2008. - 107 с.

Свидетельства об отраслевой регистрации разработок

8. Щербина А. О. Свидетельство об отраслевой регистрации разработки № 9537 от

27.11.2007 «Программа регистрации акустического сигнала Sound Processor v.l.l»

9. Щербина А. О. Свидетельство об отраслевой регистрации разработки № 9766 от

10.01.2008 «Программа для выделения геоакустических импульсов и вычисления их направлений SmartPeleng v0.7»

Публикации материалов научно-технических конференций

10. Ю. В. Марапулец, А. В. Купцов, А. О. Щербина, Г. Ю. Голотенко. Исследование возмущений геоакустической эмиссии в период подготовки землетрясений // Сборник докладов III международной конференции «Солнечно-земные связи и электромагнитные предвестники землетрясений. - г. Петропавловск-Камчатский: ИКИР ДВО РАН. - 2004. - http://ikir.kamchatka.ru/ Russian/Science/2004/l-13.pdf.

11. Ю. В. Марапулец, А. О. Щербина. Аппаратно-программный комплекс для регистрации геоакустических возмущений // Рыбохозяйственное образование и научные исследования в Камчатском регионе. Материалы научно-технической конференции профессорско-преподавательского состава и аспирантов КамчатГТУ. -г. Петропавловск-Камчатский : КамчатГТУ. - 2006. - С. 153-155.

12. Ю. В. Марапулец, А. О. Щербина. Система автоматического поиска и определения направлений геоакустических сигналов на выходе векторно-фазового комбинированного приемника // Теория и практика научных исследований в рыбо-хозяственной отрасли Камчатки. Материалы ежегодной научно-технической конференции профессорско-преподавательского состава КамчатГТУ. - г. Петропавловск-Камчатский: КамчатГТУ. - 2007. - Часть 1.-С. 129-131.

13. В. А. Гордиенко, Т. В. Гордиенко, А. В. Григорьев, А. В. Купцов, Ю. В. Марапулец, А. О. Щербина. Методы изучения высокочастотной геоакустической эмиссии комбинированным приемным модулем на базе векторного приемника // Тез. еже-

годной научной конференции МГУ «Ломоносовские чтения 2007». - Москва : МГУ,-2007.-С.136.

14. Y. Marapulets, A. Kuptsov, I. Larionov, М. Mischenko, A. Sherbina, V. Gordienko. Anisotropy of high-frequency geoacoustic emission at different stages of seismic event preparation // XXIV General Assembly Of The International Union of Geodesy and Geophysics (IUGG'2007). - 2-13 July 2007. - Perugia, Italia.

15.0. Rulenko, Y. Marapulets, M. Mischenko, A. Scherbina. Simultaneous disturbances of high-frequency geoacoustic emission and of electric field in the near-ground air // XXIV General Assembly Of The International Union of Geodesy and Geophysics (IUGG'2007). - 2-13 Juy 2007. - Perugia, Italia.

16. В. А. Гордиенко, Т. В. Гордиенко, Ю. В. Марапулец, А. О. Щербина. Результаты исследования анизотропии геоакустических сигналов на озере Микижа вектор-но-фазовыми методами // Сборник докладов IV Международной конференции «Солнечно-земные связи и предвестники землетрясений». - г. Петропавловск-Камчатский : ЖИР ДВО РАН. - 2007. - С. 140-146.

17. Л. В. Купцов, И. А. Ларионов, Ю. В. Марапулец, А. О. Щербина. Геоакустическая система мониторинга и прогнозирования сейсмической активности на полуострове Камчатка // Сборник докл. Междунар. конф. по мягким вычислениям и измерениям (SCM'2005). - СПб. - 2005. - Т.2. - С. 159-161.

18.А. В. Купцов, И. А. Ларионов, Ю. В. Марапулец, М. А. Мищенко, А. В. Шадрин, А. О. Щербина. Применение нейронных сетей для оценки возмущений геоакустической эмиссии на заключительной стадии подготовки землетрясений на Камчатке // Сборник докл. Междунар. конф. по мягким вычислениям и измерениям (SCM'2006). - СПб. - 2006. - Т.1. - С. 238-240.

19. Я. А. Ларионов, Ю. В. Марапулец, М. А. Мищенко, А. В. Шадрин, А. О. Щербина Нейронные сети на картах Кохонена в задачах оценки возмущений геоакустической эмиссии // Сборник докл. Междунар. конф. по мягким вычислениям и измерениям (SCM'2007). - СПб. - 2007. - Т. 1. - С.290-292.

20. И. А. Ларионов, Ю. В. Марапулец, М. А. Мищенко, А. О. Щербина, А. В. Шадрин. Система реального времени для обработки сигналов геоакустической эмиссии // Сб. докладов III Всероссийской научной конференции «Проектирование инженерных и научных приложений в среде Matlab». - СПб. - 2007. - С. 1435-1442.

21. Я А. Гордиенко, Т. В. Гордиенко, Ю. В. Марапулец, А. О. Щербина. Высокочастотные оперативные предвестники землетрясений // Физические проблемы экологии (Экологическая физика): сб. научн. трудов. - М. : МАКС Пресс. - 2008. -№15.- С. 70-81.

22. Б. М. Шевцов, Ю. В. Марапулец, А. О. Щербина. Особенности анизотропии высокочастотной геоакустической эмиссии на Камчатском полигоне // Физика геосфер: материалы докл. VI Всеросс. симпоз. - г. Владивосток. -2009. - С. 149-154.

Щербина Альберт Олегович

ИЗМЕНЕНИЕ НАПРАВЛЕННОСТИ ВЫСОКОЧАСТОТНОЙ ГЕОАКУСТИЧЕСКОЙ ЭМИССИИ В ПЕРИОДЫ ДЕФОРМАЦИОННЫХ ВОЗМУЩЕНИЙ

АВТОРЕФЕРАТ

Отпечатано в типографии ООО «ФОН»

Подписано к печати: 26.05.2010 г. Формат: 60x84/16 Усл. печ. л. 1,4 Тираж 100 экз. Заказ 1513-Б

ВВЕДЕНИЕ.

Глава 1. ИССЛЕДОВАНИЕ ГЕОАКУСТИЧЕСКОЙ ЭМИССИИ ВЕКТОРНО-ФАЗОВЫМИ МЕТОДАМИ.

1.1. Особенности геоакустической эмиссии.

1.2. Выбор средств регистрации геоакустической эмиссии.

1.3. Особенности применения комбинированного приемника для регистрации геоакустической эмиссии.

1.4. Анализ геоакустического сигнала.

1.5. Математическое описание геоакустического импульса.

Глава 2. ОСОБЕННОСТИ СИСТЕМЫ РЕГИСТРАЦИИ ГЕОАКУСТИЧЕСКИХ СИГНАЛОВ.

2.1. Автоматизированный метод выявления и анализа направленности геоакустического сигнала.

2.2. Общее описание комплекса.

2.3. Система регистрации.

2.4. Система обработки.

Глава 3. РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЯ

ГЕОАКУСТИЧЕСКОЙ ЭМИССИИ.

3.1. Проверка работоспособности комплекса.

3.2. Особенности регистрируемого геоакустического сигнала.

3.3. Анализ активности геоакустической эмиссии в сейсмически спокойные периоды.

3.4. Анализ активности геоакустической эмиссии в периоды возмущений.

Актуальность проблемы

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

В лабораторных условиях акустико-эмиссионный контроль позволяет изучать процессы деформации и разрушения материалов [43]. В геофизике метод акустико-эмиссионного контроля широко применяется при изучении состояния горных пород [36]. Упругие колебания могут наблюдаться в широком диапазоне длин волн (от тысяч километров сейсмических волн в земной коре до нанометров дислокационных подвижек и связанных с ними колебаний в различных средах). В геофизике принято акустическую эмиссию рассматривать в четырех частотных диапазонах: сейсмическом (от долей до 10 Гц), высокочастотном сейсмическом (10-100 Гц), сейсмоакустическом (100-1000 Гц) и акустическом (более 1 кГц).

На камчатском геодинамическом полигоне ИКИР ДВО РАН (Камчатский край, Елизовский район, п. Паратунка) проводятся исследования геоакустической эмиссии в частотном диапазоне от 0.1 Гц до 10 кГц с помощью гидрофонов, установленных в естественных и искусственных водоемах [15, 18]. В результате исследований было установлено, что в суточном временном интервале до готовящегося землетрясения, наблюдается повышение интенсивности геоакустической эмиссии. Этот эффект наиболее ярко проявляется в килогерцовом диапазоне частот и определяется деформациями пород в пунктах наблюдений на удалении первых сотен километров от эпицентров готовящихся землетрясений. Механизм возникновения описываемого объясняется следующим образом. В моменты подготовки сильных сейсмических событий под действием деформационных сил на обширной площади с радиусом порядка сотен километров вокруг его 2 эпицентра в Земной коре накапливаются напряжения. После превышения ими определенного порога наступает период деформационных возмущений с лавинообразным растрескиванием пород и проскальзыванием пластов по существующим разломам, что сопровождается характерным повышением уровня геоакустической эмиссии, проявляющимся как в увеличении амплитуды отдельных импульсов, так и в существенном увеличении их потока (числа импульсов в единицу времени).

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

В ранних работах, проводимых в ИКИР ДВО РАН, уже предпринимались попытки оценить направленность геоакустического излучения, однако конструктивные особенности использованных приемников не позволили в полной мере исследовать эти свойства [15]. В дальнейшем система геоакустических наблюдений была модернизирована. Вместо пяти направленных гидрофонов была установлена комбинированная приемная система, включающая гидрофон с круговой характеристикой направленности и трехкомпонентный векторный приемник, позволяющий регистрировать в той же точке пространства три компоненты градиента звукового давления. Используя эти данные можно достаточно легко определять пеленг акустических сигналов.

Приемная система с такими возможностями, установленная вблизи дна водоема и позволяющая регистрировать сигналы геоакустической эмиссии в частотном диапазоне от 5 до 11000 Гц, использовалась впервые [8 - 10]. Данные наблюдений, полученные в течение пятилетнего периода, позволили исследовать характеристики геоакустической эмиссии как при деформационных возмущениях, так и во время их отсутствия. Среди особенностей эмиссии в периоды подготовки землетрясений отмечаются сильные вариации интенсивности акустических сигналов и изменение их пеленга, исследованию которых и посвящена данная работа.

Предмет исследования - характеристики высокочастотной геоакустической эмиссии осадочных пород.

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

1. Разработка метода автоматизированного обнаружения геоакустических импульсов и анализа их направленности.

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

3. Проведение регулярных наблюдений, анализ полученных данных.

4. Исследование распределения интенсивности геоакустического излучения по направлениям в различные периоды сейсмической активности.

Структура диссертации

Структура диссертации определена исходя из решенных задач и полученных результатов в соответствии с поставленной целью. Диссертационная работа состоит из введения, трех глав и заключения.

ГЛАВА 3. РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЯ ГЕОАКУСТИЧЕСКОЙ ЭМИССИИ

3.1. Проверка работоспособности комплекса

После установки разработанного комплекса на пункте измерений была произведена ее калибровка и оценка работоспособности. Для этого в зимний период, в отсутствии метеорологических и техногенных помех, а так же сильной и средней сейсмотектонической активности, был совершен круговой обход приемника с радиусом 10 метров и остановками каждые 30°. На всем пути осуществлялась генерация сигнала аналогичного геоакустическому. Результаты обработки полученных данных приведены на рис. 3.1. На диаграмме отдельными точками отмечены зарегистрированные импульсы с соответствующими значениями направлений прихода относительно Востока и временем регистрации относительно начала эксперимента.

Время, часы:г.шнуты

Рис. 3.1. Результат проверки работоспособности системы регистрации и анализа направленности геоакустической эмиссии в натурных условиях.

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

74 фиксированного положения к другому. Точки, не принадлежащие этим областям, обусловлены сигналами фоновой активности геоакустической эмиссии и растрескиванием льда.

Анализ полученных данных позволил определить поправку в -8°, компенсирующую неточную ориентацию установленного приемника по сторонам света.

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

3.2. Особенности регистрируемого геоакустического сигнала

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

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

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

75 непосредственно зависит от траектории движения частиц среды. Поэтому, в случае идеальной среды, гистограмма распределения этих коэффициентов будет повторять гистограмму, полученную в результате моделирования, приведенную на рис. 1.10,а, и иметь максимум в нулевой точке. В реальных же условиях этого не происходит, абсолютный максимум гистограммы распределения эллиптичностей описывающего эллипса располагается в области, отличной от нулевой. Иллюстрация этому приведена на рис. 3.2., где представлена информация о сотнях тысяч геоакустических импульсов, зарегистрированных в период с начала 2007 года по текущий момент. Для ее построения использовались импульсы, регистрировавшиеся только в те промежутки времени, в которых отсутствовали помехи и заметные деформационные процессы. В данном случае распределение имеет абсолютный максимум в значении 0.12, а при значении 0.25 ослабляется более чем в 4 раза.

Коэффициент эллиптичности Рис. 3.2. Распределение эллиптичностей орбит движения частиц среды в геоакустических импульсах.

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

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

150-1

§ 120

0 л

1 905

0 60m

1 30

0 -1-1---f--1-1-1

0 50 100 150 200

Размер малой полуоси описывающего эллипса, отн.сд.

0 50 100 150 200

Размер малой полуоси описывающего эллипса, отн.ед. в

Рис. 3.3. Распределение малых полуосей описывающих эллипсов геоакустических импульсов зарегистрированных 22.07.2008 (а), 23.11.2008 (б) и 05.01.2009 (в). Пунктирной линией отмечен уровень шумов соответствующего периода.

В качестве примера влияния ревербераций на результаты работы метода можно рассмотреть импульс, приведенный на рис. 1.17, где видно, что направление достаточно высокого по амплитуде сигнала реверберации не совпадает с основным направлением импульса. Поэтому, при расчете малой полуоси, которая, в данном случае, будет определяться относительно направления 40° (220°), будет учтена часть сигнала реверберации, лежащая в диапазонах 40° ± 20° и 220° ± 20° (см. Главу 2). В приведенном случае уровень данной части сигнала составляет примерно 70 единиц АЦП, что согласуется с приведенными графиками (рис. 3.3).

3.3. Анализ активности геоакустической эмиссии в сейсмически спокойные периоды

Для регистрации сейсмически спокойной (фоновой) активности геоакустической эмиссии необходимо чтобы в пределах суток до и после периода регистрации отсутствовали достаточно сильные землетрясения. Кроме этого обязательно отсутствие помех как техногенного, так и метеорологического характера [15, 26, 45]. Активным или возмущенным назовем период времени, в течение которого наблюдается заметное, относительно фонового периода, повышение уровня геоакустической эмиссии. Часто такие периоды регистрируются в пределах суток до землетрясения.

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

Пример сейсмически спокойного (фонового) периода, зарегистрированного 15.11.2004 в условиях отсутствия деформационных и атмосферных возмущений, приведен на рис. 3.4.

Здесь и далее под усредненным сигналом геоакустической эмиссии Ps подразумевается 6-й диапазон системы фильтрации [26, 42], полученный выделением из геоакустического сигнала полосы 1500-6000 Гц с последующей разбивкой на односекундные интервалы и их усреднением.

0.0200.0160.0130.0100.0060.003

Ps, O(f).

Па имп./с

J" ii. JI I I .,! I -. - . . I.-1 ь.4.» I.4.

3.40 2.83 2.27 1.70 1.13 0.56

0 01 02 03 04 05 06 07 08 09 10 11 12

Время, UT

Рис. 3.4. Графики усредненного сигнала геоакустической эмиссии (сплошная) и активности геоакустической эмиссии (пунктирная), соответствующие фоновому периоду 15.11.2004. Сплошной линей обозначен усредненный сигнал Ps, пунктирной — геоакустическая активность £2(0

Как видно, в данный период отсутствуют сильные возмущения, как усредненного сигнала геоакустической эмиссии Ps, так и ее активности Q(t), среднее значение которой составляет 0.19 имп./с. На графике азимутального распределения активности геоакустической эмиссии, так же не наблюдается сильных выбросов (рис. 3.5) по каким-либо отдельным направлениям. Но при этом отмечается наличие достаточно широких областей повышенной активности, соответствующих азимутам от -20° до 60° и от 210° до 240°.

Для сравнения в том же масштабе приведены фоновые периоды зимы 2007-2008 гг. и весны-лета 2008 года (рис. 3.6). Активность геоакустической эмиссии первого периода определялась с ноября 2007 г. по февраль 2008 г. и в среднем составляет 0.63 имп./с. Второй период регистрировался с мая по июнь 2008 г., среднее значение его активности равно 0.14 имп./с.

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

Рис. 3.5. Азимутальное распределение активности геоакустической эмиссии D(a), соответствующее фоновому периоду 15.11.2004.

По множеству нерегулярно расположенных лепестков на диаграмме распределения интенсивности геоакустического излучения по направлениям (рис. 3.6) можно судить о сложном характере напряжений, а по асимметрии расположения максимумов излучения можно сделать вывод о неоднородности свойств среды вокруг точки наблюдений.

Рис. 3.6. Азимутальное распределение активности геоакустической эмиссии D(а), соответствующие фоновым периодам зимы 2007-2008 гг. (а) и весны-лета 2008 (б).

Анализ приведенных графиков позволяет сделать некоторые заключения. Во-первых, в фоновые периоды источники геоакустической эмиссии рассредоточены в некоторой степени равномерно во всем контролируемом объеме. Во-вторых, наблюдается некоторая изменчивость азимутального распределения активности геоакустической эмиссии длительного характера, что видно из сравнения ее в 2004 г. и в 2008 г. Однако две диаграммы 2008 года имеют и некоторые сходства, в частности, в обоих случаях в направлениях около 60°, 125° и 310° заметны локальные максимумы активности.

Непрерывные наблюдения, проведенные в период с августа 2007 г. по декабрь 2009 г. позволили построить график сезонных изменений активности геоакустической эмиссии Q.(t) (рис. 3.7) с усреднением в суточном временном интервале. Как видно, ее поведение хорошо аппроксимируется синусоидой, максимумы которой приходится на осень, а минимумы — на весну. Можно отметить сезонные и годовые вариации активности, а также заметную повторяемость ее локальных максимумов и минимумов, что, видимо, определяется метеорологическими факторами.

1.8 1.6 1.4 и 1.2 | 1.9

S0'8 G0.6

0.4

0.2 0

07.07 09.07 11.07 01.08 03.08 05.08 07.08 09.08 11.08 01.09 03.09 05.09 07.09 09.09 11.09

Дата, мм гг

Рис. 3.7. Сезонные изменения активности геоакустической эмиссии Q (t) (сплошная линия) и график синусоиды для сравнения (пунктирная линия).

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

Л. v\ / ~ N ч S h 4

V. \ С) / / \ \ / ч \

J \ / /г- \ \ И г' / Л \\ /J \

N / у * /[ V

V \ \ / у

N

3.4. Анализ активности геоакустической эмиссии в периоды возмущений

Как было отмечено выше, периоды возмущений характеризуются значительным повышением уровня геоакустической эмиссии, которые часто наблюдаются в преддверии сейсмических событий [15, 25]. Характеристики рассмотренных землетрясений и соответствующих им возмущений приведены в таблице.

ЗАКЛЮЧЕНИЕ

Результатом работы, проведенной в рамках диссертации, можно считать следующее:

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

2. За время исследования с 2004 по 2009 год получен представительный объем статистических данных о распределении интенсивности геоакустического излучения по направлениям. Детально рассмотрено свыше 20 случаев деформационных возмущений, предваряющих сейсмические события с энергетическим классом выше 10. Построен график многолетних сезонных изменений активности геоакустической эмиссии в сейсмически-спокойные периоды. Установлено, что минимумы активности приходятся на весну, максимумы - на осень. При этом активность геоакустической эмиссии варьирует от 0.1 до 1.2 импульсов в секунду.

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

4. По результатам анализа направленности геоакустической эмиссии перед пятнадцатью землетрясениями с энергетическим классом больше 10, произошедшими за период 2004 — 2009 гг. по азимутам 90 - 150 градусов на расстояниях до 250 км, произведена оценка ориентации оси наибольшего сжатия пород в пункте наблюдения. Установлено, что для всех случаев, диапазон ориентации оси колебался в интервале углов 105 — 140 градусов со средним значением в 123 градуса. Полученный результат совпал с генеральным направлением действия тектонических напряжений у побережья Южной Камчатки равным 125 градусам.

Дальнейшее направление работ предполагает модернизацию методов обработки данных, в том числе с вертикального канала. Запланированное на 2010 год дооснащение системы регистрации дополнительными векторными приемниками сделает возможным детальное изучение механизмов генерации геоакустического сигнала и вычисление тензоров сейсмических моментов источников. Полученные результаты могут быть использованы при создании систем акустического контроля состояния приповерхностных пород.

Кроме этого, будут улучшены методы автоматизации сбора, хранения и статистической обработки непрерывно увеличивающегося объема данных.

1. Аносов Г. И. Глубинное сейсмическое зондирование Камчатки / Г. И. Аносов, С. К. Биккенина, А. А. Попов, К. Ф. Сергеев К. Ф., В. К. Утна-син, В. И. Федорченко. -М. : Наука, 1978. 130 с.

2. Виноградов С. Д. Условия на разрыве и спектры излучаемых им волн / С. Д. Виноградов // Изв. АН СССР. Физика Земли. 1976. - № 7. - С. 20-26.

3. Виноградов С.Д. Физические процессы в очагах землетрясений / С. Д. Виноградов, К. И. Кузнецова, А. Г. Москвина и др. // М. : Наука. -1980.-С. 129-140.

4. Гордиенко В. А. Векторно-фазовые методы в акустике / В. А. Гордиен-ко. М. : ФИЗМАТ-ЛИТ, 2007. - 480 с.

5. Гамма Э. Приемы объектно-ориентированного проектирования. Паттерны проектирования = Design Patterns: Elements of Reusable Object-Oriented Software / Э. Гамма, P. Хелм, P. Джонсон, Дж. Влиссидес. — СПб. : «Питер», 2007.-366 с. ISBN 978-5-469-01136-1.

6. Гордиенко В. А. Векторно-фазовые методы в низкочастотной аэро- и сейсмоакустике / В. А. Гордиенко, Б. И. Гончаренко // Сб. трудов семинара научной школы профессора Красильникова. М. : МГУ. 2002. -С.209-229.

7. Гордиенко В. А. Геоакустическая локация областей подготовки землетрясений / В. А. Гордиенко, Т. В. Гордиенко, А. В. Купцов, И. А. Ларионов, Ю.В. Марапулец, А. Н. Рутенко, Б. М. Шевцов // ДАН. 2006. -Т. 407,№5.-С. 669-672.

8. Долгих Г. И. Деформационные и акустические предвестники землетрясений / Г. И. Долгих, А. В. Купцов, И. А. Ларионов, Ю. В. Марапулец, В. А. Швец, Б. М. Шевцов, О. Н. Широков, В. А. Чупин, С. В. Яковенко // ДАН. 2007. - Т.413, № 5. - С. 96-100.

9. Коплиен Дж. Программирование на С++. Классика CS / Дж. Коплиен. — СПб. : Питер, 2005. 479 с.

10. М.Короновский Н. В. Напряженное состояние Земной коры / Н. В. Коро-новский //Соросовский образовательный журнал. — 1997. — №1. С. 5056.

11. Купцов А. В. Изменение характера геоакустической эмиссии в связи с землетрясением на Камчатке / А. В. Купцов // Физика Земли. 2005. — № Ю.-С. 59-65.

12. Купцов А. В. Особенности геоакустической эмиссии при подготовке Камчатских землетрясений / А. В. Купцов, И. А. Ларионов, Б. М. Шевцов // Вулканология и сейсмология. 2005. - № 5. - С. 45-49.

13. Марапулец Ю. В. Методы исследования пространственной анизотропии геоакустической эмиссии / Ю. В. Марапулец, А. О. Щербина // Электронный журнал "Техническая акустика". 2008. - №14. - http://ejta.org.

14. Марапулец Ю. В. Методы исследования высокочастотной геоакустической эмиссии : Монография / Ю. В. Марапулец, М. А. Мищенко, А. В. Шадрин, А. О. Щербина. г. Петропавловск-Камчатский : КГТУ. -2008. - 107 с.

15. Райс Дж. Механика очага землетрясений / Дж. Райе. — М. : Мир, 1982. -217 с.

16. Рихтер Дж. Создание эффективных \У1п32-приложений / Дж. Рихтер. — СПб. : Питер, 2003. 752 с.

17. Свердлин Г. М. Прикладная гидроакустика: Учебное пособие / Г. М. Свердлин. — 2-е издание, переработанное и дополненное. JI. : Судостроение, 1990. - 320 с.

18. Свидетельство об отраслевой регистрации разработки № 9537 от2711.2007 «Программа регистрации акустического сигнала Sound Processor v. 1.1»

19. Свидетельство об отраслевой регистрации разработки № 9766 от1001.2008 «Программа для выделения геоакустических импульсов и вычисления их направлений SmartPeleng v0.7»

20. Скребнев Г. К. Комбинированные гидроакустические приемники / Г. К. Скребнев. СПб. : Элмор, 1997. - 200 с.

21. Славина Л.Б. Исследование по физике землетрясений / Л.Б. Славина // М. : Наука. 1976. - С. 217-236.

22. Соболев Г. А. Физика землетрясений и предвестников / Г. А. Соболев, А. В. Пономарев. М. : Наука, 2003. - 270 с.

23. Страуструп Б. Дизайн и эволюция С++ / Б. Страуструп; пер. с англ. — М. : ДМК Пресс; СПб. : Питер, 2006. 448 е.

24. Урик Роберт Дж. Основы гидроакустики / Дж. Роберт Урик; пер. с англ. Л. : Судостроение, 1978. - 448 с.

25. Шамина О. Г. Модельные исследования неоднородных и трещиноватых сред / О. Г. Шамина, В. И. Понятовская. М. : ИФЗ РАН, 1993. - 179 с.

26. Шевцов Б. М. Модели напряженно-деформированного состояния горных пород при подготовке землетрясений и их связь с геоакустическими наблюдениями / Б. М. Шевцов, А.С. Пережогин // Вычислительные технологии. 2009. - Т. 14, № 3. - С. 49-57.

27. Шевцов Б. М. О направленности приповерхностной высокочастотной геоакустической эмиссии в периоды деформационных возмущений / Шевцов Б. М., Ю. В. Марапулец, А. О. Щербина // ДАН. 2010. - Т. 430, №1. С. 119-122.

28. Щербина А.О. Аппаратно-программный комплекс мониторинга геоакустической эмиссии / А. О. Щербина, М. А. Мищенко, И. А. Ларионов // Вестник. КРАУНЦ, серия науки о Земле. — г. Петропавловск-Камчатский. 2005. - №2. - С. 128-132.

29. Christian U. G. Acoustic emission testing / U. Grosse Christian, Ohtsu Ma-sayasu. Springer-Verlag, 2010. - 399 p.

30. Gregori, G.P. Lognormality and fractality. A unique rationale for different phenomena / G. P. Gregori, G. Paparo // Phys. Rev. 2009.