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

add params for solve_network

Browse Source
This commit is contained in:
virio-andreyana 2023-05-27 12:22:53 +02:00
parent 28d39680c6
commit 1c73bb0342
4 changed files with 37 additions and 19 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
rules/solve_electricity.smk
View File

@ -6,6 +6,12 @@
rule solve_network: rule solve_network:
params: params:
solving=config["solving"], solving=config["solving"],
config_parts={
"foresight":config["foresight"],
"planning_horizons":config["scenario"]["planning_horizons"],
"co2_sequestration_potential":config["sector"].get("co2_sequestration_potential", 200),
"H2_retrofit_capacity_per_CH4":config["sector"]["H2_retrofit_capacity_per_CH4"]
},
input: input:
network=RESOURCES + "networks/elec_s{simpl}_{clusters}_ec_l{ll}_{opts}.nc", network=RESOURCES + "networks/elec_s{simpl}_{clusters}_ec_l{ll}_{opts}.nc",
output: output:

6
rules/solve_myopic.smk
View File

@ -84,6 +84,12 @@ ruleorder: add_existing_baseyear > add_brownfield
rule solve_sector_network_myopic: rule solve_sector_network_myopic:
params: params:
solving=config["solving"], solving=config["solving"],
config_parts={
"foresight":config["foresight"],
"planning_horizons":config["scenario"]["planning_horizons"],
"co2_sequestration_potential":config["sector"].get("co2_sequestration_potential", 200),
"H2_retrofit_capacity_per_CH4":config["sector"]["H2_retrofit_capacity_per_CH4"]
},
input: input:
overrides="data/override_component_attrs", overrides="data/override_component_attrs",
network=RESULTS network=RESULTS

6
rules/solve_overnight.smk
View File

@ -6,6 +6,12 @@
rule solve_sector_network: rule solve_sector_network:
params: params:
solving=config["solving"], solving=config["solving"],
config_parts={
"foresight":config["foresight"],
"planning_horizons":config["scenario"]["planning_horizons"],
"co2_sequestration_potential":config["sector"].get("co2_sequestration_potential", 200),
"H2_retrofit_capacity_per_CH4":config["sector"]["H2_retrofit_capacity_per_CH4"]
},
input: input:
overrides="data/override_component_attrs", overrides="data/override_component_attrs",
network=RESULTS network=RESULTS

38
scripts/solve_network.py
View File

@ -44,14 +44,14 @@ pypsa.pf.logger.setLevel(logging.WARNING)
from pypsa.descriptors import get_switchable_as_dense as get_as_dense from pypsa.descriptors import get_switchable_as_dense as get_as_dense
def add_land_use_constraint(n, config): def add_land_use_constraint(n, param, config):
if "m" in snakemake.wildcards.clusters: if "m" in snakemake.wildcards.clusters:
_add_land_use_constraint_m(n, config) _add_land_use_constraint_m(n, param, config)
else: else:
_add_land_use_constraint(n, config) _add_land_use_constraint(n, param)
def _add_land_use_constraint(n, config): def _add_land_use_constraint(n, param):
# warning: this will miss existing offwind which is not classed AC-DC and has carrier 'offwind' # warning: this will miss existing offwind which is not classed AC-DC and has carrier 'offwind'
for carrier in ["solar", "onwind", "offwind-ac", "offwind-dc"]: for carrier in ["solar", "onwind", "offwind-ac", "offwind-dc"]:
@ -80,10 +80,10 @@ def _add_land_use_constraint(n, config):
n.generators.p_nom_max.clip(lower=0, inplace=True) n.generators.p_nom_max.clip(lower=0, inplace=True)
def _add_land_use_constraint_m(n, config): def _add_land_use_constraint_m(n, param, config):
# if generators clustering is lower than network clustering, land_use accounting is at generators clusters # if generators clustering is lower than network clustering, land_use accounting is at generators clusters
planning_horizons = config["scenario"]["planning_horizons"] planning_horizons = param["planning_horizons"]
grouping_years = config["existing_capacities"]["grouping_years"] grouping_years = config["existing_capacities"]["grouping_years"]
current_horizon = snakemake.wildcards.planning_horizons current_horizon = snakemake.wildcards.planning_horizons
@ -141,7 +141,7 @@ def add_co2_sequestration_limit(n, limit=200):
) )
def prepare_network(n, solve_opts=None, config=None): def prepare_network(n, solve_opts=None, config=None, param=None):
if "clip_p_max_pu" in solve_opts: if "clip_p_max_pu" in solve_opts:
for df in ( for df in (
n.generators_t.p_max_pu, n.generators_t.p_max_pu,
@ -191,11 +191,11 @@ def prepare_network(n, solve_opts=None, config=None):
n.set_snapshots(n.snapshots[:nhours]) n.set_snapshots(n.snapshots[:nhours])
n.snapshot_weightings[:] = 8760.0 / nhours n.snapshot_weightings[:] = 8760.0 / nhours
if config["foresight"] == "myopic": if param["foresight"] == "myopic":
add_land_use_constraint(n, config) add_land_use_constraint(n, param, config)
if n.stores.carrier.eq("co2 stored").any(): if n.stores.carrier.eq("co2 stored").any():
limit = config["sector"].get("co2_sequestration_potential", 200) limit = param["co2_sequestration_potential"]
add_co2_sequestration_limit(n, limit=limit) add_co2_sequestration_limit(n, limit=limit)
return n return n
@ -222,7 +222,7 @@ def add_CCL_constraints(n, config):
agg_p_nom_limits: data/agg_p_nom_minmax.csv agg_p_nom_limits: data/agg_p_nom_minmax.csv
""" """
agg_p_nom_minmax = pd.read_csv( agg_p_nom_minmax = pd.read_csv(
config["electricity"]["agg_p_nom_limits"], index_col=[0, 1] config['electricity']["agg_p_nom_limits"], index_col=[0, 1]
) )
logger.info("Adding generation capacity constraints per carrier and country") logger.info("Adding generation capacity constraints per carrier and country")
p_nom = n.model["Generator-p_nom"] p_nom = n.model["Generator-p_nom"]
@ -370,7 +370,7 @@ def add_SAFE_constraints(n, config):
Which sets a reserve margin of 10% above the peak demand. Which sets a reserve margin of 10% above the peak demand.
""" """
peakdemand = n.loads_t.p_set.sum(axis=1).max() peakdemand = n.loads_t.p_set.sum(axis=1).max()
margin = 1.0 + config["electricity"]["SAFE_reservemargin"] margin = 1.0 + config['electricity']["SAFE_reservemargin"]
reserve_margin = peakdemand * margin reserve_margin = peakdemand * margin
# TODO: do not take this from the plotting config! # TODO: do not take this from the plotting config!
conv_techs = config["plotting"]["conv_techs"] conv_techs = config["plotting"]["conv_techs"]
@ -590,13 +590,13 @@ def extra_functionality(n, snapshots):
add_pipe_retrofit_constraint(n) add_pipe_retrofit_constraint(n)
def solve_network(n, config, opts="", **kwargs): def solve_network(n, config, solving_param, opts="", **kwargs):
set_of_options = config["solving"]["solver"]["options"] set_of_options = solving_param["solver"]["options"]
solver_options = ( solver_options = (
config["solving"]["solver_options"][set_of_options] if set_of_options else {} solving_param["solver_options"][set_of_options] if set_of_options else {}
) )
solver_name = config["solving"]["solver"]["name"] solver_name = solving_param["solver"]["name"]
cf_solving = config["solving"]["options"] cf_solving = solving_param["options"]
track_iterations = cf_solving.get("track_iterations", False) track_iterations = cf_solving.get("track_iterations", False)
min_iterations = cf_solving.get("min_iterations", 4) min_iterations = cf_solving.get("min_iterations", 4)
max_iterations = cf_solving.get("max_iterations", 6) max_iterations = cf_solving.get("max_iterations", 6)
@ -675,10 +675,10 @@ if __name__ == "__main__":
else: else:
n = pypsa.Network(snakemake.input.network) n = pypsa.Network(snakemake.input.network)
n = prepare_network(n, solve_opts, config=snakemake.config) n = prepare_network(n, solve_opts, config=snakemake.config, param=snakemake.params["config_parts"])
n = solve_network( n = solve_network(
n, config=snakemake.config, opts=opts, log_fn=snakemake.log.solver n, config=snakemake.config, solving_param=snakemake.params["solving"], opts=opts, log_fn=snakemake.log.solver
) )
n.meta = dict(snakemake.config, **dict(wildcards=dict(snakemake.wildcards))) n.meta = dict(snakemake.config, **dict(wildcards=dict(snakemake.wildcards)))