《基于非线性预报-校正内点法的电力系统》由会员分享，可在线阅读，更多相关《基于非线性预报-校正内点法的电力系统（7页珍藏版）》请在金锄头文库上搜索。
1、基于非线性预报基于非线性预报- -校正内点法的电力系校正内点法的电力系 统统 无功优化研究无功优化研究吴 锴（天津市电力公司城西供电分公司 天津 300113）A NONLINEAR PREDICTOR-CORRECTOR INTERIOR POINT METHOD FOR REACTIVE POWER OPTIMIZATION IN POWER SYSTEMABSTRACTABSTRACT The predictor-corrector technology is applied to nonlinear primal-dual interior point method and a lon
2、ger iteration step can be obtained by use of the improved nonlinear predictor-corrector primal-dual interior point method than pure primal-dual interior point method, so the convergence can be speeded up. Using the improved method for power system reactive power optimization a large number of inequa
3、lity constraints in the objective function can be effectively dealt with. The simulation results from four test systems, i.e., IEEE14-bus system, IEEE 30-bus system, IEEE 57-bus system and IEEE 118-bus system, show that the improved algorithm possesses fast convergence speed and good robustness.KEYK
4、EY WORDSWORDS Power system；Reactive power optimization；Predictor-corrector interior point method；Nonlinear primal-dual interior point method摘摘 要要 在非线性原-对偶内点法的基础上引入了预报-校正技术，使改进后的非线性预报校正内点法获得了较纯原-对偶内点法更大的迭代步长，从而加速了算法的收敛。应用该方法求解电力系统无功优化问题时能有效处理目标函数中的大量不等式约束。IEEE 14 节点、IEEE 30 节点、IEEE 57节点和IEEE 118 节点系统
5、的仿真结果表明，该算法收敛快、鲁棒性好。关键词关键词 电力系统；无功优化；预报-校正内点法；非线性原-对偶内点法1 1 引言引言电力系统无功优化问题在数学上是一个多变量、高维数、多约束、离散和连续变量共存的非凸、非线性规划问题，可利用逐次线性规划法、逐次二次规划法、非线性规划法、牛顿法等来求解1,2。但这些方法收敛较慢，而且对不等式约束的处理都不太理想，甚至可能出现数值稳定问题。1984 年AT&T贝尔实验室的印度数学家Karmarkar 提出了一种新的内点法，该算法的速度极快，引起了优化问题研究领域大批学者的关注。1991 年Clements 首先将内点法用于求解电力系统状态估计问题，此后，
6、内点法在电力系统优化运行中得到了广泛应用。文献3将线性内点法应用于大规模电力系统的优化问题，求解了巴西一个1832 节点和北美一个3467 节点系统的无功调度问题。文献4应用原对偶仿射尺度内点法求解无功优化问题，该方法属于线性规划法的范畴，需将目标函数和约束条件线性化。文献5,6用二次内点法求解无功优化问题，保留了目标函数的非线性特性，但约束条件仍需线性化。由于无功优化问题本身具有强非线性特性，因此无论其目标函数还是约束条件线性化后必然引起误差。文献7比较了应用线性和非线性内点法求解无功优化问题的执行效果，并指出非线性内点法的优化效果优于线性内点法。因此，近年来人们都致力于非线性内点法的研究。
7、文献8将同伦方法与非线性内点法相结合，成功求解了几种极端情况下的无功优化问题，但由于引入了同伦变量，故算法的计算用时较长，几乎为内点法的2 倍，而且只在小型系统上验证了方法的有效性。随着互联电网规模的扩大，对大型电力系统全局无功优化算法的计算速度和收敛性的要求也越来越高。本文所采用的非线性预报校正内点法将问题的求解分成预报和校正两个阶段进行，首先在预报阶段求出仿射方向，然后在校正阶段对其进行校正从而得到牛顿方向。该算法能获得较大的迭代步长，具有良好的收敛特性和鲁棒性，近年来深受人们的青睐。本文将该算法应用于电力系统的无功优化，仿真结果表明了该算法的有效性和实用性。2 2 无功优化的数学模型无功
8、优化的数学模型电力系统无功优化问题的数学模型可表示为：(1) hxhhxgtsxf)(0)(. .)(min式中，X为包括控制变量和状态变量的N维列向量；N为变量数目；目标函数为全系统网损f(x)最小。其中 niijijjijiijVVVVGxf122)cos2()(式中：n为系统节点数目，表示节点与节点 之间有线路相连；、分别为节点ijjiiVjV、的电压幅值；为节点 、之间的电导；为节点 、的电压相角差。ijijGijijij为等式约束，即潮流平衡方程0)(xg njijijijijjiLiGiiBGVVPPP10)sincos(2-5) njijijijijjiLiCiGiiBGVVQQ
9、QQ10)cossin(（ =1,2,n）i式中：、分别为节点 处的发电机有功和无功出力，若该节点未接发电机则其值为GiPGiQi零；、分别为节点 的有功和无功负荷，若该节点处无负荷则其值为零；为节LiPLiQiciQ点 处的无功补偿装置的无功输出，若该节点处没有无功补偿装置则其值为零；为节点iijG，之间的电导，为节点 ，之间的电纳。ijijBij本文考虑的不等式约束条件有：发电机节点无功出力的约束，STATCON 的电流幅值约束，电容器电纳的约束，节电电压幅值的约束，可调变压器变比的约束，各支路传输功率约束，不等式约束表示如下：hxhh)(个GiGiGiQQQgn个CiiCCiIIIIcn
10、个 (4-3)CiCiCiBBBBcnn 个iiiVVV个ijijijTTTTS个lijiijijijijjiijlPGVBGVVPP2)sincos(lS本文的主要目的是突出优化算法，因此对离散变量未做像文献9那样的特殊处理，而是仍将离散变量当作连续变量处理。3 3 预报预报校正内点法校正内点法预报校正内点法的思想是由Kojima、Mizuno和Yoshise于1989年提出的，后来Mehrotra对该方法做了大量研究并在1992年取得成功。与原对偶内点法相比该算法在每次迭代中只增加了一次前代回代计算，但可以明显减少收敛次数，优化速度明显提高。1994年Wu等10首次将该方法应用于求解电力系
11、统最优潮流问题。预报校正内点法获得成功应用的关键是：将增加了高阶信息的牛顿系统分两步来计算，第一步是预报，得到仿射方向；第二步是校正，得到校正方向，从而使牛顿方向可以高阶近似逼近中心路径。两步计算用的是同一个系数矩阵，只需进行一次因子化，与常规内点法相比只在每步迭代增加了一次前代回代计算。引入松弛变量，为全部不等式约束的个数，将式(2-1)中rrRuRlul, 0, 0r的不等式约束转化为等式约束，并将目标函数改造为障碍函数，可得：(4-4) rjrjjjulxf11)lnln()(min. .ts0)(xg0)(huxh0)(hlxh式中，为扰动因子（或称障碍参数） 。显然，式(4-4)是只
12、含等式约束的优化问题，可以直接用拉格朗日乘子法求解，其拉格朗日函数为：)()()lnln()(),(11hlxhzxgyulxfwulzyxLTTrjjrjj )(huxhwT(4-5)式中，为式(2-1)中等式约束的个数；，为拉格朗mRymrrRwRz,wzy,日乘子，亦称对偶变量，且。由 KKT(Karush-Kuhn-Tucker)一阶必要条件可知：0, 0wz(4-6)0)(0)()()(0)(00)(011xgLwzxhyxgxfLhuxhLeUWeweULhlxhLeLZeLzeLLyxxxxwuzll式中，为各元素均为 1 的维列向量；er；),.,(21rllldiagL ),
13、.,(21ruuudiagU ),.,(21rzzzdiagZ 。将式(4-6)中的各等式展成泰勒级数得到：),.,(21rwwwdiagW (4-7)()()()()()()()()()()()()()(222xgxxgyxgwzxhxfyxgwzxhxyxgwzxhxfhuxhuxxhWeUeUWeuWwUhlxhlxxhZeLeLZelZzLT xxxxxxxxxT xT x 将式(4-7)写成矩阵形式为：= =0)(0000)(0)(0)(0)(00000000)(0000000xgxgHxhxhxhIWUxhIZLT xxxxT xT xyxuwlz+ + (4-8)()()(xgL
14、huxhUWehlxhLZex 0000 WeUeZeLe式(4-8)即为预报-校正内点法的修正方程式，其中：yxgwzxhxfHxxx)()()(222)()()(wzxhyxgxfLxxxx记分别为式（1-1）等号右边的第 1、第 2 项，则式（1-1）可表示为：coafDD ,（1-2）coafDDM2在预报阶段先由线形方程（1-3）afafDM2解出仿射方向，然后计算仿射迭代步长afaf（1-4）dpaf,min式中， 、分别为原变量和对偶变量的迭代步长：pd 1 ,minmin00iiuiilii puu ll 1 ,minmin00iiwiizii dww zz式中，为安全因子，通常取=0.9995。然后计算仿射间隙：（1-)()()()(afafT afafafafT afafafwwuuzzll5）并计算仿射障碍参数（1- 2 . 0 ,)(min22 afaf afb6）式中， ；wuzlTT于是
