# For general data science
library(tidyverse)
# For working with datetime
library(lubridate)
# For pretty tables
library(gt)
knitr::opts_chunk$set(echo = TRUE)
Introduction
The Utah Division of Oil, Gas and Mining requires oil and gas operators on producing, shut-in and temporarily abandoned wells to fill a form titled Form 11
(https://oilgas.ogm.utah.gov/pub/report_forms/frm_11_2005.pdf) with information on flared and vented gas.
Data from the monthly reports (2012-present) are available at the division’s website under its Data Research Center tab here. To download the data for our selected tiem frame, click Disposition Data (2009-2014; can open in Excel),
Disposition Data (2015-2019; can open in Excel),
and Disposition Data (2020-2024; can open in Excel)
.
Load Data
We compiled and cleaned the data from the state regulatory agency and saved it as a new dataset available from the Howard Center here. (For more information about our data cleaning, please contact: howardcenter@asu.edu)
ut <- read.csv(url("https://cronkitedata.s3.amazonaws.com/howard/gaslit/ut_flaring_and_venting.csv"))
Analysis
FACT: “According to state data, oil and gas companies in Utah flared and vented 11,390,778,000 cubic feet of natural gas between 2012 and 2020.”
We first first for the years covered (2012-2020), then group by year to calculate annual flared and vented volumes. Using the gt
package, we calculate a grand total.
ut_annualtotals <- ut %>%
filter(year >= 2012 & year <= 2020) %>%
group_by(year) %>%
summarise(Total = sum(vol_vented_flared, na.rm = TRUE))
ut_annualtotals %>%
gt(rowname_col = "year") %>%
tab_stubhead(label = "Year") %>%
tab_header(title = "Flaring Volumes Reported to UDOGM",
subtitle = "Flaring volumes in Mcf (thousand cubic feet)") %>%
fmt_number(column = "Total",
decimals = 0) %>%
grand_summary_rows(fns = list("Total" = "sum"),
columns = "Total",
formatter=gt::fmt_number,
decimals=0)
Year |
Total |
2012 |
1,660,554 |
2013 |
2,424,249 |
2014 |
2,013,742 |
2015 |
987,235 |
2016 |
597,587 |
2017 |
1,216,323 |
2018 |
951,458 |
2019 |
944,164 |
2020 |
595,466 |
Total |
11,390,778 |
LS0tCnRpdGxlOiAiVXRhaCBmbGFyaW5nIGFuZCB2ZW50aW5nIHRvdGFscyIKYXV0aG9yOiAiWm9oYSBUdW5pbywgU2FyYWggSHVudCIKb3V0cHV0OiAKICBodG1sX2RvY3VtZW50OgogICAgICBjb2RlX2ZvbGRpbmc6IHNob3cKICAgICAgZGZfcHJpbnQ6IHBhZ2VkCiAgICAgIHRvYzogdHJ1ZQogICAgICB0b2NfZmxvYXQ6IHRydWUKICAgICAgY29kZV9kb3dubG9hZDogdHJ1ZQotLS0KCmBgYHtyIHNldHVwLCBtZXNzYWdlID0gRkFMU0UsIGNsYXNzLnNvdXJjZT0iZm9sZC1oaWRlIn0KCiMgRm9yIGdlbmVyYWwgZGF0YSBzY2llbmNlCmxpYnJhcnkodGlkeXZlcnNlKQoKIyBGb3Igd29ya2luZyB3aXRoIGRhdGV0aW1lCmxpYnJhcnkobHVicmlkYXRlKQoKIyBGb3IgcHJldHR5IHRhYmxlcwpsaWJyYXJ5KGd0KQoKa25pdHI6Om9wdHNfY2h1bmskc2V0KGVjaG8gPSBUUlVFKQoKYGBgCgojIEludHJvZHVjdGlvbgoKVGhlIFV0YWggRGl2aXNpb24gb2YgT2lsLCBHYXMgYW5kIE1pbmluZyByZXF1aXJlcyBvaWwgYW5kIGdhcyBvcGVyYXRvcnMgb24gcHJvZHVjaW5nLCBzaHV0LWluIGFuZCB0ZW1wb3JhcmlseSBhYmFuZG9uZWQgd2VsbHMgdG8gZmlsbCBhIGZvcm0gdGl0bGVkIGBGb3JtIDExYCAoaHR0cHM6Ly9vaWxnYXMub2dtLnV0YWguZ292L3B1Yi9yZXBvcnRfZm9ybXMvZnJtXzExXzIwMDUucGRmKSB3aXRoIGluZm9ybWF0aW9uIG9uIGZsYXJlZCBhbmQgdmVudGVkIGdhcy4gCgpEYXRhIGZyb20gdGhlIG1vbnRobHkgcmVwb3J0cyAoMjAxMi1wcmVzZW50KSBhcmUgYXZhaWxhYmxlIGF0IHRoZSBkaXZpc2lvbidzIHdlYnNpdGUgdW5kZXIgaXRzIERhdGEgUmVzZWFyY2ggQ2VudGVyIHRhYiBbaGVyZV0oaHR0cHM6Ly9vaWxnYXMub2dtLnV0YWguZ292L29pbGdhc3dlYi9kYXRhLWNlbnRlci9kYy1tYWluLnhodG1sKS4gVG8gZG93bmxvYWQgdGhlIGRhdGEgZm9yIG91ciBzZWxlY3RlZCB0aWVtIGZyYW1lLCBjbGljayBgRGlzcG9zaXRpb24gRGF0YSAoMjAwOS0yMDE0OyBjYW4gb3BlbiBpbiBFeGNlbCksYCBgRGlzcG9zaXRpb24gRGF0YSAoMjAxNS0yMDE5OyBjYW4gb3BlbiBpbiBFeGNlbCksYCBhbmQgYERpc3Bvc2l0aW9uIERhdGEgKDIwMjAtMjAyNDsgY2FuIG9wZW4gaW4gRXhjZWwpYC4gCgoKIyBMb2FkIERhdGEKCldlIGNvbXBpbGVkIGFuZCBjbGVhbmVkIHRoZSBkYXRhIGZyb20gdGhlIHN0YXRlIHJlZ3VsYXRvcnkgYWdlbmN5IGFuZCBzYXZlZCBpdCBhcyBhIG5ldyBkYXRhc2V0IGF2YWlsYWJsZSBmcm9tIHRoZSBIb3dhcmQgQ2VudGVyIFtoZXJlXShodHRwczovL2RyaXZlLmdvb2dsZS5jb20vZmlsZS9kLzFSM0ZQQW9GWkRpSnpkbk1KcGxqVkhybWFrLVRadTRkQS92aWV3P3VzcD1zaGFyaW5nKS4gKEZvciBtb3JlIGluZm9ybWF0aW9uIGFib3V0IG91ciBkYXRhIGNsZWFuaW5nLCBwbGVhc2UgY29udGFjdDogaG93YXJkY2VudGVyQGFzdS5lZHUpIAoKYGBge3IgbG9hZF9hbmRfY2xlYW5fZGF0YSwgZWNobz1UUlVFfQoKdXQgPC0gcmVhZC5jc3YodXJsKCJodHRwczovL2Nyb25raXRlZGF0YS5zMy5hbWF6b25hd3MuY29tL2hvd2FyZC9nYXNsaXQvdXRfZmxhcmluZ19hbmRfdmVudGluZy5jc3YiKSkKCmBgYAoKCiMgQW5hbHlzaXMKCkZBQ1Q6ICJBY2NvcmRpbmcgdG8gc3RhdGUgZGF0YSwgb2lsIGFuZCBnYXMgY29tcGFuaWVzIGluIFV0YWggZmxhcmVkIGFuZCB2ZW50ZWQgMTEsMzkwLDc3OCwwMDAgY3ViaWMgZmVldCBvZiBuYXR1cmFsIGdhcyBiZXR3ZWVuIDIwMTIgYW5kIDIwMjAuIgoKV2UgZmlyc3QgZmlyc3QgZm9yIHRoZSB5ZWFycyBjb3ZlcmVkICgyMDEyLTIwMjApLCB0aGVuIGdyb3VwIGJ5IHllYXIgdG8gY2FsY3VsYXRlIGFubnVhbCBmbGFyZWQgYW5kIHZlbnRlZCB2b2x1bWVzLiBVc2luZyB0aGUgYGd0YCBwYWNrYWdlLCB3ZSBjYWxjdWxhdGUgYSBncmFuZCB0b3RhbC4gCgoKYGBge3J9Cgp1dF9hbm51YWx0b3RhbHMgPC0gdXQgJT4lCiAgZmlsdGVyKHllYXIgPj0gMjAxMiAmIHllYXIgPD0gMjAyMCkgJT4lCiAgZ3JvdXBfYnkoeWVhcikgJT4lCiAgc3VtbWFyaXNlKFRvdGFsID0gc3VtKHZvbF92ZW50ZWRfZmxhcmVkLCBuYS5ybSA9IFRSVUUpKQoKCnV0X2FubnVhbHRvdGFscyAlPiUKICAgZ3Qocm93bmFtZV9jb2wgPSAieWVhciIpICU+JQogICB0YWJfc3R1YmhlYWQobGFiZWwgPSAiWWVhciIpICU+JQogICB0YWJfaGVhZGVyKHRpdGxlID0gIkZsYXJpbmcgVm9sdW1lcyBSZXBvcnRlZCB0byBVRE9HTSIsCiAgICAgICAgICAgICAgc3VidGl0bGUgPSAiRmxhcmluZyB2b2x1bWVzIGluIE1jZiAodGhvdXNhbmQgY3ViaWMgZmVldCkiKSAlPiUKICAgZm10X251bWJlcihjb2x1bW4gPSAiVG90YWwiLAogICAgICAgICAgICAgIGRlY2ltYWxzID0gMCkgJT4lCiAgIGdyYW5kX3N1bW1hcnlfcm93cyhmbnMgPSBsaXN0KCJUb3RhbCIgPSAic3VtIiksCiAgICAgICAgICAgICAgICAgICAgICBjb2x1bW5zID0gIlRvdGFsIiwKICAgICAgICAgICAgICAgICAgICAgIGZvcm1hdHRlcj1ndDo6Zm10X251bWJlciwKICAgICAgICAgICAgICAgICAgICAgIGRlY2ltYWxzPTApCgoKYGBgCgo=