基于Matlab电力变压器励磁涌流的分析和仿真

基于Matlab电力变压器励磁涌流的分析和仿真
宁德师范学院
论文(设
专业指导教师
学生
学号题
目) 毕业计
目录
1 变压器空载合闸励磁涌流产生机理 (1)
1.1 变压器励磁涌流的定义 (1)
1.2 变压器励磁涌流产生的原因 (1)
2 变压器空载合闸物理过程分析 (1)
2.1 单相变压器的涌流分析 (1)
2.2 三相变压器的涌流分析 (4)
2.3 励磁涌流的影响及抑制措
施 (5)
3 变压器励磁涌流的仿真 (5)
3.1 变压器仿真模型构建 (5)
3.3 励磁涌流仿真结果的分析 (6)
3.4 励磁涌流与短路电流比较 (8)
4 结束语........................................................................................................................ (8)
参考文献........................................................................................................................ (9)
基于Matlab的三相变压器励磁涌流仿真分析
摘要：阐述了变压器空载合闸时励磁涌流产生的机理，在单相变压器空载合闸的理论基础上，运用Matlab电气系统模块库构建仿真模型，对三相双绕组变压器空载合闸的过程进行仿真及分析。

对不同状态下的励磁涌流做进一步分析，分析结果和理论分析相吻合，验证了仿真的有效性。

关键词：变压器；Matlab；励磁涌流
1 变压器空载合闸励磁涌流产生机理
1.1 变压器励磁涌流的定义
通常在正常运行的变压器中的励磁电流非常小，大约仅有额定电流的3%~8%，而大型电力变压器的励磁
[1]涌流还不足额定电流的1%，如此小的励磁涌流并不足以破坏电力系统的稳定性。

因为变压器本身的铁芯材
料呈非线性特性，并附带磁通饱和特性，导致在空载合闸的瞬间，会产生很大的冲击电流，该值可达额定电
[2]流的3~4倍，是正常空载运行电流的几十倍甚至百倍以上。

1.2 变压器励磁涌流产生的原因
对变压器的进行空载合闸操作有两种，即：（1）电力变压器的空载投入电网运行；（2）电网发生故障要
[3]切除变压器，待故障排除后变压器的再次投入。

如图1所示，是变压器铁芯近似磁化特性曲线。

从图中可
以看出，饱和曲线的延长线与坐标纵轴相交于点S，把S点的饱和磁通量定义为?s。

在正常运行状态下，饱和磁通介于0~?s之间变化，励磁阻抗很大，一般以变压器额定电压和电流为基准的励磁阻抗Zm?100，故变压器的励磁涌流i?很小，可近似为零；但是，当变压器空载投入时，变压器铁芯磁通量大于?s时，达到瞬变磁通?x，由下图可以看出，变压器励磁涌流i?沿着磁化特性曲线将迅速增大。

它的大小与变压器等值阻抗、合闸初相角、变压器铁芯剩磁大小、变压器绕组接线方式、变压器铁芯的材质及结构等诸多因素有关。

图1 变压器铁芯磁化曲线
2 变压器空载合闸物理过程分析
2.1 单相变压器的涌流分析
电力系统中的变压器中主要是三相变压器，但分析三相变压器的励磁涌流可以在分析单相变压器励磁涌流的基础上进行。

图2是变压器接线图，二次侧开路、一次侧在t?0时刻合闸到电压为u1
的电网上，其中：u11sin(?t??)
1
式中?为变压器的合闸初始相位角。

图2 空载合闸到电网接线图
在t?0期间，变压器一次侧绕组中电流i1满足如下微分方程式：
d??i1r11sin??t??? （1）dt
其中，?是与一次侧绕组相交链的总磁通，它包括主磁通和漏磁通。

由于现代变压器电阻非常小，在上
N1式中电阻压降i1r1较小，所以在分析瞬变过程中的初始阶段暂不考虑，这样可以更清楚的看出在初始阶段电流较大的物理本质。

r1存在是使瞬态分量衰减的基本原因，因此，在研究瞬态电流衰减时，必须计算及r1的影响。

当忽略r1时，式（1）变为：
N1d??1sin??t??? （2）dt
os?(t???)C （3）1解微分方程得：??其中，C由初始条件决定。

考虑到变压器空载合闸前磁链为0，根据磁链守恒原理，有：
?|t?0???|t?0??0
得：
C?
于是式（3）变为：
????cos(?t??)]??m[cos??cos(?t??)] （4）
1
?m?1 ?N11cos? ?N1
式中?m为稳态磁通最大值。

从式（4）可以看出，磁通?的瞬变过程与合闸时刻（t?0）电压的初始相角?有关。

下面讨论两种特殊情况。

?时合闸（即u11时合闸）。

由式（4）可得：2
? ????mcos(?t?)??msin?t
（5）2（1）在电压初相角??
2
这和稳态的运行情况相同，从t?0开始，变压器一次侧电流i1在铁芯中就建立了稳态磁通?msin?t，而不发生瞬变过程，一次侧电流i1也是正常运行时的稳态空载电流i0。

（2）在电压初相角??0合闸（即u1?0时合闸）。

由式（4）可得：
????mcos?t??m??'??" （6）'"其中，???m为磁通的瞬时分量，是一个常数，因忽略了电阻r1，故无衰减；????mcos?t为磁通
的稳态分量。

与式（5）、（6）对应的磁通变化曲线如图3、图4所示。

从t?0开始经过半个周期即t?
到最大值：?时，磁通?达?
?max?2?m
即瞬变过程中磁通可达到稳态分量最大值的2倍，此时的瞬变过程最为强烈，是最不利的合闸情况。

??
图3 ?? ?空载合闸时磁通曲线图4 ??0空载合闸时磁通曲线 2
还有一种情况是：在研究瞬态电流衰减时，由于有电阻r1存在，合闸电流将逐渐衰减。

衰减快慢由时间
L1决定，L1是一次侧绕组的全电感。

一般容量较小的变压器衰减的速度快，约几个周波就达到稳r1
[4]定状态；大型变压器则衰减的比较慢，甚至要持续几十秒才达到稳定。

有必要对励磁涌流的间断角进行分析，图5是图解法描绘的变压器励磁涌流波形。

当变压器正常运行时常数T?对应图b的0~?1和?2~2?这两段，当变压器铁芯饱和时，励磁电流急剧增加，如图?1~?2段所示。

此时可以看出，励磁电流在一个周期出现了间断，设间断角为?j，则有：?j??1?(2???2) （7）令式（6）中???t??，且当???1，???2时，均有???s，代入式（6）中，得：
?1?Arccos(
且有?1??2?2?
将式（8）代入式（7）得:
3
?mcos???s) （8）?m
?j?2?1?2Arccos(?mcos???s) （9）?m
实际上，间断角大小与变压器饱和程度有关，铁芯越饱和，励磁涌流的尖顶特性越明显。

而变压器铁芯的饱和程度又与变压器的剩磁，合闸相角等有关。

发生正向饱和时，正向剩磁越大，当合闸相角接近0°时饱和最为严重，间断角越小。

同理可得，反向饱和时，反向剩磁越大，饱和最严重发生在合闸角为180°，
[5]此时间断角最小。

ΦΦ1(a) 变压器铁芯磁化曲线(b)励磁涌流波形
图5 变压器励磁涌流图
综上分析可得，单相变压器励磁涌流的特征有这几点：
（1）变压器空载合闸励磁涌流的波形之间有间断角。

（2）包含大量的高次谐波，其中主要以二次谐波为主。

[6]（3）包含大量的非周期分量，使涌流偏于时间轴的一侧。

必须注意的是，以上的分析将变压器铁芯励磁特性作了线性化的简便处理，从而忽略了非线性、磁滞特性、局部磁滞回线等，从而使复杂的非线性问题简单化。

事实上，变压器铁芯饱和后的励磁电感并非常数，它随着饱和程度的加深而逐渐变小，由此可知实际的励磁涌流波形呈现的是尖
顶波。

以下对三相变压器励磁涌流的分析也作了类似简化。

2.2 三相变压器的涌流分析
三相变压器的励磁涌流同样的与电力系统电压大小与合闸初相角、剩磁大小与方向、铁芯材料等因素关系密切。

但因为其磁路的特殊结构，不同的接线组别等各方面特定因素，导致在分析时相比单相变压器来说要复杂不少。

但它有个明显的特点，即无论空载合闸的初相角或大或小，必定都会有励磁涌流产生。

理由是三相变压器的磁通都分别滞后三相外加电压90°，即三相的合闸角度互相差120°，故不能使三相变压器铁
[7]芯中的暂态磁通全为零。

三相变压器的绕组连接方式和磁路结构有许多种，对励磁涌流的大小和波形有一定影响。

大型变压器一般是由三个单相变压器组成的变压器组，由于三个铁芯的磁路各自完全独立，所以对单相变压器的分析方法同样适合于三相变压器。

最为常见的变压器接线方式为Y/?接法，当Y侧空载合闸时，变压器
一次侧产生的励磁涌流ia，ib，ic。

变压器差动保护的二次电流相位调整通常采用星形向三角形变换来调整差流平衡，这样一来。

电流互感器二次侧流入差动保护的电流实际上是变压器一次侧的两相电流之差，如下图所示：ijA?ia?ib ijB?ib?ic ijC?ic?ia
4
h1
图6 YN，d接线三相变压器接线图
因此，要分析的励磁涌流实际上是变压器一次侧两相涌流的差值。

因
为这个特殊的关系，三相变压器励磁涌流在合闸初相角为±30°时，幅值才会出现最大值。

一次侧三相电流中两相涌流的方向相同且直流分量相差不大时，二次侧涌流中就会有一相涌流的直流分量很小，甚至为零，波形特征体现在该相涌流对称于时间轴，称为对称涌流。

对称涌流是由剩磁方向相同的两相涌流相减产生的的电流。

与之相对应的直流分量较大的涌流称为非对称涌流。

而它的来源是剩磁方向相反的两相涌流相减生成的。

大量的分析研究表明，三相变压器的空载合闸励磁涌流特点归纳如下：（1）波形肯定存在间断角。

不论是单侧性或是周期性涌流的波形都
含有不同程度的间断角，周期性涌流的间断角较小。

（2）三相变压器励磁涌流同样包含较大的二次谐波分量，其含量与
饱和磁通和合闸角直接关联。

（3）由于三相电压之间存在120°的相位差，因此导致三相励磁电流不相同。

任何情况下的空载合闸，
[8]至少有两相会出现不同程度的励磁涌流。

2.3 励磁涌流的影响及抑制措施
变压器在空载合闸或切除外部故障切除后重新投入供电时，虽说对变压器直接危害不大，但是流经变压器的电流、电压波形产生严重畸变，并让安装在变压器一次侧的过电流保护继电器动作，导致变压器合不上闸，如遇到这种情况可以再合闸一次，可以进行多次尝试，总能在恰当的时刻合闸成功。

除此之外，过大的励磁涌流还带有大量的二次谐波，严重破坏了电网的电能质量，由于电力电子器件敏感性较强，谐波分量很可能严重破坏电
力电子器件，造成电力电子器件失控或损坏，而且还会产生谐振过电压，使连在变压器附近的电气设备无法正常工作。

传统的抑制方法是保护装置采用谐波制动、放大动作电流定值等方法躲过励磁涌流，可是还没解决根本本矛盾。

目前主要有两种抑制励磁涌流的方法，一是通过变压器外部控制削减励磁涌流。

二是通过内部控制，即通过改变变压器的内部结构达到削弱励磁涌流的目的。

这些方法的实现方式分别如下：
（1）并联合闸电阻，通过合闸电阻承受冲击电流，在冲击电流衰减到一定范围后再切除合闸电阻。

（2）内插电阻法，在三相变压器的中性点处联接一个接地电阻，以承受这种不平衡电流，从而使得变压器的励磁涌流得以衰减。

（3）变压器选相分合闸技术，利用电压与磁通的相位关系，通过控制开关合闸时间来达到抑制励磁涌流的目的。

（4）改变变压器绕组的分布，通过改变变压器的结构增加暂态等效电感来抑制励磁涌流。

由于需要改变变压器的结构，它的发展受到一定限制。

[9]（5）变压器低压侧并联电容器，这种方法通过抑制变压器磁通达到饱和从而抑制了励磁涌流。

3 变压器励磁涌流的仿真
3.1 变压器仿真模型构建
如图7所示，是双侧电源双绕组的变压器的电力系统。

5
图7 双侧电源的双绕组变压器电力系统
图8 Simulink仿真模型
根据需要的从Matlab电力系统工具箱中提取相关元件建立仿真模型，对变压器进行仿真，图8就是对应的Simulink仿真模型。

对仿真模型进行简化分析，故设置变压器两侧接线方式和电压等级均一致，并令模型中的变压器铁芯饱和。

其他参数的设置如表1所示。

表1 仿真模型相关参数项目参数设定
35kV 电压源额定电压
变压器类型三相双绕组
50MVA 变压器额定容量
35kV 变压器额定电压
0.002 一次侧和二次侧绕组电阻标么值
0.08 一次侧和二次侧绕组漏电感标么值
500 励磁电阻标么值
仿真模型中的UM和UN模块所起到的作用相当于电压互感器和电流互感器，用来采集电压、电流信号并转化成分析所需要的simulink信号。

断路器QF1和QF2用于控制变压器的空载合闸，其中QF1动作时间设定为0，要消除故障点和QF2在模型中的作用，只要把动作时间设定的比仿真运行时间长即可。

[10]
为了观察励磁涌流的波形，可以通过示波器模块，如图9所示。

图9 示波器模块
3.3 励磁涌流仿真结果的分析
根据需要建立仿真模型，把变压器设定为Y/Y连接，合闸初相角设定为0°，仿真时间设为0.45s，运行仿真，得到的涌流波形如图10所示。

6
下面是把三相变压器接线改为Yn/d11连接，合闸初相角和其他参数均不变，重新运行仿真，得到图11的励磁涌流波形。

图11 Yn/d11接法的三相变压器励磁涌流
由波形图可以看出励磁涌流都不可避免的产生了间断角，间断角和励磁涌流的二次谐波有关系，间断角越大，二次谐波涌流就越大。

包含有很大成分的非周期分量，导致涌流偏于时间轴的一侧，也可以很明显的发现三相变压器的接法对励磁涌流影响甚大，就Ia来说，涌流的最大值都相差了将近2倍，，并且不论是哪一相，Y/Y接法涌流的衰减速度明显快于Yn/d11接法。

下面采用FFT Analysis（快速傅立叶分析）对变压器的涌流进行分析。

接线仍为Yn/d11连接，更改变压
7
器A相空载合闸的初相角，得到如表2所示的分析表。

表2 不同合闸初相角下的励磁涌流谐波分析表合闸初相角励磁涌
流(%) 直流(DC) 基波(Fund) 三次谐波(h3) 四次谐波(h4) 五次谐波(h5) A相100 6.04 5.41 3.43
0°B相100 7.16 8.82
C相100 12.01 8.52
A相100 7.51 8.58
60°B相100 11.06 8.03
C相100 5.40 5.64 3.35
A相100 11.03 23.11 7.84
120°B相100 4.89 5.86 3.30
C相100 7.26 23.08 8.27
A相100 4.77 5.89 3.32
180°B相100 7.26 7.97
C相100 12.02 7.79
63.29 57.81 54.84 57.49 55.41 63.49 55.19 63.64 57.79 63.68 58.44 54.43
二次谐波(h2) 45.48 50.65 55.49 51.00 53.96 44.88 54.25 44.42 50.24 44.29 48.84 55.72
23.35 23.51 23.36 23.26
22.77 23.13
从表2可以看出，不论哪一相的励磁涌流谐波含量都很大，并且以二次谐波为主。

涌流包含很大的直流分量，但这不能证明在其他初始条件下也是如此。

3.4 励磁涌流与短路电流比较
设置故障模块1，使变压器空载合闸后在0.2s~0.3s发生三相短路故障，三相变压器为Y/Y型连接，合闸初相角角为0°，其他参数不变。

运行仿真得到的波形如图12所示。

图12 励磁涌流和短路电流的区别
通过对比可以很容易的区别二者，A相励磁涌流峰值比短路电流小，B、
C两相的励磁涌流峰值比短路电流大，短路电流的衰减速度也比励磁涌流的衰减速度快。

4 结束语
在对空载合闸励磁涌流的理论分析的基础上，合理利用Matlab软件进行仿真使得分析过程变得比较简单方便，同时，得益于Matlab软件功能比较丰富，还可以对励磁涌流进行FFT分析，并且可以将励磁涌流和短路电流相比较，可以很直观的在波形上就看出二者的区别，从而避免了某些情况下对电网的错误操作。

尤
8
其是在变压器保护这一方面提供了很大的帮助，是一种不错的分析手段。

参考文献：
[1] 王展，常勇，刘春雪．变压器励磁涌流的分析及抑制方法探讨[J]．中国电力教育，2009(2)：169-171．
[2] 谢明琛，张广溢．电机学[M]．重庆：重庆大学出版社，2010，53-55．
[3] 王巨丰，何振东，李明武，杨文斌．利用分相关断策略消除电力变压器的励磁涌流[J]．广西电力，2005(5)：43-45．
[4] 谢应璞．电机学上册[M]．成都：四川大学出版社，2007，87-89．
[5] 白全林. 电力变压器励磁涌流的分析与研究[D]．陕西：西安科技大学，2012：10．
[6] 何越，熊元新，姜山，徐迎，魏徐．变压器空载合闸励磁涌流的仿真分析研究[J]．电力学报，2010，25（1)：33-42．
[7] 孙鸣，江海峰，丁明．变压器差动保护的新原理研究[J]．煤炭学
报，2003，28（3）：311-315．[8] 方愉冬. 基于小波变换的变压器保护原理探讨及实现[D]．浙江：浙江大学大学，2002：7．
[9] 陈仁云，梁娇兰．基于Matlab的抑制变压器励磁涌流仿真研究[J]．南方电网技术，2012，6（3)：86-89．
[10] 于群，曹娜．MATLAB/Simulink电力系统建模与仿真[M]．北京：机械工业出版社，2011，163-166．
9。