《武汉工程大学学报》 2018年01期
93-97
出版日期：2018-02-25
ISSN:1674-2869
CN:42-1779/TQ
直驱永磁风电机组轴系扭振的非线性控制
风能的利用近年来得到了巨大的发展，但风电大规模并网给电力系统的稳定运行带来了新的挑战，并网运行后电网的扰动或风速的变化会导致机网扭振，对风电机组的稳定运行和设备安全产生很大的影响［1-5］。在大容量风电场中，严重的扭振现象会造成风电机组损耗加剧甚至损坏，影响电力系统的安全运行［6-7］。因此，深入研究风电机组的动态特性及其与电网的相互作用，以及风机轴系的扭振问题是保证机组安全运行的重要课题。直驱永磁风电传动系统具有“柔性”，受到电网扰动时轴系出现的低频振荡，使传动轴出现应力疲劳［8-12］。国内学者已经探讨了风机的稳定性运行和轴系扭振的问题［13-16］，已有学者对非线性控制理论应用于电力系统稳定分析进行了深入的探讨［17-19］。本文基于直驱永磁风机（direct-driven permanent magnet synchronous generators，D-PMSG ）轴系扭振，建立了风机发电机组的模型，采用非线性控制策略实现协调控制来抑制轴系间扭矩振荡，仿真结果验证了该控制策略抑制轴系扭振的有效性。1 直驱永磁风电系统模型1.1 风电系统模型结构图风电机组系统采用如图1所示的系统模型，由空气动力学可知，风力机产生的机械功率为风力机从空气中捕获的风能，由（1）式表示：[Pw=0.5ρπr2V3wCp(β,γ),] [γ=rωwVw] （1）式（1）中[ρ]为空气的密度，[r]为风力机叶片半径，[Vw]为风速，[Cp]为风能利用系数，是[γ]和[β]的非线性函数，[β]为叶片的桨距角，[γ]为风力机叶尖速比，[ωw]为风力机的转速，[P=Tωw]。 1.2 系统数学模型1.2.1 机械轴系模型 风电机的轴系通常等效为两个或三个质量块［4］，这里采用两质块，数学模型为式（2）：[2Hwdωwdt=Tw-Kθ-Dωw-ωg2Hgdωgdt=-Te+Kθ+Dωw-ωgdθwgdt=ω0ωw-ωg ] （2）式（2）中[D]为阻尼系数，[K]为刚性系数，[Hw]、[Hg]为风机、发电机惯性时间常数，[θw]，[θg]分别为相对角位移，[ωw]、[ωg]为转子角速度，[Tw]，[Te] 分别为机械转矩、电磁转矩1.2.2 桨距角模型 风速通过叶片桨距角[β]的控制来实现，如下方程表示桨距角模型［8］： [dβdt=1Ts(β0-β)] （3）式（3）中[Ts]是桨距控制系统的惯性时间常数，[β0]是桨距角初始值。1.2.3 D-PMSG 模型 设dq坐标系以同步速度旋转，且q轴超前于d轴，则在dq轴坐标下的D-PMSG模型为式（4）：[di1ddt=-RaL1i1d+npωmi1q+1L1udcm1cosθ1di1qdt=-RaL1i1q-npωmi1d-1L1λ0npωm+1L1udcm1sinθ1] （4）式（4）中[u1d]、[u1q]分别为d、q轴电压，[i1d]、[i1q]分别为d、q 轴电流；[L1d]、[L1q]分别为定子直轴电感和交轴电感；[λ0]为永磁体磁链；[np]为转子的极对数；[ωm]及[Ra]分别为定子电角速度和电阻。根据能量守恒，不计变流器的能量损失，直流环节表示为式（5）：[dudcdt=32Ci1dm1cosθ1-32Ci1qm1sinθ1-32Ci2dm2cosθ2+32Ci2qm2sinθ2] （5）式（5）中，[m1,θ1,m2,θ2]分别为变换器直流环节两侧变比和导通角。1.2.4 电网侧数学模型 由结构图，取d轴与a相电压矢量重合，d轴超前q轴90°电角度，则网侧模型为［13］式（6）：[di2ddt=-RL2i2d+npωgi2q+1L2udcm2cosθ2-UdL2di2qdt=-RL2i2q-npωgi2d+1L2udcm2sinθ2] （6）式（6）中[u2d]、[u2q]分别为d、q轴电压，[i2d]、[i2q]分别为d、q轴电流；[UdL2]为电网电压，[ωg]为电网角速度，[L2]、R分别为连接电感及等值电阻。综合式（2）~式（6）构成了分析直驱永磁风电机轴系扭振的数学模型，可以写出式（7）形式［18］： [x=f(x)+j=19gj(x)u] （7）其中 [x=x1,x2,x3,x4,x5,x6,x7,x8,x9T=ωw,ωg,θwg,β,i1d,i1q,udc,i2d,i2qT] 根据对轴系扭振模型和控制目标的分析，确定3个输入量分别为：[u1=Tw,u2=m1cosθ1,u3=m2sinθ2]。2 直驱永磁风电系统轴系扭振的非线性控制 根据对风机控制做分析，将桨距角，机侧、网侧变流器作为控制对象。 考虑非线性仿射控制系统［16］：[x(t)=f(x)+g1(x)u1++gm(x)umy1(t)=h1(x)ym(t)=hm(x)] （8）式（8）中[x∈Rn]为状态量，[f(x)=[f1(x)fn(x)]T∈Rn] 、[gj(x)=[g1j(x)?gnj(x)]T∈Rn]为光滑向量，[u1um]为控制标量；[h1(x)hm(x)]为输出函数；[y1ym]为输出标量。 为保证协调控制的效果，输出函数选为多个状态量的线性组合［19］：桨距角[h1(x)=c11ωw+c12ωg+c13θwg+c14i1d]，机侧变流器[h2(x)=c23θwg+c24β+c26i1q+c27udc]，网侧变流器[h3(x)=c32ωg+c37udc+c39i2q]。将三个输出函数写成矩阵的形式：[y=[h1(x)h2(x)h3(x)]T=Cx] （9）[C=c11c12c13c140000000c23c240c26c27000c320000c370c39]根据非线性控制理论，先计算矩阵［22］：[B(x)=Lg1L0fh1(x)Lg2L0fh1(x)Lg3L0fh1(x)Lg1L1fh2(x)Lg2L1fh2(x)Lg3L1fh2(x)Lg1L2fh3(x)Lg2L2fh3(x)Lg3L2fh2(x)]（10）经过计算，系统的相对阶分别为[r1+r2+r3=5