A water level real-time tracking observation method based on working face water outflow treatment

Abstract

The application discloses a water level real-time tracking observation method based on working face water treatment, relates to the field of water level real-time tracking observation, and comprises the following steps: judging and analyzing the particle grading state and the particle stress state of the combination of cement slurry and fissure zones after grouting treatment of working face water, obtaining the combination synoptic value, analyzing the fissure state before and after grouting of tectonic fissures, obtaining the filling degree value of tectonic fissures, performing real-time monitoring and analysis on the four aquifer levels of working face treatment, obtaining the four aquifer level estimation value, performing water level comprehensive state evaluation of working face water treatment on the combination synoptic value, the filling degree value of tectonic fissures and the four aquifer level estimation value, obtaining the water level comprehensive evaluation value, receiving the water level comprehensive evaluation value, comparing and analyzing the water level comprehensive evaluation value with the set water level comprehensive evaluation threshold, and making corresponding alarms and determinations, so that the water level comprehensive state can be analyzed from multiple angles and in deep levels, and safety accidents in coal mining can be reduced.

Description

Technical Field

[0001] This invention relates to the field of real-time water level tracking and observation, specifically a method for real-time water level tracking and observation based on the treatment of effluent from a working face. Background Technology

[0002] With the development of society and economy, people's demand for energy is increasing. Coal mining is also attracting more and more attention. During the mining process, the working face is affected by the sandstone fissure water and other water-bearing elements in the roof, so grouting and other treatment measures are required. For the treated working face, the water level needs to be monitored in real time. Traditional real-time water level monitoring often relies on personnel using simple measuring instruments, which cannot guarantee the accuracy and precision of real-time monitoring.

[0003] In the real-time monitoring of water levels during the treatment of water in the working face, the determination of the comprehensive water level status is a key factor affecting the safety of the working face. Failure to determine the comprehensive water level status in a timely manner and issue alarms for abnormal conditions can lead to safety accidents in coal mining. However, the determination of the comprehensive water level status often lacks the ability to analyze it from multiple perspectives and in depth. Therefore, timely determination of the comprehensive water level status and in-depth analysis from multiple perspectives are the main technical problems to be solved. Summary of the Invention

[0004] The purpose of this invention is to provide a real-time water level tracking and observation method based on the treatment of water effluent from the working face. This method can take into account the particle size distribution and particle stress of the formed aggregate, as well as the fracture state, and also consider the real-time monitoring and analysis of water level from four perspectives. It can analyze the comprehensive water level status from multiple angles and in depth.

[0005] The objective of this invention can be achieved through the following technical solutions:

[0006] A method for real-time water level tracking and observation based on effluent treatment at the working face includes the following steps:

[0007] Step 1: After the cement grouting treatment of the working face, the combination formed by the cement grout and the crack zone is analyzed for particle size distribution and particle stress state to obtain the composite value of the combination.

[0008] Step 2: Analyze the state of the structural cracks before and after grouting to obtain the filling degree value of the structural cracks;

[0009] Step 3: Conduct real-time monitoring and analysis of the water level of the four cut surfaces at the working face to obtain an estimated water level status.

[0010] Step 4: Evaluate the overall water level status of the working face water treatment based on the combined body morphology value, the filling degree value of structural fissures, and the estimated water level value of the four water-bearing structures, and obtain the overall water level status value through graphic construction processing;

[0011] Step 5: Receive the comprehensive water level assessment value, compare and analyze it with the set comprehensive water level assessment threshold, and issue corresponding alarms and judgments.

[0012] Optionally, the method for obtaining the complex symmetry value is as follows:

[0013] Step 101: The cement grout and the fracture zone formed by the grouting treatment of the working face after water discharge are imaged using an electron microscope. The image is converted into a grayscale image. Threshold segmentation using the maximum inter-class variance method is used to separate all particles from the image background. The pixel area of ​​each particle is calculated, and the pixel area is converted to the actual area according to a set ratio coefficient to obtain the actual area sm of each particle. The equivalent diameter dz of each particle is calculated by formula. The equivalent diameter of each particle is plotted on the x-axis with the particle diameter size and the frequency of particle occurrence in each size interval as the y-axis to establish a particle size frequency distribution histogram. The frequency number in each size interval is counted and compared with the set reference peak frequency threshold for each size interval. When the frequency number of the size interval is greater than the reference peak frequency threshold of the size interval, the size interval is determined as a peak interval. The distance between each adjacent peak interval is counted and represented by the number of intervals. The total peak distance value zj is obtained by summing the distance between the intervals and the number of peaks fs. The peak value fz of the composite is obtained by calculating the total peak distance value zj and the number of peaks fs.

[0014] Step 102: Arrange the equivalent diameters of each particle in ascending order to obtain the particle diameters in the 60% and 10% ranges. Divide the particle diameters in the 60% range by the particle diameters in the 10% range to obtain the particle distribution non-uniformity coefficient value bj. Obtain the particle diameters in the 30% range, square them, and divide them by the product of the particle diameters in the 60% and 10% ranges to obtain the particle distribution curvature coefficient value ql. Calculate the particle distribution non-uniformity coefficient value bj, the particle distribution curvature coefficient value ql, and the peak value fz of the aggregate using the formula to obtain the aggregate gradation value jz.

[0015] Step 103: According to the set formula The particle contact force JF is calculated, ZF is the normal force; lF is the centrifugal force of each particle; ft represents the tangential contact force under static and dynamic friction conditions.

[0016] Step 104: Calculate the particle contact area jm, divide the particle contact force JF by the particle contact area jm, and multiply by the correction factor coefficient to obtain the particle contact stress value. Obtain the contact stress value of each particle in the composite, process the average contact stress value of each particle to obtain the average contact stress value of each particle, take the maximum contact stress value of each particle in the composite and divide it by the average contact stress value to obtain the stress distribution value of the composite.

[0017] Step 105: Multiply the gradation value jz of the composite by the weighting factor coefficient g1 to obtain product coefficient one; multiply the stress concentration anisotropy value of the composite by the weighting factor coefficient g2 to obtain product coefficient two; divide product coefficient one by product coefficient two to obtain the composite gradation value.

[0018] Optionally, step 101 specifically includes:

[0019] The cement slurry and fracture zone bond formed after grouting treatment at the working face were imaged using a scanning electron microscope. The image was converted to grayscale, and threshold segmentation using the maximum inter-class variance method was applied to separate all particles from the image background. The pixel area of ​​each particle was calculated, and the pixel area was converted to the actual area according to a set ratio coefficient to obtain the actual area sm of each particle. The formula yields the equivalent diameter dz of each particle. A particle size frequency distribution histogram is constructed, plotting the particle diameter on the x-axis and the frequency of particle occurrence within each size interval on the y-axis. The frequency count within each size interval is counted and compared with a predefined reference peak frequency threshold. If the frequency count of a size interval exceeds its reference peak frequency threshold, that interval is identified as a peak interval. The distance between adjacent peak intervals is counted, represented by the number of intervals, and summed to obtain the total peak interval value zj. The number of peaks fs is also counted. Finally, the total peak interval value zj and the number of peaks fs are normalized according to the predefined formula. The peak value fz of the composite is obtained by calculation, where a1 and a2 are natural constants, and a1>a2, and h1 and h2 are reference weight factor coefficients set for the total peak spacing and the number of peaks, respectively. The specific values ​​are determined by professionals in this field.

[0020] Optionally, step 102 specifically includes: arranging the equivalent diameters of each particle in ascending order, obtaining the particle diameters in the 60% and 10% ranges; dividing the 60% range by the 10% range to obtain the particle distribution non-uniformity coefficient value bj; obtaining the 30% range particle diameters; squaring the 30% range and dividing by the product of the 60% and 10% range particle diameters to obtain the particle distribution curvature coefficient value ql; and normalizing the particle distribution non-uniformity coefficient value bj, the particle distribution curvature coefficient value ql, and the peak value fz of the aggregate according to a set formula. The gradation value jz of the complex is obtained, where h3, h4, h5, h6, h7, h8, h9 and h10 are set correction factor coefficients to improve the accuracy of the calculation, h6>h3>h8, h7>h5>h10, and a3 is a natural constant.

[0021] Optionally, step 103 specifically includes: capturing motion images of each particle within a time interval Δt, calculating the angular displacement change Δθ between adjacent time interval points, dividing the angular displacement change Δθ by the time interval Δt to obtain the rotational angular velocity xz of each particle, calculating the centrifugal force lF of each particle using a formula, obtaining the adjacent contact diameter xd of the particles, obtaining the equivalent radius bj according to the formula bj=(xd×dz) / (xd+dz), obtaining the distance zj between the centers of two particles, obtaining the normal overlap fc according to the formula fc=xd+dz-zj, and using the formula The normal contact stiffness fj is obtained, where E is the elastic modulus of the particle, according to the established formula ZF = fj × fc. 3 / 2 The normal force ZF is obtained. The motion of the particle on the plane is simulated using discrete element method (DEM). When the particle begins to slide, the contact force at that moment is obtained using a micro-force sensor, and the static friction coefficient is derived using Coulomb's law. When the particle slides at a constant speed, the contact force at that moment is obtained using the micro-force sensor, and the kinetic friction coefficient is derived using Coulomb's law. When the particle is in a static friction state, the tangential contact force is less than or equal to the normal force ZF multiplied by the static friction coefficient. When the particle is in a kinetic friction state, the tangential contact force is equal to the normal force ZF multiplied by the kinetic friction coefficient. The tangential contact force under static and kinetic friction states is represented by Ft, according to the set formula. The particle contact force JF was calculated.

[0022] Optionally, step 104 specifically includes: according to the set formula The particle contact area jm is obtained, where υ is the Poisson's ratio of the particle and χ is the correction factor coefficient, which is used to improve the accuracy of the calculation. The particle contact force JF is divided by the particle contact area jm and multiplied by the correction factor coefficient to obtain the particle contact stress value. The contact stress value of each particle in the composite is obtained. The average contact stress value of each particle is processed to obtain the average contact stress of each particle. The maximum contact stress value of each particle in the composite is taken and divided by the average contact stress to obtain the stress anisotropy value of the composite.

[0023] Optionally, the specific steps for determining the filling degree of the structural fissure are as follows:

[0024] Step 201: Before the grouting operation, use an electromagnetic induction instrument to measure above the structural fracture area to obtain electromagnetic signals, calculate the amplitude of the signal and mark it as signal amplitude - fz1, and obtain the frequency value - pv1 by Fourier transforming the signal.

[0025] Step 202: After the grouting operation, use an electromagnetic induction instrument to measure the same location in the structural fracture area as in step 201 to obtain the electromagnetic signal, calculate the amplitude of the signal and mark it as signal amplitude two fz2, and obtain the frequency value two pv2 by Fourier transforming the signal.

[0026] Step 203: According to the set formula The degree of filling of the structural fracture is obtained as ctz, where fz max To set the ideal inductance amplitude for when the crack is completely filled, let pv max Set the ideal induction frequency so that the crack is completely filled.

[0027] Optionally, the specific steps for estimating the water level in the four-dimensional water table are as follows:

[0028] Step 301: Install four long-range hydrological observation wells in the mining area to acquire water level information in real time using water level gauges, and record the changes in water level after grouting. Establish a water level change curve with monitoring time as the x-axis and the water level of the four-crystal cut as the y-axis. Perform tangent processing on the curves at each monitoring time point, calculate the slope value of each tangent, and compare the slope value of each tangent with a set rapid change slope threshold. When a point is greater than or equal to the set rapid change slope threshold, it is identified as a rapid change marker point. This process is repeated sequentially from the start of grouting until the first rapid change marker point is identified and marked as a rapid change marker. Point 1: Compare the slope values ​​of each tangent with the set stable change slope threshold. When the point is less than or equal to the set stable change slope threshold, the point is identified as a stable change marker point. Starting from the first rapid change marker point, identify the first stable change marker point to the right and mark it as stable change marker point 1. The period from the start of grouting to the rapid change marker point 1 is defined as the slow change period of the four water level. The period from the rapid change marker point 1 to the stable change marker point 1 is defined as the rapid change period of the four water level. The period from the stable change marker point 1 to the end of monitoring is defined as the stable change period of the four water level.

[0029] Step 302: Obtain the duration of the period of slow water level change in the four-harvest area and mark it as the time lag value of water level change shz. Obtain the duration of the period of rapid water level change in the four-harvest area and add it to the time lag value of water level change shz to obtain the recovery time value of water level change hfz. Calculate the maximum and minimum water level values ​​during the stable water level change period in the four-harvest area and perform difference calculation to obtain the water level fluctuation amplitude value sbz. Obtain the slope value of the tangent at each monitoring time point during the rapid water level change period in the four-harvest area and average it to obtain the water level change rate value ssz. Calculate the water level initial state value scz by combining the time lag value of water level change shz, the recovery time value of water level change hfz, the water level fluctuation amplitude value sbz, and the water level change rate value ssz.

[0030] Step 303: By setting up several water level monitoring points during the stable water level change period of the four-crystal area, and arranging these monitoring points in the order they were set, the water level value of each monitoring point is obtained. The distance value between each monitoring point and the next monitoring point in the order is also obtained. The water level difference is obtained by subtracting the water level value of the next monitoring point from the water level value of each monitoring point. The water level difference is then divided by the distance value to obtain the water level gradient value of each monitoring point. Finally, the sum of the water level gradient values ​​of all monitoring points is divided by the number of monitoring points to obtain the water level gradient value td of the four-crystal area. i , i represents the monitoring time number, and the formula calculates the stable value wdz of the water level gradient in the four areas;

[0031] Step 304: Multiply the initial water level value scz and the stable water level gradient value wdz of the four-bearing water bodies by the weight factor coefficients of f1 and f2 respectively, and sum the products to obtain the estimated water level value of the four-bearing water bodies.

[0032] Optionally, the specific steps for alarm and judgment in step five are as follows:

[0033] When the comprehensive water level assessment value is less than the set comprehensive water level assessment threshold WT0, the corresponding coal mine working face treatment effect is unqualified and an alarm is issued;

[0034] When the comprehensive water level assessment value is greater than or equal to the set comprehensive water level assessment threshold WT0 and less than the set comprehensive water level assessment threshold WT1, the corresponding coal mine working face treatment effect is qualified.

[0035] When the comprehensive water level assessment value is greater than or equal to the set comprehensive water level assessment threshold WT1, the corresponding coal mine working face treatment effect is good.

[0036] The set comprehensive water level assessment threshold WT0 <WT1。

[0037] Optionally, in step four, obtaining the comprehensive water level assessment value through graphical construction specifically includes:

[0038] The combined heddle value, the filling degree of tectonic fissures, and the estimated water level of the four water-bearing states are normalized. An ellipse is constructed with the length of the combined heddle value as the major axis and the filling degree of tectonic fissures as the minor axis. A circle is constructed with the length of the estimated water level of the four water-bearing states as the radius. The center of the circle coincides with the right endpoint of the major axis of the ellipse. The non-overlapping area formed by the circle and the ellipse is identified and compared to obtain the comprehensive water level assessment value.

[0039] Compared with the prior art, the beneficial effects of the present invention are:

[0040] 1. This invention analyzes the particle size distribution and stress state of the cement grout and fracture zone formed after grouting treatment of the working face, and obtains the overall value of the composite. It also analyzes the fracture state before and after grouting to obtain the filling degree of the fracture. Furthermore, it monitors and analyzes the water level of the four water-bearing zones in real time to obtain the estimated water level value. This invention takes into account the particle size distribution and stress of the composite after formation, as well as the fracture state, and considers the real-time monitoring and analysis of the four water-bearing zones, enabling a multi-angle and in-depth analysis of the overall water level state.

[0041] 2. This invention assesses the comprehensive water level status of the working face by combining the overall water level value, the filling degree value of structural fissures, and the estimated water level value. The comprehensive water level status value is obtained through graphic construction and processing. The comprehensive water level status value is received and compared with the set comprehensive water level status value threshold. Corresponding alarms and judgments are made, which can promptly determine the comprehensive water level status and issue alarms for abnormal status, thereby reducing safety accidents in coal mining. Attached Figure Description

[0042] To more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. The following drawings are not drawn to scale according to the actual size, but are intended to show the main idea of ​​the present invention.

[0043] Figure 1 This is a flowchart of the method of the present invention. Detailed Implementation

[0044] The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are also within the scope of protection of the present invention.

[0045] The method for real-time water level tracking and observation based on the treatment of effluent from the working face in this invention specifically includes:

[0046] Step 1: After the cement grouting treatment of the working face, the combination formed by the cement grout and the crack zone is analyzed for particle size distribution and particle stress state to obtain the composite value of the combination.

[0047] Step 2: Analyze the state of the structural cracks before and after grouting to obtain the filling degree value of the structural cracks;

[0048] Step 3: Conduct real-time monitoring and analysis of the water level of the four cut surfaces at the working face to obtain an estimated water level status.

[0049] Step 4: Evaluate the overall water level status of the working face water treatment based on the combined body morphology value, the filling degree value of structural fissures, and the estimated water level value of the four water-bearing structures, and obtain the overall water level status value through graphic construction processing;

[0050] Step 5: Receive the comprehensive water level assessment value, compare and analyze it with the set comprehensive water level assessment threshold, and issue corresponding alarms and judgments.

[0051] The specific operational steps based on the systolic symptom values ​​are as follows:

[0052] Step 101: The cement grout and the fracture zone formed by the grouting treatment of the working face after water discharge are imaged using an electron microscope. The image is converted into a grayscale image. Threshold segmentation using the maximum inter-class variance method is used to separate all particles from the image background. The pixel area of ​​each particle is calculated, and the pixel area is converted to the actual area according to a set ratio coefficient to obtain the actual area sm of each particle. The equivalent diameter dz of each particle is calculated by formula. The equivalent diameter of each particle is plotted on the x-axis with the particle diameter size and the frequency of particle occurrence in each size interval as the y-axis to establish a particle size frequency distribution histogram. The frequency number in each size interval is counted and compared with the set reference peak frequency threshold for each size interval. When the frequency number of the size interval is greater than the reference peak frequency threshold of the size interval, the size interval is determined as a peak interval. The distance between each adjacent peak interval is counted and represented by the number of intervals. The total peak distance value zj is obtained by summing the distance between the intervals and the number of peaks fs. The peak value fz of the composite is obtained by calculating the total peak distance value zj and the number of peaks fs.

[0053] Step 102: Arrange the equivalent diameters of each particle in ascending order to obtain the particle diameters in the 60% and 10% ranges. Divide the particle diameters in the 60% range by the particle diameters in the 10% range to obtain the particle distribution non-uniformity coefficient value bj. Obtain the particle diameters in the 30% range, square them, and divide them by the product of the particle diameters in the 60% and 10% ranges to obtain the particle distribution curvature coefficient value ql. Calculate the particle distribution non-uniformity coefficient value bj, the particle distribution curvature coefficient value ql, and the peak value fz of the aggregate using the formula to obtain the aggregate gradation value jz.

[0054] Step 103: According to the set formula The particle contact force JF is calculated, ZF is the normal force; lF is the centrifugal force of each particle; Ft represents the tangential contact force under static and dynamic friction conditions.

[0055] Step 104: Calculate the particle contact area jm, divide the particle contact force JF by the particle contact area jm, and multiply by the correction factor coefficient to obtain the particle contact stress value. Obtain the contact stress value of each particle in the composite, process the average contact stress value of each particle to obtain the average contact stress value of each particle, take the maximum contact stress value of each particle in the composite and divide it by the average contact stress value to obtain the stress distribution value of the composite.

[0056] Step 105: Multiply the gradation value jz of the composite by the weighting factor coefficient g1 to obtain product coefficient one; multiply the stress concentration anisotropy value of the composite by the weighting factor coefficient g2 to obtain product coefficient two; divide product coefficient one by product coefficient two to obtain the composite gradation value.

[0057] The specific operational steps for stress concentration anisotropy in the composite are as follows:

[0058] Step 103: By capturing motion images of each particle within a time interval Δt, calculate the angular displacement change Δθ between adjacent time intervals. Divide the angular displacement change Δθ by the time interval Δt to obtain the rotational angular velocity xz of each particle. Calculate the centrifugal force lF of each particle using the formula. Obtain the adjacent contact diameter xd of the particles. According to the formula bj=(xd×dz) / (xd+dz), obtain the equivalent radius bj. Obtain the distance zj between the centers of two particles. According to the formula fc=xd+dz-zj, obtain the normal overlap fc. The normal contact stiffness fj is obtained, where E is the elastic modulus of the particle, according to the established formula ZF = fj × fc. 3 / 2 The normal force ZF is obtained. The motion of the particle on the plane is simulated using discrete element method (DEM). When the particle begins to slide, the contact force at that moment is obtained using a micro-force sensor, and the static friction coefficient is derived using Coulomb's law. When the particle slides at a constant speed, the contact force at that moment is obtained using the micro-force sensor, and the kinetic friction coefficient is derived using Coulomb's law. When the particle is in a static friction state, the tangential contact force is less than or equal to the normal force ZF multiplied by the static friction coefficient. When the particle is in a kinetic friction state, the tangential contact force is equal to the normal force ZF multiplied by the kinetic friction coefficient. The tangential contact force under static and kinetic friction states is represented by Ft. According to the established formula... The particle contact force JF was calculated.

[0059] Step 104: According to the set formula The particle contact area jm is obtained, where υ is the Poisson's ratio of the particle and χ is the correction factor coefficient, which is used to improve the accuracy of the calculation. The particle contact force JF is divided by the particle contact area jm and multiplied by the correction factor coefficient to obtain the particle contact stress value. The contact stress value of each particle in the composite is obtained. The average contact stress value of each particle is processed to obtain the average contact stress of each particle. The maximum contact stress value of each particle in the composite is taken and divided by the average contact stress to obtain the stress anisotropy value of the composite.

[0060] The specific steps for calculating the gradation value of the complex are as follows:

[0061] Step 102: Arrange the equivalent diameters of each particle in ascending order, obtaining the particle diameters in the 60% and 10% ranges. Divide the 60% particle diameters by the 10% particle diameters to obtain the particle distribution non-uniformity coefficient value bj. Obtain the particle diameters in the 30% range, square them, and divide by the product of the 60% and 10% particle diameters to obtain the particle distribution curvature coefficient value ql. Normalize the particle distribution non-uniformity coefficient value bj, the particle distribution curvature coefficient value ql, and the peak value fz of the aggregate according to the set formula. The gradation value jz of the complex is obtained, where h3, h4, h5, h6, h7, h8, h9 and h10 are set correction factor coefficients to improve the accuracy of the calculation, h6>h3>h8, h7>h5>h10, and a3 is a natural constant.

[0062] The specific steps for calculating the peak values ​​of the complex are as follows:

[0063] Step 101: The cement grout and fracture zone bond formed after grouting treatment at the working face are imaged using a scanning electron microscope. The image is converted to grayscale, and all particles are separated from the image background using threshold segmentation with the maximum inter-class variance method. The pixel area of ​​each particle is calculated, and the pixel area is converted to the actual area according to a set ratio coefficient to obtain the actual area sm of each particle. The formula yields the equivalent diameter dz of each particle. A particle size frequency distribution histogram is constructed, plotting the particle diameter on the x-axis and the frequency of particle occurrence within each size interval on the y-axis. The frequency count within each size interval is counted and compared with a predefined reference peak frequency threshold. If the frequency count of a size interval exceeds its reference peak frequency threshold, that interval is identified as a peak interval. The distance between adjacent peak intervals is counted, represented by the number of intervals, and summed to obtain the total peak interval value zj. The number of peaks fs is also counted. Finally, the total peak interval value zj and the number of peaks fs are normalized according to the predefined formula. The peak value fz of the composite is obtained by calculation, where a1 and a2 are natural constants, and a1>a2, and h1 and h2 are reference weight factor coefficients set for the total peak spacing and the number of peaks, respectively. The specific values ​​are determined by professionals in this field.

[0064] The specific steps for determining the filling degree of the fracture are as follows:

[0065] Step 201: Before the grouting operation, use an electromagnetic induction instrument to measure above the structural fracture area to obtain electromagnetic signals, calculate the amplitude of the signal and mark it as signal amplitude - fz1, and obtain the frequency value - pv1 by Fourier transforming the signal.

[0066] Step 202: After the grouting operation, use an electromagnetic induction instrument to measure the same location in the structural fracture area as in step 201 to obtain the electromagnetic signal, calculate the amplitude of the signal and mark it as signal amplitude two fz2, and obtain the frequency value two pv2 by Fourier transforming the signal.

[0067] Step 203: Analyze and obtain the filling degree value ctz of the structural fracture. Based on the established formula... The degree of filling of the structural fracture is obtained as ctz, where fz max To set the ideal inductance amplitude for when the crack is completely filled, let pv max Set the ideal induction frequency so that the crack is completely filled.

[0068] The specific steps for estimating the water level of four water-bearing structures are as follows:

[0069] Step 303: By setting up several water level monitoring points during the stable water level change period of the four-crystal area, and arranging these monitoring points in the order they were set, the water level value of each monitoring point is obtained. The distance value between each monitoring point and the next monitoring point in the order is also obtained. The water level difference is obtained by subtracting the water level value of the next monitoring point from the water level value of each monitoring point. The water level difference is then divided by the distance value to obtain the water level gradient value of each monitoring point. Finally, the sum of the water level gradient values ​​of all monitoring points is divided by the number of monitoring points to obtain the water level gradient value td of the four-crystal area. i , i represents the monitoring time number, and the formula calculates the stable value wdz of the water level gradient in the four areas;

[0070] Step 304: Multiply the initial water level value scz and the stable water level gradient value wdz of the four-bearing water bodies by the weight factor coefficients of f1 and f2 respectively, and sum the products to obtain the estimated water level value of the four-bearing water bodies.

[0071] The specific operational steps for setting the initial water level value are as follows:

[0072] Step 301: Install four long-range hydrological observation wells in the mining area to acquire water level information in real time using water level gauges, and record the changes in water level after grouting. Establish a water level change curve with monitoring time as the x-axis and the water level of the four-crystal cut as the y-axis. Perform tangent processing on the curves at each monitoring time point, calculate the slope value of each tangent, and compare the slope value of each tangent with a set rapid change slope threshold. When a point is greater than or equal to the set rapid change slope threshold, it is identified as a rapid change marker point. This process is repeated sequentially from the start of grouting until the first rapid change marker point is identified and marked as a rapid change marker. Point 1: Compare the slope values ​​of each tangent with the set stable change slope threshold. When the point is less than or equal to the set stable change slope threshold, the point is identified as a stable change marker point. Starting from the first rapid change marker point, identify the first stable change marker point to the right and mark it as stable change marker point 1. The period from the start of grouting to the rapid change marker point 1 is defined as the slow change period of the four water level. The period from the rapid change marker point 1 to the stable change marker point 1 is defined as the rapid change period of the four water level. The period from the stable change marker point 1 to the end of monitoring is defined as the stable change period of the four water level.

[0073] Step 302: Obtain the duration of the period of slow water level change in the four-harvest area and mark it as the time lag value of water level change shz. Obtain the duration of the period of rapid water level change in the four-harvest area and add it to the time lag value of water level change shz to obtain the recovery time value of water level change hfz. Calculate the maximum and minimum water level values ​​during the stable water level change period in the four-harvest area and perform difference calculation to obtain the water level fluctuation amplitude value sbz. Obtain the slope value of the tangent at each monitoring time point during the rapid water level change period in the four-harvest area and average it to obtain the water level change rate value ssz. Calculate the initial water level value scz by combining the time lag value of water level change shz, the recovery time value of water level change hfz, the water level fluctuation amplitude value sbz, and the water level change rate value ssz.

[0074] The specific steps for handling the corresponding alarms and judgments are as follows:

[0075] When the comprehensive water level assessment value is less than the set comprehensive water level assessment threshold WT0, the corresponding coal mine working face treatment effect is unqualified, and an alarm is issued, which is displayed on the screen set up at the coal mine construction entrance in text message 1; when the comprehensive water level assessment value is greater than or equal to the set comprehensive water level assessment threshold WT0 and less than the set comprehensive water level assessment threshold WT1, the corresponding coal mine working face treatment effect is qualified, and the alarm is displayed on the screen set up at the coal mine construction entrance in text message 2; when the comprehensive water level assessment value is greater than or equal to the set comprehensive water level assessment threshold WT1, the corresponding coal mine working face treatment effect is good, and the alarm is displayed on the screen set up at the coal mine construction entrance in text message 3.

[0076] Example 1:

[0077] like Figure 1 As shown, the method for real-time tracking and observation of water level based on the treatment of effluent from the working face includes the following steps:

[0078] Step 1: After the cement grouting treatment of the working face, the combination formed by the cement grout and the crack zone is analyzed for particle size distribution and particle stress state to obtain the composite value of the combination.

[0079] Step 2: Analyze the state of the structural cracks before and after grouting to obtain the filling degree value of the structural cracks;

[0080] Step 3: Conduct real-time monitoring and analysis of the water level of the four cut surfaces at the working face to obtain an estimated water level status.

[0081] Step 4: Receive the overall water level value of the combined body, the filling degree value of the structural fissures, and the estimated water level value of the four water-bearing structures to evaluate the overall water level status of the working face water treatment and obtain the overall water level status rating.

[0082] Step 5: Receive the comprehensive water level assessment value, compare and analyze it with the set comprehensive water level assessment threshold, and issue corresponding alarms and judgments.

[0083] The specific steps for calculating the systolic symmetry value are as follows:

[0084] Step 101: The cement grout and fracture zone bond formed after grouting treatment at the working face are imaged using a scanning electron microscope. The image is converted to grayscale, and all particles are separated from the image background using threshold segmentation with the maximum inter-class variance method. The pixel area of ​​each particle is calculated, and the pixel area is converted to the actual area according to a set ratio coefficient to obtain the actual area sm of each particle. The formula is used to obtain the equivalent diameter dz of each particle. A particle size frequency distribution histogram is constructed, with the particle diameter size as the x-axis and the frequency of particle occurrence in each size interval as the y-axis. The particle diameter sizes are categorized as 0-2µm, 2-4µm, 4-8µm, 8-16µm, 16-30µm, 30-50µm, 50-100µm, 100-200µm, 200-400µm, and above 400µm, where µm represents micrometers. The frequency count within each size interval is counted and compared with a set reference peak frequency threshold. If the frequency count of a size interval is greater than the reference peak frequency threshold, that size interval is determined as a peak interval. The distance between adjacent peak intervals is counted, represented by the number of intervals, and summed to obtain the total peak interval value zj. The number of peaks fs is also counted. The total peak interval value zj and the number of peaks fs are then normalized according to the set formula. The peak value fz of the aggregate is obtained by calculation, where a1 and a2 are natural constants, and a1>a2. h1 and h2 are reference weight factor coefficients set for the total peak spacing and the number of peaks, respectively. The specific values ​​are determined by professionals in this field. It should be noted that when the number of peaks is larger, the spacing between peaks is smaller, the particle gradation is worse, and the peak value of the aggregate is smaller.

[0085] Step 102: Arrange the equivalent diameters of each particle in ascending order, obtaining the particle diameters in the 60% and 10% ranges. Divide the 60% particle diameters by the 10% particle diameters to obtain the particle distribution non-uniformity coefficient value bj. Obtain the particle diameters in the 30% range, square them, and divide by the product of the 60% and 10% particle diameters to obtain the particle distribution curvature coefficient value ql. Normalize the particle distribution non-uniformity coefficient value bj, the particle distribution curvature coefficient value ql, and the peak value fz of the aggregate according to the set formula. The graded value jz of the aggregate is obtained, where h3, h4, h5, h6, h7, h8, h9, and h10 are set correction factor coefficients used to improve the accuracy of the calculation, h6>h3>h8, h7>h5>h10, and a3 is a natural constant with a specific value of 51.102. As can be seen from the formula, the curvature coefficient is between 1 and 3. The larger the value of the non-uniformity coefficient of particle distribution, the larger the graded value of the aggregate.

[0086] Step 103: By capturing motion images of each particle within a time interval Δt, calculate the angular displacement change Δθ between adjacent time intervals. Divide the angular displacement change Δθ by the time interval Δt to obtain the rotational angular velocity xz of each particle. Then, use the formula lF = λ × xz. 2 ×dz / 2, we obtain the centrifugal force lF of each particle, where λ is an adjustment factor coefficient used to improve the accuracy of the calculation. We obtain the adjacent contact diameter xd of the particles, and according to the formula bj=(xd×dz) / (xd+dz), we obtain the equivalent radius bj. We also obtain the distance zj between the centers of the two particles, and according to the formula fc=xd+dz-zj, we obtain the normal overlap fc. From the formula... The normal contact stiffness fj is obtained, where E is the elastic modulus of the particle, according to the established formula ZF = fj × fc. 3 / 2 The normal force ZF is obtained. The motion of the particle on the plane is simulated using discrete element method (DEM). When the particle begins to slide, the contact force at that moment is obtained using a micro-force sensor, and the static friction coefficient is derived using Coulomb's law. When the particle slides at a constant speed, the contact force at that moment is obtained using the micro-force sensor, and the kinetic friction coefficient is derived using Coulomb's law. When the particle is in a static friction state, the tangential contact force is less than or equal to the normal force ZF multiplied by the static friction coefficient. When the particle is in a kinetic friction state, the tangential contact force is equal to the normal force ZF multiplied by the kinetic friction coefficient. The tangential contact force under static and kinetic friction states is represented by Ft, according to the set formula. The particle contact force JF was calculated.

[0087] Step 104: According to the set formula The particle contact area jm is obtained, where υ is the Poisson's ratio of the particle and χ is the correction factor coefficient, which is used to improve the accuracy of the calculation. The particle contact force JF is divided by the particle contact area jm and multiplied by the correction factor coefficient to obtain the particle contact stress value. The contact stress value of each particle in the composite is obtained. The average contact stress value of each particle is processed to obtain the average contact stress value of each particle. The maximum contact stress value of each particle in the composite is taken and divided by the average contact stress value to obtain the stress anisotropy value of the composite.

[0088] Step 105: Multiply the gradation value jz of the composite by the weighting factor coefficient g1 to obtain product coefficient one; multiply the stress concentration anisotropy value of the composite by the weighting factor coefficient g2 to obtain product coefficient two; divide product coefficient one by product coefficient two to obtain the composite gradation value.

[0089] The specific steps for determining the degree of filling of the structural fractures are as follows:

[0090] Step 201: Before the grouting operation, use an electromagnetic induction instrument to measure above the structural fracture area to obtain electromagnetic signals, calculate the amplitude of the signal and mark it as signal amplitude - fz1, and obtain the frequency value - pv1 by Fourier transforming the signal.

[0091] Step 202: After the grouting operation, use an electromagnetic induction instrument to measure the same location in the structural fracture area as in step 201 to obtain the electromagnetic signal, calculate the amplitude of the signal and mark it as signal amplitude two fz2, and obtain the frequency value two pv2 by Fourier transforming the signal.

[0092] Step 203: According to the set formula The degree of filling of the structural fracture is obtained as ctz, where fz max To set the ideal inductance amplitude for when the crack is completely filled, let pv max Set the ideal induction frequency so that the crack is completely filled.

[0093] The specific steps for calculating the water level status of the four water-bearing sections are as follows:

[0094] Step 301: Install four long-range hydrological observation wells in the mining area to acquire water level information in real time using water level gauges, and record the changes in water level after grouting. Establish a water level change curve with monitoring time as the x-axis and the water level of the four-crystal cut as the y-axis. Perform tangent processing on the curves at each monitoring time point, calculate the slope value of each tangent, and compare the slope value of each tangent with a set rapid change slope threshold. When a point is greater than or equal to the set rapid change slope threshold, it is identified as a rapid change marker point. This process is repeated sequentially from the start of grouting until the first rapid change marker point is identified and marked as a rapid change marker. Point 1: Compare the slope values ​​of each tangent with the set stable change slope threshold. When the point is less than or equal to the set stable change slope threshold, the point is identified as a stable change marker point. Starting from the first rapid change marker point, identify the first stable change marker point to the right and mark it as stable change marker point 1. The period from the start of grouting to the rapid change marker point 1 is defined as the slow change period of the four water level. The period from the rapid change marker point 1 to the stable change marker point 1 is defined as the rapid change period of the four water level. The period from the stable change marker point 1 to the end of monitoring is defined as the stable change period of the four water level.

[0095] Step 302: Obtain the duration of the period of slow water level change in the four-catchment area, and mark it as the time lag value shz of water level change. Obtain the duration of the period of rapid water level change in the four-catchment area, and add it to the time lag value shz of water level change to obtain the recovery time value hfz of water level change. Calculate the maximum and minimum water level values ​​during the stable water level change period in the four-catchment area, and calculate the difference to obtain the water level fluctuation amplitude value sbz. Obtain the slope of the tangent at each monitoring time point during the rapid water level change period in the four-catchment area, and average it to obtain the water level change rate value ssz. Then, calculate the time lag value s of water level change. The water level lag value hz, the recovery time value hfz, the water level fluctuation amplitude value sbz, and the water level change rate value ssz are normalized. The initial water level value scz is calculated according to the set formula scz=cot4×ssz / (cot1×shz+cot2×hfz+cot3×sbz+1.202), where cot1, cot2, cot3, and cot4 are the set weight factor coefficients for the water level change time lag value, the water level change recovery time value, the water level fluctuation amplitude value, and the water level change rate value, respectively. The specific values ​​are determined by professionals in this field.

[0096] Step 303: By setting up several water level monitoring points during the stable water level change period of the four-crystal area, and arranging these monitoring points in the order they were set, the water level value of each monitoring point is obtained. The distance value between each monitoring point and the next monitoring point in the order is also obtained. The water level difference is obtained by subtracting the water level value of the next monitoring point from the water level value of each monitoring point. The water level difference is then divided by the distance value to obtain the water level gradient value of each monitoring point. Finally, the sum of the water level gradient values ​​of all monitoring points is divided by the number of monitoring points to obtain the water level gradient value td of the four-crystal area. i 'i' represents the monitoring time number, i = 1, 2, 3, ..., n1, where n1 is the maximum value of the monitoring time number, determined according to the set formula. The stable values ​​of the four water level gradients, wdz, are obtained. A1 is the average value of the four water level gradients, and η is the correction factor coefficient value, which is used to improve the accuracy of the calculation.

[0097] Step 304: Multiply the initial water level value scz and the stable water level gradient value wdz of the four-bearing water bodies by the weight factor coefficients of f1 and f2 respectively, and sum the products to obtain the estimated water level value of the four-bearing water bodies;

[0098] The specific steps for calculating the comprehensive water level assessment value are as follows:

[0099] The combined heddle value, the filling degree of structural fissures, and the four water level estimates are normalized. An ellipse is constructed with the length of the combined heddle value as the major axis and the filling degree of structural fissures as the minor axis. A circle is constructed with the length of the four water level estimates as the radius. The center of the circle coincides with the right endpoint of the major axis of the ellipse. The non-overlapping area formed by the circle and the ellipse is identified and compared to obtain the comprehensive water level assessment.

[0100] The specific steps for generating the corresponding alarms and judgments are as follows:

[0101] Step 501: Install a display screen at the coal mine construction entrance to show the overall water level status of the working face;

[0102] Step 502: When the comprehensive water level assessment value is less than the set comprehensive water level assessment threshold WT0, the corresponding coal mine working face treatment effect is unsatisfactory, and an alarm is issued. The text message "The treatment effect of this working face is poor; construction should be stopped immediately, and a safety investigation team should be established promptly" is sent to the display screen set up at the coal mine construction entrance. When the comprehensive water level assessment value is greater than or equal to the set comprehensive water level assessment threshold WT0 but less than the set comprehensive water level assessment threshold WT1, the corresponding coal mine working face treatment effect is satisfactory. The text message "The treatment effect of this working face is satisfactory; increase the frequency of water level monitoring" is sent to the display screen set up at the coal mine construction entrance. When the comprehensive water level assessment value is greater than or equal to the set comprehensive water level assessment threshold WT1, the corresponding coal mine working face treatment effect is good. The text message "The treatment effect of this working face is good" is sent to the display screen set up at the coal mine construction entrance. The set comprehensive water level assessment threshold WT0 is not specified in the original text. <WT1。

[0103] The foregoing description is illustrative of the invention and should not be construed as limiting it. Although several exemplary embodiments of the invention have been described, those skilled in the art will readily understand that many modifications can be made to the exemplary embodiments without departing from the novel teachings and advantages of the invention. Therefore, all such modifications are intended to be included within the scope of the invention as defined in the claims. It should be understood that the foregoing description is illustrative of the invention and should not be construed as limiting it to the specific embodiments disclosed, and modifications to the disclosed embodiments and other embodiments are intended to be included within the scope of the appended claims.

Claims

1. A method for real-time tracking and observation of water level based on effluent treatment at a working face, characterized in that, Includes the following steps: Step 1: Analyze the particle size distribution and stress state of the bond formed between the cement grout and the fracture zone after grouting treatment at the working face to obtain the bond's helix value. The method for obtaining the bond's helix value is as follows: Step 101: The cement grout and fracture zone bond formed after the grouting treatment of the working face are imaged using a scanning electron microscope. The image is converted to grayscale, and all particles are separated from the image background using threshold segmentation with the Otsu's method. The pixel area of ​​each particle is calculated, and the pixel area is converted to the actual area according to a set ratio coefficient to obtain the actual area of ​​each particle. The formula calculates the equivalent diameter of each particle. A particle size frequency distribution histogram is constructed by plotting the equivalent diameter of each particle on the x-axis and the frequency of particle occurrence within each size interval on the y-axis. The frequency count within each size interval is counted and compared with a predefined reference peak frequency threshold for each size interval. If the frequency count of a size interval exceeds the reference peak frequency threshold, that size interval is identified as a peak interval. The distance between adjacent peak intervals is then calculated, represented by the number of intervals, and summed to obtain the total peak interval distance. And count the number of peak values. and the total peak spacing value and peak values The peak values ​​of the complex were calculated. ; Step 102: Arrange the equivalent diameters of all particles in ascending order, and obtain the particle diameters in the 60th percentile and the 10th percentile. Divide the particle diameters in the 60th percentile by the particle diameters in the 10th percentile to obtain the particle distribution non-uniformity coefficient. Obtain the particle diameters that rank in the top 30%. Square these top 30% particle diameters and divide by the product of the top 60% and top 10% particle diameters to obtain the curvature coefficient value of the particle distribution. The coefficient of non-uniformity of particle distribution Curvature coefficient of particle distribution Peak values ​​of the complex The formula calculates the gradation value of the complex. ; Step 103: According to the set formula Calculated particle contact force , Positive pressure; The centrifugal force of each particle; It represents the tangential contact force under static and dynamic friction conditions; Step 104: Calculate the particle contact area , to increase particle contact force Divide by particle contact area Multiply by the correction factor to obtain the particle contact stress value, obtain the contact stress value of each particle in the composite, process the average contact stress value of each particle to obtain the average contact stress value of each particle, take the maximum contact stress value of each particle in the composite and divide it by the average contact stress value to obtain the stress concentration value of the composite. Step 105: Calculate the gradation values ​​of the aggregate. Multiply by the weighting factor coefficient This yields a product coefficient of 1, which is used to multiply the stress set anisotropy of the composite by the weighting factor coefficient. We obtain the product coefficient 2, and divide the product coefficient 1 by the product coefficient 2 to obtain the complex value of the complex. Step 2: Analyze the state of the structural cracks before and after grouting to obtain the filling degree value of the structural cracks; Step 3: Conduct real-time monitoring and analysis of the water level of the four cut surfaces at the working face to obtain an estimated water level status. Step 4: Evaluate the overall water level status of the working face water treatment based on the combined body morphology value, the filling degree value of structural fissures, and the estimated water level value of the four water-bearing structures, and obtain the overall water level status value through graphic construction processing; Step 5: Receive the comprehensive water level assessment value, compare and analyze it with the set comprehensive water level assessment threshold, and issue corresponding alarms and judgments.

2. The method for real-time tracking and observation of water level based on the treatment of effluent from the working face according to claim 1, characterized in that, Step 101 specifically includes: The cement slurry and fracture zone bond formed after grouting treatment at the working face were imaged using a scanning electron microscope. The image was converted to grayscale, and threshold segmentation using the maximum inter-class variance method was applied to separate all particles from the image background. The pixel area of ​​each particle was calculated, and the pixel area was converted to the actual area according to a set ratio coefficient to obtain the actual area of ​​each particle. , by setting The formula yields the equivalent diameter of each particle. A particle size frequency distribution histogram is constructed by plotting the equivalent diameter of each particle on the x-axis and the frequency of particle occurrence within each size interval on the y-axis. The frequency count within each size interval is counted and compared with a predefined reference peak frequency threshold for each size interval. If the frequency count of a size interval exceeds the reference peak frequency threshold, that size interval is identified as a peak interval. The distance between adjacent peak intervals is then calculated, represented by the number of intervals, and summed to obtain the total peak interval distance. And count the number of peak values. and the total peak spacing value and peak values Normalization is performed based on the set formula. The peak values ​​of the complex were calculated. ,in and They are the natural constants, and , and These are the reference weight factor coefficients set for the total peak spacing and the number of peaks, respectively.

3. The method for real-time tracking and observation of water level based on the treatment of effluent from the working face according to claim 1 or 2, characterized in that, Step 102 specifically includes: arranging the equivalent diameters of each particle in ascending order, obtaining the particle diameters in the 60% and 10% percentiles, and dividing the particle diameters in the 60% percentiles by the particle diameters in the 10% percentiles to obtain the particle distribution non-uniformity coefficient value. Obtain the particle diameters that rank in the top 30%. Square these top 30% particle diameters and divide by the product of the top 60% and top 10% particle diameters to obtain the curvature coefficient value of the particle distribution. The coefficient of non-uniformity of particle distribution Curvature coefficient of particle distribution Peak values ​​of the complex Normalization is performed based on the set formula. The gradation values ​​of the complex were obtained. h3, h4, h5, h6, h7, h8, h9, and h10 are the set correction factor coefficients used to improve the accuracy of the calculation, h6>h3>h8, h7>h5>h10, and a3 is a natural constant.

4. The method for real-time tracking and observation of water level based on the treatment of effluent from the working face according to claim 1 or 2, characterized in that, Step 103 specifically includes: through a time interval Internally, images of the motion of each particle are captured, and the change in angular displacement at adjacent time intervals is calculated. The change in angular displacement Divide by time interval The rotational angular velocity of each particle is obtained. The centrifugal force of each particle is calculated using the formula. Obtain the adjacent contact diameter of the particles. According to the formula The equivalent radius is obtained. And obtain the distance between the centers of the two particles. According to the formula The normal overlap is obtained. From the formula The normal contact stiffness is obtained. ,in The elastic modulus of the particles is determined according to the established formula. Receive positive pressure The motion of particles on a plane is simulated using discrete element method (DEM). When a particle begins to slide, the contact force at that moment is obtained using a micro-force sensor, and the static friction coefficient is deduced using Coulomb's law. When the particle slides at a constant speed, the contact force at that moment is obtained using a micro-force sensor, and the kinetic friction coefficient is deduced using Coulomb's law. When the particle is in a static friction state, the tangential contact force is less than or equal to the normal force. Multiplying by the static friction coefficient, when the particles are in a state of dynamic friction, the tangential contact force equals the normal force. Multiply by the coefficient of kinetic friction to convert the tangential contact force under both static and kinetic friction conditions into a coefficient of dynamic friction. This indicates that, according to the set formula The particle contact force was calculated. .

5. The method for real-time tracking and observation of water level based on the treatment of effluent from the working face according to claim 1 or 2, characterized in that, Step 104 specifically includes: according to the set formula To obtain the particle contact area Where υ is the Poisson's ratio of the particle, To correct the factor coefficients and improve the accuracy of calculations, the particle contact force is... Divide by particle contact area Multiply by the correction factor to obtain the particle contact stress value, obtain the contact stress value of each particle in the composite, process the average contact stress value of each particle to obtain the average contact stress value of each particle, take the maximum contact stress value of each particle in the composite and divide it by the average contact stress value to obtain the stress concentration value of the composite.

6. The method for real-time tracking and observation of water level based on the treatment of effluent from the working face according to claim 1 or 2, characterized in that, The specific operational steps for determining the degree of filling of the structural fissures are as follows: Step 201: Before grouting, use an electromagnetic induction instrument to measure above the structural fracture area to obtain electromagnetic signals, calculate the signal amplitude, and mark it as signal amplitude one. The signal is subjected to Fourier transform to obtain the frequency value. ; Step 202: After grouting, use an electromagnetic induction instrument to measure the electromagnetic signal at the same location as in step 201 in the structural fracture area, calculate the signal amplitude, and mark it as signal amplitude two. The signal is subjected to Fourier transform to obtain the frequency value. ; Step 203: According to the set formula The degree of filling of the structural fractures was obtained. ,in To set the ideal inductance amplitude for when the crack is completely filled, Set the ideal induction frequency so that the crack is completely filled.

7. The method for real-time tracking and observation of water level based on the treatment of effluent from the working face according to claim 1 or 2, characterized in that, The specific steps for estimating the water level in the four water-bearing states are as follows: Step 301: Install four long-range hydrological observation wells in the mining area to acquire water level information in real time using water level gauges, and record the changes in water level after grouting. Establish a water level change curve with monitoring time as the x-axis and the water level of the four-crystal cut as the y-axis. Perform tangent processing on the curves at each monitoring time point, calculate the slope value of each tangent, and compare the slope value of each tangent with a set rapid change slope threshold. When a point is greater than or equal to the set rapid change slope threshold, it is identified as a rapid change marker point. This process is repeated sequentially from the start of grouting until the first rapid change marker point is identified and marked as a rapid change marker. Point 1: Compare the slope values ​​of each tangent with the set stable change slope threshold. When the point is less than or equal to the set stable change slope threshold, the point is identified as a stable change marker point. Starting from the first rapid change marker point, identify the first stable change marker point to the right and mark it as stable change marker point 1. The period from the start of grouting to the rapid change marker point 1 is defined as the slow change period of the four water level. The period from the rapid change marker point 1 to the stable change marker point 1 is defined as the rapid change period of the four water level. The period from the stable change marker point 1 to the end of monitoring is defined as the stable change period of the four water level. Step 302: Obtain the duration of the four periods of slow water level change and mark them as the time lag value of water level change. Obtain the duration of the period of rapid water level change in the four-crystal area and the time lag value of the water level change. Add them together to get the recovery time value of the water level change. The maximum and minimum water level values ​​during the stable water level period in the four-in-one area were statistically analyzed, and the difference was calculated to obtain the water level fluctuation range. The slope values ​​of the tangents at each monitoring time point during the rapid water level change period were obtained, and the average value was processed to obtain the water level change rate value. The time lag value of water level change Recovery time value of water level change Water level fluctuation range and water level change rate value The initial water level value was obtained through calculation. ; Step 303: By setting up several water level monitoring points during the stable water level change period of the four-catchment area, and arranging these monitoring points in the order they were set, the water level value of each monitoring point is obtained. The distance value between each monitoring point and the next monitoring point in the order is also obtained. The water level difference is obtained by subtracting the water level value of the next monitoring point from the water level value of each monitoring point. The water level gradient value of each monitoring point is obtained by dividing the water level difference by the distance value. The sum of the water level gradient values ​​of each monitoring point is divided by the number of monitoring points to obtain the water level gradient value of the four-catchment area. , The monitoring time number is used to calculate the stable value of the water level gradient in the four sections. ; Step 304: Initial water level value and the stable value of the four water level gradients Multiply by respectively and The weighting factor coefficients are calculated, and the product is summed to obtain the four water level estimates.

8. The method for real-time tracking and observation of water level based on the treatment of effluent from the working face according to claim 1 or 2, characterized in that, The specific steps for alarm and judgment in step five are as follows: When the comprehensive water level assessment value is less than the set comprehensive water level assessment threshold At that time, the corresponding coal mine working face treatment effect is unsatisfactory, and an alarm is issued; When the comprehensive water level assessment value is greater than or equal to the set comprehensive water level assessment threshold The water level is below the set comprehensive assessment threshold. At that time, the corresponding coal mine working face treatment effect was satisfactory; When the comprehensive water level assessment value is greater than or equal to the set comprehensive water level assessment threshold The treatment effect on the corresponding coal mine working face is good; The set comprehensive water level assessment threshold .

9. The method for real-time tracking and observation of water level based on the treatment of effluent from the working face according to claim 1 or 2, characterized in that, In step four, obtaining the comprehensive water level assessment value through graphical construction specifically includes: The combined heddle value, the filling degree of tectonic fissures, and the estimated water level of the four water-bearing states are normalized. An ellipse is constructed with the length of the combined heddle value as the major axis and the filling degree of tectonic fissures as the minor axis. A circle is constructed with the length of the estimated water level of the four water-bearing states as the radius. The center of the circle coincides with the right endpoint of the major axis of the ellipse. The non-overlapping area formed by the circle and the ellipse is identified and compared to obtain the comprehensive water level assessment value.

Images