Four US Linearizations

EIA Field production of crude, and four Hubbert models based on different linearizations. Source: EIA for the data, models as described in the text.

I'd like to begin an investigation into how much error we might expect to make by extrapolating Hubbert linearizations forward (for background on the method, try this post). It's going to take several posts, I think; there's a bit of work to do, and I'm not sure quite what schedule I'll be able to accomplish on it. Anyway, let's start gently with a look into the situation with the US, which is a popular target for linearization, because it's one of the most mature regions (though perhaps Romania is the best of all). Anyway, the EIA happily has US production data back to 1859, and to that I added one point for 2005, by summing the monthly data, and multiplying by 12/10. Note that this data includes Alaska - I prefer to work with whole countries.

To kick us off, I started with a regression of the data between 1958 (which is the point Deffeyes chooses) and today. That gives me this picture (Model 1):

Linearization plot of EIA Field production of crude, together with linear fit for the data from 1958 to 2005. Green data is used for line fitting, plum data is not. Source: EIA for the data. 2005 data is extrapolated from data through October.

Clearly, we got a pretty good fit on recent data, and we ended up with a URR (ultimately recoverable resource) estimate of about 230gb, and a K of 5.7%. Remember that K is best thought of as the initial growth rate in the model, or equivalently, the final decline rate (with change rates gradually interpolating from one to the other over the history).

Now, recall that there's this slightly black art in the linearization process. They hardly ever fit well early on in the history, so one has to choose some point from which to start fitting. Obviously, choosing different points gives slightly different results. How sensitive are we to that? Well, if we picked 1930, which is about the earliest vaguely plausible place to try and start a straight line, we'd end up looking like this picture, which I'm calling Model 2:

Linearization plot of EIA Field production of crude, together with linear fit for the data from 1930 to 2005. Green data is used for line fitting, plum data is not. Source: EIA for the data. 2005 data is extrapolated from data through October.

Clearly, this doesn't do quite as well on the recent data. Our URR has changed to about 220gb, and our K is now 6.2%. So the URR changed by around 5%, and K changed by around 10% of it's value. Those are probably decent rough estimates of our uncertainty in these quantities with the full history known to us now.

But how would we have done in the past? Based on model 1, the 50% point would have passed in about 1976. So let's try and do a linearization analysis if we only have data up to 1976 (but sticking with the 1930 starting point). In the world case, this is roughly analagous to trying to extrapolate linearizations today. That gets us Model 3:

Linearization plot of EIA Field production of crude, together with linear fit for the data from 1930 to 1976. Green data is used for line fitting, plum data is not. Source: EIA for the data. 2005 data is extrapolated from data through October.

Clearly, we've dropped our URR down to about 190gb - off by a little more than 15% from a current estimate. Our K estimate is now 6.4% - about 11% higher than we would estimate today. Finally, let's suppose we were doing this back in 1955, and trying to scoop Hubbert. I would have come up with something like Model 4 here:

Linearization plot of EIA Field production of crude, together with linear fit for the data from 1930 to 1955. Green data is used for line fitting, plum data is not. Source: EIA for the data. 2005 data is extrapolated from data through October.

Our estimate of URR is now down to about 180gb - over 20% low. Our estimate of K has gone up to 6.6% (and truth be told, we're lucky to do that well - there's a lot of countries where we'd go hopelessly wrong trying to extrapolate from this early on in the history).

If we translate these four models back into production versus time graphs, we get the following picture:

EIA Field production of crude, and four Hubbert models based on different linearizations. Source: EIA for the data, models as described in the text.

Again, Model 1 is our best current fit. It is indeed a remarkably good fit for such a very simple model applied to what would appear to be a very complex process (I confess that I still don't fully understand why this model works as well as it seems to in practice). Model 2 is the one where we started fitting in 1930 instead of 1958. It doesn't work quite as well, but it's not bad either.

Model 3 is if we only have data to 1976. You can see pretty clearly where it goes wrong: it misses the Alaska North Slope production. Of course, if we had really been doing this in 1976, we would have known that there was a super-giant field about to come on stream much later than the others and change the picture. And of course model 4 is a little worse again - that's the one where we only have the history to 1955.

The conclusion in the US seems to be this: doing the linearization in the early days causes you to underestimate URR and overestimate K. However, it's by no means useless. The errors even in 1956 are much less than a factor of 2 - it gets us in the right ballpark. And I think it's particularly striking that even when you get the URR somewhat wrong, the shape of the downslope is still about right. Very early on in the history, you could tell with about a 10% relative error what the post peak decline rates were going to be. That's remarkable.

What I'd like to do next is write a little program to more systematically estimate the uncertainties in the linearization. I think we can also make some improvement by underweighting the early data and overweighting the later part of the history. Then I want to get into people's concerns about whether the world is fundamentally different enough from individual countries that we cannot safely use the same process. But those will have to wait for future posts.

What if it were to be found that ANWR contained 60 Billion barrels instead of 12, how would this make a theoretical model 5 look? I only say this to play devil's advocate. Let's also assume that Exxon and Halliburton could pull the oil out at ultra-efficient, western speeds, with many, many drill rigs and a new pipeline straight into the heart of Detroit. C'mon Stuart, I know you can put one more curve on the graph :)
You're right that you don't know if there's going to be advances in tech, or new oil finds.

But that would only gives us an offset of a couple of years before the peak hits. The downslope would be the same after that. I think Stuart shows us that here.

Actually, the figures for ANWR are 4.3 and 11.8 Gb (95% and 5% probabilities). This isn't going to change the analysis much.

If it contained 60 Gb, Detroit would have been receiving that oil for many years now... ;)
The people making these estimates are the USGS.  With much pardon to my friends in the USGS, I do not believe them for a second.  The USGS is not engaged in finding and developing oil and gas.  Oil companies are.  (I used to work for the USGS so I know a little bit about their capabilities or lack thereof.)

A reasonable P50 estimate for ANWR might be 2 billion barrels. The structures are small and complex.  12 billion is out of the question, and 60 billion is for people who don't know what a rock is.  No one in the industry is predicting finding another Prudhoe Bay or Kuparak River field in ANWR.

And of course no one has found a drop of oil in ANWR.  There could be nothing there worth trying to go after.  If it ever opened up, the industry might drill 5 to 10 wildcats, find some small volumes, and walk away.

This is a crucial point. I recall learning in grad school that the USGS has only drilled one hole (or not many more) and that the results of the tests were not public. The probability assessments noted above seem to be based on statistcal modeling using geological features:

"An evaluation was made of each of 10 petroleum plays (similar geologic settings). For each play, USGS constructed statistical distributions of the number and size of potential accumulations based on a probabilistic range of geologic attributes."

http://www.eia.doe.gov/pub/oil_gas/petroleum/analysis_publications/arctic_national_wildlife_refuge/h tml/analysisdiscussion.html

Does anyone know any more about the amount of information we have on ANWR resources and hopw reliable they are?

This is perhaps more to do with the analysis on global depletions rates using the same technique but it does apply to  this USA-only analysis.

The "softness" of your predicted global depletion rates compares to what you have shown here for the USA. My problem with global versus country-specific depletion rates is the fact that when the USA peaked, oil I suspect became increasingly more expensive to extract from US wells. So, would it not be the case that American oil consumers would buy up more cheaper and foreign oil even if domestic oil production could be increased? The result being companies profits drop and wells are capped or throttled?

This would have the effect of "softening" the USA depletion curve due to cheaper alternatives abroad.

Now for the global situation, at a global peak, there is no cheap seller of last resort to go to. The world cannot go off-planet to obtain cheaper crude to soften the depletion curve like the USA did.

You only have two options, pump faster and accelerate the depletion or cut back in terms of recession or alternatives.

The truth as ever will be in between these.

That would seem therefore to suggest depletion rates will be more like 7%+. It seems in the case of recently peaked Britain they have decided just to pump and pump or perhaps Peak Oil has already forced them too?

Britain's capacity is all off-shore, and the platforms require a minimum production rate to justify keeping them open.  Pumping like mad might be the strategy for maximizing ROI (they run out sooner, but they can retire the expensive platforms sooner too).
But the earlier models weren't really wrong, they just did not account for the additional North Slope oil - which of course they could not.  So Hubbert only models the "natural" behavior of fields, or aggregations of wells/fields.  But if something else happens to monkey with them, it breaks down.  This could be geopolitical, natural disaster, new technologies, or new discoveries.  So the accuracy when applied to world oil production will be dependant on having a handle on new oil discoveries and all those other factors.  OTOH, when you aggregate things, I bet you'll mostly see the effects of the few really big fields, and the effects of exceptional situations will tend to average out.  If we had the data, we could tell where things were tracking.

I bet on the whole Hubbert will be pretty true as long as we account for new fields coming on line, with my biggest concern being whether techniques used to extend the production volume have significantly changed the shape of the curve.

But we do have a good grasp of how much is left to be discovered.  The discovery data is well ahead of production.  If we tie the discovery data into production, this should be factorable in as a variance percentage.  

Actually to a degree the Hubbert model does try and account for oil that is yet to be discovered.

In Deffeyes 2005 book "Beyond Oil", he applies Hubbert linearization not only to the oil we have already used, but also to two others.

The second plot is the amount of oil we have used plus the known reserves (removing bogus reserves from OPEC nations).

The third plot is 'hits'.  Essentially a measure of oilfield discoveries, and this is an attempt to measure the amount of oil in oilfields yet to be discovered.

You can then fit straight lines to all 3 of these, and he has such graphs in one of his books.  In the summary, he has the parameters for the lines that he fit to the things.  All 3 have a Qt of 2.013 trillion barrels.

The first curve (oil that we have consumed already) has a=5.9%, peak 2005, and percentage used is 49%.

The second curve (oil that we have consumed plus reserves) has a=7.2%, peak 1978, and percentage of total of 82%.

The third curve (hits - oilfields that we know about) has a=8.1%, peak in 1964, and percentage of total of 94%.

The third curve implies that only 6% of the total world oil is in oilfields that we don't know about yet.

I have to admit that I am still a little fuzzy on the concept of 'hits', and what it means.

Twilight,
Hubbert works well for USA oil, but by your definition, it is not entirely "natural".  For geopolitical factors, many areas were placed out-of-bounds for environmental reasons. There were many innovations integrated into the production data.  Even with the major recent discovery, the North Slope, the older estimates project forward fairly well.  

For USA oil, Hubbert seems to cope well with most of the "unnatural" factors that you mention.   Going global, politics will probably throttle production more than disasters.  New discoveries probably won't be as significant for the world as for the USA.  It is likely that we will need some stupendous additional innovation just to match the innovation embedded in the USA models. This reasoning implies less global URR "loft" and global depletion rate "depression" than we see with the USA models.

Yes, I guess that's why I suspect that the Hubbert model will work, because even in this example, where we did indeed have all of those issues, it still works pretty well.  That tells me that you'd really have to "monkey" with things a lot to get to the point where the errors are really big.  I think it's just that the effects of the really big fields swamp out everything else - combining with smaller fields of similar shape just shifts the curve a bit.  As goes Ghawar, so goes the world.
Stuart, you really spend a good bit of time with excellent graphs and charts.  I believe you come up with some great results.  Linearization works wonderfully well when depletion rates behave according to Hubbert's models.  That is, normal drilling and pumping.  When you factor in huge gas-water injection variables along with horizontal and bottle-brush drilling (especially from super large fields), the results using linearization becomes fuzzy.

I would like to see you try linearization on the North Sea Fields....and see how they behave.  Collapse rates that will occur in Mexico's largest field, as well as the Saudis, the Kuwaitis, and the Russians will be quite interesting to see if they behave nicely to linearization.

We have looked at the North Sea, and the linearization method does predict the high decline rates. You can see the UK graph in this post here, where linearization calls for K=13% (once it settles down after the dual peak). Indeed we have seen very high decline rates of 6-15% in recent years in the UK. Furthermore, it calls for very high decline rates in Norway also:

width="90%">

OFF TOPIC

BUT.....

http://news.ft.com/cms/s/f39fa8e4-7e25-11da-8ef9-0000779e2340.html

I think this is the first time China has overtly expressed a desire to move away from the dollar. Please correct me if it isn't. The dollar is already moving lower - if this continues into the immediate future, it doesn't bode well for international trade.

My company is watching this closely, as dollar valuation plays heavily into our international exploration budgets.

GeoPoet, the Chinese have been sending up some less official sounding trial balloons for at least several weeks, but you're right, this is a new one on me and that is why, for the first time in my life, I invested in options on gold futures when the price of gold dipped in December. In December economist Hu Yongding, who sits on the monetary policy committee of the People's Bank of China, was quoted in the English-Language Hong Kong Standard, saying this very thing. In the same article, however, a more senior source "corrected" Mr. Hu. Was this "damage control," or was it a way of distancing a trial balloon from "official policy?" I'm no Sinologist, I have no opinion on that one.

And last week semiofficial sources started talking about increasing the percentage of gold in China's official monetary reserves. You can follow these sorts of things on the Kitco.com website, which is dedicated to information about investing in gold.

I wouldn't say this was off topic Geopoet. Such a shift by the Chinese could initiate the unwinding of huge global financial imbalances, which would have enormous implications for energy, particularly for the US. If the dollar were to lose its status as global reserve currency, oil would become much more expensive for Americans who have been insulated from currency risk as an aspect of oil price volatility. The US would be unlikely to take it lying down, which would usher in a very dangerous era indeed. Perhaps China has been speaking to Iran and Russia about pricing oil in euros?
This unwinding is all too inevitable, I am afraid. Financial systems are stable for long periods -- and suddenly they are not. On any given day, the markets focus on one number -- the employment report, what the Fed did, GNP, Microsoft earnings -- and there is not always a predictable reason why number x is the number du jour. Commentators will talk as if they know why the market focussed on the number that turned out to be the number du jour, but, trust me, we are all cluelesss. If you are a trader like me, you just go with the flow.

I bring this up to say: for several years the Chinese have accepted our essentially worthless dollars in order to keep their export-led growth going. They will keep doing that until they decide, that they won't. And when the markets wake up to that (before or after the fact, I can't say) everything will change. No one can say how it will look, but it won't be pretty. For some interesting (if economically heterodox) thinking on the subject, try this: http://www.business-in-asia.com/dollar_crisis.html.

I think of the number du jour as a post hoc rationalization for a gut feelings, in other words window dressing for unconscious herding behaviour. Most trading seems to have little to do with rational choices based on real information, and a lot to do with instinctive emotional responses to the evidence of others' actions.

The whole system seems to feed off itself, producing its own internal dynamic which is only loosely connected to reality. The inevitable result is swings of positive feedback at various degrees of trend, so to speak. There are those who see in this the application of fractals to market behaviour. Robert Prechter is the one I am most familiar with, but I also remember an article to that effect by Mandelbrot (I'll try to find the link). Its an intriguing concept which I follow out of interest, but am not intellectually wedded to. The psychological argument seems more compelling to me than the potential mathematical applications, but then my approach to life is generally more qualitative than quantitative.

Here's a link to the article I remember by Mandelbrot: http://www.elliottwave.com/education/SciAmerican/Mandelbrot_Article2.htm
Sunlight, I just read through your link (which needs to have an extra period removed from the end of the URL in order to work). Note to the-powers-that-be at TOD - this should be a separate topic as it addresses the important aspects of inflation versus deflation and global financial imbalance. The interviewee is Richard Duncan (not the same Richard Duncan associated with the Olduvai theory), whose book The Dollar Crisis I would recommend very highly.

Here's an excerpt from the interview Sunlight links to:

"Question 6:  Also in the book you note that "over-investment causes excess capacity and excess capacity causes deflation."  Could you explain to our readers exactly what deflation is, why it is to be feared and why you believe the threat of global deflation is real?

Answer:   Inflation comes about when there is too much demand relative to supply.  Deflation is caused when there is too much supply relative to demand--or, more precisely, purchasing power.  When credit is abundant, companies borrow and expand industrial capacity and aggregate supply increases.  However, the purchasing power of the public does not necessarily expand in line with supply.  During periods of very rapid credit expansion as in the 1970s and `80s in Japan, in the Asia Crisis countries during the 1990s and in China during the last 15 years, supply grows much faster than the personal income of the public.  When supply exceeds purchasing power, prices fall.  That is why Japan and China are suffering from deflation now.  The Asia Crisis countries avoided deflation by devaluing their currencies and exporting deflation abroad.  For example, the devaluation of the South Korean Won after the Asia Crisis contributed to the downward pressure on global semiconductor and steel prices.

The US current account deficit is flooding the global economy with dollar liquidity.  When the dollar earnings of the surplus nations are deposited into their domestic banking systems, those dollars, being exogenous to those banking system, act as high powered money and spark off an explosion of credit creation.  Excessive credit creation permits over-investment, which, in turn, causes excess capacity and deflation.  So long as the huge US current account deficits continue to flood the world with dollars, global deflationary pressures are very likely to continue to build, as reckless credit creation results in more industrial capacity than can be absorbed at the prevailing price level.

The reason that deflation is harmful is because it undermines corporate profitability and, thereby, leads to rising unemployment and falling purchasing power."

agreed. Very good link .I would see this discussed.Thanks.
I am in the process of preparing a presentation on this material to a discussion group of my fellow traders. I would also be happy to make it available to TOD if you think that would be helpful. It might take a bit longer. As a reader of "The Dollar Crisis" you, however, will know most of what I have to say on the topic.
Personally, I think it would be very helpful. We typically focus on supply here and voluntary or planned demand-side management, but I think this material helps to shed light on the nature of involuntary/unplanned demand reductions. The unwinding of these huge financial imbalances as the credit bubble implodes will have an enormous impact on purchasing power. If purchasing power falls off a cliff then demand, even for energy in cases where there is a structural dependency on it, falls, and price falls with it for the duration of the deflationary period.

Demand is not what we would all like to use, it's what we can afford to use, and demand will plummet with havoc in the global financial system. Economic activity would be further constrained by the very high real interest rates that result when inflation is negative, further depressing demand (a positive feedback loop). Nominal interest rates cannot go lower than zero, but zero is not low enough when the real rate of interest is the nominal rate minus negative inflation. Japan has had this problem for many years.

I think Richard Duncan's work helps to bridge the gap between the economics camp and the energy camp in order to build a realistic scenario for the future. Whether we will have inflation or deflation is a crucial distinction when it comes to planning for the future as the strategies one would employ are very different. Personally, I am convinced that we are headed towards an extremely severe global economic depression. This topic of discussion, and Sunlight's link to the Richard Duncan interview, would make a good separate thread.

absolutely.  please send it along to the editors address...
Thanks Prof G.
and I'll make an open thread...and put that in there in the first comment.
maybe off topic....but HOLY MOLY...i'm aware of the previous trial ballons, but now the cat is out of the bag. i was waiting for the chinese to start asserting ,what i'm sure they believe, is their "proper" politico-economic place in the world....and that time is now...they need commodities and they need oil to keep their feverish economy going, viz.
"It is a subtle but clear signal that they are interested in moving away from the US dollar into other currencies, and are interested in setting up some kind of strategic commodity fund, maybe just for oil, but maybe for other commodities," he said.

...he did indeed.

Very interesting! I think it is very important to try to assess the quality of the estimates resulting from the Hubbert linearization if we want to trust its predictive power.

If I may, I computed the URR fluctuations for a finer and larger range of years (assuming a minimum of 20 years between the starting year and the ending year).

We can see clearly that the ending year has a larger impact on the quality of the estimate. From 1989 up to 2005, the estimate is quite stable (between 200 and 240 Gb) regardless of the starting year.

The same for K which is quite stable:

Question: How do you find the peak year date? it's not directly available from the Hubbert linearization fit and requires a second fit (usually on the production curve).

Another intresting aspect we should look into is the impact of the discovery distribution on the Hubbert linearization. Some countries are problematic such as Nigeria for example.

Now THAT is a cool analysis. Can you do it for the World as a whole?
Here we go:



The resulting RMS error between the data and each curve:

Very nice! That's saved me some work.

To find the peak (at least of the Hubbert curve), just divide your URR estimate by 2 and then find the year where Q = URR/2.

Ok, so you take the time when the midpoint of production is reached. But in the case of model 4, you're supposed to have the data only up to 1955 far from the production midpoint!
You're right - so in those cases there's no choice but to fit in the P/t or Q/t domain (what I outlined doing is essentially a simple estimator in the Q/t domain in any case).
Very interesting! I think it is very important to try to assess the quality of the estimates resulting from the Hubbert linearization if we want to trust its predictive power.

The quality of the estimates will only be known when the actual data are gathered.  The URR estimate graph is a plot of a hypothesized function.

 Thanks for your work.It seems to me with my limited mathematical perspective,  the key to Hubbert's model is the large wells size relative to later wells/discoveries. So as I understand this, it likely throws  off Models 3 &4.  So other large well/discoveries coming on line would be one error , plus the number of large wells that use water injection would be another error to expect.
The other errors are economic & social, the latter factor being the one I expect to be too difficult to more than list possibilities,yet possibly the most significant.
An analysis Im working on is taking the world (hubbert) graph of yearly production (currently 30B Barrels) and ESTIMATING going forward what % of that 30B was needed to produce the energy. Also, I am putting in a base amount of oil consupmtion that allows the world economy to function properly without cease in trade, starting of wars etc. Obviously these will be pretty subjective, but the objective is to create 3 graphs - a) the gross world production profile (annually) b) the net production available for non-energy consumption and c)the above two graphs shown AFTER a baseline fixed amount of oil is consumed to keep the world system flowing.

The presumption being that seemingly benign decline rates of 4%, when viewed from net energy, and fixed vs marginal usage may be ginormous.

Hubbert's model I guess is only barrls, independant of net energy, so yes;with no big discoveries this factor can't be compensated for as when U.S. peaked..  Also if we knew how many of the big wells , accounting for their production, and if water injection would tell us some about decline %.
It seems that there is a relationship between K and the URR. When the K estimate increases, the URR has a tendency to be underestimated (model 3 and 4). Conversly, when K gets lower the URR gets higher (model 1 and 2). K controls the steepness of the curve and the URR the area under the curve. In summary, a flat curve leads to a higher URR and a sharp narrow curve leads to a small URR. This effect is not surprising when you notice that there is a mathematical relationship between these variables at the point of maximum production:

Pmax= k x URR/4
 

Pmax being the m