基于内点法的最优潮流计算

摘要内点法是一种能在可行域内部寻优的方法，即从初始内点出发，沿着中心路径方向在可行域内部直接走向最优解的方法。其中路径跟踪法是目前最具有发展潜力的一类内点算法，该方法鲁棒性强，对初值的选择不敏感，在目前电力系统优化问题中得到了广泛的应用。本文采用路径跟踪法进行最优求解，首先介绍了路径跟踪法的基本模型，并且结合具体算例，用编写的Matlab程序进行仿真分析，验证了该方法在最优潮流计算中的优越性能。关键词：最优潮流、内点法、路径跟踪法、仿真电力系统稳态分析课程文论第Ⅰ页共Ⅰ页目次0、引言..........................................................................................................................11、路径跟踪法的基本数学模型..................................................................................22、路径跟踪法的最优潮流求解思路..........................................................................33、具体算例及程序实现流程......................................................................................63.1、算例描述........................................................................................................63.2、程序具体实现流程........................................................................................74、运行结果及分析....................................................................................................114．1运行结果.....................................................................................................114．2结果分析......................................................................................................155、结论........................................................................................................................166、编程中遇到的问题................................................................................................17参考文献......................................................................................................................19附录..............................................................................................................................20电力系统稳态分析课程论文第1页共30页0、引言电力系统最优潮流，简称OPF（OptimalPowerFlow）。OPF问题是一个复杂的非线性规划问题，要求满足待定的电力系统运行和安全约束条件下，通过调整系统中可利用控制手段实现预定目标最优的系统稳定运行状态。针对不同的应用，OPF模型课以选择不同的控制变量、状态变量集合，不同的目标函数，以及不同的约束条件，其数学模型可描述为确定一组最优控制变量u，以使目标函数取极小值，并且满足如下等式和不等式。{𝑚𝑖𝑛𝑢𝑓(𝑥,𝑢)𝑆.𝑡.ℎ(𝑥,u)=0𝑔(𝑥,𝑢)≤0（0-1）其中𝑚𝑖𝑛𝑢𝑓(𝑥,𝑢)为优化的目标函数，可以表示系统运行成本最小、或者系统运行网损最小。𝑆.𝑡.ℎ(𝑥,u)=0为等式约束，表示满足系统稳定运行的功率平衡。𝑔(𝑥,𝑢)≤0为不等式约束，表示电源有功出力的上下界约束、节点电压上下线约束、线路传输功率上下线约束等等。电力系统最优潮流算法大致可以分为两类：经典算法和智能算法。其中经典算法主要是指以简化梯度法、牛顿法、内点法和解耦法为代表的基于线性规划和非线性规划以及解耦原则的算法，是研究最多的最优潮流算法,这类算法的特点是以一阶或二阶梯度作为寻找最优解的主要信息。智能算法主要是指遗传算法和模拟退火发等，这类算法的特点是不以梯度作为寻优信息，属于非导数的优化方法。因此经典算法的优点是能按目标函数的导数信息确定搜索方向，计算速度快，算法比较成熟，结果可信度高。缺点是对目标函数及约束条件有一定的限制，可能出现局部极小时难以收敛。而智能算法的优点是计算与导数无关，灵活性高，随机性强，缺点是算法不稳定，结果不可信，并且控制参数需凭经验给出。通过对这些常见算法的简单比较，内点法具有其优越的性能，特别是路径跟踪法，其算法收敛迅速，鲁棒性强，对初值的选择不敏感，其迭代次数与系统规模或控制变量的数目关系不大，因此本文采用该方法进行最优计算。电力系统稳态分析课程论文第2页共30页1、路径跟踪法的基本数学模型内点法最初的基本思路是希望通过寻优迭代过程始终在可行域内进行，因此，初始点应在可行域内，并在可行域的边界设置‘障碍’使迭代点接近边界时其目标函数迅速增大，从而保证迭代点均在可行域的内点。但是对于大规模实际问题而言，寻找初始点往往十分困难。为此许多学者长期以来致力于内点算法初始‘内点’条件的改进。以下介绍的路径跟踪法只要求在寻优过程中松弛变量和拉格朗日乘子满足简单的大于0或者小于0的条件，即可代替原来必须在可行域内求解的要求，使得计算过程大为简化。一般可以将最优潮流模型简化为如下的非线性优化模型。Obj.𝑚𝑖𝑛𝑢𝑓(𝑥,𝑢)（1-1）s.t.𝑆.𝑡.ℎ(𝑥,u)=0（1-2）𝑔−≤𝑔(𝑥,𝑢)≤𝑔−（1-3）其中𝑚𝑖𝑛𝑢𝑓(𝑥,𝑢)为优化的目标函数，𝑆.𝑡.ℎ(𝑥,u)=0为等式约束，𝑔(𝑥,𝑢)为不等式约束，路径跟踪内点法的基本思路是：首先将式（1-3）的不等约束变成等式约束：𝑔(𝑥,𝑢)+u=𝑔−（1-4）𝑔(𝑥,𝑢)−l=𝑔−（1-5）其中松弛变量l=[𝑙1,…,𝑙𝑟]𝑇，u=[𝑢1,…,𝑢𝑟]𝑇，应满足u0，l0这样原问题就转化为问题A：Obj.𝑚𝑖𝑛𝑢𝑓(𝑥,𝑢)S.t.()0()()hxgxuggxlg然后，把目标函数改造成障碍函数，该函数在可行域内应接近于原函数f(x)，而在边界时变得很大。一次可得带优化问题B：obj.11min.()log()log()rrrrjjfxuluu电力系统稳态分析课程论文第3页共30页s.t.()0()()hxgxuggxlg其中扰动因子或者障碍因子u0。当l或u接近边界时，以上函数将趋于无穷大，因此满足以上障碍目标函数的极小解不可能在边界上找到。这样就通过目标函数的变化把含不等式限制的优化问题A变成只含等式限制优化的问题B了，因此可以直接用拉格朗日乘子法来求解。优化问题B的拉格朗日目标函数为：11()()[()][()]log()log()rrTTTrrjjLfxyhxzgxlgwgxuguluu式中：y,z和w均为拉格朗日乘子。因此最后简化的求解问题就是求取上述表达式的极小解。2、路径跟踪法的最优潮流求解思路路径跟踪法的最优潮流求解过程就是对拉格朗日目标函数求极小值问题：11()()[()][()]log()log()rrTTTrrjjLfxyhxzgxlgwgxuguluu式中：y,z和w均为拉格朗日乘子。该问题极小值存在的必要条件是拉格朗日函数对所有变量及乘子的偏导数为0。即：11()()()(+w)0()0()0()+u00U0xxxxyzwlluuLLfxhxygxzxLLhxyLLgxlgzLLgxgwLLzLeLLZeelLLweLUWeeu（2-1）通过上述表达式可以解出：电力系统稳态分析课程论文第4页共30页2TTlzuwr（2-2）定义：TTGaplzuw，称为互补间隙。可得：2Gapr（2-3）如果x*是优化问题A的最优解，当u固定时，x(u)是优化问题B的解，那么当Gap→0，u→0时，产生的序列{x(u)}收敛至x*。建议采用：2Gapr。式中(0,1)称为中心参数，一般取0.1，在大多数场合可获得较好的收敛效果。通过偏导数为0的表达式可以可得内点法的修正方程为：111'0000000()000000000()00000'()0000()0lTzxuTwxxxTyxLLzILZLlIgxULwILuIgxLxHhxLyhx（2-4）求解方程可得到第k次迭代的修正量，于是最优解的一个新的近似解为：(1)()(1)()(1)()(1)()(1)()(1)()kkpkkpkkpkkdkkdkkdxxaxllaluuauyyayzzazwwaw（2-5）式中，pa和da为步长：0.9995min{min(,0;,0),1}iipiiiilualulu0.9995min{min(,0;,0),1}iidiiiizwazwzw（2-6）其潮流计算的流程图如下图1所示，其中初始化部分包括：电力系统稳态分析课程论文第5页共30页(1)、设置松弛变量l和u，保证[l，u]T0(2)、设置拉格朗日乘子w、y、z，满足[w0,z0，Y！=0]T(3)、设置优化问题的初值。(4)、取中心参数(0,1)，给定计算精度，迭代次数初值K=0。图1.内点法潮流计算流程图输出最优解，停止计算计算互补间隙Gap计算扰动因子计算和更新原始变量及拉格朗日乘子kKmGap求解修正方程，求出，，，，，输出“计算不收敛”初始化是是否否电力系统稳态分析课程论文第6页共30页3、具体算例及程序实现流程这部分主要有算例描述以及程序的实现流程两部分，其中算例描述主要是对系统参数以及优化问题进行说明。而程序的实现流程主要描述的是最优潮流计算中所涉及的矩阵方程的描述。3.1、算例描述该算例为王锡凡编写的《现代电力系统分析》中的3-1的例题，是以系统燃料最省为最优潮流的目标函数。选择该题目作为算例分析的原因是，该题目有比较详细的解题思路以及列写出了比较详细的迭代结果，方便对编写程序的运行结果进行比对。求如下图所示简化系统的系统燃料最省的最优潮流计算：除了由上图所提供的系统母线负荷功率数据、线路参数和变压器之路参数数据、变压器