Method and device for predicting water quality characteristics of mine water and computer equipment

Abstract

The invention relates to the technical field of prediction of mine water quality characteristics, and particularly provides a prediction method and device for mine water quality characteristics and computer equipment, and the method comprises the steps: obtaining coal mine geological drilling data to determine a water-bearing layer and a water-resisting layer; various water quality parameter concentrations of all the aquifers are obtained through measurement based on the underground water samples of all the aquifers; calculating the development height of the water flowing fractured zone based on the coal seam mining thickness; based on the development height of the water-flowing fractured zone, determining a target water-bearing stratum set uncovered by the water-flowing fractured zone, and calculating the water inflow of each target water-bearing stratum in the target water-bearing stratum set into the well so as to obtain the total water inflow of the mine; and based on the water inflow of each target aquifer, the total water inflow of the mine and various water quality parameter concentrations of each target aquifer, calculating to obtain predicted concentrations of water quality parameters included in the mine water. By utilizing the method, the mine water quality prediction precision can be improved.

Description

Technical Field

[0001] The present invention relates to the technical field of prediction of mine water quality characteristics, and in particular to a method, device and computer equipment for predicting mine water quality characteristics. Background Art

[0002] Coal is an important strategic resource in my country, accounting for 56.2% of the country's energy consumption, and it is expected that the country's energy consumption will still be dominated by coal in the future. In the past 10 years, the strategic focus of my country's coal industry has accelerated its westward shift. At present, the output of raw coal in the western region accounts for more than 50% of the country's total coal output. However, affected by the destruction of the overburden structure caused by coal seam mining, the groundwater stored in the coal seam roof aquifer enters the mine through the mining-induced water-conducting fissures to form mine water. According to incomplete statistics, the total annual water inflow in the six northwestern provinces and regions is approximately 2.9 billion m 3 , and the quality of mine water is complex, with high mineralization, high turbidity, and high fluoride mine water widely distributed. According to current requirements, mine water must not be discharged without treatment. The quality of mine water discharged from coal mines should refer to the "Surface Water Environmental Quality Standard" (GB3838-2002). It must reach the national surface water environmental quality of Class III or above before it can be discharged. Therefore, accurate prediction of mine water quality characteristics is of great significance to the construction of mine water treatment facilities and the resource utilization of mine water. In coal mining practice, coal mining produces water-conducting fissures in the rock strata, directly exposing and approaching multiple aquifers. The water quality characteristics of different aquifers are generally quite different. Therefore, the mine water formed by the convergence of groundwater from multiple aquifers and mixing underground has many sources and complex components. It is urgent to develop a prediction method for mine water quality characteristics. Summary of the invention

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

[0004] To this end, the first object of the present invention is to propose a method for predicting mine water quality characteristics to improve the accuracy of mine water quality prediction.

[0005] The second object of the present invention is to provide a device for predicting the water quality characteristics of mine water.

[0006] A third object of the present invention is to provide a computer device.

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

[0008] To achieve the above object, the present invention provides a method for predicting the water quality characteristics of mine water in a first aspect, comprising:

[0009] Obtain geological drilling data for coal mines to identify aquifers and aquicludes;

[0010] Based on the groundwater samples of each aquifer, the concentrations of various water quality parameters of all aquifers were obtained;

[0011] Calculate the height of water-conducting fracture zones based on coal seam mining thickness;

[0012] Determine the target aquifer set exposed by the water-conducting fracture zone based on the development height of the water-conducting fracture zone, calculate the water inflow of each target aquifer in the target aquifer set into the mine, and then obtain the total water inflow of the mine;

[0013] The predicted concentration of the water quality parameters included in the mine water is calculated based on the water inflow of each target aquifer, the total water inflow of the mine and the concentration of various water quality parameters of each target aquifer.

[0014] In the method of the first aspect of the present invention, if the target aquifer set includes a first target aquifer that is completely penetrated by a water-conducting fracture zone; the first water inflow from the first target aquifer into the well is calculated based on the horizontal permeability coefficient of the first target aquifer, the elevation of the natural water level of the first target aquifer to the elevation of the corresponding bottom plate, the thickness of the first target aquifer, the head reduction value of the first target aquifer, the influence radius of the first target aquifer and the reference radius of the first target aquifer.

[0015] In the method of the first aspect of the present invention, the first water inflow from any first target aquifer into the well satisfies:

[0016]

[0017] In the formula, Q α is the first water yield of the first target aquifer; K α is the horizontal permeability coefficient of the first target aquifer; H α The elevation from the natural water level of the first target aquifer to the bottom plate of the first target aquifer; M α is the thickness of the first target aquifer; S α is the hydraulic head reduction value of the first target aquifer; R α0 is the influence radius of the first target aquifer; r α0 is the reference radius of the first target aquifer.

[0018] In the method of the first aspect of the present invention, if the target aquifer set includes a second target aquifer partially penetrated by a water-conducting fracture zone; the second target aquifer includes a penetrated section and a non-penetrated section; the water inflow of the penetrated section of the second target aquifer is calculated based on the horizontal permeability coefficient of the penetrated section, the thickness of the penetrated section, the head reduction value of the penetrated section, the influence radius of the penetrated section and the reference radius of the penetrated section in the second target aquifer; the water inflow of the non-penetrated section of the second target aquifer is calculated based on the vertical permeability coefficient of the non-penetrated section in the second target aquifer, the equivalent area of ​​the affected range of the water-conducting fracture zone, the natural elevation of the water head at the top of the second target aquifer to the elevation of the corresponding bottom plate, the elevation from the top of the water-conducting fracture zone to the bottom plate of the second target aquifer and the thickness of the non-penetrated section; based on the water inflow of the penetrated section and the water inflow of the non-penetrated section, the second water inflow into the well corresponding to the second target aquifer is obtained.

[0019] In the method of the first aspect of the present invention, the water inflow of the penetration section satisfies:

[0020]

[0021] In the formula, Q β K is the water yield of the second target aquifer through the section; β M is the horizontal permeability coefficient of the penetration section; β is the thickness of the through section; S β R is the head reduction value of the through section; β0 r is the influence radius of the penetration section; β0 The reference radius for the through segment.

[0022] In the method of the first aspect of the present invention, the water inflow of the non-penetrated section satisfies:

[0023]

[0024] In the formula, Q c K is the inflow of the unpenetrated section of the second target aquifer; v is the vertical permeability coefficient of the unpenetrated section; A is the equivalent area of ​​the water-conducting fracture zone; H c H is the natural elevation of the water head at the top of the second target aquifer to the elevation of the bottom plate of the second target aquifer; d M is the elevation from the top of the water-conducting fracture zone to the bottom of the second target aquifer; c The thickness of the section not penetrated.

[0025] In the method of the first aspect of the present invention, the predicted concentration of water quality parameters included in the mine water is calculated based on the water inflow of each target aquifer, the total water inflow of the mine and the concentration of various water quality parameters of each target aquifer, including: obtaining the proportion of water inflow of each target aquifer based on the water inflow of each target aquifer and the total water inflow of the mine; and calculating the predicted concentration of water quality parameters included in the mine water based on the concentration of various water quality parameters of each target aquifer and the proportion of water inflow of the corresponding layer.

[0026] To achieve the above-mentioned purpose, the second aspect of the present invention provides a prediction device for mine water quality characteristics, comprising:

[0027] Aquitard identification module, used to obtain geological drilling data of coal mines to determine aquifers and aquitards;

[0028] The water quality parameter calculation module of each aquifer in the water quality database is used to obtain the concentration of various water quality parameters of all aquifers based on the groundwater samples of each aquifer;

[0029] A water-conducting fracture zone development height prediction module is used to calculate the water-conducting fracture zone development height based on the coal seam mining thickness;

[0030] A water inflow prediction module for each aquifer, which is used to determine the target aquifer set exposed by the water-conducting fracture zone based on the development height of the water-conducting fracture zone, calculate the water inflow of each target aquifer in the target aquifer set into the well, and then obtain the total water inflow of the mine;

[0031] The coupling prediction module is used to calculate the predicted concentration of water quality parameters included in the mine water based on the water inflow of each target aquifer, the total water inflow of the mine and the concentration of various water quality parameters of each target aquifer.

[0032] To achieve the above-mentioned purpose, the third aspect of the present invention proposes a computer device, comprising: a processor, and a memory communicatively connected to the processor; the memory stores computer execution instructions; the processor executes the computer execution instructions stored in the memory to implement the method proposed in the first aspect of the present invention.

[0033] To achieve the above-mentioned purpose, the fourth aspect of the present invention proposes a computer-readable storage medium, in which computer-executable instructions are stored. When the computer-executable instructions are executed by a processor, they are used to implement the method proposed in the first aspect of the present invention.

[0034] The method, system, computer equipment and storage medium for predicting the water quality characteristics of mine water provided by the present invention obtain geological drilling data of coal mines to determine aquifers and aquicludes; obtain the concentrations of various water quality parameters of all aquifers based on the groundwater samples of each aquifer; calculate the development height of the water-conducting fracture zone based on the mining thickness of the coal seam; determine the target aquifer set exposed by the water-conducting fracture zone based on the development height of the water-conducting fracture zone, calculate the water inflow of each target aquifer in the target aquifer set into the well, and then obtain the total water inflow of the mine; calculate the predicted concentration of water quality parameters included in the mine water based on the water inflow of each target aquifer, the total water inflow of the mine and the concentrations of various water quality parameters of each target aquifer. In this case, a water quality database containing the concentrations of various water quality parameters of all aquifers is established, taking into account the amount of water inflowing into the well from each target aquifer in the target aquifer set exposed by the water-conducting fracture zone. The concentrations of water quality parameters included in mine water are predicted by coupling the water inflow of each target aquifer and the water quality parameter concentrations of the corresponding layer in the water quality database, thereby improving the accuracy of mine water quality prediction.

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

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

[0037] Figure 1 A schematic flow chart of a method for predicting the water quality characteristics of mine water provided by an embodiment of the present invention;

[0038] Figure 2 A combined characteristic map of a coal seam, aquifer, and water-conducting fracture zone provided by an embodiment of the present invention;

[0039] Figure 3 A curve showing the relationship between the development height of the water-conducting fracture zone and the mining height of the coal seam provided in an embodiment of the present invention;

[0040] Figure 4 A block diagram of a device for predicting mine water quality characteristics provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0041] Embodiments of the present invention are described in detail below, examples of which are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended to be used to explain the present invention, and should not be construed as limiting the present invention.

[0042] The following describes a method and system for predicting water quality characteristics of mine water according to an embodiment of the present invention with reference to the accompanying drawings.

[0043] The embodiment of the present invention provides a method for predicting the water quality characteristics of mine water, so as to improve the prediction accuracy of the water quality of mine water.

[0044] Figure 1 A schematic flow chart of a method for predicting mine water quality characteristics provided by an embodiment of the present invention.

[0045] like Figure 1 As shown, the prediction method of the mine water quality characteristics includes the following steps:

[0046] Step S101, obtaining coal mine geological drilling data to determine aquifers and aquicludes.

[0047] Specifically, in step S101, the coal mine geological drilling data of the target area is obtained, and the coal seam, aquifer and impermeable layer are delineated according to the drilling column chart.

[0048] Step S102, based on the groundwater samples of each aquifer, the concentrations of various water quality parameters of all aquifers are obtained.

[0049] In step S102, a water quality database is constructed based on the groundwater samples of each aquifer determined in step S101. The water quality database includes the concentrations of various water quality parameters of all aquifers. The types of water quality parameters include K + 、Na + , Ca 2+ Mg 2+ 、SO4 2- , Cl - 、HCO3 - 、F - , SS (suspended solids), mineralization, etc.

[0050] Specifically, the water quality database construction process includes: obtaining groundwater samples from different aquifers collected during on-site exploration of coal mines, the number of groundwater samples collected from each aquifer is a set number p, and the set number p is not less than 5, and for the groundwater samples from the same aquifer, the types of water quality parameters contained in each water sample and the concentrations of various water quality parameters are determined. Among them, for any water quality parameter in the groundwater of any aquifer, the average value of the test concentration of the water quality parameter of the set number p of water samples of the aquifer is used as the required concentration of the water quality parameter of the aquifer.

[0051] Take a water quality parameter m in the water sample of aquifer α1 as an example:

[0052]

[0053] In the formula, C α1mis the concentration of groundwater quality parameter m in aquifer α1; C m1 , C m2 , C m3 , ..., C mp They are the test concentrations of water quality parameter m corresponding to the p water samples collected from aquifer α1.

[0054] Step S103, calculating the height of the water-conducting fracture zone based on the coal seam mining thickness.

[0055] Specifically, in step S103, measured data of water-conducting fracture zones in the mined coal mines in the mining area are obtained, and the height H of the water-conducting fracture zones in the statistical data is D The relationship between the mining thickness of the coal seam and the thickness M is generally the empirical formula for estimating the water-conducting fracture zone:

[0056] H D =A×M+C

[0057] Where A is the fitting parameter and C is a constant. The formula is used to obtain the height of the water-conducting fracture zone in the target area based on the coal seam mining thickness in the target area.

[0058] Step S104, based on the height of the water-conducting fracture zone, determine the target aquifer set exposed by the water-conducting fracture zone, calculate the water inflow into the mine from each target aquifer in the target aquifer set, and then obtain the total water inflow of the mine.

[0059] In step S104, based on the coal mine geological data of the target area and the height of the water-conducting fracture zone in the target area calculated using the empirical formula of the water-conducting fracture zone, the distribution of the aquifer exposed by the water-conducting fracture zone in the target area is determined, and the distribution includes the position and number of the aquifers exposed by the water-conducting fracture zone. Each aquifer exposed by the water-conducting fracture zone can be regarded as a target aquifer, and all target aquifers constitute a target aquifer set.

[0060] In step S104, since the water-conducting fracture zone exposure is divided into two situations: completely penetrated by the water-conducting fracture zone and partially penetrated by the water-conducting fracture zone, the target aquifer set includes the first target aquifer completely penetrated by the water-conducting fracture zone and / or the second target aquifer partially penetrated by the water-conducting fracture zone.

[0061] Take the target aquifer set including n first target aquifers and 1 second target aquifer as an example. Figure 2 This is a combined characteristic diagram of coal seams, aquifers, and water-conducting fracture zones provided by an embodiment of the present invention. Figure 2 As shown, the n first target aquifers are aquifer α1, aquifer α2, ..., aquifer α above the coal seam. n . 1 The second target aquifer is aquifer z. H Dis the height of the water-conducting fracture zone. M is the mining thickness of the coal seam, M α1 is the thickness of aquifer α1, M α2 is the thickness of aquifer α2, …, M αn Aquifer α n Thickness, M β is the thickness of the penetration section in aquifer z, M c is the thickness of the unpenetrated section in aquifer z.

[0062] In step S104, if the target aquifer set includes a first target aquifer that is completely penetrated by a water-conducting fracture zone; the first water inflow from the first target aquifer into the well is calculated based on the horizontal permeability coefficient of the first target aquifer, the elevation of the natural water level of the first target aquifer to the elevation of the corresponding bottom plate, the thickness of the first target aquifer, the head reduction value of the first target aquifer, the influence radius of the first target aquifer and the reference radius of the first target aquifer.

[0063] The first inflow of water from any first target aquifer into the well (i.e. into the mining space) satisfies:

[0064]

[0065] In the formula, Q α is the first water yield of the first target aquifer α, in m 3 / d;K α is the horizontal permeability coefficient of the first target aquifer, in m / d; H α The elevation from the natural water level of the first target aquifer to the bottom plate of the first target aquifer, in meters; M α is the thickness of the first target aquifer, in m; S α is the head reduction value of the first target aquifer, in m; R α0 is the influence radius of the first target aquifer, in m; r α0 is the reference radius of the first target aquifer, in meters.

[0066]

[0067] R α0 =R α +r α0

[0068]

[0069] Where η α L is the influence coefficient of the pit shape, which is generally taken as 1.12; α is the length of the coal mining face, in m; B α is the width of the coal mining face, in m.

[0070] In step S104, if the target aquifer set includes a second target aquifer partially penetrated by a water-conducting fracture zone; the second target aquifer includes a penetrated section (i.e., the portion of the second target aquifer penetrated by the water-conducting fracture zone) and an unpenetrated section (i.e., the portion of the second target aquifer not penetrated by the water-conducting fracture zone). The amount of water flowing into the mining space from the second target aquifer consists of two parts, one part is the direct inflow of the penetrated section (also called the water inflow of the penetrated section), and the other part is the groundwater overflow of the unpenetrated section (also called the water inflow of the unpenetrated section).

[0071] In step S104, the water inflow of the penetrated section of the second target aquifer is calculated based on the horizontal permeability coefficient of the penetrated section, the thickness of the penetrated section, the head reduction value of the penetrated section, the influence radius of the penetrated section and the reference radius of the penetrated section in the second target aquifer; the water inflow of the non-penetrated section of the second target aquifer is calculated based on the vertical permeability coefficient of the non-penetrated section in the second target aquifer, the equivalent area of ​​the affected range of the water-conducting fracture zone, the natural elevation of the water head at the top of the second target aquifer to the elevation of the corresponding bottom plate, the elevation from the top of the water-conducting fracture zone to the bottom plate of the second target aquifer and the thickness of the non-penetrated section; and the second water inflow into the well corresponding to the second target aquifer is obtained based on the water inflow of the penetrated section and the water inflow of the non-penetrated section.

[0072] The water inflow of the penetration section meets the following requirements:

[0073]

[0074] In the formula, Q β is the water yield of the second target aquifer, in m 3 / d;K β is the horizontal permeability coefficient of the penetration section, in m / d; M β is the thickness of the penetration section, in m; S β is the head reduction value of the penetration section (the difference between the head value at the top of the water-conducting fracture zone and the bottom elevation of the second target aquifer), in meters; R β0 is the penetration radius, in m; r β0 It is the reference radius of the penetration segment, in meters.

[0075]

[0076] R β0 =R β +R β0

[0077]

[0078] Where η β L is the influence coefficient of the pit shape, which is generally taken as 1.12; βis the length of the coal mining face, in m; B β is the width of the coal mining face, in m.

[0079] The water inflow of the unpenetrated section meets the following requirements:

[0080]

[0081] In the formula, Q c is the inflow of the unpenetrated section of the second target aquifer, in m 3 / d;K v is the vertical permeability coefficient of the unpenetrated section, in m / d; A is the equivalent area of ​​the water-conducting fracture zone, in m 2 ;H c H is the natural elevation of the water head at the top of the second target aquifer to the elevation of the bottom plate of the second target aquifer, in meters; d M is the elevation from the top of the water-conducting fracture zone to the bottom of the second target aquifer, in meters; c is the thickness of the unpenetrated section, in m.

[0082] In step S104, if the target aquifer set includes a first number n of first target aquifers and one second target aquifer, the total water inflow of the mine satisfies:

[0083] Q=Q α1 +Q α2 ++Q αn +(Q β +Q c )

[0084] Where Q is the total water inflow of the mine, Q α1 , Q α2 ,...,Q αn are the first water yields corresponding to the n first target aquifers, Q β +Q c is the second water yield of the second target aquifer.

[0085] Step S105, based on the water inflow of each target aquifer, the total water inflow of the mine and the concentration of various water quality parameters of each target aquifer, the predicted concentration of the water quality parameters included in the mine water is calculated.

[0086] In step S105, the predicted concentration of water quality parameters included in the mine water is calculated based on the water inflow of each target aquifer, the total water inflow of the mine and the concentrations of various water quality parameters of each target aquifer, including: obtaining the proportion of water inflow of each target aquifer based on the water inflow of each target aquifer and the total water inflow of the mine; and calculating the predicted concentration of water quality parameters included in the mine water based on the concentrations of various water quality parameters of each target aquifer and the proportion of water inflow of the corresponding layer.

[0087] Taking the target aquifer set including n first target aquifers and 1 second target aquifer as an example, the water yield ratio of each target aquifer satisfies:

[0088]

[0089] Based on the proportion of water inflow of each target aquifer and the concentration of various water quality parameters in the corresponding layer, the water quality characteristics of the mine water can be estimated, that is, the concentrations corresponding to the different water quality parameters contained in the mine water can be estimated.

[0090] Taking a water quality parameter m contained in mine water as an example, the predicted concentration of the water quality parameter m of mine water satisfies:

[0091] C Km =C α1m δ α1 +C α2m δ α2 +...+C αnm δ αn +C zm (δ β +δ c )

[0092] In the formula, C Km is the predicted concentration of the water quality parameter m of the mine water; C α1m is the concentration of the groundwater quality parameter m in the aquifer α1 calculated in step S102; δ α1 is the proportion of water inflow from aquifer α1; C α2m is the calculated concentration of groundwater quality parameter m in aquifer α2; δ α2 is the proportion of water inflow from aquifer α2; C αnm is the calculated aquifer α n The concentration of groundwater quality parameter m; δ αn Aquifer α n The proportion of water inflow; C zm is the calculated concentration of groundwater quality parameter m in aquifer z; δ β is the proportion of water inflow in the aquifer z penetration section; δ c is the proportion of water inflow from the unpenetrated section of aquifer z. β +δ c is the proportion of water inflow from aquifer z.

[0093] The following takes an example in which the aquifer includes two aquifers, aquifer α1 and aquifer z, where aquifer α1 is the Yan'an Formation aquifer and aquifer z is the Zhiluo Formation aquifer, and describes the prediction method of the mine water quality characteristics of the present invention in combination with specific numerical values.

[0094] 1) Collect coal mine geological drilling data in the target area, and according to the drilling column chart, delineate the 2 coal as the main mining coal seam and the aquifer as the Yan'an Formation and Zhiluo Formation aquifer.

[0095] 2) Combined with the on-site investigation of the coal mine, 5 groundwater samples were collected from the Yan'an Formation and Zhiluo Formation aquifers, and the typical water quality parameters of the water samples (such as K + +Na + , Ca 2+ Mg 2+ 、SO4 2- , Cl - 、HCO3 - The average value of the tested concentrations of the water quality parameters of the five water samples is defined as the concentration of the corresponding water quality parameters of the aquifer, as shown in Table 1.

[0096] Table 1 Concentration of various water quality parameters in aquifers

[0097]

[0098] 3) Collect measured data of water-conducting fracture zones in mined coal mines within the mining area, as shown in Table 2.

[0099] Table 2 Statistical data of measured data of water-conducting fracture zones in the mining area

[0100] Coal Mine No. Working surface Mining height(m) Height of water-conducting fracture zone (m) 1 31104 3.0 68.30 2 42101 3.3 55.02 3 22101 5.5 107.49 4 20104 5.0 135.40 5 30101 4.50 106.10 6 2304 3.5 90.55 7 122106 6.0 137.63 8 112201 5.8 162.33 9 101 2.15 35.44 10 11201 3.2 104.36 11 11203 2.67 47.00 12 12403 2.26 49.32 13 12406 4.4 81.75 14 52306 7.3 137.32 15 N15203 4.5 127.80 16 N1209 4.8 130.60 17 311101 5.72 126.00 18 113103 4.95 132.90

[0101] In Table 2, the mining height is also called the coal seam mining thickness, and the height of the water-conducting fracture zone is also called the height of the water-conducting fracture zone development. The relationship between the height of the water-conducting fracture zone development (HD) and the coal seam mining thickness (M) is statistically analyzed. Figure 3 The relationship curve between the development height of the water-conducting fracture zone and the mining height of the coal seam provided in the embodiment of the present invention is shown in FIG. Figure 3 As shown in Figure 2, the empirical formula for estimating the water-conducting fracture zone is: D =22.932×M+1.8891. Determination coefficient R 2 =0.7474.

[0102] 4) The main coal seam 2 of the target area is 4m thick, the overlying Yan'an Formation is 32m thick, and the Zhiluo Formation is 90m thick. The empirical formula is used to obtain the height of the water-conducting fracture zone in the target area based on the coal seam mining thickness (i.e., M = 4m). D目标 =22.932×4+1.8891=93.62(m). It can be judged that the water-conducting fracture zone completely communicates with the Yan'an Formation (i.e., it penetrates 32m) and partially communicates with the Zhiluo Formation (i.e., the length of the penetrated section is 61.62m, and the length of the unpenetrated section is 28.38m). That is, the target aquifer set includes 1 first target aquifer and 1 second target aquifer. The first target aquifer is the Yan'an Formation aquifer, and the second target aquifer is the Zhiluo Formation aquifer.

[0103] 5) Since the water-conducting fracture zone completely penetrates the Yan'an Formation aquifer, referring to the calculation formula of the first water yield of the first target aquifer, the first water yield of the Yan'an Formation aquifer (i.e., aquifer α1) is:

[0104]

[0105] R α10 =R α1 +r α10 =42.87+1670.76=1713.63(m)

[0106]

[0107] In the formula, Q α1 is the first inflow of the Yan'an Formation aquifer; K α1 is the horizontal permeability coefficient of the Yan'an Formation aquifer; H α1 The elevation from the natural water level of the Yan'an Formation aquifer to the bottom of the first target aquifer; M α1 is the thickness of the Yan'an Formation aquifer; S α1 is the head reduction value of the Yan'an Formation aquifer; R α10 is the influence radius of the Yan'an Formation aquifer; r α10 is the reference radius of the Yan'an Formation aquifer.

[0108] Since the water-conducting fracture zone does not completely penetrate the Zhiluo Formation aquifer (the length of the penetrated section is 61.62m, and the length of the unpenetrated section is 28.38m), the water inflow of the penetrated section (61.62m) of the water-conducting fracture zone meets the following requirements:

[0109]

[0110] R β0 =R β +r β0 =134.65+1670.76=1805.41(m)

[0111]

[0112] The water inflow of the unpenetrated section (28.38m) of the water-conducting fracture zone is:

[0113]

[0114] The total water inflow of the mine meets the following requirements:

[0115] Q=Q α1 +Q β +Q c =2234.67+12943.36+9456.25=24634.28(m3 / d)

[0116] The proportion of penetrated sections and unpenetrated sections in the Yan'an Formation aquifer and the Zhiluo Formation aquifer is as follows:

[0117]

[0118] The water quality characteristics of mine water can be estimated from this, and the water quality parameter SO4 2- The concentration is as follows:

[0119]

[0120] Similarly, the concentrations of other water quality parameters can also be predicted, as shown in Table 3 below.

[0121] Table 3 Prediction values ​​of mine water quality parameters

[0122]

[0123] In order to implement the above embodiment, the present invention also proposes a device for predicting the water quality characteristics of mine water.

[0124] Figure 4 A block diagram of a device for predicting mine water quality characteristics provided by an embodiment of the present invention.

[0125] like Figure 4 As shown, the prediction device for the mine water quality characteristics includes an aquiclude identification module, a water quality parameter calculation module for each aquifer in the water quality database, a water-conducting fracture zone development height prediction module, a water yield prediction module for each aquifer, and a coupling prediction module, wherein:

[0126] Aquitard identification module, used to obtain geological drilling data of coal mines to determine aquifers and aquitards;

[0127] The water quality parameter calculation module of each aquifer in the water quality database is used to obtain the concentration of various water quality parameters of all aquifers based on the groundwater samples of each aquifer;

[0128] A water-conducting fracture zone development height prediction module is used to calculate the water-conducting fracture zone development height based on the coal seam mining thickness;

[0129] The water inflow prediction module for each aquifer is used to determine the target aquifer set exposed by the water-conducting fracture zone based on the development height of the water-conducting fracture zone, calculate the water inflow of each target aquifer in the target aquifer set into the mine, and then obtain the total water inflow of the mine;

[0130] The coupled prediction module is used to calculate the predicted concentration of water quality parameters included in the mine water based on the water inflow of each target aquifer, the total water inflow of the mine and the concentration of various water quality parameters of each target aquifer.

[0131] Furthermore, in a possible implementation of an embodiment of the present invention, in each aquifer inflow prediction module, if the target aquifer set includes a first target aquifer that is completely penetrated by a water-conducting fracture zone; based on the horizontal permeability coefficient of the first target aquifer, the elevation of the natural water level of the first target aquifer to the elevation of the corresponding bottom plate, the thickness of the first target aquifer, the head reduction value of the first target aquifer, the influence radius of the first target aquifer and the reference radius of the first target aquifer, the first water inflow of the first target aquifer into the well is calculated.

[0132] Furthermore, in a possible implementation of an embodiment of the present invention, in each aquifer yield prediction module, if the target aquifer set includes a second target aquifer partially penetrated by a water-conducting fracture zone; the second target aquifer includes a penetrated section and a non-penetrated section; the penetrated section yield of the second target aquifer is calculated based on the horizontal permeability coefficient of the penetrated section, the thickness of the penetrated section, the head reduction value of the penetrated section, the influence radius of the penetrated section and the reference radius of the penetrated section in the second target aquifer; the non-penetrated section yield of the second target aquifer is calculated based on the vertical permeability coefficient of the non-penetrated section in the second target aquifer, the equivalent area of ​​the affected range of the water-conducting fracture zone, the natural elevation of the water head at the top of the second target aquifer to the elevation of the corresponding bottom plate, the elevation from the top of the water-conducting fracture zone to the bottom plate of the second target aquifer and the thickness of the non-penetrated section; the second yield of the corresponding second target aquifer flowing into the well is obtained based on the penetrated section yield and the non-penetrated section yield.

[0133] Furthermore, in a possible implementation of the embodiment of the present invention, the coupling prediction module includes a water yield proportion prediction module for each aquifer and a mine water quality characteristics prediction module.

[0134] Furthermore, in a possible implementation of an embodiment of the present invention, a water yield ratio prediction module for each aquifer is used to obtain the water yield ratio of each target aquifer based on the water yield of each target aquifer and the total water yield of the mine.

[0135] Furthermore, in a possible implementation of an embodiment of the present invention, a mine water quality characteristic prediction module is used to calculate the predicted concentration of water quality parameters included in the mine water based on the concentration of various water quality parameters of each target aquifer and the proportion of water inflow of the corresponding layer.

[0136] It should be noted that the aforementioned explanation of the embodiment of the method for predicting the water quality characteristics of mine water is also applicable to the device for predicting the water quality characteristics of mine water in this embodiment, and will not be repeated here.

[0137] In the embodiment of the present invention, the aquifer and aquiclude are determined by obtaining geological drilling data of the coal mine; the concentrations of various water quality parameters of all aquifers are obtained based on the groundwater samples of each aquifer; the height of the water-conducting fracture zone is calculated based on the mining thickness of the coal seam; the target aquifer set exposed by the water-conducting fracture zone is determined based on the height of the water-conducting fracture zone, and the amount of water flowing into the well from each target aquifer in the target aquifer set is calculated to obtain the total water flow of the mine; the predicted concentration of the water quality parameters included in the mine water is calculated based on the water flow of each target aquifer, the total water flow of the mine and the concentrations of various water quality parameters of each target aquifer. In this case, a water quality database containing the concentrations of various water quality parameters of all aquifers is established, taking into account the amount of water flowing into the well from each target aquifer in the target aquifer set exposed by the water-conducting fracture zone, and the concentration of the water quality parameters included in the mine water is predicted by coupling the water flow of each target aquifer and the water quality parameter concentration of the corresponding layer in the water quality database, thereby improving the accuracy of the mine water quality prediction.

[0138] The method and device of the present invention are a new type of mine water quality prediction method and device. By calculating the amount of groundwater from each aquifer entering the coal mine mining space, combined with the water quality characteristics of each aquifer, based on the mass balance principle, the water quality characteristics of mine water are predicted, and efficient and accurate prediction of mine water quality is achieved, which solves the problems of difficulty in predicting mine water quality and low prediction accuracy, and provides technical support for the construction of mine water treatment facilities and the resource utilization of mine water.

[0139] Compared with the prior art, the present invention has the following advantages: on the basis of establishing a groundwater quality database of each aquifer, the process of groundwater gushing into the well from the aquifer is fully considered, and the calculation method of the water gushing amount of aquifers with completely connected (i.e. completely penetrated) water-conducting fracture zones and the water gushing amount of aquifers with incomplete communication (i.e. incomplete penetration) of water-conducting fracture zones are distinguished; by coupling the amount of water gushing into the mining space from each aquifer and the water quality of each aquifer, a prediction method for the water quality characteristics of mine water is formed; the method is highly scientific and simple to operate, and is of great significance to the construction of mine water treatment facilities and the resource utilization of mine water.

[0140] In order to implement the above embodiments, the present invention also proposes a computer device, including: a processor, and a memory communicatively connected to the processor; the memory stores computer execution instructions; the processor executes the computer execution instructions stored in the memory to implement the method provided by the above embodiments.

[0141] In order to implement the above embodiments, the present invention further proposes a computer-readable storage medium, in which computer-executable instructions are stored. When the computer-executable instructions are executed by a processor, they are used to implement the methods provided by the above embodiments.

[0142] In order to implement the above embodiments, the present invention further provides a computer program product, including a computer program, which implements the methods provided by the above embodiments when executed by a processor.

[0143] The collection, storage, use, processing, transmission, provision and disclosure of user personal information involved in this invention shall comply with the relevant laws and regulations and shall not violate public order and good morals.

[0144] It should be noted that personal information from users should be collected for legitimate and reasonable purposes and should not be shared or sold outside of these legitimate uses. In addition, such collection / sharing should be carried out after receiving the user's informed consent, including but not limited to notifying the user to read the user agreement / user notice and sign the agreement / authorization including authorization of relevant user information before the user uses the function. In addition, any necessary steps should be taken to protect and safeguard access to such personal information data and ensure that others who have access to personal information data comply with its privacy policy and procedures.

[0145] The present invention contemplates providing embodiments where users can selectively block the use or access of personal information data. That is, the present disclosure contemplates providing hardware and / or software to prevent or block access to such personal information data. Once the personal information data is no longer needed, risks can be minimized by limiting data collection and deleting the data. In addition, when applicable, such personal information is de-identified to protect the privacy of the user.

[0146] In the description of the aforementioned embodiments, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction 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 do not necessarily refer 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, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, unless they are contradictory.

[0147] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Therefore, the features defined as "first" and "second" may explicitly or implicitly include at least one of the features. In the description of the present invention, the meaning of "plurality" is at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0148] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, segment or portion of code comprising one or more executable instructions for implementing the steps of a custom logical function or process, and the scope of the preferred embodiments of the present invention includes alternative implementations in which functions may not be performed in the order shown or discussed, including performing functions in a substantially simultaneous manner or in reverse order depending on the functions involved, which should be understood by technicians in the technical field to which the embodiments of the present invention belong.

[0149] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, device or apparatus (such as a computer-based system, a system including a processor, or other system that can fetch instructions from an instruction execution system, device or apparatus and execute the instructions), or in combination with these instruction execution systems, devices or apparatuses. For the purpose of this specification, "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by an instruction execution system, device or apparatus, or in combination with these instruction execution systems, devices or apparatuses. More specific examples of computer-readable media (a non-exhaustive list) include the following: an electrical connection with one or more wires (electronic device), a portable computer disk box (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), a fiber optic device, and a portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically, for example, by optically scanning the paper or other medium and then editing, interpreting or processing in other suitable ways if necessary, and then stored in a computer memory.

[0150] It should be understood that the various parts of the present invention can be implemented by hardware, software, firmware or a combination thereof. In the above-mentioned embodiments, a plurality of steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0151] A person skilled in the art may understand that all or part of the steps in the method for implementing the above-mentioned embodiment may be completed by instructing related hardware through a program, and the program may be stored in a computer-readable storage medium, which, when executed, includes one or a combination of the steps of the method embodiment.

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

[0153] The storage medium mentioned above may be a read-only memory, a magnetic disk or an optical disk, 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 cannot be understood as limiting the present invention. A person of ordinary skill in the art may change, modify, replace and modify the above embodiments within the scope of the present invention.

Claims

1. A method for predicting mine water quality characteristics, characterized in that: include: Obtain geological drilling data for coal mines to identify aquifers and aquicludes; Based on the groundwater samples of each aquifer, the concentrations of various water quality parameters of all aquifers were obtained; Calculate the height of water-conducting fracture zones based on coal seam mining thickness; Determine the target aquifer set exposed by the water-conducting fracture zone based on the development height of the water-conducting fracture zone, calculate the water inflow of each target aquifer in the target aquifer set into the mine, and then obtain the total water inflow of the mine; The predicted concentration of the water quality parameters included in the mine water is calculated based on the water inflow of each target aquifer, the total water inflow of the mine and the concentration of various water quality parameters of each target aquifer.

2. The method for predicting the mine water quality characteristics according to claim 1, characterized in that: If the target aquifer set includes a first target aquifer that is completely penetrated by a water-conducting fracture zone; the first water inflow from the first target aquifer into the well is calculated based on the horizontal permeability coefficient of the first target aquifer, the elevation of the natural water level of the first target aquifer to the elevation of the corresponding bottom plate, the thickness of the first target aquifer, the head reduction value of the first target aquifer, the influence radius of the first target aquifer and the reference radius of the first target aquifer.

3. The method for predicting the mine water quality characteristics according to claim 2, characterized in that: The first inflow of water from any first target aquifer into the well meets the following requirements: In the formula, Q α is the first water yield of the first target aquifer; K α is the horizontal permeability coefficient of the first target aquifer; H α The elevation from the natural water level of the first target aquifer to the bottom plate of the first target aquifer; M α is the thickness of the first target aquifer; S α is the hydraulic head reduction value of the first target aquifer; R α0 is the influence radius of the first target aquifer; r α0 is the reference radius of the first target aquifer.

4. The method for predicting the mine water quality characteristics according to claim 2, characterized in that: If the target aquifer set includes a second target aquifer partially penetrated by a water-conducting fracture zone; the second target aquifer includes a penetrated section and a non-penetrated section; the water inflow of the penetrated section of the second target aquifer is calculated based on the horizontal permeability coefficient of the penetrated section, the thickness of the penetrated section, the head reduction value of the penetrated section, the influence radius of the penetrated section and the reference radius of the penetrated section in the second target aquifer; the water inflow of the non-penetrated section of the second target aquifer is calculated based on the vertical permeability coefficient of the non-penetrated section in the second target aquifer, the equivalent area of ​​the affected range of the water-conducting fracture zone, the natural elevation of the water head at the top of the second target aquifer to the elevation of the corresponding bottom plate, the elevation from the top of the water-conducting fracture zone to the bottom plate of the second target aquifer and the thickness of the non-penetrated section; based on the water inflow of the penetrated section and the water inflow of the non-penetrated section, the second water inflow into the well corresponding to the second target aquifer is obtained.

5. The method for predicting the mine water quality characteristics according to claim 4, characterized in that: The water inflow of the penetration section meets the following requirements: In the formula, Q β K is the water yield of the second target aquifer through the section; β M is the horizontal permeability coefficient of the penetration section; β is the thickness of the through section; S β R is the head reduction value of the penetration section; β0 is the influence radius of the penetration section; γ β0 The reference radius for the through segment.

6. The method for predicting the mine water quality characteristics according to claim 4, characterized in that: The water inflow of the unpenetrated section meets the following requirements: In the formula, Q c K is the inflow of the unpenetrated section of the second target aquifer; v is the vertical permeability coefficient of the unpenetrated section; A is the equivalent area of ​​the water-conducting fracture zone; H c H is the natural elevation of the water head at the top of the second target aquifer to the elevation of the bottom plate of the second target aquifer; d M is the elevation from the top of the water-conducting fracture zone to the bottom of the second target aquifer; c The thickness of the section not penetrated.

7. The method for predicting mine water quality characteristics according to claim 1, characterized in that: The predicted concentration of the water quality parameters of the mine water is calculated based on the water inflow of each target aquifer, the total water inflow of the mine and the concentration of various water quality parameters of each target aquifer, including: Obtaining the water inflow ratio of each target aquifer based on the water inflow of each target aquifer and the total water inflow of the mine; The predicted concentration of water quality parameters included in mine water is calculated based on the concentration of various water quality parameters of each target aquifer and the proportion of water inflow of the corresponding layer.

8. A device for predicting mine water quality characteristics, characterized in that: include: Aquitard identification module, used to obtain geological drilling data of coal mines to determine aquifers and aquitards; The water quality parameter calculation module of each aquifer in the water quality database is used to obtain the concentration of various water quality parameters of all aquifers based on the groundwater samples of each aquifer; A water-conducting fracture zone development height prediction module is used to calculate the water-conducting fracture zone development height based on the coal seam mining thickness; A water inflow prediction module for each aquifer, which is used to determine the target aquifer set exposed by the water-conducting fracture zone based on the development height of the water-conducting fracture zone, calculate the water inflow of each target aquifer in the target aquifer set into the well, and then obtain the total water inflow of the mine; The coupling prediction module is used to calculate the predicted concentration of water quality parameters included in the mine water based on the water inflow of each target aquifer, the total water inflow of the mine and the concentration of various water quality parameters of each target aquifer.

9. A computer 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 7.

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 7 when executed by a processor.

Images