365 lines
18 KiB
Python
365 lines
18 KiB
Python
|
"""
|
||
|
|
Capital costs module
|
||
|
|
SPDX - License - Identifier: LGPL - 3.0 - or -later
|
||
|
|
Copyright © 2023 Project Coder Guille Gutierrez guillermo.gutierrezmorote@concordia.ca
|
||
|
|
Code contributor Pilar Monsalvete Alvarez de Uribarri pilar.monsalvete@concordia.ca
|
||
|
|
Code contributor Oriol Gavalda Torrellas oriol.gavalda@concordia.ca
|
||
|
|
"""
|
||
|
|
import math
|
||
|
|
|
||
|
|
import pandas as pd
|
||
|
|
import numpy_financial as npf
|
||
|
|
from hub.city_model_structure.building import Building
|
||
|
|
import hub.helpers.constants as cte
|
||
|
|
from scripts.costs.configuration import Configuration
|
||
|
|
from scripts.costs.constants import SKIN_RETROFIT, SKIN_RETROFIT_AND_SYSTEM_RETROFIT_AND_PV, SYSTEM_RETROFIT_AND_PV
|
||
|
|
from scripts.costs.cost_base import CostBase
|
||
|
|
|
||
|
|
|
||
|
|
class CapitalCosts(CostBase):
|
||
|
|
"""
|
||
|
|
Capital costs class
|
||
|
|
"""
|
||
|
|
|
||
|
|
def __init__(self, building: Building, configuration: Configuration):
|
||
|
|
super().__init__(building, configuration)
|
||
|
|
self._yearly_capital_costs = pd.DataFrame(
|
||
|
|
index=self._rng,
|
||
|
|
columns=[
|
||
|
|
'B2010_opaque_walls',
|
||
|
|
'B2020_transparent',
|
||
|
|
'B3010_opaque_roof',
|
||
|
|
'B1010_superstructure',
|
||
|
|
'D2010_photovoltaic_system',
|
||
|
|
'D3020_heat_and_cooling_generating_systems',
|
||
|
|
'D3040_distribution_systems',
|
||
|
|
'D3050_other_hvac_ahu',
|
||
|
|
'D3060_storage_systems',
|
||
|
|
'D40_dhw',
|
||
|
|
'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, 'B1010_superstructure'] = 0
|
||
|
|
self._yearly_capital_costs.loc[0, 'D2010_photovoltaic_system'] = 0
|
||
|
|
self._yearly_capital_costs.loc[0, 'D3020_heat_and_cooling_generating_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, 'D3060_storage_systems'] = 0
|
||
|
|
self._yearly_capital_costs.loc[0, 'D40_dhw'] = 0
|
||
|
|
# self._yearly_capital_costs.loc[0, 'D5020_lighting_and_branch_wiring'] = 0
|
||
|
|
|
||
|
|
self._yearly_capital_incomes = pd.DataFrame(
|
||
|
|
index=self._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
|
||
|
|
self._yearly_capital_costs.fillna(0, inplace=True)
|
||
|
|
self._own_capital = 1 - self._configuration.percentage_credit
|
||
|
|
self._surface_pv = 0
|
||
|
|
for roof in self._building.roofs:
|
||
|
|
self._surface_pv += roof.solid_polygon.area * roof.solar_collectors_area_reduction_factor
|
||
|
|
|
||
|
|
def calculate(self) -> tuple[pd.DataFrame, pd.DataFrame]:
|
||
|
|
if self._configuration.retrofit_scenario in (SKIN_RETROFIT, SKIN_RETROFIT_AND_SYSTEM_RETROFIT_AND_PV):
|
||
|
|
self.skin_capital_cost()
|
||
|
|
if self._configuration.retrofit_scenario in (SYSTEM_RETROFIT_AND_PV, SKIN_RETROFIT_AND_SYSTEM_RETROFIT_AND_PV):
|
||
|
|
self.energy_system_capital_cost()
|
||
|
|
|
||
|
|
self.skin_yearly_capital_costs()
|
||
|
|
self.yearly_energy_system_costs()
|
||
|
|
self.yearly_incomes()
|
||
|
|
return self._yearly_capital_costs, self._yearly_capital_incomes
|
||
|
|
|
||
|
|
def skin_capital_cost(self):
|
||
|
|
"""
|
||
|
|
calculating skin costs
|
||
|
|
:return:
|
||
|
|
"""
|
||
|
|
surface_opaque = 0
|
||
|
|
surface_transparent = 0
|
||
|
|
surface_roof = 0
|
||
|
|
surface_ground = 0
|
||
|
|
|
||
|
|
for thermal_zone in self._building.thermal_zones_from_internal_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
|
||
|
|
|
||
|
|
chapter = self._capital_costs_chapter.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('B1010_superstructure').refurbishment[0]
|
||
|
|
self._yearly_capital_costs.loc[0, 'B2010_opaque_walls'] = capital_cost_opaque * self._own_capital
|
||
|
|
self._yearly_capital_costs.loc[0, 'B2020_transparent'] = capital_cost_transparent * self._own_capital
|
||
|
|
self._yearly_capital_costs.loc[0, 'B3010_opaque_roof'] = capital_cost_roof * self._own_capital
|
||
|
|
self._yearly_capital_costs.loc[0, 'B1010_superstructure'] = capital_cost_ground * self._own_capital
|
||
|
|
capital_cost_skin = capital_cost_opaque + capital_cost_ground + capital_cost_transparent + capital_cost_roof
|
||
|
|
return capital_cost_opaque, capital_cost_transparent, capital_cost_roof, capital_cost_ground, capital_cost_skin
|
||
|
|
|
||
|
|
def skin_yearly_capital_costs(self):
|
||
|
|
skin_capital_cost = self.skin_capital_cost()
|
||
|
|
for year in range(1, self._configuration.number_of_years):
|
||
|
|
self._yearly_capital_costs.loc[year, 'B2010_opaque_walls'] = (
|
||
|
|
-npf.pmt(
|
||
|
|
self._configuration.interest_rate,
|
||
|
|
self._configuration.credit_years,
|
||
|
|
skin_capital_cost[0] * self._configuration.percentage_credit
|
||
|
|
)
|
||
|
|
)
|
||
|
|
self._yearly_capital_costs.loc[year, 'B2020_transparent'] = (
|
||
|
|
-npf.pmt(
|
||
|
|
self._configuration.interest_rate,
|
||
|
|
self._configuration.credit_years,
|
||
|
|
skin_capital_cost[1] * self._configuration.percentage_credit
|
||
|
|
)
|
||
|
|
)
|
||
|
|
self._yearly_capital_costs.loc[year, 'B3010_opaque_roof'] = (
|
||
|
|
-npf.pmt(
|
||
|
|
self._configuration.interest_rate,
|
||
|
|
self._configuration.credit_years,
|
||
|
|
skin_capital_cost[2] * self._configuration.percentage_credit
|
||
|
|
)
|
||
|
|
)
|
||
|
|
self._yearly_capital_costs.loc[year, 'B1010_superstructure'] = (
|
||
|
|
-npf.pmt(
|
||
|
|
self._configuration.interest_rate,
|
||
|
|
self._configuration.credit_years,
|
||
|
|
skin_capital_cost[3] * self._configuration.percentage_credit
|
||
|
|
)
|
||
|
|
)
|
||
|
|
|
||
|
|
def energy_system_capital_cost(self):
|
||
|
|
chapter = self._capital_costs_chapter.chapter('D_services')
|
||
|
|
energy_system_components = self.system_components()
|
||
|
|
system_components = energy_system_components[0]
|
||
|
|
component_categories = energy_system_components[1]
|
||
|
|
component_sizes = energy_system_components[-1]
|
||
|
|
capital_cost_heating_and_cooling_equipment = 0
|
||
|
|
capital_cost_domestic_hot_water_equipment = 0
|
||
|
|
capital_cost_energy_storage_equipment = 0
|
||
|
|
capital_cost_distribution_equipment = 0
|
||
|
|
capital_cost_lighting = 0
|
||
|
|
capital_cost_pv = self._surface_pv * chapter.item('D2010_photovoltaic_system').initial_investment[0]
|
||
|
|
# capital_cost_lighting = self._total_floor_area * \
|
||
|
|
# chapter.item('D5020_lighting_and_branch_wiring').initial_investment[0]
|
||
|
|
for (i, component) in enumerate(system_components):
|
||
|
|
if component_categories[i] == 'generation':
|
||
|
|
capital_cost_heating_and_cooling_equipment += chapter.item(component).initial_investment[0] * component_sizes[i]
|
||
|
|
elif component_categories[i] == 'dhw':
|
||
|
|
capital_cost_domestic_hot_water_equipment += chapter.item(component).initial_investment[0] * \
|
||
|
|
component_sizes[i]
|
||
|
|
elif component_categories[i] == 'distribution':
|
||
|
|
capital_cost_distribution_equipment += chapter.item(component).initial_investment[0] * \
|
||
|
|
component_sizes[i]
|
||
|
|
else:
|
||
|
|
capital_cost_energy_storage_equipment += chapter.item(component).initial_investment[0] * component_sizes[i]
|
||
|
|
|
||
|
|
self._yearly_capital_costs.loc[0, 'D2010_photovoltaic_system'] = capital_cost_pv
|
||
|
|
self._yearly_capital_costs.loc[0, 'D3020_heat_and_cooling_generating_systems'] = (
|
||
|
|
capital_cost_heating_and_cooling_equipment * self._own_capital)
|
||
|
|
self._yearly_capital_costs.loc[0, 'D3040_distribution_systems'] = (
|
||
|
|
capital_cost_distribution_equipment * self._own_capital)
|
||
|
|
self._yearly_capital_costs.loc[0, 'D3060_storage_systems'] = (
|
||
|
|
capital_cost_energy_storage_equipment * self._own_capital)
|
||
|
|
self._yearly_capital_costs.loc[0, 'D40_dhw'] = (
|
||
|
|
capital_cost_domestic_hot_water_equipment * self._own_capital)
|
||
|
|
# self._yearly_capital_costs.loc[0, 'D5020_lighting_and_branch_wiring'] = capital_cost_lighting * self._own_capital
|
||
|
|
capital_cost_hvac = capital_cost_heating_and_cooling_equipment + capital_cost_distribution_equipment + capital_cost_energy_storage_equipment + capital_cost_domestic_hot_water_equipment
|
||
|
|
return (capital_cost_pv, capital_cost_heating_and_cooling_equipment, capital_cost_distribution_equipment,
|
||
|
|
capital_cost_energy_storage_equipment, capital_cost_domestic_hot_water_equipment, capital_cost_lighting, capital_cost_hvac)
|
||
|
|
|
||
|
|
def yearly_energy_system_costs(self):
|
||
|
|
chapter = self._capital_costs_chapter.chapter('D_services')
|
||
|
|
system_investment_costs = self.energy_system_capital_cost()
|
||
|
|
system_components = self.system_components()[0]
|
||
|
|
component_categories = self.system_components()[1]
|
||
|
|
component_sizes = self.system_components()[2]
|
||
|
|
for year in range(1, self._configuration.number_of_years):
|
||
|
|
costs_increase = math.pow(1 + self._configuration.consumer_price_index, year)
|
||
|
|
self._yearly_capital_costs.loc[year, 'D2010_photovoltaic_system'] = (
|
||
|
|
-npf.pmt(
|
||
|
|
self._configuration.interest_rate,
|
||
|
|
self._configuration.credit_years,
|
||
|
|
system_investment_costs[0] * self._configuration.percentage_credit
|
||
|
|
)
|
||
|
|
)
|
||
|
|
self._yearly_capital_costs.loc[year, 'D3020_heat_and_cooling_generating_systems'] = (
|
||
|
|
-npf.pmt(
|
||
|
|
self._configuration.interest_rate,
|
||
|
|
self._configuration.credit_years,
|
||
|
|
system_investment_costs[1] * self._configuration.percentage_credit
|
||
|
|
)
|
||
|
|
)
|
||
|
|
self._yearly_capital_costs.loc[year, 'D3040_distribution_systems'] = (
|
||
|
|
-npf.pmt(
|
||
|
|
self._configuration.interest_rate,
|
||
|
|
self._configuration.credit_years,
|
||
|
|
system_investment_costs[2] * self._configuration.percentage_credit
|
||
|
|
)
|
||
|
|
)
|
||
|
|
self._yearly_capital_costs.loc[year, 'D3060_storage_systems'] = (
|
||
|
|
-npf.pmt(
|
||
|
|
self._configuration.interest_rate,
|
||
|
|
self._configuration.credit_years,
|
||
|
|
system_investment_costs[3] * self._configuration.percentage_credit
|
||
|
|
)
|
||
|
|
)
|
||
|
|
self._yearly_capital_costs.loc[year, 'D40_dhw'] = (
|
||
|
|
-npf.pmt(
|
||
|
|
self._configuration.interest_rate,
|
||
|
|
self._configuration.credit_years,
|
||
|
|
system_investment_costs[4] * self._configuration.percentage_credit
|
||
|
|
)
|
||
|
|
)
|
||
|
|
# self._yearly_capital_costs.loc[year, 'D5020_lighting_and_branch_wiring'] = (
|
||
|
|
# -npf.pmt(
|
||
|
|
# self._configuration.interest_rate,
|
||
|
|
# self._configuration.credit_years,
|
||
|
|
# system_investment_costs[5] * self._configuration.percentage_credit
|
||
|
|
# )
|
||
|
|
# )
|
||
|
|
# if (year % chapter.item('D5020_lighting_and_branch_wiring').lifetime) == 0:
|
||
|
|
# reposition_cost_lighting = (
|
||
|
|
# self._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._configuration.retrofit_scenario in (SYSTEM_RETROFIT_AND_PV, SKIN_RETROFIT_AND_SYSTEM_RETROFIT_AND_PV):
|
||
|
|
if (year % chapter.item('D2010_photovoltaic_system').lifetime) == 0:
|
||
|
|
self._yearly_capital_costs.loc[year, 'D2010_photovoltaic_system'] += (
|
||
|
|
self._surface_pv * chapter.item('D2010_photovoltaic_system').reposition[0] * costs_increase
|
||
|
|
)
|
||
|
|
for (i, component) in enumerate(system_components):
|
||
|
|
if (year % chapter.item(component).lifetime) == 0 and year != (self._configuration.number_of_years - 1):
|
||
|
|
if component_categories[i] == 'generation':
|
||
|
|
reposition_cost_heating_and_cooling_equipment = chapter.item(component).reposition[0] * component_sizes[i] * costs_increase
|
||
|
|
self._yearly_capital_costs.loc[year, 'D3020_heat_and_cooling_generating_systems'] += reposition_cost_heating_and_cooling_equipment
|
||
|
|
elif component_categories[i] == 'dhw':
|
||
|
|
reposition_cost_domestic_hot_water_equipment = chapter.item(component).reposition[0] * component_sizes[i] * costs_increase
|
||
|
|
self._yearly_capital_costs.loc[year, 'D40_dhw'] += reposition_cost_domestic_hot_water_equipment
|
||
|
|
elif component_categories[i] == 'distribution':
|
||
|
|
reposition_cost_distribution_equipment = chapter.item(component).reposition[0] * component_sizes[i] * costs_increase
|
||
|
|
self._yearly_capital_costs.loc[year, 'D3040_distribution_systems'] += reposition_cost_distribution_equipment
|
||
|
|
else:
|
||
|
|
reposition_cost_energy_storage_equipment = chapter.item(component).initial_investment[0] * component_sizes[i] * costs_increase
|
||
|
|
self._yearly_capital_costs.loc[year, 'D3060_storage_systems'] += reposition_cost_energy_storage_equipment
|
||
|
|
|
||
|
|
def system_components(self):
|
||
|
|
system_components = []
|
||
|
|
component_categories = []
|
||
|
|
sizes = []
|
||
|
|
energy_systems = self._building.energy_systems
|
||
|
|
for energy_system in energy_systems:
|
||
|
|
demand_types = energy_system.demand_types
|
||
|
|
generation_systems = energy_system.generation_systems
|
||
|
|
distribution_systems = energy_system.distribution_systems
|
||
|
|
for generation_system in generation_systems:
|
||
|
|
if generation_system.system_type != cte.PHOTOVOLTAIC:
|
||
|
|
heating_capacity = generation_system.nominal_heat_output or 0
|
||
|
|
cooling_capacity = generation_system.nominal_cooling_output or 0
|
||
|
|
installed_capacity = max(heating_capacity, cooling_capacity) / 1000
|
||
|
|
if cte.DOMESTIC_HOT_WATER in demand_types and cte.HEATING not in demand_types:
|
||
|
|
component_categories.append('dhw')
|
||
|
|
sizes.append(installed_capacity)
|
||
|
|
if generation_system.system_type == cte.HEAT_PUMP:
|
||
|
|
system_components.append(self.heat_pump_type(generation_system, domestic_how_water=True))
|
||
|
|
elif generation_system.system_type == cte.BOILER and generation_system.fuel_type == cte.ELECTRICITY:
|
||
|
|
system_components.append(self.boiler_type(generation_system))
|
||
|
|
else:
|
||
|
|
system_components.append('D302010_template_heat')
|
||
|
|
elif cte.HEATING or cte.COOLING in demand_types:
|
||
|
|
component_categories.append('generation')
|
||
|
|
sizes.append(installed_capacity)
|
||
|
|
if generation_system.system_type == cte.HEAT_PUMP:
|
||
|
|
item_type = self.heat_pump_type(generation_system)
|
||
|
|
system_components.append(item_type)
|
||
|
|
elif generation_system.system_type == cte.BOILER:
|
||
|
|
item_type = self.boiler_type(generation_system)
|
||
|
|
system_components.append(item_type)
|
||
|
|
else:
|
||
|
|
if cte.COOLING in demand_types and cte.HEATING not in demand_types:
|
||
|
|
system_components.append('D302090_template_cooling')
|
||
|
|
else:
|
||
|
|
system_components.append('D302010_template_heat')
|
||
|
|
|
||
|
|
if generation_system.energy_storage_systems is not None:
|
||
|
|
energy_storage_systems = generation_system.energy_storage_systems
|
||
|
|
for storage_system in energy_storage_systems:
|
||
|
|
if storage_system.type_energy_stored == 'thermal':
|
||
|
|
component_categories.append('thermal storage')
|
||
|
|
sizes.append(storage_system.volume or 0)
|
||
|
|
system_components.append('D306010_storage_tank')
|
||
|
|
if distribution_systems is not None:
|
||
|
|
for distribution_system in distribution_systems:
|
||
|
|
component_categories.append('distribution')
|
||
|
|
sizes.append(self._building.cooling_peak_load[cte.YEAR][0] / 3.6e6)
|
||
|
|
system_components.append('D3040_distribution_systems')
|
||
|
|
return system_components, component_categories, sizes
|
||
|
|
|
||
|
|
def yearly_incomes(self):
|
||
|
|
capital_cost_skin = self.skin_capital_cost()[-1]
|
||
|
|
system_investment_cost = self.energy_system_capital_cost()
|
||
|
|
capital_cost_hvac = system_investment_cost[-1]
|
||
|
|
capital_cost_pv = system_investment_cost[0]
|
||
|
|
|
||
|
|
self._yearly_capital_incomes.loc[0, 'Subsidies construction'] = (
|
||
|
|
capital_cost_skin * self._archetype.income.construction_subsidy/100
|
||
|
|
)
|
||
|
|
self._yearly_capital_incomes.loc[0, 'Subsidies HVAC'] = capital_cost_hvac * self._archetype.income.hvac_subsidy/100
|
||
|
|
self._yearly_capital_incomes.loc[0, 'Subsidies PV'] = capital_cost_pv * self._archetype.income.photovoltaic_subsidy/100
|
||
|
|
self._yearly_capital_incomes.fillna(0, inplace=True)
|
||
|
|
|
||
|
|
@staticmethod
|
||
|
|
def heat_pump_type(generation_system, domestic_how_water=False):
|
||
|
|
source_medium = generation_system.source_medium
|
||
|
|
supply_medium = generation_system.supply_medium
|
||
|
|
if domestic_how_water:
|
||
|
|
heat_pump_item = 'D4010_hot_water_heat_pump'
|
||
|
|
else:
|
||
|
|
if source_medium == cte.AIR and supply_medium == cte.WATER:
|
||
|
|
heat_pump_item = 'D302020_air_to_water_heat_pump'
|
||
|
|
elif source_medium == cte.AIR and supply_medium == cte.AIR:
|
||
|
|
heat_pump_item = 'D302050_air_to_air_heat_pump'
|
||
|
|
elif source_medium == cte.GROUND and supply_medium == cte.WATER:
|
||
|
|
heat_pump_item = 'D302030_ground_to_water_heat_pump'
|
||
|
|
elif source_medium == cte.GROUND and supply_medium == cte.AIR:
|
||
|
|
heat_pump_item = 'D302100_ground_to_air_heat_pump'
|
||
|
|
elif source_medium == cte.WATER and supply_medium == cte.WATER:
|
||
|
|
heat_pump_item = 'D302040_water_to_water_heat_pump'
|
||
|
|
elif source_medium == cte.WATER and supply_medium == cte.AIR:
|
||
|
|
heat_pump_item = 'D302110_water_to_air_heat_pump'
|
||
|
|
else:
|
||
|
|
heat_pump_item = 'D302010_template_heat'
|
||
|
|
return heat_pump_item
|
||
|
|
|
||
|
|
@staticmethod
|
||
|
|
def boiler_type(generation_system):
|
||
|
|
fuel = generation_system.fuel_type
|
||
|
|
if fuel == cte.ELECTRICITY:
|
||
|
|
boiler_item = 'D302080_electrical_boiler'
|
||
|
|
elif fuel == cte.GAS:
|
||
|
|
boiler_item = 'D302070_natural_gas_boiler'
|
||
|
|
else:
|
||
|
|
boiler_item = 'D302010_template_heat'
|
||
|
|
return boiler_item
|
||
|
|
|
||
|
|
|
||
|
|
|
||
|
|
|