劉立峰，孫贊東，楊海軍，韓劍發(fā)，敬兵，趙海濤

(1. 中國(guó)石油大學(xué)(北京) 地質(zhì)地球物理綜合研究中心，北京，102249；2. 中國(guó)石油塔里木油田分公司 勘探開(kāi)發(fā)研究院，新疆 庫(kù)爾勒，841000)

塔里木盆地碳酸鹽巖基質(zhì)孔隙度低，滲透性能差，不能成為有效的儲(chǔ)集空間，其化學(xué)成分和物理特性決定了碳酸鹽巖儲(chǔ)層更容易遭受各種因素的影響和改造。塔里木盆地經(jīng)歷過(guò)多期次的構(gòu)造升降運(yùn)動(dòng)、多次暴露地表的風(fēng)化淋濾、盆內(nèi)各種流體的作用以及大規(guī)?；鹕交顒?dòng)，使得區(qū)內(nèi)不同區(qū)域經(jīng)受了各種條件下成巖作用的改造和疊加，形成了大量的溶孔、溶洞和溶縫，極大地改善碳酸鹽巖儲(chǔ)層的儲(chǔ)集性能。因此，有效儲(chǔ)集空間類型是經(jīng)過(guò)多期多旋回構(gòu)造、溶蝕等作用改造后形成的次生溶蝕孔洞和裂縫，及它們交互耦合而成的網(wǎng)狀裂縫等，常常表現(xiàn)為縫－縫交錯(cuò)、縫－洞相連、縫－洞疊合，高產(chǎn)、穩(wěn)產(chǎn)的關(guān)鍵是鉆遇裂縫和溶蝕孔洞發(fā)育帶[1-3]。多年油氣勘探實(shí)踐證明，奧陶系碳酸鹽巖是塔里木盆地尋找大型油氣田的重要領(lǐng)域，現(xiàn)已相繼在塔河、輪南、塔中等地區(qū)都發(fā)現(xiàn)了億噸級(jí)的大油氣田，具有巨大的勘探潛力[4-5]。但是，縫洞型碳酸鹽巖儲(chǔ)層非均質(zhì)性極強(qiáng)，其外部形態(tài)和內(nèi)部結(jié)構(gòu)極不規(guī)則，縱向上和橫向上變化快，縫洞儲(chǔ)集體與致密碳酸鹽巖基質(zhì)隨機(jī)交互出現(xiàn)，儲(chǔ)集能力和滲流能力空間分布不均勻，地球物理特性十分復(fù)雜，儲(chǔ)層分布異常復(fù)雜。同時(shí)，奧陶系碳酸鹽巖地層埋藏深，地震反射信號(hào)弱，地震資料主頻和信噪比都較低，地震資料品質(zhì)不理想，運(yùn)用單因素的地質(zhì)或地球物理方法難以準(zhǔn)確預(yù)測(cè)碳酸鹽巖儲(chǔ)層的分布[6-7]。因此，對(duì)該地區(qū)奧陶系碳酸鹽巖溶洞、溶孔、裂縫發(fā)育區(qū)以及碳酸鹽巖有效儲(chǔ)層的預(yù)測(cè)，依然是當(dāng)前塔里木盆地碳酸鹽巖油氣勘探實(shí)現(xiàn)擴(kuò)大發(fā)現(xiàn)、降低鉆探風(fēng)險(xiǎn)所面臨的主要難題。

1 縫洞型碳酸鹽巖儲(chǔ)層預(yù)測(cè)技術(shù)應(yīng)用效果分析

本文選取塔里木盆地中古21井區(qū)為研究區(qū)。該工區(qū)位于塔中Ⅰ號(hào)斷裂帶的西北部，存在多個(gè)勘探目的層段，包括石炭系的東河砂巖和志留系的瀝青砂巖等，自從中古21井和中古11井在下奧陶統(tǒng)獲得高產(chǎn)工業(yè)油氣流后，從而取得了該區(qū)下奧陶統(tǒng)勘探的重大突破，展示了良好的勘探前景。但由于奧陶系碳酸鹽巖儲(chǔ)層十分復(fù)雜，且目前只有3口井鉆至下奧陶統(tǒng)，其碳酸鹽巖儲(chǔ)層的分布規(guī)律目前仍不十分清楚，制約了該區(qū)油氣勘探的擴(kuò)大發(fā)現(xiàn)。針對(duì)縫洞型碳酸鹽巖儲(chǔ)層的特點(diǎn)，充分利用多種地震信息進(jìn)行綜合預(yù)測(cè)，篩選出一套具有針對(duì)性的儲(chǔ)層預(yù)測(cè)技術(shù)系列，以期對(duì)塔里木盆地縫洞型油氣藏勘探起到有益的指導(dǎo)作用。分析結(jié)果表明，碳酸鹽巖儲(chǔ)層在地震剖面上表現(xiàn)為強(qiáng)振幅、低相干、低頻率、高衰減、低阻抗的地震響應(yīng)特征。

1.1 振幅分析技術(shù)

振幅是最基本也是最重要的地震屬性，如果不考慮地震數(shù)據(jù)采集、處理的影響，那么，在碳酸鹽巖儲(chǔ)

層中影響振幅的則是其孔、洞、縫的發(fā)育程度，能夠揭示儲(chǔ)層橫向展布規(guī)律及儲(chǔ)層內(nèi)部流體的變化[7]。碳酸鹽巖縫洞體的反射特征一般表現(xiàn)為：縫洞體頂部相對(duì)弱振幅，縫洞體底界振幅強(qiáng)，縫洞體內(nèi)部振幅較圍巖振幅強(qiáng)、相對(duì)于底界振幅弱[8]。與實(shí)際的鉆井、測(cè)井資料對(duì)比分析發(fā)現(xiàn)：當(dāng)下奧陶統(tǒng)碳酸鹽巖內(nèi)部存在縫洞時(shí)，局部孔隙度增大、速度降低，縫洞體與圍巖速度差異增大，可以在縫洞內(nèi)部或外部形成強(qiáng)振幅反射，即在弱振幅背景下強(qiáng)振幅(即弱中強(qiáng))區(qū)域縫洞儲(chǔ)層較為發(fā)育(圖1)。

圖1 中古21井區(qū)下奧陶統(tǒng)均方根振幅數(shù)據(jù)體Fig.1 Lower Ordovician RMS Amplitude data of ZG21 well area

1.2 相干體技術(shù)

相干體技術(shù)的核心就是對(duì)地震數(shù)據(jù)進(jìn)行求同存異，以突出那些不相干的數(shù)據(jù)，通過(guò)計(jì)算橫向上波形的相似性，得到三維地震相關(guān)性的估計(jì)值，用以識(shí)別細(xì)微的巖層橫向非均一性和斷裂特征[9]。雖然不能直接預(yù)測(cè)裂縫或孔洞， 但能從宏觀上預(yù)測(cè)碳酸鹽巖孔、洞、縫的發(fā)育帶。研究區(qū)目的層地層產(chǎn)狀較緩，巖性變化不大，在排除因斷裂錯(cuò)位而在橫向上出現(xiàn)的非相似性情況下，影響地震道不相關(guān)因素主要是溶蝕孔洞及裂縫發(fā)育帶和微小斷裂，即碳酸鹽巖的主要儲(chǔ)集空間。對(duì)中古 21井區(qū)進(jìn)行了基于本征算法相干體的提取，以檢測(cè)縫洞型碳酸鹽巖儲(chǔ)層裂縫發(fā)育帶和溶洞發(fā)育特征(圖2)，發(fā)現(xiàn)高產(chǎn)油氣井(如ZG21和ZG11)都處于弱相干值區(qū)域，而處于強(qiáng)相干區(qū)域井的儲(chǔ)層發(fā)育較差(如 TZ35)。

1.3 頻譜分解技術(shù)

圖2 中古21井區(qū)下奧陶統(tǒng)相干數(shù)據(jù)體Fig.2 Lower Ordovician coherence data of ZG21 well area

頻譜分解技術(shù)是一種基于頻譜分析的地震成像方法，在頻率域內(nèi)通過(guò)分析地震信號(hào)的時(shí)變特征，可以揭示儲(chǔ)層橫向之間由于巖性、物性等因素引起的微小振幅變化，刻畫(huà)碳酸鹽巖儲(chǔ)層中縫洞發(fā)育所引起的地震反射波振幅和頻率的變化異常，達(dá)到精細(xì)預(yù)測(cè)描述儲(chǔ)層的目的[10]。其中非正交Gabor-Morlet小波變換的分頻方法直接使用頻率尺度參數(shù)來(lái)控制頻譜分解，避免了常規(guī)小波變換尺度參數(shù)與頻率不直接對(duì)應(yīng)的缺陷，同時(shí)克服了常規(guī)譜分解方法時(shí)窗對(duì)計(jì)算結(jié)果的影響，增強(qiáng)了地震信號(hào)的局部特性，提高了穩(wěn)定性和分辨率[11]。本次研究采用基于非正交Gabor-Morlet小波變換的分頻方法對(duì)縫洞型碳酸鹽巖儲(chǔ)層進(jìn)行預(yù)測(cè)，用測(cè)井解釋和產(chǎn)油氣情況等單井資料，對(duì)分頻屬性進(jìn)行標(biāo)定，發(fā)現(xiàn) 22Hz分頻屬性異常值與碳酸鹽巖縫洞型儲(chǔ)層段有良好的對(duì)應(yīng)關(guān)系(圖3)，高產(chǎn)油氣井(如ZG21和 ZG11)儲(chǔ)層發(fā)育段都位于分頻屬性調(diào)諧振幅高值異常區(qū)，而干井(如TZ35)位于分頻屬性調(diào)諧振幅低值區(qū)。

圖3 中古21井區(qū)下奧陶統(tǒng)分頻數(shù)據(jù)體(22 Hz)Fig.3 Lower Ordovician frequency division data (22 Hz)of ZG21 well area

1.4 地震波吸收衰減技術(shù)

大量的正演模型和理論研究表明[12-13]：地震波在聚集了石油、天然氣的儲(chǔ)層中傳播時(shí)，高頻能量衰減要比低頻能量衰減大，即存在“高頻吸收”現(xiàn)象。在吸收衰減屬性中，地震波頻率衰減梯度對(duì)檢測(cè)流體和預(yù)測(cè)儲(chǔ)層分布最為敏感。頻率衰減梯度是指在頻譜分解基礎(chǔ)上的高頻端振幅包絡(luò)的擬合斜率，表示高頻段的地震波能量隨頻率的變化情況。利用地震波頻率衰減梯度來(lái)檢測(cè)縫洞系統(tǒng)對(duì)地震波高頻成分的吸收衰減情況(圖4)。通過(guò)對(duì)比分析發(fā)現(xiàn)：各井碳酸鹽巖縫洞儲(chǔ)層發(fā)育段與高頻率衰減梯度異常有良好的對(duì)應(yīng)關(guān)系，位于高頻率衰減梯度分布區(qū)的鉆井(如ZG21和ZG11)油氣產(chǎn)量很高，而位于低頻率衰減梯度分布區(qū)則儲(chǔ)層較差(如TZ35)，表明利用地震波吸收衰減屬性能有效地對(duì)縫洞型碳酸鹽巖儲(chǔ)層進(jìn)行預(yù)測(cè)。

圖4 中古21井區(qū)下奧陶統(tǒng)頻率衰減梯度數(shù)據(jù)體Fig.4 Lower Ordovician frequency attenuation gradient data of ZG21 well area

1.5 多井約束波阻抗反演技術(shù)

多井約束波阻抗反演通過(guò)測(cè)井與地震資料的融合，將測(cè)井資料豐富的高頻信息和完整的低頻成分加入到反演結(jié)果中，有效地補(bǔ)充了地震信號(hào)有效帶寬的不足，能得到高分辨率的地層波阻抗資料。碳酸鹽巖通常表現(xiàn)為高速度、高密度、高波阻抗，當(dāng)存在溶蝕孔洞和裂縫時(shí)，特別是被油氣水或與圍巖差異大的巖性充填后，會(huì)引起速度、密度、波阻抗降低；因此，可以利用波阻抗信息來(lái)預(yù)測(cè)碳酸鹽巖縫洞儲(chǔ)層的分布(圖 5)。高產(chǎn)工業(yè)油氣流井(如 ZG21和 ZG11)位于低阻抗區(qū)域，而干井(如TZ35)位于相對(duì)高的阻抗區(qū)域，即縫洞儲(chǔ)層發(fā)育的部位對(duì)應(yīng)的是高阻抗背景之下的弱阻抗區(qū)。

圖5 中古21井區(qū)下奧陶統(tǒng)波阻抗數(shù)據(jù)體Fig.5 Lower Ordovician impedance data of ZG21 well area

2 地震屬性優(yōu)化及綜合預(yù)測(cè)

每一種儲(chǔ)層預(yù)測(cè)技術(shù)都是從不同角度反映地震信息在縱向和橫向上的變化，它們與巖性、物性、孔隙流體性質(zhì)之間的關(guān)系非常復(fù)雜，是多種復(fù)雜地質(zhì)因素或參數(shù)的綜合反映。每種預(yù)測(cè)技術(shù)都有其自身的特殊性與局限性，存在預(yù)測(cè)的多解性和精度問(wèn)題。本文將相干體、波阻抗等“特殊屬性”都劃入“廣義地震屬性”的范疇，以便和常規(guī)屬性一起進(jìn)行優(yōu)化處理，對(duì)碳酸鹽巖儲(chǔ)層進(jìn)行綜合預(yù)測(cè)。

地震屬性優(yōu)化技術(shù)可以從眾多地震屬性中挑選出與研究目標(biāo)關(guān)系最密切、反應(yīng)最敏感、相互獨(dú)立的優(yōu)勢(shì)屬性，利用優(yōu)化算法將多種地震屬性進(jìn)行優(yōu)化處理，能明顯降低多解性，提高預(yù)測(cè)精度[14]。目前，地震屬性優(yōu)化算法很多，如聚類分析、因子分析、搜索算法、決策分析、遺傳與神經(jīng)網(wǎng)絡(luò)聯(lián)合分析等，但實(shí)際的應(yīng)用效果并不是很理想，都存在著優(yōu)化結(jié)果穩(wěn)定性和準(zhǔn)確性較差等缺點(diǎn)，特別是針對(duì)塔里木盆地這樣復(fù)雜的縫洞型碳酸鹽巖儲(chǔ)層，盲目的應(yīng)用還可能導(dǎo)致產(chǎn)生錯(cuò)誤的預(yù)測(cè)結(jié)果。

本文利用多元逐步判別分析和核主成分分析等相結(jié)合的方法對(duì)多種地震屬性進(jìn)行優(yōu)化處理。該方法可以綜合多種地震屬性的有效信息，建立地震屬性與儲(chǔ)層信息的某種線性或非線性函數(shù)關(guān)系，能實(shí)現(xiàn)地震屬性到地質(zhì)參數(shù)的轉(zhuǎn)換，同時(shí)可以將其中判別能力差或起干擾因素的地震屬性剔除，最終在判別函數(shù)中只保留數(shù)量不多而判別能力強(qiáng)的屬性參數(shù)，并且能有效地避免由于多地震屬性之間具有的相關(guān)性造成的信息重復(fù)和冗余，從而對(duì)井間未知區(qū)域的儲(chǔ)層信息進(jìn)行準(zhǔn)確的綜合判識(shí)。

具體步驟如下(見(jiàn)圖6)。

圖6 地震多屬性分析及優(yōu)化儲(chǔ)層預(yù)測(cè)工作流程Fig.6 Processing chart for seismic multi-attribute analysis and optimization of reservoir prediction

(1) 地震屬性的優(yōu)選。不同的地震屬性所代表的地質(zhì)意義不完全相同，一些屬性可能對(duì)預(yù)測(cè)結(jié)果起著干擾作用，為此必須在眾多的地震屬性中優(yōu)選那些有用的信息[15]。地震屬性優(yōu)選是指從地震屬性集合中選出合理的地震屬性子集，主要從地震屬性的物理意義、敏感性和相關(guān)性等方面考慮。在上述分析的基礎(chǔ)上，對(duì)研究區(qū)奧陶系碳酸鹽巖儲(chǔ)層提取多種地震屬性體，然后結(jié)合研究區(qū)實(shí)際地質(zhì)情況，利用測(cè)井解釋和產(chǎn)油氣情況等單井資料進(jìn)行標(biāo)定，選擇多個(gè)對(duì)碳酸鹽巖儲(chǔ)層反映敏感且具有明確物理意義的屬性，并進(jìn)行屬性間的相關(guān)分析，選擇彼此間相關(guān)性小的屬性集，最終確定8種優(yōu)勢(shì)地震屬性體，包括均方根振幅數(shù)據(jù)體、瞬時(shí)頻率數(shù)據(jù)體、相干數(shù)據(jù)體、弧長(zhǎng)數(shù)據(jù)體、能量半衰時(shí)數(shù)據(jù)體、分頻調(diào)諧能量數(shù)據(jù)體、頻率衰減梯度數(shù)據(jù)體、波阻抗數(shù)據(jù)體。

(2) 地震屬性預(yù)處理。地震屬性的提取難免會(huì)受到個(gè)別地段信噪比較低的影響和噪聲的干擾，使得參數(shù)出現(xiàn)“毛刺”、“野值”等異常數(shù)據(jù)，這樣就會(huì)干擾地震屬性參數(shù)的地質(zhì)標(biāo)定，且利用屬性參數(shù)進(jìn)行模式識(shí)別時(shí)，容易造成解釋中出現(xiàn)“假異?！盵14]。因此，對(duì)優(yōu)選的8種地震屬性數(shù)據(jù)體內(nèi)的異常數(shù)據(jù)進(jìn)行手工剔除和適當(dāng)?shù)钠交瑸V波處理，并且進(jìn)行中心標(biāo)準(zhǔn)化(零均值化)處理，使各種屬性的數(shù)值變換到某種相同的尺度之下，處理后每種地震屬性的均值為零，標(biāo)準(zhǔn)差為1。

(3) 選取訓(xùn)練樣本。根據(jù)研究區(qū)鉆井所揭示的儲(chǔ)層實(shí)際情況，將奧陶系碳酸鹽巖儲(chǔ)層分為油氣層、油氣顯示層和非儲(chǔ)層3類，將其作為建立判別模型的分類。綜合測(cè)井解釋結(jié)果及試油、巖芯等資料確定各井在碳酸鹽巖層段內(nèi)各深度點(diǎn)的儲(chǔ)層類型。由于地震數(shù)據(jù)體是時(shí)間域的，而由已知鉆井得到的儲(chǔ)層信息是深度域的，因此，需要利用高精度的合成地震記錄將其由深度域轉(zhuǎn)換成時(shí)間域，才能確定每口井在時(shí)間域內(nèi)每個(gè)采樣點(diǎn)上的儲(chǔ)層類型。并假設(shè)距井口100 m范圍內(nèi)儲(chǔ)層類型與井口位置一致，選取奧陶系碳酸鹽巖層段井旁100 m范圍內(nèi)的屬性點(diǎn)作為優(yōu)化的訓(xùn)練樣本。

(4) 地震屬性多元逐步判別。各種地震屬性對(duì)儲(chǔ)層參數(shù)的影響是不同的，有的起主要作用，有的起次要作用，有的甚至起干擾作用。多元逐步判別可以將其中判別能力最差的地震屬性剔除，最終在判別函數(shù)中只保留數(shù)量不多而判別能力強(qiáng)的屬性參數(shù)。利用多元逐步判別計(jì)算經(jīng)過(guò)標(biāo)準(zhǔn)化后的8種優(yōu)勢(shì)地震屬性體對(duì)判別模型的貢獻(xiàn)，其中頻率衰減梯度數(shù)據(jù)體對(duì)模型的貢獻(xiàn)最大，其次是波阻抗數(shù)據(jù)體、分頻調(diào)諧能量數(shù)據(jù)體、均方根振幅數(shù)據(jù)體、相干數(shù)據(jù)體、能量半衰時(shí)數(shù)據(jù)體，用上述6種地震屬性數(shù)據(jù)體建立的判別函數(shù)具有非常顯著的判別效果，而瞬時(shí)頻率數(shù)據(jù)體、弧長(zhǎng)數(shù)據(jù)體對(duì)模型的判別效果不顯著或起干擾作用，經(jīng)逐步判別將其剔除。

(5) 核主成分分析。在多數(shù)情況下，優(yōu)選后的地震屬性集的空間維數(shù)仍然較高，必然會(huì)增加預(yù)測(cè)分類的計(jì)算量；同時(shí)，優(yōu)選的地震屬性即使考慮了相關(guān)屬性的篩選，仍然不能保證屬性間絕對(duì)的不相關(guān)，存在著冗余信息，直接影響到預(yù)測(cè)結(jié)果的穩(wěn)定性，因此，必須對(duì)地震屬性空間進(jìn)行壓縮和相關(guān)性處理，以提高模式識(shí)別的速度與精度。核主成分分析可將原來(lái)的地震屬性參數(shù)轉(zhuǎn)化為彼此獨(dú)立的一組新綜合變量，同時(shí)能有效地處理多個(gè)地震屬性之間的非線性關(guān)系。

對(duì)經(jīng)過(guò)多元逐步判別所引入的6種地震屬性數(shù)據(jù)體進(jìn)行核主成分分析，將其轉(zhuǎn)變?yōu)楸舜霜?dú)立的新變量Y1，Y2，…，Y6，計(jì)算特征空間 F中協(xié)方差矩陣的特征值和累計(jì)貢獻(xiàn)率。其中前3個(gè)核主成分的累計(jì)貢獻(xiàn)率達(dá)到95%，集中了原始地震屬性的絕大部分信息，并能很好地概括各儲(chǔ)層類型之間的差異。

(6) 綜合判別分析及儲(chǔ)層綜合預(yù)測(cè)。最后根據(jù)選取的訓(xùn)練樣本，對(duì)核主成分分析得到的前3個(gè)核主成分進(jìn)行基于貝葉斯準(zhǔn)則的判別分析，建立該地區(qū)儲(chǔ)層的判別模型[16]。將研究區(qū)內(nèi)各屬性點(diǎn)代入建立的判別模型，就可以得到該地區(qū)奧陶系碳酸鹽巖儲(chǔ)層類型在三維空間中的分布情況，并對(duì)各儲(chǔ)層類型的概率進(jìn)行歸一化處理，得到多屬性優(yōu)化數(shù)據(jù)體(見(jiàn)圖7)，含油氣區(qū)用2～3表示，其值越大，表示油氣儲(chǔ)層類型的概率越大；含油氣顯示區(qū)用1～2表示，非儲(chǔ)層區(qū)用0～1表示，進(jìn)而可以對(duì)該地區(qū)奧陶系碳酸鹽巖儲(chǔ)層進(jìn)行定量預(yù)測(cè)。

圖7 中古21井區(qū)下奧陶統(tǒng)多屬性優(yōu)化數(shù)據(jù)體Fig.7 Lower Ordovician multi-attribute optimization data of ZG21 well area

建立在多種統(tǒng)計(jì)方法分析基礎(chǔ)上的地震多屬性優(yōu)化方法，綜合考慮了多種地震屬性的有效信息，以已知井?dāng)?shù)據(jù)作為訓(xùn)練樣本，并經(jīng)過(guò)多元逐步判別和核主成分分析后，對(duì)判別結(jié)果不顯著甚至是起干擾因素的信息已被剔除，并利用相對(duì)獨(dú)立的核主成分進(jìn)行綜合判別分析，保證了運(yùn)算結(jié)果的穩(wěn)定性，同時(shí)考慮了各地震屬性間的非線性關(guān)系，提高了預(yù)測(cè)結(jié)果的可信度。圖8所示為研究區(qū)各井的優(yōu)化數(shù)據(jù)體過(guò)井剖面，圖9所示為下奧陶統(tǒng)鷹山組頂界多屬性優(yōu)化數(shù)據(jù)體沿層切片。通過(guò)對(duì)比分析發(fā)現(xiàn)，多屬性優(yōu)化結(jié)果預(yù)測(cè)的儲(chǔ)層與各井碳酸鹽巖縫洞型儲(chǔ)層發(fā)育層段都有良好的對(duì)應(yīng)關(guān)系。圖8中箭頭指示各井的產(chǎn)油氣層段，說(shuō)明優(yōu)化結(jié)果能對(duì)碳酸鹽巖型儲(chǔ)層做到準(zhǔn)確定位，可以很好地反映碳酸鹽巖儲(chǔ)層的位置和空間形態(tài)。

圖8 中古21井區(qū)下奧陶統(tǒng)多屬性優(yōu)化數(shù)據(jù)體過(guò)井剖面Fig.8 Well profiles of Lower Ordovician multi-attribute optimization data of ZG21 well area

圖9 下奧陶統(tǒng)鷹山組頂界優(yōu)化數(shù)據(jù)體沿層切片F(xiàn)ig.9 Hierarchy slicing of multi-attribute optimization data of top surface of YingShan formation, Lower Ordovician

圖10 中古21井區(qū)下奧陶統(tǒng)溶蝕孔洞儲(chǔ)層三維雕刻Fig.10 3D carving Lower Ordovician solution pole reservoir of ZG21 well area

運(yùn)用三維可視化技術(shù)對(duì)多屬性優(yōu)化數(shù)據(jù)體進(jìn)行精細(xì)刻畫(huà)(圖10)，可以準(zhǔn)確預(yù)測(cè)溶蝕孔洞型碳酸鹽巖儲(chǔ)集體在三維空間的有利分布范圍，深化認(rèn)識(shí)碳酸鹽巖儲(chǔ)層的空間分布規(guī)律，進(jìn)而使該地區(qū)碳酸鹽巖儲(chǔ)層建模、油藏動(dòng)態(tài)模擬和制定開(kāi)發(fā)方案變?yōu)榭赡埽蔀樘岣咩@探成功率，井位的部署，指導(dǎo)勘探開(kāi)發(fā)提供重要依據(jù)。

3 結(jié)論

(1) 針對(duì)縫洞型碳酸鹽巖儲(chǔ)層的特點(diǎn)，通過(guò)中古21井區(qū)地震、地質(zhì)、鉆井等資料的綜合分析，形成了以振幅分析技術(shù)、基于本征算法的相干體技術(shù)、頻譜分解技術(shù)、地震波吸收衰減技術(shù)和測(cè)井約束波阻抗反演技術(shù)為核心的技術(shù)系列。

(2) 利用多元逐步判別分析和核主成分分析相結(jié)合的辦法對(duì)多種地震屬性進(jìn)行優(yōu)化求解。該方法綜合考慮了多種地震屬性的有效信息，經(jīng)過(guò)多元逐步判別和核主成分分析后，剔除了判別能力差或起干擾因素的地震屬性，最終在判別函數(shù)中只保留數(shù)量不多而判別能力強(qiáng)的屬性參數(shù)，并利用相對(duì)獨(dú)立的核主成分進(jìn)行綜合判別分析，保證了運(yùn)算結(jié)果的穩(wěn)定性；同時(shí)考慮了各地震屬性間的非線性關(guān)系，提高了預(yù)測(cè)結(jié)果的可信度，有效解決了多種地震屬性存在的多解性問(wèn)題。其優(yōu)化結(jié)果能較好地刻畫(huà)縫洞型碳酸鹽巖儲(chǔ)層的有利分布范圍，取得了良好的地質(zhì)效果，以期為地震多屬性優(yōu)化與縫洞型碳酸鹽巖儲(chǔ)層預(yù)測(cè)研究提供一種科學(xué)、有效的方法。

[1] 王世星, 曹輝蘭, 靳文芳, 等. 碳酸鹽巖縫洞系統(tǒng)地震響應(yīng)特征分析和塔中卡1區(qū)縫洞儲(chǔ)層預(yù)測(cè)[J]. 石油物探, 2005, 44(5):421-427.WANG Shi-xing, CAO Hui-lan, JIN Wen-fang, et al. Seismic response and prediction of fracture-cavity system in carbonate reservoir: A case study in the Ka-1 field[J]. Geophysical Prospecting for Petrole, 2005, 44(5): 421-427.

[2] 韓革華, 漆立新, 李宗杰, 等. 塔河油田奧陶系碳酸鹽巖縫洞型儲(chǔ)層預(yù)測(cè)技術(shù)[J]. 石油與天然氣地質(zhì), 2006, 27(6): 860-870.HAN Ge-hua, QI Li-xin, LI Zong-jie, et al. Prediction of the Ordovician fractured-vuggy carbonate reservoirs in Tahe Oilfield[J]. Oil & Gas Geology, 2006, 27(6): 860-870.

[3] 楊子川. 塔河油田碳酸鹽巖儲(chǔ)層預(yù)測(cè)技術(shù)與應(yīng)用[J]. 勘探地球物理進(jìn)展, 2004, 27(6): 432-439.YANG Zi-chuan. Carbonate reservoir prediction technology and its application in Tahe Oilfield[J]. Progress in Exploration Geophysics, 2004, 27(6): 432-439.

[4] 趙宗舉. 海相碳酸鹽巖儲(chǔ)集層類型、成藏模式及勘探思路[J].石油勘探與開(kāi)發(fā), 2008, 35(6): 692-703.ZHAO Zong-ju. Types, accumulation models and exploration concepts of marine carbonate reservoirs[J]. Petroleum Exploration and Development, 2008, 35(6): 692-703.

[5] 陳志海, 戴勇, 郎兆新. 縫洞性碳酸鹽巖油藏儲(chǔ)滲模式及其開(kāi)采特征[J]. 石油勘探與開(kāi)發(fā), 2005, 32(3): 101-105.CHEN Zhi-hai, DAI Yong, LANG Zhao-xin. Storage-percolation modes and production performance of the karst reservoirs in Tahe Oilfield[J]. Petroleum Exploration and Development, 2005,32(3): 101-105.

[6] 楊海軍, 鄔光輝, 韓建發(fā), 等. 塔里木盆地中央隆起帶奧陶系碳酸鹽巖臺(tái)緣帶油氣富集特征[J]. 石油學(xué)報(bào), 2007, 28(4):26-28.YANG Hai-jun, WU Guang-hui, HAN Jian-fa, et al.Characteristics of hydrocarbon enrichment along the Ordovician carbonate platform margin in the central uplift of Tarim Basin[J].Acta Petrolei Sinica, 2007, 28(4): 26-28.

[7] 溫志新, 王紅漫, 漆立新, 等. 塔河油田奧陶系縫洞型碳酸鹽巖儲(chǔ)層預(yù)測(cè)研究[J]. 地學(xué)前緣, 2008, 15(1): 94-100.WEN Zhi-xin, WANG Hong-man, QI Li-xin, et al. Research on the prediction of Ordovician carbonatite reservoir of fissure-cave type in the Tahe Oilfield[J]. Earth Science Frontiers, 2008, 15(1):94-100.

[8] 蔡瑞. 基于譜分解技術(shù)的縫洞型碳酸鹽巖溶洞識(shí)別方法[J].石油勘探與開(kāi)發(fā), 2005, 32(2): 82-85.CAI Rui. Carbonate cave identification by using a spectral decomposition technique[J]. Petroleum Exploration and Development, 2005, 32(2): 82-85.

[9] 孫夕平, 杜世通. 相干體技術(shù)算法研究及其在地震資料解釋中的應(yīng)用[J]. 石油大學(xué)學(xué)報(bào): 自然科學(xué)版, 2003, 27(2): 32-35.SUN Xi-ping, DU Shi-tong. Development and application of algorithm of coherency cub technique to seismic interpretation[J].Journal of the University of Petroleum: Science and Technology,2003, 27(2): 32-35.

[10] 龔洪林, 王振卿, 李錄明, 等. 應(yīng)用地震分頻技術(shù)預(yù)測(cè)碳酸鹽巖儲(chǔ)層[J]. 地球物理學(xué)進(jìn)展, 2008, 23(1): 129-135.GONG Hong-lin, WANG Zhen-qing, LI Lu-ming, et al.Predicting carbonate reservoir by applying seismic spectral decomposition technique[J]. Progress in Geophysics, 2008, 23(1):129-135.

[11] 苑書(shū)金. 非正交小波變換的頻譜分解技術(shù)在大牛地氣田儲(chǔ)層預(yù)測(cè)中的應(yīng)用[J]. 天然氣工業(yè), 2007, 27(12): 52-54.YUAN Shu-jin. Application of spectrum decomposition of nonorthogonal wavelet transforms to reservoir prediction in Daniudi gas field[J]. Natural Gas Industry, 2007, 27(12): 52-54.[12] Sams M S, Neep J P, Worthington M H, et al. The measurement of velocity dispersion and frequency-dependent intrinsic attenuation in sedimentary rocks[J]. Geophysics, 1997, 62(5):1456-1464.

[13] Sinha S, Routh P S. Spectral decomposition of seismic data with continuous-wavelet transform[J]. Geophysics, 2005, 70(6):19-25.

[14] 印興耀, 周靜毅. 地震屬性優(yōu)化方法綜述[J]. 石油地球物理勘探, 2005, 40(4): 482-489.YIN Xing-yao, ZHOU Jing-yi. Summary of optimum methods of seismic attributes[J]. Oil Geophysical Prospecting, 2005,40(4): 482-489.

[15] 馬中高, 管路平, 賀振華, 等. 利用模型正演優(yōu)選地震屬性進(jìn)行儲(chǔ)層預(yù)測(cè)[J]. 石油學(xué)報(bào), 2003, 24(6): 35-39.MA Zhong-gao, GUAN Lu-ping, HE Zheng-hua, et al. Reservoir characterization using seismic attributes optimized with rock physics modeling[J]. Acta Petrolei Sinica, 2003, 24(6): 35-39.

[16] SUN Zan-dong. Seismic methods for heavy oil reservoir monitoring and characterization[D]. Calgary: University of Calgary. Department of Geology and Geophysics, 1999:108-117.