Compare commits
No commits in common. "725d6d1e48e256929313e6198ec6471c39a2de3e" and "a230c6b61c7b9f2956e72a25e1bb2166fbf99d9c" have entirely different histories.
725d6d1e48
...
a230c6b61c
@ -1,209 +0,0 @@
|
|||||||
import pandas as pd
|
|
||||||
import numpy as np
|
|
||||||
import hub.helpers.constants as cte
|
|
||||||
from hub.exports.energy_building_exports_factory import EnergyBuildingsExportsFactory
|
|
||||||
from hub.imports.geometry_factory import GeometryFactory
|
|
||||||
from hub.imports.results_factory import ResultFactory
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
inflation_rate = 0.03
|
|
||||||
discount_rate = 0.05
|
|
||||||
period = 25
|
|
||||||
installation_cost = 0
|
|
||||||
tax_deduct= 0
|
|
||||||
incentive= 0
|
|
||||||
capacity =30
|
|
||||||
degradation_rate = 0.01
|
|
||||||
year_of_replacement_list= [12]
|
|
||||||
replacement_ratio = 0.1
|
|
||||||
maintenance_cost_ratio =0.01
|
|
||||||
dataframe_path= r'C:\Users\z_keshav\CMM_PV\data\test.csv'
|
|
||||||
Building_function = "residential"
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
def calculate_pv_system_metrics(
|
|
||||||
dataframe_path, # input from Hub
|
|
||||||
Building_function, # input from Hub
|
|
||||||
inflation_rate,
|
|
||||||
discount_rate,
|
|
||||||
period,
|
|
||||||
capacity, # input from Hub
|
|
||||||
degradation_rate,
|
|
||||||
year_of_replacement_list,
|
|
||||||
replacement_ratio,
|
|
||||||
maintenance_cost_ratio,
|
|
||||||
installation_cost=0,
|
|
||||||
tax_deduct=0,
|
|
||||||
incentive=0,
|
|
||||||
):
|
|
||||||
# Read the data
|
|
||||||
dataframe = pd.read_csv(dataframe_path)
|
|
||||||
building_hourly_consumption = dataframe['GRID_kWh'] # input from Hub
|
|
||||||
PV_hourly_generation = dataframe['PV_roofs_top_E_kWh'] # input from Hub
|
|
||||||
|
|
||||||
# Defining tariff based on building function
|
|
||||||
if Building_function == "residential": # Rate D when the maximum power demand has reached 50 kW or more
|
|
||||||
grid_current_tariff = 0.06704 # Residential tariff in $/kWh
|
|
||||||
elif Building_function == "commercial": # Rate G: General rate for small-power customers with demand ≤ 50 kW
|
|
||||||
grid_current_tariff = 0.11518 # Commercial tariff in $/kWh
|
|
||||||
|
|
||||||
# Initial Calculations for Year 1
|
|
||||||
first_year_generation_PV = PV_hourly_generation.sum()
|
|
||||||
first_year_self_consumption = np.minimum(PV_hourly_generation, building_hourly_consumption).sum()
|
|
||||||
first_year_grid_purchase = np.maximum(building_hourly_consumption - PV_hourly_generation, 0).sum()
|
|
||||||
first_year_PV_export = np.maximum(PV_hourly_generation - building_hourly_consumption, 0).sum()
|
|
||||||
|
|
||||||
# Cost per kW determination
|
|
||||||
if capacity <= 2.5:
|
|
||||||
cost_per_kW = 4000
|
|
||||||
elif 2.5 < capacity <= 5:
|
|
||||||
cost_per_kW = 3000
|
|
||||||
elif 5 < capacity <= 10:
|
|
||||||
cost_per_kW = 2500
|
|
||||||
elif 10 < capacity <= 15:
|
|
||||||
cost_per_kW = 2300
|
|
||||||
elif 15 < capacity <= 20:
|
|
||||||
cost_per_kW = 2000
|
|
||||||
elif 20 < capacity <= 10000:
|
|
||||||
cost_per_kW = 1800
|
|
||||||
else:
|
|
||||||
cost_per_kW = 1449
|
|
||||||
|
|
||||||
|
|
||||||
# Initial costs
|
|
||||||
initial_cost = capacity * cost_per_kW
|
|
||||||
|
|
||||||
# Discounted metrics initialization
|
|
||||||
discounted_generation_per_year = {}
|
|
||||||
discounted_self_consumption_per_year = {}
|
|
||||||
discounted_building_export_per_year = {}
|
|
||||||
discounted_grid_purchase_per_year = {}
|
|
||||||
discounted_total_generation = 0
|
|
||||||
discounted_total_self_consumption = 0
|
|
||||||
discounted_total_building_export = 0
|
|
||||||
discounted_total_grid_purchase = 0
|
|
||||||
discounted_annual_cost = {}
|
|
||||||
discounted_total_cost = 0
|
|
||||||
discounted_income_per_year = {}
|
|
||||||
total_discounted_income = 0
|
|
||||||
total_discounted_net_metering_income = 0
|
|
||||||
|
|
||||||
# Replacement costs calculation
|
|
||||||
replacement_cost = {
|
|
||||||
year: capacity * cost_per_kW * replacement_ratio * ((1 + inflation_rate) ** year) / (
|
|
||||||
(1 + discount_rate) ** year)
|
|
||||||
for year in year_of_replacement_list
|
|
||||||
}
|
|
||||||
|
|
||||||
# Yearly calculations
|
|
||||||
for year in range(1, period + 1):
|
|
||||||
# Apply degradation to PV generation for the current year
|
|
||||||
PV_hourly_generation_degraded = PV_hourly_generation * ((1 - degradation_rate) ** (year - 1))
|
|
||||||
|
|
||||||
# Hourly self-consumption and export considering degraded generation
|
|
||||||
building_hourly_self_consumption = np.minimum(PV_hourly_generation_degraded, building_hourly_consumption)
|
|
||||||
building_hourly_export = np.maximum(PV_hourly_generation_degraded - building_hourly_consumption, 0)
|
|
||||||
building_hourly_grid_purchase = np.maximum(building_hourly_consumption - PV_hourly_generation_degraded, 0).sum()
|
|
||||||
|
|
||||||
# Annual values
|
|
||||||
annual_self_consumption = building_hourly_self_consumption.sum()
|
|
||||||
annual_generation = PV_hourly_generation_degraded.sum()
|
|
||||||
annual_PV_export = building_hourly_export.sum()
|
|
||||||
annual_grid_purchase = building_hourly_grid_purchase.sum()
|
|
||||||
|
|
||||||
# Discounted values
|
|
||||||
discounted_generation = annual_generation / ((1 + discount_rate) ** year)
|
|
||||||
discounted_self_consumption = annual_self_consumption / ((1 + discount_rate) ** year)
|
|
||||||
discounted_building_export = annual_PV_export / ((1 + discount_rate) ** year)
|
|
||||||
discounted_grid_purchase = annual_grid_purchase / ((1 + discount_rate) ** year)
|
|
||||||
|
|
||||||
# Add to total discounted values
|
|
||||||
discounted_generation_per_year[year] = discounted_generation
|
|
||||||
discounted_self_consumption_per_year[year] = discounted_self_consumption
|
|
||||||
discounted_building_export_per_year[year] = discounted_building_export
|
|
||||||
discounted_grid_purchase_per_year[year] = discounted_grid_purchase
|
|
||||||
|
|
||||||
# Calculate total values in Life Cycle
|
|
||||||
discounted_total_generation += discounted_generation
|
|
||||||
discounted_total_self_consumption += discounted_self_consumption
|
|
||||||
discounted_total_building_export += discounted_building_export
|
|
||||||
discounted_total_grid_purchase += discounted_grid_purchase
|
|
||||||
|
|
||||||
# Annual costs
|
|
||||||
annual_opex = initial_cost * maintenance_cost_ratio * ((1 + inflation_rate) ** year) / (
|
|
||||||
(1 + discount_rate) ** year)
|
|
||||||
annual_cost = (
|
|
||||||
initial_cost if year == 1
|
|
||||||
else annual_opex + replacement_cost.get(year, 0)
|
|
||||||
)
|
|
||||||
discounted_annual_cost[year] = annual_cost
|
|
||||||
discounted_total_cost += annual_cost
|
|
||||||
|
|
||||||
# Tariff adjustment for income
|
|
||||||
inflated_grid_tariff = grid_current_tariff * ((1 + inflation_rate) ** (year - 1))
|
|
||||||
discounted_factor = ((1 + discount_rate) ** year) ** -1
|
|
||||||
|
|
||||||
# Income from self-consumption and net metering
|
|
||||||
self_consumption_income = discounted_self_consumption * inflated_grid_tariff
|
|
||||||
net_metering_income = min(annual_PV_export, first_year_grid_purchase) * inflated_grid_tariff * discounted_factor
|
|
||||||
tax_deduction_income = (
|
|
||||||
initial_cost * (1 + tax_deduct) * ((1 - tax_deduct) ** (year - 1)) * tax_deduct
|
|
||||||
)
|
|
||||||
|
|
||||||
annual_income = self_consumption_income + net_metering_income + tax_deduction_income
|
|
||||||
discounted_income_per_year[year] = annual_income
|
|
||||||
total_discounted_income += annual_income
|
|
||||||
total_discounted_net_metering_income += net_metering_income
|
|
||||||
|
|
||||||
total_discounted_income += incentive
|
|
||||||
|
|
||||||
# LCOE calculations
|
|
||||||
if discounted_total_generation == 0:
|
|
||||||
raise ValueError("Discounted generation is zero, cannot calculate LCOE.")
|
|
||||||
|
|
||||||
# To compute the LCOE for exported energy accurately,
|
|
||||||
# you should isolate the portion of the discounted income that comes only from energy exported to the grid,
|
|
||||||
# over the total discounted exported energy
|
|
||||||
# Loec of purchasing from grid is same as tariff
|
|
||||||
|
|
||||||
lcoe_pv = discounted_total_cost / discounted_total_generation
|
|
||||||
|
|
||||||
total_transaction = (
|
|
||||||
discounted_total_self_consumption +
|
|
||||||
discounted_total_building_export +
|
|
||||||
discounted_total_grid_purchase
|
|
||||||
)
|
|
||||||
|
|
||||||
# lcoe of exported electricity for net metering
|
|
||||||
lcoe_export = (
|
|
||||||
total_discounted_net_metering_income / discounted_total_building_export if discounted_total_building_export > 0 else 0)
|
|
||||||
|
|
||||||
# lcoe of the whole system combining various transactions
|
|
||||||
lcoe_system = (
|
|
||||||
(discounted_total_self_consumption / total_transaction) * lcoe_pv +
|
|
||||||
(discounted_total_grid_purchase / total_transaction) * grid_current_tariff -
|
|
||||||
(discounted_total_building_export / total_transaction) * lcoe_export
|
|
||||||
)
|
|
||||||
|
|
||||||
# NPV calculation
|
|
||||||
npv = total_discounted_income - discounted_total_cost
|
|
||||||
|
|
||||||
return {
|
|
||||||
'LCOE_PV': lcoe_pv,
|
|
||||||
'LCOE_system': lcoe_system,
|
|
||||||
'NPV': npv,
|
|
||||||
'Annual_PV_generation': first_year_generation_PV,
|
|
||||||
'Annual_building_self_consumption': first_year_self_consumption,
|
|
||||||
'Annual_grid_purchase': first_year_grid_purchase,
|
|
||||||
'Annual_PV_export': first_year_PV_export,
|
|
||||||
'Discounted_total_cost': discounted_total_cost,
|
|
||||||
'Total_discounted_income': total_discounted_income,
|
|
||||||
'Discounted_generation_per_year': discounted_generation_per_year,
|
|
||||||
'Discounted_self_consumption_per_year': discounted_self_consumption_per_year,
|
|
||||||
'Discounted_annual_cost': discounted_annual_cost,
|
|
||||||
'Discounted_income_per_year': discounted_income_per_year
|
|
||||||
}
|
|
||||||
|
|
||||||
#example
|
|
||||||
Loading…
Reference in New Issue
Block a user