重庆大学学报  2022, Vol. 45 Issue (1): 9-17  DOI: 10.11835/j.issn.1000-582X.2020.229 RIS(文献管理工具)
0

引用本文 

谭涛亮, 姜云鹏, 李震, 钱峰, 任洲洋, 刘俊磊, 颜伟. 考虑灵活性的地区电网新能源消纳能力评估方法[J]. 重庆大学学报, 2022, 45(1): 9-17. DOI: 10.11835/j.issn.1000-582X.2020.229.
TAN Taoliang, JIANG Yunpeng, LI Zhen, QIAN Feng, REN Zhouyang, LIU Junlei, YAN Wei. Assessment method of renewable energy accommodation capacity of regional power grid considering flexibility[J]. Journal of Chongqing University, 2022, 45(1): 9-17. DOI: 10.11835/j.issn.1000-582X.2020.229.

基金项目

广东电网有限责任公司科技资助项目(031700KK52180002)

作者简介

谭涛亮(1980-), 男, 高级工程师, 博士, 主要研究方向为电力调度及分析, (E-mail)30940480@qq.com

文章历史

收稿日期: 2019-12-04
考虑灵活性的地区电网新能源消纳能力评估方法
谭涛亮 1, 姜云鹏 2, 李震 1, 钱峰 3, 任洲洋 2, 刘俊磊 3, 颜伟 2     
1. 广东电网有限责任公司 阳江供电局, 广东 阳江 529500;
2. 重庆大学 输配电装备及系统安全与新技术国家重点实验室, 重庆 400044;
3. 广东电网有限责任公司 电力调度控制中心, 广州 510000
摘要: 在水电富集的地区,水电站作为一种灵活性电源用于提高电网的新能源消纳能力。提出了一种考虑灵活性的地区电网新能源消纳能力评估方法,该方法能够保证系统安全可靠,同时全额消纳新能源。首先,考虑源荷功率的相关性,采用K均值聚类技术生成典型场景;然后,通过灵活性电源的优化调度,使得系统在正常运行状态下全额消纳新能源,并在故障状态下使负荷削减量的期望值最小;最后,选取中国南方某沿海地区110 kV电网进行算例分析,研究新能源渗透率和最小负荷削减量期望值的关系,进而评估系统的新能源消纳能力,验证文中方法的有效性。
关键词: 最小负荷削减量    新能源消纳    新能源渗透率    场景生成    电源灵活性    
Assessment method of renewable energy accommodation capacity of regional power grid considering flexibility
TAN Taoliang 1, JIANG Yunpeng 2, LI Zhen 1, QIAN Feng 3, REN Zhouyang 2, LIU Junlei 3, YAN Wei 2     
1. Yangjiang Power Supply Bureau of Guangdong Power Grid Corporation, Yangjiang 529500, Guangdong, P. R. China;
2. State Key Laboratory of Power Transmission Equipment & System Security and New Technology, Chongqing University, Chongqing 400044, P. R. China;
3. Power Dispatching and Control Center, Guangdong Grid Co., Guangzhou 510000, P. R. China
Abstract: In some areas with rich water resources, hydropower station is a flexible power source which can be used to promote the renewable energy accommodation capacity of the power grid. This paper proposes a renewable energy accommodation capacity assessment method of regional power grid considering flexibility which can ensure full accommodation for renewable energy in the premise of safety and reliability. First, representative scenarios are generated using K-means clustering considering the correlation between the power sources and loads. Second, by optimizing the flexible generator dispatching, renewable energy can be fully used under normal operation states, and the expectation of load-shedding is minimized under fault conditions. Finally, with the data of a 110 kV regional power grid in a coastal area of South China, the relationship between renewable energy penetration and expectation of minimum load-shedding is captured, and the renewable energy accommodation capability of the selected network is evaluated by the proposed method, verifying the effectiveness of the method.
Keywords: minimum load-shedding    renewable energy accommodation    renewable energy penetration    scenario generation    generation flexibility    

近年来,中国新能源发电已经取得较好的环境效益和经济效益。截至2017年底,中国风电装机容量达到了1.6×08 kW,光伏装机容量达到了1.3×108 kW。然而,弃风、弃光量也分别达到了4.19×1010 kW·h和7.3×109 kW·h[1]。随着新能源装机容量的快速增长,新能源消纳成为学者关注的核心问题之一。南方地区水资源丰富,建立了一定规模的蓄水式水电站,系统具有较高的灵活性。在一定的新能源渗透率下,充分利用系统中的电源灵活性可以实现新能源的全额消纳。新能源消纳需在保证系统安全稳定运行的前提下,合理利用各种资源提高新能源的效益[2],系统不能因消纳新能源,产生大量的切负荷、线路过载等不良运行状态。但风电场、光伏电站的并网运行改变了系统的电源结构,也改变了网络的潮流分布,系统面临着线路过载和切负荷的风险[3]。文中以负荷削减量为系统安全可靠评价指标,提出了基于最小切负荷模型的新能源消纳评估方法,用于分析新能源渗透率和负荷削减量之间的关系,为统筹协调系统的安全可靠运行和新能源消纳提供参考。

当前,国内外关于新能源消纳的研究较多。文献[4-5]基于系统在正常运行状态下的网络安全约束,评估了风电消纳能力。其中,文献[4]评估了系统在特定时间段内的风电消纳能力,但没有考虑系统内灵活性电源的调节能力;文献[5]分析了系统在典型日峰时段内的风电消纳能力,但缺乏系统在极端时段(峰时段或谷时段)的消纳能力数据,不能全面反映系统的风电消纳能力。文献[6-8]考虑了区域内可灵活调度的水电资源,优化水电机组启停状态,并利用风电水电互补特性来提升新能源消纳能力,简化了水电站的运行模型,而蓄水式水电站的水库承担着防洪、灌溉和供水等其他综合性任务,文献没有考虑水库因其他任务需求而限制的机组发电能力。此外,文献[4-8]只考虑了系统在正常运行状态下的各类约束,没有对系统进行N-1安全校核,系统发生故障时可能出现线路过载、切负荷等现象,影响供电的可靠性。

为了评估系统的新能源消纳能力以及分析系统的安全可靠性,文中提出考虑灵活性的地区电网新能源消纳能力评估方法。首先,考虑源荷相关性,基于K均值聚类技术生成若干典型场景;然后,在各个场景下,通过对灵活性电源的优化调度,保证系统在正常运行状态下不出现新能源弃能、线路过载和切负荷现象;在故障情况下,根据N-1安全校核的预想故障集,优化不同新能源渗透率下的最小负荷削减量的期望值;最后,结合南方某沿海地区110 kV电网的实际数据,分析了新能源渗透率和最小负荷削减量期望值的关系,评估了系统的新能源消纳能力,验证了文中方法的有效性。

1 考虑源荷相关性的场景生成法

新能源出力和负荷功率都具有波动性和不确定性,在不同时刻可能存在着明显差异。如果新能源消纳能力的评估时刻过多,则计算量巨大。场景分析法能够分析原始数据,生成多个典型场景,新能源的消纳评估只需对生成的典型场景进行分析计算。此外,电网中的电源出力和负荷功率之间始终保持着对应关系,在对源荷功率进行分析时,选择某个时刻的电源出力就必须选择该时刻对应的负荷功率。因此,在场景生成时需要计及电源出力和负荷功率之间的相关性[9]

K均值聚类技术作为一种经典的场景分析法,既能挖掘不同时刻场景性质的异同,聚类分析生成典型场景,也能保证所生成的典型场景中新能源出力和负荷功率的相关性[10-11]。因此,基于K均值聚类技术进行场景生成。考虑到K均值聚类初始聚类中心的随机选择对聚类结果影响较大,在场景生成时采用密度参数法[12-13]优化初始聚类中心,进而提高聚类精度。基于K均值聚类技术对源荷序列进行典型场景生成主要为以下4个阶段:

第1阶段,考虑不同时段负荷水平的差异性,通过对典型日不同时刻的源荷功率数据采样分析,划分典型日为峰平谷3个时段。

第2阶段,考虑源荷功率的相关性,以采样时刻为序列分别建立峰、平、谷各个时段内的原始场景PttopPtlevPtlow,每个场景可表示为

$ \left\{ {\begin{array}{*{20}{c}} {\boldsymbol{P}_t^{{\rm{top}}} = \left[ {{\rm{ }}P_{1.t}^{{\rm{top}}}, P_{2.t}^{{\rm{top}}}, \ldots , P_{m.t}^{{\rm{top}}}} \right]\;\;\;t = 1, 2, \ldots , {T_{{\rm{top}}}}, }\\ {\boldsymbol{P}_t^{{\rm{lev}}} = {\rm{ }}\left[ {P_{1.t}^{{\rm{lev}}}, P_{2.t}^{{\rm{lev}}}, \ldots , P_{m.t}^{{\rm{lev}}}} \right]\;\;\;t = 1, 2, \ldots , {T_{{\rm{lev}}}}, }\\ {\boldsymbol{P}_t^{{\rm{low}}} = {\rm{ }}\left[ {P_{1.t}^{{\rm{low}}}, P_{2.t}^{{\rm{low}}}, \ldots , P_{m.t}^{{\rm{low}}}} \right]\;\;t = 1, 2, \ldots , {T_{{\rm{low}}}};} \end{array}} \right. $ (1)

式中:${\rm{ }}P_{x.t}^{{\rm{top}}}P_{x.t}^{{\rm{lev}}}{\rm{1}}和P_{x.t}^{{\rm{low}}}$分别为峰平谷时段内的第t个采样时刻在第x个采样点测得的源荷功率,x=1, 2, …, mm为采样点总数;TtopTlevTlow分别为峰平谷时段内的采样时刻总数。

第3阶段,确定峰、平、谷各个时段内的聚类个数KtopKlevKlow,利用密度参数法依次优化生成KtopKlevKlow个初始聚类中心;再计算各个场景与聚类中心的欧式距离,将其分配到距离最近的类中;最后,利用均值法迭代更新聚类中心直至收敛,获得KtopKlevKlow个类,KtopKlevKlow个聚类中心以及场景发生的概率。

第4阶段,考虑聚类结果应和实际系统中的源荷功率场景相吻合,分别将峰平谷各时段内生成的聚类中心同该时段内的原始场景进行比较,寻找最相似的原始场景进行替换,保证生成的典型场景的真实性。

2 基于最小负荷削减量的新能源消纳评估模型

面向中国水电资源富集的南方地区,提出利用系统中蓄水式水电站的灵活性来实现风电、光伏发电和水电全额消纳的方法。由于新能源电源的并网运行改变了电网中的电源结构,系统发生某些故障时,电网可能会出现线路过载、切负荷的情况。因此,含新能源的系统不仅要考虑新能源消纳状况,还需要计及系统在正常状态和故障状态下的安全可靠性。为此,文中以负荷削减量为系统安全可靠评价指标,建立了基于最小负荷削减量的新能源消纳评估模型。

2.1 目标函数

以系统负荷削减量的期望值F最小为优化目标,建立目标函数为

$ {\rm{min}}F = \sum\limits_{s = 1}^{{N_s}} {\left( {\frac{1}{{{N_f}}}\sum\limits_{f = 1}^{{N_f}} {\left( {\sum\limits_{i = 1}^{{N_n}} {P_{{\rm{loss}}.i}^{s.f}} } \right){\rm{ }}} } \right){\rm{ }}{p_s}{\rm{ }}} , $ (2)

式中:sfi分别为场景、故障状态和节点的编号;NsNfNn分别为场景、故障状态和节点的数量;ps为第s个场景发生的概率;${\rm{ }}P_{{\rm{loss}}.i}^{s.f}$为第s个场景下的第i个节点在第f个故障状态下的有功负荷削减量。

2.2 正常运行状态的约束条件 2.2.1 风电和光伏的消纳约束

为了防止出现弃风、弃光问题,必须使各场景下系统所消纳的风电和光伏发电功率等于风电场和光伏电站的最大出力,建立如下约束:

$ \sum\limits_{w = 1}^{{N_w}} {{W_{s.w}} + \sum\limits_{pv = 1}^{{N_{PV}}} {P{V_{s.pv}} = W_s^{{\rm{max}}} + PV_s^{{\rm{max}}}} } , $ (3)

式中:NwNPV分别为并网的风电场和光伏电站总数;wpv分别为风电场和光伏电站的编号;Ws.wPVs.pv分别为第s个场景下的第w座风电场和第pv座光伏电站的出力情况;WsmaxPVsmax分别为第s个场景下风电场和光伏电站最大出力。

2.2.2 水电的消纳约束

蓄水式水电站具有容量较大的水库,当风电场和光伏电站出力较大时,水库可以将该时段多余的发电用水量进行储存,用于其他时段的水力发电,建立的水电电量平衡约束为

$ \sum\limits_{s = 1}^{{N_{{\rm{top}}}}} {\left( {\sum\limits_{h = 1}^{{N_h}} {H_{s.h}^{{\rm{top}}}} } \right)} \Delta T_s^{{\rm{top}}} + \sum\limits_{s = 1}^{{N_{{\rm{lev}}}}} {\left( {\sum\limits_{h = 1}^{{N_h}} {H_{s.h}^{{\rm{lev}}}} } \right)} \Delta T_s^{{\rm{lev}}} + \sum\limits_{s = 1}^{{N_{{\rm{low}}}}} {\left( {\sum\limits_{h = 1}^{{N_h}} {H_{s.h}^{{\rm{low}}}} } \right)} \Delta T_s^{{\rm{low}}} = WH_T^{{\rm{max}}}, $ (4)

式中:NtopNlevNlow分别为典型日峰时段、平时段和谷时段内的场景总数;Nh为并网的水电站总数;h为水电站编号;Hs.htopHs.hlevHs.hlow分别为峰平谷各个时段内的第s个场景的第h座水电站出力情况;ΔTstop、ΔTslev和ΔTslow分别为峰平谷各个时段内的第s个场景所持续的时间;WHTmax为水电站在某个典型日的最大发电量。

2.2.3 功率平衡约束[14]
$ \left\{ {\begin{array}{*{20}{l}} {G_i^s + PV_i^s + W_i^s + H_i^s - P_{{\rm{load}}.i}^s = P_i^s,}\\ {P_i^s = B_{ii}^s\theta _i^s + \sum\limits_{j \in i,j \ne i{\rm{ }}} {B_{ij}^s\theta _j^s} ,} \end{array}} \right. $ (5)

式中:$G_i^s、PV_i^s、W_i^s和H_i^s$分别为节点i在第s个场景下的火电、光伏发电、风电和水电的有功注入功率;Pload.is为节点i在第s个场景下的有功负荷;θisθjs分别为节点i与节点j在第s个场景下的电压相角;BiisBijs分别是以$\frac{1}{{{x_{ij}}}}$为支路导纳建立起来的节点导纳矩阵的自导纳和互导纳,$B_{ii}^s = - \sum\limits_{j \in i, j \ne i} {\frac{1}{{{x_{ij}}}}, B_{ij}^s = - \frac{1}{{{x_{ij}}}}} $

2.2.4 线路传输容量约束
$ - {P_{ij.{\rm{max}}}} \le P_{ij}^s \le {P_{ij.{\rm{max}}}}, $ (6)

式中:Pij.max为支路ij之间的传输极限;Pijs为支路ij在第s个场景下所传输的有功功率,计算公式[15]

$ P_{ij}^s = \frac{{\theta _i^s - \theta _j^s}}{{{x_{ij}}}}。$ (7)
2.2.5 机组出力约束
$ G_{i.{\rm{min}}}^s \le G_i^s \le G_{i.{\rm{max}}}^s, $ (8)
$ PV_{i.{\rm{min}}}^s \le PV_i^s \le PV_{i.{\rm{max}}}^s, $ (9)
$ W_{i.{\rm{min}}}^s \le W_i^s \le W_{i.{\rm{max}}}^s, $ (10)
$ H_{i.{\rm{min}}}^s \le H_i^s \le H_{i.{\rm{max}}}^s, $ (11)

式中:Gi.maxsGi.minsPVi.maxsPVi.minsWi.maxsWi.mins以及Hi.maxsHi.mins分别为节点i所连接的火电站、光伏电站、风电场、水电站在第s个场景下的有功出力上下限。需要说明的是,当风电、光伏发电被全额消纳时,PVis=PVi.maxsWis=Wi.maxs

2.2.6 水电站发电流量约束

水库承担着一定量的防洪、灌溉和生活用水等任务,水电站需要根据水库调度部门安排发电用水量,可用水电资源应为调度部门分配的发电流量上下限,建立如下约束:

$ q_{i.{\rm{min}}}^s \le q_i^s \le q_{i.{\rm{max}}}^s, $ (12)

式中,qi.maxsqi.mins分别为节点i所连接的水电站在第s个场景时所分配的发电流量上下限。发电流量与水电出力的转化关系为[16]

$ H_i^s = 9.81{\eta _i}q_i^sh_i^s, $ (13)

式中:ηi为节点i所连接的水电站的发电效率;qishis分别为节点i所连接的水电站在第s个场景下的发电流量和发电水头平均高度。

2.2.7 水电站水量平衡约束[17]
$ V_i^{s + 1} = V_i^s + (Q_i^s - q_i^s - q_i^{s.l}), $ (14)

式中:VisVis+1分别为节点i所连接的水电站在第s个场景初和末的水库蓄水量;Qisqisqis.l分别为节点i所连接的水电站在第s个场景时的入库流量、发电流量和其他流量(灌溉、生活用水等流量)。

2.2.8 水电站库容约束[18]
$ V_{i.{\rm{min}}}^s \le V_i^s \le V_{i.{\rm{max}}}^s, $ (15)

式中,Vi.maxsVi.mins分别为节点i所连接的水电站在第s个场景时的蓄水量上下限。

2.3 故障状态的约束条件

由于系统在故障情况下可能出现切负荷,并且切负荷会造成系统潮流的改变,因此,需要修改功率平衡约束,增加负荷削减量约束,其他约束条件与式(3)(4)(6)(15)相同。

2.3.1 功率平衡约束
$ \left\{ {\begin{array}{*{20}{l}} {G_i^f + PV_i^f + W_i^f + H_i^f - P_{{\rm{load}}.i}^f + P_{{\rm{loss}}.i}^f = P_i^f,}\\ {P_i^f = B_{ii}^f\theta _i^f + \sum\limits_{j \in i,j \ne i} {B_{ij}^f\theta _j^f} ,} \end{array}} \right. $ (16)

式中:GifPVifWifHif分别为节点i在第f个故障状态下的火电、光伏、风电和水电的有功注入功率;Pload.ifPloss.if分别为节点i在第f个故障状态下的有功负荷和有功负荷削减量;θifθjf分别为节点i与节点j在第f个故障状态下的电压相角;BiifBijf分别是在第f个故障状态下,以$\frac{1}{{x_{ij}^f}}$为支路导纳建立起来的节点导纳矩阵的自导纳及互导纳,$B_{ii}^f = - \sum\limits_{j \in i, j \ne i} {\frac{1}{{x_{ij}^f}}, B_{ij}^f = - \frac{1}{{x_{ij}^f}}} {\rm{ }}$

2.3.2 节点负荷削减量约束

节点负荷削减量受限于节点的负荷功率,建立如下约束

$ 0 \le P_{{\rm{loss}}.i}^f \le P_{{\rm{load}}.i}^f。$ (17)
2.4 模型求解

由于式(2)~式(18)是典型的线性优化模型,文中采用成熟的商业软件CPLEX求解器计算,具体求解流程如图 1所示。

图 1 模型求解流程图 Fig. 1 The flow chart of model solution
3 算例分析 3.1 测试网络

为了验证文中方法的有效性,根据中国南方某沿海地区110 kV电网2018年的电源结构和网架结构等信息,仿真分析新能源渗透率与最小负荷削减量期望值的关系,评估系统的新能源消纳能力。2018年该区域有3座新能源电场,总容量为220.24 MW;1座火电厂,总容量为1 200 MW;4座蓄水式水电站,总容量为69.6 MW,年度负荷峰值为1 080.49 MW,网络N-1预想故障集中的故障数量为97个。

3.2 场景生成结果

在冬、夏和春、秋季节内各选取1个典型日来反映当地电网2018年的运行情况。在每个典型日内,以电网中的123个节点作为功率采样点,每隔5 min对各采样点进行1次功率采样获得源荷功率的原始数据。结合当地电网的实际运行情况,典型日内的峰平谷时段划分情况,如表 1所示。

表 1 典型日内的峰平谷时段划分情况 Table 1 The division of peak, flat and valley periods during the typical day

根据采样获得的原始数据,采用文中提出的场景生成方法,分别在3个典型日的峰平谷时段内各生成3个典型场景,计算各个典型场景发生的概率。以2018年春/秋季的典型日3月21日的峰时段为例,得到典型场景如图 2所示。依据2018年3月21日平时段、谷时段和冬季、夏季典型日各时段内的原始数据,同样可以生成对应的典型场景和获得各典型场景发生的概率,不再详细罗列。

图 2 2018年3月21日峰时段的典型场景生成结果 Fig. 2 The typical scenarios of peak periods on march 21st, 2018
3.3 评估结果 3.3.1 新能源渗透率和最小负荷削减量的关系

基于测试网络和典型场景数据,在正常运行状态下,实现系统内的新能源全额消纳,避免系统出现线路过载和切负荷等问题;在故障状态下,根据N-1安全校核的预想故障集,计算出各个典型场景下的最小负荷削减量,如表 2所示。

表 2 各典型场景下的最小负荷削减量 Table 2 The minimum load-shedding under each typical scenario  

进一步结合预想故障数量和各个典型场景概率,可获得最小负荷削减量的期望值。经计算可知,2018年该区域最小负荷削减量的期望值为0.288 1 MW。根据新能源渗透率的定义[19],可获得当前该区域的新能源渗透率为20.38%。调整新能源电场的装机容量,并利用文中方法对最小负荷削减量的期望值进行求解,获得不同新能源渗透率下的最小负荷削减量的期望值。需要注意的是,在改变新能源渗透率时仍需实现系统在正常运行状态下全额消纳新能源,并且不出现线路过载和切负荷等问题。根据不同新能源渗透率下的最小负荷削减量的期望值,两者之间的关系曲线,如图 3所示。

图 3 新能源渗透率和最小负荷削减量期望值的关系曲线 Fig. 3 The relation curve between renewable energy penetration and minimum load-shedding expected value

可以看出,当新能源渗透率在0~68.74%时,最小负荷削减量的期望值为固定值,与新能源是否接入、新能源渗透率是否提高均没有关系,其值始终为0.288 1 MW,因为在该渗透率区间内,电网的负荷需求主要依靠火电供给,某些线路故障后系统将失负荷,在负荷功率一定的前提下,系统失去的负荷量也是固定的。当新能源渗透率高于68.74%时,最小负荷削减量的期望值同新能源渗透率成正相关,负荷削减量将随着新能源消纳电量的增加而增大,原因是此时配电网中更多的负荷需求将依靠新能源电源供给,并且随着渗透率的增加,新能源电源输出功率增多,导致电源周围的支路更为薄弱,使得故障发生后线路越限的数量增加,造成切负荷量的增多,给系统的安全可靠运行带来更多的隐患。受限于新能源电场输电线路传输容量的限制,新能源渗透率的上限为84.58%,对应最小负荷削减量的期望值最大为0.759 2 MW。当新能源渗透率超过上限时,系统在正常运行状态下会出现线路过载的现象,无法安全运行。根据新能源渗透率上限,可以获得系统的新能源消纳能力,该区域年度新能源消纳电量上限为3.592×105 MW·h。

3.3.2 忽略电源的灵活性时新能源渗透率和最小负荷削减量的关系

为了进一步分析系统中灵活性电源的作用,在模型中不考虑蓄水式水电站的调节作用,求解不同新能源渗透率下的最小负荷削减量的期望值。图 4给出了系统忽略灵活性电源的调节作用时,新能源渗透率和最小负荷削减量期望值的关系曲线。

图 4 忽略电源灵活性时新能源渗透率和最小负荷削减量期望值的关系曲线 Fig. 4 The relation curve between renewable energy penetration and minimum load-shedding expected value ignoring generation flexibility

可以看出,如果不考虑系统内灵活性电源的调节作用,新能源渗透率高于68.74%时,最小负荷削减量的期望值将和新能源渗透率成正相关,但是新能源渗透率上限仅为78.24%。这是因为忽略系统中灵活性电源的调节作用,新能源渗透率为78.24%时,线路在某些场景下会过载,导致新能源无法全额消纳。根据新能源渗透率上限,获得该区域年度新能源消纳电量上限为3.312×105 MW·h,该值明显低于考虑电源灵活性时系统的新能源消纳能力。此外,新能源渗透率上限对应的负荷削减量期望值为0.377 8 MW,该值高于图 3中新能源渗透率78.24%时的负荷削减量期望值0.369 2 MW。这是因为预想故障发生时,系统缺少可灵活调节电源,难以安全可靠地满足负荷的供电需求,往往导致更多的负荷削减量。由此可见,利用系统内的灵活性电源可有效提高系统的新能源消纳能力,减小新能源消纳时的负荷削减量,提高系统的安全可靠性。

4 结束语

为了充分利用地区电网水电资源的灵活性以及有效评估地区电网的新能源消纳能力,文中提出了考虑灵活性的地区电网新能源消纳能力评估方法,该方法考虑了系统的正常运行状态和N-1预想故障状态。在正常运行状态下,利用区域内电源的灵活性来实现新能源的全额消纳,避免系统出现线路过载、切负荷等问题;在预想故障状态下,进行N-1安全校核,优化最小负荷削减量的期望值。结果表明,当新能源渗透率超过一定值时,提高新能源渗透率会增加最小负荷削减量的期望值,意味着为了消纳更多的新能源,系统在故障状态下将以更多的负荷削减量为代价。通过对比分析,利用系统内的电源灵活性能够有效地提升新能源消纳能力,减少最小负荷削减量,提高新能源消纳时系统的安全可靠性。

参考文献
[1]
中华人民共和国国家能源局. 2017年度我国新能源并网情况[R]. 北京: 中华人民共和国国家能源局, 2018.
National Energy Administration of PRC. China's Renewable Energy Grid-connected Condition in 2017[R]. Beijing: National Energy Administration of PRC, 2018. (in Chinese)
[2]
徐乾耀, 康重庆, 江长明, 等. 多时空尺度风电消纳体系初探[J]. 电力系统保护与控制, 2013, 41(1): 28-32.
Xu Q Y, Kang C Q, Jiang C M, et al. Preliminary analysis on wind power accommodation system from multiple temporal and spatial scale perspective[J]. Power System Protection and Control, 2013, 41(1): 28-32. (in Chinese)
[3]
邵剑. 基于风险理论的含大规模风电电力系统调度决策方法研究[D]. 武汉: 华中科技大学, 2015.
Shao J. Research on decision-making methods for large-scale wind power integrated system scheduling based on risk theory[D]. Wuhan: Huazhong University of Science and Technology, 2015. (in Chinese)
[4]
Ran L, Yuan T, Ma M, et al. Probabilistic method for evaluating wind power heating improving the accommodation capacity of renewable energy power system[C]//20182nd IEEE Conference on Energy Internet and Energy System Integration (EI2). Beijing, China, 2018.
[5]
贾文昭, 杨彬, 杨志刚, 等. 计及调峰调频能力的风电消纳能力评估[J]. 华北电力技术, 2014(5): 10-16.
Jia W Z, Yang B, Yang Z G, et al. Assessment of wind power accommodation and ability of peak regulation and frequency modulation[J]. North China Electric Power, 2014(5): 10-16. (in Chinese) DOI:10.3969/j.issn.1003-9171.2014.05.003
[6]
马建伟, 孙芊, 张景超, 等. 计及预测可信度时段差异性的含分布式可再生能源配电网调度策略[J]. 电力系统及其自动化学报, 2018, 30(5): 56-62.
Ma J W, Sun Q, Zhang J C, et al. Dispatching strategy for distribution network with distributed renewable energies considering forecasting credibility interval difference[J]. Proceedings of the CSU-EPSA, 2018, 30(5): 56-62. (in Chinese) DOI:10.3969/j.issn.1003-8930.2018.05.010
[7]
辛颂旭. 水电与风电联合运行分析[J]. 中国电力, 2013, 46(8): 85-89.
Xin S X. Combined operation analysis of hydropower and wind power[J]. Electric Power, 2013, 46(8): 85-89. (in Chinese) DOI:10.3969/j.issn.1004-9649.2013.08.017
[8]
曾雪婷, 刘天琪, 李茜, 等. 基于虚拟电源配置策略的风光水火多源互补短期优化调度[J]. 电网技术, 2016, 40(5): 1379-1386.
Zeng X T, Liu T Q, Li Q, et al. Short-term complementary optimal dispatch model of multi-source hybrid power system based on virtual power configuration strategy[J]. Power System Technology, 2016, 40(5): 1379-1386. (in Chinese)
[9]
李丹. 考虑时空相关性的源荷功率概率建模和概率预测方法[D]. 重庆: 重庆大学, 2017.
Li D. Research on the probabilistic modeling and prediction method of Source&Load power considering temporal-spatial correlation[D]. Chongqing: Chongqing University, 2017. (in Chinese)
[10]
彭春华, 于蓉, 孙惠娟. 基于K-均值聚类多场景时序特性分析的分布式电源多目标规划[J]. 电力自动化设备, 2015, 35(10): 58-65.
Peng C H, Yu R, Sun H J. Multi-objective DG planning based on K-means clustering and multi-scenario timing characteristics analysis[J]. Electric Power Automation Equipment, 2015, 35(10): 58-65. (in Chinese)
[11]
Tsekouras G J, Hatziargyriou N D, Dialynas E N. Two-Stage pattern recognition of load curves for classification of electricity customers[J]. IEEE Transactions on Power Systems, 2007, 22(3): 1120-1128. DOI:10.1109/TPWRS.2007.901287
[12]
Xu T, Chiang H, Liu G, et al. Hierarchical K-means method for clustering large-scale advanced metering infrastructure data[J]. IEEE Transactions on Power Delivery, 2017, 32(2): 609-616. DOI:10.1109/TPWRD.2015.2479941
[13]
胡锋, 朱承治, 汪志华. 基于改进K-means算法的电力负荷分类研究[J]. 电子测量技术, 2018, 41(12): 44-48.
Hu F, Zhu C Z, Wang Z H. Research on the power load classification based on improved K-means algorithm[J]. Electronic Measurement Technology, 2018, 41(12): 44-48. (in Chinese)
[14]
Stott B, Jardim J, Alsaç O. DC power flow revisited[J]. IEEE Transactions on Power Systems, 2009, 24(3): 1290-1300. DOI:10.1109/TPWRS.2009.2021235
[15]
Sauer P W. On the formulation of power distribution factors for linear load flow methods[J]. IEEE Transactions on Power Apparatus and Systems, 1981, PAS-100(2): 764-770. DOI:10.1109/TPAS.1981.316928
[16]
侯文庭, 韦化. 考虑核电可调度性的风-光-核-水-火多源协调短期优化调度[J]. 电工技术学报, 2018, 33(12): 2873-2882.
Hou W T, Wei H. A multi-source coordinated short-term dispatch model considering the dispatchability of nuclear power plants[J]. Transactions of China Electrotechnical Society, 2018, 33(12): 2873-2882. (in Chinese)
[17]
Chen Y, Liu F, Wei W, et al. Robust unit commitment for large-scale wind generation and run-off-river hydropower[J]. CSEE Journal of Power and Energy Systems, 2016, 2(4): 66-75. DOI:10.17775/CSEEJPES.2016.00051
[18]
Chen H, Xuan P, Wang Y, et al. Key technologies for integration of multitype renewable energy sources-research on multi-timeframe robust scheduling/dispatch[J]. IEEE Transactions on Smart Grid, 2016, 7(1): 471-480.
[19]
徐栎. 面向可靠性的含分布式电源配电系统相关问题研究[D]. 天津: 天津大学, 2015.
Xu L. Research on reliability issues of distribution power system with DG integration[D]. Tianjin: Tianjin University, 2015. (in Chinese)
图 1 模型求解流程图 Fig. 1 The flow chart of model solution
表 1 典型日内的峰平谷时段划分情况 Table 1 The division of peak, flat and valley periods during the typical day
图 2 2018年3月21日峰时段的典型场景生成结果 Fig. 2 The typical scenarios of peak periods on march 21st, 2018
表 2 各典型场景下的最小负荷削减量 Table 2 The minimum load-shedding under each typical scenario  
图 3 新能源渗透率和最小负荷削减量期望值的关系曲线 Fig. 3 The relation curve between renewable energy penetration and minimum load-shedding expected value
图 4 忽略电源灵活性时新能源渗透率和最小负荷削减量期望值的关系曲线 Fig. 4 The relation curve between renewable energy penetration and minimum load-shedding expected value ignoring generation flexibility
考虑灵活性的地区电网新能源消纳能力评估方法
谭涛亮 , 姜云鹏 , 李震 , 钱峰 , 任洲洋 , 刘俊磊 , 颜伟