背景技术
随着全球能源结构向清洁、低碳化转型,我国已提出“2030碳达峰、2060碳中和”的宏伟目标。配电网作为连接用户与大电网的枢纽,将面临大规模分布式电源(distributedgenerations, DG)、电动汽车、储能、电采暖等新低碳元素的接入,由此引发的系统承载能力不足、负载不均衡、电压越限等问题给配电网规划运行带来严峻挑战。
从网架升级角度,柔性互联技术是解决上述问题的有效手段,其含义为采用电力电子柔性互联设备(flexible interconnected devices, FID)升级/构建配电网联络节点/通道,利用FID的动态潮流控制能力以及故障隔离能力,实现配电网柔性闭环运行,具体包括:①实时共享互联设备的容量,包括正常运行均衡负载以及故障负荷快速转供;②动态输出无功,平抑系统电压波动;③直流母线接入数据中心等直流负荷,减少变流环节,提高能量转换效率。
配电网柔性互联技术已经历了10余年的理论研究,并开展了一些示范工程建设,主要针对高压和中压配电网,例如英国的Network Equilibrium项目,北京怀柔的三端柔性环网示范工程等,而0.4kV的低压柔性互联配电网(low-volatge flexible distributionnetwork, LVFDN)研究相对较少。相对柔性中压互联,低压FID由于技术难度与成本相对较低、应用场景多样,将有潜力率先规模化应用,形成相邻台区广泛柔性互联的“蜂窝状配电网”。
文献《考虑规模化快充负荷的低压互联配电台区风险评估》对比了配电台区中快充负荷分散接入的传统方式和台区低压侧互联的新型接入方式,对台区柔性互联风险进行了评估;文献《交直流配电台区建设改造模式研究》以低压配电台区的整体负荷均衡为目标,提出基于直流线路的低压配电台区负荷转供方法;文献《规模化电动汽车负荷的柔性台区协同经济调度》提出了考虑规模化电动汽车快充负荷的柔性台区经济调度方法;文献《低压智能柔性互连交直流混合配电网设计》设计了海宁某台区基于柔性互连的低压交直流混供配电网方案。示范工程方面,英国2014年启动了FUN-LV工程,其低压FID规模达到24个,展示了较好的经济、社会和环境效益;宁波北仑开展了国内最早的柔性台区示范工程;山东济南建成了基于台区智能融合终端的低压柔直互联系统。
上述研究对LVFDN结构、风险评估、调度方法等进行了研究,然而对于LVFDN的供电能力研究尚属空白。最大供电能力(total supply capability, TSC)是配电网规划、评估以及安全分析的经典指标,针对传统10kV中压配网,TSC已建立了从模型算法到应用的完善体系,已对中压柔性互联配电网的供电能力进行了如文献《柔性配电网的最大供电能力模型》的研究。文献《Newton C, Lang P, Terry S. Field trial results of powerelectronics in low-voltage distribution networks》也提及了N-1故障下柔性设备对负荷的柔性转带策略,但只考虑了配变通过柔性互联装置互为N-1备用的情况。这些研究对于LVFDN的TSC研究具有借鉴意义,但是在FID负荷连续分配、N-1后的负荷二次转供等难点问题上未有涉及。LVFDN的TSC模型由于中低压负荷转供的双重约束将更加复杂,现有求解方法不再适用,需要开展新的研究。
发明内容
本发明的目的在于克服现有技术的不足,提出计及低压台区柔性互联的配电网最大供电能力计算方法,通过分析台区柔性互联的配电网中-低压协同运行方式,建立考虑台区柔性互联、负荷多级转供等约束的配电网TSC模型,针对所提出模型特征进行求解得到配电网最大供电能力的计算结果。
本发明解决其技术问题是采取以下技术方案实现的:
计及低压台区柔性互联的配电网最大供电能力计算方法,包括以下步骤:
步骤1、根据台区柔性互联的配电网的结构及其特点,建立多种约束条件的配电网TSC模型;
步骤2、获取台区柔性互联的配电网参数;
步骤3、根据步骤1中构建的配电网TSC模型和步骤2中台区柔性互联的配电网参数进行解算,得到配电网最大供电能力的数值。
而且,所述步骤1中约束条件包括:配电网特性约束、运行约束和N-1安全约束。
而且,所述配电网特性约束包括:
配电网特性约束1、配电网节点功率从电网流出方向为正,注入为负;
配电网特性约束2、由于城市电网线路短并且网损小,简化馈线出口潮流包括网损,并采用直流潮流计算;
配电网特性约束3、故障集考虑主变故障退运、馈线出口线路故障退运和配变故障退运,不考虑10kV支线故障以及0.4kV线路故障;
配电网特性约束4、城市地区分布式电源渗透率低,节点净功率呈中压流向低压供电特性,若不考虑DG,则转化为纯负荷情况的供电能力计算,否则将将DG简化功率为负的负荷进行计算。
而且,所述运行约束为:配电网台区个数为n,台区的配变高压进线净功率等于配变所供用户总视在功率以及与台区互联的FID端口注入功率之和,考虑正常运行时配变所供用户总视在功率以及与台区互联的FID端口注入功率之和需要小于配变的额定容量,以及不允许配变潮流倒送的约束,则:
其中,
和
为序数;
为台区的配变所供用户总视在功率;
为与台区
互联的FID端口注入功率之和;
为配变
的额定容量;
;
为重载系数;
FID端口转移功率需满足端口容量约束为:
FID的各个端口功率之和满足基尔霍夫方程,并且忽略装置自身损耗:
其中,
为柔性装置端口容量;
为FID所有端口编号的集合;
在直流潮流的条件下,配电网潮流简化为功率平衡方程,线路
流过的功率等于其下游全部台区配变高压侧进线净功率之和,线路容量约束为:
其中,
为线路
流过的功率;
为线路
下游台区所有编号的集合;
为线路
的额定容量,
为台区
的配变所供用户总视在功率;
为与台区
互联的FID端口注入功率之和;
其中
为主变
流过的功率;
为主变
下游台区所有编号的集合;
为主变
的额定容量。
而且,所述N-1安全约束为:若故障集为
,单一故障
发生后,新拓扑条件下线路
、主变
下游台区所有编号的集合分别为
和
,线路容量约束为:
主变容量约束为:
故障发生时,负荷
通过互联的FID的端口
转出到其他端口,并按比例分配到其他互联配变:
其中,
为线路
流过的功率,
为配变
的所供用户总视在功率,
为与台区
互联的FID端口注入功率之和;
为线路
的额定容量,
为主变
流过的功率,
为线路
下游台区所有编号的集合,
为主变
下游台区所有编号的集合,
为主变
的额定容量,
为台区的配变所供用户总视在功率,
为台区
故障后经过FID的端口
向台区
转移的负荷比例系数。
而且,所述步骤1中构建的配电网TSC模型为:所有用户负荷
的和:
其中:
为TSC的数值,
为台区的配变所供用户总视在功率,建立的配电网TSC模型为非线性非凸规划模型,在达到TSC的前提下以台区用户负荷均衡为目标进行二次优化,目标函数为:
其中,
为台区负荷的方差,
为台区
配变负载率,
为配变
容量,
为
个台区配变负载率的平均值。
而且,所述步骤2中台区柔性互联的配电网参数包括:配电网结构参数、主变集合、馈线集合和配变集合。
而且,所述步骤3包括以下步骤:
步骤3.1、根据配电网结构参数,计算用户视在功率
和FID端口注入功率
;
步骤3.2、生成故障集
,其中
是单一故障,为故障集
中的故障元素;
步骤3.3、令
,由均衡转带原则,故障后配电网重构形成新拓扑;
步骤3.4、对主变集合、馈线集合、配变集合的元素依次进行N-1安全校验,同时取
;
步骤3.5、在新拓扑下,计及运行约束和N-1安全约束,并判断此时的
i是否达到最大,若达到最大,则进行步骤3.6,否则使
,并重复步骤3.5;
步骤3.6、判断
是否达到最大值,若达到最大值则由目标函数,联立约束条件,求解模型,得到配电网最大供电能力,否则使
,并返回步骤3.3。
本发明的优点和积极效果是:
本发明通过分析台区柔性互联的配电网中-低压协同运行方式,建立考虑台区柔性互联、负荷多级转供等约束的配电网TSC模型,并针对所提出模型特征进行求解得到配电网最大供电能力的计算结果。本发明考虑了LVFDN中-低压协同的灵活运行方式,能够更有效地利用负荷转供能力,进一步提升可靠性;同时建立了考虑台区柔性互联、负荷多级转供的LVFDN的TSC模型;并且针对所提出模型的非线性非凸规划特征,提出了基于分支定界算法的TSC模型求解方法,能够精确求解LVFDN的TSC数值以及变化。
具体实施方式
以下结合附图对本发明做进一步详述。
计及低压台区柔性互联的配电网最大供电能力计算方法,包括以下步骤:
步骤1、根据台区柔性互联的配电网的结构及其特点,建立多种约束条件的配电网TSC模型。
一个具备一定规模低压台区柔性互联的配电网(LVFDN)如图1所示。
正常运行时,基于FID的功率连续调节功能,互联台区能够通过FID实现负荷均衡,缓解配变重载、过载的风险;此外,基于FID的独立输出无功输出功能,还可以为台区提供电压支撑,缓解分布式光伏等接入造成的电能质量问题。
当某个低压台区的配变故障退出运行后,首先计算互联的其他台区配变的剩余容量,确定必须切除的负荷以及可以转供的负荷,然后这些可转供的负荷将通过FID快速转带至互联的其他台区配变。
当某段10kV中压馈线发生故障并隔离后,为恢复非故障区供电,配电网将同时采取两种手段(中-低压协同):一是通过中压线路开关操作转移非故障区负荷,尽可能的恢复供电;二是通过低压FID转移台区的部分负荷(互联台区中至少有一个处于非故障区)。两种手段协同运用,以损失负荷最小为目标。
当110kV或35kV变压器发生故障退出运行时,可以等效为若干条馈线失电,处理方式与馈线故障类似。
如图2所示为一个局部低压台区柔性互联配网的运行方式,在正常运行时,开关(包含负荷开关和断路器)K1、K2、K4、K6、K7闭合,K3、K5断开。此时,配变D1和D2、D3和D4分别通过FID连接,以负载率均衡为目标运行。
当配变D1故障时,D1的负荷总开关断开,D1所供的L1部分负荷(可转供部分)经FID转供至D2;
当馈线F2出口故障后,K4断开,K3闭合,此时,D3由馈线F1供电,同时考虑到馈线F1容量约束,可将L3部分负荷转供至D4(二次转供)。图2中箭头表明了故障后负荷(仅讨论有功功率)的转供路径,即负荷的电源由箭头始端的电源A变为箭头末端的电源B,蓝色箭头表示负荷一次转供路径,橙色箭头表示负荷二次转供路径。
约束条件包括:配电网特性约束、运行约束和N-1安全约束。
城市地区土地空间资源紧张,负荷面密度大,台区供电能力不足的矛盾更加突出,最有可能率先应用低压柔性互联技术。本发明针对城市场景特征,结合工程实际得到以下研究的基本假设,这些假设在现有供电能力研究中也经常采用。因此,配电网特性约束包括:
配电网特性约束1、配电网节点功率从电网流出方向为正,注入为负;
配电网特性约束2、由于城市电网线路短并且网损相对小,简化馈线出口潮流包括的网损,并采用直流潮流计算;
配电网特性约束3、故障集考虑主变故障退运、馈线出口线路故障退运和配变故障退运,不考虑10kV支线故障以及0.4kV线路故障。
配电网特性约束4、城市地区分布式电源渗透率低,节点净功率呈中压流向低压供电特性,若不考虑DG,则可以转化为纯负荷情况的供电能力计算,否则将将DG简化功率为负的负荷进行计算。
运行约束为:设配电网台区个数为n,台区的配变高压进线净功率等于配变所供用户总视在功率以及与台区互联的FID端口注入功率之和。考虑正常运行时配变所供用户总视在功率以及与台区互联的FID端口注入功率之和需要小于配变的额定容量,以及不允许配变潮流倒送的约束,有:
式中:
和
为序数;
为台区的配变所供用户总视在功率(不包含FID的注入功率);
为与台区
互联的FID端口注入功率之和,取正表示功率由配变流向FID(FID端口等效为负荷),取负则表示功率由FID注入配变(FID等效为电源);
为配变
的额定容量;
;
为重载系数,正常运行中一般取0.7~0.8,在N-1场景下可以近似认为等于1(允许短时重载运行)。
FID端口转移功率需满足端口容量约束:
式中:为柔性装置FID端口容量。
FID的各个端口功率之和满足基尔霍夫方程,并忽略装置自身损耗:
在直流潮流的条件下,配电网潮流简化为功率平衡方程,线路
流过的功率等于其下游全部台区配变高压侧进线净功率之和,线路容量约束为:
式中:
为线路
流过的功率;
为线路
下游台区所有编号的集合;
为线路
的额定容量,
为台区
的配变所供用户总视在功率;
为与台区
互联的FID端口注入功率之和。
类似的,主变i的容量约束为:
式中:
为主变
流过的功率;
为主变
下游台区所有编号的集合;
为主变
的额定容量。
N-1安全约束为:TSC定义配电网满足N-1安全准则的最大负荷供应能力。因此除了正常运行下的安全约束外,还必须考虑N-1安全约束。在LVFDN中,N-1安全约束指配电网单一元件在系统故障后退出运行,系统通过网络重构、FID功率调节等方式,至少能找到一种新的运行方式使得非故障区负荷仍能保持供电,同时满足系统安全性约束(1)—(5)。本发明考虑的N-1故障集包括单个台区配变故障、馈线出口故障和主变故障。
无论是单个台区配变故障、馈线出口故障或主变故障,系统都将发生两个变化:一是网络重构造成的拓扑变化,二是FID的端口功率变化,二者的目的都是尽可能多的保证不失负荷。在新的拓扑和FID功率分配下,系统仍要满足运行约束(1)—(5)。
记故障集为
,单一故障
发生后,新拓扑条件下线路
、主变
下游台区所有编号的集合分别为
和
,则运行约束(4)、(5)可以表示为:
现有TSC研究均未涉及台区N-1故障,主要原因是台区下低压配网线路均为辐射结构,若台区配变故障通常只能停电等待恢复,从结构上不符合N-1安全性。在LVFDN中,由于FID的存在,台区配变故障后负荷可以转移到其他配变。从另一个视角,FID将作为台区i负荷的新电源。设台区i配变发生故障,其负荷将尽可能通过互联的FID的端口i转出到其他端口,本质上是按一定比例分配到其他互联配变,这得益于FID的多端口潮流灵活分配能力,因此:
式中:
为线路
流过的功率,
为配变
的所供用户总视在功率,
为与台区
互联的FID端口注入功率之和;
为线路
的额定容量,
为主变
流过的功率,
为线路
下游台区所有编号的集合,
为主变
下游台区所有编号的集合,
为主变
的额定容量,
为台区的配变所供用户总视在功率,
为台区
故障后经过FID的端口
向台区
转移的负荷比例系数。
因此,依据TSC的定义,TSC模型的目标函数同传统供电能力模型一致,为所有用户负荷
的和。
传统配网的TSC模型中,只有用户负荷
S D,i 为变量,而在LVFDN的TSC模型中,还有FID的端口功率
为变量。TSC结果不仅对应了一组用户负荷分布,还对应了至少一组可行的FID的端口功率分布。由于在N-1后,模型的会重新优化一次FID的端口功率分布,这本质上是对柔性互联的台区负荷进行了二次转供——一次转供是中压刚性的开关重构实现,二次转供是调整FID端口功率。传统刚性配网受到开关操作灵活性和寿命的限制,一般不进行二次转供,而LVFDN可以解决这一问题,通过二次转供提升网络灵活性和可靠性。得到考虑低压台区柔性互联的TSC模型为:
由于(8)中
和
均为优化变量,同时N-1后不同拓扑、FID功率分配的方案以取“并集”的方式判定系统安全性,因此该模型为非线性非凸规划模型。
由于模型存在多解,一些解对应的负荷分布很不均衡,与实际差距较大。根据负荷均衡思想,在达到TSC的前提下以台区用户负荷均衡为目标进行二次优化,目标函数为:
其中,
为台区负荷的方差;
为
个台区配变负载率的平均值。
步骤2、获取台区柔性互联的配电网参数。
配电网参数包括:配电网结构参数、主变集合、馈线集合和配变集合。
步骤3、根据步骤1中构建的配电网TSC模型和步骤2中台区柔性互联的配电网参数进行解算,得到配电网最大供电能力的数据。
针对传统TSC模型的线性规划求解方法不再适用。如图3所示,本发明以分支定界算法为核心,并对子问题进行线性规划松弛和凸包络近似。这一方法能够快速遍历所有子问题,快速收敛至全局最优解,具有较好的鲁棒性。
步骤3.1、根据配电网结构参数,计算用户视在功率
和FID端口注入功率
;
步骤3.2、生成故障集
,其中
是单一故障,为故障集
中的故障元素;
步骤3.3、令
,由均衡转带原则,故障后配电网重构形成新拓扑;
步骤3.4、对主变集合、馈线集合、配变集合的元素依次进行N-1安全校验,同时取
;
步骤3.5、在新拓扑下,计及运行约束和N-1安全约束,并判断此时的
i是否达到最大,若达到最大,则进行步骤3.6,否则使
,并重复步骤3.5;
步骤3.6、判断
是否达到最大值,若达到最大值则由目标函数,联立约束条件,求解模型,得到配电网最大供电能力,否则使
,并返回步骤3.3。
相比于已有TSC模型和求解算法,本发明TSC模型和算法首次考虑了台区配变级别的负荷转带,并在此基础上实现负荷的二次转供。通过FID装置,对于多个备用的电源点(台区配变),负荷可以任意比例转供。
参考英国低压柔性城市网络工程(FUN-LV)的实际配电网,如附图3所示构造算例,以验证本发明所提出方法。统一取FID设备容量
为0.3MVA。将附图3简化为用点边图表示,拓扑结构如附图4所示。本发明只考虑台区、馈线出口处和主变的N-1故障,此时同线路上的低压柔性互联不影响
N-1负荷转带,故简化省去同线路上的低压互联线路。算例包括4台主变压器,11条中压馈线,5个FID(其中4个双端口FID、1个3端口FID、共11个端口)。主变T1容量为15MVA,T2、T3和T4的容量均为20MVA,馈线容量均为8MVA。配电变压器的容量均为0.6MVA。
根据本发明建立算例的TSC模型,并进行求解。在MATLAB平台中,采用Yalmip的全局优化求解器BMIBNB进行求解,在i5-8300H-8G计算机中,平均计算一组算例的TSC均衡解的时间为2.14秒。
计算得到TSC值(均衡解)为53.4MVA,对应的11条馈线负荷、以及11个参与柔性互联的台区负荷见表1。
表1 TSC水平下算例电网的各馈线、柔性台区负荷
采用经典的TSC准确性校验方法,对表1的TSC负荷分布进行N-1安全校验,结果表明,在当前TSC负荷下,配电网刚好满足N-1安全,即再以任何方式增加任意大小负荷,都至少存在一个N-1故障使得算例配网出现不安全的状态,这说明本发明所提出的LVFDN的TSC模型算法是有效的。为进一步验证,表2给出了一个略高于TSC水平的负荷分布。
表2某高于TSC水平负荷下算例电网的各馈线、柔性台区负荷
表2此时负荷之和为53.5MVA.。经N-1校验后,发现F1馈线在F3馈线出口故障和F4馈线出口故障场景下出现了容量越限情况,且越限容量即为增加负荷0.1MVA,见表3。
表3 某高于TSC水平负荷的N-1校验结果
因此,通过本发明计算出的LVFDN的TSC是准确的。
需要强调的是,本发明所述的实施例是说明性的,而不是限定性的,因此本发明包括并不限于具体实施方式中所述的实施例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,同样属于本发明保护的范围。