Theory and application of numerical simulation method of capillary force enhanced oil production

doi:10.1007/s10483-015-1917-6

Abstract

Abstract: A kind of second-order implicit upwind fractional step finite difference methods are presented for the numerical simulation of coupled systems for enhanced (chemical) oil production with capillary force in the porous media. Some techniques, e.g., the calculus of variations, the energy analysis method, the commutativity of the products of difference operators, the decomposition of high-order difference operators, and the theory of a priori estimate, are introduced. An optimal order error estimate in the l² norm is derived. The method is successfully used in the numerical simulation of the enhanced oil production in actual oilfields. The simulation results are satisfactory and interesting.

Key words: second-order implicit upwind fractional step difference, application in actual oilfield, optimal order l² estimate, three-dimensional porous coupled system on consideration capillary force, enhanced (chemical) oil production

2010 MSC Number:

O241.82

Yirang YUAN;Aijie CHENG;Danping YANG;Changfeng LI;Yunxin LIU. Theory and application of numerical simulation method of capillary force enhanced oil production. Applied Mathematics and Mechanics (English Edition), 2015, 36(3): 379-400.

TrendMD

References

[1] Ewing, R. E., Yuan, Y. R., and Li, G. Finite element for chemical-flooding simulation. Proceeding of the 7th International Conference of Finite Element Method in Flow Problems, The University of Alabama in Huntsville, Huntsville, 1264-1271 (1989)
[2] Yuan, Y. R. Theory and Application of Numerical Simulation of Energy Sources, Basis of Numerical Simulation of Chemical Production (Tertiary Oil Recovery), Science Press, Beijing, 257-304 (2013)
[3] Yuan, Y. R., Yang, D. P., and Qi, L. Q. Research on algorithms of applied software of the polymer. Proceedings on Chemical Flooding, Petroleum Industry Press, Beijing, 246-253 (1998)
[4] Douglas, J., Jr. and Russell, T. F. Numerical method for convection-dominated diffusion problems based on combining the method of characteristics with finite element or finite difference procedures. SIAM Journal on Numerical Analysis, 19, 871-885 (1982)
[5] Douglas, J., Jr. Simulation of miscible displacement in porous media by a modified method of characteristic procedure. Lecture Notes in Mathematics, Springer, Berlin/New York, 64-70 (1982)
[6] Douglas, J., Jr. Finite difference methods for two-phase incompressible flow in porous media. SIAM Journal on Numerical Analysis, 20, 681-696 (1983)
[7] Ewing, R. E., Russell, T. F., and Wheeler, M. F. Convergence analysis of an approximation of miscible displacement in porous media by mixed finite elements and a modified method of characteristics. Computer Methods in Applied Mechanics and Engineering, 47(1-2), 73-92 (1984)
[8] Douglas, J., Jr. and Roberts, J. E. Numerical method for a model for compressible miscible displacement in porous media. Mathematics of Computation, 4, 441-459 (1983)
[9] Yuan, Y. R. Time stepping along characteristics for the finite element approximation of compressible miscible displacement in porous media. Mathematica Numerica Sinica, 14, 386-406 (1992)
[10] Yuan, Y. R. Finite difference methods for a compressible miscible displacement problem in porous media. Mathematica Numerica Sinica, 15, 16-28 (1993)
[11] Ewing, R. E. The Mathematics of Reservior Simulation, Society for Industrial and Applied Mathematics, Philadelphia (1983)
[12] Yuan, Y. R. Characteristic finite difference methods for moving boundary value problem of numerical simulation of oil deposit. Science in China, 37, 1442-1453 (1994)
[13] Yuan, Y. R. The characteristic mixed-finite element method and analysis for three-dimensional moving boundary value problem. Science in China, 39, 276-288 (1996)
[14] Yuan, Y. R. Finite difference method and analysis for three-dimensional semiconductor device of heat conduction. Science in China, 39, 1140-1151 (1996)
[15] Axelsson, O. and Gustafasson, I. A modified upwind scheme for convective transport equations and the use of a conjugate gradient method for the solution of non-symmetric systems of equations. Journal of the Institute of Mathematics and Its Applications, 23, 321-337 (1979)
[16] Ewing, R. E., Lazarov, R. D., and Vassilevski, A. T. Finite difference shceme for parabolic problems on composite grids with refinement in time and space. SIAM Journal on Numerical Analysis, 31, 1605-1622 (1994)
[17] Lazarov, R. D., Mishev, I. D., and Vassilevski, P. S. Finite volume method for convection-diffusion problems. SIAM Journal on Numerical Analysis, 33, 31-55 (1996)
[18] Peaceman, D. W. Fundamantal of Numerical Reservoir Simulation, Elsevier, Amsterdam (1980)
[19] Douglas, J., Jr. and Gunn, J. E. Two high-order correct difference analogues for the equation of multidimensional heat flow. Mathematics of Computation, 17, 71-80 (1963)
[20] Douglas, J., Jr. and Gunn, J. E. A general formulation of alternating direction methods, part 1: parabolic and hyperbolic problems. Numerische Mathematik, 6, 428-453 (1964)
[21] Marchuk, G. I. Splitting and alternating direction method. Handbook of Numerical Analysis (eds.Ciarlet, P. G. and Lions, J. L.), Elsevior Science Publishers, Paris, 197-460 (1990)
[22] Yuan, Y. R. Theory and application of upwind finite difference method for moving boundary value problem of three-dimensional percolation coupled system. Science in China, 40, 103-126 (2010)
[23] Yuan, Y. R. The second-order upwind finite difference fractional steps method for moving boundary value problem of nonlinear percolation coupled system. Science in China, 42, 845-864 (2012)
[24] Shen, P. P., Liu, M. X., and Tang, L. Mathematical Problems of Petroleum Exploration and Development: Mathematical Problems of Oil-Gas Fields Development, Science Press, Beijing, 197- 264 (2002)