计及合环约束的分布式电源与配电网联合优化

2023-09-06唐海天冯银平

上海电力大学学报 2023年4期

唐海天, 冯银平

(1.上海电力大学, 上海 200090; 2.上海西门子医疗器械有限公司, 上海 201318)

双碳战略提出以来,大力发展可再生能源已成为一种趋势[1-2]。分布式电源(Distributed Generation,DG)作为可再生能源的重要组成部分,其渗透率正在逐年提高,但存在出力不稳定、消纳难等问题,对配电网的冲击日益增强[3-4]。配电网作为输电和用电的中间环节,对于供电可靠性及电能质量起到重要作用[5-6]。除了来自DG的冲击外,配电网中还存在负荷分布不均衡、网损高等问题,因此合理规划配电网,有利于提高DG的消纳,并优化线路的运行。

目前,对于配电网规划问题的解决方案主要有以下两种:一是对配电网络进行扩容,如新建变电站、开关站等[7],但这种方法成本很高,并且许多城市用地紧张,没有多余的空间进行扩容;二是进行配电网重构[8],通过改变联络开关的状态,对线路的拓扑结构进行优化重构。

在主动配电网环境下,拓扑结构因其可灵活调整的能力常被用于智能配电网的主动管理策略[9],目前已有学者提出采用配电网重构的方式优化配电网的运行[10-11],但在DG接入配置的问题上,少有学者考虑网络重构[12-15]。文献[12]建立了双层规划模型,迭代求解DG选址定容问题以及配电网最优拓扑问题;文献[13]考虑了系统负荷增长的问题,并将网络重构、DG选址定容作为应对该问题的手段;文献[14]计及负荷时序特性,建立了多场景下综合规划的模型。但以上模型均为基于传统人工智能算法求解的模型,无法直接调用商业求解器,也无法保证求解结果的全局最优性。文献[15]考虑动态重构,并在选址定容模型中引入了虚拟支路电压,但未计及运行过程中DG的时序特性,也没有对合环电流进行校验,难以保证运行的安全性以及结果的可行性。

基于以上研究,本文提出了一种计及DG时序特性的双层优化模型,通过二阶锥松弛技术转换为混合整数二阶锥规划(Mixed Integer Second-Order Cone Programming,MISOCP),并调用cplex求解器进行求解。此外,为避免产生不符合安全运行约束的不可行解,本文提出了改进的合环电流约束模型,当线路拓扑发生改变时对合环电流进行校验。

1 分布式电源与配电网联合优化模型

分布式电源与配电网联合优化模型是一个双层优化模型。上层模型为规划层,考虑线路拓扑重构得到分布式电源的最优接入容量及位置;下层为运行层,在运行中以开关开合为变量,进一步提升分布式电源的消纳。上下层模型的目标函数及约束条件相同,求解变量不同,下层模型以上层模型的结果作为输入,计及分布式电源与负荷的时序特性,联合求解最优的运行规划。

1.1 目标函数

目标函数主要基于线路网损、电压偏差和经济效益3个指标。其公式为

minW=αPloss+βΔUΣ+γCop

(1)

式中:W——综合优化目标;

α、β、γ——网损、电压偏差量、开关动作成本的权重;

Ploss——线路网损;

ΔUΣ——电压偏差量;

Cop——开关动作成本。

线路网损公式为

(2)

式中:T——一日中划分的时间段数,取24;

Nline——线路的支路数;

Rij——支路ij上的电阻;

Iij,t——t时刻支路ij上流过的电流。

电压偏差量公式为

(3)

式中:Nd——线路的节点数;

UN——节点处的额定电压;

Ui,t——t时刻节点i处的电压。

开关动作成本公式为

(4)

式中:μ——单次动作成本,取7元/次[16];

Jij——支路ij开关动作计数,若动作则取1,否则取0。

1.2 约束条件

配电网重构的过程中需要满足以下5项约束。其中,为确保可调用商业求解器进行求解,对于非线性项采用文献[17]中的二阶锥形式。

(1) 节点功率平衡约束

Sin=SLoad-SDG

(5)

式中:Sin、SLoad、SDG——节点的注入视在功率、负荷视在功率和分布式电源视在功率。

(2) 节点电压约束

(6)

式中:Umin、Umax——系统允许的最低电压和最高电压;

ui,t、uj,t——过程变量,数值等于t时刻节点i和节点j处电压的平方;

Pij,t、Qij,t——t时刻线路ij上传输的有功功率及无功功率;

Xij——线路ij上的电抗;

iij,t——过程变量,数值上等于t时刻支路ij上电流的平方。

(3) 传输功率容量约束

(7)

式中:rij,t、Tij,t——过程变量,展现潮流约束的二阶锥形式,无实际物理意义;

Gij、Bij——线路ij上的电导与电纳,若两节点之间不存在物理拓扑关系时取0;

Pmax、Qmax——线路允许通过的有功、无功功率的最大值;

Uj,t——t时刻节点j处的电压;

θij,t——t时刻线路ij两端的电压相位差。

(4) 辐射状约束配电网通常是闭环设计、开环运行,在重构的过程中不允许出现环网或孤岛,因此需要满足以下辐射状约束:

(8)

式中:ΦLine——系统各支路集合;

Fij——布尔变量,取1时表示支路ij联通;

Nroot——线路的根节点数;

u(j)——节点j的父节点的集合;

Φnode——系统中除电源节点外其他节点的集合,任意节点的父节点有且仅有1个。

(5) 开关动作约束

(9)

式中:NM——每日开关动作总次数的上限;

Nm——单个开关每日动作次数的上限。

1.3 分布式电源时序模型

分布式电源的时序特性与配电网运行密切相关,在上层优化模型确定好额定容量及接入位置后,考虑风速及光照强度,求得分布式电源的时序特性曲线,并将其叠加至线路负荷进行下层优化。

(1) 风力发电机组出力模型风速近似认为服从WeiBull分布[18],其有功出力与风速的关系式[19]如下

(10)

式中:Pw、PN——风力发电机组实际有功功率与额定有功功率;

vt——t时刻风速;

vi、vN、v0——切入风速、额定风速及切出风速。

(2) 光伏电源出力模型光照强度近似认为服从Beta分布[20],其有功出力与光照强度关系式[21]如下

Ppv=npvELIAsurη

(11)

式中:Ppv——光伏电源出力;

npv——光伏板数量;

ELI——光照强度;

Asur——光伏板面积;

η——光电转换效率。

1.4 双层优化模型的优化变量及求解方法

尽管上下层模型的目标函数及约束条件保持一致,但优化变量完全不同,上下层优化变量表达式分别如下

χc,l={xi,loc,Si|χconf,opt},i∈Φnode

(12)

(13)

式中:χc,l——上层模型优化变量;

xi,loc——分布式电源的接入位置;

Si——节点i处分布式电源的接入容量;

χconf,opt——上层模型定容时最优情况对应的拓扑结构;

χconf,t——t时刻下层模型的拓扑结构;

χconf——下层模型优化变量。

上述分布式电源与配电网联合优化模型的总体架构如图1所示。

2 合环电流约束模型

为保证供电可靠性,重构的过程中多采用合环的方式[22]。为避免产生不可行解,本文提出了改进的合环电流约束模型。

2.1 合环电流计算

合环电流的计算采用基于戴维南等效定理的叠加法[23],将网络分解成有源网络和无源网络两部分。其原理如图2所示。

图2 采用叠加法分解网络的原理示意

对于合环前的纯辐射网,可通过潮流计算得到合环前的电流,具体公式为

(14)

式中:Ij,t、Pj,t、Qj,t——t时刻合环前节点j处注入的电流、有功及无功功率。

对于合环后的等效无源环网,不考虑原有电源及负荷,仅保留线路等值阻抗及合环前两端母线的电势差,电流的计算公式为

(15)

Ijk,t=Ij,t+Ic,jk

(16)

式中:Ic,jk——线路环流,取节点j流向节点k为正方向;

Ujk,Δ——节点j和节点k之间的电压降;

Zjk,eq——戴维南等效定理求得的j、k点之间的等值阻抗;

Ijk,t——合环后节点j流向节点k的电流。

但式中Ujk,Δ、Zjk,eq均为向量形式的非线性项,为避免模型非凸,本文将对Ic,jk作出如下改进,取过程变量如下

(17)

式中:ic,jk——过程变量,数值上等于线路jk环流的平方;

Ujk,re、Ujk,im——过程变量,数值等于节点j和节点k之间电压实部之差和虚部之差;

Uj,re、Uk,re——节点j和节点k处电压的实部;

Uj,im、Uk,im——节点j和节点k处电压的虚部。

将式(15)转换为二阶锥形式,即

(18)

式中:Rjk,eq、Xjk,eq——戴维南等效定理求得的节点j和节点k点之间的等值电阻和电抗。

合环电流的最大瞬时值iM及最大有效值IM计算公式分别为

(19)

(20)

式中:Ta——衰减时间常数。

2.2 合环电流约束

(1) 稳态约束合环重构后线路上任意时刻的电流应小于线路允许通过的最大电流,即

Ijk,t≤Imax

(21)

式中:Imax——合环后线路上允许流过的最大电流。

(2) 暂态约束在合环操作中,应保证线路的冲击电流可以躲过1段速断保护,暂态电流的最大有效值躲过2段保护,稳态电流应躲过3段的整定值,可用公式表示为

(22)

式中:Iset1、Iset2、Iset3——线路1段、2段、3段的继电保护整定值。

3 算例分析

3.1 算例说明

本文以IEEE33节点模型为例,其拓扑结构如图3所示。在MATLAB环境下,调用cplex求解器对模型进行求解。线路与分布式电源的详细参数分别参考文献[24]和文献[25]。

图3 IEEE33节点配电网模型拓扑结构

为验证本文方法的有效性,设置4个算例对不同的优化方法进行对比,具体如表1所示。

表1 算例设置

此外,针对算例4,在不同场景和求解算法下进行进一步分析,研究联合优化的最优情况,并验证本文方法的优越性。

3.2 不同优化方法结果对比

不同优化方法的结果对比如表2所示。

表2 不同优化方法结果对比

与初始情况算例1相比,算例2和算例3的优化方法下网损分别降低了10.2%和15.2%,算例4的联合优化方法下网损降低了27.4%,优化效果更加明显。其优化前后的电压时空分布如图4所示。由图4可以看出,联合优化后的电压分布更加紧致,有效避免了分布式电源无序接入导致的电压越限等问题。因此可以得出结论,联合优化方法可以更好地优化配电网的运行。

配电网重构的过程中需要进行2～10次合解环操作,通过校验各次合解环电流的最大值可确保解的可行性以及过程的安全性。算例4的具体求解结果如表3所示。此处选取中午12点进行校验,仅保留每案例中合环稳态电流及冲击电流的最大值。