-
本文改进了多弛豫时间多组分伪势格子玻尔兹曼方法, 使其可以适用于极大黏度比的情况. 研究了二维直通道内的两相流动, 分析了饱和度、黏度比、润湿性和毛细管数对相对渗透率和相渗曲线的影响. 结果表明, 大黏度相的相对渗透率随着黏度比的增大而增大, 维持高流动能力的饱和度范围也随之扩大, 小黏度相与之相反; 黏性耦合效应对两相相对渗透率的影响存在极限, 即不会导致小黏度相的相对渗透率低于最小值, 也不会使大黏度相的相对渗透率高于最大值; 非润湿条件会使大黏度相的流动能力增强, 相对渗透率增加, 并强化相渗曲线的非单调性; 润湿条件会减弱大黏度相的流动能力, 相对渗透率降低; 毛细管数的增加会导致更明显的接触角滞后效应, 使两相相对渗透率降低.Immiscible two-phase flow exists in many fields such as the chemical engineering, energy engineering, and geography. It is of great significance in enhancing the oil recovery and geological carbon dioxide storage. Relative permeability is a comprehensive reflection of the effects of various factors including the saturation, contact angle, viscosity ratio, and capillary number of immiscible two-phase flow, and it is a key parameter describing such a flow and thus can directly reflect the flow capacity of each phase. As a new fluid numerical calculation method, the lattice Boltzmann method (LBM) has its own unique advantage in simulating various flows because it does not need to capture the phase interface, and thus has been widely used in recent years. In this research, the multi-relaxation time multi-component pseudopotential lattice Boltzmann method is improved to make it possible to solve the two-phase relative permeability under high viscosity ratios. It is known that the LBM is a particle-based hydrodynamic method, with the spontaneous diffusion of components in multi-component system inevitably causing a mixing between components, leading to distortion of the velocity of two-phase flow with a high viscosity ratio. The color gradient model is introduced into the pseudopotential model for the treatment of mixing, the multiple relaxation time is employed to improve the numerical stability and viscosity ratio range, and the numerical distortion caused by mixing is therefore resolved. The resulting model is validated by the use of the two-phase layered flow in a 2D infinite channel, and accurate numerical results are obtained even if the viscosity ratio reaches 1000. Based on the improved model, the two-phase flow in a two-dimensional horizontal channel is analyzed and the effects of saturation, viscosity ratio, wettability and capillary number on relative permeability and phase permeation curves are investigated. The results show that the relative permeability of the high-viscosity phase increases with the increase of the viscosity ratio, and the saturation range for maintaining high flow capacity broadens with the increase of the low-viscosity phase. There is a limit for the influence of viscous coupling effect on the relative permeability of two phases: at a certain saturation, the relative permeability of the small viscosity phase has a minimum value while that of the large viscosity phase has a maximum value, the viscous coupling effect will never make the relative permeability of the small viscosity phase to be lower than the minimum value and that of the large viscosity phase to be higher than the maximum value. Non-wetting condition plays a role in increasing the flow capacity and relative permeability of the high-viscosity phase, while the wetting condition can weaken the flow capability, relative permeability of the high-viscosity phase. When the high-viscosity phase is a non-wetting phase, the decrease of the contact area with the solid wall increases its proportion in the middle flow passage, which enables it to maintain a strong flow capability in a wide saturation range. Wall wetting tends to inhibit fluid flow and consequently reduce the relative permeability. The increase of the capillary number may cause a more obvious contact angle hysteresis effect, which reduces the relative permeability of two phases.
-
Keywords:
- two-phase flow /
- relative permeability /
- viscosity ratio /
- lattice Boltzmann method
[1] Zhao H, Ning Z, Kang Q, Chen L, Zhao T 2017 Int. Commun. Heat Mass 85 53Google Scholar
[2] Zhang D, Papadikis K, Gu S 2016 Adv. Water Resour. 95 61Google Scholar
[3] Dullien F A L, Dong M 1996 Transport Porous Med. 25 97Google Scholar
[4] Bentsen R G, Manai A A 1993 Transport Porous Med. 11 243Google Scholar
[5] Akbarabadi M, Piri M 2013 Adv. Water Resour. 52 190Google Scholar
[6] Hao L, Cheng P 2010 Int. J. Heat Mass Tran. 53 1908Google Scholar
[7] Li Z, Chen T, Ning Y, Zhang K, Guan Q 2019 Qin SPE Middle East Oil and Gas Show and Conference, Bahrain, March 18–21, 2019 SPE-194770-MS
[8] Wang J, Dong M, Asghari K 2006 SPE/DOE Symposium on Improved Oil Recovery, USA, April 22–26, 2006 SPE-99763-MS
[9] Inamuro T, Ogata T, Tajima S, Konishi N 2004 J. Comput. Phys. 198 628Google Scholar
[10] Rothman D H, Keller J M 1988 J. Stat. Phys. 52 1119Google Scholar
[11] Shan X, Chen H 1993 Phys. Rev. E 47 1815Google Scholar
[12] Swift M R, Orlandini E, Osborn W R, Yeomans J M 1996 Phys. Rev. E 54 5041Google Scholar
[13] He X, Shan X, Doolen G D 1998 Phys. Rev. E 57 R13
[14] Ramstad T, Øren P E, Bakke S 2010 SPE J. 15 917Google Scholar
[15] Liu T, Wang M 2022 Transport Porous Med. 144 111Google Scholar
[16] Li H, Pan C, Miller C T 2005 Phys. Rev. E 72 026705Google Scholar
[17] Yiotis A G, Psihogios J, Kainourgiakis M E, Papaioannou A, Stubos A K 2007 Colloid. Surface. A 300 35Google Scholar
[18] Huang H, Li Z, Liu S, Lu X Y 2009 Int. J. Numer. Meth. Fl. 61 341Google Scholar
[19] Ghassemi A, Pak A 2011 J. Petrol Sci. Eng. 77 135Google Scholar
[20] Gao C, Xu R N, Jiang P X 2015 Int. J. Numer. Meth Fl. 25 1957
[21] Dou Z, Zhou Z F 2013 Int. J. Heat Fluid Fl. 42 23Google Scholar
[22] Zhang D, Li S, Jiao S, Shang Y, Dong M 2019 Int. J. Heat Mass Tran. 134 311Google Scholar
[23] Li S, Jiang F, Wei B, Hou J, Liu H 2021 Phys. Fluids 33 063302Google Scholar
[24] Huang H, Sukop M, Lu X 2015 Multiphase Lattice Boltzmann Methods: Theory and Application (Oxford: Wiley Blackwell) pp96–99
[25] Guo Z, Zheng C, Shi B 2008 Phys. Rev. E 77 036707Google Scholar
[26] Zheng J, Chen Z, Xie C, Wang Z, Lei Z, Ju Y, Wang M 2018 Comput. Fluids 168 21Google Scholar
-
图 1 两组分表面张力验证结果 (a), (b) 方形液滴在表面张力的作用下逐渐变为圆形液滴; (c) 液滴内外的压力差与曲率半径的线性关系
Fig. 1. The verification of two-component surface tension: (a), (b) The square drops gradually become round drops under the action of surface tension; (c) linear relationship between the pressure difference inside and outside the droplet and the radius of curvature.
图 7
$ {S}_{1}=0.7 $ 时, 不同黏度比下的流动形态(a1)—(a5)和速度分布(b1)—(b5) (a1), (b1) M = 1; (a2), (b2) M = 20; (a3), (b3) M = 40; (a4), (b4) M = 100; (a5), (b5) M = 200Fig. 7. Flow pattern (a1)–(a5) and velocity distribution (b1)–(b5) for the different viscosity ratios at
$ {S}_{1}=0.7 $ : (a1), (b1) M = 1; (a2), (b2) M = 20; (a3), (b3) M = 40; (a4), (b4) M = 100; (a5), (b5) M = 200.表 1 不同黏度比下的两组分表面张力
Table 1. Surface tensions of two-component under different viscosity ratios.
黏度比M 表面张力σ 1 0.1452 20 0.1379 50 0.1330 100 0.1371 200 0.1325 -
[1] Zhao H, Ning Z, Kang Q, Chen L, Zhao T 2017 Int. Commun. Heat Mass 85 53Google Scholar
[2] Zhang D, Papadikis K, Gu S 2016 Adv. Water Resour. 95 61Google Scholar
[3] Dullien F A L, Dong M 1996 Transport Porous Med. 25 97Google Scholar
[4] Bentsen R G, Manai A A 1993 Transport Porous Med. 11 243Google Scholar
[5] Akbarabadi M, Piri M 2013 Adv. Water Resour. 52 190Google Scholar
[6] Hao L, Cheng P 2010 Int. J. Heat Mass Tran. 53 1908Google Scholar
[7] Li Z, Chen T, Ning Y, Zhang K, Guan Q 2019 Qin SPE Middle East Oil and Gas Show and Conference, Bahrain, March 18–21, 2019 SPE-194770-MS
[8] Wang J, Dong M, Asghari K 2006 SPE/DOE Symposium on Improved Oil Recovery, USA, April 22–26, 2006 SPE-99763-MS
[9] Inamuro T, Ogata T, Tajima S, Konishi N 2004 J. Comput. Phys. 198 628Google Scholar
[10] Rothman D H, Keller J M 1988 J. Stat. Phys. 52 1119Google Scholar
[11] Shan X, Chen H 1993 Phys. Rev. E 47 1815Google Scholar
[12] Swift M R, Orlandini E, Osborn W R, Yeomans J M 1996 Phys. Rev. E 54 5041Google Scholar
[13] He X, Shan X, Doolen G D 1998 Phys. Rev. E 57 R13
[14] Ramstad T, Øren P E, Bakke S 2010 SPE J. 15 917Google Scholar
[15] Liu T, Wang M 2022 Transport Porous Med. 144 111Google Scholar
[16] Li H, Pan C, Miller C T 2005 Phys. Rev. E 72 026705Google Scholar
[17] Yiotis A G, Psihogios J, Kainourgiakis M E, Papaioannou A, Stubos A K 2007 Colloid. Surface. A 300 35Google Scholar
[18] Huang H, Li Z, Liu S, Lu X Y 2009 Int. J. Numer. Meth. Fl. 61 341Google Scholar
[19] Ghassemi A, Pak A 2011 J. Petrol Sci. Eng. 77 135Google Scholar
[20] Gao C, Xu R N, Jiang P X 2015 Int. J. Numer. Meth Fl. 25 1957
[21] Dou Z, Zhou Z F 2013 Int. J. Heat Fluid Fl. 42 23Google Scholar
[22] Zhang D, Li S, Jiao S, Shang Y, Dong M 2019 Int. J. Heat Mass Tran. 134 311Google Scholar
[23] Li S, Jiang F, Wei B, Hou J, Liu H 2021 Phys. Fluids 33 063302Google Scholar
[24] Huang H, Sukop M, Lu X 2015 Multiphase Lattice Boltzmann Methods: Theory and Application (Oxford: Wiley Blackwell) pp96–99
[25] Guo Z, Zheng C, Shi B 2008 Phys. Rev. E 77 036707Google Scholar
[26] Zheng J, Chen Z, Xie C, Wang Z, Lei Z, Ju Y, Wang M 2018 Comput. Fluids 168 21Google Scholar
计量
- 文章访问数: 3958
- PDF下载量: 106
- 被引次数: 0