Power Generation by Fuel Source (in GWh)

There is available data on the webiste of Department of Energy, Philippines that shows the power capacity and generation. The data is publicly available here: https://doe.gov.ph/key-energy-statistics-dashboards/power-capacity-and-generation.

Objective

The Objective of this analysis is to see how at a high level how close the Philippines is to reducing emissions via renewable energy sources: - Solar Energy - Wind Energy - Hydropower - Geothermal Energy - Biomass Energy - Ocean Energy

This dataset shows data over the last 20 years in GWh:

pwr <-read_csv("GenerationTablePH.csv", show_col_types = FALSE)
pwr
## # A tibble: 31 × 10
##    Years Biomass  Coal Geothermal Hydro `Natural Gas` `Oil-based` Solar  Wind
##    <dbl>   <dbl> <dbl>      <dbl> <dbl>         <dbl>       <dbl> <dbl> <dbl>
##  1  1990       0  1934       5466  6062             0       12434     0     0
##  2  1991       0  1942       5758  5145             0       12804     0     0
##  3  1992       0  1791       5700  4440             0       13939     0     0
##  4  1993       0  2015       5667  5030             0       13867     0     0
##  5  1994       0  1348       6320  5862             0       16929     0     0
##  6  1995       0  2109       6135  6232             0       19078     0     0
##  7  1996       0  4855       6534  7030             0       18288     0     0
##  8  1997       0  7363       7237  6069            12       19116     0     0
##  9  1998       0  9388       8914  5084            20       18174     0     0
## 10  1999       0 11183      10594  7840            16       11799     0     0
## # ℹ 21 more rows
## # ℹ 1 more variable: `Grand Total` <dbl>

What we can see from the data is that the Philippines does not use Ocean Energy but has all the other energy types.

Next is to understand using summary statistics to create an initial analysis framework and comparing this with data from the latest available year (2020)

summary(pwr)
##      Years         Biomass            Coal         Geothermal   
##  Min.   :1990   Min.   :   0.0   Min.   : 1348   Min.   : 5466  
##  1st Qu.:1998   1st Qu.:   0.0   1st Qu.: 8376   1st Qu.: 8076  
##  Median :2005   Median :   0.0   Median :16194   Median :10242  
##  Mean   :2005   Mean   : 201.9   Mean   :20746   Mean   : 9247  
##  3rd Qu.:2012   3rd Qu.: 189.5   3rd Qu.:30173   3rd Qu.:10454  
##  Max.   :2020   Max.   :1261.0   Max.   :58176   Max.   :11626  
##      Hydro        Natural Gas      Oil-based         Solar       
##  Min.   : 4440   Min.   :    0   Min.   : 2474   Min.   :   0.0  
##  1st Qu.: 6631   1st Qu.:   14   1st Qu.: 4766   1st Qu.:   0.0  
##  Median : 7870   Median :16366   Median : 6293   Median :   1.0  
##  Mean   : 7794   Mean   :11174   Mean   : 8817   Mean   : 204.3  
##  3rd Qu.: 9260   3rd Qu.:19547   3rd Qu.:12619   3rd Qu.:   1.5  
##  Max.   :10252   Max.   :22354   Max.   :19116   Max.   :1373.0  
##       Wind         Grand Total    
##  Min.   :   0.0   Min.   : 25649  
##  1st Qu.:   0.0   1st Qu.: 40614  
##  Median :  17.0   Median : 56568  
##  Mean   : 217.2   Mean   : 58402  
##  3rd Qu.:  81.5   3rd Qu.: 74094  
##  Max.   :1153.0   Max.   :106041
subset(pwr, Years==2020)
## # A tibble: 1 × 10
##   Years Biomass  Coal Geothermal Hydro `Natural Gas` `Oil-based` Solar  Wind
##   <dbl>   <dbl> <dbl>      <dbl> <dbl>         <dbl>       <dbl> <dbl> <dbl>
## 1  2020    1261 58176      10757  7192         19497        2474  1373  1026
## # ℹ 1 more variable: `Grand Total` <dbl>

Based on what we know coming from the overall table, we can create an initial segmentation based on the data provided split up into the following metrics:

  • Non-Renewable
    • Declining: Where mean is larger than or equal to 2020 value
    • Flat or Increasing: Where mean is less than the 2020 value
  • Renewable Sources
    • Source Available from 1990: Where 1990 value is not 0
    • Newly Introduced Source: Where the 1990 value is set at 0

Given this analysis, we can start by using R to analyze the information and properly segmenting as necessary, defining the trend and the precentage of the total:

# creating multiple vectors 
energy_source <- names(pwr)[2:9] # exclude years and grand total
latest_year <- as.numeric(pwr[which(pwr[,"Years"]=="2020"),2:9]) # get the 2020 row and get the data
first_year <- as.numeric(pwr[which(pwr[,"Years"]=="1990"),2:9]) # get the 1990 row and get the data
mean <- c(mean(pwr$Biomass),mean(pwr$Coal),mean(pwr$Geothermal),mean(pwr$Hydro),mean(pwr$`Natural Gas`),mean(pwr$`Oil-based`),mean(pwr$Solar),mean(pwr$Wind)) # get the means on per column

combined_data <- data.frame(energy_source,first_year,latest_year,mean) # create a dataframe

# build categories
ren <- energy_source[c(1,3,4,7,8)] #list of renewable energy sources

combined_data$nr_inc <- (!combined_data$energy_source %in% ren) & combined_data$latest_year >= combined_data$mean # non-renewable energy sources that are increasing

combined_data$nr_dec <- (!combined_data$energy_source %in% ren) & combined_data$latest_year < combined_data$mean # non-renewable energy sources that are decreasig

combined_data$new_ren <- combined_data$first_year==0 & combined_data$energy_source %in% ren # newly introduced renewable sources of ernergy
combined_data$existing_ren <- combined_data$first_year!=0 & combined_data$energy_source %in% ren # available source

# defining the trend if increasing. Increasing is good for renewable but bad for non-renewable
combined_data$inc_trend <- combined_data$latest_year >= combined_data$mean

# lastly, get the total as a percentage of the year 2020
combined_data$total_2020 <- as.numeric(pwr[which(pwr[,"Years"]=="2020"),"Grand Total"])
combined_data$share <- combined_data$latest_year/combined_data$total_2020

We have now completed the analysis table in for the categorization, trends and share.

combined_data
##   energy_source first_year latest_year       mean nr_inc nr_dec new_ren
## 1       Biomass          0        1261   201.9032  FALSE  FALSE    TRUE
## 2          Coal       1934       58176 20746.1613   TRUE  FALSE   FALSE
## 3    Geothermal       5466       10757  9247.3871  FALSE  FALSE   FALSE
## 4         Hydro       6062        7192  7793.8710  FALSE  FALSE   FALSE
## 5   Natural Gas          0       19497 11173.6129   TRUE  FALSE   FALSE
## 6     Oil-based      12434        2474  8816.9677  FALSE   TRUE   FALSE
## 7         Solar          0        1373   204.2581  FALSE  FALSE    TRUE
## 8          Wind          0        1026   217.2258  FALSE  FALSE    TRUE
##   existing_ren inc_trend total_2020      share
## 1        FALSE      TRUE     101756 0.01239239
## 2        FALSE      TRUE     101756 0.57172059
## 3         TRUE      TRUE     101756 0.10571367
## 4         TRUE     FALSE     101756 0.07067888
## 5        FALSE      TRUE     101756 0.19160541
## 6        FALSE     FALSE     101756 0.02431306
## 7        FALSE      TRUE     101756 0.01349306
## 8        FALSE      TRUE     101756 0.01008294

We can now analyze them one by one:

Non-renewable Energy Sources

Sources where energy sources are declining

subset(combined_data,nr_dec==TRUE,c(energy_source,mean,latest_year,share))
##   energy_source     mean latest_year      share
## 6     Oil-based 8816.968        2474 0.02431306

The Philippines used to rely more heavily on oil but has declined substantially since the start. To see how this compares, we compare it to the data in 1990:

print(paste("Oil Based in 1990: ",round(as.numeric(pwr[which(pwr[,"Years"]=="1990"),"Oil-based"]) / as.numeric(pwr[which(pwr[,"Years"]=="1990"),"Grand Total"]),2)))
## [1] "Oil Based in 1990:  0.48"

In 1990, Oil was 48% of the energy production but it is now 2% which is good sign.

Sources where energy sources are increasing

subset(combined_data,nr_inc==TRUE,c(energy_source,mean,latest_year,share))
##   energy_source     mean latest_year     share
## 2          Coal 20746.16       58176 0.5717206
## 5   Natural Gas 11173.61       19497 0.1916054

The Philippines continues to rely on Coal and Natural Gas and it accounts for 76% of the energy source. We compare the same data to how it was in 1990:

print(paste("Coal in 1990: ",round(as.numeric(pwr[which(pwr[,"Years"]=="1990"),"Coal"]) / as.numeric(pwr[which(pwr[,"Years"]=="1990"),"Grand Total"]),2))) # coal
## [1] "Coal in 1990:  0.07"
print(paste("Natural Gas in 1990: ",round(as.numeric(pwr[which(pwr[,"Years"]=="1990"),"Natural Gas"]) / as.numeric(pwr[which(pwr[,"Years"]=="1990"),"Grand Total"]),2))) # natural gas
## [1] "Natural Gas in 1990:  0"

In 1990, these two sources are not primary sources back in 1990 only accounting for 7%. Unfortunately, these two energy sources were heavily invested in and now account for majority of the energy production in the country.

Renewable Sources

Sources which were previously available

subset(combined_data,existing_ren==TRUE,c(energy_source,mean,latest_year,share))
##   energy_source     mean latest_year      share
## 3    Geothermal 9247.387       10757 0.10571367
## 4         Hydro 7793.871        7192 0.07067888

17% of the energy output in the Philippines comes from Geothermal and Hyro sources. We can compare this to the 1991 data:

print(paste("Geothermal: ",round(as.numeric(pwr[which(pwr[,"Years"]=="1990"),"Geothermal"]) / as.numeric(pwr[which(pwr[,"Years"]=="1990"),"Grand Total"]),2))) # coal
## [1] "Geothermal:  0.21"
print(paste("Hydro: ",round(as.numeric(pwr[which(pwr[,"Years"]=="1990"),"Hydro"]) / as.numeric(pwr[which(pwr[,"Years"]=="1990"),"Grand Total"]),2))) # natural gas
## [1] "Hydro:  0.23"

The interesting thing about these two sources of energy is that: - Geothermals trend is increasing versus the mean, but the min and max values are not too far from one another:

paste("Min of Geothermal: ",min(pwr$Geothermal))
## [1] "Min of Geothermal:  5466"
paste("Max of Geothermal: ",max(pwr$Geothermal))
## [1] "Max of Geothermal:  11626"
  • This suggest that geothermal energy generally falls within a small range and energy output has not increased as much over time
paste("Min of Hydro: ",min(pwr$Hydro))
## [1] "Min of Hydro:  4440"
paste("Max of Hydro: ",max(pwr$Hydro))
## [1] "Max of Hydro:  10252"
  • Hydro is decreasing as well but just like geothermal energy, it falls within a small range and output has not increased over time

We can then hypothesize that apart from these two sources coming from the environment,

New Sources

subset(combined_data,new_ren==TRUE,c(energy_source,mean,latest_year,share))
##   energy_source     mean latest_year      share
## 1       Biomass 201.9032        1261 0.01239239
## 7         Solar 204.2581        1373 0.01349306
## 8          Wind 217.2258        1026 0.01008294

Biomass, Solar and Wind are increasing over time but the share versus the total generation is quite small at 3.5%.

Conclusion

There’s a long way to go for the Philippines to rely on renewable sources of energy but will require additional investment along the way to increase generation capacity. However, it’s good news that oil-based has gone down substantially and with that use case, the same could happen to coal and natural gas once significant investment is made in renewable energy.

Category Sub Category Energy Sources Strategy
Non Renewable Declining Oil-based Decline / Stop using
Increasing Coal, Natural Gas Keep Steady or decline
Renewable Existing Geothermal, Hydro Increase/Expand Capacity
New Biomass, Solar, Wind Increase/Expand Capacity