王艷妮 劉 剛 潘源敦 王獻(xiàn)勇

(1.江蘇省有色金屬地質(zhì)勘查局地質(zhì)信息中心 2.中國地質(zhì)大學(xué)(武漢)計(jì)算機(jī)學(xué)院3.教育部“構(gòu)造與油氣資源”重點(diǎn)實(shí)驗(yàn)室 4.中國石化荊州勘查技術(shù)研究中心)

特爾菲與蒙特卡羅聯(lián)合方法在油氣資源預(yù)測評(píng)價(jià)中的應(yīng)用

王艷妮1劉 剛2，3潘源敦4王獻(xiàn)勇2

(1.江蘇省有色金屬地質(zhì)勘查局地質(zhì)信息中心 2.中國地質(zhì)大學(xué)(武漢)計(jì)算機(jī)學(xué)院3.教育部“構(gòu)造與油氣資源”重點(diǎn)實(shí)驗(yàn)室 4.中國石化荊州勘查技術(shù)研究中心)

由于油氣資源成藏條件復(fù)雜、影響因素眾多，其預(yù)測評(píng)價(jià)存在很大的不確定性。傳統(tǒng)的方法估算儲(chǔ)量只能提供一個(gè)確定的儲(chǔ)量數(shù)值，不能定量說明該數(shù)值的可靠程度。為了更加精確和合理地估算油氣資源的儲(chǔ)量，將特爾菲法和蒙特卡羅法結(jié)合起來對油氣資源的儲(chǔ)量進(jìn)行預(yù)測評(píng)價(jià)，特爾菲法綜合各影響因子，蒙特卡羅法實(shí)際運(yùn)算三角分布數(shù)據(jù)、觀測數(shù)據(jù)或者經(jīng)驗(yàn)估值，計(jì)算的結(jié)果是一條儲(chǔ)量逆累積概率曲線。在使用儲(chǔ)量數(shù)值時(shí)，選用一定概率條件下的一個(gè)儲(chǔ)量范圍值，這樣更能反映人們對地下資源的認(rèn)識(shí)程度，因而對勘探開發(fā)工作更富有指導(dǎo)意義。圖3表3參6

油氣資源評(píng)價(jià) 特爾菲法 蒙特卡羅法 逆累積概率

0 引言

油氣資源定量評(píng)價(jià)是制定國家能源政策與編制油氣發(fā)展規(guī)劃的重要依據(jù)。在油氣田勘探開發(fā)工作中，對油氣儲(chǔ)量和油氣資源量的預(yù)測與估算是一項(xiàng)十分重要的工作。在油氣田勘探初期，為了取得某一區(qū)塊的勘探風(fēng)險(xiǎn)投資，通常需要提供該區(qū)塊的儲(chǔ)量和了解其可靠程度，即得到該儲(chǔ)量可能出現(xiàn)的概率。

針對油氣資源預(yù)測評(píng)價(jià)因數(shù)復(fù)雜、參數(shù)繁多的特點(diǎn)，本文提出用特爾菲法［1］對各油氣資源評(píng)價(jià)因子進(jìn)行預(yù)測估值，綜合各個(gè)專家的評(píng)價(jià)意見后，得到整個(gè)區(qū)間段上各個(gè)因子的估計(jì)值及其逆累計(jì)概率，然后將各因子預(yù)測結(jié)果參與蒙特卡羅法［2］的計(jì)算，完成油氣資源的預(yù)測評(píng)價(jià)過程，并給出評(píng)價(jià)結(jié)果的累積概率值及其分布曲線以輔助支持決策。

1 油氣資源預(yù)測評(píng)價(jià)過程

基于地質(zhì)類比模型和概率模型相結(jié)合的方式，采用特爾菲和蒙特卡羅計(jì)算方法完成油氣資源的預(yù)測評(píng)價(jià)過程。技術(shù)流程如圖1所示:

將特爾菲與概率地質(zhì)類比模型相結(jié)合完成多個(gè)地質(zhì)專家意見的綜合，以提高結(jié)果的可信度。在預(yù)測評(píng)價(jià)過程中，特爾菲法的一項(xiàng)主要工作，是每輪征詢之后的數(shù)據(jù)分析和處理。本文中主要采用三角分布抽樣方法對數(shù)據(jù)進(jìn)行處理和分析，具體計(jì)算過程如下:

圖1 油氣資源預(yù)測評(píng)價(jià)技術(shù)流程圖

將整個(gè)區(qū)間段劃分為m個(gè)小區(qū)間，求出同一油氣資源評(píng)價(jià)因子所有專家的所有估算值Qij中的極大值Qmax和極小值Qmin，作為最終評(píng)價(jià)的區(qū)間估計(jì)值的端點(diǎn)。進(jìn)而求出(m+1)個(gè)小區(qū)間的分界值Qk。

完成特爾菲法的第一輪計(jì)算后，若評(píng)價(jià)結(jié)果不滿足收斂要求，則把第一輪的評(píng)價(jià)結(jié)果發(fā)給各個(gè)專家，各個(gè)專家參考上一輪的結(jié)果修改自己的意見，進(jìn)行下一輪評(píng)價(jià)，如此反復(fù)進(jìn)行三次，當(dāng)大多數(shù)油氣資源評(píng)價(jià)因子滿足要求時(shí)，評(píng)價(jià)即告結(jié)束。如此，即可求出所有專家對某一個(gè)油氣資源評(píng)價(jià)因子的估計(jì)值及其逆累計(jì)概率。

在各專家結(jié)論的基礎(chǔ)上，基于蒙特卡羅法完成油氣資源量的計(jì)算和估計(jì)。本文中采用成因法(有機(jī)碳法)計(jì)算油氣資源量。具體如下:

式中:

Q資—油(氣)資源量，108t;

Q生—生油(氣)量，108t;

As—有效烴源巖面積，km2;

Hs—有效烴源巖厚度，m;

γ—烴源巖密度，g/cm3;

C—烴源巖殘余有機(jī)碳含量，%;

K恢復(fù)—有機(jī)碳恢復(fù)系數(shù);

K產(chǎn)烴率—有機(jī)碳產(chǎn)烴率，mg/gToc;

K運(yùn)聚—運(yùn)聚系數(shù)。

對于勘探程度較低的區(qū)域，不易直接取得各項(xiàng)計(jì)算各個(gè)因子，本文均采用概率統(tǒng)計(jì)方法模擬各個(gè)計(jì)算因子，這樣更符合此海域資料的客觀狀況。因此對計(jì)算方法進(jìn)行如下變通:

式中:

Fx—因子x的逆累計(jì)概率分布

2 實(shí)例計(jì)算

編程實(shí)現(xiàn)的軟件界面如圖2所示?？梢钥吹教峁┝硕喾N分布數(shù)據(jù)的處理方式，蒙特卡洛法計(jì)算式抽樣次數(shù)默認(rèn)為2000次，并可根據(jù)實(shí)際情況進(jìn)行參數(shù)調(diào)整。

模擬計(jì)算我國某一含油(氣)單元的天然氣資源量，其中面積和密度兩個(gè)因子是已知的確定值，將其看作常數(shù)。

各專家對研究區(qū)油氣資源評(píng)價(jià)因子—產(chǎn)氣率的評(píng)價(jià)意見匯總，如表1所示。

經(jīng)特爾菲法計(jì)算各個(gè)油氣資源評(píng)價(jià)因子在不同概率下的估算值，如表2所示。

在表2中，厚度因子是通過統(tǒng)計(jì)型分布進(jìn)行模擬其隨機(jī)變量及其函數(shù)分布的，而其余四個(gè)因子通過插值型模擬的。在計(jì)算出了各個(gè)因子的函數(shù)分布之后，進(jìn)行蒙特卡羅的連乘運(yùn)算。求天然氣資源量的總公式為:

圖2 蒙特卡洛計(jì)算法軟件界面

表1 產(chǎn)氣率專家估計(jì)值三角分布表

表2 各個(gè)因子在不同概率下的估算值

表3 不同概率下的天然氣資源量分布

圖3 研究區(qū)天然氣資源量逆累積概率分布曲線

為了進(jìn)行對比，若不用聘請專家，直接用蒙特卡羅法計(jì)算儲(chǔ)量及其逆累積概率。如圖3所示。

傳統(tǒng)的油氣資源儲(chǔ)量估值一般都是用蒙特卡羅計(jì)算法，直接求出油氣資源的儲(chǔ)量。本文中綜合各個(gè)專家對油氣資源評(píng)價(jià)因子評(píng)價(jià)意見后，再使各油氣評(píng)價(jià)因子參與蒙特卡羅計(jì)算而得到的油氣資源量的逆累積分布曲線圖。分析圖3可看出，所得到的天然氣儲(chǔ)量估值更加精確，范圍更加確切。

3 結(jié)論

特爾菲法提供了對各個(gè)油氣資源評(píng)價(jià)因子的綜合統(tǒng)計(jì)分析功能，其特點(diǎn)是:匿名性、反饋性、統(tǒng)計(jì)型，將特爾菲法和蒙特卡羅計(jì)算方法相結(jié)合，來完成油氣資源的預(yù)測評(píng)價(jià)過程，可以得到在不同概率下的油氣資源量的估計(jì)曲線。比以往的用蒙特卡羅法單獨(dú)估算油氣資源儲(chǔ)量，有更好的使用性，評(píng)價(jià)結(jié)果更加準(zhǔn)確，更有利于油氣資源開發(fā)人員決策支持。

運(yùn)用特爾菲法在選擇專家時(shí)應(yīng)該注意以下兩點(diǎn):①專家的類型，不同類型的專家由于考慮問題的側(cè)重點(diǎn)不同，對同一因素給出的預(yù)測往往是不一樣的。②專家的數(shù)量，專家太少容易片面，太多則數(shù)據(jù)整理工作量大，不易收斂;應(yīng)根據(jù)問題的性質(zhì)、規(guī)模及決策等因素確定［3］。

1 趙旭東，崔鳳文.用于石油資源評(píng)價(jià)的特爾菲法［J］.石油勘探與開發(fā)，1983，3:79－82.

2 王銀宏嚴(yán)光生翟裕生.三部式潛在礦產(chǎn)資源定量評(píng)價(jià)與蒙特卡羅模擬［J］.中國礦業(yè)，2006，15(6):4－7.

3 宋遠(yuǎn)方.層次分析法及特爾菲法在礦區(qū)規(guī)劃決策分析中的應(yīng)用［J］. 煤炭工程，1989，9:12－16.

4 文環(huán)明，肖慈，王旬，等.蒙特卡羅法在油氣儲(chǔ)量估算中的應(yīng)用［J］.成都理工學(xué)院學(xué)報(bào)，2002，29(5):487－492.

5 陳榮書.石油及天然氣地質(zhì)學(xué)［M］.武漢:中國地質(zhì)大學(xué)出版社，1994

6 查先進(jìn).信息分析與預(yù)測［M］.武漢:武漢大學(xué)出版社，2000，8.

APPLICATION OF BOTH DELPHI AND MONTE CARLO METHODS TO EVALUATION OF HYDROCARBON RESOURCES

WANG Yanni1，LIU Gang2，3，PAN Yuandun4and WANG Xianyong2(1．Information Center，Nonferrous Metals Geological Exploration Bureau of Jiangsu Province;2．College of Resources，China University of Geosciences(Wuhan);3．Key Laboratory of Tectonics and Petroleum Resources，Ministry of Education;4．Sinopec Jing Zhou Exploration Technology Research Center)

Since reservoir－forming conditions are not only complicated but also influenced by many factors，to evaluate hydrocarbon resources is very uncertain．Some traditional methods to assess reserves can only provide a certain value，which can not quantitatively demonstrate a reliability of the value．In this paper，in order to more accurately and reasonably assess reserves，both Delphi and Monte Carlo methods are used．During assessment，all influencing factors are integrated by Delphi method whereas the operations of triangular distribution data，observed data or empirical estimation are actually performed by Monte Carlo method．The assessed result is an inverse cumulative－ probability curve of reserves．When using reserves values，a range value of reserves under a given probability condition is selected，which can further reflect an understanding of hydrocarbon resources．

evaluation of hydrocarbon resources evaluation，Delphi method，Monte Carlo method，inverse cumulative－probability

王艷妮，女，1984年出生，陜西韓城人。目前就職于江蘇省有色金屬地質(zhì)勘查局地質(zhì)信息中心，從事油氣資源評(píng)價(jià)及地質(zhì)找礦信息化方面的工作。地址:(210007)江蘇省南京市白下區(qū)大光路26號(hào)華東大廈4樓。電話:13675167460。E－mail:wynhcs@qq.com

．NATURALGAS EXPLORATION＆DEVELOPMENT．v．34，no．3 ，pp．24 －26，7/25/2011

(修改回稿日期 2011－02－23 編輯 王曉清)