CN104794287B - 一种风工程尾流计算方法 - Google Patents

Abstract

本发明公开了一种风工程尾流计算方法，其特征在于，包括如下步骤：步骤一：计算尾流膨胀系数；步骤二：基于步骤一，计算Park模型下的尾流风速；步骤三：修正Park模型下的初始尾流半径；步骤四：基于步骤二和步骤三，计算修正后的Park模型尾流风速；步骤五：基于步骤四，建立新的尾流模型，将新的尾流模型命名为Park‑Gauss模型；步骤六：计算Park‑Gauss模型下的尾流风速。本发明修正后的初始尾流半径明显提高了尾流的预测精度；可以很好的模拟尾流区的风速，不仅在精度上与试验结果接近一致而且在径向分布上也更加符合真实流场，同时继承了工程模型的形式简单、容易编码、计算高效等优点。

Description

一种风工程尾流计算方法

技术领域

本发明涉及一种计算方法，具体涉及一种风工程尾流计算方法，本发明属于风力发电领域。

背景技术

随着不可再生能源短缺问题日益尖锐，作为未来可持续开发且绿色无污染的风能将成为新型能源领域的热点之一。风电场选址是风场建设前期的重要工作，其对后期风场建设和效益预测都有重要作用。在大型风场的微观选址中，受到目前计算资源和能力限制，采取全CFD方法模拟预测计算成本会很昂贵。在风能发电工程领域，一些精度可接受且形式简单的工程尾流模型具有很强的实用性。尾流效应会影响风场的功率产出和机械载荷，因此一个好的尾流模型可以更好地预测疲劳载荷和风场布局性能。Jensen和Katic等人提出了Park模型，这是个形式简单的经验模型，其尾流衰减常数k等于尾流的线性膨胀速率。Park模型假设尾流区风速在径向呈常数分布，这与真实情况不符合。经试验测量在真实的流场中，在径向尾流风速呈多项式分布或高斯型分布。1988年Ainslie提出了涡粘性模型(Eddy-viscosity model),基于轴对称坐标系下求解N-S方程，因此也被认为是一个简单的RANS模型，并假设尾流呈高斯型分布。Park模型和涡粘性模型，都低估了尾流区的风速亏损。2009年Brower等人基于考虑地表阻力引起的内部边界层方法提出了DAWM模型(Deep-Array Wake Model)，这也是基于求解RANS方程的CFD模型。2011年Ott等人提出了Fuga模型，它是一个线性RANS模型，把风力机等效成制动盘来模拟尾流。2011年Montavon等人]提出了WindModeller模型，采用k-ε湍流模型，制动盘模型模拟风力机，被应用于商业软件CFX。Xue等人提出了LES尾流模型，其对真实物体近壁面网格要求非常高，LES可以很好对尾流进行数值模拟的同时也付出了很大的计算代价。可见，现有的尾流计算方法不能完全避免精度低、形式复杂、编码难、计算低效等技术问题，现有技术尚无法解决这些问题。

发明内容

为解决现有技术的不足，本发明的目的在于提供一种风工程尾流计算方法，从而解决了现有技术不能完全避免尾流计算过程中精度低、形式复杂、编码难、计算低效的技术问题。

为了实现上述目标，本发明采用如下的技术方案：

一种风工程尾流计算方法，其特征在于，包括如下步骤：

步骤一：计算尾流膨胀系数；

步骤二：基于步骤一，计算Park模型下的尾流风速；

步骤三：修正Park模型下的初始尾流半径；

步骤四：基于步骤二和步骤三，计算修正后的Park模型尾流风速；

步骤五：基于步骤四，建立新的尾流模型，将新的尾流模型命名为Park-Gauss模型；

步骤六：计算Park-Gauss模型下的尾流风速。

前述的一种风工程尾流计算方法，其特征在于，在计算尾流膨胀系数过程中涉及到地表粗糙度对尾流的影响。

前述的一种风工程尾流计算方法，其特征在于，k＝0.5/ln(H/z₀)，其中，k为尾流膨胀系数，H为风力机轮毂高度，z₀为地表粗糙度。

前述的一种风工程尾流计算方法，其特征在于，其中，u_*为距离轮盘x处的尾流风速，r₀为初始尾流半径，a为轴向诱导因子。

前述的一种风工程尾流计算方法，其特征在于，所述步骤三中，修正Park模型下的初始尾流半径如下：

前述的一种风工程尾流计算方法，其特征在于，所述步骤四中，修正后的Park模型尾流风速的计算方式为：把步骤三中修正后的初始尾流半径代入步骤二中，得到修正后的尾流风速。

前述的一种风工程尾流计算方法，其特征在于，所述Park-Gauss模型为：其中，u为Park-Gauss模型下的尾流风速，r为径向上距离风轮中心的距离，r_x为风力机下游x位置处的尾流半径。

前述的一种风工程尾流计算方法，其特征在于，把步骤四中修正后的尾流风速代入步骤五中，得到Park-Gauss模型下的尾流风速。

本发明的有益之处在于：本发明修正后的初始尾流半径明显提高了尾流的预测精度；可以很好的模拟尾流区的风速，不仅在精度上与试验结果接近一致而且在径向分布上也更加符合真实流场，同时继承了工程模型的形式简单、容易编码、计算高效等优点。

附图说明

图1是本发明Park模型和Park-Gauss模型示意图；

图2是本发明Sexbierum风力机在下游2.5D处径向速度分布示意图；

图3是本发明Sexbierum风力机在下游5.5D处径向速度分布示意图；

图4是本发明Sexbierum风力机在下游8.0D处径向速度分布示意图；

图5是本发明Garrad Hassan风力机在下游5.0D处径向速度分布示意图；

图6是本发明Garrad Hassan风力机在下游7.5D处径向速度分布示意图；

图7是本发明Garrad Hassan风力机在下游7.5D处径向速度分布示意图；

图8是本发明整体流程示意图。

具体实施方式

以下结合附图和具体实施例对本发明作具体的介绍。

参照图8所示，本发明一种新的风工程尾流计算方法，包括如下步骤：

步骤1，计算尾流膨胀系数k，考虑地表粗糙度对尾流的影响，

k＝0.5/ln(H/z₀)

式中，H为轮毂高度，z₀为地表粗糙度。

步骤2，计算Park模型尾流风速。Park模型假设尾流是线性扩张，尾流风速随风向x的距离变化而变化，在径向方向风速呈常数分布。

式中u₀为来流风速，u_*是距离轮盘x处的尾流风速，r₀是初始尾流半径，a为轴向诱导因子。

步骤3，修正初始尾流半径，N.O.Jensen认为初始尾流半径r₀等于风轮半径r_d，考虑叶尖涡等现象的存在，真实初始半径肯定会比风轮半径r_d要大。本发明对其进行修正，采用紧靠风轮下游处的尾流半径为初始尾流半径。

步骤4，计算修正后的Park模型尾流风速，把步骤3中代入步骤2中计算尾流风速。

步骤5，建立新尾流模型Park-Gauss模型：

步骤6，计算新模型下的尾流风速，把步骤4计算的风速代入步骤5中。

下面根据各步骤一次进行展开说明。

1.计算尾流膨胀系数

根据风力机的轮毂高度z，当地地表粗糙度z₀，计算出尾流膨胀系数k

k＝0.5/ln(H/z₀) (1)

在不知预安装风机位置处的地表粗糙度时，我们可以依据Barthlmie等人建议，对于陆上风机k可取0.075，对于海上风机k可取0.05。

2.Park模型

Park模型假设尾流是线性扩张，尾流风速随风向x的距离变化而变化，在径向方向风速呈常数分布。根据来流风速u₀，a为轴向诱导因子，初始尾流半径r₀＝r_d，计算距离轮盘x处的尾流风速u_*

3.修正初始尾流半径

N.O.Jensen认为初始尾流半径r₀等于风轮半径r_d，考虑叶尖涡等现象的存在，真实初始半径肯定会比风轮半径r_d要大。本发明对其进行修正，采用紧靠风轮下游处的尾流半径为初始尾流半径。结合动量守恒原理和制动盘理论

u_*＝(1-2a)u₀ (4)

可以求出修正后的初始尾流半径为：

上式中，v₀是初始尾流半径处的风速，r_x是距离轮盘x处的尾流半径，其中r_x＝kx+r₀。(理论示意图见图1)。

4.建立Park-Gauss模型

Park尾流模型假设尾流速度只随风向x的距离变化，在径向方向风速为常数。径向速度分布呈帽子形状，因此Park模型又被成为top-hat模型。涡粘性模型假设尾流在径向呈对称且像高斯型分布。经试验测量在真实的流场中，在径向尾流风速呈多项式分布或高斯型分布。由此，本发明结合Park模型和涡粘性模型，提出了一种新的尾流模型：假设尾流区呈线性扩张，膨胀速率为k,尾流区径向风速呈高斯对称分布，使其更接近于真实尾流场，称为Park-Gauss模型。

假设Park-Gauss模型为：

式中A,B,C为待定系数，r_x为风力机下游x位置处的尾流半径，r为径向上距离风轮中心的距离。由于该Park-Gauss模型是基于Park模型和涡粘性模型假设出来的，所以满足以下几点假设：

假设当径向距离r扩展到尾流半径r_x时，在D点和E点处风速恢复到来流风速u₀。

假设两个模型的尾流半径相同，并且质量通量相等(即速度型下所包含的面积相等)。

由(7)，(8)得出：

把B＝0代入(9)得：

由于没有原函数，所以这里用幂级数展开式近似值来计算。

联立公式(10),(14)可求得A,B以及C：

再将(1.22)代入(1.13)得：

综上所述可得，Park-Gauss主要分为两步来实现：

预测步：u_*＝u₀[1-2a/(1+kx/r₀)²]

校正步：

5新模型验证

5.1用Sexbierum实测风场数据验证新模型有效性

Sexbierum风场布置了18台额定功率310kW,风轮直径D＝30.1m(近似看成D＝30m)，轮毂高度H＝35m，切入速度u_in＝5m/s,额定速度u＝14m/s,切出速度u_out＝20m/s。在风场周围安装了7台气象桅杆，用于测量风速、风向、气压、温度等。并对其中一台风力机下游2.5D、5.5D以及8.0D位置处尾流数据进行了详细的测量。来流风速u₀＝8.5m/s,风机推力系数C_T＝0.75，来流湍流强度I₀＝10％。

步骤1，计算尾流膨胀系数k，考虑地表粗糙度对尾流的影响，

k＝0.5/ln(H/z₀)

式中，H为轮毂高度，z₀为地表粗糙度。

步骤2，计算Park模型尾流风速。Park模型假设尾流是线性扩张，尾流风速随风向x的距离变化而变化，在径向方向风速呈常数分布。

式中u₀为来流风速，u_*是距离轮盘x处的尾流风速，r₀是初始尾流半径，a为轴向诱导因子。

步骤3，修正初始尾流半径，N.O.Jensen认为初始尾流半径r₀等于风轮半径r_d，考虑叶尖涡等现象的存在，真实初始半径肯定会比风轮半径r_d要大。本发明对其进行修正，采用紧靠风轮下游处的尾流半径为初始尾流半径。

步骤4，计算修正后的Park模型尾流风速，把步骤3中代入步骤2中计算尾流风速。

步骤5，建立新尾流模型Park-Gauss模型：

步骤6，计算新模型下的尾流风速，把步骤4计算的风速代入步骤5中。

通过图2-图4分析可得：

F-M2.5D-S:表示Sexbierum风场在风机下游2.5D处径向尾流实测数据；P-U2.5D-rd:表示Park模型在r₀＝r_d时风机下游2.5D处径向尾流风速；P-G2.5D-rd:表示Park-Gauss模型在r₀＝r_d时风机下游2.5D处径向尾流风速；P-U2.5D-r0:表示表示Park模型在时风机下游2.5D处径向尾流风速，其他以此类推。

通过图2-图4分析可得：

Park模型，在初始尾流半径r₀＝r_d条件下，低估了尾流效应，风速都比试验实测数据高。在风力机下游2.5D,5.5D,8.0D处，分别比试验数据高出43.9％，23.1％，15.2％。在修正后尾流半径条件下，对尾流效应低估情况有明显改善。在风力机下游2.5D,5.5D,8.0D处，分别比试验数据高出38.3％，17.2％，12.3％。

Park-Gauss模型，在初始尾流半径r₀＝r_d条件下，和Park模型一样都低估了尾流效应。但相比于后者，Park-Gauss模型对尾流预测都有很大的改善，不但径向风速趋势和试验数据近似一致，而且对尾流效应低估情况也有很明显的提高。在风力机下游2.5D,5.5D,8.0D处，分别比试验数据高出10.5％，9.5％，7.0％。在修正后尾流半径条件下，Park-Gauss模型在径向上风速走势更加接近于试验数据。在风力机下游2.5D,5.5D,8.0D处，几乎与试验数据是吻合的。

5.2用Garrad Hassan风洞试验结果验证新模型有效性

Garrad Hassan公司在1989年，开展了一系列风力机尾流的风洞实验研究[20]。试验选取的对象是风轮直径为0.27m的水平轴风力机模型。缩比尺度为1:160，对应的模拟风力机的风轮直径D＝43.2m，轮毂高度H＝50m，运行环境为地面粗糙度0.075m，2kmX2km的平坦地形。来流风速u0＝5.3m/s，叶尖速比λ1＝2.9，推力系数CT＝0.62。

步骤1，计算尾流膨胀系数k，考虑地表粗糙度对尾流的影响，

k＝0.5/ln(H/z₀)

式中，H为轮毂高度，z₀为地表粗糙度。

步骤2，计算Park模型尾流风速。Park模型假设尾流是线性扩张，尾流风速随风向x的距离变化而变化，在径向方向风速呈常数分布。

式中u₀为来流风速，u_*是距离轮盘x处的尾流风速，r₀是初始尾流半径，a为轴向诱导因子。

步骤3，修正初始尾流半径，N.O.Jensen认为初始尾流半径r₀等于风轮半径r_d，考虑叶尖涡等现象的存在，真实初始半径肯定会比风轮半径r_d要大。本发明对其进行修正，采用紧靠风轮下游处的尾流半径为初始尾流半径。

步骤4，计算修正后的Park模型尾流风速，把步骤3中代入步骤2中计算尾流风速。

步骤5，建立新尾流模型Park-Gauss模型：

步骤6，计算新模型下的尾流风速，把步骤4计算的风速代入步骤5中。

经过分析图5-图7可知：

Park模型，在初始尾流半径r₀＝r_d条件下，低估了尾流效应，风速都比试验实测数据高。在风力机下游5.0D,7.5D,10D处，分别比试验数据高出19.7％，10.0％，2.1％。在修正后尾流半径条件下，对尾流效应低估情况有明显改善。在风力机下游5.0D,7.5D,10D处，分别比试验数据高出17.5％，8.1％，0.6％。

Park-Gauss模型，在初始尾流半径r₀＝r_d条件下，Park-Gauss模型对尾流预测都有很大的改善，不但径向风速趋势和试验数据近似一致，而且对尾流效应低估情况也有很明显的提高。在风机下游5.0D，7.5D处稍稍高于试验数据，在10D处略低于试验数据。在修正后尾流半径条件下，Park-Gauss模型在径向上风速走势更加接近于试验数据。在风力机下游5.0D,7.5D,几乎与试验数据是吻合的，在10D处略低于试验数据。

可见，本发明经分别与风场实测数据和风洞试验结果比较，得出修正后的初始尾流半径明显提高了尾流的预测精度；新提出的Park-Gauss模型可以很好的模拟尾流区的风速，不仅在精度上与试验结果接近一致而且在径向分布上也更加符合真实流场，同时继承了工程模型的形式简单、容易编码、计算高效等优点。

以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解，上述实施例不以任何形式限制本发明，凡采用等同替换或等效变换的方式所获得的技术方案，均落在本发明的保护范围内。

Claims (2)

1.一种风工程尾流计算方法，其特征在于，包括如下步骤：

步骤一：计算尾流膨胀系数，，k＝0.5/ln(H/z₀)，其中，k为尾流膨胀系数，H为风力机轮毂高度，z₀为地表粗糙度；

步骤二：基于步骤一，计算Park模型下的尾流风速，其中，u_*为距离轮盘x处的尾流风速，r₀为初始尾流半径，a为轴向诱导因子，u₀为来流风速；

步骤三：修正Park模型下的初始尾流半径，修正Park模型下的初始尾流半径如下：r_d为分轮半径；

步骤四：基于步骤二和步骤三，计算修正后的Park模型尾流风速，修正后的Park模型尾流风速的计算方式为：把步骤三中修正后的初始尾流半径代入步骤二中，得到修正后的尾流风速；

步骤五：基于步骤四，建立新的尾流模型，将新的尾流模型命名为Park-Gauss模型，所述Park-Gauss模型为：其中，u为Park-Gauss模型下的尾流风速，r为径向上距离风轮中心的距离，r_x为风力机下游x位置处的尾流半径；

步骤六：计算Park-Gauss模型下的尾流风速，把步骤四中修正后的尾流风速代入步骤五中，得到Park-Gauss模型下的尾流风速。

2.根据权利要求1所述的一种风工程尾流计算方法，其特征在于，在计算尾流膨胀系数过程中涉及到地表粗糙度对尾流的影响。