Claims
- 1. A method of determining a parameter of interest of subsurface formations containing a sand and a shale surrounding a borehole, the method comprising:(a) conveying an electromagnetic logging tool into the borehole and using at least one electromagnetic transmitter and one electromagnetic receiver on the tool to obtain measurements related to a horizontal and vertical resistivity of the formation; (b) obtaining a plurality of layer intervals for the subsurface; (c) deriving a horizontal and vertical resistivity at a plurality of depths from said measurements using an anisotropic whole space model; (d) obtaining from the derived horizontal and vertical resistivities a layered model; (e) deriving a shoulder bed corrections based upon the layered model and the whole space model; (f) applying the shoulder bed corrections to said measurements to obtain shoulder bed corrected measurements; (g) determining a difference between said shoulder bed corrected measurements at the center of each said plurality of layers to a synthetic response of the whole space model at said centers; and (h) iteratively repeating steps (c)-(g) to give an updated layered model.
- 2. The method of claim 1 wherein said electromagnetic logging tool is conveyed one one of (i) a wireline, (ii) a drillstring, and (iii) coiled tubing.
- 3. The method of claim 1 wherein the parameter of interest is at least one of (i) a horizontal resistivity of one of said layers, (ii) a vertical resistivity of one of said layers, (iii) a fluid saturation of one of said layers.
- 4. The method of claim 1 wherein deriving said horizontal and vertical resistivity at said plurality of depths further comprises performing a whole space inversion.
- 5. The method of claim 4 wherein performing a whole space inversion further comprises computing a synthetic tool response.
- 6. The method of claim 4 wherein performing a whole space inversion further comprises updating estimates of said horizontal and vertical resistivity using a Lanczos procedure.
- 7. The method of claim 1 wherein obtaining said plurality of layer intervals further comprises determining a first and a second derivative with respect to depth of said measurements.
- 8. The method of claim 1 wherein obtaining said plurality of layer intervals further comprises:(i) Fourier transforming said signals into the spatial frequency domain; (ii) low pass filtering said Fourier transformed signals at a cutoff about equal to an axial resolution of said induction signals; (iii) calculating a central first derivative of said filtered Fourier transformed signals; (iv) calculating an inverse Fourier transform of said central first derivative; (v) selecting roots of said inverse Fourier transformed central first derivative; and (vi) testing localized properties of said inverse Fourier transformed central first derivative within a selected number of data sample points of said roots, thereby providing indications of formation layer boundaries at axial positions most likely to be true ones of said formation layer boundaries.
- 9. The method of claim 1 wherein deriving shoulder bed corrections further comprises:(i) determining a 1-D layered synthetic response is calculated from the layered model using a synthetic tool response; (ii) determining a whole space response using said horizontal and vertical resistivities; and (iii) obtaining said shoulder bed correction as a difference between said 1-D layered synthetic response and said whole space response.
- 10. The method of claim 1 wherein step (h) is carried out until said difference is less than a predetermined threshold.
- 11. The method of claim 1 wherein step (h) is carried out a prespecified number of iterations.
CROSS-REFERENCES TO RELATED APPLICATIONS
This application claims priority from U.S. Provisional Patent Application Ser. No. 60/164,055 filed on Nov. 8, 1999.
US Referenced Citations (9)
Number |
Name |
Date |
Kind |
5446654 |
Chemali et al. |
Aug 1995 |
A |
5656930 |
Hagiwara |
Aug 1997 |
A |
5854991 |
Gupta et al. |
Dec 1998 |
A |
5867806 |
Strickland et al. |
Feb 1999 |
A |
5883515 |
Strack et al. |
Mar 1999 |
A |
5999883 |
Gupta et al. |
Dec 1999 |
A |
5999884 |
Kriegshauser et al. |
Dec 1999 |
A |
6098019 |
Hakvoort et al. |
Aug 2000 |
A |
6304086 |
Minerbo et al. |
Oct 2001 |
B1 |
Non-Patent Literature Citations (1)
Entry |
Gerald W. Hohmann et al.; Inversion of Controlled-Source Electromagnetic Data, Investigations in Geophysics, Society of Exploration Geophysicists, vol. 1, Theory, 1987, pp. 469-503. |
Provisional Applications (1)
|
Number |
Date |
Country |
|
60/164055 |
Nov 1999 |
US |