CC BY
1
Вестник КРАУНЦ. Физ.-мат. науки. 2024. Т. 47. №2. C.58-74. ISSN 2079-6641
ИНФОРМАЦИОННЫЕ И ВЫЧИСЛИТЕЛЬНЫЕ ТЕХНОЛОГИИ " https://doi.org/10.26117/2079-6641-2024-47-2-58-74 Научная статья
Полный текст на русском языке УДК 550.8; 004.8; 004.93
Поиск признаков изменения состояния приповерхностного слоя литосферы на основе результата анализа изображений, отражающих динамику характеристик сигнала геоакустической эмиссии
Ю. И. Сенкевич*
Институт космофизических исследований и распространения радиоволн ДВО РАН, 684034, Камчатский край, Елизовский район, с. Паратунка, ул. Мира, 7, Россия
Аннотация. Физическое состояние литосферы, включая ее поверхностные слои, определенным образом характеризует вероятность опасности возникновения катастрофических сейсмических событий, указывает на степень угрозы жизнедеятельности людей и опасности экономического ущерба. Одним из направлений оценки состояния литосферы является результат анализа особенностей изменчивости сигнала геоакустической эмиссии в точке наблюдения. Представлены результаты исследований, проводимых в институте космофизических исследований и распространения радиоволн Дальневосточного отделения РАН, позволяющие идентифицировать динамику характеристик сигнала геоакустической эмиссии, связанную с изменением состоянием приповерхностного слоя литосферы. Распознавание и оценка изменчивости характеристик сигналов геоакустической эмиссии на заданном временном масштабе наблюдений осуществлялась с использованием нейросетевого подхода. Разработана методика классификации наблюдаемых аномалий динамики характеристик сигнала геоакустической эмиссии.
Ключевые слова: состояние литосферы, геоакустическая эмиссия, распознавание образов, динамика характеристик сигнала, нейронные сети
Получение: 05.07.2024; Исправление: 04.08.2024; Принятие: 22.08.2024; Публикация онлайн: 25.08.2024
Для цитирования. Сенкевич Ю. И. Поиск признаков изменения состояния приповерхностного слоя литосферы на основе результата анализа изображений, отражающих динамику характеристик сигнала геоакустической эмиссии // Вестник КРАУНЦ. Физ.-мат. науки. 2024. Т. 47. № 2. C. 58-74. EDN: HGQSBR. https://doi.org/10.26117/2079-6641-2024-47-2-58-74.
Финансирование. Работа выполнена за счет Государственного задания институте космофизических исследований и распространения радиоволн Дальневосточного отделения РАН (рег. № темы 124012300245-2). Конкурирующие интересы. Конфликтов интересов в отношении авторства и публикации нет. Авторский вклад и ответственность. Автор участвовал в написании статьи и полностью несет ответственность за предоставление окончательной версии статьи в печать.
* Корреспонденция: А E-mail: senkevich@ikir.ru
Контент публикуется на условиях Creative Commons Attribution 4.0 International License © Сенкевич Ю. И., 2024
© ИКИР ДВО РАН, 2024 (оригинал-макет, дизайн, составление)
Vestnik KRAUNG. Fiz.-Mat. nauki. 2024. vol. 47. no. 2. P. 58-74. ISSN 2079-6641
INFORMATION AND COMPUTING TECHNOLOGIES " https://doi.org/10.26117/2079-6641-2024-47-2-58-74 Research Article Full text in Russian MSG 00A69, 05G75
Search signs changes in the state of the near-surface layer of the lithosphere based on the image analysis results reflecting the characteristics dynamics of the geoacoustic emission signal
Y. I. Senkevich*
Institute of Cosmophysical Research and Radio Wave Propagation FEB RAS, 684034, v. Paratunka, Mirnaya st., 7, Russia
Abstract. The lithosphere physical state, including its surface layers, in a certain way characterizes the likelihood of the risk of catastrophic seismic events, indicates the degree of threat to human life and the risk of economic damage. One of the directions for assessing the lithosphere state is the analyzing result the characteristics variability of the geoacoustic emission signal at the observation point. The studies results carried out at the Institute of Cosmophysical Research and Radio Wave Propagation Far Eastern Branch of the Russian Academy of Sciences are presented. This make it possible to identify the dynamic characteristics of the geoacoustic emission signal associated with changes in the near-surface of the lithosphere. Recognition and assessment of the characteristics variability of geoacoustic emission signals on a given time scale of observations was carried out using a neural network approach. A methodology has been developed for classifying observed anomalies in the dynamics of geoacoustic emission signal characteristics.
Key words: lithosphere state, geoacoustic emission, pattern recognition, of signal characteristics dynamics, neural networks.
Received: 05.07.2024; Revised: 04.08.2024; Accepted: 22.08.2024; First online: 25.08.2024
For citation. Senkevich Y. I. Search signs changes in the state of the near-surface layer of the lithosphere based on the image analysis results reflecting the characteristics dynamics of the geoacoustic emission signal. Vestnik KRAUNC. Fiz.-mat. nauki. 2024, 47: 2,58-74. EDN: HGQSBR. https://doi.org/10.26117/2079-6641-2024-47-2-58-74. Funding. The work was supported by the Institute of Cosmophysical Research and Radio Wave Propagation FAB RAS State Task (subject registration No. 124012300245-2).
Competing interests. There are no conflicts of interest regarding authorship and publication.
Contribution and Responsibility. The author participated in the writing of the article and is fully responsible for submitting the final version of the article to the press.
* Correspondence: A E-mail: senkevich@ikir.ru
The content is published under the terms of the Creative Commons Attribution 4.0 International License © Senkevich Y.I., 2024
© Institute of Cosmophysical Research and Radio Wave Propagation, 2024 (original layout, design, compilation)
Введение
Решение задачи поиска предсказания землетрясений и цунами на основании изучения изменений, наблюдаемых физических полей Земли, сегодня активно развивающаяся область исследований в геофизике [1-3]. В частности, изучение аномалий акустических полей в верхних слоях литосферы одно из направлений в этой области, весьма перспективно [4-6]. Исследовательской базой этого направления служит анализ динамических характеристик измеряемых параметров акустических шумов, образующихся под действием сил, порождаемых внутренними процессами Земли, которые проецируются на ее поверхность в виде продольных и поперечных волн. Возникающая при этом деформация в почвенном слое вызывает образование высокочастотного акустического сигнала, который получил название геоакустической эмиссии (ГАЭ) [7]. Сигналы ГАЭ опосредовано отражают процессы, порождающие сейсмические явления и поэтому в определенной степени могут служить базой для их описания. Преобразованные датчиками в электрические сигналы, ГАЭ представляет собой импульсный поток в диапазоне частот от единиц герц до сотен килогерц. Поскольку природа известных сейсмических явлений описывается нелинейными динамическими процессами [8], то генерируемые при этом поверхностные сейсмоакустические сигналы носят нестационарный стохастический характер. Этот факт накладывает известные ограничения на стремление получить устойчивые признаки возникающих сейсмических явлений с использованием статистических наблюдений, которые описаны моделями нелинейной динамики [9-11]. Кроме того, на результаты анализа наблюдений существенное влияние оказывают климатические и погодные условия, накладывающие на сигнал дополнительный акустический шум в точке наблюдения.
В настоящее время в Институте космофизических исследований и распространения радиоволн дальневосточного отделения Российской академии наук (ИКИР ДВО РАН) накоплена база измерений ГАЭ. На основании проводимых теоретических исследований в период с 2001 года по настоящее время сформулирована и проверяется гипотеза о существовании связи обнаруженной изменчивости характеристик сигналов ГАЭ с фиксируемыми изменениями состояния приповерхностного слоя литосферы [12-15]. Поиск подтверждающих фактов, указывающих на правомерность этой гипотезы, является основной целью проводимых исследований. Предполагается, что в силу влияния сейсмических волн на деформационные процессы в неоднородной среде литосферного слоя, в нем происходит растрескивание среды, порождающее акустические импульсы. Это явление отражается на интенсивности генерации потока импульсов акустической эмиссии определенным образом, связанным с характером сейсмических волн. В результате сейсмическое воздействие меняет состояние слоя, что отражается в динамике характеристик сигналов ГАЭ (интенсивность потока, амплитудно-фазовые показатели импульсов, скважность, частота следования, частота модуляции импульсов). Наблюдая изменение характеристик, можно попытаться выявить определенные закономерности изменений этих показателей, и связать их
с определяющим фактором, в интересуемом случае - с параметрами сейсмических волн. Изменения характеристик сигнала ГАЭ с 2001 года фиксируются в Лаборатории акустических исследований ИКИР ДВО РАН в ходе наблюдений, как изменение статистических распределений амплитуды, межимпульсных интервалов и концентрации импульсов во времени. В общем случае задача исследования формулируется, как стремление качественно описать наблюдаемую изменчивость свойств среды под действием деформаций, вызванных воздействием внешних факторов. Изменения свойств среды фиксируются опосредовано через обнаружение аномального поведения динамики характеристик сигналов ГАЭ.
Пусть за некоторое время в результате измерений получаем вектор значений определенного параметра сигнала ГАЭ (например, серии межимпульсных интервалов, амплитуд, значений частотной или амплитудной модуляции, или их производных) X = (Х1, Х2,..., Хи). Обозначим место вектора, как точку в и-мерном фазовом пространстве. Получим п последовательных измерений и разместим их в том же фазовом пространстве. В результате получим фазовую траекторию, описывающую динамику выбранного параметра за время Т = пД"Ь. Очевидно, что траектория однозначно описывается матрицей вида
fXi X2 u Xn
U=4
. (1)
Х1,2 Х2,2 ■ ■ ■ Хп,2
,Х1,1 Х2,1 ••• Хп,1
Будем называть (1) отображением динамики выбранного параметра сигнала ГАЭ с дискретностью в масштабе измерения Т, или матрицей-отображений динамики параметра, или просто — отображением. В ходе наблюдений в фазовом пространстве складывается непрерывная картина динамики характеристик ГАЭ, которую можно представить как последовательное соединение фрагментов траектории вида (1). Согласно законам нелинейной динамики, состояния открытых систем отображаются в фазовом пространстве некоторыми петлями сложных траекторий [16-18]. При достаточно длительном наблюдении за изменениями состоянием системы, следует ожидать появление самоподобных траекторий, имеющими сходное описание. Этот факт создает прецедент для поиска закономерности смены состояния системы и, как следствие, возможность обнаружения закономерностей в сигнале на уровне описания динамики его характеристик через отображения. Степень сходства описания можно определить через сравнение отображений вида (1). Однако в данном случае при выборе метода сравнения отображений следует учитывать специфику их построения, связанную с порядком их составления из векторов параметров.
Ц={ аьа2,аэ, ••• ап} , (2)
где аг = (Хг,)}, ) = 1,и. То есть, существует определяющая порядок построения отображения временная ось, и поэтому корректное сравнение должно осуществляться синхронным попарным сопоставлением векторов сравниваемых
отображений вдоль временной оси. Примем следующее определение степени сходства отображений: два равно размерных отображения А и В подобны до степени Н, равного сумме коэффициентов корреляции Пирсона попарно выбираемых векторов сравниваемых матриц щ и Ъ вида (2):
Н(А, В) = ^ т(аьЪ0 = ^
г=1
г=1
СОУ (аьЪ0
(3)
где ста, стЪ — дисперсии векторов а и Ъ, соответственно. Из свойства корреляции вытекает, что максимальное значение степени сходства отображения шириной в п векторов равно Нтах = |пх Н|. Если Н(А, В) = —Нтах, то каждая из матриц является негативным отображением по отношению к другой. Если Н(А, В) = 0, это означает, что отображения полностью не совпадают.
Разбивая эпоху наблюдения сигнала ГАЭ на последовательные фрагменты равной длительности Т, сформируем множество отображений (ит), отражающих динамику выбранного параметра сигнала ГАЭ за эпоху анализа длительностью Т х М. Сравнивая между собой попарно отображения полученного множества фрагментов по критерию (3) получим симметричную матрицу значений тк(к = 1, М) степеней совпадений множества отображений П (4):
1
Т1,2
Т2,1 1
тп—1,1
тк,1
т1,к—2 т2,к—2 Т1,к— 1 Т2,к—1 т1,к т2,к
Тп—1,2 Тк,2
тк—1,к—2 тк,к—2
1 тк,к—1 тк—1,к 1
(4)
Зададимся значением порога сходства отображений Л0, и применим правило отбора сравниваемых значений 5 построчно:
V =
1, Л > Ло 0
(5)
В результате получим бинарную матрицу, где любая строка, содержащая хотя бы одну <1> вне главной диагонали, будет обозначать кластер, объединяющий все отображения с номерами элементов этой строки, в которых стоят <1>. Представленная последовательность вычислений по существу представляет алгоритм кластеризации, в котором критерием близости выступает корреляционная связь отображений вида (3). Формирование матрицы (4) позволяет выделить кластеры отображений с характерными особенностями динамики характеристик ГАЭ для их дальнейшего сопоставления с зафиксированными ранее сейсмическими событиями.
Таким образом, в настоящем исследовании задача обнаружения смены состояния среды фиксируется на основании распознавания отображений динамики характеристик наблюдаемого сигнала ГАЭ вида (4). Предшествующий
распознаванию отображений этап их разбиения на классы носит трудоемкий методический и вычислительный характер. Среди эффективных алгоритмов кластеризации в современной обработке отображений выделяются инструменты, разработанные в рамках искусственного интеллекта. Достижения в этой области на основе структур нейронных сетей показали высокую эффективность их внедрения в практику цифровой обработки сигналов и распознавания образов. Что и явилось мотивацией к использованию нейронных сетей в настоящем исследовании на этапах препроцессинга сигналов, формирования тензоров отображений и их кластеризации.
Методика эксперимента
Особенности приема сигнала ГАЭ вблизи поверхности почвенного слоя земли, описанные в ряде более ранних работ [19,20], состоят в том, что расчет вычисления статистических распределений амплитудно-частотных характеристик импульсов заменен на распределение характеристик, полученных по результатам сравнения структуры этих импульсов [21]. В проводимых нами исследованиях на этапе предобработки импульсного сигнала ГАЭ амплитудно-фазовая структура каждого составляющего сигнал импульса претерпевает инвариантное преобразование, позволяющее описать его структуру в виде бинарной матрицы. Это преобразование позволяет провести сравнительный анализ структур последовательности импульсов в сигнале ГАЭ и их кластеризацию по степени сходства матриц. В результате на фрагменте сигнала заданной длительности удается описать и определить численно разнообразие структур импульсов, а далее сформировать множество, составляющее классы - алфавиты. Каждый элемент алфавита связан со статистикой его встречаемости на фрагменте сигнала и представляет объект, включающий амплитудно-фазовую характеристику импульса. Далее выполняется ранжирование элементов алфавита по значению размерности символов и степени сходства их матриц [21]. Весь анализируемый сигнал разбивается на равные фрагменты, для которых вычисляются отдельные алфавиты и выстраивается их распределение по времени. Каждое полученное распределение элементов алфавита выкладывается на двумерную сетку, где по одной оси вносятся элементы алфавита по их размерности, а по другой выстраиваются алфавиты, полученные на последовательно нарезанных фрагментах сигнала. В итоге формируется трехмерное отображение изменчивости состава алфавитов по времени, где третья ось представлена значениями частоты встречаемости элементов определенного алфавита для каждого фрагмента сигнала. Таким образом сформированная матрица представляет отображение и позволяет наблюдать динамику характеристик сигнала ГАЭ. В данном случае в символы алфавита "упакованы"амплитудно-фазовые характеристики сигнала ГАЭ, которые последовательно отражены в алфавитах от одного временного фрагмента к другому.
Для проверки высказанного предположения о существовании связи сейсмического воздействия на характеристики ГАЭ в ходе проводимого исследования была разработана методика, контроля сезонных изменений состояния приповерхностного слоя литосферы. Методика включала следующие основные этапы: предварительную подготовку сигналов ГАЭ с использованием авторского метода структурно-лингвистической обработки; формирование трехмерных отображений, отражающих динамику характеристик сигнала ГАЭ на заданных временных интервалах и подготовку этих отображений для передачи в нейронную сеть; кластеризация подготовленных отображений и выделение классов аномалий динамики характеристик ГАЭ; обучение нейронной сети на основании выделенных групп отображений, проверка результата обучения и сопоставление результатов классификации с сезонным состоянием поверхностного слоя литосферы в районе измерений.
Первый этап методики предусматривает подготовку данных и начинается с применения структурно-лингвистической обработки временных файлов, содержащих сигналы ГАЭ, которые были получены в ходе наблюдений на Камчатском геодинамическом полигоне в рамках выполнения государственного задания ИКИР ДВО РАН с 2016 по 2021 годы [22]. Каждый файл это отрезок сигнала, представленный 15 минутным его фрагментом. Длительность фрагмента принята в лаборатории акустических исследований ИКИР ДВО РАН за условный стандарт — минимально достаточный блок для обработки данных выбранными методами. Файл содержит от нескольких десятков до нескольких тысяч импульсов ГАЭ. В результате структурно-лингвистической обработки из каждого файла выделяется алфавит. После сортировки по значениям размерности элементов каждый алфавит представляется как вектор. Далее составляется ряд векторов, полученных из последовательно следующих друг за другом фрагментов сигнала. В памяти компьютера формируется блок размерностью 512 на 2880 точек, который визуально может быть представлен трехмерным графиком, составленным из 2880 векторов. Пример трех сформированных отображений представлен на рис. 1.
Рис. 1. Контроль результата работы функции предобработки отображений [Figure 1. The controls result of the image preprocessing function]
На отображениях снизу вверх размещены элементы алфавитов по значениям их размерностей, слева направо приведена временная шкала, цветом отражены значения частоты появления элементов в алфавитах. Выбор ряда из 2880 векторов определен средне месячным периодом анализа наблюдений (4 файла * 24 часа * 30 суток) [24].
Отображение динамики характеристик ГАЭ в течении месяца позволяет оценить сезонную изменчивость состояния верхних слоев литосферы, достаточную для демонстрации работоспособности подхода, заложенного в представляемой методике.
Общее число подготовленных таким образом данных за период наблюдений с 2017 по 2021 годы составил 60 графиков (12 месяцев * 5 лет). Практика исследований показала, что внешнее представление полученных графиков выглядит как отображение со сложным цветным узором. Визуальный анализ полученного материала практически затруднителен в виду большого разнообразия трудно подающихся описанию цветных паттернов в поле отображениях. В обход трудностей зрительного анализа было принято собрать тензор данных для обработки отображений нейросетью. Путем интерполяции подготовленные отображения преобразуются в растровое представление в формате Portable Network Graphics (PNG), масштабируются до размеров 64х64 пикселя. Формат новых отображений выбирается исходя из подгонки их к доступным параметрам входного слоя используемой архитектуры нейронной сети. Размер получившегося массива отображений недостаточен для качественной обработки данных нейронной сетью (60 отображений). По этой причине использовался метод генерации предыскаженных клонов отображений на базе существующих отображений. Варианты предыскажений отображений ограничены из-за вероятности потерять полезную информацию. которая привязана к масштабу и векторному построению исходных отображений. Использовались преобразования и трансформации, которые только смещали отображение по оси времени, а также производилась линейная корректировка яркости и контрастности, поскольку такие изменения допустимы. В этом состояла одна из особенностей обработки исходной базы отображений. Для каждого исходного подготовленного отображения была создана серия из 50 дубликатов отображений с предыскажениями, заложенными в генераторе аугментации. В результате выполненной предобработки итоговое число нового массива отображений составило 3000. Исходные 60 неискаженных отображений используются только на заключительном этапе распознавания отображений по классам.
Второй этап обработки данных направлен непосредственно на определение особенностей динамики характеристик сигналов ГАЭ путем выделение специфики паттернов полученных отображений, передаваемой изменениями палитры рисунка, и дальнейшем разбиении отображений на группы по степени сходства пространства признаков применением алгоритмов кластерного анализа.
Отображения относятся к образам так называемых «тепловых карт». Распознавание таких отображений имеет свои особенности, связанные с
отсутствием на них определенных объектов. Место объектов занимают градиенты цветности или уровней серого, которые несут основную информационную нагрузку. Для выполнения кластеризации было принято решение провести воспользоваться алгоритмом Ллойда — (метод K-средних), использующий минимум-дисперсионный подход. Для возможности сравнения распределений по кластерам принудительное разбиение осуществлялось от 2 до 6 кластеров. Результат выполнения кластеризации представлен в табл. 1.
Таблица 1
Результаты работы алгоритма кластеризации [The results of the clustering
algorithm]
Заданное число кластеров Полученное число элементов в кластерах
2 929 2071
3 808 907 1285
4 800 863 523 814
5 516 863 271 788 562
6 861 336 561 436 538 268
Для оценки эффективности проведенного кластерного анализа использовался Коэффициент силуэта, позволяющий измерить компактность и разделение кластеров. Он дает оценку в диапазоне от -1 до 1, где более высокое значение указывает на лучшее качество кластера. Положительное значение указывает на то, что образцы хорошо соответствуют выделенным кластерам, а отрицательное значение предполагает, что они, возможно, были отнесены к неправильным кластерам. Показатели оценки кластеризации представлены в табл. 2.
Таблица 2
Результаты оценки эффективности кластерного анализа [The results of evaluating the effectiveness of cluster analysis]
Число полученных кластеров Значение Контурного коэффициента
2 0,722281
3 0,953454
4 0,859461
5 0,756904
6 0,745786
По результатам второго этапа выполняемой методики можно сделать вывод, что разделение на группы для всех вариантов числа заданных кластеров прошло успешно. При этом явно выделяется максимум показателя качества кластеризации для значения 3. Число элементов в указанном случае равномерно распределено по кластерам.
Следующий этап представляемой методики состоит в обучении нейронной сети на основании выделенных групп отображений по лучшему результату оценки качества кластеризации. Затем выполняется проверка результата обучения с использованием выборки 60 исходных неискаженных отображений на предмет качества классификации используемой нейронной сети. В качестве обучающей выборки данных был выбран результат кластеризации с 3 выделенными группами, файлам которых были присвоены соответствующие метки со значениями классов 0, 1, 2.
Для проведения эксперимента была разработана нейронная сеть с архитектурой типа КезЫе^ Для эффективной работы нейросети были проведены многочисленные эксперименты с изменением глобальных параметров, численностью сверточных слоев, числа остаточных связей, числа фильтров. Найденная оптимальная архитектура нейросети, на которой были получены лучшие результаты обучения и классификации представлена на рис. 2.
Layer (type) Output Shape Paiam # Connected to
img [InputLayer] [[More G4, 64, 3]] и []
conv2d_75 [Conv2D) (None, 62, 62, 32] S96 ['img[0][В]']
conv2d_76 (ConvZD) (None, 60. 60, 32] 9248 [ conv2d_75[0]ff]']
max_pooling2d_12 {MaxPooli ng2D) (None, 20, 20, 33) 0 [■conv2d_76[0][0] ]
conv2d_77 (Conv2D) (None, 20, 20, 32] 924S ['max_poolïng2d_12[0][0]1]
canv2d_78 (ConvJD) (None, 20, 20, 32] 924В fconv2d,77[0] [0] 1 ]
add_2G (Add) (None„ 20. 20. 32] 0 [1canv2d_7B[0][011, nia)c_paoling2d_i2[0] [0]1 ]
conv2d_79 (Conv2D) (None„ 20. 20. 32] 9246 ["add_20[0][0]'J
conv2d_80 (Conv2D) (None, 20, 20, 32] 9248 ['conv2d_79[0][0]1]
add_21 (Add) (None, 20, 20, 33] [1conv2d_80(0][0]1, 1add_20[0][0]']
conv2d_81 (Conv2D) (None, IS, ie, 33] 9248 [ 'add_21[0][0]1]
globally erage_pooling2d_l 2 (GlobalAveiagePoolingiD) (None, 32] 0 [■conv2d_Sl[0][0]1]
dense_24 (Dense] (None, 256) 344 S ['global average__pooling2d 12[ 0][0]']
dropout_l2 [Dropout) (None, 256) 0 [ dense_24[0][0] J
dense_25 (Dense] (None„ 3] 771 ["dropout_12[01[0]']
Рис. 2. Описание архитектуры нейросети [Figure 2. Description of the neural network architecture]
Результаты
Результаты тестирования нейросети следующие: ошибка классификации на тестовой выборке составила 0,3426; точность классификации на тестовой выборке составила 0,9133. Для всех классов нейросеть показала уверенное распознавание отображений по классам обучения с малыми значениями ошибок классификации уже на 40 эпохе обучения. Нормированная матрица ошибок представлена на рис. 3.
Верные классы 0 0,89 0,07 0,04
1 0,01 0,98 0,01
2 0,1 0,03 0,87
0 1 2
Предсказанные классы
Рис. 3. Матрица ошибок классификации [Figure 3. Classification error matrix]
Нейросеть верно распознала 89% отображений в Классе 0, 98% отображений в Классе 1 и 87% отображений в Классе 2. Таким образом средняя точность распознавания составила 91 %. Графики обучения модели на рис. 4 и 5 показывают ровную тенденцию без всплесков переобучения.
Рис. 4. Графики, показывающие долю верных ответов в процессе обучения модели на 64 эпохах
[Figure 4. Graphs showing the proportion of correct answers in the learning process of the model at 64 epochs]
Рис. 5. Графики ошибки в процессе обучения модели на 64 эпохах [Figure 5. Error graphs in the learning process of the model at 64 epochs]
Заключительный этап методики решает задачу сопоставления результатов классификации с измененным состоянием поверхностного слоя литосферы в районе измерений. Этот этап осуществляется визуальной оценкой особенностей отображений, которые в ходе выполнения методики были разделены на кластеры, а затем был собран классификатор состояний литосферы на основе нейросети.
Рис. б. Примеры отображений, принадлежащих определенным классам [Figure 6. Examples of images belonging to different classes]
На рис. 6 представлены ряды по 10 экземпляров отображений, представляющих определенный класс. Можно видеть достаточно заметные отличия паттернов отображений для каждого из классов. Однако сделать убедительные выводы по состоянию поверхностного слоя литосферы на основании только наблюдаемой динамики характеристик сигнала ГАЭ по паттернам отображений и времени проявления паттернов оказалось недостаточно. Здесь в явном виде помимо
нелинейности самих сейсмических явлений, формирующих специфические паттерны на отображениях, значительно влияют погодные и сезонные природные влияния в виде осадков, температурных градиентов и перепадов атмосферного давления, которые были ранее зафиксированы в наших исследованиях (например,
[24]).
Заключение
Проведенное исследование показало, что разработанная методика позволяет выйти на новый уровень оценки динамики характеристик сигналов ГАЭ, представленных в виде отображений. Анализ отображений с использованием алгоритмов кластеризации и классификации, построенных на базе нейронных сетей, дает возможность уверенно обнаружить общие и отличительные черты и закономерности, скрытые в импульсном потоке сигналов ГАЭ, несмотря на то, что исходные обрабатываемые отображения не имеют строго обозначенных контурных очертаний, а это в значительной степени затрудняет формирование пространства отличительных признаков, на основании которых формируются кластеры. К сожалению малое число исходных отображений не позволяет провести достаточно глубокий анализ, как с позиции применяемых методов предобработки данных, на основе использования нейросетевого подхода, так и с позиции применения различных методов кластеризации. Однако разработанная методика, и разработанные в ходе ее реализации компьютерные программы, позволяют проводить дальнейшие эксперименты в автоматическом режиме на всех этапах от формирования датасета и до вывода результатов распознавания и классификации. Разработанный классификатор на базе архитектуры успешно справляется
с поставленной технической целью и обеспечивает перспективу развития архитектуры в сторону ее адаптации для решения задачи поиска признаков изменчивости состояния приповерхностного слоя литосферы.
Аббревиатуры
ГАЭ Геоакустическая эмиссия Список литературы
1. Богомолов Л.М., Сычева Н. А. Прогноз землетрясений в XXI веке: предыстория и концепции, Геосистемы переходных зон, 2022. Т. 6, №3, С. 145-181 Б01: 10.30730/^.2022.6.3.145-164.164-182.
2. Родкин М.В., Липеровская Е. В. Проблемы и новый подход к прогнозу сильных землетрясений / Проблемы комплексного геофизического мониторинга сейсмоактивных регионов, Девятая Всероссийская научно-техническая конференция с международным участием, 24-30 сентября 2023 г. Тезисы докладов.. Петропавловск-Камчатский, Камчатский филиал ФИЦ ЕГС РАН, 2023, С. 24.
3. Шебалин П.Н., Гвишиани А. Д., Дзебоев Б. А., Скоркина А. А. Почему необходимы новые подходы к оценке сейсмической опасности?, Доклады Российской академии наук. Науки о Земле, 2022. Т. 507, №1, С. 91-97 Б01: 10.31857/Б2686739722601466.
4. Хайретдинов М.С., Ковалевский В. В., Воскобойникова Г. М., Седухина Г. Ф., Якименко А. А. Сейсмоакустические волны сейсмовибраторов в системе литосфера-атмосфера, Вестник НЯЦ РК, 2018. Т. 2, С. 44-48.
5. Спивак А. А., Рыбнов Ю. С. Акустические эффекты сильных землетрясений, Физика Земли, 2021, С. 41-50 DOI: 10.31857/S0002333721010099.
6. Муратов П. В., Руленко О. П., Марапулец Ю.В., Солодчуку А. А. Электрический и акустический отклик приповерхностных осадочных пород на прохождение сейсмических волн от землетрясений, Вестник КРАУНЦ. Физ.-мат. науки, 2018. №5(25), С. 62-73 DOI: 10.18454/2079-6641-2018-25-5-62-73.
7. Гапеев М. И., Марапулец Ю. В. Моделирование зон относительных сдвиговых деформаций перед сильными землетрясениями на Камчатке, произошедшими в период 2018-2021 гг, Вестник КРАУНЦ. Физ.-мат.. науки, 2021. Т. 37, №4, С. 53-66 DOI: 10.26117/2079-6641-2021-37-4-53-66.
8. Соломин С. В. Физика земли, Учебное пособие для студентов, обучающихся по специальности 011200 - геофизика., Т. 2 (Сейсмология). Саратов: Из-во «Научная книга», 2008. 42 с.
9. Самсонов В. Б. и др. Нелинейная динамика Земли: сферы и структуры самоорганизации, монография: под ред. В. Б. Самсонова. Саратов: Совместный учеб.-науч. центр ин-та аграрных проблем РАН и Саратовского гос. ун-та им. Н. Г. Чернышевского "Социальная география и региональный анализ", 2005. 217 с.
10. Мазуров Б. Т. Математическое моделирование при исследовании геодинамики, монография:. Новосибирск: Агентство «Сибпринт», 2019.360 с.
11. Сычев В.Н., Богомолов Л.М. Потенциал методов нелинейной динамики для анализа геофизических рядов и сейсмичности / Труды конференции Солнечно-земные связи и физика предвестников землетрясений VI международная конференция 9-13 сентября 2013 г., Тезисы докладов. Паратунка, Камчатский край, Институт космофизических исследовангий и распространения радиоволн ДВО РАН, 2013.
12. Марапулец Ю.В., Сенкевич Ю.И., Луковенкова О.О., Солодчук А.А., Ларионов И.А., Мищенко М.А., Малкин Е.И., Щербина А.О., Гапеев М.И. Комплексный анализ акустических и электромагнитных сигналов для оценки уровня сейсмической опасности. Владивосток: Дальнаука, 2020.120 с.
13. Сенкевич Ю.И. и др. Методика выделения информативных признаков в сигналах геоакустической эмиссии, Труды СПИИРАН, 2019. Т. 18, №5, С. 1066-1092 DOI: 10.15622/sp.2019.18.5.1066-1092.
14. Сенкевич Ю. И., Мищенко М. А., Солодчук А. А., Луковенкова О. О., Гапеев М. И. Стохастическое моделирование сигнала геоакустической эмиссии, Цифровая обработка сигналов, 2023. №1, С. 27-34.
15. Игнатьев Ю.Г. Математические модели теоретической физики с примерами решения задач в СКМ Maple., Издание второе, исправленное и дополненное. Казань: Казанский университет, 2023. 248 с.
16. Кузнецов А. П. Динамические системы и бифуркации. Саратов: Издательский центр «Наука», 2015.168 с.
17. Малевич Н.Е. Нелинейная динамика, учеб. пособие. Самара: Изд-во Самар. гос. аэрокосм ун-та, 2007.160 с.
18. Сенкевич Ю. И., Мищенко М. А. Методика оценки состояния приповерхностных осадочных пород по результатам наблюдений динамических характеристик геоакустической эмиссии, Вестник КРАУНЦ. Физ.-мат. науки. Т. 45, №4, С. 109-121 DOI: 10.26117/2079-66412023-45-4-109-121.
19. Сенкевич Ю.И., Марапулец Ю.В., Луковенкова О. О., Солодчук А. А. Методика выделения информативных признаков в сигналах геоакустической эмиссии, Труды СПИИРАН. №18(5), С. 1066-1092 DOI: 10.15622/sp.2019.18.5.1066-1092.
20. Сенкевич Ю.И., Луковенкова О. О., Солодчук А. А. Методика формирования Реестра геофизических сигналов на примере сигналов геоакустической эмиссии, Геосистемы переходных зон, 2018. Т. 2, №4, С. 409-418 DOI: 10.30730/2541-8912.2018.2.4.409-418.
21. Сенкевич Ю.И. и др. Способ обнаружения комплексного предвестника землетрясений, пат. № 2758582. Рос. Федерация G01V 11/00, G01V 1/00, G01V 3/12. ФГБУН ИКИР ДВО РАН. № 2020138668; 26.11.2020, Бюл. № 31.
22. Мищенко М. А., Ларионов И. А., Щербина А. О. База данных «Реестр сейсмоакустических сигналов по наблюдениям в долине реки Карымшина, св-во о регистрации баз данных 2023620001: ФГБУН ИКИР ДВО РАН от 09.01.2023.
23. Senkevich Yu.I. Search for Hidden Patterns in Acoustic and Electromagnetic Pulse Signals, IEEE, 2020 DOI: 10.1109/SCM50615.2020.9198754.
24. Марапулец Ю.В., Солодчук А. А. Суточный ход температуры как причина возникновения периодической вариации высокочастотной геоакустической эмиссии, Метеорология и гидрология, 2016. Т. 4, С. 37-44.
Информация об авторе
1 Сенкевич Юрий Игоревич^ - доктор технических наук, доцент, ведущий научный сотрудник, Лаборатории акустических исследований, Институт космофизических исследований и распространения радиоволн ДВО РАН, Паратунка, Россия, © СЖСГО 0000-0003-0875-6112.
References
[1] Bogomolov L. M., Sycheva N.A. Prognoz zemletryaseniy v XXI veke: predystoriya i kontseptsii. Geosistemy perekhodnykh zon, 2022, no. 6(3) pp. 145-182 DOI: 10.30730/gtrz.2022.6.3.145-164.164-182 (In Russian).
[2] Rodkin M.V., Liperovskaya E.V. Problemy i novyy podkhod k prognozu sil'nykh zemletryaseniy. Problemy kompleksnogo geofizicheskogo monitoringa seysmoaktivnykh regionov. Devyataya Vserossiyskaya nauchno-tekhnicheskaya konferentsiya s mezhdunarodnym uchastiem. 24 - 30 September 2023. g. Tezisy dokladov. Petropavlovsk-Kamchatskiy: Kamchatskiy filial FITs EGS RAN, 2023. 24 p. (In Russian).
[3] Shebalin P. N., Gvishiani A. D., Dzeboev B. A., Skorkina A. A. Why are new approaches to seismic hazard assessment required? Doklady Rossiyskoy akademii nauk. Nauki o Zemle, 2022. vol. 507, no 1, pp. 91-97. DOI: 10.31857/S2686739722601466. (In Russian).
[4] Khayretdinov M.S., et al. Seismic-acoustic waves of seismic vibrators in the lithosphere-athmosphere system. Vestnik NYaTs RK, 2018. vol. 2, pp. 44-48. (In Russian).
[5] Spivak A. A., Rybnov Y. S. Acoustic effects of strong earthquakes. Izvestiya, Physics of the Solid Earth. 2021. vol. 57. pp. 37-45.
[6] Muratov P. V., Rulenko O.P., Marapulets Yu. V., Solodchukv A. A. Elektricheskiy i akusticheskiy otklik pripoverkhnostnykh osadochnykh porod na prokhozhdenie seysmicheskikh voln ot zemletryaseniy. Vestnik KRAUNC. Fiz.-mat. nauki, 2018. no. 5(25), pp. 62-73, DOI: 10.18454/2079-6641-2018-25-5-62-73 (In Russian).
[7] Gapeev M.I., Marapulets Yu. V. Modelirovanie zon otnositel'nykh sdvigovykh deformatsiy pered sil'nymi zemletryaseniyami na Kamchatke, proizoshedshimi v period 2018-2021 gg. Vestnik KRAUNC. Fiz.-mat. nauki, 2021. vol. 37, no. 4, pp. 53-66, DOI: 10.26117/20796641-2021-37-4-53-66 (In Russian).
[8] Solomin S. V. Fizika zemli, Uchebnoe posobie dlya studentov, obuchayushchikhsyapo spetsial'nosti 011200 - geofizika., vol. 2 (Seysmologiya). Saratov: Iz-vo «Nauchnaya kniga», 2008. 42 p. (In Russian).
[9] Samsonov V. B., et al. Nelineynaya dinamika Zemli: sfery i struktury samoorganizatsii, monografiya: pod red. V. B. Samsonova. Saratov: Sovmestnyy ucheb.-nauch. tsentr in-ta agrarnykh problem RAN i Saratovskogo gos. un-ta im. N. G. Chernyshevskogo "Sotsial'naya geografiya i regional'nyy analiz 2005. 217 p. (In Russian).
[10] Mazurov B.T. Matematicheskoe modelirovanie pri issledovanii geodinamiki [Mathematical modeling in the study of geodynamics]. Novosibirsk: Agentstvo «Sibprint», 2019. 360 p. (In Russian).
[11] Sychev V. N., Bogomolov L.M. Potentsial metodov nelineynoy dinamiki dlya analiza geofizicheskikh ryadov i seysmichnosti. Trudy konferentsii Solnechno-zemnye svyazi i fizika predvestnikov zemletryaseniy VI mezhdunarodnaya konferentsiya 9-13 September 2013. 2013. (In Russian).
[12] Marapulets Yu. V., et al. Kompleksnyy analiz akusticheskikh i elektromagnitnykh signalov dlya otsenki urovnya seysmicheskoy opasnosti [Comprehensive analysis of acoustic and electromagnetic signals for assessing the level of seismic hazard]. Vladivostok: Dal'nauka. 2020. 120 p. (In Russian).
[13] Senkevich Yu. I., et al. Technique of informative features selection in geoacoustic emission signals. SPIIRAS Proceedings, 2019. vol. 18, no. 5, pp. 1066-1092 DOI: 10.15622 /sp.2019.18.5.1066-1092. (In Russian).
[14] Senkevich Yu. I. Stochastic modeling of geoacoustic emission pulse signal. Digital signal processing. 2023. no. 1. pp. 27-34. (In Russian).
[15] Ignat'ev Yu.G. Matematicheskie modeli teoreticheskoy fiziki s primerami resheniya zadach v SKM Maple [Mathematical models of theoretical physics with examples of solving problems in SCM Maple]. Izdanie vtoroe, ispravlennoe i dopolnennoe. Kazan': Kazanskiy universitet, 2023. 248 p. (In Russian).
[16] Kuznetsov A.P. Dinamicheskie sistemy i bifurkatsii []. Saratov: Izdatel'skiy tsentr «Nauka», 2015. 168 p. (In Russian).
[17] Malevich N.E. Malevich N., E. Nonlinear dynamics [Nonlinear dynamics]. Samara: SamGAU, 2007. 160 p. (In Russian).
[18] Senkevich Yu. I., Mishchenko M.A. Metodika otsenki sostoyaniya pripoverkhnostnykh osadochnykh porod po rezul'tatam nablyudeniy dinamicheskikh kharakteristik geoakusticheskoy emissii. Vestnik KRAUNTs. Fiz.-mat. nauki. 2023. vol. 45, no. 4, pp. 109-121, DOI: 10.26117/2079-6641-2023-45-4-109-121 (In Russian).
[19] Senkevich, Yu.I., Marapulets, Yu.V., Lukovenkova, O.O., Solodchuk, A.A. Metodika vydeleniya informativnykh priznakov v signalakh geoakusticheskoy emissii. Trudy SPIIRAN. no. 18(5), pp. 1066-1092 DOI: 10.15622/sp.2019.18.5.1066-1092 (In Russian).
[20] Senkevich Yu.I., Lukovenkova O.O., Solodchuk A.A. Metodika formirovaniya Reestra geofizicheskikh signalov na primere signalov geoakusticheskoy emissii. Geosistemy perekhodnykh zon. vol. 2. no. 4, pp. 409-418 DOI: 10.30730/2541-8912.2018.2.4.409-418 (In Russian).
[21] Senkevich Yu.I., Marapulets Yu.V., Lukovenkova O.O., Solodchuk A.A., Mishchenko M.A., Malkin E.I., Gapeev M.I. Sposob obnaruzheniya kompleksnogo predvestnika zemletryaseniy, patent 2758582, Ros.Federatsiya G01V 11/00, G01V 1/00, G01V 3/12 (FGBUN IKIR DVO RAN. № 2020138668; 26.11.2020, Byul. № 31. 9 s.) (In Russian).
[22] Mishchenko M.A., Larionov I.A., Shcherbina A.O. Baza dannykh «Reestr seysmoakusticheskikh signalov po nablyudeniyam v doline reki Karymshina, FGBUN IKIR DVO RAN ot 09.01.2023, registration certificate Data Base 2023620001. (In Russian).
[23] Senkevich Yu.I. Search for Hidden Patterns in Acoustic and Electromagnetic Pulse Signals, 2020 XXIII International Conference on Soft Computing and Measurements (SCM). IEEE, 2020, 2020 DOI: 10.1109/SCM50615.2020.9198754 (In Russian).
[24] Marapulets Yu.V., Solodchuk A.A. Sutochnyy khod temperatury kak prichina vozniknoveniyaperiodicheskoy variatsii vysokochastotnoy geoakusticheskoy emissii. Meteorologiya i gidrologiya. 2016. vol. 4, pp. 37-44. (In Russian).
Information about the author
Senevich Yury Igorevich & - D. Sci. (Tech)),
Docent, Leading Researcher, Laboratory of Acoustic
Research, Institute of Cosmophysical and Radio
Wave Propagation FEB RAS, Paratunka, Russia, ©ORCID 0000-0003-0875-6112.
