CN112019118B - 一种直驱风电变流器无定子电压测量矢量控制方法 - Google Patents
一种直驱风电变流器无定子电压测量矢量控制方法 Download PDFInfo
- Publication number
- CN112019118B CN112019118B CN202010864025.1A CN202010864025A CN112019118B CN 112019118 B CN112019118 B CN 112019118B CN 202010864025 A CN202010864025 A CN 202010864025A CN 112019118 B CN112019118 B CN 112019118B
- Authority
- CN
- China
- Prior art keywords
- stator
- generator
- current
- flux linkage
- axis
- 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
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/22—Current control, e.g. using a current control loop
-
- 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/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
一种直驱风电变流器无定子电压测量矢量控制方法,包含定子电流计算环节、初始磁链计算环节、磁链锁相环节、参考电流计算环节、电流闭环控制环节五部分。定子电流计算环节采集三相定子电流Isa、Isb、Isc,转化为αβ坐标系上的定子电流Isα、Isβ和DQ坐标上的定子电流Isd、Isq;初始磁链计算环节通过输出7个ΔT0时间的零矢量,分8个阶段确定允许起动时刻的发电机角频率ω0和磁链角度初值θψ0,输出机侧变流器允许起动标识ST_Per;磁链锁相环节在ST_Per=1后在ω0、θψ0的基础上进行锁相,输出发电机角频率ω和磁链角度θψd;参考电流计算环节输出DQ坐标系上的定子参考电流IsdRef、IsqRef;电流闭环控制环节采用PI调节器对DQ坐标系上的定子电流进行闭环控制,输出αβ坐标系上的定子控制电压Usca、Uscβ和DQ坐标系上的定子控制电压Uscd、Uscq。
Description
技术领域
本发明涉及一种直驱风电变流器无定子电压测量矢量控制方法。
背景技术
近年来随着绿色能源的大力推广与发展,我国风电事业得到了飞速的发展,累计装机总量和年新增装机总量均位居世界第一。随着陆上低风速风电场的开发和海上风电的快速发展,陆上及海上风电采用大兆瓦级风电机组已达成共识,功率等级基本从4MW起步,一般采用基于永磁同步发电机的直驱技术路线,变流器采用三相电压型交-直-交双向变流器。受制于产业链影响,目前大功率海上直驱风电机组电压等级大多还是沿用陆上风电的690V低压技术路线,存在电缆过多引起的系列问题,越来越多的整机厂尝试采用1.14kV或3.3kV中压技术路线,对应的变流器功率单元为二极管箝位三电平拓扑。
直驱风电变流器分为机侧变流器和网侧变流器两部分,其中机侧变流器一般采用磁链定向矢量控制方式,需要定子电压、直流母线电压和输出电流信息。直流母线电压测量用于保证直流电压的稳定,交流电流测量提供电流反馈信号,实现过流保护,两者一般不宜省去。由于机侧变流器直接或者通过小感值的du/dt电感间接连接到发电机定子,输出为基波频率不固定的高频脉宽波,定子电压无法通过常规PT进行测量,需要采用高成本的传感器隔离方案或低成本的电阻分压非隔离方案。690V低压技术路线中一般采用电阻分压非隔离方案,这种方案在电压等级越来越高时安全隐患逐渐凸显。为了降低变流器成本和规避中压技术路线相关的安全隐患,最好方法是采用无定子电压测量的变流器控制技术。
专利201510755713.3“一种无电压传感器的三相变流器模型预测控制方法”,公开了一种基于二次规划求解算法寻优的无电压传感器的变流器模型预测控制方法,可以减小计算量、减少参数整定的难度,但控制对象仅是频率固定的电网,同时未考虑控制方法的磁链初值问题;论文《无电网电压传感器三相PWM整流器虚拟电网磁链定向矢量控制研究》(赵仁德,贺益康等)提出一种基于电网虚拟磁链的无电网电压传感器三相PWM整流器矢量控制技术,在控制方法运行之前通过输出连续几个PWM周期的零矢量获取控制方法中的虚拟磁链初值,仅适用于频率固定的电网,且应用在大兆瓦海上直驱风电变流器时连续输出几个PWM周期的零矢量会引起严重的过电流问题。
发明内容
本发明克服现有技术的缺点,提出一种针对直驱风电变流器的无定子电压测量矢量控制方法。
直驱风电变流器分为网侧变流器和机侧变流器两部分,其中机侧变流器包含二极管钳位三电平拓扑组成的主电路和控制单元,所述的矢量控制方法在机侧变流器控制单元中实现。控制单元与上位机通讯得到有功功率给定值PGiven,采集主电路中的三相定子电流Isa、Isb、 Isc和直流母线电压Udc,输出12路PWM调制信号SPWM到主电路中的电力电子器件。三相定子电流Isa、Isb、Isc以从机侧变流器到发电机定子方向为正方向,忽略发电机定子电阻Rs,认为定子空载电动势角度超前发电机磁链角度
所述控制方法包含定子电流计算环节、初始磁链计算环节、磁链锁相环节、参考电流计算环节、电流闭环控制环节五部分。
所述的定子电流计算环节采集三相定子电流Isa、Isb、Isc,将三相定子电流Isa、Isb、Isc分别转化为α轴定子电流Isα、β轴定子电流Isβ、D轴定子电流Isd和Q轴定子电流 Isq。所述的初始磁链计算环节通过输出7个ΔT0时间的零矢量,分8个阶段确定机侧变流器允许起动时刻,确定允许起动时刻的发电机角频率ω0和发电机磁链角度初值θψ0,在允许起动时刻输出机侧变流器允许起动标识ST_Per;磁链锁相环节在收到机侧变流器允许起动标识ST_Per后,在允许起动时刻的发电机角频率ω0和发电机磁链角度初值θψ0的基础上进行锁相,输出发电机角频率ω和发电机磁链角度θψ。所述的参考电流计算环节输出D 轴定子参考电流IsdRef和Q轴定子参考电流IsqRef。所述的电流闭环控制环节采用PI调节器对DQ坐标系上的定子电流进行闭环控制,输出α轴定子控制电压Usca、β轴定子控制电压 Uscβ、D轴定子控制电压Uscd和Q轴定子控制电压Uscq。
定子电流计算环节采集发电机三相定子电流Isa、Isb、Isc,输入磁链锁相环节输出的发电机磁链角度θψ,输出α轴定子电流Isα、β轴定子电流Isβ、D轴定子电流Isd、Q 轴定子电流Isq。初始磁链计算环节输入机侧变流器起动命令CMD_ST,输入定子电流计算环节输出的α轴定子电流Isα、β轴定子电流Isβ,输出机侧变流器允许起动标识ST_Per、允许起动时刻的发电机角频率ω0、发电机磁链角度初值θψ0。磁链锁相环节输入电流闭环控制环节输出的D轴定子控制电压Uscd、Q轴定子控制电压Uscq,输入定子电流计算环节输出的D轴定子电流Isd、Q轴定子电流Isq,输入初始磁链计算环节输出的机侧变流器允许起动标识ST_Per、允许起动时刻的发电机角频率ω0、发电机磁链角度初值θψ0,输出发电机磁链角度θψ、发电机角频率ω。参考电流计算环节输入有功功率给定值PGiven,输入磁链锁相环节输出发电机交频率ω,输入电流闭环控制环节输出的D轴定子控制电压Uscd、 Q轴定子控制电压Uscq,输出D轴定子参考电流IdRef、Q轴定子参考电流IqRef。电流闭环控制环节输入参考电流计算环节输出的D轴定子参考电流IdRef、Q轴定子参考电流IqRef,输入定子电流计算环节输出的D轴定子电流Isd、Q轴定子电流Isq,输入磁链锁相环节输出的发电机角频率ω、发电机磁链角度θψ,输出D轴定子控制电压Uscd、Q轴定子控制电压 Uscq、α轴定子控制电压Uscα和β轴定子控制电压Uscβ。
电流闭环控制环节输出的D轴定子控制电压Uscd和Q轴定子控制电压Uscq用于下一个控制周期的磁链锁相和参考电流计算,电流闭环控制环节输出的α轴定子控制电压Usca和β轴定子控制电压Uscβ经PWM调制策略后输出12路PWM调制信号SPWM,用于控制主电路三相逆变桥中的电力电子器件开通与关断,实现本发明所述的矢量控制。
以下对各个环节的功能和实现方法进行详细说明。
(1)定子电流计算环节
在定子电流计算环节中,采集的三相定子电流Isa、Isb、Isc经过3s/2s坐标变换得到αβ坐标系上的α轴定子电流Isα和β轴定子电流Isβ;α轴定子电流Isα和β轴定子电流Isβ经过基于发电机磁链角度θψ的2s/2r坐标变换得到DQ坐标系上的D轴定子电流Isd和Q轴定子电流Isq;
(2)初始磁链计算环节
在初始磁链计算环节中,设定参数包含发电机磁链ψSet、零矢量作用时间ΔT0,变流器输出最高频率fMax。采用短时间输出零矢量测量发电机定子电流的方式获取定子空载电动势角度,当采用零矢量将发电机定子短时短路时,发电机在αβ坐标系上的数学方程如公式 (1)所示。
公式(1)中,Isα为α轴定子电流,Isβ为β轴定子电流,Esα为α轴定子空载电动势、 Esβ为β轴定子空载电动势,Ld为发电机直轴电感,ΔT0为零矢量作用时间。
定子空载电动势角度计算公式为θe=a tan 2(Esβ,Esα),将公式(1)代入,可求得机侧变流器输出零矢量时刻的定子空载电动势角度θe,如公式(2)所示。
θe=a tan 2(-Isβ,-Isα) (2)
当机侧变流器起动命令CMD_ST从0变为1时,按照8个阶段完成磁链初值计算流程,不同阶段用Flux_Step表示。
阶段1(Flux_Step=0):检测到机侧变流器起动命令CMD_ST从0变1时,第1次输出ΔT0时间的零矢量,调用定子电流计算环节输出α轴定子电流Isα和β轴定子电流Isβ,根据公式(2)计算第1次输出零矢量时刻的定子空载电动势角度θe1;同时根据公式(3) 计算第2、3、4次输出零矢量的时间间隔TSet1,将定时器Tcnt清0,将Flux_Step设置为1;
公式(3)中,fMax为变流器输出最高频率。
阶段2(Flux_Step=1):定时器计时到TSet1后,第2次输出ΔT0时间的零矢量,调用定子电流计算环节输出α轴定子电流Isα和β轴定子电流Isβ,根据公式(2)计算第2次输出零矢量时刻的定子空载电动势角度θe2;同时将定时器Tcnt清0,将Flux_Step设置为2;
阶段3(Flux_Step=2):定时器计时到TSet1后,第3次输出ΔT0时间的零矢量,调用定子电流计算环节输出α轴定子电流Isα和β轴定子电流Isβ,根据公式(2)计算第3次输出零矢量时刻的定子空载电动势角度θe3;同时将定时器Tcnt清0,将Flux_Step设置为3;
阶段4(Flux_Step=3):定时器计时到TSet1后,第4次输出ΔT0时间的零矢量,调用定子电流计算环节输出α轴定子电流Isα和β轴定子电流Isβ,根据公式(2)计算第4 次输出零矢量时刻的定子空载电动势角度θe4。分别计算每个TSet1间隔的定子空载电动势角度差值,如公式(4)所示。
对3次空载电动势角度差值求取平均值,根据平均值计算1~4阶段的发电机频率f1,如公式(5)所示。
计算第5次输出零矢量的时间间隔TSet2,同时将定时器Tcnt清0,Flux_Step设置为4。 TSet2计算方法如公式(6)所示。
公式(6)中,f1为初始磁链计算1~4阶段的发电机频率,θe4为第4次输出零矢量时对应的定子空载电动势角度。
阶段5(Flux_Step=4):定时器计时到TSet2后,第5次输出ΔT0时间的零矢量,调用定子电流计算环节输出α轴定子电流Isα和β轴定子电流Isβ,根据公式(2)计算第5次输出零矢量时刻的定子空载电动势角度θe5;在发电机A相空载电动势最大值处确认空载电动势角度,定子空载电动势角度θe5角度应该在0附近,误差小于认为角度估算正确;同时采用公式(7)计算第6次输出零矢量的时间间隔TSet3,将定时器Tcnt清0,将Flux_Step 设置为5。
阶段6(Flux_Step=5):定时器计时到TSet3后,第6次输出ΔT0时间的零矢量,调用定子电流计算环节输出α轴定子电流Isα和β轴定子电流Isβ,根据公式(2)计算第6次输出零矢量时刻的定子空载电动势角度θe6,间隔确认空载电动势角度,f1>0时定子空载电动势角度θe6角度应该在附近,f1<0时定子空载电动势角度θe6角度应该在附近,误差小于认为角度估算正确;同时定时器Tcnt清0,将Flux_Step设置为6;
阶段7(Flux_Step=6):定时器计时到TSet3后,第7次输出ΔT0时间的零矢量,调用定子电流计算环节输出α轴定子电流Isα和β轴定子电流Isβ,根据公式(2)计算第7次输出零矢量时刻的定子空载电动势角度θe7,间隔再次确认空载电动势角度,f1>0时定子空载电动势角度θe7角度应该在附近,f1<0时定子空载电动势角度θe7角度应该在附近,误差小于认为角度估算正确;计算每个TSet3间隔的定子空载电动势角度差值,如公式(8)所示。
对2次定子空载电动势角度差值求取平均值,根据平均值计算5~7阶段的发电机频率f2,如公式(9)所示。
计算到机侧变流器起动允许标识有效的时间间隔TSet4,同时定时器Tcnt清0,将Flux_Step设置为7。TSet4计算方法如公式(10)所示。
公式(10)中,f2为初始磁链计算5~7阶段的发电机频率,θe7为第7次输出零矢量时对应的定子空载电动势角度。
阶段8(Flux_Step=7):定时器计时到TSet4后,机侧变流器允许起动标识ST_Per置1,Flux_Step设置为8。采用公式(11)计算允许起动时刻的发电机角频率ω0,采用公式(12)计算发电机磁链角度初值θψ0。
ω0=2πf2 (11)
在公式(11)(12)中,f2为初始磁链计算第5~7阶段的发电机频率。在定子空载电动势角度为0时的机侧变流器允许起动时刻,如果发电机正转(频率为正),发电机磁链角度初值θψ0为如果发电机反转(频率为负),发电机磁链角度初值θψ0为
(3)磁链锁相环节
磁链锁相环节在机侧变流器允许起动标识ST_Per置1后开始工作。在磁链锁相环节中,设定参数包括发电机交轴电感Lq,发电机磁链ψSet,DQ轴定子电压耦合系数Kcoff。
输入的Q轴定子参考电流Isq乘以发电机角频率ω,再乘以发电机交轴电感Lq得到D轴电枢反应电压Usdar,输入的D轴定子控制电压Uscd加上D轴电枢反应电压Usdar得到D 轴定子空载电动势Esd,对D轴定子空载电动势Esd求取积分得到发电机D轴磁链ψsd,发电机D轴磁链ψsd的积分初值为发电机磁链ψSet;
输入的Q轴定子控制电压Uscq减去D轴定子控制电压Uscd、DQ轴定子电压耦合系数Kcoff、发电机角频率ω三者的乘积得到Q轴定子空载电动势Esq,Q轴定子空载电动势Esq除以发电机D轴磁链ψsd得到发电机角频率ω,由于发电机角频率ω求解涉及到除法运算,在机侧变流器允许起动时刻无法得到Q轴定子空载电动势Esq信息,为此将除法运算屏蔽,将输出的发电机角频率ω设置为初始磁链计算环节输出的允许起动时刻的发电机角频率ω0;
对发电机角频率ω求取积分得到发电机磁链角度θψ,发电机磁链角度θψ的积分初值为初始磁链计算环节输出的发电机磁链角度初值θψ0;
(4)参考电流计算环节
参考电流计算环节在机侧变流器允许起动标识ST_Per置1后开始工作,参考电流计算环节设定参数包括发电机磁链ψSet。
定子参考电压UsRef和定子电压有效值UsRms相减得到定子电压误差UsErr,定子电压误差UsErr经PI调节器输出D轴定子参考电流IsdRef;根据输入的有功功率给定值PGiven、D 轴定子控制电压Uscd、Q轴定子控制电压Uscq和D轴定子参考电流IsdRef计算Q轴定子参考电流IqRef,计算公式为
(5)电流闭环控制环节
电流闭环控制环节在机侧变流器允许起动标识ST_Per置1后开始工作。电流闭环控制环节设定参数包含发电机直轴电感Ld、发电机交轴电感Lq、发电机磁链ψSet。
D轴定子参考电流IsdRef减去D轴定子电流Isd得到D轴定子电流误差IsdErr,D轴定子电流误差IsdErr经PI调节器输出D轴调节电压UsdPI;Q轴定子参考电流IsqRef减去Q 轴定子电流Isq得到Q轴定子电流误差IsqErr,Q轴定子电流误差IsqErr经PI调节器输出Q 轴调节电压UsqPI;
求取D轴定子控制电压Uscd和Q轴定子控制电压Uscq,计算方法如公式(13)所示。
D轴定子控制电压Uscd和Q轴定子控制电压Uscq经过基于发电机磁链角度θψ的2r/2s坐标变换,得到α轴定子控制电压Uscα和β轴定子控制电压Uscβ。
本发明的优点在于,无定子电压测量实现发电机矢量控制,省去定子电压测量既可节约成本又可以有效避免高电压测量带来的安全隐患。在机侧变流器允许起动之前,通过输出7 个很短时间ΔT0的零矢量来判定发电机角频率和发电机磁链角度,可以在无定子电压测量的情况下实现基本无电流冲击起动,有效避免了机侧变流器起动暂态过程中的过电流。在初始磁链计算环节输出的零矢量作用时间ΔT0在几十us量级,保证发电机磁链角度初值获取准确的前提下可以将磁链测试电流控制在额定电流的10%以内,有效避免了发电机磁链角度初值获取输出零矢量过程中的过电流和对发电机运行的影响。
附图说明
图1直驱风电机组及其变流器电气原理图;
图2机侧变流器主电路和控制单元原理图;
图3无定子电压测量矢量控制方法及各环节信号传递框图;
图4初始磁链计算环节流程图;
图5磁链锁相环节原理图;
图6本发明机侧变流器典型电压电流波形。
具体实施方式
以下结合附图和具体实施方式进一步说明本发明。
如图1所示,直驱风电机组包含永磁同步发电机(PMSG)、直驱风电变流器和电网三部分。
其中直驱风电变流器分为网侧变流器110和机侧变流器130两部分,网侧变流器110的交流侧连接电网100,机侧变流器130的交流侧连接到发电机140的定子,网侧变流器110 和机侧变流器130通过直流母线120连接在一起,实现发电机140和电网100的柔性连接。
风电变流器发出的电能并入电网100,需要满足电网并网导则和电网适应性要求,且电网频率固定,电压谐波含量较低,采用常规的PT即可对电网电压进行低成本测量。由于PT 的作用将电网的高电压和控制用的低电压进行电气隔离有效降低了安全隐患,网侧变流器 110的交流电压测量不能省略也没有省略的必要,本发明主要研究机侧变流器130相关的无定子电压测量的矢量控制方法。
图1所示的机侧变流器130电气原理如图2所示。机侧变流器130包含主电路200和控制单元210两部分。主电路200包含直流支撑电容C1、C2,如图中201,包含S11、S12、 S13、S14、S21、S22、S23、S24、S31、S32、S33、S34共12支IGBT和D11、D12、 D21、D22、D31、D32共6支二极管组成的二极管箝位三电平三相逆变桥;三相逆变桥输出连接到永磁同步发电机定子,如图2中202;主电路PN两端连接到直流母线,如图中206;本发明所述的矢量控制方法在控制单元210中实现,控制单元采集主电路中的直流母线电压 Udc和三相定子电流Isa、Isb、Isc,如3图中203、204;接收上位机发送的有功功率给定值PGiven和机侧变流器起动命令CMD_ST,如图2中220、221;输出12路PWM调制信号 SPWM到主电路中三相逆变桥电力电子器件,如图2中205。
如图3所示,本发明所述的无定子电压测量矢量控制方法包含定子电流计算环节300、初始磁链计算环节310、磁链锁相环节320、参考电流计算环节330、电流闭环控制环节340 五部分。
定子电流计算环节300采集发电机三相定子电流Isa、Isb、Isc,如图3中301、302、303,输入磁链锁相环节输出的发电机磁链角度θψ,如图3中304,输出α轴定子电流Isα、β轴定子电流Isβ、D轴定子电流Isd、Q轴定子电流Isq,如图3中305、306、307、308;
初始磁链计算环节310输入机侧变流器起动命令CMD_ST,如图3中311,输入定子电流计算环节输出的α轴定子电流Isα、β轴定子电流Isβ,如图3中312、313,输出机侧变流器允许起动标识ST_Per、允许起动时刻的发电机角频率ω0、发电机磁链角度初值θψ0,如图3中314、315、316;
磁链锁相环节320输入电流闭环控制环节输出的D轴定子控制电压Uscd、Q轴定子控制电压Uscq,如图3中325、326,输入定子电流计算环节输出的Q轴定子电流Isq,如图3中 324,输入初始磁链计算环节输出的机侧变流器允许起动标识ST_Per、允许起动时刻的发电机角频率ω0、发电机磁链角度初值θψ0,如图3中323、321、322,输出发电机磁链角度θψ、发电机角频率ω,如图3中328、327;
参考电流计算环节330输入机侧变流器允许起动标识ST_Per,如图3中332,输入有功功率给定值PGiven,如图3中335,输入磁链锁相环节输出的发电机角频率ω,如图3 中331,输入电流闭环控制环节输出的D轴定子控制电压Uscd、Q轴定子控制电压Uscq,如图3中333、334,输出D轴定子参考电流IdRef、Q轴定子参考电流IqRef,如图3中336、 337;
电流闭环控制环节340输入参考电流计算环节输出的D轴定子参考电流IdRef、Q轴定子参考电流IqRef,如图3中336、337,输入定子电流计算环节输出的D轴定子电流Isd、Q轴定子电流Isq,如图3中341、342,输入磁链锁相环节输出的发电机角频率ω、发电机磁链角度θψ,如图3中327、328;输出D轴定子控制电压Uscd、Q轴定子控制电压Uscq,如图3中345、346,输出α轴定子控制电压Uscα和β轴定子控制电压Uscβ,如图3中343、 344。
在图3所示的定子电流计算环节300中,采集三相定子电流Isa、Isb、Isc,经过3s/2s坐标变换得到αβ坐标系上的α轴定子电流Isα和β轴定子电流Isβ;α轴定子电流Isα和β轴定子电流Isβ经过基于发电机磁链角度θψ的2s/2r坐标变换得到DQ坐标系上的D轴定子电流Isd和Q轴定子电流Isq。
本发明所述的无定子电压测量矢量控制方法包含以下步骤:
步骤1:控制单元上电初始化,从存储器中读取运行需要的参数,包含:发电机磁链ψSet、发电机直轴电感Ld、发电机交轴电感Lq、变流器输出最高频率fMax、DQ轴定子电压耦合系数Kcoff、零矢量作用时间ΔT0。
步骤2:图3所示的初始磁链计算环节310检测到机侧变流器起动命令CMD_ST=1时进行初始磁链计算,如图4中400、401,初始磁链计算分为8个阶段,用Flux_Step区分,如图4中402。
阶段1(Flux_Step=0):检测到机侧变流器起动命令CMD_ST从0变1时,第1次输出ΔT0时间的零矢量,调用定子电流计算环节300输出α轴定子电流Isα和β轴定子电流 Isβ,根据公式(2)计算第1次输出零矢量时刻的定子空载电动势角度θe1;同时根据公式 (3)计算第2、3、4次输出零矢量的时间间隔Tset1,将定时器Tcnt清0,如图4中410、411;将Flux_Step设置为1,切换到阶段2,如图4中412。
阶段2(Flux_Step=1):定时器计时到TSet1后,第2次输出ΔT0时间的零矢量,调用定子电流计算环节300输出α轴定子电流Isα和β轴定子电流Isβ,根据公式(2)计算第2次输出零矢量时刻的定子空载电动势角度θe2;同时将定时器Tcnt清0,如图4中415、 416;将Flux_Step设置为2,切换到阶段3,如图4中417;
阶段3(Flux_Step=2):定时器计时到TSet1后,第3次输出ΔT0时间的零矢量,调用定子电流计算环节300输出α轴定子电流Isα和β轴定子电流Isβ,根据公式(2)计算第 3次输出零矢量时刻的定子空载电动势角度θe3;同时将定时器Tcnt清0,如图4中420、421;将Flux_Step设置为3,切换到阶段4,如图4中422;
阶段4(Flux_Step=3):定时器计时到TSet1后,第4次输出ΔT0时间的零矢量,调用定子电流计算环节300输出α轴定子电流Isα和β轴定子电流Isβ,采用公式(2)计算第4次输出零矢量时刻的定子空载电动势角度θe4,采用公式(4)分别计算每个TSet1间隔的定子空载电动势角度差值Δθe1、Δθe2、Δθe3,采用公式(5)对3次空载电动势角度差值求取平均值,根据平均值计算1~4阶段的发电机频率f1。采用公式(6)计算第5次输出零矢量的时间间隔TSet2,同时将定时器Tcnt清0,如图4中425、426所示,将Flux_Step设置为4,切换到阶段5,如图4中427;
阶段5(Flux_Step=4):定时器计时到TSet2后,第5次输出ΔT0时间的零矢量,调用定子电流计算环节300输出α轴定子电流Isα和β轴定子电流Isβ,根据公式(2)计算第 5次输出零矢量时刻的定子空载电动势角度θe5;在发电机A相空载电动势最大值处确认空载电动势角度,定子空载电动势角度θe5角度应该在0附近,误差小于认为角度估算正确;同时采用公式(7)计算第6次输出零矢量的时间间隔TSet3,将定时器Tcnt清0,如图4中430、 431;将Flux_Step设置为5,切换到阶段6,如图4中432;
阶段6(Flux_Step=5):定时器计时到TSet3后,第6次输出ΔT0时间的零矢量,调用定子电流计算环节300输出α轴定子电流Isα和β轴定子电流Isβ,根据公式(2)计算第 6次输出零矢量时刻的定子空载电动势角度θe6,间隔确认空载电动势角度,f1>0时定子空载电动势角度θe6角度应该在附近,f1<0时定子空载电动势角度θe6角度应该在附近,误差小于认为角度估算正确;同时定时器Tcnt清0,如图4中435、436;将Flux_Step 设置为6,切换到阶段7,如图4中437;
阶段7(Flux_Step=6):定时器计时到TSet3后,第7次输出ΔT0时间的零矢量,,调用定子电流计算环节300输出α轴定子电流Isα和β轴定子电流Isβ,采用公式(2)计算第 7次输出零矢量时刻的定子空载电动势角度θe7,间隔再次确认空载电动势角度,f1>0时定子空载电动势角度θe7角度应该在附近,f1<0时定子空载电动势角度θe7角度应该在附近,误差小于认为角度估算正确;采用公式(8)计算每个Tset3间隔的定子空载电动势角度差值。采用公式(9)对2次定子空载电动势角度差值求取平均值,根据平均值计算5~7阶段的发电机频率f2;采用公式(10)计算到机侧变流器起动允许标识有效的时间间隔TSet4,同时定时器Tcnt清0,如图4中440、441所示;将Flux_Step设置为7,切换到阶段8,如图4中442;
阶段8(Flux_Step=7):定时器计时到TSet4后,机侧变流器允许起动标识ST_Per置1,采用公式(11)计算允许起动时刻的发电机角频率ω0,采用公式(12)计算发电机磁链角度初值θψ0,如图4中445、446所示;将Flux_Step设置为8,如图4中447。机侧变流器起动之前的初始磁链计算结束,准备机侧变流器起动。
步骤3:图3所示的磁链锁相环节320检测到机侧变流器允许起动标识ST_Per=1时开始磁链锁相,具体流程如图5所示。
输入的Q轴定子参考电流Isq乘以发电机角频率ω,再乘以发电机交轴电感Lq得到D轴电枢反应电压Usdar,如图5中500、501、502,输入的D轴定子控制电压Uscd加上D轴电枢反应电压Usdar得到D轴定子空载电动势Esd,如图5中511,对D轴定子空载电动势Esd求取积分得到发电机D轴磁链ψsd,如图5中512,发电机D轴磁链ψsd的积分初值为发电机磁链ψSet;
输入的Q轴定子控制电压Uscq减去D轴定子控制电压Uscd、DQ轴定子电压耦合系数Kcoff、发电机角频率ω三者的乘积,得到Q轴定子空载电动势Esq,如图5中520、521、 522;Q轴定子空载电动势Esq除以发电机D轴磁链ψsd得到发电机角频率ω,由于发电机角频率ω求解涉及到除法运算,在机侧变流器允许起动时刻无法得到Q轴定子空载电动势 Esq信息,此时将除法运算屏蔽,将输出的发电机角频率ω设置为初始磁链计算环节输出的允许起动时刻的发电机角频率ω0,如图5中523;对发电机角频率ω求取积分得到发电机磁链角度θψ,发电机磁链角度θψ的积分初值为初始磁链计算环节输出的发电机磁链角度初值θψ0,如图5中530;
步骤4:图3中所示的参考电流计算环节330在机侧变流器允许起动标识ST_Per置1后开始计算参考电流。在参考电流计算环节330中,输入的发电机角频率ω乘以发电机磁链ψSet得到定子参考电压UsRef;根据输入的D轴定子控制电压Uscd和Q轴定子控制电压 Uscq求取定子电压有效值UsRms,计算公式为
定子参考电压UsRef和定子电压有效值UsRms相减得到定子电压误差UsErr,定子电压误差UsErr经PI调节器输出D轴定子参考电流IsdRef;根据输入的有功功率给定值PGiven、D 轴定子控制电压Uscd、Q轴定子控制电压Uscq和D轴定子参考电流IsdRef计算Q轴定子参考电流IqRef,计算公式为
步骤5:图3中所示的电流闭环控制环节340在机侧变流器允许起动标识ST_Per置1后开始对定子电流进行闭环控制。在电流闭环控制环节340中,D轴定子参考电流IsdRef减去D轴定子电流Isd得到D轴定子电流误差IsdErr,D轴定子电流误差IsdErr经PI调节器输出D轴调节电压UsdPI;Q轴定子参考电流IsqRef减去Q轴定子电流Isq得到Q轴定子电流误差IsqErr,Q轴定子电流误差IsqErr经PI调节器输出Q轴调节电压UsqPI;采用公式(13) 求取D轴定子控制电压Uscd和Q轴定子控制电压Uscq;
D轴定子控制电压Uscd和Q轴定子控制电压Uscq经过基于发电机磁链角度θψ的2r/2s坐标变换后得到α轴定子控制电压Uscα和β轴定子控制电压Uscβ。
步骤6:电流闭环控制环节输出的α轴定子控制电压Uscα和β轴定子控制电压Uscβ经PWM调制策略后,输出12路PWM信号SPWM到主电路三相逆变桥,控制主电路中IGBT开通与关断,实现所述的矢量控制功能。
采用本发明所述的无定子电压测量矢量控制方法研制的5MW/3.3kV变流器起动加载波形如图6所示,(a)为三相定子电压(kV),(b)为三相定子电流(kA),(c)为有功功率(MW), (d)为三相电压放大图(kV),(e)为三相电流放大图(kA)。在图6中,在0.15s处接收到机侧变流器起动命令CMD_ST,在t1、t2、t3、t4时刻间隔Tset1输出ΔT0=50us的4个零矢量计算1~4阶段的发电机频率f1,在t5、t6、t7时刻间隔Tset3输出ΔT0=50us的3个零矢量计算4~7阶段的发电机频率f2,在t8时刻输出机侧变流器允许起动标识ST_Per,机侧变流器开始起动。5MW/3.3kV变流器额定峰值电流1360A,输出零矢量时定子电流峰值最大50A,机侧变流器起动时定子电流峰值最大120A,小于额定电流的10%。
Claims (3)
1.一种直驱风电变流器无定子电压测量矢量控制方法,其特征在于,所述的矢量控制方法包含定子电流计算环节、初始磁链计算环节、磁链锁相环节、参考电流计算环节和电流闭环控制环节五部分;所述的定子电流计算环节采集三相定子电流Isa、Isb、Isc,转化为α轴定子电流Isα、β轴定子电流Isβ、D轴定子电流Isd和Q轴定子电流Isq;所述的初始磁链计算环节通过输出7个△T0时间的零矢量,△T0为零矢量作用时间,分8个阶段确定机侧变流器允许起动时刻,确定允许起动时刻的发电机角频率ω0和发电机磁链角度初值θψ0,在允许起动时刻输出机侧变流器允许起动标识ST_Per;所述的磁链锁相环节在收到机侧变流器允许起动标识ST_Per后,在允许起动时刻的发电机角频率ω0和发电机磁链角度初值θψ0的基础上进行锁相,输出发电机角频率ω和发电机磁链角度θψ;所述的参考电流计算环节输出D轴定子参考电流IsdRef和Q轴定子参考电流IsqRef;所述的电流闭环控制环节采用PI调节器对DQ坐标系上的定子电流进行闭环控制,输出α轴定子控制电压Usca、β轴定子控制电压Uscβ、D轴定子控制电压Uscd和Q轴定子控制电压Uscq;所述的电流闭环控制环节输出的D轴定子控制电压Uscd和Q轴定子控制电压Uscq用于下一个控制周期的磁链锁相和参考电流计算,电流闭环控制环节输出的α轴定子控制电压Usca和β轴定子控制电压Uscβ经PWM调制策略后输出12路PWM调制信号SPWM,用于控制主电路三相逆变桥中的电力电子器件开通与关断,实现所述的矢量控制;
所述的磁链锁相环节在机侧变流器允许起动标识ST_Per=1时开始磁链锁相,输入的Q轴定子电流Isq乘以发电机角频率ω,再乘以发电机交轴电感Lq得到D轴电枢反应电压Usdar,输入的D轴定子控制电压Uscd加上D轴电枢反应电压Usdar得到D轴定子空载电动势Esd;对D轴定子空载电动势Esd求取积分得到发电机D轴磁链ψsd,发电机D轴磁链ψsd的积分初值为发电机磁链ψSet;输入的Q轴定子控制电压Uscq减去D轴定子控制电压Uscd、DQ轴定子电压耦合系数Kcoff与发电机角频率ω三者的乘积,得到Q轴定子空载电动势Esq;Q轴定子空载电动势Esq除以发电机D轴磁链ψsd得到发电机角频率ω,由于发电机角频率ω求解涉及到除法运算,在机侧变流器允许起动时刻无法得到Q轴定子空载电动势Esq信息,为此将除法运算屏蔽,将输出的发电机角频率ω设置为初始磁链计算环节输出的允许起动时刻的发电机角频率ω0;对发电机角频率ω求取积分得到发电机磁链角度θψ,发电机磁链角度θψ的积分初值为初始磁链计算环节输出的发电机磁链角度初值θψ0。
2.根据权利要求1所述的直驱风电变流器无定子电压测量矢量控制方法,其特征在于:所述的控制方法接收有功功率给定值PGiven,采集主电路中的三相定子电流Isa、Isb、Isc和直流母线电压Udc,采集信号不包含三相定子电压;所述矢量控制方法的定子电流计算环节采集发电机三相定子电流Isa、Isb、Isc,输入磁链锁相环节输出的发电机磁链角度θψ,输出α轴定子电流Isα、β轴定子电流Isβ、D轴定子电流Isd、Q轴定子电流Isq;初始磁链计算环节输入机侧变流器起动命令CMD_ST,输入定子电流计算环节输出的α轴定子电流Isα、β轴定子电流Isβ,输出机侧变流器允许起动标识ST_Per、允许起动时刻的发电机角频率ω0、发电机磁链角度初值θψ0;磁链锁相环节输入电流闭环控制环节输出的D轴定子控制电压Uscd和Q轴定子控制电压Uscq,输入定子电流计算环节输出的D轴定子电流Isd和Q轴定子电流Isq,输入初始磁链计算环节输出的机侧变流器允许起动标识ST_Per、允许起动时刻的发电机角频率ω0、发电机磁链角度初值θψ0,输出发电机磁链角度θψ、发电机角频率ω;参考电流计算环节输入有功功率给定值PGiven,输入磁链锁相环节输出的发电机角频率ω,输入电流闭环控制环节输出的D轴定子控制电压Uscd和Q轴定子控制电压Uscq,输出D轴定子参考电流IdRef、Q轴定子参考电流IqRef;电流闭环控制环节输入参考电流计算环节输出的D轴定子参考电流IdRef和Q轴定子参考电流IqRef,输入定子电流计算环节输出的D轴定子电流Isd和Q轴定子电流Isq,输入磁链锁相环节输出的发电机角频率ω和发电机磁链角度θψ,输出D轴定子控制电压Uscd、Q轴定子控制电压Uscq、α轴定子控制电压Uscα和β轴定子控制电压Uscβ。
3.根据权利要求1或2所述的直驱风电变流器无定子电压测量矢量控制方法,其特征在于:所述的初始磁链计算环节通过输出7个△T0时间的零矢量,分8个阶段确定机侧变流器允许起动的时刻,确定允许起动时刻的发电机角频率ω0和发电机磁链初始角度θψ0,在允许起动时刻输出机侧变流器允许起动标识ST_Per;
前4个阶段间隔TSet1输出4个△T0时间的零矢量,测量4个零矢量时刻的定子电流采用公式(2)计算定子空载电动势角度;
θe=a tan 2(-Isβ,-Isα) (2)
公式(2)中,Isα为α轴定子电流,Isβ为β轴定子电流,atan2为反正切运算符号,θe为定子空载电动势角度,1~4阶段零矢量作用时刻的空载电动势角度分别为θe1、θe2、θe3、θe4;
TSet1计算方法如公式(3)所示,其中fMax为设定的变流器输出最高频率;
根据1~4阶段零矢量作用时刻的空载电动势角度θe1、θe2、θe3、θe4,采用公式(4)计算TSet1时间间隔的定子空载电动势角度差值△θe1、△θe2、△θe3,
采用公式(5)对3个空载电动势角度差值求取平均值,平均值除以时间间隔TSet1,再除以2π得到1~4阶段的发电机频率f1;
在第4阶段采用公式(6)计算到下一个定子空载电动势角度为0的时间间隔TSet2;其中,f1为初始磁链计算1~4阶段的发电机频率,θe4为第4次输出零矢量时对应的定子空载电动势角度;
在5、6、7阶段分别在0、位置确认定子空载电动势角度,同时在第7阶段计算到下一个定子空载电动势角度为0的时间间隔TSet4,TSet4计算方法如公式(10)所示,其中,f2为初始磁链计算5~7阶段的发电机频率,θe7为第7次输出零矢量时对应的定子空载电动势角度;
在第8阶段定子空载电动势角度为0的时刻输出机侧变流器允许起动标识ST_Per,同时输出允许起动时刻的发电机角频率ω0和发电机磁链角度初值θψ0到磁链锁相环节,计算方法如公式(11)、(12)所示;
ω0=2πf2 (11)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010864025.1A CN112019118B (zh) | 2020-08-25 | 2020-08-25 | 一种直驱风电变流器无定子电压测量矢量控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010864025.1A CN112019118B (zh) | 2020-08-25 | 2020-08-25 | 一种直驱风电变流器无定子电压测量矢量控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112019118A CN112019118A (zh) | 2020-12-01 |
CN112019118B true CN112019118B (zh) | 2022-04-15 |
Family
ID=73503531
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010864025.1A Active CN112019118B (zh) | 2020-08-25 | 2020-08-25 | 一种直驱风电变流器无定子电压测量矢量控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112019118B (zh) |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0653000B2 (ja) * | 1984-09-14 | 1994-07-06 | いすゞ自動車株式会社 | リラクタンス発電機制御装置 |
JP4425193B2 (ja) * | 2005-08-16 | 2010-03-03 | 三洋電機株式会社 | モータの位置センサレス制御装置 |
CN101795007B (zh) * | 2010-03-18 | 2012-05-16 | 清华大学 | 基于同步化模型的风电双馈电机定子功率控制方法 |
CN101977006B (zh) * | 2010-10-28 | 2012-05-23 | 天津大学 | 电网故障情况下双馈风力发电机的控制方法 |
CN102437813B (zh) * | 2011-12-26 | 2014-04-09 | 中国东方电气集团有限公司 | 一种基于无速度传感器的永磁同步电机的转子角度、转速估计方法 |
JP5644820B2 (ja) * | 2012-08-17 | 2014-12-24 | 株式会社安川電機 | モータ制御装置 |
-
2020
- 2020-08-25 CN CN202010864025.1A patent/CN112019118B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112019118A (zh) | 2020-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9252601B2 (en) | Method for controlling a power converter in a wind turbine generator | |
CN102914740B (zh) | 快速辨识异步电机参数的方法 | |
JP3982232B2 (ja) | 同期発電機のセンサレス制御装置と制御方法 | |
EP2398140B1 (en) | Adjustable speed drive lifetime improvement method | |
CN107154761B (zh) | 发电机全控整流***和控制方法 | |
Ahmed et al. | A novel stand-alone induction generator system for AC and DC power applications | |
Jou et al. | Direct power control of a DFIG in wind turbines to improve dynamic responses | |
CN105119536A (zh) | 一种低成本电机驱动器拓扑及其控制方法 | |
CN108039843B (zh) | 一种dfig-dc***的定子频率控制方法 | |
CN110875696A (zh) | 用于大功率永磁直驱***的牵引变流器及其控制方法 | |
CN107947659B (zh) | 一种dfig-dc***的定子电流正弦化控制方法 | |
Xiang et al. | Switching frequency dynamic control for DFIG wind turbine performance improvement around synchronous speed | |
Pan et al. | Series compensated open-winding PM generator wind generation system | |
CN109494799B (zh) | 基于开绕组结构的永磁同步发电机直接并网***及其控制方法 | |
CN105099320B (zh) | 永磁直驱风力发电机输出有功功率控制方法及装置 | |
CN112019118B (zh) | 一种直驱风电变流器无定子电压测量矢量控制方法 | |
CN109088428B (zh) | 高电压穿越装置、方法,包括该装置的变流*** | |
WO2010018424A1 (en) | A power regulation scheme for a converter feeding into an asymmetrical grid | |
Lee et al. | Output voltage control of PWM inverters for stand-alone wind power generation systems using feedback linearization | |
CN109842316A (zh) | Dc-ac变换器的双前馈控制方法 | |
CN115579901A (zh) | 一种lcc-hvdc受端电网静态电压稳定性分析方法及*** | |
CN113300417B (zh) | 一种增强双馈风机同步稳定性的控制方法及*** | |
US11387757B2 (en) | Inductance-based estimation of rotor magnet temperature | |
Abolhassani et al. | Integrated doubly-fed electric alternator/active filter (IDEA), a viable power quality solution, for wind energy conversion systems | |
Fan | Mathematical modelling of grid connected fixed-pitch variable-speed permanent magnet synchronous generators for wind turbines |
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 |