基于改进LBM的气液自发渗吸过程中动态润湿效应模拟  

SIMULATION OF DYNAMIC WETTING EFFECT DURING GAS-LIQUID SPONTANEOUS IMBIBITION BASED ON MODIFIED LBM

作　　者：张晟庭 李靖 陈掌星 张涛[3] 吴克柳[1] 冯东 毕剑飞 朱上 Zhang Shengting;Li Jing;Chen Zhangxing;Zhang Tao;Wu Keliu;Feng Dong;Bi Jianfei;Zhu Shang(State Key Laboratory of Petroleum Resources and Exploration,China University of Petroleum(Beijing),Beijing 102249,China;Department of Chemistry and Petroleum Engineering,University of Calgary,Alberta T2N1N4,Canada;State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation,Southwest Petroleum University,Chengdu 610500,China)

机构地区：[1]中国石油大学(北京)油气资源与探测国家重点实验室,北京102249 [2]卡尔加里大学化学与石油工程系,加拿大阿尔伯塔T2N1N4 [3]西南石油大学油气藏地质及开发工程国家重点实验室,成都610500

出　　处：《力学学报》2023年第2期355-368,共14页Chinese Journal of Theoretical and Applied Mechanics

基　　金：国家自然科学基金(52104051,52174041,52204058);博士后创新人才支持计划(BX20220350);中国博士后科学基金(2022M713459);中国石油大学(北京)科研基金(2462021QNXZ002)资助项目。

摘　　要：微通道内气液自发渗吸是广泛发生在自然界及诸多工业领域的物理现象,而动态接触角是影响整个渗吸过程的关键因素.针对该问题,本文使用改进的伪势多相流格子玻尔兹曼方法(LBM),直接捕捉微通道内气液自发渗吸过程中的实时接触角,并分析接触角的动态变化特性及其对渗吸长度的影响.首先,本文在原始的伪势多相流LBM的基础上耦合Peng-Robinson(PR)状态方程,改进流体−流体作用力以及流−固作用力格式,并采用精确差分方法将外力添加至LBM框架中.然后,通过校准模型的热力学一致性,模拟测试界面张力,静态平衡接触角等界面现象验证了模型的准确性.最后,基于建立的模拟方法,在水平方向上模拟微通道内气液自发渗吸过程.结果表明:渗吸过程中的接触角呈现动态变化特征,在渗吸初期,因受到惯性力的影响存在较大波动;随着渗吸距离的增大,其逐渐减小并趋近于静态平衡接触角.渗吸过程中的接触角与微通道尺寸及静态接触角有关,随着微通道宽度增大,实时的动态接触角与静态接触角相差大;随着静态接触角增大,实时的动态接触角与静态接触角的相差增大.此外,忽略动态接触角的Lucas-Washburn(LW)方程所预测的弯液面位置与模拟结果存在一定偏差,利用模拟得到实时动态接触角数据可以直接用于校正LW方程,校正后的LW方程预测的弯液面位置与模拟结果基本一致.Gas-liquid spontaneous imbibition in microchannels is a widely occurring physical phenomenon in nature and many industrial fields.The dynamic contact angle is the key factor affecting the whole gas-liquid imbibition process.In this work,we use a modified pseudopotential multiphase flow lattice Boltzmann method(LBM)to capture the realtime contact angle during gas-liquid spontaneous imbibition in microchannels and analyze the dynamic characteristics of the contact angle and its effects on the imbibition length.Firstly,we coupled the Peng-Robinson(PR)equation of state to the original pseudopotential multiphase flow LBM,improved the fluid-fluid interaction force and fluid-solid interaction force formats,and added the external forces to the LBM framework by using the exact difference method.Then,the accuracy of the model was verified by calibrating the thermodynamic consistency of the model and simulating interfacial phenomena such as interfacial tension and static equilibrium contact angles.Finally,based on the established simulation method,the spontaneous gas-liquid percolation process in the microchannel is simulated in the horizontal direction.The results show that the contact angle in the imbibition process is dynamic and varies greatly in the early stage of imbibition due to the inertia force.With the further increase of the imbibition distance,it gradually decreases and tends to the static equilibrium contact angle.The contact angle in the imbibition process is related to the microchannel size and the static contact angle.As the width of the microchannel increases,the difference between the dynamic contact angle and the static contact angle in real-time increases;as the static contact angle increases,the difference between the dynamic contact angle and the static contact angle in real-time increases.In addition,the Lucas-Washburn(LW)equation,which ignores the dynamic contact angle,predicts the position of the meniscus is different from the simulated results.The realtime dynamic contact angle data obtained from

关 键 词：气液自发渗吸 格子玻尔兹曼方法 动态接触角 微通道 Lucas-Washburn 方程 

分 类 号：TE319[石油与天然气工程—油气田开发工程]