CN106897486B - 考虑湍流强度影响的抛物线形风电机组尾流模型计算方法 - Google Patents
考虑湍流强度影响的抛物线形风电机组尾流模型计算方法 Download PDFInfo
- Publication number
- CN106897486B CN106897486B CN201710021657.XA CN201710021657A CN106897486B CN 106897486 B CN106897486 B CN 106897486B CN 201710021657 A CN201710021657 A CN 201710021657A CN 106897486 B CN106897486 B CN 106897486B
- Authority
- CN
- China
- Prior art keywords
- wake
- wind speed
- wind
- turbulence intensity
- wind wheel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 15
- 238000012821 model calculation Methods 0.000 title abstract description 5
- 238000004364 calculation method Methods 0.000 claims abstract description 20
- 238000011084 recovery Methods 0.000 claims abstract description 20
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 239000012467 final product Substances 0.000 claims description 2
- 230000004907 flux Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 description 5
- 230000007613 environmental effect Effects 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Wind Motors (AREA)
Abstract
本发明属于风电场微观选址技术领域，尤其涉及一种考虑湍流强度影响的抛物线形风电机组尾流模型计算方法，包括如下步骤：S1：假设尾流区域半径线性增长，且紧邻风轮后方尾流区域半径与风轮相同，求得风轮后轴向距离x位置处尾流区域半径r；S2：假设风轮后轴向距离x位置处风速沿径向均匀分布，根据质量守恒得到均匀风速v*；S3：根据质量守恒和均匀风速v*，得到轮毂高度水平面上任意点风速v；S4：考虑尾流区域湍流强度对尾流恢复系数k的影响，参考叶片旋转附加湍流强度计算经验公式，得到尾流恢复系数k的计算公式；S5：将尾流恢复系数k代入风速v计算公式，得到风电机组尾流区风速分布的经验尾流模型。
Description
技术领域
本发明属于风电场微观选址技术领域，尤其涉及一种考虑湍流强度影响的抛物线形风电机组尾流模型计算方法。
背景技术
能源短缺和环境污染问题日益严重，风能作为一种可再生清洁能源得到越来越多地应用。陆上和海上风电场逐渐建成投产，许多问题也随之而来，其中，下游风电机组发电量受上游风电机组尾流效应影响的问题极为突出。风电机组尾流效应主要导致尾流区风速降低，造成尾流区内风电机组发电量损失。风电场微观选址和风电场功率预测都需要考虑尾流效应的影响。
工程应用中，需要利用结构简单、计算时间短且计算精度较高的尾流区风速计算模型来量化尾流效应对发电量的影响。准确的尾流模型建立将有助于风电场的微观选址从而提高风电场经济效益以实现风资源的高效利用。商业软件采用的经验尾流模型主要有Jensen、Larsen、Ainslie模型等，其中，Jensen模型结构最简单，对于尾流区风速沿径向均匀分布的假设与实际相差较大，且未考虑环境湍流强度对尾流恢复的影响，低估了尾流区域风速恢复情况；Larsen模型假设风速沿径向呈高斯分布，但高估了风速恢复情况，尤其在近尾流区精确度较差；Ainslie模型采用涡漩粘性理论求解N-S方程得到尾流区风速分布，计算精度较高但计算时间长。
发明内容
针对上述问题，本发明提出了一种考虑湍流强度影响的抛物线形风电机组尾流模型计算方法。
一种考虑湍流强度影响的抛物线形风电机组尾流模型计算方法，包括如下步骤：
步骤1：假设尾流区域半径线性增长，且紧邻风轮后方尾流区域半径与风轮相同，求得风轮后轴向距离x位置处尾流区域半径r；
步骤2：假设风轮后轴向距离x位置处风速沿径向均匀分布，根据质量守恒得到均匀风速v*；
步骤3：根据质量守恒和均匀风速v*，得到轮毂高度水平面上任意点风速v；
步骤4：考虑尾流区域湍流强度对尾流恢复系数k的影响，用抛物线函数描述风速沿径向变化规律，参考叶片附加湍流强度计算经验公式，得到尾流恢复系数k的计算公式；
步骤5：将尾流恢复系数k代入风速v计算公式，得到风电机组尾流区风速分布的经验尾流模型。
优选的，所述经验尾流模型为抛物线尾流模型。
优选的，所述尾流区域半径r＝r0+kx：其中，r0为风轮半径。
优选的，所述尾流区域湍流强度包括环境湍流强度和叶片旋转附加湍流强度。
一种采用上述方法计算风电机组尾流区风速分布的经验尾流模型的方法，包括以下步骤：
步骤1：确定参考坐标系，以风轮中心为坐标原点，风轮旋转轴为x轴，径向为y轴，得到轮毂高度水平面各点位置坐标；
步骤2：根据来流风速，对照机组推力系数随风速变化曲线得到机组该工况下推力系数；
步骤3：将各输入参数代入抛物线尾流模型，得到轮毂高度水平面任意位置处风速值。
优选的，所述x轴方向平行于来流方向，所述y轴方向垂直于来流方向。
本发明的有益效果在于：
本发明在适当假设前提下，根据质量守恒，同时考虑到风速沿径向变化规律及环境湍流强度和叶片旋转附加湍流强度对尾流区风速恢复速率的影响，推导得到一种计算风电机组尾流区风速分布的经验尾流模型——抛物线尾流模型。
本发明同时考虑到风速沿径向变化规律及环境湍流强度和叶片旋转附加湍流强度对尾流区风速恢复速率的影响，建立一个可以快速准确计算风电机组尾流区风速分布的经验模型，为风电场微观选址计算尾流效应的影响提供参考。
附图说明
图1是控制体的动量守恒图；
图2是本发明的推导流程图；
图3是风电场布局及参考坐标系；
图4是风速值与测量值对比图。
具体实施方式
下面结合附图，对实施例作详细说明。
一种计算风电机组尾流区风速分布的经验尾流模型的方法，该方法包括风速计算公式和湍流强度计算公式，其中，输入参数有风电机组推力系数Ct、风轮直径D、机组后轴向距离x、径向距离y、环境湍流强度I0、来流风速v0，输出参数有叶片旋转附加湍流强度Ia、风速v。
一种计算风电机组尾流区风速分布的经验尾流模型的方法，如图2所示，包括以下步骤：
步骤1：假设尾流区域半径线性增长，且紧邻风轮后方尾流区域半径与风轮相同，求得风轮后轴向距离x位置处尾流区域半径r。
步骤2：假设风轮后轴向距离x位置处，风速沿径向均匀分布，根据质量守恒，推导得到均匀风速v*。
步骤3：考虑风速沿径向呈抛物线状变化规律，根据质量守恒和均匀风速v*，推导得到轮毂高度水平面上任意点风速v。
步骤4：考虑环境湍流强度和叶片旋转附加湍流强度对尾流恢复系数k的影响，参考叶片旋转附加湍流强度计算经验公式，得到尾流恢复系数k的计算公式。
步骤5：将参数k代入风速v计算公式，得到风电机组尾流区风速分布的经验尾流模型——抛物线尾流模型。
进一步，所述步骤1包括：
步骤11：引入尾流恢复系数k和x位置处尾流区域半径r，风轮半径为r0。
r＝r0+kx (1)
进一步，所述步骤2包括：
步骤21：以图1所示控制体为例，应用质量守恒，
其中，v1为紧邻风轮后尾流区风速，v0为来流风速，v*为x位置处均匀风速。
步骤22：由风轮动量理论得知，
v1＝(1-2a)v0 (3)
其中，a为轴向诱导因子，与推力系数Ct有如下关系
Ct＝4a(1-a) (4)
步骤23：联立公式(1)—(4)可得，
进一步，所述步骤3包括：
步骤31：引入风轮直径D＝2r0，风轮后x距离处无量纲半径r1。
步骤32：假设风轮后轴向距离为x处，风速沿径向(y轴方向)呈抛物线状分布
v＝Ay2+B (7)
当y＝r1时风速恢复到来流风速
Ar1 2+B＝v0 (8)
步骤33：根据质量守恒，风轮后轴向距离为x处，均匀风速与抛物线状分布有相同的通量
步骤34：将公式(6)—(8)代入公式(9)求解得
进一步，所述步骤4包括：
步骤41：尾流区域湍流强度包括两部分，环境湍流强度I0和叶片旋转附加湍流强度Ia
步骤42：叶片旋转附加湍流强度Ia由经验公式(12)求得
步骤43：尾流恢复系数k由经验公式(13)求得，其中，κ为经验常数取值0.4，可根据具体湍流强度变化情况进行修正。
k≈κ×I (13)
进一步，所述步骤5包括：
步骤51：将尾流恢复系数k代入公式(10)，得到一种风电机组尾流区风速分布的经验尾流模型——抛物线尾流模型，如公式(14)所示。
实施例1
一种计算风电机组尾流区风速分布的经验尾流模型的应用，包括以下步骤：
步骤1：确定参考坐标系，以风轮中心为坐标原点，风轮旋转轴为x轴(平行于来流方向)，径向(垂直于来流方向)为y轴，得到轮毂高度水平面各点位置坐标。
步骤2：根据来流风速，对照机组推力系数随风速变化曲线得到机组该工况下推力系数。
步骤3：将各输入参数代入抛物线尾流模型，可计算得到轮毂高度水平面任意位置处风速值。
实施例2
Sexbierum风场是位于荷兰北部的平坦地形的陆上风电场，总装机5.4MW，拥有18台额定功率为310kW的机组。风电机组风轮直径30m，轮毂高度30m，切入风速、额定风速、切出风速分别为5m/s、14m/s、20m/s。风场机组的布局如图3所示。正北方向为180°，T18与T27连线方向为231°，在T18后方2.5D、5.5D、8D位置处分别立有测风塔进行尾流区风速分布测量。
应用本发明进行风电场中单台风电机组尾流区风速计算，包括以下步骤：
步骤1：确定参考坐标系，以风轮中心为坐标原点，风轮旋转轴为x轴(平行于来流方向—231°)，径向(垂直于来流方向)为y轴，得到轮毂高度水平面各点位置坐标，如图3。
步骤2：该算例中来流风速为8.5m/s，对应机组推力系数曲线得知，风轮推力系数为0.75。
步骤3：将各输入参数代入抛物线尾流模型，可计算得到轮毂高度水平面任意位置处风速值。
进一步，所述步骤3包括：
步骤31：分析测风塔测风数据得出轮毂高度处环境湍流强度I0为10％。
步骤32：输入风轮半径风轮直径30m,半径15m。
步骤33：确定计算域范围，x和y的取值范围，-2＜x＜2，0＜y＜300。
步骤34：将计算域范围进行网格划分，得到离散点，分别代入抛物线尾流模型计算得到个点风速。
提取机组后2.5D、5.5D、8D位置处风速计算值与测量值进行对比，如图4所示，抛物线尾流模型能较准确地预测单台机组尾流区风速分布，可用于风电场微观选址、风电功率短期预测等工程应用。
此实施例仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。
Claims (3)
1.一种考虑湍流强度影响的抛物线形风电机组尾流模型计算方法，其特征在于，包括如下步骤：
步骤1：假设尾流区域半径线性增长，且紧邻风轮后方尾流区域半径与风轮相同，求得风轮后轴向距离x位置处尾流区域半径r；
步骤2：假设风轮后轴向距离x位置处风速沿径向均匀分布，根据质量守恒得到均匀风速v*；
步骤3：根据质量守恒和均匀风速v*，以抛物线函数描述风速沿径向分布规律，得到轮毂高度水平面上任意点风速v；
步骤4：考虑尾流区域湍流强度对尾流恢复系数k的影响，参考叶片旋转附加湍流强度计算经验公式，得到尾流恢复系数k的计算公式；
步骤5：将尾流恢复系数k代入风速v计算公式，得到风电机组尾流区风速分布的经验尾流模型；
所述步骤3包括：
步骤31：引入风轮直径D＝2r0，风轮后x距离处无量纲半径r1；
步骤32：假设风轮后轴向距离为x处，风速沿径向呈抛物线状分布
v＝Ay2+B (7)
当y＝r1时风速恢复到来流风速
Ar1 2+B＝v0 (8)
步骤33：根据质量守恒，风轮后轴向距离为x处，均匀风速与抛物线状分布有相同的通量
步骤34：将公式(6)—(8)代入公式(9)求解得
所述步骤4包括：
步骤41：尾流区域湍流强度包括两部分，环境湍流强度I0和叶片旋转附加湍流强度Ia
步骤42：叶片旋转附加湍流强度Ia由经验公式(12)求得
步骤43：尾流恢复系数k由经验公式(13)求得，其中，κ为经验常数取值0.4，可根据具体湍流强度变化情况进行修正；
k≈κ×I (13)
所述步骤5包括：
步骤51：将尾流恢复系数k代入公式(10)，得到一种风电机组尾流区风速分布的经验尾流模型——抛物线尾流模型，如公式(14)所示；
2.一种采用权利要求1所述的方法计算风电机组尾流区风速分布的方法，其特征在于，包括以下步骤：
步骤1：确定参考坐标系，以风轮中心为坐标原点，风轮旋转轴为x轴，径向为y轴，得到轮毂高度水平面各点位置坐标；
步骤2：根据来流风速，对照机组推力系数随风速变化曲线得到机组当前工况下推力系数；
步骤3：将各输入参数代入抛物线尾流模型，得到轮毂高度水平面任意位置处风速值。
3.根据权利要求2所述方法，其特征在于，所述x轴方向平行于来流方向，所述y轴方向垂直于来流方向。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710021657.XA CN106897486B (zh) | 2017-01-12 | 2017-01-12 | 考虑湍流强度影响的抛物线形风电机组尾流模型计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710021657.XA CN106897486B (zh) | 2017-01-12 | 2017-01-12 | 考虑湍流强度影响的抛物线形风电机组尾流模型计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106897486A CN106897486A (zh) | 2017-06-27 |
| CN106897486B true CN106897486B (zh) | 2020-07-07 |
Family
ID=59198702
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710021657.XA Active CN106897486B (zh) | 2017-01-12 | 2017-01-12 | 考虑湍流强度影响的抛物线形风电机组尾流模型计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106897486B (zh) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108108562B (zh) * | 2017-12-27 | 2021-02-19 | 华北电力大学 | 一种基于高斯分布的风电机组尾流的解析建模方法 |
| CN109325248A (zh) * | 2018-02-11 | 2019-02-12 | 同济大学 | 风电场内立体式风速分布关联图建立方法 |
| CN108536907B (zh) * | 2018-03-01 | 2021-11-30 | 华北电力大学 | 一种基于简化动量定理的风电机组远场尾流解析建模方法 |
| CN108491617B (zh) * | 2018-03-19 | 2022-04-19 | 华北电力大学 | 一种基于尾流边界膨胀的bp模型的简化模型 |
| CN108509718B (zh) * | 2018-03-29 | 2022-02-25 | 华北电力大学 | 一种基于质量守恒的远场尾流二维解析模型 |
| CN110516328A (zh) * | 2019-08-14 | 2019-11-29 | 中国科学院寒区旱区环境与工程研究所 | 一种建立风电机出力仿真模型的方法 |
| CN111651870B (zh) * | 2020-05-20 | 2023-07-25 | 中国能源建设集团广东省电力设计研究院有限公司 | 一种基于改进Jensen模型的风力机尾流计算方法及系统 |
| CN112001131B (zh) * | 2020-07-28 | 2022-11-11 | 河海大学 | 一种改进风力机尾流叠加计算方法 |
| CN111985097B (zh) * | 2020-08-13 | 2023-01-06 | 大唐三亚未来能源研究所有限公司 | 一种考虑海浪波高影响的海上风电机组尾流计算方法 |
| CN112347611A (zh) * | 2020-10-15 | 2021-02-09 | 华北电力大学 | 一种风力机远场尾流流向湍流度计算方法 |
| CN112199908B (zh) * | 2020-10-16 | 2022-06-28 | 中国船舶重工集团海装风电股份有限公司 | 一种基于流体力学的风电机组偏航控制尾流模型修正方法 |
| CN112434252B (zh) * | 2020-11-27 | 2022-12-09 | 中国华能集团清洁能源技术研究院有限公司 | 一种考虑风电场局地环境因素的尾流计算方法 |
| CN113153636B (zh) * | 2021-04-14 | 2022-11-22 | 中国科学院力学研究所 | 一种促进风力发电机尾流恢复的流动控制装置和方法 |
| CN114021488B (zh) * | 2021-10-20 | 2023-02-28 | 中国大唐集团新能源科学技术研究院有限公司 | 风轮前风能估算方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101852172A (zh) * | 2010-03-09 | 2010-10-06 | 山东科技大学 | 一种考虑风电场内尾流影响的风电机组输入风速计算方法 |
| CN103884485A (zh) * | 2014-02-25 | 2014-06-25 | 国家电网公司 | 一种基于多尾流模型的风机尾流分析方法 |
| CN103942608A (zh) * | 2014-03-28 | 2014-07-23 | 华北电力大学 | 一种基于尾流模型的风电场优化调度方法 |
-
2017
- 2017-01-12 CN CN201710021657.XA patent/CN106897486B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101852172A (zh) * | 2010-03-09 | 2010-10-06 | 山东科技大学 | 一种考虑风电场内尾流影响的风电机组输入风速计算方法 |
| CN103884485A (zh) * | 2014-02-25 | 2014-06-25 | 国家电网公司 | 一种基于多尾流模型的风机尾流分析方法 |
| CN103942608A (zh) * | 2014-03-28 | 2014-07-23 | 华北电力大学 | 一种基于尾流模型的风电场优化调度方法 |
Non-Patent Citations (2)
| Title |
|---|
| Comparison and validation of wake models based on field measurements with lidar;Li Li 等;《IET 5th Renewable Power Generation Conference》;20160922;第1-6页 * |
| WIND TURBINE WAKE VERTICAL DISTRIBUTIONS CONSIDERING DIFFERENT INFLOW SHEAR INDICES;Yansong Cui 等;《IET 4th Renewable Power Generation Conference》;20151018;第1-6页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106897486A (zh) | 2017-06-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106897486B (zh) | 考虑湍流强度影响的抛物线形风电机组尾流模型计算方法 | |
| CN109376389B (zh) | 一种基于2D_k Jensen模型的三维尾流数值模拟方法 | |
| Bossanyi | Combining induction control and wake steering for wind farm energy and fatigue loads optimisation | |
| Sanchez Gomez et al. | The effect of wind direction shear on turbine performance in a wind farm in central Iowa | |
| He et al. | A novel three-dimensional wake model based on anisotropic Gaussian distribution for wind turbine wakes | |
| CN104794347B (zh) | 一种风电机组机舱传递函数分区拟合方法 | |
| Schepers et al. | Analysis of wake measurements from the ECN wind turbine test site wieringermeer, EWTW | |
| CN103884485A (zh) | 一种基于多尾流模型的风机尾流分析方法 | |
| CN106548414B (zh) | 一种海上风电场发电量计算方法 | |
| CN103244354B (zh) | 风力发电机组功率曲线自适应优化方法 | |
| CN103886179B (zh) | 一种基于尾流效应分群的风电场风机聚合方法 | |
| CN104675629B (zh) | 一种变速风力发电机组的最大风能捕获方法 | |
| Hassan et al. | Predicting the amount of energy generated by aWind turbine based on the weather data | |
| Mazon et al. | An assessment of the sea breeze energy potential using small wind turbines in peri-urban coastal areas | |
| CN106951977B (zh) | 一种基于尾流效应的风速预测模型的构建方法 | |
| CN115544884A (zh) | 一种基于数据驱动的大型风电场尾流快速计算方法及系统 | |
| CN105023099B (zh) | 一种考虑湍流强度的风力发电机出力评估方法 | |
| Hu et al. | An experimental investigation on the effects of turbine rotation directions on the wake interference of wind turbines | |
| Maeda et al. | Investigation of wake characteristics of a Horizontal Axis Wind Turbine in vertical axis direction with field experiments | |
| CN106919730B (zh) | 一种采用风速衰减因子的风电场尾流确定方法 | |
| Meister et al. | Simulation of a 5MW wind turbine in an atmospheric boundary layer | |
| Guggeri et al. | Simulation of a 7.7 MW onshore wind farm with the Actuator Line Model | |
| CN106779202B (zh) | 一种考虑空气湿度的风电功率预测方法 | |
| CN109753759A (zh) | 一种基于等效功率的风轮等效风速计算方法 | |
| Filipowicz et al. | Study of building integrated wind turbines operation on the example of Center of Energy AGH |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |