A Rigorous Compressible Streamline Formulation for Two and Three-Phase Black-Oil Simulation
- Hao Cheng (Texas A&M U.) | Ichiro Osako (Texas A&M U.) | Akhil Datta-Gupta (Texas A&M U.) | Michael J. King (BP America Inc)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- December 2006
- Document Type
- Journal Paper
- 407 - 417
- 2006. Society of Petroleum Engineers
- 5.2.1 Phase Behavior and PVT Measurements, 5.1.5 Geologic Modeling, 5.4.2 Gas Injection Methods, 5.5 Reservoir Simulation, 4.3.4 Scale, 5.4.1 Waterflooding, 5.6.5 Tracers, 5.3.3 Particle Transportation, 5.3.1 Flow in Porous Media, 4.1.2 Separation and Treating, 5.1 Reservoir Characterisation, 4.6 Natural Gas, 7.6.2 Data Integration, 5.5.7 Streamline Simulation, 3.3.6 Integrated Modeling, 5.3.2 Multiphase Flow
- 3 in the last 30 days
- 773 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
Streamline simulators have received increased attention in the petroleum industry because of their ability to effectively handle multimillion-cell detailed geologic models and large simulation models. The efficiency of streamline simulation has relied primarily on the decoupling of the 3D saturation equation into 1D equations along streamlines using the streamline time of flight as the spatial coordinate. Until now, this decoupling has been strictly valid for incompressible flow. Applications to compressible flow have generally lacked strong theoretical foundations, and very often yielded mixed or unsatisfactory results.
In this paper, for the first time we generalize streamline models to compressible flow using a rigorous formulation while retaining much of its favorable characteristics. Our new formulation is based on three major elements and requires only minor modifications to existing streamline models. First, we introduce an "effective density?? for the total fluids along the streamlines. This density captures the changes in the fluid volume with pressure and can be conveniently and efficiently traced along streamlines. Thus, we simultaneously compute time of flight and volume changes along streamlines. Second, we incorporate a density-dependent source term in the streamline saturation equation to account for compressibility effects. Third, the effective density, fluid volumes, and the time-of-flight information are used to incorporate cross-streamline effects through use of pressure updates and remapping of saturations. Our proposed approach preserves the 1D nature of the saturation calculations and all the associated advantages of the streamline approach. The saturation calculations are fully decoupled from the underlying grid and can be carried out using large timesteps without grid-based stability limits.
We demonstrate the validity and practical utility of our approach using synthetic and field examples and comparison with a commercial finite-difference simulator. A comparison of the number of pressure solutions and the CFL numbers for the streamline and finite-difference simulation indicates that our proposed compressible streamline approach is likely to offer substantial computational advantage.
Streamline simulators have become increasingly popular for high-resolution reservoir simulation using multimillion-cell geologic models. For incompressible or slightly compressible flow and under convection-dominated conditions, streamline models are well known to outperform conventional finite-difference simulation in terms of computational speed. Streamline models can also result in improved accuracy because of subgrid resolution and reduced numerical dispersion and grid-orientation effects (King and Datta-Gupta 1998; Datta-Gupta 2000). To a large extent, the efficiency of the current streamline formulation is a consequence of the incompressibility assumption that allows us to easily and effectively decouple the pressure and saturation calculations during flow simulation. This decoupling has been greatly facilitated by the introduction of the streamline time of flight coordinate (Datta-Gupta and King 1995). Specifically, utilizing the time of flight as the spatial coordinate, the multidimensional saturation calculations are reduced to a series of 1D solutions along streamlines. These 1D solutions can be carried out independently and using relatively large timesteps, as they are not impacted by the underlying geologic grid-based stability limitations. This is the primary advantage of streamline simulation. In addition, for heterogeneity-dominated flow and adverse mobility ratio conditions, the streamlines need to be updated infrequently, leading to further savings in computation time (King and Datta-Gupta 1998; Datta-Gupta 2000).
|File Size||2 MB||Number of Pages||11|
Bear, J. 1973. Dynamics of Fluid in Porous Media. New York: DoverPublications.
Bratvedt, F., Gimse, T., and Tegnander, C. 1996. Streamline Computations forPorous Media Flow Including Gravity. Transport in Porous Media 25(1): 63-78. DOI: http://dx.doi.org/10.1007/BF00141262.
Coats, K.H. 2001. IMPESStability: The CFL Limit. Paper SPE 66345 presented at the SPE ReservoirSimulation Symposium, Houston, 11-14 February. DOI:http://dx.doi.org/10.2118/66345-MS.
Crane, M., Bratvedt, F., Bradvedt, K., Childs, P., and Olufsen, R. 2000. A Fully Compositional StreamlineSimulator. Paper SPE 63156 presented at the SPE Annual Technical Conferenceand Exhibition, Dallas, 1-4 October. DOI:http://dx.doi.org/10.2118/63156-MS.
Datta-Gupta, A. 2000. Streamline Simulation: A TechnologyUpdate. JPT 52 (12): 65-73, 84. SPE-65604-PA. DOI:http://dx.doi.org/10.2118/65604-PA.
Datta-Gupta, A. and King, M.J. 1995. A Semianalytic Approach to Tracer FlowModeling in Heterogeneous Permeable Media. Advances in Water Resources18 (1): 9-24. DOI: http://dx.doi.org/10.1016/0309-1708(94)00021-V.
He, Z., Yoon, S., and Datta-Gupta, A. 2002. Streamline-Based Production DataIntegration Under Changing Field Conditions. SPEJ 7 (4): 423-436.SPE-71333-PA. DOI: http://dx.doi.org/10.2118/71333-PA.
Ingebrigtsen, L., Bratvedt, F., and Berge, J. 1999. A Streamline-Based Approach toSolution of Three-Phase Flow. Paper SPE 51904 presented at the 1999 SPEReservoir Simulation Symposium, Houston, 14-17 February. DOI:http://dx.doi.org/10.2118/51904-MS.
Jessen, K. and Orr, F.M. Jr. 2002. Compositional StreamlineSimulation. Paper SPE 77379 presented at the SPE Annual TechnicalConference and Exhibition, San Antonio, Texas, 29 September-2 October. DOI:http://dx.doi.org/10.2118/77379-MS.
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:http://dx.doi.org/10.2118/92873-MS.
King, M.J. and Datta-Gupta, A. 1998. Streamline simulation: A currentperspective. In Situ 22 (1): 91-140.
Martin, J.C. and Wegner, R.E. 1979. Numerical Solution of Multiphase, 2DIncompressible Flow Using Stream-Tube Relationships. SPEJ 19(5): 313-323. SPE-7140-PA. DOI: http://dx.doi.org/10.2118/7140-PA.
Osako, I., Datta-Gupta, A., and King, M.J. 2004. Timestep Selection During StreamlineSimulation Through Transverse Flux Correction. SPEJ 9(4):450-462. SPE-79688-PA. DOI: http://dx.doi.org/10.2118/79688-PA.
Pollock, D.W. 1988. Semianalytical Computation of Path Lines forFinite-Difference Models. Ground Water 26 (6): 743-750. DOI: http://dx.doi.org/10.1111/j.1745-6584.1988.tb00425.x.
Ponting, D.K. 1998. HybridStreamline Methods. Paper SPE 39756 presented at the SPE Asia PacificConference on Integrated Modeling for Asset Management, Kuala Lumpur, 23-24March. DOI: http://dx.doi.org/10.2118/39756-MS.
Thiele, M.R., Batycky, R.P., and Blunt, M.J. 1997. A Streamline-Based 3D Field-ScaleCompositional Reservoir Simulator. Paper SPE 38889 presented at the SPEAnnual Technical Conference and Exhibition, San Antonio, Texas, 5-8 October.DOI: http://dx.doi.org/10.2118/38889-MS.