CN117052471A - Monitoring method and device for grouting effect of coal mine ground area - Google Patents
Monitoring method and device for grouting effect of coal mine ground area Download PDFInfo
- Publication number
- CN117052471A CN117052471A CN202311179475.7A CN202311179475A CN117052471A CN 117052471 A CN117052471 A CN 117052471A CN 202311179475 A CN202311179475 A CN 202311179475A CN 117052471 A CN117052471 A CN 117052471A
- Authority
- CN
- China
- Prior art keywords
- grouting
- measuring point
- ground
- coal mine
- frequency domain
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000003245 coal Substances 0.000 title claims abstract description 107
- 230000000694 effects Effects 0.000 title claims abstract description 77
- 238000000034 method Methods 0.000 title claims abstract description 64
- 238000012544 monitoring process Methods 0.000 title claims abstract description 37
- 238000001514 detection method Methods 0.000 claims abstract description 146
- 230000008859 change Effects 0.000 claims abstract description 35
- 238000004590 computer program Methods 0.000 claims description 21
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 13
- 239000002002 slurry Substances 0.000 claims description 11
- 230000007704 transition Effects 0.000 claims description 11
- 238000012795 verification Methods 0.000 claims description 11
- 230000006870 function Effects 0.000 claims description 10
- 238000010183 spectrum analysis Methods 0.000 claims description 9
- 238000003860 storage Methods 0.000 claims description 8
- 238000012937 correction Methods 0.000 claims description 7
- 238000007689 inspection Methods 0.000 claims description 6
- 238000012806 monitoring device Methods 0.000 claims description 5
- 230000035699 permeability Effects 0.000 claims description 5
- 230000006698 induction Effects 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 4
- 230000001939 inductive effect Effects 0.000 claims description 3
- 238000007569 slipcasting Methods 0.000 claims 1
- 238000009826 distribution Methods 0.000 description 13
- 238000005553 drilling Methods 0.000 description 12
- 238000005259 measurement Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 10
- 238000001914 filtration Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 238000009792 diffusion process Methods 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000010276 construction Methods 0.000 description 4
- 238000009499 grossing Methods 0.000 description 4
- 238000005065 mining Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000005358 geomagnetic field Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004164 analytical calibration Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 239000008398 formation water Substances 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
- E21F17/18—Special adaptations of signalling or alarm devices
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
The application provides a method and a device for monitoring grouting effect of a coal mine ground area. The method comprises the following steps: setting a measuring line and a detection frequency in a coal mine ground area to be detected, wherein the measuring line comprises a plurality of measuring points; the method comprises the steps of utilizing a ground-air frequency domain electromagnetic detection method to obtain ground-air frequency domain electromagnetic detection data of each measuring point in a coal mine ground area to be detected in different grouting time periods under detection frequency, wherein the grouting time periods comprise a pre-grouting time period, a grouting middle time period and a post-grouting time period, and the ground-air frequency domain electromagnetic detection data comprise voltage, current, frequency and waveforms; for each period, calculating apparent resistivity of the corresponding measuring point based on the earth-air frequency domain electromagnetic detection data of each measuring point in the period; and determining the grouting effect of the ground area of the coal mine to be detected according to the change rate of the apparent resistivity of each measuring point in different time periods. The application can solve the monitoring hysteresis of the grouting effect of the coal mine ground area in the prior art, and realize the dynamic monitoring of the grouting effect of the coal mine ground area.
Description
Technical Field
The application relates to the technical field of exploration and treatment of coal mine water damage ground areas, in particular to a method and a device for monitoring grouting effect of coal mine ground areas.
Background
The coal mine safety production is greatly threatened by water damage, and 3 serious water damage accidents and 7 serious water damage accidents occur in the coal mine in the whole country in the year 2020. In order to prevent water inrush accidents in coal mining, regional treatment techniques are commonly used in coal mining.
Because the floor grouting reinforcement or the drainage of the aquifer has the following defects, on one hand, the grouting effect cannot be evaluated correctly, and on the other hand, the underground water resource is damaged or wasted. At present, the grouting effect can be judged to a certain extent by detecting the grouting effect by means of the traditional drilling method, however, the drilling method has certain hysteresis, dynamic monitoring of the grouting effect cannot be realized, the treatment result of many mines is not accurately evaluated after the mines are treated in areas, the drilling method has higher cost and longer working period, the cost is not beneficial to saving, and the mining is efficient.
Therefore, there is a need for a method that can dynamically and accurately monitor grouting effects and save costs and efficiently produce.
Disclosure of Invention
The application provides a method and a device for monitoring grouting effect of a coal mine ground area, which are used for solving the problem of hysteresis in monitoring grouting effect of the coal mine ground area in the prior art.
In a first aspect, the application provides a method for monitoring grouting effect of a coal mine ground area, which comprises the following steps:
setting a measuring line and a detection frequency in a coal mine ground area to be detected, wherein the measuring line comprises a plurality of measuring points;
acquiring ground-air frequency domain electromagnetic detection data of each measuring point in the ground area of the coal mine to be detected in different grouting time periods under the detection frequency by using a ground-air frequency domain electromagnetic detection method, wherein the grouting time periods comprise a pre-grouting time period, a grouting middle time period and a post-grouting time period, and the ground-air frequency domain electromagnetic detection data comprise voltage, current, frequency and waveforms;
for each period, calculating apparent resistivity of the corresponding measuring point based on the earth-air frequency domain electromagnetic detection data of each measuring point in the period;
and determining the grouting effect of the coal mine ground area to be detected according to the change rate of apparent resistivity of each measuring point in different time periods.
In a second aspect, the application provides a device for monitoring grouting effect of a coal mine ground area, comprising:
the setting module is used for setting a measuring line and a detection frequency in a coal mine ground area to be detected, wherein the measuring line comprises a plurality of measuring points;
the acquisition module is used for acquiring ground-air frequency domain electromagnetic detection data of each measuring point in the ground area of the coal mine to be detected in different grouting time periods by using a ground-air frequency domain electromagnetic detection method, wherein the grouting time periods comprise a pre-grouting time period, a grouting middle time period and a post-grouting time period, and the ground-air frequency domain electromagnetic detection data comprise voltage, current, frequency and waveforms;
The calculation module is used for calculating the apparent resistivity of the corresponding measuring point according to the ground-air frequency domain electromagnetic detection data of each measuring point in each time period;
and the determining module is used for determining the grouting effect of the coal mine ground area to be detected according to the change rate of the apparent resistivity of each measuring point in different time periods.
In a third aspect, the present application provides a terminal comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor implementing the steps of the method according to the first aspect or any one of the possible implementations of the first aspect when the computer program is executed.
In a fourth aspect, the present application provides a computer readable storage medium storing a computer program which when executed by a processor implements the steps of the method of the first aspect or any one of the possible implementations of the first aspect.
The application provides a method and a device for monitoring grouting effect of a coal mine ground area, which are characterized in that ground-air frequency domain electromagnetic detection data of each measuring point in the coal mine ground area to be detected in a period before grouting, a period in grouting and a period after grouting are respectively obtained under the same detection frequency by using a ground-air frequency domain electromagnetic detection method, apparent resistivity of the corresponding measuring point is calculated according to the ground-air frequency domain electromagnetic detection data of each measuring point in each period, and grouting effect of the coal mine ground area to be detected is determined according to the change rate of apparent resistivity of each measuring point in the period before grouting, the period in grouting and the period after grouting. The application can solve the monitoring hysteresis of the grouting effect of the coal mine ground area in the prior art by using the earth-air frequency domain electromagnetic detection method, realizes the dynamic monitoring of the grouting effect of the coal mine ground area, and can avoid the monitoring of the grouting effect of the coal mine ground area by using the drilling method, thereby having low cost, saving cost and improving the exploitation efficiency.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are needed in the embodiments or the description of the prior art will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present application, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of an implementation of a method for monitoring grouting effect of a coal mine ground area, which is provided by an embodiment of the application;
FIG. 2 is a schematic distribution diagram of apparent resistivity during post-grouting time intervals according to an embodiment of the present application;
FIG. 3 is a schematic diagram of the distribution of resistivity after diffusion of grouting formation slurry during a post-grouting period provided by an embodiment of the present application;
FIG. 4 is a schematic distribution diagram of apparent resistivity during a pre-grouting period according to an embodiment of the present application;
FIG. 5 is a schematic diagram showing the distribution of apparent resistivity over time during grouting according to an embodiment of the present application;
fig. 6 is a schematic structural diagram of a monitoring device for grouting effect in a coal mine ground area according to an embodiment of the application;
fig. 7 is a schematic diagram of a terminal according to an embodiment of the present application.
Detailed Description
In the following description, for purposes of explanation and not limitation, specific details are set forth such as the particular system architecture, techniques, etc., in order to provide a thorough understanding of the embodiments of the present application. It will be apparent, however, to one skilled in the art that the present application may be practiced in other embodiments that depart from these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present application with unnecessary detail.
For the purpose of making the objects, technical solutions and advantages of the present application more apparent, the following description will be made by way of specific embodiments with reference to the accompanying drawings.
Fig. 1 is a flowchart of an implementation of a method for monitoring grouting effect in a coal mine ground area according to an embodiment of the present application, which is described in detail below:
in step 101, a measuring line and a detecting frequency are set for a coal mine ground area to be detected, wherein the measuring line comprises a plurality of measuring points.
In the embodiment of the application, the coal mine ground area to be detected on the coal mine ground is determined by detection purposes according to the detection scheme design. The actual topography survey is carried out on the coal mine ground area to be detected, the coordinates of the mine square are converted into WGS84 coordinates, and the GPS navigation system is used for positioning. And correcting the GPS navigation system by using the known coordinate datum points before positioning, selecting an optimal survey line scheme, and determining the electrode pit position. Setting a survey line according to the optimal survey line scheme, and arranging the survey line according to the optimal survey line scheme. The network measurement density is arranged according to practical conditions, and minimum values are set for line distances among the measuring lines and point distances of a plurality of measuring points in the measuring lines, for example, the minimum values are 10 meters, namely, the point distances are 10 meters, and the line distances are 10 meters. The purpose of setting up the survey line is accurate finding out the exploration scope inflection point on the austenitic map.
Illustratively, the GPS coordinates of a known coordinate reference point, such as a certain point coordinate X, are queried in the mine: 4203613.6049, y:19668524.1705, the corresponding GPS coordinates are X:37.9489251604, y:112.9177971243, at least two points are needed to determine accuracy. If the GPS coordinates of the accurate coordinate reference points are not found on the mine, the site detection is needed.
It should be noted that the line is required to be parallel to the line between the AB poles of the emission source, the included angle is not more than 15 °, and the line is selected for the purpose of determining the route and the emission source.
Because the general observation frequency of the earth space frequency domain electromagnetic detection method is in the frequency range of 1Hz to 1000Hz, in the embodiment of the application, the full frequency domain detection is firstly carried out, the data processing is carried out, the imaging is carried out, the analysis of the drilling data is compared, and the detection frequency domain which is most suitable for the actual situation is determined. And setting corresponding detection frequency for the coal mine ground area to be detected according to the determined detection frequency domain which is most suitable for reality. The embodiment of the application adopts the same detection frequency, so that the detection precision can be effectively improved, the data quantity generated by multiple detection frequencies can be reduced, the interpretation speed can be improved, the influence of other physical conditions can be reduced, and the contrast effect can be enhanced.
When the detection is performed before the treatment of the coal mine ground area to be detected, namely, the detection is performed in a period before grouting, the detection frequency set in the coal mine ground area to be detected is adjusted according to the grouting quantity of the corresponding layer of the known area treatment, and only the detection frequency is selected, so that the data quantity and the workload of data processing are reduced, and the detection efficiency is improved.
In one possible implementation, after the survey line is set in the coal mine ground area to be tested, a flight route is formulated.
Specific: chargeable equipment such as a power battery, a ground station, a receiving station and the like is fully charged, an overhaul machine maintains various equipment, records are made, experimental flight is carried out before detection, a flight plan is made, a flight route and an airplane landing point are designed, and the to-be-flown route is led into the ground station of the unmanned aerial vehicle for use by the unmanned aerial vehicle. And according to the flying route, the unmanned flying speed, flying height, transmitting frequency of the transmitter, transmitting current, receiving and transmitting distance and gain of the receiver are set.
Wherein, in the embodiment of the application, the survey line is substantially identical to the flying route. The difference is that: the lines are all from one direction to another, for example, line L1 is from north to south 1200 meters, and line L2 is from north to south 1200 meters. The flight direction and the survey line are not necessarily the same, and the survey line flies in opposite directions, for example, the survey line L1 direction is from north to south by 1200 meters, and the survey line L2 direction is from south to north by 1200 meters. The purpose of this setting is to facilitate mission setting for a flight frame, which may include several flight lines, so that it is necessary to separately design the survey line and the flight lines.
In step 102, using an earth space frequency domain electromagnetic detection method, obtaining earth space frequency domain electromagnetic detection data of each measuring point in a coal mine ground area to be detected in different grouting time periods under detection frequency, wherein the grouting time periods comprise a pre-grouting time period, a grouting middle time period and a grouting post-time period, and the earth space frequency domain electromagnetic detection data comprise voltage, current, frequency and waveforms.
The ground-air frequency domain electromagnetic detection method (Semi-Airborne Electromagnetic Method, SAEM) adopts a ground emission mode, and an air measurement response magnetic field working mode, and has the advantages of wide detection range, strong adaptability to complex terrains, low detection cost and high detection efficiency. The ground-air frequency domain electromagnetic detection method has the advantages of large detection, is sensitive to low-low resistance bodies, is not limited by topography, can realize dynamic monitoring before, during and after grouting, realizes evaluation of slurry diffusion and grouting effect, and has very important practical significance for preventing water burst accidents in coal mines.
In the embodiment of the application, by using an earth-air frequency domain electromagnetic detection method, under the detection frequency set in step 101, earth-air frequency domain electromagnetic detection data of each measuring point in a coal mine ground area to be detected in a period before grouting, a period in grouting and a period after grouting are obtained, wherein the earth-air frequency domain electromagnetic detection data comprises transmitting end data and receiving end data, the transmitting end data comprises voltage, current, frequency and different waveforms, and the receiving end data comprises a GPS data file and recorded flight track data which can be received.
The method comprises the steps of detecting before regional treatment, acquiring the electromagnetic detection data of the ground-air frequency domain of each measuring point in the ground region of the coal mine to be detected in the period before grouting under the detection frequency determined in the period before grouting, and knowing the water-rich condition of the stratum in the period before grouting.
And (3) detecting in the regional treatment, and carrying out second detection on the coal mine ground region to be detected under the same detection frequency to obtain ground-air frequency domain electromagnetic detection data of the grouting period.
And detecting after the region is treated, and detecting the ground region of the coal mine to be detected for the third time through detecting after the treatment of the ground-air electromagnetic region, so as to obtain ground-air frequency domain electromagnetic detection data of the grouting period.
Because only a unique drilling method is used for checking grouting effect in the past grouting engineering, and the application of the ground-to-air frequency domain electromagnetic detection method can realize the detection with no damage, continuous multiple stages and low cost on the basis of no drilling. In the hydraulic Cheng Ouyu control of coal mine prevention and treatment, the application adopts an earth-air frequency domain electromagnetic detection method to realize dynamic monitoring of grouting effect and auxiliary verification hole construction, the horizontal branch holes are generally at 60m intervals, and if the diffusion radius is large, the distance between the branch holes is increased, so that the grouting speed is improved. The diffusion radius is small, the distance between branch holes is reduced, and the safety of the coal mine is improved, so that the safe construction and efficient mining of the coal mine are realized.
In one possible implementation, after step 102, the method may further include:
judging whether the impedance phase corresponding to the ground-air frequency domain electromagnetic detection data is larger than a preset impedance phase or not;
if the impedance phase corresponding to the ground space frequency domain electromagnetic detection data is larger than the preset impedance phase, calculating the mean square relative error of the impedance phase;
if the impedance phase corresponding to the ground space frequency domain electromagnetic detection data is smaller than or equal to the preset impedance phase, calculating the mean square error of the impedance phase;
if the mean square relative error or the mean square error is not larger than the preset working error, determining that the acquired ground-to-air frequency domain electromagnetic detection data accords with the standard;
wherein, the mean square relative error calculation formula is:
wherein M is i The mean square relative error, n is the number of frequency points which are checked, observed and participated in statistics, m i The Carnikom resistivity of the ith frequency point;
the mean square error calculation formula is:
wherein ε i In the form of a mean square error,original observed impedance phase data for the ith frequency bin,/>Observing impedance phase data for the inspection of the ith frequency point;
accordingly, step 103 may include:
for each period, the apparent resistivity of the corresponding measuring point is calculated based on the electromagnetic detection data of the ground-air frequency domain, which accords with the standard, of each measuring point in the period.
According to the embodiment of the application, after the detection work is finished, the three-level detection secondary acceptance check is strictly executed on the detection and acceptance of the acquired ground-air frequency domain electromagnetic detection data, so that the high quality and reliability of the ground-air frequency domain electromagnetic detection data are ensured.
The ground-air electromagnetic field work quality inspection line is arranged, and accurate and reliable ground-air frequency domain electromagnetic detection data is guaranteed. A certain number of quality inspection routes are arranged in the working area and are basically and uniformly distributed in the whole working area and the whole construction process, so that the working quality of the whole area can be comprehensively reflected. The number of the inspection lines is greater than 5% of the total number of the lines; the observation of the quality check points follows the observation mode of 'two by two different', namely, the observation mode is carried out by different operators at different times by using the same instrument set in the same measuring line.
Specifically, it is determined whether the impedance phase corresponding to the electromagnetic detection data in the ground-air frequency domain obtained in step 102 is greater than a preset impedance phase, and if so, the mean square relative error corresponding to the impedance phase is calculated by using a mean square relative error calculation formula. If the impedance phase is not greater than the preset impedance phase, calculating the relative corresponding mean square error of the impedance by using a mean square error calculation formula. And then the quality of the electromagnetic detection data in the ground-to-air frequency domain is measured by utilizing the relative mean square error or the judgment result of the magnitude of the relative mean square error and the preset working error. Wherein, the preset impedance phase is set to be 200mrad according to the technical specification of the controllable source audio magnetotelluric method referenced by the data quality measurement method. The quality and the reliability of the ground-air frequency domain electromagnetic detection data are ensured.
When the quality check result of the same point shows one of the following situations, the quality disqualification of the earth-space frequency domain electromagnetic detection data of the point should be determined:
(1) The mean square relative error or the frequency point number of which the mean square error exceeds the preset working error is larger than 1/3 of the total number of the frequency points of the physical point;
(2) The mean square relative error or the frequency point number of which the mean square error exceeds 2 times of the preset working error is larger than 5% of the total number of the frequency points of the physical point;
(3) The mean square relative error or the observed value of the mean square error exceeding the preset working error continuously appears on 3 adjacent frequency points;
(4) The mean square relative error or mean square error is greater than the preset working error.
Correspondingly, after reliability and quality of the ground-air frequency domain electromagnetic detection data are judged and the obtained ground-air frequency domain electromagnetic detection data are determined to accord with the standard, the corresponding apparent resistivity is calculated based on the ground-air frequency domain electromagnetic detection data, which accord with the standard, of each measuring point in each period.
In step 103, for each period, apparent resistivity of the corresponding measurement point is calculated based on the earth-air frequency domain electromagnetic detection data of the respective measurement points in the period.
Wherein, the apparent resistivity is reflected by the formation deposition rule, and when the formation deposition distribution is uniform and is not damaged or changed, the apparent resistivity is uniformly changed in a layering way; if the influence of a plurality of factors such as sandstone gaps, formation water and the like is encountered, the apparent resistivity value is relatively smaller than that of a stratum with even sediment distribution.
In the embodiment of the application, the apparent resistivity of each measuring point in each period is calculated according to the earth-air frequency domain electromagnetic detection data acquired in the step 102.
In one possible implementation, before step 103, the method may further include:
sequentially carrying out denoising, spectrum analysis, calibration, normalization, posture correction and leveling treatment on the ground-to-air frequency domain electromagnetic detection data;
accordingly, step 103 may include:
for each period, the apparent resistivity of the corresponding measuring point is calculated based on the preprocessed ground-air frequency domain electromagnetic detection data of each measuring point in the period.
Specifically, because the earth-air frequency domain electromagnetic detection data is influenced by the change of the air coil gesture, a low-frequency geomagnetic field component and a horizontal magnetic field component generated by the same-frequency transmission signal are easily introduced, so that the quality of the earth-air frequency domain electromagnetic detection data is obviously reduced. The field has more obvious influence on the remote source receiving point, because the signal is attenuated rapidly along with the increase of the receiving distance, the signal is weak at the remote source, and the gesture noise can cause the obvious reduction of the signal-to-noise ratio of the data. Therefore, the application needs to perform preprocessing such as noise elimination and filtering on the ground-air frequency domain electromagnetic detection data before the data inversion imaging and after the acquired ground-air frequency domain electromagnetic detection data is determined to meet the standard. The method comprises the steps of removing a low-frequency baseline through wavelet denoising, filtering, superposing, performing spectrum analysis (FFT), calibrating and normalizing, correcting the posture, leveling and the like, and preprocessing ground-air frequency domain electromagnetic detection data meeting the standard.
The pretreatment process comprises the following steps:
(1) Denoising includes wavelet denoising to remove low frequency baselines, filtering, and superposition:
wavelet denoising removes low frequency baselines: aiming at low-frequency noise introduced by low-frequency geomagnetic field and coil low-frequency motion, a wavelet transformation method is adopted to remove a low-frequency baseline, so that the influence of baseline drift in an effective signal is reduced.
And (3) filtering: the out-of-band noise is filtered out by filtering the data before the spectrum analysis, so that the leakage site is prevented when the data is subjected to truncated data spectrum analysis.
And (3) superposition: and in the sampling time, multi-period superposition is carried out on the data, and non-whole period random noise in the signal is removed.
(2) Spectral analysis: and carrying out spectrum analysis operation on the processed time sequence to obtain the amplitude and the phase of each frequency point signal.
(3) Calibrating and normalizing: and (3) performing instrument calibration and magnetic sensor calibration on the data after the spectrum analysis to obtain the corresponding amplitude and phase of the measured magnetic field at each frequency point, and performing current normalization on the amplitude and phase to obtain a normalized magnetic field.
(4) Posture correction: and carrying out posture correction on the calibrated data, and removing the same-frequency posture noise in the data through the data and the correction factors.
(5) Leveling: and carrying out space moving average filtering on the data of different measuring lines, different measuring points and different frequencies. So as to realize the data is too flat and ensure the measurement consistency of adjacent measuring lines.
Correspondingly, for each period, the apparent resistivity of the corresponding measuring point is calculated based on the preprocessed earth-air frequency domain electromagnetic detection data of each measuring point in the period.
In one possible implementation, step 103 may include:
the apparent resistivity is calculated by a first formula:
wherein ρ is 1 For apparent resistivity, σ 1 For uniform earth conductivity, f is the target frequency point transmitting frequency, k 1 B is the wavenumber of the uniform ground z Mu, the induction field at Rx 0 Is vacuum magnetic permeability, I is current amplitude in an electric dipole source, L is wire length, lambda is integral variable, mu 1 Epsilon for uniform earth permeability 1 For uniform earth dielectric constant, z is the height of the receiving inductive coil placement, J 1 Is a first order bessel function of the first class.
Specifically, the earth-air frequency domain electromagnetic detection method is a method for obtaining the underground apparent resistivity distribution by measuring a single vertical magnetic field component in the air. In the uniform half-space model, the vertical magnetic field response at a certain point is shown in formula (1):
the apparent resistivity of each point can be obtained through magnetic field fitting iteration by referring to the concept of the apparent resistivity of the whole area in the wide area electromagnetic method and the apparent depth is defined by a skin depth formula.
Specific: when the power supply is arranged on the ground, the length of the wire is L, the current amplitude in the electric dipole source is I, and the receiving induction coil is placed in the air with the height of |z| (z)<0) Induced magnetic field B at Rx z Is shown in the formula (2):
wherein mu 0 ＝4π×10 -7 H/m，r TE For reflection coefficient related to formation resistivity and thickness +.>k 0 Is wave number in air, lambda is integral variable, J 1 Is a first order bessel function of the first class.
By measuring magnetic field signals at specific positions in the air under different frequencies, detection of resistivity structures with different depths in an underground space can be realized.
And (3) referring to the concept of the apparent resistivity of the whole region, obtaining the apparent resistivity of each point through magnetic field fitting iteration, wherein a specific calculation formula refers to a first formula.
In one possible implementation manner, the embodiment of the application also needs to edit and smooth the calculated apparent resistivity result and perform the processing calculation of the simulated section mapping. The method comprises the following steps:
(1) Editing and smoothing: editing and smoothing is a method for manually eliminating noise, and is to manually correct or eliminate special frequency points which are different from actual conditions, such as abrupt jump points, abnormal points which are generated in low-frequency conditions, and the like, in the original data on the basis of keeping the acquired original data so as to achieve the aim of conforming to the general change trend of resistivity.
(2) Spatial filtering: in order to smooth the curve, the curve is required to be subjected to spatial filtering smoothing treatment, and the embodiment of the application adopts a spatial filtering algorithm of '3 points are transversely taken and 3 points are longitudinally taken'.
(3) Drawing a simulated section: and carrying out boundary smoothing, noise elimination treatment and static effect correction on the obtained ground-air frequency domain electromagnetic detection data, and carrying out space filtering to obtain the resistivity pseudo sections with one-to-one correspondence of frequency and resistivity.
In step 104, the grouting effect of the coal mine ground area to be measured is determined according to the change rate of apparent resistivity of each measuring point in different time periods.
In the embodiment of the application, the grouting effect of the coal mine ground area to be detected is determined according to the change rate of apparent resistivity of each measuring point in the period before grouting, the period in grouting and the period after grouting. The grouting effect comprises the grouting effect meeting the standard requirement and not meeting the standard requirement.
In one possible implementation, step 104 may include:
judging whether the change rate of apparent resistivity of each measuring point is increased along with the transition of grouting time period or not;
if the change rate of apparent resistivity of each measuring point is increased along with the lapse of grouting time period, determining that the grouting effect of the ground area of the coal mine to be detected meets the standard requirement.
Specifically, for each measuring point, judging whether the change rate of apparent resistivity of the measuring point is increased along with the transition of the grouting period, if so, determining that the grouting effect of the ground area of the coal mine to be detected meets the standard requirement.
The earth-space frequency domain electromagnetic detection method is sensitive to underground barriers, low in cost and free of terrain limitation. The resistivity of the grouting material is higher than that of water, if the resistivity is high, the grouting effect is good, and if the resistivity is not proportional to the grouting amount, the grouting material can be leaked. The grouting is generally carried out by judging whether the grouting is finished or not through pressure, and is generally several megapascals, a pressurized water test is carried out after the grouting, if the total pressure is not less than 2-3 times of the maximum hydrostatic pressure of the injected aquifer, and the grouting can be carried out after the grouting is stabilized for more than 30 minutes, the hole can be sealed, and the grouting is considered to be finished by the drilling, and the grouting effect reaches the standard.
In one possible implementation, after determining whether the rate of change of apparent resistivity of each measurement point increases with the passage of the grouting period, the method may further include:
if a measuring point exists, the change rate of apparent resistivity of which does not increase along with the transition of the grouting period, determining that the grouting effect of the ground area of the coal mine to be detected does not meet the standard requirement.
Specifically, if the change rate of apparent resistivity of some measuring points does not increase along with the lapse of grouting time, determining that the grouting effect of the measuring points in the ground area of the coal mine to be detected does not reach the standard requirement.
In one possible implementation manner, after determining that the grouting effect of the coal mine ground area to be tested does not meet the standard requirement, the method may further include:
outputting position verification hole indication information of a target measuring point, wherein the target measuring point is a measuring point in which the change rate of apparent resistivity in all the measuring points does not increase along with the transition of grouting time period;
and acquiring the water-rich quantity of the verification hole at the position of the target measuring point by the user, and determining that the grouting quantity of the target measuring point is less or grouting slurry is lost according to the water-rich quantity.
Specifically, during detection after regional treatment, the corresponding apparent resistivity is obtained according to the earth-air frequency domain electromagnetic detection data of the post-grouting period, and the distribution schematic diagram of the apparent resistivity condition of the specific post-grouting period is shown in fig. 2. And comparing the apparent resistivity with the state of the period before grouting and the apparent resistivity value change of the grouting layer in the ground area treatment, if a relatively low-resistance area exists, judging that the area is a grouting weak area with small grouting amount or grouting slurry leakage, and outputting position verification hole indication information of a target measuring point through a display device. And guiding the construction of the verification hole according to the position verification hole indication information. And the drilling verification is carried out on the area with lower apparent resistivity, if the area is not rich in water, grouting slurry overflows, the grouting slurry theory is more than that of the area with lower apparent resistivity in the area with higher apparent resistivity, and the overall grouting effect of the area treatment of the area can be judged to reach the standard requirement. The distribution diagram of resistivity after grouting stratum slurry diffusion in a specific period after grouting in a coal mine ground area to be detected is shown in fig. 3.
In a possible implementation manner, the embodiment of the application can also obtain apparent resistivity according to the earth-space frequency domain electromagnetic detection data of the period before grouting when the region is detected before treatment, and can know the water-rich condition of the stratum in the period before grouting. A schematic distribution of apparent resistivity for a specific pre-grouting period is shown in fig. 4.
In a possible implementation manner, the embodiment of the application can also obtain the apparent resistivity equivalent curve section chart according to the earth air frequency domain electromagnetic detection data of the grouting middle period when the grouting is performed in the regional treatment, and can compare the change of the apparent resistivity curve because the resistivity of the regional treatment slurry is higher than that of the stratum in the period before grouting, thereby obtaining the diffusion condition of the slurry in the regional treatment, and judging whether the grouting effect of the regional treatment of the ground is achieved, further adjusting the interval of the horizontal radial holes, reducing the workload or improving the mine recovery safety. A schematic distribution of apparent resistivity for a specific post-grouting period is shown in fig. 5.
In a possible implementation manner, the embodiment of the application also utilizes the grouting amount in the regional treatment process to verify the attribute of the coal mine ground regional treatment layer to be detected, evaluates the data detected before regional treatment by combining the information of the known drilling holes, and obtains the attribute of the grouting layer in the coal mine ground region to be detected according to the evaluated data; according to the attribute of the grouting horizon of the coal mine ground area to be detected, calculating faults similar to the attribute of the fault in the recoverable working face in the fault of the coal mine ground area to be detected, performing linear fitting on the calculated faults, and predicting the faults in the adjacent working faces according to the roadway disclosure information and the attribute of the fault linearly fitted in the coal mine ground area to be detected. The application can carry out advanced remote prediction on the adjacent working surface, and is convenient for optimizing the reasonable arrangement of the adjacent working surface in advance.
The application determines the grouting effect condition of the coal mine ground area to be detected according to the detected resistivity distribution of the coal mine ground area to be detected based on the change relation between the change rate of apparent resistivity and the change relation between the time interval before grouting, the time interval during grouting and the time interval after grouting of the coal mine ground area to be detected, and carries out corresponding evaluation.
The application provides a monitoring method of grouting effect of coal mine ground area, which is characterized in that the earth-air frequency domain electromagnetic detection method is utilized to respectively obtain the earth-air frequency domain electromagnetic detection data of each measuring point in the coal mine ground area to be detected in the period before grouting, the period in grouting and the period after grouting under the same detection frequency, the apparent resistivity of the corresponding measuring point is calculated according to the earth-air frequency domain electromagnetic detection data of each measuring point in each period, and the grouting effect of the coal mine ground area to be detected is determined according to the change rate of the apparent resistivity of each measuring point in the period before grouting, the period in grouting and the period after grouting. The application can solve the monitoring hysteresis of the grouting effect of the coal mine ground area in the prior art by using the earth-air frequency domain electromagnetic detection method, realizes the dynamic monitoring of the grouting effect of the coal mine ground area, and can avoid the monitoring of the grouting effect of the coal mine ground area by using the drilling method, thereby having low cost, saving cost and improving the exploitation efficiency.
It should be understood that the sequence number of each step in the foregoing embodiment does not mean that the execution sequence of each process should be determined by the function and the internal logic, and should not limit the implementation process of the embodiment of the present application.
The following are device embodiments of the application, for details not described in detail therein, reference may be made to the corresponding method embodiments described above.
Fig. 6 is a schematic structural diagram of a monitoring device for grouting effect in a coal mine ground area according to an embodiment of the present application, and for convenience of explanation, only a portion relevant to the embodiment of the present application is shown, which is described in detail below:
as shown in fig. 6, the monitoring device 6 for grouting effect in the ground area of the coal mine comprises:
the setting module 61 is configured to set a measurement line and a detection frequency for a coal mine ground area to be detected, where the measurement line includes a plurality of measurement points;
the acquisition module 62 is configured to acquire, by using a ground-air frequency domain electromagnetic detection method, ground-air frequency domain electromagnetic detection data of each measurement point in a ground area of a coal mine to be detected at a detection frequency in different grouting periods, where the grouting periods include a period before grouting, a period during grouting and a period after grouting, and the ground-air frequency domain electromagnetic detection data includes voltage, current, frequency and waveform;
A calculation module 63, configured to calculate, for each period, apparent resistivity of a corresponding measurement point based on earth-air frequency domain electromagnetic detection data of each measurement point in the period;
the determining module 64 is configured to determine a grouting effect of the ground area of the coal mine to be measured according to the change rate of apparent resistivity of each measuring point in different time periods.
The application provides a monitoring device for grouting effect of coal mine ground area, which is characterized in that the earth-air frequency domain electromagnetic detection data of each measuring point in the coal mine ground area to be detected in the time period before grouting, the time period in grouting and the time period after grouting are respectively obtained under the same detection frequency by utilizing an earth-air frequency domain electromagnetic detection method, the apparent resistivity of the corresponding measuring point is calculated according to the earth-air frequency domain electromagnetic detection data of each measuring point in each time period, and the grouting effect of the coal mine ground area to be detected is determined according to the change rate of the apparent resistivity of each measuring point in the time period before grouting, the time period in grouting and the time period after grouting. The application can solve the monitoring hysteresis of the grouting effect of the coal mine ground area in the prior art by using the earth-air frequency domain electromagnetic detection method, realizes the dynamic monitoring of the grouting effect of the coal mine ground area, and can avoid the monitoring of the grouting effect of the coal mine ground area by using the drilling method, thereby having low cost, saving cost and improving the exploitation efficiency.
In one possible implementation, the computing module may be specifically configured to:
the apparent resistivity is calculated by a first formula:
wherein p is 1 For apparent resistivity, σ 1 For uniform earth conductivity, f is the target frequency point transmitting frequency, k 1 B is the wavenumber of the uniform ground z Mu, the induction field at Rx 0 Is vacuum magnetic permeability, I is current amplitude in an electric dipole source, L is wire length, lambda is integral variable, mu 1 Epsilon for uniform earth permeability 1 For uniform earth dielectric constant, z is the height of the receiving inductive coil placement, J 1 Is a first order bessel function of the first class.
In one possible implementation manner, the determining module may specifically include:
the judging module is used for judging whether the change rate of the apparent resistivity of each measuring point is increased along with the transition of the grouting period;
and the standard determining module is used for determining that the grouting effect of the ground area of the coal mine to be detected meets the standard requirement if the change rate of the apparent resistivity of each measuring point is increased along with the transition of the grouting period.
In one possible implementation manner, after the determining module, the apparatus may further include:
and the substandard determining module is used for determining that the grouting effect of the ground area of the coal mine to be detected does not reach the standard requirement if a measuring point exists, wherein the change rate of the apparent resistivity of the measuring point does not increase along with the transition of the grouting period.
In one possible implementation, after the failure to reach the standard determination module, the apparatus may be further configured to:
outputting position verification hole indication information of a target measuring point, wherein the target measuring point is a measuring point in which the change rate of apparent resistivity in all the measuring points does not increase along with the transition of grouting time period;
and acquiring the water-rich quantity of the verification hole at the position of the target measuring point by the user, and determining that the grouting quantity of the target measuring point is less or grouting slurry is lost according to the water-rich quantity.
In one possible implementation, after the acquisition module, the apparatus may be further configured to:
judging whether the impedance phase corresponding to the ground-air frequency domain electromagnetic detection data is larger than a preset impedance phase or not;
if the impedance phase corresponding to the ground space frequency domain electromagnetic detection data is larger than the preset impedance phase, calculating the mean square relative error of the impedance phase;
if the impedance phase corresponding to the ground space frequency domain electromagnetic detection data is smaller than or equal to the preset impedance phase, calculating the mean square error of the impedance phase;
if the mean square relative error or the mean square error is not larger than the preset working error, determining that the acquired ground-to-air frequency domain electromagnetic detection data accords with the standard;
wherein, the mean square relative error calculation formula is:
wherein M is i The mean square relative error, n is the number of frequency points which are checked, observed and participated in statistics, m i The Carnikom resistivity of the ith frequency point;
the mean square error calculation formula is:
wherein ε i In the form of a mean square error,original observed impedance phase data for the ith frequency bin,/>Observing impedance phase data for the inspection of the ith frequency point;
accordingly, the computing module may include:
for each period, the apparent resistivity of the corresponding measuring point is calculated based on the electromagnetic detection data of the ground-air frequency domain, which accords with the standard, of each measuring point in the period.
In one possible implementation, before the computing module, the apparatus may further include:
the preprocessing module is used for sequentially carrying out denoising, spectrum analysis, calibration, normalization, attitude correction and leveling treatment on the ground-to-air frequency domain electromagnetic detection data;
correspondingly, the computing module specifically may include:
for each period, the apparent resistivity of the corresponding measuring point is calculated based on the preprocessed ground-air frequency domain electromagnetic detection data of each measuring point in the period.
Fig. 7 is a schematic diagram of a terminal according to an embodiment of the present application. As shown in fig. 7, the terminal 7 of this embodiment includes: a processor 70, a memory 71, and a computer program 72 stored in the memory 71 and executable on the processor 70. The processor 70, when executing the computer program 72, performs the steps of the above-described embodiments of the method for monitoring grouting effects in the ground area of each coal mine, for example, steps 101 to 104 shown in fig. 1. Alternatively, the processor 70 may implement the functions of the modules in the above-described apparatus embodiments, such as the functions of the modules shown in fig. 6, when executing the computer program 72.
By way of example, the computer program 72 may be partitioned into one or more modules that are stored in the memory 71 and executed by the processor 70 to complete the present application. The one or more modules/units may be a series of computer program instruction segments capable of performing the specified functions, which instruction segments are used for describing the execution of the computer program 72 in the terminal 7. For example, the computer program 72 may be partitioned into the modules shown in FIG. 6.
The terminal 7 may be a computing device such as a desktop computer, a notebook computer, a palm computer, a cloud server, etc. The terminal 7 may include, but is not limited to, a processor 70, a memory 71. It will be appreciated by those skilled in the art that fig. 7 is merely an example of the terminal 7 and is not limiting of the terminal 7, and may include more or fewer components than shown, or may combine some components, or different components, e.g., the terminal may further include input and output devices, network access devices, buses, etc.
The processor 70 may be a central processing unit (Central Processing Unit, CPU), or may be another general purpose processor, a digital signal processor (Digital Signal Processor, DSP), an application specific integrated circuit (Application Specific Integrated Circuit, ASIC), a Field-programmable gate array (Field-Programmable Gate Array, FPGA) or other programmable logic device, a discrete gate or transistor logic device, a discrete hardware component, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory 71 may be an internal storage unit of the terminal 7, such as a hard disk or a memory of the terminal 7. The memory 71 may be an external storage device of the terminal 7, such as a plug-in hard disk, smart Media Card (SMC), a Secure Digital (SD Card, flash Card) or the like, which are provided on the terminal 7. Further, the memory 71 may also include both an internal storage unit and an external storage device of the terminal 7. The memory 71 is used for storing the computer program as well as other programs and data required by the terminal. The memory 71 may also be used for temporarily storing data that has been output or is to be output.
It will be apparent to those skilled in the art that, for convenience and brevity of description, only the above-described division of the functional units and modules is illustrated, and in practical application, the above-described functional distribution may be performed by different functional units and modules according to needs, i.e. the internal structure of the apparatus is divided into different functional units or modules to perform all or part of the above-described functions. The functional units and modules in the embodiment may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit, where the integrated units may be implemented in a form of hardware or a form of a software functional unit. In addition, the specific names of the functional units and modules are only for distinguishing from each other, and are not used for limiting the protection scope of the present application. The specific working process of the units and modules in the above system may refer to the corresponding process in the foregoing method embodiment, which is not described herein again.
In the foregoing embodiments, the descriptions of the embodiments are emphasized, and in part, not described or illustrated in any particular embodiment, reference is made to the related descriptions of other embodiments.
Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus/terminal and method may be implemented in other manners. For example, the apparatus/terminal embodiments described above are merely illustrative, e.g., the division of the modules or units is merely a logical function division, and there may be additional divisions when actually implemented, e.g., multiple units or components may be combined or integrated into another system, or some features may be omitted or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed may be an indirect coupling or communication connection via interfaces, devices or units, which may be in electrical, mechanical or other forms.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in the embodiments of the present application may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units may be implemented in hardware or in software functional units.
The integrated modules/units, if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the present application may implement all or part of the flow of the method of the above embodiment, or may be implemented by instructing related hardware by a computer program, where the computer program may be stored in a computer readable storage medium, and the computer program may implement the steps of the method embodiment for monitoring grouting effects of each coal mine ground area when executed by a processor. Wherein the computer program comprises computer program code which may be in source code form, object code form, executable file or some intermediate form etc. The computer readable medium may include: any entity or device capable of carrying the computer program code, a recording medium, a U disk, a removable hard disk, a magnetic disk, an optical disk, a computer Memory, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), an electrical carrier signal, a telecommunications signal, a software distribution medium, and so forth. It should be noted that the computer readable medium may include content that is subject to appropriate increases and decreases as required by jurisdictions in which such content is subject to legislation and patent practice, such as in certain jurisdictions in which such content is not included as electrical carrier signals and telecommunication signals.
The above embodiments are only for illustrating the technical solution of the present application, and not for limiting the same; although the application has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present application, and are intended to be included in the scope of the present application.
Claims (10)
1. The method for monitoring the grouting effect of the coal mine ground area is characterized by comprising the following steps of:
setting a measuring line and a detection frequency in a coal mine ground area to be detected, wherein the measuring line comprises a plurality of measuring points;
acquiring ground-air frequency domain electromagnetic detection data of each measuring point in the ground area of the coal mine to be detected in different grouting time periods under the detection frequency by using a ground-air frequency domain electromagnetic detection method, wherein the grouting time periods comprise a pre-grouting time period, a grouting middle time period and a post-grouting time period, and the ground-air frequency domain electromagnetic detection data comprise voltage, current, frequency and waveforms;
for each period, calculating apparent resistivity of the corresponding measuring point based on the earth-air frequency domain electromagnetic detection data of each measuring point in the period;
And determining the grouting effect of the coal mine ground area to be detected according to the change rate of apparent resistivity of each measuring point in different time periods.
2. The method for monitoring grouting effect of coal mine ground area according to claim 1, wherein for each period, calculating apparent resistivity of the corresponding measuring point based on the ground-air frequency domain electromagnetic detection data of each measuring point in the period comprises:
calculating the apparent resistivity by a first formula:
wherein ρ is 1 Sigma for the apparent resistivity 1 For uniform earth conductivity, f is the target frequency point transmitting frequency, k 1 B is the wavenumber of the uniform ground z Mu, the induction field at Rx 0 Is vacuum magnetic conductivity, I is current amplitude in the electric dipole source, L is wire length, lambda is productFractional variable, mu 1 Epsilon for uniform earth permeability 1 For uniform earth dielectric constant, z is the height of the receiving inductive coil placement, J 1 Is a first order bessel function of the first class.
3. The method for monitoring grouting effect of coal mine ground area according to claim 1, wherein the determining grouting effect of the coal mine ground area according to the change rate of apparent resistivity of each measuring point in different time periods comprises:
Judging whether the change rate of apparent resistivity of each measuring point is increased along with the transition of grouting time period or not;
if the change rate of apparent resistivity of each measuring point is increased along with the lapse of grouting time period, determining that the grouting effect of the coal mine ground area to be detected meets the standard requirement.
4. A method of monitoring grouting effects in a coal mine floor area as claimed in claim 3, wherein after said separately determining whether the rate of change of apparent resistivity at each station increases with the passage of grouting time, the method further comprises:
if a measuring point exists, the change rate of apparent resistivity of which does not increase along with the transition of the grouting period, determining that the grouting effect of the coal mine ground area to be detected does not meet the standard requirement.
5. The method for monitoring grouting effects of a coal mine floor area according to claim 4, wherein after the determining that the grouting effects of the coal mine floor area to be tested do not meet the standard requirements, the method further comprises:
outputting position verification hole indication information of a target measuring point, wherein the target measuring point is a measuring point in which the change rate of apparent resistivity in all measuring points does not increase along with the transition of grouting time period;
And acquiring the water-rich quantity of the verification hole at the position of the target measuring point by a user, and determining that the grouting quantity of the target measuring point is small or grouting slurry is lost according to the water-rich quantity.
6. The method for monitoring grouting effect of coal mine ground area according to claim 1, wherein after the method for utilizing the ground-air frequency domain electromagnetic detection method obtains the ground-air frequency domain electromagnetic detection data of each measuring point in the coal mine ground area to be detected in different grouting periods at the detection frequency, the method further comprises:
judging whether the impedance phase corresponding to the ground-air frequency domain electromagnetic detection data is larger than a preset impedance phase or not;
if the impedance phase corresponding to the ground-air frequency domain electromagnetic detection data is larger than the preset impedance phase, calculating the mean square relative error of the impedance phase;
if the impedance phase corresponding to the ground-air frequency domain electromagnetic detection data is smaller than or equal to the preset impedance phase, calculating the mean square error of the impedance phase;
if the mean square relative error or the mean square error is not larger than a preset working error, determining that the acquired ground-air frequency domain electromagnetic detection data meets a standard;
wherein, the mean square relative error calculation formula is:
Wherein M is i For the mean square relative error, n is the number of frequency points which are checked, observed and participated in statistics, m i The Carnikom resistivity of the ith frequency point;
the mean square error calculation formula is as follows:
wherein ε i For the mean square error of the signal to be processed,is the firstOriginal observed impedance phase data of i frequency points, < >>Observing impedance phase data for the inspection of the ith frequency point;
correspondingly, for each period, calculating the apparent resistivity of the corresponding measuring point based on the earth-air frequency domain electromagnetic detection data of each measuring point in the period comprises the following steps:
for each period, the apparent resistivity of the corresponding measuring point is calculated based on the electromagnetic detection data of the ground-air frequency domain, which accords with the standard, of each measuring point in the period.
7. The method for monitoring grouting effects of a coal mine floor area according to claim 1, wherein before the calculating, for each period, the apparent resistivity of the corresponding measuring point based on the earth-air frequency domain electromagnetic detection data of the respective measuring point in the period, the method further comprises:
sequentially carrying out denoising, spectrum analysis, calibration, normalization, attitude correction and leveling treatment on the ground-air frequency domain electromagnetic detection data;
correspondingly, for each period, calculating the apparent resistivity of the corresponding measuring point based on the earth-air frequency domain electromagnetic detection data of each measuring point in the period comprises the following steps:
For each period, the apparent resistivity of the corresponding measuring point is calculated based on the preprocessed ground-air frequency domain electromagnetic detection data of each measuring point in the period.
8. Monitoring device of colliery ground area slip casting effect, its characterized in that includes:
the setting module is used for setting a measuring line and a detection frequency in a coal mine ground area to be detected, wherein the measuring line comprises a plurality of measuring points;
the acquisition module is used for acquiring ground-air frequency domain electromagnetic detection data of each measuring point in the ground area of the coal mine to be detected in different grouting time periods by using a ground-air frequency domain electromagnetic detection method, wherein the grouting time periods comprise a pre-grouting time period, a grouting middle time period and a post-grouting time period, and the ground-air frequency domain electromagnetic detection data comprise voltage, current, frequency and waveforms;
the calculation module is used for calculating the apparent resistivity of the corresponding measuring point according to the ground-air frequency domain electromagnetic detection data of each measuring point in each time period;
and the determining module is used for determining the grouting effect of the coal mine ground area to be detected according to the change rate of the apparent resistivity of each measuring point in different time periods.
9. A terminal comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the computer program, performs the steps of the method of monitoring grouting effects in a coal mine floor area as claimed in any one of claims 1 to 7.
10. A computer readable storage medium storing a computer program, wherein the computer program when executed by a processor performs the steps of the method for monitoring grouting effects of a coal mine floor area as claimed in any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311179475.7A CN117052471A (en) | 2023-09-13 | 2023-09-13 | Monitoring method and device for grouting effect of coal mine ground area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311179475.7A CN117052471A (en) | 2023-09-13 | 2023-09-13 | Monitoring method and device for grouting effect of coal mine ground area |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117052471A true CN117052471A (en) | 2023-11-14 |
Family
ID=88666382
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311179475.7A Pending CN117052471A (en) | 2023-09-13 | 2023-09-13 | Monitoring method and device for grouting effect of coal mine ground area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117052471A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117556198A (en) * | 2024-01-11 | 2024-02-13 | 中国地质大学（北京） | Apparent resistivity denoising calculation method based on full-waveform well ground method |
-
2023
- 2023-09-13 CN CN202311179475.7A patent/CN117052471A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117556198A (en) * | 2024-01-11 | 2024-02-13 | 中国地质大学（北京） | Apparent resistivity denoising calculation method based on full-waveform well ground method |
| CN117556198B (en) * | 2024-01-11 | 2024-04-02 | 中国地质大学（北京） | Apparent resistivity denoising calculation method based on full-waveform well ground method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Abd Allah et al. | Three-dimensional resistivity characterization of a coastal area: Application of Grounded Electrical-Source Airborne Transient Electromagnetic (GREATEM) survey data from Kujukuri Beach, Japan | |
| CN104656157B (en) | A kind of method and device in identification shale gas dessert area | |
| CN105116453A (en) | Transient electromagnetic exploration method and device for natural gas hydrate in permafrost zone | |
| CN107085240B (en) | Slope magnetofluid detection system and method | |
| CN102590874B (en) | Method for detecting ground surface crack of upland coal-mining subsidence paddy field | |
| Allah et al. | Three-dimensional resistivity modelling of grounded electrical-source airborne transient electromagnetic (GREATEM) survey data from the Nojima Fault, Awaji Island, south-east Japan | |
| CN103995301A (en) | Method and device for evaluating total organic carbon content in shale gas reservoir | |
| CN117052471A (en) | Monitoring method and device for grouting effect of coal mine ground area | |
| CN111045091B (en) | Recognition and positioning method for ancient river channel covered by basalt | |
| CN104714254A (en) | Coal mine water-accumulating gob detection method and device based on transient electromagnetic method | |
| Weymer et al. | Review of electromagnetic induction for mapping barrier island framework geology | |
| CN103176216A (en) | Pipeline detection method and borehole antenna | |
| Keating | Improved use of the local wavenumber in potential-field interpretation | |
| CN111735494A (en) | Method for monitoring permeation increasing process of low-permeability polluted site | |
| CN105220694A (en) | A kind of engineering jackstone amount evaluation method | |
| Reid et al. | An overview of the SkyTEM airborne EM system with Australian examples | |
| Lane et al. | Filling in the gaps? validation and integration of airborne EM data with surface and subsurface observations for catchment management? an example from Bendigo, Victoria, Australia | |
| RU2363965C1 (en) | Method designed to monitor local irregularities and geodynamic zones of geological section top part (gst) | |
| CN104749645B (en) | A kind of method and apparatus of coal seam substrate interface undulation form detection in ash top difficult to understand | |
| RU2678535C1 (en) | Method of monitoring state of diaphragm from boron-cutting alumina-cement concrete piles in earth dam by electrical tomography | |
| CN112462432A (en) | High-precision detection method and device for coal field fire area cavity | |
| Shan et al. | Portable audio magnetotellurics-experimental measurements and joint inversion with radiomagnetotelluric data from Gotland, Sweden | |
| Magawata et al. | Geo-electric assessment of kali failed dam project aliero, North Western Nigeria | |
| CN114488346B (en) | Underground space abnormal body detection method, device, equipment and medium | |
| CN115182362B (en) | Side slope time-sharing partition accurate targeting protection method and device, electronic equipment and medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |