Conditional Statistical Moment Equations for Dynamic Data Integration in Heterogeneous Reservoirs
- Liyong Li (Chevron Corp.) | Hamdi A. Tchelepi (Stanford University)
- Document ID
- Society of Petroleum Engineers
- SPE Reservoir Evaluation & Engineering
- Publication Date
- June 2006
- Document Type
- Journal Paper
- 280 - 288
- 2006. Society of Petroleum Engineers
- 7.6.2 Data Integration, 4.3.4 Scale, 5.1 Reservoir Characterisation, 4.5 Offshore Facilities and Subsea Systems, 5.3.2 Multiphase Flow, 5.5 Reservoir Simulation, 5.1.8 Seismic Modelling, 5.1.5 Geologic Modeling, 5.3.1 Flow in Porous Media
- 0 in the last 30 days
- 316 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
An inversion method for the integration of dynamic (pressure) data directly into statistical moment equations (SMEs) is presented. The method is demonstrated for incompressible flow in heterogeneous reservoirs. In addition to information about the mean, variance, and correlation structure of the permeability, few permeability measurements are assumed available. Moreover, few measurements of the dependent variable are available. The first two statistical moments of the dependent variable (pressure) are conditioned on all available information directly. An iterative inversion scheme is used to integrate the pressure data into the conditional statistical moment equations (CSMEs). That is, the available information is used to condition, or improve the estimates of, the first two moments of permeability, pressure, and velocity directly. This is different from Monte Carlo (MC) -based geostatistical inversion techniques, where conditioning on dynamic data is performed for one realization of the permeability field at a time. In the MC approach, estimates of the prediction uncertainty are obtained from statistical post-processing of a large number of inversions, one per realization.
Several examples of flow in heterogeneous domains in a quarter-five-spot setting are used to demonstrate the CSME-based method. We found that as the number of pressure measurements increases, the conditional mean pressure becomes more spatially variable, while the conditional pressure variance gets smaller. Iteration of the CSME inversion loop is necessary only when the number of pressure measurements is large. Use of the CSME simulator to assess the value of information in terms of its impact on prediction uncertainty is also presented.
The properties of natural geologic formations (e.g., permeability) rarely display uniformity or smoothness. Instead, they usually show significant variability and complex patterns of correlation. The detailed spatial distributions of reservoir properties, such as permeability, are needed to make performance predictions using numerical reservoir simulation. Unfortunately, only limited data are available for the construction of these detailed reservoir-description models. Consequently, our incomplete knowledge (uncertainty) about the property distributions in these highly complex natural geologic systems means that significant uncertainty accompanies predictions of reservoir flow performance.
To deal with the problem of characterizing reservoir properties that exhibit such variability and complexity of spatial correlation patterns when only limited data are available, a probabilistic framework is commonly used. In this framework, the reservoir properties (e.g., permeability) are assumed to be a random space function. As a result, flow-related properties such as pressure, velocity, and saturations are random functions. We assume that the available information about the permeability field includes a few measurements in addition to the spatial correlation structure, which we take here as the two-point covariance. This incomplete knowledge (uncertainty) about the detailed spatial distribution of permeability is the only source of uncertainty in our problem. Uncertainty about the detailed distribution of the permeability field in the reservoir leads to uncertainty in the computed predictions of the flow field (e.g., pressure).
|File Size||2 MB||Number of Pages||9|
Dagan, G. 1985. Stochastic Modeling of Groundwater Flow by Unconditional andConditional Probabilities, the Inverse Problem. Water Resour. Res. 21 (1):65-72.
Dietrich, C.R. and Newsam, G.N. 1989. A Stability Analysis of theGeostatistical Approach to Aquifer Transmissivity Identification.Stochastic Hydro. Hydraul. 3: 293-316.
Gomez-Hernandez, J.J., Sahuquillo, A., and Capilla, J.E. 1997. StochasticSimulation of Transmissivity Fields Conditional to Both Transmissivity andPiezometric Data—I. Theory. J. of Hydrology 203: 162-174.
Hernandez, A.F., Neuman, S.P., Guadagnini, A., and Carrera, J. 2003. Conditioning Mean SteadyState Flow on Hydraulic Head and Conductivity Through GeostatisticalInversion. Stochastic Envir. Res. and Risk Assessment 17: 329-338.
Kitanidis, P. and Vomvoris, E. 1983. A Geostatistical Approach to TheInverse Problem in Groundwater Modelling (Steady State) and One-DimensionalSimulations. Water Resour. Res. 19 (3): 677-690.
Li, L. and Tchelepi, H.A. 2003. Conditional Stochastic MomentEquations for Uncertainty Analysis of Flow in Heterogeneous Reservoirs.SPEJ 8 (4): 392-400. SPE-87337-PA.
Li, L. and Tchelepi, H. 2004. Statistical Moment Equations for Flow inComposite Heterogeneous Porous Media. Paper presented at the 9th EuropeanConference on the Mathematics of Oil Recovery, Cannes, France, 30 August-2September.
Li, L., Tchelepi, H., and Zhang, D. 2003. Perturbation-BasedMoment Equation Approach for Flow in Heterogeneous Porous Media: ApplicabilityRange and Analysis of High-Order Terms. J. of Computational Physics 188:296-317.
Oliver, D.S., Cunha, L., and Reynolds, A.C. 1997. Markov Chain Monte CarloMethods for Conditioning a Permeability Field to Pressure Data. MathematicalGeology 29: 61-91.
Sun, N.-Z. and Yeh, W.W.-G. 1992. A Stochastic Inverse Solution forTransient Groundwater Flow: Parameter Identification and Reliability Analysis.Water Resour. Res. 28 (12): 3269-3280.
Wen, X.-H., Deutsch, C.V., and Cullick, A.S. 2003. Inversion of DynamicProduction Data for Permeability: Fast Streamline-Based Computation ofSensitivity Coefficients of Fractional Flow Rate. J. of Hydrology 281:296-312.
Xue, G. and Datta-Gupta, A. 1997. Structure Preserving Inversion: AnEffective Approach to Conditioning Stochastic Reservoir Models to DynamicData. Paper SPE 38727 presented at the SPE Annual Technical Conference andExhibition, San Antonio, Texas, 5-8 October.
Yeh, T.-C.J., Jin, M., and Hanna, S. 1996. An Iterative Stochastic InverseMethod: Conditional Effective Transmissivity and Hydraulic Head Fields. WaterResour. Res. 32 (1): 85-92.