Main of costs_workflow (with hardcoded values for demands and power) operating. Lifecyclecosts first version

life_cycle_costs.py

@@ -5,8 +5,10 @@ Copyright © 2022 Project Author Pilar Monsalvete Alvarez de Uribarri pilar_mons
 Project contributor 2023 Author Oriol Gavaldà Torrellas oriol.gavalda@concordia.ca
 """
 import sys
+import math
 from hub.hub_logger import logger
 from hub.helpers.dictionaries import Dictionaries
+from hub.helpers import constants as cte
 
 
 class LifeCycleCosts:
@@ -30,6 +32,13 @@ class LifeCycleCosts:
     building = self._building
     surface_opaque = 0
     surface_transparent = 0
+    surface_roof = 0
+    surface_ground = 0
+    factor_pv=0.5
+    factor_heating_power=0.1 #kW/m2
+    factor_cooling_power=0.1 #kW/m2
+    total_floor_area=0
+
     try:
       function = Dictionaries().hub_function_to_montreal_custom_costs_function[building.function]
       archetype = self._search_archetype(self._cost_catalog, function)
@@ -42,20 +51,43 @@ class LifeCycleCosts:
 
     for internal_zone in building.internal_zones:
         for thermal_zone in internal_zone.thermal_zones:
+          total_floor_area+=thermal_zone.total_floor_area
+          print(total_floor_area)
           for thermal_boundary in thermal_zone.thermal_boundaries:
-            surface_opaque += thermal_boundary.opaque_area * (1-thermal_boundary.window_ratio)
+            if thermal_boundary.type == 'Ground':
+              surface_ground += thermal_boundary.opaque_area
+            elif thermal_boundary.type == 'Roof':
+              surface_roof += thermal_boundary.opaque_area
+            elif thermal_boundary.type == 'Wall':
+              surface_opaque += thermal_boundary.opaque_area * (1-thermal_boundary.window_ratio)
             surface_transparent += thermal_boundary.opaque_area * thermal_boundary.window_ratio
+    print(f'total floor area {total_floor_area}')
 
-    chapters = archetype.capital_cost.general_chapters
+    chapters = archetype.capital_cost
-
+    capital_cost_skin=0
+    capital_cost_services=0
     if self._retrofitting_scenario == 1 or self._retrofitting_scenario == 3:
       chapter = chapters.chapter('B_shell')
       capital_cost_opaque = surface_opaque * chapter.item('B2010_opaque_walls').refurbishment[0]
       capital_cost_transparent = surface_transparent * chapter.item('B2020_transparent').refurbishment[0]
+      capital_cost_roof = surface_roof * chapter.item('B3010_opaque_roof').refurbishment[0]
+      capital_cost_ground = surface_ground * chapter.item('B10_superstructure').refurbishment[0]
+      capital_cost_skin= capital_cost_opaque+capital_cost_transparent+capital_cost_roof+capital_cost_ground
+      print (f'capital cost skin {capital_cost_skin}')
+    if self._retrofitting_scenario == 2 or self._retrofitting_scenario == 3:
+      chapter = chapters.chapter('D_services')
+      capital_cost_pv = surface_roof * factor_pv * chapter.item('D301010_photovoltaic_system').initial_investment[0]
+      capital_cost_heating_equipment = total_floor_area * factor_heating_power * chapter.item('D3020_heat_generating_systems').initial_investment[0]
+      capital_cost_cooling_equipment = total_floor_area * factor_cooling_power * chapter.item('D3030_cooling_generation_systems').initial_investment[0]
+      capital_cost_distribution_equipment = total_floor_area * factor_cooling_power * chapter.item('D3040_distribution_systems').initial_investment[0]
+      capital_cost_other_hvac_ahu = total_floor_area * factor_cooling_power * chapter.item('D3080_other_hvac_ahu').initial_investment[0]
+      capital_cost_lighting = total_floor_area * factor_pv * chapter.item('D5020_lighting_and_branch_wiring').initial_investment[0]
+      capital_cost_services=capital_cost_pv+capital_cost_heating_equipment+capital_cost_cooling_equipment+capital_cost_distribution_equipment+capital_cost_other_hvac_ahu+capital_cost_lighting
 
-      print(f'cost_opaque {capital_cost_opaque}')
+    capital_cost_subtotal = capital_cost_skin + capital_cost_services
-      print(f'cost_transparent {capital_cost_transparent}')
+    capital_cost_total= capital_cost_subtotal*(1+chapters.design_allowance)*(1+chapters.overhead_and_profit)
 
+    return capital_cost_total
   @staticmethod
   def _search_archetype(costs_catalog, function):
     costs_archetypes = costs_catalog.entries('archetypes').archetypes
@@ -64,32 +96,137 @@ class LifeCycleCosts:
         return building_archetype
     raise KeyError('archetype not found')
 
-
-'''total_capital_costs = self._capital_costs_at_year_0
-    #for year in range(1, self._number_of_years + 1):
-    #  costs_increase = math.pow(1 + self._consumer_price_index, year) / math.pow(1 + self._discount_rate, year)
-    #  for item in self._items:
-    #    total_capital_costs += item.reposition_costs[year] * costs_increase
-    #return total_capital_costs
-
   def calculate_end_of_life_costs(self):
+    building = self._building
+    total_floor_area = 0
+
+    try:
+      function = Dictionaries().hub_function_to_montreal_custom_costs_function[building.function]
+      archetype = self._search_archetype(self._cost_catalog, function)
+    except KeyError:
+      logger.error(f'Building {building.name} has unknown costs archetype for building function: '
+                   f'{building.function}\n')
+      sys.stderr.write(f'Building {building.name} has unknown costs archetype for building function: '
+                       f'{building.function}\n')
+      return
+
+    for internal_zone in building.internal_zones:
+      for thermal_zone in internal_zone.thermal_zones:
+        total_floor_area += thermal_zone.total_floor_area
+        print(total_floor_area)
+
     price_increase = 0
     for year in range(1, self._number_of_years + 1):
       price_increase += math.pow(1 + self._consumer_price_index, year) / math.pow(1 + self._discount_rate, year)
-    return self._end_of_life_cost * price_increase
+    price_increase_average = price_increase/self._number_of_years
+
+    return total_floor_area*archetype.end_of_life_cost*price_increase_average
+
 
   def calculate_total_operational_costs(self):
     total_operational_costs = 0
+    peak_cost=0
+    monthly_cost=0
+    variable_cost=0
+
+    building = self._building
+    total_floor_area = 0
+
+    try:
+      function = Dictionaries().hub_function_to_montreal_custom_costs_function[building.function]
+      archetype = self._search_archetype(self._cost_catalog, function)
+    except KeyError:
+      logger.error(f'Building {building.name} has unknown costs archetype for building function: '
+                   f'{building.function}\n')
+      sys.stderr.write(f'Building {building.name} has unknown costs archetype for building function: '
+                       f'{building.function}\n')
+      return
+
+    for internal_zone in building.internal_zones:
+      for thermal_zone in internal_zone.thermal_zones:
+        total_floor_area += thermal_zone.total_floor_area
+
+    if self._retrofitting_scenario==1 or self._retrofitting_scenario==3:
+      specific_heating_demand=50
+    else:
+      specific_heating_demand=190
+
+    heating_demand= specific_heating_demand * total_floor_area
+    cooling_demand= 10 * total_floor_area
+
+    if self._retrofitting_scenario==2 or self._retrofitting_scenario==3:
+      heating_SCOP=1
+      cooling_SEER=2.8
+    else:
+      heating_SCOP = 2.5
+      cooling_SEER = 4
+
+    electricity_heating=heating_demand/heating_SCOP
+    electricity_cooling=cooling_demand/cooling_SEER
+    electricity_lighting=11* total_floor_area
+    electricity_plug_loads=19*total_floor_area
+    domestic_hot_water_demand= 50* total_floor_area
+
+    total_electricity_consumption= electricity_cooling+electricity_heating+electricity_lighting+domestic_hot_water_demand+electricity_plug_loads
+
+    peak_electricity_demand= 0.1*total_floor_area
+
+    operational_cost_year_0=total_electricity_consumption*archetype.operational_cost.fuels[0].variable[0]
+    peak_cost_year_0=peak_electricity_demand*archetype.operational_cost.fuels[0].fixed_power*12
+    monthly_cost_year_0=archetype.operational_cost.fuels[0].fixed_monthly*12*(total_floor_area/100)
+    print(f'operational_cost_year_0 {operational_cost_year_0}')
+    print(f'peak_cost_year_0 {peak_cost_year_0}')
+    print(f'monthly_cost_year_0 {monthly_cost_year_0}')
+
     for year in range(1, self._number_of_years + 1):
-      for fuel in self._fuels:
+      peak_cost += operational_cost_year_0 \
-        total_operational_costs += fuel.operational_cost \
+                                   * math.pow(1 + self._consumer_price_index, year) / math.pow(1 + self._discount_rate, year)
-                                   * math.pow(1 + fuel.energy_price_index, year) / math.pow(1 + self._discount_rate, year)
+      monthly_cost +=peak_cost_year_0 \
+                                   * math.pow(1 + self._consumer_price_index, year) / math.pow(1 + self._discount_rate, year)
+      variable_cost +=monthly_cost_year_0 \
+                                   * math.pow(1 + self._consumer_price_index, year) / math.pow(1 + self._discount_rate, year)
+      total_operational_costs = peak_cost + monthly_cost + variable_cost
+
     return total_operational_costs
 
+
   def calculate_total_maintenance_costs(self):
-    total_maintenance_costs = 0
+
+    building = self._building
+    total_floor_area = 0
+
+    factor_pv = 0.5
+    factor_heating_power = 0.1  # kW/m2
+    factor_cooling_power = 0.1  # kW/m2
+    surface_roof = 0
+    maintenance_pv = 0
+    maintenance_heating = 0
+    maintenance_cooling = 0
+
+    try:
+      function = Dictionaries().hub_function_to_montreal_custom_costs_function[building.function]
+      archetype = self._search_archetype(self._cost_catalog, function)
+    except KeyError:
+      logger.error(f'Building {building.name} has unknown costs archetype for building function: '
+                   f'{building.function}\n')
+      sys.stderr.write(f'Building {building.name} has unknown costs archetype for building function: '
+                       f'{building.function}\n')
+      return
+
+    for internal_zone in building.internal_zones:
+        for thermal_zone in internal_zone.thermal_zones:
+          total_floor_area+=thermal_zone.total_floor_area
+          for thermal_boundary in thermal_zone.thermal_boundaries:
+            if thermal_boundary.type == 'Roof':
+              surface_roof += thermal_boundary.opaque_area
+    surface_pv = surface_roof*factor_pv
+    maintenance_pv_0=surface_pv*archetype.operational_cost.maintenance_pv
+    maintenance_heating_0=total_floor_area*factor_heating_power*archetype.operational_cost.maintenance_heating
+    maintenance_cooling_0=total_floor_area*factor_cooling_power*archetype.operational_cost.maintenance_cooling
     for year in range(1, self._number_of_years + 1):
       costs_increase = math.pow(1 + self._consumer_price_index, year) / math.pow(1 + self._discount_rate, year)
-      for concept in self._concepts:
+      maintenance_pv +=maintenance_pv_0*costs_increase
-        total_maintenance_costs += concept.mantainance_costs * costs_increase
+      maintenance_heating +=maintenance_heating_0*costs_increase
-    return total_maintenance_costs'''
+      maintenance_cooling +=maintenance_cooling_0*costs_increase
+    total_maintenance_costs = maintenance_pv + maintenance_heating + maintenance_cooling
+    return total_maintenance_costs

main.py

@@ -14,7 +14,7 @@ from hub.imports.geometry_factory import GeometryFactory
 from life_cycle_costs import LifeCycleCosts
 from hub.catalog_factories.costs_catalog_factory import CostCatalogFactory
 
-file_path = (Path(__file__).parent.parent / 'costs_workflow' / 'input_files' / 'Citylayers_neighbours.geojson')
+file_path = (Path(__file__).parent.parent / 'costs_workflow' / 'input_files' / 'selected_building_2864.geojson')
 out_path = (Path(__file__).parent.parent / 'costs_workflow' / 'out_files')
 files = glob.glob(f'{out_path}/*')
 for file in files:
@@ -32,17 +32,22 @@ city = GeometryFactory('geojson',
 print(f'city created from {file_path}')
 ConstructionFactory('nrcan', city).enrich()
 print('enrich constructions... done')
-number_of_years = 40
+number_of_years = 30
-consumer_price_index = 0.1
+consumer_price_index = 0.04
 discount_rate = 0.06
-retrofitting_scenario = 1
+#retrofitting_scenario = 2
 catalog = CostCatalogFactory('montreal_custom').catalog
 
-for building in city.buildings:
+for i in range(1,4) :
-  lcc = LifeCycleCosts(building, catalog, number_of_years, consumer_price_index, discount_rate, retrofitting_scenario)
+    for building in city.buildings:
-  total_capital_costs = lcc.calculate_capital_costs()
+      lcc = LifeCycleCosts(building, catalog, number_of_years, consumer_price_index, discount_rate, i)
-  # end_of_life_costs = lcc.calculate_end_of_life_costs()
+      total_capital_costs = lcc.calculate_capital_costs()
-  # total_operational_costs = lcc.calculate_total_operational_costs()
+      print(f'total capital costs scenario {i} are {total_capital_costs}')
-  # total_maintenance_costs = lcc.calculate_total_maintenance_costs()
+      end_of_life_costs = lcc.calculate_end_of_life_costs()
-  # life_cycle_costs = total_capital_costs + end_of_life_costs + total_operational_costs + total_maintenance_costs
+      print(f'end_of_life_costs scenario {i} are {end_of_life_costs}')
-
+      total_operational_costs = lcc.calculate_total_operational_costs()
+      print(f'total_operational_costs scenario {i} are {total_operational_costs}')
+      total_maintenance_costs = lcc.calculate_total_maintenance_costs()
+      print(f'total_maintenance_costs scenario {i} are {total_maintenance_costs}')
+      life_cycle_costs = total_capital_costs + end_of_life_costs + total_operational_costs + total_maintenance_costs
+      print(f'life_cycle_costs scenario {i} are {life_cycle_costs}')