diff --git a/config/config.default.yaml b/config/config.default.yaml index 63cdcabff..246a15f47 100644 --- a/config/config.default.yaml +++ b/config/config.default.yaml @@ -605,7 +605,6 @@ sector: hydrogen_turbine: false SMR: true SMR_cc: true - regional_methanol_demand: false regional_oil_demand: false regional_coal_demand: false regional_co2_sequestration_potential: @@ -638,6 +637,19 @@ sector: # - onshore # more than 50 km from sea - nearshore # within 50 km of sea # - offshore + methanol: + regional_methanol_demand: false + methanol_reforming: false + methanol_reforming_cc: false + methanol_to_kerosene: false + methanol_to_olefins: false + methanol_to_power: + ccgt: false + ccgt_cc: false + ocgt: false + allam: false + biomass_to_methanol: false + biomass_to_methanol_cc: false ammonia: false min_part_load_fischer_tropsch: 0.5 min_part_load_methanolisation: 0.3 @@ -803,7 +815,7 @@ industry: # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#costs costs: year: 2030 - version: v0.9.1 + version: v0.9.2 social_discountrate: 0.02 fill_values: FOM: 0 @@ -1254,8 +1266,19 @@ plotting: liquid: '#25c49a' kerosene for aviation: '#a1ffe6' naphtha for industry: '#57ebc4' - methanolisation: '#83d6d5' - methanol: '#468c8b' + methanol-to-kerosene: '#C98468' + methanol-to-olefins/aromatics: '#FFA07A' + Methanol steam reforming: '#FFBF00' + Methanol steam reforming CC: '#A2EA8A' + methanolisation: '#00FFBF' + biomass-to-methanol: '#EAD28A' + biomass-to-methanol CC: '#EADBAD' + allam methanol: '#B98F76' + CCGT methanol: '#B98F76' + CCGT methanol CC: '#B98F76' + OCGT methanol: '#B98F76' + methanol: '#FF7B00' + methanol transport: '#FF7B00' shipping methanol: '#468c8b' industry methanol: '#468c8b' # co2 diff --git a/doc/configtables/sector.csv b/doc/configtables/sector.csv index 7de2a173d..000927ced 100644 --- a/doc/configtables/sector.csv +++ b/doc/configtables/sector.csv @@ -1,182 +1,193 @@ -,Unit,Values,Description -transport,--,"{true, false}",Flag to include transport sector. -heating,--,"{true, false}",Flag to include heating sector. -biomass,--,"{true, false}",Flag to include biomass sector. -industry,--,"{true, false}",Flag to include industry sector. -agriculture,--,"{true, false}",Flag to include agriculture sector. -fossil_fuels,--,"{true, false}","Flag to include imports of fossil fuels ( [""coal"", ""gas"", ""oil"", ""lignite""])" -district_heating,--,,`prepare_sector_network.py `_ --- potential,--,float,maximum fraction of urban demand which can be supplied by district heating --- progress,--,Dictionary with planning horizons as keys., Increase of today's district heating demand to potential maximum district heating share. Progress = 0 means today's district heating share. Progress = 1 means maximum fraction of urban demand is supplied by district heating --- district_heating_loss,--,float,Share increase in district heat demand in urban central due to heat losses --- supply_temperature_approximation,,, --- -- max_forward_temperature,°C,Dictionary with country codes as keys. One key must be 'default'., Max. forward temperature in district heating (if ambient temperature lower-or-equal `lower_threshold_ambient_temperature`) --- -- min_forward_temperature,°C,Dictionary with country codes as keys. One key must be 'default'., Min. forward temperature in district heating (if ambient temperature higher-or-equal `upper_threshold_ambient_temperature`) --- -- return_temperature,°C,Dictionary with country codes as keys. One key must be 'default'.,Return temperature in district heating. Must be lower than forward temperature --- -- lower_threshold_ambient_temperature,°C,float, Assume `max_forward_temperature` if ambient temperature is below this threshold --- -- upper_threshold_ambient_temperature,°C,float, Assume `min_forward_temperature` if ambient temperature is above this threshold --- -- rolling_window_ambient_temperature, h, int, Rolling window size for averaging ambient temperature when approximating supply temperature --- heat_source_cooling,K,float,Cooling of heat source for heat pumps --- heat_pump_cop_approximation,,, --- -- refrigerant,--,"{ammonia, isobutane}",Heat pump refrigerant assumed for COP approximation --- -- heat_exchanger_pinch_point_temperature_difference,K,float,Heat pump pinch point temperature difference in heat exchangers assumed for approximation. --- -- isentropic_compressor_efficiency,--,float,Isentropic efficiency of heat pump compressor assumed for approximation. Must be between 0 and 1. --- -- heat_loss,--,float,Heat pump heat loss assumed for approximation. Must be between 0 and 1. --- heat_pump_sources,--,, --- -- urban central,--,List of heat sources for heat pumps in urban central heating, --- -- urban decentral,--,List of heat sources for heat pumps in urban decentral heating, --- -- rural,--,List of heat sources for heat pumps in rural heating, -cluster_heat_buses,--,"{true, false}",Cluster residential and service heat buses in `prepare_sector_network.py `_ to one to save memory. -,,, -bev_dsm_restriction _value,--,float,Adds a lower state of charge (SOC) limit for battery electric vehicles (BEV) to manage its own energy demand (DSM). Located in `build_transport_demand.py `_. Set to 0 for no restriction on BEV DSM -bev_dsm_restriction _time,--,float,Time at which SOC of BEV has to be dsm_restriction_value -transport_heating _deadband_upper,°C,float,"The maximum temperature in the vehicle. At higher temperatures, the energy required for cooling in the vehicle increases." -transport_heating _deadband_lower,°C,float,"The minimum temperature in the vehicle. At lower temperatures, the energy required for heating in the vehicle increases." -,,, -ICE_lower_degree_factor,--,float,Share increase in energy demand in internal combustion engine (ICE) for each degree difference between the cold environment and the minimum temperature. -ICE_upper_degree_factor,--,float,Share increase in energy demand in internal combustion engine (ICE) for each degree difference between the hot environment and the maximum temperature. -EV_lower_degree_factor,--,float,Share increase in energy demand in electric vehicles (EV) for each degree difference between the cold environment and the minimum temperature. -EV_upper_degree_factor,--,float,Share increase in energy demand in electric vehicles (EV) for each degree difference between the hot environment and the maximum temperature. -bev_dsm,--,"{true, false}",Add the option for battery electric vehicles (BEV) to participate in demand-side management (DSM) -,,, -bev_availability,--,float,The share for battery electric vehicles (BEV) that are able to do demand side management (DSM) -bev_energy,--,float,The average size of battery electric vehicles (BEV) in MWh -bev_charge_efficiency,--,float,Battery electric vehicles (BEV) charge and discharge efficiency -bev_charge_rate,MWh,float,The power consumption for one electric vehicle (EV) in MWh. Value derived from 3-phase charger with 11 kW. -bev_avail_max,--,float,The maximum share plugged-in availability for passenger electric vehicles. -bev_avail_mean,--,float,The average share plugged-in availability for passenger electric vehicles. -v2g,--,"{true, false}",Allows feed-in to grid from EV battery -land_transport_fuel_cell _share,--,Dictionary with planning horizons as keys.,The share of vehicles that uses fuel cells in a given year -land_transport_electric _share,--,Dictionary with planning horizons as keys.,The share of vehicles that uses electric vehicles (EV) in a given year -land_transport_ice _share,--,Dictionary with planning horizons as keys.,The share of vehicles that uses internal combustion engines (ICE) in a given year. What is not EV or FCEV is oil-fuelled ICE. -transport_electric_efficiency,MWh/100km,float,The conversion efficiencies of electric vehicles in transport -transport_fuel_cell_efficiency,MWh/100km,float,The H2 conversion efficiencies of fuel cells in transport -transport_ice_efficiency,MWh/100km,float,The oil conversion efficiencies of internal combustion engine (ICE) in transport -agriculture_machinery _electric_share,--,float,The share for agricultural machinery that uses electricity -agriculture_machinery _oil_share,--,float,The share for agricultural machinery that uses oil -agriculture_machinery _fuel_efficiency,--,float,The efficiency of electric-powered machinery in the conversion of electricity to meet agricultural needs. -agriculture_machinery _electric_efficiency,--,float,The efficiency of oil-powered machinery in the conversion of oil to meet agricultural needs. -Mwh_MeOH_per_MWh_H2,LHV,float,"The energy amount of the produced methanol per energy amount of hydrogen. From `DECHEMA (2017) `_, page 64." -MWh_MeOH_per_tCO2,LHV,float,"The energy amount of the produced methanol per ton of CO2. From `DECHEMA (2017) `_, page 66." -MWh_MeOH_per_MWh_e,LHV,float,"The energy amount of the produced methanol per energy amount of electricity. From `DECHEMA (2017) `_, page 64." -shipping_hydrogen _liquefaction,--,"{true, false}",Whether to include liquefaction costs for hydrogen demand in shipping. -,,, -shipping_hydrogen_share,--,Dictionary with planning horizons as keys.,The share of ships powered by hydrogen in a given year -shipping_methanol_share,--,Dictionary with planning horizons as keys.,The share of ships powered by methanol in a given year -shipping_oil_share,--,Dictionary with planning horizons as keys.,The share of ships powered by oil in a given year -shipping_methanol _efficiency,--,float,The efficiency of methanol-powered ships in the conversion of methanol to meet shipping needs (propulsion). The efficiency increase from oil can be 10-15% higher according to the `IEA `_ -,,, -shipping_oil_efficiency,--,float,The efficiency of oil-powered ships in the conversion of oil to meet shipping needs (propulsion). Base value derived from 2011 -aviation_demand_factor,--,float,The proportion of demand for aviation compared to today's consumption -HVC_demand_factor,--,float,The proportion of demand for high-value chemicals compared to today's consumption -,,, -time_dep_hp_cop,--,"{true, false}",Consider the time dependent coefficient of performance (COP) of the heat pump -heat_pump_sink_T,°C,float,The temperature heat sink used in heat pumps based on DTU / large area radiators. The value is conservatively high to cover hot water and space heating in poorly-insulated buildings -reduce_space_heat _exogenously,--,"{true, false}",Influence on space heating demand by a certain factor (applied before losses in district heating). -reduce_space_heat _exogenously_factor,--,Dictionary with planning horizons as keys.,"A positive factor can mean renovation or demolition of a building. If the factor is negative, it can mean an increase in floor area, increased thermal comfort, population growth. The default factors are determined by the `Eurocalc Homes and buildings decarbonization scenario `_" -retrofitting,,, --- retro_endogen,--,"{true, false}",Add retrofitting as an endogenous system which co-optimise space heat savings. --- cost_factor,--,float,Weight costs for building renovation --- interest_rate,--,float,The interest rate for investment in building components --- annualise_cost,--,"{true, false}",Annualise the investment costs of retrofitting --- tax_weighting,--,"{true, false}",Weight the costs of retrofitting depending on taxes in countries --- construction_index,--,"{true, false}",Weight the costs of retrofitting depending on labour/material costs per country -tes,--,"{true, false}",Add option for storing thermal energy in large water pits associated with district heating systems and individual thermal energy storage (TES) -tes_tau,,,The time constant used to calculate the decay of thermal energy in thermal energy storage (TES): 1- :math:`e^{-1/24τ}`. --- decentral,days,float,The time constant in decentralized thermal energy storage (TES) --- central,days,float,The time constant in centralized thermal energy storage (TES) -boilers,--,"{true, false}",Add option for transforming gas into heat using gas boilers -resistive_heaters,--,"{true, false}",Add option for transforming electricity into heat using resistive heaters (independently from gas boilers) -oil_boilers,--,"{true, false}",Add option for transforming oil into heat using boilers -biomass_boiler,--,"{true, false}",Add option for transforming biomass into heat using boilers -overdimension_heat_generators,,,Add option for overdimensioning heating systems by a certain factor. This allows them to cover heat demand peaks e.g. 10% higher than those in the data with a setting of 1.1. --- decentral,--,float,The factor for overdimensioning (increasing CAPEX) decentral heating systems --- central,--,float,The factor for overdimensioning (increasing CAPEX) central heating systems -chp,--,"{true, false}",Add option for using Combined Heat and Power (CHP) -micro_chp,--,"{true, false}",Add option for using Combined Heat and Power (CHP) for decentral areas. -solar_thermal,--,"{true, false}",Add option for using solar thermal to generate heat. -solar_cf_correction,--,float,The correction factor for the value provided by the solar thermal profile calculations -marginal_cost_storage,currency/MWh ,float,The marginal cost of discharging batteries in distributed grids -methanation,--,"{true, false}",Add option for transforming hydrogen and CO2 into methane using methanation. -coal_cc,--,"{true, false}",Add option for coal CHPs with carbon capture -dac,--,"{true, false}",Add option for Direct Air Capture (DAC) -co2_vent,--,"{true, false}",Add option for vent out CO2 from storages to the atmosphere. -allam_cycle,--,"{true, false}",Add option to include `Allam cycle gas power plants `_ -hydrogen_fuel_cell,--,"{true, false}",Add option to include hydrogen fuel cell for re-electrification. Assuming OCGT technology costs -hydrogen_turbine,--,"{true, false}",Add option to include hydrogen turbine for re-electrification. Assuming OCGT technology costs -SMR,--,"{true, false}",Add option for transforming natural gas into hydrogen and CO2 using Steam Methane Reforming (SMR) -SMR CC,--,"{true, false}",Add option for transforming natural gas into hydrogen and CO2 using Steam Methane Reforming (SMR) and Carbon Capture (CC) -regional_methanol_demand,--,"{true, false}",Spatially resolve methanol demand. Set to true if regional CO2 constraints needed. -regional_oil_demand,--,"{true, false}",Spatially resolve oil demand. Set to true if regional CO2 constraints needed. -regional_co2 _sequestration_potential,,, --- enable,--,"{true, false}",Add option for regionally-resolved geological carbon dioxide sequestration potentials based on `CO2StoP `_. --- attribute,--,string or list,Name (or list of names) of the attribute(s) for the sequestration potential --- include_onshore,--,"{true, false}",Add options for including onshore sequestration potentials --- min_size,Gt ,float,Any sites with lower potential than this value will be excluded --- max_size,Gt ,float,The maximum sequestration potential for any one site. --- years_of_storage,years,float,The years until potential exhausted at optimised annual rate -co2_sequestration_potential,--,Dictionary with planning horizons as keys.,The potential of sequestering CO2 in Europe per year and investment period -co2_sequestration_cost,currency/tCO2,float,The cost of sequestering a ton of CO2 -co2_sequestration_lifetime,years,int,The lifetime of a CO2 sequestration site -co2_spatial,--,"{true, false}","Add option to spatially resolve carrier representing stored carbon dioxide. This allows for more detailed modelling of CCUTS, e.g. regarding the capturing of industrial process emissions, usage as feedstock for electrofuels, transport of carbon dioxide, and geological sequestration sites." -,,, -co2network,--,"{true, false}",Add option for planning a new carbon dioxide transmission network -co2_network_cost_factor,p.u.,float,The cost factor for the capital cost of the carbon dioxide transmission network -,,, -cc_fraction,--,float,The default fraction of CO2 captured with post-combustion capture -hydrogen_underground _storage,--,"{true, false}",Add options for storing hydrogen underground. Storage potential depends regionally. -hydrogen_underground _storage_locations,,"{onshore, nearshore, offshore}","The location where hydrogen underground storage can be located. Onshore, nearshore, offshore means it must be located more than 50 km away from the sea, within 50 km of the sea, or within the sea itself respectively." -,,, -ammonia,--,"{true, false, regional}","Add ammonia as a carrrier. It can be either true (copperplated NH3), false (no NH3 carrier) or ""regional"" (regionalised NH3 without network)" -min_part_load_fischer _tropsch,per unit of p_nom ,float,The minimum unit dispatch (``p_min_pu``) for the Fischer-Tropsch process -min_part_load _methanolisation,per unit of p_nom ,float,The minimum unit dispatch (``p_min_pu``) for the methanolisation process -,,, -use_fischer_tropsch _waste_heat,--,"{true, false}",Add option for using waste heat of Fischer Tropsch in district heating networks -use_fuel_cell_waste_heat,--,"{true, false}",Add option for using waste heat of fuel cells in district heating networks -use_electrolysis_waste _heat,--,"{true, false}",Add option for using waste heat of electrolysis in district heating networks -electricity_transmission _grid,--,"{true, false}",Switch for enabling/disabling the electricity transmission grid. -electricity_distribution _grid,--,"{true, false}",Add a simplified representation of the exchange capacity between transmission and distribution grid level through a link. -electricity_distribution _grid_cost_factor,,,Multiplies the investment cost of the electricity distribution grid -,,, -electricity_grid _connection,--,"{true, false}",Add the cost of electricity grid connection for onshore wind and solar -transmission_efficiency,,,Section to specify transmission losses or compression energy demands of bidirectional links. Splits them into two capacity-linked unidirectional links. --- {carrier},--,str,The carrier of the link. --- -- efficiency_static,p.u.,float,Length-independent transmission efficiency. --- -- efficiency_per_1000km,p.u. per 1000 km,float,Length-dependent transmission efficiency ($\eta^{\text{length}}$) --- -- compression_per_1000km,p.u. per 1000 km,float,Length-dependent electricity demand for compression ($\eta \cdot \text{length}$) implemented as multi-link to local electricity bus. -H2_network,--,"{true, false}",Add option for new hydrogen pipelines -gas_network,--,"{true, false}","Add existing natural gas infrastructure, incl. LNG terminals, production and entry-points. The existing gas network is added with a lossless transport model. A length-weighted `k-edge augmentation algorithm `_ can be run to add new candidate gas pipelines such that all regions of the model can be connected to the gas network. When activated, all the gas demands are regionally disaggregated as well." -H2_retrofit,--,"{true, false}",Add option for retrofiting existing pipelines to transport hydrogen. -H2_retrofit_capacity _per_CH4,--,float,"The ratio for H2 capacity per original CH4 capacity of retrofitted pipelines. The `European Hydrogen Backbone (April, 2020) p.15 `_ 60% of original natural gas capacity could be used in cost-optimal case as H2 capacity." -gas_network_connectivity _upgrade ,--,float,The number of desired edge connectivity (k) in the length-weighted `k-edge augmentation algorithm `_ used for the gas network -gas_distribution_grid,--,"{true, false}",Add a gas distribution grid -gas_distribution_grid _cost_factor,,,Multiplier for the investment cost of the gas distribution grid -,,, -biomass_spatial,--,"{true, false}",Add option for resolving biomass demand regionally -biomass_transport,--,"{true, false}",Add option for transporting solid biomass between nodes -biogas_upgrading_cc,--,"{true, false}",Add option to capture CO2 from biomass upgrading -conventional_generation,,,Add a more detailed description of conventional carriers. Any power generation requires the consumption of fuel from nodes representing that fuel. -biomass_to_liquid,--,"{true, false}",Add option for transforming solid biomass into liquid fuel with the same properties as oil -biosng,--,"{true, false}",Add option for transforming solid biomass into synthesis gas with the same properties as natural gas -bioH2,--,"{true, false}",Add option for transforming solid biomass into hydrogen with carbon capture -municipal_solid_waste,--,"{true, false}",Add option for municipal solid waste -limit_max_growth,,, --- enable,--,"{true, false}",Add option to limit the maximum growth of a carrier --- factor,p.u.,float,The maximum growth factor of a carrier (e.g. 1.3 allows 30% larger than max historic growth) --- max_growth,,, --- -- {carrier},GW,float,The historic maximum growth of a carrier --- max_relative_growth,,, --- -- {carrier},p.u.,float,The historic maximum relative growth of a carrier -,,, -enhanced_geothermal,,, --- enable,--,"{true, false}",Add option to include Enhanced Geothermal Systems --- flexible,--,"{true, false}",Add option for flexible operation (see Ricks et al. 2024) --- max_hours,--,int,The maximum hours the reservoir can be charged under flexible operation --- max_boost,--,float,The maximum boost in power output under flexible operation --- var_cf,--,"{true, false}",Add option for variable capacity factor (see Ricks et al. 2024) --- sustainability_factor,--,float,Share of sourced heat that is replenished by the earth's core (see details in `build_egs_potentials.py `_) -solid_biomass_import,,, --- enable,--,"{true, false}",Add option to include solid biomass imports --- price,currency/MWh,float,Price for importing solid biomass --- max_amount,Twh,float,Maximum solid biomass import potential --- upstream_emissions_factor,p.u.,float,Upstream emissions of solid biomass imports +,Unit,Values,Description +transport,--,"{true, false}",Flag to include transport sector. +heating,--,"{true, false}",Flag to include heating sector. +biomass,--,"{true, false}",Flag to include biomass sector. +industry,--,"{true, false}",Flag to include industry sector. +agriculture,--,"{true, false}",Flag to include agriculture sector. +fossil_fuels,--,"{true, false}","Flag to include imports of fossil fuels ( [""coal"", ""gas"", ""oil"", ""lignite""])" +district_heating,--,,`prepare_sector_network.py `_ +-- potential,--,float,maximum fraction of urban demand which can be supplied by district heating +-- progress,--,Dictionary with planning horizons as keys., Increase of today's district heating demand to potential maximum district heating share. Progress = 0 means today's district heating share. Progress = 1 means maximum fraction of urban demand is supplied by district heating +-- district_heating_loss,--,float,Share increase in district heat demand in urban central due to heat losses +-- supply_temperature_approximation,,, +-- -- max_forward_temperature,°C,Dictionary with country codes as keys. One key must be 'default'., Max. forward temperature in district heating (if ambient temperature lower-or-equal `lower_threshold_ambient_temperature`) +-- -- min_forward_temperature,°C,Dictionary with country codes as keys. One key must be 'default'., Min. forward temperature in district heating (if ambient temperature higher-or-equal `upper_threshold_ambient_temperature`) +-- -- return_temperature,°C,Dictionary with country codes as keys. One key must be 'default'.,Return temperature in district heating. Must be lower than forward temperature +-- -- lower_threshold_ambient_temperature,°C,float, Assume `max_forward_temperature` if ambient temperature is below this threshold +-- -- upper_threshold_ambient_temperature,°C,float, Assume `min_forward_temperature` if ambient temperature is above this threshold +-- -- rolling_window_ambient_temperature, h, int, Rolling window size for averaging ambient temperature when approximating supply temperature +-- heat_source_cooling,K,float,Cooling of heat source for heat pumps +-- heat_pump_cop_approximation,,, +-- -- refrigerant,--,"{ammonia, isobutane}",Heat pump refrigerant assumed for COP approximation +-- -- heat_exchanger_pinch_point_temperature_difference,K,float,Heat pump pinch point temperature difference in heat exchangers assumed for approximation. +-- -- isentropic_compressor_efficiency,--,float,Isentropic efficiency of heat pump compressor assumed for approximation. Must be between 0 and 1. +-- -- heat_loss,--,float,Heat pump heat loss assumed for approximation. Must be between 0 and 1. +-- heat_pump_sources,--,, +-- -- urban central,--,List of heat sources for heat pumps in urban central heating, +-- -- urban decentral,--,List of heat sources for heat pumps in urban decentral heating, +-- -- rural,--,List of heat sources for heat pumps in rural heating, +cluster_heat_buses,--,"{true, false}",Cluster residential and service heat buses in `prepare_sector_network.py `_ to one to save memory. +,,, +bev_dsm_restriction _value,--,float,Adds a lower state of charge (SOC) limit for battery electric vehicles (BEV) to manage its own energy demand (DSM). Located in `build_transport_demand.py `_. Set to 0 for no restriction on BEV DSM +bev_dsm_restriction _time,--,float,Time at which SOC of BEV has to be dsm_restriction_value +transport_heating _deadband_upper,°C,float,"The maximum temperature in the vehicle. At higher temperatures, the energy required for cooling in the vehicle increases." +transport_heating _deadband_lower,°C,float,"The minimum temperature in the vehicle. At lower temperatures, the energy required for heating in the vehicle increases." +,,, +ICE_lower_degree_factor,--,float,Share increase in energy demand in internal combustion engine (ICE) for each degree difference between the cold environment and the minimum temperature. +ICE_upper_degree_factor,--,float,Share increase in energy demand in internal combustion engine (ICE) for each degree difference between the hot environment and the maximum temperature. +EV_lower_degree_factor,--,float,Share increase in energy demand in electric vehicles (EV) for each degree difference between the cold environment and the minimum temperature. +EV_upper_degree_factor,--,float,Share increase in energy demand in electric vehicles (EV) for each degree difference between the hot environment and the maximum temperature. +bev_dsm,--,"{true, false}",Add the option for battery electric vehicles (BEV) to participate in demand-side management (DSM) +,,, +bev_availability,--,float,The share for battery electric vehicles (BEV) that are able to do demand side management (DSM) +bev_energy,--,float,The average size of battery electric vehicles (BEV) in MWh +bev_charge_efficiency,--,float,Battery electric vehicles (BEV) charge and discharge efficiency +bev_charge_rate,MWh,float,The power consumption for one electric vehicle (EV) in MWh. Value derived from 3-phase charger with 11 kW. +bev_avail_max,--,float,The maximum share plugged-in availability for passenger electric vehicles. +bev_avail_mean,--,float,The average share plugged-in availability for passenger electric vehicles. +v2g,--,"{true, false}",Allows feed-in to grid from EV battery +land_transport_fuel_cell _share,--,Dictionary with planning horizons as keys.,The share of vehicles that uses fuel cells in a given year +land_transport_electric _share,--,Dictionary with planning horizons as keys.,The share of vehicles that uses electric vehicles (EV) in a given year +land_transport_ice _share,--,Dictionary with planning horizons as keys.,The share of vehicles that uses internal combustion engines (ICE) in a given year. What is not EV or FCEV is oil-fuelled ICE. +transport_electric_efficiency,MWh/100km,float,The conversion efficiencies of electric vehicles in transport +transport_fuel_cell_efficiency,MWh/100km,float,The H2 conversion efficiencies of fuel cells in transport +transport_ice_efficiency,MWh/100km,float,The oil conversion efficiencies of internal combustion engine (ICE) in transport +agriculture_machinery _electric_share,--,float,The share for agricultural machinery that uses electricity +agriculture_machinery _oil_share,--,float,The share for agricultural machinery that uses oil +agriculture_machinery _fuel_efficiency,--,float,The efficiency of electric-powered machinery in the conversion of electricity to meet agricultural needs. +agriculture_machinery _electric_efficiency,--,float,The efficiency of oil-powered machinery in the conversion of oil to meet agricultural needs. +Mwh_MeOH_per_MWh_H2,LHV,float,"The energy amount of the produced methanol per energy amount of hydrogen. From `DECHEMA (2017) `_, page 64." +MWh_MeOH_per_tCO2,LHV,float,"The energy amount of the produced methanol per ton of CO2. From `DECHEMA (2017) `_, page 66." +MWh_MeOH_per_MWh_e,LHV,float,"The energy amount of the produced methanol per energy amount of electricity. From `DECHEMA (2017) `_, page 64." +shipping_hydrogen _liquefaction,--,"{true, false}",Whether to include liquefaction costs for hydrogen demand in shipping. +,,, +shipping_hydrogen_share,--,Dictionary with planning horizons as keys.,The share of ships powered by hydrogen in a given year +shipping_methanol_share,--,Dictionary with planning horizons as keys.,The share of ships powered by methanol in a given year +shipping_oil_share,--,Dictionary with planning horizons as keys.,The share of ships powered by oil in a given year +shipping_methanol _efficiency,--,float,The efficiency of methanol-powered ships in the conversion of methanol to meet shipping needs (propulsion). The efficiency increase from oil can be 10-15% higher according to the `IEA `_ +,,, +shipping_oil_efficiency,--,float,The efficiency of oil-powered ships in the conversion of oil to meet shipping needs (propulsion). Base value derived from 2011 +aviation_demand_factor,--,float,The proportion of demand for aviation compared to today's consumption +HVC_demand_factor,--,float,The proportion of demand for high-value chemicals compared to today's consumption +,,, +time_dep_hp_cop,--,"{true, false}",Consider the time dependent coefficient of performance (COP) of the heat pump +heat_pump_sink_T,°C,float,The temperature heat sink used in heat pumps based on DTU / large area radiators. The value is conservatively high to cover hot water and space heating in poorly-insulated buildings +reduce_space_heat _exogenously,--,"{true, false}",Influence on space heating demand by a certain factor (applied before losses in district heating). +reduce_space_heat _exogenously_factor,--,Dictionary with planning horizons as keys.,"A positive factor can mean renovation or demolition of a building. If the factor is negative, it can mean an increase in floor area, increased thermal comfort, population growth. The default factors are determined by the `Eurocalc Homes and buildings decarbonization scenario `_" +retrofitting,,, +-- retro_endogen,--,"{true, false}",Add retrofitting as an endogenous system which co-optimise space heat savings. +-- cost_factor,--,float,Weight costs for building renovation +-- interest_rate,--,float,The interest rate for investment in building components +-- annualise_cost,--,"{true, false}",Annualise the investment costs of retrofitting +-- tax_weighting,--,"{true, false}",Weight the costs of retrofitting depending on taxes in countries +-- construction_index,--,"{true, false}",Weight the costs of retrofitting depending on labour/material costs per country +tes,--,"{true, false}",Add option for storing thermal energy in large water pits associated with district heating systems and individual thermal energy storage (TES) +tes_tau,,,The time constant used to calculate the decay of thermal energy in thermal energy storage (TES): 1- :math:`e^{-1/24τ}`. +-- decentral,days,float,The time constant in decentralized thermal energy storage (TES) +-- central,days,float,The time constant in centralized thermal energy storage (TES) +boilers,--,"{true, false}",Add option for transforming gas into heat using gas boilers +resistive_heaters,--,"{true, false}",Add option for transforming electricity into heat using resistive heaters (independently from gas boilers) +oil_boilers,--,"{true, false}",Add option for transforming oil into heat using boilers +biomass_boiler,--,"{true, false}",Add option for transforming biomass into heat using boilers +overdimension_heat_generators,,,Add option for overdimensioning heating systems by a certain factor. This allows them to cover heat demand peaks e.g. 10% higher than those in the data with a setting of 1.1. +-- decentral,--,float,The factor for overdimensioning (increasing CAPEX) decentral heating systems +-- central,--,float,The factor for overdimensioning (increasing CAPEX) central heating systems +chp,--,"{true, false}",Add option for using Combined Heat and Power (CHP) +micro_chp,--,"{true, false}",Add option for using Combined Heat and Power (CHP) for decentral areas. +solar_thermal,--,"{true, false}",Add option for using solar thermal to generate heat. +solar_cf_correction,--,float,The correction factor for the value provided by the solar thermal profile calculations +marginal_cost_storage,"currency/MWh ",float,The marginal cost of discharging batteries in distributed grids +methanation,--,"{true, false}",Add option for transforming hydrogen and CO2 into methane using methanation. +coal_cc,--,"{true, false}",Add option for coal CHPs with carbon capture +dac,--,"{true, false}",Add option for Direct Air Capture (DAC) +co2_vent,--,"{true, false}",Add option for vent out CO2 from storages to the atmosphere. +allam_cycle,--,"{true, false}",Add option to include `Allam cycle gas power plants `_ +hydrogen_fuel_cell,--,"{true, false}",Add option to include hydrogen fuel cell for re-electrification. Assuming OCGT technology costs +hydrogen_turbine,--,"{true, false}",Add option to include hydrogen turbine for re-electrification. Assuming OCGT technology costs +SMR,--,"{true, false}",Add option for transforming natural gas into hydrogen and CO2 using Steam Methane Reforming (SMR) +SMR CC,--,"{true, false}",Add option for transforming natural gas into hydrogen and CO2 using Steam Methane Reforming (SMR) and Carbon Capture (CC) +regional_oil_demand,--,"{true, false}",Spatially resolve oil demand. Set to true if regional CO2 constraints needed. +regional_co2 _sequestration_potential,,, +-- enable,--,"{true, false}",Add option for regionally-resolved geological carbon dioxide sequestration potentials based on `CO2StoP `_. +-- attribute,--,string or list,Name (or list of names) of the attribute(s) for the sequestration potential +-- include_onshore,--,"{true, false}",Add options for including onshore sequestration potentials +-- min_size,Gt ,float,Any sites with lower potential than this value will be excluded +-- max_size,Gt ,float,The maximum sequestration potential for any one site. +-- years_of_storage,years,float,The years until potential exhausted at optimised annual rate +co2_sequestration_potential,--,Dictionary with planning horizons as keys.,The potential of sequestering CO2 in Europe per year and investment period +co2_sequestration_cost,currency/tCO2,float,The cost of sequestering a ton of CO2 +co2_sequestration_lifetime,years,int,The lifetime of a CO2 sequestration site +co2_spatial,--,"{true, false}","Add option to spatially resolve carrier representing stored carbon dioxide. This allows for more detailed modelling of CCUTS, e.g. regarding the capturing of industrial process emissions, usage as feedstock for electrofuels, transport of carbon dioxide, and geological sequestration sites." +,,, +co2network,--,"{true, false}",Add option for planning a new carbon dioxide transmission network +co2_network_cost_factor,p.u.,float,The cost factor for the capital cost of the carbon dioxide transmission network +,,, +cc_fraction,--,float,The default fraction of CO2 captured with post-combustion capture +hydrogen_underground _storage,--,"{true, false}",Add options for storing hydrogen underground. Storage potential depends regionally. +hydrogen_underground _storage_locations,,"{onshore, nearshore, offshore}","The location where hydrogen underground storage can be located. Onshore, nearshore, offshore means it must be located more than 50 km away from the sea, within 50 km of the sea, or within the sea itself respectively." +,,, +methanol,--,--,Add methanol as carrrier and add enabled methnol technologies +-- regional_methanol_demand,--,"{true, false}",Spatially resolve methanol demand. Set to true if regional CO2 constraints needed. +-- methanol_reforming,--,"{true, false}"," Add methanol reforming" +-- methanol_reforming_cc,--,"{true, false}"," Add methanol reforming with carbon capture" +-- methanol_to_kerosene,--,"{true, false}"," Add methanol to kerosene" +-- methanol_to_olefins,--,"{true, false}"," Add methanol to olefins" +-- methanol_to_power,--,--," Add different methanol to power technologies" +-- -- ccgt,--,"{true, false}"," Add combined cycle gas turbine (CCGT) using methanol" +-- -- ccgt_cc,--,"{true, false}"," Add combined cycle gas turbine (CCGT) with carbon capture using methanol" +-- -- ocgt,--,"{true, false}"," Add open cycle gas turbine (OCGT) using methanol" +-- -- allam,--,"{true, false}"," Add Allam cycle gas power plants using methanol" +,,, +ammonia,--,"{true, false, regional}","Add ammonia as a carrrier. It can be either true (copperplated NH3), false (no NH3 carrier) or ""regional"" (regionalised NH3 without network)" +min_part_load_fischer _tropsch,per unit of p_nom ,float,The minimum unit dispatch (``p_min_pu``) for the Fischer-Tropsch process +min_part_load _methanolisation,per unit of p_nom ,float,The minimum unit dispatch (``p_min_pu``) for the methanolisation process +,,, +use_fischer_tropsch _waste_heat,--,"{true, false}",Add option for using waste heat of Fischer Tropsch in district heating networks +use_fuel_cell_waste_heat,--,"{true, false}",Add option for using waste heat of fuel cells in district heating networks +use_electrolysis_waste _heat,--,"{true, false}",Add option for using waste heat of electrolysis in district heating networks +electricity_transmission _grid,--,"{true, false}",Switch for enabling/disabling the electricity transmission grid. +electricity_distribution _grid,--,"{true, false}",Add a simplified representation of the exchange capacity between transmission and distribution grid level through a link. +electricity_distribution _grid_cost_factor,,,Multiplies the investment cost of the electricity distribution grid +,,, +electricity_grid _connection,--,"{true, false}",Add the cost of electricity grid connection for onshore wind and solar +transmission_efficiency,,,Section to specify transmission losses or compression energy demands of bidirectional links. Splits them into two capacity-linked unidirectional links. +-- {carrier},--,str,The carrier of the link. +-- -- efficiency_static,p.u.,float,Length-independent transmission efficiency. +-- -- efficiency_per_1000km,p.u. per 1000 km,float,Length-dependent transmission efficiency ($\eta^{\text{length}}$) +-- -- compression_per_1000km,p.u. per 1000 km,float,Length-dependent electricity demand for compression ($\eta \cdot \text{length}$) implemented as multi-link to local electricity bus. +H2_network,--,"{true, false}",Add option for new hydrogen pipelines +gas_network,--,"{true, false}","Add existing natural gas infrastructure, incl. LNG terminals, production and entry-points. The existing gas network is added with a lossless transport model. A length-weighted `k-edge augmentation algorithm `_ can be run to add new candidate gas pipelines such that all regions of the model can be connected to the gas network. When activated, all the gas demands are regionally disaggregated as well." +H2_retrofit,--,"{true, false}",Add option for retrofiting existing pipelines to transport hydrogen. +H2_retrofit_capacity _per_CH4,--,float,"The ratio for H2 capacity per original CH4 capacity of retrofitted pipelines. The `European Hydrogen Backbone (April, 2020) p.15 `_ 60% of original natural gas capacity could be used in cost-optimal case as H2 capacity." +"gas_network_connectivity _upgrade ",--,float,The number of desired edge connectivity (k) in the length-weighted `k-edge augmentation algorithm `_ used for the gas network +gas_distribution_grid,--,"{true, false}",Add a gas distribution grid +gas_distribution_grid _cost_factor,,,Multiplier for the investment cost of the gas distribution grid +,,, +biomass_spatial,--,"{true, false}",Add option for resolving biomass demand regionally +biomass_transport,--,"{true, false}",Add option for transporting solid biomass between nodes +biogas_upgrading_cc,--,"{true, false}",Add option to capture CO2 from biomass upgrading +conventional_generation,,,Add a more detailed description of conventional carriers. Any power generation requires the consumption of fuel from nodes representing that fuel. +biomass_to_liquid,--,"{true, false}",Add option for transforming solid biomass into liquid fuel with the same properties as oil +biosng,--,"{true, false}",Add option for transforming solid biomass into synthesis gas with the same properties as natural gas +bioH2,--,"{true, false}",Add option for transforming solid biomass into hydrogen with carbon capture +municipal_solid_waste,--,"{true, false}",Add option for municipal solid waste +limit_max_growth,,, +-- enable,--,"{true, false}",Add option to limit the maximum growth of a carrier +-- factor,p.u.,float,The maximum growth factor of a carrier (e.g. 1.3 allows 30% larger than max historic growth) +-- max_growth,,, +-- -- {carrier},GW,float,The historic maximum growth of a carrier +-- max_relative_growth,,, +-- -- {carrier},p.u.,float,The historic maximum relative growth of a carrier +,,, +enhanced_geothermal,,, +-- enable,--,"{true, false}",Add option to include Enhanced Geothermal Systems +-- flexible,--,"{true, false}",Add option for flexible operation (see Ricks et al. 2024) +-- max_hours,--,int,The maximum hours the reservoir can be charged under flexible operation +-- max_boost,--,float,The maximum boost in power output under flexible operation +-- var_cf,--,"{true, false}",Add option for variable capacity factor (see Ricks et al. 2024) +-- sustainability_factor,--,float,Share of sourced heat that is replenished by the earth's core (see details in `build_egs_potentials.py `_) +solid_biomass_import,,, +-- enable,--,"{true, false}",Add option to include solid biomass imports +-- price,currency/MWh,float,Price for importing solid biomass +-- max_amount,Twh,float,Maximum solid biomass import potential +-- upstream_emissions_factor,p.u.,float,Upstream emissions of solid biomass imports diff --git a/doc/release_notes.rst b/doc/release_notes.rst index 7340fe316..a3411e87f 100644 --- a/doc/release_notes.rst +++ b/doc/release_notes.rst @@ -10,9 +10,13 @@ Release Notes .. Upcoming Release + +* Add technology options for methanol, like electricity production from methanol, biomass to methanol, methanol to kerosene, ... + * Change the heating demand from final energy which includes losses in legacy equipment to thermal energy service based on JRC-IDEES. Efficiencies of existing heating capacities are lowered according to the conversion of final energy to thermal energy service. For overnight scenarios or future planning horizon this change leads to a reduction in heat supply. * Updated district heating supply temperatures based on `Euroheat's DHC Market Outlook 2024`__ and `AGFW-Hauptbericht 2022 `__. `min_forward_temperature` and `return_temperature` (not given by Euroheat) are extrapolated based on German values. + * Made the overdimensioning factor for heating systems specific for central/decentral heating, defaults to no overdimensionining for central heating and no changes to decentral heating compared to previous version. * bugfix: The carrier of stores was silently overwritten by their bus_carrier as a side effect when building the co2 constraints diff --git a/rules/build_sector.smk b/rules/build_sector.smk index 60db6c144..b3cfc4afa 100755 --- a/rules/build_sector.smk +++ b/rules/build_sector.smk @@ -1088,6 +1088,9 @@ rule prepare_sector_network: hourly_heat_demand_total=resources( "hourly_heat_demand_total_elec_s{simpl}_{clusters}.nc" ), + industrial_production=resources( + "industrial_production_elec_s{simpl}_{clusters}_{planning_horizons}.csv" + ), district_heat_share=resources( "district_heat_share_elec_s{simpl}_{clusters}_{planning_horizons}.csv" ), diff --git a/scripts/prepare_sector_network.py b/scripts/prepare_sector_network.py index 4ddcac327..4a9c1be6b 100755 --- a/scripts/prepare_sector_network.py +++ b/scripts/prepare_sector_network.py @@ -157,7 +157,7 @@ def define_spatial(nodes, options): spatial.methanol.nodes = ["EU methanol"] spatial.methanol.locations = ["EU"] - if options["regional_methanol_demand"]: + if options["methanol"]["regional_methanol_demand"]: spatial.methanol.demand_locations = nodes spatial.methanol.industry = nodes + " industry methanol" spatial.methanol.shipping = nodes + " shipping methanol" @@ -819,11 +819,315 @@ def add_allam(n, costs): carrier="allam", p_nom_extendable=True, # TODO: add costs to technology-data - capital_cost=0.6 * 1.5e6 * 0.1, # efficiency * EUR/MW * annuity - marginal_cost=2, - efficiency=0.6, + capital_cost=costs.at["allam", "fixed"] * costs.at["allam", "efficiency"], + marginal_cost=costs.at["allam", "VOM"] * costs.at["allam", "efficiency"], + efficiency=costs.at["allam", "efficiency"], efficiency2=costs.at["gas", "CO2 intensity"], - lifetime=30.0, + lifetime=costs.at["allam", "lifetime"], + ) + + +def add_biomass_to_methanol(n, costs): + + n.madd( + "Link", + spatial.biomass.nodes, + suffix=" biomass-to-methanol", + bus0=spatial.biomass.nodes, + bus1=spatial.methanol.nodes, + bus2="co2 atmosphere", + carrier="biomass-to-methanol", + lifetime=costs.at["biomass-to-methanol", "lifetime"], + efficiency=costs.at["biomass-to-methanol", "efficiency"], + efficiency2=-costs.at["solid biomass", "CO2 intensity"] + + costs.at["biomass-to-methanol", "CO2 stored"], + p_nom_extendable=True, + capital_cost=costs.at["biomass-to-methanol", "fixed"] + / costs.at["biomass-to-methanol", "efficiency"], + marginal_cost=costs.loc["biomass-to-methanol", "VOM"] + / costs.at["biomass-to-methanol", "efficiency"], + ) + + +def add_biomass_to_methanol_cc(n, costs): + + n.madd( + "Link", + spatial.biomass.nodes, + suffix=" biomass-to-methanol CC", + bus0=spatial.biomass.nodes, + bus1=spatial.methanol.nodes, + bus2="co2 atmosphere", + bus3=spatial.co2.nodes, + carrier="biomass-to-methanol CC", + lifetime=costs.at["biomass-to-methanol", "lifetime"], + efficiency=costs.at["biomass-to-methanol", "efficiency"], + efficiency2=-costs.at["solid biomass", "CO2 intensity"] + + costs.at["biomass-to-methanol", "CO2 stored"] + * (1 - costs.at["biomass-to-methanol", "capture rate"]), + efficiency3=costs.at["biomass-to-methanol", "CO2 stored"] + * costs.at["biomass-to-methanol", "capture rate"], + p_nom_extendable=True, + capital_cost=costs.at["biomass-to-methanol", "fixed"] + / costs.at["biomass-to-methanol", "efficiency"] + + costs.at["biomass CHP capture", "fixed"] + * costs.at["biomass-to-methanol", "CO2 stored"], + marginal_cost=costs.loc["biomass-to-methanol", "VOM"] + / costs.at["biomass-to-methanol", "efficiency"], + ) + + +def add_methanol_to_power(n, costs, types={}): + # TODO: add costs to technology-data + + nodes = pop_layout.index + + if types["allam"]: + logger.info("Adding Allam cycle methanol power plants.") + + n.madd( + "Link", + nodes, + suffix=" allam methanol", + bus0=spatial.methanol.nodes, + bus1=nodes, + bus2=spatial.co2.df.loc[nodes, "nodes"].values, + bus3="co2 atmosphere", + carrier="allam methanol", + p_nom_extendable=True, + capital_cost=costs.at["allam", "fixed"] * costs.at["allam", "efficiency"], + marginal_cost=costs.at["allam", "VOM"] * costs.at["allam", "efficiency"], + efficiency=costs.at["allam", "efficiency"], + efficiency2=0.98 * costs.at["methanolisation", "carbondioxide-input"], + efficiency3=0.02 * costs.at["methanolisation", "carbondioxide-input"], + lifetime=25, + ) + + if types["ccgt"]: + logger.info("Adding methanol CCGT power plants.") + + # efficiency * EUR/MW * (annuity + FOM) + capital_cost = costs.at["CCGT", "efficiency"] * costs.at["CCGT", "fixed"] + + n.madd( + "Link", + nodes, + suffix=" CCGT methanol", + bus0=spatial.methanol.nodes, + bus1=nodes, + bus2="co2 atmosphere", + carrier="CCGT methanol", + p_nom_extendable=True, + capital_cost=capital_cost, + marginal_cost=costs.at["CCGT", "VOM"], + efficiency=costs.at["CCGT", "efficiency"], + efficiency2=costs.at["methanolisation", "carbondioxide-input"], + lifetime=costs.at["CCGT", "lifetime"], + ) + + if types["ccgt_cc"]: + logger.info( + "Adding methanol CCGT power plants with post-combustion carbon capture." + ) + + # TODO consider efficiency changes / energy inputs for CC + + # efficiency * EUR/MW * (annuity + FOM) + capital_cost = costs.at["CCGT", "efficiency"] * costs.at["CCGT", "fixed"] + + capital_cost_cc = ( + capital_cost + + costs.at["cement capture", "fixed"] + * costs.at["methanolisation", "carbondioxide-input"] + ) + + n.madd( + "Link", + nodes, + suffix=" CCGT methanol CC", + bus0=spatial.methanol.nodes, + bus1=nodes, + bus2=spatial.co2.df.loc[nodes, "nodes"].values, + bus3="co2 atmosphere", + carrier="CCGT methanol CC", + p_nom_extendable=True, + capital_cost=capital_cost_cc, + marginal_cost=costs.at["CCGT", "VOM"], + efficiency=costs.at["CCGT", "efficiency"], + efficiency2=costs.at["cement capture", "capture_rate"] + * costs.at["methanolisation", "carbondioxide-input"], + efficiency3=(1 - costs.at["cement capture", "capture_rate"]) + * costs.at["methanolisation", "carbondioxide-input"], + lifetime=costs.at["CCGT", "lifetime"], + ) + + if types["ocgt"]: + logger.info("Adding methanol OCGT power plants.") + + n.madd( + "Link", + nodes, + suffix=" OCGT methanol", + bus0=spatial.methanol.nodes, + bus1=nodes, + bus2="co2 atmosphere", + carrier="OCGT methanol", + p_nom_extendable=True, + capital_cost=costs.at["OCGT", "fixed"] * costs.at["OCGT", "efficiency"], + marginal_cost=costs.at["OCGT", "VOM"] * costs.at["OCGT", "efficiency"], + efficiency=costs.at["OCGT", "efficiency"], + efficiency2=costs.at["methanolisation", "carbondioxide-input"], + lifetime=costs.at["OCGT", "lifetime"], + ) + + +def add_methanol_to_olefins(n, costs): + nodes = spatial.nodes + nhours = n.snapshot_weightings.generators.sum() + nyears = nhours / 8760 + + tech = "methanol-to-olefins/aromatics" + + logger.info(f"Adding {tech}.") + + demand_factor = options["HVC_demand_factor"] + + industrial_production = ( + pd.read_csv(snakemake.input.industrial_production, index_col=0) + * 1e3 + * nyears # kt/a -> t/a + ) + + p_nom_max = ( + demand_factor + * industrial_production.loc[nodes, "HVC"] + / nhours + * costs.at[tech, "methanol-input"] + ) + + co2_release = ( + costs.at[tech, "carbondioxide-output"] / costs.at[tech, "methanol-input"] + + costs.at["methanolisation", "carbondioxide-input"] + ) + + n.madd( + "Link", + nodes, + suffix=f" {tech}", + carrier=tech, + capital_cost=costs.at[tech, "fixed"] / costs.at[tech, "methanol-input"], + marginal_cost=costs.at[tech, "VOM"] / costs.at[tech, "methanol-input"], + p_nom_extendable=True, + bus0=spatial.methanol.nodes, + bus1=spatial.oil.naphtha, + bus2=nodes, + bus3="co2 atmosphere", + p_min_pu=1, + p_nom_max=p_nom_max.values, + efficiency=1 / costs.at[tech, "methanol-input"], + efficiency2=-costs.at[tech, "electricity-input"] + / costs.at[tech, "methanol-input"], + efficiency3=co2_release, + ) + + +def add_methanol_to_kerosene(n, costs): + nodes = pop_layout.index + nhours = n.snapshot_weightings.generators.sum() + + demand_factor = options["aviation_demand_factor"] + + tech = "methanol-to-kerosene" + + logger.info(f"Adding {tech}.") + + all_aviation = ["total international aviation", "total domestic aviation"] + + p_nom_max = ( + demand_factor + * pop_weighted_energy_totals.loc[nodes, all_aviation].sum(axis=1) + * 1e6 + / nhours + * costs.at[tech, "methanol-input"] + ) + + capital_cost = costs.at[tech, "fixed"] / costs.at[tech, "methanol-input"] + + n.madd( + "Link", + spatial.h2.locations, + suffix=f" {tech}", + carrier=tech, + capital_cost=capital_cost, + bus0=spatial.methanol.nodes, + bus1=spatial.oil.kerosene, + bus2=spatial.h2.nodes, + efficiency=costs.at[tech, "methanol-input"], + efficiency2=-costs.at[tech, "hydrogen-input"] + / costs.at[tech, "methanol-input"], + p_nom_extendable=True, + p_min_pu=1, + p_nom_max=p_nom_max.values, + ) + + +def add_methanol_reforming(n, costs): + logger.info("Adding methanol steam reforming.") + + tech = "Methanol steam reforming" + + capital_cost = costs.at[tech, "fixed"] / costs.at[tech, "methanol-input"] + + n.madd( + "Link", + spatial.h2.locations, + suffix=f" {tech}", + bus0=spatial.methanol.nodes, + bus1=spatial.h2.nodes, + bus2="co2 atmosphere", + p_nom_extendable=True, + capital_cost=capital_cost, + efficiency=1 / costs.at[tech, "methanol-input"], + efficiency2=costs.at["methanolisation", "carbondioxide-input"], + carrier=tech, + lifetime=costs.at[tech, "lifetime"], + ) + + +def add_methanol_reforming_cc(n, costs): + logger.info("Adding methanol steam reforming with carbon capture.") + + tech = "Methanol steam reforming" + + # TODO: heat release and electricity demand for process and carbon capture + # but the energy demands for carbon capture have not yet been added for other CC processes + # 10.1016/j.rser.2020.110171: 0.129 kWh_e/kWh_H2, -0.09 kWh_heat/kWh_H2 + + capital_cost = costs.at[tech, "fixed"] / costs.at[tech, "methanol-input"] + + capital_cost_cc = ( + capital_cost + + costs.at["cement capture", "fixed"] + * costs.at["methanolisation", "carbondioxide-input"] + ) + + n.madd( + "Link", + spatial.h2.locations, + suffix=f" {tech} CC", + bus0=spatial.methanol.nodes, + bus1=spatial.h2.nodes, + bus2="co2 atmosphere", + bus3=spatial.co2.nodes, + p_nom_extendable=True, + capital_cost=capital_cost_cc, + efficiency=1 / costs.at[tech, "methanol-input"], + efficiency2=(1 - costs.at["cement capture", "capture_rate"]) + * costs.at["methanolisation", "carbondioxide-input"], + efficiency3=costs.at["cement capture", "capture_rate"] + * costs.at["methanolisation", "carbondioxide-input"], + carrier=f"{tech} CC", + lifetime=costs.at[tech, "lifetime"], ) @@ -2309,6 +2613,32 @@ def add_heat(n: pypsa.Network, costs: pd.DataFrame, cop: xr.DataArray): ) +def add_methanol(n, costs): + + methanol_options = options["methanol"] + if not any(methanol_options.values()): + return + + logger.info("Add methanol") + add_carrier_buses(n, "methanol") + + if options["biomass"]: + if methanol_options["biomass_to_methanol"]: + add_biomass_to_methanol(n, costs) + + if methanol_options["biomass_to_methanol"]: + add_biomass_to_methanol_cc(n, costs) + + if methanol_options["methanol_to_power"]: + add_methanol_to_power(n, costs, types=methanol_options["methanol_to_power"]) + + if methanol_options["methanol_reforming"]: + add_methanol_reforming(n, costs) + + if methanol_options["methanol_reforming_cc"]: + add_methanol_reforming_cc(n, costs) + + def add_biomass(n, costs): logger.info("Add biomass") @@ -3129,7 +3459,7 @@ def add_industry(n, costs): / nhours ) - if not options["regional_methanol_demand"]: + if not options["methanol"]["regional_methanol_demand"]: p_set_methanol = p_set_methanol.sum() n.madd( @@ -3252,7 +3582,7 @@ def add_industry(n, costs): * efficiency ) - if not options["regional_methanol_demand"]: + if not options["methanol"]["regional_methanol_demand"]: p_set_methanol_shipping = p_set_methanol_shipping.sum() n.madd( @@ -3526,6 +3856,9 @@ def add_industry(n, costs): efficiency3=process_co2_per_naphtha, ) + if options["methanol"]["methanol_to_olefins"]: + add_methanol_to_olefins(n, costs) + # aviation demand_factor = options.get("aviation_demand_factor", 1) if demand_factor != 1: @@ -3570,6 +3903,9 @@ def add_industry(n, costs): efficiency2=costs.at["oil", "CO2 intensity"], ) + if options["methanol"]["methanol_to_kerosene"]: + add_methanol_to_kerosene(n, costs) + # TODO simplify bus expression n.madd( "Load", @@ -4413,6 +4749,9 @@ def add_enhanced_geothermal(n, egs_potentials, egs_overlap, costs): if options["ammonia"]: add_ammonia(n, costs) + if options["methanol"]: + add_methanol(n, costs) + if options["industry"]: add_industry(n, costs)