Streamline Tracing on General Triangular or Quadrilateral Grids
- Sebastien Francois Matringe (Stanford University) | Ruben Juanes (Massachusetts Inst. of Tech.) | Hamdi A. Tchelepi (Stanford University)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- June 2007
- Document Type
- Journal Paper
- 217 - 233
- 2007. Society of Petroleum Engineers
- 5.5 Reservoir Simulation, 5.3.2 Multiphase Flow, 1.2.3 Rock properties, 4.3.4 Scale, 5.1 Reservoir Characterisation, 5.5.7 Streamline Simulation, 4.3.1 Hydrates, 5.5.8 History Matching
- 2 in the last 30 days
- 475 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
Streamline methods have received renewed interest over the past decade as an attractive alternative to traditional finite-difference (FD) simulation. They have been applied successfully to a wide range of problems including production optimization, history matching, and upscaling. Streamline methods are also being extended to provide an efficient and accurate tool for compositional reservoir simulation. One of the key components in a streamline method is the streamline tracing algorithm. Traditionally, streamlines have been traced on regular Cartesian grids using Pollock's method. Several extensions to distorted or unstructured rectangular, triangular, and polygonal grids have been proposed. All of these formulations are, however, low-order schemes.
Here, we propose a unified formulation for high-order streamline tracing on unstructured quadrilateral and triangular grids, based on the use of the stream function. Starting from the theory of mixed finite-element methods (FEMs), we identify several classes of velocity spaces that induce a stream function and are therefore suitable for streamline tracing. In doing so, we provide a theoretical justification for the low-order methods currently in use, and we show how to extend them to achieve high-order accuracy. Consequently, our streamline tracing algorithm is semi-analytical: within each gridblock, the streamline is traced exactly. We give a detailed description of the implementation of the algorithm, and we provide a comparison of low- and high-order tracing methods by means of representative numerical simulations on 2D heterogeneous media.
Streamline simulation is now accepted as a practical tool for reservoir simulation. It represents a fast alternative to the classical FD or finite-volume (FV) methods. However, streamline simulation is still a young technology and does not offer the same capabilities as more traditional methods. Here, we investigate the extension of the streamline method to simulate problems on unstructured or highly distorted grids with full tensor permeability fields.
In streamline simulation, the flow problem (pressure equation) and the transport problem (saturation equations) are solved sequentially in an operator-splitting fashion. The transport problem is solved along the streamlines using a 1D formulation of the transport equation expressed in terms of the time-of-flight variable (Bradvedt et al. 1993; Batycky et al. 1997; King and Datta-Gupta 1998). A background simulation grid is used to solve the flow problem and trace the streamlines. Therefore, extension of the streamline method to general triangular or quadrilateral grids hinges on the ability to: (1) properly discretize the pressure equation, and (2) accurately trace the streamlines on these advanced grids.
These two problems are linked. The key link between discretization and streamline tracing resides in the velocity field description. To each discretization corresponds a particular form of velocity field, and the streamline tracing algorithm has to be adapted to each type of velocity field.
|File Size||1 MB||Number of Pages||7|
Aavatsmark, I. 2002. An Introduction to Multipoint Flux Approximation forQuadrilateral Grids. Comput. Geosci. 6: 405-432. DOI: http://dx.doi.org/10.1023/A:1021291114475.
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: 10.2118/36726-PA.
Bratvedt, F. et al. 1993. Frontline and Frontsim, two full scale, two-phase,black-oil reservoir simulators based on front tracking. Surv. Math. Ind.3: 185-215.
Brezzi, F., Douglas, J., and Marini, L.D. 1985. Two families of mixed finiteelements for second order elliptic problems. Numer. Math. 47:217-235.
Brezzi, F. and Fortin, M. 1991. Mixed and Hybrid Finite Element Methods. NewYork City: Springer-Verlag.
Cordes, C. and Kinzelbach, W. 1992. Continuous Groundwater Velocity Fieldsand Path Lines in Linear, Bilinear, and Trilinear Finite Elements. WaterResour. Res. 28 (11): 2903-2911.
Durlofsky, L.J. 1991. Accuracy of Mixed and Control Volume Finite ElementApproximations to Darcy Velocity and Related Quantities. Water Resour.Res. 30 (4): 329-348.
Fortin, M and Glowinski, R. 1983. Augmented Lagrangian Methods:Applications to the Numerical Solution of Boundary-Value Problems.Amsterdam: North-Holland Co.
Hægland, H. 2003. Streamline Tracing on Irregular Girds. Masters Thesis.University of Bergen, Bergen, Norway.
Jimenez, E., Sabir, K., Datta-Gupta, A., and King, M.J. 2005. Spatial Error and Convergence inStreamline Simulation. Paper SPE 92873 presented at the SPE ReservoirSimulation Symposium, The Woodlands, Texas, 31 January-2 February. DOI:10.2118/92873-MS.
Juanes, R. and Matringe, S.F. Unified Formulation of Velocity Fields forStreamline Tracing on Two-Dimensional Unstructured Grids. Submitted forpublication.
King, M.J. and Datta-Gupta, A. 1998. Streamline Simulation: A CurrentPerspective. In Situ 22 (1): 91-140.
Matringe, S.F. and Gerritsen, M.G. 2004. On Accurate Tracing ofStreamlines. Paper SPE 89920 presented at the SPE Annual TechnicalConference and Exhibition, Houston, 26-29 September. DOI: 10.2118/89920-MS.
Mosé, R., Siegel, P., Ackerer, P., and Chavent, G. 1994. Application of theMixed Hybrid Finite Element Approximation in a Groundwater Flow Model: Luxuryor Necessity? Water Resour. Res. 3 (11): 3001-3012.
Pollock, D.W. 1988. Semianalytical Computation of Path Lines forFinite-Difference Models. Ground Water 26 (6): 743-750.
Prevost, M., Edwards, M.G., and Blunt, M.J. 2002. Streamline Tracing on CurvilinearStructured and Unstructured Grids. SPEJ 7 (2): 139-148.SPE-78663-PA. DOI: 10.2118/78663-PA.
Raviart, P.A. and Thomas, J.M. 1977. A Mixed Finite Element Method forSecond Order Elliptic Problems. InMathematical Aspects of the FiniteElement Method, Lecture Notes in Mathematics 606. New York:Springer-Verlag.