Моделирование направления осей главных напряжений горных пород при подготовке землетрясений
Authors: М.И. Гапеев , А.А. Солодчук 
Published: November 2024

Основная гипотеза возникновения предсейсмических аномалий направленности геоакустической эмиссии состоит в том, что на заключительной стадии подготовки землетрясения формируются постоянные направления осей главных напряжений пород в точке наблюдений. Направления этих осей в свою очередь определяют преимущественную ориентацию акустических источников. Для подтверждения этой гипотезы произведено моделирование ориентации осей главных напряжений, обусловленных процессом подготовки землетрясений. Расчеты основаны на модели, построенной в рамках линейной теории упругости: земная кора рассматривается в виде однородного изотропного упругого полупространства, силовое воздействие в очаге готовящегося землетрясения — в виде комбинации двойных пар сил. Учтена потенциальная энергия упругих деформаций, накопленная в процессе подготовки землетрясения. В работе использован каталог механики очагов землетрясений «The Global Centroid-Moment-Tensor Catalog». Из него были выбраны параметры сейсмических событий, произошедших вблизи полуострова Камчатка с 1976 по 2020 годы. Так как ориентация акустического излучения зависит от азимутального направления на эпицентры землетрясений, то все рассматриваемые сейсмические события были разделены на три группы методом K-средних по пространственному расположению их эпицентров. Результаты моделирования сопоставлены с экспериментальными оценками направления осей главных напряжений в пункте наблюдений «Микижа» (52,99° с. ш., 158,22° в. д.). Оценки были получены ранее по аномалиям направленности геоакустической эмиссии. Показано, что гистограммы распределения направлений осей главных напряжений согласуются с результатами экспериментальных оценок для двух групп землетрясений. Модальные интервалы и локальные максимумы гистограмм распределения попадают в диапазоны оценок ориентации осей главных напряжений от 290° до 320° и от 20° до 50° соответственно.

The process of rock stress-strain state change causes acoustic radiation, which is called rock acoustic emission or geoacoustic emission. The relation between the earthquake preparation process and rock acoustic emission variations, which are called pre-seismic anomalies, has been stated in a number of investigations. The general mechanism of occurrences of these anomalies is associated with the fact that formation of a preparing earthquake source causes changes in the stress-strain state of rocks, surrounding it. One of the kinds of anomalies, occurring at the final stage of earthquake preparation, is the appearance of clearly expressed direction of acoustic activity. The main hypothesis of occurrence of this phenomenon is that a preparing earthquake source impact causes formation of constant direction of principle stress axes at an observation point. In their turn, the direction of these axes determines the primary orientation of acoustic radiation sources. To confirm this hypothesis, axis orientations of the main stresses, determined by the earthquake preparation process, were modeled. The estimates are based on the model constructed within the framework of elasticity linear theory where the Earth crust is considered in the form a homogeneous isotropic elastic half-space and the force impact at a preparing earthquake source is considered in the form of a combination of forces double pairs. Elastic deformation potential energy, accumulated during the earthquake preparation process, is taken into account. In the paper, we used the data from the catalog of earthquake source mechanics «The Global Centroid-Moment-Tensor Catalog» on seismic events occurred near Kamchatka peninsula from 1976 until 2020. As long as the acoustic radiation direction depends on the azimuthal direction to earthquake epicenter, all the considered seismic events were divided into three groups by the method of K-averages according to spatial locations of their epicenters. The modeling results were compared with experimental estimates of the main stress axis directions at Mikizha observation site (52.99° N, 158.22° E). The estimates were earlier obtained based on the geoacoustic emission directivity anomalies. It was shown that histograms of main stress axis direction distributions agree with the results of estimates for two groups of earthquakes. Modal intervals of distribution histograms fall within the range of experimental estimates from 290° to 320° and from 20° to 50° accordingly.