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

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

Info

Publication number
CN110224393A
CN110224393A CN201910438343.9A CN201910438343A CN110224393A CN 110224393 A CN110224393 A CN 110224393A CN 201910438343 A CN201910438343 A CN 201910438343A CN 110224393 A CN110224393 A CN 110224393A
Authority
CN
China
Prior art keywords
power
typical scene
node
typical
new energy
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
Application number
CN201910438343.9A
Other languages
English (en)
Other versions
CN110224393B (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

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48—Controlling the sharing of the in-phase component
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, 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。
每个典型场景表示如下：
式中，和分别为高峰、平衡和低谷时段内的第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如下所示：
式中，s为典型场景的编号。Ns为典型场景的总数。ns为典型场景s所表示的典型场景数量。N为典型日中的典型场景总数。Ts为典型场景s所属的典型日在一年中所代表的天数。
3)对电力系统的水电站进行优化调度，使电力系统在正常运行状态下满足正常运行约束条件。
所述正常运行约束条件主要包括风电与光伏发电的消纳约束、水电的消纳约束、正常运行功率平衡约束、线路传输容量约束、机组出力约束、水电站库容约束、水电站发电流量约束和水电站水量平衡约束。
风电与光伏发电的消纳约束如下所示：
式中，Nw和NPV分别为电力系统并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。和分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的出力情况。和分别为第s个典型场景下风电场和光伏电站总的最大出力。
水电的消纳约束如下所示：
式中，Nstop、Nslev和Nslow分别为典型日内高峰平衡低谷时段内的典型场景总数。Nh为并网的水电站总数。sh为水电站的编号。和分别为高峰平衡低谷各个时段内的第s个典型场景下的第sh座水电站的出力情况。 和分别为高峰平衡低谷各个时段内的第s个典型场景所持续的时间。为水电站在任意典型日的最大发电量。
正常运行功率平衡约束如下所示：
其中，和分别为节点i在第s个典型场景下的火电、光伏、风电和水电的有功注入功率。和分别为节点i在第s个典型场景下的有功负荷和注入功率。和分别为节点i与节点j在第s个典型场景下的电压相角。和分别是以为支路导纳建立起来的节点导纳矩阵的自导纳和互导纳。
正常运行时第s个典型场景的节点i的自导纳如下所示：
正常运行时第s个典型场景的节点i和节点j互导纳如下所示：
线路传输容量约束如下所示：
式中，Pij.max为支路ij之间的传输极限。为支路ij在第s个典型场景下所传输的有功功率。
支路ij在第s个典型场景下所传输的有功功率如下所示：
机组出力约束分别如公式10至公式13所示：
式中，和分别为节点i所连火电机组在第s个典型场景下的有功出力上下限。和分别为节点i所连光伏电站在第s个典型场景下的有功出力上下限。和分别为节点i所连风电场在第s个典型场景下的有功出力上下限。和分别为节点i所连接的水电站在第s个典型场景时的有功出力上下限。
水电站发电流量约束如下所示：
式中，和分别为节点i所连接的水电站在第s个典型场景时所分配的发电流量上下限。
发电流量和水电出力的转化关系如下所示：
式中，ηi为节点i所连接的水电站的发电效率。为节点i所连接的水电站在第s个典型场景下的发电流量。为节点i所连接的水电站在第s个典型场景下的发电水头平均高度。
水电站水量平衡约束如下所示：
其中，和分别为节点i所连接的水电站在第s个典型场景初时和末时的水库蓄水量。和分别为节点i所连接的水电站在第s个典型场景时的入库流量、发电流量和其它流量。所述其它流量主要包括灌溉流量和生活用水流量。
水电站库容约束如下所示：
式中，和分别为节点i所连接的水电站在第s个典型场景时的蓄水量上下限。为节点i所连接的水电站在第s个典型场景时的蓄水量。
4)基于N-1预想故障集，建立最小切负荷模型。将电力系统数据输入最小切负荷模型中，得到最小负荷削减量的期望值，使电力系统在N-1预想故障状态下满足故障状态约束条件。
所述故障状态约束条件主要包括风电和光伏发电的消纳约束、水电的消纳约束、故障状态功率平衡约束、线路传输容量约束、机组出力约束、故障状态节点负荷削减量约束、水电站库容约束、水电站发电流量约束、水电站水量平衡约束。
故障状态功率平衡约束如下所示：
式中，和分别为节点i在第f个故障状态下的火电、光伏、风电和水电机组的有功注入功率。为节点i在第f个故障状态下的有功负荷。为节点i在第f个故障状态下的有功负荷削减量。和分别为节点i与节点j在第f个故障状态下的电压相角。和分别是在第f个故障状态下以为支路导纳建立起来的节点导纳矩阵的自导纳及互导纳。
第f个故障状态下节点i的自导纳如下所示：
第f个故障状态下节点i和节点j的互导纳如下所示：
故障状态节点负荷削减量约束如下所示：
5)改变电力系统中新能源电场的装机容量，并重复步骤3和步骤4，计算不同新能源渗透率下的最小负荷削减量的期望值，从而计算得到新能源渗透率上限和新能源消纳能力。
最小切负荷模型如下所示：
式中，s、f和i分别为典型场景、N-1预想故障状态和电力网络节点的编号。Ns、Nf和Nn分别为典型场景、N-1预想故障状态和电力网络节点的数量。ps为第s个场景发生的概率。为第s个场景下的第i个节点在第f个故障状态下的有功负荷削减量。
6)根据不同新能源渗透率下的最小负荷削减量的期望值，建立新能源渗透率和最小负荷削减量期望值的关系曲线，从而判断新能源消纳时系统的安全稳定性。电力系统安全稳定的判断标准为：新能源渗透率ξ＜ξmax。
建立新能源渗透率和最小负荷削减量期望值关系的主要步骤如下：
6.1)改变风电场和光伏电站的装机容量，从而改变电力系统中的新能源渗透率。新能源渗透率ξ如下所示：
式中，Nw和NPV分别为并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。和分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的装机容量。PLD.max为电力系统中负荷功率的最大值。
6.2)利用最小切负荷模型计算不同新能源渗透率下的最小负荷削减量的期望值。
6.3)当电力系统无法满足步骤3或步骤4中的任一约束条件时，则输出电力系统可接受的新能源渗透率上限ξmax。电力系统的新能源消纳能力M如下所示：
式中，M为年度新能源消纳电量。Nw和NPV分别为并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。和分别为第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。
每个典型场景表示如下：
式中，和分别为高峰、平衡和低谷时段内的第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如下所示：
式中，s为典型场景的编号。Ns为典型场景的总数，Ns个典型场景是由所选取的典型日中的原始场景通过K均值聚类技术生成的，系统中的所有原始场景均是通过Ns个典型场景进行表示。ns为典型场景s所表示的典型场景数量。N为典型日中的典型场景总数。Ts为典型场景s所属的典型日在一年中所代表的天数。
3)对电力系统的灵活性电源蓄水式水电站进行优化调度，使电力系统在正常运行状态下满足正常运行约束条件。
所述正常运行约束条件主要包括风电与光伏发电的消纳约束、水电的消纳约束、正常运行功率平衡约束、线路传输容量约束、机组出力约束、水电站库容约束、水电站发电流量约束和水电站水量平衡约束。
风电与光伏发电的消纳约束如下所示：
式中，Nw和NPV分别为电力系统并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。和分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的出力情况。和分别为第s个典型场景下风电场和光伏电站总的最大出力。为了防止出现弃风弃光问题，必须使各场景下系统所消纳的风电和光伏发电等于风电场和光伏电站的最大出力。
水电的消纳约束如下所示：
式中，Nstop、Nslev和Nslow分别为典型日内高峰平衡低谷时段内的典型场景总数。Nh为并网的水电站总数。sh为水电站的编号。和分别为高峰平衡低谷各个时段内的第s个典型场景下的第sh座水电站的出力情况。和分别为高峰平衡低谷各个时段内的第s个典型场景所持续的时间。为水电站在任意典型日的最大发电量。蓄水式水电站具有较大容量的水库，当风电场和光伏电站出力较大时，水库可以将该场景下的多余发电用水量进行储存，用于其它场景时的水力发电。
正常运行功率平衡约束如下所示：
其中，和分别为节点i在第s个典型场景下的火电、光伏、风电和水电的有功注入功率。和分别为节点i在第s个典型场景下的有功负荷和注入功率。和分别为节点i与节点j在第s个典型场景下的电压相角。和分别是以为支路导纳建立起来的节点导纳矩阵的自导纳和互导纳。
正常运行时第s个典型场景的节点i的自导纳如下所示：
正常运行时第s个典型场景的节点i和节点j互导纳如下所示：
线路传输容量约束如下所示：
式中，Pij.max为支路ij之间的传输极限。为支路ij在第s个典型场景下所传输的有功功率。
支路ij在第s个典型场景下所传输的有功功率如下所示：
机组出力约束分别如公式10至公式13所示：
式中，和分别为节点i所连火电机组在第s个典型场景下的有功出力上下限。和分别为节点i所连光伏电站在第s个典型场景下的有功出力上下限。和分别为节点i所连风电场在第s个典型场景下的有功出力上下限。和分别为节点i所连接的水电站在第s个典型场景时的有功出力上下限。
水电站发电流量约束如下所示：
式中，和分别为节点i所连接的水电站在第s个典型场景时所分配的发电流量上下限。水库承担着一定量的防洪、农业灌溉和生活用水等任务，水电站需要根据水库调度部门调度计划进行安排发电用水量，可用水电资源应为调度部门分配的发电流量上下限。
发电流量和水电出力的转化关系如下所示：
式中，ηi为节点i所连接的水电站的发电效率。为节点i所连接的水电站在第s个典型场景下的发电流量。为节点i所连接的水电站在第s个典型场景下的发电水头平均高度。
水电站水量平衡约束如下所示：
其中，和分别为节点i所连接的水电站在第s个典型场景初时和末时的水库蓄水量。和分别为节点i所连接的水电站在第s个典型场景时的入库流量、发电流量和其它流量。所述其它流量主要包括灌溉流量和生活用水流量。
水电站库容约束如下所示：
式中，和分别为节点i所连接的水电站在第s个典型场景时的蓄水量上下限。为节点i所连接的水电站在第s个典型场景时的蓄水量。
4)基于N-1预想故障集，建立最小切负荷模型。将电力系统数据输入最小切负荷模型中，得到最小负荷削减量的期望值，使电力系统在N-1预想故障状态下满足故障状态约束条件。
所述故障状态约束条件主要包括风电和光伏发电的消纳约束、水电的消纳约束、故障状态功率平衡约束、线路传输容量约束、机组出力约束、故障状态节点负荷削减量约束、水电站库容约束、水电站发电流量约束、水电站水量平衡约束。
故障状态功率平衡约束如下所示：
式中，和分别为节点i在第f个故障状态下的火电、光伏、风电和水电机组的有功注入功率。为节点i在第f个故障状态下的有功负荷。为节点i在第f个故障状态下的有功负荷削减量。和分别为节点i与节点j在第f个故障状态下的电压相角。和分别是在第f个故障状态下以为支路导纳建立起来的节点导纳矩阵的自导纳及互导纳。
第f个故障状态下节点i的自导纳如下所示：
式中，为第f个故障状态下支路导纳。
第f个故障状态下节点i和节点j的互导纳如下所示：
故障状态节点负荷削减量约束如下所示：
5)改变电力系统中新能源电场的装机容量，并重复步骤3和步骤4，计算不同新能源渗透率下的最小负荷削减量的期望值，从而计算得到新能源渗透率上限和新能源消纳能力。
最小切负荷模型如下所示：
式中，s、f和i分别为典型场景、N-1预想故障状态和电力网络节点的编号。Ns、Nf和Nn分别为典型场景、N-1预想故障状态和电力网络节点的数量。ps为第s个场景发生的概率。为第s个场景下的第i个节点在第f个故障状态下的有功负荷削减量。
6)根据不同新能源渗透率下的最小负荷削减量的期望值，建立新能源渗透率和最小负荷削减量期望值的关系曲线，从而判断新能源消纳时系统的安全稳定性。电力系统安全稳定的判断标准为：新能源渗透率ξ＜ξmax。
建立新能源渗透率和最小负荷削减量期望值关系曲线的主要步骤如下：
6.1)结合典型场景中的负荷功率，改变风电场和光伏电站的装机容量，从而改变电力系统中的新能源渗透率。新能源渗透率ξ如下所示：
式中，Nw和NPV分别为并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。和分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的装机容量。PLD.max为电力系统中负荷功率的最大值。
6.2)利用最小切负荷模型计算不同新能源渗透率下的最小负荷削减量的期望值。
6.3)当电力系统无法满足步骤3或步骤4中的任一约束条件时，则输出电力系统可接受的新能源渗透率上限ξmax。电力系统的新能源消纳能力M如下所示：
式中，M为年度新能源消纳电量。Nw和NPV分别为并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。和分别为第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。
每个典型场景表示如下：
式中，和分别为高峰、平衡和低谷时段内的第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如下所示：
式中，s为典型场景的编号。Ns为典型场景的总数。ns为典型场景s所表示的典型场景数量。N为典型日中的典型场景总数。Ts为典型场景s所属的典型日在一年中所代表的天数。
实施例4：
一种基于最小切负荷模型的新能源消纳评估方法，主要步骤参见实施例2，其中，正常运行约束条件主要包括风电与光伏发电的消纳约束、水电的消纳约束、正常运行功率平衡约束、线路传输容量约束、机组出力约束、水电站库容约束、水电站发电流量约束和水电站水量平衡约束。
风电与光伏发电的消纳约束如下所示：
式中，Nw和NPV分别为电力系统并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。和分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的出力情况。和分别为第s个典型场景下风电场和光伏电站总的最大出力。
水电的消纳约束如下所示：
式中，Nstop、Nslev和Nslow分别为典型日内高峰平衡低谷时段内的典型场景总数。Nh为并网的水电站总数。sh为水电站的编号。和分别为高峰平衡低谷各个时段内的第s个典型场景下的第sh座水电站的出力情况。和分别为高峰平衡低谷各个时段内的第s个典型场景所持续的时间。为水电站在任意典型日的最大发电量。
正常运行功率平衡约束如下所示：
其中，和分别为节点i在第s个典型场景下的火电、光伏、风电和水电的有功注入功率。和分别为节点i在第s个典型场景下的有功负荷和注入功率。和分别为节点i与节点j在第s个典型场景下的电压相角。和分别是以为支路导纳建立起来的节点导纳矩阵的自导纳和互导纳。
正常运行时第s个典型场景的节点i的自导纳如下所示：
正常运行时第s个典型场景的节点i和节点j互导纳如下所示：
线路传输容量约束如下所示：
式中，Pij.max为支路ij之间的传输极限。为支路ij在第s个典型场景下所传输的有功功率。
支路ij在第s个典型场景下所传输的有功功率如下所示：
机组出力约束分别如公式8至公式11所示：
式中，和分别为节点i所连火电机组在第s个典型场景下的有功出力上下限。和分别为节点i所连光伏电站在第s个典型场景下的有功出力上下限。和分别为节点i所连风电场在第s个典型场景下的有功出力上下限。和分别为节点i所连接的水电站在第s个典型场景时的有功出力上下限。
水电站发电流量约束如下所示：
式中，和分别为节点i所连接的水电站在第s个典型场景时所分配的发电流量上下限。
发电流量和水电出力的转化关系如下所示：
式中，ηi为节点i所连接的水电站的发电效率。为节点i所连接的水电站在第s个典型场景下的发电流量。为节点i所连接的水电站在第s个典型场景下的发电水头平均高度。
水电站水量平衡约束如下所示：
其中，和分别为节点i所连接的水电站在第s个典型场景初时和末时的水库蓄水量。和分别为节点i所连接的水电站在第s个典型场景时的入库流量、发电流量和其它流量。所述其它流量主要包括灌溉流量和生活用水流量。
水电站库容约束如下所示：
式中，和分别为节点i所连接的水电站在第s个典型场景时的蓄水量上下限。为节点i所连接的水电站在第s个典型场景时的蓄水量。
实施例5：
一种基于最小切负荷模型的新能源消纳评估方法，主要步骤参见实施例2，其中，故障状态约束条件主要包括风电和光伏发电的消纳约束、水电的消纳约束、故障状态功率平衡约束、线路传输容量约束、机组出力约束、故障状态节点负荷削减量约束、水电站库容约束、水电站发电流量约束、水电站水量平衡约束。
故障状态功率平衡约束如下所示：
式中，和分别为节点i在第f个故障状态下的火电、光伏、风电和水电机组的有功注入功率。为节点i在第f个故障状态下的有功负荷。为节点i在第f个故障状态下的有功负荷削减量。和分别为节点i与节点j在第f个故障状态下的电压相角。和分别是在第f个故障状态下以为支路导纳建立起来的节点导纳矩阵的自导纳及互导纳。
第f个故障状态下节点i的自导纳如下所示：
第f个故障状态下节点i和节点j的互导纳如下所示：
故障状态节点负荷削减量约束如下所示：
风电与光伏发电的消纳约束如下所示：
式中，Nw和NPV分别为电力系统并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。和分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的出力情况。和分别为第s个典型场景下风电场和光伏电站总的最大出力。
水电的消纳约束如下所示：
式中，Nstop、Nslev和Nslow分别为典型日内高峰平衡低谷时段内的典型场景总数。Nh为并网的水电站总数。sh为水电站的编号。和分别为高峰平衡低谷各个时段内的第s个典型场景下的第sh座水电站的出力情况。和分别为高峰平衡低谷各个时段内的第s个典型场景所持续的时间。为水电站在任意典型日的最大发电量。
线路传输容量约束如下所示：
式中，Pij.max为支路ij之间的传输极限。为支路ij在第s个典型场景下所传输的有功功率。
支路ij在第s个典型场景下所传输的有功功率如下所示：
机组出力约束分别如公式9至公式12所示：
式中，和分别为节点i所连火电机组在第s个典型场景下的有功出力上下限。和分别为节点i所连光伏电站在第s个典型场景下的有功出力上下限。和分别为节点i所连风电场在第s个典型场景下的有功出力上下限。和分别为节点i所连接的水电站在第s个典型场景时的有功出力上下限。
水电站发电流量约束如下所示：
式中，和分别为节点i所连接的水电站在第s个典型场景时所分配的发电流量上下限。
发电流量和水电出力的转化关系如下所示：
式中，ηi为节点i所连接的水电站的发电效率。为节点i所连接的水电站在第s个典型场景下的发电流量。为节点i所连接的水电站在第s个典型场景下的发电水头平均高度。
水电站水量平衡约束如下所示：
其中，和分别为节点i所连接的水电站在第s个典型场景初时和末时的水库蓄水量。和分别为节点i所连接的水电站在第s个典型场景时的入库流量、发电流量和其它流量。所述其它流量主要包括灌溉流量和生活用水流量。
水电站库容约束如下所示：
式中，和分别为节点i所连接的水电站在第s个典型场景时的蓄水量上下限。为节点i所连接的水电站在第s个典型场景时的蓄水量。
实施例6：
一种基于最小切负荷模型的新能源消纳评估方法，主要步骤参见实施例2，其中，最小切负荷模型如下所示：
式中，s、f和i分别为典型场景、N-1预想故障状态和电力网络节点的编号。Ns、Nf和Nn分别为典型场景、N-1预想故障状态和电力网络节点的数量。ps为第s个场景发生的概率。为第s个场景下的第i个节点在第f个故障状态下的有功负荷削减量。
实施例7：
一种基于最小切负荷模型的新能源消纳评估方法，主要步骤参见实施例2，其中，建立新能源渗透率和最小负荷削减量期望值关系的主要步骤如下：
6.1)改变风电场和光伏电站的装机容量，从而改变电力系统中的新能源渗透率。新能源渗透率ξ如下所示：
式中，Nw和NPV分别为并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。和分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的装机容量。PLD.max为电力系统中负荷功率的最大值。
6.2)利用最小切负荷模型计算不同新能源渗透率下的最小负荷削减量的期望值。
6.3)当电力系统无法满足步骤3或步骤4中的任一约束条件时，则输出电力系统可接受的新能源渗透率上限ξmax。电力系统的新能源消纳能力M如下所示：
式中，M为年度新能源消纳电量。Nw和NPV分别为并网的风电场和光伏电站总数。sw和sPV分别为风电场和光伏电站的编号。和分别为第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日高峰时段的典型场景发生的概率
3)对电力系统的水电站进行优化调度，使电力系统在正常运行状态下满足正常运行约束条件。
4)基于N-1预想故障集，建立最小切负荷模型。将电力系统数据输入最小切负荷模型中，利用CPLEX求解器求解得到最小负荷削减量的期望值，使电力系统在N-1预想故障状态下满足故障状态约束条件。
基于测试网络和典型场景数据，在正常运行状态下，实现系统内的新能源全额消纳，并且避免出现线路过载和切负荷等问题；在故障状态下，根据N-1安全校核的预想故障集，计算出各个典型场景下的最小负荷削减量，计算结果如表3所示。
表3各典型场景下的最小负荷削减量(MW)
进一步结合预想故障数量和各个典型场景概率，计算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；
每个典型场景表示如下：
式中，和分别为高峰、平衡和低谷时段内的第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如下所示：
式中，s为典型场景的编号；Ns为典型场景的总数；ns为典型场景s所表示的典型场景数量；N为典型日中的典型场景总数；Ts为典型场景s所属的典型日在一年中所代表的天数。
3.根据权利要求1所述的一种基于最小切负荷模型的新能源消纳评估方法，其特征在于：所述正常运行约束条件主要包括风电与光伏发电的消纳约束、水电的消纳约束、正常运行功率平衡约束、线路传输容量约束、机组出力约束、水电站库容约束、水电站发电流量约束和水电站水量平衡约束；
风电与光伏发电的消纳约束如下所示：
式中，Nw和NPV分别为电力系统并网的风电场和光伏电站总数；sw和sPV分别为风电场和光伏电站的编号；和分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的出力情况；和分别为第s个典型场景下风电场和光伏电站总的最大出力；
水电的消纳约束如下所示：
式中，Nstop、Nslev和Nslow分别为典型日内高峰平衡低谷时段内的典型场景总数；Nh为并网的水电站总数；sh为水电站的编号；和分别为高峰平衡低谷各个时段内的第s个典型场景下的第sh座水电站的出力情况；ΔTs top、ΔTs lev和ΔTs low分别为高峰平衡低谷各个时段内的第s个典型场景所持续的时间；为水电站在任意典型日的最大发电量；
正常运行功率平衡约束如下所示：
其中，PVi s、Wi s和分别为节点i在第s个典型场景下的火电、光伏、风电和水电的有功注入功率；和Pi s分别为节点i在第s个典型场景下的有功负荷和注入功率；和分别为节点i与节点j在第s个典型场景下的电压相角；和分别是以为支路导纳建立起来的节点导纳矩阵的自导纳和互导纳；
正常运行时第s个典型场景的节点i的自导纳如下所示：
正常运行时第s个典型场景的节点i和节点j互导纳如下所示：
线路传输容量约束如下所示：
式中，Pij.max为支路ij之间的传输极限；为支路ij在第s个典型场景下所传输的有功功率；
支路ij在第s个典型场景下所传输的有功功率如下所示：
机组出力约束分别如公式10至公式13所示：
式中，和分别为节点i所连火电机组在第s个典型场景下的有功出力上下限；和分别为节点i所连光伏电站在第s个典型场景下的有功出力上下限；和分别为节点i所连风电场在第s个典型场景下的有功出力上下限；和分别为节点i所连接的水电站在第s个典型场景时的有功出力上下限；
水电站发电流量约束如下所示：
式中，和分别为节点i所连接的水电站在第s个典型场景时所分配的发电流量上下限；
发电流量和水电出力的转化关系如下所示：
式中，ηi为节点i所连接的水电站的发电效率；为节点i所连接的水电站在第s个典型场景下的发电流量；为节点i所连接的水电站在第s个典型场景下的发电水头平均高度；
水电站水量平衡约束如下所示：
其中，Vi s.t和Vi s.t+1分别为节点i所连接的水电站在第s个典型场景初时和末时的水库蓄水量；和分别为节点i所连接的水电站在第s个典型场景时的入库流量、发电流量和其它流量；所述其它流量主要包括灌溉流量和生活用水流量；
水电站库容约束如下所示：
式中，和分别为节点i所连接的水电站在第s个典型场景时的蓄水量上下限；Vi s为节点i所连接的水电站在第s个典型场景时的蓄水量。
4.根据权利要求1或2所述的一种基于最小切负荷模型的新能源消纳评估方法，其特征在于，最小切负荷模型如下所示：
式中，s、f和i分别为典型场景、N-1预想故障状态和电力网络节点的编号；Ns、Nf和Nn分别为典型场景、N-1预想故障状态和电力网络节点的数量；ps为第s个场景发生的概率；为第s个场景下的第i个节点在第f个故障状态下的有功负荷削减量。
5.根据权利要求1或3所述的一种基于最小切负荷模型的新能源消纳评估方法，其特征在于：所述故障状态约束条件主要包括风电和光伏发电的消纳约束、水电的消纳约束、故障状态功率平衡约束、线路传输容量约束、机组出力约束、故障状态节点负荷削减量约束、水电站库容约束、水电站发电流量约束、水电站水量平衡约束；
故障状态功率平衡约束如下所示：
式中，PVi f、Wi f和分别为节点i在第f个故障状态下的火电、光伏、风电和水电机组的有功注入功率；为节点i在第f个故障状态下的有功负荷；为节点i在第f个故障状态下的有功负荷削减量；和分别为节点i与节点j在第f个故障状态下的电压相角；和分别是在第f个故障状态下以为支路导纳建立起来的节点导纳矩阵的自导纳及互导纳；
第f个故障状态下节点i的自导纳如下所示：
第f个故障状态下节点i和节点j的互导纳如下所示：
故障状态节点负荷削减量约束如下所示：
6.根据权利要求1所述的一种基于最小切负荷模型的新能源消纳评估方法，其特征在于，建立新能源渗透率和最小负荷削减量期望值关系的主要步骤如下：
1)改变风电场和光伏电站的装机容量，从而改变电力系统中的新能源渗透率；新能源渗透率ξ如下所示：
式中，Nw和NPV分别为并网的风电场和光伏电站总数；sw和sPV分别为风电场和光伏电站的编号；和分别为第s个典型场景下的第sw座风电场和第sPV座光伏电站的装机容量；PLD.max为电力系统中负荷功率的最大值；
2)利用最小切负荷模型计算不同新能源渗透率下的最小负荷削减量的期望值；
3)当电力系统无法满足步骤3或步骤4中的任一约束条件时，则输出电力系统可接受的新能源渗透率上限ξmax；电力系统的新能源消纳能力M如下所示：
式中，M为年度新能源消纳电量；Nw和NPV分别为并网的风电场和光伏电站总数；sw和sPV分别为风电场和光伏电站的编号；和分别为第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 true CN110224393A (zh) 2019-09-10
CN110224393B 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)

Cited By (12)

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

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
HANDREA BERNANDO TAMBUNAN 等: ""Intermittent Renewable Energy Source (IRES) Penetration Level into Bangka Power System"", 《2018 INTERNATIONAL CONFERENCE ON SMART GREEN TECHNOLOGY IN ELECTRICAL AND INFORMATION SYSTEMS (ICSGTEIS)》 *
彭春华 等: ""基于K- 均值聚类多场景时序特性分析的分布式电源多目标规划"", 《电力自动化设备》 *
马建伟 等: ""计及预测可信度时段差异性的含分布式可再生能源配电网调度策略"", 《电力系统及其自动化学报》 *

Cited By (19)

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

Also Published As

Publication number Publication date
CN110224393B (zh) 2020-08-25

Similar Documents

Publication Publication Date Title
CN110224393A (zh) 一种基于最小切负荷模型的新能源消纳评估方法
Mc Garrigle et al. How much wind energy will be curtailed on the 2020 Irish power system?
Cheng et al. China׳ s small hydropower and its dispatching management
CN102427244B (zh) 大规模光伏风电信息接入系统
CN110909954B (zh) 一种最大化可再生能源利用的多阶段电源规划方法
CN103577901A (zh) 潮间带风电接入电网的方法
CN104578183B (zh) 提高电力系统低碳化水平的联络线送电计划优化方法
Kizito et al. Optimal distributed generator placement in utility-based microgrids during a large-scale grid disturbance
CN105244922B (zh) 限风情况下双馈型风电场群有功出力优化分配方法
CN111799794B (zh) 一种考虑输电阻塞盈余的输电网扩展规划方法
CN108155674A (zh) 考虑不确定分布特性的水火电联合调度方法及系统
CN117175543A (zh) 一种可调节负荷的配电网规划策略优化方法及系统
Gauché et al. A concentrating solar power value proposition for South Africa
CN110460090B (zh) 一种水风光微电网电源容量配置方法
Peng et al. Research on inter-regional renewable energy accommodation assessment method based on time series production simulation
Zhang et al. Research on active distribution network structure planning for multi distributed generation access
CN112184016A (zh) 复杂电网下便捷判断水光互补一体化光伏规模的方法
Chi et al. A case study of developing an intelligent management system for energy internet
CN106485593A (zh) 一种电网侧新能源发电受限原因分析方法
CN108736521A (zh) 一种具有电网接入控制功能的风电场能量管理系统
CN104993526A (zh) 多类型电源联合运行系统装机配比计算方法
Cui et al. Review of transmission planning with large-scale wind power integration
Kang et al. New Measurement Method of Distributed Power Generation
CN113162101B (zh) 一种含互斥约束的交直流受端电网优化规划方法及系统
CN110880793B (zh) 一种基于日流量的小水电微电网电源容量配置方法

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.