电力系统潮流分析
电力系统最优潮流分析

电力系统最优潮流分析电力系统是现代社会中最重要的系统工程之一,为社会生产和人民生活提供了绝大部分能量。
电能的生产需要耗费大量的燃料,而目前电能在输送、分配和消费过程中存在着大量的损耗。
因此如何采取适当措施节约能源,提高整个电力系统的运行效率,优化系统的运行方式,是国内外许多学者一直关注与研究的热点。
电力系统的最优化运行是指在确保电力系统安全运行、满足用户用电需求的前提下,如何通过调度系统中各发电机组或发电厂的运行,从而使系统发电所需的总费用或所消耗的总燃料达到最小的运筹决策问题。
数学上可将此问题描述为非线性规划或混合非线性规划问题。
最优潮流问题是指在满足必须的系统运行和安全约束条件下,通过调整系统中可利用控制手段实现预定目标最优的系统稳定运行状态。
同经典的经济调度法相比,最优潮流具有全面规划、统筹考虑等优点,它可将安全运行和最优经济运行等问题进行综合考虑,通过统一的数学模型来描述,从而将电力系统对经济性、安全性以及电能质量等方面的要求统一起来。
最优潮流问题的提出把电力系统的最优运行理论提高到一个新的高度,受到了国内外学者高度重视。
最优潮流已在电力系统中的安全运行、电网规划、经济调度、阻塞管理、可靠性分析以及能量管理系统等方面得到了广泛应用,成为了电力系统网络运行分析和优化中不可或缺的工具。
一、最优潮流问题研究的意义最优潮流可将电力系统可靠性与电能质量量化成相应的经济指标,并最终达到优化资源配置、降低成本、提高服务质量的目的。
因此最优潮流研究具有传统潮流计算无法比拟的意义,主要体现在以下两个方面。
一方面,通过最优潮流计算可指导系统调度员的操作,保证系统在经济、安全、可靠的状态下运行。
具体表现为:第一,当所求问题以目标函数、控制变量和约束条件的形式固定下来后,就一定可以求出唯一最优解,并且该结果不受人为因素的影响。
第二,最优潮流的寻优过程可以自动识别界约束,在解逐渐趋于最优的过程中可得到网络传输瓶颈信息,从而可以指导电网扩容与规划。
潮流分析

与电路原理的分析方法有何区别? 与电路原理的分析方法有何区别?
模型和方法的特殊性: 模型和方法的特殊性:潮流模型和方法研究
条件已经变了:复电流复功率,非线性、 条件已经变了:复电流复功率,非线性、复杂性 电力系统特性:功率平衡(时时处处) 电力系统特性:功率平衡(时时处处) 数学基础: 数学基础:非线性方程组的数值解
把元件抽象成支路,研究支路之间的联结关系。 把元件抽象成支路,研究支路之间的联结关系。 并联) (串、并联) 元件之间的联结关系是网络拓扑, 元件之间的联结关系是网络拓扑,它与元件特性无 关,仅取决于运行方式。 仅取决于运行方式。
网络元件特性约束
元件特性约束与元件联结关系无关
基本原理
13
网络方程
网络矩阵是描述网络的数学工具。 网络矩阵是描述网络的数学工具。
以姆欧继电器为例
区外/反向接地故障由于过渡电阻的影响, 区外 反向接地故障由于过渡电阻的影响,引起故障 反向接地故障由于过渡电阻的影响 电压相量落入动作区= 误动 电压相量落入动作区=>误动 保护安装在受端, 保护安装在受端,正向区外经小电阻接地故障 保护安装在受端,背后母线经小电阻接地故障 保护安装在受端, 经小电阻接地两相故障, 经小电阻接地两相故障,超前相整定点外超越
21
潮流解的工程约束
22
PQ和PV节点互相转化 和 节点互相转化
结点的无功功率越界, 如PU结点的无功功率越界,此时 结点转化 结点的无功功率越界 此时PU结点转化 结点。 为PQ结点。 结点
对于PQ结点,若计算结果Ui越界,需修正该结 对于 结点,若计算结果 越界, 结点 点的无功功率使结点电压维持在允许范围内, 点的无功功率使结点电压维持在允许范围内, 实质上即使该结点转化为PU结点。 实质上即使该结点转化为 结点。 结点
九节点电力系统分析潮流程序设计

九节点电力系统分析潮流程序设计下载提示:该文档是本店铺精心编制而成的,希望大家下载后,能够帮助大家解决实际问题。
文档下载后可定制修改,请根据实际需要进行调整和使用,谢谢!本店铺为大家提供各种类型的实用资料,如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等,想了解不同资料格式和写法,敬请关注!Download tips: This document is carefully compiled by this editor. I hope that after you download it, it can help you solve practical problems. The document can be customized and modified after downloading, please adjust and use it according to actual needs, thank you! In addition, this shop provides you with various types of practical materials, such as educational essays, diary appreciation, sentence excerpts, ancient poems, classic articles, topic composition, work summary, word parsing, copy excerpts, other materials and so on, want to know different data formats and writing methods, please pay attention!1.引言随着电力系统规模的不断扩大和电力市场的逐渐开放,电力系统潮流分析成为了电力系统规划、运行和市场交易中不可或缺的一部分。
简单电力系统潮流分析

电力系统分析(刘学军主编)机械工业出版社
第3章简单电力系统的潮流分析
(2) 变压器并联支路的损耗
2 2 SYT P j Q Y U U G jB U G jB U YT YT T 1 1 T T 1 T T 1
(3-11)
PY 2 QY 2
1 2 GU 2 (3-12) 2 1 2 BU 2 2
线路首端;
1 Y 1 2 2 SY 1 U1 U1 GU1 j BU1 2 2 2 PY 1 j QY 1 1 PY 1 GU12 2 1 QY 1 BU12 2
图3-2电力线路的π形等值电路
电力系统分析(刘学军主编)机械工业出版社
第3章简单电力系统的潮流分析
1.电力线路阻抗中功率损耗 (1)串联支路的单相功率损耗 设已知线路末端流出功率 S ,末端电压 U 2 ,则串联支 2 路的单相功率损耗为:
S2 2 P22 Q2 S R jX Z 2 2 U2 U2 2 2 P22 Q2 P22 Q2 R j X 2 2 U2 U2 P2 j Q2 2 2 P2 Q2 P2 R 2 U2 2 P22 Q2 Q2 X 2 U2
第3章简单电力系统的潮流分析
电力系统分析教材配套课件
第3章简单电力系统分析
电力系统分析(刘学军主编)机械工业出版社
第3章简单电力系统的潮流分析
3.1电力网络的电压降落和功率损耗 3.2输电线路的运行特性 3.3开式电网的潮流分析
电力系统潮流计算方法分析

电力系统潮流计算方法分析1.黎曼法是最简单和最直接的计算方法。
该方法直接利用电力系统的基本方程式,即功率平衡方程式和节点电压方程式来计算潮流分布。
然而,黎曼法需要利用复杂的矩阵方程式来解决系统中节点电压的计算,计算量大且计算速度较慢,对大型复杂系统不适用。
2.高斯-赛德尔法是一种迭代法,将电网中的节电清设置为未知数,并采用全局迭代求解。
该方法通过迭代计算不断逼近潮流分布,直到满足系统中所有节点的电压和功率平衡方程为止。
高斯-赛德尔法具有迭代次数多、耗时较长的缺点,但计算稳定可靠,对于小型系统具有较好的适用性。
3.牛顿-拉夫逊法是一种基于牛顿迭代思想的高效潮流计算方法。
该方法通过利用电力系统中的雅可比矩阵,将潮流计算问题转化为解非线性方程组的问题。
牛顿-拉夫逊法的迭代速度和稳定性较高,适用于大型复杂系统的潮流计算。
综上所述,电力系统潮流计算方法可以选择黎曼法、高斯-赛德尔法和牛顿-拉夫逊法等不同的算法进行计算。
选择合适的计算方法应根据系统的规模、复杂度以及计算时间要求来综合考虑。
实际应用中,通常会根据具体情况采用不同的方法进行潮流计算,以获得准确和高效的结果。
同时,随着电力系统的发展和智能化技术的应用,也出现了一些基于机器学习和深度学习的潮流计算方法。
这些方法利用大数据和智能算法,通过学习和分析系统历史数据,能够更好地预测和计算系统潮流分布,提高计算效率和准确性。
这些方法在未来的电力系统潮流计算中具有潜力和广阔的应用前景。
总结起来,电力系统潮流计算是电力系统分析和规划的重要工作,不同的计算方法有不同的优劣势,合理选择计算方法对于准确评估系统稳定性和可靠性至关重要。
随着技术的进步和应用的发展,电力系统潮流计算方法也在不断演化和改进,以满足电力系统智能化和可持续发展的需求。
最新电力系统分析潮流计算实验报告

电力系统分析潮流计算实验报告姓名:XXXXXX 学号:XXXXXXXXXX 班级:XXXXXXXX一、实验目的掌握潮流计算计算机算法的方法,熟悉MATLAB的程序调试方法。
二、实验准备根据课程内容,熟悉MATLAB软件的使用方法,自行学习MATLAB程序的基础语法,并根据所学知识编写潮流计算牛顿拉夫逊法(或PQ分解法) 的计算程序,用相应的算例在MATLAB上进行计算、调试和验证。
三、实验要求每人一组,在实验课时内,调试和修改运行程序,用算例计算输出潮流结果。
四、程序流程五、实验程序%本程序的功能是用牛拉法进行潮流计算%原理介绍详见鞠平著《电气工程》%默认数据为鞠平著《电气工程》例8.4所示数据%B1是支路参数矩阵%第一列和第二列是节点编号。
节点编号由小到大编写%对于含有变压器的支路,第一列为低压侧节点编号,第二列为高压侧节点编号%第三列为支路的串列阻抗参数,含变压器支路此值为变压器短路电抗%第四列为支路的对地导纳参数,含变压器支路此值不代入计算%第五烈为含变压器支路的变压器的变比,变压器非标准电压比%第六列为变压器是否是否含有变压器的参数,其中“1”为含有变压器,“0”为不含有变压器%B2为节点参数矩阵%第一列为节点注入发电功率参数%第二列为节点负荷功率参数%第三列为节点电压参数%第四列%第五列%第六列为节点类型参数,“1”为平衡节点,“2”为PQ节点,“3”为PV节点参数%X为节点号和对地参数矩阵%第一列为节点编号%第二列为节点对地参数%默认算例% n=4;% n1=4;% isb=4;% pr=0.00001;% B1=[1 2 0.1667i 0 0.8864 1;1 3 0.1302+0.2479i 0.0258i 1 0;1 4 0.1736+0.3306i 0.0344i 1 0;3 4 0.2603+0.4959i 0.0518i 1 0];% B2=[0 0 1 0 0 2;0 -0.5-0.3i 1 0 0 2;0.2 0 1.05 0 0 3;0 -0.15-0.1i 1.05 0 0 1];% X=[1 0;2 0.05i;3 0;4 0];clear;clc;num=input('是否采用默认数据?(1-默认数据;2-手动输入)');if num==1n=4;n1=4;isb=4;pr=0.00001;B1=[1 2 0.1667i 0 0.8864 1;1 3 0.1302+0.2479i 0.0258i 1 0;1 4 0.1736+0.3306i 0.0344i 1 0;3 4 0.2603+0.4959i 0.0518i 1 0];B2=[0 0 1 0 0 2;0 -0.5-0.3i 1 0 0 2;0.2 0 1.05 0 0 3;0 -0.15-0.1i 1.05 0 0 1];X=[1 0;2 0.05i;3 0;4 0];elsen=input('请输入节点数:n=');n1=input('请输入支路数:n1=');isb=input('请输入平衡节点号:isb=');pr=input('请输入误差精度:pr=');B1=input('请输入支路参数:B1=');B2=input('请输入节点参数:B2=');X=input('节点号和对地参数:X=');endTimes=1; %迭代次数%创建节点导纳矩阵Y=zeros(n);for i=1:n1if B1(i,6)==0 %不含变压器的支路p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1/B1(i,3);Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4);Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4);else %含有变压器的支路p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-B1(i,5)/B1(i,3);Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+B1(i,5)/B1(i,3)+(1-B1(i,5))/B1(i,3);Y(q,q)=Y(q,q)+B1(i,5)/B1(i,3)+(B1(i,5)*(B1(i,5)-1))/B1(i,3);endendfor i=1:n1Y(i,i)=Y(i,i)+X(i,2); %计及补偿电容电纳enddisp('导纳矩阵为:');disp(Y); %显示导纳矩阵%初始化OrgS、DetaSOrgS=zeros(2*n-2,1);DetaS=zeros(2*n-2,1);%创建OrgS,用于存储初始功率参数h=0;j=0;for i=1:n %对PQ节点的处理if i~=isb&B2(i,6)==2 %不是平衡点&是PQ点h=h+1;for j=1:n%公式8-74%Pi=ei*(Gij*ej-Bij*fj)+fi*(Gij*fj+Bij*ej)%Qi=fi*(Gij*ej-Bij*fj)-ei*(Gij*fj+Bij*ej)OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real (Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j ,3))+imag(Y(i,j))*real(B2(j,3)));endendendfor i=1:n %对PV节点的处理,注意这时不可再将h初始化为0if i~=isb&B2(i,6)==3 %不是平衡点&是PV点h=h+1;for j=1:n%公式8-75-a%Pi=ei*(Gij*ej-Bij*fj)+fi*(Gij*fj+Bij*ej)%Qi=fi*(Gij*ej-Bij*fj)-ei*(Gij*fj+Bij*ej)OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real (Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j ))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendend%创建PVU 用于存储PV节点的初始电压PVU=zeros(n-h-1,1);t=0;for i=1:nif B2(i,6)==3t=t+1;PVU(t,1)=B2(i,3);endend%创建DetaS,用于存储有功功率、无功功率和电压幅值的不平衡量h=0;for i=1:n %对PQ节点的处理if i~=isb&B2(i,6)==2h=h+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1); %delPiDetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1); %delQiendendt=0;for i=1:n %对PV节点的处理,注意这时不可再将h初始化为0if i~=isb&B2(i,6)==3h=h+1;t=t+1;DetaS(2*h-1,1)=real(B2(i,1))-OrgS(2*h-1,1); %delPiDetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2; %delUi endend% DetaS%创建I,用于存储节点电流参数i=zeros(n-1,1);h=0;for i=1:nif i~=isbh=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));%conj求共轭endend%创建Jacbi(雅可比矩阵)Jacbi=zeros(2*n-2);h=0;k=0;for i=1:n %对PQ节点的处理if B2(i,6)==2h=h+1;for j=1:nif j~=isbk=k+1;if i==j %对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));else %非对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);endif k==(n-1) %将用于内循环的指针置于初始值,以确保雅可比矩阵换行k=0;endendendendendk=0;for i=1:n %对PV节点的处理if B2(i,6)==3h=h+1;for j=1:nif j~=isbk=k+1;if i==j %对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=2*imag(B2(i,3));Jacbi(2*h,2*k)=2*real(B2(i,3));else %非对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;endif k==(n-1) %将用于内循环的指针置于初始值,以确保雅可比矩阵换行k=0;endendendendenddisp('初始雅可比矩阵为:');disp(Jacbi);%求解修正方程,获取节点电压的不平衡量DetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS; %inv矩阵求逆% DetaU%修正节点电压j=0;for i=1:n %对PQ节点处理if B2(i,6)==2j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendfor i=1:n %对PV节点的处理if B2(i,6)==3j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endend% B2%开始循环**********************************************************************while abs(max(DetaU))>prOrgS=zeros(2*n-2,1);h=0;j=0;for i=1:nif i~=isb&B2(i,6)==2h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real (Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j ))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendendfor i=1:nif i~=isb&B2(i,6)==3h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real (Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j ))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendend% OrgS%创建DetaSh=0;for i=1:nif i~=isb&B2(i,6)==2h=h+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);endendt=0;for i=1:nif i~=isb&B2(i,6)==3h=h+1;t=t+1;% DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h-1,1)=real(B2(i,1))-OrgS(2*h-1,1);DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2;endend% DetaS%创建Ii=zeros(n-1,1);h=0;for i=1:nif i~=isbh=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));endend% I%创建JacbiJacbi=zeros(2*n-2);h=0;k=0;for i=1:nif B2(i,6)==2h=h+1;for j=1:nif j~=isbk=k+1;if i==jJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));elseJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);endif k==(n-1)k=0;endendendendendk=0;for i=1:nif B2(i,6)==3h=h+1;for j=1:nif j~=isbk=k+1;if i==jJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=2*imag(B2(i,3));Jacbi(2*h,2*k)=2*real(B2(i,3));elseJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;endif k==(n-1)k=0;endendendend% JacbiDetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS;% DetaU%修正节点电压j=0;for i=1:nif B2(i,6)==2j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendfor i=1:nif B2(i,6)==3j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endend% B2Times=Times+1; %迭代次数加1enddisp('迭代次数为:');disp(Times);disp('收敛时电压修正量为::');disp(DetaU);for k=1:nE(k)=B2(k,3);e(k)=real(E(k));f(k)=imag(E(k));V(k)=sqrt(e(k)^2+f(k)^2);sida(k)=atan(f(k)./e(k))*180./pi;end%=============== 计算各输出量=========================== disp('各节点的实际电压标幺值E为(节点号从小到大排列):'); disp(E); %显示各节点的实际电压标幺值E用复数表示disp('-----------------------------------------------------')disp('各节点的电压大小V为(节点号从小到大排列):');disp(V); %显示各节点的电压大小V的模值disp('-----------------------------------------------------');disp('各节点的电压相角sida为(节点号从小到大排列):');disp(sida); %显示各节点的电压相for p=1:nfor q=1:nC(p)=C(p)+conj(Y(p,q))*conj(E(q)); %计算各节点的注入电流的共轭值endS(p)=E(p)*C(p); %计算各节点的功率S = 电压X 注入电流的共轭值enddisp('各节点的功率S为(节点号从小到大排列):');disp(S); %显示各节点的注入功率Sline=zeros(n1,5);disp('-----------------------------------------------------');disp('各条支路的首端功率Si为(顺序同您输入B1时一致):');for i=1:n1p=B1(i,1);q=B1(i,2);Sline(i,1)=B1(i,1);Sline(i,2)=B1(i,2);if B1(i,6)==0Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));Siz(i)=Si(p,q);elseSi(p,q)=E(p)*(conj(E(p))*((1-B1(i,5))/B1(i,3))+(conj(E(p))-conj(E(q)))*(B1(i,5)/B1(i,3)));Siz(i)=Si(p,q);endSSi(p,q)=Si(p,q);Sline(i,3)=Siz(i);ZF=['S(',num2str(p),',',num2str(q),')=',num2str(SSi(p,q))];disp(ZF);enddisp('-----------------------------------------------------');disp('各条支路的末端功率Sj为(顺序同您输入B1时一致):');for i=1:n1p=B1(i,1);q=B1(i,2);if B1(i,6)==0Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));Sjy(i)=Sj(q,p);elseSj(q,p)=E(q)*(conj(E(q))*((B1(i,5)*(B1(i,5)-1))/B1(i,3))+(conj(E(q))-conj(E(p)))*(B1(i,5)/B1(i,3)));Sjy(i)=Sj(q,p);endSSj(q,p)=Sj(q,p);Sline(i,4)=Sjy(i);ZF=['S(',num2str(q),',',num2str(p),')=',num2str(SSj(q,p))];disp(ZF);enddisp('-----------------------------------------------------');disp('各条支路的功率损耗DS为(顺序同您输入B1时一致):');for i=1:n1p=B1(i,1);q=B1(i,2);DS(i)=Si(p,q)+Sj(q,p);DDS(i)=DS(i);Sline(i,5)=DS(i);ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DDS(i))];disp(ZF);enddisp('-----------------------------------------------------');disp('各支路首端编号末端编号首端功率末端功率线路损耗');disp(Sline);六、运行结果及其分析是否采用默认数据?(1-默认数据;2-手动输入)1导纳矩阵为:2.9056 -11.5015i 0.0000 + 5.3173i -1.6606 +3.1617i -1.2450 + 2.3710i0.0000 + 5.3173i 0.0000 - 4.6633i 0.0000 + 0.0000i 0.0000 + 0.0000i-1.6606 + 3.1617i 0.0000 + 0.0000i 2.4904 - 4.7039i -0.8298 + 1.5809i-1.2450 + 2.3710i 0.0000 + 0.0000i -0.8298 + 1.5809i 2.0749 - 3.9089i初始雅可比矩阵为:11.1267 2.7603 -5.3173 0 -3.1617 -1.6606-3.0509 11.8762 0 -5.3173 1.6606 -3.1617-5.3173 0 5.3173 0 0 00 -5.3173 0 4.0092 0 0-3.3198 -1.7436 0 0 4.8217 2.69800 0 0 0 0 2.1000迭代次数为:4收敛时电压修正量为::1.0e-05 *0.0349-0.2445-0.0101-0.5713-0.0931-0.0073各节点的实际电压标幺值E为(节点号从小到大排列):0.9673 - 0.0655i 1.0252 - 0.1666i 1.0495 - 0.0337i 1.0500 + 0.0000i -----------------------------------------------------各节点的电压大小V为(节点号从小到大排列):0.9695 1.0387 1.0500 1.0500-----------------------------------------------------各节点的电压相角sida为(节点号从小到大排列):-3.8734 -9.2315 -1.8419 0各节点的功率S为(节点号从小到大排列):-0.0000 + 0.0000i -0.5000 - 0.3000i 0.2000 + 0.1969i 0.3277 + 0.0443i -----------------------------------------------------各条支路的首端功率Si为(顺序同您输入B1时一致):S(1,2)=-0.5-0.30713iS(1,3)=-0.24266-0.197iS(1,4)=-0.25734-0.11013iS(3,4)=-0.055551+0.0017528i-----------------------------------------------------各条支路的末端功率Sj为(顺序同您输入B1时一致):S(2,1)=0.5+0.24606iS(3,1)=0.25555+0.1952iS(4,1)=0.2712+0.1014iS(4,3)=0.056496-0.057061i-----------------------------------------------------各条支路的功率损耗DS为(顺序同您输入B1时一致):DS(1,2)=0-0.06107iDS(1,3)=0.012892-0.0018014iDS(1,4)=0.013863-0.0087295iDS(3,4)=0.00094545-0.055308i-----------------------------------------------------各支路首端编号末端编号首端功率末端功率线路损耗1.0000 + 0.0000i2.0000 + 0.0000i -0.5000 - 0.3071i 0.5000 + 0.2461i 0.0000 - 0.0611i 1.0000 + 0.0000i3.0000 + 0.0000i -0.2427 - 0.1970i 0.2556 + 0.1952i 0.0129 - 0.0018i 1.0000 + 0.0000i4.0000 + 0.0000i -0.2573 - 0.1101i 0.2712 + 0.1014i 0.0139 - 0.0087i3.0000 + 0.0000i4.0000 + 0.0000i -0.0556 + 0.0018i 0.0565 - 0.0571i 0.0009 - 0.0553i七、实验体会及感悟通过这次实验,首先让我对matlab软件有了初步的了解,对它强大的矩阵运算能力有了更深的体会,同时掌握了设置断点和断点调试的一般方法,结合课本上的程序流程图和参考资料上的例子单步跟踪调试,再一次的熟悉了牛顿拉夫逊法潮流计算的一般方法和步骤,对计算机计算潮流计算有了更进一步的认识,在学习潮流计算时,虽然依次学习了节点导纳矩阵,功率方程、雅可比矩阵,但不能将它们联系起来,更不知道其中的原委,通过程序的编写,知道了其中的联系,也知道了每个方程、矩阵在计算中的作用。
简单电力系统分析潮流计算

简单电力系统分析潮流计算电力系统潮流计算是电力系统分析中的一项重要任务。
其目的是通过计算各个节点的电压、电流、有功功率、无功功率等参数,来确定系统中各个元件的运行状态和互相之间的相互影响。
本文将介绍电力系统潮流计算的基本原理、计算方法以及应用。
潮流计算的基本原理是基于电力系统的节点电压和支路功率之间的网络方程。
通过对节点电压进行迭代计算,直到满足所有支路功率平衡方程为止,得到系统的运行状态。
潮流计算的基本问题可以表示为以下方程组:P_i = V_i * (G_i * cos(θ_i - θ_j ) + B_i * sin(θ_i -θ_j )) - V_j * (G_i * cos(θ_i - θ_j ) - B_i * sin(θ_i -θ_j )) (1)Q_i = V_i * (G_i * sin(θ_i - θ_j ) - B_i * cos(θ_i -θ_j )) - V_j * (G_i * sin(θ_i - θ_j ) + B_i * cos(θ_i -θ_j )) (2)其中,P_i为节点i的有功功率注入;Q_i为节点i的无功功率注入;V_i和θ_i分别为节点i的电压幅值和相角;V_j和θ_j分别为节点j的电压幅值和相角;G_i和B_i分别为支路i的导纳的实部和虚部。
对于一个电力系统,如果知道了节点注入功率和线路的导纳,就可以通过潮流计算求解出各节点的电压和功率。
这是一种不断迭代的过程,直到系统达到平衡状态。
潮流计算的方法有多种,常见的有高斯-赛德尔迭代法、牛顿-拉夫逊迭代法等。
其中,高斯-赛德尔迭代法是最常用的一种方法。
高斯-赛德尔迭代法的思想是从已知节点开始,逐步更新其他节点的电压值,直到所有节点的电压值收敛为止。
具体步骤如下:1.初始化所有节点电压的初始值;2.根据已知节点的注入功率和节点电压,计算其他节点的电压值;3.判断节点电压是否收敛,如果收敛则结束计算,否则继续迭代;4.更新未收敛节点的电压值,返回步骤2高斯-赛德尔迭代法的优点是简单有效,但其收敛速度较慢。
电力系统分析第04章复杂电力系统潮流计算

电力系统分析第04章复杂电力系统潮流计算潮流计算是电力系统分析的一个重要工具,用于计算电力系统中各节点的电压幅值和相角,以及各支路的功率潮流分布情况。
复杂电力系统潮流计算主要包括节点潮流计算和线路潮流计算两部分。
节点潮流计算是指计算电力系统各节点的电压幅值和相角。
节点潮流计算的基本原理是根据节点复功率方程和节点电流平衡方程,建立节点潮流计算的数学模型。
该模型可以用于计算电力系统中各节点的电压幅值和相角,并找出潮流计算过程中出现的问题。
线路潮流计算是指计算电力系统中各支路的功率潮流分布情况。
线路潮流计算的基本原理是根据支路潮流方程,建立线路潮流计算的数学模型。
该模型可以用于计算电力系统中各支路的功率潮流,包括有功功率、无功功率和视在功率等。
在复杂电力系统潮流计算中,需要考虑以下几个方面。
首先,需要确定电力系统的潮流计算方法,常用的有直接法、迭代法和改进迭代法。
直接法适用于小型电力系统,计算速度较快,但对于大型电力系统不太适用。
迭代法采用不断迭代的方式计算潮流,适用于大型电力系统,计算精度较高。
改进迭代法是对迭代法的改进,可以提高计算速度和精度。
其次,需要确定电力系统的节点类型。
电力系统中的节点可以分为平衡节点、PQ节点、PV节点和参考节点。
平衡节点的有功功率和无功功率都为零,用于维持整个系统的功率平衡。
PQ节点的有功功率和无功功率是已知的,需要通过潮流计算来确定该节点的电压幅值和相角。
PV节点的有功功率是已知的,需要通过潮流计算来确定该节点的无功功率和电压幅值。
参考节点是一个已知电压值的节点,作为其他节点电压相角的参考点。
最后,需要考虑电力系统潮流计算的收敛性和稳定性。
收敛性是指潮流计算的结果是否能够收敛到一个稳定的值。
如果潮流计算不能收敛,则需要调整潮流计算的参数或算法,以提高收敛性。
稳定性是指潮流计算结果对电力系统的扰动是否具有稳定的响应。
如果潮流计算结果不稳定,则需要进一步分析系统的动态行为,以寻找稳定的解决方案。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
潮流计算的意义 (1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。 (2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。 (3)正常检修及下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。 (4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。 总结为在和规划方案的研究中,都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时,为了电力系统的运行状态,也需要进行大量而快速的潮流计算。因此,潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时,采用离线潮流计算;在的实时监控中,则采用在线潮流计算。 潮流计算的发展史 利用电子计算机进行潮流计算从20世纪50年代中期就已经开始。此后,潮流计算曾采用了各种不同的方法,这些方法的发展主要是围绕着对潮流计算的一些基本要求进行的。对潮流计算的要求可以归纳为下面几点: (1)算法的可靠性或收敛性 (2)计算速度和内存占用量 (3)计算的方便性和灵活性 电力系统潮流计算属于稳态分析范畴,不涉及系统元件的动态特性和过渡过程。因此其数学模型不包含微分方程,是一组高阶非线性方程。非线性代数方程组的解法离不开迭代,因此,潮流计算方法首先要求它是能可靠的收敛,并给出正确答案。随着电力系统规模的不断扩大,潮流问题的方程式阶数越来越高,目前已达到几千阶甚至上万阶,对这样规模的方程式并不是采用任何数学方法都能保证给出正确答案的。这种情况促使电力系统的研究人员不断寻求新的更可靠的计算方法。 在用数字计算机求解电力系统潮流问题的开始阶段,人们普遍采用以节点导纳为基础的高斯-赛德尔迭代法(一下简称导纳法)。这个方法的原理比较简单,要求的数字计算机的内存量也比较小,适应当时的制作水平和电力系统理论水平,于是电力系统计算人员转向以阻抗矩阵为主的逐次代入法(以下简称阻抗法)。 20世纪60年代初,数字计算机已经发展到第二代,计算机的内存和计算速度发生了很大的飞跃,从而为阻抗法的采用创造了条件。阻抗矩阵是满矩阵,阻抗法要求计算机储存表征系统接线和参数的阻抗矩阵。这就需要较大的内存量。而且阻抗法每迭代一次都要求顺次取阻抗矩阵中的每一个元素进行计算,因此,每次迭代的计算量很大。 阻抗法改善了电力系统潮流计算问题的收敛性,解决了导纳法无法解决的一些系统的潮流计算,在当时获得了广泛的应用,曾为我国电力系统设计、运行和研究作出了很大的贡献。但是,阻抗法的主要缺点就是占用计算机的内存很大,每次迭代的计算量很大。当系统不断扩大时,这些缺点就更加突出。为了克服阻抗法在内存和速度方面的缺点,后来发展了以阻抗矩阵为基础的分块阻抗法。这个方法把一个大系统分割为几个小的地区系统,在计算机内只需存储各个地区系统的阻抗矩阵及它们之间的联络线的阻抗,这样不仅大幅度的节省了内存容量,同时也提高了节省速度。 克服阻抗法缺点的另一途径是采用牛顿-拉夫逊法(以下简称)。牛顿法是数学中求解非线性方程式的典型方法,有较好的收敛性。解决电力系统潮流计算问题是以导纳矩阵为基础的,因此,只要在迭代过程中尽可能保持方程式的稀疏性,就可以大大提高牛顿潮流程序的计算效率。自从20世纪60年代中期采用了最佳顺序消去法以后,牛顿法在收敛性、内存要求、计算速度方面都超过了阻抗法,成为直到目前仍被广泛采用的方法。 在牛顿法的基础上,根据电力系统的特点,抓住主要矛盾,对纯数学的牛顿法进行了改造,得到了P-Q分解法。P-Q分解法在计算速度方面有显着的提高,迅速得到了推广。 牛顿法的特点是将非线性方程线性化。20世纪70年代后期,有人提出采用更精确的模型,即将泰勒级数的高阶项也包括进来,希望以此提高算法的性能,这便产生了保留非线性的潮流算法。另外,为了解决病态潮流计算,出现了将潮流计算表示为一个无约束非线性规划问题的模型,即非线性规划潮流算法。 近20多年来,潮流算法的研究仍然非常活跃,但是大多数研究都是围绕改进牛顿法和P-Q分解法进行的。此外,随着人工智能理论的发展,遗传算法、、模糊算法也逐渐被引入潮流计算。但是,到目前为止这些新的模型和算法还不能取代牛顿法和P-Q分解法的地位。由于电力系统规模的不断扩大,对计算速度的要求不断提高,计算机的并行计算技术也将在潮流计算中得到广泛的应用,成为重要的研究领域。 潮流计算的发展趋势 通过几十年的发展,潮流算法日趋成熟。近几年,对潮流算法的研究仍然是如何改善传统的潮流算法,即高斯-塞德尔法、牛顿法和快速解耦法。牛顿法,由于其在求解非线性潮流方程时采用的是逐次线性化的方法,为了进一步提高算法的收敛性和计算速度,人们考虑采用将泰勒级数的高阶项或非线性项也考虑进来,于是产生了二阶潮流算法。后来又提出了根据直角坐标形式的潮流方程是一个二次代数方程的特点,提出了采用直角坐标的保留非线性快速潮流算法。 对于保留非线性算法典型论文有: 1.文献[保留非线性的电力系统概率潮流计算]提出了它在电力系统概率潮流计算中的应用。该文献提出了一种新的概率潮流计算方法,它保留了潮流方程的非线性,又利用了P-Q解耦方法,因而数学模型精度较高,且保留了P-Q解耦的优点,有利于大电网的随机潮流计算,用提出的方法对一个典型的系统进行了计算,其数值用MonteCarlo随机模拟作了验证,得到了满意的结果。 2.文献[基于系统分割的保留非线性的快速P-Q解耦潮流计算法]分析研究了保留非线性的P-Q解耦快速潮流计算法。该文献提出了一种新的状态估计算法,既保留了量测方程非线性又利用了快速P-Q分解方法,因此数学模型精度高且保留了快速P-Q分解的优点,提高了状态估计的计算精度和速度.采用系统分割方法将大系统分割为多个小系统,分别对每个小系统进行状态估计,然后对各小系统的状态估计结果进行协调,得到整个系统具有同一参考节点的状态估计结果,这样可大大提高状态估计的计算速度,有利于进行大电网的状态估计.在18节点系统上进行的数字仿真实验验证了该方法的有效性。岩本伸一等提出了一种保留非线性的快速潮流计算法,但用的是,因而没法利用P-Q解耦。为了更有利于大电网的潮流计算,将此原理推广用于P-Q解耦。这样,既利用了保留非线性的快速算法,在迭代中使用常数,又保留了P-Q解耦的优点。 对于一些病态系统,应用非线性潮流计算方法往往会造成计算过程的振荡或者不收敛,从数学上讲,非线性的潮流计算方程组本来就是无解的。这样,人们提出来了将潮流方程构造成一个函数,求此函数的最小值问题,称之为非线性规划潮流的计算方法。优点是原理上保证了计算过程永远不会发散。如果将数学规划原理和牛顿潮流算法有机结合一起就是最优乘子法。另外,为了优化系统的运行,从所有以上的可行潮流解中挑选出满足一定指标要求的一个最佳方案就是最优潮流问题。最优潮流是一种同时考虑经济性和安全性的分析优化问题。OPF 在电力系统的安全运行、经济调度、可靠性分析、能量管理以及电力定价等方面得到了广泛的应用。 最优潮流方面的典型论文有: 1.文献[电力系统最优潮流新算法的研究]以NCP 方法为基础,提出了一种新的求解最优潮流算法——投影渐近半光滑牛顿型算法。该文献以NCP方法为基础,提出了一种新的求解OPF算法——投影渐近半光滑牛顿型算法。针对电力系统的特点,本文的研究工作如下: 1.建立了与OPF问题的KKT系统等价的带界约束的半光滑方程系统。与已有的NCP方法相比,新的模型由于无需考虑界约束对应的对偶变量(乘子变量),降低了问题的维数,从而适用于解大规模的电力系统问题。 2.基于建立的新模型,本文提出了一类新的Newton型算法,该算法一方面保持界约束的相容性,另一方面有较好的全局与局部超线性收敛性,同时,算法结构简单,易于实现。 3.考虑到电力系统固有的弱耦合特性,受传统解耦最优潮流方法的启示,在所提出的新Newton型方法的基础上,本文又设计了一类分解方法。新方法基于解耦——校正的策略实现算法,不仅充分利用了系统的弱耦合特性,同时保证分解算法在理论上的收敛性。 4.根据所提出的两种算法,用标准的IEEE电力测试系统进行数值实验,并与已有的其他方法进行比较。结果显示新算法具有良好的收敛性和计算效果,在电力系统的规划与运行方面将有广阔的应用前景。 2.文献[基于可信域内点法的最优潮流问题研究]介绍了OPF内点法具有收敛性强、多项式时间复杂性等优点,是极具潜力的优秀算法之一。 电力系统不断发展,使得OPF算法跻身于极其困难、非凸的大规模非线性规划行列。可信域和线性搜索方法是保证最优化算法全局收敛性能的两类技术,将内点法和可信域、线性搜索方法有机结合,构造新的优化算法,是数学规划领域的研究热点。 此方面的典型文献有: 1.文献[电力市场环境下基于最优潮流的输电容量充裕度研究]首先以最优潮流为工具,选取系统中的关键线路作为系统输电容量充裕度的研究对象,从电网运行的安全性、可靠性的角度系统地研究了稳定限额对输电容量充裕度的影响,指出稳定限额因子与影子