From 886b7c5f485e6bc1642f7361c9c99b9014df2428 Mon Sep 17 00:00:00 2001
From: Tom Brown <tom@nworbmot.org>
Date: Sun, 28 Jul 2019 13:10:21 +0200
Subject: [PATCH] Include micro-CHP using PEMFC

---
 config.yaml                       |  6 +++---
 data/costs.csv                    |  5 +++++
 scripts/prepare_sector_network.py | 31 ++++++++++++++++++++++++++++++-
 3 files changed, 38 insertions(+), 4 deletions(-)

diff --git a/config.yaml b/config.yaml
index 0aec2ff5..9e3069f7 100644
--- a/config.yaml
+++ b/config.yaml
@@ -2,15 +2,15 @@ logging_level: INFO
 
 results_dir: 'results/'
 summary_dir: results
-run: '190720-new_industry'
+run: '190726-microchp'
 
 scenario:
   sectors: [E] # ,E+EV,E+BEV,E+BEV+V2G] # [ E+EV, E+BEV, E+BEV+V2G ]
   simpl: ['']
-  lv: [1.0,1.25,1.5,2.0]#, 1.125, 1.25, 1.5, 2.0]# or opt
+  lv: [1.0]#, 1.125, 1.25, 1.5, 2.0]# or opt
   clusters: [128] #[90, 128, 181] #[45, 64, 90, 128, 181, 256] #, 362] # (2**np.r_[5.5:9:.5]).astype(int) minimum is 37
   opts: ['']  #for pypsa-eur
-  sector_opts: [Co2L0-3H-T-H-B-I]#,Co2L0p1-3H-T-H-B-I,Co2L0p25-3H-T-H-B-I,Co2L0p5-3H-T-H-B-I]#[Co2L0-3H-T-H-B-I-onwind0-solar3,Co2L0-3H-T-H-B-I-onwind0p125-solar3,Co2L0-3H-T-H-B-I-onwind0p25-solar3,Co2L0-3H-T-H-B-I-onwind0p50-solar3,Co2L0-3H-T-H-B-I-solar3]#,Co2L0-3H-T-H-B-I-onwind0p25-solar3]#,Co2L0p05-3H-T-H-B-I,Co2L0p10-3H-T-H-B-I,Co2L0p20-3H-T-H-B-I,Co2L0p30-3H-T-H-B-I,Co2L0p50-3H-T-H-B-I]#[Co2L-3H-T-H,Co2L0p10-3H-T-H,Co2L0-3H-T-H,Co2L0p20-3H-T-H]   #Co2L-3H-T-H,Co2L0p10-3H-T-H,Co2L0p20-3H-T-HCo2L-3H-T-H,Co2L0p10-3H-T-H,Co2L0p30-3H-T-H,Co2L0p50-3H-T-H] #Co2L-3H,Co2L-3H-T,, LC-FL, LC-T, Ep-T, Co2L-T]
+  sector_opts: [Co2L0-3H-T-H-B-I,Co2L0-3H-T-H-B-I-nodistrict]#,Co2L0p1-3H-T-H-B-I,Co2L0p25-3H-T-H-B-I,Co2L0p5-3H-T-H-B-I]#[Co2L0-3H-T-H-B-I-onwind0-solar3,Co2L0-3H-T-H-B-I-onwind0p125-solar3,Co2L0-3H-T-H-B-I-onwind0p25-solar3,Co2L0-3H-T-H-B-I-onwind0p50-solar3,Co2L0-3H-T-H-B-I-solar3]#,Co2L0-3H-T-H-B-I-onwind0p25-solar3]#,Co2L0p05-3H-T-H-B-I,Co2L0p10-3H-T-H-B-I,Co2L0p20-3H-T-H-B-I,Co2L0p30-3H-T-H-B-I,Co2L0p50-3H-T-H-B-I]#[Co2L-3H-T-H,Co2L0p10-3H-T-H,Co2L0-3H-T-H,Co2L0p20-3H-T-H]   #Co2L-3H-T-H,Co2L0p10-3H-T-H,Co2L0p20-3H-T-HCo2L-3H-T-H,Co2L0p10-3H-T-H,Co2L0p30-3H-T-H,Co2L0p50-3H-T-H] #Co2L-3H,Co2L-3H-T,, LC-FL, LC-T, Ep-T, Co2L-T]
   # Co2L will give default (5%); Co2L0p25 will give 25% CO2 emissions; Co2Lm0p05 will give 5% negative emissions
 
 
diff --git a/data/costs.csv b/data/costs.csv
index c866b09e..b414c0ec 100644
--- a/data/costs.csv
+++ b/data/costs.csv
@@ -175,6 +175,11 @@ decentral CHP,2030,discount rate,0.04,per unit,Palzer thesis
 central CHP,2030,lifetime,25,years,HP
 central CHP,2030,investment,650,EUR/kWel,HP
 central CHP,2030,FOM,3,%/year,HP
+micro CHP,2030,lifetime,20,years,DEA for PEMFC with methane (for unit consuming 2kW CH4)
+micro CHP,2030,investment,4500,EUR/kWCH4,DEA for PEMFC with methane (for unit consuming 2kW CH4)
+micro CHP,2030,FOM,6,%/year,DEA for PEMFC with methane (for unit consuming 2kW CH4)
+micro CHP,2030,efficiency,0.351,per unit,DEA for PEMFC with methane (for unit consuming 2kW CH4)
+micro CHP,2030,efficiency-heat,0.609,per unit,DEA for PEMFC with methane (for unit consuming 2kW CH4)
 decentral solar thermal,2030,discount rate,0.04,per unit,Palzer thesis
 decentral solar thermal,2030,FOM,1.3,%/year,HP
 decentral solar thermal,2030,investment,270000,EUR/1000m2,HP
diff --git a/scripts/prepare_sector_network.py b/scripts/prepare_sector_network.py
index 130622f3..d450e7d7 100644
--- a/scripts/prepare_sector_network.py
+++ b/scripts/prepare_sector_network.py
@@ -978,7 +978,7 @@ def add_heat(network):
                      efficiency=costs.at['central resistive heater','efficiency'])
 
         network.madd("Link",
-                     rural + " gas boiler",
+                     rural + " rural gas boiler",
                      p_nom_extendable=True,
                      bus0=["EU gas"]*len(rural),
                      bus1=rural + " rural heat",
@@ -1010,6 +1010,35 @@ def add_heat(network):
                      efficiency2=costs.at['gas','CO2 intensity'],
                      efficiency=costs.at['central gas boiler','efficiency'])
 
+        network.madd("Link",
+                     rural + " rural micro CHP",
+                     p_nom_extendable=True,
+                     bus0=["EU gas"]*len(rural),
+                     bus1=rural,
+                     bus2=rural + " rural heat",
+                     bus3="co2 atmosphere",
+                     carrier="rural micro CHP",
+                     efficiency=costs.at['micro CHP','efficiency'],
+                     efficiency2=costs.at['micro CHP','efficiency-heat'],
+                     efficiency3=costs.at['gas','CO2 intensity'],
+                     capital_cost=costs.at['micro CHP','fixed'])
+
+
+        network.madd("Link",
+                     urban_decentral + " urban decentral micro CHP",
+                     p_nom_extendable=True,
+                     bus0=["EU gas"]*len(urban_decentral),
+                     bus1=urban_decentral,
+                     bus2=urban_decentral + " urban decentral heat",
+                     bus3="co2 atmosphere",
+                     carrier="urban decentral micro CHP",
+                     efficiency=costs.at['micro CHP','efficiency'],
+                     efficiency2=costs.at['micro CHP','efficiency-heat'],
+                     efficiency3=costs.at['gas','CO2 intensity'],
+                     capital_cost=costs.at['micro CHP','fixed'])
+
+
+
     if options["chp"]:
 
         #additional bus, to which we can also connect biomass