配电网潮流及线路损耗计算

1、解题思路

采用前推回代法求解。辐射状配电网的显著特征是从任意给定母线到源点有唯一路径，前推回代法就是充分利用了配电网络的这一特征，沿这些唯一的供电路径修正电压和电流。在前推过程中，首先根据配电检测终端量测的节点各相负荷的有功、无功功率计算支路电流；然后从各负荷支路开始向潮流的前方直到源点根据基尔霍夫电流定律计算各支路的电流分布；最后求出源节点流出的三相电流。在回代过程中，由已知电源电压和所求得的三相电流，从源节点向各负荷节点根据基尔霍夫电压定律计算系统所有节点的三相电压。每次迭代须对负荷电流做修正，即按求得的各负荷节点电压修正配电变压器的损耗，由修正后的变压器损耗和测量的负荷功率修正负荷支路电流。经过反复迭代和修正，直到两次迭代的各节点三相电压差均小于设定误差。最后，利用前推回代潮流计算收敛后的电流、电压值和等效线路的阻抗、配电变压器的阻抗和导纳值算出相应的各支路损耗和各变压器损耗值。

2、程序流程

根据基于节点(支路)分层的广度优先搜索编号法对网络进行节点和支路编号；

设置电压初值，本程序电压初值设为线电压10kV；

计算网络三相等效参数；

计算负荷电流；

前推计算支路电流；

回代计算节点电压；

收敛判定。如果收敛则进行下一步，否则把计算出的各节点电压值替代迭代前的各节点电压值代入步骤4)重新计算负荷电流；

计算支路损耗；

计算变压器损耗；

计算馈线总损耗

3、程序

clear;clc

format long

Z=[0,1,2.047*0.0601+i*2.047*0.0898,0,0,0,0,0,0,0

1,2,0.084*0.125+i*0.084*0.4,0,0,0,0,0,0,0

2,3,0.018*0.125+i*0.018*0.4,0,0,0,0,0,0,0

2,4,0.113*0.387+i*0.113*0.119,0,0,0,0,0,0,0

3,5,0.101*0.125+i*0.101*0.4,0,0,0,0,0,0,0

3,6,0.032*0.387+i*0.032*0.119,2.6*10^2*10^3/200^2,4*10^2*10/200,0.48*10^-3/10^2,1.3*200/10^2*10^-5,20.595+i*9.75,26.1225+i*10.2,19.9575+i*8.4

4,7,0.113*0.387+i*0.113*0.119,0,0,0,0,44.025+i*48.75,0+i*0,67.95+i*7.5

5,8,0.031*0.125+i*0.031*0.4,0,0,0,0,0,0,0

5,9,0.037*0.387+i*0.037*0.119,3.65*15^2*10^3/315^2,4*15^2*10/315,0.67*10^-3/15^2,1.1*315/15^2*10^-5,12.3875+i*2.75,11.6375+i*0.75,17.775+i*2.75

6,10,0.009*0.125+i*0.009*0.4,0,0,0,0,0,0,0

8,11,0.019*0.125+i*0.019*0.4,0,0,0,0,0,0,0

8,12,0.052*0.387+i*0.052*0.119,3.65*10^2*10^3/315^2,4*10^2*10/315,0.48*10^-3/10^2,1.1*315/10^2*10^-5,19.1125+i*6.5,19.0625+i*6,18.4+i*4.75

10,13,0.085*0.387+i*0.085*0.119,1.5*10^2*10^3/100^2,4*10^2*10/100,0.29*10^-3/10^2,1.6*100/10^2*10^-5,1.44375+i*1.95,0+i*0,0+i*0

10,14,0,3.65*10^2*10^3/315^2,4*10^2*10/315,0.48*10^-3/10^2,1.1*315/10^2*10^-5,9.096+i*3.3,7.311+i*1.8,15.234+i*6.6

10,15,0,3.65*10^2*10^3/315^2,4*10^2*10/315,0.67*10^-3/10^2,1.1*315/10^2*10^-5,33.453+i*2.7,24.675+i*3.6,25.206+i*2.7

11,16,0.048*0.125+i*0.048*0.4,0,0,0,0,0,0,0

11,17,0.103*0.153+i*0.103*0.102,3.65*10^2*10^3/315^2,4*10^2*10/315,0.67*10