Scientific journal
International Journal of Applied and fundamental research
ISSN 1996-3955
ИФ РИНЦ = 0,593

GEOPHYSICAL MONITORING OF THE GEODYNAMIC PROCESSES IN THE EARTH CRUST

Grib N.N. 1 Trofimenko S.V. 1 Trofimenko S.V. 2 Grib G.V. 1 Imaev V.S. 1 Imaev V.S. 3
1 Technical Institute branch of North-Eastern Federal University named after M.K. Ammosov
2 Institute of Tectonics and Geophysics
3 Institute of the Earth’s crust Siberian Branch Russian Academy of Sciences
2457 KB
Geophysical monitoring of the geological environment is one of the main directions in the geodynamic research which makes it possible to follow the processes in the Earth Crust. For prediction purposes it is interesting to know the statistics of the distribution of electromagnetic radiation anomalies which manifests itself in the marked periods of seismic processes. The task is to divide the sources of electromagnetic radiation into exogenous, endogenous and of anthropogenic origin. The law-based change in the distribution statistics of electromagnetic anomalies may be the result of the superposition of two processes. The total signal contains the effect from the close and far away earthquakes due to the wake-up of surface structures in the zones of dynamic influences of activated fissures and also radiation from anthropogenic seismicity which leads to the temporary surface breakup of blocks and, as a consequence, to the decrease of electromagnetic radiation.
Monitoring
geological environment
pulsed electromagnetic radiation
seismic activity

Геофизический мониторинг геологической среды – одно из основных направлений геодинамических исследований, который позволяет дистанционно отслеживать процессы в земной коре. Теоретические разработки, модельные оценки параметров и возможных источников электромагнитных вариаций, лабораторные и полевые эксперименты по изучению сейсмоэлектромагнитных явлений позволили установить физическую природу электромагнитного излучения (ЭМИ) горных пород в их естественном залегании в условиях изменяющегося деформационного процесса земных недр [1, 2, 3].

Цель работы: получить некоторую ясность в статистика распределений аномалий ЭМИ и связи их с геодинамическими процессами в земной коре.

Материалы и методы исследования

Для прогностических целей представляет несомненный интерес статистика распределений аномалий ЭМИ, реализованная на выделенных периодах сейсмического процесса. Фактически данная задача сводится к задаче разделения источников электромагнитного излучения экзогенной, эндогенной и техногенной природы.

Предварительная подготовка данных заключалась в том, что в каждой из трех компонент сигнала: «север-юг» gri1.wmf, «восток-запад» gri2.wmf и «электрическая» gri3.wmf методом фильтрации выделялись импульсные помехи (gri4.wmf), по критерию превышения среднестатистических значений максимумов суточных аномалий ЭМИ за годичный цикл измерений (gri5.wmf). В результате данной процедуры из исходной базы данных формировался файл данных, удовлетворяющий условиям:

gri6.wmf, gri7.wmf, gri8.wmf (1)

где gri9.wmf gri10.wmf gri11.wmf среднестатистические значения максимумов суточных аномалий ЭМИ за годичный цикл измерений соответствующих компонент.

Распределения аномалий ЭМИ gri12.wmf оценивалось по выборочным данным по двум пунктам регистрации электромагнитных излучений («Нерюнгри» – «Иенгра») на годичных реализациях. Средняя длина массива составила: для часовой статистики – более 100000 значений; для годовой статистики – 5602 значения.

На рис. 1 показаны результаты статистических построений в двух пунктах регистрации электромагнитных излучений («Нерюнгри» – «Иенгра»).

grib1.tif

Рис. 1. Изменение частоты появления импульсной составляющей ЭМИ в течение недели в двух пунктах регистрации «Нерюнгри-Иенгра»

Для семидневной статистики для обоих пунктов регистрации «Нерюнгри» и «Иенгра» характерно наличие максимума интенсивности в воскресенье. Причем для пункта «Нерюнгри» данный эффект наблюдается только для «С-Ю» компоненты с двукратным превышением амплитуды. Для пункта «Нерюнгри» изменение частоты появления импульсных ЭМИ представлено выпуклой функцией с минимумом активности от вторника до четверга. В пункте «Иенгра» на среду-четверг приходится максимум активности с наличием линейного тренда. Данные по двум пунктам не коррелированны (Kкор<0.2), хотя максимум в воскресение характерен для обоих пунктов регистрации.

Расположение пунктов регистрации относительно угольного разреза и систем активных разломов показано на рис. 2.

Результаты исследования и их обсуждение

На рис. 3 показана статистика распределения взрывов в Нерюнгринском угольном разрезе, расположенном в 6 – 15 км на СЗ от г. Нерюнгри (рис. 3, а) и землетрясений в радиусе 100 км от карьера (рис. 3, б). Для того чтобы данные по взрывам и землетрясениям случайным образом не перемешивались, к рассмотрению принимались взрывы, для которых разность времен пробега ts-tp=(1-2)c, а для землетрясений 2с< ts-tp <10c.

grib2.tif

Рис. 2. Расположение пунктов регистрации ЭМИ в системе активных разломов Южной Якутии

1, 2-активные разломы.: К-Кабактинский; Н-Нижне-Нерюнгринский; Б-Беркакитский; С-Суннагино-Ларбинский; Ю-Я-Южно-Якутский; В-Верхне-Гонамский; С-С (с)-Северо-Становой (северная ветвь); С-С (о)-Северо-Становой (основная ветвь); Ю-С (с)-Южно-Становой (северная ветвь); Ю-С (о)-Южно-Становой (основная ветвь); 3-эпицентры землетрясений по каталогу ГС РАН; 4-населенные пункты; 5 – технологическая площадка Нерюнгринского угольного разреза; 6 – радиус корреляции техногенной и природной сейсмичности; 7 – пункты регистрации ЭМИ

Для сравнения, на рис. 4 представлена статистика распределения землетрясений по данным из каталога ГС РАН за весь период инструментальных наблюдений (рис. 4, а) и по всем землетрясениям, зарегистрированным НСС в период 2006 – 2010 г.г. (рис. 4, б). Коэффициент корреляции для массивов данных «взрывы-все землетрясения» составляет 0.92, что и отражают формы графиков на (рис. 3, а и 4, б). Для графиков на (рис. 3, б и 4, а) коэффициент корреляции составляет 0.7, а для всех остальных случаев близок к нулю. Таким образом, вблизи карьера динамика сейсмичности в низкоэнергетической области (Е<107Дж) в значительной степени формируется под действием техногенного фактора. Для землетрясений вне радиуса корреляции распределение землетрясений на семидневных реализациях отличается равномерной статистикой с минимумом активности, приходящимся на пятницу (рис. 3, б и 4, а). Для данной статистики характерно наличие локального минимума, приходящегося на пятницу, и линейного отрицательного тренда. Причем, для землетрясений с энергией E >108Дж (рис. 4, а) угловой коэффициент равен 1.25, а для низкоэнергетических классов E<107Дж – 1.57. т.е. на максимум высвобождаемой упругой энергии в ближней зоне за счет взрывов и землетрясений (рис. 3, а и 4, б) может быть причиной снижения фоновой сейсмической активности.

grib3.tiff

Рис. 3. Распределение взрывов и землетрясений за пределами радиуса корреляции (3-10с) по данным Нерюнгринской сейсмостанции за 2006-2010 гг.

grib4.tiff

Рис. 4. Распределение землетрясений (345 событий) по данным каталога ГС РАН за 1962-2005 гг. (а) и всех землетрясений в ближней зоне (827 событий в пределах радиуса корреляции) по данным Нерюнгринской сейсмостанции за 2006-2010 гг. (б)

Закономерное изменение в статистике распределений выбросов по дням недели может быть следствием наложения двух процессов.

Во-первых, по аналогии со статистикой распределений слабых землетрясений (рис. 3, б) форма кривых подобна наблюдаемому эффекту «выходных дней» [6, 10], хотя в данном локальном случае наблюдается обратный эффект «выходных дней» т.е. на минимум сейсмической активности приходится максимум импульсных ЭМИ.

С другой стороны, известно, что распределения штормовых микросейсмических шумов и микро-землетрясений коррелируют с лунно-солнечными приливными вариациями силы тяжести. Эта избыточная сейсмическая энергия может трансформироваться в энергию электромагнитного излучения, что и наблюдается в статистике распределений регулярной составляющей gri13.wmf импульсных аномалий ЭМИ [4, 8, 9]. В контексте данного предположения аномалии ЭМИ gri14.wmf – это аналог штормовых микросейсм.

Дополнительно к изложенным причинам характерного проявления статистики распределений выбросов gri15.wmf ЭМИ суммарный сигнал может содержать эффект от ближних и дальних землетрясений за счет активизации поверхностных структур в зонах динамического влияния активизированных разломов при прохождении упругих волн. Согласно принципу суперпозиции электромагнитных волн оба эффекта могут проявляться одновременно, вследствие чего разделение суммарного сигнала в модели (1) весьма проблематично [5].

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

Например, для представленной статистики распределений «ЭМИ – сейсмичность» характерна в среднем противофазное распределение числа событий на семидневных реализациях. То-есть, увеличение числа взрывов и землетрясений в радиусе корреляции приводит к уменьшению числа выбросов (импульсов) ЭМИ gri16.wmf в пункте «Нерюнгри», в противоположность с изменением регулярной (периодической) составляющей, где зависимость прямая. Для пункта «Иенгра» уменьшение числа событий в течение недели (рис. 3, б и 4, а) приводит к положительному тренду gri17.wmf (рис. 1).

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

В блоковой структуре земной коры [7] динамика gri18.wmf изменений контролируется кинематикой взаимодействия блоков и микроблоков и пропорциональна скорости диссипативных процессов энергии упругих деформаций. Причем периодические компоненты аномалий ЭМИ будут проявляться только при наличии подвижных систем активных разломов. При уменьшении сейсмической активности gri19.wmf (числа землетрясений за определенный интервал времени) увеличивается относительная деформация, как шовных зон, так и самих блоков. В этом случае, при сохранении постоянства горизонтальной составляющей тензора напряжений возможно образование новых микротрещин или обновление уже имеющихся при достижении предельной деформации разрушения для данного уровня иерархии блоков.

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

Выводы

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

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

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

gri21.wmf, (2)

где gri22.wmf – функция предвестника сейсмической активности (gri23.wmf и gri24.wmf импульсная (выбросы) и регулярная составляющие сигнала ЭМИ); gri25.wmf – прогнозная пространственно – временная функция параметров землетрясения. Причем априори функция источника сигнала не известна, т.е. для неё не известен тип модели, определяемой по (2). Поэтому, на современном этапе развития цифровых технологий геофизического мониторинга, одной из важнейших задач становится поиск закономерностей типа (2) и разработка физических моделей на основе геофизического эксперимента.

Работа выполнена при финансовой поддержке Минобрнауки (ЗАДАНИЕ № 5.1771.2014К).