Battery electric vehicle

AUTHORS.rst

@@ -13,4 +13,5 @@ Authors
 * Sarah Berendes
 * Marie-Claire Gering
 * Julian Endres
+* Sabine Haas
 * Felix Maurer

CHANGELOG.rst

@@ -9,6 +9,9 @@ Unreleased
 
 Features
 
+* Add tests for BEV facades developed in #94 `#142 <https://github.com/oemof/oemof-tabular/pull/142>`_
+
+
 Fixes
 
 

src/oemof/tabular/_facade.py

@@ -151,7 +151,7 @@ def __init__(self, *args, **kwargs):
             add_subnodes, sender=self
         )
 
-    def _nominal_value(self):
+    def _nominal_value(self, value=None):
         """Returns None if self.expandable ist True otherwise it returns
         the capacity
         """
@@ -168,18 +168,23 @@ def _nominal_value(self):
                     "to_from": self.to_from_capacity,
                 }
             else:
-                return self.capacity
+                if value:
+                    return value
+                else:
+                    return self.capacity
 
-    def _investment(self):
+    def _investment(self, bev=False):
         if not self.expandable:
             self.investment = None
             return self.investment
+
         if self.capacity_cost is None:
             msg = (
                 "If you set `expandable`to True you need to set "
                 "attribute `capacity_cost` of component {}!"
             )
             raise ValueError(msg.format(self.label))
+
         # If storage component
         if isinstance(self, GenericStorage):
             # If invest costs/MWH are given
@@ -207,7 +212,7 @@ def _investment(self):
                     age=getattr(self, "age", 0),
                     fixed_costs=getattr(self, "fixed_costs", None),
                 )
-        # If other component than storage
+        # If other component than storage or Bev
         else:
             self.investment = Investment(
                 ep_costs=self.capacity_cost,

src/oemof/tabular/constraint_facades.py

@@ -2,6 +2,38 @@
 from dataclasses import dataclass
 
 from oemof.solph.constraints.integral_limit import generic_integral_limit
+from pyomo.environ import Constraint
+
+from oemof.tabular.facades import Bev
+
+
+def var2str(var):
+    return "_".join(
+        [i.label if not isinstance(i, int) else str(i) for i in var]
+    )
+
+
+def get_bev_label(var):
+    return var[1].label.split("-storage")[0]
+
+
+def get_period(model, year, constraint_type=None):
+    if model.es.periods:
+        years = [period.year.min() for period in model.es.periods]
+        for period_index, period_year in enumerate(years):
+            if period_year == year:
+                return period_index
+        raise ValueError(
+            f"'{constraint_type}' constraint facade:\n"
+            f"Year {year} is not in model.PERIODS."
+        )
+    elif year == model.es.timeindex.year.min():
+        return 0
+    else:
+        raise ValueError(
+            f"'{constraint_type}' constraint facade:\n"
+            f"Year {year} is not in model.timeindex."
+        )
 
 
 class ConstraintFacade(abc.ABC):
@@ -39,4 +71,150 @@ def build_constraint(self, model):
         )
 
 
-CONSTRAINT_TYPE_MAP = {"generic_integral_limit": GenericIntegralLimit}
+@dataclass
+class BevShareMob(ConstraintFacade):
+    # TODO: rework docstring
+    """This constraint is only feasible if the definition of one vehicle is the
+    same (charging_capacity/storage_capacity) for all three bev technologies"""
+    name: str
+    type: str
+    year: int
+    label: str
+    share_mob_flex_G2V: (int, float)
+    share_mob_flex_V2G: (int, float)
+    share_mob_inflex: (int, float)
+
+    @staticmethod
+    def map_share2vars(model, share_mob, period):
+        invest_vars = []
+        for node in model.es.nodes:
+            if isinstance(node, Bev):
+                invest_vars.extend(
+                    [
+                        inv
+                        for inv in model.InvestmentFlowBlock.invest
+                        if node.electricity_bus.label == inv[0].label
+                        and f"{node.facade_label}-storage" in inv[1].label
+                        and period == inv[2]
+                    ]
+                )
+
+        share_mob = {
+            inv_var: value
+            for key, value in share_mob.items()
+            for inv_var in invest_vars
+            if key in inv_var[1].label
+        }
+        return invest_vars, share_mob
+
+    @staticmethod
+    def convert_share(share):
+        if 0 <= share <= 1:
+            return share
+        elif 0 <= share <= 100:
+            return share / 100
+        else:
+            raise ValueError(f"Mob share: {share} not in [0,1] or [0,100]")
+
+    def build_constraint(self, model):
+        period = get_period(model, self.year, self.__class__.__name__)
+
+        if period > len(model.PERIODS):
+            raise ValueError(
+                f"Period {period} not in model.PERIODS {model.PERIODS}"
+            )
+
+        share_mob = {
+            f"{self.label}-G2V-storage": self.convert_share(
+                self.share_mob_flex_G2V
+            ),
+            f"{self.label}-V2G-storage": self.convert_share(
+                self.share_mob_flex_V2G
+            ),
+            f"{self.label}-inflex-storage": self.convert_share(
+                self.share_mob_inflex
+            ),
+        }
+
+        invest_vars, share_mob = self.map_share2vars(
+            model=model, share_mob=share_mob, period=period
+        )
+
+        def investment_constraints_rule(InvestmentFlowBlock):
+            return InvestmentFlowBlock.invest[inv_var] == share_mob[
+                inv_var
+            ] * sum(InvestmentFlowBlock.invest[iv] for iv in invest_vars)
+
+        for inv_var in invest_vars:
+            name = f"mob_share_{get_bev_label(inv_var)}_{period}"
+            setattr(
+                model.InvestmentFlowBlock,
+                name,
+                Constraint(rule=investment_constraints_rule),
+            )
+
+
+# TODO maybe move inside facade
+@dataclass
+class BevEqualInvest(ConstraintFacade):
+    name: str
+    type: str
+    year: int
+
+    @staticmethod
+    def double_with_offset(lst):
+        result = []
+        for i in range(len(lst) - 1):
+            result.append((lst[i], lst[i + 1]))
+        return result
+
+    @staticmethod
+    def get_bev_invest_vars(model, period):
+        all_invest_vars = {}
+        for node in model.es.nodes:
+            if isinstance(node, Bev):
+                invest_vars_bev = list(
+                    set(
+                        inv
+                        for inv in model.InvestmentFlowBlock.invest
+                        if inv[2] == period
+                        for edge in inv[:2]
+                        if node.facade_label in edge.label
+                    )
+                )
+                all_invest_vars[node.facade_label] = invest_vars_bev
+        return all_invest_vars
+
+    def build_constraint(self, model):
+        # TODO add checks
+        period = get_period(model, self.year, self.__class__.__name__)
+        if period > len(model.PERIODS):
+            raise ValueError(
+                f"Period {period} not in model.PERIODS {model.PERIODS}"
+            )
+
+        invest_vars = self.get_bev_invest_vars(model, period)
+
+        def equate_variables_rule(InvestmentFlowBlock):
+            return (
+                InvestmentFlowBlock.invest[var1]
+                == InvestmentFlowBlock.invest[var2]
+            )
+
+        for bev_label, invest_vars in invest_vars.items():
+            for i, (var1, var2) in enumerate(
+                self.double_with_offset(invest_vars)
+            ):
+                name = f"{bev_label}_equal_invest_{i}_({period})"
+                setattr(
+                    model.InvestmentFlowBlock,
+                    name,
+                    Constraint(rule=equate_variables_rule),
+                )
+
+
+CONSTRAINT_TYPE_MAP = {
+    "generic_integral_limit": GenericIntegralLimit,
+    "bev_equal_invest": BevEqualInvest,
+    "bev_share_mob": BevShareMob,
+}

src/oemof/tabular/examples/datapackages/private_transport/README.md

@@ -0,0 +1,5 @@
+# Emission constraint example for oemof-tabular
+
+Run `scripts/infer.py` from the datapackage root directory to add the
+meta data file `datapackage.json` after updating the resources of the
+datapackage.

src/oemof/tabular/examples/datapackages/private_transport/data/constraints/bev_mob_share.csv

@@ -0,0 +1,3 @@
+name;type;year;label;share_mob_flex_G2V;share_mob_flex_V2G;share_mob_inflex
+bev_share_mob_2020;bev_share_mob;2020;BEV;10;20;70
+bev_share_mob_2030;bev_share_mob;2030;BEV;20;30;50

src/oemof/tabular/examples/datapackages/private_transport/data/constraints/bev_total_invest.csv

@@ -0,0 +1,3 @@
+name;type;year
+bev_equal_invest;bev_equal_invest;2020
+bev_equal_invest;bev_equal_invest;2030

src/oemof/tabular/examples/datapackages/private_transport/data/elements/bus.csv

@@ -0,0 +1,2 @@
+name;type;balanced
+bus-electricity;bus;true

src/oemof/tabular/examples/datapackages/private_transport/data/elements/dispatchable.csv

@@ -0,0 +1,4 @@
+name;type;carrier;tech;capacity;bus;marginal_cost;profile;output_parameters
+gas;dispatchable;gas;gt;1000;bus-electricity;40;1;{"custom_attributes": {"emission_factor": 10}}
+coal;dispatchable;coal;st;1000;bus-electricity;40;1;{"custom_attributes": {"emission_factor": 20}}
+lignite;dispatchable;lignite;st;500;bus-electricity;20;1;{"custom_attributes": {"emission_factor": 30}}

src/oemof/tabular/examples/datapackages/private_transport/data/elements/excess.csv

@@ -0,0 +1,2 @@
+name;type;bus
+electricity-excess;excess;bus-electricity

src/oemof/tabular/examples/datapackages/private_transport/data/elements/load.csv

@@ -0,0 +1,2 @@
+name;amount;profile;type;bus
+electricity-demand;5000;electricity-load-profile;load;bus-electricity

src/oemof/tabular/examples/datapackages/private_transport/data/elements/storage.csv

@@ -0,0 +1,2 @@
+name,carrier,tech,storage_capacity,capacity,capacity_cost,storage_capacity_initial,type,bus
+el-storage,lithium,battery,100,10,10,0.5,storage,bus-electricity

src/oemof/tabular/examples/datapackages/private_transport/data/elements/volatile.csv

@@ -0,0 +1,2 @@
+name;type;carrier;tech;capacity;capacity_cost;bus;marginal_cost;profile;output_parameters
+wind;volatile;wind;onshore;50;;bus-electricity;0;wind-profile;{}

src/oemof/tabular/examples/datapackages/private_transport/data/sequences/load_profile.csv

@@ -0,0 +1,4 @@
+timeindex,electricity-load-profile
+2011-01-01T00:00:00Z,0.000745659236
+2011-01-01T01:00:00Z,0.000709651546
+2011-01-01T02:00:00Z,0.00068564642

src/oemof/tabular/examples/datapackages/private_transport/data/sequences/volatile_profile.csv

@@ -0,0 +1,4 @@
+timeindex,wind-profile
+2011-01-01T00:00:00Z,0.147532
+2011-01-01T01:00:00Z,0.184181
+2011-01-01T02:00:00Z,0.223937