开头最近在做配电网分布式电源规划这块的项目碰到一个挺典型的难题传统分布式电源优化配置模型里光伏电站大多被当成一个恒定功率因数的PQ节点来处理也就是说只考虑了有功出力无功这块基本按固定功率因数折算一下就算完事。但实际跑仿真的时候你会发现这样做出来的配置方案拿到现场往往不大好用特别是馈线末端电压波动、重载时段无功不足导致电压越限这些场景下固定的无功假设和真实运行状态偏差很大。原因是光伏电站现在的无功响应能力已经和以前完全不一样了。现代光伏逆变器一般都具备快速无功调节功能响应时间在几十毫秒到几百毫秒这个量级完全可以参与配电网的电压调节。如果优化配置模型里不考虑这个快速无功响应特性等于白白浪费了光伏电站本身自带的电压支撑能力配置出来的容量、选址很难说是最优的。所以这篇文章就围绕“计及光伏电站快速无功响应特性的分布式电源优化配置方法”展开记录一下我用Matlab实现这个模型的完整思路、数学推导、代码框架和踩坑经验。内容覆盖配电网潮流计算、光伏逆变器无功容量建模、无功电压灵敏度分析、粒子群优化算法求解这几个模块适合正在做分布式电源规划、微电网优化配置、电压治理相关课题的研究生和工程人员参考。代码思路可以直接迁移到IEEE 33节点等标准算例上也可以改吧改吧用到实际馈线上。1. 整体设计思路为什么要在配置阶段就考虑无功响应1.1 传统配置模型的盲区分布式电源优化配置本质上是回答三个问题装在哪里、装多大、什么运行策略。传统方法里光伏电站运行策略这一环往往被简化处理。最常见的做法是把光伏出力曲线简化为有功出力曲线无功按固定功率因数0.98甚至1.0折算代入潮流计算后验证节点电压是否越限、网损是否降低然后靠优化算法去找最优配置方案。这种简化在光伏渗透率不高、配电网裕度较大的年代问题不大。但现在分布式光伏渗透率越来越高一个台区装个几百千瓦甚至兆瓦级的光伏都很常见问题就来了光伏大发时段往往也是负荷较轻的时段有功功率就地消纳不掉就要倒送倒送功率在馈线上产生电压抬升严重的时候会导致末端电压超过1.07 pu。这种情况下光伏本身有没有无功调节能力直接决定了电压能不能压得住。固定功率因数模型默认光伏无功输出等于有功乘以某个系数完全没办法刻画光伏逆变器在电压越限时快速增发无功、压低电压的过程。做优化配置的时候就会倾向于把光伏容量压得很小甚至得出“这个台区不能再装光伏了”的结论——这其实不是光伏装多了而是模型没把光伏的无功潜力利用起来。1.2 快速无功响应特性到底是什么光伏逆变器的无功调节能力来自电力电子变换器的控制灵活性。现代并网逆变器在控制层面通常支持PQ控制、PV控制、VSG控制等多种模式其中和配电网电压治理最直接相关的是无功功率控制Q控制和无功-电压下垂控制Q-V droop control。所谓“快速无功响应”有两个关键量化指标。第一是响应速度逆变器电流内环的时间常数一般在几十毫秒量级无功功率从零增长到额定值通常只需要100到300毫秒第二是调节范围只要逆变器不过热、不过调制在视在功率容量范围内可以灵活分配有功和无功。一般光伏逆变器的无功容量边界近似为Q_max sqrt(S_rated^2 - P^2)其中S_rated是逆变器额定视在功率P是当前有功出力。这个关系非常重要它说明一个关键事实光伏在满发状态下无功调节能力最弱在欠发状态下无功调节能力最强。而恰恰是光伏满发的午间时段配电网最容易出现高电压这种“越需要无功的时候无功容量反而越小”的张力就是优化配置模型必须显式刻画的核心矛盾。如果配置阶段不把这个无功容量随有功变化的特性建模进去等到实际运行的时候才发现并联的无功补偿装置不够用那就只能额外投入SVG、电容器组成本高不说响应速度还比不上逆变器的快速无功调节。1.3 把无功响应纳入配置模型的可行路径在配置阶段考虑光伏无功响应其实不需要做动态仿真那是另外一个时间尺度的问题。我的做法是把无功响应特性映射到静态潮流模型里核心思路是在潮流计算中把光伏电站建模为具备可变无功输出能力的节点其无功输出范围由逆变器容量边界给定具体的无功输出值由电压偏差和优化目标联合决定。具体来说配置模型采用双层迭代框架。外层是容量和选址优化内层是考虑光伏无功调节后的配电网潮流计算。内层潮流计算时光伏节点根据本节点的电压偏差调整无功输出无功调节量受逆变器无功容量边界约束。这样潮流收敛后得到的节点电压、网损就是真实考虑了光伏快速无功响应后的结果外层优化算法基于这个结果去评价配置方案的优劣。这个方案的好处是模型复杂度可控计算量在可接受范围内同时能抓住光伏无功响应配置方案影响的主要矛盾。相比动态仿真它不需要搭建逆变器控制模型和电磁暂态仿真相比固定功率因数模型它又能真实反映无功调节能力对电压的支撑作用。属于性价比很高的一种建模思路。2. 数学模型构建目标函数、约束条件与光伏无功建模2.1 目标函数设计分布式电源优化配置的目标函数需要综合考虑经济性和技术性不能只看一个指标。我的模型采用多目标加权的方式三个子目标分别是年度综合费用最小、系统有功网损最小、电压偏差最小。年度综合费用包括分布式电源的年投资等年值费用和运行维护费用。等年值计算的处理方法是用资金回收系数把初始投资摊到每一年计算公式为C_inv C_cap * P_DG * r * (1r)^n / ((1r)^n - 1)其中C_cap是单位容量投资成本元/kWP_DG是安装容量r是贴现率n是项目寿命。运维费用按安装容量和年发电量的一定比例估算。系统有功网损通过潮流计算得到这没什么好说的直接对每条支路的I^2R求和。电压偏差指标用的是所有节点电压偏差的平方和这个指标能比较敏锐地反映电压越限和电压质量恶化程度。之所以不用最大电压偏差而用平方和是因为平方和能体现整体电压水平避免单个节点电压被平均掉。三个目标量纲不一致直接加权求和没有意义。我的处理办法是采用标幺化归一每个子目标先除以各自在“不装分布式电源”基准情况下的值然后再加权。这样三个子目标的数值都在1.0附近波动权重设置就相对直观。2.2 潮流计算中的约束条件潮流计算要满足的约束包括节点功率平衡约束、节点电压上下限约束、支路电流上限约束、分布式电源总装机容量约束。节点功率平衡约束就是经典的有功和无功功率平衡方程对于每个节点iP_i - P_load_i P_DG_i - P_loss_i? 不对更准确的说法是注入功率等于负荷功率减去分布式电源出力然后根据节点导纳矩阵计算注入电流和功率。无功环节同样处理只不过光伏节点在潮流迭代过程中无功是变量受无功容量边界约束。电压约束方面配电网节点的电压一般要求控制在0.93到1.07 pu范围内。分布式电源接入后尤其是集中接入末端节点时电压抬升是一个很现实的问题。模型里需要显式校验每个节点的电压幅值如果不满足约束就要通过优化算法去调整配置方案。还有一个容易被忽略的约束是分布式电源总容量不能超过变电站主变容量的某个比例比如25%或30%。这个约束的目的是防止分布式电源渗透率过高带来的调峰和反送问题在工程实际中基本是硬性要求。2.3 光伏电站快速无功响应建模这是整个模型的核心部分。我把光伏电站的无功输出建模为电压偏差的分段函数并考虑逆变器容量边界限制。首先确定光伏节点的无功调节能力。设光伏当前有功出力为P_pv逆变器额定视在功率为S_rated那么无功调节能力的上限和下限为Q_max sqrt(S_rated^2 - P_pv^2) Q_min -sqrt(S_rated^2 - P_pv^2)注意这里无功的正负号约定吸收无功为正容性还是发出无功为正感性不同文献定义可能不一样自己代码里统一就行关键是边界公式别搞错。然后定义无功-电压响应关系。一种常用的建模方式是线性下垂控制电压偏差大于死区时光伏按比例增发无功Q_pv K_q * (V_i - V_ref)其中V_ref是参考电压一般取1.0 puK_q是下垂系数V_i是光伏接入节点的实际电压。计算得到的Q_pv还需要经过容量边界约束if Q_pv Q_max Q_pv Q_max; elseif Q_pv Q_min Q_pv Q_min; end但这里有一个迭代收敛的问题。无功输出改变后节点电压会变电压变了无功又该变所以必须做迭代。我的做法是在潮流计算内部每迭代一次就根据最新电压更新光伏无功直到相邻两次迭代的无功差值小于阈值。对于大多数配电网算例这种迭代方式十几步就能收敛不会显著拖慢优化速度。有意思的是下垂系数K_q的取值对结果影响很大。K_q取值太大电压稍微偏一点无功就顶满了实际上等于把光伏变成了恒无功源会削弱电压调节的平滑性K_q取值太小无功响应太迟钝电压越限了都拉不回来。我的经验是K_q取额定容量的1.0到2.0倍之间比较合理这时候电压偏差0.05 pu就能触发50%左右的无功输出既不会太激进也不会太迟钝。3. Matlab代码实现总体架构与关键模块拆解3.1 代码整体结构整个Matlab工程我分了五个核心文件每个文件职责单一方便调试和复用。第一个是主程序文件负责参数初始化、调用优化算法、汇总结果和输出图表。第二个是配电网数据文件包括IEEE 33节点系统的支路参数、节点负荷、基准电压和基准功率。第三个是光伏无功响应计算函数输入当前节点电压、光伏有功出力和逆变器参数输出无功输出值。第四个是潮流计算函数采用前推回代法或牛顿-拉夫逊法内部调用光伏无功响应函数做迭代。第五个是粒子群优化算法函数负责容量选址寻优。这种模块化设计最大的好处是方便替换。比如你今天用IEEE 33节点明天想换IEEE 69节点只需要改配电网数据文件就行其他模块逻辑完全不用动。再比如你不想用粒子群想换成遗传算法或者灰狼算法只需要把优化算法函数换掉目标函数和潮流计算都能复用。3.2 光伏无功响应函数的实现细节无功响应函数虽然逻辑很简单但有几个细节决定成败。先看核心代码。function Q_pv pvReactiveResponse(V_pu, P_pv_pu, S_rated_pu, K_q, V_ref, V_dead) % 光伏电站快速无功响应模型 % 输入: % V_pu - 光伏接入节点当前电压 (pu) % P_pv_pu - 光伏当前有功出力 (pu, 基准容量为系统基准功率) % S_rated_pu - 逆变器额定视在功率 (pu) % K_q - 无功-电压下垂系数 % V_ref - 参考电压 (pu) % V_dead - 电压死区 (pu) % 输出: % Q_pv - 光伏无功输出 (pu, 感性为正) % Step 1: 计算逆变器可用无功容量边界 Q_cap sqrt(max(S_rated_pu^2 - P_pv_pu^2, 0)); % Step 2: 计算电压偏差, 引入死区避免频繁动作 dV V_pu - V_ref; if abs(dV) V_dead Q_pv 0; elseif dV 0 Q_pv K_q * (dV - V_dead); else Q_pv K_q * (dV V_dead); end % Step 3: 无功容量边界约束 Q_pv max(min(Q_pv, Q_cap), -Q_cap); end有几个细节要特别说明。首先sqrt(max(S_rated^2 - P_pv^2, 0))处理了一个边界情况当P_pv无限接近S_rated时括号内可能因为浮点数误差出现微小负数不加max保护直接开根号会得到复数运行直接报错。这个坑我踩过一次排查半天才发现是浮点数精度问题。其次是死区的设计。实际工程中逆变器不可能电压稍微波动一点就频繁调节无功否则功率器件寿命受影响。所以我在响应模型里加了死区电压偏差在死区范围内时无功为零。这个设计在优化配置阶段的意义是避免模型过度乐观——如果假设光伏在所有场景下都完美调节无功配置出来的方案可能过于激进。另外无功输出采用“感性为正”的约定这在Matlab的潮流计算里需要和负荷的无功约定保持一致。我习惯把负荷的无功用正数表示感性负荷吸收无功光伏吸收无功时也是正数光伏发出无功时是负数。这个约定如果不统一潮流计算结果会完全错乱。3.3 潮流计算模块的迭代逻辑配电网潮流计算我用的是前推回代法这种方法对辐射状配电网收敛性好编程简单计算速度快非常适合嵌入到优化算法里反复调用。和牛顿-拉夫逊法相比前推回代不需要求雅可比矩阵代码量少很多而且对配电网这种R/X比值较大的网络反而更稳定。含光伏无功响应的潮流计算核心逻辑是先忽略无功响应做一次基础潮流然后根据电压偏差更新光伏无功再重新算潮流如此迭代。用伪代码描述如下% 初始化光伏无功 Q_pv zeros(n_pv, 1); for iter 1:max_iter % 前推回代潮流计算 [V_mag, V_angle, P_branch, Q_branch] backwardForwardSweep(...); % 更新光伏无功出力 Q_pv_new zeros(n_pv, 1); for k 1:n_pv Q_pv_new(k) pvReactiveResponse(V_mag(idx_pv(k)), ...); end % 检查收敛 if max(abs(Q_pv_new - Q_pv)) tol break; end Q_pv Q_pv_new; end收敛判据用的是相邻两次迭代的光伏无功变化量。这个收敛判据比用电压变化量更直接因为整个迭代的目的就是让光伏无功和节点电压自洽。收敛精度我设的是0.0001 pu实际运行中通常5到15次内迭代就能收敛。3.4 粒子群算法求解配置问题外层优化我选粒子群算法PSO主要原因是它参数少、实现简单、对非线性约束的适配性好。PSO的决策变量设计是这样的如果有N个候选安装节点每个节点的安装容量范围是0到上限值那么决策变量就是N维向量每个维度代表对应节点的安装容量。这样选址和定容一次性解决容量为0就代表该节点不安装。% PSO参数 n_particles 100; n_iter 200; w 0.9; % 惯性权重 c1 1.5; % 个体学习因子 c2 1.5; % 社会学习因子 % 初始化种群 particles rand(n_particles, n_vars) .* upper_bound; velocities zeros(n_particles, n_vars);PSO速度更新公式是标准的velocities w * velocities c1 * rand * (pbest - particles) c2 * rand * (gbest - particles); particles particles velocities;需要注意的几个问题。首先是粒子速度的限幅如果速度过大粒子可能跳出搜索空间产生不合理的配置方案。我一般把速度上限设为决策变量范围的10%到20%。其次是约束处理如果某个粒子的决策变量超过上限直接截断到上限如果总容量超过约束按比例缩减各节点容量。目标函数计算是性能瓶颈每次粒子位置更新都要重新计算潮流。IEEE 33节点系统配20个粒子迭代100次总共需要2000次潮流计算每次潮流计算里含光伏无功迭代约10次整体算下来几秒钟就能完成。但如果换到IEEE 123节点这种大系统计算时间会显著上升这时候可以考虑用雅可比矩阵稀疏化或者并行计算工具箱加速。4. 算例分析与对比到底带来了多少改善4.1 测试系统与场景设置我用IEEE 33节点标准测试系统做验证。这个系统有33个节点、32条支路基准电压12.66 kV基准功率10 MVA系统总负荷约3715 kW和2300 kvar。按照惯例我把节点0设为平衡节点变电站母线其余节点都是PQ节点。光伏候选安装节点选了5个分别是节点8、节点12、节点18、节点24、节点30。这几个节点的位置有代表性节点18是馈线末端电压灵敏度最高适合观察光伏无功对末端电压的支撑效果节点8和节点12在馈线中段节点24和节点30在另一条分支上。每个节点的光伏安装容量上限设为500 kW。光伏出力曲线我简化处理取四个典型场景光伏大发出力为额定容量的90%对应夏季正午、光伏中等出力60%对应春秋季上午、光伏低出力20%对应早晚、光伏出力为零夜间。负荷也按典型日负荷曲线取对应时段的数值。这样四个场景基本能覆盖光伏运行的主要区间。4.2 两种建模方式的对比结果为了验证考虑快速无功响应的效果我设置了两个对比方案。方案一用传统固定功率因数模型功率因数取0.98方案二用本文的快速无功响应模型下垂系数取1.5死区取0.01 pu。两个方案都用同样的PSO参数跑优化保证对比的公平性。优化结果差异非常明显。方案一的优化结果是总共安装约1350 kW光伏方案二总共安装约1580 kW光伏。方案二的装机容量比方案一多出17%。这个差异的核心原因是方案一要保持固定0.98的功率因数光伏逆变器始终留出一定无功容量不参与有功发电而且这个预留比例是固定的系统不需要时就浪费了方案二的光伏逆变器在电压正常时不发无功把全部容量用于有功输出电压需要时才动态增发无功所以同一个节点能装下更大的光伏容量。网损方面方案一优化后的系统网损是约85 kW方案二约72 kW降幅15%左右。电压质量方面方案二的最大电压偏差为0.041 pu方案一为0.052 pu。这些指标都说明考虑无功响应后光伏配置方案在技术和经济指标上都有提升。4.3 下垂系数灵敏度分析我还做了一组下垂系数的灵敏度分析这个结果对实际工程有直接参考意义。固定光伏配置不变把K_q从0.5逐步调到3.0观察系统最大电压偏差和总网损的变化。结果显示K_q从0.5增大到1.5时最大电压偏差明显下降从0.062 pu降到0.041 puK_q继续从1.5增大到3.0电压改善效果变得非常平缓基本稳定在0.040到0.042 pu之间。这说明下垂系数存在一个经济合理的区间超过这个区间后继续加大投入意义不大。这也解释了为什么在建模过程中K_q不能随便设一个值就完事。我的建议是在做正式优化前先跑一遍灵敏度分析摸清K_q在什么范围内对结果影响显著然后在这个范围内选择。如果你完全依赖优化算法去自动寻优K_q和配置容量会大幅增加决策变量维度计算时间翻好几倍而且容易陷入局部最优。5. 常见问题与调试经验速查5.1 潮流计算不收敛这是最常见的问题。光伏无功响应模型引入了迭代环节如果无功更新逻辑和潮流算法配合不好很容易出现振荡不收敛。第一种情况是下垂系数太大导致振荡。电压偏差大时无功输出也大无功改变电压相位和幅值可能反过来让电压偏差更大形成正反馈振荡。解决办法是减小K_q或者无功更新时加入阻尼因子比如新Q等于旧Q乘以0.5加计算值乘以0.5做一阶低通滤波。第二种情况是光伏无功初值设置不合理。如果初值设置得离真实解太远前几次迭代可能步子迈太大。我习惯把光伏无功初值设为0让迭代逐步逼近。第三种情况是系统本身处于临界状态比如负荷特别重、电压已经低到接近下限。这时候光伏无功已经顶满了但电压还是越限潮流就会在越限点附近反复跳。解决办法是检查是不是某些场景下需要增加无功补偿装置或者提升变电站母线电压。5.2 IEEE 33节点系统数据格式容易出错IEEE 33节点系统的经典数据格式是每个节点一行内容包括支路编号、首端节点、末端节点、支路电阻和电抗。很多人下载的数据文件里首端节点编号是从0开始的0代表变电站母线但Matlab数组索引是从1开始的两者差1很容易搞错。我的做法是读取数据后统一加1把节点编号映射到Matlab索引同时单独存一个原始编号数组用于结果展示。还有一个坑是IEEE 33节点的单位——经典数据里阻抗单位是欧姆功率单位是kW和kvar潮流计算时如果不做标幺化换算结果会非常离谱。必须先把所有电气量换算到标幺制下再计算。5.3 粒子群算法的随机性问题PSO是随机优化算法每次运行结果可能不一样。如果代码里不固定随机种子同一个算例跑十次可能得到十个不同的最优解。这不是代码错了是算法本身的特性。在科研或工程报告中需要稳定结果时我通常在程序开头固定随机种子rng(42);这样每次运行的结果都能复现。但要提醒的是固定随机种子意味着只有一个局部最优解不一定代表全局最优。建议的做法是多跑几次不同的随机种子对比结果后再下结论。另外一个PSO的坑是收敛速度不均。前期粒子分散探索全局搜索能力强后期粒子聚集局部搜索能力强。但参数设置不当可能导致前期收敛太快过早陷入局部最优或者后期收敛太慢计算资源浪费。我的经验是惯性权重w从0.9线性递减到0.4配合速度限幅大多数配电网优化问题几十代就能收敛到稳定值。5.4 代码逻辑与算例复现我们现在来把上面拆解的各个模块串起来用一段可直接运行的主程序骨架展示整体调用逻辑方便你对照自己的代码做修改。%% 主程序: 计及光伏快速无功响应的分布式电源优化配置 clear; clc; close all; % 固定随机种子, 保证结果可复现 rng(42); %% 1. 加载IEEE 33节点系统数据 [bus, branch] loadIEEE33bus(); baseMVA 10; % 基准功率 MVA baseKV 12.66; % 基准电压 kV Zbase baseKV^2 / baseMVA; % 阻抗基准值 Ohm %% 2. 设置光伏参数 pvNodes [9, 13, 19, 25, 31]; % 注意: 原始编号8,12,18,24,30加1 pvNodes pvNodes - 1; % 如果bus矩阵索引从0开始则调整 n_pv length(pvNodes); S_rated_max 500; % 单节点最大容量 kW S_rated_vec zeros(n_pv, 1); % 加载光伏出力曲线 P_pv_curve [0.9, 0.6, 0.2, 0.0]; load_curve [0.8, 1.0, 0.7, 0.4]; %% 3. 粒子群优化参数 n_vars n_pv; n_particles 50; n_iter 100; ub repmat(S_rated_max / baseMVA, 1, n_vars); % 标幺化 lb zeros(1, n_vars); %% 4. 初始化粒子种群 particles repmat(lb, n_particles, 1) rand(n_particles, n_vars) .* (ub - lb); velocities zeros(n_particles, n_vars); pbest particles; pbest_fit inf(n_particles, 1); gbest particles(1, :); gbest_fit inf; %% 5. 粒子群主循环 for iter 1:n_iter w 0.9 - 0.5 * (iter / n_iter); % 线性递减惯性权重 for i 1:n_particles % 计算目标函数值 fitness objectiveFunction(particles(i, :), bus, branch, ... pvNodes, P_pv_curve, load_curve, baseMVA); % 更新个体最优 if fitness pbest_fit(i) pbest_fit(i) fitness; pbest(i, :) particles(i, :); end % 更新全局最优 if fitness gbest_fit gbest_fit fitness; gbest particles(i, :); end end % 更新速度和位置 for i 1:n_particles velocities(i, :) w * velocities(i, :) ... c1 * rand(1, n_vars) .* (pbest(i, :) - particles(i, :)) ... c2 * rand(1, n_vars) .* (gbest - particles(i, :)); % 速度限幅 v_max 0.2 * (ub - lb); velocities(i, :) max(min(velocities(i, :), v_max), -v_max); % 位置更新与边界约束 particles(i, :) particles(i, :) velocities(i, :); particles(i, :) max(min(particles(i, :), ub), lb); end fprintf(迭代 %d/%d, 最优适应度: %.4f\n, iter, n_iter, gbest_fit); end %% 6. 结果展示 figure; bar(gbest * baseMVA); xlabel(候选安装节点); ylabel(配置容量 (kW)); title(分布式电源优化配置结果); grid on;目标函数里最核心的部分是把每个场景的潮流算一遍然后把网损、电压偏差和投资成本加权求和。每个场景的光伏有功出力需要按光伏出力的标幺值和配置容量相乘得到然后对每个光伏节点调用pvReactiveResponse更新无功。5.5 从33节点到更大规模系统的迁移这套代码迁移到IEEE 69节点、123节点甚至实际馈线时有几个地方需要特别注意。首先是阻抗数据格式差异。IEEE 69节点的数据格式和33节点类似但123节点系统是三相不平衡系统前推回代法要扩展成三相前推回代光伏无功响应策略也要从单相扩展到三相独立控制。这已经不是小改动了需要重新写潮流模块。其次是计算效率问题。实际馈线可能有几百个节点候选安装节点如果也设置几十个PSO的决策变量维度就很高计算量呈指数级增长。我的建议是做两阶段优化第一阶段用连续变量在全部候选节点上粗优化第二阶段在结果基础上选几个关键节点做精细优化。这个思路类似多分辨率搜索能大幅减少计算时间。再有一个是数据校验。实际系统的负荷数据可能有错误或者缺失数据质量直接决定优化结果靠不靠谱。我在处理实际馈线项目时第一步永远是用没有分布式电源的基态潮流校验系统模型看算出来的网损和电压分布和实际量测对不对得上对不上就先调数据不急着跑优化。5.6 光伏无功响应模型的适用边界最后把话说明白一点这个模型不是万能的它有自己的适用边界。我在建模时做了简化假设光伏逆变器在稳态潮流层面表现为电压-无功下垂特性没有考虑动态响应过程和时间常数。如果研究对象是暂态电压稳定性、次同步振荡这类问题动态模型是必须的这个静态模型只能做配置阶段的初步评估。另外配电网里已经有其他无功补偿设备的情况下光伏无功响应策略需要和它们做协调。比如配电变压器自带的有载调压分接头OLTC动作速度慢但调节范围大光伏逆变器响应速度快但调节范围有限两者协调好了能做到全局优化。这个协调策略在配置模型里可以考虑为联合优化但模型复杂度会上一个台阶等以后有专门的需求再详细展开。在现阶段这个模型解决的核心问题就是在规划层面回答清楚“光伏无功响应能力到底值多少钱、省多少投资”这个关键问题。实测下来把光伏的无功特性建进去之后同样的电网骨架能多接纳不少分布式光伏网损和电压指标还更漂亮。对于现在各地光伏渗透率快速提升的局面这个建模思路有明确的工程参考价值。
