The present invention relates to a method of predicting the pressure sensitivity of seismic velocity within reservoir rocks. More generally, the invention relates to rock physics and the modeling of static and dynamic reservoir properties. The invention also relates to an heuristic approach for interpreting seismic activity within cemented sandstone reservoirs.
The general problem addressed by the present invention is how to predict the composition of a rock formation and, in particular, whether (and to what extent) it is saturated with oil, from seismic velocity measurements. In order to interpret seismic surveys it is necessary to establish relationships between the measured velocities and the intrinsic rock properties.
It is known that seismic compression (p) wave velocity is strongly dependent on an effective rock pressure. The effective pressure is the difference between the confining pressure (of the overlying rock column) and the pore pressure (which may be equal to, greater than or less than the hydrostatic pressure).
In general, velocity rises with increasing confining pressure and levels off (to a terminal velocity) when the effective pressure is high. This effect is thought to be due to crack closure: at low effective pressure cracks are open and easily closed by an increase in pressure (resulting in a small bulk modulus, K, and low velocity); as the effective pressure increases the cracks are all closed (resulting in an increase in K and in velocity).
Static rock physics modeling may be used to generate 3-Dimensional (3D) data plots of rock properties at a particular instance in time. Dynamic rock physics modeling, on the other hand, provides tools for estimating the evolution of rock properties over time. This is also referred to as 4-Dimensional (4D) modeling, where the fourth dimension represents time.
Rock physics models for fluid and stress dependency in reservoir rocks have been found to be essential for the quantification and interpretation of 4D seismic signatures during reservoir depletion and injection. However, our ability to predict the sensitivity of seismic data to pressure from first principles is poor.
The current state of the art requires that we calibrate the pressure dependence of seismic velocity with core measurements. A major challenge is the fact that consolidated rocks often break up during coring, and hence the pressure sensitivity is likely to be over-predicted in the laboratory relative to in-situ conditions. For unconsolidated sands, acquisition of core samples is not usually feasible due to the friable nature of the sediments.
One physical model that has been applied to predict pressure sensitivity in unconsolidated granular media is the Hertz-Mindlin contact theory (as described in, for example, Avseth et al., 2005, “Quantitative Seismic Interpretation; Applying Rock Physics Tools to Reduce Interpretation Risk”, Cambridge University Press). Several other empirical models have also been suggested (e.g. Bachrach and Avseth, 2008, “Rock physics modeling of unconsolidated sands: Accounting for non-uniform contacts and heterogeneous stress fields in the effective media approximation with applications to hydrocarbon exploration”, Geophysics, 73, E197-E209) which have fitting parameters that correlate with the micro-crack intensity, soft porosity and aspect ratio of the rock, and feasibility studies can be undertaken based on assumptions about these parameters. However, these models are not easily applied to moderately consolidated sandstones with contact cement, where crack parameters and aspect ratios are difficult to quantify.
In order to predict the seismic velocities of a rock at a given point in time, knowing only the porosity, mineralogic composition, and the elastic moduli of the mineral constituents, we can at best predict the upper and lower bounds of the seismic velocities. However, if we know the geometric details of how the mineral grains and pores are arranged relative to each other, we can predict more exact seismic properties using static rock physics modelling. There are several models that account for the microstructure and texture of rocks and these, in principle, allow us to go the other way: to predict the type of rock and microstructure from seismic velocities. This rock physics diagnostic technique was introduced by Dvorkin and Nur in 1996 as a means to infer rock microstructure from velocity-porosity relations. This technique is conducted by adjusting an effective-medium theoretical model curve to a trend in the seismic data, assuming that the microstructure of the rock matches that used in the model.
A main short-coming with the Dvorkin-Nur contact-cement model is that is does not include pressure sensitivity. Instead, it is assumed that the cemented grain contacts immediately loose pressure sensitivity as the cementation process initiates. However, from in-situ observations, we know that cemented reservoirs can have significant pressure sensitivity. This could either be related to fractures not captured by the microstructural scale model, or by a patchy cementation where some grain contacts are cemented and others are loose.
As with the Hertz-Mindlin contact theory for loose granular media, the Dvorkin-Nur contact-cement model is also found to often overpredict shear stiffness in cemented sandstones. This could be related to non-uniform grain contacts and tangential slip at loose contacts, associated heterogeneous stress chains, and/or relative roll and torsion not taken into account in the contact theory. A reduced shear factor (Ft) has been introduced to honour this “reduced shear effect” in the contact theory and this varies between 0 and 1 representing the boundary conditions between no-friction (Walton smooth contact theory) and no-slip (Walton rough or Hertz-Mindlin contact theory) conditions, and for loose sands this parameter can be estimated directly from dry rock using Poisson's ratio. For cemented sandstones, this parameter is a pure fitting parameter, yet it has been found to correlate with the degree of cementation.
More details of the models discussed above (e.g. the Hertz-Mindlin contact theory model, the Walton smooth pressure sensitive model, the Hashin-Shtrikman model and the Dvorkin-Nur contact-cement model) are provided in “The Rock Physics Handbook: Tools for Seismic Analysis of Porous Media” by Gary Mavko, Tapan Mukerji and Jack Dvorkin.
One way to heuristically quantify the pore stiffness of a rock is to measure the distance between an upper and a lower elastic bound at a given pressure. Marion and Nur introduced the Bounding Average Method as a relative measure of pore stiffness of a rock and this is illustrated in
A similar approach has been suggested to quantify the degree of consolidation, and to define a weight function, W, depending on where the rock data (e.g. sandstone data) is plotted between an upper and a lower bound in the elastic moduli versus porosity domain and this is illustrated in
A similar weight factor can also be calculated from the dry shear modulus data in
It is therefore an aim of the present invention to provide a method of predicting the pressure sensitivity of seismic velocity within reservoir rocks, which ameliorates at least some of the afore-mentioned problems.
It is proposed here to employ a method of predicting the pressure sensitivity of seismic velocity within reservoir rocks, and which comprises the following steps:
Thus, embodiments of the present invention provide a method which takes into account the level of cementation of the rock and which provides a suitable model for interpreting seismic velocity data (obtained from or modelled for a particular rock formation) and which includes pressure sensitivity in the model to an appropriate degree. An advantage of the present method is that the pressure sensitivity of partially cemented rock is accounted for.
It is noted that the terms ‘friable sand’, ‘loose sand’, ‘unconsolidated sand’ and ‘unconsolidated rock’ are used interchangeably throughout this specification. In addition, the terms ‘uncemented’, ‘partially cemented’ and ‘cemented’ have been used in a manner synonymous with the terms ‘unconsolidated’, ‘partially consolidated’ and ‘consolidated’, respectively.
The method of the present invention may be considered a ‘hybrid’ model as it combines existing models for unconcolidated and consolidated sands, respectively, into a model with can be used for partially consolidated sands.
The method is particularly suitable for use in relation to sandstone rock formations, although it may be applied to other rock types also.
The dry rock moduli and porosity may be obtained based upon well log data, i.e. data measured at specific well locations. Alternatively, the dry rock moduli and porosity may be obtained from a theoretical model of the rock geology.
Although data from core samples may be used to provide the dry rock moduli and porosity, it is an advantage of the present invention that the input data can be obtained without requiring such core samples.
In specific embodiments, well log data may be employed to calibrate one or more of the models.
The method may be used to quantify the sensitivity of seismic velocity within rock to pressure.
The method may be used to determine a pressure or pressure change within rock using seismic survey results.
The method may involve creating cross-plots of dry rock moduli versus porosity, including elastic bounds for different degrees of consolidation.
The effective rock moduli and corresponding seismic velocities as a function of pressure for a partially cemented rock may be estimated by a weighted average of the friable sand and cemented models.
The method may be used to quantify the sensitivity of seismic velocity to pressure in cemented sandstones, without the need for core measurements.
The method may be used during mapping of reservoir pressure from 3-D and 4-D seismic data.
The method may comprise the step of using the predicted pressure sensitivity to interpret seismic velocity data and thereby predict the composition of a rock formation.
All three models may be required in order to determine the degree of cementation, however, depending on the degree of cementation, one or more of the models may be required in order to determine the sensitivity of seismic velocity to pressure. Accordingly, if it is determined that the rock under consideration comprises only one of friable sand, partially cemented rock or cemented rock, the relevant model may be employed on the entire dataset. However, if it is determined that the rock under consideration comprises portions of more than one of friable sand, partially cemented rock or cemented rock, the relevant models may only be employed on the relevant portions of the dataset.
The first model for rock comprising friable sand may comprise or be based upon the Hertz-Mindlin model for unconsolidated sands. In certain embodiments, the first model may comprise the Walton Smooth contact theory model.
The second model for rock comprising partially cemented rock may comprise a modified contact model that is pressure sensitive. More specifically, the second model may comprise the Walton smooth pressure sensitive model or the Hertz-Mindlin model (defining a Hashin-Shtrikman soft bound) in combination with the Dvorkin-Nur contact cement model or the Constant Cement model (defining a Hashin-Shtrikman stiff bound).
The third model for rock comprising cemented rock may comprise or be based upon the Dvorkin-Nur contact cement model for consolidated sands or the Constant Cement model.
The degree of cementation may be determined by modelling upper and lower elastic bounds based on the porosity of the rock and then establishing the weighting function to account for the degree of cementation of the rock. Elastic bounds may be determined separately for bulk modulus and shear modulus data. Thus, different weighting functions (WK and WG) may be obtained in relation to the bulk modulus and the shear modulus data.
The lower (soft) bound may be determined using the Hertz-Mindlin model or the Walton Smooth model in combination with the Hashin-Shtrikman model to determine the relationship between the elastic moduli and porosity for unconsolidated sands at a given (in situ) pressure.
The upper (stiff) bound may be determined using the Dvorkin-Nur contact cement model or Constant Cement model in combination with the Hashin-Shtrikman model to determine the relationship between the elastic moduli and porosity for consolidated sands (i.e. where all grains are taken to be cemented). In practice, the upper bound is determined at the degree of cementation where the rock is substantially non-compressible and therefore the pressure sensitivity is determined to be negligible. Accordingly, the upper bound may be defined at a hypothetical maximum pressure where no further stress sensitivity will be seen.
In certain embodiments, the upper bound may be obtained by matching the Hertz-Mindlin model (or Walton smooth pressure sensitive model) with the lower bound Hashin-Shtrikman model at a very high pressure such that the upper bound superimposes onto a Constant Cement model (or Dvorkin-Nur contact cement model) for the cemented rock.
In particular embodiments, the degree of cementation at which the rock is considered substantially non-compressible and is therefore considered cemented rock may be 10%. In other embodiments, the degree of cementation at which no further stress sensitivity is observed may be 8%, 12% or another value obtained through further modelling and/or experimentation.
The Applicants envisage conducting further studies to determine uncertainties related to the upper and lower bounds and, if possible, to find ways to improve on the accuracy of the upper and lower bounds.
The weighting function may be linear and may vary between 0 representing no cementation and 1 representing the degree of cementation at which the rock is substantially non-compressible and therefore all grain contacts are taken to be cemented. The bulk modulus weighting function, WK, may be calculated from Equation (2) below, where Kdry is the pressure sensitive dry bulk modulus (which has been modelled or observed) at porosity (P0), Ksoft is the pressure sensitive soft (i.e. lower bound) bulk modulus at the same porosity (P0), and Kstiff is the pressure insensitive stiff (i.e. upper bound) bulk modulus at this porosity value.
Similarly, the shear modulus weighting function, WG, may be calculated from Equation (3) below, where Gdry is the pressure sensitive dry shear modulus (which has been modelled or observed) at porosity (P0), Gsoft is the pressure sensitive soft (i.e. lower bound) shear modulus at the same porosity (P0), and Gstiff is the pressure insensitive stiff (i.e. upper bound) shear modulus at this porosity value.
Thus, in the above embodiments, the linear weighting functions are used to interpolate between the soft and stiff bounds. In other embodiments, the weighting functions may be non-linear and may be obtained using other methods. For example, a two-step Hashin-Shtrikman modeling approach may be employed whereby a first interpolation is performed between uncemented and cemented end members at a high porosity, followed by a second interpolation between the high porosity and low porosity (mineral point) end members.
The weighting functions allow us to estimate vertical pressure sensitivity in partially cemented structures. The pressure dependence of the dry bulk and shear elastic moduli may be obtained using equations (4) and (5) below.
K
dry(Peff)=(1−WK)·Ksoft(Peff)+WK·Kstiff (4)
G
dry(Peff)=(1−WG)·Gsoft(Peff)+WG·Gstiff (5)
From these dry elastic moduli we may use known techniques to calculate the expected seismic velocities and acoustic impedances at various pressures (taking into account the effect of pore fluid using Gassmann's (1951) equations). For example, we may determine the stress sensitivity on these parameters for rock saturated with gas, oil or brine. We may then map well log data against the modelled data and correlate the results so as to establish the rock properties giving rise to the well log data. This therefore serves to indicate the likely rock structure producing the recorded seismic data. The method may also be used to determine how the properties of the rock have changed over time by comparing the results of seismic data obtained at one time to that obtained at another time.
The upper bound may be modelled first since it is both data and pressure independent. The lower bound may be modelled second as it is data independent but pressure dependent. An initial value for the pressure dependence may be input into the model for the lower bound from well log data or a geological model (e.g. where a cement volume is assumed). Well log data or modelled data may then be applied between bounds and the weighting functions and resulting pressure dependence calculated as described above before the data is converted into velocity plots for comparison with measured seismic data.
In some embodiments, regression modelling may be employed on a simulated dataset so as to derive equations for calculating seismic velocities directly from porosity, effective pressure and cement volume values.
In accordance with a second aspect of the invention, there is provided a computer system configured to carry out the above method.
In accordance with a third aspect of the invention, there is provided a computer program, comprising computer readable code which, when run on a computer system causes the computer system to carry out the above method.
In accordance with a fourth aspect of the invention, there is provided a computer program product comprising a computer readable medium and a computer program according to the third aspect of the invention, wherein the computer program is stored on the computer readable medium.
Embodiments of the invention will now be described with reference to the accompanying drawings, in which:
Note that the selected Gullfaks data plots along the more stress sensitive curves than the Statfjord data. This is likely related to the degree of consolidation.
An embodiment of the present invention there is provided a method of predicting the pressure sensitivity of seismic velocity within sandstone reservoir rocks, which comprises: defining the degree of cementation of rock as at least one of friable sand, partially cemented rock comprising a degree of cementation up to a level at which the rock is substantially non-compressible (in this example, up to 10% cementation), and cemented rock comprising a degree of cementation at which the rock is substantially non-compressible (in this example, 10% cementation and above). For rock comprising friable sand we define a first model specifying a dependence of seismic velocity upon pressure. For rock comprising partially cemented rock we define a second model specifying a dependence of seismic velocity upon pressure and a weighting function accounting for a degree of cementation of the rock. For rock comprising cemented rock we define a third model demonstrating an insensitivity of seismic velocity to pressure. For a given dry rock moduli and porosity, we determine a degree of cementation, select the appropriate model, and use the selected model to predict the sensitivity of seismic velocity to pressure.
In this embodiment, the first model employed for rock comprising friable sand is the Hertz-Mindlin contact theory model. The second model employed for rock comprising partially cemented rock comprises the Hertz-Mindlin contact theory model (defining a Hashin-Shtrikman soft bound) in combination with the Dvorkin-Nur contact cement model (defining a Hashin-Shtrikman stiff bound). The third model employed for rock comprising cemented rock is the Dvorkin-Nur contact cement model for consolidated sands.
The degree of cementation is determined by modelling upper and lower elastic bounds based on the porosity of the rock and then establishing the weighting function to account for the degree of cementation of the rock. The elastic bounds are determined separately for bulk modulus and shear modulus data such that different weighting functions (WK and WG) are obtained in relation to the bulk modulus and the shear modulus data. It should be noted that the weighting functions will be different for the bulk and shear moduli because the relative location of the elastic bounds will be affected by the reduced tangential shear stiffness mentioned above.
The lower (soft) bound is determined using the Hertz-Mindlin contact theory model in combination with the Hashin-Shtrikman model to determine the relationship between the elastic moduli and porosity for unconsolidated sands at a given (in situ) pressure.
The upper (stiff) bound is determined using the Dvorkin-Nur contact cement model in combination with the Hashin-Shtrikman model to determine the relationship between the elastic moduli and porosity for consolidated sands (i.e. having at least 10% cementation such that the effect of pressure is taken to be equivalent to that where all grains are cemented).
In this particular embodiment, the upper bound is estimated by matching the Hertz-Mindlin model with the lower bound Hashin-Shtrikman model at a very high pressure such that the upper bound superimposes onto a Constant Cement model for 10% cementation. This approach provides a soft and stiff bound with the same basic shape, which gives a more stable and realistic weighting function estimation for a given porosity value.
The weighting functions in the present embodiment are calculated in accordance with equations (2) and (3) above. The pressure sensitivity of the dry bulk and shear elastic moduli is then derived from the weighting functions in accordance with equations (4) and (5) above. Thus, by combining the Hertz-Mindlin contact theory model for unconsolidated sands with a stiff contact cement model, we can obtain a modified contact model for heterogeneous contacts that is pressure sensitive via the fraction of unconsolidated grain contacts.
Any sandstone data point (e.g. from well log data) can then be inverted for the weighting factors WK, WG allowing us to estimate stress curves for each data point.
In a particular example, the applicants simulated synthetic data for a wide range of porosities and cement volumes using the static rock physics models described above. The cement volume was varied between 0 and 10%, as it was assumed that if the cement volume was higher than this there will be no stress sensitivity at the grain contacts. Porosity was varied between 0 and 0.4 (i.e. 40%) and a synthetic elastic moduli dataset was created covering all possible combinations of porosity and cement volume within these ranges and the results are shown in
Next, the applicants estimated the soft and stiff bounds that enclose the simulated dataset in the moduli-porosity domain. These bounds were modeled by combining Hertz-Mindlin and lower bound Hashin-Shtrikman as described above. The soft bound is the unconsolidated sand model where the reference effective stress (P0) is set to 20 MPa in this example. This represents the effective stress at around 2 km burial depth, which is the depth we expect quartz cementation to initiate in the North Sea. Any data point falling on this soft bound should represent unconsolidated sands where all grain contacts are stress-sensitive.
The stiff bound is defined by increasing the effective stress in the Hertz-Mindlin model so that it mimics the 10% constant cement model. In this example, we find that an effective stress of 600 MPa must be selected in order to get this match. For any practical reason, this stiff bound should therefore be considered what happens when all grain contacts are closed, and there is no stress sensitivity in the sandstone data falling on this bound. A linear weighting function is then defined between the soft and the stiff bounds (e.g. using Equations 2 and 3), both for bulk modulus and shear modulus versus porosity. This weighting function will define the stress sensitivity of the cemented sandstone. In this example, the soft bound is defined with a reduced shear factor Ft=0. The reduced shear factor for the stiff bound is set to 0.5. It has been demonstrated that this parameter is depth dependent and likely associated with degree of diagenesis. This parameter may therefore be further updated in an iterative scheme to fit a calibration data set (not demonstrated here).
Next, we derive the effective bulk and shear modulus as a function of effective stress for the cemented sandstones, depending on the estimated weight factors (i.e. consolidation degree).
From these dry elastic moduli the applicants used known techniques to calculate the expected seismic velocities and acoustic impedances at various pressures (taking into account the effect of pore fluid using Gassmann's (1951) equations). Accordingly,
We then map well log data against the modelled data and correlate the results so as to establish the rock properties giving rise to the well log data.
Simplified Regression Models and Dynamic Rock Physics Templates
In another embodiment, regression modelling is employed on a simulated dataset so as to derive equations for calculating seismic velocities directly from porosity, effective pressure and cement volume values.
In this particular example, a nonlinear regression is performed on the simulated dataset for porosities ranging from 0.20 to 0.40. Strictly speaking, contact theory is only valid at relatively high porosities (of greater than approximately 0.20) and pressure sensitivity at lower porosities should be quantified using inclusion models (i.e. aspect ratios). Moreover, the applicants have found that regression easily becomes unstable if we include the whole porosity range since the shape of the velocity-porosity trends are very different at lower porosities relative to higher porosities.
In this example, the applicants chose a mathematical formulation similar to the one suggested by Eberhart-Phillips et al. (1989). However, slight modifications were made to obtain a satisfactory fit between the regression formulae and the simulated data.
Firstly, regression was performed on the “plane” representing well-cemented sandstones, with cement volumes ranging from 0.08-0.1. The resulting dry velocities were found to be given by equations (6) and (7) below as a function of porosity and effective pressure:
Vp
stiff=4992−10171·φ+9548·φ2+123·Peff0.2777 (6)
Vs
stiff=3013−5513·φ+3519·φ2+106·Peff0.2851 (7)
Next, regression on the “plane” representing unconsolidated sands (i.e. no cement) was performed. The simulated data here represents sands with only Walton smooth contacts (i.e. Hertz-Mindlin with a reduced shear factor Ft=0). The resulting dry velocities were found to be given by equations (8) and (9) below as a function of porosity and effective pressure:
Vp
soft=1204−6069·φ+5788·φ2+exp(7.377)·Peff0.1556 (8)
Vs
soft=769−3799·φ+3562·φ2+exp(6.839)·Peff0.1582 (9)
Finally, we define the dry velocities as a function of porosity, effective pressure and cement volume to be a weighted average of the soft and the stiff velocities, with respect to cement volume in accordance with equation (10):
where n(Peff)=3.5+2·Peff/ 20e6, reflecting that the weighting average will change with pressure. The same formulation is then applied to determine effective shear wave velocities also.
In accordance with the above, the applicants have established a heuristic approach to estimate fluid (using existing Gassman's theory) and pressure sensitivity in cemented sandstones using non-uniform contact theory combined with modified Hashin-Shtrikman elastic bounds. Embodiments of the invention expand on existing static models of cemented sandstones to account for stress sensitivity using elastic bounds in the porosity-moduli domain, where we define a soft bound to be stress sensitive (c.f., Hertz-Mindlin contact theory) and a stiff bound to be insensitive to stress (c.f., Dvorkin-Nur contact cement model). Based on the location of a data point (well log data or inverted seismic data) between these bounds, we are able to quantify expected pressure and fluid sensitivity in elastic and seismic properties (including moduli, velocities, acoustic impedance and Vp/Vs) of cemented sandstones. We have also established regression formulas that can be used to estimate dynamic rock physics templates for reservoir sandstones where the input parameters are confined to cement volume, porosity and effective pressure. This approach can be applied to predict the effect of pressure changes for example during 4-D monitoring analysis.
In the above embodiments, we assume that the cemented rock consists of a binary mixture of cemented and uncemented grain contacts, or “patchy cementation” (as shown in
A rough “fudge” parameter during application of contact theory is the so-called “slip-factor” or reduced shear factor. It is often seen that Hertz-Mindlin as well as Dvorkin-Nur overpredicts shear stiffness compared to measurements. In loose sands, it is shown that the Walton smooth contact theory (no friction) often gives the best fit. With increasing consolidation and/or pressure, the Hertz-Mindlin model gradually becomes more suited. Further studies will be conducted to better understand the physics behind the slip-factor and, if possible, obtain a better understanding of how this parameter varies as a function of, for example, burial depth and pressure.
Embodiments of the present method may assume clean, homogenous and isotropic reservoir sandstones. Presence of clay in the rock frame can be accounted for in the existing workflow. However, interbedded sand-shale sequences will affect pressure sensitivity in a more complex way. It is possible that depletion of reservoir sands will cause pore pressure increase in interbedded shales. Very little work has been done to model or document the effect of heterogeneity on pressure sensitivity and the plan is to investigate this in more detail.
It will be understood that various modifications may be made to the above described embodiments without departing from the scope of the present invention, as defined in the accompanying claims.
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/EP11/65887 | 9/13/2011 | WO | 00 | 5/8/2013 |
Number | Date | Country | |
---|---|---|---|
61382780 | Sep 2010 | US |