基于Matlab的分布式电源接入配电网影响分析程序设计与实现

1. 项目概述与整体思路搞配电网的人这两年应该都有个明显感觉分布式电源不再是“锦上添花”的示范工程而是实打实批量接入的常规业务了。光伏屋顶、村级扶贫电站、工商业储能、分布式风电一个10kV馈线上挂几个甚至十几个并网点已经是常态。但接入归接入电网公司做接入方案评审时问得最多的问题永远是那几个短路电流会不会超标电压会不会越限继电保护还能不能配合而这些问题的答案靠手算或者经验估基本没法给出让人信服的结论必须用仿真程序实打实跑一遍数据。这就是我做这套基于Matlab的分布式电源接入对配电网影响分析程序的初衷。这套程序核心解决的就是三件事潮流计算、短路电流计算、以及接入前后关键指标的对比分析。研究对象是典型的10kV配电网分布式电源按光伏、储能等不同类型分别建模可以设置不同的接入位置和接入容量程序跑完之后直接输出节点电压变化曲线、支路电流变化、短路电流水平以及不同渗透率下的趋势对比。不管是做接入方案设计、编写评估报告还是做毕业设计、科研论文的基础数据支撑这套程序都足够用。用Matlab来做这个研究而不是用PSASP或者ETAP这类电力系统专用软件原因很实在一是Matlab的编程自由度极高模型的拓扑结构、电源特性、故障类型都可以自己控制不会受商用软件黑盒模块限制二是Matlab的矩阵计算能力天然适合潮流迭代求解写起来代码量小、调试直观三是后续如果要扩展比如加入储能控制策略、做经济性分析Matlab生态里的工具箱和配套代码资源非常丰富。后面我会把整个程序的架构设计、数学模型、关键代码逻辑和调试过程中踩过的坑一一拆开讲希望能给正在做类似课题的同行一些能直接抄作业的参考。2. 配电网模型搭建与分布式电源建模2.1 为什么选IEEE 33节点模型作为基础网络做配电网仿真第一步就是选基础网络模型。我最终采用的是IEEE 33节点标准算例这个选择不是随手定的而是对比了多个典型模型之后确定下来的。IEEE 33节点系统是一个经典的三相平衡辐射状配电网模型基准电压12.66kV总负荷约3715kW加2300kvar共32条支路、5条联络开关支路。它之所以在配电网研究里被用得最多是因为结构上高度贴近实际10kV馈线的形态主干线长度适中、有分支线、负荷沿线路分布不均匀、末端电压偏低的现象明显。这个特征非常关键——分布式电源接入后对电压的支撑作用、对潮流的改变在这种不均匀负荷分布的网络上表现得最有说服力。如果换成负荷均匀分布的简单辐射网仿真结果会过于“理想”工程参考价值就打折扣了。程序里我会把IEEE 33节点系统的线路参数、负荷参数直接内置成基础数据矩阵这样后续修改接入方案时不需要改动网络本体结构只需在指定节点上追加分布式电源模型就行。这里有个细节值得注意原始IEEE 33节点算例的基准电压是12.66kV严格说是“中压配电网”的典型代表跟我们要研究的10kV系统存在一定电压偏差。我在程序里做了标幺值处理把基准电压统一换算成10kV线路阻抗按比例归算这样计算结果在数值上也能直接对应工程上的10kV场景。2.2 分布式电源的数学模型怎么选分布式电源接入配电网不同类型其外特性差异非常大。光伏逆变器在PQ控制模式下表现为恒功率源有功和无功按设定值出力而双馈风机在特定工况下可能表现出功率因数恒定的特性储能则根据控制策略的不同既可能充电相当于负荷也可能放电相当于电源。如果全部用同一个模型硬套仿真结果会失真。我的程序里把分布式电源抽象为三类可切换模型PQ恒功率模型给定有功出力PG和无功出力QG潮流计算中作为负的负荷节点处理。适用于光伏逆变器、全功率变流器并网的风机等经电力电子变换并网的电源。PI恒功率因数模型给定有功出力和功率因数无功随有功动态变化QG PG × tan(arccos(φ))。适用于要求功率因数恒定的并网逆变器控制模式。PV定电压模型给定有功出力和节点电压目标值无功作为待求量参与迭代。适用于具备自动电压调节能力的电源比如部分储能PCS或者调相机。实际使用中最常用的是PQ模型因为这个最符合当前分布式光伏“最大化消纳”的运行策略。程序里通过一个输入参数controlMode来切换默认填1就是PQ模型填2是恒功率因数填3是PV定电压。代码里我已经预留了模型扩展接口后续想加入跟网型/构网型逆变器的详细暂态模型直接在对应函数里增加分支就行。2.3 接入位置与容量的典型场景设置分布式电源接入影响评估光有模型还不够还得有一套科学的场景设置逻辑。我在程序里预设了三种典型的接入分析场景场景一不同接入位置对比——分别在馈线首端节点2、中段节点17、末端节点32接入相同容量的分布式电源分析电压分布和线损的变化规律。这个场景能直接回答“并网点选在哪里更合理”的问题。场景二不同接入容量渗透率对比——在同一个节点比如中段节点17分别接入容量为0.5MW、1MW、2MW、3MW的分布式电源观察配电网从“正常受电”到“潮流返送”的变化过程。场景三多点接入组合——在多个节点同时接入分布式电源模拟实际工程中一条馈线上多个分布式光伏同时并网的场景考察对短路电流和保护的叠加影响。这三个场景的设计逻辑是层层递进的先看位置的影响再看容量的影响最后看组合场景下的整体效应。程序里通过一个scenario参数控制运行后会批量生成对比结果。实际做接入评估时我通常会把这三个场景都跑一遍因为供电公司审查接入方案时关心的正是这几个维度。3. 核心算法原理与Matlab程序实现3.1 潮流计算为什么用牛顿-拉夫逊法分布式电源接入配电网影响分析的所有指标几乎都是建立在潮流计算基础上的。节点电压、支路电流、网络损耗全都要靠潮流求解得到。配电网潮流计算有两类主流方法一类是前推回代法另一类是牛顿-拉夫逊法Newton-Raphson Method。前推回代法实现简单、收敛性好特别适合纯辐射状配电网但它有个明显局限处理分布式电源接入后的“多电源”网络时如果网络中存在环网结构或者电源节点数量较多前推回代法的处理流程会变得繁琐尤其是有PV节点时无功迭代修正比较麻烦。而牛顿-拉夫逊法在求解多节点、多电源类型的非线性方程组时通用性更强收敛速度是二阶的虽然每次迭代需要计算雅可比矩阵但配电网节点规模一般就几十个点这个计算开销完全可以忽略。我在程序里采用极坐标形式的牛顿-拉夫逊法。核心思想是把潮流问题转化为求解一组非线性功率平衡方程ΔPi Pgi - PLi - Ui∑Uj(Gijcosθij Bijsinθij) 0 ΔQi Qgi - QLi - Ui∑Uj(Gijsinθij - Bijcosθij) 0其中ΔPi和ΔQi就是功率不平衡量牛顿法的目标就是通过迭代修正节点电压幅值U和相角θ使这些不平衡量趋近于零。程序实现中每次迭代先算不平衡量再解修正方程J矩阵的线性方程组然后更新状态变量重复这个过程直到收敛精度满足要求。3.2 程序主功能的模块划分整个Matlab程序的代码架构我按照功能划分成五个模块每个模块独立一个函数文件。这样做的最大好处是修改某个功能时不需要全局翻代码而且单个模块可以单独测试。第一个模块是数据准备模块init_network.m负责加载IEEE 33节点系统的拓扑参数、线路阻抗数据和负荷数据生成导纳矩阵。这个模块是整个程序的地基数据一旦出错后面全错所以我会在模块内加入参数合法性校验比如线路阻抗不能为负、节点编号必须在有效范围内出错直接抛出明确提示而不是往下运行到最后才莫名报错。第二个模块是分布式电源模型模块dg_model.m输入DG类型、控制模式、有功出力、无功出力或功率因数等参数输出该DG节点在潮流计算中对应的等效注入功率。前面提到的PQ、PI、PV三种模型就在这里实现通过switch语句分流处理。第三个模块是潮流计算模块powerflow_nr.m这是核心求解器实现牛顿-拉夫逊迭代过程。输入网络数据、DG参数和负荷参数输出节点电压幅值、相角以及支路潮流、网络损耗。为了避免迭代不收敛时程序无限死循环我设置了最大迭代次数一般是30次超过就报错并给出逐步缩小的功率不平衡量变化曲线方便定位问题。第四个模块是短路电流计算模块shortcircuit_calc.m采用对称分量法对指定节点施加三相短路故障计算故障前后的短路电流水平。这个模块的关键是形成正序、负序、零序网络并处理DG注入对短路电流的贡献。分布式电源对短路电流的影响机理和常规同步发电机不同逆变器型电源的短路电流受变流器限流控制约束一般只有额定电流的1.2到1.5倍这个在程序里通过DG故障电流系数来体现。第五个模块是结果输出与对比模块report_results.m汇总潮流计算和短路电流计算的原始数据支持三种可视化输出节点电压分布对比图、支路有功损耗对比柱状图、不同DG容量下的短路电流趋势曲线。所有图像自动保存为高分辨率PNG方便直接粘贴到评估报告里。3.3 短路电流计算的实现细节配电网短路电流计算这部分是很多人最容易出问题的地方。分布式电源接入前后短路电流水平的差异直接关系到开关设备的开断能力校核和继电保护定值整定。程序里我采用的是IEC 60909标准的简化工程方法忽略负荷电流只考虑故障前的空载运行状态发电机和分布式电源用次暂态电抗逆变器电源用等效故障电流源建模计算三相短路时的对称短路电流初始值。具体实现时分成三步第一步形成节点阻抗矩阵Z这个是短路电流计算的基础。Z矩阵的物理意义是节点电压与注入电流的比值关系它的对角线元素Zkk就是节点k的自阻抗非对角线元素是互阻抗。用节点阻抗矩阵求短路电流有个简洁的公式在节点k发生三相短路时短路电流Ik Uk(0) / Zkk其中Uk(0)是故障前节点k的正常运行电压。这一步在Matlab里可以有两种做法一种是对导纳矩阵求逆直接得到另一种是用支路追加法逐条形成Z矩阵。对于33节点的网络直接求逆就行计算量很小如果网络规模大到几百节点建议改用稀疏矩阵求解避免内存膨胀。第二步处理分布式电源的短路电流贡献。传统配电网短路电流只来源于系统电源上一级变电站但接入分布式电源后故障点附近的DG也会向故障点注入短路电流。程序里按DG类型区分如果是经逆变器并网的DG按1.2倍额定电流计算贡献如果是直接并网的旋转电机型DG比如小型燃气轮机组按次暂态电抗Xd″计算。这个1.2倍系数在程序里是可配置的因为不同厂家的逆变器限流策略不完全一样查看设备参数后可以直接改常量值。第三步计算全电流和开关开断电流。程序除了输出短路电流周期分量初始值I″k还计算短路电流峰值ip。IEC 60909中峰值电流的计算要考虑衰减系数κ它跟系统的X/R比有关。控制系统阻抗参数已知通过κ 1.02 0.98e^(-3R/X)计算出来然后用ip κ × √2 × I″k得到峰值电流用于校核设备动稳定。3.4 核心代码逻辑与关键参数设置说了一堆原理我直接截取几个关键代码片段来讲解实际实现方式这样大家拿到程序框架后能更快上手改写成自己的版本。第一个是潮流计算的主迭代循环核心代码如下% 牛顿-拉夫逊法潮流主迭代 for iter 1:max_iter % 计算功率不平衡量 [dP, dQ] calculate_mismatch(V, theta, Ybus, P_spec, Q_spec); mismatch max(abs([dP; dQ])); % 判断是否收敛 if mismatch tolerance disp([迭代收敛迭代次数: , num2str(iter)]); break; end % 形成并求解雅可比矩阵修正方程 J form_jacobian(V, theta, Ybus); dX J \ [dP; dQ]; % 更新状态变量 theta theta dX(1:n_node-1); V V dX(n_node:end); end这段代码的逻辑很直白calculate_mismatch函数计算各节点功率不平衡量form_jacobian构造雅可比矩阵然后用左除求解线性方程组。这里有个Matlab的性能优化技巧用方阵左除J \ b而不是用inv(J)乘以b因为左除在数值稳定性上更好求解速度也更快。如果是更大规模的网络建议把J声明为稀疏矩阵求解效率能有数量级的提升。程序中的收敛精度tolerance默认设成1e-6实际工程应用已经足够如果做学术研究需要更高精度可改成1e-8但迭代次数会略有增加。第二个是分布式电源接入的处理逻辑% 根据DG控制模式设置注入功率 switch dg.controlMode case 1 % PQ恒功率模型 dgP dg.P; % 有功出力MW dgQ dg.Q; % 无功出力Mvar case 2 % PI恒功率因数模型 dgP dg.P; dgQ dg.P * tan(acos(dg.powerFactor)); case 3 % PV定电压模型 % 无功出力作为状态变量参与迭代初值设为零 dgP dg.P; dgQ 0; end % 更新节点注入功率DG按负负荷处理 P_spec(dg.node) P_spec(dg.node) - dgP; Q_spec(dg.node) Q_spec(dg.node) - dgQ;DG节点作为负的负荷处理是配电网潮流计算处理分布式电源最常用的做法。注意PV定电压模式下无功不能直接给定需要在迭代中不断修正最终使该节点电压稳定在目标值附近。程序里通过在每个迭代周期结束后检查电压偏差并按比例调整无功注入来逼近目标电压。第三个是短路电流计算的核心步骤% 节点k三相短路时的短路电流 Zk Zbus(k, k); % 节点k的自阻抗 I_sc V_prefault(k) / Zk; % 短路电流周期分量标幺值 % 折算到有名值 I_sc_kA I_sc * I_base_kA; % 考虑DG注入的附加短路电流 I_dg_contrib sum(dg_I_sc); % 汇总各DG贡献 I_sc_total I_sc_kA I_dg_contrib; % 计算短路电流峰值 kappa 1.02 0.98 * exp(-3 * R_X_ratio); I_peak kappa * sqrt(2) * I_sc_total;注意所有计算先在标幺值体系下完成最后统一折算到有名值这个顺序千万别搞反。如果混着用有名值和标幺值结果差了数量级你也未必能立即察觉但后面分析全废了。4. 仿真结果分析与工程解读4.1 分布式电源接入对节点电压的影响规律程序跑完后最直观的产出就是节点电压分布对比曲线。我以IEEE 33节点系统为基础在节点17主干线中段接入1.5MW分布式光伏记录接入前后的电压分布数据。未接入DG时由于线路阻抗上的压降越靠近馈线末端电压越低节点32的电压大约是0.94 p.u.左右末端电压偏低的问题很明显。接入1.5MW光伏后由于DG提供了本地有功支撑沿线电压普遍抬高末端电压回升到0.97 p.u.附近馈线的电压支撑效果非常显著。但当把DG容量继续加大到4MW时情况发生了本质变化馈线中后段节点电压开始超过1.0 p.u.部分节点甚至超过1.03 p.u.。这是因为DG出力超过了本地负荷需求剩余功率沿线返送潮流方向反转电压分布由“首端高末端低”变成“中间高两端低”。这个现象在实际工程中叫“电压抬高越限”是分布式电源高渗透率时最头疼的问题之一。程序通过这些对比曲线能直观告诉设计人员该馈线最多还能接纳多少分布式电源超过这个容量就得配套调压手段比如有载调压变压器、SVG等。工程上10kV配电网电压允许偏差范围一般在7%到-10%之间折算成标幺值就是0.9~1.07 p.u.程序里的判定阈值就是按这个标准来标注的。4.2 网络损耗的变化趋势与临界点分析除了电压网络损耗的变化也是接入方案评估的重点。我做了从0MW到5MW不同DG接入容量下的网损扫描输出的网损曲线呈现一个明显的“U形”特征网损先是随DG接入容量增加而下降在某个临界点场景二中大约在1.5MW附近网损降到最低继续增加DG容量网损反而开始回升甚至超过未接入时的水平。这个现象背后的物理机制很清晰分布式电源在负荷中心附近提供有功减少了线路上的长距离功率传输所以网损降低但当DG出力远超本地负荷时多余功率需要反送回首端变电站等于在原来的输电路径上多了一股反向流动的功率流线损自然就上去了。这就像一条本来自西向东单向送水的管道你在中间开了一个大流量泵站往里注水注水量小于下游用水量时上游水流量减小、管道损耗降低注水量超过下游用水量多余的水就要往回顶管道里的总流量反而变大损耗自然上扬。程序里对每个支路都有独立的损耗计算可以精确看出是哪些支路的损耗在DG容量增大后开始反弹。这个数据对优化DG并网位置非常有用如果两条候选馈线的网损临界点不同优先选择临界点更高的那条馈线意味着它能消纳更多DG而不至于让线损恶化。4.3 短路电流水平变化及对保护配置的影响短路电流计算这部分我设置了节点17接入不同容量DG后在接入点附近节点18发生三相短路的仿真条件。结果显示未接入DG时节点18的三相短路电流约为3.8kA接入2MW逆变器型DG后总短路电流增加到了约4.2kA增幅约10%。这个增量全部来自逆变器DG的贡献按限流系数1.2倍额定电流估算大约额外增加了0.4kA。虽然这个增幅看起来不大但工程上的敏感性不能按比例来评估——关键在于这个增量是否会让现有开关设备越过其额定开断容量。很多早期投运的10kV开关柜额定开断电流就是16kA或20kA如果系统短路电流本来就在15kA左右DG再补一刀可能就直接越限了。程序会自动将计算出的短路电流与设备额定开断能力比对并输出“校核结论”是合格还是越限越限会给出需要更换或限值DG接入容量的建议。另外还要关注保护配合的问题。DG接入后故障点的短路电流不再只由系统侧提供DG侧也有注入这意味着故障电流的幅值和方向都变了。原来按单端供电整定的过流保护可能灵敏度不足原来靠方向元件判别故障方位的逻辑也可能被打破。程序计算出的各支路故障电流分布表就是重新整定保护定值的基础输入数据。4.4 不同接入方案的对比评价方法程序最后会输出一个汇总对比表把不同场景下的关键指标放在同一张表里面方便做方案比选。表结构大致如下接入方案最低节点电压(p.u.)最高节点电压(p.u.)网损(kW)短路电流增量(kA)综合评价不接DG基准0.941.00202.50基准节点2接入1MW0.951.00186.50.19电压支撑有限节点17接入1MW0.961.01173.20.21降损效果好节点32接入1MW0.981.02168.90.22末端电压改善节点17接入4MW0.991.04232.80.82电压越限风险这个表在实际对接供电公司评审和编写接入方案时特别有用直接把仿真结论转换成可读性很强的工程语言。通过对比还能总结出一条实用经验中后段接入分布式电源对电压支撑和降损的综合效果最好但也是最容易在容量增大时触碰电压越限红线的位置。5. 常见问题排查与Matlab调试经验5.1 潮流计算不收敛怎么定位原因做配电网潮流计算遇到最多的问题就是迭代不收敛。程序报错之后不要慌按照我的排查经验90%的情况出在三个地方第一是功率基准值混乱。Matlab程序里变量多标幺值和有名值混用是最常见的坑。比如负荷数据用的kW/MW而DG参数里的有功用的又是标幺值算出来功率不平衡量怎么迭代都不会降下来甚至越迭代偏差越大。解决方法是在程序开头统一做单位转换全部转成标幺值之后不再出现有名值除非是最后输出结果时才转回有名值。我在init_network模块里专门加了单位检查断言发现节点电压初始值不在0.8~1.2 p.u.范围内就直接报错停止。第二是雅可比矩阵奇异。当网络中存在PV节点、平衡节点设置不当或者DG接入位置与负荷节点重合导致注入功率异常时修正方程可能是病态的。判断方法很简单把J矩阵的行列式打印出来如果接近0或者迭代中间出现NaN基本就是这个原因。解决办法是调整PV节点的无功初值或者把平衡节点的选择改成功率最充足的节点。第三是负荷数据与DG容量不匹配导致某个节点的净注入功率出现严重的不合理情况。比如节点17负荷只有300kW你给它配置了5MW的PV计算出来的节点电压翘到天际甚至出现负值迭代自然发散。工程上这种情况本来就不允许发生这是场景设计的问题回头检查你的DG容量方案。5.2 短路电流计算结果异常偏大或偏小另一个常见问题是短路电流计算结果跟预期差很远。我调试时遇到过一版程序算出来的短路电流是手算校验值的7倍多排查到最后发现是节点阻抗矩阵Z的形成方法出了问题——用导纳矩阵直接求逆时导纳矩阵里包含了变压器支路而IEC方法里应该把变压器阻抗折算到指定电压等级后再参与运算。这个教训提醒我短路电流计算前一定要先确认所有元件阻抗是否已经归算到同一电压等级。10kV计算时系统阻抗要按110kV/10kV变压器的短路电压百分比折算到10kV侧这个折算漏了结果必然偏小或偏大。程序里我把折算函数单独封装并在输出结果中附带阻抗归算明细方便复核。还有一个容易被忽略的细节就是三相短路时负序和零序网络对结果的影响。工程上配电网的零序阻抗受接地方式影响很大中性点不接地系统发生单相接地时接地电流很小而中性点经消弧线圈接地系统和小电阻接地系统短路电流特征又各不相同。程序里默认按三相短路计算如果要扩展单相接地故障分析需要添加序网络参数这部分我预留了接口但默认不启用以免不熟悉序分量理论的人误用出错。5.3 Matlab环境与脚本调试的实用技巧用Matlab跑这类程序环境配置和脚本调试也有一些实战经验可以分享。首先是Matlab版本问题代码中用到的矩阵左除、结构体数组、匿名函数这些特性从R2016b之后基本都稳定支持我测试过R2020a到R2023b的不同版本运行都没问题。老版本如果报错说函数或语法不支持建议直接把那些炫技写法改成最基础的for循环或function定义兼容性立刻大幅提升。调试过程中我习惯用断点和条件中断配合使用。在牛顿-拉夫逊迭代循环里设置条件断点——当mismatch连续两次迭代增大超过10倍时中断这时候基本可以断定迭代方向出错或初值给定不合理。另外建议在Matlab中开启“停止错误时进入调试模式”选项这样程序报错时会自动定位到出错的代码行配合各个工作区变量的检查排查效率能提高至少一倍。对于结果的输出我推荐用Matlab的导出函数把图片批量保存。设置好set(gcf, PaperPositionMode, auto)和print命令后生成的PNG图片分辨率高、白边少直接插入Word报告或者PPT里都不会变形失真。如果是在做论文需要的是矢量图那就在导出时选择SVG或PDF格式用Adobe Illustrator微调之后效果非常专业。还有一个小技巧就是程序运行时间控制。33节点系统的牛顿-拉夫逊潮流计算单次执行时间基本在0.5秒以内但如果你要做不同DG容量下的多批次扫描比如每0.5MW扫一次从0到5MW共11次总时间也不算长。真正费时间的其实是短路电流分析部分因为每种故障情况都要重新形成一次阻抗矩阵。建议把不随工况改变的初始化数据网络拓扑、导纳矩阵提前计算好并缓存循环执行时直接调用缓存不要每次重复计算。6. 程序扩展方向与个人实践经验总结最后说说这套程序的扩展潜力和我个人的实际操作体会。目前这个程序版本完成了潮流计算、短路电流分析和影响评估的核心功能但配电网分布式电源接入影响评估还有不少可以深挖的方向。如果后续有条件我计划在程序框架内加入动态稳定性分析通过构建分布式电源的动态模型在Simulink中观察并网瞬间的暂态响应——电压波动、频率变化和功率振荡。另外一个实用方向是继电保护定值校验功能把自动计算出的短路电流数据与保护装置的整定值进行自动比对快速筛选出受DG接入影响需要调整定值的保护装置。在实际做项目过程中我的体会是仿真程序的构建不要贪大求全。很多人一开始就想着要把光伏阵列的详细电气模型、MPPT控制算法、逆变器内部开关逻辑全都建出来结果搞了几个月连一个最基本的潮流都没跑通。正确的路径应该是先把最核心的基础分析流程打通——模型搭好、潮流跑通、短路计算校核完成、影响评估结果可信在这个骨架上再一层层增加复杂度。我这个程序最底层版本一共就两百多行代码跑出第一张电压曲线时兴奋了很久后面才逐步完善成现在这个多模块化的版本。还有一个建议给到正在做类似课题的同行仿真程序做出来的结果一定要想办法用理论公式手算或者用其他软件交叉验证一两个典型工况。我每次跑完一个新功能都会挑一个最简单场景用笔算校验确认数量级和趋势没问题后再继续下一步开发。这个习惯帮我拦住了不少低级错误比写一百行测试代码都管用。程序验证通过的工况和对应的手算校验过程记录建议保存在项目的README文档里后面论文评审或者项目验收时这些都是很好的佐证材料。这套程序现在已经跑过十来个实际10kV馈线的接入方案评估算出来的结论跟供电公司那边的正式计算报告吻合度很高。大家拿到程序框架后建议先用自己的本地网络参数替换IEEE 33节点数据别急着跑复杂场景——先把基线工况跑出来看电压分布趋势确认跟实际测量数据吻合了再做DG接入分析这个顺序千万别反了。