基于MATLAB的电力系统潮流计算开题报告

1、1 文献综述电力系统潮流计算是对复杂电力系统正常和故障条件下稳态运行状态的计算。其目的是求取电力系统在给定运行方式下的节点电压和功率分布，用以检查系统各元件是否过负荷、各点电压是否满足要求、功率分布和分配是否合理以及功率损耗等，是电力系统计算分析中的一种最根本的计算1。潮流计算是电力系统的各种计算的根底，同时它又是研究电力系统的一项重要分析功能，是进行故障计算，继电保护鉴定，平安分析的工具。电力系统潮流计算是计算系统动态稳定和静态稳定的根底。在电力系统规划设计和现有电力系统运行方式的研究中，都需要利用电力系统潮流计算来定量的比较供电方案或运行方式的合理性、可靠性和经济性1。对于正在规划的电力系

2、统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置整定计算、电力系统故障计算和稳定计算等提供原始数据。潮流计算的目的在于:确定是电力系统的运行方式；检查系统中的各元件是否过压或过载；为电力系统继电保护的整定提供依据；为电力系统的稳定计算提供初值，为电力系统规划和经济运行提供分析的根底。因此，电力系统潮流计算是电力系统中一项最根本的计算，既具有一定的独立性，又是研究其他问题的根底1。1.2 国内外开展现状利用电子计算机进行潮流计算从20世纪50年代中期就已经开始。此后，潮流计算曾采用了各种不同的方法，这些方法的开展主要是围绕着对潮流计算的一些根本要求进

3、行的。对潮流计算的要求可以归纳为下面几点： 1算法的可靠性或收敛性 2计算速度和内存占用量 3计算的方便性和灵活性 电力系统潮流计算属于稳态分析范畴，不涉及系统元件的动态特性和过渡过程。因此其数学模型不包含微分方程，是一组高阶非线性方程。非线性代数方程组的解法离不开迭代，因此，潮流计算方法首先要求它是能可靠的收敛，并给出正确答案。随着电力系统规模的不断扩大，潮流问题的方程式阶数越来越高，目前已到达几千阶甚至上万阶，对这样规模的方程式并不是采用任何数学方法都能保证给出正确答案的。这种情况促使电力系统的研究人员不断寻求新的更可靠的计算方法2。 1.2.1 高斯-赛德尔迭代法在用数字计算机求解电力系

4、统潮流问题的开始阶段，人们普遍采用以节点导纳矩阵为根底的高斯-赛德尔迭代法一下简称导纳法。这个方法的原理比较简单，要求的数字计算机的内存量也比较小，适应当时的电子数字计算机制作水平和电力系统理论水平，于是电力系统计算人员转向以阻抗矩阵为主的逐次代入法以下简称阻抗法。 20世纪60年代初，数字计算机已经开展到第二代，计算机的内存和计算速度发生了很大的飞跃，从而为阻抗法的采用创造了条件。阻抗矩阵是满矩阵，阻抗法要求计算机储存表征系统接线和参数的阻抗矩阵。这就需要较大的内存量。而且阻抗法每迭代一次都要求顺次取阻抗矩阵中的每一个元素进行计算，因此，每次迭代的计算量很大3。 阻抗法改善了电力系统潮流计算

5、问题的收敛性，解决了导纳法无法解决的一些系统的潮流计算，在当时获得了广泛的应用，曾为我国电力系统设计、运行和研究作出了很大的奉献。但是，阻抗法的主要缺点就是占用计算机的内存很大，每次迭代的计算量很大。当系统不断扩大时，这些缺点就更加突出。为了克服阻抗法在内存和速度方面的缺点，后来开展了以阻抗矩阵为根底的分块阻抗法。这个方法把一个大系统分割为几个小的地区系统，在计算机内只需存储各个地区系统的阻抗矩阵及它们之间的联络线的阻抗，这样不仅大幅度的节省了内存容量，同时也提高了计算速度4。1.2.2 牛顿-拉夫逊法和P-Q分解法克服阻抗法缺点的另一途径是采用牛顿-拉夫逊法以下简称牛拉法。牛拉法是数学中求解

6、非线性方程式的典型方法，有较好的收敛性。解决电力系统潮流计算问题是以导纳矩阵为根底的，因此，只要在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿潮流程序的计算效率。自从20世纪60年代中期采用了最正确顺序消去法以后，牛拉法在收敛性、内存要求、计算速度方面都超过了阻抗法，成为直到目前仍被广泛采用的方法。 在牛拉法的根底上，根据电力系统的特点，抓住主要矛盾，对纯数学的牛拉法进行了改造，得到了P-Q分解法。P-Q分解法在计算速度方面有显著的提高，迅速得到了推广5。 牛拉法的特点是将非线性方程线性化。20世纪70年代后期，有人提出采用更精确的模型，即将泰勒级数的高阶项也包括进来，希望以

7、此提高算法的性能，这便产生了保存非线性的潮流算法。另外，为了解决病态潮流计算，出现了将潮流计算表示为一个无约束非线性规划问题的模型，即非线性规划潮流算法6。 近20多年来，潮流算法的研究仍然非常活泼，但是大多数研究都是围绕改良牛拉法和P-Q分解法进行的。此外，随着人工智能理论的开展，遗传算法、人工神经网络、模糊算法也逐渐被引入潮流计算。但是，到目前为止这些新的模型和算法还不能取代牛拉法和P-Q分解法的地位。由于电力系统规模的不断扩大，对计算速度的要求不断提高，计算机的并行计算技术也将在潮流计算中得到广泛的应用，成为重要的研究领域7。通过几十年的开展，潮流算法日趋成熟。近几年，对潮流算法的研究仍

8、然是如何改善传统的潮流算法，即高斯-塞德尔法、牛拉法和快速解耦法。牛拉法，由于其在求解非线性潮流方程时采用的是逐次线性化的方法，为了进一步提高算法的收敛性和计算速度，人们考虑采用将泰勒级数的高阶项或非线性项也考虑进来，于是产生了二阶潮流算法。后来又提出了根据直角坐标形式的潮流方程是一个二次代数方程的特点，提出了采用直角坐标的保存非线性快速潮流算法8。1.2.3 基于MATLAB 的电力系统潮流计算开展前景MATLAB 自1980 年问世以来，以其学习简单、使用方便以及其它高级语言所无可比较的强大的矩阵处理功能越来越受到世人的关注。目前，它已成为国际控制界最流行、使用最广泛的语言了。它的强大的矩

9、阵处理功能给电力系统的分析、计算带来许多方便。在处理潮流计算时，其计算机软件的速度已无法满足大电网模拟和实时控制的仿真要求，而高效的潮流问题相关软件的研究已成为大规模电力系统仿真计算的关键。随着计算机技术的不断开展和成熟,对MATLAB 潮流计算的研究为快速、详细地解决大电网的计算问题开辟了新思路。MATLAB 语言允许用户以数学形式的语言编写程序, 其比BASIC 语言和FORT RAN 等更为接近书写的数学表达格式, 且程序易于调试。在计算要求相同的情况下, 使用MATLAB 编程, 工作量将会大为减少9。基于MATLAB 的电力系统潮流计算使计算机在计算、分析、研究复杂的电力系统潮流分布

10、问题上又前进了一步。矩阵输入、输出格式简单, 与数学书写格式相似; 以双精度类型进行数据的存储和运算, 数据精确度高，能进行潮流计算中的各种矩阵运算, 包括求逆、求积和矩阵L R 分解等, 其程序的编写也因MATLAB 提供了许多功能函数而变得简单易行。另外, MATLAB稀疏矩阵技术的引入, 使电力系统潮流计算由传统方法转变为优化算法成为可能10。参考文献1 陈衍.电力系统稳态分析M.北京：中国电力出版社，2007.2 李维波. MATLAB 在电气工程中的应用M. 北京: 中国电力出版社，2007.3 武晓朦，张飞廷. 电力系统的P- Q 分解法潮流计算J.现代电子技术,2002 ,142

11、(11): 105-106 4 拉夫逊的潮流就算方法J.许昌学院学报. 2006.3,25(2):27-305 杨帆.电力系统潮流计算程序设计J.山西冶金，2007，1062：42-446 刘军,刘学军. MATLAB 在电力系统分析中的应用J. 电力系统及其自动化学报, 2000.4,12(2): 23-257 李家坤,刘姣姣. MATLAB仿真技术在电力电子教学中的应用J. 长江工程职业技术学院学院2021.3,26(1): 75-778 连小洲，焦洁. MATLAB在电力系统计算中的应用J. 江西电力职业技术学院学报，2005.12,18(4): 10-119 张宁,江红梅,张渭. 基于

12、MATLAB的电力系统潮流计算J.西北农林科技大学学报(自然科学版)，2004.12,32(12): 124-12610 周卫星，张颖.基于Matlab的电力系统潮流计算J.科技咨询导报，2007，10：70-7111 徐一哲，沈瑞寒.基于Matlab的电力系统潮流分析J.中外企业家，2021,5(下)：206-2082 本课题要研究或解决的问题和拟采用的研究手段及途径u 掌握潮流计算的手工算法,开式网络的潮流计算，简单闭合环网的潮流计算，配电网的潮流计算等。u 了解目前流行的潮流计算计算机算法，即高斯-赛德尔迭代法，牛顿拉夫逊法和P-Q分解法。u 结合要求编写程序进行潮流计算。u 结合MAT

13、LAB软件进行系统仿真。u 利用MATLAB软件，在(或给定) 某些运行参数的情况下,计算出系统中的全部参数,包括各母线电压的大小和相位、各发电机和负荷的功率及电流、以及各个变压器和线路等元件所通过的功率、电流和损耗。简单辐射网络的潮流计算。首先选定一个四节点或者五节点的区域电网，画出电路图。同一端的功率和电压，求另一端功率和电压。方法为从功率、电压端，齐头并进逐段求解功率和电压，如下列图所示求功率求电压求功率求电压求功率求电压U4S4计算起点1234始端电压、末端功率，求始端功率、末端电压以此居多；或末端电压、始端功率，求末端功率、始端电压。求解方法总结为“一来、二去共两步来逼近需求解的网络

14、功率和电压分布。一来即：设所有未知电压节点的电压为线路额定电压，从功率端开始逐段求功率，直到推得电压点得功率；二去即：从电压点开始，用推得的功率和电压点的电压，往回逐段向未知电压点求电压。在计算中，上述过程一般只需要做一次。但当一次“来、去完毕后，此电压与初始假设电压相差较大时，可再一次假设未知电压节点的电压值为刚刚计算得到的节点4电压值，继续进行“来、去计算，直到前后两次同一点的电压值相差不大。对节点导纳矩阵的形成和修改加以介绍，包括节点电压方程的意义，节点导纳矩阵的特点，以及在原有网络上增加节点或者添加删除支路对节点导纳矩阵的影响。对高斯-赛德尔迭代法计算的主要步骤和原理进行简单介绍，重点

15、介绍牛顿拉夫逊法直角坐标的计算过程，包括牛拉法的推导过程，原理和意义，然后给出修正方程式，形成雅克比矩阵，再介绍牛拉法的计算步骤。高斯-赛德尔迭代法计算的主要步骤1根据网络结构、参数，形成节点导纳矩阵；2迭代计算各节点电压。对PQ节点：a设某节点为平衡节点。b设各节点电压的初始值。c据初始值电压、 (i2、3、n)及的节点注入功率、进行第一次迭代，求、即有d第二次迭代。e第次迭代。f第k+1次迭代后，当所有节点电压都满足 (精度)时，说明迭代收敛，那么第k+1次的结果即为所求。3计算功率。牛顿拉夫逊法直角坐标的计算步骤1. 输入原始数据和信息：y、C、Pis、Qis、Uis、约束条件2. 形成

16、节点导纳矩阵YB=CTyC3. 设置各节点电压初值ei(0), fi(0)4. 将初始值代入(4-38)求不平衡量rPi(0), rQi(0), rUi2(0)5. 计算雅可比矩阵各元素Hij、Lij、Nij、Jij、Rij、Sij6. 解修正方程(4-37) ，求r ei(k), r fi(k)7. 求节点电压新值ei(k+1) =ei(k) -r ei(k), fi(k+1) = fi(k) -r fi(k)8. 判断是否收敛：Max| r fi(k) |, Max| r ei(k) |9. 重复迭代第4、5、6、7步，直到满足第8步的条件10. 求平衡节点的功率和PV节点的Qi及各支路的功率接下来对潮流计算的计算计算法必须满足的约束条件加以介绍，在论文的计算过程中，得出结果之后用约束条件进行校验，如果不能满足要求，那么要修改程序或者修改某些变量的给定值。2.2.3 用MATLAB