\documentclass[twoside]{article}
\usepackage[silent,pdf,hyper]{elxshort}
\usepackage{color}
\usepackage{hyperref}
\journalid{ES}{19}{4}{19}
\citnumber{}{2}
\paperid{2019ES000663}
\papertype{RAR}
\cpright{}{2019}
\lefthead{manevich et al.}
\righthead{geoecological aspects}
\received{27}{September}{2018}
\accepted{3}{April}{2019}
\published{16}{July}{2019}
\definecolor{brown}{RGB}{157,23,4}

\authaddr{A. I. Manevich,
Geophysical Center of RAS, 3 Molodezhnaya St., 119296 Moscow. (ai.manevich@yandex.ru)}


\begin{document}

\title{Geoecological aspects of stress-strain state modeling results of Leninsky coal deposit (Kuzbass, Russia)}

\author{1}{Manevich}{A. I.}{}{1}{2}{}
\author{2}{Kolikov}{K. S.}{}{2}{}{}
\author{3}{Egorova}{E. A.}{}{2}{}{}

\affil{1}{Geophysical Center RAS}{Moscow}{}{Russia}
\affil{2}{Mining Institute of NUST ``MISiS''}{Moscow}{}{Russia}



\keywd{Leninsky coal deposit}
\keywd{mine}
\keywd{Kuzbass}
\keywd{modeling}
\keywd{stress-strain state}
\keywd{coal mining}

\abstract{This paper describes the analysis of numerical simulation results of the stress-strain state of Leninsky coal deposit area (Kuzbass) before and after allotment mining of mining fields specified in this article.
       The analysis of the influence of coal mining on the stress-strain state of rocks is performed taking into account the fault tectonics of the area. The particularities of the influence of regional stress field changes on geoecological safety are considered. The geoecological consequences of the exploitation of coal seams of the Leninsky coal deposit were analyzed.}

\section{1. Introduction}

       Leninsky coal deposit geographically belongs to the North-Western part of the Kuznetsk coal basin, inclining to the Kuznetsk's Alatau. The Kuznetsk coal basin itself is a large synclinorium, which has the shape of an irregular quadrangle in the plan, the long axis of which is oriented in the NW--SE direction
[\itc{Nedra,} \reflink{Nedr67}{1967}].
 In geomorphological terms, the region is an intermountain basin. It is bounded on the east by the Kuznetsk Alatau, on the south by mountain Shoria, and on the south-west and west by Salair. Against the background of Altai-Sayan folded region, represents a set of mountain ranges, clearly expressed in relief in the form of mountain ranges
[\itc{Nedra,} \reflink{Nedr67}{1967}].
 Modern geodynamic activity is expressed in the growth of new elevations due to the reduction of depressions resulted from the general uplift of the territory
 [\itc{Panina and Zaitsev,} \reflink{Pani12}{2012}].
  This is confirmed by geomorphological data, by the results of re-leveling, by observations using global navigation satellite systems (GNSS), as well as by spontaneous gas emission in the territory of the Leninsky coal deposit (LCD)
[\itc{Batugina and Petukhov,} \reflink{Batu88}{1988};
\itc{Kartashev and Yesikov,} \reflink{Kart83}{1983};
\itc{Novikov et al.,} \reflink{Novi13}{2013};
\itc{Panina and Zaitsev,} \reflink{Pani12}{2012}].

       A large amount of geological, tectonic and seismological data about the state of this area has been accumulated over the long period of LCD operation
[\itc{Adushkin,} \reflink{Adus15}{2015};
\itc{Emanov,} \reflink{Eman09}{2009};
\itc{Lazarevich and Polyakov,} \reflink{Laza09}{2009};
\itc{Novikov et al.,} \reflink{Novi13}{2013};
\itc{Panina and Zaitsev,} \reflink{Pani12}{2012};
\itc{Yakovlev et al.,} \reflink{Yako13}{2013}].
Large-scale mining operations cause massive anthropogenic changes in the geological environment, which lead to a negative impact on the geological environment, water objects, the ground surface and the biosphere. In this regard, the research of the stress-strain state of the Kuzbass areas is of particular relevance. And in this paper additional data are presented that could be useful in the geoecologycal assessment of the area.


 \begin{sloppypar}
 This paper aims to analyze the stress-strain state of the area of the LCD. For this, the SSS of a fragment of the geological environment of the studied area was modeled. In our work, we used the map of fault tectonics of Central Kuzbass
[\itc{Novikov et al.,} \reflink{Novi13}{2013};
\itc{Yakovlev et al.,} \reflink{Yako13}{2013}]
with drawn planned mine fields
%[Database, 2018]
(Database of State Geological Maps, All-Russian Scientific Research Geological Institute named  A. P. Karpinsky of the  RAS, % Russian Academy of Sciences,
 Electronic resource, URL:
 \href{http://webmapget.vsegei.ru/index.html}{http://webmapget.vsegei.ru/index.html},
 reference date: 10.09.2018)
as a geometric model. To determine the tectonic component of stresses, we used the data on the direction of the axes of the main tectonic stresses. The data was based on the results of morphotectonic studies of the Kuznetsk basin
[\itc{Novikov et al.,} \reflink{Novi13}{2013};
\itc{Panina and Zaitsev,} \reflink{Pani12}{2012}].
\end{sloppypar}


       The modeling of the stress-strain state of LCD is made using the calculation method of stress-strain state of heterogeneous block arrays with system of tectonic faults. In this case the fault refers to the volumetric area of the dispersion of rocks with an elastic modulus significantly lower than in the surrounding rocks
[\itc{Morozov et al.,} \reflink{Moro12}{2012},  \reflink{Moro18}{2018}].

\section{2. Structural and Tectonic Model of Leninsky Coal Deposit}

       From geological point of view, LCD is located in complicated tectonic conditions. The coal seams of this deposit are formed by Leninskaya syncline, which is a wide fold with asymmetric wings and mildly sloping downfold. The syncline downfold is sinking from the Central part of the field to the North-West (2--3\deg) and South-East (3--5\deg). The South-Western wing of the fold is complicated by Leninskaya anticline with the angle of inclination of wings varying from 20 to 40\deg \
[\itc{Popov,} \reflink{Popo76}{1976}].

       The area of Leninsky coal deposit, as per the results of morphotectonic studies
[\itc{Novikov et al.,} \reflink{Novi13}{2013}]
refers to the boundary between two neotectonic blocks of the Prisalair and Kemerovo sub-districts. Usually, there is no exposure of the latest disturbance zones, since they are covered by the deluvial aprons of piedmont foothills in case of differentiated movements along them or by alluvial sediments. This is only possible if there are no significant vertical movements
[\itc{Novikov et al.,} \reflink{Novi13}{2013}].
The thickness of covering complex is from 20 to 40~m
[\itc{Novikov et al.,} \reflink{Novi13}{2013}].
The heights of neighbour neotectonic blocks within the modelling area are comparable and differ by 20--30~m
[\itc{Novikov et al.,} \reflink{Novi13}{2013};
\itc{Panina and Zaitsev,} \reflink{Pani12}{2012}].


\begin{figure}[t]                     % Fig 1
\figurewidth{20pc}
\setimage{}{}{20pc}{}{2019es000663-f01}
\bottomcaption{Scheme of fault tectonics of the Leninsky coal deposit based on data of
[\itc{Lazarevich and Polyakov,} \reflink{Laza09}{2009};
\itc{Novikov et al.,} \reflink{Novi13}{2013};
\itc{Yakovlev et al.,} \reflink{Yako13}{2013}]
and borders of main mine fields
(Database of State Geological Maps, All-Russian Scientific Research Geological Institute named  A. P. Karpinsky of the  RAS, % Russian Academy of Sciences,
 Electronic resource, URL:
 \href{http://webmapget.vsegei.ru/index.html}{http://webmapget.vsegei.ru/index.html},
 reference date: 10.09.2018):
1 -- tectonic faults;
2 -- agglomeration of Leninsk-Kuznetsky city;
3 -- Leninsk-Kuznetsky city;
4 -- mine fields of LCD scheduled to mining, which refer to the area of the supposed dispersion of rocks.}
\vspace*{\baselineskip}
\end{figure}


       Thus, the structural-tectonic scheme of tectonic fault was used to simulate the stress-strain state (SSS) of the Leninsky coal deposit area
(\figref{1}),           % (Fig. 1),
which was verified by geophysical studies, the results of which are presented in
[\itc{Lazarevich and Polyakov,} \reflink{Laza09}{2009};
\itc{Novikov et al.,} \reflink{Novi13}{2013};
\itc{Yakovlev et al.,} \reflink{Yako13}{2013}].
The SSS of a fragment of the geological environment of $150\times 140$~km was simulated. The depth of the simulated layer corresponds to the depth of the mine fields planned for the development (or already developed) and is within 500--1000~m from the ground surface. The simulated layer is composed of elastic and fragile rocks: coal, mudstone, siltstone and sandstone. Within the framework of the elastic formulation of the problem, the model of SSS of an elastic isotropic layer with Young's modulus
$E = 8\times 10^3$~MPa disturbed by a system of tectonic faults with Young's modulus lower by two orders and the total Poisson's ratio
$\mu = 0.25$ is simulated
%[Coal, 2003;
[\itc{Geoinformark,} \reflink{Geoi03}{2003};
\itc{Morozov et al.,} \reflink{Moro18}{2018}].
 We assumed physical and mechanical parameters of the model tentatively, based on the fact that the simulation results are aimed to obtain qualitative rather than quantitative estimates of SSS of the area.


\fitopen{1.1}
       The fault in the model corresponds to mechanically dispersed environment with the width of the fracture zone up to 500 m
 [\itc{Morozov et al.,} \reflink{Moro12}{2012}].
 The dip of fault was assumed as vertical at interval depth of intensive mining works.The external field of tectonic stresses is specified by the dominant compression stress $\sigma_{yy}$ with the submeridional orientation axis, which is consistent with the materials of morphotectonic
  [ \itc{Novikov et al.,} \reflink{Novi13}{2013};
  \itc{Panina and Zaitsev,} \reflink{Pani12}{2012}]
  and seismological
 [\itc{Emanov,} \reflink{Eman09}{2009};
 \itc{Ovsyuchenko et al.,} \reflink{Ovsy10}{2010};
 \itc{Zhalkovsky et al.,} \reflink{Zhal95}{1995}]
 research.


 \begin{sloppypar}
 The model of the block diagram after mining of fields is as follows. The boundaries of the main mine fields planned for mining and located in the range of 500--1000~m from the ground surface
 %[Database, 2018]
 (Database of State Geological Maps, All-Russian Scientific Research Geological Institute named  A.~P.~Karpinsky of the  RAS, % Russian Academy of Sciences,
 Electronic resource, URL:
 \href{http://webmapget.vsegei.ru/index.html}{http://webmapget.vsegei.ru/index.html},
 reference date: 10.09.2018)
 were drawn on the scheme of fault tectonics. Excavation, drive working and all large-scale underground mining within the boundaries of mine fields in general form the volume of dispersed rocks. First of all, this is due to the technology of coal face operations. LCD used the rock pressure control by the complete collapse of the roof of coal seams. Excavation of coal seams at different depths with the collapse of the roof rocks causes the appearance of large zones, with reduced physical and mechanical properties in comparison with the host rocks
 [\itc{Batugin et al.,} \reflink{Batu16}{2016};                      %%%
 \itc{Bykadorov et al.,} \reflink{Byka18}{2018};         % ???
 \itc{Tarasenko et al.,} \reflink{Tara17}{2017}].
 This is due to the physical appearance of voids and the filling of these voids with water.
 \end{sloppypar}


       We analyzed fields of the stress tensor components $\sigma_{xx}$, $\sigma_{yy}$, $\tau_{xy}$ and fields of the stress intensity $\sigma_i$.  Region includes mine fields of Leninsky coal deposit and agglomeration of % Leninsk-Kuzneckiy
 Leninsk-Kuznetsky
 city before and after mining of the fields.
\fitclose

\section{3. Method of Stress-Strain State Modeling}

	In this study we used the methodology and software to calculate stress-strain state of block heterogeneous massifs disturbed by a system of tectonic faults in the external field of tectonic stresses
[\itc{Morozov et al.,} \reflink{Moro12}{2012},  \reflink{Moro18}{2018}].
The medium inside the study area is isotropic-elastic, disturbed by the system of tectonic faults and zones of mined out areas.

 \begin{sloppypar}
        In modelling process, the fault and zones of worked-out coal seams correspond to isotropic-elastic medium with a modulus of elasticity twice lower than elastic modulus of hosting rocks.
 \end{sloppypar}

       The stress-strain state in the massif is calculated by the finite element method as displacements based on four-square isoparametric elements.
[\itc{Morozov et al.,} \reflink{Moro12}{2012}].
The choice of the finite element method as a calculation approach stems from physical considerations that yield conditions for minimizing the total potential strain energy through-out the massif; the efficiency of mathematical calculations is also a factor.

       The relations between stresses and strains are taken according to the generalized Hooke elasticity law
[\itc{Baklashov,} \reflink{Bakl04}{2004}]:

\begin{eqnarray*}
\sigma_{ij} = D_{ij} \times \varepsilon_{ij}
\end{eqnarray*}
where $\sigma_{ij}$  are the components  of the stress tensor:

\begin{eqnarray*}
\sigma_{ij} = \left\{
 \begin{array}{cc}
\sigma_{xx} &\tau_{xy} \\
\tau_{xy} & \sigma_{yy}
\end{array}
\right\}
\end{eqnarray*}
$\varepsilon_{ij}$ are the components of elastic strain tensor:

\begin{eqnarray*}
\varepsilon_{ij} = \left\{
 \begin{array}{cc}
\varepsilon_{xx} &0.5 \varepsilon_{xy} \\
0.5 \varepsilon_{xy} & \varepsilon_{yy}
\end{array}
\right\}
\end{eqnarray*}
$D(\mu, E)$  is the matrix of elastic coefficients for isotropic-elastic material. Where $\mu$ is Poisson's ratio; $E$ is Young's modulus,~MPa.

       The stress intensity is calculated by the following formula
\eqref{1}
       [\itc{Baklashov,} \reflink{Bakl04}{2004}]:
 \begin{equation}\label{1}
\sigma_{i} =
(\sigma_{xx}^2 +\sigma_{yy}^2 - \sigma_{xx} \times \sigma_{yy} +
3 \times \sigma_{xy}^2)^{1/2}
\end{equation}                     % (4)

       The change of SSS of the studied area is analysed by the magnitude of difference of stress intensity before and after mining of fields:
\begin{equation}\label{2}
\Delta \sigma_{i} =
\langle\sigma_{i} \rangle_I -\langle\sigma_{i} \rangle_{II}
\end{equation}                     % (5)
where $\langle\sigma_{i} \rangle_I$ and  $\langle\sigma_{i} \rangle_{II}$ is the intensity of stresses at an arbitrary point before and after mining, respectively.

\begin{figure}[t]                     % Fig 2
\figurewidth{20pc}
\setimage{}{}{20pc}{}{2019es000663-f02}
\caption{Map of stress intensity $\sigma_i$ prior to the mining of fields:
1 -- tectonic faults;
2~--~agglomeration of Leninsk-Kuznetsky city;
3 -- Leninsk-Kuznetsky city.}
\end{figure}


\begin{figure}[t]                     % Fig 3
\figurewidth{20pc}
\setimage{}{}{20pc}{}{2019es000663-f03}
\caption{The map of shear stress $\tau_{xy}$ prior to the mining of fields:
1 -- tectonic faults;
2  -- agglomeration of Leninsk-Kuznetsky city;
3  -- Leninsk-Kuznetsky city.}
\end{figure}


\begin{figure}[t]                     % Fig 4
\figurewidth{20pc}
\setimage{}{}{20pc}{}{2019es000663-f04}
\caption{Map of stress intensity $\sigma_i$ after mining of fields:
1 -- tectonic faults;
2  -- agglomeration of Leninsk-Kuznetsky city;
3 -- Leninsk-Kuznetsky city;
4 -- mine fields of LCD scheduled to mining which refer to the area of the supposed dispersion of rocks.}
\end{figure}

       To model the geological massif SSS in the coordinate system $y, x $ (coinciding with the direction North--South and West--East), we establish the conditions of the external field of tectonic stresses. We used stress values $\sigma_{xx}= 10$~MPa and
$\sigma_{yy} = 30$~MPa based on the different results of measurements of the main stresses of the Earth's crust in various regions
[\itc{Selin,} \reflink{Seli08}{2008};
 \itc{Zubkov,} \reflink{Zubk16}{2016}].   %%% ??? +
 The selected ratio includes the horizontal tectonic component and the lateral repulse stress of vertical stresses $\sigma_{zz}$.

\section{4. Results of Numerical Simulation of SSS of LCD}


 \begin{sloppypar}
 \figref{2}          % Figure 2
represents the map of the LCD stress intensity before the mining works. Background values $\sigma_i$ in initial massif are in the range of 10--20~MPa. Near the agglomeration of Leninsk-Kuznetsky city and Leninsky deposit, the zone of increased stress intensity $> 35$~MPa elongated in the NNE orientation is distinguished. One of the coal mine fields belong to this zone.
\end{sloppypar}



\figref{3}          %  Figure 3
shows map of the shear stress ($\tau_{xy}$). The studied area contains the zone of increased shear stresses. The size and orientation of this zone corresponds to the zone of increased stress intensity. The background values of the shear stresses $\tau_{xy}$ in the study area are in the range of
$\pm 5$~MPa
	
       At the next iteration, the SSS of the studied area was calculated for the case when all the planned mine fields in the simulated depth range of 500--100~m are worked out.
\figref{4}          % Figures 4
and
\figref{5}          % 5
show maps of stress intensity and shear stresses.

%\def\bottomfraction{.8}
%\def\textfraction{.1}

\begin{figure}[t]                     % Fig 5
\figurewidth{20pc}
\setimage{}{}{20pc}{}{2019es000663-f05}
\caption{The map of shear stress $\tau_{xy}$ after mining of fields:
1 -- tectonic faults;
2 -- the agglomeration of Leninsk-Kuznetsky city;
3 -- Leninsk-Kuznetsky city;
4 -- mine fields of LCD scheduled to mining which refer to the area of the supposed dispersion of rocks.}
\end{figure}

       After the end of coal face operations in the area of city agglomeration of Leninsk-Kuznetsky and Leninsky deposit, the area with increased stress
(\figref{4}
and
\figref{5})          %          % (Fig. 4 and 5)
is  significantly expanded in the North-West and South-East direction, covering an area of about 2.5 times more before the mining works in the fields. Also, a new area with increased value of stress intensity and shear stress at the end of the mine field in the South-East direction is being formed; this area is currently being worked by two mines: mine named after Kirov and mine ``Polysaevskaya''.

\figref{6}          %        Figure 6
shows the map of the stress intensity difference $\Delta \sigma_i$ calculated by using the
formula
\eqref{2}.           % (5).


	The range of released stress intensity 0--5~MPa occupies a significant area and refers to the background. The maximum values of the stress intensity difference $\Delta \sigma_i$  reach 15--20~MPa. Mining works form areas with high stress intensity which refer to zones of increased energy saturation.

	It should be noted that there is uncertainty in the interpretation of the simulation results. During the long period of deposit coal mining the stress relaxation may occur. This relaxation is caused by the induced seismic phenomena and stress redistribution as a result of slow movements.

\section{5. Geoecological Aspects of Simulation Results}
	
       Analysis of the results of SSS modelling based on the system approach
[\itc{Gvishiani et al.,} \reflink{Gvis18}{2018}]
showed that prior to the mining works in the area there was a local zone with heterogeneity of stress fields in the massif. This zone is located between the endings of two tectonic structures limiting the Leninskaya anticline from the Eastern side and oriented to the NNE direction. Then, during the process of development of mining and working-out of coal-bearing seams of LCD, new zones of local heterogeneity of stress field are formed. These zones are confined to the places of dispersion of the geological environment arising from the mining of coal seams. The area of LCD, at the time of mining, includes three zones of significantly changed stress intensity. The regions of negative values limiting the region from the West and the East indicate the increased energy saturation of the hosting medium, and the region with positive values $\Delta \sigma_i$  indicates the reduced energy saturation. The alternation of such zones at one geological formation can induce the realization of the energy of potential deformations, and serve as a trigger effect for the dynamic events of the host rock mass
[\itc{Adushkin,} \reflink{Adus15}{2015};
\itc{Emanov et al.,} \reflink{Eman18}{2018};
\itc{Petukhov and Batugina,} \reflink{Petu99}{1999}].

       According to the latest data, the problems of safety assessment of mining operations in underground coal mining do not lose their relevance
[\itc{Batugin,} \reflink{Batu18}{2018};
\itc{Litvinov et al.,} \reflink{Litv17}{2017}].
The processes of methane release and the formation of conditions for the occurrence of gas-dynamic phenomena in the coal massif is caused by several factors
[\itc{Kolikov et al.,} \reflink{Koli16}{2016};
\itc{Litvinov et al.,} \reflink{Litv17}{2017}]:
a) the increase of the depth of mining operations and the rate of coal seams mining (the rate of advance of the working face), the complication of mining and geological conditions significantly affect the processes of methane release and the formation of conditions for the occurrence of gas-dynamic phenomena in the coal massif. The developed models of the stress-strain state of the LCD area can be used as boundary conditions for bigger scale models. A small-scale study of the SSS of the coal massif is necessary for qualitative assessment of changes in the relative permeability of the coal massif, and further search for potentially dangerous zones of formation of gas-dynamic phenomena
[\itc{Kolikov et al.,} \reflink{Koli17}{2017};
 \itc{Petukhov and Batugina,} \reflink{Petu99}{1999}].

\begin{figure}[t]                     % Fig 6
\figurewidth{20pc}
\setimage{}{}{20pc}{}{2019es000663-f06}
\caption{The difference between the stress intensity $ \Delta \sigma_i$ before and after mining of fields:
1 -- tectonic faults;
2 -- the agglomeration of Leninsk-Kuznetsky city;
3 -- Leninsk-Kuznetsky city.}
\end{figure}

       The stress state can have a significant impact on geo-ecological safety at the liquidation of mines and reclamation of district. In the study
 [\itc{Morozov et al.,} \reflink{Moro12}{2012}]
 it was shown that the release of stress contributes to the flow of water fluids from the more intense surrounding massifs into the weakened zone, thereby forming a horizontal gradient of the filtration. Leninsky coal deposit is crossed almost in the middle by the river Inya, which belongs to the III category under the conditions of safety of mining operations under water objects. The third group includes water objects that intersect in the mined areas with geological disturbances, or where the working occurs
[\itc{Gosgortechnadzor,} \reflink{Gosg98}{1998}].
Inya River is located in the area of Leninskaya syncline. The safety depth of mining under this river is 200~m, according to the calculation for the objects of category III
[\itc{Gosgortechnadzor,} \reflink{Gosg98}{1998}].

       The change of the hydrological regime is an environmental hazard factor. The change in the level of groundwater, the creation of depression craters or dewatering in rivers, lakes and other open water bodies can lead to flooding of the ground surface, to settling of the ground surface and to climate change. The mechanism activating the displacement of the Earth surface after the end of mining works (in the liquidation step) has also another nature. It is due to the decrease in the strength properties of wet rocks during the flooding of worked out area by mine water  [\itc{Batugin and Klimanova,} \reflink{Batu00}{2000}].

       During the coal seams excavation the technology of roof fall is used. The filling of  worked out area by water may occur due to the drainage of the water horizons entering into the zone of water-conducting fractures. Therefore, the large inflows are possible not only under the river or floodplain, but also at a considerable distance from them. In a number of cases mentioned in
 [\itc{Gvirtsman et al.,} \reflink{Gvir77}{1977}],
 there was no increase in the inflow of water into the mine workings and there were no records about formation of cracks on the ground surface extending under the waterbody. It follows that with sufficient depth of mining, the cracks that are forming on the ground surface, do not communicate with the zone of water-conducted cracks over the worked-out space.
	
       However, the presence of hydraulic connection of the Inya River and some layers of LCD can be estimated. The following cases of increases of inflows into the mined-out space were recorded at undermining of flood-lands of Inya River
[\itc{Gvirtsman et al.,} \reflink{Gvir77}{1977}]:
  \begin{itemize}
 \item
 Boldyrevsky seam, the depth of the working is 150~m, the increase of the inflow was 40~m$^3$/h;
\item
 Tomachevsky seam, the depth of the working is 20--120~m, the increase in the inflow was 100~m$^3$/h;
\item
 Polenov seam, the depth of the working is 120--200~m, the increase in the inflow have not been observed.
 \end{itemize}

\section{6. Summary}
	
       The development of mining operations at the mining of LCD coal seams alters geological and geomechanical conditions in the host rocks, leads to the activation of negative and dangerous geodynamic processes and phenomena: induced earthquakes, rockbursts and tectonic strokes, displacement of the ground surface, creation of new surfaces of weakening in the geological massif, change of physical and mechanical properties of rocks, etc.

       This is due to several interrelated and mutually influencing factors: a) the movement of a large volume of rocks from the geological environment to the ground surface (as a consequence, the physical appearance of a large amount of worked-out space); b) changes in the stress-strain state of the LCD area; c) changes in the hydrological regime of the area; d) changes in the physical and mechanical properties of the geological environment.

       The analysis of the results of modelling showed the presence of local zones with elevated stresses in the geological massif. These stresses relate to the tectonic structures and areas of mining of coal seams. The area of the Leninsky coal deposit at the time of completion of mining works potentially includes three areas of released stress intensity, the presence of which must be taken into account when assessing the safety of mining operations. Also, given the presence of hydraulic connection of Inya River with the deposit it is possible to judge about the change of the hydrological mode of the area on the basis of stress-strain state modelling. All these aspects should be taken into account in operation, liquidation and reclamation of the LCD area.

       In conclusion, it should be noted that the obtained values of the SSS characteristics should be considered only as a qualitative illustration. Despite this, the results of numerical modelling of the SSS of the Leninsky coal deposit can be used as data for qualitative assessment of the impact of mine operation on the host geological environment for the conditions of Kuzbass. The obtained data could be useful in the formulation of geological and geophysical observations and geodynamic studies of the area.

\acknowledgment{This work was conducted in the framework of budgetary funding of GC RAS, adopted by the Ministry of Science and Higher Education of the Russian Federation.}

\begin{references}

\reference{Adus15}
\reftype{5}
\refauthf{Adushkin}{V. V.}{}
\refyear{2015}
\reftitle{Trigger seismicity of the Kuznetsk Basin, Trigger effects in geosystems.
Theses of the reports of the III All-Russian seminar-meeting}
\refpubl{Institute of Dynamics of Geospheres RAS}
\reflocat{Moscow}
\refpages{8--28}
\endreference{}

\reference{Bakl04}
\reftype{3}
\refauthf{Baklashov}{I. V.}{}
\refyear{2004}
\refbook{Geomechanics. Volume 1. Fundamentals of Geomechanics}
\refpubl{MSMU}
\reflocat{Moscow}
\refpages{208}
\endreference{}

\reference{Batu18}
\reftype{2}
\refauthf{Batugin}{A.}{}
\refyear{2018}
\reftitle{Critically Stressed Areas of Earth's Crust as Medium for Man-Caused Hazards}
 \refjour{E3S Web Conf., VII International Scientific Conference ``Problems of Complex Development of Georesources''}
\refvol{56}
\refpages{02007}
\refdoi{10.1051/e3sconf/20185602007}
\endreference{}

\reference{Batu00}
\reftype{1}
\refauthf{Batugin}{A. S.}{}
\refauth{Klimanova}{V. G.}{}
\refyear{2000}
\reftitle{Flooding of the liquidated mines as a possible cause of technogenic earthquakes in mining regions}
 \refjour{Mining Information and Analytical Bulletin} % (scientific and technical journal)}
\refvol{2}
\refpages{58--61}
\endreference{}

\reference{Batu16}
\reftype{1}
\refauthf{Batugin}{A. S.}{}
\refgroup{et al.}
\refyear{2016}
\reftitle{Tectonophysical model of fault tectonic rock burst with wing sliding Journal of Liaoning Technical University}
\refjour{Natural Science}
\refvol{35}
\refpages{561--565}
\endreference{}


\reference{Batu88}
\reftype{3}
\refauthf{Batugina}{I. M.}{}
\refauth{Petukhov}{I. M.}{}
\refyear{1988}
\refbook{Geodynamic Zoning of Deposits in the Design and Operation of Mines}
\refpubl{Nedra}
\reflocat{Moscow}
\refpages{166}
\endreference{}


\reference{Byka18}   %%% ???
\reftype{2}
\refauthf{Bykadorov}{A. I.}{}
\refauth{Chernukha}{A. V.}{}
\refauth{Svirko}{S. V.}{}
\refyear{2018}
\reftitle{Aspects of liquidation of mines of the Prokopyevsko-Kiselevsky coal deposit}
\refjour{Russian Coal Journal}
\refnumb{2}
\refpages{88--94}
\refdoi{10.18796/0041-5790-2018-2-88-94}
\endreference{}

\reference{Eman09}
\reftype{1}
\refauthf{Emanov}{ A. F.}{}
\refauth{Emanov}{ A. A. }{}
\refauth{Leskova}{ E. V. }{}
\refauth{Fateev}{A. V.}{}
\refauth{Semin}{A. Yu.}{}
\refyear{2009}
\reftitle{Seismic activation at coil mining in Kuzbass}
\refjour{Physical Mesomechanics}
\refvol{12}
\refnumb{1}
\refpages{37--43}
\endreference{}

\reference{Eman18}
\reftype{2}
\refauthf{Emanov}{ A. A.}{}
\refauth{Emanov}{ A. F. }{}
\refauth{Fateev}{ A. V. }{}
\refauth{Leskova}{E. V.}{}
\refyear{2018}
\reftitle{Simultaneous Impact of Open-Pit and Underground Mining on the Subsurface and Induced Seismicity}
\refjour{Seismic Instruments}
\refvol{54}
\refnumb{4}
\refpages{479--487}
\refdoi{10.3103/S0747923918040035}
\endreference{}


\reference{Geoi03}
\reftype{3}
\refauthf{Geoinformark}{}{}
\refyear{2003}
\refbook{Coal Base of Russia,
Volume 2. Coal Basins and Deposits of West Siberia (Kuznetsky, Gorlovsky, West-Siberian basins, the Deposits of the Altai Territory and the Altai Republic)}
\refpubl{Geoinformark}
\reflocat{Moscow}
\refpages{604}
\endreference{}

\reference{Gosg98}
\reftype{3}
\refauthf{Gosgortechnadzor}{}{}
\refyear{1998}
\refbook{Resolution of the Gosgortechnadzor of Russia of 16.03.1998  \#13 ``On the approval of the rules for the protection of structures and natural objects from the harmful influence of underground mining in coal deposits''}
\refpubl{Gosgortechnadzor}
\reflocat{Russia}
\endreference{}


\reference{Gvir77}
\reftype{3}
\refauthf{Gvirtsman}{B. Ya.}{}
\refauth{Katsnelson}{N. N.}{}
\refauth{Boshenyatov}{E. V.}{}
\refyear{1977}
\refbook{Safe Excavation of Coal Under Water Objects}
\refpubl{Nedra}
\reflocat{Moscow}
\refpages{175}
\endreference{}

\reference{Gvis18}
\reftype{4}
\refauthf{Gvishiani}{A. D.}{}
\refauth{Weisberg}{L. A.}{}
\refauth{Tatarinov}{V. N.}{}
\refauth{Manevich}{A. I.}{}
\refyear{2018}
\reftitle{System analysis in mining sciences and decreasing environment damage}
\refbook{Proceedings of the conference dedicated to the memory of Academician A. V. Kryazhimsky ``System Analysis: Modeling and Management''}
\refpubl{Steklov Institute of Mathematics RAS}
\reflocat{Moscow}
\refpages{43--45}
\refdoi{10.4213/proc20588}
\endreference{}

\reference{Kart83}
\reftype{4}
\refauthf{Kartashev}{M. P.}{}
\refauth{Yesikov}{N. P.}{}
\refyear{1983}
\reftitle{The study of deformations of the earth's surface on the Kuzbass polygon}
\refbook{Results of researches in the international geophysical projects ``Modern Movements and Deformations of the Earth's Crust on Geodynamic Polygons''}
\refpubl{Nauka}
\reflocat{Moscow}
\refpages{102--105}
\endreference{}

\reference{Koli17}
\reftype{2}
\refauthf{Kolikov}{K.\! S.}{}
\refauth{Egorova}{E.\! A.}{}
\refauth{Mazanik}{E. V.}{}
\refyear{2017}
\reftitle{Geoinformation support of the forecast of dangerous natural processes during the mining of coal seams}
\refjour{Geoinformatics Research Papers: Proceedings of the Geophisical Center RAS}
\refvol{5}
\refnumb{2}
\refpages{3--12}
\refdoi{10.2205/2017BS045}
\endreference{}

\reference{Koli16}
\reftype{1}
\refauthf{Kolikov}{K. S.}{}
\refauth{Nikitin}{S. G.}{}
\refauth{Manevich}{A. I.}{}
\refyear{2016}
\reftitle{Analytical evaluation of the forecast methan emission recommended safety regulations}
\refjour{Safety of Labor in Industry}
\refnumb{8}
\refpages{34--39}
\endreference{}

\reference{Laza09}
\reftype{4}
\refauthf{Lazarevich}{T. I.}{}
\refauth{Polyakov}{A. N.}{}
\refyear{2009}
\reftitle{Mining monitoring of seismic and geodynamic safety of Kuzbass}
 \refbook{Mining Geomechanics and Mine Surveying: A Collection of Scientific Papers}
\refpubl{VNIMI}
\reflocat{SPb}
\refpages{252}
\endreference{}

\reference{Litv17}
\reftype{2}
\refauthf{Litvinov}{A. R.}{}
\refauth{Kolikov}{K. S.}{}
\refauth{Ishhneli}{O. G.}{}
\refyear{2017}
\reftitle{Accident and injury at coal industry enterprises in 2010--2015}
\refjour{Bulletin of the Scientific Center for the Safety of Work in the Coal Industry}
\refnumb{2}
\refpages{6--17}
\refdoi{10.26631/arc2-2017-6-17}
\endreference{}

\enlargethispage{-18\baselineskip}

\reference{Moro12}
\reftype{2}
\refauthf{Morozov}{V. N.}{}
\refauth{Kolesnikov}{I. Yu.}{}
\refauth{Tatarinov}{V. N.}{}
\refyear{2012}
\reftitle{Modeling the Hazard Levels of Stress-Strain State in Structural Blocks in Nizhnekanskii Granitoid Massif for Selecting Nuclear Waste Disposal Sites}
\refjour{Water Resources}
\refvol{39}
\refnumb{7}
\refpages{756--769}
\refdoi{10.1134/S009780781207007X}
\endreference{}

\reference{Moro18}
\reftype{2}
\refauthf{Morozov}{V. N.}{}
\refauth{Tatarinov}{V. N.}{}
\refauth{Kolesnikov}{I. Yu.}{}
\refauth{Manevich}{A. I.}{}
\refyear{2018}
\reftitle{Modeling the Stress-Strain State in the Epicentral Zone of a Strong Earthquake in Iran (December 26, 2003, $Mw = 6.6$)}
\refjour{Izvestiya, Physics of the Solid Earth}
\refvol{54}
\refnumb{4}
\refpages{602--611}
\refdoi{10.1134/S1069351318040080}
\endreference{}

\reference{Nedr67}
\reftype{3}
\refauthf{Nedra}{}{}
%Geology of the USSR
\refyear{1967}
\refbook{Geology of the USSR, Volume 14. Geological Description. Western Siberia}
\refpubl{Nedra}
\reflocat{Moscow}
\refpages{664}
\endreference{}

\reference{Novi13}
\reftype{2}
\refauthf{Novikov}{ I. S.}{}
\refauth{Cherkas}{ O. V. }{}
\refauth{Mamedov}{ G. M. }{}
\refgroup{et al.}  %%\refauth{Simonov}{Yu. G.}{} \refauth{Simonova}{T.\! Yu.}{} \refauth{Nastavko}{V.\! G.}{}
\refyear{2013}
\reftitle{Stages of activation and tectonic divisibility of the Kuznetsk coal basin (Southern Siberia)}
\refjour{Geology and Geophysics}
\refvol{54}
\refnumb{3}
\refpages{424--437}
\refdoi{10.1016/j.rgg.2013.02.007}
\endreference{}

\reference{Ovsy10}
\reftype{1}
\refauthf{Ovsyuchenko}{A. N.}{}
\refauth{Rogozhin}{E. A.}{}
\refauth{Novikov}{S. S.}{}
\refgroup{et al.}  %%\refauth{Marakhanov}{A. V.}{} \refauth{Larkov}{A. S.}{} \refauth{Akbiev}{R. T.}{} \refauth{Mogushkov}{I. M.}{}
\refyear{2010}
\reftitle{Paleogeological and tectonic studies of earthquake-prone areas in the south of the Kuznetsk Basin, Seismic resistant construction}
\refjour{Safety of Structures}
\refnumb{6}
\refpages{35--45}
\endreference{}

\reference{Pani12}
\reftype{2}
\refauthf{Panina}{L. V.}{}
\refauth{Zaitsev}{V. A.}{}
\refyear{2012}
\reftitle{Neotectonics and geodynamics of the Kuznetsk cavity}
\refjour{Bulletin of the Moscow University Series 4. Geology}
\refnumb{6}
\refpages{13--20}
\refdoi{10.3103/S014587521206004X}
\endreference{}

\reference{Petu99}
\reftype{3}
\refauthf{Petukhov}{I. M.}{}
\refauth{Batugina}{I. M.}{}
\refyear{1999}
\refbook{Geodynamics of Soils (Second edition, revised and enlarged)}
\refpubl{Nedra Communications LTL}
\reflocat{Moscow}
\refpages{287}
\endreference{}

\reference{Popo76}
\reftype{1}
\refauthf{Popov}{Yu. N.}{}
\refyear{1976}
 \reftitle{Some features of the tectonics of mine fields in the North-Eastern part of the Leninsky district of Kuzbass}
\refjour{Izvestiya of the Tomsk Polytechnic Institute}
\refvol{236}
\refpages{119--128}
\endreference{}


\reference{Seli08}
\reftype{1}
\refauthf{Selin}{K. V.}{}
\refyear{2008}
\reftitle{Measurement of the initial stresses of the array with depth in the regions of the world}
\refjour{Mining Information and Analytical Bulletin}
\refvol{10}
\refpages{297--302}
\endreference{}

\reference{Tara17}
\reftype{2}
\refauthf{Tarasenko}{ I. A.}{}
\refauth{Zinkov}{ A. V. }{}
\refauth{Chudaev}{ O. V. }{}
\refauth{Vetoshkina}{A. V.}{}
\refauth{Holodilov}{I. I.}{}
\refyear{2017}
 \reftitle{Technogenic effect of liquidation of coal mines on earth's entrails: hydrogeochemical aspect}
\refjour{IOP Conference Series: Earth and Environmental Science}
\refvol{87}
\refnumb{4}
\refdoi{10.1088/1755-1315/87/4/042022}
\endreference{}

\reference{Yako13}
\reftype{2}
\refauthf{Yakovlev}{D.\! V.}{}
\refauth{Lazarevich}{T.\! I.}{}
\refauth{Tsirel'}{S.\! V.}{}
\refyear{2013}
 \reftitle{Natural and induced seismic activity in Kuzbass}
\refjour{Journal of Mining Science}
\refvol{49}
\refnumb{6}
\refpages{862--872}
\refdoi{10.1134/S1062739149060038}
\endreference{}

\reference{Zhal95}
\reftype{1}
\refauthf{Zhalkovsky}{N. D.}{}
\refauth{Kuchai}{O. A.}{}
\refauth{Muchnaya}{V. I.}{}
\refyear{1995}
\reftitle{Seismicity and some characteristics of the stress state of the earth's crust of the Altai-Sayan region}
\refjour{Geology and Geophysics}
\refvol{36}
\refnumb{10}
\refpages{20--30}
\endreference{}

\reference{Zubk16}
\reftype{1}
\refauthf{Zubkov}{A. V.}{}
\refyear{2016}
\reftitle{The low of natural stress formation of the Earth's crust}
\refjour{Lithosphere}
\refnumb{5}
\refpages{146--151}
\endreference{}

\vspace*{1ex}

\end{references}

\end{document}



