Mine water disaster grouting treatment effect evaluation method and system

Abstract

The present application relates to the technical field of mine safety evaluation, in particular to a mine water disaster grouting treatment effect evaluation method and system.The present application can accurately identify and delineate the most active water level change area responding to grouting by quantitatively extracting multiple hydrological dynamic characteristics such as water level slope change rate, peak rebound interval and continuous decline time length, and then accurately spatially superimposes and compares the area with the actual action area of grouting agent determined by analyzing the dynamic relationship between grouting flow rate and permeability, thereby directly verifying whether the plugging material accurately covers the key water conducting channel, solving the problem that the prior art cannot determine the pertinence of the treatment measures only by the change of macroscopic parameters, and at the same time, the continuous time stability of residual water pressure in the plugging area is analyzed to determine the long-term reliability of the plugging body, which makes up for the defects of traditional methods that only pay attention to the change of water inflow and ignore the quality and durability of the plugging body.

Description

Technical Field

[0001] This invention relates to the field of mine safety assessment technology, and in particular to a method and system for evaluating the effectiveness of grouting treatment for mine water hazards. Background Technology

[0002] The field of mine safety assessment technology involves the identification, monitoring, assessment, and pre-control of various safety risk factors during mine operation, including mine gas control, mine water hazard prevention, surrounding rock stability analysis, mine pressure manifestation monitoring, and disaster early warning.

[0003] Among them, the evaluation method for the treatment effect of grouting in mine water hazards refers to the means of analyzing the treatment effect by observing hydrogeological parameters such as water level changes, water inflow changes and grouting material diffusion range during the underground mining process of coal mines or metal mines to prevent or control water inrush disasters.

[0004] Current technologies for evaluating the effectiveness of grouting treatment for mine water hazards primarily rely on direct observation of macroscopic hydrogeological parameters such as water inflow and water level. This approach has significant limitations. Due to the lack of in-depth analysis of the dynamic changes in hydrogeological parameters, the evaluation process cannot accurately pinpoint the core area where the hydraulic response caused by grouting is most intense. When mine hydrogeological conditions are complex, simply observing a reduction in overall water inflow or a drop in water level in a certain area may erroneously interpret a large-scale, weak hydraulic impact as a successful local treatment, thus neglecting the true assessment of the sealing effect on key water outlets or main water diversion channels, leading to potential water inrush risks. Furthermore, existing technologies typically focus only on the final results after treatment, such as the diffusion range of the grouting material. However, they lack effective process monitoring and analysis methods to assess whether the diffusion behavior of the grouting agent during injection is abnormal, such as whether there is significant loss or failure to bond as expected. At the same time, evaluation methods emphasize short-term effects and lack reliable criteria for judging whether the formed sealing body possesses long-term stability and compressive strength. Ultimately, the evaluation conclusions tend to be qualitative descriptions with strong subjectivity, making it difficult to conduct refined and standardized comparisons of the advantages and disadvantages of different grouting processes. This hinders the continuous improvement and optimization of treatment technologies. Summary of the Invention

[0005] To address the technical problems existing in the prior art, embodiments of the present invention provide a method for evaluating the effectiveness of grouting treatment for mine water hazards, comprising the following steps:

[0006] To achieve the above objectives, the present invention adopts the following technical solution: a method for evaluating the effectiveness of grouting treatment for mine water hazards, comprising the following steps:

[0007] S1: Obtain water level change information in the corresponding grouting area before and after each grouting method, extract water level change features, identify the grouting area with the most significant water level change and draw the corresponding boundary curve, compare the maximum offset amplitude and density of the boundary curve, and determine the water level change area.

[0008] S2: Obtain the residual water pressure in the sealing area corresponding to the grouting area, analyze the continuity of water pressure changes over a continuous time period to determine the sealing status and obtain the sealing effect judgment result;

[0009] S3: Obtain the grouting flow rate and permeability information of the boundary of the stable grouting area in the sealing effect judgment result, analyze the diffusion behavior of the grouting agent outside the boundary of the target area, and obtain the grouting agent diffusion analysis result;

[0010] S4: Based on the grouting agent diffusion analysis results, determine whether the grouting agent is distributed around the water level change area, identify the grouting agent's action area, and obtain effective sealing coverage information;

[0011] S5: Obtain the regional boundary displacement data of the water level change area under each grouting method, and combine it with the grouting agent action area of ​​the effective sealing coverage information to generate the grouting treatment effect evaluation result.

[0012] As a further aspect of the present invention, step S1 specifically comprises:

[0013] S101: Obtain water level change information in the corresponding grouting area before and after each grouting method, extract water level change features from them, collect and integrate them, and establish a set of water level feature parameters for measuring points.

[0014] S102: Compare the characteristic parameters of each measuring point in the set of water level characteristic parameters of the measuring points with the corresponding benchmark values ​​one by one, filter the set of measuring points whose parameter values ​​exceed the benchmark values, connect the outermost measuring points of the target measuring point set and draw the boundary curve of the active water level area.

[0015] S103: Calculate the maximum offset amplitude and density corresponding to each measuring point on the boundary curve of the active water level area, and compare them with the offset amplitude threshold and density threshold to determine the water level change area.

[0016] As a further aspect of the present invention, the water level change characteristics include the rate of change of water level slope, the peak rebound interval, and the duration of continuous decline.

[0017] As a further aspect of the present invention, step S2 specifically comprises:

[0018] S201: Obtain the residual water pressure in the sealing area corresponding to the grouting area, analyze the trend of water pressure change over time, and calculate the water pressure change amplitude at each monitoring point;

[0019] S202: Determine whether the water pressure change amplitude at each monitoring point is continuously maintained within the water pressure stability discrimination interval and reaches the shortest maintenance period, determine the continuity of water pressure change, and generate the water pressure stabilization period of the blocked area.

[0020] S203: Substitute the water pressure stabilization cycle of the blocked area into the preset blocking state classification rules, classify and determine the blocking state of the current blocked area, and generate a blocking effect judgment result.

[0021] As a further aspect of the present invention, the preset blocking state division rule is determined based on the length of the water pressure stabilization period in the blocking area.

[0022] As a further aspect of the present invention, step S3 specifically comprises:

[0023] S301: Obtain the grouting flow rate and permeability information of the boundary of the stable grouting area in the sealing effect judgment result, and establish a boundary diffusion correlation parameter set;

[0024] S302: Analyze the variation trend of the grouting flow velocity of the boundary diffusion correlation parameter over a continuous time, and combine the variation trend with the permeability at the same time point to determine the synchronous variation law of the two, and obtain the analysis results of the dynamic relationship between flow velocity and permeability.

[0025] S303: Call the flow velocity-permeability dynamic relationship analysis results. When the flow velocity is continuously decreasing and the permeability is not increasing, it is determined that the diffusion behavior of the grouting agent is abnormal, and the grouting agent diffusion analysis results are generated.

[0026] As a further aspect of the present invention, step S4 specifically comprises:

[0027] S401: Obtain the grouting agent action area from the grouting agent diffusion analysis results, and call the water level change area. Align and overlay the geospatial data of these two areas to generate a spatial overlay map of the grouting range and water level change area.

[0028] S402: Based on the spatial overlay diagram of the grouting range and the water level change zone, calculate the spatial distance between the outer boundary of the grouting agent action area and the outer boundary of the water level change zone, determine whether the grouting agent is distributed outside the water level change zone, and obtain the spatial deviation between the grouting boundary and the water level boundary.

[0029] S403: Based on the spatial deviation between the grouting boundary and the water level boundary, the calculated deviation is compared with the allowable deviation value for effective sealing coverage to identify whether the grouting agent's action area completely covers the water level change area, thereby obtaining the effective sealing coverage information.

[0030] As a further aspect of the present invention, step S5 specifically comprises:

[0031] S501: Obtain the regional boundary displacement data of the water level change area under each grouting method, combine it with the grouting agent action area in the effective coverage information of the sealing, and integrate it with the sealing effect judgment result to establish a sub-method governance efficiency dataset;

[0032] S502: Extract the regional boundary displacement data and grouting agent action range of each grouting method in the sub-method governance efficiency dataset, calculate the maximum offset amplitude of the boundary displacement and the degree of overlap of the action range, and use it as the boundary dynamic control conformity.

[0033] S503: Invoke the boundary dynamic control conformity and combine it with the sub-method governance efficiency dataset to comprehensively evaluate the boundary control, sealing status and diffusion behavior of each grouting method, and generate grouting governance effect evaluation results to assist in selecting the most suitable grouting method to achieve the best sealing effect and water level change control.

[0034] The evaluation system for the effectiveness of grouting treatment in mine water hazards includes:

[0035] The boundary recognition module obtains water level change information in the corresponding grouting area before and after each grouting method, extracts water level change features, identifies the grouting area with the most significant water level change and draws the corresponding boundary curve, compares the maximum offset amplitude and density of the boundary curve, and determines the water level change area.

[0036] The sealing effect judgment module obtains the residual water pressure in the sealing area corresponding to the grouting area, analyzes the continuity of water pressure changes over a continuous time, determines the sealing status, and obtains the sealing effect judgment result.

[0037] The diffusion behavior analysis module obtains the grouting flow rate and permeability information of the boundary of the stable grouting area in the sealing effect judgment result, analyzes the diffusion behavior of the grouting agent outside the boundary of the target area, and obtains the grouting agent diffusion analysis result.

[0038] The sealing and coverage assessment module, based on the grouting agent diffusion analysis results, determines whether the grouting agent is distributed outside the water level change area, identifies the grouting agent's effective area, and obtains sealing and coverage information.

[0039] The treatment effect evaluation module acquires the regional boundary displacement data of the water level change area under each grouting method, and generates the grouting treatment effect evaluation result by combining the grouting agent action area with the effective coverage information of the sealing.

[0040] Compared with the prior art, the advantages and positive effects of the present invention are as follows:

[0041] This invention, by quantitatively extracting multiple hydrological dynamic features such as the rate of change of water level slope, peak rebound interval, and duration of continuous decline, can accurately identify and delineate the water level change area most actively responding to grouting. This area is then precisely spatially superimposed and compared with the actual action area of ​​the grouting agent determined by analyzing the dynamic relationship between grouting flow rate and permeability. This directly verifies whether the sealing material accurately covers the key water-conducting channels, solving the problem that existing technologies cannot determine the targeted nature of treatment measures based solely on changes in macroscopic parameters. Furthermore, by analyzing the continuous-time stability of residual water pressure within the sealing area, the long-term stability of the sealing body can be assessed. The reliability of the grouting process is assessed, which overcomes the shortcomings of traditional methods that only focus on changes in water inflow while ignoring the quality and durability of the sealing body. In addition, the analysis of the synchronous change law of grouting flow rate and permeability can identify abnormal diffusion behavior during the grouting process, providing a new dimension for evaluating the effectiveness of the treatment process. Finally, the boundary control compliance, sealing stability and diffusion behavior are comprehensively evaluated to form a multi-dimensional, full-process quantitative evaluation system. This provides a scientific basis for objectively comparing the advantages and disadvantages of different grouting methods and optimizing treatment schemes, significantly improving the accuracy and comprehensiveness of the evaluation of the grouting treatment effect of mine water hazards. Attached Figure Description

[0042] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0043] Figure 1 This is a schematic diagram of the steps of the present invention;

[0044] Figure 2 This is a detailed schematic diagram of S1 of the present invention;

[0045] Figure 3 This is a detailed schematic diagram of S2 of the present invention;

[0046] Figure 4 This is a detailed schematic diagram of S3 of the present invention;

[0047] Figure 5 This is a detailed schematic diagram of S4 of the present invention;

[0048] Figure 6 This is a detailed schematic diagram of S5 of the present invention;

[0049] Figure 7 This is a system module diagram of the present invention. Detailed Implementation

[0050] The technical solution of the present invention will now be described with reference to the accompanying drawings.

[0051] In embodiments of the present invention, words such as "exemplarily," "for example," etc., are used to indicate that something is an example, illustration, or description. Any embodiment or design described as "exemplary" in the present invention should not be construed as being more preferred or advantageous than other embodiments or designs. Specifically, the use of the word "exemplary" is intended to present the concept in a concrete manner. Furthermore, in embodiments of the present invention, the meaning expressed by "and / or" can be both, or either one.

[0052] In the embodiments of this invention, the terms "image" and "picture" may sometimes be used interchangeably. It should be noted that, without emphasizing the distinction between them, they convey the same meaning. Similarly, the terms "of," "corresponding (relevant)," and "corresponding" may sometimes be used interchangeably. It should be noted that, without emphasizing the distinction between them, they convey the same meaning.

[0053] In this embodiment of the invention, sometimes a subscript such as W1 may be written in a non-subscript form such as W1. When the difference is not emphasized, the meaning they express is the same.

[0054] To make the technical problems, technical solutions and advantages of the present invention clearer, a detailed description will be given below in conjunction with the accompanying drawings and specific embodiments.

[0055] Please see Figure 1 This invention provides a method for evaluating the effectiveness of grouting treatment for mine water hazards, comprising the following steps:

[0056] S1: Obtain water level change information in the corresponding grouting area before and after each grouting method, extract water level change features, identify the grouting area with the most significant water level change and draw the corresponding boundary curve, compare the maximum offset amplitude and density of the boundary curve, and determine the water level change area.

[0057] S2: Obtain the residual water pressure in the sealing area corresponding to the grouting area, analyze the continuity of water pressure changes over a continuous time period to determine the sealing status and obtain the sealing effect judgment result;

[0058] S3: Obtain the grouting flow rate and permeability information of the boundary of the stable grouting area in the sealing effect judgment result, analyze the diffusion behavior of the grouting agent outside the boundary of the target area, and obtain the grouting agent diffusion analysis results;

[0059] S4: Based on the grouting agent diffusion analysis results, determine whether the grouting agent is distributed outside the water level change area, identify the grouting agent's action area, and obtain information on the effective coverage of the sealing.

[0060] S5: Obtain the regional boundary displacement data of the water level change area under each grouting method, and combine it with the grouting agent action area of ​​the effective sealing coverage information to generate the grouting treatment effect evaluation result.

[0061] Please see Figure 2 Step S1 is as follows:

[0062] S101: Obtain water level change information in the corresponding grouting area before and after each grouting method, extract water level change features from them, collect and integrate them, and establish a set of water level feature parameters for measuring points. Water level change features include water level slope change rate, peak rebound interval and continuous decline time length.

[0063] First, data on the continuous changes in water level elevation over time before and after grouting were obtained from multiple water level monitoring points pre-deployed in the rock mass surrounding the underground tunnel excavation face. Based on this data, the water level change characteristics of each monitoring point were extracted and calculated. The analysis here aims to reveal the sealing effect and influence range of the grout on the water channels in the rock mass fissures. Specifically, the process of obtaining the water level slope change rate involves calculating the rate of rise and fall of the water level per unit hour, and then calculating the change in these rates between adjacent hours. Drastic changes in this value reflect the dynamic process of the groundwater flow path being rapidly cut off or diverted by the grout. The process of obtaining the peak rebound interval is as follows... First, determine the time when the water level drops to its lowest point due to grouting and drainage. Then, record the time interval from that lowest point when the pressurized water head first rebounds due to surrounding replenishment. This interval reflects the hydraulic sealing of the blocked area. The process of obtaining the continuous drop time is to record the cumulative total time from the start of the grouting effect, during which the water level drops continuously without any rebound or leveling off. This time directly reflects the effective time of a single grouting operation. The characteristic values ​​obtained from the above three calculations are combined to establish a set of unique characteristic parameters for each monitoring point. Finally, the parameters of all monitoring points are integrated to establish a set of water level characteristic parameters for the monitoring points.

[0064] S102: Compare the characteristic parameters of each measuring point in the water level characteristic parameter set with the corresponding benchmark value one by one, filter the set of measuring points whose parameter values ​​exceed the benchmark value, connect the outermost measuring points of the target measuring point set and draw the boundary curve of the active water level area.

[0065] The system retrieves the water level characteristic parameter set of each measuring point and compares each parameter with a preset benchmark value. These benchmark values ​​are determined based on statistical analysis of data from dozens of curtain grouting projects conducted under similar aquifer geological conditions. The 85th percentile of the distribution of each characteristic parameter in the historical data is selected as the benchmark. This aims to screen out the monitoring points that are most sensitive and drastic to the grouting response. Specifically, the benchmark value for the rate of change of water level slope is set to 0.25 meters per hour, the benchmark value for the peak rebound interval is 2.5 hours, and the benchmark value for the duration of continuous decline is 4.5 hours. During the comparison, if any characteristic parameter value of a measuring point exceeds the set benchmark value range, such as if its calculated duration of continuous decline exceeds 4.5 hours, then the measuring point is screened out. All the screened measuring points together constitute the target measuring point set. Finally, the geospatial coordinates of these target measuring points are obtained, and the outermost measuring points are connected by lines on the drawing to create a closed curve, which is the boundary curve of the active water level area.

[0066] S103: Calculate the maximum offset amplitude and density corresponding to each measuring point on the boundary curve of the active water level area, compare it with the offset amplitude threshold and density threshold, and determine the water level change area.

[0067] For the drawn boundary curve of the active water level area, two key indicators need to be calculated: maximum offset amplitude and density. The maximum offset amplitude is calculated by overlaying and comparing the actual boundary curve with the theoretical influence boundary line planned in the engineering design stage on a map, and finding the maximum vertical distance between the two lines. The density is calculated by counting the total number of target measuring points falling within the area enclosed by the boundary curve and then dividing it by the total area of ​​the area. Subsequently, these two calculated indicators are compared with the corresponding thresholds. The offset amplitude threshold and density threshold are set with reference to ten similar foundation pit dewatering or tunnel water blocking projects that have been evaluated as successful. 1.2 times the statistical average value of the corresponding parameters in these successful cases is taken as the judgment standard to ensure the rigor of the evaluation. After specific calculation, the offset amplitude threshold is set to 1.8 meters and the density threshold is 0.018 measuring points per square meter. When the calculated maximum offset amplitude and density both exceed their respective thresholds, the area is finally determined to be an area with significant hydraulic changes due to grouting, that is, a water level change area.

[0068] Please see Figure 3 Step S2 is as follows:

[0069] S201: Obtain the residual water pressure in the sealing area corresponding to the grouting area, analyze the trend of water pressure change over time, and calculate the water pressure change amplitude at each monitoring point;

[0070] After the grout solidifies and forms a complete underground sealing curtain, residual water pressure readings are continuously obtained from water pressure monitoring points located on the back side (water-repellent side) of the curtain. Subsequently, the trend of water pressure values ​​recorded at each monitoring point over time is analyzed. The purpose of this analysis is to verify whether the sealing curtain has formed a long-term effective hydraulic barrier. The specific analysis process involves plotting the trend of water pressure changes over time at each monitoring point, observing whether it changes from violent fluctuations to a gradual change after grouting is completed, and calculating the water pressure change amplitude within a fixed monitoring period. The calculation method for this amplitude is to find the highest and lowest water pressure values ​​recorded in a complete monitoring period, such as 24 hours, and the difference between the two is the water pressure change amplitude for that period.

[0071] S202: Determine whether the water pressure change amplitude at each monitoring point is continuously maintained within the water pressure stability discrimination interval and reaches the shortest maintenance period, determine the continuity of water pressure change, and generate the water pressure stabilization cycle of the blocked area.

[0072] For each monitoring point, the calculated water pressure variation amplitude is used to determine whether it has reached a stable state. This determination includes two criteria: the water pressure stability discrimination interval and the shortest maintenance period. The water pressure stability discrimination interval is set based on the results of indoor pressure fatigue tests on the grouting materials (such as cement-water glass). The tests show that when the water pressure fluctuation amplitude of the solidified body is less than 0.015 MPa, its internal microcracks will not expand, and it can maintain long-term structural stability and impermeability. Therefore, the stability interval is set to 0 to 0.015 MPa. The shortest maintenance period is set in accordance with the safety observation period requirements for adjacent processes (such as tunnel secondary lining) in the building safety code of the project location, and is set to 48 hours. The determination process is to check whether the water pressure variation amplitude of a monitoring point falls within the stability discrimination interval for multiple consecutive cycles, and whether the duration of this state reaches the minimum maintenance period of 48 hours. If the conditions are met, the water pressure variation at that point is determined to be continuous, and the total time to reach stability is recorded as the water pressure stabilization period at that point.

[0073] S203: Substitute the water pressure stabilization cycle of the blocked area into the preset blocking state classification rules, classify and determine the blocking state of the current blocked area, and generate the blocking effect judgment result. The preset blocking state classification rules are determined according to the length of the water pressure stabilization cycle of the blocked area.

[0074] The water pressure stabilization period of the entire blocked area (usually the slowest response, i.e., the shortest stabilization period, among all monitoring points) is substituted into a set of preset blocking status classification rules. These rules are based on a retrospective analysis of the effects of twenty past groundwater blocking projects of different risk levels. The average stabilization time of historically rated "excellent" projects and the average stabilization time of successful projects rated "qualified" are used as the dividing line. Specifically, the rules are as follows: if the stabilization period exceeds 72 hours, the blocking status is classified as "highly efficient and stable"; if the period is between 48 and 72 hours, it is classified as "basically stable"; if the period is less than 48 hours, it is classified as "underly stable". By comparing the calculated stabilization period with these rules, the blocking status of the current blocked area is determined, and the final blocking effect judgment result is generated.

[0075] Please see Figure 4 Step S3 is as follows:

[0076] S301: Obtain the grouting flow rate and permeability information of the boundary of the stable grouting area in the sealing effect judgment result, and establish a boundary diffusion correlation parameter set;

[0077] Based on the sealing effect assessment, grouting areas classified as "highly efficient and stable" or "basically stable" were selected. Grouting flow velocity and soil permeability information measured through in-situ water injection tests were retrieved from the boundaries of these areas during the grouting process. These data, presented in pairs as time points, flow velocity values, and permeability values, were collected and organized to establish a set of related parameters for analyzing boundary diffusion behavior.

[0078] S302: Analyze the variation trend of the concentrated grouting flow velocity, a boundary diffusion correlation parameter, over a continuous time period, and combine the variation trend with the permeability at the same time point to determine the synchronous variation law of the two, and obtain the analysis results of the dynamic relationship between flow velocity and permeability.

[0079] By calling the boundary diffusion correlation parameter set, the trend of grouting flow rate over a continuous time period is first analyzed. This analysis aims to determine whether the grout is effectively filling the voids in the soil and rock. This is usually manifested as a decrease in flow rate over time under constant pressure. Then, the flow rate at each time point is combined with the permeability at the same time to determine whether there is a synchronous change pattern between the two. For example, if the grouting flow rate decreases while the permeability at the corresponding location also decreases synchronously within a certain time period, this relationship is recorded as positive synchronization. This indicates that the grout is solidifying at the expected location and reducing its water conductivity. Conversely, if the flow rate decreases while the permeability remains unchanged or increases, it is recorded as asynchronous or negative synchronization. By making such point-by-point judgments on the data of the entire grouting process, a complete dynamic analysis result of the flow rate-permeability relationship is finally obtained.

[0080] S303: Call the dynamic relationship analysis results of flow velocity and permeability. When it shows that the flow velocity is continuously decreasing and the permeability is not increasing, it is determined that the diffusion behavior of the grouting agent is abnormal and the diffusion analysis results of the grouting agent are generated.

[0081] The dynamic relationship analysis results of flow velocity and permeability are retrieved and reviewed according to a preset anomaly judgment rule. The rule is as follows: if the analysis results show a continuous phenomenon reflecting a decrease in flow velocity but no decrease in permeability, and the duration of this phenomenon exceeds the initial setting time of the grouting material (for example, the initial setting time of cement grout is set to 15 minutes), it is judged as an anomaly. This rule is set to exclude instantaneous data changes caused by normal fluctuations in pumping pressure, and to ensure that the captured data is the actual diffusion problem of grout that may be lost or not set as expected. During execution, if the dynamic relationship analysis results of a certain area are found to meet the above anomaly rule, it is judged that the diffusion behavior of the grouting agent at that location is abnormal, and this conclusion is recorded, and finally the grouting agent diffusion analysis results are generated.

[0082] Please see Figure 5 Step S4 is as follows:

[0083] S401: Obtain the grouting agent action area from the grouting agent diffusion analysis results, and call the water level change area. Align and overlay the geospatial data of these two areas to generate a spatial overlay map of the grouting range and water level change area.

[0084] The area determined to have normal diffusion behavior in the grouting agent diffusion analysis results is the actual effective area of ​​the grouting agent. At the same time, the water level change area determined in the previous steps is called, and the geospatial data of the two areas are accurately aligned in coordinate and scaled on the same map view. Then, an overlay analysis is performed to generate a spatial overlay map that can intuitively show the relative position and overlap of the two areas.

[0085] S402: Based on the spatial overlay diagram of the grouting range and the water level change zone, calculate the spatial distance between the outer boundary of the grouting agent's action area and the outer boundary of the water level change zone, determine whether the grouting agent is distributed outside the water level change zone, and obtain the spatial deviation between the grouting boundary and the water level boundary.

[0086] On the spatial overlay map, the spatial distance between the outer boundary of the grouting agent's action area and the outer boundary of the water level change area is calculated. Specifically, a series of dense sampling points are selected from the boundary of the grouting agent's action area, and then the shortest distance from each sampling point to the boundary of the water level change area is calculated. Through this process, it can be determined whether the distribution range of the grouting agent exceeds the periphery of the water level change area, or whether it fails to completely cover certain areas. All measured distance values ​​that indicate that the grouting range is smaller than the water level change range are collected to form the spatial deviation between the grouting boundary and the water level boundary.

[0087] S403: Based on the spatial deviation between the grouting boundary and the water level boundary, compare the calculated deviation with the allowable deviation value of the effective coverage of the sealing, identify whether the grouting agent action area completely covers the water level change area, and obtain the effective coverage information of the sealing.

[0088] The maximum value of the spatial deviation is compared with a preset allowable deviation value for effective sealing coverage. This allowable value is set based on the risk assessment of groundwater bypass and related engineering safety redundancy requirements, and is combined with multiple computer numerical simulations. The simulation results show that when the coverage deviation of the grouting range from the target area is less than 0.5 meters, the probability of subsequent bypass seepage leading to project failure is less than one in ten thousand. Therefore, the allowable deviation value is set at 0.5 meters. If the calculated maximum deviation exceeds this allowable value, it is identified that the grouting agent's action area has not completely covered the water level change area, and information on the effective sealing coverage is generated accordingly.

[0089] Please see Figure 6 Step S5 is as follows:

[0090] S501: Obtain the regional boundary displacement data of the water level change area under each grouting method, combine it with the grouting agent action area in the effective coverage information of the sealing, and integrate it with the sealing effect judgment results to establish a sub-method governance efficiency dataset;

[0091] For each grouting method used in the project, such as high-pressure fracturing grouting and low-pressure permeation grouting, the boundary displacement data of the corresponding water level change area, the effective coverage information of the sealing, and the final sealing effect judgment results are collected. These data obtained from different steps are correlated and integrated according to the grouting method to construct a dataset that can comprehensively reflect the treatment effectiveness of different methods.

[0092] S502: Extract the regional boundary displacement data and grouting agent action range of each grouting method in the sub-method governance efficiency dataset, calculate the maximum offset amplitude of the boundary displacement and the degree of overlap of the action range, and use it as the boundary dynamic control conformity.

[0093] From the grouting effectiveness dataset, the boundary displacement data (i.e., maximum offset amplitude) and the grouting agent's effective range (overlap calculated through spatial overlay diagrams) corresponding to each grouting method are extracted. These two indicators are then combined to calculate the boundary dynamic control conformity. The evaluation rules for this conformity are as follows: when the maximum offset amplitude is less than 1.0 meter and the overlap of the effective range is higher than 98%, it is rated as "high"; when the maximum offset amplitude is between 1.0 meter and 2.0 meter or the overlap is between 90% and 98%, it is rated as "medium"; when the maximum offset amplitude is greater than 2.0 meter or the overlap is less than 90%, it is rated as "low". The numerical definition of this rule is determined by conducting data regression analysis on a large number of completed excellent and qualified engineering cases to find the parameter range most closely related to the final treatment effect level.

[0094] S503: Call the boundary dynamic control fit degree and combine it with the sub-mode governance efficiency dataset to comprehensively evaluate the boundary control, sealing status and diffusion behavior of each grouting method, and generate grouting governance effect evaluation results to assist in selecting the most suitable grouting method to achieve the best sealing effect and water level change control.

[0095] The system invokes the assessed boundary dynamic control conformity and, combined with the blocking status and diffusion behavior information (indirectly reflected by blocking coverage information) from the sub-method governance effectiveness dataset, comprehensively evaluates each grouting method. The evaluation process employs a weighted scoring method, converting qualitative descriptions such as "high," "medium," "low," "highly efficient and stable," and "basically stable" into numerical values ​​(e.g., high = 3, medium = 2, low = 1), and assigning weights to different indicators. The weighting is based on multiple rounds of evaluation conducted by 20 senior engineering and technical experts in this field. A questionnaire survey was conducted, and experts ranked and scored the various indicators according to their importance in determining the success or failure of the entire grouting project. The results showed that the weight of boundary control and sealing status was 0.4, while diffusion behavior, as a process indicator, had a weight of 0.2. The scores of each indicator were multiplied by their weights and then summed to obtain a final comprehensive score. Based on the preset score range (for example, a score between 2.5 and 3.0 is excellent, between 1.5 and 2.5 is good, and between 1.0 and 1.5 is average), the final evaluation results of the grouting treatment effect for each grouting method were generated.

[0096] Please see Figure 7 The evaluation system for the effectiveness of grouting treatment in mine water hazards includes:

[0097] The boundary recognition module obtains water level change information in the corresponding grouting area before and after each grouting method, extracts water level change features, identifies the grouting area with the most significant water level change and draws the corresponding boundary curve, compares the maximum offset amplitude and density of the boundary curve, and determines the water level change area.

[0098] The sealing effect judgment module obtains the residual water pressure in the sealing area corresponding to the grouting area, analyzes the continuity of water pressure changes over a continuous time, determines the sealing status, and obtains the sealing effect judgment result.

[0099] The diffusion behavior analysis module obtains the grouting flow rate and permeability information of the boundary of the stable grouting area in the sealing effect judgment result, analyzes the diffusion behavior of the grouting agent outside the boundary of the target area, and obtains the diffusion analysis results of the grouting agent.

[0100] The sealing and coverage assessment module, based on the grout diffusion analysis results, determines whether the grout is distributed outside the water level change area, identifies the grout's effective coverage area, and obtains sealing and coverage information.

[0101] The treatment effect evaluation module acquires the regional boundary displacement data of the water level change area under each grouting method, and combines it with the grouting agent action area information of the effective sealing coverage information to generate the grouting treatment effect evaluation result.

[0102] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A method for evaluating the effectiveness of grouting treatment for mine water hazards, characterized in that, Includes the following steps: S1: Obtain water level change information in the corresponding grouting area before and after each grouting method, extract water level change features, identify the grouting area with the most significant water level change and draw the corresponding boundary curve, compare the maximum offset amplitude and density of the boundary curve, and determine the water level change area. S101: Obtain water level change information in the corresponding grouting area before and after each grouting method, extract water level change features from them, collect and integrate them, and establish a set of water level feature parameters for measuring points. S102: Compare the characteristic parameters of each measuring point in the set of water level characteristic parameters of the measuring points with the corresponding benchmark values ​​one by one, filter the set of measuring points whose parameter values ​​exceed the benchmark values, connect the outermost measuring points of the target measuring point set and draw the boundary curve of the active water level area. S103: Calculate the maximum offset amplitude and density corresponding to each measuring point on the boundary curve of the active water level area, and compare them with the offset amplitude threshold and density threshold to determine the water level change area; the water level change characteristics include the water level slope change rate, peak rebound interval and continuous decline time length; S2: Obtain the residual water pressure in the sealing area corresponding to the grouting area, analyze the continuity of water pressure changes over a continuous time period to determine the sealing status and obtain the sealing effect judgment result; S3: Obtain the grouting flow rate and permeability information of the boundary of the stable grouting area in the sealing effect judgment result, analyze the diffusion behavior of the grouting agent outside the boundary of the target area, and obtain the grouting agent diffusion analysis result; S301: Obtain the grouting flow rate and permeability information of the boundary of the stable grouting area in the sealing effect judgment result, and establish a boundary diffusion correlation parameter set; S302: Analyze the variation trend of the grouting flow velocity of the boundary diffusion correlation parameter over a continuous time, and combine the variation trend with the permeability at the same time point to determine the synchronous variation law of the two, and obtain the analysis results of the dynamic relationship between flow velocity and permeability. S303: Call the flow velocity-permeability dynamic relationship analysis results. When the flow velocity is continuously decreasing and the permeability is not increasing, it is determined that the diffusion behavior of the grouting agent is abnormal, and the grouting agent diffusion analysis results are generated. S4: Based on the grouting agent diffusion analysis results, determine whether the grouting agent is distributed around the water level change area, identify the grouting agent's action area, and obtain effective sealing coverage information; S5: Obtain the regional boundary displacement data of the water level change area under each grouting method, and combine it with the grouting agent action area of ​​the effective sealing coverage information to generate the grouting treatment effect evaluation result.

2. The method for evaluating the effectiveness of mine water hazard grouting treatment according to claim 1, characterized in that, Step S2 is as follows: S201: Obtain the residual water pressure in the sealing area corresponding to the grouting area, analyze the trend of water pressure change over time, and calculate the water pressure change amplitude at each monitoring point; S202: Determine whether the water pressure change amplitude at each monitoring point is continuously maintained within the water pressure stability discrimination interval and reaches the shortest maintenance period, determine the continuity of water pressure change, and generate the water pressure stabilization period of the blocked area. S203: Substitute the water pressure stabilization cycle of the blocked area into the preset blocking state classification rules, classify and determine the blocking state of the current blocked area, and generate a blocking effect judgment result.

3. The method for evaluating the effectiveness of mine water hazard grouting treatment according to claim 2, characterized in that, The preset blocking state division rules are determined based on the length of the water pressure stabilization period in the blocking area.

4. The method for evaluating the effectiveness of mine water hazard grouting treatment according to claim 1, characterized in that, Step S4 is as follows: S401: Obtain the grouting agent action area from the grouting agent diffusion analysis results, and call the water level change area. Align and overlay the geospatial data of these two areas to generate a spatial overlay map of the grouting range and water level change area. S402: Based on the spatial overlay diagram of the grouting range and the water level change zone, calculate the spatial distance between the outer boundary of the grouting agent action area and the outer boundary of the water level change zone, determine whether the grouting agent is distributed outside the water level change zone, and obtain the spatial deviation between the grouting boundary and the water level boundary. S403: Based on the spatial deviation between the grouting boundary and the water level boundary, the calculated deviation is compared with the allowable deviation value for effective sealing coverage to identify whether the grouting agent's action area completely covers the water level change area, thereby obtaining the effective sealing coverage information.

5. The method for evaluating the effectiveness of mine water hazard grouting treatment according to claim 1, characterized in that, Step S5 is as follows: S501: Obtain the regional boundary displacement data of the water level change area under each grouting method, combine it with the grouting agent action area in the effective coverage information of the sealing, and integrate it with the sealing effect judgment result to establish a sub-method governance efficiency dataset; S502: Extract the regional boundary displacement data and grouting agent action range of each grouting method in the sub-method governance efficiency dataset, calculate the maximum offset amplitude of the boundary displacement and the degree of overlap of the action range, and use it as the boundary dynamic control conformity. S503: Invoke the boundary dynamic control conformity and combine it with the sub-method governance efficiency dataset to comprehensively evaluate the boundary control, sealing status and diffusion behavior of each grouting method, and generate grouting governance effect evaluation results to assist in selecting the most suitable grouting method to achieve the best sealing effect and water level change control.

6. A system for evaluating the effectiveness of grouting treatment for mine water hazards, characterized in that, The system is used to implement the method for evaluating the effect of mine water hazard grouting treatment as described in any one of claims 1-5, and the system includes: The boundary recognition module obtains water level change information in the corresponding grouting area before and after each grouting method, extracts water level change features, identifies the grouting area with the most significant water level change and draws the corresponding boundary curve, compares the maximum offset amplitude and density of the boundary curve, and determines the water level change area. The sealing effect judgment module obtains the residual water pressure in the sealing area corresponding to the grouting area, analyzes the continuity of water pressure changes over a continuous time, determines the sealing status, and obtains the sealing effect judgment result. The diffusion behavior analysis module obtains the grouting flow rate and permeability information of the boundary of the stable grouting area in the sealing effect judgment result, analyzes the diffusion behavior of the grouting agent outside the boundary of the target area, and obtains the grouting agent diffusion analysis result. The sealing and coverage assessment module, based on the grouting agent diffusion analysis results, determines whether the grouting agent is distributed outside the water level change area, identifies the grouting agent's effective area, and obtains sealing and coverage information. The treatment effect evaluation module acquires the regional boundary displacement data of the water level change area under each grouting method, and generates the grouting treatment effect evaluation result by combining the grouting agent action area with the effective coverage information of the sealing.

Images