CN106547996B - 一种风电场虚拟测风塔构建方法 - Google Patents
一种风电场虚拟测风塔构建方法 Download PDFInfo
- Publication number
- CN106547996B CN106547996B CN201611119467.3A CN201611119467A CN106547996B CN 106547996 B CN106547996 B CN 106547996B CN 201611119467 A CN201611119467 A CN 201611119467A CN 106547996 B CN106547996 B CN 106547996B
- Authority
- CN
- China
- Prior art keywords
- wind power
- power plant
- wind
- data
- moment
- 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
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
本发明公开了一种风电场虚拟测风塔构建方法,它包括:1采集风电场n台风机的基础数据和历史数据;2采集风电场n台风机实时数据;3根据风电场n台风机的历史数据和实时数据信息,模拟计算第ξ个时刻整个风电场轮毂高度的风速和风向;4计算第ξ个时刻整个风电场轮毂高度的温度Tsrtξ;5、计算第ξ个时刻风电场空气密度、湿度和气压;6根据计算得到的数据,形成风电场第ξ时刻的虚拟测风塔数据;解决了由于高原山区的特殊地理位置和气象特征下,几个测风塔或风机测风仪数据无法全面代表复杂地形的整个风电场,以及覆冰等灾害频发造成测风塔倒塌、测风仪器损坏或短时间工作失灵造成测风塔数据质量差风电场的风电功率预测精度较低等问题。
Description
技术领域
本发明属于风力发电技术领域,尤其涉及一种风电场虚拟测风塔构建方法。
背景技术
高原山区风电场的风机排布,相比平原风电场要复杂得多。高原山区风电场的风机排布,没有统一的规律,其排布不但要满足平原风电场的间距和排距原则,还需要根据地形地貌的具体情况,尽量排布在山脊区域。于是,风机之间的距离没有规律,一般都比平原风电场的风机间距大,高原山区的风电场所占面积也大。此时,风机尾流的影响往往不是影响风电场出力的重要的因素,而风电场内风速风向的微观变化,对风电场出力影响很大。
由于高原山区风电场每一个风电机组所处的位置都有一定的特殊性,其湍流与回风的影响都不一样,导致高原山区同一测风塔不同高度间风速相关性好,而同一风电场各测风塔间的风速相关性不好,使得测风塔不具有足够的代表性,任何一个或几个测风塔或风机测风仪数据都无法全面代表复杂地形的整个山地风电场,而且由于高原山区的特殊地理位置和气象特征使得覆冰灾害频发,已立的测风塔因覆冰或年久失修而损坏或倒塔,或者覆冰造成测风仪器损坏或短时间工作失灵,测风数据明显异常、不可信。由于以上问题极易造成测风塔数据质量差导致风电场的风电功率预测精度较低。
发明内容:
本发明要解决的技术问题:提供一种风电场虚拟测风塔构建方法,以解决现有技术高原山区风电场的高原山区的特殊地理位置和气象特征下,任何一个或几个测风塔或风机测风仪数据都无法全面代表复杂地形的整个风电场,以及覆冰等灾害频发造成测风塔倒塌、测风仪器损坏或短时间工作失灵,造成测风塔数据质量差导致风电场的风电功率预测精度较低等技术问题。
本发明技术方案:
一种风电场虚拟测风塔构建方法,它包括:
步骤1、采集风电场n台风机的基础数据和历史数据;
步骤2、采集风电场n台风机实时数据,实时数据包括风电场所在区域内所有风机的风机测风仪测得风速数据Wssrt和风向数据Wdsrt、所有风机的风机温度传感器测得温度数据Tsrt、所有风机的实时出力数据Pwtrt、整个风电场的实时出力数据Pwfrt;
步骤3、根据风电场n台风机的历史数据和实时数据信息,模拟计算第ξ个时刻整个风电场轮毂高度的风速Wssrtξ和风向Wdsrtξ;
步骤4、根据风电场n台风机的历史数据和实时数据信息,模拟计算第ξ个时刻整个风电场轮毂高度的温度Tsrtξ;
步骤5、根据步骤3、4的计算结果,计算第ξ个时刻风电场空气密度、湿度和气压;
步骤6、汇总步骤3、4、5计算得到的数据,形成风电场第ξ时刻的虚拟测风塔数据。
步骤1所述基础数据包括风电场区域内所有风机的轮毂高度、历史数据包括风电场所在区域内所有风机的历史风速数据Wssh和风向数据Wdsh、所有风机的历史温度数据Tsh、所有风机的历史出力数据Pwth和整个风电场的历史出力数据Pwfh。
步骤3所述根据风电场n台风机的历史数据和实时数据信息,模拟计算第ξ个时刻整个风电场轮毂高度的风速Wssrtξ和风向Wdsrtξ的方法包括:
步骤3.1、利用风电场所有风电机组的历史风速风向数据,以及整个风电场的历史出力数据,建立所有风电机组风速风向‐风电场风电功率数据样本库;
步骤3.2、建立所有风电机组风速风向‐风电场风电功率模型;
步骤3.3、将3.1建立的所有风电机组风速风向‐风电场风电功率数据样本库作为训练集,对风电机组风速风向‐风电场风电功率模型进行训练,确定模型中n台风电机组的风速风向数据的权值;其中,第i台风机的风速风向数据的权值记做Wqzi,则有Wqz1+Wqz2+...Wqzi+...+Wqzn=1;且有1≤i≤n;
步骤3.4、计算第ξ个时刻整个风电场轮毂高度的风速Wssrtξ和风向Wdsrtξ
第ξ个时刻整个风电场轮毂高度的风速Wssrtξ的计算公式为:
第ξ个时刻整个风电场轮毂高度的风向Wdsrtξ的表达式为:
步骤4所述的根据风电场n台风机的历史数据和实时数据信息,模拟计算第ξ个时刻整个风电场轮毂高度的温度Tsrtξ的计算方法包括:
步骤4.1、利用风电场所有风电机组的历史温度数据,以及所有风机的历史出力数据和整个风电场的历史出力数据,建立所有风电机组温度‐风电场风电功率数据样本库;
步骤4.2、建立所有风电机组温度‐风电场风电功率模型;
步骤4.3、将4.1建立的所有风电机组温度‐风电场风电功率数据样本库作为训练集,对风电机组温度‐风电场风电功率模型进行训练,确定模型中所有风电机组的温度数据的权值;
其中,第i台风机的温度数据的权值记做Tqzi,则有Tqz1+Tqz2+...Tqzi+...+Tqzn=1;
且有1≤i≤n;
步骤4.4、模拟计算第ξ个时刻整个风电场轮毂高度的温度Tsrtξ
第ξ个时刻整个风电场轮毂高度的温度Tsrtξ的表达式为:
步骤5所述根据步骤3、4的计算结果,计算第ξ个时刻风电场空气密度、湿度和气压的计算方法包括:
步骤5.1、根据风电场第ξ个时刻的输出功率Pwfrtξ与风速Wssrtξ,计算第ξ个时刻的空气密度ρsrtξ
第ξ个时刻的空气密度ρsrtξ的计算表达式为:
其中F为风轮叶片旋转一周的扫掠面积;
步骤5.2、根据大气压力和空气密度计算公式,以及空气湿度公式,计算得到湿度和气压
第ξ个时刻的气压presrtξ的计算表达式为:
presrtξ=ρsrtξ×(273.15+Tsrtξ)×R
式中,R为气体常数,其值为287;
第ξ个时刻的湿度Hsrtξ的计算表达式为:
本发明的有益效果:
针对本发明的技术问题,本发明提出了一种风电场虚拟测风塔构建方法,其主要是采集目标风电场整场的风电机组测风系统以及测温系统的观测数据,利用熵权法综合评价模型确定各台风机的权重系数来构建全场的虚拟测风塔,并给出虚拟测风塔的所有实时物理量的计算值,包括各高层的风速风向、温度湿度与气压等实时信息;解决了由于高原山区的特殊地理位置和气象特征下,任何一个或几个测风塔或风机测风仪数据都无法全面代表复杂地形的整个风电场,以及覆冰等灾害频发造成测风塔倒塌、测风仪器损坏或短时间工作失灵,造成测风塔数据质量差导致风电场的风电功率预测精度较低等问题。
具体实施方式:
一种风电场虚拟测风塔构建方法,它包括:
步骤1、采集风电场n台风机的基础数据和历史数据;步骤1所述基础数据包括风电场区域内所有风机的轮毂高度、历史数据包括风电场所在区域内所有风机的历史风速数据Wssh和风向数据Wdsh、所有风机的历史温度数据Tsh、所有风机的历史出力数据Pwth和整个风电场的历史出力数据Pwfh。
步骤2、采集风电场n台风机实时数据,实时数据包括风电场所在区域内所有风机的风机测风仪测得风速数据Wssrt和风向数据Wdsrt、所有风机的风机温度传感器测得温度数据Tsrt、所有风机的实时出力数据Pwtrt、整个风电场的实时出力数据Pwfrt;
步骤3、根据风电场n台风机的历史数据和实时数据信息,模拟计算第ξ个时刻整个风电场轮毂高度的风速Wssrtξ和风向Wdsrtξ;
步骤3.1、利用风电场所有风电机组的历史风速风向数据,以及整个风电场的历史出力数据,建立所有风电机组风速风向‐风电场风电功率数据样本库;
步骤3.2、建立所有风电机组风速风向‐风电场风电功率模型;可采用熵权法综合评价模型,建立所有风电机组风速风向‐风电场风电功率模型;
步骤3.3、将3.1建立的所有风电机组风速风向‐风电场风电功率数据样本库作为训练集,对风电机组风速风向‐风电场风电功率模型进行训练,确定模型中n台风电机组的风速风向数据的权值;其中,第i台风机的风速风向数据的权值记做Wqzi,则有Wqz1+Wqz2+...Wqzi+...+Wqzn=1;且有1≤i≤n;
步骤3.4、模拟计算第ξ个时刻整个风电场轮毂高度的风速Wssrtξ和风向Wdsrtξ
第ξ个时刻整个风电场轮毂高度的风速Wssrtξ的计算公式为:
式中:Wssrtiξ为第ξ个时刻第i台风机轮毂高度的风速。
第ξ个时刻整个风电场轮毂高度的风向Wdsrtξ的表达式为:
式中:Wdsrtiξ为第ξ个时刻第i台风机的轮毂高度的风向。
步骤4、根据风电场n台风机的历史数据和实时数据信息,模拟计算第ξ个时刻整个风电场轮毂高度的温度Tsrtξ;
步骤4所述的根据风电场n台风机的历史数据和实时数据信息,模拟计算第ξ个时刻整个风电场轮毂高度的温度Tsrtξ的计算方法包括:
步骤4.1、利用风电场所有风电机组的历史温度数据,以及所有风机的历史出力数据和整个风电场的历史出力数据,建立所有风电机组温度‐风电场风电功率数据样本库;
步骤4.2、建立所有风电机组温度‐风电场风电功率模型;本发明采用熵权法综合评价模型,建立所有风电机组温度‐风电场风电功率模型;
步骤4.3、将4.1建立的所有风电机组温度‐风电场风电功率数据样本库作为训练集,对风电机组温度‐风电场风电功率模型进行训练,确定模型中所有风电机组的温度数据的权值;
其中,第i台风机的温度数据的权值记做Tqzi,则有Tqz1+Tqz2+...Tqzi+...+Tqzn=1;
且有1≤i≤n;
步骤4.4、模拟计算第ξ个时刻整个风电场轮毂高度的温度Tsrtξ
第ξ个时刻整个风电场轮毂高度的温度Tsrtξ的表达式为:
式中:Tsrtiξ为第ξ个时刻第i台风机轮毂高度的温度。
步骤5、根据步骤3、4的计算结果,计算第ξ个时刻风电场空气密度、湿度和气压;
步骤5所述根据步骤3、4的计算结果,计算第ξ个时刻风电场空气密度、湿度和气压的计算方法包括:
步骤5.1、根据风电场第ξ个时刻的输出功率Pwfrtξ与风速Wssrtξ,计算第ξ个时刻的空气密度ρsrtξ
第ξ个时刻的空气密度ρsrtξ的计算表达式为:
其中F为风轮叶片旋转一周的扫掠面积;
步骤5.2、根据大气压力和空气密度计算公式,以及空气湿度经验公式,计算得到湿度和气压
第ξ个时刻的气压presrtξ的计算表达式为:
presrtξ=ρsrtξ×(273.15+Tsrtξ)×R
式中,R为气体常数,其值为287;
第ξ个时刻的湿度Hsrtξ的计算表达式为:
步骤6、汇总步骤3、4、5计算得到的数据,形成风电场第ξ时刻的虚拟测风塔数据。其中风电场第ξ时刻的虚拟测风塔数据的表达式为:
VAT={Wssrtξ、Wdsrtξ、Tsrtξ、ρsrtξ、Hsrtξ、presrtξ}。
Claims (2)
1.一种风电场虚拟测风塔构建方法,它包括:
步骤1、采集风电场n台风机的基础数据和历史数据;
步骤2、采集风电场n台风机实时数据,实时数据包括风电场所在区域内所有风机的风机测风仪测得风速数据Wssrt和风向数据Wdsrt、所有风机的风机温度传感器测得温度数据Tsrt、所有风机的实时出力数据Pwtrt、整个风电场的实时出力数据Pwfrt;
步骤3、根据风电场n台风机的历史数据和实时数据信息,模拟计算第ξ个时刻整个风电场轮毂高度的风速Wssrtξ和风向Wdsrtξ;
步骤3所述根据风电场n台风机的历史数据和实时数据信息,模拟计算第ξ个时刻整个风电场轮毂高度的风速Wssrtξ和风向Wdsrtξ的方法包括:
步骤3.1、利用风电场所有风电机组的历史风速风向数据,以及整个风电场的历史出力数据,建立所有风电机组风速风向-风电场风电功率数据样本库;
步骤3.2、建立所有风电机组风速风向-风电场风电功率模型;
步骤3.3、将3.1建立的所有风电机组风速风向-风电场风电功率数据样本库作为训练集,对风电机组风速风向-风电场风电功率模型进行训练,确定模型中n台风电机组的风速风向数据的权值;其中,第i台风机的风速风向数据的权值为Wqzi,则有Wqz1+Wqz2+...Wqzi+...+Wqzn=1;且有1≤i≤n;
步骤3.4、计算第ξ个时刻整个风电场轮毂高度的风速Wssrtξ和风向Wdsrtξ
第ξ个时刻整个风电场轮毂高度的风速Wssrtξ的计算公式为:
第ξ个时刻整个风电场轮毂高度的风向Wdsrtξ的表达式为:
步骤4、根据风电场n台风机的历史数据和实时数据信息,模拟计算第ξ个时刻整个风电场轮毂高度的温度Tsrtξ;
步骤4所述的根据风电场n台风机的历史数据和实时数据信息,模拟计算第ξ个时刻整个风电场轮毂高度的温度Tsrtξ的计算方法包括:
步骤4.1、利用风电场所有风电机组的历史温度数据,以及所有风机的历史出力数据和整个风电场的历史出力数据,建立所有风电机组温度-风电场风电功率数据样本库;
步骤4.2、建立所有风电机组温度-风电场风电功率模型;
步骤4.3、将4.1建立的所有风电机组温度-风电场风电功率数据样本库作为训练集,对风电机组温度-风电场风电功率模型进行训练,确定模型中所有风电机组的温度数据的权值;
其中,第i台风机的温度数据的权值记做Tqzi,则有Tqz1+Tqz2+...Tqzi+...+Tqzn=1;
且有1≤i≤n;
步骤4.4、模拟计算第ξ个时刻整个风电场轮毂高度的温度Tsrtξ第ξ个时刻整个风电场轮毂高度的温度Tsrtξ的表达式为:
步骤5、根据步骤3、4的计算结果,计算第ξ个时刻风电场空气密度、湿度和气压;
步骤5所述根据步骤3、4的计算结果,计算第ξ个时刻风电场空气密度、湿度和气压的计算方法包括:
步骤5.1、根据风电场第ξ个时刻的输出功率Pwfrtξ与风速Wssrtξ,计算第ξ个时刻的空气密度ρsrtξ
第ξ个时刻的空气密度ρsrtξ的计算表达式为:
其中F为风轮叶片旋转一周的扫掠面积;
步骤5.2、根据大气压力和空气密度计算公式,以及空气湿度公式,计算得到湿度和气压
第ξ个时刻的气压presrtξ的计算表达式为:
presrtξ=ρsrtξ×(273.15+Tsrtξ)×R
式中,R为气体常数,其值为287;
第ξ个时刻的湿度Hsrtξ的计算表达式为:
步骤6、汇总步骤3、4、5计算得到的数据,形成风电场第ξ时刻的虚拟测风塔数据。
2.根据权利要求1所述的风电场虚拟测风塔构建方法,其特征在于:
步骤1所述基础数据包括风电场区域内所有风机的轮毂高度、历史数据包括风电场所在区域内所有风机的历史风速数据Wssh和风向数据Wdsh、所有风机的历史温度数据Tsh、所有风机的历史出力数据Pwth和整个风电场的历史出力数据Pwfh。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611119467.3A CN106547996B (zh) | 2016-12-08 | 2016-12-08 | 一种风电场虚拟测风塔构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611119467.3A CN106547996B (zh) | 2016-12-08 | 2016-12-08 | 一种风电场虚拟测风塔构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106547996A CN106547996A (zh) | 2017-03-29 |
CN106547996B true CN106547996B (zh) | 2019-10-29 |
Family
ID=58396698
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611119467.3A Active CN106547996B (zh) | 2016-12-08 | 2016-12-08 | 一种风电场虚拟测风塔构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106547996B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108665102B (zh) * | 2018-05-11 | 2022-08-05 | 中国船舶重工集团海装风电股份有限公司 | 基于中尺度数据预测风电场实时发电量的方法 |
CN109784563B (zh) * | 2019-01-18 | 2023-05-23 | 南方电网科学研究院有限责任公司 | 一种基于虚拟测风塔技术的超短期功率预测方法 |
CN111090932B (zh) * | 2019-12-10 | 2023-06-09 | 华能威宁风力发电有限公司 | 一种适应中等复杂地形的在运风电场场地标定方法 |
CN111932400B (zh) * | 2020-08-06 | 2021-04-30 | 南通大学 | 风/储一体功率调度计划优化实现方法 |
CN114139775A (zh) * | 2021-11-15 | 2022-03-04 | 明阳智慧能源集团股份公司 | 一种基于虚拟测风塔的风电功率预测系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104504464A (zh) * | 2014-12-11 | 2015-04-08 | 国家电网公司 | 基于风区风带风速规律的风功率预测方法 |
CN105844361A (zh) * | 2016-03-29 | 2016-08-10 | 北京金风科创风电设备有限公司 | 风电机组风功率预测方法、解缆方法和装置 |
-
2016
- 2016-12-08 CN CN201611119467.3A patent/CN106547996B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104504464A (zh) * | 2014-12-11 | 2015-04-08 | 国家电网公司 | 基于风区风带风速规律的风功率预测方法 |
CN105844361A (zh) * | 2016-03-29 | 2016-08-10 | 北京金风科创风电设备有限公司 | 风电机组风功率预测方法、解缆方法和装置 |
Non-Patent Citations (2)
Title |
---|
基于人工神经网络的风电场短期功率预测;彭怀午等;《太阳能学报》;20110831;第32卷(第8期);第1247-1249页 * |
风电功率预测方法研究;范强等;《贵州电力技术》;20161031;第19卷(第10期);第13-16页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106547996A (zh) | 2017-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106547996B (zh) | 一种风电场虚拟测风塔构建方法 | |
CN106779208B (zh) | 一种基于虚拟测风塔技术的风电超短期功率预测方法 | |
Wagner et al. | The influence of the wind speed profile on wind turbine performance measurements | |
Iungo et al. | Volumetric lidar scanning of wind turbine wakes under convective and neutral atmospheric stability regimes | |
CN104699936B (zh) | 基于计算流体力学短期风速预测的风电场扇区管理方法 | |
Li et al. | LiDAR technology for wind energy potential assessment: Demonstration and validation at a site around Lake Erie | |
CN109958588A (zh) | 结冰预测方法、装置、存储介质、模型生成方法及装置 | |
CN106815456A (zh) | 一种风电机组功率特性评价方法 | |
CN105911467B (zh) | 复杂地形下的风电机组功率曲线考核评估方法 | |
Hau et al. | The wind resource | |
Wharton et al. | Measuring turbine inflow with vertically-profiling lidar in complex terrain | |
Radünz et al. | The variability of wind resources in complex terrain and its relationship with atmospheric stability | |
Rodrigo et al. | Iea-task 31 wakebench: Towards a protocol for wind farm flow model evaluation. part 1: Flow-over-terrain models | |
CN115358606B (zh) | 一种平缓地形下在役风电场能效评估方法及系统 | |
CN108062722B (zh) | 基于风速变异系数的山地风电场模型风机机械功率计算方法 | |
CN105335617B (zh) | 风电场尾流效应评估方法及装置 | |
CN109784563A (zh) | 一种基于虚拟测风塔技术的超短期功率预测方法 | |
CN106815773A (zh) | 一种风电场功率特性评估方法 | |
CN109636019B (zh) | 基于神经网络算法的测风塔布置方案确定方法 | |
Antoniou et al. | Influence of wind characteristics on turbine performance | |
Li et al. | Wind velocity prediction at wind turbine hub height based on CFD model | |
CN115310375A (zh) | 基于流体动力学模型的风电机组机舱传递函数拟合方法 | |
CN214944743U (zh) | 一种海上风电机组载荷测试数据采集系统 | |
Björck et al. | Dynamic stall and 3D effects | |
Gasch et al. | The wind |
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 |