[发明专利]基于相关分析求取风电场稳态等效风速与发电功率的方法

权利要求书

1.一种基于相关分析求取风电场稳态等效风速与发电功率的方法，其特征在于具有以下步骤：

1).风电场运行数据的预处理：

首先对从风电场直接采集到的“生数据”进行预处理，进行膨胀和腐蚀两种基本运算；设f(m)为定义在Zⁿ上的离散函数，即f：Zⁿ→Z，结构元素B为Zⁿ上的有限子集，即B关于原点的对称集合为B^s＝{-b：b∈B}，B关于点m平移集合为B_m＝{b+m：b∈B，m∈Zⁿ}，则f(m)关于B的膨胀和腐蚀运算分别为

(f⊕Bs)(m)=maxb∈Bx{f(b)}]]>

(fΘBs)(m)=minb∈Bx{f(b)}]]>

由腐蚀和膨胀运算可得到f(m)关于B的开运算和闭运算分别为

(fob)=(fΘb)⊕b]]>

(f·b)=(f⊕b)Θb]]>

这里符号ο和·分别代表开运算和闭运算；

选用扁平结构元素，在其定义域内取零；直接采用开、闭运算的均值构成混合滤波器，然后通过选取合适宽度和幅值的结构元素，该混合滤波器可以对连续干扰及随机背景噪声干扰取得很好的抑制效果，满足分析要求；

2).求取等效风速：

(1)形成风速矩阵

设风电场内运行风机的数量为n，取出某一主风向下m个时刻滤波后每台风机感受到的风速，根据统计学要求，m≥(2～3)n，形成风速矩阵如下

V=v11v12···v1nv21v22···v2nMvm1vm2···vmn]]>

对风速矩阵标准化；风速样本数据标准化的实质是将样本变换为平均为0，方差为1的标准化数据；即对每一个风速分量作标准化变换，变换公式为：

xij=vij-vj‾Sj(i=1,2,L,n;j=1,2,L,m)]]>

其中：——样本均值，vj‾=1nΣk=1nvkj;]]>

S_j——样本标准差，Sj=1n-1Σk=1n(vkj-vj‾);]]>

x_ij——标准化后的第i时刻第j台风机的风速数据；标准化后的风速矩阵可表示为

X=x11x12···x1nx21x22···x2nMxm1xm2···xmn]]>

(2)计算风电场所有机组风速间的相关矩阵

对于风电场内的n台风机，所有风机间的相关系数所构成的矩阵就是相关矩阵R，相关矩阵R中的每一个元素由相应的相关系数所表示；

R=XXT=1r12···r1nr211···r2nMrn1rn2···1]]>

其中：r_ij——相关系数，rij=1n-1Σk=1nxkixkj;]]>

(3)找出相关矩阵的特征值和特征向量

由相关矩阵R，求解特征方程：

|R-λI|＝0

其中：I——单位矩阵，为n阶方阵；

λ——R的特征值；

通过求解特征方程，可得到n个特征值λ_i(i＝1，2，L，n)，和对应于每一个特征值的特征向量ξ_i＝(ξ_i1，ξ_i2，L，ξ_in)(i＝1，2，L，n)；且有λ1≥λ2≥λ3≥L≥λn≥0与之对应的特征向量相互正交；

再将求解出来的特征值λ_i(i＝1，2，L，n)代入齐次代数方程组

|λ_iI-R|X＝0

展开为

λi-1-r12···-r1n-r21λi-1···-r2nMM···M-rn1-rn2···λi-1ξ1iξ2iMξni=0]]>

从而可以求出相应的特征向量，且ξ12+ξ22+L+ξn2=1;]]>

(4)求取等效风速

根据得到的n个特征向量，可以把n台风机感受风速的主要成分表示为：

F1=ξ11x1+ξ12x2+L+ξ1nxnF2=ξ21x1+ξ22x2+L+ξ2nxnLFn=ξn1x1+ξn2x2+L+ξnnxn]]>

可写为通式

F_i＝ξ_i1x₁+ξ_i2x₂+L+ξ_inx_n(i＝1，2，L，n)

以上求得的主成分相互正交，且每一个主成分的载荷系数之平方和等于对应的特征根λ；由于λ₁≥λ₂≥λ₃≥L≥λ_n≥0，各主要成分对应的方差是逐次递减的；若用p_i表示第i个主成分的方差解释率，即

pi=λiΣi=1nλi(i=1,2,L,n)]]>

那么，前q台风机的累积方差解释率可以表示为

Σj=1qpi=Σj=1qλjΣi=1nλi(i=1,2,L,n;j=1,2,L,q)]]>

选取特征值大于1且累积贡献率超过总方差80％的因子个数；假设满足要求的因子个数为t，则整个风电场的等效风速可表示为：

v_i(eq)＝ξ_i1x₁+ξ_i2x₂+L+ξ_inx_n(i＝1，2，L，t)

3).求取风电场发电功率：

在求出风电场等效风速的基础上，本发明采用支持向量机的方法求出等效风速对应的风电场发电功率；

支持向量机回归问题是通过非线性映射Φ：(m₀≥n₀，m₀、n₀分别表示空间维数)，把输入空间的样本X映射到一个高维特征空间，然后在该空间中做线性回归；对于给定训练数据集其中输入数据X_k∈R^N，输出数据Y_k∈R，支持向量机对应的函数回归估计为

Y(X)＝ωΦ(X)+b

式中：ω、Φ(X)为m₀维向量；b为偏置量，并且ω和b可以通过下式来确定：

minω,b,ξ,ξ*12ωT-ω+CΣi=1N(ξi+ξi*)]]>

Y_i-[ω^T-Φ(X_i)]+b≤

s.t.ϵ+ξi[ωT-Φ(Xi)]+b-Yi≤ϵ+ξi*]]>

(ξi,ξi*≥0,i=1,2,···,N)]]>

利用Langrange函数和Wolfe的对偶理论，并利用核技巧在高维空间求解上式中的ω，其中核函数的选取有多项式函数K(X_kX_l)＝[X_kX_l-c₀]^d，c₀≥0，c₀为一个参数；径向基函数K(X_k，X_l)＝exp(-||X_k-X_l||/σ²)；Sigoid函数K(X_k，X_l)＝tanh[k(X_k-X_l)+v]，k＞0，v＜0等；最终得到ω表达式为

ω=Σi=1N(αi-αi*)Φ(Xi)]]>

根据Karush-Kuhn-Tucker(KKT)条件可得到系数b，相应回归函数为：

Y(X)=Σi=1N(αi-αi*)K(Xk,X)+b]]>

式中不为零的α_i，对应的向量称为支持向量；得到支持向量后，即可求得回归函数Y(X)；

把上面求出的等效风速数据进行分类，第1类数据用于构建风电场等值模型，第2类数据用于模型的优化，第3类数据用于模型的验证；将第1类数据和相应的已有发电功率经过支持向量回归机的计算，可以得到一条输出曲线

P_发电功率＝f(v_1(eq)，v_2(eq)，L，v_t(eq))

其中：p_发电功率——风电场的发电功率；

v_i(eq)——表示第i个主要成分的等效风速，(i＝1，2，L，t)；

通过第2类数据采用最小二乘法对模型进行优化，最后再通过第3类数据对模型进行验证；