CN108631363B - 基于主从式控制的逆变器光伏虚拟同步控制方法 - Google Patents
基于主从式控制的逆变器光伏虚拟同步控制方法 Download PDFInfo
- Publication number
- CN108631363B CN108631363B CN201810462242.0A CN201810462242A CN108631363B CN 108631363 B CN108631363 B CN 108631363B CN 201810462242 A CN201810462242 A CN 201810462242A CN 108631363 B CN108631363 B CN 108631363B
- Authority
- CN
- China
- Prior art keywords
- inverter
- power
- photovoltaic
- voltage
- pvi
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 82
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 51
- 238000004891 communication Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims description 29
- 230000001939 inductive effect Effects 0.000 claims description 28
- 230000009466 transformation Effects 0.000 claims description 17
- 238000005070 sampling Methods 0.000 claims description 15
- 230000005611 electricity Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 4
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000007665 sagging Methods 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 abstract 1
- 230000001965 increasing effect Effects 0.000 description 9
- 241000196324 Embryophyta Species 0.000 description 6
- 238000013016 damping Methods 0.000 description 5
- 238000004146 energy storage Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000011217 control strategy Methods 0.000 description 4
- 230000003247 decreasing effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- SJWPTBFNZAZFSH-UHFFFAOYSA-N pmpp Chemical compound C1CCSC2=NC=NC3=C2N=CN3CCCN2C(=O)N(C)C(=O)C1=C2 SJWPTBFNZAZFSH-UHFFFAOYSA-N 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- H02J3/385—
-
- 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/40—Synchronising a generator for connection to a network or to another generator
- H02J3/42—Synchronising a generator for connection to a network or to another generator with automatic parallel connection when synchronisation is achieved
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Inverter Devices (AREA)
- Control Of Electrical Variables (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了基于主从式控制的逆变器光伏虚拟同步控制方法。该方法将光伏电站的光伏逆变器设置为主从控制模式，主逆变器执行最大功率追踪控制方法，从逆变器执行光伏虚拟同步控制方法，在主逆变器获得最大功率点的功率和电压通过通信发送给从逆变器，使从逆变器调整功率追踪的步长和调频参数。本发明改进原有光伏虚拟同步控制方法，使从逆变器始终处于调频状态，对电网频率有更好的调节效果。
Description
技术领域
本发明属于光伏逆变器并网控制领域，具体涉及一种基于主从式控制的逆变器光伏虚拟同步控制方法。
背景技术
随着分布式发电技术逐步成熟，光伏、风电等新能源发电应用越来越广泛，但是分布式发电系统是通过逆变器等电力电子装置输出电能与电网相连，虽然暂态响应速度快，但其不具备惯性，也无法支撑电网的一次调频。大规模新能源的并网时电网渗透率不断增加，传统同步发电机的装机比例不断下降，这导致电力系统旋转惯量和系统阻尼的降低，使电网鲁棒性变弱。而且，光伏发电系统多采用最大功率跟踪方式并网，其“只管发电，不管电网”的特性，也导致电网能量过剩，加剧了电网频率的恶化。
近年来国内外学者提出了虚拟同步发电机(VSG)的控制方法，借鉴传统电力系统中对电网天然友好的同步发电机的运行特点和控制方式，使光伏并网系统模拟同步发电机的惯性和阻尼特性，有助于其实现友好并网，虚拟同步技术借鉴其运行特性以及控制方式，模拟传统同步发电机的转动惯性和阻尼的外特性，实现友好并网。现今高渗透率电网的功率供需关系不再要求光伏发电系统时刻处于最大功率追踪状态，也迫切希望其能够参与到电网调频调压的过程中来，如何依托灵活可控的电力电子技术提高光伏系统的可靠性，使其成为共同维护电网安全稳定的“参与者”，具有非常重要的意义。
现有光伏虚拟同步技术改进传统光伏逆变器的控制策略，其控制独立，相互不干扰，但是由于为保证系统稳定，需隔一段时间进行最大功率追踪来确定备用功率点，然后在备用功率点处调频，以确保输出功率不会超过最大功率，但是由于不断在调频状态切换，最大功率追踪过程不能调频，影响调频效果。
目前，对于光伏虚拟同步发电机技术，已有多篇学术论文进行分析并提出解决方案，例如：
1、题为(“基于虚拟同步发电机的光伏并网系统及仿真分析”，毛弋，向海燕，康伦，《湖南师范大学自然科学学报》，2013年第36卷第3期27-32页)文章中将虚拟同步控制策略应用于光伏并网系统，使光伏逆变器具有同步发电机的外特性，但文章对系统的稳定性问题并未做深入研究，同时也不能使光伏逆变器的出力灵活可控。
2、题为“A consensus-based frequency control for photovoltaic connectedvirtual synchronous generators in microgrid”，Guo Y,Chen L,《IEEE InternationalConference on Power System Technology》2016:1-6.(“一种基于一致性频率控制的光伏微网虚拟同步发电机”，郭岩，陈来军，IEEE电力系统技术国际会议，2016年第1-6页)文章将虚拟同步技术和光伏逆变器的最大功率追踪技术结合，使光伏逆变器在没有储能的情况下具有调频调压功能。但是文章存在以下不足：该文并未对最大功率追踪方法进行改进，使系统稳定性依靠限幅PI控制器来调节，但是参数不易整定。
3、题为(“考虑源端动态特性的光伏虚拟同步机多模式运行控制”，郑天文,陈来军,刘炜,《中国电机工程学报》，2017年第37卷第2期第454-463页)文章中增加了虚拟同步发电机的电磁部分和机电特性，在运行特性上和传统同步发电机更接近。但此方法存在对光伏虚拟同步发电机的稳定性未作过多研究，针对在功率和电压上额外增加一个PI控制器，控制器的参数不易整定，且当光伏功率不足时，并网状态减小电压幅值容易造成逆变器脱网。
发明内容
本发明目的是基于拉格朗日插值预测的光伏虚拟同步控制方法需离线测试最大功率点，且单个无储能光伏虚拟同步机需要间隔一段时间进行最大功率追踪以刷新备用功率点及调频系数，此期间不能调频的问题，提供一种基于主从式控制的逆变器光伏虚拟同步控制方法，本方法针对现有较强通信连接的光伏电站，对光伏电站部分光伏逆变器加以改进，使整个光伏电站具有调频调压的功能。
为实现上述目的，提出一种基于主从式控制的逆变器光伏虚拟同步控制方法，该方法涉及的逆变器台数为n，n为整数且n≥2，将逆变器分为主逆变器和从逆变器，所有逆变器之间采用通讯方式连接，主逆变器编号为#1，从逆变器编号为#i，i∈[2,n]且i为正整数，主逆变器采用最大功率追踪控制方法，从逆变器采用光伏虚拟同步控制方法；
从逆变器采用的光伏虚拟同步控制方法的步骤如下：
步骤1、采样主逆变器#1光伏电池板最大功率点的功率PMPP、光伏电池板最大功率点的电压UMPP和开路电压UOC，并经通讯下发给从逆变器；
步骤2、采样从逆变器#i输出相电压Uoai,Uobi,Uoci，并经输出相电压坐标变换方程得到输出电压dq轴分量Uodi,Uoqi，采样从逆变器#i网侧电感电流Ioai,Iobi,Ioci，并经网侧电感电流坐标变换方程得到网侧电感电流dq轴分量Iodi,Ioqi，采样从逆变器#i桥臂电感电流ILai,ILbi,ILci，并经桥臂电感电流坐标变换方程得到桥臂电感电流dq轴分量ILdi,ILqi，其中d轴为有功轴，q轴为无功轴；通过锁相环获得电网频率ωg，采样从逆变器#i在k时刻的直流侧光伏输出电流Ipvi(k)和从逆变器#i在k时刻的直流侧光伏输出电压Upvi(k)；
步骤3、设置电网角频率基准值ωref和从逆变器#i上层有功功率指令Pref0i，并根据步骤1得到的最大功率点的功率PMPP和步骤2得到的电网频率ωg，通过功率叠加方程获得功率指令其表达式为：
其中，Ji为从逆变器#i的虚拟转动惯量，s为拉普拉斯算子，ki为从逆变器#i的预设定平均有功功率输出值占当前最大功率点功率PMPP的百分比，π为圆周率；
步骤4、根据步骤3中得到的功率指令通过拉格朗日插值预测法得到期望参考功率指令
步骤5、根据步骤2中得到的输出电压dq轴分量Uodi,Uoqi和网侧电感电流dq轴分量Iodi,Ioqi，通过功率计算方程获得平均有功功率Poi；
步骤6、根据步骤2中得到的从逆变器#i在k时刻的直流侧光伏输出电流Ipvi(k)和从逆变器#i在k时刻的直流侧光伏输出电压Upvi(k)、步骤4中得到的期望参考功率指令和步骤5中得到的平均有功功率Poi，通过基于多项式拟合的功率跟踪法获得直流侧参考电压Urefi；
步骤7、根据步骤2中得到的从逆变器#i在k时刻的直流侧光伏输出电压Upvi(k)和步骤6中得到的直流侧参考电压Urefi，经过电压环获得光伏板参考功率Prefi；
步骤8、根据步骤7中得到的光伏板参考功率Prefi和步骤5中得到的平均有功功率Poi，经过功率外环得到d轴有功功率指令Pdrefi；设置指令电压Edref和从逆变器#i上层无功功率指令Qref0i，并根据步骤2中得到的输出电压d轴分量Uodi，通过一次调压方程获得无功功率指令Qdrefi；
步骤9、根据步骤2中得到的输出电压d轴分量Uodi和步骤8中得到的d轴有功功率指令Pdrefi、无功功率指令Qdrefi，经过电流计算方法得到电流d轴指令Icdrefi和电流q轴指令Icqrefi；
步骤10、将步骤9中得到的电流d轴指令Icdrefi和步骤2中得到的桥臂电感电流d轴分量ILdi，经过d轴电流闭环控制方程，得到d轴输出信号Uidi；将步骤9中得到的电流q轴指令Icqrefi与步骤2中得到的桥臂电感电流q轴分量ILqi，经过q轴电流闭环控制方程，得到q轴输出信号Uiqi；
步骤11、将步骤2中得到的输出电压d轴分量Uodi和输出电压q轴分量Uoqi分别加上步骤10中得到的d轴输出信号Uidi和q轴输出信号Uiqi，得到dq坐标系下的调制波Umdi和Umqi，其表达式分别为：
Umdi＝Uodi+Uidi
Umqi＝Uoqi+Uiqi；
步骤12、将步骤11中得到的dq坐标系下的调制波Umdi和Umqi经坐标反变换方程得到从逆变器桥臂电压的三相调制波Umai,Umbi,Umci，经PWM调制后作为IGBT电路的驱动信号。
优选的，步骤2中所述所述电网频率ωg的计算公式为：
其中，ω0为公共耦合点电压的额定角频率，Kp_pll为锁相环比例积分调节器的比例调节系数，Ki_pll为锁相环比例积分调节器的积分调节系数，s为拉普拉斯算子。
优选的，步骤4中所述拉格朗日插值预测法包含以下步骤：
步骤4.1，对步骤3中得到的功率指令进行采样处理，采样时间为拉格朗日插值时间T，并按时间先后顺序记录三次采样获得的功率指令：
其中，为从逆变器#i在k-2时刻的功率指令，为从逆变器#i在k-1时刻的功率指令，为从逆变器#i在k时刻的功率指令；
步骤4.2，计算期望参考功率指令其计算公式如下：
其中，t为预测的时间间隔，t取值为t＝2T，T为拉格朗日插值时间，T取值为T＝5Ts，Ts为从逆变器#i的ADC采样时间。
优选的，步骤5中的通过功率计算方程获得平均有功功率Poi的计算公式为：
其中Tf为一阶低通滤波器的时间常数，s为拉普拉斯算子。
优选的，步骤6所述基于多项式拟合的功率跟踪法包括以下步骤：
步骤6.1，求电压步长Ustep，其计算公式为：
当时，取Ustep＝0；
当时，取Ustep＝Uthreshold_high；
当时，取
其中，为未限幅电压步长，Ustep为电压步长，Uthreshold_low为低电压步长阈值，Uthreshold_high为高电压步长阈值；
步骤6.2，计算从逆变器#i在k时刻光伏电池板的功率Ppvi(k)，其计算公式为：
Ppvi(k)＝Upvi(k)·Ipvi(k)
其中，k为时刻标志；
步骤6.3，计算从逆变器#i的符号标志flagi，其计算公式为：
flagi＝sign(Ppvi(k)-Ppvi(k-1))×sign(Upvi(k)-Upvi(k-1))
其中，flagi为从逆变器#i的符号标志，Ppvi(k)为从逆变器#i在k时刻光伏电池板的功率，Ppvi(k-1)为从逆变器#i在k-1时刻光伏电池板的功率，Upvi(k)为从逆变器#i在k时刻直流侧光伏输出电压，Upvi(k-1)为从逆变器#i在k-1时刻直流侧光伏输出电压，sign为符号函数数学运算符，其含义如下：
其中，x为自变量；
步骤6.4，求直流侧参考电压Urefi；
当期望参考功率指令大于平均有功功率Poi时，执行Urefi＝Upvi(k)+Ustep×flagi；
当期望参考功率指令小于等于平均有功功率Poi时，如果flagi大于等于零，执行Urefi＝Upvi(k)+Ustep，如果flagi小于零，执行Urefi＝Upvi(k)-Ustep×flagi。
优选的，步骤7中的经过电压环获得光伏板参考功率Prefi的表达式为：
Prefi＝(Upvi(k)-Urefi)Gdc(s)
其中，Gdc(s)为直流侧电压闭环比例积分调节器，其表达式为：
Gdc(s)＝kdcki+kdcpi/s
其中，kdcki为从逆变器#i直流侧电压闭环比例调节器系数，kdcpi为从逆变器#i直流侧电压闭环积分调节器系数，s为拉普拉斯算子。
优选的，步骤8所述功率外环和一次调压方程的表达式分别为：
Pdrefi＝(Prefi-Poi)Gp(s)
其中，ni为从逆变器#i无功功率下垂系数，Gp(s)为功率闭环比例积分调节器，其表达式为：
Gp(s)＝kpki+kppi/s
其中，kpki为从逆变器#i功率闭环比例调节器系数，kppi为从逆变器#i功率闭环积分调节器系数，s为拉普拉斯算子。
优选的，步骤9所述电流计算方法得到电流d轴指令Icdrefi和电流q轴指令Icqrefi的表达式分别为：
优选的，步骤10所述d轴电流闭环控制方程和q轴电流闭环控制方程的表达式分别为：
Uidi＝(Icdrefi-ILdi)GI(s)
Uiqi＝(Icqrefi-ILqi)GI(s)
其中，Uidi为d轴输出信号，Uiqi为q轴输出信号，GI(s)为电流闭环比例调节器，其表达式为：
GI(s)＝kIki+kIpi/s
其中，kIki为从逆变器#i电流闭环比例调节器系数,kIpi为从逆变器#i电流闭环积分调节器系数。
本发明公开的基于主从式控制的逆变器光伏虚拟同步控制方法，与现有的光伏并网逆变器相比，其有益效果体现在：
基于拉格朗日插值预测的光伏虚拟同步控制方法需离线测试最大功率点，且单个无储能光伏虚拟同步机需要间隔一段时间进行最大功率追踪以刷新备用功率点及调频系数，此期间不能调频的问题，
1、本控制方法不需要离线测试最大功率点，通过主逆变器获知，减少工作量，并减小误差；
2、本控制方法解决了单个无储能光伏虚拟同步机需要间隔一段时间进行最大功率追踪以刷新备用功率点及调频系数，此期间不能调频的问题；
3、本控制方法仅对现有的光伏逆变器的控制策略进行改进，无需增加额外储能设备，节约成本；
4、本控制方法使光伏逆变器具有更好的调频效果，并实现系统稳定安全可靠的运行，实现友好并网；
5、本控制方法基于多项式拟合的功率追踪，实现了对功率的快速精确的追踪，提高了光伏发电系统的快速性和稳态精度，对同步发电机的转动惯量和阻尼模拟更精确；
6、本控制方法容易对光伏电站进行实现改造，在光伏电站预留一部分逆变器，对控制策略加以改进，使整个电站具有调频调压能力；
7、本控制方法实现了光伏逆变器功率输出灵活可控，使光伏逆变器兼顾PQ并网优势；
8、本控制方法依托现有光伏电站较强的通信能力，节约成本。
附图说明
图1为本发明实施例光伏电站并网结构图。
图2为本发明实施例从逆变器控制结构框图。
图3为本发明实施例从逆变器光伏虚拟同步机多项式拟合示意图。
图4为本发明实施例从逆变器加入采用拉格朗日插值预测法的光伏虚拟同步发电机前后增减负荷电网频率变化波形。
图5为本发明实施例从逆变器加入采用拉格朗日插值预测法的光伏虚拟同步发电机后增减负荷输出有功功率波形。
具体实施方式
下面结合附图对本实施例进行具体的描述。
图1为本发明实施例光伏电站并网结构图。所有逆变器具体参数如下：额定输出线电压为380V/50Hz，直流侧滤波电容Cin＝15mF，桥臂侧滤波电感Lf＝0.06mH，交流测滤波电容Cf＝300uF，网侧滤波电感值为Lg＝0.02mH，逆变器额定容量为500KVA。每个逆变器与电网联线阻抗均为ZL＝0.001+j0.001Ω。本实施例中m＝1，i＝1。
图2为本发明实施例从逆变器控制结构框图。本发明控制方法的步骤如下：
一种基于主从式控制的逆变器光伏虚拟同步控制方法，该方法涉及的逆变器的台数为n，n为整数且n≥2，将逆变器分为主逆变器和从逆变器，所有逆变器之间采用通讯方式连接，主逆变器编号为#1，从逆变器编号为#i，i∈[2,n]且i为正整数，主逆变器采用最大功率追踪控制方法，从逆变器采用光伏虚拟同步控制方法。
从逆变器采用的光伏虚拟同步控制方法的步骤如下：
步骤1、采样主逆变器#1光伏电池板最大功率点的功率PMPP、光伏电池板最大功率点的电压UMPP和开路电压UOC，并经通讯下发给从逆变器；
本实施例中，主逆变器处于温度T＝25℃和光照强度Wd＝1000W/m2下，PMPP＝500KW，UMPP＝645.4V，Uoc＝839.2V。
步骤2、采样从逆变器#i输出相电压Uoai,Uobi,Uoci，并经输出相电压坐标变换方程得到输出电压dq轴分量Uodi,Uoqi，采样从逆变器#i网侧电感电流Ioai,Iobi,Ioci，并经网侧电感电流坐标变换方程得到网侧电感电流dq轴分量Iodi,Ioqi，采样从逆变器#i桥臂电感电流ILai,ILbi,ILci，并经桥臂电感电流坐标变换方程得到桥臂电感电流dq轴分量ILdi,ILqi，其中d轴为有功轴，q轴为无功轴；通过锁相环获得电网频率ωg，采样从逆变器#i在k时刻的直流侧光伏输出电流Ipvi(k)和从逆变器#i在k时刻的直流侧光伏输出电压Upvi(k)；
所述电网频率ωg的计算公式为：
其中，ω0为公共耦合点电压的额定角频率，Kp_pll为锁相环PI调节器的比例调节系数，Ki_pll为锁相环PI调节器的积分调节系数，s为拉普拉斯算子。本实施例中，锁相环中ω0＝100πrad/s，Kp_pll＝0.5，Ki_pll＝1。
所述输出相电压坐标变换方程的表达式为：
所述网侧电感电流坐标变换方程的表达式为：
所述桥臂电感电流坐标变换方程的表达式为：
以上六个公式中，θ为d轴和q轴的相位差，
步骤3、设置电网角频率基准值ωref和从逆变器#i上层有功功率指令Pref0i，并根据步骤1得到的最大功率点的功率PMPP和步骤2得到的电网频率ωg，通过功率叠加方程获得功率指令其表达式为：
其中，Ji为从逆变器#i的虚拟转动惯量，s为拉普拉斯算子，ki为从逆变器#i的预设定平均有功功率输出值占当前最大功率点功率PMPP的百分比，π为圆周率，本实施例中Pref0i＝0W，为防止系统超调过大且充分利用逆变器容量，Ji＝20kg·m2，ki＝100％。
步骤4、根据步骤3中得到的功率指令通过拉格朗日插值预测法得到期望参考功率指令
所述拉格朗日插值预测法包含以下步骤：
步骤4.1，对步骤3中得到的功率指令进行采样处理，采样时间为拉格朗日插值时间T，并按时间先后顺序记录三次采样获得的功率指令：
其中，为从逆变器#i在k-2时刻的功率指令，为从逆变器#i在k-1时刻的功率指令，为从逆变器#i在k时刻的功率指令；
步骤4.2，计算期望参考功率指令其计算公式如下：
其中，t为预测的时间间隔，t取值为t＝2T，T为拉格朗日插值时间，T取值为T＝5Ts，Ts为从逆变器#i的ADC采样时间。本实施例Ts＝0.2ms，T＝1ms，t＝2ms。
步骤5、根据步骤2中得到的输出电压dq轴分量Uodi,Uoqi和网侧电感电流dq轴分量Iodi,Ioqi，通过功率计算方程获得平均有功功率Poi。
所述平均有功功率Poi的计算公式为：
其中Tf为一阶低通滤波器的时间常数，本实例中Tf＝1e-4s。
步骤6、根据步骤2中得到的从逆变器#i在k时刻的直流侧光伏输出电流Ipvi(k)和从逆变器#i在k时刻的直流侧光伏输出电压Upvi(k)、步骤4中得到的期望参考功率指令和步骤5中得到的平均有功功率Poi，通过基于多项式拟合的功率跟踪法获得直流侧参考电压Urefi。
所述基于多项式拟合的功率跟踪法包括两部分，第一部分为通过多项式拟合方法求直流侧参考电压Urefi步长的大小，第二部分为通过功率跟踪法求直流侧参考电压Urefi的方向，具体分为以下5个步骤，其中，6.1属于第一部分，6.2—6.4属于第二部分。
步骤6.1，求电压步长Ustep，其计算公式为：
当时，取Ustep＝0；
当时，取Ustep＝Uthreshold_high；
当时，取
其中，为未限幅电压步长，Ustep为电压步长，Uthreshold_low为低电压步长阈值，Uthreshold_high为高电压步长阈值；
步骤6.2，计算从逆变器#i在k时刻光伏电池板的功率Ppvi(k)，其计算公式为：
Ppvi(k)＝Upvi(k)·Ipvi(k)
其中，k为时刻标志；
步骤6.3，计算从逆变器#i的符号标志flagi，其计算公式为：
flagi＝sign(Ppvi(k)-Ppvi(k-1))×sign(Upvi(k)-Upvi(k-1))
其中，flagi为从逆变器#i的符号标志，Ppvi(k)为从逆变器#i在k时刻光伏电池板的功率，Ppvi(k-1)为从逆变器#i在k-1时刻光伏电池板的功率，Upvi(k)为从逆变器#i在k时刻直流侧光伏输出电压，Upvi(k-1)为从逆变器#i在k-1时刻直流侧光伏输出电压，sign为符号函数数学运算符，其含义如下：
其中，x为自变量；
步骤6.4，求直流侧参考电压Urefi；
当期望参考功率指令大于平均有功功率Poi时，执行Urefi＝Upvi(k)+Ustep×flagi；
当期望参考功率指令小于等于平均有功功率Poi时，如果flagi大于等于零，执行Urefi＝Upvi(k)+Ustep，如果flagi小于零，执行Urefi＝Upvi(k)-Ustep×flagi。
电压步长Ustep的大小根据多项式曲线确定，多项式拟合曲线如图3所示，多项式曲线为二次曲线的一部分，二次曲线的顶点为光伏曲线的最大功率点(UMPP，PMPP)，且二次曲线经过光伏曲线的(Uoc，0)点，本实施例中，Uthreshold_low＝0.3V，Uthreshold_high＝100V，从主逆变器中获得PMPP＝500KW，UMPP＝645.4V，Uoc＝839.2V。
步骤7、根据步骤2中得到的从逆变器#i在k时刻的直流侧光伏输出电压Upvi(k)和步骤6中得到的直流侧参考电压Urefi，经过电压环获得光伏板参考功率Prefi。
所述经过电压环获得光伏板参考功率Prefi的表达式为：
Prefi＝(Upvi(k)-Urefi)Gdc(s)
其中，Gdc(s)为直流侧电压闭环比例积分调节器，其表达式为：
Gdc(s)＝kdcki+kdcpi/s
其中，kdcki为从逆变器#i直流侧电压闭环比例调节器系数，kdcpi为从逆变器#i直流侧电压闭环积分调节器系数，本实例中kdcki＝4550，kdcpi＝1.05e6。
步骤8、根据步骤7中得到的光伏板参考功率Prefi和步骤5中得到的平均有功功率Poi，经过功率外环得到d轴有功功率指令Pdrefi；设置指令电压Edref和从逆变器#i上层无功功率指令Qref0i，并根据步骤2中得到的输出电压d轴分量Uodi，通过一次调压方程获得无功功率指令Qdrefi。
所述功率外环和一次调压方程的表达式分别为：
Pdrefi＝(Prefi-Poi)Gp(s)
其中，ni为从逆变器#i无功功率下垂系数，Gp(s)为功率闭环比例积分调节器，其表达式为：
Gp(s)＝kpki+kppi/s
其中，kpki为从逆变器#i功率闭环比例调节器系数，kppi为从逆变器#i功率闭环积分调节器系数。本实例中取ni＝1/4000，kpki＝0.7，kppi＝1200，Qref0i＝0var。
步骤9、根据步骤2中得到的输出电压d轴分量Uodi和步骤8中得到的d轴有功功率指令Pdrefi、无功功率指令Qdrefi，经过电流计算方法得到电流d轴指令Icdrefi和电流q轴指令Icqrefi。
所述电流计算方法得到电流d轴指令Icdrefi和电流q轴指令Icqrefi的表达式分别为：
步骤10、将步骤9中得到的电流d轴指令Icdrefi和步骤2中得到的桥臂电感电流d轴分量ILdi，经过d轴电流闭环控制方程，得到d轴输出信号Uidi；将步骤9中得到的电流q轴指令Icqrefi与步骤2中得到的桥臂电感电流q轴分量ILqi，经过q轴电流闭环控制方程，得到q轴输出信号Uiqi。
所述d轴电流闭环控制方程和q轴电流闭环控制方程的表达式分别为：
其中，Uidi为d轴输出信号，Uiqi为q轴输出信号，GI(s)为电流闭环比例调节器，其表达式为：
GI(s)＝kIki+kIpi/s
其中，kIki为从逆变器#i电流闭环比例调节器系数,kIpi为从逆变器#i电流闭环积分调节器系数。本实例中取kIki＝40,kIpi＝1。
步骤11、将步骤2中得到的输出电压d轴分量Uodi和输出电压q轴分量Uoqi分别加上步骤10中得到的d轴输出信号Uidi和q轴输出信号Uiqi，得到dq坐标系下的调制波Umdi和Umqi，其表达式分别为：
步骤12、将步骤11中得到的dq坐标系下的调制波Umdi和Umqi经坐标反变换方程得到从逆变器桥臂电压的三相调制波Umai,Umbi,Umci，经PWM调制后作为IGBT电路的驱动信号。
所述坐标反变换方程的表达式为：
Umai＝Umdicosθ+Umqisinθ
其中，θ为d轴和q轴的相位差。
本实施例中发明适用于传统集中式三相光伏逆变器。以下所示为500KW三相光伏逆变器系统加入拉格朗日插值预测法的仿真波形。
主逆变器采用最大功率追踪控制方法，从逆变器采用光伏虚拟同步控制方法，0s时并网运行，0.25s时突加400KW公共阻性负载，在1.00s突降200KW阻性负载。
图4为增加拉格朗日插值预测法的光伏虚拟同步机前后增减负荷电网频率变化波形，加入拉格朗日插值预测法前，突增负载电网频率跌落较快，加入拉格朗日插值预测法后频率跌落变缓，突降负载时，相较未加入拉格朗日插值预测法，电网频率更迟缓恢复，因此加入拉格朗日插值预测法后增加了电网的惯性和阻尼。
图5为加入采用拉格朗日插值预测法的光伏虚拟同步发电机后增减负荷输出有功功率波形，加入采用拉格朗日插值预测法的光伏虚拟同步发电机后，平均有功功率超前于功率指令，在动态过程中有更好的调频效果。
Claims (8)
1.一种基于主从式控制的逆变器光伏虚拟同步控制方法，本控制方法所涉及的逆变器台数为n，n为整数且n≥2，其特征在于，将逆变器分为主逆变器和从逆变器，所有逆变器之间采用通讯方式连接，主逆变器编号为#1，从逆变器编号为#i，i∈[2,n]且i为正整数，主逆变器采用最大功率追踪控制方法，从逆变器采用光伏虚拟同步控制方法；
从逆变器采用的光伏虚拟同步控制方法的步骤如下：
步骤1、采样主逆变器#1光伏电池板最大功率点的功率PMPP及该光伏电池板最大功率点的电压UMPP和开路电压UOC，并经通讯下发给从逆变器；
步骤2、采样从逆变器#i输出相电压Uoai,Uobi,Uoci，并经输出相电压坐标变换方程得到输出电压dq轴分量Uodi,Uoqi；采样从逆变器#i网侧电感电流Ioai,Iobi,Ioci，并经网侧电感电流坐标变换方程得到网侧电感电流dq轴分量Iodi,Ioqi；采样从逆变器#i桥臂电感电流ILai,ILbi,ILci，并经桥臂电感电流坐标变换方程得到桥臂电感电流dq轴分量ILdi,ILqi，其中d轴为有功轴，q轴为无功轴；通过锁相环获得电网频率ωg，采样从逆变器#i在k时刻的直流侧光伏输出电流Ipvi(k)和从逆变器#i在k时刻的直流侧光伏输出电压Upvi(k)；
步骤3、设置电网角频率基准值ωref和从逆变器#i上层有功功率指令Pref0i，并根据步骤1得到的最大功率点的功率PMPP和步骤2得到的电网频率ωg，通过功率叠加方程获得功率指令其表达式为：
其中，Ji为从逆变器#i的虚拟转动惯量，s为拉普拉斯算子，ki为从逆变器#i的预设定平均有功功率输出值占当前最大功率点功率PMPP的百分比，π为圆周率；
步骤4、根据步骤3中得到的功率指令通过拉格朗日插值预测法得到期望参考功率指令
步骤5、根据步骤2中得到的输出电压dq轴分量Uodi,Uoqi和网侧电感电流dq轴分量Iodi,Ioqi，通过功率计算方程获得平均有功功率Poi；
步骤6、根据步骤2中得到的从逆变器#i在k时刻的直流侧光伏输出电流Ipvi(k)和从逆变器#i在k时刻的直流侧光伏输出电压Upvi(k)、步骤4中得到的期望参考功率指令和步骤5中得到的平均有功功率Poi，通过基于多项式拟合的功率跟踪法获得直流侧参考电压Urefi，其基于多项式拟合的功率跟踪法包括以下步骤：
步骤6.1，求电压步长Ustep，其计算公式为：
当时，取Ustep＝0；
当时，取Ustep＝Uthreshold_high；
当时，取
其中，为未限幅电压步长，Ustep为电压步长，Uthreshold_low为低电压步长阈值，Uthreshold_high为高电压步长阈值；
步骤6.2，计算从逆变器#i在k时刻光伏电池板的功率Ppvi(k)，其计算公式为：
Ppvi(k)＝Upvi(k)·Ipvi(k)
其中，k为时刻标志；
步骤6.3，计算从逆变器#i的符号标志flagi，其计算公式为：
flagi＝sign(Ppvi(k)-Ppvi(k-1))×sign(Upvi(k)-Upvi(k-1))
其中，flagi为从逆变器#i的符号标志，Ppvi(k)为从逆变器#i在k时刻光伏电池板的功率，Ppvi(k-1)为从逆变器#i在k-1时刻光伏电池板的功率，Upvi(k)为从逆变器#i在k时刻直流侧光伏输出电压，Upvi(k-1)为从逆变器#i在k-1时刻直流侧光伏输出电压，sign为符号函数数学运算符，其含义如下：
其中，x为自变量；
步骤6.4，求直流侧参考电压Urefi；
当期望参考功率指令大于平均有功功率Poi时，执行Urefi＝Upvi(k)+Ustep×flagi；
当期望参考功率指令小于等于平均有功功率Poi时，如果flagi大于等于零，执行Urefi＝Upvi(k)+Ustep，如果flagi小于零，执行Urefi＝Upvi(k)-Ustep×flagi；
步骤7、根据步骤2中得到的从逆变器#i在k时刻的直流侧光伏输出电压Upvi(k)和步骤6中得到的直流侧参考电压Urefi，经过电压环获得光伏板参考功率Prefi；
步骤8、根据步骤7中得到的光伏板参考功率Prefi和步骤5中得到的平均有功功率Poi，经过功率外环得到d轴有功功率指令Pdrefi；设置指令电压Edref和从逆变器#i上层无功功率指令Qref0i，并根据步骤2中得到的输出电压d轴分量Uodi，通过一次调压方程获得无功功率指令Qdrefi；
步骤9、根据步骤2中得到的输出电压d轴分量Uodi和步骤8中得到的d轴有功功率指令Pdrefi、无功功率指令Qdrefi，经过电流计算方法得到电流d轴指令Icdrefi和电流q轴指令Icqrefi；
步骤10、将步骤9中得到的电流d轴指令Icdrefi和步骤2中得到的桥臂电感电流d轴分量ILdi，经过d轴电流闭环控制方程，得到d轴输出信号Uidi；将步骤9中得到的电流q轴指令Icqrefi与步骤2中得到的桥臂电感电流q轴分量ILqi，经过q轴电流闭环控制方程，得到q轴输出信号Uiqi；
步骤11、将步骤2中得到的输出电压d轴分量Uodi和输出电压q轴分量Uoqi分别加上步骤10中得到的d轴输出信号Uidi和q轴输出信号Uiqi，得到dq坐标系下的调制波Umdi和Umqi，其表达式分别为：
步骤12、将步骤11中得到的dq坐标系下的调制波Umdi和Umqi经坐标反变换方程得到从逆变器桥臂电压的三相调制波Umai,Umbi,Umci，经PWM调制后作为IGBT电路的驱动信号。
2.根据权利要求1所述的基于主从式控制的逆变器光伏虚拟同步控制方法，其特征在于，步骤2所述电网频率ωg的计算公式为：
其中，ω0为公共耦合点电压的额定角频率，Kp_pll为锁相环比例积分调节器的比例调节系数，Ki_pll为锁相环比例积分调节器的积分调节系数，s为拉普拉斯算子。
3.根据权利要求1所述的基于主从式控制的逆变器光伏虚拟同步控制方法，其特征在于，步骤4中所述拉格朗日插值预测法包含以下步骤：
步骤4.1，对步骤3中得到的功率指令进行采样处理，采样时间为拉格朗日插值时间T，并按时间先后顺序记录三次采样获得的功率指令：
其中，为从逆变器#i在k-2时刻的功率指令，为从逆变器#i在k-1时刻的功率指令，为从逆变器#i在k时刻的功率指令；
步骤4.2，计算期望参考功率指令其计算公式如下：
其中，t为预测的时间间隔，t取值为t＝2T，T为拉格朗日插值时间，T取值为T＝5Ts，Ts为从逆变器#i的ADC采样时间。
4.根据权利要求1所述的基于主从式控制的逆变器光伏虚拟同步控制方法，其特征在于，步骤5中的通过功率计算方程获得平均有功功率Poi的计算公式为：
其中Tf为一阶低通滤波器的时间常数，s为拉普拉斯算子。
5.根据权利要求1所述的基于主从式控制的逆变器光伏虚拟同步控制方法，其特征在于，步骤7中的经过电压环获得光伏板参考功率Prefi的表达式为：
Prefi＝(Upvi(k)-Urefi)Gdc(s)
其中，Gdc(s)为直流侧电压闭环比例积分调节器，其表达式为：
Gdc(s)＝kdcki+kdcpi/s
其中，kdcki为从逆变器#i直流侧电压闭环比例调节器系数，kdcpi为从逆变器#i直流侧电压闭环积分调节器系数，s为拉普拉斯算子。
6.根据权利要求1所述的基于主从式控制的逆变器光伏虚拟同步控制方法，其特征在于，步骤8所述功率外环和一次调压方程的表达式分别为：
Pdrefi＝(Prefi-Poi)Gp(s)
其中，ni为从逆变器#i无功功率下垂系数，Gp(s)为功率闭环比例积分调节器，其表达式为：
Gp(s)＝kpki+kppi/s
其中，kpki为从逆变器#i功率闭环比例调节器系数，kppi为从逆变器#i功率闭环积分调节器系数，s为拉普拉斯算子。
7.根据权利要求1所述的基于主从式控制的逆变器光伏虚拟同步控制方法，其特征在于，步骤9所述电流计算方法得到电流d轴指令Icdrefi和电流q轴指令Icqrefi的表达式分别为：
8.根据权利要求1所述的基于主从式控制的逆变器光伏虚拟同步控制方法，其特征在于，步骤10所述d轴电流闭环控制方程和q轴电流闭环控制方程的表达式分别为：
Uidi＝(Icdrefi-ILdi)GI(s)
Uiqi＝(Icqrefi-ILqi)GI(s)
其中，Uidi为d轴输出信号，Uiqi为q轴输出信号，GI(s)为电流闭环比例调节器，其表达式为：
GI(s)＝kIki+kIpi/s
其中，kIki为从逆变器#i电流闭环比例调节器系数,kIpi为从逆变器#i电流闭环积分调节器系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810462242.0A CN108631363B (zh) | 2018-05-15 | 2018-05-15 | 基于主从式控制的逆变器光伏虚拟同步控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810462242.0A CN108631363B (zh) | 2018-05-15 | 2018-05-15 | 基于主从式控制的逆变器光伏虚拟同步控制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108631363A CN108631363A (zh) | 2018-10-09 |
| CN108631363B true CN108631363B (zh) | 2019-10-11 |
Family
ID=63693353
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810462242.0A Active CN108631363B (zh) | 2018-05-15 | 2018-05-15 | 基于主从式控制的逆变器光伏虚拟同步控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108631363B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110266056B (zh) * | 2019-06-25 | 2022-03-25 | 国网湖南省电力有限公司 | 基于光储独立并网的虚拟同步发电控制方法 |
| CN114123218B (zh) * | 2021-11-24 | 2024-02-02 | 国网山东省电力公司电力科学研究院 | 分布式光伏接入末端配电网的过电压分层分区抑制方法 |
| CN115021283B (zh) * | 2022-07-06 | 2024-05-03 | 西安交通大学 | 一种无储能光伏电压型控制方法及系统 |
| CN116905614A (zh) * | 2023-06-13 | 2023-10-20 | 埃梯梯智慧水务科技有限公司 | 基于光伏发电的无负压叠压供水及污水处理方法与系统 |
| CN117175665B (zh) * | 2023-11-02 | 2024-01-26 | 深圳海辰储能科技有限公司 | Pcs多机并联的发波同步系统、方法及相关设备 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106655272B (zh) * | 2017-01-16 | 2018-12-04 | 湖南大学 | 抑制故障瞬时冲击电流型虚拟同步逆变器及其控制方法 |
-
2018
- 2018-05-15 CN CN201810462242.0A patent/CN108631363B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108631363A (zh) | 2018-10-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108631363B (zh) | 基于主从式控制的逆变器光伏虚拟同步控制方法 | |
| Arul et al. | Control strategies for a hybrid renewable energy system: A review | |
| WO2018153222A1 (zh) | 一种基于内模控制的微电网并离网平滑切换控制方法 | |
| CN111817326B (zh) | 一种交流微网孤岛模式下的分散式储能soc控制及整合方法 | |
| Zhang et al. | Active power reserve photovoltaic virtual synchronization control technology | |
| CN108667072B (zh) | 基于有功备用的光伏并网逆变器控制方法 | |
| AU2012222318A1 (en) | Synchronous power controller for a generating system based on static power converters | |
| CN108512247B (zh) | 基于拉格朗日插值预测的光伏虚拟同步控制方法 | |
| CN109888829B (zh) | 基于改进感性下垂控制的光伏微网系统离并网无缝切换系统 | |
| CN108173292B (zh) | 基于功率追踪的光伏虚拟同步控制方法 | |
| CN112436545B (zh) | 孤岛/并网双模式下提升微电网运行稳定性的控制方法 | |
| CN106130062A (zh) | 一种两级式光伏并网发电系统及中间直流侧电压稳定方法 | |
| CN115579944B (zh) | 一种具有自限流保护能力的构网型储能控制系统及方法 | |
| Lee et al. | Inertia-free stand-alone microgrid—Part I: Analysis on synchronized GPS time-based control and operation | |
| Alcala et al. | A three-phase back-to-back converter for reactive power compensation, current harmonic filtering and active power compensation | |
| Zhao et al. | Control Strategies of Microgrid at Micro-source Level and System Level. | |
| Duan et al. | Controller design and stability analysis of electric spring with switchable smart load | |
| Ashabani | Synchronous converter and synchronous-VSC-state of art of universal control strategies for smart grid integration | |
| Wang et al. | Strategy for the seamless mode transfer of an inverter in a master-slave control independent microgrid | |
| Wang et al. | Comparison and Analysis of grid-Forming Control Methods of Wind Turbine Converter under Load Disturbance | |
| Małaczek et al. | Control of low voltage microgrid in autonomous operation mode | |
| Liu et al. | Island operation control strategy for BESS in parallel based on droop control | |
| Yuqiang et al. | Operation control strategy of improved master-slave control for microgrid | |
| Tarraso et al. | Grid forming control for power converters based on an inertial phase locked loop (IPLL) | |
| CN113394825B (zh) | 一种极弱电网下并网逆变器正交功率同步控制方法 |
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 |