History Matching for Determination of Fracture Permeability and Capillary Pressure
- Tao Gang (PetroTel) | Mohan G. Kelkar (U. of Tulsa)
- Document ID
- Society of Petroleum Engineers
- SPE Reservoir Evaluation & Engineering
- Publication Date
- October 2008
- Document Type
- Journal Paper
- 813 - 822
- 2008. Society of Petroleum Engineers
- 1.10.1 Drill string components and drilling tools (tubulars, jars, subs, stabilisers, reamers, etc), 5.8.6 Naturally Fractured Reservoir, 5.3.1 Flow in Porous Media, 5.6.4 Drillstem/Well Testing, 5.5.5 Evaluation of uncertainties, 5.1.5 Geologic Modeling, 5.3.2 Multiphase Flow, 3.3.2 Borehole Imaging and Wellbore Seismic, 5.5 Reservoir Simulation, 4.3.4 Scale, 5.5.8 History Matching, 5.1 Reservoir Characterisation, 5.4.2 Gas Injection Methods, 5.8.7 Carbonate Reservoir
- 0 in the last 30 days
- 916 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 5.00|
|SPE Non-Member Price:||USD 35.00|
For naturally fractured reservoirs, it is very difficult to quantify future prediction without proper fracture properties, including such factors as fracture permeability, fracture orientation, and fracture length, and capillary pressure. Capillary pressure and fracture permeability can interact in a way that can produce nonunique results. Evaluation of core data may provide us with some information about imbibition capillary pressure curves; however, quantification of fracture permeability is determined largely by matching production data.
This paper presents an integrated approach to history matching naturally fractured reservoirs by adjusting the fracture permeability of individual fractures and water/oil capillary pressure curves. This paper is focused mainly on history-matching procedure and not on geological-model construction. By minimizing an objective function to history matching production data, we generate estimates of fracture permeability of individual fractures and water/oil capillary pressure curves. All implementations were incorporated into a commercial simulator and iterated in the automatic-history-matching scheme. The adjoint method and an efficient direct solver were used to reduce central-processing-unit (CPU) time for calculating the sensitivity-coefficient matrix. A 2D synthetic case was used, with the fracture distribution from a Middle East reservoir, to validate this method.
Naturally fractured reservoirs have received much more attention in the last decade because of low oil recovery in many of them. Fluid flow in naturally fractured reservoirs is controlled largely by the distribution, orientation, and interconnectivity of the fracture system (Bourbiaux et al. 2005). A great deal of effort has been made in characterization of the geological fracture-network system based on the analysis, interpolation, and extrapolation of the fracture information acquired in wells, derived from seismic data, and obtained from outcrop analog data. As a result, there has been a large amount of work concerned with diagnosing and characterizing fractured reservoirs and developing viable techniques to use fractures effectively (Barenblatt et al. 1960; Bourbiaux et al. 1998, 2002; Dershowitz et al. 2000; Hu and Jenni 2005; Thomas et al. 1983, 1991; Gilman and Kazemi 1983). The reservoir characterization, modeling, and simulation of naturally fractured reservoirs present unique challenges compared with conventional reservoirs, because the fracture, the matrix, and the mass transfer between matrix and fracture have to be characterized properly. In this paper, production data were incorporated into the static model to characterize the fracture permeability of individual fractures, which is difficult to quantify properly without production data.
The capillary pressure in naturally fractured reservoirs plays a much more important role than that in conventional reservoirs. Capillary forces in fractured reservoirs are an important component of driving mechanisms, while the dynamic role of the capillary forces in a conventional reservoir is much more limited. In a fractured reservoir, capillary forces may contribute to the displacement process through imbibition, or they may oppose it through drainage (Reiss 1981).
It is essential to represent capillary pressure properly. Generally there are four different methods of measuring capillary pressure curves, including porous plate, centrifuge, air mercury injection, and water-vapor desorption. These methods are based on the availability of core data and suffer limitation because of scales over which the data are collected. It would be very useful if capillary pressure curves could be extracted from production data through the history-matching process.
Generally, geological models for naturally fractured reservoirs derived from static data alone cannot reproduce the field production history, which might be ascribed to the insufficient consideration of fracture effects on flow, insufficient dynamic characterization of the distribution of fracture systems, and insufficient consideration of the interaction between matrix and fracture. In other words, the hydrodynamic properties of a fracture-network system and capillary pressure curves need to be characterized properly. In this paper, production data were used to characterize fracture permeability and to estimate water/oil capillary pressure curves, which mainly control hydrodynamic properties of naturally fractured reservoirs.
Reconciling geologic models to the dynamic response of the reservoir is critical to building reliable reservoir models, which are very important in reservoir management and key economic decision making. Classical history-matching procedures whereby the reservoir parameters are adjusted manually by trial-and-error procedure can be tedious and time consuming. In recent years, several techniques have been developed for integrating production data to improve reservoir models, but only a few of these papers have been focused on naturally fractured reservoirs (Oliver 1994; Gao and Reynolds 2006; Cheng et al. 2004; Landa and Horne 1997; Cui and Kelkar 2005; Gang 2006; Gang and Kelkar 2006). The history matching of naturally fractured reservoirs requires that both the matrix and the fracture system be quantified properly. The field experience suggests that it is very difficult to quantify the fracture system properly without production data, while the matrix system might be characterized reasonably with static data. No research work so far has focused on estimating capillary pressure curves and fracture permeability by use of production data.
The history-matching process generally requires a optimization algorithm to find the minimum of an objective function that defines the mismatch between the observed and calculated data. Gradient-based algorithms are used widely to minimize the objective function in automatic history matching. For this algorithm, with an initial estimate of model parameters--in this case, the fracture permeability and water/oil capillary pressure parameters--the optimization algorithm will calculate the change of the model parameters with information from the computation of gradient and the Hessian of the objective function. The sensitivity-coefficient matrix, the partial derivatives of production data with respect to model parameters, is needed to evaluate the gradient and the Hessian.
In this paper we propose an approach to history matching fractured reservoirs through adjusting the fracture permeability of individual fractures and the capillary pressure parameters. First, with the initial estimate of model parameters, the adjoint method was used to calculate the sensitivity of production data with respect to model parameters, grid permeability, and capillary pressure. An efficient direct solver was used to reduce the CPU time, considering the special linear-algebra problem in this process. Second, the sensitivity of production data with respect to fracture permeability was computed through the chain rule. Third, changes of fracture permeability and capillary pressure parameters were computed by use of an optimization algorithm with information from the computation of the gradient and the Hessian. Fourth, the fracture permeability and capillary pressure curves were updated, and then the grid permeability was updated with new fracture permeability. Assuming that the relation between fracture permeability and grid permeability is known is a limitation of this paper. Though previous researchers have made assumptions regarding this relationship (da Silva 1989; Thomas et al. 1991; Cui and Kelkar 2005), we cannot confirm its validity independently. Fifth, the updated model was simulated to predict the reservoir response. The above process was repeated until satisfactory history-matching results were obtained. Fig. 1 shows the flow chart of this history-matching process. This approach avoids the time-consuming trial-and-error procedure associated with manual history matching. Through adjusting the fracture permeability instead of grid permeability directly, the geological consistency is preserved. The estimation of capillary pressure curves is very useful for quantifying reservoir performance when limited core data are available in naturally fractured reservoirs.
|File Size||2 MB||Number of Pages||10|
Barenblatt, G.I., Zheltov, I.P,. and Kochina, I.N. 1960. Basic Concepts in theTheory of Seepage of Homogeneous Liquids in Fissured Rocks. J. ofApplied Mathematics and Mechanics 24 (5): 1286-1303.DOI:10.1016/0021-8928(60)90107-6.
Bentsen, R.G. and Anli, J. 1977. Using Parameter Estimation Techniquesto Convert Centrifuge Data into a Capillary-Pressure Curve. SPEJ17 (1): 57-64; Trans., AIME, 263. SPE-5026-PA. DOI:10.2118/5026-PA.
Bourbiaux, B., Basquet, R., Cacas, M.C., Deniel, J.M., and Sarda, S. 2002.An Integrated Workflow to Accountfor Multi-Scale Fractures in Reservoir Simulation Models: Implementation andBenefits. Paper SPE 78489 presented at the Abu Dhabi InternationalPetroleum Exhibition and Conference, Abu Dhabi, UAE 13-16 October. DOI:10.2118/78489-MS.
Bourbiaux, B., Basquet, R., Daniel, J.M., Hu, L.Y., Jenni, S., Lange, A.,and Rasolofosaon, P. 2005. Fractured Reservoirs Modeling: A Review of theChallenges and Some Recent Solutions. First Break 23 (September):33-40.
Bourbiaux, B., Cacas, M.C., Sarda, S., and Sabathier, J.C. 1998. A Rapid andEfficient Methodology to Convert Fractured Reservoir Images Into aDual-Porosity Model. Revue de l'institut Français du Pétrole 53(6): 785-799.
Brooks, R.H. and Corey, A.T. 1964. Hydraulic Properties of Porous Media.Hydrology Paper No. 3, Colorado State University, Fort Collins, Colorado,22-27.
Brooks, R.H. and Corey, A.T. 1966. Properties of Porous Media AffectingFluid Flow. J. of Irrigation and Drainage Engineering 92 (IR2):61-68.
Cheng, H., Wen, X.-H., Milliken. W.J., Datta-Gupta, A. 2004. Field Experiences with Assisted andAutomatic History Matching Using Streamline Models. Paper SPE 89857presented at the SPE Annual Technical Conference and Exhibition, Houston, 26-29September. DOI 10.2118/89857-MS.
Cui, H. and Kelkar, M. 2005. Automatic History Matching ofNaturally Fractured Reservoirs and a Case Study. Paper SPE 94037 presentedat the SPE Western Regional Meeting, Irvine, California, USA, 30 March-1 April.DOI: 10.2118/94037-MS.
da Silva, F.V. 1989. Primaryand Enhanced Recovery of Ekofisk Field: A Single- and Double-Porosity NumericalSimulation Study. Paper SPE 19840 presented at the SPE Annual TechnicalConference and Exhibition, San Antonio, Texas, USA, 8-11 Oct. DOI:10.2118/19840-MS.
Dershowitz, B., LaPointe, P. Eiben, T., and Wei, L. 2000. Integration of Discrete FeatureNetwork Methods With Conventional Simulator Approaches. SPEREE3 (2): 165-170. SPE-62498-PA. DOI: 10.2118/62498-PA.
ECLIPSE Reference Manual 2004A. 2004. Houston:Schlumberger/GeoQuest.
Fletcher, R. 1987. Practical Methods of Optimization. New York City:John Wiley & Sons.
Gang, T. 2006. Data Integration and Model Improvement for NaturallyFractured Reservoirs. PhD dissertation, University of Tulsa, Tulsa,Oklahoma.
Gang, T. and Kelkar, M. 2006. Efficient History Matching inNaturally Fractured Reservoirs. Paper SPE 99578 presented at the SPE/DOESymposium on Improved Oil Recovery, Tulsa, 22-26 April. DOI:10.2118/99578-MS.
Gang, T. and Kelkar, M. 2007. A More General Capillary PressureCurve and Its Estimation from Production Data. Paper SPE 108180 presentedat the Rocky Mountain Oil & Gas Technology Symposium, Denver, 16-18 April.DOI: 10.2118/108180-MS.
Gao, G. and Reynolds, A.C. 2006. An Improved Implementation of theLBFGS Algorithm for Automatic History Matching. SPEJ 11 (1),5-17. SPE-90058-PA. DOI: 10.2118/90058-PA.
Gilman, J.R. and Kazemi, H. 1983. Improvements in Simulation ofNaturally Fractured Reservoirs. SPEJ 23 (4): 695-707.SPE-10511-PA. DOI: 10.2118/10511-PA.
Hu, L.Y. and Jenni, S. 2005. History Matching of Object-BasedStochastic Reservoir Models. SPEJ 10 (3): 312-323.SPE-81503-PA. DOI: 10.2118/81503-PA.
Huang, D.D., Honarpour, M.M., and Al-Hussainy, R. 1997. An Improved Modelfor Relative Permeability and Capillary Pressure Incorporating Wettability.Paper SCA 9718 presented at the Society of Core Analysts InternationalSymposium, Calgary, 7-10 September.
Landa, J.L. and Horne, R.N. 1997. A Procedure to Integrate Well TestData, Reservoir Performance History and 4-D Seismic Information into aReservoir Description. Paper SPE 38653 presented at the SPE AnnualTechnical Conference and Exhibition, San Antonio, Texas, 5-8 October. DOI:10.2118/38653-MS.
Li, K. 2004. GeneralizedCapillary Pressure and Relative Permeability Model Inferred from FractalCharacterization of Porous Media. Paper SPE 89874 presented at the SPEAnnual Technical Conference and Exhibition, Houston, 26-29 September. DOI:10.2118/89874-MS.
Oliver, D.S. 1994. Incorporation of Transient Pressure Data Into ReservoirCharacterization. In Situ 18 (3): 243-275.
Parker, J.C., Lenhard, R.J., and Kuppusamy, T. 1987. A Parametric Model forConstitutive Properties Governing Multiphase Flow in Porous Media. WaterResources Research 23 (4): 618-624. DOI:10.1029/WR023i004p00618.
Pooladi-Darvish, M. and Firoozabadi, A. 2000. Cocurrent and CountercurrentImbibition in a Water-Wet Matrix Block. SPEJ 5 (1): 3-11.SPE-38443-PA. DOI: 10.2118/38443-PA.
Reiss, L.H. 1981. The Reservoir Engineering Aspects of FracturedFormations. Houston: Gulf Publishing Co.
Skjaeveland, S.M., Siqveland, L.M., Kjosavik, A., Hammervold Thomas, W.L.,and Virnovsky, G.A. 2000. Capillary Pressure Correlation forMixed-Wet Reservoirs. SPEREE 3 (1): 60-67. SPE-60900-PA. DOI:10.2118/60900-PA.
Thomas, L.K., Dixon, T.N., and Pierson, R. 1991. Ekofisk Nitrogen Injection.SPEFE 6 (2): 151-160. SPE-19839-PA. DOI: 10.2118/19839-PA.
Thomas, L.K., Dixon, T.N., and Pierson, R.G. 1983. Fractured Reservoir Simulation.SPEJ 23 (1): 42-54. SPE-9305-PA. DOI: 10.2118/9305-PA.