Crude-Oil Boom and Coming Bust in Colorado

L. David Roper
9 October, 2016
World Fossil Fuels Depletion



The technique of fractionating ("fracking") shale/dolomite formations. A "boom" in nonrenewable resource extraction from the Earth is defined as a time period in which extraction is occurring very fast in a given area; thus, many workers come in from outside the area to man the drilling rigs, to build housing for the oil workers and to provide other services for the increased population.

This article shows mathematically that the current Colorado crude-oil boom will become a bust within a decade. A "bust" in nonrenewable resource extraction from the Earth is defined to begin at the time when extraction of the resource peaks and then falls to negligible amounts over a time period.

Crude-Oil Extraction Data

The U.S. Energy Information Administration gives monthly and annual crude-oil extraction data for Oklahoma since 1981. Earlier data are available back to 1905.

Those data are fitted by a depletion function, the Verhulst function, in this study to determine when the extraction will peak.

The data and the fits to the data are given in a later section.

Crude-Oil Reserves

A reliable estimate of reserves is needed to fit extraction data by a function for projecting into the future for a nonrenewable resource Here is a good definition of reserves of a nonrenewable resource.

The U.S. Energy Information Administration gives reserves estimates from 2006 to 2011 for crude-oil extraction in Colorado , which are shown here by black dots for years 2005 to 2014:

The curve is a fit to the 4 last data using the Verhulst function described in the next section, assuming that the curve will be symmetrical. Since the reserves estimates have been rising each year since 2005, the fit is done to get a rough estimate of the peak value of the reserves estimate in the future, which is ~4.5 x 109 barrels.

Go to the top.

Verhulst-Function Fit to Crude-Oil Extraction Data

The depletion function that is used in this article is the Verhulst function:

The asymmetry parameter, n, must be greater than 0.


The maximum of P(t) occurs at , which yields the peak value  .

For the symmetric case (n=1):  and .

For a depletion situation for which there are N peaks the depletion function is:


When a peak is symmetrical, the Verhulst function simplifies to


The depletion curve for crude-oil extraction for Oklahoma has essentially three peaks; so N=3 in the depletion equation. The fit to the extraction data using the reserves peak, ~5.2 x 109 barrels, given in a previous section is:

The final large peak in the red curve is assumed to be symmetric. Asymmetry would shift the peak backward for n > 1 (depletion curve skewed forward with a smaller peak) and forward for n < 1 (depletion curve skewed backward with a larger peak). The blue curve is for n=5, a value that is close to what J. David Hughes finds in his detailed micro-analysis of extraction of tight oil by fracking in the U.S.

The area under the curve is equal to the amount already extracted (~2.4 x 109 barrels) plus the reserves value (~4.5 x 109 barrels).

Even if this very high reserves value (~4.5 x 109 barrels) is correct, the final peak occurs at year ~2021. The onset of the bust could be extended out to later years by imposing environmental regulations and/or taxes on the extraction of natural gas, thereby reducing the extraction rate.

Go to the top.

L. David Roper interdisciplinary studies
World Fossil Fuels Depletion

L. David Roper,
9 October, 2016