Automatic Dip Picking From Wellbore Azimuthal Image Logs

• US 20150098627A1
• Filed: 08/29/2014
• Published: 04/09/2015
• Est. Priority Date: 10/08/2013
• Status: Active Grant

• Current Assignee: Shin-Ju Chu Ye
• Original Assignee: Shin-Ju Chu Ye
• Inventors: Ye, Shin-Ju Chu
• Granted Patent: US 9,366,135 B2
• Time in Patent Office: Days
• US Class Current: 382/109
• CPC Class Codes: E21B 49/00 Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
• CPC Class Codes: G01V 1/50 Analysing data

Abstract

Method for automatically determining a planar event such as a bedding plane from azimuthal wellbore log data (501). The image log may optionally first be filtered using a smoothing filter (502). Then, the image log is inputted to a gradient (first-derivative) filter (503). The image log is also used to find the optimal local sinusoidal trend (or orientation) on the image by determining the best sine wave at each depth among all azimuthal curves from the image using a minimum-variance method (504). Finally, bed boundaries are located based on highest contrast (506) calculated along the sinusoidal lines on the gradient image (505).

Claims

1. A method for determining a bedding surface from an azimuthal image made from a signal provided by a logging tool placed in a cylindrical wellbore, said method comprising:

determining, using a computer and based on minimum variance, a best sine wave in the image at each of a plurality of selected depths;

computing contrast at each selected depth along the best sine wave in the image; and

locating one or more bedding surface boundaries at one or more depths with highest contrasts.

2. The method of claim 1, further comprising filtering the image with a gradient filter before computing contrast.

3. The method of claim 1, wherein determining a best sine wave at a selected depth comprises:

locating a depth window in the image centered at the selected depth, then dividing the window into a plurality of discrete depth increments;

at each depth increment, computing a variance of values in the image along sine waves for a plurality of different amplitude and phase values for the sine waves;

for each combination of amplitude and phase values, summing the variances of all sine waves within the depth window; and

using the amplitude and phase values giving minimum summed variance to determine the best sine wave.

4. The method of claim 1, further comprising applying a smoothing filter to the image before the determining a best sine wave.

5. The method of claim 1, wherein the variance is computed using image data at each azimuthal bin derived from the logging tool signal.

6. The method of claim 5, further comprising using interpolated data between available azimuthal bins to provide more variance information.

7. The method of claim 1, wherein the contrast is computed for image data at each azimuthal bin provided by the logging tool, then summed over azimuthal bins.

8. The method of claim 7, further comprising using interpolated data between available azimuthal bins to provide more contrast information.

9. The method of claim 1, wherein the contrast is based directly on the gradient, or is based upon a figure of merit that depends upon the gradient.

10. The method of claim 1, wherein the azimuthal logging tool is one of a group consisting of a gamma-ray tool, a density tool, a neutron tool, an acoustic tool, and a resistivity tool.

11. The method of claim 1, wherein the method is fully automated, i.e. all actions in the method are performed using a computer.

12. A non-transitory computer usable medium having a computer readable program code embodied therein, said computer readable program code adapted to be executed to implement a method for determining a bedding surface from an azimuthal image made from a signal provided by a logging tool placed in a cylindrical wellbore, said method comprising:

determining, based on minimum variance, a best sine wave in the image at each of a plurality of selected depths;

computing contrast at each selected depth along the best sine wave in the image; and

locating one or more bedding surface boundaries at one or more depths with highest contrasts.

13. The non-transitory computer usable medium of claim 12, wherein said method further comprises filtering the image with a gradient filter before computing contrast.

14. The non-transitory computer usable medium of claim 12, wherein determining a best sine wave at a selected depth comprises:

locating a depth window in the image centered at the selected depth, then dividing the window into a plurality of discrete depth increments;

at each depth increment, computing a variance of values in the image along sine waves for a plurality of different amplitude and phase values for the sine waves;

for each combination of amplitude and phase values, summing the variances of all sine waves within the depth window; and

using the amplitude and phase values giving minimum summed variance to determine the best sine wave.

15. The non-transitory computer usable medium of claim 12, wherein said method further comprises applying a smoothing filter to the image before the determining a best sine wave.