CN111709112A - 一种海上风电运行模拟方法、装置及存储介质 - Google Patents
一种海上风电运行模拟方法、装置及存储介质 Download PDFInfo
- Publication number
- CN111709112A CN111709112A CN202010370948.1A CN202010370948A CN111709112A CN 111709112 A CN111709112 A CN 111709112A CN 202010370948 A CN202010370948 A CN 202010370948A CN 111709112 A CN111709112 A CN 111709112A
- Authority
- CN
- China
- Prior art keywords
- wind speed
- wind power
- corrected
- electricity generation
- power output
- 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
- 238000004088 simulation Methods 0.000 title claims abstract description 65
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000005611 electricity Effects 0.000 claims abstract description 84
- 230000000694 effects Effects 0.000 claims abstract description 83
- 238000012795 verification Methods 0.000 claims abstract description 27
- 238000006243 chemical reaction Methods 0.000 claims description 27
- 238000013461 design Methods 0.000 claims description 27
- 238000004364 calculation method Methods 0.000 claims description 19
- 238000012937 correction Methods 0.000 claims description 17
- 238000005315 distribution function Methods 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 10
- 230000005653 Brownian motion process Effects 0.000 claims description 7
- 238000005537 brownian motion Methods 0.000 claims description 7
- 230000006870 function Effects 0.000 description 22
- 238000010248 power generation Methods 0.000 description 16
- 238000010586 diagram Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000003203 everyday effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000001932 seasonal effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 230000002354 daily effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012546 transfer 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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Wind Motors (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及电力系统规划技术领域,公开了一种海上风电运行模拟方法、装置及存储介质,所述方法包括:S1:获取海上风速模拟值;S2:对所述海上风速模拟值进行修正,获得修正风速;S3:根据所述修正风速获得风电出力;S4:根据当前的尾流效应系数对所述风电出力进行修正,获得修正风电出力;S5:根据所述修正风电出力获得发电利用小时数;S6:对所述发电利用小时数进行校验;S7:当所述发电利用小时数不满足校验要求时,对所述当前的尾流效应系数进行修正,并返回步骤S4。本发明提供的一种海上风电运行模拟方法、装置及存储介质,能够使风电运行模拟的发电利用小时数接近小时数设计值,保证海上风电场出力曲线模拟值的准确性。
Description
技术领域
本发明涉及电力系统规划技术领域,特别是涉及一种海上风电运行模拟方法、装置及存储介质。
背景技术
海上风电具有资源丰富、发电利用小时数相对较高、技术相对高端的特点,是新能源发展的前沿领域,是我国沿海地区可再生能源中最具规模化发展潜力的领域。
随着系统中海上风电装机规模和比重的不断增大,海上风电出力对电力系统的影响将越来越显著,海上风电出力特性分析工作越来越重要。海上风电出力随机性、间歇性很强,受海域、气候等影响很大。风电出力除了幅值上的随机性以外,还存在明显的时间相依性,包括风电出力发电与负荷的相关性,风电场发电出力的自相关性等,这给开展海上风电出力特性分析工作带来挑战。如何有效模拟出未来规划年份的海上风电出力曲线及出力特性,是开展海上风电运行模拟工作的核心内容。
目前,海上风电运行模拟方法主要分为两大类:(1)模拟风速序列,进而转换成风电出力序列;(2)通过对历史风电出力曲线进行聚类分析,得出典型风电出力曲线。
上述方法侧重于通过建模方法模拟风电出力曲线,在模拟时重点关注如何在模型方法中考虑随机特性、波动特性、时间相依性、空间相关性等风电特性。发电利用小时数指标是风电场规划时评价风场资源好坏的重要指标,也是衡量风电运行模拟结果准确性的重要指标。现有的模拟方法较少关注风电出力特性中的发电利用小时数这一指标,未有效利用风电场规划设计时的发电利用小时数信息。
发明内容
本发明实施例所要解决的技术问题是:提供一种海上风电运行模拟方法、装置及存储介质,利用风电场规划设计时的发电利用小时数信息,构建考虑尾流效应系数的海上风电出力表达式,通过修正尾流效应系数的方式,使得模拟的发电利用小时数接近发电利用小时数设计值,保证海上风电场出力曲线模拟值的准确性。
为了解决上述技术问题,第一方面,本发明实施例提供一种海上风电运行模拟方法,所述方法包括:
S1:获取海上风速模拟值;
S2:对所述海上风速模拟值进行修正,获得修正风速;
S3:根据所述修正风速获得风电出力;
S4:根据当前的尾流效应系数对所述风电出力进行修正,获得修正风电出力;
S5:根据所述修正风电出力获得发电利用小时数;
S6:对所述发电利用小时数进行校验;
S7:当所述发电利用小时数不满足校验要求时,对所述当前的尾流效应系数进行修正,并返回步骤S4。
作为一个优选方案,所述获取海上风速模拟值,具体包括:
所述获取海上风速模拟值,具体包括:
获取历史风速序列以及满足标准布朗运动的第一随机变量、第二随机变量;
根据所述历史风速序列获得历史风速平均值、衰减系数以及风速Weibull分布函数;其中,所述风速Weibull分布函数的表达式为:c为Weibull分布尺度参数,k为Weibull分布形状参数,vt为t时刻的海上风速模拟值;
根据所述历史风速平均值、所述衰减系数、所述风速Weibull分布函数获得风速转换函数;其中,所述风速转换函数的表达式为:θ为所述衰减系数,μ为所述历史风速平均值,l为定义域左端点,y为定义域上的变量;
根据所述第一随机变量、所述第二随机变量、所述历史风速平均值、所述衰减系数以及所述转换函数获取所述海上风速模拟值;其中,所述海上风速模拟值的表达式为:vt+1为t+1时刻的海上风速模拟值,Wt为所述第一随机变量,Wt+1为所述第二随机变量。
作为一个优选方案,所述对所述海上风速模拟值进行修正,获得修正风速,具体包括:
根据所述历史风速序列获得海上风电场的月规律因子、小时规律因子;
根据所述月规律因子以及所述小时规律因子对所述海上风速模拟值进行修改正,获得所述修正风速;其中,所述修正风速的表达式为:vt *=vtkhkm,vt *为所述修正风速,km为所述月规律因子,kh为所述小时规律因子。
作为一个优选方案,所述根据所述修正风速获得风电出力,具体包括:
获取海上风电场的额定功率;
获取风速与风功率转换函数;其中,所述风速与风功率转换函数的表达式为:
其中,vci为切入风速,vr为额定风速,vco为切出风速;
根据所述额定功率、所述风速与风功率转换函数获得风电出力;其中,所述风电出力的表达式为:Pt=f(t)Pr,Pt为t时刻的风电出力,Pr为所述额定功率。
作为一个优选方案,所述根据当前的尾流效应系数对所述风电出力进行修正,获得修正风电出力,具体包括:
获取海上风电场的可用机组台数;
获取海上风电场当前的尾流效应系数;
根据所述风电出力、所述可用机组台数以及所述当前的尾流效应系数,获得修正风电出力;其中,所述修正出力的表达式为:Pt*=n(1-η)Pt,Pt*为t时刻的修正风电出力,n为所述可用机组台数,η为所述当前的尾流效应系数。
作为一个优选方案,所述根据所述修正风电出力获得发电利用小时数,具体为:
作为一个优选方案,所述对所述发电利用小时数进行校验,具体包括:
将所述发电利用小时数与小时数设计值进行对比;
当所述发电利用小时数与所述小时数设计值的差值不大于预设范围时,判定所述发电利用小时数满足校验要求;
当所述发电利用小时数与所述小时数设计值的差值大于预设范围时,判定所述发电利用小时数不满足校验要求。
作为一个优选方案,所述对所述当前的尾流效应系数进行修正,具体包括:
将η=η'作为当前的尾流效应系数。
为了解决上述技术问题,第二方面,本发明实施例提供一种海上风电运行模拟装置,所述装置包括:
风速获取模块,用于获取海上风速模拟值;
风速修正模块,用于对所述海上风速模拟值进行修正,获得修正风速;
出力计算模块,用于根据所述修正风速获得风电出力;
出力修正模块,用于根据当前的尾流效应系数对所述风电出力进行修正,获得修正风电出力;
小时数计算模块,用于根据所述修正风电出力获得发电利用小时数;
小时数校验模块,用于对所述发电利用小时数进行校验;
尾流效应系数修正模块,用于进行尾流效应系数修正。
为了解决上述技术问题,第三方面,本发明实施例提供一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被执行时实现如第一方面任一项所述的海上风电运行模拟方法。
与现有技术相比,本发明实施例提供的一种海上风电运行模拟方法、装置及存储介质,其有益效果在于:在研究海上风电运行模拟方法时重点关注发电利用小时数指标,有效利用风电场规划设计时的发电利用小时数信息,提出了一种考虑发电利用小时数的海上风电运行模拟方法,将影响发电利用小时数模拟值的尾流效应系数与风电出力关联起来,构建考虑尾流效应系数的海上风电出力表达式,给出了尾流效应系数的调整方法,能够做到通过调整尾流效应系数的方式使模拟的发电利用小时数接近发电利用小时数设计值,保证海上风电场出力曲线模拟值的准确性;而且该方法原理简单,计算速度快,易于实施。
附图说明
为了更清楚地说明本发明实施例的技术特征,下面将对本发明实施例中所需要使用的附图做简单地介绍,显而易见地,下面所描述的附图仅仅是本发明的一些实施例,对于本领域技术人员来说,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明提供的一种海上风电运行模拟方法的一个优选实施例的流程示意图;
图2是本发明提供的一种海上风电运行模拟装置的一个优选实施例的结构示意图。
具体实施方式
为了对本发明的技术特征、目的、效果有更加清楚的理解,下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例仅用于说明本发明,但是不用来限制本发明的保护范围。基于本发明的实施例,本领域技术人员在没有付出创造性劳动的前提下所获得的其他实施例,都应属于本发明的保护范围。
在本发明的描述中,应当理解的是,本文中的编号本身,例如“第一”、“第二”等,仅用于区分所描述的对象,不具有顺序或者技术含义,不能理解为规定或者暗示所描述的对象的重要性。
图1所示为本发明提供的一种海上风电运行模拟方法的一个优选实施例的流程示意图。
如图1所示,所述方法包括:
S1:获取海上风速模拟值;
S2:对所述海上风速模拟值进行修正,获得修正风速;
S3:根据所述修正风速获得风电出力;
S4:根据当前的尾流效应系数对所述风电出力进行修正,获得修正风电出力;
S5:根据所述修正风电出力获得发电利用小时数;
S6:对所述发电利用小时数进行校验;
S7:当所述发电利用小时数不满足校验要求时,对所述当前的尾流效应系数进行修正,并返回步骤S4。
需要说明的是,每个海上风电系统都会包括多个风电场,每一个风电场又会包括多个机组,为了便于说明,本发明实施例仅选取其中一个风电场来实施本发明,但本发明实施例的保护范围不限于此。
具体实施时,首先对海上风电场进行运行模拟,获取海上风速模拟值,由于风速在不同季节和每日内不同时刻会呈现一定的大小规律,为反映风电场的月平均出力水平和小时平均出力水平,需要对所述风速模拟值进行修正,获取修正风速,并根据所述修正风速计算获得风电出力,考虑到风场的尾流效应与风电场风机出力可靠性,需要根据当前的尾流效应系数对所述风电出力进行修正,获得修正风电出力,并根据所述修正风电出力获得风电场的发电利用小时数,将所述发电利用小时数与风电场规划设计时的发电利用小时设计值进行对比做校验,当所述发电小时数满足校验要求时,则输出海上风电计算结果,当所述发电利用小时数不满足校验要求时,则需要对所述当前的尾流效应系数进行修正,并重新计算修正出力,直至所述发电利用小时数满足校验要求为止。
本发明实施例提供的一种海上风电运行模拟方法,在研究海上风电运行模拟时重点关注发电利用小时数指标,有效利用风电场规划设计时的发电利用小时数信息,将影响发电利用小时数模拟值的尾流效应系数与风电出力关联起来,构建考虑尾流效应系数的海上风电出力表达式,能够通过调整尾流效应系数的方式使模拟的发电利用小时数接近发电利用小时数设计值,保证海上风电场出力曲线模拟值的准确性。
在一个优选方案中,所述获取海上风速模拟值,具体包括:
获取历史风速序列以及满足标准布朗运动的第一随机变量、第二随机变量;
根据所述历史风速序列获得历史风速平均值、衰减系数以及风速Weibull分布函数;其中,所述风速Weibull分布函数的表达式为:c为Weibull分布尺度参数,k为Weibull分布形状参数,vt为t时刻的海上风速模拟值;
根据所述第一随机变量、所述第二随机变量、所述历史风速平均值、所述衰减系数以及所述转换函数获取所述海上风速模拟值;其中,所述海上风速模拟值的表达式为:vt+1为t+1时刻的海上风速模拟值,Wt为所述第一随机变量,Wt+1为所述第二随机变量。
需要说明的是,为海上风速模拟核心公式,式中:t为时段编号,1≤t≤T,T为时段数量,一般T=8760,vt为t时刻的海上风速模拟值,vt+1为t+1时刻的海上风速模拟值,θ为衰减系数,μ为风速平均值,w(vt)为风速有关的转换函数,Wt为满足标准布朗运动的第一随机变量,Wt+1为满足标准布朗运动的第二随机变量。
具体而言,其中各个参数的获取方式如下:
1)获取历史风速序列,一般为全年小时级风速序列;
2)获取满足标准布朗运动的第一随机变量、第二随机变量;
3)通过计算历史风速平均值得到历史风速平均值;
4)通过计算历史风速自相关系数以后拟合得到衰减系数;
5)通过计算历史风速的Weibull分布尺度参数c、形状参数k,得到风速Weibull分布函数;
6)根据所述衰减系数、所述风速Weibull分布函数计算得到风速转换函数。
最后,将上述计算得到的各参数带入海上风速模拟公式获取所述海上风速模拟值。
在一个优选方案中,所述对所述海上风速模拟值进行修正,获得修正风速,具体包括:
根据所述历史风速序列获得海上风电场的月规律因子、小时规律因子;
根据所述月规律因子以及所述小时规律因子对所述海上风速模拟值进行修改正,获得所述修正风速;其中,所述修正风速的表达式为:vt *=vtkhkm,kh为所述小时规律因子,km为所述月规律因子。
由于风速在不同季节和每日内不同时刻会呈现一定的大小规律,为反映风电场的月平均出力水平和小时平均出力水平,需要对所述海上风速模拟值进行修正,获得修正风速,进一步得到满足海上风电场季节特性和日特性的修正风速曲线。
具体而言,满足海上风电场季节特性的月规律因子km、日特征的小时规律因子kh是根据所述历史风速序列计算获得的,获得所述月规律因子、所述小时规律因子后,再对所述海上风速模拟值进行修正获得所述修正风速。
在一个优选方案中,所述根据所述修正风速获得风电出力,具体包括:
获取海上风电场的额定功率;
获取风速与风功率转换函数;其中,所述风速与风功率转换函数的表达式为:
其中,vci为切入风速,vr为额定风速,vco为切出风速;
根据所述额定功率、所述风速与风功率转换函数获得风电出力;其中,所述风电出力的表达式为:Pt=f(t)Pr,Pr为所述额定功率。
具体而言,需要预先根据风机的出力特性曲线计算出风速与风功率转换函数,在计算风电出力时,还要获取海上风电场的额定功率、额定风速以及风速模拟值的切入风速、切出风速,将所述额定功率与所述风速与风功率转换函数相乘得到所述风电出力。
在一个优选方案中,所述根据当前的尾流效应系数对所述风电出力进行修正,获得修正风电出力,具体包括:
获取海上风电场的可用机组台数;
获取海上风电场当前的尾流效应系数;
根据所述风电出力、所述可用机组台数以及所述当前的尾流效应系数,获得修正风电出力;其中,所述修正出力的表达式为:Pt*=n(1-η)Pt,Pt*为t时刻的修正风电出力,n为所述可用机组台数,η为所述当前的尾流效应系数。
具体实施时,考虑风场尾流效应与风电场风机出力可靠性,需要先获取可用机组台数以及尾流效应系数,并进一步求出整个风电场的修正风电出力。
在一个优选方案中,所述根据所述修正风电出力获得发电利用小时数,具体为:
具体而言,首先对一个发电时段的所有时刻的修正出力进行求和计算获得总修正出力,并将所述总风电出力除以额定功率获得发电利用小时数。
在一个优选方案中,所述对所述发电利用小时数进行校验,具体包括:
将所述发电利用小时数与小时数设计值进行对比;
当所述发电利用小时数与所述小时数设计值的差值不大于预设范围时,判定所述发电利用小时数满足校验要求;
当所述发电利用小时数与所述小时数设计值的差值大于预设范围时,判定所述发电利用小时数不满足校验要求。
需要说明的是,所述预设范围可以根据风电运行模拟的实际需求进行设置,将其设置得较大,则可以减少计算量,将其设置得较少,则模拟得到的出力曲线越准确,本发明实施例对此不做限制。
具体实施时,将所述发电利用小时数与风电场规划时的发电利用小时数设计值进行对比,若满足校验要求,则输出此时的海上风电计算结果,并结束计算;若不满足校验要求,则要对尾流效应系数进行修正。
在一个优选方案中,所述对所述当前的尾流效应系数进行修正,具体包括:
将η=η'作为当前的尾流效应系数。
具体而言,当所述发电利用小时数不满足校验要求时,对尾流效应系数进行修正调整,获得修正尾流效应系数,并基于所述修正尾流效应系数重新进行所述发电利用小时数的计算,直至所述发电利用小时数满足校验要求为止,并输出海上风电计算结果,结束计算。
综上,本发明实施例提供一种海上风电运行模拟方法,在研究海上风电运行模拟时重点关注发电利用小时数指标,有效利用了风电场规划设计时的发电利用小时数信息,将影响发电利用小时数模拟值的尾流效应系数与风电出力关联起来,构建考虑尾流效应系数的海上风电出力表达式,能够通过调整尾流效应系数的方式使模拟的发电利用小时数接近发电利用小时数设计值,保证海上风电场出力曲线模拟值的准确性;而且尾流效应系数的调整方法原理简单,计算速度快,易于实施。
应当理解,本发明实现上述海上风电运行模拟方法中的全部或部分流程,也可以通过计算机程序来指令相关的硬件来完成,计算机程序可存储于一计算机可读存储介质中,该计算机程序在被处理器执行时,可实现上述海上风电运行模拟方法的步骤。其中,计算机程序包括计算机程序代码,计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。计算机可读介质可以包括:能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、电载波信号、电信信号以及软件分发介质等。需要说明的是,计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减,例如在某些司法管辖区,根据立法和专利实践,计算机可读介质不包括电载波信号和电信信号。
图2所述为本发明提供的一种海上风电运行模拟装置的一个优选实施例的结构示意图。
如图2所示,所述装置包括:
风速获取模块,用于获取海上风速模拟值;
风速修正模块,用于对所述海上风速模拟值进行修正,获得修正风速;
出力计算模块,用于根据所述修正风速获得风电出力;
出力修正模块,用于根据当前的尾流效应系数对所述风电出力进行修正,获得修正风电出力;
小时数计算模块,用于根据所述修正风电出力获得发电利用小时数;
小时数校验模块,用于对所述发电利用小时数进行校验;
尾流效应系数修正模块,用于进行尾流效应系数修正。
具体实施时,所述风速获取模块首先对海上风电场进行模拟,获取海上风速模拟值,由于风速在不同季节和每日内不同时刻会呈现一定的大小规律,为反映风电场的月平均出力水平和小时平均出力水平,需要在所述风速修正模块对所述风速模拟值进行修正,获取修正风速,所述出力计算模块根据所述修正风速计算获得风电出力,考虑到风场的尾流效应与风电场风机出力可靠性,需要所述出力修正模块根据当前的尾流效应系数对所述风电出力进行修正,获得修正风电出力,所述小时数计算模块根据所述修正风电出力获得风电场的发电利用小时数,所述小时数校验模块将所述发电利用小时数与风电场规划设计时的发电利用小时设计值进行对比做校验,当所述发电小时数满足校验要求时,则所述装置输出海上风电计算结果,当所述发电利用小时数不满足校验要求时,则需要所述尾流效应系数修正模块对所述当前的尾流效应系数进行修正,并将修正后的尾流效应系数发送至所述出力模块修正重新计算修正出力,并进一步在所述小时数计算模块计算发电利用小时数,直至所述发电利用小时数在所述小时数校验模块满足校验要求为止。
本发明实施提供的一种海上风电运行模拟装置,能够通过调整尾流效应系数的方式使模拟的发电利用小时数接近发电利用小时数设计值,保证海上风电场出力曲线模拟值的准确性。
优选地,所述风速获取模块具体用于:
获取历史风速序列以及满足标准布朗运动的第一随机变量、第二随机变量;
根据所述历史风速序列获得历史风速平均值、衰减系数以及风速Weibull分布函数;其中,所述风速Weibull分布函数的表达式为:c为Weibull分布尺度参数,k为Weibull分布形状参数,vt为t时刻的海上风速模拟值;
根据所述第一随机变量、所述第二随机变量、所述历史风速平均值、所述衰减系数以及所述转换函数获取所述海上风速模拟值;其中,所述海上风速模拟值的表达式为:vt+1为t+1时刻的海上风速模拟值,Wt为所述第一随机变量,Wt+1为所述第二随机变量。
优选地,所述风速修正模块具体用于:
根据所述历史风速序列获得海上风电场的月规律因子、小时规律因子;
根据所述月规律因子以及所述小时规律因子对所述海上风速模拟值进行修改正,获得所述修正风速;其中,所述修正风速的表达式为:vt *=vtkhkm,kh为所述小时规律因子,km为所述月规律因子。
优选地,所述出力计算模块具体用于:
获取海上风电场的额定功率;
获取风速与风功率转换函数;其中,所述风速与风功率转换函数的表达式为:
其中,vci为切入风速,vr为额定风速,vco为切出风速;
根据所述额定功率、所述风速与风功率转换函数获得风电出力;其中,所述风电出力的表达式为:Pt=f(t)Pr,Pr为所述额定功率。
优选地,所述出力修正模块具体用于:
获取海上风电场的可用机组台数;
获取海上风电场当前的尾流效应系数;
根据所述风电出力、所述可用机组台数以及所述当前的尾流效应系数,获得修正风电出力;其中,所述修正出力的表达式为:Pt*=n(1-η)Pt,Pt*为t时刻的修正风电出力,n为所述可用机组台数,η为所述当前的尾流效应系数。
优选地,所述发电利用小时数获取模块具体用于:
优选地,所述发电利用小时数校验模块具体用于:
将所述发电利用小时数与小时数设计值进行对比;
当所述发电利用小时数与所述小时数设计值的差值不大于预设范围时,判定所述发电利用小时数满足校验要求;
当所述发电利用小时数与所述小时数设计值的差值大于预设范围时,判定所述发电利用小时数不满足校验要求。
优选地,所述尾流效应系数修正模块具体用于:
将η=η'作为当前的尾流效应系数。
综上,本发明实施例提供一种海上风电运行模拟装置,有效利用了风电场规划设计时的发电利用小时数信息,构建考虑尾流效应系数的海上风电出力表达式,能够通过调整尾流效应系数的方式使模拟的发电利用小时数接近发电利用小时数设计值,保证海上风电场出力曲线模拟值的准确性;而且尾流效应系数的调整方法原理简单,计算速度快,易于实施。
以上所述,仅是本发明的优选实施方式,但本发明的保护范围并不局限于此,应当指出,对于本领域技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干等效的明显变形和/或同等替换,这些明显变形和/或同等替换也应视为本发明的保护范围。
Claims (10)
1.一种海上风电运行模拟方法,其特征在于,所述方法包括:
S1:获取海上风速模拟值;
S2:对所述海上风速模拟值进行修正,获得修正风速;
S3:根据所述修正风速获得风电出力;
S4:根据当前的尾流效应系数对所述风电出力进行修正,获得修正风电出力;
S5:根据所述修正风电出力获得发电利用小时数;
S6:对所述发电利用小时数进行校验;
S7:当所述发电利用小时数不满足校验要求时,对所述当前的尾流效应系数进行修正,并返回步骤S4。
2.根据权利要求1所述的海上风电运行模拟方法,其特征在于,所述获取海上风速模拟值,具体包括:
获取历史风速序列以及满足标准布朗运动的第一随机变量、第二随机变量;
根据所述历史风速序列获得历史风速平均值、衰减系数以及风速Weibull分布函数;其中,所述风速Weibull分布函数的表达式为:c为Weibull分布尺度参数,k为Weibull分布形状参数,vt为t时刻的海上风速模拟值;
根据所述历史风速平均值、所述衰减系数、所述风速Weibull分布函数获得风速转换函数;其中,所述风速转换函数的表达式为:θ为所述衰减系数,μ为所述历史风速平均值,l为定义域左端点,y为定义域上的变量;
3.根据权利要求2所述的海上风电运行模拟方法,其特征在于,所述对所述海上风速模拟值进行修正,获得修正风速,具体包括:
根据所述历史风速序列获得海上风电场的月规律因子、小时规律因子;
根据所述月规律因子以及所述小时规律因子对所述海上风速模拟值进行修改正,获得所述修正风速;其中,所述修正风速的表达式为:vt *=vtkhkm,vt *为所述修正风速,km为所述月规律因子,kh为所述小时规律因子。
5.根据权利要求4所述的海上风电运行模拟方法,其特征在于,所述根据当前的尾流效应系数对所述风电出力进行修正,获得修正风电出力,具体包括:
获取海上风电场的可用机组台数;
获取海上风电场当前的尾流效应系数;
根据所述风电出力、所述可用机组台数以及所述当前的尾流效应系数,获得修正风电出力;其中,所述修正出力的表达式为:Pt*=n(1-η)Pt,Pt*为t时刻的修正风电出力,n为所述可用机组台数,η为所述当前的尾流效应系数。
7.根据权利要求6所述的海上风电运行模拟方法,其特征在于,所述对所述发电利用小时数进行校验,具体包括:
将所述发电利用小时数与小时数设计值进行对比;
当所述发电利用小时数与所述小时数设计值的差值不大于预设范围时,判定所述发电利用小时数满足校验要求;
当所述发电利用小时数与所述小时数设计值的差值大于预设范围时,判定所述发电利用小时数不满足校验要求。
9.一种海上风电运行模拟装置,其特征在于,所述装置包括:
风速获取模块,用于获取海上风速模拟值;
风速修正模块,用于对所述海上风速模拟值进行修正,获得修正风速;
出力计算模块,用于根据所述修正风速获得风电出力;
出力修正模块,用于根据当前的尾流效应系数对所述风电出力进行修正,获得修正风电出力;
小时数计算模块,用于根据所述修正风电出力获得发电利用小时数;
小时数校验模块,用于对所述发电利用小时数进行校验;
尾流效应系数修正模块,用于进行尾流效应系数修正。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被执行时实现如权利要求1至8任一项所述的海上风电运行模拟方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010370948.1A CN111709112B (zh) | 2020-04-30 | 2020-04-30 | 一种海上风电运行模拟方法、装置及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010370948.1A CN111709112B (zh) | 2020-04-30 | 2020-04-30 | 一种海上风电运行模拟方法、装置及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111709112A true CN111709112A (zh) | 2020-09-25 |
CN111709112B CN111709112B (zh) | 2023-05-16 |
Family
ID=72536820
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010370948.1A Active CN111709112B (zh) | 2020-04-30 | 2020-04-30 | 一种海上风电运行模拟方法、装置及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111709112B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112883577A (zh) * | 2021-02-26 | 2021-06-01 | 广东电网有限责任公司 | 一种海上风电场出力典型场景生成方法及存储介质 |
CN113326636A (zh) * | 2021-06-30 | 2021-08-31 | 广东电网有限责任公司 | 一种远海海上风电场时序出力曲线快速生成方法及系统 |
CN113904338A (zh) * | 2021-11-19 | 2022-01-07 | 广东电网有限责任公司 | 风电并网系统计及频率特性的概率潮流计算方法及系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103617308A (zh) * | 2013-10-30 | 2014-03-05 | 河海大学 | 一种风电场频域等效模型的构建方法 |
CN104133989A (zh) * | 2014-07-15 | 2014-11-05 | 华北电力大学 | 计及覆冰损失的风电场时序输出功率计算方法 |
CN104794325A (zh) * | 2015-03-10 | 2015-07-22 | 国家电网公司 | 一种基于随机差分方程的集群风电场出力时序模拟方法 |
US20150240789A1 (en) * | 2014-02-25 | 2015-08-27 | State Grid Corporation Of China | Method of analyzing wake flow of wind turbine based on multiple wake flow models |
CN106951977A (zh) * | 2017-02-08 | 2017-07-14 | 中能电力科技开发有限公司 | 一种基于尾流效应的风速预测模型的构建方法 |
CN108039737A (zh) * | 2017-12-29 | 2018-05-15 | 国网能源研究院有限公司 | 一种源网荷协调运行模拟系统 |
CN109614718A (zh) * | 2018-12-14 | 2019-04-12 | 国网山东省电力公司经济技术研究院 | 一种风电出力仿真模拟模型与方法 |
CN111062617A (zh) * | 2019-12-18 | 2020-04-24 | 广东电网有限责任公司电网规划研究中心 | 一种海上风电出力特性分析方法及系统 |
-
2020
- 2020-04-30 CN CN202010370948.1A patent/CN111709112B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103617308A (zh) * | 2013-10-30 | 2014-03-05 | 河海大学 | 一种风电场频域等效模型的构建方法 |
US20150240789A1 (en) * | 2014-02-25 | 2015-08-27 | State Grid Corporation Of China | Method of analyzing wake flow of wind turbine based on multiple wake flow models |
CN104133989A (zh) * | 2014-07-15 | 2014-11-05 | 华北电力大学 | 计及覆冰损失的风电场时序输出功率计算方法 |
CN104794325A (zh) * | 2015-03-10 | 2015-07-22 | 国家电网公司 | 一种基于随机差分方程的集群风电场出力时序模拟方法 |
CN106951977A (zh) * | 2017-02-08 | 2017-07-14 | 中能电力科技开发有限公司 | 一种基于尾流效应的风速预测模型的构建方法 |
CN108039737A (zh) * | 2017-12-29 | 2018-05-15 | 国网能源研究院有限公司 | 一种源网荷协调运行模拟系统 |
CN109614718A (zh) * | 2018-12-14 | 2019-04-12 | 国网山东省电力公司经济技术研究院 | 一种风电出力仿真模拟模型与方法 |
CN111062617A (zh) * | 2019-12-18 | 2020-04-24 | 广东电网有限责任公司电网规划研究中心 | 一种海上风电出力特性分析方法及系统 |
Non-Patent Citations (2)
Title |
---|
陈启超: "促进可再生能源消纳的多能互补时序运行模拟技术" * |
陈颖;周海;王文鹏;曹潇;丁杰;: "风电场输出功率超短期预测结果分析与改进" * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112883577A (zh) * | 2021-02-26 | 2021-06-01 | 广东电网有限责任公司 | 一种海上风电场出力典型场景生成方法及存储介质 |
CN113326636A (zh) * | 2021-06-30 | 2021-08-31 | 广东电网有限责任公司 | 一种远海海上风电场时序出力曲线快速生成方法及系统 |
CN113904338A (zh) * | 2021-11-19 | 2022-01-07 | 广东电网有限责任公司 | 风电并网系统计及频率特性的概率潮流计算方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111709112B (zh) | 2023-05-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Karki et al. | A simplified wind power generation model for reliability evaluation | |
Olauson et al. | Modelling the Swedish wind power production using MERRA reanalysis data | |
CN111709112B (zh) | 一种海上风电运行模拟方法、装置及存储介质 | |
Sorensen et al. | Power fluctuations from large wind farms | |
WO2017054537A1 (zh) | 一种长时间尺度光伏出力时间序列建模方法及装置 | |
Karki et al. | Wind power simulation model for reliability evaluation | |
Wangdee et al. | Considering load-carrying capability and wind speed correlation of WECS in generation adequacy assessment | |
Nuño et al. | On the simulation of aggregated solar PV forecast errors | |
CN104091077B (zh) | 一种风电场光伏电站联合出力模拟方法 | |
Todeschini et al. | Medium-term variability of the UK's combined tidal energy resource for a net-zero carbon grid | |
Reddy et al. | Hybrid approach for short term wind power forecasting | |
Kennedy et al. | A probabilistic model for simulating long-term wind-power output | |
Olauson et al. | Scenarios and time series of future wind power production in Sweden | |
Padhee et al. | Analyzing effects of seasonal variations in wind generation and load on voltage profiles | |
CN107992970B (zh) | 风力发电机组的输出功率预测方法和设备 | |
Salles et al. | Risk analysis methodologies for financial evaluation of wind energy power generation projects in the Brazilian system | |
Lang et al. | Verification of wind power forecasts provided in real-time to the Irish Transmission System Operator | |
Karki et al. | Transmission system adequacy evaluation considering wind power | |
CN111259552A (zh) | 一种大规模风电场概率型时序出力模拟方法和装置 | |
Fernandez‐Bernal et al. | Wind speed generation for dynamic analysis | |
Wannakam et al. | Evaluation of Wind Energy Production Using Weibull Distribution and Artificial Neural Networks | |
Pang et al. | Research on wind power output reconstruction technology based on multiple time scales | |
Khalkho et al. | Reliability Evaluation of Nondispatchable Energy Sources in Generation Planning: A Wind Electrical System Case Study | |
Hoevenaars | Temporal resolution in time series and probabilistic models of renewable power systems | |
Ulm et al. | Wind and Solar energy optimal integration |
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 |