adding average LCO-ccsinje and adjustment cost bugfixes #492

tabeado · 2023-12-08T16:44:17Z

added average levelized cost of ccsinje and cost from electricity demand for ccsinje in report.LCOE.
bugfix of average adjustment cost in report.LCOE
bugfix of "capital cost w/ adjCost" for ccsinje, dac, and some transport techs in report.Technology

tabeado · 2023-12-08T16:45:16Z

cchrisgong · 2023-12-08T17:35:17Z

R/reportLCOE.R

- v_directteinv_wadj <- readGDX(gdx,name=c("o_avgAdjCostInv"),field="l",format="first_found")[,ttot,]
+ v_directteinv <- readGDX(gdx,name=c("v_costInvTeDir","vm_costInvTeDir","v_directteinv"),field="l",format="first_found")[,ttot,] ## Total direct Investment Cost in Timestep
+ #v_directteinv_wadj <- readGDX(gdx,name=c("o_avgAdjCostInv"),field="l",format="first_found")[,ttot,] ## average adjustment cost per unit installed in period
+ v_adjteinv <- readGDX(gdx,name=c("v_costInvTeAdj"),field="l",format="first_found")[,ttot,] ## total adjustment cost in period


maybe the R variable can also be called "v_costInvTeAdj"? normally I think the convention is to go with the REMIND variable name unless it is not possible.. (like v_directteinv above this)

I agree in principle, but wanted to stick with the convention used in this part of the script. "v_costInvTeDir" turns into "v_directteinv", and I just applied the same logic. What did you mean saying "(like v_directteinv above this)"?

cchrisgong · 2023-12-08T17:46:09Z

R/reportLCOE.R

+
+    # calculate secondary Fuel cost for ccsinje 
+      te_annual_secFuel_cost <- new.magpie(getRegions(te_inv_annuity),ttot_from2005, "ccsinje" , fill=0)  
+      te_annual_secFuel_cost[,,"ccsinje"] <- setNames(-pm_prodCouple[,,"ccsinje"] * Fuel.Price[,,"seel"] * vm_co2CCS[,,"ccsinje.1"], "ccsinje") 


I am probably missing something, but why negative here? (-pm_prodCouple) oh I see nevermind it is because CCS consumes energy, would be nice to write this more explicitly (like below)

could also write a short documentation of the meanings of the three things multiplied here, I read it as
negative sign (so pm_prodCouple turns positive) * energy demand of CO2 injection (pm_prodCouple) * secondary electricity price * amount of CO2 captured (with units like below)

It says in line 249 that negative values mean own consumption, positive values mean coupled product. As I want positive cost, I need to multiply.
I can adjust the line with the units if you want to have a more detailed documentation.

cchrisgong · 2023-12-08T17:52:53Z

R/reportLCOE.R

+    # calculate secondary Fuel cost for ccsinje 
+      te_annual_secFuel_cost <- new.magpie(getRegions(te_inv_annuity),ttot_from2005, "ccsinje" , fill=0)  
+      te_annual_secFuel_cost[,,"ccsinje"] <- setNames(-pm_prodCouple[,,"ccsinje"] * Fuel.Price[,,"seel"] * vm_co2CCS[,,"ccsinje.1"], "ccsinje") 
+          # units: pm_prodCouple[,,"ccsinje"] = [TWa/GtC]; Fuel.Price = [USD2005/GtC]; vm_co2CCS= [GtC]; te_annual_secFuel_cost = [USD2005]


fuel price has unit USD2005/TWa, right?

Good catch, will adjust

cchrisgong · 2023-12-08T17:56:27Z

R/reportLCOE.R

-                                                     te_annual_OMF_cost[,,"ccsinje"])
- # distribute ccs injection cost over techs
+                                                     te_annual_OMF_cost[,,"ccsinje"] + te_annual_secFuel_cost[,,"ccsinje"])
+ # distribute ccs injection cost over techs => this does not include adjustment cost incurred! => what is the reasoning behind it?


don't quite understand this comment, do you mean why is adjustment cost not added as a OMF cost? What is the unit of total_ccsInj_cost?

Btw I don't quite get the formula here; shouldn't the total cost = annuitized (investment cost + other fixed costs (OMF, adjCost)) + running cost (energy cost, OMV cost, CO2 tax)... It looks like OMF and energy cost are here, but the others(?)

I put this in the commit accidentally. The point is that only ccs injection cost without adjustment cost are distributed across technologies. I just wondered why there is no version that includes adjustment cost as these can be very substantial

The general formula for total cost should be (and is):
annuitized (investment cost + adjCost) + OMF cost (don't need to be annuitized) + running cost (energy cost, OMV cost, CO2 tax).
The ccs cost only has investment cost + OMF cost + now also energy cost
Adjustment cost are missing, that was precisely the point of my comment

cchrisgong · 2023-12-08T17:59:39Z

R/reportLCOE.R

+              #### Carbon Transport and storage ("ccsinje")
+              setNames(te_annual_inv_cost[,ttot_from2005,"ccsinje"]/
+                         vm_co2CCS_tCO2[,,"ccsinje.1"],
+                       paste0("LCOx|average|","cco2_stored|","ccsinje","|supply-side", "|Investment Cost")),


instead of calling it LCOx, maybe something more concrete like "LCOCS" Levelized cost of carbon/Co2 storage?

why do we give "cco2_stored" in the variable name? I thought usually for readability we have something more like "Captured CO2|Injected CO2 |Supply side|..." etc. not sure what the convention is

I called it "cco2_stored" because the convention for the other items in this dimension is pe2se$all_te. In my understanding, "cco2" would be the remind element corresponding to that. However, it is not about cco2 (all captured co2) but only the part of it that is stored, thus the additional "_stored"

about LCOx vs LCOCS: I think that is a general question of how we want to call "LCO other than E" in this script. There are plans to also add LCO of CDR options. Do we want to give all of them their own abbreviation? Or do we want to have a generic LCOx and then filter according to the technologies, as it is the case with LCOE?

cchrisgong · 2023-12-08T18:01:16Z

R/reportLCOE.R


 # calculate standing system LCOE
 # divide total cost of standing system in that time step by total generation (before curtailment) in that time step
-# exception: grid and storage cost are calculate by dividing by generation after curtailment
-# convert from USD2005/MWh to USD2015/MWh (*1.2)
+# exceptions: grid and storage cost are calculate by dividing by generation after curtailment


maybe add another exception here that for CCS, it is divided by CO2 (or C?) stored

cchrisgong

Maybe would also be nice to have a trial run (or output LCOE csv available - or its screenshot), so a CCS expert can comment on the costs (not an expert here :P) maybe Jess or Anne?

tabeado · 2023-12-12T09:34:45Z

You can compare results here: /p/tmp/tabeado/remind2_clean/reportLCOE_adjustment

cchrisgong

Thanks for the nice edits! I approve this change :) Happy holiday~ :3

tabeado added 6 commits December 5, 2023 17:59

bug-fix average AdjCost + add average LCO-ccsinje

af9a3d7

Merge branch 'master' of https://github.com/tabeado/remind2

0a75957

add electricity cost for ccsinje

061fbdd

bugfix capital cost w adj cost

8d845f2

Merge branch 'master' of https://github.com/tabeado/remind2

f9c005c

increment version

9acad08

tabeado requested a review from cchrisgong December 8, 2023 16:44

cchrisgong reviewed Dec 8, 2023

View reviewed changes

cchrisgong requested changes Dec 8, 2023

View reviewed changes

tabeado added 2 commits December 19, 2023 18:42

cleanup, documentation, rename invest cost variables

b2dbc3d

Merge branch 'master' of https://github.com/pik-piam/remind2

78efc63

cchrisgong approved these changes Dec 20, 2023

View reviewed changes

tabeado added 2 commits January 4, 2024 13:59

Merge branch 'master' of https://github.com/pik-piam/remind2

2164f40

increment version

d875e21

tabeado merged commit bb3bd38 into pik-piam:master Jan 4, 2024
4 checks passed

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

adding average LCO-ccsinje and adjustment cost bugfixes #492

adding average LCO-ccsinje and adjustment cost bugfixes #492

tabeado commented Dec 8, 2023

tabeado commented Dec 8, 2023

cchrisgong Dec 8, 2023

tabeado Dec 12, 2023

cchrisgong Dec 8, 2023 •

edited

Loading

cchrisgong Dec 8, 2023

tabeado Dec 12, 2023

cchrisgong Dec 8, 2023

tabeado Dec 12, 2023

cchrisgong Dec 8, 2023

tabeado Dec 12, 2023

tabeado Dec 12, 2023

cchrisgong Dec 8, 2023

cchrisgong Dec 8, 2023 •

edited

Loading

tabeado Dec 12, 2023

tabeado Dec 12, 2023

cchrisgong Dec 8, 2023

cchrisgong left a comment

tabeado commented Dec 12, 2023 •

edited

Loading

cchrisgong left a comment

adding average LCO-ccsinje and adjustment cost bugfixes #492

adding average LCO-ccsinje and adjustment cost bugfixes #492

Conversation

tabeado commented Dec 8, 2023

tabeado commented Dec 8, 2023

Choose a reason for hiding this comment

Choose a reason for hiding this comment

cchrisgong Dec 8, 2023 • edited Loading

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

cchrisgong Dec 8, 2023 • edited Loading

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

cchrisgong left a comment

Choose a reason for hiding this comment

tabeado commented Dec 12, 2023 • edited Loading

cchrisgong left a comment

Choose a reason for hiding this comment

cchrisgong Dec 8, 2023 •

edited

Loading

cchrisgong Dec 8, 2023 •

edited

Loading

tabeado commented Dec 12, 2023 •

edited

Loading