运用MATLAB软件进行潮流计算论文
- 格式:doc
- 大小:1.37 MB
- 文档页数:47
基于matlab的电力系统潮流计算仿真分析本文旨在介绍电力系统潮流计算仿真分析的背景和目的,并简要概述本文的主要内容和结构安排。
潮流计算是电力系统运行中的重要环节,通过计算电力系统中各节点的电压和功率分布情况,可以帮助分析系统的运行状态、调控能力以及潜在的问题。
随着电力系统的规模不断扩大和复杂性的增加,利用计算机进行潮流计算仿真分析已成为一种必要且有效的方法。
而matlab作为一种功能强大的科学计算软件,被广泛应用于电力系统的潮流计算仿真分析。
本研究的目的是基于matlab,开展电力系统潮流计算仿真分析,以探究系统运行状态、发现潜在的问题,并提出相应的优化方案。
通过仿真分析,可以评估系统的稳定性、安全性和可靠性,为电力系统运行与规划提供重要的参考依据。
本文主要包括以下内容:研究背景和意义:介绍电力系统潮流计算仿真分析的背景和其在电力系统运行中的重要性。
相关理论与方法:介绍电力系统潮流计算的基本理论和常用的计算方法,以及matlab在电力系统仿真中的应用。
模型构建与数据处理:详细阐述潮流计算仿真中的模型构建过程,以及对系统数据的处理和准备。
仿真结果与分析:展示仿真计算得到的结果,并进行相应的分析和讨论。
优化方案提出与评估:根据仿真结果,提出相应的优化方案,并进行评估和比较。
结论与展望:总结全文的研究内容和结论,并展望未来进一步的研究方向。
通过本文的研究和分析,我们将深入了解电力系统潮流计算仿真分析的原理和方法,为电力系统的优化和运行提供有效的技术支持。
本部分将介绍电力系统的组成,包括发电机组、输电网和配电网等,以及相关概念和术语,为后续的潮流计算仿真分析奠定基础。
潮流计算是电力系统中重要的分析方法,用于计算系统中各节点的电压幅值和相角,以及线路和设备的功率潮流分布。
潮流计算的基本原理是建立节点潮流方程和数学模型,通过求解这些方程来得到系统的潮流状态。
节点潮流方程节点潮流方程描述了电力系统中各节点的电压和功率之间的关系。
基于MATLAB进行潮流计算本文介绍了基于MATLAB软件的潮流计算方法。
电力系统潮流计算方法分为手算潮流和计算机潮流计算两类。
手算潮流主要适用于规模较小的辐射型电力潮流计算,而计算机潮流计算有两种途径:编程实现网络方程的迭代求解和借助电力系统分析仿真软件搭建系统模型完成潮流计算。
MATLAB具有强大的矩阵运算功能和电力系统仿真平台,可以为实现潮流计算提供更便捷的手段。
本文采用极坐标形式牛顿─拉夫逊法进行潮流计算,为其他形式的潮流计算提供借鉴。
Abstract: The power flow n method can be divided into two categories: hand n of tidal current and computer power flow XXX simplified equivalent circuits。
making it XXX: programming XXX ns。
or using power system XXX system model for power flow n。
MATLAB are has strong matrix ns and its power system XXX-Raphson method of power flow n in polar coordinates with MATLAB are。
and can serve as a reference for other forms of power flow n.1.电力系统中的牛顿法潮流计算是一种常用的电力系统分析方法。
该方法基于节点电压的相等条件和潮流方程的等式条件,通过迭代求解电压和相位的不平衡量,最终得到各节点的电压、相位和功率等参数。
2.牛顿法潮流计算的步骤包括输入系统原始数据、形成节点导纳矩阵、给定各节点电压初值、计算功率偏差向量、判断收敛条件、计算雅克比矩阵、解修正方程、计算节点电压和相位的修正值、迭代计算直至满足收敛条件、计算各节点功率等参数并输出计算结果。
基于MATLAB的电⼒系统潮流计算_毕业设计论⽂基于MATLAB的电⼒系统潮流计算摘要潮流计算是电⼒系统最基本、最常⽤的计算。
根据系统给定的运⾏条件、⽹络接线及元件参数,通过潮流计算可以确定各母线的电压(幅值及相⾓),各元件中流过的功率、整个系统的功率损耗等。
潮流计算是实现电⼒系统安全经济发供电的必要⼿段和重要⼯作环节。
因此潮流计算在电⼒系统的规划设计、⽣产运⾏、调度管理及科学研究中都有着⼴泛的应⽤。
本次设计的主要⽬的就是⾯向⼀般的电⼒⽹络,形成节点导纳矩阵,确定合适的算法,编写通⽤的计算程序,得到计算结果。
设计中主要介绍了⽜顿拉夫逊和PQ分解两种算法,PQ分解法虽然在结构上⽐⽜顿法更加简化,但是针对⼀般⽹络现代计算机在存储空间及计算速度上已经⼗分强⼤,鉴于对⽜顿法的熟悉与其算法的直观性,本次设计在编程时采⽤了⽜顿拉夫逊法的直⾓坐标形式。
解⽅程的过程利⽤Matlab的强⼤计算功能,编写M语⾔,合理设置变量,实现通⽤计算功能。
关键词: 电⼒系统,潮流计算,⽜顿—拉夫逊法,Matlab。
AbstractPower system load flow calculation is the most basic and commonly used calculations. Given according to the system operating conditions, the network connection and device parameters can be determined by power flow calculation of the bus voltage (magnitude and phase angle), the power flowing through the components, overall system power consumption and so on. Flow calculation is to achieve economic development of power system supply the necessary means and important part of the work. Therefore flow calculation in power system planning and design, production and operation, scheduling management, and scientific research have a wide range of applications.The main purpose of this design is for the general electricity network, the formation of the node admittance matrix, determine the appropriate method, the preparation of general-purpose computer program to get results. Introduces the design and the PQ decomposition Newton Raphson two algorithms, PQ decomposition although the structure is more streamlined than the Newton method, but for the general network of modern computer storage space and computing speed has been very strong, in view of the Newton Familiar with its intuitive algorithm, this design in programming using Newton Raphson polar form. The process of solving equations using matlab powerful computing capabilities, the preparation of M language, a reasonable set variables, to achieve general-purpose computing functions.Keywords: power system, power flow calculation, Newton - Raphson method, Matlab.⽬录摘要 (I)Abstract ................................................................................................................................................ II ⽬录.................................................................................................................................................... I II 1 引⾔ .. (1)1.1 潮流计算⽬的 (1)1.2 潮流计算意义 (1)1.3 潮流计算发展史 (1)1.4基于MATLAB 的电⼒系统潮流计算发展前景 (2)2简单电⼒系统潮流计算 (4)2.1简单辐射⽹络的潮流计算 (4)2.1.1简单⽀路的潮流分布和电压降落 (4)2.1.2 辐射型⽹络的⼿⼯潮流计算⽅法 (6)2.2 简单环⽹的潮流计算 (7)2.2.1两端电压相等 (7)2.2.2两端电压不相等 (8)3 复杂电⼒系统潮流计算的计算机算法 (10)3.1电⼒⽹络⽅程及等值电路 (10)3.2节点导纳矩阵形成及修改 (11)3.3节点的分类 (14)3.3.1 PQ节点 (14)3.3.2 PV节点 (14)3.3.3 平衡节点 (14)3.4潮流计算的约束条件 (15)3.5⽜顿-拉夫逊法(直⾓坐标) (15)3.5.1⽜顿-拉夫逊法的推导过程 (15)3.5.2潮流计算时的修正⽅程(直⾓坐标) (17) 3.5.3雅可⽐矩阵的特点: (19)3.5.4⽜顿-拉夫逊法计算步骤 (19)3.6 P-Q分解法潮流计算 (20)3.6.1 P-Q分解法潮流计算概述 (20)3.6.2 P-Q分解法的潮流计算步骤 (20)3.6.3 P-Q分解法的特点 (21)4 Matlab概述 (22)4.1Matlab简介 (22)4.2 Matlab中的变量 (22)4.3 Matlab编程 (23)4.3.1矩阵的输⼊ (23)4.3.2矩阵的运算 (24)4.3.3 MatLab的控制流 (24)5 ⽜顿法潮流计算程序设计及实例 (26)5.1⼿算 (26)5.2计算机算法的数据输⼊ (29)5.3潮流计算程序 (30)5.3 计算结果分析 (36)结论 (37)参考⽂献 (38)附录A 程序流程图 (39)附录B Matlab仿真 (40)致谢 (1)1 引⾔1.1 潮流计算⽬的电⼒系统潮流计算是研究电⼒系统稳态运⾏情况的⼀种基本电⽓计算。
1 引言1.1 本课题的目的和意义电力系统潮流计算是对复杂电力系统正常和故障条件下稳态运行状态的计算。
其目的是求取电力系统在给定运行方式下的节点电压和功率分布,用以检查系统各元件是否过负荷、各点电压是否满足要求、功率分布和分配是否合理以及功率损耗等,是电力系统计算分析中的一种最基本的计算[1]。
潮流计算是电力系统的各种计算的基础,同时它又是研究电力系统的一项重要分析功能,是进行故障计算,继电保护鉴定,安全分析的工具。
电力系统潮流计算是计算系统动态稳定和静态稳定的基础。
在电力系统规划设计和现有电力系统运行方式的研究中,都需要利用电力系统潮流计算来定量的比较供电方案或运行方式的合理性、可靠性和经济性[1]。
对于正在规划的电力系统,通过潮流计算,可以为选择电网供电方案和电气设备提供依据.潮流计算还可以为继电保护和自动装置整定计算、电力系统故障计算和稳定计算等提供原始数据。
潮流计算的目的在于:确定是电力系统的运行方式;检查系统中的各元件是否过压或过载;为电力系统继电保护的整定提供依据;为电力系统的稳定计算提供初值,为电力系统规划和经济运行提供分析的基础。
因此,电力系统潮流计算是电力系统中一项最基本的计算,既具有一定的独立性,又是研究其他问题的基础[1].1。
2 国内外发展现状利用电子计算机进行潮流计算从20世纪50年代中期就已经开始。
此后,潮流计算曾采用了各种不同的方法,这些方法的发展主要是围绕着对潮流计算的一些基本要求进行的。
对潮流计算的要求可以归纳为下面几点:(1)算法的可靠性或收敛性(2)计算速度和内存占用量(3)计算的方便性和灵活性电力系统潮流计算属于稳态分析范畴,不涉及系统元件的动态特性和过渡过程。
因此其数学模型不包含微分方程,是一组高阶非线性方程。
非线性代数方程组的解法离不开迭代,因此,潮流计算方法首先要求它是能可靠的收敛,并给出正确答案。
随着电力系统规模的不断扩大,潮流问题的方程式阶数越来越高,目前已达到几千阶甚至上万阶,对这样规模的方程式并不是采用任何数学方法都能保证给出正确答案的。
基于MATLAB的电力系统潮流计算【摘要】本文基于MATLAB,通过建立电力系统潮流计算的数学模型,实现了节点电压和支路功率的准确计算。
在正文部分中,详细介绍了数学模型的建立过程,节点电压和支路功率的计算方法,以及在MATLAB中的实现步骤。
通过对计算结果进行分析,发现了电力系统中可能存在的问题和优化方向。
在结论部分中,总结了本文的研究意义和实际应用价值,并展望了未来的改进方向。
该研究对电力系统的稳定运行和优化管理具有重要意义,为实际工程应用提供了有力支持。
通过本文的研究,可以更好地了解和应用MATLAB在电力系统潮流计算中的优势,推动电力系统领域的进一步发展。
【关键词】电力系统、潮流计算、MATLAB、数学模型、节点电压、支路功率、实现、结果分析、研究意义、改进展望、实际应用。
1. 引言1.1 概述电力系统潮流计算是电力系统分析中重要的一部分,通过对电力系统中节点间的电压、功率等参数进行计算,可以有效地评估系统的运行状况。
随着电力系统规模的不断扩大和电力负荷的增加,潮流计算的准确性和效率显得尤为重要。
基于MATLAB的电力系统潮流计算,可以帮助工程师更方便地进行系统分析和优化。
本文旨在探讨基于MATLAB的电力系统潮流计算方法,从数学模型建立开始,详细介绍节点电压计算和支路功率计算的过程,然后通过MATLAB编程实现这些计算。
将对计算结果进行分析,探讨其在电力系统优化中的应用前景。
通过本文的研究,将有助于深入理解电力系统潮流计算的原理和实现方法,为电力系统规划和运行提供更准确、高效的分析工具。
本文还将探讨MATLAB在电力系统领域的实际应用价值,为未来的研究和工程实践提供参考。
1.2 研究背景电力系统潮流计算是电力系统分析中十分重要的一个环节,它主要用于分析电力系统中各节点和支路上的电压、电流以及功率等参数。
通过潮流计算,可以帮助电力系统运行人员了解系统当前的负荷情况,优化系统运行,提高系统运行的效率和稳定性。
基于MATLAB进行潮流计算学生:王仕龙2011148213指导老师:李咸善摘要:电力系统潮流计算方法有两类,即手算潮流和计算机潮流计算。
手算潮流主要借助于形成简化的等值电路来实现,这种方法尤其适用于规模不大的辐射型电力潮流计算。
计算机潮流计算的实现有两种途径:其一是编程实现网络方程的迭代求解;其二是借助与电力系统分析仿真软件,搭建系统模型来完成潮流计算。
MATLAB具有强大的矩阵运算功能,同时其具有电力系统仿真平台也为直观地实现潮流计算提供了更便捷的手段[1]。
本文是基于MATLAB软件,采用极坐标形式牛顿─拉夫逊法进行潮流计算,为其他形式的潮流计算有借鉴的作用。
关键词:电力系统;计算机潮流计算;MATLAB ;牛顿─拉夫逊法Abstract:The power flow calculation method has two kinds,which are the hand calculation of tidal current and computer power flow calculation.Hand calculation tidal current is mainly realized by means of the formation of simplified equivalent circuit.This method is especially suitable for small scale radiation power flow calculation.There are two ways to realize the computer power flow calculation.The first one is through the programming iteration for solving network equation,the second one is with the help of analysis of power system simulation software to build the system model complete the power flow calculation.The software of MATLAB has strong matrix function,.At the same time,It’s power system simulation platform provides a more convenient means to realize power flow calculation intuitively[1].This paper is based on the software of MATLAB to calculate the power flow calculation by adopting the form of Newton-Raphson method of power flow calculation of polar coordinates.And it can be the role of reference of other forms of power flow calculation.Key words: power system computer; power flow calculation;MATLAB;Newton-Raphson1.计算原理电力系统潮流是指系统中所有运行参数的总体,包括各个母线电压的大小和相位,各个发电机和负荷的功率及电流,以及各个变压器和线路等元件所通过的功率,电流和其中的损耗。
电力系统潮流计算是电力系统分析和运行中的重要环节。
潮流计算主要用于确定电力系统中各个节点的电压、功率和潮流方向,以便进行功率平衡、电压稳定和线路负荷等方面的分析和评估。
MATLAB作为一种强大的数学建模和仿真工具,被广泛应用于电力系统潮流计算的研究和实际工程中。
本文将介绍MATLAB在电力系统潮流计算中的应用,包括算法原理、建模步骤和实例分析等内容。
一、潮流计算的基本原理潮流计算是指在给定电力网拓扑结构、负荷信息和发电机功率的情况下,通过迭代计算求解节点电压的复数值,以确定各节点的电压幅值和相角,进而计算各支路和各节点上的有功和无功功率。
潮流计算的基本原理是基于电力系统的潮流方程和节点功率平衡等基本理论,通过建立节点电压的复数方程组,利用迭代计算方法求解该方程组,从而得到节点的电压和功率信息。
二、MATLAB在潮流计算中的应用MATLAB作为一种功能强大的数学建模和仿真工具,具有丰富的数学计算函数和图形显示功能,适合于电力系统潮流计算的建模和仿真。
在MATLAB环境下,可以利用其矩阵运算、方程求解和数据可视化等功能,实现电力系统潮流计算的数学模型和算法的实现。
下面将介绍MATLAB在电力系统潮流计算中的具体应用步骤。
1. 建立电力系统潮流计算的数学模型在MATLAB环境下,首先需要建立电力系统潮流计算的数学模型,包括节点电压方程、支路潮流方程、节点功率平衡方程等。
利用MATLAB的矩阵运算和符号计算工具,可以将电力系统的节点和支路参数、负荷信息、发电机功率等数据表示为矩阵形式,建立电力系统潮流计算的数学模型。
2. 编写潮流计算的求解算法在建立电力系统潮流计算的数学模型后,需要编写潮流计算的求解算法。
在MATLAB环境下,可以利用其丰富的数学计算函数和优化工具,实现潮流计算的迭代求解算法,包括高斯-赛德尔迭代法、牛顿-拉夫逊迭代法等。
通过编写求解算法,可以实现电力系统潮流计算的数值求解过程。
3. 进行潮流计算的仿真实验在完成潮流计算的求解算法后,可以利用MATLAB进行潮流计算的仿真实验。
MATLAB在潮流计算演示教学中的应用MATLAB在潮流计算演示教学中的应用潮流计算是电力系统分析的重要内容之一,它可以求解电网中各节点电压、电流、功率、潮流等参数。
在电力系统教学中,潮流计算是一个重要的课程内容,也是学生理论知识与实践运用的桥梁。
因此,如何提高潮流计算教学的效果是学术界和工业界共同关注的问题。
MATLAB是一款功能强大的计算软件,广泛应用于工程、科学和教育领域,在潮流计算教学中也有重要的应用。
一、MATLAB在潮流计算仿真中的应用MATLAB可以进行各种电力系统仿真,包括潮流计算仿真。
它可以计算无功补偿、阻抗、励磁等各种因素对电网的影响,以及计算故障流、负荷流、电压降等等电网的各种参数。
通过使用MATLAB,学生可以模拟各种操作和状态,了解电网在不同负载状态下的响应,并研究电路参数对电网的影响。
二、MATLAB在潮流计算可视化中的应用MATLAB提供了强大的图形功能,可以用来可视化潮流计算。
学生可以绘制节点电压、相角、功率等参数的图示,以及各个节点之间电流流向的关系图。
在课堂上,教师可以通过这些可视化图像演示潮流计算的过程,提高学生的理解能力。
三、MATLAB在潮流计算分析中的应用MATLAB还可以用于对潮流计算结果进行分析,发现其中的规律和特点。
学生可以借助MATLAB进行统计分析、回归分析、协方差分析等方法,对潮流计算的结果进行深入探究。
这些分析结果有助于学生深入理解电网潮流的本质,提高其对电力系统的认知。
四、MATLAB在潮流计算算法优化中的应用MATLAB还可以用于潮流计算算法的优化,参与设计各种算法的控制语句和流程,包括线性规划、非线性规划、最小二乘法等多种算法。
学生可以在课程中了解各种算法的基本原理和具体操作流程,进而掌握MATLAB以及其他编程工具的使用方法,提高自主研究和创新能力。
综上所述,MATLAB在潮流计算教学中发挥了重要的作用,为学生提供了一个高效和灵活的学习平台。
摘要本文运用MATLAB软件进行潮流计算,对给定题目进行分析计算,再应用DDRTS软件,构建系统图进行仿真,最终得到合理的系统潮流。
潮流计算是电力系统最基本最常用的计算。
根据系统给定的运行条件,网络接线及元件参数,通过潮流计算可以确定各母线的电压幅值和相角,各元件流过的功率,整个系统的功率损耗。
潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。
因此,潮流计算在电力系统的规划计算,生产运行,调度管理及科学计算中都有着广泛的应用。
首先,画出系统的等效电路图,在计算出各元件参数的基础上,应用牛顿—拉夫逊Newton-Raphson法以及MATLAB软件进行计算对给定系统图进行了四种不同负荷下的潮流计算,经过调节均得到符合电压限制及功率限制的潮流分布。
其次,牛顿—拉夫逊Newton-Raphson法具有较好的收敛性,上述计算过程经过四到五次迭代后均能收敛。
根据运算结果,分析各支路损耗和系统总损耗。
最后,应用DDRTS软件,构建系统图,对给定负荷重新进行分析,潮流计算后的结果也能满足相应的参数要求。
关键词:牛顿-拉夫逊法MATLAB DDRTS 潮流计算目录1.摘要 (2)2.题目原始资料 (2)3.题目分析 (5)4.题目求解 (6)1)根据题意要求画出等值电路 (6)2)读程序画出拉夫逊法的流程图 (7)3)变电所负荷为题目所给数据进行求解 (8)4)编写程序并运行 (10)5)具体调压调损耗过程 (10)1.改变变压器变比调压 (10)2.改变发电机机端电压调压 (12)3.负荷按照一定比例变化的潮流计算分析 (15)4.轮流断开支路双回线中的一条的潮流计算 (19)5.仿真并比较 (26)6.设计心得 (28)7.参考文献 (29)一、 题目原始资料:1.系统图:两个发电厂分别通过变压器和输电线路与四个变电所相连。
变电所1变电所2母线2、发电厂资料:母线1和2为发电厂高压母线,发电厂一总装机容量为(300MW ),母线3为机压母线,机压母线上装机容量为(100MW),最大负荷和最小负荷分别为50MW和20MW;发电厂二总装机容量为(200MW )。
3、变电所资料:(一)变电所1、2、3、4低压母线的电压等级分别为:35KV 10KV 35KV 10KV(二)变电所的负荷分别为:50MW 50MW 40MW 70MW(三)每个变电所的功率因数均为cosφ=0.85;(四)变电所1和变电所3分别配有两台容量为75MVA的变压器,短路损耗414KW,短路电压(%)=16.7;变电所2和变电所4分别配有两台容量为63MVA的变压器,短路损耗为245KW,短路电压(%)=10.5;4、输电线路资料:发电厂和变电所之间的输电线路的电压等级及长度标于图中,单位长度的电阻为Ω.0,单位长度的电抗为Ω0.402,单位长度的电纳为S-617*2.78。
10二、课程设计容及要求(一)课设容:1.对给定的网络查找潮流计算所需的各元件等值参数,画出等值电路图。
2.输入各支路数据,各节点数据利用给定的程序进行在变电所在某一负荷情况下的潮流计算,并对计算结果进行分析。
3.跟随变电所负荷按一定比例发生变化,进行潮流计算分析。
1)4个变电所的负荷同时以2%的比例增大;2)4个变电所的负荷同时以2%的比例下降3)1和4号变电所的负荷同时以2%的比例下降,而2和3号变电所的负荷同时以2%的比例上升;4.在不同的负荷情况下,分析潮流计算的结果,如果各母线电压不满足要求,进行电压的调整。
(变电所低压母线电压10KV要求调整围在9.5-10.5之间;电压35KV要求调整围在35-36之间)5.轮流断开支路双回线中的一条,分析潮流的分布。
(几条支路断几次)6.利用DDRTS软件,进行绘制系统图进行上述各种情况潮流的分析,并进行结果的比较。
7.最终形成课程设计成品说明书。
(二) 课设要求:1.在读懂程序的基础上画出潮流计算基本流程图2.通过输入数据,进行潮流计算输出结果3.对不同的负荷变化,分析潮流分布,写出分析说明。
4.对不同的负荷变化,进行潮流的调节控制,并说明调节控制的方法,并列表表示调节控制的参数变化。
5.打印利用DDRTS 进行潮流分析绘制的系统图,以及潮流分布图。
三、 题目分析1.节点设置及分类根据系统图可知此系统为两端供电网路,将母线1,2设为节点1,2,将变电所1、2、3、4的高低压侧分别设为节点3、4、5、6、7、8、9、10。
并且,将节点1设为平衡节点,将节点2设为PV 节点,其余节点设为PQ 节点。
2.参数求取根据题目原始资料,计算发电厂、变压器及线路的参数。
(1)运用下列公式计算变压器参数:采用变压器参数为折算至高压侧的数值,其变比K<1,其中,SUP R NNkT221000=SU U XNNKT1002%=X R ZTT T+=计算变压器分接头变比变压器有5个抽头,电压调节围为N U ±2*2.5%, N U 对应的分接头开始时设变压器高压侧接主接头。
5.382315.38%)51(2201=+⨯=k 1123111%)51(2201=+⨯=k5.385.2255.38%)5.21(2202=+⨯=k 115.22511%)5.21(2202=+⨯=k5.382205.3812203=⨯=k 112201112203=⨯=k5.385.2145.38%)5.21(2204=-⨯=k 115.21411%)5.21(2204=-⨯=k5.382095.38%)51(2205=-⨯=k 1120911%)51(2205=-⨯=k(2)计算线路参数L jx r jX R Z )(+=+= (3)变电所负荷分别为:变电所1 L S =50+j30.98722 变电所2 L S =50+j30.98722 变电所3 L S =40+j24.78977 变电所4 L S =70+j43.38210 将参数整理如下表:利用牛顿拉夫逊法进行求解,用MATLAB软件编程,可以求解系统潮流分布根据题目的不同要求对参数进行调整,通过调节变压器变比和发电厂的电压,求解出合理的潮流分布,最后用DDRTS进行潮流分析,将两者进行比较。
四、题目求解(一)根据题意要求画出等值电路(二)读程序画出牛顿拉夫逊法的流程图(三)变电所负荷为题目所给数据进行求解B1矩阵:1、支路首端号;2、末端号;3、支路阻抗;4、线路对地电纳(或变压器导纳);5、支路的变比;6、支路首端处于K侧为1,1侧为0;7、线路/变压器标识(0/1)变压器参数当支路首端处于K侧标识为1时归算至末端侧,0归算至首端侧B2矩阵:1、该节点发电机功率;2、该节点负荷功率;3、PQ节点电压初始值,或平衡节点及PV节点电压的给定值4、节点所接无功补偿并联电容(感)的电纳5、节点分类标号:1为平衡节点(应为1号节点);2为PQ节点;3为PV节点;B1=[ 1 3 8.5+j20.1 j0.000556 1 01 5 4.25+j10.05 j0.000278 1 01 7 7.65+j18.09 j0.0005004 1 02 9 8.5+j20.1 j0.000556 1 03 4 0.05455+j1.65024 0 220/38.5 15 6 0.0037345+j0.10083 0 220/11 17 8 0.05455+j1.65024 0 220/38.5 17 9 6.8+j16.08 j0.0004448 1 09 10 0.0037345+j0.10083 0 220/11 1]B2=[ 0 0 231 0 1200 0 231 0 30 0 220 0 20 50+j30.98722 35 0 20 0 220 0 20 50+j30.98722 10 0 20 0 220 0 20 40+j24.78977 35 0 20 0 220 0 20 70+j43.38210 10 0 2]由于各节点电压不应超过各自的允许围,由题目要求知,变压器低压侧电压有名值变化围如下表:同时保证发电厂发出的有功功率要大于变电所有功功率,并且两者差值越小表明线路有功功率损耗越小,越满足实际要求。
(四)编写程序并运行由于程序较复杂,故将其放入附录中,具体程序和结果见附录一由附录一的结果可以观察到各个节点电压有名值分别为:由此观察到节点4、6、8、10的电压都较正常围偏高,因此调节变压器分接头和发电厂电压,此系统进行多次调整,最终得到合理结果。
(五)具体调压调损耗过程1.改变变压器变比调压由于调节的次数很多,过程较为复杂,只选取其中有代表性的两组数据如下所述。
调节方法:调节结果:由电压和功率损耗的有名值可以观察到,最后一次的电压和有功损耗与题目要求较接近,但仍超出要求的围。
第二次调变比的具体的潮流分布如下:各节点的功率S为(节点号从小到大排列):Columns 1 through 50.2224 + 0.9261i 2.0000 - 0.3935i -0.0000 + 0.0000i -0.5000 - 0.3099i 0.0000 + 0.0000iColumns 6 through 10-0.5000 - 0.3099i -0.0000 - 0.0000i -0.4000 - 0.2479i 0.0000 + 0.0000i -0.7000 - 0.4338i得到的图像如下:2.改变发电机机端电压调压由于调节变压器分接头达到极限仍未达到题目要求的围,所以接下来调节发电机机端电压。
此系统经过多次调节后达到要求,选取三个调节过程,对应的调节方法和结果见下表:调节方法:调节方法电厂1电压厂2电压分接头1 分接头2 分接头3 分接头438.5/2311/231 38.5/231 11/231 未调整230 230138.5/23第一次229 22911/231 38.5/231 11/231138.5/23第二次227 22711/231 38.5/231 11/2311调节结果:由电压和功率损耗的有名值可以观察到,最后一次的电压符合题目要求,可认为是合理的。
具体的潮流分布如下表:各节点的功率S为(节点号从小到大排列):Columns 1 through 50.2281 + 0.9732i 2.0000 - 0.3651i -0.0000 - 0.0000i -0.5000 - 0.3099i 0.0000 + 0.0000iColumns 6 through 10-0.5000 - 0.3099i 0.0000 + 0.0000i -0.4000 - 0.2479i -0.0000 + 0.0000i -0.7000 - 0.4338i得到的图像如下:与DDRTS比较:潮流计算方式节点4 节点6 节点8 节点10 DDRTS 35.2307 10.274 35.6356 10.0043 MATLAB 35.2401 10.2766 35.6513 10.00173.负荷按照一定比例变化的潮流计算的分析⑴4个变电所的负荷同时以2%比例增大;为了使电压符合题目要求,因此将K4变比由11/231调为11/225.5各节点的电压模值大小V为(节点号从小到大排列):Columns 1 through 6225.0000 225.0000 220.9222 35.1842 222.7473 10.2668Columns 7 through 10221.3415 35.6094 220.1434 10.2557各节点的电压相角sida为(节点号从小到大排列):Columns 1 through 60 9.9672 -0.9616 -4.6100 -0.4506 -3.0950 Columns 7 through 102.2891 -0.5886 5.1016 1.4397各节点的功率S为(节点号从小到大排列):1.0e+002 * Columns 1 through 30.2696 + 0.9911i 2.0000 - 0.3541i -0.0000 + 0.0000i Columns 4 through 6-0.5100 - 0.3161i 0.0000 + 0.0000i -0.5100 - 0.3161i Columns 7 through 9-0.0000 + 0.0000i -0.4080 - 0.2529i 0.0000 + 0.0000i Column 10-0.7140 - 0.4425i节点电压000 000 222 842 473 668 425 094 434 557各节点相角0 9.9672-0.9616-4.6100-0.4506-3.09502.2891-0.58865.10161.4397各节点功率2.2502.2160 +0.3894i2.2089 -0.0371i0.3507 -0.0283i2.2274 -0.0175i0.1025 -0.0055i2.2116 +0.0884i0.3561 -0.0037i2.1927 +0.1958i0.1025 +0.0026i与DDRTS比较:潮流计算方式节点4 节点6 节点8 节点10 DDRTS 35.2156 10.2718 35.6194 10.2548 MATLAB 35.1842 10.2668 35.6094 10.2557⑵四个变电所的负荷同时以2%的比例下降;节点1 节点2 节点3 节点4 节点5 节点6节点7 节点8节点9 节点10各节点电压225.0000225.0000221.157235.2956222.855910.2863221.525135.6951220.368810.0186各节点相0 10.1203-0.9302-4.4206-0.4347-2.96932.3717-0.38415.23561.5499各节点功率0.1866 +0.9540i2.0000 -0.3785i0.0000 +0.0000i-0.4900 -0.3037i0.0000 +0.0000i-0.4900 -0.3037i0.0000 -0.0000i-0.3920 -0.2429i0.0000 +0.0000i-0.6860 -0.4251i潮流计算方式节点4 节点6 节点8 节点10 DDRTS 35.3107 10.2884 35.7508 10.0843 MATLAB 35.2956 10.2863 35.6951 10.0186⑶变电所1、4负荷同时以2%的比例下降,变电所2、3负荷同时以2%的比例上升:节点1 节点2 节点3 节点4节点5 节点6节点7 节点8节点9 节点10各节点电225.0000225.0000221.157235.2956222.747310.2668221.422435.6234220.311610.0157各节点相角0 10.0861-0.9302-4.4206-4.4206-3.09502.3451-0.53045.20621.5185各节点功率0.2226 +0.9755i2.0000 -0.3723i0.0000 +0.0000i-0.4900 -0.3037i0.0000 +0.0000i-0.5100 -0.3161i0 -0.4080 -0.2529i-0.0000 +0.0000i-0.6860 -0.4251i潮流计算方式节点4 节点6 节点8 节点10 DDRTS 35.3107 10.2718 35.6789 10.0815 MATLAB 35.2956 10.2668 35.6234 10.01574.轮流断开支路双回线中的一条的潮流计算⑴断开①③支路后,参数变化及其调节方法:为了使电压符合题目要求,因此将K1变比由38.5/231调为38.5/225.5 支路矩阵B1第一行发生变化:B1=[ 1 3 17+j40.2 j0.000278 1 0]各节点的电压模值大小V为(节点号从小到大排列):Columns 1 through 6225.0000 225.0000 215.4963 35.1896 222.8017 10.2766Columns 7 through 10221.4275 35.6513 220.2453 10.0017各节点的电压相角sida为(节点号从小到大排列):Columns 1 through 60 10.0441 -1.7973 -5.3763 -0.4427 -3.0321Columns 7 through 102.3309 -0.4859 5.1697 1.4003各节点的功率S为(节点号从小到大排列):1.0e+002 *Columns 1 through 30.2352 + 1.1318i 2.0000 - 0.3651i 0Columns 4 through 6-0.5000 - 0.3099i -0.0000 - 0.0000i -0.5000 - 0.3099iColumns 7 through 90.0000 + 0.0000i -0.4000 - 0.2479i -0.0000 + 0.0000iColumn 10-0.7000 - 0.4338i节点1 节点2 节点3 节点4 节点5 节点6节点7 节点8节点9 节点10各节点电压225.0000225.0000215.496335.1896222.801710.2766221.427535.6513220.245310.0017各节点相角0 10.0441-1.7973-5.3763-0.4427-3.03212.3309-0.48595.16971.4003各节点功率0.2352 +1.1318i2.0000 -0.3651i0.0000 +0.0000i-0.5000 -0.3099i-0.0000 -0.0000i-0.5000 -0.3099i0.0000 +0.0000i-0.4000 -0.2479i-0.0000 +0.0000i-0.7000 -0.4338i与DDRTS比较:潮流计算方式节点4 节点6 节点8 节点10 DDRTS 35.261 10.2864 35.692 10.0059 MATLAB 35.1896 10.2766 35.6513 10.0017⑵断开①⑤支路后,参数变化及其调节方法:为了使电压符合题目要求,因此将K2变比由11/231调为11/225.5 支路矩阵B第二行发生变化:B1=[ 1 5 8.5+j20.1 j0.000139 1 0]节点1 节点2 节点3 节点4 节点5 节点6节点7 节点8节点9 节点10各节点电压225.0000225.0000221.040035.2401220.226010.4144221.427535.6513220.245310.0017各节点相角0 10.0441-0.9459-4.5151-0.8618-3.38532.3309-0.48595.16971.4003各0.231 2.000-0.00-0.50-0.00-0.500.000-0.40-0.00-0.70节点功率4 +1.0506i0 -0.3651i00 -0.0000i00 -0.3099i00 +0.0000i00 -0.3099i0 +0.0000i00 -0.2479i00 +0.0000i00 -0.4338i潮流计算方式节点4 节点6 节点8 节点10 DDRTS 35.2966 10.204 35.692 10.0059 MATLAB 35.2401 10.4144 35.6513 10.0017⑶断开①⑦支路后,参数变化及其调节方法:为了使电压符合题目要求,因此将K4变比由11/231调为11/225.5支路矩阵B第三行发生变化:B1=[ 1 7 15.3+j36.18 j0.0002502 1 0]节点1 节点2 节点3 节点4节点5 节点6节点7 节点8节点9 节点10各节225.0000225.0000221.040035.2401222.801710.2766218.934035.2193218.870710.2007点电压各节点相角0 12.0967-0.9459-4.5151-0.4427-3.03214.57141.68747.33973.7093各节点功率0.2411 +0.9845i2.0000 -0.2160i0.0000 -0.0000i-0.5000 -0.3099i-0.0000 -0.0000i-0.5000 -0.3099i0.0000 -0.0000i-0.4000 -0.2479i-0.0000 +0.0000i-0.7000 -0.4338i潮流计算方式节点4 节点6 节点8 节点10 DDRTS 35.2958 10.2072 35.2153 10.1955 MATLAB 35.2401 10.2766 35.2193 10.2007(4)断开⑦⑨支路后,参数变化及其调节方法:为了使电压符合题目要求,因此将K4变比由11/231调为11/225.5 支路矩阵B第八行发生变化:B1=[7 9 13.6+j32.16 j0.0002224 1 0]节点1 节点2 节点3 节点4 节点5 节点6节点7 节点8节点9 节点10各节点电压225.0000225.0000221.040035.2401222.801610.2766220.972835.5726219.121210.2136各节点相角0 12.7172-0.9459-4.5151-0.4427-3.03212.324-0.50487.93874.3171各节点功率0.2531 +1.0190i2.0000 -0.2432i-0.0000 -0.0000i-0.5000 -0.3099i0.0000 +0.0000i-0.5000 -0.3099i-0.0000 -0.0000i-0.4000 -0.2479i-0.0000 +0.0000i-0.7000 -0.4338i(5)断开②⑨支路后,参数变化及其调节方法:为了使电压符合题目要求,因此将K4变比由11/231调为11/225.5 支路矩阵B第四行发生变化:B1=[2 9 17+j40.2 j0.000278 1 0]与DDRTS比较:潮流计算方式节点4 节点6 节点8 节点10 DDRTS 35.2966 10.2076 35.3482 10.0644 MATLAB 35.2401 10.2766 35.302 10.0794 五、仿真并比较应用DDRTS软件,以系统变电所为正常负荷情况构建系统图,如下:进行潮流计算,得到潮流图如下:将DDRTS仿真得到的潮流图中各个点电压的有名值与MATLAB计算得到的各节点电压进行比较:由上述比较可观察到: 应用仿真和MATLAB计算的到的结果大致相同,可以认为此潮流计算是正确的。