Inversion of fractures with combination of production performance and in-situ stress analysis data

Abstract Successful identification of the fractures in fractured reservoirs is important to guarantee an effective development. Considering that the production performance of fractured reservoirs contains important information of the distribution of fractures, the production performance can be applied for the inversion of fractures. However, the inversion of fractures is difficult to achieve because it is an inverse problem with an inherent defect of multiplicity of solutions. In order to alleviate the defect, we estimate the possible distribution ranges of fractures which can be obtained by analyzing the correlation of fractures and in-situ stress, and the ranges are used as constraint conditions when the production performance is applied for obtaining the accurate distribution of fractures. Firstly, we choose the geometric parameters of fractures such as midpoint coordinate, azimuth and extension length of fractures as inversion parameters and use production data as inversion indexes. Secondly, we simulate the flow behavior in fractured reservoirs based on the Discrete Fracture Matrix (DFM) module of Matlab Reservoir Simulation Toolbox (MRST) to explicitly describe the effect of fractures on the flow behavior of fluid. Thirdly, the possible distribution ranges of fractures which can be obtained by analyzing the correlation between fractures and in-situ stress based on Griffith failure criterion are used as the constraint conditions of inversion parameters. Finally, Simultaneous Perturbation Stochastic Approximation (SPSA) algorithm is adopted to minimize the inversion objective function to obtain the accurate distribution of fractures. Theoretical cases verify that the method is effective for the accurate inversion of fractures while the inversion results of more fractures become worse because more fractures make the sensitivities of production performance to individual fractures decrease greatly.

[1]  Michael Andrew Christie,et al.  History Matching and Uncertainty Quantification: Multiobjective Particle Swarm Optimisation Approach , 2011 .

[2]  Jan M. Nordbotten,et al.  An efficient multi-point flux approximation method for Discrete Fracture-Matrix simulations , 2012, J. Comput. Phys..

[3]  Dean S. Oliver,et al.  Conditioning Geostatistical Models to Two-Phase Production Data , 1999 .

[4]  Mohan Kelkar,et al.  Efficient History Matching in Naturally Fractured Reservoirs , 2006 .

[5]  Rainer Helmig,et al.  Comparison of cell- and vertex-centered discretization methods for flow in a two-dimensional discrete-fracture–matrix system , 2009 .

[6]  Li Chun-yan,et al.  Upper Ordovician paleo tectonic stress field simulating and fracture distribution forecasting in Tazhong area of Tarim Basin,Xinjiang,China , 2011 .

[7]  Albert C. Reynolds,et al.  An Improved Implementation of the LBFGS Algorithm for Automatic History Matching , 2006 .

[8]  Stein Krogstad,et al.  Open-source MATLAB implementation of consistent discretisations on complex grids , 2012, Computational Geosciences.

[9]  Jose Rodrigues,et al.  Calculating derivatives for automatic history matching , 2006 .

[10]  Abbas Firoozabadi,et al.  Numerical Simulation of Water Injection in 2D Fractured Media Using Discrete-Fracture Model , 2001 .

[11]  G. I. Barenblatt,et al.  Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata] , 1960 .

[12]  Ossama Abdelkhalik,et al.  Dynamic Penalty Function Evolution Algorithms for History Matching of Oil and Gas Reservoir Models , 2012 .

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

[14]  Jiamin Jiang,et al.  Numerical study of complex fracture geometries for unconventional gas reservoirs using a discrete fracture-matrix model , 2015 .

[15]  A. Reynolds,et al.  Large-scale history matching with quadratic interpolation models , 2013, Computational Geosciences.

[16]  Hui Zhao,et al.  A two-stage efficient history matching procedure of non-Gaussian fields , 2016 .

[17]  N. Kalogerakis,et al.  Reliability of Horizontal Well Performance On a Field Scale Through Automatic History Matching , 1995 .

[18]  Li Yang,et al.  Theoretical research on production optimization of oil reservoirs , 2010 .

[19]  Hanqiao Jiang,et al.  The investigation of fracture aperture effect on shale gas transport using discrete fracture model , 2014 .

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

[21]  K. Aziz,et al.  Matrix-fracture transfer shape factors for dual-porosity simulators , 1995 .

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

[23]  Satomi Suzuki,et al.  History Matching of Naturally Fractured Reservoirs Using Elastic Stress Simulation and Probability Perturbation Method , 2007 .

[24]  Albert C. Reynolds,et al.  A Stochastic Optimization Algorithm for Automatic History Matching , 2007 .

[25]  Akhil Datta-Gupta,et al.  Stochastic Reservoir Modeling Using Simulated Annealing and Genetic Algorithm , 1995 .

[26]  Wang Bi-feng.  Quantitative relationship between crustal stress and parameters of tectonic fracture , 2010 .

[27]  J. Noorishad,et al.  An upstream finite element method for solution of transient transport equation in fractured porous media , 1982 .