From fc5bfcb5082d3cfff7ac29215b980c9920b6a694 Mon Sep 17 00:00:00 2001 From: Maryam Majidi Date: Fri, 12 Nov 2021 19:23:59 +0100 Subject: [PATCH] delete pycodes --- test/pycodes/cigre_timeseries_15min.json | 1 - test/pycodes/create_net.py | 124 ----------------------- test/pycodes/storage_test_case.py | 70 ------------- 3 files changed, 195 deletions(-) delete mode 100644 test/pycodes/cigre_timeseries_15min.json delete mode 100644 test/pycodes/create_net.py delete mode 100644 test/pycodes/storage_test_case.py diff --git a/test/pycodes/cigre_timeseries_15min.json b/test/pycodes/cigre_timeseries_15min.json deleted file mode 100644 index 31e0641..0000000 --- a/test/pycodes/cigre_timeseries_15min.json +++ /dev/null @@ -1 +0,0 @@ -{"timestep":{"0":0,"1":1,"2":2,"3":3,"4":4,"5":5,"6":6,"7":7,"8":8,"9":9,"10":10,"11":11,"12":12,"13":13,"14":14,"15":15,"16":16,"17":17,"18":18,"19":19,"20":20,"21":21,"22":22,"23":23,"24":24,"25":25,"26":26,"27":27,"28":28,"29":29,"30":30,"31":31,"32":32,"33":33,"34":34,"35":35,"36":36,"37":37,"38":38,"39":39,"40":40,"41":41,"42":42,"43":43,"44":44,"45":45,"46":46,"47":47,"48":48,"49":49,"50":50,"51":51,"52":52,"53":53,"54":54,"55":55,"56":56,"57":57,"58":58,"59":59,"60":60,"61":61,"62":62,"63":63,"64":64,"65":65,"66":66,"67":67,"68":68,"69":69,"70":70,"71":71,"72":72,"73":73,"74":74,"75":75,"76":76,"77":77,"78":78,"79":79,"80":80,"81":81,"82":82,"83":83,"84":84,"85":85,"86":86,"87":87,"88":88,"89":89,"90":90,"91":91,"92":92,"93":93,"94":94,"95":95},"pv":{"0":0.0061428571,"1":0.0058772894,"2":0.0055054945,"3":0.0052930403,"4":0.0049212454,"5":0.0046190476,"6":0.0042857143,"7":0.0039285714,"8":0.0036739927,"9":0.0033919414,"10":0.003029304,"11":0.0027472527,"12":0.0025714286,"13":0.0022857143,"14":0.0019952381,"15":0.0017333333,"16":0.0014908425,"17":0.0012600733,"18":0.0010555556,"19":0.0008571429,"20":0.0006984127,"21":0.0005311355,"22":0.0003553114,"23":0.0002222222,"24":0.0001318681,"25":0.0000659341,"26":0.000010989,"27":0.0,"28":0.0000612245,"29":0.000244898,"30":0.0003673469,"31":0.0006530612,"32":0.0010340136,"33":0.0017063492,"34":0.0024365079,"35":0.0220952381,"36":0.0538095238,"37":0.1234285714,"38":0.1466285714,"39":0.3191666667,"40":0.4307142857,"41":0.5313333333,"42":0.5887142857,"43":0.6539365079,"44":0.5895714286,"45":0.5693015873,"46":0.4239333333,"47":0.3119166667,"48":0.335,"49":0.4423142857,"50":0.4685238095,"51":0.4832619048,"52":0.4542952381,"53":0.4389428571,"54":0.2424603175,"55":0.148247619,"56":0.1802619048,"57":0.1397619048,"58":0.1086380952,"59":0.1067509158,"60":0.1093150183,"61":0.0926785714,"62":0.0563988095,"63":0.0353142857,"64":0.0254761905,"65":0.0232698413,"66":0.0222142857,"67":0.0211587302,"68":0.0203650794,"69":0.0195714286,"70":0.0189761905,"71":0.0175873016,"72":0.0168333333,"73":0.0162380952,"74":0.015452381,"75":0.0148333333,"76":0.0141836735,"77":0.0138435374,"78":0.0131031746,"79":0.0125277778,"80":0.0118359788,"81":0.0113703704,"82":0.0106709957,"83":0.0100865801,"84":0.0097359307,"85":0.0090857143,"86":0.0086857143,"87":0.008152381,"88":0.0073445378,"89":0.0067955182,"90":0.0063837535,"91":0.0058347339,"92":0.0054229692,"93":0.0047037037,"94":0.0041428571,"95":0.0036666667},"wind":{"0":0.6226,"1":0.6235846154,"2":0.6250564103,"3":0.6258974359,"4":0.6273692308,"5":0.6451384615,"6":0.6851282051,"7":0.7132666667,"8":0.7194974359,"9":0.6861128205,"10":0.6508,"11":0.6307333333,"12":0.6394,"13":0.6826,"14":0.7283733333,"15":0.77724,"16":0.8079589744,"17":0.7493384615,"18":0.6996555556,"19":0.6606,"20":0.6293555556,"21":0.5971948718,"22":0.5641179487,"23":0.5404555556,"24":0.5353589744,"25":0.5531128205,"26":0.5679076923,"27":0.5868666667,"28":0.6072285714,"29":0.5883142857,"30":0.5757047619,"31":0.5522571429,"32":0.5317904762,"33":0.5513666667,"34":0.5723,"35":0.6030333333,"36":0.6208,"37":0.6304666667,"38":0.6375866667,"39":0.6457,"40":0.6276,"41":0.5918888889,"42":0.5655666667,"43":0.5334222222,"44":0.5332666667,"45":0.5890444444,"46":0.7842666667,"47":0.9233333333,"48":1.0052166667,"49":1.0118133333,"50":0.9797333333,"51":0.9418333333,"52":0.90148,"53":0.8783333333,"54":0.8383555556,"55":0.81672,"56":0.7973,"57":0.7880444444,"58":0.7884,"59":0.7939538462,"60":0.8028564103,"61":0.80865,"62":0.8126083333,"63":0.81268,"64":0.8094333333,"65":0.75,"66":0.7076666667,"67":0.6527777778,"68":0.6118888889,"69":0.571,"70":0.5403333333,"71":0.4687777778,"72":0.48895,"73":0.5566444444,"74":0.6610666667,"75":0.7616,"76":0.754,"77":0.6961333333,"78":0.582,"79":0.5136666667,"80":0.4635407407,"81":0.5363703704,"82":0.6475151515,"83":0.7433030303,"84":0.8007757576,"85":0.8503866667,"86":0.8732266667,"87":0.90368,"88":0.932172549,"89":0.9282509804,"90":0.9253098039,"91":0.9213882353,"92":0.9184470588,"93":0.886837037,"94":0.8583393939,"95":0.838369697},"residential":{"0":0.256022,"1":0.248835,"2":0.2119733333,"3":0.201825,"4":0.209164,"5":0.1986,"6":0.18905,"7":0.186924,"8":0.1870066667,"9":0.1836575,"10":0.181048,"11":0.1855666667,"12":0.1822822222,"13":0.192295,"14":0.19936,"15":0.1822875,"16":0.186546,"17":0.192088,"18":0.203122,"19":0.22803,"20":0.2446333333,"21":0.240875,"22":0.26249,"23":0.3318133333,"24":0.43635,"25":0.501545,"26":0.52612,"27":0.53715,"28":0.614585,"29":0.6310525,"30":0.6089275,"31":0.60786,"32":0.624715,"33":0.670312,"34":0.69094,"35":0.64048,"36":0.6392,"37":0.636695,"38":0.652735,"39":0.65687,"40":0.678112,"41":0.6487575,"42":0.611575,"43":0.62915,"44":0.60457,"45":0.6140533333,"46":0.702965,"47":0.726462,"48":0.7271,"49":0.7306,"50":0.71867,"51":0.6943366667,"52":0.6500266667,"53":0.632864,"54":0.5709366667,"55":0.5649275,"56":0.565415,"57":0.5238728571,"58":0.49257,"59":0.49336,"60":0.494898,"61":0.48185,"62":0.46305,"63":0.44926,"64":0.45825,"65":0.5083566667,"66":0.51753,"67":0.5749233333,"68":0.63139,"69":0.695165,"70":0.72372,"71":0.776485,"72":0.81294,"73":0.8595785714,"74":0.86886,"75":0.902426,"76":0.9147775,"77":0.9045125,"78":0.884404,"79":0.8297333333,"80":0.8209328571,"81":0.7957233333,"82":0.738732,"83":0.70009,"84":0.6890433333,"85":0.67773,"86":0.64621,"87":0.613804,"88":0.5403266667,"89":0.4990371429,"90":0.4692085714,"91":0.43566,"92":0.41087,"93":0.344715,"94":0.304026,"95":0.2863281818}} \ No newline at end of file diff --git a/test/pycodes/create_net.py b/test/pycodes/create_net.py deleted file mode 100644 index 2edaab5..0000000 --- a/test/pycodes/create_net.py +++ /dev/null @@ -1,124 +0,0 @@ -import os -import json -import pathlib -import tempfile -import pandas as pd -import numpy as np -import simbench as sb -import pandapower as pp -import pandapower.networks as pn -from pandapower.converter.powermodels.to_pm import convert_pp_to_pm, init_ne_line, dump_pm_json, convert_to_pm_structure -from pandapower.auxiliary import _add_ppc_options, _add_opf_options -from pandapower.opf.pm_storage import add_storage_opf_settings - -jul_path = pathlib.PurePath(pathlib.Path.home(), ".julia") - -if not os.path.exists(jul_path): - raise KeyError("julia failed. Check julia install!") - -dev_path = pathlib.PurePath(jul_path, "dev", "PandaModels") - -if os.path.exists(dev_path): - Warning("PandaModels is in development mode.") - json_path = pathlib.PurePath(dev_path, "test", "data") -else: - json_path = tempfile.gettempdir() - -types = ["pm", "powerflow", "opf", "custom"] - -net = {type: pp.create_empty_network() - for type in types} - -net["tnep"] = pn.create_cigre_network_mv() -net["ots"] = pn.case5() - -min_vm_pu = 0.95 -max_vm_pu = 1.05 - -for type in types: - #create buses - bus1 = pp.create_bus(net[type], vn_kv=220., geodata=(5,9), min_vm_pu=min_vm_pu, max_vm_pu=max_vm_pu) - bus2 = pp.create_bus(net[type], vn_kv=110., geodata=(6,10), min_vm_pu=min_vm_pu, max_vm_pu=max_vm_pu) - bus3 = pp.create_bus(net[type], vn_kv=110., geodata=(10,9), min_vm_pu=min_vm_pu, max_vm_pu=max_vm_pu) - bus4 = pp.create_bus(net[type], vn_kv=110., geodata=(8,8), min_vm_pu=min_vm_pu, max_vm_pu=max_vm_pu) - bus5 = pp.create_bus(net[type], vn_kv=110., geodata=(6,8), min_vm_pu=min_vm_pu, max_vm_pu=max_vm_pu) - - #create 220/110/110 kV 3W-transformer - pp.create_transformer3w_from_parameters(net[type], bus1, bus2, bus5, vn_hv_kv=220, vn_mv_kv=110, - vn_lv_kv=110, vk_hv_percent=10., vk_mv_percent=10., - vk_lv_percent=10., vkr_hv_percent=0.5, - vkr_mv_percent=0.5, vkr_lv_percent=0.5, pfe_kw=10, - i0_percent=0.1, shift_mv_degree=0, shift_lv_degree=0, - sn_hv_mva=100, sn_mv_mva=50, sn_lv_mva=50) - - #create 110 kV lines - l1 = pp.create_line(net[type], bus2, bus3, length_km=70., std_type='149-AL1/24-ST1A 110.0') - l2 = pp.create_line(net[type], bus3, bus4, length_km=50., std_type='149-AL1/24-ST1A 110.0') - l3 = pp.create_line(net[type], bus4, bus2, length_km=40., std_type='149-AL1/24-ST1A 110.0') - l4 = pp.create_line(net[type], bus4, bus5, length_km=30., std_type='149-AL1/24-ST1A 110.0') - - #create loads - pp.create_load(net[type], bus2, p_mw=60) - pp.create_load(net[type], bus3, p_mw=70) - pp.create_load(net[type], bus4, p_mw=10) - - #create generators - g1 = pp.create_gen(net[type], bus1, p_mw=40, min_p_mw=0, max_p_mw=200, vm_pu=1.01, slack=True) - pp.create_poly_cost(net[type], g1, 'gen', cp1_eur_per_mw=1) - - g2 = pp.create_gen(net[type], bus3, p_mw=40, min_p_mw=0, max_p_mw=200, vm_pu=1.01) - pp.create_poly_cost(net[type], g2, 'gen', cp1_eur_per_mw=3) - - g3 = pp.create_gen(net[type], bus4, p_mw=50, min_p_mw=0, max_p_mw=200, vm_pu=1.01) - pp.create_poly_cost(net[type], g3, 'gen', cp1_eur_per_mw=3) - pp.runpp(net[type]) - - -net["tnep"]["bus"].loc[:, "min_vm_pu"] = min_vm_pu -net["tnep"]["bus"].loc[:, "max_vm_pu"] = max_vm_pu -net["tnep"]["line"].loc[:, "max_loading_percent"] = 80. -net["tnep"]["line"] = pd.concat([net["tnep"]["line"]] * 2, ignore_index=True) -net["tnep"]["line"].loc[max(net["tnep"]["line"].index) + 1:, "in_service"] = False - -new_lines = net["tnep"]["line"].loc[max(net["tnep"]["line"].index) + 1:].index -init_ne_line(net["tnep"], new_lines, construction_costs=np.ones(len(new_lines))) - -pp.runpp(net["tnep"]) - -test_pm_json = os.path.join(json_path, "test_pm.json") -test_powerflow_opf_json = os.path.join(json_path, "test_pf.json") -test_powermodels_json = os.path.join(json_path, "test_opf.json") -test_custom_json = os.path.join(json_path, "test_custom.json") -test_ots_json = os.path.join(json_path, "test_ots.json") -test_tnep_json = os.path.join(json_path, "test_tnep.json") - - -test_pm = convert_pp_to_pm(net["pm"], pm_file_path=test_pm_json, correct_pm_network_data=True, calculate_voltage_angles=True, ac=False, - trafo_model="t", delta=1e-8, trafo3w_losses="hv", check_connectivity=True, - pp_to_pm_callback=None, pm_model="DCPPowerModel", pm_solver="ipopt", - pm_mip_solver="cbc", pm_nl_solver="ipopt") - -test_powerflow = convert_pp_to_pm(net["powerflow"], pm_file_path=test_powerflow_opf_json, correct_pm_network_data=True, calculate_voltage_angles=True, ac=True, - trafo_model="t", delta=1e-8, trafo3w_losses="hv", check_connectivity=True, - pp_to_pm_callback=None, pm_model="ACPPowerModel", pm_solver="ipopt", - pm_mip_solver="cbc", pm_nl_solver="ipopt") - -test_powermodels = convert_pp_to_pm(net["opf"], pm_file_path=test_powermodels_json, correct_pm_network_data=True, calculate_voltage_angles=True, ac=True, - trafo_model="t", delta=1e-8, trafo3w_losses="hv", check_connectivity=True, - pp_to_pm_callback=None, pm_model="ACPPowerModel", pm_solver="ipopt", - pm_mip_solver="cbc", pm_nl_solver="ipopt") - -test_custom = convert_pp_to_pm(net["custom"], pm_file_path=test_custom_json, correct_pm_network_data=True, calculate_voltage_angles=True, ac=False, - trafo_model="t", delta=1e-8, trafo3w_losses="hv", check_connectivity=True, - pp_to_pm_callback=None, pm_model="ACPPowerModel", pm_solver="ipopt", - pm_mip_solver="cbc", pm_nl_solver="ipopt") - -test_ots = convert_pp_to_pm(net["ots"], pm_file_path=test_ots_json, correct_pm_network_data=True, calculate_voltage_angles=True, ac=True, - trafo_model="t", delta=1e-8, trafo3w_losses="hv", check_connectivity=True, - pp_to_pm_callback=None, pm_model="DCPPowerModel", pm_solver="juniper", - pm_mip_solver="cbc", pm_nl_solver="ipopt") - -test_tnep = convert_pp_to_pm(net["tnep"], pm_file_path=test_tnep_json, correct_pm_network_data=True, calculate_voltage_angles=True, ac=True, - trafo_model="t", delta=1e-8, trafo3w_losses="hv", check_connectivity=True, - pp_to_pm_callback=None, pm_model="ACPPowerModel", pm_solver="juniper", - pm_mip_solver="cbc", pm_nl_solver="ipopt") diff --git a/test/pycodes/storage_test_case.py b/test/pycodes/storage_test_case.py deleted file mode 100644 index 9acd081..0000000 --- a/test/pycodes/storage_test_case.py +++ /dev/null @@ -1,70 +0,0 @@ -from pandapower.converter.powermodels.to_pm import convert_pp_to_pm -from timeseries.run_profile_cython.run_cython import profiles_from_ds -import json -import os -import tempfile -import matplotlib.pyplot as plt -import numpy as np -import pandas as pd -import pandapower as pp -import pandapower.networks as nw -from pandapower.control.controller.const_control import ConstControl -from pandapower.timeseries.data_sources.frame_data import DFData - - - -net = nw.create_cigre_network_mv("pv_wind") -# set some limits -min_vm_pu = 0.95 -max_vm_pu = 1.05 - -net["bus"].loc[:, "min_vm_pu"] = min_vm_pu -net["bus"].loc[:, "max_vm_pu"] = max_vm_pu - -net["line"].loc[:, "max_loading_percent"] = 100. - -# close all switches -net.switch.loc[:, "closed"] = True -# add storage to bus 10 -pp.create_storage(net, 10, p_mw=0.5, max_e_mwh=.2, soc_percent=0., q_mvar=0., controllable=True) - -# net.load["controllable"] = False -# net.sgen["controllable"] = False - -net["load"].loc[:, "type"] = "residential" -net.sgen.loc[:, "type"] = "pv" -net.sgen.loc[8, "type"] = "wind" - -input_file = r"C:\Users\x230\.julia\dev\PandaModels\test\pycodes\cigre_timeseries_15min.json" -time_series = pd.read_json(input_file) - - -time_series.sort_index(inplace=True) - -load_df = pd.DataFrame(time_series["residential"]).dot([net["load"]["p_mw"].values]) - -ConstControl(net, "load", "p_mw", element_index=net.load.index[net.load.type == "residential"], - profile_name=net.load.index[net.load.type == "residential"], data_source=DFData(load_df)) - -pv_df = pd.DataFrame(time_series["pv"]).dot([net["sgen"].loc[:7, "p_mw"].values]) - -ConstControl(net, "sgen", "p_mw", element_index=net.sgen.index[net.sgen.type == "pv"], - profile_name=net.sgen.index[net.sgen.type == "pv"], data_source=DFData(pv_df)) - -wind_df = pd.DataFrame(time_series["wind"]).dot([net["sgen"].loc[8, "p_mw"]]) - -ConstControl(net, "sgen", "p_mw", element_index=net.sgen.index[net.sgen.type == "wind"], - profile_name=net.sgen.index[net.sgen.type == "wind"], data_source=DFData(time_series["wind"])) - - -storage_results = pp.runpm_storage_opf(net, n_timesteps=96, time_elapsed=0.25) - -# time_series.sort_values(by="timestep", inplace=True) -# time_series["residential"].plot() -# time_series["pv"].plot() -# time_series["wind"].plot() -# plt.xlabel("time step") -# plt.xlabel("relative value") -# plt.legend() -# plt.grid() -# plt.show() \ No newline at end of file