Working face effluent treatment effect evaluation method

Abstract

The invention relates to the technical field of grouting engineering, in particular to a working face effluent treatment effect evaluation method. The invention provides a comprehensive grouting project evaluation method, and the method achieves the deep understanding and quantitative evaluation of the grouting project effect through systematic steps and fine data analysis. Through detailed data definition and collection, the hydrogeological structure data of the grouting drilling hole and the stratum can be accurately obtained; meanwhile, the generated hydrogeological structure data matrix and the generated accumulated grouting quantity matrix provide convenience for further weighting processing and matrix merging; through matrix merging and center matrix calculation, comprehensive analysis of the grouting amount and the hydrogeological structure data is achieved, and the correlation between different hydrogeological structure data, the correlation between the grouting amount and the hydrogeological structure data and the correlation between the grouting amount of different grouting drill holes are quantitatively analyzed.

Description

Technical Field

[0001] The invention relates to the technical field of grouting engineering, in particular to a method for evaluating the effect of water discharge treatment on a working face. Background Art

[0002] Water in the working face is often the result of a combination of geological and hydrological factors. In underground projects such as coal mines and tunnels, water in the working face can originate from water in geological structures such as water-bearing rock strata, faults, collapse columns, and tectonic fracture zones. Water in the working face can also flow through existing goafs. Furthermore, high mine pressure or hydrostatic pressure can cause expansion of the tunnel floor, damage surrounding rock, and cracks, which can then flow through water-bearing strata and cause water intrusion.

[0003] Grouting is a technique that involves injecting slurry into the ground to reinforce it and prevent water leaks. Grouting is commonly used in underground projects such as coal mines and tunnels to address water intrusion at the working face. Once injected, the grouting material fills cracks and voids within the formation, forming a water barrier that prevents groundwater from seeping through.

[0004] Grouting, as an important geological engineering technique, is widely used in mining, tunnel construction, subway construction, water conservancy projects, and other fields. It is primarily used to reinforce strata, waterproof and plug leaks, and improve foundation bearing capacity. However, the effectiveness of grouting is often affected by a variety of factors, such as the hydrogeological conditions of the stratum, the selection and proportion of grouting materials, and grouting process parameters. Therefore, accurately evaluating grouting projects to ensure their effectiveness is achieved is a pressing issue in engineering practice.

[0005] In view of the above problems, it is necessary to propose a method for evaluating the effect of water discharge control on working faces. Summary of the Invention

[0006] The purpose of the present invention is to provide a method for evaluating the effect of water discharge treatment on a working face. The method achieves an in-depth understanding and quantitative evaluation of the effect of grouting engineering through systematic steps and sophisticated data analysis. Through detailed data definition and collection, the present invention can accurately obtain the hydrogeological structural data of the grouting borehole and the formation, providing a basis for subsequent analysis.

[0007] The purpose of the present invention can be achieved through the following technical solutions: A method for evaluating the effect of water discharge treatment at a working face comprises the following steps: Step 1: Hydrological data monitoring and analysis: Hydrological data includes working face hydrogeological structure data and working face natural hydrological condition data. Working face hydrogeological structure data includes mining area aquifer water level, borehole unit water inflow, and water depth in the open mining area. Natural hydrological condition data includes daily precipitation, surface runoff, and groundwater level. Hydrological data are used to generate the working face hydrogeological structure matrix and working face hydrological condition matrix, and extract the hydrogeological synergy index, the contribution index of each grouting borehole, and the grouting synergy degree. Step 2: Grouting Project Operation Evaluation: Real-time monitoring of the water level elevation of each monitoring well in the fourth aquifer in the working face, generating a water level elevation duration curve, and calculating groundwater loss parameters, treatment efficiency parameters, and treatment effect parameters. Simultaneously, real-time monitoring of the total grouting flow and average pore pressure of the grouting holes, generating a grouting flow and pore pressure duration curve, and calculating the material contribution parameter and pore pressure control parameter. Step 3: Economic benefit evaluation: Obtain the total footage of each grouting borehole and the footage data at different stages, as well as the final cumulative injection volume of the sealing material and the grouting loss volume, and calculate the economic evaluation parameters; Step 4: Environmental impact assessment: Obtain the leakage data of chemical substances treated at the water surface, including the COD and BOD values ​​of the fourth aquifer, and calculate the environmental impact parameters; Step 5. Comprehensive evaluation and report generation: Record all evaluation data, including the hydrogeological synergy index, the contribution index of each grouting borehole, the grouting synergy, the groundwater loss parameter, the treatment efficiency parameter, the treatment effect parameter, the material contribution parameter, the pore pressure control parameter, the economic evaluation parameter and the environmental impact parameter, generate an evaluation report, and highlight the evaluation data that exceeds the threshold.

[0008] Optionally, the step 1 specifically includes: S101, data definition and collection: number each grouting borehole, and collect the hydrogeological structure data of the stratum where the working face is located, including the water level of the aquifer in the mining area, the unit water inflow of the aquifer borehole, the water accumulation depth of the underground empty mining area, the rock stratum fracture ratio of the working face, the maximum fracture spacing of the rock stratum in the working face, and the minimum fracture spacing of the rock stratum in the working face. At the same time, the cumulative total injection volume of each grouting borehole is calculated; the hydrogeological structure data and the cumulative total injection volume of each grouting borehole are merged to generate a hydrogeological structure data matrix and a cumulative grouting volume matrix; S102, weighted processing of hydrogeological structure data: Calculating a locally weighted hydrogeological structure data matrix based on the elevation of the working face and the ground elevation, combined with the hydrogeological structure data matrix, daily precipitation, surface runoff, and groundwater level; S103, matrix merging and central matrix calculation: merging the cumulative grouting volume matrix and the locally weighted hydrogeological structure data matrix to obtain a hydrological comprehensive matrix, and calculating the central matrix; S104, covariance matrix calculation and feature extraction: calculate the covariance matrix of the central matrix, and extract the overall hydrogeological synergy index, the contribution index of each grouting borehole, and the grouting synergy degree from the covariance matrix of the central matrix.

[0009] Optionally, the S101 specifically includes: Number each grouting borehole, the number symbol is ; n is the total number of grouting holes; the data observation period is numbered, and the number symbol is ; ;in is the number of days included in the entire data observation period; Aquifer water level in the mining area Aquifer borehole unit water yield , Depth of water accumulation in underground mining area , rock fracture rate of working face , Maximum crack spacing of rock strata at the working face Minimum crack spacing between working face rock layers , all in mm; the fracture rate of the rock formation at the working face is a dimensionless parameter; The specific process of merging the hydrogeological structure data and the cumulative total injection volume of each grouting borehole to generate a data matrix is ​​as follows: Generate hydrogeological structure data matrix; ; Count the cumulative total injection volume of each grouting borehole i on different recording dates t , unit is t, generate cumulative grouting amount matrix .

[0010] Optionally, the specific process of S102 is: Obtain the elevation H of the working face and the ground elevation H0; collect the natural hydrological condition data of the working face environment during each data observation period, including daily precipitation , surface runoff and groundwater level ; Among them, the elevation of the working surface and ground elevation The unit is ; The unit of daily precipitation is ; The unit of surface runoff is ; The unit of groundwater level is ; By formula Calculate daily natural hydrological condition parameters ,in and All are preset weight factors; generate hydrological condition matrix ; Through the formula

[0011] Calculate the locally weighted hydrogeological structure data matrix V, where | is the augmented matrix dividing line.

[0012] Optionally, the specific process of S103 is: Merged cumulative grouting volume matrix and the locally weighted hydrogeological structure data matrix , and obtain the hydrological comprehensive matrix ,in is the augmented matrix dividing line; By formula: Calculate the center matrix ,in and is the weighted hydrogeological structure data matrix and cumulative grouting volume matrix The column mean vector of .

[0013] Optionally, the specific process of S104 is as follows: By formula calculate The covariance matrix of dimension ,in: ; where · is the matrix multiplication operator; is the central matrix The transposed matrix of is the geological structure synergy matrix, that is, the covariance matrix between different hydrogeological structure data; is the grouting volume-geological structure submatrix, that is, the covariance matrix between the grouting volume of different grouting boreholes and the hydrogeological structure data; is the grouting amount coordination sub-matrix, that is, the covariance matrix between the grouting amounts of different grouting boreholes; By formula Calculate the overall hydrogeological synergy index , each grouting drilling hole Contribution index and grouting coordination .

[0014] Optionally, the specific process of step 2 is as follows: During the grouting treatment process at the water outlet, the water level elevation in each monitoring well in the fourth aquifer is monitored in real time to generate a water level elevation duration curve for the fourth aquifer; feature extraction is performed on the water level elevation duration curve for the fourth aquifer, and the average water level elevation in each monitoring well is calculated to obtain an average water level curve; Get the average thickness of the fourth aquifer , in m, to obtain the time when the loss occurs Average water level elevation in the average water level curve , in m; get the time when the loss ends Average water level elevation in the average water level curve , in m; get the preset lead time before the loss occurs When the water level in the average water level curve is , unit is m; By formula Calculation of groundwater loss parameters , Governance efficiency parameters and governance effect parameters ,in 、 and All are preset influence coefficients; During the grouting treatment process at the water surface, the total grouting flow of all grouting holes is monitored in real time , the unit is and the average pore pressure of all grouting holes , the unit is MPa, generate the grouting flow rate and pore pressure duration curve; Extract features from the grouting flow rate and pore pressure duration curve to obtain the analytical formula for the grouting flow rate ; Get the maximum average pore pressure Time of maximum average pore pressure generation , through the formula ; Calculating material contribution parameters and pore pressure control parameters ,in is the preset influence coefficient, is the preset pore pressure threshold; is the cumulative grouting volume, in units of The unit of the maximum average pore pressure is Pa, and the unit of the maximum average pore pressure generation time is seconds; the material contribution parameter The unit is .

[0015] Optionally, the specific process of calculating the economic evaluation parameters in step 3 is as follows: Get each grouting borehole Total footage , one-step progress , Second opening footage and three-opening footage , the unit is m; Get each grouting borehole The final cumulative injection volume of sealing materials, including cement injection volume , unit is t; fine aggregate injection amount , unit is t; and water glass injection volume , unit is t; get each grouting borehole Grouting loss , the unit is ; By formula Calculation of economic evaluation parameters ,in 、 and is the preset weight factor.

[0016] Optionally, the specific process of calculating the environmental impact parameters in step 4 is: Obtain data on chemical seepage from surface water treatment, including the quaternary aquifer value and value , through the formula Calculation of environmental impact parameters ;in and The unit is ; is the cumulative grouting volume, in units of .

[0017] Compared with the prior art, the present invention has the following beneficial effects: 1. The present invention provides a comprehensive grouting engineering evaluation method, which achieves an in-depth understanding and quantitative evaluation of the grouting engineering effect through systematic steps and detailed data analysis; through detailed data definition and collection, the present invention can accurately obtain the hydrogeological structural data of the grouting borehole and the stratum, providing a basis for subsequent analysis. At the same time, the generated hydrogeological structural data matrix and cumulative grouting volume matrix provide convenience for further weighted processing and matrix merging; the present invention takes into account the impact of natural hydrological conditions on hydrogeological structural data, and through weighted processing, makes the data closer to the actual situation. At the same time, through matrix merging and central matrix calculation, a comprehensive analysis of grouting volume and hydrogeological structural data is achieved; 2. By calculating the covariance matrix, the present invention can quantitatively analyze the correlations between different hydrogeological structure data, between grouting volume and hydrogeological structure data, and between grouting volumes in different grouting boreholes. This helps to gain a deeper understanding of the improvement of hydrogeological structures by grouting projects and the mutual influence between grouting boreholes.

[0018] 3. The present invention can evaluate the operation of the grouting project in real time by monitoring the water level elevation of the fourth aquifer and the total grouting flow and average pore pressure of the grouting holes in real time. At the same time, by calculating the groundwater loss parameters, treatment efficiency parameters, treatment effect parameters, material contribution parameters and pore pressure control parameters, the effect and efficiency of the grouting project can be accurately evaluated. In addition, the present invention comprehensively considers the total footage of the grouting borehole and the cumulative injection volume of the sealing material, as well as the grouting loss volume. By calculating the economic evaluation parameters, the economic benefits of the grouting project can be comprehensively evaluated. By obtaining the leakage data of the chemical substances treated at the water surface and calculating the environmental impact parameters, the potential impact of the grouting project on the environment can be evaluated. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] In order to facilitate understanding by those skilled in the art, the present invention is further described below with reference to the accompanying drawings: Figure 1 is a flow chart of the method of the present invention; Figure 2 This is a schematic diagram of a time course curve of the water level elevation of the fourth aquifer according to the present invention; Figure 3 It is a schematic diagram of the grouting flow rate and pore pressure duration curve of the present invention. DETAILED DESCRIPTION

[0020] The technical solutions of the present invention will be clearly and completely described below in conjunction with the embodiments. Obviously, the embodiments described are only some of the embodiments of the present invention, not all of them. All other embodiments obtained by ordinary technicians in this field based on the embodiments of the present invention without making any creative efforts shall fall within the scope of protection of the present invention.

[0021] It should be understood that the terms “include” and “comprising” used in the specification and claims of the present disclosure indicate the presence of described features, integers, steps, operations, elements and / or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or collections thereof.

[0022] It should also be understood that the terms used in this disclosure are for the purpose of describing particular embodiments only and are not intended to limit the disclosure. As used in this disclosure and the claims, the singular forms "a," "an," and "the" are intended to include the plural forms unless the context clearly indicates otherwise. It should also be further understood that the term "and / or" used in this disclosure and the claims refers to any and all possible combinations of one or more of the associated listed items, and includes these combinations; See also Figure 1 As shown, the working face water discharge treatment effect evaluation method of the present invention comprises the following steps: Step 1: Hydrological data monitoring and analysis; Collect, define, and process hydrological data for the working face area. This involves systematically collecting and weighting data on the mining area's hydrogeological structure and natural hydrological conditions. Data on each grouting borehole is collected and defined, along with hydrogeological data on the working face, including aquifer levels, borehole specific water yield, and water accumulation depth in the open-cut area. Natural hydrological data is also collected, including daily precipitation, surface runoff, and groundwater levels. Weighted processing of this data generates a hydrogeological structure and hydrological condition matrix, providing the foundation for grouting monitoring data analysis and model calculations in step two.

[0023] S101, data definition and collection; Number each grouting borehole, the number symbol is ; n is the total number of grouting holes; the data observation period is numbered, and the number symbol is ; ;in is the number of days included in the entire data observation period.

[0024] During each data observation period, hydrogeological structural data of the stratum where the working face is located are collected, including: Aquifer water level in the mining area Aquifer borehole unit water yield , Depth of water accumulation in underground mining area , rock fracture rate of working face , Maximum crack spacing of rock strata at the working face Minimum crack spacing between working face rock layers , all in mm; the fracture rate of the rock formation at the working face is a dimensionless parameter; The units of the water level of the aquifer in the mining area, the unit water yield of the aquifer borehole, the depth of water accumulation in the underground mining area, the maximum crack spacing of the rock formation at the working face, and the minimum crack spacing of the rock formation at the working face are all millimeters (mm). Among them, the fracture rate of the working face rock formation is a dimensionless parameter.

[0025] Generate hydrogeological structure data matrix .

[0026] Count the cumulative total injection volume of each grouting borehole i on different recording dates t , unit is t, generate cumulative grouting amount matrix .

[0027] S102, weighted processing of hydrogeological structure data; Obtain the elevation H of the working face and the ground elevation H0. Collect the natural hydrological condition data of the working face environment during each data observation period, including daily precipitation. , surface runoff and groundwater level .

[0028] Among them, the elevation of the working surface and ground elevation The unit is ; The unit of daily precipitation is ; The unit of surface runoff is ; The unit of groundwater level is ; By formula Calculate daily natural hydrological condition parameters ,in and All are preset weight factors; generate hydrological condition matrix .

[0029] By formula

[0030] Calculate the locally weighted hydrogeological structure data matrix V, where | is the augmented matrix dividing line.

[0031] It should be noted that the purpose of calculating local weights is to quantify the additional impact of environmental factors on hydrogeological structural data. Grouting, the primary method for managing water discharge from the working face, fills cracks in the coal seam at the working face, forcing water that would have seeped into the working face to remain in the original aquifer or flow into other underground runoff, thereby improving the hydrogeological structure. However, while quantifying changes in hydrogeological structure, changes in environmental factors such as precipitation and underground mining cannot be ignored.

[0032] S103, matrix merging and center matrix calculation; Merged cumulative grouting volume matrix and the locally weighted hydrogeological structure data matrix , and obtain the hydrological comprehensive matrix , where | is the augmented matrix dividing line.

[0033] Furthermore, through the formula: Calculate the center matrix ,in and is the weighted hydrogeological structure data matrix and cumulative grouting volume matrix The column mean vector of .

[0034] S104, covariance matrix calculation and feature extraction; By formula calculate The covariance matrix of dimension ,in: ; where · is the matrix multiplication operator; is the central matrix The transposed matrix of is the geological structure synergy matrix, that is, the covariance matrix between different hydrogeological structure data; is the grouting volume-geological structure submatrix, that is, the covariance matrix between the grouting volume of different grouting boreholes and the hydrogeological structure data; is the grouting amount coordination sub-matrix, that is, the covariance matrix between the grouting amounts of different grouting boreholes; It should be noted that the covariance calculated from the covariance matrix reflects the correlation between any two dimensions, namely, the correlation between any hydrogeological structure data and the grouting volume of the grouting borehole, the correlation between different hydrogeological structure data, and the relationship between different grouting volumes of the grouting borehole. The magnitude and sign of the value reflect the correlation between the two dimensions, including positive correlation, negative correlation, or no correlation. A larger covariance value indicates a more significant positive correlation between the two dimensions; conversely, a smaller covariance value indicates a more significant negative correlation between the two dimensions; and a covariance value closer to 0 indicates a weaker correlation between the two dimensions.

[0035] It should be further explained that the quantitative correlation between hydrogeological structure data and grouting volume of grouting boreholes reflects the contribution of the grouting volume of each grouting borehole to the improvement of the overall hydrogeological structure; the quantitative correlation between different hydrogeological structure data reflects whether the overall hydrogeological structure is improved; the quantitative relationship between the grouting volumes of different grouting boreholes reflects whether there is a positive or negative mutual influence between the grouting volumes of these two grouting boreholes.

[0036] Furthermore, through the formula Calculate the overall hydrogeological synergy index , each grouting drilling hole Contribution index and grouting coordination .

[0037] Step 2: Grouting project operation assessment; See also Figure 2 As shown (for example), during the grouting treatment process at the water outlet, the water level elevation in each monitoring well in the fourth aquifer is monitored in real time to generate a water level elevation curve for the fourth aquifer. Feature extraction is performed on the water level elevation curve for the fourth aquifer, and the average water level elevation in each monitoring well is calculated to obtain the average water level curve.

[0038] Get the average thickness of the fourth aquifer , in m, to obtain the time when the loss occurs Average water level elevation in the average water level curve , in m; get the time when the loss ends Average water level elevation in the average water level curve , in m; get the preset lead time before the loss occurs When the water level in the average water level curve is , unit is m; By formula Calculation of groundwater loss parameters , Governance efficiency parameters and governance effect parameters ,in 、 and All are preset influence coefficients; See also Figure 3 As shown (for example), during the grouting treatment process at the water surface, the total grouting flow of all grouting holes is monitored in real time. and the average pore pressure of all grouting holes , generate grouting flow and pore pressure duration curves.

[0039] Extract features from the grouting flow rate and pore pressure duration curve to obtain the analytical formula for the grouting flow rate ; Get the maximum average pore pressure Time of maximum average pore pressure generation , through the formula ; Calculating material contribution parameters and pore pressure control parameters ,in is the preset influence coefficient, is the preset pore pressure threshold; is the cumulative grouting volume, in units of The unit of the maximum average pore pressure is Pa, and the unit of the maximum average pore pressure generation time is seconds; the material contribution parameter The unit is .

[0040] Step 3: Economic benefit evaluation; Get each grouting borehole Total footage , one-step progress , Second opening footage and three-opening footage ;The unit is m; Get each grouting borehole The final cumulative injection volume of sealing materials, including cement injection volume , unit is t; fine aggregate injection amount , unit is t; and water glass injection volume , unit is t; get each grouting borehole Grouting loss , the unit is ; By formula Calculation of economic evaluation parameters ,in 、 and is the preset weight factor.

[0041] Step 4: Environmental impact assessment; Obtain data on chemical seepage from surface water treatment, including the quaternary aquifer value and value , through the formula Calculation of environmental impact parameters ;in and The unit is ; is the cumulative grouting volume, in units of .

[0042] Step 5: Comprehensive evaluation and report generation; Record all evaluation data obtained through calculation, including hydrogeological synergy index , each grouting drilling hole Contribution index , grouting coordination , groundwater loss parameters , Governance efficiency parameters , governance effect parameters , material contribution parameters , pore pressure control parameters , economic evaluation parameters and economic evaluation parameters Generate an evaluation report and highlight the evaluation data that exceeds the threshold.

[0043] The preferred embodiments of the present invention disclosed above are intended only to help illustrate the present invention. These preferred embodiments do not exhaustively describe all details, nor do they limit the present invention to specific embodiments. Obviously, many modifications and variations are possible based on the contents of this specification. These embodiments are selected and described in detail in this specification to better explain the principles and practical applications of the present invention, thereby enabling those skilled in the art to better understand and utilize the present invention. The present invention is limited only by the claims and their full scope and equivalents.

Claims

1. A method for evaluating the effect of water discharge treatment on a working face, characterized in that: The following steps are involved: Step 1: Hydrological data monitoring and analysis: Hydrological data includes working face hydrogeological structure data and working face natural hydrological condition data. Working face hydrogeological structure data includes mining area aquifer water level, borehole unit water inflow, and water depth in the open mining area. Natural hydrological condition data includes daily precipitation, surface runoff, and groundwater level. Hydrological data are used to generate the working face hydrogeological structure matrix and working face hydrological condition matrix, and extract the hydrogeological synergy index, the contribution index of each grouting borehole, and the grouting synergy degree. Step 2: Grouting Project Operation Evaluation: Real-time monitoring of the water level elevation of each monitoring well in the fourth aquifer in the working face, generating a water level elevation duration curve, and calculating groundwater loss parameters, treatment efficiency parameters, and treatment effect parameters. Simultaneously, real-time monitoring of the total grouting flow and average pore pressure of the grouting holes, generating a grouting flow and pore pressure duration curve, and calculating the material contribution parameter and pore pressure control parameter. Step 3: Economic benefit evaluation: Obtain the total footage of each grouting borehole and the footage data at different stages, as well as the final cumulative injection volume of the sealing material and the grouting loss volume, and calculate the economic evaluation parameters; Step 4: Environmental impact assessment: Obtain the leakage data of chemical substances treated at the water surface, including the COD and BOD values ​​of the fourth aquifer, and calculate the environmental impact parameters; Step 5. Comprehensive evaluation and report generation: Record all evaluation data, including the hydrogeological synergy index, the contribution index of each grouting borehole, the grouting synergy, the groundwater loss parameter, the treatment efficiency parameter, the treatment effect parameter, the material contribution parameter, the pore pressure control parameter, the economic evaluation parameter and the environmental impact parameter, generate an evaluation report, and highlight the evaluation data that exceeds the threshold.

2. The method for evaluating the effect of water discharge treatment on a working face according to claim 1, characterized in that: The step 1 specifically includes: S101, data definition and collection: number each grouting borehole, and collect the hydrogeological structure data of the stratum where the working face is located, including the water level of the aquifer in the mining area, the unit water inflow of the aquifer borehole, the water accumulation depth of the underground empty mining area, the rock stratum fracture ratio of the working face, the maximum fracture spacing of the rock stratum in the working face, and the minimum fracture spacing of the rock stratum in the working face. At the same time, the cumulative total injection volume of each grouting borehole is calculated; the hydrogeological structure data and the cumulative total injection volume of each grouting borehole are merged to generate a hydrogeological structure data matrix and a cumulative grouting volume matrix; S102, weighted processing of hydrogeological structure data: Calculating a locally weighted hydrogeological structure data matrix based on the elevation of the working face and the ground elevation, combined with the hydrogeological structure data matrix, daily precipitation, surface runoff, and groundwater level; S103, matrix merging and center matrix calculation: merging the cumulative grouting volume matrix and the locally weighted hydrogeological structure data matrix to obtain a hydrological comprehensive matrix, and calculating the center matrix; S104, covariance matrix calculation and feature extraction: calculate the covariance matrix of the central matrix, and extract the overall hydrogeological synergy index, the contribution index of each grouting borehole, and the grouting synergy degree from the covariance matrix of the central matrix.

3. The method for evaluating the effect of water discharge treatment on a working face according to claim 2, characterized in that: The S101 specifically includes: Number each grouting borehole, the number symbol is ; n is the total number of grouting holes; the data observation period is numbered, and the number symbol is ; Where T is the number of days included in the entire data observation period; Aquifer water level in the mining area Aquifer borehole unit water yield , Depth of water accumulation in underground mining area , rock fracture rate of working face , Maximum crack spacing of rock strata at the working face Minimum crack spacing between working face rock layers , all in mm; the fracture rate of the rock formation at the working face is a dimensionless parameter; The specific process of merging the hydrogeological structure data and the cumulative total injection volume of each grouting borehole to generate a data matrix is ​​as follows: Generate hydrogeological structure data matrix; ； Count the cumulative total injection volume of each grouting borehole i on different recording dates t , unit is t, generate cumulative grouting amount matrix .

4. The working face water discharge treatment effect evaluation method according to claim 2 or 3, characterized in that: The specific process of S102 is as follows: Obtain the elevation H of the working face and the ground elevation H0; collect the natural hydrological condition data of the working face environment during each data observation period, including daily precipitation , surface runoff and groundwater level ; Among them, the elevation of the working surface and ground elevation The unit is ; The unit of daily precipitation is ; The unit of surface runoff is ; The unit of groundwater level is ; By formula Calculate daily natural hydrological condition parameters ,in and All are preset weight factors; generate hydrological condition matrix ; Through the formula Calculate the locally weighted hydrogeological structure data matrix V, where | is the augmented matrix dividing line.

5. The method for evaluating the effect of water discharge treatment on a working face according to claim 2 or 3, characterized in that: The specific process of S103 is as follows: Merged cumulative grouting volume matrix and the locally weighted hydrogeological structure data matrix , and obtain the hydrological comprehensive matrix ,in is the augmented matrix dividing line; By formula: Calculate the center matrix ,in and is the weighted hydrogeological structure data matrix and cumulative grouting volume matrix The column mean vector of .

6. The working face water discharge treatment effect evaluation method according to claim 2 or 3, characterized in that: The specific process of S104 is as follows: By formula calculate The covariance matrix of the dimension ,in: ； where · is the matrix multiplication operator; is the central matrix The transposed matrix of is the geological structure synergy matrix, that is, the covariance matrix between different hydrogeological structure data; is the grouting volume-geological structure submatrix, that is, the covariance matrix between the grouting volume of different grouting boreholes and the hydrogeological structure data; is the grouting amount coordination sub-matrix, that is, the covariance matrix between the grouting amounts of different grouting boreholes; By formula Calculate the overall hydrogeological synergy index , each grouting drilling hole Contribution index and grouting coordination .

7. The method for evaluating the effect of water discharge treatment on a working face according to any one of claims 1 to 3, characterized in that: The specific process of step 2 is as follows: During the grouting treatment process at the water outlet, the water level elevation in each monitoring well in the fourth aquifer is monitored in real time to generate a water level elevation duration curve for the fourth aquifer; feature extraction is performed on the water level elevation duration curve for the fourth aquifer, and the average water level elevation in each monitoring well is calculated to obtain an average water level curve; Get the average thickness of the fourth aquifer , in m, to obtain the time when the loss occurs Average water level elevation in the average water level curve , unit is m; Get the time when the leakage ends Average water level elevation in the average water level curve , in m; get the preset lead time before the loss occurs When the water level in the average water level curve is , unit is m; By formula Calculation of groundwater loss parameters , Governance efficiency parameters and governance effect parameters ,in 、 and All are preset influence coefficients; During the grouting treatment process at the water surface, the total grouting flow of all grouting holes is monitored in real time , the unit is and the average pore pressure of all grouting holes , the unit is MPa, generate the grouting flow rate and pore pressure duration curve; Extract features from the grouting flow rate and pore pressure duration curve to obtain the analytical formula for the grouting flow rate ; Get the maximum average pore pressure Time of maximum average pore pressure generation , through the formula ; Calculating material contribution parameters and pore pressure control parameters ,in is the preset influence coefficient, is the preset pore pressure threshold; is the cumulative grouting volume, in units of The unit of the maximum average pore pressure is Pa, and the unit of the maximum average pore pressure generation time is seconds; the material contribution parameter The unit is .

8. The method for evaluating the effect of water discharge treatment on a working face according to any one of claims 1 to 3, characterized in that: The specific process of calculating the economic evaluation parameters in step 3 is as follows: Get each grouting borehole Total footage , one-step progress , Second opening footage and three-opening footage , the unit is m; Get each grouting borehole The final cumulative injection volume of sealing materials, including cement injection volume , unit is t; amount of fine aggregate injected , unit is t; and water glass injection volume , unit is t; get each grouting borehole Grouting loss , the unit is ; By formula Calculation of economic evaluation parameters ,in 、 and is the preset weight factor.

9. The method for evaluating the effect of water discharge treatment on a working face according to any one of claims 1 to 3, characterized in that: The specific process of calculating the environmental impact parameters in step 4 is as follows: Obtain data on chemical seepage from surface water treatment, including the quaternary aquifer value and value , through the formula Calculation of environmental impact parameters ;in and The unit is ; is the cumulative grouting volume, in units of .

Images