A program to create permeability fields that honor single-phase flow rate and pressure data

Abstract Accurate prediction (or simulation) of reservoir performance or contaminant transport in groundwater requires a realistic geological model representative of the reservoir/aquifer heterogeneity. Geostatistics provides tools for constructing such complex geological models constrained by different types of available (hard and soft) data and providing an assessment of related uncertainty. Permeability and flow data are nonlinearly related through the flow equations. Derivation of permeability models that honor flow response data is typically an inverse problem. This paper presents a FORTRAN program for generating permeability fields conditional to multiple-well single-phase flow rate and pressure data through an iterative inverse technique, called the sequential self-calibration (SSC) method. The SSC method is geostatistically-based, that is, it generates multiple equiprobable realizations that honor the input geostatistics of permeability and match pressure data for the given flow rate, under the given boundary conditions. The unique aspects of SSC are: (1) the master point concept that reduces the amount of computation, (2) a propagation mechanism based on kriging that accounts for spatial correlations of perturbations and (3) a fast method for computing all sensitivity coefficients within a single flow simulation run. Results from running the SSC code using an example data set are presented.

[1]  D. W. Peaceman.  Fundamentals of numerical reservoir simulation , 1977 .

[2]  Akhil Datta-Gupta,et al.  On the Sensitivity and Spatial Resolution of Transient Pressure and Tracer Data For Heterogeneity Characterization , 1997 .

[3]  A. Sahuquillo,et al.  Stochastic simulation of transmissivity fields conditional to both transmissivity and piezometric data—I. Theory , 1997 .

[4]  Clayton V. Deutsch,et al.  GSLIB: Geostatistical Software Library and User's Guide , 1993 .

[5]  Dean S. Oliver,et al.  Reparameterization Techniques for Generating Reservoir Descriptions Conditioned to Variograms and Well-Test Pressure Data , 1996 .

[6]  Dean S. Oliver,et al.  Multiple Realizations of the Permeability Field From Well Test Data , 1996 .

[7]  D. Oliver,et al.  Markov chain Monte Carlo methods for conditioning a permeability field to pressure data , 1997 .

[8]  Clayton V. Deutsch,et al.  ANNEALING TECHNIQUES APPLIED TO RESERVOIR MODELING AND THE INTEGRATION OF GEOLOGICAL AND ENGINEERING (WELL TEST) DATA , 1992 .

[9]  Clayton V. Deutsch,et al.  Challenges in reservoir forecasting , 1996 .

[10]  Andrés Sahuquillo,et al.  Stochastic simulation of transmissivity fields conditional to both transmissivity and piezometric data 2. Demonstration on a synthetic aquifer , 1997 .

[11]  J. Gómez-Hernández,et al.  Joint Sequential Simulation of MultiGaussian Fields , 1993 .

[12]  A. Journel.  Fundamentals of geostatistics in five lessons , 1991 .

[13]  A. S. Cullick,et al.  High Resolution Reservoir Models Integrating Multiple-Well Production Data , 1997 .

[14]  R. M. Srivastava,et al.  Integrating Seismic Data in Reservoir Modeling: The Collocated Cokriging Alternative , 1992 .

[15]  Roland N. Horne,et al.  A Procedure to Integrate Well Test Data, Reservoir Performance History and 4-D Seismic Information into a Reservoir Description , 1997 .

[16]  F. Roggero.  Direct Selection of Stochastic Model Realizations Constrained to Historical Data , 1997 .

[17]  M. Marietta,et al.  Pilot Point Methodology for Automated Calibration of an Ensemble of conditionally Simulated Transmissivity Fields: 1. Theory and Computational Experiments , 1995 .

[18]  Dean S. Oliver,et al.  Three-dimensional reservoir description from multiwell pressure data and prior information , 1997 .

[19]  J. J. Gómez-Hernández,et al.  Probabilistic assessment of travel times in groundwater modeling , 1994 .

[20]  Hua Zhu,et al.  Formatting and Integrating Soft Data: Stochastic Imaging via the Markov-Bayes Algorithm , 1993 .

[21]  Xian-Huan Wen,et al.  Significance of conditioning to piezometric head data for predictions of mass transport in groundwater modeling , 1996 .

[22]  A. Tarantola.  Inverse problem theory : methods for data fitting and model parameter estimation , 1987 .

[23]  Andrés Sahuquillo,et al.  Stochastic simulation of transmissivity fields conditional to both transmissivity and piezometric head data—3. Application to the Culebra formation at the Waste Isolation Pilot Plan (WIPP), New Mexico, USA , 1998 .

[24]  N. Sun.  Inverse problems in groundwater modeling , 1994 .