WO2012036783A1 - System and method for sweet zone identification in shale gas reservoirs - Google Patents

System and method for sweet zone identification in shale gas reservoirs Download PDF

Info

Publication number
WO2012036783A1
WO2012036783A1 PCT/US2011/044132 US2011044132W WO2012036783A1 WO 2012036783 A1 WO2012036783 A1 WO 2012036783A1 US 2011044132 W US2011044132 W US 2011044132W WO 2012036783 A1 WO2012036783 A1 WO 2012036783A1
Authority
WO
WIPO (PCT)
Prior art keywords
data
neutron
density
radioactivity
porosity
Prior art date
Application number
PCT/US2011/044132
Other languages
French (fr)
Inventor
Chengbing Liu
Original Assignee
Chevron U.S.A. Inc.
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 Chevron U.S.A. Inc. filed Critical Chevron U.S.A. Inc.
Priority to CA2809969A priority Critical patent/CA2809969C/en
Priority to CN2011800439070A priority patent/CN103098062A/en
Priority to EP11825591.8A priority patent/EP2616978A1/en
Priority to AU2011302598A priority patent/AU2011302598B2/en
Priority to JP2013528197A priority patent/JP2013542412A/en
Priority to EA201390369A priority patent/EA201390369A1/en
Priority to BR112013005708A priority patent/BR112013005708A2/en
Publication of WO2012036783A1 publication Critical patent/WO2012036783A1/en

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing 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
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/10Locating fluid leaks, intrusions or movements

Definitions

  • the present invention relates generally to methods and systems for identification of the sweet zone in shale gas reservoirs and more particularly to combining types of well log information to identify the sweet zone.
  • a computer implemented method for automatically identifying a hydrocarbon (such as kerogen, gas, oil) rich zone in a well bore includes obtaining well log data including neutron data, density data, radioactivity data, and resistivity data representative of physical characteristics of a formation surrounding the well bore and computing an apparent neutron porosity and an apparent density porosity based on the neutron data and density data.
  • a normalized neutron-density separation is computed based on the computed apparent neutron porosity and the computed apparent density porosity and a baseline of normal shale is determined for each data type.
  • a computer system for automatically identifying a hydrocarbon rich zone in a well bore includes a computer readable medium having computer readable well log data stored thereon, the well log data including neutron data, density data, radioactivity data, and resistivity data representative of physical characteristics of a formation surrounding the well bore.
  • a processor of the computer system is configured and arranged to compute an apparent neutron porosity and an apparent density porosity based on the neutron data and density data, to compute a normalized neutron-density separation based on the computed apparent neutron porosity and the computed apparent density porosity, to compute a baseline of normal shale for the neutron data, density data, radioactivity data and resistivity data, and to compute the presence or absence of a hydrocarbon rich zone based on the computed normalized neutron-density separation, the radioactivity data, the resistivity data, and the determined baselines.
  • the computations outlined above are done at each depth level logged in the well.
  • Figure 1 is a flowchart illustrating a method in accordance with an embodiment of the invention
  • Figure 2 is an example of a set of well logs showing a determined sweet zone indicator and sweet zone quality index in accordance with an embodiment of the invention.
  • Figure 3 schematically illustrates a system for performing a method in accordance with an embodiment of the invention.
  • One method of characterizing a formation is to make measurements of characteristics along a borehole penetrating the formation, either during or after drilling operations, i.e., well logging.
  • Well logging includes a number of techniques including resistivity/conductivity measurements, ultrasound, NMR, neutron, density, uranium concentration and radiation scattering, for example.
  • Borehole data of this type is often used to replace or supplement the collection of cores for direct inspection.
  • logged borehole data is analyzed by human interpreters in order to characterize a subsurface geological formation to allow decisions to be made regarding potential of the well or to determine information about the nature of the surrounding geologic area.
  • the inventors have determined that by combining information from a variety of well logs, a quantitative approach may be pursued to identify formations or portions of formations that are likely to be rich in organic material and therefore likely to offer potential in hydrocarbon production, without requiring human interpretation.
  • well log data is obtained.
  • the well log data comprises neutron, density, uranium concentration and resistivity data.
  • uranium concentration is replaced by gamma ray data, each being a type of radioactivity data.
  • the well log data may be acquired by any of a variety of well logging techniques, or may be existing well log data stored locally or remotely from a computer system on which the method is executed. In a particular example and not by way of limitation, the well log data may be from a shale formation.
  • Equation 1 sets out the calculation for PHIT D:
  • Equation 1 P M is the density of the rock matrix (where the matrix is selected to be a calcite matrix or other appropriate matrix, depending on the geology of the shale formation), P B is the bulk density of the rock, and P F is the density of fluid in the rock (where the fluid may be selected to be water).
  • this Equation will produce a value of 0.0 where the ratio (PM - PB)/ (PM - PF) is negative, 1.0 when the ratio is greater than one, and the value of the ratio where the ratio is between zero and one. That is, it calculates a porosity value that is bounded by zero and one.
  • Equation 2 an apparent neutron porosity
  • Equation 2 TMPH is the neutron porosity reading of the rock, TNPM is the neutron porosity of the matrix and TNPF is the neutron porosity of the fluid.
  • this Equation produces a value equal to the ratio (TNPH - TNPM)/ (TNPF - TNPM) for values between zero and one, and is bounded by zero and one for all other values of the ratio.
  • VWSH NDS normalized neutron-density separation
  • VWSH NDS max(min([(PHIT_N - PHIT D) - (PHIT N - PHIT_D) min ]/ [(PHIT N - PHIT_D) ns - (PHIT N - PHIT_D) min ], 1.0), -1.0) Eqn. 3
  • Equation 3 the newly introduced quantity (PHIT N - PHIT_D) ns is the neutron-density separation for normal shales, while (PHIT N - PHIT_D) m i n represents a minimum value of the neutron-density separation.
  • (PHIT N - PHIT_D) m i n is taken to be zero and that portion of the numerator and denominator is eliminated. This equation produces values between minus one and one, although in most cases the values are between zero and one.
  • a baseline value for each of the quantities is determined. For an embodiment using neutron, density, uranium concentration and resistivity data, baselines are determined for each of these. For embodiments in which gamma ray data replaces uranium concentration data, a baseline for gamma ray log readings is determined.
  • step 20a the values determined in the preceding steps are used to generate a sweet zone indicator (RNR) in accordance with the if statement in Equation 4.
  • VWSH NDS NSBSL is the normalized neutron-density separation baseline for normal shales
  • URAN is a uranium concentration
  • URAN NSBSL is baseline uranium concentration for normal shales
  • RD is a resistivity value of the log data
  • RD NSBSL is a baseline resistivity for normal shales
  • FVBSL, FUBSL and FRBSL are adjustment factors for the respective baselines.
  • the baseline for each type of log may be a constant, or may vary with depth and thus be represented by a curve or trendline, depending on the geological or borehole conditions.
  • a shale interval is chosen to determine the baseline value or curve.
  • the respective adjustment factors, FVBSL, FUBSL and FRBSL are selected to reduce measurement noise and also to reduce high frequency variations in the actual geological structure, thereby improving reliability of the indicator. In an embodiment, these are determined by Monte Carlo experimentation.
  • the adjustment factors may also be adjusted in accordance with the experience of a user based on local geological conditions, analogues, and data quality and/or data provenance.
  • Equation 4 is replaced by Equation 5.
  • steps 12 and 14 could be performed in any order.
  • step 18 the baseline determination for each type of well log performed in step 18 could, in principle, be performed in advance of any of the other calculations, and after all calculations except those of step 20, which depend on the results of step 18.
  • Equation 4 or 5 Evaluation of either Equation 4 or 5 will return a value of one or zero, indicating presence or absence of a sweet zone respectively.
  • the indicator may then be used as a basis for determining a depth to initiate a horizontal drilling operation, or otherwise to guide production drilling decisions.
  • Figure 2 illustrates a number of well logs and derivative products in accordance with an embodiment of the invention.
  • the first column shows radiation data derived from gamma ray measurements.
  • the space between the two curves in the central portion of the log is indicative of uranium and represents a difference between spectral gamma radiation (the right hand curve) and computed gamma radiation (left hand curve).
  • Some additional curves along the left-hand side of the trace are not relevant to the method described herein.
  • the second column shows depth of the well.
  • the third column shows resistivity data for a number of different depths of investigation.
  • the fourth column shows neutron and density data and the fifth shows uranium data.
  • the indicator may be supplemented with a quality index that quantifies the quality of the identified sweet zone. This is illustrated in Figure 2, in the sixth column where 30 indicates a region in which the sweet zone indicator is one and 32 is a curve indicating the quality index within the zone 30.
  • the region 30 corresponds to the shaded region in column 5 where normalized neutron-density separation is less than its baseline and the intersection of that shaded region with the shaded region in column 6 where uranium concentration is above its respective baseline.
  • resistivity is above its baseline substantially throughout the region in which normalized neutron- density separation is less than its baseline.
  • a sweet zone quality index may be calculated based on the data used to determine the sweet zone indicator.
  • quality indexes are calculated for each of the data types, then those calculated quality indexes are used to compute an overall quality index that allows for comparison between or among various formations.
  • SQI NDS min(max([VWSH_NDS_NSBSL - VWSH_NDS]/[VWSH_NDS_NSBSL -
  • SQI_RD min(max([logio(RD)-logio(RD_NSBSL)]/[logio(RD max )-logio(RD_NSBSL)],0), 1)
  • SQI min(max([SQI_NDS ⁇ (W nds /(W nds + W gr + W rd )) + SQI_GR ⁇ (W gr /(W nds + Wgr +
  • Equation 10 and 11 the choice between Equation 10 and 11 will depend on availability of uranium data. Where uranium data is not available, gamma ray data is used in accordance with Equation 11. Otherwise, Equation 10 is generally preferable.
  • the W quantities are respective weighting factors, and the default value is 1.
  • the respective weighting factor has a subscript of nds when referring to the neutron-density separation data, uran when referring to the uranium data, gr when referring to the gamma ray data, and rd when referring to the resistivity data. An operator may elect to weight the quantities differently, based on the observed geological conditions, data quality and/or provenance, or other factors.
  • Equations 6 through 11 are various measures of the sweet zone quality index based on individual well logs and normalizing constants.
  • SQI NDS refers to the sweet zone quality index from the neutron- density separation data
  • VWSH_NDS min is the minimum value of VWSH_NDS_NSBSL (with a default of zero).
  • SWI URAN refers to the sweet zone quality index from the uranium concentration data
  • URAN max refers to the maximum of the uranium concentration data (default value of 10 in ppm).
  • SQI GR refers to the sweet zone quality index from gamma ray data
  • GR max refers to the maximum of the gamma ray data (default value of 200 in API units).
  • SQI RD refers to the sweet zone quality index from resistivity data
  • RD max refers to the maximum of the resistivity data (default value of 100 in ohm-meter units).
  • SQI refers to the sweet gas quality index which is a combination of previous determined parameters from Equations 6 and 9, and either Equation 7 or Equation 8 depending on whether uranium concentration data is available.
  • the foregoing methods may be implemented in a computer system and computer executable instructions for performing the method may be stored on a tangible computer readable medium.
  • the system includes a data storage device or memory 202.
  • the stored data may be made available to a processor 204, such as a programmable general purpose computer.
  • the processor 204 may include interface components such as a display 206 and a graphical user interface 208, and is used to implement the above-described transforms in accordance with embodiments of the invention.
  • the graphical user interface may be used both to display data and processed data products and to allow the user to select among options for implementing aspects of the method.
  • Data may be transferred to the system 200 via a bus 210 either directly from a data acquisition device, or from an intermediate storage or processing facility (not shown).

Abstract

A computer system and computer implemented method for automatically identifying a hydrocarbon (such as kerogen, gas, oil) rich zone in a well bore includes obtaining well log data comprising neutron data, density data, radioactivity data, and resistivity data representative of physical characteristics of a formation surrounding the well bore and computing an apparent neutron porosity and an apparent density porosity based on the neutron data and density data. A normalized neutron-density separation is computed based on the computed apparent neutron porosity and the computed apparent density porosity and a baseline of the formation is determined for each data type. Using the computed normalized neutron-density separation, the radioactivity data, the resistivity data, and the determined baselines, the presence or absence of a hydrocarbon rich zone is determined. A quality index may further be derived from the data.

Description

SYSTEM AND METHOD FOR SWEET ZONE IDENTIFICATION IN SHALE GAS
RESERVOIRS
FIELD OF THE INVENTION
[0001] The present invention relates generally to methods and systems for identification of the sweet zone in shale gas reservoirs and more particularly to combining types of well log information to identify the sweet zone.
BACKGROUND OF THE INVENTION
[0002] Quick identification of the kerogen-rich sweet zone in wells, mapping the sweet zone areas, and placement of the horizontal holes within the sweet zone is one of the most important tasks in shale gas exploration and development. As shale gas plays have become more important to the oil and gas industry, methods of identifying kerogen-rich zones have gained in importance. In many cases, the existing methods are applicable only to the specific formation in which they have been applied, and do not have general relevance to new areas of exploration and development.
SUMMARY OF THE INVENTION
[0003] According to one implementation of the present invention, a computer implemented method for automatically identifying a hydrocarbon (such as kerogen, gas, oil) rich zone in a well bore includes obtaining well log data including neutron data, density data, radioactivity data, and resistivity data representative of physical characteristics of a formation surrounding the well bore and computing an apparent neutron porosity and an apparent density porosity based on the neutron data and density data. A normalized neutron-density separation is computed based on the computed apparent neutron porosity and the computed apparent density porosity and a baseline of normal shale is determined for each data type. Using the computed normalized neutron-density separation, the radioactivity data, the resistivity data, and the determined baselines, the presence or absence of a hydrocarbon rich zone is determined. A quality index may further be derived from the data. The computation of the presence or absence of a hydrocarbon rich zone and quality index is done at each depth level logged in the well. [0004] In an embodiment, a computer system for automatically identifying a hydrocarbon rich zone in a well bore includes a computer readable medium having computer readable well log data stored thereon, the well log data including neutron data, density data, radioactivity data, and resistivity data representative of physical characteristics of a formation surrounding the well bore. A processor of the computer system is configured and arranged to compute an apparent neutron porosity and an apparent density porosity based on the neutron data and density data, to compute a normalized neutron-density separation based on the computed apparent neutron porosity and the computed apparent density porosity, to compute a baseline of normal shale for the neutron data, density data, radioactivity data and resistivity data, and to compute the presence or absence of a hydrocarbon rich zone based on the computed normalized neutron-density separation, the radioactivity data, the resistivity data, and the determined baselines. The computations outlined above are done at each depth level logged in the well.
[0005] The above summary section is provided to introduce a selection of concepts in a simplified form that are further described below in the detailed description section. The summary is not intended to identify key features or essential features of the claimed subject matter, nor is it intended to be used to limit the scope of the claimed subject matter. Furthermore, the claimed subject matter is not limited to implementations that solve any or all disadvantages noted in any part of this disclosure.
BRIEF DESCRIPTION OF THE DRAWINGS
[0006] These and other features of the present invention will become better understood with regard to the following description, pending claims and accompanying drawings where:
[0007] Figure 1 is a flowchart illustrating a method in accordance with an embodiment of the invention;
[0008] Figure 2 is an example of a set of well logs showing a determined sweet zone indicator and sweet zone quality index in accordance with an embodiment of the invention; and
[0009] Figure 3 schematically illustrates a system for performing a method in accordance with an embodiment of the invention. DETAILED DESCRIPTION OF THE INVENTION
[0010] It is useful to evaluate subterranean formations to determine whether they are likely to contain significant amounts of organic matter, and therefore act as a good source of hydrocarbon resources. One method of characterizing a formation is to make measurements of characteristics along a borehole penetrating the formation, either during or after drilling operations, i.e., well logging. Well logging includes a number of techniques including resistivity/conductivity measurements, ultrasound, NMR, neutron, density, uranium concentration and radiation scattering, for example. Borehole data of this type is often used to replace or supplement the collection of cores for direct inspection. Conventionally, logged borehole data is analyzed by human interpreters in order to characterize a subsurface geological formation to allow decisions to be made regarding potential of the well or to determine information about the nature of the surrounding geologic area.
[0011] The inventors have determined that by combining information from a variety of well logs, a quantitative approach may be pursued to identify formations or portions of formations that are likely to be rich in organic material and therefore likely to offer potential in hydrocarbon production, without requiring human interpretation.
[0012] In this regard, a method in accordance with the present invention is illustrated in the flowchart of Figure 1. At step 10, well log data is obtained. In an embodiment, the well log data comprises neutron, density, uranium concentration and resistivity data. In another embodiment, uranium concentration is replaced by gamma ray data, each being a type of radioactivity data. As will be appreciated, the well log data may be acquired by any of a variety of well logging techniques, or may be existing well log data stored locally or remotely from a computer system on which the method is executed. In a particular example and not by way of limitation, the well log data may be from a shale formation.
[0013] From the density well log data, an apparent density porosity (PHIT D) is calculated in step 14. In this regard, Equation 1 sets out the calculation for PHIT D:
PHIT D = min(max(((pM - pB)/ (PM - PF)), 0.0), 1.0) Eqn. 1
[0014] With respect to Equation 1 , PM is the density of the rock matrix (where the matrix is selected to be a calcite matrix or other appropriate matrix, depending on the geology of the shale formation), PB is the bulk density of the rock, and PF is the density of fluid in the rock (where the fluid may be selected to be water). As will be appreciated, this Equation will produce a value of 0.0 where the ratio (PM - PB)/ (PM - PF) is negative, 1.0 when the ratio is greater than one, and the value of the ratio where the ratio is between zero and one. That is, it calculates a porosity value that is bounded by zero and one.
[0015] At step 12, an apparent neutron porosity (PHIT N) is calculated in accordance with Equation 2:
PHIT N = min(max(((TNPH - TNPM)/ (TNPF - TNPM)), 0.0), 1.0) Eqn. 2
[0016] In Equation 2, TMPH is the neutron porosity reading of the rock, TNPM is the neutron porosity of the matrix and TNPF is the neutron porosity of the fluid. Similarly to Equation 1 , this Equation produces a value equal to the ratio (TNPH - TNPM)/ (TNPF - TNPM) for values between zero and one, and is bounded by zero and one for all other values of the ratio.
[0017] Using the results of Equations 1 and 2, a value for normalized neutron-density separation (VWSH NDS) may be calculated (step 16) in accordance with Equation 3 :
VWSH NDS = max(min([(PHIT_N - PHIT D) - (PHIT N - PHIT_D)min]/ [(PHIT N - PHIT_D)ns - (PHIT N - PHIT_D)min], 1.0), -1.0) Eqn. 3
[0018] In Equation 3, the newly introduced quantity (PHIT N - PHIT_D)ns is the neutron-density separation for normal shales, while (PHIT N - PHIT_D)min represents a minimum value of the neutron-density separation. In an embodiment, (PHIT N - PHIT_D)min is taken to be zero and that portion of the numerator and denominator is eliminated. This equation produces values between minus one and one, although in most cases the values are between zero and one.
[0019] At step 18, a baseline value for each of the quantities is determined. For an embodiment using neutron, density, uranium concentration and resistivity data, baselines are determined for each of these. For embodiments in which gamma ray data replaces uranium concentration data, a baseline for gamma ray log readings is determined.
[0020] At step 20a, the values determined in the preceding steps are used to generate a sweet zone indicator (RNR) in accordance with the if statement in Equation 4. RNR = 1 if (VWSH NDS < VWSH NDS NSBSL · FVBSL and URAN > URAN NSBSL · FUBSL and RD > RD NSBSL · FRBSL) else RNR = 0 Eqn. 4
[0021] In Equation 4, by way of example and not limitation, VWSH NDS NSBSL is the normalized neutron-density separation baseline for normal shales, URAN is a uranium concentration, URAN NSBSL is baseline uranium concentration for normal shales, RD is a resistivity value of the log data, RD NSBSL is a baseline resistivity for normal shales and, FVBSL, FUBSL and FRBSL are adjustment factors for the respective baselines. Thus, if neutron-density separation is less than an adjusted baseline, and uranium and resistivity are above their respective adjusted baselines, then the indicator takes the value one, otherwise it takes the value zero.
[0022] As will be appreciated, the baseline for each type of log may be a constant, or may vary with depth and thus be represented by a curve or trendline, depending on the geological or borehole conditions. Typically, a shale interval is chosen to determine the baseline value or curve. The respective adjustment factors, FVBSL, FUBSL and FRBSL are selected to reduce measurement noise and also to reduce high frequency variations in the actual geological structure, thereby improving reliability of the indicator. In an embodiment, these are determined by Monte Carlo experimentation. The adjustment factors may also be adjusted in accordance with the experience of a user based on local geological conditions, analogues, and data quality and/or data provenance.
[0023] In alternate step 20b, for the case where uranium logs are replaced with gamma ray logs, Equation 4 is replaced by Equation 5.
RNR = 1 if (VWSH NDS < VWSH NDS NSBSL · FVBSL and GR > GR NSBSL · FGBSL and RD > RD NSBSL · FRBSL) else RNR = 0 Eqn. 5
[0024] The newly introduced quantities in Equation 5 are GR, which indicates gamma ray data, GR NSBSL which is the gamma ray baseline for normal shale and FGBSL, the adjustment factor for the gamma ray baseline. That is, for Equation 5, gamma ray data replaces the uranium data of Equation 4, but the equations otherwise operate in accordance with common principles. [0025] In general, the adjustment factors are selected to be close to one, and in an embodiment are limited to a range between 0.5 and 1.5. In a particular embodiment, (VSBSL, FUBSL, FRBSL, FGBSL) = (0.6, 0.99, 0.99, 0.99).
[0026] As will be appreciated, steps 12 and 14 could be performed in any order.
Likewise, the baseline determination for each type of well log performed in step 18 could, in principle, be performed in advance of any of the other calculations, and after all calculations except those of step 20, which depend on the results of step 18.
[0027] Evaluation of either Equation 4 or 5 will return a value of one or zero, indicating presence or absence of a sweet zone respectively. The indicator may then be used as a basis for determining a depth to initiate a horizontal drilling operation, or otherwise to guide production drilling decisions.
[0028] Figure 2 illustrates a number of well logs and derivative products in accordance with an embodiment of the invention. The first column shows radiation data derived from gamma ray measurements. The space between the two curves in the central portion of the log is indicative of uranium and represents a difference between spectral gamma radiation (the right hand curve) and computed gamma radiation (left hand curve). Some additional curves along the left-hand side of the trace are not relevant to the method described herein.
[0029] The second column shows depth of the well. The third column shows resistivity data for a number of different depths of investigation. The fourth column shows neutron and density data and the fifth shows uranium data.
[0030] In an embodiment, the indicator may be supplemented with a quality index that quantifies the quality of the identified sweet zone. This is illustrated in Figure 2, in the sixth column where 30 indicates a region in which the sweet zone indicator is one and 32 is a curve indicating the quality index within the zone 30.
[0031] As will be appreciated, in accordance with Equation 4, above, the region 30 corresponds to the shaded region in column 5 where normalized neutron-density separation is less than its baseline and the intersection of that shaded region with the shaded region in column 6 where uranium concentration is above its respective baseline. In the illustrated example, resistivity is above its baseline substantially throughout the region in which normalized neutron- density separation is less than its baseline.
[0032] In an embodiment, a sweet zone quality index may be calculated based on the data used to determine the sweet zone indicator. In particular, quality indexes are calculated for each of the data types, then those calculated quality indexes are used to compute an overall quality index that allows for comparison between or among various formations.
SQI NDS = min(max([VWSH_NDS_NSBSL - VWSH_NDS]/[VWSH_NDS_NSBSL -
Figure imgf000008_0001
SQI URAN = min(max(([URAN - URAN_NSBSL]/[URANmax-URAN_NSBSL]), 0), 1)
Eqn. 7
SQI GR = min(max([GR - GR_NSBSL]/[GRmax - GR NSBSL], 0), 1) Eqn. 8
SQI_RD = min(max([logio(RD)-logio(RD_NSBSL)]/[logio(RDmax)-logio(RD_NSBSL)],0), 1)
Eqn. 9
SQI = min(max([SQI_NDS · (Wnds/(Wnds + Wuran + Wrd)) + SQIJJRAN · (W^W^ + Wuran + Wrd)) + SQI_RD · (Wrd/(Wnds + + Wrd))], 0), 1)
Eqn. 10
Or
SQI = min(max([SQI_NDS · (Wnds/(Wnds + Wgr + Wrd)) + SQI_GR · (Wgr/(Wnds + Wgr +
Wrd)) + SQI_RD · (Wrd/(Wnds + Wgr + Wrd))], 0), 1) Eqn. 11
[0033] As will be appreciated, the choice between Equation 10 and 11 will depend on availability of uranium data. Where uranium data is not available, gamma ray data is used in accordance with Equation 11. Otherwise, Equation 10 is generally preferable. In Equations 10 and 11, the W quantities are respective weighting factors, and the default value is 1. The respective weighting factor has a subscript of nds when referring to the neutron-density separation data, uran when referring to the uranium data, gr when referring to the gamma ray data, and rd when referring to the resistivity data. An operator may elect to weight the quantities differently, based on the observed geological conditions, data quality and/or provenance, or other factors.
[0034] The newly introduced quantities in Equations 6 through 11 are various measures of the sweet zone quality index based on individual well logs and normalizing constants. In Equation 6, SQI NDS refers to the sweet zone quality index from the neutron- density separation data, and VWSH_NDSmin is the minimum value of VWSH_NDS_NSBSL (with a default of zero). In Equation 7, SWI URAN refers to the sweet zone quality index from the uranium concentration data, and URANmax refers to the maximum of the uranium concentration data (default value of 10 in ppm). In Equation 8, SQI GR refers to the sweet zone quality index from gamma ray data, and GRmax refers to the maximum of the gamma ray data (default value of 200 in API units). In Equation 9, SQI RD refers to the sweet zone quality index from resistivity data, and RDmax refers to the maximum of the resistivity data (default value of 100 in ohm-meter units). In Equations 10 and 11, SQI refers to the sweet gas quality index which is a combination of previous determined parameters from Equations 6 and 9, and either Equation 7 or Equation 8 depending on whether uranium concentration data is available.
[0035] In an embodiment, the foregoing methods may be implemented in a computer system and computer executable instructions for performing the method may be stored on a tangible computer readable medium.
[0036] A system 200 for performing the method is schematically illustrated in Figure
3. The system includes a data storage device or memory 202. The stored data may be made available to a processor 204, such as a programmable general purpose computer. The processor 204 may include interface components such as a display 206 and a graphical user interface 208, and is used to implement the above-described transforms in accordance with embodiments of the invention. The graphical user interface may be used both to display data and processed data products and to allow the user to select among options for implementing aspects of the method. Data may be transferred to the system 200 via a bus 210 either directly from a data acquisition device, or from an intermediate storage or processing facility (not shown). [0037] While in the foregoing specification this invention has been described in relation to certain preferred embodiments thereof, and many details have been set forth for purpose of illustration, it will be apparent to those skilled in the art that the invention is susceptible to alteration and that certain other details described herein can vary considerably without departing from the basic principles of the invention. In addition, it should be appreciated that structural features or method steps shown or described in any one embodiment herein can be used in other embodiments as well.

Claims

WHAT IS CLAIMED IS:
1. A computer implemented method of automatically identifying a hydrocarbon rich zone in a well bore, comprising:
obtaining well log data comprising neutron data, density data, radioactivity data, and resistivity data representative of physical characteristics of a formation surrounding the well bore;
computing an apparent neutron porosity and an apparent density porosity based on the neutron data and density data;
computing a normalized neutron-density separation based on the computed apparent neutron porosity and the computed apparent density porosity;
determining baselines of the formation for each of the neutron data, density data, radioactivity data and resistivity data;
using the computed normalized neutron-density separation, the radioactivity data, the resistivity data, and the determined baselines to determine the presence or absence of a hydrocarbon rich zone.
2. A method as in claim 1, wherein the radioactivity data is selected from the group consisting of uranium log data and gamma ray data.
3. A method as in claim 2, wherein, if uranium log data is available, it is preferentially selected over gamma ray data for use as the radioactivity data and if uranium log data is unavailable, gamma ray data is used as the radioactivity data.
4. A method as in claim 1 , wherein the using the computed normalized neutron- density separation, the radioactivity data, and the resistivity data to determine the presence or absence of a hydrocarbon rich zone comprises comparing each to a respective baseline and, based on the comparison, determining the presence or absence of the hydrocarbon rich zone.
5. A method as in claim 4, wherein the comparing comprises determining whether normalized neutron-density separation is less than its respective baseline, whether the radioactivity data is greater than its respective baseline and whether the resistivity data is greater than its respective baseline, such that when all three conditions are true, a hydrocarbon rich zone is determined to be present.
6. A method as in claim 4, wherein prior to the comparing, each of the respective baselines is multiplied by a respective adjustment factor.
7. A method as in claim 6, wherein each of the respective correction factors is between 0.5 and 1.5.
8. A method as in claim 1, further comprising, visually displaying an indicator of the presence or absence of a hydrocarbon rich zone.
9. A method as in claim 1, further comprising determining a sweet zone quality index by:
determining a respective quality index for each of the normalized neutron-density separation, the radioactivity data and the resistivity data; and
determining the sweet zone quality index by performing a weighted averaging of the respective quality indexes.
10. A method as in claim 1, wherein the formation comprises a shale gas reservoir.
11. A system configured and arranged to automatically identify a hydrocarbon rich zone in a well bore, comprising:
a computer readable medium having computer readable well log data stored thereon, the well log data comprising neutron data, density data, radioactivity data, and resistivity data representative of physical characteristics of a formation surrounding the well bore;
a processor, configured and arranged to:
compute an apparent neutron porosity and an apparent density porosity based on the neutron and density data;
compute a normalized neutron-density separation based on the computed apparent neutron porosity and the computed apparent density porosity; compute a baseline of the formation for the neutron, density, radioactivity and resistivity; and
compute the presence or absence of a hydrocarbon rich zone based on the computed normalized neutron-density separation, the radioactivity data, the resistivity data, and the determined baselines.
12. A system as in claim 11, further comprising a display, configured and arranged to produce a visual display of an indicator of the computed presence or absence of a hydrocarbon rich zone.
13. A system as in claim 11, wherein the processor is further configured and arranged to select uranium log data as the radioactivity data when uranium log data is available, or to select gamma ray data as the radioactivity data when uranium log data is unavailable.
PCT/US2011/044132 2010-09-13 2011-07-15 System and method for sweet zone identification in shale gas reservoirs WO2012036783A1 (en)

Priority Applications (7)

Application Number Priority Date Filing Date Title
CA2809969A CA2809969C (en) 2010-09-13 2011-07-15 System and method for sweet zone identification in shale gas reservoirs
CN2011800439070A CN103098062A (en) 2010-09-13 2011-07-15 System and method for sweet zone identification in shale gas reservoirs
EP11825591.8A EP2616978A1 (en) 2010-09-13 2011-07-15 System and method for sweet zone identification in shale gas reservoirs
AU2011302598A AU2011302598B2 (en) 2010-09-13 2011-07-15 System and method for sweet zone identification in shale gas reservoirs
JP2013528197A JP2013542412A (en) 2010-09-13 2011-07-15 System and method for sweet zone identification in a shale gas reservoir
EA201390369A EA201390369A1 (en) 2010-09-13 2011-07-15 SYSTEM AND METHOD FOR IDENTIFICATION OF A LOW-GAS GRID ZONE IN SLAST GAS PLASTERS
BR112013005708A BR112013005708A2 (en) 2010-09-13 2011-07-15 system and method for sweet zone identification in shale gas reservoirs

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US12/880,436 US8626447B2 (en) 2010-09-13 2010-09-13 System and method for sweet zone identification in shale gas reservoirs
US12/880,436 2010-09-13

Publications (1)

Publication Number Publication Date
WO2012036783A1 true WO2012036783A1 (en) 2012-03-22

Family

ID=45807523

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2011/044132 WO2012036783A1 (en) 2010-09-13 2011-07-15 System and method for sweet zone identification in shale gas reservoirs

Country Status (9)

Country Link
US (1) US8626447B2 (en)
EP (1) EP2616978A1 (en)
JP (1) JP2013542412A (en)
CN (1) CN103098062A (en)
AU (1) AU2011302598B2 (en)
BR (1) BR112013005708A2 (en)
CA (1) CA2809969C (en)
EA (1) EA201390369A1 (en)
WO (1) WO2012036783A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11360004B2 (en) 2017-09-15 2022-06-14 Korea Institute Of Geoscience And Mineral Resourues Shale gas extracting device and extracting method therefor

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8332155B2 (en) * 2010-09-13 2012-12-11 Chevron U.S.A. Inc. System and method for hydrocarbon gas pay zone characterization in a subterranean reservoir
EP2707828A4 (en) * 2011-05-10 2016-06-01 Chevron Usa Inc System and method for hydrocarbon pay zone definition in a subterranean reservoir
US9291690B2 (en) * 2012-06-22 2016-03-22 Chevron U.S.A. Inc. System and method for determining molecular structures in geological formations
CN104573339B (en) * 2014-12-24 2018-03-09 中国石油大学(北京) The geologic parameter of shale gas reservoir determines method and apparatus
US10209393B2 (en) 2015-01-23 2019-02-19 Halliburton Energy Services, Inc. Method to correct and pulsed neutron fan based interpretation for shale effects
CN106761728B (en) * 2017-02-14 2019-10-01 中国石油大学(北京) A kind of recognition methods of the advantageous interval of marine facies shale formation
KR101985497B1 (en) * 2017-06-09 2019-09-04 한국지질자원연구원 Method for estimating water saturation rate of tight gas reservoir composed of shale
CN112180443B (en) * 2019-07-04 2024-03-01 中国石油天然气集团有限公司 Shale gas two-dimensional seismic dessert area optimization method and device
CN111487176B (en) * 2020-05-13 2022-06-10 南京宏创地质勘查技术服务有限公司 Method for calculating porosity occupied by liquid hydrocarbon in shale oil system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030070480A1 (en) * 2001-10-11 2003-04-17 Herron Michael M. Real time petrophysical evaluation system
US20040099804A1 (en) * 2001-01-23 2004-05-27 Keyu Liu Oil reservoirs
US20080114547A1 (en) * 2004-04-30 2008-05-15 Schlumberger Technology Corporation Method and System for Determining Hydrocarbon Properties
WO2009015252A2 (en) * 2007-07-26 2009-01-29 Schlumberger Canada Limited System and method for estimating formation characteristics in a well

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3638484A (en) 1968-11-05 1972-02-01 Schlumberger Technology Corp Methods of processing well logging data
US4233839A (en) * 1979-01-15 1980-11-18 Schlumberger Technology Corporation Apparatus and method for determining characteristics of subsurface formations
US4686364A (en) 1985-07-19 1987-08-11 Schlumberger Technology Corporation In situ determination of total carbon and evaluation of source rock therefrom
US4916616A (en) 1986-12-08 1990-04-10 Bp Exploration, Inc. Self-consistent log interpretation method
FR2710987B1 (en) * 1993-10-06 1996-01-05 Schlumberger Services Petrol Combined logging device.
JPH10227868A (en) * 1997-02-12 1998-08-25 Anadrill Internatl Sa Method and apparatus for measurement of density of stratum
BR0111521A (en) * 2000-06-06 2004-04-20 Halliburton Energy Serv Inc Real time method for maintaining the stability of an underground formation
US6686736B2 (en) 2000-08-30 2004-02-03 Baker Hughes Incorporated Combined characterization and inversion of reservoir parameters from nuclear, NMR and resistivity measurements
JP2005127983A (en) * 2003-09-30 2005-05-19 Mitsubishi Heavy Ind Ltd Valuation methods of implant, underground resources, underground waste, underground cache, and geological structure, and in-building monitoring method, using hard ray or gamma ray
US7587373B2 (en) * 2005-06-24 2009-09-08 Halliburton Energy Services, Inc. Neural network based well log synthesis with reduced usage of radioisotopic sources
US7538547B2 (en) 2006-12-26 2009-05-26 Schlumberger Technology Corporation Method and apparatus for integrating NMR data and conventional log data
WO2009121131A1 (en) * 2008-03-31 2009-10-08 Southern Innovation International Pty Ltd Method and apparatus for borehole logging
US8573298B2 (en) 2008-04-07 2013-11-05 Baker Hughes Incorporated Method for petrophysical evaluation of shale gas reservoirs
CN101749012B (en) * 2008-12-08 2012-12-12 中国石油天然气集团公司 Method of determining oil reservoir exploitation level
CN101787884B (en) * 2010-01-28 2013-03-13 中国石油集团川庆钻探工程有限公司 Method for judging fluid type of reservoir through acoustic porosity-neutron porosity differential

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040099804A1 (en) * 2001-01-23 2004-05-27 Keyu Liu Oil reservoirs
US20030070480A1 (en) * 2001-10-11 2003-04-17 Herron Michael M. Real time petrophysical evaluation system
US20080114547A1 (en) * 2004-04-30 2008-05-15 Schlumberger Technology Corporation Method and System for Determining Hydrocarbon Properties
WO2009015252A2 (en) * 2007-07-26 2009-01-29 Schlumberger Canada Limited System and method for estimating formation characteristics in a well

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11360004B2 (en) 2017-09-15 2022-06-14 Korea Institute Of Geoscience And Mineral Resourues Shale gas extracting device and extracting method therefor

Also Published As

Publication number Publication date
JP2013542412A (en) 2013-11-21
AU2011302598B2 (en) 2015-04-16
BR112013005708A2 (en) 2016-05-10
US20120065887A1 (en) 2012-03-15
US8626447B2 (en) 2014-01-07
AU2011302598A1 (en) 2013-03-21
CN103098062A (en) 2013-05-08
CA2809969A1 (en) 2012-03-22
CA2809969C (en) 2019-01-15
EA201390369A1 (en) 2013-07-30
EP2616978A1 (en) 2013-07-24

Similar Documents

Publication Publication Date Title
US8626447B2 (en) System and method for sweet zone identification in shale gas reservoirs
CN108713089B (en) Estimating formation properties based on borehole fluid and drilling logs
CA3070868C (en) Resource density screening tool
US8332155B2 (en) System and method for hydrocarbon gas pay zone characterization in a subterranean reservoir
US10393641B2 (en) Methods of determining cementation exponent and saturation exponent in porous media from dielectric dispersion data
EP1800149A2 (en) Improving the accuracy of shaly sand formation evaluation
US20220082014A1 (en) Estimation of fluid saturation of a formation from integration of multiple well logs
AU2014290779B2 (en) System and method for estimating porosity distribution in subterranean reservoirs
Davies Permeability Modelling of a Sandstone Reservoir in Parts of the Niger Delta
US10013748B2 (en) Downhole imaging techniques using lamination analysis
Li et al. Petrophysical Rock Typing Based on Pore Geometry Improves Permeability and Bound Fluid Volume Estimation in Heterogeneous Sandstone Formations
Taylor Lithostratigraphic and Petrophysical Analysis of the Middle Devonian Marcellus Shale at the Mamont Prospect, Westmoreland County, Pennsylvania
CN114086938A (en) Gas saturation prediction method for heterogeneous sandstone reservoir
CN113311502A (en) Method and device for identifying conventional oil layer and shale oil layer in shale layer system
Samotorova et al. Exploring Petrophysical Uncertainties Thanks to Stochastic Well Log Data Interpretation: Termokarstovoye Field Case
CN115075800A (en) Hydrocarbon reservoir identification method, apparatus, medium and product
CN112443321A (en) Reservoir quality index-based oil layer identification method and device
Abrar et al. Flow Units Characterization and Application of Matrix Exponential Value (‘M-Value’) to Identify Secondary Porosity in Carbonates: A Successful Story to Identify By-Passed Oil and Increase Production
Tesch et al. Evaluation of methods to quantify pore pressure uncertainty from geophysical well logs
Sutiyono Reservoir Water Saturation and Permeability Modeling in the Pangkah Field

Legal Events

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

Ref document number: 201180043907.0

Country of ref document: CN

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

Ref document number: 11825591

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2011825591

Country of ref document: EP

ENP Entry into the national phase

Ref document number: 2809969

Country of ref document: CA

ENP Entry into the national phase

Ref document number: 2013528197

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2011302598

Country of ref document: AU

Date of ref document: 20110715

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 201390369

Country of ref document: EA

REG Reference to national code

Ref country code: BR

Ref legal event code: B01A

Ref document number: 112013005708

Country of ref document: BR

ENP Entry into the national phase

Ref document number: 112013005708

Country of ref document: BR

Kind code of ref document: A2

Effective date: 20130308