ISRN Mechanical Engineering

Volume 2014 (2014), Article ID 529748, 37 pages

http://dx.doi.org/10.1155/2014/529748

## A New Theory to Forecast the Price of Nonrenewable Energy Resources with Mass and Energy-Capital Conservation Equations

Department of Industrial Engineering, University of Rome “Tor Vergata,” Via del Politecnico 1, 00133 Rome, Italy

Received 31 December 2013; Accepted 1 April 2014; Published 5 June 2014

Academic Editors: K. Abhary and Y. Demirel

Copyright © 2014 Fabio Gori. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The mass and energy-capital conservation equations are employed to study the time evolution of mass and price of nonrenewable energy resources, extracted and sold to the market, in case of no-accumulation and no-depletion, that is, when the resources are extracted and sold to the market at the same mass flow rate. The Hotelling rule for nonrenewable resources, that is, an exponential increase of the price at the rate of the current interest multiplied the time, is shown to be a special case of the general energy-capital conservation equation when the mass flow rate of extracted resources is unity. The mass and energy-capital conservation equations are solved jointly to investigate the time evolution of the extracted resources.

#### 1. Introduction

The price evolution of nonrenewable energy resources is a very important problem in an economy based on nonrenewable energy resources. The economics of exhaustible resources was investigated by Hotelling [1], and reviewed in [2], with the conclusion that the price of a resource increases exponentially with the product of time and interest rate of capital.

After the oil crisis of 1973 the problem of the energy resources became a very important issue and several energy models and forecasts have been developed and employed. During the decade of 1970s, among the others, McCall [3] presented a linear program for the world oil industry which was employed by an oil company. Tewksbury [4] described a practical general approach to forecast foreign oil prices to the year 2000. Pindyck [5] described a new version of an econometric policy model of the natural gas industry. Adelman [6] forecasted a decade of rising real oil prices for the US. DuMoulin and Newland [7] came to the conclusion that a rapid take-off in demand is likely to be followed by a rapid increase in the price of oil.

During the 1980s, among the others, Pearce [8] presented an overall look at the world energy demand forecasts with the conclusion that the judgment, rather than explicit modelling, is used to suggest a world crude oil scenario in which oil prices rise at 10% per year. French [9] discussed the markets of world oil, natural gas, coal, electricity, and the energy demand, with the conclusion that from mid-decade onward energy prices will most likely increase more rapidly than the overall price level. Roberts [10] concluded his analysis that through the mid-1980s, real crude prices should be flat to down and thereafter increase at a rate equal to or less than inflation. Campbell and Hubbard [11] showed that forecasts of prices—including oil, natural gas, equipment and drilling costs, and money—affect the evaluation of all projects. Clark [12] focused on generalized price formation in the U.S. gas industry. David and Carson [13] focused on the banks which use the long term oil supply and demand situation as the underlying support for making oil and gas price projections. Lorentsen and Roland [14] outlined a very open model, which can easily impose different assumptions about oil supply and demand in order to indicate a range of feasible projections for oil prices. Curlee [15] concluded that overall assessment of forecasts and recent oil market trends suggests that prices will remain constant in real terms for the remainder of the 1980s, to increase by 2-3% during the 1990s and beyond. Anon [16] studied the U.S. demand for petroleum products in 1987, based on the assumption that world oil prices will average 13 US $/bbl through the third quarter of 1986 and rise gradually to 18 US $/bbl by the end of 1987. Dielwart and Coles [17] discussed the current market uncertainties and their impact on Canadian oil and natural gas prices. Yohe [18] presented an analytical technique based on an adaptive expectations model of incorporating current information into long-term forecast. Sekera [19] presented three price forecasting analyses with forecasts under each theory and the recommendation that all long term forecasts be based on fundamental/economic analysis, while short term forecasting needs to utilize other methods. Gately [20] reviewed recent and current opinions about prospects for prices in the world market for crude oil with two alternative views. Holden et al. [21] examined the accuracy and properties of forecasts by the OECD for 24 countries and 8 variables. Dougherty [22] presented some rules of thumb that may help to understand the evolution of oil prices, to identify factors to consider when deciding on the credibility of an oil price forecast, rather than to prophesize future oil price. Amano [23] developed an annual, small-scale econometric model of the world oil market to analyze oil market conditions and oil prices for the period 1986–1991.

During the 1990s, among the others, Garces [24] explored the effects of natural gas storage on short-term gas price forecasts, with preliminary results which indicate that by using natural gas storage as a causal variable, the final gas price forecasts can decrease by as much as 5% over the forecast period. Angelier [25] concluded his analytical framework with the forecast that, in the year 2000, oil prices will not be significantly different from those of 1990s. Bopp and Lady [26] concluded their analysis that when actual prices are employed, future prices did correctly anticipate the observed seasonal pattern. Abramson and Finizza [27] developed a system that forecasts crude oil prices via Monte Carlo analyses of the network. Charleson and Weber [28] forecasted energy consumption in Western Australia with Bayesian vector autoregressions to 2010. Fesharaki [29] believed that for the next three to five years oil prices will remain at lower levels than generally predicted, because price increases, even when they occur, will not be sustainable for very long. Huntington [30] reviewed forecasts of oil prices over the 1980s that were made in 1980, identifying the sources of errors due to such factors as exogenous GNP assumptions, resource supply conditions outside the cartel, and demand adjustments to price changes. Moosa and Al-Loughani [31], on the basis of monthly observations on spot and future prices of the West Texas Intermediate (WTI) crude oil, suggested that future prices are neither unbiased nor efficient forecasters of spot prices. Abramson [32] discussed a knowledge-based system that models the crude oil market as a belief network and uses scenario generation and Monte-Carlo analysis to forecast oil prices. Skov [33] reviewed some of the more significant forecasts of energy supply, demand, and oil prices made in the last 25 years and identifies the basic premises used and why they led to incorrect forecasts. Abramson and Finizza [34] described a case study in the use of inherently probabilistic belief network models to produce probabilistic forecasts of average annual oil prices. Santini [35] investigated the statistical model of oil prices examining issues and trends related to both the U.S. and world oil supply. Williams [36], in order to forecast the future, examined historical trends and presents the new idea that there is an inverse relationship between crude oil prices volatility and the strength of the relationship between the owners of crude oil (host governments) and the oil companies that develop, refine, market, and sell the products. Chaudhuri and Daniel [37] demonstrated that the nonstationary behavior of US dollar real exchange rates, over the post-Bretton Woods era, is due to the nonstationary behavior of real oil prices. Pindyck [38] discussed how price forecasts also influence investment decisions and the choice of products to produce. Silvapulle and Moosa [39] examined the relationship between the spot and futures prices of WTI crude oil using a sample of daily data with the result that both spot and futures markets react simultaneously to new information.

During the 2000s, among the others, Alba and Bourdaire [40] consider the actual situation as a “cohabitation” between oil and the other energies with the oil price, extremely volatile, reflecting the trial and error adjustment of the market share left to the other energies. The Centre for Global Energy Studies [41] covers several topics including oilfields set to begin production in 2000, non-OPEC production, demand in western OECD countries and the Asia Pacific region, and oil price forecasts. According to the US energy information [42] global crude oil prices in 2000 will rise to 24 US $/bbl, up by 2 US $/bbl from its earlier estimate, while spot prices for West Texas Intermediate crude will average more than 22 US $/bbl through 2001. Hare [43] reported that according to forecasts by the Energy Information Administration, increases in natural gas prices will continue all the way to the summer and winter of 2000, and that high oil prices will still be seen in 2001. According to the managing director of Shell’s exploration and production [44] crude oil prices will average in the mid-teens of US $/bbl over the next 10 years. The Energy Department of Energy Information Administration [45] forecasts the monthly average price of crude oil imported to the USA to stay above 28 US $/bbl for the rest of 2000. Walde and Dos Santos [46] identified factors that have to be considered in speculating about the future evolution of oil prices. Birky et al. [47] examined the potentiality of unconventional natural gas resources with all other fossil fuels combined with the conclusion of a widespread production and use of these unconventional reserves. The trend of world crude oil price in 1999 has been reviewed and analyzed by Yang [48] and the prospects on world crude oil price in 2000 were predicted. MacAvoy and Moshkin [49] used a model framework consisting of a simultaneous equations system for production from reserves and for demands for production in the residential, commercial, and industrial sectors, with a significant negative trend in the long-term price of natural gas. Morana [50] showed how the GARCH properties of oil price changes can be employed to forecast the oil price distribution over short-term horizons, with a semiparametric methodology based on the bootstrap approach. Tang and Hammoudeh [51] investigated the behavior of the nonlinear model, based on the target zone theory, which has very good forecasting ability when the oil price approaches the upper or the lower limit of the band. Alvarez-Ramirez et al. [52] studied daily records of international crude oil prices using multifractal analysis methods, where rescaled range Hurst analysis provides evidence that the crude oil market is a persistent process with long-run memory effects, demonstrating that the crude oil market is consistent with the random-walk assumption only at time scales on the order of days to weeks. Ye et al. [53] presented a short-term monthly forecasting model of West Texas Intermediate crude oil spot price using OECD petroleum inventory levels, which are a measure of the balance, or imbalance, between petroleum production and demand, and thus provide a good market barometer of crude oil price change. Sadorsky [54] used an ARMAX-ARCH model to estimate the conditional expected returns of petroleum futures prices under time-varying risk with results from a small forecasting experiment which indicates that the out-of-sample forecasts from an ARMAX-ARCH model generally outperform a random walk for all forecast horizons. Fong and See [55] examined the temporal behaviour of volatility of daily returns on crude oil futures using a generalised regime switching model that allows for abrupt changes in mean and variance, GARCH dynamics, basis-driven time-varying transition probabilities, and conditional leptokurtosis. Taal et al. [56] gave a summary of the most common methods used for cost estimation of heat exchange equipment in the process industry and the sources of energy price projections, showing the relevance of the choice of the right method and the most reliable source of energy price forecast used when choosing between alternative retrofit projects or when trying to determine the viability of a retrofit project. Cortazar and Schwartz [57] developed a parsimonious three-factor model of the term structure of oil future prices that can be easily estimated from available future price data. Cabedo and Moya [58] proposed to use value at risk (VaR) for oil price risk quantification, providing an estimation for the maximum oil price change associated with a likelihood level, which fits the continuous oil price movements well and provides an efficient risk quantification. Burg et al. [59] developed an econometric modeling of the world oil market suggesting that world oil prices are likely to fall in the latter part of 2004 and in 2005 as global demand pressures ease. Burg et al. [60] made a market forecast for various energy sources, including oil and gas and coal, with the conclusion that prices are forecasted to decline in 2005 to average 38 US $/bbl WTI. Mirmirani and Li [61] applied VAR and ANN techniques to make ex-post forecast of U.S. oil price movements. Abosedra and Bagesthani [62] evaluated the predictive accuracy of 1-, 3-, 6-, 9-, and 12-month-ahead crude oil future prices for the period from January 1991 until December 2001. Gori and Takanen [63] forecasted the energy demand in a specific country, where the analysis of the electricity demand was focused, for the first time, on the energy consumption and the possible substitution among the different energy resources by using a modified form of the econometric model EDM (Energy Demand Model). Feng et al. [64] predicted the oil price fluctuation by ARFIMA model, which takes the long memory feature into consideration, showing that ARFIMA model is better than ARMA model. Movassagh and Modjtahedi [65] tested the fair-game efficient-markets hypothesis for the natural gas future prices over the period 1990 through 2003. Ye et al. [66] presented a short-term forecasting model of monthly West Texas Intermediate crude oil spot prices using readily available OECD industrial petroleum inventory levels. Yousefi et al. [67] illustrated an application of wavelets as a possible vehicle for investigating the issue of market efficiency in future markets for oil. Sun and Lai [68] presented a forecast model of time series oil price about NYMEX with BP neural networks, analyzing the oil price fluctuation trend. Sadorsky [69] used several different univariate and multivariate statistical models to estimate forecasts of daily volatility in petroleum future price returns. Ye et al. [70] showed the effect that surplus crude oil production capacity has on short-term crude oil prices. Moshiri and Foroutan [71] modelled and forecasted daily crude oil future prices from 1983 to 2003, listed in NYMEX, applying ARIMA and GARCH models, tested for chaos using embedding dimension, BDS(L), Lyapunov exponent, and neural networks tests, and set up a nonlinear and flexible ANN model to forecast the series. The United States Energy Information Administration (EIA) [72] expects worldwide oil demand to increase from 80 million bbl/day in 2003 to 98 million bbl/day in 2015 and to 118 million bbl/day in 2030. The latest forecast reference case calls for crude oil prices to climb from 31 US $/bbl in 2003 to 57 US $/bbl in 2030. Ghouri [73] analyzed qualitatively and quantitatively the relationship between US monthly ending oil stocks position with that of West Texas Intermediate (WTI) oil prices from February 1995 to July 2004, concluding that, if other things are held constant, WTI is inversely related to the petroleum products (PPP), combined petroleum products and crude oil (CPPP), crude oil alone (crude), total oil stocks including petroleum products, crude oil and strategic petroleum reserves SPR (total), total gasoline (TGO), and total distillate (TDO). Ye et al. [74] incorporated low- and high-inventory variables in a single equation model to forecast short-run WTI crude oil prices enhancing the model fit and forecast ability. The approach of using mass and energy-capital conservation equations to investigate the price evolution with time throughout the use of economic parameters was proposed in [75, 76] by the present author. The Hotelling rule was generalized in [75] with the conclusion that the price of the extracted resources increases exponentially with the product of the time and the difference between the inflation rate and the extraction rate of the resources, PIFE, “price increase factor of extracted resources.” A further generalization was done in [76] with the introduction of the price of the selling resources, which depends on PIFE and the new parameter PIFS, “price increase factor of selling resources,” that is, the difference between the prime rate of interest and the extraction rate of the resources. The political events and many complicated factors which happened in the last decades have made oil prices highly nonlinear and even chaotic, with irregularities and random oscillations, linked to the day-to-day events, giving reliability to forecast oil price and consumption in really short terms only, by using adaptive neural fuzzy inference systems (ANFIS), [77]. P. K. Narayan and S. Narayan [78] examined the volatility of the crude oil price using daily data for the period from 1991 to 2006, in various subsamples in order to judge the robustness of the results, implying that the behaviour of the oil prices tend to change over short periods of time. Knetsch [79] developed an oil price forecasting technique, based on the present value model of rational commodity pricing, suggesting the shifting of the forecasting problem to the marginal convenience yield, which can be derived from the cost-of-carry relationship. Deutsche Bank [80] reported that light sweet crude oil prices will average 80 US $/bbl on the New York Mercantile Exchange in 2008 while natural gas prices will average 7.75 US $/MMbtu and perhaps stay put there. It also reported that crude oil in this decade is likely to average nearly 55 US $/bbl. It seems that crude oil markets are following the reciprocal of the declining value of the US dollar and ignoring deteriorating oil fundamentals. MacAskie and Jablonowski [81] developed a decision-analytic model to value commodity price forecasts in the presence of futures markets and applied the method to a data set on crude oil prices. Askari and Krichene [82] observed markets expecting oil prices to remain volatile and jumpy and, with higher probabilities, to rise, rather than fall, above the expected mean. Xu et al. [83] investigated with ILS approach the forecast of the relationship between commodity inventory levels and crude oil spot prices effectively, with the conclusion that their empirical study suggests that both the ILS method and the confidence interval method can produce comparable quality forecasts. Yu et al. [84] proposed empirical mode decomposition (EMD) based on neural network ensemble learning paradigm for world crude oil spot price forecasting, decomposing the original crude oil spot price series into a finite, and often small, number of intrinsic mode functions (IMFs). Fan et al. [85] applied pattern matching technique to multistep prediction of crude oil prices and proposed a new approach: generalized pattern matching based on genetic algorithm (GPMGA), which can be used to forecast future crude oil price based on historical observations. This approach can detect the most similar pattern in contemporary crude oil prices from the historical data. Kang et al. [86] investigated the efficacy of a volatility model for three crude oil markets—Brent, Dubai, and West Texas Intermediate (WTI)—with regard to its ability to forecast and identify volatility stylized facts, in particular volatility persistence or long memory with the conclusion that CGARCH and FIGARCH models are useful for modeling and forecasting persistence in the volatility of crude oil prices. Hamilton [87] examined the factors responsible for changes in crude oil prices, reviewing the statistical behavior of oil prices, in relation to the predictions of theory, and looking in detail at key features of petroleum demand and supply, discussing the role of commodity speculation, OPEC, and resource depletion. Cuaresma et al. [88], using a simple unobserved components model, showed that explicitly modelling asymmetric cycles on crude oil prices improves the forecast ability of univariate time series models of the oil price. Ye et al. [89] predicted crude oil prices from 1992 through early 2004 by using OECD’s relative inventories and OPEC’s excess production capacity improving forecasts for the post-Gulf War I time period over models without the ratchet mechanism. Cheong [90] investigated the time-varying volatility of two major crude oil markets, the West Texas Intermediate (WTI) and the Europe Brent, with a flexible autoregressive conditional heteroskedasticity (ARCH) model to take into account the stylized volatility facts such as clustering volatility, asymmetric news impact, and long memory volatility. Ghaffari and Zare [91] developed a method based on soft computing approaches to predict the daily variation of the crude oil price of the West Texas Intermediate (WTI), with comparison to the actual daily variation of the oil price, and the difference is implemented to activate the learning algorithms. The energy supply curve, that is, price versus consumption, of nonrenewable energy resources was constructed in [92] where new parameters were introduced to forecast the price evolution of nonrenewable resources. The present theory has been reviewed and applied to the period from 1966 until 2006 in [93].

During the 2010s, among the others, de Souza e Silva et al. [94] investigated the usefulness of a nonlinear time series model, known as hidden Markov model (HMM), to predict future crude oil price movements, developing a forecasting methodology that consists of employing wavelet analysis to remove high frequency price movements, assumed as noise, and using the probability distribution of the price return accumulated over the next days to infer future price trends. He et al. [95] presented two forecasting models, one based on a vector error correction mechanism and the other based on a transfer function framework with the range taken as a driver variable, for forecasting the daily highs and lows, showing that both of these models offer significant advantages over the naïve random walk and univariate ARIMA models in terms of out-of-sample forecast accuracy. Miller and Ni [96] examined how future real GDP growth relates to changes in the forecasted long-term average of discounted real oil prices and to changes in unanticipated fluctuations of real oil prices around the forecasts. Forecasts were conducted using a state-space oil market model, in which global real economic activity and real oil prices share a common stochastic trend. Yaziz et al. [97] obtained the daily West Texas Intermediate (WTI) crude oil prices data from Energy Information Administration (EIA) from the 2nd of January, 1986, to the 30th of September, 2009, by using the Box-Jenkins methodology and generalized autoregressive conditional heteroscedasticity (GARCH) approach in analyzing the crude oil prices, which is able to capture the volatility by the nonconstant of conditional variance. Prat and Uctum [98] suggested a mixed expectation model, defined as a linear combination of these traditional processes, interpreted as the aggregation of individual mixing behavior and of heterogeneous groups of agents using these simple processes, consistent with the economically rational expectations theory. It was shown that the target oil price included in the regressive component of this model depends on the long-run marginal cost of crude oil production and on the short term macroeconomic fundamentals whose effects are subject to structural changes. Jammazi and Aloui [99] combined the dynamic properties of the multilayer back propagation neural network and the recent Harr A trous wavelet decomposition, implementing a hybrid model HTW-MPNN to achieve prominent prediction of crude oil price, by using three variants of activation function, namely, sigmoid, bipolar sigmoid, and hyperbolic tangent in order to test the model’s flexibility. Mingming and Jinliang [100] constructed a multiple wavelet recurrent neural network (MWRNN) simulation model, in which trend and random components of crude oil and gold prices were considered to capture multiscale data characteristics, while a real neural network (RNN) was utilized to forecast crude oil prices at different scales, showing that the model has high prediction accuracy. Baumeister and Kilian [101] constructed a monthly real-time dataset consisting of vintages for the period from January, 1991, to December, 2010, that is suitable for generating forecasts of the real price of oil from a variety of models and documenting that revisions of the data typically represent news, and introducing backcasting and nowcasting techniques to fill gaps in the real-time data. It was shown that the real-time forecasts of the real oil price can be more accurate than the nochange forecast at horizons up to 1 year. Azadeh [102] presented a flexible algorithm based on artificial neural network (ANN) and fuzzy regression (FR) to cope with optimum long-term oil price forecasting in noisy, uncertain, and complex environments, incorporating the oil supply, crude oil distillation capacity, oil consumption of non-OECD, U.S. refinery capacity, and surplus capacity as economic indicators. Wang et al. [103] investigated the interacting impact between the crude oil prices and the stock market indices in China and analyzed the corresponding statistical behaviors introducing a jump stochastic time effective neural network model and applying it to forecast the fluctuations of the time series for the crude oil prices and the stock indices, and studying the corresponding statistical properties by comparison. Hu et al. [104] attempted to accurately forecast prices of the crude oil futures by adopting three popular neural networks methods including the multilayer perceptron, the Elman recurrent neural network (ERNN), and the recurrent fuzzy neural network (RFNN), with the conclusion that learning performance can be improved by increasing the training time and that the RFNN has the best predictive power and the MLP has the worst one among the three underlying neural networks. Ruelke et al. [105] derived internal consistency restrictions on short, medium, and long-term oil price forecasts, by analysing whether oil price forecasts extracted from the Survey of Professional Forecasters (SPF), conducted by the European Central Bank (ECB), satisfy these internal consistency restrictions, finding that neither short-term forecast is consistent with medium-term forecasts nor that medium-term forecasts are consistent with long-term forecasts. A new category of cases, that is, the negative inflation rate, has been introduced, within the present theory, in [106], where some preliminary results have been presented. Pierdzioch et al. [107] found that the loss function of a sample of oil price forecasters is asymmetric in the forecast error, indicating that the loss oil price forecasters incurred when their forecasts exceeded the price of oil tended to be larger than the loss they incurred when their forecast fell short of the price of oil. Accounting for the asymmetry of the loss function does not necessarily make forecasts look rational. Shin et al. [108] proposed a study to exploit the method of representing the network between the time-series entities and to then employ SSL to forecast the upward and downward movement of oil prices by using one-month-ahead monthly crude oil price predictions between January, 1992, and June, 2008. Alquist [109] addressed some of the key questions that arise in forecasting the price of crude oil. Xiong et al. [110] proposed a revised hybrid model built upon empirical mode decomposition (EMD) based on the feed-forward neural network (FNN) modeling framework incorporating the slope-based method (SBM), which is capable of capturing the complex dynamic of crude oil prices. The results obtained in this study indicate that the proposed EMD-SBM-FNN model, using the MIMO strategy, is the best in terms of prediction accuracy with accredited computational load. Chang and Lai [111] attempted to develop an integrated model to forecast the cycle of energy prices with respect to the level of economic activity, showing that the oil price cycle and economic activities have bidirectional causality in the short run, and that the upwards (downwards) cycle of oil prices is accompanied by expansion (contraction) of economic activities, and vice versa, with a comovement trend in the long run. Conceptually, the model developed in this work is useful with regard to forecasting the level of economic activities using the oil price cycle, as most economic activities depend on energy. Salehnia et al. [112] used the gamma test for the first time as a mathematically nonparametric nonlinear smooth modeling tool to choose the best input combination before calibrating and testing models and developing several nonlinear models with the aid of the gamma test, including regression models, local linear regression (LLR), dynamic local linear regression (DLLR) and artificial neural networks (ANN) models. Shin et al. [113] exploited the method of representing the network between the time-series entities to employ semisupervised learning (SSL) to forecast the upward and downward movement of oil prices, by using one-month-ahead monthly crude oil price predictions between January, 1992, and June, 2008. Yang et al. [114] forecasted the international crude oil price by using the grey system theory and creating a MATLAB program to achieve it. The present theory has been confirmed in case of negative inflation rate [115], projecting the forecast to the following five months of 2013 [116], with a very recent verification [117].

#### 2. Mass Conservation Equation of Extracted Resources

Assume the mass of nonrenewable energy resources is in the reservoir, , reported in Figure 1. The mass flow rate of extraction from the reservoir is, , with dimensions (mass/annum).

The mass conservation equation of the extracted resources is which, after integration, becomes where is the mass of resources in the reservoir at the initial time .

Assuming the extraction rate, , at the time , given by the mass flow rate of extraction with constant in the time interval is

The evolution of the mass flow rate of extraction, given by (4) and presented in Figure 2, depends on the value of . is constant with time if , increases if , and decreases if .

The mass of nonrenewable resources in the reservoir becomes, dividing by , and its time evolution is presented in Figure 3.

In case of a constant extraction rate, , that is, , (1) gives that is, the mass of resources in the reservoir decreases linearly with time. Equation (6) gives the time when , that is, the time of exhaustion of the nonrenewable resources.

If, for example, the ratio is equal to 50 years, the time to exhaust the resources is years. If , the mass of resources in the reservoir is exhausted faster than for , that is, before 50 years. The time , which gives , is given by (5) as Assuming years, the number of years for the exhaustion of the resources depends on the value of . From (7) it can be found that(i) years for (1/annum);(ii) years for (1/annum).

If , (5) can be written as which, for an infinite time , gives The extraction policy to exhaust the nonrenewable resources in an infinite time is to have which gives the critical extraction rate, , For years the critical extraction rate is (1/annum).

If , the time to exhaust the resources is given by which depends on . For years the time is(i) years for (1/annum);(ii) years for (1/annum).

The owner of the nonrenewable resources does not find convenient to have at an infinite time, , because the resources could remain nonextracted.

If , the nonrenewable resources will never be extracted and a certain amount of resources will remain in the reservoir as nonextracted. For example, if (1/annum) the amount of resources which will remain nonextracted is years.

#### 3. Hotelling Rule

Assume the nonrenewable resource can be extracted without cost. Denote the initial price per unit of resource at the initial time and suppose the interest rate, , is constant in the time interval from to . The owner of the resources unit has two options: to keep the unit of resources in the reservoir for the time interval with the real expectation that the price at will be higher, saying , or to sell the unit of resources and put the money, obtained from the sale, in a bank for the same time interval to get an interest return equal to . Under competitive conditions the owner is indifferent between these two options; then, the price of the unit at , , must be equal to the price at the time , increased of the interest return . In other words or, taking the limit as , it can be obtained After integration, (15) becomes also called Hotelling rule, [2], in honour of Hotelling [1], where, , the interest rate of capital, is assumed constant in the time period .

Equation (16) states that the price of nonrenewable resources can only increase exponentially with time if the interest rate is positive while only decreasing in case of negative inflation rate. The price evolution is independent of any other variable present in the energy game of nonrenewable resources.

#### 4. Energy-Capital Conservation Equation of Extracted Resources

The energy-capital conservation equation can be written in unsteady state, similarly to the energy conservation equation of thermal engineering, for the nonextracted resources of the reservoir , Figure 1, with interest rate (1/annum), as The left term expresses the time variation of the capital , due to the resources of the reservoir , while is the capital flow rate entering into the reservoir , which derives from the mass flow rate of extraction, .

The capital is written as where is the current price of a unity of extracted resources and is the mass flow rate of extraction.

The right term of (17), that is, the capital flow rate entering into the reservoir, is written as where is the interest rate of the nonextracted resources.

Equation (17) then becomes which, for a mass flow rate of extraction being equal to 1, that is, , it gives (15), that is, the Hotelling rule before integration. In conclusion, the Hotelling rule can be obtained from the energy-capital conservation equation with a mass flow rate of extraction being equal to .

Differentiation of the capital conservation equation, (20), and division by the capital gives

Equation (21) can be used, with (3), to study the price evolution of the extracted resources with interest rate . Equation (21) becomes, indicating as for the extracted resources and using (3), which, after integration, gives where, that is, the difference between the interest rate of nonextracted resources, , and the extraction rate, , is a new parameter called “price increase factor of extracted resources,” PIFE. The price evolution with time of the extracted resources, given by (23), is shown in Figure 4.

The price of extracted resources evolves with time according to the value of PIFE, . If the extraction rate is equal to the interest rate of nonextracted resources, , the price increase factor, PIFE is equal to zero, , and the price remains constant with time. If the extraction rate is higher than , the price increase factor, PIFE is negative, , and the price decreases with time. If the extraction rate is lower than , the price increase factor, PIFE is positive, , and the price increases with time. In case of a conservative policy of extraction; that is, , the price increase factor, is greater than zero, , and the price increases with time.

#### 5. Mass and Energy-Capital Conservation Equations of Selling Resources

Assume that nonextracted and extracted energy resources are in the two reservoirs, and , of Figure 5.

The nonextracted resources are in the reservoir with mass and interest rate . The mass flow rate of extracted resources from the reservoir is while the capital flow rate, relative to the mass flow rate , is . The resources of the reservoir are extracted at the price of the extracted resources, , and allocated in the reservoir , where is the mass of extracted resources. The resources of the reservoir are sold at the price . The mass flow rate of the selling resources from the reservoir is while the capital flow rate, relative to the mass flow rate , is .

The mass conservation equation of the extracted resources in the reservoir is while the relative energy-capital conservation equation is which becomes and finally

In case of no-accumulation and no-depletion of the extracted resources, , the right term of (25) is zero; that is, the mass of extracted resources in the reservoir is constant with time. The mass flow rates of extracted and selling resources have the same extraction rate from the reservoirs and ; that is, which gives, after integration, and the time evolution is similar to that of Figure 2.

Equation (28) becomes or

The substitution of (23) into (32) and the solution of the relative differential equation give the time evolution of the price of sold resources where is called the “Price Increase Factor of Selling resources,” PIFS, and the “critical initial price of selling resources,” CIPS.

The price of selling resources, (33), has an extreme (maximum or minimum) for the time given by The extreme is a maximum if otherwise it is a minimum.

The time of the extreme is zero if is equal to which is called “critical initial price extreme of selling resources,” CIPES.

The time of the maximum is greater than zero if or For the time of the maximum is for , or for .

For , that is, , is given by where CIPS, that is, , is not defined.

The price of selling resources has an extreme for the time given by

The extreme is a maximum if and a minimum if . The time of the maximum is zero if the initial price is

The possible cases can be classified into six categories:(i)Category 0—;(ii)Category 1: , that is, ;(iii)Category 2: , that is, ;(iv)Category 3: , that is, ;(v)Category 4: , that is, or ;(vi)Category 5: .On their side, each category presents the following possible cases.

*Category 0 (**).* If the interest rate of nonextracted resources is , (i.e., ), the price evolution of selling resources is
where is the PIF of sold resources*. *

Figure 6 presents the price evolution of selling resources according to the relation between extraction rate, , and interest rate of extracted resources .

*Case 0-A (**, (i.e., **)).* The price increases with time at different rates, according to the values of . The rate of increase is higher for , , intermediate for , , and lower for , .

*Case 0-B (**, (i.e.,* *, **)).* The price is constant with time at .

*Case 0-C (**, (i.e.,* *, **)).* The price is decreasing with time towards .

In conclusion, for the criterium to evaluate the price evolution of selling resources, , is . If the price increases, if the price remains constant, and if the price decreases towards .

*Category 1 (**, that is, **).* The price evolutions of are reported in Figure 7 according to the relation between and .

*Case 1-A (**, that is, **).* The price of selling resource increases with time if the initial price is , that is, equal to or greater than .

*Case 1-B (** and ** or ** and **).* The price of selling resource increases temporarily with time up to , given by (19), and then decreases for if the initial price has a value comprised between CIPS and CIPES, that is, , or, for , if the initial price is greater than CIPES, that is, .

*Case 1-C (** and **, or ** and **, or **).* The price of selling resource decreases with time for if the initial price is smaller than CIPS and CIPES, that is, , or, for , if the initial price is smaller than CIPES, that is, , or, for for every initial price because .

*Category 2 (**, that is, **).* The time evolution is presented in Figure 8 according to and or the relation between and .

*Case 2-A (**, that is, **).* The price of selling resource increases with time if the initial price is , that is, equal to or greater than , .

*Case 2-B (**, that is, **).* The price of selling resource decreases temporarily with time up to , given by (19), and then increases if the initial price, , has a value comprised between CIPS and CIPES, that is, .

*Case 2-C (**, that is, **).* In all other cases the price of selling resource decreases with time including *,* or , or .

*Category 3 (**, that is, **).* The price evolution of selling resources, , is dependent on the relation between and , as reported in Figure 9.

*Case 3-A (**, that is, **).* The price of selling resource increases with time if the initial price is , that is, greater than , .

*Case 3-B (**, that is, **).* The price of selling resource remains constant with time if the initial price is , that is, equal to , *.*

*Case 3-C (**, that is, **, or **, that is, **).* The price of selling resource, for , decreases with time if the initial price is , that is, smaller than CIPS and CIPES, . In the other cases the price of the selling resource decreases with time for every value of .

*Category 4 (**, that is, **).* The time evolutions are presented in Figure 10.

*Case 4-A (**, that is, **).* The price cannot increase with time.

*Case 4-B (**, that is, **).* The price of selling resource increases temporarily with time up to , given by (2), and then decreases if the initial price is , that is, greater than .

*Case 4-C (**, that is, **, or **, that is, **, or **, that is, **).* The price of selling resource decreases with time if the initial price is , that is, equal to or smaller than . In all the other cases the price of selling resource decreases with time for every value of .

*Category 5 (**).* A new category of cases has been introduced for negative inflation rate, , which was registered from March to October, 2009, after the economic crisis of 2008. This new category is interesting because the Hotelling rule forecasts a decrease of the oil price if the interest rate is assumed to be equal to the inflation rate, that is, negative, while the conclusions of the present approach are different.

The time evolutions are presented in Figure 11.

*Case 5-A (**, that is, **).* The price of selling resource increases with time for every value of the initial price because , since CIPS and CIPES are negative, or .

*Case 5-B (**, that is, **).* The price of selling resource increases temporarily with time up to and then decreases if the initial price, , is , that is, smaller than .

*Case 5-C (**, that is, **).* The price of selling resource decreases with time if the initial price, , is , that is, greater than *.*

#### 6. Price Difference between Selling and Extracted Resources

The difference between selling and extracted resources is defined as where is the critical initial price difference, CIPD, .

The price difference has an extreme (maximum or minimum) for

The extreme is a maximum if otherwise it is a minimum.

The time of the extreme is zero if is equal to where is the critical initial extreme of the initial price difference, CIEIPD.

The time of the extreme is greater than zero if For the time of the extreme is for , and for .

The time is when .

The possible cases can be classified in five categories: Category 0D: ; Category 1D: , that is, ; Category 2D: , that is, ; Category 3D: , that is, ; Category 4D: , that is, or .

On their side, each category presents the following possible cases.

*Category 0D (**).* The difference between selling and extracted resources, , is function of the initial difference and the extraction rate, , as compared to . It can be remarked that for the critical initial price and differences of sold resources are , , and .

The price difference is then given by

The time of the extreme is given by which is obtained for and the extreme is a maximum for

The price difference is at the time

Four initial conditions are investigated:(i)if the initial price difference is or ;(ii)if the initial price difference is ;(iii)if the initial price difference is ;(iv)if the initial price difference is or .

The time evolutions are presented in Figure 12.

*Case 0D-A (**, (i.e., **, **, **)*

*Case 0D-A1*. If the initial , increases exponentially towards with a higher rate for , , an intermediate one for , , and a smaller one for , *.*

*Case 0D-A2*. If the initial is the evolution is similar to the previous case but now increases from a negative value and all the curves pass through the same point at . If the initial , increases exponentially with the minimum at .

*Case 0D-A3*. If the initial is and , decreases initially to a minimum and then increases exponentially as in the previous cases.

*Case 0D-B (**, (i.e., **, ** is not defined and **))*

*Case 0D-B1*. If the initial , tends towards as asymptote.

*Case 0D-C (**, (i.e., **, **))*

*Case 0D-C1*. If the initial value , the difference decreases towards , with the maximum at for .

*Case 0D-C2*. If the initial value , the difference increases up to the time when reaches a maximum value, , and then decreases towards .

*Case 0D-C3*. If the initial , increases towards 0*.*

*Case 0D-C4*. If the initial value , decreases towards for , tends to 0 for , and remains constant for *.*

In conclusion, for , the criteria to evaluate the evolution of are the extraction rate and the initial value . If , increases towards if the initial is . If , tends towards the asymptotic value , for . If , tends towards if .

*Category 1D (**, that is, **).* The price evolution of the difference between selling and extracted resources, , is investigated according to the initial value and the extraction rate , as compared to .

*Case 1D-A (**, (i.e., **, **))*. Three relations between and can be investigated:(i), that is, ;(ii), that is, ;(iii), that is, .

If , that is, , the following four initial conditions are investigated:(i)the initial price difference is (or );(ii)the initial price difference is ;(iii)the initial price difference is ;(iv)the initial price difference is , or ;

and the price evolutions are presented in Figure 13.

*Case 1D-A1*. If the difference increases exponentially towards with a higher rate for , an intermediate one for , and a smaller one for .

*Case 1D-A2*. If the difference increases exponentially towards as in the previous case and the curves pass through the same point at *. *If the difference increases exponentially with the minimum at .

*Case 1D-A3*. If , decreases initially down to a minimum and then increases exponentially as in the previous cases.

If , that is, , the following cases are discussed.

*Case 1D-A1*. If the difference increases exponentially towards with a higher rate for , an intermediate one for , and a smaller one for .

If , that is, , the following cases are discussed.

*Case 1D-A1*. If the difference increases exponentially towards with a higher rate for , an intermediate one for , and a smaller one for .

*Case 1D-B (**). * remains constant at with time.

*Case 1D-C1 (**, (i.e., **, **, **, and **)).* If is ,* Δp* decreases towards *. *If is , increases up to a maximum and then decreases towards . If , decreases towards .

*Case 1D-C2 (**, (i.e., **, **, **)).* The difference decreases towards for every initial .

In conclusion, for , (i.e., ) the criteria to evaluate the price difference evolution are the extraction rate , the initial condition, and the relation between and . For and , increases towards if . For and , increases towards if . For and , increases towards if , while increases up to a maximum and then decreases if *.* For , decreases towards *.*

*Category 2D (**, that is, **).* The price evolution of the difference is investigated for according to the initial value and the relation between and .

The price evolutions are presented in Figure 14.

*Case 2D-A (**, (i.e., **)).* If , that is, , the following situations are discussed.

*Case 2D-A1*. If the difference increases towards .

If , that is, , the following situations are discussed.

*Case 2D-A1*. If the difference increases towards . If is the difference increases towards after a minimum at .

If , that is, , the following situations are discussed.

*Case 2D-A1*. If the difference increases towards .

*Case 2D-B (**, (i.e., **, **, ** for **, ** for **, and ** for **))*

*Case 2D-B1*. If , the difference remains constant at . If the difference tends to .

*Case 2D-B2*. If , the difference decreases towards . For , remains constant at .

*Case 2D-C (**, (i.e., **, ** for **, ** for **, and ** for **))*

*Case 2D-C1*. If is the difference decreases towards after a maximum at .

*Case 2D-C2*. If , the difference decreases towards *.* If the difference decreases towards . If the difference decreases towards *.*

For , tends towards 0, if , increases up to a maximum and then tends towards +0, if , and tends towards , if . For , decreases to if , remains constant at if , and increases towards −0 if . For , decreases towards 0 for every .

*Case 2D-D (**, (i.e., **, **)).* For every the difference decreases down to a minimum and then increases towards −0.

In conclusion, for , (i.e., ) the criteria to evaluate the evolution of are the extraction rate and the initial condition. For , increases towards if , or or , decreases towards if or . For , tends towards (). If , tends towards for every initial except for and when tends towards 0 after a maximum.

*Category 3D (**, that is, **).* The price evolution of the difference is investigated for , (i.e., ) according to the initial value and the relation between and .

The price evolutions are reported in Figure 15.

*Case 3D-1 (**, that is, **)*

*Case 3D-1A*. If the difference increases exponentially with time towards .

*Case 3D-1B*. If the difference remains constant at the initial value .

*Case 3D-1C*. If the difference decreases towards .

*Case 3D-2 (**, (i.e., **)).* The difference decreases towards for every initial *.*

In conclusion, for , the criteria to evaluate the evolution of are the relation between and and the initial condition . For , increases towards if , remains constant if , and decreases towards if . For , tends towards the asymptotic value .

*Category 4D (**, that is, **).* The price difference is given by

The critical initial price difference and the critical initial extreme price difference are not defined and the only critical price defined is the critical initial price extreme of the sold resources , already defined in the second part of (38) or (45).

The price difference has an extreme for the time given by The extreme is a maximum if and a minimum if .

The time of maximum is zero if and is greater than zero if

The difference is zero for the time equal to

The price evolution of the difference is investigated according to the initial value and the extraction rate .

The price evolutions are reported in Figure 16.

*Case 4D-1 (**, (i.e., **))*

*Case 4D-1A*. The price cannot increase with time.

*Case 4D-1B*. If the difference initially increases up to the maximum time and then decreases towards .

*Case 4D-1C*. If the difference decreases towards with the maximum at . If , the difference decreases towards .

*Case 4D-2 (**, (i.e., **, **))*

*Case 4D-2A*. The price cannot increase with time.

*Case 4D-2B*. If , the difference decreases with time, reaches a minimum at , and then increases towards −0. If the difference increases towards 0.

*Case 4D-3 (**, (i.e., **, **)).* The difference decreases linearly towards for every initial .

*Discussion*. The discussion on the price evolutions of sold resources, , and price difference, , in case of no-accumulation and no-depletion of the extracted resources, is carried out summarizing the results in the following Tables 1–8.

Table 1 presents the price trends of , , and , for nonrenewable resources without interest rate, = 0, that is, , , , and with the following conclusions:(i)the price of extracted resources increases towards if , (i.e., ), remains constant at if , (i.e., ), and decreases towards if ;(ii)the price of sold resources increases towards if , (i.e., ), remains constant at if , (i.e., ), and decreases towards if , (i.e., );(iii)if the extraction rate , (i.e., ), the trends of the price difference depend on the initial value, . If , increases towards , if , increases towards after a minimum, and if , decreases towards ;(iv)if the extraction rate is , (i.e., ), the price difference increases towards if while tending towards if . It can be noticed that cannot be smaller than because cannot be smaller than 0;(v)if the extraction rate , (i.e., ), the price difference increases towards if while tends towards 0 if ;(vi)if the extraction rate , (i.e., ), the price difference tends towards .;(vii)if the extraction rate , (i.e., ), the price difference tends towards for any but with a maximum at if .

Table 2 presents the price trends of , , and , for , that is, *, *with the following conclusions:(i)the price increases towards , because ;(ii)the price increases towards if , (i.e., ), and ;(iii)the price initially increases, reaches a maximum at , and then decreases towards in the following cases: and , and , and and ;(iv)the price, , decreases towards , in all the other cases;(v)the price difference has different trends for according to the relation between and . If , the price difference increases towards if , increases towards after a minimum if , and decreases towards if . If , the price difference increases towards if , remains constant at if , and decreases to if . If , the price difference increases towards if , decreases towards , after a maximum, if , and decreases to if ;(vi)the price difference for initially increases up to a maximum at and then decreases towards if ;(vii)the price difference decreases towards in all other cases.

Table 3 presents the price trends of , , and , for , that is, , with the following conclusions:(i)the price is constant at , because ;(ii)the price increases towards , for , and ;(iii)the price is constant at for and or for and every ;(iv)the price decreases towards for and or decreases linearly towards for and every ;(v)the price difference increases towards for , and ;(vi)the price difference is constant at , for , and or for and every ;(vii)the price difference decreases towards , for , and or linearly to for and every .

Table 4 presents the price trends of , , and for , that is, , with the following conclusions:(i)the price decreases towards , because ;(ii)the price increases towards for (i.e., ) and or after a minimum for ;(iii)the price tends to for ;(iv)the price decreases towards for and ;(v)the price decreases towards 0 without or with a minimum, in the other two cases;(vi)the price difference increases towards for and , if , if , if and . If a minimum is present if ;(vii)the price difference tends to if for every ;(viii)the price difference decreases towards for and , or for while passes through a maximum if and for ;(ix)the price difference remains constant at for , and ;(x)the price difference tends to 0 in the other three cases.

Table 5 resumes the trends of for , that is, , .

The price increases towards if , remains constant at if , and decreases to if .

Table 6 resumes the trends of for , that is, . The trends can be divided according to the initial conditions which are grouped into four categories.

The greater value between and is indicated as while the smaller one as . Consider the following:(i)for , increases to if , tends to if , and tends to if ;(ii)for , increases to if , remains constant at if , tends to 0 if , and , tends to if ;(iii)for , increases to after a minimum if , and tends to after a maximum if ;(iv)for , decreases to if , remains constant at if , and tends to 0 if .

Table 7 resumes the trends of for different from 0. The trends can be divided into four classes, A-B-C-D, according to the relation between and , and . In each class the trends are depending on the value of or , and the initial value . The greater between and is indicated as and the smaller one as .

In class A, , the trends of are to increase towards for if , for if , and for if and if after a minimum. A temporary increase, followed by a decrease towards , is present if