刘士利 罗英楠 刘宗烨 邓金锁 蔡国伟

（1. 现代电力系统仿真控制与绿色电能新技术教育部重点实验室（东北电力大学） 吉林 132012 2. 先进输电技术国家重点实验室（全球能源互联网研究院有限公司） 北京 102209）

0 引言

随着能源枯竭和环境污染问题越发严峻，目前全球海上风电开发呈现出深远海化和集群化的发展趋势[1]，离海岸距离超过100km 的深远海域具有更丰富的风能资源，是未来海上风电探寻的重点方向[2]。德国等部分欧洲国家的深远海风电场已经投运[3]，中国已经进行了基础前瞻性探索[4-6]，其中，大规模风电跨海远距离输送是深远海域风能开发涉及的关键问题之一。

海上风电输送技术难度较大和成本较高，一般采用电缆线路。采用直流输电方式虽然可以减小损耗提升载流量，但直流电压下空间电荷积累效应给电缆绝缘带来了严峻挑战，在一定程度上限制了输送容量的提升。若采用工频交流输电，则由电容效应导致的无功电流和电磁感应引起的护套感应电流将显著影响电缆的输送容量、传输距离和经济性。而介于工频和直流之间的低频乃至分频输电方式可能在海上风电并网中得到广泛应用[7-9]。因此电缆在低频输电方式下的损耗和载流特性是深远海风电高效输送的关键因素之一。

目前，关于电缆护套感应电压、环流损耗的研究多集中于工频输电方式，采用的方法有IEC 标准法及有限元方法。文献[10]分析了多回路电力电缆不同敷设方式对环流的影响，指出了不同相序组合下的环流分布规律；文献[11]基于光电复合海缆的有限元模型，提出了光纤温度变化模型；文献[12]建立了光电复合海缆的热路模型；文献[13-15]采用有限元多物理场仿真方法，分析了不同敷设条件下海底电缆的温度和载流特性；文献[16]通过建立环形纹和螺纹护套的参数方程，提出了线芯与波纹护套的互感和等效直径方法误差的解析公式；文献[17]采用有限元法研究了护套材料磁导率、体积电阻率对护套感应电压的影响；文献[18]基于电路模型推导了双回路敷设电缆护套环流的计算方程，提出了减小护套环流的技术方案；部分学者还研究了电缆及周围媒质的温度场分布和热特性[19-21]。而关于电缆在低频输电方式下的损耗和载流特性研究较少。文献[22]基于系统技术经济性，讨论了分频输电方式下交流电缆载流量、传输极限等问题；文献[23]研究了城市电缆线路在低频输电方式下的运行特性。综上所述，关于三芯海底电缆在低频输电方式下的损耗特性，鲜有文献报道。

本文采用有限元方法，建立了三芯铠装电缆的电磁-热多物理场仿真模型，研究了其在低频运行方式下的损耗和载流特性，分析了缆芯电流在低频输电方式下的不均匀分布特性，研究了输电频率、缆芯半径等因素对护套感应电压、损耗因子及载流量的影响，并分析了IEC 标准计算误差，计算结果为提高深远海电缆线路的经济性提供了数据支撑。

1 三芯铠装电缆护套感应电压解析计算

运行过程中电缆护套损耗由护套上的感应电压和电流引起，因此在分析电缆护套损耗之前应首先明确护套感应电压的计算。相较于单芯电缆，三芯铠装电缆能大幅度节省输电走廊，提高传输容量，因此海上风电传输往往采用三芯光纤复合海底电缆，其结构示意图如图1 所示。

图1 三芯铠装电缆结构示意图Fig.1 Structure of three-core armored cable

根据电磁场理论，缆芯电流激发的交变磁场将在电缆护套和铠装中产生感应电压，设三芯电缆流过三相平衡电流，其复数形式为

式中，I为三相电流有效值。

按IEC 标准，当导体P 和A、B、C 三相缆芯平行布置时，如图2 所示，则三相电流产生的与P导体交链的磁通可表示为式中，RP为导体P 的几何半径；S1、S2、S3分别为三相缆芯与导体P 之间的距离。

图2 电缆护套感应电压等效分析模型Fig.2 Equivalent analysis model for induced voltage of cable sheath

对于三芯铠装海底电缆，护套可视为上述导体P 的特例，其与缆芯的间距等于护套平均半径，此时，护套感应电压计算式为

式中，d为护套直径的几何平均值；S为缆芯轴间距离。

上述护套感应电压计算公式是基于电缆简化模型推导得出，模型首先假设电流均匀分布，在处理某些环节时进一步假设电流集中于缆芯中央。而实际电流则由于趋肤效应和邻近效应的影响，在缆芯中并非均匀分布；且随着传输容量的上升，缆芯半径越大、结构越复杂，电流不均匀分布越明显，解析法计算就愈发困难。因此，对于三芯铠装海底电缆，采用上述解析公式进行计算，精度难以满足要求，有必要建立精确的有限元场模型进行分析。

2 三芯电缆电磁特性有限元分析

2.1 电磁-热耦合分析模型

电磁-热多物理场耦合模型以电磁场的麦克斯韦方程组和热传导的傅里叶定律为基础，电缆中电磁参数满足

式中， ˙H为磁场强度矢量相量（A/m）； ˙E为电场强度矢量相量（V/m）；σ为电导率（S/m）；ω为角频率（rad/s）；ε为介电常数（F/m）；μ为磁导率（H/m）。

引入相量形式的矢量磁位˙A，上述方程可以转化为频域下的二阶偏微分方程，即

由上述方程得出矢量磁位˙A之后，可按式（6）求得磁感应强度。

式中，Φ为与金属护套交链的磁通。

电缆实际运行过程中，缆芯和金属护套的电导率均受温度影响，其与温度的函数关系为

基于电磁场和温度场双向耦合的电缆损耗分析流程如图3 所示。

图3 电磁-热耦合分析原理Fig.3 Principle of electromagnetic-thermal coupling analysis

2.2 缆芯电流分布及护套感应电压

三芯电缆的几何尺寸见表1，敷设条件如图4所示。通过电磁-热多物理场仿真计算，三芯铠装电缆的磁场分布如图5 所示。

图4 电缆敷设条件Fig.4 Cable laying conditions

图5 电缆磁场分布Fig.5 Magnetic field distribution of cable

表1 电缆物理参数Tab.1 Cable physical parameters

当输电频率为50Hz 时，缆芯和护套中的电流分布如图6 所示。可以看出，缆芯电流明显呈不均匀分布，在三根缆芯相互靠近的一侧电流密度较大，在远离的一侧较小。缆芯电流的不均匀分布由趋肤效应和邻近效应引起。缆芯电流分布随频率的变化趋势如图7 所示。可以看出，频率越高，电流不均匀度越明显；且电流密度平均值位置和缆芯中心偏离较大，该偏离距离正是导致IEC 标准计算护套感应电压存在误差的原因之一。可以预见，电缆结构越复杂，其误差越大。

图6 缆芯与护套的电流密度Fig.6 Current density of cable core and sheath

图7 缆芯电流径向分布Fig.7 Radial distribution of core current

不同频率时，采用有限元法计算的护套感应电压如图8 所示，可以看出，采用电流均匀分布假设时，计算结果偏大，且随着频率增加，二者之间的偏差逐渐增大，分频时误差为8.4%，工频时的误差为13%。

图8 电缆护套感应电压Fig.8 Induced voltage of cable sheath

2.3 护套损耗

电缆护套中的感应电流会产生额外的损耗，从而导致电缆温度进一步上升，因此，准确计算护套损耗对于合理确定三芯海底电缆载流量，提升海上风电的经济性具有重要意义。根据IEC 标准，在不考虑铠装的情况下，护套中的总损耗由环流损耗和涡流损耗构成，损耗因子可表示为

式中，X为电缆护套电抗（Ω/m）；Rs为护套电阻（Ω/m）；R为缆芯电阻（Ω/m）；S为三相导体间的缆芯轴间距离(m)。

式中，ts为护套厚度（mm）；Ds为电缆护套的外径（mm）；ρ0为护套材料的电阻率（Ω·m）。

上述公式在计算λ1时没有考虑邻近效应，这同样导致公式在计算三芯海底电缆损耗因子时精度降低。输电频率一方面影响缆芯电流分布、缆芯电阻及缆芯欧姆损耗，另一方面还影响护套感应电压及护套损耗，而IEC 标准在计算λ1过程中做出的某些假设并未计及上述影响。此外，若考虑铠装，1λ′可由式（14）给出，此时电缆结构更加复杂，IEC 标准的计算误差会更大。因此，需要构建有限元模型进行详细讨论。

采用有限元分析可以得出，输电频率及缆芯半径对缆芯电阻的影响如图9 所示，可以看出，由于趋肤效应和邻近效应的存在，随着频率的增加，电阻明显增加，从而将导致缆芯损耗增加。

图9 输电频率及缆芯半径对电阻的影响Fig.9 Influence of transmission frequency and core radius on resistance

输电频率对损耗因子的影响如图10 所示，可以看出，采用有限元计算时，如采用电流均匀分布模型，其计算结果和IEC 标准计算结果基本相同，且大于实际值（有限元采用电流不均匀分布模型计算结果），这一方面验证了IEC 标准的假设条件，另一方面也揭示了标准的局限性。随着频率的上升，趋肤效应和邻近效应越发明显，电流分布不均匀导致的损耗比率λ1也明显变化，当频率由分频50/3Hz 增加到工频50Hz 时，λ1误差由14%增加到18.4%。

图10 频率对λ1 的影响Fig.10 Influence of frequency on λ1

除了频率，缆芯半径同样影响损耗因子λ1，不同频率时缆芯半径对1λ的影响如图11 和图12 所示。在工频输电方式下，IEC 计算结果显示，λ1随缆芯半径增加而逐渐增加，并存在饱和趋势；而有限元结果则显示，当缆芯半径小于24mm 时，λ1随缆芯增加呈增长趋势，当缆芯半径大于24mm 时，λ1则随缆芯增加呈现下降趋势，此时IEC 标准法误差将进一步增大，究其原因在于，缆芯半径越大，电流密度平均值位置和缆芯中心偏离越大，导致IEC 标准的电流集中分布假设和电缆实际电流分布情况相差较远，从而导致IEC 计算结果误差进一步增大。分频输电方式下，缆芯半径对λ1的影响与工频情况下类似。总之，电流的不均匀分布对损耗有十分明显的影响，因此在计算电缆损耗时，应将电流分布情况考虑在内。

图11 50Hz 时缆芯半径对λ1 的影响Fig.11 Influence of core radius on λ1 at 50Hz

图12 17Hz 时缆芯半径对λ1 的影响Fig.12 Influence of core radius on λ1 at 17Hz

当同时计及护套损耗和铠装损耗时，采用有限元计算的电缆各部分损耗如图13 所示，可以看出，随着频率的上升，电缆整体电磁损耗、缆芯损耗、铠装及护套损耗均呈上升趋势。电缆总损耗比率如图14 所示。

图13 电缆各部分损耗Fig.13 Loss of each part of the cable

图14 电缆损耗比率Fig.14 Cable loss ratio

2.4 载流量

基于上述分析，充分考虑缆芯电流不均匀分布，采用电磁-热耦合模型计算得到的三芯铠装海底电缆的载流量如图15 所示。

图15 电缆载流量Fig.15 Ampacity of the cable

可以看出，相同缆芯截面积时，电缆工频载流量低于分频载流量；和有限元结果相比，IEC 计算结果偏大，随着缆芯截面积增大，IEC 计算误差有饱和趋势，工频时误差最大值约为15%，分频时误差最大值约为10%。图16 为三芯电缆满载时的温度场分布，此时缆芯温度为90℃，温度由内向外逐渐降低。

图16 电缆温度场分布Fig.16 Distribution of cable temperature

3 结论

本文基于电磁-热耦合原理研究了三芯铠装海底电缆在低频输电方式下的损耗特性，得出如下结论：

1）缆芯电流密度平均值点偏离缆芯中心，偏离距离随输电频率的增大而增大，此偏离距离是采用IEC 标准计算感应电压、损耗因子和载流量时存在误差的重要原因之一。

2）采用IEC 标准计算的三芯铠装电缆护套损耗因子λ1偏大，输电频率越高，缆芯半径越大，λ1误差越大，分频及工频输电时误差分别达到14%和18.4%；另外，当缆芯半径增大时，实际损耗因子λ1变化规律呈现饱和甚至减小的趋势，此时（半径为24mm）由于IEC 标准计算误差进一步增大，应慎重采用。

3）采用IEC 标准计算的电缆载流量偏小，该误差随输电频率和缆芯截面积的增加而增大，当截面积取2 000mm2时，分频及工频下的载流量误差分别达10%和15%左右。