Sub-surface imaging of dielectric structures and voids via narrowband electromagnetic resonance scattering

Patent metadata

Field	Value
Publication number	US-10976461-B2
Application number	US-201816162197-A
Country	US
Kind code	B2
Filing date	Oct 16, 2018
Priority date	Oct 17, 2017
Publication date	Apr 13, 2021
Grant date	Apr 13, 2021

Abstract

Subsurface imaging with information of shape, volume, and dielectric properties is achieved with low frequencies and a ramp waveform. The low frequencies have a lower attenuation compared to the penetration losses of radar frequencies. The technique operates at wavelengths which are comparable to the object or void being imaged, and can be applied to detect and image underground aquifers, magma chambers, man-made tunnels and other underground structures.

First claim

What is claimed is:

1. A method comprising:
generating electromagnetic waves in a wavelength range, the electromagnetic waves having a plurality of orthogonal linear polarizations;
transmitting, by a transmitter, the electromagnetic waves to an object to be imaged; and
detecting, by a receiver and at each polarization of the plurality of orthogonal linear polarizations, electromagnetic waves scattered by the object,
wherein:
the object is located below ground and is embedded in a lossy dielectric material,
the transmitter and the receiver are located above ground, and
the wavelength range comprises wavelengths between 0.2 and 3 times a perimeter or lateral dimension of the object,
the method further comprising:
applying a constant threshold to determine if the object is present or not; and
determining the constant threshold based on thermal noise and scattering losses of the lossy dielectric material.

2. The method of claim 1 , wherein the electromagnetic waves scattered by the object are resonant electromagnetic waves.

3. The method of claim 1 , further comprising calculating resonant wavelengths of the object, and wherein the generating the electromagnetic waves comprises generating the electromagnetic waves at the resonant wavelengths.

4. The method of claim 1 , wherein the object is void.

5. The method of claim 1 , further comprising taking a ratio of two polarizations of the plurality of orthogonal linear polarizations, thereby cancelling attenuation losses.

6. The method of claim 1 , wherein transmitting the electromagnetic waves comprises sweeping a wavelength over the wavelength range.

7. The method of claim 6 , wherein the sweeping the wavelength is linear.

8. The method of claim 7 , wherein the sweeping the wavelength is based on a wavelength step size, and a frequency step size corresponding to the wavelength step size is smaller than a speed of light divided by two times a depth of the object from ground level.

9. The method of claim 1 , further comprising applying a ramp filter in a frequency domain to obtain a spectral domain ramp response.

10. The method of claim 9 , further comprising multiplying the spectral domain ramp response by (jω) −2 , where j is an imaginary unit, and ω is a radial frequency.

11. The method of claim 10 , further comprising:
applying an inverse Fourier transform to the spectral domain ramp response, thereby obtaining a time domain ramp response; and
calculating a surface area for the object as a function of depth.

12. The method of claim 11 , wherein the applying the inverse Fourier transform and the calculating the surface area are for each polarization of the plurality of orthogonal linear polarizations.

13. The method of claim 1 , wherein the transmitting the electromagnetic waves is carried out at a plurality of angles, the angles being between a range axis and a unit cross-product of two polarizations of the plurality of orthogonal linear polarizations scattered by the object.

14. The method of claim 2 , further comprising detecting a bulk relative dielectric permittivity of the object by measuring higher order modes of the resonant electromagnetic waves.

15. The method of claim 14 , wherein the detecting the bulk relative dielectric permittivity comprises calculating a simulated response by the object, and minimizing an error between the simulated response and the electromagnetic waves scattered by the object.

16. The method of claim 1 , further comprising detecting a shape, location, dielectric permittivity, and dielectric conductivity of the body, based on the electromagnetic waves scattered by the object.

17. The method of claim 1 , wherein the wavelength range corresponds to frequencies below 1 MHz.

18. A method comprising:
generating electromagnetic waves in a wavelength range, the electromagnetic waves having a plurality of orthogonal linear polarizations;
transmitting, by a transmitter, the electromagnetic waves to an object to be imaged; and
detecting, by a receiver and at each polarization of the plurality of orthogonal linear polarizations, electromagnetic waves scattered by the object,
wherein:
the object is located below ground,
the transmitter and the receiver are located above ground,

the wavelength range comprises wavelengths between 0.2 and 3 times a perimeter or lateral dimension of the object,
transmitting the electromagnetic waves comprises sweeping a wavelength over the wavelength range, and
the sweeping the wavelength is linear and based on a wavelength step size, a frequency step size corresponding to the wavelength step size being smaller than a speed of light divided by two times a depth of the object from ground level.

19. A method comprising:
generating electromagnetic waves in a wavelength range, the electromagnetic waves having a plurality of orthogonal linear polarizations;
transmitting, by a transmitter, the electromagnetic waves to an object to be imaged; and
detecting, by a receiver and at each polarization of the plurality of orthogonal linear polarizations, electromagnetic waves scattered by the object,
wherein:
the object is located below ground,
the transmitter and the receiver are located above ground, and

the wavelength range comprises wavelengths between 0.2 and 3 times a perimeter or lateral dimension of the object,
the method further comprising:
applying a ramp filter in a frequency domain to obtain a spectral domain ramp response;
multiplying the spectral domain ramp response by (jω) −2 , where j is an imaginary unit, and ω is a radial frequency;
applying an inverse Fourier transform to the spectral domain ramp response, thereby obtaining a time domain ramp response; and
calculating a surface area for the object as a function of depth.

20. The method of claim 19 , wherein the applying the inverse Fourier transform and the calculating the surface area are for each polarization of the plurality of orthogonal linear polarizations.

Assignees

• California Inst Of Techn

Inventors

• Arumugam Darmindra D

Classifications

• G01S13/90

  using synthetic aperture techniques {, e.g. synthetic aperture radar [SAR] techniques} · CPC title

• G01S7/025

  involving the transmission of linearly polarised waves · CPC title

• G01S13/885

  for ground probing (prospecting or detecting using electromagnetic waves G01V3/12) · CPC title

• G01V3/12 Primary

  operating with electromagnetic waves {(operating with millimetre waves G01V8/005)} · CPC title

• G01N2021/4792

  Polarisation of scatter light · CPC title