A detection method, device and electronic equipment for the water-rich situation of coal and rock strata in a roadway

Abstract

The present invention provides a method, apparatus and electronic device for detecting the water-rich situation of roadway coal and rock strata. The method includes: transmitting a transient electromagnetic signal to the roadway coal and rock strata; receiving the transient electromagnetic response reflected by the roadway coal and rock strata; correcting the transient electromagnetic response based on a preset correction coefficient to obtain a corrected transient electromagnetic response; the correction coefficient is calculated based on the transient electromagnetic response under the condition of no interference and the transient electromagnetic response under the condition of interference; determining the water-rich situation of the roadway coal and rock strata based on the corrected transient electromagnetic response. The present invention can improve the accuracy of detecting the water-rich situation of roadway coal and rock strata with metals present.

Description

Technical Field

[0001] The present invention relates to the technical field of geological exploration, and particularly to a method, device and electronic equipment for detecting the water-rich situation of coal and rock strata in roadways. Background Art

[0002] Hidden water-conducting structures in front of and on the sidewalls of roadway headings are the main hidden dangers of mine water inrush, which will cause serious harm to the safe production of coal mines. Therefore, it is necessary to find out the water-rich situation of coal and rock strata within a certain range in front of the roadway heading in advance during the roadway driving process.

[0003] The transient electromagnetic method is a commonly used method for detecting the water-rich situation of coal and rock strata in roadways. The transient electromagnetic method can detect water-conducting structures with different positions and forms in the surrounding space in coal mine roadways. With the advantages of flexible construction, large detection range, small volume effect, strong directionality, high resolution, and sensitivity to low-resistance areas, the transient electromagnetic method has become the main method for coal mine water hazard detection.

[0004] During the actual detection process, metal is one of the main factors affecting the detection accuracy of the transient electromagnetic method. Due to geological reasons and production needs, there are many metal objects installed in coal mine roadways, such as metal mesh sheets, bolts, anchor cables for support, and I-beams in the shed support section; tracks and belt racks for transportation; cables and switches for power supply; roadheaders for production, etc. These metal objects will have a certain impact on the transient electromagnetic detection results. Compared with geological bodies, metals can generate induction electromagnetic fields with greater intensity, which will significantly enhance the intensity of the induced electric field received by the coil. Moreover, the closer the metal is, the larger the metal size, the greater the induced electromotive force, the smaller the apparent resistivity, and the smaller the detection depth, resulting in a significant decrease in the measured apparent resistivity compared to the roadway environment without metal interference, and the induced electromotive force attenuation curve is smoother. During the actual underground transient electromagnetic detection process, metals will not only absorb and attenuate the primary field emitted, but also absorb and attenuate the secondary field again during the return process after turning off the emission current. Since the resistivity of metals is relatively low, the attenuation rate of induced current in them is much slower than that in high-resistance media, which will also cause the early response curve to change slowly. Therefore, when there are metals in the coal and rock strata of the roadway, there are interference problems in the transient electromagnetic method, resulting in inaccurate detection of the water-rich situation of the coal and rock strata in the roadway. Summary of the Invention

[0005] The present invention provides a method, device and electronic equipment for detecting the water-rich situation of coal and rock strata in roadways, which can improve the accuracy of detecting the water-rich situation of coal and rock strata in roadways with metals.

[0006] In a first aspect, the present invention provides a method for detecting the water-rich condition of roadway coal and rock strata, including: transmitting a transient electromagnetic signal to the roadway coal and rock strata; receiving the transient electromagnetic response reflected by the roadway coal and rock strata; correcting the transient electromagnetic response based on a preset correction coefficient to obtain a corrected transient electromagnetic response; the correction coefficient is calculated based on the transient electromagnetic response under the condition of no interference and the transient electromagnetic response under the condition of interference; determining the water-rich condition of the roadway coal and rock strata based on the corrected transient electromagnetic response.

[0007] The present invention provides a method for detecting the water-rich condition of roadway coal and rock strata. By correcting the transient electromagnetic response reflected by the roadway coal and rock strata and determining the water-rich condition of the roadway coal and rock strata based on the corrected transient electromagnetic response. Since the correction coefficient is calculated based on the transient electromagnetic response under the condition of no interference and the transient electromagnetic response under the condition of interference, the water-rich condition determined based on the corrected transient electromagnetic response is more in line with the water-rich condition under the condition of no interference, reducing the interference of metal on the transient electromagnetic method and improving the accuracy of detecting the water-rich condition of roadway coal and rock strata with metal.

[0008] In a possible implementation manner, determining the water-rich condition of the roadway coal and rock strata based on the corrected transient electromagnetic response includes: determining the apparent resistivity of the roadway coal and rock strata based on the corrected transient electromagnetic response; determining the water-rich condition of the roadway coal and rock strata based on the apparent resistivity.

[0009] In a possible implementation manner, determining the water-rich condition of the roadway coal and rock strata based on the apparent resistivity includes: if the apparent resistivity is greater than the set resistivity, it is determined that the water content of the roadway coal and rock strata is less; if the apparent resistivity is less than or equal to the set resistivity, it is determined that the water content of the roadway coal and rock strata is greater.

[0010] In a possible implementation manner, correcting the transient electromagnetic response based on a preset correction coefficient to obtain a corrected transient electromagnetic response includes: determining the corrected transient electromagnetic response based on the following formula;

[0011] D C = k * D U ;

[0012] where k is the preset correction coefficient, D C is the corrected transient electromagnetic response, D0 is the simulated response under the condition of no interference, D R is the simulated response under the condition of interference, D U is the transient electromagnetic response.

[0013] In a possible implementation manner, before correcting the transient electromagnetic response based on a preset correction coefficient to obtain the corrected transient electromagnetic response, it further includes: based on the transient electromagnetic response, searching a correction coefficient database to determine an interference response that matches the transient electromagnetic response; the interference response is a transient electromagnetic response obtained by simulating the situation where there is interference with the roadway coal rock stratum; determining the correction coefficient corresponding to the interference response as the preset correction coefficient.

[0014] In a possible implementation manner, before correcting the transient electromagnetic response based on a preset correction coefficient to obtain the corrected transient electromagnetic response, it further includes: Step 1: Based on the geological characteristics of the roadway coal rock stratum, construct a non-interference roadway coal rock stratum model and an interference roadway coal rock stratum model; Step 2: Perform finite volume division on the non-interference roadway coal rock stratum model and divide it into M volume elements to obtain a non-interference model; Step 3: Perform finite volume division on the interference roadway coal rock stratum model and divide it into M volume elements to obtain an interference model; where the M volume elements of the non-interference model correspond one by one to the M volume elements of the interference model; Step 4: Simulate the emission source to emit transient electromagnetic signals to the non-interference model and the interference model; simulate to obtain the non-interference response of the non-interference model and the interference response of the interference model; Step 5: Determine the correction coefficient based on the non-interference response and the interference response; Step 6: Change the position of the emission source, the geological characteristics of the roadway coal rock stratum, and the interference information in the interference roadway coal rock stratum model, and repeat Steps 1 to 4 to obtain a correction coefficient database; where the correction coefficient database includes multiple correction coefficients and the interference response corresponding to each correction coefficient.

[0015] In a possible implementation manner, determining the correction coefficient based on the non-interference response and the interference response includes: based on the response data at each time point in the interference response, draw an interference response curve; fit the curve slope in the starting time period of the interference curve to obtain a starting slope; starting from the starting moment, if the error between the slope at a certain moment in the interference response curve and the starting slope time is greater than a set value, then determine this moment as the anchor point moment; determine the correction coefficient based on the data of the non-interference response and the interference response at each moment before the anchor point moment.

[0016] In a second aspect, an embodiment of the present invention provides a detection device for the water-rich situation of a roadway coal rock stratum, which is characterized in that it includes: a communication module, configured to emit a transient electromagnetic signal to the roadway coal rock stratum; receive the transient electromagnetic response reflected by the roadway coal rock stratum; a processing module, configured to correct the transient electromagnetic response based on a preset correction coefficient to obtain a corrected transient electromagnetic response; the correction coefficient is calculated based on the transient electromagnetic response in a non-interference situation and the transient electromagnetic response in an interference situation; determine the water-rich situation of the roadway coal rock stratum based on the corrected transient electromagnetic response.

[0017] In a possible implementation, the processing module is specifically configured to determine the apparent resistivity of the roadway coal and rock strata based on the corrected transient electromagnetic response, and determine the water-rich condition of the roadway coal and rock strata based on the apparent resistivity.

[0018] In a possible implementation, the processing module is specifically configured to determine that the water content in the roadway coal and rock strata is low if the apparent resistivity is greater than the set resistivity, and determine that the water content in the roadway coal and rock strata is high if the apparent resistivity is less than or equal to the set resistivity.

[0019] In a possible implementation, the processing module is specifically configured to determine the corrected transient electromagnetic response based on the following formula:

[0020] D C = k * D U ;

[0021] where k is a preset correction coefficient, D C is the corrected transient electromagnetic response, D0 is the simulated response under no interference, D R is the simulated response under interference, and D U is the transient electromagnetic response.

[0022] In a possible implementation, the processing module is further configured to search the correction coefficient database based on the transient electromagnetic response to determine the interference response matching the transient electromagnetic response; the interference response is the transient electromagnetic response obtained by simulating the situation where there is interference in the roadway coal and rock strata; and determine the correction coefficient corresponding to the interference response as the preset correction coefficient.

[0023] In a possible implementation, the processing module is further configured to perform the following steps: Step 1: Based on the geological characteristics of the roadway coal and rock strata, construct an undisturbed roadway coal and rock strata model and a disturbed roadway coal and rock strata model; Step 2: Perform finite volume division on the undisturbed roadway coal and rock strata model, and divide it into M volume elements to obtain an undisturbed model; Step 3: Perform finite volume division on the disturbed roadway coal and rock strata model, and divide it into M volume elements to obtain a disturbed model; where the M volume elements of the undisturbed model correspond one-to-one to the M volume elements of the disturbed model; Step 4: Simulate the emission source to emit transient electromagnetic signals to the undisturbed model and the disturbed model; simulate and obtain the undisturbed response of the undisturbed model and the disturbed response of the disturbed model; Step 5: Based on the undisturbed response and the disturbed response, determine the correction coefficient; Step 6: Change the position of the emission source, the geological characteristics of the roadway coal and rock strata, and the interference information in the disturbed roadway coal and rock strata model, and repeat Steps 1 to 4 to obtain a correction coefficient database; where the correction coefficient database includes multiple correction coefficients and the disturbed response corresponding to each correction coefficient.

[0024] In a possible implementation, the processing module is specifically configured to draw a disturbed response curve based on the response data at each time point in the disturbed response; fit the curve slope in the starting time period of the disturbed curve to obtain the starting slope; starting from the starting moment, if the error between the slope of a certain moment in the disturbed response curve and the starting slope time is greater than the set value, then determine that moment as the anchor point moment; based on the undisturbed response data and the disturbed response data at each moment before the anchor point moment, determine the correction coefficient.

[0025] In a third aspect, an embodiment of the present invention provides an electronic device, characterized in that the electronic device includes a memory and a processor, the memory stores a computer program, and the processor is configured to call and run the computer program stored in the memory to execute the steps of the method as described in the first aspect and any possible implementation manner in the first aspect.

[0026] In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium, the computer-readable storage medium stores a computer program, characterized in that when the computer program is executed by a processor, it implements the steps of the method as described in the first aspect and any possible implementation manner in the first aspect.

[0027] The technical effects brought by any implementation manner in the second to fourth aspects above can refer to the technical effects brought by the corresponding implementation manner in the first aspect, which will not be elaborated here. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] To more clearly illustrate the technical solutions in the embodiments of the present invention, the following will briefly introduce the accompanying drawings required for the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0029] Figure 1 is a schematic flowchart of a method for detecting the water-rich situation of roadway coal and rock strata provided by an embodiment of the present invention;

[0030] Figure 2 is a schematic diagram of a transient electromagnetic response provided by an embodiment of the present invention;

[0031] Figure 3 is a schematic diagram of a transient electromagnetic response before and after correction provided by an embodiment of the present invention;

[0032] Figure 4 is a schematic structural diagram of a device for detecting the water-rich situation of roadway coal and rock strata provided by an embodiment of the present invention;

[0033] Figure 5 is a schematic structural diagram of an electronic device provided by an embodiment of the present invention. Detailed implementation manners

[0034] In the following description, for the purpose of illustration rather than limitation, specific details such as specific system structures and technologies are put forward in order to thoroughly understand the embodiments of the present invention. However, those skilled in the art should clearly understand that the present invention can also be implemented in other embodiments without these specific details. In other cases, the detailed descriptions of well-known systems, devices, circuits, and methods are omitted to avoid unnecessary details from interfering with the description of the present invention.

[0035] In the description of the present invention, unless otherwise specified, " / " means "or". For example, A / B can represent A or B. The "and / or" herein is only a description of the association relationship of the associated objects, indicating that three relationships can exist. For example, A and / or B can represent: A exists alone, A and B exist simultaneously, and B exists alone. In addition, "at least one" and "multiple" mean two or more. The words such as "first" and "second" do not limit the quantity and execution order, and the words such as "first" and "second" do not necessarily limit to be different.

[0036] In the embodiments of the present application, words such as "exemplary" or "for example" are used to represent examples, illustrations, or explanations. Any embodiment or design solution described as "exemplary" or "for example" in the embodiments of the present application should not be construed as being more preferred or having more advantages than other embodiments or design solutions. Rather, the use of words such as "exemplary" or "for example" is intended to present relevant concepts in a specific manner for easy understanding.

[0037] In addition, the terms "comprising" and "having" and any variations thereof mentioned in the description of the present application are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or device that includes a series of steps or modules is not limited to the listed steps or modules, but may optionally further include other unlisted steps or modules, or may optionally further include other steps or modules inherent to these processes, methods, products, or devices.

[0038] To make the objectives, technical solutions, and advantages of the present invention clearer, the following will be described through specific embodiments in conjunction with the drawings of the present invention.

[0039] As described in the background art, currently when there is metal in the roadway coal and rock strata, there are interference problems in the transient electromagnetic method, resulting in inaccurate detection of the water-rich situation in the roadway coal and rock strata.

[0040] To solve the above technical problems, as Figure 1 shown, an embodiment of the present invention provides a method for detecting the water-rich situation in the roadway coal and rock strata. The execution subject is a detection device for the water-rich situation in the roadway coal and rock strata, and this detection method includes steps S101 - S104.

[0041] S101. Transmit a transient electromagnetic signal to the roadway coal and rock strata.

[0042] As a possible implementation, the detection device may transmit a transient electromagnetic signal to the roadway coal and rock strata at a set measuring point.

[0043] S102. Receive the transient electromagnetic response reflected by the roadway coal and rock strata.

[0044] S103. Based on a preset correction coefficient, correct the transient electromagnetic response to obtain a corrected transient electromagnetic response.

[0045] In the embodiments of the present application, the correction coefficient is calculated based on the transient electromagnetic response under the condition of no interference and the transient electromagnetic response under the condition of interference.

[0046] As a possible implementation, the detection device may determine the corrected transient electromagnetic response based on the following formula.

[0047] D C = k * D U ;

[0048] where k is a pre-set correction coefficient. D C is the corrected transient electromagnetic response, D0 is the simulated response without interference, and D R is the simulated response under interference, and D U is the transient electromagnetic response.

[0049] Optionally, the correction coefficient can be a pre-set fixed value. Alternatively, the detection device can also determine a correction coefficient that matches the received transient electromagnetic response based on the received transient electromagnetic response.

[0050] As a possible implementation, before step S103, the detection device can also search a correction coefficient database based on the transient electromagnetic response to determine an interference response that matches the transient electromagnetic response; and determine the correction coefficient corresponding to the interference response as the pre-set correction coefficient. Among them, the interference response is the transient electromagnetic response obtained by simulating the situation where there is interference to the roadway coal and rock strata.

[0051] S104. Determine the water-rich situation of the roadway coal and rock strata based on the corrected transient electromagnetic response.

[0052] As a possible implementation, the detection device can determine the apparent resistivity of the roadway coal and rock strata based on the corrected transient electromagnetic response; and determine the water-rich situation of the roadway coal and rock strata based on the apparent resistivity.

[0053] Exemplarily, if the apparent resistivity is greater than the set resistivity, the detection device can determine that the water content in the roadway coal and rock strata is less.

[0054] Another exemplarily, if the apparent resistivity is less than or equal to the set resistivity, the detection device can determine that the water content in the roadway coal and rock strata is greater.

[0055] The present invention provides a method for detecting the water-rich situation of roadway coal and rock strata. By correcting the transient electromagnetic response reflected by the roadway coal and rock strata and determining the water-rich situation of the roadway coal and rock strata based on the corrected transient electromagnetic response. Since the correction coefficient is calculated based on the transient electromagnetic response without interference and the transient electromagnetic response under interference, the water-rich situation determined based on the corrected transient electromagnetic response is more in line with the water-rich situation without interference, reducing the interference of metal on the transient electromagnetic method and improving the accuracy of detecting the water-rich situation of roadway coal and rock strata with metal.

[0056] Optionally, for the method for detecting the water-rich situation of roadway coal and rock strata provided in the embodiments of the present invention, before step S103, it further includes steps one to six.

[0057] Step 1: Based on the geological characteristics of the roadway coal and rock strata, construct an undisturbed roadway coal and rock strata model and a disturbed roadway coal and rock strata model.

[0058] Step 2: Conduct a finite volume division on the undisturbed roadway coal and rock strata model, and divide it into M volume elements to obtain an undisturbed model.

[0059] Step 3: Conduct a finite volume division on the disturbed roadway coal and rock strata model, and divide it into M volume elements to obtain a disturbed model.

[0060] Among them, the M volume elements of the undisturbed model correspond one by one to the M volume elements of the disturbed model.

[0061] Step 4: Simulate the emission of transient electromagnetic signals from the emission source to the undisturbed model and the disturbed model, and simulate to obtain the undisturbed response of the undisturbed model and the disturbed response of the disturbed model.

[0062] It should be noted that the embodiments of the present application adopt a two-dimensional transient electromagnetic forward calculation method for simulation, use finite volume difference simulation, divide non-uniform grids, process the emission source as an analytical solution and add it as an initial condition at t>0, and use the Dirichlet boundary condition as the boundary condition.

[0063] In some embodiments, the finite volume difference simulation method of the transient electromagnetic response can use the Maxwell's equations under quasi-static conditions as the control equation for the induced electromagnetic field in the time domain. The Maxwell's equations under quasi-static conditions are as follows.

[0064]

[0065]

[0066]

[0067]

[0068] b = μh

[0069] j = σe;

[0070] Among them, e is the electric field strength, with the unit of V / m, b is the magnetic flux density, with the unit of T, h is the magnetic field strength, with the unit of A / m, j is the current density flux, with the unit of A / m2, s m is the magnetic source, s e is the electric source, σ is the conductivity, with the unit of S / m, μ is the magnetic permeability, with the unit of H / m, C is a discrete curl operator that maps the variables on the edges to the faces, C edges→faces ; is the inner product matrix of the face element variables, and maps the reciprocal of the magnetic permeability μ -1Average from the grid center to the cell faces is the inner product matrix of the variables on the edges, which averages the conductivity σ from the grid center to the edges s m and s e are the magnetic flux or electric field components obtained by numerical integration, defined on the cell faces and edges respectively

[0071] It should be noted that is Faraday's law, which describes how a time-varying magnetic flux generates a rotating electric field is Ampere's law, which describes how an electric current generates a rotating magnetic field is the constitutive relation, which describes how the field and the flux are converted

[0072] In numerical simulation, a discrete grid needs to be established and the governing equations need to be discretized. The discrete method used in the embodiments of the present application is the finite volume method, which divides the model into a tensor grid and defines physical quantities at different positions of the grid cells. Specifically, physical properties (σ, μ) and scalar potential are defined at the cell center vector fields (e, h) are defined on the cell edges; vector fluxes (b, j) are defined on the cell faces. Based on this, discrete differential operators are constructed, including curl gradient divergence and inner product averaging operators, including averaging the variables at the cell center to the 6 cell faces averaging the variables at the cell center to the 12 cell edges. The finite volume method on the staggered grid discretizes and solves the transient electromagnetic problem, and two sets of unknowns need to be considered. One is the field defined on the cell edges, and the other is the flux defined on the grid cell faces. For the transient electromagnetic inverse problem in the time domain, the electric field e is discretized on the cell edges, the magnetic flux density b is discretized on the cell faces, the physical property parameters conductivity σ and permeability μ are discretized at the cell center, and the external current source s e is discretized on the cell edges

[0073] In some embodiments, for the discretization in the time domain, the detection device can discretize the Maxwell equations in the time domain based on the following Euler difference formula to obtain the discrete form of the Maxwell equations

[0074]

[0075]

[0076] b(0) = b0

[0077] Based on the Maxwell's equations under quasi-static conditions and the discrete forms of the Maxwell's equations, the discrete forms of the Maxwell's equations at each moment can be obtained as follows.

[0078]

[0079]

[0080] where, Δt k = t k+1 - t k is the step size, and the superscripts k, k + 1 are time indices. The matrices of all time steps are assembled together to form a lower diagonal matrix. The matrix on the main diagonal is represented by and the matrix on the sub-diagonal is represented by . The lower diagonal matrix is as follows.

[0081]

[0082] It should be noted that the initial condition u0 of the Maxwell's equations depends on the type of the emission source and the emission waveform, and can be calculated by analytical or numerical methods. For example, if the emission source is a loop source and the emission waveform is a step current, the magnetic flux b0 can be obtained by calculating the analytical solution of the loop source in the full space and used as the initial value. If there is no analytical solution for the reflection waveform, set u0 = 0 and numerically obtain the waveform value of s m at subsequent moments, obtain the electric field on the edge at each moment, and then calculate other physical fields (h, j, b).

[0083] It should be noted that in the selection of boundary conditions, a sufficiently large grid space is adopted. At positions far from the emission source and the anomaly, the secondary field decays greatly. Therefore, the Dirichlet boundary condition is adopted, and the tangential component of the electric field and the normal component of the magnetic field on the boundary are set to zero.

[0084] For the electromagnetic field diffusion equation, the initial values of the electric field and the magnetic field need to be given, that is, the electric field value at t = t0, and the magnetic field value at. A magnetic dipole source in a homogeneous full space is used as the emission source, and the initial transient electromagnetic field is used as the emission source. When the current of the magnetic dipole source in the full space is disconnected, the generated electromagnetic field can be expressed by the following formula.

[0085]

[0086]

[0087]

[0088]

[0089]

[0090]

[0091] Among them, M = Idl is the emission magnetic moment, I is the emission current, and dl is the length of the current source.

[0092] Step Five: Determine the correction coefficient based on the interference-free response and the interference-present response.

[0093] As a possible implementation, the detection device can determine the correction coefficient based on Steps A1 - A4.

[0094] A1. Based on the response data at each time point in the interference-present response, plot the interference-present response curve.

[0095] A2. Fit the curve slope in the starting time period of the interference-present curve to obtain the starting slope.

[0096] A3. Starting from the starting moment, if the error between the slope at a certain moment in the interference-present response curve and the starting slope time is greater than the set value, then determine this moment as the anchor point moment.

[0097] A4. Based on the data of the interference-free response and the interference-present response at each moment before the anchor point moment, determine the correction coefficient.

[0098] It should be noted that for the correction preprocessing of transient electromagnetic measured data, there are various metal bodies in the underground environment, which will cause a gap between the inverted results and the actual situation. Therefore, for the unprocessed actual data, we need to perform correction processing first to correct the horizontal differences caused by human factors or environmental factors.

[0099] By establishing a model without considering low-resistance interferences such as rails and support, and a model of a roadway containing low-resistance bodies such as rails, support, anchor nets, etc. The responses without interference and with interference are simulated. The feature is that when the transient electromagnetic signal induced by metal interference can generate a stronger induced electromagnetic field relative to the geological body, and the closer the metal is, the larger its size, the greater the induced electromotive force, the smaller the apparent resistivity, and the smaller the detection depth. In the forward simulation, the uniform background conductivity value is set, and the background response participates in the correction. By adjusting the background conductivity parameters in the interference-free model and the interference-present model, the problem of a large reduction in the resistivity of the inverted data caused by metal interference is eliminated.

[0100] In this correction process, two forward simulations are mainly used to obtain the interference-free and interference-present data. If the interference-present data is close to the actual data, the interference-present data D R , and the data obtained by simulating the situation without interference is D0. Assume that the uncorrected actual data is D U , and the corrected data is D C, the expression for data correction is: D C =(D0 / D R )D U .

[0101] As Figure 2 shown, the first curve 201 is a simulated response with a high degree of coincidence between the actual data in the early stage and the late stage, and the second curve 202 is the transient electromagnetic response actually collected. When the transient electromagnetic response actually collected has a high degree of coincidence with the simulated response in the early stage and the late stage, it can be considered that the simulated response can be used for correction; if the coincidence between the simulated data with interference and the actual data is good, the correction coefficient can be constructed using this data and the data without interference, and the interference in the early stage of the actual data can be corrected. Then find a time node when the data change trends of each measuring point are basically the same; use the data of a certain measuring point at this time node as the anchor point, obtain the correction coefficient from the data at this time point, and then use this correction coefficient to correct the responses of the measuring points on each survey line. Or the data collected by the mine transient electromagnetic instrument from the underground roadway site, with the ordinate being the induced electromotive force dB / dt (V / m2) and the abscissa being the time t (s), select a piece of data with a disturbing body as the reference for correction, obtain the correction coefficient and apply it to the correction of the actual data.

[0102] As Figure 3 shown, a schematic diagram of the transient electromagnetic response before and after correction provided by the embodiment of the present application. The second curve 202 is the transient electromagnetic response actually collected, and the third curve 301 is the transient electromagnetic response after correction. By comparing the transient electromagnetic response after correction with that before correction, it can be seen that the characteristic of the early response curve being relatively flat due to metal low-resistance interference has been corrected, and the attenuation characteristic of the entire response curve is consistent with the response curve without interference. At the same time, the low-resistance response characteristic in the transient electromagnetic response before correction has not disappeared due to correction, and the low-resistance abnormal response characteristic in Figure 1 is still well preserved.

[0103] It should be noted that a uniform background conductivity parameter is set in the constructed model. By adjusting the background conductivity parameter, the problems of a large reduction in the resistivity of the inversion data and a reduction in the detection distance caused by metal interference can be eliminated. Applying the correction method to the inversion processing of actual data can make the obtained resistivity profile more consistent with the actual formation resistivity, and the position of the low-resistance abnormal body has a good correspondence with the low-resistance response characteristic in the actual data.

[0104] Step Six: Change the position of the emission source, the geological characteristics of the roadway coal seam, and the interference information in the roadway coal seam model with interference, and repeat Steps One to Four to obtain a correction coefficient database; where the correction coefficient database includes multiple correction coefficients and the interference responses corresponding to each correction coefficient.

[0105] In this way, the embodiments of the present application can adopt finite volume difference simulation to forward calculate the transient electromagnetic response with interference and the transient electromagnetic response without interference, thereby determining the correction coefficient, realizing the correction of the transient electromagnetic response, reducing the interference of metal on the transient electromagnetic method, and improving the accuracy of detecting the water-rich situation of the roadway coal and rock strata with metal.

[0106] In addition, the embodiments of the present application can also change the position of the emission source, the geological characteristics of the roadway coal and rock strata, and the interference information in the roadway coal and rock strata model with interference to determine the correction coefficients under different working conditions, making the correction coefficients more accurate, thereby further improving the accuracy of detecting the water-rich situation of the roadway coal and rock strata with metal.

[0107] It should be understood that the magnitudes of the sequence numbers of the steps in the above embodiments do not mean the order of execution. The order of execution of each process should be determined according to its function and internal logic, and should not constitute any limitation to the implementation process of the embodiments of the present invention.

[0108] The following is an apparatus embodiment of the present invention. For details not described in detail, reference can be made to the corresponding method embodiments above.

[0109] Figure 4 The structural schematic diagram of a detection apparatus for the water-rich situation of a roadway coal and rock strata provided by an embodiment of the present invention is shown. The detection apparatus 400 includes a communication module 401 and a processing module 402.

[0110] The communication module 401 is configured to transmit a transient electromagnetic signal to the roadway coal and rock strata; receive the transient electromagnetic response reflected by the roadway coal and rock strata.

[0111] The processing module 402 is configured to correct the transient electromagnetic response based on a preset correction coefficient to obtain a corrected transient electromagnetic response; the correction coefficient is calculated based on the transient electromagnetic response without interference and the transient electromagnetic response with interference; determine the water-rich situation of the roadway coal and rock strata based on the corrected transient electromagnetic response.

[0112] In a possible implementation manner, the processing module 402 is specifically configured to determine the apparent resistivity of the roadway coal and rock strata based on the corrected transient electromagnetic response; determine the water-rich situation of the roadway coal and rock strata based on the apparent resistivity.

[0113] In a possible implementation manner, the processing module 402 is specifically configured to determine that the water content in the roadway coal and rock strata is less if the apparent resistivity is greater than the set resistivity; determine that the water content in the roadway coal and rock strata is greater if the apparent resistivity is less than or equal to the set resistivity.

[0114] In a possible implementation, the processing module 402 is specifically configured to determine the corrected transient electromagnetic response based on the following formula:

[0115] D C = k * D U ;

[0116] where k is a preset correction coefficient, D C is the corrected transient electromagnetic response, D0 is the simulated response under no interference, D R is the simulated response under interference, and D U is the transient electromagnetic response.

[0117] In a possible implementation, the processing module 402 is further configured to, based on the transient electromagnetic response, search the correction coefficient database to determine the interference response that matches the transient electromagnetic response; the interference response is the transient electromagnetic response obtained by simulating the case where there is interference with the roadway coal and rock strata; and determine the correction coefficient corresponding to the interference response as the preset correction coefficient.

[0118] In a possible implementation, the processing module 402 is further configured to perform the following steps: Step 1: Based on the geological characteristics of the roadway coal and rock strata, construct a non-interference roadway coal and rock strata model and an interference roadway coal and rock strata model; Step 2: Perform finite volume division on the non-interference roadway coal and rock strata model and divide it into M volume elements to obtain a non-interference model; Step 3: Perform finite volume division on the interference roadway coal and rock strata model and divide it into M volume elements to obtain an interference model; where the M volume elements of the non-interference model correspond one-to-one with the M volume elements of the interference model; Step 4: Simulate the transient electromagnetic signal emitted by the emission source to the non-interference model and the interference model; simulate and obtain the non-interference response of the non-interference model and the interference response of the interference model; Step 5: Determine the correction coefficient based on the non-interference response and the interference response; Step 6: Change the position of the emission source, the geological characteristics of the roadway coal and rock strata, and the interference information in the interference roadway coal and rock strata model, and repeat Steps 1 to 4 to obtain a correction coefficient database; where the correction coefficient database includes multiple correction coefficients and the interference response corresponding to each correction coefficient.

[0119] In a possible implementation, the processing module 402 is specifically configured to, based on the response data at each time point in the interference response, draw an interference response curve; fit the curve slope in the starting time period of the interference curve to obtain a starting slope; starting from the starting moment, if the error between the slope at a certain moment in the interference response curve and the starting slope time is greater than a set value, then determine that moment as the anchor point moment; and determine the correction coefficient based on the data of the non-interference response and the interference response at each moment before the anchor point moment.

[0120] Figure 5 is a schematic structural diagram of an electronic device provided by an embodiment of the present invention. As Figure 5 shown, the electronic device 500 of this embodiment includes: a processor 501, a memory 502, and a computer program 503 stored in the memory 502 and executable on the processor 501. When the processor 501 executes the computer program 503, the steps in the above method embodiments are implemented, for example Figure 1 the steps 101 to 104 shown. Alternatively, when the processor 501 executes the computer program 503, the functions of each module / unit in the above device embodiments are implemented. For example, Figure 4 the functions of the communication module 401 and the processing module 402 shown.

[0121] Exemplarily, the computer program 503 can be divided into one or more modules / units. The one or more modules / units are stored in the memory 502 and executed by the processor 501 to complete the present invention. The one or more modules / units can be a series of computer program instruction segments capable of performing specific functions, and these instruction segments are used to describe the execution process of the computer program 503 in the electronic device 500. For example, the computer program 503 can be divided into Figure 4 the communication module 401 and the processing module 402 shown.

[0122] The so-called processor 501 may be a central processing unit (CPU), or may also be other general-purpose processors, digital signal processors (DSPs), application specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.

[0123] The memory 502 may be an internal storage unit of the electronic device 500, such as a hard disk or memory of the electronic device 500. The memory 502 may also be an external storage device of the electronic device 500, such as a plug-in hard disk equipped on the electronic device 500, a Smart Media Card (SMC), a Secure Digital (SD) card, a Flash Card, etc. Further, the memory 502 may also include both the internal storage unit of the electronic device 500 and an external storage device. The memory 502 is used to store the computer program and other programs and data required by the terminal. The memory 502 may also be used to temporarily store data that has been output or is to be output.

[0124] Those skilled in the art can clearly understand that, for the convenience and brevity of description, only the above division of each functional unit and module is used as an example. In actual applications, the above functions can be assigned to different functional units and modules according to needs, that is, the internal structure of the device is divided into different functional units or modules to complete all or part of the functions described above. Each functional unit and module in the embodiments can be integrated into a processing unit, or each unit can exist physically alone, or two or more units can be integrated into one unit. The above integrated units can be implemented in the form of hardware or in the form of software functional units. In addition, the specific names of each functional unit and module are only for the convenience of mutual distinction and do not limit the protection scope of the present application. The specific working processes of the units and modules in the above system can refer to the corresponding processes in the foregoing method embodiments and will not be elaborated herein.

[0125] In the above embodiments, the descriptions of the various embodiments have their own emphases. For the parts not detailed or recorded in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0126] Those of ordinary skill in the art can realize that the units and algorithm steps of each example described in combination with the embodiments disclosed herein can be implemented by electronic hardware, or by a combination of computer software and electronic hardware. Whether these functions are executed in hardware or software depends on the specific application and design constraints of the technical solution. Professional technicians can use different methods to implement the described functions for each specific application, but such implementation should not be considered to exceed the scope of the present invention.

[0127] In the embodiments provided by the present invention, it should be understood that the disclosed device / terminal and method can be implemented in other ways. For example, the device / terminal embodiments described above are merely illustrative. For example, the division of the modules or units is only a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed couplings or direct couplings or communication connections between each other can be through some interfaces. The indirect couplings or communication connections of the devices or units can be in electrical, mechanical or other forms.

[0128] The units described as separate components may or may not be physically separated. The components displayed as units may or may not be physical units, that is, they can be located in one place or distributed to multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0129] In addition, in each embodiment of the present invention, the functional units can be integrated in a processing unit, or each unit can exist physically alone, or two or more units can be integrated in one unit. The above integrated units can be implemented in the form of hardware or in the form of software functional units.

[0130] If the integrated module / unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, to implement all or part of the processes in the above method embodiments of the present invention, it can also be completed by a computer program instructing relevant hardware. The computer program can be stored in a computer-readable storage medium. When the computer program is executed by a processor, the steps of the above method embodiments can be implemented. Among them, the computer program includes computer program code, and the computer program code can be in the form of source code, object code, executable file or some intermediate form. The computer-readable medium can include: any entity or device that can carry the computer program code, recording medium, USB flash drive, mobile hard disk, magnetic disk, optical disc, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signal, telecommunication signal, and software distribution medium, etc.

[0131] The above-described embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: they can still modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements on some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the various embodiments of the present invention, and should all be included within the protection scope of the present invention.

Claims

1. A detection method for the water-rich situation of coal and rock strata in a roadway, characterized in that, Including: Transmitting a transient electromagnetic signal to the coal and rock strata in the roadway; Receiving the transient electromagnetic response reflected by the coal and rock strata in the roadway; Based on a preset correction coefficient, correcting the transient electromagnetic response to obtain a corrected transient electromagnetic response; the correction coefficient is calculated based on the transient electromagnetic response under non-interference conditions and the transient electromagnetic response under interference conditions; Based on the corrected transient electromagnetic response, determining the water-rich condition of the coal and rock strata in the roadway; Before the step of correcting the transient electromagnetic response based on a preset correction coefficient to obtain a corrected transient electromagnetic response, it further includes: Step 1: Based on the geological characteristics of the coal and rock strata in the roadway, constructing a non-interference coal and rock strata model of the roadway and an interference coal and rock strata model of the roadway; Step 2: Performing finite volume division on the non-interference coal and rock strata model of the roadway, dividing it into M volume elements to obtain a non-interference model; Step 3: Performing finite volume division on the interference coal and rock strata model of the roadway, dividing it into M volume elements to obtain an interference model; where, the M volume elements of the non-interference model correspond one by one to the M volume elements of the interference model; Step 4: Simulating the transient electromagnetic signal transmitted by the emission source to the non-interference model and the interference model; simulating to obtain the non-interference response of the non-interference model and the interference response of the interference model; Step 5: Based on the non-interference response and the interference response, determining the correction coefficient; Step 6: Changing the position of the emission source, the geological characteristics of the coal and rock strata in the roadway, and the interference information in the interference coal and rock strata model of the roadway, repeating Steps 1 to 4 to obtain a correction coefficient database; where, the correction coefficient database includes multiple correction coefficients and the interference response corresponding to each correction coefficient; The determining the correction coefficient based on the non-interference response and the interference response includes: Based on the response data at each time point in the interference response, drawing an interference response curve; fitting the curve slope in the starting time period of the interference curve to obtain a starting slope; starting from the starting moment, if the error between the slope of a certain moment in the interference response curve and the starting slope is greater than a set value, then determining this moment as the anchor point moment; based on the data of the non-interference response and the interference response at each moment before the anchor point moment, determining the correction coefficient.

2. The detection method for the water-rich situation of roadway coal and rock strata according to claim 1, characterized in that The determining the water-rich condition of the coal and rock strata in the roadway based on the corrected transient electromagnetic response includes: Based on the corrected transient electromagnetic response, determining the apparent resistivity of the coal and rock strata in the roadway; Based on the apparent resistivity, determining the water-rich condition of the coal and rock strata in the roadway.

3. The detection method for the water-rich situation of roadway coal and rock strata according to claim 2, characterized in that, The determining the water-rich condition of the coal and rock strata in the roadway based on the apparent resistivity includes: If the apparent resistivity is greater than the set resistivity, it is determined that the water content in the coal and rock strata in the roadway is less; If the apparent resistivity is less than or equal to the set resistivity, it is determined that the water content in the coal and rock strata in the roadway is greater.

4. The detection method for the water-rich situation of roadway coal and rock strata according to claim 1, characterized in that The correcting the transient electromagnetic response based on a preset correction coefficient to obtain a corrected transient electromagnetic response includes: Based on the following formula, determining the corrected transient electromagnetic response; D C = k * D U ; where k is the preset correction coefficient, D C is the corrected transient electromagnetic response, D0 is the simulated response under no interference, and D R is the simulated response under interference, and D U is the transient electromagnetic response.

5. The detection method for the water-rich situation of roadway coal and rock strata according to claim 1, characterized in that, Before correcting the transient electromagnetic response based on a preset correction factor to obtain the corrected transient electromagnetic response, it further includes: Based on the transient electromagnetic response, searching a correction factor database to determine an interference response that matches the transient electromagnetic response; the interference response is a transient electromagnetic response obtained by simulating a situation where there is interference with the roadway coal and rock strata; Determining the correction factor corresponding to the interference response as the preset correction factor.

6. A detection device for the water-rich situation of roadway coal and rock strata, characterized in that, It includes: A communication module for transmitting a transient electromagnetic signal to the roadway coal and rock strata; receiving the transient electromagnetic response reflected by the roadway coal and rock strata; A processing module for correcting the transient electromagnetic response based on a preset correction factor to obtain the corrected transient electromagnetic response; the correction factor is calculated based on the transient electromagnetic response under the non-interference condition and the transient electromagnetic response under the interference condition; determining the water-rich situation of the roadway coal and rock strata based on the corrected transient electromagnetic response; The processing module is further configured to perform the following steps: Step 1: Based on the geological characteristics of the roadway coal and rock strata, construct a non-interference roadway coal and rock strata model and an interference roadway coal and rock strata model; Step 2: Perform finite volume division on the non-interference roadway coal and rock strata model, and divide it into M volume elements to obtain a non-interference model; Step 3: Perform finite volume division on the interference roadway coal and rock strata model, and divide it into M volume elements to obtain an interference model; where, the M volume elements of the non-interference model correspond one by one to the M volume elements of the interference model; Step 4: Simulate the emission source to emit a transient electromagnetic signal to the non-interference model and the interference model; simulate to obtain the non-interference response of the non-interference model and the interference response of the interference model; Step 5: Determine the correction factor based on the non-interference response and the interference response; Step 6: Change the position of the emission source, the geological characteristics of the roadway coal and rock strata, and the interference information in the interference roadway coal and rock strata model, and repeat Steps 1 to 4 to obtain a correction factor database; where, the correction factor database includes multiple correction factors and the interference response corresponding to each correction factor; The processing module is specifically configured to draw an interference response curve based on the response data at each time point in the interference response; fit the curve slope in the starting time period of the interference curve to obtain the starting slope; starting from the starting moment, if the error between the slope of a certain moment in the interference response curve and the starting slope is greater than the set value, then determine that moment as the anchor point moment; determine the correction factor based on the data of the non-interference response and the interference response at each moment before the anchor point moment.

7. An electronic device, characterized in that, The electronic device includes a memory and a processor. The memory stores a computer program, and the processor is configured to call and run the computer program stored in the memory to execute the steps of the method according to any one of claims 1 to 5.

8. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the steps of the method according to any one of claims 1 to 5 as above.

Images