CN104794259B - 基于测风塔相互验证的风电场上网电量偏差计算方法 - Google Patents
基于测风塔相互验证的风电场上网电量偏差计算方法 Download PDFInfo
- Publication number
- CN104794259B CN104794259B CN201510107717.0A CN201510107717A CN104794259B CN 104794259 B CN104794259 B CN 104794259B CN 201510107717 A CN201510107717 A CN 201510107717A CN 104794259 B CN104794259 B CN 104794259B
- Authority
- CN
- China
- Prior art keywords
- data
- anemometer tower
- power plant
- wind power
- 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.)
- Active
Links
Landscapes
- Wind Motors (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于测风塔相互验证的风电场上网电量偏差计算方法,其中步骤一:勘测风电场场区地形资料;步骤二:勘测风电场场区粗糙度;步骤三:测风塔完整有效性分析;步骤四:获取各点实测数据;步骤五:利用步骤一、步骤二得到的风电场场区地形资料数据CON、风电场场区粗糙度数据ROU及其它辅助资料(如机组参数等),采用Meteodyn WT进行风资源模拟分析;分别采用其中一个测风塔模拟推算其它测风塔风速;步骤六:通过实测数据与推测数据计算计算偏差值;步骤七:计算测风塔相互验证加权平均误差值,其中权值为各塔的有效完整率,得到整个场区的综合误差,再通过总误差ER计算上网电量误差。
Description
技术领域
本发明涉及风电场上网电量折减分析,更具体地说它是一种基于测风塔相互验证的风电场上网电量偏差计算方法。
背景技术
风电场年上网电量的准确估算多少直接关系到风电场的效益水平和风险程度,决定着企业投资方向和投资决策。
风电场年上网电量通常依照《国家发展和改革委员会关于对中国风电发电量折减问题的说明》的方法来计算分析:风电场发电量折减指对影响风电场实际出力的因素进行逐个分析,得出各因素所引起的发电量减少的数值。折减因素包括但不限于风电机组尾流的影响、空气密度、低温等气候条件、控制和湍流、风电机组可利用率、风电场功率曲线保证值、场用电、线损等。
(1)尾流折减
尾流折减是指风电机组间由于相互影响而降低的发电量,反映了风电场风电机组的排列效率,通常采用相关的尾流模型进行分析计算。
(2)空气密度折减
空气密度折减指用于发电量计算的风电机组功率曲线(平均风速/功率曲线)对应的空气密度与场区实际空气密度不一致时,风电机组受此影响而降低的发电量,视具体情况计算而定。
(3)控制和湍流折减
控制和湍流折减指由于风电机组受风电场内湍流等风况影响而 降低的发电量。其典型折减数值可在5%左右。
(4)叶片污染折减
叶片污染折减是指由于叶片表面污染,影响了其对风能的捕获能力而降低的发电量,其典型折减数值可在6%以内。
(5)风电机组可利用率折减
风电机组可利用率折减是指由于风电机组本身质量的问题,降低了正常工作时间而减少的发电量。由于国内风电机组整机制造厂商设计运行经验不足,机组技术成熟度较低,制造工艺相对落后,售后服务跟不上等原因,风电机组可利用率常达不到设计值,其典型的折减数值为5%~10%。
(6)风电机组功率曲线保证率折减
风电机组功率曲线保证率折减是指风电机组实际功率曲线达不到设计值,从而降低的发电量。风电机组厂商提供的功率曲线保证率一般为95%,此项折减的典型值为5%。
(7)场用电、线损等折减
场用电、线损等折减指风电机组发出的电力,在集电线路中的损耗和场内自用的损耗,其数值可以根据风电场设计方案估算得出,典型值一般在3%~10%之间,根据风电场实际情况而定。
(8)软件计算误差折减
软件计算误差折减是指由于发电量计算软件对风电场适应性差异,可能造成理论发电量估算的不确定性,需要进行折减,其典型折减数值一般在5%~10%之间。
(9)电网频率波动与限电等折减
电网频率波动与限电等折减指由于电网的频率波动、为保障电网安全而暂时限电等影响,风电场会由于暂时脱网而影响风电场发电量上网,其典型数值一般在3%~5%之间。
(10)大规模风电场尾流折减
大规模风电场尾流折减指由于大规模开发风电场,影响了小区域 内的风况特征,使得风能无法及时得到恢复,从而降低了风电场的整体出力,对于此种影响的分析方法和折减的具体数值,现阶段国内外仍处于探索中。
由于中国地域辽阔,各地区自然条件差异较大,以上折减系数需要根据风电场实际情况进行适当调整。
目前国内风电场的年上网电量采用的综合折减法计算公式如下所示:
AEPnet=AEPgross×Nt
Nt=∑Ni
式中:AEPnet为年上网电量,AEPgross为年理论发电量,Nt为综合折减系数,Ni为第i个折减因素的折减系数。
利用风资源评估软件计算风电场扣除尾流后的年理论发电,再考虑风电场各种折减因素的影响程度最终确定。风电场年上网电量计算折减因素众多,如尾流、空气密度、叶片污染等。这些影响因素大部分可以依照经验值进行计算分析,而风资源评估软件的计算误差无法通过经验值来合理估算,若风资源评估软件分析误差仅采用常规经验值来分析,带来的误差将影响风电场年上网电量估算的准确度,直接影响着企业对本风电开发项目的投资决策。
复杂地形条件下的风电场,由于风资源评估软件无法准确模拟出场区风能资源分布,由此引起的评估误差更大,造成项目投资的不确定性更大。如何合理确定风资源评估软件引起的误差,准确计算风电场的年上网电量成为风电场投资开发的关键。目前关于风资源评估软件带来的年上网电量偏差没有实际有效的解决办法。
发明内容
本发明的目的在于提供一种克服现有的风电场的年上网电量综合折减系数中的软件分析误差无法合理确定的基于测风塔相互验证的风电场上网电量偏差计算方法。
本发明是按如下的技术方案实现的:
步骤一:勘测风电场场区地形资料
依据风电场区域地形资料,结合现场各机位实际踏勘,确定各机位地形的情况,再对场区局部复杂地形进行核实,最终调整地形资料的偏差,得到整个风电场场区地形资料数据CON;
步骤二:勘测风电场场区粗糙度
依据风电场场区地理影像资料,分析各机位和场区周边植被、建筑等情况,调整风电场场区地表粗糙因子情况,得到整个风电场场区粗糙度数据ROU;
步骤三:测风塔完整有效性分析
根据测风塔实测数据,结合现场踏勘,分析判断各测风塔的完整有效率;测风塔Ti(其中i=1:n)的有效完整率为Vi(其中i=1:n);测风塔有效完整率的计算方式为:
Vi=(YCNi-QCNi-UVNi)/YCNi (1)
其中YCNi(其中i=1:n)为Ti数据的应测数量;QCNi(其中i=1:n)为Ti数据的为缺测数量;UCNi(其中i=1:n)为Ti数据的为无效数据数量;
步骤四:获取各点实测数据,测风塔Ti的实测风速为Mi(其中i=1:n);
步骤五:利用步骤一、步骤二得到的风电场场区地形资料数据CON、风电场场区粗糙度数据ROU及其它辅助数据(如机组参数等),采用Meteodyn WT进行风资源模拟分析;分别采用其中一个测风塔模拟推算其它测风塔风速,其中用Ti的风速Mi推测测风塔Tj的推测风速S(i,j)(其中i,j=n)的计算方式为:
S(i,j)=Meteodyn WT[Mi,j,CON,ROU,辅助数据] (2)
步骤六:通过实测数据与推测数据计算计算偏差值,用Ti的风速Mi推测测风塔Tj风速的偏差Ems(i,j)(其中i,j=n)的计算方式为:
Ems(i,j)=|(Mj-S(i,j))|/Mi*100% (3)
步骤七:计算测风塔相互验证加权平均误差值:
计算得到各测风塔间的相互模拟误差后,取各测风塔推算偏差的加权平均误差值,其中权值为各塔的有效完整率V作为对整个场区模拟误差,用以调整风电场年上网电量的计算值;
测风塔相互验证加权平均误差值的详细计算过程如下所示,其中用Ti的实测风速Mi推测数据的误差Eri(i=1:n)为各分项误差;ER为测风塔相互验证加权平均得到的综合误差;
其中Eri的计算方式如下:
Eri=(Ems(i,1)+Ems(i,2)+…Ems(i,i-1)+Ems(i,i+1)…Ems(i,n))/n-1,i=1:n (4)
整个场区的综合误差的计算方式如下:
ER=(V1*Er1+V2*Er2+…+Vn*Ern)/n,i=1:n (5)
再通过总误差ER计算上网电量误差。
附图说明
图1为T1测风塔推算的场区风速分布情况。
图2为T2测风塔推算的场区风速分布情况。
图3为T3测风塔推算的场区风速分布情况。
图4为运用软件推测其他测风塔风速流程图。
图5为T1测风塔模拟其它测风塔的结果。
图6为T2测风塔模拟其它测风塔的结果。
图7为T3测风塔模拟其它测风塔的结果。
图8为各测风塔模拟的分项误差值
具体实施方式
以下结合附图和实施例对本发明作进一步地详细描述,但该实施例不应该理解为对本发明的限制,仅作举例而已,同时通过说明本发明的优点将变得更加清楚和容易理解。
基于测风塔相互验证的风电场上网电量偏差计算方法,以某风电场为例,该风电场设计安装24台2.0MW风力发电机组,总装机容量48MW,机组轮毂高度80m。风电场的设计年上网电量为8891万kWh,风电场年等效负荷小时数为1852h,容量系数为0.211。
步骤一:风电场场区地形资料(CON)
采用风电场场址内1:10000地形图,现场实测1:2000地形图。各机位详细踏勘,对地形图未能准确表达场区实际地形的区域进行纠偏处理,以保证地形图能基本反映风电场实际情况。
步骤二:风电场场区粗糙度(ROU)
采用风电场场址卫星影像图制作场区的粗糙度文件,由于影像图无法准确表达场区内植被和建筑物的高度,现场实际踏勘确定植被和建筑物高度,修正场区粗糙度文件。
步骤三:测风塔完整有效性分析
整个风电场场区内设立有三个测风塔,分别位于风电场的北部、中部和南部区域,塔编号分别为T1、T2和T3,测风塔风向和风能分布存在一定差异。根据测风塔数据进行分析,判断各测风塔的完整有效性Vi(其中i=1,2,3)。
Vi=(YCNi-QCNi-UVNi)/YCNi (1)
其中YCNi(其中i=1,2,3)为Ti数据的应测数量;QCNi(其中i=1,2,3)为Ti数据的为缺测数量;UCNi(其中i=1,2,3)为Ti数据的为无效数据数量;
测风塔有效完整率表
测风塔 | 有效完整率 |
T1 | 0.971 |
T2 | 1.000 |
T3 | 0.908 |
步骤四:获取各点实测数据,测风塔Ti的实测风速为Mi(其中i=1,2,3);
实测测得M1=6.48,M2=6.52,M3=6.29;
步骤五:利用步骤一、步骤二得到的风电场场区地形资料数据CON、风电场场区粗糙度数据ROU及其它辅助数据(如机组参数等),采用Meteodyn WT进行风资源模拟分析;分别采用其中一个测风塔模拟推算其它测风塔风速,其中用Ti的风速Mi推测测风塔Tj的推测风速S(i,j)(其中i,j=n)的计算方式为:
S(i,j)=Meteodyn WT[Mi,j,CON,ROU,辅助数据] (2)
测风塔相互模拟计算结果如下表:
各测风塔相互模拟的平均风速表
步骤六:通过实测数据与推测数据计算计算偏差值,用Ti的风速Mi推测测风塔Tj风速的偏差Ems(i,j)(其中i,j=n)的计算方式为:
Ems(i,j)=|(Mj-S(i,j))|/Mi*100% (3)
各测风塔间相互模拟计算结果偏差值分别为:
Ems(1,2)=|6.52-6.840|/6.52=4.908%
Ems(1,3)=|6.29-6.510|/6.29=3.497%
Ems(2,1)=|6.48-6.579|/6.48=1.533%
Ems(2,3)=|6.29-6.348|/6.29=0.926%
Ems(3,1)=|6.48-6.630|/6.48=2.314%
Ems(3,2)=|6.52-6.590|/6.52=1.073%
表4 各测风塔相互模拟的误差表
步骤七:计算测风塔相互验证加权平均误差值:
计算得到各测风塔间的相互模拟误差后,取各测风塔推算偏差的加权平均误差值,其中权值为各塔的有效完整率V作为对整个场区模拟误差,用以调整风电场年上网电量的计算值;
测风塔相互验证加权平均误差值的详细计算过程如下所示,其中用Ti的实测风速Mi推测数据的误差Eri(i=1:n)为各分项误差;ER为测风塔相互验证加权平均得到的综合误差;
其中Eri的计算方式如下:
Eri=(Ems(i,1)+Ems(i,2)+…Ems(i,i-1)+Ems(i,i+1)…Ems(i,n))/n-1,i=1:n (4)
整个场区的综合误差的计算方式如下:
ER=(V1*Er1+V2*Er2+…+Vn*Ern)/n,i=1:n (5)
再通过总误差ER计算上网电量误差。
现将该风电场数据带入计算结果如下:
Er1=Ems(1,2)+Ems(1,3)=(4.908%+3.497%)/2=4.203%
Er2=Ems(2,1)+Ems(2,3)=(1.533%+0.926%)/2=1.230%
Er3=Ems(3,1)+Ems(3,2)=(2.314%+1.073%)/2=1.230%
表5 测风塔相互验证加权平均误差计算结果表
依照前述步骤(1-4),计算测风塔相互验证加权平均误差值:
ER=(V1*Er1+V2*Er2+…+Vn*Ern)/n=(V1*Er1+V2*Er2+V3*Er3)/3
=(0.971*4.203%+1.000*1.230%+0.908*1.694%)/3=2.283%
采用基于测风塔相互验证降低复杂风电场年发电量计算偏差综合计算方法计算的测风塔相互验证加权平均误差值作为整个场区模拟误差,用以调整风电场年上网电量的设计值。通过各测风塔相互推算分析,软件误差为2.283%,用以作为年上网电量软件误差偏差值,其它因素偏差纠正值采用经验值的取值。风电场实际运行满一年后,对比设计计算结果,风电场年上网量偏差仅为1.68%(该风电场实际发电量未考虑第一年机组故障,当年电网接入问题等因素,若考虑影响发电量的情况,实际发电量值与设计值将更加一致),说明该风电场计算效果非常理想。
需要说明的是:对于所属领域的技术人员来说,在不改变本发明原理的前提下还可以对本发明作出若干的改变或变形,这同样属于本发明的保护范围。
Claims (1)
1.基于测风塔相互验证的风电场上网电量偏差计算方法,其特征在于,包括如下步骤:
步骤一:勘测风电场场区地形资料
依据风电场区域地形资料,结合现场各机位实际踏勘,确定各机位地形的情况,再对场区局部复杂地形进行核实,最终调整地形资料的偏差,得到整个风电场场区地形资料数据CON;
步骤二:勘测风电场场区粗糙度
依据风电场场区地理影像资料,分析各机位和场区周边植被、建筑情况,调整风电场场区地表粗糙因子情况,得到整个风电场场区粗糙度数据ROU;
步骤三:测风塔完整有效性分析
根据测风塔实测数据,结合现场踏勘,分析判断各测风塔的完整有效率;测风塔Ti,其中i=1:n的有效完整率为Vi,其中i=1:n;测风塔有效完整率的计算方式为:
Vi=(YCNi-QCNi-UVNi)/YCNi (1)
其中YCNi,其中i=1:n,为Ti数据的应测数量;QCNi,其中i=1:n,为Ti数据的缺测数量;UVNi,其中i=1:n,为Ti数据的无效数据数量;
步骤四:获取各点实测数据,测风塔Ti的实测风速为Mi,其中i=1:n;
步骤五:利用步骤一、步骤二得到的风电场场区地形资料数据CON、风电场场区粗糙度数据ROU及辅助数据,采用MeteodynWT进行风资源模拟分析;分别采用其中一个测风塔模拟推算其它测风塔风速,其中用Ti的风速Mi推测测风塔Tj的推测风速S(i,j),其中i,j=n,的计算方式为:
S(i,j)=MeteodynWT[Mi,j,CON,ROU,辅助数据] (2);
步骤六:通过实测数据与推测数据计算计算偏差值,用Ti的风速 Mi推测测风塔Tj风速的偏差Ems(i,j),其中i,j=n,的计算方式为:
Ems(i,j)=|(Mj-S(i,j))|/Mi*100% (3)
步骤七:计算测风塔相互验证加权平均误差值:
计算得到各测风塔间的相互模拟误差后,取各测风塔推算偏差的加权平均误差值,其中权值为各塔的有效完整率V作为对整个场区模拟误差,用以调整风电场年上网电量的计算值;
测风塔相互验证加权平均误差值的详细计算过程如下所示,其中用Ti的实测风速Mi推测数据的误差Eri,i=1:n,为各分项误差;ER为测风塔相互验证加权平均得到的综合误差;
其中Eri的计算方式如下:
Eri=(Ems(i,1)+Ems(i,2)+…Ems(i,i-1)+Ems(i,i+1)…Ems(i,n))/n-1,i=1:n (4)
整个场区的综合误差的计算方式如下:
ER=(V1*Er1+V2*Er2+…+Vn*Ern)/n,i=1:n (5)
再通过总误差ER计算上网电量误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510107717.0A CN104794259B (zh) | 2015-03-12 | 2015-03-12 | 基于测风塔相互验证的风电场上网电量偏差计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510107717.0A CN104794259B (zh) | 2015-03-12 | 2015-03-12 | 基于测风塔相互验证的风电场上网电量偏差计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104794259A CN104794259A (zh) | 2015-07-22 |
CN104794259B true CN104794259B (zh) | 2017-10-10 |
Family
ID=53559050
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510107717.0A Active CN104794259B (zh) | 2015-03-12 | 2015-03-12 | 基于测风塔相互验证的风电场上网电量偏差计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104794259B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106296462B (zh) * | 2016-08-18 | 2019-05-31 | 中国能源建设集团江苏省电力设计院有限公司 | 一种基于双测风塔数据的既有风电场粗糙度值确定方法 |
CN108536881B (zh) | 2017-03-06 | 2020-10-16 | 新疆金风科技股份有限公司 | 用于计算风电场发电量的方法和设备 |
CN107194141B (zh) * | 2017-03-24 | 2020-04-24 | 中国农业大学 | 一种区域风能资源精细化评估方法 |
CN107153894B (zh) | 2017-06-02 | 2018-11-27 | 北京金风科创风电设备有限公司 | 一种风电场的预测风速校正方法及装置 |
CN107958335A (zh) * | 2017-12-06 | 2018-04-24 | 北京天润新能投资有限公司 | 一种基于开发地图的风电场装机规模及上网电量估算系统及估算方法 |
CN110942218B (zh) * | 2018-09-25 | 2023-10-10 | 北京天润新能投资有限公司 | 计算风电场上网电量的方法、装置以及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103034763A (zh) * | 2012-12-12 | 2013-04-10 | 上海市电力公司 | 大容量风电场风机聚合模型建立和仿真方法 |
CN103219725A (zh) * | 2013-01-18 | 2013-07-24 | 云南电网公司大理供电局 | 一种基于实时运行数据的风电场等值建模方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6948910B2 (en) * | 2002-07-12 | 2005-09-27 | Polacsek Ronald R | Spiral-based axial flow devices |
US20100217550A1 (en) * | 2009-02-26 | 2010-08-26 | Jason Crabtree | System and method for electric grid utilization and optimization |
US8873813B2 (en) * | 2012-09-17 | 2014-10-28 | Z Advanced Computing, Inc. | Application of Z-webs and Z-factors to analytics, search engine, learning, recognition, natural language, and other utilities |
-
2015
- 2015-03-12 CN CN201510107717.0A patent/CN104794259B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103034763A (zh) * | 2012-12-12 | 2013-04-10 | 上海市电力公司 | 大容量风电场风机聚合模型建立和仿真方法 |
CN103219725A (zh) * | 2013-01-18 | 2013-07-24 | 云南电网公司大理供电局 | 一种基于实时运行数据的风电场等值建模方法 |
Non-Patent Citations (3)
Title |
---|
基于机舱风速数据的;姜文玲等;《电网技术》;20140331;第38卷(第3期);第647-652页 * |
输电线路对风电场无功电能计量偏差的影响分析;李勇等;《内蒙古电力技术》;20121231;第30卷(第6期);第122-124页 * |
风电场的发电可靠性模型及其应用;陈树勇等;《中国电机工程学报》;20000331;第20卷(第3期);第26-29页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104794259A (zh) | 2015-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104794259B (zh) | 基于测风塔相互验证的风电场上网电量偏差计算方法 | |
US11408399B2 (en) | Forecasting output power of wind turbine in wind farm | |
CN102628876B (zh) | 一种包含上下游效应实时监测的超短期预测方法 | |
CN102663251B (zh) | 基于计算流体力学模型的风电场功率物理预测方法 | |
CN102567636B (zh) | 基于相关向量机的高精度风电场功率区间预测方法 | |
CN104319807B (zh) | 一种基于Copula函数获取多风电场容量可信度的方法 | |
Abbes et al. | Wind resource estimation and wind park design in El-Kef region, Tunisia | |
CN102945508A (zh) | 一种基于模型校正的风电功率预测预报系统及方法 | |
Che et al. | A wind power forecasting system based on the weather research and forecasting model and Kalman filtering over a wind-farm in Japan | |
CN109784563B (zh) | 一种基于虚拟测风塔技术的超短期功率预测方法 | |
CN106321368A (zh) | 基于运营风机scada数据的海上风电场尾流损失测量方法 | |
CN113496311A (zh) | 光伏电站发电功率预测方法及系统 | |
CN106815773A (zh) | 一种风电场功率特性评估方法 | |
CN107679361A (zh) | 一种测风数据处理方法 | |
CN107313901A (zh) | 一种风电机组发电量损失的评估方法及装置 | |
CN115358606B (zh) | 一种平缓地形下在役风电场能效评估方法及系统 | |
Srinivas et al. | Multi-site downscaling of maximum and minimum daily temperature using support vector machine | |
CN102664409A (zh) | 一种基于实测数据风电场风电功率实时预测的计算方法 | |
Bo et al. | Hybrid PSO-BP neural network approach for wind power forecasting | |
Louassa et al. | Effects of local ambient air temperatures on wind park performance: case of the Kaberten wind park | |
CN105262148B (zh) | 计及风电出力特性的规划年电力平衡方法 | |
Ma et al. | Ensemble transform with 3D rescaling initialization method | |
Rahmani et al. | Wind power assessment and site matching of wind turbines in Lootak of Zabol | |
CN116993030A (zh) | 变化条件下水库压咸补淡调度方法及系统 | |
CN115586591A (zh) | 风暴潮的自动预报方法以及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |