电力系统故障运行潮流分析
搞电力系统的人都知道，系统故障时的潮流分析比正常工况刺激多了。就像你正吃着火锅唱着歌，突然变压器炸了，这时候要快速判断全网哪个节点电压会崩、哪条线路会过载，这时候故障潮流计算就是你的救命稻草。
先看个实际场景：某35节点系统发生三相短路，我们需要快速计算故障后的电压分布。用Python+PYPOWER搞这事挺方便，先加载系统数据：
import pypower.api as pp from pypower.runpf import runpf case_data = pp.case35() # 正常状态潮流计算 result = runpf(case_data)[0]这段代码跑的是常规潮流计算，但故障分析得加点料。比如在节点15接入短路阻抗，这时候的处理关键在修改导纳矩阵：
def add_fault(case, bus_idx, Zf): # 获取原始导纳矩阵 Ybus, _, _ = pp.makeYbus(case) # 在故障节点追加接地导纳 Ybus[bus_idx, bus_idx] += 1/Zf return Ybus # 节点15发生金属性接地短路（Zf=0.0001j） fault_Ybus = add_fault(case_data, 14, 0.0001j) # 注意Python索引从0开始这里有个坑：直接修改导纳矩阵后，传统牛顿-拉夫逊法可能会发散。实战中更常用的是叠加法——先算正常潮流，再叠加上故障分量。看这段故障电流计算：
V_normal = result['bus'][:, 7] + 1j*result['bus'][:, 8] # 正常电压 Zbus = np.linalg.inv(fault_Ybus) # 节点阻抗矩阵 If = V_normal[14] / (Zbus[14,14] + 0.0001j) # 故障点电流拿到故障电流后，重点来了——全网电压跌了多少？用这个公式暴力计算：
delta_V = Zbus[:,14] * If V_fault = V_normal - delta_V这时候你可能会骂娘："节点22电压掉到0.7pu了！" 赶紧用matplotlib画个电压分布热力图：
plt.figure(figsize=(10,6)) plt.scatter(case_data['bus'][:,1], case_data['bus'][:,2], c=np.abs(V_fault), cmap='RdYlGn') plt.colorbar(label='Voltage (p.u.)') plt.title('Post-fault Voltage Distribution') plt.xlabel('X Coordinate') plt.ylabel('Y Coordinate')!电压分布热力图
看到没？故障点附近一片飘红，这时候要赶紧拉闸保命。但实际系统可比这复杂得多，比如考虑发电机暂态电抗、负荷动态特性啥的。不过这个代码框架已经能帮你抓住主要矛盾了。
最后说个冷知识：故障计算时别迷信PQ节点模型。系统崩溃前两秒，那些恒功率负荷其实早变成恒阻抗特性了，这时候用ZIP负荷模型更靠谱。下次有机会再唠这个坑。
