一种基于时序蒙特卡罗模拟的新能源接纳能力评估方法
【专利摘要】本发明是一种基于时序蒙特卡罗模拟的新能源接纳能力评估方法。本发明提出了接纳能力分析模型，可分析新能源接纳能力的限制因素，如调峰调频，线路潮流约束等，根据不同限制因素导致的新能源弃电量，可识别系统的薄弱环节。本发明建立了基于时序蒙特卡罗模拟的新能源接纳能力评估框架，考虑了运行方式变化，负荷波动，新能源波动，水电出力安排，元件故障，机组检修等系统中随时间变化的因素，可为规划人员提供新能源弃电量，新能源发电量，新能源接纳能力值等评价指标，本发明按照充分利用水电水量的原则对水电出力进行安排，充分发挥水电的调峰作用，相比之前的评估方法更加全面。
【专利说明】一种基于时序蒙特卡罗模拟的新能源接纳能力评估方法
【技术领域】
[0001]本发明属于新能源接纳能力评估领域，尤其是一种从规划角度评估主网对大规模新能源的接纳能力的基于时序蒙特卡罗模拟的新能源接纳能力评估方法。
【背景技术】
[0002]发展清洁的可再生新源是解决环境问题和能源制约问题的有效解决方法之一，也是人类社会实现可持续发展的必经之路。近几年，新能源发电发展迅速，装机容量快速上升。因为新能源发电绿色环保且可再生，所以电网应尽可能的消纳新能源的功率和电量。
[0003]新能源普遍具有能量密度低、分布广、间歇性和波动性强等特点，这给电网的规划运行了带来巨大挑战。对主网来说，大规模的风、光电厂接入系统会引起调频问题、调峰问题、电力电量消纳问题和系统潮流大幅变化等问题。因此有必要研究如何评估电网对新能源的接纳能力。
[0004]目前对新能源接纳能力的评估方法主要包括基于电力系统调峰能力的接纳能力评估、基于经济调度的接纳评估、基于稳定的系统接纳能力评估等。但这些方法只考虑了部分影响因素，对系统运行方式时序变化没有计及，因而有待进一步完善。
【发明内容】
[0005]针对现有评估方法的不足，本发明的目的在于提出一种基于时序蒙特卡罗模拟的新能源接纳能力评估方法。本发明是一种考虑系统运行方式时序变化的新能源接纳能力评估方法，可分析设备检修，元件故障，负荷波动，新能源出力的随机性，网架结构及调峰调频约束对新能源接纳能力的影响，可分析不同约束对系统接纳新能源的约束作用，识别系统接纳新能源的薄弱环节，并提出新能源接纳能力的评估指标，为电力系统规划提供指导。
[0006]本发明的基于时序蒙特卡罗模拟的新能源接纳能力评估方法，包括以下步骤:
[0007]I)输入负荷、线路、发电机组信息等原始数据，定义最大模拟年份Nmax与收敛标准；
[0008]2)模拟年次数K=I;
[0009]3)根据等备用原则安排机组检修，抽样得到以下元件的时序状态，包括时序负荷、火电机组和线路的状态持续时间、新能源时序出力，按照充分利用水电水量的原则安排水电的时序出力；
[0010]4)模拟小时数T = O;
[0011]5)对时刻T的系统状态进行评估，求解接纳能力分析模型，如果模型存在最优解，则统计可接纳的新能源出力与对应成本，进入步骤7)，否则进入步骤6)；
[0012]6)判断限制新能源接纳的因素；
[0013]7)如果 T < 8760，T = T+1，进入步骤 5)，否则 K = K+1 ；
[0014]8)如果K等于最大模拟年份Nmax,而且每月的新能源接纳能力已满足收敛条件，进入步骤9)，否则继续进入下一年的状态评估并进入步骤3)；[0015]9)统计并输出所有新能源接纳能力的计算指标。
[0016]本发明建立了基于时序蒙特卡罗模拟的新能源接纳能力评估框架，考虑了运行方式变化，负荷波动，新能源波动，水电出力安排，元件故障，机组检修等系统中随时间变化的因素，可为规划人员提供新能源弃电量，新能源发电量，新能源接纳能力值等评价指标，相比之前的评估方法更加全面。本发明按照充分利用水电水量的原则对水电出力进行安排，充分发挥水电的调峰作用。本发明提出了接纳能力分析模型，可分析新能源接纳能力的限制因素，如调峰调频，线路潮流约束等，根据不同限制因素导致的新能源弃电量，可识别系统的薄弱环节。
【专利附图】
【附图说明】
[0017]图1基于时序蒙特卡罗模拟的新能源接纳能力评估框架。
[0018]图2水电位置安排示意图。
[0019]图3水电出力安排过程。
[0020]图4系统状态图。 [0021]图5新能源接纳能力分布柱状图。
[0022]图6新能源接纳能力累积概率分布。
【具体实施方式】
[0023]本发明是一种基于时序蒙特卡罗模拟的新能源接纳能力评估方法，包括状态选择，状态评估，结果统计三部分，具体分为以下步骤:
[0024]?状态选择:对机组检修进行安排，形成系统不同元件的时序出力模型，包括时序负荷出力，线路运行/故障状态序列，新能源时序出力，水电机组时序出力，火电机组开机状态。其中水电出力按照充分利用水电水量的原则进行安排。
[0025]?状态评估:在不同元件时序出力的基础上组合而成系统状态序列，对不同的系统时序状态执行新能源接纳能力分析模型进行评估，目标为运行成本最小，变量为火电机组出力，约束包括调峰约束、调频约束、线路潮流约束、功率平衡约束。如果评估的结果是系统不能接纳现有新能源出力，则依次去除不同的约束以判断系统限制新能源接纳的因素。
[0026]?结果统计:统计新能源接纳能力的评估指标，包括不同限制因素导致的新能源弃电量，新能源发电量，新能源接纳能力值等。
[0027]本发明新能源接纳能力的评估流程如图1所示，具体步骤如下:
[0028]I)输入负荷、线路、发电机组信息等原始数据，定义最大模拟年份Nmax与收敛标准；
[0029]2)模拟年次数K=I;
[0030]3)根据等备用原则安排机组检修，抽样得到以下元件的时序状态，包括时序负荷、火电机组和线路的状态持续时间、新能源时序出力，按照充分利用水电水量的原则安排水电的时序出力。水电出力的安排过程如下:
[0031]本发明中将水电机组的出力分为两部分，一部分为强迫出力，即水电站为保证下游供水，所要求的一定水流量下水电机组的出力，另一部分为水电站的可调出力，第一部分出力可承担系统的基荷，对第二部分出力采用二分搜索法寻找其在系统时序负荷曲线上的工作位置和工作容量，为体现水电对新能源接纳的作用，此处使用的时序负荷曲线为实际时序负荷曲线与新能源的时序出力曲线的差值，并将其定义为等效时序负荷曲线。每次安排完一台水电机组出力后，则对等效时序负荷曲线进行修正，从系统时序负荷曲线上扣除上一台水电机组的工作出力。每台机组详细的安排过程介绍如下。
[0032]假设我们已经按一定的顺序在系统等效时序负荷曲线图上安排了前i_l个水电机组的工作位置，并且从系统时序负荷曲线上扣除了前i_l个水电机组工作出力后系统的修正时序负荷曲线如图2所示，其中:Ρ^系统修正时序负荷，Ptax-负荷的最大值，P^a-水电机组在系统修正时序负荷曲线上的工作位置，P0ut-水电机组实际最大出力。
[0033]若已知第i个水电机组的可调出力和一段时期内的可调电量分别为PH,pEy，则采用二分法确定水电机组i在系统修正时序负荷曲线上的工作位置的步骤如下:
[0034]a)令水电机组i在系统修正时序负荷曲线上的工作位置Pkxa，上限Xi2及下限Xil的初值为
[0035]PLoca = PLmax，Xj2 = PL—，Xii = O (11)
[0036]其中Ptaax为系统某段时间内的最大负荷。
[0037]其工作出力的初值PQut，上限Pi2及下限Pil为
[0038]P0ut — PH, i，Pi2 — Ph, i，?ii — O (12)
[0039]b)计算修正时序负荷曲线图上Pkjca到P^-Ptxit之间的面积Ei,即水电机组i工作于修正时序负荷曲线上Xil到Xi2位置时的发电量，若Ei > Eh,i时转到步骤c)，若Ei < Eh,i;则转到步骤d)；
[0040]c)若Ei与Eh,i之间的差值小于一定范围，则认为结果收敛，转到步骤e)，否则若Ei> Eh, i，则令 Pi2 = P0ut，P0ut = (Pn+Pi2) /2，若 Ei < Eh, i，则 Pil = P0ut，P0ut = (Pn+Pi2) /2，计算修正时序负荷曲线图上P^a到P^a-Ptxit之间的面积Ei,回到步骤c)；
[0041]d)若Ei与Eh, i之间的差值小于一定范围，则认为结果收敛，转到步骤e)，若Ei >Bh, i，则 Xil = PLocaJ Ploc3 = (XiI+Xi2) /2，若 Ei < Eh, i，则 Xi2 = PLoca, PLoca = (Xn+Xi2) /2，计算修正时序负荷曲线图上P^a到Ρ_-Ρ_之间的面积Ei,回到步骤d)；
[0042]e)得到水电机组相应的工作位置及时序水电机组出力。
[0043]按照此流程可逐个安排水电机组的出力，考虑到水电电量的预测准确性问题，水电的出力过程安排是以月或周为单位。
[0044]4)模拟小时数T = O;
[0045]5)进行状态评估，求解接纳能力分析模型，如果模型存在最优解，则统计可接纳的新能源出力与对应成本，进入步骤7)，否则进入步骤6)。
[0046]接纳能力分析模型的目标函数为系统运行成本和新能源弃电量最小，即
[0047]
【权利要求】
1.一种基于时序蒙特卡罗模拟的新能源接纳能力评估方法，其特征在于包括以下步骤: 1)输入负荷、线路、发电机组信息原始数据，定义最大模拟年份Nmax与收敛标准； 2)模拟年次数K=I; 3)根据等备用原则安排机组检修，抽样得到以下元件的时序状态，包括时序负荷、火电机组和线路的状态持续时间、新能源时序出力，按照充分利用水电水量的原则安排水电的时序出力； 4)模拟小时数T= O5 5)对时刻T的系统状态进行评估，求解接纳能力分析模型，如果模型存在最优解，则统计可接纳的新能源出力与对应成本，进入步骤7)，否则进入步骤6)； 6)判断限制新能源接纳的因素； 7)如果T< 8760，T = T+1，进入步骤5)，否则K = K+1 ； 8)如果K等于最大模拟年份Nmax,而且每月的新能源接纳能力已满足收敛条件，进入步骤9)，否则继续进入下一年的状态评估并进入步骤3)； 9)统计并输出所有新能源接纳能力的计算指标。
2.根据权利要求1所述基于时序蒙特卡罗模拟的新能源接纳能力评估方法，其特征在于考虑电力系统运行方式的时序变化，能分析设备检修，元件故障，负荷波动，新能源出力的随机性，网架结构及调峰调频约束对新能源接纳能力的影响。
3.根据权利要求1所述基于时序蒙特卡罗模拟的新能源接纳能力评估方法，其特征在于将水电机组的出力分为两部分，一部分为强迫出力，用来承担系统的基荷；另一部分为水电站的可调出力，采用二分搜索法寻找其在系统时序负荷曲线上的工作位置和工作容量，为体现水电对新能源接纳的作用，此处使用的时序负荷曲线为实际时序负荷曲线与新能源的时序出力曲线的差值，并将其定义为等效时序负荷曲线，每次安排完一台水电机组出力后，则对等效时序负荷曲线进行修正，从系统时序负荷曲线上扣除上一台水电机组的工作出力。
4.根据权利要求1所述基于时序蒙特卡罗模拟的新能源接纳能力评估方法，其特征在于对系统进行状态评估时，求解接纳能力分析模型，接纳能力分析模型的目标函数为系统运行成本最小和新能源弃电量最小，约束包括:节点功率平衡约束，调峰约束，调频约束，线路潮流约束，新能源弃电量的约束。
5.根据权利要求4所述基于时序蒙特卡罗模拟的新能源接纳能力评估方法，其特征在于上述系统运行成本和切新能源量最小的目标函数为:
N
1*1 式中: N-火电机组数； Pi, t-火电机组i在t时段的出力； Cpii(Piit)-火电机组i的运行费用；一般机组运行费用取为功率的二次函数形式，如下式所示:
6.根据权利要求1所述基于时序蒙特卡罗模拟的新能源接纳能力评估方法，其特征在于接纳能力分析模型的最优解保证同时取得成本最小与新能源出力切除功率最小，如果切除功率为0，模型的最优解即为系统最小的运行成本；如果切除功率大于0，假设存在一个更大的运行成本与更小的切除功率；在切除功率更小时，火电承担的负荷更小，将导致运行成本的降低，因此与假设矛盾，所以接纳能力分析模型的最优解能同时保证取得成本与新能源切除功率最小。
7.根据权利要求1所述基于时序蒙特卡罗模拟的新能源接纳能力评估方法，其特征在于:判断限制新能源接纳的因素，包括线路潮流约束，爬坡率约束，调峰能力约束；判定的方法是依次消除对应的约束方程，然后对接纳能力分析模型求解，若在去除某项约束后，模型存在最优解，则 可判定这项约束为当前系统状态下新能源接纳能力的主要约束原因。
8.根据权利要求1所述基于时序蒙特卡罗模拟的新能源接纳能力评估方法，其特征在于提出评估新能源接纳能力的指标，包括新能源接纳能力值的分布，不同原因导致的新能源弃电量，系统接纳新能源的发电量。
【文档编号】H02J3/28GK104037790SQ201410219466
【公开日】2014年9月10日 申请日期:2014年5月22日 优先权日:2014年5月22日
【发明者】金小明, 别朝红, 李蓉蓉, 宁光涛, 孙灿, 杨柳, 胡源, 张东辉, 高玉洁 申请人:南方电网科学研究院有限责任公司, 西安交通大学, 海南电网公司