Hubbert linearization

Last updated

The Hubbert linearization is a way to plot production data to estimate two important parameters of a Hubbert curve, the approximated production rate of a nonrenewable resource following a logistic distribution:

Contents

The linearization technique was introduced by Marion King Hubbert in his 1982 review paper. [1] The Hubbert curve [2] is the first derivative of a logistic function, which has been used for modeling the depletion of crude oil in particular, the depletion of finite mineral resources in general [3] and also population growth patterns. [4]

Example of a Hubbert Linearization on the US Lower-48 crude oil production. HubbertLin US Lower48.svg
Example of a Hubbert Linearization on the US Lower-48 crude oil production.

Principle

The first step of the Hubbert linearization consists of plotting the yearly production data (P in bbl/y) as a fraction of the cumulative production (Q in bbl) on the vertical axis and the cumulative production on the horizontal axis. This representation exploits the linear property of the logistic differential equation:

with

We can rewrite (1) as the following:

The above relation is a line equation in the P/Q versus Q plane. Consequently, a linear regression on the data points gives us an estimate of the line slope calculated by -k/URR and intercept from which we can derive the Hubbert curve parameters:

Examples

Global oil production

The geologist Kenneth S. Deffeyes applied this technique in 2005 to make a prediction about the peak of overall oil production at the end of the same year, which has since been found to be premature. [5] He did not make a distinction between "conventional" and "non-conventional" oil produced by fracturing, aka tight oil, which has continued further growth in oil production. However, since 2005 conventional oil production has not grown anymore. [6]

US oil production

The charts below gives an example of the application of the Hubbert Linearization technique in the case of the US Lower-48 oil production. The fit of a line using the data points from 1956 to 2005 (in green) gives a URR of 199 Gb and a logistic growth rate of 6%.

Norway oil production

The Norwegian Hubbert linearization estimates an URR = 30 Gb and a logistic growth rate of k = 17%.

Alternative techniques

Second Hubbert linearization

The Hubbert linearization principle can be extended to the first derivatives of the production rate [7] by computing the derivative of (2):

The left term, the rate of change of production per current production, is often called the decline rate. The decline curve is a line that starts at +k, crosses zero at URR/2 and ends at −k. Thus, we can derive the Hubbert curve parameters:

Hubbert parabola

This representation was proposed by Roberto Canogar [8] and applied to the oil depletion problem. It is equation (2) multiplied by Q.

The parabola starts from the origin (0,0) and passes through (URR,0). Data points until t are used by the least squares fitting method to find an estimate for URR.

Logit transform

David Rutledge applied the logit transform for the analysis of coal production data, [9] which often has a worse signal-to-noise ratio than the production data for hydrocarbons. The integrative nature of cumulation serves as a low pass, filtering noise effects. The production data is fitted to the logistic curve after transformation using e(t) as normalized exhaustion parameter going from 0 to 1.

The value of URR is varied so that the linearized logit gives a best fit with a maximal coefficient of determination .


Related Research Articles

<span class="mw-page-title-main">Logistic function</span> S-shaped curve

A logistic function or logistic curve is a common S-shaped curve with the equation

The Hubbert curve is an approximation of the production rate of a resource over time. It is a symmetric logistic distribution curve, often confused with the "normal" gaussian function. It first appeared in "Nuclear Energy and the Fossil Fuels," geologist M. King Hubbert's 1956 presentation to the American Petroleum Institute, as an idealized symmetric curve, during his tenure at the Shell Oil Company. It has gained a high degree of popularity in the scientific community for predicting the depletion of various natural resources. The curve is the main component of Hubbert peak theory, which has led to the rise of peak oil concerns. Basing his calculations on the peak of oil well discovery in 1948, Hubbert used his model in 1956 to create a curve which predicted that oil production in the contiguous United States would peak around 1970.

<span class="mw-page-title-main">Logistic regression</span> Statistical model for a binary dependent variable

In statistics, the logistic model is a statistical model that models the log-odds of an event as a linear combination of one or more independent variables. In regression analysis, logistic regression is estimating the parameters of a logistic model. Formally, in binary logistic regression there is a single binary dependent variable, coded by an indicator variable, where the two values are labeled "0" and "1", while the independent variables can each be a binary variable or a continuous variable. The corresponding probability of the value labeled "1" can vary between 0 and 1, hence the labeling; the function that converts log-odds to probability is the logistic function, hence the name. The unit of measurement for the log-odds scale is called a logit, from logistic unit, hence the alternative names. See § Background and § Definition for formal mathematics, and § Example for a worked example.

<span class="mw-page-title-main">M. King Hubbert</span> American geologist (1903–1989)

Marion King Hubbert was an American geologist and geophysicist. He worked at the Shell research lab in Houston, Texas. He made several important contributions to geology, geophysics, and petroleum geology, most notably the Hubbert curve and Hubbert peak theory, with important political ramifications. He was often referred to as "M. King Hubbert" or "King Hubbert".

<span class="mw-page-title-main">Curve fitting</span> Process of constructing a curve that has the best fit to a series of data points

Curve fitting is the process of constructing a curve, or mathematical function, that has the best fit to a series of data points, possibly subject to constraints. Curve fitting can involve either interpolation, where an exact fit to the data is required, or smoothing, in which a "smooth" function is constructed that approximately fits the data. A related topic is regression analysis, which focuses more on questions of statistical inference such as how much uncertainty is present in a curve that is fitted to data observed with random errors. Fitted curves can be used as an aid for data visualization, to infer values of a function where no data are available, and to summarize the relationships among two or more variables. Extrapolation refers to the use of a fitted curve beyond the range of the observed data, and is subject to a degree of uncertainty since it may reflect the method used to construct the curve as much as it reflects the observed data.

<span class="mw-page-title-main">Hubbert peak theory</span> One of the primary theories on peak oil

The Hubbert peak theory says that for any given geographical area, from an individual oil-producing region to the planet as a whole, the rate of petroleum production tends to follow a bell-shaped curve. It is one of the primary theories on peak oil.

In statistics, a generalized linear model (GLM) is a flexible generalization of ordinary linear regression. The GLM generalizes linear regression by allowing the linear model to be related to the response variable via a link function and by allowing the magnitude of the variance of each measurement to be a function of its predicted value.

<span class="mw-page-title-main">Regression analysis</span> Set of statistical processes for estimating the relationships among variables

In statistical modeling, regression analysis is a set of statistical processes for estimating the relationships between a dependent variable and one or more independent variables. The most common form of regression analysis is linear regression, in which one finds the line that most closely fits the data according to a specific mathematical criterion. For example, the method of ordinary least squares computes the unique line that minimizes the sum of squared differences between the true data and that line. For specific mathematical reasons, this allows the researcher to estimate the conditional expectation of the dependent variable when the independent variables take on a given set of values. Less common forms of regression use slightly different procedures to estimate alternative location parameters or estimate the conditional expectation across a broader collection of non-linear models.

<span class="mw-page-title-main">Peak oil</span> Point in time when the maximum rate of petroleum extraction is reached

Peak oil is the theorized point in time when the maximum rate of global oil production will occur, after which oil production will begin an irreversible decline. The primary concern of peak oil is that global transportation heavily relies upon the use of gasoline and diesel fuel. Switching transportation to electric vehicles, biofuels, or more fuel-efficient forms of travel may help reduce oil demand.

<span class="mw-page-title-main">Log–log plot</span> 2D graphic with logarithmic scales on both axes

In science and engineering, a log–log graph or log–log plot is a two-dimensional graph of numerical data that uses logarithmic scales on both the horizontal and vertical axes. Power functions – relationships of the form – appear as straight lines in a log–log graph, with the exponent corresponding to the slope, and the coefficient corresponding to the intercept. Thus these graphs are very useful for recognizing these relationships and estimating parameters. Any base can be used for the logarithm, though most commonly base 10 are used.

Oil depletion is the decline in oil production of a well, oil field, or geographic area. The Hubbert peak theory makes predictions of production rates based on prior discovery rates and anticipated production rates. Hubbert curves predict that the production curves of non-renewing resources approximate a bell curve. Thus, according to this theory, when the peak of production is passed, production rates enter an irreversible decline.

<span class="mw-page-title-main">Dose–response relationship</span> Measure of organism response to stimulus

The dose–response relationship, or exposure–response relationship, describes the magnitude of the response of an organism, as a function of exposure to a stimulus or stressor after a certain exposure time. Dose–response relationships can be described by dose–response curves. This is explained further in the following sections. A stimulus response function or stimulus response curve is defined more broadly as the response from any type of stimulus, not limited to chemicals.

Latent growth modeling is a statistical technique used in the structural equation modeling (SEM) framework to estimate growth trajectories. It is a longitudinal analysis technique to estimate growth over a period of time. It is widely used in the field of psychology, behavioral science, education and social science. It is also called latent growth curve analysis. The latent growth model was derived from theories of SEM. General purpose SEM software, such as OpenMx, lavaan, AMOS, Mplus, LISREL, or EQS among others may be used to estimate growth trajectories.

Jean H. Laherrère is a French petroleum engineer and consultant, best known as the co-author of an influential 1998 Scientific American article entitled The End of Cheap Oil.

Peak coal is the peak consumption or production of coal by a human community. Peak coal can be driven by peak demand or peak supply. Historically, it was widely believed that the supply-side would eventually drive peak coal due to the depletion of coal reserves. However, since the increasing global efforts to limit climate change, peak coal has been driven by demand. This is due in large part to the rapid expansion of natural gas and renewable energy. As of 2024 over 40% of all energy sector CO2 emissions are from coal, and many countries have pledged to phase-out coal.

<span class="mw-page-title-main">Predicting the timing of peak oil</span>

Peak oil is the point at which oil production, sometimes including unconventional oil sources, hits its maximum. Predicting the timing of peak oil involves estimation of future production from existing oil fields as well as future discoveries. The most influential production model is Hubbert peak theory, first proposed in the 1950s. The effect of peak oil on the world economy remains controversial.

<span class="mw-page-title-main">Sustainability measurement</span> Quantitative basis for the informed management of sustainability

Sustainability measurement is a set of frameworks or indicators used to measure how sustainable something is. This includes processes, products, services and businesses. Sustainability is difficult to quantify. It may even be impossible to measure as there is no fixed definition. To measure sustainability, frameworks and indicators consider environmental, social and economic domains. The metrics vary by use case and are still evolving. They include indicators, benchmarks and audits. They include sustainability standards and certification systems like Fairtrade and Organic. They also involve indices and accounting. They can include assessment, appraisal and other reporting systems. The metrics are used over a wide range of spatial and temporal scales. For organizations, sustainability measures include corporate sustainability reporting and Triple Bottom Line accounting. For countries, they include estimates of the quality of sustainability governance or quality of life measures, or environmental assessments like the Environmental Sustainability Index and Environmental Performance Index. Some methods let us track sustainable development. These include the UN Human Development Index and ecological footprints.

Peak minerals marks the point in time when the largest production of a mineral will occur in an area, with production declining in subsequent years. While most mineral resources will not be exhausted in the near future, global extraction and production has become more challenging. Miners have found ways over time to extract deeper and lower grade ores with lower production costs. More than anything else, declining average ore grades are indicative of ongoing technological shifts that have enabled inclusion of more 'complex' processing – in social and environmental terms as well as economic – and structural changes in the minerals exploration industry and these have been accompanied by significant increases in identified Mineral Reserves.

In statistics, the class of vector generalized linear models (VGLMs) was proposed to enlarge the scope of models catered for by generalized linear models (GLMs). In particular, VGLMs allow for response variables outside the classical exponential family and for more than one parameter. Each parameter can be transformed by a link function. The VGLM framework is also large enough to naturally accommodate multiple responses; these are several independent responses each coming from a particular statistical distribution with possibly different parameter values.

A quantile-parameterized distribution (QPD) is a probability distributions that is directly parameterized by data. They were created to meet the need for easy-to-use continuous probability distributions flexible enough to represent a wide range of uncertainties, such as those commonly encountered in business, economics, engineering, and science. Because QPDs are directly parameterized by data, they have the practical advantage of avoiding the intermediate step of parameter estimation, a time-consuming process that typically requires non-linear iterative methods to estimate probability-distribution parameters from data. Some QPDs have virtually unlimited shape flexibility and closed-form moments as well.

References

  1. Hubbert, M. King (1982). "Techniques of Prediction as Applied to the Production of Oil and Gas". In Gass, Saul I. (ed.). Oil and Gas Supply Modeling (PDF) (proceedings of a symposium held at the Department of Commerce, Washington, DC, June 18–20, 1980). NBS Special Publication 631. Washington (DC): National Bureau of Standards. pp. 16–141.
  2. Claerbout, Jon; Muir, Francis (2008). "Hubbert math" (PDF). Retrieved 2020-06-08.
  3. Roper, David. "Where Have All the Metals Gone?" (PDF). Archived from the original (PDF) on 2007-09-28.
  4. Roper, David. "Projection of World Population". Archived from the original on 2007-02-18.
  5. Deffeyes, Kenneth (February 24, 2005). Beyond Oil - The view from Hubbert's peak. Hill and Wang. ISBN   978-0-8090-2956-3.
  6. Andruleit; Franke; Ladage; Lutz; Pein; Rebscher; Schauer; Schmidt; von Goerne (August 2019). BGR Energy Study 2018 (PDF). Data and Developments Concerning German and Global Energy Supplies. Vol. 22. Hannover: Federal Institute for Geosciences and Natural Resources (BGR). p. 43. Although global conventional crude oil production has stagnated since 2005, it still retains a share of around 79 % of total production, and will therefore continue to play a most significant role in the long term in supplying liquid hydrocarbons (Fig. 3-4).
  7. Sam Foucher (2006-08-18). "A Different Way to Perform the Hubbert Linearization". The Oil Drum . Retrieved 2020-06-08.
  8. Canogar, Roberto (2006-09-06). "The Hubbert Parabola". GraphOilogy. Archived from the original on 2011-07-08. Retrieved 2007-03-07.
  9. Rutledge, David (2011-01-01). "Estimating long-term world coal production with logit and probit transforms". International Journal of Coal Geology. 85 (1). Elsevier: 23–33. Bibcode:2011IJCG...85...23R. doi: 10.1016/j.coal.2010.10.012 .