简介这是一份面向电力系统专业学生与工程技术人员的《高等电力系统分析》课件聚焦阻抗模拟这一核心方法重点讲解如何通过 Ward 等值、阻抗/导纳矩阵以及规范阻抗形式来模拟不对称故障。内容从基础概念入手覆盖 A 相接地、BC 相间短路、BC 两相接地等典型故障场景并给出等值阻抗通用公式与支路内部故障的规范化处理思路对学习故障分析、保护配置和系统等效建模很有帮助。课件为单个 PPT 文件大小仅 359KB便于下载和课堂演示作者为 wxg520cxl。目前已有 75 人浏览学习适合作为电力系统课程复习、考研备考或教学备课的补充资料。资源虽小但知识点密度较高包含 Ward 等值分解、导纳阵构造、10-6/10^6 数值近似等关键技巧可帮助读者快速梳理阻抗模拟的完整脉络。1. 阻抗模拟课件电力系统分析里最容易被忽略的一层功夫初学电力系统分析时很多人把不对称故障计算理解成背公式记住A相接地、BC相间短路、两相接地的边界条件套对称分量法就结束。真正动手做暂态仿真、保护定值校核之后才会发现工程程序里根本不这么干而是把所有故障统一表达成阻抗/导纳网络交矩阵求解器去处理。这套《高等电力系统分析》课件正好覆盖这个工程缺口。它把Ward等值、不对称故障规范阻抗形式、支路内部故障新增节点再消去这三件事串成一条完整的计算链课件里反复出现的“10⁻⁶代表短路、10⁶代表开路”恰恰是故障模拟的数值基础。适合正在做故障计算、保护整定或者临时需要把教材公式改写成程序的人这份课件不是给你补概念而是告诉你故障计算在程序里到底长什么样。2. Ward等值从分块方程组到等值网络的完整推导2.1 分块方程A X B保留系统与外部系统的第一次拆分课件里Ward等值的起点是一个标准的线性方程组形式是A X B。这个A不是最终要用的节点导纳阵而是与网络方程同构的系数矩阵你要算的电压或者电流状态量都在X里。把节点按“保留/消去”分成两组之后方程写成块状[ \begin{bmatrix} A_{11} A_{12} \ A_{21} A_{22} \end{bmatrix} \begin{bmatrix} X_1 \ X_2 \end{bmatrix}\begin{bmatrix} B_1 \ B_2 \end{bmatrix} ]其中X₁是保留节点研究区域加边界节点的状态量X₂是外部节点的状态量。第一行展开是A₁₁X₁ A₁₂X₂ B₁移项后得到课件中的第一个关键公式X₁ A₁₁⁻¹B₁ − A₁₁⁻¹A₁₂X₂到这里先别急着往代码里搬两个量的物理含义值得停下来想清楚。A₁₁⁻¹B₁这一项含义是外部系统完全去掉、只保留系统自身时的解课件里一般记为X₁的零阶近似A₁₁⁻¹A₁₂这一项反映的是外部系统通过边界节点对保留系统的耦合影响记作W₁₂。于是整个式子变成X₁ X₁⁽⁰⁾ − W₁₂X₂这是Ward等值在阻抗域最精炼的表达保留节点的响应等于本系统自身响应减去外部系统经边界注入的修正量。工程上做外部系统等值真正要保留的就是这个W₁₂和等值后的边界注入。2.2 消去外部节点等值导纳的舒尔补形式接着看第二行A₂₁X₁ A₂₂X₂ B₂。把上一节X₁的表达式代进去整理之后可以得到课件中给出的外部节点解X₂ (A₂₂ − A₂₁A₁₁⁻¹A₁₂)⁻¹(B₂ − A₂₁A₁₁⁻¹B₁)括号里那个A₂₂ − A₂₁A₁₁⁻¹A₁₂就是数值线性代数里的舒尔补。高斯消去时先消X₁剩下的系数矩阵就是它。在Ward等值里它的角色是消去外部系统之后边界节点看到的等值网络参数。如果A是节点导纳阵那这一项就是边界节点的等值导纳修正如果A是回路阻抗阵那这一项就是等值阻抗。课件后面所有等值阻抗通用公式本质上都在反复用这个舒尔补结构。再用分块规模算一下开销设全系统n个节点保留r个节点外部n−r个。直接解全系统是O(n³)量级做Ward等值后求A₁₁⁻¹B₁和W₁₂的开销大约在O(r³)加上一次性O(r²(n−r))。外部节点越多等值省的计算量越可观这正是多馈入系统、大规模互联电网用Ward等值的原因。但注意这个等式成立的前提是A₁₁可逆也就是保留节点之间不能出现纯开路导致的对角为零这一点在导纳阵装配时就要预防。2.3 数值算例两个节点手算验证等值前后响应一致空讲公式容易变成黑匣子我习惯用一个最小算例把手推过程走一遍。假设一个两节点网络系数矩阵取A [[2, −1], [−1, 3]]右端项B [1, 0]把节点1当保留节点、节点2当外部节点。先直接解原方程组由第二行可得X₂ X₁/3代入第一行2X₁ − X₁/3 1解出X₁ 0.6完整解是X₁ 0.6、X₂ 0.2。再看Ward等值怎么做。这里A₁₁ [2]A₁₂ [−1]A₂₁ [−1]A₂₂ [3]。外部节点消去后的等值系数是A₁₁′ A₁₁ − A₁₂A₂₂⁻¹A₂₁ 2 − (−1)×(1/3)×(−1) 5/3右端项B₁′ B₁ − A₁₂A₂₂⁻¹B₂ 1。于是等值后只解一个一元方程X₁ 1/(5/3) 0.6和全系统直接解得出的X₁完全一致。这个算例虽小但说明了Ward等值消去外部节点不改变保留节点的电压/电流响应代价完全体现在边界参数修正上这也是“等值”这两个字的数学依据。实际程序里你可能会用更大的矩阵但验证思路完全一样先全系统解一遍再等值解一遍两个结果对比误差在数值精度内等值就算通过。2.4 Ward等值的适用边界什么系统适合做等值多次用下来我的体会是Ward等值不是万能的边界条件比较苛刻。第一外部系统最好以恒阻抗线性元件为主负荷如果当成恒功率注入运行点一变等值注入就失真。第二研究区域和外部区域之间的联络线越少、越集中等值精度越高边界节点应该选在电气距离近、电压支撑强的母线上。第三如果外部系统里含有需要重点观察的机组或故障点别把它们消掉否则故障后暂态过程根本等值不出来。课件里的Ward等值部分适合用来理解原理和做潮流层面的化简真正做暂态稳定分析时还是建议把研究区域附近的边界尽量划大一点宁多留几个节点也不要贪图计算速度把边界划得过分干净。3. 不对称故障的规范阻抗形式三种故障的装配与10⁻⁶/10⁶约定3.1 三种典型不对称故障的阻抗接线表课件第二部分直接用阻抗接线的方式画出了三种不对称故障A相接地、BC相间短路、BC两相接地。A、B、C三相各自通过阻抗Z_A、Z_B、Z_C连接到中性点NZN是中性点对地阻抗。三种故障归一化之后区别只在于三相各自的阻抗值怎么取故障类型ZAZBZC说明A相金属性接地10⁻⁶10⁶10⁶A相经零阻抗接地B、C相开路A相经阻抗接地ZN10⁶10⁶ZN为过渡电阻常见于弧光接地BC相间短路不接地10⁶ZbcZbcB、C相间直接连接对地开路BC相间接地短路10⁶ZNZNB、C经阻抗接地A相开路这张表的核心是“规范”两个字。无论故障发生在哪一相、是否带过渡电阻都用同一套三相阻抗装配流程程序不需要为每一种故障单独写分支。实际建模时先按表格把每相的阻抗值确定下来再换算成导纳填进矩阵顺手把故障类型名也留在注释里。A相经阻抗接地那一行的ZN如果取0就是金属性接地如果取实际过渡电阻弧光接地、树枝碰线这类带电阻故障也能直接模拟不需要改程序结构。3.2 10⁻⁶与10⁶短路与开路的数值化约定课件里专门强调了“采用10⁻⁶代表0、采用10⁶代表∞”这个约定在数值计算里几乎是必须的原因有两个。第一个原因是避免除零。理想短路阻抗是0直接放在分母上做导纳运算导纳就是无穷大矩阵元素直接爆掉理想开路的阻抗是无穷大导纳是0对角元素为0又会把矩阵推成奇异。10⁻⁶和10⁶这两个数量级是工程上通用的替代值短路导纳高达10⁶开路导纳低到10⁻⁶数值上既不溢出、又不奇异计算结果和理想值之间的误差在6个数量级以下对工程精度完全够。第二个原因是让所有故障统一进一个装配流程。不管你是短路还是开路在程序里都是同一个填表函数只是那个位置的阻抗值写10⁻⁶还是10⁶的区别。这个约定的坑在于方向容易记反10⁻⁶是给短路支路的10⁶是给开路支路的。我见过不止一次同事把故障相的阻抗写成10⁶结果A相接地故障算出来三相电流几乎为零误以为保护拒动其实是短路支路被填成了开路。后面第5章排查章节会再展开这里先记住方向。3.3 从阻抗装配到三相导纳阵程序里真正写入的矩阵阻抗表只是人看的程序算的是导纳阵。以一个单节点三相网络为例节点导纳阵按相展开排布规则是对角元素相加、非对角元素取负对角元素Yₐₐ′是A相所有并联导纳之和包括A相对地导纳、A相上的并联电抗器等非对角元素−Yₐᵦ是A、B相之间的互导纳取负号是因为互导纳的电流方向和自导纳相反。以A相金属性接地为例装配时A相的对地阻抗填10⁻⁶导纳是10⁶B、C相对地阻抗填10⁶导纳是10⁻⁶。装配完成之后的Y矩阵代入I YU解出来就是各相的电压和电流。程序里往往是先算I YU再对需要求电压的节点做U Y⁻¹I两步都绕不开这个装配好的矩阵。很多初学的人会在这一步卡住为什么不先画故障网络再手推公式答案很简单手推公式只适用于教材里的三种对称故障场景一旦出现带过渡电阻、非线性电弧或者多故障点并存手推的边界条件就崩了而导纳阵装配完全不受影响你只需要改对应位置的阻抗值。课件里“用规范的阻抗形式模拟不对称故障”这一整节落脚点就是把这个矩阵装配流程讲透后面所有等值阻抗计算都建立在它之上。4. 等值阻抗通用公式与支路内部故障新增节点、消去节点的一步到位4.1 三相导纳阵求逆等值阻抗通用公式的落地课件里等值阻抗通用公式的推导建立在三相导纳阵的基础上。前面装配出的Y阵求逆之后得到的就是等值阻抗阵Z Y⁻¹。对3×3矩阵程序里可以直接调矩阵求逆手算时用伴随矩阵法Z adj(Y) / det(Y)其中det(Y)对于三相系统可以展开成三项乘积的差。具体展开式不一定要背程序也不会让你手算行列式但有两个物理解读值得记住。对角元素Zₐₐ表示A相看进去的等值自阻抗也就是A相注入单位电流时A相的电压响应非对角元素Zₐᵦ表示A、B两相之间的互阻抗它的存在说明三相网络在故障后不再解耦这正是对称分量法在不对称故障时要逐序分解的原因。课件里把这个求逆过程称作等值阻抗通用公式想表达的是任一三相网络只要Y阵装配正确等值阻抗就由Y⁻¹唯一确定不需要针对故障类型手动构造等值电路。程序实现时直接调numpy.linalg.inv或者Matlab的inv()但验证结果时还是要回到伴随矩阵的物理解读上来防止矩阵奇异时程序返回垃圾结果你还不知道。4.2 支路内部故障的规范化处理新增节点把故障位置变成普通节点课件最后一部分处理的是线路内部故障。一段线路i–j中间距离i点不远处f点发生故障常见做法是把f点新增为节点原线路被拆成两段i–f段导纳Y_L1、f–j段导纳Y_L2f点经故障阻抗Z_F连接到地。节点导纳阵从2×2扩成3×3排布如下节点ifjiYii Y_L1−Y_L10f−Y_L1Y_L1 Y_L2 Y_F−Y_L2j0−Y_L2Yjj Y_L2其中Y_F 1/Z_F是故障点对地导纳金属性接地时Z_F取10⁻⁶Y_F就是10⁶。这里Y_L1和Y_L2不是随便取的如果线路i–j总导纳是Y_L故障点距i端比例为p那么Y_L1 Y_L/p、Y_L2 Y_L/(1−p)看着反直觉但物理上是对的线路段越短导纳越大。这个p就是故障位置参数程序里只要改p就能扫出沿线任意位置故障的电气量曲线这在保护定值校核里非常实用。4.3 消去故障节点后的网络合并把故障拍平回原来的节点规模f节点在物理上是真实的故障点但网络求解时它没有注入电流属于无源节点。把f节点从方程组里消去实际上是对f所在行做高斯消去消去之后i、j两点之间出现新的转移导纳原网络节点规模恢复。按3×3矩阵手动消去第二行第二列可以得到i–j转移导纳变为−Y_L1·Y_L2 / (Y_L1 Y_L2 Y_F)i点对地导纳修正Yii Y_L1 − Y_L1² / (Y_L1 Y_L2 Y_F)j点对地导纳修正Yjj Y_L2 − Y_L2² / (Y_L1 Y_L2 Y_F)物理解读很直观故障点引入的导纳Y_F通过分压作用落到i、j两端改变了原线路的等值对地导纳和互导纳。如果Y_F特别大金属性接地i–j转移导纳趋近于−Y_L1·Y_L2/(Y_L1Y_L2)相当于线路被故障点“劈”成了两段并联如果Y_F接近开路修正项趋近于0原线路保持不变。这套做法最大的好处是规范化不管是节点故障还是支路内部任意位置的故障最终都走同一条新增节点→装配导纳阵→消去无源节点→求逆流程故障位置参数化之后程序里只需改p和Y_F不需要重写整个网络方程。5. 阻抗模拟常见问题排查四个翻车点与对应处理5.1 10⁻⁶和10⁶的方向填反故障直接算成开路现象单相接地故障算出来的故障电流几乎为0三相电压全部接近正常保护逻辑怎么校都校不上。原因程序里把故障相的阻抗填成了10⁶金属性短路被模拟成了开路故障电流自然出不来。解决装配时统一先写阻抗表再转导纳阵给短路和开路各定义一个常量比如Z_SHORT 1e-6、Z_OPEN 1e6注释写明“短路用1e-6、开路用1e6”不要在装配函数里直接写裸数字。所有故障类型都从这个常量表取值方向就不会错。5.2 开路相对角导纳填零矩阵奇异求逆失败现象Y阵求逆报singular或者算出来的电压出现inf/nan程序直接崩。原因开路相的对地导纳填了0对角元素为零矩阵不可逆。这是最容易忽略的细节短路相用了10⁻⁶觉得开路相“不接东西”就直接写0结果矩阵奇异。解决开路相对地导纳也必须有值用10⁻⁶量级替代0。数值上是“漏电”的开路物理效果接近理想开路但矩阵保持非奇异。这就是课件里10⁻⁶约定最重要的工程含义不只短路要避零开路也要兜底。5.3 Ward等值的外部负荷用恒阻抗近似运行点一变误差失控现象潮流算例里等值前后的损耗对不上暂态仿真时边界母线电压偏差超过工程容忍度。原因等值时把外部系统负荷当恒阻抗吸收进去但实际负荷是恒功率特性运行点变了恒阻抗模型的注入就不对了等值精度随之劣化。解决把外部系统中对边界影响大的发电机节点保留下来或者用诺顿等值注入替代纯阻抗等值再进一步的做法是给边界节点配上可调的等值注入随运行点修正。别为了省计算量把边界划得太小留两三个关键机组比堆几百个远端节点有用得多。5.4 课件公式符号和教材对不上代码也跟着写错现象照着课件公式写程序解出来的电压符号反了或者矩阵维度对不上debug半天找不到根源。原因课件用A X B的分块形式表述Ward等值符号体系和常用的节点导纳阵Y U I不一样A可能是导纳阵也可能是阻抗阵B的符号约定各套教材不一致直接照抄容易抄错。解决动手写代码前先用最小两节点算例手推一遍把A₁₁、A₁₂、B₁逐个对上再落到代码。我一般会把课件每个公式在笔记本上对照教材重写一遍标出符号含义表再开写这一步能省掉后面大量debug时间尤其是涉及分块矩阵维度时提前标清楚每一块的尺寸比事后猜快得多。6. 进阶验证用一个最小Python脚本把课件公式跑通6.1 最小验证脚本课件看到这个程度最有效的验证方式是把A相接地故障重新实现一遍。我用Python演示最小流程装配三相导纳阵、求逆得等值阻抗、解出故障电流三步对应课件前三节的核心内容。import numpy as np Z_SHORT 1e-6 # 短路阻抗模拟0 Z_OPEN 1e6 # 开路阻抗模拟无穷大 # A相金属性接地A相短路接地B/C相开路 z_a Z_SHORT z_b Z_OPEN z_c Z_OPEN # 三相单节点导纳阵只含对地导纳无相间耦合 y np.array([ [1/z_a, 0, 0], [0, 1/z_b, 0], [0, 0, 1/z_c] ], dtypecomplex) # 等值阻抗阵 导纳阵求逆 z_eq np.linalg.inv(y) # 施加A相单位电压B/C相无注入求解故障电流 u np.array([10j, 0, 0]) i y u print(等值自阻抗 Z_AA , abs(z_eq[0, 0])) print(故障电流 I_A , abs(i[0]), 标幺值1e6量级说明金属性短路)这段代码对应课件里A相接地那一页的完整计算链。参数说明Z_SHORT是金属性短路的数值替身Z_OPEN是开路的数值替身y矩阵的对角元素分别对应A、B、C相对地导纳求逆后z_eq的对角元素就是各相等值自阻抗故障电流直接用导纳阵乘电压得到。对于BC相间短路只需要把z_b和z_c都改成Z_SHORTy矩阵改成含相间耦合的形式其余逻辑完全不用动。6.2 把课件公式转成自己的笔记与汇报材料最后说一个课件整理的技巧。这份PPT里公式密集直接在屏幕上反复翻页效率很低。我习惯把每一页关键公式截图导出成高清图片再插进自己的OneNote或汇报文档里旁边用手写批注写清参数含义。导出时记得把页面分辨率调高别用默认压缩否则公式边缘全是锯齿放大就糊了。整理完再对照第5章那几个坑做一轮自查基本就能把课件的计算链吃透。从那以后我每拿到一套矩阵公式密集的课件都强制走同一个流程先按故障类型填阻抗表再核对短路/开路量级方向最后用最小脚本复算一遍。这套流程走下来再回自己的仿真工程你会明显感觉到故障计算不是背公式而是搭矩阵、填参数、求逆三步走。希望帮到你。本文还有配套的精品资源点击获取 从Ward等值到阻抗模拟：电力系统故障计算的工程实现
简介这是一份面向电力系统专业学生与工程技术人员的《高等电力系统分析》课件聚焦阻抗模拟这一核心方法重点讲解如何通过 Ward 等值、阻抗/导纳矩阵以及规范阻抗形式来模拟不对称故障。内容从基础概念入手覆盖 A 相接地、BC 相间短路、BC 两相接地等典型故障场景并给出等值阻抗通用公式与支路内部故障的规范化处理思路对学习故障分析、保护配置和系统等效建模很有帮助。课件为单个 PPT 文件大小仅 359KB便于下载和课堂演示作者为 wxg520cxl。目前已有 75 人浏览学习适合作为电力系统课程复习、考研备考或教学备课的补充资料。资源虽小但知识点密度较高包含 Ward 等值分解、导纳阵构造、10-6/10^6 数值近似等关键技巧可帮助读者快速梳理阻抗模拟的完整脉络。1. 阻抗模拟课件电力系统分析里最容易被忽略的一层功夫初学电力系统分析时很多人把不对称故障计算理解成背公式记住A相接地、BC相间短路、两相接地的边界条件套对称分量法就结束。真正动手做暂态仿真、保护定值校核之后才会发现工程程序里根本不这么干而是把所有故障统一表达成阻抗/导纳网络交矩阵求解器去处理。这套《高等电力系统分析》课件正好覆盖这个工程缺口。它把Ward等值、不对称故障规范阻抗形式、支路内部故障新增节点再消去这三件事串成一条完整的计算链课件里反复出现的“10⁻⁶代表短路、10⁶代表开路”恰恰是故障模拟的数值基础。适合正在做故障计算、保护整定或者临时需要把教材公式改写成程序的人这份课件不是给你补概念而是告诉你故障计算在程序里到底长什么样。2. Ward等值从分块方程组到等值网络的完整推导2.1 分块方程A X B保留系统与外部系统的第一次拆分课件里Ward等值的起点是一个标准的线性方程组形式是A X B。这个A不是最终要用的节点导纳阵而是与网络方程同构的系数矩阵你要算的电压或者电流状态量都在X里。把节点按“保留/消去”分成两组之后方程写成块状[ \begin{bmatrix} A_{11} A_{12} \ A_{21} A_{22} \end{bmatrix} \begin{bmatrix} X_1 \ X_2 \end{bmatrix}\begin{bmatrix} B_1 \ B_2 \end{bmatrix} ]其中X₁是保留节点研究区域加边界节点的状态量X₂是外部节点的状态量。第一行展开是A₁₁X₁ A₁₂X₂ B₁移项后得到课件中的第一个关键公式X₁ A₁₁⁻¹B₁ − A₁₁⁻¹A₁₂X₂到这里先别急着往代码里搬两个量的物理含义值得停下来想清楚。A₁₁⁻¹B₁这一项含义是外部系统完全去掉、只保留系统自身时的解课件里一般记为X₁的零阶近似A₁₁⁻¹A₁₂这一项反映的是外部系统通过边界节点对保留系统的耦合影响记作W₁₂。于是整个式子变成X₁ X₁⁽⁰⁾ − W₁₂X₂这是Ward等值在阻抗域最精炼的表达保留节点的响应等于本系统自身响应减去外部系统经边界注入的修正量。工程上做外部系统等值真正要保留的就是这个W₁₂和等值后的边界注入。2.2 消去外部节点等值导纳的舒尔补形式接着看第二行A₂₁X₁ A₂₂X₂ B₂。把上一节X₁的表达式代进去整理之后可以得到课件中给出的外部节点解X₂ (A₂₂ − A₂₁A₁₁⁻¹A₁₂)⁻¹(B₂ − A₂₁A₁₁⁻¹B₁)括号里那个A₂₂ − A₂₁A₁₁⁻¹A₁₂就是数值线性代数里的舒尔补。高斯消去时先消X₁剩下的系数矩阵就是它。在Ward等值里它的角色是消去外部系统之后边界节点看到的等值网络参数。如果A是节点导纳阵那这一项就是边界节点的等值导纳修正如果A是回路阻抗阵那这一项就是等值阻抗。课件后面所有等值阻抗通用公式本质上都在反复用这个舒尔补结构。再用分块规模算一下开销设全系统n个节点保留r个节点外部n−r个。直接解全系统是O(n³)量级做Ward等值后求A₁₁⁻¹B₁和W₁₂的开销大约在O(r³)加上一次性O(r²(n−r))。外部节点越多等值省的计算量越可观这正是多馈入系统、大规模互联电网用Ward等值的原因。但注意这个等式成立的前提是A₁₁可逆也就是保留节点之间不能出现纯开路导致的对角为零这一点在导纳阵装配时就要预防。2.3 数值算例两个节点手算验证等值前后响应一致空讲公式容易变成黑匣子我习惯用一个最小算例把手推过程走一遍。假设一个两节点网络系数矩阵取A [[2, −1], [−1, 3]]右端项B [1, 0]把节点1当保留节点、节点2当外部节点。先直接解原方程组由第二行可得X₂ X₁/3代入第一行2X₁ − X₁/3 1解出X₁ 0.6完整解是X₁ 0.6、X₂ 0.2。再看Ward等值怎么做。这里A₁₁ [2]A₁₂ [−1]A₂₁ [−1]A₂₂ [3]。外部节点消去后的等值系数是A₁₁′ A₁₁ − A₁₂A₂₂⁻¹A₂₁ 2 − (−1)×(1/3)×(−1) 5/3右端项B₁′ B₁ − A₁₂A₂₂⁻¹B₂ 1。于是等值后只解一个一元方程X₁ 1/(5/3) 0.6和全系统直接解得出的X₁完全一致。这个算例虽小但说明了Ward等值消去外部节点不改变保留节点的电压/电流响应代价完全体现在边界参数修正上这也是“等值”这两个字的数学依据。实际程序里你可能会用更大的矩阵但验证思路完全一样先全系统解一遍再等值解一遍两个结果对比误差在数值精度内等值就算通过。2.4 Ward等值的适用边界什么系统适合做等值多次用下来我的体会是Ward等值不是万能的边界条件比较苛刻。第一外部系统最好以恒阻抗线性元件为主负荷如果当成恒功率注入运行点一变等值注入就失真。第二研究区域和外部区域之间的联络线越少、越集中等值精度越高边界节点应该选在电气距离近、电压支撑强的母线上。第三如果外部系统里含有需要重点观察的机组或故障点别把它们消掉否则故障后暂态过程根本等值不出来。课件里的Ward等值部分适合用来理解原理和做潮流层面的化简真正做暂态稳定分析时还是建议把研究区域附近的边界尽量划大一点宁多留几个节点也不要贪图计算速度把边界划得过分干净。3. 不对称故障的规范阻抗形式三种故障的装配与10⁻⁶/10⁶约定3.1 三种典型不对称故障的阻抗接线表课件第二部分直接用阻抗接线的方式画出了三种不对称故障A相接地、BC相间短路、BC两相接地。A、B、C三相各自通过阻抗Z_A、Z_B、Z_C连接到中性点NZN是中性点对地阻抗。三种故障归一化之后区别只在于三相各自的阻抗值怎么取故障类型ZAZBZC说明A相金属性接地10⁻⁶10⁶10⁶A相经零阻抗接地B、C相开路A相经阻抗接地ZN10⁶10⁶ZN为过渡电阻常见于弧光接地BC相间短路不接地10⁶ZbcZbcB、C相间直接连接对地开路BC相间接地短路10⁶ZNZNB、C经阻抗接地A相开路这张表的核心是“规范”两个字。无论故障发生在哪一相、是否带过渡电阻都用同一套三相阻抗装配流程程序不需要为每一种故障单独写分支。实际建模时先按表格把每相的阻抗值确定下来再换算成导纳填进矩阵顺手把故障类型名也留在注释里。A相经阻抗接地那一行的ZN如果取0就是金属性接地如果取实际过渡电阻弧光接地、树枝碰线这类带电阻故障也能直接模拟不需要改程序结构。3.2 10⁻⁶与10⁶短路与开路的数值化约定课件里专门强调了“采用10⁻⁶代表0、采用10⁶代表∞”这个约定在数值计算里几乎是必须的原因有两个。第一个原因是避免除零。理想短路阻抗是0直接放在分母上做导纳运算导纳就是无穷大矩阵元素直接爆掉理想开路的阻抗是无穷大导纳是0对角元素为0又会把矩阵推成奇异。10⁻⁶和10⁶这两个数量级是工程上通用的替代值短路导纳高达10⁶开路导纳低到10⁻⁶数值上既不溢出、又不奇异计算结果和理想值之间的误差在6个数量级以下对工程精度完全够。第二个原因是让所有故障统一进一个装配流程。不管你是短路还是开路在程序里都是同一个填表函数只是那个位置的阻抗值写10⁻⁶还是10⁶的区别。这个约定的坑在于方向容易记反10⁻⁶是给短路支路的10⁶是给开路支路的。我见过不止一次同事把故障相的阻抗写成10⁶结果A相接地故障算出来三相电流几乎为零误以为保护拒动其实是短路支路被填成了开路。后面第5章排查章节会再展开这里先记住方向。3.3 从阻抗装配到三相导纳阵程序里真正写入的矩阵阻抗表只是人看的程序算的是导纳阵。以一个单节点三相网络为例节点导纳阵按相展开排布规则是对角元素相加、非对角元素取负对角元素Yₐₐ′是A相所有并联导纳之和包括A相对地导纳、A相上的并联电抗器等非对角元素−Yₐᵦ是A、B相之间的互导纳取负号是因为互导纳的电流方向和自导纳相反。以A相金属性接地为例装配时A相的对地阻抗填10⁻⁶导纳是10⁶B、C相对地阻抗填10⁶导纳是10⁻⁶。装配完成之后的Y矩阵代入I YU解出来就是各相的电压和电流。程序里往往是先算I YU再对需要求电压的节点做U Y⁻¹I两步都绕不开这个装配好的矩阵。很多初学的人会在这一步卡住为什么不先画故障网络再手推公式答案很简单手推公式只适用于教材里的三种对称故障场景一旦出现带过渡电阻、非线性电弧或者多故障点并存手推的边界条件就崩了而导纳阵装配完全不受影响你只需要改对应位置的阻抗值。课件里“用规范的阻抗形式模拟不对称故障”这一整节落脚点就是把这个矩阵装配流程讲透后面所有等值阻抗计算都建立在它之上。4. 等值阻抗通用公式与支路内部故障新增节点、消去节点的一步到位4.1 三相导纳阵求逆等值阻抗通用公式的落地课件里等值阻抗通用公式的推导建立在三相导纳阵的基础上。前面装配出的Y阵求逆之后得到的就是等值阻抗阵Z Y⁻¹。对3×3矩阵程序里可以直接调矩阵求逆手算时用伴随矩阵法Z adj(Y) / det(Y)其中det(Y)对于三相系统可以展开成三项乘积的差。具体展开式不一定要背程序也不会让你手算行列式但有两个物理解读值得记住。对角元素Zₐₐ表示A相看进去的等值自阻抗也就是A相注入单位电流时A相的电压响应非对角元素Zₐᵦ表示A、B两相之间的互阻抗它的存在说明三相网络在故障后不再解耦这正是对称分量法在不对称故障时要逐序分解的原因。课件里把这个求逆过程称作等值阻抗通用公式想表达的是任一三相网络只要Y阵装配正确等值阻抗就由Y⁻¹唯一确定不需要针对故障类型手动构造等值电路。程序实现时直接调numpy.linalg.inv或者Matlab的inv()但验证结果时还是要回到伴随矩阵的物理解读上来防止矩阵奇异时程序返回垃圾结果你还不知道。4.2 支路内部故障的规范化处理新增节点把故障位置变成普通节点课件最后一部分处理的是线路内部故障。一段线路i–j中间距离i点不远处f点发生故障常见做法是把f点新增为节点原线路被拆成两段i–f段导纳Y_L1、f–j段导纳Y_L2f点经故障阻抗Z_F连接到地。节点导纳阵从2×2扩成3×3排布如下节点ifjiYii Y_L1−Y_L10f−Y_L1Y_L1 Y_L2 Y_F−Y_L2j0−Y_L2Yjj Y_L2其中Y_F 1/Z_F是故障点对地导纳金属性接地时Z_F取10⁻⁶Y_F就是10⁶。这里Y_L1和Y_L2不是随便取的如果线路i–j总导纳是Y_L故障点距i端比例为p那么Y_L1 Y_L/p、Y_L2 Y_L/(1−p)看着反直觉但物理上是对的线路段越短导纳越大。这个p就是故障位置参数程序里只要改p就能扫出沿线任意位置故障的电气量曲线这在保护定值校核里非常实用。4.3 消去故障节点后的网络合并把故障拍平回原来的节点规模f节点在物理上是真实的故障点但网络求解时它没有注入电流属于无源节点。把f节点从方程组里消去实际上是对f所在行做高斯消去消去之后i、j两点之间出现新的转移导纳原网络节点规模恢复。按3×3矩阵手动消去第二行第二列可以得到i–j转移导纳变为−Y_L1·Y_L2 / (Y_L1 Y_L2 Y_F)i点对地导纳修正Yii Y_L1 − Y_L1² / (Y_L1 Y_L2 Y_F)j点对地导纳修正Yjj Y_L2 − Y_L2² / (Y_L1 Y_L2 Y_F)物理解读很直观故障点引入的导纳Y_F通过分压作用落到i、j两端改变了原线路的等值对地导纳和互导纳。如果Y_F特别大金属性接地i–j转移导纳趋近于−Y_L1·Y_L2/(Y_L1Y_L2)相当于线路被故障点“劈”成了两段并联如果Y_F接近开路修正项趋近于0原线路保持不变。这套做法最大的好处是规范化不管是节点故障还是支路内部任意位置的故障最终都走同一条新增节点→装配导纳阵→消去无源节点→求逆流程故障位置参数化之后程序里只需改p和Y_F不需要重写整个网络方程。5. 阻抗模拟常见问题排查四个翻车点与对应处理5.1 10⁻⁶和10⁶的方向填反故障直接算成开路现象单相接地故障算出来的故障电流几乎为0三相电压全部接近正常保护逻辑怎么校都校不上。原因程序里把故障相的阻抗填成了10⁶金属性短路被模拟成了开路故障电流自然出不来。解决装配时统一先写阻抗表再转导纳阵给短路和开路各定义一个常量比如Z_SHORT 1e-6、Z_OPEN 1e6注释写明“短路用1e-6、开路用1e6”不要在装配函数里直接写裸数字。所有故障类型都从这个常量表取值方向就不会错。5.2 开路相对角导纳填零矩阵奇异求逆失败现象Y阵求逆报singular或者算出来的电压出现inf/nan程序直接崩。原因开路相的对地导纳填了0对角元素为零矩阵不可逆。这是最容易忽略的细节短路相用了10⁻⁶觉得开路相“不接东西”就直接写0结果矩阵奇异。解决开路相对地导纳也必须有值用10⁻⁶量级替代0。数值上是“漏电”的开路物理效果接近理想开路但矩阵保持非奇异。这就是课件里10⁻⁶约定最重要的工程含义不只短路要避零开路也要兜底。5.3 Ward等值的外部负荷用恒阻抗近似运行点一变误差失控现象潮流算例里等值前后的损耗对不上暂态仿真时边界母线电压偏差超过工程容忍度。原因等值时把外部系统负荷当恒阻抗吸收进去但实际负荷是恒功率特性运行点变了恒阻抗模型的注入就不对了等值精度随之劣化。解决把外部系统中对边界影响大的发电机节点保留下来或者用诺顿等值注入替代纯阻抗等值再进一步的做法是给边界节点配上可调的等值注入随运行点修正。别为了省计算量把边界划得太小留两三个关键机组比堆几百个远端节点有用得多。5.4 课件公式符号和教材对不上代码也跟着写错现象照着课件公式写程序解出来的电压符号反了或者矩阵维度对不上debug半天找不到根源。原因课件用A X B的分块形式表述Ward等值符号体系和常用的节点导纳阵Y U I不一样A可能是导纳阵也可能是阻抗阵B的符号约定各套教材不一致直接照抄容易抄错。解决动手写代码前先用最小两节点算例手推一遍把A₁₁、A₁₂、B₁逐个对上再落到代码。我一般会把课件每个公式在笔记本上对照教材重写一遍标出符号含义表再开写这一步能省掉后面大量debug时间尤其是涉及分块矩阵维度时提前标清楚每一块的尺寸比事后猜快得多。6. 进阶验证用一个最小Python脚本把课件公式跑通6.1 最小验证脚本课件看到这个程度最有效的验证方式是把A相接地故障重新实现一遍。我用Python演示最小流程装配三相导纳阵、求逆得等值阻抗、解出故障电流三步对应课件前三节的核心内容。import numpy as np Z_SHORT 1e-6 # 短路阻抗模拟0 Z_OPEN 1e6 # 开路阻抗模拟无穷大 # A相金属性接地A相短路接地B/C相开路 z_a Z_SHORT z_b Z_OPEN z_c Z_OPEN # 三相单节点导纳阵只含对地导纳无相间耦合 y np.array([ [1/z_a, 0, 0], [0, 1/z_b, 0], [0, 0, 1/z_c] ], dtypecomplex) # 等值阻抗阵 导纳阵求逆 z_eq np.linalg.inv(y) # 施加A相单位电压B/C相无注入求解故障电流 u np.array([10j, 0, 0]) i y u print(等值自阻抗 Z_AA , abs(z_eq[0, 0])) print(故障电流 I_A , abs(i[0]), 标幺值1e6量级说明金属性短路)这段代码对应课件里A相接地那一页的完整计算链。参数说明Z_SHORT是金属性短路的数值替身Z_OPEN是开路的数值替身y矩阵的对角元素分别对应A、B、C相对地导纳求逆后z_eq的对角元素就是各相等值自阻抗故障电流直接用导纳阵乘电压得到。对于BC相间短路只需要把z_b和z_c都改成Z_SHORTy矩阵改成含相间耦合的形式其余逻辑完全不用动。6.2 把课件公式转成自己的笔记与汇报材料最后说一个课件整理的技巧。这份PPT里公式密集直接在屏幕上反复翻页效率很低。我习惯把每一页关键公式截图导出成高清图片再插进自己的OneNote或汇报文档里旁边用手写批注写清参数含义。导出时记得把页面分辨率调高别用默认压缩否则公式边缘全是锯齿放大就糊了。整理完再对照第5章那几个坑做一轮自查基本就能把课件的计算链吃透。从那以后我每拿到一套矩阵公式密集的课件都强制走同一个流程先按故障类型填阻抗表再核对短路/开路量级方向最后用最小脚本复算一遍。这套流程走下来再回自己的仿真工程你会明显感觉到故障计算不是背公式而是搭矩阵、填参数、求逆三步走。希望帮到你。本文还有配套的精品资源点击获取