costs_workflow/costs/life_cycle_costs.py

"""
Life cycle costs module
"""

from hub.city_model_structure.building import Building
from configuration import Configuration


class LifeCycleCosts:
  """
  Life cycle costs class
  """
  def __init__(self, building: Building, configuration: Configuration):
    self._building = building
    self._configuration = configuration
    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
    for archetype in self._configuration.cost_catalog.entries('archetypes').archetype:
      if str(building.function) == str(archetype.function):
        self._archetype = archetype
        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

    total_floor_area = self._building.floor_area

    for internal_zone in 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

    chapters = archetype.capital_cost

    peak_heating = building.heating_peak_load[cte.YEAR].values[0]/1000
    peak_cooling = building.cooling_peak_load[cte.YEAR].values[0]/1000
    # todo: change area pv when the variable exists
    roof_area = 0
    for roof in building.roofs:
      roof_area += roof.solid_polygon.area
    surface_pv = roof_area * 0.5

    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.loc[0, 'Subsidies construction'] = 0
    self._yearly_capital_incomes.loc[0, 'Subsidies HVAC'] = 0
    self._yearly_capital_incomes.loc[0, 'Subsidies PV'] = 0

    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