简介面向电力系统专业学生与电气工程师的教学PPT系统讲解对称分量法在不对称故障分析中的应用覆盖正序、负序、零序分量的基本概念、序阻抗计算、带负荷发电机的序网络建模以及单相接地、两相短路、两相接地短路等典型故障的分析流程。资源为单个PPT文件大小5.58MB内容精炼紧凑适合配合课程学习或自修使用目前已有348人学习下载。PPT还结合MATLAB函数abc2sc、sc2abc以及lgfault、llfault、dlgfault进行对称分量转换与故障计算同时介绍利用节点阻抗矩阵求解系统响应的思路帮助读者将抽象理论转化为可操作的算法步骤。通过掌握本资源读者可深入理解不对称故障的序网络分析方法提升电力系统暂态分析与故障计算能力。1. 对称分量法不是玄学为什么三相不对称能拆成三个对称网络拿一个实际工况开场10kV 馈线发生单相接地故障相电流飙到几千安另外两相电流相位和幅值全变了。这时候你若还按第九章的三相对称短路去算方程直接没法解因为三相电路已经不对称到连单相等效电路都画不出来。对称分量法解决的就是这件事——它由 Fortescue 在 1918 年提出核心是把一组不对称的三相相量强行分解成三组各自对称的分量正序abc 顺序互差 120°、负序acb 顺序互差 120°、零序三相同相位。这样原本纠缠在一起的三相网络就被拆成三个独立的单相序网络每个序网络单独用单相电路求解最后再叠加回相域。对电力系统继电保护整定、故障分析、负序电流监测来说这是绕不开的基本功。这篇内容基于一份经典的电力系统分析课程 PPTCHAPTER 10我把里面的矩阵推导、序阻抗结论和 MATLAB 故障计算函数串起来讲清楚怎么从 abc 域落到 012 域再一步步算到故障电流和电压适合正在学电力系统分析、或者在用 MATLAB 做故障仿真的工程师。2. 从 abc 到 012对称分量变换矩阵与 MATLAB 的 abc2sc/sc2abc 实操2.1 旋转因子 a 和对称分量变换矩阵的由来定义旋转因子 a 1∠120° -0.5 j0.866。a² 1∠240°a³ 1∠0°。注意性质 1 a a² 0这个性质在后面推导零序分量时非常关键。对于三相不对称电流 Ia、Ib、Ic用 a 相的正序、负序、零序分量表示时满足Ib a²·Ia1 a·Ia2 Ia0Ic a·Ia1 a²·Ia2 Ia0写成矩阵形式就是| Ia | | 1 1 1 | | Ia0 | | Ib | | 1 a² a | | Ia1 | | Ic | | 1 a a² | | Ia2 |中间这个矩阵就是对称分量转换矩阵 A。它的逆矩阵 A⁻¹ 是A⁻¹ (1/3) * | 1 1 1 | | 1 a a² | | 1 a² a |逆矩阵用来从相量反解序分量也就是 abc2sc 函数的底层逻辑。注意系数 1/3这意味着零序分量等于三相电流之和的三分之一。所以当中性点不接地、三相电流之和为 0 时零序电流天然不存在而单相接地故障时故障相电流很大三相电流之和不再是 0零序分量就会出现。2.2 MATLAB 函数 sctm、abc2sc、sc2abc、rec2pol、pol2rec 的参数与用法这套 PPT 里给出的 MATLAB 函数本质是围绕 A 矩阵封装的操作。sctm 定义 A 矩阵输入到工作区后可以直接用 inv 求逆。实际分析中用得最多的是两个转换函数abc2sc(Iabc)输入不对称三相相量输出零序、正序、负序分量直角坐标或极坐标。sc2abc(I012)输入序分量输出三相相量。rec2pol(F)直角坐标转极坐标返回 [幅值, 相角]。pol2rec(F)极坐标转直角坐标。phasor(F)绘制相量图。这些函数支持两种输入格式一种是 (3×1) 的直角坐标复数数组另一种是 (3×2) 的矩阵第一列为幅值第二列为相角度。我一般直接用极坐标输入这样故障类型一眼能看出来。例 10.1 求一组不对称电流的对称分量Ia 1.6∠25°Ib 1.0∠180°Ic 0.9∠132°。直接跑% 例10.1: 不对称电流 - 对称分量 Iabc [1.6 25; 1.0 180; 0.9 132]; % 极坐标输入列1幅值列2相角 I012 abc2sc(Iabc); % 得到直角坐标形式的序分量 I012p rec2pol(I012) % 转成极坐标查看幅值和相角运行结果I012p 0.4512 96.4529 0.9435 -0.0550 0.6024 22.3157三行分别对应零序、正序、负序。这里零序幅值 0.4512 不等于 0说明三相电流之和不为 0系统存在接地回路或不对称。正序分量接近 0.9435∠0°负序也有 0.6024表示不对称程度较高。观察相角可以快速判断故障类型如果是单相接地通常负序与正序幅值关系明显两相短路时零序为 0。例 10.2 是反向变换已知序分量 V012 [0.6∠90°; 1.0∠30°; 0.8∠-30°]求三相电压。% 例10.2: 对称分量 - 三相电压 V012 [0.6 90; 1.0 30; 0.8 -30]; Vabc sc2abc(V012); % 转换到相域 Vabcp rec2pol(Vabc)结果Vabcp 1.7088 24.1825 0.4000 90.0000 1.7088 155.8175可以看到 b 相电压只有 0.4∠90°而 a 相和 c 相幅值相等、相位不对称这就是典型的不对称电压分布。用phasor函数可以直接画出不对称相量图和三个序分量的相量图用来验证边界条件对不对。2.3 零序分量的物理含义中性点接地是分水岭从公式上看零序分量是三相电流的算术平均。实际系统里零序电流只有在有接地回路时才可能流通。比如中性点不接地的 Y 型负荷三相电流之和必然为 0零序电流不存在而中性点经阻抗 Zn 接地时零序电流可以流经 Zn且此时流过 Zn 的电流是 3·I0所以在序网络中中性点阻抗要按 3Zn 计。这个 3 倍关系在变压器零序等效电路里也是核心后面第 4 章会用到。3. 序阻抗输电线、发电机、变压器各自的正负零序差别在哪3.1 Y 型连接负荷的序阻抗矩阵推导PPT 里从 Y 型三相对称负荷入手中性点经 Zn 接地每相自阻抗 Zs相间互感 Zm。列 KVL 方程后得到相域阻抗矩阵 Zabc再利用变换矩阵 A 做相似变换Z012 A⁻¹ · Zabc · A经过矩阵乘法利用 1aa²0 消去互感和中性点阻抗的耦合项最后得到序分量阻抗值零序 Z0Zs 3Zn正序 Z1Zs - Zm负序 Z2Zs - Zm注意这里默认忽略了互感 Zm 对零序的影响。实际上如果考虑互感零序阻抗还包含 -2Zm 项即 Z0 Zs 3Zn - 2Zm。忽略互感后Z012 变成对角阵说明对称分量法最大的好处就是在对称元件上正序、负序、零序三个序网络完全解耦可以单独计算。我自己习惯在 MATLAB 里用例子验证这一步% 求Y型负荷的序阻抗矩阵 Zs j*12; % 自阻抗 Zm j*4; % 互感阻抗 Zn 0; % 中性点阻抗不接地时无穷大先设0看效果 Za [Zs Zm Zm; Zm Zs Zm; Zm Zm Zs]; % Zabc if Zn ~ inf Zn_actual Zn; else Zn_actual 1e6; % 近似开路 end % 实际加中性点阻抗Zm 位置不变对角线加 Zn但相间互感已含 % 更严谨的做法是构造完整 Zabc下面使用常见形式 Zabc [ZsZn_actual Zn_actual Zn_actual; ... Zn_actual ZsZn_actual Zn_actual; ... Zn_actual Zn_actual ZsZn_actual]; % 定义旋转因子 a exp(1j*2*pi/3); A [1 1 1; 1 a^2 a; 1 a a^2]; Z012 inv(A)*Zabc*A; disp(Z012)这段代码里有个容易踩的坑如果直接在 Zabc 对角线加 Zn而互感项不加 Zn那算出来不对。正确做法是中性点阻抗通过 KCL 等效到每相时对角线加 Zn非对角线也加 Zn因为任意两相电流之和会在 Zn 上产生压降影响第三相。上面代码中我用了含 Zn 的全矩阵这样 Z012 的 (1,1) 位置才会正确出现 Zs 3Zn。如果中性点不接地实际上 Zn 趋于无穷序网络里零序阻抗就变成无穷大对应零序网络开路。3.2 输电线和发电机负序等于正序是常态零序看大地回路输电线路是静止元件正序和负序电流流过的几何排列完全一样所以 Z1 Z2。零序则不同三相零序电流同相位必须通过大地和架空地线构成回路。PPT 里用等边三角形排列的导线推导单位长度零序电抗可以写成X0 X1 3·Xn其中 Xn 是考虑大地回路影响的等效电抗包含 Carson 公式里的土壤电阻率和等值深度项。所以输电线路的零序阻抗通常比正序阻抗大得多经验值是 X0 ≈ 3·X1 左右具体要看大地电阻率。同步发电机是旋转设备正序电抗随故障瞬间不同而变化次暂态电抗、暂态电抗、稳态电抗都可以作为正序电抗取决于你算的是故障初瞬还是稳态。负序电抗的物理图像是负序电流产生的旋转磁场相对转子以两倍同步速反转励磁绕组对负序磁场没有贡献阻尼绕组只在交轴起作用所以工程上近似取负序电抗等于正序次暂态电抗 X2 ≈ Xd。零序电抗更小因为零序电流产生的气隙合成磁动势为零没有电枢反应只有漏磁通所以 X0 ≈ 漏电抗。这里给一个对比表方便记忆元件正序 Z1负序 Z2零序 Z0输电线路与几何排列相关等于 Z1大于 Z1含大地回路同步发电机Xd/Xd/Xd约等于 Xd约等于漏抗变压器漏阻抗等于 Z1取决于接线和中性点接地等于漏阻抗或开路3.3 变压器零序等效电路绕组接法和中性点接地共同决定通路变压器是静态设备漏阻抗与相序无关所以 Z1 Z2 Z0 漏阻抗。但是变压器能否流通零序电流完全由绕组连接方式决定。PPT 图 10.6 给了几种典型接线Y-Y 两侧中性点都接地两侧都有零序通路零序等值电路就是每相漏阻抗。Y-Y 仅一次侧中性点接地二次侧中性点不接地二次侧零序电流无法流通一次侧也被迫开路等效为无穷大阻抗。Y-△ 一次侧中性点接地一次侧 Y 接地形成零序回路△ 侧内部环流但出线上没有零序电流等效电路是一次侧漏阻抗串联△ 侧短路对零序而言。Y-△ 一次侧中性点不接地无通路等效开路。△-△零序电流只能在△内部环流外部线路无零序通路。做故障计算时第一步就是根据变压器接线画零序网络。比如系统里常见的 Yd11 变压器高压侧中性点接地那么接地故障时零序网络高压侧可以流通低压侧△绕组只是提供内部环流不会把零序传到低压线路。4. 不对称故障怎么算单相接地、两相短路、两相接地短路的序网连接与函数调用4.1 序网络连接规则串联与并联的物理本质对称分量法把不对称故障的边界条件转换到序域后正、负、零序网络会在故障点按特定方式连接。以故障点作为端口三个序网络等效成三个戴维南电路每相只有一个串联阻抗从故障点看进去。然后根据故障类型单相接地A 相接地边界条件为 Ua0IbIc0。变换到序域后得到 I0I1I2且 U0U1U20。这意味着三个序网络在故障点串联因为流过同一电流电压相加为零。两相短路B、C 相短路边界条件为 Ib-IcUa 任意Ia0。得到 I00I1-I2U1U2。正序和负序网络在故障点并联零序不参与。两相接地短路B、C 相接地边界条件为 UbUc0Ia0。得到 I00 不成立而是 U0U1U2且 I0I1I20。三个序网络在故障点并联。记忆方法接地故障通常涉及零序网络两相短路没有接地零序电流为零。单相接地是串联两相接地是并联两相短路是正负序并联。4.2 lgfault、llfault、dlgfault 函数参数拆解PPT 给出三个故障计算函数分别对应单相接地、两相短路、两相接地。它们的输入参数是% 单相接地 lgfault(zdata0, zbus0, zdata1, zbus1, zdata2, zbus2, V) % 两相短路 llfault(zdata1, zbus1, zdata2, zbus2, V) % 两相接地短路 dlgfault(zdata0, zbus0, zdata1, zbus1, zdata2, zbus2, V)参数含义zdata0、zdata1、zdata2零序、正序、负序网络的支路数据格式通常是 [起始节点 终止节点 阻抗值]用于构成节点阻抗矩阵。zbus0、zbus1、zbus2对应的节点阻抗矩阵可以由 zdata 用zbus函数生成也可以自己用inv求节点导纳矩阵的逆。V故障前节点电压向量通常是所有节点电压均为 1∠0°标幺值也可以考虑负荷时用潮流结果。函数内部会先把网络变换到故障端口等效阻抗再按故障类型连接序网最后输出故障点各序电流电压、各节点电压和各支路电流。实际使用时我一般先把系统参数按序网分别整理成三个支路数据表再调用zbus生成矩阵。下面是一个伪代码流程% 构造简单三节点系统的序网数据 zdata1 [1 2 0.1j; 2 3 0.2j; 1 3 0.3j]; % 正序支路 zdata2 zdata1; % 负序网络与正序相同静态元件 zdata0 [1 2 0.3j; 2 3 0.5j]; % 零序只有接地回路支路 zbus1 zbus(zdata1); % 生成正序节点阻抗矩阵 zbus2 zbus(zdata2); zbus0 zbus(zdata0); V [1; 1; 1]; % 故障前各节点电压标幺值 % 计算单相接地故障 lgfault(zdata0, zbus0, zdata1, zbus1, zdata2, zbus2, V);注意这里zdata里的阻抗是各元件的序阻抗不是相阻抗。比如变压器正序阻抗用漏抗零序阻抗要根据接线方式决定是漏抗还是无穷大无穷大就不画该支路。发电机零序阻抗很小但只有中性点接地时才接入零序网络。4.3 例 10.3网孔法 vs 对称分量法同一个问题的两条路PPT 例 10.3 是一个对称三相系统相电压 100VY 型对称负荷中性点不接地每相串联电抗 Zs j12Ω相间互感 Zm j4Ω。要求线电流。这个题目的妙处在于负荷中性点不接地意味着零序电流为 0但相间互感让三相之间耦合直接解网孔方程很麻烦而对称分量法刚好解耦。网孔法代码PPT 给的思路% 例10.3 (a) 网孔法 Zs j*12; Zm j*4; VL 100*sqrt(3); % 两个网孔方程加 KCL Z [(Zs-Zm) -(Zs-Zm) 0; ... 0 (Zs-Zm) -(Zs-Zm); ... 1 1 1]; V [VL*cos(pi/6)j*VL*sin(pi/6); ... VL*cos(-pi/2)j*VL*sin(-pi/2); 0]; I inv(Z)*V;对称分量法则简单得多因为负荷对称且中性点不接地正序阻抗 Z1 Zs - Zm j8Ω负序阻抗 Z2 与正序相同零序阻抗为无穷大。正序电压取线电压的相量形式比如 Ua 100∠0°那么正序电流I1 Ua / Z1 100 / (j8) -j12.5 A负序电流为零电源对称零序电流为零。线电流就是正序电流。这个结果和网孔法求出来一致但计算量差了一个矩阵求逆。对比下来就能理解为什么工程上偏爱对称分量法——正常运行时系统三相对称故障是对称性被破坏的唯一来源把故障用序网拆开后大量正常工作下的对称关系可以直接复用。5. 对称分量法排错与进阶从相量图验证结果到节点阻抗矩阵程序的参数传递5.1 用 phasor 画图检查转换是否合理MATLAB 工具包里phasor函数能同时画出原始三相相量和分解后的零、正、负序相量图。我自己调试故障计算时习惯先跑一个已知的对称系统如果不加故障abc2sc出来的负序和零序应当接近 0。如果负序或零序不为 0说明数据输入有问题比如三相电压幅值不相等、相序错乱或者坐标转换时相角单位弄错一个是度一个是弧度。rec2pol和pol2rec会自动处理单位但手工输入时一定要确认第二列是度。另外画图时注意 x 轴和 y 轴比例零序分量很多时候很小直接看数值比看图更准确。5.2 节点阻抗矩阵在不对称故障程序里的角色PPT 10.8 节专门讲用节点阻抗矩阵分析不对称故障。节点阻抗矩阵 Zbus 的物理意义是Zbus 的第 k 列元素是在节点 k 注入单位电流时各节点电压。对不对称故障我们只需要在故障点 f 处的三个序网络等效阻抗Zff0、Zff1、Zff2它们就是各序节点阻抗矩阵在故障点对角线元素。因为故障点可能不是发电机节点也不一定在支路端点上所以程序内部要先计算故障点所在支路的位置然后用支路参数把故障点等效为新增节点再对 Zbus 做修正。修正方法有两种一是支路追加法把故障点当作新节点用zbus原函数扩展矩阵二是通过故障支路阻抗和两端节点阻抗利用公式直接算出 Zff。PPT 里lgfault等函数接收 zdata 和 zbus 两个参数目的就是允许程序内部根据 zdata 判断故障点位置再结合 zbus 计算端口等值阻抗。给一个故障点端口阻抗计算的常见套路% 假设故障在支路 i-j 的中点支路阻抗 zij % 已知 zbus 的 i、j 列元素 Zii zbus(i,i); Zjj zbus(j,j); Zij zbus(i,j); % 故障点距 i 的比例系数 alpha0~1 alpha 0.5; % 新增节点 f 的自阻抗 Zff Zii alpha^2*(Zjj - 2*Zii Zij) 2*alpha*(Zij - Zii) alpha*(1-alpha)*zij; % 简化式实际按支路追加法推导注意这个公式是经过简化的真实推导要结合支路追加公式。关键思想是故障点不是天然节点必须把支路分割再合成新的等值阻抗。这个细节在课程设计或工程计算里很容易被忽略很多同学直接拿原 zbus 的对角线用算出来故障电流偏大因为没有计入故障点到两端母线之间的支路阻抗分压。5.3 程序参数传递的三个容易错的地方检查自己或学生的程序时我发现三个高频错误。第一个是zbus的参数顺序zdata 里支路阻抗应该用标幺值且所有序网的基准容量和基准电压必须一致。如果正序网络用 100MVA 基准零序网络用 50MVA 基准算出来 Zff 完全乱套。第二个是变压器零序回路的处理中性点经电抗接地的变压器零序等值电路里要串入 3Zn而不是 Zn很多人漏了 3 倍系数。第三个是故障类型对应的序网连接方式两个函数虽然在参数上一样但内部连接逻辑不同不要试图通过改参数把llfault算成两相接地。判断方法很简单看故障后零序电流如果零序电流为 0说明是两相短路如果零序电流不为 0 且正序负序相等那是单相接地如果三个序电流都不为 0 且电压相等那是两相接地。调试时建议先在简单两节点系统上验证比如一个电源经线路带负荷做两相短路手算预期结果与程序输出对比。手算公式两相短路时故障相电流幅值 1.732 × E / (Z1 Z2)其中 E 是故障前相电压。如果程序输出一致再放心去算复杂网络这样能快速定位问题是出在序网参数还是连接逻辑上。本文还有配套的精品资源点击获取