def add_test_impedance(net):
b1, b2, ln = add_grid_connection(net, zone="test_impedance")
b3 = pp.create_bus(net, vn_kv=220., zone="test_impedance")
rij = 0.02
xij = 0.01
rji = 0.03
xji = 0.005
s = 2.
pl = 1
ql = 0.5
pp.create_impedance(net,
b2,
b3,
rft_pu=rij,
xft_pu=xij,
rtf_pu=rji,
xtf_pu=xji,
sn_mva=s,
index=pp.get_free_id(net.impedance) + 1)
pp.create_impedance(net,
b2,
b3,
rft_pu=rij,
xft_pu=xij,
rtf_pu=rji,
xtf_pu=xji,
sn_mva=s,
index=pp.get_free_id(net.impedance) + 1,
in_service=False)
pp.create_load(net, b3, p_mw=pl, q_mvar=ql)
net.last_added_case = "test_impedance"
return net
def test_impedance_line_replacement():
# create test net
net1 = pp.create_empty_network(sn_kva=1.1e3)
pp.create_buses(net1, 2, 10)
pp.create_ext_grid(net1, 0)
pp.create_impedance(net1, 0, 1, 0.1, 0.1, 8.7)
pp.create_load(net1, 1, 7, 2)
# validate loadflow results
pp.runpp(net1)
net2 = copy.deepcopy(net1)
pp.replace_impedance_by_line(net2)
pp.runpp(net2)
assert pp.nets_equal(net1, net2,
exclude_elms={"line",
"impedance"}) # Todo: exclude_elms
cols = [
"p_from_kw", "q_from_kvar", "p_to_kw", "q_to_kvar", "pl_kw", "ql_kvar",
"i_from_ka", "i_to_ka"
]
assert np.allclose(net1.res_impedance[cols].values,
net2.res_line[cols].values)
net3 = copy.deepcopy(net2)
pp.replace_line_by_impedance(net3)
pp.runpp(net3)
assert pp.nets_equal(net2, net3, exclude_elms={"line", "impedance"})
assert np.allclose(net3.res_impedance[cols].values,
net2.res_line[cols].values)
def impedance_net():
net = pp.create_empty_network()
b1 = pp.create_bus(net, 220)
b2 = pp.create_bus(net, 30)
pp.create_ext_grid(net,
b1,
s_sc_max_mva=100.,
s_sc_min_mva=40.,
rx_min=0.1,
rx_max=0.1)
pp.create_ext_grid(net,
b2,
s_sc_max_mva=100.,
s_sc_min_mva=40.,
rx_min=0.1,
rx_max=0.1)
pp.create_impedance(net,
b1,
b2,
rft_pu=0.01,
xft_pu=0.02,
rtf_pu=0.05,
xtf_pu=0.01,
sn_kva=1e3)
return net
def test_create_replacement_switch_for_branch():
net = pp.create_empty_network()
bus0 = pp.create_bus(net, vn_kv=0.4)
bus1 = pp.create_bus(net, vn_kv=0.4)
bus2 = pp.create_bus(net, vn_kv=0.4)
bus3 = pp.create_bus(net, vn_kv=0.4)
pp.create_ext_grid(net, bus0, vm_pu=0.4)
line0 = pp.create_line(net,
bus0,
bus1,
length_km=1,
std_type="NAYY 4x50 SE")
line1 = pp.create_line(net,
bus2,
bus3,
length_km=1,
std_type="NAYY 4x50 SE")
impedance0 = pp.create_impedance(net, bus1, bus2, 0.01, 0.01, sn_kva=1e5)
impedance1 = pp.create_impedance(net, bus1, bus2, 0.01, 0.01, sn_kva=1e5)
pp.create_load(net, bus2, 1)
pp.runpp(net)
# look that the switch is created properly
tb.create_replacement_switch_for_branch(net, 'line', line0)
tb.create_replacement_switch_for_branch(net, 'impedance', impedance0)
net.line.in_service.at[line0] = False
net.impedance.in_service.at[impedance0] = False
assert 'REPLACEMENT_line_0' in net.switch.name.values
assert 'REPLACEMENT_impedance_0' in net.switch.name.values
assert net.switch.closed.at[0]
assert net.switch.closed.at[1]
pp.runpp(net)
# look that the switch is created with the correct closed status
net.line.in_service.at[line1] = False
net.impedance.in_service.at[impedance1] = False
tb.create_replacement_switch_for_branch(net, 'line', line1)
tb.create_replacement_switch_for_branch(net, 'impedance', impedance1)
assert 'REPLACEMENT_line_1' in net.switch.name.values
assert 'REPLACEMENT_impedance_1' in net.switch.name.values
assert ~net.switch.closed.at[2]
assert ~net.switch.closed.at[3]
def test_xward_buses():
"""
Issue: xward elements create dummy buses for the load flow, that are cleaned up afterwards.
However, if the load flow does not converge, those buses end up staying in the net and don't get
removed. This can potentially lead to thousands of dummy buses in net.
"""
net = pp.create_empty_network()
bus_sl = pp.create_bus(net, 110, name='ExtGrid')
pp.create_ext_grid(net, bus_sl, vm_pu=1)
bus_x = pp.create_bus(net, 110, name='XWARD')
pp.create_xward(net, bus_x, 0, 0, 0, 0, 0, 10, 1.1)
iid = pp.create_impedance(net, bus_sl, bus_x, 0.2, 0.2, 1e3)
bus_num1 = len(net.bus)
pp.runpp(net)
bus_num2 = len(net.bus)
assert bus_num1 == bus_num2
# now - make sure that the loadflow doesn't converge:
net.impedance.at[iid, 'rft_pu'] = 1
pp.create_load(net, bus_x, 1e6, 0)
with pytest.raises(LoadflowNotConverged):
# here the load flow doesn't converge and there is an extra bus in net
pp.runpp(net)
bus_num3 = len(net.bus)
assert bus_num3 == bus_num1
def net():
net = pp.create_empty_network()
bus0 = pp.create_bus(net, vn_kv=0.4)
bus1 = pp.create_bus(net, vn_kv=0.4)
bus2 = pp.create_bus(net, vn_kv=0.4)
bus3 = pp.create_bus(net, vn_kv=0.4)
bus4 = pp.create_bus(net, vn_kv=0.4)
bus5 = pp.create_bus(net, vn_kv=0.4)
bus6 = pp.create_bus(net, vn_kv=0.4)
pp.create_ext_grid(net, bus0, vm_pu=0.4)
line0 = pp.create_line(net,
bus0,
bus1,
length_km=0,
std_type="NAYY 4x50 SE")
line1 = pp.create_line_from_parameters(net,
bus2,
bus3,
length_km=1,
r_ohm_per_km=0,
x_ohm_per_km=0.1,
c_nf_per_km=0,
max_i_ka=1)
line2 = pp.create_line_from_parameters(net,
bus3,
bus4,
length_km=1,
r_ohm_per_km=0,
x_ohm_per_km=0,
c_nf_per_km=0,
max_i_ka=1)
impedance0 = pp.create_impedance(net, bus1, bus2, 0.01, 0.01, sn_kva=1e5)
impedance1 = pp.create_impedance(net, bus4, bus5, 0, 0, sn_kva=1e5)
impedance2 = pp.create_impedance(net,
bus5,
bus6,
0,
0,
rtf_pu=0.1,
sn_kva=1e5)
return net
def test_unsupplied_buses_with_switches():
net = pp.create_empty_network()
pp.create_buses(net, 8, 20)
pp.create_buses(net, 5, 0.4)
pp.create_ext_grid(net, 0)
pp.create_line(net, 0, 1, 1.2, "NA2XS2Y 1x185 RM/25 12/20 kV")
pp.create_switch(net, 0, 0, "l", closed=True)
pp.create_switch(net, 1, 0, "l", closed=False)
pp.create_line(net, 0, 2, 1.2, "NA2XS2Y 1x185 RM/25 12/20 kV")
pp.create_switch(net, 0, 1, "l", closed=False)
pp.create_switch(net, 2, 1, "l", closed=True)
pp.create_line(net, 0, 3, 1.2, "NA2XS2Y 1x185 RM/25 12/20 kV")
pp.create_switch(net, 0, 2, "l", closed=False)
pp.create_switch(net, 3, 2, "l", closed=False)
pp.create_line(net, 0, 4, 1.2, "NA2XS2Y 1x185 RM/25 12/20 kV")
pp.create_switch(net, 0, 3, "l", closed=True)
pp.create_switch(net, 4, 3, "l", closed=True)
pp.create_line(net, 0, 5, 1.2, "NA2XS2Y 1x185 RM/25 12/20 kV")
pp.create_switch(net, 0, 6, "b", closed=True)
pp.create_switch(net, 0, 7, "b", closed=False)
pp.create_transformer(net, 0, 8, "0.63 MVA 20/0.4 kV")
pp.create_switch(net, 0, 0, "t", closed=True)
pp.create_switch(net, 8, 0, "t", closed=False)
pp.create_transformer(net, 0, 9, "0.63 MVA 20/0.4 kV")
pp.create_switch(net, 0, 1, "t", closed=False)
pp.create_switch(net, 9, 1, "t", closed=True)
pp.create_transformer(net, 0, 10, "0.63 MVA 20/0.4 kV")
pp.create_switch(net, 0, 2, "t", closed=False)
pp.create_switch(net, 10, 2, "t", closed=False)
pp.create_transformer(net, 0, 11, "0.63 MVA 20/0.4 kV")
pp.create_switch(net, 0, 3, "t", closed=True)
pp.create_switch(net, 11, 3, "t", closed=True)
pp.create_transformer(net, 0, 12, "0.63 MVA 20/0.4 kV")
pp.create_buses(net, 2, 20)
pp.create_impedance(net, 0, 13, 1, 1, 10)
pp.create_impedance(net, 0, 14, 1, 1, 10, in_service=False)
ub = top.unsupplied_buses(net)
assert ub == {1, 2, 3, 7, 8, 9, 10, 14}
ub = top.unsupplied_buses(net, respect_switches=False)
assert ub == {14}
def add_test_impedance(net):
b1, b2, ln = add_grid_connection(net, zone="test_impedance")
b3 = pp.create_bus(net, vn_kv=220., zone="test_impedance")
r = 0.02
x = 0.01
s = 2000
pl = 1000
ql = 500
pp.create_impedance(net,
b2,
b3,
r_pu=r,
x_pu=x,
sn_kva=s,
index=pp.get_free_id(net.impedance) + 1)
pp.create_load(net, b3, p_kw=pl, q_kvar=ql)
net.last_added_case = "test_impedance"
return net
def example_multivoltage():
"""
Returns the multivoltage example network from the pandapower tutorials.
OUTPUT:
net - multivoltage example network
EXAMPLE:
>>> import pandapower.networks
>>> net = pandapower.networks.example_multivoltage()
"""
net = pp.create_empty_network()
# --- Busses
# HV
# Double busbar
pp.create_bus(net, name='Double Busbar 1', vn_kv=380, type='b')
pp.create_bus(net, name='Double Busbar 2', vn_kv=380, type='b')
for i in range(10):
pp.create_bus(net, name='Bus DB T%s' % i, vn_kv=380, type='n')
for i in range(1, 5):
pp.create_bus(net, name='Bus DB %s' % i, vn_kv=380, type='n')
# Single busbar
pp.create_bus(net, name='Single Busbar', vn_kv=110, type='b')
for i in range(1, 6):
pp.create_bus(net, name='Bus SB %s' % i, vn_kv=110, type='n')
for i in range(1, 6):
for j in [1, 2]:
pp.create_bus(net,
name='Bus SB T%s.%s' % (i, j),
vn_kv=110,
type='n')
# Remaining
for i in range(1, 5):
pp.create_bus(net, name='Bus HV%s' % i, vn_kv=110, type='n')
# MV
pp.create_bus(net, name='Bus MV0 20kV', vn_kv=20, type='n')
for i in range(8):
pp.create_bus(net, name='Bus MV%s' % i, vn_kv=10, type='n')
# LV
pp.create_bus(net, name='Bus LV0', vn_kv=0.4, type='n')
for i in range(1, 6):
pp.create_bus(net, name='Bus LV1.%s' % i, vn_kv=0.4, type='m')
for i in range(1, 5):
pp.create_bus(net, name='Bus LV2.%s' % i, vn_kv=0.4, type='m')
pp.create_bus(net, name='Bus LV2.2.1', vn_kv=0.4, type='m')
pp.create_bus(net, name='Bus LV2.2.2', vn_kv=0.4, type='m')
# --- Lines
# HV
hv_lines = pd.DataFrame()
hv_lines['line_name'] = ['HV Line%s' % i for i in range(1, 7)]
hv_lines['from_bus'] = [
'Bus SB 2', 'Bus HV1', 'Bus HV2', 'Bus HV1', 'Bus HV3', 'Bus SB 3'
]
hv_lines['to_bus'] = [
'Bus HV1', 'Bus HV2', 'Bus HV4', 'Bus HV4', 'Bus HV4', 'Bus HV3'
]
hv_lines['std_type'] = '184-AL1/30-ST1A 110.0'
hv_lines['length'] = [30, 20, 30, 15, 25, 30]
hv_lines['parallel'] = [1, 1, 1, 1, 1, 2]
for _, hv_line in hv_lines.iterrows():
from_bus = pp.get_element_index(net, "bus", hv_line.from_bus)
to_bus = pp.get_element_index(net, "bus", hv_line.to_bus)
pp.create_line(net,
from_bus,
to_bus,
length_km=hv_line.length,
std_type=hv_line.std_type,
name=hv_line.line_name,
parallel=hv_line.parallel)
# MV
mv_lines = pd.DataFrame()
mv_lines['line_name'] = ['MV Line%s' % i for i in range(1, 9)]
mv_lines['from_bus'] = ['Bus MV%s' % i for i in list(range(7)) + [0]]
mv_lines['to_bus'] = ['Bus MV%s' % i for i in list(range(1, 8)) + [7]]
mv_lines['length'] = 1.5
mv_lines['std_type'] = 'NA2XS2Y 1x185 RM/25 12/20 kV'
for _, mv_line in mv_lines.iterrows():
from_bus = pp.get_element_index(net, "bus", mv_line.from_bus)
to_bus = pp.get_element_index(net, "bus", mv_line.to_bus)
pp.create_line(net,
from_bus,
to_bus,
length_km=mv_line.length,
std_type=mv_line.std_type,
name=mv_line.line_name)
# LV
lv_lines = pd.DataFrame()
lv_line_idx = [
'1.1', '1.2', '1.3', '1.4', '1.6', '2.1', '2.2', '2.3', '2.4', '2.2.1',
'2.2.2'
]
lv_lines['line_name'] = ['LV Line%s' % i for i in lv_line_idx]
lv_line_idx = [
'0', '1.1', '1.2', '1.3', '1.4', '0', '2.1', '2.2', '2.3', '2.2',
'2.2.1'
]
lv_lines['from_bus'] = ['Bus LV%s' % i for i in lv_line_idx]
lv_line_idx = [
'1.1', '1.2', '1.3', '1.4', '1.5', '2.1', '2.2', '2.3', '2.4', '2.2.1',
'2.2.2'
]
lv_lines['to_bus'] = ['Bus LV%s' % i for i in lv_line_idx]
lv_lines['length'] = [0.08] * 5 + [0.12] * 6
lv_lines['std_type'] = ['NAYY 4x120 SE'] * 7 + ['15-AL1/3-ST1A 0.4'] * 4
for _, lv_line in lv_lines.iterrows():
from_bus = pp.get_element_index(net, "bus", lv_line.from_bus)
to_bus = pp.get_element_index(net, "bus", lv_line.to_bus)
pp.create_line(net,
from_bus,
to_bus,
length_km=lv_line.length,
std_type=lv_line.std_type,
name=lv_line.line_name)
# --- Transformer
hv_bus = pp.get_element_index(net, "bus", "Bus DB 2")
lv_bus = pp.get_element_index(net, "bus", "Bus SB 1")
pp.create_transformer_from_parameters(net,
hv_bus,
lv_bus,
sn_kva=300000,
vn_hv_kv=380,
vn_lv_kv=110,
vscr_percent=0.06,
vsc_percent=8,
pfe_kw=0,
i0_percent=0,
tp_pos=0,
shift_degree=0,
name='EHV-HV-Trafo')
hv_bus = pp.get_element_index(net, "bus", "Bus MV4")
lv_bus = pp.get_element_index(net, "bus", "Bus LV0")
pp.create_transformer_from_parameters(net,
hv_bus,
lv_bus,
sn_kva=400,
vn_hv_kv=10,
vn_lv_kv=0.4,
vscr_percent=1.325,
vsc_percent=4,
pfe_kw=0.95,
i0_percent=0.2375,
tp_side="hv",
tp_mid=0,
tp_min=-2,
tp_max=2,
tp_st_percent=2.5,
tp_pos=0,
shift_degree=150,
name='MV-LV-Trafo')
# Trafo3w
hv_bus = pp.get_element_index(net, "bus", "Bus HV2")
mv_bus = pp.get_element_index(net, "bus", "Bus MV0 20kV")
lv_bus = pp.get_element_index(net, "bus", "Bus MV0")
pp.create_transformer3w_from_parameters(net,
hv_bus,
mv_bus,
lv_bus,
vn_hv_kv=110,
vn_mv_kv=20,
vn_lv_kv=10,
sn_hv_kva=40000,
sn_mv_kva=15000,
sn_lv_kva=25000,
vsc_hv_percent=10.1,
vsc_mv_percent=10.1,
vsc_lv_percent=10.1,
vscr_hv_percent=0.266667,
vscr_mv_percent=0.033333,
vscr_lv_percent=0.04,
pfe_kw=0,
i0_percent=0,
shift_mv_degree=30,
shift_lv_degree=30,
tp_side="hv",
tp_mid=0,
tp_min=-8,
tp_max=8,
tp_st_percent=1.25,
tp_pos=0,
name='HV-MV-MV-Trafo')
# --- Static generators
# HV
pp.create_sgen(net,
pp.get_element_index(net, "bus", 'Bus SB 5'),
p_kw=-20000,
q_kvar=-4000,
sn_kva=45000,
type='WP',
name='Wind Park')
# MV
mv_sgens = pd.DataFrame()
mv_sgens['sgen_name'] = [
'Biogas plant', 'Further MV Generator', 'Industry Generator', 'PV Park'
]
mv_sgens['bus'] = ['Bus MV6', 'Bus MV0', 'Bus MV0 20kV', 'Bus MV5']
mv_sgens['p'] = [-500, -500, -15000, -2000]
mv_sgens['q'] = [0, -50, -3000, -100]
mv_sgens['sn'] = [750, 1000, 20000, 5000]
mv_sgens['type'] = ['SGEN', 'SGEN', 'SGEN', 'PV']
for _, sgen in mv_sgens.iterrows():
bus_idx = pp.get_element_index(net, "bus", sgen.bus)
pp.create_sgen(net,
bus_idx,
p_kw=sgen.p,
q_kvar=sgen.q,
sn_kva=sgen.sn,
type=sgen.type,
name=sgen.sgen_name)
# LV
lv_sgens = pd.DataFrame()
lv_sgens['sgen_name'] = ['PV'] + ['PV(%s)' % i for i in range(1, 6)]
lv_sgens['bus'] = [
'Bus LV%s' % i for i in ['1.1', '1.3', '2.3', '2.4', '2.2.1', '2.2.2']
]
lv_sgens['p'] = [-6, -5, -5, -5, -5, -5]
lv_sgens['q'] = 0
lv_sgens['sn'] = [12, 10, 10, 10, 10, 10]
lv_sgens['type'] = 'PV'
for _, sgen in lv_sgens.iterrows():
bus_idx = pp.get_element_index(net, "bus", sgen.bus)
pp.create_sgen(net,
bus_idx,
p_kw=sgen.p,
q_kvar=sgen.q,
sn_kva=sgen.sn,
type=sgen.type,
name=sgen.sgen_name)
# --- Loads
# HV
hv_loads = pd.DataFrame()
hv_loads['load_name'] = ['MV Net %s' % i for i in range(5)]
hv_loads['bus'] = ['Bus SB 4', 'Bus HV1', 'Bus HV2', 'Bus HV3', 'Bus HV4']
hv_loads['p'] = 38000
hv_loads['q'] = 6000
for _, load in hv_loads.iterrows():
bus_idx = pp.get_element_index(net, "bus", load.bus)
pp.create_load(net,
bus_idx,
p_kw=load.p,
q_kvar=load.q,
name=load.load_name)
# MV
mv_loads = pd.DataFrame()
mv_loads['load_name'] = ['Further MV-Rings', 'Industry Load'
] + ['LV Net %s' % i for i in [1, 2, 3, 5, 6, 7]]
mv_loads['bus'] = ['Bus MV0', 'Bus MV0 20kV'
] + ['Bus MV%s' % i for i in [1, 2, 3, 5, 6, 7]]
mv_loads['p'] = [6000, 18000, 400, 400, 400, 400, 400, 400]
mv_loads['q'] = [2000, 4000, 100, 60, 60, 60, 60, 60]
for _, load in mv_loads.iterrows():
bus_idx = pp.get_element_index(net, "bus", load.bus)
pp.create_load(net,
bus_idx,
p_kw=load.p,
q_kvar=load.q,
name=load.load_name)
# LV
lv_loads = pd.DataFrame()
idx = ['', '(1)', '(2)', '(3)', '(4)', '(5)']
lv_loads['load_name'] = ['Further LV-Feeders Load'] + [
'Residential Load%s' % i for i in idx[0:5]
] + ['Rural Load%s' % i for i in idx[0:6]]
lv_loads['bus'] = [
'Bus LV%s' % i for i in [
'0', '1.1', '1.2', '1.3', '1.4', '1.5', '2.1', '2.2', '2.3', '2.4',
'2.2.1', '2.2.2'
]
]
lv_loads['p'] = [100] + [10] * 11
lv_loads['q'] = [10] + [3] * 11
for _, load in lv_loads.iterrows():
bus_idx = pp.get_element_index(net, "bus", load.bus)
pp.create_load(net,
bus_idx,
p_kw=load.p,
q_kvar=load.q,
name=load.load_name)
# --- Other
# Shunt
pp.create_shunt(net,
pp.get_element_index(net, "bus", 'Bus HV1'),
p_kw=0,
q_kvar=-960,
name='Shunt')
# ExtGrids
pp.create_ext_grid(net,
pp.get_element_index(net, "bus", 'Double Busbar 1'),
vm_pu=1.03,
va_degree=0,
name='External grid',
s_sc_max_mva=10000,
rx_max=0.1,
rx_min=0.1)
# Gen
pp.create_gen(net,
pp.get_element_index(net, "bus", 'Bus HV4'),
vm_pu=1.03,
p_kw=-1e5,
name='Gas turbine')
# Impedance
pp.create_impedance(net,
pp.get_element_index(net, "bus", 'Bus HV3'),
pp.get_element_index(net, "bus", 'Bus HV1'),
rft_pu=0.074873,
xft_pu=0.198872,
sn_kva=100000,
name='Impedance')
# xwards
pp.create_xward(net,
pp.get_element_index(net, "bus", 'Bus HV3'),
ps_kw=23942,
qs_kvar=-12241.87,
pz_kw=2814.571,
qz_kvar=0,
r_ohm=0,
x_ohm=12.18951,
vm_pu=1.02616,
name='XWard 1')
pp.create_xward(net,
pp.get_element_index(net, "bus", 'Bus HV1'),
ps_kw=3776,
qs_kvar=-7769.979,
pz_kw=9174.917,
qz_kvar=0,
r_ohm=0,
x_ohm=50.56217,
vm_pu=1.024001,
name='XWard 2')
# --- Switches
# HV
# Bus-bus switches
hv_bus_sw = pd.DataFrame()
hv_bus_sw['bus_name'] = ['DB DS%s' % i for i in range(14)] + \
['DB CB%s' % i for i in range(5)] + \
['SB DS%s.%s' % (i, j) for i in range(1, 6) for j in range(1, 3)] + \
['SB CB%s' % i for i in range(1, 6)]
hv_bus_sw['from_bus'] = ['Double Busbar %s' % i for i in [2, 1, 2, 1, 2, 1, 2, 1, 2, 1]] + \
['Bus DB T%s' % i for i in [2, 4, 6, 8, 0, 3, 5, 7, 9]] + \
['Bus SB T1.1', 'Single Busbar', 'Bus SB T2.1', 'Single Busbar',
'Bus SB T3.1', 'Single Busbar', 'Bus SB T4.1', 'Single Busbar',
'Bus SB T5.1', 'Single Busbar'] + \
['Bus SB T%s.2' % i for i in range(1, 6)]
hv_bus_sw['to_bus'] = ['Bus DB %s' % i for i in
['T0', 'T1', 'T3', 'T3', 'T5', 'T5', 'T7', 'T7', 'T9', 'T9',
'1', '2', '3', '4', 'T1', 'T2', 'T4', 'T6', 'T8']] + \
['Bus SB %s' % i for i in
['1', 'T1.2', '2', 'T2.2', '3', 'T3.2', '4', 'T4.2', '5', 'T5.2']] + \
['Bus SB T%s.1' % i for i in range(1, 6)]
hv_bus_sw['type'] = ['DS'] * 14 + ['CB'] * 5 + ['DS'] * 10 + ['CB'] * 5
hv_bus_sw['et'] = 'b'
hv_bus_sw['closed'] = [
bool(i)
for i in [1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1] +
[1] * 15
]
for _, switch in hv_bus_sw.iterrows():
from_bus = pp.get_element_index(net, "bus", switch.from_bus)
to_bus = pp.get_element_index(net, "bus", switch.to_bus)
pp.create_switch(net,
from_bus,
to_bus,
et=switch.et,
closed=switch.closed,
type=switch.type,
name=switch.bus_name)
# Bus-Line switches
hv_buses = net.bus[(net.bus.vn_kv == 380) | (net.bus.vn_kv == 110)].index
hv_ls = net.line[(net.line.from_bus.isin(hv_buses))
& (net.line.to_bus.isin(hv_buses))]
for _, line in hv_ls.iterrows():
for bus in [line.from_bus, line.to_bus]:
pp.create_switch(net,
bus,
line.name,
et='l',
closed=True,
type='LBS',
name='Switch %s - %s' %
(net.bus.name.at[bus], line['name']))
# MV
# Bus-line switches
mv_buses = net.bus[(net.bus.vn_kv == 10) | (net.bus.vn_kv == 20)].index
mv_ls = net.line[(net.line.from_bus.isin(mv_buses))
& (net.line.to_bus.isin(mv_buses))]
for _, line in mv_ls.iterrows():
for bus in [line.from_bus, line.to_bus]:
pp.create_switch(net,
bus,
line.name,
et='l',
closed=True,
type='LBS',
name='Switch %s - %s' %
(net.bus.name.at[bus], line['name']))
open_switch_id = net.switch[(
net.switch.name == 'Switch Bus MV5 - MV Line5')].index
net.switch.closed.loc[open_switch_id] = False
# LV
# Bus-line switches
lv_buses = net.bus[net.bus.vn_kv == 0.4].index
lv_ls = net.line[(net.line.from_bus.isin(lv_buses))
& (net.line.to_bus.isin(lv_buses))]
for _, line in lv_ls.iterrows():
for bus in [line.from_bus, line.to_bus]:
pp.create_switch(net,
bus,
line.name,
et='l',
closed=True,
type='LBS',
name='Switch %s - %s' %
(net.bus.name.at[bus], line['name']))
# Trafoswitches
# HV
pp.create_switch(net,
pp.get_element_index(net, "bus", 'Bus DB 2'),
pp.get_element_index(net, "trafo", 'EHV-HV-Trafo'),
et='t',
closed=True,
type='LBS',
name='Switch DB2 - EHV-HV-Trafo')
pp.create_switch(net,
pp.get_element_index(net, "bus", 'Bus SB 1'),
pp.get_element_index(net, "trafo", 'EHV-HV-Trafo'),
et='t',
closed=True,
type='LBS',
name='Switch SB1 - EHV-HV-Trafo')
# LV
pp.create_switch(net,
pp.get_element_index(net, "bus", 'Bus MV4'),
pp.get_element_index(net, "trafo", 'MV-LV-Trafo'),
et='t',
closed=True,
type='LBS',
name='Switch MV4 - MV-LV-Trafo')
pp.create_switch(net,
pp.get_element_index(net, "bus", 'Bus LV0'),
pp.get_element_index(net, "trafo", 'MV-LV-Trafo'),
et='t',
closed=True,
type='LBS',
name='Switch LV0 - MV-LV-Trafo')
# --- Powerflow
# run power flow and generate result tables
pp.runpp(net, init='dc', calculate_voltage_angles=True, Numba=False)
return net
def test_replace_zero_branches_with_switches():
net = pp.create_empty_network()
bus0 = pp.create_bus(net, vn_kv=0.4)
bus1 = pp.create_bus(net, vn_kv=0.4)
bus2 = pp.create_bus(net, vn_kv=0.4)
bus3 = pp.create_bus(net, vn_kv=0.4)
bus4 = pp.create_bus(net, vn_kv=0.4)
bus5 = pp.create_bus(net, vn_kv=0.4)
bus6 = pp.create_bus(net, vn_kv=0.4)
pp.create_ext_grid(net, bus0, vm_pu=0.4)
line0 = pp.create_line(net,
bus0,
bus1,
length_km=0,
std_type="NAYY 4x50 SE")
line1 = pp.create_line_from_parameters(net,
bus2,
bus3,
length_km=1,
r_ohm_per_km=0,
x_ohm_per_km=0.1,
c_nf_per_km=0,
max_i_ka=1)
line2 = pp.create_line_from_parameters(net,
bus3,
bus4,
length_km=1,
r_ohm_per_km=0,
x_ohm_per_km=0,
c_nf_per_km=0,
max_i_ka=1)
impedance0 = pp.create_impedance(net, bus1, bus2, 0.01, 0.01, sn_kva=1e5)
impedance1 = pp.create_impedance(net, bus4, bus5, 0, 0, sn_kva=1e5)
impedance2 = pp.create_impedance(net,
bus5,
bus6,
0,
0,
rtf_pu=0.1,
sn_kva=1e5)
# test for line with zero length
tb.replace_zero_branches_with_switches(net,
elements=('line', ),
zero_length=True,
zero_impedance=False,
in_service_only=True)
assert 'REPLACEMENT_line_0' in net.switch.name.values
assert ~net.line.in_service.at[0]
# test for in_service_only
tb.replace_zero_branches_with_switches(net,
elements=('line', ),
zero_length=True,
zero_impedance=False,
in_service_only=True)
assert len(net.switch.loc[net.switch.name == 'REPLACEMENT_line_0']) == 1
tb.replace_zero_branches_with_switches(net,
elements=('line', ),
zero_length=True,
zero_impedance=False,
in_service_only=False)
assert len(net.switch.loc[net.switch.name == 'REPLACEMENT_line_0']) == 2
assert ~net.switch.closed.at[1]
# test for line with zero impedance
tb.replace_zero_branches_with_switches(net,
elements=('line', ),
zero_length=False,
zero_impedance=True,
in_service_only=True)
assert 'REPLACEMENT_line_1' not in net.switch.name.values
assert 'REPLACEMENT_line_2' in net.switch.name.values
# test for impedance
tb.replace_zero_branches_with_switches(net,
elements=('impedance', ),
zero_length=False,
zero_impedance=True,
in_service_only=True)
assert 'REPLACEMENT_impedance_0' not in net.switch.name.values
assert 'REPLACEMENT_impedance_1' in net.switch.name.values
assert 'REPLACEMENT_impedance_2' not in net.switch.name.values
# now run everything altogether
net.switch.drop(net.switch.index, inplace=True)
net.line.loc[:, 'in_service'] = True
net.impedance.loc[:, 'in_service'] = True
tb.replace_zero_branches_with_switches(net,
elements=('impedance', 'line'),
zero_length=True,
zero_impedance=True,
in_service_only=True)
assert 'REPLACEMENT_impedance_1' in net.switch.name.values
assert 'REPLACEMENT_line_0' in net.switch.name.values
assert 'REPLACEMENT_line_2' in net.switch.name.values
assert ~net.line.in_service.at[0]
assert net.line.in_service.at[1]
assert ~net.line.in_service.at[2]
assert 'REPLACEMENT_impedance_0' not in net.switch.name.values
assert 'REPLACEMENT_impedance_2' not in net.switch.name.values
assert 'REPLACEMENT_line_1' not in net.switch.name.values
assert net.impedance.in_service.at[0]
assert ~net.impedance.in_service.at[1]
assert net.impedance.in_service.at[2]
net.sgen
pp.create_shunt(net,
pp.get_element_index(net, "bus", 'Bus HV1'),
p_kw=0,
q_kvar=-960,
name='Shunt')
# show shunt table
net.shunt
# Impedance
pp.create_impedance(net,
pp.get_element_index(net, "bus", 'Bus HV3'),
pp.get_element_index(net, "bus", 'Bus HV1'),
rft_pu=0.074873,
xft_pu=0.198872,
sn_kva=100000,
name='Impedance')
# show impedance table
net.impedance
# xwards
pp.create_xward(net,
pp.get_element_index(net, "bus", 'Bus HV3'),
ps_kw=23942,
qs_kvar=-12241.87,
pz_kw=2814.571,
qz_kvar=0,
r_ohm=0,
