CN107026460B - 基于新能源大规模并网下频率偏差系数的确定方法 - Google Patents

基于新能源大规模并网下频率偏差系数的确定方法 Download PDF

Info

Publication number
CN107026460B
CN107026460B CN201710289167.8A CN201710289167A CN107026460B CN 107026460 B CN107026460 B CN 107026460B CN 201710289167 A CN201710289167 A CN 201710289167A CN 107026460 B CN107026460 B CN 107026460B
Authority
CN
China
Prior art keywords
load
frequency deviation
interconnected
variance
power
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.)
Expired - Fee Related
Application number
CN201710289167.8A
Other languages
English (en)
Other versions
CN107026460A (zh
Inventor
冯劲松
李卫东
刘娆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201710289167.8A priority Critical patent/CN107026460B/zh
Publication of CN107026460A publication Critical patent/CN107026460A/zh
Application granted granted Critical
Publication of CN107026460B publication Critical patent/CN107026460B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06Electricity, gas or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Abstract

本发明公开了新能源大规模并网下频率偏差系数的确定方法,包括以下步骤:步骤1:采集控制区的负荷统计特性和风电出力特性;步骤2:根据CPS标准推导计算出电网的二次调频责任;步骤3:计算各区域的负荷和风电出力在分钟级别的统计方差,按照方差分布比例来重新分配各区域的频率偏差系数;步骤4:通过仿真比较改进之后的各控制区域的CPS标准值情况,得到新的频率偏差系数分配方法。本发明公开方法中所提供的新的频率偏差系数分配方法能能公平、合理的分配各控制区域的调频责任,有利于系统在现有基础上最大限度的消纳风电。

Description

基于新能源大规模并网下频率偏差系数的确定方法
技术领域
本发明属于电力系统二次调频技术领域,涉及一种基于新能源大规模并网下频率偏差系数的确定方法。涉及专利分类号H02发电、变电或配电H02J供电或配电的电路装置或系统;电能存储系统H02J3/00交流干线或交流配电网络的电路装置。
背景技术
随着环境问题的严峻以及可再生能源的发展,可再生能源将大规模地并入到电力系统之中,这对电力系统的稳定运行是一个挑战。在保证电网安全运行的前提下,为了能够更多的消纳可再生能源,提出了区域合并消纳风电的策略。在该策略下如何合理公平的分配确定区域的调频责任,即频率偏差系数,成为了一个难题。
而我国目前普遍实行的是按照区域规模来分配一个固定频率偏差系数,其中用该区域的年负荷最大功率来衡量区域规模。但这种方式无法时刻保证区域的频率偏差系数始终处于合理的区间,对电网的频率控制不利。另外,这种方法有着两个前提,即各区域日负荷曲线形态相似和各区负荷扰动分量构成相似。在风电并网之后,由于风电出力的波动性和随机性,不可避免的对接入风电区域的负荷波动产生影响。所以在可再生能源大规模并网下,按照以往的频率偏差系数的确定方法有着诸多不利之处。
发明内容
根据现有频率偏差系数确定方法的问题,本发明根据CPS标准中对区域调频效果的要求,提出按照分钟级负荷波动方差的大小来分配频率偏差系数的方法。具体的技术方案是:一种新能源大规模并网下频率偏差系数的确定方法,包括以下步骤:
步骤1:采集至少包括当前控制区域的负荷统计特性和风电出力特性的电网数据;根据负荷、风电曲线,得到此区域波动特性曲线;
步骤2:依据控制性能评价标准CPS(Control Performance Standard),计算出由多个控制区域互联之后形成的联合区域的二次调频责任Bs
步骤3:计算多个控制区域的负荷统计特性和风电出力特性在分钟级别的统计均方差Si,按照方差分布比例,重新分配各区域的频率偏差系数Bi
作为优选的实施方式,所述依据CPS标准,计算出电网的二次调频责任Bs过程如下:
—叠加多个所述控制区域的区域控制偏差ACEi,形成联合区域的总区域控制偏差ACES,忽略联络线交换功率在计划和测量中的误差和错误,将所述总区域控制偏差表示为所述二次调频责任和△f的乘积;
—根据互联区域中有功出力应尽量接近整个电网的负荷扰动量ΔPLs的原则,得出如下公式:
Figure GDA0002583519260000021
根据均方差公式RMS2(X)=σ2(X)+AVG2(X)=AVG(X2),对CPS标准中,作为指标要求的,Δf一分钟长期统计特性小于一目标值作为对频率控制,进行变换,得到如下表达式
Figure GDA0002583519260000022
将所述二次调频责任Bs代入上式,得到
Figure GDA0002583519260000023
化简并整理得到
Figure GDA0002583519260000024
—根据电网规定的ε1值和电网ΔPLs的历史数据,根据上述算出互联电力系统总的二次调频责任Bs
其中:Δf1为互联系统Δf1min的平均值,ε1为互联电力系统在一年的时间段内频差的1min平均值的均方根值,ACES为互联系统总的区域控制偏差;ACEi为区域i的区域控制偏差;ΔPtiei为区域i的联络线实际交换功率偏差;Bi为区域的频率偏差系数;Bs为互联电力系统总的频率偏差系数;ΔPLs为系统负荷扰动;
Figure GDA0002583519260000025
为互联系统一分钟ACE的平均值;
Figure GDA0002583519260000026
为互联系统一分钟负荷扰动的平均值。
作为优选的实施方式,所述步骤3具体过程如下:
根据方差公式:
Figure GDA0002583519260000031
其中S为样本方差;xi为样本数据点;M为样本平均值;n为样本数量;
根据滚动平均法分离出负荷和风电在分钟级别的负荷曲线,每分钟取500 个点计算出负荷及风电的方差;
Figure GDA0002583519260000032
式中:LFt是经滚动平均法后t时刻的负荷峰值;2M是样本数量;
分钟级别负荷分量的幅值LM为Lt及其对应的LFt的差值,将得到的LM数据代入式方差公式中即可得到方差;
互联电力系统区域i的方差为Si,则风电大规模并网后的频率偏差系数为:
Figure GDA0002583519260000033
本发明公开了新能源大规模并网下频率偏差系数的确定方法,包括以下步骤:步骤1:采集控制区的负荷统计特性和风电出力特性;步骤2:根据CPS标准推导计算出电网的二次调频责任;步骤3:计算各区域的负荷和风电出力在分钟级别的统计方差,按照方差分布比例来重新分配各区域的频率偏差系数;步骤4:通过仿真比较改进之后的各控制区域的CPS标准值情况,得到新的频率偏差系数分配方法。本发明公开方法中所提供的新的频率偏差系数分配方法能能公平、合理的分配各控制区域的调频责任,有利于系统在现有基础上最大限度的消纳风电。
附图说明
为了更清楚的说明本发明的实施例或现有技术的技术方案,下面将对实施例或现有技术描述中所需要使用的附图做一简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明频率偏差系数分配算法图;
图2为本发明单区域系统发电机负荷构成模型图;
图3为本发明CPS标准计算模型图;
图4为本发明的算法流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明的技术方案进行清楚、完整地描述:
如图1-4所示:一种新能源大规模并网下频率偏差系数的确定方法,包括以下步骤:
步骤1:采集至少包括当前控制区域的负荷统计特性和风电出力特性的电网数据;根据负荷、风电曲线,得到此区域波动特性曲线;
步骤2:依据控制性能评价标准(Control Performance Standard—CPS),计算出当前联合区域(多个区域互联之后的大区域)二次调频责任Bs
步骤3:计算多个控制区域的负荷统计特性和风电出力特性在分钟级别的统计均方差Si,按照方差分布比例,重新分配各区域的频率偏差系数Bi
所述依据CPS标准,计算出电网的二次调频责任Bs过程如下:
根据我国电网采用的对Δf的一分钟长期统计特性小于一目标值作为对频率控制的指标要求:
Figure GDA0002583519260000041
如果将互联区域中所有区域的ACEi相加,得到区域总的ACEs为:
ACEs=∑ACEi=∑ΔPtiei+(∑Bi)Δf (2.2)
忽略联络线交换功率在计划和测量中的误差和错误,∑ΔPtiei=0,
公式(2.2)化为:
ACEs≈Bs·Δf (2.3)
为了达到规定的频率控制目标,互联区域中有功出力应尽量接近整个电网的负荷扰动量ΔPLs,即ACEs≈-ΔPLs,所以:
Figure GDA0002583519260000042
根据均方差公式:
RMS2(X)=σ2(X)+AVG2(X)=AVG(X2) (2.5)
根据式(2.5)化简式(2.1)可得:
Figure GDA0002583519260000051
将固定值Bs代入得
Figure GDA0002583519260000052
根据式(2.1)得:
Figure GDA0002583519260000053
最终整理得到:
Figure GDA0002583519260000054
—根据电网规定的ε1值和电网ΔPLs的历史数据,按照公式(2.9)即可算出互联电力系统总的二次调频责任Bs
其中:Δf1为互联系统Δf1min的平均值,ε1为互联电力系统在一年的时间段内频差的1min平均值的均方根值,ACES为互联系统总的区域控制偏差;ACEi为区域i的区域控制偏差;ΔPtiei为区域i的联络线实际交换功率偏差;Bi为区域的频率偏差系数;Bs为互联电力系统总的频率偏差系数;ΔPLs为系统负荷扰动;
Figure GDA0002583519260000055
为互联系统一分钟ACE的平均值;
Figure GDA0002583519260000056
为互联系统一分钟负荷扰动的平均值。
所述步骤3具体过程如下:
根据方差公式:
Figure GDA0002583519260000057
其中S为样本方差;xi为样本数据点;M为样本平均值;n为样本数量;
根据滚动平均法分离出负荷和风电在分钟级别的负荷曲线,每分钟取500 个点计算出负荷及风电的方差;
Figure GDA0002583519260000061
式中:LFt是经滚动平均法后t时刻的负荷峰值;2M是样本数量;
分钟级别负荷分量的幅值LM为Lt及其对应的LFt的差值,将得到的LM数据代入式(3.1)中即可得到方差;
互联电力系统区域i的方差为Si,则风电大规模并网后的频率偏差系数为:
Figure GDA0002583519260000062
实施例:
以两区域电力系统仿真结果,说明本发明的实施过程,通过仿真比较改进之后的各控制区域的CPS标准值情况,得到新的频率偏差系数分配方法。
仿真模拟典型日负荷、风电曲线,利用两区域互联电力系统仿真模型,得到24小时两区域的CPS评价标准值。其中,区域1是采用按照区域规模分配B 系数的方法,区域2是利用负荷波动方差来分配B系数的方法。
表一 两种方法的CPS标准值
区域1 区域2
CPS1 178.5 195.8
CPS2 90 100
本文所提出的在新能源大规模并网下的频率偏差系数确定方法,符合新能源并网形势下的调频责任分配原则,较原有的按照区域规模分配B系数的方法,区域的CPS考核标准值更满足要求,同时也更公平。
以上所述,仅为本发明的一种具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。

Claims (2)

1.一种新能源大规模并网下频率偏差系数的确定方法,其特征在于包括以下步骤:
步骤1:采集至少包括当前控制区域的负荷统计特性和风电出力特性的电网数据;根据负荷、风电曲线,得到此区域波动特性曲线;
步骤2:依据控制性能评价标准CPS,计算出由多个控制区域互联之后形成的联合区域的二次调频责任Bs
步骤3:计算多个控制区域的负荷统计特性和风电出力特性在分钟级别的统计均方差Si,按照方差分布比例,重新分配各区域的频率偏差系数Bi;所述依据CPS标准,计算出电网的二次调频责任Bs过程如下:
—叠加多个所述控制区域的区域控制偏差ACEi,形成联合区域的总区域控制偏差ACES,忽略联络线交换功率在计划和测量中的误差和错误,将所述总区域控制偏差表示为所述二次调频责任和频率偏差△f的乘积;
—根据互联区域中有功出力应尽量接近整个电网的负荷扰动量ΔPLs的原则,得出如下公式:
Figure FDA0002561171020000011
根据均方差公式RMS2(X)=σ2(X)+AVG2(X)=AVG(X2),对CPS标准中,作为指标要求的,Δf一分钟长期统计特性小于一目标值作为对频率控制,进行变换,得到如下表达式
Figure FDA0002561171020000012
将所述二次调频责任Bs带入上式,得到
Figure FDA0002561171020000013
化简并整理得到
Figure FDA0002561171020000014
—根据电网规定的ε1值和电网ΔPLs的历史数据,根据上述算出互联电力系统总的二次调频责任Bs
其中:Δf1为互联系统Δf 1min的平均值,ε1为互联电力系统在一年的时间段内频差的1min平均值的均方根值,ACES为互联系统总的区域控制偏差;ACEi为区域i的区域控制偏差;ΔPtiei为区域i的联络线实际交换功率偏差;Bi为区域的频率偏差系数;Bs为互联电力系统总的频率偏差系数;ΔPLs为系统负荷扰动;
Figure FDA0002561171020000021
为互联系统一分钟ACE的平均值;
Figure FDA0002561171020000022
为互联系统一分钟负荷扰动的平均值。
2.根据权利要求1所述的新能源大规模并网下频率偏差系数的确定方法,其特征在于所述步骤3具体过程如下:
根据方差公式:
Figure FDA0002561171020000023
其中S为样本方差;xi为样本数据点;M为样本平均值;n为样本数量;
根据滚动平均法分离出负荷和风电在分钟级别的负荷曲线,每分钟取500个点计算出负荷及风电的方差;
Figure FDA0002561171020000024
式中:LFt是经滚动平均法后t时刻的负荷峰值;2M是样本数量;
分钟级别负荷分量的幅值LM为Lt及其对应的LFt的差值,将得到的LM数据代入式方差公式中即可得到方差;
互联电力系统区域i的方差为Si,则风电大规模并网后的频率偏差系数为:
Figure FDA0002561171020000025
CN201710289167.8A 2017-04-27 2017-04-27 基于新能源大规模并网下频率偏差系数的确定方法 Expired - Fee Related CN107026460B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710289167.8A CN107026460B (zh) 2017-04-27 2017-04-27 基于新能源大规模并网下频率偏差系数的确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710289167.8A CN107026460B (zh) 2017-04-27 2017-04-27 基于新能源大规模并网下频率偏差系数的确定方法

Publications (2)

Publication Number Publication Date
CN107026460A CN107026460A (zh) 2017-08-08
CN107026460B true CN107026460B (zh) 2020-09-29

Family

ID=59527348

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710289167.8A Expired - Fee Related CN107026460B (zh) 2017-04-27 2017-04-27 基于新能源大规模并网下频率偏差系数的确定方法

Country Status (1)

Country Link
CN (1) CN107026460B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108258722B (zh) * 2018-01-09 2022-08-09 国网辽宁省电力有限公司 一种提升电网新能源接纳能力的功频技术
CN108777481B (zh) * 2018-05-25 2020-09-08 国家电网公司西北分部 一种电网控制区频率偏差系数设置方法
CN108879721B (zh) * 2018-06-15 2021-11-02 国网河南省电力公司电力科学研究院 一种基于风电消纳的频率控制方法
CN109245090B (zh) * 2018-09-06 2021-04-20 大连理工大学 一种大功率缺失下频率最低点预测的解析模型的建模方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105207253A (zh) * 2015-09-15 2015-12-30 重庆大学 考虑风电及频率不确定性的agc随机动态优化调度方法
CN106208039A (zh) * 2016-02-25 2016-12-07 大连理工大学 一种源网荷互动运行控制性能评价标准评价方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105207253A (zh) * 2015-09-15 2015-12-30 重庆大学 考虑风电及频率不确定性的agc随机动态优化调度方法
CN106208039A (zh) * 2016-02-25 2016-12-07 大连理工大学 一种源网荷互动运行控制性能评价标准评价方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Controlling of Frequency Deviations in Interconnected Power Systems using Smart Techniques;Komal Wadhwa等;《1st IEEE International Conference on Power Electronics, Intelligent Control and Energy Systems (ICPEICES-2016)》;20161231;第1-6页 *

Also Published As

Publication number Publication date
CN107026460A (zh) 2017-08-08

Similar Documents

Publication Publication Date Title
CN107026460B (zh) 基于新能源大规模并网下频率偏差系数的确定方法
KR102304537B1 (ko) 송전 주파수 응답
Gevorgian Highly accurate method for real-time active power reserve estimation for utility-scale photovoltaic power plants
CN104778367B (zh) 基于单一状态断面的广域戴维南等值参数在线计算方法
KR20140105506A (ko) 전력 생산 용량의 적응 및 에너지 저장 유닛 사이즈의 결정
Jin et al. Equivalent modeling of wind energy conversion considering overall effect of pitch angle controllers in wind farm
CN102590685B (zh) 一种配电网电流匹配状态估计方法
CN104680253B (zh) 基于可靠性与经济性相协调的配电网规划投资优化方法
CN110299762B (zh) 基于pmu准实时数据的主动配电网抗差估计方法
CN110543696B (zh) 一种未建模小机组参与电力市场出清和安全校核的方法
CN106786608B (zh) 一种适用于分布式电源接入的不确定谐波潮流计算方法
CN107516902B (zh) 计及市场运作规则的多种可中断负荷实时协调控制方法
CN103311925B (zh) 一种复杂电网自组织临界状态的辨识方法
CN103927597A (zh) 基于自回归滑动平均模型的风电功率超短期预测方法
CN109636033B (zh) 一种分布式光伏项目自发自用综合电价预测方法
CN109149566A (zh) 一种大功率缺失下频率最低点预测的仿真模型的建模方法
Ding et al. A dynamic period partition method for time-of-use pricing with high-penetration renewable energy
Eltohamy et al. Technical investigation for power system flexibility
Tazi et al. Demand and energy management in smart grid: Techniques and implementation
CN109245090A (zh) 一种大功率缺失下频率最低点预测的解析模型的建模方法
CN105281373A (zh) 不同风电出力水平下电网弃风电量计算方法和系统
Lave et al. Full-scale demonstration of distribution system parameter estimation to improve low-voltage circuit models
Eltohamy et al. Power system flexibility metrics review with high penetration of variable renewable generation
Wu et al. Short-term wind power prediction based on time series analysis model
Ba et al. Power energy output prediction of small wind urban for decision making

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200929

Termination date: 20210427

CF01 Termination of patent right due to non-payment of annual fee