alex/CIFS_Storage_Demand_Model
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Added start of analysis work

Browse Source
This commit is contained in:
Alex 2026-01-23 14:09:14 -07:00
parent 16c2e872db
commit b5ccf98aac
1 changed files with 20 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
Analysis.r Normal file
View File

@ -0,0 +1,20 @@
library(tidyverse)
library(RcppRoll)
RES <- readRDS("Results/Storage_Values_by_Facility_and_Variable_Discounts.Rds")
RES[[1]]
names(RES)
read_csv("Data/Curie_Spent_Fuel_Site_Totals.csv")
CDATA <- RES$Per_5 %>% left_join(read_csv("Data/Curie_Spent_Fuel_Site_Totals.csv")) %>% mutate(Revenue=parse_number(Revenue),Year=parse_number(Year),Marginal_Value=Revenue/Total_Tons) %>% select(Facility,Year,Revenue,Marginal_Value,Total_Tons)
CDATA <- CDATA %>% group_by(Year) %>% arrange(Year,desc(Marginal_Value)) %>% mutate(CIFS_total=cumsum(Total_Tons)) %>% ungroup
CDATA <- CDATA %>% group_by(Year) %>% mutate(CIFS_20=CIFS_total<=20000,CIFS_40=CIFS_total<=40000,CIFS_80=CIFS_total<=80000,CIFS_150=CIFS_total<=150000) %>% ungroup
SUMMARY <- rbind(CDATA %>% filter(CIFS_20) %>% group_by(Year) %>% summarize(Revenue=sum(Revenue),Marginal_Value=min(Marginal_Value),Total_tons=sum(Total_Tons),CIFS_Size=20000),
CDATA %>% filter(CIFS_40) %>% group_by(Year) %>% summarize(Revenue=sum(Revenue),Marginal_Value=min(Marginal_Value),Total_tons=sum(Total_Tons),CIFS_Size=40000),
CDATA %>% filter(CIFS_80) %>% group_by(Year) %>% summarize(Revenue=sum(Revenue),Marginal_Value=min(Marginal_Value),Total_tons=sum(Total_Tons),CIFS_Size=80000),
CDATA %>% filter(CIFS_150) %>% group_by(Year) %>% summarize(Revenue=sum(Revenue),Marginal_Value=min(Marginal_Value),Total_tons=sum(Total_Tons),CIFS_Size=150000))
###Graph the waiting value (value added by higher need, vs cost of waiting)
COST <- 500
SUMMARY <- SUMMARY %>% group_by(CIFS_Size) %>% mutate(Waiting_Rev=lead(Revenue,1)-Revenue,Waiting_Rev=ifelse(is.na(Waiting_Rev),0,Waiting_Rev),Waiting_Cost=0.05*(Revenue-COST),Waiting_Profit=Waiting_Rev-Waiting_Cost)
ggsave(ggplot(SUMMARY ,aes(x=Year,y=Waiting_Profit,color=as.factor(CIFS_Size)))+geom_point()+ scale_x_continuous(breaks = seq(1960,2083, by=5))+theme_bw(),filename="test.png",width=20)
ggsave(ggplot(SUMMARY ,aes(x=Year,y=Revenue,color=as.factor(CIFS_Size)))+geom_point()+geom_smooth()+ scale_x_continuous(breaks = seq(1960,2083, by=5)),filename="test.png",width=20)