基于相场方法的孔隙尺度油水两相流体流动模拟

计算物理 ›› 2020, Vol. 37 ›› Issue (4): 439-447.

基于相场方法的孔隙尺度油水两相流体流动模拟

冯其红¹, 赵蕴昌¹, 王森¹, 张以根², 孙业恒², 史树彬³

1. 中国石油大学(华东)石油工程学院, 山东青岛 266580;
2. 中国石化胜利油田分公司勘探开发研究院, 山东东营 257000;
3. 中国石化胜利油田分公司石油工程技术研究院, 山东东营 257000

收稿日期:2019-03-11 修回日期:2019-05-28 出版日期:2020-07-25 发布日期:2020-07-25
通讯作者: 王森(1986-),副教授,主要从事多尺度多相流体流动与传热研究,E-mail:wangsena1@126.com
作者简介:冯其红(1969-),教授,博士生导师,研究方向:油气田开发工程,E-mail:fengqihong@126.com
基金资助:
国家重大科技专项（2016ZX05011-001）、国家自然科学基金（51474233、U1762213）及中石化"十条龙"项目"胜利整装油田特高含水期深度堵调技术"（217007）资助

Pore-scale Oil-Water Two-phase Flow Simulation Based on Phase Field Method

FENG Qihong¹, ZHAO Yunchang¹, WANG Sen¹, ZHANG Yigen², SUN Yeheng², SHI Shubin³

1. School of Petroleum Engineering, China University of Petroleum(East China), Qingdao, Shandong 266580, China;
2. Exploration and Development Research Institute, SINOPEC Shengli Oilfield Company, Dongying, Shandong 257000, China;
3. Petroleum Engineering Technology Research Institute, SINOPEC Shengli Oilfield Company, Dongying, Shandong 257000, China

Received:2019-03-11 Revised:2019-05-28 Online:2020-07-25 Published:2020-07-25

摘要/Abstract

摘要： 基于真实岩心颗粒粒径分布，利用过程法构建疏松砂岩油藏的三维孔隙结构模型，利用相场方法建立两相流体流动数学模型并利用有限元方法进行求解，研究驱替速度、流体性质、润湿性对剩余油分布以及采出程度的影响.结果表明：驱替速度的增大和油水粘度比的减小会导致较大的毛管数，进而有利于采出程度的提高；就润湿性而言，水湿条件下毛管力是水驱油的动力，而在油湿条件下是阻力，因此水湿岩心采出程度更高.同时，从孔隙尺度对油水渗流机理及剩余油分布机理进行揭示，结果表明：由于多孔介质的复杂孔隙结构，流体在流经不同孔隙时呈现不同的流动特征，进而对油水两相流整体的压力分布、流速分布造成重要影响.

关键词: 多相流, 剩余油, 多孔介质, 相场方法, 提高采收率

Abstract: Based on a typical particle size distribution of conventional sandstone, we establish a three-dimensional porous media of unconsolidated sandstone reservoir with process method. By using phase field method(PFM), a two-phase flow model is established on the basis of two-dimensional porous media model, which is solved with finite element method. Influences of displacement velocity, fluid properties and wettability on remaining oil distribution and recovery are discussed. It shows that displacement speed increase and inviscosity ratio decrease lead to a higher capillary number, thereby facilitating enhancing oil recovery. As for wettability, capillary force is the driving force for displacement under-water wet condition and it is a resistance under oil-wet condition,which results in a higher oil recovery in water-wet core.Meanwhile, this study reveals mechanism of two-phase flow and remaining oil distribution in pore scale. It can be seen that due to complex pore structure, fluid exhibits various characteristics as flowing through pores,which have an important impact on pressure distribution and fluid velocity distribution of oil-water two-phase flow.

Key words: multiphase flow, remaining oil, porous media, phase field method, enhanced oil recovery

中图分类号:

TE312

冯其红, 赵蕴昌, 王森, 张以根, 孙业恒, 史树彬. 基于相场方法的孔隙尺度油水两相流体流动模拟[J]. 计算物理, 2020, 37(4): 439-447.

FENG Qihong, ZHAO Yunchang, WANG Sen, ZHANG Yigen, SUN Yeheng, SHI Shubin. Pore-scale Oil-Water Two-phase Flow Simulation Based on Phase Field Method[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2020, 37(4): 439-447.

参考文献

[1] 闫伟超, 孙建孟. 微观剩余油研究现状分析[J]. 地球物理学进展, 2016, 31(5):2198-2211.
[2] 胡雪涛, 李允. 随机网络模拟研究微观剩余油分布[J]. 石油学报, 2000, 21(4):46-51.
[3] 赵福麟. EOR原理[M]. 石油工业大学出版社, 2001:10-11.
[4] 马炳杰. 注水停止对剩余油分布状态影响的实验研究[J]. 科学技术与工程, 2016, 16(12):254-259.
[5] 李洪玺, 刘全稳, 王健, 等. 弱凝胶流动规律及其提高采收率机理新认识[J]. 油气地质与采收率, 2006, 13(4):85-87.
[6] 许长福, 刘红现, 钱根宝, 等. 克拉玛依砾岩储集层微观水驱油机理[J]. 石油勘探与开发, 2011, 38(6):725-732.
[7] WANG S, FENG Q, DONG Y L, et al. A dynamic pore-scale network model for two-phase imbibition[J]. Journal of Natural Gas Science and Engineering, 2015, 26:118-129.
[8] SUCCI S. The lattice Boltzmann equation:For fluid dynamics and beyond[M]. Oxford University Press, 2001.
[9] LI J, HO M T, WU L, et al. On the unintentional rarefaction effect in LBM modeling of intrinsic permeability[J]. Advances in Geo-Energy Research, 2018, 2(4):404-409.
[10] HIRT C W, NICHOLS B D. Volume of fluid (VOF) method for the dynamics of free boundaries[J]. Journal of Computational Physics, 1981, 39(1):201-225.
[11] AMIRI H A A, HAMOUDA A A. Evaluation of level set and phase field methods in modeling two phase flow with viscosity contrast through dual-permeability porous medium[J]. International Journal of Multiphase Flow, 2013, 52:22-34.
[12] HATIBOGLU C U, BABADAGLI T. Pore-scale studies of spontaneous imbibition into oil-saturated porous media[J]. Physical Review E, 2008, 77(6):066311.
[13] BLUNT M J. Flow in porous media-pore-network models and multiphase flow[J]. Current Opinion in Colloid & Interface Science, 2001, 6(3):197-207.
[14] CAI J, WEI W, HU X, et al. Fractal characterization of dynamic fracture network extension in porous media[J]. Fractals, 2017, 25(2):1750023.
[15] AMIRI H A A, HAMOUDA A A. Pore-scale modeling of non-isothermal two phase flow in 2D porous media:Influences of viscosity, capillarity, wettability and heterogeneity[J]. International Journal of Multiphase Flow, 2014, 61:14-27.
[16] GAO Y J, WANG J F, LUO Z R, et al. Nano-twin structure simulation with phase field crystal method[J]. Chinese Journal of Computational Physics, 2013, 30(4):577-581.
[17] GAO Y J, YANG R L, HUANG L L, et al. Phase field crystal simulation in nano-scale for crack extension with pre-deformation[J]. Chinese Journal of Computational Physics, 2017, 34(4):453-460.
[18] ZHU G, YAO J, ZHANG L, et al. Investigation of the dynamic contact angle using a direct numerical simulation method[J]. Langmuir, 2016, 32(45):11736-11744.
[19] ROKHFOROUZ M R, AKHLAGHIAMIRIH A. Phase-field simulation of counter-current spontaneous imbibition in a fractured heterogeneous porous medium[J]. Physics of Fluids, 2017, 29(6):062104.
[20] ZHANG N, YAO J, HUANG Z Q, et al. Locally conservative Galerkin numerical simulation for two-phase flow in porous media[J]. Chinese Journal of Computational Physics, 2013, 30(5):667-674.
[21] BAKKE S, OREN P E. 3-D pore-scale modelling of sandstones and flow simulations in the pore networks[J]. SPE Journal, 1997, 2(2):136-149.
[22] REN P E, BAKKE S. Process based reconstruction of sandstones and prediction of transport properties[J]. Transport in Porous Media, 2002, 46(2):311-343.
[23] BADALASSI V E, CENICEROS H D, BANERJEE S. Computation of multiphase systems with phase field models[J]. Journal of Computational Physics, 2003, 190(2):371-397.
[24] WASHBURN E W. The dynamics of capillary flow[J]. Physical Review, 1921, 17(3):273.
[25] LIU H, VALOCCHI A J, WERTH C, et al. Pore-scale simulation of liquid CO₂ displacement of water using a two-phase lattice Boltzmann model[J]. Advances in Water Resources, 2014, 73:144-158.
[26] 姜汉桥, 姚军, 姜瑞忠. 油藏工程原理与方法[M]. 东营:石油大学出版社, 2000:109-117.
[27] MENG Q, LIU H, WANG J, et al. Effect of wetting-phase viscosity on cocurrent spontaneous imbibition[J]. Energy & Fuels, 2016, 30(2):835-843.