diff --git a/life_cycle_costs.py b/life_cycle_costs.py index cd259f4..5887367 100644 --- a/life_cycle_costs.py +++ b/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}') - print(f'cost_transparent {capital_cost_transparent}') + capital_cost_subtotal = capital_cost_skin + capital_cost_services + 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: - total_operational_costs += fuel.operational_cost \ - * math.pow(1 + fuel.energy_price_index, year) / math.pow(1 + self._discount_rate, year) + peak_cost += operational_cost_year_0 \ + * math.pow(1 + self._consumer_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: - total_maintenance_costs += concept.mantainance_costs * costs_increase - return total_maintenance_costs''' + maintenance_pv +=maintenance_pv_0*costs_increase + maintenance_heating +=maintenance_heating_0*costs_increase + maintenance_cooling +=maintenance_cooling_0*costs_increase + total_maintenance_costs = maintenance_pv + maintenance_heating + maintenance_cooling + return total_maintenance_costs diff --git a/main.py b/main.py index b0a94d9..3835539 100644 --- a/main.py +++ b/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 -consumer_price_index = 0.1 +number_of_years = 30 +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: - lcc = LifeCycleCosts(building, catalog, number_of_years, consumer_price_index, discount_rate, retrofitting_scenario) - total_capital_costs = lcc.calculate_capital_costs() - # end_of_life_costs = lcc.calculate_end_of_life_costs() - # total_operational_costs = lcc.calculate_total_operational_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 - +for i in range(1,4) : + for building in city.buildings: + lcc = LifeCycleCosts(building, catalog, number_of_years, consumer_price_index, discount_rate, i) + total_capital_costs = lcc.calculate_capital_costs() + print(f'total capital costs scenario {i} are {total_capital_costs}') + end_of_life_costs = lcc.calculate_end_of_life_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}')