分布式电源接入配电网影响仿真是这几年电网规划、新能源并网、电力系统方向研究生和一线工程人员绕不开的课题。我一个项目里用 Matlab/Simulink 把 10kV 馈线模型搭起来逐一对比分布式电源在不同位置、不同容量下的电压分布、短路电流、谐波和谐波畸变率然后把影响量级量化出来。整个过程从建模到参数计算再到结果提取踩了不少坑。这篇博文把完整的技术路线和实操细节写出来特别是那些仿真文档里不会写的经验希望能帮你少走弯路让复现这个题目时心里更有底。1. 分布式电源接入配电网影响研究的核心思路拆解1.1 为什么“接入一个电源”会引发连锁反应传统配电网是单电源辐射网潮流方向固定从变电站母线流向用户。这时候电压沿馈线逐步降低保护按单端电源整定短路电流计算也简单。分布式电源Distributed GenerationDG接入后馈线上出现多个电源点潮流方向可能倒转故障时DG会向短路点注入电流。任何一个环节被打破整体评估方法就得重写。很多人一开始容易低估这件事的复杂度觉得“不就是多并一台光伏逆变器嘛”。实际工程里分布式电源接入影响评估涉及四个层面电压水平、短路电流水平、谐波电能质量、继电保护配合。这四个层面互相关联不能孤立看。比如电压抬升可能触发逆变器过压脱网脱网后又导致馈线末端电压过低DG容量增加抬高短路电流又可能超出开关设备开断能力甚至影响保护灵敏度。从研究角度最需要先想清楚的是评估目标是什么。如果是电网规划侧重点在电压和短路电流如果是并网方案设计谐波和孤岛响应更关键如果是继电保护整定重点看故障电流大小和方向变化。目标不同模型复杂度、仿真步长、调节方式都不同。这个项目定制化程度很高但分析框架是通用的。1.2 影响机理前先分清“源”的种类分布式电源不是单一概念光伏、储能、风电、燃机、燃料电池都属于分布式电源但它们的并网接口和故障特性差异很大。最根本的区别在于旋转型电源同步发电机、异步风机和逆变器型电源光伏、储能、部分风机。旋转型电源接在配电网里故障时能提供持续短路电流相当于给系统增加了一个额外的强电源支路短路电流贡献不容忽视。逆变器型电源情况不同逆变器本身有过流限幅控制输出电流一般限制在额定电流的 1.2~1.5 倍而且控制响应速度很快。这意味着光伏并网点的短路电流提升幅度远没有同步机型DG那么大但谐波含量和电压波动问题更突出。项目中我做的是 10kV 配电网DG 以光伏为主考虑到未来储能接入也兼顾了逆变器型电源的通用特性。模型处理时把光伏分成两部分直流侧光伏阵列 MPPT DC/DC和交流侧三相逆变器 滤波 并网控制。如果只做潮流和短路分析可以简化成 P、Q 注入源如果要做谐波和动态分析必须用开关级模型。这个取舍后面再展开。2. Matlab 仿真方案设计从潮流计算到电磁暂态的路线选择2.1 用 m 脚本还是 Simulink 模块取决于你要什么精度Matlab 做配电网研究两条路线并行不悖一是纯 m 脚本编写潮流计算、短路电流计算程序优点是透明可控、方便批量扫描二是用 Simulink 搭建模型做动态仿真优点是逆变器控制、保护动作这些细节可以真实还原。我这批研究是混合路线。涉及多节点电压分布、DG容量和位置扫描用 m 脚本写牛顿-拉夫逊潮流涉及逆变器控制策略对电压闪变和谐波的影响用 SimulinkSimscape Electrical 库里的三相元器件搭时间域仿真。两个结果互相校验很多问题互为镜像。举个例子分析 DG 接入对电压分布的影响时我写了一个 m 脚本内置 IEEE 33 节点配电网参数潮流计算采用牛顿-拉夫逊法。改变 DG 接入节点编号和注入功率循环跑几十种场景几分钟出所有节点的电压曲线。这个场景用 Simulink 跑就要搭几十个模型效率低太多。反过来分析逆变器谐波频谱时需要看开关频率附近的谐波分布这就必须在 Simulink 里搭 PWM 逆变器模型用 1e-5~1e-6 秒级步长做暂时变仿真再用 Powergui 的 FFT 工具分析电流波形谐波成分。2.2 模型颗粒度稳态、准稳态、电磁暂态的三档选择模型颗粒度决定了仿真结果的可信度和计算成本。我总结了三档第一档是纯稳态模型。线路用阻抗负荷用恒定功率DG 用 PQ 或 PV 节点代表。用在潮流计算、静态电压分析。优点是快适合几十节点网络反复测算缺点是看不到器件级动态。第二档是准稳态模型。在稳态基础上加入 DG 出力波动、负荷时序曲线用连续潮流或逐时仿真看一天 24 小时的电压变化。用于评估 DG 波动性给配电网带来的电压越限风险。第三档是电磁暂态模型。在 Simulink 中搭三相电路逆变器用带 IGBT 和 PWM 控制的开关级模型电缆、变压器用集中参数或分布参数模型。用于分析谐波、电压暂降、孤岛过程、保护动作行为。计算量最大但结果最真实。顺序上先做第一二档找到问题节点和越限场景再用第三档精确考察最严重的场景。千万别一开始就上电磁暂态模型调试成本会让你崩溃。2.3 整体仿真框架别让脚本和模型变成一团乱麻这类研究涉及大量场景组合很容易陷入“改一个参数就得重新跑一遍”的泥潭。我的做法是建三层结构数据层把配电网参数、DG参数、负荷曲线集中放在独立脚本或 .mat 文件里。检修时只改数据不动算法代码。 算法层潮流计算、短路计算、谐波分析各自封装成 function。比如潮流程序写成一个my_NR_powerflow(branch, bus, DGbus, DGpq)函数输入电网拓扑和 DG 注入返回节点电压和支路功率。 场景层循环脚本遍历 DG 接入位置、容量、渗透率等变量调用算法层函数批量收集结果。这个结构初看“绕”实操时是真省事。改一个 DG 接入位置不用翻代码只需把场景矩阵对应的行改掉。后面对比结果出图也只需要在场景层统一处理。3. 核心模型搭建与参数计算配电网友和目标 DG 建模细节3.1 基准配电网模型从零搭一套可复用的 10kV 馈线配电网仿真研究里IEEE 33 节点系统是天生的“标准试验场”几乎每篇同类文献都引用它。它包含 33 个节点、32 条支路总负荷约 3.7MW 2.3Mvar额定电压 12.66kV常用于 10kV 系统近似。更重要的是网上到处有支路阻抗、节点负荷的公开参数方便复现和横向对比。我先按 Linux 环境还是 Windows 无所谓M文件里用两个矩阵bus和branch存储系统数据。bus存储节点编号、负荷有功、无功branch存储首端节点、末端节点、电阻、电抗。这个工作量不大但必须反复检查一个参数抄错后面全乱。配电网建模几个要点线路阻抗需要区分单位阻抗和总阻抗。IEEE 33 节点原始参数是每条支路的总电阻、总电抗单位Ω不要乘线路长度。变压器和联络开关先简化初始网络不用闭合联络开关保持单辐射结构。负荷模型静态分析用恒定功率就够动态分析建议加一定比例的恒阻抗成分比如 30%防止电压过低时仿真发散。节点电压基准值取 10kV或 12.66kV按网络参数定功率基准值取 1MVA。潮流计算收敛后记录各节点电压标幺值。这一步是整个项目的地基地基不平后面盖多高都危险。3.2 光伏 DG 模型的三个层次PQ源、单机逆变器和等效机群光伏逆变器并网模型按项目需求分三个层次我建议按研究目标选好再做避免过度建模。层次一PQ 注入源。用潮流计算节点类型实现在 DG 接入节点把注入有功、无功设成常数。适用于电压分布、渗透率评估。光伏逆变器一般在单位功率因数运行无功设为 0但若研究无功支撑能力可把无功设为负值发出无功或正值吸收无功。层次二完整单机逆变器模型。Simulink 里用 PV 阵列模块根据光照、温度计算 I-V 特性、Boost 变换器带 MPPT 控制、三相两电平逆变器带 PQ 控制或电压电流双闭环。适用于谐波分析、动态响应研究。层次三等效机群模型。几十台户用光伏不可能逐台建模用一台容量等于总和的大逆变器等效。适用于调度和稳控研究。我这里重点用了层次一和层次二。层次一解决“接入位置/容量→电压分布→短路电流”这个主链条层次二专门验证最恶劣场景下的谐波和逆变器响应。3.3 关键参数计算那些不能“拍脑袋”的数参数计算是容易糊弄但绝不能糊弄的环节。我把几个关键参数的计算方法和取值依据列出来DG 容量选择常见做法是按渗透率DG 总容量 / 系统总负荷扫描。比如系统总负荷 3.7MW选择 10%、20%、30%、40%、50% 渗透率对应 DG 总容量 0.37MW、0.74MW、1.11MW、1.48MW、1.85MW。这个扫描能直观画出渗透率对电压和短路电流的影响曲线。逆变器额定容量光伏阵列直流侧峰值功率和逆变器交流侧额定容量之间需要留余量。工程上 DC/AC 容配比取 1.1~1.3即 500kW 光伏阵列配 400~450kVA 逆变器。仿真里如果只做并网研究为简化可设逆变器容量等于光伏峰值功率乘以 0.9 左右。注意逆变器容量影响短路电流限幅值不能随便取。并网变压器10kV 配电网里分布式电源通常经过 0.4kV/10kV 升压变压器并网。变压器阻抗电压 uk% 典型值取 4%~6%仿真里设 4.5%。这个数影响短路电流计算也影响谐波谐振特性不能省略。滤波器参数逆变器输出一般加 L 或 LCL 滤波器。L 的取值在 0.05~0.15 pu以逆变器容量为基准LCL 的谐振频率控制在 10 倍基波和 0.5 倍开关频率之间。我用的开关频率 5kHzLCL 谐振频率设计在 1kHz 左右避免在 50Hz 附近产生不利谐振。这些数据不是凭空定的设计每个参数前我都会先想它影响哪一项结果、数量级大概是多少。想不清的先做敏感性分析画曲线看趋势再定。4. 结果分析DG接入对配电网影响的量化评估4.1 电压分布接入位置比想象中更敏感我跑了 IEEE 33 节点基准模型无 DG 时馈线末端节点电压约 0.93 pu标幺值在合格范围但偏低。接入 0.5MW 光伏后如果接在靠近变电站的节点 2电压提升极小几乎看不到变化接入馈线中段的节点 8末端电压抬升到 0.95 pu接入接近末端的节点 18末端电压抬升到 0.98 pu但光伏接入点附近电压反而可能超过 1.03 pu。接入节点DG容量(MW)接入点电压(pu)末端电压(pu)电压最大偏移无DG00.970.93-7%节点20.50.990.94-6%节点80.51.020.952%/-5%节点180.51.050.975%/-3%节点181.01.081.028%数据说明两件事DG 接入位置靠近馈线末端电压抬升效果越显著容量超过一定值后局部电压可能越上限10kV 系统允许偏差一般为 7%/-10%即 10.7kV~9.0kV。规划和设计时要特别注意不能只看“有没有提升”还要看“提升过头没有”。分析电压影响还有一个维度就是 DG 出力波动。我在 Simulink 里给光照加了一天 24 小时曲线中午光伏大发时节点电压冲到最高傍晚出力下降后电压回落。一天内电压波动范围如果超过 2.5%对敏感负荷就是隐患。4.2 短路电流逆变器限幅特性给继电保护带来的新挑战配电网短路电流计算传统方法假定故障后系统电压恒定为 1 pu各电源的短路电流贡献按电抗换算。DG 接入后故障点附近新增电源支路短路电流水平整体上升。我用三相短路仿真看了几个典型故障点故障位置无DG短路电流(kA)接入0.5MW DG后(kA)接入1MW DG后(kA)节点83.23.64.1节点182.12.83.5DG 容量增加短路电流贡献呈近似线性增长但幅度有限因为逆变器限幅。旋转型 DG 就不一样同步电机短路电流可能达到额定电流的 6~10 倍这才是对开关设备开断能力的真正威胁。分析时一定要区分 DG 类型不能全部按恒压源处理。短路电流的工程意义在于校核馈线开关短路开断电流是否足够保护定值是否需要调整DG 下游故障上游保护能否正确动作计算得到的短路电流超出设计值就要考虑限流措施如串电抗或调整运行方式。4.3 谐波与电能质量多台逆变器叠加后的“脏电”风险单台逆变器并网电流 THD 一般控制在 3%~5% 以下满足国标并不难。但多台逆变器集中在同一条馈线谐波叠加效应不容忽视。我在 Simulink 里把 10 台 100kW 逆变器分散接入 15 个节点馈线首端测得的电流 THD 从单台的 2.8% 上升到 4.7%特定次谐波5 次、7 次放大明显。谐波分析要区分两个位置的指标节点电压 THD 和支路电流 THD。国标对公用电网谐波电压有明确限值10kV 通常要求 THD≤4%对谐波电流按接入容量分配限值。仿真结果用来判断是否需要加装滤波装置或者调整逆变器控制策略比如引入谐波补偿环。这里特别提醒一下谐波谐振问题。配电网中存在电缆电容、变压器电感和并联电容器在特定频率可能形成并联谐振。当逆变器谐波注入频率接近谐振频率时谐波电流被放大电压波形畸变加剧。仿真时用阻抗扫描扫频功能画出节点阻抗随频率变化的曲线找到谐振峰判断风险频率是否落在逆变器主要谐波频段内。4.4 孤岛保护与重合闸配合仿真的“隐藏雷区”DG接入对配电网保护的影响很多人做电压和短路电流分析时会忽略孤岛问题。简单说配电网故障导致主网跳闸后如果 DG 仍在发电部分负荷由 DG 独立供电就形成非计划孤岛。此时电压和频率失控检修人员安全受到威胁所以并网规程强制要求防孤岛保护。仿真里我做了两种孤岛场景一是主网断路器跳开DG 带本馈线负荷继续运行观察电压频率偏移二是无防孤岛策略和带主动式防孤岛策略主动注入扰动的对比。结果显示被动式防孤岛依靠电压、频率检测在 DG 出力与负荷功率基本匹配时存在检测盲区主动式防孤岛策略能有效缩短检测时间。这个结论对工程实践有直接参考价值也提醒做仿真研究时别只盯着电压短路电流保护配合必须纳入分析框架。5. 常见问题与调试经验那些文档里不写的坑5.1 潮流计算不收敛多半不是算法问题牛顿-拉夫逊潮流不收敛初学者第一反应是改迭代次数、调容差但真正原因往往出在数据和初始值上。我踩过的几个典型坑负荷数据单位不统一。IEEE 33 节点原始数据里有功单位是 kW无功是 kvar如果在程序中按 MW、Mvar 处理需要除以 1000。单位错了负荷太大潮流直接发散。PV 节点无功越限。在 m 脚本里把 DG 设成 PV 节点定电压定有功但 DG 无功出力超出逆变器容量限制时程序迭代会震荡。解决办法是加无功越限约束当无功越限时自动转成 PQ 节点并记录越限状态。初始电压给得太偏。配电β网负荷重的节点初始电压给 1.0 pu迭代几步就飞了。我的习惯是初始电压全部给 1.0但电容和电抗比值大的线路会敏感可以先试 0.98 或直接用平值启动再逐步加载。5.2 Simulink 仿真发散先查数值解算器Simulink 里搭完并网逆变器一跑就发散经验是八成出在三个地方解算器步长过大。电力电子模型开关频率 5kHz步长至少取开关周期的 1/20也就是 1e-5 秒量级。用变步长解算器时把最大步长限制在 1e-5否则容易错失开关动作时刻。控制环 PI 参数不合。逆变器电流环 PI 参数直接照抄文献容易翻车因为不同容量、不同滤波参数下电流环增益差异很大。我的做法是先断开并网开关单独调内环让电流跟踪给定值正常再闭合并网开关调外环。直流侧电压源内阻过高。PV 阵列模块输出特性是电流源型若并联电容太小MPPT 扰动时直流电压波动大容易触发保护退出。直流母线电容按经验取 3000~5000uF/MW至少保证直流电压纹波小于 2%。5.3 结果可视化与对比让数据“会说话”仿真做完一堆数据不会展示等于白做。我用几个固定的可视化套路出图快且能直接放到报告里电压分布图横轴节点编号纵轴电压标幺值把无 DG 和不同 DG 接入场景画成一组对比曲线越限范围用阴影标注。这张图最能直观体现 DG 的影响。短路电流柱状图按故障点分组无 DG、有 DG 两组染色对比误差棒可省但数据标签必须标清楚。谐波频谱图用 FFT 分析得到电流频谱把 50Hz 基波幅值归一到 100%各次谐波用百分比柱状图展示。50Hz 和 2.5k~5k 开关频率附近会出两簇谱峰标注出来就是很好的结论。扫描曲线改变 DG 渗透率横轴渗透率、纵轴电压偏差或短路电流得到趋势线一眼看出安全边界在哪里。出图时注意坐标轴范围和标注清晰调好线宽推荐 2pt 以上、对比色推荐色盲友好的蓝红和图例位置出图后缩小到 1/4 尺寸依然清晰再导出。6. 这个题目的扩展方向与我的实操建议项目做到后面你会发现配电网友 DG 影响分析有大量值得深入的方向。如果时间允许我建议按以下顺序扩展先把主场景跑扎实。无 DG 基线、不同容量、不同位置、不同渗透率的电压和短路电流矩阵全部整理出来。这是毕业论文或技术报告的核心支撑。再补动态场景。光伏出力波动、负荷波动下电压越限的概率性分析。可以用某地全年的光照和负荷数据做时序仿真统计越限小时数。这个结论对实际工程最有用。然后做保护配合。用 Simulink 搭馈线保护过流/方向过流模拟 DG 接入前后保护动作行为分析保护误动、拒动的边界条件。这个方向发论文很容易出亮点。谐波治理和逆变器控制优化也是热门方向。在 LCL 滤波器参数设计、谐波补偿环、虚拟同步机控制策略里选一个仿真对比改进前后的 THD 和电压支撑效果。我个人的实操体会是做这类研究模型的“简洁正确”远比“复杂炫技”重要。很多初学者一上来就搭全三相开关级模型跑一个场景要几十分钟结果还不收敛最后项目卡壳。我的建议是先用潮流算法把静态问题摸透再按需上电磁暂态模型验证关键场景。分层验证、逐步复杂化效率最高。最后分享一个小技巧在整个项目开始前花一天时间把所有仿真结果的图表命名规范、存储目录结构定好。比如fig_contour/voltage_by_DG_position.png、data_export/case_18_500kW.mat。等做完整批扫描你会发现这个习惯帮你省下至少三天整理数据的时间。项目文件不混乱研究结论才站得住脚。