CN109245165A - 三相级联h桥光伏逆变器直流侧电压波动抑制方法 - Google Patents
三相级联h桥光伏逆变器直流侧电压波动抑制方法 Download PDFInfo
- Publication number
- CN109245165A CN109245165A CN201811254704.6A CN201811254704A CN109245165A CN 109245165 A CN109245165 A CN 109245165A CN 201811254704 A CN201811254704 A CN 201811254704A CN 109245165 A CN109245165 A CN 109245165A
- Authority
- CN
- China
- Prior art keywords
- phase
- voltage
- grid
- inverter
- current
- 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 22
- 238000002347 injection Methods 0.000 claims abstract description 21
- 239000007924 injection Substances 0.000 claims abstract description 21
- 230000001276 controlling effect Effects 0.000 claims description 12
- 230000001105 regulatory effect Effects 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 9
- 230000001360 synchronised effect Effects 0.000 claims description 9
- 238000005192 partition Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 239000003990 capacitor Substances 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 9
- 230000005611 electricity Effects 0.000 description 6
- 238000005286 illumination Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 2
- 230000002401 inhibitory effect Effects 0.000 description 2
- 230000005764 inhibitory process Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000004886 process control Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000011664 signaling Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 230000003068 static 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/24—Arrangements for preventing or reducing oscillations of power in networks
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Inverter Devices (AREA)
Abstract
本发明公开了一种三相级联H桥光伏逆变器直流侧电压波动抑制方法,其目的是为了解决级联H桥光伏逆变器直流侧电压二倍频波动大的问题。该方法包括如下步骤:(1)总直流侧电压控制,用来实现各H桥单元直流侧电压跟踪其最大功率点电压,并得到系统有功电流指令值;(2)网侧电流解耦控制,能够实现有功电流和无功电流的独立控制,同时产生逆变器原始调制信号;(3)抑制直流侧电压波动控制,通过注入三次谐波修正三相原始调制信号以减小H桥直流侧电压波动。该控制方法能够实现级联H桥光伏逆变器稳定运行的条件下明显减小直流侧电压波动,降低电容成本,且具有算法简单,工程易实现等优点。
Description
技术领域
本发明涉及一种三相级联H桥光伏逆变器直流侧电压波动抑制方法,特别涉及一种基于三次谐波注入的三相级联H桥光伏逆变器H桥直流侧电压波动抑制方法,属于级联型光伏逆变器直流侧电压控制技术领域。
背景技术
光伏并网发电由于提供清洁能源,且环境友好而备受关注。面对如何提高光伏系统效率、降低发电成本等问题,级联H桥多电平逆变器由于其模块化易拓展、系统效率高、并网电流总谐波失真(THD)小等优势而成为研究的热点。
然而,级联型H桥光伏逆变器各H桥模块直流侧相互独立,且三相逆变器结构中每相由单相逆变器组成,交流侧有功功率的二倍频波动会导致各H桥模块直流侧电压出现二倍频纹波。为此,级联型H桥光伏逆变器各H桥模块直流侧均配置较大容量的电容吸收二倍频波动能量实现稳压功能,增大了系统成本和逆变器的体积、重量,同时降低了系统运行的稳定性。因此,对级联型H桥光伏逆变器采取一定的直流侧电压波动抑制控制具有突出的工程意义。
目前,国内外学者针对级联型H桥光伏逆变器直流侧电压波动抑制方法研究较多。如2016年IEEE文献“Cascaded H-Bridge Multilevel PV Topology for Alleviation ofPer-Phase Power Imbalances and Reduction of Second Harmonic Voltage Ripple”Townsend C D,Yu Y,Konstantinou G,《IEEE Transactions on Power Electronics》,2016,31(8),5574-5586(“抑制相间功率不平衡和降低直流侧电压波动的级联H桥多电平光伏拓扑”,《IEEE学报-电力电子期刊》2016年第31卷第8期5574-5586页)提出了一种通过控制直流侧DC/DC变换器降低级联H桥直流侧电压波动的控制策略,但该控制方法基于其提出的带直流侧DC/DC变换器的H桥拓扑,使得硬件实现较复杂,同时电能经2次转换,降低了系统效率。
如文献《电工技术学报》2014年第29卷第2期46-54页“丁明,陈中,程旭东.级联储能变换器直流链纹波电流的抑制策略”提出了基于直流侧有源滤波器的级联H桥变流器拓扑,并利用直流侧有源滤波器抑制H桥直流侧电流和电压波动,但改方法需要在H桥直流侧并联有源滤波器,将承担全部直流侧电压应力,使得硬件实现较复杂,增加了装置体积和成本,同时控制系统的动态响应性能也会受到影响。
综上所述,现有降低级联型H桥光伏逆变器直流侧电压波动的控制策略主要存在如下问题:
(1)现有技术在H桥直流侧串联DC/DC变换器或并联有源滤波器,改变了传统的H桥逆变器拓扑,不利于H桥逆变器模块化结构特点的实现。
(2)现有技术均利用其他电路补偿或吸收H桥直流侧电压波动能量来降低直流侧电压波动,没有针对H桥结构本身的控制来降低其直流侧电压波动。加入的其他电路结构增加了系统的体积和成本,不易于工程实现,同时增加了一级能量转换,降低了系统的效率。
(3)现有技术采用有源滤波器抑制直流侧电压波动,会影响系统的动态响应性能。
发明内容
本发明要解决的问题就是克服上述方案的局限性,针对三相级联H桥光伏逆变器H桥直流侧电压二倍频波动大这一问题,提出一种基于三次谐波注入的三相级联H桥光伏逆变器H桥直流侧电压波动抑制方法。该方法不需要改变传统级联H桥光伏逆变器拓扑,只需改进控制算法,能明显降低各H桥单元直流侧电压波动,从而减小直流侧电容值。
为解决本发明的技术问题,本发明提供了一种一种三相级联H桥光伏逆变器直流侧电压波动抑制方法,所述的三相级联H桥光伏逆变器包括ABC三相,每相由N个带有光伏组件的H桥单元和电感LS组成,本控制方法包括总直流侧电压控制、网侧电流控制和抑制直流侧电压波动控制,步骤如下:
步骤1,总直流侧电压控制
步骤1.1,对三相中的每个H桥单元的直流侧电压采样并依次经过100Hz陷波器和200Hz陷波器滤波,得到三相中每相的N个H桥单元的直流侧电压实际值并记为VPVAi,VPVBi,VPVCi,i=1,2,3...N,ABC分别表示逆变器的三相电路,即A相,B相,C相;采样三相中每相的N个H桥单元的直流侧电流实际值并记为IPVAi,IPVBi,IPVCi,i=1,2,3...N,ABC分别表示逆变器的三相电路,即A相,B相,C相;采样三相的电网电压实际值并记为Vgrid_A,Vgrid_B,Vgrid_C,采样三相的电网电流实际值分别为记为Igrid_A,Igrid_B,Igrid_C,ABC分别表示逆变器的三相电路,即A相,B相,C相;
步骤1.2,通过对三相中每个H桥单元直流侧进行最大功率点跟踪控制,得到三相中每相的N个H桥单元的直流侧电压指令值并记为VPVAi *,VPVBi *,VPVCi *,i=1,2,3...N,ABC分别表示逆变器的三相电路,即A相,B相,C相;
步骤1.3,根据步骤1.1得到的三相中每相的N个H桥单元的直流侧电压实际值VPVAi,VPVBi,VPVCi和步骤1.2得到的三相中每相的N个H桥单元的直流侧电压指令值VPVAi *,VPVBi *,VPVCi *,通过电压调节器,计算得到并网有功电流的指令值Id *,其计算式为:
其中,KVPAi,KVPBi,KVPCi为电压调节器比例系数,i=1,2,3...N;KVIAi,KVIBi,KVICi为电压调节器积分系数,i=1,2,3...N,s为拉普拉斯算子;
步骤2,网侧电流控制
步骤2.1,对步骤1.1中采样的三相电网电压实际值Vgrid_A,Vgrid_B,Vgrid_C进行锁相获得电网电压的相位θ和角频率ω;通过同步旋转坐标变换将步骤1.1中采样的三相电网电压实际值Vgrid_A,Vgrid_B,Vgrid_C转换成旋转坐标系下的电网电压有功分量Vd和电网电压无功分量Vq;通过同步旋转坐标变换将步骤1.1中采样的三相电网电流实际值Igrid_A,Igrid_B,Igrid_C转换成旋转坐标系下的电网电流有功分量Id和电网电流无功分量Iq;
电网电压有功分量Vd和电网电压无功分量Vq的计算式分别为:
电网电流有功分量Id和电网电流无功分量Iq的计算式分别为:
步骤2.2,设逆变器并网无功电流指令值Iq *为0,根据步骤1.3得到的并网有功电流的指令值Id *和步骤2.1得到的电网电流有功分量Id和电网电流无功分量Iq,分别通过有功电流调节器和无功电流调节器,计算得到d轴PI调节值Ed和q轴PI调节值Eq,其计算式分别为:
其中,KiP为有功电流调节器比例系数,KiI为无功电流调节器积分系数,s为拉普拉斯算子;
步骤2.3,根据步骤2.1得到的电网电压有功分量Vd、电网电压无功分量Vq、电网电流有功分量Id、电网电流无功分量Iq、电网电压角频率ω和步骤2.2中得到的d轴PI调节值Ed和q轴PI调节值Eq,计算得到d轴电压控制值Ud和q轴电压控制值Uq,其计算式分别为:
其中,LS为滤波电感;
步骤2.4,根据步骤2.1得到的电网电压的相位θ,将步骤2.3得到的d轴电压控制值Ud和q轴电压控制值Uq经过同步旋转坐标系逆变换得到自然坐标系下的逆变器A相调制波信号VrA、逆变器B相调制波信号VrB、逆变器C相调制波信号VrC,其计算式分别为:
步骤3,抑制直流侧电压波动控制
步骤3.1,对步骤2.4中计算得到的逆变器A相调制波信号VrA进行锁相,得到逆变器A相调制波信号的幅值VrAm和逆变器A相调制波信号的相位θA;
步骤3.2,由步骤3.1中计算得到的逆变器A相调制波信号的幅值VrAm和步骤1.1中得到的A相N个H桥单元的直流侧电压实际值VPVAi计算A相调制波的调制度mkA0,其计算式为:
步骤3.3,由步骤3.2计算得到的A相调制波的调制度mkA0计算三次谐波的调制度mk3。其拟合计算式为:
其中Fj为拟合计算式系数,j=0,1,2,3,4,5,6;
步骤3.4,由步骤3.3计算得到的三次谐波的调制度mk3和步骤1.1中得到的A相H桥单元的直流侧电压实际值VPVAi,计算注入三次谐波的幅值V3rd,其计算式为:
步骤3.5,根据步骤3.1得到的逆变器A相调制波信号的相位θA和步骤3.4得到的注入三次谐波的幅值V3rd计算注入的三次谐波信号u3rd,由此三次谐波信号抑制H桥直流侧电压波动,其计算式为:
u3rd=V3rdcos(3θA)
步骤3.6,由步骤2.4得到的逆变器A相调制波信号VrA、逆变器B相调制波信号VrB、逆变器C相调制波信号VrC和步骤3.5得到的注入的三次谐波信号u3rd,计算注入三次谐波后逆变器三相调制波信号VrA1,VrB1,VrC1,其计算式为:
步骤3.7,计算三相中每相的N个H桥单元的功率值pAi,pBi,pCi和三相中每相总功率值pA,pB,pC,i=1,2,3...N,ABC分别表示逆变器的三相电路,即A相,B相,C相,其计算式分别为:
步骤3.8,根据步骤3.7计算得到的三相中每相的N个H桥单元的功率值pAi,pBi,pCi和三相中每相总功率值pA,pB,pC计算三相中每相的N个H桥单元功率分配系数kAi,kBi,kCi,i=1,2,3...N,ABC分别表示逆变器的三相电路,即A相,B相,C相,其计算式分别为:
其中,pavg为三相平均功率,
步骤3.9,根据步骤1.1得到的三相中每相的N个H桥单元的直流侧电压实际值VPVAi,VPVBi,VPVCi,步骤3.6得到的注入三次谐波后逆变器三相调制波信号VrA1,VrB1,VrC1和步骤3.8得到的三相中每相的N个H桥单元功率分配系数kAi,kBi,kCi,计算三相中每相的N个H桥单元最终的调制信号mAi,mBi,mCi,其计算式分别为:
其中,i=1,2,3...N,ABC分别表示逆变器的三相电路,即A相,B相,C相。
与现有技术相比,本发明公开的一种基于三次谐波注入的三相级联H桥光伏逆变器H桥直流侧电压波动抑制方法,在各H桥单元以最大功率跟踪运行条件下实现了减小直流侧电压波动的效果,其有益效果具体体现在:
1、本发明提出的方法通过在三相调制波信号中注入三次谐波信号,只需改进原有控制算法,实现简单。
2、本发明提出的方法可以明显减小直流侧电压波动。
3、本发明提出的方法在保证逆变器稳定运行条件下基波调制度最大值由传统的1增大至提高了H桥直流侧电压利用率。
附图说明
图1是本发明三相级联H桥光伏逆变器主电路拓扑框图。
图2是本发明三相级联H桥光伏逆变器总控制结构框图。
图3是本发明实施例1三相电网电压实际值锁相控制结构框图。
图4是本发明实施例1逆变器A相调制波信号锁相控制结构框图。
图5是额定温度、额定光照条件下三相级联H桥光伏逆变器每相第一个H桥单元调制信号波形。
图6是额定温度、额定光照条件下三相级联H桥光伏逆变器每相第一个H桥单元直流侧电压波形。
图7是额定温度、光照均为200W/m2条件下三相级联H桥光伏逆变器每相第一个H桥单元调制信号波形。
图8是额定温度、光照均为200W/m2条件下三相级联H桥光伏逆变器每相第一个H桥单元直流侧电压波形。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,下面结合附图及实施例,对本发明作进一步清楚、完整地描述。
图1为本发明实施例三相级联H桥光伏逆变器拓扑结构,如图中所示,所述的三相级联H桥光伏逆变器包括ABC三相,每相由N个带有光伏组件的H桥单元和电感LS组成,具体的,A相包含N个H桥单元为H桥A1、H桥A2…H桥AN,每个H桥单元直流侧分别与光伏电池板PVA1、PVA2…PVAN连接;B相包含N个H桥单元为H桥B1、H桥B2…H桥BN,每个H桥单元直流侧分别与光伏电池板PVB1、PVB2…PVBN连接;C相包含N个H桥单元为H桥C1、H桥C2…H桥CN,每个H桥单元直流侧分别与光伏电池板PVC1、PVC2…PVCN连接;光伏电池板工作条件为在额定温度25℃,额定光照强度1000W/m2下的最大功率点电压为32.5V,每块光伏电池板通过14.1mF电容与每个H桥单元相连,级联系统通过1.5mH电感LS连接到电网,三相的电网电压实际值分别为Vgrid_A、Vgrid_B、Vgrid_C;三相的电网电流实际值分别为Igrid_A、Igrid_B、Igrid_C。
本发明的控制框图如图2所示,包括总直流侧电压控制、网侧电流解耦控制和抑制直流侧电压波动控制三部分。
实施例1
步骤1.1,对三相中的每个H桥单元的直流侧电压采样并依次经过100Hz陷波器和200Hz陷波器滤波,得到三相中每相的N个H桥单元的直流侧电压实际值并记为VPVAi,VPVBi,VPVCi,i=1,2,3...N,ABC分别表示逆变器的三相电路,即A相,B相,C相;采样三相中每相的N个H桥单元的直流侧电流实际值并记为IPVAi,IPVBi,IPVCi,i=1,2,3...N,ABC分别表示逆变器的三相电路,即A相,B相,C相;采样三相的电网电压实际值并记为Vgrid_A,Vgrid_B,Vgrid_C,采样三相的电网电流实际值分别为记为Igrid_A,Igrid_B,Igrid_C,ABC分别表示逆变器的三相电路,即A相,B相,C相。
本实施例中,以每相三个H桥单元为例,每个H桥单元初始时的直流侧电压实际值为VPVA1=VPVA2=VPVA3=VPVB1=VPVB2=VPVB3=VPVC1=VPVC2=VPVC3=32.5V。
步骤1.2,通过对三相中每个H桥单元直流侧进行最大功率点跟踪控制,得到三相中每相的N个H桥单元的直流侧电压指令值并记为VPVAi *,VPVBi *,VPVCi *,i=1,2,3...N,ABC分别表示逆变器的三相电路,即A相,B相,C相。
最大功率点跟踪控制过程为:根据步骤1.1得到的三相中每相的N个H桥单元的直流侧电压实际值VPVAi,VPVBi,VPVCi和三相中每相的N个H桥单元的直流侧电流实际值IPVAi,IPVBi,IPVCi,求得当前每个H桥单元直流侧连接的光伏电池板的输出功率PPVAi,PPVBi,PPVCi,i=1,2,3…N,其中ABC分别表示逆变器的三相电路,即A相,B相,C相,其计算式为:
本实施例中每隔0.3s按照上式计算一次当前每个H桥单元直流侧连接的光伏电池板的输出功率,比较相邻两次的计算值,直到取得每个H桥单元输出功率的最大值,该输出功率对应的直流侧电压实际值即为H桥单元直流侧电压指令值。
本实施例中,初始时刻t=0s时,各H桥单元均工作在额定温度T=25℃,额定光照强度EA1=EA2=EA3=EB1=EB2=EB3=EC1=EC2=EC3=1000W/m2的条件下,得到每个H桥单元的直流侧电压指令值VPVA1 *=VPVA2 *=VPVA3 *=VPVB1 *=VPVB2 *=VPVB3 *=VPVC1 *=VPVC2 *=VPVC3 *=32.5V。
步骤1.3,根据步骤1.1得到的三相中每相的N个H桥单元的直流侧电压实际值VPVAi,VPVBi,VPVCi和步骤1.2得到的三相中每相的N个H桥单元的直流侧电压指令值VPVAi *,VPVBi *,VPVCi *,通过电压调节器,计算得到并网有功电流的指令值Id *,其计算式为:
其中,KVPAi,KVPBi,KVPCi为电压调节器比例系数,i=1,2,3...N;KVIAi,KVIBi,KVICi为电压调节器积分系数,i=1,2,3...N,s为拉普拉斯算子。电压调节器比例系数和电压调节器积分系数按照常规并网逆变器进行设计,本实施例中,KVPAi=KVPBi=KVPCi=5,KVIAi=KVIBi=KVICi=200。
步骤2,网侧电流控制
步骤2.1,对步骤1.1中采样的三相电网电压实际值Vgrid_A,Vgrid_B,Vgrid_C进行锁相获得电网电压的相位θ和角频率ω;通过同步旋转坐标变换将步骤1.1中采样的三相电网电压实际值Vgrid_A,Vgrid_B,Vgrid_C转换成旋转坐标系下的电网电压有功分量Vd和电网电压无功分量Vq;通过同步旋转坐标变换将步骤1.1中采样的三相电网电流实际值Igrid_A,Igrid_B,Igrid_C转换成旋转坐标系下的电网电流有功分量Id和电网电流无功分量Iq。
(1)三相电网电压锁相环锁相过程控制框图如图3所示。
设锁相前电网电压的相位为θ0,T3s2r为三相静止坐标系至旋转坐标系的坐标变换矩阵,计算式为:
将三相电网电压实际值Vgrid_A,Vgrid_B,Vgrid_C通过坐标变换转换为旋转坐标系下的d轴分量Vd0和q轴分量Vq0,其计算式为:
将q轴分量Vq0输入电网电压锁相环调节器得到电网电压锁相环调节器输出ω0,其计算式为:
KfP为电网电压锁相环调节器比例系数,KfI为电网电压锁相环调节器积分系数,s为拉普拉斯算子。本实施例中,KfP=10,KfI=15。
将电网电压锁相环调节器输出ω0与电网电压额定角频率100π叠加后进行积分,再与2π取余,得到电网电压的相位θ0,其计算式为:
其中,mod为取余运算,s为拉普拉斯算子。
重复上述步骤,当电网电压锁相环调节器输出ω0趋于0时,取假设的电网电压的相位θ0的值作为电网电压的相位θ。电网电压的角频率为:
ω=ω0+100π
(2)电网电压有功分量Vd和电网电压无功分量Vq的计算式分别为:
(3)电网电流有功分量Id和电网电流无功分量Iq的计算式分别为:
步骤2.2,设逆变器并网无功电流指令值Iq *为0,根据步骤1.3得到的并网有功电流的指令值Id *和步骤2.1得到的电网电流有功分量Id和电网电流无功分量Iq,分别通过有功电流调节器和无功电流调节器,计算得到d轴PI调节值Ed和q轴PI调节值Eq,其计算式分别为:
其中,KiP为有功电流调节器比例系数,KiI为无功电流调节器积分系数,s为拉普拉斯算子。有功电流调节器比例系数和无功电流调节器积分系数按照常规并网逆变器进行设计,本实施例中,KiP=100,KiI=400。
步骤2.3,根据步骤2.1得到的电网电压有功分量Vd、电网电压无功分量Vq、电网电流有功分量Id、电网电流无功分量Iq、电网电压角频率ω和步骤2.2中得到的d轴PI调节值Ed和q轴PI调节值Eq,计算得到d轴电压控制值Ud和q轴电压控制值Uq,其计算式分别为:
其中,LS为滤波电感。
步骤2.4,根据步骤2.1得到的电网电压的相位θ,将步骤2.3得到的d轴电压控制值Ud和q轴电压控制值Uq经过同步旋转坐标系逆变换得到自然坐标系下的逆变器A相调制波信号VrA、逆变器B相调制波信号VrB、逆变器C相调制波信号VrC,其计算式分别为:
步骤3,抑制直流侧电压波动控制
步骤3.1,对步骤2.4中计算得到的逆变器A相调制波信号VrA进行锁相,得到逆变器A相调制波信号的幅值VrAm和逆变器A相调制波信号的相位θA。
逆变器A相调制波信号VrA的锁相过程控制框图如图4所示。
逆变器A相调制波信号VrA通过全通滤波器后转换为正交信号VrA_β,其计算式为:
其中,s为拉普拉斯算子。
设锁相前逆变器A相调制波信号VrA相位为θA0,图中T2s2r为两相静止坐标系至旋转坐标系的坐标变换矩阵,计算式为:
将逆变器A相调制波信号VrA和正交信号VrA_β通过坐标变换转换为旋转坐标系下的d轴分量Vd_A和q轴分量Vq_A,其计算式为:
将q轴分量Vq_A输入逆变器A相调制波信号锁相环调节器得到锁相环调节器输出ωA0,其计算式为:
KfPA为逆变器A相调制波信号锁相环调节器比例系数,KfIA为逆变器A相调制波信号锁相环调节器积分系数,s为拉普拉斯算子。本实施例中,KfPA=10,KfIA=15。
将逆变器A相调制波信号锁相环调节器输出ωA0与电网电压额定角频率100π叠加后进行积分,再与2π取余,得到逆变器A相调制波信号的相位θA0,其计算式为:
其中,mod为取余运算,s为拉普拉斯算子。
重复上述步骤,当逆变器A相调制波信号锁相环调节器输出ωA0趋于0时,取假设的逆变器A相调制波信号的相位θA0的值作为逆变器A相调制波信号的相位θA。逆变器A相调制波信号的幅值VrAm为:
步骤3.2,由步骤3.1中计算得到的逆变器A相调制波信号的幅值VrAm和步骤1.1中得到的A相N个H桥单元的直流侧电压实际值VPVAi计算A相调制波的调制度mkA0,其计算式为:
步骤3.3,由步骤3.2计算得到的A相调制波的调制度mkA0计算三次谐波的调制度mk3。其拟合计算式为:
其中Fj为拟合计算式系数,j=0,1,2,3,4,5,6,本实施例中取:F0=138.336146884827,F1=-1.010.352587537018,F2=3.035.937825227597,F3=-4830.832644846181,F4=4295.930456145725,F5=-2024.482695947416,F6=395.079987072500。
步骤3.4,由步骤3.3计算得到的三次谐波的调制度mk3和步骤1.1中得到的A相H桥单元的直流侧电压实际值VPVAi,计算注入三次谐波的幅值V3rd,其计算式为:
步骤3.5,根据步骤3.1得到的逆变器A相调制波信号的相位θA和步骤3.4得到的注入三次谐波的幅值V3rd计算注入的三次谐波信号u3rd,由此三次谐波信号抑制H桥直流侧电压波动,其计算式为:
u3rd=V3rdcos(3θA)
本实施例中t=0.15s注入三次谐波。
步骤3.6,由步骤2.4得到的逆变器A相调制波信号VrA、逆变器B相调制波信号VrB、逆变器C相调制波信号VrC和步骤3.5得到的注入的三次谐波信号u3rd,计算注入三次谐波后逆变器三相调制波信号VrA1,VrB1,VrC1,其计算式为:
步骤3.7,计算三相中每相的N个H桥单元的功率值pAi,pBi,pCi和三相中每相总功率值pA,pB,pC,i=1,2,3...N,ABC分别表示逆变器的三相电路,即A相,B相,C相,其计算式分别为:
步骤3.8,根据步骤3.7计算得到的三相中每相的N个H桥单元的功率值pAi,pBi,pCi和三相中每相总功率值pA,pB,pC计算三相中每相的N个H桥单元功率分配系数kAi,kBi,kCi,i=1,2,3...N,ABC分别表示逆变器的三相电路,即A相,B相,C相,其计算式分别为:
其中,pavg为三相平均功率,
步骤3.9,根据步骤1.1得到的三相中每相的N个H桥单元的直流侧电压实际值VPVAi,VPVBi,VPVCi,步骤3.6得到的注入三次谐波后逆变器三相调制波信号VrA1,VrB1,VrC1和步骤3.8得到的三相中每相的N个H桥单元功率分配系数kAi,kBi,kCi,计算三相中每相的N个H桥单元最终的调制信号mAi,mBi,mCi,其计算式分别为:
其中,i=1,2,3...N,ABC分别表示逆变器的三相电路,即A相,B相,C相。
图5为三相级联H桥光伏逆变器每相第一个H桥单元调制信号波形图。t=0.15s时刻注入三次谐波抑制H桥直流侧电压波动。从图中可以看出,注入三次谐波前三相调制信号为对称正弦波,调制度均为0.6,小于0.6495,则注入三次谐波幅值与基频调制波相同以最大程度减小直流侧电压波动。注入三次谐波后三相调制信号为对称鞍型波,且其幅值均小于1。
图6为三相级联H桥光伏逆变器每相第一个H桥单元直流侧电压波形图。t=0.15s时刻注入三次谐波抑制H桥直流侧电压波动。从图中可以看出,抑制前直流侧电压波动最大值为△VPVAmax1=△VPVBmax1=△VPVCmax1=2.0V,抑制后直流侧电压波动最大值为△VPVAmax2=△VPVBmax2=△VPVCmax2=1.35V,三相直流侧电压电压波动均降低32.5%。
实施例2
本实施例中,以每相三个H桥单元为例,每个H桥单元初始时的直流侧电压实际值为VPVA1=VPVA2=VPVA3=VPVB1=VPVB2=VPVB3=VPVC1=VPVC2=VPVC3=30.5V。初始时刻t=0s时,各H桥单元均工作在额定温度T=25℃,光照强度EA1=EA2=EA3=EB1=EB2=EB3=EC1=EC2=EC3=200W/m2的条件下,得到每个H桥单元的直流侧电压指令值VPVA1 *=VPVA2 *=VPVA3 *=VPVB1 *=VPVB2 *=VPVB3 *=VPVC1 *=VPVC2 *=VPVC3 *=30.5V。t=0.15s注入三次谐波。
本实施例中,电压环调节器参数为:KVPAi=KVPBi=KVPCi=5,KVIAi=KVIBi=KVICi=200,电流环调节器参数为:KiP=100,KiI=400。
图7为三相级联H桥光伏逆变器每相第一个H桥单元调制信号波形图。t=0.15s时刻注入三次谐波抑制H桥直流侧电压波动。从图中可以看出,注入三次谐波前三相调制信号为对称正弦波,调制度为0.68,大于0.6495。注入三次谐波后三相调制信号为对称鞍型波,且其幅值均等于1,保证H桥单元不过调制的条件下注入最大的三次谐波以最大程度减小直流侧电压波动。
图8为三相级联H桥光伏逆变器每相第一个H桥单元直流侧电压波形图。t=0.15s时刻注入三次谐波抑制H桥直流侧电压波动。从图中可以看出,抑制前直流侧电压波动最大值为△VPVAmax1=△VPVBmax1=△VPVCmax1=1.56V,抑制后直流侧电压波动最大值为△VPVAmax2=△VPVBmax2=△VPVCmax2=1.05V,三相直流侧电压电压波动均降低32.69%。
Claims (1)
1.一种三相级联H桥光伏逆变器直流侧电压波动抑制方法,所述的三相级联H桥光伏逆变器包括A、B、C三相,每相由N个带有光伏组件的H桥单元和电感LS组成,其特征在于,本控制方法包括总直流侧电压控制、网侧电流控制和抑制直流侧电压波动控制,步骤如下:
步骤1,总直流侧电压控制
步骤1.1,对三相中的每个H桥单元的直流侧电压采样并依次经过100Hz陷波器和200Hz陷波器滤波,得到每相的N个H桥单元的直流侧电压实际值并记为VPVAi,VPVBi,VPVCi,i=1,2,3...N,A、B、C分别表示逆变器的三相电路,即A相,B相,C相;采样三相中每相的N个H桥单元的直流侧电流实际值并记为IPVAi,IPVBi,IPVCi,i=1,2,3...N;采样三相的电网电压实际值并记为Vgrid_A,Vgrid_B,Vgrid_C,采样三相的电网电流实际值分别为记为Igrid_A,Igrid_B,Igrid_C;
步骤1.2,通过对三相中每个H桥单元直流侧进行最大功率点跟踪控制,得到三相中每相的N个H桥单元的直流侧电压指令值并记为VPVAi *,VPVBi *,VPVCi *;
步骤1.3,根据步骤1.1得到的三相中每相的N个H桥单元的直流侧电压实际值VPVAi,VPVBi,VPVCi和步骤1.2得到的三相中每相的N个H桥单元的直流侧电压指令值VPVAi *,VPVBi *,VPVCi *,通过电压调节器,计算得到并网有功电流的指令值Id *,其计算式为:
其中,KVPAi,KVPBi,KVPCi为电压调节器比例系数,i=1,2,3...N;KVIAi,KVIBi,KVICi为电压调节器积分系数,i=1,2,3...N,s为拉普拉斯算子;
步骤2,网侧电流控制
步骤2.1,对步骤1.1中采样的三相电网电压实际值Vgrid_A,Vgrid_B,Vgrid_C进行锁相获得电网电压的相位θ和角频率ω;通过同步旋转坐标变换将步骤1.1中采样的三相电网电压实际值Vgrid_A,Vgrid_B,Vgrid_C转换成旋转坐标系下的电网电压有功分量Vd和电网电压无功分量Vq;通过同步旋转坐标变换将步骤1.1中采样的三相电网电流实际值Igrid_A,Igrid_B,Igrid_C转换成旋转坐标系下的电网电流有功分量Id和电网电流无功分量Iq;
电网电压有功分量Vd和电网电压无功分量Vq的计算式分别为:
电网电流有功分量Id和电网电流无功分量Iq的计算式分别为:
步骤2.2,设逆变器并网无功电流指令值Iq *为0,根据步骤1.3得到的并网有功电流的指令值Id *和步骤2.1得到的电网电流有功分量Id和电网电流无功分量Iq,分别通过有功电流调节器和无功电流调节器,计算得到d轴PI调节值Ed和q轴PI调节值Eq,其计算式分别为:
其中,KiP为有功电流调节器比例系数,KiI为无功电流调节器积分系数,s为拉普拉斯算子;
步骤2.3,根据步骤2.1得到的电网电压有功分量Vd、电网电压无功分量Vq、电网电流有功分量Id、电网电流无功分量Iq、电网电压角频率ω和步骤2.2中得到的d轴PI调节值Ed和q轴PI调节值Eq,计算得到d轴电压控制值Ud和q轴电压控制值Uq,其计算式分别为:
其中,LS为滤波电感;
步骤2.4,根据步骤2.1得到的电网电压的相位θ,将步骤2.3得到的d轴电压控制值Ud和q轴电压控制值Uq经过同步旋转坐标系逆变换得到自然坐标系下的逆变器A相调制波信号VrA、逆变器B相调制波信号VrB、逆变器C相调制波信号VrC,其计算式分别为:
步骤3,抑制直流侧电压波动控制
步骤3.1,对步骤2.4中计算得到的逆变器A相调制波信号VrA进行锁相,得到逆变器A相调制波信号的幅值VrAm和逆变器A相调制波信号的相位θA;
步骤3.2,由步骤3.1中计算得到的逆变器A相调制波信号的幅值VrAm和步骤1.1中得到的A相N个H桥单元的直流侧电压实际值VPVAi计算A相调制波的调制度mkA0,其计算式为:
步骤3.3,由步骤3.2计算得到的A相调制波的调制度mkA0计算三次谐波的调制度mk3,其拟合计算式为:
其中Fj为拟合计算式系数,j=0,1,2,3,4,5,6;
步骤3.4,由步骤3.3计算得到的三次谐波的调制度mk3和步骤1.1中得到的A相H桥单元的直流侧电压实际值VPVAi,计算注入三次谐波的幅值V3rd,其计算式为:
步骤3.5,根据步骤3.1得到的逆变器A相调制波信号的相位θA和步骤3.4得到的注入三次谐波的幅值V3rd计算注入的三次谐波信号u3rd,由此三次谐波信号抑制H桥直流侧电压波动,其计算式为:
u3rd=V3rdcos(3θA)
步骤3.6,由步骤2.4得到的逆变器A相调制波信号VrA、逆变器B相调制波信号VrB、逆变器C相调制波信号VrC和步骤3.5得到的注入的三次谐波信号u3rd,计算注入三次谐波后逆变器三相调制波信号VrA1,VrB1,VrC1,其计算式为:
步骤3.7,计算三相中每相的N个H桥单元的功率值pAi,pBi,pCi和三相中每相总功率值pA,pB,pC,i=1,2,3...N,A、B、C分别表示逆变器的三相电路,即A相,B相,C相,其计算式分别为:
步骤3.8,根据步骤3.7计算得到的三相中每相的N个H桥单元的功率值pAi,pBi,pCi和三相中每相总功率值pA,pB,pC计算三相中每相的N个H桥单元功率分配系数kAi,kBi,kCi,i=1,2,3...N,A、B、C分别表示逆变器的三相电路,即A相,B相,C相,其计算式分别为:
其中,pavg为三相平均功率,
步骤3.9,根据步骤1.1得到的三相中每相的N个H桥单元的直流侧电压实际值VPVAi,VPVBi,VPVCi,步骤3.6得到的注入三次谐波后逆变器三相调制波信号VrA1,VrB1,VrC1和步骤3.8得到的三相中每相的N个H桥单元功率分配系数kAi,kBi,kCi,计算三相中每相的N个H桥单元最终的调制信号mAi,mBi,mCi,其计算式分别为:
其中,i=1,2,3...N,A、B、C分别表示逆变器的三相电路,即A相,B相,C相。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811254704.6A CN109245165B (zh) | 2018-10-26 | 2018-10-26 | 三相级联h桥光伏逆变器直流侧电压波动抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811254704.6A CN109245165B (zh) | 2018-10-26 | 2018-10-26 | 三相级联h桥光伏逆变器直流侧电压波动抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109245165A true CN109245165A (zh) | 2019-01-18 |
CN109245165B CN109245165B (zh) | 2020-08-28 |
Family
ID=65082097
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811254704.6A Active CN109245165B (zh) | 2018-10-26 | 2018-10-26 | 三相级联h桥光伏逆变器直流侧电压波动抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109245165B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110086196A (zh) * | 2019-04-01 | 2019-08-02 | 合肥工业大学 | 一种弱电网下单相级联h桥并网逆变器控制方法 |
CN110098637A (zh) * | 2019-05-23 | 2019-08-06 | 合肥工业大学 | 级联h桥光伏并网逆变器双模式平滑切换控制方法 |
CN111416382A (zh) * | 2020-03-09 | 2020-07-14 | 上海发电设备成套设计研究院有限责任公司 | 一种两级式三相级联光伏并网逆变器的控制方法 |
CN112564171A (zh) * | 2020-12-14 | 2021-03-26 | 青岛大学 | 一种级联h桥光伏并网逆变器调制波的配置策略 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130114311A1 (en) * | 2009-06-29 | 2013-05-09 | Ideal Power Converters, Inc. | Power Transfer Devices, Methods, and Systems with Crowbar Switch Shunting Energy-Transfer Reactance |
WO2017072593A1 (en) * | 2015-10-26 | 2017-05-04 | Rakesh Goel | System and method for improved performance and synchronization of a grid tie inverter |
CN107026474A (zh) * | 2017-05-12 | 2017-08-08 | 合肥工业大学 | 减小级联h桥逆变器直流电压波动的功率均衡控制方法 |
CN107733269A (zh) * | 2017-10-12 | 2018-02-23 | 合肥工业大学 | 扩大级联h桥型光伏逆变器运行范围的方波补偿控制方法 |
-
2018
- 2018-10-26 CN CN201811254704.6A patent/CN109245165B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130114311A1 (en) * | 2009-06-29 | 2013-05-09 | Ideal Power Converters, Inc. | Power Transfer Devices, Methods, and Systems with Crowbar Switch Shunting Energy-Transfer Reactance |
WO2017072593A1 (en) * | 2015-10-26 | 2017-05-04 | Rakesh Goel | System and method for improved performance and synchronization of a grid tie inverter |
CN107026474A (zh) * | 2017-05-12 | 2017-08-08 | 合肥工业大学 | 减小级联h桥逆变器直流电压波动的功率均衡控制方法 |
CN107733269A (zh) * | 2017-10-12 | 2018-02-23 | 合肥工业大学 | 扩大级联h桥型光伏逆变器运行范围的方波补偿控制方法 |
Non-Patent Citations (1)
Title |
---|
Y. YU, G. KONSTANTINOU, B. HREDZAK, AND V. AGELIDIS: "Power Balance of Cascaded H-Bridge Multilevel", 《IEEE TRANSACTIONS ON POWER ELECTRONICS》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110086196A (zh) * | 2019-04-01 | 2019-08-02 | 合肥工业大学 | 一种弱电网下单相级联h桥并网逆变器控制方法 |
CN110086196B (zh) * | 2019-04-01 | 2020-08-28 | 合肥工业大学 | 一种弱电网下单相级联h桥并网逆变器控制方法 |
CN110098637A (zh) * | 2019-05-23 | 2019-08-06 | 合肥工业大学 | 级联h桥光伏并网逆变器双模式平滑切换控制方法 |
CN110098637B (zh) * | 2019-05-23 | 2020-11-24 | 合肥工业大学 | 级联h桥光伏并网逆变器双模式平滑切换控制方法 |
CN111416382A (zh) * | 2020-03-09 | 2020-07-14 | 上海发电设备成套设计研究院有限责任公司 | 一种两级式三相级联光伏并网逆变器的控制方法 |
CN111416382B (zh) * | 2020-03-09 | 2021-07-09 | 上海发电设备成套设计研究院有限责任公司 | 一种两级式三相级联光伏并网逆变器的控制方法 |
CN112564171A (zh) * | 2020-12-14 | 2021-03-26 | 青岛大学 | 一种级联h桥光伏并网逆变器调制波的配置策略 |
Also Published As
Publication number | Publication date |
---|---|
CN109245165B (zh) | 2020-08-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107528491B (zh) | 一种级联h桥多电平变换器及其控制方法 | |
CN109245165A (zh) | 三相级联h桥光伏逆变器直流侧电压波动抑制方法 | |
CN111769591B (zh) | 基于双分裂变压器的多逆变器系统双模式组合控制方法 | |
Hu et al. | Impedance characteristic analysis and stability improvement method for DFIG system within PLL bandwidth based on different reference frames | |
CN104362668A (zh) | 一种电压不平衡/谐波畸变时双馈风力发电机的控制方法 | |
Rastogi et al. | Performance investigation of two-level reduced-switch D-STATCOM in grid-tied solar-PV array with stepped P&O MPPT algorithm and modified SRF strategy | |
CN103475033A (zh) | 无锁相环节三相lcl型并网逆变器电流控制方法及系统 | |
Ni et al. | Fault-tolerant operation of DFIG-WT with four-switch three-phase grid-side converter by using simplified SVPWM technique and compensation schemes | |
CN110086196A (zh) | 一种弱电网下单相级联h桥并网逆变器控制方法 | |
CN113629763B (zh) | 非理想电网下中压直挂储能变流器电流控制方法及系统 | |
CN113991715A (zh) | 非理想电网下中压直挂不对称混合储能系统控制方法 | |
Sundarabalan et al. | Unified Power Quality Conditioner for Enhancement of Power Quality and Hybrid Power Generation Injection to Grid | |
CN104426160A (zh) | 光伏并网逆变器兼容低电压穿越的正负序解耦控制方法 | |
Aiswarya et al. | A comparative performance analysis of PV grid interface STATCOM control algorithms | |
Seo et al. | Harmonics and reactive power compensation method by grid-connected Photovoltaic generation system | |
CN110071524A (zh) | 单相级联h桥光伏逆变器虚拟同步控制方法 | |
Bhanu et al. | SVPWM: torque level controlling of wind turbine system using fuzzy and ABC-DQ transformation | |
CN110071523A (zh) | 基于单元有功备用的级联型逆变器虚拟同步控制方法 | |
CN113013898A (zh) | 基于远端电网锁相的并网逆变器次同步振荡抑制方法 | |
Mora et al. | Active power oscillation elimination in 4-leg grid-connected converters under unbalanced network conditions | |
Surasmi et al. | Grid integrated solar photovoltaic converter control for power quality enhancement | |
CN110912130A (zh) | 一种双交流母线并网变换器的电路结构及其谐波补偿方法 | |
Mao et al. | Control of a cascaded STATCOM with battery energy storage system under unbalanced and distorted grid voltage conditions | |
Reddy et al. | Hybrid renewable energy sources based four leg inverter for power quality improvement | |
Preetha et al. | Power quality enhancement using Distributed Generation inverters with active power control |
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 |