In the 2016 revision of the standard, the contribution of full converter elements like wind parks of photovoltaic or wind power plants, which was previously neglected, has been integrated.
In the paper Short Circuit Calculation with Fullsize Converters According to IEC 60909, the process is explained and an example is given. This example is used here to validate the correct implementation of the 2016 revision of the standard in pandapower.
The following network is introduced as example in the paper:
This network can be modeled in pandapower with the line parameters given in the paper as:
import pandapower as pp
import pandapower.shortcircuit as sc
import numpy as np
net = pp.create_empty_network()
for i in range(4):
pp.create_bus(net, vn_kv=110., index=i+1)
pp.create_ext_grid(net, 1, s_sc_max_mva=20*110*np.sqrt(3), rx_max=0.1)
pp.create_std_type(net, {"r_ohm_per_km": 0.120, "x_ohm_per_km": 0.393, "c_nf_per_km": 0.19,
"max_i_ka": 0.6}, "example_type")
for fb, tb, length in [(1, 2, 100), (1, 3, 50), (2, 3, 50), (3, 4, 25)]:
pp.create_line(net, from_bus=fb, to_bus=tb, length_km=length, std_type="example_type")
for b, p in [(2, 100), (3, 50), (4, 50)]:
pp.create_sgen(net, b, p_mw=p, sn_mva=p, k=1.2)
First, we disable the static generators that represent the wind parks and run a short circuit calculation without the wind park contribution:
net.sgen.in_service = False
sc.calc_sc(net, ip=True)
The result given in the paper without wind park contribution for a fault at bus 2 is Ikss = 2.9133 kA and ip = 5.9746 kA. Checking the results table shows that pandapower gives the same results:
net.res_bus_sc
ikss_ka | ip_ka | |
---|---|---|
1 | 20.000000 | 49.384390 |
2 | 2.913118 | 5.974240 |
3 | 3.706618 | 7.658568 |
4 | 2.399375 | 4.897167 |
Minor deviations in the results can be explained since the values in the paper are obtained from rounded intermediate results in the paper.
We now activate the three wind parks and run the short circuit calculation again:
net.sgen.in_service = True
sc.calc_sc(net, ip=True)
The result given in the paper with wind park contribution for a fault at bus 2 is Ikss = 3.9034 kA and ip = 7.3746 kA. Checking the results table shows that pandapower results yields the same result:
net.res_bus_sc
ikss_ka | ip_ka | |
---|---|---|
1 | 21.259673 | 51.165837 |
2 | 3.903126 | 7.374321 |
3 | 4.794773 | 9.197452 |
4 | 3.214790 | 6.050337 |