changed calculation of capex from investment -> project development c…

…osts were lower than needed be
rl-institut · Jun 12, 2019 · bb75e1f · bb75e1f
1 parent 58ed992
commit bb75e1f
Show file tree

Hide file tree

Showing 5 changed files with 233 additions and 81 deletions.
diff --git a/D_process_input.py b/D_process_input.py
@@ -65,6 +65,7 @@ def economic_values(experiment):
                           'project']
 
         for item in component_list:
+            print(item)
             # --------------------------------------------------#
             # CAPEX without opex/a                              #
             # --------------------------------------------------#
@@ -85,7 +86,8 @@ def economic_values(experiment):
             # costs. Includes project_cost_annuity, distribution_grid_cost_annuity, maingrid_extension_cost_annuity for
             # consistency eventhough these are not used in optimization.
             experiment.update({item + '_cost_annuity': experiment[item + '_cost_annuity'] / 365*experiment['evaluated_days']})
-
+            if item == 'project':
+                print(experiment[item + '_cost_investment'], experiment[item + '_cost_capex'], experiment[item + '_cost_annuity']*365/experiment['evaluated_days'])
         return experiment
 
     def add_timeseries(experiment_s):

diff --git a/Z_economic_functions.py b/Z_economic_functions.py
@@ -12,8 +12,10 @@ def crf(project_life, wacc):
 
     def capex_from_investment(investment_t0, lifetime, project_life, wacc, tax):
         # [quantity, investment, installation, weight, lifetime, om, first_investment]
-        number_of_investments = int(round(project_life / lifetime+0.5))
-
+        if project_life == lifetime:
+            number_of_investments = 1
+        else:
+            number_of_investments = int(round(project_life / lifetime + 0.5))
         # costs with quantity and import tax at t=0
         first_time_investment = investment_t0 * (1+tax)
 
@@ -28,6 +30,8 @@ def capex_from_investment(investment_t0, lifetime, project_life, wacc, tax):
 
         # Substraction of component value at end of life with last replacement (= number_of_investments - 1)
         if number_of_investments * lifetime > project_life:
+            print(number_of_investments * lifetime)
+            print(investment_t0)
             last_investment = first_time_investment / ((1 + wacc) ** ((number_of_investments - 1) * lifetime))
             linear_depreciation_last_investment = last_investment / lifetime
             capex = capex -  linear_depreciation_last_investment * (number_of_investments * lifetime - project_life)

diff --git a/Z_output_functions.py b/Z_output_functions.py
@@ -57,9 +57,11 @@ def overall_results_title(settings, number_of_project_sites, sensitivity_array_d
 
         title_overall_results = pd.concat([title_overall_results, pd.DataFrame(columns=[
             'capacity_pv_kWp',
-            'capacity_wind_kW',
             'capacity_storage_kWh',
             'power_storage_kW',
+            'capacity_rectifier_ac_dc_kW',
+            'capacity_inverter_kW',
+            'capacity_wind_kW',
             'capacity_genset_kW',
             'capacity_pcoupling_kW',
             'consumption_fuel_annual_l',