Multirate-Transfer Dual-Porosity Modeling of Gravity Drainage and Imbibition
- Ginevra Di Donato (Imperial College) | Huiyun Lu (Imperial College) | Zohreh Tavassoli (Imperial College) | Martin Julian Blunt (Imperial College)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- March 2007
- Document Type
- Journal Paper
- 77 - 88
- 2007. Society of Petroleum Engineers
- 5.5.3 Scaling Methods, 5.5 Reservoir Simulation, 4.1.5 Processing Equipment, 5.4.1 Waterflooding, 4.3.4 Scale, 5.8.6 Naturally Fractured Reservoir, 5.1.2 Faults and Fracture Characterisation, 5.1 Reservoir Characterisation, 5.5.7 Streamline Simulation, 5.4.2 Gas Injection Methods, 1.6.9 Coring, Fishing, 5.3.4 Reduction of Residual Oil Saturation, 5.2.1 Phase Behavior and PVT Measurements, 4.1.2 Separation and Treating
- 3 in the last 30 days
- 997 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
We develop a physically motivated approach to modeling displacement processes in fractured reservoirs. To find matrix/fracture transfer functions in a dual-porosity model, we use analytical expressions for the average recovery as a function of time for gas gravity drainage and countercurrent imbibition. For capillary-controlled displacement, the recovery tends to its ultimate value with an approximately exponential decay (Barenblatt et al. 1990). When gravity dominates, the approach to ultimate recovery is slower and varies as a power law with time (Hagoort 1980). We apply transfer functions based on these expressions for core-scale recovery in field-scale simulation.
To account for heterogeneity in wettability, matrix permeability, and fracture geometry within a single gridblock, we propose a multirate model (Ponting 2004). We allow the matrix to be composed of a series of separate domains in communication with different fracture sets with different rate constants in the transfer function.
We use this methodology to simulate recovery in a Chinese oil field to assess the efficiency of different injection processes. We use a streamline-based formulation that elegantly allows the transfer between fracture and matrix to be accommodated as source terms in the 1D transport equations along streamlines that capture the flow in the fractures (Di Donato et al. 2003; Di Donato and Blunt 2004; Huang et al. 2004).
This approach contrasts with the current Darcy-like formulation for fracture/matrix transfer based on a shape factor (Gilman and Kazemi 1983) that may not give the correct average behavior (Di Donato et al. 2003; Di Donato and Blunt 2004; Huang et al. 2004). Furthermore, we show that recovery is exceptionally sensitive to parameters that describe the physics of the displacement process, highlighting the need to make careful core-scale measurements of recovery.
Di Donato et al.(2003) and Di Donato and Blunt (2004) proposed a dual-porosity streamline-based model for simulating flow in fractured reservoirs. Conceptually, the reservoir is composed of two domains: a flowing region with high permeability that represents the fracture network and a stagnant region with low permeability that represents the matrix (Barenblatt et al. 1960; Warren and Root 1963). The streamlines capture flow in the flowing regions, while transfer from fracture to matrix is accommodated as source/sink terms in the transport equations along streamlines. Di Donato et al. (2003) applied this methodology to study capillary-controlled transfer between fracture and matrix and demonstrated that using streamlines allowed multimillion-cell models to be run using standard computing resources. They showed that the run time could be orders of magnitude smaller than equivalent conventional grid-based simulation (Huang et al. 2004). This streamline approach has been applied by other authors (Al-Huthali and Datta-Gupta 2004) who have extended the method to include gravitational effects, gas displacement, and dual-permeability simulation, where there is also flow in the matrix. Thiele et al. (2004) have described a commercial implementation of a streamline dual-porosity model based on the work of Di Donato et al. that efficiently solves the 1D transport equations along streamlines.
|File Size||1 MB||Number of Pages||12|
Al-Huthali, A. and Datta-Gupta, A. 2004. Streamline simulation ofcounter-current imbibition in naturally fractured reservoirs. J. ofPetroleum Science and Engineering 43 (3-4): 271-300.
Aronofsky, J.S., Masse, L., and Natanson, S.G. 1958. A Model for the Mechanism of OilRecovery From the Porous Matrix Due to Water Invasion in FracturedReservoirs. Trans., AIME 213. SPE-932-G. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=932-G.
Barenblatt, G.I., Zheltov, I.P., and Kochina, I.N. 1960. Basic concepts inthe theory of seepage of homogeneous liquids in fissured rocks. J. Appl.Math. Mech., Engl. Transl. 24: 1286-1303.
Barenblatt, G.I., Entov, V.M., and Ryzhik, V.M. 1990. Theory of FluidFlows Through Natural Rocks. Dordrecht: Kluwer Academic Publishers.
Batycky, R.P., Blunt, M.J., and Thiele, M.R. 1997. A 3D Field-Scale Streamline-BasedReservoir Simulator. SPERE 12 (4): 246-254. SPE-36726-PA.DOI: http://www.spe.org/elibrary/servlet/spepreview?id=36726-PA.
Behbahani, H. and Blunt, M.J. 2005. Analysis of Imbibition in Mixed-WetRocks Using Pore-Scale Modeling. SPEJ 10 (4): 466-474.SPE-90132-PA. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=90132-PA.
Belayneh, M. and Cosgrove, J.W. 2004. Fracture pattern variations around amajor fold and their implications regarding fracture prediction using limiteddata: an example from the Bristol Channel Basin, UK. In The initiation,Propagation and Arrest of Joints and Other Fractures, ed.J.W.Cosgrove and T.Engelder. London: Geological Soc., Special Publications231: 89-102.
Civan, F. and Rasmussen, M.L. 2001. Asymptotic Analytical Solutions forImbibition Waterfloods in Fractured Reservoirs. SPEJ 6(2):171-181. SPE-71312-PA. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=71312-PA.
Civan, F., Wang, W., and Gupta, A. 1999. Effect of Wettability and Matrix toFracture Transfer on the Waterflooding in Fractured Reservoirs. Paper SPE52197 presented at the SPE Mid-Continent Operations Symposium, Oklahoma City,Oklahoma, 28-31 March. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=52197-MS.
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, 28 August-2 September.
Di Donato, G. and Blunt, M.J. 2004. Streamline-based dual-porositysimulation of reactive transport and flow in fractured reservoirs. WaterResources Research 40 (4). DOI: http://dx.doi.org/10.1029/2003WR002772.
Di Donato, G., Huang, W., and Blunt, M.J. 2003. Streamline-Based Dual-PorositySimulation of Fractured Reservoirs. Paper SPE 84036 presented at the SPEAnnual Technical Conference and Exhibition, Denver, Colorado, 5-8 October. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=84036-MS.
Gelhar, L.W. and Collins, M.A. 1971. General analysis of longitudinaldispersion in nonuniform flow. Water Resources Research 7 (6):1511-1521.
Gilman J.R. and Kazemi, H. 1983. Improvements in Simulation ofNaturally Fractured Reservoirs. SPEJ 23 (5): 695-707.SPE-10511-PA. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=10511-PA.
Hagoort, J. 1980. Oil Recoveryby Gravity Drainage. SPEJ 20 (3): 139-150. SPE-7424-PA. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=7424-PA.
Horie, T., Firoozabadi, A., and Ishimoto, K. 1990. Laboratory Studies of CapillaryInteraction in Fracture/Matrix Systems. SPERE 5 (3): 353-360.SPE-18282-PA. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=18282-PA.
Huang, W., Di Donato, G., and Blunt, M.J. 2004. Comparison ofstreamline-based and grid-based dual porosity simulation. J. of PetroleumScience and Engineering 43: 129-137.
Kazemi, H., Gilman, J.R., and Elsharkawy, A.A. 1992. Analytical and Numerical Solution ofOil Recovery From Fractured Reservoirs With Empirical Transfer Functions.SPERE 7 (2):219-227. SPE-19849-PA. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=19849-PA.
Matthäi, S.K., Mezentsev, A., and Belayneh, M. 2005. Control-Volume Finite-ElementTwo-Phase Flow Experiments With Fractured Rock Represented by Unstructured 3DHybrid (CVFE) Meshes. Paper SPE 93341 presented at the SPE ReservoirSimulation Symposium, Houston, 31 January-2 February. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=93341-MS.
Moreno, J., Kazemi, H., and Gilman, J.R. 2004. Streamline Simulation ofCountercurrent Water-Oil and Gas-Oil Flow in Naturally Fractured Dual-PorosityReservoirs. Paper SPE 89880 presented at the SPE Annual TechnicalConference and Exhibition, Houston, 26-29 September. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=89880-MS.
Morrow, M.R and Mason, G. 2001. Recovery of oil by spontaneous imbibition.Current Opinion in Colloid and Interface Science 6: 321-337.
Ponting, D. 2004. Characterization and Modeling of Fractured Reservoirs:Flow Simulation. Paper presented at the European Conference on the Mathematicsof Oil Recovery, Cannes, France, 28 August-2 September.
Por, G.J., Boerrigter, P., Maas, J.G., and de Vries, A. 1989. A Fractured Reservoir SimulatorCapable of Modeling Block-Block Interaction. Paper SPE 19807 presented atthe SPE Annual Technical Conference and Exhibition, San Antonio, Texas, 8-11October. DOI: http://www.spe.org/elibrary/servlet/spepreview?id=19807-MS.
Tavassoli, Z., Zimmerman, R.W., and Blunt, M.J. 2005a. Analytic Analysis forOil Recovery During Counter-Current Imbibition in Strongly Water-Wet Systems.Transport in Porous Media 58: 173-189. DOI: http://dx.doi.org/10.1007/s11242-004-5474-4.
Tavassoli, Z., Zimmerman, R.W., and Blunt, M.J. 2005b. Analysis ofcounter-current imbibition with gravity in weakly water-wet systems. J. ofPetroleum Science and Engineering 48: 94-104.
Terez, I.E. and Firoozabadi, A. 1999. Water Injection in Water-WetFractured Porous Media: Experiments and a New Model With ModifiedBuckley-Leverett Theory. SPEJ (1999) 4(2): 135-141.SPE-56854-PA. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=56854-PA.
Thiele, M.R., Batycky R.P., Iding, M., and Blunt, M.J. 2004. Extension ofStreamline-Based Dual Porosity Flow Simulation to Realistic Geology. Paperpresented at the European Conference on the Mathematics of Oil Recovery,Cannes, France, 28 August-2 September.
Warren, J.E. and Root, P.J 1963. The Behavior of Naturally FracturedReservoirs. SPEJ 3: 245-255; Trans., AIME,228. SPE-426-PA. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=426-PA.
Zhang, X., Morrow, N.R., and Ma, S. 1996. Experimental Verification of aModified Scaling Group for Spontaneous Imbibition. SPERE 11(4): 280-285. SPE-30762-PA. DOI:http://www.spe.org/elibrary/servlet/spepreview?id=30762-PA.
Zhou, D., Jia, L., Kamath, J., and Kovscek, A.R. 2002. Scaling ofcounter-current imbibition processes in low-permeability porous media. J. ofPetroleum Science and Engineering 33: 61-74.