Direct current method inversion-based apparent resistivity data processing method and system
By applying the apparent resistivity data processing method based on DC inversion in coal mines, the problem of difficulty in monitoring the resistivity field changes in the prior art is solved, and accurate identification and dynamic monitoring of geological medium changes during coal seam excavation and mining is achieved, which improves safety and reduces costs.
Patent Information
- Application Number
- CN202510126854.2
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-01-27
- Publication Date
- 2025-04-11
AI Technical Summary
The existing DC power monitoring methods lack effective inversion analysis methods in coal mines, and it is difficult to monitor the changes in the resistivity field in a timely manner, which increases the difficulty and risk of excavation and mining work.
A method and system for processing apparent resistivity data based on DC inversion is provided. By inversion of the four-dimensional electric field data during the excavation and mining of the target coal seam, normalization process is performed, normalization gradient and timing change rate are calculated, and the geological medium change boundary and timing change situation are determined.
It realizes accurate identification of broken areas and instable areas during excavation and mining, dynamically reflects the distribution of dangerous areas, reduces the testing workload and construction costs, and improves the safety during excavation and mining.
Smart Images
Figure CN119805593A_ABST
Abstract
Description
Technical Field
[0001] The invention relates to the technical field of mine safety assessment, and in particular to a method and system for processing apparent resistivity data based on direct current electrical inversion. Background Art
[0002] During the excavation and mining of the mine, the geological structure is complex, and the instability of rocks and coal seams may lead to dangerous events such as gas outbursts and rock collapse. Traditional geological exploration methods such as drilling and acoustic wave detection have limitations and cannot effectively monitor changes during the coal mining process. The direct current method has been used as an effective geological exploration method due to its high sensitivity and wide range of applications. It can measure the resistivity field of underground media, thereby reflecting the structure and changes of geological bodies.
[0003] The existing DC electrical monitoring method lacks effective inversion analysis means in coal mine work, making it difficult to monitor the changes in the resistivity field in a timely manner, increasing the difficulty and risk of tunneling and mining work. Summary of the invention
[0004] In order to solve the above technical problems existing in the prior art, an embodiment of the present invention provides a method and system for processing apparent resistivity data based on direct current inversion. The technical solution is as follows:
[0005] On the one hand, a method for processing apparent resistivity data based on direct current electrical inversion is provided, comprising: performing apparent resistivity cloud map inversion on four-dimensional electric field data during excavation and mining of a target coal seam to obtain inverted data; performing normalization processing on the inverted data to obtain normalized resistivity distribution at each position of the target coal seam; determining a normalized λ distribution cloud map of the normalized resistivity distribution based on a natural neighbor interpolation method; respectively calculating a normalized gradient and a time series change rate of the normalized λ distribution cloud map; and determining a geological medium change boundary and a time series change of the target coal seam based on the normalized gradient and the time series change amount.
[0006] Furthermore, the method also includes: obtaining four-dimensional electric field data during the excavation and mining of the target coal seam by electrical equipment.
[0007] Furthermore, the inverted data is normalized to obtain the normalized resistivity distribution of each position of the target coal seam, including:
[0008]
[0009] In the formula, ρ i is the apparent resistivity value of the i-th point in the inverted data; ρ min , maxare the minimum and maximum apparent resistivity values in the test area, respectively, and λ is the normalized resistivity value.
[0010] Further, calculating the normalized gradient of the normalized λ distribution cloud map includes:
[0011] λ G (x,y,z)=λ(x,y,z)Grad
[0012] In the formula, λ G is the normalized gradient value, λ is the normalized resistivity value, and (x, y, z) is the three-dimensional space coordinate.
[0013] Further, calculating the time series change rate of the normalized λ distribution cloud diagram includes:
[0014]
[0015] Where (x, y, z, t) is the four-dimensional space-time coordinate, λ is the normalized resistivity value, is the timing change rate.
[0016] Furthermore, based on the normalized gradient and the time series change amount, the geological medium change boundary and the time series change of the target coal seam are determined, including: determining the position where the normalized gradient exceeds a preset threshold as the geological medium change boundary of the target coal seam; if the time series change rate is less than 0, determining that the structure of the target coal seam is in a compaction process; if the time series change amount is greater than 0, determining that the structure of the target coal seam is in a slow destruction stage.
[0017] On the other hand, a system for processing apparent resistivity data based on direct current electrical inversion is provided, comprising: an inversion module, a normalization module, a first determination module, a calculation module and a second determination module; wherein the inversion module is used to perform apparent resistivity cloud map inversion on four-dimensional electric field data during excavation and mining of a target coal seam to obtain inverted data; the normalization module is used to perform normalization processing on the inverted data to obtain the normalized resistivity distribution of each position of the target coal seam; the first determination module is used to determine the normalized λ distribution cloud map of the normalized resistivity distribution based on the natural neighbor interpolation method; the calculation module is used to respectively calculate the normalized gradient and time series change rate of the normalized λ distribution cloud map; the second determination module is used to determine the geological medium change boundary and time series change of the target coal seam based on the normalized gradient and the time series change amount, respectively.
[0018] Furthermore, the second determination module is also used to: determine the position where the normalized gradient exceeds a preset threshold as the geological medium change boundary of the target coal seam; if the time series change rate is less than 0, determine that the structure of the target coal seam is in a compaction process; if the time series change amount is greater than 0, determine that the structure of the target coal seam is in a slow destruction stage.
[0019] On the other hand, an electronic device is provided, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the method provided in the embodiment of the present invention when executing the computer program.
[0020] On the other hand, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores computer instructions, and when the computer instructions are executed by a processor, the method provided in the embodiment of the present invention is implemented.
[0021] The present invention provides a method and system for processing apparent resistivity data based on direct current electrical inversion. By testing the apparent resistivity of coal seam excavation and mining areas, the normalized electric field distribution and the gradient distribution characteristics of the gradient in the excavation and mining process are determined, the broken area and the unstable area in the excavation and mining process are identified, and the potential dangerous area is determined, thereby reducing the testing workload and construction cost, and providing a scientific and accurate technical means for processing the apparent resistivity obtained by direct current electrical inversion. BRIEF DESCRIPTION OF THE DRAWINGS
[0022] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.
[0023] Figure 1 It is a flow chart of a method for processing apparent resistivity data based on direct current electrical inversion provided by an embodiment of the present invention;
[0024] Figure 2 This is an overall schematic diagram of a direct current method laying, installation and testing provided by an embodiment of the present invention;
[0025] Figure 3 It is a cross-sectional schematic diagram of a direct current method laying and installation test provided by an embodiment of the present invention;
[0026] Figure 4 It is a schematic diagram of an apparent resistivity data processing system based on direct current electrical inversion provided by an embodiment of the present invention.
[0027] Description of reference numerals:
[0028] 1. Rock stratum, 2. Coal seam area, 3. Mudstone, 4. Bottom extraction tunnel, 5. DC electrical instrument, 6. Electrode, 7. Survey line. DETAILED DESCRIPTION
[0029] The technical solution of the present invention is described below in conjunction with the accompanying drawings.
[0030] In the embodiments of the present invention, words such as "exemplarily" and "for example" are used to indicate examples, illustrations or explanations. Any embodiment or design described as "example" in the present invention should not be interpreted as being more preferred or more advantageous than other embodiments or designs. Specifically, the use of the word "example" is intended to present the concept in a specific way. In addition, in the embodiments of the present invention, the meaning expressed by "and / or" can be both, or it can be either of the two.
[0031] In order to make the technical problems, technical solutions and advantages to be solved by the present invention more clear, a detailed description will be given below with reference to the accompanying drawings and specific embodiments.
[0032] Embodiment 1
[0033] Figure 1 FIG. 1 is a flow chart of a method for processing apparent resistivity data based on direct current inversion according to an embodiment of the present invention. Figure 1 As shown, the method specifically comprises the following steps:
[0034] Step S102, performing apparent resistivity cloud map inversion on the four-dimensional electric field data during the excavation and mining of the target coal seam to obtain inverted data.
[0035] Specifically, the embodiment of the present invention obtains four-dimensional electric field data of the target coal seam during excavation and mining through electrical equipment. The inverted data includes apparent resistivity data under four-dimensional changes (x, y, z, t); x, y, z represent three-dimensional spatial coordinates, and t represents time information.
[0036] The inverted data is then exported and the required data type is set to a .dat format file with five columns: x, y, z, t, and R; where R represents the apparent resistivity in four-dimensional space coordinates.
[0037] Step S104, normalizing the inverted data to obtain the normalized resistivity distribution of each position of the target coal seam.
[0038] Specifically, the normalization calculation formula includes:
[0039]
[0040] In the formula, ρ i is the apparent resistivity value of the i-th point in the inverted data; ρmin , max are the minimum and maximum apparent resistivity values in the test area, respectively, and λ is the normalized resistivity value.
[0041] Step S106, determining a normalized λ distribution cloud diagram of the normalized resistivity distribution based on the natural neighbor interpolation method.
[0042] Step S108, respectively calculating the normalized gradient and the time series change rate of the normalized λ distribution cloud map.
[0043] Specifically, the normalized gradient of the normalized λ distribution cloud map is calculated, including:
[0044] λ G (x,y,z)=λ(x,y,z)Grad
[0045] In the formula, λ G is the normalized gradient value, (x, y, z) is the three-dimensional space coordinate.
[0046] Specifically, the time series change rate of the normalized lambda distribution cloud graph is calculated, including:
[0047]
[0048] Where (x, y, z, t) is the four-dimensional space-time coordinate, is the timing rate of change.
[0049] Step S110, determining the geological medium change boundary and the time series change of the target coal seam based on the normalized gradient and the time series change amount respectively.
[0050] Specifically, step S110 also includes the following steps:
[0051] Step S1101, the position where the normalized gradient exceeds the preset threshold is determined as the geological medium change boundary of the target coal seam. For example, the preset threshold is 0.2, that is, the normalized gradient of the geological medium change interface of the target coal seam is not less than 0.2.
[0052] Step S1102: if the time series change rate is less than 0, it is determined that the structure of the target coal seam is in a compaction process.
[0053] Step S1103: If the time series variation is greater than 0, it is determined that the structure of the target coal seam is in a slow destruction stage, wherein a greater time series variation indicates a greater destruction speed.
[0054] From the above description, it can be seen that the embodiment of the present invention provides a method for processing apparent resistivity data based on direct current inversion, which has the following technical effects compared with the prior art:
[0055] 1. The present invention arranges DC electrical measurement lines at the excavation and mining working faces and utilizes normalization processing and gradient calculation methods to avoid the need for full-process drilling or cumbersome on-site testing in traditional methods, thereby significantly reducing the testing workload and construction costs.
[0056] 2. The present invention uses the direct current method to continuously test the entire area to be tested, obtains high-resolution apparent resistivity field and its time series changes, breaks through the limitation of the existing method of "using one point to represent the whole", achieves comprehensive coverage of risks in the area, and provides a more intuitive and accurate risk assessment.
[0057] 3. The present invention combines normalized resistivity value and gradient calculation to accurately identify broken areas and unstable areas during excavation and mining, and dynamically reflects the distribution of dangerous areas during excavation and mining through time series change analysis. Compared with traditional electrical exploration, it has stronger scientificity and accuracy.
[0058] 4. The present invention realizes the identification and early warning functions of dangerous areas during excavation and mining by dynamically monitoring the normalized resistivity and its gradient, thereby greatly improving the safety during excavation and mining.
[0059] Embodiment 2
[0060] Figure 2 This is an overall schematic diagram of a direct current method laying, installation and testing provided according to an embodiment of the present invention. Figure 3 , a direct current method for laying and installing a test section according to an embodiment of the present invention is provided. , a method for processing apparent resistivity data based on direct current method inversion according to an embodiment of the present invention is provided in such a Figure 2 and Figure 3 The application steps in the coal seam shown are as follows:
[0061] S201, based on the characteristics of the coal seam area 2, select the bottom extraction tunnel 4 below the coal seam as the test site.
[0062] S202: A DC electrical method measuring line is laid out at the selected bottom extraction tunnel test location, 48 electrodes 6 are arranged, the distance between adjacent electrodes is 3 m, the drilling depth is 0.2 m, and the electrodes 6 are connected to the DC electrical method instrument 5 through the measuring line 7;
[0063] S203, start the direct current resistivity instrument 5, test the apparent resistivity of the area, and obtain four-dimensional electric field resistivity data during coal seam excavation and mining;
[0064] S204, performing apparent resistivity cloud map inversion on the measured data, exporting the inverted data and setting the required data type;
[0065] S205, normalizing the initial data and calculating the normalized apparent resistivity distribution of each four-dimensional coordinate;
[0066] S206, using the natural neighbor method to determine the normalized λ distribution cloud map;
[0067] S207, calculating a normalized gradient on the normalized cloud image to determine a boundary of geological medium change;
[0068] S208, calculating the time series change rate of the normalized cloud graph to determine the time series change situation;
[0069] S209, calculated rate of change When , it indicates that the structure is in the process of compaction; When , it indicates that the structure is in the slow destruction stage, and the larger the value, the faster the destruction speed.
[0070] In summary, the present invention takes the difference in conductivity of different media as the physical property basis, and preliminarily determines the coal and rock conditions by conducting DC electrical testing in the coal seam bottom extraction tunnel. Then, the initial measured values are normalized and gradient calculated to identify the broken areas and unstable areas during the excavation and mining process, and determine the potential dangerous areas during the excavation and mining process, and evaluate their danger.
[0071] Embodiment 3
[0072] Figure 4 FIG. 1 is a schematic diagram of an apparent resistivity data processing system based on direct current inversion according to an embodiment of the present invention. Figure 3 As shown, the system includes: an inversion module 10 , a normalization module 20 , a first determination module 30 , a calculation module 40 and a second determination module 50 .
[0073] Specifically, the inversion module 10 is used to perform apparent resistivity cloud map inversion on the four-dimensional electric field data during the excavation and mining of the target coal seam to obtain inverted data;
[0074] A normalization module 20 is used to perform normalization processing on the inverted data to obtain the normalized resistivity distribution of each position of the target coal seam;
[0075] A first determination module 30 is used to determine a normalized λ distribution cloud map of the normalized resistivity distribution based on a natural neighbor interpolation method;
[0076] A calculation module 40, used to respectively calculate the normalized gradient and the time series change rate of the normalized λ distribution cloud map;
[0077] The second determination module 50 is used to determine the geological medium change boundary and the time series change of the target coal seam based on the normalized gradient and the time series change amount respectively.
[0078] Specifically, the second determining module 50 is further configured to:
[0079] The position where the normalized gradient exceeds a preset threshold is determined as the geological medium change boundary of the target coal seam;
[0080] If the time series change rate is less than 0, it is determined that the structure of the target coal seam is in the process of compaction;
[0081] If the time series variation is greater than 0, it is determined that the structure of the target coal seam is in the slow destruction stage.
[0082] The present invention further provides an electronic device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the method provided in the embodiment of the present invention when executing the computer program.
[0083] The present invention also provides a computer-readable storage medium, which stores computer instructions. When the computer instructions are executed by a processor, the method provided in the embodiment of the present invention is implemented.
[0084] It will be apparent to those skilled in the art that the invention is not limited to the details of the exemplary embodiments described above and that the invention can be implemented in other specific forms without departing from the spirit or essential features of the invention. Therefore, the embodiments should be considered exemplary and non-limiting in all respects, and the scope of the invention is defined by the appended claims rather than the foregoing description, and it is intended that all variations falling within the meaning and scope of the equivalent elements of the claims be included in the invention. Any reference numeral in a claim should not be considered as limiting the claim to which it relates.
[0085] In addition, it should be understood that although the present specification is described according to implementation modes, not every implementation mode contains only one independent technical solution. This description of the specification is only for the sake of clarity. Those skilled in the art should regard the specification as a whole. The technical solutions in each embodiment may also be appropriately combined to form other implementation modes that can be understood by those skilled in the art.
Claims
1. A method for processing apparent resistivity data based on direct current inversion, characterized in that: include: Perform apparent resistivity cloud map inversion on the four-dimensional electric field data during the excavation and mining of the target coal seam to obtain the inverted data; Normalizing the inverted data to obtain normalized resistivity distribution at each position of the target coal seam; Determining a normalized λ distribution cloud diagram of the normalized resistivity distribution based on a natural neighbor interpolation method; Respectively calculating the normalized gradient and the time series change rate of the normalized λ distribution cloud map; Based on the normalized gradient and the time series variation, the geological medium variation boundary and the time series variation of the target coal seam are determined respectively.
2. The method according to claim 1, characterized in that: The method further comprises: acquiring four-dimensional electric field data of the target coal seam during excavation and mining by means of electrical equipment.
3. The method according to claim 1, characterized in that: The inverted data is normalized to obtain the normalized resistivity distribution of each position of the target coal seam, including: In the formula, ρ i is the apparent resistivity value of the i-th point in the inverted data; ρ min , max are the minimum and maximum apparent resistivity values in the test area, respectively, and λ is the normalized resistivity value.
4. The method according to claim 1, characterized in that: Calculating the normalized gradient of the normalized lambda distribution cloud map includes: λ G (x,y,z)=λ(x,y,z)Degree In the formula, λ G is the normalized gradient value, λ is the normalized resistivity value, and (x, y, z) is the three-dimensional space coordinate.
5. The method according to claim 1, characterized in that: Calculating the time series change rate of the normalized lambda distribution cloud diagram includes: Where (x, y, z, t) is the four-dimensional space-time coordinate, λ is the normalized resistivity value, is the timing change rate.
6. The method according to claim 1, characterized in that: Based on the normalized gradient and the time series variation, respectively, determining the geological medium variation boundary and time series variation of the target coal seam includes: Determine the position where the normalized gradient exceeds a preset threshold as the geological medium change boundary of the target coal seam; If the time series change rate is less than 0, it is determined that the structure of the target coal seam is in a compaction process; If the time series variation is greater than 0, it is determined that the structure of the target coal seam is in a slow destruction stage.
7. An apparent resistivity data processing system based on direct current electrical inversion, characterized in that: include: Inversion module, normalization module, first determination module, calculation module and second determination module; wherein, The inversion module is used to perform apparent resistivity cloud map inversion on the four-dimensional electric field data during the excavation and mining of the target coal seam to obtain inverted data; The normalization module is used to perform normalization processing on the inverted data to obtain the normalized resistivity distribution of each position of the target coal seam; The first determining module is used to determine a normalized λ distribution cloud map of the normalized resistivity distribution based on a natural neighbor interpolation method; The calculation module is used to respectively calculate the normalized gradient and the time series change rate of the normalized λ distribution cloud map; The second determination module is used to determine the geological medium change boundary and the time series change of the target coal seam based on the normalized gradient and the time series change amount respectively.
8. The system according to claim 7, characterized in that: The second determining module is further used for: Determine the position where the normalized gradient exceeds a preset threshold as the geological medium change boundary of the target coal seam; If the time series change rate is less than 0, it is determined that the structure of the target coal seam is in a compaction process; If the time series variation is greater than 0, it is determined that the structure of the target coal seam is in a slow destruction stage.
9. An electronic device, characterized in that: include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the method according to any one of claims 1 to 6 when executing the computer program.
10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and when the computer instructions are executed by a processor, the method according to any one of claims 1 to 6 is implemented.