Вестник КРАУНЦ. Физ.-мат. науки. 2021. Т. 34. №1. C. 137-149. ISSN 2079-6641
УДК 534.23:537.87:550.34 Научная статья
Аномалии перед камчатскими землетрясениями (ML > 4.75) в сигналах электромагнитного излучения и геоакустической
эмиссии в 2013 году
О. О. Луковенкова, Е. И. Малкин, М. А. Мищенко, А. А. Солодчук
Институт космофизических исследований и распространения радиоволн ДВО РАН, 684034, Камчатский край, с. Паратунка, ул. Мирная, 7 E-mail: micle@ikir.ru
В работе рассматриваются предсейсмические аномалии сигналов электромагнитного излучения в диапазоне очень низких частот (ОНЧ) и геоакустической эмиссии, зарегистрированных на Камчатке в 2013 году. Выявлены признаки возникновения аномальных возмущений в исследуемых сигналах и особенности их проявления. Проведен совместный анализ предсейсмических аномалий в сигналах электромагнитного ОНЧ-излучения и геоакустической эмиссии методом наложения эпох. Представлены оценки вероятности возникновения таких аномалий в рассматриваемых геофизических сигналах.
Ключевые слова: электромагнитное излучение, геоакустическая эмиссия, предсейсмическая аномалия, метод наложения эпох
DOI: 10.26117/2079-6641-2021-34-1-137-149
Поступила в редакцию: 17.12.2020 В окончательном варианте: 05.03.2021
Для цитирования. Луковенкова О. О., Малкин Е. И., Мищенко М.А., Солодчук А. А. Аномалии перед камчатскими землетрясениями (ML > 4.75) в сигналах электромагнитного излучения и геоакустической эмиссии в 2013 году // Вестник КРАУНЦ. Физ.-мат. науки. 2021. Т. 34. № 1. C. 137-149. DOI: 10.26117/2079-6641-2021-34-1-137-149
Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.Org/licenses/by/4.0/deed.ru)
© Луковенкова О. О. и др., 2021
Введение
Результаты ранее проведенных исследований [1, 2, 3, 4] показывают, что перед сильными землетрясениями на полуострове Камчатка наблюдаются аномалии электромагнитного ОНЧ-излучения (ЭИ) и геоакустической эмиссии (ГАЭ). В указанных выше работах аномалии в электромагнитных и геоакустических сигналах рассматривались независимо друг от друга. Авторы настоящей работы полагают, что предсейсмические аномалии электромагнитного ОНЧ-излучения и геоакустической эмиссии могут быть связаны между собой. Данная работа посвящена исследованию предсейсмических аномалий сигналов ЭИ и ГАЭ и проверке гипотезы об их связи.
Финансирование. Исследование выполнялось без финансовой поддержки фондов.
Признаки аномалий
Под аномалией будем понимать существенное изменение параметров принимаемого сигнала. В исследуемых геофизических сигналах аномалии проявляются по-разному.
Аномалии в электромагнитных сигналах. Источником естественного электромагнитного излучения атмосферно-грозового происхождения являются молниевые разряды. Электромагнитный импульс, генерируемый молниевым разрядом (далее атмосферик), представляет собой периодическое колебание с крутым передним фронтом. Для оценки плотности потока электромагнитной энергии атмосферика вычисляется вектор Умова-Пойнтинга 8, который равен векторному произведению
8=Е х Н = ЩН8 - Е8Иу)+](ЕхН8 - Е8Нх)+к(ЕхИу - ЕуИх),
где Е и Н — векторы напряжённости электрического и магнитного полей соответственно; Ех, Еу, Ег, Нх, Ну, Н — проекции векторов Е и Н на координатные оси.
Расчет 8 выполняется для фрагментов сигнала длительностью 0.01-0.1 с, содержащих атмосферики. По полученным значениям 8 составляется временной ряд из среднесекундных значений модуля вектора Умова-Пойнтинга I(г).
При распространении атмосфериков в волноводе Земля — Ионосфера наблюдается их затухание. Наиболее сильное затухание возникает при отражении атмосфериков от дневного ^-слоя ионосферы. Таким образом, появление ^-слоя приводит к появлению суточных периодических колебаний во временном ряде I(г).
Рис. 1. Пример ряда среднесекундных значений модуля вектора Умова-Пойнтинга I (г)
На рис. 1 изображен пример ряда I(г), вычисленного для электромагнитного сигнала, зарегистрированного в условиях хорошей погоды. Серыми точками обозначены рассчитанные значения I(г), а черной линией — ряд I(г), полученный наложением суточных интервалов за год и нормированием на количество суток.
Результаты анализа I(г) показывают, что иногда в ряде возникают нарушения нормальной суточной периодичности. Первым типом аномалий электромагнитного
сигнала являются сигналы с резким изменением периода суточного хода I(t), данные аномалии обозначаются EDD (рис. 2а).
Вторым типом аномалий являются сигналы с нестандартными значениями дисперсии I(t) относительно I(t). Аномалии этого типа — сигналы со значением дисперсии I(t) более 0.8, обозначаются EDS (рис. 2б).
Рис. 2. Примеры аномалий электромагнитного сигнала: а — аномалия типа EDD, б — аномалия типа EDS
Для выявления аномалий в сигналах электромагнитного излучения выделены следующие признаки: 1) нестандартная величина периода суточного хода ряда среднесекундных значений модуля вектора Умова-Пойнтинга I(t) (у обычных сигналов равна 1440 мин); 2) превышение значения дисперсии I(t) (значение дисперсии более 0.8).
Аномалии в геоакустических сигналах. Сигналы геоакустической эмиссии складываются из последовательности релаксационных импульсов различной амплитуды и длительности, с крутым передним фронтом и частотой заполнения от единиц до десятков килогерц [5]. Результаты анализа геоакустических сигналов показывают, что в фоновый период наблюдаются незначительные по амплитуде геоакустические импульсы с частотой следования до единиц импульсов в секунду (рис. 3).
Рис. 3. Пример сигнала ГАЭ, зарегистрированного в сейсмически спокойный период: а — шум, б — импульс
В период подготовки землетрясения во время роста напряжений и скорости деформирования пород наблюдается увеличение как амплитуды импульсов, так и их количества в единицу времени, при этом частота следования импульсов достигает десятков и даже сотен в секунду (рис. 4).
Рис. 4. Пример серии импульсов сигнала ГАЭ, зарегистрированного в период подготовки землетрясения
Аномальные возмущения присутствуют и в обработанных накопленных сигналах ГАЭ. Обработка сигналов ГАЭ включает цифровую фильтрацию в семи частотных поддиапазонах: 0.1-10 Гц, 30-60 Гц, 70-200 Гц, 200-600 Гц, 600-2000 Гц, 20006500 Гц и более 6500 Гц, детектирование, интегрирование и синхронную запись отсчетов с частотой 1 Гц.
В периоды подготовки землетрясений аномалии в накопленных сигналах ГАЭ проявляются либо в виде характерного увеличения среднего уровня сигнала в
несколько раз (обозначаются ЛБ, рис. 5), либо в виде групп квазипериодически появляющихся импульсов небольшой длительности (обозначаются Л/, рис. 6).
Рис. 5. Пример обработанного накопленного сигнала ГАЭ. На графиках а, б, в, г, д, е, ж — акустическое давление Ра, накопленное в частотных поддиапазонах: 0.1-10 Гц, 30-60 Гц, 70-200 Гц, 200-600 Гц, 600-2000 Гц, 20006500 Гц, больше 6500 Гц соответственно. Прямоугольником выделены предсейсмические аномалии сплошного типа ЛБ (графики д, е, ж)
С учетом вышесказанного, были определены следующие признаки аномалий в геоакустических сигналах: 1) по амплитуде — амплитуда накопленного сигнала более чем в четыре раза превышает фоновые значения; 2) по длительности — аномалия длится более 15 минут; 3) по типу — аномалия в накопленном сигнале ГАЭ имеет сплошной (ЛБ) или импульсный характер (Л/).
Рис. 6. Пример обработанного накопленного сигнала ГАЭ. На графиках а, б, в, г, д, е, ж — акустическое давление Ра, накопленное в частотных поддиапазонах: 0.1-10 Гц, 30-60 Гц, 70-200 Гц, 200-600 Гц, 600-2000 Гц, 20006500 Гц, больше 6500 Гц соответственно. Прямоугольником выделены предсейсмические аномалии импульсного типа Л1 (графики е и ж). Стрелкой указано время регистрации серии импульсов, изображенной на рис. 4
Совместный анализ предсейсмических аномалий ЭИ и ГАЭ
Анализ проводился на основе данных наблюдений за 2013 год. Была сформирована выборка аномалий электромагнитных и геоакустических сигналов, зарегистрированных в этот период (всего 794 аномалии: 158 аномалий геоакустических сигналов, зарегистрированных в пункте наблюдений "Микижа",
и 636 аномалий электромагнитных сигналов, зарегистрированный в пункте наблюдений "Карымшина"). Для каждой аномалии приводились ее тип (EDD, EDS, AI, AS), время начала и длительность (с точностью до минут). За период года наблюдений на основе данных об аномалиях были сформированы два формализованных по времени ряда с разрешением 1 мин (далее ряды аномалий). Каждому элементу этих рядов присваивалось значение 1 (наличие аномалии) или 0 (отсутствие аномалии). Пример построения ряда аномалий приведен на рис. 7.
а)
Дата, UT Тип Начало аномалии, UT Длительность, мин
02.08.2013 AI 1:12 36
02.08.2013 AI 15:40 41
03.08.2013 AI 0:13 16
03.08.2013 AI 1:45 33
Значение ряда 1h
оН
02.08.2013 00:00
02.08.2013 12:00
б)
03.08.2013 00:00
Дата и время, UT
Рис. 7. Пример построения ряда аномалий: а — выборка, б — формализованный ряд аномалий
В исследовании использовались данные о землетрясениях из регионального сейсмического каталога (Камчатский филиал ФИЦ ЕГС РАН, Каталог землетрясений Камчатки и Командорских островов (1962 г. - наст. вр.), http://sdis.emsd.ru/info/earthquakes/catalogue.php). Выбирались землетрясения с локальной магнитудой Мь не менее 4.75 (всего 200 землетрясений). Землетрясения были разделены на четыре группы: Q\, й2, йэ, й4. В группу Q\ были включены землетрясения, для которых выполнялось условие Яй/Яр < 1, где Яр = ю04Э*Мь — радиус зоны Добровольского [6], Яй — расстояние от пункта наблюдений "Микижа" до эпицентра землетрясения; в группу й2 — землетрясения, для которых 1 < Яй/Яп < 2; в йэ — 2 < Яй/Яп < Э; в й4 — Яй/Яп > Э. Зоны, в которых фиксируются землетрясения из перечисленных групп, отмечены на карте Камчатского региона, изображенной на рис. 8.
С учетом высокой скорости распространения электромагнитных сигналов в волноводе Земля — Ионосфера расстоянием между пунктами наблюдений "Микижа" и "Карымшина" (приблизительно 20 км) при расчете дистанции до эпицентра землетрясения можно пренебречь. Поэтому эпицентральное расстояние для землетрясений рассчитывалось относительно пункта наблюдений "Микижа".
Для оценки времени появления предсейсмического эффекта отдельно для каждого из рассматриваемых рядов аномалий использовался метод наложения эпох [7]. Для каждого землетрясения выбирались фрагменты рядов аномалий длиной 50 суток (25 суток до и 25 суток после землетрясения), которые впоследствии усреднялись.
Рис. 8. Зоны, в которых фиксируются землетрясения из групп Ql, Q2, 2з, ^4, на карте Камчатского региона
Результаты анализа методом наложения эпох показали, что перед землетрясениями группы Q4 аномалии в электромагнитных и геоакустических сигналах отсутствуют, поэтому они были исключены из выборки, и дальнейший анализ проводился с использованием землетрясений из групп Ql, Q2 и Qз (всего 151 землетрясение).
Установлено, что предсейсмические аномалии в сигналах ЭИ и ГАЭ наблюдаются перед землетрясениями в разные периоды времени. На рис. 9а, 10а представлены усредненные ряды аномалий сигналов ЭИ и ГАЭ, построенные для временных периодов 25 суток до и 25 суток после землетрясений. В сигналах ЭИ предсейсмические аномалии наблюдаются за 9-13 дней перед землетрясениями, максимальное значение усредненного ряда аномалий фиксируется за 10.6 дней (рис. 9а).
В сигналах ГАЭ предсейсмические аномалии регистрируются за 5-21 день до землетрясения (рис. 10а), максимальное значение усредненного ряда аномалий — за 7 дней. При этом график ряда имеет более сложную форму и содержит несколько пиков, основной из которых проявляется на интервале 5-9 дней до землетрясения.
Для проверки корректности полученных временных оценок был проведен вычислительный эксперимент. Был сгенерирован искусственный ряд землетрясений, расположенных случайным образом на временной шкале исследуемого года наблюдений.
Методом наложения эпох оценивалось время появления предсейсмических аномалий в окрестности "искусственных" землетрясений. Установлено, что в таком случае предсейсмические аномалии отсутствуют (рис. 9б, 10б). Таким образом показано, что эффект возникновения предсейсмических аномалий в сигналах ЭИ и ГАЭ неслучаен.
Совместный анализ рядов аномалий сигналов ЭИ и ГАЭ (рис. 9а, 10а), показывает, что предсейсмические аномалии в электромагнитных и геоакустических сигналах наблюдаются со смещением по времени, поэтому была произведена оценка
взаимной корреляционной функции (ВКФ) усредненных рядов. Максимальное значение ВКФ составило 0.025 при уровне статистической значимости 0.05, что свидетельствует об отсутствии линейной статистической связи.
Рис. 9. Накопленные ряды аномалий электромагнитного ОНЧ-излучения в окрестности а — реальных и б — "искусственных" землетрясений
Рис. 10. Накопленные ряды аномалий геоакустической эмиссии в окрестности а реальных и б — "искусственных" землетрясений
Оценка вероятности совместного проявления предсейсмических аномалий геоакустических и электромагнитных сигналов
Для оценки вероятности возникновения совместных предсейсмических аномалий в сигналах ЭИ и ГАЭ, рассматривались аномалии, зарегистрированные в течение 15 суток перед землетрясениями. Такой временной интервал соответствует периоду появления предвестников в исследуемых полях. Для поиска скрытых закономерностей использовались методы интеллектуального анализа (поиск последовательностей и ассоциативных правил).
Таблица 1
Значения поддержки и вероятностей для правил размера 1 и 2
Правило Поддержка, к Вероятность, Р
(общее число транзакций 151)
Размер правила 1
AI ^ Q 38 25.17%
AS ^ Q 9 5.96%
EDD ^ Q 14 9.27%
EDS ^ Q 55 36.42%
Размер правила 2
AI - EDD ^ Q 7 4.64%
AI - EDS ^ Q 30 19.87%
AS - EDD ^ Q 2 1.32%
AS - EDS ^ Q 5 3.31%
EDD - AI ^ Q 10 6.62%
EDD - AS ^ Q 0 0.00%
EDS - AI ^ Q 31 20.53%
EDS - AS ^ Q 6 3.97%
Оценка вероятности возникновения предсейсмических аномалий выполняется следующим образом:
1) Правилом размера п будем называть условие "аномалии А1, ..., Ап предшествуют землетрясению й" (А1 — ... — Ап ^ й).
2) Для каждого землетрясения строится цепочка аномалий, зарегистрированных не более, чем за 15 суток до землетрясения. Если в 15-дневном временном интервале перед землетрясением й наблюдается другое землетрясение й*, то во избежание многократного учета одного и того же правила, каждая аномалия соотносится с ближайшим по времени землетрясением. Таким образом, 15-дневный интервал перед землетрясением й сокращается до времени наступления землетрясения й*. Построенные цепочки аномалий образуют множество транзакций для процедуры поиска ассоциативных правил.
3) Формируется набор правил размера 1 (правила типа А ^ й) и предполагается, что поддержка к каждого из них равна 0.
4) Если в 15-дневном временном интервале перед й наблюдаются аномалии А1, ..., Ап (в той же временной последовательности), то поддержка правила А1 — ... — Ап ^ й увеличивается на 1.
5) Для каждого правила рассчитывается вероятность P, равная отношению поддержки k к общему количеству транзакций.
6) Правила размера n + 1 формируются путем объединения правил размера n с уровнем поддержки, превышающим заданный.
По результатам применения метода к сформированным выборкам аномалий (794 аномалии) и землетрясений (151 землетрясение) были получены следующие оценки (табл. 1): вероятность, что в 15-дневный период перед землетрясением будет зарегистрирована аномалия EDD сигнала ЭИ составляет около 9%, аномалия EDS сигнала ЭИ — около 36%, импульсная аномалия AI сигнала ГАЭ — около 25%, сплошная аномалия AS сигнала ГАЭ — около 6%. Вероятность того, что перед землетрясением будут зарегистрированы и импульсная аномалия AI сигнала ГАЭ, и аномалия EDS сигнала ЭИ в среднем составляет около 20%.
Заключение
В результате исследования определены признаки возникновения предсейсмических аномалий в электромагнитных и геоакустических сигналах, согласно которым отобраны предсейсмические аномальные проявления сигналов ЭИ и ГАЭ, зарегистрированных в 2013 году. В результате совместного анализа рядов аномалий сигналов ЭИ и ГАЭ установлено отсутствие линейной статистической связи между ними. Аномалии в сигналах обоих типов вызваны динамикой сейсмотектонического процесса и между собой могут быть связаны опосредованно. Установлено, что вероятность возникновения предсейсмических аномалий в рассматриваемых геофизических сигналах составляет не более 36%. Дальнейшие исследования предполагают использование описанных в работе методов для анализа рядов данных, охватывающих больший временной период.
Конкурирующие интересы. Авторы заявляют, что конфликтов интересов в отношении авторства и публикации нет.
Авторский вклад и ответственность. Все авторы участвовали в написании статьи и полностью несут ответственность за предоставление окончательной версии статьи в печать. Окончательная версия рукописи была одобрена всеми авторами.
Список литература/References
[1] Купцов А. В., Марапулец Ю.В., Шевцов Б.М., "Анализ изменений геоакустической эмиссии в процессе подготовки сильных землетрясений на Камчатке", Исследовано в России, 7 (2004), 2809-2818. [Kupcov A.V., Marapulets Yu.V., Shevtsov B.M., "Analiz izmenenij geoakusti-cheskoj ehmissii v processe podgotovki sil'nyh zemletryasenij na Kam-chatke", Issledovano v Rossii, 2004, №7, 2809-2818 (Russian))].
[2] Мищенко М. А., "Статистика возникновения предсейсмических аномалий в геоакустической эмиссии и атмосферном электрическом поле", Вестник КРАУНЦ. Физ.-мат. науки, 2016, № 3(14), 47-52. [Mishchenko M. A., "Statistika vozniknoveniya predseysmich-eskikh anomaliy v geoakusticheskoy emissii i atmosfernom elektricheskom pole", Vestnik KRAUNTs. Fiz.-mat. nauki, 14:3 (2016), 47:52 (in Russian)].
[3] Marapulets Yu., Solodchuk A., Shcherbina A., "Changes of geoacoustic emission directivity at "Mikizha" site associated with earthquakes in Kamchatka", E3S Web of Conferences, 11 (2016), 00014.
[4] Дружин Г. И. и др., "Акустические и электромагнитные излучения перед землетрясением на Камчатке", ДАН, 472:5 (2017), 584-589. [Druzhin G. I. et al., "Acoustic and electromagnetic emissions preceding the earthquake in Kamchatka", Doklady Earth Sciences, 472:2 (2017), 215-219].
[5] Марапулец Ю.В., Шевцов Б.М., Мезомасштабная акустическая эмиссия, Дальна-ука, Владивосток, 2012, 126 с. [Marapulets Yu.V., Shevtsov B.M., Mezomasshtabnaya akusticheskaya emissiya, Dal'nauka, Vladivostok, 2012 (in Russian), 126 pp.]
[6] Добровольский И. П., Математическая теория подготовки и прогноза тектонического землетрясения, Физматлит, М., 2009, 240 с. [Dobrovol'skiy I.P., Matem-aticheskaya teoriya podgotovki i prognoza tektonicheskogo zemletryaseniya, Fizmatlit, Moskva, 2009 (in Russian), 240 pp.]
[7] Мищенко М. А., "Статистический анализ возмущений геоакустической эмиссии, предшествующих сильным землетрясениям на Камчатке", Вестник КРАУНЦ. Физ.-мат. науки, 2011, №1(2), 56-64. [Mishchenko M.A., "Statisticheskiy analiz vozmushcheniy geoakusticheskoy emissii, predshestvuyushchikh sil'nym zemletryaseniyam na Kamchatke", Vestnik KRAUNTs. Fiz.-mat. nauki, 2011, №1(2) (in Russian), 56-64 pp.]
Vestnik KRAUNC. Fiz.-Mat. Nauki. 2021. vol. 34. no. 1. P. 137-149. TSSN 2079-6641
MSC 86-05 Research Article
Anomalies in electromagnetic and geoacoustic emission signals before Kamchatka earthquakes (ML > 4.75) in 2013
O. O. Lukovenkova, E. I. Malkin, M. A. Mishchenko, A. A. Solodchuk
Institute of Cosmophysical Research and Radio Wave Propagation, FEB RAS, 684034 Kamchatka region, Paratunka, Mirnaya str., 7, Russia E-mail: micle@ikir.ru
The paper considers pre-seismic anomalies in the signals of electromagnetic radia-tion in a very low frequency range and geoacoustic emission recorded in Kamchatka in 2013. Signs of anomaly occurrence in the studied signals and peculiarities of their manifestation were revealed. Using the epoch folding method, joint analysis of the detected pre-seismic electro-magnetic and geoacoustic anomalies was carried out. The estimates of occurrence probability for such amomalies of the considered geophysical fields are presented.
Key words: electromagnetic radiation, geoacoustic emission, pre-seismic anomaly, epoch folding method
DOT: 10.26117/2079-6641-2021-34-1-137-149
Original article submitted: 17.12.2020 Revision submitted: 05.03.2021
For citation. Lukovenkova O. O., Malkin E. T., Mishchenko M. A., Solodchuk A. A. Anomalies in electromagnetic and geoacoustic emission signals before Kamchatka earthquakes (ML > 4.75) in 2013. Vestnik KRAUNC. Fiz.-mat. nauki. 2021,34: 1,137-149. DOT: 10.26117/2079-6641-202134-1-137-149
Competing interests. The authors declare that there are no conflicts of interest regarding authorship and publication.
Contribution and Responsibility. All authors contributed to this article. Authors are solely responsible for providing the final version of the article in print. The final version of the manuscript was approved by all authors.
The content is published under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/deed.ru)
© Lukovenkova O.O. et al., 2021
Funding. The study was carried out without financial support from foundations.