""" Life cycle costs module """ import pandas as pd from hub.city_model_structure.building import Building import hub.helpers.constants as cte from configuration import Configuration class LifeCycleCosts: """ Life cycle costs class """ def __init__(self, building: Building, configuration: Configuration, retrofit_scenario): self._building = building self._configuration = configuration self._retrofit_scenario = retrofit_scenario self._total_floor_area = 0 for internal_zone in building.internal_zones: for thermal_zone in internal_zone.thermal_zones: self._total_floor_area += thermal_zone.total_floor_area self._archetype = None self._capital_costs_chapter = None rng = range(configuration.number_of_years) self._yearly_capital_costs = pd.DataFrame( index=rng, columns=[ 'B2010_opaque_walls', 'B2020_transparent', 'B3010_opaque_roof', 'B10_superstructure', 'D301010_photovoltaic_system', 'D3020_heat_generating_systems', 'D3030_cooling_generation_systems', 'D3040_distribution_systems', 'D3080_other_hvac_ahu', 'D5020_lighting_and_branch_wiring' ], dtype='float' ) self._yearly_capital_costs.loc[0, 'B2010_opaque_walls'] = 0 self._yearly_capital_costs.loc[0]['B2020_transparent'] = 0 self._yearly_capital_costs.loc[0, 'B3010_opaque_roof'] = 0 self._yearly_capital_costs.loc[0]['B10_superstructure'] = 0 self._yearly_capital_costs.loc[0, 'D3020_heat_generating_systems'] = 0 self._yearly_capital_costs.loc[0, 'D3030_cooling_generation_systems'] = 0 self._yearly_capital_costs.loc[0, 'D3040_distribution_systems'] = 0 self._yearly_capital_costs.loc[0, 'D3080_other_hvac_ahu'] = 0 self._yearly_capital_costs.loc[0, 'D5020_lighting_and_branch_wiring'] = 0 self._yearly_capital_incomes = pd.DataFrame( index=rng, columns=[ 'Subsidies construction', 'Subsidies HVAC', 'Subsidies PV' ], dtype='float' ) self._yearly_capital_incomes.loc[0, 'Subsidies construction'] = 0 self._yearly_capital_incomes.loc[0, 'Subsidies HVAC'] = 0 self._yearly_capital_incomes.loc[0, 'Subsidies PV'] = 0 for archetype in self._configuration.cost_catalog.entries('archetypes').archetype: if str(building.function) == str(archetype.function): self._archetype = archetype self._capital_costs_chapter = self._archetype.capital_cost break if not self._archetype: raise KeyError('archetype not found') def calculate_capital_costs(self): """ Calculate capital cost :return: pd.DataFrame """ surface_opaque = 0 surface_transparent = 0 surface_roof = 0 surface_ground = 0 capital_cost_pv = 0 capital_cost_opaque = 0 capital_cost_ground = 0 capital_cost_transparent = 0 capital_cost_roof = 0 capital_cost_heating_equipment = 0 capital_cost_cooling_equipment = 0 capital_cost_distribution_equipment = 0 capital_cost_other_hvac_ahu = 0 capital_cost_lighting = 0 for internal_zone in self._building.internal_zones: for thermal_zone in internal_zone.thermal_zones: for thermal_boundary in thermal_zone.thermal_boundaries: 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 peak_heating = self._building.heating_peak_load[cte.YEAR].values[0] / 1000 peak_cooling = self._building.cooling_peak_load[cte.YEAR].values[0] / 1000 surface_pv = 0 for roof in self._building.roofs: surface_pv += roof.solid_polygon.area * roof.solar_collectors_area_reduction_factor self._yearly_capital_costs.fillna(0, inplace=True) if self._retrofitting_scenario in (SKIN_RETROFIT, SKIN_RETROFIT_AND_SYSTEM_RETROFIT_AND_PV): 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] self._yearly_capital_costs.loc[0, 'B2010_opaque_walls'] = capital_cost_opaque * (1-PERCENTAGE_CREDIT) self._yearly_capital_costs.loc[0]['B2020_transparent'] = capital_cost_transparent * (1-PERCENTAGE_CREDIT) self._yearly_capital_costs.loc[0, 'B3010_opaque_roof'] = capital_cost_roof * (1-PERCENTAGE_CREDIT) self._yearly_capital_costs.loc[0]['B10_superstructure'] = capital_cost_ground * (1-PERCENTAGE_CREDIT) if self._retrofitting_scenario in (SYSTEM_RETROFIT_AND_PV, SKIN_RETROFIT_AND_SYSTEM_RETROFIT_AND_PV): chapter = chapters.chapter('D_services') capital_cost_pv = surface_pv * chapter.item('D301010_photovoltaic_system').initial_investment[0] self._yearly_capital_costs.loc[0]['D301010_photovoltaic_system'] = capital_cost_pv capital_cost_heating_equipment = ( peak_heating * chapter.item('D3020_heat_generating_systems').initial_investment[0] ) capital_cost_cooling_equipment = ( peak_cooling * chapter.item('D3030_cooling_generation_systems').initial_investment[0] ) capital_cost_distribution_equipment = ( peak_cooling * chapter.item('D3040_distribution_systems').initial_investment[0] ) capital_cost_other_hvac_ahu = peak_cooling * chapter.item('D3080_other_hvac_ahu').initial_investment[0] capital_cost_lighting = total_floor_area * chapter.item('D5020_lighting_and_branch_wiring').initial_investment[0] self._yearly_capital_costs.loc[0, 'D3020_heat_generating_systems'] = capital_cost_heating_equipment * (1-PERCENTAGE_CREDIT) self._yearly_capital_costs.loc[0, 'D3030_cooling_generation_systems'] = capital_cost_cooling_equipment * (1-PERCENTAGE_CREDIT) self._yearly_capital_costs.loc[0, 'D3040_distribution_systems'] = capital_cost_distribution_equipment * (1-PERCENTAGE_CREDIT) self._yearly_capital_costs.loc[0, 'D3080_other_hvac_ahu'] = capital_cost_other_hvac_ahu * (1-PERCENTAGE_CREDIT) self._yearly_capital_costs.loc[0, 'D5020_lighting_and_branch_wiring'] = capital_cost_lighting * (1-PERCENTAGE_CREDIT) for year in range(1, self._number_of_years): chapter = chapters.chapter('D_services') costs_increase = math.pow(1 + self._consumer_price_index, year) self._yearly_capital_costs.loc[year, 'B2010_opaque_walls'] = -npf.pmt(INTEREST_RATE, CREDIT_YEARS, capital_cost_opaque * (PERCENTAGE_CREDIT)) self._yearly_capital_costs.loc[year, 'B2020_transparent'] = -npf.pmt(INTEREST_RATE, CREDIT_YEARS, capital_cost_transparent * (PERCENTAGE_CREDIT) ) self._yearly_capital_costs.loc[year, 'B3010_opaque_roof'] = -npf.pmt(INTEREST_RATE, CREDIT_YEARS,capital_cost_roof * (PERCENTAGE_CREDIT)) self._yearly_capital_costs.loc[year, 'B10_superstructure'] = -npf.pmt(INTEREST_RATE, CREDIT_YEARS, capital_cost_ground * (PERCENTAGE_CREDIT)) self._yearly_capital_costs.loc[year, 'D3020_heat_generating_systems'] = -npf.pmt(INTEREST_RATE,CREDIT_YEARS, capital_cost_heating_equipment * (PERCENTAGE_CREDIT)) self._yearly_capital_costs.loc[year, 'D3030_cooling_generation_systems'] = -npf.pmt(INTEREST_RATE, CREDIT_YEARS, capital_cost_cooling_equipment * (PERCENTAGE_CREDIT)) self._yearly_capital_costs.loc[year, 'D3040_distribution_systems'] = -npf.pmt(INTEREST_RATE, CREDIT_YEARS, capital_cost_distribution_equipment * (PERCENTAGE_CREDIT)) self._yearly_capital_costs.loc[year, 'D3080_other_hvac_ahu'] = -npf.pmt(INTEREST_RATE, CREDIT_YEARS, capital_cost_other_hvac_ahu * (PERCENTAGE_CREDIT)) self._yearly_capital_costs.loc[year, 'D5020_lighting_and_branch_wiring'] = -npf.pmt(INTEREST_RATE, CREDIT_YEARS, capital_cost_lighting * (PERCENTAGE_CREDIT)) if (year % chapter.item('D3020_heat_generating_systems').lifetime) == 0: reposition_cost_heating_equipment = peak_heating * chapter.item('D3020_heat_generating_systems').reposition[0] \ * costs_increase self._yearly_capital_costs.loc[year, 'D3020_heat_generating_systems'] += reposition_cost_heating_equipment if (year % chapter.item('D3030_cooling_generation_systems').lifetime) == 0: reposition_cost_cooling_equipment = peak_cooling \ * chapter.item('D3030_cooling_generation_systems').reposition[0] \ * costs_increase self._yearly_capital_costs.loc[year, 'D3030_cooling_generation_systems'] += reposition_cost_cooling_equipment if (year % chapter.item('D3080_other_hvac_ahu').lifetime) == 0: reposition_cost_hvac_ahu = peak_cooling * chapter.item('D3080_other_hvac_ahu').reposition[0] * costs_increase self._yearly_capital_costs.loc[year, 'D3080_other_hvac_ahu'] = reposition_cost_hvac_ahu if (year % chapter.item('D5020_lighting_and_branch_wiring').lifetime) == 0: reposition_cost_lighting = total_floor_area * chapter.item('D5020_lighting_and_branch_wiring').reposition[0] \ * costs_increase self._yearly_capital_costs.loc[year, 'D5020_lighting_and_branch_wiring'] += reposition_cost_lighting if self._retrofitting_scenario in (SYSTEM_RETROFIT_AND_PV, SKIN_RETROFIT_AND_SYSTEM_RETROFIT_AND_PV): if (year % chapter.item('D301010_photovoltaic_system').lifetime) == 0: self._yearly_capital_costs.loc[year]['D301010_photovoltaic_system'] += surface_pv \ * chapter.item( 'D301010_photovoltaic_system').reposition[0] * costs_increase capital_cost_skin = capital_cost_opaque + capital_cost_ground + capital_cost_transparent + capital_cost_roof capital_cost_hvac = ( capital_cost_heating_equipment + capital_cost_cooling_equipment + capital_cost_distribution_equipment + capital_cost_other_hvac_ahu + capital_cost_lighting ) self._yearly_capital_incomes.loc[0, 'Subsidies construction'] = ( capital_cost_skin * archetype.income.construction_subsidy/100 ) self._yearly_capital_incomes.loc[0, 'Subsidies HVAC'] = capital_cost_hvac * archetype.income.hvac_subsidy/100 self._yearly_capital_incomes.loc[0, 'Subsidies PV'] = capital_cost_pv * archetype.income.photovoltaic_subsidy/100 self._yearly_capital_incomes.fillna(0, inplace=True) return self._yearly_capital_costs, self._yearly_capital_incomes