引言
注水开发是油田保持地层压力、提高采收率的最有效且应用最广的方法, 但地层在经历长时间流体冲刷后易在注采井之间形成优势通道, 注入水并未能按照预期的驱替路径流动, 导致波及范围小、驱替效率低[1-2]. 问题的关键在于传统水驱主要利用压力场实现注入水在地下的渗流与波及, 而压力场调控难度大. 因此, 如何发展绿色智能驱油技术, 进一步提高原油采收率成为亟待解决的难题[3-4].
离子液体是一类可调控、多功能的绿色环保材料, 本质上是一种在室温或室温附近温度下呈现液态的盐, 又称为室温离子液体[5]. 与常见的盐一样, 离子液体仅由阳离子和阴离子组成, 通过静电相互作用结合在一起. 而与固态盐形成晶体结构不同, 离子液体由能够自由移动的离子组成, 如图1所示, 从而使盐具有流动的特性[6]. 离子液体具有良好的离子导电性与电化学稳定性, 主要应用于电化学领域[7-8], 同时离子液体具有多样性, 许多功能型离子液体正在研发中, 在生物医学[9]、绿色化学[10-11]以及光电子[12]等多个学科领域具有广阔的应用前景.
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/8//lxxb2021-156-1.jpg)
class="figure_img
figure_type1 bbb " id="Figure1" />
图
1
离子液体宏观与微观状态示意图
Figure
1.
Diagram of the macroscopic and microscopic states of ionic liquids
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
全尺寸图片
幻灯片
最近, 离子液体在提高采收率领域的应用研究受到****们关注. Dahbag等[13]在离子液体驱油实验的基础上, 利用CMG-STARS软件中的表面活性剂模型对离子液体驱油过程进行数值模拟, 实验接触角测量和模拟的相对渗透率曲线变化表明: 离子液体提高采收率机理主要是促进岩石润湿性向更亲水转变. Nabipour等[14]分别研究了地层盐水和海水对离子液体驱中界面张力和润湿性变化的影响, 结果表明: 离子液体不仅可以承受25 g/L的含盐度, 而且还能更有效的降低油水界面张力. Velusamy等[15]研究了6种常见咪唑型离子液体在高含盐和不含盐条件下的界面张力变化, 接触角测量结果表明: 在室温和高含盐条件(NaCl浓度为200 g/L)下离子液体能够使油水界面张力降低49.7%. Manshad等[16]通过悬滴法和接触角测试证明了离子液体可以作为一种新型表面活性剂来改变岩石润湿性和降低油水界面张力, 岩心驱替实验表明: 与水驱相比, 原油采收率提高了13%. Pillai等[17]研究了[C8mim]BF4, [C10mim]BF4, [C12mim]BF4这3种合成离子液体对界面张力和油湿岩石润湿性的影响, 实验发现这3种离子液体不仅能在高矿化度和高温条件下保持稳定, 还能有效提高岩石润湿性能并且显著降低界面张力. BinDahbag等[18]利用驱油实验的采收率和压降数据, 对岩心盐水驱和离子液体驱过程进行了拟合, 发现在高矿化度条件下采用离子液体驱获得的采收率比盐水驱更高. 可见, 离子液体在提高采收率方面的研究思路主要是将其视为传统表面活性剂的替代品, 仅仅关注了离子液体优异的化学特性. 然而离子液体同样表现出良好的物理场响应特征, 有望把材料的响应性和驱动性有机地结合为一体, 通过多物理场刺激来调控材料的运移路径, 从而实现离子液体智能驱油, 达到提高采收率的目的.
磁场作为一种有效的驱动和控制手段, 由于不需要与样品或试剂的直接接触即可实现对分析物的操控, 同时具有驱动力强、可控性高等优势, 被广泛用于生物医学、微流控等研究领域[19-20]. 在提高采收率领域, Kothari等[21]利用表面活性剂与磁性纳米颗粒复配制备了铁磁流体, 并总结了铁磁流体驱油机理: 表面活性剂的加入将大大降低界面张力, 同时铁磁流体具有与磁场相互作用的能力, 外加磁场可以提高驱替液的流动能力, 表现出比表面活性剂驱油更高的波及效率. 姚军等[22]开展了磁场?渗流场耦合作用下的铁磁流体多孔介质流动数值模拟, 研究了利用磁场控制的铁磁流体驱油问题, 结果表明: 利用磁场控制铁磁流体的驱替路径可以提高波及范围从而提高采收率, 为智能驱油提供了新的思路. 所使用的铁磁流体是一类特殊的智能纳米材料, 不仅表现出一般流体的性质, 同时可以受到磁场的控制[23]. 由于离子液体功能的可定制性, 利用特殊的阳离子和阴离子可以制备与铁磁流体一样带有磁性的离子液体. 磁性离子液体是一类具有磁性的功能化离子液体, 不仅保留了一般离子液体的蒸气压低、液程宽、良好的溶解能力等优点, 而且能够在磁场作用下产生一定的磁化强度, 并对外界磁场产生宏观响应[24]. Hayashi和Hamaguchi[25]在2004年首次提出了磁性离子液体的概念, 他们利用[bmim]Cl和FeCl3为原料合成了一种离子液体([bmim]FeCl4), 该离子液体能够被磁铁吸引, 表现出很强的顺磁性. 图2中的照片显示了磁性离子液体([bmim]FeCl4)在磁场下的响应, 无磁铁作用时(图2(a)), 离子液体层在水层下扩散, 当磁铁靠近时, 离子液体层被磁铁吸引(图2(b)), 最后在磁场力的作用下被扭曲形成凹曲线型(图2(c)).
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/8//lxxb2021-156-2.jpg)
class="figure_img
figure_type1 bbb " id="Figure2" />
图
2
离子液体([bmim]FeCl4)的磁场响应[25]
Figure
2.
Magnetic field response of ionic liquid [bmim]FeCl4[25]
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
全尺寸图片
幻灯片
虽然磁场对磁性流体具备很好的控制能力, 但磁性离子液体是一种特殊的离子液体, 大多数离子液体并不具备磁场响应, 因此磁场会限制离子液体的应用范围. 相比磁场而言, 离子液体的电磁场响应更加值得关注, 有望应用于调控水驱油路径, 研究的关键在于揭示多孔介质中离子液体在电磁场作用下的流动机理. 目前尚未有关于电磁场作用下离子液体在多孔介质中的相关研究. 本文在分析离子液体在毛细管中电磁场响应机理的基础上, 建立了电磁场?渗流场耦合作用下离子液体多孔介质流动模型, 通过理论推导与数值分析, 研究电磁场?渗流场耦合作用下离子液体在多孔介质中的流动规律, 并从离子液体优选、离子液体智能驱油、电致热量影响3个方面进行讨论.
1.
电磁场?渗流场耦合流动模型
1.1
毛细管中的洛伦兹力
在离子液体中有3种离子存在形态: 离子对(呈电中性)、离子团(数量很少)以及能够自由移动的阴、阳离子, 正是由于这些自由移动的阴、阳离子存在, 离子液体表现出优异的导电性能, 同时自由移动的阴、阳离子也是离子液体对电磁场产生响应的根本原因[26]. 下面分析离子液体在毛细管中的电磁场响应机理, 首先建立离子液体在电磁场作用下的毛细管模型, 如图3所示, 毛细管处在相互垂直的均匀电场(
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M1.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M3.jpg)
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/8//lxxb2021-156-3.jpg)
class="figure_img
figure_type1 bbb " id="Figure3" />
图
3
离子液体在电磁场作用下的毛细管模型
Figure
3.
Capillary model of ionic liquids under electromagnetic field
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
全尺寸图片
幻灯片
在电磁场作用下, 毛细管模型中任意一个自由离子会同时受到电场力和洛伦兹力的作用, 其受力情况可用如下公式表示[27]
$${{boldsymbol{F}}_i} = {n_i}e{boldsymbol{E}} + {n_i}e{{boldsymbol{v}}_i} times {boldsymbol{B}}$$ ![]() | (1) |
式中:
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M4.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M10.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M7.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M5.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M6.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M10.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M7.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M5.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M8.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M9.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M10.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M7.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M5.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M11.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M12.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M13.jpg)
由于离子液体内部阴、阳离子所受电场力的方向相反, 并且离子液体整体呈电中性, 则式(1)右边第一项对离子液体内部所有阴、阳离子求和后等于零, 即宏观上离子液体整体所受电场力为零. 式(1)右边第二项为运动的离子在磁场作用下产生的洛伦兹力, 其大小与在垂直于磁场的平面内的速度大小成正比, 离子速度在
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M14.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M15.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M16.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M17.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M18.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M19.jpg)
离子液体中不同速度分量的离子在宏观上受到的总洛伦兹力也不同. 沿毛细管轴向移动而具有
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M21.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M20.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M21.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M20.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M22.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M23.jpg)
$${{boldsymbol{v}}_{iy pm }} = {mu _ pm }{boldsymbol{E}}$$ ![]() | (2) |
式中
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M24.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M25.jpg)
如图3所示, 在毛细管模型中取一单位体积模型, 假设该单位体积模型中离散的阴、阳离子的粒子数密度分别为
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M26.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M27.jpg)
m{V}}}$
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M28.jpg)
$${{boldsymbol{F}}_{ m{V}}}{ m{ = }}sumlimits_i {{n_i}e{{boldsymbol{v}}_{iy}} times {boldsymbol{B}}} = {boldsymbol{E}} times {boldsymbol{B}}left( {{mu _ - }{m_ - }{n_ - }e + {mu _ + }{m_ + }{n_ + }e} ight)$$ ![]() | (3) |
由于离子液体中的离子团很少, 根据电中性原则, 阴、阳离子所带的电荷量
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M29.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M30.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M31.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M32.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M33.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M33-1630562005284.jpg)
$${{boldsymbol{F}}_{ m{V}}}{ m{ = }}{boldsymbol{E}} times {boldsymbol{B}}({mu _ - } + {mu _ + }){{lambda mq} / 2}$$ ![]() | (4) |
假设离子液体的密度为
ho $
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M34.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M35.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M36.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M37.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M38.jpg)
$$sigma { m{ = }}({mu _ - } + {mu _ + }){{lambda mq} / 2}$$ ![]() | (5) |
则单位质量的离子液体总洛伦兹力
m{L}}}$
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M39.jpg)
$${{boldsymbol{F}}_{ m{L}}}{ m{ = }}{{sigma {boldsymbol{E}} times {boldsymbol{B}}} / ho }$$ ![]() | (6) |
1.2
离子液体运动方程
下面从理论上对电磁场?渗流场耦合作用下离子液体在毛细管中的运动方程进行推导, 假设毛细管内离子液体的密度为
ho $
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M40.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M41.jpg)
$$frac{{partial {boldsymbol{v}}}}{{partial t}} + {boldsymbol{v}} cdot nabla {boldsymbol{v}} = {boldsymbol{f}} - frac{1}{ ho }nabla p + frac{mu }{ ho }{nabla ^2}{boldsymbol{v}}$$ ![]() | (7) |
式中,
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M42.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M43.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M44.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M45.jpg)
$$frac{{partial {boldsymbol{v}}}}{{partial t}} + {boldsymbol{v}} cdot nabla {boldsymbol{v}} = frac{1}{ ho }sigma {boldsymbol{E}} times {boldsymbol{B}} + {boldsymbol{g}} - frac{1}{ ho }nabla p + frac{mu }{ ho }{nabla ^2}{boldsymbol{v}}$$ ![]() | (8) |
假设毛细管直径均匀, 则离子液体在毛细管中流速只有沿轴向的分量, 这样运动方程中的对流项
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M46.jpg)
$$frac{{partial {boldsymbol{v}}}}{{partial t}} = frac{1}{ ho }sigma {boldsymbol{E}} times {boldsymbol{B}}{ m{ + }}{boldsymbol{g}} - frac{1}{ ho }nabla p + frac{mu }{ ho }{nabla ^2}{boldsymbol{v}}$$ ![]() | (9) |
1.3
多孔介质一维流动模型
为了研究电磁场?渗流场耦合作用下离子液体在多孔介质中的一维流动规律, 建立基于“直毛细管”的多孔介质一维流动模型, 如图4所示, 该模型由一束平行的等径毛细管所构成[30].
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/8//lxxb2021-156-4.jpg)
class="figure_img
figure_type1 bbb " id="Figure4" />
图
4
基于“直毛细管”的多孔介质一维流动模型
Figure
4.
One-dimensional flow model of porous media based on “straight capillary”
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
全尺寸图片
幻灯片
假设离子液体沿
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M47.jpg)
$$frac{{partial v}}{{partial t}} = frac{1}{ ho }sigma EB - frac{1}{ ho }frac{{partial p}}{{partial x}} + frac{mu }{ ho }{nabla ^2}v$$ ![]() | (10) |
模型求解选用圆柱坐标系
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M48.jpg)
$$frac{{partial v}}{{partial t}} + frac{1}{ ho }frac{{partial p}}{{partial x}} - frac{1}{ ho }sigma EB = frac{mu }{ ho }left[ {frac{1}{r}frac{partial }{{partial r}}left( {rfrac{{partial v}}{{partial r}}} ight){ m{ + }}frac{1}{{{r^2}}}frac{{{partial ^2}v}}{{partial {theta ^2}}} + frac{{{partial ^2}v}}{{partial {x^2}}}} ight]$$ ![]() | (11) |
当毛细管两端压差恒定且外加电磁场稳定时, 离子液体做定常流动, 有
m{ = }}0$
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M49.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M50.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M51.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M52.jpg)
$$frac{{{ m{d}}p}}{{{ m{d}}x}} - sigma EB = mu left( {frac{{{partial ^2}v}}{{partial {r^2}}} + frac{1}{r}frac{{partial v}}{{partial r}}} ight)$$ ![]() | (12) |
设毛细管半径为
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M53.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M54.jpg)
$$v = frac{1}{{4mu }}left( {sigma EB - frac{{{ m{d}}p}}{{{ m{d}}x}}} ight)left( {r_0^2 - {r^2}} ight)$$ ![]() | (13) |
通过厚度为
m{d}}r$
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M56.jpg)
m{d}}Q = v{
m{d}}A$
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M57.jpg)
m{d}}A$
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M58.jpg)
$$Q = intlimits_A {v{ m{d}}A} = int_0^{{r_0}} {v cdot 2{text{π}} r{ m{d}}r} $$ ![]() | (14) |
最后将式(13)代入式(14)积分可得电磁场?渗流场耦合作用下离子液体在毛细管中Poiseuille流动的流量方程
$$Q = frac{{{text{π}} r_0^4}}{{8mu }}left( {sigma EB - frac{{{ m{d}}p}}{{{ m{d}}x}}} ight)$$ ![]() | (15) |
1.4
多孔介质三维流动模型
图4所示的一维流动模型只能描述离子液体在一个方向上的流动规律, 为了弥补这一缺陷将模型作如下修正: 在3个互相正交的方向上各放置三分之一的毛细管[30], 建立基于“毛细管组”的多孔介质三维流动模型, 如图5所示.
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/8//lxxb2021-156-5.jpg)
class="figure_img
figure_type1 bbb " id="Figure5" />
图
5
基于“毛细管组”的多孔介质三维流动模型
Figure
5.
Three-dimensional flow model of porous media based on “capillary group”
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
全尺寸图片
幻灯片
假设离子液体可以在
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M59.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M60.jpg)
m{L}}}$
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M63.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M64.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M65.jpg)
$$left. {begin{array}{*{20}{l}} {dfrac{{partial {v_x}}}{{partial t}} = dfrac{1}{ ho }sigma EBcos alpha sin beta - dfrac{1}{ ho }dfrac{{partial p}}{{partial x}} + dfrac{mu }{ ho }{nabla ^2}{v_x}} {dfrac{{partial {v_y}}}{{partial t}} = dfrac{1}{ ho }sigma EBsin alpha sin beta - dfrac{1}{ ho }dfrac{{partial p}}{{partial y}} + dfrac{mu }{ ho }{nabla ^2}{v_y}} {dfrac{{partial {v_{textit{z}}}}}{{partial t}} = dfrac{1}{ ho }sigma EBcos beta - g - dfrac{1}{ ho }dfrac{{partial p}}{{partial {textit{z}}}} + dfrac{mu }{ ho }{nabla ^2}{v_{textit{z}}}} end{array}} ight}$$ ![]() | (16) |
类似一维模型, 在圆柱坐标系下采用无滑移边界求解并积分后可得离子液体3个方向上的体积流量
$$left. {begin{array}{*{20}{c}} {{Q_x} = dfrac{{{text{π}} r_0^4}}{{8mu }}left( {sigma EBcos alpha sin beta - dfrac{{partial p}}{{partial x}}} ight)} {{Q_y} = dfrac{{{text{π}} r_0^4}}{{8mu }}left( {sigma EBsin alpha sin beta - dfrac{{partial p}}{{partial y}}} ight)} {{Q_{textit{z}}} = dfrac{{{text{π}} r_0^4}}{{8mu }}left( {sigma EBcos beta - dfrac{{partial p}}{{partial {textit{z}}}} - ho g} ight)} end{array}} ight}$$ ![]() | (17) |
由于
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M66.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M67.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M68.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M69.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M70.jpg)
$$left. {begin{array}{*{20}{l}} {{{bar v}_x} = dfrac{{r_0^2}}{{8mu }}left( {sigma EBcos alpha sin beta - dfrac{{partial p}}{{partial x}}} ight)} {{{bar v}_y} = dfrac{{r_0^2}}{{8mu }}left( {sigma EBsin alpha sin beta - dfrac{{partial p}}{{partial y}}} ight)} {{{bar v}_{textit{z}}} = dfrac{{r_0^2}}{{8mu }}left( {sigma EBcos beta - dfrac{{partial p}}{{partial {textit{z}}}} - ho g} ight)} end{array}} ight}$$ ![]() | (18) |
2.
数值分析
下面以疏松砂岩油藏经长期注水开发后形成的高渗通道为研究对象, 对电磁场?渗流场耦合作用下离子液体在多孔介质中的流动规律进行数值分析. 在数值计算过程中, 根据高渗通道孔隙半径分布特点, 将多孔介质流动模型中的毛细管半径
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M71.jpg)
2.1
洛伦兹力
离子液体在电磁场作用下产生响应的本质是洛伦兹力的存在. 图6为单位体积离子液体的洛伦兹力与电磁场强度的三维关系图, 其中离子液体的电导率为0.5 S/m. 从图中可以看出, 当离子液体的电导率一定时, 单位体积离子液体的洛伦兹力与电磁场强度有关; 当磁感应强度一定时, 洛伦兹力随着电场强度的增大而增大, 当电场强度一定时, 洛伦兹力随着磁感应强度的增大而增大, 即电场强度与磁感应强度的乘积越大, 单位体积离子液体的洛伦兹力越大.
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/8//lxxb2021-156-6.jpg)
class="figure_img
figure_type1 bbb " id="Figure6" />
图
6
单位体积离子液体的洛伦兹力与电磁场强度的关系
Figure
6.
The relationship between Lorentz force per unit volume of ionic liquids and electromagnetic field strength
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
全尺寸图片
幻灯片
2.2
毛细管流量
图7为毛细管流量与离子液体电导率、黏度的三维关系图, 其中电场强度为5 × 104 V/m, 磁感应强度为0.5 T, 压力梯度为5.0 kPa/m. 从图中可以看出, 当电磁场强度和压力梯度一定时, 毛细管流量与离子液体的电导率与黏度有关; 当黏度一定时, 毛细管流量随着离子液体电导率的增大而增大, 当电导率一定时, 毛细管流量随着离子液体黏度的增大而减小, 即离子液体的电导率与黏度的比值越大, 离子液体的毛细管流量也越大.
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/8//lxxb2021-156-7.jpg)
class="figure_img
figure_type1 bbb " id="Figure7" />
图
7
毛细管流量与离子液体电导率、黏度的关系
Figure
7.
The relationship between capillary flow rate and conductivity and viscosity of ionic liquid
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
全尺寸图片
幻灯片
图8为离子液体毛细管流量与电磁场强度的三维关系图, 其中电导率为0.5 S/m, 黏度为100 mPa·s, 压力梯度为5.0 kPa/m. 从图中可以看出, 当离子液体电导率、黏度和压力梯度一定时, 离子液体的毛细管流量只与电磁场强度有关, 当磁感应强度一定时, 离子液体的毛细管流量随着电场强度的增大而增大, 当电场强度一定时, 离子液体的毛细管流量随着磁感应强度的增大而增大, 即电场强度与磁感应强度的乘积越大, 离子液体的毛细管流量也越大.
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/8//lxxb2021-156-8.jpg)
class="figure_img
figure_type1 bbb " id="Figure8" />
图
8
离子液体的毛细管流量与电磁场强度的关系
Figure
8.
The relationship between capillary flow rate of ionic liquids and electromagnetic field strength
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
全尺寸图片
幻灯片
下面用电场强度与磁感应强度的乘积来表征电磁场强度, 即
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M73.jpg)
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/8//lxxb2021-156-9.jpg)
class="figure_img
figure_type1 bbb " id="Figure9" />
图
9
毛细管流量与压力梯度、电磁场强度的三维关系
Figure
9.
Three-dimensional relationship between capillary flow rate, pressure gradient and electromagnetic field strength
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
全尺寸图片
幻灯片
2.3
平均线性流速
下面利用三维流动模型分析电磁场?渗流场耦合作用下离子液体在多孔介质中的流动特征. 图10(a)和图10(b)分别为离子液体在
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M78.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M74.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M79.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M75.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M76.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M77.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M78.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M74.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M79.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M75.jpg)
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/8//lxxb2021-156-10.jpg)
class="figure_img
figure_type1 bbb " id="Figure10" />
图
10
离子液体在(a) x方向和(b) y方向平均线性流速与洛伦兹力方向参数的三维关系
Figure
10.
Three-dimensional relationship between the average linear velocity of ionic liquids in the (a) x-direction and (b) y-direction and the direction parameter of Lorentz force
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
全尺寸图片
幻灯片
由于离子液体在
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M86.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M84.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M83.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M80.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M82.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M86.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M84.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M83.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M80.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M86.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M84.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M83.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M80.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M87.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M85.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M86.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M84.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M83.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M80.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M87.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M85.jpg)
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/8//lxxb2021-156-11.jpg)
class="figure_img
figure_type1 bbb " id="Figure11" />
图
11
离子液体在z方向的平均线性流速与洛伦兹力方向参数、密度的三维关系
Figure
11.
Three-dimensional relationship between the average linear velocity of ionic liquids in the z-direction, the direction parameter of Lorentz force and density
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
全尺寸图片
幻灯片
3.
讨论
3.1
离子液体优选
根据文献[31]中的公式, 可以求出5种离子液体在温度为25 °C时的电导率和黏度数据, 并计算相应电导率与黏度的比值, 如表1所示.
下面利用基于“直毛细管”简化的多孔介质一维流动模型来分析表1中4种离子液体在电磁场?渗流场耦合作用下的流动特点. 图12为不同离子液体的毛细管流量与压力梯度的关系, 其中电磁场强度为2.5 × 104 V/m·T. 由图中分析可知, [bmim][(CF3SO2)2N]的毛细管流量明显高于其他离子液体, 主要是因为其电导率与黏度的比值最大; 而[bmim][BF4]的电导率与黏度的比值最小, 因此其毛细管流量最小; [bmim][CF3SO3]与[bmim][BF4]两种离子液体对应的直线几乎重合, 两者的毛细管流量相近, 结合表1中数据分析可知, 主要是因为这两种离子液体的电导率与黏度的比值相当. 由局部放大图可知, 在较小压力梯度时, 因[bmim][BF4]的电导率与黏度比值略大于[bmim][CF3SO3], 其毛细管流量要大于[bmim][CF3SO3]; 随着压力梯度不断增大, 由于[bmim][CF3SO3]的黏度比[bmim][BF4]小, 其毛细管流量逐渐高于[bmim][BF4]. 因此, 在优选离子液体时, 应选择电导率与黏度的比值较大的离子液体作为流体介质, 若离子液体电导率与黏度的比值相近时, 应选择黏度较小的离子液体, 从而保证离子液体在多孔介质中更容易流动.
表
1
五种离子液体的电导率和黏度
Table
1.
Conductivity and viscosity of five ionic liquids
table_type2 ">
Ionic liquids | [bmim][(C2F5SO2)2N] | [bmim][(CF3SO2)2N] | [bmim][CF3SO3] | [bmim][PF6] | [bmim][BF4] |
$sigma $/(S·m?1) | 0.14 | 0.39 | 0.29 | 0.15 | 0.35 |
$mu $/(mPa·s) | 108.40 | 49.58 | 83.64 | 258.74 | 99.42 |
${sigma cdot mu^{-1} }$/(S·m?1·Pa?1·s?1) | 1.29 | 7.87 | 3.47 | 0.58 | 3.52 |
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
导出CSV
|显示表格
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/8//lxxb2021-156-12.jpg)
class="figure_img
figure_type1 bbb " id="Figure12" />
图
12
不同离子液体的毛细管流量与压力梯度的关系
Figure
12.
The relationship between capillary flow and pressure gradient of different ionic liquids
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
全尺寸图片
幻灯片
3.2
离子液体智能驱油
离子液体具有良好的电磁场响应, 电磁场的可调控性为离子液体的智能驱油提供了思路. 基于“毛细管组”简化的多孔介质三维流动模型, 从理论上说明了通过改变角度参数
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M91.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M92.jpg)
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/8//lxxb2021-156-13.jpg)
class="figure_img
figure_type1 bbb " id="Figure13" />
图
13
在平面上存在优势通道时(a)传统水驱和(b)离子液体驱的流动路径示意
Figure
13.
Diagram of the flow paths of (a) traditional waterflooding and (b) ionic liquids flooding when there is a dominant channel laterally
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
全尺寸图片
幻灯片
图14为在纵向上存在高渗层时传统水驱和离子液体驱的流动路径示意图, 由于高渗层的存在以及重力作用, 传统水驱时注入水在压力差的驱动下主要沿着高渗层流动(图14(a)), 进入低渗层中的注入水较少; 离子液体驱时, 由于竖直向上的洛伦兹力的存在(14(b)), 原本进入高渗层的离子液体会向低渗层中流动, 进而改变离子液体驱替方向. 可见, 无论是平面上存在优势通道, 还是纵向上存在高渗层, 通过电磁场的调控可以改变洛伦兹力的方向, 从而控制离子液体朝低波及区域驱替, 达到提高波及范围而提高采出程度目的.
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/8//lxxb2021-156-14.jpg)
class="figure_img
figure_type1 bbb " id="Figure14" />
图
14
在纵向上存在高渗层时(a)传统水驱和(b)离子液体驱的流动路径示意
Figure
14.
Diagram of the flow paths of (a) traditional waterflooding and (b) ionic liquids flooding when there is a high permeability layer vertically
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
全尺寸图片
幻灯片
3.3
电致热量影响
离子液体在电磁场作用下会产生电流, 电流会产生热量即电致热量, 随着热量不断增加, 电磁场产生的热量会使离子液体的温度升高, 而离子液体是温度敏感材料, 其电导率和黏度也会随着温度的变化发生相应的改变, 从而影响电磁场?渗流场耦合作用下离子液体在多孔介质中的流动. 根据文献[31]中的公式, 可以求出[bmim][PF6]和[bmim][BF4]这两种常见离子液体在不同温度
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M93.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M94.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M95.jpg)
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/lxxb2021-156_M96.jpg)
表
2
离子液体[bmim][PF6]在不同温度下的电导率和黏度
Table
2.
Conductivity and viscosity of ionic liquid [bmim][PF6] at different temperatures
table_type1 ">
T/°C | 25 | 35 | 45 | 55 | 65 |
$sigma $/(S·m?1) | 0.15 | 0.25 | 0.39 | 0.57 | 0.80 |
$mu $/(mPa·s) | 258.74 | 140.04 | 84.17 | 54.81 | 38.00 |
${sigma cdot mu^{-1} }$/(S·m?1·Pa?1·s?1) | 0.58 | 1.79 | 4.63 | 10.40 | 21.05 |
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
导出CSV
|显示表格
表
3
离子液体[bmim][BF4]在不同温度下的电导率和黏度
Table
3.
Conductivity and viscosity of ionic liquid [bmim][BF4] at different temperatures
table_type1 ">
T/°C | 25 | 35 | 45 | 55 | 65 |
$sigma $/(S·m?1) | 0.35 | 0.55 | 0.81 | 1.13 | 1.51 |
$mu $/(mPa·s) | 99.42 | 60.29 | 39.41 | 27.34 | 19.89 |
${sigma cdot mu^{-1} }$/(S·m?1·Pa?1·s?1) | 3.52 | 9.12 | 20.55 | 41.33 | 75.92 |
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
导出CSV
|显示表格
![](https://lxxb.cstam.org.cn/fileLXXB/journal/article/lxxb/2021/8//lxxb2021-156-15.jpg)
class="figure_img
figure_type1 bbb " id="Figure15" />
图
15
两种常见离子液体的毛细管流量与温度的关系
Figure
15.
The relationship between capillary flow and temperature of two common ionic liquids
![](https://lxxb.cstam.org.cn/cn/article/doi/10.6052/style/web/images/public/download.png)
全尺寸图片
幻灯片
4.
结论与认识
在分析离子液体在毛细管中电磁场响应机理的基础上, 建立了电磁场?渗流场耦合作用下离子液体多孔介质流动模型. 通过理论推导与数值分析得到以下结论与认识:
(1)电磁场?渗流场耦合作用下毛细管流量大小主要由离子液体电导率与黏度的比值(内因)、电磁场强度与压力梯度(外因)两方面因素决定;
(2)电磁场通过洛伦兹力在离子液体上产生一个电磁驱动压强, 形成一个类似压力梯度的电磁驱动等效压力梯度, 从而改变离子液体的流量, 当电磁场强度为2.0 × 104 V/m·T时, 电磁场在电导率为0.5 S/m的离子液体上可形成10 kPa/m电磁驱动等效压力梯度;
(3)通过调整电磁场方向即可控制离子液体在多孔介质中的流动方向, 解决利用压力差难以控制流动路径的难题, 为离子液体智能驱油提供理论依据;
(4)离子液体在电磁场作用产生的电致热量, 不仅会影响离子液体的电导率和黏度, 也会降低原油黏度、改变岩石润湿性. 因此, 电磁场产生的热效应对离子液体的流动能力及驱油效率的影响有待进一步探索.