双馈风力发电的矢量控制策略-电机及其系统分析与仿真

1双馈风力发电的矢量控制策略双馈电机在结构上与三相绕线式异步电机类似，其定子和转子均安放三相对称绕组，都可以与电网进行能量的交换。其定子绕组直接接入工频电网，转子绕组通过双馈变流器与电网连接。转子绕组连接于一个频率、相位、幅值均可调的三相电源激励，转子通入励磁电流产生旋转磁场，再加上转子的转速在气隙中产生一个同步旋转磁场。通过控制输入转子绕组的电流，不仅可以保证电机定子侧输出的电压和频率保持与电网电压一致，而且还可以调节双馈电机定子侧的功率因数。稳态运行时，双馈变流器根据所检测的电机转速调节流入双馈电机转子绕组的励磁电流频率以保证定转子电流所产生的旋转磁场在空间上保持相对静止，实现定子侧感应电势的频率与电网频率相同，以实现双馈型风力发电系统的变速恒频运行。双馈风力发电的系统原理图如图1所示。图1双馈风力发电系统原理图双馈变换器目前的多采用两电平双PWM变换器，其结构图如图2所示。图2两电平双PWM变换器21双馈发电机的数学模型1.1三相坐标轴系下数学模型定子绕组采用发电机惯例，定子电流流出为正，转子绕组采用电动机惯例，转子电流流入为正。则双馈发电机在三相静止坐标轴系下的模型为图3所示：图3三相坐标轴系下双馈发电机模型针对此模型可以得到三相坐标轴系下电压方程、磁链方程、运动方程和转矩方程为：电压方程：（1）转子侧电压方程：（2）定子侧电压方程：3ua1、ub1、uc1、ua2、ub2、uc2分别表示定转子电压，下标为1表示为定子侧，为2表示转子侧；ψa1、ψb1、ψc1、ψa2、ψb2、ψc2表示定、转子侧磁链；ia1、ib1、ic1、ia2、ib2、ic2为定子，转子相电流；R1、R2为定子，转子绕组的等效电阻。（3）磁链方程：其中Lm1是与定子绕组交链的最大互感磁通对应的定子电感；Lm2是与转子绕组交链的最大互感磁通对应的转子互感；Ll1，Ll2分别为定，转子漏电感；θr为转子的位置角。（4）运动方程：TL为风力机提供的拖动力矩；Dg为与转速成正比的阻转距阻尼系数；Kg为扭转弹性转矩系数；Jg为发电机的转动惯量。（5）转矩方程：41.2同步旋转坐标系dq下的数学模型将一个三相静止坐标系下的的矢量，通过变换用一个两相静止坐标系或两相旋转坐标系里的矢量表示，在变换时可采取幅值不变或者功率不变的原则。坐标关系如图4所示轴系为静止两相坐标，dq为同步旋转坐标轴系。图4坐标关系（1）定子绕组电压方程（2）转子绕组电压方程（3）转子磁链方程（4）定子磁链方程5（5）转矩方程（6）运动方程与三相静止坐标系下一致（7）定子有功功率、无功功率方程2转子侧矢量控制常用的转子侧矢量控制方式有转子磁链定向矢量控制、气隙磁场定向矢量控制、定子磁场定向矢量控制、电压定向矢量控制。其中定子磁场定向矢量控制使用较多。忽略定子绕组电阻后，则定子电压方程可以简化为u1表示相电压瞬时值，1表示定子磁链。定子磁链方向为同步旋转坐标系d轴方向时，则定子磁链具有如下关系此外对于基于定子电压和转子转速的U—ω定子磁链检测法，有6对于基于定子电压和同步转速的U—ω1法，在取定子磁链定向后，若忽略定子的电阻，则定子电压矢量和定子磁链矢量之间相位差90°，幅值相差一个同步速ω1的倍数。磁基于场定向的交流励磁变速恒频风力发电机的功率调节及转矩调节则可得其中7双馈发电机的功率调节包括有功功率P1和无功功率Q1的调节。交流励磁电机转子侧电流在q轴上的分量iq2与定子侧输出的有功功率P1成线性比例关系，交流励磁电机转子侧电流在d轴上的分量id2与定子侧输出的无功功率Q1成线性比例关系，即可以通过调节交流励磁电机转子侧电流（id2，iq2）来直接调节交流励磁电机定子侧输出的有功功率、无功功率（P1、Q1）。双馈励磁电机dq轴电压（ud2，uq2）和dq轴电流（id2，iq2）也存在线性关系，那么此刻只需要通过调节转子侧电压值（ud2、uq2）即可以间接的调节系统所产生的有功功率P1、无功功率Q1，并使之达到给定值P1*和Q1*。而转子侧电压值（ud2、uq2）经过旋转变换后得到相应的发电机转子三相电压控制指令（ua2*、ub2*、uc2*），并将它们作为调制波与三角载波比较以产生SPWM脉冲去控制主电路开关管IGBT的通断，实现有功功率P1，无功功率Q1的独立调节控制。交流励磁电机的电磁转矩Te与其转子侧电流的q轴上的分量iq2成线性比例关系，即在电机的运行过程中，可以通过对电机转子侧电流iq2的调节就可以达到对电机电磁转矩直接控制的目的。交流励磁电机转子dq轴电压（ud2，uq2）和转子dq轴电流（id2，iq2）也存在线性关系，那么此刻只需要调节转子侧电压值（ud2，uq2）即可以调节系统所产生的电磁转矩Te、无功功率Qs，并使之达到给定值Te*和Qs*。转子侧电压值（ud2，uq2）经过旋转变换后得到相应的发电机转子三相电压控制指令（ua2*，ub2*，uc2*），并将它们作为调制波与三角载波比较以产生SPWM脉冲去控制主电路开关管IGBT的通断，实现变速恒频、电磁转矩Te、无功功率Qs的独立调节控制。83网侧变换器控制网侧变换器主电路如图5所示图5网侧变换器主电路稳定工作状态时，母线上的直流电压恒定，网侧变换器的三相桥按正弦调制规律驱动。调节网侧变换器输出交流电压的幅值和相位就能控制电感电流的大小以及电流与电网电压的相位角，从而使该变换器运行于单位功率因数整流运行、单位功率因数逆变运行和非单位功率因数运行。设三相电网电压平衡，则Sk（k＝a、b、c）分别为三相桥臂的开关函数。其中：Sk=1：对应桥臂上管导通，下管关断；Sk=0：对应桥臂下管导通，上管关断。，9在同步旋转坐标轴系下，有Sd、Sq是开关函数Sk（k＝1、2、3）变换到d-q坐标系中的d、q轴相应的开关函数。设则：令则：上式中Ed、Eq为电网电压扰动的影响。104双馈风力发电矢量控制Matlab模型双馈风力发电的整体Matlab模型如图6所示，网侧变换器控制的Matlab模型如图7所示，转子控制Matlab模型如图8所示图6双馈风力发电整体Matlab模型9MWWindFarm(6x1.5MW)WindFarm-DFIGDetailedModelDiscrete,Ts=5e-006s.Scope-K-MWABCNGroundingTransformer[Q_pu][wr][Vdc][P_pu]Iabc_B575Vabc_B575[wr][Vdc]Iabc_B25Vabc_B25[P_pu][Q_pu]Wind(m/s)Qref_pumABCDFIGWindTurbine015ABCabcB575(575V)ABCabcB25(25kV)ABCabcB120(120kV)ABCABC30kmline3.3ohmsABCABC2500MVAX0/X1=3ABCabc25kV/575V6*1.75MVAABCabc120kV/25kV47MVANABC120kVVabc_B575(pu)Iabc_B575(pu)P(MW)Q(Mvar)Vdc(V)wr(pu)Vabc_B25(pu)Iabc_B25(pu)Vdc_Vwr_pu(IGspeed)P_puQ_pu11图7网侧变换器控制的Matlab模型WindTurbinecontrols-Grid-sideconvertercontrolsystemWhereVnom=RMSph-phnominalvoltagem=Vnom*2*sqrt(2/3)/VdcRequiredvalueformodulationindexmtoobtain1pugeneratedvoltagebytheconverter:VdcregulatorCurrentregulatorModulationindex&phaseUctrlgridconverter1Uctrl_grid_convdq0sin_cosabcdq0_to_abcTransformationVdc_nomVdc_ref(V)0cossinPolartoCartesianVnom*2*sqrt(2/3)K[Vdc][w][id_ref][iq_ref][Vdc][id_ref][w][iq_ref][iq_ref][w][id_ref]PIPIDemuxDemuxDemuxDemuxR_RLR_RLL_RLL_RLCartesiantoPolar1Avoiddivisionbyzero0-Mod_index_max_grid-K--pu6Theta_PLL5Vdc4Vdqs3w2Idq_gc1Iq_refIq_refId_refId_refIq_refvq'vd'VqIdIqIqVdmVdq0Vd_ctrl_grid_convVq_ctrl_grid_conv12图8转子控制的Matlab模型varregulatorVoltageregulatorIqr_reflimitIdrrefIqrrefWindTurbinecontrols-Rotor-sideconvertercontrolsystemElectromagnetictorque(Tem)controllerCurrentregulatorWhereVnom_r=RMSph-phnominalvoltagem=Vnom_r*2*sqrt(2/3)/VdcRequiredvalueformodulationindexmtoobtain1pugeneratedvoltagebytheconverter:Modulationindex&phaseNote:Thevar/voltcontrolsystemisbasedonthereport:DynamicModelingofGE1.5and3.6WindTurbine-Generators;Preparedby:NicholasW.Miller,WilliamW.Price,JuanJ.Sanchez-Gasca;October27,2003,Version3.0GE-PowerSystemsEnergyConsultingCopyright2002GeneralElectricCompany,U.S.A.Refertofigure4.3onpage4.4ofthisreporttoseethecontrolsystem.(Lls+Lm)Lm/Ls(Lm+Lls)/RsRs^2(Lm+Lls)^2LmLm*(Lm+Lls)/RsLlr+LmRsLm/RsLlr+LmVnom_r*2*sqrt(2/3)K1Irotor_max^2MaxMinInKiHoldResetResetValueOutDiscrete-Time(Trapezoidal)IntegratorINT[Ids_estim][Iqr_ref][Vdqs][Iqs_estim][Idr_ref][Idqs][w_wr][w]-K-[Ids_estim][w_wr][Idr_ref][w_wr][Iqs_estim][w_wr][w][w][w][Idr_ref][Iqr_ref][Iqr_ref][w][w][w][Idr_ref][Vdqs][Iqr_ref][w_wr][Idr_ref][Idqs]KTsz-1KTsz-1xoPIDemuxDemuxDemuxboolean1Ki_volt01CartesiantoPolarAvo