A Compositional Streamline Formulation With Compressibility Effects Conference Paper uri icon


  • Streamline simulators are gaining increasing acceptance because of their computational advantage, intuitive appeal and ability to visualize the flow patterns in 3D. Much of the current industry experience has been limited to black oil simulation, specifically two phase water-oil displacements under incompressible or slightly compressible conditions. However, given the favorable computational scaling properties of streamline models, the potential advantage for compositional simulation can be even more compelling. Although several papers have discussed compositional streamline formulation, they all suffer from a major limitation, particularly for compressible flow. The decoupling of the compositional pressure and transport equations, which form the basis for streamline simulation, has been accomplished under severely restrictive assumptions. In this paper, we examine the implications of these assumptions on the accuracy of compositional streamline simulation using a novel and rigorous treatment of compressibility. Our compositional streamline formulation builds on the recent work of Cheng et al.1 and uses the concept of effective density to redefine the bi-streamfunctions and decouple the 3D conservation equations to ID transport equations for overall compositions along streamlines. The effective density accounts for volume changes with pressure and can be conveniently traced along streamlines. The streamline equations are solved using a third order TVD scheme to minimize numerical dispersion. The phase compositions and saturations are obtained via thermodynamic flash calculations. We propose a novel 'optimal' gridding strategy for efficient solution of the ID compositional equations along streamlines. The approach is based on a bias-variance trade-off of 'slowness' along streamlines and adaptively coarsens the grid for the ID solution without compromising accuracy. This significantly reduces the number of flash calculations during the ID solutions. We compare the results from our proposed formulation with finite difference simulation and also, existing compositional streamline formulations in the literature to highlight the importance of the rigorous treatment of compressibility effects. Finally, we examine the scaling behavior of the streamline and finite difference models for compositional simulation to identify potential computational advantages for large-scale field applications. Copyright 2007, Society of Petroleum Engineers Inc.

author list (cited authors)

  • Osako, I., & Datta-Gupta, A.

citation count

  • 9

publication date

  • January 2007


  • SPE  Publisher