A method and system for calculating the amount of water in the interlayer of a sandy conglomerate aquifer

Abstract

The present application provides a kind of sandy conglomerate aquifer interlayer water content calculation method and system, wherein the method comprises determining the interlayer position based on the parameter information of the rock stratum above the coal seam and goaf;For the area between the coal seam and the ground, the 3DEC software is used to simulate the deformation and failure of rock stratum to obtain the target three-dimensional stratum model;Based on the interlayer water accumulation and the target three-dimensional stratum model, the upper and lower boundary surfaces in the interlayer position are obtained;The interlayer water content is calculated based on the upper and lower boundary surfaces.The method of the present application can more accurately calculate the sandy conglomerate interlayer water content.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of mine water disaster prevention, and particularly relates to a method and system for calculating water content in a sandstone aquifer bed separation. BACKGROUND

[0002] Mine water inrush is one of the three major dynamic disasters in coal mines, which seriously affects the safety production of coal mines. In the southern Ordos Basin of China, Jurassic coal seams are mainly mined, which are threatened by the sandstone aquifer bed of Cretaceous Zhidan Group in the roof. The coal seams in this area have the characteristics of large burial depth, large mining thickness, and large mining disturbance. The main water-filled aquifer is the sandstone aquifer bed of Cretaceous Luohe and Yijun groups. After mining, the conglomerate layer is broken, the pore and fracture are increased, the permeability coefficient and water abundance of the aquifer are greatly improved, and the water hazard threat of the working face is increased. The thick sandstone and gravel layer is easy to form a "cavity" type bed separation between the lower soft mudstone. With the development of sandstone and gravel layer fissure, it is connected with the lower "cavity", which improves the water storage capacity of the bed separation, redistributes the underground water, and forms a local water-rich area. In addition to the static water load of the bed separation itself, when the bed separation is filled with water, it will also transfer part of the load of the overlying rock. Under the overall load, the lower water-resisting protective layer may be unstable and broken, resulting in water inrush. The special properties of the sandstone aquifer bed make the bed separation developed in different ways, and the bed separation water inrush has large instantaneous water quantity and short duration, which causes large pressure on the working face drainage system. Therefore, the calculation of the water content in the sandstone aquifer bed separation is of great significance to the safety production of mines. At present, there is little research on the development characteristics of the sandstone aquifer bed separation, especially the calculation method of the water content in the sandstone bed separation. SUMMARY

[0003] The present application aims to at least partially solve one of the problems in the related art.

[0004] To this end, the first object of the present application is to provide a method for calculating the water content in the sandstone aquifer bed separation, so as to more accurately calculate the water content in the sandstone bed separation.

[0005] The second object of the present application is to provide a system for calculating the water content in the sandstone aquifer bed separation.

[0006] The third object of the present application is to provide an electronic device.

[0007] The fourth object of the present application is to provide a computer readable storage medium.

[0008] To achieve the above objects, the first aspect of the present application provides a method for calculating the water content in the sandstone aquifer bed separation, comprising:

[0009] determining the position of the bed separation based on the parameter information of the rock layer above the coal seam and the goaf;

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

[0011] Based on the separation water accumulation condition and the target three-dimensional stratum model, the upper and lower boundary curved surfaces in the separation position are obtained;

[0012] Based on the upper and lower boundary curved surfaces, the separation water accumulation amount is calculated.

[0013] In the method of the first aspect of the present application, the separation position is determined based on the parameter information of the rock stratum above the coal seam and the goaf, including: based on the parameter information of the rock stratum above the coal seam, using the key layer discriminant to determine whether the sand and gravel aquifer contains a key layer; based on the parameter information of the rock stratum of the water flowing fractured zone, it is determined whether the goaf meets the spatial conditions for developing separation; if the spatial conditions for developing separation are met and the sand and gravel aquifer contains a key layer, the separation position is determined in combination with the softening coefficient and the water resistance of the rock stratum.

[0014] In the method of the first aspect of the present application, the upper and lower boundary curved surfaces in the separation position are obtained based on the separation water accumulation condition and the target three-dimensional stratum model, including: based on the upper and lower boundary layers of the separation position in the target three-dimensional stratum model, the theoretical upper boundary curved surface and the theoretical lower boundary curved surface are fitted; for the separation water accumulation condition, the corrected upper boundary curved surface and the corrected lower boundary curved surface are calculated; if the separation is full of water, the theoretical upper boundary curved surface is the target upper boundary curved surface, and the corrected lower boundary curved surface is the target lower boundary curved surface; if the separation is not full of water, the corrected upper boundary curved surface is the target upper boundary curved surface, and the theoretical lower boundary curved surface is the target lower boundary curved surface.

[0015] In the method of the first aspect of the present application, the separation water accumulation amount is calculated based on the upper and lower boundary curved surfaces, including: based on the target upper boundary curved surface and the target lower boundary curved surface, the total volume between the upper and lower boundary layers in the separation position is calculated; based on the target three-dimensional stratum model, the total volume of the gravel blocks between the upper and lower boundary layers is calculated; the total volume between the layers is subtracted from the total volume of the gravel blocks to obtain the separation water accumulation amount.

[0016] To achieve the above purpose, the second aspect of the present application provides a sand and gravel aquifer separation water accumulation amount calculation system, which comprises:

[0017] The determination module is used to determine the separation position based on the parameter information of the rock stratum above the coal seam and the goaf;

[0018] The modeling module is used to simulate the deformation and failure of the rock stratum for the area between the coal seam and the ground to obtain a target three-dimensional stratum model;

[0019] a fitting module configured to fit an upper boundary surface and a lower boundary surface in the position of the separation layer based on the water-accumulation condition of the separation layer and the target three-dimensional stratum model;

[0020] a calculating module configured to calculate the water-accumulation volume of the separation layer based on the upper boundary surface and the lower boundary surface.

[0021] In the system of the second aspect of the present application, the determining module is specifically configured to: determine whether the sandstone aquifer contains a key layer based on the parameter information of the rock layer above the coal seam using a key layer discriminant; determine whether the goaf satisfies the spatial condition for developing the separation layer based on the parameter information of the rock layer of the water-conducting fractured zone; and determine the position of the separation layer in combination with the softening coefficient and the water-resisting capacity of the rock layer if the goaf satisfies the spatial condition for developing the separation layer and the sandstone aquifer contains the key layer.

[0022] In the system of the second aspect of the present application, the fitting module is specifically configured to: fit the theoretical upper boundary surface and the theoretical lower boundary surface based on the upper and lower boundary layers in the position of the separation layer in the target three-dimensional stratum model; calculate the corrected upper boundary surface and the corrected lower boundary surface for the water-accumulation condition of 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; and 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.

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

[0024] To achieve the above object, the third aspect of the present application provides an electronic device, which comprises a processor and a memory connected with the processor in communication; the memory stores computer-executed instructions; and the processor executes the computer-executed instructions stored in the memory to implement the method provided in the first aspect of the present application.

[0025] To achieve the above object, the fourth aspect of the present application provides a computer-readable storage medium, which stores computer-executed instructions; and the computer-executed instructions are executed by a processor to implement the method provided in the first aspect of the present application.

[0026] The sandstone aquifer parting layer water accumulation amount calculation method, system, electronic device and storage medium provided by the present application determine the parting layer 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, the 3DEC software is used to simulate the rock layer deformation and failure to obtain a target three-dimensional stratum model; the upper and lower boundary surfaces in the parting layer position are obtained based on the parting layer water accumulation situation and the target three-dimensional stratum model; and the parting layer water accumulation amount is calculated based on the upper and lower boundary surfaces. In this case, compared with the parting layer water accumulation amount calculated without considering the water accumulation in the prior art, the present application obtains the required upper and lower boundary surfaces by comprehensively considering the parting layer water accumulation situation and the target three-dimensional stratum model, and then calculates the parting layer water accumulation amount, so that the sandstone parting layer water accumulation amount can be more accurately calculated.

[0027] Additional aspects and advantages of the present application will be set forth in part in the description that follows, and in part will become apparent to those skilled in the art upon examination of the following and / or can be learned by practice of the present application. BRIEF DESCRIPTION OF DRAWINGS

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

[0029] Figure 1 A flowchart of a sandstone aquifer parting layer water accumulation amount calculation method provided by an embodiment of the present application;

[0030] Figure 2 A rock layer schematic diagram provided by an embodiment of the present application;

[0031] Figure 3 A Python fitted rock layer subsidence surface diagram provided by an embodiment of the present application;

[0032] Figure 4 A schematic diagram of a method for calculating parting layer volume using Python provided by an embodiment of the present application;

[0033] Figure 5 A block diagram of a sandstone aquifer parting layer water accumulation amount calculation system provided by an embodiment of the present application. DETAILED DESCRIPTION

[0034] Embodiments of the present application are described in detail below, examples of which are shown in the accompanying drawings, in which the same or similar reference numerals represent the same or similar elements or elements having the same or similar functions throughout. The embodiments described below by reference to the accompanying drawings are exemplary and are intended to explain the present application, and cannot be understood as limiting the present application.

[0035] The sandstone aquifer parting layer water accumulation amount calculation method and system of the embodiments of the present application are described below with reference to the accompanying drawings.

[0036] The water in bed separation water mine drainage capacity needs to be determined according to the predicted water in bed separation volume, the water in bed separation water accidents of the existing mines threatened by sandstone and gravel are mainly referenced to predict the water inrush volume, the water in bed separation volume empirical formula derived based on the thin plate small deflection theory does not consider the influence of accumulated water on the deflection of the rock stratum, and the calculation result is often greatly deviated from the actual situation.

[0037] Based on this, the sandstone and gravel aquifer water in bed separation accumulated water volume calculation method is provided to more accurately calculate the sandstone and gravel water in bed separation accumulated water volume.

[0038] Figure 1 A flowchart of the sandstone and gravel aquifer water in bed separation accumulated water volume calculation method provided by the embodiment of the present application is shown.

[0039] As shown in Figure 1 , the sandstone and gravel aquifer water in bed separation accumulated water volume calculation method includes the following steps:

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

[0041] In step S101, the water in bed separation position is determined based on the parameter information of the rock stratum above the coal seam and the goaf, including: determining whether the sandstone and gravel aquifer range contains a key layer based on the parameter information of the rock stratum above the coal seam by using a key layer discriminant; determining whether the goaf satisfies the space condition for developing water in bed separation based on the parameter information of the rock stratum of the water flowing fractured zone; if the space condition for developing water in bed separation is satisfied and the sandstone and gravel aquifer range contains a key layer, the water in bed separation position is determined in combination with the softening coefficient and the water resistance 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 gravity of each rock stratum between the coal seam and the ground, etc. The parameter information of the rock stratum of the water flowing fractured zone includes the original thickness and dilatancy 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 rock stratum type, the water abundance of the aquifer, the water level of the aquifer, etc.

[0044] Specifically, Figure 2 A rock stratum diagram is provided by the embodiment of the present application. As shown in Figure 2 , the underground coal seam to the ground includes multiple rock strata. After the coal seam is mined in the advancing direction of the working face, the mined part forms a goaf, the roof rock stratum is affected by the mining pressure, and gradually forms a caving zone, and a fractured zone and a curved subsidence zone are formed above the caving zone. Among them, the caving zone and the fractured zone form a water flowing fractured zone.

[0045] The position of the delamination of the sand and gravel aquifer threatening water inrush is determined based on the parameter information of the coal seam and the rock strata above the goaf. The position of the delamination is determined mainly from two conditions. The first condition is to determine whether the sand and gravel aquifer range between the coal seam and the ground has a delamination developed stratum structure, that is, to determine whether the sand and gravel aquifer range contains a key stratum. The second condition is to determine whether the goaf satisfies the space condition of the delamination development. If the two conditions are met at the same time, it is indicated that the delamination is generated in the sand and gravel aquifer range.

[0046] For the first condition, all the key strata between the coal seam and the ground are determined by a key stratum discriminant, wherein the key stratum discriminant satisfies:

[0047]

[0048] In the formula, E n+1 is the elastic modulus of the n+1 layer of rock strata above the coal seam, h n+1 is the thickness of the n+1 layer of rock strata above the coal seam, γ i is the specific gravity of the i layer of rock strata above the coal seam, h i is the thickness of the i layer of rock strata above the coal seam, represents the cumulative sum of the product of the specific gravity and the thickness of the first layer to the n layer of rock strata above the coal seam. γ n+1 is the specific gravity of the n+1 layer of rock strata above the coal seam, E i is the elastic modulus of the i layer of rock strata above the coal seam, h i is the thickness of the i layer of rock strata above the coal seam. Starting from the first layer of rock strata above the coal seam, when the formula (1) is satisfied for the first time, the n+1 layer can be determined as the key stratum, which is the first key stratum above the coal seam. Starting from the first key stratum, the formula (2) is used 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, it is determined whether the sand and gravel aquifer range contains a key stratum based on the sand and gravel aquifer range. If yes, the key stratum contained in the sand and gravel aquifer range is screened out.

[0050] For the second condition, a development discriminant is used to determine whether the goaf satisfies the space condition of the delamination development, wherein the development discriminant 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 dilatancy coefficient of the i layer of rock strata in the water-conducting fractured zone after the coal seam is mined, d iis the original thickness of the i-th stratum in the water flowing fractured zone after coal mining.∑(K i -1)d i The sum of the layers is the number of all strata in the water flowing fractured zone after coal mining. If both conditions (1) and (2) are met, it means that the separation occurs in the sand and gravel aquifer.

[0053] After determining that the separation occurs in the sand and gravel aquifer, the separation range is analyzed by the rock mechanics characteristics of the strata. Considering that the upper and lower boundary layers of the separation should be hard rock and soft rock respectively, the soft rock should have a certain water resistance capacity, and the sand and gravel aquifer is mainly composed of interbedded sandstone and conglomerate, therefore the development of the sand and gravel aquifer separation is usually a regional separation spanning multiple strata. The target key layer is determined according to the key layers in the sand and gravel aquifer, which is the key layer closest to the ground among the key layers in the sand and gravel aquifer. From the target key layer, the layers are analyzed layer by layer downward until the softening coefficient and water resistance capacity of a stratum meet the requirements, and the stratum with the softening coefficient and water resistance capacity meeting the requirements is taken as the lower boundary layer of the separation. The softening coefficient and water resistance capacity meeting the requirements means that the softening coefficient is less than or equal to a set coefficient threshold (for example, 0.75 MPa), and the relative water resistance (i.e. the stratum is a mudstone type). The target key layer is taken as the upper boundary layer of the separation.

[0054] The water abundance of the aquifer is determined according to the unit water inflow q of the borehole, and the separation water accumulation of the sand and gravel aquifer is determined based on the water abundance of the aquifer, the water level of the aquifer and the position of the separation. The water accumulation is divided into two cases: full water accumulation and incomplete water accumulation. When the water abundance of the aquifer is strong or extremely strong (for example, the unit water inflow of the borehole meets 1.0 L / (s·m)<q≤5.0 L / (s·m) for strong water abundance; q>5.0 L / (s·m) for extremely strong water abundance), the water level of the aquifer is higher than the determined position of the separation (here it means that the water level of the aquifer is higher than the upper boundary layer of the separation), and the water level is continuously and slowly decreasing, it can be determined that the separation has accumulated full water, otherwise the separation has not accumulated full water.

[0055] In step S102, for the region between the coal seam and the ground, the 3DEC software is used to simulate the deformation and failure of the strata to obtain the target three-dimensional stratum model.

[0056] Specifically, in step S102, for the area between the coal seam and the ground, an initial three-dimensional stratum model of the overburden rock of the area is established by using 3DEC (3DimensionDistinct Element Code) software, wherein the initial three-dimensional stratum model is established according to the stratum data actually exposed by drilling in the vertical direction, a random generated Voronoi polygon is used to simulate a gravel block with a random shape and size, the size of the gravel block is limited by setting a range of edge length, the inclination width of the initial three-dimensional stratum model in the horizontal direction should be consistent with the actual engineering background and more than a reserved set length (for example, 30 m) of surrounding rock in front and back, in order to accurately simulate the law of separation to closure, the length of the initial three-dimensional stratum model in the strike direction should not be less than twice the key layer breakage distance in the range of the sandy gravel aquifer (if there are multiple key layers, the maximum value of the breakage distance is taken), wherein the calculation formula of the key layer breakage distance is:

[0057]

[0058] In the formula, l is the key layer breakage distance, h is the calculated thickness of the key layer, R t is the calculated tensile strength of the key layer, and q is the calculated load of the overburden rock layer of the key layer.

[0059] Considering that the initial three-dimensional stratum model is a simulation model for the rock stratum above the unexplored coal seam, in order to better simulate the deformation and failure of the overburden rock after mining, the mechanical behavior of the rock is described by using the Mohr-Coulomb model, the mechanical behavior of the rock joint is described by using the Coulomb slip criterion, and the yield fracture and collapse are performed. Thus, the deformation and failure of the overburden rock are simulated by simulating excavation to obtain a target three-dimensional stratum model.

[0060] In step S103, the upper and lower boundary curved surfaces in the separation position are obtained based on the separation water accumulation condition and the target three-dimensional stratum model.

[0061] In step S103, the upper and lower boundary curved surfaces in the separation position are obtained based on the separation water accumulation condition and the target three-dimensional stratum model, including: fitting the upper and lower boundary layers in the separation position in the target three-dimensional stratum model to obtain a theoretical upper boundary curved surface and a theoretical lower boundary curved surface; calculating a corrected upper boundary curved surface and a corrected lower boundary curved surface for the separation water accumulation condition; if the separation is full of water, the theoretical upper boundary curved surface is the target upper boundary curved surface, and the corrected lower boundary curved surface is the target lower boundary curved surface; if the separation is not full of water, the corrected upper boundary curved surface is the target upper boundary curved surface, and the theoretical lower boundary curved surface is the target lower boundary curved surface.

[0062] Specifically, the process of obtaining the theoretical upper boundary curved surface and the theoretical lower boundary curved 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 position of the separation layer upper and lower boundary layer determined in step S101, use conditional statements to extract the id (identification number) of the gravel blocks of the upper and lower boundary layer, and store them in the list list1, list2 in python, extract the id of all gravel blocks in the separation layer range (i.e. between the upper and lower boundary layers) and store them in list3.

[0064] 12) Execute the excavation command to simulate the deformation and failure of the overburden rock, obtain the target three-dimensional stratum model, and update list1, list2 and list3 based on the id of the gravel blocks in the target three-dimensional stratum model. Traverse list1 and list2 respectively to obtain the position coordinates of each gravel block after the stratum moves due to excavation, and fit the settlement surface w1(x,y) and the settlement surface w2(x,y) according to the three-dimensional coordinates. The settlement surface w1(x,y) is the theoretical lower boundary surface, and the settlement surface w2(x,y) is the theoretical upper boundary surface. Thus, numerical planning is carried out in 3DEC software using Python, which can realize the surface fitting of the separation layer upper and lower boundary layers determined in step S101. Figure 3 The Python fitted stratum settlement surface diagram provided by the embodiment of the present application. Figure 3 A stratum settlement surface diagram of the entire sandy gravel aquifer fitted based on the target three-dimensional stratum model is shown. w is the settlement surface of any layer in the sandy gravel aquifer.

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

[0066] 21) Since the influence of water accumulation in the separation layer on the settlement of the lower boundary of the separation layer is not considered in the target three-dimensional stratum model, the influence of water accumulation in the separation layer on the settlement of the lower boundary of the separation layer is considered, and for the separation layer upper and lower boundary layers above the goaf, the stratum deflection micro-surface equation under the separation layer water is solved to obtain the stratum deflection surface, and the lower boundary of the separation layer is corrected, so the corrected lower boundary surface satisfies:

[0067]

[0068] In the formula, Where q is the transverse load on the lower boundary layer of the separation layer above the goaf, m is the number of series, a and b are the length along the working face strike and the width along the inclination of the lower boundary layer of the separation layer above the goaf, D is the bending stiffness, p is the density of groundwater, g is the acceleration of gravity, and (x, y) is the coordinates of the gravel block.

[0069] 22) During the working face mining process, according to the mining progress and the water level change of the sand and gravel aquifer, the current separation water accumulation situation is dynamically analyzed, if the separation is not filled with water, the sand and gravel aquifer water level is used as the upper boundary, that is, the sand and gravel aquifer water level is used to obtain the corrected upper boundary surface.

[0070] The final required target upper boundary surface and target lower boundary surface are determined considering the separation water accumulation situation.

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

[0072] If the separation is not filled with water: the upper boundary is the aquifer water level, 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 theoretical lower boundary surface.

[0073] Step S104, the separation water accumulation amount is calculated based on the upper and lower boundary surfaces.

[0074] In step S104, the separation water accumulation amount is calculated based on the upper and lower boundary surfaces, including: calculating the total volume of the interlayer between the upper and lower boundaries in the separation position based on the target upper boundary surface and the target lower boundary surface; calculating the total volume of the gravel blocks in the interlayer between the upper and lower boundaries based on the target three-dimensional stratum model; the total volume of the interlayer is obtained by subtracting the total volume of the gravel blocks from the total volume of the interlayer.

[0075] The separation water accumulation amount satisfies:

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

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

[0078] Figure 4 is a Python calculation method for separation volume provided by the embodiment of the present application. As Figure 4As shown, the curved subsidence zone is formed on the fracture zone, the upper and lower boundary layers of the separation layer position are determined by using step S101, and the water level of the aquifer above the upper boundary of the separation layer affects the calculation of the separation layer water accumulation. The fitting subsidence surface is obtained by using Python for the upper boundary of the separation layer and the lower boundary of the separation layer, the target upper boundary surface and the target lower boundary surface are determined in combination with the separation layer water accumulation, the total volume of the gravel block between the upper and lower boundary layers is read, and the separation layer water accumulation is obtained in combination with the total volume and the total volume of the gravel block.

[0079] In order to realize the above-mentioned embodiment, the present application further provides a sandstone aquifer separation layer water accumulation calculation system.

[0080] Figure 5 A block diagram of a sandstone aquifer separation layer water accumulation calculation system provided by an embodiment of the present application.

[0081] As Figure 5 shown, the sandstone aquifer separation layer water accumulation calculation system 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 configured to determine the separation layer position based on the parameter information of the rock layer above the coal seam and the goaf.

[0083] The modeling module 12 is configured to simulate the rock layer deformation and failure for the region between the coal seam and the ground by using 3DEC software to obtain a target three-dimensional stratum model.

[0084] The fitting module 13 is configured to obtain the upper and lower boundary surfaces in the separation layer position based on the separation layer water accumulation and the target three-dimensional stratum model.

[0085] The calculation module 14 is configured to calculate the separation layer water accumulation based on the upper and lower boundary surfaces.

[0086] Further, in a possible implementation manner of the embodiment of the present application, the determination module 11 is specifically configured to: determine whether the sandstone aquifer range contains a key layer based on the parameter information of the rock layer above the coal seam by using a key layer discriminant; determine whether the goaf satisfies the spatial condition for developing the separation layer based on the parameter information of the rock layer of the water flowing fracture zone; and if the spatial condition for developing the separation layer is satisfied and the sandstone aquifer range contains the key layer, determine the separation layer position in combination with the softening coefficient and the water-resisting capacity of the rock layer.

[0087] Further, in a possible implementation manner of the embodiment of the present application, the fitting module 13 is specifically configured to: fit the theoretical upper boundary surface and the theoretical lower boundary surface based on the upper and lower boundary layers of the off-layer position in the target three-dimensional stratum model; calculate the corrected upper boundary surface and the corrected lower boundary surface for the off-layer water accumulation condition; if the off-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 off-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 application, the calculation module 14 is specifically configured to: calculate the total volume between the upper and lower boundary layers in the off-layer position based on the target upper boundary surface and the target lower boundary surface; calculate the total volume of the gravel block between the upper and lower boundary layers based on the target three-dimensional stratum model; and obtain the off-layer water accumulation volume by subtracting the total volume of the gravel block from the total volume between the upper and lower boundary layers.

[0089] It should be noted that the foregoing explanation of the sand and gravel aquifer off-layer water accumulation volume calculation method embodiment is also applicable to the sand and gravel aquifer off-layer water accumulation volume calculation system of the embodiment, which will not be described here.

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

[0091] Compared with the prior art, the sand and gravel aquifer off-layer water accumulation volume calculation method of the present application can better target the geological structure and rock mechanical characteristics, and has strong applicability, which has the following advantages:

[0092] The construction of the drainage system of the mine threatened by off-layer water needs to consider the off-layer gushing water volume, and under the influence of mining, the sand and gravel aquifer develops pores and fissures, and is connected with the possible off-layer "cavity" to form a good water storage space. At this time, the traditional off-layer volume calculation formula cannot accurately evaluate the volume of the off-layer water, and the off-layer volume calculation method based on 3DEC-Python of the present application can accurately calculate the fissure volume of the gravel.

[0093] The water in the separation layer is a dynamic process, and the application considers the dynamic development of the separation layer, calculates the water in the separation layer before mining of the working face, provides a basis for water prevention and control design of the working face, and dynamically corrects the water in the separation layer according to the mining progress and the water level change of the aquifer during the mining process, so that the mine can make real-time adjustment on the drainage system of the working face. The water in the separation layer calculation method is reliable, and can provide protection for safe mining of the mine.

[0094] In order to realize the above-mentioned embodiments, the application further provides an electronic device, comprising: a processor and a memory connected with the processor; the memory stores computer execution instructions; and the processor executes the computer execution instructions stored in the memory to realize the method provided by the above-mentioned embodiments.

[0095] In order to realize the above-mentioned embodiments, the application further provides a computer readable storage medium, the computer readable storage medium stores computer execution instructions, and the computer execution instructions are executed by the processor to realize the method provided by the above-mentioned embodiments.

[0096] In order to realize the above-mentioned embodiments, the application further provides a computer program product, comprising a computer program, the computer program is executed by the processor to realize the method provided by the above-mentioned embodiments.

[0097] In the foregoing embodiment description, the description of the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" means that the specific features, structures, materials or characteristics described in combination with the embodiment or example are included in at least one embodiment or example of the application. In the description, the illustrative description of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any suitable manner in any one or more embodiments or examples. In addition, the person skilled in the art can combine and combine the different embodiments or examples described in the description and the features of the different embodiments or examples without contradiction.

[0098] In addition, the terms "first", "second" are only for descriptive purposes, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of indicated technical features. Therefore, the features defined with "first", "second" can explicitly or implicitly include at least one feature. In the description of the application, the meaning of "multiple" is at least two, for example, two, three, etc., unless otherwise specifically limited.

[0099] Any process or method description in the flowchart or otherwise herein can be understood as representing a module, segment, or portion of code comprising one or more executable instructions for implementing custom logic functions or processes, and the scope of preferred embodiments of the invention includes additional implementations in which functions may be performed not in the order shown or discussed, including substantially simultaneously or in reverse order depending on the functions involved, as should be understood by those skilled in the art to which embodiments of the invention pertain.

[0100] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (such as a computer-based system, a processor-included system, or other system that can fetch and execute instructions from, an instruction execution system, apparatus, or device). For the purposes of this specification, "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transmit programs for use by, or in conjunction with, an instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of computer-readable media include: an electrical connection having one or more wires (electronic device), a portable computer disk drive (magnetic device), random access memory (RAM), read-only memory (ROM), erasable and editable read-only memory (EPROM or flash memory), fiber optic devices, and portable optical disc read-only memory (CDROM). Alternatively, the computer-readable medium may be paper or other suitable media on which the program can be printed, since the program can be obtained electronically, for example, by optically scanning the paper or other medium, followed by editing, interpreting, or otherwise processing as necessary, 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 memory and executed by a suitable instruction execution system. For example, if implemented in hardware as in another embodiment, it can be implemented using any of the following techniques known in the art, or a combination thereof: discrete logic circuits having logic gates for implementing logical functions on data signals, application-specific integrated circuits (ASICs) having suitable combinational logic gates, programmable gate arrays (PGAs), field-programmable gate arrays (FPGAs), etc.

[0102] Those skilled in the art of the present technology can understand that all or part of the steps carried out by the above-mentioned embodiment method can be completed by programs instructing related hardware, and the programs can be stored in a computer readable storage medium. When the program is executed, it includes one of the steps of the method embodiment or a combination thereof.

[0103] In addition, each functional unit in each embodiment of the present application can be integrated into one processing module, or each unit can exist physically alone, or two or more units can be integrated into one module. The integrated module can be realized in the form of hardware or in the form of a software functional module. When the integrated module is realized 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 storage medium mentioned above can be a read-only memory, a magnetic disk or an optical disk, etc. Although the embodiments of the present application have been shown and described above, it should be understood that the above-mentioned embodiments are exemplary and cannot be understood as limiting the present application, and those skilled in the art can make changes, modifications, replacements and variations to the above-mentioned embodiments within the scope of the present application.

Claims

1. A method for calculating the water accumulation in a sandstone or conglomerate aquifer, characterized in that, include: The location of the separation is determined based on the parameter information of the rock strata above the coal seam and the goaf; For the area between the coal seam and the surface, 3DEC software is used to simulate the deformation and failure of the rock strata to obtain the target three-dimensional stratum model; Based on the delamination water accumulation and the target 3D stratigraphic model, the upper and lower boundary surfaces at the delamination location are obtained, including: The theoretical boundary surface and the theoretical lower boundary surface are obtained by fitting the upper and lower boundary layers at the delamination locations in the target three-dimensional stratigraphic model. For cases of water accumulation due to delamination, the upper boundary surface and the lower boundary surface are calculated and corrected. If the delamination is filled with water, then theoretically the boundary surface is the target upper boundary surface, and the modified lower boundary surface is the target lower boundary surface; if the delamination is not filled with water, then the modified upper boundary surface is the target upper boundary surface, and the theoretical lower boundary surface is the target lower boundary surface. The calculation of the corrected upper boundary surface and the corrected lower boundary surface for the case of delamination and water accumulation includes: Considering the impact of water accumulation on the settlement of the lower boundary of the delamination layer, for the upper and lower boundary layers of the delamination layer in the goaf, the rock strata deflection micro-surface equation under the water accumulation in the delamination layer is solved to obtain the rock strata deflection surface. The lower boundary of the delamination layer is then corrected, and the corrected lower boundary surface satisfies the following: in, , q is the lateral load on the lower boundary layer of the goaf, m is the number of stages of development, a and b are the length and dip width of the lower boundary layer of the goaf along the working face, D is the bending stiffness, ρ is the density of groundwater, g is the gravitational acceleration, (x,y) are the coordinates of the gravel block, and h is the thickness of the key layer. The amount of water accumulation due to delamination is calculated based on the upper and lower boundary surfaces, including: Calculate the total interlayer volume between the upper and lower boundary layers at the delamination location based on the target upper boundary surface and the target lower boundary surface; Calculate the total volume of gravel between the upper and lower boundary layers based on the target three-dimensional stratigraphic model; The total interlayer volume minus the total gravel volume yields the delamination water accumulation, wherein the delamination water accumulation satisfies the following conditions: Among them, V f This refers to the amount of water accumulated outside the stratum. For the target upper boundary surface, Let V be the target lower boundary surface, and V be the total volume of the gravel block.

2. The method for calculating the amount of water accumulated in a sandstone or conglomerate aquifer after separation according to claim 1, characterized in that, The determination of the delamination location based on parameter information of the rock strata above the coal seam and goaf includes: Based on the parameter information of the rock strata above the coal seam, the key layer discriminant is used to determine whether there is a key layer within the sandstone and conglomerate aquifer. Determine whether the goaf meets the spatial conditions for delamination development based on the parameter information of the rock strata in the water-conducting fracture zone; If the spatial conditions for the development of abscission are met and a key layer is present within the aquifer of the sandstone and conglomerate, the location of the abscission can be determined by combining the softening coefficient and water-retaining capacity of the rock strata.

3. A system for calculating the water accumulation in a sandstone or conglomerate aquifer, characterized in that, include: The determination module is used to determine the location of the separation based on the parameter information of the rock strata above the coal seam and the goaf; The modeling module is used to simulate rock deformation and failure in the area between the coal seam and the ground using 3DEC software to obtain a target three-dimensional geological model. The fitting module is used to obtain the upper and lower boundary surfaces at the delamination location based on the delamination water accumulation situation and the target three-dimensional stratigraphic model. The fitting module is specifically used to: fit the upper and lower boundary layers at the delamination location in the target three-dimensional stratigraphic model to obtain the theoretical boundary surface and the theoretical lower boundary surface; calculate the corrected upper boundary surface and the corrected lower boundary surface for the delamination water accumulation situation; if the delamination 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 delamination 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. The calculation of the corrected upper boundary surface and the corrected lower boundary surface for the case of delamination and water accumulation includes: Considering the impact of water accumulation on the settlement of the lower boundary of the delamination layer, for the upper and lower boundary layers of the delamination layer in the goaf, the rock strata deflection micro-surface equation under the water accumulation in the delamination layer is solved to obtain the rock strata deflection surface. The lower boundary of the delamination layer is then corrected, and the corrected lower boundary surface satisfies the following: in, , q represents the lateral load on the lower boundary layer of the goaf, m is the order of expansion, a and b are the length and dip width of the lower boundary layer of the goaf along the working face, D is the bending stiffness, ρ is the density of groundwater, g is the acceleration due to gravity, (x, y) are the coordinates of the gravel block, and h is the thickness of the key layer. The calculation module is used to calculate the amount of water accumulation outside the layer based on the upper and lower boundary surfaces; The calculation module is specifically used for: calculating the total interlayer volume between the upper and lower boundary layers at the delamination location based on the target upper boundary surface and the target lower boundary surface; calculating the total volume of gravel between the upper and lower boundary layers based on the target three-dimensional stratigraphic model; and subtracting the total volume of gravel from the total interlayer volume to obtain the delamination water accumulation, wherein the delamination water accumulation satisfies: Among them, V f This refers to the amount of water accumulated outside the stratum. For the target upper boundary surface, Let V be the target lower boundary surface, and V be the total volume of the gravel block.

4. The system for calculating the water accumulation in sandstone and conglomerate aquifers according to claim 3, characterized in that, The determination module is specifically used to: determine whether a key layer is present within the sandstone aquifer based on the parameter information of the rock strata above the coal seam using the key layer discriminant; determine whether the goaf meets the spatial conditions for the development of delamination based on the parameter information of the rock strata in the water-conducting fracture zone; if the spatial conditions for the development of delamination are met and a key layer is present within the sandstone aquifer, then determine the location of the delamination by combining the softening coefficient and water-retaining capacity of the rock strata.

5. An electronic device, characterized in that, include: A processor, and a memory communicatively connected to the processor; The memory stores computer-executed instructions; The processor executes computer execution instructions stored in the memory to implement the method as described in any one of claims 1-2.

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

Images