AUTOMATIC METHOD FOR THREE-DIMENSIONAL STRUCTURAL INTERPRETATION OF BOREHOLE IMAGES ACQUIRED IN HIGH-ANGLE AND HORIZONTAL WELLS

• US 20160370480A1
• Filed: 02/27/2015
• Published: 12/22/2016
• Est. Priority Date: 02/28/2014
• Status: Active Grant

• Current Assignee: Schlumberger Technology Corporation (Schlumberger NV)
• Original Assignee: Schlumberger Technology Corporation (Schlumberger NV)
• Inventors: Picco, Vittorio, Shetty, Sushil, Rasmus, John, Morriss, Christopher Edward, Ito, Koji, Asif, Shahzad
• Granted Patent: US 10,466,375 B2
• Time in Patent Office: Days
• US Class Current: 1/1
• CPC Class Codes: G01V 1/30 Analysis G01V1/50 takes precedence
• CPC Class Codes: G01V 1/34 Displaying seismic recordings or visualisation of seismic data or attributes
• CPC Class Codes: G01V 1/50 Analysing data
• CPC Class Codes: G01V 3/18 specially adapted for well-logging
• CPC Class Codes: G01V 5/04 specially adapted for well-logging

Abstract

A method for automatic interpretation of bulls-eye and sinusoidal features observed in LWD images is disclosed. In some embodiments, the method includes an automatic workflow for extracting smooth contours from images that demarcate boundaries of structural features, followed by projection of the contours to three-dimensional (3D) point clouds in the well coordinate system for structural interpretation. The method may characterize both sinusoidal features and bulls-eye features, taking into account variations of formation dip/azimuth, or well inclination/azimuth, on the topology of a structural feature. The disclosed method may be sufficiently fast for use in real-time analysis and interpretation, or to provide constraints for physics-based data inversion processing.

Claims

1. A method comprising:

accepting as input to a processor an azimuthally substantially continuous wellbore image obtained using a well logging tool moved along a wellbore penetrating a subsurface formation;

in the processor, extracting contours from the substantially azimuthally continuous image;

in the processor, grouping the extracted contours into clusters corresponding to a single transition zone; and

in the processor, mapping the clustered contours having a measured depth extent that is greater than a length-scale over which a dip of the subsurface formation varies to a three-dimensional space corresponding to a coordinate system associated with the well; and

in the processor, calculating a dip and azimuth of a layer boundary of the subsurface formation from the mapped clustered contours and at least one of storing and displaying the calculated dip and azimuth.

2. The method of claim 1, wherein the contours comprise at least one of open contours, closes contours and combinations thereof.

3. The method of claim 1 further comprising, in the processor, prior to grouping the extracted contours into the clusters, filtering the extracted contours to remove spurious contours resulting from the extraction.

4. The method of claim 3, wherein the filtering comprises computing a Fourier series approximation for each extracted contour, and filtering out contours for which a quality of fit parameter is less than a selected threshold.

5. The method of claim 1 further comprising, in the processor, for extracted contours having a measured depth extent that is less than the length-scale over which the dip of the subterranean formation varies, estimating relative dip and apparent azimuth based on a first harmonic approximation of a contour.

6. The method of claim 1, wherein the well logging tool comprises at least one of a resistivity tool, a density tool and a photoelectric factor tool.

7. The method of claim 1, wherein the extracting the contours is performed using at least one of a marching-squares algorithm, a square tracing algorithm, a Moore-Neighbor algorithm, a radial sweep algorithm, a Theo Pavlidis'"'"' algorithm, an asymptotic decider algorithm, a cell-by-cell algorithm, and combinations thereof.

8. The method of claim 1, wherein the well logging tool comprises a logging-while-drilling tool.

9. The method of claim 1, wherein the well logging tool comprises a wireline tool, and wherein the azimuthally substantially continuous borehole image is obtained by processing azimuthally non-continuous borehole images using an in-painting algorithm.

10. A method for determining dip and azimuth of a subsurface formation, comprising:

moving a well logging tool along an interior of a wellbore penetrating the subsurface formation to generate a substantially azimuthally continuous image from measurements made by the well logging tool;

in a processor, extracting contours from the substantially azimuthally continuous image;

in the processor, grouping the extracted contours into clusters corresponding to a single transition zone;

in the processor, mapping the clustered contours having a measured depth extent that is greater than a length-scale over which a dip of the subsurface formation varies to a three-dimensional space corresponding to a coordinate system associated with the well; and

in the processor, calculating a dip and azimuth of a layer boundary of the subsurface formation from the mapped clustered contours and at least one of storing and displaying the calculated dip and azimuth.

11. The method of claim 10, wherein the contours comprise at least one of open contours, closes contours and combinations thereof.

12. The method of claim 10 further comprising, in the processor, prior to grouping the extracted contours into the clusters, filtering the extracted contours to remove spurious contours resulting from the extraction.

13. The method of claim 12, wherein the filtering comprises computing a Fourier series approximation for each extracted contour, and filtering out contours for which a quality of fit parameter is less than a selected threshold.

14. The method of claim 10 further comprising, in the processor, for extracted contours having a measured depth extent that is less than the length-scale over which the dip of the subterranean formation varies, estimating relative dip and apparent azimuth based on a first harmonic approximation of a contour.

15. The method of claim 10, wherein the well logging tool comprises at least one of a resistivity tool, a density tool and a photoelectric factor tool.

16. The method of claim 10, wherein the extracting the contours is performed using at least one of a marching-squares algorithm, a square tracing algorithm, a Moore-Neighbor algorithm, a radial sweep algorithm, a Theo Pavlidis'"'"' algorithm, an asymptotic decider algorithm, a cell-by-cell algorithm, and combinations thereof.

17. The method of claim 10, wherein the well logging tool comprises a logging-while-drilling tool.

18. The method of claim 17 wherein the moving the well logging tool comprises operating the well logging tool while drilling the wellbore

19. The method of claim 10, wherein the well logging tool comprises a wireline tool, and wherein the azimuthally substantially continuous borehole image is obtained by processing azimuthally non-continuous borehole images using an in-painting algorithm.

20. A system for determining dip and azimuth of a subsurface formation, comprising:

a processor programmed to accept as input thereto an azimuthally substantially continuous wellbore image obtained using a well logging tool moved along a wellbore penetrating a subsurface formation;

the processor programmed to extract contours from the substantially azimuthally continuous image;

the processor programmed to group the extracted contours into clusters corresponding to a single transition zone;

the processor programmed to map the clustered contours having a measured depth extent that is greater than a length-scale over which a dip of the subsurface formation varies to a three-dimensional space corresponding to a coordinate system associated with the well; and

the processor programmed to calculate a dip and azimuth of a layer boundary of the subsurface formation from the mapped clustered contours and at least one of storing and displaying the calculated dip and azimuth.