Hi all,
I just noticed that the sector_opts wildcard for changing the CO2 sequestration limit seems to not be working.
I have found the reason, which is that the e_nom_opt of the co2 sequestration store is set to np.inf(), so multiplying it by any factor obviously does nothing.
n.madd("Store",
spatial.co2.nodes,
e_nom_extendable=True,
e_nom_max=np.inf,
capital_cost=options['co2_sequestration_cost'],
carrier="co2 stored",
bus=spatial.co2.nodes
)
The actual limit is applied in the script solve_network.py, in the function add_co2_sequestration_limit(), but it expects the wildcard to be formatted as seq2000000 or seq2e8 (for example) instead of co2 stored+e2
However, one line in solve_network is incorrect. It is now:
opts = snakemake.wildcards.opts.split('-')
But should be
opts = snakemake.wildcards.sector_opts.split('-')
I'm just writing this message 1) in case anyone else has this problem, and 2) to make sure the developers are aware that the these issues exist. I can write a pull request to make the change in solve_network.py, but I'm not sure if there are different plans to handle sequestration.
Also, I'm wondering if there's a specific justification for 200Mt as the default sequestration potential value. I was copied on an email from Dr. Brown a few months ago that said it was in order to compensate for process emissions, but it seems like 150Mt is the usual amount.
Best,
Paul
Thanks,
Paul