CN106548414A - 一种海上风电场发电量计算方法 - Google Patents
一种海上风电场发电量计算方法 Download PDFInfo
- Publication number
- CN106548414A CN106548414A CN201610995568.0A CN201610995568A CN106548414A CN 106548414 A CN106548414 A CN 106548414A CN 201610995568 A CN201610995568 A CN 201610995568A CN 106548414 A CN106548414 A CN 106548414A
- Authority
- CN
- China
- Prior art keywords
- wind
- electric field
- marine
- wind electric
- marine wind
- 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.)
- Granted
Links
- 230000005684 electric field Effects 0.000 title claims abstract description 117
- 238000000205 computational method Methods 0.000 title claims abstract description 17
- 230000003746 surface roughness Effects 0.000 claims abstract description 29
- 230000000694 effects Effects 0.000 claims abstract description 11
- 238000000034 method Methods 0.000 claims abstract description 10
- 238000005755 formation reaction Methods 0.000 claims description 18
- 230000015572 biosynthetic process Effects 0.000 claims description 15
- 230000007935 neutral effect Effects 0.000 claims description 11
- 230000005611 electricity Effects 0.000 claims description 7
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 238000005315 distribution function Methods 0.000 claims description 6
- 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 3
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003653 coastal water Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种海上风电场发电量计算方法,包括以下步骤:收集测风数据;求解海面粗糙度;利用理查森数法求解理查森数并判定海上风电场热稳定性;求解莫宁长度;对风轮廓模型进行修正并求解修正后的风速数据;对Jensen尾流模型进行修正并对海上风电场不同热稳定性赋予对应的不同的尾流耗散系数中的常数;求解单台机组尾流影响下的下游机组风轮前的风速;求解多台机组尾流影响下的下游机组风轮前的风速;求解海上风电场不同热稳定性下单台机组输出功率;分别计算海上风电场不同热稳定性下单台机组年总发电量,以及海上风电场年总发电量。充分考虑近海风电场不同区域、不同时间风能分布的差异性,实现快速准确地计算海上风电场发电量。
Description
技术领域
本发明涉及一种发电量计算方法,特别是涉及一种海上风电场发电量计算方法,属于海上风电场风资源计算技术领域。
背景技术
我国海上风能资源丰富,靠近负荷中心,但近年来海上风电产业发展缓慢。随着技术的进步、政策的不断落实,大规模的海上风电工程相继开工建设,海上风电将迎来快速发展时期。
随着近海风电产业的发展,提高海上风能评估准确性有着重要的意义。海面粗糙度和大气热稳定性是海上风能影响的主要因素,海面粗糙度又称为海面空气动力学粗糙长度,其概念是从陆面上对数风廓线理论中的粗糙度延伸应用到海面而得来的,定义为风速等于零的高度,目前海上风能资源计算评估中还没有充分考虑粗糙度的变化,粗糙度取为同一值,影响了风能计算的准确性,造成风能计算误差。
大气热稳定性是指受温度分布状况影响大气垂直方向运动受到抑制或加强的程度。热稳定性对海上风电场空气动力场有重要的影响,但目前针对风电场微观尺度空气动力场的研究大部分是在中性大气条件下开展的,很少考虑风电场的热稳定性,降低了风电场空气动力场的预测精度。
发明内容
本发明的主要目的在于,克服现有技术中的不足,提供一种海上风电场发电量计算方法,充分考虑近海风电场不同区域、不同时间风能分布的差异性,实现快速准确地计算海上风电场发电量。
为了达到上述目的,本发明所采用的技术方案是:
一种海上风电场发电量计算方法,包括以下步骤:
1)收集海上风电场所建的测风塔的一年的测风数据,所述测风数据包括同一测风年内测风塔所测的不同高度的风速、风向及温度数据,其中不同高度包括距离海面10m高度以及测风塔轮毂高度H;
2)利用海面10m高度的风速数据以及海面粗糙度计算公式,求解海面粗糙度z0;
3)利用理查森数法求解理查森数Ri,并根据理查森数Ri判定海上风电场热稳定性;
4)根据理查森数Ri求解莫宁长度L;
5)利用海面粗糙度z0以及莫宁长度L对风轮廓模型进行修正,并求解修正后的风速数据V0;
6)对Jensen尾流模型进行修正,并对海上风电场不同热稳定性赋予对应的不同的尾流耗散系数k中的尾流耗散常数kw;
7)求解单台机组尾流影响下的下游机组风轮前的风速vx;
8)求解多台机组尾流影响下的下游机组风轮前的风速vj(t);
9)根据风电场风力发电机组的功率曲线以及步骤8)求解得到的多台机组影响下的下游机组风轮前的风速vj(t),求解海上风电场不同热稳定性下单台机组输出功率E(v);
10)分别计算海上风电场不同热稳定性下单台机组年总发电量Ej,以及海上风电场年总发电量E。
本发明进一步设置为:所述步骤2)中的海面10m高度的海面粗糙度计算公式为,
其中,z0为海面粗糙度,s为冯卡门常数、取值为0.35,U10为海面10m高度处风速。
本发明进一步设置为:所述步骤3)中的理查森数Ri的计算公式为,
ΔT=T2-T1 (4)
Δu=u2-u1 (5)
其中,g为重力加速度,I为z1和z2的算术平方根,z1和z2分别为上、下两个气层的高度,T1和T2分别为上、下两个气层的温度,T为上、下两个气层的温度T1和T2的平均值,ΔT为上、下两个气层的温差值,u1和u2分别为上、下两个气层的速度,Δu为上、下两个气层的速度差;
所述步骤3)中的判定海上风电场热稳定性为,
当理查森数Ri的数值范围为Ri>0.2,则判定海上风电场热稳定性为稳定;
当理查森数Ri的数值范围为-0.6<Ri<=0.2,则判定海上风电场热稳定性为中性;
当理查森数Ri的数值范围为-2.5<Ri<=-0.6,则判定海上风电场热稳定性为不稳定;
当理查森数Ri的数值范围为Ri<=-2.5,则判定海上风电场热稳定性为极不稳定。
本发明进一步设置为:所述步骤4)中的莫宁长度L的计算公式为,
其中,L为莫宁长度,I为z1和z2的算术平方根,Ri为理查森数。
本发明进一步设置为:所述步骤5)中的风轮廓模型修正为,
其中,u*为摩擦速度,s为冯卡门常数,取值为0.35;z为垂直方向上风轮廓对应的高度值,z0为海面粗糙度,ψm为风速的通用函数,L为莫宁长度;
摩擦速度u*的求解公式为u* 2=C10U10 2,其中C10为阻力系数,C10可以根据吴京公式求解,U10为海面10m高度处风速;
其中,风速的通用函数ψm的计算公式为,
当海上风电场热稳定性为稳定时,
当海上风电场热稳定性为中性时,
当海上风电场热稳定性为不稳定或极不稳定时,
y=[1-(16·z/L)]1/4 (11)
其中,y为公式(10)中单列出的公共项;
所述步骤5)中的求解修正后的风速数据V0为,根据式(7)及风速的通用函数ψm的计算公式,计算考虑了海上风电场热稳定性和海面粗糙度后的轮毂高度z=H处的风速数据u(z)。
本发明进一步设置为:所述步骤6)对Jensen尾流模型进行修正,对海上风电场不同热稳定性赋予对应的不同的尾流耗散系数k中的尾流耗散常数kw,具体为,
6-1)根据步骤3)判定的海上风电场热稳定性,将步骤5)求解得到的修正后的风速数据V0对应分为稳定、中性、不稳定和极不稳定共四组;
6-2)计算尾流耗散系数k,k=kw(σG+σ0)/v0;
其中,kw为尾流耗散常数,σG和σ0分别为风电机组产生的湍流和自然湍流的均方差,v0为自然风速;
6-3)对海上风电场不同热稳定性赋予对应的不同的尾流耗散系数k中的尾流耗散常数kw,包括,
当海上风电场热稳定性为稳定时,赋予尾流耗散系数k中的尾流耗散常数kw值为0.098;
当海上风电场热稳定性为中性时,赋予尾流耗散系数k中的尾流耗散常数kw值为0.048;
当海上风电场热稳定性为不稳定时,赋予尾流耗散系数k中的尾流耗散常数kw值为0.051;
当海上风电场热稳定性为极不稳定时,赋予尾流耗散系数k中的尾流耗散常数kw值为0.044。
本发明进一步设置为:所述步骤7)求解单台机组尾流影响下的下游机组风轮前的风速vx,具体为,
7-1)根据动量理论有,
其中,ρ为空气密度,R和Rw分别为叶轮半径和尾流半径,vx为受尾流影响的风速,vT为通过叶片的风速;
7-2)根据推力系数公式求解得到自然风速v0、通过叶片的风速vT和风电机组的推力系数CT具有如下关系,
vT=v0(1-CT)1/2 (15)
7-3)单台机组尾流影响下的下游机组风轮前的风速vx的计算公式为,
其中,X为两个风电机组的间距。
本发明进一步设置为:所述步骤8)求解多台机组尾流影响下的下游机组风轮前的风速vj(t),其计算公式为,
其中,vj(t)为作用在任意一台机组上的风速,vj0(t)为没有经过任何塔影影响作用在第j台风力发电机组上的风速、即自由流风速,vmj(t)为考虑机组间尾流效应时第m台风力发电机组作用在第j台风力发电机组上的尾流风速,表示在第j台风力发电机组处第m台风力发电机组的投影面积与第j台风力发电机组面积Arot-j的比,n为风力发电机组总台数,t表示时刻。
本发明进一步设置为:所述步骤10)分别计算海上风电场不同热稳定性下单台机组年总发电量Ej,以及海上风电场年总发电量E,具体为,
10-1)采用威布尔分布描述平均风速变化,
则平均风速的概率密度函数f(v)为,
平均风速的累积分布函数F(v)为,
F(v)=1-exp(-(v/c)p) (19)
其中,p为决定分布范围的形状参数,c为决定位置的尺度参数,v为测风数据中的风速;
10-2)将步骤5)修正后的风速数据V0带入公式(18)中的v,分别求解海上风电场不同热稳定性下的平均风速的概率密度函数;
将步骤5)修正后的风速数据V0带入公式(19)中的v,分别求解海上风电场不同热稳定性下的平均风速的累积分布函数;
10-3)计算单台机组年总发电量Ej,
其中,p(θ)为对应角度θ的风向频率,vin为风电机组的切入风速,vout为风电机组的切出风速,N(v)为相应风速等级出现的全年累计小时数,E(v)为速度v通过风电场风力发电机组的功率曲线得到的单台机组输出功率,f(v)为通过威布尔分布的平均风速的概率密度函数得到的频率;
将步骤5)修正后的风速数据V0带入公式(20)中的v,分别求解海上风电场不同热稳定性下的单台机组年总发电量Ej;
10-4)计算海上风电场年总发电量E,
其中,n为海上风电场的风电机组总台数。
与现有技术相比,本发明具有的有益效果是:
通过测风数据的收集,选择随风速变化的海面粗糙度计算公式,利用理查森数法判定海上风电场热稳定性,充分考虑了海面粗糙度随时间变化特征以及不同热稳定性情况,不仅对风轮廓模型进行修正,还对Jensen尾流模型进行修正,给出不同热稳定性下合适的尾流耗散系数中的尾流耗散常数,并利用修正后的风速数据求解海上风电场年总发电量。不仅体现出近海风电场不同区域、不同时间风能分布的差异性,使得计算得到的风速和风能更可靠,而且大大提高了既有海上风电场及其临近区域风资源分布规律模拟结果的准确性,为后续既有海上风电场临近场区发电量估算和选址工作奠定基础,可对近海风电场微观选址、短期风功率预测等有一定的指导意义,在工程中有较好的应用前景。
上述内容仅是本发明技术方案的概述,为了更清楚的了解本发明的技术手段,下面结合附图对本发明作进一步的描述。
附图说明
图1为本发明一种海上风电场发电量计算方法的流程图;
图2为本发明的不同热稳定性下风轮廓示意图;
图3为本发明修正的Jensen尾流模型;
图4为本发明的风机功率特性曲线。
具体实施方式
下面结合说明书附图,对本发明作进一步的说明。
本发明提供一种海上风电场发电量计算方法,如图1所示,包括以下步骤:
1)收集海上风电场所建的测风塔的一年的测风数据,所述测风数据包括同一测风年内测风塔所测的不同高度的风速、风向及温度数据,其中不同高度包括距离海面10m高度以及测风塔轮毂高度H;这里的海上风电场包括既有海上风电场或其临近区域。
2)利用海面10m高度的风速数据及海面粗糙度计算公式,求解海面粗糙度z0;
考虑海面粗糙度随时间变化特征,海面粗糙度受风况和海况的影响,风速越大粗糙度也越大,所以海面粗糙度不是同一值,海面10m高度的海面粗糙度计算公式为,
其中,z0为海面粗糙度,s为冯卡门常数、取值为0.35,U10为海面10m高度处风速。
3)利用理查森数法求解理查森数Ri,并根据理查森数Ri判定海上风电场热稳定性;
考虑海上风电场热稳定性随时间变化特征,以及海上风电场本身特性,理查森数法综合了湍流激发的热力因子和动力因子的作用,可更多的反应湍流状况信息,能够准确判定不同边界条件下热稳定性,所以采用梯度的理查森数方法,其理查森数Ri的计算公式为,
ΔT=T2-T1 (4)
Δu=u2-u1 (5)
其中,g为重力加速度,I为z1和z2的算术平方根,z1和z2分别为上、下两个气层的高度,T1和T2分别为上、下两个气层的温度,T为上、下两个气层的温度T1和T2的平均值,ΔT为上、下两个气层的温差值,u1和u2分别为上、下两个气层的速度,Δu为上、下两个气层的速度差。
根据计算得到的理查森数Ri来判定海上风电场热稳定性,如表1所示,即
当理查森数Ri的数值范围为Ri>0.2,则判定海上风电场热稳定性为稳定;
当理查森数Ri的数值范围为-0.6<Ri<=0.2,则判定海上风电场热稳定性为中性;
当理查森数Ri的数值范围为-2.5<Ri<=-0.6,则判定海上风电场热稳定性为不稳定;
当理查森数Ri的数值范围为Ri<=-2.5,则判定海上风电场热稳定性为极不稳定。
Ri数值范围 | 稳定度情况 |
Ri<=-2.5 | 极不稳定 |
-2.5<Ri<=-0.6 | 不稳定 |
-0.6<Ri<=0.2 | 中性 |
Ri>0.2 | 稳定 |
表1
4)根据理查森数Ri求解莫宁长度L;
莫宁长度L的计算公式为,
其中,L为莫宁长度,I为高度z1和z2的算术平方根,Ri为理查森数。
5)利用海面粗糙度z0以及莫宁长度L对风轮廓模型进行修正,并求解修正后的风速数据V0;
考虑了海面粗糙度以及海上风电场热稳定性相互作用,如图2所示,风轮廓模型修正为,
其中,u*为摩擦速度,s为冯卡门常数、取值为0.35,z为垂直方向上风轮廓对应的高度值,z0为海面粗糙度,ψm为风速的通用函数,L为莫宁长度;
摩擦速度u*的求解公式为u* 2=C10U10 2,其中C10为阻力系数,C10可以根据吴京公式求解,U10为海面10m高度处风速;
其中,风速的通用函数ψm的计算公式为,
当海上风电场热稳定性为稳定时,
当海上风电场热稳定性为中性时,
当海上风电场热稳定性为不稳定或极不稳定时,
y=[1-(16·z/L)]1/4 (11)
其中,y为公式(10)中单列出的公共项;
根据式(7)及风速的通用函数ψm的计算公式,计算考虑了海上风电场热稳定性和海面粗糙度后的轮毂高度z=H处的风速数据u(z)即为修正后的风速数据V0。
6)对Jensen尾流模型进行修正,并对海上风电场不同热稳定性赋予对应的不同的尾流耗散系数k中的尾流耗散常数kw;
6-1)根据步骤3)判定的海上风电场热稳定性,将步骤5)求解得到的修正后的风速数据V0对应分为稳定、中性、不稳定和极不稳定共四组;
6-2)考虑海上风电场不同热稳定性下,风电场风力发电机组之间尾流耗散特性不同,修正Jensen尾流模型,如图3所示;
计算尾流耗散系数k,k=kw(σG+σ0)/v0;
其中,kw为常数,σG和σ0分别为风电机组产生的湍流和自然湍流的均方差,v0为自然风速;
6-3)对海上风电场不同热稳定性赋予对应的不同的尾流耗散系数k中的尾流耗散常数kw,如表2所示,包括,
当海上风电场热稳定性为稳定时,赋予尾流耗散系数k中的尾流耗散常数kw值为0.098;
当海上风电场热稳定性为中性时,赋予尾流耗散系数k中的尾流耗散常数kw值为0.048;
当海上风电场热稳定性为不稳定时,赋予尾流耗散系数k中的尾流耗散常数kw值为0.051;
当海上风电场热稳定性为极不稳定时,赋予尾流耗散系数k中的尾流耗散常数kw值为0.044。
表2
7)求解单台机组尾流影响下的下游机组风轮前的风速vx;
7-1)根据动量理论有,
其中,ρ为空气密度,R和Rw分别为叶轮半径和尾流半径,vx为受尾流影响的风速,vT为通过叶片的风速;
7-2)根据推力系数公式求解得到自然风速v0、通过叶片的风速vT和风电机组的推力系数CT具有如下关系,
vT=v0(1-CT)1/2 (15)
7-3)单台机组尾流影响下的下游机组风轮前的风速vx的计算公式为,
其中,X为两个风电机组的间距。
8)求解多台机组尾流影响下的下游机组风轮前的风速vj(t);
根据动量守恒定律,其计算公式为,
其中,vj(t)为作用在任意一台机组上的风速,vj0(t)为没有经过任何塔影影响作用在第j台风力发电机组上的风速、即自由流风速,vmj(t)为考虑机组间尾流效应时第m台风力发电机组作用在第j台风力发电机组上的尾流风速,表示在第j台风力发电机组处第m台风力发电机组的投影面积与第j台风力发电机组面积Arot-j的比,n为风力发电机组总台数,t表示时刻。
9)根据风电场风力发电机组的功率曲线以及步骤8)求解得到的多台机组影响下的下游机组风轮前的风速vj(t),求解海上风电场不同热稳定性下单台机组输出功率E(v),如图4所示。
10)分别计算海上风电场不同热稳定性下单台机组年总发电量Ej,以及海上风电场年总发电量E;
10-1)采用威布尔分布描述平均风速变化,
则平均风速的概率密度函数f(v)为,
平均风速的累积分布函数F(v)为,
F(v)=1-exp(-(v/c)p) (19)
其中,p为决定分布范围的形状参数,c为决定位置的尺度参数,v为测风数据中的风速;
10-2)将步骤5)修正后的风速数据V0带入公式(18)中的v,分别求解海上风电场不同热稳定性下的平均风速的概率密度函数;
将步骤5)修正后的风速数据V0带入公式(19)中的v,分别求解海上风电场不同热稳定性下的平均风速的累积分布函数;
10-3)计算单台机组年总发电量Ej,
其中,p(θ)为对应角度θ的风向频率,vin为风电机组的切入风速,vout为风电机组的切出风速,N(v)为相应风速等级出现的全年累计小时数,E(v)为速度v通过风电场风力发电机组的功率曲线得到的单台机组输出功率,f(v)为通过威布尔分布的平均风速的概率密度函数得到的频率;
将步骤5)修正后的风速数据V0带入公式(20)中的v,分别求解海上风电场不同热稳定性下的单台机组年总发电量Ej;
10-4)计算海上风电场年总发电量E,
其中,n为海上风电场的的风力发电机组总台数。
风速分布是描述风速概率分布对时间的函数,即概率密度,数学上用来描述概率密度分布的函数很多,常用的为威布尔分布和瑞利分布;本发明实施例采用的是两参数的威布尔分布,其中形状参数和尺度参数均为正值。整个海上风电场的年总发电量可以看成是各单台风力机不同热稳定性下的年总发电量的代数和,所以先求4种热稳定性下单台风电机组的年总发电量,再计算所有机组年发电量总和,从而实现对整个海上风电场年总发电量的估算。
以上显示和描述了本发明的基本原理、主要特征及优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (9)
1.一种海上风电场发电量计算方法,其特征在于,包括以下步骤:
1)收集海上风电场所建的测风塔的一年的测风数据,所述测风数据包括同一测风年内测风塔所测的不同高度的风速、风向及温度数据,其中不同高度包括距离海面10m高度以及测风塔轮毂高度H;
2)利用海面10m高度的风速数据以及海面粗糙度计算公式,求解海面粗糙度z0;
3)利用理查森数法求解理查森数Ri,并根据理查森数Ri判定海上风电场热稳定性;
4)根据理查森数Ri求解莫宁长度L;
5)利用海面粗糙度z0以及莫宁长度L对风轮廓模型进行修正,并求解修正后的风速数据V0;
6)对Jensen尾流模型进行修正,并对海上风电场不同热稳定性赋予对应的不同的尾流耗散系数k中的尾流耗散常数kw;
7)求解单台机组尾流影响下的下游机组风轮前的风速vx;
8)求解多台机组尾流影响下的下游机组风轮前的风速vj(t);
9)根据风电场风力发电机组的功率曲线以及步骤8)求解得到的多台机组影响下的下游机组风轮前的风速vj(t),求解海上风电场不同热稳定性下单台机组输出功率E(v);
10)分别计算海上风电场不同热稳定性下单台机组年总发电量Ej,以及海上风电场年总发电量E。
2.根据权利要求1所述的一种海上风电场发电量计算方法,其特征在于:所述步骤2)中的海面10m高度的海面粗糙度计算公式为,
其中,z0为海面粗糙度,s为冯卡门常数、取值为0.35,U10为海面10m高度处风速。
3.根据权利要求1所述的一种海上风电场发电量计算方法,其特征在于:所述步骤3)中的理查森数Ri的计算公式为,
ΔT=T2-T1 (4)
Δu=u2-u1 (5)
其中,g为重力加速度,I为z1和z2的算术平方根,z1和z2分别为上、下两个气层的高度,T1和T2分别为上、下两个气层的温度,T为上、下两个气层的温度T1和T2的平均值,ΔT为上、下两个气层的温差值,u1和u2分别为上、下两个气层的速度,Δu为上、下两个气层的速度差;
所述步骤3)中的判定海上风电场热稳定性为,
当理查森数Ri的数值范围为Ri>0.2,则判定海上风电场热稳定性为稳定;
当理查森数Ri的数值范围为-0.6<Ri<=0.2,则判定海上风电场热稳定性为中性;
当理查森数Ri的数值范围为-2.5<Ri<=-0.6,则判定海上风电场热稳定性为不稳定;
当理查森数Ri的数值范围为Ri<=-2.5,则判定海上风电场热稳定性为极不稳定。
4.根据权利要求3所述的一种海上风电场发电量计算方法,其特征在于:所述步骤4)中的莫宁长度L的计算公式为,
其中,L为莫宁长度,I为z1和z2的算术平方根,Ri为理查森数。
5.根据权利要求4所述的一种海上风电场发电量计算方法,其特征在于:所述步骤5)中的风轮廓模型修正为,
其中,u*为摩擦速度,s为冯卡门常数,取值为0.35;z为垂直方向上风轮廓对应的高度值,z0为海面粗糙度,ψm为风速的通用函数,L为莫宁长度;
摩擦速度u*的求解公式为u* 2=C10U10 2,其中C10为阻力系数,C10可以根据吴京公式求解,U10为海面10m高度处风速;
其中,风速的通用函数ψm的计算公式为,
当海上风电场热稳定性为稳定时,
当海上风电场热稳定性为中性时,
当海上风电场热稳定性为不稳定或极不稳定时,
y=[1-(16·z/L)]1/4 (11)
其中,y为公式(10)中单列出的公共项;
所述步骤5)中的求解修正后的风速数据V0为,根据式(7)及风速的通用函数ψm的计算公式,计算考虑了海上风电场热稳定性和海面粗糙度后的轮毂高度z=H处的风速数据u(z)。
6.根据权利要求5所述的一种海上风电场发电量计算方法,其特征在于:所述步骤6)对Jensen尾流模型进行修正,对海上风电场不同热稳定性赋予对应的不同的尾流耗散系数k中的尾流耗散常数kw,具体为,
6-1)根据步骤3)判定的海上风电场热稳定性,将步骤5)求解得到的修正后的风速数据V0对应分为稳定、中性、不稳定和极不稳定共四组;
6-2)计算尾流耗散系数k,k=kw(σG+σ0)/v0;
其中,kw为尾流耗散常数,σG和σ0分别为风电机组产生的湍流和自然湍流的均方差,v0为自然风速;
6-3)对海上风电场不同热稳定性赋予对应的不同的尾流耗散系数k中的尾流耗散常数kw,包括,
当海上风电场热稳定性为稳定时,赋予尾流耗散系数k中的尾流耗散常数kw值为0.098;
当海上风电场热稳定性为中性时,赋予尾流耗散系数k中的尾流耗散常数kw值为0.048;
当海上风电场热稳定性为不稳定时,赋予尾流耗散系数k中的尾流耗散常数kw值为0.051;
当海上风电场热稳定性为极不稳定时,赋予尾流耗散系数k中的尾流耗散常数kw值为0.044。
7.根据权利要求6所述的一种海上风电场发电量计算方法,其特征在于:所述步骤7)求解单台机组尾流影响下的下游机组风轮前的风速vx,具体为,
7-1)根据动量理论有,
其中,ρ为空气密度,R和Rw分别为叶轮半径和尾流半径,vx为受尾流影响的风速,vT为通过叶片的风速;
7-2)根据推力系数公式求解得到自然风速v0、通过叶片的风速vT和风电机组的推力系数CT具有如下关系,
vT=v0(1-CT)1/2 (15)
7-3)单台机组尾流影响下的下游机组风轮前的风速vx的计算公式为,
其中,X为两个风电机组的间距。
8.根据权利要求7所述的一种海上风电场发电量计算方法,其特征在于:所述步骤8)求解多台机组尾流影响下的下游机组风轮前的风速vj(t),其计算公式为,
其中,vj(t)为作用在任意一台机组上的风速,vj0(t)为没有经过任何塔影影响作用在第j台风力发电机组上的风速、即自由流风速,vmj(t)为考虑机组间尾流效应时第m台风力发电机组作用在第j台风力发电机组上的尾流风速,表示在第j台风力发电机组处第m台风力发电机组的投影面积与第j台风力发电机组面积Arot-j的比,n为风力发电机组总台数,t表示时刻。
9.根据权利要求8所述的一种海上风电场发电量计算方法,其特征在于:所述步骤10)分别计算海上风电场不同热稳定性下单台机组年总发电量Ej,以及海上风电场年总发电量E,具体为,
10-1)采用威布尔分布描述平均风速变化,
则平均风速的概率密度函数f(v)为,
平均风速的累积分布函数F(v)为,
F(v)=1-exp(-(v/c)p) (19)
其中,p为决定分布范围的形状参数,c为决定位置的尺度参数,v为测风数据中的风速;
10-2)将步骤5)修正后的风速数据V0带入公式(18)中的v,分别求解海上风电场不同热稳定性下的平均风速的概率密度函数;
将步骤5)修正后的风速数据V0带入公式(19)中的v,分别求解海上风电场不同热稳定性下的平均风速的累积分布函数;
10-3)计算单台机组年总发电量Ej,
其中,p(θ)为对应角度θ的风向频率,vin为风电机组的切入风速,vout为风电机组的切出风速,N(v)为相应风速等级出现的全年累计小时数,E(v)为速度v通过风电场风力发电机组的功率曲线得到的单台机组输出功率,f(v)为通过威布尔分布的平均风速的概率密度函数得到的频率;
将步骤5)修正后的风速数据V0带入公式(20)中的v,分别求解海上风电场不同热稳定性下的单台机组年总发电量Ej;
10-4)计算海上风电场年总发电量E,
其中,n为海上风电场的风电机组总台数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610995568.0A CN106548414B (zh) | 2016-11-11 | 2016-11-11 | 一种海上风电场发电量计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610995568.0A CN106548414B (zh) | 2016-11-11 | 2016-11-11 | 一种海上风电场发电量计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106548414A true CN106548414A (zh) | 2017-03-29 |
CN106548414B CN106548414B (zh) | 2020-03-10 |
Family
ID=58395710
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610995568.0A Active CN106548414B (zh) | 2016-11-11 | 2016-11-11 | 一种海上风电场发电量计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106548414B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107944153A (zh) * | 2017-11-28 | 2018-04-20 | 华北电力大学 | 风轮、风力发电机组、发电系统及风轮半径的确定方法 |
CN108537411A (zh) * | 2018-03-16 | 2018-09-14 | 中国能源建设集团江苏省电力设计院有限公司 | 一种低风速风电场风电机组粗糙度及风资源分析计算方法 |
CN108629461A (zh) * | 2018-05-14 | 2018-10-09 | 华北电力大学 | 一种基于简化动量定理的近场尾流预测模型 |
CN108959174A (zh) * | 2018-07-27 | 2018-12-07 | 中国大唐集团新能源科学技术研究院有限公司 | 一种风电系统发电量的计算方法 |
CN109190212A (zh) * | 2018-08-20 | 2019-01-11 | 明阳智慧能源集团股份公司 | 复杂地形陆上风电场智能化自动布机方法及其实现系统 |
CN110543649A (zh) * | 2018-05-29 | 2019-12-06 | 北京金风科创风电设备有限公司 | 基于快速评估流体模型和尾流模型的风机排布方法和装置 |
CN111985097A (zh) * | 2020-08-13 | 2020-11-24 | 中国大唐集团未来能源科技创新中心有限公司 | 一种考虑海浪波高影响的海上风电机组尾流计算方法 |
CN112241612A (zh) * | 2020-09-15 | 2021-01-19 | 浙江运达风电股份有限公司 | 一种考虑大气热稳定性的风资源评估综合计算外推方法 |
CN113204731A (zh) * | 2021-04-28 | 2021-08-03 | 华北电力大学 | 一种基于中性等效风速的海上风廓线模型及其使用方法 |
CN117849797A (zh) * | 2024-01-29 | 2024-04-09 | 中国船舶集团风电发展有限公司 | 根据合成孔径雷达反演海表风速方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104077435A (zh) * | 2014-06-20 | 2014-10-01 | 内蒙古电力勘测设计院有限责任公司 | 用于风电场设计与优化的方法和系统 |
CN105825438A (zh) * | 2016-04-06 | 2016-08-03 | 中国能源建设集团江苏省电力设计院有限公司 | 一种基于海面粗糙度的海上风电场风能评估方法 |
-
2016
- 2016-11-11 CN CN201610995568.0A patent/CN106548414B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104077435A (zh) * | 2014-06-20 | 2014-10-01 | 内蒙古电力勘测设计院有限责任公司 | 用于风电场设计与优化的方法和系统 |
CN105825438A (zh) * | 2016-04-06 | 2016-08-03 | 中国能源建设集团江苏省电力设计院有限公司 | 一种基于海面粗糙度的海上风电场风能评估方法 |
Non-Patent Citations (1)
Title |
---|
ALFREDO PENA 等: "Modeling large offshore wind farms under different atmospheric stability regimes with the Park wake model", 《RENEWABLE ENERGY》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107944153A (zh) * | 2017-11-28 | 2018-04-20 | 华北电力大学 | 风轮、风力发电机组、发电系统及风轮半径的确定方法 |
CN108537411B (zh) * | 2018-03-16 | 2021-08-13 | 中国能源建设集团江苏省电力设计院有限公司 | 一种低风速风电场风电机组粗糙度及风资源分析计算方法 |
CN108537411A (zh) * | 2018-03-16 | 2018-09-14 | 中国能源建设集团江苏省电力设计院有限公司 | 一种低风速风电场风电机组粗糙度及风资源分析计算方法 |
CN108629461A (zh) * | 2018-05-14 | 2018-10-09 | 华北电力大学 | 一种基于简化动量定理的近场尾流预测模型 |
CN110543649B (zh) * | 2018-05-29 | 2023-04-07 | 北京金风科创风电设备有限公司 | 基于快速评估流体模型和尾流模型的风机排布方法和装置 |
CN110543649A (zh) * | 2018-05-29 | 2019-12-06 | 北京金风科创风电设备有限公司 | 基于快速评估流体模型和尾流模型的风机排布方法和装置 |
CN108959174A (zh) * | 2018-07-27 | 2018-12-07 | 中国大唐集团新能源科学技术研究院有限公司 | 一种风电系统发电量的计算方法 |
CN109190212A (zh) * | 2018-08-20 | 2019-01-11 | 明阳智慧能源集团股份公司 | 复杂地形陆上风电场智能化自动布机方法及其实现系统 |
CN109190212B (zh) * | 2018-08-20 | 2023-05-26 | 明阳智慧能源集团股份公司 | 复杂地形陆上风电场智能化自动布机方法及其实现系统 |
CN111985097A (zh) * | 2020-08-13 | 2020-11-24 | 中国大唐集团未来能源科技创新中心有限公司 | 一种考虑海浪波高影响的海上风电机组尾流计算方法 |
CN112241612A (zh) * | 2020-09-15 | 2021-01-19 | 浙江运达风电股份有限公司 | 一种考虑大气热稳定性的风资源评估综合计算外推方法 |
CN113204731A (zh) * | 2021-04-28 | 2021-08-03 | 华北电力大学 | 一种基于中性等效风速的海上风廓线模型及其使用方法 |
CN113204731B (zh) * | 2021-04-28 | 2023-09-29 | 华北电力大学 | 一种基于中性等效风速的海上风廓线模型的使用方法 |
CN117849797A (zh) * | 2024-01-29 | 2024-04-09 | 中国船舶集团风电发展有限公司 | 根据合成孔径雷达反演海表风速方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106548414B (zh) | 2020-03-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106548414A (zh) | 一种海上风电场发电量计算方法 | |
Hezaveh et al. | Increasing the power production of vertical-axis wind-turbine farms using synergistic clustering | |
Wu et al. | Modeling turbine wakes and power losses within a wind farm using LES: An application to the Horns Rev offshore wind farm | |
CN102663251B (zh) | 基于计算流体力学模型的风电场功率物理预测方法 | |
Han et al. | Atmospheric stability and topography effects on wind turbine performance and wake properties in complex terrain | |
Jeong et al. | Blockage corrections for wind tunnel tests conducted on a Darrieus wind turbine | |
Frandsen | Turbulence and turbulence-generated structural loading in wind turbine clusters | |
Saleem et al. | Aerodynamic analysis of an airborne wind turbine with three different aerofoil-based buoyant shells using steady RANS simulations | |
Vermeer et al. | Wind turbine wake aerodynamics | |
Battisti et al. | Aerodynamic measurements on a vertical axis wind turbine in a large scale wind tunnel | |
Krogstad et al. | “Blind Test 3” calculations of the performance and wake development behind two in-line and offset model wind turbines | |
CN103996074B (zh) | 基于cfd和改进pso的复杂地形风电场微观选址方法 | |
Lee et al. | Observing and simulating wind-turbine wakes during the evening transition | |
CN103745024B (zh) | 基于三维尾流模型修正风电机组尾部风速功率特性评估法 | |
CN107194097A (zh) | 基于风电场气动模拟和风速风向数据的分析方法 | |
CN101949363A (zh) | 计及风电场输入风速和风向随机波动的风电机组分组方法 | |
CN103617308B (zh) | 一种风电场频域等效模型的构建方法 | |
Gupta | Development of a time-accurate viscous Lagrangian vortex wake model for wind turbine applications | |
Pierella et al. | Wind tunnel investigation on the effect of the turbine tower on wind turbines wake symmetry | |
van der Laan et al. | Improved energy production of multi-rotor wind farms | |
Krüger et al. | Validation of a coupled atmospheric-aeroelastic model system for wind turbine power and load calculations | |
Zhang et al. | Effects of blade number on the aerodynamic performance and wake characteristics of a small horizontal-axis wind turbine | |
Wu et al. | Investigation of multiblade wind-turbine wakes in turbulent boundary layer | |
Taschner et al. | On the performance of the helix wind farm control approach in the conventionally neutral atmospheric boundary layer | |
Abedi | Assessment of flow characteristics over complex terrain covered by the heterogeneous forest at slightly varying mean flow directions:(A case study of a Swedish wind farm) |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |