WO2016209265A1 - Identifying formation layer boundaries on well log measurements - Google Patents

Identifying formation layer boundaries on well log measurements Download PDF

Info

Publication number
WO2016209265A1
WO2016209265A1 PCT/US2015/037992 US2015037992W WO2016209265A1 WO 2016209265 A1 WO2016209265 A1 WO 2016209265A1 US 2015037992 W US2015037992 W US 2015037992W WO 2016209265 A1 WO2016209265 A1 WO 2016209265A1
Authority
WO
WIPO (PCT)
Prior art keywords
well log
log measurements
derivative
smoothing
formation layer
Prior art date
Application number
PCT/US2015/037992
Other languages
French (fr)
Inventor
Mehdi Eftekhari FAR
John Andrew Quirein
Natasa MEKIC
Original Assignee
Landmark Graphics Corporation
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Landmark Graphics Corporation filed Critical Landmark Graphics Corporation
Priority to PCT/US2015/037992 priority Critical patent/WO2016209265A1/en
Priority to US15/026,564 priority patent/US9874094B2/en
Priority to FR1654724A priority patent/FR3037992B1/en
Publication of WO2016209265A1 publication Critical patent/WO2016209265A1/en

Links

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V5/00—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
    • G01V5/04—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging

Definitions

  • the present disclosure generally relates to systems and methods for identifying formation layer boundaries on well log measurements. More particularly, the present disclosure relates to systems and methods for identifying formation layer boundaries on well log measurements using a second derivative of the last iteratively smoothed well log measurements and a fluctuation index to determine the extent of smoothing.
  • FIG. 1 is a flow diagram illustrating one embodiment of a method for implementing the present disclosure.
  • FIG. 2 is a well log illustrating synthetic (simulated) resistivity well log measurements used as exemplary input for step 102 in FIG. 1.
  • FIG. 3 is a first derivative well log illustrating a first derivative of the resistivity well log measurements in FIG. 2 and a first derivative of the last smoothed resistivity well log measurements in FIG. 2, which may be used in step 104 of FIG. 1.
  • FIG. 4 is a second derivative well log illustrating a second derivative of the resistivity well log measurements in FIG. 2 and a second derivative of the last smoothed resistivity well log measurements in FIG. 2 at each sampled point, which may be used in step 114 of FIG. 1.
  • FIGS. 5A-5B are two well logs illustrating real resistivity well log measurements at a predetermined depth in a well bore and a comparison of the formation layer boundaries displayed in step 118 of FIG. 1 (FIG. 5A) with the formation layer boundaries displayed according to the well-known variance based (FIG. 5B).
  • FIGS. 6A-6B are two well logs illustrating real resistivity well log measurements at another predetermined depth in the same well bore and a comparison of the formation layer boundaries displayed in step 118 of FIG. 1 (FIG. 6 A) with the formation layer boundaries displayed according to the well-known variance based (FIG. 6B).
  • FIGS. 7A-7B are two well logs illustrating real resistivity well log measurements at still another predetermined depth in the same well bore and a comparison of the formation layer boundaries displayed in step 118 of FIG. 1 (FIG. 7A) with the formation layer boundaries displayed according to the well-known variance based (FIG. 7B).
  • FIGS. 8A-8D are well logs illustrating real gamma ray well log measurements at a predetermined depth in a well bore and a comparison of the formation layer boundaries displayed in step 118 of FIG. 1 at different resolutions.
  • FIGS. 9A-9D are well logs illustrating real gamma ray well log measurements at a predetermined depth in a well bore and a comparison of the formation layer boundaries displayed in step 118 of FIG. 1 at different resolutions.
  • FIG. 10 is a block diagram illustrating one embodiment of a computer system for implementing the present disclosure.
  • the present disclosure overcomes one or more deficiencies in the prior art by providing systems and methods for identifying formation layer boundaries on well log measurements using a second derivative of the last iteratively smoothed well log measurements and a fluctuation index to determine the extent of smoothing.
  • the present disclosure includes a method for identifying formation layer boundaries on well log measurements, which comprises: a) computing a first derivative using one of a first derivative of the well log measurements and last smoothed well log measurements; b) computing a fluctuation index by dividing the first derivative by the number of data points in the well log measurements; c) smoothing one of the well log measurements and the last smoothed well log measurements using a computer processor; d) repeating steps a) - c) until the fluctuation index is less than a predetermined threshold; e) computing a second derivative using one of a second derivative of the well log measurements and the last smoothed well log measurements at each sampled point; f) identifying each second derivative with a zero value, which represents an inflection point; and g) displaying a formation layer boundary on one of the well log measurements and the last smoothed well log measurements at each identified inflection point.
  • the present disclosure includes a non-transitory program carrier device tangibly carrying computer executable instructions for identifying formation layer boundaries on well log measurements, the instructions being executable to implement: a) computing a first derivative using one of a first derivative of the well log measurements and last smoothed well log measurements; b) computing a fluctuation index by dividing the first derivative by the number of data points in the well log measurements; c) smoothing one of the well log measurements and the last smoothed well log measurements; d) repeating steps a) - c) until the fluctuation index is less than a predetermined threshold; e) computing a second derivative using one of a second derivative of the well log measurements and the last smoothed well log measurements at each sampled point; f) identifying each second derivative with a zero value, which represents an inflection point; and g) displaying a formation layer boundary on one of the well log measurements and the last smoothed well log measurements at each identified inflection point.
  • the present disclosure includes a non-transitory program carrier device tangibly carrying computer executable instructions for identifying formation layer boundaries on well log measurements, the instructions being executable to implement: a) computing a first derivative using one of a first derivative of the well log measurements and last smoothed well log measurements; b) computing a fluctuation index by dividing the first derivative by the number of data points in the well log measurements; c) smoothing one of the well log measurements and the last smoothed well log measurements; d) repeating steps a) - c) until the fluctuation index is less than a predetermined threshold, wherein the predetermined threshold represents an optimal formation layer boundary resolution; e) computing a second derivative using one of a second derivative of the well log measurements and the last smoothed well log measurements at each sampled point; and f) identifying each second derivative with a zero value, which represents an inflection point.
  • FIG. 1 a flow diagram illustrates one embodiment of a method 100 for implementing the present disclosure.
  • the method 100 enables the display of a layer boundary on the last smoothed well log measurements where each respective second derivative value is equal to zero.
  • the method 100 is based on a guided up-scaling of the data to be used for formation layer boundary detection.
  • the up-scaling is adjusted automatically according to data resolution.
  • step 102 well log measurements from a well bore and a predetermined formation layer boundary resolution are input using the client interface and/or the video interface described further in reference to FIG. 10.
  • the well log measurements may comprise, for example, resistivity or gamma ray measurements.
  • the predetermined formation layer boundary resolution is a parameter representing the thiclcness (resolution) of the formation layer boundaries displayed in step 118.
  • a well log illustrates synthetic resistivity well log measurements that may be used as exemplary input.
  • a first derivative (N) is computed using a first derivative of the well log measurements from step 102 or the last smoothed well log measurements from step 112, and techniques well known in the art.
  • the first derivative (N) represents the number of peaks and troughs, which are also referred to as the number of zero crossings, in the first derivative of the well log measurements from step 102 or the last smoothed well log measurements from step 112.
  • a first derivative well log illustrates a first derivative of the resistivity well log measurements in FIG. 2 and a first derivative of the last smoothed resistivity well log measurements in FIG. 2 as a result of step 112.
  • a fluctuation index (FI) is computed by dividing the first derivative (N) from step 104 by the number of data points (M) in the well log measurements from step 102.
  • step 108 the method 100 determines if FI from step 106 is less than a predetermined threshold representing the optimal thickness (resolution) of the formation layer boundaries displayed in step 118. If FI is less than a predetermined threshold, then the method 100 proceeds to step 114. Otherwise, the method 100 proceeds to step 110.
  • a smoothing window size and a number of smoothing iterations are defined using FI from step 106 a predetermined scaling factor (SF) and a predetermined constant (C).
  • the smoothing window size is equal to FI * SF and the number of smoothing iterations is equal to FI * C. The larger the smoothing window, the better (more smooth) the results.
  • step 112 smoothing is performed on the well log measurements from step 102 or the last smoothed well log measurements from step 112 using the smoothing window size and number of smoothing iterations from step 110, and smoothing techniques well-known in the art.
  • the method returns to step 104 and repeats steps 104-112 until FI is less than the predetermined threshold.
  • a second derivative ( ⁇ ') is computed using a second derivative of the well log measurements from step 102 or the last smoothed well log measurements from step 112 at each sampled point, and techniques well-known in the art.
  • the second derivative ( ⁇ ') thus, represents a value at each sampled point in the second derivative of the well log measurements from step 102 or the last smoothed well log measurements from step 112. If the second derivative ( ⁇ ') value is zero, then an inflection point is represented by the corresponding sampled point.
  • a second derivative well log illustrates a second derivative of the resistivity well log measurements in FIG. 2 and a second derivative of the last smoothed resistivity well log measurements in FIG. 2 (as a result of step 112) at each sampled point.
  • each second derivative ( ⁇ ') from step 114 with a zero value i.e. an inflection point
  • ⁇ ' a zero value
  • step 118 a formation layer boundary is displayed on the well log measurements from step 102 or the last smoothed well log measurements form step 112 at each inflection point identified in step 116 using the video interface described further in reference to FIG. 10.
  • each solid horizontal line represents a formation layer boundary displayed according to the method 100.
  • FIGS. 5A-5B two well logs illustrate real resistivity well log measurements (dashed lines) at a predetermined depth in a well bore and a comparison of the formation layer boundaries (solid horizontal lines) displayed in step 118 of FIG. 1 (FIG. 5A) with the formation layer boundaries (solid horizontal lines) displayed according to the well-known variance based (FIG. 5B).
  • FIGS. 6A-6B two well logs illustrate real resistivity well log measurements (dashed lines) at another predetermined depth in the same well bore and a comparison of the formation layer boundaries (solid horizontal lines) displayed in step 118 of FIG. 1 (FIG. 6A) with the formation layer boundaries (solid horizontal lines) displayed according to the well- known variance based (FIG. 6B).
  • FIGS. 7A-7B two well logs illustrate real resistivity well log measurements (dashed lines) at still another predetermined depth in the same well bore and a comparison of the formation layer boundaries (solid horizontal lines) displayed in step 118 of FIG. 1 (FIG. 7A) with the formation layer boundaries (solid horizontal lines) displayed according to the well-known variance based (FIG. 7B).
  • the method 100 is capable of identifying more actual formation layer boundaries on well log measurements than the popular variance based method.
  • the well logs illustrate real gamma ray well log measurements (dashed lines) at a predetermined depth in a well bore and a comparison of the formation layer boundaries (solid horizontal lines) displayed in step 118 of FIG. 1 at different resolutions.
  • the well logs illustrate real gamma ray well log measurements (dashed lines) at a predetermined depth in a well bore and a comparison of the formation layer boundaries (solid horizontal lines) displayed in step 118 of FIG. 1 at different resolutions.
  • the predetermined formation layer boundary resolution parameter in step 102 is capable of adjusting the thickness (resolution) of the formation layer boundaries displayed in step 118.
  • the method 100 is thus, useful for interpretation and processing of different petrophysical and geophysical data that may lead to more efficient drilling operations through the adjustment of the same based on the location of the formation layer boundaries.
  • the present disclosure may be implemented through a computer-executable program of instructions, such as program modules, generally referred to as software applications or application programs executed by a computer.
  • the software may include, for example, routines, programs, objects, components and data structures that perform particular tasks or implement particular abstract data types.
  • the software forms an interface to allow a computer to react according to a source of input.
  • DecisionSpace ® which is a commercial software application marketed by Landmark Graphics Corporation, may be used as an interface application to implement the present disclosure.
  • the software may also cooperate with other code segments to initiate a variety of tasks in response to data received in conjunction with the source of the received data.
  • the software may be stored and/or carried on any variety of memory such as CD- ROM, magnetic disk, bubble memory and semiconductor memory (e.g. various types of RAM or ROM). Furthermore, the software and its results may be transmitted over a variety of carrier media such as optical fiber, metallic wire and/or through any of a variety of networks, such as the Internet.
  • memory such as CD- ROM, magnetic disk, bubble memory and semiconductor memory (e.g. various types of RAM or ROM).
  • the software and its results may be transmitted over a variety of carrier media such as optical fiber, metallic wire and/or through any of a variety of networks, such as the Internet.
  • the disclosure may be practiced with a variety of computer-system configurations, including hand-held devices, multiprocessor systems, microprocessor-based or programmable-consumer electronics, minicomputers, mainframe computers, and the like. Any number of computer-systems and computer networks are acceptable for use with the present disclosure.
  • the disclosure may be practiced in distributed-computing environments where tasks are performed by remote- processing devices that are linked through a communications network.
  • program modules may be located in both local and remote computer- storage media including memory storage devices.
  • the present disclosure may therefore, be implemented in connection with various hardware, software or a combination thereof, in a computer system or other processing system.
  • FIG. 10 a block diagram illustrates one embodiment of a system for implementing the present disclosure on a computer.
  • the system includes a computing unit, sometimes referred to as a computing system, which contains memory, application programs, a client interface, a video interface, and a processing unit.
  • the computing unit is only one example of a suitable computing environment and is not intended to suggest any limitation as to the scope of use or functionality of the disclosure.
  • the memory primarily stores the application programs, which may also be described as program modules containing computer-executable instructions, executed by the computing unit for implementing the present disclosure described herein and illustrated in FIGS. 1-9.
  • the memory therefore, includes a formation layer boundary identification module, which enables steps 104-118 described in reference to FIG. 1.
  • the formation layer boundary identification module may integrate functionality from the remaining application programs illustrated in FIG. 10.
  • DecisionSpace ⁇ may be used as an interface application to perform step 102 in FIG. 1.
  • DecisionSpace ® may be used as interface application, other interface applications may be used, instead, or the formation layer boundary identification module may be used as a stand-alone application.
  • the computing unit typically includes a variety of computer readable media.
  • computer readable media may comprise computer storage media and communication media.
  • the computing system memory may include computer storage media in the form of volatile and/or nonvolatile memory such as a read only memory (ROM) and random access memory (RAM).
  • ROM read only memory
  • RAM random access memory
  • a basic input/output system (BIOS) containing the basic routines that help to transfer information between elements within the computing unit, such as during start-up, is typically stored in ROM.
  • the RAM typically contains data and/or program modules that are immediately accessible to, and/or presently being operated on, the processing unit.
  • the computing unit includes an operating system, application programs, other program modules, and program data.
  • the components shown in the memoiy may also be included in other removable/nonremovable, volatile/nonvolatile computer storage media or they may be implemented in the computing unit through an application program interface ("API") or cloud computing, which may reside on a separate computing unit connected through a computer system or network,
  • API application program interface
  • cloud computing which may reside on a separate computing unit connected through a computer system or network
  • a hard disk drive may read from or write to nonremovable, nonvolatile magnetic media
  • a magnetic disk drive may read from or write to a removable, nonvolatile magnetic disk
  • an optical disk drive may read from or write to a removable, nonvolatile optical disk such as a CD ROM or other optical media.
  • removable/nonremovable, volatile/nonvolatile computer storage media that can be used in the exemplary operating environment may include, but are not limited to, magnetic tape cassettes, flash memory cards, digital versatile disks, digital video tape, solid state RAM, solid state ROM, and the like.
  • the drives and their associated computer storage media discussed above provide storage of computer readable instructions, data structures, program modules and other data for the computing unit.
  • a client may enter commands and information into the computing unit through the client interface, which may be input devices such as a keyboard and pointing device, commonly referred to as a mouse, trackball or touch pad. Input devices may include a microphone, joystick, satellite dish, scanner, or the like. These and other input devices are often connected to the processing unit through the client interface that is coupled to a system bus, but may be connected by other interface and bus structures, such as a parallel port or a universal serial bus (USB).
  • USB universal serial bus
  • a monitor or other type of display device may be connected to the system bus via an interface, such as a video interface.
  • a graphical user interface may also be used with the video interface to receive instructions from the client interface and transmit instructions to the processing unit.
  • computers may also include other peripheral output devices such as speakers and printer, which may be connected through an output peripheral interface.

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • Geophysics (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
  • Measurement Of Radiation (AREA)

Abstract

Systems and methods for identifying formation layer boundaries on well log measurements using a second derivative of the last iteratively smoothed well log measurements and a fluctuation index to determine the extent of smoothing.

Description

IDENTIFYING FORM ATION LAYER BOUNDARIES
ON WELL LOG MEASUREMENTS
FIELD OF THE DISCLOSURE
[0001] The present disclosure generally relates to systems and methods for identifying formation layer boundaries on well log measurements. More particularly, the present disclosure relates to systems and methods for identifying formation layer boundaries on well log measurements using a second derivative of the last iteratively smoothed well log measurements and a fluctuation index to determine the extent of smoothing.
BACKGROUND
[0002] Detection of formation layer boundaries on well log measurements is required for interpretation, modeling and inversion tasks. Existing methods for automatic formation boundary detection, which are mostly based on data variance and frequency analysis, are not robust and effective. Conventional formation layer boundary detection on well log measurements therefore, is often manually performed and time consuming. These disadvantages are compounded in situations involving multi-sensor data and multi-well data.
BRIEF DESCRIPTION OF THE DRAWINGS
[0003] The present disclosure is described below with references to the accompanying drawings in which like elements are referenced with like reference numerals, and in which:
[0004] FIG. 1 is a flow diagram illustrating one embodiment of a method for implementing the present disclosure.
[0005] FIG. 2 is a well log illustrating synthetic (simulated) resistivity well log measurements used as exemplary input for step 102 in FIG. 1.
[0006] FIG. 3 is a first derivative well log illustrating a first derivative of the resistivity well log measurements in FIG. 2 and a first derivative of the last smoothed resistivity well log measurements in FIG. 2, which may be used in step 104 of FIG. 1.
[0007] FIG. 4 is a second derivative well log illustrating a second derivative of the resistivity well log measurements in FIG. 2 and a second derivative of the last smoothed resistivity well log measurements in FIG. 2 at each sampled point, which may be used in step 114 of FIG. 1.
[0008] FIGS. 5A-5B are two well logs illustrating real resistivity well log measurements at a predetermined depth in a well bore and a comparison of the formation layer boundaries displayed in step 118 of FIG. 1 (FIG. 5A) with the formation layer boundaries displayed according to the well-known variance based (FIG. 5B).
[0009] FIGS. 6A-6B are two well logs illustrating real resistivity well log measurements at another predetermined depth in the same well bore and a comparison of the formation layer boundaries displayed in step 118 of FIG. 1 (FIG. 6 A) with the formation layer boundaries displayed according to the well-known variance based (FIG. 6B).
[0010] FIGS. 7A-7B are two well logs illustrating real resistivity well log measurements at still another predetermined depth in the same well bore and a comparison of the formation layer boundaries displayed in step 118 of FIG. 1 (FIG. 7A) with the formation layer boundaries displayed according to the well-known variance based (FIG. 7B).
[0011] FIGS. 8A-8D are well logs illustrating real gamma ray well log measurements at a predetermined depth in a well bore and a comparison of the formation layer boundaries displayed in step 118 of FIG. 1 at different resolutions.
[0012] FIGS. 9A-9D are well logs illustrating real gamma ray well log measurements at a predetermined depth in a well bore and a comparison of the formation layer boundaries displayed in step 118 of FIG. 1 at different resolutions.
[0013] FIG. 10 is a block diagram illustrating one embodiment of a computer system for implementing the present disclosure.
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
[0014] The present disclosure overcomes one or more deficiencies in the prior art by providing systems and methods for identifying formation layer boundaries on well log measurements using a second derivative of the last iteratively smoothed well log measurements and a fluctuation index to determine the extent of smoothing.
[0015] In one embodiment, the present disclosure includes a method for identifying formation layer boundaries on well log measurements, which comprises: a) computing a first derivative using one of a first derivative of the well log measurements and last smoothed well log measurements; b) computing a fluctuation index by dividing the first derivative by the number of data points in the well log measurements; c) smoothing one of the well log measurements and the last smoothed well log measurements using a computer processor; d) repeating steps a) - c) until the fluctuation index is less than a predetermined threshold; e) computing a second derivative using one of a second derivative of the well log measurements and the last smoothed well log measurements at each sampled point; f) identifying each second derivative with a zero value, which represents an inflection point; and g) displaying a formation layer boundary on one of the well log measurements and the last smoothed well log measurements at each identified inflection point.
[0016] In another embodiment, the present disclosure includes a non-transitory program carrier device tangibly carrying computer executable instructions for identifying formation layer boundaries on well log measurements, the instructions being executable to implement: a) computing a first derivative using one of a first derivative of the well log measurements and last smoothed well log measurements; b) computing a fluctuation index by dividing the first derivative by the number of data points in the well log measurements; c) smoothing one of the well log measurements and the last smoothed well log measurements; d) repeating steps a) - c) until the fluctuation index is less than a predetermined threshold; e) computing a second derivative using one of a second derivative of the well log measurements and the last smoothed well log measurements at each sampled point; f) identifying each second derivative with a zero value, which represents an inflection point; and g) displaying a formation layer boundary on one of the well log measurements and the last smoothed well log measurements at each identified inflection point.
[0017] In yet another embodiment, the present disclosure includes a non-transitory program carrier device tangibly carrying computer executable instructions for identifying formation layer boundaries on well log measurements, the instructions being executable to implement: a) computing a first derivative using one of a first derivative of the well log measurements and last smoothed well log measurements; b) computing a fluctuation index by dividing the first derivative by the number of data points in the well log measurements; c) smoothing one of the well log measurements and the last smoothed well log measurements; d) repeating steps a) - c) until the fluctuation index is less than a predetermined threshold, wherein the predetermined threshold represents an optimal formation layer boundary resolution; e) computing a second derivative using one of a second derivative of the well log measurements and the last smoothed well log measurements at each sampled point; and f) identifying each second derivative with a zero value, which represents an inflection point.
[0018] The subject matter of the present disclosure is described with specificity, however, the description itself is not intended to limit the scope of the disclosure. The subject matter thus, might also be embodied in other ways, to include different structures, steps and/or combinations similar to those described herein, in conjunction with other present or future technologies. Moreover, although the term "step" may be used herein to describe different elements of methods employed, the term should not be interpreted as implying any particular order among or between various steps herein disclosed unless otherwise expressly limited by the description to a particular order. While the present disclosure may be applied in the oil and gas industry, it is not limited thereto and may also be applied in other industries (e.g. drilling water wells) to achieve similar results.
Method Description
[0019] Referring now to FIG. 1, a flow diagram illustrates one embodiment of a method 100 for implementing the present disclosure. The method 100 enables the display of a layer boundary on the last smoothed well log measurements where each respective second derivative value is equal to zero. In this manner, the method 100 is based on a guided up-scaling of the data to be used for formation layer boundary detection. In order for the method 100 to be automatic, independent of the data resolution and input parameters, the up-scaling is adjusted automatically according to data resolution.
[0020] In step 102, well log measurements from a well bore and a predetermined formation layer boundary resolution are input using the client interface and/or the video interface described further in reference to FIG. 10. The well log measurements may comprise, for example, resistivity or gamma ray measurements. And, the predetermined formation layer boundary resolution is a parameter representing the thiclcness (resolution) of the formation layer boundaries displayed in step 118. In FIG. 2, a well log illustrates synthetic resistivity well log measurements that may be used as exemplary input.
[0021] In step 104, a first derivative (N) is computed using a first derivative of the well log measurements from step 102 or the last smoothed well log measurements from step 112, and techniques well known in the art. The first derivative (N) represents the number of peaks and troughs, which are also referred to as the number of zero crossings, in the first derivative of the well log measurements from step 102 or the last smoothed well log measurements from step 112. In FIG. 3, a first derivative well log illustrates a first derivative of the resistivity well log measurements in FIG. 2 and a first derivative of the last smoothed resistivity well log measurements in FIG. 2 as a result of step 112.
[0022] In step 106, a fluctuation index (FI) is computed by dividing the first derivative (N) from step 104 by the number of data points (M) in the well log measurements from step 102.
[0023] In step 108, the method 100 determines if FI from step 106 is less than a predetermined threshold representing the optimal thickness (resolution) of the formation layer boundaries displayed in step 118. If FI is less than a predetermined threshold, then the method 100 proceeds to step 114. Otherwise, the method 100 proceeds to step 110.
[0024] In step 110, a smoothing window size and a number of smoothing iterations are defined using FI from step 106 a predetermined scaling factor (SF) and a predetermined constant (C). The smoothing window size is equal to FI * SF and the number of smoothing iterations is equal to FI * C. The larger the smoothing window, the better (more smooth) the results.
[0025] In step 112, smoothing is performed on the well log measurements from step 102 or the last smoothed well log measurements from step 112 using the smoothing window size and number of smoothing iterations from step 110, and smoothing techniques well-known in the art. The method returns to step 104 and repeats steps 104-112 until FI is less than the predetermined threshold.
[0026] In step 114, a second derivative (Ν') is computed using a second derivative of the well log measurements from step 102 or the last smoothed well log measurements from step 112 at each sampled point, and techniques well-known in the art. The second derivative (Ν') thus, represents a value at each sampled point in the second derivative of the well log measurements from step 102 or the last smoothed well log measurements from step 112. If the second derivative (Ν') value is zero, then an inflection point is represented by the corresponding sampled point. In FIG. 4, a second derivative well log illustrates a second derivative of the resistivity well log measurements in FIG. 2 and a second derivative of the last smoothed resistivity well log measurements in FIG. 2 (as a result of step 112) at each sampled point.
[0027] In step 116, each second derivative (Ν') from step 114 with a zero value (i.e. an inflection point) is identified. Each inflection point represents a formation layer boundary.
[0028] In step 118, a formation layer boundary is displayed on the well log measurements from step 102 or the last smoothed well log measurements form step 112 at each inflection point identified in step 116 using the video interface described further in reference to FIG. 10. In FIG. 4, each solid horizontal line represents a formation layer boundary displayed according to the method 100. If the second derivative is computed in step 114 using the well log measurements from step 102, then many unwanted thin formation layer boundaries may be detected compared to if the second derivative is computed using the last smoothed well log measurements from step 112, which is more favorable for interpretation and modeling. The smoothing performed in step 112 thus, can remove small fluctuations from the data resulting in a reasonable number of inflection points identified in step 116. A good measure of the fluctuations is the number of peaks and troughs in the well log measurements from step 102, normalized by the number of data points.
[0029] In FIGS. 5A-5B, two well logs illustrate real resistivity well log measurements (dashed lines) at a predetermined depth in a well bore and a comparison of the formation layer boundaries (solid horizontal lines) displayed in step 118 of FIG. 1 (FIG. 5A) with the formation layer boundaries (solid horizontal lines) displayed according to the well-known variance based (FIG. 5B). In FIGS. 6A-6B, two well logs illustrate real resistivity well log measurements (dashed lines) at another predetermined depth in the same well bore and a comparison of the formation layer boundaries (solid horizontal lines) displayed in step 118 of FIG. 1 (FIG. 6A) with the formation layer boundaries (solid horizontal lines) displayed according to the well- known variance based (FIG. 6B). In FIGS. 7A-7B, two well logs illustrate real resistivity well log measurements (dashed lines) at still another predetermined depth in the same well bore and a comparison of the formation layer boundaries (solid horizontal lines) displayed in step 118 of FIG. 1 (FIG. 7A) with the formation layer boundaries (solid horizontal lines) displayed according to the well-known variance based (FIG. 7B). As demonstrated by these figures, the method 100 is capable of identifying more actual formation layer boundaries on well log measurements than the popular variance based method.
[0030] In FIGS. 8A-8D, the well logs illustrate real gamma ray well log measurements (dashed lines) at a predetermined depth in a well bore and a comparison of the formation layer boundaries (solid horizontal lines) displayed in step 118 of FIG. 1 at different resolutions. In FIGS. 9A-9D, the well logs illustrate real gamma ray well log measurements (dashed lines) at a predetermined depth in a well bore and a comparison of the formation layer boundaries (solid horizontal lines) displayed in step 118 of FIG. 1 at different resolutions. As demonstrated by these figures, the predetermined formation layer boundary resolution parameter in step 102 is capable of adjusting the thickness (resolution) of the formation layer boundaries displayed in step 118. The method 100 is thus, useful for interpretation and processing of different petrophysical and geophysical data that may lead to more efficient drilling operations through the adjustment of the same based on the location of the formation layer boundaries.
System Description
[0031] The present disclosure may be implemented through a computer-executable program of instructions, such as program modules, generally referred to as software applications or application programs executed by a computer. The software may include, for example, routines, programs, objects, components and data structures that perform particular tasks or implement particular abstract data types. The software forms an interface to allow a computer to react according to a source of input. DecisionSpace®, which is a commercial software application marketed by Landmark Graphics Corporation, may be used as an interface application to implement the present disclosure. The software may also cooperate with other code segments to initiate a variety of tasks in response to data received in conjunction with the source of the received data. The software may be stored and/or carried on any variety of memory such as CD- ROM, magnetic disk, bubble memory and semiconductor memory (e.g. various types of RAM or ROM). Furthermore, the software and its results may be transmitted over a variety of carrier media such as optical fiber, metallic wire and/or through any of a variety of networks, such as the Internet.
[0032] Moreover, those skilled in the art will appreciate that the disclosure may be practiced with a variety of computer-system configurations, including hand-held devices, multiprocessor systems, microprocessor-based or programmable-consumer electronics, minicomputers, mainframe computers, and the like. Any number of computer-systems and computer networks are acceptable for use with the present disclosure. The disclosure may be practiced in distributed-computing environments where tasks are performed by remote- processing devices that are linked through a communications network. In a distributed- computing environment, program modules may be located in both local and remote computer- storage media including memory storage devices. The present disclosure may therefore, be implemented in connection with various hardware, software or a combination thereof, in a computer system or other processing system.
[0033] Referring now to FIG. 10, a block diagram illustrates one embodiment of a system for implementing the present disclosure on a computer. The system includes a computing unit, sometimes referred to as a computing system, which contains memory, application programs, a client interface, a video interface, and a processing unit. The computing unit is only one example of a suitable computing environment and is not intended to suggest any limitation as to the scope of use or functionality of the disclosure.
[0034] The memory primarily stores the application programs, which may also be described as program modules containing computer-executable instructions, executed by the computing unit for implementing the present disclosure described herein and illustrated in FIGS. 1-9. The memory therefore, includes a formation layer boundary identification module, which enables steps 104-118 described in reference to FIG. 1. The formation layer boundary identification module may integrate functionality from the remaining application programs illustrated in FIG. 10. In particular, DecisionSpace© may be used as an interface application to perform step 102 in FIG. 1. Although DecisionSpace® may be used as interface application, other interface applications may be used, instead, or the formation layer boundary identification module may be used as a stand-alone application.
[0035] Although the computing unit is shown as having a generalized memory, the computing unit typically includes a variety of computer readable media. By way of example, and not limitation, computer readable media may comprise computer storage media and communication media. The computing system memory may include computer storage media in the form of volatile and/or nonvolatile memory such as a read only memory (ROM) and random access memory (RAM). A basic input/output system (BIOS), containing the basic routines that help to transfer information between elements within the computing unit, such as during start-up, is typically stored in ROM. The RAM typically contains data and/or program modules that are immediately accessible to, and/or presently being operated on, the processing unit. By way of example, and not limitation, the computing unit includes an operating system, application programs, other program modules, and program data.
[0036] The components shown in the memoiy may also be included in other removable/nonremovable, volatile/nonvolatile computer storage media or they may be implemented in the computing unit through an application program interface ("API") or cloud computing, which may reside on a separate computing unit connected through a computer system or network, For example only, a hard disk drive may read from or write to nonremovable, nonvolatile magnetic media, a magnetic disk drive may read from or write to a removable, nonvolatile magnetic disk, and an optical disk drive may read from or write to a removable, nonvolatile optical disk such as a CD ROM or other optical media. Other removable/nonremovable, volatile/nonvolatile computer storage media that can be used in the exemplary operating environment may include, but are not limited to, magnetic tape cassettes, flash memory cards, digital versatile disks, digital video tape, solid state RAM, solid state ROM, and the like. The drives and their associated computer storage media discussed above provide storage of computer readable instructions, data structures, program modules and other data for the computing unit.
[0037] A client may enter commands and information into the computing unit through the client interface, which may be input devices such as a keyboard and pointing device, commonly referred to as a mouse, trackball or touch pad. Input devices may include a microphone, joystick, satellite dish, scanner, or the like. These and other input devices are often connected to the processing unit through the client interface that is coupled to a system bus, but may be connected by other interface and bus structures, such as a parallel port or a universal serial bus (USB).
[0038] A monitor or other type of display device may be connected to the system bus via an interface, such as a video interface. A graphical user interface ("GUI") may also be used with the video interface to receive instructions from the client interface and transmit instructions to the processing unit. In addition to the monitor, computers may also include other peripheral output devices such as speakers and printer, which may be connected through an output peripheral interface.
[0039] Although many other internal components of the computing unit are not shown, those of ordinary skill in the art will appreciate that such components and their interconnection are well-known.
[0040] While the present disclosure has been described in connection with presently preferred embodiments, it will be understood by those skilled in the art that it is not intended to limit the disclosure to those embodiments. It is therefore, contemplated that various alternative embodiments and modifications may be made to the disclosed embodiments without departing from the spirit and scope of the disclosure defined by the appended claims and equivalents thereof.

Claims

1. A method for identifying formation layer boundaries on well log measurements, which comprises: a) computing a first derivative using one of a first derivative of the well log measurements and last smoothed well log measurements; b) computing a fluctuation index by dividing the first derivative by the number of data points in the well log measurements; c) smoothing one of the well log measurements and the last smoothed well log measurements using a computer processor; d) repeating steps a) - c) until the fluctuation index is less than a predetermined threshold; e) computing a second derivative using one of a second derivative of the well log measurements and the last smoothed well log measurements at each sampled point; f) identifying each second derivative with a zero value, which represents an inflection point; and g) displaying a formation layer boundary on one of the well log measurements and the last smoothed well log measurements at each identified inflection point.
2. The method of claim 1, further comprising adjusting one or more drilling operations based on each displayed formation layer boundary.
3. The method of claim 1, wherein the well log measurements comprise one of resistivity and gamma ray measurements.
4. The method of claim 1, wherein the predetermined threshold represents an optimal formation layer boundary resolution.
5. The method of claim 1, wherein the smoothing is based on a smoothing window size and a number of smoothing iterations.
6. The method of claim 5, wherein the smoothing window size is determined by multiplying the fluctuation index by a predetermined scaling factor.
7. The method of claim 5, wherein the number of smoothing iterations is determined by multiplying the fluctuation index by a predetermined constant.
8. The method, of claim 1, further comprising adjusting the predetermined threshold.
9. A non-transitory program carrier device tangibly carrying computer executable instructions for identifying formation layer boundaries on well log measurements, the instructions being executable to implement: a) computing a first derivative using one of a first derivative of the well log measurements and last smoothed well log measurements; b) computing a fluctuation index by dividing the first derivative by the number of data points in the well log measurements; c) smoothing one of the well log measurements and the last smoothed well log measurements; d) repeating steps a) - c) until the fluctuation index is less than a predetermined threshold; e) computing a second derivative using one of a second derivative of the well log measurements and the last smoothed well log measurements at each sampled point; f) identifying each second derivative with a zero value, which represents an inflection point; and g) displaying a formation layer boundary on one of the well log measurements and the last smoothed well log measurements at each identified inflection point.
10. The program carrier device of claim 9, further comprising adjusting one or more drilling operations based on each displayed formation layer boundary.
1 1. The program carrier device of claim 9, wherein the well log measurements comprise one of resistivity and gamma ray measurements.
12. The program carrier device of claim 9, wherein the predetermined threshold represents an optimal formation layer boundary resolution.
13. The program carrier device of claim 9, wherein the smoothing is based on a smoothing window size and a number of smoothing iterations.
14. The program carrier device of claim 13, wherein the smoothing window size is determined by multiplying the fluctuation index by a predetermined scaling factor.
15. The program carrier device of claim 13, wherein the number of smoothing iterations is determined by multiplying the fluctuation index by a predetermined constant.
16. The program carrier device of claim 9 further comprising adjusting the predetermined threshold.
17. A non-transitory program carrier device tangibly carrying computer executable instructions for identifying formation layer boundaries on well log measurements, the instructions being executable to implement: a) computing a first derivative using one of a first derivative of the well log measurements and last smoothed well log measurements; b) computing a fluctuation index by dividing the first derivative by the number of data points in the well log measurements; c) smoothing one of the well log measurements and the last smoothed well log measurements; d) repeating steps a) - c) until the fluctuation index is less than a predetermined threshold, wherein the predetermined threshold represents an optimal formation layer boundary resolution; e) computing a second derivative using one of a second derivative of the well log measurements and the last smoothed well log measurements at each sampled point; and f) identifying each second derivative with a zero value, which represents an inflection point.
18. The program carrier device of claim 17, wherein the smoothing is based on a smoothing window size and a number of smoothing iterations.
19. The program carrier device of claim 18, wherein the smoothing window size is determined by multiplying the fluctuation index by a predetermined scaling factor.
20. The program carrier device of claim 18, wherein the number of smoothing iterations is determined by multiplying the fluctuation index by a predetermined constant.
PCT/US2015/037992 2014-07-25 2015-06-26 Identifying formation layer boundaries on well log measurements WO2016209265A1 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
PCT/US2015/037992 WO2016209265A1 (en) 2015-06-26 2015-06-26 Identifying formation layer boundaries on well log measurements
US15/026,564 US9874094B2 (en) 2014-07-25 2015-06-26 Identifying formation layer boundaries on well log measurements
FR1654724A FR3037992B1 (en) 2015-06-26 2016-05-26 IDENTIFICATION OF TRAINING LAYER LIMITS ON WELL DIAGRAM MEASUREMENTS

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/US2015/037992 WO2016209265A1 (en) 2015-06-26 2015-06-26 Identifying formation layer boundaries on well log measurements

Publications (1)

Publication Number Publication Date
WO2016209265A1 true WO2016209265A1 (en) 2016-12-29

Family

ID=57538939

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2015/037992 WO2016209265A1 (en) 2014-07-25 2015-06-26 Identifying formation layer boundaries on well log measurements

Country Status (2)

Country Link
FR (1) FR3037992B1 (en)
WO (1) WO2016209265A1 (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4585939A (en) * 1983-10-05 1986-04-29 Halliburton Company Multi-function natural gamma ray logging system
US6466872B1 (en) * 1999-11-08 2002-10-15 Baker Hughes Incorporated Method for determination of apparent resistivities of anisotropic reservoirs
US6856909B2 (en) * 2001-03-20 2005-02-15 Shell Oil Company Detecting a boundary in a formation
US20060186887A1 (en) * 2005-02-22 2006-08-24 Strack Kurt M Method for identifying subsurface features from marine transient controlled source electromagnetic surveys
US20080015780A1 (en) * 2006-03-30 2008-01-17 Council Of Scientific And Industrial Research Non-linear inversion technique for interpretation of geophysical data using analytically computed first and second order derivatives

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4585939A (en) * 1983-10-05 1986-04-29 Halliburton Company Multi-function natural gamma ray logging system
US6466872B1 (en) * 1999-11-08 2002-10-15 Baker Hughes Incorporated Method for determination of apparent resistivities of anisotropic reservoirs
US6856909B2 (en) * 2001-03-20 2005-02-15 Shell Oil Company Detecting a boundary in a formation
US20060186887A1 (en) * 2005-02-22 2006-08-24 Strack Kurt M Method for identifying subsurface features from marine transient controlled source electromagnetic surveys
US20080015780A1 (en) * 2006-03-30 2008-01-17 Council Of Scientific And Industrial Research Non-linear inversion technique for interpretation of geophysical data using analytically computed first and second order derivatives

Also Published As

Publication number Publication date
FR3037992B1 (en) 2019-03-22
FR3037992A1 (en) 2016-12-30

Similar Documents

Publication Publication Date Title
US10385658B2 (en) In-situ wellbore, core and cuttings information system
EP2973429B1 (en) Basin-to-reservoir modeling
AU2013392090B2 (en) Systems and methods for optimizing existing wells and designing new wells based on the distribution of average effective fracture lengths
US10197692B2 (en) Validation of depth-depth curves using time-seismic depth interval velocity
FR3038338A1 (en) CORRECTION OF DEVIATION AND DISPERSION EFFECTS ON ACOUSTICAL DIAGRAM MEASUREMENTS OF WELLS DUE TO STRATIFIED FORMATIONS
US9874094B2 (en) Identifying formation layer boundaries on well log measurements
US11549355B2 (en) Avoiding geological formation boundaries during drilling operations
NO20180862A1 (en) Transferring logging data from an offset well location to a target well location
WO2016209265A1 (en) Identifying formation layer boundaries on well log measurements
US11085292B2 (en) Avoiding geological bodies that are hazards to drilling operations
CA2796915C (en) Systems and methods for computing a default 3d variogram model
US20180217286A1 (en) Selecting logging data for petrophysical modelling and completion optimization
US20190004197A1 (en) Updating models of complex geological sequences field of the disclosure

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 15026564

Country of ref document: US

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 15896551

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205 DATED 06.03.2018)

122 Ep: pct application non-entry in european phase

Ref document number: 15896551

Country of ref document: EP

Kind code of ref document: A1