陸 江，鄧孝亮

(中海石油(中國)有限公司湛江分公司，廣東 湛江 524057)

0 引 言

油氣資源評價是油氣勘探的一個重要環(huán)節(jié)，也是勘探規(guī)劃編制和勘探部署的基礎(chǔ)，因此，對油氣資源的準(zhǔn)確預(yù)測就顯得格外重要。油氣資源分布受多種因素的影響，其評價方法也很多，歸納起來，主要可以分為成因法、類比法和統(tǒng)計法三大類[1-4]。成因法在計算資源量時主要依據(jù)油氣生、運、聚、散的原理，其優(yōu)點是考慮了油氣生成、運移、聚集和保存的地質(zhì)過程，便于把握區(qū)域整體的勘探潛力和方向。但其也有明顯的缺點:一是油氣的生排烴機理目前還存在爭論，不同理論模型計算的結(jié)果差別較大；此外，排聚系數(shù)的定量評價及參數(shù)的選取也具有較大困難[3-8]。類比法主要根據(jù)未知區(qū)與高勘探程度地區(qū)(刻度區(qū))的相似性來選取體積豐度、面積豐度等類比參數(shù)，進行資源量的計算，具有思路簡單、評價快速等優(yōu)點。但一般盆地的地質(zhì)條件具有較大的差異，即很難找到完全類似的類比盆地或區(qū)域來建立刻度區(qū)與未知區(qū)的定量關(guān)系。因此，該方法僅適用于大區(qū)域勘探早期的概略資源評估[7-10]。統(tǒng)計法是利用統(tǒng)計學(xué)的基本原理和方法，根據(jù)高勘探程度地區(qū)油氣田的分布規(guī)律和儲量發(fā)現(xiàn)規(guī)律建立數(shù)學(xué)模型進行油氣資源評價的方法。其優(yōu)點是主要應(yīng)用研究區(qū)的資料，并且不考慮油氣生、排、運、聚、散的過程，回避了其他方法資源評價過程中參數(shù)無法定量評價的問題，同時能獲得評價區(qū)的油氣藏數(shù)量及規(guī)模，結(jié)果較為直觀[3-4,9-11]。在盆地的勘探部署中，剩余油氣資源是其基礎(chǔ)，對其進行準(zhǔn)確的預(yù)測是勘探研究的關(guān)鍵，并且成因法等在計算資源量時參數(shù)的選取多受主觀因素的影響，且其資源量多為靜態(tài)的，很難直接運用于勘探部署中。因此，為了能夠準(zhǔn)確預(yù)測勘探過程剩余油氣資源，基于勘探過程的統(tǒng)計法即是勘探部署研究的最優(yōu)選擇。

1 南海北部灣盆地北部坳陷石油地質(zhì)概況

南海北部灣盆地北部坳陷是南海大陸架重要的產(chǎn)油區(qū)，也是增儲上產(chǎn)的重要領(lǐng)域。坳陷位于中國南海大陸架中部海域，面積約為4.6×104km2，可劃分為多個凹陷和凸起，其中潿西南和海中凹陷是主要的生烴凹陷。

北部坳陷是一個以新生代沉積為主的裂谷盆地，其形成、演化過程復(fù)雜，經(jīng)歷了斷陷期—斷拗期—拗?jǐn)嗥凇窒萜?個階段。對應(yīng)于構(gòu)造演化在縱向上形成了“斷陷湖—拗?jǐn)嗪獮I海潮坪—淺?！彼拇蟪练e層序疊置，主要發(fā)育古近系、新近系以及第四系地層。其中古近系發(fā)育的湖相泥為主要的烴源巖，其上的新近系地層是主要的含油層系。

北部坳陷經(jīng)歷了30多年的勘探，在潿西南凹陷及周緣發(fā)現(xiàn)了一批大中型油田，且勘探程度較高。截至2015年年底，共鉆探井100多口，已探明原油地質(zhì)儲量超過1×108m3。雖然以前針對北部坳陷開展過多次的資源評價工作，但主要局限于成因法(主要是盆地模擬技術(shù))，幾乎沒有涉及統(tǒng)計法，因此，選用勘探效益法和分形模型2種非線性統(tǒng)計法從不同角度對其資源潛力進行評價。

2 勘探效益法和分形模型基本原理

2.1 勘探效益法

一個盆地(或地區(qū))隨著勘探的進展和投入的增加，其油氣資源將會不斷地被發(fā)現(xiàn)(探明)，并且隨著時間的推進，探明儲量將會在達到1個(或幾個)峰值后逐漸降低[1]。這種勘探工作量(指勘探過程中的鉆井?dāng)?shù)或物探工作量等)與探明儲量的關(guān)系可以表示為該探區(qū)的勘探效益，與該探區(qū)油藏規(guī)模與分布、油氣豐度以及勘探程度等因素有關(guān)。勘探效益法即利用累計探明儲量與累計勘探投入(可簡化為探井?dāng)?shù)或進尺)建立統(tǒng)計關(guān)系對待探明資源量(或待探明儲量)進行預(yù)測的方法[4,12]。由于探區(qū)勘探效益的好壞主要受勘探技術(shù)、地質(zhì)認(rèn)識程度及勘探程度的影響，而經(jīng)濟條件變化對其沒有太大的直接影響，因此，該方法用單位勘探工作量而不是單位時間或單位投入資金的探明儲量來度量勘探效益。另外，在預(yù)測時可以根據(jù)勘探效益的變化來確定飽和鉆探(即經(jīng)濟下限)時的勘探工作量及其探明儲量。

勘探效益法主要適用于勘探程度較高的地區(qū)。儲量的發(fā)現(xiàn)(探明)一般隨著勘探工作量(探井?dāng)?shù)或進尺)的增加而呈下降趨勢，其探明儲量可用下式表示：

ΔQ=(dQ/dh)Δh

(1)

(2)

式中：Q為累計發(fā)現(xiàn)量(探明儲量)，104m3；Qλ為盆地(探區(qū))當(dāng)前累計發(fā)現(xiàn)量(探明儲量)，104m3；h為累計探井?dāng)?shù)(累計進尺)，口；hλ為累計發(fā)現(xiàn)量達到Qλ時所對應(yīng)的累計探井，口；hν為累計資源量達到Qν時所對應(yīng)的累計探井，口；dQ/dh為單位探井(進尺)發(fā)現(xiàn)率，104m3/口。

在對歷史資料進行擬合過程中，單位探井發(fā)現(xiàn)率y(即dQ/dh)呈下降趨勢，這種趨勢主要有對數(shù)和指數(shù)2種模型。

指數(shù)下降曲線模型可表述為:

y=ea+bh

(3)

對數(shù)下降曲線模型可表述為:

y=a+blnh

(4)

式中：y為單位探井發(fā)現(xiàn)率，104m3/口；a、b為該盆地(探區(qū))的經(jīng)驗常數(shù)，可由歷史數(shù)據(jù)擬合得出。

2.2 分形模型

分形是數(shù)學(xué)家Mandelbrot于1967年在研究英國海岸長度時首次提出的，并創(chuàng)立了分形幾何學(xué)，形成了研究事物分形性質(zhì)及應(yīng)用的理論[13-14]。分形是指組成整體的各個部分以某種方式與整體相似，這種自相似性既可以是完全相同，也可以只是統(tǒng)計相似，并且可以用分形維數(shù)來定量描述[15-16]。

分形理論應(yīng)用很廣，在油氣勘探開發(fā)中，主要用來進行地震信號的濾波、定量解釋斷裂特征、儲層孔滲與非均質(zhì)性的表征、油氣的運聚過程及資源的定量評價等。根據(jù)國內(nèi)外油氣勘探開發(fā)經(jīng)驗，在一個盆地或地區(qū)的勘探過程中，大的油氣藏一般在早期被發(fā)現(xiàn)，中后期主要發(fā)現(xiàn)中小型油氣藏，且小型油氣藏在發(fā)現(xiàn)數(shù)量上占絕大部分。并且，油氣藏的累計數(shù)量與對應(yīng)的儲量規(guī)模之間存在分形關(guān)系(一般是冪函數(shù)關(guān)系)[17-19]。這種分形關(guān)系可以表示為：

(5)

式中：r為油氣田儲量規(guī)模，104m3；N(r)為儲量規(guī)模大于r的油氣藏的累計個數(shù)，個；D為分形維數(shù)；C為該地區(qū)的經(jīng)驗常數(shù)。

在以油氣田累計數(shù)量N(r)為縱坐標(biāo)，儲量規(guī)模r為橫坐標(biāo)的雙對數(shù)坐標(biāo)系下，直線段即為該地區(qū)的分形曲線。其斜率為油氣田儲量規(guī)模的分形維數(shù)，維數(shù)越小，代表油氣田規(guī)模由大變小的速度越慢，即油氣田的規(guī)模越集中。分形曲線與縱坐標(biāo)的交點代表該盆地最終所能發(fā)現(xiàn)油氣藏的總數(shù)，而拐點的實際意義是目前待發(fā)現(xiàn)油藏的最大儲量級別及對應(yīng)的累計個數(shù)。因此，可以利用分形方法來較好地預(yù)測未發(fā)現(xiàn)油氣田的數(shù)量及其儲量[15]。

3 實例應(yīng)用

3.1 勘探效益法預(yù)測

北部坳陷的儲量發(fā)現(xiàn)過程具有較大的波動性，其中，最高峰在1997至1998年，探明儲量超過2 500×104m3，其他年份探明儲量較低，多低于500×104m3。由于年度新增探明儲量受多種因素影響，因此，為了更好地反映儲量發(fā)現(xiàn)規(guī)律，并降低隨機性的干擾，從1985年起，將發(fā)現(xiàn)的儲量按每5 a進行平均，擬合出儲量發(fā)現(xiàn)趨勢曲線(圖1)。由圖1可以看出，曲線主要表現(xiàn)為先平穩(wěn)增長，然后進入發(fā)現(xiàn)高峰期，最后逐漸降低的過程，總體趨勢滿足勘探效益法預(yù)測油氣資源的條件。

1996年以前有多年的勘探停滯期且總體探明儲量較少，而1997和1998年是儲量發(fā)現(xiàn)的高峰期，其后開始下降。為了研究該坳陷的勘探前景，在探明儲量發(fā)現(xiàn)的平穩(wěn)期對其勘探效益進行分析，以此預(yù)測該坳陷的待探明儲量。首先對2000年及以后各年的累計探井?dāng)?shù)和探明儲量進行統(tǒng)計，并計算出歷年單井探明儲量(圖2)；然后利用單井探明儲量的變化趨勢對今后的探明儲量進行預(yù)測。通過對累計探井?dāng)?shù)與單井探明儲量進行擬合，可以得到單位發(fā)現(xiàn)量的對數(shù)下降模型(圖3)，求得參數(shù)a、b分別為550.71和-90.73。

圖1 南海北部灣盆地北部坳陷探明儲量及趨勢擬合曲線

圖2 南海北部灣盆地北部坳陷原油勘探效益統(tǒng)計曲線

圖3 南海北部灣盆地北部坳陷單井發(fā)現(xiàn)率預(yù)測趨勢

假定最低單井探明儲量取50×104m3/口，則按照預(yù)測趨勢計算出累計探井?dāng)?shù)為249口，計算出其可探明儲量為 25 173.5×104m3。即南海北部灣盆地北部坳陷原油預(yù)計總可探明儲量約為2.5×108m3，預(yù)計待探明儲量約為1.2×108m3。

3.2 分形模型預(yù)測

前已述及，一個盆地油氣藏累計數(shù)量與儲量級別之間具有分形關(guān)系，因此，首先根據(jù)已發(fā)現(xiàn)油藏的探明儲量將油藏的儲量規(guī)模劃分為9個級別(表1)。

表1 南海北部灣盆地北部坳陷油藏儲量規(guī)模分級

在已發(fā)現(xiàn)油藏儲量分級的基礎(chǔ)上，統(tǒng)計各級別以上累計油藏數(shù)量。然后以盆地油藏儲量規(guī)模r作為橫坐標(biāo)、油藏累計數(shù)量N(r≥r1)作為縱坐標(biāo)，在雙對數(shù)坐標(biāo)下研究儲量級別與累計油藏數(shù)量的關(guān)系(圖4)。從圖4可以看出，研究區(qū)在已發(fā)現(xiàn)油藏中，儲量規(guī)模較大的油藏其累計數(shù)量與儲量規(guī)模之間在雙對數(shù)坐標(biāo)下存在較好的線性關(guān)系，即分形關(guān)系(圖4中直線為分形曲線)，該直線的斜率即為該地區(qū)油藏儲量分布的分形維數(shù)(1.161)。圖4中分形曲線向上延伸，規(guī)模較小的油藏儲量分布逐漸偏離了分形曲線，其分離的點探明儲量為800.0×104m3，即認(rèn)為研究區(qū)探明儲量大于800.0×104m3的油藏大部分均已發(fā)現(xiàn)。已發(fā)現(xiàn)油藏分布與分形曲線圍成的區(qū)域為未發(fā)現(xiàn)油藏數(shù)量及規(guī)模。經(jīng)計算，南海北部北部灣盆地北部坳陷儲量規(guī)模大于50×104m3的油藏還有120個，總可探明儲量為3.0×108m3，未探明儲量為1.6×108m3。

綜上所述，南海北部灣盆地北部坳陷已發(fā)現(xiàn)油藏數(shù)量及儲量規(guī)模符合分形分布，可以運用勘探效益法和分形模型對未探明儲量及未探明油藏的規(guī)模及數(shù)量進行預(yù)測?？碧叫б娣A(yù)測其原油預(yù)計總可探明儲量為2.5×108m3，預(yù)計待探明儲量約為1.2×108m3；分形模型預(yù)測南海北部灣盆地北部坳陷儲量規(guī)模大于50×104m3的油藏還有120個，預(yù)計總可探明儲量為3.0×108m3，預(yù)計待探明儲量1.6×108m3。2種方法預(yù)測結(jié)果雖有差異，但也合理，分析認(rèn)為主要有以下原因：①下限的選擇標(biāo)準(zhǔn)不同，即單井探明儲量50.0×104m3與油藏規(guī)模50.0×104m3并不完全對等；②南海北部灣盆地北部坳陷近年勘探效益偏低，單井探明儲量下降變快，使得勘探效益法預(yù)測結(jié)果偏低，而分形模型則主要受已發(fā)現(xiàn)油藏數(shù)量與規(guī)模的影響，與各年的勘探成效關(guān)系不大。

圖4累計油藏數(shù)量與儲量規(guī)模關(guān)系曲線

4 結(jié) 論

(1) 勘探效益法和分形模型等非線性方法在資源預(yù)測過程中不需要考慮油氣資源生、排、運、聚的過程，減小了評價過程中參數(shù)選擇時主觀因素，得出的結(jié)果更為客觀，在南海油氣資源的預(yù)測中具有很好的效果。

(2) 通過非線性方法進行分析得出，南海北部灣盆地北部坳陷預(yù)計待探明原油約為1.4×108m3，仍具有較大的勘探潛力，為勘探部署與規(guī)劃提供了地質(zhì)依據(jù)。

(3) 油氣資源的評價影響因素較為復(fù)雜，運用單一的方法進行預(yù)測可能由于方法的適用性、地質(zhì)認(rèn)識以及勘探現(xiàn)狀等的影響而使結(jié)果誤差較大，因此，在實際工作中應(yīng)該盡量選用多種方法從不同的角度進行預(yù)測，用以相互驗證，減小誤差，提高預(yù)測結(jié)果的可行性。

[1] 周總瑛.統(tǒng)計法在石油資源定量評價中的應(yīng)用——以江漢盆地潛江凹陷為例[J].石油實驗地質(zhì)，2007，29(2)：207-211.

[2] WHITE D A，GEHMAN H M. Methods of estimating oil and gas resources[J]. AAPG Bulletin，1979，63(12)：2183-2192.

[3] 張林曄，李政，孔樣星，等.成熟探區(qū)油氣資源評價方法研究——以渤海灣盆地牛莊洼陷為例[J].天然氣地球科學(xué)，2014，25(4)：477-489.

[4] 祝厚勤，龐雄奇，姜振學(xué)，等.勘探效益法預(yù)測油氣資源的原理及應(yīng)用[J].地質(zhì)科技情報，2006，25(5)：75-79.

[5] 廖前進，于學(xué)敏，何詠梅，等.統(tǒng)計法在中—高成熟探區(qū)資源評價中的適應(yīng)性研究[J].石油學(xué)報，2005，26(增刊)：102-105.

[6] 張寬.中國近海油氣資源評價述評及評價方法探討[J].中國海上油氣(地質(zhì))，2001，15(4)：229-235.

[7] 趙文智，胡素云，沈成喜，等.油氣資源評價的總體思路和方法體系[J].石油學(xué)報，2005，26(增刊)：12-18.

[8] 趙文智，胡素云.油氣資源-儲量分類體系探討[J].石油學(xué)報，2005，26(增刊)：7-11.

[9] 祝厚勤，劉平蘭，龐雄奇.勘探效益法在蘇北盆地金湖凹陷油氣資源評價中的應(yīng)用[J].天然氣地球科學(xué)，2008，19(3)：362-366.

[10] 張寬，宮少波，胡根成.中國近海第三輪油氣資源評價方法述評[J].中國海上油氣，2004，16(4)：218-221.

[11] 郭元嶺，高磊，趙樂強，等.濟陽坳陷石油勘探效果分析[J].河南石油，2001，15(3)：15 -17.

[12] 張厚和，陳蓉.油氣儲量增長趨勢預(yù)測方法及其在中國近海油氣評價中的應(yīng)用[J].天然氣地球科學(xué)，2007，18(5)：684-688.

[13] 郭秋麟，謝紅兵，米石云，等.油氣資源分布的分形特征及應(yīng)用[J].石油學(xué)報，2009，30(3)：379-385.

[14] MANDELBROT B B.How long is the coast of Britain？Statisticalself-similarity and fractional dimension [J ]. Science，1967，156(3775)：636-638.

[15] 陳新，王緒龍，靳濤.石油儲量分布的分形特征及其預(yù)測[J].新疆地質(zhì)，2001，19(4)：297-299.

[16] 劉曉東，徐景禎，王興濤.分形方法預(yù)測氣田數(shù)量及其儲量[J].石油學(xué)報，2000，21(2)：42-44.

[17] 曾怡.分形法預(yù)測油氣儲量與資源量[J].石油實驗地質(zhì)，1998，20(2)：152-154.

[18] POON D C，MCCORMACK M，THIMM H F.The application of fractalgeostatistics to oil and gas property evaluation and reserve estimates[J].The Journal of Canadian Petroleum Technology，1993，32(10)：24-27.

[19] 宋寧，王鐵冠，劉東鷹，等.分形方法在蘇北盆地金湖凹陷石油資源評價中的應(yīng)用[J].地質(zhì)科學(xué)，2006，41(4)：578-585.