Forecasting Underground Water Dynamics within the Technogenic Environment of a Mine Field. Case Study

The objective is to analyze the dynamics of the underground water of a mine field based on the study of the geofiltration process of the rock mass disturbed by mining to achieve safe extraction operations as well as subsurface territories at the stage of the mining enterprise closure. Numerical modeling, based on a finite difference method under the conditions of multifactority and definite uncertainty of processes of transformation of technogenic environment of a mine field, helps solve a problem concerning underground water dynamics forecasting. A hydrodynamic model of the M.I. Stashkov mine was developed while solving option series of epignosis problems in terms of the chronology of mine field stoping. The abovementioned made it possible to identify regularities of the history of filtration, the capacity parameters of rock mass and the expansion of areas of heightened hydraulic conductivity as well as to evaluate qualitatively the water balance components of a carbonic watered formation and an overlying one. The stage of mining closure helped obtain the forecasting hydrodynamic solutions. The efficiency of measures, concerning reduction of water ingress into mine workings and the mitigation of surface ecological effects of mine flooding was evaluated quantitatively. It was determined that implementation of the water control procedures makes it possible to perform a 10–38% decrease in water ingress. In this context, they may be applied both independently and simultaneously. In terms of mine closure and flooding, a period of complete underground water recovery takes three years; in the process, surface zones of potential waterlogging and swamping are developed within the floodplain of Samara River, located at the territory of Western Donbas (Ukraine). The scientific novelty is to define regularities of hydraulic conductivity transformation of the rock mass of a mine field starting from the mine working roof fall, up to its compaction during the mine operation period. To do that, nonstationary identification problems were solved, using numerical modeling. The abovementioned makes it possible to improve the reliability of hydrodynamic prognoses and develop technological schemes to control water at the state of the mine closure.


Introduction
Researchers consider a process of underground coal mining as one of the major technogenesis types that is comparable with geogenesis in terms of intensity and metamorphosis depth [1,2]. Large-scale transformations of rock mass factor lead to changes in the state and characteristics of both the subsurface [3][4][5] and the whole environment [6,7]. Underground water is the basic dynamic subsurface component helping to perform energy, mass, and heat transfer [8,9] as well as to transfer common groundwater contaminants [10][11][12][13]. Water can vary the properties of the geosystem and environment [14], and activate exogenic processes and even endogenic ones [15,16]. Forecasting of hydrodynamic changes in the underworked rock mass and prevention of their negative effects are relevant scientific and practical problems [17][18][19][20][21].
Geomechanical processes, arising in the rock mass in the process of mining with complete gob caving, vary the physicomechanical characteristics of the rock mass [22][23][24][25]. Technogenic fissuring progresses within the caving areas, resulting in the increased hydraulic conductivity as well as in the increased capacity properties of the rock mass that experiences permanent temporal variations during the whole period of the deposit operation [26][27][28]. In terms of a part of the Donetsk coal basin (Central Donbas, Ukraine) and steep pitch, the authors of papers [29,30] identified that under the influence of changes in geomechanical processes in rock mass, stipulated by its repeated underworking, the hydraulic conductivity varies according to hyperbolic dependence with stabilization at 10 −7 -10 −6 m 2 level.
In this context, the filtration coefficient of rocks experiences almost a 100-times increase; the gravitation capacity increases up to 14 times; and the infiltration supply increases up to 3.5 times. Within the areas of tensile strain concentration, rock disintegration increases 1.2-1.4 times. At the same time, the compression strain is 1.6 times more than the tensile one, resulting in rock solidification and hydraulic conductivity bed formation [31].
Studies, described by paper [32], showed how changes in the geomechanical situation of the rock mass influence its state. It was identified that the major rock mass changes take place at a distance from the mine working periphery of up to five times its radii. An increase in the distance between the fissures from 0.5 m up to 1.5 m decreases the filtration coefficient in the mine working neighborhood by 25-50%.
Paper [33] evaluated the formation of a secondary aquifer according to the results of observations about the movement and compaction of caved rocks. The compaction of the caved rocks and their cementation comprises three periods; mainly, the processes depend upon the load being developed after the rock moves to a day surface. The period of the rock movement to a day surface is 2 to 6 months after the stope displacement. In this context, the authors believed that new structures, having other filtration properties, are formed. However, they did not specify their quantitative characteristics.
The majority of the studies concern the analysis of the filtration parameters of the rock mass at the stage of active displacements and the formation of zones of water conducted fissures [34][35][36]. Moreover, they differ by orders of magnitude in terms of various data [37]. Hence, they are not common depending upon specific mining-geological and technological conditions of coal extraction [38][39][40]. The problem of the transformation of filtration and capacity characteristics of mined-out areas remains uncertain over 40-60 years of a coal deposit operation; nevertheless, their influence on the reliability of hydrogeological forecasts is determinative [41].
Under the conditions, complicated by numerous factors with many unknowns, digital mathematical modeling of geofiltration is perhaps the only justified method to analyze the rock mass of a mine field disturbed technogenically.
The idea is also supported by the authors of paper [42], who dealt with the development of the hydrodynamic model of the Central Donbas area where mine fields are connected hydraulically. Surface waterlogging areas, being almost 33% of the territory, were identified in terms of the results of prognosis modeling of successive and simultaneous mine floods. A differential approach was proposed to select the drainage measures protecting the mine field surfaces against waterlogging.
A numerical model of a technogenically disturbed rock mass within the Kuznetsk coal field was developed by the authors of paper [43] to evaluate the amount of underground water in the zone of drainage impact of underground mine workings. At the same time, the geofiltration model developed by the authors of paper [44] is the basis to evaluate the capacity amount of underground water in the flooded mine as well as to forecast changes in the hydrodynamic mode while using the water as a heat carrier with the help of a thermal power module.
In the context of mining regions, the technogenic mode of underground water is characterized by the breach of hydrodynamical, hydrochemical, and thermal behavior; it is notable for its specific dynamisms and dependence upon mining and geological factors [45][46][47][48].
It should be mentioned that each mining region in Ukraine and in the world has similar features of technogenic changes in the underground and surface hydrosphere [49][50][51][52][53]. They are manifested as the formation of subaerial subsidence troughs and waterlogging areas, quantitative and qualitative troubles in the functioning of the drinking water abstraction, and ground water contamination due to infiltration losses from containment ponds [54,55]. At the same time, solving the problem of safe mining under conditions of complex and multi-factor geofiltration processes is still relevant [56,57], especially during the closure of a mining enterprise [41]. Hence, the research objective is to analyze the dynamics of the underground water of a technogenic medium of a mine field to achieve the safety of extraction operations and near-surface territories at the stage of closure of a mining enterprise.
To succeed, the following problems should be solved: - The substantiation of a scheme of a geofiltration numerical model of a mine field, taking into consideration both natural and technogenic factors of underground water mode formation; - The identification of basic factors forming water ingress in a mine as well as the regularities of the transformation of filtration parameters of the disturbed rock mass; - The determination of the influence of the mine working flooding on the ecological state of the neighboring territories and substantiation of engineering measures of water control.
While solving nonstationary identification problems using a numerical model and determining regularities of the hydraulic conductivity transformation of the mine field rock mass starting from the moment of the mine working roof fall, up to its solidification along the mine operation period, we can be sure in reasoning the scientific novelty of the paper, as it will help improve the reliability of the hydrodynamic prognoses and develop scientifically substantiated technological schemes of water control at the mine closure stage.

Methodology
Scientific and technical approaches to model geofiltration within coal deposits were substantiated in terms of the hydrodynamically open type of the M.I. Stashkov mine (Donetsk coal basin "Donbas", Ukraine) where coal seams have an active hydraulic connection with the thick formation of waterflooded Meso-Cenozoic rocks. Figure 1 shows the geographic location of the site under study. The hydrodynamic situation of a part of the deposit is complicated by coal occurrence in inundable river portions where the surface watercourse may be an additional formation source of inflows into mine workings and the floodplain is a zone of potential waterlogging in the process of a mine closure and the recovery of the underground water level. The algorithm to model numerically a technogenic mode of underground water is based upon differential filtration Equation (1), solved by means of numerical iteration methods with the help of a system of finite-difference equations: where H is the unknown pressure function; W is the feeding (discharge) throughout the aquifer; T x , T y are the conductivity of the aquifer in terms of linear coordinates x and y, respectively; µ* is the coefficient of the elastic water yield; and t is time. The finite-difference approximation of Equation (1) is based upon the spatial and time-dependent discrete representation of a geofiltration flow in terms of which continuous filtration field H(x,y,z) is replaced by an imaginary network area characterized by H(x j ,y i ,z k ) pressure values in each node of the network per time moments with ∆t interval: where Φ is the filtration resistance of the flow between analytical blocks expressed through ∆x, ∆y, block geometry and conductivity values. Papers [60][61][62] specify general methods of the development of the numerical, mathematical geofiltration model involving several stages, among which are the following: -Systemization of data of mining and geological conditions; - The model identification in terms of the physical and dynamical analogy of the object; -Forecasting solutions (generally, integrity of the latter is identified by means of the schematization quality of the conditions as well as the model verification).  To solve differential equations of geofiltration and to model the filtration processes, the authors applied the MIFp software developed in Ukraine by a branch laboratory of hydrogeological studies of the Dnipro Department of the Ukrainian Geological Prospecting Institute. The software is protected by copyright; unfortunately, it cannot be sent through the internet. In turn, it is used to solve various hydrogeological problems starting from the evaluation of underground water reserves up to the assessment of the disturbed mode of underground water within mining regions. The program was tested while solving hydrogeological problems of mine closure in the Central Donbas area and hydroecological problems of the Kryvyi Rih iron-ore field (Ukraine). The results of early studies, published in [41,57,63], were obtained using corresponding software; their validity was confirmed practically.
The software is meant for numerical solutions of differential equations describing processes of plane and spatial filtration as well as mass transfer within nonhomogeneous stratified formation, involving the following: interaction of underground and surface water; migration through separating layers; infiltration feeding; evaporation depending upon the ground water depth; water surface entering; operation of water wells and pumping wells; time-dependent changes in boundary conditions and parameters; transition of a pressure filtration mode into a pressureless one and vice versa; layer thinning; water conductivity dependence of pressureless levels upon the underground water level; and rock anisotropy. The algorithm advantage is the possibility to represent the thinning (outcrop under the overburden rocks) of coal seams, being the determining factor while forming the water inflow into the mine workings.
The representation of the aquiferous complex of overburden deposits as well as natural and technogenic coal systems has become an important stage in the process of mine-field model development. The abovementioned has made it possible to use the model as a permanent one for each operating stage of a mining enterprise while responding promptly to rapid changes in a hydrogeological situation and solving engineering problems to optimize water control measures.

Analysis of Natural and Technogenic Conditions, and Substantiation of Hydrodynamic Scheme of a Mine Field Model
The analyzed mine field is within the Dnipropetrovsk Region (Ukraine) in the Samara River basin; its area is almost 83.6 km 2 .
The mine field territory is represented by slightly hilled plain separated by a river and ravine network. The Samara River crosses the mine field, dividing it into two uneven shares: large-a right-bank part being 30-80 m higher than the floodplain, and small-a left bank part represented by the floodplain where absolute marks are +62 . . . +72 m.
The structure section consists of two hydrogeological levels. The upper one involves Meso-Cenozoic overburden deposits represented by loose sandy and clay Paleogene rocks (Bucha, Kyiv, and Kharkiv measures); Neogene rocks (Sarmatian measures); and Quaternary rocks (Figure 3). The formation thickness varies from 40 m in the Samara River valley up to 140 m within a watershed.
Meso-Cenozoic deposits do not include continuous impermeable layers. Hence, all water-bearing strata are interconnected, forming a unified aquifer system. Under the natural conditions, its supply was performed at the expense of atmospheric precipitation infiltration; the water drained into the Samara River and its feeders.
The lower hydrogeological level is represented by a coal measure of the Mississippian period. The thickness of the layer is up to 410 m; it is the key productive formation in Western Donbas. Underground water is confined to coal, sandstone, and limestone seams, alternated with impermeable argillite and aleurite. Water of the deposits is fissure-stratal and pressured. The filtration coefficient of the sandstone and coal varies in a broad range; its average value is 0.74 m/day. The natural hydraulic conductivity of the rocks decreases in a structured section; depths deeper than 250-300 m unwatered coal and sandstone are prospected.
In the process of the geofiltration model development, the mine field is approximated with the help of grid with 77 × 43 cells ( Figure 4). As for the actual conditions, the geometry of the mine field is 15.4 × 8.6 km. Each grid point represents the seam hypsometry as well as the filtration and capacity parameters of both aquifers and separating formations.  The most watered sandstone with 30-50 m thickness occurs within the roof of the C 5 + C 5 B coal seam. The sandstone outcrops under Meso-Cenozoic deposits. In terms of the feature complex, the mineral is characterized as alluvial or channel. Since the mine is of a hydrodynamically open type, in addition to capacity reserves, the water inflow is formed owing to the involved resources at the expense of the flow of overburden deposit aquifers as well as the surface water inflow. The representation of coal seam mining within productive carbon formation to involve the actual hydraulic connection between the upper and lower hydrogeological levels was performed with the help of mine field contouring in terms of C 10 , C 8 , C 7 , C 6 , C 5 , C 4 , and C 1-3 coal layers ( Figure 5). The abovementioned has made it possible to define extra feeding within the outcropping areas of the coal seams as well as the paleochannel sandstone under Meso-Cenozoic deposits.
Hence, within a vertical section, the mine field model is represented by means of 8 layers, where upper layer 1 demonstrates the aquifer system of overburden Bucha-Kyiv deposits, and 2-8 layers show the mined-out coal seams and paleochannel sandstone C 10 , C 8 , C 7 , C 6 , C 5 , C 4 , and C 1-3 , outcropping directly under layer 1 (Figure 6a).  The external hydrodynamic boundary of analytical layer 1 of Bucha-Kyiv deposits (Figure 4) demonstrates the natural feeding system as well as the underground flow into the Samara River (Derichlet's boundary condition, H = const). In terms of flow lines, Newmann's boundary condition with zero loss (Q = 0) is specified.
In the plan, hydrodynamic boundaries of coal seams (analytical layers 2-8) coincide with the mine field boundary determined in terms of zones of tectonic disturbances (i.e., Pozdovzhnyi, Petrovskyi, Bohdanivskyi, and Nikolskyi faults). They are impermeable barriers for underground water filtration (Q = 0). The boundary of the coal seam outcrop under the watered formation (Figure 6b) is specified by Newton's boundary condition (H = f (Q)), representing the interconnection between the levels and losses resulting from mine drainage.
To represent the natural mode of ground water, the boundaries are specified along the flow lines (i.e., perpendicularly to hydroisohypses); at the stage of hydrodynamic changes, they are not waterproof barriers. The river has water-permeable underflow deposits; hence, its impact on the water inflow into mine workings is recorded when mining operations approach the coal seam outcrop under the floodplain. As for the role of tectonic displacements within the carbonous formation, practices of their crossing by mine workings mean that they are waterproof barriers for underground water.
Internal boundaries of the model are represented by mine workings, being also a drainage contour for sandstone occurring within the immediate roof or within the water fissure zone. In this context, sandstone is considered part of the mine working, i.e., as a drainage periphery with the specified level of underground water decrease at a mark of the coal seam floor. The model scheme has helped to involve the hydraulic connection between the carbonous formation and overburden deposits as well as a prerequisite concerning the availability of the paleochannel (i.e., alluvial) sandstone within the coal seam roof.

Identification of the Model
The most important stage to develop a mine field model involves solving identification or epignosis problems in terms of the natural hydrodynamic situation and in terms of the situation disturbed by mining. Then the solutions are applied to specify the boundary conditions, filtration parameters, and capacity parameters, using the underground water level as well as water rate inflows.
Finally, the identification results of the mine field model are evaluated in terms of the object functioning data, i.e., functional correspondence is determined when the model reproduces the processes for which the actual behavior was recorded under the natural conditions. The development process of the mine field model reproduces the chronology of coal seam mining during the periods in which were recorded the most typical changes in water inflows (Figure 7). The verification process of the mine field model identifies the reliability of the forecasted solutions as well as the water control efficiency at the stage of the deposit operation and at the stage of mine closure.

Results and Discussion
A process of natural hydrodynamic pattern reproduction, using the developed mine field model for the mine preoperational period helped to support the reliability of the actual underground water data as well as to balance the numerical model (Table 1). It was determined that under natural conditions, the underground water of the Bucha-Kyiv deposit flows into the Samara River. The fact is supported by the groundwater contour near the river as well as by the moving direction of underground flows (Figure 4). The Samara intake operation with 2000-3200 m 3 /day of water discharge is a development pressure source for the period.
The research results, connected with the model response to fluid dynamic excitation by means of the introduction of Derichlet's boundary condition, according to the mining schedule, showed that in terms of each period, the underground water level of overburden deposits corresponds to actual ones (Q act ) as well as the recorded water inflows (Q rec ) ( Figure 8).
The average error of the solution is not more than 10% in terms of the water inflow rate; it is 7% in terms of the level. In this context, the mismatch of the total model balance is not more than 0.01%. Hence, in accordance with mining schedules, the authors characterize periods of coal seam mining where the most typical changes in water inflows were recorded.
The results of the monitoring observations in terms of the 3756a well, represented in Figure 9, are the control data for the changes in the underground water level in the Bucha aquifer in the process of coal seam mining.  Changes in the direction of filtration flows show that from a zone of underground water drainage, the river transforms into a feeding zone of the Bucha-Kyiv aquifer ( Figure 10). A decrease in the underground water level in the Bucha deposits achieves 13.5 m, relative to the preoperational state. Period 2 (years of 1990-1994) is characterized by maximum inflows into the mine as well as the mining of C 8 , C 6 , C 5 , and C 4 coal seams. The key task of the research stage is to evaluate quantitatively the factors and formation reasons of abnormally high water inflows in terms of the C 5 seam.
During the period in which the total mine inflow is 1788 m 3 /h, mining takes place within the zone of paleochannel sandstone (C 8 and C 5 seams) with water conductivity, being comparable with that of the Bucha deposits and in the immediate neighborhood of the Samara River bed (at a 200 m distance). Analysis of balance components (Table 2) demonstrates that 69% (i.e., 570 m 3 /h) of water inflows within the C 5 seam is leakage from the Bucha aquifer and river water underrun. Hence, during the period, 40% of the water inflow into the mine is water filtration from the Samara River. Table 2. Components of water balance in terms of C 5 seam.

Expandable Components of the Balance Water Inflow Value, m 3 /h
Mine working-underground water interconnection (inclusive of leakage from Bucha aquifer at the expense of filtration from the Samara River)

(570)
Flow through the level roof 0.9 Flow through the floor 53.7 Consequently, the conical depression displacement in analytical layer one toward the river and underground water, lowering down to 31 m, correspond to the data of monitoring observations in the year of 1994 ( Figure 11). It should be noted that during the time frame, the water intake in the Samara inlet experienced more than a twofold decrease (i.e., from 3200 down to 1800 m 3 /day).
Period 3 (years of 1994-1998) is characterized by water inflow lowering down to mine workings (i.e., from 1788 to 1454 m 3 /h). C 8 , C 7 , C 6 , C 5 , and C 4 coal seams are mined. Water inflow within the main haulage roads of C 10 seam is recorded. For the period, the decrease in the underground water level of the Bucha-Kyiv deposits is 29 m. In 1995, hydraulic discharge of the Samara water inlet was 180 m 3 /day. It was shut down in 1998.
During Periods 4-6, a decrease in the total mine water inflow is lasting (Figure 7) since the mining zone distances from the area where coal seams outcrop under the overburden deposits.
The modeling of the actual dynamics of water inflows and their levels at each stage of the mine field extraction has helped to determine a number of regularities concerning the formation of the underground technogenic mode: 1.
If natural water conductivity of coal seams is 0.5-1.5 m 2 /day, then actual water inflow values may be achieved when transit zones of hyperpermeability are available. Within a mine field, they are represented by means of sandstone formations with up to 50-60 m thickness and up to 30 m 2 /day water conductivity. The sandstone occurs in the coal seam roof and outcrops under the watered overburden deposits. According to lithological and facial analysis, it is alluvial sandstone of the Early Carboniferous paleochannels; regularities of its spreading over the area correspond to the heightened water ingress during mining.

2.
A total of 70% of the water inflow into a mine are provided from Meso-Cenozoic deposits; 30% are the capacity reserves of carbonous formation. As a result, the formation of a conical depression up to 30 m takes place within the Bucha-Kyiv aquifer as well as changes in directions of filtration flows in the neighborhood of the Samara River. Under the disturbed conditions, the river stops being a drainage zone transforming into a feeding zone for overburden deposits. Thus, it becomes an extra source for water mine workings.

3.
Movement of the conical depression within the Bucha level performs synchronously by putting into operation a new coal seam toward a zone of outcrop under overburden deposits.

4.
Maximum involvement of the Samara River flow in the mine drainage was recorded while extracting the C 5 seam within the river plain flood (i.e., at a 200 m distance from its bed) and while approaching safe mining boundaries. Hence, a maximum of 69% (i.e., 570 m 3 /h) of the water inflow within the seam depends upon the river water leakage. The index drops down to 17% if the mining distances from the coal seam outcrop under the Bucha-Kyiv deposits and the thickness of the underflow decreases within the eastern part of the mine filed.

5.
A value of the water inflow into a mine cannot correlate with the increased mining area. Within the area mined out with the roof caving, carbonous formation thickness is a time-variable value which increases 10-15 times in the process of rock displacement. The value decreases after 5-10 years and stops performing a drainage function after 15-20 years ( Figure 12). It is possible to determine such transformations for the specific mining and geological conditions while solving long-term (i.e., during 15-20 years) time series of non-stationary identification problems, being a mandatory stage for the methods developing adequate models of mine fields, and obtaining durable, high-integrity forecasts. Relying upon the research results, concerning the temporal transformation of filtration properties of the disturbed rock mass, the authors propose three possible variants of activities to be performed during the phaseout stage to support the drainage mode in the M.I. Stashkov mine as well as a variant to protect a plain flood against waterlogging, while constructing an irrigation intake in a productive formation of paleochannel sandstone with increased hydraulic conductivity. The dynamics of rock mass flooding differ from the dynamics of level recovery in a mine shaft. Within the stopes of C 10 , C 8 , and C 6 seams, flooding lasts during the first year since they are already mined out and levels are recovered partially. The waterlogging of mine shafts and the recovery of the level within C 5 and C 4 seams lasts for three years ( Figure 13). The Bucha-Kyiv aquifer response to drainage-off in the mine is manifested as a level recovery over three years with the formation of waterlogging zones along the Samara River bed. The influence of the underground water level recovery on the operation of neighboring mines is seen in the increased water inflow, being 50 m 3 /h. Similar to the preoperational period, the plain flood of the Samara River, with almost a 9.75 km 2 area, remains a zone of potential waterlogging and swamping where the absolute height is less than +75 m ( Figure 14). It should be mentioned that the areas were not underworked. They are beyond the zone of geomechanical influence of mining.

2.
Variant 2. Mine flooding with preservation of the available underground system of the main mine drainage. Figure 14. Fragment of the mine field section: 1, 2, 3-projections of C 5 seam on a day surface (plane bed, seam periphery, and influence area, respectively); 4-boundary of potential waterlogging; 5-zone of water wells.
Preservation of the main mine drainage at a 225 m level (elevation is −115 m) helps lower the underground water level from the overburden deposits by 4 m within the zone of potential (i.e., preoperational) waterlogging. However, in the longer term, the variant involves preservation of the available system for drainage of the mineralized mine water with its further export to the Samara River in the amounts of 5869.2 thousand cubic meters, annually. The fact cannot be considered to be an environmentally friendly measure.

3.
Variant 3. The mine abandonment with drainage mode support and by-pit equipping with submersible pumps.
To control the environmental safety of a near-surface aquifer system (i.e., to prevent waterlogging and changes in the underground water as well as in the surface one), it is required to keep underground water marks lower than the foot of overburden deposits within the mining area. To keep the underground mine water level (i.e., underground water of carbonous deposits) within the undermining zone with +20 m elevation being the minimum level of the overburden deposit floor, drainage using submersible pumps is considered in terms of following subvariants.
Mine water elevation levels in a mine shaft are the following: The response of the overburden deposit aquifer to the operation of the submersible pumps shows that the drainage influence radius is quite smaller, compared with variant 2; within the potential waterlogging area, the level lowering is 0.1 to 1.2 m. In this context, the decrease in the potential waterlogging area is only 15%.
The use of submersible pumps with a drainage rate relative to 225 m level (elevation is −115 m) or 300 m level (elevation is −190 m) should involve free water migration from the overlying levels. If the drainage variants are applied, then the situation in the Samara River plain flood is similar to that formed at the final stage of a mine field extraction. However, the abovementioned involves long-term drainage and export of the mineralized mine water to the river.

4.
Variant 4. Alternative to protect the plain flood against waterlogging by means of water storage facility construction within the productive strata of paleochannel sandstone being of heightened hydraulic conductivity.
Substantiation of the variant as well as the water well location follows the solution of the identification problems and practices of C 5 seam mining in the immediate neighborhood of its outcrop under overburden deposits and the Samara River bed. During 1994-1995, abnormal water inflows into the mine were recorded (i.e., 1817 m 3 /h; 822 m 3 /h of them within C 5 seam); the level of the Bucha deposits achieved 30 m. In this context, the paleochannel sandstone, occurring in the main roof of the C 5 seam with 14-22 m thickness, served as an extra zone for the transition of the underground water of overburden deposits.
Defining the reasons and factors of abnormal water inflow manifestation makes it possible to use the disturbed zone of the underflow sandstone within a −55 to −50 m elevation range to reduce natural waterlogging areas after mine abandonment. Figure 14 demonstrates the well location in the plan.
Sufficient lowering of the level within a potential waterlogging area (i.e., up to 2 m) is recorded if wells with a total of 150 m 3 /h or 3600 m 3 /day operate. According to the dynamics of mine field flooding, its operation should be foreseen after 6 months from the onset of the mine field flooding. The average mineralization of the water to be pumped out will be almost 2 g/dm 3 .
The variant is the most acceptable, ecologically, with synchronous level mode monitoring, which will help use the industrial water for internal needs of the mine office as well as increase the reclaimed land area up to 90%.

Conclusions
As a result of the carried out complex of scientific studies, the paper helped to solve topical problems related to the scientific substantiation of a scheme of a geofiltration model of a mine field, taking into consideration both natural and technogenic factors of underground water mode formation. This was possible owing to solving nonstationary identification problems and determining regularities of the rock mass hydraulic conductivity transformation (involving a numerical model), being of scientific and applied value.
It is possible to forecast underground water dynamics within the technogenic environment of a mine field if a permanent numerical model is developed. The model is identified using reconstruction of the whole hydrodynamic history of deposit mining in terms of characteristic time periods.
The determined regularities of changes in the technogenic hydraulic conductivity of a mine field of rock mass increase the validity of long-term hydrodynamic forecasts, making it possible to develop technological procedures for water control. The hydrodynamic model of the M.I. Stashkov mine field is an integral part of the geofiltration model of a share of the Donetsk coal basin (Western Donbas) for solving the current problems of water control as well as the environmentally friendly process of mining phaseout.
According to the forecast analysis, in the context of the recovery of the underground water level of carbonous and overburden deposits, areas of potential waterlogging and swamping are formed in the Samara River plain food, which existed during the preoperational period. The areas were not underworked; moreover, they are beyond the zone of geomechanical influence of the mining operations.
Alternatives to the support drainage mode in the M.I. Stashkov mine were considered as well as an option to protect the floodplain against waterlogging by means of water storage facility construction within the productive strata of the paleochannel sandstone, being of heightened hydraulic conductivity. Three variations to support mine drainage were proposed as engineering measures to protect the river plain flood against waterlogging. The use of submersible pumps for water pumping out from the mine shaft over 140 m level is an inefficient procedure from the viewpoint of a decrease in waterlogging areas since they are within the mine field periphery. Drainage support at the 225 m level makes it possible to approach the hydrogeological situation, being typical for the deposit mining period; however, it involves long-term drainage of the mineralized mine water and its discharge to the Samara River.
The substantiation of water storage facility construction has become possible, owing to solving epignosis problems for a period of abnormally high inflows into the mine. The abovementioned helps use the productive formation of paleochannel sandstone within the C 5 seam roof as a transit zone of heightened hydraulic conductivity to protect the plain flood against waterlogging. The results of the forecasting hydrodynamic solutions concerning the mine flooding outcomes are the foundations to develop its abandonment project.
A hydrodynamic model of the M.I. Stashkov mine field will become a constituent part to develop a regional hydrodynamic Western Donbas model. It is planned to use the regional model for studying and considering the regional water balance control alternatives in terms of the mining closure to minimize environmental risks and substantiate the geocirculation of the heat supply systems with the use of capacity resources of the flooded mine fields.