A high-precision cross-hole geological radar detection method

Abstract

The application discloses a high-precision cross-hole geological radar detection method, comprising the following steps: acquiring radar single-hole self-sending and self-receiving data of a single excitation source of two detection holes respectively, and acquiring radar double-hole transmission data of the two detection holes; generating an initial geological model by using the radar single-hole self-sending and self-receiving data; converting the radar double-hole transmission data into radar single-hole reflection wave data by virtual source mapping; carrying out forward modeling on the radar double-hole transmission data based on the initial geological model to obtain a forward wave field; carrying out wave field reverse continuation on the radar single-hole reflection wave data based on the initial geological model to obtain a reverse-time wave field; imaging the forward wave field and the reverse-time wave field at the same time to obtain imaging results, and then superimposing the imaging results of the excitation sources at different depths of the same detection hole to obtain wave field migration imaging. Compared with the prior art, the method improves the imaging precision.

Description

Technical Field

[0001] This invention relates to the field of geophysical exploration technology, specifically to a high-precision cross-hole ground-penetrating radar exploration method. Background Technology

[0002] Cross-hole ground-penetrating radar (GPR) is a high-resolution geophysical exploration method widely used in engineering geological exploration, underground structural defect detection, karst geological body identification, and tunnel advanced geological prediction. Its basic principle involves placing transmitting and receiving antennas in two parallel boreholes. By studying the propagation characteristics of electromagnetic waves in the underground medium (such as travel time, amplitude, and waveform), the electrical parameter distribution of the medium between the boreholes is inverted, thereby achieving detailed imaging of the geological structure.

[0003] However, existing cross-hole radar detection systems typically employ a "one-transmit, one-receive" working mode: only the transmitting antenna is placed in one borehole, and only the receiving antenna is placed in another borehole. This mode has the following technical defects: (1) Low data acquisition efficiency: Each detection can only acquire transmission data from one transmitting position to multiple receiving positions, requiring multiple movements of the transmitting antenna to complete the scanning of the entire profile. (2) Single detection perspective: It can only acquire electromagnetic wave information penetrating the medium from one direction, lacking multi-angle constraints, resulting in limited uniqueness and resolution of the inversion results. (3) Inability to acquire reflection information: Traditional transmission imaging mainly utilizes the travel time of the first arrival wave, making it difficult to fully utilize the fine structural information carried by reflected waves, scattered waves, etc. (4) Insufficient equipment utilization: The radar probe in each borehole has a single function; the transmitting probe can only transmit, and the receiving probe can only receive, resulting in a waste of hardware resources. Summary of the Invention

[0004] The purpose of this invention is to overcome the above-mentioned technical deficiencies and propose a high-precision cross-hole ground-penetrating radar detection method to solve the technical problem of insufficient detection accuracy in the existing "one-transmitter-one-receiver" detection mode.

[0005] To achieve the above-mentioned technical objectives, the present invention adopts the following technical solution: In a first aspect, the present invention provides a high-precision cross-hole ground-penetrating radar detection method, comprising the following steps: The radar single-aperture spontaneous transmission and reception data and the radar dual-aperture transmission data of the two detector apertures were acquired separately from the single excitation source of the two detector apertures. An initial geological model is generated based on the radar single-hole self-transmitted and self-received data. The radar dual-aperture transmission data is converted into radar single-aperture reflection wave data through virtual source mapping. Based on the initial geological model, a forward modeling simulation was performed on the radar dual-aperture transmission data to obtain the forward propagation wavefield. Based on the initial geological model, the wavefield is reverse-extended on the radar single-hole reflection wave data to obtain the inverse time wavefield. The forward propagation wave field and the reverse propagation wave field at the same time are imaged. After obtaining the imaging results, the imaging results of excitation sources at different depths of the same probe hole are superimposed to obtain wave field migration imaging.

[0006] In some embodiments, acquiring radar single-aperture self-transmitted and self-received data from a single excitation source of each of the two detection apertures includes: The reflection signals of the excitation source at different depths of the first detection aperture are collected respectively to obtain the first radar single aperture self-transmission and self-reception data of the excitation source. The reflected signals of the excitation source at different depths of the second detection aperture are collected to obtain the second radar single aperture self-transmission and self-reception data of the excitation source.

[0007] In some embodiments, the method for acquiring the radar dual-aperture transmission data of the two detection apertures includes: When the excitation source is located at a preset depth of the first detection hole, transmission signals at different depths within the second detection hole are collected to obtain the radar dual-hole transmission data.

[0008] In some embodiments, in the virtual source mapping, the virtual reflected wave field between points a and b in the same probe aperture is: R(a,b,t)= , Where p(a,s,t) and p(b,s,t) represent the wave field information at points a and b, respectively, s represents the position of the virtual excitation source, and t represents time.

[0009] In some embodiments, the step of performing forward modeling on the radar dual-aperture transmission data based on the initial geological model to obtain the forward propagation wavefield includes: Numerical simulations were conducted based on the initial geological model to calculate the wave field of the electromagnetic wave emitted by the excitation source propagating in the positive direction of the time axis. The wave field information at each time step was stored and denoted as S(χ,t,s). i ); χ represents the location of the receiving source, t represents time, s i This indicates the position of the i-th excitation source.

[0010] In some embodiments, the step of performing wavefield reverse extrapolation on the radar single-aperture reflection wave data based on the initial geological model to obtain a reverse-time wavefield includes: The radar single-aperture reflected wave data is used as the excitation source to propagate the wave field in reverse along the time axis, and the wave field information at each time step is stored and denoted as R(χ,t,s). i ).

[0011] In some embodiments, the wavefield migration imaging is denoted as .

[0012] Secondly, the present invention also provides a high-precision cross-hole ground-penetrating radar detection device, comprising: The data acquisition module acquires radar single-aperture spontaneous transmission and reception data from a single excitation source of the two detection apertures, and radar dual-aperture transmission data from the two detection apertures. The model generation module uses the radar single-hole self-transmitted and self-received data to generate an initial geological model; The virtual source mapping module converts the radar dual-aperture transmission data into radar single-aperture reflection wave data through virtual source mapping. The forward modeling module performs forward modeling on the radar dual-aperture transmission data based on the initial geological model to obtain the forward propagation wavefield. The reverse continuation module performs wavefield reverse continuation on the radar single-aperture reflection wave data based on the initial geological model to obtain the inverse time wavefield. The offset imaging module images the forward propagation wave field and the reverse propagation wave field at the same time to obtain the imaging results. Then, the imaging results of excitation sources at different depths of the same probe hole are superimposed to obtain wave field offset imaging.

[0013] Thirdly, the present invention also provides an electronic device, including a memory and a processor, wherein, The memory is used to store programs; The processor, coupled to the memory, is used to execute the program stored in the memory to implement the steps in the high-precision cross-hole ground-penetrating radar detection method described above.

[0014] Fourthly, the present invention also provides a readable storage medium for storing a computer-readable program or instruction, which, when executed by a processor, enables the implementation of the steps in the high-precision cross-hole ground-penetrating radar detection method described above.

[0015] Compared with existing technologies, the high-precision cross-hole ground-penetrating radar detection method provided by this invention improves data acquisition efficiency, expands the detection angle, and improves azimuth resolution through single-hole directional self-transmitting and self-receiving mode and dual-hole joint transmission mode; it also improves data utilization, thereby improving the detection accuracy of external anomalies; and electromagnetic wave imaging is not limited by tilt angle, improving multi-wave field utilization and imaging accuracy. Attached Figure Description

[0016] Figure 1 This is a flowchart of the high-precision cross-hole ground-penetrating radar detection method provided in the embodiments of the present invention; Figure 2 This is a schematic diagram of the structure of the radar probe provided in an embodiment of the present invention; Figure 3 A schematic diagram of a cross-hole ground-penetrating radar directional detection device according to another embodiment of the present disclosure is shown; Figure 4 This is a schematic block diagram of an electronic device provided in an embodiment of the present invention. Detailed Implementation

[0017] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the invention.

[0018] To address the technical problem of insufficient detection accuracy in existing detection methods, this invention provides a high-precision cross-hole ground-penetrating radar detection method that can achieve high-precision detection.

[0019] Please see Figure 1 , Figure 1 This is a flowchart of a high-precision cross-hole ground-penetrating radar (GPR) detection method provided in an embodiment of the present invention. The high-precision cross-hole GPR detection method includes the following steps: S101. Acquire radar single-aperture spontaneous transmission and reception data from a single excitation source of each of the two probe apertures, and acquire radar dual-aperture transmission data from the two probe apertures.

[0020] S102. Generate an initial geological model using the radar single-hole self-transmitted and self-received data.

[0021] S103. Convert the radar dual-aperture transmission data into radar single-aperture reflection wave data through virtual source mapping.

[0022] S104. Based on the initial geological model, perform forward modeling on the radar dual-aperture transmission data to obtain the forward propagation wavefield.

[0023] S105. Based on the initial geological model, the radar single-aperture reflection wave data is subjected to wavefield reverse extension to obtain the inverse time wavefield.

[0024] S106. Image the forward propagation wave field and the reverse propagation wave field at the same time to obtain the imaging results. Then, superimpose the imaging results of excitation sources at different depths of the same probe hole to obtain wave field migration imaging.

[0025] In some embodiments, in step S1, a radar probe is used to collect data; this radar probe has both transmitting and receiving functions. For example... Figure 2 As shown, Figure 2This is a schematic diagram of the radar probe provided in an embodiment of the present invention. The radar probe includes a battery 1, a storage module 2, a gyroscope 3, a main control module 4, a low-frequency antenna 5, a high-frequency antenna 6, a housing 7, and a rubber anti-collision head 8. The battery 1, storage module 2, gyroscope 3, main control module 4, low-frequency antenna 5, and high-frequency antenna 6 are all housed within the housing 7. The battery 1 is connected to the storage module 2, gyroscope 3, main control module 4, low-frequency antenna 5, and high-frequency antenna 6, respectively, providing them with power. The storage module 2, gyroscope 3, low-frequency antenna 5, and high-frequency antenna 6 are respectively connected to the main control module 4 via optical fiber communication, exchanging and transmitting information with the main control module 4. The rubber anti-collision head 8 is located at one end of the housing 7 and serves as a shock absorber when the radar probe is placed into the detection hole.

[0026] In this specific embodiment, battery 1 uses a high-power battery for continuous power supply, and storage module 2 uses a large-capacity internal memory card, enabling long-term wireless data acquisition. Gyroscope 3 uses a high-precision gyroscope to record the device's orientation and synchronizes with low-frequency antenna 5 and high-frequency antenna 6 to achieve electromagnetic wave reception orientation recording.

[0027] Data is acquired using a combination of low-frequency antenna 5 and high-frequency antenna 6. One set of high-frequency antennas 6 is used for fine-grained detection, primarily targeting minute anomalies or structural details, while one set of low-frequency antennas 5 is used for long-range detection, mainly to extend the detection range and detect large anomalies. Each set of low-frequency antennas 5 or high-frequency antennas 6 has four identical butterfly antennas, which can emit electromagnetic waves in four directions radially outward from the radar. The rubber anti-collision head 8 is made of high-density rubber material to ensure that the equipment is not damaged during mobile measurement.

[0028] This radar probe has a single-aperture self-transmitting and self-receiving mode and a dual-aperture combined mode.

[0029] Single-hole self-transmitting and self-receiving mode: The main control module 4 controls the low-frequency antenna 5 and the high-frequency antenna 6 to transmit electromagnetic wave signals in four directions of the detection hole through optical fiber communication. The electromagnetic waves propagate in the strata outside the detection hole, and after being reflected back by the geological body outside the detection hole, they are received by the low-frequency antenna 5 and the high-frequency antenna 6. The main control module 4 synchronizes with the gyroscope through optical fiber communication and records the electromagnetic wave signal and azimuth angle θ and time t at the receiving point χ at this time, which is recorded as p(χ,t,θ) and transmitted to the storage module 2 for storage through optical fiber communication.

[0030] Dual-aperture combined mode: The main control module 4 of the first radar probe controls the low-frequency antenna 5 and the high-frequency antenna 6 to generate electromagnetic waves, while the main control module 4 of the second radar probe controls the low-frequency antenna 5 and the high-frequency antenna 6 to receive electromagnetic waves. The electromagnetic waves propagate in the strata outside the probe hole, are transmitted through the geological body outside the probe hole, and are received by the second radar probe. The main control module 4 of the second radar probe controls the storage module 2 to record the received electromagnetic wave signals through time synchronization.

[0031] Therefore, in step S101, the step of acquiring radar single-aperture spontaneous transmission and reception data of a single excitation source from each of the two detection apertures includes acquiring the reflection signals of the excitation source at different depths in the first detection aperture to obtain the first radar single-aperture spontaneous transmission and reception data of the excitation source. Similarly, acquiring the reflection signals of the excitation source at different depths in the second detection aperture to obtain the second radar single-aperture spontaneous transmission and reception data of the excitation source. The acquisition method involves setting the radar probe to single-aperture spontaneous transmission and reception mode, placing it in the detection aperture, and acquiring data sequentially from top to bottom or bottom to top at preset distance intervals. The data acquired from the two detection apertures are recorded as p1(χ,t,θ) and p2(χ,t,θ), respectively.

[0032] In step S101, the step of acquiring radar dual-aperture transmission data of the two detection apertures includes setting the radar probe to dual-aperture combined mode, placing the first radar probe as the excitation source into the first detection aperture at a preset depth, and using the second radar probe as the receiving source to collect transmission signals at different depths of the second detection aperture, so as to obtain radar dual-aperture transmission data of the excitation source at that depth.

[0033] In a preferred embodiment, the first radar probe can be placed at the same depth within the second detection aperture as an excitation source, and the second radar probe can be used as a receiving source to collect transmission signals at different depths within the first detection aperture, thus obtaining another set of radar dual-aperture transmission data at that depth. Superimposing the two sets of data yields the final radar dual-aperture transmission data, which improves data acquisition accuracy.

[0034] In some embodiments, in step S102, based on the collected single-aperture radar data p1(χ,t,θ) and p2(χ,t,θ), the approximate distribution of geological anomalies around the detection hole can be determined through the single-aperture radar self-transmission and self-reception profile. The location of the geological anomalies can be determined by the correspondence between the angles recorded by the gyroscopes 3 inside the two radar probes and the electromagnetic wave signals, thereby initially determining the initial model M0.

[0035] In reverse time migration imaging algorithms, the initial model is a key factor affecting imaging performance. To obtain better imaging results, the initial model should be as close to the real model as possible. If the initial model differs significantly from the real model, the imaging effect may even be worse than using a uniform initial model. Therefore, during data acquisition, a single-hole self-transmitting and self-receiving mode is first used to obtain the basic conditions of the hole perimeter, pile body, and pile perimeter, thereby determining the initial model M0.

[0036] Since reverse time migration (RTM) provides the best imaging effect for reflected waves, and since both reflected and transmitted waves exist simultaneously in cross-hole radar (CT) detection, its application in CT imaging generally only yields images of the interface perpendicular to the borehole, not those parallel to it. This severely limits the application of RTM in borehole radar imaging. Therefore, it is necessary to perform a mirror transformation on the CT data to make it equivalent to reflected wave borehole imaging. Specifically, step S103 converts the radar dual-aperture transmission data into radar single-aperture reflected wave data through virtual source mapping.

[0037] Assume there is an excitation source c in the first probe hole, and two receiving points a and b in the second probe hole. Numerous diffractors exist between the first and second probe holes. Let p(a,c,t) and p(b,c,t) represent the wavefield records received at points a and b respectively after excitation by source c. To improve imaging quality, we use virtual source mapping to convert the wavefield records into a record emitted at point b, reflected by geological bodies, and then received at point b, denoted as R(a,c,t).

[0038] The wave field records at points a and b can be represented as follows: p(a,c,t)=e(t)*g ac (t) and p(b,c,t)=e(t)*g bc (t).

[0039] Where e(t) is the time excitation function of the excitation source, t represents time, and g ac (t) and g bc (t) is the Green's function at points a and b when excited at point c, and represents the impulse response between points ac and bc, respectively; * indicates convolution operation.

[0040] The cross-correlation between the wave field records at points a and b can be expressed as:

[0041] The above equation shows that the impulse response between points a and b is contained in g. ab In (t), g ab (t) represents the Green's function at point a when the excitation source is located at point b, and the actual response R(a,b,t) and Green's function g. ab The difference between (t) and (t) is only e(t)*e(-t).

[0042] When there are many excitation sources:

[0043] In this context, the first term on the left, the delta function δ(ab)δ(t), has almost no effect, and the second term R(a,b,-t) is only meaningful when t<0. Therefore, R(a,b,t)=-p s (a,t)*p s (b,-t).

[0044] The above equation shows that the wave field with point b as the virtual excitation source and point a as the receiving source, as well as the wave field with point a as the virtual excitation source and point b as the receiving source, can both be obtained by cross-correlation of the receiving records at points a and b. Therefore, the reflected wave field with point b as the virtual excitation source and point a as the receiving source in the second probe aperture can be obtained by the propagation response between the excitation source c in the first probe aperture and the two receiving points a and b in the second probe aperture.

[0045]

[0046] The above equation shows that the virtual reflected wave field R(a,b,t) between points a and b in the same borehole can be obtained by cross-correlation calculation of the wave field records p(a,s,t) and p(b,s,-t) of the two points and summing them over all source points.

[0047] In some embodiments, step S104 is based on the initial geological model M O The steps for performing forward modeling on radar dual-aperture transmission data to obtain the forward propagation wavefield include, based on the initial geological model M... O Numerical simulations were conducted to calculate the wave field of the electromagnetic wave emitted by the excitation source propagating along the positive time axis, and the wave field information at each time step was stored and denoted as S(χ,t,s). i ); χ represents the location of the receiving source, t represents time, s i This indicates the position of the i-th excitation source.

[0048] Based on the above embodiments, step S105 involves performing a reverse wavefield extrapolation on the radar single-aperture reflection wave data based on the initial geological model to obtain a reverse-time wavefield. This includes using the radar single-aperture reflection wave data as an excitation source to propagate the wavefield in reverse along the time axis, and storing the wavefield information for each time step as R(χ,t,s). i ).

[0049] S106. Imaging the forward-propagating wavefield and the reverse-propagating wavefield at the same time to obtain the imaging result. This imaging result is calculated based on the radar dual-aperture transmission data acquired under the set depth conditions in step S101. By modifying the depth conditions of the excitation source when acquiring the radar dual-aperture transmission data in step S101, imaging results under other depth conditions can be obtained. Similarly, imaging results under different depth conditions for the same probe aperture are acquired sequentially, and the various imaging results are superimposed to obtain the desired wavefield migration imaging. This is denoted as... .

[0050] Please see Figure 3 , Figure 3 A schematic diagram of a cross-hole ground-penetrating radar directional detection device according to another embodiment of the present disclosure is shown. The high-precision cross-hole ground-penetrating radar detection device 300 includes a data acquisition module 310, a model generation module 320, a virtual source mapping module 330, a forward modeling module 340, a reverse continuation module 350, and a migration imaging module 360.

[0051] The data acquisition module 310 acquires radar single-aperture self-transmitting and self-receiving data from a single excitation source of the two detection apertures, and radar dual-aperture transmission data from the two detection apertures.

[0052] The model generation module 320 uses radar single-hole self-transmitted and self-received data to generate an initial geological model.

[0053] The virtual source mapping module 330 converts radar dual-aperture transmission data into radar single-aperture reflection wave data through virtual source mapping.

[0054] The forward modeling module 340 performs forward modeling on radar dual-aperture transmission data based on the initial geological model to obtain the forward propagation wavefield.

[0055] The reverse continuation module 350 performs wavefield reverse continuation on radar single-hole reflection wave data based on the initial geological model to obtain the reverse time wavefield.

[0056] The migration imaging module 360 ​​images the forward propagation wave field and the reverse propagation wave field at the same time to obtain the imaging results. Then, the imaging results of excitation sources at different depths of the same probe hole are superimposed to obtain wave field migration imaging.

[0057] The aforementioned high-precision cross-hole ground-penetrating radar detection device 300 can be implemented as a computer program, which can be used in, for example... Figure 4 It runs on the electronic device shown. Figure 4 This is a schematic block diagram of an electronic device provided in an embodiment of the present invention. The electronic device 500 is a host computer or a server.

[0058] The electronic device 500 includes a processor 502, a memory, and a network interface 505 connected via a device bus 501, wherein the memory may include a storage medium 503 and internal memory 504.

[0059] The storage medium 503 may store an operating system 5031 and a computer program 5032. When the computer program 5032 is executed, it enables the processor 502 to perform a high-precision cross-hole ground-penetrating radar detection method.

[0060] The processor 502 provides computing and control capabilities to support the operation of the entire electronic device 500.

[0061] The internal memory 504 provides an environment for the operation of the computer program 5032 in the storage medium 503. When the computer program 5032 is executed by the processor 502, the processor 502 can perform a high-precision cross-hole ground-penetrating radar detection method.

[0062] This network interface 505 is used for network communication, such as providing data transmission. Those skilled in the art will understand that... Figure 4 The structure shown is merely a block diagram of a portion of the structure related to the present invention and does not constitute a limitation on the electronic device 500 to which the present invention is applied. The specific electronic device 500 may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0063] The processor 502 is used to run a computer program 5032 stored in a memory to implement the high-precision cross-hole ground-penetrating radar detection method disclosed in the embodiments of the present invention.

[0064] Those skilled in the art will understand that Figure 4 The embodiments of the computer device shown do not constitute a limitation on the specific configuration of the computer device. In other embodiments, the computer device may include more or fewer components than shown, or combine certain components, or have different component arrangements. For example, in some embodiments, the computer device may include only memory and a processor. In such embodiments, the structure and function of the memory and processor are consistent with those described above, and will not be repeated here.

[0065] It should be understood that, in this embodiment of the invention, the processor 502 may be a Central Processing Unit (CPU), or it may be other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or any conventional processor.

[0066] In another embodiment of the present invention, a computer-readable storage medium is provided. This computer-readable storage medium may be a non-volatile computer-readable storage medium or a volatile computer-readable storage medium. The computer-readable storage medium stores a computer program, wherein when executed by a processor, the computer program implements the high-precision cross-hole ground-penetrating radar detection method disclosed in the embodiments of the present invention.

[0067] Those skilled in the art will readily 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. 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, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in terms of function in the foregoing description. 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 implementation should not be considered beyond the scope of this invention.

[0068] In the 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. Units with the same function may be grouped into one unit. For example, multiple units or components may be combined or integrated into another device, or some features may be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed may be indirect coupling or communication connection through some interfaces, devices, or units, or may be electrical, mechanical, or other forms of connection.

[0069] 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 the embodiments of the present invention, depending on actual needs.

[0070] Furthermore, 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. The integrated unit can be implemented in hardware or as a software functional unit.

[0071] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or 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 an electronic device (which may be a personal computer, a backend server, or a network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present 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), magnetic disks, or optical disks.

[0072] The specific embodiments of the present invention described above do not constitute a limitation on the scope of protection of the present invention. Any other corresponding changes and modifications made in accordance with the technical concept of the present invention should be included within the scope of protection of the claims of the present invention.

Claims

1. A high-precision cross-hole ground-penetrating radar detection method, characterized in that, The steps include the following: The radar single-aperture spontaneous transmission and reception data and the radar dual-aperture transmission data of the two detector apertures were acquired separately from the single excitation source of the two detector apertures. An initial geological model is generated based on the radar single-hole self-transmitted and self-received data. The radar dual-aperture transmission data is converted into radar single-aperture reflection wave data through virtual source mapping. Based on the initial geological model, a forward modeling simulation was performed on the radar dual-aperture transmission data to obtain the forward propagation wavefield. Based on the initial geological model, the wavefield is reverse-extended on the radar single-hole reflection wave data to obtain the inverse time wavefield. The forward propagation wave field and the reverse propagation wave field at the same time are imaged. After obtaining the imaging results, the imaging results of excitation sources at different depths of the same probe hole are superimposed to obtain wave field migration imaging.

2. The high-precision cross-hole ground-penetrating radar detection method according to claim 1, characterized in that, The acquisition of radar single-aperture self-transmitted and self-received data from a single excitation source of each of the two detection apertures includes: The reflection signals of the excitation source at different depths of the first detection aperture are collected respectively to obtain the first radar single aperture self-transmission and self-reception data of the excitation source. The reflected signals of the excitation source at different depths of the second detection aperture are collected to obtain the second radar single aperture self-transmission and self-reception data of the excitation source.

3. The high-precision cross-hole ground-penetrating radar detection method according to claim 2, characterized in that, The method for acquiring the radar dual-aperture transmission data of the two detection apertures includes: When the excitation source is located at a preset depth of the first detection hole, transmission signals at different depths within the second detection hole are collected to obtain the radar dual-hole transmission data.

4. The high-precision cross-hole ground-penetrating radar detection method according to claim 3, characterized in that, In the virtual source mapping, the virtual reflected wave field between points a and b in the same detection aperture is: R(a,b,t)= , Where p(a,s,t) and p(b,s,t) represent the wave field information at points a and b, respectively, s represents the location of the virtual excitation source, and t represents time.

5. The high-precision cross-hole ground-penetrating radar detection method according to claim 3, characterized in that, The forward modeling of the radar dual-aperture transmission data based on the initial geological model yields the forward propagation wavefield, including: Numerical simulations were conducted based on the initial geological model to calculate the wave field of the electromagnetic wave emitted by the excitation source propagating in the positive direction of the time axis. The wave field information at each time step was stored and denoted as S(χ,t,s). i ); χ represents the location of the receiving source, t represents time, s i This indicates the position of the i-th excitation source.

6. The high-precision cross-hole ground-penetrating radar detection method according to claim 5, characterized in that, The process of performing wavefield reverse extrapolation on the radar single-aperture reflection wave data based on the initial geological model to obtain the inverse time wavefield includes: The radar single-aperture reflected wave data is used as the excitation source to propagate the wave field in reverse along the time axis, and the wave field information at each time step is stored and denoted as R(χ,t,s). i ).

7. The high-precision cross-hole ground-penetrating radar detection method according to claim 6, characterized in that, The wavefield migration imaging is denoted as .

8. A high-precision cross-hole ground-penetrating radar detection device, characterized in that, include: The data acquisition module acquires radar single-aperture spontaneous transmission and reception data from a single excitation source of the two detection apertures, and radar dual-aperture transmission data from the two detection apertures. The model generation module uses the radar single-hole self-transmitted and self-received data to generate an initial geological model; The virtual source mapping module converts the radar dual-aperture transmission data into radar single-aperture reflection wave data through virtual source mapping. The forward modeling module performs forward modeling on the radar dual-aperture transmission data based on the initial geological model to obtain the forward propagation wavefield. The reverse continuation module performs wavefield reverse continuation on the radar single-aperture reflection wave data based on the initial geological model to obtain the inverse time wavefield. The offset imaging module images the forward propagation wave field and the reverse propagation wave field at the same time to obtain the imaging results. Then, the imaging results of excitation sources at different depths of the same probe hole are superimposed to obtain wave field offset imaging.

9. An electronic device, characterized in that, Including memory and processor, among which, The memory is used to store programs; The processor, coupled to the memory, is used to execute the program stored in the memory to implement the steps in the high-precision cross-hole ground-penetrating radar detection method according to any one of claims 1 to 7.

10. A readable storage medium, characterized in that, Used to store computer-readable programs or instructions, which, when executed by a processor, can implement the steps in the high-precision cross-hole ground-penetrating radar detection method according to any one of claims 1 to 7.

Images