基于MATAB的电力系统潮流计算介绍杨浩(东华大学电气信息学院，上海)摘要:简单介绍了电力系统潮流计算的几种传统方法，重点讲述了在ATB开发环境下，应用电力系统仿真工具箱MATPWR实现对电力网络的潮流求解，并且给出了基本算例。关键词:潮流计算;ATPOWItrducton ope sytemflow caultion basd n MTLABAbtact: e everatrdionalmthd foersysm ow calculation a introduced iefly. ypicaly, h ower flow aculatioo pwr ntwork i implmented by MTPOWE， asliotoolbox ofper syste ine MATLAdevelpig enviroment. A n elemntay cacaioexaple is givnintis thesis.eyords: e flow aculain; TPOWER0引言电力系统潮流计算是研究电力系统稳态运行情况的一种基本电气计算。它的任务是根据给定的运行条件和网路结构确定整个系统的运行状态,如各母线上的电压（幅值及相角)、网络中的功率分布以及功率损耗等。它在电网规划、编制年运行方式、电力系统检修，以及事故预判、设备退出应用方面都有重要意义。在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时,为了实时监控电力系统的运行状态,也需要进行大量而快速的潮流计算。因此,潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。利用电子计算机进行潮流计算从世纪50年代中期就已经开始。此后,潮流计算曾采用了各种不同的方法，这些方法的发展主要是围绕着对潮流计算的一些基本要求进行的。对潮流计算的要求可以归纳为下面几点:算法的可靠性或收敛性；计算速度和内存占用量 ；计算的方便性和灵活性。电力系统潮流计算属于稳态分析范畴,不涉及系统元件的动态特性和过渡过程。因此其数学模型不包含微分方程，是一组高阶非线性方程。非线性代数方程组的解法离不开迭代，因此，潮流计算方法首先要求它是能可靠的收敛。随着电力系统规模的不断扩大，潮流问题的方程式阶数越来越高,目前已达到几千阶甚至上万阶,对这样规模的方程式并不是采用任何数学方法都能保证给出正确答案的。这种情况促使电力系统的研究人员不断寻求新的更可的计算方法。用数字计算机求解电力系统潮流的历史,经历了以节点导纳矩阵为基础的高斯-赛德尔迭代法，求解非线性方程式的典型方法牛顿拉夫逊法，以及其改进方法快速解耦Q分解法。近20多年来,潮流算法的研究仍然非常活跃,但是大多数研究都是围绕改进牛顿法和-Q分解法进行的。此外，随着人工智能理论的发展,遗传算法、人工神经网络、模糊算法也逐渐被引入潮流计算。但是，到目前为止这些新的模型和算法还不能取代牛顿法和PQ分解法的地位。本文将主要阐述牛顿法和P分解法的MATLAB实现,进而分析其性能的可靠性。1 百度百科，电力系统潮流计算，2011/11/5访问;2 陈珩，电力系统稳态分析，中国电力出版社，2007;3 徐劲松、宁玉琳、杨永峰，基于Matlab的电力系统PQ分解法潮流计算研究，电气传动自动化，2011，33(2):1018;4 于群、曹娜，MATLAB/Simulink电力系统建模与仿真，机械工业出版社，2011.电力系统潮流计算1.1电力系统潮流计算模型电力系统主要由发电机、变压器、输电线路和负荷构成,它们都是具有物理意义的电气元件。而电力系统分析中,需要对这些元件进行等值数学模型的建立,才能实现稳态下量的计算分析。而电力系统元件根据各自的性能均可以建立R、C等值线性模型,变压器可以等效成型网络，输电线路可以等效成型网络等，如下图。 变压器型等值网络输电线路型等值网络针对电力系统等值网络可以列些节点电压方程:式中为导纳矩阵，为节点注入电流矩阵，为节点电压列向量。在工程实践中,通常已知的既不是节点电压，也不是节点电流，而是各节点的功率和。所以一般建立的是功率方程:展开式为:上式含有n个非线性复数方程,是潮流计算问题的基本方程,若用极坐标形式表示,即: (）这样电力系统潮流计算的问题，便成为了根据非线性功率方程解各节点有功功率、无功功率、电压幅值和电压相角的问题。电力系统中的节点一般分为三种类型：节点:已知的是节点注入功率、（包括负荷功率和电源功率）,待求的是节点电压幅值和电压相角，系统中大部分是这类节点;节点:已知节点的负荷有功功率和电源有功功率，以及负荷的无功功率、电压幅值,待求的是电源的无功功率和电压相角；平衡节点：潮流计算时,一般只设一个平衡节点，已知负荷功率和，节点电压幅值和电压相角，待求的是电源功率、。1.2电力系统潮流计算方法高斯赛德尔迭代法: 高斯赛德尔迭代法是在电力系统中应用最早的一种潮流计算方法，初值选定后，计算每个变量的估计值，每一步迭代过程都利用前一步已求得的最新估计值,直接代入求解下一个变量的估计值，而不是等待一次迭代过程结束后再利用。一个迭代过程结束后得到的每个变量的估计值,再重新代入公式计算变量估计值，如此反复迭代，不断得到新的估计值，直到得出满意的结果为止。这种方法编程简单、占用内存少,但算法收敛性差、计算时间长，后来被牛顿法所取代。牛顿拉夫逊法：牛顿拉夫逊法是求解非线性代数方程有效的迭代计算方法，已经成为求解电力系统潮流问题应用最为广泛的一种方法。而牛顿型潮流计算的核心问题是修正方程式的建立和求解。将式中的电压相量改写成复数形式:并将实部、虚部分别对应： （）此外，由于系统中还有电压大小给定的节点，还应补充一组方程式：将（2)式改写成残差形式： （3）对上式中的、分别求导,便可以得到修正方程组的雅可比矩阵: (4)应用牛顿拉夫逊法计算潮流的基本步骤有:1. 形成节点导纳矩阵;2. 设各节点电压的初值、;3. 将各节点电压的初值代入式（3)求修正方程式中的不平衡量、以及;4. 将各节点电压的初值代入式（4)求修正方程式的系数矩阵雅可比矩阵的各个元素、以及、；5. 解修正方程式,求各节点电压的变化量,即修正量、;6. 计算各节点电压的新值，即修正后值：；7. 运用各节点电压的新值自第三步开始进入下一次迭代；8. 计算平衡节点功率和线路功率。目前,牛顿法因为具有收敛快、内存较少的优点而普遍应用，但是雅克比矩阵计算量大、编程复杂。PQ分解法：针对牛顿法的缺点,在其基础上又提出了快速解耦P分解法。其基本思想是：把节点功率表示为电压向量的极坐标方程式，抓住主要矛盾,把有功功率误差作为修正电压向量角度的依据,把无功功率误差作为修正电压幅值的依据，把有功功率和无功功率迭代分开进行。它密切地结合了电力系统的固有特点，无论是内存占用量还是计算速度方面都比牛顿-拉夫逊法有了较大的改进。简单、快速、节省内存和收敛可靠成为该算法的突出优点,是当前国内外优先使用的算法。与牛顿法相比，分解法具有如下特点：1. 以一个(n-1)阶和一个（m1)阶系数矩阵、替代原有的(+m2）阶系数矩阵,提高了计算速度,降低了对贮存容量的要求;2. 以迭代过程中保持不变的系数矩阵、替代起变化的系数矩阵，显著地提高了计算速度;3. 以对称的系数矩阵、替代不对称的系数矩阵,使求逆等运算量和所需的储存容量都大为减少。2基于MTLA的潮流分析MATPOR简介MAPOER是一个基于MALAB的M文件的组建包,用来解决电力潮流和优化潮流的问题，是由卡奈尔大学电气学院电力系统工程研究中心RAYD.ZIMMENRan，ARLOS.Murlo和甘德强在ROET TOMS的指导下开发出来的。它的特点是简单、易懂且代码公开。APOER所用的所有数据文件均为MLA的M文件或者AT文件,他们用来定义和返回变量：baseMV，bu,rn,gen，areas和gencost。变量baseV 是标量,其他的都是矩阵。矩阵的每一行都对应于一个单一的母线,线路或者发电机组。列的数据类似于标准的IEE和TI列的数据格式。使用MATPWE首先需要对电力网络结构的数学模型,建立相应的MATPR数据文件,其中包括母线(Bus)数据矩阵、发电机（Geratr)数据矩阵以及支路（Bran）数据矩阵，实现对电力系统具体网络,到数学图论模型,再到数字化处理模型的建立。MAPOWER可以采用潮流计算的四种方法:牛顿法、快速解耦（X)、快速解耦（X）、高斯赛德尔法,通过改变变量“PF_AG”进行选择。此外,MTWER还可以实现直流潮流、最优潮流的计算仿真,在这里主要通过具体实例来介绍复杂网络的潮流计算。22具体算例分析(IE -butest stm）对上图所示的IEEE 9bus st systm,按照MTPOWER的数据文件格式编写数据M文件如下:untion mpc=case9eee%MATOER Cas Fm：Verson2pc.version=2；%- PrFl Data-% syte MV basmpc.baseMA=100;%us data%s_ te d Gs B areaVm V seK oemVnmc.us 3 0 0 00 1 00 .1 0.9; 20 0 00 1 1 0 10 1 1.1 0.; 3 20 0 0 1 0 100 1 1.10.9; 41 0 0 0 1 0 100 .19; 5 1 15 50 001 0 100 11.1 9； 6 9 3000 1 1 0 00 1 1.1 0.; 7 0 0 0 0 1 1 0 100 1 1.1 0; 8 1 0 35 0 0 1 1 010 1. 09； 9 0 0 1 100 1 .1 0; ;%geneatora%bus Pg Q Qmax Qmin Vgas Statu Pmax minmpc.ge= 10 0 30 -30 1 1001 25 0; 2 13 300 00 1 10 1 20 1; 3 0 00 -300