Diffraction tomography systems and methods with fixed detector arrays

• US 4,594,662 A
• Filed: 09/19/1983
• Issued: 06/10/1986
• Est. Priority Date: 11/12/1982
• Status: Expired due to Term

• Current Assignee: Schlumberger Technology Corporation (SLB)
• Original Assignee: Schlumberger Technology Corporation (SLB)
• Inventors: Devaney, Anthony J.
• Primary Examiner(s): NOT, DEFINED
• Assistant Examiner(s): Jablon, Clark A.
• Application Number:
  US06/533,391
• Time in Patent Office: 995 Days
• Field of Search: 364/414, 364/420, 364/421, 364/422, 364/400, 73/602, 250/253, 250/256, 250/269, 128/660, 367/14, 378/4, 378/901
• US Class Current: 700/90
• CPC Class Codes: G01S 15/8977 using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
• CPC Class Codes: G06T 11/006 Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
• CPC Class Codes: G06T 2211/421 Filtered back projection [FBP]
• CPC Class Codes: Y10S 128/916 Ultrasound 3-D imaging
• CPC Class Codes: Y10S 378/901 Computer tomography program or processor

Abstract

Diffraction tomography systems and methods are disclosed for the reconstruction of physical properties of two- and three-dimensional objects from information collected by detector arrays which are fixed in space. Improved preprocessing procedures and an optimal filtered backpropagation technique are disclosed for both geophysical and ultrasonic diffraction tomography systems and methods. A novel method for reconstructing properties of three-dimensional objects and slices thereof in real time is also disclosed. An optimal filtered backpropagation technique for the backscatter mode is further disclosed.

Claims

1. A method for the tomographic reconstruction of a subsurface geological formation which has diffracted at least one wave of energy interacting with it, comprising the steps of:

(a) obtaining one or more signals which are a function of at least one of the phase and amplitude of the diffracted propagating wave after said wave has interacted with said formation;

(b) preprocessing said obtained signals, wherein said preprocessing includes the step of scaling and phase correcting signals; and

(c) converting said preprocessed signals by means of a filtered backpropagation technique into an array representing the partial reconstruction profile of the formation.

2. A method according to claim 1 further comprising:

(d) displaying said reconstruction profile of the formation.

3. A method according to claim 1 wherein said preprocessing step further comprises:

converting said one or more obtained signals from analog to digital form;

transforming said digital signal over time with a Fourier transform;

normalizing said transformed signals;

complex phase converting said normalized signals;

phase unwrapping said complex phase converted signals for input into said scaling and phase correcting step.

4. A method for the tomographic reconstruction of a subsurface geological formation comprising:

(a) directing at a phase tilt one or more waves of energy toward said geological formation;

(b) detecting as a function of time, the one or more fields produced by said one or more waves;

(c) producing one or more signals which are functions of said one or more detected fields;

(d) preprocessing said produced signals wherein said preprocessing includes the step of scaling and phase correcting signals; and

(e) filtering said scaled and phase corrected preprocessed signals according to a filtered backpropagation technique to generate a two-dimensional array representing the partial reconstruction of said geological formation.

5. A method according to claim 4 further comprising:

(f) interpolating said partial reconstruction onto a master array;

(g) changing said phase tilt; and

(h) repeating steps a-g at different view angles to obtain a two-dimensional reconstruction of said formation.

6. A method according to claim 5 further comprising:

(i) displaying said two-dimensional reconstruction of said formation.

7. A method according to claim 5 wherein said preprocessing step further comprises:

converting said one or more produced signals from analog to digital form;

transforming said digital signal over time with a Fourier transform;

normalizing said transformed signals;

complex phase converting said normalized signals;

phase unwrapping said complex phase converted signals for input into said scaling and phase correcting step.

8. A method according to claim 4 wherein:

said preprocessed signals are filtered and backpropagated according to a filtered backpropagation technique

space="preserve" listing-type="equation">H.sub.q,m =h.sub.m e.sup.iks.sbsp.o.sup..r e.sup.iγ

.sbsp.m.sup.(z.sbsp.q-.sup.l.sbsp.o.sup.)

wherein H_q,m is a filtered backpropagation technique for object depth z_q, and m indexes sample values of preprocessed signals;

h_m is a filter which equals the absolute value of ##EQU15## where M is the number of detectors detecting said fields;

k is the wavenumber associated with the dominant frequency of said waves;

z_q is the distance of the investigated object depth q from the axis of rotation of said object; and

.sub. o is the distance between the axis of rotation of said object and the detector array. s₀.r is a scalar product between the unit propagation vector s₀ and the position vector r; and

θ

_n is the phase tilt of said directed waves of energy.

9. A method according to claim 8 wherein:

said filter h_m is changed to equal ##EQU16## when -π

≦

θ

_o ≦

-π

/2, wherein S is a step function.

10. A method according to claim 9 wherein:

said directing and detecting steps are accomplished by sources and receivers located along an identical line.

11. A method according to claim 4 wherein:

said directing and said detecting steps are accomplished by sources and receivers located along an identical line.

12. A method for the reconstruction of a three-dimensional object which has diffracted waves of energy propagating through it at at least one view angle on a plane for at least one planar orientation, comprising the steps of:

(a) obtaining signals which are a function of at least one of the phase and amplitude of each propagating wave after said wave has interacted with said object;

(b) preprocessing and separating said obtained signals to obtain sets of preprocessed signals, each set representing the signals obtained from the propagation of an energy wave at a separate view angle, said preprocessing and separating including the step of scaling and phase correcting signals.

(c) converting said sets of preprocessed signals by means of a filtered backpropagation technique into arrays representing the partial reconstruction of the the object; and

(d) interpolating each partial reconstruction onto a master array to obtain a two-dimensional reconstruction of a planar projection of the object.

13. A method according to claim 12 further comprising:

(e) additionally separating said obtained signals into second sets of signals, each second set of signals representing the signals obtained from the propagation of an energy wave at a separate planar orientation, whereby by said converting step, the two-dimensional planar projected reconstruction is for each planar orientation, and a three-dimensional reconstruction of the object can be obtained from a plurality of said two-dimensional planar projection reconstructions.

14. A method for tomographic reconstruction of a three-dimensional object comprising:

(a) directing at a view angle on a plane one or more waves of energy toward an object being examined;

(b) detecting as a function of time the one or more fields produced by said one or more waves;

(c) producing one or more signals which are functions of said one or more detected fields;

(d) preprocessing said produced signals wherein said preprocessing includes the step of scaling and phase correcting signals; and

(e) filtering said preprocessed signals according to a filtered backpropagation technique to generate a two-dimensional array representing the partial reconstruction of the object.

15. A method according to claim 14 further comprising:

(f) adding said partial reconstruction array onto a master array;

(g) changing said view angle on said plane; and

(h) repeating steps a-g at different view angles on said plane to obtain a two-dimensional reconstruction of a planar projection of the object.

16. A method according to claim 15 further comprising:

(i) changing said plane on which said waves of energy are directed at a view angle; and

(j) repeating steps a-i at a sufficient number of different planar orientations to obtain a plurality of planar projections of the three-dimensional object profile from which a full three-dimensional reconstruction of the object profile can be obtained.

17. A method according to claim 16 further comprising:

displaying said planar projection of the three-dimensional reconstruction of the object.

18. A method according to claim 14 wherein said preprocessing step further comprises:

converting said produced signal from analog to digital form;

transforming said digital signals over time with a Fourier transform;

normalizing said transformed signals;

complex phase converting said normalized signals;

phase unwrapping said complex phase converted signals; and

projecting a two-dimensional array of complex phase converted signals onto the line formed by the intersection of the detecting plane with the plane of rotation of the directing-detecting system to produce a one-dimensional array for input into said filtering step.

19. A method according to claim 18 wherein:

said directing and said detecting steps are accomplished by two-dimensional arrays of sources and receivers located in the same plane.

20. A method according to claim 14 wherein

said preprocessed signals are filtered and backpropagated according to a filtered backpropagation technique

21. A method for tomographic reconstruction of a three-dimensional object comprising:

(a) directing at a view angle on a plane one or more waves of energy toward an object being examined;

(b) detecting as a function of time the one or more fields produced by said one or more waves;

(c) producing one or more signals which are functions of said one or more detected fields;

(d) preprocessing said produced signals wherein said preprocessing includes,

converting said produced signal from analog to digital form;

transforming said digital signals over time with a Fourier transform;

normalizing said transformed signals;

complex phase converting said normalized signals;

phase unwrapping said complex phase converted signals scaling and phase correcting said phase unwrapped signals;

(e) filtering said preprocessed signals according to a filtered backpropagation technique;

(f) changing said view angle on said plane;

(g) changing said plane on which said waves of energy are directed at a view angle; and

(h) repeating steps a-g at a sufficient number of different planar orientations such that a full three-dimensional reconstruction of the object profile can be obtained.

22. A method according to claim 21 wherein:

said preprocessed signals are filtered and backpropagated according to a filtered backpropagation technique

space="preserve" listing-type="equation">H.sub.q,m,n =h.sub.m,n (S.sub.o)e.sup.-iks.sbsp.o.sup..r e.sup.iγ

.sbsp.m,n.sup.(z-l.sbsp.o.sup.)

wherein H_q,m,n is a filtered backpropagation technique for object depth z_q, and n and m index sample values of preprocessed signals;

h_m,n is a deblurring filter;

##EQU18## where M and N are the number of detector rows and columns used in detecting said wave;

k is the wavenumber associated with the frequency of said directed waves;

z_q is the distance of the investigated object depth q from the axis of rotation of said object; and

.sub. o is the distance between the axis of rotation of said object and the detector array.

s₀ r is a scalar product between the unit propagation vector s₀ =sin θ

₀ cos φ

₀ sin θ

₀ sin φ

₀ cos θ

₀) and the position vector r, wherein s₀ has two degrees of freedom;

θ

_n is the view angle; and

φ

₀ is the planar orientation.

23. A method according to claim 22 wherein:

a planar section of said object is reconstructed at a given depth by setting z_q of said filtered backpropagation technique equal to said given depth.

24. A method according to claim 21 wherein:

said directing and said detecting steps are accomplished by two-dimensional arrays of sources and receivers located in the same plane.

25. A system for the tomographic reconstruction of a subsurface geological formation which has diffracted at least one wave of energy interacting with it comprising:

(a) obtaining means for obtaining one or more signals which are a function of at least one of the phase and amplitude of the diffracted propagating wave after said wave has interacted with said formation;

(b) preprocessing means for preprocessing said signals obtained by said obtaining means, wherein said preprocessing means includes means for scaling and phase correcting signals; and

(c) filtering means for converting said preprocessed signals by means of a filtered backpropagation technique into an array representing the partial reconstruction profile of the formation.

26. A system for the tomographic reconstruction of a subsurface geological formation comprising:

source means for directing at a phase-tilt one or more waves of energy toward said geological formation;

(b) detecting means for detecting as a function of time, the one or more fields produced by said one or more waves and for producing one or more signals which are functions of said one or more fields wherein said detecting means is located on the same line as said source means;

(c) preprocessing means for preprocessing said produced signals wherein said preprocessing means includes means for scaling and phase correcting signals; and

(d) filtering means for filtering said scaled and phase corrected preprocessed signals according to a filtered backpropagation technique to generate a two-dimensional array representing the partial reconstruction of said geological formation.

27. A system for tomographic reconstruction of a three-dimensional object comprising:

(a) source means for directing at a view angle on a plane one or more waves of energy directed toward an object being examined;

(b) detecting means for detecting as a function of time the one or more fields produced by said one or more waves, and producing one or more signals which are functions of said one or more fields;

(c) converting means for correcting said produced signals from analog to digital form;

(d) transforming means for Fourier transforming said digital signals over time;

(e) normalizing means for normalizing said transformed signals;

(f) complex phase converting means for complex phase converting said normalized signals;

(g) phase unwrapping means for phase unwrapping said complex phase converted signals to produce a two-dimensional array;

(h) scaling and phase converting means for scaling and phase converting said phase unwrapped signals;

(i) filtering means for filtering said scaled and phase converted signals according to a filtered backpropagation technique

space="preserve" listing-type="equation">H.sub.q,m,n =h.sup.m,n (S.sub.o)e.sub.-iks o.sup..r e.sup.iγ

.sbsp.m,n.sup.(z-l.sbsp.o.sup.)

wherein H_q,m,n is a filtered backpropagation technique for object depth z_q, and n and m index sample values of preprocessed signals;

h_m,n is a deblurring filter;

##EQU19## where M and N are the number of detector rows and columns used in detecting said wave;

k is the wavenumber associated with the frequency of said directed waves;

z_q is the distance of the investigated object depth q from the axis of rotation of said object; and

.sub. o is the distance between the axis of rotation of said object and the detector array.

s₀.r is a scalar product between the unit propagation vector s₀ =sin θ

₀ cos φ

₀ sin θ

₀ sin φ

₀ cos θ

₀) and the position vector r;

θ

_n is the view angle; and

φ

_n is the planar orientation.

28. A system according to claim 27 wherein:

said source means and detecting means are co-planar.