Method and system for calculating separated layer water accumulation amount of glutenite aquifer

Abstract

The invention provides a glutenite aquifer separation layer water accumulation amount calculation method and system. The method comprises the steps that the separation layer position is determined based on parameter information of a coal seam and a rock stratum above a goaf; for an area from a coal seam to the ground, simulating a rock stratum deformation and damage condition by using 3DEC software to obtain a target three-dimensional stratum model; obtaining upper and lower boundary curved surfaces in the separation position based on the separation water accumulation condition and the target three-dimensional stratum model; and calculating the separation water accumulation based on the upper and lower boundary curved surfaces. By means of the method, the glutenite separation layer water accumulation amount can be calculated more accurately.

Description

Technical Field

[0001] The present invention relates to the technical field of mine water disaster prevention and control, and particularly relates to a method and a system for calculating the accumulated water volume of the separated layer in a sandy conglomerate aquifer. Background Art

[0002] Mine water inrush is one of the three major dynamic disasters in coal mines, seriously affecting the safe production of coal mines. In the southern part of the Ordos Basin in China, Jurassic coal seams are mainly mined, which are threatened by the sandy conglomerate aquifer of the Cretaceous Zhidan Group in the roof. The coal seams in this area are characterized by large burial depth and high mining thickness, and the mining disturbance is large; the main water-inrush aquifers are the sandy conglomerate aquifers of the Cretaceous Luohe Formation and Yijun Formation. After mining, the gravel strata are broken, and the pore fissures are increased, greatly improving the permeability coefficient and water-richness of the aquifer, and increasing the threat of water disasters in the working face; a "cavity"-type separated layer is easily formed between the thick sandy conglomerate and the underlying soft argillaceous rock. With the development of fissures in the sandy conglomerate layer, it is connected to the "cavity" below, improving the water storage capacity of the separated layer, redistributing the groundwater, and forming a locally water-rich area. In addition to its own hydrostatic load, when the separated layer is filled with water, the separated layer water will also transfer a part of the load of the overlying strata. Under the action of the overall load, it may cause the instability and fracture of the lower water-resisting protective layer, resulting in water inrush. The special properties of the sandy conglomerate aquifer make the occurrence states of the separated layers developed in it different, and there are also different water-inrush modes. At the same time, the instantaneous water volume of the separated layer water inrush is large and the duration is short, causing great pressure on the working face drainage system. Therefore, the calculation of the accumulated water volume of the separated layer in the sandy conglomerate aquifer is of great significance for the safe production of mines. At present, there is little research on the development characteristics of the separated layer in the sandy conglomerate aquifer, especially the lack of a calculation method for the accumulated water volume of the sandy conglomerate separated layer. Summary of the Invention

[0003] The present invention aims to solve at least one of the technical problems in the related technologies to some extent.

[0004] To this end, the first object of the present invention is to propose a method for calculating the accumulated water volume of the separated layer in a sandy conglomerate aquifer to calculate the accumulated water volume of the sandy conglomerate separated layer more accurately.

[0005] The second object of the present invention is to propose a system for calculating the accumulated water volume of the separated layer in a sandy conglomerate aquifer.

[0006] The third object of the present invention is to propose an electronic device.

[0007] The fourth object of the present invention is to propose a computer-readable storage medium.

[0008] To achieve the above object, the first aspect of the present invention proposes a method for calculating the accumulated water volume of the separated layer in a sandy conglomerate aquifer, including:

[0009] Determining the position of the separated layer based on the parameter information of the strata above the coal seam and the goaf;

[0010] For the area between the coal seam and the ground surface, use 3DEC software to simulate the deformation and failure of rock strata to obtain a target three-dimensional strata model;

[0011] Based on the separated seam water accumulation situation and the target three-dimensional strata model, obtain the upper and lower boundary surfaces in the separated seam position;

[0012] Calculate the separated seam water accumulation volume based on the upper and lower boundary surfaces.

[0013] In the method of the first aspect of the present invention, the determination of the separated seam position based on the parameter information of the rock strata above the coal seam and the goaf includes: based on the parameter information of the rock strata above the coal seam, use the key stratum discriminant formula to determine whether there is a key stratum within the range of the sandy conglomerate aquifer; based on the parameter information of the rock strata in the water-conducting fracture zone, determine whether the goaf meets the spatial conditions for the development of separated seams; if the spatial conditions for the development of separated seams are met and there is a key stratum within the range of the sandy conglomerate aquifer, then determine the separated seam position in combination with the softening coefficient and water-resisting ability of the rock strata.

[0014] In the method of the first aspect of the present invention, the obtaining of the upper and lower boundary surfaces in the separated seam position based on the separated seam water accumulation situation and the target three-dimensional strata model includes: fitting the upper and lower boundary layers in the separated seam position in the target three-dimensional strata model to obtain a theoretical upper boundary surface and a theoretical lower boundary surface; for the separated seam water accumulation situation, calculate a corrected upper boundary surface and a corrected lower boundary surface; if the separated seam is full of water, the theoretical upper boundary surface is the target upper boundary surface, and the corrected lower boundary surface is the target lower boundary surface; if the separated seam is not full of water, the corrected upper boundary surface is the target upper boundary surface, and the theoretical lower boundary surface is the target lower boundary surface.

[0015] In the method of the first aspect of the present invention, the calculation of the separated seam water accumulation volume based on the upper and lower boundary surfaces includes: calculating the total volume between the upper and lower boundary layers in the separated seam position based on the target upper boundary surface and the target lower boundary surface; calculating the total volume of gravel blocks between the upper and lower boundary layers based on the target three-dimensional strata model; subtracting the total volume of gravel blocks from the total volume between the layers to obtain the separated seam water accumulation volume.

[0016] To achieve the above object, the second aspect of the present invention proposes a calculation system for the separated seam water accumulation volume of a sandy conglomerate aquifer, including:

[0017] A determination module for determining the separated seam position based on the parameter information of the rock strata above the coal seam and the goaf;

[0018] A modeling module for simulating the deformation and failure of rock strata for the area between the coal seam and the ground surface using 3DEC software to obtain a target three-dimensional strata model;

[0019] A fitting module, configured to obtain upper and lower boundary surfaces in the separated layer position based on the separated layer water accumulation situation and the target three-dimensional stratum model;

[0020] A calculation module, configured to calculate the separated layer water accumulation volume based on the upper and lower boundary surfaces.

[0021] In the system according to the second aspect of the present invention, the determination module is specifically configured to: based on the parameter information of the rock strata above the coal seam, use the key stratum discriminant formula to determine whether there is a key stratum within the range of the sandy conglomerate aquifer; based on the parameter information of the rock strata in the water-conducting fissure zone, determine whether the goaf meets the spatial conditions for the development of separated layers; if the spatial conditions for the development of separated layers are met and there is a key stratum within the range of the sandy conglomerate aquifer, then determine the separated layer position by combining the softening coefficient and water-resisting ability of the rock strata.

[0022] In the system according to the second aspect of the present invention, the fitting module is specifically configured to: fit the upper and lower boundary layers in the separated layer position in the target three-dimensional stratum model to obtain a theoretical upper boundary surface and a theoretical lower boundary surface; for the separated layer water accumulation situation, calculate a corrected upper boundary surface and a corrected lower boundary surface; if the separated layer is full of water, the theoretical upper boundary surface is the target upper boundary surface, and the corrected lower boundary surface is the target lower boundary surface; if the separated layer is not full of water, the corrected upper boundary surface is the target upper boundary surface, and the theoretical lower boundary surface is the target lower boundary surface.

[0023] In the system according to the second aspect of the present invention, the calculation module is specifically configured to: calculate the total volume between the upper and lower boundary layers in the separated layer position based on the target upper boundary surface and the target lower boundary surface; calculate the total volume of gravel blocks between the upper and lower boundary layers based on the target three-dimensional stratum model; subtract the total volume of gravel blocks from the total volume between the layers to obtain the separated layer water accumulation volume.

[0024] To achieve the above object, a third aspect of the present invention provides an electronic device, including: a processor, and a memory communicatively connected to the processor; the memory stores computer-executable instructions; the processor executes the computer-executable instructions stored in the memory to implement the method provided in the first aspect of the present invention.

[0025] To achieve the above object, a fourth aspect of the present invention provides a computer-readable storage medium, in which computer-executable instructions are stored, and when the computer-executable instructions are executed by a processor, they are used to implement the method provided in the first aspect of the present invention.

[0026] The method, system, electronic device and storage medium for calculating the separated layer water accumulation amount in a glutenite aquifer provided by the present invention determine the separated layer position based on the parameter information of the rock layers above the coal seam and the goaf; for the area between the coal seam and the ground surface, the 3DEC software is used to simulate the deformation and failure of the rock layers to obtain the target three-dimensional stratum model; based on the separated layer water accumulation situation and the target three-dimensional stratum model, the upper and lower boundary surfaces in the separated layer position are obtained; and the separated layer water accumulation amount is calculated based on the upper and lower boundary surfaces. In this case, compared with the water accumulation amount calculated without considering water accumulation in the prior art, the present invention comprehensively obtains the required upper and lower boundary surfaces based on the separated layer water accumulation situation and the target three-dimensional stratum model, and then calculates the separated layer water accumulation amount, which can calculate the separated layer water accumulation amount in the glutenite more accurately.

[0027] Additional aspects and advantages of the present invention will be given in part in the following description, become apparent in part from the following description, or be learned through the practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] The above and / or additional aspects and advantages of the present invention will become apparent and be readily understood from the following description of the embodiments in conjunction with the accompanying drawings, in which:

[0029] Figure 1 is a schematic flow chart of a method for calculating the separated layer water accumulation amount in a glutenite aquifer provided by an embodiment of the present invention;

[0030] Figure 2 is a schematic diagram of a rock layer provided by an embodiment of the present invention;

[0031] Figure 3 is a graph of the rock layer settlement surface fitted by Python provided by an embodiment of the present invention;

[0032] Figure 4 is a schematic diagram of a method for calculating the separated layer volume using Python provided by an embodiment of the present invention;

[0033] Figure 5 is a block diagram of a system for calculating the separated layer water accumulation amount in a glutenite aquifer provided by an embodiment of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0034] The embodiments of the present invention will be described in detail below. The examples of the embodiments are shown in the accompanying drawings, in which the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below by referring to the accompanying drawings are exemplary and are intended to explain the present invention, and should not be construed as a limitation of the present invention.

[0035] The method and system for calculating the separated layer water accumulation amount in a glutenite aquifer according to the embodiments of the present invention will be described below with reference to the accompanying drawings.

[0036] The drainage capacity of separated seam water in a mine needs to be determined according to the predicted volume of separated seam water. Existing mines threatened by separated seam water in sandy conglomerate mainly refer to the predicted water inrush volume of separated seam water in this mine or adjacent mines. In the past, an empirical formula for the volume of separated seam was derived based on the thin plate small deflection theory, which did not consider the influence of accumulated water on the deflection of the rock stratum, and the calculation results often deviated greatly from the actual situation.

[0037] Based on this, the embodiment of the present invention provides a method for calculating the accumulated water volume of separated seam in sandy conglomerate aquifer to calculate the accumulated water volume of separated seam in sandy conglomerate more accurately.

[0038] Figure 1 It is a schematic flow chart of a method for calculating the accumulated water volume of separated seam in sandy conglomerate aquifer provided by the embodiment of the present invention.

[0039] As Figure 1 shown, the method for calculating the accumulated water volume of separated seam in sandy conglomerate aquifer includes the following steps:

[0040] Step S101, determining the separated seam position based on the parameter information of the rock stratum above the coal seam and the goaf.

[0041] In step S101, determining the separated seam position based on the parameter information of the rock stratum above the coal seam and the goaf includes: based on the parameter information of the rock stratum above the coal seam, using the key stratum discriminant formula to determine whether there is a key stratum within the range of the sandy conglomerate aquifer; based on the parameter information of the rock stratum in the water-conducting fissure zone, determining whether the goaf meets the spatial conditions for the development of separated seam; if the spatial conditions for the development of separated seam are met and there is a key stratum within the range of the sandy conglomerate aquifer, then determining the separated seam position in combination with the softening coefficient and water-resisting ability of the rock stratum.

[0042] In step S101, the parameter information of the rock stratum above the coal seam includes the elastic modulus, thickness, and specific weight of each rock stratum between the coal seam and the ground surface, etc. The parameter information of the rock stratum in the water-conducting fissure zone includes the original thickness and swelling coefficient of each rock stratum, etc.

[0043] In step S101, the parameter information of the rock stratum above the coal seam and the goaf also includes the softening coefficient of the rock stratum, the type of rock stratum, the water-richness of the aquifer, the water level of the aquifer, etc.

[0044] Specifically, Figure 2 It is a schematic diagram of the rock stratum provided by the embodiment of the present invention. As Figure 2 shown, there are multiple layers of rock strata between the underground coal seam and the ground surface. After the coal seam is mined along the working face advancing direction, the mined part forms a goaf, and the roof rock stratum is affected by the mining pressure and gradually forms a caving zone. Above the caving zone, a fissure zone and a bending subsidence zone are formed accordingly. Among them, the caving zone and the fissure zone form the water-conducting fissure zone.

[0045] Determine the separation position of the sandy gravel aquifer that threatens water inrush to the working face based on the parameter information of the strata above the coal seam and the goaf. The determination of the separation position is mainly determined from two conditions: The first condition is to determine whether there is a stratum structure with developed separation within the range of the sandy gravel aquifer between the coal seam and the ground, that is, to determine whether there are key strata within the range of the sandy gravel aquifer; The second condition is to determine whether the goaf meets the spatial conditions for the development of separation. If both conditions are met, it indicates that separation occurs within the range of the sandy gravel aquifer.

[0046] For the first condition, first identify all the key strata between the coal seam and the ground through the key stratum discriminant formula, where the key stratum discriminant formula satisfies:

[0047]

[0048] In the formula, E n+1 is the elastic modulus of the (n + 1)-th layer of rock above the coal seam, h n+1 is the thickness of the (n + 1)-th layer of rock above the coal seam, γ i is the specific weight of the i-th layer of rock above the coal seam, h i is the thickness of the i-th layer of rock above the coal seam, represents the cumulative sum of the products of the specific weights and thicknesses of the first to n-th layers of rock above the coal seam. γ n+1 is the specific weight of the (n + 1)-th layer of rock above the coal seam, E i is the elastic modulus of the i-th layer of rock above the coal seam, h i is the thickness of the i-th layer of rock above the coal seam. Starting from the first layer of rock above the coal seam, when the formula (1) is first satisfied, the (n + 1)-th layer can be determined as a key stratum, which is the first key stratum above the coal seam. Starting from the first key stratum, re-use the formula (2) to determine the second key stratum, and so on to determine all the key strata between the coal seam and the ground.

[0049] After determining all the key strata between the coal seam and the ground, based on the range of the sandy gravel aquifer, determine whether there are key strata within the range of the sandy gravel aquifer. If so, screen out the key strata contained within the range of the sandy gravel aquifer.

[0050] For the second condition, use the development discriminant formula to determine whether the goaf meets the spatial conditions for the development of separation, where the development discriminant formula satisfies:

[0051] M - ∑(K i - 1)d i > 0 (2)

[0052] In the formula, M is the total thickness of the mined coal seam, K i is the swelling coefficient of the i-th layer of rock in the water-conducting fractured zone after coal mining, d iis the original thickness of the i-th layer of rock in the water-conducting fissure zone after coal seam mining. ∑(K i -1)d i The summation layer number targeted is all the rock layer numbers in the water-conducting fissure zone after coal seam mining. If these two conditions of Equation (1) and Equation (2) are simultaneously satisfied, it indicates that bedding separation occurs within the range of the glutenite aquifer.

[0053] After determining that bedding separation occurs within the range of the glutenite aquifer, analyze the range of bedding separation through the rock mechanical characteristics of the rock layers. Considering that the upper and lower boundary layers of bedding separation should be hard rock and soft rock respectively, and the soft rock needs to have a certain water-resisting ability, and the glutenite aquifer is mainly composed of an interbedded lithological combination of sandstone and conglomerate, so the development of bedding separation in the glutenite aquifer is usually a regional bedding separation spanning multiple rock layers. Determine the target key layer according to the key layer within the range of the glutenite aquifer. The target key layer is the key layer closest to the ground among the key layers within the range of the glutenite aquifer. Analyze layer by layer from the target key layer downwards until a rock layer whose softening coefficient and water-resisting ability meet the requirements is found, and use this rock layer whose softening coefficient and water-resisting ability meet the requirements as the lower boundary layer of bedding separation. The softening coefficient and water-resisting ability meeting the requirements mean that the softening coefficient is less than or equal to the set coefficient threshold (such as 0.75 MPa), and it is relatively water-resistant (that is, this rock layer is a mudstone-like rock layer). Use the target key layer as the upper boundary layer of bedding separation.

[0054] Determine the water-richness of the aquifer according to the unit water inflow q of the borehole. Based on the water-richness of the aquifer, the water level of the aquifer and the position of bedding separation, judge the water accumulation situation of the bedding separation in the glutenite aquifer. The water accumulation situation is divided into two situations: full of water and not full of water. Among them, when the water-richness of the aquifer is strong or extremely strong (for example, the unit water inflow of the borehole satisfies 1.0 L / (s·m) < q ≤ 5.0 L / (s·m) for strong water-richness; satisfies q > 5.0 L / (s·m) for extremely strong water-richness), the water level of the aquifer is higher than the determined position of bedding separation (here it means the water level of the aquifer is higher than the upper boundary layer of bedding separation) and the water level shows a continuous and slow decline, it can be judged that the bedding separation is full of water, otherwise the bedding separation is not full of water.

[0055] Step S102, for the area between the coal seam and the ground surface, use 3DEC software to simulate the deformation and failure of the rock layers to obtain the target three-dimensional strata model.

[0056] Specifically, in step S102, for the area between the coal seam and the ground surface, an initial three-dimensional strata model of the overlying rock in this area is established using 3DEC (3Dimension Distinct Element Code) software. Vertically, it is established based on the actual strata data revealed by boreholes. Randomly generated Thiessen polygons are used to simulate gravel blocks of random shapes and sizes, and the side length range is set to limit the size of the gravel blocks. Horizontally, the dip width of the initial three-dimensional strata model should be consistent with the actual engineering background, and the surrounding rock with a set length (e.g., 30 m) or more should be reserved at the front and back. To accurately simulate the law of separation generation to closure, the strike length of the initial three-dimensional strata model should not be less than twice the breakage distance of the key stratum within the range of the sandy conglomerate aquifer (if there are multiple key strata, the maximum value of the breakage moment is taken). The calculation formula for the breakage distance of the key stratum is as follows:

[0057]

[0058] In the formula, l is the breakage distance of the key stratum, h is the thickness of the key stratum being calculated, R t is the tensile strength of the key stratum being calculated, and q is the load of the overlying strata of the key stratum being calculated.

[0059] Considering that the initial three-dimensional strata model is a simulation model of the rock strata above the unmined coal seam, to better simulate the deformation and failure of the overlying rock after mining, the Mohr-Coulomb model is used to describe the mechanical behavior of rocks, and the Coulomb slip criterion is used to describe the mechanical behavior of rock joints to perform yield fracture and caving. Thus, the deformation and failure of the overlying rock are simulated through simulated excavation to obtain the target three-dimensional strata model.

[0060] Step S103, obtaining the upper and lower boundary surfaces in the separation position based on the separation water accumulation situation and the target three-dimensional strata model.

[0061] In step S103, obtaining the upper and lower boundary surfaces in the separation position based on the separation water accumulation situation and the target three-dimensional strata model includes: fitting the upper and lower boundary layers in the separation position of the target three-dimensional strata model to obtain the theoretical upper boundary surface and the theoretical lower boundary surface; calculating the corrected upper boundary surface and the corrected lower boundary surface for the separation water accumulation situation; if the separation is full of water, the theoretical upper boundary surface is the target upper boundary surface, and the corrected lower boundary surface is the target lower boundary surface; if the separation is not full of water, the corrected upper boundary surface is the target upper boundary surface, and the theoretical lower boundary surface is the target lower boundary surface.

[0062] Specifically, the process of obtaining the theoretical upper boundary surface and the theoretical lower boundary surface includes:

[0063] 11) In the initial equilibrium state of the 3DEC model (i.e., for the initial three-dimensional stratum model), traverse all the gravel blocks in the model. According to the positions of the separated layer's upper and lower boundary layers determined in step S101, use conditional statements to extract the IDs (identification numbers) of the gravel blocks in the upper and lower boundary layers, and store them in tables list1 and list2 in Python in the form of a list. Extract the IDs of the gravel blocks in all the rock layers within the separated layer (i.e., between the upper and lower boundary layers), and store them in table list3.

[0064] 12) Execute the excavation command to simulate the deformation and failure of the overlying rock to obtain the target three-dimensional stratum model, and update list1, table list2, and table list3 based on the IDs of the gravel blocks in the target three-dimensional stratum model. Traverse table list1 and table list2 respectively, obtain the position coordinates of each gravel block after the excavation causes the rock layer to move, and fit the settlement surface w 1 (x, y) and the settlement surface w 2 (x, y). Among them, the settlement surface w 1 (x, y) is the theoretical lower boundary surface, and the settlement surface w 2 (x, y) is the theoretical upper boundary surface. Thus, using Python for numerical programming in 3DEC software can achieve surface fitting for the upper and lower boundary layers of the separated layer determined in step S101. Figure 3 This is the rock layer settlement surface diagram fitted by Python provided by the embodiment of the present invention. Figure 3 It shows a rock layer settlement surface diagram of the entire sandy conglomerate aquifer fitted based on the target three-dimensional stratum model. w is the settlement surface of any layer within the sandy conglomerate aquifer.

[0065] The process of obtaining the corrected upper boundary surface and the corrected lower boundary surface includes:

[0066] 21) Since the influence of separated layer water accumulation on the settlement of the lower boundary of the separated layer is not considered in the target three-dimensional stratum model, consider the influence of water accumulation on the settlement of the lower boundary of the separated layer. For the upper and lower boundary layers of the separated layer above the goaf, obtain the rock layer deflection surface by disassembling the rock layer deflection micro-surface equation under separated layer water accumulation, and correct the lower boundary of the separated layer. Therefore, the corrected lower boundary surface satisfies:

[0067]

[0068] In the formula, where q is the lateral load on the lower boundary layer of the separated layer above the goaf, the series expanded by m, a and b are the lengths along the working face strike and the dip widths of the lower boundary layer of the separated layer above the goaf, D is the flexural rigidity, ρ is the density of groundwater, g is the acceleration due to gravity. (x, y) are the coordinates of the gravel block.

[0069] 22) During the face mining process, based on the mining progress and the water level change of the sandy conglomerate aquifer, dynamically analyze the current separated layer water accumulation situation. If the separated layer is not filled with water, the water level of the sandy conglomerate aquifer needs to be used as the upper boundary, that is, the corrected upper boundary surface is obtained by using the water level of the sandy conglomerate aquifer.

[0070] Determine the final required target upper boundary surface and target lower boundary surface considering the separated layer water accumulation situation.

[0071] Specifically, if the separated layer is filled with water: the required upper boundary is the upper boundary determined by the separated layer discrimination in step S101, so the target upper boundary surface is the theoretically fitted upper boundary surface in Python; the difference in the lower boundary is the rock layer deflection surface obtained by dissociating the rock layer deflection micro-surface equation under separated layer water accumulation, that is, the target lower boundary surface is the corrected lower boundary surface.

[0072] If the separated layer is not filled with water: the upper boundary is the water level of the aquifer, that is, the target upper boundary surface is the corrected upper boundary surface; the lower boundary is the lower surface fitted in Python, that is, the target lower boundary surface is the theoretically lower boundary surface.

[0073] Step S104, calculate the separated layer water accumulation volume based on the upper and lower boundary surfaces.

[0074] In step S104, calculating the separated layer water accumulation volume based on the upper and lower boundary surfaces includes: calculating the total volume between the upper and lower boundary layers in the separated layer position based on the target upper boundary surface and the target lower boundary surface; calculating the total volume of gravel blocks between the upper and lower boundary layers based on the target three-dimensional strata model; subtracting the total volume of gravel blocks from the total volume between the layers to obtain the separated layer water accumulation volume.

[0075] The separated layer water accumulation volume satisfies:

[0076] V f =∫∫w 1,目标 -w 2,目标 dxdy-V (5)

[0077] In the formula, V f is the separated layer water accumulation volume, w 1,目标 is the target upper boundary surface, w 2,目标 is the target lower boundary surface. V is the total volume of gravel blocks, which can be obtained by traversing list3 to get the volume of each gravel block. The integration range is the range of the mined coal seam.

[0078] Figure 4 This is a schematic diagram of the method for calculating the separated layer volume using Python provided by the embodiment of the present invention. As Figure 4As shown in the figure, a bending subsidence zone is formed on the fissure zone. The upper and lower boundary layers of the separated zone position are determined by using step S101. The water level of the aquifer above the upper boundary of the separated zone affects the calculation of the subsequent accumulated water volume in the separated zone. For the upper boundary and the lower boundary of the separated zone, the fitted subsidence surfaces are obtained by using Python respectively. Combining the water accumulation situation in the separated zone, the target upper boundary surface and the target lower boundary surface are determined. The total volume of gravel blocks between the upper and lower boundary layers is read, and the accumulated water volume in the separated zone is obtained by combining the total volume between the layers and the total volume of gravel blocks.

[0079] To implement the above embodiments, the present invention also proposes a system for calculating the accumulated water volume in the separated zone of a glutenite aquifer.

[0080] Figure 5 It is a block diagram of a system for calculating the accumulated water volume in the separated zone of a glutenite aquifer provided by an embodiment of the present invention.

[0081] As Figure 5 shown, the system for calculating the accumulated water volume in the separated zone of the glutenite aquifer includes a determination module 11, a modeling module 12, a fitting module 13 and a calculation module 14, wherein:

[0082] The determination module 11 is used to determine the separated zone position based on the parameter information of the rock strata above the coal seam and the goaf.

[0083] The modeling module 12 is used to simulate the deformation and failure of the rock strata by using 3DEC software for the area between the coal seam and the ground surface to obtain a target three-dimensional strata model.

[0084] The fitting module 13 is used to obtain the upper and lower boundary surfaces in the separated zone position based on the water accumulation situation in the separated zone and the target three-dimensional strata model.

[0085] The calculation module 14 is used to calculate the accumulated water volume in the separated zone based on the upper and lower boundary surfaces.

[0086] Further, in a possible implementation manner of the embodiment of the present invention, the determination module 11 is specifically used for: based on the parameter information of the rock strata above the coal seam, using the key stratum discriminant formula to determine whether there is a key stratum within the range of the glutenite aquifer; based on the parameter information of the rock strata in the water-conducting fissure zone, determining whether the goaf meets the spatial conditions for the development of a separated zone; if the spatial conditions for the development of a separated zone are met and there is a key stratum within the range of the glutenite aquifer, then the separated zone position is determined by combining the softening coefficient and the water-resisting ability of the rock strata.

[0087] Further, in a possible implementation manner of the embodiment of the present invention, the fitting module 13 is specifically configured to: fit the upper and lower boundary layers at the separation layer position in the target three-dimensional strata model to obtain a theoretical upper boundary surface and a theoretical lower boundary surface; calculate a corrected upper boundary surface and a corrected lower boundary surface for the case of water accumulation in the separation layer; if the separation layer is full of water, the theoretical upper boundary surface is the target upper boundary surface, and the corrected lower boundary surface is the target lower boundary surface; if the separation layer is not full of water, the corrected upper boundary surface is the target upper boundary surface, and the theoretical lower boundary surface is the target lower boundary surface.

[0088] Further, in a possible implementation manner of the embodiment of the present invention, the calculation module 14 is specifically configured to: calculate the total volume between the upper and lower boundary layers in the separation layer position based on the target upper boundary surface and the target lower boundary surface; calculate the total volume of gravel blocks between the upper and lower boundary layers based on the target three-dimensional strata model; subtract the total volume of gravel blocks from the total volume between the layers to obtain the water accumulation volume in the separation layer.

[0089] It should be noted that the foregoing explanation of the embodiment of the calculation method for the water accumulation volume in the separation layer of the glutenite aquifer also applies to the calculation system for the water accumulation volume in the separation layer of the glutenite aquifer in this embodiment, and will not be elaborated here.

[0090] In the embodiment of the present invention, the separation layer position is determined based on the parameter information of the coal seam and the rock strata above the goaf; for the area between the coal seam and the ground surface, the 3DEC software is used to simulate the deformation and failure of the rock strata to obtain the target three-dimensional strata model; the upper and lower boundary surfaces in the separation layer position are obtained based on the water accumulation situation in the separation layer and the target three-dimensional strata model; the water accumulation volume in the separation layer is calculated based on the upper and lower boundary surfaces. In this case, compared with the water accumulation volume calculated without considering water accumulation in the prior art, the present invention comprehensively considers the water accumulation situation in the separation layer and the target three-dimensional strata model to obtain the required upper and lower boundary surfaces, and then calculates the water accumulation volume in the separation layer, which can more accurately calculate the water accumulation volume in the glutenite separation layer.

[0091] Compared with the prior art, the calculation method for the water accumulation volume in the separation layer of the glutenite aquifer of the present invention can better target the geological structure and rock mechanics characteristics, and has strong applicability. It has the following advantages:

[0092] The construction of the drainage system in mines threatened by separated layer water needs to consider the sudden inflow water volume of the separated layer. Under the influence of mining, the pores and fissures of the glutenite aquifer are developed and connected to the possible "cavity" formed by the separated layer, forming a good water storage space. At this time, the traditional separated layer volume calculation formula cannot accurately evaluate the volume of the separated layer water. The separated layer volume calculation method based on 3DEC-Python of the present invention can accurately calculate the fissure volume of the conglomerate.

[0093] The water accumulation in the separated strata is a dynamic process. Considering the dynamic development of the separated strata, the present invention calculates the water accumulation amount in the separated strata before the working face is mined, providing a basis for the water control design of the working face. During the mining process, according to the mining progress and the change of the water level in the aquifer, the water accumulation amount is dynamically corrected, facilitating the mine to make real-time adjustments to the drainage system of the working face. The above method for calculating the water accumulation amount in the separated strata is reliable and can provide guarantee for the safe mining of the mine.

[0094] To implement the above embodiments, the present invention also provides an electronic device, including: a processor, and a memory communicatively connected to the processor; the memory stores computer-executable instructions; the processor executes the computer-executable instructions stored in the memory to implement the method provided in the foregoing embodiments.

[0095] To implement the above embodiments, the present invention also provides a computer-readable storage medium storing computer-executable instructions, which are used to implement the method provided in the foregoing embodiments when executed by a processor.

[0096] To implement the above embodiments, the present invention also provides a computer program product including a computer program, which implements the method provided in the foregoing embodiments when executed by a processor.

[0097] In the description of the foregoing embodiments, the descriptions referring to terms such as "one embodiment", "some embodiments", "example", "specific example", or "some examples" mean that the specific features, structures, materials, or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms are not necessarily directed to the same embodiment or example. Moreover, the specific features, structures, materials, or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, without contradiction, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of different embodiments or examples.

[0098] In addition, the terms "first" and "second" are only used for descriptive purposes and cannot be construed as indicating or implying relative importance or implicitly indicating the quantity of the indicated technical features. Thus, the features defined with "first" and "second" may explicitly or implicitly include at least one of the features. In the description of the present invention, "a plurality" means at least two, such as two, three, etc., unless otherwise specifically defined.

[0099] Any process or method description represented in a flowchart or otherwise described herein can be understood to represent a module, segment, or portion of code including one or more executable instructions for implementing a customized logic function or process. The scope of the preferred embodiments of the present invention includes additional implementations, where functions may be executed in a substantially simultaneous manner or in a reverse order according to the functions involved, rather than in the order shown or discussed, which should be understood by those skilled in the art to which the embodiments of the present invention pertain.

[0100] The logic and / or steps represented in a flowchart or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing a logical function, and can be embodied specifically in any computer-readable medium for use by or in connection with an instruction execution system, apparatus, or device, such as a computer-based system, a system including a processor, or other systems that can fetch and execute instructions from the instruction execution system, apparatus, or device. For the purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of the computer-readable medium include the following: an electrical connection having one or more wires (electronic device), a portable computer diskette (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). Additionally, the computer-readable medium can even be paper or other suitable medium on which the program can be printed, as the program can be obtained electronically, for example, by optically scanning the paper or other medium, followed by editing, interpretation, or otherwise processing as appropriate, and then stored in a computer memory.

[0101] It should be understood that various parts of the present invention can be implemented in hardware, software, firmware, or a combination thereof. In the above embodiments, multiple steps or methods can be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, any one or a combination of the following techniques well known in the art can be used: discrete logic circuits having logic gate circuits for implementing logical functions on data signals, application specific integrated circuits having appropriate combinational logic gate circuits, programmable gate arrays (PGAs), field programmable gate arrays (FPGAs), etc.

[0102] Those of ordinary skill in the art can understand that all or part of the steps carried out in implementing the method of the above embodiments can be completed by instructing relevant hardware through a program, and the program can be stored in a computer-readable storage medium. When the program is executed, it includes one or a combination of the steps of the method embodiments.

[0103] In addition, each functional unit in various embodiments of the present invention can be integrated into a processing module, or each unit can exist physically alone, or two or more units can be integrated into one module. The above-mentioned integrated module can be implemented in the form of hardware or in the form of a software functional module. When the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it can also be stored in a computer-readable storage medium.

[0104] The above-mentioned storage medium can be a read-only memory, a magnetic disk, an optical disc, etc. Although the embodiments of the present invention have been shown and described above, it can be understood that the above embodiments are exemplary and should not be construed as limiting the present invention. Those of ordinary skill in the art can make changes, modifications, substitutions, and variations to the above embodiments within the scope of the present invention.

Claims

1. A method for calculating the water accumulation in the stratum of a sandy conglomerate aquifer, characterized in that: include: Determine the separation position based on the parameter information of the coal seam and the rock layer above the goaf; For the area between the coal seam and the ground, 3DEC software is used to simulate the deformation and damage of the rock formation to obtain the target three-dimensional formation model; Based on the water accumulation in the stratum and the target three-dimensional stratum model, upper and lower boundary surfaces in the stratum separation position are obtained; The amount of water accumulated in the separation layer is calculated based on the upper and lower boundary surfaces.

2. The method for calculating the amount of water accumulation in the stratum of a sandy conglomerate aquifer according to claim 1, characterized in that: The step of determining the separation position based on the parameter information of the coal seam and the rock layer above the goaf includes: Based on the parameter information of the rock layer above the coal seam, the key layer discriminant formula is used to determine whether the gravel aquifer contains a key layer. Based on the parameter information of the rock formation in the water-conducting fracture zone, it is determined whether the goaf meets the spatial conditions for developing separation layers; If the spatial conditions for the development of abscission are met and the gravel aquifer contains a key layer, the location of the abscission layer is determined in combination with the softening coefficient and water-isolating capacity of the rock formation.

3. The method for calculating the amount of water accumulation in the stratum of a sandy conglomerate aquifer according to claim 1, characterized in that: The step of obtaining upper and lower boundary surfaces in the delamination position based on the delamination water accumulation condition and the target three-dimensional stratum model includes: Based on the upper and lower boundary layer fitting of the separation position in the target three-dimensional formation model, a theoretical boundary surface and a theoretical lower boundary surface are obtained; According to the water accumulation in the separation layer, the corrected upper boundary surface and the corrected lower boundary surface are calculated; If the detached layer is full of water, the theoretical boundary surface is the target upper boundary surface, and the modified lower boundary surface is the target lower boundary surface; if the detached layer is not full of water, the modified upper boundary surface is the target upper boundary surface, and the theoretical lower boundary surface is the target lower boundary surface.

4. The method for calculating the amount of water accumulation in the stratum of a sandy conglomerate aquifer according to claim 3, characterized in that: The calculating of the amount of water accumulation in the stratum based on the upper and lower boundary surfaces comprises: Calculating the total interlayer volume between the upper and lower boundary layers in the separation position based on the target upper boundary curved surface and the target lower boundary curved surface; Calculate the total volume of gravel blocks between the upper and lower boundary layers based on the target three-dimensional stratum model; The total volume of the interlayer minus the total volume of the gravel blocks gives the amount of water accumulated in the abscission layer.

5. A system for calculating water accumulation in a sandy conglomerate aquifer, characterized in that: include: A determination module, used for determining the separation position based on parameter information of the coal seam and the rock layer above the goaf; Modeling module, used to simulate the deformation and damage of rock formations in the area between coal seams and the ground using 3DEC software to obtain the target three-dimensional formation model; A fitting module, used for obtaining upper and lower boundary surfaces in the delamination position based on the delamination water accumulation condition and the target three-dimensional formation model; A calculation module is used to calculate the amount of water accumulation in the stratum based on the upper and lower boundary surfaces.

6. The system for calculating the amount of water accumulation in the stratum of sandy conglomerate aquifer according to claim 5, characterized in that: The determination module is specifically used for: based on the parameter information of the rock layer above the coal seam, using the key layer discriminant formula to determine whether the gravel aquifer contains a key layer; based on the parameter information of the rock layer in the water-conducting fracture zone, determining whether the goaf meets the spatial conditions for the development of strata; if the spatial conditions for the development of strata are met and the gravel aquifer contains a key layer, determining the stratum separation position in combination with the softening coefficient and water-isolating capacity of the rock layer.

7. The system for calculating the amount of water accumulation in the stratum of sandy conglomerate aquifer according to claim 5, characterized in that: The fitting module is specifically used to: obtain a theoretical boundary surface and a theoretical lower boundary surface based on the upper and lower boundary layer fitting of the abscission layer position in the target three-dimensional formation model; calculate a corrected upper boundary surface and a corrected lower boundary surface according to the water accumulation in the abscission layer; if the abscission layer is full of water, the theoretical boundary surface is the target upper boundary surface, and the corrected lower boundary surface is the target lower boundary surface; if the abscission layer is not full of water, the corrected upper boundary surface is the target upper boundary surface, and the theoretical lower boundary surface is the target lower boundary surface.

8. The system for calculating water accumulation in the stratum of sandy conglomerate aquifer according to claim 7, characterized in that: The calculation module is specifically used to: calculate the total interlayer volume between the upper and lower boundary layers in the delamination position based on the target upper boundary surface and the target lower boundary surface; calculate the total volume of gravel blocks between the upper and lower boundary layers based on the target three-dimensional formation model; and obtain the delamination water volume by subtracting the total volume of gravel blocks from the total interlayer volume.

9. An electronic device, characterized in that: include: A processor, and a memory communicatively connected to the processor; The memory stores computer-executable instructions; The processor executes the computer-executable instructions stored in the memory to implement the method according to any one of claims 1 to 4.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer-executable instructions, which are used to implement the method according to any one of claims 1 to 4 when executed by a processor.

Images