📚 课程目录第4章谐波潮流计算基础

4. 谐波潮流计算基础

谐波源定位不能只靠“看波形”“比幅值”，它需要一个能描述谐波在电网中如何流动、在哪里被吸收、在哪里被放大的计算框架。谐波潮流计算正是这个框架：把工频网络在某一谐波次数下重新建模，求出各节点谐波电压和各支路谐波电流，再据此判断谐波源的位置与责任。本节从建模、导纳矩阵、求解方法和定位用途四个层面展开。

一、谐波网络建模：把工频元件“翻译”成谐波域元件

谐波潮流计算的基本假设是：在 h 次谐波下，系统仍可视为线性网络，但各元件的参数随频率变化。因此建模的核心是给出每个元件在 h 次谐波下的等值阻抗或导纳。

同步发电机：用 h 次谐波下的负序电抗近似，Z_G,h ≈ jhX_d″，X_d″ 为次暂态电抗。发电机通常作为谐波电压源或背景谐波源处理。
变压器：忽略励磁支路时，Z_T,h ≈ R_T + jhX_T；若考虑绕组间电容和铁芯饱和，需增加并联支路，但定位计算中常先取简化模型。
输电线路：用分布参数或 π 型等值，Z_L,h = R + jhX_L，并联电容导纳 Y_C,h = jhωC。长线路在 h 较大时电容效应明显，会改变谐波流向。
负荷：这是建模中最不确定的部分。常用恒阻抗模型 Z_load,h = R_load + jhX_load，也可用恒功率、恒电流或综合模型。负荷的谐波阻抗直接影响谐波在母线上的分配。
谐波源：整流器、变频器、电弧炉等用 h 次谐波电流源 I_h 表示，注入节点。其幅值和相位是定位计算的关键输入。
并联电容器组：Z_C,h = 1/(jhωC)，在某个谐波次数附近可能与系统电感发生并联谐振，导致谐波电流放大。

建模时先确定参考节点（通常取大地或某一平衡节点），规定所有电流从节点流向参考方向为正，全文及图示保持该约定。节点编号、支路编号和导纳矩阵的行列顺序必须一致，否则后续求解会错位。

二、节点导纳矩阵：谐波网络的“骨架”

对 h 次谐波，把每个元件的等值导纳填入节点导纳矩阵 Y_h。矩阵元素规则为：

对角元素 Y_ii,h = 与节点 i 相连的所有支路导纳之和，加上节点 i 对参考节点的并联导纳；
非对角元素 Y_ij,h = −节点 i 与节点 j 之间支路导纳之和（互导纳取负）。

以三节点系统为例，节点 1 为谐波源注入节点，节点 2 为母线负荷，节点 3 为并联电容器。设支路 1-2 导纳为 y_12,h，支路 2-3 导纳为 y_23,h，节点 3 对地电容导纳为 y_C,h，则：

Y_h = [ y12,h        -y12,h         0      ]
      [ -y12,h   y12,h+y23,h     -y23,h   ]
      [ 0          -y23,h     y23,h+yC,h ]

这个矩阵是稀疏的，实际系统中可用稀疏存储。它的物理意义是：节点电压向量 U_h 与注入电流向量 I_h 满足 I_h = Y_h · U_h。谐波源定位中，导纳矩阵决定了谐波电流在各节点间的分配比例，是判断“谁注入、谁吸收”的基础。

三、谐波潮流求解方法

谐波潮流本质上是求解线性方程组 Y_hU_h = I_h。常用方法有三类：

直接法：对每个谐波次数 h，直接对 Y_h 求逆或做 LU 分解，得到 U_h = Y_h⁻¹I_h。适用于节点数不大、谐波次数不多的场景，结果精确。
迭代法：当负荷模型含非线性或需要与基波潮流交替求解时，用高斯-赛德尔或牛顿类迭代。每次迭代更新节点注入电流，直到电压变化小于收敛阈值，例如 max|ΔU_h| < 1×10⁻⁵ p.u.。
解耦法：利用各次谐波之间近似独立的特点，对每个 h 分别建立 Y_h 并求解，不跨谐波次数耦合。这是工程中最常用的方式，计算量小，便于定位分析。

求解后得到每个节点的 h 次谐波电压 U_i,h 和每条支路的 h 次谐波电流 I_ij,h = (U_i,h − U_j,h)·y_ij,h。这些量就是定位判据的原始数据。

四、在谐波源定位中的用途

谐波潮流计算结果可直接用于三类定位判断：

谐波流向判断：若某支路 h 次谐波有功功率方向从节点 i 流向节点 j，说明谐波源更可能在 i 侧。定义 P_h = Re(U_i,h · I_ij,h^*)，P_h > 0 表示 i 侧注入，P_h < 0 表示 j 侧注入。
谐波责任划分：在公共连接点 PCC 处，用系统侧等值谐波阻抗和用户侧等值谐波阻抗，把 PCC 谐波电压分解为系统背景贡献和用户注入贡献，比值即为责任分摊依据。
谐振与放大识别：通过扫描不同 h 下的 Y_h 条件数或节点电压放大倍数，找出并联谐振点。若某节点在 h=5 附近电压放大倍数超过 3，说明该处电容器与系统电感可能谐振，定位时需优先排查该节点附近的谐波源。

定位计算中，谐波源电流的相位往往比幅值更难获取。工程上可先用幅值排序锁定可疑节点，再用相位关系或有功方向做二次确认，避免单凭幅值误判。

把上述流程落到一个具体算例：三节点系统，基波下支路 1-2 电抗 X₁₂ = 0.1 p.u.，支路 2-3 电抗 X₂₃ = 0.15 p.u.，节点 3 并联电容器 C 在 h=5 时容抗 X_C,5 = 0.08 p.u.。取 h=5，则 y₁₂,₅ = 1/(j5×0.1) = −j2.0 p.u.，y₂₃,₅ = 1/(j5×0.15) = −j1.333 p.u.，y_C,5 = j1/0.08 = j12.5 p.u.。节点 1 注入 5 次谐波电流 I₁,₅ = 1.0∠0° p.u.，节点 2、3 无注入。代入 Y₅U₅ = I₅ 求解，可得到 U₁,₅、U₂,₅、U₃,₅，再算支路电流和 P₅。若 P₅ 在支路 1-2 上为正，则谐波源在节点 1 侧；若节点 3 电压幅值明显高于节点 2，则需检查电容器是否引起 5 次谐波放大。这些数值属于教学示例假设，实际工程应以实测参数和现场数据为准。

谐波潮流计算的结果对负荷模型和电容器参数非常敏感。定位前应至少核对：变压器接线组别、电容器实际投切组数、线路长度与型号、谐波源电流频谱。参数偏差可能导致谐波流向判断反转。

总结来说，谐波网络建模提供元件级参数，导纳矩阵提供网络级结构，谐波潮流求解提供节点和支路级结果，三者共同支撑谐波源定位与责任划分。掌握这条链路，才能在面对实测谐波数据时，从“看到异常”走到“找到源头、分清责任”。