Method for reconstructing complex wave attributes from limited view measurements

• US 7,330,792 B2
• Filed: 09/11/2006
• Issued: 02/12/2008
• Est. Priority Date: 07/28/2003
• Status: Expired due to Fees

• Current Assignee: Alan Witten
• Original Assignee: Alan Witten
• Inventors: Witten, Alan
• Primary Examiner(s): McElheny, Jr.; Donald E
• Application Number:
  US11/518,821
• Publication Number:
  US 20070035438A1
• Time in Patent Office: 519 Days
• Field of Search: 702/16, 702/17
• US Class Current: 702/16
• CPC Class Codes: G01N 29/0672 by acoustic tomography medical tomography A61B8/13
• CPC Class Codes: G01N 29/07 by measuring propagation velocity or propagation time of acoustic waves
• CPC Class Codes: G01N 29/11 by measuring attenuation of acoustic waves
• CPC Class Codes: G01N 29/46 by spectral analysis, e.g. Fourier analysis or wavelet analysis spectral signal processing per se G06F17/14
• CPC Class Codes: G01V 1/28 Processing seismic data, e.g. for interpretation or for event detection G01V1/48 takes precedence
• CPC Class Codes: G06F 17/141 Discrete Fourier transforms
• CPC Class Codes: G06F 17/17 Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method G06F17/18 takes precedence ; interpolation for numerical control G05B19/18

Abstract

A method is disclosed for reconstructing complex wave attributes from limited view measurements of a scattering object. The method involves the analytic continuation of the Fourier transform of the object function into the area in which there is an absence of K-space coverage by requiring objects to be an even function. (It is assumed that physical objects are even functions, and it is this assumption that allows analytic continuation.) When the object function is not centered at the origin, the measurements are shifted to the origin prior to determining the analytic continuation and returned to their original location following analytic continuation.

Claims

1. An apparatus for reconstructing and displaying wave speed and attenuation described by an object function O from limited view ground penetrating radar measurements of a subterranean object with associated wavevector K, comprising:

a ground penetrating radar device capable of acquiring said limited view ground penetrating radar measurements of said subterranean object;

a processor that computes Fourier transformed measurements from said limited view ground penetrating radar measurements, determines a Fourier transformed object function Õ

from said object function O, detennines an analytic relationship between said Fourier transformed object function Õ and

said Fourier transformed measurements, analytically extends said Fourier transform Õ

by specifying that Õ

(K)=Õ

(−

K), and reconstructs said wave speed and attenuation by inverting said analytically extended Fourier transform of Õ

; and

an output device that displays said reconstructed wave speed and attenuation.

2. The apparatus of claim 1 wherein said ground penetrating radar device uses a reflection measurement geometry where a single transmitting antenna and a single receiving antenna with a fixed separation are moved in unison above said subterranean object.

3. The apparatus of claim 1 wherein said output device displays said wave speed and displays said attenuation so that said displayed wave speed and said displayed attenuation have identical shapes and different contrasts.

4. An apparatus for reconstructing and displaying wave speed and attenuation described by an object function O from limited view ground penetrating radar measurements of a subterranean object with associated wavevector K:

a ground penetrating radar device capable of acquiring said limited view ground penetrating radar measurements of said subterranean object;

a processor that computes Fourier transformed measurements from said limited view ground penetrating radar measurements, determines a midpoint of said object, creates shifted Fourier transformed measurements by shifting said Fourier transformed measurements so that said midpoint is located at the origin, determines an analytic relationship between said object function O and said shifted Fourier transformed measurements, determines the Fourier transformed object function Õ

of said object function O from said Fourier transformed measurements using said analytic relationship, analytically extends said Fourier transform Õ

by specifying that Õ

(K)=Õ

(−

K), determines shifted wave speed and attenuation by inverting said analytically extended Fourier transform of Õ

, and reconstructs said wave speed and attenuation by shifting said shifted wave speed and attenuation back to said midpoint; and

an output device that displays said reconstructed wave speed and attenuation.

5. The apparatus of claim 4 wherein said ground penetrating radar device uses a reflection measurement geometry where a single transmitting antenna and a single receiving antenna with a fixed separation are moved in unison above said subterranean object.

6. The apparatus of claim 4 wherein said output device displays said wave speed and displays said attenuation so that said displayed wave speed and said displayed attenuation have identical shapes and different contrasts.

7. An apparatus for reconstructing and displaying wave speed and attenuation described by an object function O from limited view seismic measurements of a geologic structure with associated wavevector K, comprising:

a seismic device capable of acquiring said limited view seismic measurements of said geologic structure;

a processor that computes Fourier transformed measurements from said limited view seismic measurements, determines a Fourier transformed object function Õ

from said object function O, detennines an analytic relationship between said Fourier transformed object function Õ and

said Fourier transformed measurements, analytically extends said Fourier transform Õ

by specifying that Õ

(K)=Õ

(−

K), and reconstructs said wave speed and attenuation by inverting said analytically extended Fourier transform of Õ

; and

an output device that displays said reconstructed wave speed and attenuation.

8. The apparatus of claim 7 wherein said output device displays said wave speed and displays said attenuation so that said displayed wave speed and said displayed attenuation have identical shapes and different contrasts.

9. An apparatus for reconstructing and displaying wave speed and attenuation described by an object function O from limited view seismic measurements of a geologic structure with associated wavevector K:

a seismic device capable of acquiring said limited view seismic measurements of said geologic structure;

a processor that computes Fourier transformed measurements from said limited view seismic measurements, determines a midpoint of said object, creates shifted Fourier transformed measurements by shifting said Fourier transformed measurements so that said midpoint is located at the origin, determines an analytic relationship between said object function O and said shifted Fourier transformed measurements, determines the Fourier transformed object function Õ

of said object function O from said Fourier transformed measurements using said analytic relationship, analytically extends said Fourier transform Õ

by specifying that Õ

(K)=Õ

(−

K), determines shifted wave speed and attenuation by inverting said analytically extended Fourier transform of Õ

, and reconstructs said wave speed and attenuation by shifting said shifted wave speed and attenuation back to said midpoint; and

an output device that displays said reconstructed wave speed and attenuation.

10. The apparatus of claim 9 wherein said output device displays said wave speed and displays said attenuation so that said displayed wave speed and said displayed attenuation have identical shapes and different contrasts.

11. An apparatus for reconstructing and displaying wave speed and attenuation described by an object function O from limited view seismic measurements of a geologic structure with associated wavevector K:

a seismic device capable of acquiring said limited view seismic measurements of said geologic structure;

a processor that defines said geologic structure as a linear superposition of a plurality of layers, mirrors said limited view seismic measurements for each of said plurality of layers so that each has a common midpoint, computes Fourier transformed measurements for each of said plurality of layers from said limited view seismic measurements, creates shifted Fourier transformed measurements for each of said plurality of layers by shifting said Fourier transformed measurements so that said midpoint is located at the origin, determines an analytic relationship between said object function O and said shifted Fourier transformed measurements for each of said plurality of layers, determines the Fourier transformed object function Õ

of said object function O from said Fourier transformed measurements using said analytic relationship for each of said plurality of layers, analytically extends said Fourier transform Õ

by specifying that Õ

(K)=Õ

(−

K) for each of said plurality of layers, determines shifted wave speed and attenuation for each of said plurality of layers by inverting said analytically extended Fourier transform of Õ

, and reconstructs said wave speed and attenuation by shifting said shifted wave speed and attenuation for each of said plurality of layers back to said midpoint; and

an output device that displays said reconstructed wave speed and attenuation.

12. The apparatus of claim 11 wherein said output device displays said wave speed and displays said attenuation so that said displayed wave speed and said displayed attenuation have identical shapes and different contrasts.