CN110224393B - 一种基于最小切负荷模型的新能源消纳评估方法 - Google Patents

一种基于最小切负荷模型的新能源消纳评估方法 Download PDF

Info

Publication number
CN110224393B
CN110224393B CN201910438343.9A CN201910438343A CN110224393B CN 110224393 B CN110224393 B CN 110224393B CN 201910438343 A CN201910438343 A CN 201910438343A CN 110224393 B CN110224393 B CN 110224393B
Authority
CN
China
Prior art keywords
typical
power
node
new energy
scene
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
Application number
CN201910438343.9A
Other languages
English (en)
Other versions
CN110224393A (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.)
Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd
Yangjiang Power Supply Bureau of Guangdong Power Grid Co Ltd
Original Assignee
Chongqing University
Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd
Yangjiang Power Supply Bureau of Guangdong Power Grid Co Ltd
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 Chongqing University, Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd, Yangjiang Power Supply Bureau of Guangdong Power Grid Co Ltd filed Critical Chongqing University
Priority to CN201910438343.9A priority Critical patent/CN110224393B/zh
Publication of CN110224393A publication Critical patent/CN110224393A/zh
Application granted granted Critical
Publication of CN110224393B publication Critical patent/CN110224393B/zh
Active 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
    • 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/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48Controlling the sharing of the in-phase component
    • 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]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开一种基于最小切负荷模型的新能源消纳评估方法,主要步骤为:1)获取T年内电力系统数据,并建立电力系统N‑1预想故障集。2)建立电力系统运行典型场景。3)对电力系统的水电站进行优化调度,使电力系统在正常运行状态下满足正常运行约束条件。4)将电力系统数据输入最小切负荷模型中,得到最小负荷削减量的期望值,使电力系统在N‑1预想故障状态下满足故障状态约束条件。5)计算得到新能源渗透率上限和新能源消纳能力。6)建立新能源渗透率和最小负荷削减量期望值的关系曲线。本发明适用于水电富集地区,本发明充分考虑了系统中的电源灵活性,提高新能源消纳水平,准确评估了系统的新能源消纳能力。

Description

一种基于最小切负荷模型的新能源消纳评估方法
技术领域
本发明涉及新能源电力系统技术领域,具体是一种基于最小切负荷模型的新能源消纳评估方法。
背景技术
近年来,新能源发电已经显示出巨大的环境效益和经济效益,我国新能源迅猛发展。在我国南方水资源丰富地区,建立了一定规模的蓄水式水电站,使得系统具有较高的灵活性,可以充分利用系统中的电源灵活性来实现一定规模新能源的全额消纳。然而,新能源的消纳是指在保证系统安全稳定运行的前提下,合理利用各种资源使得系统能够高效利用新能源。可见,新能源消纳仍是以系统的安全稳定运行为大前提,系统在消纳新能源时,不应产生大量的切负荷、线路过载等不良运行状态。但新能源电源的并网运行改变了电力系统的电源结构,网络的潮流分布也显著改变,使得电力系统可能面临着线路过载和切负荷的风险。
电力系统中电源的灵活性,能够有效提升新能源消纳能力,特别是在我国南方水电富集的地区,利用灵活性电源蓄水式水电站的调节作用可以实现区域内新能源的全额消纳。然而,现有的大多数关于新能源消纳评估的研究简化了水电站的运行模型,蓄水式水电站的水库可能承担着防洪、灌溉和供水等其它综合性任务,研究中没有考虑水库因其它任务需求而限制机组的发电能力,因此无法准确评估系统的新能源消纳能力。此外,目前大多数研究只考虑了系统在正常运行状态下,并没有对系统进行N-1安全校验。在实际运行中,新能源接入的系统在N-1预想故障状态可能出现切负荷现象,影响供电的可靠性。目前缺乏一种综合考虑系统的正常运行状态和N-1预想故障状态进行新能源消纳评估的方法。
发明内容
本发明的目的是解决现有技术中存在的问题。
为实现本发明目的而采用的技术方案是这样的,一种基于最小切负荷模型的新能源消纳评估方法,主要包括以下步骤:
1)获取T年内电力系统数据,并建立电力系统N-1预想故障集。
2)利用K均值聚类方法对电力系统数据进行处理,建立电力系统运行典型场景,并分别计算每个典型场景发生的概率。
建立电力系统运行典型场景的主要步骤如下:
2.1)基于T年内电力系统数据,选出T年内电力系统在每个季节的典型日,并设置典型日的功率采样点和采样时刻。
在每个采样时刻,对所有典型日的功率采样点进行功率数据采样,建立原始数据集,并将典型日划分为高峰、平衡和低谷三个时段。
2.2)根据采样点、采样时刻和原始数据集,建立典型日内高峰、平衡和低谷三个时段内的典型场景。典型日内的典型场景总数N=Ttop+Tlev+Tlow
每个典型场景表示如下:
Figure BDA0002071261830000021
式中,
Figure BDA0002071261830000022
Figure BDA0002071261830000023
分别为高峰、平衡和低谷时段内的第t个采样时刻测得的源荷功率。x=1,2,…,m,m为采样点总数。Ttop、Tlev和Tlow分别为高峰、平衡和低谷时段内的采样时刻总数。
2.3)利用K均值法对典型场景进行聚类,主要步骤如下:
2.3.1)根据典型日高峰、平衡和低谷三个时段内的典型场景,确定高峰、平衡和低谷三个时段内的聚类个数Ktop、Klev和Klow。利用密度参数法优化生成Ktop、Klev和Klow个初始聚类中心。
2.3.2)分别计算每个典型场景和所有聚类中心的欧式距离,并将所述典型场景分配到距离最近的类中。利用均值法迭代更新聚类中心直至收敛,从而在高峰时段建立Ktop个类和Ktop个聚类中心,在平衡时段建立Klev个类和Klev个聚类中心,在低谷时段建立Klow个类和Klow个聚类中心。每个类中的典型场景数量即为每个聚类中心表示的典型场景数量。
2.4)对每个聚类中心和所述聚类中心所在类中的所有典型场景中的源荷功率数据进行比较,找出与聚类中心差距最小的典型场景,并利用与聚类中心差距最小的典型场景替换聚类中心。每个替换后的聚类中心即为典型日的典型场景。
2.5)第s个典型场景发生的概率Ps如下所示:
Figure BDA0002071261830000024
式中,s为典型场景的编号。Ns为典型场景的总数。ns为典型场景s所表示的典型场景数量。N为典型日中的典型场景总数。Ts为典型场景s所属的典型日在一年中所代表的天数。
3)对电力系统的水电站进行优化调度,使电力系统在正常运行状态下满足正常运行约束条件。
所述正常运行约束条件主要包括风电与光伏发电的消纳约束、水电的消纳约束、正常运行功率平衡约束、线路传输容量约束、机组出力约束、水电站库容约束、水电站发电流量约束和水电站水量平衡约束。
风电与光伏发电的消纳约束如下所示:
Figure BDA0002071261830000025
式中,Nw和NPV分别为电力系统并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。
Figure BDA0002071261830000031
Figure BDA0002071261830000032
分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的出力情况。
Figure BDA0002071261830000033
Figure BDA0002071261830000034
分别为第s个典型场景下风电场和光伏电站总的最大出力。
水电的消纳约束如下所示:
Figure BDA0002071261830000035
式中,Nstop、Nslev和Nslow分别为典型日内高峰平衡低谷时段内的典型场景总数。Nh为并网的水电站总数。sh为水电站的编号。
Figure BDA0002071261830000036
Figure BDA0002071261830000037
分别为高峰平衡低谷各个时段内的第s个典型场景下的第sh座水电站的出力情况。
Figure BDA0002071261830000038
Figure BDA0002071261830000039
Figure BDA00020712618300000310
分别为高峰平衡低谷各个时段内的第s个典型场景所持续的时间。
Figure BDA00020712618300000311
为水电站在任意典型日的最大发电量。
正常运行功率平衡约束如下所示:
其中,
Figure BDA00020712618300000313
Figure BDA00020712618300000314
分别为节点i在第s个典型场景下的火电、光伏、风电和水电的有功注入功率。
Figure BDA00020712618300000315
Figure BDA00020712618300000316
分别为节点i在第s个典型场景下的有功负荷和注入功率。
Figure BDA00020712618300000317
Figure BDA00020712618300000318
分别为节点i与节点j在第s个典型场景下的电压相角。
Figure BDA00020712618300000319
Figure BDA00020712618300000320
分别是以
Figure BDA00020712618300000321
为支路导纳建立起来的节点导纳矩阵的自导纳和互导纳。
正常运行时第s个典型场景的节点i的自导纳
Figure BDA00020712618300000322
如下所示:
Figure BDA00020712618300000323
正常运行时第s个典型场景的节点i和节点j互导纳
Figure BDA00020712618300000324
如下所示:
Figure BDA00020712618300000325
线路传输容量约束如下所示:
Figure BDA0002071261830000041
式中,Pij.max为支路ij之间的传输极限。
Figure BDA0002071261830000042
为支路ij在第s个典型场景下所传输的有功功率。
支路ij在第s个典型场景下所传输的有功功率
Figure BDA0002071261830000043
如下所示:
Figure BDA0002071261830000044
机组出力约束分别如公式10至公式13所示:
Figure BDA0002071261830000045
Figure BDA0002071261830000046
Figure BDA0002071261830000047
Figure BDA0002071261830000048
式中,
Figure BDA0002071261830000049
Figure BDA00020712618300000410
分别为节点i所连火电机组在第s个典型场景下的有功出力上下限。
Figure BDA00020712618300000411
Figure BDA00020712618300000412
分别为节点i所连光伏电站在第s个典型场景下的有功出力上下限。
Figure BDA00020712618300000413
Figure BDA00020712618300000414
分别为节点i所连风电场在第s个典型场景下的有功出力上下限。
Figure BDA00020712618300000415
Figure BDA00020712618300000416
分别为节点i所连接的水电站在第s个典型场景时的有功出力上下限。
水电站发电流量约束如下所示:
Figure BDA00020712618300000417
式中,
Figure BDA00020712618300000418
Figure BDA00020712618300000419
分别为节点i所连接的水电站在第s个典型场景时所分配的发电流量上下限。
发电流量
Figure BDA00020712618300000420
和水电出力
Figure BDA00020712618300000421
的转化关系如下所示:
Figure BDA00020712618300000422
式中,ηi为节点i所连接的水电站的发电效率。
Figure BDA00020712618300000423
为节点i所连接的水电站在第s个典型场景下的发电流量。
Figure BDA00020712618300000424
为节点i所连接的水电站在第s个典型场景下的发电水头平均高度。
水电站水量平衡约束如下所示:
Figure BDA0002071261830000051
其中,
Figure BDA0002071261830000052
Figure BDA0002071261830000053
分别为节点i所连接的水电站在第s个典型场景初时和末时的水库蓄水量。
Figure BDA0002071261830000054
Figure BDA0002071261830000055
分别为节点i所连接的水电站在第s个典型场景时的入库流量、发电流量和其它流量。所述其它流量主要包括灌溉流量和生活用水流量。
水电站库容约束如下所示:
Figure BDA0002071261830000056
式中,
Figure BDA0002071261830000057
Figure BDA0002071261830000058
分别为节点i所连接的水电站在第s个典型场景时的蓄水量上下限。
Figure BDA0002071261830000059
为节点i所连接的水电站在第s个典型场景时的蓄水量。
4)基于N-1预想故障集,建立最小切负荷模型。将电力系统数据输入最小切负荷模型中,得到最小负荷削减量的期望值,使电力系统在N-1预想故障状态下满足故障状态约束条件。
所述故障状态约束条件主要包括风电和光伏发电的消纳约束、水电的消纳约束、故障状态功率平衡约束、线路传输容量约束、机组出力约束、故障状态节点负荷削减量约束、水电站库容约束、水电站发电流量约束、水电站水量平衡约束。
故障状态功率平衡约束如下所示:
Figure BDA00020712618300000510
式中,
Figure BDA00020712618300000511
Figure BDA00020712618300000512
分别为节点i在第f个故障状态下的火电、光伏、风电和水电机组的有功注入功率。
Figure BDA00020712618300000513
为节点i在第f个故障状态下的有功负荷。
Figure BDA00020712618300000514
为节点i在第f个故障状态下的有功负荷削减量。
Figure BDA00020712618300000515
Figure BDA00020712618300000516
分别为节点i与节点j在第f个故障状态下的电压相角。
Figure BDA00020712618300000517
Figure BDA00020712618300000518
分别是在第f个故障状态下以
Figure BDA00020712618300000519
为支路导纳建立起来的节点导纳矩阵的自导纳及互导纳。
第f个故障状态下节点i的自导纳
Figure BDA00020712618300000520
如下所示:
Figure BDA0002071261830000061
第f个故障状态下节点i和节点j的互导纳
Figure BDA0002071261830000062
如下所示:
Figure BDA0002071261830000063
故障状态节点负荷削减量约束如下所示:
Figure BDA0002071261830000064
5)改变电力系统中新能源电场的装机容量,并重复步骤3和步骤4,计算不同新能源渗透率下的最小负荷削减量的期望值,从而计算得到新能源渗透率上限和新能源消纳能力。
最小切负荷模型如下所示:
Figure BDA0002071261830000065
式中,s、f和i分别为典型场景、N-1预想故障状态和电力网络节点的编号。Ns、Nf和Nn分别为典型场景、N-1预想故障状态和电力网络节点的数量。ps为第s个场景发生的概率。
Figure BDA0002071261830000066
为第s个场景下的第i个节点在第f个故障状态下的有功负荷削减量。
6)根据不同新能源渗透率下的最小负荷削减量的期望值,建立新能源渗透率和最小负荷削减量期望值的关系曲线,从而判断新能源消纳时系统的安全稳定性。电力系统安全稳定的判断标准为:新能源渗透率ξ<ξmax
建立新能源渗透率和最小负荷削减量期望值关系的主要步骤如下:
6.1)改变风电场和光伏电站的装机容量,从而改变电力系统中的新能源渗透率。新能源渗透率ξ如下所示:
Figure BDA0002071261830000067
式中,Nw和NPV分别为并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。
Figure BDA0002071261830000068
Figure BDA0002071261830000069
分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的装机容量。PLD.max为电力系统中负荷功率的最大值。
6.2)利用最小切负荷模型计算不同新能源渗透率下的最小负荷削减量的期望值。
6.3)当电力系统无法满足步骤3或步骤4中的任一约束条件时,则输出电力系统可接受的新能源渗透率上限ξmax。电力系统的新能源消纳能力M如下所示:
Figure BDA0002071261830000071
式中,M为年度新能源消纳电量。Nw和NPV分别为并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。
Figure BDA0002071261830000072
Figure BDA0002071261830000073
分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的出力情况。Ns为典型场景的数量。ps为第s个典型场景发生的概率。
本发明的技术效果是毋庸置疑的。本发明以负荷削减量为系统安全稳定评价指标,提出了基于最小切负荷模型的新能源消纳评估方法,分析新能源渗透率和负荷削减量之间的关系,从而为统筹协调系统的安全稳定运行和新能源消纳提供指导意见。
本发明的评估方法可引导新能源的合理发展;充分考虑了系统中的电源灵活性,并且考虑了系统的正常运行状态和N-1预想故障状态,分析新能源渗透率和最小负荷削减量的关系,可为保障电力系统安全稳定运行的实际工作提供指导。
本发明适用于水电富集地区,通过灵活性电源蓄水式水电站的优化调度,实现区域内新能源的全额消纳,解决弃风、弃光、弃水问题,该方法充分考虑了系统中的电源灵活性,提高新能源消纳水平,准确评估了系统的新能源消纳能力。
本发明综合考虑了系统的正常运行状态和N-1预想故障状态,优化求解了系统在N-1预想故障状态下的最小负荷削减量的期望值,分析了新能源渗透率和最小负荷削减量期望值的关系以及新能源消纳时系统的安全稳定性,为引导新能源的合理发展以及保障系统安全稳定运行奠定基础。
附图说明
图1是新能源消纳评估方法的流程示意图;
图2是2018年典型日3月21日峰时段生成的典型场景;
图3是新能源渗透率和最小负荷削减量期望值的关系曲线。
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实施例1:
参见图1,一种基于最小切负荷模型的新能源消纳评估方法,主要包括以下步骤:
1)获取T年内电力系统数据,并建立电力系统N-1预想故障集。N-1预想故障集是根据电网的实际运行方式建立的。N-1预想故障集即无法满足电力系统N-1准则时所出现的故障。N-1准则是指在正常运行方式下的电力系统中任一元件(如线路、发电机、变压器等)无故障或因故障断开,电力系统应能保持稳定运行和正常供电,其他元件不过载且系统电压和频率均在允许范围内。
2)考虑源荷功率的相关性,利用K均值聚类方法对电力系统数据进行处理,建立电力系统运行典型场景,并分别计算每个典型场景发生的概率。
建立电力系统运行典型场景的主要步骤如下:
2.1)基于T年内电力系统数据,选出T年内电力系统在每个季节的典型日,并设置典型日的功率采样点和采样时刻。
在每个采样时刻,对所有典型日的功率采样点进行功率数据采样,建立原始数据集,并将典型日划分为高峰、平衡和低谷三个时段。
2.2)根据采样点、采样时刻和原始数据集,以采样时刻为序列建立典型日内高峰、平衡和低谷三个时段内的典型场景。典型日内的典型场景总数N=Ttop+Tlev+Tlow
每个典型场景表示如下:
Figure BDA0002071261830000081
式中,
Figure BDA0002071261830000082
Figure BDA0002071261830000083
分别为高峰、平衡和低谷时段内的第t个采样时刻测得的源荷功率。x=1,2,…,m,m为采样点总数。Ttop、Tlev和Tlow分别为高峰、平衡和低谷时段内的采样时刻总数。
2.3)利用K均值法对典型场景进行聚类,主要步骤如下:
2.3.1)根据典型日高峰、平衡和低谷三个时段内的典型场景,确定高峰、平衡和低谷三个时段内的聚类个数Ktop、Klev和Klow。利用密度参数法优化生成Ktop、Klev和Klow个初始聚类中心。
2.3.2)分别计算每个典型场景和所有聚类中心的欧式距离,并将所述典型场景分配到距离最近的类中。利用均值法迭代更新聚类中心直至收敛,从而在高峰时段建立Ktop个类和Ktop个聚类中心,在平衡时段建立Klev个类和Klev个聚类中心,在低谷时段建立Klow个类和Klow个聚类中心。每个类中的典型场景数量即为每个聚类中心表示的典型场景数量。
2.4)对每个聚类中心和所述聚类中心所在类中的所有典型场景中的源荷功率数据进行比较,找出与聚类中心差距最小的典型场景,也即与聚类中心的源荷功率数据差值最小的典型场景,并利用与聚类中心差距最小的典型场景替换聚类中心,从而保证所生成的典型场景来源于系统中实际的源荷功率场景。每个替换后的聚类中心即为典型日的典型场景。
2.5)第s个典型场景发生的概率Ps如下所示:
Figure BDA0002071261830000091
式中,s为典型场景的编号。Ns为典型场景的总数,Ns个典型场景是由所选取的典型日中的原始场景通过K均值聚类技术生成的,系统中的所有原始场景均是通过Ns个典型场景进行表示。ns为典型场景s所表示的典型场景数量。N为典型日中的典型场景总数。Ts为典型场景s所属的典型日在一年中所代表的天数。
3)对电力系统的灵活性电源蓄水式水电站进行优化调度,使电力系统在正常运行状态下满足正常运行约束条件。
所述正常运行约束条件主要包括风电与光伏发电的消纳约束、水电的消纳约束、正常运行功率平衡约束、线路传输容量约束、机组出力约束、水电站库容约束、水电站发电流量约束和水电站水量平衡约束。
风电与光伏发电的消纳约束如下所示:
Figure BDA0002071261830000092
式中,Nw和NPV分别为电力系统并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。
Figure BDA0002071261830000093
Figure BDA0002071261830000094
分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的出力情况。
Figure BDA0002071261830000095
Figure BDA0002071261830000096
分别为第s个典型场景下风电场和光伏电站总的最大出力。为了防止出现弃风弃光问题,必须使各场景下系统所消纳的风电和光伏发电等于风电场和光伏电站的最大出力。
水电的消纳约束如下所示:
Figure BDA0002071261830000097
式中,Nstop、Nslev和Nslow分别为典型日内高峰平衡低谷时段内的典型场景总数。Nh为并网的水电站总数。sh为水电站的编号。
Figure BDA0002071261830000098
Figure BDA0002071261830000099
分别为高峰平衡低谷各个时段内的第s个典型场景下的第sh座水电站的出力情况。
Figure BDA00020712618300000910
Figure BDA00020712618300000911
分别为高峰平衡低谷各个时段内的第s个典型场景所持续的时间。
Figure BDA00020712618300000912
为水电站在任意典型日的最大发电量。蓄水式水电站具有较大容量的水库,当风电场和光伏电站出力较大时,水库可以将该场景下的多余发电用水量进行储存,用于其它场景时的水力发电。
正常运行功率平衡约束如下所示:
Figure BDA0002071261830000101
其中,
Figure BDA0002071261830000102
Figure BDA0002071261830000103
分别为节点i在第s个典型场景下的火电、光伏、风电和水电的有功注入功率。
Figure BDA0002071261830000104
Figure BDA0002071261830000105
分别为节点i在第s个典型场景下的有功负荷和注入功率。
Figure BDA0002071261830000106
Figure BDA0002071261830000107
分别为节点i与节点j在第s个典型场景下的电压相角。
Figure BDA0002071261830000108
Figure BDA0002071261830000109
分别是以
Figure BDA00020712618300001010
为支路导纳建立起来的节点导纳矩阵的自导纳和互导纳。
正常运行时第s个典型场景的节点i的自导纳
Figure BDA00020712618300001011
如下所示:
Figure BDA00020712618300001012
正常运行时第s个典型场景的节点i和节点j互导纳
Figure BDA00020712618300001013
如下所示:
Figure BDA00020712618300001014
线路传输容量约束如下所示:
Figure BDA00020712618300001015
式中,Pij.max为支路ij之间的传输极限。
Figure BDA00020712618300001016
为支路ij在第s个典型场景下所传输的有功功率。
支路ij在第s个典型场景下所传输的有功功率
Figure BDA00020712618300001017
如下所示:
Figure BDA00020712618300001018
机组出力约束分别如公式10至公式13所示:
Figure BDA00020712618300001019
Figure BDA00020712618300001020
Figure BDA00020712618300001021
Figure BDA0002071261830000111
式中,
Figure BDA0002071261830000112
Figure BDA0002071261830000113
分别为节点i所连火电机组在第s个典型场景下的有功出力上下限。
Figure BDA0002071261830000114
Figure BDA0002071261830000115
分别为节点i所连光伏电站在第s个典型场景下的有功出力上下限。
Figure BDA0002071261830000116
Figure BDA0002071261830000117
分别为节点i所连风电场在第s个典型场景下的有功出力上下限。
Figure BDA0002071261830000118
Figure BDA0002071261830000119
分别为节点i所连接的水电站在第s个典型场景时的有功出力上下限。
水电站发电流量约束如下所示:
Figure BDA00020712618300001110
式中,
Figure BDA00020712618300001111
Figure BDA00020712618300001112
分别为节点i所连接的水电站在第s个典型场景时所分配的发电流量上下限。水库承担着一定量的防洪、农业灌溉和生活用水等任务,水电站需要根据水库调度部门调度计划进行安排发电用水量,可用水电资源应为调度部门分配的发电流量上下限。
发电流量
Figure BDA00020712618300001113
和水电出力
Figure BDA00020712618300001114
的转化关系如下所示:
Figure BDA00020712618300001115
式中,ηi为节点i所连接的水电站的发电效率。
Figure BDA00020712618300001116
为节点i所连接的水电站在第s个典型场景下的发电流量。
Figure BDA00020712618300001117
为节点i所连接的水电站在第s个典型场景下的发电水头平均高度。
水电站水量平衡约束如下所示:
Figure BDA00020712618300001118
其中,
Figure BDA00020712618300001119
Figure BDA00020712618300001120
分别为节点i所连接的水电站在第s个典型场景初时和末时的水库蓄水量。
Figure BDA00020712618300001121
Figure BDA00020712618300001122
分别为节点i所连接的水电站在第s个典型场景时的入库流量、发电流量和其它流量。所述其它流量主要包括灌溉流量和生活用水流量。
水电站库容约束如下所示:
Figure BDA00020712618300001123
式中,
Figure BDA00020712618300001124
Figure BDA00020712618300001125
分别为节点i所连接的水电站在第s个典型场景时的蓄水量上下限。
Figure BDA0002071261830000121
为节点i所连接的水电站在第s个典型场景时的蓄水量。
4)基于N-1预想故障集,建立最小切负荷模型。将电力系统数据输入最小切负荷模型中,得到最小负荷削减量的期望值,使电力系统在N-1预想故障状态下满足故障状态约束条件。
所述故障状态约束条件主要包括风电和光伏发电的消纳约束、水电的消纳约束、故障状态功率平衡约束、线路传输容量约束、机组出力约束、故障状态节点负荷削减量约束、水电站库容约束、水电站发电流量约束、水电站水量平衡约束。
故障状态功率平衡约束如下所示:
Figure BDA0002071261830000122
式中,
Figure BDA0002071261830000123
Figure BDA0002071261830000124
分别为节点i在第f个故障状态下的火电、光伏、风电和水电机组的有功注入功率。
Figure BDA0002071261830000125
为节点i在第f个故障状态下的有功负荷。
Figure BDA0002071261830000126
为节点i在第f个故障状态下的有功负荷削减量。
Figure BDA0002071261830000127
Figure BDA0002071261830000128
分别为节点i与节点j在第f个故障状态下的电压相角。
Figure BDA0002071261830000129
Figure BDA00020712618300001210
分别是在第f个故障状态下以
Figure BDA00020712618300001211
为支路导纳建立起来的节点导纳矩阵的自导纳及互导纳。
第f个故障状态下节点i的自导纳
Figure BDA00020712618300001212
如下所示:
Figure BDA00020712618300001213
式中,
Figure BDA00020712618300001214
为第f个故障状态下支路导纳。
第f个故障状态下节点i和节点j的互导纳
Figure BDA00020712618300001215
如下所示:
Figure BDA00020712618300001216
故障状态节点负荷削减量约束如下所示:
Figure BDA00020712618300001217
5)改变电力系统中新能源电场的装机容量,并重复步骤3和步骤4,计算不同新能源渗透率下的最小负荷削减量的期望值,从而计算得到新能源渗透率上限和新能源消纳能力。
最小切负荷模型如下所示:
Figure BDA0002071261830000131
式中,s、f和i分别为典型场景、N-1预想故障状态和电力网络节点的编号。Ns、Nf和Nn分别为典型场景、N-1预想故障状态和电力网络节点的数量。ps为第s个场景发生的概率。
Figure BDA0002071261830000132
为第s个场景下的第i个节点在第f个故障状态下的有功负荷削减量。
6)根据不同新能源渗透率下的最小负荷削减量的期望值,建立新能源渗透率和最小负荷削减量期望值的关系曲线,从而判断新能源消纳时系统的安全稳定性。电力系统安全稳定的判断标准为:新能源渗透率ξ<ξmax
建立新能源渗透率和最小负荷削减量期望值关系曲线的主要步骤如下:
6.1)结合典型场景中的负荷功率,改变风电场和光伏电站的装机容量,从而改变电力系统中的新能源渗透率。新能源渗透率ξ如下所示:
Figure BDA0002071261830000133
式中,Nw和NPV分别为并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。
Figure BDA0002071261830000134
Figure BDA0002071261830000135
分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的装机容量。PLD.max为电力系统中负荷功率的最大值。
6.2)利用最小切负荷模型计算不同新能源渗透率下的最小负荷削减量的期望值。
6.3)当电力系统无法满足步骤3或步骤4中的任一约束条件时,则输出电力系统可接受的新能源渗透率上限ξmax。电力系统的新能源消纳能力M如下所示:
Figure BDA0002071261830000136
式中,M为年度新能源消纳电量。Nw和NPV分别为并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。
Figure BDA0002071261830000137
Figure BDA0002071261830000138
分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的出力情况。Ns为典型场景的数量。ps为第s个典型场景发生的概率。
实施例2:
一种基于最小切负荷模型的新能源消纳评估方法,主要包括以下步骤:
1)获取T年内电力系统数据,并建立电力系统N-1预想故障集。
2)利用K均值聚类方法对电力系统数据进行处理,建立电力系统运行典型场景,并分别计算每个典型场景发生的概率。
3)对电力系统的水电站进行优化调度,使电力系统在正常运行状态下满足正常运行约束条件。
4)基于N-1预想故障集,建立最小切负荷模型;将电力系统数据输入最小切负荷模型中,得到最小负荷削减量的期望值,使电力系统在N-1预想故障状态下满足故障状态约束条件。
5)改变电力系统中新能源电场的装机容量,并重复步骤3和步骤4,计算不同新能源渗透率下的最小负荷削减量的期望值,从而计算得到新能源渗透率上限和新能源消纳能力。
6)根据不同新能源渗透率下的最小负荷削减量的期望值,建立新能源渗透率和最小负荷削减量期望值的关系曲线,从而判断新能源消纳时系统的安全稳定性。
实施例3:
一种基于最小切负荷模型的新能源消纳评估方法,主要步骤参见实施例2,其中,建立电力系统运行典型场景的主要步骤如下:
1)基于T年内电力系统数据,选出T年内电力系统在每个季节的典型日,并设置典型日的功率采样点和采样时刻。
在每个采样时刻,对所有典型日的功率采样点进行功率数据采样,建立原始数据集,并将典型日划分为高峰、平衡和低谷三个时段。
2)根据采样点、采样时刻和原始数据集,建立典型日内高峰、平衡和低谷三个时段内的典型场景。典型日内的典型场景总数N=Ttop+Tlev+Tlow
每个典型场景表示如下:
Figure BDA0002071261830000141
式中,
Figure BDA0002071261830000142
Figure BDA0002071261830000143
分别为高峰、平衡和低谷时段内的第t个采样时刻测得的源荷功率。x=1,2,…,m,m为采样点总数。Ttop、Tlev和Tlow分别为高峰、平衡和低谷时段内的采样时刻总数。
3)利用K均值法对典型场景进行聚类,主要步骤如下:
3.1)根据典型日高峰、平衡和低谷三个时段内的典型场景,确定高峰、平衡和低谷三个时段内的聚类个数Ktop、Klev和Klow。利用密度参数法优化生成Ktop、Klev和Klow个初始聚类中心。
3.2)分别计算每个典型场景和所有聚类中心的欧式距离,并将所述典型场景分配到距离最近的类中。利用均值法迭代更新聚类中心直至收敛,从而在高峰时段建立Ktop个类和Ktop个聚类中心,在平衡时段建立Klev个类和Klev个聚类中心,在低谷时段建立Klow个类和Klow个聚类中心。每个类中的典型场景数量即为每个聚类中心表示的典型场景数量。
4)对每个聚类中心和所述聚类中心所在类中的所有典型场景中的源荷功率数据进行比较,找出与聚类中心差距最小的典型场景,并利用与聚类中心差距最小的典型场景替换聚类中心。每个替换后的聚类中心即为典型日的典型场景。
5)第s个典型场景发生的概率Ps如下所示:
Figure BDA0002071261830000151
式中,s为典型场景的编号。Ns为典型场景的总数。ns为典型场景s所表示的典型场景数量。N为典型日中的典型场景总数。Ts为典型场景s所属的典型日在一年中所代表的天数。
实施例4:
一种基于最小切负荷模型的新能源消纳评估方法,主要步骤参见实施例2,其中,正常运行约束条件主要包括风电与光伏发电的消纳约束、水电的消纳约束、正常运行功率平衡约束、线路传输容量约束、机组出力约束、水电站库容约束、水电站发电流量约束和水电站水量平衡约束。
风电与光伏发电的消纳约束如下所示:
Figure BDA0002071261830000152
式中,Nw和NPV分别为电力系统并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。
Figure BDA0002071261830000153
Figure BDA0002071261830000154
分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的出力情况。
Figure BDA0002071261830000155
Figure BDA0002071261830000156
分别为第s个典型场景下风电场和光伏电站总的最大出力。
水电的消纳约束如下所示:
Figure BDA0002071261830000157
式中,Nstop、Nslev和Nslow分别为典型日内高峰平衡低谷时段内的典型场景总数。Nh为并网的水电站总数。sh为水电站的编号。
Figure BDA0002071261830000158
Figure BDA0002071261830000159
分别为高峰平衡低谷各个时段内的第s个典型场景下的第sh座水电站的出力情况。
Figure BDA00020712618300001510
Figure BDA00020712618300001511
分别为高峰平衡低谷各个时段内的第s个典型场景所持续的时间。
Figure BDA0002071261830000161
为水电站在任意典型日的最大发电量。
正常运行功率平衡约束如下所示:
Figure BDA0002071261830000162
其中,
Figure BDA0002071261830000163
Figure BDA0002071261830000164
分别为节点i在第s个典型场景下的火电、光伏、风电和水电的有功注入功率。
Figure BDA0002071261830000165
Figure BDA0002071261830000166
分别为节点i在第s个典型场景下的有功负荷和注入功率。
Figure BDA0002071261830000167
Figure BDA0002071261830000168
分别为节点i与节点j在第s个典型场景下的电压相角。
Figure BDA0002071261830000169
Figure BDA00020712618300001610
分别是以
Figure BDA00020712618300001611
为支路导纳建立起来的节点导纳矩阵的自导纳和互导纳。
正常运行时第s个典型场景的节点i的自导纳
Figure BDA00020712618300001612
如下所示:
Figure BDA00020712618300001613
正常运行时第s个典型场景的节点i和节点j互导纳
Figure BDA00020712618300001614
如下所示:
Figure BDA00020712618300001615
线路传输容量约束如下所示:
Figure BDA00020712618300001616
式中,Pij.max为支路ij之间的传输极限。
Figure BDA00020712618300001617
为支路ij在第s个典型场景下所传输的有功功率。
支路ij在第s个典型场景下所传输的有功功率
Figure BDA00020712618300001618
如下所示:
Figure BDA00020712618300001619
机组出力约束分别如公式8至公式11所示:
Figure BDA00020712618300001620
Figure BDA00020712618300001621
Figure BDA00020712618300001622
Figure BDA0002071261830000171
式中,
Figure BDA0002071261830000172
Figure BDA0002071261830000173
分别为节点i所连火电机组在第s个典型场景下的有功出力上下限。
Figure BDA0002071261830000174
Figure BDA0002071261830000175
分别为节点i所连光伏电站在第s个典型场景下的有功出力上下限。
Figure BDA0002071261830000176
Figure BDA0002071261830000177
分别为节点i所连风电场在第s个典型场景下的有功出力上下限。
Figure BDA0002071261830000178
Figure BDA0002071261830000179
分别为节点i所连接的水电站在第s个典型场景时的有功出力上下限。
水电站发电流量约束如下所示:
Figure BDA00020712618300001710
式中,
Figure BDA00020712618300001711
Figure BDA00020712618300001712
分别为节点i所连接的水电站在第s个典型场景时所分配的发电流量上下限。
发电流量
Figure BDA00020712618300001713
和水电出力
Figure BDA00020712618300001714
的转化关系如下所示:
Figure BDA00020712618300001715
式中,ηi为节点i所连接的水电站的发电效率。
Figure BDA00020712618300001716
为节点i所连接的水电站在第s个典型场景下的发电流量。
Figure BDA00020712618300001717
为节点i所连接的水电站在第s个典型场景下的发电水头平均高度。
水电站水量平衡约束如下所示:
Figure BDA00020712618300001718
其中,
Figure BDA00020712618300001719
Figure BDA00020712618300001720
分别为节点i所连接的水电站在第s个典型场景初时和末时的水库蓄水量。
Figure BDA00020712618300001721
Figure BDA00020712618300001722
分别为节点i所连接的水电站在第s个典型场景时的入库流量、发电流量和其它流量。所述其它流量主要包括灌溉流量和生活用水流量。
水电站库容约束如下所示:
Figure BDA00020712618300001723
式中,
Figure BDA00020712618300001724
Figure BDA00020712618300001725
分别为节点i所连接的水电站在第s个典型场景时的蓄水量上下限。
Figure BDA00020712618300001726
为节点i所连接的水电站在第s个典型场景时的蓄水量。
实施例5:
一种基于最小切负荷模型的新能源消纳评估方法,主要步骤参见实施例2,其中,故障状态约束条件主要包括风电和光伏发电的消纳约束、水电的消纳约束、故障状态功率平衡约束、线路传输容量约束、机组出力约束、故障状态节点负荷削减量约束、水电站库容约束、水电站发电流量约束、水电站水量平衡约束。
故障状态功率平衡约束如下所示:
Figure BDA0002071261830000181
式中,
Figure BDA0002071261830000182
Figure BDA0002071261830000183
分别为节点i在第f个故障状态下的火电、光伏、风电和水电机组的有功注入功率。
Figure BDA0002071261830000184
为节点i在第f个故障状态下的有功负荷。
Figure BDA0002071261830000185
为节点i在第f个故障状态下的有功负荷削减量。
Figure BDA0002071261830000186
Figure BDA0002071261830000187
分别为节点i与节点j在第f个故障状态下的电压相角。
Figure BDA0002071261830000188
Figure BDA0002071261830000189
分别是在第f个故障状态下以
Figure BDA00020712618300001810
为支路导纳建立起来的节点导纳矩阵的自导纳及互导纳。
第f个故障状态下节点i的自导纳
Figure BDA00020712618300001811
如下所示:
Figure BDA00020712618300001812
第f个故障状态下节点i和节点j的互导纳
Figure BDA00020712618300001813
如下所示:
Figure BDA00020712618300001814
故障状态节点负荷削减量约束如下所示:
Figure BDA00020712618300001815
风电与光伏发电的消纳约束如下所示:
Figure BDA00020712618300001816
式中,Nw和NPV分别为电力系统并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。
Figure BDA00020712618300001817
Figure BDA00020712618300001818
分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的出力情况。
Figure BDA00020712618300001819
Figure BDA00020712618300001820
分别为第s个典型场景下风电场和光伏电站总的最大出力。
水电的消纳约束如下所示:
Figure BDA0002071261830000191
式中,Nstop、Nslev和Nslow分别为典型日内高峰平衡低谷时段内的典型场景总数。Nh为并网的水电站总数。sh为水电站的编号。
Figure BDA0002071261830000192
Figure BDA0002071261830000193
分别为高峰平衡低谷各个时段内的第s个典型场景下的第sh座水电站的出力情况。
Figure BDA0002071261830000194
Figure BDA0002071261830000195
分别为高峰平衡低谷各个时段内的第s个典型场景所持续的时间。
Figure BDA0002071261830000196
为水电站在任意典型日的最大发电量。
线路传输容量约束如下所示:
Figure BDA0002071261830000197
式中,Pij.max为支路ij之间的传输极限。
Figure BDA0002071261830000198
为支路ij在第s个典型场景下所传输的有功功率。
支路ij在第s个典型场景下所传输的有功功率
Figure BDA0002071261830000199
如下所示:
Figure BDA00020712618300001910
机组出力约束分别如公式9至公式12所示:
Figure BDA00020712618300001911
Figure BDA00020712618300001912
Figure BDA00020712618300001913
Figure BDA00020712618300001914
式中,
Figure BDA00020712618300001915
Figure BDA00020712618300001916
分别为节点i所连火电机组在第s个典型场景下的有功出力上下限。
Figure BDA00020712618300001917
Figure BDA00020712618300001918
分别为节点i所连光伏电站在第s个典型场景下的有功出力上下限。
Figure BDA00020712618300001919
Figure BDA00020712618300001920
分别为节点i所连风电场在第s个典型场景下的有功出力上下限。
Figure BDA00020712618300001921
Figure BDA00020712618300001922
分别为节点i所连接的水电站在第s个典型场景时的有功出力上下限。
水电站发电流量约束如下所示:
Figure BDA00020712618300001923
式中,
Figure BDA0002071261830000201
Figure BDA0002071261830000202
分别为节点i所连接的水电站在第s个典型场景时所分配的发电流量上下限。
发电流量
Figure BDA0002071261830000203
和水电出力
Figure BDA0002071261830000204
的转化关系如下所示:
Figure BDA0002071261830000205
式中,ηi为节点i所连接的水电站的发电效率。
Figure BDA0002071261830000206
为节点i所连接的水电站在第s个典型场景下的发电流量。
Figure BDA0002071261830000207
为节点i所连接的水电站在第s个典型场景下的发电水头平均高度。
水电站水量平衡约束如下所示:
Figure BDA0002071261830000208
其中,
Figure BDA0002071261830000209
Figure BDA00020712618300002010
分别为节点i所连接的水电站在第s个典型场景初时和末时的水库蓄水量。
Figure BDA00020712618300002011
Figure BDA00020712618300002012
分别为节点i所连接的水电站在第s个典型场景时的入库流量、发电流量和其它流量。所述其它流量主要包括灌溉流量和生活用水流量。
水电站库容约束如下所示:
Figure BDA00020712618300002013
式中,
Figure BDA00020712618300002014
Figure BDA00020712618300002015
分别为节点i所连接的水电站在第s个典型场景时的蓄水量上下限。
Figure BDA00020712618300002016
为节点i所连接的水电站在第s个典型场景时的蓄水量。
实施例6:
一种基于最小切负荷模型的新能源消纳评估方法,主要步骤参见实施例2,其中,最小切负荷模型如下所示:
Figure BDA00020712618300002017
式中,s、f和i分别为典型场景、N-1预想故障状态和电力网络节点的编号。Ns、Nf和Nn分别为典型场景、N-1预想故障状态和电力网络节点的数量。ps为第s个场景发生的概率。
Figure BDA00020712618300002018
为第s个场景下的第i个节点在第f个故障状态下的有功负荷削减量。
实施例7:
一种基于最小切负荷模型的新能源消纳评估方法,主要步骤参见实施例2,其中,建立新能源渗透率和最小负荷削减量期望值关系的主要步骤如下:
6.1)改变风电场和光伏电站的装机容量,从而改变电力系统中的新能源渗透率。新能源渗透率ξ如下所示:
Figure BDA0002071261830000211
式中,Nw和NPV分别为并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。
Figure BDA0002071261830000212
Figure BDA0002071261830000213
分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的装机容量。PLD.max为电力系统中负荷功率的最大值。
6.2)利用最小切负荷模型计算不同新能源渗透率下的最小负荷削减量的期望值。
6.3)当电力系统无法满足步骤3或步骤4中的任一约束条件时,则输出电力系统可接受的新能源渗透率上限ξmax。电力系统的新能源消纳能力M如下所示:
式中,M为年度新能源消纳电量。Nw和NPV分别为并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。
Figure BDA0002071261830000215
Figure BDA0002071261830000216
分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的出力情况。Ns为典型场景的数量。ps为第s个典型场景发生的概率。
实施例8:
一种验证基于最小切负荷模型的新能源消纳评估方法的实验,主要步骤如下:
1)获取T年内电力系统数据,并建立电力系统N-1预想故障集。以我国南方某沿海地区110kV电网2018年的电源结构、网架结构、实际运行方式等信息作为实施例的原始数据,2018年该地有3座新能源电场,总容量达到220.24MW;1座火电厂,总容量达到1200MW,4座蓄水式水电站,总容量达到69.6MW,2018年该区域负荷的最大功率为1080.49MW,网络N-1预想故障集中的故障数量为97个。
2)利用K均值聚类方法对电力系统数据进行处理,建立电力系统运行典型场景,并分别计算每个典型场景发生的概率。
在冬、夏和春/秋季节内各选取一个典型日来反映当地电网2018年的运行情况。其中,冬季选择了12月21日,夏季选择了8月5日,春/秋季选择了3月21日。
选择该区域电网中的126个节点作为功率采样点,每隔5min对各采样点进行一次功率采样获得源荷功率的原始数据,结合该区域电力系统的实际运行情况,典型日内的高峰、平衡和低谷时段划分情况如表1所示。
表1典型日内的高峰、平衡和低谷时段划分情况
时段 持续时间
高峰时段 10:00-12:00、15:00-20:00
平衡时段 7:00-10:00、12:00-15:00、20:00-0:00
低谷时段 0:00-7:00
根据采样获得的原始数据,采用K均值聚类技术分别在3个典型日的高峰、平衡和低谷时段内均生成3个典型场景,并计算各个典型场景发生的概率。
图2为2018年春/秋季典型日3月21日高峰时段的典型场景的生成结果;表2为2018年春/秋季典型日3月21日高峰时段的典型场景发生的概率;依据2018年3月21日平衡时段、低谷时段和冬季、夏季典型日各时段内的原始场景,同样可以生成对应的典型场景和获得各典型场景发生的概率。
表2 2018年3月21日高峰时段的典型场景发生的概率
Figure BDA0002071261830000221
3)对电力系统的水电站进行优化调度,使电力系统在正常运行状态下满足正常运行约束条件。
4)基于N-1预想故障集,建立最小切负荷模型。将电力系统数据输入最小切负荷模型中,利用CPLEX求解器求解得到最小负荷削减量的期望值,使电力系统在N-1预想故障状态下满足故障状态约束条件。
基于测试网络和典型场景数据,在正常运行状态下,实现系统内的新能源全额消纳,并且避免出现线路过载和切负荷等问题;在故障状态下,根据N-1安全校核的预想故障集,计算出各个典型场景下的最小负荷削减量,计算结果如表3所示。
表3各典型场景下的最小负荷削减量(MW)
Figure BDA0002071261830000231
进一步结合预想故障数量和各个典型场景概率,计算2018年该区域最小负荷削减量的期望值为0.2881MW;根据新能源电场的装机容量以及2018年该区域最大的负荷功率,计算2018年该区域的新能源渗透率为20.38%,新能源消纳电量为415601.62MW·h。
5)改变电力系统中新能源电场的装机容量,并重复步骤3和步骤4,计算不同新能源渗透率下的最小负荷削减量的期望值,从而计算得到新能源渗透率上限和新能源消纳能力。
调整新能源电场的装机容量来改变系统中的新能源渗透率,求解不同新能源渗透率下的最小负荷削减量的期望值。其中,在改变新能源渗透率时仍需使系统满足正常运行状态下和预想故障状态下的约束条件,当不满足时,即可获得系统的新能源渗透率上限,经计算系统的新能源渗透率上限为84.58%。
6)根据不同新能源渗透率下的最小负荷削减量的期望值,建立新能源渗透率和最小负荷削减量期望值的关系曲线,分析两者之间的关系,进而分析新能源消纳时系统的安全稳定性。
根据不同新能源渗透率下的最小负荷削减量的期望值,建立新能源渗透率和最小负荷削减量期望值的关系曲线如图3所示。
从图3可以看出,当新能源渗透率小于68.74%时,最小负荷削减量的期望值为固定值,即使新能源渗透率为0,没有新能源接入时,最小负荷削减量的期望值仍是0.2881MW,说明新能源渗透率在0-68.74%区间内,系统最小负荷削减量的期望值与新能源是否接入、新能源渗透率是否提高均没有关系;当新能源渗透率高于68.74%时,最小负荷削减量的期望值同新能源渗透率成正相关,系统最小负荷削减量的期望值将随着新能源消纳电量的增加而增大,进而给系统的安全稳定运行带来更多的隐患;受限于新能源电场输电线路传输容量的限制,新能源渗透率的上限为84.58%,对应最小负荷削减量的期望值最大为0.7592MW;当新能源渗透率超过上限时,系统在正常运行状态下会出现线路过载的现象,无法安全稳定运行。
根据新能源渗透率的上限,评估系统的新能源消纳能力,经计算,系统年度新能源消纳电量最大为3592193.13MW·h。

Claims (7)

1.一种基于最小切负荷模型的新能源消纳评估方法,其特征在于,主要包括以下步骤:
1)获取T年内电力系统数据,并建立电力系统N-1预想故障集;
2)利用K均值聚类方法对电力系统数据进行处理,建立电力系统运行典型场景,并分别计算每个典型场景发生的概率;
3)对电力系统的水电站进行优化调度,使电力系统在正常运行状态下满足正常运行约束条件;
4)基于N-1预想故障集,建立最小切负荷模型;将电力系统数据输入最小切负荷模型中,得到最小负荷削减量的期望值,使电力系统在N-1预想故障状态下满足故障状态约束条件;
5)改变电力系统中新能源电场的装机容量,并重复步骤3和步骤4,计算不同新能源渗透率下的最小负荷削减量的期望值,从而计算得到新能源渗透率上限和新能源消纳能力;
6)根据不同新能源渗透率下的最小负荷削减量的期望值,建立新能源渗透率和最小负荷削减量期望值的关系曲线,从而判断新能源消纳时系统的安全稳定性。
2.根据权利要求1所述的一种基于最小切负荷模型的新能源消纳评估方法,其特征在于,建立电力系统运行典型场景的主要步骤如下:
1)基于T年内电力系统数据,选出T年内电力系统在每个季节的典型日,并设置典型日的功率采样点和采样时刻;
在每个采样时刻,对所有典型日的功率采样点进行功率数据采样,建立原始数据集,并将典型日划分为高峰、平衡和低谷三个时段;
2)根据采样点、采样时刻和原始数据集,建立典型日内高峰、平衡和低谷三个时段内的典型场景;典型日内的典型场景总数N=Ttop+Tlev+Tlow
每个典型场景表示如下:
Figure FDA0002554113500000011
式中,
Figure FDA0002554113500000012
Figure FDA0002554113500000013
分别为高峰、平衡和低谷时段内的第t个采样时刻测得的源荷功率;x=1,2,…,m,m为采样点总数;Ttop、Tlev和Tlow分别为高峰、平衡和低谷时段内的采样时刻总数;
3)利用K均值法对典型场景进行聚类,主要步骤如下:
3.1)根据典型日高峰、平衡和低谷三个时段内的典型场景,确定高峰、平衡和低谷三个时段内的聚类个数Ktop、Klev和Klow;利用密度参数法优化生成Ktop、Klev和Klow个初始聚类中心;
3.2)分别计算每个典型场景和所有聚类中心的欧式距离,并将所述典型场景分配到距离最近的类中;利用均值法迭代更新聚类中心直至收敛,从而在高峰时段建立Ktop个类和Ktop个聚类中心,在平衡时段建立Klev个类和Klev个聚类中心,在低谷时段建立Klow个类和Klow个聚类中心;每个类中的典型场景数量即为每个聚类中心表示的典型场景数量;
4)对每个聚类中心和所述聚类中心所在类中的所有典型场景中的源荷功率数据进行比较,找出与聚类中心差距最小的典型场景,并利用与聚类中心差距最小的典型场景替换聚类中心;每个替换后的聚类中心即为典型日的典型场景;
5)第s个典型场景发生的概率Ps如下所示:
Figure FDA0002554113500000021
式中,s为典型场景的编号;Ns为典型场景的总数;ns为典型场景s所表示的典型场景数量;N为典型日中的典型场景总数;Ts为典型场景s所属的典型日在一年中所代表的天数。
3.根据权利要求1所述的一种基于最小切负荷模型的新能源消纳评估方法,其特征在于:所述正常运行约束条件主要包括风电与光伏发电的消纳约束、水电的消纳约束、正常运行功率平衡约束、线路传输容量约束、机组出力约束、水电站库容约束、水电站发电流量约束和水电站水量平衡约束;
风电与光伏发电的消纳约束如下所示:
Figure FDA0002554113500000022
式中,Nw和NPV分别为电力系统并网的风电场和光伏电站总数;sw和sPV分别为风电场和光伏电站的编号;
Figure FDA0002554113500000023
Figure FDA0002554113500000024
分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的出力情况;
Figure FDA0002554113500000025
Figure FDA0002554113500000026
分别为第s个典型场景下风电场和光伏电站总的最大出力;
水电的消纳约束如下所示:
Figure FDA0002554113500000027
式中,Nstop、Nslev和Nslow分别为典型日内高峰平衡低谷时段内的典型场景总数;Nh为并网的水电站总数;sh为水电站的编号;
Figure FDA0002554113500000028
Figure FDA0002554113500000029
分别为高峰平衡低谷各个时段内的第s个典型场景下的第sh座水电站的出力情况;
Figure FDA00025541135000000210
Figure FDA00025541135000000211
分别为高峰平衡低谷各个时段内的第s个典型场景所持续的时间;
Figure FDA00025541135000000212
为水电站在任意典型日的最大发电量;
正常运行功率平衡约束如下所示:
Figure FDA0002554113500000031
其中,
Figure FDA0002554113500000032
Wi s
Figure FDA0002554113500000033
分别为节点i在第s个典型场景下的火电、光伏、风电和水电的有功注入功率;
Figure FDA0002554113500000034
和Pi s分别为节点i在第s个典型场景下的有功负荷和注入功率;
Figure FDA0002554113500000035
Figure FDA0002554113500000036
分别为节点i与节点j在第s个典型场景下的电压相角;
Figure FDA0002554113500000037
Figure FDA0002554113500000038
分别是以
Figure FDA0002554113500000039
为支路导纳建立起来的节点导纳矩阵的自导纳和互导纳;
正常运行时第s个典型场景的节点i的自导纳
Figure FDA00025541135000000310
如下所示:
Figure FDA00025541135000000311
正常运行时第s个典型场景的节点i和节点j互导纳
Figure FDA00025541135000000312
如下所示:
Figure FDA00025541135000000313
线路传输容量约束如下所示:
Figure FDA00025541135000000314
式中,Pij.max为支路ij之间的传输极限;
Figure FDA00025541135000000315
为支路ij在第s个典型场景下所传输的有功功率;
支路ij在第s个典型场景下所传输的有功功率
Figure FDA00025541135000000316
如下所示:
Figure FDA00025541135000000317
机组出力约束分别如公式10至公式13所示:
Figure FDA00025541135000000318
Figure FDA00025541135000000319
Figure FDA00025541135000000320
Figure FDA00025541135000000321
式中,
Figure FDA0002554113500000041
Figure FDA0002554113500000042
分别为节点i所连火电机组在第s个典型场景下的有功出力上下限;
Figure FDA0002554113500000043
Figure FDA0002554113500000044
分别为节点i所连光伏电站在第s个典型场景下的有功出力上下限;
Figure FDA0002554113500000045
Figure FDA0002554113500000046
分别为节点i所连风电场在第s个典型场景下的有功出力上下限;
Figure FDA0002554113500000047
Figure FDA0002554113500000048
分别为节点i所连接的水电站在第s个典型场景时的有功出力上下限;
水电站发电流量约束如下所示:
Figure FDA0002554113500000049
式中,
Figure FDA00025541135000000410
Figure FDA00025541135000000411
分别为节点i所连接的水电站在第s个典型场景时所分配的发电流量上下限;
发电流量
Figure FDA00025541135000000412
和水电出力
Figure FDA00025541135000000413
的转化关系如下所示:
Figure FDA00025541135000000414
式中,ηi为节点i所连接的水电站的发电效率;
Figure FDA00025541135000000415
为节点i所连接的水电站在第s个典型场景下的发电流量;
Figure FDA00025541135000000416
为节点i所连接的水电站在第s个典型场景下的发电水头平均高度;
水电站水量平衡约束如下所示:
Figure FDA00025541135000000417
其中,Vi s.t和Vi s.t+1分别为节点i所连接的水电站在第s个典型场景初时和末时的水库蓄水量;
Figure FDA00025541135000000418
Figure FDA00025541135000000419
分别为节点i所连接的水电站在第s个典型场景时的入库流量、发电流量和其它流量;所述其它流量主要包括灌溉流量和生活用水流量;
水电站库容约束如下所示:
Figure FDA00025541135000000420
式中,
Figure FDA00025541135000000421
Figure FDA00025541135000000422
分别为节点i所连接的水电站在第s个典型场景时的蓄水量上下限;Vi s为节点i所连接的水电站在第s个典型场景时的蓄水量。
4.根据权利要求1或2所述的一种基于最小切负荷模型的新能源消纳评估方法,其特征在于,最小切负荷模型如下所示:
Figure FDA0002554113500000051
式中,s、f和i分别为典型场景、N-1预想故障状态和电力网络节点的编号;Ns、Nf和Nn分别为典型场景、N-1预想故障状态和电力网络节点的数量;ps为第s个场景发生的概率;
Figure FDA0002554113500000052
为第s个场景下的第i个节点在第f个故障状态下的有功负荷削减量。
5.根据权利要求1或3所述的一种基于最小切负荷模型的新能源消纳评估方法,其特征在于:所述故障状态约束条件主要包括风电和光伏发电的消纳约束、水电的消纳约束、故障状态功率平衡约束、线路传输容量约束、机组出力约束、故障状态节点负荷削减量约束、水电站库容约束、水电站发电流量约束、水电站水量平衡约束;
故障状态功率平衡约束如下所示:
Figure FDA0002554113500000053
式中,
Figure FDA0002554113500000054
Figure FDA0002554113500000055
分别为节点i在第f个故障状态下的火电、光伏、风电和水电机组的有功注入功率;
Figure FDA0002554113500000056
为节点i在第f个故障状态下的有功负荷;
Figure FDA0002554113500000057
为节点i在第f个故障状态下的有功负荷削减量;
Figure FDA0002554113500000058
Figure FDA0002554113500000059
分别为节点i与节点j在第f个故障状态下的电压相角;
Figure FDA00025541135000000510
Figure FDA00025541135000000511
分别是在第f个故障状态下以
Figure FDA00025541135000000512
为支路导纳建立起来的节点导纳矩阵的自导纳及互导纳;
第f个故障状态下节点i的自导纳
Figure FDA00025541135000000513
如下所示:
Figure FDA00025541135000000514
第f个故障状态下节点i和节点j的互导纳
Figure FDA00025541135000000515
如下所示:
Figure FDA00025541135000000516
故障状态节点负荷削减量约束如下所示:
Figure FDA00025541135000000517
6.根据权利要求1所述的一种基于最小切负荷模型的新能源消纳评估方法,其特征在于,建立新能源渗透率和最小负荷削减量期望值关系的主要步骤如下:
1)改变风电场和光伏电站的装机容量,从而改变电力系统中的新能源渗透率;新能源渗透率ξ如下所示:
Figure FDA0002554113500000061
式中,Nw和NPV分别为并网的风电场和光伏电站总数;sw和sPV分别为风电场和光伏电站的编号;
Figure FDA0002554113500000062
Figure FDA0002554113500000063
分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的装机容量;PLD.max为电力系统中负荷功率的最大值;
2)利用最小切负荷模型计算不同新能源渗透率下的最小负荷削减量的期望值;
3)当电力系统无法满足步骤3或步骤4中的任一约束条件时,则输出电力系统可接受的新能源渗透率上限ξmax;电力系统的新能源消纳能力M如下所示:
Figure FDA0002554113500000064
式中,M为年度新能源消纳电量;Nw和NPV分别为并网的风电场和光伏电站总数;sw和sPV分别为风电场和光伏电站的编号;
Figure FDA0002554113500000065
Figure FDA0002554113500000066
分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的出力情况;Ns为典型场景的数量;ps为第s个典型场景发生的概率。
7.根据权利要求1所述的一种基于最小切负荷模型的新能源消纳评估方法,其特征在于,电力系统安全稳定的判断标准为:新能源渗透率ξ<ξmax
CN201910438343.9A 2019-05-24 2019-05-24 一种基于最小切负荷模型的新能源消纳评估方法 Active CN110224393B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910438343.9A CN110224393B (zh) 2019-05-24 2019-05-24 一种基于最小切负荷模型的新能源消纳评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910438343.9A CN110224393B (zh) 2019-05-24 2019-05-24 一种基于最小切负荷模型的新能源消纳评估方法

Publications (2)

Publication Number Publication Date
CN110224393A CN110224393A (zh) 2019-09-10
CN110224393B true CN110224393B (zh) 2020-08-25

Family

ID=67818095

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910438343.9A Active CN110224393B (zh) 2019-05-24 2019-05-24 一种基于最小切负荷模型的新能源消纳评估方法

Country Status (1)

Country Link
CN (1) CN110224393B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110729767A (zh) * 2019-10-09 2020-01-24 广东电网有限责任公司阳江供电局 一种含水电的地区电网风光容量优化配置方法
CN110717694B (zh) * 2019-10-28 2023-02-03 国网河南省电力公司电力科学研究院 基于新能源消纳期望值的储能配置随机决策方法及装置
CN111030166A (zh) * 2019-11-29 2020-04-17 国网甘肃省电力公司建设分公司 一种基于预想故障下的可再生能源出力调控运行方法
CN111709109B (zh) * 2020-04-28 2022-06-07 中国能源建设集团江苏省电力设计院有限公司 一种考虑源荷时序相关性的光伏消纳容量计算方法及装置
CN112149309B (zh) * 2020-09-30 2023-07-21 贵州电网有限责任公司 一种考虑温度-负荷相关性的调温负荷曲线拟合方法
CN112417643B (zh) * 2020-10-13 2023-05-30 国网山东省电力公司电力科学研究院 基于送风机电流的火电机组最大出力实时评估方法及系统
CN113315171B (zh) * 2021-05-18 2022-10-18 南方电网能源发展研究院有限责任公司 新能源消纳能力评估方法、装置、计算机设备和存储介质
CN113098075B (zh) * 2021-05-24 2022-07-12 广东电网有限责任公司 光伏消纳能力的分析方法、装置、电子设备及存储介质
CN113779768B (zh) * 2021-08-16 2024-05-17 西安交通大学 一种需求响应模型构建方法、电子设备及存储介质
CN114925893B (zh) * 2022-05-09 2024-06-04 重庆大学 一种考虑源荷时序模拟与水库月间协调的跨省跨区全年购电策略分层优化方法
CN115641152B (zh) * 2022-12-07 2023-05-12 国网浙江省电力有限公司宁波供电公司 综合能源系统故障识别方法、装置、设备和存储介质
CN115912384B (zh) * 2022-12-21 2024-05-28 国网湖北省电力有限公司经济技术研究院 改善特高压直流受端含高比例新能源系统的电压稳定方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107292502A (zh) * 2017-06-13 2017-10-24 南京软核科技有限公司 一种配电网可靠性评估方法
CN107808216A (zh) * 2017-10-24 2018-03-16 重庆大学 电‑气‑热互联系统弃风弃光和电气热负荷削减综合最小优化模型构建方法
CN108183512A (zh) * 2018-02-23 2018-06-19 南方电网科学研究院有限责任公司 一种接入新能源的电力系统的可靠性评估方法
CN109449923A (zh) * 2018-10-16 2019-03-08 南方电网科学研究院有限责任公司 一种有源配电系统运行灵活性的量化分析方法及相关产品
CN109636671A (zh) * 2018-12-18 2019-04-16 中南大学 一种考虑多层次市场联动型的风电消纳优化策略

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107292502A (zh) * 2017-06-13 2017-10-24 南京软核科技有限公司 一种配电网可靠性评估方法
CN107808216A (zh) * 2017-10-24 2018-03-16 重庆大学 电‑气‑热互联系统弃风弃光和电气热负荷削减综合最小优化模型构建方法
CN108183512A (zh) * 2018-02-23 2018-06-19 南方电网科学研究院有限责任公司 一种接入新能源的电力系统的可靠性评估方法
CN109449923A (zh) * 2018-10-16 2019-03-08 南方电网科学研究院有限责任公司 一种有源配电系统运行灵活性的量化分析方法及相关产品
CN109636671A (zh) * 2018-12-18 2019-04-16 中南大学 一种考虑多层次市场联动型的风电消纳优化策略

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"Intermittent Renewable Energy Source (IRES) Penetration Level into Bangka Power System";Handrea Bernando Tambunan 等;《2018 International Conference on Smart Green Technology in Electrical and Information Systems (ICSGTEIS)》;20190509;全文 *
"基于K- 均值聚类多场景时序特性分析的分布式电源多目标规划";彭春华 等;《电力自动化设备》;20151031;第35卷(第10期);全文 *
"计及预测可信度时段差异性的含分布式可再生能源配电网调度策略";马建伟 等;《电力系统及其自动化学报》;20180531;第30卷(第5期);全文 *

Also Published As

Publication number Publication date
CN110224393A (zh) 2019-09-10

Similar Documents

Publication Publication Date Title
CN110224393B (zh) 一种基于最小切负荷模型的新能源消纳评估方法
CN106874630B (zh) 基于电量消纳的区域电网新能源开发潜力评估方法
CN106786799B (zh) 一种直流联络线功率阶梯化发电计划优化方法
CN110909954B (zh) 一种最大化可再生能源利用的多阶段电源规划方法
CN113489003B (zh) 一种考虑风光水一体化互补运行的源网协调规划方法
CN110994606B (zh) 一种基于复杂适应系统理论的多能源电源容量配置方法
CN112531786A (zh) 区域电网电力外送方法
CN105244922B (zh) 限风情况下双馈型风电场群有功出力优化分配方法
CN110460090B (zh) 一种水风光微电网电源容量配置方法
CN117175543A (zh) 一种可调节负荷的配电网规划策略优化方法及系统
CN112836957B (zh) 一种考虑供电可靠性的区域综合能源系统规划方法
CN108039728B (zh) 提升电网清洁能源消纳能力的实时调度方法
Wu et al. Techno-economic analysis of contingency reserve allocation scheme for combined UHV DC and AC receiving-end power system
CN106485593B (zh) 一种电网侧新能源发电受限原因分析方法
Crosara et al. Generation adequacy in the Nordic and Baltic area: The potential of flexible residential electric heating
CN109447456A (zh) 计及短期需求响应的风电并网电力系统规划方法
Yang et al. Locational marginal pricing-oriented transmission network expansion planning with consideration of line capacity sensitivity
Li et al. Study on grid planning method considering multiple energy access
Zhang Study on the Effects of Different Measures in Promoting Renewable Energy Consumption
Mazumder et al. Integration of solar power in smart grid, perspective: Chittagong (Bangladesh)
Cai et al. Maintenance Schedule of Power Transmission Equipment Considering Large-scale Renewable Energy Access
CN117763960A (zh) 一种偏远地区微电网运行评估方法及系统
Farhat et al. Hive visualizer: evaluate local impact of future scenarios with high penetration of electromobility and renewable energy
Bai et al. Reliability Assessment of Hydro Dominant Systems with Daily and Monthly Inflow Management
Hu et al. Analysis on flexibility resources of Western Inner Mongolia power grid

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210621

Address after: No. 110, Mojiang Road, Yangjiang City, Guangdong Province, 529500

Patentee after: YANGJIANG POWER SUPPLY BUREAU OF GUANGDONG POWER GRID Corp.

Patentee after: ELECTRIC POWER DISPATCHING CONTROL CENTER OF GUANGDONG POWER GRID Co.,Ltd.

Address before: No. 110, Mojiang Road, Jiangcheng District, Yangjiang City, Guangdong Province, 529500

Patentee before: YANGJIANG POWER SUPPLY BUREAU OF GUANGDONG POWER GRID Corp.

Patentee before: Chongqing University

Patentee before: ELECTRIC POWER DISPATCHING CONTROL CENTER OF GUANGDONG POWER GRID Co.,Ltd.