电力系统分析潮流计算讲解

1、电力系统分析潮流计算报告寸录配电网概述3 1配电网的分类3. 3.L2配电网运行的特点及要求L3配电网潮流计算的意义-计算原理及计算流程2.1前推回代法计算原理2. 2前推回代法计算流程2. 3主程序清单:9.2. 4输入文件清单:4.1. 12. 5计算结果清单：1.2.三,前推回代法计算流程图3参考文献配电网才既述1配电网的分类在电力网中重要起分配电能作用的网络就称为配电网；配电网按电压等级来分类，可分为高压配电网(35- UOKV),中压配电网(610KV,苏州有20KV的),低压配电网(220/380V )；在负载率较大的特大型城市，220KV电网也有配电功能。按供电区的功能来分类，可

2、分为城市配电网，农村配电网和工厂配电网等。在城市电网系统中，主网是指110KV及其以上电压等级的电网，主要 起连接 区域高压(220KV及以上)电网的作用。配电网是指35KV及其以下电压等级的电网，作用是给城市里各个配电站和各类用电负荷供给电源。从投资角度看，我国与国外先进国家的发电、输电、配电投资比率差异很刚从重发大，国外基本上是电网投资大于电厂投资，输电投资小于配电投资。我 电轻供电状态中转变过来，而在供电投资中，输电投资大于配电投资。从我国城网 改造之后，将逐渐从输电投资转入配电建设为主。本文是基于前推回代法的配电网潮流分析计算的研究，研究是是以根节点为10kV的电压等级的配电网。配电网

3、运行的特点及要求配电系统相对于输电系统来说，由于电压等级低、供电范围小，但与用户直接 相连，是供电部门对用户服务的窗口，因而决定了配电网运行有如下特点和基 本要 求：10kv中压配电网在运行中，负荷节点数多，一般无表计实时记录负荷，无法应用现在传统潮流程序进行配电网的计算分析，要求建立新的数学 模型和计算方法。（2）随着铁道电气化和用户电子设备的大量使用，配电网运行中有大量的 谐波源、三相电压不平衡、电压闪变等污染，要求准确测量与计算配电网中的 谐波 分布，从而采取有效措施抑制配电网运行中的谐波危害。（3）由于环保条件日趋严格的制约，要求配电网运行能制定不影响城市绿 化、防火、防爆、防噪音等技

4、术和组织措施，以便减少配电网运行对环境的污 染。（4）随着用户对供电可靠性和电压质量指标的提高，还靠人工操作已无 法适应，要求现代配电网运行不断提高自动化、智能水平。由于“电能”作为商品将进入市场竞争，要求各电力公司采用需求侧管理和用 户电力技术，以降低配电网运行的线损和年运行费用，提高运行的经济性，从 而降 低配电成本，并积极协助用户搞好优化用电计划、节约用电，推行战略节电 和战略 负荷开拓等积极措施，进一步提高对用户的服务质量和降低供电企业的成 本，达到 双方受益的目的。配电网潮流计算的意义配电网潮流计算是配电网网络分析的基础，配电网的网络重构、故障处 理、无功优化和状态估计等都需要用到配

5、电网的潮流数据。由于配电网结构特点 都 是开环运行的，配电网呈辐射状，配电线路的电阻电抗比（R/X）较大，利用常规方法进行潮流计算会导致算法不收敛，而前推回代法是线性收敛的，解决了潮流计算收敛难的问题。近年来,开发配电管理系统（DMS）成为人们研究的热点。而配电网潮流计算作为DMS的高级应用软件之一，更是整个问题研究和分析的基础。二计算原理及计算流程2. 1前推回代法计算原理前推回代法已知配电网的始端电压和末端负荷，以馈线为基本计算单位。最初假 设全网电压都为额定电压，根据负荷功率由末端向始端逐段推算，仅计算各元件中的功率损耗而不计算节点电压，求得各支路上的电流和功率损 耗，并据此获得始端功率

6、，这是回代过程；再根据给定的始端电压和求得的始端功 率，由始端向末端逐段推算电压降落，求得各节点电压，这是前推过程。如 此重复 上述过程，直至各个节点的功率偏差满足允许条件为o图2-1所示的网络结构即为典型的辐射状配电网结构。图21辐射状配电系统首先要搜索节点矢系，确定拓扑结构表。为了配合算法和避免复杂的网络编号，采用以下原始数据输入结构，不用形成节点导纳矩阵，就可以H动搜索 节点矢系，确定网络的拓扑结构。节点结构体：节点号节点有功节点无功支路结构体：支路首端节点号支路末端节点号支路电阻支路电抗根据线路首末 节点，就可以确定每个节点连接的节点及其矣系，从而可以形成整体的树状的矢系 结构。为了形

7、成层次矢系，确定节点计算顺序，要利用网络拓扑结构，经过多次按 层遍历的广度优先搜索，形成层次尖系，确定前推后代潮流算法的节点计算顺序。 具体方法如下：搜索末梢节点作为第一层节点；搜索末梢节点的父节点作为第二层节继续搜索第二层节点的父节点作为第三层节点，这样反的搜索下去，直到搜索到某层节点的父节点全部是根节点时停止搜索；删除在后面层次中有重复的前面层次中的节点，形成真正的层次尖系，确定潮流计算的节点顺序0前推回代法基于支路电流进行，首先假定各节点的电压幅值为1,幅角为0,具 体计算步骤为：1)从第一层节点开始，根据基尔霍夫电流定律，求支路上的电流：li 厂 Sj/Uj 二 Pj-jQj/Uj2-

8、)SU式中是节点j的功率，w是节点j的电压。2)从第二层开始逐层计算非末梢节点的注入电流，根据基尔霍夫电流定律应 等于(34 )式与该节点流出电流之和：111 lij(k)二 Sj/Uj(k) + ?Ijk(k)-(2-2)3)由步骤1)和2)可求出所有支路的支电流，再利用已知的根节点电压，从根节点向后顺次求得各个负荷节点的电压 其中i为父节点为子节点，乙为i、Zij4)计算各节点的电压幅值修正量:DUj(l)二 uj(i)- u5)计算节点电压修正量的最大j间支删脾抗。Uj二Ui(i)- 1-(0)2-4)max (dUj(k)(S-3)r a /习、n 仃 XRUI其中k为迭代次数，若最大

9、电压修正量小于阈值，则跳出循环，输出电压计算结果；否则重复步骤（1） s （6）,直到满足（35 ）式的条件为7）在得到各个节点的电压电流后，就可以计算线路潮流S:5uJ a2-6）总之，（2-1 ）（2-3 ）式为一组递归方程，对树进行前向遍历，从树的末梢节点出发，利用已知的负荷功率，逐一计算（2-1 ）（2-2）式，即可求得根节点处的电流，再从根节点出发，对树进行后向遍历，用（2-3）式可 求各 节点电压。这样就完成一次前推后代的计算，迭代重复进行，直至满足收敛标准为 止。2.2前推回代法计算流程要看懂前推回代法计算程序，报告叙述计算原理及计算流程。绘制计算流程 框图。确定前推回代支路次序

10、（广度优先，或深度优先），编写前推回代计算输 入文件。进行潮流计算。下列为节点配电网结构图及系统支路参数和系统负荷参数表。6210 11 1213!节点配电网结构图表1系统支路参数支路R(Q-)X(Q.)B/2 (T K) (S)23.3673.6850.0232. 3562.5410.03s41. 1451.280.04. 5245.040.02R0. 8561. 140.0672. 7452. 9650.0283. 7434. 2510.0Q2.2372. 7560.0qm4. 144. 6960.02 111.3281. 7630.01T122. 4362. 8660.04 133. 5

11、213.9660.0表2系统负荷参数节点编号节点类型节点初始电压(kV)Pi ( MV A)Qi (MVA)1根节点10.4002中间节点10.00.03420. 03013中间节点10.00.06930. 06424中间节点10.00.08450. 07635叶节点10.00.02950. 02616中间节点10.00. 04740. 04097叶节点10.00. 11760. 09578中间节点10.00.09460. 08579叶节点10.00.09160. 085910叶节点10.00. 02710. 022911中间节点10.00.06960. 064312叶节点10.00. 067

12、60. 057913叶节点10. 00.02980. 02423主程序清单：PQ, FT, RX二 easel 13() ；%调用数据文件NN 二 size(PQ, 1)；% 节点数NB 二 size (FT, 1) ；%支路数数V二PQ(： ,1)；%v初始电压相量maxd=lk=lwhile maxd0. 0001PQ2二PQ ；%每一次迭代各节点的注入有功和无功相同PL 二 0. 0：for i=l： NBkf二FT(i, 1)；%前推始节点号kt二FT(i,2) ；%前推终节点号x 二(PQ2(kf,2 厂 2+PQ2(kf,3 厂 2)/V (kf)/V (kf) ； %计算沿线电流

13、平方 APQl(i, l)Z2PQ2(kf,2)+RX(i, l)*x； % 计算支路首端有功 /MW为计算沿支路的无功损耗/MvhRX(i, 1 厂 RPQ1 (i, 2)=PQ2 (kf, 3)+RX(i, 2)*x；RX(i,2)XPQ2(kt,2尸PQ2(kt,2)+PQl(i, 1) ； %用PQ1去修正支路末端节点的 有功P单位MW PQ2 (kt,3)= PQ2(kt,3)+PQl(i,2) ； %用PQ1去修正支路末端节点的有功Q单位MvarPL 二 PL+RX(i, l)*x；endangle。尸0. 0；fbii 二1kf二FT(i,2)；%回代始节点号毗二FT (i, 1

14、)；%回代终节点号%计算支路电压损耗的纵%计算支路电压损耗的横0/二，咎* Ttr 口士占草社3 m dv (PQ1 (i, 1)*RX (i, 1) +PQ1 (i, 2)*RX(i, 2) )/V (kf);分量dvldv2 Zl(PQl(i, l)*RX(i,2)-PQl(i,2)*RX(i, l)A(kf)；分量 dv2V2 (kt) =sqil(V(kf)-dvl) 2+dv2 2)；angle(kt)=angle(kf)+atand(dv2/(V(kf)-dvl)； end maxd=abs(V2(2)-V (2)；V2(l)二 V(l)；fbii二 3： 1： NNif abs

15、(V2 (i)-V(i)maxd；maxd=abs(V2(i)-V (i)；endendmaxd k=k+lV=V2%节点电压模计算结果kvanglePL%节点电压角度计算结果单位 外网损单位MWend clear24 输入文件 j 青单：function PQ, FT, RX=casell3 () PQ=%节点电压有功无功10.4 0010.0 0. 0342 0. 030110.0 0. 0693 0. 064210.0 0. 0845 0. 076310.0 0. 0295 0. 026110.0 0. 0474 0. 040910.0 0. 1176 0. 095710.0 0. 09

16、46 0. 085710.0 0. 0916 0. 085910.0 0. 0271 0. 022910.0 0. 0696 0. 064310.0 0. 0676 0. 057910.0 0. 0298 0. 0242FT=%首端末端10 35412 112；RX=%R4. 144. 5243. 5211452.4361.3287450. 8562.2377432.3563.367):X696043.9661.282. 8661.7632. 965147564.2512. 541 3.6852. 5计算结果清单maxd = 1k 二 1maxd =0.1780k=2maxd =0.1787k

17、=3maxd =0.1793k =4maxd =0.1226k=5maxd =0.0010k=6maxd =2.602 le-004 k =7maxd =6.1046e-005k=8PQ1 =0. 02720. 02960. 02300. 02620. 02990. 14440. 02430. 12730. 06780. 13790. 11830. 16610. 09200. 18910.05810. 12310. 09640.13790. 08640. 17500. 38520. 34450.8116V=10 - 40000.72789. 8795 9. 699/ 9. 6652 9. 6377 9. 8492 9. 7872 9. 7326 9. 6870 9. 6764 9. 6578 9. 6235 9. 6444angle =00. 3011 0. 4011 0. 4250 0. 4438 0. 3431 0. 3942 0. 38990. 4266 0. 4209 0. 4498 0. 4823 0. 4452PL =0.0488三前推回代法计算流程图经过之前的理论知识，我们对MATLAE的基本指令有了初步的了解，我们就可以 画出程序编程框图，按照程序编程框图，我们可以通过MATLAE软件进行编程，实现 我们所求的前推回代法计算。程序