Numerical simulation of water-gas two-phase displacement process in unsaturated granite residual soil

CAI Peichen; QUE Yun; LI Xian

doi:10.16030/j.cnki.issn.1000-3665.202010017

2021 Vol. 48, No. 6

Article Contents

Article navigation > Hydrogeology & Engineering Geology > 2021 > 48(6): 54-63

Next Article Previous Article

CAI Peichen, QUE Yun, LI Xian. Numerical simulation of water-gas two-phase displacement process in unsaturated granite residual soil[J]. Hydrogeology & Engineering Geology, 2021, 48(6): 54-63. doi: 10.16030/j.cnki.issn.1000-3665.202010017

Citation:

CAI Peichen, QUE Yun, LI Xian. Numerical simulation of water-gas two-phase displacement process in unsaturated granite residual soil[J]. Hydrogeology & Engineering Geology, 2021, 48(6): 54-63. doi: 10.16030/j.cnki.issn.1000-3665.202010017

Numerical simulation of water-gas two-phase displacement process in unsaturated granite residual soil

School of Civil Engineering, Fuzhou University, Fuzhou, Fujian　350108, China

More Information

Corresponding author: QUE Yun, queyun_2001@fzu.edu.cn

Received Date: 12 October 2020

Revised Date: 05 March 2021

Publish Date: 15 November 2021

Abstract

Abstract

Infiltration of current in unsaturated soil is essentially a two-phase flow problem of water displacing air in the process of infiltration. The accuracy of traditional two-phase flow research methods cannot meet the needs of engineering, and it is not conducive to repetitive research. In contrast, the numerical simulation method of multiphase flow at the meso-scale can better simulate the whole dynamic process of water flooding. Mao Huan, Qren and others have achieved great results in the field of pore meso-scale research, but on the one hand, most of them focus on the study of rock multiphase flow, and there are some differences between the research object and the actual pore structure. On the other hand, the widely used pore network model method cannot directly show the change of particle velocity at any time and cannot present the phase interface movement state. In view of this, in order to reveal the dynamic percolation mechanism of water-gas two-phase displacement of unsaturated granite residual soil, this paper selected undisturbed granite residual soil in Fuzhou as the research object and studied the dynamic characteristics of two-phase displacement of undisturbed soil samples by using industrial CT scanning images and Level Set method. The results show that the Level Set method can properly capture the interface position between two immiscible fluids for meso-scale water-gas two-phase displacement simulation. The water-gas two-phase displacement process has the characteristics of large pore preferential flow, and the ‘low around’ phenomenon can easily appear in the higher ground of porosity roundness in general. The displacement rate is mainly controlled by the tortuosity of the channel, displacement speed is relatively high in the straight and wide channel. There is an obvious phenomenon of ‘preferential passage’, and its seepage time is positively correlated with the first rapid and then slow characteristics, and the maximum and minimum growth rate are 10.77 % and 1.90 %, respectively. The velocity distribution of the pore cross-section is related to the pore structure, and the phenomena of ‘reflux’ and ‘flow around’ cause the displacement velocity to drop sharply, whose decreasing degree can reach 21.62 %. The maximum displacement resistance appears at the hole wall, and the narrower the hole, the greater the resistance. The displacement efficiency is directly proportional to the displacement pressure difference, and the initial pressure growth effect is significant (up to 25.49 %, and only 1.47 % later). The research results can provide a theoretical basis for the study of the water migration mechanism of porous slopes, and also can enrich the theoretical basis of rainfall-induced landslides and are helpful in preventing natural disasters.
- unsaturated granite; residual soil /
- CT scan slice /
- porous media /
- Level Set /
- water-gas two-phase displacement

FullText(HTML)

References

[1]	陈嘉伟, 高游, 付俊杰, 等. 不同类型黏土的强度特性及其预测[J]. 水文地质工程地质,2020,47(3):101 − 106. [CHEN Jiawei, GAO You, FU Junjie, et al. Strength of different clayey soils and its prediction[J]. Hydrogeology & Engineering Geology,2020,47(3):101 − 106. (in Chinese with English abstract) Google Scholar
[2]	郑晓磊, 刘志峰, 王晓宏, 等. 二维非均匀多孔介质中不可压两相驱替的有限分析算法[J]. 计算物理,2015,32(5):586 − 594. [ZHENG Xiaolei, LIU Zhifeng, WANG Xiaohong, et al. Finite analytic numerical method for two-Phase incompressible flow in 2D teterogeneous porous media[J]. Chinese Journal of Computational Physics,2015,32(5):586 − 594. (in Chinese with English abstract) doi: 10.3969/j.issn.1001-246X.2015.05.011 CrossRef Google Scholar
[3]	乔文静, 叶淑君, 吴吉春. 非均质孔隙介质中两相流的光透法应用研究[J]. 水文地质工程地质,2015,42(2):112 − 119. [QIAO Wenjing, YE Shujun, WU Jichun. A study of the two-phase flow in heterogeneous porous media with the light transmission method[J]. Hydrogeology & Engineering Geology,2015,42(2):112 − 119. (in Chinese with English abstract) Google Scholar
[4]	乔梁. 泄洪洞突扩突跌掺气水流特性的数值模拟研究[D]. 昆明: 昆明理工大学, 2014. Google Scholar QIAO Liang. Numerical simulation study on aerated water flow characteristics of sudden expansion and sudden fall of spillway tunnel[D]. Kunming: Kunming University of Science and Technology, 2014. (in Chinese with English abstract) Google Scholar
[5]	于明旭, 朱维耀, 宋洪庆. 低渗透储层可视化微观渗流模型研制[J]. 辽宁工程技术大学学报(自然科学版),2013,32(12):1646 − 1650. [YU Mingxu, ZHU Weiyao, SONG Hongqing. Development of microscopic visualization flow model of low-permeability reservoir[J]. Journal of Liaoning Technical University(Natural Science),2013,32(12):1646 − 1650. (in Chinese with English abstract) doi: 10.3969/j.issn.1008-0562.2013.12.014 CrossRef Google Scholar
[6]	王伟, 赵永攀, 江绍静, 等. 鄂尔多斯盆地特低渗油藏CO₂非混相驱实验研究[J]. 西安石油大学学报(自然科学版),2017,32(6):87 − 92. [WANG Wei, ZHAO Yongpan, JIANG Shaojing, et al. Experimental study on CO₂immiscible flooding in ultra-low permeability reservoirs, Ordos Basin[J]. Journal of Xi’an Shiyou University(Natural Science Edition),2017,32(6):87 − 92. (in Chinese with English abstract) Google Scholar
[7]	吕伟峰, 冷振鹏, 张祖波, 等. 应用CT扫描技术研究低渗透岩心水驱油机理[J]. 油气地质与采收率,2013,20(2):87 − 90. [LYU Weifeng, LENG Zhenpeng, ZHANG Zubo, et al. Study on waterflooding mechanism in low-permeability cores using CT scan technology[J]. Petroleum Geology and Recovery Efficiency,2013,20(2):87 − 90. (in Chinese with English abstract) doi: 10.3969/j.issn.1009-9603.2013.02.022 CrossRef Google Scholar
[8]	毛欢. 基于格子Boltzmann方法的多孔介质内两相流流动特性的研究[D]. 东营: 中国石油大学(华东), 2017. Google Scholar MAO Huan. Study on characteristics of two-phase flow in porous media with Lattice Boltzmann method[D]. Dongying: China University of Petroleum (Huadong), 2017. (in Chinese with English abstract) Google Scholar
[9]	ØREN P E, BAKKE S. Process based reconstruction of sandstones and prediction of transport properties[J]. Transport in Porous Media,2002,46(2/3):311 − 343. Google Scholar
[10]	陈民锋, 姜汉桥. 基于孔隙网络模型的微观水驱油驱替特征变化规律研究[J]. 石油天然气学报(江汉石油学院学报),2006,28(5):91 − 95. [CHEN Minfeng, JIANG Hanqiao. Law of characteristic variations of microscopic water displacement based on pore network model[J]. Journal of Oil and Gas Technology,2006,28(5):91 − 95. (in Chinese with English abstract) Google Scholar
[11]	朱光普, 姚军. 基于相场模型的表面活性剂两相流动模拟研究[C] //中国力学学会. 第十届全国流体力学学术会议论文摘要集. 杭州: 中国力学学会, 2018: 212-212. Google Scholar ZHU Guangpu, YAO Jun. Research on two-phase flow simulation of surfactant based on phase field model[C] // Chinese Society of Mechanics. Collection of abstracts of the 10th National Conference on fluid mechanics. Hangzhou: Chinese Society of Mechanics, 2018: 212-212. (in Chinese) Google Scholar
[12]	OSHER S, SETHIAN J A. Fronts propagating with curvature-dependent speed: Algorithms based on Hamilton-Jacobi formulations[J]. Journal of Computational Physics,1988,79(1):12 − 49. doi: 10.1016/0021-9991(88)90002-2 CrossRef Google Scholar
[13]	王琳琳, 田辉, 李国君. 基于Level Set方法对油水和气水两相界面的数值模拟[J]. 应用力学学报,2010,27(2):298 − 302. [WANG Linlin, TIAN Hui, LI Guojun. Numerical simulation of oil-water and air-water two-phase flow based on level set methods[J]. Chinese Journal of Applied Mechanics,2010,27(2):298 − 302. (in Chinese with English abstract) Google Scholar
[14]	高亚军, 姜汉桥, 王硕亮, 等. 基于Level Set有限元方法的微观水驱油数值模拟[J]. 石油地质与工程,2016,30(5):91 − 96. [GAO Yajun, JIANG Hanqiao, WANG Shuoliang, et al. Numerical simulation of microscopic water-oil displacement based on level set finite element method[J]. Petroleum Geology and Engineering,2016,30(5):91 − 96. (in Chinese with English abstract) doi: 10.3969/j.issn.1673-8217.2016.05.023 CrossRef Google Scholar
[15]	LIU J J, SONG R. Investigation of water and CO₂ flooding using pore-scale reconstructed model based on micro-CT images of Berea sandstone core[J]. Progress in Computational Fluid Dynamics,2015,15(5):317 − 326. doi: 10.1504/PCFD.2015.072013 CrossRef Google Scholar
[16]	赵晓磊, 齐秋菊, 郭春超. COMSOL与MATLAB联合仿真的实现方法[J]. 中国新技术新产品,2014(24):17 − 19. [ZHAO Xiaolei, QI Qiuju, GUO Chunchao. The realization method of COMSOL and MATLAB co-simulation[J]. New Technology & New Products of China,2014(24):17 − 19. (in Chinese with English abstract) doi: 10.3969/j.issn.1673-9957.2014.24.016 CrossRef Google Scholar
[17]	GUNDE A C, BERA B, MITRA S K. Investigation of water and CO₂ (carbon dioxide) flooding using micro-CT (micro-computed tomography) images of Berea sandstone core using finite element simulations[J]. Energy,2010,35(12):5209 − 5216. doi: 10.1016/j.energy.2010.07.045 CrossRef Google Scholar
[18]	杨永飞, 尹振, 姚军, 等. 多孔介质中水气交替注入微观渗流模拟[J]. 地球科学,2013,38(4):853 − 858. [YANG Yongfei, YI Zhen, YAO Jun, et al. Pore-scale simulation of microcosmic flow during water-alternating-gas (WAG) in porous media[J]. Earth Science,2013,38(4):853 − 858. (in Chinese with English abstract) Google Scholar
[19]	张石峰, 张冀生, 于炯. 两相流模型数值解[J]. 新疆工学院学报,1996,17(1):45 − 51. [ZHANG Shifeng, ZHANG Jisheng, YU Jiong. Numerical solution of two-phase flow model[J]. Journal of Xinjiang Institute of Technology,1996,17(1):45 − 51. (in Chinese with English abstract) Google Scholar
[20]	张丽萍, 陈儒章, 邬燕虹, 等. 风化花岗岩坡地土壤剖面大孔隙特性的空间分布[J]. 土壤学报,2018,55(3):620 − 632. [ZHANG Liping, CHEN Ruzhang, WU Yanhong, et al. Spatial distribution of macro-pore properties in soil profile on a slope of weathering granite[J]. Acta Pedologica Sinica,2018,55(3):620 − 632. (in Chinese with English abstract) doi: 10.11766/trxb201709070266 CrossRef Google Scholar
[21]	冯其红, 赵蕴昌, 王森, 等. 基于相场方法的孔隙尺度油水两相流体流动模拟[J]. 计算物理,2020,37(4):439 − 447. [FENG Qihong, ZHANG Yunchang, WANG Sen, et al. Pore-scale oil-water two-phase flow simulation based on phase pield method[J]. Chinese Journal of Computational Physics,2020,37(4):439 − 447. (in Chinese with English abstract) Google Scholar
[22]	吴丰, 姚聪, 丛林林, 等. 岩石气水两相渗流的玻璃刻蚀驱替实验与有限元数值模拟对比[J]. 岩性油气藏,2019,31(4):121 − 132. [WU Feng, YAO Cong, CONG Linlin, et al. Comparison of glass etching displacement experiment and finite element numerical simulation for gas-water two-phase seepage in rocks[J]. Lithologic Reservoirs,2019,31(4):121 − 132. (in Chinese with English abstract) Google Scholar
[23]	张鹏伟, 胡黎明, MEEGODA Jay, 等. 基于岩土介质三维孔隙结构的两相流模型[J]. 岩土工程学报,2020,42(1):37 − 45. [ZHANG Pengwei, HU Liming, MEEGODA Jay, et al. Two-phase flow model based on 3D pore structure of geomaterials[J]. Chinese Journal of Geotechnical Engineering,2020,42(1):37 − 45. (in Chinese with English abstract) Google Scholar