简介面向电力系统工程师与研究人员这份压缩包提供了一套基于 MATLAB 的谐波潮流计算工具重点解决非线性负荷增多背景下电网谐波分析与治理问题。程序通过谐波解耦算法将复杂多变量系统分解为独立子问题降低计算规模并提升效率适用于教学验证与工程估算场景。包内共 9 个文件其中 1 个 m 脚本为主程序8 个 txt 数据文件分别存放不同规模节点系统的参数如节点信息、线路与负荷数据便于直接代入测试。整个压缩包仅 15KB轻量易用适合快速上手和二次开发。目前已有 586 人浏览学习。借助该资源读者可以掌握谐波潮流计算的主流程理解矩阵对角化解耦思路在电力系统中的应用并基于自带数据自行扩展更多节点案例为谐波治理和电能质量分析提供参考。1. 谐波潮流计算拿到 harmonic.rar 先别急着跑搞懂解耦才算入门做电力谐波分析的人手里大概率攒过不少以harmonic命名的压缩包。harmonic.rar这类资源里通常装的是谐波潮流计算的源码或算例核心就两件事一是把电网里的非线性负荷产生的谐波电流算清楚二是把这些电流注入网络后各节点电压的畸变情况算出来。谐波潮流计算不是基波潮流的简单叠加它和常规牛拉法最大的区别在于基波只有一个频率谐波有 3、5、7、9 次一堆频率而且不同频率在网络里的阻抗特性完全不同。谐波解耦则是这一堆计算里最关键的数学前提——没有它每次谐波都要跟基波联立求解规模直接爆炸。这篇笔记就围绕谐波潮流计算与谐波解耦展开讲清楚原理、步骤、参数和坑。2. 为什么要解耦谐波潮流计算的数学模型与求解思路2.1 谐波源建模恒流源注入是工程默认诺顿等效是进阶选项做谐波潮流计算第一步不是写代码而是想清楚谐波源怎么表示。电力系统里典型的谐波源是整流器、变频器、电弧炉这一类非线性设备。工程上最常见的处理方式是恒流源模型——把非线性负荷当成一个向系统注入特定次数谐波电流的电流源注入幅值和相角由设备的典型频谱给定。比如 6 脉动整流器的特征谐波是 6k±1 次即 5、7、11、13 次幅值大致按 1/h 衰减这个衰减规律可以直接写进参数表。恒流源的好处是简单、稳定、收敛快坏处是忽略了非线性负荷端电压变化时注入电流也会变化的物理事实。谐波电压畸变严重时恒流源模型会高估谐波水平。如果要做更精细的评估就把谐波源改成诺顿等效电路一个谐波电流源并联一个谐波阻抗。这个谐波阻抗取负荷在基波下的等值阻抗乘以次数 h即 Z_h h·Z_1。这个近似在 50Hz 到 2500Hz 范围内基本够用但对容性负荷和高次谐振场景会失真。我一般先跑恒流源如果某些母线谐波电压异常高再换诺顿模型复核这样能快速定位到底是源的问题还是网络的问题。2.2 谐波解耦的核心逻辑为什么各次谐波可以分开独立求解谐波解耦这四个字字面上指的是把基波潮流和各次谐波潮流之间的耦合解开各自独立求解。为什么能解耦因为电力系统的稳态响应是线性的——只要系统里的元件参数不随频率混叠不同频率的电压电流在傅里叶变换下是正交的不会互相产生能量交换。基波潮流决定各节点的基波电压幅值和相角这个结果用来修正各次谐波导纳阵里的负荷阻抗各次谐波注入电流只在各自频率的网络阻抗上产生压降和别的频率没有直接关系。用数学语言表达常规潮流求解的是 F(x)0 的非线性方程组谐波潮流加进来后理论上应该是 F(x, I_h)0 的增广方程组。但工程实现很少这么干常见做法是把求解拆成两阶段第一阶段跑一次基波潮流牛拉法或 PQ 分解法得到各节点电压 U_1 和各支路功率第二阶段对每一个关注的谐波次数 h单独构建 h 次谐波导纳阵 Y_h求解 Y_h·U_h I_h。这个两阶段法的名字就叫谐波解耦潮流计算也是 harmonic.rar 这类代码包最常见的内部结构。2.3 谐波导纳阵的构建频率相关参数是最大变量解耦之后工作量全落在构建各次谐波导纳阵上。电力元件在谐波下的阻抗特性和基波完全不一样这是谐波潮流计算里最容易被忽略也最容易出错的地方。输电线路的谐波阻抗按频率线性变化感抗 X_Lh h·X_L1容抗 X_Ch X_C1/h这个线性关系对所有集中参数模型都成立。变压器漏抗同样线性放大但铁损和涡流损耗会随频率上升在计算中常等效为并联电阻阻值大致和 h² 成正比。发电机在谐波下用次暂态电抗 X_d 的 h 倍近似。负荷最麻烦——恒功率负荷在基波潮流里是 PjQ 的功率约束但到了谐波潮流里没有“功率”这个概念只能近似成恒阻抗Z_load_h h·(U_1²/S_load)*。这个近似在电动机占比较高的工业负荷里偏差比较大电动机的等效阻抗频率特性更接近感抗加电阻而不是单纯的线性放大。所以做谐波潮流计算前先看清楚负荷模型是恒阻抗、恒电流还是恒功率这个选择直接决定了结果的可用性。import numpy as np # 构建 5 次谐波导纳阵示例3 节点系统 # Y_bus_h 的每个元素是节点导纳按频率修正后组装 # 支路参数电阻 R, 电抗 X1(基波), 电纳 B1(基波) # h 5 R np.array([0.01, 0.015, 0.008]) # 支路电阻 pu X1 np.array([0.08, 0.10, 0.06]) # 基波电抗 pu B1 np.array([0.02, 0.025, 0.015]) # 基波充电电纳 pu h 5 # 按频率修正感抗放大 h 倍容纳放大 h 倍 Xh h * X1 Bh h * B1 Zh R 1j * Xh # 谐波阻抗 # 导纳 Yh 1 / Zh Yshunt_h 1j * Bh / 2 # 线路对地充电导纳拆两半 # 3 节点系统的节点导纳阵组装简单拓扑1-2, 2-3, 1-3 Ybus_h np.zeros((3, 3), dtypecomplex) for k, (i, j) in enumerate([(0, 1), (1, 2), (0, 2)]): Ybus_h[i, i] Yh[k] Yshunt_h[k] Ybus_h[j, j] Yh[k] Yshunt_h[k] Ybus_h[i, j] - Yh[k] Ybus_h[j, i] - Yh[k] # 注入 5 次谐波电流pu幅值按典型整流器频谱给定 I_h np.array([0, 0.03 - 0.01j, 0], dtypecomplex) # 求解节点谐波电压 U_h np.linalg.solve(Ybus_h, I_h) print(5次谐波节点电压(pu):, np.abs(U_h))这段代码展示了谐波解耦计算最核心的一步把基波导纳阵按频率修正后重新组装然后每个谐波次数独立求解一次线性方程。之所以能这样做是因为前面的谐波解耦已经切断了各次谐波之间的耦合。Xh h * X1和Bh h * B1这两行是整个计算里最关键的两处频率修正漏掉任何一处高次谐波计算结果都会失真。组装导纳阵的过程和基波潮流里的 Y_bus 组装完全一致区别仅在于支路参数按 h 次频率取值。实际工程中很少只用 3 节点算例真实系统的 Y_h 可能上千阶但结构完全一样。稀疏矩阵存储和 LU 分解在 Python 里用 scipy.sparse 就能解决计算耗时几乎可以忽略。谐波潮流计算真正耗时的其实不是矩阵求解而是数据准备——把系统的每个元件参数按谐波频率换算好这个环节出错率最高。3. 用 Harmonic 算例包跑通最小系统目录结构、数据文件与执行流程3.1 算例包目录结构与数据文件定位拿到 harmonic.rar 解压后先别急着双击运行脚本。我习惯先看目录搞清每个文件的角色。典型的谐波潮流代码包会包含这几类内容系统数据文件母线、支路、变压器、负荷、谐波源数据文件各节点的注入频谱、主计算脚本、结果输出脚本。母线数据文件是文本格式常见后缀是.csv或.txt里面每行代表一条母线列依次是母线编号、电压等级、基波电压幅值、相角、有功负荷、无功负荷。支路文件则列出发端母线、受端母线、电阻、电抗、充电电纳、变比。谐波源文件一般按母线编号给出各次谐波注入电流的幅值和相角。数据文件是最容易翻车的地方。很多算例包自带的系统数据来自 IEEE 测试系统比如 IEEE 14 节点或 30 节点这些数据默认是标幺值pu基准容量通常在文件头注释里写着 100 MVA。如果直接把这些数据喂给程序却不统一基准容量谐波计算结果会偏离实际好几个数量级。所以拿到数据文件的第一件事是确认基准容量并且看清楚每个母线的电压等级后续的 THD 百分比计算依赖额定电压作为基准。3.2 算例参数设置注入频谱和负荷模型是决定结果的两个开关参数设置是谐波潮流计算里最能拉开新手和老手差距的环节。注入频谱参数决定谐波源强度通常按设备类型查标准谱6 脉动整流器注入 5、7、11、13 次幅值随次数衰减12 脉动整流器注入 11、13 次为主5、7 次大幅削弱。判断注入参数是否合理的经验法则是5 次谐波电流幅值占基波电流的比例通常在 15% 到 25% 之间超过 30% 说明谐波源模型取得偏保守结果会有较大裕量。IEEE 519 标准里对注入限值有明确建议照着取不会出大错。负荷模型选项直接决定导纳阵里对角元素的取值。前面提过负荷在谐波下的阻抗近似为 h 倍基波阻抗但工程代码里通常提供一个开关LOAD_MODEL Z代表恒阻抗LOAD_MODEL PQ代表把功率换算成阻抗。这两种方式对低次谐波影响不大但到了 23 次以上恒阻抗模型给出的谐波电压往往偏高因为纯线性放大导致高次感抗过大、容抗过小负荷吸收谐波的能力被低估了。算例包里如果默认是恒功率模型建议手动改成恒阻抗试一次对比两次结果的差异再决定最终采用的模型。# 谐波潮流主流程伪代码结构以 harmonic 类工具包常见流程为参考 import numpy as np from scipy.sparse import lil_matrix from scipy.sparse.linalg import spsolve # 1. 读入系统数据 # buses: array of [bus_id, base_kv, P_load, Q_load] # branches: array of [from_bus, to_bus, R, X, B, tap] # harmonic_source: dict {bus_id: {h: [mag, angle]}} # 2. 基波潮流求解牛顿-拉夫逊此处调用外部函数 # U1, theta1 run_power_flow(buses, branches) # 得到基波电压用于后续负荷阻抗修正 # 3. 对关注的谐波次数依次循环 for h in [5, 7, 11, 13, 23, 25]: # 3.1 按频率修正支路阻抗 # Z_h R 1j * h * X # B_h h * B # 3.2 组装谐波导纳阵 Y_h Y_h lil_matrix((n_bus, n_bus), dtypecomplex) # 3.3 构建谐波电流注入向量 I_h np.zeros(n_bus, dtypecomplex) for bus, spec in harmonic_source.items(): if h in spec: mag, ang spec[h] I_h[bus] mag * np.exp(1j * np.deg2rad(ang)) # 3.4 求解 Y_h * U_h I_h U_h spsolve(Y_h.tocsr(), I_h) # 3.5 保存该次谐波电压折算成有名值 U_h_volts U_h * base_kv * 1000 / np.sqrt(3)这段主流程把谐波潮流计算的骨架展示得很清楚。注意 3.3 和 3.1 的顺序——谐波电流注入的幅值和相角是给定参数不随基波潮流结果变化而导纳阵里的负荷阻抗要用基波电压幅值修正所以基波潮流必须先算完。谐波解耦在此处的体现就是基波潮流和每次谐波各算各的循环之间只共享系统参数不共享迭代过程。参数说明里有几个容易被忽略的点。base_kv是母线额定线电压把标幺值谐波电压折算成伏特时用的是相电压形式。如果算例包里出现U_h_volts U_h * base_kv * 1000而不除根号三说明作者用的可能是相电压基准这时候把结果和标准对比会偏差 1.732 倍。还有I_h向量只在有谐波源的母线处非零其余母线全是零很多谐波潮流结果里出现了一个没有谐波源的母线谐波电压畸变率反而比源母线还高正是因为并联谐振把它放大成了局部高电压。3.3 结果文件解读THD、各次谐波电压与向量组算完谐波潮流之后结果文件里最重要的一个指标是电压总畸变率 THD。计算公式是 THD sqrt(Σ_{h≥2} U_h²) / U_1 × 100%其中 U_1 是基波电压有效值U_h 是各次谐波电压有效值。IEEE 519 和国标 GB/T 14549 对公用连接点的 THD 限值有明确规定比如 110kV 及以上电压等级 THD 限值是 2%380V 等级是 5%。谐波潮流计算的最终输出就是要拿 THD 和标准对比判断哪些母线超标、超在哪次。除了 THD结果文件里通常还有各次谐波电压的向量组数据。向量组的价值在于判断谐波传播方向和谐振位置。举例如果 11 次谐波电压在某些母线幅值异常放大而其他母线正常把各母线 11 次谐波电压的幅值画成柱状图峰值母线往往就是并联谐振中心。4. 谐波潮流计算避坑指南5 个高频翻车现场与排查办法4.1 基波潮流不收敛谐波结果全部作废现象基波潮流迭代多次不收敛弹出奇异矩阵告警但谐波计算照常输出了一组数据。原因谐波解耦算法要求基波潮流先收敛因为谐波导纳阵里的负荷阻抗修正依赖基波电压幅值。基波不收敛时程序通常用上一次迭代的电压值硬算结果就是谐波电压全部失真。还有一个隐蔽诱因是数据文件里负荷功率和发电机出力不平衡总有功差值超过了设定的收敛裕度。解决先单独跑基波潮流并观察收敛曲线。如果确实不收敛检查注入总功率是否为零、母线电压初始值是否接近 1.0pu、无功补偿是否配置合理。基波收敛是一切谐波计算的地基这一步化验血泪经验直接跳过基波不收敛去调谐波参数等于开车不扶方向盘只看后视镜。4.2 谐波解耦的前提被三相不平衡打破现象三相系统算出的各次谐波电压严重偏离实测值尤其是 3 次谐波。原因谐波解耦的理论前提是系统三相对称各次谐波正序、负序、零序分量相互独立。但实际电网存在三相不平衡负序分量在基波潮流阶段就会引入不对称导致 3 次谐波不再被 Yyn0 接线的变压器隔绝。当负序不平衡度超过 2% 时解耦误差会明显放大。解决判断系统是否适合解耦算法。对三相不平衡度较高的配电网建议改用三相谐波潮流模型或者在数据准备阶段对不平衡母线做等效归算把负序分量折算进谐波源注入里而不是直接忽略。4.3 变压器接线对零序谐波的隔离被忽视现象计算结果显示某条母线的 3 次谐波电压异常高但实际测量并没有这么严重。原因3、9、15 次等零序谐波在 Yyn0 接线的变压器星形中性点接地-星形中性点接地高压侧和低压侧之间没有通路零序谐波电流被挡在低压侧。很多谐波潮流代码包的变压器模型默认用漏抗加理想变压器表示没有区分接线组别对零序分量的阻断作用。解决检查变压器模型是否包含接线组别信息。计算零序谐波分量时Yyn0 接线的变压器要在铁芯模型里增加励磁支路或直接把零序谐波设为不穿越。Dyn11 接线的变压器给 3 次谐波提供了流通路径计算模型里要把三角形绕组产生的环流效应考虑进去。这条是谐波潮流计算里最容易被忽略的边界条件没有之一。4.4 并联谐振点判断失误只算到 13 次就停了现象某母线电压 THD 虽然没超标但 23 次谐波电压单项畸变率已经逼近限值。原因计算只覆盖到 13 次或 25 次就收工了。并联谐振通常发生在系统容性元件和感性元件阻抗相等的频率点。并联电容器组和变压器漏抗构成经典谐振回路谐振频率一般在 10 到 30 次之间。如果只算到 13 次可能会错过发生在 17 或 23 次的谐振峰。解决先把关注频段扩展到 50 次2500Hz做一次频率扫描找出谐振峰所在次数再决定谐波潮流的最高计算次数。频率扫描的做法是对每个目标次数 h 计算母线看进去的等值阻抗幅值 |Z_h|画出 |Z_h|-h 曲线峰值处就是并联谐振点。只要在谐振频率附近加密计算密度就能避免 THD 评估时漏掉高次分量。4.5 负荷阻抗频率模型的过度简化现象同一个系统用恒阻抗负荷模型和用恒功率换算模型算出的 THD 差异超过 30%。原因谐波潮流里负荷表达没有统一约定恒阻抗模型的感抗随频率线性放大在高次谐波下负荷几乎不吸收谐波电流谐波电压自然偏高恒功率换算模型则把负荷在基波电压下算出的阻抗直接乘以 h两者在高次频段的差异被放大。解决对比两个模型的差异取其平均值或按保守原则取偏大值。判断哪个更合理的标准是负荷构成——电动机为主的工业负荷更接近恒阻抗特性商业和居民负荷更接近恒功率特性。5. 谐波潮流算完不算完用叠加反演和扫参验证结果可信度谐波潮流计算输出一堆数据但第一步要做的不是直接拿去写报告而是把计算结果反演回时域波形做一次合理性验证。做法很简单把基波电压和各次谐波电压按相角叠加还原出母线电压的时域波形。如果叠加后的波形顶底不对称、明显畸变说明各次谐波相角之间的一致性有问题。谐波潮流计算里相角是实数输出各次谐波相角混合后应当还原出可辨识的整流纹波特征。如果反演波形看起来杂乱无章大概率是某个谐波源的相角基准没对齐或者正序负序分量被错误合并。这个土办法在新手阶段帮我抓出过很多次谐波源数据文件里的相角错误。第二个验证手段是对注入电流频谱做灵敏度扫参。电网里非线性负荷实测频谱很难完全准确厂家给的典型值往往和实际有偏差。固定网络参数不变把某条母线的 5 次谐波电流幅值从 50% 扫到 150%观察所有母线 THD 的变化曲线。如果某条母线 THD 对注入幅值极度敏感说明该母线正好落在谐振点附近需要进一步做频率扫描确认谐振频率。如果所有母线 THD 对注入变化都不敏感说明系统对该次谐波有天然抑制能力谐波源模型精度要求可以放宽。这个扫参过程能帮你判断结果的可信边界——一个对参数极度敏感的系统数据文件里任何一个谐波源参数不准都会导致结论失真。我个人的习惯是每跑完一个算例就在结果文件里附加一栏“最大灵敏母线”记录扫参过程中 THD 变化率最大的母线编号。下次做治理方案时直接看这栏就知道电容补偿支路或 APF 应该优先挂在哪条母线。谐波潮流计算本质上是工程判断工具结果的意义不在于数字精确到小数点后三位而在于帮你定位谐振风险最高的位置和频次。用三个手段——时域反演、灵敏度扫参、频率扫描——交叉验证之后计算结论才具备落地的说服力。这套验证流程陪我处理过不少谐波治理项目希望帮到你。本文还有配套的精品资源点击获取 谐波潮流计算与谐波解耦：从原理到工程实践的完整指南
简介面向电力系统工程师与研究人员这份压缩包提供了一套基于 MATLAB 的谐波潮流计算工具重点解决非线性负荷增多背景下电网谐波分析与治理问题。程序通过谐波解耦算法将复杂多变量系统分解为独立子问题降低计算规模并提升效率适用于教学验证与工程估算场景。包内共 9 个文件其中 1 个 m 脚本为主程序8 个 txt 数据文件分别存放不同规模节点系统的参数如节点信息、线路与负荷数据便于直接代入测试。整个压缩包仅 15KB轻量易用适合快速上手和二次开发。目前已有 586 人浏览学习。借助该资源读者可以掌握谐波潮流计算的主流程理解矩阵对角化解耦思路在电力系统中的应用并基于自带数据自行扩展更多节点案例为谐波治理和电能质量分析提供参考。1. 谐波潮流计算拿到 harmonic.rar 先别急着跑搞懂解耦才算入门做电力谐波分析的人手里大概率攒过不少以harmonic命名的压缩包。harmonic.rar这类资源里通常装的是谐波潮流计算的源码或算例核心就两件事一是把电网里的非线性负荷产生的谐波电流算清楚二是把这些电流注入网络后各节点电压的畸变情况算出来。谐波潮流计算不是基波潮流的简单叠加它和常规牛拉法最大的区别在于基波只有一个频率谐波有 3、5、7、9 次一堆频率而且不同频率在网络里的阻抗特性完全不同。谐波解耦则是这一堆计算里最关键的数学前提——没有它每次谐波都要跟基波联立求解规模直接爆炸。这篇笔记就围绕谐波潮流计算与谐波解耦展开讲清楚原理、步骤、参数和坑。2. 为什么要解耦谐波潮流计算的数学模型与求解思路2.1 谐波源建模恒流源注入是工程默认诺顿等效是进阶选项做谐波潮流计算第一步不是写代码而是想清楚谐波源怎么表示。电力系统里典型的谐波源是整流器、变频器、电弧炉这一类非线性设备。工程上最常见的处理方式是恒流源模型——把非线性负荷当成一个向系统注入特定次数谐波电流的电流源注入幅值和相角由设备的典型频谱给定。比如 6 脉动整流器的特征谐波是 6k±1 次即 5、7、11、13 次幅值大致按 1/h 衰减这个衰减规律可以直接写进参数表。恒流源的好处是简单、稳定、收敛快坏处是忽略了非线性负荷端电压变化时注入电流也会变化的物理事实。谐波电压畸变严重时恒流源模型会高估谐波水平。如果要做更精细的评估就把谐波源改成诺顿等效电路一个谐波电流源并联一个谐波阻抗。这个谐波阻抗取负荷在基波下的等值阻抗乘以次数 h即 Z_h h·Z_1。这个近似在 50Hz 到 2500Hz 范围内基本够用但对容性负荷和高次谐振场景会失真。我一般先跑恒流源如果某些母线谐波电压异常高再换诺顿模型复核这样能快速定位到底是源的问题还是网络的问题。2.2 谐波解耦的核心逻辑为什么各次谐波可以分开独立求解谐波解耦这四个字字面上指的是把基波潮流和各次谐波潮流之间的耦合解开各自独立求解。为什么能解耦因为电力系统的稳态响应是线性的——只要系统里的元件参数不随频率混叠不同频率的电压电流在傅里叶变换下是正交的不会互相产生能量交换。基波潮流决定各节点的基波电压幅值和相角这个结果用来修正各次谐波导纳阵里的负荷阻抗各次谐波注入电流只在各自频率的网络阻抗上产生压降和别的频率没有直接关系。用数学语言表达常规潮流求解的是 F(x)0 的非线性方程组谐波潮流加进来后理论上应该是 F(x, I_h)0 的增广方程组。但工程实现很少这么干常见做法是把求解拆成两阶段第一阶段跑一次基波潮流牛拉法或 PQ 分解法得到各节点电压 U_1 和各支路功率第二阶段对每一个关注的谐波次数 h单独构建 h 次谐波导纳阵 Y_h求解 Y_h·U_h I_h。这个两阶段法的名字就叫谐波解耦潮流计算也是 harmonic.rar 这类代码包最常见的内部结构。2.3 谐波导纳阵的构建频率相关参数是最大变量解耦之后工作量全落在构建各次谐波导纳阵上。电力元件在谐波下的阻抗特性和基波完全不一样这是谐波潮流计算里最容易被忽略也最容易出错的地方。输电线路的谐波阻抗按频率线性变化感抗 X_Lh h·X_L1容抗 X_Ch X_C1/h这个线性关系对所有集中参数模型都成立。变压器漏抗同样线性放大但铁损和涡流损耗会随频率上升在计算中常等效为并联电阻阻值大致和 h² 成正比。发电机在谐波下用次暂态电抗 X_d 的 h 倍近似。负荷最麻烦——恒功率负荷在基波潮流里是 PjQ 的功率约束但到了谐波潮流里没有“功率”这个概念只能近似成恒阻抗Z_load_h h·(U_1²/S_load)*。这个近似在电动机占比较高的工业负荷里偏差比较大电动机的等效阻抗频率特性更接近感抗加电阻而不是单纯的线性放大。所以做谐波潮流计算前先看清楚负荷模型是恒阻抗、恒电流还是恒功率这个选择直接决定了结果的可用性。import numpy as np # 构建 5 次谐波导纳阵示例3 节点系统 # Y_bus_h 的每个元素是节点导纳按频率修正后组装 # 支路参数电阻 R, 电抗 X1(基波), 电纳 B1(基波) # h 5 R np.array([0.01, 0.015, 0.008]) # 支路电阻 pu X1 np.array([0.08, 0.10, 0.06]) # 基波电抗 pu B1 np.array([0.02, 0.025, 0.015]) # 基波充电电纳 pu h 5 # 按频率修正感抗放大 h 倍容纳放大 h 倍 Xh h * X1 Bh h * B1 Zh R 1j * Xh # 谐波阻抗 # 导纳 Yh 1 / Zh Yshunt_h 1j * Bh / 2 # 线路对地充电导纳拆两半 # 3 节点系统的节点导纳阵组装简单拓扑1-2, 2-3, 1-3 Ybus_h np.zeros((3, 3), dtypecomplex) for k, (i, j) in enumerate([(0, 1), (1, 2), (0, 2)]): Ybus_h[i, i] Yh[k] Yshunt_h[k] Ybus_h[j, j] Yh[k] Yshunt_h[k] Ybus_h[i, j] - Yh[k] Ybus_h[j, i] - Yh[k] # 注入 5 次谐波电流pu幅值按典型整流器频谱给定 I_h np.array([0, 0.03 - 0.01j, 0], dtypecomplex) # 求解节点谐波电压 U_h np.linalg.solve(Ybus_h, I_h) print(5次谐波节点电压(pu):, np.abs(U_h))这段代码展示了谐波解耦计算最核心的一步把基波导纳阵按频率修正后重新组装然后每个谐波次数独立求解一次线性方程。之所以能这样做是因为前面的谐波解耦已经切断了各次谐波之间的耦合。Xh h * X1和Bh h * B1这两行是整个计算里最关键的两处频率修正漏掉任何一处高次谐波计算结果都会失真。组装导纳阵的过程和基波潮流里的 Y_bus 组装完全一致区别仅在于支路参数按 h 次频率取值。实际工程中很少只用 3 节点算例真实系统的 Y_h 可能上千阶但结构完全一样。稀疏矩阵存储和 LU 分解在 Python 里用 scipy.sparse 就能解决计算耗时几乎可以忽略。谐波潮流计算真正耗时的其实不是矩阵求解而是数据准备——把系统的每个元件参数按谐波频率换算好这个环节出错率最高。3. 用 Harmonic 算例包跑通最小系统目录结构、数据文件与执行流程3.1 算例包目录结构与数据文件定位拿到 harmonic.rar 解压后先别急着双击运行脚本。我习惯先看目录搞清每个文件的角色。典型的谐波潮流代码包会包含这几类内容系统数据文件母线、支路、变压器、负荷、谐波源数据文件各节点的注入频谱、主计算脚本、结果输出脚本。母线数据文件是文本格式常见后缀是.csv或.txt里面每行代表一条母线列依次是母线编号、电压等级、基波电压幅值、相角、有功负荷、无功负荷。支路文件则列出发端母线、受端母线、电阻、电抗、充电电纳、变比。谐波源文件一般按母线编号给出各次谐波注入电流的幅值和相角。数据文件是最容易翻车的地方。很多算例包自带的系统数据来自 IEEE 测试系统比如 IEEE 14 节点或 30 节点这些数据默认是标幺值pu基准容量通常在文件头注释里写着 100 MVA。如果直接把这些数据喂给程序却不统一基准容量谐波计算结果会偏离实际好几个数量级。所以拿到数据文件的第一件事是确认基准容量并且看清楚每个母线的电压等级后续的 THD 百分比计算依赖额定电压作为基准。3.2 算例参数设置注入频谱和负荷模型是决定结果的两个开关参数设置是谐波潮流计算里最能拉开新手和老手差距的环节。注入频谱参数决定谐波源强度通常按设备类型查标准谱6 脉动整流器注入 5、7、11、13 次幅值随次数衰减12 脉动整流器注入 11、13 次为主5、7 次大幅削弱。判断注入参数是否合理的经验法则是5 次谐波电流幅值占基波电流的比例通常在 15% 到 25% 之间超过 30% 说明谐波源模型取得偏保守结果会有较大裕量。IEEE 519 标准里对注入限值有明确建议照着取不会出大错。负荷模型选项直接决定导纳阵里对角元素的取值。前面提过负荷在谐波下的阻抗近似为 h 倍基波阻抗但工程代码里通常提供一个开关LOAD_MODEL Z代表恒阻抗LOAD_MODEL PQ代表把功率换算成阻抗。这两种方式对低次谐波影响不大但到了 23 次以上恒阻抗模型给出的谐波电压往往偏高因为纯线性放大导致高次感抗过大、容抗过小负荷吸收谐波的能力被低估了。算例包里如果默认是恒功率模型建议手动改成恒阻抗试一次对比两次结果的差异再决定最终采用的模型。# 谐波潮流主流程伪代码结构以 harmonic 类工具包常见流程为参考 import numpy as np from scipy.sparse import lil_matrix from scipy.sparse.linalg import spsolve # 1. 读入系统数据 # buses: array of [bus_id, base_kv, P_load, Q_load] # branches: array of [from_bus, to_bus, R, X, B, tap] # harmonic_source: dict {bus_id: {h: [mag, angle]}} # 2. 基波潮流求解牛顿-拉夫逊此处调用外部函数 # U1, theta1 run_power_flow(buses, branches) # 得到基波电压用于后续负荷阻抗修正 # 3. 对关注的谐波次数依次循环 for h in [5, 7, 11, 13, 23, 25]: # 3.1 按频率修正支路阻抗 # Z_h R 1j * h * X # B_h h * B # 3.2 组装谐波导纳阵 Y_h Y_h lil_matrix((n_bus, n_bus), dtypecomplex) # 3.3 构建谐波电流注入向量 I_h np.zeros(n_bus, dtypecomplex) for bus, spec in harmonic_source.items(): if h in spec: mag, ang spec[h] I_h[bus] mag * np.exp(1j * np.deg2rad(ang)) # 3.4 求解 Y_h * U_h I_h U_h spsolve(Y_h.tocsr(), I_h) # 3.5 保存该次谐波电压折算成有名值 U_h_volts U_h * base_kv * 1000 / np.sqrt(3)这段主流程把谐波潮流计算的骨架展示得很清楚。注意 3.3 和 3.1 的顺序——谐波电流注入的幅值和相角是给定参数不随基波潮流结果变化而导纳阵里的负荷阻抗要用基波电压幅值修正所以基波潮流必须先算完。谐波解耦在此处的体现就是基波潮流和每次谐波各算各的循环之间只共享系统参数不共享迭代过程。参数说明里有几个容易被忽略的点。base_kv是母线额定线电压把标幺值谐波电压折算成伏特时用的是相电压形式。如果算例包里出现U_h_volts U_h * base_kv * 1000而不除根号三说明作者用的可能是相电压基准这时候把结果和标准对比会偏差 1.732 倍。还有I_h向量只在有谐波源的母线处非零其余母线全是零很多谐波潮流结果里出现了一个没有谐波源的母线谐波电压畸变率反而比源母线还高正是因为并联谐振把它放大成了局部高电压。3.3 结果文件解读THD、各次谐波电压与向量组算完谐波潮流之后结果文件里最重要的一个指标是电压总畸变率 THD。计算公式是 THD sqrt(Σ_{h≥2} U_h²) / U_1 × 100%其中 U_1 是基波电压有效值U_h 是各次谐波电压有效值。IEEE 519 和国标 GB/T 14549 对公用连接点的 THD 限值有明确规定比如 110kV 及以上电压等级 THD 限值是 2%380V 等级是 5%。谐波潮流计算的最终输出就是要拿 THD 和标准对比判断哪些母线超标、超在哪次。除了 THD结果文件里通常还有各次谐波电压的向量组数据。向量组的价值在于判断谐波传播方向和谐振位置。举例如果 11 次谐波电压在某些母线幅值异常放大而其他母线正常把各母线 11 次谐波电压的幅值画成柱状图峰值母线往往就是并联谐振中心。4. 谐波潮流计算避坑指南5 个高频翻车现场与排查办法4.1 基波潮流不收敛谐波结果全部作废现象基波潮流迭代多次不收敛弹出奇异矩阵告警但谐波计算照常输出了一组数据。原因谐波解耦算法要求基波潮流先收敛因为谐波导纳阵里的负荷阻抗修正依赖基波电压幅值。基波不收敛时程序通常用上一次迭代的电压值硬算结果就是谐波电压全部失真。还有一个隐蔽诱因是数据文件里负荷功率和发电机出力不平衡总有功差值超过了设定的收敛裕度。解决先单独跑基波潮流并观察收敛曲线。如果确实不收敛检查注入总功率是否为零、母线电压初始值是否接近 1.0pu、无功补偿是否配置合理。基波收敛是一切谐波计算的地基这一步化验血泪经验直接跳过基波不收敛去调谐波参数等于开车不扶方向盘只看后视镜。4.2 谐波解耦的前提被三相不平衡打破现象三相系统算出的各次谐波电压严重偏离实测值尤其是 3 次谐波。原因谐波解耦的理论前提是系统三相对称各次谐波正序、负序、零序分量相互独立。但实际电网存在三相不平衡负序分量在基波潮流阶段就会引入不对称导致 3 次谐波不再被 Yyn0 接线的变压器隔绝。当负序不平衡度超过 2% 时解耦误差会明显放大。解决判断系统是否适合解耦算法。对三相不平衡度较高的配电网建议改用三相谐波潮流模型或者在数据准备阶段对不平衡母线做等效归算把负序分量折算进谐波源注入里而不是直接忽略。4.3 变压器接线对零序谐波的隔离被忽视现象计算结果显示某条母线的 3 次谐波电压异常高但实际测量并没有这么严重。原因3、9、15 次等零序谐波在 Yyn0 接线的变压器星形中性点接地-星形中性点接地高压侧和低压侧之间没有通路零序谐波电流被挡在低压侧。很多谐波潮流代码包的变压器模型默认用漏抗加理想变压器表示没有区分接线组别对零序分量的阻断作用。解决检查变压器模型是否包含接线组别信息。计算零序谐波分量时Yyn0 接线的变压器要在铁芯模型里增加励磁支路或直接把零序谐波设为不穿越。Dyn11 接线的变压器给 3 次谐波提供了流通路径计算模型里要把三角形绕组产生的环流效应考虑进去。这条是谐波潮流计算里最容易被忽略的边界条件没有之一。4.4 并联谐振点判断失误只算到 13 次就停了现象某母线电压 THD 虽然没超标但 23 次谐波电压单项畸变率已经逼近限值。原因计算只覆盖到 13 次或 25 次就收工了。并联谐振通常发生在系统容性元件和感性元件阻抗相等的频率点。并联电容器组和变压器漏抗构成经典谐振回路谐振频率一般在 10 到 30 次之间。如果只算到 13 次可能会错过发生在 17 或 23 次的谐振峰。解决先把关注频段扩展到 50 次2500Hz做一次频率扫描找出谐振峰所在次数再决定谐波潮流的最高计算次数。频率扫描的做法是对每个目标次数 h 计算母线看进去的等值阻抗幅值 |Z_h|画出 |Z_h|-h 曲线峰值处就是并联谐振点。只要在谐振频率附近加密计算密度就能避免 THD 评估时漏掉高次分量。4.5 负荷阻抗频率模型的过度简化现象同一个系统用恒阻抗负荷模型和用恒功率换算模型算出的 THD 差异超过 30%。原因谐波潮流里负荷表达没有统一约定恒阻抗模型的感抗随频率线性放大在高次谐波下负荷几乎不吸收谐波电流谐波电压自然偏高恒功率换算模型则把负荷在基波电压下算出的阻抗直接乘以 h两者在高次频段的差异被放大。解决对比两个模型的差异取其平均值或按保守原则取偏大值。判断哪个更合理的标准是负荷构成——电动机为主的工业负荷更接近恒阻抗特性商业和居民负荷更接近恒功率特性。5. 谐波潮流算完不算完用叠加反演和扫参验证结果可信度谐波潮流计算输出一堆数据但第一步要做的不是直接拿去写报告而是把计算结果反演回时域波形做一次合理性验证。做法很简单把基波电压和各次谐波电压按相角叠加还原出母线电压的时域波形。如果叠加后的波形顶底不对称、明显畸变说明各次谐波相角之间的一致性有问题。谐波潮流计算里相角是实数输出各次谐波相角混合后应当还原出可辨识的整流纹波特征。如果反演波形看起来杂乱无章大概率是某个谐波源的相角基准没对齐或者正序负序分量被错误合并。这个土办法在新手阶段帮我抓出过很多次谐波源数据文件里的相角错误。第二个验证手段是对注入电流频谱做灵敏度扫参。电网里非线性负荷实测频谱很难完全准确厂家给的典型值往往和实际有偏差。固定网络参数不变把某条母线的 5 次谐波电流幅值从 50% 扫到 150%观察所有母线 THD 的变化曲线。如果某条母线 THD 对注入幅值极度敏感说明该母线正好落在谐振点附近需要进一步做频率扫描确认谐振频率。如果所有母线 THD 对注入变化都不敏感说明系统对该次谐波有天然抑制能力谐波源模型精度要求可以放宽。这个扫参过程能帮你判断结果的可信边界——一个对参数极度敏感的系统数据文件里任何一个谐波源参数不准都会导致结论失真。我个人的习惯是每跑完一个算例就在结果文件里附加一栏“最大灵敏母线”记录扫参过程中 THD 变化率最大的母线编号。下次做治理方案时直接看这栏就知道电容补偿支路或 APF 应该优先挂在哪条母线。谐波潮流计算本质上是工程判断工具结果的意义不在于数字精确到小数点后三位而在于帮你定位谐振风险最高的位置和频次。用三个手段——时域反演、灵敏度扫参、频率扫描——交叉验证之后计算结论才具备落地的说服力。这套验证流程陪我处理过不少谐波治理项目希望帮到你。本文还有配套的精品资源点击获取