复杂山地风电场后评估技术及方法探析.docxVIP

【摘要】近十年来国内装机容量增长迅猛，部分风电场实际发电利用小时数明显低于设计水平或厂家保证的水平。特别针对国内的复杂山地风电场，发电量不达设计值有两大主因，第一是设计阶段的风流场建模错误，第二是运行阶段的风机性能劣化。本文主要结合SCADA数据分析与风电场建模方法，准确判断发电量低的根本原因，创新性地实现功率曲线量化对标，并探索风电机组性能监督方法和发电量提升手段。

【关键词】?流场整定，功率曲线量化对标，SCADA，数据交叉对比，理论发电小时数提升

1流场整定

首先，我们需要明确哪些因素可能导致复杂山地风电场设计阶段的风流数据模拟产生错误：1）测风时间不足，2）测风塔选址不具备代表性，3）风流场建模模型选择不合理，4）微观选址点位变动后未能重新进行流场建模。因此，在进行后评估工作之前，需要先进行流场整定工作，其目的在于获得风电场真实的风流特性图谱，进而得到每台机组真实的理论发电量数值，从而为SCADA功率曲线导出的理论发电量对标奠定基础。为了更好地考虑复杂山地对风流场的影响，我们采用借助计算流体力学（CFD）技术计算流体的特性。描述流体运动的方程是质量和动量守恒方程（Navier-Stokes方程）。三维瞬时Navier-Stokes方程的一般形式（笛卡尔张量形式）由方程（1）给出：??????

???

（1）连续性方程为方程（2）：?

（2）瞬时风速是一个随机变量，通常将瞬时风速分解为平均风速（独立于时间变化）与随机分量（随时间变化很快）。因此瞬时风速可以写成如下方程（3）：

（3）

把方程（3）代入方程（1）与方程（2），并对时间进行积分，就可以得到如下方程（4）：

（4）本文采用计算流体力学（CFD）软件对相关项目进行评估及分析。它求解雷诺平均方程与质量平均方程。对于不可压缩稳态流体，方程（2）将变为方程（5），方程（4）将变为方程（6）如下所示：

（5）

（6）

湍流部分则借助湍流粘性的概念进行参数化。此方法建立在雷诺应力项与平均风速梯度成正比的假设之上，如方程（7）所示：

（7）湍流粘性等于湍流长度尺度与速度尺度的乘积，而后两者都是湍流的特征尺度。湍流的速度尺度由湍流动能的平方根给出，如方程（8）所示：（8）

湍流动能则是由传输方程（9）给出，如下：

（9）

湍流动能产生项和耗散项分别表达如方程（10）与方程（11）：（10）（11）边界条件是由软件自动生成，例如入口边界条件在垂直方向上被划分为三个不同区域以充分考虑大气边界层的结构。第一层从地表到高度hs，为表面层，风廓线特征遵循对数法则。过渡层，高度从hs到hbl，为Eckman层。在更高的区域中，地理风条件被应用。湍流动能是一个常量，根据入口处定义的粗糙度进行计算。计算区域的底层的边界条件为具有粗糙度赋值的墙壁类型，粗糙度的影响是由处在障碍物范围内的网格中的吸收井项来模拟的。对称条件被应用于计算区域的侧面，其含义是没有风流和标量通量穿过边界。顶面边界条件为Neumann条件：所有变量在Z方向上的导数均为零。出口边界条件为压力条件。

2?流场整定过程中的网格生成机制

网格自动生成技术采用结构化网格对整个空间进行离散化，并对特定区域（即绘图区域）处进行网格加密，网格单元维度的扩展以及长宽比受到控制，以避免收敛的不稳定性。特别值得注意的是，本文采用的CFD计算方法会根据不同的来风方向重新生成网格，可以使网格与来风方向正交化本文采用的网格分辨率为水平方向25米，垂直方向5米。在数值求解方面，MIGAL-S求解器可以在整个计算区域上在每个迭代循环上同时求解风速分量和压强，也就是通常所说的耦合求解器。这种方法要求更强大的存储容量，但具有显著加速收敛的优点。全部的离散方程由一个矩阵系统来表示，作为ILU(0)类型的不完全LU分解的对象。为了改进求解的稳健性，进行了一个预处理来增大主对角线上的项。一个GMRES类型的预解算子也被用于该过程。此外，MIGAL-S求解器使用了多重网格求解技术。这种方式实现在不同疏密的网格中陆续求解方程（从最细致的到最粗的网格），这样使得收敛最慢的误差的低频分量收敛更快。针对所求解问题的类型，测试表明V型循环比W型循环在遍历网格时能够更快地收敛。而疏密程度不同的网格的层数则自动由MIGAL-S求解器进行计算。

图1网格生成效果图（风向=正北360度）

3流场整定结果

通过对比主导风向的整定后的流场图谱可以清晰地看到与可研阶段在风速和发电量上的差异。

图2流场整定前后

风速图谱的差异利用整定后的图谱，结合专业CFD软件，可以重新计算每台机组的理论发电小时数。就一个国内20MW风电场而言，其各台机组实际上网电量均远远低于可研值，本文对其进行了流场整定，整定前后的结果如表1所示：

非常明