Acoustic Monitoring of Anomaly Stressed Zones, Determination Their Positions, Surfaces, Evaluation of Catastrophic Risk

Abstract Problem statement. Self-organization is not a universal property of matter, it exists under certain internal and external conditions and this is not associated with a special class of substances. The study of the morphology and dynamics of migration of anomalous zones associated with increased stresses is of particular importance in the development of deep deposits, complicated by dynamic phenomena in the form of rock impacts. Applied method and design: An important tool for this study is geophysical explora-tion. To describe the geological environment in the form of an array of rocks with its natural and technogenic heterogeneity, one should use its more adequate description, which is a discrete model of the medium in the form of a piecewise heterogeneous block medium with embedded heterogeneities of a lower rank than the block size. This nesting can be traced several times, i.e. changing the scale of the study; we see that heterogeneities of a lower rank now appear in the form of blocks for heterogeneities of the next rank. A simple averaging of the measured geophysical parameters can lead to distorted ideas about the structure of the medium and its evolution. Typical results: We have analyzed the morphology of the structural features of disintegration zones before a strong dynamic phenomenon. The introduction of the pro-posed integrated passive and active geophysical monitoring into the mining system, aimed at studying the transient processes of the redistribution of stress-strain and phase states, can help prevent catastrophic dynamic manife-stations during the development of deep-located deposits. Concluding note (Practical value/implications): Active geophysical monitoring methods hierarchical elastic structure located in the J-th layer of the N-layer elastic medium are constructed.

changing the scale of the study; we see that heterogeneities of a lower rank now appear in the form of blocks for heterogeneities of the next rank. A simple averaging of the measured geophysical parameters can lead to distorted ideas about the structure of the medium and its evolution. Typical results: We have analyzed the morphology of the structural features of disintegration zones before a strong dynamic phenomenon. The introduction of the proposed integrated passive and active geophysical monitoring into the mining system, aimed at studying the transient processes of the redistribution of stress-strain and phase states, can help prevent catastrophic dynamic manifestations during the development of deep-located deposits. Concluding note (Practical value/implications): Active geophysical monitoring methods should be tuned to a model of a hierarchical heterogeneous environment. Iterative algorithms for 2-D modeling and interpretation for sound diffraction and a linearly polarized transverse elastic wave on the inclusion with a

Introduction
The formation of structures during irreversible processes is associated with a qualitative leap when threshold (critical) parameters are reached. Self-organization is a supercritical phenomenon when the parameters of the system exceed their critical values. When a system deviates strongly from an equilibrium state, its variables satisfy non-linear equations. Nonlinearity is an important and general feature of processes that occur far from equilibrium. Moreover, the supercritical return of entropy is possible only if there is an unusual, special internal structure of the system [1]. This means that self-organization is not a universal property of matter, it exists under certain internal and external conditions and this is not associated with a special class of substances. So, there are two classes of irreversible processes: 1) destruction of the structure near the equilibrium position, this is a universal property of systems under arbitrary conditions; 2) the appearance of structures far from the equilibrium position under the conditions that the system is open and has nonlinear internal dynamics, and its external parameters have supercritical values. I. Prigogine called them dissipative structures [2]. The study of the morphology and dynamics of migration of these zones (named as disintegration zones) is of particular importance in the development of deep deposits, complicated by dynamic phenomena in the form of mountain impacts. An important tool for this study is geophysical exploration. As shown in [3], to describe the geological environment in the form of an array of rocks with its natural and technogenic heterogeneity, one should use its more adequate description, which is a discrete model of the medium in the form of a piecewise inhomogeneous block medium with embedded heterogeneities of lower rank than the block size. This nesting can be traced several times, i.e. changing the scale of the study; we see that heterogeneities of a lower rank now appear in the form of blocks for heterogeneities of the next rank. A simple averaging of the measured geophysical parameters can lead to distorted ideas about the structure of the medium and its evolution [3]. One of the fundamental mining problems, which are traditionally referred to the problems of geomechanics, is the development of theoretical and experimental methods for studying the structure and condition of rock masses in order to predict and prevent catastrophic phenomena during mining. This problem is compounded by the fact that the rock massif is under direct or indirect technogenic impact, which leads to significant unsteadiness of both the structure and the state of the mass [4]. When conducting mining operations in highly stressed rock massifs, technogenic seismicity is manifested, the forecasting and prevention issues of which receive a lot of attention in all countries with a developed mining industry. An important role here belongs to the short-term forecast; the methodology for selecting criteria for it is still a problem, both in mining and in seismology [5]. In the framework of the Institute of Mining of the Siberian Branch of Russian Academy of sciences, important results have been achieved on the study of the state of a rock massif in the framework of nonlinear geomechanics [6] using geophysical methods with a resolution to detect the nucleation and decay of self-organizing structures [7]. In [7], studies were conducted aimed at developing criteria for the spatial-temporal complex active and passive seismic and electromagnetic monitoring to prevent destructive dynamic phenomena based on six-year seismological monitoring data carried out by the mountain impact service at the Tashtagol underground mine and the experience gained use of the system of induction electromagnetic spatial-temporal monitoring developed at the Institute of geophysics of Ural Branch of Russian Academy of sciences for arrays of various material composition before and after mass explosions.

Theory, Research Methods
We analyzed the morphology of the structural features of disintegration zones before a strong dynamic phenomenon. During the next cycle of electromagnetic observations at the Tashtagol mine in August 2007 On August 9, a mountain impact occurred with an energy of log E = 6.9 in the pillar, located in the orth 3 alignment at a level of 16 m below the soil of horizon −280, Figure 1.
Three days before a mountain shock, in the orth 3 ( Figure 1, Figure 2), in the geoelectric sections of the soil, sub vertical discrete structures are found in which disintegration zones are combined. These structures appeared in the resonance mode at different frequencies and only at one frequency for each of the unit vectors. We discovered the same phenomenon earlier in one day at Estyuninskaya mine and SUBR mine 15 [8]. The appearance of these structures of sub vertical morphology is a precursor of a strong dynamic phenomenon; however, to determine the place and magnitude of an event, it is necessary to have information on the state of the unit arrays and their belonging to the corresponding ranks on the stability of the array, as was done in [9].
At present, theoretical results on modeling the electromagnetic and seismic fields in a layered medium with hierarchical inclusions are in demand. Modeling algorithms were constructed in the electromagnetic case for 3D heterogeneity, in the seismic case for 2D heterogeneity [10] [11]. It is shown that with an increase in the degree of hierarchy of the medium, the degree of spatial non-linearity of the distribution of components of the seismic and electromagnetic fields increases, which corresponds to the detailed monitoring experiments in shock hazardous mines of the Tashtagolsky mine and SUBR. The constructed theory  demonstrated how complicated the process of complexion methods using an electromagnetic and seismic field is to study the response of a medium with a hierarchical structure. This problem is linked with the formulation and solution of the inverse problem for the propagation of electromagnetic and seismic fields in such complex media. In [12] [13], the problem of constructing an algorithm for solving the inverse problem using the equation of the theoretical inverse problem for the 2D Helmholtz equation was considered. Explicit equations of the theoretical inverse problem are written out for the cases of scattering of an electromagnetic field (E and H polarization) and scattering of a linearly polarized elastic wave in a layered conducting and elastic medium with a hierarchical conducting or elastic inclusion, which are the basis for determining the contours of misaligned inclusions l-rank of hierarchical structure. Obviously, when solving the inverse problem, monitoring systems configured to study the hierarchical structure of the environment should be used as the initial monitoring data. On the other hand, the more complex the medium, the each wave field brings its own information about its internal structure; therefore, the interpretation of the seismic and electromagnetic fields must be carried out separately, without mixing these databases.

Two-Dimensional Sound Diffraction Modeling Abnormally Stressed Heterogeneity of a Hierarchical Type, Located in an N-Layer Elastic Medium
In [14], an algorithm for modeling sound diffraction by a two-dimensional elastic hierarchical inclusion located in the J-th layer of an N-layer medium is described ( ) wave number for a longitudinal wave, in the above expression, the index ji means belonging to the properties of the medium inside the heterogeneity, ja is outside the heterogeneity, λ is the Lame constant; σ is the density of the medium; ω is the circular frequency; gradϕ = u vector of displacements; φ 0 is the potential of a normal seismic field in a layered medium in the absence of heterogeneity: We assume that the density of the hierarchical inclusion for all ranks l and the enclosing layer are the same, and the elastic parameters of the hierarchical inclusion for all ranks differ from the elastic parameters of the enclosing medium, then the system of Equations (3) can be rewritten in the form:

Elastic Shear Wave Diffraction Modeling on an Abnormally Stressed Heterogeneity of a Hierarchical Type, Located in an N-Layer Elastic Medium
Similarly to (1), the same process is written for modeling the propagation of an elastic transverse wave in an N-layer medium with a two-dimensional hierar-chical structure of an arbitrary section morphology using the integral relations written in [15].
which coincides with the corresponding expression for the normal field in the work [15]. It should be noted that the structure of Equations (4) coincides with the general case when the hierarchical heterogeneity has not only elastic parameters different from the parameters of the enclosing medium, but also density parameters at all ranks differ from the density parameters of the enclosing layer. The difference between these problems lies only in the values of the wave number. Thus, the response of the medium associated with the longitudinal wave is more sensitive to the region of elastic heterogeneities in the array. This should be taken into account when assessing the status of a complex organized geological environment. In [13], the problem of constructing an algorithm for solving the inverse problem using the equation of the theoretical inverse problem for the 2-D Helmholtz equation was considered. An explicit equation of the theoretical inverse problem is obtained for the cases of scattering of a linearly polarized elastic wave in a layered elastic medium with a hierarchical elastic inclusion, the density of which for all ranks is equal to the density of the containing layer. An iterative algorithm for determining the contours of misaligned inclusions of the k-th rank in a hierarchical structure is constructed with the sequential use of the solution of the direct problem of calculating the elastic field of k-1 rank. With an increase in the degree of hierarchy of the structure of the medium, the degree of spatial non-linearity of the distribution of the components of the seismic field increases, which implies the exclusion of linearization methods for creating interpretation methods. This problem is linked with the solution of the inverse problem for the propagation of the seismic field in such complex environments using explicit equations of the theoretical inverse problem. For the first time, an equation was written for determining the surface of an abnormally stressed inclusion in a hierarchical layered-block medium according to acoustic monitoring data. In practice, using this algorithm, according to acoustic monitoring data, we can localize the region of a possible source of a rock shock or an earthquake and estimate the degree of anomalous elastic stresses.

Discussion of Results and Conclusions
Comparing expressions (1) and (2) we can draw the following conclusions. When constructing a complex seismic-gravity model without taking into account the anomalous influence of the stress-strain state inside the inclusion, the analysis of the anomalous acoustic effect using the data on the propagation of a longitudinal wave shows that it is also more sensitive to the form of inclusion, compared with the acoustic effect on the propagation of a transverse wave. However, it follows from these expressions that the influence of the elastic parameters in the host medium in the seismic model cannot be neglected, and they, when interpreted, affect the values of the desired anomalous densities. If these values are used in constructing the gravitational density model, then these density values will not reflect the material composition of the analyzed medium. When constructing an anomalously stressed geomechanical model without taking into account the anomalous effect of density heterogeneities within the inclusion, an analysis of the anomalous acoustic effect using data on the propagation of a shear wave shows that it is also more sensitive to the form of inclusion, compared with the acoustic effect on the propagation of a longitudinal wave. However, it follows from these expressions that the influence of the density parameters in the host medium in the seismic model cannot be neglected, and when interpreted, they affect the values of the desired anomalous elastic parameters that cause the anomalous stress state. If these values are used in the construction of the geomechanical model, then these values of the elastic parameters will not reflect the stress state of the analyzed medium. It is shown that with an increase in the degree of hierarchy of the medium, the degree of spatial non-linearity of the distribution of components of the seismic and electromagnetic fields increases, which corresponds to the detailed monitoring experiments in shock-hazardous mines of the Tashtagolsky mine and SUBR. The constructed theory demonstrated how complicated the process of complexion methods using an electromagnetic and seismic field to study the response of a medium with a hierarchical structure. This problem is linked with the formulation and solution of the inverse problem for the propagation of electromagnetic and seismic fields in such complex media. In [13], the problem of constructing an algorithm for solving the inverse problem using the equation of the theoretical inverse problem for the 2D Helmholtz equation was considered. Using the theory of solving the inverse problem, one can trace the migration of zones of abnormal stresses, their increase or decrease due to cyclic explosive influences during mining of the array, and also assess the possible risk of high-energy dynamic phenomena in the array.