1. 为什么突然想做风光并网的仿真建模先说个背景。前阵子我接手了一个新能源消纳评估的活儿客户那边给的资料里有一句很扎眼的话某区域电网在风电出力高峰期PCC公共连接点电压偏差一度到了额定值的±8%谐波畸变率也逼近国标限值。当时我脑子里第一反应不是这电网真脆弱而是这数据是怎么量化出来的——因为光伏和风电的出力波动性、间歇性加上大量电力电子变流器接入电能质量问题是必然会发生的但难就难在把它建模、仿真、量化出来而不是停留在可能有影响这种定性结论上。所以这个项目其实是我自己给自己挖的一个坑搭建一套风光并网仿真模型用数据说话把并网对电压偏差、频率波动、谐波畸变、闪变这几个维度的影响全部跑一遍。用到的工具是MATLAB/Simulink配合Simscape Electrical以前叫SimPowerSystems模块库这是目前做电力系统仿真最主流的方案之一。整个模型分为三块风光发电单元、并网变流器及控制策略、电网侧等值模型。下文会把我从零开始搭模型的完整过程、参数整定的依据、实测中踩的坑全部摊开来讲适合正在做新能源并网方向毕业设计或工程预研的同行参考。有一点先讲清楚仿真建模不是把物理世界1:1搬进软件里而是在精度和计算成本之间找平衡。你不可能把每台风机的桨距角调节、每个光伏板的MPPT动态都建得极其精细那样仿真步长会被拖到微秒级一次暂态仿真跑几个小时毫无意义。真正有价值的做法是分时间尺度建模——电磁暂态、机电暂态、稳态潮流各取所需。我最终选定的建模策略是风光单元用等效模型风速模型双馈风机DFIG简化机电模型、光伏阵列工程数学模型Boost电路并网侧用两电平或三电平VSC变流器电网侧用一个含多馈线负荷的等值网络。这套组合既能反映风光出力波动对电网电压、频率的影响又不会让仿真计算量爆炸。下面把每个环节的细节展开。2. 风光发电单元建模稳态精度和动态响应怎么兼顾2.1 光伏阵列工程数学模型是性价比之王光伏仿真模型五花八门最精确的是基于物理机理的二极管五参数模型但这个模型对光照和温度的变化极其敏感在Simulink里如果每个光伏组件都用五参数模型加微分方程求解整个系统仿真速度会慢到怀疑人生。我用的方法是工程数学模型原理来自光伏电池的I-V特性方程I I_ph - I_0 [exp(q(V IR_s)/(nKT)) - 1] - (V IR_s)/R_sh这个方程本身不复杂但里面的I_ph光生电流、I_0反向饱和电流、R_s串联电阻、R_sh并联电阻、n二极管理想因子在物理模型里都需要精确参数。工程数学模型的做法是利用厂商提供的四个标准参数做近似短路电流I_sc、开路电压V_oc、最大功率点电流I_mp、最大功率点电压V_mp。通过这四个参数反推I_ph、I_0等中间量然后在给定光照S和温度T下修正得到不同工况下的I-V曲线。具体在Simulink里我建立了一个Mask封装的PV阵列子系统输入是光照强度W/m²和电池温度℃输出是直流电压电流。关键参数如下参数数值备注标准光照S_ref1000 W/m²标准测试条件标准温度T_ref25 ℃标准测试条件组件串联数N_s20每串组件数并联组串数N_p10组串并联数最大功率点电压V_mp35.2 V单板厂商规格书给定最大功率点电流I_mp9.94 A单板厂商规格书给定开路电压V_oc45.8 V单板厂商规格书给定短路电流I_sc10.58 A单板厂商规格书给定串并联之后整个光伏阵列的额定容量大约是20 × 10 × 35.2 × 9.94 ≈ 70 kW。这个容量规模是我故意选的——单个70kW的光伏单元接入380V低压配电网正好可以模拟分布式光伏入户的场景。如果需要做更大规模的光伏电站把N_s和N_p按比例放大就行工程数学模型的好处就是参数化缩放非常方便。这里有个关键点温度对光伏输出的影响比很多人想象中大得多。标准测试条件STC是25℃但实际运行时光伏板表面温度经常到50-60℃每升高1℃开路电压大约下降0.35%最大功率点功率大约下降0.45%。所以仿真时不能只设一个固定温度一定要加入温度变化输入否则对并网电压的影响评估会偏乐观。2.2 双馈风机DFIG电磁暂态模型可以简化到什么程度风电部分我选的是双馈异步发电机DFIG这是目前陆上风电的主流机型。DFIG的完整模型包括空气动力学模型桨距角控制、传动链模型两质量块、发电机电磁模型五阶或三阶微分方程、转子侧变流器RSC和网侧变流器GSC及其控制。如果全部建出来Simulink里光这一个风机模型就要几十个模块不太现实。我的折中方案空气动力学部分简化为风速-机械功率查表模型发电机部分用三阶暂态模型桨距角控制保留但变流器用平均模型Average Model替代开关器件模型。说下为什么能这么简化风速-机械功率查表模型风机的气动特性可以用功率系数C_p(λ, β)来描述λ是叶尖速比β是桨距角。C_p曲线可以直接从风机厂商的技术文档里获得经验公式。机械功率P_m 0.5ρπR²V³C_p代入风速V就能算出来。查表模型的好处是跳过复杂的空气动力学微分方程在仿真步长较大的情况下也能保证稳定性。三阶发电机模型DFIG的定子磁链动态比转子磁链动态快很多可以忽略定子磁链的暂态只保留转子磁链的三个状态变量d轴、q轴两个分量加上转速作为状态方程。这样模型从五阶降到三阶计算量大幅减少但转子侧动态保留完整足以反映风机对电网电压波动的响应。变流器平均模型把IGBT开关器件替换为受控电压源和受控电流源用控制信号的平均值驱动。这样仿真时间步长可以从微秒级放大到毫秒级甚至更大。缺点是丢失了开关频率附近的谐波特征——这意味着谐波分析只能针对电网侧其他非线性源比如光伏逆变器风机变流器自身的开关谐波需要额外单独做。你可能会问那DFIG的谐波影响岂不是被我人为丢掉了确实是但这是有意为之的取舍。DFIG变流器产生的谐波主要集中在开关频率附近通常2-10kHz而IEEE 519标准里关注的谐波频段也是到5kHz。如果用平均模型这部分谐波确实是盲区。我的处理办法是在变流器出口额外叠加一个基于开关频率的谐波电流源模型把谐波特征注入回去这样可以兼顾计算效率和谐波分析的完整性。风机的额定参数如下额定功率1.5MW额定电压690V额定风速11.5m/s切入风速3m/s切出风速25m/s。一台1.5MW风机配合前面70kW光伏并到一个10kV馈线上这种容量配比在真实的农村配电网里非常典型。2.3 并网接口为什么逆变器模型比风机电机模型更影响电能质量光伏和风电最终的并网点都是通过电力电子变流器接入电网的。对电能质量来说逆变器模型几乎决定了谐波和电压稳定性的仿真结果。所以在建模时变流器的拓扑和控制策略是重头戏。我用的拓扑是光伏侧Boost升压三相两电平VSC逆变器风电侧DFIG的转子侧和网侧各一个VSC。这里只说最常用的两电平VSC控制策略是经典的双闭环——外环电压/功率控制内环电流控制。以光伏并网逆变器为例外环直流母线电压控制在参考值比如700V输出d轴电流参考值i_d*q轴电流参考值i_q*设为0单位功率因数运行或者根据无功需求调节。内环PI控制器跟踪i_d和i_q参考值输出调制电压。调制策略SPWM或SVPWM。我用的SVPWM直流电压利用率比SPWM高约15%。在Simulink里搭建时控制环时间常数很关键。外环PI的时间常数通常取20-50ms内环PI取1-5ms。如果内外环时间常数靠得太近控制会互相干扰很容易出现发散。我在调参前先按内环比外环快10倍的原则粗调再根据阶跃响应的超调量和调节时间细调最后确定内环PIKp 2.5Ki 150外环PIKp 0.8Ki 20这套参数在70kW光伏并网仿真中表现良好直流母线电压超调小于5%调节时间约50ms。变流器模型我做了两套版本对比一套是详细的IGBT开关模型一套是平均模型。仿真结果显示对于电压偏差和频率波动的评估两套模型结果基本一致但对于谐波畸变率的评估详细模型能捕捉到5kHz以内的开关谐波平均模型则完全看不到。结论是如果你主要研究电压稳定性和频率响应平均模型足够如果研究谐波至少要保留一个变流器的详细开关模型。3. 电网侧等值建模与电能质量评价指标体系3.1 配电网等值从无穷大电源到多馈线负荷的取舍电网侧模型的精细程度直接决定了仿真结论的可靠性。常见做法是直接用Three-Phase Source模块模拟无穷大电源加一个短路阻抗这样最简单但对于电能质量研究来说极度不推荐。原因很简单无穷大电源的电压是刚性固定的无论你注入多少谐波或无功波动PCC电压几乎不变——你根本观察不到风光并网导致的电压偏差和谐波变化。我搭的是一个10kV多馈线等值配电网模型包含以下部分上级电网等值110kV/10kV变压器容量31.5MVA短路阻抗10.5%两条10kV馈线各带不同性质的负荷馈线1线性负荷为主居民照明、电机功率因数0.9馈线2含非线性负荷整流器、变频器用于对比风光并网是否加剧谐波风光并网点PCC设在馈线1的中段变压器和线路模型用PI等值电路等值阻抗值的计算是这一步的核心。线路和变压器的阻抗决定了PCC电压对功率波动的敏感度可以用短路比SCR来衡量SCR S_sc / S_new其中S_sc是PCC点的短路容量S_new是风光装机容量。我计算了一下S_cc 上级变压器容量 / 变压器阻抗 31.5MVA / 0.105 ≈ 300MVA风光总装机 1.5MW 70kW ≈ 1.57MWSCR ≈ 300 / 1.57 ≈ 191这个SCR偏高意味着并网影响会被弱化。为了模拟真实的弱电网场景这也是我认为最有工程意义的研究场景我把PCC的短路容量人为调低——在PCC点串联一个电抗器等效把SCR降到10左右。这个操作在仿真建模中有争议因为它改变了网络的物理结构。但在工程评估里接入点短路容量低导致的弱电网问题恰恰是分布式电源并网最常遇到的场景。SCR从191降到10之后同样的风光出力波动PCC电压波动幅度放大了近20倍。这个对比做完之后我才真正理解为什么业内总在强调弱电网下的新能源并网是需要专门研究的问题。3.2 扰动场景设置别只测稳态动态过程才是电能质量问题的重灾区仿真建模的核心价值在于能够人为制造各种极端场景这在真实电网中几乎不可能等到。我设计了以下几类扰动场景场景一光照骤降云层遮挡设置20s时光照从1000W/m²骤降到400W/m²持续5s后恢复。模拟光伏出力快速下降时PCC电压的跌落深度和恢复特性。场景二风速阶跃阵风设置30s时风速从8m/s阶跃到12m/s模拟风机出力突然增加对PCC电压的影响。12m/s接近额定风速风机出力接近额定值压升最明显。场景三三相短路故障在PCC点下游1km处设置一个三相短路故障持续0.1s后切除。这个场景主要看风光机组在故障期间和故障清除后的动态响应——是否脱网、恢复时间多长。场景四非线性负荷冲击在馈线2投入一个200kW的六脉波整流器负荷观察风光并网是否让谐波畸变率进一步恶化。每个场景仿真时长为60sSimulink仿真步长设为50微秒在风光单元出力变化较平缓的时间段可以放大步长但为了方便整体分析我全程用了固定步长。这带来的问题是仿真时间较长——60s仿真大概需要实际运行15分钟左右还可以接受。3.3 电能质量指标怎么量化不只THD还有这些电能质量评价不是感觉电压不稳就行需要精确到指标。我在Simulink里通过傅里叶分析模块和数学运算模块实时计算以下指标指标计算方法评价标准电压偏差(U_actual - U_nominal) / U_nominal × 100%GB/T 12325-2008±7%以内频率偏差f_actual - 50HzGB/T 15945-2008±0.2Hz以内电压总谐波畸变率THD√(ΣU_h²) / U_1 × 100%GB/T 14549-199310kV以下≤5%电压闪变Pst短期闪变严重度GB/T 12326-2008Pst≤1三相不平衡度负序分量/正序分量 × 100%GB/T 15543-2008≤2%短时≤4%实时计算这些指标有一个好处可以看到同一个扰动场景下不同指标的响应时序比如光照骤降时首先受影响的是电压偏差几百毫秒后谐波畸变率才跟着变化这种时序关系只有仿真能看得清楚。提示在Simulink里做FFT分析时采样窗口长度一定要和生产厂家电能质量分析仪的设置保持一致通常10个周波或12个周波否则你计算出的THD与实测值没有可比性。我这边的设置是10个周波0.2s50Hz对应IEC 61000-4-7标准。4. 仿真结果分析风光并网到底把电网污染成了什么样4.1 电压偏差弱电网下光伏的午间高电压现象先看最直观的电压偏差结果。在SCR191强电网条件下光伏满发70kW时PCC电压升高只有不到0.3%完全在标准范围内。但把SCR降到10之后同样的光伏出力导致PCC电压升高了约4.7%——虽然还在±7%限值内但已经超过了一半。这说明分布式光伏的电压影响和电网强度强相关。更有意思的是光照骤降场景光照从1000降到400W/m²后PCC电压从1.047pu跌到0.985pu电压跌落幅度达6.2%。这个跌落速度很快约0.5s内完成如果没有动态无功补偿这个电压波动足够让附近的敏感负荷感受到明显的闪变。风电部分的电压影响更复杂。风速阶跃从8m/s到12m/s时DFIG有功出力从约0.4MW跳到1.4MWPCC电压升高了约2.8%。但由于DFIG的变流器具备无功调节能力可以在电压升高时吸收无功容性运行把电压升高抑制到1.5%以内。这就是DFIG相对于光伏的天然优势——光伏逆变器在单位功率因数控制下没有任何无功调节裕度。4.2 谐波畸变率非线性负荷和逆变器开关谐波的叠加效应谐波分析是这次仿真中最有意思的部分。单独开光伏逆变器详细开关模型开关频率5kHz时PCC电压THD约为2.1%单独开非线性负荷时THD约为3.2%两者同时运行时THD达到了6.8%——超过了5%的标准限值。这个结果非常有工程意义风光并网和谐波源负荷之间的谐波叠加不是简单的算术加法。2.1%加3.2%等于5.3%但实测仿真结果是6.8%说明两者之间存在谐波放大效应。原因是光伏逆变器的输出阻抗和非线性负荷的谐波电流在PCC点产生了谐振放大特定次数的谐波主要是5次、7次被放大了。我在Simulink里用FFT分析工具看了各次谐波含量谐波次数光伏单独运行THD非线性负荷单独运行THD两者同时运行THD5次0.8%1.9%3.8%7次0.5%0.7%1.5%11次0.4%0.2%0.7%13次0.2%0.1%0.5%总THD2.1%3.2%6.8%5次谐波的放大比例最夸张从2.7%的代数和变成了3.8%。查了一下原因PCC点到变压器之间的线路阻抗和光伏逆变器的输出滤波电容在250Hz附近形成了并联谐振点导致5次谐波电流放大。解决思路有两个一是调整逆变器滤波器的参数把谐振点移到非特征谐波频率二是在PCC点加装无源滤波器5次单调谐滤波器。我在模型里加了5次滤波器后总THD降到了4.1%效果明显。4.3 频率波动风光占比多少才会威胁频率稳定频率稳定性是电力系统最重要的指标。仿真结果显示在我这套模型里风光出力波动引起的频率变化很小——最大偏差不到0.03Hz远低于0.2Hz的限值。原因很简单我的模型里电网仍然是一个以同步发电机为主的大电网风光占比只有约5%系统惯量足够大。但如果把风光占比提高到30%以上通过扩大风光容量实现情况就完全不同了。我做了个扩展仿真把风电容量提升到10MW、光伏提升到2MW总占比约30%左右频率最大偏差放大到0.12Hz接近限值的一半。更严重的是在风电出力骤降0.5s内损失5MW的情况下频率变化率达到0.15Hz/s这个变化率已经接近一些低频减载装置的启动阈值。这个结论印证了业内的担心高比例新能源接入后系统惯量下降频率稳定性成为比电压稳定性更棘手的问题。我的模型里还保留了同步发电机的惯性环节如果完全用电力电子变流器带孤岛负荷频率波动会大到不可接受——需要储能系统的虚拟惯量支撑。这是我下一阶段想做的扩展方向。4.4 暂态响应故障穿越能力决定脱网风险三相短路故障场景的结果最揪心。故障发生后光伏逆变器在故障开始后25ms检测到电压跌落保护动作60ms内脱网DFIG风机在故障期间由于Crowbar电路投入转子侧变流器被旁路机组靠定子直接并网运行故障清除后约200ms恢复正常励磁故障切除后PCC电压恢复过程中出现了约0.4s的振荡衰减这是风光机组和电网之间的动态交互从电能质量角度看光伏的快速脱网其实是双刃剑。一方面保护了逆变器本体另一方面脱网造成的光伏出力瞬间消失加剧了电压跌落——故障切除后电压恢复不完全长期处于0.95pu左右直到光伏重新并网后才回归1.0pu。现代并网标准如GB/T 19964-2016对光伏和风电的故障穿越能力提出了明确要求低压穿越期间光伏逆变器需要保持并网并向电网注入无功电流。我这套模型里用的还是传统保护策略下一步改进方向是加入低电压穿越LVRT控制策略让光伏在电压跌落到0.2pu时仍保持并网运行200ms。目前国内新投运的新能源电站都已经具备这个功能仿真模型也必须跟上这个标准。5. 建模过程中踩过的坑和实用技巧5.1 Simulink仿真速度的隐形杀手代数环和变步长我第一次搭完模型运行仿真时间长得无法忍受——60s仿真跑了快两个小时。排查了很久罪魁祸首是代数环。代数环的产生原因逆变器控制环节电流内环和主电路LCL滤波器之间存在瞬时耦合Simulink在计算时形成了无法直接求解的代数依赖被迫在每个仿真步长内进行迭代求解速度自然就慢了。解决办法在控制回路中增加Memory模块或Unit Delay模块打破代数环会引入一个步长的延迟对稳定性影响很小。在电流采样信号后加一个Rate Transition模块把控制频率和主电路仿真频率解耦。改用变步长求解器如ode23tb在动态变化剧烈时自动缩小步长平稳时放大步长整体速度提升5到10倍。改了之后仿真时间从2小时降到15分钟效果立竿见影。5.2 参数标幺化vs有名值别让小数值吃光仿真精度风电DFIG模型的转子侧参数和定子侧参数数量级差异巨大——定子电阻可能只有0.005pu转子漏抗只有0.1pu左右。Simulink默认用double精度计算但这些小数值在积分计算中很容易被截断误差吃掉导致模型出现高频数值振荡。我后来把所有电气参数统一改成标幺值pu体系发电机模型内部的电压、电流、阻抗都归一化到基准值上。这样所有参数数量级都在0.01到10之间数值稳定性大幅提升。同时Simulink的三相变压器、异步电机等内置模块默认都支持pu参数输入改动成本很低。5.3 风光波动场景别用阶跃函数硬模拟光照从1000W/m²骤降到400W/m²如果用Step模块直接跳变PCC电压波形会出现严重的振荡甚至数值发散因为实际光照变化是有一个过程的——云层遮挡是渐进的风速变化也是渐变的。正确做法是用Ramp模块加Saturation限幅把骤降过程模拟成0.5-1s的斜坡变化或者用Uniform Random Number模块叠加Sine Wave模拟更真实的光照波动。我甚至用真实的气象数据5分钟间隔的光照和风速测量值做插值后作为输入结果比任何理想化输入都真实。这个方法强烈推荐尤其是你做的是影响评估性质的分析输入的可信度直接决定结论的可信度。5.4 逆变器功率开关的仿真步长不要迷信越小越好逆变器开关模型详细IGBT模型对仿真步长非常敏感。步长太大2μs会导致开关动作瞬间的数值震荡甚至发散步长太小0.1μs则仿真时间爆炸。我的做法是分模块设置步长逆变器开关电路用0.5μs步长电网等值部分用50μs步长通过Simulink的Rate Transition模块连接不同步长的子模块。这样做的方式能兼顾逆变器谐波分析的精度和整体仿真的速度。5.5 结果验证仿真和理论计算的交叉验证仿真模型做出来之后一定要做验证否则就是个黑箱——你自己都不知道结果对不对。我的验证方法稳态潮流对比用MATPOWER或手算IEEE节点潮流对比仿真模型在稳态下的PCC电压和有功无功分布误差控制在1%以内。短路电流对比在PCC点设置三相短路短路电流的仿真值和理论计算值基于短路容量和等值阻抗对比误差小于3%。谐波传播对比用单一谐波源注入对比仿真测量的谐波电压传播和解析计算的阻抗频率扫描结果谐振点位置误差在5%以内。只有这三个层面都通过了我才有信心把仿真结果作为电能质量影响评估的依据。6. 从仿真到工程落地一些真实的边界和后续方向仿真做完了但我必须泼一盆冷水仿真的结论再漂亮也只能作为工程判断的参考不能替代真实测量。因为仿真模型里永远有简化——风速模型是统计模型而不是实际风场的湍流特征光照模型没有考虑光伏板表面的尘埃遮挡电网等值模型丢掉了下级馈线上几百个用户的随机负荷波动。这些都是仿真的固有边界。但反过来仿真的价值也是不可替代的你可以安全地在模型里做短路故障、做极端天气、做满负荷冲击这些在真实电网上几乎不可能做实验。我这次的仿真结果最终用在了两方面一是为客户提供光伏并网工程的可研报告数据支撑二是为后续的储能容量配置和电压治理方案设计提供基线数据。关于后续方向我想做以下扩展加入储能系统的虚拟同步机控制VSG改善高占比新能源下的频率稳定性把光伏逆变器的LVRT控制策略落实进去满足最新的并网标准要求用实测数据驱动模型校准让仿真结果和真实电能质量监测数据的误差进一步缩小。最后分享一个个人心得建模这个事儿最大的敌人不是数学而是你对系统和场景的理解深度。仿真模型本质上是你对物理系统认知的具象化——你对风光并网的理解有多深你的模型就能建得多准。不要在软件操作上花太多时间多花时间理解设备的控制逻辑、电网的运行特性这才是模型质量的真正决定因素。