Using differential evolution for compositional history-matching of a tight gas condensate well in the Montney Formation in western Canada

Abstract Production data analysis for low-permeability unconventional reservoirs is a challenging task, particularly for cases where multi-phase flow occurs within the reservoir. Analytical models developed to account for multi-phase flow typically require calculation of pseudo variables, which in turn require knowledge of relative permeability and fluid data. In the presence of sparse sampling, the analytical models often do not provide satisfactory results when there are so many unknown parameters. In such situations, numerical models are better suited, using a history matching framework to assist with reservoir and fluid characterization. In this work, we implement an assisted history-matching routine to characterize reservoir fluids and extract reservoir and hydraulic fracture properties for a hydraulically-fractured horizontal well completed in a tight gas condensate reservoir within the Montney Formation in western Alberta, Canada. The initial water distribution (e.g. movable water profile in the reservoir), in situ fluid (e.g. initial hydrocarbon composition with C 7+ properties) and reservoir properties (e.g. permeability in the matrix and around the fracture, and pressure dependent fracture permeability) are described in terms of 20 unknown parameters, which creates a high-dimensional inverse problem. We use the Differential Evolution algorithm, which is a powerful population-based optimization algorithm, and employ numerical compositional simulations to match pressure, water and hydrocarbon rates, and surface compositions of the produced fluids. Application of this optimization routine results in a good match to all measured data. The DE algorithm is repeated for an extra run to check for the existence of other non-unique solutions. The history-match results helped determine parameters for well/reservoir description and develop a compositional fluid model based on the measured separator composition data. The collected samples for both DE runs, along with one thousand extra samples from quasi-random sequence sampling design, provide a pool of data with invaluable information that are used to perform the global sensitivity analysis and to rank the contribution of each descriptive parameter on the variances of the reservoir outputs. In this way, the value of production data and surface compositions for the characterization of reservoir and fluid is quantified. This work aims to provide a practical and simple workflow for analysis of unconventional reservoirs where the direct analytical approaches cannot be applied.

[1]  Donald R. Jones,et al.  Efficient Global Optimization of Expensive Black-Box Functions , 1998, J. Glob. Optim..

[2]  Dan Simon,et al.  Evolutionary Optimization Algorithms , 2013 .

[3]  Jun Xia,et al.  An efficient integrated approach for global sensitivity analysis of hydrological model parameters , 2013, Environ. Model. Softw..

[4]  British Columbia Oil,et al.  Energy Briefing Note: The Ultimate Potential for Unconventional Petroleum from the Montney Formation of British Columbia and Alberta , 2013 .

[5]  Keith Atwood,et al.  An Integrated Approach for Understanding Oil and Gas Reserves Potential in Eagle Ford Shale Formation , 2011 .

[6]  Kendall E. Atkinson An introduction to numerical analysis , 1978 .

[7]  Pascal Longuemare Geomechanics in reservoir simulation , 2002 .

[8]  Michael J. Mayerhofer,et al.  What Is Stimulated Reservoir Volume , 2010 .

[9]  Yasin Hajizadeh,et al.  Sequential Bayesian Optimization Coupled with Differential Evolution for Geological Well Testing , 2014 .

[10]  Michael J. Mayerhofer,et al.  Integration of Microseismic-Fracture-Mapping Results With Numerical Fracture Network Production Modeling in the Barnett Shale , 2006 .

[11]  Volker Schwieger Sensitivity analysis as a general tool for model optimisation – examples for trajectory estimation , 2007 .

[12]  C. Cipolla,et al.  Reservoir Modeling in Shale-Gas Reservoirs , 2010 .

[13]  Shahab D. Mohaghegh,et al.  Fast Track Analysis of Shale Numerical Models , 2012 .

[14]  Lucas W. Bazan,et al.  Improving Production in the Eagle Ford Shale With Fracture Modeling, Increased Fracture Conductivity, and Optimized Stage and Cluster Spacing Along the Horizontal Wellbore , 2010, TGC 2010.

[15]  Anthony J. Culyer,et al.  Box-And-Whisker Plot , 2014, The SAGE Encyclopedia of Research Design.

[16]  Turgay Ertekin,et al.  Advanced Well Structures: An Artificial Intelligence Approach to Field Deployment and Performance Prediction , 2014 .

[17]  Louis Mattar,et al.  Practical Analytical Model To Simulate Production of Horizontal Wells With Branch Fractures , 2012 .

[18]  Michael Nikolaou,et al.  Analysis of Production History for Unconventional Gas Reservoirs With Statistical Methods , 2013 .

[19]  R. Storn,et al.  Differential Evolution - A simple and efficient adaptive scheme for global optimization over continuous spaces , 2004 .

[20]  Akhil Datta-Gupta,et al.  Improved Characterization and Performance Assessment of Shale Gas Wells by Integrating Stimulated Reservoir Volume and Production Data , 2011 .

[21]  I. Sobol On the distribution of points in a cube and the approximate evaluation of integrals , 1967 .

[22]  I. Jolliffe Principal Component Analysis , 2002 .

[23]  L. Klinkenberg The Permeability Of Porous Media To Liquids And Gases , 2012 .

[24]  Dongxiao Zhang,et al.  Ensemble Based Characterization and History Matching of Naturally Fractured Tight/Shale Gas Reservoirs , 2010 .

[25]  Kamy Sepehrnoori,et al.  Simulation of gas desorption and geomechanics effects for unconventional gas reservoirs , 2014 .

[26]  Matthew Mauldon,et al.  Determining Size and Shape of Fractures from Trace Data using a Differential Evolution Algorithm , 2006 .

[27]  Gary B. Lamont,et al.  Evolutionary Algorithms for Solving Multi-Objective Problems (Genetic and Evolutionary Computation) , 2006 .

[28]  Yasin Hajizadeh,et al.  Application of Differential Evolution as a New Method for Automatic History Matching , 2009 .

[29]  A. Makurat,et al.  The influence of natural fractures, faults and earth stresses on reservoir performance - geomechanical analysis by numerical modelling , 1994 .

[30]  T. Brubaker,et al.  Nonlinear Parameter Estimation , 1979 .

[31]  Curtis Hays Whitson,et al.  Modeling Water Flow in Hydraulically-Fractured Shale Wells , 2013 .

[32]  Donald L. Katz,et al.  Surface Tension of Hydrocarbons , 1939 .

[33]  Donald R. Jones,et al.  A Taxonomy of Global Optimization Methods Based on Response Surfaces , 2001, J. Glob. Optim..

[34]  Zhangxin Chen,et al.  Differential Evolution for Assisted History Matching Process: SAGD Case Study , 2013 .

[35]  Eamonn J. Keogh,et al.  Curse of Dimensionality , 2010, Encyclopedia of Machine Learning.

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

[37]  Abdollah Orangi,et al.  Unconventional Shale Oil and Gas-Condensate Reservoir Production, Impact of Rock, Fluid, and Hydraulic Fractures , 2011 .

[38]  Sherilyn Williams-Stroud,et al.  The use of microseismicity to understand subsurface-fracture systems and to increase the effectiveness of completions: Eagle Ford Shale, Texas , 2013 .

[39]  Aaron John White,et al.  Estimating Long Term Well Performance in the Montney Shale Gas Reservoir , 2014 .

[40]  Andrea Saltelli,et al.  Sensitivity Analysis for Importance Assessment , 2002, Risk analysis : an official publication of the Society for Risk Analysis.

[41]  Christopher R. Clarkson,et al.  Analysis of production data in shale gas reservoirs: Rigorous corrections for fluid and flow properties , 2012 .

[42]  Yasin Hajizadeh,et al.  Population-based sampling methods for geological well testing , 2015, Computational Geosciences.

[43]  Louis Mattar,et al.  Analyzing Production Data From Unconventional Gas Reservoirs With Linear Flow and Apparent Skin , 2010 .

[44]  Turgay Ertekin,et al.  Optimization Of The Design Of Transverse Hydraulic Fractures In Horizontal Wells Placed In Dual Porosity Tight Gas Reservoirs , 2011 .

[45]  Ted Urbancic,et al.  Microseismic Imaging of Hydraulic Fracture Complexity in the Barnett Shale , 2002 .

[46]  Antonin Settari,et al.  A Coupled Reservoir and Geomechanical Simulation System , 1998 .

[47]  James P. Brill,et al.  A Study of Two-Phase Flow in Inclined Pipes , 1973 .

[48]  Louis Mattar,et al.  Analytical Model for History Matching and Forecasting Production in Multifrac Composite Systems , 2012 .

[49]  T. Hastie,et al.  Using multivariate adaptive regression splines to predict the distributions of New Zealand ’ s freshwater diadromous fish , 2005 .

[50]  John R. Franchi,et al.  Integrated Flow Modeling , 2000 .

[51]  Elyezer Lolon,et al.  Modeling Well Performance in Shale-Gas Reservoirs , 2009 .

[52]  Aleksander Juell,et al.  Optimized Well Modeling of Liquid-Rich Shale Reservoirs , 2013 .

[53]  Riccardo Poli,et al.  Particle swarm optimization , 1995, Swarm Intelligence.

[54]  Christopher R. Clarkson,et al.  Production data analysis of tight gas condensate reservoirs , 2015 .

[55]  Christopher R. Clarkson,et al.  Production Data Analysis of Overpressured Liquid-Rich Shale Reservoirs: Effect of Degree of Undersaturation , 2014 .

[56]  J. Friedman Multivariate adaptive regression splines , 1990 .

[57]  Carlos A. Coello Coello,et al.  Multi-objective Optimization Using Differential Evolution: A Survey of the State-of-the-Art , 2008 .

[58]  F. Graziani,et al.  A Practical Global Sensitivity Analysis Methodology for Multi-Physics Applications , 2008 .

[59]  John R. Fanchi Shared earth modeling , 2002 .

[60]  Abeeb A. Awotunde,et al.  Efficient Drilling Time Optimization with Differential Evolution , 2014 .

[61]  James N. Albright,et al.  Acoustic Emissions as a Tool for Hydraulic Fracture Location: Experience at the Fenton Hill Hot Dry Rock Site , 1982 .

[62]  Luis F. Ayala,et al.  Analysis of Capillary Pressure and Relative Permeability Effects on the Productivity of Naturally Fractured Gas-Condensate Reservoirs , 2009 .

[63]  Max Kuhn,et al.  Applied Predictive Modeling , 2013 .

[64]  R. Storn,et al.  Differential Evolution: A Practical Approach to Global Optimization (Natural Computing Series) , 2005 .

[65]  Christopher R. Clarkson,et al.  A New Analytical Method for Analyzing Linear Flow in Tight/Shale Gas Reservoirs: Constant-Flowing-Pressure Boundary Condition , 2012 .

[66]  Wataru Sakamoto,et al.  MARS: selecting basis functions and knots with an empirical Bayes method , 2007, Comput. Stat..

[67]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .

[68]  Amos Nur,et al.  Pore pressure profiles in fractured and compliant rocks1 , 1994 .

[69]  John E. Walsh,et al.  Assessing the response of area burned to changing climate in western boreal North America using a Multivariate Adaptive Regression Splines (MARS) approach , 2009 .

[70]  Irving Langmuir THE CONSTITUTION AND FUNDAMENTAL PROPERTIES OF SOLIDS AND LIQUIDS. PART I. SOLIDS. , 1916 .

[71]  William John Lee,et al.  Gas Reserves Estimation in Resource Plays , 2010 .

[72]  Michael J. Mayerhofer,et al.  Imaging Seismic Deformation Induced by Hydraulic Fracture Complexity , 2009 .

[73]  Christopher A. Wright,et al.  Optimizing Horizontal Completion Techniques in the Barnett Shale Using Microseismic Fracture Mapping , 2004 .

[74]  Christopher A. Wright,et al.  Analysis and Prediction of Microseismicity Induced by Hydraulic Fracturing , 2004 .

[75]  Christopher K. I. Williams,et al.  Gaussian Processes for Machine Learning (Adaptive Computation and Machine Learning) , 2005 .

[76]  Melanie Mitchell,et al.  An introduction to genetic algorithms , 1996 .