[发明专利]一种故障风机损失发电量的估计方法

摘要

本发明公开了一种故障风机损失发电量的估计方法，包括如下步骤1)根据风机的地理信息，建立风机的无向图模型；2)风机发生故障停运时根据风机的无向图模型，试图获取地理位置距离故障风机最近的两台正常运行风机；3)如果获取成功，利用地理位置距离故障风机最近的两台正常运行风机风速风向，对故障风机故障期间所损失的发电量进行估计；4)如果获取失败，或者所获取的风机与故障风机之间的距离大于测风塔与故障风机之间的距离，则利用测风塔风速风向，对故障风机故障期间所损失的发电量进行估计。本发明的估计方法，提供了对风机故障期间所损失的发电量的合理估计，可作为风电场对风机故障造成的经济损失进行量化评估时的数据依据。

主权项

一种故障风机损失发电量的估计方法，其特征在于，包括以下步骤：1)根据风机的地理信息，建立风机的无向图模型；2)风机发生故障停运时根据风机的无向图模型，获取地理位置距离故障风机最近的两台正常运行风机；3)如果获取成功，利用地理位置距离故障风机最近的两台正常运行风机风速风向，对故障风机故障期间所损失的发电量进行估计；4)如果获取失败，或者所获取的风机与故障风机之间的距离大于测风塔与故障风机之间的距离，则利用测风塔风速风向，对故障风机故障期间所损失的发电量进行估计；在所述步骤3)中，风机发生故障停运时根据风机的无向图模型，获取地理位置距离故障风机最近的两台正常运行风机，获取方法如下：设风电场共有n台风机，D(ni,nj)表示两台风机ni与nj之间的物理距离，D(ni,T)表示风机与测风塔T之间的物理距离，其中1≤i≤n，1≤j≤n；增设两台虚拟风机nn+1与nn+2，任意一台风机与该两台虚拟风机之间的物理距离为无穷大，即D(ni,nn+1)＝∞，D(ni,nn+2)＝∞，其中1≤i≤n，∞表示无穷大；假设风机A在时间段t∈[t1,t2]发生故障停运，首先将距离风机A最近的两台正常运行的风机B和风机C设为增设的虚拟风机：B＝nn+1，C＝nn+2；遍历n台风机，进行n轮循环，其中对风机ni(1≤i≤n)进行如下操作：21)如果ni为风机A本身，则结束本轮循环并开始下一轮循环，否则进行步骤22)；22)如果风机ni在时间段t∈[t1,t2]内正常运行，并且D(ni,A)＜D(B,A)，则置B为B＝ni，结束本轮循环并开始下一轮循环，否则进行步骤23)；23)如果风机ni在时间段t∈[t1,t2]内正常运行，并且D(ni,A)＜D(C,A)，则置C为C＝ni，结束本轮循环并开始下一轮循环；循环结束后，如果B≠nn+1，D(A,B)＜D(A,T)且C≠nn+2，D(A,C)＜D(A,T),则表明获取成功，否则获取失败；在所述步骤3)中，如果获取成功，利用地理位置距离故障风机最近的两台正常运行风机风速，对故障风机故障期间所损失的发电量进行估计，估计方法如下：设风机B和风机C的风速分别为v1(t)和v2(t)，风向为d1(t)和d2(t)；假设在时间段[t1,t2]内d1(t)和d2(t)与南北方向的夹角分别为θ1(t)和θ2(t)，其中θ1,θ2∈[0°,90°]，则风机B和风机C的风速在南北方向的分量vHB(t)、vHC(t)和东西方向的分量vLB(t)、vLC(t)分别为：风机A和风机B、C之间的距离D(A,B)、D(A,C)在南北方向的分量DH和东西方向的分量DL分别为：风机A的风速在南北方向的估计分量和东西方向的估计分量为风机B和C在两个方向上的分量根据距离成比例的叠加：v^HA(t)=DH(A,B)vHB(t)+DH(A,C)vHC(t)DH(A,B)+DH(A,C)v^LA(t)=DL(A,B)vLB(t)+DL(A,C)vLC(t)DL(A,B)+DL(A,C),t∈[t1,t2]---(5)]]>将式(1)、(2)、(3)、(4)代入式(5)可得:v^HA(t)=D(A,B)cos(θ1(t))v1(t)cos(θ1(t))+D(A,C)cos(θ2(t))v2(t)cos(θ2(t))D(A,B)cos(θ1(t))+D(A,C)cos(θ2(t))v^LA(t)=D(A,B)sin(θ1(t))v1(t)sin(θ1(t))+D(A,C)sin(θ2(t))v2(t)sin(θ2(t))D(A,B)sin(θ1(t))+D(A,C)sin(θ2(t))---(6)]]>得到风机A的风速在南北方向和东西方向两个估计分量后，可得风机A的风速估计值如下:v^(t)=(v^HA(t))2+(v^LA(t))2,t∈[t1,t2]---(7)]]>对风机A的实际输出功率特性曲线分段使用最小二乘法进行曲线拟合，获得风机A的输出功率特性曲线表达式P(v)；风机A在时间段t∈[t1,t2]内的发电量估计值如下为：K^=∫t1t2P(v^(t))dt=∫t1t2P((v^HA(t))2+(v^LA(t))2dt---(8)]]>将式(6)代入式(8)，并对此复合函数积分求解，便可得风机A在时间段t∈[t1,t2]内的发电量估计值。