Abstract
Flow of fluids and transport of solutes in porous media are subjects of wide interest in several fields of applications: reservoir engineering, subsurface hydrology, chemical engineering, etc. In this paper we will study two-phase flow in a model consisting of two different types of sediments. Here, the absolute permeability, the relative permeabilities and the capillary pressure are discontinuous functions in space. This leads to interior boundary value problems at the interface between the sediments. The saturation Sw will be discontinuous or experience large gradients at the interface. A new solution procedure for such problems will be presented. The method combines the modified method of characteristics with a weak formulation where the basis functions are discontinuous at the interior boundary. The modified method of characteristics will provide a good first approximation for the jump in the discontinuous basis functions, which leads to a fast converging iterative solution scheme for the complete problem.
The method has been implemented in a two-dimensional simulator, and results from numerical experiments will be presented.
Similar content being viewed by others
References
O. Axelsson and V.A. Barker, Finite Element Solution of Boundary Value Problems(Academic Press, London, 1984).
J.W. Barrett and K.W. Morton, Approximate symmetrization and Petrov-Galerkin methods for diffusion-convection problems, Comput. Methods Appl. Mech. Engrg. 45 (1984) 97-122.
H. Bertin, A. Graue and T. Eilertsen, Two dimensional imaging of viscous and capillary effects in a dipping stratified heterogeneous media, in: International Symposium of the Society of Core Analysts, Stavanger, Norway (1994).
Ø Bøe, Finite-element methods for porous media flow, Ph.D. thesis, Department of Applied Mathematics, University of Bergen (1990).
K.A. Børresen and A. Graue, Numerical simulations and 2d-imaging experiments of flow in tilted cross beds, in: Proceedings SCA Symposium(September 1996).
K.A. Børresen and A. Graue, Comparing simulations to 2d-imaging experiments of fluid flow in a dipping crossbedded reservoir model, in: 71th Annual Technical Conference and Exhibition of SPE, SPE paper No. 36731 (October 1996).
K.A. Børresen and A. Graue, Fluid flow across capillary heterogeneities, experiments and simulations, submitted to Transport in Porous Media.
F. Bratvedt, K. Bratvedt, C.F. Buchholz, L. Holden, H. Holden and N.H. Risebro, A new front-tracking method for reservoir simulation, SPE Reservoir Engineering (1992) pp. 107-116.
T. Bu and L.B. Håøy, On the importance of correct inclusion of capillary pressure in reservoir simulation, in: Proceedings 8th European IOR Symposium(1995).
G. Chavent and J. Jaffre, Mathematical Models and Finite Elements for Reservoir Simulation(North-Holland, Amsterdam, 1987).
G. Chavent, J. Jaffre and J.E. Roberts, Generalized cell-centered finite volume methods: application to two-phase flow in porous media, in: Computational Science for the 21st Century, eds. M.O. Bristau et al. (Wiley, New York, 1997) pp. 231-241.
H.K. Dahle, Adaptive characteristic operator splitting techniques for convection-dominated diffusion problems in one and two space dimensions, Ph.D. thesis, Department of Mathematics, University of Bergen (1988).
H.K. Dahle, M.S. Espedal and O. Sævareid, Characteristic, local grid refinement techniques for reservoir flow problems, Internat. J. Numer. Methods Engrg. 34 (1992) 1051-1069.
J.Jr. Douglas and T.F. Russell, Numerical methods for convection-dominated diffusion problems based on combining the method of characteristics with finite element or finite difference procedures, SIAM J. Numer. Anal. 19 (1982) 871-885.
C.J. van Duijn, J. Molenaar and M.J. de Neef, The effect of the capillary forces on immiscible two-phase flow in heterogeneous porous media, Transport in Porous Media 21 (1995) 71-93.
T. Eilertsen, K.A. Børresen, H. Bertin and A. Graue, Dimensional imaging experiments and numerical simulation of cross layered reservoir models, submitted to J. Petroleum Science Engrg.
B.G. Ersland, On numerical methods for including the effect of capillary presure forces on two-phase, immiscible flow in a layered porous media, Ph.D. thesis, Department of Applied Mathematics, University of Bergen (1996).
B.G. Ersland and M. Espedal, Domain decomposition method for heterogeneous reservoir flow, in: Domain Decomposition Methods in Science and Engineering, eds. R. Glowinsky et al. (Wiley, New York, 1997) pp. 477-484.
M.S. Espedal and R.E. Ewing, Characteristic Petrov-Galerkin subdomain methods for two-phase immiscible flow, Comput. Methods Appl. Mech. Engrg. 64 (1987) 113-135.
M. Espedal, K.J. Hersvik and B.G. Ersland, Domain decomposition method for heterogeneous reservoir flow, in: Proceedings 10th International Conference on Domain Decomposition Methods, Boulder, USA (1997).
T. Gimse and N.H. Risebro, Solution of the cauchy problem for a conservation law with a discontinuous flux function, SIAM J. Math. Anal. 23 (1992) 635-648.
T. Gimse and N.H. Risebro, A note on reservoir simulation for heterogeneous porous media, Transport in Porous Media 10 (1993) 257-270.
A. Graue, Imaging the effects of capillary heterogeneities on local saturation development in long-core floods, SPE paper No. 21101, Society of Petroleum Engineers (1990).
R. Hansen and M.S. Espedal, On the numerical solution of nonlinear flow models with gravity, Internat. J. Numer. Methods Engrg. 38 (1995) 2017-2032.
P.W. Hemker, A numerical study of stiff two-point boundary problems, Ph.D. thesis, Mathematisch Centrum, Amsterdam (1977).
J. Jaffre, Flux calculation at the interface between two rock types for two-phase flow in porous media, Transport in Porous Media 21 (1995) 195-207.
K.H. Karlsen, K. Brusdal, H.K. Dahle, S. Evje and K.A. Lie, The corrected operator splitting approach applied to nonlinear advection-diffusion problem, to appear in Comput. Methods Appl. Mech. Engrg.
K.H. Karlsen, K.-A. Lie, N.H. Risebro and J. Fryen, A front-tracking approach to a two-phase fluid-flow model with capillary forces, In Situ 22(1) (special issue on reservoir simulation) (1998) 59-89.
D.S. Koetsier, Mechanisms affecting capillary entrapment in heterogeneous porous media, Ph.D. thesis, Reservoir Engineering, Mining and Petroleum Engineering, Delft University of Technology (1993).
T.F.M. Kortekaas, Water/oil displacement characteristics in crossbedded reservoir zones, Society of Petroleum Engineers J. (1985) pp. 917-926.
J.I. Kristiansen, M. Mikkelsen and K. Esbensen, A modified Leverett approach and pls-regression for integrated formation evaluation, in: Proceedings 33rd Annual Symposium of Society of Professional Well Log Analysis(1992).
K.W. Morton, Numerical solution of convection-diffusion problems, Appl. Math. Math. Comput. 12 (Chapman & Hall, 1996).
R. Nybø, To-fase reservoarmodelller med diskontinuitet i absolutt og relativ permeabilitet, Thesis for the cand. scient. degree, Department of Applied Mathematics, University of Bergen (1997).
O. Sævareid, On local grid refinement techniques for reservoir flow problems, Ph.D. thesis, Department of Applied Mathematics, University of Bergen (1990).
X.-C. Tai, T.O.W. Johansen, H.K. Dahle and M. Espedal, A characteristic domain splitting method, in: Domain Decomposition Methods in Science and Engineering, eds. R. Glowinsky et al. (Wiley, New York, 1997) pp. 317-324.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Ersland, B.G., Espedal, M.S. & Nybø, R. Numerical methods for flow in a porous medium with internal boundaries. Computational Geosciences 2, 217–240 (1998). https://doi.org/10.1023/A:1011554320427
Issue Date:
DOI: https://doi.org/10.1023/A:1011554320427