王迪，張益明，牛聰，黃饒，韓利

(1.中海油研究總院有限責(zé)任公司，北京 100028；2.中國(guó)海洋石油國(guó)際有限公司，北京 100028)

0 引言

深水濁積砂巖儲(chǔ)層是當(dāng)前世界油氣勘探開(kāi)發(fā)的熱點(diǎn)目標(biāo)，以西非的尼日爾三角洲盆地為典型代表[1-4]。相比陸地和淺水沉積，深水區(qū)油氣勘探開(kāi)發(fā)成本高且風(fēng)險(xiǎn)大，因此，如何準(zhǔn)確預(yù)測(cè)優(yōu)質(zhì)儲(chǔ)層并開(kāi)展烴類(lèi)檢測(cè)對(duì)深水勘探尤為重要。疊前反演是進(jìn)行烴類(lèi)檢測(cè)的有效地球物理技術(shù)手段，其關(guān)鍵是首先優(yōu)選出能夠反映儲(chǔ)層流體性質(zhì)的敏感彈性參數(shù)，然后通過(guò)反演獲得該參數(shù)進(jìn)行流體識(shí)別。早在20世紀(jì)80年代，Ostrander[5]開(kāi)展了基于AVO技術(shù)的流體識(shí)別研究。隨后，Smith等[6]提出了流體因子的概念。Goodway等[7]提出了λρ-μρ交會(huì)分析的烴類(lèi)檢測(cè)方法。Russell等[8]利用縱波阻抗與橫波阻抗構(gòu)建了新的流體指示因子。張玉潔等[9]基于噴流效應(yīng)對(duì)Russell流體因子進(jìn)行了推廣和應(yīng)用。姜仁等[10]利用Russell流體因子開(kāi)展了致密儲(chǔ)層的含氣性檢測(cè)。鄭靜靜等[11]對(duì)各類(lèi)流體因子進(jìn)行了分類(lèi)與系統(tǒng)總結(jié)。張廣智等[12]研究了角度流體因子提取方法，能夠更準(zhǔn)確地識(shí)別目的層處的異常，同時(shí)提高剖面的信噪比。謝玉洪等[13]探索了“暗點(diǎn)”型油氣藏的流體指示因子，可以放大II類(lèi)AVO異常的差異。

通常情況下，流體因子是通過(guò)疊前反演縱波阻抗和橫波阻抗等參數(shù)間接計(jì)算獲得，容易出現(xiàn)累計(jì)誤差，為此諸多學(xué)者提出了流體因子的直接反演方法。宗兆云等[14]提出了包含拉梅系數(shù)的彈性阻抗方程，可直接反演拉梅參數(shù)。印興耀等[15]提出了包含Gassmann流體項(xiàng)的彈性阻抗方程，可直接反演高精度的流體項(xiàng)。李紅梅[16]應(yīng)用疊前彈性參數(shù)直接反演技術(shù)，提取了對(duì)儲(chǔ)層流體敏感的λρ、泊松比和高靈敏度流體因子等參數(shù)，對(duì)KD地區(qū)淺層河道砂巖油藏進(jìn)行了識(shí)別。在流體因子敏感性評(píng)價(jià)方面，楊培杰等[17]建立了一種敏感因子定量分析方法。桂金詠等[18]提出了面向?qū)嶋H儲(chǔ)層的流體因子優(yōu)選方法。

在實(shí)際應(yīng)用中，流體因子往往受到諸多因素、特別是孔隙度的影響，由于忽略孔隙度而出現(xiàn)流體識(shí)別假象的案例越來(lái)越多。張世鑫[19]研究發(fā)現(xiàn)復(fù)雜儲(chǔ)層受孔隙度影響會(huì)產(chǎn)生流體識(shí)別假象，如何將孔隙度的影響從流體因子中剔除顯得至關(guān)重要。李英等[20]通過(guò)流體替代、孔隙度替代等正演模擬顯示，在孔隙度波動(dòng)范圍較大的區(qū)域，流體變化導(dǎo)致的地震異常響應(yīng)往往會(huì)被孔隙度影響所掩蓋，此時(shí)“亮點(diǎn)”的響應(yīng)往往表征的只是優(yōu)質(zhì)儲(chǔ)層而非含油氣層。為了壓制孔隙度的影響，諸多學(xué)者[21-23]指出流體模量只與孔隙內(nèi)流體類(lèi)型有關(guān)而與孔隙度無(wú)關(guān)，具有更好的流體識(shí)別效果，開(kāi)展了流體模量反演方法的研究。由于流體模量反演在實(shí)際應(yīng)用中受限于地震資料品質(zhì)、干巖石模量估算精度等問(wèn)題，結(jié)果存在不確定性。常規(guī)反演方法仍然是目前業(yè)界主流的技術(shù)手段，如何優(yōu)選出靶區(qū)真正體現(xiàn)流體性質(zhì)差異的敏感流體因子、提高常規(guī)反演流體識(shí)別精度顯得至關(guān)重要。為此，本文在綜合分析儲(chǔ)層地震響應(yīng)特征的基礎(chǔ)上，提出了一種新的流體因子敏感性定量分析方法，篩選的流體因子具有“對(duì)流體性質(zhì)敏感度高、對(duì)孔隙度敏感度低”的特征，能夠最大限度地突出流體異常，壓制孔隙度的影響。實(shí)際應(yīng)用結(jié)果表明，該流體因子能夠有效區(qū)分油層和高孔隙水層，與已鉆井吻合率更高。

1 區(qū)塊概況

尼日爾三角洲盆地位于非洲大陸西部，面積達(dá)12×104km2，是當(dāng)今世界熱門(mén)的深水勘探區(qū)域。尼日爾三角洲是始新世以來(lái)，尼日爾河和貝努埃河沿著比達(dá)盆地和貝努埃槽向南匯合流入大西洋，其所攜帶的大量碎屑物質(zhì)和有機(jī)質(zhì)形成巨大的海退型前積三角洲沉積。尼日爾三角洲深水扇則是由尼日爾河和尼日爾三角洲提供物源，在深海平原形成的大型深水扇。本文研究的S區(qū)塊位于尼日利亞海域(圖1)，北距陸上哈科特港200 km，水深為1 300～1 450 m，主力儲(chǔ)層為中新統(tǒng)R1180層系發(fā)育的水道化朵葉沉積體，砂體單層厚度大，物性較好，孔隙度介于22%～30%之間，儲(chǔ)層條件優(yōu)越。研究區(qū)內(nèi)共有3口已鉆井，鉆探結(jié)果如表1所示。其中，W1和W2井鉆遇油層，W3井鉆遇水層。

圖1 研究區(qū)位置示意

表1 研究區(qū)已鉆井結(jié)果

2 儲(chǔ)層地震響應(yīng)特征及影響因素

2.1 疊后振幅響應(yīng)特征及影響因素

圖2所示為過(guò)W1、W2和W3井的連井地震剖面，從圖中可以看到，儲(chǔ)層均表現(xiàn)為振幅“亮點(diǎn)”特征。其中，W3井振幅最強(qiáng)，鉆遇水層，W1和W2井振幅稍弱，卻鉆遇油層。根據(jù)表1可知，W1井砂巖孔隙度為26.0%，W2井砂巖孔隙度為23.4%，W3井砂巖孔隙度為29.7%，推測(cè)物性可能是導(dǎo)致水層出現(xiàn)振幅“亮點(diǎn)”的主要原因。

圖2 過(guò)W1、W2和W3井的連井地震剖面

進(jìn)一步開(kāi)展疊后振幅的影響因素分析，包括流體性質(zhì)、孔隙度和儲(chǔ)層厚度(圖3)。可以看出，當(dāng)儲(chǔ)層孔隙度均為26%時(shí)，油層和水層的振幅存在明顯差異，含油砂巖振幅相比含水砂巖要強(qiáng)。受調(diào)諧作用的影響，油層和水層的振幅均在調(diào)諧厚度處達(dá)到最大值?？紤]到已鉆井揭示的砂巖厚度主要集中在20～25 m，在該范圍內(nèi)厚度變化對(duì)疊后振幅的影響相對(duì)較弱，不是該區(qū)疊后振幅的主要影響因素。當(dāng)儲(chǔ)層孔隙度發(fā)生變化時(shí)，孔隙度為30%的水層其振幅與孔隙度為26%的油層振幅基本相當(dāng)。即高孔隙含水砂巖也會(huì)形成強(qiáng)振幅“亮點(diǎn)”特征，與油層表現(xiàn)出的特征相似，從而導(dǎo)致利用疊后振幅預(yù)測(cè)流體性質(zhì)存在多解性。

圖3 疊后振幅影響因素分析

2.2 疊前AVO響應(yīng)特征及影響因素

圖4所示為W1、W2和W3井旁道集，可以看到，油層和水層表現(xiàn)為遠(yuǎn)道振幅增強(qiáng)的II-III類(lèi)AVO異常，異常強(qiáng)度表現(xiàn)為W3>W1>W2井，利用常規(guī)P*G屬性很難準(zhǔn)確識(shí)別流體性質(zhì)。

圖4 W1、W2和W3井旁地震道集

以W1井鉆遇油層為例，開(kāi)展流體和孔隙度替代分析，研究?jī)烧邔?duì)疊前地震響應(yīng)特征的影響。首先，保持原狀流體性質(zhì)不變，孔隙度Por由22%增加到34%，間隔3%。然后，將原狀流體替換為水層，再做一次孔隙度替換，正演道集結(jié)果如圖5所示，截距—梯度交會(huì)結(jié)果如圖6所示。

圖5 油層(a)及水層(b)AVO特征隨孔隙度變化

據(jù)圖6可知：在相同孔隙度條件下，油和水的截距、梯度均有差異，此時(shí)單獨(dú)依靠截距或梯度均可以區(qū)分油水。在不同孔隙度條件下，油和水的截距、梯度存在一定重疊。例如，孔隙度為31%的水層，其梯度與孔隙度為34%的油層相當(dāng)，其截距與孔隙度為28%的油層相當(dāng)。在截距—梯度交會(huì)圖上，高孔隙含油和含水砂巖均位于第三象限的位置，從而導(dǎo)致利用常規(guī)AVO屬性不能有效區(qū)分油水。

圖6 不同流體、孔隙度的截距—梯度交會(huì)

3 流體因子敏感性定量評(píng)價(jià)

地下介質(zhì)充填不同流體類(lèi)型時(shí)會(huì)表現(xiàn)出不同的巖石物理特征，從而為通過(guò)疊前反演實(shí)現(xiàn)儲(chǔ)層流體檢測(cè)提供了可能。其中，最為關(guān)鍵的環(huán)節(jié)在于篩選出真正反映目標(biāo)區(qū)流體性質(zhì)差異的敏感彈性參數(shù)，即流體敏感因子。常用的流體因子選擇方法是通過(guò)測(cè)井曲線計(jì)算目標(biāo)層段油氣層與水層的流體因子平均值之差來(lái)衡量流體識(shí)別的能力，并沒(méi)有考慮孔隙度對(duì)流體識(shí)別效果的影響。然而孔隙度是導(dǎo)致本區(qū)烴類(lèi)檢測(cè)多解性的關(guān)鍵因素，適用于研究區(qū)的流體因子必須對(duì)流體性質(zhì)變化敏感的同時(shí)對(duì)儲(chǔ)層孔隙度變化不敏感。針對(duì)該實(shí)際情況，本文提出了一種基于流體替代和孔隙度替代的流體因子敏感性定量評(píng)價(jià)方法，具體實(shí)現(xiàn)步驟如下：

1)保持原狀地層孔隙度不變，采用Gassmann方程進(jìn)行流體替代，根據(jù)式(1)計(jì)算流體敏感系數(shù)A，這是目前進(jìn)行流體因子篩選的常規(guī)方法，并沒(méi)有考慮孔隙度的影響：

(1)

其中，Pw代表含水狀態(tài)下的彈性參數(shù)值，Po代表含油狀態(tài)下的彈性參數(shù)值。

2)保持原狀地層流體性質(zhì)不變，進(jìn)行孔隙度替代，根據(jù)式(2)計(jì)算孔隙度敏感系數(shù)B：

(2)

其中，P?1代表孔隙度為?1時(shí)的彈性參數(shù)值，P?2代表孔隙度為?2的彈性參數(shù)值。

3)為了突出流體、壓制孔隙度的影響，定義評(píng)價(jià)因子C，其與參數(shù)A和B的關(guān)系如式(3)：

(3)

C值越大代表對(duì)流體敏感性越高、對(duì)孔隙度敏感性越低。當(dāng)A=B時(shí)，C=0；當(dāng)AB時(shí)，C>0，極限值為1。根據(jù)C值大小可優(yōu)選出研究區(qū)的最佳敏感流體因子。

以W1井鉆遇油層為例，開(kāi)展流體和孔隙度替代分析，流體性質(zhì)由原狀油層替換成水層，孔隙度由原狀27%替換到31%，獲得不同彈性參數(shù)的值，然后計(jì)算對(duì)應(yīng)的敏感性系數(shù)A、B和評(píng)價(jià)因子C，結(jié)果如表2所示。對(duì)表中數(shù)據(jù)進(jìn)行直方圖統(tǒng)計(jì)分析，從圖7a中可以看出，λρ、λ/μ、PI和f這4個(gè)參數(shù)對(duì)流體性質(zhì)均較為敏感，其中λρ敏感性最高，按照常規(guī)分析會(huì)選擇λρ作為研究區(qū)的敏感流體因子。但是根據(jù)圖7b可知，λρ、PI和f對(duì)孔隙度也具有較高的敏感性，容易受到孔隙度變化的影響，導(dǎo)致流體檢測(cè)出現(xiàn)多解性，而參數(shù)λ/μ對(duì)孔隙度敏感性低。通過(guò)圖7c的評(píng)價(jià)因子對(duì)比，可更加直觀地發(fā)現(xiàn)，壓制孔隙度影響以后，相比其他參數(shù)，λ/μ具有最高的流體敏感性。所以，λ/μ是研究區(qū)去除物性影響后的最佳流體因子。

表2 流體和孔隙度替代后彈性參數(shù)值及敏感系數(shù)計(jì)算結(jié)果

圖7 不同彈性參數(shù)的流體敏感系數(shù)A(a)、孔隙度敏感系數(shù)B(b)和評(píng)價(jià)因子C(c)

4 應(yīng)用效果

圖8所示為基于常規(guī)方法篩選的流體敏感因子λρ與利用本文提出方法篩選的流體敏感因子λ/μ對(duì)流體識(shí)別的效果對(duì)比。通過(guò)圖8a可以看出，當(dāng)油層和水層孔隙度相當(dāng)或者油層孔隙度大于水層時(shí)，λρ能夠有效區(qū)分油和水；隨著水層孔隙度增大，高孔隙度水層與相對(duì)低孔隙度油層的λρ存在重疊，流體識(shí)別效果變差。通過(guò)圖8b可以看到，λ/μ能夠有效區(qū)分不同孔隙度條件下的油層和水層，當(dāng)λ/μ<1.2時(shí)表現(xiàn)為油層，進(jìn)一步證實(shí)了λ/μ具有良好的烴類(lèi)檢測(cè)能力。

圖8 λρ(a)及λ/μ(b)流體識(shí)別效果對(duì)比

圖9所示為疊前反演λρ和λ/μ進(jìn)行烴類(lèi)檢測(cè)的效果對(duì)比。從圖9a中可以看出，λρ容易受孔隙度的影響，雖然W1和W2井油層有異常，但是W3井水層異常更為明顯，烴類(lèi)檢測(cè)結(jié)果與實(shí)鉆結(jié)果不符。從圖9b中可以看出，利用彈性參數(shù)λ/μ開(kāi)展烴類(lèi)檢測(cè)，W1和W2井油層表現(xiàn)為低λ/μ值，W3井水層表現(xiàn)為相對(duì)高λ/μ值，高孔水砂巖的異常假象得到壓制，油水得到有效區(qū)分。

圖9 過(guò)W1、W2、W3井的λρ反演剖面(a)和流體因子λ/μ反演剖面(b)

利用流體因子λ/μ開(kāi)展R1180層油氣平面展布規(guī)律預(yù)測(cè)，結(jié)果如圖10所示。從圖10a中可以觀察到研究區(qū)發(fā)育東、西兩支水道化朵葉沉積體，W1井鉆遇西支朵葉，W2和W3井鉆遇東支朵葉，朵葉砂體在平面上表現(xiàn)為振幅“亮點(diǎn)”異常，西支朵葉振幅級(jí)別要強(qiáng)于東支朵葉。從λ/μ反演結(jié)果(圖10b)可以看到，W3井水層對(duì)應(yīng)的強(qiáng)振幅異常假象得到了明顯壓制，W1和W2井油層異常得到了有效保留。西支朵葉大面積的強(qiáng)振幅區(qū)域基本都是高孔含水砂巖的響應(yīng)，只有局部位置發(fā)育油氣層，從而規(guī)避了高孔含水砂巖帶來(lái)的勘探風(fēng)險(xiǎn)，有效指導(dǎo)了西支朵葉的下一步鉆探部署。

圖10 R1180層均方根振幅(a)和流體因子λ/μ反演結(jié)果(b)

5 結(jié)論和認(rèn)識(shí)

1)復(fù)雜儲(chǔ)層的流體檢測(cè)受孔隙度影響較大，如何將孔隙度的影響從流體因子中剔除是獲取準(zhǔn)確流體識(shí)別效果的關(guān)鍵。本文針對(duì)疊前反演敏感流體參數(shù)優(yōu)選問(wèn)題，提出了基于流體和孔隙度替代的流體因子敏感性定量分析方法，篩選出的流體因子具有“對(duì)流體性質(zhì)敏感度高、對(duì)孔隙度敏感度低”的特征，能夠最大限度地突出流體異常，提高流體識(shí)別精度。該方法和思路可以推廣應(yīng)用到其他類(lèi)似區(qū)塊，通過(guò)優(yōu)選適用于區(qū)塊自身的流體敏感因子開(kāi)展烴類(lèi)檢測(cè)工作。

2)將該方法應(yīng)用于尼日爾三角洲盆地S區(qū)塊，優(yōu)選λ/μ作為最佳流體敏感因子，實(shí)際應(yīng)用效果證實(shí)λ/μ能夠有效壓制高孔含水砂巖在地震剖面上表現(xiàn)出的強(qiáng)振幅“亮點(diǎn)”假象，提高烴類(lèi)檢測(cè)成功率。