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

add other co2 potentials

Browse Source
This commit is contained in:
lisazeyen 2024-04-10 12:06:53 +02:00
parent e1502c95dd
commit 9d9ab98279
2 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
config/config.default.yaml
View File

@ -515,7 +515,7 @@ sector:
regional_coal_demand: false regional_coal_demand: false
regional_co2_sequestration_potential: regional_co2_sequestration_potential:
enable: false enable: false
attribute: 'conservative estimate Mt' attribute: ['conservative estimate Mt', 'conservative estimate GAS Mt', 'conservative estimate OIL Mt', 'conservative estimate aquifer Mt']
include_onshore: false include_onshore: false
min_size: 3 min_size: 3
max_size: 25 max_size: 25

6
scripts/build_sequestration_potentials.py
View File

@ -23,13 +23,13 @@ def area(gdf):
def allocate_sequestration_potential( def allocate_sequestration_potential(
gdf, regions, attr="conservative estimate Mt", threshold=3 gdf, regions, attr="conservative estimate Mt", threshold=3
): ):
gdf = gdf.loc[gdf[attr] > threshold, [attr, "geometry"]] gdf = gdf.loc[gdf[attr].sum(axis=1) > threshold, attr + ["geometry"]]
gdf["area_sqkm"] = area(gdf) gdf["area_sqkm"] = area(gdf)
overlay = gpd.overlay(regions, gdf, keep_geom_type=True) overlay = gpd.overlay(regions, gdf, keep_geom_type=True)
overlay["share"] = area(overlay) / overlay["area_sqkm"] overlay["share"] = area(overlay) / overlay["area_sqkm"]
adjust_cols = overlay.columns.difference({"name", "area_sqkm", "geometry", "share"}) adjust_cols = overlay.columns.difference({"name", "area_sqkm", "geometry", "share"})
overlay[adjust_cols] = overlay[adjust_cols].multiply(overlay["share"], axis=0) overlay[adjust_cols] = overlay[adjust_cols].multiply(overlay["share"], axis=0)
return overlay.dissolve("name", aggfunc="sum")[attr] return overlay.dissolve("name", aggfunc="sum")[attr].sum(axis=1)
if __name__ == "__main__": if __name__ == "__main__":
@ -37,7 +37,7 @@ if __name__ == "__main__":
from _helpers import mock_snakemake from _helpers import mock_snakemake
snakemake = mock_snakemake( snakemake = mock_snakemake(
"build_sequestration_potentials", simpl="", clusters="181" "build_sequestration_potentials", simpl="", clusters="128"
) )
set_scenario_config(snakemake) set_scenario_config(snakemake)