CN104574216A - 一种基于wams数据的风电出力特性分析方法 - Google Patents
一种基于wams数据的风电出力特性分析方法 Download PDFInfo
- Publication number
- CN104574216A CN104574216A CN201510031408.XA CN201510031408A CN104574216A CN 104574216 A CN104574216 A CN 104574216A CN 201510031408 A CN201510031408 A CN 201510031408A CN 104574216 A CN104574216 A CN 104574216A
- Authority
- CN
- China
- Prior art keywords
- wind power
- wind
- power output
- load
- electricity generation
- 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.)
- Pending
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 7
- 230000005611 electricity Effects 0.000 claims description 93
- 230000000875 corresponding effect Effects 0.000 claims description 12
- 238000010248 power generation Methods 0.000 claims description 7
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims description 6
- 230000001174 ascending effect Effects 0.000 claims description 3
- 230000002596 correlated effect Effects 0.000 claims description 3
- 238000010219 correlation analysis Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 4
- 230000001360 synchronised effect Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000005684 electric field Effects 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000029087 digestion Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000009897 systematic effect Effects 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- 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)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Wind Motors (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明提供一种基于WAMS数据的风电出力特性分析方法,包括以下步骤:确定风电出力特性;分析风电出力特性;定义风电场群,并对其进行统计分析。本发明提供的基于WAMS数据的风电出力特性分析方法,对风电出力特性指标进行量化,从而能够为分析大规模风电接入对电网的影响提供技术参考。同时,通过对大量风电运行监测数据的积累分析,也有助于协助调度运行人员把握风电运行规律,为更加科学合理地消纳风电和保障电网安全稳定运行提供分析基础。
Description
技术领域
本发明涉及一种分析方法,具体涉及一种基于WAMS数据的风电出力特性分析方法。
背景技术
随着风电等新能源装机容量的不断持续增长,其逐渐成为我国能源战略的重要补充部分。然而,由于风力发电具有随机性、波动性、间歇性等固有特性,风电接入必然会给电网运行带来一定影响,尤其是大规模风电接入条件下其影响更大。目前,风电出力特性的研究大多为定性分析,缺乏系统性的风电出力特性指标体系建立,由此导致难以深入掌握大规模风电的运行规律和特性。
广域测量系统(WideAreaMeasurementSystem,WAMS)采用同步相角测量技术,通过逐步布局全网关键测点的同步相角测量单元(PMU),实现对全网同步相角及电网主要数据的实时高速率采集。WAMS通过电力调度数据网络实时传送到广域监测主站系统,从而提供对电网正常运行与事故扰动情况下的实时监测与分析计算,并及时获得并掌握电网运行的动态过程。PMU能够以数百Hz的速率采集电流、电压信息,通过计算获得测点的功率、相位、功角等信息,并以每秒几十帧的频率向主站发送。PMU通过全球定位系统(GPS)对时,能够保证全网数据的同步性,时标信息与数据同时存储并发送到主站。因此,WAMS能够实现在同一时间参考轴下获取大型互联系统实时动态信息和稳态信息,为电力系统区间动态监视、分析甚至控制提供前提条件,使电力系统监视从稳态阶段提高到动态阶段。
发明内容
为了更好地把握大规模风电的运行规律和特性,提高风电消纳能力和新能源运行管理水平,保障电网安全稳定运行,本发明提供一种基于WAMS数据的风电出力特性分析方法,对风电出力特性指标进行量化,从而能够为分析大规模风电接入对电网的影响提供技术参考。同时,通过对大量风电运行监测数据的积累分析,也有助于协助调度运行人员把握风电运行规律,为更加科学合理地消纳风电和保障电网安全稳定运行提供分析基础。
为了实现上述发明目的,本发明采取如下技术方案:
本发明提供一种基于WAMS数据的风电出力特性分析方法,所述方法包括以下步骤:
步骤1:确定风电出力特性;
步骤2:分析风电出力特性;
步骤3:定义风电场群,并对其进行统计分析。
所述步骤1中,风电出力特性主要包括风电出力极值特性、风电功率峰谷特性、风电功率短时波动特性、风电调峰贡献、日反调峰特性、偏差相关性、风电负荷率概率分布、风电穿透率和风电利用小时数。
所述步骤2中,对风电出力极值特性、风电功率峰谷特性、风电功率短时波动特性、风电调峰贡献、日反调峰特性、偏差相关性、风电负荷率概率分布、风电穿透率和风电利用小时数分别进行分析。
所述风电出力极值特性主要包括风电出力最大值、风电出力最小值及风电出力均值;
1)所述风电出力最大值主要包括日出力最大值、月出力最大值和任意时间段的出力最大值;
2)所述风电出力最小值主要包括日出力最小值、月出力最小值和任意时间段的出力最小值;
3)所述风电出力均值主要包括日平均出力、月平均出力和任意时间段的平均出力;其中,任意时间段的平均出力用Pavg表示,有:
其中,Pi为日有功出力,T为任意时间段。
所述风电功率峰谷特性主要包括日风电功率最大值、日风电功率最小值和日风电功率最大峰谷差;日风电功率最大峰谷差用Pmpv表示,有:
Pmpv=Pmax-Pmin (2)
其中,Pmax为日风电功率最大值,Pmin为日风电功率最小值。
所述风电功率短时波动特性主要包括风电功率短时波动量、风电功率最大上升量和风电功率下降速率;
所述风电功率短时波动量为间隔时段n内风电功率最大值与风电功率最小值之差,表示为:
ΔPn=Pn max-Pn min (3)
其中,ΔPn为风电功率短时波动量,Pn max为间隔时段n内风电功率最大值,Pn min为间隔时段n内风电功率最小值,n取1min、5min或15min。
所述风电调峰贡献为实际风电负荷峰谷差与净风电负荷峰谷差之差,净风电负荷等于实际风电负荷减去风电发电出力;
风电调峰贡献用ΔPpc表示,有:
ΔPpc=ΔPsjfh-ΔPjfh (4)
其中,ΔPsjfh为实际风电负荷峰谷差,ΔPjfh为净风电负荷峰谷差。
所述日反调峰特性是指风电出力峰谷与电网负荷峰谷不匹配,风电出力最大值出现在用电低谷期,或风电出力最小值出现在用电高峰期;
设Pi为日负荷序列,其中i=1,2,…,N,N表示负荷采样点总数;Pave为日负荷均值,Δk为负荷峰谷判断阈值倍数,取0.1;用电高峰期和用电低谷期分别用Tpeak和Tvalley表示,且满足Tpeak∈[Tp1,Tp2],Tvalley∈[Tv1,Tv2],其中,Tp1和Tp2分别为用电高峰期负荷段开始时刻和结束时刻,Tv1和Tv2分别为用电低谷期负荷段开始时刻和结束时刻;具体有:
1)当负荷曲线上连续3个负荷采样点的值大于等于Pave*(1+Δk)时,即:
Pi≥Pave*(1+Δk)i=k1,k1+1,k1+2 (5)
则负荷曲线上负荷采样点k1的值Pk1对应的时刻即为Tp1;其中,k1,k1+1,k1+2为连续3个负荷采样点;
当负荷曲线连续满足式(5)后,出现3个负荷采样点的值小于等于Pave*(1+Δk)时,即:
Pi≤Pave*(1+Δk)i=k2,k2+1,k2+2 (6)
则负荷曲线上负荷采样点k2的值Pk2对应的时刻即为Tp2;其中,k2,k2+1,k2+2为连续3个负荷采样点;
2)当负荷曲线上连续3个负荷采样点的值小于等于Pave*(1-Δk)时,即:
Pi≤Pave*(1-Δk)i=k3,k3+1,k3+2 (7)
则负荷曲线上负荷采样点k3的值Pk3对应的时刻即为Tv1;其中,k3,k3+1,k3+2为连续3个负荷采样点;
当负荷曲线连续满足式(7)后,出现3个负荷采样点的值大于等于Pave*(1-Δk)时,即:
Pi≥Pave*(1-Δk)i=k4,k4+1,k4+2 (8)
则负荷曲线上负荷采样点k4的值Pk4对应的时刻即为Tv2;其中,k4,k4+1,k4+2为连续3个负荷采样点;
风电出力高峰期和风电出力低谷期分别用Twpeak和Twvalley表示,同理,得到风电出力高峰期开始时刻Twp1和结束时刻Twp2,以及风电出力低谷期开始时刻Twv1和结束时刻Twv2;
分析用电高峰期、用电低谷期、风电出力高峰期和风电出力低谷期的一致性,若风电出力最大值出现在用电低谷期,或风电出力最小值出现在用电高峰期,则认为出现反调峰。
所述偏差相关性采用相关系数表示,假定数据序列X和Y分别表示风电出力与频率偏差,则两者的相关系数表示为:
其中,和分别表示风电出力均值与频率偏差均值;r表示风电出力与频率偏差的相关系数,取值在-1和+1之间,当r的绝对值愈接近1时,表明风电出力与频率偏差关系愈紧密;如果r接近于0,则表明风电出力与频率偏差关系不紧密。
所述风电负荷率概率分布是将风电负荷率划分为11个范围段,分别为5%及以下、5%~10%、10%~20%、20%~30%、30%~40%、40%~50%、50%~60%、60%~70%、70%~80%、80%~90%、90%~100%;风电负荷率为风电功率值与风电可用容量的比值,风电负荷率概率为落在各范围段的发电时间长度与总时间长度的比值,有:
Rrate=Woput/Luput×100% (10)
其中,Rrate为风电负荷率,Woput为风电功率值,Luput为风电可用容量。
所述风电穿透率为风电场装机容量占系统总负荷的比例,有:
Rwpp=Wcap/Lmax×100% (11)
其中,Rwpp为风电穿透率,Wcap为风电场装机容量,Lmax为系统总负荷。
所述风电利用小时数统计周期内风电发电量与风电场装机容量之比,有:
Hutime=Woput/Wcap (12)
其中,Hutime为风电利用小时数,Woput为风电发电量,Wcap为风电场装机容量。
所述步骤3包括以下步骤:
步骤3-1:根据风电场的地理特性和电气特性,将相关特征的风电场组合定义为风电场群;
步骤3-2:对风电场群进行功率和电量的统计、分析和查询,以及风电功率波动变化与电网频率变化的相关性分析,实现对风电场的分层、分区深度统计分析。
与现有技术相比,本发明的有益效果在于:
(1)从风电出力随机波动性、峰谷特性和调峰特性等方面出发,提供风电出力特性的分析方法,主要包括风电出力极值特性、风电功率峰谷特性、风电功率短时波动特性、风电调峰贡献、日反调峰特性、偏差相关性、风电负荷率概率分布、风电穿透率和风电利用小时数;
(2)通过整合电网风电实时数据与非实时数据,量化分析了风电场及用户自定义风电场群的运行特性,可对任意时间范围内的风电出力特性指标进行全面量化分析;有助于调度运行及研究人员对风电运行规律的把握,以及电网对风电的接纳能力和电网安全运行水平的提高,对大规模风电接入后电网运行具有指导作用。
附图说明
图1是本发明实施例中风电出力特性分析数据流程图;
图2是本发明实施例中风电出力特性示意图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
本发明通过分析风电出力特性,为运行人员全面把控风电场及其集群的出力特性提供技术依据和参考。通过分析风电出力随机波动性、峰谷特性和调峰特性,主要包括风电出力极值特性、风电功率最大峰谷差、风电功率短时波动、风电调峰贡献和反调峰、指标偏差相关性、风电负荷率概率分布、风电穿透率和风电利用小时等特性。通过对风电出力特性的分析,可以达到量化地描述风电出力特性和分析其对电网影响程度。根据风电场地理特性、电气特性对风电场进行组合后,实现对风电场集群出力特性的分层、分区深度挖掘。依据各风电出力特性的发展变化趋势,可准确判断大规模风电接入对电网的影响程度,能够为风电场运行管理人员提供决策支持。
进行风电场出力特性分析时,风电出力数据主要从D5000平台数据库获取,包括PMU数据和RTU数据等,获取到风电出力数据后在应用服务器上计算各指标值。风电出力特性分析数据流程如图1所示。
基于WAMS数据的风电出力特性分析方法包括以下步骤:
步骤1:确定风电出力特性;
步骤2:分析风电出力特性;
步骤3:定义风电场群,并对其进行统计分析。
所述步骤1中,(如图2)风电出力特性主要包括风电出力极值特性、风电功率峰谷特性、风电功率短时波动特性、风电调峰贡献、日反调峰特性、偏差相关性、风电负荷率概率分布、风电穿透率和风电利用小时数。
所述步骤2中,对风电出力极值特性、风电功率峰谷特性、风电功率短时波动特性、风电调峰贡献、日反调峰特性、偏差相关性、风电负荷率概率分布、风电穿透率和风电利用小时数分别进行分析。
所述风电出力极值特性主要包括风电出力最大值、风电出力最小值及风电出力均值;
1)所述风电出力最大值主要包括日出力最大值、月出力最大值和任意时间段的出力最大值;
2)所述风电出力最小值主要包括日出力最小值、月出力最小值和任意时间段的出力最小值;
3)所述风电出力均值主要包括日平均出力、月平均出力和任意时间段的平均出力;其中,任意时间段的平均出力用Pavg表示,有:
其中,Pi为日有功出力,T为任意时间段。
所述风电功率峰谷特性主要包括日风电功率最大值、日风电功率最小值和日风电功率最大峰谷差;日风电功率最大峰谷差用Pmpv表示,有:
Pmpv=Pmax-Pmin (2)
其中,Pmax为日风电功率最大值,Pmin为日风电功率最小值。
风电功率峰谷差可在一定程度上反映风电出力的波动特性,日风电功率最大峰谷差反映风电功率一天的最大波动量。
所述风电功率短时波动特性是影响系统调频设置的一个重要因素,当风电场容量在系统中所占的比例较大时,其输出功率的随机波动性对电网频率的影响会比较显著,将会影响到电网的电能质量和系统中其他一些频率敏感负荷的正常工作。同时,风电场出力频繁的波动和变化也会引起电网电压波动和闪变。因此,风电功率最大峰谷差和短时波动指标的分析与研究对于分析风力发电的波动性和随机性,以及风电对电网的影响非常重要。
所述风电功率短时波动特性主要包括风电功率短时波动量、风电功率最大上升量和风电功率下降速率;
所述风电功率短时波动量为间隔时段n内风电功率最大值与风电功率最小值之差,表示为:
ΔPn=Pn max-Pn min (3)
其中,ΔPn为风电功率短时波动量,Pn max为间隔时段n内风电功率最大值,Pn min为间隔时段n内风电功率最小值,n取1min、5min或15min。
调峰的目的就是要满足电网实际负荷的变化,通过比较实际负荷、风电功率、净负荷的变化曲线,可以分析风电对电网风电调峰贡献的大小,即风电出力的变化对常规机组调峰的影响。
所述风电调峰贡献为实际风电负荷峰谷差与净风电负荷峰谷差之差,净风电负荷等于实际风电负荷减去风电发电出力;净负荷曲线主要反映电网中除风电之外的其它常规机组的实际发电曲线。
风电调峰贡献用ΔPpc表示,有:
ΔPpc=ΔPsjfh-ΔPjfh (4)
其中,ΔPsjfh为实际风电负荷峰谷差,ΔPjfh为净风电负荷峰谷差。
所述日反调峰特性是指风电出力峰谷与电网负荷峰谷不匹配,风电出力最大值出现在用电低谷期,或风电出力最小值出现在用电高峰期。
由于社会生活和生产周期性、气象因素等影响,不同地区用电负荷峰谷出现的时间段存在一定差异,若以固定时间段作为负荷峰谷出现的时段必然存在问题。为准确确定负荷高峰期和低谷期,在此根据负荷曲线的实际变化动态计算其峰谷,再根据风电出力的峰谷准确地分析是否发生日内反调峰。
设Pi为日负荷序列,其中i=1,2,…,N,N表示负荷采样点总数;Pave为日负荷均值,Δk为负荷峰谷判断阈值倍数,取0.1;用电高峰期和用电低谷期分别用Tpeak和Tvalley表示,且满足Tpeak∈[Tp1,Tp2],Tvalley∈[Tv1,Tv2],其中,Tp1和Tp2分别为用电高峰期负荷段开始时刻和结束时刻,Tv1和Tv2分别为用电低谷期负荷段开始时刻和结束时刻;具体有:
1)当负荷曲线上连续3个负荷采样点的值大于等于Pave*(1+Δk)时,即:
Pi≥Pave*(1+Δk)i=k1,k1+1,k1+2 (5)
则负荷曲线上负荷采样点k1的值Pk1对应的时刻即为Tp1;其中,k1,k1+1,k1+2为连续3个负荷采样点;
当负荷曲线连续满足式(5)后,出现3个负荷采样点的值小于等于Pave*(1+Δk)时,即:
Pi≤Pave*(1+Δk)i=k2,k2+1,k2+2 (6)
则负荷曲线上负荷采样点k2的值Pk2对应的时刻即为Tp2;其中,k2,k2+1,k2+2为连续3个负荷采样点;
2)当负荷曲线上连续3个负荷采样点的值小于等于Pave*(1-Δk)时,即:
Pi≤Pave*(1-Δk)i=k3,k3+1,k3+2 (7)
则负荷曲线上负荷采样点k3的值Pk3对应的时刻即为Tv1;其中,k3,k3+1,k3+2为连续3个负荷采样点;
当负荷曲线连续满足式(7)后,出现3个负荷采样点的值大于等于Pave*(1-Δk)时,即:
Pi≥Pave*(1-Δk)i=k4,k4+1,k4+2 (8)
则负荷曲线上负荷采样点k4的值Pk4对应的时刻即为Tv2;其中,k4,k4+1,k4+2为连续3个负荷采样点;
风电出力高峰期和风电出力低谷期分别用Twpeak和Twvalley表示,同理,得到风电出力高峰期开始时刻Twp1和结束时刻Twp2,以及风电出力低谷期开始时刻Twv1和结束时刻Twv2;
分析用电高峰期、用电低谷期、风电出力高峰期和风电出力低谷期的一致性,若风电出力最大值出现在用电低谷期,或风电出力最小值出现在用电高峰期,则认为出现反调峰。
所述偏差相关性采用相关系数表示,假定数据序列X和Y分别表示风电出力与频率偏差,则两者的相关系数表示为:
其中,和分别表示风电出力均值与频率偏差均值;r表示风电出力与频率偏差的相关系数,取值在-1和+1之间,当r的绝对值愈接近1时,表明风电出力与频率偏差关系愈紧密;如果r接近于0,则表明风电出力与频率偏差关系不紧密。
所述风电负荷率概率分布是将风电负荷率划分为11个范围段,分别为5%及以下、5%~10%、10%~20%、20%~30%、30%~40%、40%~50%、50%~60%、60%~70%、70%~80%、80%~90%、90%~100%;风电负荷率为风电功率值与风电可用容量的比值,风电负荷率概率为落在各范围段的发电时间长度与总时间长度的比值,有:
Rrate=Woput/Luput×100% (10)
其中,Rrate为风电负荷率,Woput为风电功率值,Luput为风电可用容量。
所述风电穿透率为风电场装机容量占系统总负荷的比例,风电场穿透功率极限定义为接入系统的风电场最大装机容量与系统最大负荷的百分比。有:
Rwpp=Wcap/Lmax×100% (11)
其中,Rwpp为风电穿透率,Wcap为风电场装机容量,Lmax为系统总负荷。
所述风电利用小时数统计周期内风电发电量与风电场装机容量之比,利用小时数可以折算到年利用小时数,这个指标是对全年利用小时的估计,对风电产业有指导意义,可以加强对风电的管理。有:
Hutime=Woput/Wcap (12)
其中,Hutime为风电利用小时数,Woput为风电发电量,Wcap为风电场装机容量。
所述步骤3包括以下步骤:
步骤3-1:根据风电场的地理特性和电气特性,将相关特征的风电场组合定义为风电场群;
步骤3-2:对风电场群进行功率和电量的统计、分析和查询,以及风电功率波动变化与电网频率变化的相关性分析,实现对风电场的分层、分区深度统计分析。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,所属领域的普通技术人员参照上述实施例依然可以对本发明的具体实施方式进行修改或者等同替换,这些未脱离本发明精神和范围的任何修改或者等同替换,均在申请待批的本发明的权利要求保护范围之内。
Claims (13)
1.一种基于WAMS数据的风电出力特性分析方法,其特征在于:所述方法包括以下步骤:
步骤1:确定风电出力特性;
步骤2:分析风电出力特性;
步骤3:定义风电场群,并对其进行统计分析。
2.根据权利要求1所述的基于WAMS数据的风电出力特性分析方法,其特征在于:所述步骤1中,风电出力特性主要包括风电出力极值特性、风电功率峰谷特性、风电功率短时波动特性、风电调峰贡献、日反调峰特性、偏差相关性、风电负荷率概率分布、风电穿透率和风电利用小时数。
3.根据权利要求2所述的基于WAMS数据的风电出力特性分析方法,其特征在于:所述步骤2中,对风电出力极值特性、风电功率峰谷特性、风电功率短时波动特性、风电调峰贡献、日反调峰特性、偏差相关性、风电负荷率概率分布、风电穿透率和风电利用小时数分别进行分析。
4.根据权利要求3所述的基于WAMS数据的风电出力特性分析方法,其特征在于:所述风电出力极值特性主要包括风电出力最大值、风电出力最小值及风电出力均值;
1)所述风电出力最大值主要包括日出力最大值、月出力最大值和任意时间段的出力最大值;
2)所述风电出力最小值主要包括日出力最小值、月出力最小值和任意时间段的出力最小值;
3)所述风电出力均值主要包括日平均出力、月平均出力和任意时间段的平均出力;其中,任意时间段的平均出力用Pavg表示,有:
其中,Pi为日有功出力,T为任意时间段。
5.根据权利要求3所述的基于WAMS数据的风电出力特性分析方法,其特征在于:所述风电功率峰谷特性主要包括日风电功率最大值、日风电功率最小值和日风电功率最大峰谷差;日风电功率最大峰谷差用Pmpv表示,有:
Pmpv=Pmax-Pmin (2)
其中,Pmax为日风电功率最大值,Pmin为日风电功率最小值。
6.根据权利要求3所述的基于WAMS数据的风电出力特性分析方法,其特征在于:所述风电功率短时波动特性主要包括风电功率短时波动量、风电功率最大上升量和风电功率下降速率;
所述风电功率短时波动量为间隔时段n内风电功率最大值与风电功率最小值之差,表示为:
ΔPn=Pnmax-Pnmin (3)
其中,ΔPn为风电功率短时波动量,Pnmax为间隔时段n内风电功率最大值,Pnmin为间隔时段n内风电功率最小值,n取1min、5min或15min。
7.根据权利要求3所述的基于WAMS数据的风电出力特性分析方法,其特征在于:所述风电调峰贡献为实际风电负荷峰谷差与净风电负荷峰谷差之差,净风电负荷等于实际风电负荷减去风电发电出力;
风电调峰贡献用ΔPpc表示,有:
ΔPpc=ΔPsjfh-ΔPjfh (4)
其中,ΔPsjfh为实际风电负荷峰谷差,ΔPjfh为净风电负荷峰谷差。
8.根据权利要求3所述的基于WAMS数据的风电出力特性分析方法,其特征在于:所述日反调峰特性是指风电出力峰谷与电网负荷峰谷不匹配,风电出力最大值出现在用电低谷期,或风电出力最小值出现在用电高峰期;
设Pi为日负荷序列,其中i=1,2,…,N,N表示负荷采样点总数;Pave为日负荷均值,Δk为负荷峰谷判断阈值倍数,取0.1;用电高峰期和用电低谷期分别用Tpeak和Tvalley表示,且满足Tpeak∈[Tp1,Tp2],Tvalley∈[Tv1,Tv2],其中,Tp1和Tp2分别为用电高峰期负荷段开始时刻和结束时刻,Tv1和Tv2分别为用电低谷期负荷段开始时刻和结束时刻;具体有:
1)当负荷曲线上连续3个负荷采样点的值大于等于Pave*(1+Δk)时,即:
Pi≥Pave*(1+Δk) i=k1,k1+1,k1+2 (5)
则负荷曲线上负荷采样点k1的值Pk1对应的时刻即为Tp1;其中,k1,k1+1,k1+2为连续3个负荷采样点;
当负荷曲线连续满足式(5)后,出现3个负荷采样点的值小于等于Pave*(1+Δk)时,即:
Pi≤Pave*(1+Δk) i=k2,k2+1,k2+2 (6)
则负荷曲线上负荷采样点k2的值Pk2对应的时刻即为Tp2;其中,k2,k2+1,k2+2为连续3个负荷采样点;
2)当负荷曲线上连续3个负荷采样点的值小于等于Pave*(1-Δk)时,即:
Pi≤Pave*(1-Δk) i=k3,k3+1,k3+2 (7)
则负荷曲线上负荷采样点k3的值Pk3对应的时刻即为Tv1;其中,k3,k3+1,k3+2为连续3个负荷采样点;
当负荷曲线连续满足式(7)后,出现3个负荷采样点的值大于等于Pave*(1-Δk)时,即:
Pi≥Pave*(1-Δk) i=k4,k4+1,k4+2 (8)
则负荷曲线上负荷采样点k4的值Pk4对应的时刻即为Tv2;其中,k4,k4+1,k4+2为连续3个负荷采样点;
风电出力高峰期和风电出力低谷期分别用Twpeak和Twvalley表示,同理,得到风电出力高峰期开始时刻Twp1和结束时刻Twp2,以及风电出力低谷期开始时刻Twv1和结束时刻Twv2;
分析用电高峰期、用电低谷期、风电出力高峰期和风电出力低谷期的一致性,若风电出力最大值出现在用电低谷期,或风电出力最小值出现在用电高峰期,则认为出现反调峰。
9.根据权利要求3所述的基于WAMS数据的风电出力特性分析方法,其特征在于:所述偏差相关性采用相关系数表示,假定数据序列X和Y分别表示风电出力与频率偏差,则两者的相关系数表示为:
其中,和分别表示风电出力均值与频率偏差均值;r表示风电出力与频率偏差的相关系数,取值在-1和+1之间,当r的绝对值愈接近1时,表明风电出力与频率偏差关系愈紧密;如果r接近于0,则表明风电出力与频率偏差关系不紧密。
10.根据权利要求3所述的基于WAMS数据的风电出力特性分析方法,其特征在于:所述风电负荷率概率分布是将风电负荷率划分为11个范围段,分别为5%及以下、5%~10%、10%~20%、20%~30%、30%~40%、40%~50%、50%~60%、60%~70%、70%~80%、80%~90%、90%~100%;风电负荷率为风电功率值与风电可用容量的比值,风电负荷率概率为落在各范围段的发电时间长度与总时间长度的比值,有:
Rrate=Woput/Luput×100% (10)
其中,Rrate为风电负荷率,Woput为风电功率值,Luput为风电可用容量。
11.根据权利要求3所述的基于WAMS数据的风电出力特性分析方法,其特征在于:所述风电穿透率为风电场装机容量占系统总负荷的比例,有:
Rwpp=Wcap/Lmax×100% (11)
其中,Rwpp为风电穿透率,Wcap为风电场装机容量,Lmax为系统总负荷。
12.根据权利要求3所述的基于WAMS数据的风电出力特性分析方法,其特征在于:所述风电利用小时数统计周期内风电发电量与风电场装机容量之比,有:
Hutime=Woput/Wcap (12)
其中,Hutime为风电利用小时数,Woput为风电发电量,Wcap为风电场装机容量。
13.根据权利要求1所述的基于WAMS数据的风电出力特性分析方法,其特征在于:所述步骤3包括以下步骤:
步骤3-1:根据风电场的地理特性和电气特性,将相关特征的风电场组合定义为风电场群;
步骤3-2:对风电场群进行功率和电量的统计、分析和查询,以及风电功率波动变化与电网频率变化的相关性分析,实现对风电场的分层、分区深度统计分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510031408.XA CN104574216A (zh) | 2015-01-22 | 2015-01-22 | 一种基于wams数据的风电出力特性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510031408.XA CN104574216A (zh) | 2015-01-22 | 2015-01-22 | 一种基于wams数据的风电出力特性分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104574216A true CN104574216A (zh) | 2015-04-29 |
Family
ID=53090206
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510031408.XA Pending CN104574216A (zh) | 2015-01-22 | 2015-01-22 | 一种基于wams数据的风电出力特性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104574216A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112485589A (zh) * | 2020-11-17 | 2021-03-12 | 国网安徽省电力有限公司电力科学研究院 | 基于分布式实时数据仓库的台区停电及故障实时感知方法 |
CN112769156A (zh) * | 2020-12-28 | 2021-05-07 | 南昌大学 | 一种计及大规模海上风电并网的源网荷储协调运行方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100271222A1 (en) * | 2007-02-12 | 2010-10-28 | Shawn Kerrigan | Automatic system information determination of distributed renewable energy systems |
CN102496962A (zh) * | 2011-12-31 | 2012-06-13 | 清华大学 | 调峰调频约束下电力系统可消纳风电的辨识与控制方法 |
CN102682222A (zh) * | 2012-05-23 | 2012-09-19 | 甘肃省电力公司电力科学研究院 | 一种基于风电波动规律的连续潮流计算方法 |
CN103778340A (zh) * | 2014-01-27 | 2014-05-07 | 国家电网公司 | 大规模新能源发电特性的统计方法 |
CN103854066A (zh) * | 2012-12-05 | 2014-06-11 | 国网能源研究院 | 多能源电力远距离输送的数据优化处理方法及系统 |
-
2015
- 2015-01-22 CN CN201510031408.XA patent/CN104574216A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100271222A1 (en) * | 2007-02-12 | 2010-10-28 | Shawn Kerrigan | Automatic system information determination of distributed renewable energy systems |
CN102496962A (zh) * | 2011-12-31 | 2012-06-13 | 清华大学 | 调峰调频约束下电力系统可消纳风电的辨识与控制方法 |
CN102682222A (zh) * | 2012-05-23 | 2012-09-19 | 甘肃省电力公司电力科学研究院 | 一种基于风电波动规律的连续潮流计算方法 |
CN103854066A (zh) * | 2012-12-05 | 2014-06-11 | 国网能源研究院 | 多能源电力远距离输送的数据优化处理方法及系统 |
CN103778340A (zh) * | 2014-01-27 | 2014-05-07 | 国家电网公司 | 大规模新能源发电特性的统计方法 |
Non-Patent Citations (3)
Title |
---|
张圣楠: "华北电网风力发电数据挖掘分析系统研究与应用", 《华东电力》 * |
李剑楠 等: "多时空尺度风电统计特性评价指标体系及其应用", 《中国电机工程学报》 * |
胡媛媛: "风电运行相关性随机性及波动性分析", 《仪器仪表与分析监测》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112485589A (zh) * | 2020-11-17 | 2021-03-12 | 国网安徽省电力有限公司电力科学研究院 | 基于分布式实时数据仓库的台区停电及故障实时感知方法 |
CN112485589B (zh) * | 2020-11-17 | 2021-11-12 | 国网安徽省电力有限公司电力科学研究院 | 基于分布式实时数据仓库的台区停电及故障实时感知方法 |
CN112769156A (zh) * | 2020-12-28 | 2021-05-07 | 南昌大学 | 一种计及大规模海上风电并网的源网荷储协调运行方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104882905B (zh) | 一种考虑暂态安全约束的新能源接纳能力评估方法 | |
CN105787606A (zh) | 一种基于超短期负荷预测的电力调度在线趋势预警系统 | |
Terciyanlı et al. | Enhanced nationwide wind-electric power monitoring and forecast system | |
Sun et al. | Research on ultra-short-term wind power prediction considering source relevance | |
CN103679282B (zh) | 风电功率爬坡的预测方法 | |
CN103235981A (zh) | 一种风电电能质量趋势预测方法 | |
Mabel et al. | Adequacy evaluation of wind power generation systems | |
CN102831321A (zh) | 一种基于蒙特卡洛方法的风电场风险评估方法 | |
CN104037790B (zh) | 一种基于时序蒙特卡罗模拟的新能源接纳能力评估方法 | |
Shang et al. | Equivalent modeling of active distribution network considering the spatial uncertainty of renewable energy resources | |
CN103065202A (zh) | 一种基于组合核函数的风电场超短期风速预测方法 | |
CN104463356A (zh) | 一种基于多维信息人工神经网络算法的光伏发电功率预测方法 | |
Meenakshi et al. | Power Quality Monitoring of a Photovoltaic System Through IoT | |
CN104156889A (zh) | 一种基于wams数据的风电场性能评估系统及其评估方法 | |
CN106786608B (zh) | 一种适用于分布式电源接入的不确定谐波潮流计算方法 | |
CN109462223A (zh) | 一种基于大数据的电能质量与可靠性分析方法 | |
Kaplan et al. | A novel method based on Weibull distribution for short-term wind speed prediction | |
CN101819243B (zh) | 基于广域信息的静态电压失稳预测方法 | |
Deceglie et al. | Quantifying year-to-year variations in solar panel soiling from PV energy-production data | |
Dong et al. | Improving the accuracy of wind speed statistical analysis and wind energy utilization in the Ningxia Autonomous Region, China | |
CN116266331A (zh) | 电网公司区域内的电力能源结构监测方法及装置 | |
CN104574216A (zh) | 一种基于wams数据的风电出力特性分析方法 | |
Eltohamy et al. | Technical investigation for power system flexibility | |
CN104657787B (zh) | 一种风功率时间序列组合预测方法 | |
Hosseini Dolatabadi et al. | Enhancing photovoltaic farm capacity estimation: A comprehensive analysis with a novel approach |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150429 |
|
RJ01 | Rejection of invention patent application after publication |