邱革非, 杨昊天, 何 超, 刘铠铭, 何虹辉
(昆明理工大学 电力工程学院,云南 昆明 650504)
近年来，随着社会经济的飞速发展，环境和能源问题逐渐受到重视。以光伏为代表的绿色能源得到了大规模推进[1]。大量户用光伏发电系统接入配电网，导致配电网中普遍存在不对称负荷以及不对称线路参数[2]，光伏系统的非全相运行，加剧了配电网的三相不平衡特性。对户用光伏发电系统接入配电网后系统运行的研究已有成果。文献[3]考虑负荷与风机的不确定性，针对上级购电成本和网损最小建立了配电网最优潮流的鲁棒模型。文献[4]采用二次规划模型，以网损和弃光最小为目标函数，解决了光伏渗透率较高情况下的最优潮流问题。文献[5]在主动配电网最优潮流计算中采用了二阶锥松弛技术来对模型中的非凸约束进行线性化处理，将优化模型转变为标准的二阶锥规划(SOCP)问题，采用不同算法对求解结果进行误差分析，试验表明松弛误差满足了计算需求。但以上研究均未考虑配电网三相不平衡的特性，现有的文献大多在配电网三相平衡的运行条件下对一相进行潮流分析。随着配电网三相不平衡的情况加剧，已有研究显然不能够正确反映当下系统中各相复杂的情况。低压配电网的三相不平衡运行会增加系统损耗，重负荷所在相使母线电压降低，从而影响用户电能质量[6]。而面对不平衡相的不正常运行状况进行无功补偿和继电保护装置的安装均依赖于系统的潮流参数。可见，建立一个针对三相不平衡的低压配电网系统进行最优潮流快速求解的平台尤为重要。
求解方法上，大量智慧算法的提出使一些学者面对各类线性和非线性模型时，直接采用智慧算法求解。文献[7]针对电力系统最优潮流计算问题，提出了一种基于Q学习和纵横交叉搜索的粒子群算法，使用状态-组合动作链的方法解决了Q学习方法中维度灾难的问题，采用纵横交叉算子提高了粒子群算法寻优过程中的收敛速度。文献[8]采用樽海鞘群算法针对电压偏移、网损、发电成本、稳定度等多目标进行求解，研究了种群中领导者选取最优个体的比例对求解结果的影响。但是智慧算法在可行域内迭代求解的过程中很容易陷入局部最优解，同时在最优解附近的反复迭代也降低了求解效率。
综上所述,本文从构建求解电网单相潮流的节点导纳矩阵出发，以支路潮流模型为基础建立了三相不平衡系统各相最优潮流求解模型，考虑了包括分布式电源、储能、离散及连续无功补偿装置等控制单元对不平衡电网各相的补偿情况，采用SOCP方法将原非凸非线性模型松弛为典型的凸优化模型，提高求解速度。最后，基于Gurobi求解器得到该模型的全局最优解，仿真结果验证了该方法的正确性和可行性。
1.1 配电网三相四线制网络拓扑结构
我国的配电系统大多采用三相四线制拓扑结构[9]，节点l与节点m之间的配电线路如图1所示，节点间均采用中性点接地的星形接法，线路的首段中性线作为该模型唯一的参考节点。各相线路自身有自阻抗，线路之间的耦合关系用互阻抗表示，线路与用户设备相连形成闭合回路。
图1 配电网三相四线制线路模型
1.2 线路模型
根据上述配电网的网络拓扑结构，考虑到各相之间的耦合关系，2个节点间的关系可以用4×4的节点阻抗矩阵Zl,m来表示[9]：
(1)
式中：Zgg(g取a,b,c,n)为串联阻抗矩阵的对角元素，作为三相线路和中线的自阻抗；
Zgh(h取a,b,c,n)为串联阻抗矩阵中非对角元素(g≠h)，作为三相线路和中线的互阻抗。
对应地，该阻抗矩阵的逆矩阵(Y=Z-1)作为这段线路的导纳矩阵Yl,m。
含有m个节点的配电网节点导纳矩阵Y可以表示为
(2)
1.3 低压配电网负载模型
为了使潮流计算的结果能够正确反映系统中的真实情况，本文采用了考虑负荷静态特性的多项式ZIP模型(Z为恒阻抗模型，I为恒电流模型，P为恒功率模型)计及用户的电压特性或静态特性。可求得节点上s相的注入电流为[10]
(3)
1.4 可控协调装置模型
该模型中的储能装置(ESS)、以光伏为代表的分布式电源(DG)逆变器、分组投切电容器(CB)、静止无功补偿器(SVC)、有载调压开关(OLTC)作为可调有功、无功源。在保证配电网供需平衡的基础上，通过调节无功功率输入和输出以改变线路上电压的大小。白天光伏发电处于高峰期而用户用电处于低谷期，通过吸收无功功率降低电网过电压；
晚上光伏发电处于低谷期而用户用电处于高峰期时发出无功功率提升电网电压。
1.4.1 连续无功补偿装置
本文拟采用SVC作为连续无功补偿装，其功率极限如下：
(4)
1.4.2 离散无功补偿装置
本文采用CB组作为离散无功补偿装置,其运行条件如下：
(5)
(6)
(7)
(8)
1.4.3 OLTC
OLTC的运行条件如下，式(9)分别为高、低侧电压与变比的关系；
式(10)为变比与OLTC抽头的位置与变比之间的关系；
式(11)为OLTC抽头的最大位置限制；
式(12)为调度周期内OLTC抽头的调节限制:
Ul,t=nlm,tUj,t
(9)
nlm,t=nlm,0+Klm,tΔnl,m
(10)
-Klm,min≤Klm,t≤Klm,max
(11)
(12)
1.4.4 光伏逆变器
本文主要考虑户用光伏发电，集聚后的光伏系统逆变器容量需满足以下关系式：
(13)
1.4.5 ESS
储能系统的运行条件如下：
(14)
1.5 配电网的潮流算法公式
配电网各节点各相的电压向量方程：
Y×V(t)=Iinj(t)
(15)
式中：V(t)为t时刻时各相电压值形成的N阶向量；
Iinj(t)为t时刻各相向节点注入的电流形成的N阶向量。
为求得各节点中各相电压值，将式(15)变形为
V(t)=Y-1×Iinj(t)
(16)
2.1 目标函数
网络损耗是判断配电网经济指标的重要参考依据，以一天24 h为周期，配电网中的有功网损最小作为目标，建立三相四线制配电网的最优潮流模型[11]。
(17)
式中：Ploss为该配电系统24 h各支路有功损耗之和；
E为低压配电网支路集合；
T为全天各时段总数；
rlm为支路l-m的电阻；
Ilm,t为在t时段内支路l-m的电流。
2.2 运行约束条件
2.2.1 潮流约束
以经典Disflow潮流模型为例，选取结构为辐射状的低压配电网络中某一相支路在t时刻下的运行状态建立支路潮流模型如图2所示。
图2 配电网单相支路潮流模型
该支路潮流模型应满足的约束条件为[12]
(18)
(19)
(20)
(21)
式中：i、j为节点的编号；
Pj,t、Qj,t分别为t时刻节点j的有功注入功率和无功注入功率；
Pij,t、Qij,t分别为t时刻支路i-j的首端有功、无功功率；
Pjk,t、Qjk,t分别为t时刻支路j-k的首端有功、无功功率；
k为以j节点为父节点的所有节点集合；
Ui,t、Uj,t分别为节点i、节点j在t时刻的电压。
2.2.2 支路电流约束
(22)
2.2.3 电压约束
三相低压配电网每个节点中任一相的电压幅值均包含限值以保证电网的安全运行：
(23)
所建立的最优潮流优化模型是一个含有混合整数变量的非凸非线性模型，且无法保证解的收敛性，该形式下的模型很难求得最优解。为有效、快速地得到最优解，本文采用二阶锥方法对该模型进行松弛变换。
二阶锥的标准形式为[3]
(24)
式中：x∈Rn为n阶向量；
Ai∈Rm*n、bi∈Rm、ci∈Rn、di∈R均为已知常数。
对于支路(i,j)∈E有支路潮流约束如下：
Vi-Vj=zijIij,∀(i,j)∈E
(25)
(26)
将式(25)代入式(26)可以得到：
(27)
(28)
(29)
(30)
(31)
经过以上一系列的松弛变换，原本性质为非凸非线性、NP-hard的模型转化成了如下的经典SOCP模型：
(32)
为了提高计算速度、简化求解难度，二阶锥松弛对约束部分放松了条件，其松弛过程必然会产生误差，通过下式定义松弛误差[13]：
(33)
文献[14]对二阶锥松弛变换在数学模型上的有效性有更为深入研究，文中表明了在特定的条件下，通过二阶锥变换的模型依然具有精确性。本文在算例部分对松弛误差进行了求解分析，验证了二阶锥松弛对模型处理的有效性。
该模型已具有全局最优性，借助成熟的Gurobi算法包可进行求解。本文在Python平台下通过pyomo建模工具,基于二阶锥松弛，对上述大量户用光伏与负载接入导致三相不平衡的低压配电网容易发生单相电压越限的情况进行建模分析。调用Gurobi算法包，求得以有功损耗最小为目标、储能系统及无功补偿装置等为可调设备下全局的最优解。
4.1 参数设置
本文采用IEEE 33节点网络进行算例仿真分析，IEEE 33节点如图3所示，相关参数见文献[15]。
图3 IEEE 33节点
考虑光伏发电作为分布式电源，假设100户住宅用户分布在该模型内，户用光伏发电集中在用户侧，每户均安装光伏发电，参数见文献[16]。为体现算例中各相负载的不平衡特性，令算例中c相配置的用户数明显大于a、b两相的用户数。节点15和节点32装设储能设备配合DG实现削峰填谷；
节点6和节点16中的三相分别装设单相无功补偿设备CB，其每组容量为50 kvar，共10组；
在节点6、节点16、节点32三相各相均安装单相无功补偿设备SVC，补偿范围为-0.1～0.3 Mvar；
节点33上装设有载调压变压器，OLTC的分接头调整步长为0.01，且单日最大调整次数为5次；
节点电压运行范围为0.98～1.13 p.u.。利用本文模型对IEEE 33节点配电系统一天24 h的运行工况进行优化。
4.2 结果分析
调用Gurobi求解器对松弛后的模型进行求解，得到各时段有功损耗之和为691.25 kW，系统优化前有功损耗为1 036.88 kW，优化后的有功损耗降低到优化前的66.7%。该系统节点6上各相电压幅值以及整个配电系统各时段优化前后的损耗对比图如图4和图5所示。
图4 IEEE 33节点配电网控制条件下节点6上a、b、c三相电压
图5 优化前后有功损耗对比
由图5可以看出，系统各时段的有功损耗曲线与负荷需求走势相似，系统的有功损耗与负荷成正比关系。本文的配电系统模型可以通过控制协调以单相无功补偿装置为代表的可控单元，对运行中的三相不平衡系统进行优化，从而减少有功损耗，提高配网运行的经济性。
节点6上各相无功补偿装置的补偿情况如图6～图8所示，OLTC变化情况如图9所示，各单位可控装置在运行条件范围内，均跟随负荷变化而调整了补偿量。
图6 节点6各时段a相无功补偿装置出力情况
图7 节点6各时段b相无功补偿装置出力情况
图8 节点6各时段c相无功补偿装置出力情况
图9 各时段OLTC的变比情况
从图6～图8中可以看出，节点6上各相的无功补偿装置为了防止无功不足导致的母线电压降低进行无功功率补偿，其各相上的设备均能基于单相潮流参数对系统进行优化。由于c相上的负载比a、b相上大，其无功补偿量也相对较大。同时，补偿集中在15～20 h之间,这是由于傍晚系统中的户用光伏正处于发电功率的低谷期，而居民用电却处于高峰期。光伏发电量与负荷用量时序的不匹配导致电压容易越下限，无功补偿装置需要输出足够的无功功率，从而保证节点电压的稳定性。
ESS的充放电功率以及SOC在各时段变化量如表1所示。
表1 ESS在各时段功率输出与SOC
ESS的输出功率大于0为充电状态，小于0为放电状态。由表1中数据可知，ESS在白天负荷低谷期充电，在晚间的负荷高峰期放电，具有较好的削峰填谷的作用。
为了验证二阶锥松弛处理模型的准确性，按照式(33)的计算式，求得多时段各支路的松弛误差，如图10所示。从图中可看到松弛误差在10-6量级，满足该情况下的运行要求，证明该方法具有较高的可行性。
图10 各支路多时段下的松弛误差
本文针对低压配电网三相不平衡运行情况，同时考虑分布式电源和单相无功补偿装置等可控设备，以配电网有功损耗最小为目标，建立了多时段配电网最优潮流模型；
并采用二阶锥松弛方法将混合整数的非凸非线性模型转换为易于求解的线性模型；
通过IEEE 33节点配电系统案例分析，验证了该方法的可行性与有效性，为进一步研究三相不平衡配电网调度与规划问题提供了参考。