Probabilistic history matching using discrete Latin Hypercube sampling and nonparametric density estimation

Abstract This paper describes a new iterative procedure for probabilistic history matching using a discrete Latin Hypercube (DLHC) sampling method and nonparametric density estimation. The iterative procedure consists of selecting a set of models based on the history matching quality (normalized misfit) to generate histograms. The histograms are smoothed and used to estimate marginal probability densities of reservoir attributes. Three selection methods are evaluated. One of them is based on a global objective function (GOF) and the others are based on a local objective function (LOF), which is composed of influenced reservoir responses identified with the aim of a correlation matrix. The methodology was successfully applied to the UNISIM-I-H benchmark case, which is a reservoir model based on Namorado field, Campos basin, Brazil. Eight iterations with 450 combinations for each one were adequate to address the problem studied in this paper. To demonstrate the robustness of the proposed method and the consistency of the results, the iterative process was repeated 10 times and the discrepancy among the results was very small. The proposed method exhibited good convergence along the iterations, reducing the variability of the objective function (average normalized misfit) by approximately 90% from the first to the last iteration. Robustness, efficiency and facility of implementation are the key features of the proposed methodology.

[1]  Denis José Schiozer,et al.  UNISIM-I: Synthetic Model for Reservoir Development and Management Applications , 2015 .

[2]  Denis José Schiozer,et al.  A Methodology To Reduce Uncertainty Constrained to Observed Data , 2009 .

[3]  Takashi Goda,et al.  History matching with iterative Latin hypercube samplings and parameterization of reservoir heterogeneity , 2014 .

[4]  Denis José Schiozer,et al.  A new methodology to reduce uncertainties in reservoir simulation models using observed data and sampling techniques , 2010 .

[5]  Michael Andrew Christie,et al.  Population MCMC methods for history matching and uncertainty quantification , 2010, Computational Geosciences.

[6]  Tahir HusainT. Husain,et al.  Latin hypercube sampling for uncertainty analysis in multiphase modelling , 2008 .

[7]  Albert C. Reynolds,et al.  Combining the Ensemble Kalman Filter With Markov-Chain Monte Carlo for Improved History Matching and Uncertainty Characterization , 2012 .

[8]  J. Caers,et al.  History Matching With an Uncertain Geological Scenario , 2006 .

[9]  Dean S. Oliver,et al.  A Hybrid Markov Chain Monte Carlo Method for Generating Permeability Fields Conditioned to Multiwell Pressure Data and Prior Information , 1998 .

[10]  David Corne,et al.  An Adaptive Evolutionary Algorithm for History-Matching , 2013 .

[11]  Jonathan Carter,et al.  Using genetic algorithms for reservoir characterisation , 2001 .

[12]  Denis José Schiozer,et al.  A methodology to evaluate and reduce reservoir uncertainties using multivariate distribution , 2015 .

[13]  Hannes Kazianka,et al.  Bayesian spatial modeling and interpolation using copulas , 2011, Comput. Geosci..

[14]  B. Bobée,et al.  A distribution function based bandwidth selection method for kernel quantile estimation , 2001 .

[15]  P. A. Slotte,et al.  Response Surface Methodology Approach for History Matching and Uncertainty Assessment of Reservoir Simulation Models , 2008 .

[16]  Denis José Schiozer,et al.  A new procedure to reduce uncertainties in reservoir models using statistical inference and observed data , 2013 .

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

[18]  Albert C. Reynolds,et al.  Investigation of the sampling performance of ensemble-based methods with a simple reservoir model , 2013, Computational Geosciences.

[19]  Charles H. Langmuir,et al.  Calculation of phase equilibrium in mineral-melt systems , 1990 .

[20]  Á. Yustres,et al.  A review of Markov Chain Monte Carlo and information theory tools for inverse problems in subsurface flow , 2011, Computational Geosciences.

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

[22]  Alexandre A. Emerick.  Analysis of the performance of ensemble-based assimilation of production and seismic data , 2016 .

[23]  K. Osadetz,et al.  Using copulas for implementation of variable dependencies in petroleum resource assessment: Example from Beaufort-Mackenzie Basin, Canada , 2012 .

[24]  Naomi S. Altman,et al.  Bandwidth selection for kernel distribution function estimation , 1995 .

[25]  Arnaud Lange,et al.  Assisted history matching for the characterization of fractured reservoirs , 2009 .

[26]  Andrew R. Conn,et al.  Combining Adjoint Calculations and Quasi-Newton Methods for Automatic History Matching , 2006 .

[27]  T. LaForce,et al.  Bayesian Reservoir History Matching Considering Model and Parameter Uncertainties , 2012, Mathematical Geosciences.

[28]  Denis José Schiozer,et al.  Bayesian history matching using artificial neural network and Markov Chain Monte Carlo , 2014 .

[29]  D. Busby,et al.  Dynamic Data Assimilation by MCMC and Sequential Design of Experiments , 2008 .

[30]  Michael Andrew Christie,et al.  Comparison of Stochastic Sampling Algorithms for Uncertainty Quantification , 2010 .

[31]  Mike Christie,et al.  Uncertainty quantification for porous media flows , 2006, J. Comput. Phys..

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