基于目标级联分析法的输电网结构优化

张 旋，张玉敏，吉兴全，杨子震，刘 健，杨 建

（1.山东科技大学 电气与自动化工程学院，山东 青岛 266590；2.国网青岛供电公司，山东 青岛 266002）

0 引言

在能源广域互联的背景下，电网互联有利于电能远距离传输，提高能源利用效率，进而为大容量、高效能机组的优化管理提供支持。然而，我国的能源与负荷呈逆向分布，源荷不平衡矛盾较为突出，而各类间歇性可再生能源发电占比日益增加，源网矛盾逐渐突出［1］，输电网运行过程中时常出现网络阻塞等问题［2］，使多区域互联电网间的调控能力面临严峻考验，因此增强电力系统区域间的互联和非同调效应具有重要意义。

输电网中的非同调现象，在不额外增加输电网建设投资的前提下，可以根据电网的实际工作条件，通过开断线路和电磁环网开环运行进行输电网结构优化，可在一定程度上消除网络阻塞，从而缓解了源、网、荷之间的冲突，提高了系统的灵活性和经济性。基于此，输电网结构优化（optimal transmission switching，OTS）成为一种通过改变电力系统中部分线路的开闭状态以改变潮流分布的有效方法［3］。目前已有部分学者对输电网的结构优化进行了研究。文献［4］提出了考虑结构优化的含可再生能源的OTS模型，通过OTS降低了系统运行成本及故障后网络阻塞造成的损失成本。文献［5］结合可再生能源，构建了计及网络结构优化和储能配置的两阶段随机优化模型，提高了电力系统的经济性。文献［6］提出了基于广义短路比灵敏度分析的OTS方法，其有助于协调系统强度与短路电流水平之间的矛盾。文献［7］在矩估计理论的基础上构建OTS模型，使得特高压接入下的电力系统安全性有所提高。由于考虑OTS的经济调度模型是典型的混合整数规划问题，模型中引入了大量的离散变量，使得模型求解难度增加，上述研究虽然验证了OTS能够有效提升系统的经济性与可靠性，但没有对提升模型的求解速度进行研究。

为此，诸多学者围绕模型的求解算法展开了深入研究。文献［8］通过使用交替方向乘子法（alternating direction method of multipliers，ADMM）实现交流系统的分布式最优潮流（optimal power flow，OPF）计算。由于ADMM源自增广拉格朗日乘子法和邻近点算法，它并没有消除一阶算法固有的缺点［9］。文献［10］采用分布式内点法对电网OPF进行分散式求解。为实现各个区域的同步迭代，文献［11］采用同步交替方向乘子法解决多子系统的并行协调优化问题，进而确定整体最优解。上述优化算法能有效提升模型的求解速度，但未考虑到电网规模庞大、受分区管理、信息不共享等因素影响，难以建立起统一的互联电网模型。

目标级联分析（analytical target cascading，ATC）法可用于加速难以用集中方式解决大规模优化问题的处理过程，并用于管理具有多个独立控制实体的系统。文献［12］利用ATC法寻找日前调度的最优发电计划，在每次迭代中，只有边界母线的电压和相角需要在主问题和子问题之间共享。因此，每次迭代的数据交换量低且通信简单。文献［13］利用ATC法搭建了双层调度框架，提出了一种面向主动配电网和虚拟微电网的线性化交互式调度模型，解决了由于上下层之间存在交互变量而导致的耦合问题。上述研究为本文开展输电网结构优化提供了良好的借鉴作用。

为提升互联电网整体运行水平和模型的计算效率，本文基于输电系统中的非同调现象，提出了一种基于ATC法的计及可再生能源的OTS模型。针对区域互联的输电网结构模型，通过母线撕裂法对互联电网进行切分，在切分后的各分区内进行独立的经济调度，利用ATC法进行解耦和求解，有效降低了求解难度。基于算例分析验证了所提模型和方法的有效性。

1 ATC法

ATC法首先将上层系统中的耦合变量传递给子系统，然后对各个子系统分别进行求解，从而获得全局 最 优 解［14］。ATC法 与ADMM和 辅 助 问 题 原 理（auxiliary problem principle，APP）法类似［15］，APP、ADMM应用了对偶概念并引入了一组罚函数，然后将原优化问题分解成若干个子问题；ATC法将整个系统分解成若干个子系统，然后应用约束松弛的概念对各问题进行并行求解［16］。本文选取区域间联络线交换功率和电压相角作为共享变量，通过设置目标变量和响应变量对相邻层间的共享变量进行建模，形成各独立区域相关的目标函数和约束条件，其中目标变量为上层（系统层）传送至下层（子系统层）的共享变量，响应变量为下层传送至上层的共享变量，具体如式（1）所示。

在式（1）所示的目标函数中，第2、3项是由子系统响应和耦合变量表示的偏差，子系统可表示为式（2）。

式中：gsub，i(·)和hsub，i(·)分别为子系统i的不等式约束和等式约束集合；xsub，i为子系统i的设计变量。

2 基于ATC法的输电网结构优化模型

考虑到可再生能源与电网结构优化的经济调度是一个复杂的非线性优化问题，基于母线撕裂法和ATC法构建OTS数学模型。为不失一般性，本文通过母线撕裂法将系统分成区域A—C这3个部分，分别建立相应的OTS模型。本节以区域A为例，构建基于ATC法的输电网多区域经济调度模型，同理可得区域B、C的数学模型。

2.1 目标函数

目标函数为系统常规机组的发电成本最小化，如式（3）所示。

2.2 约束条件

1）发电机出力约束。

2）爬坡约束。

3）支路功率约束。4）电压相角约束。

5）节点功率平衡约束。

6）区域耦合约束。

式中：ΔA为区域A的相邻区域集合；ΓA，B为区域A和B的联络支路集合；(i，j)∈ΓA，B表示区域A内节点i和区域B内节点j分别为联络支路的首、末节点；Pij和θij分别为联络支路的功率和相角差。

7）支路开断数量约束。

式中：JA为区域A中最大允许的线路开断数目。

8）风电和光伏出力约束。

2.3 模型简化

以图1所示的上、下两层系统结构为例，对本文所提模型进行分解。图中：上层系统为区域A，下层系统为区域B和区域C；t为从上层系统发送到下层系统的共享变量组成的向量，将其作为目标向量；r为下层系统向上层系统传输的共享变量组成的向量，将其作为响应变量。由图1（a）可见，为追求系统总运行成本的最小化，将上下层系统决策量进行迭代优化，上层系统向下分流，下层系统向上反馈，从而实现互联电网的整体优化。为便于描述ATC算法的求解过程，本文给出模型的矩阵表达式为：

图1 上、下层系统的分解Fig.1 Decomposition of upper- and lower-layer system

式中：gA、hA分别为上层系统中区域A需满足的不等式和等式约束；x为上层系统中除联络支路以外的决策变量，x、t1、t2、…、tn为上层系统的区域变量；、分别为第n个下层系统（例如图1（a）区域B）的不等式、等式约束条件；y为第n个下层系统的变量需满足的本地约束，y、r1、r2、…、rn、t′1、t′2、…、t′n为第n个下层系统区域变量；c表示上层系统区域变量与第n个下层系统区域变量的耦合约束。

3 模型转化与求解

3.1 模型解耦

为便于形成各区域相关的目标函数及约束条件，设置2组变量集合。根据ATC解算思想，将约束式（19）分别在上层和下层系统的模型中进行求解。这里利用罚函数ζ来松弛耦合约束，如式（20）所示。

式中：⊙表示Hadamard积；λ和μ分别为罚函数ζ的一次项和二次项的乘子向量。

上下层系统的优化模型可表示为：

通过式（20）将上、下层系统间的耦合约束进行松弛处理后，上层系统和下层系统只剩下区域变量中需要满足的本地约束和本地决策变量，由此上、下两层系统实现了解耦。

3.2 解耦的上层系统与下层系统优化模型

在基于ATC法的互联电网结构优化调度中，上层系统的目标函数为：

式中：λn、μn分别为拉格朗日一次项、二次项的乘子，通过不断更新拉格朗日乘子使目标函数满足收敛条件，以至趋于最优；为rn优化后的新变量。

上层的区域A在求解自身模型时，将虚拟负荷tn优化后的值以参数的形式传递给下层系统。上层系统在更新时，需要将虚拟负荷联动优化，并对每个区域的虚拟发电机进行优化。

由式（20）可见，与上层系统相同，在下层系统进行独立优化时，需要将虚拟发电机rn和虚拟负荷tn结合进行优化，通过引入罚函数对耦合约束进行松弛，并加入到下层系统的目标函数中。以区域B为例，下层系统的目标函数表示为：

通过不断更新乘子使目标函数满足收敛条件，以趋于最优，从而上、下层系统可独立进行更新、求解。

3.3 收敛判据与乘子更新原则

优化调度算法的收敛判据为：式（25）表示在第k次迭代中，下层系统虚拟发电机rnk与上层系统虚拟负荷tnk作为耦合变量，其差值需要满足所要求的精度ε1；式（26）表示分解后的上层系统和下层系统的总体效益在相邻2次迭代中需要满足精度ε2。若不能同时满足式（25）和式（26），则依据式（27）更新乘子。

式中：β的作用是为了加快收敛速度，一般取值为2≤β≤3；λ和μ的初值一般选取较小的常数，本文取值为λ=μ=1。

3.4 基于ATC法的OTS优化求解流程

基于ATC法的OTS优化求解流程如附录A图A1所示，具体步骤如下。

1）设定ATC法的最大迭代次数为k′，并置当前迭代次数k=1，输入常规机组参数，设置各个优化变量以及罚函数乘子等初值。

2）求解下层系统。根据下层系统模型，即式（24）和式（18）对每个区域进行并行优化，并将求解得到的功率传递给上层系统。

3）求解上层系统。根据上层系统模型，即式（23）和式（17）进行优化，并将求解得到的功率，传递给下层系统。

4）判断是否同时满足式（25）和式（26），若同时满足，则终止迭代过程并输出最优调度结果，否则根据式（27）更新乘子，置迭代次数k=k+1并返回步骤2）。如果k＞k′，则说明算法不收敛，终止计算。

4 算例分析

选择改进的IEEE 14节点系统及IEEE 118节点测试系统对模型进行仿真分析。计算机配置为intel i7-6300处理器、8 GB内存，采用GAMS软件进行编程，并选用CPLEX工具对模型进行求解。

4.1 IEEE 14节点系统

IEEE 14节点系统包含5台发电机、20条线路，机组和线路参数见附录A表A1—A3，风电场参数和光伏发电参数分别见附录A表A4和表A5。考虑到光伏与风力发电具有较大的间歇性与随机性，故选取其出力上、下限中的随机变量进行测试。通过母线撕裂法将该系统分解成三区域系统，如图2所示。设置区域A和区域B的开断输电线数为1。

图2 三区域互联系统Fig.2 Three regional interconnection system

4.1.1 ATC法有效性分析

为了验证所采用ATC法的求解优势，将ATC、APP法所得结果进行对比，如表1与图3所示。

表1 不同算法的求解结果Table 1 Results of different algorithms

图3 算法结果对比Fig.3 Comparison of algorithm results

由表1和图3可知：在运行成本方面，采用APP法求解的运行成本为 $ 18 212.23，而ATC法求得的运行成本为 $ 18 394.46，由于初始参数的设置与网络拓扑的改变，使得ATC法的运行成本略高于APP法；在收敛性能上，ATC法的迭代次数和迭代时间分别为34次和7.71 s，APP法在迭代次数和收敛时间上的表现劣于ATC法，说明了本文所采用的ATC法收敛性能较好。APP法的收敛性能取决于其参数和电网分区，由于OTS问题的复杂性，很难得到充分而有效的解；而ATC法可扩展性强，对系统的OTS优化问题具有较好的适用性。

4.1.2 输电网结构优化及可再生能源接入影响分析

为了证明OTS的效果和可再生能源发电对系统运行的影响，设置如下4种方案并进行对比分析：方案1，不考虑OTS及可再生能源发电；方案2，不考虑OTS，考虑可再生能源发电；方案3，考虑OTS，不考虑可再生能源发电；方案4，同时考虑OTS及可再生能源发电。

1）不考虑输电网结构优化的影响分析。

当不考虑OTS时，为了验证可再生能源接入对输电网运行的影响，对比分析了方案1、2下各区域运行成本和机组出力，结果如表2所示。

表2 方案1和方案2的运行结果对比Table 2 Operating result comparison between Scheme 1 and Scheme 2

由表2分析可知：在方案2中，由于考虑了可再生能源发电的接入，极大地释放了机组GA1的发电能力，使得经济性较好的机组GA1出力增加，经济性最差的机组GA2出力减少；采用方案1时，各区域成本分别为 $ 8 476.62、$ 7 363.38及 $ 4 363.69，而采用方案2时，各区域运行成本与采用方案1时相比分别降低了19.97 %、1.39 %和-0.13 %，其中区域A成本降低幅度较大，而区域B、C的成本变化幅度较小，这说明下层区域B、C能够凭借自身的发电资源实现分布自治，从而极大地减少了上层区域A的功率救济，使得区域A的经济性得到明显提高。

2）输电网结构优化影响分析。

为了证明可再生能源接入时考虑OTS的效果，方案2和方案4的机组出力和运行成本如表3所示，输电线路运行情况如表4所示。

表4 考虑OTS的输电线路运行情况Table 4 Transmission line operation condition considering OTS

由表3分析可知：与不考虑OTS的方案2相比，采用方案4时，机组GA1出力增加了43.27 %，经济性最差的机组GA2出力减少了35.27 %，这说明考虑OTS可实现源出力模式与网架结构的配合，使得发电资源合理分配；在经济性上，采用方案4时，各区域运行成本分别为 $ 5 814.34、$ 6 086.17和 $ 4 261.58，与采用方案2时相比分别下降了14.29 %、16.18 % 和2.47 %，由此说明将电网结构视为“动态可变”的经济调度方式能够达到缓解网络阻塞的目的，提升了系统的运行经济性。

表3 方案2和方案4的运行结果对比Table 3 Operating result comparison between Scheme 2 and Scheme 4

由表4分析可知：当考虑OTS时，系统可通过负荷模式对输电线路状态进行优化调节，调节过程中区域A内支路1-101处于停运状态，从而降低了传输阻塞程度，输电线路1-2的传输功率得以大幅度增加，进一步表明了风、光的接入及开断输电线路的措施同时应用于电网调度中，可有效提升系统运行的经济性，更大程度地实现源网协同调度；方案3和方案4的开断线路分别为2-101、202-2和1-101、202-2，这说明在考虑OTS后，输电线路的运行状态可以依据相应负荷模式进行适当调整，可有效提高电网运行的灵活性。

4.2 IEEE 118节点系统

IEEE 118节点测试系统结构图如附录B图B1所示，系统发电机组特性数据、输电线路数据以及系统负荷数据见文献［18］。

4.2.1 算法性能对比

为证明ATC法在大规模系统中的适用性，在IEEE 118节点测试系统中比较了ATC法和APP法，结果如图4所示。

图4 不同收敛精度的算法优化结果Fig.4 Optimization results of algorithms with different convergence accuracy

由图4可知：随着收敛精度的提高，ATC法趋于最优的速度最快，展现了良好的收敛性能；收敛精度的提高导致迭代次数也随之增加；相比于APP法，在任一精度要求下，ATC法始终保持良好的收敛性能和优化结果，这说明ATC法对于大规模系统仍具有较好的适用性。

4.2.2 优化结果对比

针对IEEE 118节点系统，分别采用上述4种方案对模型进行测试分析，所得的结果如表5所示，具体分析如下。

表5 4种方案的运行成本对比Table 5 Comparison of operating costs of four cases

方案4的总成本相较于方案1—3分别下降了18.99 %、17.65 %、2.50 %；三区域运行成本分别为$ 47 311.34、$ 63 384.13、$ 52 117.87，与方案1—3相比都有不同程度的降低。这说明在保证系统安全运行的前提下，通过电网结构调度，提升了可再生能源的接纳水平，使系统运行的经济性得以提高。方案4中区域A与区域B的运行成本较方案1—3的下降程度高于区域C，这是由于区域A、B阻塞较严重，故节省成本占比较为明显。综上所述，将风电、光伏系统及开断输电线路同时应用于电网中，实现系统区域间的协同合作，能够降低发电成本，提高能源利用效率，从而提升系统整体运行的灵活性和经济性。

5 结论

本文对多区域电力系统调度优化问题进行研究，提出了一种基于ATC法的输电网结构优化模型。算例分析所得结论如下。

1）所提模型能够有效提高互联电网整体运行的经济性，使电网灵活运行的能力得以充分发挥，高效利用全网资源，实现互联电网整体的统一优化。

2）所提模型增加了表示输电线路是否运行的离散变量，使其能够合理决策出系统所需的输电网架构。将电网结构优化、可再生能源与ATC相结合，极大地释放了网架结构的传输能力，提高了模型和方法的通用性。

3）本文所采用的ATC法收敛速度快，应用于多区域电网分布式求解时具有明显的效率优势，对较大规模的系统也具有较好的适应性。

后续研究可考虑如何提高模型的计算效率与计算精度。

附录见本刊网络版（http：//www.epae.cn）。