CN105048444A - 一种基于测风塔测风数据的风电场弃风电量确定方法 - Google Patents

Abstract

本发明涉及一种基于测风塔测风数据的风电场弃风电量确定方法，该方法以测风塔测风数据为基础，考虑地形、地貌及风电机组尾流等影响，计算风电机组轮毂高度处的风速，再根据功率曲线获得风电机组理论功率，对风电场所有风电机组的理论功率求和得到风电场理论功率，最后将风电场理论功率和实际功率分别对时间积分，得到理论电量和实际电量，两者相减得到风电场弃风电量，实现对风电场弃风电量的精确计算。该方法普遍适用于各种风电场，成功解决了样板机法计算弃风电量精确度不高的问题。

Description

一种基于测风塔测风数据的风电场弃风电量确定方法

技术领域

本发明涉及一种风力发配电技术，具体讲涉及一种基于测风塔测风数据的风电场弃风电量确定方法。

背景技术

大型风电基地覆盖面积广，一般包含多个风电场或风电场群，由于电网输送通道受限、调峰能力不足等因素会导致一定程度的弃风。正确科学的认识弃风、以合理的方式计算弃风电量将不仅有助于风电产业的健康、平稳发展，还能提升电网调度运行水平、推动风电规划与电网规划的协调发展、提高清洁能源利用率。

目前，风电行业对弃风电量的评估主要采用的是样板机方法，即在弃风时段保留不超过风机总台数10％的样板机正常运行，根据样板机的实际功率计算风电场的理论功率。然而，受地形、地貌等因素的影响，风电场每台风机的出力并非简单的线性关系，样板风机方法计算精度受到一定程度的制约。因此，需要提供一种以测风塔测风数据为基础，综合考虑地形、地貌、尾流等影响，计算风电场的理论电量，并与实际电量相比的方法，以实现对弃风电量的精确计算。

发明内容

针对现有技术的不足，本发明的目的是提供一种基于测风塔测风数据的风电场弃风电量确定方法，该方法以测风塔测风数据为基础，考虑地形、地貌及风电机组尾流等影响，计算风电机组轮毂高度处的风速，再根据功率曲线获得风电机组理论功率，对风电场所有风电机组的理论功率求和得到风电场理论功率，最后将风电场理论功率和实际功率分别对时间积分，得到理论电量和实际电量，两者相减得到风电场弃风电量，实现对风电场弃风电量的精确计算。

本发明的目的是采用下述技术方案实现的：

本发明提供一种基于测风塔测风数据的风电场弃风电量确定方法，其改进之处在于，所述方法包括下述步骤：

(1)确定测风数据在地形、地貌和风机尾流影响下的变化；

(2)确定风电机组的理论功率；

(3)确定风电场理论功率；

(4)从能量管理系统中获得风电场实际功率；

(5)确定风电场弃风电量。

进一步地，所述步骤(1)包括：

①确定地形对风速的影响：地形扰动下的受扰边界层包括内层、外层和中间层；

②确定地貌对风速的影响；

③确定风电机组尾流对风速的影响。

进一步地，所述①中，风电机组边界层外层风速的扰动的求解包括：

边界层外层在地形扰动下的流场变化按照势流理论求解，将地形变化看作对未受扰流场的小扰动,由势流理论得：

u'＝▽χ1)；

其中：u'为变化地形对未受扰流场的扰动；χ(r,φ,z)为柱坐标下的势函数，▽为哈密顿算子；r、φ、z分别表示柱坐标系中的三个坐标变量；

以风电机组所在位置为坐标原点，则边界层外层在地形扰动下的流场变化转换为求解以下定解问题：

其中：u₀为上风向未受扰水平风速矢量；h(r,φ)为地形高度函数；R为研究区域半径，R＝10km为10km外的地形变化将不再影响风电机组位置处的流场；L为地形扰动在垂直方向上的长度；

χ(r,φ,z)的通解：

$\mathrm{\χ}(r,\mathrm{\φ},z)=\underset{n=0}{\overset{\∞}{\mathrm{\Σ}}}\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}{J}_n\left({\mathrm{\α}}_{j}r\right)(A_{\mathrm{nj}}\cos \mathrm{n\φ}+B_{\mathrm{nj}}\sin \mathrm{n\φ})\sinh {\mathrm{\α}}_j(L_j-z)---3);$

式中,J_n(α_jr)为n阶贝塞尔函数；

由边界条件和贝塞尔函数的正交性以及表达式1)得：

${\stackrel{\→}{u}}_{\mathrm{WT}}^{\′}=\frac{1}{2R}\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}\sinh \left[\frac{c_j^1}{R}(L_j-z){]c}_j^1\right[(A_{1j}+B_{1j})e_r(B_{1j}-A_{1j})e_{\mathrm{\φ}}]---4);$

其中：为地形变化对风电机组位置处边界层外层流场的扰动，为一阶贝塞尔函数的第j个零点；e_r、e_φ分别为径向与方位角方向的单位向量；系数A_1j、B_1j由下式5)确定：

$A_{1j}+{\mathrm{iB}}_{1j}=-\frac{2{\∫}_0^R{\∫}_0^{2\mathrm{\π}}{\stackrel{\→}{u}}_0\·\▿h(r,\mathrm{\φ})e^{\mathrm{i\φ}}r{J}_1\left(\frac{c_j^1}{R}r\right)\mathrm{d\φdr}}{\mathrm{\π}R\cosh \left(1\right)c_j^1{\left[J_2\left(c_j^1\right)\right]}^2}---5);$

其中：▽h(r,φ)包含地形变化信息；

根据表达式4)和5)，在给定风电机组地理位置以及地形变化信息后，获得地形变化对风电机组位置边界层外层流场的扰动。

进一步地，所述①中，对风电机组边界层内层风速的扰动求解，包括：

边界层内层流场扰动随高度按对数风廓线变化，在边界层内层层顶扰动达到最大并大于势流解，内层流场对于同一高度势流解的修正值为：

$\mathrm{\Δ}{\stackrel{\→}{u}}_j\left(z\right)=U_0\left(z\right)\frac{U_0^2\left(L_j\right)}{U_0^2\left(z_j^{\′}\right)}\▿{\mathrm{\χ}}_j---6);$

式中，U₀(z)为上风向未受扰风矢量在高度z的风速，L_j为地形扰动在垂直方向上的长度尺度，z'_j＝max(z,l_j)，其中l_j为边界层内层高度，l_j＜＜L_j，l_j由下式确定：z_oj为对应l_j的相对粗糙度，上风向为均一地形时，z_0j＝z₀，上风向为非均一地形时： $z_{0j}=z_{0n}\×{\left(\frac{z_{0n+1}}{z_{0n}}\right)}^{w_n},$ 其中： $W_n=\exp (-\frac{x_n}{D}),D={5L}_j,$ x_n为第n个粗糙度变化与风电机组的距离；▽为哈密顿算子。

进一步地，所述①包括：

边界层中间层的范围为l_j≤z≤4l_j，边界层中间层流场扰动为：

${\stackrel{\→}{u}}_{m_j}\left(z\right)=k_{\mathrm{wf}}[\mathrm{\Δ}{\stackrel{\→}{u}}_{l_j}\left(z\right)+\▿{\mathrm{\χ}}_{l_j}]+(1-k_{\mathrm{wf}})\▿{\mathrm{\χ}}_{{4l}_j}---7);$

其中：k_wf为加权因子，Δu_j(l_j)为边界层内层顶势流解修正值；分别为z＝l_j与z＝4l_j的势流解；l_j为边界层内层高度，l_j＜＜L_j，l_j由下式确定：z_oj为对应l_j的相对粗糙度，上风向为均一地形时，z_0j＝z₀，上风向为非均一地形时： $z_{0j}=z_{0n}\×{\left(\frac{z_{0n+1}}{z_{0n}}\right)}^{w_n},$ 其中： $W_n=\exp (-\frac{x_n}{D}),D={5L}_j,$ x_n为第n个粗糙度变化与风电机组的距离。

进一步地，所述②包括：

流经变化粗糙度的下风向风廓线描述为：

$u\left(z\right)=\left\{\begin{array}{cc}{u}^{\′}\frac{\ln (z/z_{01})}{\ln (0.3h/z_{01})}& z\≥0.3h\\ u^{\′\′}+(u^{\′}-u^{\′\′})\frac{\ln (z/0.09h)}{\ln (0.3/0.09)}& 0.09h\≤z\≤0.3h\\ u^{\′\′}\frac{\ln (z/z_{02})}{\ln (0.09h/z_{02})}& z\≤0.09h\end{array}\right.---8);$

其中：z₀₂为风电机组位置的粗糙度，z₀₁为距离风电机组位置最近的上风向粗糙度，u_*2、u_*1分别为对应z₀₁、z₀₂的摩擦速度，κ＝0.4为卡曼常数，u_*1为z₀₁对应的摩擦速度，h为边界层内层高度，由下式确定：

$\frac{h}{z_0^{\′}}(\ln \frac{h}{z_0^{;}}-1)=0.9\frac{x}{z_0^{\′}}---9);$

其中：z'₀＝max(z₀₁,z₀₂)，x为粗糙度变化位置与风电机组位置的距离；表示研究区域内粗糙度的最大值；

粗糙度变化扰动下，摩擦速度之间的关系，如下式：

$\frac{u_{*n+1}}{u_{*n}}=\frac{\ln (h/z_{0n})}{\ln (h/z_{0n+1})}---10);$

式中，z_0n、z_0n+1分别为上风向粗糙度与距离最近的下风向粗糙度，u_*n、u_*n+1为对应z_0n、z_0n+1的摩擦速度；

粗糙度变化位置距离风电机组位置越远其影响亦越弱，加入距离权重因子表示距离的作用，得：

$z_{0\mathrm{neffe}}=z_{0n+1}\×{\left(\frac{Z_{0n}}{z_{0n+1}}\right)}^{w_n}---11);$

式中：z_0effe为等效粗糙度，为第n个粗糙度的距离权重因子，D＝10km，即认为10km外的粗糙度变化将不再对风电机组位置的风廓线产生影响。

进一步地，所述③中，确定风电机组尾流对风速的影响包括：确定风电机组尾流模型；

尾流模型被称为Larsen尾流模型，假定下风向不同位置的风速衰减具有相似性，并且风速只发生中等程度的衰减，则通过下式计算下风向L＝x处的尾流影响区域半径：

$\left\{\begin{array}{c}{R}_w={\left[\frac{35}{2p}\right]}^{\frac{1}{5}}{\left[3c_1^2\right]}^{\frac{1}{5}}{\left[C_T\mathrm{Ax}\right]}^{\frac{1}{3}}\\ c_1=l{\left(C_T\mathrm{Ax}\right)}^{-\frac{1}{3}}\end{array}\right.---12);$

其中：c₁为无量纲混合长；l为普朗特混合长，A为风力机扫风面积，C_T为风电机组推力系数；

Larsen尾流模型最终的风速衰减表达式为：

$\mathrm{DU}=-\frac{U_{\mathrm{WT}}}{9}{\left(C_T{\mathrm{Ax}}^{-2}\right)}^{\frac{1}{3}}[R_w^{\frac{3}{2}}{\left(3c_1^2{C}_T\mathrm{Ax}\right)}^{-\frac{1}{2}}-\left(\frac{{35}^{\frac{3}{10}}}{2\mathrm{\π}}{\left(3c_1^2\right)}^{-\frac{1}{5}}\right){]}^2---13);$

式中：U_WT为风电机组轮毂高度的平均风速；R_w由式12)确定；

进一步地，所述步骤(2)中，根据功率曲线，得到风电机组理论功率；功率曲线由风机制造商提供；

所述步骤(3)中，对风电场所有风电机组理论功率求和得到风电场理论功率。

进一步地，所述步骤(5)中，将风电场理论功率和风电场实际功率分别对时间积分，得到风电场理论电量和风电场实际电量，两者相减得到风电场弃风电量；

进一步地，利用分段二次插值法求出设定时段内各风电机组理论功率，对时间积分得到风电场理论电量：其中P_T为风电场理论功率，t₀、t₁分别为起始时间和结束时间；

风电场理论电量和风电场实际电量相减得到风电场弃风电量，表达式如下：E_c＝E_T-E_M；其中：E_c为风电场弃风电量；E_T为风电场理论电量；E_M为风电场实际电量。

与现有技术比，本发明达到的有益效果是：

本发明提供的风电场弃风电量确定方法基于测风塔测风数据，采用了基于势流理论的粗糙度变化模型与地形变化模型，可以有效反映风电场局地效应对大气边界层风的扰动作用。解析求解大气运动方程有利于缩短计算时间，满足弃风电量计算的时效要求，且降低了对计算机计算能力的要求。模型涉及参数较少，鲁棒性与工程实用性强。该方法普遍适用于各种风电场，成功解决了样板机法计算弃风电量精确度不高的问题。

附图说明

图1是本发明提供的基于测风塔测风数据的风电场弃风电量确定方法的流程图；

图2是本发明提供的粗糙度变化下的边界层内层发展示意图；

图3是本发明提供的典型功率曲线示意图。

具体实施方式

下面结合附图对本发明的具体实施方式作进一步的详细说明。

本发明提供的基于测风塔测风数据的风电场弃风电量确定方法的流程图如图1所示，包括下述步骤：

(1)确定测风数据在地形、地貌和风机尾流影响下的变化；所述测风数据包括风速，确定测风数据在地形、地貌和风机尾流影响下的变化包括：

①确定地形对风速的影响：地形扰动下的受扰边界层包括内层、外层和中间层；

一、对风电机组边界层外层风速的扰动求解，包括：

边界层外层在地形扰动下的流场变化按照势流理论求解，将地形变化看作对未受扰流场的小扰动,由势流理论得：

u'＝▽χ1)；

其中：u'为变化地形对未受扰流场的扰动；χ(r,φ,z)为柱坐标下的势函数，▽为哈密顿算子；r、φ、z分别表示柱坐标系中的三个坐标变量；

以风电机组所在位置为坐标原点，则边界层外层在地形扰动下的流场变化转换为求解以下定解问题：

其中：u₀为上风向未受扰水平风速矢量；h(r,φ)为地形高度函数；R为研究区域半径，R＝10km,即认为10km外的地形变化将不再影响风电机组位置处的流场；L为地形扰动在垂直方向上的长度；

χ(r,φ,z)的通解：

$\mathrm{\χ}(r,\mathrm{\φ},z)=\underset{n=0}{\overset{\∞}{\mathrm{\Σ}}}\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}{J}_n\left({\mathrm{\α}}_{j}r\right)(A_{\mathrm{nj}}\cos \mathrm{n\φ}+B_{\mathrm{nj}}\sin \mathrm{n\φ})\sinh {\mathrm{\α}}_j(L_j-z)---3);$

式中,J_n(α_jr)为n阶贝塞尔函数；

由边界条件和贝塞尔函数的正交性以及表达式1)得：

${\stackrel{\→}{u}}_{\mathrm{WT}}^{\′}=\frac{1}{2R}\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}\sinh \left[\frac{c_j^1}{R}(L_j-z){]c}_j^1\right[(A_{1j}+B_{1j})e_r(B_{1j}-A_{1j})e_{\mathrm{\φ}}]---4);$

其中：为地形变化对风电机组位置处边界层外层流场的扰动，为一阶贝塞尔函数的第j个零点；e_r、e_φ分别为径向与方位角方向的单位向量；系数A_1j、B_1j由下式5)确定：

$A_{1j}+{\mathrm{iB}}_{1j}=-\frac{2{\∫}_0^R{\∫}_0^{2\mathrm{\π}}{\stackrel{\→}{u}}_0\·\▿h(r,\mathrm{\φ})e^{\mathrm{i\φ}}r{J}_1\left(\frac{c_j^1}{R}r\right)\mathrm{d\φdr}}{\mathrm{\π}R\cosh \left(1\right)c_j^1{\left[J_2\left(c_j^1\right)\right]}^2}---5);$

其中：▽h(r,φ)包含地形变化信息；

根据表达式4)和5)，在给定风电机组地理位置以及地形变化信息后，获得地形变化对风电机组位置边界层外层流场的扰动。

二、对风电机组边界层内层风速的扰动求解，包括：

边界层内层流场扰动随高度按对数风廓线变化，在边界层内层层顶扰动达到最大并大于势流解，内层流场对于同一高度势流解的修正值为：

$\mathrm{\Δ}{\stackrel{\→}{u}}_j\left(z\right)=U_0\left(z\right)\frac{U_0^2\left(L_j\right)}{U_0^2\left(z_j^{\′}\right)}\▿{\mathrm{\χ}}_j---6);$

式中，U₀(z)为上风向未受扰风矢量在高度z的风速，L_j为地形扰动在垂直方向上的长度尺度，z'_j＝max(z,l_j)，其中l_j为边界层内层高度，l_j＜＜L_j，l_j由下式确定：z_oj为对应l_j的相对粗糙度，上风向为均一地形时，z_0j＝z₀，上风向为非均一地形时： $z_{0j}=z_{0n}\×{\left(\frac{z_{0n+1}}{z_{0n}}\right)}^{w_n},$ 其中： $W_n=\exp (-\frac{x_n}{D}),D={5L}_j,$ x_n为第n个粗糙度变化与风电机组的距离。

三、对风电机组边界层中间层风速的扰动求解，包括：

边界层中间层的范围为l_j≤z≤4l_j，边界层中间层流场扰动为：

${\stackrel{\→}{u}}_{m_j}\left(z\right)=k_{\mathrm{wf}}[\mathrm{\Δ}{\stackrel{\→}{u}}_{l_j}\left(z\right)+\▿{\mathrm{\χ}}_{l_j}]+(1-k_{\mathrm{wf}})\▿{\mathrm{\χ}}_{{4l}_j}---7);$

其中：k_wf为加权因子，Δu_j(l_j)为边界层内层层顶势流解修正值；； 分别为z＝l_j与z＝4l_j的势流解；l_j为边界层内层高度，l_j＜＜L_j，l_j由下式确定：z_oj为对应l_j的相对粗糙度，上风向为均一地形时，z_0j＝z₀，上风向为非均一地形时： $z_{0j}=z_{0n}\×{\left(\frac{z_{0n+1}}{z_{0n}}\right)}^{w_n},$ 其中： $W_n=\exp (-\frac{x_n}{D}),D={5L}_j,$ x_n为第n个粗糙度变化与风电机组的距离。

②确定地貌对风速的影响，包括：

流经变化粗糙度的下风向风廓线描述为：

$u\left(z\right)=\left\{\begin{array}{cc}{u}^{\′}\frac{\ln (z/z_{01})}{\ln (0.3h/z_{01})}& z\≥0.3h\\ u^{\′\′}+(u^{\′}-u^{\′\′})\frac{\ln (z/0.09h)}{\ln (0.3/0.09)}& 0.09h\≤z\≤0.3h\\ u^{\′\′}\frac{\ln (z/z_{02})}{\ln (0.09h/z_{02})}& z\≤0.09h\end{array}\right.---8);$

其中：z₀₂为风电机组位置的粗糙度，z₀₁为距离风电机组位置最近的(风电机组位置附近一般有多个粗糙度区域，例如图2中，风机附近有一片树林和一块粗糙度区域，这里那块区域的粗糙度就是距离风电机组最近的粗糙度)上风向粗糙度， u_*2、u_*1分别为对应z₀₁、z₀₂的摩擦速度，κ＝0.4为卡曼常数，u_*1为z₀₁对应的摩擦速度，h为边界层内层高度，由下式确定：

$\frac{h}{z_0^{\′}}(\ln \frac{h}{z_0^{;}}-1)=0.9\frac{x}{z_0^{\′}}---9);$

其中：z'₀＝max(z₀₁,z₀₂)，x为粗糙度变化位置与风电机组位置的距离；表示研究区域内粗糙度的最大值；

粗糙度变化扰动下，摩擦速度之间的关系，如下式：

$\frac{u_{*n+1}}{u_{*n}}=\frac{\ln (h/z_{0n})}{\ln (h/z_{0n+1})}---10);$

式中，z_0n、z_0n+1分别为上风向粗糙度与距离最近的下风向粗糙度，u_*n、u_*n+1为对应z_0n、z_0n+1的摩擦速度；

粗糙度变化位置距离风电机组位置越远其影响亦越弱，加入距离权重因子表示距离的作用，得：

$z_{0\mathrm{neffe}}=z_{0n+1}\×{\left(\frac{Z_{0n}}{z_{0n+1}}\right)}^{w_n}---11);$

式中：z_0effe为等效粗糙度，为第n个粗糙度的距离权重因子，D＝10km，即认为10km外的粗糙度变化将不再对风电机组位置的风廓线产生影响。

③确定风电机组尾流对风速的影响，包括：确定风电机组尾流模型；

尾流模型被称为Larsen尾流模型，假定下风向不同位置的风速衰减具有相似性，并且风速只发生中等程度的衰减，则通过下式计算下风向L＝x处的尾流影响区域半径：

$\left\{\begin{array}{c}{R}_w={\left[\frac{35}{2p}\right]}^{\frac{1}{5}}{\left[3c_1^2\right]}^{\frac{1}{5}}{\left[C_T\mathrm{Ax}\right]}^{\frac{1}{3}}\\ c_1=l{\left(C_T\mathrm{Ax}\right)}^{-\frac{1}{3}}\end{array}\right.---12);$

其中：c₁为无量纲混合长；l为普朗特混合长，A为风力机扫风面积，C_T为风电机组推力系数；

Larsen尾流模型最终的风速衰减表达式为：

$\mathrm{DU}=-\frac{U_{\mathrm{WT}}}{9}{\left(C_T{\mathrm{Ax}}^{-2}\right)}^{\frac{1}{3}}[R_w^{\frac{3}{2}}{\left(3c_1^2{C}_T\mathrm{Ax}\right)}^{-\frac{1}{2}}-\left(\frac{{35}^{\frac{3}{10}}}{2\mathrm{\π}}{\left(3c_1^2\right)}^{-\frac{1}{5}}\right){]}^2---13);$ 式中：U_WT为风电机组轮毂高度的平均风速；R_w由式12)确定。确定风电机组的理论功率：根据功率曲线，得到风电机组理论功率；典型功率曲线示意图如图3所示。

(3)确定风电场理论功率：对风电场所有风电机组理论功率求和得到风电场理论功率。

(4)从能量管理系统中获得风电场实际功率；

(5)确定风电场弃风电量：将风电场理论功率和风电场实际功率分别对时间积分，得到风电场理论电量和风电场实际电量，两者相减得到风电场弃风电量，具体为：

利用分段二次插值法求出设定时段内各风电机组理论功率，对时间积分得到风电场理论电量：其中P_T为风电场理论功率，t₀、t₁分别为起始时间和结束时间；

风电场理论电量和风电场实际电量相减得到风电场弃风电量，表达式如下：E_c＝E_T-E_M；其中：E_c为风电场弃风电量；E_T为风电场理论电量；E_M为风电场实际电量。

最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求范围当中。

Claims (10)

1.一种基于测风塔测风数据的风电场弃风电量确定方法，其特征在于，所述方法包括下述步骤：

(1)确定测风数据在地形、地貌和风机尾流影响下的变化；

(2)确定风电机组的理论功率；

(3)确定风电场理论功率；

(4)从能量管理系统中获得风电场实际功率；

(5)确定风电场弃风电量。

2.如权利要求1所述的风电场弃风电量确定方法，其特征在于，所述步骤(1)包括：

①确定地形对风速的影响：地形扰动下的受扰边界层包括内层、外层和中间层；

②确定地貌对风速的影响；

③确定风电机组尾流对风速的影响。

3.如权利要求2所述的风电场弃风电量确定方法，其特征在于，所述①中，风电机组边界层外层风速的扰动的求解包括：

边界层外层在地形扰动下的流场变化按照势流理论求解，将地形变化看作对未受扰流场的小扰动,由势流理论得：

u'＝▽χ1)；

其中：u'为变化地形对未受扰流场的扰动；χ(r,φ,z)为柱坐标下的势函数，▽为哈密顿算子；r、φ、z分别表示柱坐标系中的三个坐标变量；

以风电机组所在位置为坐标原点，则边界层外层在地形扰动下的流场变化转换为求解以下定解问题：

其中：u₀为上风向未受扰水平风速矢量；h(r,φ)为地形高度函数；R为研究区域半径，R＝10km为10km外的地形变化将不再影响风电机组位置处的流场；L为地形扰动在垂直方向上的长度；

χ(r,φ,z)的通解：

式中,J_n(α_jr)为n阶贝塞尔函数；

由边界条件和贝塞尔函数的正交性以及表达式1)得：

其中：为地形变化对风电机组位置处边界层外层流场的扰动，为一阶贝塞尔函数的第j个零点；e_r、e_φ分别为径向与方位角方向的单位向量；系数A_1j、B_1j由下式5)确定：

其中：▽h(r,φ)包含地形变化信息；

根据表达式4)和5)，在给定风电机组地理位置以及地形变化信息后，获得地形变化对风电机组位置边界层外层流场的扰动。

4.如权利要求2所述的风电场弃风电量确定方法，其特征在于，所述①中，对风电机组边界层内层风速的扰动求解，包括：

边界层内层流场扰动随高度按对数风廓线变化，在边界层内层层顶扰动达到最大并大于势流解，内层流场对于同一高度势流解的修正值为：

式中，U₀(z)为上风向未受扰风矢量在高度z的风速，L_j为地形扰动在垂直方向上的长度尺度，z'_j＝max(z,l_j)，其中l_j为边界层内层高度，l_j＜＜L_j，l_j由下式确定：z_oj为对应l_j的相对粗糙度，上风向为均一地形时，z_0j＝z₀，上风向为非均一地形时：其中：x_n为第n个粗糙度变化与风电机组的距离；▽为哈密顿算子。

5.如权利要求2所述的风电场弃风电量确定方法，其特征在于，所述①包括：

边界层中间层的范围为l_j≤z≤4l_j，边界层中间层流场扰动为：

其中：k_wf为加权因子，Δu_j(l_j)为边界层内层顶势流解修正值；分别为z＝l_j与z＝4l_j的势流解；l_j为边界层内层高度，l_j＜＜L_j，l_j由下式确定：z_oj为对应l_j的相对粗糙度，上风向为均一地形时，z_0j＝z₀，上风向为非均一地形时：其中：x_n为第n个粗糙度变化与风电机组的距离。

6.如权利要求2所述的风电场弃风电量确定方法，其特征在于，所述②包括：

流经变化粗糙度的下风向风廓线描述为：

其中：z₀₂为风电机组位置的粗糙度，z₀₁为距离风电机组位置最近的上风向粗糙度，u_*2、u_*1分别为对应z₀₁、z₀₂的摩擦速度，κ＝0.4为卡曼常数，u_*1为z₀₁对应的摩擦速度，h为边界层内层高度，由下式确定：

其中：z'₀＝max(z₀₁,z₀₂)，x为粗糙度变化位置与风电机组位置的距离；表示研究区域内粗糙度的最大值；

粗糙度变化扰动下，摩擦速度之间的关系，如下式：

式中，z_0n、z_0n+1分别为上风向粗糙度与距离最近的下风向粗糙度，u_*n、u_*n+1为对应z_0n、z_0n+1的摩擦速度；

粗糙度变化位置距离风电机组位置越远其影响亦越弱，加入距离权重因子表示距离的作用，得：

式中：z_0effe为等效粗糙度，为第n个粗糙度的距离权重因子，D＝10km，即认为10km外的粗糙度变化将不再对风电机组位置的风廓线产生影响。

7.如权利要求2所述的风电场弃风电量确定方法，其特征在于，所述③中，确定风电机组尾流对风速的影响包括：确定风电机组尾流模型；

尾流模型被称为Larsen尾流模型，假定下风向不同位置的风速衰减具有相似性，并且风速只发生中等程度的衰减，则通过下式计算下风向L＝x处的尾流影响区域半径：

其中：c₁为无量纲混合长；l为普朗特混合长，A为风力机扫风面积，C_T为风电机组推力系数；

Larsen尾流模型最终的风速衰减表达式为：

式中：U_WT为风电机组轮毂高度的平均风速；R_w由式12)确定。

8.如权利要求1所述的风电场弃风电量确定方法，其特征在于，所述步骤(2)中，根据功率曲线，得到风电机组理论功率；功率曲线由风机制造商提供；

所述步骤(3)中，对风电场所有风电机组理论功率求和得到风电场理论功率。

9.如权利要求1所述的风电场弃风电量确定方法，其特征在于，所述步骤(5)中，将风电场理论功率和风电场实际功率分别对时间积分，得到风电场理论电量和风电场实际电量，两者相减得到风电场弃风电量。

10.如权利要求9所述的风电场弃风电量确定方法，其特征在于，利用分段二次插值法求出设定时段内各风电机组理论功率，对时间积分得到风电场理论电量：其中P_T为风电场理论功率，t₀、t₁分别为起始时间和结束时间；

风电场理论电量和风电场实际电量相减得到风电场弃风电量，表达式如下：E_c＝E_T-E_M；其中：E_c为风电场弃风电量；E_T为风电场理论电量；E_M为风电场实际电量。