[0053] 本发明提出的一种风力发电机转子侧变换器的功率控制方法，该方法应用于风力发电机，通过本发明公开的功率控制方法，可以高效稳定的实现对风力发电机转子侧变换器的功率控制。
[0054] 在本实施例中，无论是该功率控制方法还是功率控制系统都是基于风力发电机以实现的。但本发明阐述的重点为本功率控制方法的执行步骤以及实现该方法的功率控制系统的组成。由此在本实施例中，针对风力发电机的结构仅简略描述了和本发明有关的几个部分，而有关风力发电机结构的描述也不应作为对本发明保护范围的限制。
[0055] 另外，本实施例中出现的“连接”一词，应做广义理解。在没有特定说明的情况下，应根据方案的实际情况，可以理解为有线连接，也可以作为无线连接；在涉及到数据交换时，还可理解为通信连接以达到模块间数据传输的目的。
[0056] 首先，下文将结合附图对本实施例公开的一种风力发电机转子侧变换器的功率控制方法做详细描述。该功率控制方法基于功率控制系统以实现，功率控制系统包括电信号采集装置及控制单元，该系统的详细结构将在说明完本方法之后在做具体描述。
[0057] 请结合图1、图2，图2中示出了控制框图，图中定子101、转子102、转子侧变换器2作为简化的风力发电机的构成，现有技术风力发电机上还包括风叶、齿轮箱等必要部件，因其并不参与到本发明公开的技术方案中，因此没有示出。转子侧变换器2为逆变器，其包含3组(共6个)开关管，开关管组成共有8种驱动状态，本发明的目的在于通过下文描述的方法得到最优的驱动状态，实现风力发电机转子侧变换器的功率控制。
[0058] 本功率控制方法的具体步骤包括：
[0059] 步骤A：对功率控制系统设定离散采样周期为Ts及采样时刻k，并初始化k＝0；
[0060] 步骤B：电信号采集装置采集第k时刻的定子101电流矢量信号Is(k)和转子102电流矢量信号Ir(k)，经转子磁链和定子磁链计算模块6，得到第k时刻的转子磁链矢量信号Ψr(k)和第k时刻的定子磁链矢量信号Ψs(k)。
[0061] 在本实施例中电信号采集装置可以为各种型号的电压传感器或电流传感器，可通过电性连接，亦或是其他非接触式测量的感应器，如：霍尔式传感器。
[0062] 步骤C：基于第k时刻的转子磁链矢量信号Ψr(k)、第k时刻的定子磁链矢量信号Ψs(k)、经过定子有功功率和定子无功功率计算模块7，得到风力发电机第k时刻的定子有功功率Ps(k)和第k时刻的定子无功功率Qs(k)。
[0063] 步骤D：电信号采集装置采集转子侧变换器2直流侧第k时刻的直流母线电压信号Udc(k)，定子有功功率和定子无功功率预测模块9依据直流母线电压信号Udc(k)生成第k时s1s2s3刻的8个转子电压矢量参考值Ur (k)，
[0064] Urs1s2s3(k)＝(2/3)(s1+s2a+s3a2)Udc(k)
[0065] 其中， j为虚部单位。s1＝1或0，s2＝1或0，s3＝1或0，s1s2s3可以为“000；001；010；011；100；101；110；111”8个组合，代表转子侧变换器2开关管8s1s2s3
个不同的驱动信号状态。由此，8个状态下的Ur (k)的值分别为：
[0066] U000r(k)＝0；
[0067] U001r(k)＝(2/3)ej4π/3Udc(k)；
[0068] U010r(k)＝(2/3)ej2π/3Udc(k)；
[0069] U011r(k)＝(2/3)ejπUdc(k)；
[0070] U100r(k)＝2/3Udc(k)；
[0071] U101r(k)＝(2/3)ej5π/3Udc(k)；
[0072] U110r(k)＝(2/3)ejπ/3Udc(k)；
[0073] U111r(k)＝0。
[0074] 基于第k时刻的8个转子电压矢量参考值Urs1s2s3(k)、第k时刻的定子有功功率Ps(k)、第k时刻的定子无功功率Qs(k)、第k时刻的转子磁链矢量信号Ψr(k)和第k时刻的定子磁链矢量信号Ψs(k)，经过定子有功功率和定子无功功率预测模块9，得到第k+1时刻的8个s1s2s3 s1s2s3定子有功功率预测值Ps (k+1)和第k+1时刻的8个定子无功功率预测值Qs (k+1)。
[0075] 步骤E：在功率控制系统中设定风力发电机的有功功率给定指令Ps，风力发电机的无功功率给定指令Qs，有功功率代价函数约束参数CP，无功功率代价函数约束参数CQ，给定代价函数约束参数的调整参数a1和a2。
[0076] 步骤F：经代价函数约束条件判断模块10，判断并筛选出符合约束条件的转子电压s1s2s3矢量参考值Ur (k)，代价函数约束条件判断模块10的约束判断条件为：
[0077] 如果|Ps‑Pss1s2s3(k+1)|≤CP，且|Qs‑Qss1s2s3(k+1)|≤CQ，则风力发电机第k时刻的转s1s2s3 s1s2s3子电压矢量参考值Ur (k)满足代价函数约束条件，当8个电压矢量参考值Ur (k)均不s1s2s3
满足约束条件时，执行步骤G；当有1个及1个以上的电压矢量参考值Ur (k)满足约束条件时，执行步骤H。
[0078] 步骤G：将CP+a1赋值给CP，将CQ+a2赋值给CQ，并返回步骤F执行。
[0079] 步骤H：基于满足约束条件的转子电压矢量参考值Urs1s2s3(k)，经优化模块11，得到最优第k+1时刻的转子电压矢量参考值Vr(k)。即：在所有满足约束条件的转子电压矢量参s1s2s3考值Ur (k)中优化得出其中的最优值，并将最优值作为第k+1时刻的转子电压矢量参考s1s2s3
值Vr(k)，优化模块的优化方法为：如果其中的一个转子电压矢量参考值Ur (k)，使得两s1s2s3 s1s2s3
个误差值“Ps与Ps (k+1)之间误差”以及“Qs与Qs (k+1)之间误差”之和最小，则该转s1s2s3
子电压矢量参考值Ur (k)为最优第k+1时刻的转子电压矢量参考值Vr(k)。
[0080] 步骤I：基于最优第k+1时刻的转子电压矢量参考值Vr(k)，经SVPWM模块12，得到风力发电机转子侧变换器开关管第k时刻的驱动信号，进而驱动风力发电机转子侧变换器2工作。本步骤的目的在于：Vr(k)作为一个矢量参考值，变换器2是无法直接识别Vr(k)并根据Vr(k)更改开关管驱动状态的。SVPWM模块12的作用就是将Vr(k)转化成变换器2的驱动信号。
[0081] 步骤J：将k+1赋值给k,等待下一采样周期,返回步骤B执行。
[0082] 进一步地，在步骤B中，转子磁链和定子磁链计算模块6的计算方式为：
[0083] Ψr(k)＝LrIr(k)+LmIs(k)，
[0084] Ψs(k)＝LsIs(k)+LmIr(k)，
[0085] 其中：Lr为转子绕组自感；Lm为定转子互感，Ls为定子绕组自感。
[0086] 进一步地，在步骤C中，定子有功功率和定子无功功率计算模块7的计算方式为：
[0087] Ps(k)＝1.5ωsσLmIm{[Ψr(k)]*Ψs(k)}
[0088] Qs(k)＝1.5ωsσ[Lr|Ψs(k)|2‑LmRe{[Ψr(k)]*Ψs(k)}]
[0089] 其中：ωs为同步旋转角频率；σ为漏磁系数，且σ＝1/[LsLr‑(Lm)2]；*为向量的共轭。
[0090] 进一步地，在步骤D中，经过所述定子有功功率和定子无功功率预测模块9，得到第s1s2s3k+1时刻的8个定子有功功率预测值Ps (k+1)和第k+1时刻的8个定子无功功率预测值s1s2s3
Qs (k+1)。
[0091] Pss1s2s3(k+1)的计算公式为：
[0092] Pss1s2s3(k+1)＝Ps(k)+1.5TsωsσLm[Im{[Urs1s2s3(k)]*Ψs(k)}+ωgRe{Ψs(k)[Ψr*(k)]}] (I)
[0093] 由此，具体地8个驱动状态下的有功功率预测值为：
[0094] Ps000(k+1)＝Ps(k)+1.5TsωsσLm[Im{[U000r(k)]*Ψs(k)}+ωgRe{Ψs(k)[Ψr*(k)]}]，
[0095] Ps001(k+1)＝Ps(k)+1.5TsωsσLm[Im{[U001r(k)]*Ψs(k)}+ωgRe{Ψs(k)[Ψr*(k)]}]，
[0096] Ps010(k+1)＝Ps(k)+1.5TsωsσLm[Im{[U010r(k)]*Ψs(k)}+ωgRe{Ψs(k)[Ψr*(k)]}]，
[0097] Ps011(k+1)＝Ps(k)+1.5TsωsσLm[Im{[U011r(k)]*Ψs(k)}+ωgRe{Ψs(k)[Ψr*(k)]}]，
[0098] Ps100(k+1)＝Ps(k)+1.5TsωsσLm[Im{[U100r(k)]*Ψs(k)}+ωgRe{Ψs(k)[Ψr*(k)]}]，
[0099] Ps101(k+1)＝Ps(k)+1.5TsωsσLm[Im{[U101r(k)]*Ψs(k)}+ωgRe{Ψs(k)[Ψr*(k)]}]，
[0100] Ps110(k+1)＝Ps(k)+1.5TsωsσLm[Im{[U110r(k)]*Ψs(k)}+ωgRe{Ψs(k)[Ψr*(k)]}]，
[0101] Ps111(k+1)＝Ps(k)+1.5TsωsσLm[Im{[U111r(k)]*Ψs(k)}+ωgRe{Ψs(k)[Ψr*(k)]}]，
[0102] Qss1s2s3(k+1)的计算公式为：
[0103] Qss1s2s3(k+1)＝Qs(k)‑1.5TsωsσLm[Re{[Urs1s2s3(k)]*Ψs(k)}‑ωgRe{[Ψr(k)]*Ψs(k)}] (II)
[0104] 由此，具体地8个驱动状态下的有功功率预测值为：
[0105] Qs000(k+1)＝Qs(k)‑1.5TsωsσLm[Re{[U000r(k)]*Ψs(k)}‑ωgRe{[Ψr(k)]*Ψs(k)}]，
[0106] Qs001(k+1)＝Qs(k)‑1.5TsωsσLm[Re{[U001r(k)]*Ψs(k)}‑ωgRe{[Ψr(k)]*Ψs(k)}]，
[0107] Qs010(k+1)＝Qs(k)‑1.5TsωsσLm[Re{[U010r(k)]*Ψs(k)}‑ωgRe{[Ψr(k)]*Ψs(k)}]，
[0108] Qs011(k+1)＝Qs(k)‑1.5TsωsσLm[Re{[U011r(k)]*Ψs(k)}‑ωgRe{[Ψr(k)]*Ψs(k)}]，
[0109] Qs100(k+1)＝Qs(k)‑1.5TsωsσLm[Re{[U100r(k)]*Ψs(k)}‑ωgRe{[Ψr(k)]*Ψs(k)}]，
[0110] Qs101(k+1)＝Qs(k)‑1.5TsωsσLm[Re{[U101r(k)]*Ψs(k)}‑ωgRe{[Ψr(k)]*Ψs(k)}]，
[0111] Qs110(k+1)＝Qs(k)‑1.5TsωsσLm[Re{[U110r(k)]*Ψs(k)}‑ωgRe{[Ψr(k)]*Ψs(k)}]，
[0112] Qs111(k+1)＝Qs(k)‑1.5TsωsσLm[Re{[U111r(k)]*Ψs(k)}‑ωgRe{[Ψr(k)]*Ψs(k)}]；
[0113] 在上述公式(I)和(II)中，ωg为风力发电机角频率。
[0114] 进一步地，在步骤H中，优化模块11的优化方法为：如果其中的一个转子电压矢量s1s2s3 s1s2s3 s1s2s3参考值Ur (k)，使得|Ps‑Ps (k+1)|+|Qs‑Qs (k+1)|的值最小，则该转子电压矢量s1s2s3
参考值Ur (k)为最优第k+1时刻的转子电压矢量参考值Vr(k)；
[0115] 结合步骤E，该优化方法的含义为，在若干满足条件的驱动状态下，实时的有功功s1s2s3 s1s2s3率Ps (k+1)与设定的有功功率Ps之间的误差为|Ps‑Ps (k+1)|，实时的无功功率s1s2s3 s1s2s3
Qs (k+1)与设定的无功功率Qs之间的误差为|Qs‑Qs (k+1)|，选取其中一个驱动状态使得这两个误差的总和是最小的，则该驱动状态即为最优的驱动状态。
[0116] 显然，其采用绝对值的算法是为了将误差值取正数。由此还有其他若干的算法，另s1s2s3 s1s2s3外的一种优化方法为：如果其中的一个转子电压矢量参考值Ur (k)，使得[Ps‑Ps (k+
2 s1s2s3 2 s1s2s3
1)]+[Qs‑Qs (k+1)]的值最小，则该转子电压矢量参考值Ur (k)为最优第k+1时刻的转子电压矢量参考值Vr(k)。本算法采用误差值的平方以保证误差值为正数。诸如此类的方法还有很多，在此不做赘述。
[0117] 本实施例提供了一组仿真数据以更清楚的展示本方法的效果，仿真的结果请参照图3、图4以及图5。
[0118] 本实施例在仿真环境下设置的基本参数如下：
[0119] 风力发电机1的额定功率Pn为1500kW，额定频率f为50Hz，电网3线电压为380V，极对数p为2，定子101电阻Rs为0.46Ω，转子102电阻Rr为0.18Ω，定子绕组自感Ls为0.006H，转子绕组自感Lr为0.212H，定转子互感Lm为0.033H。上述为风力发电机1相关的参数设置，符合实际情况即可。
[0120] 本实施例公开的仿真实验的运行时间共2s，图3、图4、图5的横坐标均表示仿真实验的时间。
[0121] 仿真环境下的风速如图3所示，风速在7.4m/s至15.4m/s之间变化，平均风速为11.3m/s。
[0122] 步骤A中的离散采样周期Ts设定为0.00001s；
[0123] 步骤E中的有功功率代价函数约束参数CP、无功功率代价函数约束参数CQ可自由设定为一个常量，本实施例中CP＝16.5kW；CQ＝16.5kVar。给定代价函数约束参数的调整参数a1和a2可自由设定，本实施例中a1＝0.5kW，a2＝0.5kVar。
[0124] 另一组需要设定的数据为有功功率给定指令Ps和无功功率给定指令Qs，两个指令5 5
可以被设定为任意的常量。如:Ps＝‑10×10W；Qs＝2.5×10 Var。但是，为了观察风力发电机1的动态和稳态响应，优选地，风力发电机的有功功率给定指令Ps和风力发电机的无功功率给定指令Qs为阶跃式时变函数。
[0125] 具体地，本实施例中，Ps和Qs的阶跃式时变函数是这样设定的：
[0126] 风力发电机1的有功功率给定指令Ps开始被设定为0W；在时间为0.4s时，有功功率5 5
给定指令Ps从0W跳变到‑5×10 W；在时间为1s时，有功功率给定指令Ps从‑5×10 W跳变到‑
5 5 5
15×10W；在时间为1.6s时，有功功率给定指令Ps从‑15×10W跳变到‑5×10W。
[0127] 风力发电机1的无功功率给定指令Qs开始被设定为‑5×105Var；在时间为0.6s时，5
无功功率给定指令Qs从‑5×10Var跳变到0Var；在时间为1s时，无功功率给定指令Qs从0Var
5 5
跳变到2.5×10 Var；在时间为1.6s时，无功功率给定指令Qs从2.5×10 Var跳变到5×
5
10Var。
[0128] 图4、图5给出了本实施例中仿真实验的结果，可以看出，定子实时有功功率(图4)和定子实时无功功率(图5)虽在图形上有微小幅度的波动，但总体上是沿有功功率给定指令Ps和无功功率给定指令Qs设定的阶跃式时变函数而变化的。
[0129] 综上，本功率控制方法可高效地实现对风力发电机转子侧变换器的功率控制，且具有较好的动态响应。
[0130] 本实施例还公开了一种风力发电机转子侧变换器的功率控制系统，该功率控制系统用以实现上述功率控制方法。具体地，功率控制系统包括电信号采集装置及控制单元。
[0131] 电信号采集装置分别与风力发电机的定子101、转子102及转子侧变换器2相连接，以采集定子101、转子102及转子侧变换器2的电压和/或电流信号。
[0132] 控制单元包括依次连接的转子磁链和定子磁链计算模块6、定子有功功率和定子无功功率计算模块7、定子有功功率和定子无功功率预测模块9、代价函数约束条件判断模块10、优化模块11、SVPWM模块12。控制单元中各模块之间采用通信连接以传输数据。控制单元与电信号采集装置之间亦采用通信连接。
[0133] 该系统详细的连接方式请参照图2，以下给以说明：电信号采集装置子电流矢量信号Is(k)和转子电流矢量信号Ir(k)并发送至转子磁链和定子磁链计算模块6，经计算得出转子磁链矢量信号Ψr(k)和定子磁链矢量信号Ψs(k)。
[0134] 转子磁链和定子磁链计算模块6将Ψr(k)和Ψs(k)分别发送至定子有功功率和定子无功功率计算模块7和定子有功功率和定子无功功率预测模块9，定子有功功率和定子无功功率计算模块7根据Ψr(k)和Ψs(k)计算得到第k时刻的定子有功功率Ps(k)和第k时刻的定子无功功率Qs(k)并将Ps(k)、Qs(k)发送至定子有功功率和定子无功功率预测模块9。
[0135] 电信号采集装置采集转子侧变换器2的电压Udc(k)并发送至定子有功功率和定子s1s2s3无功功率预测模块9，定子有功功率和定子无功功率预测模块9计算得出Ps (k+1)、s1s2s3 s1s2s3
Qs (k+1)，同时生成对应驱动状态下的转子电压矢量参考值Ur (k)。并将上述所得数据发送至代价函数约束条件判断模块10。
[0136] 代价函数约束条件判断模块10；判断并筛选出符合约束条件的转子电压矢量参考s1s2s3 s1s2s3值Ur (k)；并将满足条件的若干Ur (k)发送至优化模块11；优化模块11优化得出最优s1s2s3
的Ur (k)作为第k+1时刻的转子电压矢量参考值Vr(k)，并将Vr(k)发送至SVPWM模块12，SVPWM模块12将Vr(k)转化成对应的转子侧变换器2的驱动信号。
[0137] 上述为各模块件的信息传输关系，信息传输的过程以代表各模块之间的连接关系。
[0138] 进一步地，电信号采集装置包括：用于采集风力发电机定子电流量的定子电流传感器4；用于采集风力发电机转子电流量的转子电流传感器5；以及用于采集转子侧变换器电压量的直流母线电压传感器8。
[0139] 进一步地，定子电流传感器4、转子电流传感器5为霍尔电流传感器，直流母线电压传感器8为霍尔电压传感器，采用非接触式的传感器可有效简化电路连接。
[0140] 本实施例仅为本发明一种较优的实施方式，尽管上述已较为详尽的给出了具体的技术方案，但不应构成对本发明保护范围的限制。尤其是仿真过程中设置的一些参数，不应理解为对本发明技术方案的限制。方法权利要求中的部分步骤，在可实现的基础上存在顺序变换的可能。本领域技术人员在可以实现的基础上，使用本发明的原理，但改变了部分步骤的顺序也应属于本发明的保护范围。
