A Pilot Point Guided Pattern Matching Approach to Integrate Dynamic Data into Geological Modeling

Abstract Methods based on multiple-point statistics (MPS) have been routinely used to characterize complex geological formations in the last decade. These methods use the available static data (for example, measured conductivities) for conditioning. Integrating dynamic data (for example, measured transient piezometric head data) into the same framework is challenging because of the complex non-linear relationship between the dynamic response and geology. The Ensemble PATtern (EnPAT) search method was recently developed as a promising technique to handle this problem. In this approach, a pattern is postulated to be composed of both parameter and state variables, and then, parameter values are sequentially (point-wise) simulated by directly sampling the matched pattern from an ensemble of training images of both geologic parameters and state variables. As a consequence, the updated ensemble of realizations of the geological parameters preserve curvilinear structures (i.e., non-multiGaussanity) as well as the complex relationship between static and dynamic data. Moreover, the uncertainty of flow and transport predictions can be assessed using the updated ensemble of geological models. In this work, we further modify the EnPAT method by introducing the pilot-point concept into the algorithm. More specifically, the parameter values at a set of randomly selected pilot point locations are simulated by the pattern searching procedure, and then a faster MPS method is used to complete the simulation by conditioning to the previously simulated pilot point values. This pilot point guided MPS implementation results in lower computational cost and more accurate inference of the parameter field. In addition, in some situations where there is sparsity of measured geologic static data, the EnPAT algorithm is extended to work only with the dynamic data. We employed a synthetic example to demonstrate the effectiveness of pilot points in the implementation of EnPAT, and also the capability of dynamic data to identify complex geologic structures when measured conductivity data are not available.

[1]  D. Oliver,et al.  Recent progress on reservoir history matching: a review , 2011 .

[2]  Xian-Huan Wen,et al.  Numerical modeling of macrodispersion in heterogeneous media: a comparison of multi-Gaussian and non-multi-Gaussian models , 1998 .

[3]  Charles F. Harvey,et al.  When good statistical models of aquifer heterogeneity go bad: A comparison of flow, dispersion, and mass transfer in connected and multivariate Gaussian hydraulic conductivity fields , 2003 .

[4]  J. Gómez-Hernández,et al.  To be or not to be multi-Gaussian? A reflection on stochastic hydrogeology , 1998 .

[5]  W. B. Whalley,et al.  The use of fractals and pseudofractals in the analysis of two-dimensional outlines: Review and further exploration , 1989 .

[6]  Jef Caers,et al.  Geostatistical History Matching Under Training-Image Based Geological Model Constraints , 2002 .

[7]  Wolfgang Nowak,et al.  Parameter Estimation by Ensemble Kalman Filters with Transformed Data , 2010 .

[8]  Andres Alcolea,et al.  Inverse problem in hydrogeology , 2005 .

[9]  Haiyan Zhou,et al.  A pattern‐search‐based inverse method , 2012 .

[10]  G. Mariéthoz,et al.  Bayesian inverse problem and optimization with iterative spatial resampling , 2010 .

[11]  L. Y. Hu,et al.  Multiple‐point geostatistics for modeling subsurface heterogeneity: A comprehensive review , 2008 .

[12]  Wen H. Chen,et al.  Generalization of the Ensemble Kalman Filter Using Kernels for Nongaussian Random Fields , 2009 .

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

[14]  Philippe Renard,et al.  Issues in characterizing heterogeneity and connectivity in non-multiGaussian media , 2008 .

[15]  Jianlin Fu,et al.  A Blocking Markov Chain Monte Carlo Method for Inverse Stochastic Hydrogeological Modeling , 2008 .

[16]  P. Kitanidis.  Quasi‐Linear Geostatistical Theory for Inversing , 1995 .

[17]  Behnam Jafarpour,et al.  A Probability Conditioning Method (PCM) for Nonlinear Flow Data Integration into Multipoint Statistical Facies Simulation , 2011 .

[18]  R. M. Srivastava,et al.  Multivariate Geostatistics: Beyond Bivariate Moments , 1993 .

[19]  H. H. Franssen,et al.  A comparison of seven methods for the inverse modelling of groundwater flow. Application to the characterisation of well catchments , 2009 .

[20]  Y. Rubin.  Applied Stochastic Hydrogeology , 2003 .

[21]  Liangping Li,et al.  An approach to handling non-Gaussianity of parameters and state variables in ensemble Kalman filtering , 2011 .

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

[23]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

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

[25]  Sebastien Strebelle,et al.  Conditional Simulation of Complex Geological Structures Using Multiple-Point Statistics , 2002 .

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

[27]  Philippe Renard,et al.  Connectivity metrics for subsurface flow and transport , 2013 .

[28]  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 .

[29]  M. Boucher,et al.  Interpretation of Interference Tests in a Well Field Using Geostatistical Techniques to Fit the Permeability Distribution in a Reservoir Model , 1984 .

[30]  Liangping Li,et al.  Inverse methods in hydrogeology: Evolution and recent trends , 2014 .

[31]  Y. Liu,et al.  Updating multipoint simulations using the ensemble Kalman filter , 2013, Comput. Geosci..

[32]  Sanjay Srinivasan,et al.  Efficient Conditional Simulation of Spatial Patterns Using a Pattern-Growth Algorithm , 2012 .

[33]  J. Jaime Gómez-Hernández,et al.  Pattern Recognition in a Bimodal Aquifer Using the Normal-Score Ensemble Kalman Filter , 2012, Mathematical Geosciences.

[34]  Dean S. Oliver,et al.  THE ENSEMBLE KALMAN FILTER IN RESERVOIR ENGINEERING-A REVIEW , 2009 .

[35]  Knud Skou Cordua,et al.  A Frequency Matching Method: Solving Inverse Problems by Use of Geologically Realistic Prior Information , 2012, Mathematical Geosciences.

[36]  Andrés Sahuquillo,et al.  Coupled inverse modelling of groundwater flow and mass transport and the worth of concentration data , 2003 .

[37]  Gregoire Mariethoz,et al.  The Direct Sampling method to perform multiple‐point geostatistical simulations , 2010 .

[38]  Xian-Huan Wen,et al.  Real-Time Reservoir Model Updating Using Ensemble Kalman Filter , 2005 .

[39]  Dean S. Oliver,et al.  The Ensemble Kalman Filter for Continuous Updating of Reservoir Simulation Models , 2006 .

[40]  A. Lavenue,et al.  Application of a coupled adjoint sensitivity and kriging approach to calibrate a groundwater flow model , 1992 .

[41]  Xian-Huan Wen,et al.  Simultaneous Integration of Pressure, Water Cut,1 and 4-D Seismic Data In Geostatistical Reservoir Modeling , 2006 .

[42]  Philippe Renard,et al.  A practical guide to performing multiple-point statistical simulations with the Direct Sampling algorithm , 2013, Comput. Geosci..

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

[44]  Harrie-Jan Hendricks Franssen,et al.  Modeling transient groundwater flow by coupling ensemble Kalman filtering and upscaling , 2012 .

[45]  J. Jaime Gómez-Hernández,et al.  Simultaneous Estimation of Geologic and Reservoir State Variables Within an Ensemble-Based Multiple-Point Statistic Framework , 2013, Mathematical Geosciences.

[46]  A. S. Cullick,et al.  A program to create permeability fields that honor single-phase flow rate and pressure data , 1999 .

[47]  Véronique Gervais,et al.  History matching of petroleum reservoir models by the Ensemble Kalman Filter and parameterization methods , 2013, Comput. Geosci..

[48]  A Blocking Markov Chain Monte Carlo Method for Inverse Stochastic Hydrogeological Modeling , 2009 .

[49]  S. P. Neuman,et al.  Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 1. Maximum Likelihood Method Incorporating Prior Information , 1986 .

[50]  A. J. Desbarats,et al.  Subsurface Flow and Transport: A Stochastic Approach , 1998 .

[51]  Philippe Renard,et al.  Blocking Moving Window algorithm: Conditioning multiple‐point simulations to hydrogeological data , 2010 .

[52]  D. A. Zimmerman,et al.  A comparison of seven geostatistically based inverse approaches to estimate transmissivities for modeling advective transport by groundwater flow , 1998 .

[53]  Yan Chen,et al.  Data assimilation for transient flow in geologic formations via ensemble Kalman filter , 2006 .

[54]  Wolfgang Nowak,et al.  Best unbiased ensemble linearization and the quasi‐linear Kalman ensemble generator , 2009 .

[55]  Eulogio Pardo-Igúzquiza,et al.  CONNEC3D: a computer program for connectivity analysis of 3D random set models☆ , 2003 .

[56]  A. Journel.  Geostatistics for Conditional Simulation of Ore Bodies , 1974 .

[57]  Graham E. Fogg,et al.  Groundwater Flow and Sand Body Interconnectedness in a Thick, Multiple-Aquifer System , 1986 .

[58]  Alexander Y. Sun,et al.  Sequential updating of multimodal hydrogeologic parameter fields using localization and clustering techniques , 2009 .

[59]  L. Hu.  Gradual Deformation and Iterative Calibration of Gaussian-Related Stochastic Models , 2000 .

[60]  András Bárdossy,et al.  Effects of non‐Gaussian copula‐based hydraulic conductivity fields on macrodispersion , 2012 .

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