Multidimensional Upwinding for Multiphase Transport in Porous Media
- Jeremy E. Kozdon (Stanford University) | Bradley Mallison (Chevron) | Margot Gerritsen (Stanford University) | Wen H. Chen (Chevron)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- June 2011
- Document Type
- Journal Paper
- 263 - 272
- 2011. Society of Petroleum Engineers
- 5.4.1 Waterflooding, 5.5 Reservoir Simulation
- multidimensional upstream methods, multiphase flow, monotone, volume relaxed methods, grid orientation effect
- 1 in the last 30 days
- 409 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 12.00|
|SPE Non-Member Price:||USD 35.00|
Multidimensional transport for reservoir simulation is typically solved by applying 1D numerical methods in each spatial-coordinate direction. This approach is simple, but the disadvantage is that numerical errors become highly correlated with the underlying computational grid. In many real-field applications, this can result in strong sensitivity to grid design not only for the computed saturation/composition fields but also for critical integrated data such as breakthrough times. Therefore, to increase robustness of simulators, especially for adverse-mobility-ratio flows that arise in a variety of enhanced-oil-recovery (EOR) processes, it is of much interest to design truly multidimensional schemes for transport that remove, or at least strongly reduce, the sensitivity to grid design.
We present a new upstream-biased truly multidimensional family of schemes for multiphase transport capable of handling countercurrent flow arising from gravity. The proposed family of schemes has four attractive properties: applicability within a variety of simulation formulations with varying levels of implicitness, extensibility to general grid topologies, compatibility with any finite-volume flow discretization, and provable stability (monotonicity) for multiphase transport. The family is sufficiently expressive to include several previously developed multidimensional schemes, such as the narrow scheme, in a manner appropriate for general-purpose reservoir simulation.
A number of waterflooding problems in homogeneous and heterogeneous media demonstrate the robustness of the method as well as reduced transverse (cross-wind) diffusion and grid-orientation effects.
|File Size||640 KB||Number of Pages||10|
Aavatsmark, I. 2002. AnIntroduction to Multipoint Flux Approximations for Quadrilateral Grids.Computational Geosciences 6 (3-4): 405-432. doi:10.1023/A:1021291114475.
Ács, G., Doleschall, S., and Farkas, E. 1985. General Purpose CompositionalModel. SPE J. 25 (4): 543-553. SPE-10515-PA. doi:10.2118/10515-PA.
Aziz, K. and Settari, A. 1979. Petroleum Reservoir Simulation.Dordrecht, The Netherlands: Kluwer Academic Publishers.
Colella, P. 1990. Multidimensional UpwindMethods for Hyperbolic Conservation Laws. Journal of ComputationalPhysics 87 (1): 171-200. doi:10.1016/0021-9991(90)90233-Q.
Corey, A.T. 1954. The Interrelation Between Gas and Oil RelativePermeabilities. Producers Monthly 19 (November): 38-41.
Crandall, M. and Majda, A. 1980. Monotone DifferenceApproximations for Scalar Conservation Laws. Mathematics ofComputation 34 (149): 1-21. doi:10.1090/S0025-5718-1980-0551288-3.
Deutsch, C.V. and Journel, A.G. 1998. GSLIB: Geostatistical SoftwareLibrary and User's Guide, second edition. Oxford, UK: Oxford UniversityPress.
Edwards, M.G. 2003. HigherDimensional Wave-Oriented Upwind Schemes With Minimal Cross-Wind Diffusion.Paper SPE 79689 presented at the SPE Reservoir Simulation Symposium, Houston,3-5 February. doi: 10.2118/79689-MS.
Edwards, M.G. and Rogers, C.F. 1998. Finite Volume DiscretizationWith Imposed Flux Continuity for the General Tensor Pressure Equation.Computational Geosciences 2 (4): 259-290. doi:10.1023/A:1011510505406.
Hurtado, F. et al. 2007. AQuadrilateral Element-Based Finite-Volume Formulation for the Simulation ofComplex Reservoirs. Paper SPE 107444 presented at the Latin American &Caribbean Petroleum Engineering Conference, Buenos Aires, 15-18 April. doi:10.2118/107444-MS.
Keilegavlen, E. 2009. Robust Control Volume Methods for Reservoir Simulationon Challenging Grids. The University of Bergen, Bergen, Norway.
Keilegavlen, E., Kozdon, J., and Mallison, B. 2010. MonotoneMulti-Dimensional Upstream Weighting on General Grids. Presented at the 12thEuropean Conference on the Mathematics of Oil Recovery (ECMOR XII), Oxford, UK,6-9 September.
Koren, B. 1991. Low-Diffusion Rotated Upwind Schemes, Multigrid andDefect Correction for Steady, Multi-Dimensional Euler Flows, No. 98,265-276. Berlin, Germany: International Series of Numerical Mathematics,Springer Basel.
Kozdon, J., Mallison, B., and Gerritsen, M. 2010. Multidimensional UpstreamWeighting for Multiphase Transport in Porous Media. ComputationalGeosciences. doi: 10.1007/s10596-010-9211-5.
Kozdon, J., Mallison, B., and Gerritsen, M. 2009. Robust Multi-D TransportSchemes With Reduced Grid Orientation Effects. Transport in PorousMedia 78 (1): 47-75. doi: 10.1007/s11242-008-9281-1.
Lamine, S. and Edwards, M.G. 2010a. Higher Order Multidimensional UpwindConvection Schemes for Flow in Porous Media on Structured and UnstructuredQuadrilateral Grids. SIAM J. Sci. Comput. 32 (3):1119-1139. doi: 10.1137/080727750.
Lamine, S. and Edwards, M.G. 2010b. Multidimensional UpwindConvection Schemes for Flow in Porous Media on Structured and UnstructuredQuadrilateral Grids. Journal of Computational and AppliedMathematics 234 (7): 2106-2117. doi:10.1016/j.cam.2009.08.069.
LeVeque, R.J. 1997. WavePropagation Algorithms for Multidimensional Hyperbolic Systems. J.Comput. Phys. 131 (2): 327-353. doi:10.1006/jcph.1996.5603.
Riaz, A., Pankiewitz, C., and Meiburg, E. 2004. Linear Stability of RadialDisplacements in Porous Media: Influence of Velocity-Induced Dispersion andConcentration-Dependent Diffusion. Physics of Fluids 16(10): 3592-3598. doi: 10.1063/1.1775431.
Roe, P.L. and Sidilkover, D. 1992. Optimum Positive Linear Schemes forAdvection in Two and Three Dimensions. SIAM J. Numer. Anal. 29 (6): 1542-1568. doi: 10.1137/0729089.
Schneider, G.E. and Raw, M.J. 1987. Control VolumeFinite-Element Method for Heat Transfer and Fluid Flow Using ColocatedVariables--1. Computational Procedure. Numerical Heat Transfer, Part A:Applications 11 (4): 363-390. doi:10.1080/10407788708913560.
Shubin, G.R. and Bell, J.B. 1984. An Analysis of the GridOrientation Effect in Numerical Simulation of Miscible Displacement.Computer Methods in Applied Mechanics and Engineering 47(1-2): 47-71. doi: 10.1016/0045-7825(84)90047-1.
Tan, C.T. and Homsy, G.M. 1987. Stability of Miscible Displacementsin Porous Media: Radial Source Flow. Physics of Fluids 30 (5): 1239-1245. doi: 10.1063/1.866289.
Todd, M.R., O'Dell, P.M., and Hirasaki, G.J. 1972. Methods for Increased Accuracy inNumerical Reservoir Simulators. SPE J. 12 (6): 515-530.SPE-3516-PA. doi: 10.2118/3516-PA.
van Ransbeeck, P. and Hirsch, C. 1997. A General Analysis of 2D/3DMultidimensional Upwind Convection Schemes, Vol. 57, 251-304. New York:Notes on Numerical Fluid Mechanics (NNFM), Springer.
Yanosik, J.L. and McCracken, T.A. 1979. A Nine-Point, Finite-DifferenceReservoir Simulator for Realistic Prediction of Adverse Mobility RatioDisplacements. SPE J. 19 (4): 253-262; Trans.,AIME, 267. SPE-5734-PA. doi: 10.2118/5734-PA.
Younis, R.M., Tchelepi, H.A., and Aziz, K. 2010. Adaptively LocalizedContinuation-Newton Method--Nonlinear Solvers That Converge All the Time.SPE J. 15 (2): 526-544. SPE-119147-PA. doi:10.2118/119147-PA.