Peak Oil News

 

  Login or Register
 
Menu
 News
 Search
 Topics
 Stories Archive
 Submit News
 Discussions
 Code of Conduct
 Forums
 Forums Search
 Last 24 Hours
 PO 24hrs
 Peak Blog
 Resources
 About Us
 Downloads
 Web Links
 PeakWiki
 PeakPortal
 Focus Search
 Peak TV
 Peak Oil Boston
 Members
 Your Account
 Members List
 Ignore List
 JOIN!
 Private Messages
 
google
 
PeakSpeak
NICKNAME

Download TeamSpeak
What is PeakSpeak?
Peak Oil on IRC
 
Photo Album
Submit Photo
Peakoil.com is You!


member photos
 
Light Sweet Crude Oil
 
Member Quotes
I want my mommy!

Buggy

Suggest Quote

 
aspo08
 
ICM
Cisco & Net App Training
 
Peak Oil News: Forums

Peakoil.com :: View topic - Peak prediction based on the Riccati equation: 2007!
 Forum FAQForum FAQ   SearchSearch   UsergroupsUsergroups   ProfileProfile   Log in to check your private messagesLog in to check your private messages   Log inLog in 

Peak prediction based on the Riccati equation: 2007!
Goto page 1, 2, 3  Next
 
Post new topic   Reply to topic   Printer-friendly version    Peakoil.com Forum Index -> Depletion Modeling
View previous topic :: View next topic  
Author Message
khebab
Moderator
Moderator


Joined: Sep 27, 2004
Posts: 935
Location: Canada

PostPosted: Fri Jun 24, 2005 6:37 am    Post subject: Peak prediction based on the Riccati equation: 2007! Add User to Ignore List Reply with quote

A very interesting scientific paper that is supposed to come out this year:

GUSEO, R., DALLA VALLE, A. and GUIDOLIN, M. (2005). World Oil Depletion Models: Price Effects Compared with Strategic or Technological Interventions

An earlier similar paper is also available:

GUSEO, R. e DALLA VALLE, A. (2004). Oil and Gas Depletion: Diffusion Models and Forecasting under Strategic Intervention, Atti della XLII Riunione Scientifica della Società Italiana di Statistica, Bari 9-11/6/2004, Vol. Sessioni Spontanee, 733-36, CLEUP, Padova.

Most of us are familiar with the curve fitting techniques used by Laherrère, Hubbert, etc. in order to predict the position of the peak of oil production. This approach has been widely criticized in particular by Michael Lynch because of the apparent arabitrary curve model employed and the large uncertainty about the URR (Ultimate Recoverable Ressource):
Quote:
Obviously, pure Hubbert models have no room for economics since the basic differential equation is the classical logistic. The under estimation of Hubbert model is a known aspect but no analysis is given in order to introduce more flexible tools capable, for instance, to track an intervention or a systematic perturbation. Lynch supposes that the origin of the problem relies on the strong assumption of a static URR (Ultimate Recoverable Resource). The impact of technology in expanding resources is sometimes reasonable but a lot of systematic interventions may affect normal evolution. For instance, the forced production of PrudhoeBay in Alaska was a positive shock, or, on the contrary, the negative effect ofaccidents in pipelines (Piper Alpha disaster) in the North Sea depressed regular production. These accidents led to examine restructuring of UK offshore safety legislation followed by major changes to the Petroleum Reven Tax. These radical interventions are not recognized by a simple Hubbert model


Quote:
Evidence is given to a common effect of that shocks, i.e., a transition to low-energy tools (computers, electronics, robotics, bio-engineering, etc.) as compared with traditional industrial products, steel-car, chemicals, etc.. The role of the "well-informed spectators" of the oil game, i.e., journals, specialized consultants, institutes, is important as they "massage the message" and, clearly, influence oil price-settings. This result contributes to enhance our work which highlights the surprising weak relationship between observed prices, oil demand and supply. The decade post 1986 is relatively stable with some exceptions that confirm paradoxical behaviours of the quantity-price relationship.


In this paper oil production is considered as a diffusion model based on a complex differential equation called generalized Bass models (GBM). Solving this equation leads to the famous Riccati formula. This approach do not require the knowledge of the URR and is able to model price shocks like for instance the 70's oil shock!

Quote:
In the sequel we may consider oil as a production directly controlled by correlated diffusion processes, i.e., individual mobility, heating, electricity production, chemical applications, etc., so that crude oil production may be thought as a parallel diffusion process. Such production is influenced by regulatory or strategic interventions of political and economical nature which are usually
persistent and conceptually deterministic.


What's amazing with this new model is that shocks are explicitely modeled:
Quote:
The subtle and shifty deviation originated in the post{war period beyond 1951 has generated an incremental consumption of oil not fully balanced by the feedbacks of 1970's. This deviation from the Hubbert model (exactly centered on the world oil data until 1950) is recovered by the GBM model with three shocks with an abrupt contraction of the diffusion process depletion component. The contraction is quite severe. The present model forecasts an URR of 1524 Gbo, a peak positioned in 2007 with a production of 76.34 mb/d, a depletion of 90% URR in 2019 with a production of 55.33 mb/d and a 95% saturation during 2023 with a production of 36.13 mb/d.


The prediction results are quite pessimistic:

Quote:
The subtle and shifty deviation originated in the post{war period beyond 1951 has generated an incremental consumption of oil not fully balanced by the feedbacks of 1970's. This deviation from the Hubbert model (exactly centered on the world oil data until 1950) is recovered by the GBM model with three shocks with an abrupt contraction of the diffusion process depletion component. The contraction is quite severe. The present model forecasts an URR of 1524 Gbo, a peak positioned in 2007 with a production of 76.34 mb/d, a depletion of 90% URR in 2019 with a production of 55.33 mb/d and a 95% saturation during 2023 with a production of 36.13 mb/d.



_________________
______________________________________
http://GraphOilogy.blogspot.com


Last edited by khebab on Fri Jun 24, 2005 10:37 am; edited 1 time in total
Back to top
View user's profile Send private message Visit poster's website
Jack
Dark Lord
Dark Lord


Joined: Aug 11, 2004
Posts: 5085

PostPosted: Fri Jun 24, 2005 6:42 am    Post subject: Add User to Ignore List Reply with quote

Outstanding! Thanks for posting this.
_________________
Dieoff. Fun to watch. Better with hot buttered popcorn! new_popcornsmiley
Back to top
View user's profile Send private message
Barbara
Light Sweet Crude
Light Sweet Crude


Joined: May 26, 2004
Posts: 1190
Location: Zoorope

PostPosted: Fri Jun 24, 2005 6:47 am    Post subject: Add User to Ignore List Reply with quote

Never heard of this on Aspoitalia.
I post it there, let's see what they say.

Thanks, and I'll keep you updated! Smile
_________________
**no english mothertongue**
--------
Objects in the rear view mirror
are closer than they appear.
Back to top
View user's profile Send private message
khebab
Moderator
Moderator


Joined: Sep 27, 2004
Posts: 935
Location: Canada

PostPosted: Fri Jun 24, 2005 7:05 am    Post subject: Add User to Ignore List Reply with quote

It's a very exciting work that use a mathematical economic model that should please Michael Lynch Smile . Note that this paper has been reviewed by Campbell.
_________________
______________________________________
http://GraphOilogy.blogspot.com
Back to top
View user's profile Send private message Visit poster's website
pup55
Expert
Expert


Joined: May 26, 2004
Posts: 3692

PostPosted: Fri Jun 24, 2005 7:45 am    Post subject: Add User to Ignore List Reply with quote

Interesting that he considers the demand spike in the 50's to be a "shock". I wish they had explained the variables in equations (1) and (3) a little better, so I guess we will have to look up the references for more detail.

Also, his fit of the verhulst curve to the historical data is a little different from ours. It looked like he minimized his error prior to 1950, with the result that the "shock" of the 60's of really high production and the subsequent "corrections" are all well above the verhulst curve.

Excellent post, Khebab.
Back to top
View user's profile Send private message
khebab
Moderator
Moderator


Joined: Sep 27, 2004
Posts: 935
Location: Canada

PostPosted: Fri Jun 24, 2005 8:04 am    Post subject: Add User to Ignore List Reply with quote

pup55 wrote:
Interesting that he considers the demand spike in the 50's to be a "shock". I wish they had explained the variables in equations (1) and (3) a little better, so I guess we will have to look up the references for more detail.

I agree, the variable a not very well defined and quite esoteric!
Code:

m:  potential market or carrying capacity
z'(t):   instantaneous adoptions (daily world production?)
D(t):   oil price
x(t):  intervention fonction or common intervention tool (function of oil price?)

_________________
______________________________________
http://GraphOilogy.blogspot.com
Back to top
View user's profile Send private message Visit poster's website
MicroHydro
Light Sweet Crude
Light Sweet Crude


Joined: Apr 10, 2005
Posts: 1286

PostPosted: Fri Jun 24, 2005 9:50 am    Post subject: Add User to Ignore List Reply with quote

Their "2007 peak" is about 8 millon bpd below current production. What is up with that?
_________________
"The world is changed... I feel it in the water... I feel it in the earth... I smell it in the air... Much that once was, is lost..." - Galadriel
Back to top
View user's profile Send private message
killJOY
Fission
Fission


Joined: Feb 21, 2005
Posts: 2478
Location: ^NNE^

PostPosted: Fri Jun 24, 2005 10:11 am    Post subject: Add User to Ignore List Reply with quote

People that try to pass off sentences like this:
Quote:
Evidence is given to a common effect of that shocks, i.e., a transition to low-energy tools (computers, electronics, robotics, bio{engineering, etc.) as compared with traditional industrial products, steel{car, chemicals, etc.. The role of the "well-informed spectators" of the oil game, i.e., journals, specialized consultants, institutes, is important as they "massage the message" and, clearly, influence oil price-settings
don't deserve the mf-ing time of day.
_________________
"By the time individuals discover that remaining resources will not be adequate for the next generation, the next generation has already been born. " David Price
Back to top
View user's profile Send private message
khebab
Moderator
Moderator


Joined: Sep 27, 2004
Posts: 935
Location: Canada

PostPosted: Fri Jun 24, 2005 10:35 am    Post subject: Add User to Ignore List Reply with quote

MicroHydro wrote:
Their "2007 peak" is about 8 millon bpd below current production. What is up with that?

That's a good question, probably because they used the 2003 BP energy review their prediction may be a little bit too pessimistic.
_________________
______________________________________
http://GraphOilogy.blogspot.com
Back to top
View user's profile Send private message Visit poster's website
Dan1195
Heavy Crude
Heavy Crude


Joined: Mar 19, 2005
Posts: 259

PostPosted: Sat Jun 25, 2005 5:38 pm    Post subject: Add User to Ignore List Reply with quote

Why do they have a URR of ~1500 gbo? Someone explain why this is not too low?
Back to top
View user's profile Send private message
khebab
Moderator
Moderator


Joined: Sep 27, 2004
Posts: 935
Location: Canada

PostPosted: Sat Jun 25, 2005 6:22 pm    Post subject: Add User to Ignore List Reply with quote

Dan1195 wrote:
Why do they have a URR of ~1500 gbo? Someone explain why this is not too low?

From what I understand the URR is estimated from the non linear fit between the Riccati model and the actual data, it's not a prior assumption.
_________________
______________________________________
http://GraphOilogy.blogspot.com
Back to top
View user's profile Send private message Visit poster's website
ohanian
Light Sweet Crude
Light Sweet Crude


Joined: Oct 17, 2004
Posts: 1087

PostPosted: Sun Jun 26, 2005 7:44 pm    Post subject: GBM tutorial Add User to Ignore List Reply with quote

tutorial

http://home.uchicago.edu/~zhuwang/gsb37201/tutorial1.pdf
Back to top
View user's profile Send private message
khebab
Moderator
Moderator


Joined: Sep 27, 2004
Posts: 935
Location: Canada

PostPosted: Sun Jun 26, 2005 8:38 pm    Post subject: Re: GBM tutorial Add User to Ignore List Reply with quote

ohanian wrote:
tutorial

http://home.uchicago.edu/~zhuwang/gsb37201/tutorial1.pdf

Ohanian, Thank you for the link!

There is also this document:

GUSEO, R. (2004). Interventi strategici e aspetti competitivi nel ciclo di vita di innovazioni; Strategic Interventions and Competitive Aspects in Innovation Life Cycle ; Working Paper Series, N. 11, Department of Statistical Sciences, University of Padua, Italy;
Too bad it's in italian but it gives some additional mathematical details.
_________________
______________________________________
http://GraphOilogy.blogspot.com
Back to top
View user's profile Send private message Visit poster's website
khebab
Moderator
Moderator


Joined: Sep 27, 2004
Posts: 935
Location: Canada

PostPosted: Sun Jun 26, 2005 11:13 pm    Post subject: Add User to Ignore List Reply with quote

Ok, I managed to reproduce the article results. Apparently, the parameter estimation is performed by a nonlinear least square algorithm (ex: Marquardt). The parameter values proposed by the author are the following (see Table 3 p. 27):
Code:
m = 1524;
p = 0.00010439; q = 0.063497;
c1 = -0.3021860; b1 = 0.05674; a1 = 80.50;
c2 = 0.0717753; b2 = 0.07187; a2 = 51.07;
c3 = -0.2272032; b3 = 0.07098; a3 = 74.60;

Wich leads to the following curve:


Figure 2

The data used here are from the 2005 BP review provided by pup55 (see thread Updated Verhulst model). We can observe that there is a significant discrepancy between the fit and the data between 1900 and 1965 which is not present on the figure in the article. Obviously, there is a problem in the data, maybe pup55 has the answer.

The oil production noted z'(t) is given by the following GBM differential equation:
Code:
z'(t)= (p + qz(t)/m)(m - z(t))x(t)

where z(t) is the cumulative oil production given by the Riccati equation. m is the equivalent of the URR so (m - z(t)) is the quantity of oil left at the time t (proven reserve). x(t) is called the control model and has a great influence on the resulting curve. In particular, x(t) is used to model discontinuities in the curve like the 70s oil shock. If x(t)= 1, the GBM becomes the simpler Bass Model. In the GBM, x(t) is used to model the effect of pricing and investment (see the link provided by Ohanian above for more details).

I used the following control model (3 exponential shocks) noted gbm3e in the article:
Code:

x(t)= 1 + E1 + E2 + E3
with Ei=ci x exp(bi x (t - ai))
the signs of ci and bi are important and influence greatly the shape of the curve.

I performed a new least square estimation using the Marquardt algorithm. I excluded the data from 1900 to 1965 because of the unexplained poor fit for these points:
Code:
F(x): 2.209514e-006
Nb of Iterations: 500
Parameters:
URR=   0.84 p= 0.000021 q= 0.096700
Control model
Shock 1: c1= -0.163277 b1= -0.015254 q1= 80.484536
Shock 2: c2= 0.001622 b2= 0.125247 q2= 51.070223
Shock 3: c3= -0.125364 b3= 0.061326 q3= 75.000000

The new URR went down to 0.84 * 1524= 1280 Gb, and the peak production date is now in 2005:


Figure 3

Other control model functions are possible in particular prices fluctuations can be injected into the model.
_________________
______________________________________
http://GraphOilogy.blogspot.com
Back to top
View user's profile Send private message Visit poster's website
ohanian
Light Sweet Crude
Light Sweet Crude


Joined: Oct 17, 2004
Posts: 1087

PostPosted: Mon Jun 27, 2005 4:58 am    Post subject: Add User to Ignore List Reply with quote

eh? What happens when m=2500 ?

When will it peak then?
Back to top
View user's profile Send private message
Display posts from previous:   
Post new topic   Reply to topic   Printer-friendly version    Peakoil.com Forum Index -> Depletion Modeling All times are GMT - 6 Hours
Goto page 1, 2, 3  Next
Page 1 of 3

 
Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum

Atom News FeedRSS 1.0 News FeedRSS 2.0 News FeedRSS Forums Feed