Compositional Modeling of Discrete-Fractured Media Without Transfer Functions by the Discontinuous Galerkin and Mixed Methods
- Hussein Hoteit (Reservoir Engr. Research Inst.) | Abbas Firoozabadi (Reservoir Engr. Research Inst.)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- September 2006
- Document Type
- Journal Paper
- 341 - 352
- 2006. Society of Petroleum Engineers
- 5.4.2 Gas Injection Methods, 5.3.1 Flow in Porous Media, 4.1.2 Separation and Treating, 1.2.3 Rock properties, 5.5.1 Simulator Development, 5.8.8 Gas-condensate reservoirs, 5.2.1 Phase Behavior and PVT Measurements, 5.5 Reservoir Simulation, 5.8.6 Naturally Fractured Reservoir, 5.3.2 Multiphase Flow, 4.3.4 Scale
- 1 in the last 30 days
- 744 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 12.00|
|SPE Non-Member Price:||USD 35.00|
In a recent work, we introduced a numerical approach that combines the mixed-finite-element (MFE) and the discontinuous Galerkin (DG) methods for compositional modeling in homogeneous and heterogeneous porous media. In this work, we extend our numerical approach to 2D fractured media. We use the discrete-fracture model (crossflow equilibrium) to approximate the two-phase flow with mass transfer in fractured media. The discrete-fracture model is numerically superior to the single-porosity model and overcomes limitations of the dual-porosity model including the use of a shape factor. The MFE method is used to solve the pressure equation where the concept of total velocity is invoked. The DG method associated with a slope limiter is used to approximate the species-balance equations. The cell-based finite-volume schemes that are adapted to a discrete-fracture model have deficiency in computing the fracture/fracture fluxes across three and higher intersecting-fracture branches. In our work, the problem is solved definitively because of the MFE formulation. Several numerical examples in fractured media are presented to demonstrate the superiority of our approach to the classical finite-difference method.
Compositional modeling in fractured media has broad applications in CO2, nitrogen, and hydrocarbon-gas injection, and recycling in gas condensate reservoirs. In addition to species transfer, the compressibility effects should be also considered for such applications. Heterogeneities and fractures add complexity to the fluid-flow modeling. Several conceptually different models have been proposed in the literature for the simulation of flow and transport in fractured porous media.
The single-porosity approach uses an explicit computational representation for fractures (Ghorayeb and Firoozabadi 2000; Rivière et al. 2000). It allows the geological parameters to vary sharply between the matrix and the fractures. However, the high contrast and different length scales in the matrix and fractures make the approach unpractical because of the ill conditionality of the matrix appearing in the numerical computations (Ghorayeb and Firoozabadi 2000).The small control volumes in the fracture grids also add a severe restriction on the timestep size because of the Courant-Freidricks-Levy (CFL) condition if an explicit temporal scheme is used.
|File Size||4 MB||Number of Pages||12|
Ács, G., Doleschall, S., and Farkas, É. 1985. General Purpose CompositionalModel. SPEJ 25 (4): 543-553. SPE-10515-PA.
Arbogast, T., Douglas, J., and Hornung, U. 1990. Derivation of the doubleporosity model of single phase via homogenization theory. SIAM J.Math. Anal. 21 (4): 823-836.
Baca, R., Arnett, R., and Langford, D. 1984. Modeling fluid flow infractured porous rock masses by finite element techniques. Intl. J. Num.Meth. Fluids 4: 337-348.
Bastian, P., Chen, Z., Ewing, R.E., Helmig, R., Jakobs, H., andReichenberger, V. 2000. Numerical solution of multiphase flow in fracturedporous media. InNumerical Treatment of Multiphase Flows inPorous Media. Z. Chen, R.E. Ewing, and Z.C. Shi (eds.). Berlin:Springer-Verlag.
Bourgeat, A. 1984. Homogenized behavior of diphasic flow in naturallyfissured reservoir with uniform fractures. Comp. Methods in AppliedMechanics and Engineering 47: 205-217.
Brezzi, F. and Fortin, M. 1991. Mixed and Hybrid Finite ElementMethod. New York: Springer-Verlag.
Chavent, G. and Cockburn, B. 1989. The local projection P0P1 -discontinuous Galerkin finite element method for scalarconservation laws. M2AN 23 (4): 565-592.
Chavent, G. and Jaffré, J. 1986. Mathematical Models and Finite Elementsfor Reservoir Simulation: Studies in Mathematics and its Applications.North Holland, Amsterdam.
Chavent, G. and Roberts, J-E. 1991. A unified physical presentation ofmixed, mixed-hybrid finite element method and standard finite differenceapproximations for the determination of velocities in water flow problems.Adv. Water Resour. 14 (6): 329-348.
Chen, Z. and Ewing, R. 1997a. Comparison of various formulations of thethree-phase flow in porous media. J. Comp. Phys. 132:362-373.
Chen, Z. and Ewing, R. 1997b. From single-phase tocompositional flow: applicability of mixed finite elements. Transport inPorous Media 27: 225-242.
Cockburn, B. and Shu, C. 1989. TVB Runge Kutta local projectiondiscontinuous Galerkin finite element method for conservative laws II: Generalframe-work. Math. Comp. 52: 411-435.
Cockburn, B. and Shu, C. 1998. The Runge-Kutta Discontinuous Galerkin Methodfor Conservative Laws V: Multidimentional Systems. J. Comput. Phys.141: 199-224.
Cockburn, B. and Shu, C. 2001. Runge-Kutta discontinuous Galerkin method forconvection-dominated problems. J. Scientific Computing 16 (3):173-261.
Durlofsky, L.J. and Chien, M.C.H. 1993. Development of a MixedFinite-Element-Based Compositional Reservoir Simulator. Paper SPE 25253presented at the SPE Symposium on Reservoir Simulation, New Orleans, 28February-3 March.
Ewing, R.E. and Heinemann, R.F. 1983. Incorporation of Mixed Finite ElementMethods in Compositional Simulation for Reduction of Numerical Dispersion.Paper SPE 12267 presented at the SPE Reservoir Simulation Symposium, SanFrancisco, 15-18 November.
Ghorayeb, K. and Firoozabadi, A. 2000. Numerical Study of Natural Convectionand Diffusion in Fractured Porous Media. SPEJ 5 (1): 12-20.SPE-51347-PA.
Granet, S., Fabrie, P., Lemmonier, P., and Quitard, M. 1998. A single phasesimulation of fractured reservoir using a discrete representation of fractures.Paper presented at the European Conference on the Mathematics of Oil Recovery,Peebles, Scotland, U.K., 8-11 September.
Granet, S., Fabrie, P., Lemonnier, P., and Quitard, M. 2001. A two-phaseflow simulation of a fractured reservoir using a new fissure element method.J. Petroleum Science and Engineering 32 (18): 35-52.
Hoteit, H. and Firoozabadi, A. 2005. Multicomponent fluid flow bydiscontinuous Galerkin and mixed methods in unfractured and fracturedmedia. Water Resour. Res. 41 (11): W11412,doi:10.1029/2005WR004339.
Hoteit, H. and Firoozabadi, A. 2006. Compositional Modeling by theCombined Discontinuous Galerkin and Mixed Methods. SPEJ 11(1): 19-34. SPE-90276-PA.
Hoteit, H., Ackerer, P., Mosé, R., Erhel, J., and Philippe, B. 2004. NewTwo-Dimensional Slope Limiters for Discontinuous Galerkin Methods on ArbitraryMeshes.Int. J. Numer. Meth. Eng. 61 (14): 2566-2593.
Karimi-Fard, M. and Firoozabadi, A. 2003. Numerical Simulation of WaterInjection in Fractured Media Using the Discrete Fractured Model and theGalerkin Method. SPEREE 6 (2): 117-126. SPE-71615-PA.
Karimi-Fard, M., Durlofsky, L.J., and Aziz, K. 2004. An Efficient Discrete-Fracture ModelApplicable for General-Purpose Reservoir Simulators. SPEJ 9(2): 227-236. SPE-79699-PA.
Kazemi, H. 1969. PressureTransient Analysis of Naturally Fractured Reservoirs With Uniform FractureDistribution. SPEJ 9 (12): 451-462; Trans., AIME,246. SPE-2156-PA.
Kim, J.G. and Deo, M.D. 1999. Comparison of the Performance of aDiscrete Fracture Multiphase Model With Those Using Conventional Methods.Paper SPE 51928 presented at the SPE Reservoir Simulation Symposium, Houston,14-17 February.
Kim, J.G. and Deo, M.D. 2000. Finite element discretefracture model for multiphase flow in porous media," AIChE J.46 (6): 1120-1130.
Mallison, B., Gerritsen, M., Jessen, K., and Orr F.M. 2005. High-Order Upwind Schemes forTwo-Phase, Multicomponent Flow. SPEJ 10 (3): 291-311.SPE-79691-PA.
Monteagudo, J. and Firoozabadi, A. 2004. Control-volume method fornumerical simulation of two-phase immiscible flow in 2D and 3Ddiscrete-fracture media. Water Resour. Res. 7: 1-20.doi:10.1029/2003WR002996.
Noorishad, J. and Mehran, M. 1982. An upstream finite element method for thesolution of transient transport equation in fractured porous media. WaterResour. Res. 18 (3):588-596.
Peddibhotla, S., Datta-Gupta, A., and Xue, G. 1997. Multiphase Streamline Modeling inThree Dimensions: Further Generalizations and a Field Application. PaperSPE 38003 presented at the SPE Reservoir Simulation Symposium, Dallas, 8-11June.
Putti, M., Yeh, W., and Mulder, W. 1990. A triangular finite volume approachwith high-resolution upwind terms for the solution of groundwater transportequations. Water Resour. Res. 26 (12): 2865-2880.
Raviart, P. and Thomas, J. 1977. A mixed hybrid finite element method forthe second order elliptic problem. Lectures Notes in Mathematics 606.New York: Springer-Verlag, 292-315.
Rivière, B., Wheeler, M.F., and Banas, K. 2000. Part II. Discontinuous Galerkin MethodApplied to Single Phase Flow in Porous Media. ComputationalGeosciences (2000) 4 (4): 337-341.
Rubin, B. and Blunt, M.J. 1991. Higher-Order Implicit Flux-LimitingSchemes for Black Oil Simulation. Paper SPE 21222 presented at the SPESymposium on Reservoir Simulation, Anaheim, California, 17-20 February.
Tan, C. and Firoozabadi, A. 1995. Theoretical analysis of miscibledisplacement in fractured porous media: I. Theory. J. Can. Petrol.Technol. 34 (2): 17-27.
Thomas, L.K., Dixon, T.N., and Pierson, R.G. 1983. Fractured Reservoir Simulation.SPEJ 23 (1): 42-54. SPE-9305-PA.
Warren, J.E. and Root, P.J. 1963. The Behavior of Naturally FracturedReservoirs. SPEJ 3 (11): 245-255; Trans., AIME,228. SPE-426-PA.
Watts, J. 1986. A CompositionalFormulation to the Pressure and Saturation Equations. SPERE 1(3): 243-252; Trans., AIME, 281. SPE-12244-PA.