History matching with iterative Latin hypercube samplings and parameterization of reservoir heterogeneity

Abstract History matching can be formulated as a global minimization of the difference between time-series observations and numerical results. Existence of a number of unknown parameters, however, makes the dimensionality of history matching intractably high. This study addresses two issues involved in solving history matching with a feasible number of simulation runs. One is the computational effort required for searching an optimal solution, the other the ill-posedness owing to reservoir heterogeneity. A new population-based search algorithm named iterative Latin hypercube samplings is proposed for the former and we would show the superior convergence of our proposed algorithm over those of other famous population-based search algorithms for a broad class of functions. As for the latter, parameterization of reservoir heterogeneity using orthonormal basis functions is considered, which can significantly reduce the number of unknown parameters to be optimized. Numerical example would reveal that our approach of history matching is efficient and of practical use.

[1]  P. N. Suganthan,et al.  Differential Evolution: A Survey of the State-of-the-Art , 2011, IEEE Transactions on Evolutionary Computation.

[2]  Thomas J. Santner,et al.  The Design and Analysis of Computer Experiments , 2003, Springer Series in Statistics.

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

[4]  K. Pruess,et al.  TOUGH2 User's Guide Version 2 , 1999 .

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

[6]  Runze Li,et al.  Design and Modeling for Computer Experiments , 2005 .

[7]  A. Owen A Central Limit Theorem for Latin Hypercube Sampling , 1992 .

[8]  M. Stein Large sample properties of simulations using latin hypercube sampling , 1987 .

[9]  Hans O. Jahns,et al.  A Rapid Method for Obtaining a Two-Dimensional Reservoir Description From Well Pressure Response Data , 1966 .

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

[11]  R. Schulze-Riegert,et al.  Evolutionary Algorithms Applied to History Matching of Complex Reservoirs , 2002 .

[12]  P. N. Suganthan,et al.  Differential Evolution Algorithm With Strategy Adaptation for Global Numerical Optimization , 2009, IEEE Transactions on Evolutionary Computation.

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

[14]  C. W. Harper,et al.  A FORTRAN IV program for comparing ranking algorithms in quantitative biostratigraphy , 1984 .

[15]  J. Caers,et al.  The Probability Perturbation Method: A New Look at Bayesian Inverse Modeling , 2006 .

[16]  L. Durlofsky,et al.  Efficient real-time reservoir management using adjoint-based optimal control and model updating , 2006 .

[17]  P. Jacquard,et al.  Permeability Distribution From Field Pressure Data , 1965 .

[18]  James Kennedy,et al.  Particle swarm optimization , 2002, Proceedings of ICNN'95 - International Conference on Neural Networks.

[19]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[20]  Andries Petrus Engelbrecht,et al.  A study of particle swarm optimization particle trajectories , 2006, Inf. Sci..

[21]  J. Walsh A Closed Set of Normal Orthogonal Functions , 1923 .

[22]  Ingrid Daubechies,et al.  Ten Lectures on Wavelets , 1992 .

[23]  D. E. Goldberg,et al.  Genetic Algorithms in Search, Optimization & Machine Learning , 1989 .

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

[25]  A. Owen Monte Carlo Variance of Scrambled Net Quadrature , 1997 .

[26]  N. Ahmed,et al.  Discrete Cosine Transform , 1996 .

[27]  Martin Cohn Walsh Functions, Sequency, and Gray Codes , 1971 .

[28]  Roland N. Horne,et al.  A Multiresolution Approach to Reservoir Parameter Estimation Using Wavelet Analysis , 2000 .

[29]  George Kingsley Zipf,et al.  Human behavior and the principle of least effort , 1949 .

[30]  Steven A. Orszag,et al.  CBMS-NSF REGIONAL CONFERENCE SERIES IN APPLIED MATHEMATICS , 1978 .

[31]  Michael Andrew Christie,et al.  Prediction under uncertainty in reservoir modeling , 2002 .

[32]  M. D. McKay,et al.  A comparison of three methods for selecting values of input variables in the analysis of output from a computer code , 2000 .

[33]  B. Jafarpour,et al.  History matching with an ensemble Kalman filter and discrete cosine parameterization , 2008 .

[34]  Rainer Storn,et al.  Differential Evolution – A Simple and Efficient Heuristic for global Optimization over Continuous Spaces , 1997, J. Glob. Optim..