高频电磁加热稠油储层温度分布及其影响因素分析

2020-02-21王正旭高德利

石油钻探技术 2020年1期

王正旭，高德利 2

（1. 石油工程教育部重点实验室（中国石油大学（北京）），北京 102249；2. 油气资源与工程国家重点实验室（中国石油大学（北京）），北京 102249）

随着常规油气资源被持续开发、剩余资源日益减少，国民经济发展对油气需求不断增长，以及勘探开发技术不断提高，地质储量丰富的稠油资源将是未来的接替能源之一，具有广阔的开发利用前景。目前采用的稠油蒸汽热采技术存在诸多缺陷：蒸汽在管线及非储层区域热损失大；深层、薄层及裂缝性储层的原油动用难度大[1-3]；地面蒸汽生产设备消耗大量天然气并导致温室气体过度排放，最终导致开采成本过高和环境污染严重[4-5]。为此，国内外开展了诸多稠油开采新技术研究[6-7]，电磁加热法是其中一种有效的稠油开采技术[8-10]。

1956 年，H. W. Ritchey[11]提出了电磁加热法稠油开采技术；随后，多位学者对该技术进行了广泛的理论研究[12-16]。近年来，因高频电磁加热具有加热快和可以体积加热、内外同时加热的特点，也符合环保要求[17-18]，基于高频波段（大于100 kHz[19]）的电磁加热法得到了高度关注[20-24]。一些学者利用数值模拟方法分析了高频电磁加热过程中损耗多孔媒质的热传递过程，比较了电磁加热和其他加热法的加热性能，并计算了电磁加热用于超稠油开采的温度分布结果[25-28]，但忽略了电磁加热过程中储层性质一直处于动态变化的情况[2]，导致计算结果出现偏差，不能准确反映电磁加热过程。因此，为了准确计算电磁加热过程中稠油储层的温度分布，笔者采用喇叭天线来辐射高频电磁波，建立了考虑稠油储层性质变化的温度分布数学模型，通过求解该模型来分析稠油储层温度分布的影响因素，以期为高频电磁加热技术的现场应用提供理论依据。

1 电磁加热机理

稠油储层是损耗媒质，电磁波在其内传播时会产生电磁损耗，最终电磁能在储层内转化为热能，实现原油增温降黏的目的。稠油储层内主要含有油、水和岩石（见图1），由于各组分存在性质差异，它们对电磁波的吸收能力也有所不同。

一般情况下，水的吸收能力较强，这是由于水分子是极性分子，置于外加电场后，每个水分子的正负电荷受到外电场力的作用，形成偶极子；当外电场按照一定频率发生变化时，偶极子被反复极化，分子的剧烈运动引起分子间摩擦并最终产生热能[19]。因此，当电磁波穿过储层内部时，水吸收较多电磁能，并在储层内部通过热传导作用加热其周围的油和岩石。

2 电磁加热模型

2.1 物理模型

为研究高频电磁加热过程中稠油储层温度分布的影响因素，建立了三维几何模型，如图2 所示。该模型包括喇叭天线、腔体及稠油储层。其中，喇叭天线的作用是辐射一定频率和功率的电磁波，腔体将电磁波从喇叭天线传输至稠油储层。图2 中，喇叭天线尺寸：La=200.0 mm，Lh=54.6 mm，Wh=109.0 mm，Wa=200.0 mm，Hh=50.0 mm 和Ha=300.0 mm；腔体尺寸：Lc=388.0 mm，Wc=470.0 mm，Hc=200.0 mm；稠油储层尺寸：Lr=388.0 mm，Wr=470.0 mm，Hr=450.0 mm。坐标原点位于腔体和储层交界面的中心位置，喇叭天线及腔体内部充满空气，并假设其壁面为理想导体，稠油储层均匀分布且具有各向同性。

图 2 电磁加热稠油储层三维几何模型Fig.2 Three-dimensional geometric model of a heavy oil reservoir under electromagnetic heating

图 1 电磁加热稠油储层示意Fig.1 Schematic of heavy oil reservoir under electromagnetic heating

2.2 储层性质

高频电磁加热过程中，储层的性质会发生动态变化，为此，分别分析了储层电导率、相对介电常数随电磁波频率的变化规律和导热系数、比热容随温度变化的特性。

电磁波频率对储层电导率和相对介电常数的影响规律如图3 所示[29]。

图 3 电磁波频率对储层电导率和相对介电常数的影响Fig.3 Changing laws of electric conductivity and relative permittivity of reservoir with frequency

通过最小二乘拟合法，可以得到储层电导率和相对介电常数关于电磁波频率的函数表达式，其相对误差分别为3.3%和4.2%。

储层电导率的计算公式为：

式中：σ 为储层的电导率，S/m；f 为电磁波频率，Hz。相对介电常数的计算公式为：

式中：εr为储层的相对介电常数。

储层导热系数和比热容与温度的关系可表示为[2]：

式中：k 为导热系数，W/（m·℃）；T 为温度，℃；c 为比热容，J/（kg·℃）。

此外，储层密度ρ 为恒定值1 950 kg/m3，且初始温度T0设为25 ℃。

2.3 电磁场与电磁热源

电磁波在损耗媒质中传播时，呈现逐渐衰减的特点（见图4），其中E 和H 分别为电场强度和磁场强度，且电磁波沿z1方向传播[25]。

因此，根据损耗媒质中的麦克斯韦方程[30]，可推导出该情况下的齐次亥姆霍兹方程[31]：

图 4 电磁波衰减示意Fig. 4 Schematic of electromagnetic waves attenuation

式中：E 为媒质中电场强度，V/m；kc为损耗媒质中的波数，rad/m。

kc为复数，也称为复数波，其表达式为[31]：

式中：ω 为角频率，rad/s；μ 为损耗媒质的磁导率，H/m；εc为复介电常数。

其中，复介电常数和角频率的表达式分别为：

式中：ε 为损耗媒质的介电常数，F/m；j 为虚数单位；σ 为媒质的电导率，S/m。

根据式（6）和式（8），可知真空中的波数k0为：

式中：ε0为真空介电常数，F/m；μ0为真空磁导率，H/m。将式（6）—式（9）代入式（5），得：

式中：μr为相对磁导率，即媒质磁导率μ 与真空磁导率μ0之比，由于假定储层是无磁的，μr=1.0；εr为相对介电常数，即损耗媒质介电常数ε 与真空介电常数ε0之比。

此外，对于电磁强度E，有如下关系：

在无源的麦克斯韦表达式中， ∇·E为0，再将式（11）代入式（10），得到电场在储层内传播时的计算表达式：

式（12）中，储层的相对介电常数和电导率是关于电磁波频率的函数。

在储层内，电磁波单位体积损耗功率的表达式为[2]：

式中：Q 为电磁波单位体积的损耗功率，W/m3。

式中：t 为加热时间，s；ρ 为储层密度，kg/m3；T 为储层被加热时间t 后的温度，℃；q 为内热源强度，W/m3，表示单位体积单位时间内所吸收的热量。

利用式（12）求得储层内电场，然后代入式（13），计算得到的Q 项等于式（14）中的内热源强度q[32]，即Q=q。其中，储层不同温度下的导热系数和比热容分别利用式（3）和式（4）进行计算。最后，求解式（14）可得到储层的温度分布。

初始温度设为恒定值，即：

此外，稠油储层外表面均是绝热的，根据傅里叶定律，该表面上温度梯度为0，边界条件表示为：

式中：s 为储层外表面；n 为s 面的法向量。

采用多物理场模拟软件COMSOL 求解上述数学模型，为了提高计算精度，对物理模型交界区域进行了网格加密处理（见图5）。

图 5 电磁加热稠油储层物理网格划分Fig.5 Grid partition of the heavy oil reservoir under electromagnetic heating

3 电场和温度场分布

电磁波频率为1.6 GHz、功率为1.5 kW 时，稠油储层加热6 h 后，模拟喇叭天线和稠油储层的电场分布，结果见图6。从图6 可知：电场强度在喇叭天线端口处较大且电磁波沿着z 轴方向向储层传播；在电磁波入射的储层面上，与y 轴平行的两边附近区域的电场强度较大，且电场在电磁波入射的储层面内沿y 轴对称分布。

图 6 电场分布模拟结果Fig.6 Calculation results of electric field distribution

储层温度分布随时间的变化情况如图7 所示。从图7 可以看出：电磁波入射面的储层温度也呈现沿y 轴对称分布的特点，这与图6（b）所示的电场分布规律类似；加热初期，在电磁波入射的储层面内，平行于y 轴的两边附近区域的温度较高，随着加热时间增长，电磁波入射储层面中心区域的温度逐渐升高，接近上述两边附近的温度，且沿加热深度逐步扩大。这是由于在热传导的作用下，加热时间延长，导致加热范围逐步扩大。

4 影响因素分析

高频电磁加热过程中，稠油储层温度分布的影响因素主要包括电磁波功率和频率，储层的导热系数、比热容、相对介电常数和电导率。

图 7 储层温度随时间的变化规律Fig.7 Changing relationship of reservoir temperature with time

图 8 不同电磁波功率下储层温度分布三维图Fig.8 Three-dimensional map of reservoir temperature distribution under different electromagnetic wave powers

4.1 电磁波功率和频率

在加热时间恒为6 h 时，电磁波功率对储层温度分布的影响见图8。从图8 可以看出，随着电磁波功率增大，储层加热深度增大，电磁波入射储层面的温度分布也趋于均匀。由此可知，电磁波功率增大，能使更多的电磁能转化为储层内的热能，储层的高温区域增大。

电磁波频率对储层温度分布的影响研究分2 种情况进行：1）基于实际情况，考虑频率变化引起储层电导率和相对介电常数的变化，其变化规律符合式（1）和式（2）；2）储层电导率和相对介电常数不受频率影响，始终为0.08 S/m 和12.5。

根据以上2 种情况，一方面分析电磁波频率对储层温度分布的影响规律，另一方面对2 种情况下的储层温度分布结果进行比较和分析。电磁波频率分别选取1.4，1.6 和1.8 GHz，对应的储层电导率分别为0.080，0.082 和0.085 S/m，相对介电常数分别为12.5，12.9 和13.2。

在情况2）的条件下，储层沿深度的温度分布受电磁波频率的影响规律如图9 所示。

从图9 可以看出，电磁波频率增大，能够使入射波储层面区域的温度急剧升高，但随着加热深度增加，温度急剧下降。因此，增大电磁波频率，有助于快速升温，但不利于增大加热距离。

通过求解数学模型，分别得到了2 种情况下储层温度沿深度的变化情况，见图10（图10 中，绿色曲线表示电磁波频率从1.4 GHz 变为1.6 GHz 或者1.8 GHz、而电导率和相对介电常数恒为1.4 GHz 时的恒定值（即0.080 S/m 和12.5）的计算结果，蓝色曲线表示根据式（1）和式（2）计算该频率下对应电导率和相对介电常数时的结果）。由图10 可知：在不考虑储层电导率和相对介电常数随电磁波频率变化的情况下，储层温度计算结果偏高，即储层性质不变时的温度计算结果要大于现场实际加热过程中的储层温度。

图 11 储层导热系数和比热容随温度变化和恒定时的温度分布对比Fig.11 A comparison of temperature distributions based on the constant and variable thermal conductivity and specific heat of reservoir with temperature

图 10 储层电导率和相对介电常数随电磁波频率变化和恒定时的温度分布对比Fig.10 A comparison of temperature distributions based on the constant and variable electrical conductivity and relative permittivity of reservoir with temperature

4.2 储层导热系数和比热容

在加热时间恒为6 h 的条件下，分2 种情况进行研究：1）基于实际情况，考虑温度变化引起储层导热系数和比热容改变，其变化规律依据式（3）和式（4）计算分析；2）储层导热系数和比热容不受温度影响，其值恒为0.2 W/（m·℃）和800 J/（kg·℃）。

通过求解数学模型，分别得到了2 种情况下的储层温度分布情况，如图11 所示。

从图11 可以看出：情况1）下的储层温度分布曲线与情况2）下的储层温度分布曲线存在交点；在该交点横坐标左侧的区域内，情况2）下的储层温度较高；在该交点横坐标右侧的一段区域内，情况1）下的储层温度较高。这说明在不考虑储层导热系数和比热容随温度变化的情况下，接近波源储层区域温度的计算结果高于储层实际温度，远离波源一定距离储层区域内温度的计算结果低于储层实际温度；随着储层深度进一步增大，2 条曲线同时下降，并趋于重合。以上温度分布曲线与文献[2]中温度分布曲线基本吻合，验证了以上结果的正确性，同时也验证了数学模型的有效性。

4.3 储层相对介电常数和电导率

研究不同稠油储层的相对介电常数和电导率对温度分布的影响规律时，电磁波频率恒为1.6 GHz，储层相对介电常数及电导率的取值范围可参考文献[29]。假设储层的相对介电常数分别为8.5，12.5和16.5，计算不同相对介电常数下的温度分布曲线，结果见图12。由图12 可知：不同相对介电常数下的储层温度随深度增大而降低，3 条温度分布曲线随着深度增大逐渐趋于重合；储层相对介电常数从8.5 增至16.5 时，储层深度小于0.15 m 时其温度均升高，表明储层相对介电常数越大，对电磁波加热开采稠油越有利。

假设储层电导率分别为0.01，0.08 和0.15 S/m，计算不同电导率下的温度分布曲线，结果如图13 所示。从图13 可以看出：随着电导率增大，沿深度方向的储层温度也明显升高，这与图12 中的温度分布曲线变化规律相似。以上研究表明，稠油储层相对介电常数和电导率在一定范围内增大，可以提高储层温度。