p_monsalvete/costs_workflow
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases 1 Wiki Activity

Merge remote-tracking branch 'origin/main'

Browse Source 
# Conflicts:
#	life_cycle_costs.py
#	main.py
This commit is contained in:
Oriol Gavalda 2023-04-25 12:02:05 -04:00
commit e9df29a323
2 changed files with 52 additions and 64 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

85
life_cycle_costs.py
View File

@ -4,22 +4,21 @@ SPDX - License - Identifier: LGPL - 3.0 - or -later
Copyright © 2022 Project Author Pilar Monsalvete Alvarez de Uribarri pilar_monsalvete@concordia.ca Copyright © 2022 Project Author Pilar Monsalvete Alvarez de Uribarri pilar_monsalvete@concordia.ca
Project contributor 2023 Author Oriol Gavaldà Torrellas oriol.gavalda@concordia.ca Project contributor 2023 Author Oriol Gavaldà Torrellas oriol.gavalda@concordia.ca
""" """
import math import sys
import hub.catalog_factories.costs_catalog_factory as CostCatalogFactory from hub.hub_logger import logger
from hub.catalog_factories.cost.montreal_custom_catalog import MontrealCustomCatalog
from hub.helpers.dictionaries import Dictionaries from hub.helpers.dictionaries import Dictionaries
class LifeCycleCosts: class LifeCycleCosts:
# todo: this should be (city, costs_catalog) or similar # todo: this should be (city, costs_catalog) or similar
def __init__(self, city, number_of_years, consumer_price_index, discount_rate, def __init__(self, building, catalog, number_of_years, consumer_price_index, discount_rate,
retrofitting_scenario): retrofitting_scenario):
self._city = city self._building = building
self._number_of_years = number_of_years self._number_of_years = number_of_years
self._consumer_price_index = consumer_price_index self._consumer_price_index = consumer_price_index
self._discount_rate = discount_rate self._discount_rate = discount_rate
self._cost_catalog=CostCatalogFactory('montreal_custom').catalog self._cost_catalog = catalog
self._end_of_life_cost = 0 self._end_of_life_cost = 0
self._capital_costs_at_year_0 = 0 self._capital_costs_at_year_0 = 0
self._items = 0 self._items = 0
@ -28,35 +27,43 @@ class LifeCycleCosts:
self._retrofitting_scenario = retrofitting_scenario self._retrofitting_scenario = retrofitting_scenario
def calculate_capital_costs(self): def calculate_capital_costs(self):
surface_opaque=0 building = self._building
surface_transparent=0 surface_opaque = 0
city = self._city surface_transparent = 0
for building in city.buildings: try:
try: function = Dictionaries().hub_function_to_montreal_custom_costs_function[building.function]
function = 'residential' #Dictionaries().hub_function_to_nrel_construction_function[building.function] archetype = self._search_archetype(self._cost_catalog, function)
archetype = self._search_archetype(self._cost_catalog, function) except KeyError:
except KeyError: logger.error(f'Building {building.name} has unknown costs archetype for building function: '
logger.error(f'Building {building.name} has unknown construction archetype for building function: ' f'{building.function}\n')
f'{building.function} and building year of construction: {building.year_of_construction} ' sys.stderr.write(f'Building {building.name} has unknown costs archetype for building function: '
f'and climate zone reference norm {self._climate_zone}\n') f'{building.function}\n')
sys.stderr.write(f'Building {building.name} has unknown construction archetype for building function: ' return
f'{building.function} and building year of construction: {building.year_of_construction} '
f'and climate zone reference norm {self._climate_zone}\n') for internal_zone in building.internal_zones:
# if building has no thermal zones defined from geometry, and the building will be divided in storeys, for thermal_zone in internal_zone.thermal_zones:
# one thermal zone per storey is assigned for thermal_boundary in thermal_zone.thermal_boundaries:
for internal_zone in building.internal_zones: surface_opaque += thermal_boundary.opaque_area * (1-thermal_boundary.window_ratio)
if internal_zone.thermal_zones is None: surface_transparent += thermal_boundary.opaque_area * thermal_boundary.window_ratio
continue
for thermal_zone in internal_zone.thermal_zones: chapters = archetype.capital_cost.general_chapters
for thermal_boundary in thermal_zone.thermal_boundaries:
surface_opaque += thermal_boundary.opaque_area*(1-thermal_boundary.window_ratio) if self._retrofitting_scenario == 1 or self._retrofitting_scenario == 3:
surface_transparent += thermal_boundary.opaque_area*(thermal_boundary.window_ratio) chapter = chapters.chapter('B_shell')
print(f'surface_opaque {surface_opaque}') capital_cost_opaque = surface_opaque * chapter.item('B2010_opaque_walls').refurbishment[0]
print(f'surface_transparent {surface_transparent}') capital_cost_transparent = surface_transparent * chapter.item('B2020_transparent').refurbishment[0]
capital_cost_op = surface_opaque*self._cost_catalog.archetypes.archetype[0].capital_cost.general_chapters[0].items[1].refurbishment[0]
capital_cost_transp= surface_transparent*self._cost_catalog.archetypes.archetype[0].capital_cost.general_chapters[0].items[2].refurbishment[0] print(f'cost_opaque {capital_cost_opaque}')
print(f'cost_opaque {capital_cost_op}') print(f'cost_transparent {capital_cost_transparent}')
print(f'cost_transparent {capital_cost_transp}')
@staticmethod
def _search_archetype(costs_catalog, function):
costs_archetypes = costs_catalog.entries('archetypes').archetypes
for building_archetype in costs_archetypes:
if str(function) == str(building_archetype.function):
return building_archetype
raise KeyError('archetype not found')
'''total_capital_costs = self._capital_costs_at_year_0 '''total_capital_costs = self._capital_costs_at_year_0
#for year in range(1, self._number_of_years + 1): #for year in range(1, self._number_of_years + 1):
@ -86,11 +93,3 @@ class LifeCycleCosts:
for concept in self._concepts: for concept in self._concepts:
total_maintenance_costs += concept.mantainance_costs * costs_increase total_maintenance_costs += concept.mantainance_costs * costs_increase
return total_maintenance_costs''' return total_maintenance_costs'''
@staticmethod
def _search_archetype(costs_catalog, function):
costs_archetypes = costs_catalog.entries('archetypes')
for building_archetype in costs_archetypes:
if (str(function) == str(building_archetype.function)):
return building_archetype
raise KeyError('archetype not found')

31
main.py
View File

@ -6,20 +6,13 @@ Copyright © 2022 Project Author Pilar Monsalvete Álvarez de Uribarri pilar.mon
import glob import glob
import os import os
import sys
from pathlib import Path
import csv
from hub.imports.construction_factory import ConstructionFactory from hub.imports.construction_factory import ConstructionFactory
from hub.imports.usage_factory import UsageFactory
from hub.exports.exports_factory import ExportsFactory
from hub.exports.energy_building_exports_factory import EnergyBuildingsExportsFactory
from hub.helpers.dictionaries import Dictionaries from hub.helpers.dictionaries import Dictionaries
from pathlib import Path from pathlib import Path
from hub.imports.geometry_factory import GeometryFactory from hub.imports.geometry_factory import GeometryFactory
from life_cycle_costs import LifeCycleCosts from life_cycle_costs import LifeCycleCosts
import hub.catalog_factories.costs_catalog_factory as CostCatalogFactory from hub.catalog_factories.costs_catalog_factory import CostCatalogFactory
from hub.catalog_factories.cost.montreal_custom_catalog import MontrealCustomCatalog
file_path = (Path(__file__).parent.parent / 'costs_workflow' / 'input_files' / 'Citylayers_neighbours.geojson') file_path = (Path(__file__).parent.parent / 'costs_workflow' / 'input_files' / 'Citylayers_neighbours.geojson')
out_path = (Path(__file__).parent.parent / 'costs_workflow' / 'out_files') out_path = (Path(__file__).parent.parent / 'costs_workflow' / 'out_files')
@ -41,19 +34,15 @@ ConstructionFactory('nrcan', city).enrich()
print('enrich constructions... done') print('enrich constructions... done')
number_of_years = 40 number_of_years = 40
consumer_price_index = 0.1 consumer_price_index = 0.1
discount_rate=0.06 discount_rate = 0.06
retrofitting_scenario = 1 retrofitting_scenario = 1
catalog = CostCatalogFactory('montreal_custom').catalog
lcc = LifeCycleCosts(city, 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, retrofitting_scenario)
# 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() # end_of_life_costs = lcc.calculate_end_of_life_costs()
# total_maintenance_costs = lcc.calculate_total_maintenance_costs() # total_operational_costs = lcc.calculate_total_operational_costs()
# life_cycle_costs = total_capital_costs + end_of_life_costs + total_operational_costs + total_maintenance_costs # total_maintenance_costs = lcc.calculate_total_maintenance_costs()
# life_cycle_costs = total_capital_costs + end_of_life_costs + total_operational_costs + total_maintenance_costs
print(f'Building name: {building.name}')
print(f'Capital costs: {total_capital_costs}')
# print(f'End of life costs: {end_of_life_costs}')
# print(f'Operational costs: {total_operational_costs}')
# print(f'Maintenance costs: {total_maintenance_costs}')
# print(f'Life cycle costs: {life_cycle_costs}')