黎斌林,申 維

(中國(guó)地質(zhì)大學(xué)(北京)地球科學(xué)與資源學(xué)院,北京 100083)

石油峰值(Peak Oil)即石油產(chǎn)量的“鐘形”(Bell-Shaped)規(guī)律，不僅反應(yīng)了資源的物理特征(Cleveland and Kaufmann,1991[1];Moroney and Berg,1999[2]),也反應(yīng)了影響油氣開(kāi)采的地質(zhì)因素。從石油峰值的內(nèi)涵上考量，石油峰值并不意味著石油枯竭，而是產(chǎn)量達(dá)到峰值后，逐步開(kāi)始下降。(Robert Rapier，2011[3])

針對(duì)石油峰值(Peak Oil)的爭(zhēng)論在全球的學(xué)術(shù)圈中從未休止。一些石油地質(zhì)學(xué)家都是石油峰值理論的擁躉者。如愛(ài)爾蘭地質(zhì)學(xué)家Campbell,法國(guó)著名的石油地質(zhì)學(xué)家Jean Laherrere等。許多都為經(jīng)濟(jì)學(xué)領(lǐng)域的學(xué)者更傾向于石油峰值論的反對(duì)者,如美國(guó)經(jīng)濟(jì)學(xué)家 Michael Lynch(2003)[4]，反對(duì)理由集中在：價(jià)格會(huì)調(diào)節(jié)人類(lèi)的能源消費(fèi)結(jié)構(gòu)。石油價(jià)格的上漲或者下降會(huì)調(diào)節(jié)人們使用其他的替代能源或者開(kāi)采更高難度的石油。技術(shù)進(jìn)步因素也會(huì)促使非常規(guī)石油向常規(guī)石油轉(zhuǎn)變,這些均決定了石油峰值理論只是一個(gè)過(guò)于悲觀的論題。也有一些學(xué)者也持有其他的觀點(diǎn)，如Murray(2012)[5]認(rèn)為石油產(chǎn)量的峰值點(diǎn)已經(jīng)過(guò)去了。Kerr(2012)[6]認(rèn)為美國(guó)及全球石油產(chǎn)量持續(xù)的增長(zhǎng)是由于對(duì)于非傳統(tǒng)石油的開(kāi)發(fā)。Mill(2012)[7]認(rèn)為不存在石油短缺的問(wèn)題，但卻認(rèn)為石油需求的峰值會(huì)晚于供給的峰值。

對(duì)全球的石油峰值預(yù)測(cè)所使用的模型主要有三種：第一種是Hubbert及修正后的Hubbert模型,如A1-Fattah and Startzman (2000)[8])；第二種方式是圖解模型。如Aleklett,2010[9]使用WOCAP模型預(yù)測(cè)石油峰值，Bakhtiari(2004)[10]所使用的圖解模型。第三種方式是Long et al.(2004)[11]提出的，石油的產(chǎn)量下降是從儲(chǔ)產(chǎn)比10開(kāi)始的。但更多的學(xué)者都認(rèn)為Hubbert模型仍然是預(yù)測(cè)石油峰值最核心的理論[8]。

本文中所指的石油為常石油，包括了原油、凝析油、液態(tài)天然氣，也包括了部分經(jīng)濟(jì)可采的委內(nèi)瑞拉重油資源，美國(guó)頁(yè)巖油，加拿大油砂等。

1 石油峰值預(yù)測(cè)

1.1 單循環(huán)Hubbert預(yù)測(cè)

Hubbert在1956年的論文Nuclear Energy and The Fossil Fuels首次提出了Hubbert模型，模型中的增長(zhǎng)曲線(xiàn)的核心基礎(chǔ)如下。

1) 對(duì)于任何的稀缺性資源的生產(chǎn)曲線(xiàn)，從一開(kāi)始曲線(xiàn)上的兩點(diǎn)已知，也就是在t=0和t=∞，初始生產(chǎn)率以及資源耗盡，產(chǎn)量均為0。意味著對(duì)任何有限性的資源而言，在經(jīng)歷了一個(gè)或者幾個(gè)極大值(maxima)產(chǎn)量后，最終產(chǎn)量將為0。

2) 第二個(gè)考量是來(lái)自積分，也就是，如果存在一個(gè)單值函數(shù)，y=f(x)。那么就存在

陳元千(1996)對(duì)Hubbert模型進(jìn)行了推導(dǎo)的Hubbert模型可以表述為[12]式(1)～(4)。

(1)

(2)

(3)

(4)

式中：NP為石油的累計(jì)產(chǎn)量，NR為可采儲(chǔ)量，t為開(kāi)發(fā)時(shí)間，Q為石油產(chǎn)量，NPM是累計(jì)可采儲(chǔ)量，tMAX是石油最大產(chǎn)量所發(fā)生的時(shí)間。

最終可采資源量是累計(jì)產(chǎn)量，剩余儲(chǔ)量，儲(chǔ)量增長(zhǎng)以及待發(fā)現(xiàn)的資源量之和[13]。本文選取美國(guó)地質(zhì)調(diào)查局(USGS)對(duì)待發(fā)現(xiàn)資源量，儲(chǔ)量增長(zhǎng)，BP的儲(chǔ)量與產(chǎn)量數(shù)據(jù),估算全球的最終可采資源量(URR)為4700Gb。

產(chǎn)量數(shù)據(jù)選取1860～2012年全球石油產(chǎn)量的數(shù)據(jù)，1860～1960年的數(shù)據(jù)來(lái)自于“Oil economist handbook”，1960～2012年的數(shù)據(jù)來(lái)自于BP統(tǒng)計(jì)，共153個(gè)產(chǎn)量數(shù)據(jù)。核心的b值的確定依據(jù)線(xiàn)性Hubbert方法確定，如圖1所示，橫坐標(biāo)為累計(jì)產(chǎn)量，縱坐標(biāo)為產(chǎn)量與累計(jì)產(chǎn)量之比，選擇擬合程度最高的數(shù)據(jù)，確定b值的確定結(jié)果是0.046。

依據(jù)上述模型和前述最終可采資源量以及累計(jì)產(chǎn)量，選取兩個(gè)不同的最終可采資源量(4700Gb和3000Gb)，累計(jì)產(chǎn)量取1300Gb，分別確定相應(yīng)的a值，依據(jù)單循環(huán)Hubbert模型，如圖2所示，未來(lái)的全球常規(guī)石油峰值將會(huì)在2034年出現(xiàn)，峰值產(chǎn)量為53.46Gb/a，與當(dāng)前的實(shí)際產(chǎn)量相比(2012年全球石油產(chǎn)量)，增長(zhǎng)41.2%。如最終可采資源量為3000Gb,依據(jù)Hubbert模型，全球石油峰值會(huì)提前至2019年，相應(yīng)的最大產(chǎn)量也減少至34.12Gb。

Laherrere及一些石油峰值理論的專(zhuān)家在對(duì)全球石油峰值的預(yù)測(cè)中，廣泛的應(yīng)用分層曲線(xiàn)(creaming curves)圖來(lái)反映石油峰值。因此，論文中對(duì)于全球石油峰值的預(yù)測(cè)都采用分層曲線(xiàn)圖來(lái)圖示全球石油產(chǎn)量(圖2)。

圖1 線(xiàn)性Hubbert

圖2 單循環(huán)Hubbert模型對(duì)全球石油產(chǎn)量的預(yù)測(cè)

1.2 多循環(huán)Hubbert模型

單循環(huán)Hubbert模型之所以成功的預(yù)測(cè)了美國(guó)的石油峰值，是因?yàn)槊绹?guó)的石油市場(chǎng)符合一些必要條件：如較為穩(wěn)定的能源市場(chǎng)；消費(fèi)量呈指數(shù)增長(zhǎng)等。但現(xiàn)實(shí)是，由于油價(jià)，資源稀缺性加劇等因素的制約，峰值之后的產(chǎn)量遞減率會(huì)高于峰值前產(chǎn)量增長(zhǎng)率，對(duì)稱(chēng)鐘形分布的單循環(huán)Hubbert曲線(xiàn)與實(shí)際不符，石油產(chǎn)量曲線(xiàn)并不會(huì)呈現(xiàn)對(duì)稱(chēng)鐘形分布。因此Hubbert模型朝多循環(huán)模型演進(jìn)，A1-Fattah與Startzman提出了多循環(huán)Hubbert模型(A1-Fattah and Startzman,2000[8])。

(5)

式中k是生產(chǎn)循環(huán)數(shù)，qmax和tmax分別為每一個(gè)循環(huán)中的生產(chǎn)率和相應(yīng)的時(shí)間，a是常數(shù)。

每一個(gè)生產(chǎn)循環(huán)都有其自身的最終可采儲(chǔ)量Q∞,公式為

(6)

多個(gè)循環(huán)的最終可采儲(chǔ)量可以表示為

(7)

那么其邏輯方程，也就是代表著多循環(huán)的累計(jì)產(chǎn)量為

(8)

未來(lái)的可采石油量就為

NFR=Npa-Q

(9)

在對(duì)多循環(huán)Hubbert模型的求解中，首先根據(jù)對(duì)每個(gè)國(guó)家的歷史石油產(chǎn)量數(shù)據(jù)，確定循環(huán)數(shù)，然后使用模型使用非線(xiàn)性方法來(lái)計(jì)算產(chǎn)量。依據(jù)全球石油產(chǎn)量歷史(1860～2012年)，確定兩個(gè)循環(huán)數(shù)，第一個(gè)峰值年為1979年，相應(yīng)的最終可采資源量依據(jù)前述公式計(jì)算近似為1000Gb。那第二個(gè)循環(huán)的最終可采儲(chǔ)量相應(yīng)的為3700Gb。依據(jù)多循環(huán)HUBBERT模型所計(jì)算出的結(jié)果發(fā)現(xiàn)，全球常規(guī)石油的峰值會(huì)出現(xiàn)在2037年，最大石油產(chǎn)量值為42.87Gb(圖2),與單循環(huán)Hubbert 模型所計(jì)算出的53.46Gb相比，峰值產(chǎn)量減少了近20%。相應(yīng)的峰值年推后了3年(圖3)。

圖3 多循環(huán)Hubbert模型對(duì)全球石油產(chǎn)量預(yù)測(cè)

2 石油峰值預(yù)測(cè)對(duì)比研究

最終可采資源量(URR)值的不同是許多學(xué)者采用了同樣的Hubbert模型，但最終預(yù)測(cè)值差異較大的最主要的原因。因此URR值的確定對(duì)峰值預(yù)測(cè)至關(guān)重要。有一些學(xué)者運(yùn)用線(xiàn)性趨勢(shì)推斷最終可采資源量，這是不科學(xué)的。本文根據(jù)美國(guó)地質(zhì)調(diào)查局(USGS)2012年對(duì)于全球儲(chǔ)量增長(zhǎng)的評(píng)估，待發(fā)現(xiàn)資源量的評(píng)估，BP的剩余儲(chǔ)量及累計(jì)產(chǎn)量的數(shù)據(jù)，依據(jù)AlHusseini(2006)[13]對(duì)于最終可采資源量的定義，估算全球的最終可采資源量(URR)在4700Gb左右。Chih-Ming Tien(2011)[14]使用數(shù)學(xué)模型和以往的文獻(xiàn)綜述表明，全球的最終可采資源量在2.5Gb～5Gb(109桶)之間。全球的最終可采資源量是一個(gè)動(dòng)態(tài)值，某些傳統(tǒng)的非常規(guī)石油資源隨著技術(shù)進(jìn)步，價(jià)格的上升逐步變成經(jīng)濟(jì)可采資源，最終可采資源量會(huì)出現(xiàn)增加。

文章采用兩個(gè)不同的URR值4700Gb與3000Gb，使用單循環(huán)Hubbert模型，所預(yù)測(cè)峰值產(chǎn)量相差23.46Gb,相應(yīng)的峰值年也推遲了20年(表1)。進(jìn)一步證明了盡管學(xué)者采用同樣的Hubbert模型，最終可采資源量的取值才是峰值預(yù)測(cè)出現(xiàn)較大差異的最主要的原因。

表1 石油峰值預(yù)測(cè)值/Gb(109桶)

相同的最終可采資源量4700Gb,運(yùn)用多循環(huán)Hubbert模型所計(jì)算出的峰值產(chǎn)量與單循環(huán)Hubbert模型相比減少了近20%，而且峰值年也推后了3年 (表1) 。同比對(duì)比研究多循環(huán)Hubbert模型與單循環(huán)Hubbert模型相比，具有以下的優(yōu)點(diǎn)(圖2)。

1) 多循環(huán)Hubbert模型所預(yù)測(cè)的全球石油產(chǎn)量峰值之后的產(chǎn)量遞減率高于峰值前產(chǎn)量增長(zhǎng)率，克服了單循環(huán)Hubbert曲線(xiàn)呈對(duì)稱(chēng)鐘形分布的缺陷。

2) 多循環(huán)Hubbert模型更加符合歷史上的由于技術(shù)、經(jīng)濟(jì)和政治等影響的產(chǎn)量數(shù)據(jù)的波動(dòng)，擬合程度遠(yuǎn)高于單循環(huán)Hubbert模型。而且單循環(huán)Hubbert產(chǎn)量曲線(xiàn)由于本身數(shù)學(xué)模型的原因，預(yù)測(cè)模型會(huì)夸大峰值產(chǎn)量。

3) 多循環(huán)Hubbert模型可以很方便很容易地使用歷史產(chǎn)量數(shù)據(jù)，而且當(dāng)新產(chǎn)量數(shù)據(jù)出現(xiàn)時(shí)，可以很輕易的更新估算。并且方程有三個(gè)未知參數(shù)，相較傳統(tǒng)的方程的四個(gè)參數(shù)，更方便計(jì)算。

在對(duì)石油峰值預(yù)測(cè)的學(xué)者中，有一些學(xué)者較為悲觀，認(rèn)為石油峰值已經(jīng)過(guò)去，或者即將來(lái)臨。如Bakhtiaria(2004)使用WOCAP模型預(yù)測(cè)了全球的石油峰值在2006～2007年。Fisher(2008)認(rèn)為全球的峰值是在20世紀(jì)70年代。Duncan(2000)預(yù)測(cè)石油的峰值在2006年，相應(yīng)的峰值年產(chǎn)量為32Gb左右。Deffeyes(2005)應(yīng)用Hubbert線(xiàn)性模型預(yù)測(cè)石油峰值在2005年[15]。Laherrere在2003年的預(yù)測(cè)中，最終可采資源量以2150Gb為計(jì)算值，最終預(yù)測(cè)的峰值年在2010年左右，峰值產(chǎn)量為27.4Gb，這與2010年的實(shí)際產(chǎn)量相差了2.8億桶。2013年Laherrere的預(yù)測(cè)延續(xù)了其之前的預(yù)測(cè)，認(rèn)為我們當(dāng)前正處于石油峰值期[16]。

全球更為主流的石油峰值預(yù)測(cè)認(rèn)為石油峰值會(huì)來(lái)臨，但由于不斷提升的地質(zhì)知識(shí)、技術(shù)以及石油價(jià)格的上升，全球過(guò)去所謂的非常規(guī)石油逐步成為經(jīng)濟(jì)可采石油，導(dǎo)致全球的石油峰值大大延遲。PFC Energy(2004)認(rèn)為石油峰值會(huì)發(fā)生在2018年之后。IEA(2006,2008,2012)認(rèn)為全球的石油峰值將會(huì)推遲到2030年之后，殼牌(Shell,2012)和?？松梨?Exxon Mobile，2012)的展望報(bào)告也都認(rèn)為石油峰值會(huì)推遲到2030年之后。

文章所采用的多循環(huán)Hubbert模型對(duì)于全球石油峰值的預(yù)測(cè)也符合主流機(jī)構(gòu)的看法，全球的石油峰值延至2030年之后，即石油產(chǎn)量在2037年左右逐步達(dá)到峰值，年產(chǎn)量相較當(dāng)前石油產(chǎn)量，增長(zhǎng)26%左右，年均增長(zhǎng)率在1.5%左右，這也非常符合當(dāng)前1.6%的年增長(zhǎng)率(2006～2011年)。單循環(huán)Hubbert生產(chǎn)曲線(xiàn)與之相比，峰值年的預(yù)測(cè)較合理，但明顯夸大未來(lái)峰值產(chǎn)量。

3 結(jié)論與討論

1) 全球的石油峰值會(huì)到來(lái)，只是因?yàn)榻陙?lái)技術(shù)進(jìn)步，油價(jià)上升，地質(zhì)領(lǐng)域知識(shí)的進(jìn)步，傳統(tǒng)非常規(guī)石油逐步經(jīng)濟(jì)可采使得全球的石油峰值推遲發(fā)生。運(yùn)用多循環(huán)Hubbert模型對(duì)于全球石油峰值的預(yù)測(cè)也符合主流機(jī)構(gòu)的看法，石油產(chǎn)量在2037年左右逐步達(dá)到峰值，峰值產(chǎn)量相較當(dāng)前石油產(chǎn)量，增長(zhǎng)26%左右，為42Gb左右。

2) 最終可采資源量(URR)依據(jù)美國(guó)地質(zhì)調(diào)查局的評(píng)估等數(shù)據(jù)確定為4700Gb，使用單循環(huán)Hubbert模型并與URR3000Gb對(duì)比研究，進(jìn)一步證明了盡管學(xué)者采用同樣的Hubbert模型，最終可采資源量的取值才是峰值預(yù)測(cè)出現(xiàn)較大差異的最主要的原因。

表2 不同機(jī)構(gòu)和學(xué)者的石油峰值預(yù)測(cè)

3) 與單循環(huán)Hubbert模型相比，多循環(huán)Hubbert產(chǎn)量曲線(xiàn)峰值后的產(chǎn)量遞減率高于峰值前產(chǎn)量增長(zhǎng)率，更符合稀缺性資源的產(chǎn)量特征和歷史產(chǎn)量數(shù)據(jù)的波動(dòng)。最終預(yù)測(cè)值對(duì)比研究也證實(shí)一定程度上克服了其夸大未來(lái)峰值產(chǎn)量的缺陷。

文章由于篇幅所限，對(duì)最終可采資源量的確定過(guò)程未能詳述。

[1] Cleveland C J,R K Kaufmann.Forecasting ultimate oil recovery and its rate of production:Incorporating economic forces into the models of M.King Hubbert[J].Energy,1999,12(2):17-46.

[2] Moroney J R,Berg M D.An integrated model of oil production[J].Energy,1999,20(1)：105-124.

[3] Robert Rapier.Five Misconceptions About Peak Oil.[EB/OL].http:∥www.theoildrum.com/node/8585,2010.

[4] Lynch,M.C.The New Pessimism about Petroleum Resources:Debunking the ‘Hubbert Model’(and Hubbert modelers)[J].Minerals and Energy,2003,18(1):21-32.

[5] Murray and King.Oil’s Tipping Point has passed[J].Nature,2012(5):433-435.

[6] Kerr.Technology is turning US oil around but not the World’s [J].Science,2012,335(7):522-523.

[7] Mill.The world has plenty of oil.[EB/OL].http:∥www.europeanenergyreview.eu/site/pagina.php?id_mailing=268&toegang=8f121ce07d74717e0b1f21d122e04521&id=3641.2012.

[8] Al-Fattah,S M，Startzman,R A.Forecasting World Natural Gas Supply[R].SPE62580，2000：1-7.

[9] Aleklett,et al.The peak of the oil age-analyzing the world oil production reference scenario in world energy outlook 2008[J].Energy Policy,2010,38(5):1398-1414.

[10] Bakhtiari A M.World oil production capacity model suggests output peak by 2006-2007.[EB/OL].http:∥www.energybulletin.net/node/147.2004.

[11] Long,et al.Long-Term World Oil Supply Scenarios:The Future Is Neither as Bleak or Rosy as Some Assert.[EB/OL].http:∥www.eia.doe.gov/pub/oil_gas/petroleum/feature_articles/2004/worldoilsupply/oilsupply04.html.2004.

[12] 陳元千.廣義翁氏預(yù)測(cè)模型的推導(dǎo)與應(yīng)用[J].天然氣工業(yè)，1996,16(2):22- 26.

[13] AlHusseini,The debate over Hubbert’s Peak:a review[J].GeoArabia,2006,11(2):181-210.

[14] Chih-Ming Tien,et al.Quantifying the Uncertainty in Estimates of Ultimately Recoverable World Conventional Oil Resources[J].SPE Economics & Management.2011,3(2):79-92.

[15] Pedro de Almeida,The Peak of Oil Production-Timings and Market Recognition[J].Energy Policy.2009,37(7):1267-127.

[16] Laherrere,Oil & gas production forecasts 1900-2100[EB/OL].http:∥aspofrance.viabloga.com/files/JL_2013_oilgasprodforecasts.pdf.