多绕组变压器复合短路阻抗与环流计算*

2024-01-15李世军夏湘滨杨佳意邝孝威陈懿媛郭盈

电气工程学报 2023年4期

李世军黄超夏湘滨杨佳意邝孝威张彬陈懿媛郭盈

(1.湖南工程学院电气与信息工程学院湘潭 411101；2.湖南华夏特变股份有限公司湘潭 411101；3.邵阳学院电气工程学院邵阳 422004；4.湖南德意电气有限公司湘潭 411101)

1 引言

多绕组变压器能够在电力工程领域应用广泛[1-5]，得益于它独特的优势和特点。因为它有多绕组且容量大，所以能输出多数值电压，能较大程度节约设备成本、减小占地面积和空间。但是它给变压器的设计和制造也带来了一些技术难题，主要体现在复合短路阻抗计算以及复杂的环流计算，特别是针对多对一和多对多绕组的情况。

常用的复合短路阻抗计算方法有直接计算法、能量法以及场路耦合法等，这些方法各自都有各自的优缺点。① 直接计算法。文献[6]针对某一种多绕组变压器，提出了漏磁路计算方法，但是此方法的应用范围不宽，难以达到广泛通用性。② 能量法。文献[7]先根据磁场能量计算两两绕组短路时的短路阻抗，然后通过矩阵变换计算出变压器的复合短路阻抗，这种方法在求解多个单对单绕组的短路阻抗时，要在有限元软件中进行多次计算，所以增加了求解难度。③ 场路耦合法。在建立变压器模型阶段，需要输入铁心饱和曲线，由于变压器铁心B-H曲线的非线性，使得单次场路计算时间很长；需要每次设定不同的短路工况，来求解某一绕组或某些绕组的复合短路阻抗，总体的工作量与工作时间很长[8-10]。

短路环流计算方面也主要采用场路耦合法，同样是求解时间较长，不同的短路工况下，无法一次求解出所有短路工况下的复合短路阻抗[11-12]。

本文设计了一种基于多绕组变压器简化有限元模型的计算方法，将多绕组变压器进行归类，只要在简化变压器磁场有限元模型得出的电感矩阵基础上，通过输入多绕组变压器的匝数、容量、额定电压参数等，就可以在软件中计算出多绕组变压器的复合短路阻抗及其短路情况下各绕组的电流环流。计算软件基于Matlab GUI 实现，只需要安装个别插件就可以在多种操作系统中运行，具有良好的兼容性；计算方法和软件实用方便，界面操作简单，利于变压器设计人员计算使用。

2 简化变压器磁场有限元模型的复合阻抗及环流计算方法

通过有限元仿真软件可以建立简化的变压器磁场有限元模型，并通过简化的磁场有限元模型推导多绕组变压器两两绕组的短路阻抗。其方法如下：① 将多绕组多匝数变压器简化为多绕组单匝数变压器有限元模型，即把多匝数的绕组看成单匝绕组，建立三维有限元模型；② 通过简化磁场模型得到两两绕组短路阻抗，计算变压器导纳矩阵；③ 根据导纳矩阵与短路工况，计算复合短路阻抗与短路环流。

2.1 简化变压器磁场计算两两短路阻抗

假设原变压器模型的i号线圈匝数为Ni，电流为Ii。首先建立多绕组变压器简化磁场有限元模型，将原来的多匝线圈按照单匝线圈处理，相应的电流由Ii变为NIi，很容易得到单匝线圈的变压器绕组电感矩阵M为

式中，lii为简化的单匝线圈自感；lij为单匝线圈互感。

当i、j号绕组短路处于短路情况时，可根据电感的定义得知

式中，ψii为单匝电感与本身交链的磁链，而在非简化的多匝线圈情况下，当多绕组发生交链时

联立式(1)～(3)可得

同理可得

由短路阻抗的定义可知

将式(4)和式(5)代入式(6)可知

又由于变压器绕组的电阻很小，可以忽略不计，可得

这样就可以通过简化的变压器有限元模型，求解两两绕组的短路阻抗[13-17]。

2.2 求解变压器导纳矩阵

多绕组变压器绕组排布如图1 所示。

图1 多绕组变压器绕组排布

其中，绕组1～m为并联高压绕组，m+1 至n为独立的低压绕组，有不同的电压等级；各绕组匝数分别为N1，N2，…，Nn。

由文献[12]可得

多绕组变压器电压平衡方程为

式中，当i=j时，Xσij为绕组的自漏电抗；i≠j时为i、j两绕组间互感漏电抗。又知

将式(11)、(12)代入式(10)并进行变换可得

将式(9)、(14)联立解得

对式(11)进行逆矩阵变换，联立式(12)可解得

2.3 根据短路工况求解复合短路阻抗与短路环流

多绕组变压器大多采用高压并联供电，低压输电的模式[18-21]。在这种情况下，如绕组编号为i的低压绕组发生短路故障，则其满足短路工作状况的边界条件为：① 各高压绕组电压相等；② 各低压绕组电压为0；③ 短路绕组i电流为额定电流；④ 除i号绕组外，其余低压绕组电流为0。结合短路工况与式(15)，可解得在此工况下的复合短路阻抗与各绕组短路电流

另外还有两种情况：① 高压绕组开路，低压绕组供电；② 高压绕组并联，各低压绕组同时都短路，一个或者多个其他绕组同时短路的情况，均可以列写相应的短路端口条件，代入式(17)求得相应的复合短路阻抗与短路环流[13,17]。

3 复合短路阻抗及环流计算软件与计算实例

根据第2 节的推导，可以得到简化变压器电感矩阵与原变压器两两短路阻抗的关系。通过一系列的矩阵变换与计算，结合各种短路工况下的端口条件，就可以计算多绕组变压器复合短路阻抗及短路环流。中间的算法很复杂，本文通过Matlab GUI 界面开发程序，编制软件来实现相关的计算。

3.1 软件计算流程

按照第2 节的计算公式，可以得到软件计算流程，如图2 所示。软件流程如下：① 通过电感矩阵与绕组数据，计算变压器的导纳矩阵；② 通过供电和短路绕组编号，判断短路工况，得出短路下的端口条件；③ 结合导纳矩阵与端口条件，计算变压器复合短路阻抗与短路环流；④ 输出数据。

图2 软件流程图

3.2 输入输出数据及软件界面

输入输出数据位于软件设计的最上层，其中需要输入的数据如下：绕组个数、绕组编号、绕组匝数、额定电压、绕组容量以及单匝线圈的电感矩阵。其中电感矩阵只需要按照提供的模板进行导入，使得数据提取更加方便快捷。编制的软件界面如图3所示。

图3 软件界面图

多绕组变压器可能存在多种低压，当需要计算的多绕组变压器只存在某一种或几种低压的情况时，需要在剩余的低压绕组个数上填0；计算多个绕组短路工况下的短路阻抗时，只需要输入多个短路绕组的编号，用空格或者逗号隔开，就可以进行相应的计算。填写相关的数据之后，点击“添加电感矩阵并计算”按钮，就可以得到复合短路阻抗与列表形式的各绕组电流。

3.3 实例验证

选取某一8 绕组变压器作为验证对象，其主变压器参数如表1 所示。

表1 变压器参数

其中有4 个绕组为并联的高压绕组，另外4 个绕组为独立的低压绕组。对多绕组变压器复合短路阻抗与环流计算方法进行验证。首先在有限元仿真软件中建立这一多绕组变压器的简化磁场模型，如图4 所示，得到简化单匝电感矩阵m，并计算其电感矩阵M如下

图4 多绕组变压器的简化磁场模型图

再将电感矩阵、变压器基本参数、短路绕组编号作为输入数据导入计算软件；最后得到多绕组变压器的复合短路阻抗与短路环流，如表2所示。

表2 复合短路阻抗和短路环流的计算值与实际测量值以及误差

由表2 可知：① 采用简化变压器磁场有限元模型计算得到的变压器复合短路阻抗计算值与实测值之间的误差，没有超过工程需求的5%，满足工程要求的误差范围；② 不同短路工况情况下，存在不同的短路环流，当变压器高压绕组对某一低压绕组短路时，各高压绕组内的电流是不同的，这主要是由绕组排布不同导致的耦合度不同造成的。

由于实测的多绕组变压器只存在高压与低压两种电压等级，此时只需要中压绕组个数为0，软件可自动识别为无中压绕组。当高压绕组并非串联情况时，可对某一或某多个高压绕组施加电压，填写至“供电绕组编号”框内，即可完成计算。

根据文献[8]中提出的场路耦合计算方法，以及文献[7]中提出的数学计算方法，得到复合短路阻抗与短路环流计算结果如表3 所示。由表3 可知，场路耦合计算方法、数学计算方法与本文所提方法计算结果基本一致，由此验证了本文计算方法的正确性，但文献[8]中场路耦合计算需要结合铁心B-H曲线进行非线性计算，总体计算时间冗长。文献[7]中数学计算方法需要计算两两绕组的短路阻抗，并利用短路阻抗矩阵进行矩阵转换来得到计算结果，计算过程较为复杂，可由矩阵计算软件编制计算程序完成。