变压器地磁感应电流-无功功率动态关系分析

2021-05-12王泽忠黄天超

电工技术学报 2021年9期

王泽忠黄天超

（华北电力大学高电压与电磁兼容北京市重点实验室北京 102206）

0 引言

太阳活动引起的地磁场剧烈变化称为磁暴。时变磁场在地面感应出的电场通过输电线、中性点接地变压器和大地构成的回路产生地磁感应电流（Geomagnetically Induction Current, GIC），其频率一般为 0.000 1～0.01Hz[1]。由其引起的变压器铁心半波饱和，会导致谐波、局部温升、无功损耗和振动、噪声增大等有害的干扰[2-8]。其中谐波和无功功率的波动会影响电力设备的安全稳定运行及整个系统的电压调控，因此危害最大[9]。近些年随着输电线路电压等级越来越高，导线电阻越来越小，变压器遭受的GIC 越来越大[10]。因此精确预测在地磁暴期间的变压器无功功率损耗就显得十分重要。

GIC 对不同类型变压器无功损耗影响大小以及两者之间的直接关系（以下简称为GIC-Q），国内外已有大量的研究[11-17]。目前普遍认为，单相变压器由于磁路结构的特点受GIC 影响导致的无功功率波动最大[18]；虽然谐波电流受变压器励磁曲线影响，但在实际测得的GIC 数值范围内，基波无功功率与GIC 之间呈线性关系[19]。然而上述研究直接把GIC等效为直流电流，没有考虑GIC 的准直流特性，尤其是在GIC 由一个正峰值向一个负峰值变化的阶段，其变化率较大，此时GIC 的这种准直流特性不可忽略。文献[20]首次提出不考虑GIC 的准直流特性，可能会导致对GIC 风险评估的降低。因此IEEE在Std C57.163—2015 标准中，将一个地磁暴事件分割为多个小区间，每个区间对应一个直流电流，以此来考虑GIC 的准直流特性对温升计算的影响，但这种方法并没有涉及对铁心动态过程的分析。文献[21]认为对于含有三角形联结绕组的变压器，不考虑GIC 的准直流特性会影响其无功功率的计算精度，并通过理论计算与实测数据的对比，证明了结论的正确性。但其没有考虑实测GIC 数据的精度对计算结果的影响，也没有针对具体变压器进行分析。

因此本文通过搭建磁暴期间变压器铁心准直流零序等效电路模型，分析研究了GIC 的准直流特性对于无功损耗的影响，相比于将GIC 当作纯直流来研究，无功损耗的计算精度明显提高。通过使用不同步长的GIC 数据，对一台1 000kV 交流变压器和一台800kV 换流变压器进行了仿真计算，证明了若直接使用原始GIC 秒数据进行无功损耗计算会造成较大误差，必须对原始数据进行插值处理，并进一步总结了GIC 准直流特性对无功损耗的影响规律。最后通过改变变压器参数，确定了各参数对无功损耗影响的灵敏度大小。研究结论可为准确预测地磁暴期间变压器无功损耗及评估其遭受GIC 侵害大小提供理论依据。

1 静态无功功率理论计算

由于变压器铁心使用的硅钢片在饱和前后微分磁导率相差的非常大，因此可以用图1 所示的两条线性直线表示其i-ϕ关系，即忽略曲线在拐点附近的过渡过程。当不考虑GIC 的动态特性时，其在某一小时间段内可以完全等效为直流电流，因此GIC的静态无功功率理论推导可以当作是直流偏磁情况来处理。

图1 直流偏磁时变压器铁心工作状态Fig.1 Transformer core operating condition under dc bias

图1 中，Lair-core、Lu分别为铁心饱和、非饱和微分电感，ibias为直流磁通对应的直流电流分量，iu为额定励磁电流峰值。当GICIGIC流入变压器时，铁心中便会产生相应的直流磁通量dcΦ，加上峰值为acΦ的额定交流磁通量，此时铁心中总的磁通量最大值mΦ将会超过膝点磁通kΦ使得铁心进入饱和状态。图中α定义为过磁通角，直流电流越大该角数值越大，直流为零时，α≈0 。图1 中各个磁通的关系如下。

式中，ω为工频角频率。

由图1 可知在GIC 的作用下，一个周期内励磁电流的表达式为

式中

对励磁电流进行傅里叶分解，可得其直流含量IGIC和基频分量有效值I1的表达式分别为

当变压器遭受直流侵害时铁心饱和，励磁电感迅速减小，此时单台变压器损耗的无功功率Q为正常时的无功损耗uQ加上励磁电感消耗的无功mQ，对应I1的两部分。假设变压器一次电压在GIC 入侵期间波形保持不变为额定基频电压，且与励磁电流的相位相差90°，即忽略铁心磁滞、涡流损耗。因此直流偏磁状态下，单台变压器的无功损耗为

式中，U1、Um分别为额定相电压的有效值和峰值。由式（7）结合式（5）和式（2）即可得到任意GICIGIC下，单台变压器消耗的无功损耗。

2 GIC-Q 动态关系分析

在磁暴期间，假设A、B、C 三相系统的参数都是一样的，那么流经变压器每一相中的GIC 无论大小还是方向都是完全相同的，又因为其具有准直流特性，因此可以将GIC 看成是一组零序电流。直流输电工程一般为12 脉动结构，因此每一个单元包含一台阀侧为三角形联结的换流变压器，而交流变压器为了补偿系统无功损耗，一般也含有三角形联结的第三绕组。由于三角形联结绕组的存在，具有准直流特性的GIC 将会在三角形联结绕组中感应出循环电流，从而抵消铁心中的准直流磁通。因此在磁暴期间，流入变压器中的GIC，不是完全作用到铁心中并产生相应的无功损耗。图2 为变压器准直流零序等效电路，表示了在磁暴期间变压器铁心的动态特性。

图2 变压器准直流零序电路Fig.2 Transformer quasi-DC zero sequence circuit diagram

图 2 中mi为实际作用到铁心产生直流磁通的电流，0R为三角形联结绕组的零序等效电阻，Ri为三角形联结绕组中的电流。因为GIC 的频率很低，因此变压器绕组漏电抗产生的电压降很小，可以忽略。变压器铁心的准直流特性用电感Lmdc表示。由图1 可知直流磁通Φdc的表达式为

式中，ω=2πf，f为系统工频。当变压器铁心工作在线性区域时，即铁心磁通最大值mΦ小于膝点磁通kΦ，直流电流与直流磁通的关系可由非饱和微分电感uL表示，所以此时准直流电感Lmdc=uL。当变压器铁心逐渐进入饱和状态时，即mΦ＞kΦ时，实际作用到铁心的直流电流mi由两部分组成，即

由定义可知，此时准直流电感Lmdc为

因此铁心饱和时的准直流电感Lmdc在数值上可以看成是电感uL和电感Ldc2的并联。

综上所述准直流电感Lmdc的表达式为

式中，ik为变压器铁心达到膝点磁通所需的直流电流，其与kΦ的关系为

此时准直流电感Lmdc随电流mi变化关系如图3 所示。图3 中数值由表1 中交流变压器参数计算得到。

图3 准直流电感与直流电流关系Fig.3 Diagram of relationship between quasi-DC inductance and DC current

表1 变压器参数Tab.1 Transformer parameters

根据图2 可以列出微分方程为

化简可得

运用向后差分法，即可求得磁暴期间变压器三角形联结绕组中的电流iR以及实际作用于铁心的直流电流mi，进而计算出对应的无功损耗。

3 仿真计算

为了对比分析，本节分别针对1 000kV 交流变压器和800kV 换流变压器进行了仿真计算。两种变压器的参数见表1，其中交流变压器为三绕组自耦变压器，第三绕组为三角形联结，而换流变压器为双绕组结构，阀侧为三角形联结绕组。GIC 采用2004年11 月9 号广东岭澳核电站的实测秒数据，该电流的峰值为75A，是我国目前为止测量到的最大数值。为了提高计算效率，仿真时间只截取了包含GIC 峰值在内的5min 这一时段，如图4 所示。

图4 实测GIC 数据Fig.4 Measured GIC data

由于ibias与IGIC相比很小，uQ与Q相比也很小。因此式（5）和式（7）可化简为

图5 为函数sinα−αcosα和α−0.5sin(2α)的关系图，由图可知当α＜0.846rad，即48.5°时，函数α−0.5sin(2α) ≈1.9(sinα−αcosα)，此时IGIC和Q的关系可以表示为

图5 函数关系图Fig.5 Functional diagram

由于使得α＞48.5°所对应的GIC 数值很大，实际情况中并不存在，因此可以认为在图4 所示的GIC 下，GIC 与无功损耗为线性关系。

图6 为不考虑GIC 准直流特性时，由式（2）、式（5）、式（7）计算得到的换流变压器无功损耗数值。通过对比图4 和图6，也验证了GIC 与无功损耗的线性关系。又因为在考虑GIC 准直流特性时，实际作用于铁心的直流电流mi为流入变压器中的实测GICiGIC减去三角形联结绕组中的电流Ri。因此三角形联结绕组中电流越大，变压器无功损耗受GIC 准直流特性影响变化量也越大。

图6 不考虑GIC 准直流特性时换流变压器无功损耗Fig.6 The reactive power loss of converter transformer with no quasi-DC characteristic

3.1 GIC 测量精度对仿真结果的影响

因为实测数据是每秒记录一次，为了避免GIC数据精度对计算结果造成影响，对GIC 秒数据进行线性插值，每两个数据点之间分别插入50 个和10个点，对应的仿真步长为0.02s 和0.1s。

图7 为分别使用步长为0.02s、1s、0.1s 的GIC数据计算的换流变压器无功损耗变化图（变化量为考虑GIC 准直流特性和将GIC 当作纯直流计算得到的无功损耗差值）。由图可知若直接使用GIC 秒数据计算，变压器的无功损耗与将GIC 当作纯直流计算相比变化较小，而使用经过插值处理的GIC 数据计算得到的无功变化却较大。从式（16）可以看出，准直流电感mL与iGIC变化率越大，三角形联结绕组中的电流Ri越大，无功损耗与不考虑GIC 准直流特性时相比变化就越大。而图3 则表明当作用于铁心的直流电流mi小于ki时电感mL最大，之后迅速减小，又因为变压器非饱和微分电感uL非常大，由式（14）计算出的ik只为0.365 A，因此若实测GIC数据采样间隔太大，如某一时刻GIC 为1A，1s 后测得GIC 为-2A，则计算时将不会考虑准直流电感为uL的情况，这会对无功计算造成较大的误差。

另外使用步长为0.02s 和0.1s 的GIC 数据计算得到的无功变化量基本一样，这说明使用步长为0.1s 的GIC 数据可以得到准确的计算结果，后续的计算都使用此数据。

图7 换流变压器无功损耗变化量Fig.7 Variation of reactive power loss in converter transformer

3.2 GIC 准直流特性对无功影响规律

通过对比图7 和图4 还可以看出，在GIC 流向发生转变所对应的波形过零点阶段，无功损耗发生了较大变化，而其余时段无功损耗与不考虑GIC 准直流特性时相比基本没有变化。这是因为在GIC 流向发生转变所对应的波形过零点阶段，准直流电感和GIC 变化率与其他时段相比都是最大的，所以在此阶段三角形联结绕组中电流激增，造成无功损耗较大的变化。另外当GIC 数值由正变负时，即流向变压器电流由流入变为流出，由于iGIC变化率为负，导致计算出的三角形联结绕组中的电流iR为负值，因此在这种情况下，实际作用于铁心的直流电流mi大于流过变压器的GIC，造成无功损耗的增大。当GIC 流向由流出变压器变为流入时，结论则相反。

3.3 两种类型变压器差异分析

由于换流变压器和交流变压器的额定电压等级不一致，为了对比分析两者无功受GIC 准直流特性影响的大小，将交流变压器额定电压归算到换流变压器的额定电压，计算的无功变化量如图8 所示。对比图7 和图8 可知，两种变压器无功损耗有着相同的变化规律，都是在GIC 流向发生转变所对应的过零点阶段发生较大变化。但是交流变压器的无功变化却小于换流变压器的。这是因为换流变压器正常工作时，由于换流阀触发延迟角不平衡等原因，阀侧绕组中存在较小的直流电流，因此在设计换流变压器时会将其工作磁通设计得较低，所以其额定磁通比交流变压器小。而两者铁心硅钢片使用的型号分别为27ZH095 和27ZH100，非常接近。再加上都是单相四柱结构，所以其膝点磁通、非饱和微分电感、饱和微分电感相差不大。表1 列出了具体数值。因此结合式（16）就可以得出如下结论：在其余参数相同的情况下，额定磁通较低的换流变压器相比于交流变压器，无功损耗更容易遭受GIC 准直流特性影响。

图8 交流变压器无功损耗变化量Fig.8 Variation of reactive power loss in AC transformer

4 参数灵敏度分析

本节针对仿真参数对无功损耗的影响进行了灵敏度分析。由表1 和图7、图8 可知，变压器额定磁通越大其无功损耗受GIC 准直流特性影响就越小，且在额定磁通增加了11.6%的情况下（以表1 中换流变压器的参数为基准），无功损耗变化量的最大值减小了39 %。通过第3 节的分析可知，膝点磁通对变压器无功损耗的影响机理与额定磁通一样，因此两者的灵敏度分析一致。

图9 仿真参数灵敏度分析图Fig.9 Sensitivity analysis diagram of simulation parameters

图9 为饱和微分电感Lair-core、非饱和微分电感Lu、三角形联结绕组电阻R0的灵敏度分析图，图中的变量均是以相应的换流变压器参数为比较对象。对比图7、图9 可知，当饱和微分电感变为原来的两倍时，无功损耗变化量数值基本不变。而非饱和微分电感越小，三角形联结绕组电阻值越大，无功损耗受GIC 准直流特性影响变化量越小。具体变化率为：非饱和微分电感减小为原来的一半，无功损耗变化量最大值减小了11.7%；三角形联结绕组电阻值增加为原来的两倍，无功损耗变化量最大值减小了32.75%。

综上所述变压器额定磁通和膝点磁通对其无功损耗受GIC 准直流特性影响最大，三角形联结绕组电阻和非饱和微分电感的影响依次减小，饱和微分电感对其几乎没有影响。