搞配电网仿真的人这两年应该都有同一个感觉手头的潮流计算工具越来越“算不准”了。以前算一个纯负荷配电网前推回代跑一遍节点电压、网损清清楚楚。但现在网里接了光伏、风机、储能变流器同样的系统、同样的负荷实测电压谐波畸变率动不动就冲到5%以上而传统潮流程序根本算不出这些分量。问题就出在——传统潮流只算基波而DG接入后注入的谐波电流已经在网络里形成了不可忽略的谐波潮流。这篇文章就围绕“含DG的33节点配电网谐波潮流计算”这个课题展开用最经典的IEEE 33节点系统做载体聊清楚三个问题谐波潮流和基波潮流到底差在哪、DG的谐波源怎么建模、以及从零搭建一个谐波潮流仿真程序时哪些环节最容易翻车。适合正在做配电网谐波分析、分布式电源并网影响评估、或者写相关方向毕业论文的同行参考。1. 基波潮流能算清楚为什么还要做谐波潮流1.1 含DG配电网的谐波问题从哪来先把谐波的来源理清楚。配电网里新增的DG大体分两类旋转电机类小型水轮、燃气轮机直连和电力电子接口类光伏逆变器、双馈风机变流器、储能PCS。如今实际工程里占绝对主导的是后者而电力电子接口恰恰就是配电网里最主要的谐波源。逆变器产生谐波的机理不复杂核心就是开关动作。以光伏并网逆变器为例IGBT按PWM载波频率高频通断输出电压是一系列宽度不同的脉冲方波傅里叶展开后除了基波还含有大量的高次分量。谐波主要分两段低频段以5、7、11、13次为主这是三相桥式结构固有的特征次谐波一般和死区效应、调制策略有关高频段则集中在开关频率及其整数倍附近频率动辄几千赫兹到几十千赫兹对电能质量的影响方式和低频段不同。传统潮流程序把这些全部忽略只解50Hz下的电压电流分布。这在DG渗透率低的年代问题不大因为谐波电流幅值很小。但现在单台光伏逆变器容量做到几十千瓦甚至上百千瓦一个台区下多台逆变器并联谐波电流叠加后完全能推高公共连接点的电压畸变率这就必须把谐波潮流纳入计算。1.2 潮流计算在高次频域的盲区在哪儿搞清楚谐波潮流的本质一句话总结在网络的每个节点上把基波、5次、7次、11次各次谐波分量的电压电流关系都解出来。基波潮流是解一组非线性方程因为功率与电压的平方成正比而谐波潮流不同各次谐波下的网络是线性网络谐波源注入的谐波电流可以看作已知量求解过程就是解一组复数线性方程。通常的做法是在非线性元件处建模为谐波电流源或诺顿等效电路然后各次谐波独立求解。这带来一个关键好处基波潮流和谐波潮流可以解耦。先求基波得到每个节点的基波电压和DG的基波工作点再根据DG的谐波发射特性计算它注入的谐波电流最后在网络导纳矩阵上叠加求解。但这里有一个传统潮流完全没覆盖的问题网络元件在不同频率下的阻抗特性。线路电抗随频率线性增长变压器漏抗也随谐波次数线性增加负荷的谐波阻抗更复杂。如果沿用基波潮流里的恒定阻抗模型谐波计算结果会严重失真。这就是为什么不能简单把基波潮流程序扩展几个频点就完事。2. 含DG配电网谐波源建模别把逆变器当成理想电流源2.1 逆变器谐波发射特性与常见等效模型做谐波潮流计算第一步就是给DG建一个合理的谐波源模型。这里我见过不少初学者踩坑直接把逆变器的谐波含量按典型值设一个百分比然后当恒定电流源注入网络。这种粗暴做法在渗透率低的时候勉强能看但只要DG接入点离母线远一点或者系统里有无功补偿电容结果就完全不可信。实际工程中逆变器的谐波发射特性与运行工作点强相关。同样一台逆变器输出功率从20%升到80%谐波电流幅值和频谱占比都在变。原因是开关时刻的死区效应、控制环路的非线性都会随输出电流变化。所以严格的做法是用诺顿等效电路一个谐波电流源并联一个谐波阻抗谐波电流幅值根据基波工作点查表或拟合得到。具体建模时有几种常用方案建模方案基本原理适用场景恒流源模型忽略谐波阻抗直接注入固定幅值相位的谐波电流谐波源远离谐振频率初步估算诺顿等效模型谐波电流源并联阻抗谐波阻抗用逆变器输出阻抗近似考虑DG与网络阻抗交互的场合时域-频域转换在Simulink里搭详细开关模型采稳态波形后FFT分解谐波频谱复杂、需要高精度验证从我的经验看做33节点系统级的谐波潮流研究诺顿等效模型是性价比最高的选择。它比恒流源多了一个关键信息DG接入点的谐波阻抗。当网络阻抗与DG谐波阻抗发生交互时实际注入网络的谐波电流会被放大或衰减恒流源模型完全无法描述这种现象。2.2 网络各元件的谐波阻抗参数怎么取谐波潮流和基波潮流最大的差别就在各个元件的频率特性上。33节点系统里主要涉及线路、负荷和DG接口三类元件它们的谐波模型我逐一说明。线路的谐波阻抗是最容易处理的也是误差最容易失控的地方。理论上架空线和电缆都有分布参数效应但33节点属于中压配电网线路长度短、电压等级低用集中参数π等效已经足够。关键在于电阻的频率特性集肤效应导致谐波频率下电阻增大。工程近似公式是R(h) R(1) × (1 0.646h² / (192 h²))。5次谐波时电阻约增加5%25次时增加约17%不能忽略。电抗则简单X(h) hX(1)。负荷的谐波阻抗是另一个容易踩坑的地方。负荷本身的物理构成五花八门照明、电机、电子设备混杂。谐波潮流里面不可能逐个建模通常用并联R-L等效即认为负荷谐波阻抗的电阻和电抗都与频率相关。更稳妥的做法是参考CIGRE和IEEE工作组推荐的负荷模型恒阻抗模型(ZIP中的Z部分)在谐波频率下保持阻抗不变这也意味着谐波电流与谐波电压成正比网络呈线性特性。在33节点算例里负荷取恒阻抗模型是基准做法实测吻合度不错。变压器在33节点系统里通常位于首端电源侧。谐波频率下变压器的励磁支路感抗随频率升高但漏抗才是主导漏抗与谐波次数近似成正比。如果研究范围包含变压器后的部分要在导纳矩阵中把变压器的谐波阻抗折算进去。由于33节点系统经常假定变压器出口就是平衡节点这个效应容易被忽略但只要DG容量大变压器压降上的谐波分量就不能不管。2.3 DG接入位置与容量对谐波传播的影响谐波源建模完成后紧接着要考虑DG接入位置和容量的影响。很多人以为谐波畸变和DG容量成正比接入点越远畸变越大实际上这个直觉并不完全对。配电网的谐波传播有个很有意思的现象同一套网络拓扑DG接在不同节点同一容量的谐波注入在公共连接点产生的电压畸变可能相差好几倍。原因是谐波电流在网络上传播时会遇到线路电感和并联电容无功补偿、电缆对地电容形成的谐振回路。以33节点系统为例末端节点(如节点18)接入大容量DG时谐波电流要流经更长的馈线才能到达上一级母线线路阻抗大谐波电压降明显但DG接入点如果恰好在某个无功补偿电容器附近并联谐振频率落在5次或7次附近谐波电流就会被放大局部电压畸变率飙升。所以在选择DG接入位置时不能只看基波潮流的电压支撑作用必须做谐波扫描画一下从DG接入点看进去的谐波阻抗频率特性曲线看有没有谐振峰落在特征谐波频率上。容量方面DG容量增大不仅导致注入谐波电流绝对值增大更复杂的是谐波电流的相位也会随输出功率变化。两台容量相同的逆变器如果控制策略不同谐波电流在并网点可能抵消也可能叠加。这就是后面第5部分要讨论的多DG叠加效应问题。3. IEEE 33节点系统改造从基础参数到DG接入方案3.1 33节点系统的标准参数与模型搭建要点IEEE 33节点系统是配电网研究领域最经典的测试馈线由美国电气电子工程师学会提出基准电压12.66千伏基准功率10兆伏安总共33个节点、32条支路额定总负荷约3.715兆瓦加2.3兆乏。系统呈辐射状结构首端节点1连接上级变电站被当作平衡节点处理其余节点都是PQ节点。这个系统的数据在各类论文里已经很标准唯一需要注意的是各版本数据略有差异。不同文献中支路阻抗数据单位不同有名值/标幺值搭建模型前务必统一归算到标幺制。我的做法是取基准容量10兆伏安、基准电压12.66千伏把所有支路阻抗、负荷功率全部折算成标幺值再输入到程序中。这样后续做谐波分析导纳矩阵的迭代逻辑更统一。系统还有5个联络开关分别位于支路8-21、9-15、12-22、18-33、25-29平时处于断开状态。做基波潮流的标准算例时联络开关保持断开这是最基本的辐射状运行方式。如果研究DG接入后的故障重构或优化调度这些联络开关才需要操作。谐波潮流计算一般也不需要考虑联络开关闭合的情况但有一点要特别注意如果闭合了联络开关网络从辐射状变成弱环状前推回代法的迭代公式要改成回路电流法或者用牛顿法不能直接套用辐射状的公式。3.2 DG接入节点与容量的合理配置思路在33节点系统里接DG第一个问题就是接在哪儿。我建议先做基波潮流观察各节点电压分布情况。原始33节点系统在满负荷下末端节点电压已经偏低节点18的电压幅值大概在0.9标幺值附近。这种弱末端正是DG接入的典型场景既能支撑电压又能降低网损。我的算例里选择三个不同位置接入DG分别代表三种典型场景接入节点场景代表接入容量说明节点18馈线末端400千瓦电压支撑效果最明显但谐波传播路径最长节点8馈线中段300千瓦谐波流向复杂可能影响上游节点节点25分支线末端200千瓦分支线阻抗大局部畸变可能突出容量选择遵循一个原则单点DG容量不超过该节点峰值负荷的60%-80%。这个比例是从工程实际沿用过来的因为DG容量过大时除了谐波问题还有倒送功率导致馈线电压越限的风险。在谐波计算里容量还直接影响逆变器工作点进而影响谐波电流注入量所以不能随意放大。DG的类型我建议统一用光伏逆变器模型原因很简单光伏是当前配电网渗透率最高的DG类型谐波特性数据5、7、11、13次典型含量容易获取结果有工程参考价值。如果要研究风机变流器谐波频谱会随风速波动而动态变化计算复杂度上了一个台阶不适合作为入门算例。4. 谐波潮流计算的两种主流算法与我的选型建议4.1 频域解耦法的实现流程目前工程上最常用、也最容易实现的谐波潮流算法是频域解耦法。这个名字听起来复杂核心思想却很朴素基波和各次谐波相互独立分步求解。第一步用前推回代法计算基波潮流得到所有节点的基波电压幅值和相角。这一步和普通潮流计算没有任何区别前推回代在辐射状配电网里收敛快、内存占用小两三次迭代就能达到很高的精度。第二步根据基波潮流结果确定各DG的运行工作点然后从谐波频谱数据中查取或计算注入电流。这里要注意的是DG注入的谐波电流大小取决于DG输出功率和无功功率所以必须先有基波潮流结果不能凭空设定谐波电流幅值。第三步对每一个关注的谐波次数分别构建谐波导纳矩阵。节点导纳矩阵的构建方法和基波一样只是每个元件的导纳值按谐波频率重新计算。我自己写实现时直接生成一个函数输入谐波次数输出33乘33的复数导纳矩阵所有线路、负荷、变压器都按频率折算。第四步解复数线性方程组。这一部分的代码如下核心就是矩阵左除% 谐波次数列表配电网重点关注奇次谐波 H_list [5, 7, 11, 13]; V_harmonic zeros(33, length(H_list)); for k 1:length(H_list) h H_list(k); % 构建h次谐波节点导纳矩阵 Y_h buildY_harmonic(h); % 获取DG在h次谐波下的注入电流列向量33节点 I_h getDGharmonicCurrent(h); % 求解 h次谐波节点电压 V_harmonic(:, k) Y_h \ I_h; end % 计算各节点总谐波畸变率THD% V1_abs abs(V_base); % 基波电压幅值 THD sqrt(sum(abs(V_harmonic).^2, 2)) ./ V1_abs * 100;第五步计算总谐波畸变率和其他指标分析结果。这一步最常被忽略的是谐波电压的相角信息。打印报告时当然只关心幅值但在分析多DG谐波叠加时不同DG谐波电流的相位关系直接决定叠加结果所以解方程时务必保留复数结果不要一开始就取模。4.2 前推回代与统一牛顿法的对比除了解耦法谐波潮流还有一类统一求解法把基波和各次谐波方程联立用牛顿法同时求解。两类方法各有适用场景我放在一个表里对比对比维度频域解耦法统一牛顿法实现难度低复用基波潮流程序即可高需要推导统一雅可比矩阵计算速度快各次谐波独立求解慢方程组规模大收敛性稳定谐波部分为线性方程对初值敏感容易发散精度忽略谐波对基波的影响考虑基波与谐波耦合适用场景中低压配电网、DG渗透率适中含非线性负荷比例高、需要高精度边界工程项目的个人建议很明确33节点这种规模的中压配电网用解耦法完全够用而且调试成本低得多。统一牛顿法引入基波和谐波的耦合关系精度确实更高但这个耦合效应在中压配电网里通常很弱——谐波电压畸变率在5%以下时谐波对基波潮流的修正量几乎可以忽略。除非研究的是谐波畸变率高达20%以上的极端工况否则解耦法已经满足工程精度。4.3 迭代收敛性处理的实操经验实际的谐波潮流代码里最容易翻车的不是算法本身而是迭代过程不收敛。我在这部分踩过几个具体的坑分享出来供参考。第一个坑是谐波电流源幅值固定导致的收敛振荡。如果DG的谐波电流模型没有考虑节点电压对其的影响在节点阻抗较大、谐波电压较高时计算出来的谐波电压可能比基波电压还大显然不符合物理实际。解决方法是给谐波电流源加上电压修正系数典型的做法是计算一次谐波电压后按该电压修正注入电流进行两到三次迭代修正for iter 1:3 V_harmonic(:, k) Y_h \ I_h; % 根据谐波电压修正DG谐波电流考虑逆变器直流侧电压限制 I_h modifyHarmonicCurrent(I_h, V_harmonic(:, k)); end第二个坑是基波潮流不收敛导致谐波无法计算。33节点系统在极端DG渗透率下可能倒送功率前推回代法在这种工况下收敛变慢。解决办法是给迭代变量加一个阻尼系数把电压修正量限制在合理范围内或者改用牛顿法求解基波潮流。第三个坑是频率选择。谐波次数取到多少才算完整基础算例取5、7、11、13四次即可这是三相桥式逆变器最常见的低频特征谐波。做精细分析时可以扩展到25次、甚至50次但要注意高频段集中参数模型误差增大结果参考价值下降。更实际的处理方式是先算基波计算节点电压畸变率如果某次谐波占总畸变的比值低于1%就可以忽略。5. 仿真结果怎么看电压分布、THD与谐振风险5.1 典型算例结果数据解读我按第3部分的接入方案跑了一组算例节点18接入400千瓦DG节点8接入300千瓦DG逆变器谐波电流按六脉波典型频谱设定基波潮流用前推回代法谐波频域解耦法求解。结果能用但也正好暴露了一些容易被忽视的问题。先看基波电压分布。接入DG前节点18末端电压约0.90标幺值接入400千瓦DG后抬升到0.95标幺值左右电压支撑效果立竿见影。节点8的300千瓦DG对首端电压影响很小但明显改变了局部馈线的功率流向。再看谐波结果以节点18处的DG为例它注入网络的5次谐波电流在流回上级电网的过程中经过节点17到节点18的长馈线线路上产生明显的谐波电压降。这导致一个现象DG接入点附近节点的5次谐波电压畸变率反而不如线路中间某个节点高。谐波电压最大的位置不在谐波源所在节点而在谐波源与系统阻抗之间某个网络阻抗匹配的位置。这个结论对谐波治理方案的制定很关键——单点治理比如在DG出口加滤波器可能无法降低最严重节点的畸变率。5.2 谐波放大的识别与定位实际算下来节点18接入DG的场景有一个典型的谐波放大现象7次谐波电压畸变率明显高于5次。初看很反常因为逆变器注入的5次谐波电流通常大于7次。查了谐波阻抗频率特性后发现节点18看进去的7次谐波阻抗明显大于5次谐波阻抗同样的注入电流产生了更大的电压降。这说明谐波潮流计算必须关注网络阻抗的频域特性。在33节点系统里电缆或架空线的对地电容不大但如果有并联无功补偿电容器组就必须做谐波阻抗频率扫描。我用一个简单函数从DG接入点看进去计算从1次到50次的等效阻抗幅值Zin zeros(1, 50); for h 1:50 Y_h buildY_harmonic(h); % 从节点n看进去的等效阻抗 Zin(h) 1 / Y_h(n, n); end plot(1:50, abs(Zin));如果某个频率下阻抗幅值出现尖峰就说明该频率存在并联谐振风险。当谐振频率与逆变器特征谐波5、7、11、13次重合时即使注入电流很小谐波电压也会被显著放大。这种情况下最有效的解决方法是调整DG接入位置或改变无功补偿容量改滤波器参数反而可能把谐振点推向另一个特征谐波频率。5.3 多DG协同接入时的叠加效应多台DG同时接入时谐波叠加问题在33节点算例里表现得非常直观。单个DG接入时各节点THD都低于3%似乎一切正常但节点8和节点18同时接入后部分节点THD不是简单的算术相加有些节点5次谐波反而减弱有些节点7次谐波明显增强。原因在于谐波电流的相位。逆变器的谐波电流相角取决于基波电流幅值、锁相环角度和控制策略两台DG的基波工作点不同谐波电流相位也不同。当它们的5次谐波电流相位接近同向叠加后电流增大相位相差较大时叠加后电流反而减小。这个现象给了工程上一个重要启示在规划多个DG接入时单纯看各台逆变器的谐波发射指标是不够的必须做联合仿真。不同厂家、不同控制参数的逆变器谐波相位特征差异很大实际工程中偶尔会遇到后接入的DG让全网谐波降低的现象原理就在于此。6. 工程应用中的几个坑和排错思路6.1 谐波次数选取与精度取舍谐波潮流计算中谐波次数的选取是一个工程权衡问题。取少了漏掉关键谐波分量取多了计算量上升但精度未必提高。我建议按以下优先级选取奇数特征次谐波5、7、11、13、17、19必取这是逆变器最显著的分量偶次谐波2、4等在故障态或不对称工况下才明显正常对称工况可以忽略三倍频谐波3、9、15等零序分量需要考虑变压器接线方式若变压器为三角形接线则不会向高压侧传播可以忽略。对33节点算例来说一个非常实用的方法是分两次计算。第一次只算5、7、11、13次观察THD结果如果某个节点的THD已经接近标准限值再扩展到25次做精细核算。标准限值参考IEEE 519或GB/T 14549低压和中压电网电压THD限值一般为5%单次谐波通常限制在3%左右。6.2 负荷与线路参数的敏感性谐波潮流计算结果对负荷模型的敏感性极高。我在试验中换了两种负荷模型恒阻抗模型和恒功率模型分别计算结果THD差异最高达到30%以上。这背后是物理机制的区别恒阻抗模型下负荷阻抗随谐波频率变化等同于固定阻抗值谐波电流随谐波电压线性变化恒功率模型下谐波电流计算要迭代修正谐波会改变基波工作点。我对做算例研究的建议是主结果采用恒阻抗模型然后补充一组恒功率模型的敏感性分析说明模型选择对结论的影响。这样写论文或做工程报告更有说服力也避免了拍脑袋选模型的质疑。线路参数方面最容易出问题的是忽略集肤效应导致的电阻频率特性。我特意试过在5次谐波时忽略集肤效应THD结果大约偏低5%-10%在25次谐波时忽略误差可以扩大到20%以上。所以高频谐波分析必须计入电阻频率修正。6.3 从算例走向实际配电网的注意事项33节点系统是一个理想化很强的测试馈线从它推广到实际配电网时有几个经验值得提前说明。第一实际配电网的三相不平衡度不可忽略。33节点算例通常做三相平衡假设但台区配电网天生三相不平衡负序电流会产生3次及倍频零序谐波计算必须扩展到三相模型。单相DG如户用光伏逆变器接入不同相序谐波传播特性差异很大。第二实际逆变器出厂时都内置了滤波器和并网标准控制谐波发射水平远低于早期文献里的六脉波理论值。在算例里可以设一个比较高的谐波含量来分析最恶劣工况但实际评估时还是要以设备实测数据或厂家提供的谐波频谱为准。第三从算例走向实勘时最关键的验证手段是实测。在DG接入点装一台电能质量分析仪连续记录一周以上的谐波数据会看到谐波随日照、风速、负荷变化的时变规律。理论仿真算的是稳态工况而实际谐波是时变的两者对比时必须保持相同的运行工况。在33节点系统上把谐波潮流计算做扎实其实就是在练一套可迁移的方法论DG谐波源建模、网络谐波阻抗构建、频域求解、结果分析。这套流程换个实际馈线换个DG类型逻辑完全一样差别只在于参数和规模。我自己的体会是初学阶段在33节点系统上把每一个环节都手动实现一遍比直接套用商业软件更有价值。因为商业软件把导纳矩阵的构建过程封装成了黑盒子一旦计算结果异常你根本没法判断是算法问题、参数问题还是模型问题。手动实现虽然前期费时间但换来的是对每个数值背后的物理意义都心里有数。这个账怎么算都划算。