研究了利用迦利金有限元方法模拟油藏中油水两相渗流问题，提供一种油藏数值模拟方法。研究方法从数学模型出发，严格推导了油藏油水两相渗流方程的有限元弱形式，在井边界处理及上游权选择等非线性处理方面给出了可行的策略。采用商业性有限元代码生成平台FEPG，生成六面体八结点有限元软件代码，模拟三维两相油藏的油水渗流规律。为验证有限元方法及软件的正确性，与基于有限差分的商业软件进行对比，结果表明有限元方法模拟计算结果有效可靠，同时网格弥散现象有所控制。
This paper simulates the oil-water two-phase flow of oil and water based on the Galerkin finite element method, as an advanced reservoir simulation method. The finite element mathematical model in the weak form is derived for the oil-water two-phase flow, Then, a feasible strategy is adopted to deal with the nonlinearity terms in the well boundary processing and the upstream weight selection. Using the finite element code generation platform FEPG, the finite element software code based on hexahedral eight nodes is generated for the simulation of the three-dimensional two-phase oil-water percolation. To validate the finite element method and the software, the results are compared with those obtained by the commercial software based on the finite difference, and it is shown that the finite element method is effective and reliable.
[1] Sebastian G. Combining finite volume and finite element methods to simulate fluid flow in geologic media[J]. ANZIAM Joural, 2002, 44: 180-201.
[2] Durlofsky L J. A triangle based mixed finite element finite volume technique for modeling two phase through porous media[J]. Journal of ComputationalPhysics, 1993, 105:252-266.
[3] Aarnes J E, Kippe V, Knut-Andreas L. Mixed multiscale finite elements and streamline methods for reservoir simulation of large geomodels[J]. AdvancesinWaterResources, 2005, 28(3):257-271.
[4] Kozdon J, Mallison B, Gerritsen M, et al. Multi-D upwinding for multi phase transport in porous media[C]//SPE Reservoir Simulation Symposium. Woodlands, Texas:SocietyofPetroleumEngineers, 2009.
[5] Efendiev Y, Ginting V, Hou T, et al. Accurate multiscale finite element methods for two-phase flow simulations[J]. Journal of Computational Physics, 2006220:155-174.
[6] Chen Z X. Finite element methods for the black oil model in petroleum reservoirs[M].Minneapolis:IMAPreprintSeries, 1994.
[7] 刘振宇. 有限元法在油藏渗流中的理论和应用[D]. 安达: 大庆石油学院, 2003. Liu Zhenyu. The theory and application of finite element method for the flow inreservoir[D].Anda:DaqingPetroleumInstitute, 2003.
[8] 王树平, 李治平, 罗勇, 等. 水平井射孔孔眼分布方式优化分析[J]. 钻采工艺, 2007(2):39-45. Wang Shuping, Li Zhiping, Luo Yong, et al. Drilling & Production Technology, 2007(2):39-45.
[9] 董平川. 油气储层流固耦合理论、数值模拟及其应用[J]. 岩石力学与工程学报, 1999, 18(4):491-491. Dong Pingchuan. Chinese Journal of Rock Mechanics and Engineering, 1999, 18(4):491-491.
[10] 沈国华. 有限元数值模拟方法在构造裂缝预测中的应用[J]. 油气地质与采收率, 2008(4):24-29. Shen Guohua. Petroleum Geology and Recovery Efficiency, 2008 (4): 24-29.
[11] 宋永东, 戴俊生. 储层构造裂缝预测研究[J]. 油气地质与采收率, 2007(6):9-13. Song Yongdong, Dai Junsheng. Petroleum Geology and Recovery Efficiency, 2007(6):9-13.
[12] 张烈辉, 李成勇, 刘启国, 等. 边界元理论在缝洞型非均质油藏渗流研究中的应用[J].石油与天然气地质, 2007(4):528-534. Zhang Liehui, Li Chengyong, Liu Qiguo, et al. Oil & Gas Geology, 2007(4): 528-534.
[13] 梁国平.有限元语言[M].北京:科学出版社, 2009. Liang Guoping.Finite element language[M].Beijing:Science Press, 2009.
