Joint Inversion Method for Radio Wave Perspective and Borehole Radar Transmission Data in Working Face

Abstract

The present invention relates to the technical field of coal mine exploration, and specifically relates to a method for joint inversion of radio wave penetration and trough wave seismic transmission data in a working face. According to the collected data, the respective energy attenuation coefficients of radio wave penetration and trough wave seismic transmission in the current coal mining working face are calculated to obtain the corresponding electromagnetic wave attenuation coefficient and trough wave attenuation coefficient respectively; these two energy attenuation coefficients are normalized into attribute data corresponding to each discrete grid of the working face; according to the preset anomaly judgment threshold, the attribute data is judged for anomalies, and the attribute data with an energy attenuation coefficient greater than the anomaly judgment threshold is respectively screened out; according to the two types of attribute data, the spatial position relationship of the corresponding geological anomalies is analyzed, and through calculation, the intersection and union corresponding to the distribution area of the corresponding geological anomalies are generated. This solution can improve the accuracy and resolution of the corresponding exploration during geological anomaly exploration.

Description

Technical Field

[0001] The present invention relates to the technical field of coal mine exploration, and particularly relates to a joint inversion method for radio wave penetration and trough wave seismic transmission data of a working face. Background Art

[0002] Coal mining tunneling has great ambiguity and uncertainty, bringing great risks to mine operations. When tunneling without knowing the information to be mined, underground accidents and complex situations (such as water inrush, gas outburst, etc.) occur frequently, and more seriously, it may even lead to casualties of personnel, causing huge economic losses and affecting the efficiency of the mine. Therefore, it is particularly important to detect and locate the geological abnormal positions of the coal mining working face.

[0003] When locating geological anomalies, existing methods include radio wave penetration method and trough wave seismic transmission method. Among them, trough wave seismic exploration has good spatial resolution for geological anomalies such as faults and collapse columns based on seismic wave transmission method, and can accurately reflect the strike information of faults and the structural development interface; radio wave penetration is based on the attenuation characteristics of the electromagnetic wave field, and has high reliability for the structural development area and water-rich geological anomalies, but its vertical resolution is relatively low.

[0004] Of course, after completing the location of geological anomalies in the coal mining working face, it is also very important to monitor the mining equipment during the mining activities of operators and mining equipment, especially when carrying out mining activities near geological anomalies.

[0005] Based on this, a joint inversion method for radio wave penetration and trough wave geological projection method of a coal mining working face is needed, which can improve the accuracy and resolution of the corresponding exploration when exploring geological anomalies. Summary of the Invention

[0006] The present invention aims to provide a joint inversion method for radio wave penetration and trough wave seismic transmission data of a working face, which can improve the accuracy and resolution of the corresponding exploration when exploring geological anomalies.

[0007] To achieve the above object, the present invention adopts the following technical scheme: A joint inversion method for radio wave penetration and trough wave seismic transmission data of a working face, comprising the following steps:

[0008] A data preprocessing step, according to the collected data, calculate the respective energy attenuation coefficients of radio wave penetration and trough wave seismic transmission methods for the current coal mining working face, and respectively obtain the corresponding electromagnetic wave attenuation coefficient and trough wave attenuation coefficient;

[0009] A data normalization step, normalize these two energy attenuation coefficients into attribute data corresponding to each discrete working face grid;

[0010] Anomaly screening step: According to the preset anomaly judgment threshold, perform anomaly judgment on the attribute data, and respectively screen out the attribute data whose energy attenuation coefficient is greater than the anomaly judgment threshold.

[0011] Spatial cross-fusion inversion step: According to two types of attribute data, analyze the spatial position relationship of the corresponding geological anomalies, and through calculation, generate the intersection and union corresponding to the distribution area of the corresponding geological anomalies.

[0012] The principle and advantages of this solution are as follows: In this solution, first, the geological conditions of the current coal mining face are detected by using the radio wave penetration method and the trough wave seismic transmission method. Then, according to these two methods, the respective corresponding energy attenuation coefficients are calculated. Both of these methods independently complete the data collection and calculation, and the data processing process is independent, and the results have a certain similarity in structure. In order to perform comprehensive analysis later, the data will be normalized before data processing to make the standards of these two types of data consistent. Then, by setting a reasonable anomaly judgment threshold, the distribution range of geological anomalies is screened out, the spatial position relationship of these two types of screened geological anomalies is analyzed, and through calculating the intersection and union of the anomaly distributions, the precise positioning of the geological structure is achieved.

[0013] This solution is a joint inversion method developed based on the similarity of the spatial distribution of different physical property parameters. This method does not participate in the process of geophysical inversion, but only constrains the similarity of the geometric structure of the target body during the inversion process, that is, the radio wave penetration method and the trough wave seismic transmission method are similar in structure. The joint inversion can be carried out on the basis of each method fitting its own data, and the inversion results are spatially fused to obtain more accurate detection results.

[0014] At the same time, compared with the processing results of the existing single inversion method or the independently interpreted results, this solution does not need to make any assumptions about various models. The spatial distributions of various physical property parameters can be consistent or inconsistent. The determination of geological anomalies is completed through the set operation of the anomaly areas. The algorithm is simple to implement, the effect is prominent, and the obtained detection results are more accurate and efficient, that is, the accuracy and resolution of the corresponding detection and interpretation are increased.

[0015] Regarding the current exploration of hidden geological structure anomalies in coal mining faces, only two geophysical exploration methods, namely the radio wave penetration method and the trough wave seismic exploration method, are effectively applied underground. Due to the complexity of geological conditions, the limitations of single-method implementation are obvious. To save costs, only single-method detection is often carried out, which can meet the needs of conventional mining. For the higher requirements of current intelligent mining for the control accuracy of geological structures, in order to achieve more beneficial effects, the conventional joint inversion method must obtain the same geological and geophysical model based on the mutual relationship between the physical and geometric parameters of the geological body, and a common physical interface must be ensured. However, for inversions based on different physical properties, their physical bases are different, and complex and cumbersome mathematical conversions of physical parameters are required, and it is difficult to ensure the reliability of the converted data, making it difficult to achieve. The spatial intersection method proposed in this solution is based on the spatial distribution of the physical properties of underground geological bodies, aiming to find the common spatial distribution characteristics among various models, without making any assumptions about various models, only constraining the similarity of the inversion geometric structure, and performing spatial fusion on the basis of each method fitting its own data to obtain more accurate detection results.

[0016] Preferably, as an improvement, it further includes an image display step. According to the intersection and union corresponding to the distribution area of geological anomalies, image reconstruction is performed to generate a corresponding abnormal distribution image of the working face and display it graphically.

[0017] Beneficial effects: The visualization of geological anomaly points is realized through graphical display, making the observation of geological anomaly points more intuitive and convenient.

[0018] Preferably, as an improvement, the normalization formula in the data normalization step is:

[0019]

[0020]

[0021] Where: β 1 is the electromagnetic wave attenuation coefficient, β 2 is the trough wave attenuation coefficient, k 1 i and k2 j are the corresponding attribute data respectively.

[0022] Preferably, as an improvement, it further includes:

[0023] Anomaly point exclusion step. According to the corresponding abnormal distribution image of the working face, anomaly investigation is carried out on each anomaly point in the image, and each anomaly point is classified according to the corresponding anomaly cause. The levels include the first level, the second level, and the third level.

[0024] Beneficial effects: After obtaining each abnormal point under the corresponding coal mine, through the investigation and grading of the abnormal points, the situation of the abnormal points can be well judged, so that early warnings can be given for some relatively dangerous abnormal points.

[0025] Preferably, as an improvement, it further includes:

[0026] An early warning step, according to the abnormal distribution image of the corresponding working face, judge whether the position information of the current mining equipment is close to the abnormal points in the abnormal distribution image of the working face and the levels of the corresponding abnormal points, and give corresponding alarms according to the corresponding distances and levels.

[0027] Beneficial effects: In this solution, when the operator enters the coal mine, the position information of the mining equipment will be located in real time to ensure that the operator does not approach the abnormal points, thus avoiding accidents.

[0028] Preferably, as an improvement, the early warning step specifically includes:

[0029] A position positioning step, detect the position information of the current mining equipment in real time and generate the corresponding position information;

[0030] A distance judgment step, using the abnormal distribution image of the working face, judge whether the distance from the position information of the mining equipment to each abnormal point is less than a preset first threshold; if so, proceed to the next step;

[0031] A level identification step, identify the abnormal points whose distances from the position information to each abnormal point are less than the preset first threshold, and identify the levels of the corresponding abnormal points;

[0032] An early warning step, if the level of the corresponding abnormal point is the first level, immediately send an alarm to the operator and send a corresponding alarm signal to the remote server; if the level of the corresponding abnormal point does not exist at the first level but exists at the second level, send an alarm signal to the operator; if the level of the corresponding abnormal point is only the third level, send a warning signal to the operator.

[0033] Beneficial effects: In this solution, the basis for judging whether the position where the operator is located is safe is not only the distance from the corresponding abnormal point, but also the level of the corresponding abnormal point. The higher the level, the greater the possibility of danger. When approaching an abnormal point with a high level, not only an alarm is given to the operator, but also an alarm is given to the staff at the remote end. This can timely dissuade the operator and enable the staff at the remote end to timely understand the position of the operator, greatly improving the property safety of the people.

[0034] Preferably, as an improvement, it further includes: a mining plan formulation step, which identifies the corresponding signal type according to the received signal and formulates a corresponding mining plan.

[0035] Beneficial effects: After the operator receives the signal, different signals mean different situations faced during mining. By identifying the signal type, the corresponding risk level can be accurately known. In this way, when the operator uses the mining equipment for mining, a corresponding mining plan can be formulated in a timely manner, which greatly improves the mining efficiency and safety. Brief Description of the Drawings

[0036] Figure 1 It is a flowchart of the joint inversion method of radio wave perspective and trough wave seismic transmission data in the first embodiment of the present invention for the working face. Detailed Embodiments

[0037] The following is further detailed through specific embodiments:

[0038] The embodiment is basically as shown in the attached Figure 1 The joint inversion method of radio wave perspective and trough wave seismic transmission data for the working face includes the following steps:

[0039] Data preprocessing step: Calculate the respective energy attenuation coefficients of radio wave perspective and trough wave seismic transmission methods for the current coal mining working face according to the collected data, and obtain the corresponding electromagnetic wave attenuation coefficient and trough wave attenuation coefficient respectively;

[0040] Data normalization step: Normalize these two energy attenuation coefficients into attribute data corresponding to each discrete working face grid; In this embodiment, the electromagnetic wave attenuation coefficient ( β 1) and the trough wave attenuation coefficient ( β 2), where β 1 and β 2 have the same number of rows and are the same as the number of discrete grids.

[0041] The normalization formula in the data normalization step is:

[0042]

[0043]

[0044] Where: β 1 is the electromagnetic wave attenuation coefficient, β 2 is the trough wave attenuation coefficient, k 1 i and k2 j are the corresponding attribute data respectively.

[0045] Anomaly screening step: According to the preset anomaly judgment threshold, perform anomaly judgment on the attribute data, and respectively screen out the attribute data with the energy attenuation coefficient greater than the anomaly judgment threshold; in this embodiment, for the setting of the anomaly judgment threshold, it is mainly determined according to the known or presumed geological structure.

[0046] Spatial cross-fusion inversion step: According to two types of attribute data, analyze the spatial position relationship of the corresponding geological anomalies, and through calculation, generate the intersection and union corresponding to the distribution area of the corresponding geological anomalies.

[0047] Specifically,

[0048]

[0049] In the formula, k 1 th and k 2 th are respectively the anomaly judgment thresholds, that is, the area with a higher attenuation coefficient is identified as the anomaly area; Ω1 and Ω2 are the anomaly distribution ranges of two different fusion methods of exploration results.

[0050] Image display step: According to the intersection and union corresponding to the distribution area of the geological anomalies, perform image reconstruction, generate the corresponding working face anomaly distribution image and display it graphically.

[0051] Anomaly point exclusion step: According to the corresponding working face anomaly distribution image, conduct anomaly investigation on each anomaly point in the image, and classify each anomaly point according to the corresponding anomaly cause.

[0052] Warning step: According to the corresponding working face anomaly distribution image, judge whether the position information of the current mining equipment is close to the anomaly points in the working face anomaly distribution image and the levels of the corresponding anomaly points, and give corresponding alarms according to the corresponding distances and levels

[0053] The warning step specifically includes:

[0054] Position positioning step: Real-time detect the position information of the current mining equipment and generate the corresponding position information;

[0055] Distance judgment step: Use the working face anomaly distribution image to judge whether the distance from the position information of the mining equipment to each anomaly point is less than the preset first threshold; if so, proceed to the next step;

[0056] Level identification step: Identify the anomaly points whose distance from the position information to each anomaly point is less than the preset first threshold, and identify the levels of the corresponding anomaly points;

[0057] Warning step: If the level of the corresponding abnormal point is at the first level, an alarm will be immediately sent to the operator, and a corresponding alarm signal will be sent to the remote server; if the level of the corresponding abnormal point does not exist at the first level but exists at the second level, an alarm signal will be sent to the operator; if the level of the corresponding abnormal point is only at the third level, a warning signal will be sent to the operator.

[0058] In this embodiment, after determining the geological abnormal points, when the corresponding operator uses the mining equipment to carry out mining activities on this coal mining face, in order to ensure the safety of the operator and the mining equipment, during the mining of the corresponding mining equipment, the position information of the mining equipment will be monitored in real time, that is, the position of the operator will be located in real time. By comparing the position information of the mining equipment with the positions of each abnormal point in the working abnormal distribution image, the corresponding abnormal points with a distance less than the preset first threshold are matched, and the level of the abnormal point is identified. The critical degree of the position where the operator is located is judged by using the dual judgment of distance and level, and then different warnings are given according to different critical degrees. For example, when Party A is operating the mining equipment of the coal mine and the mining equipment is continuously mining into the coal seam, the corresponding mining equipment will get closer and closer to each abnormal point. When the mining equipment approaches the corresponding abnormal points A, B, and C, if it is detected that the distances between the mining equipment and these three abnormal points are all less than the preset first threshold, then when identifying the levels of these three abnormal points, the levels of all three abnormal points will be identified; if only the distance of one abnormal point is less than the preset first threshold, only the level of one point will be identified. When it is identified that the level of the corresponding abnormal point exists at the first level, this indicates that the position where the operator is located is extremely dangerous. Therefore, an alarm will be sent to the operator immediately, and an alarm signal will also be sent to the remote staff synchronously to notify the remote staff, so as to avoid the problem that the remote staff do not know the position of the operator after an accident occurs, resulting in untimely and inaccurate rescue.

[0059] It further includes: a mining plan formulation step, which identifies the corresponding signal type according to the received signal and formulates a corresponding mining plan. In this embodiment, when the operator receives the corresponding signal, the signal type of the signal will be identified, so that the operator can formulate the corresponding mining plan in the first time, so that corresponding countermeasures can be taken in time when the mining equipment mines to the position of the corresponding abnormal point, avoiding accidents and greatly improving the safety and reliability of the mining work.

[0060] The above are only embodiments of the present invention, and common knowledge such as specific technical solutions and / or characteristics in the solutions are not described in detail herein. It should be noted that for those skilled in the art, without departing from the technical solution of the present invention, several deformations and improvements can be made, which should also be regarded as the protection scope of the present invention, and these will not affect the implementation effect of the present invention and the practicability of the patent. The protection scope required by this application shall be subject to the content of its claims, and the specific implementation manners and the like recorded in the specification can be used to interpret the content of the claims.

Claims

1. Joint inversion method for radio wave perspective and trough wave seismic transmission data of working face, characterized in that: It includes the following steps: Data preprocessing step: According to the collected data, calculate the respective energy attenuation coefficients of radio wave perspective and trough wave seismic transmission method for the current coal mining working face, and obtain the corresponding electromagnetic wave attenuation coefficient and trough wave attenuation coefficient respectively; Data normalization step: Normalize these two energy attenuation coefficients into attribute data corresponding to the discrete working face grid one by one; Abnormality screening step: According to the preset abnormality judgment threshold, judge the attribute data for abnormalities, and respectively screen out the attribute data whose energy attenuation coefficient is greater than the abnormality judgment threshold; Spatial cross-fusion inversion step: According to the two attribute data, analyze the spatial position relationship of the corresponding geological abnormalities, and generate the intersection and union corresponding to the distribution area of the corresponding geological abnormalities through calculation; It also includes an image display step: According to the intersection and union corresponding to the distribution area of the geological abnormalities, perform image reconstruction, generate the corresponding working face abnormality distribution image and display it graphically; Abnormality point investigation step: According to the corresponding working face abnormality distribution image, investigate each abnormality point in the image, and classify each abnormality point according to the on-site geological information and the abnormality value range. The levels include the first level, the second level and the third level; Early warning step: According to the corresponding working face abnormality distribution image, judge whether the position information of the current mining equipment is close to the abnormality points in the working face abnormality distribution image and the levels of the corresponding abnormality points, and give corresponding alarms according to the corresponding distance and level.

2. The joint inversion method for radio wave perspective and trough wave seismic transmission data of working face according to claim 1, characterized in that: The normalization formula in the data normalization step is: Wherein: β 1 is the electromagnetic wave attenuation coefficient, β 2 is the trough wave attenuation coefficient, k 1 i and k2 j are the corresponding attribute data respectively.

3. The joint inversion method for radio wave perspective and trough wave seismic transmission data of working face according to claim 1, characterized in that: The early warning step specifically includes: Position positioning step: Real-time detect the position information of the current mining equipment and generate the corresponding position information; Distance judgment step: Use the working face abnormality distribution image to judge whether the distance from the mining equipment position information to each abnormality point is less than the preset first threshold; if so, proceed to the next step; Level identification step: Identify the abnormality points whose distance from the position information to each abnormality point is less than the preset first threshold, and identify the levels of the corresponding abnormality points; Early warning step: If the level of the corresponding abnormality point is the first level, immediately send an alarm to the operator and send a corresponding alarm signal to the remote server; if the level of the corresponding abnormality point does not exist in the first level but exists in the second level, send an alarm signal to the operator; if the level of the corresponding abnormality point is only the third level, send a warning signal to the operator.

4. The joint inversion method for radio wave perspective and trough wave seismic transmission data of working face according to claim 3, characterized in that: It also includes: Mining plan formulation step: According to the received signal, identify the corresponding signal type and formulate the corresponding mining plan.

Images