多元線性回歸法在DX北帶砂礫巖儲層孔隙度定量預(yù)測中的應(yīng)用

杜 波 于正軍 韓建軍

(1.中國石化勝利油田公司物探研究院 2.山東省東營市住建委)

DX北帶受基底陳南大斷裂的控制，斷坡陡傾，呈溝梁相間的格局，發(fā)育了YA16、YA18和YN93三大古沖溝體系。溝系廣泛發(fā)育近岸水下扇、扇三角洲、深水濁積扇、滑塌濁積扇等多種類型砂礫巖扇體，這些砂礫巖扇體沿陡坡帶呈裙帶狀分布，縱向上多期疊置，橫向上疊合連片，具有溝扇對應(yīng)、大溝對大扇的特點。扇三角洲前緣亞相、近岸水下扇扇中亞相和前緣滑塌扇的砂礫巖體具有較好的儲層物性，能夠形成良好儲層。陡坡帶砂礫巖體埋藏深、鄰近多個生油洼陷，具有良好的油源條件。目前的勘探認(rèn)為，沙四上亞段和沙三段具有構(gòu)造背景的砂礫巖具有最好的成藏條件。進(jìn)行儲層孔隙度預(yù)測對進(jìn)一步推進(jìn)砂礫巖勘探具有重要意義。

1 技術(shù)流程

根據(jù)砂礫巖體的沉積特點，并考慮到實際操作的可行性，砂礫巖體儲層孔隙度平面預(yù)測的技術(shù)流程包括：①在礫巖體期次劃分的基礎(chǔ)上，進(jìn)行精細(xì)合成地震記錄標(biāo)定以確定地震分層與井分層的對應(yīng)關(guān)系；②從井資料出發(fā)，統(tǒng)計各期次砂礫巖體的孔隙度，同時依托地震數(shù)據(jù)體，提取各期次砂礫巖體的多種地震屬性，通過從屬性物理意義和相關(guān)性回歸分析確定所用地震屬性的類型和個數(shù)；③建立地震屬性與孔隙度之間的函數(shù)關(guān)系，通過函數(shù)關(guān)系把地震屬性轉(zhuǎn)換為孔隙度的一種表達(dá)方式，獲得期次砂礫巖體儲層孔隙度的平面預(yù)測結(jié)果(圖1)。

圖1 砂礫巖體儲層孔隙度平面預(yù)測技術(shù)流程圖

2 應(yīng)用條件

嚴(yán)密控制數(shù)據(jù)參數(shù)質(zhì)量，充分滿足多元回歸法應(yīng)用條件是獲得良好預(yù)測效果的基礎(chǔ)[1-3]。

2.1 砂礫巖體期次儲層孔隙度統(tǒng)計

研究工區(qū)內(nèi)有50余口探井，大部分都有聲波測井曲線，但可收集到孔隙度成果曲線的較少。為了準(zhǔn)確求得孔隙度平面預(yù)測的函數(shù)模型，需要有盡可能多的井點孔隙度數(shù)據(jù)，因此，可以利用聲波測井曲線間接獲取更多井點孔隙度數(shù)據(jù)。選取具有孔隙度曲線和聲波曲線的YA22井、YN937井、YA162井、YN920井、YN935井中的典型砂礫巖體層段進(jìn)行統(tǒng)計擬合，求得各井砂礫巖體聲波速度與孔隙度的函數(shù)關(guān)系式。統(tǒng)計擬合表明砂礫巖體孔隙度與聲波速度呈較簡單的線性關(guān)系。由于測井年代和測井條件的不同，各井?dāng)M合函數(shù)系數(shù)有所差異，為了統(tǒng)一儲層孔隙度的統(tǒng)計對比條件，獲得區(qū)域上普遍適用的函數(shù)關(guān)系式，在直角坐標(biāo)系中用作圖法對各函數(shù)式進(jìn)一步擬合獲得新的函數(shù)式：

y=0.4604x-21.643

(1)

式中：

y—孔隙度；

x—聲波速度。

以上函數(shù)的數(shù)據(jù)源涵蓋了YA16井、YA18井之間兩大古沖溝的砂礫巖主扇體的扇根、扇中和扇間部位，深度為1000m～4500m，基本覆蓋了現(xiàn)今礫巖勘探深度范圍，因而數(shù)據(jù)統(tǒng)計具有典型代表性。用(1)式可方便快捷地從砂礫巖聲波測井中求取更多孔隙度數(shù)據(jù)以滿足預(yù)測的需要。

2.2 地震層位解釋和屬性提取

陡坡帶砂礫巖扇體沉積復(fù)雜，并受后期構(gòu)造改造，上下層位之間產(chǎn)狀迥異，單個層位控制時窗難以準(zhǔn)確獲取砂礫巖期次的地震屬性信息，因此使用上下兩個層位界面來界定提取屬性的時窗是正確的選擇。地震期次層位解釋時根據(jù)砂礫巖的沉積特點和提取屬性的要求，在陡坡凸起方向沿基巖面解釋至上下層位齊平，在洼陷方向擴大解釋范圍，以利于宏觀比較、提高預(yù)測結(jié)果與地質(zhì)規(guī)律的吻合程度。據(jù)此，DX北帶自上而下解釋砂礫巖扇體12期次，其中沙三段4個期次，沙四上3個期次，沙四下5個期次(圖2)。以每一期次上下兩個層位界面來界定屬性時窗，然后利用SATK軟件分別提取各個期次振幅能量、頻率、相位、波形等方面20多種地震屬性。

圖2 DX北帶砂礫巖體期次解釋成果剖面圖

2.3 地震屬性優(yōu)選

(1)地震屬性交匯相關(guān)性分析

地震屬性相關(guān)性分析是確定屬性間的獨立性以及地震屬性對儲層參數(shù)反映敏感程度的重要工具，它的主要任務(wù)是計算地震屬性參數(shù)間以及地震屬性參數(shù)與儲層參數(shù)間的相關(guān)性，并根據(jù)相關(guān)性大小確定相對獨立的、能夠反映儲層參數(shù)本質(zhì)特征的地震屬性參數(shù)[4-6]。

DX北帶第8期次砂礫巖儲層孔隙度與22種常用地震屬性的相關(guān)性分析認(rèn)為，零點相位數(shù)、負(fù)振幅總和、強度總和等屬性與砂礫巖孔隙度的相關(guān)性較大，而平均振幅、主頻率、瞬時相位等屬性與砂礫巖孔隙度相關(guān)性很小。該期次有12個井點孔隙度值，相關(guān)性大于0.4的屬性有8種，根據(jù)統(tǒng)計學(xué)定義認(rèn)為其具有顯著性意義，表明這些地震屬性與砂礫巖孔隙度之間存在著某種內(nèi)在的聯(lián)系。

進(jìn)一步考察此8種地震屬性之間的獨立性。交匯相關(guān)性分析(表1)表明，帶寬比率偏差與帶寬比率除偏、負(fù)振幅總和與正振幅總和之間具有很好的相關(guān)性，相關(guān)系數(shù)分別為0.996、0.702，表明這兩組屬性之間有很大的相似性，選擇其一即可。

表1 8種地震屬性交匯相關(guān)系數(shù)統(tǒng)計表

(2)地震屬性物理意義分析

雖然交匯相關(guān)性分析大大地減小了待選地震屬性范圍，但由于勘探中獲得孔隙度數(shù)據(jù)有限，以及孔隙度統(tǒng)計、地震屬性提取中不可避免地存在誤差，使得某些本身無物理意義或者不相干的地震屬性也顯示較好的相關(guān)性，而使用這些屬性做孔隙度預(yù)測顯然不合適。純數(shù)學(xué)的相關(guān)性檢驗方法不能消除這種誤差的影響[4-7]，因此，在待選屬性不多的情況下，通過孔隙度預(yù)測原理和地震屬性物理意義的分析來確定最優(yōu)屬性組合是最好的選擇。

地震屬性中振幅能量與儲層孔隙度關(guān)系是最密切的，但是地震振幅能量還與儲層厚度有關(guān)。根據(jù)Widess調(diào)諧原理，當(dāng)?shù)貙雍穸融呄蛴谡{(diào)諧厚度1/4波長(λ)時具有最強的振幅能量，理論上調(diào)諧振幅能量等于儲層上下界面振幅能量的疊加。因此，基于砂礫巖體期次的儲層孔隙度平面預(yù)測需要主要考慮的問題是如何確定儲層孔隙度和厚度對期次振幅能量的貢獻(xiàn)量。為此，取泥巖速度110μs/ft、標(biāo)準(zhǔn)砂礫巖體速度60μs/ft及30Hz雷克子波建立了不同厚度和速度的正極性一維地震正演模擬模型(圖3)。從圖3中可以看出，儲層上界面呈正極性反射是砂礫巖體的反射，儲層下界面負(fù)極性反射是泥巖的反射；當(dāng)薄層儲層厚度為(1/4～1/8)λ(即42m～21m)時，上下界面的地震反射波產(chǎn)生疊加而有最大的振幅能量，這就是儲層厚度的調(diào)諧效應(yīng)，因其影響，厚度低至λ/32(5m)的薄儲層仍有較強的能量反射。

注：泥巖速度110μs/ft

另一方面，儲層速度的提高會導(dǎo)致地層反射系數(shù)的變大，從而使地震反射增強。假設(shè)厚層儲層速度(波阻抗)需要提高到x倍時與薄層調(diào)諧振幅的能量相當(dāng)，則有

(p2-p1)/(p2+p1) ：(xp2-p1)/(xp2+p1)=1∶2

x=3/(2-p2/p1)>1.5

(2)

式中：

p1、p2—波阻抗，p1=ρ1v1，p2=ρ2v2。

很明顯，厚層儲層速度提高到1.5倍時與薄層調(diào)諧振幅的能量相當(dāng)。從圖3中也可以看出，厚層儲層速度提高到40μs/ft時與60μs/ft的薄層調(diào)諧振幅的能量相當(dāng)，證實了理論推導(dǎo)的正確性。砂礫巖體的聲波速度為50μs/ft ～80μs/ft，上述的速度變化幅度并不明顯，所以，單憑速度提高對增強地震反射的貢獻(xiàn)是有限的。此外，砂礫巖體具有強烈的非均質(zhì)性，大套的礫巖體實際上是由許多物理性質(zhì)不同的小層構(gòu)成，從這個意義上講，砂礫巖體期次的振幅能量主要來自多個薄層儲層產(chǎn)生的調(diào)諧效應(yīng)。因此，利用振幅能量預(yù)測砂礫巖體儲層孔隙度應(yīng)選擇那些能反映期次內(nèi)反射界面?zhèn)€數(shù)以及組合起來能消除調(diào)諧效應(yīng)的屬性，即強度總和、負(fù)振幅總和、零點相位數(shù)。其中，強度總和反映了時窗內(nèi)振幅絕對值的總和，零點相位數(shù)表示地震波零相位的個數(shù)，反映了地層反射界面?zhèn)€數(shù)。在正極性地震正演模擬中，這3種屬性變量能完整地描述砂礫巖體期次振幅能量變化關(guān)系，通過回歸擬合求得的函數(shù)式可以有效消除砂礫巖體期次內(nèi)薄層調(diào)諧效應(yīng)的影響。

2.4 數(shù)據(jù)異常值分析

由于預(yù)測方法中需要的各種數(shù)據(jù)不可避免地存在誤差，有時這些誤差比較大，會嚴(yán)重影響預(yù)測效果，因此，在擬合回歸函數(shù)時應(yīng)當(dāng)盡可能消除這些誤差的影響[7]。通過離散數(shù)據(jù)交匯趨勢性分析可以識別某些偶然誤差的存在。

圖4是第8期次砂礫巖體孔隙度與強度總和、負(fù)振幅總和、零點相位數(shù)三種地震屬性參數(shù)的交匯圖，在誤差校正之前孔隙度與這3種地震屬性的相關(guān)性系數(shù)分別為53.1%、54%和61.9%，通過觀察交匯圖的樣點分布情況可以看出，10、11兩點(空心圓)遠(yuǎn)離樣點總體趨勢，因此有理由相信這兩個點是存在較大誤差的異常值，可以通過校正統(tǒng)計值或者予以剔除來消除其對預(yù)測的影響。剔除兩個異常點后，孔隙度與地震屬性的相關(guān)性分別提高到70.7%、70.9%和73.1%(圖4)。

圖4 第8期次砂礫巖體孔隙度與3種地震屬性相關(guān)性的散點交匯圖

3 孔隙度預(yù)測效果

在確定地震屬性最優(yōu)組合并剔除異常值提高相關(guān)性之后，利用LPM軟件可自動求得多元回歸方程。第8期次砂礫巖體儲層孔隙度(Ф)的多元回歸方程為

Ф= 1.98991+0.000110721x1+0.00015185x2

-0.294354x3

(3)

式中：

x1—強度總和；

x2—負(fù)振幅總和；

x3—零點相位數(shù)。

各砂礫巖體期次深度、地層反射界面?zhèn)€數(shù)不同，其擬合函數(shù)的各屬性變量組合形式不同，表現(xiàn)為多元回歸方程系數(shù)的變化，如第5期次多元回歸方程為：

Ф=-7.02295+0.0000784052x1+0.000085249x2

+0.675982x3

(4)

孔隙度平面預(yù)測結(jié)果(圖5)清晰地反映了DX北帶砂礫巖扇體的展布特征：砂礫巖扇體沿溝谷分布，在古沖溝前方扇中的部位具有較高的孔隙度，而在扇根、扇端以及古沖溝交匯的部位孔隙度較低，與砂礫巖扇體地質(zhì)沉積規(guī)律是相吻合的；孔隙度預(yù)測值與實測值誤差為-1.272～1.010，平均絕對誤差僅為0.351(表2)，說明預(yù)測結(jié)果也是比較準(zhǔn)確的。

4 結(jié)論與認(rèn)識

(1)多元線性回歸法預(yù)測砂礫巖體儲層平面孔隙度比較實用，易實現(xiàn)。從孔隙度統(tǒng)計、層位解釋、屬性優(yōu)化分析等方面充分滿足預(yù)測方法的應(yīng)用條件是獲得良好預(yù)測效果的基礎(chǔ)。

表2 DX北帶兩個期次砂礫巖體孔隙度預(yù)測值與實測值誤差

(2)多元線性回歸預(yù)測儲層孔隙度的關(guān)鍵問題是地震屬性的優(yōu)選，需要選擇反映儲層參數(shù)敏感的、個數(shù)最少、相互獨立的地震屬性。地震屬性相關(guān)性和物理意義分析結(jié)果表明，強度總和、負(fù)振幅總和、零點相位數(shù)是最優(yōu)地震屬性組合。

(3)受多種因素影響，數(shù)據(jù)樣點不可避免地存在奇異值。通過離散數(shù)據(jù)交匯趨勢性分析可有效識別并去除異常值的影響，從而有助于提高孔隙度預(yù)測精度。

(4)使用多元線性回歸法預(yù)測的DX北帶砂礫巖體儲層平面孔隙度特征符合地質(zhì)規(guī)律，預(yù)測精度較高。

1 樂有喜，王永剛，張軍華，等.儲層參數(shù)平面分布預(yù)測方法評價[J].地質(zhì)與勘探，2001，37((5)：56-60.

2 張麗艷.砂礫巖儲層孔隙度和滲透率預(yù)測方法[J].測井技術(shù)，2005，29(3)：212-215.

3 顧偉欣，周紅.基于逐步回歸分析的孔隙度預(yù)測方法[J].石油地質(zhì)與工程，2008，22(1)：37-39.

4 周宗良，肖建玲，張楓，等.地震屬性的優(yōu)化處理及儲層厚度的定量解釋[J].新疆地質(zhì)，2002，20(3)：262-265.

5 王治國，尹成，雷小蘭，等.預(yù)測砂巖孔隙度的地震多屬性優(yōu)化模式對比[J].石油地球物理勘探，2011，46(3)：442-448.

6 張新亮，何麗菁，吳俊.基于地震多屬性的孔隙度預(yù)測——以川東A氣田為例[J].新疆石油地質(zhì)，2011，(4)：385-386.

7 于志鈞，趙旭東.石油數(shù)學(xué)地質(zhì)[M].北京：石油工業(yè)出版社，1986.