Next Article in Journal
Optimization of Liquid-State Anaerobic Digestion by Defining the Optimal Composition of a Complex Mixture of Substrates Using a Simplex Centroid Design
Next Article in Special Issue
Numerical Groundwater Model to Assess the Fate of Nitrates in the Coastal Aquifer of Arborea (Sardinia, Italy)
Previous Article in Journal
LithoSFR Model for Mapping Groundwater Potential Zones Using Remote Sensing and GIS
Previous Article in Special Issue
Application of Variable Weight Theory in the Suitability Evaluation of Regional Shallow Geothermal Energy Development
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Dynamic Programming Models for Improvement of Technological Approaches to Combat Negative Water Leakage in the Underground Space

by
Sérgio Lousada
1,2,3,4,5,*,
Svitlana Delehan
6 and
Andrii Khorolskyi
7
1
Department of Civil Engineering and Geology (DECG), Faculty of Exact Sciences and Engineering (FCEE), University of Madeira (UMa), 9000-082 Funchal, Portugal
2
CITUR-Madeira-Research Centre for Tourism Development and Innovation, 9000-082 Funchal, Portugal
3
VALORIZA-Research Centre for Endogenous Resource Valorization, Polytechnic Institute of Portalegre (IPP), 7300-110 Portalegre, Portugal
4
Research Group on Environment and Spatial Planning (MAOT), University of Extremadura, 06071 Badajoz, Spain
5
RISCO—Civil Engineering Department, University of Aveiro, 3810-193 Aveiro, Portugal
6
Centre for Interdisciplinary Research, Uzhhorod National University, 88000 Uzhhorod, Ukraine
7
Laboratory of Mining Problems, Branch for Physics of Mining Processes, National Academy of Sciences of Ukraine, 49005 Dnipro, Ukraine
*
Author to whom correspondence should be addressed.
Water 2024, 16(14), 1952; https://doi.org/10.3390/w16141952
Submission received: 15 May 2024 / Revised: 16 June 2024 / Accepted: 3 July 2024 / Published: 10 July 2024
(This article belongs to the Special Issue Water-Related Geoenvironmental Issues, 2nd Edition)

Abstract

:
The article solves an urgent problem, which is to develop a new approach to finding solutions to improve technological methods to combat negative water leakage in underround spaces. We propose the use of dynamic programming methods to select the optimal technology to secure such spaces. In accordance with the algorithm proposed in this paper, the problem was broken into a number of stages. At each stage, an optimal solution was sought (organisation of transport, delivery of materials to the destination, selection of materials, etc.). Thus, we applied a decomposition approach that allowed us to take into account the variety of parameters that affect the efficiency of the process. All these stages and their corresponding technological solutions were formalised by building network models. In these network models, vertices corresponded to solutions, and the distances between vertices (edges) corresponded to the value of the optimisation parameter. Thus, the shortest route from the initial to the final vertex corresponded to the optimal technological solution to combat negative water leakage in underground spaces. Based on the systematisation of data on technologies to combat water inflow into underground spaces, basic and refined models were developed. These models allowed us to take into account the risks associated with water breakthroughs into underground spaces. To minimise the risks, additional measures to combat water inflows are envisaged. In the practical part of this study, the results of the selection of a method with which to control water inflows are presented. This method involves the use of anchoring to reduce water filtration. According to the results of field observations, no water breakthroughs into the underground space were recorded.

1. Introduction

Throughout its history, mankind has undertaken the construction of structures in underground spaces. At first, it was caves, then tunnels, vaults, labyrinths. Today, underground structures are used not only in production activities, but also in public construction (for example, subways). Regardless of the stage of development of science and technology, there was always a problem in dealing with negative water manifestations. Negative water effects include water filtration, the soaking of rocks surrounding the contour of an underground structure, and a high probability of water inflow into underground spaces. All these phenomena can be called unfavourable hydrogeological conditions. The negative phenomena associated with water infiltration seriously complicate the aforementioned processes, significantly increasing the cost and slowing the pace of the construction of underground structures, and sometimes making it impossible altogether. To solve these problems, mankind has created technological methods and means that make it impossible for water to enter the underground space. Such means include lining underground spaces with different types of materials, applying surfactants to minimise water filtration, etc. However, not every solution is optimal. Perhaps the most reliable technology reduces the risk of water leaks, but it is significantly more expensive or requires more human resources. That is why we set ourselves the task of developing decision-making approaches that would allow us to obtain the most rational solution. By “rationality”, we mean obtaining a solution regarding the technology intended to prevent water seepage into underground spaces that will meet two conditions: first, it will be optimal from the point of view of minimising the optimisation parameter (for example, the cost of constructing a support); second, it will cause the least harm to the environment, and, if possible, will allow for the disposal of production waste. One way or another, in the first stage, it all comes down to developing decision-making tools.
In this article, we propose a new methodology and decision-making tools for the design of measures to minimise the negative impact of groundwater in underground spaces.
However, first of all, it is necessary to analyse existing approaches to decision-making. There are several classifications of decision-making methods and tools in production. These classifications can be divided into three types: classification based on economic and mathematical models (this classification was popular until the 1980s), classification based on the results of decisions, i.e., quantitative or qualitative characteristics (this classification is still popular, but it became most popular in the 1990s), and classification based on methods of obtaining decisions and the duration of forecasting (2010—present).
The authors of [1] described economic and mathematical models, as well as the scope of the models. It is worth noting that the difference in the use of mathematical modelling methods lies solely in the limitations and types of parameter changes, as well as the requirements for the resulting solution. As already mentioned, in [2,3], a retrospective analysis of decision-making methods in production was carried out. According to these works, classifications can be divided according to the type of economic and mathematical model used [4], with the main differences being the nature of the parameter change and the design time. For example, if the process includes several stages that change over time, dynamic models are used. Another classification approach is based on the type of decision obtained [5], i.e., will we get a quantitative or qualitative answer as a result of the decision? For example, in the organisation of transport deliveries, a quantitative solution is appropriate, because we will get an answer to the question “Which of the alternatives in the transport chain has the lowest delivery cost, route length, etc”. When applying hierarchical analysis methods, a qualitative characteristic is obtained, that is, we will get an answer to the question “Which of the obtained solutions is the most satisfactory for our needs?” For example, when solving the problem of concentrator location, the connection with the enterprise, location, community interests, etc. are taken into account. The last classification type [6] takes into account the duration of the forecast (long-term or short-term), as well as the process of obtaining a solution (software, modelling on equivalent materials, simulation modelling, finite element method, etc.).
The following economic and mathematical methods are currently known: linear programming [7], quadratic programming [8], integer programming [9] (for example, when choosing mechanisation means by combinatorial optimisation, the number of equipment types should be an integer), stochastic programming [10] (in which nonlinear changes in parameters are possible) dynamic programming [11] (in which the process changes in time and includes several stages), geometric programming [12], inventory management [13] (in which the problem of allocation is solved), game theory [14] (when it is necessary to solve the problem of ensuring interests), network models and optimisation algorithms on network models and graphs [15] (which search for the shortest route corresponding to the optimal solution), models based on establishing dependencies and predicting changes in parameters [16], service theory [17], study of reliability as the ability of an object to perform its functions [18] (for example, the productivity of a technological chain is the ability of an object to reproduce a given level of productivity), forecasting and building correlation models [19], and finally, the application of simulation modelling methods using electrical circuits, information systems, etc. [20].
It is worth noting that the vast majority of the described models and methods [7,8,9,10,11,12,13,14,15,16,17,18,19,20] are used in short-term planning. In addition, these models contain quantitative rather than qualitative dimensions of the decision. However, the decision-making process is influenced not only by economic, technological, and social factors, but also by qualitative ones. If the objective function is, for example, the unit cost or transportation costs [21], the decision will be optimal in terms of a certain parameter, but there is a high probability that it will not meet the ‘quality’ characteristic, i.e., a set of features that characterises the object (e.g., for a mineral, the degree of depletion, ash content, etc.; and for an enterprise, infrastructure, location, communication with suppliers and consumers, etc.) [22]. It is also necessary to take into account production risks associated with unfavourable production conditions [23,24]. Water-related risks are also considered to be risks, but in addition to complicating production, the quality of water resources deteriorates [25]. All of this affects the quality of life of the population in regions where underground spaces are developed [26,27]. As noted, the methods [7,8,9,10,11,12,13,14,15,16,17,18,19,20] provide an answer about quantitative characteristics, not qualitative ones. If the transportation problem is solved [28,29], the solution will be optimal and quantitative, for example, in terms of specific transportation costs, the length of the transportation chain, etc. However, this does not take into account the wishes of suppliers, consumers, etc. That is why a classification system has been developed that takes into account the types of solutions and models used [30,31]. These models are divided into information and optimisation models. Optimisation models provide a quantitative answer, while information models reflect the qualitative characteristics of expert groups. Informational models include mathematical programming methods [32,33,34,35,36,37,38], as well as game theory and decision-making criteria under conditions of uncertainty (Wald, Laplace, Hurwitz, Savage, Maximin). Information models are based on the AHP hierarchy analysis method [39] developed by T. Saaty; they comprise variations of the AHP method [40,41,42,43,44,45,46,47,48,49,50,51], a method that involves the use of fuzzy sets [52,53]. In these methods, hierarchies are built and priorities are set. As a result, we get a high-quality solution that takes into account the wishes of the expert group. The disadvantages of these methods include situations in which it is impossible to build a connectivity matrix, i.e., when one of the parameters significantly outweighs the others (for example, the cost of purchasing equipment) [54]. Another disadvantage is the need to use expensive software [55,56]. In [57], a separate group of approaches based on simulation modelling methods is identified [58]. Simulation modelling is used to model groups of production parameters by breaking larger problems into local ones and applying appropriate algorithms to find the optimal solution for each group [59,60,61].
Based on the analysis of recent works [62,63,64,65,66,67,68], the following generalisation can be made, which forms trends in the development of approaches to process design: the optimality criterion is “quality”, but it has different quantitative indicators; the method of breaking down larger problems into local ones should be applied; but only after finding the optimal solution in the previous stage should one proceed to finding a solution in the next stage, in other words, applying the principle of ‘Belman optimality’. The following questions need to be answered for the task of preventing water leakage in underground spaces:
What technology should be chosen? The chosen technology should not only prevent water from entering the underground space but should also be inexpensive. The cost can be expressed not only in monetary terms, but also in terms of the amount of resources required. Of course, the use of concrete-based materials is cheaper than technologies that use composite materials, but the former are more resource-intensive and cause more environmental damage.
How to minimise the negative impact on the environment? Or, which technology is environmentally preferable? Here, we analyse not only the amount of resources required to prevent water spills, but also the environmental damage caused by the materials and technologies used. Nowadays, it is no longer enough to analyse and develop economic strategies; environmental strategies should be developed as well [69,70].
Is it possible to use production waste to prevent water leakage? This possibility can be applied not only in the production process, when the technology involves leaving materials in the underground space, but also on the basis of the comprehensive utilisation of production waste.
The answers to these questions are contained in the present study. We have proposed not only approaches to decision-making, but also paid considerable attention to the possibility of minimising the technogenic impact on the environment.

2. Materials and Methods

2.1. Study Area

The approaches presented in this paper and the proposed dynamic programming model are valid for countries where minerals are extracted underground at considerable depths, i.e., when there is a need to secure the mined space. An additional condition is the likelihood of water inflows into such underground spaces, i.e., the presence of waterlogged layers in the rock mass. The above technological operations are described for conditions with unstable or medium stability of roof rocks. The rocks of the workings’ base are prone to uplift. In the absence of unfavourable hydrogeological conditions, as well as when there is no need to secure the excavated space (the rocks are stable), the above models are not applicable. For the conditions of Eastern Europe and countries where mining is carried out under unfavourable mining and geological conditions, this approach is appropriate. In addition, it is worth noting that the main area of application is the extraction of seam minerals (coal, potash, manganese, etc.). For iron ore deposits, this approach can be applied only in capital workings when there is a need to secure the mined space.
The proposed methodology can be applied to other processes that take place in underground spaces. These are tasks that are not only related to securing the mined space to prevent negative manifestations of water; dynamic programming models can be used to develop transport chains and emergency response plans. It is worth noting that we have presented the application of the dynamic programming method in terms of minimising the cost of construction of workings (the forward problem), but it is also possible to apply and solve the reverse problem. If a mineral is being extracted and there is a need to extract the maximum amount of it, i.e., to maximise the cost of 1 tonne of mineral, then the inverse problem is solved.
When used in regions with different hydrogeological characteristics, existing technologies for the control of water inflows should be taken into account. For example, in Ukraine, Polesia, and the Czech Republic, anchoring is a very common method in underground spaces; it ensures the monolithicity of the waterlogged rock mass. At the same time, for the conditions in the Republic of Kazakhstan, where the water content is lower, the use of frame supports is more traditional. Also, if such projects are not related to the extraction of minerals, for example, in the construction of subways or underground storage facilities, concrete supports are traditionally used. The only differences are in the existing technologies to control water inflows into underground spaces, in accordance with the country’s conditions.
Dynamic programming has several advantages:
Reduced execution time, i.e., by breaking down a task into smaller subtasks to avoid repetitive calculations, thereby reducing the overall execution time.
The efficient use of resources, i.e., the memorisation of results avoids repeated calculations and reduces resource requirements.
  • The optimal solution: dynamic programming usually provides an optimal solution to a problem.
Limitations of dynamic programming:
  • Dynamic programming has certain limitations.
  • Dependence on the structure of the problem: to use dynamic programming, the problem must have an optimal structure, which may not always be achieved.
  • The need for a large amount of memory: the memorisation of results can require a significant amount of memory, especially for large problems.
The time complexity of the Bellman-Ford algorithm depends on the number of vertices in the graph. To reduce the number of calculations, you can use the time-to-memory ratio, and with more memory, the execution time decreases.
Part of the data calculated by the algorithm about the distances to the vertices of the graph can be stored in memory. Then, if you need to find a path from the same node, the path can be loaded from memory instead of having to be calculated again.
The data calculated by the Bellman-Ford algorithm can be partially stored, which will allow you to choose the optimal time–memory ratio for a particular system. Unlike the Dijkstra algorithm, the Bellman-Ford algorithm allows negative edge weights to occur in the graph. In addition, the complexity of the Bellman-Ford algorithm is O(VE) (V is the number of vertices, E is the number of edges), while the complexity of the Dijkstra algorithm in the classical implementation is O(V2E). Thus, the Bellman-Ford algorithm was chosen to implement data transmission in a multi-parameter system.

2.2. Dynamic Programming Model

To solve this problem, we used dynamic programming methods.
Dynamic programming models can be applied in less controlled environments, for example, to design mining, processing, and beneficiation processes. In this case, not only quantitative indicators should be addressed (the reduction of production costs and expenses, the maximisation of profits), but also qualitative indicators should be taken into account. The main purpose of dynamic programming models is to make decisions in multi-stage processes, i.e., to find the optimal/extreme (largest or smallest) value of the parameter to be optimised. In this case, “multistage” means a certain process, technology, process structure, which contains a certain sequence and the possibility of its division into stages. The term “dynamic” indicates that the design process is time-consuming. However, it is worth noting that dynamic programming models, due to software complexities in terms of implementation, should only be used for long, multi-stage processes, such as the consolidation of excavated spaces under complex hydrogeological conditions. This is because at the decision-making stage, a “compromise” decision is made, i.e., either to reduce the cost of mining or to not apply the measures, in which case the risks of water inflow and subsequent water filtration will increase significantly. The amount of environmental damage is much greater than the economic effect. Thus, dynamic programming models are used in several cases:
-
When the process is multi-stage and all stages are interconnected;
-
when the process is time-consuming;
-
when the efficiency of the process is influenced not only by quantitative but also by qualitative indicators.
In particular, the decision-making process can be considered using the following functional equation [71]
x ( t ) t = f ( x ( t ) , u ( t ) , t ) ; x ( τ ) = y ; τ t t 1 ,
  • τ—an arbitrary parameter;
  • y—is an arbitrary vector.
The direct immersion invariant method allows us to consider the production process in time t with a limited amount of resources. All possible solutions are considered, and each subsequent solution (alternative) must be no worse than the previous one (Figure 1).
To make a decision, we will present the process of finding the final solution as a set of procedures for finding optimal solutions at each stage (Figure 2). For this purpose, we will present Equation (1) in the form
x ( t + 1 ) = f ( x ( t ) , u ( t ) , t ) ; x ( τ ) = y ; τ t t 1 ,
Figure 2. Schematic of the computational procedure.
Figure 2. Schematic of the computational procedure.
Water 16 01952 g002
The process of finding a solution consists of examining possible solution vectors x for compliance with optimal solution domain F0. If the value of vector x at time t is FF0, you need to memorise the result and postpone it (memorise) on phase grid x; otherwise, the search should be interrupted. After that, based on responses x, x1, x2, we calculated the value of the optimisation parameter.
By using graphs and their corresponding network models, information can be presented in a compact and informative way. Dynamic programming methods can be used to select measures to combat water leakage in underground spaces. The idea is as follows: these stages and corresponding decisions are represented as vertices in a graph, and the distances between the vertices are the values of the optimisation parameter (e.g., cost). The shortest route will correspond to the optimal solution. It should be noted that the model is divided into stages and contains only real connections. That is, if there is a solution, we calculate the value of the optimisation parameter; if not, we do not eat the vertices or take the distance to be infinity.
The use of dynamic programming has a number of advantages:
Firstly, when making a decision, a multi-parameter problem is broken down into smaller ones, i.e., a decomposition approach is applied. The essence of this approach is as follows: at each stage, one “priority” control parameter is selected. It is the optimisation of this parameter that allows us to identify the optimal solution at this stage, which shapes the overall efficiency at the subsequent stages. For example, when choosing a technology for preventing water inflows into underground spaces, a number of tasks should be solved: to ensure the stability of the rocks surrounding the contour of the underground structure (to do this, we reduce filtration by using anchoring [72,73,74]); to ensure insulation against water ingress; to organise the delivery of materials, etc. However, in parallel, we solve the problems of material delivery, work organisation, and waste disposal—all these issues are important, but within a certain time frame. In other words, we take the general technology of combating negative water leakage into the underground space and break it down into stages: the delivery of materials to the warehouse, the delivery of materials underground, the delivery of materials to the place of construction of the support, work to minimise negative water leakage, etc. In each stage, we analyse alternative solutions. Without finding the optimal solution at a certain stage, we do not proceed to the next stage.
Dynamic programming models are based on the principle of Belman optimality [75]. According to Belman’s optimality principle, “If a control is optimal, then whatever the initial state of the system and the control of the system at an initial time, the subsequent control is optimal with respect to the state that the system will take as a result of the initial control”. Thus, the resulting solution will be optimal in terms of a certain optimisation parameter.
The search for a solution can be carried out in both forward and backward order. For the task of choosing the optimal technological solution to combat water inflows, the direct task is to minimise the cost of preventing water seepage into the underground space. The inverse problem is to minimise the waste that will be generated. In other words, the inverse problem involves an input condition that determines the quality of the solution.
After describing the principles of dynamic programming, we should focus on the decision-making procedure for improving technological solutions to combat negative water leakage in underground spaces.

2.3. The Development of a Model for Finding an Optimal Technological Solution in Complex Hydrogeological Conditions

Groundwater has a negative impact on the condition of underground structures and technological processes. In accordance with our objectives, we have established the relationship between hydrogeological factors and technological processes (Table 1).
Table 1 shows that the problem of choosing a technological solution comprises justifying not only the safest and most reliable method, but also the most resource-saving one. In addition to cost, the time factor should be taken into account.
Regardless of the chosen optimisation parameter (cost of construction of 1 m of workings, costs of ensuring stability of 1 m3 of underground space, etc.), this indicator depends not only on the cost of resources but also on the consistency of technological processes at each stage of construction and operation (Table 2).
Regardless of the type of technological solutions to combat water leakage in underground spaces, the cost components shown in Table 2 are identical, but different at intermediate stages, i.e., when comparing spray concrete and anchor supports, the amount of materials for the construction of support, time costs, labour intensity are analysed, but the level of mechanisation and technological stages are different.
As noted earlier, regardless of the type of support that prevents water inflows into underground structures, the stages of construction and operation of structures are identical, so the life cycle can be represented as a network model that takes into account alternative options (Figure 3).
Each vertex, from 1 to 17, corresponds to a separate solution, except for vertex 1, the start, and vertex 17, the end. These vertices are needed to connect the intermediate vertices (2–16) that correspond to the process solutions. In this case, the connections between the vertices are mutually exclusive, i.e., if the fastening is sprayed concrete (2) or anchored (3), there is no possibility of reuse (15), and there are no dismantling costs. The model can also be considered in reverse order: if it is not possible to dismantle the support (16) from the excavated space, then anchor supports (3) are not used.
As mentioned earlier, the model is structured by levels. To make an optimal decision, the shortest route from the first (1) to the last vertex (17) should be found. For this purpose, dynamic programming algorithms can be applied to network models and graphs; in our case, the Bellman-Ford algorithm [76]. For the algorithm to work, condition (3) must be met—the process must be uninterrupted in time. That is, to find the shortest route from node 1 to node j, some edge (arc) must be the final one.
f i = min i ( i , j ) ( f i + t i j ) ,
Additional conditions are reflected in [77]. According to these requirements, only realistic alternatives should be reflected in the model, and the parameters should be managed centrally; without finding an optimal solution in one stage, it is impossible to proceed to the search for an optimal solution in the next.
In this case, depending on the stage of construction of the support, the parameters to be taken into account will be different.
The developed approach allows us to take into account hydrogeological parameters. This is done in the following way: knowing the degree of watering of the massif and its hydrogeological characteristics, it is possible to provide measures to improve the stability of the workings. For example, knowing that there is a risk of water breakthrough, it is necessary to use anchoring, which will reduce the water permeability of the workings and cover the cost of making grooves for drainage. Additional technological stages can be combined with existing ones. To do this, the network model (Figure 3) can be adjusted for the degree of hydrogeological risk, and a certain cost increase factor (shaded vertices) can be added to the intermediate stage (wavy line).
Conversely, when considering the model in the opposite case, having an idea of the support used (e.g., anchoring), it is possible to weaken the connections in the model (Figure 4), i.e., the probability of water inflows into the workings will be reduced, which will reduce the cost of auxiliary operations.
Thus, at the beginning of the decision-making process, a network model will be proposed with which to choose the optimal technological solution. In the absence of information on unfavourable hydrogeological factors, the search for a technological solution is carried out according to the parameter “costs of water inflow control per 1 m3 of underground space”, etc. If it is known that the massif is waterlogged, then certain links in the network model should be strengthened, which will result in an increase in costs, while others should be weakened.
In the following section, we present the results of studies on the use of dynamic programming methods to select optimal technological solutions to combat water leakage in underground spaces.

3. Results

3.1. Systematisation of the Stages of Construction of Support to Combat Negative Water Leakage: Basic Model

We considered the conditions in Western Donbas (Ukraine) as the object of this study. The Western Donbas geological and industrial region is characterised by complex mining, geological, and hydrogeological conditions. The structure of the region’s productive strata contains a large number of aquifers, which leads to wetting and collapse of the surrounding rocks, threatens safety, and increases the complexity of excavations. The effectiveness of existing methods to combat water inflows is insufficient, as they bring about large financial costs and do not take into account changes in the stress–strain state of rocks around the underground space and its impact on the permeability of the massif [78,79].
Our goal was to prevent water seepage into underground spaces in the conditions of the Samarska coal mine. Underground workings are used to transport coal and auxiliary materials. However, a complicating factor is the location of the mine in the area of the Bogdanovsky dump geological disturbance. This results in frequent water breakthroughs into the underground space.
Surface and groundwaters are located within the mine area. The main water inflows come from the widespread Buchak aquifer, which is composed of lightly clayey sands. Water from overburden deposits does not contribute to the watering of the mine workings. In the coal deposits, the water-bearing layers are sandstones and limestones, the thickness of which varies from 7–10 m to 100.7 m. The rocks of the upper part of the strata have increased water-bearing capacity. The water-bearing capacity of rocks in tectonic fault zones does not differ from the water-bearing capacity of faulted water-bearing rocks. However, this does not exclude the temporary inflow of water to mine workings when they are close to tectonic faults.
The average annual water inflow to the mine is 41 m3/h, and the salinity of mine water is 20–35 g/L. The water is corrosive, with a large amount of solid sediment, and has a pH of 6.5–6.8. The water hardness ranges from 21 to 40 mg equivalent per litre.
Therefore, our task was to choose the most rational technology for the construction of workings to prevent water leakage in the underground space.
Regardless of the support material (frame-load-bearing, sprayed concrete, anchoring, etc.), we found that there were common stages that involved solving a number of local problems, and we will focus on them below (Table 3).
Table 3 shows that regardless of the type of support, construction technology, or materials used, all stages are identical (I–VII). This allowed us to formalise the task of selecting the technology for the construction of workings. To do this, it was necessary to represent the cycle of roofing the excavated space as a network model. The vertices of the network model corresponded to certain stages, and the distances between the vertices were the values of the optimisation parameter (e.g., cost, time, quantity of materials, etc.).
Let us depict a network model of the process of the working face construction in Figure 5.
As previously mentioned, the vertices are taken as technology options, and the distance between the vertices is the value of the optimisation parameter. Table 4 provides an explanation of the network model shown in Figure 5.
Then, to find a solution that corresponds to the optimal technology for the construction of workings, it came down to finding the shortest route from vertex 1 (the start vertex required to execute the algorithm) to vertex 25 (the end vertex required to complete the algorithm). With the help of optimisation algorithms, the route could be found automatically. In Figure 5, for clarity, the optimal route is highlighted with a thickened line; it is the shortest and optimal route. As noted earlier, only actual links were analysed.
Based on the results of the study, we proposed the use of anchoring [80,81]. Based on studies of the geomechanical parameters of the rock mass, it was found that the use of anchoring would reduce the value of the heterogeneity by 25–30% and the permeability in the roof of the mine workings by 25–35% when tampering with the coal aquifer. Anchoring prevents the development of fissures, preserves the host rocks in a natural, monolithic state, and increases the stability of the mine workings, even if there is a waterlogged coal seam in the roof. A significant reduction in the size of the filtration area and a decrease in permeability within it leads to a decrease in the intensity of the filtration movement of fluid and to the prevention or reduction of water inflow into the anchored workings, which makes it possible to use anchoring as a technological method to reduce water inflow into the mine workings [82,83,84]. Figure 6 shows the results of the implemented technological solution to combat the negative effects of water in underground spaces.
Figure 7 shows a close-up of the technological solution. Also, Figure 7 shows that there are no signs of water in the area of the geological fault according to the observations.
The considered basic model allowed us to choose the optimal technological solutions to combat water inflows in the underground space. However, there are situations when hydrogeological conditions are unfavourable and there are risks associated with water breakthroughs into underground spaces. Therefore, technological solutions should be provided to minimise these risks.

3.2. Model for Finding Optimal Technological Solutions to Combat Water Inflows in Difficult Hydrogeological Conditions

The procedure is identical to the one described in Section 3.1. However, measures related to the prevention of negative phenomena associated with water inflows are taken into account.
Table 5 shows the main stages of the support construction and describes possible hydrogeological factors that affect the efficiency of the process. Based on the analysis of these factors, additional technological stages are envisaged to increase the stability of the rock mass. Knowledge of the technologies make it possible to develop a network model of alternatives to support work under difficult hydrogeological conditions.
Table 5 shows that all the stages are universal, i.e., they do not depend on the technology and type of support used. This allowed us to find the optimal technology for controlling water inflows into the underground space based on dynamic programming.
Let us depict in Figure 8 a network model of the process of building a working face. It should be noted that Figure 8a shows the basic variant and Figure 8b a variant taking into account complex hydrogeological conditions. We note that the consideration of complex hydrogeological conditions can be taken into account in different ways: Figure 8b shows a simplified version, when the network model is corrected (loaded) by a certain coefficient equal to the probability of occurrence of an unfavourable hydrogeological factor.
For example, in the baseline case, the preparatory work technology (7–10) corresponds to a cost value of 10 c.u. However, we know that there are additional costs for the application of surfactants and possibly re-drilling holes; these additional costs are 50% for this type of support, so the value of the optimisation parameter for this technology should be multiplied by 1.5.
As noted earlier, the vertices are taken as technology options, and the distance between the vertices is the value of the optimisation parameter.
It should be noted that the accounting of amendments can be applied for express analysis, even with limited resources, but the design process can be carried out more accurately by introducing additional stages to the basic model (Figure 8a) to combat the negative impact of various factors, e.g., strengthening the massif, the construction of drainage grooves, etc. Figure 9 shows a parameter optimisation model which takes into account complex hydrogeological conditions. It is worth noting that by neutralising the negative impact of water on the massif condition, it is possible to minimise costs in subsequent stages, i.e., preliminary measures (technology option) affect the efficiency of the process as a whole.
As shown in Figure 9, the number of main stages is identical, but in difficult hydrogeological conditions, it is necessary to provide intermediate stages which are directly related to measures to strengthen the state of the rock mass. If the provided technology is able to control water inflows, then an additional stage is not provided; otherwise, it should be provided (Table 6).
Table 6 shows that to find the optimal solution, we need to find the distance from node 1 to node 34. The optimal solution is the one that corresponds to the shortest path in the network model (Figure 9). A prerequisite is that there is a connection between the technological stages.
Thus, if it is known that the array is waterlogged, then certain links in the network model should be strengthened, which will result in an increase in costs, while others should be weakened.

4. Discussion

In the present study, we proposed the use of dynamic programming methods to improve technological solutions to combat water inflows into underground spaces. Based on our analysis of existing approaches to finding optimal solutions, we found that, currently, the most advanced methods are those that break down the overall problem into partial ones, whereby it is possible to take into account a variety of factors and increase the accuracy of the solution.
Dynamic programming methods can be used to obtain optimal solutions in each stage, which contributes to the maximum probability of achieving the goal. In accordance with the proposed solution, we propose an algorithm with which to find the optimal technological solution to improve the methods of combating negative water leakage in underground spaces. According to this algorithm, it is first necessary to divide the general task of substantiating the technology for combating water leakage into stages that correspond to the search for a specific solution (transport, fastening, dismantling, etc.). After that, alternative solutions are presented in the form of a network model. Optimisation algorithms are used to find the shortest route in the network model that corresponds to the optimal solution.
If there are risks associated with unfavourable hydrogeological conditions, the proposed model should be adjusted. The number of basic stages is identical, but intermediate stages are introduced, which include measures related to the minimisation of hydrogeological risks. However, the construction of a refined network model involves certain time costs. To prevent this, we propose the use of a “simplified” network model, which is identical to the basic one but which makes allowances for negative hydrogeological factors (e.g., water breakthrough into the underground space).
Our approach is universal. It can be applied not only in industrial underground structures but also in civil engineering (subways, tunnels). As an example of this approach, we analysed the conditions in an industrial enterprise—a coal mine. However, this was no coincidence, as coal mines often face risks of water leakage for which there are a large number of alternative solutions, and the structures have a long service life (from 20 to 100 years). This variability allowed us to investigate the effectiveness of the proposed solutions. For the conditions of an enterprise that operates within a geological fault, the use of anchoring is proposed, as this reduces the permeability of roof rocks by 25–30%. According to the results of field observations, no water seepage into the underground space was recorded.
In conclusion, it is worth focusing on further areas of research. Our approaches are basic. We should additionally focus on the search for technologies to combat water inflows that involve the use of industrial waste.
The above models and algorithms can be implemented in software. Therefore, the task in the short term is to create a decision support system in the form of software. This will create a bank of design solutions and significantly increase the dimensionality of the tasks.

5. Conclusions

The utilisation of dynamic programming models represents a significant advancement in the domain of combating negative water leakage in underground spaces. Through our study, we have proposed the application of dynamic programming methods to enhance existing technological solutions, aiming to address water inflows effectively.
To operationalise this approach, we have formulated an algorithm for identifying the optimal technological solution to combat negative water leakage. This algorithm involves segmenting the task into stages corresponding to specific solution components and utilising optimisation algorithms to identify the shortest route in the network model, corresponding to the optimal solution. In scenarios with risks associated with unfavourable hydrogeological conditions, our model can be adjusted to incorporate intermediate stages focused on risk mitigation measures. While constructing refined network models incurs time costs, the proposal of a “simplified” version accounts for negative hydrogeological factors, offering a practical alternative without compromising effectiveness.
By analysing the conditions of a coal mine as an illustrative example, we have showcased the effectiveness of our approach. For instance, in enterprises operating within geological fault zones, anchoring techniques have proven to be effective in reducing roof rock permeability, consequently preventing water seepage into the underground space.
In conclusion, our study underscores the potential of dynamic programming models to improve technological approaches for addressing negative water leakage in underground spaces. By leveraging these advancements and continuing research efforts, we can enhance resilience against water-related challenges within various engineering contexts.

Author Contributions

Conceptualization, S.L., S.D. and A.K.; methodology, A.K.; software, S.L. and S.D.; validation, S.D. and A.K.; formal analysis, S.L. and S.D.; investigation, S.L. and S.D.; resources, S.D. and A.K.; data curation, S.L. and A.K.; writing—original draft preparation, S.L. and A.K.; writing—review and editing, S.L. and A.K.; visualisation, S.D. and A.K.; supervision, S.L., S.D. and A.K.; project administration, S.L. and S.D.; funding acquisition, S.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Nate, S.; Bilan, Y.; Cherevatskyi, D.; Kharlamova, G.; Lyakh, O.; Wosiak, A. The Impact of Energy Consumption on the Three Pillars of Sustainable Development. Energies 2021, 14, 1372. [Google Scholar] [CrossRef]
  2. Hrinov, V.; Khorolskyi, A. Improving the Process of Coal Extraction Based on the Parameter Optimization of Mining Equipment. E3S Web Conf. 2018, 60, 00017. [Google Scholar] [CrossRef]
  3. Shepherd, N.G.; Rudd, J.M. The Influence of Context on the Strategic Decision-Making Process: A Review of the Literature. Int. J. Manag. Rev. 2014, 16, 340–364. [Google Scholar] [CrossRef]
  4. Stankevich, A. Explaining the Consumer Decision-Making Process: Critical Literature Review. J. Int. Bus. Res. Mark. 2017, 2, 7–14. [Google Scholar] [CrossRef]
  5. Khorolskyi, A.; Hrinov, V.; Mamaikin, O.; Demchenko, Y. Models and Methods to Make Decisions While Mining Production Scheduling. Min. Miner. Depos. 2019, 13, 53–62. [Google Scholar] [CrossRef]
  6. Li, X.; Chen, G.; Dong, Z.Y. Co-optimisation Model for the Long-term Design and Decision Making in Community Level Cloud Energy Storage System. IET Renew. Power Gener. 2020, 14, 3518–3525. [Google Scholar] [CrossRef]
  7. Huang, S.; Li, G.; Ben-Awuah, E.; Afum, B.O.; Hu, N. A Robust Mixed Integer Linear Programming Framework for Underground Cut-and-Fill Mining Production Scheduling. Int. J. Min. Reclam. Environ. 2020, 34, 397–414. [Google Scholar] [CrossRef]
  8. Gill, P.E.; Wong, E. Methods for Convex and General Quadratic Programming. Math. Prog. Comp. 2015, 7, 71–112. [Google Scholar] [CrossRef]
  9. Little, J.; Topal, E. Strategies to Assist in Obtaining an Optimal Solution for an Underground Mine Planning Problem Using Mixed Integer Programming. Int. J. Min. Miner. Eng. 2011, 3, 152. [Google Scholar] [CrossRef]
  10. MacNeil, J.A.L.; Dimitrakopoulos, R.G. A Stochastic Optimization Formulation for the Transition from Open Pit to Underground Mining. Optim. Eng. 2017, 18, 793–813. [Google Scholar] [CrossRef]
  11. Khorolskyi, A.; Hrinov, V.; Kaliushenko, O. Network Models for Searching for Optimal Economic and Environmental Strategies for Field Development. Procedia Environ. Sci. Eng. Manag. 2019, 6, 463–471. [Google Scholar]
  12. Balezentis, T.; Streimikiene, D.; Siksnelyte-Butkiene, I. Energy Storage Selection for Sustainable Energy Development: The Multi-Criteria Utility Analysis Based on the Ideal Solutions and Integer Geometric Programming for Coordination Degree. Environ. Impact Assess. Rev. 2021, 91, 106675. [Google Scholar] [CrossRef]
  13. Hill, A.; Brickey, A.J.; Cipriano, I.; Goycoolea, M.; Newman, A. Optimization Strategies for Resource-Constrained Project Scheduling Problems in Underground Mining. INFORMS J. Comput. 2022, 34, 3042–3058. [Google Scholar] [CrossRef]
  14. Özyurt, M.C.; Karadogan, A. A New Model Based on Artificial Neural Networks and Game Theory for the Selection of Underground Mining Method. J. Min. Sci. 2020, 56, 66–78. [Google Scholar] [CrossRef]
  15. Khorolskyi, A.; Hrinov, V.; Mamaikin, O.; Fomychova, L. Research into Optimization Model for Balancing the Technological Flows at Mining Enterprises. E3S Web Conf. 2020, 201, 01030. [Google Scholar] [CrossRef]
  16. Petlovanyi, M.; Sai, K.; Malashkevych, D.; Popovych, V.; Khorolskyi, A. Influence of Waste Rock Dump Placement on the Geomechanical State of Underground Mine Workings. IOP Conf. Ser. Earth Environ. Sci. 2023, 1156, 012007. [Google Scholar] [CrossRef]
  17. Kwinta, A.; Gradka, R. Analysis of the Damage Influence Range Generated by Underground Mining. Int. J. Rock Mech. Min. Sci. 2020, 128, 104263. [Google Scholar] [CrossRef]
  18. Vayenas, N.; Peng, S. Reliability Analysis of Underground Mining Equipment Using Genetic Algorithms: A Case Study of Two Mine Hoists. J. Qual. Maint. Eng. 2014, 20, 32–50. [Google Scholar] [CrossRef]
  19. Thoft-Christensen, P.; Murotsu, Y. Application of Structural Systems Reliability Theory; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  20. Petropoulos, F.; Apiletti, D.; Assimakopoulos, V.; Babai, M.Z.; Barrow, D.K.; Ben Taieb, S.; Bergmeir, C.; Bessa, R.J.; Bijak, J.; Boylan, J.E.; et al. Forecasting: Theory and Practice. Int. J. Forecast. 2022, 38, 705–871. [Google Scholar] [CrossRef]
  21. Tsai, S.-B.; Xue, Y.; Zhang, J.; Chen, Q.; Liu, Y.; Zhou, J.; Dong, W. Models for Forecasting Growth Trends in Renewable Energy. Renew. Sustain. Energy Rev. 2017, 77, 1169–1178. [Google Scholar] [CrossRef]
  22. Daganzo, C. Multinomial Probit: The Theory and Its Application to Demand Forecasting; Elsevier: Amsterdam, The Netherlands, 2014. [Google Scholar]
  23. Lousada, S.A.N. Hydraulics: Theory—Volume I; Universidade Da Madeira: Funchal, Portugal, 2018; p. 348. ISBN 978-989-8805-48-5. [Google Scholar]
  24. Lane, S.N. Natural Flood Management. WIREs Water 2017, 4, e1211. [Google Scholar] [CrossRef]
  25. Blagojević, B.; Vasilevska, L.; Anđelković, D.; Bogojević, A.; Lousada, S. A Framework for Assessing Nature-Based Urban Stormwater Management Solutions: A Preliminary Spatial Analysis Approach Applied to Southeast Serbia. Water 2023, 15, 3604. [Google Scholar] [CrossRef]
  26. Jenkins, K.; Surminski, S.; Hall, J.; Crick, F. Assessing Surface Water Flood Risk and Management Strategies under Future Climate Change: Insights from an Agent-Based Model. Sci. Total Environ. 2017, 595, 159–168. [Google Scholar] [CrossRef] [PubMed]
  27. Lousada, S. Geoinformatics in Support of Urban Politics and the Development of Civil Engineering; IGI Global: Hershey, PA, USA, 2023. [Google Scholar]
  28. Li, W.; De Vriend, H.J.; Wang, Z.; Van Maren, D.S. Morphological Modeling Using a Fully Coupled, Total Variation Diminishing Upwind-Biased Centered Scheme: Efficient and Accurate Coupled Morphological Modeling. Water Resour. Res. 2013, 49, 3547–3565. [Google Scholar] [CrossRef]
  29. Le, M.-H.; Cordier, S.; Lucas, C.; Cerdan, O. A Faster Numerical Scheme for a Coupled System Modeling Soil Erosion and Sediment Transport. Water Resour. Res. 2015, 51, 987–1005. [Google Scholar] [CrossRef]
  30. Gunantara, N. A Review of Multi-Objective Optimization: Methods and Its Applications. Cogent Eng. 2018, 5, 1502242. [Google Scholar] [CrossRef]
  31. Bottou, L.; Curtis, F.E.; Nocedal, J. Optimization Methods for Large-Scale Machine Learning. SIAM Rev. 2018, 60, 223–311. [Google Scholar] [CrossRef]
  32. Ball, M.O. Heuristics Based on Mathematical Programming. Surv. Oper. Res. Manag. Sci. 2011, 16, 21–38. [Google Scholar] [CrossRef]
  33. Pishvaee, M.S.; Khalaf, M.F. Novel Robust Fuzzy Mathematical Programming Methods. Appl. Math. Model. 2016, 40, 407–418. [Google Scholar] [CrossRef]
  34. Grossmann, I.E. Advances in Mathematical Programming Models for Enterprise-Wide Optimization. Comput. Chem. Eng. 2012, 47, 2–18. [Google Scholar] [CrossRef]
  35. Xu, G.; Wan, S.; Wang, F.; Dong, J.; Zeng, Y. Mathematical Programming Methods for Consistency and Consensus in Group Decision Making with Intuitionistic Fuzzy Preference Relations. Knowl.-Based Syst. 2016, 98, 30–43. [Google Scholar] [CrossRef]
  36. Herskovits, J.; Mappa, P.; Goulart, E.; Mota Soares, C.M. Mathematical Programming Models and Algorithms for Engineering Design Optimization. Comput. Methods Appl. Mech. Eng. 2005, 194, 3244–3268. [Google Scholar] [CrossRef]
  37. D’Ambrosio, C.; Lodi, A.; Wiese, S.; Bragalli, C. Mathematical Programming Techniques in Water Network Optimization. Eur. J. Oper. Res. 2015, 243, 774–788. [Google Scholar] [CrossRef]
  38. Yang, L.; Liu, S.; Tsoka, S.; Papageorgiou, L.G. Mathematical Programming for Piecewise Linear Regression Analysis. Expert Syst. Appl. 2016, 44, 156–167. [Google Scholar] [CrossRef]
  39. Saaty, T.L.; Vargas, L.G. The Analytic Network Process. In Decision Making with the Analytic Network Process; International Series in Operations Research & Management Science; Springer: Boston, MA, USA, 2013; Volume 195, pp. 1–40. [Google Scholar] [CrossRef]
  40. Brans, J.-P.; De Smet, Y. PROMETHEE Methods. In Multiple Criteria Decision Analysis; Greco, S., Ehrgott, M., Figueira, J.R., Eds.; International Series in Operations Research & Management Science; Springer: New York, NY, USA, 2016; Volume 233, pp. 187–219. [Google Scholar] [CrossRef]
  41. Corrente, S.; Figueira, J.R.; Greco, S. The SMAA-PROMETHEE Method. Eur. J. Oper. Res. 2014, 239, 514–522. [Google Scholar] [CrossRef]
  42. Singh, A.; Gupta, A.; Mehra, A. Best Criteria Selection Based PROMETHEE II Method. Opsearch 2021, 58, 160–180. [Google Scholar] [CrossRef]
  43. Ishizaka, A.; Nemery, P. Selecting the Best Statistical Distribution with PROMETHEE and GAIA. Comput. Ind. Eng. 2011, 61, 958–969. [Google Scholar] [CrossRef]
  44. Figueira, J.R.; Mousseau, V.; Roy, B. ELECTRE Methods. In Multiple Criteria Decision Analysis; Greco, S., Ehrgott, M., Figueira, J.R., Eds.; International Series in Operations Research & Management Science; Springer: New York, NY, USA, 2016; Volume 233, pp. 155–185. [Google Scholar] [CrossRef]
  45. Llamazares, B. An Analysis of the Generalized TODIM Method. Eur. J. Oper. Res. 2018, 269, 1041–1049. [Google Scholar] [CrossRef]
  46. Alali, F.; Tolga, A.C. Portfolio Allocation with the TODIM Method. Expert Syst. Appl. 2019, 124, 341–348. [Google Scholar] [CrossRef]
  47. Taherdoost, H.; Madanchian, M. An Effective Compromising Ranking Technique for Decision Making. Macro Manag. Public Policies 2023, 5, 27–33. [Google Scholar] [CrossRef]
  48. Ayağ, Z.; Özdemir, R.G. A Fuzzy AHP Approach to Evaluating Machine Tool Alternatives. J. Intell. Manuf. 2006, 17, 179–190. [Google Scholar] [CrossRef]
  49. Srdjevic, B.; Medeiros, Y.D.P. Fuzzy AHP Assessment of Water Management Plans. Water Resour. Manag. 2008, 22, 877–894. [Google Scholar] [CrossRef]
  50. Sahoo, S.; Dhar, A.; Kar, A. Environmental Vulnerability Assessment Using Grey Analytic Hierarchy Process Based Model. Environ. Impact Assess. Rev. 2016, 56, 145–154. [Google Scholar] [CrossRef]
  51. Roozbahani, A.; Ghased, H.; Shahedany, M.H. Inter-Basin Water Transfer Planning with Grey COPRAS and Fuzzy COPRAS Techniques: A Case Study in Iranian Central Plateau. Sci. Total Environ. 2020, 726, 138499. [Google Scholar] [CrossRef]
  52. Liu, Y.; Wang, C.; Chun, Y.; Yang, L.; Chen, W.; Ding, J. A Novel Method in Surface Water Quality Assessment Based on Improved Variable Fuzzy Set Pair Analysis. Int. J. Environ. Res. Public Health 2019, 16, 4314. [Google Scholar] [CrossRef]
  53. Yan, F.; Liu, L.; Zhang, Y.; Chen, M.; Chen, N. The Research of Dynamic Variable Fuzzy Set Assessment Model in Water Quality Evaluation. Water Resour. Manag. 2016, 30, 63–78. [Google Scholar] [CrossRef]
  54. Saaty, T.L. Decision Making with the Analytic Hierarchy Process. Int. J. Serv. Sci. 2008, 1, 83. [Google Scholar] [CrossRef]
  55. Bildirici, M.E.; Lousada, S.; Yılmaz Genç, S. Terrorism, Freshwater, and Environmental Pollution: Evidence of Afghanistan, Burkina Faso, Iraq, Arab Republic of Egypt, Cameroon, Mali, Mozambique, Niger, Nigeria, Somalia, Syrian Arab Republic, and Pakistan. Water 2022, 14, 2684. [Google Scholar] [CrossRef]
  56. Lousada, S.A.N. Hydrology, Water Resources and Environment: Theorical Classes; Universidade Da Madeira: Funchal, Portugal, 2019; pp. 1–220. [Google Scholar]
  57. Paravarzar, S.; Pourrahimian, Y.; Nasab, H.A.; Emery, X. Short-Term Underground Mine Planning: A Review. Int. J. Min. Miner. Eng. 2021, 12, 1. [Google Scholar] [CrossRef]
  58. Wang, K.; Davies, E.G.R. Municipal Water Planning and Management with an End-Use Based Simulation Model. Environ. Model. Softw. 2018, 101, 204–217. [Google Scholar] [CrossRef]
  59. Babamiri, O.; Marofi, S. A Multi-Objective Simulation–Optimization Approach for Water Resource Planning of Reservoir–River Systems Based on a Coupled Quantity–Quality Model. Environ. Earth Sci. 2021, 80, 389. [Google Scholar] [CrossRef]
  60. Singh, A. Simulation–Optimization Modeling for Conjunctive Water Use Management. Agric. Water Manag. 2014, 141, 23–29. [Google Scholar] [CrossRef]
  61. Grundmann, J.; Schütze, N.; Schmitz, G.H.; Al-Shaqsi, S. Towards an Integrated Arid Zone Water Management Using Simulation-Based Optimisation. Environ. Earth Sci. 2012, 65, 1381–1394. [Google Scholar] [CrossRef]
  62. Reinhart, R.; Dang, T.; Hand, E.; Papachristos, C.; Alexis, K. Learning-Based Path Planning for Autonomous Exploration of Subterranean Environments. In Proceedings of the 2020 IEEE International Conference on Robotics and Automation (ICRA), Paris, France, 31 May–31 August 2020; pp. 1215–1221. [Google Scholar] [CrossRef]
  63. Li, S.; Huang, Q.; Hu, B.; Pan, J.; Chen, J.; Yang, J.; Zhou, X.; Wang, X.; Yu, H. Mining Method Optimization of Difficult-to-Mine Complicated Orebody Using Pythagorean Fuzzy Sets and TOPSIS Method. Sustainability 2023, 15, 3692. [Google Scholar] [CrossRef]
  64. Emdini Gliwan, S.; Crowe, K. A Network Flow Model for Operational Planning in an Underground Gold Mine. Mining 2022, 2, 712–724. [Google Scholar] [CrossRef]
  65. Fazlali, A.; Shourian, M. A Demand Management Based Crop and Irrigation Planning Using the Simulation-Optimization Approach. Water Resour. Manag. 2018, 32, 67–81. [Google Scholar] [CrossRef]
  66. Delehan, S.; Melehanych, H.; Khorolskyi, A. The Traditions and Technologies of Ecological Construction in Portugal. In Proceedings of the 4th International Conference on Advances in Environmental Engineering, Ostrava, Czech Republic, 20–22 November 2023; MDPI: Basel, Switzerland, 2023; p. 23. [Google Scholar] [CrossRef]
  67. Bhatt, D.; Mall, R.K. Surface Water Resources, Climate Change and Simulation Modeling. Aquat. Procedia 2015, 4, 730–738. [Google Scholar] [CrossRef]
  68. Soleimanian, E.; Afshar, A.; Molajou, A. A Review on Water Simulation Models for the WEF Nexus: Development Perspective. Environ. Sci. Pollut. Res. 2022, 29, 79769–79785. [Google Scholar] [CrossRef] [PubMed]
  69. Salli, S.; Pochepov, V.; Mamaykin, O. Theoretical Aspects of the Potential Technological Schemes Evaluation and Their Susceptibility to Innovations. In Progressive Technologies of Coal, Coalbed Methane, and Ores Mining; CRC Press: London, UK, 2014; pp. 479–484. [Google Scholar] [CrossRef]
  70. Guo, D.; Hou, H.; Long, J.; Guo, X.; Xu, H. Underestimated Environmental Benefits of Tailings Resource Utilization: Evidence from a Life Cycle Perspective. Environ. Impact Assess. Rev. 2022, 96, 106832. [Google Scholar] [CrossRef]
  71. Eddy, S.R. What Is Dynamic Programming? Nat. Biotechnol. 2004, 22, 909–910. [Google Scholar] [CrossRef]
  72. Fu, M.; Liu, S.; Jia, H.; He, D. Experimental Study of an Orientation and Resin-Lifting Device for Improving the Performance of Resin-Anchored Roof Bolts. Rock Mech. Rock Eng. 2020, 53, 211–231. [Google Scholar] [CrossRef]
  73. Bondarenko, V.I.; Kovalevska, I.A.; Symanovych, H.A.; Sachko, R.M.; Sheka, I.V. Integrated Research into the Stress-Strain State Anomalies, Formed and Developed in the Mass under Conditions of High Advance Velocities of Stope Faces. IOP Conf. Ser. Earth Environ. Sci. 2023, 1254, 012062. [Google Scholar] [CrossRef]
  74. Bondarenko, V.I.; Kovalevska, I.A.; Podkopaiev, S.V.; Sheka, I.V.; Tsivka, Y.S. Substantiating Arched Support Made of Composite Materials (Carbon Fiber-Reinforced Plastic) for Mine Workings in Coal Mines. IOP Conf. Ser. Earth Environ. Sci. 2022, 1049, 012026. [Google Scholar] [CrossRef]
  75. Bellman, R.E. Dynamic Programming; Princeton University Press: Princeton, NJ, USA, 2010. [Google Scholar] [CrossRef]
  76. He, G.; Zhang, X.; Sun, Y.; Luo, G.; Chen, L. Contour Line Simplification Method Based on the Two-level Bellman–Ford Algorithm. Trans. GIS 2021, 25, 396–423. [Google Scholar] [CrossRef]
  77. Hemalatha, S.; Valsalal, P. Identification of Optimal Path in Power System Network Using Bellman Ford Algorithm. Model. Simul. Eng. 2012, 2012, 913485. [Google Scholar] [CrossRef]
  78. Małkowski, P.; Niedbalski, Z.; Majcherczyk, T.; Bednarek, Ł. Underground Monitoring as the Best Way of Roadways Support Design Validation in a Long Time Period. Min. Miner. Depos. 2020, 14, 1–14. [Google Scholar] [CrossRef]
  79. Krukovskyi, O.; Bulich, Y.; Kurnosov, S.; Yanzhula, O.; Demin, V. Substantiating the Parameters for Selecting a Pillar Width to Protect Permanent Mine Workings at Great Depths. IOP Conf. Ser. Earth Environ. Sci. 2022, 970, 012049. [Google Scholar] [CrossRef]
  80. Małkowski, P.; Niedbalski, Z.; Balarabe, T. A Statistical Analysis of Geomechanical Data and Its Effect on Rock Mass Numerical Modeling: A Case Study. Int. J. Coal Sci. Technol. 2021, 8, 312–323. [Google Scholar] [CrossRef]
  81. Sdvyzhkova, O.; Patyńska, R. Effect of Increasing Mining Rate on Longwall Coal Mining—Western Donbass Case Study. Stud. Geotech. Mech. 2016, 38, 91–98. [Google Scholar] [CrossRef]
  82. Skipochka, S.; Krukovskyi, O.; Palamarchuk, T.; Prokhorets, L. On the Methodology for Considering Scale Effect of Rock Strength. Min. Miner. Depos. 2020, 14, 24–30. [Google Scholar] [CrossRef]
  83. Kochmar, I.; Karabyn, V. Water Extracts from Waste Rocks of the Coal Industry of Chernvonograd Mining Area (Ukraine): Problems of Environmental Safety and Civil Protection. Ecol. Eng. Environ. Technol. 2023, 24, 247–255. [Google Scholar] [CrossRef]
  84. Huang, Z.; Cheng, Y.; Zhang, D.; Yan, D.; Shen, Y. Seismic Fragility and Resilience Assessment of Shallowly Buried Large-Section Underground Civil Defense Structure in Soft Soils: Framework and Application. Tunn. Undergr. Space Technol. 2024, 146, 105640. [Google Scholar] [CrossRef]
Figure 1. Diagram of the procedure for finding an optimal solution.
Figure 1. Diagram of the procedure for finding an optimal solution.
Water 16 01952 g001
Figure 3. Network model for finding the optimal type of fastening that will prevent water inflow into the underground space.
Figure 3. Network model for finding the optimal type of fastening that will prevent water inflow into the underground space.
Water 16 01952 g003
Figure 4. Network model for finding the optimal type of support with consideration of hydrogeological factors.
Figure 4. Network model for finding the optimal type of support with consideration of hydrogeological factors.
Water 16 01952 g004
Figure 5. Network model for optimising the parameters of mine construction.
Figure 5. Network model for optimising the parameters of mine construction.
Water 16 01952 g005
Figure 6. Absence of water leakage in the underground space of the Samarska mine in the transition zone of the high-amplitude geological fault Bogdanovsky discharge (based on the results of the recommended technology).
Figure 6. Absence of water leakage in the underground space of the Samarska mine in the transition zone of the high-amplitude geological fault Bogdanovsky discharge (based on the results of the recommended technology).
Water 16 01952 g006
Figure 7. The state of the contour rocks in the underground space of the Samarska mine in the transition zone of the high-amplitude geological fault Bogdanovsky discharge.
Figure 7. The state of the contour rocks in the underground space of the Samarska mine in the transition zone of the high-amplitude geological fault Bogdanovsky discharge.
Water 16 01952 g007
Figure 8. Network model for the optimisation of mine construction parameters: (a) basic model; (b) adjusted for complex hydrogeological conditions.
Figure 8. Network model for the optimisation of mine construction parameters: (a) basic model; (b) adjusted for complex hydrogeological conditions.
Water 16 01952 g008
Figure 9. Network model for optimising the parameters of excavation construction in complex hydrogeological conditions.
Figure 9. Network model for optimising the parameters of excavation construction in complex hydrogeological conditions.
Water 16 01952 g009
Table 1. Relationship between hydrogeological and technological factors in the construction and operation of workings.
Table 1. Relationship between hydrogeological and technological factors in the construction and operation of workings.
Negative Impact of GroundwaterHydrogeological FactorTechnological Factor and Risk
reduced productivity specific water inflowreduced equipment performance
the degree of watering of rocks and water permeabilityincreased construction time and costs
groundwater head, aquifer thickness, water permeabilityflooding of workings, additional costs for drainage
reduced stability of productsfiltration, degree of water saturationadditional costs for the support of the workings, measures to combat the heaving of the sole of the workings
reduced labour safetywater inflow, water permeability, filtrationdestruction of temporary or permanent fasteners
deterioration in the operation of transport equipmentwater saturationcomplicated transport of materials, increased construction time
Table 2. Costs for the construction and operation of the support.
Table 2. Costs for the construction and operation of the support.
Cost Component of the Total CostConstruction and Operation PhaseRisks Associated with Water Inflows
materials from which the fasteners will be made to prevent water leakageinitial stageNone
transport and warehousing costsloading, transport and unloading of fastenersare present, it is necessary to take into account water permeability
tools, time costspreparatory stages (drilling holes for anchors, etc.), transport of materialsare present, it is necessary to have an idea of the thickness of the aquifers of the massif
materials, labour intensityconstruction of temporary supportare present, it is necessary to analyse the filtering
maintenance costs, labour intensityconstruction of a permanent support and further operationare present, the water inflow into the workings is analysed
Table 3. Formalisation of the task of erecting a support in the underground space to prevent water leakage.
Table 3. Formalisation of the task of erecting a support in the underground space to prevent water leakage.
StageStage NameObjectivesThe Essence of the Task Solution
Iorganisation of loading and unloading operations on the surfacebuild a rational structure of the cycle for the organisation of surface loading operationsa balance should be found between resource flows to speed up cargo operations
IIdelivery of fastening materials in the underground space, taking into account the cost of storing materials in the underground spacea rational structure of the transport chain should be built with the condition of maximum preservation of the quality and quantity of materials for the construction of workingsa balance should be struck between transport operations, suppliers (warehouses) and delivery to the site of construction
IIIcarrying out preparatory work required for the construction of temporary support (drilling holes, constructing grooves, stripping, etc.)minimising time spent on preparatory operations and reducing costsit is necessary to select drilling, assembly and loading equipment for these operations
IVerection of a temporary fixtureorganisation of the technology for the construction of workings and temporary consolidation of the excavated spaceorganisation of work on securing the excavated space. main requirements: ease of dismantling, minimisation of the amount of fastening materials, possibility of removing waste generated as a result of previous operations (I–III)
Vcarrying out work on the construction of a permanent supportminimisation of the amount of materials required to securely fix the produced spaceorganisation of works on the construction of permanent support. it is necessary to minimise equipment downtime and reduce the consumption of fastening materials. at the same time, it is necessary to ensure the appropriate reliability of the support
VIcleaning the underground space from the rock mass formed as a result of the construction of the supportminimisation of manual labour, reduction of wasteit is necessary to organise a technology that will keep the maximum amount of waste in the underground space. at the same time, the number of human resources involved in cleaning the produced space should be minimised
VIIdismantling of fasteners in the underground spaceorganisation of technology with the least time and financial coststhe technology of dismantling the support in the underground space should be organised. safe working conditions should be ensured. the number of human resources involved in these operations should be minimal
Table 4. Network formalisation of the problem of mine construction.
Table 4. Network formalisation of the problem of mine construction.
Stage Designation in Figure 5Stage NameStarting PointFinal Peak (Peaks at Intermediate Stages)Interpretation
Iloading of fasteners on the surface12–41—the starting point; 2–4 options for transport technology
IIdelivery of fasteners in the underground space, taking into account warehouse costs25–82—the best option after the 1st stage; 5–8 tops corresponding to the transport technologies
IIIpreparatory work69–116—optimal solution after 2 stages; 9–11 variants of borehole drilling technologies
IVconstruction of temporary support912–159—optimal technology after three stages; 12–15 variants of temporary support construction technology
Vconstruction of permanent support1216–1712—optimal technology after four stages; 16—fastening technologies
VIremoval and transport of waste generated during construction1618–2016—optimal solution after 5 stages; 18–20 technologies for transporting production waste
VIIremoval of fastening materials1821–2418—optimal technology after six stages; 21–24 dismantling technologies
VIIIcompletion242521—optimal technology after completion of seven stages; 25—completion of the cycle for selecting the optimal technology
Table 5. Stages of roof support construction, taking into account hydrogeological conditions.
Table 5. Stages of roof support construction, taking into account hydrogeological conditions.
StageStage NameComplicating Hydrogeological ConditionsProcedure for Taking into Account the Hydrogeological Factor at the Design Stage
FactorHydrogeological
Parameter
Iloading of fasteners on the surfaceNoneNonenone
IIdelivery of fasteners in products, taking into account warehouse costsmaterial sticking to the delivery meansfiltration coefficientit is necessary to allow for additional downtime for cleaning the delivery vehicles (as a result, the speed of the workings is recalculated), or an additional stage for strengthening the rock mass
caking of materials
IIIpreparatory work (drilling holes, etc.)destruction of boreholes
filtration coefficient, degree of watering of rocksit is necessary to include additional stages for the treatment of the contact surfaces of the array, the application of surface-active substances (surfactants)
unsatisfactory working conditions (waterlogged workings)location and characteristics of waterstops, degree of watering of rocks
IVerection of temporary supportunsatisfactory characteristics of the solution, destruction of boreholeswater absorption coefficient, air absorption coefficientadditional funds should be provided for strengthening the rock mass, additional costs for drainage ditches
Vconstruction of permanent supportunsatisfactory characteristics of the solution, unfavourable working conditionsilocation and characteristics of water stops, filtration coefficient, hydrostatic and hydrodynamic headadditional measures to strengthen the rock mass, application of surfactants, injection of strengthening solutions into boreholes, drainage grooves
VIexcavation and transport of rock mass in the course of workunfavourable working conditions, rock mass stickingfiltration coefficient, degree of watering of rocksconstruction of drainage grooves, additional time spent on cleaning delivery vehicles
VIIremoval
fasteners
destruction of the fastening material, inability to remove the structure (due to destruction)chemical aggressiveness of water, location and characteristics of karsts and quicksandaccounting for the time spent on removing deformed support elements, additional measures to minimise water breakthroughs into the workings (injecting the massif with solutions, resins, etc.)
Table 6. Network formalisation of the problem of construction of workings in difficult hydrogeological conditions.
Table 6. Network formalisation of the problem of construction of workings in difficult hydrogeological conditions.
The Stage Designation in Figure 9Stage NameStarting PointFinal Peak (Peaks at Intermediate Stages)Interpretation
Iloading of fasteners on the surface13–51—the peak to start from; 3–5 options for transport technology
1–2, 2–52—availability of cleaning equipment, 5—transport technology
IIdelivery of fasteners in products, taking into account warehouse costs37–103—optimal technology after the first stage; 6—measures for sealing materials; 7–10 transport technologies
IIIpreparatory work (drilling holes, etc.)812–148—optimal technology after two stages; 11—measures for applying surfactants; 12–14 drilling technologies
IVconstruction of temporary support1216–1912—optimal technology after three stages; 15—additional costs for processing the massif; 16–19 options for the technology of erecting temporary support
Vconstruction of permanent support1621–2216—optimal technology after four stages; 20—additional costs for the construction of grooves;
21—22 fastening technologies
VIexcavation and transport of rock mass in the course of work2125–2721—optimal technology after five stages; 23—additional measures for cleaning vehicles,
24—measures for transporting rock mass; 25–27 transport technologies
VIIfastener extraction2529–3225—optimal technology after six stages; 28—additional measures for soil blasting; 29–32 dismantling technologies
VIIICompletion293429—optimal technology after seven stages; 33—additional measures to reinforce disturbed areas; 34—completion
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lousada, S.; Delehan, S.; Khorolskyi, A. Application of Dynamic Programming Models for Improvement of Technological Approaches to Combat Negative Water Leakage in the Underground Space. Water 2024, 16, 1952. https://doi.org/10.3390/w16141952

AMA Style

Lousada S, Delehan S, Khorolskyi A. Application of Dynamic Programming Models for Improvement of Technological Approaches to Combat Negative Water Leakage in the Underground Space. Water. 2024; 16(14):1952. https://doi.org/10.3390/w16141952

Chicago/Turabian Style

Lousada, Sérgio, Svitlana Delehan, and Andrii Khorolskyi. 2024. "Application of Dynamic Programming Models for Improvement of Technological Approaches to Combat Negative Water Leakage in the Underground Space" Water 16, no. 14: 1952. https://doi.org/10.3390/w16141952

APA Style

Lousada, S., Delehan, S., & Khorolskyi, A. (2024). Application of Dynamic Programming Models for Improvement of Technological Approaches to Combat Negative Water Leakage in the Underground Space. Water, 16(14), 1952. https://doi.org/10.3390/w16141952

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop