US 4,399,525 A · Method for interpreting well log records to yield indications of gas/oil in an earth formation such as a sandstone, limestone, or dolostone

Abstract

The present invention indicates that acoustic characteristics associated with a first well penetrating a gas/oil-containing strata such as a sandstone, limestone or dolostone, can be normalized (and favorably compared) with similar characteristics of a second, contiguous well of known hydrocarbon potential, such characteristics being calculated and displayed by means of a machine-implemented data processing method in which well logging and geologic data are fed thereto to calculate such characteristics without the need for shear-wave velocities. In more detail, in accordance with the invention brine-saturated bulk and shear moduli, (i.e., Kw* and Gw*) of a sandstone, limestone or dolostone penetrated by the second well can be predicted as a function of, say, brine-saturated P-wave modulus (Pw*) alone (independent of shear-wave velocity). In that way, resulting acoustic values including distinct velocity and amplitude values as a function of a saturation operator can ultimately be provided. Such values, when compared to actual, in-well, field-generated characteristics of the first well, are surprisingly accurate predictors of the amount of gas/oil maturation in and around the first well. The method has particular accuracy in designating gas zones within formations of interest due to use of the complete Benedict-Webb-Rubin non-ideal gas law in pseudo-reduced form.

Claims21

Dependent claims

System for resolving event characteristics of an acoustic log of a first well so as to determine gas/oil fractions in a selected zone of an earth formation by creating synthetic gas-indicating characteristics from petrophysical parameters corresponding to said selected zone of interest associated with a second well of substantially zero hydrocarbon potential, comprising:

(a) means for obtaining said acoustic log for said first well;
(b) means for determining actual compressional velocity and density values of said zone of interest based on field information associated with said second well in which selection of a zonal lithology from a group comprising sandstone, limestone and dolostone occurs;
(c) a digital programmable computer for:
(i) predicting brine-saturated bulk and shear moduli (Kw* and Gw*) using empirical equations containing coefficients and said P-wave modulus value (Pw*), said coefficients varying in empirical fashion as a function of at least porosity, pressure and temperature;
(ii) estimating a frame modulus value (Ka) of said zone of interest using a relationship involving Kw*, Kb, Km and .phi., where Kw* is brine-saturated bulk density of the aggregate, Kb is the bulk modulus of brine, and Km is the bulk modulus of the solid material making up the zone of interest and .phi. is porosity;
(iii) predicting bulk modulus K* values with the gas/oil present in said zone of interest in differing amounts using a relationship involving Ka of step (ii), Km, .phi., Kf, where Kf is the bulk modulus of a liquid-gas mixture of a selected oil/gas gravity and saturation;
(iv) assuming Gw* is equal to Ggas/oil*, predicting P-wave modulus values of the zone with gas/oil present (Pgas/oil*) based on the relationship
P.sub.gas/oil* =K.sub.gas* +(4/3)G.sub.w*
(v) predicting acoustic velocities within said zone from said P-wave modulus values (Pgas/oil*) of step
(iv) as a function of said selected gas/oil gravity and saturation;
and
(d) comparitor means comparing said predicted velocities or derivations thereof with actual velocities or derivations thereof, of (a) associated with said acoustic log of said first well to predict presence and amounts of gas/oil fractions in said zone of interest at said first well.

9 more elements

System for resolving event characteristics of a sonic log of a first well so as to determine gas/oil fractions in a selected zone of an earth formation by creating synthetic gas-indicating characteristics from brine-saturated petrophysical parameters corresponding to said selected zone of interest including a programmable digital computer having a CPU-control means, program and data memories connected to said CPU-control means and display-storage means also connected to said CPU-control means, comprising:

(a) means for obtaining said acoustic log for said first well;
(b) means for determining actual compressional velocity and density values of said zone of interest based on field information associated with said second well in which selection of a zonal lithology from a group comprising sandstone, limestone and dolostone, occurs;
(c) said CPU-control means of said digital computer operating on data within said data memory under control of program statements initially held in said program memory but sequentially fetched to said CPU-control means so as:
(i) from brine-saturated bulk and shear moduli (Kw* and Gw*), calculated by empirical equations including a true P-wave modulus value (Pw*) corresponding to said zone of interest based on said actual velocity and density values of (a) and estimated frame modulus values (Ka) of said zone of interest, to predict bulk modulus values K* with gas/oil present in said zone of interest in differing amounts;
(ii) to predict P-wave modulus values of the zone with gas/oil present (Pgas/oil*) based on the relationship
P.sub.gas/oil* =K*.sub.gas/oil +(4/3)G*.sub.w
(iii) to predict sonic velocities within said zone (Pgas/oil*) from said P-wave modulus values of step (ii) as a function of gas/oil saturation;
(iv) to compare said predicted sonic velocities or derivations thereof with actual velocities or similar derivatives thereof, of (a) associated with said log of said first well to predict amounts of gas/oil fractions in said zone of interest at said first well.

5 more elements

System for resolving event characteristics of a sonic log of a first well so as to determine gas/oil fractions in a selected zone of an earth formation by creating synthetic gas/oil-indicating characteristics from brine-saturated petrophysical parameters corresponding to said selected zone of interest but associated with a second well of substantially zero hydrocarbon potential, comprising:

(a) means for obtaining said acoustic log for said first well;
(b) means for determining actual compressional velocity and density values of said zone of interest based on field information associated with said second well in which selection of a zonal lithology from a group comprising sandstone, limestone and dolostone, occurs;
(c) a programmable digital computer for:
(i) from calculated normalized brine-saturated petrophysical parameters involving data of (b) predicting compressional velocities of said zone with gas/oil in differing amounts present in said zone, in accordance with the relationship ##EQU23## where rho-b is the bulk density, .phi. is the porosity, G* is the shear modulus at the saturations of interest, K.sub.f equals the bulk modulus for the liquid-gas mixture, K.sub.a is the frame modulus, K.sub.m is the bulk modulus of a hypothetical rock, i.e., same composition but perfectly cemented with zero porosity;
(ii) from sonic velocities of said zone calculated from (i), predicting synthetic sonic response information with gas/oil present in said zone;
and
(iii) comparing said predicted sonic response information with gas/oil present with said actual field-acquired sonic logging response data of said sonic log of (A) so as to determine both presence and amounts of gas/oil within said zone of interest at said first well.

4 more elements

Method for resolving event characteristics of an acoustic log of a first well so as to determine gas/oil fractions in a selected zone of an earth formation by creating synthetic gas-indicating characteristics from petrophysical parameters associated with a second well of substantially zero hydrocarbon potential, comprising:

(a) obtaining said acoustic log for said first well by field measurement;
(b) determining actual compressional velocity and density values of said zone of interest based on field information associated with said second well in which selection of a zonal lithology from a group comprising sandstone, limestone and dolostone, occurs;
(c) predicting brine-saturated bulk and shear moduli (Kw* and Gw*) for said lithology of step (b) using empirical equations containing coefficients and a P-wave modulus value (Pw*), said coefficients varying in empirical fashion as a function of at least porosity, pressure and temperature, said P-wave modulus value (Pw*) being a true P-wave modulus value which corresponds to said zone of interest, based on data of (b) associated with said second well;
(d) estimating a frame modulus value (Ka) of said lithology of said zone of interest using a relationship involving Kw*, Kb, Km and .phi., where Kw* is brine-saturated bulk density of the aggregate, Kb is the bulk modulus of brine, Km is the bulk modulus of the solid material making up said formation at said zone of interest and .phi. is porosity;
(e) predicting a bulk modulus value K* with the gas/oil present at said zone of interest using a relationship involving Ka of step (d), Km, .phi., Kf, where Kf is the bulk modulus of a gas/oil/brine mixture of a selected gas/oil gravity and saturation;
(f) assuming Gw* is equal to Ggas/oil*, predicting a P-wave modulus value of the zone with gas/oil present (Pgas/oil*) based on the relationship
P.sub.gas/oil* =K*.sub.gas/oil +(4/3)G*.sub.w
(g) predicting at an acoustic velocity at said zone from said P-wave modulus values (Pgas/oil*) of step (f);
(h) repeating steps (e)-(f) for a series of gas/oil-brine mixtures of different gas/oil gravities and saturations to provide a suite of synthetic acoustic velocities at said zone of interest as the function of a plurality of gas/oil gravities and saturations;
and
(i) comparing said suite of estimated velocities with velocity information associated with said actual acoustic log of (a) of said first well to accurately predict presence and amounts of gas/oil fractions at said zone of interest.

8 more elements

Method for resolving event characteristics of an acoustic log of a first well so as to determine gas/oil fractions at a selected zone of an earth formation by creating synthetic gas-indicating characteristics from petrophysical parameters corresponding to said selected zone of interest of substantially zero hydrocarbon potential, comprising the steps of:

(a) obtaining said acoustic log for said first well by field measurement;
(b) determining actual compressional velocity and density values of said zone of interest based on field information associated with said zone in and about its area of zero hydrocarbon potential in which selection of a zonal lithology from a group comprising sandstone, limestone and dolostone, occurs;
(c) from brine-saturated bulk and shear moduli (Kw* and Gw*), calculated by empirical equations including a true P-wave modulus value (Pw*) which corresponds to said zone of interest based on data of (b) and estimated frame modulus values (Ka) at said zone of interest, predicting bulk modulus K* with the gas/oil present at said zone of interest in said first well in varying amounts;
(d) predicting P-wave modulus values at said zone with gas/oil present (Pgas/oil*) in said varying amounts based on the relationship
P.sub.gas/oil* =K*.sub.gas/oil +(4/3)G*.sub.w
(e) predicting acoustic velocities at said zone in said first well from step (d) as a function of gas/oil saturation;
(f) indexing predicted and actual velocity values of (a) as a function of different gas/oil saturation and gravities;
and
(g) determining a best match of all indexed values so as to predict amounts of gas/oil at said zone of interest.

6 more elements

Method for resolving events of an acoustic log of a first well associated with at least a selected zone of an earth formation so as to determine gas/oil fractions therein by creating synthetic gas/oil-indicating characteristics from brine-saturated petrophysical parameters corresponding to said selected zone of interest at a second well of substantially zero hydrocarbon potential, comprising the steps of:

(a) obtaining said acoustic log for said first well by field measurement;
(b) determining actual compressional velocity and density values of said zone of interest based on field information associated with said second well in which selection of a zonal lithology from a group comprising sandstone, limestone and dolostone, occurs;
(c) from brine-saturated bulk and shear moduli (Kw*) and Gw*) calculated from empirical equations including a true brine-saturated P-wave modulus value (Pw*) which corresponds to said zone of interest based on data of (b) associated with said second well, and estimated frame modulus values (Ka) of said zone of interest, predicting a plurality of bulk moduli (K*) with gas/oil present in said zone of interest in preselected amounts;
(d) predicting a multiplicity of P-wave moduli of said zone with gas/oil present (Pgas/oil*) in said preselected amounts based on the relationship
P.sub.gas/oil* =K*.sub.gas/oil +(4/3)G*.sub.w
where Pgas/oil* is a multiplicity of predicted P-wave moduli with gas/oil present in said preselected amounts;
Kgas/oil* is a plurality of predicted bulk moduli of said zone with gas/oil present in said preselected amounts;
Gw* is the brine-saturated shear modulus of said zone;
(e) predicting synthetic acoustic values from acoustic velocities within said zone calculated from step (d) as a function of different gas/oil gravities and saturations;
and
(f) comparing said synthetic values indexed to different gas/oil gravities and/or saturations with actual field-produced values produced from said acoustic log of (a) to predict amounts of gas/oil in said zone of interest at said first well.

8 more elements

Method for resolving event characteristics of an acoustic log of a first well so as to determine gas/oil fractions in a selected zone of an earth formation by creating synthetic gas/oil-indicating characteristics from brine-saturated petrophysical parameters corresponding to said selected zone of interest located at a second well of substantially zero hydrocarbon potential, comprising the steps of:

(a) obtaining said acoustic log for said first well by field measurement;
(b) determining actual compressional velocity and density values of said zone of interest based on field information associated with said second well in which selection of a zonal lithology from a group comprising sandstone, limestone and dolostone, occurs;
(c) from calculated normalized brine-saturated petrophysical parameters involving a true brine-saturated P-wave modulus value obtained from (b) corresponding to said zone of interest at said second well, predicting velocities of said zone with gas/oil in differing amounts present in said zone, in accordance with the relationship ##EQU28## where rho-b is the bulk density, .phi. is the porosity, G* is the shear modulus at the saturations of interest, K.sub.f equals the bulk modulus for the liquid-gas mixture, K.sub.a is the frame modulus, K.sub.m is the bulk modulus of a hypothetical rock, i.e., same composition but perfectly cemented with zero porosity;
and
(d) comparing said predicted velocity information with actual field-acquired data associated with said acoustic log of (a) associated with said first well, so as to determine both presence and amounts of gas/oil within said zone of interest at said first well.

2 more elements

Dependencies are inferred from the claim text.

Method for interpreting well log records to yield indications of gas/oil in an earth formation such as a sandstone, limestone, or dolostone

Life of this patent

Method for interpreting well log records to yield indications of gas/oil in an earth formation such as a sandstone, limestone, or dolostone

Life of this patent

Abstract

Claims21