Streamline Tracing on General Triangular or Quadrilateral Grids

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 lowand high-order tracing methods by means of representative numerical simulations on 2D heterogeneous media.

[1]  D. W. Pollock Semianalytical Computation of Path Lines for Finite‐Difference Models , 1988 .

[2]  Michael G. Edwards,et al.  Streamline Tracing on Curvilinear Structured and Unstructured Grids , 2002 .

[3]  W. Kinzelbach,et al.  Continuous Groundwater Velocity Fields and Path Lines in Linear, Bilinear, and Trilinear Finite Elements , 1992 .

[4]  Margot Gerritsen,et al.  On Accurate Tracing of Streamlines , 2004 .

[5]  M. Fortin,et al.  Augmented Lagrangian methods : applications to the numerical solution of boundary-value problems , 1983 .

[6]  I. Aavatsmark,et al.  An Introduction to Multipoint Flux Approximations for Quadrilateral Grids , 2002 .

[7]  L. Durlofsky Accuracy of mixed and control volume finite element approximations to Darcy velocity and related quantities , 1994 .

[8]  L. D. Marini,et al.  Two families of mixed finite elements for second order elliptic problems , 1985 .

[9]  Martin J. Blunt,et al.  A 3D Field-Scale Streamline-Based Reservoir Simulator , 1997 .

[10]  G. Chavent,et al.  Application of the mixed hybrid finite element approximation in a groundwater flow model: Luxury or necessity? , 1994 .

[11]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[12]  Akhil Datta-Gupta,et al.  Spatial Error and Convergence in Streamline Simulation , 2007 .

[13]  Michael G. Edwards,et al.  Streamline Tracing on Curvilinear Structured and Unstructured Grids , 2001 .

[14]  H. Hægland,et al.  Streamline Tracing on Irregular Grids , 2003 .

[15]  Philippe G. Ciarlet,et al.  7. A Mixed Finite Element Method , 2002 .