Method of determining planar events from borehole or core images

• US 7,236,887 B2
• Filed: 02/22/2005
• Issued: 06/26/2007
• Est. Priority Date: 02/26/2004
• Status: Active Grant

• Current Assignee: Schlumberger Technology Corporation (SLB)
• Original Assignee: Schlumberger Technology Corporation (SLB)
• Inventors: Etchecopar, Arnaud, Onda, Hitoshi, Mathieu, Gilles, Wang, Yinyu
• Primary Examiner(s): MCELHENY JR, DONALD E
• Application Number:
  US11/063,945
• Publication Number:
  US 20050192753A1
• Time in Patent Office: 854 Days
• Field of Search: 702/6, 702/7, 702/10, 324/339, 324/343
• US Class Current: 702/10
• CPC Class Codes: G01V 9/00 Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00

Abstract

The present invention relates to a user-guided dip picking method from electrical borehole images using Hough Transform. It includes defining a top and bottom reference curve, applying the Hough Transform to the borehole image limited by the top and bottom reference dips so that the Hough Transform space is restrained in the neighborhood of the reference dips. A 1D average trace representing the extrema of the 3D Hough parameter space for each depth is then computed and the extrema exceeding a given threshold are search in this 1D trace. The position in the 3D Hough space of an extremum represents a dip event.

Claims

1. A method for detecting planar events from an image using a Hough transform technique comprising the steps of:

a. Storing a borehole image into a computer based system;

said computer based system performing the following tasks;

b. Invite a human user to specify an upper and a bottom reference curve that border the pattern corresponding to a series of essentially co-planar events to delimit a subset of the 3D Hough accumulator array corresponding to the whole borehole image, said reference curves being used to define a range of reference planar event parameters and a depth interval;

c. Computing the accumulator matrix for said subset;

d. Computing a 1D average trace representing the extrema of the subset of the Hough transform accumulator for each depth;

e. Searching in the 1D average trace the extrema beyond a given threshold;

f. For each extremum in the 1D average trace, getting the horizontal coordinates of the point of the 3D Hough accumulator corresponding to this extremum at depth Zi, said horizontal coordinates representing one planar event to detect; and

g. Detecting the planar events from the image from all the planar events detected from step (f);

h. Storing the detected planar events from the image in the computer based system for retrieval by the human user.

3. The method according to claim 1, wherein said image is a borehole image or a core image from a subterranean well.

4. The method according to claim 1, wherein said geometric pattern are represented by a sinusoid.

5. The method according to claim 1, wherein said reference curves are sinusoids.

6. The method according to claim 1, wherein said reference curves are stored in a database.

7. The method according to claim 1, wherein said reference curves are computed by reference to planar events already identified.

8. The method according to claim 1, wherein said reference curves are provided by the user.

9. The method according to claim 8, wherein said reference curves are traced by the user over the wellbore image.

10. The method according to claim 1, wherein the 1D average trace is obtained by computing H_average, the mean value of H(X0, Y0, Z0) where H(X0, Y0, Z0) is the 3D-accumulator matrix of the Hough transform, centralizing H(X0, Y0, Z0) in the value domain Hc(X0, Y0, Z0)=H(X0, Y0, Z0)−

H_average and, for each depth Z0, finding the maximum of absolute value V(Z0) of Hc.

11. The method according to claim 1, wherein the accumulator matrix of the Hough transform is computed by computing the median value along the sinusoidal features of the image pixel measurements.

2. A method for detecting multiple set of planar events from an image using a Hough transform technique comprising the steps of:

a. Storing a borehole image into a computer based system;

said computer based system performing the following tasks;

b. Invite a human user to specify at least three reference curves to border the multiple pattern corresponding to the multiple set of a series of essentially co-planar events;

c. Sorting the reference curves along depth;

d. For each consecutive pair of reference curves, said reference curves being used to define a range of reference planar event parameters and a depth interval;

 

i. Computing the accumulator matrix for the subset of the 3D Hough accumulator delimited by said reference curves;

 

ii. Computing a 1D average trace representing the extrema of the subset of the Hough transform accumulator for each depth;

 

iii. Searching in the 1D average trace the extrema beyond a given threshold;

 

iv. For each extremum in the 1D average trace, getting the horizontal coordinates of the point of the 3D Hough accumulator corresponding to this extremum at depth Zi, said horizontal coordinates representing one planar event to detect;

 

v. Detecting the planar events from the image from all the planar events detected from step (iv);

 

vi. Storing the detected planar events from the image in the computer based system for retrieval by the human user.

12. The method according to claim 2, wherein said image is a borehole image or a core image from a subterranean well.

13. The method according to claim 2, wherein said geometric pattern are represented by a sinusoid.

14. The method according to claim 2, wherein said reference curves are sinusoids.

15. The method according to claim 2, wherein said reference curves are stored in a database.

16. The method according to claim 2, wherein said reference curves are computed by reference to planar events already identified.

17. The method according to claim 2, wherein said reference curves are provided by the user.

18. The method according to claim 17, wherein said reference curves are traced by the user over the wellbore image.

19. The method according to claim 2, wherein the 1D average trace is obtained by computing H_average, the mean value of H(X0, Y0, Z0) where H(X0, Y0, Z0) is the 3D-accumulator matrix of the Hough transform, centralizing H(X0, Y0, Z0) in the value domain Hc(X0, Y0, Z0)=H(X0, Y0, Z0)−

H_average and, for each depth Z0, finding the maximum of absolute value V(Z0) of Hc.

20. The method according to claim 2, wherein the accumulator matrix of the Hough transform is computed by computing the median value along the sinusoidal features of the image pixel measurements.