Zhao Jianlin, Zhang Guangqing, Wu Keliu, Qin Feifei, Fei Linlin, Derome Dominique, Carmeliet Jan
College of Petroleum Engineering, China University of Petroleum-Beijing, Beijing 102249, China.
State Key Laboratory of Petroleum Resources and Engineering, China University of Petroleum-Beijing, Beijing 102249, China.
Langmuir. 2024 Apr 9;40(14):7364-7374. doi: 10.1021/acs.langmuir.3c03534. Epub 2024 Mar 27.
Wetting films can develop in the corners of pore structures during imbibition in a strongly wetting porous medium, which may significantly influence the two-phase flow dynamics. Due to the large difference in scales between main meniscus and corner film, accurate and efficient modeling of the dynamics of corner film remains elusive. In this work, we develop a novel two-pressure dynamic pore network model incorporating the interacting capillary bundle model to analyze the competition between the main meniscus and corner film flow in real porous media. A pore network with four-point star-shaped pore bodies and throat bonds is extracted from the real porous medium based on the pore shape factor and pore cross-sectional area, which is then decomposed into several layers of sub-pore networks, where the first layer of sub-pore network simulates the main meniscus flow while the upper layers characterize the corner film flow. The two-phase flow conductance of throat bonds for different layers of sub-pore networks are determined by high-resolution two-phase lattice Boltzmann modeling, thus inherently considering the viscous coupling effect. In addition, two artificial neural network models are developed to predict the two phases' flow conductance based on the shape of the throat cross section and the fluid properties. The accuracy of the developed model is validated with a lattice Boltzmann simulation of imbibition in a strongly wetting square tube. Then the model is used to simulate imbibition in a strongly wetting sandstone porous medium, and the competition between the main meniscus and the corner film flow is analyzed. The results show that with decreasing capillary number and viscosity ratio between wetting and nonwetting fluids, the development of the wetting corner film becomes more significant.
在强润湿性多孔介质的吸渗过程中，湿润膜会在孔隙结构的角落处形成，这可能会显著影响两相流动力学。由于主流体弯月面和角落膜之间的尺度差异很大，准确而有效地模拟角落膜的动力学仍然难以实现。在这项工作中，我们开发了一种新颖的双压力动态孔隙网络模型，该模型结合了相互作用的毛细管束模型，以分析实际多孔介质中主流体弯月面和角落膜流动之间的竞争。基于孔隙形状因子和孔隙横截面积，从实际多孔介质中提取出具有四点星形孔隙体和喉道键的孔隙网络，然后将其分解为几层子孔隙网络，其中第一层子孔隙网络模拟主流体弯月面流动，而上面几层表征角落膜流动。通过高分辨率的两相格子玻尔兹曼模型确定不同层子孔隙网络喉道键的两相流通量，从而内在地考虑了粘性耦合效应。此外，还开发了两个人工神经网络模型，以基于喉道横截面形状和流体性质预测两相的通量。通过对强润湿性方管中吸渗过程的格子玻尔兹曼模拟验证了所开发模型的准确性。然后，该模型被用于模拟强润湿性砂岩多孔介质中的吸渗过程，并分析了主流体弯月面和角落膜流动之间的竞争。结果表明，随着毛细数以及润湿性流体与非润湿性流体之间粘度比的降低，湿润角落膜的发展变得更加显著。
