NSW Electricity Network

From MoxyWiki

Contents

General Description

This is a general instance of the Electricity Network. It represents the electricity network within New South Wales and has been taken from the File:Line Diagram of the NSW Energy Network.pdf

Properties of NSW Electricity Network

The distributed efficiency of the NSW Electricity Network is 0.96info.png
 0.96

Last calculated 72 weeks ago
. This refers to the efficiency of distributed generation using the electricity transmission and distribution network. [1]. The centralized efficiency of the NSW Electricity Network is 0.945info.png
 0.945

Last calculated 72 weeks ago
[2]. This refers to the efficiency of centralized generation using the electricity transmission and distribution network.[1]. The inverse of this number is conventionally referred to as the Loss Factor, and is used for converting from the electricity consumed in the state and the electricity generated. The initial cost of the hardware required to support additional distributed solar capacity, the transmission and distribution cost per Watt for the NSW Electricity Network, for every Watt of additional installed distributed generation the cost is AU$1.30 AUD 2009/wattinfo.png
'1.30 AUD 2009/watt'

Last calculated 72 weeks ago
per Watt. [1]. The Operation cost as a percentage of the capital cost for the NSW Electricity Market is $0.02info.png
 0.02

Last calculated 72 weeks ago
. This is the percentage of the total installed cost of the Transmission and Distribution equipment installed specifically to cater for the new solar capacity that will have to be paid in operation and maintenance costs per year. [1]. The replacement period for the new transmission and distribution equipment installed to support the additional solar capacity is estimated to need to be replaced after approximately 25 yrinfo.png
'25 yr'

Last calculated 72 weeks ago
. [1]. The average emissions per MWh of the NSW Electricity Network is 0.906t/MWhinfo.png
'0.906t/MWh'

Last calculated 72 weeks ago
tonnes of CO2e/MWh [2]. This value reflects the amount of CO2 equivalent emissions that are generated in the conventional, centralized NSW power generation network for every MWh of generated power. Note that this is for generated power, not consumed power.

Robustness

The Robustness of NSW Electricity Grid has been analysed by comparing a number of different outage scenarios. The Load Flow analysis gives an indication of how the network will respond under different conditions.

Network

The NSW Electricity Network is made up of the following buses.

Bus Id RealPowerDemand ReactivPowerDemand ShuntSusceptance VoltageMagnitude RealPowerOutput MaximumReactivePowerOutput MinimumReactivePowerOutput VoltageMagnitudeSetpoint
Armidale Electrical Bus 2 230 MWinfo.png
'230 MW'

Last calculated 65 weeks ago
10 MVAinfo.png
'10 MVA'

Last calculated 65 weeks ago
-0.623 MVAinfo.png
'-0.623 MVA'

Last calculated 65 weeks ago
Avon Electrical Bus 29 0 MWinfo.png
'0 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Bayswater Powerstation 9 0 MWinfo.png
'0 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
1.057info.png
1.057

Last calculated 65 weeks ago
0 MWinfo.png
'0 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
1.057info.png
1.057

Last calculated 65 weeks ago
Beaconsfield West Electrical Bus 27 590 MWinfo.png
'590 MW'

Last calculated 65 weeks ago
170 MVAinfo.png
'170 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Canberra Electrical Bus 33 0 MWinfo.png
'0 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
1info.png
1

Last calculated 73 weeks ago
230 MWinfo.png
'230 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
1info.png
1

Last calculated 73 weeks ago
Dapto Electrical Bus 28 460 MWinfo.png
'460 MW'

Last calculated 65 weeks ago
40 MVAinfo.png
'40 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Dederang Electrical Bus 40 940 MWinfo.png
'940 MW'

Last calculated 65 weeks ago
-350 MVAinfo.png
'-350 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Dumaresq Substation 1 0 MWinfo.png
'0 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
1.035info.png
1.035

Last calculated 65 weeks ago
1000 MWinfo.png
'1000 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
-300 MVAinfo.png
'-300 MVA'

Last calculated 65 weeks ago
1.035info.png
1.035

Last calculated 65 weeks ago
Eraring Powerstation 10 190 MWinfo.png
'190 MW'

Last calculated 65 weeks ago
-142 MVAinfo.png
'-142 MVA'

Last calculated 65 weeks ago
1.03info.png
1.03

Last calculated 65 weeks ago
1200 MWinfo.png
'1200 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
1.03info.png
1.03

Last calculated 65 weeks ago
Ingleburn Electrical Bus 23 113 MWinfo.png
'113 MW'

Last calculated 65 weeks ago
71 MVAinfo.png
'71 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Jindera Electrical Bus 38 122 MWinfo.png
'122 MW'

Last calculated 65 weeks ago
86 MVAinfo.png
'86 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Kangaroo Valley Electrical Bus 32 521 MWinfo.png
'521 MW'

Last calculated 65 weeks ago
270 MVAinfo.png
'270 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Kemps Creek Electrical Bus 11 0 MWinfo.png
'0 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
1info.png
1

Last calculated 73 weeks ago
1390 MWinfo.png
'1390 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
1info.png
1

Last calculated 73 weeks ago
Liddell Powerstation 6 40 MWinfo.png
'40 MW'

Last calculated 65 weeks ago
90 MVAinfo.png
'90 MVA'

Last calculated 65 weeks ago
1.05info.png
1.05

Last calculated 65 weeks ago
1850 MWinfo.png
'1850 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
1.05info.png
1.05

Last calculated 65 weeks ago
Lismore Electrical Bus 3 140 MWinfo.png
'140 MW'

Last calculated 65 weeks ago
36 MVAinfo.png
'36 MVA'

Last calculated 65 weeks ago
-0.044 MVAinfo.png
'-0.044 MVA'

Last calculated 65 weeks ago
Liverpool Electrical Bus 24 240 MWinfo.png
'240 MW'

Last calculated 65 weeks ago
38 MVAinfo.png
'38 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Lower Tumut Electrical Bus 34 2 MWinfo.png
'2 MW'

Last calculated 65 weeks ago
-300 MVAinfo.png
'-300 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Marulan Electrical Bus 30 82 MWinfo.png
'82 MW'

Last calculated 65 weeks ago
20 MVAinfo.png
'20 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Mount Piper Powerstation 20 120 MWinfo.png
'120 MW'

Last calculated 65 weeks ago
16 MVAinfo.png
'16 MVA'

Last calculated 65 weeks ago
1.03info.png
1.03

Last calculated 65 weeks ago
570 MWinfo.png
'570 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
1.03info.png
1.03

Last calculated 65 weeks ago
Munmorah Powerstation 13 21 MWinfo.png
'21 MW'

Last calculated 65 weeks ago
46 MVAinfo.png
'46 MVA'

Last calculated 65 weeks ago
1.02info.png
1.02

Last calculated 65 weeks ago
470 MWinfo.png
'470 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
1.02info.png
1.02

Last calculated 65 weeks ago
Murray Electrical Bus 35 0 MWinfo.png
'0 MW'

Last calculated 65 weeks ago
-13.5 MVAinfo.png
'-13.5 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Muswellbrook Electrical Bus 5 20 MWinfo.png
'20 MW'

Last calculated 65 weeks ago
120 MVAinfo.png
'120 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Newcastle Electrical Bus 7 1000 MWinfo.png
'1000 MW'

Last calculated 65 weeks ago
290 MVAinfo.png
'290 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Regentville Electrical Bus 19 150 MWinfo.png
'150 MW'

Last calculated 65 weeks ago
90 MVAinfo.png
'90 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Sterland Electrical Bus 14 0 MWinfo.png
'0 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
1info.png
1

Last calculated 73 weeks ago
Sydney East Electrical Bus 17 570 MWinfo.png
'570 MW'

Last calculated 65 weeks ago
50 MVAinfo.png
'50 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Sydney North Electrical Bus 16 1070 MWinfo.png
'1070 MW'

Last calculated 65 weeks ago
-226 MVAinfo.png
'-226 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Sydney South Electrical Bus 26 903 MWinfo.png
'903 MW'

Last calculated 65 weeks ago
-61 MVAinfo.png
'-61 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Sydney West Electrical Bus 18 1070 MWinfo.png
'1070 MW'

Last calculated 65 weeks ago
210 MVAinfo.png
'210 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Tamworth Electrical Bus 4 130 MWinfo.png
'130 MW'

Last calculated 65 weeks ago
17 MVAinfo.png
'17 MVA'

Last calculated 65 weeks ago
0.817 MVAinfo.png
'0.817 MVA'

Last calculated 65 weeks ago
Tomago Aluminium Smelter 8 740 MWinfo.png
'740 MW'

Last calculated 65 weeks ago
280 MVAinfo.png
'280 MVA'

Last calculated 65 weeks ago
1.37 MVAinfo.png
'1.37 MVA'

Last calculated 65 weeks ago
Tuggerah Electrical Bus 15 230 MWinfo.png
'230 MW'

Last calculated 65 weeks ago
80 MVAinfo.png
'80 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Upper Tumut Electrical Bus 36 0 MWinfo.png
'0 MW'

Last calculated 65 weeks ago
-14.4 MVAinfo.png
'-14.4 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Vales Point Powerstation 12 33 MWinfo.png
'33 MW'

Last calculated 65 weeks ago
17 MVAinfo.png
'17 MVA'

Last calculated 65 weeks ago
1.02info.png
1.02

Last calculated 65 weeks ago
1200 MWinfo.png
'1200 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
1.02info.png
1.02

Last calculated 65 weeks ago
Vineyard Electrical Bus 25 103 MWinfo.png
'103 MW'

Last calculated 65 weeks ago
69 MVAinfo.png
'69 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Wagga Wagga Electrical Bus 37 250 MWinfo.png
'250 MW'

Last calculated 65 weeks ago
-131.7 MVAinfo.png
'-131.7 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Wallerawang Powerstation 22 150 MWinfo.png
'150 MW'

Last calculated 65 weeks ago
18 MVAinfo.png
'18 MVA'

Last calculated 65 weeks ago
1.03info.png
1.03

Last calculated 65 weeks ago
890 MWinfo.png
'890 MW'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
1.03info.png
1.03

Last calculated 65 weeks ago
Wellington Electrical Bus 21 200 MWinfo.png
'200 MW'

Last calculated 72 weeks ago
35 MVAinfo.png
'35 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Wodonga Electrical Bus 39 40 MWinfo.png
'40 MW'

Last calculated 65 weeks ago
36 MVAinfo.png
'36 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Yass Electrical Bus 31 216 MWinfo.png
'216 MW'

Last calculated 65 weeks ago
16 MVAinfo.png
'16 MVA'

Last calculated 65 weeks ago
0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago


and the following transmission lines.

Line ShortName FromBus FromId ToBus ToId Resistance Reactance TotalLineChargingSusceptance
01X Electrical Transmission Line '01X' Kangaroo Valley Electrical Bus 32 Yass Electrical Bus 31 0.001692info.png
0.001692

Last calculated 65 weeks ago
0.01301info.png
0.01301

Last calculated 65 weeks ago
0.1595info.png
0.1595

Last calculated 65 weeks ago
01 Electrical Transmission Line '01' Kangaroo Valley Electrical Bus 32 Upper Tumut Electrical Bus 36 0.003818info.png
0.003818

Last calculated 65 weeks ago
0.03245info.png
0.03245

Last calculated 65 weeks ago
0.3617info.png
0.3617

Last calculated 65 weeks ago
02 Electrical Transmission Line '02' Yass Electrical Bus 31 Upper Tumut Electrical Bus 36 0.005888info.png
0.005888

Last calculated 65 weeks ago
0.04686info.png
0.04686

Last calculated 65 weeks ago
0.5542info.png
0.5542

Last calculated 65 weeks ago
03 Electrical Transmission Line '03' Yass Electrical Bus 31 Lower Tumut Electrical Bus 34 0.005166info.png
0.005166

Last calculated 65 weeks ago
0.03951info.png
0.03951

Last calculated 65 weeks ago
0.4901info.png
0.4901

Last calculated 65 weeks ago
051 Electrical Transmission Line '051' Lower Tumut Electrical Bus 34 Wagga Wagga Electrical Bus 37 0.004022info.png
0.004022

Last calculated 65 weeks ago
0.03071info.png
0.03071

Last calculated 65 weeks ago
0.3822info.png
0.3822

Last calculated 65 weeks ago
07 Electrical Transmission Line '07' Kangaroo Valley Electrical Bus 32 Lower Tumut Electrical Bus 34 0.003972info.png
0.003972

Last calculated 65 weeks ago
0.03046info.png
0.03046

Last calculated 65 weeks ago
0.3735info.png
0.3735

Last calculated 65 weeks ago
10 Electrical Transmission Line '10' Avon Electrical Bus 29 Dapto Electrical Bus 28 0.000311info.png
0.000311

Last calculated 65 weeks ago
0.003298info.png
0.003298

Last calculated 65 weeks ago
0.04159info.png
0.04159

Last calculated 65 weeks ago
11 Electrical Transmission Line '11' Sydney South Electrical Bus 26 Dapto Electrical Bus 28 0.00195info.png
0.00195

Last calculated 65 weeks ago
0.02061info.png
0.02061

Last calculated 65 weeks ago
0.2602info.png
0.2602

Last calculated 65 weeks ago
12 Electrical Transmission Line '12' Liverpool Electrical Bus 24 Sydney South Electrical Bus 26 0.000522info.png
0.000522

Last calculated 65 weeks ago
0.005267info.png
0.005267

Last calculated 65 weeks ago
0.07325info.png
0.07325

Last calculated 65 weeks ago
13 Electrical Transmission Line '13' Kemps Creek Electrical Bus 11 Sydney South Electrical Bus 26 0.00069info.png
0.00069

Last calculated 65 weeks ago
0.006991info.png
0.006991

Last calculated 65 weeks ago
0.09616info.png
0.09616

Last calculated 65 weeks ago
14 Electrical Transmission Line '14' Sydney North Electrical Bus 16 Kemps Creek Electrical Bus 11 0.001431info.png
0.001431

Last calculated 65 weeks ago
0.01504info.png
0.01504

Last calculated 65 weeks ago
0.1914info.png
0.1914

Last calculated 65 weeks ago
16 Electrical Transmission Line '16' Avon Electrical Bus 29 Marulan Electrical Bus 30 0.002838info.png
0.002838

Last calculated 65 weeks ago
0.02178info.png
0.02178

Last calculated 65 weeks ago
0.2673info.png
0.2673

Last calculated 65 weeks ago
18 Electrical Transmission Line '18' Canberra Electrical Bus 33 Dapto Electrical Bus 28 0.001713info.png
0.001713

Last calculated 65 weeks ago
0.01317info.png
0.01317

Last calculated 65 weeks ago
0.1615info.png
0.1615

Last calculated 65 weeks ago
20 Electrical Transmission Line '20' Sydney North Electrical Bus 16 Sydney West Electrical Bus 18 0.000954info.png
0.000954

Last calculated 65 weeks ago
0.01007info.png
0.01007

Last calculated 65 weeks ago
0.1273info.png
0.1273

Last calculated 65 weeks ago
21A Electrical Transmission Line '21A' Munmorah Powerstation 13 Sterland Electrical Bus 14 0.000759info.png
0.000759

Last calculated 65 weeks ago
0.008023info.png
0.008023

Last calculated 65 weeks ago
0.1012info.png
0.1012

Last calculated 65 weeks ago
21B Electrical Transmission Line '21B' Sterland Electrical Bus 14 Tuggerah Electrical Bus 15 0.000398info.png
0.000398

Last calculated 65 weeks ago
0.004073info.png
0.004073

Last calculated 65 weeks ago
0.05128info.png
0.05128

Last calculated 65 weeks ago
21C Electrical Transmission Line '21C' Sterland Electrical Bus 14 Sydney North Electrical Bus 16 0.001462info.png
0.001462

Last calculated 65 weeks ago
0.01549info.png
0.01549

Last calculated 65 weeks ago
0.1954info.png
0.1954

Last calculated 65 weeks ago
22 Electrical Transmission Line '22' Vales Point Powerstation 12 Sydney North Electrical Bus 16 0.002463info.png
0.002463

Last calculated 65 weeks ago
0.02611info.png
0.02611

Last calculated 65 weeks ago
0.3292info.png
0.3292

Last calculated 65 weeks ago
23 Electrical Transmission Line '23' Vales Point Powerstation 12 Munmorah Powerstation 13 0.000208info.png
0.000208

Last calculated 65 weeks ago
0.002174info.png
0.002174

Last calculated 65 weeks ago
0.0274info.png
0.0274

Last calculated 65 weeks ago
24 Electrical Transmission Line '24' Newcastle Electrical Bus 7 Vales Point Powerstation 12 0.00106info.png
0.00106

Last calculated 65 weeks ago
0.01013info.png
0.01013

Last calculated 65 weeks ago
0.148info.png
0.148

Last calculated 65 weeks ago
25 Electrical Transmission Line '25' Eraring Powerstation 10 Vineyard Electrical Bus 25 0.003147info.png
0.003147

Last calculated 65 weeks ago
0.03074info.png
0.03074

Last calculated 65 weeks ago
0.4496info.png
0.4496

Last calculated 65 weeks ago
26 Electrical Transmission Line '26' Munmorah Powerstation 13 Sydney West Electrical Bus 18 0.003555info.png
0.003555

Last calculated 65 weeks ago
0.035info.png
0.035

Last calculated 65 weeks ago
0.5077info.png
0.5077

Last calculated 65 weeks ago
27 Electrical Transmission Line '27' Sydney North Electrical Bus 16 Sydney East Electrical Bus 17 0.000631info.png
0.000631

Last calculated 65 weeks ago
0.006431info.png
0.006431

Last calculated 65 weeks ago
0.08611info.png
0.08611

Last calculated 65 weeks ago
28 Electrical Transmission Line '28' Sydney North Electrical Bus 16 Sydney East Electrical Bus 17 0.000647info.png
0.000647

Last calculated 65 weeks ago
0.006769info.png
0.006769

Last calculated 65 weeks ago
0.08545info.png
0.08545

Last calculated 65 weeks ago
29 Electrical Transmission Line '29' Sydney West Electrical Bus 18 Vineyard Electrical Bus 25 0.000591info.png
0.000591

Last calculated 65 weeks ago
0.005796info.png
0.005796

Last calculated 65 weeks ago
0.08479info.png
0.08479

Last calculated 65 weeks ago
30 Electrical Transmission Line '30' Sydney West Electrical Bus 18 Liverpool Electrical Bus 24 0.00048info.png
0.00048

Last calculated 65 weeks ago
0.005035info.png
0.005035

Last calculated 65 weeks ago
0.06366info.png
0.06366

Last calculated 65 weeks ago
31 Electrical Transmission Line '31' Bayswater Powerstation 9 Regentville Electrical Bus 19 0.004892info.png
0.004892

Last calculated 65 weeks ago
0.04799info.png
0.04799

Last calculated 65 weeks ago
0.7021info.png
0.7021

Last calculated 65 weeks ago
32 Electrical Transmission Line '32' Bayswater Powerstation 9 Sydney West Electrical Bus 18 0.00541info.png
0.00541

Last calculated 65 weeks ago
0.05293info.png
0.05293

Last calculated 65 weeks ago
0.7742info.png
0.7742

Last calculated 65 weeks ago
33 Electrical Transmission Line '33' Liddell Powerstation 6 Bayswater Powerstation 9 0.000176info.png
0.000176

Last calculated 65 weeks ago
0.001696info.png
0.001696

Last calculated 65 weeks ago
0.02482info.png
0.02482

Last calculated 65 weeks ago
34 Electrical Transmission Line '34' Liddell Powerstation 6 Bayswater Powerstation 9 0.000176info.png
0.000176

Last calculated 65 weeks ago
0.001696info.png
0.001696

Last calculated 65 weeks ago
0.02482info.png
0.02482

Last calculated 65 weeks ago
35 Electrical Transmission Line '35' Mount Piper Powerstation 20 Marulan Electrical Bus 30 0.002841info.png
0.002841

Last calculated 65 weeks ago
0.0375info.png
0.0375

Last calculated 65 weeks ago
0.8098info.png
0.8098

Last calculated 65 weeks ago
36 Electrical Transmission Line '36' Mount Piper Powerstation 20 Marulan Electrical Bus 30 0.00284info.png
0.00284

Last calculated 65 weeks ago
0.0375info.png
0.0375

Last calculated 65 weeks ago
0.8097info.png
0.8097

Last calculated 65 weeks ago
37 Electrical Transmission Line '37' Avon Electrical Bus 29 Kemps Creek Electrical Bus 11 0.001821info.png
0.001821

Last calculated 65 weeks ago
0.01929info.png
0.01929

Last calculated 65 weeks ago
0.2521info.png
0.2521

Last calculated 65 weeks ago
38 Electrical Transmission Line '38' Sydney West Electrical Bus 18 Regentville Electrical Bus 19 0.000501info.png
0.000501

Last calculated 65 weeks ago
0.004912info.png
0.004912

Last calculated 65 weeks ago
0.07186info.png
0.07186

Last calculated 65 weeks ago
39 Electrical Transmission Line '39' Sydney West Electrical Bus 18 Yass Electrical Bus 31 0.009566info.png
0.009566

Last calculated 65 weeks ago
0.07349info.png
0.07349

Last calculated 65 weeks ago
0.9011info.png
0.9011

Last calculated 65 weeks ago
41 Electrical Transmission Line '41' Sydney South Electrical Bus 26 Beaconsfield West Electrical Bus 27 0.000396info.png
0.000396

Last calculated 65 weeks ago
0.003412info.png
0.003412

Last calculated 65 weeks ago
2.479info.png
2.479

Last calculated 65 weeks ago
4 Electrical Transmission Line '4' Marulan Electrical Bus 30 Yass Electrical Bus 31 0.00458info.png
0.00458

Last calculated 65 weeks ago
0.0351info.png
0.0351

Last calculated 65 weeks ago
0.4311info.png
0.4311

Last calculated 65 weeks ago
5 Electrical Transmission Line '5' Marulan Electrical Bus 30 Yass Electrical Bus 31 0.004735info.png
0.004735

Last calculated 65 weeks ago
0.03642info.png
0.03642

Last calculated 65 weeks ago
0.4465info.png
0.4465

Last calculated 65 weeks ago
62 Electrical Transmission Line '62' Wagga Wagga Electrical Bus 37 Jindera Electrical Bus 38 0.003996info.png
0.003996

Last calculated 65 weeks ago
0.03084info.png
0.03084

Last calculated 65 weeks ago
0.3756info.png
0.3756

Last calculated 65 weeks ago
64 Electrical Transmission Line '64' Lower Tumut Electrical Bus 34 Upper Tumut Electrical Bus 36 0.001607info.png
0.001607

Last calculated 65 weeks ago
0.01295info.png
0.01295

Last calculated 65 weeks ago
0.1508info.png
0.1508

Last calculated 65 weeks ago
65 Electrical Transmission Line '65' Murray Electrical Bus 35 Upper Tumut Electrical Bus 36 0.001797info.png
0.001797

Last calculated 65 weeks ago
0.01464info.png
0.01464

Last calculated 65 weeks ago
0.1679info.png
0.1679

Last calculated 65 weeks ago
66 Electrical Transmission Line '66' Lower Tumut Electrical Bus 34 Murray Electrical Bus 35 0.00298info.png
0.00298

Last calculated 65 weeks ago
0.0263info.png
0.0263

Last calculated 65 weeks ago
0.2778info.png
0.2778

Last calculated 65 weeks ago
67 Electrical Transmission Line '67' Murray Electrical Bus 35 Dederang Electrical Bus 40 0.00441info.png
0.00441

Last calculated 65 weeks ago
0.0343info.png
0.0343

Last calculated 65 weeks ago
0.451info.png
0.451

Last calculated 65 weeks ago
68 Electrical Transmission Line '68' Murray Electrical Bus 35 Dederang Electrical Bus 40 0.00438info.png
0.00438

Last calculated 65 weeks ago
0.0342info.png
0.0342

Last calculated 65 weeks ago
0.45info.png
0.45

Last calculated 65 weeks ago
6 Electrical Transmission Line '6' Canberra Electrical Bus 33 Kangaroo Valley Electrical Bus 32 0.007606info.png
0.007606

Last calculated 65 weeks ago
0.05851info.png
0.05851

Last calculated 65 weeks ago
0.7172info.png
0.7172

Last calculated 65 weeks ago
71 Electrical Transmission Line '71' Mount Piper Powerstation 20 Wallerawang Powerstation 22 0.000232info.png
0.000232

Last calculated 65 weeks ago
0.002224info.png
0.002224

Last calculated 65 weeks ago
0.03231info.png
0.03231

Last calculated 65 weeks ago
72 Electrical Transmission Line '72' Mount Piper Powerstation 20 Wellington Electrical Bus 21 0.006822info.png
0.006822

Last calculated 65 weeks ago
0.05214info.png
0.05214

Last calculated 65 weeks ago
0.6582info.png
0.6582

Last calculated 65 weeks ago
73 Electrical Transmission Line '73' Bayswater Powerstation 9 Mount Piper Powerstation 20 0.003977info.png
0.003977

Last calculated 65 weeks ago
0.05241info.png
0.05241

Last calculated 65 weeks ago
1.131info.png
1.131

Last calculated 65 weeks ago
74 Electrical Transmission Line '74' Bayswater Powerstation 9 Wallerawang Powerstation 22 0.004215info.png
0.004215

Last calculated 65 weeks ago
0.05465info.png
0.05465

Last calculated 65 weeks ago
1.164info.png
1.164

Last calculated 65 weeks ago
76 Electrical Transmission Line '76' Sydney South Electrical Bus 26 Wallerawang Powerstation 22 0.004063info.png
0.004063

Last calculated 65 weeks ago
0.03998info.png
0.03998

Last calculated 65 weeks ago
0.5879info.png
0.5879

Last calculated 65 weeks ago
77 Electrical Transmission Line '77' Wallerawang Powerstation 22 Ingleburn Electrical Bus 23 0.003455info.png
0.003455

Last calculated 65 weeks ago
0.03408info.png
0.03408

Last calculated 65 weeks ago
0.4989info.png
0.4989

Last calculated 65 weeks ago
78 Electrical Transmission Line '78' Sydney South Electrical Bus 26 Ingleburn Electrical Bus 23 6e-04info.png
0.0006

Last calculated 65 weeks ago
0.005834info.png
0.005834

Last calculated 65 weeks ago
0.08791info.png
0.08791

Last calculated 65 weeks ago
81 Electrical Transmission Line '81' Liddell Powerstation 6 Newcastle Electrical Bus 7 0.002914info.png
0.002914

Last calculated 65 weeks ago
0.03058info.png
0.03058

Last calculated 65 weeks ago
0.3856info.png
0.3856

Last calculated 65 weeks ago
82 Electrical Transmission Line '82' Liddell Powerstation 6 Tomago Aluminium Smelter 8 0.003321info.png
0.003321

Last calculated 65 weeks ago
0.03439info.png
0.03439

Last calculated 65 weeks ago
0.4456info.png
0.4456

Last calculated 65 weeks ago
83 Electrical Transmission Line '83' Muswellbrook Electrical Bus 5 Liddell Powerstation 6 0.000704info.png
0.000704

Last calculated 65 weeks ago
0.005355info.png
0.005355

Last calculated 65 weeks ago
0.06807info.png
0.06807

Last calculated 65 weeks ago
84 Electrical Transmission Line '84' Tamworth Electrical Bus 4 Liddell Powerstation 6 0.005571info.png
0.005571

Last calculated 65 weeks ago
0.04298info.png
0.04298

Last calculated 65 weeks ago
0.5282info.png
0.5282

Last calculated 65 weeks ago
85 Electrical Transmission Line '85' Armidale Electrical Bus 2 Tamworth Electrical Bus 4 0.005511info.png
0.005511

Last calculated 65 weeks ago
0.03368info.png
0.03368

Last calculated 65 weeks ago
0.4255info.png
0.4255

Last calculated 65 weeks ago
86 Electrical Transmission Line '86' Armidale Electrical Bus 2 Tamworth Electrical Bus 4 0.004159info.png
0.004159

Last calculated 65 weeks ago
0.03198info.png
0.03198

Last calculated 65 weeks ago
0.3921info.png
0.3921

Last calculated 65 weeks ago
88 Electrical Transmission Line '88' Tamworth Electrical Bus 4 Muswellbrook Electrical Bus 5 0.005086info.png
0.005086

Last calculated 65 weeks ago
0.03913info.png
0.03913

Last calculated 65 weeks ago
0.4798info.png
0.4798

Last calculated 65 weeks ago
89 Electrical Transmission Line '89' Armidale Electrical Bus 2 Lismore Electrical Bus 3 0.01222info.png
0.01222

Last calculated 65 weeks ago
0.09428info.png
0.09428

Last calculated 65 weeks ago
1.149info.png
1.149

Last calculated 65 weeks ago
8C Electrical Transmission Line '8C' Dumaresq Substation 1 Armidale Electrical Bus 2 0.003762info.png
0.003762

Last calculated 65 weeks ago
0.04763info.png
0.04763

Last calculated 65 weeks ago
0.7243info.png
0.7243

Last calculated 65 weeks ago
8E Electrical Transmission Line '8E' Dumaresq Substation 1 Armidale Electrical Bus 2 0.003798info.png
0.003798

Last calculated 65 weeks ago
0.04755info.png
0.04755

Last calculated 65 weeks ago
0.7243info.png
0.7243

Last calculated 65 weeks ago
8 Electrical Transmission Line '8' Marulan Electrical Bus 30 Dapto Electrical Bus 28 0.002849info.png
0.002849

Last calculated 65 weeks ago
0.02182info.png
0.02182

Last calculated 65 weeks ago
0.2681info.png
0.2681

Last calculated 65 weeks ago
92 Electrical Transmission Line '92' Newcastle Electrical Bus 7 Vales Point Powerstation 12 0.001345info.png
0.001345

Last calculated 65 weeks ago
0.01389info.png
0.01389

Last calculated 65 weeks ago
0.1797info.png
0.1797

Last calculated 65 weeks ago
93 Electrical Transmission Line '93' Newcastle Electrical Bus 7 Eraring Powerstation 10 0.000608info.png
0.000608

Last calculated 65 weeks ago
0.005983info.png
0.005983

Last calculated 65 weeks ago
0.08332info.png
0.08332

Last calculated 65 weeks ago
94 Electrical Transmission Line '94' Newcastle Electrical Bus 7 Tomago Aluminium Smelter 8 0.000717info.png
0.000717

Last calculated 65 weeks ago
0.006888info.png
0.006888

Last calculated 65 weeks ago
0.0993info.png
0.0993

Last calculated 65 weeks ago
99B Electrical Transmission Line '99B' Jindera Electrical Bus 38 Wodonga Electrical Bus 39 0.00164info.png
0.00164

Last calculated 65 weeks ago
0.01259info.png
0.01259

Last calculated 65 weeks ago
0.1646info.png
0.1646

Last calculated 65 weeks ago
99X Electrical Transmission Line '99X' Wodonga Electrical Bus 39 Dederang Electrical Bus 40 0.0016info.png
0.0016

Last calculated 65 weeks ago
0.0123info.png
0.0123

Last calculated 65 weeks ago
0.16226info.png
0.16226

Last calculated 65 weeks ago
… further results


A map of all the Electrical Busses.


Photos
Videos
Wikipedia
Webcams


A number of assumptions have been made about the network:

Property Value
Shunt Conductance (MW demanded at V = 1.0 p.u.) 0 MWinfo.png
'0 MW'

Last calculated 65 weeks ago
Bus Area 1info.png
1

Last calculated 73 weeks ago
Voltage Angle 0info.png
0

Last calculated 73 weeks ago
Base Voltage 330 kVinfo.png
'330 kV'

Last calculated 65 weeks ago
Loss Zone 1info.png
1

Last calculated 73 weeks ago
Maximum Voltage Magnitude (p.u.) 1.1info.png
1.1

Last calculated 65 weeks ago
Minimum Voltage Magnitude (p.u.) 0.9info.png
0.9

Last calculated 72 weeks ago
Reactive Power Output 0 MVAinfo.png
'0 MVA'

Last calculated 65 weeks ago
Base MVA 100 MVAinfo.png
'100 MVA'

Last calculated 65 weeks ago
Maximum Real Power Output 9999 MWinfo.png
'9999 MW'

Last calculated 65 weeks ago
Minimum Real Power Output 9999 MWinfo.png
'9999 MW'

Last calculated 65 weeks ago

Network Description

The following Octave[3] code selects the buses, generators, branches and network defaults and forms up matricies that can be processed by the Matpower [4] Load Flow analysis tool. It then displays the outputs of the computation on a map.


%Graphical Controls
%P_N = 4; %Number of plots (ONE OR FOUR)
%SCENARIO_NUMBERS = [ 1, 3,     4,     5];                 %DEFINE SCENARIOS TO COMPARE(between 0 and 5) IN MULTIPLOT. For single plot, first scenario is 
%V_P = [-33.75,150.5,8]; %Define map viewpoint (lat, long, zoom number). Zoom number MUST be between 6 & displayed.

L_P = [25,100]; % Legend Position ([0,0] = TOP LEFT)

for PLOTS = 1:P_N;

mpc.version = '2';

%Get the NSW network defaults baseMVA
%MachineBaseMva ?MaximumRealPowerOutput ?MinimumRealPowerOutput ?MaximumVoltageMagnitude ?MinimumVoltageMagnitude ?BusArea ?VoltageAngle ?BaseVoltage ?ShuntConductance
nswDefaults = sparql("SELECT ?MachineBaseMva ?MaximumRealPowerOutput ?MinimumRealPowerOutput ?MaximumVoltageMagnitude ?MinimumVoltageMagnitude ?BusArea ?VoltageAngle ?BaseVoltage ?ShuntConductance WHERE{
a:NSW_Electricity_Network prop:MachineBaseMva ?MachineBaseMva.
a:NSW_Electricity_Network prop:MaximumRealPowerOutput ?MaximumRealPowerOutput.
a:NSW_Electricity_Network prop:MinimumRealPowerOutput ?MinimumRealPowerOutput.
a:NSW_Electricity_Network prop:MaximumVoltageMagnitude ?MaximumVoltageMagnitude.
a:NSW_Electricity_Network prop:MinimumVoltageMagnitude ?MinimumVoltageMagnitude.
a:NSW_Electricity_Network prop:BusArea ?BusArea.
a:NSW_Electricity_Network prop:VoltageAngle ?VoltageAngle.
a:NSW_Electricity_Network prop:BaseVoltage ?BaseVoltage.
a:NSW_Electricity_Network prop:ShuntConductance ?ShuntConductance.
}");

%Parse the baseMVA
mpc.baseMVA = sscanf (char(nswDefaults(1,1)), "%f MVA");


%?bus ?id ?BusType ?realPowerDemand ?reactivPowerDemand ?shuntSusceptance ?VoltageMagnitude ?RealPowerOutput ?MaximumReactivePowerOutput ?MinimumReactivePowerOutput ?VoltageMagnitudeSetpoint ?GeoLocation ?ShortName
bus = sparql("select ?bus ?id ?BusType ?realPowerDemand ?reactivPowerDemand ?shuntSusceptance ?VoltageMagnitude ?RealPowerOutput 
  ?MaximumReactivePowerOutput ?MinimumReactivePowerOutput ?VoltageMagnitudeSetpoint ?GeoLocation ?ShortName WHERE{
?bus rdf:type cat:Electrical_Bus.
OPTIONAL{?bus prop:BusIdentifier ?id.}
OPTIONAL{?bus prop:BusType ?BusType.}
OPTIONAL{?bus prop:RealPowerDemand ?realPowerDemand.}
OPTIONAL{?bus prop:ReactivPowerDemand ?reactivPowerDemand.}
OPTIONAL{?bus prop:ShuntSusceptance ?shuntSusceptance.}
OPTIONAL{?bus prop:VoltageMagnitude ?VoltageMagnitude.}
OPTIONAL{?bus prop:RealPowerOutput ?RealPowerOutput.}
OPTIONAL{?bus prop:MaximumReactivePowerOutput ?MaximumReactivePowerOutput.}
OPTIONAL{?bus prop:MinimumReactivePowerOutput ?MinimumReactivePowerOutput.}
OPTIONAL{?bus prop:VoltageMagnitudeSetpoint ?VoltageMagnitudeSetpoint.}
OPTIONAL{?bus prop:GeoLocation ?GeoLocation.}
OPTIONAL{?bus prop:ShortName ?ShortName.}
}
ORDER BY ?bus");

[m,n] =size(bus);
B_n = n;
busCtr = 0;
genCtr = 0;

%----------------------------------------------------------------------------------------------------------------%
%Build up the mpc.bus tablefrom bus. Take defaults from the nswDefaults variable [Bus Data (mpc.bus) matrix definition ]

for i = 1:n
BUS_I(i) = cell2mat(bus(2,i)); %BUS_I  (column 1) Bus number
BUS_TYPE(i) = cell2mat(bus(3,i)); %BUS_TYPE  (colomn 2) Bus type (1 = PQ, 2 = PV, 3 = ref, 4 = isolated)
PD(i) =sscanf (char(bus(4,i)), "%d MW"); %PD  (column 3) Real Power Demand
QD(i) =sscanf (char(bus(5,i)), "%d MVA"); %QD  (column 4) Reactive Power Demand
GS(i) =  sscanf (char(nswDefaults(9,1)), "%f MW"); %GS (Column 5)  Shunt Conductance
if (length (char(bus(6,i)))<4), BS(i) = 0;  %BS (Column 6)  Shunt Susceptance
else BS(i) = sscanf (char(bus(6,i)), "%f MVA"); 
endif;
BUS_AREA(i) = cell2mat(nswDefaults(6,1)); %BUS_AREA (Column 7)  Area Number
if (length (char(bus(7,i)))<1), VM(i) = 1; 
else VM(i) = cell2mat(bus(7,i));  %VM (Column 8) Volatge Magnitude
endif;
VA(i) = cell2mat(nswDefaults(7,1)); %VA (Column 9) Volatge Angle
BASE_KV(i) = sscanf (char(nswDefaults(8,1)), "%f kV"); %BASE_KV (Column 10) Base Volatge (kV)
ZONE(i) = 1; %ZONE  (Column 11) Loss Zone
VMAX(i) = cell2mat(nswDefaults(4,1)); %VMAX (Coulmn 12)  Maximum voltage magnitude
VMIN(i) = cell2mat(nswDefaults(5,1)); %VMIN (Coulmn 13)  Minimum voltage magnitude
endfor

mpc.bus = [BUS_I', BUS_TYPE', PD',QD',GS',BS',BUS_AREA',VM',VA', BASE_KV',ZONE',VMAX',VMIN'];

%----------------------------------------------------------------------------------------------------------------%
%Build up mpc.gen table from bus. Take defaults from the nswDefaults variable

for j = 1:n

if cell2mat(bus(3,j)) ==1, GEN_BUS(j)=0; %GEN_BUS (Column 1) Bus number
else GEN_BUS(j)=cell2mat(bus(2,j)); 
endif
if (length (char(bus(8,j)))<4), PG(j) = 0;  %PG (Column 2) Real Power Output
else PG(j) = sscanf (char(bus(8,j)), "%f MW"); 
endif;
QG(j) = 0; %QG (Column 3)  Reactive Power Output  !HARD CODED
if (length (char(bus(9,j)))<4), QMAX(j) = 0;  %QMAX (Column 4) Maximum Reactive Power Output
else QMAX(j) = sscanf (char(bus(9,j)), "%f MVA"); 
endif;
if (length (char(bus(10,j)))<4), QMIN(j) = 0;  %QMIN (Column 5) Minimum Reactive Power Output
else QMIN(j) = sscanf (char(bus(10,j)), "%f MVA"); 
endif;
if (length (char(bus(11,j)))<1), VG(j) = 0;  %VG (Column 6) Voltage Magnitude
else VG(j) = cell2mat(bus(11,j)); 
endif;
MBASE(j) =  sscanf (char(nswDefaults(1,1)), "%f MVA"); %MBASE (Column 7) Total MVA base of machine !HARDCODED
GEN_STATUS(j) = 1; %GEN_STATUS (Column 8) Machine Status !HARDCODED
PMAX(j) = 9999; %PMAX (Column 9) Maximum reaal power output !HARDCODED
PMIN(j) = 9999; %PMIN (Column 10) Minumum real power output !HARDCODED
endfor

mpc.gen = [GEN_BUS', PG',QG', QMAX', QMIN', VG',MBASE',GEN_STATUS',PMAX',PMIN'];

a = mpc.gen;
[s, i] = sort (a (:, 1)); 
mpc.gen = a (i, :);
len = length (mpc.gen)-length (mpc.gen(:,1)(mpc.gen(:,1)~=0)); %Delete zero leading rows
mpc.gen(1:len,:)=[];

%----------------------------------------------------------------------------------------------------------------%
%Get the line data
%?line ?ShortName ?FromBus ?FromId ?FromGeoLocation ?ToBus ?ToId ?ToGeoLocation ?Resistance ?Reactance ?TotalLineChargingSusceptance
branch = sparql("select ?line ?ShortName ?FromBus ?FromId ?FromGeoLocation ?ToBus ?ToId ?ToGeoLocation ?Resistance ?Reactance ?TotalLineChargingSusceptance WHERE{
?line rdf:type cat:Electrical_Transmission.
?line prop:ShortName ?ShortName.
?line prop:FromBus ?FromBus.
?FromBus prop:BusIdentifier ?FromId.
?FromBus prop:GeoLocation ?FromGeoLocation.
?line prop:ToBus ?ToBus.
?ToBus prop:BusIdentifier ?ToId.
?ToBus prop:GeoLocation ?ToGeoLocation.
?line prop:Resistance ?Resistance.
?line prop:Reactance ?Reactance.
?line prop:TotalLineChargingSusceptance  ?TotalLineChargingSusceptance .
} ORDER BY ?line");

%build up the mpc.branch matrix
[m,n] =size(branch);
for i = 1:n

FBUS(i) = cell2mat(branch(4,i)); %F_Bus From Bus  (column 1)
TBUS(i) = cell2mat(branch(7,i)); %T_BUS To Bus  (column 2)
BR_R(i) = cell2mat(branch(9,i)); %BR_R Resistance  (column 3)
BR_X(i) = cell2mat(branch(10,i)); %BR_X Reactance  (column 4)
BR_B(i) = cell2mat(branch(11,i)); %BR_B Total line charging suseptance  (column 5)
RATE_A(i) = 1000; %RATE_A MVA rating A (column 6) !HARDCODED
RATE_B(i) = 1000; %RATE_B MVA rating B (column 7) !HARDCODED
RATE_C(i) = 1000; %RATE_C MVA rating C (column 8) !HARDCODED
TAP(i) = 0; %TAP Transformer off Nominal turns ratio (column 9) !HARDCODED
SHIFT(i) = 0; %SHIFT Transformer phase shift angles (column 10) !HARDCODED
BR_STATUS(i) = 1; %BR_STATUS Initial branch status (column 11) !HARDCODED
ANGMIN(i) = -360; %ANGMIN Minimum angle difference (column 12) !HARDCODED
ANGMAX(i) = 360; %ANDMAX Maximum angle difference (column 13) !HARDCODED

endfor

mpc.branch = [FBUS',TBUS',BR_R',BR_X',BR_B', RATE_A', RATE_B', RATE_C', TAP',SHIFT',BR_STATUS',ANGMIN',ANGMAX'];
mpc;

%----------------------------------------------------------------------------------------------------------------%
%Get the scenario data

%----------------------------------------------------------------------------------------------------------------%
%Evaluate the relevant case

SCENARIO = [ 1, 0.6, 1.05, 1.05, 1.05, 1.05; %SCENARIO MATRIX
1, 0.6, 1.05, 1.05, 1.05, 1.05;
1, 1, 1, 0, 1, 1;
1, 1, 1, 1, 0, 1;
1, 1, 1, 1, 1, 0 ];

SN = SCENARIO_NUMBERS(1,PLOTS)+1;

% Scenario 0 - No modification
% Scenario 1 - Base case (light load)
% Scenario 2 - Heavy load 
% Scenario 3 - Heavy loading plus outage of the lines 12-13 and 12-16 (Vales - Sydney North Vales - Munmorah)
% Scenario 4 - Heavy loading plus outage of generator #6 (Liddell)
% Scenario 5 - Heavy loading plus outage of the lines 9-19 and 9-20 (Bayswater - Regentville and Bayswater - Mt Piper)

mpc.bus(:,[3:4]) = SCENARIO(1,SN)*mpc.bus(:,[3:4]);
mpc.gen(:,[2]) = SCENARIO(2,SN)*mpc.gen(:,[2]);
mpc.branch([18,19],[11]) = SCENARIO(3,SN)*mpc.branch([18,19],[11]);
mpc.gen([2],[8]) = SCENARIO(4,SN)*mpc.gen([2],[8]);
mpc.branch([27,48],[11]) = SCENARIO(5,SN)*mpc.branch([27,48],[11]);
diary off;
results = runpf(mpc, mpoption( "OUT_ALL",0));
diary on;
%----------------------------------------------------------------------------------------------------------------%
%Apparent Power Calculation

PP_IN = diag(results.branch(:,14)*(results.branch(:,14))'); %Test Apparent power in is under limit
QQ_IN = diag(results.branch(:,15)*(results.branch(:,15))');
APPARENT_POWER_IN = sqrt (PP_IN+QQ_IN);

PP_OUT = diag(results.branch(:,16)*(results.branch(:,16))'); %Test Apparent power out is under limit
QQ_OUT = diag(results.branch(:,17)*(results.branch(:,17))');
APPARENT_POWER_OUT = sqrt (PP_OUT+QQ_OUT);

for t = 1: n
if APPARENT_POWER_IN(t,1) < 1000, APPARENT_POWER_TEST(t,1) =1; else APPARENT_POWER_TEST(t,1) = 0; endif; %Ensure BOTH power in AND power out is under limit. 
if APPARENT_POWER_OUT(t,1) < 1000, APPARENT_POWER_TEST(t,2) =1; else APPARENT_POWER_TEST(t,2) = 0; endif;
APPARENT_POWER_TEST(t,3) = APPARENT_POWER_TEST(t,2) + APPARENT_POWER_TEST(t,1);
endfor

%----------------------------------------------------------------------------------------------------------------%
%PLOTTING THE RESULTS
%----------------------------------------------------------------------------------------------------------------%


%Set Viewpoint characteristics (i.e. Lat, Long, zoom) ZOOM LEVEL MUST BE BETWEEN 6 & 8 



%V_P = [-33.75, 150.5, 8]; %Reasonable co-ordinates for zoom level 8 (Sydney level)
%V_P = [-33, 149, 7]; %Reasonable co-ordinates for zoom level 7 (North Coast level)
%V_P = [-33, 148, 6]; %Reasonable co-ordinates for zoom level 6 (NSW level)


%----------------------------------------------------------------------------------------------------------------%
%Get Image from Google Maps API

LAT = num2str(V_P(1,1));
LONG = num2str(V_P(1,2));
ZOOM = num2str(V_P(1,3));

CENTER = strcat (LAT,",",LONG);

mapbgf = urlread(strcat("http://maps.google.com/maps/api/staticmap?center=",CENTER,"&zoom=",ZOOM,"&size=640x640&sensor=false&style=feature:all",char(124),"visibility:off&style=feature:water",char(124),"saturation:-100",char(124),"lightness:-50",char(124),"visibility:on&style=feature:administrative.province",char(124),"element:geometry",char(124),"visibility:on&style=feature:landscape",char(124),"element:geometry",char(124),"lightness:100",char(124),"visibility:on"));
bgfid=fopen("static.png","w");
fwrite(bgfid,mapbgf,"char");
fclose(bgfid);

I = imread("static.png");

if P_N == 4, PLOT_GRID = 2; else PLOT_GRID = 1; endif
subplot(PLOT_GRID,PLOT_GRID,PLOTS); imshow(I)
hold on

%----------------------------------------------------------------------------------------------------------------%

%Zoom Level transform variables

T  = [ 6.05, 2.97, 1.452; %Transform matrix 
6.9, 3.5, 1.75;
45.87, 90.9, 183.15;
53.145, 108.67, 221.23 ];

Z_n = V_P(1,3) - 5; %Zoom Number (1,2 or 3, where 1 = "zoom level 6", 2 = "zoom level 7", 3 = "zoom level 8")


Y = -V_P (1,1) - T(1,Z_n); %translation in y
X = V_P(1,2) - T(2,Z_n); %translation in x
XD = T(3,Z_n); %dilation factor x
YD = T(4,Z_n); %dilation factor y


%----------------------------------------------------------------------------------------------------------------%

%Color scheme


%Ok lines- 4284d3 – 66,132,211, width 1 
LINES = [11,97,164]/256;
%Err lines - FFB540 – 255,181,64, width 2
ERRLINES = [255,191,0]/256;
%Out lines – 000000 – 0,0,0 – width 1, could we make these dashed?
OUTLINES = [255,49,0]/256;
%Ok bus - 04346C – 4, 52,108 – Radius 1 - Dot
BUS_COL = [3,62,107]/256;
%Err bus - A66500 – 166,101,0 – Radius 2 – Concentric circles
ERR_BUS_COL = [166,124,0]/256;

%----------------------------------------------------------------------------------------------------------------%
subplot(PLOT_GRID,PLOT_GRID,PLOTS);

%Plot Branches
for i = 1:n
c = [sscanf(char(branch(5,i)), '%f,'),sscanf(char(branch(8,i)), '%f,')];
if ((APPARENT_POWER_IN(i,1)) == 0),
plot (((c(2,:)-X)*XD),((Y+c(1,:))*-YD),"color",OUTLINES,"linewidth",2,"linestyle","--")
elseif (abs (APPARENT_POWER_TEST(i,3)) < 2), 
plot (((c(2,:)-X)*XD),((Y+c(1,:))*-YD),"color",ERRLINES,"linewidth",2)
else plot (((c(2,:)-X)*XD),((Y+c(1,:))*-YD),"color",LINES,"linewidth",1)
endif
endfor


%Plot Buses
for i = 1:B_n
GL = sscanf(char(bus(12,i)), '%f,');
if results.bus(i,8) > 1.1
for k = 1:3 plot (((GL(2,1)-X)*XD),((Y+GL(1,1))*-YD), "o","color",ERR_BUS_COL, "markersize",(2*k)) endfor
elseif results.bus(i,8) < 0.9 
for k = 1:3 plot (((GL(2,1)-X)*XD),((Y+GL(1,1))*-YD), "o","color",ERR_BUS_COL,"markersize",(2*k)) endfor
else for k = 1:3 plot (((GL(2,1)-X)*XD),((Y+GL(1,1))*-YD), "o","color",BUS_COL,"markersize",k) endfor
endif
endfor

%----------------------------------------------------------------------------------------------------------------%

%Plot Legend



L_L = 3; % Line length in legend
if P_N == 4, L_S = 20; else L_S = 10; endif; % Legend Spacing  
if P_N == 4, L_TS = 30; else L_TS = 15; endif; % Legend Title spacing

%Legend images
for k = 1:3 plot (L_P(1,1), L_P(1,2)+L_TS, "o","color",ERR_BUS_COL,"markersize",(2*k)) endfor
for k = 1:3 plot (L_P(1,1), L_P(1,2)+L_TS+L_S, "o","color",BUS_COL,"markersize",k) endfor
plot ([(L_P(1,1)-L_L),(L_P(1,1)+L_L)],[L_P(1,2)+L_TS+2*L_S,L_P(1,2)+L_TS+2*L_S],"color",OUTLINES,"linewidth",2)
plot ([(L_P(1,1)-L_L),(L_P(1,1)+L_L)],[L_P(1,2)+L_TS+3*L_S,L_P(1,2)+L_TS+3*L_S],"color",ERRLINES,"linewidth",2)
plot ([(L_P(1,1)-L_L),(L_P(1,1)+L_L)],[L_P(1,2)+L_TS+4*L_S,L_P(1,2)+L_TS+4*L_S],"color",LINES,"linewidth",1)

%Legend text
text (L_P(1,1), L_P(1,2), "LEGEND")
text (L_P(1,1)+10, L_P(1,2)+L_TS, "Bus Voltage Outside Limits")
text (L_P(1,1)+10, L_P(1,2)+L_TS+L_S, "Bus OK")
text (L_P(1,1)+10, L_P(1,2)+L_TS+2*L_S, "Line Outage")
text (L_P(1,1)+10, L_P(1,2)+L_TS+3*L_S, "Line Overload")
text (L_P(1,1)+10, L_P(1,2)+L_TS+4*L_S, "Line Ok")

%TITLE

SCEN_TITLE = strcat ("Scenario",{" "},num2str(SN-1));
text (300,0, SCEN_TITLE)

%----------------------------------------------------------------------------------------------------------------%
endfor

hold off

References

Facts about NSW Electricity NetworkRDF feed
BaseVoltage330 kVinfo.png
'330 kV'
  +
BusArea1info.png
1
  +
CapitalCost1.30 AUD 2009/wattinfo.png
'1.30 AUD 2009/watt'
  +
CentralizedEfficiency0.945info.png
 0.945
  +
CiteMatthew Sullivan  +, and Calculating Demand for NGACs  +
DistributedEfficiency0.96info.png
 0.96
  +
Emissions0.906t/MWhinfo.png
'0.906t/MWh'
  +
LossZone1info.png
1
  +
MachineBaseMva100 MVAinfo.png
'100 MVA'
  +
MaximumRealPowerOutput9999 MWinfo.png
'9999 MW'
  +
MaximumVoltageMagnitude1.1info.png
1.1
  +
MinimumRealPowerOutput9999 MWinfo.png
'9999 MW'
  +
MinimumVoltageMagnitude0.9info.png
0.9
  +
Octave

%Graphical Controls
%P_N = 4; %Number of plots (ONE OR FOUR)
%SCENARIO_NUMBERS = [ 1, 3,     4,     5];                 %DEFINE SCENARIOS TO COMPARE(between 0 and 5) IN MULTIPLOT. For single plot, first scenario is 
%V_P = [-33.75,150.5,8]; %Define map viewpoint (lat, long, zoom number). Zoom number MUST be between 6 & displayed.

L_P = [25,100]; % Legend Position ([0,0] = TOP LEFT)

for PLOTS = 1:P_N;

mpc.version = '2';

%Get the NSW network defaults baseMVA
%MachineBaseMva ?MaximumRealPowerOutput ?MinimumRealPowerOutput ?MaximumVoltageMagnitude ?MinimumVoltageMagnitude ?BusArea ?VoltageAngle ?BaseVoltage ?ShuntConductance
nswDefaults = sparql("SELECT ?MachineBaseMva ?MaximumRealPowerOutput ?MinimumRealPowerOutput ?MaximumVoltageMagnitude ?MinimumVoltageMagnitude ?BusArea ?VoltageAngle ?BaseVoltage ?ShuntConductance WHERE{
a:NSW_Electricity_Network prop:MachineBaseMva ?MachineBaseMva.
a:NSW_Electricity_Network prop:MaximumRealPowerOutput ?MaximumRealPowerOutput.
a:NSW_Electricity_Network prop:MinimumRealPowerOutput ?MinimumRealPowerOutput.
a:NSW_Electricity_Network prop:MaximumVoltageMagnitude ?MaximumVoltageMagnitude.
a:NSW_Electricity_Network prop:MinimumVoltageMagnitude ?MinimumVoltageMagnitude.
a:NSW_Electricity_Network prop:BusArea ?BusArea.
a:NSW_Electricity_Network prop:VoltageAngle ?VoltageAngle.
a:NSW_Electricity_Network prop:BaseVoltage ?BaseVoltage.
a:NSW_Electricity_Network prop:ShuntConductance ?ShuntConductance.
}");

%Parse the baseMVA
mpc.baseMVA = sscanf (char(nswDefaults(1,1)), "%f MVA");


%?bus ?id ?BusType ?realPowerDemand ?reactivPowerDemand ?shuntSusceptance ?VoltageMagnitude ?RealPowerOutput ?MaximumReactivePowerOutput ?MinimumReactivePowerOutput ?VoltageMagnitudeSetpoint ?GeoLocation ?ShortName
bus = sparql("select ?bus ?id ?BusType ?realPowerDemand ?reactivPowerDemand ?shuntSusceptance ?VoltageMagnitude ?RealPowerOutput 
  ?MaximumReactivePowerOutput ?MinimumReactivePowerOutput ?VoltageMagnitudeSetpoint ?GeoLocation ?ShortName WHERE{
?bus rdf:type cat:Electrical_Bus.
OPTIONAL{?bus prop:BusIdentifier ?id.}
OPTIONAL{?bus prop:BusType ?BusType.}
OPTIONAL{?bus prop:RealPowerDemand ?realPowerDemand.}
OPTIONAL{?bus prop:ReactivPowerDemand ?reactivPowerDemand.}
OPTIONAL{?bus prop:ShuntSusceptance ?shuntSusceptance.}
OPTIONAL{?bus prop:VoltageMagnitude ?VoltageMagnitude.}
OPTIONAL{?bus prop:RealPowerOutput ?RealPowerOutput.}
OPTIONAL{?bus prop:MaximumReactivePowerOutput ?MaximumReactivePowerOutput.}
OPTIONAL{?bus prop:MinimumReactivePowerOutput ?MinimumReactivePowerOutput.}
OPTIONAL{?bus prop:VoltageMagnitudeSetpoint ?VoltageMagnitudeSetpoint.}
OPTIONAL{?bus prop:GeoLocation ?GeoLocation.}
OPTIONAL{?bus prop:ShortName ?ShortName.}
}
ORDER BY ?bus");

[m,n] =size(bus);
B_n = n;
busCtr = 0;
genCtr = 0;

%----------------------------------------------------------------------------------------------------------------%
%Build up the mpc.bus tablefrom bus. Take defaults from the nswDefaults variable [Bus Data (mpc.bus) matrix definition ]

for i = 1:n
BUS_I(i) = cell2mat(bus(2,i)); %BUS_I  (column 1) Bus number
BUS_TYPE(i) = cell2mat(bus(3,i)); %BUS_TYPE  (colomn 2) Bus type (1 = PQ, 2 = PV, 3 = ref, 4 = isolated)
PD(i) =sscanf (char(bus(4,i)), "%d MW"); %PD  (column 3) Real Power Demand
QD(i) =sscanf (char(bus(5,i)), "%d MVA"); %QD  (column 4) Reactive Power Demand
GS(i) =  sscanf (char(nswDefaults(9,1)), "%f MW"); %GS (Column 5)  Shunt Conductance
if (length (char(bus(6,i)))<4), BS(i) = 0;  %BS (Column 6)  Shunt Susceptance
else BS(i) = sscanf (char(bus(6,i)), "%f MVA"); 
endif;
BUS_AREA(i) = cell2mat(nswDefaults(6,1)); %BUS_AREA (Column 7)  Area Number
if (length (char(bus(7,i)))<1), VM(i) = 1; 
else VM(i) = cell2mat(bus(7,i));  %VM (Column 8) Volatge Magnitude
endif;
VA(i) = cell2mat(nswDefaults(7,1)); %VA (Column 9) Volatge Angle
BASE_KV(i) = sscanf (char(nswDefaults(8,1)), "%f kV"); %BASE_KV (Column 10) Base Volatge (kV)
ZONE(i) = 1; %ZONE  (Column 11) Loss Zone
VMAX(i) = cell2mat(nswDefaults(4,1)); %VMAX (Coulmn 12)  Maximum voltage magnitude
VMIN(i) = cell2mat(nswDefaults(5,1)); %VMIN (Coulmn 13)  Minimum voltage magnitude
endfor

mpc.bus = [BUS_I', BUS_TYPE', PD',QD',GS',BS',BUS_AREA',VM',VA', BASE_KV',ZONE',VMAX',VMIN'];

%----------------------------------------------------------------------------------------------------------------%
%Build up mpc.gen table from bus. Take defaults from the nswDefaults variable

for j = 1:n

if cell2mat(bus(3,j)) ==1, GEN_BUS(j)=0; %GEN_BUS (Column 1) Bus number
else GEN_BUS(j)=cell2mat(bus(2,j)); 
endif
if (length (char(bus(8,j)))<4), PG(j) = 0;  %PG (Column 2) Real Power Output
else PG(j) = sscanf (char(bus(8,j)), "%f MW"); 
endif;
QG(j) = 0; %QG (Column 3)  Reactive Power Output  !HARD CODED
if (length (char(bus(9,j)))<4), QMAX(j) = 0;  %QMAX (Column 4) Maximum Reactive Power Output
else QMAX(j) = sscanf (char(bus(9,j)), "%f MVA"); 
endif;
if (length (char(bus(10,j)))<4), QMIN(j) = 0;  %QMIN (Column 5) Minimum Reactive Power Output
else QMIN(j) = sscanf (char(bus(10,j)), "%f MVA"); 
endif;
if (length (char(bus(11,j)))<1), VG(j) = 0;  %VG (Column 6) Voltage Magnitude
else VG(j) = cell2mat(bus(11,j)); 
endif;
MBASE(j) =  sscanf (char(nswDefaults(1,1)), "%f MVA"); %MBASE (Column 7) Total MVA base of machine !HARDCODED
GEN_STATUS(j) = 1; %GEN_STATUS (Column 8) Machine Status !HARDCODED
PMAX(j) = 9999; %PMAX (Column 9) Maximum reaal power output !HARDCODED
PMIN(j) = 9999; %PMIN (Column 10) Minumum real power output !HARDCODED
endfor

mpc.gen = [GEN_BUS', PG',QG', QMAX', QMIN', VG',MBASE',GEN_STATUS',PMAX',PMIN'];

a = mpc.gen;
[s, i] = sort (a (:, 1)); 
mpc.gen = a (i, :);
len = length (mpc.gen)-length (mpc.gen(:,1)(mpc.gen(:,1)~=0)); %Delete zero leading rows
mpc.gen(1:len,:)=[];

%----------------------------------------------------------------------------------------------------------------%
%Get the line data
%?line ?ShortName ?FromBus ?FromId ?FromGeoLocation ?ToBus ?ToId ?ToGeoLocation ?Resistance ?Reactance ?TotalLineChargingSusceptance
branch = sparql("select ?line ?ShortName ?FromBus ?FromId ?FromGeoLocation ?ToBus ?ToId ?ToGeoLocation ?Resistance ?Reactance ?TotalLineChargingSusceptance WHERE{
?line rdf:type cat:Electrical_Transmission.
?line prop:ShortName ?ShortName.
?line prop:FromBus ?FromBus.
?FromBus prop:BusIdentifier ?FromId.
?FromBus prop:GeoLocation ?FromGeoLocation.
?line prop:ToBus ?ToBus.
?ToBus prop:BusIdentifier ?ToId.
?ToBus prop:GeoLocation ?ToGeoLocation.
?line prop:Resistance ?Resistance.
?line prop:Reactance ?Reactance.
?line prop:TotalLineChargingSusceptance  ?TotalLineChargingSusceptance .
} ORDER BY ?line");

%build up the mpc.branch matrix
[m,n] =size(branch);
for i = 1:n

FBUS(i) = cell2mat(branch(4,i)); %F_Bus From Bus  (column 1)
TBUS(i) = cell2mat(branch(7,i)); %T_BUS To Bus  (column 2)
BR_R(i) = cell2mat(branch(9,i)); %BR_R Resistance  (column 3)
BR_X(i) = cell2mat(branch(10,i)); %BR_X Reactance  (column 4)
BR_B(i) = cell2mat(branch(11,i)); %BR_B Total line charging suseptance  (column 5)
RATE_A(i) = 1000; %RATE_A MVA rating A (column 6) !HARDCODED
RATE_B(i) = 1000; %RATE_B MVA rating B (column 7) !HARDCODED
RATE_C(i) = 1000; %RATE_C MVA rating C (column 8) !HARDCODED
TAP(i) = 0; %TAP Transformer off Nominal turns ratio (column 9) !HARDCODED
SHIFT(i) = 0; %SHIFT Transformer phase shift angles (column 10) !HARDCODED
BR_STATUS(i) = 1; %BR_STATUS Initial branch status (column 11) !HARDCODED
ANGMIN(i) = -360; %ANGMIN Minimum angle difference (column 12) !HARDCODED
ANGMAX(i) = 360; %ANDMAX Maximum angle difference (column 13) !HARDCODED

endfor

mpc.branch = [FBUS',TBUS',BR_R',BR_X',BR_B', RATE_A', RATE_B', RATE_C', TAP',SHIFT',BR_STATUS',ANGMIN',ANGMAX'];
mpc;

%----------------------------------------------------------------------------------------------------------------%
%Get the scenario data

%----------------------------------------------------------------------------------------------------------------%
%Evaluate the relevant case

SCENARIO = [ 1, 0.6, 1.05, 1.05, 1.05, 1.05; %SCENARIO MATRIX
1, 0.6, 1.05, 1.05, 1.05, 1.05;
1, 1, 1, 0, 1, 1;
1, 1, 1, 1, 0, 1;
1, 1, 1, 1, 1, 0 ];

SN = SCENARIO_NUMBERS(1,PLOTS)+1;

% Scenario 0 - No modification
% Scenario 1 - Base case (light load)
% Scenario 2 - Heavy load 
% Scenario 3 - Heavy loading plus outage of the lines 12-13 and 12-16 (Vales - Sydney North Vales - Munmorah)
% Scenario 4 - Heavy loading plus outage of generator #6 (Liddell)
% Scenario 5 - Heavy loading plus outage of the lines 9-19 and 9-20 (Bayswater - Regentville and Bayswater - Mt Piper)

mpc.bus(:,[3:4]) = SCENARIO(1,SN)*mpc.bus(:,[3:4]);
mpc.gen(:,[2]) = SCENARIO(2,SN)*mpc.gen(:,[2]);
mpc.branch([18,19],[11]) = SCENARIO(3,SN)*mpc.branch([18,19],[11]);
mpc.gen([2],[8]) = SCENARIO(4,SN)*mpc.gen([2],[8]);
mpc.branch([27,48],[11]) = SCENARIO(5,SN)*mpc.branch([27,48],[11]);
diary off;
results = runpf(mpc, mpoption( "OUT_ALL",0));
diary on;
%----------------------------------------------------------------------------------------------------------------%
%Apparent Power Calculation

PP_IN = diag(results.branch(:,14)*(results.branch(:,14))'); %Test Apparent power in is under limit
QQ_IN = diag(results.branch(:,15)*(results.branch(:,15))');
APPARENT_POWER_IN = sqrt (PP_IN+QQ_IN);

PP_OUT = diag(results.branch(:,16)*(results.branch(:,16))'); %Test Apparent power out is under limit
QQ_OUT = diag(results.branch(:,17)*(results.branch(:,17))');
APPARENT_POWER_OUT = sqrt (PP_OUT+QQ_OUT);

for t = 1: n
if APPARENT_POWER_IN(t,1) < 1000, APPARENT_POWER_TEST(t,1) =1; else APPARENT_POWER_TEST(t,1) = 0; endif; %Ensure BOTH power in AND power out is under limit. 
if APPARENT_POWER_OUT(t,1) < 1000, APPARENT_POWER_TEST(t,2) =1; else APPARENT_POWER_TEST(t,2) = 0; endif;
APPARENT_POWER_TEST(t,3) = APPARENT_POWER_TEST(t,2) + APPARENT_POWER_TEST(t,1);
endfor

%----------------------------------------------------------------------------------------------------------------%
%PLOTTING THE RESULTS
%----------------------------------------------------------------------------------------------------------------%


%Set Viewpoint characteristics (i.e. Lat, Long, zoom) ZOOM LEVEL MUST BE BETWEEN 6 & 8 



%V_P = [-33.75, 150.5, 8]; %Reasonable co-ordinates for zoom level 8 (Sydney level)
%V_P = [-33, 149, 7]; %Reasonable co-ordinates for zoom level 7 (North Coast level)
%V_P = [-33, 148, 6]; %Reasonable co-ordinates for zoom level 6 (NSW level)


%----------------------------------------------------------------------------------------------------------------%
%Get Image from Google Maps API

LAT = num2str(V_P(1,1));
LONG = num2str(V_P(1,2));
ZOOM = num2str(V_P(1,3));

CENTER = strcat (LAT,",",LONG);

mapbgf = urlread(strcat("http://maps.google.com/maps/api/staticmap?center=",CENTER,"&zoom=",ZOOM,"&size=640x640&sensor=false&style=feature:all",char(124),"visibility:off&style=feature:water",char(124),"saturation:-100",char(124),"lightness:-50",char(124),"visibility:on&style=feature:administrative.province",char(124),"element:geometry",char(124),"visibility:on&style=feature:landscape",char(124),"element:geometry",char(124),"lightness:100",char(124),"visibility:on"));
bgfid=fopen("static.png","w");
fwrite(bgfid,mapbgf,"char");
fclose(bgfid);

I = imread("static.png");

if P_N == 4, PLOT_GRID = 2; else PLOT_GRID = 1; endif
subplot(PLOT_GRID,PLOT_GRID,PLOTS); imshow(I)
hold on

%----------------------------------------------------------------------------------------------------------------%

%Zoom Level transform variables

T  = [ 6.05, 2.97, 1.452; %Transform matrix 
6.9, 3.5, 1.75;
45.87, 90.9, 183.15;
53.145, 108.67, 221.23 ];

Z_n = V_P(1,3) - 5; %Zoom Number (1,2 or 3, where 1 = "zoom level 6", 2 = "zoom level 7", 3 = "zoom level 8")


Y = -V_P (1,1) - T(1,Z_n); %translation in y
X = V_P(1,2) - T(2,Z_n); %translation in x
XD = T(3,Z_n); %dilation factor x
YD = T(4,Z_n); %dilation factor y


%----------------------------------------------------------------------------------------------------------------%

%Color scheme


%Ok lines- 4284d3 – 66,132,211, width 1 
LINES = [11,97,164]/256;
%Err lines - FFB540 – 255,181,64, width 2
ERRLINES = [255,191,0]/256;
%Out lines – 000000 – 0,0,0 – width 1, could we make these dashed?
OUTLINES = [255,49,0]/256;
%Ok bus - 04346C – 4, 52,108 – Radius 1 - Dot
BUS_COL = [3,62,107]/256;
%Err bus - A66500 – 166,101,0 – Radius 2 – Concentric circles
ERR_BUS_COL = [166,124,0]/256;

%----------------------------------------------------------------------------------------------------------------%
subplot(PLOT_GRID,PLOT_GRID,PLOTS);

%Plot Branches
for i = 1:n
c = [sscanf(char(branch(5,i)), '%f,'),sscanf(char(branch(8,i)), '%f,')];
if ((APPARENT_POWER_IN(i,1)) == 0),
plot (((c(2,:)-X)*XD),((Y+c(1,:))*-YD),"color",OUTLINES,"linewidth",2,"linestyle","--")
elseif (abs (APPARENT_POWER_TEST(i,3)) < 2), 
plot (((c(2,:)-X)*XD),((Y+c(1,:))*-YD),"color",ERRLINES,"linewidth",2)
else plot (((c(2,:)-X)*XD),((Y+c(1,:))*-YD),"color",LINES,"linewidth",1)
endif
endfor


%Plot Buses
for i = 1:B_n
GL = sscanf(char(bus(12,i)), '%f,');
if results.bus(i,8) > 1.1
for k = 1:3 plot (((GL(2,1)-X)*XD),((Y+GL(1,1))*-YD), "o","color",ERR_BUS_COL, "markersize",(2*k)) endfor
elseif results.bus(i,8) < 0.9 
for k = 1:3 plot (((GL(2,1)-X)*XD),((Y+GL(1,1))*-YD), "o","color",ERR_BUS_COL,"markersize",(2*k)) endfor
else for k = 1:3 plot (((GL(2,1)-X)*XD),((Y+GL(1,1))*-YD), "o","color",BUS_COL,"markersize",k) endfor
endif
endfor

%----------------------------------------------------------------------------------------------------------------%

%Plot Legend



L_L = 3; % Line length in legend
if P_N == 4, L_S = 20; else L_S = 10; endif; % Legend Spacing  
if P_N == 4, L_TS = 30; else L_TS = 15; endif; % Legend Title spacing

%Legend images
for k = 1:3 plot (L_P(1,1), L_P(1,2)+L_TS, "o","color",ERR_BUS_COL,"markersize",(2*k)) endfor
for k = 1:3 plot (L_P(1,1), L_P(1,2)+L_TS+L_S, "o","color",BUS_COL,"markersize",k) endfor
plot ([(L_P(1,1)-L_L),(L_P(1,1)+L_L)],[L_P(1,2)+L_TS+2*L_S,L_P(1,2)+L_TS+2*L_S],"color",OUTLINES,"linewidth",2)
plot ([(L_P(1,1)-L_L),(L_P(1,1)+L_L)],[L_P(1,2)+L_TS+3*L_S,L_P(1,2)+L_TS+3*L_S],"color",ERRLINES,"linewidth",2)
plot ([(L_P(1,1)-L_L),(L_P(1,1)+L_L)],[L_P(1,2)+L_TS+4*L_S,L_P(1,2)+L_TS+4*L_S],"color",LINES,"linewidth",1)

%Legend text
text (L_P(1,1), L_P(1,2), "LEGEND")
text (L_P(1,1)+10, L_P(1,2)+L_TS, "Bus Voltage Outside Limits")
text (L_P(1,1)+10, L_P(1,2)+L_TS+L_S, "Bus OK")
text (L_P(1,1)+10, L_P(1,2)+L_TS+2*L_S, "Line Outage")
text (L_P(1,1)+10, L_P(1,2)+L_TS+3*L_S, "Line Overload")
text (L_P(1,1)+10, L_P(1,2)+L_TS+4*L_S, "Line Ok")

%TITLE

SCEN_TITLE = strcat ("Scenario",{" "},num2str(SN-1));
text (300,0, SCEN_TITLE)

%----------------------------------------------------------------------------------------------------------------%
endfor

hold off
OperationCostAsPercentageOfFirstCost0.02info.png
 0.02
  +
ReactivePowerOutput0 MVAinfo.png
'0 MVA'
  +
RegionNew South Wales  +
ReplacementPeriod25 yrinfo.png
'25 yr'
  +
ShuntConductance0 MWinfo.png
'0 MW'
  +
VoltageAngle0info.png
0
  +