bjoern/pypsa-eur
1
0
Fork 0
Code Issues Pull Requests Actions 1 Packages Projects Releases Wiki Activity

adjust s_max_pu clipping

Browse Source
This commit is contained in:
Philipp Glaum 2022-02-21 09:18:20 +01:00
parent 86942b0afc
commit 7e692f01e3
1 changed files with 7 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
scripts/add_electricity.py
View File

@ -526,13 +526,13 @@ def estimate_renewable_capacities(n, tech_map):
.where(lambda s: s>0.1, 0.)) # only capacities above 100kW
n.generators.loc[tech_i, 'p_nom_min'] = n.generators.loc[tech_i, 'p_nom']
def attach_line_rating(n, fn, s_max_pu_factor):
def attach_line_rating(n, fn, s_max_pu_factor,dlr_factor):
s_max = xr.open_dataarray(fn).to_pandas().transpose()
n.lines_t.s_max_pu = s_max / n.lines.s_nom[s_max.columns] #only considers overhead lines
# account for maximal voltage angles of maximally 30 degree.
x = n.lines.type.map(n.line_types["x_per_length"])*n.lines.length/(n.lines.v_nom**2)
s_max_pu_cap = (np.pi / (6 * x * n.lines.s_nom)).clip(lower=1) # need to clip here as cap values are below 1 -> would mean the line cannot be operated at actual given pessimistic ampacity
n.lines_t.s_max_pu = n.lines_t.s_max_pu.clip(upper=s_max_pu_cap, lower=1, axis=1)*s_max_pu_factor
x_pu = n.lines.type.map(n.line_types["x_per_length"])*n.lines.length/(n.lines.v_nom**2)
s_max_pu_cap = (np.pi / (6 * x_pu * n.lines.s_nom)).clip(lower=1) # need to clip here as cap values are below 1 -> would mean the line cannot be operated at actual given pessimistic ampacity
n.lines_t.s_max_pu = n.lines_t.s_max_pu.clip(upper=s_max_pu_cap, lower=1, axis=1)*s_max_pu_factor*dlr_factor
def add_nice_carrier_names(n, config):
carrier_i = n.carriers.index
@ -584,8 +584,9 @@ if __name__ == "__main__":
update_p_nom_max(n)
if snakemake.config["lines"]["line_rating"]:
s_max_pu_factor = snakemake.config["lines"]["s_max_pu"]
attach_line_rating(n, snakemake.input.line_rating, s_max_pu_factor)
s_max_pu_factor = snakemake.config["lines"]["s_max_pu"] #factor mainly used for N-1 security
dlr_factor=snakemake.config["lines"]["dlr_factor"] #factor due to overestimation of the wind speed in hourly averaged wind data
attach_line_rating(n, snakemake.input.line_rating, s_max_pu_factor, dlr_factor)
add_nice_carrier_names(n, snakemake.config)