History Matching of Naturally Fractured Reservoirs Using Elastic Stress Simulation and Probability Perturbation Method
- Satomi Suzuki (Stanford University) | Colin Daly (Schlumberger) | Jef Karel Caers (Stanford University) | Dietmar Mueller (Shell Intl.)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- March 2007
- Document Type
- Journal Paper
- 118 - 129
- 2007. Society of Petroleum Engineers
- 1.6.9 Coring, Fishing, 5.5.8 History Matching, 5.1 Reservoir Characterisation, 5.1.5 Geologic Modeling, 4.1.2 Separation and Treating, 4.3.4 Scale, 5.6.9 Production Forecasting, 5.8.6 Naturally Fractured Reservoir, 4.1.5 Processing Equipment, 5.6.4 Drillstem/Well Testing, 1.2.2 Geomechanics, 5.5 Reservoir Simulation, 5.1.2 Faults and Fracture Characterisation, 1.10 Drilling Equipment, 5.1.7 Seismic Processing and Interpretation
- 2 in the last 30 days
- 754 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
The application of elastic stress simulation for fracture modeling provides a more realistic description of fracture distribution than conventional statistical and geostatistical techniques, allowing the integration of geomechanical data and models into reservoir characterization. The geomechanical prediction of the fracture distribution accounts for the propagation of fracture caused by stress perturbation associated with faults. However, the challenge lies in estimating the past remote stress conditions which induced structural deformation and fracturing, the limited applicability of the elasticity assumption, and the latent uncertainty in the structural geometry of faults. The integration of historical production data and well-test permeability into geomechanical fracture modeling is a practical way to reduce such uncertainty. We propose to combine geostatistical algorithms for history matching with geomechanical elastic simulation models for developing an integrated yet efficient fracture modeling tool.
This paper presents an integrated approach to history matching of naturally fractured reservoirs which includes (1) fracture trend prediction through elastic stress simulation; (2) geostatistical population of fracture density based on a fracture trend model; (3) fracture permeability modeling integrating fracture density, matrix permeability and well-test permeability; and (4) numerical flow simulation and history matching. All of these implementations are incorporated into a single forward modeling process and iterated in the automatic history-matching scheme. To obtain a history match on a reservoir model, we jointly perturb the large-scale fracture trend and local-scale geostatistical fluctuations of fracture densities rather than perturbing permeability calibrated from fractures. This strategy enables us to preserve the geological/geomechanical consistency throughout the history-matching process. The geomechanically simulated fracture trend model is calibrated to both production data and the reservoir geological structure (faults and horizons) by searching for the optimum remote stress condition for elastic stress-field simulation. The latter is achieved by matching the actually observed structural deformation with the simulated one. The smaller-scale fluctuation of fracture density is simultaneously history matched through the probability perturbation method of Caers (Caers 2003; Hoffman and Caers 2005; Caers 2007). The methodology is presented on a synthetic reservoir application.
The modeling of the density and pattern of fracture distributions can take different approaches depending on the origin and the type of fracture sets and on the ultimate reservoir engineering questions raised. In this paper, we focus on the modeling of shear fractures which are generated by structural deformation accompanied with fault slip. Recently, an application of the elastic stress simulation has been proposed for predicting the pattern of shear/tensile fractures or the pattern of secondary faults and shown promising results (Bourne and Willemse 2001; Maerten et al. 2002; Bourne et al. 2001). The elastic simulation numerically simulates the structural deformation of the reservoir by solving linear elasticity equations under given boundary conditions, and simultaneously calculates the corresponding stress/strain tensor fields (Bourne and Willemse 2001; Maerten et al. 2002; Bourne et al. 2001; Daly and Mueller 2004; Roxar FracPerm Reference Manual 2005). The boundary conditions consist of (1) location/geometry of fault surface, (2) stress conditions or displacement conditions on the fault surfaces, and (3) the remote loads applied to the structure at the time of structural deformation accompanied with fault slippage. First, satisfying boundary conditions and by minimizing strain energy, the linear elastic equations are solved to obtain a structural deformation field which is expressed by the displacement vector. Next, strain field is computed from the displacement gradient based on the definition of strain. Finally, under the assumption of elasticity, stress is calculated from strain by means of Hook's law.
|File Size||2 MB||Number of Pages||12|
Bernabe, Y. 1995. The Transport Properties of Networks of Cracks and Pores.J. Geophysical Research 100 (b3): 4231-4241.DOI: http://dx.doi.org/10.1029/94JB02986.
Bourne S.J., Rijkels, L., Stephenson, B.J., and Willemse, E.J.M. 2001.Predictive Modeling of Naturally Fractured Reservoirs Using Geomechanics andFlow Simulation. GeoArabia6 (1): 27-42.
Bourne, S.J. and Willemse, E.J.M. 2001. Elastic Stress Control on thePattern of Tensile Fracturing Around a Small Fault Network at Nash Point, UK.J. of Structural Geology 23: 1753-1770.
Caers, J. 2003. HistoryMatching Under Training-Image-Based Geological Constraints. SPEJ 8(3): 218-226. SPE-74716-PA. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=74716-PA.
Caers, J. 2007. Comparing the Gradual Deformation with the ProbabilityPerturbation for Solving Inverse Problem. Mathematical Geology 39(2).
Crouch, S.L. and Starfield, A.M. 1983. Boundary Element Methods in SolidMechanics. London: Allen and Unwin.
Daly, C. and Mueller, D. 2004. Characterization and Modeling of FracturedReservoirs: Static Model. Paper presented at the European Conference on theMathematics of Oil Recovery, Cannes, France, 29 August-2 September.
Heffer, K.J., King, P.R., and Jones, A.D.W. 1999. Fracture Modeling as Part ofIntegrated Reservoir Description. Paper SPE 53347 presented at the SPEMiddle East Oil Show and Conference, Bahrain, 20-23 February. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=53347-MS.
Hoffman, B. T. and Caers, J. 2005. Regional Probability Perturbations forHistory Matching. J. of Petroleum Science and Engineering 46:53-71.
Hu, L-Y. and Le Ravalec-Dupin, M. 2004. An Improved Gradual DeformationMethod for Reconciling Random and Gradient Searches in StochasticOptimizations. Mathematical Geology 36 (6): 703-720.
Hu, L-Y., Blanc, G., and Noetinger, B. 2001. Gradual Deformation andIterative Calibration of Sequential Simulations. Mathematical Geology33: 475-489. DOI: http://dx.doi.org/10.1023/A:1011088913233.
Jaeger, J.C. and Cook, N.G.W. 1979. Fundamentals of Rock Mechanics.London: Chapman and Hall.
Journel, A.G. 1994. Modeling Uncertainty: Some Conceptual Thoughts. InGeostatistics for the Next Century, Dimitrakoponlos R et.al.Dordrecht, The Netherlands: Kluwer. 30-43.
King, G. and Cocco, M. 2001. Fault Interaction by Elastic Stress Changes;New Clues from Earthquake Sequences. Advances in Geophysics 44:1-46.
Maerten, L., Gillespie, P., and Pollard, D. 2002. Effects of Local StressPerturbation on Secondary Fault Development. J. of Structural Geology24: 145-153.
Press, W.H., Flannery, B.P., Teukolsky, S.A., Vetterling, W.T. 1992.Numerical Recipes in C. 2nd edition. Cambridge U. Press.
Roxar. 2005. FracPerm Reference Manual.
Zoback, M.D. and Townend, J. 2001. Implications of Hydrostatic PorePressures and High Crustal Strength for the Deformation of IntraplateLithosphere. Tectonophysics 336 (1-4): 19-30. DOI: http://dx.doi.org/10.1016/S0040-1951(01)00091-9.