潮流计算常用线路参数
电力系统潮流计算(九节点)
辽宁工程技术大学电力系统分析课程设计设计题目9节点电力网络潮流计算指导教师院（系、部）专业班级学号姓名日期电气工程系课程设计标准评分模板电力系统分析课程设计任务书9节点系统单线图如下：基本数据如下：母线名基准电压区域号电压上限电压下限发电 1 16.5000 2 18.1500 14.8500 发电 2 18.000 1 19.800 16.2000 发电 3 13.8000 1 15.1800 12.4200 GEN1-230 230.000 2 0.0000 0.0000 GEN2-230 230.000 1 0.0000 0.0000 GEN3-230 230.000 1 0.0000 0.0000 STNA-230 230.000 2 0.0000 0.0000 STNB-230 230.000 2 0.0000 0.0000 STNC-230 230.000 1 0.0000 0.0000数据组I 侧母线J 侧母线编号所属区域单位正序电阻正序电抗正序充电电纳的1/2常规GEN1-230 STNA-230 1 I侧标么0.010000 0.085000 0.088000 常规STNA-230 GEN2-230 2 I侧标么0.032000 0.161000 0.153000 常规GEN2-230 STNC-230 3 I侧标么0.008500 0.072000 0.074500 常规STNC-230 GEN3-230 4 I侧标么0.011900 0.100800 0.104500 常规GEN3-230 STNB-230 5 I侧标么0.039000 0.170000 0.179000 常规STNB-230 GEN1-230 6 I侧标么0.017000 0.092000 0.079000表3 两绕组变压器数据负荷数据电网12-1班数据目录1 PSASP软件简介 (1)1.1 PSASP平台的主要功能和特点 (7)1.2 PSASP的平台组成 (8)2 牛顿拉夫逊潮流计算简介 (9)2.1 牛顿—拉夫逊法概要 (9)2.2 直角坐标下的牛顿—拉夫逊潮流计算 (11)2.3 牛顿—拉夫逊潮流计算的方法 (6)3 九节点系统单线图及元件数据 (8)3.1 九节点系统单线图 (8)3.2 系统各项元件的数据 (9)4 潮流计算的结果 (11)4.1 潮流计算后的单线图 (17)4.2 潮流计算结果输出表格 (18)5 结论 (22)6 参考文献 (17)1 PSASP软件简介“电力系统分析综合程序”（Power System Analysis Software Package,PSASP）是一套历史悠久、功能强大、使用方便的电力系统分析程序，是高度集成和开发具有我国自主知识产权的大型软件包。
含分布式电源的配电网潮流计算
含分布式电源的配电网潮流计算一、本文概述随着可再生能源的快速发展和广泛应用，分布式电源（Distributed Generation，DG）在配电网中的渗透率逐年提高。
分布式电源包括风力发电、光伏发电、微型燃气轮机等，它们具有位置灵活、规模适中、与环境兼容性强等特点，是智能电网的重要组成部分。
然而，分布式电源的接入对配电网的潮流分布、电压质量、系统稳定性等方面都产生了显著影响。
因此，准确进行含分布式电源的配电网潮流计算，对于保障配电网安全、经济运行具有重要意义。
本文旨在探讨含分布式电源的配电网潮流计算方法。
本文将对分布式电源的类型、特性及其在配电网中的应用进行简要介绍。
将重点分析分布式电源接入对配电网潮流计算的影响，包括电源位置、容量、出力特性等因素。
在此基础上，本文将提出一种适用于含分布式电源的配电网潮流计算模型和方法，并对其准确性、有效性进行验证。
本文还将对含分布式电源的配电网潮流计算在实际工程中的应用前景进行讨论。
通过本文的研究，旨在为配电网规划、运行和管理人员提供一套有效的潮流计算工具和方法，以应对分布式电源大量接入带来的挑战。
本文的研究成果也有助于推动智能电网、可再生能源等领域的技术进步和应用发展。
二、分布式电源建模在配电网潮流计算中，分布式电源（Distributed Generation，DG）的建模是至关重要的一步。
分布式电源通常包括风能、太阳能、小水电、生物质能等多种类型，它们的接入位置和容量对配电网的潮流分布、电压质量、系统稳定性等方面都有显著影响。
建模过程中，首先需要明确分布式电源的类型和特性。
例如，对于光伏电源，其输出功率受到光照强度、温度等自然条件的影响，具有随机性和波动性；而对于风力发电，其输出功率则受到风速、风向、湍流强度等因素的影响，同样具有不确定性。
因此，在建模时需要考虑这些不确定性因素，以更准确地描述分布式电源的实际运行状况。
需要根据分布式电源的具体接入方式和位置，建立相应的数学模型。
潮流计算
U U % 1 N 100 U 1.N
UN
U U % 2 N 100 U 2.N
UN
第三章 输电线路运行特性及简单电力系统潮流估算
2. 电力线路的功率损耗
下图为电力线路的П型等值电路，其中Z=R+jX，Y=G+jB
为电力线路每相阻抗和导纳。一般情况下，线路一端的功率和 电压是已知的，要求计算另一端的功率和电压。
*
dU
U1
U 2
S2
*
(R
jX )
U2
(3 2)
第三章 输电线路运行特性及简单电力系统潮流估算
*
从而得：U1
U 2
S2
*
(R
jX )
பைடு நூலகம்U2
取U 2为参考向量，即令U 2 U200 ,则式（3 - 2）变成
*
dU 2
S2
*
U2
(R
jX )
P2
jQ2 U2
(R
jX )
P2R Q2 X j P2 X Q2R
3)计
算变压器始
端电压U1
(以U
为参考向量
2
）
U1
U 2
P2 RT Q2 X T U2
j
P2 X T Q2RT U2
4)计算变压器并联支路的功率损耗
S~YT PYT jQYT YTU1 *U1 (GT jBT )*U12 (GT jBT )U12
5)计算变压器的输入功率
~ S1
~
S1
S 1
~
.
y1
U1
~
~
S S ' '
1
Z
YZ
2
电网的潮流计算
（2）容纳
电力网线路的导线间，在电压作用下形成 电场，导线上电荷与电压的比值称为线路 里的容纳b＝WC。 若线路长度为L（KM），则每相线路容纳 B＝bL，每相线路容性电流Ic按下式计算： Ic＝UxB（KA） Ux――线路相电压
单位:西门子(S)
LGJ－ LGJ－ LGJ－ LGJ－ 120 150 185 240
0.27 0.21 0.17 0.13
LGJ－ LGJ－ 300 400
0.107 0.08
感 0.446 0.435 0.429 0.422 0.415 0.407 0.399 0.391 抗x
r、x
x r
S
LGJ导线r、x与S的关系曲线
2、一端电源供电网的潮流计算。
（4）电导
线路的电导是指对应于相间绝缘介质损耗 功率的参数，对于架空线路，相间绝缘介 质损耗主要为电晕损耗与绝缘子泄漏损耗； 对于电缆线路，主要为绝缘介质的极化与 泄漏损耗。
架空线路由绝缘子漏电所产生的有功损耗 很小，可以忽略不计，因此，架空线路相 间绝缘介质损耗主要有电晕所决定，电晕 所损耗的有功功率一般用参数电导表征。
计算出电力系统在正常及各种可能的故 障运行方式下的潮流分布，各节点电压及 元件中的功率损耗，对于电力系统设计、 运行都是十分必要的。
目前，计算机已广泛应运于电力系统的运 行、设计和科学研究各个方面。自1956年 成功地运用它计算潮流分布以来，几乎所 有主要的电力系统计算都已使用计算机。
1、电力网的功率损耗与电压计算。
（1）施工时不要磨损导线，要保持导线及金属元 件表面光滑，以防电场不均匀。
（2）增大导线半径，减小导线表面附近的电场强 度，可采用分裂导线、扩径导线、空心导线等。
电力系统潮流计算
3.2.1 节点电压方程与节点导纳矩阵和阻抗矩阵
将节点电压法应用于电力系统潮流计算，变量为节点电压与节
点注入电流。通常以大地作为电压幅值的参考（U0 = 0），以
系统中某一指定母线的电压角度作为电压相角的参考，以支路
导纳作为电力网的参数进行计算。节点注入电流规定为流向网
络为正，流出为负。
Pmax P
表征年有功负荷曲线特点的两个指标
0
年最大负荷利用小时数 Tmax
t Tmax 8760
根据年负荷曲线，可求得全年所需电能：
8760
A 0
Pdt MWh
定义年最大负荷（最大值 Pmax）利用小时： Tmax
A Pmax
h
Tmax 越大，负荷曲线越平坦
负荷曲线为一水平线时， Tmax 达到最大值8760 (h)
2
1 ZT1
2
Zl
T2
34
3
ZT2 4
YT3
Yl /2
YT2
已知末端功率和电压， 计算网上潮流分布。
1 ZT1 2 Zl
3 ZT2 4
已知始端功率和电压， 计算网上潮流分布。
Y20
Y30
已知末端功率和始端电 压，计算网上的潮流。
不管哪种情况，先作等值电路
3.1.3 辐射形网络的分析计算
1）已知末端功率、电压 利用前面的方法，从末端逐级 往上推算，直至求得各要求的量。
Pm(t)
损耗称年电能损耗，是电网运行经
济性的指标。
Pmi
1）年电能损耗的准确计算方法
已知各负荷的年有功和无功负荷曲线 时，理论上可准确计算年电能损耗。
8760小时分为 n 段，第 i 时段时间为 Dti (h)，全网功率损耗为DPi (MW)，则 全网年电能损耗为
潮流计算 例题
例题：如图1所示的简单电力网中，已知变压器的参数为S N =31.5MV A ，0S S 031kW,190kW,%=10.5,%=0.7P P U I ∆=∆=；线路单位长度的参数为61110.21/km,0.416/km, 2.7410S/km r x b -=Ω=Ω=⨯。
如图所示的简单电力网中，当线路首端电压U A =120kV 时，试求：（1）线路和变压器的电压损耗；（2）变压器运行在额定变比时的低压侧电压及电压偏移。
说明：以上计算忽略电压降落的横分量。
图1解：如题画等值电路图如下：线路参数为：0.21408.40.4164016.64l l l l R rl X x l ==⨯=Ω==⨯=Ω变压器参数为Ω=⨯⨯⨯=⨯∆=317.210)105.31(110190103232322N N S T S U P R Ω=⨯⨯⨯⨯⨯=33.4010105.311001105.1010100%33232N N S T S U U X （1） 变压器的功率损耗和励磁功率为222T S 02N 2515()19031193.760.1937631.5S P P P kW MW S +∆=∆+∆=⨯+==222S 0N T N %%10.5(2515)0.731.5 3.0538var 10010010031.5100U S I S Q M S ⨯+⨯∆=+=+=⨯1点处线路的充电功率var 66308.01104074.22121222M lU b Q N l B =⨯⨯⨯==计算L S 2 为：MVAj j Q Q Q j P P S B T LD T LD L 39.1719.25)66308.00538.315(19376.025)(22+=-+++=-∆++∆+=线路阻抗中的功率损耗为：MW R U Q P P l L L l 65044.0104021.011017390251901032223222222=⨯⨯⨯+=⨯+=∆-- v a r 2885.11040416.011017390251901032223222222M X U Q P Q l L L l =⨯⨯⨯+=⨯+=∆-- 计算功率1S 为 M V Aj j Q Q j P P S l L l I 68.1884.25)2885.139.17(65044.019.25)(221+=+++=∆++∆+=线路电压损耗（忽略电压降落的横分量） 1125.848.418.6816.64 4.40120=l l l A PR Q X U kV U +⨯+⨯∆== 1点电压为：1120 4.40115.60-A l U U U kV =∆=-=计算功率TS 2 为 M V Aj j Q Q j P P S T LD T LD T 833.17163.25)833.215(163.025)(2+=+++='∆++'∆+= 变压器电压损耗 22125.163 2.3217.83340.33 6.73115.60=T T T T T P R Q X U kV U +⨯+⨯∆== （2） 变压器低压侧折算到高压侧的电压为21115.60 6.73108.87=-T U U U kV '∆=-= 变压器低压侧的实际电压 22108.8710.8910=U U kV k '== 电压偏移为2210.8910%100%8.9%10N N U U m U --=⨯==。
BPA潮流计算解析
节点数据 (B卡 )：发电机节点、负荷节点、中枢母线 节点
支路数据 ：线路参数(L卡)、变压器支路参数(T卡)
§2.1 潮流文件编制
(1) 一级控制语句 一级控制语句用“(……)”表示，表示潮流作业的开始和结束。 (2) 二级控制语句 二级控制语句用“/……\”表示，.1 潮流文件编制 §2.2 地理接线图绘制 §2.3 潮流计算深入分析
§2 BPA潮流计算
BPA程序简介 程序规模
– 1) 10000个节点; – 2) 15200条支路，其中包括300台带负荷调
压（LTC）变压器，70台移相器; – 3) 40条两端直流线路和25条多端直流线路; – 4) 120个功率交换区域，所有的区域共可包
根据国调中心调研情况，全国30个网、省（市）调 中 ， 日 常 采 用 BPA 进 行 潮 流 计 算 的 有 12 家 、 日 常 采 用 PASAP进行潮流计算的有19家,日常采用PSS/E进行潮流 计算的有5家（其中有部分调度机构根据需要采用了两 种以上工具进行潮流计算） 。华北、华东电网内采用 BPA较多，其它区域采用PSASP较多。
BPA潮流计算
2013年3月
§1 概述
安全稳定计算分析（《国家电网安全稳定 计算规范》） ➢潮流分析
➢无功电压分析 ➢静态安全分析 ➢静态稳定计算分析
➢ 暂态稳定计算分析 ➢ 小扰动动态稳定分析
➢短路电流安全校核 ➢电压稳定计算分析 ➢频率稳定计算分析 ➢再同步稳定计算分析
§1 概述
电力系统静态安全分析
§1 概述
BPA软件的应用介绍
➢ BPA软件的应用 潮流计算、暂态稳定计算、最优潮流、短路电
4、电力系统潮流计算
由前刚求得的S1和已 知的首端电压U1，可 求得阻抗上压降：
又
第四章
开式电力网的潮流计算 可求得该段 阻抗环节末 端电压为： 推广：
任何复杂电网的等值电路，均可分解成多个阻 抗环节； 求解复杂电网的潮流，可按上述一段阻抗环节 的功率、电压算法，逐个环节重复、递推计算。
电 力 系 潮 流 计 算
练习--课本例题4-1及补充习题-1
电 力 系 潮 流 计 算
潮流计算——电力系统分析（稳态分析）中一 种最基本的计算。 稳态分析计算—— 不考虑发电机的参数，将机端母线作为系 统的边界点。
潮流的基本概念
在发电机母线上功率被注入网络； 而在变（配）电站的母线上接入负荷； 两者之间，功率在网络中流动。 对于这种流动的功率，电力生产部门称之为 潮
电 力 系 潮 流 计 算
将电压和功率由末端向始端交替推进 ；
对于110KV及以下网络，可略去电压降落 的横分量，从而使计算简化； 计算中须注意到变压器参数及电压的归算。
第四章
开式电力网的潮流计算
练一练：[补充-2］
电 力 系 潮 流 计 算
一电力线路长100km，末端接有一台容量为 20MVA、变比为110/38.5KV的降压变压器； 变压器低压侧负荷为15+j11.25MVA。 正常运行时负荷要求电压为36KV，试求线路 始端母线应具备多大的功率和电压才满足要 求。 2 T 3 1
（4-25）
当s=2，即已知环节末端量来求阻抗损耗的计 算式，即为：
（4-26）
（4-27）
第四章
（一）电力网的功率损耗
即电容功率 )
电 力 系 潮 流 线路首端的无功损耗： 计 算
电力系统分析潮流计算最终完整版
电力系统分析潮流计算实验报告姓名：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软件有了初步的了解，对它强大的矩阵运算能力有了更深的体会，同时掌握了设置断点和断点调试的一般方法，结合课本上的程序流程图和参考资料上的例子单步跟踪调试，再一次的熟悉了牛顿拉夫逊法潮流计算的一般方法和步骤，对计算机计算潮流计算有了更进一步的认识，在学习潮流计算时，虽然依次学习了节点导纳矩阵，功率方程、雅可比矩阵，但不能将它们联系起来，更不知道其中的原委，通过程序的编写，知道了其中的联系，也知道了每个方程、矩阵在计算中的作用。
3.2电力网络潮流计算的手算解法要点
3.2电⼒⽹络潮流计算的⼿算解法要点3.2 电⼒⽹络潮流计算的⼿算解法3.2.1 电压降落及功率损耗计算1.电⼒线路上功率损耗与电压降落的计算电压是电能质量的指标之⼀，电⼒⽹络在运⾏过程中必须把某些母线上的电压保持在⼀定范围内，以满⾜⽤户电⽓设备的电压处于额定电压附近的允许范围内。
电⼒系统计算中常⽤功率⽽不⽤电流，这是因为实际系统中的电源、负荷常以功率形式给出，⽽电流是未知的。
当电流（功率）在电⼒⽹络中的各个元件上流过时，将产⽣电压降落，直接影响⽤户端的电压质量。
因此，电压降落的计算为分析电⼒⽹运⾏状态所必需。
电压降落即为该⽀路⾸末两端电压的相量差。
对如图3.3所⽰系统，已知末端相电压及功率求线路功率损耗及电压降落，设末端电压为，末端功率为，则线路末端导纳⽀路的功率损耗为（3-8）则阻抗末端的功率为阻抗⽀路中损耗的功率为，（3-9）阻抗⽀路始端的功率，线路始端导纳⽀路的功率损耗，（3-10）线路⾸端功率，从式（3-8）-（3-10）可知，线路阻抗⽀路有功功率和⽆功功率损耗均为正值，⽽导纳⽀路的⽆功功率损耗为负值，表⽰线路阻抗既损耗有功功率⼜损耗⽆功功率，导纳⽀路实际上是发出⽆功功率的（⼜称充电功率），充当⽆功功率源的作⽤，也就是说，当线路轻载运⾏时，线路只消耗很少的⽆功功率，甚⾄会发出⽆功功率。
⾼压线路在轻载运⾏时发出的⽆功功率，对⽆功缺乏的系统可能是有益的，但对于超⾼压输电线路是不利的，当线路输送的⽆功功率⼩于线路的充电功率时，线路始端电压可能会低于末端电压，或者说末端电压⾼于始端电压，若末端电压升⾼可能会导致绝缘的损坏，是应加以避免的，⼀般为了防⽌末端电压的升⾼，线路末端常连接有并联电抗器在轻载或空载时抵消充电功率，避免出现线路电压过⾼。
从以上推导不难看出，要想求出始端导纳⽀路的功率损耗及，必须先求出始端电压。
设与实轴重合，即，如图3-4所⽰。
图3-3 电⼒线路的电压和功率图3-4 利⽤末端电压计算始端电压则由(3-11)令则有（3-12）从⽽得出功率⾓在⼀般电⼒系统中，远远⼤于δU，也即电压降落的横分量的值δU对电压U1的⼤⼩影响很⼩，可以忽略不计，所以同理，也可以从始端电压、始端功率求取电压降落及末端电压和末端功率的计算公式。
电力系统分析潮流计算例题
电力系统的潮流计算西安交通大学自动化学院 2012。
103．1 电网结构如图3-11所示，其额定电压为10KV 。
已知各节点的负荷功率及参数： MVAj S )2.03.0(2+=，MVAj S )3.05.0(3+=，MVA j S )15.02.0(4+=Ω+=)4.22.1(12j Z ,Ω+=)0.20.1(23j Z ,Ω+=)0.35.1(24j Z试求电压和功率分布.解：（1）先假设各节点电压均为额定电压，求线路始端功率。
0068.00034.0)21(103.05.0)(22223232232323j j jX R V Q P S N +=++=++=∆0019.00009.0)35.1(1015.02.0)(22224242242424j j jX R V Q P S N +=++=++=∆则: 3068.05034.023323j S S S +=∆+=1519.02009.024424j S S S +=∆+=6587.00043.122423'12j S S S S +=++=又0346.00173.0)4.22.1(106587.00043.1)(22212122'12'1212j j jX R V Q P S N +=++=++=∆故： 6933.00216.112'1212j S S S +=∆+=（2） 再用已知的线路始端电压kVV 5.101=及上述求得的线路始端功率12S ，求出线路各点电压.kVV X Q R P V 2752.05.104.26933.02.10216.1)(11212121212=⨯+⨯=+=∆ kV V V V 2248.101212=∆-≈kVV V V kV V X Q R P V 1508.100740.0)(242422424242424=∆-≈⇒=+=∆kVV V V kV V X Q R P V 1156.101092.0)(232322323232323=∆-≈⇒=+=∆（3）根据上述求得的线路各点电压，重新计算各线路的功率损耗和线路始端功率.0066.00033.0)21(12.103.05.022223j j S +=++=∆ 0018.00009.0)35.1(15.1015.02.022224j j S +=++=∆ 故 3066.05033.023323j S S S +=∆+=1518.02009.024424j S S S +=∆+=则 6584.00042.122423'12j S S S S +=++=又0331.00166.0)4.22.1(22.106584.00042.122212j j S +=++=∆ 从而可得线路始端功率 6915.00208.112j S +=这个结果与第（1）步所得计算结果之差小于0。
潮流计算
Sb SG STc S0c jQB 2 jQB3
1 b Tb 2 c Tc 3
A
d Td
SLDb
G
SG
SL D d
14
二、两级电压的开式电力网计算 计算方法一：包含理想变压器，计算时，经过理 想变压器功率保持不变，两侧电压之比等于实际 变比k。 T b d c L-1 L-2 SLD A
V1 arctg V1 V1
4
网络元件的功率损耗
功率损耗包括：电阻和等值电抗上的损耗 对地等值导纳上产生的损耗
V1S1 , I1 S ' I
jQB1
B j 2
R jX
S '', I S 2 , I 2 V2
jQB 2
B j 2
线路
VS1 , I1
线路
S0 (GT jBT )V 2
I0% S0 P0 jQ0 P0 j SN 100
开式网络的电压和功率分布计算
一、已知供电点电压和负荷点功率时的计算方法 已知末端的功率和电压：从末端开始依次计算出 电压降落和功率损耗。 已知电源点的电压和负荷的功率：采取近似的方 法通过叠代计算求得满足一定精度的结果
X2 k2 X2
T
A
A
B2 B2 / k 2 d c L-2 SLD
R'2+ j X'2 j B'2/2
16
R1+ jX1
j B1/2 j B1/2
b ΔS0
Z'T
c' j B'2/2
d'
SLD
二、两级电压的开式电力网计算 计算方法三：用π型等值电路代表变压器
电力线路的参数与等值电路以及潮流计算的简单介绍
U 2 Z
P Q j 2
2
U* 2
两边同乘 e3 j30 U U
U 1 U 2 U U 2 Z
e P Q 3 ( j30 2 j
)
2
U* 2
U 2 Z
3(P2 jQ2)
3e
j30
U* 2
U 2
Z
P2 jQ2 U* 2
**
U
U 1 U 2 Z
P2 jQ2 电压降落 U2
基本概念
二、电压降落、电压损耗、电压偏移
目的：对于一条线路（变压器）有负荷流过时，首末端电压不等，造
成电压 损耗，可以推导已知端的S和U时求另一端的S和U
u 1
I
u 2
R jX
S 2 P2 jQ2
1、已知U2及S2求U1
I
S 2 U 2
*
P Q j 2
2
U* 2
U 1 U 2 U U 2 Z I
电力线路的参数与等值电路
一.单位长度电力线路的参数
1、电阻 r1=ρ/ s
ρ电阻率
单位：Ω•mm2/km 铜：18.8 铝: 31.3
与温度有关
S 截面积 mm2
一般是查表 rt=r20(1+α(t-20))
钢线电阻：导磁集肤、磁滞效应交流电阻> 直流电阻，和电流有关查手册
电力线路的参数与等值电路
以U2为参考电压
U
(R
jX ) P2 jQ2 U2
I2
U 1 U U' U 2 U'
P2 R Q2 X j P2 X Q2 R U' j U'
U2
U2
纵分量 横分量
2、已知U1及S1求U2
电气工程基础第五章
频率的变化对电力系统的正常运行也是十分有害的。频率下降会使发电厂的许多重要设备如给水泵、循环水泵、风机等的出力下降，造成水压、风力不足，使整个发电厂的有功出力减少，导致频率进一步下降，如不采取必要措施，就会产生所谓“频率崩溃”的恶性循环；频率的变化可能会使汽轮机的叶片产生共振，降低叶片寿命，严重时会产生裂纹甚至断片，造成重大事故。另外，频率的下降，会使异步电动机和变压器的励磁电流增大，无功损耗增加，给电力系统的无功平衡和电压调整增加困难。
一次调频为有差调频，频率不会恢复到初始值。一次调频时，系统中所有发电机组均参与。
二次调频的调节范围比一次调频大，可将频率恢复到偏移的允许范围或初始值。二次调频时，仅系统中特定的调频电厂（主调频厂与辅助调频厂）参与。
5-12如何选择主调频厂？
解：主调频厂一般应按下列条件选择：
①具有足够的调节容量和范围；
解：（1）两台变压器变比均为110/11时，各变压器通过的负荷功率；
（2）要使变压器T2满载运行，应如何调整变压器的变比？
(1)变比相同，可等效为电压相等的两端供电网，变压器的负载功率为：
（2）T2满载运行时，通过的负荷功率为
循环功率为
由
得
由
得到
或
5-9我国规定频率的额定值是多少？允许偏移值是多少？系统低频运行对用户和系统有什么危害？
在电力系统的设计和运行中都要用到潮流计算的结果，例如电力网规划设计时，要根据潮流计算的结果选择导线截面和电气设备，确定电力网主结线方案，计算网络的电能损耗和运行费用，进行方案的经济比较；电力系统运行时，要根据潮流计算的结果制定检修计划，校验电能质量，采取调频和调压措施，确定最佳运行方式，整定继电保护和自动装置。
电力系统潮流计算
电力系统潮流计算The final revision was on November 23, 2020电力系统课程设计题目: 电力系统潮流计算院系名称：电气工程学院专业班级：电气F1206班学生姓名：学号：指导教师：张孝远12节点的分类 (5)3 计算方法简介 (6)牛顿—拉夫逊法原理 (6)牛顿—拉夫逊法概要 (6)牛顿法的框图及求解过程 (8)MATLAB简介 (9)4 潮流分布计算 (10)系统的一次接线图 (10)参数计算 (10)丰大及枯大下地潮流分布情况 (14)该地区变压器的有功潮流分布数据 (15)重、过载负荷元件统计表 (17)5 设计心得 (17)参考文献 (18)附录:程序 (19)原始资料一、系统接线图见附件1。
二、系统中包含发电厂、变电站、及其间的联络线路。
500kV变电站以外的系统以一个等值发电机代替。
各元件的参数见附件2。
设计任务1、手动画出该系统的电气一次接线图，建立实际网络和模拟网络之间的联系。
2、根据已有资料，先手算出各元件的参数，后再用Matlab表格核算出各元件的参数。
3、潮流计算1）对两种不同运行方式进行潮流计算，注意110kV电网开环运行。
2）注意将电压调整到合理的范围110kV母线电压控制在106kV～117kV之间;220kV母线电压控制在220 kV～242kV之间。
附件一：72水电站2水电站1303x40C20+8B 2x8A2x31.5D4x7.5水电站5E2x1090+120H12.5+31.5FG1x31.5水电站324L2x150火电厂1x50M110kV线路220kV线路课程设计地理接线示意图110kV变电站220kV变电站牵引站火电厂水电站500kV变电站附件二：1、变压器：两个220kV变电站均采用参数一致的三绕组变压器，具体参数如下。
220kV变电站参数表110kV及以下的变电站的变压器省略，即可将负荷直接挂在110kV母线上。
潮流计算
武汉理工大学《电力系统分析》课程设计说明书节点数：4 支路数：4 计算精度：0.00010支路1：0.0200+j0.08001┠—————□—————┨3支路2：0.0400+j0.12001┠—————□—————┨4支路3：0.0500+j0.14002┠—————□—————┨4支路4：0.0400+j0.12003┠—————□—————┨4节点1：PQ节点，S(1)=-0.6000-j0.2500节点2：PQ节点，S(2)=-0.8000-j0.3500节点3：PV节点，P(3)=0.4000 V(3)=0.9500节点4：平衡节点，U(4)=1.0000∠0.0000运用matlab软件对选定课设题目进行潮流计算。
潮流计算是电力系统课程中必须掌握也是非常重要的计算。
潮流计算是指对电力系统正常运行状况的分析和计算。
在已知系统条件情况下，给定一些初始条件，进而计算出系统运行的电压和功率等；潮流计算方法很多：高斯-塞德尔法、牛顿-拉夫逊法、PQ分解法、直流潮流法等。
通过潮流计算，可以确定各母线的电压幅值和相角,各元件流过的功率和整个系统的功率损耗。
潮流计算是实现安全经济发供电的必要手段和重要工作环节。
因此潮流计算在电力系统的规划计算，生产运行，调度管理及科学计算中都有广泛的运用。
本课程设计采用PQ分解法进行电力系统分析的潮流计算程序的编制与调试，获得电力系统中各节点电压，为进一步进行电力系统分析作准备。
关键词：matlab 潮流计算PQ分解法1.题目原始数据及其化简 (1)2.PQ分解法 (2)2.1PQ分解法基本思想 (2)2.2 PQ分解法潮流计算基本步骤 (5)3编程及运行 (6)3.1 PQ分解法潮流计算程序框图 (6)3.2源程序代码 (7)3.3运行程序及结果分析： (16)4.小结 (18)5.参考文献 (19)1.题目原始数据及其化简原始数据：节点数：4 支路数：4 计算精度：0.00010 支路1：0.0200+j0.08001┠—————□—————┨3支路2：0.0400+j0.12001┠—————□—————┨4支路3：0.0500+j0.14002┠—————□—————┨4支路4：0.0400+j0.12003┠—————□—————┨4节点1：PQ节点，S(1)=-0.6000-j0.2500节点2：PQ节点，S(2)=-0.8000-j0.3500节点3：PV节点，P(3)=0.4000 V(3)=0.9500节点4：平衡节点，U(4)=1.0000∠0.0000根据原始数据所画电路简化图如图1：1 34 2图1电路简化图2.PQ 分解法2.1PQ 分解法基本思想PQ 分解法是从改进和简化牛顿法潮流程序的基础上提出来的，它的基本思想是：把节点功率表示为电压向量的极坐标形式，以有功功率误差作为修正电压向量角度的依据，以无功功率误差作为修正电压幅值的依据，这样，n-1+m 阶的方程式便分解为一个n-1阶和一个m 阶的方程，这两组方程分别进行轮流迭代，这就是所谓的有功-无功功率分解法。
第3章作业答案电力系统潮流计算(已修订)
第三章 电力系统的潮流计算3-1 电力系统潮流计算就是对给定的系统运行条件确定系统的运行状态。
系统运行条件是指发电机组发出的有功功率和无功功率（或极端电压），负荷的有 功功率和无功功率等。
运行状态是指系统中所有母线（或称节点）电压的幅值和 相位，所有线路的功率分布和功率损耗等。
3-2 电压降落是指元件首末端两点电压的相量差。
电压损耗是两点间电压绝对值之差。
当两点电压之间的相角差不大时， 可以近似地认为电压损耗等于电压降落的纵分量。
电压偏移是指网络中某点的实际电压同网络该处的额定电压之差。
电压 偏移可以用kV 表示，也可以用额定电压的百分数表示。
电压偏移=%100⨯-NNV V V 功率损耗包括电流通过元件的电阻和等值电抗时产生的功率损耗和电压 施加于元件的对地等值导纳时产生的损耗。
输电效率是是线路末端输出的有功功率2P 与线路首端输入的有功功率1P 之比。
输电效率=%10012⨯P P 3-3 网络元件的电压降落可以表示为()•••••+=+=-2221V V I jX R V V δ∆式中，•2V ∆和•2V δ分别称为电压降落的纵分量和横分量。
从电压降落的公式可见，不论从元件的哪一端计算，电压降落的纵、横分量计算公式的结构都是一样的，元件两端的电压幅值差主要有电压降落的纵分量决定，电压的相角差则由横分量决定。
在高压输电线路中，电抗要远远大于电阻，即R X 〉〉，作为极端的情况，令0=R ，便得V QX V /=∆，V PX V /=δ上式说明，在纯电抗元件中，电压降落的纵分量是因传送无功功率而产生的，而电压降落的横分量则是因为传送有功功率产生的。
换句话说，元件两端存在电压幅值差是传送无功功率的条件，存在电压相角差则是传送有功功率的条件。
3-4 求解已知首端电压和末端功率潮流计算问题的思路是，将该问题转化成已知同侧电压和功率的潮流计算问题。
首先假设所有未知点的节点电压均为额定电压，从线路末端开始，按照已知末端电压和末端潮流计算的方法，逐段向前计算功率损耗和功率分布，直至线路首端。