Research on Numerical Simulations of Heavy Oil Thermal Recovery by Steam Injection

The works on the developments of the numerical scheme for multi-phase flows in porous media and the adaptive mesh refinement in numerical simulations of heavy oil thermal recovery by steam injection are introduced. First, a finite analytic method was developed to deal with a challenging problem in numerical simulations in the case of strong geological heterogeneity, which cause the upscaling permeability to have a rapid change cross the grid interfaces and therefore the nodal flow effects will lead the flow fingering to the high permeability region. With the traditional numerical scheme, refining the coarse grid enough is the only manner to describe the flow pattern accurately. Second, for the incompressible two-phase flows, a numerical scheme, which could avoid the grid orientation effects efficiently, was proposed. Next, the adaptive mesh refinement technique was applied to the numerical simulations for the processes of thermal recovery by steam injection, where the different heterogeneous cases of the reservoir were considered, like the reservoir with the permeability variations, the different rock-types, the fractured porous media, the complex faulted reservoir or complex boundary reservoirs. The proposed AMR technique is fast and can give good accuracy. At last, on the basis of above theoretical research, a software package for the reservoir simulations, where the AMR technique was included, was developed.

[1]  Benoit Noetinger,et al.  Calculation of Internodal Transmissivities in Finite Difference Models of Flow in Heterogeneous Porous Media , 1995 .

[2]  J. L. Yanosik,et al.  A Nine-Point, Finite-Difference Reservoir Simulator for Realistic Prediction of Adverse Mobility Ratio Displacements , 1979 .

[3]  K. Aziz,et al.  The Grid Orientation Effect in Reservoir Simulation , 1991 .

[4]  Michel Quintard,et al.  Adaptive Mesh Refinement for One-Dimensional Three-Phase Flow with Phase Change in Porous Media , 2006 .

[5]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

[6]  M. R. Todd,et al.  Methods for Increased Accuracy in Numerical Reservoir Simulators , 1972 .

[7]  J. Trangenstein Multi-scale iterative techniques and adaptive mesh refinement for flow in porous media , 2002 .

[8]  Robert W. Zimmerman,et al.  Accuracy of the Renormalization Method for Computing Effective Conductivities of Heterogeneous Media , 2001 .

[9]  Hassan Mahani,et al.  Generation of Voronoi Grid Based on Vorticity for Coarse-Scale Modeling of Flow in Heterogeneous Formations , 2010 .

[10]  Richard D. Hornung,et al.  Adaptive Mesh Refinement and Multilevel Iteration for Flow in Porous Media , 1997 .

[11]  Zhi-Feng Liu,et al.  Finite analytic numerical method for two-dimensional fluid flow in heterogeneous porous media , 2013, J. Comput. Phys..

[12]  Zhuoxin Bi,et al.  Large Multi-Scale Iterative Techniques and Adaptive Mesh Refinement for Miscible Displacement Simulation , 2002 .

[13]  Rainer Grauer,et al.  Adaptive Mesh Refinement for Singular Current Sheets in Incompressible Magnetohydrodynamic Flows , 1997 .

[14]  K. H. Coats,et al.  Effects of grid type and difference scheme on pattern steamflood simulation results , 1982 .

[15]  Haishan Luo,et al.  Oil Displacement for One-Dimensional Three-Phase Flow with Phase Change in Fractured Media , 2009 .

[16]  Pierre Haldenwang,et al.  Dynamically adapted mesh refinement for combustion front tracking , 2002 .

[17]  Stanley Osher,et al.  Minimization of Grid Orientation Effects Through Use of Higher Order Finite Difference Methods , 1993 .

[18]  Michel Quintard,et al.  Adaptive Mesh Refinement for One-Dimensional Three-Phase Flows in Heterogeneous Fractured Porous Media , 2008 .

[19]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

[20]  C. Brand,et al.  Modeling reservoir geometry with irregular grids , 1991 .

[21]  Gautam S. Shiralkar,et al.  A general formulation for simulating physical dispersion and a new nine-point scheme , 1991 .

[22]  Richard A. Dawe,et al.  Experimental displacement patterns in a 2 × 2 quadrant block with permeability and wettability heterogeneities—problems for numerical modelling , 2008 .