CN101521481A - 一种双馈异步风力发电系统不对称协调直接功率控制方法 - Google Patents

一种双馈异步风力发电系统不对称协调直接功率控制方法 Download PDF

Info

Publication number
CN101521481A
CN101521481A CN200910097499A CN200910097499A CN101521481A CN 101521481 A CN101521481 A CN 101521481A CN 200910097499 A CN200910097499 A CN 200910097499A CN 200910097499 A CN200910097499 A CN 200910097499A CN 101521481 A CN101521481 A CN 101521481A
Authority
CN
China
Prior art keywords
signal
voltage
stator
phase
gsc
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
Application number
CN200910097499A
Other languages
English (en)
Other versions
CN101521481B (zh
Inventor
周鹏
贺益康
章玮
孙丹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN2009100974991A priority Critical patent/CN101521481B/zh
Publication of CN101521481A publication Critical patent/CN101521481A/zh
Application granted granted Critical
Publication of CN101521481B publication Critical patent/CN101521481B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Eletrric Generators (AREA)

Abstract

本发明公开了一种双馈异步风力发电系统不对称协调直接功率控制方法。通过一个比例调节器和两个谐振调节器组成的比例谐振调节器调节双馈异步发电机的定子电磁、无功功率以及网侧变换器的有功、无功功率,从而实现在无需进行任何正、负序电量分解的情况下对双馈异步发电机和网侧变换器的运行状况进行有效控制。采用本发明可实现双馈异步风力发电系统在电网不对称故障情况下的整体最优化控制效果,消除电网电压不对称引起的双馈异步风力发电系统总有功、无功功率、电磁转矩、直流母线电压波动,有效提高该系统在电网电压故障条件下的运行控制性能,确保供电电能质量和电力系统的运行稳定性及安全。

Description

一种双馈异步风力发电系统不对称协调直接功率控制方法
技术领域
本发明涉及风力发电机的控制方法,特别是一种双馈异步风力发电系统不对称协调直接功率控制方法。
背景技术
目前我国的风电技术大多停留在理想电网条件下的运行控制,由于实际电网经常有各类对称、不对称故障发生,因此必须开展电网故障下的运行控制研究并提出相应控制技术。现代大型风力发电系统主要有双馈异步发电机(DFIG)和永磁同步发电机两种类型,其中DFIG系统是当前的主流机型。近年来国际上DFIG机组控制技术的研究多集中在电网对称故障下的低电压穿越运行控制,但电网不对称故障更为频繁、几率更大,因此,DFIG故障运行研究已从对称故障向不对称故障延伸。在DFIG控制系统中若未曾考虑电网电压的不对称,很小的不对称电压将造成定子电流的高度不对称,致使定子绕组产生不对称发热,发电机产生转矩脉动,导致输入电网的有功、无功功率发生振荡。若风电机组相对电网容量足够大,这种缺乏不对称电网电压控制能力的风电机组不得不从电网中解列,以防引发后续的更大电网故障。但从电网安全角度又要求风电机组能承受最大达2%的稳态和相对较大瞬态不对称电压而不退出电网,这就要求风电机组能实现电网电压不对称故障穿越运行。目前,国内、外已经兴起了对这种不对称电网电压条件下DFIG发电机及相关励磁变频器控制方法与实施方案的研究。检索到DFIG风电机组不对称电网条件下运行控制的相关专利和研究文章有:
I.胡家兵,贺益康等.不对称电网电压条件下双馈异步风力发电系统的建模与控制.电力系统自动化,2007,31(14):47-56.
II.L Xu,and Y.Wang,“Dynamic Modeling and Control of DFIG Based WindTurbines under Unbalanced Network Conditions,”IEEE Trans.Power System,Vol.22,No.1,pp.314-323,Feb.2007.
III.Lie Xu,“Coordinated control of DFIG′s rotor and grid side convertersduring network unbalance,”IEEE Trans.Power Electronics,vol.23,no.3,pp.1041-1049,May.2008.
IV.Abo-Khalil,A.G.,Dong-Choon Lee,and Jeong-Ik Jang,“Control ofback-to-back PWM converters for DFIG wind turbine systems under unbalancedgrid voltage,”in Proc.of ISIE 2007,pp.2637-2642,Jun.2007.
V.CARTWRIGHT P,XU L.System controller for e.g.wind powered doublyfed induction generator attached to wind turbine,has grid imbalance detector whichcontrols current to cancel imbalance in grid served by generators[Patent].PatentNumber:GB2420456-A.Date:20060524.Application Number:GB025662.Date:20041123.
不对称电网电压条件下,上述文献提出的方法都是基于对称分量理论的矢量控制方法。这些方法的核心思想是将转子电流分解为正序和负序分量,通过分别控制转子电流的正序和负序分量来间接控制DFIG 3的输出功率和电磁转矩,其原理可用图1来说明。DFIG 3的转子侧变换器1采用双比例积分调节器19分别对转子正、负序电流作独立控制;但为实现对正、负序转子电流的分别调节,必须首先获得反馈转子电流的正、负序分量,其处理过程是:利用二个三相电流霍尔传感器4分别采集DFIG3的三相定子电流Isabc和转子电流信号Irabc,利用一个三相电压霍尔传感器5采集三相定子电压信号Usabc;采集得到的三相定、转子电流信号Isabc和Irabc,定子电压信号Usabc分别经过静止三相到二相坐标变换模块6,得到包含正、负序分量的定子电压综合矢量Usαβ和定、转子电流综合矢量Isαβ和Irαβ;其中Usαβ、Isαβ分别通过正、反转同步速旋转坐标变换模块13,14,得到在电网电压不对称条件下正、反转同步速旋转坐标系中含有直流量与两倍频2ωs交流量之和的电压综合矢量
Figure A200910097499D00061
;Irαβ分别通过正、反转转差速旋转坐标变换模块11,12,得到转子电流综合矢量
Figure A200910097499D00062
Figure A200910097499D0006151032QIETU
;然后采用2ωs频率陷波器16(或低通滤波器、1/4电网电压基波周期延时等方法)来滤除
Figure A200910097499D0006151632QIETU
Figure A200910097499D0006151636QIETU
中2ωs频率的交流成分,从而获得其正、负序分量
Figure A200910097499D00064
(直流量)。在此基础上,通过定子磁链观测器18获取转子电流指令值计算模块17和反馈补偿解耦模块15所需的定子磁链分量
Figure A200910097499D00065
,根据电网电压不对称条件下DFIG3不同的控制目标由转子电流指令值计算模块17计算获得转子电流指令
Figure A200910097499D00066
Figure A200910097499D00067
并与转子电流反馈信号
Figure A200910097499D00068
比较获得误差信号,然后分别在正、反转同步速旋转坐标系中采用比例积分器19对误差信号作比例-积分调节,调节得到信号经反馈补偿解耦模块15补偿解耦获得正、反转同步速旋转坐标系中的正、负序转子电压参考值
Figure A200910097499D00069
分别通过不同的旋转坐标变换模块9,10转换得到定子坐标系中的正、负序转子电压参考值,并相加后得到空间矢量脉宽调制SVPWM模块2的参考信号
Figure A200910097499D000610
,经过SVPWM模块2调制获得转子侧变换器1的开关信号以控制DFIG 3运行,实现不对称电网电压条件下DFIG 3正、负序转子电流在正、反转同步旋转坐标系中的独立闭环控制,达到所要求的控制目标。此外,该方法采用软件锁相环8电路对电网电压的频率和相位进行检测,转子位置采用编码器7测定,再通过微分器35计算速度,为定、转子电压、电流采集信号实现正、反转旋转坐标变换提供依据。
由上述分析过程可见,电网电压不对称条件下传统DFIG 3控制方法的实质是将不对称系统分解成正、负序对称分量系统后,再分别在正、反转同步旋转坐标系中实现d、q轴解耦控制。虽然转子正、负序电流在正、反转同步旋转坐标系中各自表现为直流量,分别采用两个PI调节器即可实现无静差独立跟踪控制,但控制实施的前提是已实现对采集转子电流的正、负序分离。图1所示传统控制方法中正、负序分离普遍采用了2ωs频率陷波器16(或低通滤波器、1/4电网电压基波周期延时等方法),分离中除引入延时外,控制系统带宽将受到影响,会造成动态跟踪误差,动态控制效果不理想。更有甚者,该电路无法区分电网电压是平衡还是不对称,如果DFIG 3运行在严格电网电压平衡状态下,控制系统仍将采用陷波器来分离转子变量,这将给系统正常控制带来了不必要的延时,严重影响了系统的动态控制性能。此外,由于传统DFIG 3控制方法仅有转子电流正序d、q轴分量和负序d、q轴分量四个可控量,因此只能在控制定子有功、无功功率平均值之外,再有选择地消除定子有功或者无功功率中的二倍频振荡,而不能同时消除定子有功、无功功率以及电磁转矩中的二倍频振荡,更无法顾及定子有功、无功功率以及电磁转矩中的基频振荡。
文献III和IV中提出利用GSC 28来补偿DFIG 3定子有功功率中的二倍频振荡分量,但是GSC 28的控制中同样需要采用正、负序分离和双比例积分调节器19,而且没有考虑GSC 28输出有功、无功功率中的基频频振荡,使得控制系统结构复杂、动态响应慢、不能达到预期的控制效果。
综上所述,亟需探索一种无需正负序分解、又能消除电网电压不对称引起的双馈异步风力发电系统总有功、无功功率、电磁转矩、直流母线电压波动的控制方法,以适应电网平衡与不对称条件下DFIG风电机组的运行控制。
发明内容
本发明的目的是提供一种双馈异步风力发电系统不对称协调直接功率控制方法,该方法无需进行任何正、负序分解,免除了正、负序分解操作而引入控制延时,并且能消除电网电压不对称引起的双馈异步风力发电系统总有功、无功功率、电磁转矩、直流母线电压波动,从而有效提高双馈异步风电系统在电网电压故障条件下的运行控制性能,确保供电电能质量和电力系统的运行稳定性及安全。
本发明的技术解决方案,一种双馈异步风力发电系统不对称协调直接功率控制方法,包括双馈异步风力发电系统中转子侧变换器RSC的不对称直接功率控制和网侧变换器GSC的不对称直接功率控制;通过对RSC的不对称直接功率控制消除双馈异步发电机DFIG定子电磁功率和无功功率中的基频ωs和两倍频2ωs振荡,保持DFIG的电磁转矩恒定;通过对GSC的不对称直接功率控制消除直流母线电压和GSC输出无功功率中的基频ωs和两倍频2ωs振荡,保持双馈异步风力发电系统向电网输出的总有功、无功功率恒定;
所述的RSC不对称直接功率控制,包括以下步骤;
(i)利用三相电压霍尔传感器采集双馈异步感应发电机DFIG三相定子电压信号Usabc;利用第一三相电流霍尔传感器采集三相定子电流Isabc,利用第二三相电流霍尔传感器采集转子三相电流信号Irabc
(ii)采集得到的三相定子电压信号Usabc经过不对称锁相环检测得到电网或者定子电压中正序分量的角频率ωs和相位θs;与此同时采用编码器检测DFIG的转子位置θr,再经过微分器计算转速ωr;并通过减法器计算得到滑差角频率ωslip=ωsr
(iii)将采集得到的三相定子电压信号Usabc,三相定子电流信号Isabc分别经过第一静止三相到二相坐标变换模块,得到定子静止坐标系中包含正、负序分量的定子电压综合矢量Usαβ,以及包含正、负序和零序分量的定子电流综合矢量Isαβ
(iv)将得到的定子静止坐标系中定子电压综合矢量Usαβ、定子电流综合矢量Isαβ经过定子有功、无功功率计算模块得到定子输出的瞬时有功功率信号Ps和无功功率信号Qs
(v)将采集得到的三相转子电流信号Irabc通过正转转子速旋转坐标变换模块,得到转子速旋转坐标系中含正、负序和零序分量的转子电流综合矢量Irαβ
(vi)将得到的定子电流综合矢量Isαβ和转子电流综合矢量Irαβ经过定子电磁功率计算模块得到定子电磁功率信号Pes
(vii)将定子电磁功率信号Pes和无功功率信号Qs与其参考电磁功率信号
Figure A200910097499D00081
和无功功率信号
Figure A200910097499D00082
经过减法器计算得到电磁功率误差信号ΔPes和无功功率误差信号ΔQs
(viii)将得到的电磁功率误差信号ΔPes和无功功率误差信号ΔQs通过第一比例谐振调节器调节;调节后的输出信号
Figure A200910097499D00091
经过第一反馈补偿解耦模块完成转子速旋转坐标系中交-直轴间的交叉解耦和动态反馈补偿,获取转子速旋转坐标系中的转子电压参考值
Figure A200910097499D00092
(ix)转子电压参考值
Figure A200910097499D00093
经过转差速旋转坐标变换模块,获得定子坐标系中转子电压参考值
Figure A200910097499D00094
作为空间矢量第一脉宽调制模块调制所需的转子电压参考信号,该信号经过空间矢量脉宽调制后获得控制DFIG运行的转子侧变换器的开关信号Sa1,Sb1,Sc1
所述的GSC不对称直接功率控制,包括以下步骤:
(i)利用单相电压霍尔传感器采集直流母线电压信号Vdc;利用第三三相电流霍尔传感器采集GSC输出的流过滤波电感的三相电流信号Igabc
(ii)将采集得到的GSC三相输出电流信号Igabc经过第二静止三相到二相坐标变换模块,得到静止坐标系中包含正、负序和零序分量的GSC电流综合矢量Igαβ
(iii)将得到的静止坐标系中定子电压综合矢量Usαβ、GSC电流综合矢量Igαβ经过GSC有功、无功功率计算模块得到GSC输出的瞬时有功功率信号Pg和无功功率信号Qg
(iv)将直流母线电压参考信号
Figure A200910097499D0009152156QIETU
与采集得到的直流母线电压信号Vdc经过减法器计算得到直流母线电压误差信号,利用第一比例积分调节器对得到的误差信号作比例-积分调节,调节输出信号经过GSC有功功率参考信号补偿模块得到GSC有功功率参考信号
Figure A200910097499D00095
(v)将GSC输出的有功功率信号Pg和无功功率信号Qg与其参考有功功率信号和无功功率信号
Figure A200910097499D00097
经过减法器计算得到GSC输出有功误差信号ΔPg和无功功率误差信号ΔQg
(vi)将得到的有功误差信号ΔPg和无功功率误差信号ΔQg通过第二比例谐振调节器调节;调节后的输出信号经过第二反馈补偿解耦模块完成同步速旋转坐标系中交-直轴间的交叉解耦和动态反馈补偿,获取同步速旋转坐标系中的GSC输出电压参考值
Figure A200910097499D00098
(vii)GSC输出电压参考值
Figure A200910097499D00099
经过同步速旋转坐标变换模块,获得的静止坐标系中GSC输出电压参考值
Figure A200910097499D000910
作为第二脉宽调制模块调制所需的GSC输出电压参考信号,该信号经过空间矢量脉宽调制后获得控制GSC运行的开关信号Sa2,Sb2,Sc2
上述的RSC不对称直接功率控制中所述的采用不对称锁相环检测电网或者定子电压中正序分量的角频率ωs和相位θs,步骤如下:
(i)利用锁相环的反馈相位信号
Figure A200910097499D0010152247QIETU
对三相定子电压信号Usabc进行正转坐标变换,得到正转坐标系中含有直流量与两倍频2ωs交流量之和的电压综合矢量Usdq
(ii)将得到的正转坐标系中电压综合矢量Usdq的q轴分量Usq经过第二比例积分调节器得到三相定子电压正序分量的频率ωs
(iii)将得到的频率信号ωs经过积分器积分得到电压正序分量的相位信号θs
(iv)Usq经过第一两倍频2ωs谐振调节器调节后的输出信号与电压正序分量的相位信号θs相加,得到锁相环的反馈相位信号
Figure A200910097499D0010152247QIETU
上述的RSC不对称直接功率控制和GSC不对称直接功率控制中所述的第一比例谐振调节器和第二比例谐振调节器分别包括一个比例控制器和两个角频率分别为基频ωs和两倍频2ωs的谐振调节器,其中比例控制器调节功率的平均值,角频率为ωs的谐振调节器可抑制功率中的基频ωs振荡,角频率为2ωs的第二倍频2ωs谐振调节器抑制功率中的两倍频2ωs振荡。
本发明提出的控制方法比传统的正、负序双d、q解耦控制方法大为简化,消除了电流内环控制环节,可直接对输出功率实施控制,由于无需进行正、负序分解,不会因此引入分解延时,两个谐振调节器可抑制输出功率的振荡,从而获得稳定的输出、较小的稳态误差以及较好的动态响应特性。而且,通过对RSC和GSC的协调,可以消除电网电压不对称引起的双馈异步风力发电系统总有功、无功功率、电磁转矩、直流母线电压波动,实现该系统在电网不对称情况下的整体最优化控制效果。采用该方法可在不对称电网电压条件下实现发电系统的增强控制目标,有效提高该类风电系统电网故障下的不间断运行能力。
本发明方法适用于除DFIG风电系统之外的其他采用高频开关自关断器件构成的各类形式PWM控制的三相或单相逆变装置在平衡与不对称电网电压条件下的有效控制,如太阳能、燃料电池发电系统的并网逆变装置,柔性输电系统的电力电子逆变装置即以电力调速传动中的双馈电动机或发电机变流装置的有效控制。
附图说明
图1是不对称电网电压条件下双馈异步风力发电机传统控制方法的原理图。
图2是双馈异步风力发电机组的原理图。
图3是本发明的双馈异步风力发电系统不对称协调直接功率控制方法的原理图。
图4是本发明的不对称软件锁相环的原理图。
图5是本发明的比例谐振调节器的原理图。
图6为电网电压瞬态不对称条件下的仿真效果图,图(A)为未采用本发明方法,图(B)采用本发明方法。图(A)和图(B)中,(a)DFIG定子三相电压(pu);(b)DFIG定子三相电流(pu);(c)DFIG转子三相电流(pu);(d)DFIG定子输出有功功率(pu);(e)GSC输出有功功率(pu);(f)系统输出总有功功率(pu);(g)DFIG定子输出无功功率(pu);(h)GSC输出无功功率(pu);(i)系统输出总无功功率(pu);(j)DFIG电磁转矩(pu);(k)直流母线电压(V);
具体实施方式
下面结合附图对本发明进一步说明。
图2是双馈异步风力发电机组的原理图,包括风轮机、齿轮箱、双馈异步发电机DFIG 3、转子侧变换器RSC 1、网侧侧变换器GSC 28、滤波电感31、直流母线电容29和升压变压器。
双馈异步发电机DFIG 3通过升压变压器连接到电网、转子侧变换器RSC 1的输入端连接到双馈异步发电机DFIG 3的转子绕组、网侧侧变换器GSC 28的输入端通过滤波电感31连接到升压变压器,直流母线电容29连接在转子侧变换器RSC 1和网侧侧变换器GSC 28的输出端之间。
图3是本发明提出的一种双馈异步风力发电系统不对称协调直接功率控制方法的原理图。以一台2MW商用变速恒频DFIG风电系统为例,采用本发明提出电网不对称时双馈异步风力发电系统的协调直接功率控制方法,其中左侧虚线框内是RSC 1的不对称直接功率控制,右侧虚线框内是GSC 28的不对称直接功率控制。通过对RSC 1的不对称直接功率控制消除DFIG 3定子电磁功率和无功功率中的基频ωs和两倍频2ωs振荡,保持双馈异步发电机DFIG 3的电磁转矩恒定;通过对GSC 28的不对称直接功率控制消除直流母线电压和GSC 28输出无功功率中的基频ωs和两倍频2ωs振荡,保持双馈异步风力发电系统向电网输出的总有功、无功功率恒定;
所述的RSC 1不对称直接功率控制,包括以下步骤:
(i)利用三相电压霍尔传感器5-1采集双馈异步感应发电机DFIG 3三相定子电压信号Usabc;利用第一三相电流霍尔传感器4-1采集三相定子电流Isabc,利用第二三相电流霍尔传感器4-2采集和转子电流信号Irabc
(ii)采集得到的三相定子电压信号Usabc经过不对称锁相环8检测得到电网或者定子电压中正序分量的角频率ωs和相位θs;与此同时采用编码器7检测DFIG的转子位置θr,再经过微分器35计算转速ωr;并通过减法器计算得到滑差角频率ωslip=ωsr
(iii)将采集得到的三相定子电压信号Usabc,三相定子电流信号Isabc分别经过第一静止三相到二相坐标变换模块6-1,得到定子静止坐标系中包含正、负序分量的定子电压综合矢量Usαβ,以及包含正、负序和零序分量的定子电流综合矢量Isαβ;以定子电压为例,静止三相到二相坐标变换如下式表达
U sα U sβ = 2 3 1 - 1 2 - 1 2 0 3 2 3 2 U sa U sb U sc ;
(iv)将得到的定子静止坐标系中定子电压综合矢量Usαβ、定子电流综合矢量Isαβ经过定子有功、无功功率计算模块22得到定子输出的瞬时有功功率信号Ps和无功功率信号Qs;其计算方法如下式表达
P s + jQ s = U sαβ × I sαβ ^ = ( U sα I sα + U sβ I sβ ) + j ( U sβ I sβ - U sα I sβ )
(v)将采集得到的三相转子电流信号Irabc通过正转转子速旋转坐标变换模块26,得到转子速旋转坐标系中含正、负序和零序分量的转子电流综合矢量Irαβ;正转转子速旋转坐标变换如下式表达
I rα I rβ = 2 3 cos ( - θ r ) sin ( - θ r ) - sin ( - θ r ) cos ( - θ r ) 1 - 1 2 - 1 2 0 3 2 3 2 I ra I rb I rc
(vi)将得到的定子电流综合矢量Isαβ和转子电流综合矢量Irαβ经过定子电磁功率计算模块(27)得到定子电磁功率信号Pes;定子电磁功率Pes计算方法如下式表达
P es = ω s T e = 3 2 ω s L m ( I sβ I rα - I sα I rβ )
其中,Te为DFIG 3的电磁转矩。
(vii)将定子电磁功率信号Pes和无功功率信号Qs与其参考电磁功率信号和无功功率信号
Figure A200910097499D0012152544QIETU
经过减法器计算得到电磁功率误差信号ΔPes和无功功率误差信号ΔQs
(viii)将得到的电磁功率误差信号ΔPes和无功功率误差信号ΔQs通过第一比例谐振调节器21-1调节;比例谐振控制器21-1的频域表达式为
C PR ( s ) = k p + k r 1 s s 2 + 2 ω c 1 s + ω s 2 + k r 2 s s 2 + 2 ω c 2 s + ( 2 ω s ) 2
其中,kp,kr1,kr2为比例、谐振调节器的系数。调节后的输出信号
Figure A200910097499D00132
经过第一反馈补偿解耦模块15完成转子速旋转坐标系中交-直轴间的交叉解耦和动态反馈补偿,获取转子速旋转坐标系中的转子电压参考值
Figure A200910097499D00133
可用下式表达
U rd + * = V rd + - ω slip ( Q s k σ U s - L r L m U sd ω s )
U rq + * = V rq + - ω slip P es k σ U s
其中
U rd + * = - C PR ( s ) ( P es * - P s ) = - [ k p + k r 1 s s 2 + 2 ω c 1 s + ω s 2 + k r 2 s s 2 + 2 ω c 2 s + ( 2 ω s ) 2 ] ( P es * - P es )
U rq + * = C PR ( s ) ( Q s * - Q s ) = [ k p + k r 1 s s 2 + 2 ω c 1 s + ω s 2 + k r 2 s s 2 + 2 ω c 2 s + ( 2 ω s ) 2 ] ( Q s * - Q s )
(ix)转子电压参考值
Figure A200910097499D00138
经过转差速旋转坐标变换模块9,获得定子坐标系中转子电压参考值
Figure A200910097499D00139
作为空间矢量第一脉宽调制模块2-1调制所需的转子电压参考信号,该信号经过空间矢量脉宽调制后获得控制DFIG运行的RSC 1的开关信号Sa1,Sb1,Sc1。其中转差速旋转坐标变换9可用下式表示
U rα * U rβ * = cos ( θ s - θ r ) sin ( θ s - θ r ) - sin ( θ s - θ r ) cos ( θ s - θ r ) U rd + * U rq + *
所述的GSC(28)不对称直接功率控制,包括以下步骤:
(i)利用单相电压霍尔传感器36采集直流母线电压信号Vdc;利用第三三相电流霍尔传感器4-3采集GSC 28输出的流过滤波电感31的三相电流信号Igabc
(ii)将采集得到的GSC 28三相输出电流信号Igabc经过第二静止三相到二相坐标变换模块6-2,得到静止坐标系中包含正、负序和零序分量的GSC 28电流综合矢量Igαβ
(iii)将得到的静止坐标系中定子电压综合矢量Usαβ、GSC 28电流综合矢量Igαβ经过GSC 28有功、无功功率计算模块34得到GSC 28输出的瞬时有功功率信号Pg和无功功率信号Qg;其计算方法如下式表达
P g + jQ g = 3 2 U sαβ × I gαβ ^ = ( U sα I gα + U sβ I gβ ) + j ( U sβ I gβ - U sα I gβ )
(iv)将直流母线电压参考信号
Figure A200910097499D00142
与采集得到的直流母线电压信号Vdc经过减法器计算得到直流母线电压误差信号,利用第一比例积分调节器19对得到的误差信号作比例-积分调节,调节输出信号经过GSC 28有功功率参考信号补偿模块32得到GSC28有功功率参考信号
Figure A200910097499D00143
;其计算方法如下式表达
P g * = C PI ( s ) ( V dc * - V dc ) + ω r T e - P s
其中比例积分器19的频域表达式CPI(s)为
CPI(s)=kp+ki/s
其中,kp,ki为比例、积分调节器的系数。
(v)将GSC 28输出的有功功率信号Pg和无功功率信号Qg与其参考有功功率信号
Figure A200910097499D00145
和无功功率信号
Figure A200910097499D00146
经过减法器计算得到GSC 28输出有功误差信号ΔPg和无功功率误差信号ΔQg
(vi)将得到的有功误差信号ΔPg和无功功率误差信号ΔQg通过第二比例谐振调节器21-2调节;调节后的输出信号经过第二反馈补偿解耦模块33完成同步速旋转坐标系中交-直轴间的交叉解耦和动态反馈补偿,获取同步速旋转坐标系中的GSC 28输出电压参考值可用下式表达
U cd + * = - C PR ( s ) ( P g * - P g ) - 2 3 ω s L g U sd Q g + U sd
U cq + * = C PR ( s ) ( Q g * - Q g ) - 2 3 ω s L g U sd P g
(vii)GSC 28输出电压参考值
Figure A200910097499D001410
经过同步速旋转坐标变换模块14,获得的静止坐标系中GSC 28输出电压参考值
Figure A200910097499D001411
作为空间矢量脉宽调制模块2调制所需的GSC输出电压参考信号,该信号经过第二脉宽调制模块2-2调制后获得控制GSC 28运行的开关信号Sa2,Sb2,Sc2。其中同步速旋转坐标变换模块14如下式表达
U cα * U cβ * = cos θ s sin θ s - sin θ s cos θ s U cd + * U cq + *
参照图4,上述的RSC 1不对称直接功率控制中所述的采用不对称锁相环8检测电网或者定子电压中正序分量的角频率ωs和相位θs,步骤如下:
(i)利用锁相环的反馈相位信号对三相定子电压信号Usabc进行正转坐标变换23,得到正转坐标系中含有直流量与两倍频2ωs交流量之和的电压综合矢量Usdq;正转坐标变换如下式表达
U sd U sq = 2 3 cos θ s ′ sin θ s ′ - sin θ s ′ cos θ s ′ 1 - 1 2 - 1 2 0 3 2 3 2 U sa U sb U sc
(ii)将得到的正转坐标系中电压综合矢量Usdq的q轴分量Usq经过第二比例积分调节器(37)得到三相定子电压正序分量的频率ωs
(iii)将得到的频率信号ωs经过积分器20积分得到电压正序分量的相位信号θs
(iv)Usq经过第一两倍频2ωs谐振调节器24-1调节后的输出信号与电压正序分量的相位信号θs相加,得到锁相环的反馈相位信号。两倍频2ωs谐振调节器24-1的频域表达式为
C R ( s ) = k r s s 2 + 2 ω c 2 s + ( 2 ω s ) 2
其中,kr为谐振调节器的系数。
参照图5,上述的RSC 1不对称直接功率控制和GSC 28不对称直接功率控制中所述的第一比例谐振调节器21-1和第二比例谐振调节器21-2分别包括一个比例控制器26和两个角频率分别为基频ωs和两倍频2ωs的谐振调节器25、24-2,其中比例控制器26调节功率的平均值,角频率为ωs的谐振调节器25可抑制功率中的基频ωs振荡,角频率为2ωs的第二倍频2ωs谐振调节器24-2抑制功率中的两倍频2ωs振荡。
参照图6(A),若不采用本发明方法,则在电压不对称条件下(0.1-0.3s),DFIG的定子有功、无功功率,电磁转矩,GSC输出有功、无功功率,系统输出总有功、无功功率,以及直流母线电压之中都出现明显的基频ωs和两倍频2ωs振荡,其中基频ωs振荡由电压突变瞬间定子磁链中的直流分量引起,而两倍频2ωs振荡则由电压中的负序分量引起;电压恢复正常(0.3s)之后,由于定子磁链中的直流分量仍然存在,DFIG的定子有功、无功功率,电磁转矩,GSC输出有功、无功功率,系统输出总有功、无功功率,以及直流母线电压之中的基频ωs振荡仍然存在。
参照图6(B),采用本发明方法之后,DFIG的定子无功功率、GSC输出无功功率、系统输出总无功功率、DFIG电磁转矩以及直流母线电压之中的基频ωs和两倍频2ωs振荡被很快抑制;但是定子有功功率中的基频ωs和两倍频2ωs振荡反而增大,此时通过控制GSC输出的有功功率对定子有功功率中的振荡分量进行补偿,从而使得整个风电系统向电网输出的总有功功率保持稳定;电压恢复正常(0.3s)之后,上述所有电量中的振荡分量都在短时间(0.03s)内减小至零,即定子磁链中的直流分量对系统的影响被完全抑制。通过图6(A)和图6(B)的对比,可见采用本发明的电网不对称时双馈异步风力发电系统的协调直接功率控制方法之后,实现了双馈异步风力发电系统在电网不对称故障情况下的整体最优化控制目标。
综上所述,本发明公开的一种双馈异步风力发电系统不对称协调直接功率控制方法无需任何正、负序分解,结构简单,动态响应块,稳态性能好;在电网电压不对称的情况下,可以实现双馈异步风力发电系统向电网输出稳定的有功、功率;同时,在系统内部,可以抑制DFIG和直流母线电压的振荡,从而避免风电机组的机械系统和直流母线电容受到损坏。本方法可增强电网不对称故障情况下对风电机组的控制能力,实现了风电机组在电网故障下的穿越运行。

Claims (3)

1.一种双馈异步风力发电系统不对称协调直接功率控制方法,其特征在于包括双馈异步风力发电系统中转子侧变换器RSC(1)的不对称直接功率控制和网侧变换器GSC(28)的不对称直接功率控制;通过对RSC(1)的不对称直接功率控制消除DFIG(3)定子电磁功率和无功功率中的基频ωs和两倍频2ωs振荡,保持双馈异步发电机DFIG(3)的电磁转矩恒定;通过对GSC(28)的不对称直接功率控制消除直流母线电压和GSC(28)输出无功功率中的基频ωs和两倍频2ωs振荡,保持双馈异步风力发电系统向电网输出的总有功、无功功率恒定;
所述的RSC(1)不对称直接功率控制,包括以下步骤:
(i)利用三相电压霍尔传感器(5-1)采集双馈异步感应发电机DFIG(3)三相定子电压信号Usabc;利用第一三相电流霍尔传感器(4-1)采集三相定子电流Isabc,利用第二三相电流霍尔传感器(4-2)采集三相转子电流信号Irabc
(ii)采集得到的三相定子电压信号Usabc经过不对称锁相环(8)检测得到电网或者定子电压中正序分量的角频率ωs和相位θs;与此同时采用编码器(7)检测DFIG的转子位置θr,再经过微分器(35)计算转速ωr;并通过减法器计算得到滑差角频率ωslip=ωsr
(iii)将采集得到的三相定子电压信号Usabc,三相定子电流信号Isabc分别经过第一静止三相到二相坐标变换模块(6-1),得到定子静止坐标系中包含正、负序分量的定子电压综合矢量Usαβ,以及包含正、负序和零序分量的定子电流综合矢量Isαβ
(iv)将得到的定子静止坐标系中定子电压综合矢量Usαβ、定子电流综合矢量Isαβ经过定子有功、无功功率计算模块(22)得到定子输出的瞬时有功功率信号Ps和无功功率信号Qs
(v)将采集得到的三相转子电流信号Irabc通过正转转子速旋转坐标变换模块(26),得到转子速旋转坐标系中含正、负序和零序分量的转子电流综合矢量Irαβ
(vi)将得到的定子电流综合矢量Isαβ和转子电流综合矢量Irαβ经过定子电磁功率计算模块(27)得到定子电磁功率信号Pes
(vii)将定子电磁功率信号Pes和无功功率信号Qs与其参考电磁功率信号和无功功率信号
Figure A200910097499C0002151910QIETU
经过减法器计算得到电磁功率误差信号ΔPes和无功功率误差信号ΔQs
(viii)将得到的电磁功率误差信号ΔPes和无功功率误差信号ΔQs通过第一比例谐振调节器(21-1)调节;调节后的输出信号
Figure A200910097499C00031
经过第一反馈补偿解耦模块(15)完成转子速旋转坐标系中交-直轴间的交叉解耦和动态反馈补偿,获取转子速旋转坐标系中的转子电压参考值
Figure A200910097499C00032
(ix)转子电压参考值经过转差速旋转坐标变换模块(9),获得定子坐标系中转子电压参考值
Figure A200910097499C0003151940QIETU
作为空间矢量第一脉宽调制模块(2-1)调制所需的转子电压参考信号,该信号经过空间矢量脉宽调制后获得控制DFIG运行的转子侧变换器(1)的开关信号Sa1,Sb1,Sc1
所述的GSC(28)不对称直接功率控制,包括以下步骤:
(i)利用单相电压霍尔传感器(36)采集直流母线电压信号Vdc;利用第三三相电流霍尔传感器(4-3)采集GSC(28)输出的流过滤波电感(31)的三相电流信号Igabc
(ii)将采集得到的GSC(28)三相输出电流信号Igabc经过第二静止三相到二相坐标变换模块(6-2),得到静止坐标系中包含正、负序和零序分量的GSC(28)电流综合矢量Igαβ
(iii)将得到的静止坐标系中定子电压综合矢量Usαβ和GSC(28)电流综合矢量Igαβ经过GSC(28)有功、无功功率计算模块(34)得到GSC(28)输出的瞬时有功功率信号Pg和无功功率信号Qg
(iv)将直流母线电压参考信号
Figure A200910097499C0003152002QIETU
与采集得到的直流母线电压信号Vdc经过减法器计算得到直流母线电压误差信号,利用第一比例积分调节器(19)对得到的误差信号作比例-积分调节,调节输出信号经过GSC(28)有功功率参考信号补偿模块(32)得到GSC(28)有功功率参考信号
Figure A200910097499C0003152029QIETU
(v)将GSC(28)输出的有功功率信号Pg和无功功率信号Qg与其参考有功功率信号
Figure A200910097499C0003152132QIETU
和无功功率信号
Figure A200910097499C0003152140QIETU
经过减法器计算得到GSC(28)输出有功误差信号ΔPg和无功功率误差信号ΔQg
(vi)将得到的有功功率误差信号ΔPg和无功功率误差信号ΔQg通过第二比例谐振调节器(21-2)调节;调节后的输出信号经过第二反馈补偿解耦模块(33)完成同步速旋转坐标系中交-直轴间的交叉解耦和动态反馈补偿,获取同步速旋转坐标系中的GSC(28)输出电压参考值
Figure A200910097499C0003152156QIETU
(vii)GSC(28)输出电压参考值
Figure A200910097499C0003152206QIETU
经过同步速旋转坐标变换模块(14),获得的静止坐标系中GSC(28)输出电压参考值作为第二脉宽调制模块(2-2)调制所需的GSC输出电压参考信号,该信号经过空间矢量脉宽调制后获得控制GSC(28)运行的开关信号Sa2,Sb2,Sc2
2.根据权利要求1所述的一种双馈异步风力发电系统不对称协调直接功率控制方法,其特征在于RSC(1)不对称直接功率控制中所述的采用不对称锁相环(8)检测电网或者定子电压中正序分量的角频率ωs和相位θs,步骤如下:
(i)利用锁相环的反馈相位信号
Figure A200910097499C0004152247QIETU
对三相定子电压信号Usabc进行正转坐标变换(23),得到正转坐标系中含有直流量与两倍频2ωs交流量之和的电压综合矢量Usdq
(ii)将得到的正转坐标系中电压综合矢量Usdq的q轴分量Usq经过第二比例积分调节器(37)得到三相定子电压正序分量的频率ωs
(iii)将得到的频率信号ωs经过积分器(20)积分得到电压正序分量的相位信号θs
(iv)Usq经过第一两倍频2ωs谐振调节器(24-1)调节后的输出信号与电压正序分量的相位信号θs相加,得到锁相环的反馈相位信号
Figure A200910097499C0004152315QIETU
3.根据权利要求1所述的一种双馈异步风力发电系统不对称协调直接功率控制方法,其特征在于RSC(1)不对称直接功率控制和GSC(28)不对称直接功率控制中所述的第一比例谐振调节器(21-1)和第二比例谐振调节器(21-2)分别包括一个比例控制器(26)和两个角频率分别为基频ωs和两倍频2ωs的谐振调节器(25、24-2),其中比例控制器(26)调节功率的平均值,角频率为ωs的谐振调节器(25)可抑制功率中的基频ωs振荡,角频率为2ωs的第二倍频2ωs谐振调节器(24-2)抑制功率中的两倍频2ωs振荡。
CN2009100974991A 2009-04-07 2009-04-07 一种双馈异步风力发电系统不对称协调直接功率控制方法 Expired - Fee Related CN101521481B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100974991A CN101521481B (zh) 2009-04-07 2009-04-07 一种双馈异步风力发电系统不对称协调直接功率控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100974991A CN101521481B (zh) 2009-04-07 2009-04-07 一种双馈异步风力发电系统不对称协调直接功率控制方法

Publications (2)

Publication Number Publication Date
CN101521481A true CN101521481A (zh) 2009-09-02
CN101521481B CN101521481B (zh) 2011-01-12

Family

ID=41081893

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100974991A Expired - Fee Related CN101521481B (zh) 2009-04-07 2009-04-07 一种双馈异步风力发电系统不对称协调直接功率控制方法

Country Status (1)

Country Link
CN (1) CN101521481B (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102158174A (zh) * 2010-01-22 2011-08-17 Abb研究有限公司 旋转电机的控制器
CN103281025A (zh) * 2013-05-17 2013-09-04 浙江大学 一种基于谐振滑模的dfig系统控制方法
CN103401498A (zh) * 2010-02-26 2013-11-20 株式会社电装 用于旋转电机的控制设备
CN104362668A (zh) * 2014-10-16 2015-02-18 中国人民解放军装甲兵工程学院 一种电压不平衡/谐波畸变时双馈风力发电机的控制方法
US9093924B2 (en) 2010-08-18 2015-07-28 Vestas Wind Systems A/S Method of controlling a grid side converter of a wind turbine and system suitable therefore
CN107124126A (zh) * 2017-03-31 2017-09-01 中国电力科学研究院 一种双馈感应发电机的无锁相环电流控制方法及装置
CN109510548A (zh) * 2018-09-29 2019-03-22 北方工业大学 一种双馈电机柔性功率控制方法及装置
CN109752654A (zh) * 2017-11-07 2019-05-14 西门子歌美飒可再生能源公司 诊断定子的绕组集合
CN110412464A (zh) * 2019-08-05 2019-11-05 黑龙江科技大学 新型研究水轮发电机负序和零序分量电磁功率的方法
CN111404191A (zh) * 2020-04-19 2020-07-10 中国电建集团华东勘测设计研究院有限公司 基于谐振控制器的mmc换流站低电压穿越方法及控制系统
CN112087166A (zh) * 2020-09-17 2020-12-15 华中科技大学 交直流混合型双馈异步全电船舶电力推进系统及控制方法
CN113162486A (zh) * 2021-04-13 2021-07-23 中国人民解放军国防科技大学 一种双边空芯直线同步电机牵引-导向解耦控制方法
CN113193808A (zh) * 2021-04-25 2021-07-30 华中科技大学 一种容错型双馈异步全电船舶电力传动系统的控制方法
CN113826318A (zh) * 2019-05-22 2021-12-21 西门子歌美飒可再生能源创新与技术有限公司 用于风力涡轮机中的双馈感应发电机的转子侧变流器的控制系统和方法

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104811115B (zh) * 2015-04-15 2017-03-22 哈尔滨工业大学 基于准比例谐振控制的永磁同步电机电感参数辨识系统及方法
CN110034582A (zh) * 2019-04-24 2019-07-19 国电联合动力技术有限公司 风力发电变流器功率不平衡控制系统及方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU3148893A (en) * 1991-11-27 1993-06-28 U.S. Windpower, Inc. Variable speed wind turbine with reduced power fluctuation and a static var mode of operation
EP1499009B1 (en) * 2003-07-15 2007-10-31 Gamesa Innovation & Technology, S.L. Unipersonal Control and protection of a doubly-fed induction generator system
GB2420456A (en) * 2004-11-23 2006-05-24 Areva T & D Uk Ltd Generator control having grid imbalance detector
US7425771B2 (en) * 2006-03-17 2008-09-16 Ingeteam S.A. Variable speed wind turbine having an exciter machine and a power converter not connected to the grid
GB2451463B (en) * 2007-07-28 2012-07-25 Converteam Technology Ltd Control methods for VSC active rectifier/inverters under unbalanced operating conditions

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102158174A (zh) * 2010-01-22 2011-08-17 Abb研究有限公司 旋转电机的控制器
CN102158174B (zh) * 2010-01-22 2015-11-25 Abb研究有限公司 旋转电机的控制器
CN103401498A (zh) * 2010-02-26 2013-11-20 株式会社电装 用于旋转电机的控制设备
US9093924B2 (en) 2010-08-18 2015-07-28 Vestas Wind Systems A/S Method of controlling a grid side converter of a wind turbine and system suitable therefore
CN103281025A (zh) * 2013-05-17 2013-09-04 浙江大学 一种基于谐振滑模的dfig系统控制方法
CN103281025B (zh) * 2013-05-17 2015-04-22 浙江大学 一种基于谐振滑模的dfig系统控制方法
CN104362668A (zh) * 2014-10-16 2015-02-18 中国人民解放军装甲兵工程学院 一种电压不平衡/谐波畸变时双馈风力发电机的控制方法
CN104362668B (zh) * 2014-10-16 2016-08-24 中国人民解放军装甲兵工程学院 一种电压不平衡/谐波畸变时双馈风力发电机的控制方法
CN107124126A (zh) * 2017-03-31 2017-09-01 中国电力科学研究院 一种双馈感应发电机的无锁相环电流控制方法及装置
CN107124126B (zh) * 2017-03-31 2021-09-03 中国电力科学研究院 一种双馈感应发电机的无锁相环电流控制方法及装置
CN109752654A (zh) * 2017-11-07 2019-05-14 西门子歌美飒可再生能源公司 诊断定子的绕组集合
CN109510548B (zh) * 2018-09-29 2020-05-19 北方工业大学 一种双馈电机柔性功率控制方法及装置
CN109510548A (zh) * 2018-09-29 2019-03-22 北方工业大学 一种双馈电机柔性功率控制方法及装置
CN113826318A (zh) * 2019-05-22 2021-12-21 西门子歌美飒可再生能源创新与技术有限公司 用于风力涡轮机中的双馈感应发电机的转子侧变流器的控制系统和方法
CN113826318B (zh) * 2019-05-22 2024-07-02 西门子歌美飒可再生能源创新与技术有限公司 用于风力涡轮机中的转子侧变流器的控制系统和方法
CN110412464A (zh) * 2019-08-05 2019-11-05 黑龙江科技大学 新型研究水轮发电机负序和零序分量电磁功率的方法
CN111404191A (zh) * 2020-04-19 2020-07-10 中国电建集团华东勘测设计研究院有限公司 基于谐振控制器的mmc换流站低电压穿越方法及控制系统
CN111404191B (zh) * 2020-04-19 2021-10-22 中国电建集团华东勘测设计研究院有限公司 基于谐振控制器的mmc换流站低电压穿越方法及系统
CN112087166A (zh) * 2020-09-17 2020-12-15 华中科技大学 交直流混合型双馈异步全电船舶电力推进系统及控制方法
CN112087166B (zh) * 2020-09-17 2022-05-27 华中科技大学 交直流混合型双馈异步全电船舶电力推进系统及控制方法
CN113162486A (zh) * 2021-04-13 2021-07-23 中国人民解放军国防科技大学 一种双边空芯直线同步电机牵引-导向解耦控制方法
CN113193808A (zh) * 2021-04-25 2021-07-30 华中科技大学 一种容错型双馈异步全电船舶电力传动系统的控制方法

Also Published As

Publication number Publication date
CN101521481B (zh) 2011-01-12

Similar Documents

Publication Publication Date Title
CN101521481B (zh) 一种双馈异步风力发电系统不对称协调直接功率控制方法
CN101505131B (zh) 一种双馈异步风力发电机的不对称直接功率控制方法
CN100527595C (zh) 交流励磁双馈异步风力发电机转子电流无延时控制方法
CN101534065B (zh) 一种并网三相电压源变换器的不对称直接功率控制方法
CN108683198B (zh) 双馈风力发电机组的电压控制型虚拟同步方法
CN107425539B (zh) 电网不对称故障下双馈风电机组的增强低电压穿越控制方法
CN101141111B (zh) 一种双馈异步风力发电机转子电流无延时控制方法
CN101453187B (zh) 一种电网不对称故障时风电机组控制基准信号的检测方法
CN101478283B (zh) 不平衡电网电压下双馈异步风力发电机转子侧逆变器控制方法
CN102082543B (zh) 不平衡电网电压下双馈风力发电机控制方法
CN101741096A (zh) 并网型变速恒频双馈感应风力发电机转子电流无延时控制方法
CN103117699B (zh) 一种基于双矢量谐振调节双馈异步风力发电机的控制方法
CN103117700B (zh) 一种不平衡电网下基于谐振反馈的dfig控制方法
CN101977006B (zh) 电网故障情况下双馈风力发电机的控制方法
CN105958534B (zh) 一种双馈风电系统不对称高电压故障穿越控制方法
CN107658911B (zh) 电网不对称故障下永磁直驱风电机组增强低电压穿越控制方法
CN101141110A (zh) 变速恒频双馈异步风力发电机转子电流无延时控制方法
CN104362668A (zh) 一种电压不平衡/谐波畸变时双馈风力发电机的控制方法
CN101557190A (zh) 不平衡电网电压下双馈异步风力发电机网侧整流器控制方法
CN103219736B (zh) 一种飞轮储能单元抑制永磁直驱风力发电系统直流侧2倍频波动的控制方法
CN108667080B (zh) 一种不平衡电网电压下的虚拟同步机有功平衡控制方法
CN105024607A (zh) 一种不平衡电网电压下基于矩阵变换器励磁的dfig控制方法
CN103346718A (zh) 电压不平衡下双转子永磁风力发电系统的控制方法
CN105024608A (zh) 一种不平衡电网下基于矩阵变换器的pmsg比例积分谐振控制方法
CN112165106A (zh) 基于双二阶广义积分器的dfig次同步振荡抑制方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110112

Termination date: 20110407