def test_pmu_with_trafo3w():
net = pp.create_empty_network()
bus_slack = pp.create_bus(net, vn_kv=110)
pp.create_ext_grid(net, bus=bus_slack)
bus_20_1 = pp.create_bus(net, vn_kv=20, name="b")
pp.create_sgen(net, bus=bus_20_1, p_mw=0.03, q_mvar=0.02)
bus_10_1 = pp.create_bus(net, vn_kv=10)
pp.create_sgen(net, bus=bus_10_1, p_mw=0.02, q_mvar=0.02)
bus_10_2 = pp.create_bus(net, vn_kv=10)
pp.create_load(net, bus=bus_10_2, p_mw=0.06, q_mvar=0.01)
pp.create_line(net,
from_bus=bus_10_1,
to_bus=bus_10_2,
std_type="149-AL1/24-ST1A 10.0",
length_km=2)
pp.create_transformer3w(net,
bus_slack,
bus_20_1,
bus_10_1,
std_type="63/25/38 MVA 110/20/10 kV")
pp.runpp(net)
add_virtual_pmu_meas_from_loadflow(net, with_random_error=False)
estimate(net, algorithm="lp", maximum_iterations=10)
pp.runpp(net)
assert np.allclose(net.res_bus.vm_pu, net.res_bus_est.vm_pu, atol=1e-2)
assert np.allclose(net.res_bus.va_degree,
net.res_bus_est.va_degree,
atol=1e-1)
def test_pmu_with_trafo3w():
net = pp.create_empty_network()
bus_slack = pp.create_bus(net, vn_kv=110)
pp.create_ext_grid(net, bus=bus_slack)
bus_20_1 = pp.create_bus(net, vn_kv=20, name="b")
pp.create_sgen(net, bus=bus_20_1, p_mw=0.03, q_mvar=0.02)
bus_10_1 = pp.create_bus(net, vn_kv=10)
pp.create_sgen(net, bus=bus_10_1, p_mw=0.02, q_mvar=0.02)
bus_10_2 = pp.create_bus(net, vn_kv=10)
pp.create_load(net, bus=bus_10_2, p_mw=0.06, q_mvar=0.01)
pp.create_line(net,
from_bus=bus_10_1,
to_bus=bus_10_2,
std_type="149-AL1/24-ST1A 10.0",
length_km=2)
pp.create_transformer3w(net,
bus_slack,
bus_20_1,
bus_10_1,
std_type="63/25/38 MVA 110/20/10 kV")
pp.runpp(net)
add_virtual_pmu_meas_from_loadflow(net, with_random_error=False)
run_se_lp_verify(net)
def test_pmu_case14():
net = nw.case9()
pp.runpp(net)
add_virtual_pmu_meas_from_loadflow(net)
run_se_lp_verify(net)
def test_pmu_case14():
net = nw.case14()
pp.runpp(net)
add_virtual_pmu_meas_from_loadflow(net)
estimate(net, algorithm="lp", maximum_iterations=20)
assert np.allclose(net.res_bus.vm_pu, net.res_bus_est.vm_pu, atol=1e-2)
assert np.allclose(net.res_bus.va_degree,
net.res_bus_est.va_degree,
atol=1e-1)