CN111123264A - Ground penetrating radar reverse time migration imaging method considering medium frequency dispersion and attenuation compensation - Google Patents
Ground penetrating radar reverse time migration imaging method considering medium frequency dispersion and attenuation compensation Download PDFInfo
- Publication number
- CN111123264A CN111123264A CN201911408517.3A CN201911408517A CN111123264A CN 111123264 A CN111123264 A CN 111123264A CN 201911408517 A CN201911408517 A CN 201911408517A CN 111123264 A CN111123264 A CN 111123264A
- Authority
- CN
- China
- Prior art keywords
- transmitting
- electromagnetic wave
- wave field
- zero
- antenna
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 40
- 230000005012 migration Effects 0.000 title claims abstract description 33
- 238000013508 migration Methods 0.000 title claims abstract description 33
- 239000006185 dispersion Substances 0.000 title claims abstract description 29
- 230000000149 penetrating effect Effects 0.000 title claims abstract description 25
- 230000005672 electromagnetic field Effects 0.000 claims abstract description 11
- 230000000704 physical effect Effects 0.000 claims abstract description 11
- 238000000034 method Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 14
- 230000008859 change Effects 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 4
- 230000005284 excitation Effects 0.000 claims description 3
- 230000001902 propagating effect Effects 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims 1
- 230000035699 permeability Effects 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 8
- 239000002609 medium Substances 0.000 description 20
- 230000005540 biological transmission Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 239000004567 concrete Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 239000002612 dispersion medium Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to the technical field of detection, and provides a reverse time migration imaging method of a ground penetrating radar by considering medium frequency dispersion and attenuation compensation, which comprises the following steps: acquiring a transmitting signal of an ith transmitting antenna on a component to be detected and an echo signal of a receiving antenna corresponding to the ith transmitting antenna; constructing an unstructured grid scattering three-dimensional offset velocity model according to the size and physical properties of a component to be detected; calculating forward electromagnetic wave field values S of the transmitting points of all the ith transmitting antennas by using an unstructured grid discrete three-dimensional offset velocity modeli(x1,y1,z1,t1) Counter-propagating electromagnetic field value R of receiving point of corresponding receiving antennai(x2,y2,z2,t2) Zero-delay cross-correlation value ofSum I (x, y, z). The reverse time migration imaging method of the ground penetrating radar with the medium frequency dispersion and attenuation compensation considered, provided by the embodiment of the invention, can effectively compensate the frequency dispersion and attenuation of electromagnetic waves generated when the electromagnetic waves are transmitted in an underground medium, and improves the imaging accuracy.
Description
Technical Field
The invention relates to the technical field of detection, in particular to a reverse time migration imaging method of a ground penetrating radar by considering medium frequency dispersion and attenuation compensation.
Background
The ground penetrating radar technology has been widely applied to the fields of engineering detection, environmental investigation, geological exploration and the like, and with the increasing complexity of detection objects and the continuous improvement of detection requirements, the high-precision ground penetrating radar imaging technology has become an urgent need in practical engineering application. In the working frequency range of the ground penetrating radar, the relative dielectric constant of a detection object such as soil, concrete and rock is a complex function of frequency change, the speed and the attenuation coefficient of electromagnetic waves are changed along with the frequency when the electromagnetic waves propagate in the frequency range, the energy attenuation is fast, and dispersion and distortion are often accompanied. In addition, with increasingly refined ground penetrating radar engineering detection and continuously improved detection frequency, the influence of the frequency dispersion characteristic of the underground medium on electromagnetic wave propagation is continuously increased, and the influence is urgently required to be considered in high-precision imaging of actually measured data.
In order to solve the technical problem, a method for imaging under the condition of a ground penetrating radar reverse time migration velocity model is provided. However, in the prior art, a conventional reverse time migration algorithm is generally adopted, the underground medium is generally assumed to be a non-dispersive medium for the calculation of the source point electromagnetic wave field and the receiving point backward transmission electromagnetic wave field, the dispersion and attenuation characteristics of the source point electromagnetic wave field and the receiving point backward transmission electromagnetic wave field generated when the source point electromagnetic wave field and the receiving point backward transmission electromagnetic wave field propagate in the actual underground medium are ignored, and the actual measurement ground penetrating radar data migration imaging is not clear enough. For the conventional ground penetrating radar reverse time migration imaging algorithm, the schematic diagram is shown in FIG. 1, xTFor transmitting antenna position, xRAnd for the position of a receiving antenna, an S dotted line represents a forward electromagnetic wave field of a source point and an R solid line represents a backward electromagnetic wave field of a receiving point, and the imaging of the target body is realized by calculating the cross correlation between the forward electromagnetic wave field and the backward electromagnetic wave field.
Disclosure of Invention
The embodiment of the invention provides a reverse time migration imaging method of a ground penetrating radar considering medium frequency dispersion and attenuation compensation, which solves the problem that the conventional reverse time migration algorithm is adopted in the prior art, the calculation of a source point electromagnetic wave field and a receiving point reverse transmission electromagnetic wave field usually assumes that an underground medium is a non-frequency dispersion medium, and the data migration imaging of the actually measured ground penetrating radar is not clear enough.
The embodiment of the invention is a reverse time migration imaging method of a ground penetrating radar with medium frequency dispersion and attenuation compensation considered, which comprises the following steps:
acquiring a transmitting signal of an ith transmitting antenna on a component to be detected and an echo signal of a receiving antenna corresponding to the ith transmitting antenna, wherein i is more than or equal to 1, and the transmitting signal comprises transmitting time t1Said echo signal comprising an echo time t2The transmitting antenna is used for transmitting electromagnetic wave signals to the component to be detected, and the receiving antenna is used for receiving echo signals returned by the multi-offset-range ground penetrating radar;
constructing an unstructured grid scattering three-dimensional offset velocity model according to the size and physical properties of the component to be detected;
calculating forward electromagnetic wave field values S of the transmitting points of all the ith transmitting antenna by using the unstructured grid discrete three-dimensional offset velocity modeli(x1,y1,z1,t1) The value R of the counter electromagnetic wave field of the receiving point of the corresponding receiving antennai(x2,y2,z2,t2) Obtaining a reverse time migration imaging result by the sum I (x, y, z) of the zero-delay cross correlation values; wherein,
x1、y1、z1respectively representing the three-dimensional spatial coordinates of the emission points, t1For positively propagating electromagnetic wave time, x2、y2、 z2Three-dimensional spatial coordinates, t, respectively representing the receiving points2For the time of back propagation of electromagnetic waves, x, y and z respectively represent three-dimensional space coordinates of the zero-delay cross-correlation value sum I (x, y and z).
Optionally, the method further includes: and carrying out spatial high-pass filtering on the reverse time migration imaging result to obtain a target reverse time migration imaging result.
Optionally, the step of constructing an unstructured grid scattering three-dimensional offset velocity model according to the size and the physical properties of the component to be measured specifically includes:
constructing a three-dimensional offset speed model according to the size and the physical properties of the component to be detected;
adopting dense mesh subdivision in the region with intense electromagnetic field change in the component to be measured, and dividing the electromagnetic field by using a dense mesh subdivision method
And (4) subdividing the area with gentle change by adopting a sparse grid to obtain an unstructured grid scattering three-dimensional offset speed model.
Optionally, the forward electromagnetic wave field value S of the emission point of all the ith transmitting antenna is calculated by using the unstructured grid scattering three-dimensional offset velocity modeli(x1,y1,z1,t1) The value R of the counter electromagnetic wave field of the receiving point of the corresponding receiving antennai(x2,y2,z2,t2) The step of obtaining the inverse time offset imaging result specifically includes:
calculating forward electromagnetic wave field values S of the transmitting points of all the ith transmitting antenna by using the unstructured grid discrete three-dimensional offset velocity modeli(x1,y1,z1,t1) The value R of the counter electromagnetic wave field of the receiving point of the corresponding receiving antennai(x2,y2,z2,t2) Zero-delay cross-correlation value ofi(x3,y3, z3) The calculation formula is as follows:
Ii(x3,y3,z3)=Si(x1,y1,z1,t1)×Ri(x2,y2,z2,t2);
all the zero-delay cross-correlation values Ii(x3,y3,z3) And summing to obtain a zero-delay cross-correlation value sum I (x, y, z) and further obtain an inverse time migration imaging result, wherein the calculation formula of the zero-delay cross-correlation value sum I (x, y, z) is as follows:
wherein N is the number of transmitting antennas, x3,y3,z3Respectively representing zero-delay cross-correlation values Ii(x3,y3,z3) Three-dimensional space coordinates of (a).
Optionally, the forward electromagnetic wave field value S of the emission point of all the ith transmitting antenna is calculated by using the unstructured grid scattering three-dimensional offset velocity modeli(x1,y1,z1,t1) The value R of the counter electromagnetic wave field of the receiving point of the corresponding receiving antennai(x2,y2,z2,t2) Zero-delay cross-correlation value ofi(x3,y3,z3) The method specifically comprises the following steps:
calculating forward electromagnetic wave field value S of the transmitting point of the ith transmitting antenna by using the unstructured grid scattering three-dimensional offset velocity modeli(x1,y1,z1,t1) The calculation formula is as follows:
calculating the reverse electromagnetic wave field value R of the receiving point of the corresponding receiving antenna by using the unstructured grid discrete three-dimensional offset velocity modeli(x2,y2,z2,t2) The calculation formula is as follows:
wherein,for the generation of frequency dispersion and attenuation terms, epsilon, for the counter-propagating electromagnetic waves0Which is the dielectric constant in a vacuum, is, the relative constants when the frequency tends to be 0 and infinity are respectively, tau is the relaxation time of the medium, J is the excitation source, sigma is the resistivity, and mu is the magnetic conductivity.
The reverse time migration imaging method of the ground penetrating radar with the medium frequency dispersion and attenuation compensation considered, provided by the embodiment of the invention, can effectively compensate the frequency dispersion and attenuation of electromagnetic waves generated when the electromagnetic waves are transmitted in an underground medium, and improves the imaging accuracy.
Drawings
FIG. 1 is a schematic diagram of a conventional reverse time migration imaging algorithm of a ground penetrating radar according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a ground penetrating radar reverse time migration imaging algorithm considering medium frequency dispersion and attenuation according to an embodiment of the present invention;
fig. 3 is a flowchart of a georadar reverse time migration algorithm with medium frequency dispersion and attenuation compensation taken into consideration according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
The reverse time migration imaging method of the ground penetrating radar with the medium frequency dispersion and attenuation compensation considered, provided by the embodiment of the invention, can effectively compensate the frequency dispersion and attenuation of electromagnetic waves generated when the electromagnetic waves are transmitted in an underground medium, and improves the imaging accuracy.
Example one
Referring to fig. 2 and fig. 3, the method for compensating the reverse time migration imaging of the ground penetrating radar by considering the medium dispersion and the attenuation includes the following steps: the first step is as follows: obtaining the transmission signal of the ith transmitting antenna on the component to be tested, and the ith transmitting antennaThe echo signals received by the receiving antenna corresponding to the transmitting antenna, wherein i is more than or equal to 1, and the transmitting signals comprise transmitting time t1The echo signal including an echo time t2The transmitting antenna is used for transmitting electromagnetic wave signals to the component to be measured, and the receiving antenna is used for receiving echo signals returned by the multi-offset-range ground penetrating radar.
In the embodiment of the present invention, it can be understood that a user may set the number of transmitting antennas, that is, the number of receiving antennas corresponding to the number, on a component to be tested according to a requirement, which is not limited specifically.
The second step is that: and constructing an unstructured grid scattering three-dimensional offset velocity model according to the size and physical properties of the component to be measured.
In the embodiment of the present invention, the member to be measured may be concrete, soil, or the like.
As an embodiment of the present invention, the physical properties of the member to be measured include: relative permittivity, conductivity, etc.
In an embodiment of the present invention, the step of constructing the unstructured grid scattering three-dimensional offset velocity model according to the size and the physical property of the component to be measured specifically includes:
constructing a three-dimensional offset speed model according to the size and physical properties of a component to be measured;
and adopting dense mesh subdivision in the region with severe electromagnetic field change in the component to be measured, and adopting sparse mesh subdivision in the region with mild electromagnetic field change to obtain an unstructured mesh dispersion three-dimensional offset speed model.
The third step: calculating forward electromagnetic wave field values S of the transmitting points of all the ith transmitting antennas by using an unstructured grid discrete three-dimensional offset velocity modeli(x1,y1,z1,t1) Counter-propagating electromagnetic field value R of receiving point of corresponding receiving antennai(x2,y2,z2,t2) Obtaining a reverse time migration imaging result by the sum I (x, y, z) of the zero-delay cross correlation values; wherein x is1、y1、z1Respectively representing the three-dimensional spatial coordinates of the emission points, t1For positively propagating electromagnetic wave time, x2、y2、z2Three-dimensional spatial coordinates, t, respectively representing the receiving points2For the time of back propagation of electromagnetic waves, x, y and z respectively represent three-dimensional space coordinates of the zero-delay cross-correlation value sum I (x, y and z).
In the embodiment of the invention, the forward electromagnetic wave field value S of the emission point of all the ith transmitting antennas is calculated by using the unstructured grid discrete three-dimensional offset velocity modeli(x1,y1,z1,t1) Counter-propagating electromagnetic field value R of receiving point of corresponding receiving antennai(x2,y2,z2,t2) The step of obtaining the inverse time offset imaging result specifically includes:
calculating forward electromagnetic wave field values S of the transmitting points of all the ith transmitting antennas by using an unstructured grid discrete three-dimensional offset velocity modeli(x1,y1,z1,t1) Counter-propagating electromagnetic field value R of receiving point of corresponding receiving antennai(x2,y2,z2,t2) Zero-delay cross-correlation value ofi(x3,y3,z3) The calculation formula is as follows:
Ii(x3,y3,z3)=Si(x1,y1,z1,t1)×Ri(x2,y2,z2,t2);
cross-correlating all zero-delay cross-correlation values Ii(x3,y3,z3) And summing to obtain a zero-delay cross-correlation value sum I (x, y, z) and further obtain an inverse time migration imaging result, wherein the calculation formula of the zero-delay cross-correlation value sum I (x, y, z) is as follows:
wherein N is the number of transmitting antennas, x3,y3,z3Respectively representing zero-delay cross-correlation values Ii(x3,y3,z3) Three-dimensional space coordinates of (a).
In the embodiment of the present invention, the sum I (x, y, z) of the zero-delay cross-correlation values is the superimposed data of the zero-delay cross-correlation values of the multiple transmission sources.
In the embodiment of the invention, the forward electromagnetic wave field value S of the emission point of all the ith transmitting antennas is calculated by using the unstructured grid discrete three-dimensional offset velocity modeli(x1,y1,z1,t1) Counter-propagating electromagnetic field value R of receiving point of corresponding receiving antennai(x2,y2,z2,t2) Zero-delay cross-correlation value ofi(x3,y3,z3) The method specifically comprises the following steps:
calculating forward electromagnetic wave field value S of the emitting point of the ith emitting antenna by using the unstructured grid scattering three-dimensional offset velocity modeli(x1,y1,z1,t1) The calculation formula is as follows:
calculating the counter-propagating electromagnetic wave field value R of the receiving point of the corresponding receiving antenna by using the unstructured grid discrete three-dimensional offset velocity modeli(x2,y2,z2,t2) The calculation formula is as follows:
wherein,for the generation of frequency dispersion and attenuation terms, epsilon, for the counter-propagating electromagnetic waves0Which is the dielectric constant in a vacuum, is, the relative constants when the frequency tends to be 0 and infinity are respectively, tau is the relaxation time of the medium, J is the excitation source, sigma is the resistivity, and mu is the magnetic conductivity.
As an embodiment of the present invention, the method further includes, after the third step:
the fourth step: and carrying out spatial high-pass filtering on the reverse time migration imaging result to obtain a target reverse time migration imaging result.
In the embodiment of the present invention, the spatial high-pass filtering uses filter coefficients as follows:
the reverse time migration imaging method of the ground penetrating radar with the medium frequency dispersion and attenuation compensation considered, provided by the embodiment of the invention, can effectively compensate the frequency dispersion and attenuation of electromagnetic waves generated when the electromagnetic waves are transmitted in an underground medium, and improves the imaging accuracy.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.
Claims (5)
1. A method for compensating reverse time migration imaging of a ground penetrating radar by considering medium dispersion and attenuation is characterized by comprising the following steps:
acquiring a transmitting signal of an ith transmitting antenna on a component to be detected and an echo signal of a receiving antenna corresponding to the ith transmitting antenna, wherein i is more than or equal to 1, and the transmitting signal comprises transmitting time t1Said echo signal comprising an echo time t2The transmitting antenna is used for transmitting electromagnetic wave signals to the component to be tested, and the receiving antenna is used for receiving echo signals returned by the multi-offset-range ground penetrating radarNumber;
constructing an unstructured grid scattering three-dimensional offset velocity model according to the size and physical properties of the component to be detected;
calculating forward electromagnetic wave field values S of the transmitting points of all the ith transmitting antenna by using the unstructured grid discrete three-dimensional offset velocity modeli(x1,y1,z1,t1) The value R of the counter electromagnetic wave field of the receiving point of the corresponding receiving antennai(x2,y2,z2,t2) Obtaining a reverse time migration imaging result by the sum I (x, y, z) of the zero-delay cross correlation values; wherein,
x1、y1、z1respectively representing the three-dimensional spatial coordinates of the emission points, t1For positively propagating electromagnetic wave time, x2、y2、z2Three-dimensional spatial coordinates, t, respectively representing the receiving points2For the time of back propagation of electromagnetic waves, x, y and z respectively represent three-dimensional space coordinates of the zero-delay cross-correlation value sum I (x, y and z).
2. The method of claim 1, wherein the method further comprises:
and carrying out spatial high-pass filtering on the reverse time migration imaging result to obtain a target reverse time migration imaging result.
3. The method of claim 1, wherein the step of constructing the unstructured grid discretized three-dimensional offset velocity model based on the dimensions and physical properties of the component to be tested specifically comprises:
constructing a three-dimensional offset speed model according to the size and the physical properties of the component to be detected;
and adopting dense mesh subdivision in the region with severe electromagnetic field change in the component to be tested, and adopting sparse mesh subdivision in the region with mild electromagnetic field change to obtain an unstructured mesh dispersion three-dimensional offset speed model.
4. The method of claim 1, wherein the method comprisesCalculating forward electromagnetic wave field values S of the transmitting points of all the ith transmitting antenna by using the unstructured grid discrete three-dimensional offset velocity modeli(x1,y1,z1,t1) The value R of the counter electromagnetic wave field of the receiving point of the corresponding receiving antennai(x2,y2,z2,t2) The step of obtaining the inverse time offset imaging result specifically includes:
calculating forward electromagnetic wave field values S of the transmitting points of all the ith transmitting antenna by using the unstructured grid discrete three-dimensional offset velocity modeli(x1,y1,z1,t1) The value R of the counter electromagnetic wave field of the receiving point of the corresponding receiving antennai(x2,y2,z2,t2) Zero-delay cross-correlation value ofi(x3,y3,z3) The calculation formula is as follows:
Ii(x3,y3,z3)=Si(x1,y1,z1,t1)×Ri(x2,y2,z2,t2);
all the zero-delay cross-correlation values Ii(x3,y3,z3) And summing to obtain a zero-delay cross-correlation value sum I (x, y, z) and further obtain an inverse time migration imaging result, wherein the calculation formula of the zero-delay cross-correlation value sum I (x, y, z) is as follows:
5. The method of claim 4, wherein said utilizing said non-junctionStructuring a grid scattering three-dimensional offset velocity model, and calculating forward electromagnetic wave field values S of the transmitting points of all the ith transmitting antennasi(x1,y1,z1,t1) The value R of the counter electromagnetic wave field of the receiving point of the corresponding receiving antennai(x2,y2,z2,t2) Zero-delay cross-correlation value ofi(x3,y3,z3) The method specifically comprises the following steps:
calculating forward electromagnetic wave field value S of the transmitting point of the ith transmitting antenna by using the unstructured grid scattering three-dimensional offset velocity modeli(x1,y1,z1,t1) The calculation formula is as follows:
calculating the reverse electromagnetic wave field value R of the receiving point of the corresponding receiving antenna by using the unstructured grid discrete three-dimensional offset velocity modeli(x2,y2,z2,t2) The calculation formula is as follows:
wherein,for the generation of frequency dispersion and attenuation terms, epsilon, for the counter-propagating electromagnetic waves0Which is the dielectric constant in a vacuum, is, are respectively provided withIs the relative constant when the frequency approaches 0 and infinity, tau is the relaxation time of the medium, J is the excitation source, sigma is the resistivity, and mu is the permeability.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911408517.3A CN111123264A (en) | 2019-12-31 | 2019-12-31 | Ground penetrating radar reverse time migration imaging method considering medium frequency dispersion and attenuation compensation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911408517.3A CN111123264A (en) | 2019-12-31 | 2019-12-31 | Ground penetrating radar reverse time migration imaging method considering medium frequency dispersion and attenuation compensation |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111123264A true CN111123264A (en) | 2020-05-08 |
Family
ID=70506196
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911408517.3A Pending CN111123264A (en) | 2019-12-31 | 2019-12-31 | Ground penetrating radar reverse time migration imaging method considering medium frequency dispersion and attenuation compensation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111123264A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111562574A (en) * | 2020-05-22 | 2020-08-21 | 中国科学院空天信息创新研究院 | MIMO ground penetrating radar three-dimensional imaging method based on backward projection |
CN114200526A (en) * | 2020-09-02 | 2022-03-18 | 中国石油天然气股份有限公司 | Electromagnetic exploration method and device |
CN115469376A (en) * | 2022-11-14 | 2022-12-13 | 四川省冶勘设计集团有限公司 | Tunnel advanced geological forecast comprehensive detection method and system based on non-homologous data |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5381151A (en) * | 1994-02-02 | 1995-01-10 | Grumman Aerospace Corporation | Signal processing for ultra-wideband impulse radar |
US20030090405A1 (en) * | 2001-11-02 | 2003-05-15 | Sol Rauch | Spread spectrum radar with leak compensation at baseband |
US20120062408A1 (en) * | 2010-06-15 | 2012-03-15 | Igor Bausov | Unsynchronized radio imaging |
CN105298489A (en) * | 2015-12-03 | 2016-02-03 | 中国石油大学(华东) | Continuous measurement method of dielectric constant frequency dispersion characteristic at wide spectrum of well-close stratum |
CN106707277A (en) * | 2017-01-16 | 2017-05-24 | 中国地质大学(北京) | Method for high-precision reverse time migration imaging based on severe relief surface ground penetrating radar data |
CN109387834A (en) * | 2018-08-28 | 2019-02-26 | 国家能源投资集团有限责任公司 | Groundwater reservoir dam body fault localization method and electronic equipment |
CN109581516A (en) * | 2018-11-29 | 2019-04-05 | 桂林理工大学 | Bent wave zone statistic adaptive threshold Coherent Noise in GPR Record denoising method and system |
CN109917379A (en) * | 2019-01-31 | 2019-06-21 | 上海卫星工程研究所 | Targets in ocean detection method, system, device and medium |
CN110568487A (en) * | 2019-09-19 | 2019-12-13 | 中国科学技术大学 | Active fault structure imaging method based on natural seismic waveform |
-
2019
- 2019-12-31 CN CN201911408517.3A patent/CN111123264A/en active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5381151A (en) * | 1994-02-02 | 1995-01-10 | Grumman Aerospace Corporation | Signal processing for ultra-wideband impulse radar |
US20030090405A1 (en) * | 2001-11-02 | 2003-05-15 | Sol Rauch | Spread spectrum radar with leak compensation at baseband |
US20120062408A1 (en) * | 2010-06-15 | 2012-03-15 | Igor Bausov | Unsynchronized radio imaging |
CN105298489A (en) * | 2015-12-03 | 2016-02-03 | 中国石油大学(华东) | Continuous measurement method of dielectric constant frequency dispersion characteristic at wide spectrum of well-close stratum |
CN106707277A (en) * | 2017-01-16 | 2017-05-24 | 中国地质大学(北京) | Method for high-precision reverse time migration imaging based on severe relief surface ground penetrating radar data |
CN109387834A (en) * | 2018-08-28 | 2019-02-26 | 国家能源投资集团有限责任公司 | Groundwater reservoir dam body fault localization method and electronic equipment |
CN109581516A (en) * | 2018-11-29 | 2019-04-05 | 桂林理工大学 | Bent wave zone statistic adaptive threshold Coherent Noise in GPR Record denoising method and system |
CN109917379A (en) * | 2019-01-31 | 2019-06-21 | 上海卫星工程研究所 | Targets in ocean detection method, system, device and medium |
CN110568487A (en) * | 2019-09-19 | 2019-12-13 | 中国科学技术大学 | Active fault structure imaging method based on natural seismic waveform |
Non-Patent Citations (4)
Title |
---|
"基于归一化互相关成像条件的GPR逆时偏移成像" * |
S. KURODA等: "A Full Waveforminversion Algorithm for Interpreting Crosshole Radar Data" * |
王敏玲等: "基于激发振幅成像条件的探地雷达逆时偏移成像" * |
赵超等: "频散介质探地雷达正演与高速公路实测资料解释" * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111562574A (en) * | 2020-05-22 | 2020-08-21 | 中国科学院空天信息创新研究院 | MIMO ground penetrating radar three-dimensional imaging method based on backward projection |
CN114200526A (en) * | 2020-09-02 | 2022-03-18 | 中国石油天然气股份有限公司 | Electromagnetic exploration method and device |
CN114200526B (en) * | 2020-09-02 | 2023-08-22 | 中国石油天然气股份有限公司 | Electromagnetic exploration method and device |
CN115469376A (en) * | 2022-11-14 | 2022-12-13 | 四川省冶勘设计集团有限公司 | Tunnel advanced geological forecast comprehensive detection method and system based on non-homologous data |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111123264A (en) | Ground penetrating radar reverse time migration imaging method considering medium frequency dispersion and attenuation compensation | |
Cui et al. | Inverse scattering of two-dimensional dielectric objects buried in a lossy earth using the distorted Born iterative method | |
Liu et al. | Frequency-domain reverse-time migration of ground penetrating radar based on layered medium green's functions | |
Sakamoto | A fast algorithm for 3-dimensional imaging with UWB pulse radar systems | |
Shaari et al. | Effects of antenna-target polarization and target-medium dielectric contrast on GPR signal from non-metal pipes using FDTD simulation | |
Zhou et al. | Subsurface cavity imaging by crosshole borehole radar measurements | |
CN108845317A (en) | A kind of frequency domain reverse-time migration algorithm based on layered medium Green's function | |
Zhou et al. | GPR imaging with RM algorithm in layered mediums | |
Xing et al. | The penetrating depth analysis of Lunar Penetrating Radar onboard Chang’e-3 rover | |
Baba et al. | The effect of localized ionospheric perturbations on subionospheric VLF propagation on the basis of finite element method | |
Symes et al. | Solution of an acoustic transmission inverse problem by extended inversion | |
Patsia et al. | Background removal, velocity estimation, and reverse-time migration: a complete GPR processing pipeline based on machine learning | |
Song et al. | High-frequency electromagnetic impedance method for subsurface imaging | |
Zhu et al. | 3D unstructured spectral element method for frequency-domain airborne EM forward modeling based on Coulomb gauge | |
On et al. | Detection of an Object Bottoming at Seabed by the Reflected Signal Modeling | |
Serrano et al. | Alternative cost function for full waveform inversion of GPR data | |
Purwandani et al. | Accuracy improvement in through the wall radar based on deconvolution and delay estimation | |
Yahia et al. | Time delay estimation, using correlation approaches applied to seismic time picking | |
Busch et al. | Combined effective wavelet estimation and full-waveform inversion of GPR data | |
Sandmeier | Non-destructive testing of concrete with electromagnetic and acoustic–elastic waves: data analysis | |
Zhang et al. | A space marching inversion algorithm for pulsed borehole radar in the time-domain | |
Xu et al. | Statistical method to localize buried landmines from GPR array measurement | |
Ge et al. | Time reverse locating of noise source based on isotime point imaging | |
Ebrahimy et al. | FDTD numerical GPR stratigraphy modeling and processing and a case study with GPR data | |
Huyen et al. | Positioning the Adjacent Buried Objects Using UWB Technology Combine with Levenberg Marquardt Algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20200508 |
|
WD01 | Invention patent application deemed withdrawn after publication |