Abnormal body positioning method and system combining borehole transient electromagnetism and geological radar

Abstract

The invention discloses an anomalous body positioning method and system combining borehole transient electromagnetic and geological radar, and relates to the technical field of mine safety detection, and the method comprises the steps: synchronously collecting the borehole transient electromagnetic data and geological radar data of an ore body to be analyzed, and generating a borehole transient electromagnetic imaging graph and a borehole geological radar imaging graph; two-dimensional coordinates of a target abnormal point in the borehole geological radar imaging graph are extracted, and the corresponding target apparent resistivity is determined; drawing an apparent resistivity contour map based on the target apparent resistivity in a target radial sectioning map of the borehole transient electromagnetic imaging map, and determining a plurality of candidate three-dimensional coordinates of a target abnormal point; at least one target feature point is selected from the borehole geological radar imaging graph, and the feature apparent resistivity corresponding to the target feature point is determined; and drawing a feature circle based on the target feature points, and screening the plurality of candidate three-dimensional coordinates to obtain coordinates of the geological anomalous body. According to the invention, the technical problem of insufficient positioning precision and reliability in the prior art is solved.

Description

Technical Field

[0001] This invention relates to the field of mine safety detection technology, and in particular to a method and system for locating anomalies by combining borehole transient electromagnetic and ground-penetrating radar. Background Technology

[0002] The precise location of geological anomalies is a core technological bottleneck in mineral resource exploration, geological disaster prevention and control, and major underground engineering construction, directly impacting project safety and resource development efficiency. From an industry development perspective, as easily exploitable shallow resources become increasingly depleted, mineral exploration has extended to depths exceeding 1000 meters. Deep geological conditions are more complex, and anomalies are more concealed, making it difficult for traditional technologies to penetrate complex strata for accurate identification. Furthermore, frequent extreme weather events and intensified human engineering activities result in geological disasters such as landslides, mine flooding, and karst cave collapses causing over 10 billion yuan in direct economic losses annually. Therefore, the early and accurate identification of geological anomalies has become crucial for disaster prevention and control.

[0003] Currently, both ground-penetrating radar (GPR) and borehole transient electromagnetic detection (BEM) technologies have inherent limitations when used individually. GPR achieves high resolution and rapid imaging using high-frequency electromagnetic waves, but these waves attenuate drastically in water-rich strata and highly conductive rock layers, limiting its detection depth to typically no more than 30 meters, and only 10-15 meters in complex strata. Furthermore, it can only provide two-dimensional profile information and cannot achieve three-dimensional calibration. For example, in tunnel advance detection, side caves are easily misidentified as obstacles directly ahead. BEM, on the other hand, can reach depths of hundreds of meters and is sensitive to low-resistivity anomalies, but its lateral resolution is only 5-10 meters, making it difficult to identify small anomalies. Moreover, due to the heterogeneity of the underground medium, multiple suspected anomalies are easily generated, making it impossible to directly determine the location of a unique anomaly. This often results in borehole failures during mineral exploration.

[0004] Existing joint detection technologies have failed to overcome core bottlenecks: First, data fusion is low, with asynchronous acquisition of the two types of data, inconsistent coordinate systems, and a lack of standardized coordinate transformation; second, they rely on experience-based judgment, simply overlaying anomalous areas without establishing quantitative models, resulting in strong subjectivity; third, the screening mechanism is simplistic, easily leading to multiple suspected anomaly coordinates, making it impossible to directly determine the location of a unique geological anomaly and eliminate the ambiguity of anomaly points. These problems result in insufficient positioning accuracy and reliability, failing to meet the precise requirements of deep exploration and major engineering projects. Summary of the Invention

[0005] To address the aforementioned technical problems in existing technologies, embodiments of the present invention provide a method and system for anomaly localization combining borehole transient electromagnetic and ground-penetrating radar. The technical solution is as follows: On one hand, a method for anomaly localization combining borehole transient electromagnetic (BEM) and ground-penetrating radar (GPR) is provided. The method includes: simultaneously acquiring BEM and GPR data of the ore body to be analyzed and preprocessing them to generate a BEM image and a GPR image, respectively; extracting two-dimensional coordinates of target anomalies from the GPR image; the two-dimensional coordinates include the drill rod advance direction coordinates and the detection depth coordinates; determining the target apparent resistivity at the target anomaly based on the BEM distance; the BEM distance is the distance between the BEM and the measured point; drawing an apparent resistivity contour map based on the target apparent resistivity in the target radial section of the BEM image, and then using the apparent resistivity contour map and the depth of the target anomaly. The depth coordinates are used to determine multiple candidate three-dimensional coordinates of the target anomaly point; the target radial section is a radial section of the borehole transient electromagnetic imaging image based on the drill rod advance direction coordinates of the target anomaly point; at least one target feature point is selected in the borehole ground-penetrating radar imaging image, and the characteristic apparent resistivity corresponding to the target feature point is determined based on the transient electromagnetic distance of the target feature point; the target feature point and the target anomaly point have the same coordinates in the drill rod advance direction; a feature circle is drawn in the target radial section image based on the coordinates of the target feature point in the detection depth direction, and the multiple candidate three-dimensional coordinates are filtered based on the characteristic apparent resistivity and the feature circle to obtain the geological anomaly body coordinates of the ore body to be analyzed; wherein, the number of filtering times is consistent with the number of target feature points selected.

[0006] Optionally, the transient electromagnetic data and ground-penetrating radar data of the ore body to be analyzed are collected simultaneously and preprocessed, including: setting transient electromagnetic and ground-penetrating radar at preset intervals in the tunnel borehole of the ore body to be analyzed; collecting the transient electromagnetic data and ground-penetrating radar data of the ore body to be analyzed based on the transient electromagnetic and ground-penetrating radar respectively, and preprocessing them.

[0007] Optionally, determining the target apparent resistivity at the target anomaly point based on the transient electromagnetic distance of the target anomaly point includes: determining the target apparent resistivity at the target anomaly point based on the transient electromagnetic distance of the target anomaly point and the correspondence between the borehole transient electromagnetic distance and apparent resistivity; wherein, the correspondence between the borehole transient electromagnetic distance and apparent resistivity includes:

[0008] In the formula, s is the transient electromagnetic distance during drilling. K For compensation coefficient, t This is the delay time after the transmitting current is turned off. ρ Resistivity μ 0 represents the permeability of free space.

[0009] Optionally, determining multiple candidate three-dimensional coordinates of the target anomaly point based on the apparent resistivity contour map and the detection depth coordinates of the target anomaly point includes: drawing a target circle based on the detection depth coordinates of the target anomaly point in the target radial section map; determining the three-dimensional coordinates of the intersection of the apparent resistivity contour map and the target circle in the borehole transient electromagnetic imaging map as multiple candidate three-dimensional coordinates; wherein, the equation of the target circle includes The multiple candidate three-dimensional coordinates are satisfy , Let be the two-dimensional coordinates of the target anomaly point, Y be the coordinates of the drill pipe advancing direction, and r0 be the coordinates of the detection depth.

[0010] Optionally, selecting at least one target feature point in the borehole ground-penetrating radar image and determining the characteristic apparent resistivity at the target feature point based on the transient electromagnetic distance of the target feature point includes: selecting at least one target feature point in the borehole ground-penetrating radar image; the two-dimensional coordinates of the target feature point... r p To select parameters, r max A preset upper limit threshold is set; based on the two-dimensional coordinates of the target feature point, the transient electromagnetic distance of the target feature point is calculated; the formula for calculating the transient electromagnetic distance of the target feature point includes: , l The spacing between transient electromagnetic and ground-penetrating radars installed in the borehole of the ore body to be analyzed is determined; based on the transient electromagnetic distance of the target feature point and the correspondence between the transient electromagnetic distance of the borehole and the apparent resistivity, the characteristic apparent resistivity at the target feature point is determined.

[0011] Optionally, the geological anomaly coordinates of the ore body to be analyzed are obtained by filtering the plurality of candidate three-dimensional coordinates based on the characteristic apparent resistivity and the characteristic circle, including: starting from the origin of the coordinates of the target radial section and connecting it to the plurality of candidate three-dimensional coordinates respectively and extending it until it intersects the characteristic circle at a plurality of corresponding characteristic intersection points; calculating the apparent resistivity corresponding to each characteristic intersection point based on the coordinates of the characteristic intersection points; calculating the error between the apparent resistivity corresponding to each characteristic intersection point and the characteristic apparent resistivity, and removing the candidate three-dimensional coordinates with errors greater than a preset error value to obtain the filtered three-dimensional coordinates; if the filtered candidate three-dimensional coordinates are unique, then the filtered three-dimensional coordinates are determined as the geological anomaly coordinates of the ore body to be analyzed.

[0012] Optionally, if the candidate three-dimensional coordinates after screening are not unique, other feature points are selected and screened repeatedly until the candidate three-dimensional coordinates after screening are unique, and the unique candidate three-dimensional coordinates are determined as the geological anomaly coordinates of the ore body to be analyzed.

[0013] On the other hand, a combined borehole transient electromagnetic and ground-penetrating radar (GPR) anomaly localization system is also provided to implement the combined borehole transient electromagnetic and GPR anomaly localization method provided in the embodiments of the present invention. The system includes: an acquisition module, an extraction module, a determination module, a candidate module, a selection module, and a filtering module. The acquisition module is used to simultaneously acquire borehole transient electromagnetic data and GPR data of the ore body to be analyzed and preprocess them to generate a borehole transient electromagnetic imaging map and a borehole GPR imaging map, respectively. The extraction module is used to extract the two-dimensional coordinates of the target anomaly points in the borehole GPR imaging map; the two-dimensional coordinates include the drill rod advance direction coordinates and the detection depth coordinates. The determination module is used to determine the target apparent resistivity corresponding to the target anomaly point based on the transient electromagnetic distance of the target anomaly point; the transient electromagnetic distance is the distance between the transient electromagnetic field and the point to be measured. The candidate module is used to select the target anomaly in the radial section of the borehole transient electromagnetic imaging map. A resistivity contour map based on the target resistivity is drawn, and multiple candidate three-dimensional coordinates of the target anomaly are determined based on the resistivity contour map and the detection depth coordinates of the target anomaly. The target radial section is a radial section of the borehole transient electromagnetic imaging map based on the drill rod advance direction coordinates of the target anomaly. The selection module is used to select at least one target feature point in the borehole ground-penetrating radar imaging map, and determine the characteristic resistivity corresponding to the target feature point based on the transient electromagnetic distance of the target feature point. The target feature point and the target anomaly have the same coordinates in the drill rod advance direction. The filtering module is used to draw a feature circle based on the target feature point's coordinates in the detection depth direction in the target radial section map, and filter the multiple candidate three-dimensional coordinates based on the characteristic resistivity and the feature circle to obtain the geological anomaly coordinates of the ore body to be analyzed. The number of filtering times is consistent with the number of target feature points selected.

[0014] On the other hand, an electronic device is also provided, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the method provided in the embodiments of the present invention.

[0015] On the other hand, a computer-readable storage medium is also provided, wherein program code is stored in the computer-readable storage medium, and the program code can be called by a processor to execute the method provided in the embodiments of the present invention.

[0016] This invention provides a method and system for locating geological anomalies by combining borehole transient electromagnetic and ground-penetrating radar. The method uses a combined detection mode of borehole transient electromagnetic and ground-penetrating radar to make up for the deficiencies of single detection technologies. It also eliminates positioning errors through a multi-round coordinate screening mechanism, which significantly improves the accuracy and reliability of locating geological anomalies. This alleviates the technical problems of insufficient positioning accuracy and reliability in existing technologies, which cannot meet the precision requirements of deep exploration and major engineering projects. Attached Figure Description

[0017] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0018] Figure 1 This is a flowchart of an anomaly localization method combining borehole transient electromagnetic and ground-penetrating radar, provided by an embodiment of the present invention; Figure 2 This is a schematic diagram of a cooperative detection system provided in an embodiment of the present invention; Figure 3 This is a borehole ground-penetrating radar imaging image provided in an embodiment of the present invention; Figure 4 This is a schematic diagram illustrating the extraction of candidate three-dimensional coordinates of target anomaly points according to an embodiment of the present invention; Figure 5 This is a schematic diagram of feature point selection in a borehole ground-penetrating radar imaging image provided by an embodiment of the present invention; Figure 6 This is a schematic diagram illustrating the screening of candidate three-dimensional coordinates provided in an embodiment of the present invention; Figure 7 This is a schematic diagram of an anomaly localization system combining borehole transient electromagnetic and ground-penetrating radar, provided in an embodiment of the present invention.

[0019] In the diagram: 1. Drilling rig, 2. Controller, 3. Tunnel, 4. Borehole, 5. Drill rod, 6. Transient electromagnetic, 7. Ground penetrating radar. Detailed Implementation

[0020] The technical solution of the present invention will now be described with reference to the accompanying drawings.

[0021] In embodiments of the present invention, words such as "exemplarily," "for example," etc., are used to indicate that something is an example, illustration, or description. Any embodiment or design described as "exemplary" in the present invention should not be construed as being more preferred or advantageous than other embodiments or designs. Specifically, the use of the word "exemplary" is intended to present the concept in a concrete manner. Furthermore, in embodiments of the present invention, the meaning expressed by "and / or" can be both, or either one.

[0022] To make the technical problems, technical solutions and advantages of the present invention clearer, a detailed description will be given below in conjunction with the accompanying drawings and specific embodiments.

[0023] Figure 1 This is a flowchart of an anomaly localization method combining borehole transient electromagnetic and ground-penetrating radar, provided by an embodiment of the present invention. Figure 1 As shown, the method specifically includes the following steps: Step S102: Simultaneously collect borehole transient electromagnetic data and ground-penetrating radar data of the ore body to be analyzed and preprocess them to generate borehole transient electromagnetic imaging map and borehole ground-penetrating radar imaging map, respectively.

[0024] Step S104: Extract the two-dimensional coordinates of the target anomaly points in the borehole ground-penetrating radar image; the two-dimensional coordinates include the coordinates of the drill rod advance direction and the detection depth.

[0025] Step S106: Based on the transient electromagnetic distance of the target anomaly point, determine the target apparent resistivity corresponding to the target anomaly point; the transient electromagnetic distance is the distance between the transient electromagnetic field and the point to be measured.

[0026] Step S108: Draw an apparent resistivity contour map based on the apparent resistivity of the target in the target radial section map of the borehole transient electromagnetic imaging image, and determine multiple candidate three-dimensional coordinates of the target anomaly point based on the apparent resistivity contour map and the detection depth coordinates of the target anomaly point; the target radial section map is a radial section map of the drill rod advance direction coordinates of the target anomaly point in the borehole transient electromagnetic imaging image.

[0027] Step S110: Select at least one target feature point in the borehole ground-penetrating radar image, and determine the characteristic apparent resistivity at the target feature point based on the transient electromagnetic distance of the target feature point; the target feature point and the target anomaly point have the same coordinates in the drill rod advance direction.

[0028] Step S112: Draw a feature circle based on the coordinates of the target feature points in the detection depth direction in the radial section map of the target, and filter multiple candidate three-dimensional coordinates based on the feature resistivity and the feature circle to obtain the coordinates of the geological anomaly of the ore body to be analyzed; wherein, the number of filtering times is consistent with the number of target feature points selected.

[0029] Specifically, step S102 further includes the following steps: Step S1021: In the tunnel borehole of the ore body to be analyzed, transient electromagnetic and ground-penetrating radar are installed at preset intervals.

[0030] Step S1022: Based on transient electromagnetic data and ground-penetrating radar (GPR) data, borehole transient electromagnetic data and GPR data of the ore body to be analyzed are acquired and preprocessed. Finally, borehole transient electromagnetic imaging map and borehole GPR imaging map are generated based on the borehole transient electromagnetic data and GPR data, respectively.

[0031] Figure 2 This is a schematic diagram of a cooperative detection system provided according to an embodiment of the present invention. Figure 2 As shown, the collaborative detection system used in this embodiment of the invention includes a drilling rig 1, a controller 2, a drill rod 5, a transient electromagnetic system 6, and a ground-penetrating radar 7. The drilling rig 1 drills boreholes 4 within the tunnel 3, and the drill rods 5 within the boreholes 4 are spaced at preset intervals. l Transient electromagnetic radar 6 and ground-penetrating radar 7 are respectively installed.

[0032] Preferably, the transient electromagnetic 6 is located behind the ground-penetrating radar 7.

[0033] Preferably, in order to achieve optimal detection results while preventing mutual interference between the transient electromagnetic radar 6 and the ground-penetrating radar 7, a preset distance is set. l =3m.

[0034] In an optional embodiment of this invention, a three-dimensional coordinate system for borehole detection is established, defining the drill rod advance direction as the Y-axis, the detection depth direction as the Z-axis, and the direction perpendicular to the YOZ plane as the X-axis. In the aforementioned three-dimensional coordinate system for borehole detection, the coordinates of anomalies in the borehole transient electromagnetic imaging image are defined as (x, y, z), and the coordinates of anomalies in the borehole ground-penetrating radar imaging image are defined as (y, r), where... .

[0035] In step S104, for the geological anomaly to be located, the target anomaly point a is identified in the two-dimensional image of the borehole ground-penetrating radar image, and the two-dimensional coordinates of the target anomaly point in the borehole ground-penetrating radar image are recorded. Where Y is the coordinate of the drill pipe advancing direction, and r0 is the coordinate of the detection depth.

[0036] Specifically, step S106 further includes: determining the target apparent resistivity at the target anomaly point based on the correspondence between the transient electromagnetic distance of the target anomaly point and the transient electromagnetic distance of the borehole and the apparent resistivity.

[0037] Specifically, based on the detection depth coordinate r0 of target anomaly point a, the distance between the transient electromagnetic field and the target anomaly point is obtained, i.e., the transient electromagnetic distance of the target anomaly point: .

[0038] Specifically, the correspondence between borehole transient electromagnetic distance and apparent resistivity includes:

[0039] In the formula, s is the transient electromagnetic distance during drilling. K For compensation coefficient, t This is the delay time after the transmitting current is turned off. ρ Resistivity μ 0 represents the permeability of free space.

[0040] Finally, based on the above correspondence, the target apparent resistivity was determined to be... .

[0041] Specifically, step S108 further includes the following steps: Step S1081: In the target radial section diagram, draw a target circle based on the detection depth coordinates of the target anomaly points. The equation of the target circle includes... .

[0042] Step S1082: Determine multiple candidate three-dimensional coordinates for the three-dimensional intersection of the apparent resistivity contour map and the target circle in the borehole transient electromagnetic imaging map.

[0043] Specifically, in the three-dimensional image of the borehole transient electromagnetic imaging map y = Y At this location, a two-dimensional apparent resistivity contour map is obtained by radial sectioning, and the apparent resistivity is then selected. ρ = ρ a contour map and target circle The set of intersection points is used as multiple candidate 3D coordinates, denoted as point. ,satisfy .

[0044] Specifically, step S110 further includes the following steps: Step S1101: Select at least one target feature point in the borehole ground-penetrating radar image.

[0045] Step S1102: Based on the two-dimensional coordinates of the target feature point, calculate the transient electromagnetic distance of the target feature point; the formula for calculating the transient electromagnetic distance of the target feature point includes: , l The spacing between transient electromagnetic and ground-penetrating radars installed in the boreholes of the ore body to be analyzed.

[0046] Step S1103: Based on the correspondence between the transient electromagnetic distance of the target feature point and the transient electromagnetic distance of the borehole and the apparent resistivity, determine the characteristic apparent resistivity at the target feature point.

[0047] In this embodiment of the invention, feature points are introduced to determine and filter multiple candidate three-dimensional coordinates. Specifically, a target feature point is defined as any point in the geological imaging map, wherein the target feature point and the target anomaly point have the same Y-axis coordinate, and the coordinate value r is any coordinate value not exceeding a preset upper limit threshold, i.e., the two-dimensional coordinate of the target feature point. r p To select parameters, r max This is a preset upper limit threshold.

[0048] Then calculate the distance between the target feature point and the transient electromagnetic field: Based on the correspondence between borehole transient electromagnetic distance and apparent resistivity, the characteristic apparent resistivity of the target feature point p in the borehole transient electromagnetic field is determined. .

[0049] Specifically, step S112 further includes the following steps: Step S1121: Starting from the origin of the target radial section, connect and extend the coordinates of multiple candidate three-dimensional coordinates until they intersect with the feature circle at multiple corresponding feature intersection points.

[0050] Step S1122: Calculate the apparent resistivity at each feature intersection point based on the coordinates of the feature intersection points.

[0051] Step S1123: Calculate the error between the apparent resistivity and the feature apparent resistivity at each feature intersection point, and remove the candidate three-dimensional coordinates with errors greater than the preset error value to obtain the filtered three-dimensional coordinates.

[0052] Step S1124: If the candidate three-dimensional coordinates after screening are unique, then the three-dimensional coordinates after screening are determined as the geological anomaly coordinates of the ore body to be analyzed.

[0053] Step S1125: If the candidate three-dimensional coordinates after screening are not unique, then the other selected feature points are screened repeatedly until the candidate three-dimensional coordinates after screening are unique, and the unique candidate three-dimensional coordinates are determined as the geological anomaly coordinates of the ore body to be analyzed.

[0054] Specifically, feature circles are drawn in the radial section diagram of the target. Starting from the origin and the i-th candidate 3D coordinates a i Connect and extend a ray until it intersects the feature circle. Intersect at characteristic intersection point The apparent resistivity at the intersection of features is recorded. ; Calculate the error between the apparent resistivity and the characteristic apparent resistivity at each characteristic intersection point. When the error Retain the corresponding candidate 3D coordinates, otherwise discard them.

[0055] If, after filtering the candidate 3D coordinates, a unique outlier coordinate cannot be obtained, the feature point p is selected again to perform a second and third filtering on the outlier after the first filtering, until the final unique and accurate coordinates of the outlier are obtained.

[0056] When i=k, a unique and accurate coordinate of the outlier is obtained, and the point is output. These are the precise coordinates of the geological anomaly.

[0057] The method provided by the embodiments of the present invention will be further described below with reference to specific examples.

[0058] Figure 3 This is a borehole ground-penetrating radar imaging image provided according to an embodiment of the present invention. For example... Figure 3 As shown, the coordinates of target anomaly point a are identified in the borehole ground-penetrating radar image, where the drill rod advance direction Y = 33.2m and r0 = 16.75m. Its two-dimensional coordinates in the borehole ground-penetrating radar image are denoted as follows: The distance between the transient electromagnetic field and the target anomaly point a was calculated. Then, the apparent resistivity at target anomaly point a is obtained by using the correspondence between the transient electromagnetic distance of the borehole and the apparent resistivity. .

[0059] Figure 4 This is a schematic diagram illustrating the extraction of candidate three-dimensional coordinates of target anomaly points according to an embodiment of the present invention. Figure 4 As shown in Figure (a), the radial section of the target is drawn. Transient electromagnetic contour maps of boreholes at that time; such as Figure 4 As shown in Figure (b), draw the target circle. This leads to two intersection points: and , as two candidate three-dimensional coordinates.

[0060] Figure 5 This is a schematic diagram illustrating feature point selection in a borehole ground-penetrating radar imaging image according to an embodiment of the present invention. Figure 5 As shown, select target feature points. The distance between the borehole transient electromagnetic field and the target feature point p is calculated. Based on the correspondence between the transient electromagnetic distance of the borehole and the apparent resistivity, the apparent resistivity corresponding to the target feature point p is obtained as follows: .

[0061] Figure 6This is a schematic diagram illustrating the filtering of candidate three-dimensional coordinates according to an embodiment of the present invention. Figure 6 As shown, draw the feature circle. Candidate 3D coordinates of the origin and the target anomaly point a i Connect and extend a ray until it intersects the feature circle. Intersect at point and points Record the intersection points respectively. Apparent resistivity Intersection Apparent resistivity .

[0062] Then calculate the feature points respectively. and The corresponding apparent resistivity error is: , .

[0063] After feature point screening, geological anomalies a The precise and unique coordinates are a The coordinates corresponding to 1, i.e. .

[0064] As described above, the embodiments of the present invention provide an anomaly localization method combining borehole transient electromagnetic and ground-penetrating radar, overcoming the shortcomings of the prior art and solving problems such as low fusion of borehole transient electromagnetic and ground-penetrating radar data, simple superposition of anomaly areas obtained from the two types of detection without establishing a quantitative model, and the inability to directly determine the location of a unique anomaly by easily obtaining the coordinates of multiple anomaly points. By jointly detecting borehole transient electromagnetic and ground-penetrating radar, the method achieves unique and accurate three-dimensional positioning of geological anomalies. This method has the advantages of complementary technologies, unique anomalies, and high accuracy.

[0065] Figure 7 This is a schematic diagram of an anomaly localization system combining borehole transient electromagnetic and ground-penetrating radar, provided according to an embodiment of the present invention. Figure 7 As shown, the system includes: acquisition module 10, extraction module 20, determination module 30, candidate module 40, selection module 50 and filtering module 60.

[0066] Specifically, the acquisition module 10 is used to simultaneously acquire borehole transient electromagnetic data and ground-penetrating radar data of the ore body to be analyzed and preprocess them to generate borehole transient electromagnetic imaging map and borehole ground-penetrating radar imaging map, respectively. Extraction module 20 is used to extract the two-dimensional coordinates of target anomaly points in the borehole ground-penetrating radar image; the two-dimensional coordinates include the coordinates of the drill rod advance direction and the detection depth. The determination module 30 is used to determine the target apparent resistivity at the target anomaly point based on the transient electromagnetic distance; the transient electromagnetic distance is the distance between the transient electromagnetic field and the point to be measured. Candidate module 40 is used to draw an apparent resistivity contour map based on the apparent resistivity of the target in the target radial section map of the borehole transient electromagnetic imaging image, and to determine multiple candidate three-dimensional coordinates of the target anomaly based on the apparent resistivity contour map and the detection depth coordinates of the target anomaly; the target radial section map is a radial section map based on the drill rod advance direction coordinates of the target anomaly in the borehole transient electromagnetic imaging image. Module 50 is selected to select at least one target feature point in the borehole ground-penetrating radar image and determine the characteristic apparent resistivity at the target feature point based on the transient electromagnetic distance of the target feature point; the target feature point and the target anomaly point have the same coordinates in the drill pipe advance direction; The filtering module 60 is used to draw a feature circle based on the coordinates of the target feature points in the detection depth direction in the radial section map of the target, and to filter multiple candidate three-dimensional coordinates based on the feature resistivity and the feature circle to obtain the coordinates of the geological anomaly of the ore body to be analyzed; wherein, the number of filtering times is consistent with the number of target feature points selected.

[0067] Specifically, candidate module 40 is also used for: In the radial section view of the target, draw a target circle based on the detection depth coordinates of the target anomaly points; The three-dimensional coordinates of the intersection of the apparent resistivity contour map and the target circle in the borehole transient electromagnetic imaging map are determined as multiple candidate three-dimensional coordinates. The equation of the target circle includes... Multiple candidate 3D coordinates are satisfy , denoted as , where Y represents the coordinates of the drill pipe advance direction, and r0 represents the detection depth coordinates.

[0068] Specifically, module 50 is also used for: Select at least one target feature point in the borehole ground-penetrating radar image; the two-dimensional coordinates of the target feature point. r p To select parameters, r max The preset upper limit threshold; The transient electromagnetic distance of the target feature point is calculated based on its two-dimensional coordinates. The formula for calculating the transient electromagnetic distance of the target feature point includes: , l The spacing between transient electromagnetic and ground-penetrating radars installed in the boreholes of the ore body to be analyzed; Based on the correspondence between the transient electromagnetic distance of the target feature point and the transient electromagnetic distance of the borehole and the apparent resistivity, the characteristic apparent resistivity at the target feature point is determined.

[0069] Specifically, the filtering module 60 is also used for: Starting from the origin of the radial section of the target, connect and extend the coordinates of multiple candidate 3D coordinates until they intersect with the feature circle at multiple corresponding feature intersection points; Based on the coordinates of the feature intersections, calculate the apparent resistivity at each feature intersection. Calculate the error between the apparent resistivity and the feature apparent resistivity at each feature intersection point, and remove the candidate 3D coordinates with errors greater than the preset error value to obtain the filtered 3D coordinates. If the candidate three-dimensional coordinates are unique after screening, then the selected three-dimensional coordinates will be determined as the geological anomaly coordinates of the ore body to be analyzed. If the candidate 3D coordinates after screening are not unique, the other feature points are selected and screened repeatedly until the candidate 3D coordinates after screening are unique. The unique candidate 3D coordinates are then determined as the geological anomaly coordinates of the ore body to be analyzed.

[0070] The present invention also 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 executes the computer program to implement the method provided in the embodiments of the present invention.

[0071] The present invention also provides a computer-readable storage medium storing program code, which can be invoked by a processor to execute the method provided in the embodiments of the present invention.

[0072] Those skilled in the art will recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementations should not be considered beyond the scope of this invention.

[0073] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working processes of the devices, apparatuses, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0074] In the several embodiments provided by this invention, it should be understood that the disclosed devices, apparatuses, and methods can be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative; for instance, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another device, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be through some interfaces; the indirect coupling or communication connection between devices or units may be electrical, mechanical, or other forms.

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

[0076] In addition, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0077] If the aforementioned functions are implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this invention, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0078] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A method for locating anomalies using a combination of borehole transient electromagnetic and ground-penetrating radar, characterized in that, The method includes: Simultaneously collect and preprocess the transient electromagnetic data and ground-penetrating radar data of the ore body to be analyzed from the borehole, and generate the transient electromagnetic imaging map and the ground-penetrating radar imaging map of the borehole respectively. Extract the two-dimensional coordinates of the target anomaly points from the borehole ground-penetrating radar image; the two-dimensional coordinates include the coordinates of the drill rod advance direction and the detection depth. Based on the transient electromagnetic distance of the target anomaly point, the target apparent resistivity corresponding to the target anomaly point is determined; the transient electromagnetic distance is the distance between the transient electromagnetic field and the point to be measured. In the target radial section of the borehole transient electromagnetic imaging image, an apparent resistivity contour map based on the apparent resistivity of the target is drawn, and multiple candidate three-dimensional coordinates of the target anomaly point are determined based on the apparent resistivity contour map and the detection depth coordinates of the target anomaly point; the target radial section is a radial section of the borehole transient electromagnetic imaging image based on the drill rod advance direction coordinates of the target anomaly point; At least one target feature point is selected in the borehole ground-penetrating radar image, and the characteristic apparent resistivity at the target feature point is determined based on the transient electromagnetic distance of the target feature point; the target feature point and the target anomaly point have the same coordinates in the drill pipe advance direction; In the radial section of the target, a feature circle is drawn based on the coordinates of the target feature points in the direction of the detection depth. The multiple candidate three-dimensional coordinates are then filtered based on the feature resistivity and the feature circle to obtain the geological anomaly coordinates of the ore body to be analyzed. The number of filtering times is consistent with the number of target feature points selected.

2. The method according to claim 1, characterized in that, Simultaneously acquire and preprocess borehole transient electromagnetic data and ground-penetrating radar data of the ore body to be analyzed, including: In the tunnel borehole of the ore body to be analyzed, transient electromagnetic and ground-penetrating radar are installed at preset intervals; Based on the transient electromagnetic data and the ground-penetrating radar data, respectively, borehole transient electromagnetic data and ground-penetrating radar data of the ore body to be analyzed are collected and preprocessed.

3. The method according to claim 1, characterized in that, Based on the transient electromagnetic distance of the target anomaly point, the target apparent resistivity at the target anomaly point is determined, including: Based on the transient electromagnetic distance of the target anomaly point and the correspondence between the borehole transient electromagnetic distance and apparent resistivity, the target apparent resistivity at the target anomaly point is determined. The relationship between the transient electromagnetic distance of the borehole and the apparent resistivity includes: In the formula, s is the transient electromagnetic distance during drilling. K For compensation coefficient, t This is the delay time after the transmitting current is turned off. ρ Resistivity μ 0 represents the permeability of free space.

4. The method according to claim 1, characterized in that, Based on the apparent resistivity contour map and the detection depth coordinates of the target anomaly, multiple candidate three-dimensional coordinates of the target anomaly are determined, including: In the radial section view of the target, a target circle is drawn based on the detection depth coordinates of the target anomaly points; The intersection of the apparent resistivity contour map and the target circle is used to determine multiple candidate three-dimensional coordinates in the borehole transient electromagnetic imaging map. The equation of the target circle includes... The multiple candidate three-dimensional coordinates are satisfy , Let be the two-dimensional coordinates of the target anomaly point, Y be the coordinates of the drill pipe advancing direction, and r0 be the coordinates of the detection depth.

5. The method according to claim 4, characterized in that, Selecting at least one target feature point in the borehole ground-penetrating radar image, and determining the characteristic apparent resistivity at the target feature point based on the transient electromagnetic distance of the target feature point, including: At least one target feature point is selected in the borehole ground-penetrating radar image; the two-dimensional coordinates of the target feature point. r p To select parameters, r max The preset upper limit threshold; Based on the two-dimensional coordinates of the target feature point, the transient electromagnetic distance of the target feature point is calculated; the formula for calculating the transient electromagnetic distance of the target feature point includes: , l The spacing between the transient electromagnetic and ground-penetrating radars installed in the boreholes of the ore body to be analyzed; Based on the correspondence between the transient electromagnetic distance of the target feature point and the transient electromagnetic distance of the borehole and the apparent resistivity, the characteristic apparent resistivity at the target feature point is determined.

6. The method according to claim 1, characterized in that, Based on the characteristic apparent resistivity and the characteristic circle, the multiple candidate three-dimensional coordinates are filtered to obtain the geological anomaly coordinates of the ore body to be analyzed, including: Starting from the origin of the radial section of the target, connect and extend the points to the multiple candidate 3D coordinates until they intersect the feature circle at multiple corresponding feature intersection points; Based on the coordinates of the feature intersections, calculate the apparent resistivity at each feature intersection. The error between the apparent resistivity at each feature intersection point and the feature apparent resistivity is calculated, and the candidate three-dimensional coordinates with errors greater than a preset error value are removed to obtain the filtered three-dimensional coordinates. If the candidate three-dimensional coordinates after screening are unique, then the three-dimensional coordinates after screening are determined as the geological anomaly coordinates of the ore body to be analyzed.

7. The method according to claim 6, characterized in that, If the candidate three-dimensional coordinates after the screening are not unique, other feature points are selected and screened repeatedly until the candidate three-dimensional coordinates after screening are unique, and the unique candidate three-dimensional coordinates are determined as the geological anomaly coordinates of the ore body to be analyzed.

8. An anomaly localization system combining borehole transient electromagnetic and ground-penetrating radar, characterized in that, A method for locating anomalies using a combination of borehole transient electromagnetic and ground-penetrating radar as described in any one of claims 1-7; the system comprises: an acquisition module, an extraction module, a determination module, a candidate module, a selection module, and a filtering module; wherein... The acquisition module is used to simultaneously acquire borehole transient electromagnetic data and ground-penetrating radar data of the ore body to be analyzed and preprocess them to generate borehole transient electromagnetic imaging map and borehole ground-penetrating radar imaging map, respectively. The extraction module is used to extract the two-dimensional coordinates of target anomalies in the borehole ground-penetrating radar image; the two-dimensional coordinates include the coordinates of the drill rod advance direction and the detection depth. The determining module is used to determine the target apparent resistivity at the target anomaly point based on the transient electromagnetic distance of the target anomaly point; the transient electromagnetic distance is the distance between the transient electromagnetic field and the point to be measured. The candidate module is used to draw an apparent resistivity contour map based on the apparent resistivity of the target in the target radial section map of the borehole transient electromagnetic imaging image, and to determine multiple candidate three-dimensional coordinates of the target anomaly point based on the apparent resistivity contour map and the detection depth coordinates of the target anomaly point; the target radial section map is a radial section map of the borehole transient electromagnetic imaging image based on the drill rod advance direction coordinates of the target anomaly point; The selection module is used to select at least one target feature point in the borehole ground-penetrating radar image, and determine the characteristic apparent resistivity at the target feature point based on the transient electromagnetic distance of the target feature point; the target feature point and the target anomaly point have the same coordinates in the drill pipe advance direction; The filtering module is used to draw a feature circle based on the coordinates of the target feature points in the detection depth direction in the radial section of the target, and to filter the multiple candidate three-dimensional coordinates based on the feature resistivity and the feature circle to obtain the geological anomaly coordinates of the ore body to be analyzed; wherein, the number of filtering times is consistent with the number of target feature points selected.

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, when executing the computer program, implements the method as claimed in any one of claims 1-7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium contains program code that can be invoked by a processor to execute the method as described in any one of claims 1 to 7.

Images