CN108649560A - 高渗透率分布式光伏发电集群实时仿真建模方法 - Google Patents
高渗透率分布式光伏发电集群实时仿真建模方法 Download PDFInfo
- Publication number
- CN108649560A CN108649560A CN201810392773.7A CN201810392773A CN108649560A CN 108649560 A CN108649560 A CN 108649560A CN 201810392773 A CN201810392773 A CN 201810392773A CN 108649560 A CN108649560 A CN 108649560A
- Authority
- CN
- China
- Prior art keywords
- model
- photovoltaic
- real
- load
- time simulation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000010248 power generation Methods 0.000 title claims abstract description 26
- 230000035699 permeability Effects 0.000 title claims abstract description 22
- 238000005094 computer simulation Methods 0.000 title claims abstract description 17
- 238000004088 simulation Methods 0.000 claims abstract description 79
- 238000004146 energy storage Methods 0.000 claims abstract description 52
- 230000001360 synchronised effect Effects 0.000 claims abstract description 17
- 239000003990 capacitor Substances 0.000 claims description 17
- 238000010521 absorption reaction Methods 0.000 claims description 12
- 230000005284 excitation Effects 0.000 claims description 11
- 230000009466 transformation Effects 0.000 claims description 8
- 230000005611 electricity Effects 0.000 claims description 5
- 238000009825 accumulation Methods 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 2
- 240000002853 Nelumbo nucifera Species 0.000 abstract 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 abstract 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000029087 digestion Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003064 k means clustering Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000006116 polymerization reaction Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- H02J3/383—
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
Abstract
本发明提供了一种高渗透率分布式光伏发电集群实时仿真建模方法,该方法在RT‑LAB实时仿真平台上构建包含源侧子系统、网侧子系统、负荷子系统和储能子系统模型进行仿真。方法包括:在RT‑LAB仿真平台中搭建了基于FPGA小步长仿真的逆变器详细模型;使用受控电压源模型代替IGBT电压源型换流器搭建了逆变器等效简化模型;搭建能够反映电网频率特性的同步发电机组模型、反映综合负荷特性的ZIP负荷模型、反映系统稳定性的储能电池模型。该方法不仅考虑了电网频率波动和源、网、荷、储四个子系统之间的协调性,并在仿真规模大的情况下确保了仿真精确度、降低了模型的复杂度和仿真计算时间。
Description
技术领域
本发明涉及高渗透率分布式光伏发电集群实时仿真建模方法,属于分布式能源并网建模与仿真的技术领域。
背景技术
光伏发电作为一种典型的分布式电源,已经成为人类利用可再生能源发电的一个重要途径。因此,对于光伏系统发电过程的仿真建模研究具有重大意义。但是电力系统频率的稳定性要求对可再生能源发电有着较高的要求,太阳能发电具有间歇性、波动性以及不可调度性等特点,因此大规模分布式光伏发电集群会存在并网难、电网消纳能力弱的问题。同时,电网渗透率逐渐增大使得大量电力电子变换器引入到电力系统中,从而大量的非线性负载也加入到了供电系统中,对电力系统造成严重的污染,对电网的规划设计、调度运行等方面带来了新的问题,因此对高渗透率分布式光伏发电系统进行仿真建模亟待研究。
目前,已有多篇国内外学术论文进行分析并提出解决方案,例如:
文献1《Remon D,Cantarellas A M,Rodriguez P.Equivalent Model of LargeScale Synchronous Photovoltaic Power Plants[J].IEEE Transactions,2016,52(6):5029-5040.》(《大型同步光伏电站的等效模型》——2016年IEEE电力电子期刊);
文献2《Ma Z,Zheng J,Zhu S.Online Clustering Modeling of Large-scalePhotovoltaic Power Plants[C]//IEEE Power&Energy Society General Meeting.IEEE,2015:1-5》(《大型光伏电站的在线聚类建模》——2015年IEEE PESGM会议论文集);
文献3《顾伟,李培鑫,盛万兴.一种分布式光伏电站集群的动态等值建模方法:,CN107565601A[P].2018.》。
文献1通过使用同步功率控制器控制等效光伏发电单元来获得大型同步光伏电站的等效模型,但其针对的是具有虚拟同步特性的光伏电站,对普通光伏电站不适用。文献2通过使用逆变器控制参数的特征距离作为聚类指标并利用K-means聚类算法,提出了一种大型光伏电站在线聚类建模仿真方法,但真实逆变器控制参数不易获取,并且此方法使光伏电站模型不断更新对等效造成一定难度。文献3选择能够准确描述光伏电站动态特性的电气量波形作为聚类指标,再根据聚类结果计算同簇光伏电站的聚合等值参数并进行配电网的等值化简。
综上所述,上述文献存在以下的不足:
1)没有考虑分布式光伏发电系统对电力系统频率的影响,没有对可以反映频率波动的电网环境进行模拟。
2)光伏发电集群的等效模型会随着动态聚类结果的不同而更新,对建模和仿真造成了一定的难度。
3)对大规模的电站集群进行仿真时,Simulink/PSCAD等仿真软件不能同时满足对仿真规模和仿真精确度的要求。
发明内容
本发明的目的是要解决高渗透率分布式光伏发电集群仿真规模大复杂度高的问题。针对现有技术的不足,本发明提供了一种高渗透率分布式光伏发电集群实时仿真建模方法。所采用的技术方案为:包括在RT-LAB实时仿真平台中搭建基于FPGA小步长仿真的逆变器详细模型;使用受控电压源模型代替IGBT电压源型换流器搭建逆变器等效简化模型;搭建能够反映电网频率特性的同步发电机组模型、反映综合负荷特性的ZIP负荷模型、反映系统稳定性的储能电池模型;并按照仿真要求将各子系统组合建模。
本发明的目的是这样实现的,本发明提供了一种高渗透率分布式光伏发电集群实时仿真建模方法,其特征在于,所述实时仿真建模方法包括基于RT-LAB实时仿真平台构建源侧实时仿真子系统、网侧实时仿真子系统、负荷实时仿真子系统和储能实时仿真子系统;其中,源侧实时仿真子系统由光伏并网逆变器模型组成,所述的光伏并网逆变器模型包括光伏并网逆变器详细模型和光伏并网逆变器等效简化模型;网侧实时仿真子系统由同步发电机组模型组成;负荷实时仿真子系统由ZIP负荷模型组成;
所述的高渗透率分布式光伏发电集群包括3个光伏发电区域,分别记为区域1、区域2和区域3,区域1、区域2、区域3和负荷实时仿真子系统彼此并联到公共母线后与网侧实时仿真子系统串联;
所述区域1、区域2和区域3中,每个区域均包括7个光伏电站和1个储能实时仿真子系统,7个光伏电站和1个储能实时仿真子系统依次串联;
所述光伏电站中均包含两台并联的光伏并网逆变器,区域1的光伏并网逆变器均使用光伏并网逆变器详细模型,区域2和区域3的光伏并网逆变器模型均使用光伏并网逆变器等效简化模型。
优选地,所述光伏并网逆变器详细模型的拓扑结构由光伏阵列、电容器、逆变器开关部分、三相滤波器和控制系统组成,其中,光伏阵列与电容器并联构成光伏并网逆变器详细模型直流侧部分,经过逆变器开关部分和三相滤波器产生所需交流电,控制系统产生控制信号控制所述逆变器开关部分;按照如下方法在RT-LAB实时仿真平台中构建光伏并网逆变器详细模型:
将光伏阵列、电容器、控制系统建立在CPU上进行微秒级步长的仿真,将逆变器开关部分和三相滤波器建立在FPGA上进行纳秒级步长的仿真。
优选地,所述光伏并网逆变器等效简化模型的拓扑结构由光伏阵列、电容器、受控电流源、受控电压源、三相滤波器和控制系统组成,其中,光伏阵列、电容器和受控电流源并联构成光伏并网逆变器等效简化模型直流侧部分,光伏并网逆变器等效简化模型交流侧部分通过受控电压源和三相滤波器产生所需交流电,控制系统产生控制信号控制受控电流源和受控电压源;按照如下步骤在RT-LAB实时仿真平台中构建光伏并网逆变器等效简化模型:
步骤3.1,采样光伏并网逆变器等效简化模型控制系统的电流内环输出信号Ua、Ub,采样光伏并网逆变器等效简化模型直流侧电压Udc,采样光伏并网逆变器等效简化模型并网点相电压Uabc和相电流Iabc,对光伏并网逆变器等效简化模型并网点相电压Uabc使用abc-dq坐标变换模块得到光伏并网逆变器等效简化模型并网点相电压dq轴分量ed、eq;
步骤3.2,将步骤3.1中得到的电流内环输出信号Ua、Ub,加入光伏并网逆变器等效简化模型并网点相电压dq轴分量ed、eq得到光伏并网逆变器等效简化模型dq坐标系下的电压信号Udref和Uqref,然后使用dq-abc坐标变换模块得到光伏并网逆变器等效简化模型三相电压信号Uaref、Ubref和Ucref,将Uaref、Ubref和Ucref作为光伏并网逆变器等效简化模型交流侧受控电压源的控制信号;
步骤3.3,根据光伏并网逆变器等效简化模型并网点相电压Uabc和相电流Iabc,按如下公式计算光伏并网逆变器等效简化模型交流侧功率输出值Pac,
Pac=3×Uabc×Iabc×cosφ
式中,cosφ是功率因数,令光伏并网逆变器等效简化模型在单位功率因数运行,cosφ=1;
步骤3.4,根据步骤3.3得到的光伏并网逆变器等效简化模型交流侧功率输出值Pac和步骤3.1得到的光伏并网逆变器等效简化模型直流侧电压Udc,计算得到光伏并网逆变器等效简化模型直流侧部分受控电流源的控制信号ih,ih=Pac÷Udc。
优选地,所述同步发电机组模型由发电机模型、调速器模型、原动机模型和励磁系统组成,调速器模型与原动机模型串联输出机械功率给发电机模型,同时励磁系统提供励磁电压给发电机模型;按照如下方法在RT-LAB实时仿真平台中构建同步发电机组模型:
步骤4.1,先给定参考频率ωref,采样发电机模型的输出频率ω,并计算得到频率偏差Δ,Δ=ωref-ω;然后将模型频率偏差Δ依次经过调速器模型和原动机模型,产生机械功率Pm作为发电机模型的输入;
步骤4.2,先给定参考电压Uref,采样发电机模型的端电压Ut,并计算得到电压偏差ε,ε=Uref-Ut;然后将电压偏差ε进入电压调节器进行放大,输出电压UR作为励磁系统的励磁电压。
优选地,所述ZIP负荷模型按照如下方法在RT-LAB实时仿真平台中构建:
步骤5.1,给定负荷有功功率指令值P0,无功功率指令值Q0,基准电压值U0,采样网侧相电压U和负荷频率偏差Δf,计算得到负荷吸收的有功功率PZIP和负荷吸收的无功功率QZIP,其计算式为:
PZIP=P0[p1(U/U0)2+p2(U/U0)+p3](1+kpfΔf)
QZIP=Q0[q1(U/U0)2+q2(U/U0)+q3](1+kqfΔf)
式中:p1是有功功率中恒阻抗负荷占总负荷的百分比,p2是有功功率中恒电流负荷占总负荷的百分比,p3是有功功率中恒功率负荷占总负荷的百分比,q1是无功功率中恒阻抗负荷占总负荷的百分比,q2是无功功率中恒电流负荷占总负荷的百分比,q3是无功功率中恒功率负荷占总负荷的百分比,kpf是有功的频率偏差影响因子,kqf是无功的频率偏差影响因子;
步骤5.2,根据步骤5.1得到的负荷吸收的有功功率PZIP和负荷吸收的无功功率QZIP计算负荷电流dq轴分量Id_ref,Iq_ref,其计算式为:
式中,Ud是网侧相电压U的d轴分量,Uq是网侧相电压U的q轴分量;
步骤5.3,使用dq-abc坐标变换模块将步骤5.2得到的负荷电流dq轴分量Id_ref、Iq_ref转变成负荷三相电流Ia,Ib,Ic。
优选地,所述储能实时仿真子系统包括储能电池模型和储能变流器模型,其中储能电池模型按照如下方法在RT-LAB实时仿真平台中构建:
给定储能电池放电电流ib,使用安时积分法按照下式计算储能电池的荷电状态,
式中soc0为储能电池初始荷电状态值,当电池充满电时为1,soc为储能电池当前荷电状态值,Qb为储能电池额定容量,α为充放电效率,t为放电时间;
通过下式求得储能电池等效电势Eb,
式中E0为蓄电池内电势,K为储能电池等效电势Eb与soc关系曲线的拟合参数1,A为储能电池Eb与soc关系曲线的拟合参数2,B为储能电池Eb与soc关系曲线的拟合参数3。
相对于现有技术,本发明的有益效果如下:
1)考虑了分布式光伏发电系统对电网频率的影响,搭建了同步发电机组模型来对电网环境进行仿真模拟。
2)考虑了源侧子系统、网侧子系统、负荷子系统和储能子系统之间的协调性。
3)基于RT-LAB实时仿真平台,采用了FPGA小步长仿真构建了逆变器详细模型,同时提出等值方法构建逆变器简化模型。既满足仿真精确度,又达到了降低模型复杂度和减小仿真计算时间的效果。
附图说明
图1为本发明高渗透率分布式光伏发电集群框架图;
图2为光伏并网逆变器详细模型的结构示意图;
图3为光伏并网逆变器等效简化模型的结构示意图;
图4为光伏并网逆变器详细模型和光伏并网逆变器等效简化模型的波形对比图;
图5为各种负荷扰动情况下有无调速器的系统频率变化曲线图;
图6为ZIP负荷模型实现框图;
具体实施方式
本发明提出了一种高渗透率分布式光伏发电集群实时仿真建模方法,所述实时仿真建模方法包括基于RT-LAB实时仿真平台构建源侧实时仿真子系统、网侧实时仿真子系统、负荷实时仿真子系统和储能实时仿真子系统;其中,源侧实时仿真子系统由光伏并网逆变器模型组成,所述的光伏并网逆变器模型包括光伏并网逆变器详细模型和光伏并网逆变器等效简化模型;网侧实时仿真子系统由同步发电机组模型组成;负荷实时仿真子系统由ZIP负荷模型组成。
所述的高渗透率分布式光伏发电集群包括3个光伏发电区域,分别记为区域1、区域2和区域3,区域1、区域2、区域3和负荷实时仿真子系统经并联到公共母线后与网侧实时仿真子系统串联。
所述区域1、区域2和区域3中,每个区域均包括7个光伏电站和1个储能实时仿真子系统,7个光伏电站和1个储能实时仿真子系统依次串联。
所述光伏电站中均包含两台并联的光伏并网逆变器,区域1的光伏并网逆变器均使用光伏并网逆变器详细模型,区域2和区域3的光伏并网逆变器模型均使用光伏并网逆变器等效简化模型。
图1为本发明高渗透率分布式光伏发电集群框架图,本实施例中的具体参数如下。在每个光伏电站中,包括两台500KW的光伏逆变器及一台连接方式为Yn11-Yn11-D、变比为0.315/0.315/35KV的双分裂变压器。另外,每个储能实时仿真子系统的额定功率为1MW。
下面分四个部分说明四个实时仿真子系统的建模方法。
(一)源侧实时仿真子系统
源侧实时仿真子系统由光伏并网逆变器模型组成,光伏并网逆变器模型包括光伏并网逆变器详细模型和光伏并网逆变器等效简化模型。
(1)光伏并网逆变器详细模型
所述光伏并网逆变器详细模型的拓扑结构由光伏阵列、电容器、逆变器开关部分、三相滤波器和控制系统组成,其中,光伏阵列与电容器并联构成光伏并网逆变器详细模型直流侧部分,经过逆变器开关部分和三相滤波器产生所需交流电,控制系统产生控制信号控制所述逆变器开关部分。
按照如下方法在RT-LAB实时仿真平台中构建光伏并网逆变器详细模型:
将光伏阵列、电容器、控制系统建立在CPU(中央处理器)上进行微秒级步长的仿真,将逆变器开关部分和三相滤波器建立在FPGA(现场可编程门阵列)上进行纳秒级步长的仿真。CPU模型提供直流侧电压信号、PWM调制信号给FPGA模型,在FPGA模型中采样交流侧电流信号作为反馈量传给CPU模型的控制系统形成闭环,具体结构见图2。
在本实施例中,光伏并网逆变器模型额定功率为500KW,逆变器基波频率ω0为100πrad/s,使用LCL滤波,逆变器侧电感L1为0.1mH,电网侧电感L2为0.05mH,滤波电容C1为200μF,阻尼电阻Rd为0.07Ω,直流侧电压初始值Udc设置为650V,直流侧电容器C为5.5μF。
(2)光伏并网逆变器等效简化模型
光伏并网逆变器等效简化模型的拓扑结构见图3。光伏并网逆变器等效简化模型的拓扑结构由光伏阵列、电容器、受控电流源、受控电压源、三相滤波器和控制系统组成,其中,光伏阵列、电容器和受控电流源并联构成光伏并网逆变器等效简化模型直流侧部分,光伏并网逆变器等效简化模型交流侧部分通过受控电压源和三相滤波器产生所需交流电,控制系统产生控制信号控制受控电流源和受控电流源。
按照如下步骤在RT-LAB实时仿真平台中构建光伏并网逆变器等效简化模型:
步骤3.1,采样光伏并网逆变器等效简化模型控制系统的电流内环输出信号Ua、Ub,采样光伏并网逆变器等效简化模型直流侧电压Udc,采样光伏并网逆变器等效简化模型并网点相电压Uabc和相电流Iabc,对光伏并网逆变器等效简化模型并网点相电压Uabc使用abc-dq坐标变换模块得到光伏并网逆变器等效简化模型并网点相电压dq轴分量ed、eq;
步骤3.2,将步骤3.1中得到的电流内环输出信号Ua、Ub,加入光伏并网逆变器等效简化模型并网点相电压dq轴分量ed、eq得到光伏并网逆变器等效简化模型dq坐标系下的电压信号Udref和Uqref,然后使用dq-abc坐标变换模块得到光伏并网逆变器等效简化模型三相电压信号Uaref、Ubref和Ucref,将Uaref、Ubref和Ucref作为光伏并网逆变器等效简化模型交流侧受控电压源的控制信号;
步骤3.3,根据光伏并网逆变器等效简化模型并网点相电压Uabc和相电流Iabc,按如下公式计算光伏并网逆变器等效简化模型交流侧功率输出值Pac,
Pac=3×Uabc×Iabc×cosφ
式中,cosφ是功率因数,令光伏并网逆变器等效简化模型在单位功率因数运行,cosφ=1;
步骤3.4,根据步骤3.3得到的光伏并网逆变器等效简化模型交流侧功率输出值Pac和步骤3.1得到的光伏并网逆变器等效简化模型直流侧电压Udc,计算得到光伏并网逆变器等效简化模型直流侧部分受控电流源的控制信号ih,ih=Pac÷Udc。
图4为逆变器FPGA详细模型和受控源等效简化模型的波形对比图,可以看到其动态响应一致。假设等效简化模型和详细模型的输出功率分别为a(i)和b(i),详细模型输出功率与时间轴围成的面积为S,等效简化模型和详细模型的输出功率与时间围成的面积之差为ΔS,则定义误差度:其中i表示不断变化的时间量。可计算求得误差度为1.24%,值很小,说明等效简化模型能很好的吻合详细模型。
(二)网侧实时仿真子系统
网侧实时仿真子系统由同步发电机组模型组成。同步发电机组模型由发电机模型、调速器模型、原动机模型和励磁系统组成,调速器模型与原动机模型串联输出机械功率给发电机模型,同时励磁系统提供励磁电压给发电机模型。
按照如下方法在RT-LAB实时仿真平台中构建同步发电机组模型:
步骤4.1,先给定参考频率ωref,采样发电机模型的输出频率ω,并计算得到频率偏差Δ,Δ=ωref-ω;然后将模型频率偏差Δ依次经过调速器模型和原动机模型,产生机械功率Pm作为发电机模型的输入;
步骤4.2,先给定参考电压Uref,采样发电机模型的端电压Ut,并计算得到电压偏差ε,ε=Uref-Ut;然后将电压偏差ε进入电压调节器进行放大,输出电压UR作为励磁系统的励磁电压。
在本实施例中,发电机组额定动率设为42MW,额定电压为110KV,调速器比例增益设置为0.5。
图5为各种负荷扰动情况下有无调速器的系统频率变化曲线图,可以看到负荷波动越大,系统频率变化越大。同时在发电机组没有一次调频功能的情况下频率波动远远大于有一次调频功能的电力系统,甚至在负荷变化5%的情况下频率波动已经超过了±0.5Hz的允许波动界限。可以看出用搭建的同步发电机组模型模拟电网可以反映电网的频率波动情况。
(三)负荷实时仿真子系统
负荷实时仿真子系统由ZIP负荷模型构成。图6为ZIP负荷模型实现框图,由该图可见,ZIP负荷模型按照如下方法在RT-LAB实时仿真平台中构建:
步骤5.1,给定负荷有功功率指令值P0,无功功率指令值Q0,基准电压值U0,采样网侧相电压U和负荷频率偏差Δf,计算得到负荷吸收的有功功率PZIP和负荷吸收的无功功率QZIP,其计算式为:
PZIP=P0[p1(U/U0)20p2(U/U0)0p3](10kpfΔf)
QZIP=Q0[q1(U/U0)2+q2(U/U0)+q3](10kqfΔf.
式中:p1是有功功率中恒阻抗负荷占总负荷的百分比,p2是有功功率中恒电流负荷占总负荷的百分比,p3是有功功率中恒功率负荷占总负荷的百分比,q1是无功功率中恒阻抗负荷占总负荷的百分比,q2是无功功率中恒电流负荷占总负荷的百分比,q3是无功功率中恒功率负荷占总负荷的百分比,kpf是有功的频率偏差影响因子,kqf是无功的频率偏差影响因子;
步骤5.2,根据步骤5.1得到的负荷吸收的有功功率PZIP和负荷吸收的无功功率QZIP计算负荷电流dq轴分量Id_ref,Iq_ref,其计算式为:
式中,Ud是网侧相电压U的d轴分量,Uq是网侧相电压U的q轴分量;
步骤5.3,使用dq-abc坐标变换模块将步骤5.2得到的负荷电流dq轴分量Id_ref、Iq_ref转变成负荷三相电流Ia,Ib,Ic。
在本实施例中,ZIP负荷模型的额定功率设置为42MW,p1=q1=0.3,p2=q2=0.3,p3=q3=0.4,kpf=1,kqf=-1。
(四)储能实时仿真子系统
所述储能实时仿真子系统包括储能电池模型和储能变流器模型,其中储能电池模型按照如下方法在RT-LAB实时仿真平台中构建:
给定储能电池放电电流ib,使用安时积分法按照下式计算储能电池的荷电状态,
式中soc0为储能电池初始荷电状态值,当电池充满电时为1,soc为储能电池当前荷电状态值,Qb为储能电池额定容量,α为充放电效率,t为放电时间;
通过下式求得储能电池等效电势Eb,
式中E0为蓄电池内电势,K为储能电池等效电势Eb与soc关系曲线的拟合参数1,A为储能电池Eb与soc关系曲线的拟合参数2,B为储能电池Eb与soc关系曲线的拟合参数3。
在本实施例中,E0=660V,R=2.35Ω,拟合参数1=0.0188,拟合参数2=0.0144,拟合参数3=2.3077。
将储能电池模型连接储能变流器模型构成储能实时仿真子系统。
综上可得,该建模方法实现简单,在仿真精确度要求高和仿真规模大的条件下,提供了一种切实可行的高渗透率光伏发电集群实时仿真建模方法。
Claims (6)
1.一种高渗透率分布式光伏发电集群实时仿真建模方法,其特征在于,所述实时仿真建模方法包括基于RT-LAB实时仿真平台构建源侧实时仿真子系统、网侧实时仿真子系统、负荷实时仿真子系统和储能实时仿真子系统;其中,源侧实时仿真子系统由光伏并网逆变器模型组成,所述的光伏并网逆变器模型包括光伏并网逆变器详细模型和光伏并网逆变器等效简化模型;网侧实时仿真子系统由同步发电机组模型组成;负荷实时仿真子系统由ZIP负荷模型组成;
所述的高渗透率分布式光伏发电集群包括3个光伏发电区域,分别记为区域1、区域2和区域3,区域1、区域2、区域3和负荷实时仿真子系统彼此并联到公共母线后与网侧实时仿真子系统串联;
所述区域1、区域2和区域3中,每个区域均包括7个光伏电站和1个储能实时仿真子系统,7个光伏电站和1个储能实时仿真子系统依次串联;
所述光伏电站中均包含两台并联的光伏并网逆变器,区域1的光伏并网逆变器均使用光伏并网逆变器详细模型,区域2和区域3的光伏并网逆变器模型均使用光伏并网逆变器等效简化模型。
2.根据权利要求1所述一种高渗透率分布式光伏发电集群实时仿真建模方法,其特征在于,所述光伏并网逆变器详细模型的拓扑结构由光伏阵列、电容器、逆变器开关部分、三相滤波器和控制系统组成,其中,光伏阵列与电容器并联构成光伏并网逆变器详细模型直流侧部分,经过逆变器开关部分和三相滤波器产生所需交流电,控制系统产生控制信号控制所述逆变器开关部分;按照如下方法在RT-LAB实时仿真平台中构建光伏并网逆变器详细模型:
将光伏阵列、电容器、控制系统建立在CPU上进行微秒级步长的仿真,将逆变器开关部分和三相滤波器建立在FPGA上进行纳秒级步长的仿真。
3.根据权利要求1所述一种高渗透率分布式光伏发电集群实时仿真建模方法,其特征在于,所述光伏并网逆变器等效简化模型的拓扑结构由光伏阵列、电容器、受控电流源、受控电压源、三相滤波器和控制系统组成,其中,光伏阵列、电容器和受控电流源并联构成光伏并网逆变器等效简化模型直流侧部分,光伏并网逆变器等效简化模型交流侧部分通过受控电压源和三相滤波器产生所需交流电,控制系统产生控制信号控制受控电流源和受控电压源;按照如下步骤在RT-LAB实时仿真平台中构建光伏并网逆变器等效简化模型:
步骤3.1,采样光伏并网逆变器等效简化模型控制系统的电流内环输出信号Ua、Ub,采样光伏并网逆变器等效简化模型直流侧电压Udc,采样光伏并网逆变器等效简化模型并网点相电压Uabc和相电流Iabc,对光伏并网逆变器等效简化模型并网点相电压Uabc使用abc-dq坐标变换模块得到光伏并网逆变器等效简化模型并网点相电压dq轴分量ed、eq;
步骤3.2,将步骤3.1中得到的电流内环输出信号Ua、Ub,加入光伏并网逆变器等效简化模型并网点相电压dq轴分量ed、eq得到光伏并网逆变器等效简化模型dq坐标系下的电压信号Udref和Uqref,然后使用dq-abc坐标变换模块得到光伏并网逆变器等效简化模型三相电压信号Uaref、Ubref和Ucref,将Uaref、Ubref和Ucref作为光伏并网逆变器等效简化模型交流侧受控电压源的控制信号;
步骤3.3,根据光伏并网逆变器等效简化模型并网点相电压Uabc和相电流Iabc,按如下公式计算光伏并网逆变器等效简化模型交流侧功率输出值Pac,
Pac=3×Uabc×Iabc×cosφ
式中,cosφ是功率因数,令光伏并网逆变器等效简化模型在单位功率因数运行,cosφ=1;
步骤3.4,根据步骤3.3得到的光伏并网逆变器等效简化模型交流侧功率输出值Pac和步骤3.1得到的光伏并网逆变器等效简化模型直流侧电压Udc,计算得到光伏并网逆变器等效简化模型直流侧部分受控电流源的控制信号ih,ih=Pac÷Udc。
4.根据权利要求1所述一种高渗透率分布式光伏发电集群实时仿真建模方法,其特征在于,所述同步发电机组模型由发电机模型、调速器模型、原动机模型和励磁系统组成,调速器模型与原动机模型串联输出机械功率给发电机模型,同时励磁系统提供励磁电压给发电机模型;按照如下方法在RT-LAB实时仿真平台中构建同步发电机组模型:
步骤4.1,先给定参考频率ωref,采样发电机模型的输出频率ω,并计算得到频率偏差Δ,Δ=ωref-ω;然后将模型频率偏差Δ依次经过调速器模型和原动机模型,产生机械功率Pm作为发电机模型的输入;
步骤4.2,先给定参考电压Uref,采样发电机模型的端电压Ut,并计算得到电压偏差ε,ε=Uref-Ut;然后将电压偏差ε进入电压调节器进行放大,输出电压UR作为励磁系统的励磁电压。
5.根据权利要求1所述一种高渗透率分布式光伏发电集群实时仿真建模方法,其特征在于,所述ZIP负荷模型按照如下方法在RT-LAB实时仿真平台中构建:
步骤5.1,给定负荷有功功率指令值P0,无功功率指令值Q0,基准电压值U0,采样网侧相电压U和负荷频率偏差Δf,计算得到负荷吸收的有功功率PZIP和负荷吸收的无功功率QZIP,其计算式为:
PZIP=P0[p1(U/U0)2+p2(U/U0)+p3](1+kpfΔf)
QZIP=Q0[q1(U/U0)2+q2(U/U0)+q3](1+kqfΔf)
式中:p1是有功功率中恒阻抗负荷占总负荷的百分比,p2是有功功率中恒电流负荷占总负荷的百分比,p3是有功功率中恒功率负荷占总负荷的百分比,q1是无功功率中恒阻抗负荷占总负荷的百分比,q2是无功功率中恒电流负荷占总负荷的百分比,q3是无功功率中恒功率负荷占总负荷的百分比,kpf是有功的频率偏差影响因子,kqf是无功的频率偏差影响因子;
步骤5.2,根据步骤5.1得到的负荷吸收的有功功率PZIP和负荷吸收的无功功率QZIP计算负荷电流dq轴分量Id_ref,Iq_ref,其计算式为:
式中,Ud是网侧相电压U的d轴分量,Uq是网侧相电压U的q轴分量;
步骤5.3,使用dq-abc坐标变换模块将步骤5.2得到的负荷电流dq轴分量Id_ref、Iq_ref转变成负荷三相电流Ia,Ib,Ic。
6.根据权利要求1所述一种高渗透率分布式光伏发电集群实时仿真建模方法,其特征在于,所述储能实时仿真子系统包括储能电池模型和储能变流器模型,其中储能电池模型按照如下方法在RT-LAB实时仿真平台中构建:
给定储能电池放电电流ib,使用安时积分法按照下式计算储能电池的荷电状态,
式中soc0为储能电池初始荷电状态值,当电池充满电时为1,soc为储能电池当前荷电状态值,Qb为储能电池额定容量,α为充放电效率,t为放电时间;
通过下式求得储能电池等效电势Eb,
式中E0为蓄电池内电势,K为储能电池等效电势Eb与soc关系曲线的拟合参数1,A为储能电池Eb与soc关系曲线的拟合参数2,B为储能电池Eb与soc关系曲线的拟合参数3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810392773.7A CN108649560B (zh) | 2018-04-27 | 2018-04-27 | 高渗透率分布式光伏发电集群实时仿真建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810392773.7A CN108649560B (zh) | 2018-04-27 | 2018-04-27 | 高渗透率分布式光伏发电集群实时仿真建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108649560A true CN108649560A (zh) | 2018-10-12 |
CN108649560B CN108649560B (zh) | 2020-06-23 |
Family
ID=63747868
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810392773.7A Active CN108649560B (zh) | 2018-04-27 | 2018-04-27 | 高渗透率分布式光伏发电集群实时仿真建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108649560B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109167395A (zh) * | 2018-11-22 | 2019-01-08 | 国网宁夏电力有限公司电力科学研究院 | 基于adpss的光伏发电系统的暂态等值电势辨识方法 |
CN111666658A (zh) * | 2020-05-15 | 2020-09-15 | 中国能源建设集团甘肃省电力设计院有限公司 | 大型光伏发电系统中多台光伏发电装置简化模型建模方法 |
CN111817333A (zh) * | 2020-06-03 | 2020-10-23 | 中国电力科学研究院有限公司 | 一种基于柔性直流输电系统特性的交流过电压计算方法及系统 |
CN113258607A (zh) * | 2021-06-17 | 2021-08-13 | 昆明理工大学 | 基于不同渗透下的光伏发电系统动态离散等值模型建立方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103973203A (zh) * | 2014-02-18 | 2014-08-06 | 国家电网公司 | 一种适用于安全稳定分析的大型光伏电站在线等值建模方法 |
CN105938578A (zh) * | 2016-04-13 | 2016-09-14 | 合肥工业大学 | 一种基于聚类分析的大规模光伏电站等值建模方法 |
CN106571644A (zh) * | 2016-11-04 | 2017-04-19 | 国网新疆电力公司经济技术研究院 | 一种调节汽轮机参数以提高光伏渗透率的仿真方法 |
CN107579541A (zh) * | 2017-08-31 | 2018-01-12 | 南京理工大学 | 一种基于模式分析的光伏电站对多机系统低频振荡的抑制方法 |
US20180102651A1 (en) * | 2016-05-26 | 2018-04-12 | Ge Energy Power Conversion Technology Ltd | Solar inverter grid emulation mode |
-
2018
- 2018-04-27 CN CN201810392773.7A patent/CN108649560B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103973203A (zh) * | 2014-02-18 | 2014-08-06 | 国家电网公司 | 一种适用于安全稳定分析的大型光伏电站在线等值建模方法 |
CN105938578A (zh) * | 2016-04-13 | 2016-09-14 | 合肥工业大学 | 一种基于聚类分析的大规模光伏电站等值建模方法 |
US20180102651A1 (en) * | 2016-05-26 | 2018-04-12 | Ge Energy Power Conversion Technology Ltd | Solar inverter grid emulation mode |
CN106571644A (zh) * | 2016-11-04 | 2017-04-19 | 国网新疆电力公司经济技术研究院 | 一种调节汽轮机参数以提高光伏渗透率的仿真方法 |
CN107579541A (zh) * | 2017-08-31 | 2018-01-12 | 南京理工大学 | 一种基于模式分析的光伏电站对多机系统低频振荡的抑制方法 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109167395A (zh) * | 2018-11-22 | 2019-01-08 | 国网宁夏电力有限公司电力科学研究院 | 基于adpss的光伏发电系统的暂态等值电势辨识方法 |
CN111666658A (zh) * | 2020-05-15 | 2020-09-15 | 中国能源建设集团甘肃省电力设计院有限公司 | 大型光伏发电系统中多台光伏发电装置简化模型建模方法 |
CN111666658B (zh) * | 2020-05-15 | 2023-06-23 | 中国能源建设集团甘肃省电力设计院有限公司 | 大型光伏发电系统中多台光伏发电装置简化模型建模方法 |
CN111817333A (zh) * | 2020-06-03 | 2020-10-23 | 中国电力科学研究院有限公司 | 一种基于柔性直流输电系统特性的交流过电压计算方法及系统 |
CN111817333B (zh) * | 2020-06-03 | 2022-08-12 | 中国电力科学研究院有限公司 | 一种基于柔性直流输电系统特性的交流过电压计算方法及系统 |
CN113258607A (zh) * | 2021-06-17 | 2021-08-13 | 昆明理工大学 | 基于不同渗透下的光伏发电系统动态离散等值模型建立方法 |
CN113258607B (zh) * | 2021-06-17 | 2022-09-20 | 昆明理工大学 | 基于不同渗透下的光伏发电系统动态离散等值模型建立方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108649560B (zh) | 2020-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Miao et al. | Investigation of microgrids with both inverter interfaced and direct AC-connected distributed energy resources | |
CN108649560A (zh) | 高渗透率分布式光伏发电集群实时仿真建模方法 | |
CN104659804A (zh) | 含有混合储能的微电网及其控制方法 | |
CN106549417A (zh) | 一种光伏‑储能系统的虚拟同步发电机控制方法及装置 | |
CN102842917B (zh) | 一种通用的并网式光伏发电系统机电暂态模型 | |
CN103124080B (zh) | 一种光伏发电系统模型的建模方法 | |
CN105608244B (zh) | 基于pscad/emtdc的分布式光伏发电系统 | |
CN104932285A (zh) | 一种光伏发电系统的等效建模方法 | |
CN104979849A (zh) | 一种适用于用户侧微电网的并网运行控制方法 | |
Akbari et al. | Voltage control of a hybrid ac/dc microgrid in grid-connected operation mode | |
CN107026468A (zh) | 一种基于自抗扰控制技术的微网实验系统 | |
CN103840695A (zh) | 一种光伏并网逆变器控制参数的辨识方法 | |
CN105680473A (zh) | 面向光伏发电系统的通用机电暂态信息物理融合建模方法 | |
Wu et al. | Influence of induction motor starting on microgrid | |
Jian et al. | VSG current balance control strategy under unbalanced grid voltage | |
CN203218889U (zh) | 一种通用的并网式光伏发电系统机电暂态模型 | |
CN108365620A (zh) | 一种储能系统基于等值模型进行仿真的方法及装置 | |
Yi et al. | Modeling and simulation of DC micro-grid based on virtual motor control | |
CN201378744Y (zh) | 一种三相光伏并网逆变控制装置 | |
CN111293713A (zh) | 一种风电场风机变流器仿真模型的建立方法及仿真方法 | |
Narkhede et al. | Challenges, modeling simulation and performance analysis of virtual power plant in Indian context | |
Ramabhotla et al. | RSCAD/RTDS | |
Lin et al. | Converter controller design methods for wind turbine systems in the DC micro grid | |
Shi et al. | Design of a real-time simulation platform to support large-scaled renewable energy access to power grid | |
Xie et al. | An Advanced DC Voltage Control for Virtual Synchronous Generator |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |