一种基于惩罚凹凸优化技术的电力系统最优潮流控制方法
【专利摘要】本发明公开了一种基于惩罚凹凸优化技术的电力系统最优潮流控制方法，包括以下步骤：首先获得电力系统网络参数，确定最优潮流控制问题的数学模型；接着利用惩罚思想和泰勒展开，得到最优潮流问题当前的近似凸问题；然后迭代地求解最优潮流控制近似凸问题得到各母线注入功率值；最后根据计算得到的母线注入功率值完成电力系统最优潮流控制。本发明利用惩罚凹凸优化技术设计电力系统最优潮流，能够在保证潮流方程可行性的前提下使得系统性能损耗最小。
【专利说明】
一种基于惩罚凹凸优化技术的电力系统最优潮流控制方法
技术领域
[0001] 本发明涉及电力系统技术领域，具体涉及一种基于惩罚凹凸优化思想的电力系统 最优潮流控制方法。
【背景技术】
[0002] 电力系统优化问题，包括规划、调度、运行于控制，其目标是系统安全性与经济性 的平衡和折中。作为其中最重要的问题之一，最优潮流(〇Ptima 1 Power FIow，0PF)控制是 指电力系统的结构参数和负荷情况都已给定时，调节可利用的控制变量(如发电机输出功 率、可调变压器抽头等)来找到能满足所有运行约束条件的，并使系统的某一性能指标(如 发电成本或网络损耗)达到最小值时的潮流控制。近年来，随着智能电网、分布式发电技术、 分布式电能存储技术的迅猛发展，在满足电力系统安全性的前提下，尽可能地提高经济性， 合理利用资源配置和现有设备以减少能源消耗的最优潮流控制这一经典问题又成为了研 究热点。
[0003] 从20世纪60年代以来，最优潮流作为电力系统运行和分析的强有力工具，一直倍 受关注。经过近50年的发展，众多最优化方法被相继引入该领域，如：线性规划、二次规划、 非线性规划以及牛顿法和解耦法等。但最优潮流是一个典型的非线性优化问题，且由于约 束的复杂性使得其计算复杂，难度较大。当前，文献[M.Farivar and S.H.Low，"Branch flow model: Relaxations and convexif i cat ion (parts I，II)，''IEEE Trans . Power Syst ·，vo 1 · 28，no · 3，pp · 2554-2572，2013 ]中提出了利用凸松弛方法--S0CP松弛求解最 优潮流问题，然而其只在一定条件下证明了该松弛是紧的。对于非凸的最优潮流问题，凸松 弛方法甚至都无法保证得到问题的可行解。因此，本发明提出基于惩罚凹凸优化的最优潮 流控制方法。
【发明内容】
[0004] 本发明的目的在于针对现有技术的不足，提供一种基于凹凸优化思想的电力系统 最优潮流控制方法，本发明方法考虑了分布式发电中逆变器的控制问题，可以保证在迭代 过程中电力系统性能损耗一直单调递减，能够实现在达到电力系统各约束要求的同时降低 电力系统性能损耗的目的，完成电力系统最优潮流控制。具体包括以下步骤：
[0005] 步骤1:获得电力系统网络参数:母线集合N和去除根母线后的集合Ν+ = Ν\{0};电 网支路集合ξ;母线之间支路的阻抗2^，'外/，./)￡。各母线节点注入功率81的约束集合5 1;母 线电压幅度值的平方Vl的下限ρ和上限丐；
[0006] 步骤2 :初始化迭代次数k = 0，最大迭代次数Kmax，收敛精度tol ;设定初始点 确定惩罚系数hh谦示从母线i到母线j的电流幅度值的平方;示母线i 至Ijj间的连线发送端的潮流；
[0007] 步骤3:利用泰勒展开，得到最优潮流问题的近似凸问题：
[0011] 其中|a| j和Re(a)分别表示复数a的幅度、共辄和实部;表示母线i处注 入功率所引起的系统性能损耗；
[0012] 步骤4:求解问题(P1)得到当前的计算结果判断是否达到收敛 精度：
或者是否达到最大迭代次数:k>Kmax;若是，输出母线的注入功率，计 算系统性能损耗，执行步骤5;否则令迭代次数k = k+l，重复步骤3和4。
[0013] 步骤5:根据计算得到的母线注入功率完成最优潮流控制。
[0014] 进一步的，所述步骤4中求解问题(P1)得到计算结果(^+1，矣+1,1^+1，^ +1；)的解决方 法为内点算法。
[0015] 本发明有益效果:本发明方法首先构造最优潮流对应规划问题，接着利用惩罚思 想和泰勒展开，得到最优潮流问题当前的近似凸问题;然后迭代地求解最优潮流控制近似 凸问题得到各母线注入功率值;最后根据计算得到的母线注入功率值完成电力系统最优潮 流控制。本发明利用惩罚凹凸优化技术设计电力系统最优潮流，能够在保证潮流方程可行 性的前提下使得系统性能损耗最小。
【附图说明】
[0016] 图1是本发明实施例采用该方法的系统模型图。
[0017] 图2是本发明实施例采用该方法的具体流程图。
[0018] 图3是本发明实施例的目标值与迭代次数的关系图。
[0019] 图4是本发明实施例系统约束可行性指标与迭代次数关系图。
【具体实施方式】
[0020] 为了使本发明的目的和效果更加清楚，下面结合附图对本发明方法的具体实施方 式进行详细说明。
[0021] 如图1所示，考虑辐射状的配电网，其由母线和连接母线的连线组成。该网络中的 根节点为变电站母线(为方便描述，下面称为根母线），其与输电网络相连。根母线使用固定 的电压，同时将从传输网络中的接收到的电力分配到其他母线。本发明定义该根母线为母 线0,其他母线为1，. . .，n;另外，令N: = . . .，n}表示电网中所有母线，定义N+: =N\{0}; (i，j)表示母线i和母线j相连，方向为i - j，且母线j在母线i与母线〇的唯一路径上。令ξ表 示网络中所有支路的集合，对任意（i，j)e|，表示有向支路i-j。
[0022] 对于任意母线ieN，令^表示母线i处的电压幅度值的平方。如上所述，变电站母 线的电压为固定值νο。定义Si = pi+iqi表示在母线i处的注入功率，其中pi、qi分别表示注入 的有功功率和无功功率。另外，定义P,为母线i到母线〇之间的唯一路径，对于辐射型网络，P, 是唯一的。对于任意连线（i，j) e ξ，令1 i」表示从母线i到母线j的电流幅度值的平方，z i j = rij+ixij表示母线i，j之间连线的阻抗;令Sij = Pij+iQij表示母线i到j间的连线发送端的潮 流(或称功率流），其中Pi」和Qi」分别表示有功功率流和无功功率流。另外，对于复数a e C，用 :反表不a的共辄。
[0023] 给定网络拓扑(Ν，ξ)、阻抗z以及变电站母线电压vo时，那么其他电网参数(s，S，v， l，so)可以通过福射网络的支流模型(branch flow model)表示如下：
[0028] 公式(la)和(lb)是功率平衡方程，公式(lc)和(Id)是欧姆公式的恒等变换。
[0029] 本发明考虑以下几种配电网可控设备:分布式发电机、逆变器、可控负载，比如电 动车辆、智能家电、并联电容器。在实际应用中，电网通过控制并联电容器和逆变器注入的 无功功率来调节电压。在设定注入功率s后，通过公式(1)可以确定其它电参数(S，v，l，so)。
[0030] 根据可控设备的不同类型，电网中母线ieN+的注入功率Sl具有不同的约束集合 Si，即：
[0031] SieSi，ieN+ (2)
[0032] 根据设备类型定义集合31为：
[0033] ①若Si代表一个额定容量为虿的并联电容器，那么萬=卜eC；| Kc⑷=(),M.vi = 0("彳丨, 若81代表一个最大发电量为Λ的太阳能电板，其通过一个容量为5；的逆变器与电网连接，那 么在=e !:C 丨 0 瓦，卜| ;
[0034] ②若Si代表一个功率因子为η、有功功率消耗在区间[-]?，-? ]连续变化的可调负 载，
[0035] 注意，Si可以表示多个上述设备总的注入功率。
[0036] 另外，需要将母线i的电压幅度值的平方Vl控制在预先设定的电压下限值^和电压 上限值巧之间，即需满足
[0037] <vri e.N+
[0038] 在功率流约束、电压约束、注入功率约束的条件下，最优潮流问题可描述如下：
[0047] 其中目标函数中-駿表示母线i注入功率所导致的系统性能损耗。若对于任 意ieN，有fi(x) = x:
即表示在电网中的总功率损耗。
[0048]
的非凸约束，上述最优潮流问题为非凸优化问题， 很难求角军。文南犬[M · Fari var and S · Η · Low，"Branch flow model ： Relaxations and convexification(parts I，II)，''IEEE Trans.Power Syst·，vol·28，no·3，pp·2554-2572， 2013]中提出了利用S0CP松弛方法求解最优潮流问题，然而只在一定条件下证明了该松弛 是紧的。对于一股的最优潮流问题，凸松弛方法甚至都无法保证得到问题的可行解。因此， 本发明提出基于惩罚凹凸优化的最优潮流控制方法。
[0049] 以非凸约束(3d)为例，首先将(3d)转化为两个不等式约束：
[0051]其中前者为凸约束，后者为非凸约束。对于后者，通过引入惩罚项，将其移至目标 函数中，得到惩罚问题：
[0053] 可以证明，当惩罚参数β大于某个门限值时。问题(4)可以通过凹凸优化进行求解。 具体地，给定Sij和Vi的当前将问题(4)中的惩罚项通过泰勒展开进行线 性逼近，即：
[0055] 可以得到如下凸问题，
[0057] 迭代地求解上述凸问题直至算法收敛，可以得到最优潮流控制结果。
[0058] 图2给出了上述基于凹凸优化技术的电力系统最优潮流控制方法的流程图。具体 地，可以描述如下：
[0059] -种基于惩罚凹凸优化技术的电力系统最优潮流控制方法，该方法包括如下步 骤：
[0060] 步骤1:获得电力系统网络参数:母线集合N和去除根母线后的集合Ν+ = Ν\{0};电 网支路集合ξ;母线之间支路的阻抗各母线节点注入功率Sl的约束集合S1;母 线电压幅度值的平方 Vl的下限Xi和上限丐;
[0061] 步骤2 :初始化迭代次数k = 0，最大迭代次数Kmax，收敛精度tol;设定初始点 ^\5<>?,/; )};确定惩罚系数比11请示从母线理」母线」的电流幅度值的平方;51谦示母线1 到j间的连线发送端的潮流；
[0062] 步骤3:利用泰勒展开，得到最优潮流问题的近似凸问题：
[0066] 其中|a|、5和Re(a)分别表示复数a的幅度、共辄和实部;/:徽4潑表示母线i处 注入功率所引起的系统性能损耗；
[0067] 步骤4:求解问题(P1)得到当前的计算结果[^+1，4+1,4+1,/^}，判断是否达到收敛 精度：
或者是否达到最大迭代次数:k>Kmax;若是，输出母线的注入功率，计 算系统性能损耗，执行步骤5;否则令迭代次数k = k+l，重复步骤3和4。
[0068] 步骤5:根据计算得到的母线注入功率完成最优潮流控制。
[0069] 进一步地，所述步骤4中求解问题(P1)得到计算结果的解决方 法为内点算法。
[0070] 下面通过具体实例对本发明的技术方案进行进一步阐述。实验中，采用SCE-47和 SCE-56电网系统进行算法验证。具体地，使用下面的实验参数：
[0071] ①设定功率损耗最小化为目标，且变电站母线的电压Vo为一个单位的基准电压 值；
[0072]②对于注入功率约束边界的设定，在电网中任意母线i e N+处可能存在多个设备， 如并联电容、可调负载、太阳能电板等；假设电网中总共存在Di个设备并将其编号为1， 2,. . .，Di;对于d=l，2,. . .，Di，sid表示设备d的注入功率。
[0073]若设备d是一个负载，且已知有功功率消耗p和无功功率消耗q，那么此时sld = -p-i · q;若已知负载d的视在功率峰值Speak，那么Sid = -speak exp(j9)。其中，0=arccos(〇.9)， 此时，负载的注入功率Sld也即一个常数；
[0074] 若设备d是一个容量为f的电容，那么有= ;
[0075] 若设备d是一个容量为歹的光伏电板，那么1??(心.)2〇,|~|<?;
[0076] 根据以上设定，此时母线i总的注入功率
[0077] 其他参数设定如下：根据电网实际情况确定各基准容量值、初始化母线总数N，电 阻r^、电抗Xlj以及在母线处相应设备类型的容量或有功功率消耗值，同时设定变电站节点 的功率为一个单位的基准功率值，令迭代总次数K max = 20，Vmax等于1 . 1单元电压基值，Vmin等 于0.9单元电压基值，收敛精度tol = 0.001，惩罚参数β = 0.001，初始化母线电压上限 5 =ν^、电压下限}^#^另外，本实施例中
为约束可行性指标，其值 接近0时说明已到可行。
[0078]图3、4是通过Matlab对所设计方法的仿真验证结果图。
[0079]图3分别给出了在SCE-47母线系统和SCE-56母线系统中应用本发明方法的计算结 果。为了比较，图中也给出了S0CP松弛方法的结果。从图中可以看出，本发明方法能够快速 收敛，在满足电力系统潮流方程、功率和电压约束的同时使系统总功率消耗随着迭代次数 不断减小直至收敛，而基于凹凸优化技术的最优潮流控制方法达到了和基于S0CP的最优潮 流控制方法几乎相同的目标值，说明本发明技术方案能够实现最优潮流控制。
[0080] 图4分别给出了在SCE-47母线系统和SCE-56母线系统中应用本发明方法后，模型 中约束可行性指标与迭代次数关系图。从图中可以看出，随着迭代次数的增加，约束条件逐 渐得到满足，且在迭代第2次后，本发明方法即可满足约束可行性指标。
[0081] 本发明不仅局限于上述【具体实施方式】，本领域一股技术人员根据本发明公开的内 容，可以采用其它多种具体实施方案实施本发明。因此，凡是采用本发明的设计结构和思 路，做一些简单的变化或更改的设计，都落入本发明保护范围。
【主权项】
1. 一种基于惩罚凹凸优化技术的电力系统最优潮流控制方法，其特征在于，该方法包 括如下步骤： 步骤1:获得电力系统网络参数:母线集合N和去除根母线后的集合Ν+=Ν\{0};电网支路 集合ξ;母线之间支路的阻抗zu，V(/，./) 各母线节点注入功率Sl的约束集合S 1;母线电压 幅度值的平方Vi的下限Xi和上限巧;夂·π in始化迭代次数k = 〇，最大迭代次数Kmax，收敛精度t〇1 ;设定初始点 淀惩罚系数β;?υ表示从母线i到母线j的电流幅度值的平方;S1^示母线i 到j间的连线发送端的潮流； 步骤3:利用泰勒展开，得到最优潮流问题的近似凸问题：SieSi,ieN+ V1 < v?. <Vrie Ni 其中|a|、石和Re(a)分别表示复数a的幅度、共辄和实部；表示母线i处注入功 率所引起的系统性能损耗； 步骤4:求解问题(Pl)得到当前的计算结j，判断是否达到收敛精度：SVo/或者是否达到最大迭代次数:k>Kmax;若是，输出母线的注入功率，计算系 统性能损耗，执行步骤5;否则令迭代次数k = k+l，重复步骤3和4。 步骤5:根据计算得到的母线注入功率完成最优潮流控制。2. 根据权利要求1所述的一种基于惩罚凹凸优化技术的电力系统最优潮流控制方法，其特征在于，所述步骤4中求解问题(Pl)得到计算结果 ?解决方法为内点 算法。
【文档编号】H02J3/06GK106026105SQ201610586607
【公开日】2016年10月12日
【申请日】2016年7月22日
【发明人】陈志勇, 赵明杰, 史清江, 徐伟强
【申请人】浙江理工大学