# For general data science
library(tidyverse)

# For data cleaning
library(janitor)

# For working with datetime
library(lubridate)

# For pretty tables
library(knitr)
library(gt)

knitr::opts_chunk$set(echo = TRUE)

Oil and gas companies in Oklahoma are not required to report volumes of gas flared or vented. However, they are required to request approval from the Oklahoma Corporation Commission (OCC) to flare or vent in excess of 50,000 cubic feet per day.

To request approval operators fill out Form 1022, the Application to Flare or Vent Gas. At the bottom of the form, there are fields for the OCC to record the approved volume per day and the number of days.

The Howard Center built a database to total approved flaring and venting volumes in Oklahoma between 2012 and 2020 by extracting data from more than 800 requests submitted to the OCC. In creating the database, the Howard Center encountered duplicate flaring data due to the same application being submitted multiple times, as well as situations where there were applications for single wells that also existed on applications for leases covering multiple wells. These duplicates were filtered out of the final data set.

We downloaded the Form 1022s from the commission’s imaging website and then opened each scanned copy to manually enter the estimated and approved flaring and venting volumes. Applications were either for a single well or a lease containing multiple wells.

Some operators submitted applications to flare from a lease and then submitted multiple applications for each well on that lease. We removed these applications to avoid doubling the volumes requested and approved.

There were also applications that were not required because the volumes requested were below 50,000 cubic feet per day. We clearly delineate these as being “not required” in our data set. We still include the approved volumes (if any) from these requests when calculating our totals to best reflect how much flaring and venting the commission authorized between 2012 and 2020.

Because many of the applications were filled out by hand, we identified several data entry errors. We consulted the commission about these. Applications that were dismissed or void were excluded from the yearly totals.

Load Data

The dataset is available from the Howard Center here.

ok <- read.csv(url("https://cronkitedata.s3.amazonaws.com/howard/gaslit/ok_flaring_and_venting.csv"))

Analysis

FACT: According to the state permit applications, the commission approved the flaring and venting of 4,838,390,000 cubic feet of natural gas from 2012 to 2020.

To calculate this total, we filtered out dismissed and denied applications. Then, we extracted the year from the approved_date variable. Finally, we grouped by year and aggregated annual volumes to arrive at a grand total.

ok %>%
  filter(outcome == "approved" | outcome == "no permit required") %>%
  mutate(effective_date = as.Date(effective_date, format = "%m/%d/%Y"),
         year = year(effective_date)) %>%
  group_by(year) %>%
  summarise(Total = sum(as.numeric(total_volume_approved), na.rm = TRUE)) %>%
  gt(rowname_col = "year") %>%
  tab_stubhead(label = "Year") %>%
  tab_header(title = "Flaring and Venting Volumes Approved by OCC",
             subtitle = "Flaring and venting volumes in Mcf (thousand cubic feet)") %>%
  fmt_number(column = "Total",
             decimals = 2) %>%
  grand_summary_rows(fns = list("Total" = "sum"), 
                           columns=("Total"), 
                           formatter=gt::fmt_number, 
                           decimals=2)
Flaring and Venting Volumes Approved by OCC
Flaring and venting volumes in Mcf (thousand cubic feet)
Year Total
2012 192,050.00
2013 1,374,885.00
2014 1,084,015.00
2015 238,540.00
2016 98,670.00
2017 280,882.00
2018 664,967.00
2019 552,996.00
2020 349,385.00
Total 4,836,390.00
LS0tCnRpdGxlOiAiT2tsYWhvbWEgRmxhcmluZyBhbmQgVmVudGluZyIKYXV0aG9yOiAiTW9sbGllIEphbWlzb24iCm91dHB1dDoKICAgIGh0bWxfZG9jdW1lbnQ6CiAgICAgIGNvZGVfZm9sZGluZzogc2hvdwogICAgICBkZl9wcmludDogcGFnZWQKICAgICAgdG9jOiB0cnVlCiAgICAgIHRvY19mbG9hdDogdHJ1ZQogICAgICBjb2RlX2Rvd25sb2FkOiB0cnVlCi0tLQoKYGBge3Igc2V0dXAsIG1lc3NhZ2U9RkFMU0UsIGNsYXNzLnNvdXJjZT0iZm9sZC1oaWRlIn0KCiMgRm9yIGdlbmVyYWwgZGF0YSBzY2llbmNlCmxpYnJhcnkodGlkeXZlcnNlKQoKIyBGb3IgZGF0YSBjbGVhbmluZwpsaWJyYXJ5KGphbml0b3IpCgojIEZvciB3b3JraW5nIHdpdGggZGF0ZXRpbWUKbGlicmFyeShsdWJyaWRhdGUpCgojIEZvciBwcmV0dHkgdGFibGVzCmxpYnJhcnkoa25pdHIpCmxpYnJhcnkoZ3QpCgprbml0cjo6b3B0c19jaHVuayRzZXQoZWNobyA9IFRSVUUpCgpgYGAKCk9pbCBhbmQgZ2FzIGNvbXBhbmllcyBpbiBPa2xhaG9tYSBhcmUgbm90IHJlcXVpcmVkIHRvIHJlcG9ydCB2b2x1bWVzIG9mIGdhcyBmbGFyZWQgb3IgdmVudGVkLiBIb3dldmVyLCB0aGV5IGFyZSByZXF1aXJlZCB0byByZXF1ZXN0IGFwcHJvdmFsIGZyb20gdGhlIE9rbGFob21hIENvcnBvcmF0aW9uIENvbW1pc3Npb24gKE9DQykgdG8gZmxhcmUgb3IgdmVudCBpbiBleGNlc3Mgb2YgNTAsMDAwIGN1YmljIGZlZXQgcGVyIGRheS4gIAoKVG8gcmVxdWVzdCBhcHByb3ZhbCBvcGVyYXRvcnMgZmlsbCBvdXQgRm9ybSAxMDIyLCB0aGUgQXBwbGljYXRpb24gdG8gRmxhcmUgb3IgVmVudCBHYXMuIEF0IHRoZSBib3R0b20gb2YgdGhlIGZvcm0sIHRoZXJlIGFyZSBmaWVsZHMgZm9yIHRoZSBPQ0MgdG8gcmVjb3JkIHRoZSBhcHByb3ZlZCB2b2x1bWUgcGVyIGRheSBhbmQgdGhlIG51bWJlciBvZiBkYXlzLgoKVGhlIEhvd2FyZCBDZW50ZXIgYnVpbHQgYSBkYXRhYmFzZSB0byB0b3RhbCBhcHByb3ZlZCBmbGFyaW5nIGFuZCB2ZW50aW5nIHZvbHVtZXMgaW4gT2tsYWhvbWEgYmV0d2VlbiAyMDEyIGFuZCAyMDIwIGJ5IGV4dHJhY3RpbmcgZGF0YSBmcm9tIG1vcmUgdGhhbiA4MDAgcmVxdWVzdHMgc3VibWl0dGVkIHRvIHRoZSBPQ0MuIEluIGNyZWF0aW5nIHRoZSBkYXRhYmFzZSwgdGhlIEhvd2FyZCBDZW50ZXIgZW5jb3VudGVyZWQgZHVwbGljYXRlIGZsYXJpbmcgZGF0YSBkdWUgdG8gdGhlIHNhbWUgYXBwbGljYXRpb24gYmVpbmcgc3VibWl0dGVkIG11bHRpcGxlIHRpbWVzLCBhcyB3ZWxsIGFzIHNpdHVhdGlvbnMgd2hlcmUgdGhlcmUgd2VyZSBhcHBsaWNhdGlvbnMgZm9yIHNpbmdsZSB3ZWxscyB0aGF0IGFsc28gZXhpc3RlZCBvbiBhcHBsaWNhdGlvbnMgZm9yIGxlYXNlcyBjb3ZlcmluZyBtdWx0aXBsZSB3ZWxscy4gVGhlc2UgZHVwbGljYXRlcyB3ZXJlIGZpbHRlcmVkIG91dCBvZiB0aGUgZmluYWwgZGF0YSBzZXQuCgpXZSBkb3dubG9hZGVkIHRoZSBGb3JtIDEwMjJzIGZyb20gdGhlIGNvbW1pc3Npb27igJlzIGltYWdpbmcgW3dlYnNpdGVdKGh0dHBzOi8vaW1hZ2luZy5vY2Mub2suZ292L2ltYWdpbmcvT0dXZWxsUmVjb3Jkcy5hc3B4KSBhbmQgdGhlbiBvcGVuZWQgZWFjaCBzY2FubmVkIGNvcHkgdG8gbWFudWFsbHkgZW50ZXIgdGhlIGVzdGltYXRlZCBhbmQgYXBwcm92ZWQgZmxhcmluZyBhbmQgdmVudGluZyB2b2x1bWVzLiBBcHBsaWNhdGlvbnMgd2VyZSBlaXRoZXIgZm9yIGEgc2luZ2xlIHdlbGwgb3IgYSBsZWFzZSBjb250YWluaW5nIG11bHRpcGxlIHdlbGxzLgoKU29tZSBvcGVyYXRvcnMgc3VibWl0dGVkIGFwcGxpY2F0aW9ucyB0byBmbGFyZSBmcm9tIGEgbGVhc2UgYW5kIHRoZW4gc3VibWl0dGVkIG11bHRpcGxlIGFwcGxpY2F0aW9ucyBmb3IgZWFjaCB3ZWxsIG9uIHRoYXQgbGVhc2UuIFdlIHJlbW92ZWQgdGhlc2UgIGFwcGxpY2F0aW9ucyB0byBhdm9pZCBkb3VibGluZyB0aGUgdm9sdW1lcyByZXF1ZXN0ZWQgYW5kIGFwcHJvdmVkLiAKClRoZXJlIHdlcmUgYWxzbyBhcHBsaWNhdGlvbnMgdGhhdCB3ZXJlIG5vdCByZXF1aXJlZCBiZWNhdXNlIHRoZSB2b2x1bWVzIHJlcXVlc3RlZCB3ZXJlIGJlbG93IDUwLDAwMCBjdWJpYyBmZWV0IHBlciBkYXkuIFdlIGNsZWFybHkgZGVsaW5lYXRlIHRoZXNlIGFzIGJlaW5nIOKAnG5vdCByZXF1aXJlZOKAnSBpbiBvdXIgZGF0YSBzZXQuIFdlIHN0aWxsIGluY2x1ZGUgdGhlIGFwcHJvdmVkIHZvbHVtZXMgKGlmIGFueSkgZnJvbSB0aGVzZSByZXF1ZXN0cyB3aGVuIGNhbGN1bGF0aW5nIG91ciB0b3RhbHMgdG8gYmVzdCByZWZsZWN0IGhvdyBtdWNoIGZsYXJpbmcgYW5kIHZlbnRpbmcgdGhlIGNvbW1pc3Npb24gYXV0aG9yaXplZCBiZXR3ZWVuIDIwMTIgYW5kIDIwMjAuIAoKQmVjYXVzZSBtYW55IG9mIHRoZSBhcHBsaWNhdGlvbnMgd2VyZSBmaWxsZWQgb3V0IGJ5IGhhbmQsIHdlIGlkZW50aWZpZWQgc2V2ZXJhbCBkYXRhIGVudHJ5IGVycm9ycy4gV2UgY29uc3VsdGVkIHRoZSBjb21taXNzaW9uIGFib3V0IHRoZXNlLiBBcHBsaWNhdGlvbnMgdGhhdCB3ZXJlIGRpc21pc3NlZCBvciB2b2lkIHdlcmUgZXhjbHVkZWQgZnJvbSB0aGUgeWVhcmx5IHRvdGFscy4gCgojIExvYWQgRGF0YQoKVGhlIGRhdGFzZXQgaXMgYXZhaWxhYmxlIGZyb20gdGhlIEhvd2FyZCBDZW50ZXIgW2hlcmVdKGh0dHBzOi8vZHJpdmUuZ29vZ2xlLmNvbS9maWxlL2QvMVhxejVKTHJzVWE1T0pINTVNRzVKR19idEMwMlQ2bENNL3ZpZXc/dXNwPXNoYXJpbmcpLiAKCmBgYHtyfQoKb2sgPC0gcmVhZC5jc3YodXJsKCJodHRwczovL2Nyb25raXRlZGF0YS5zMy5hbWF6b25hd3MuY29tL2hvd2FyZC9nYXNsaXQvb2tfZmxhcmluZ19hbmRfdmVudGluZy5jc3YiKSkKCmBgYAoKIyBBbmFseXNpcwoKRkFDVDogQWNjb3JkaW5nIHRvIHRoZSBzdGF0ZSBwZXJtaXQgYXBwbGljYXRpb25zLCB0aGUgY29tbWlzc2lvbiBhcHByb3ZlZCB0aGUgZmxhcmluZyBhbmQgdmVudGluZyBvZiA0LDgzOCwzOTAsMDAwIGN1YmljIGZlZXQgb2YgbmF0dXJhbCBnYXMgZnJvbSAyMDEyIHRvIDIwMjAuCgpUbyBjYWxjdWxhdGUgdGhpcyB0b3RhbCwgd2UgZmlsdGVyZWQgb3V0IGRpc21pc3NlZCBhbmQgZGVuaWVkIGFwcGxpY2F0aW9ucy4gVGhlbiwgd2UgZXh0cmFjdGVkIHRoZSB5ZWFyIGZyb20gdGhlIGBhcHByb3ZlZF9kYXRlYCB2YXJpYWJsZS4gRmluYWxseSwgIHdlIGdyb3VwZWQgYnkgeWVhciBhbmQgYWdncmVnYXRlZCBhbm51YWwgdm9sdW1lcyB0byBhcnJpdmUgYXQgYSBncmFuZCB0b3RhbC4gCgoKYGBge3J9CgpvayAlPiUKICBmaWx0ZXIob3V0Y29tZSA9PSAiYXBwcm92ZWQiIHwgb3V0Y29tZSA9PSAibm8gcGVybWl0IHJlcXVpcmVkIikgJT4lCiAgbXV0YXRlKGVmZmVjdGl2ZV9kYXRlID0gYXMuRGF0ZShlZmZlY3RpdmVfZGF0ZSwgZm9ybWF0ID0gIiVtLyVkLyVZIiksCiAgICAgICAgIHllYXIgPSB5ZWFyKGVmZmVjdGl2ZV9kYXRlKSkgJT4lCiAgZ3JvdXBfYnkoeWVhcikgJT4lCiAgc3VtbWFyaXNlKFRvdGFsID0gc3VtKGFzLm51bWVyaWModG90YWxfdm9sdW1lX2FwcHJvdmVkKSwgbmEucm0gPSBUUlVFKSkgJT4lCiAgZ3Qocm93bmFtZV9jb2wgPSAieWVhciIpICU+JQogIHRhYl9zdHViaGVhZChsYWJlbCA9ICJZZWFyIikgJT4lCiAgdGFiX2hlYWRlcih0aXRsZSA9ICJGbGFyaW5nIGFuZCBWZW50aW5nIFZvbHVtZXMgQXBwcm92ZWQgYnkgT0NDIiwKICAgICAgICAgICAgIHN1YnRpdGxlID0gIkZsYXJpbmcgYW5kIHZlbnRpbmcgdm9sdW1lcyBpbiBNY2YgKHRob3VzYW5kIGN1YmljIGZlZXQpIikgJT4lCiAgZm10X251bWJlcihjb2x1bW4gPSAiVG90YWwiLAogICAgICAgICAgICAgZGVjaW1hbHMgPSAyKSAlPiUKICBncmFuZF9zdW1tYXJ5X3Jvd3MoZm5zID0gbGlzdCgiVG90YWwiID0gInN1bSIpLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgY29sdW1ucz0oIlRvdGFsIiksIAogICAgICAgICAgICAgICAgICAgICAgICAgICBmb3JtYXR0ZXI9Z3Q6OmZtdF9udW1iZXIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICBkZWNpbWFscz0yKQoKCgpgYGAK