基于前推回代潮流计算(MATLAB)

1、本文主要运用前推回代潮流法计算菲宝线配电网络的线损，该网络图如下所示：p.号表3-3架空线路参数首节点尾节点导线型号长度(km)单位电阻9/km)单位电抗(QAm)12JKLGJ/Q-10kV-2400.0840.1250.125JKLGJ/Q-10kV-2400.0180.1250.167JKLGJ/Q-10kV-2400.0090.1250.159JKLGJ/Q-10kV-2400.1010.1250.1911JKLGJ/Q-10kV-2400.0310.1250.11112JKLGJ/Q-10kV-2400.0190.1250.11215JKLGJ/Q-10kV-2400.0480.12

2、50.11516JKLGJ/Q-10kV-2400.0480.1250.11617LJ2400.0540.1320.3051718LJ2400.0380.1320.305表3-4电缆线路参数首节点尾节点电缆型号长度(km)单位电阻(Q,/km)单位电抗(Q/km)01Y.JV22-&7/15KV-3*300MM22.0470.06010.089823YJV22-10KV-3*50MM20.1130.3870.11934YJV22-10KV-3*50MM20.1130.3870.11956YJV22-10KV-3*50MM20.0320.3870.119nonaorn007niin表3-5变压器

3、参数节点i4V5变压器名称型号联结组标号短路损耗(kW)短路电压百分比(幼空载损耗(kW)空载电流百分比僅)终端安装处。4日祥#1变S9-315/10Y,vnO3.6540.671.106菲宝织造配变S9-200/10Y,ynO2.640.48L3125大堡衲变S9-315/10Y,ynO3.6540.67L1124大堡村13#变S11-M-315/10Y,vnO3.6540.481.118建邦1#变S9-100/10Y,ynO1.540.29L6110凯欣1#变S9-315/10Y,ynO3.6540.67L1113宏升1#变S11-M-315/10Y,vnO3.6540.481.1114豪

4、宝14#变S9-M-315/10Y,ynO3.6540.671.1120贝瑞特服装5#变S9-315/10Y,ynO3.6540.67L1121贝瑞特服装6#变S9-315/10Y,ynO3.6540.671.1122艺彩布业配变S9-315/10Y,ynO3.6540.671.1123艺彩布业2#变S9-315/10Y,ynO3.6540.67L11a0表示配变抄表终端安装在变压器的高压侧，1代表在低压侧在本文中，主要的解答过程如下：一、节点和支路编号为了编程的方便，在此源点编号为1，其他节点随机编号，而支路的编号与该支路终点的编号取值一样。菲宝线配网的编号图如下所示：图二节点及支路编号二、

5、给节点分层在给节点分层之前，举例说明节点分层的方法。在图三中，先找出网络中的一颗树干（树干是随机的），同时去掉目标树干上的末梢节点。先将树干上的节点分层（树干上有n个节点就分成年n层），再找树干上各节点对应的子节点，将这些子节点分在同一层，依此类推，最后将末梢点分在同一层。在图一中，可将该网络分成六层，分层情况如下：第一层：1第二层：2第三层：3第四层：4第五层：6，8第六层：5，7，9，10在菲宝线潮流计算中，先建立原始数据矩阵Raw_data,第1列、第二列分别为支路始节点和末节点；第3-5列为该支路对应的线路参数，分别对应线路的长度、线路的单位电阻和单位阻抗（由于线路小于100km，线路

6、的导纳约等于零）；第5列表示支路的末节点是否装有变压器，若无置为0，若有且配网抄表终端安装在变压器的低压侧置为1，若配网抄表终端安装在变压器的高压侧置为2；第7-12列为对应变压器参数，分别为该变压器的额定容量、额定电压、短路损耗、短路电压百分比、空载损耗和空载电流百分比；第13-14列分别为该变压器所接负额的有功功率和无功功率。原始数据矩阵Raw_data，如图四所示田Raw_data123斗5678910111213141122.04700.06010.08980000000002230.08400.12500.400000000000033斗0.01800.12500.400000000

7、0000斗450.10100.12500.40000000000005560.03100.12500.40000000000006670.01900.12500.40000000000007780.04800.12500.40000000000008890.04800.12500.400000000000099100.05100.13200.30500000000001010110.04500.26800.11301315103.6500斗0.67001.100018.462519.50001111120001315103.6500斗0.67001.100055.67509.250012121

8、30001315103.6500斗0.67001.100014.340081313140001315103.6500斗0.67001.100010.41253143150.11300.38700.11900000000001515160.11300.38700.11902315103.6500斗0.67001.100044.025048.750016斗170.03200.38700.11901200102.6000斗0.48001.300020.59509.75001717180.00900.12500.40001315103.6500斗0.67001.100033.45302.7000181

9、8190001315103.6500斗0.48001.10009.09603.30001919200.08500.38700.11901100101.5000斗0.29001.60001.44381.9500205210.03700.38700.11901315103.6500斗0.67001.100012.38752.7500216220.05200.38700.11901315103.6500斗0.48001.100019.11256.5000227230.10300.15300.10201315103.6500斗0.67001.100021.04255.2500图四原始数据矩阵Raw_data图五节点分层矩阵layer三、计算过程第一步：给节点电压赋初值，计算出变压器的总注入功率。第二步：计算末梢点的注入电流，由外层向根节点推进，计算各节点的注入电流I(该电流就是支路电流)，也就是回推求解各支路电流。第三步：前推求解各节点电压U。U(c+1)=U(c+1)-I(c)Z第四步：求出各节点电压的修正量&U,同时计算出mXQU).。U(c+l)=U(c+1)-U