趙 虎,孫 勇,趙容容,羅 鑫,安虹伊,陳思锜,徐 姁
(1.西南石油大學天然氣地質(zhì)四川省重點實驗室,四川成都610500;2.西南石油大學地球科學與技術(shù)學院,四川成都610500;3.中國石油西南油氣田分公司勘探事業(yè)部,四川成都610041;4.中石油西南油氣田分公司重慶氣礦,重慶400707)
自1976年在建南構(gòu)造發(fā)現(xiàn)第一個上二疊統(tǒng)長興組生物礁氣藏,四川盆地逐漸發(fā)現(xiàn)了大量的生物礁氣藏[1-4]。特別是圍繞開江-梁平海槽不斷深化研究,陸續(xù)在海槽兩側(cè)發(fā)現(xiàn)了普光、龍崗、元壩、鐵山等多個中大型生物礁氣藏[5-8]。四川盆地長興組生物礁氣藏勘探主要集中在臺緣相帶,而對于分布范圍廣、面積小、非均質(zhì)性強的臺內(nèi)礁研究相對較少[9-10]。針對生物礁儲層的地震地質(zhì)特征,前人研究認為,古地貌是控制生物礁儲層發(fā)育規(guī)模的最關(guān)鍵因素。目前古地貌恢復(fù)技術(shù)應(yīng)用較為成熟,如曾韜等[11]選用殘厚法并結(jié)合現(xiàn)代巖溶學和巖溶動力學理論,恢復(fù)了元壩地區(qū)茅口組頂面巖溶古地貌,為該區(qū)下一步儲層勘探提供了方向;胡偉光等[12]將古地貌恢復(fù)技術(shù)與地震相分析技術(shù)、波形分類技術(shù)、分頻處理技術(shù)和相干體識別技術(shù)相結(jié)合進行生物礁儲層預(yù)測,取得了較好的效果。吳勇等[13]結(jié)合古地貌分析技術(shù),提出了一種利用古坡度和相對坡度平面異常特征對生物礁邊界進行識別的方法,并對川東北羅頂寨地區(qū)臺緣型生物礁邊界進行了預(yù)測。研究表明,識別的生物礁邊界與鉆井結(jié)果吻合率高,生物礁邊界局部特征更加精細。而針對規(guī)模小、分布散的生物礁儲層預(yù)測技術(shù),前人做了大量研究,如林昌榮等[14]提出了一種新型的局部指數(shù)擬合異常地震參數(shù)提取方法,準確地預(yù)測了地層含油氣位置。疊前預(yù)測方法在生物礁氣藏勘探中也有較為廣泛的應(yīng)用,如周路等[15]在大貓坪地區(qū)長興組生物礁儲層預(yù)測中,針對生物礁體巖性組合復(fù)雜,儲集層空間分布復(fù)雜、厚度小、非均質(zhì)性強、埋藏深度大的特點,采用疊前參數(shù)反演、波形聚類和頻譜分解技術(shù),刻畫生物礁邊界,預(yù)測了生物礁分布,得到了更加精細的生物礁內(nèi)部特征[14-16]。關(guān)于生物礁內(nèi)部刻畫方面,蒲勇[17]采用正演模擬方法建立了生物礁的地震反射特征,結(jié)果顯示,其整體特征呈現(xiàn)為丘狀或杏仁狀隆起,內(nèi)部多雜亂反射,能量時強時弱,同相軸上超現(xiàn)象明顯,為礁灘體的刻畫提供了重要的依據(jù);姚軍等[18]以長興組地層展布-礁灘體展布-礁灘體儲層分布逐級遞進思路準確預(yù)測了遂寧高能帶礁灘儲層分布;王浩等[19]從點、線、面、體多個維度剖析生物礁形態(tài),刻畫了生物礁空間展布。顯然,前人都認為對于非均質(zhì)性強的生物礁儲層的內(nèi)幕刻畫,需利用單井相識別、模型正演、古地貌分析、地震屬性優(yōu)選和儲層定量預(yù)測等技術(shù)進行綜合研究[20]。五百梯地區(qū)經(jīng)過多年勘探,認為長興組臺內(nèi)同樣存在礁灘儲層發(fā)育的有利區(qū),但臺內(nèi)礁勘探一直未取得突破,難點在于臺內(nèi)礁具有分布面積小、演化規(guī)律不明顯、非均質(zhì)性強等特點,儲層預(yù)測存在多解性。因此,如何提高儲層預(yù)測精度是重中之重[21-22]。
鑒于此,本文提出從地層、沉積相、古地貌及地震反演等方面入手,首先,在地震資料處理時拓寬頻帶,突出礁灘體的地震異常特征,同時結(jié)合長興組不同時期古地貌演化特征,明確生物礁的橫向展布特征;其次,利用地質(zhì)分析等技術(shù),明確儲層縱向發(fā)育位置,圈定礁灘體優(yōu)質(zhì)儲層發(fā)育區(qū)域。最后,結(jié)合地震屬性和儲層反演等方法,對礁灘儲層進行綜合評價,形成臺內(nèi)礁的預(yù)測模式和方法,降低礁灘儲層預(yù)測的多解性。
研究區(qū)位于四川盆地東部開江縣、開縣境內(nèi),開江—梁平海槽東側(cè),區(qū)域構(gòu)造位于印支期開江古隆起的東北斜坡上,大天池構(gòu)造帶北傾末端斷下盤的潛伏構(gòu)造,構(gòu)造形態(tài)為一短軸狀背斜,剖面形態(tài)為箱狀。工區(qū)面積為800余km2,工區(qū)內(nèi)在長興組生物礁儲層獲氣井18口,工業(yè)氣井16口(圖1)。長興組自下而上分為3段,其中,長一段厚度為40~180m,巖性以灰?guī)r夾泥頁巖為主;長二段厚度為20~145m,巖性以顆?;?guī)r、泥晶灰?guī)r、礁灰(云)巖為主;長三段厚度為10~110m,巖性以泥晶灰?guī)r、顆?;?guī)r為主,含少量礁灰(云巖)。地層厚度變化大,斜坡相帶長興組整體厚度為91~200m;臺緣帶整體厚度較大,為167~245m;臺地帶厚度為150~160m,局部存在加厚區(qū),厚度為305m(圖2)。
圖1 研究區(qū)位置
圖2 長興組斜坡-臺緣-臺地帶東西向連井測井曲線
研究發(fā)現(xiàn),區(qū)內(nèi)發(fā)育有臺緣礁、臺內(nèi)二排礁和三排礁,但受地震分辨率的影響,臺內(nèi)二、三排礁寬度較小、內(nèi)部刻畫困難,針對這一問題,本文利用基于連續(xù)小波變換的拓頻方法,對地震數(shù)據(jù)進行了拓低頻處理,提高同相軸連續(xù)性。由于波傳播過程中不同頻帶信號都有所損失,導(dǎo)致頻帶變窄,對此,本文提出將地震信號分解到連續(xù)小波域,優(yōu)選小尺度小波信號(低頻信號),并進行適當?shù)难a償,具體思路如下。
首先,對地震數(shù)據(jù)進行連續(xù)小波變換(CWT)處理,變換公式如下:
(1)
式中:W(τ,s)為小波系數(shù);f(t)為原始信號;s為縮放參數(shù),控制小波函數(shù)的伸縮;τ為位置參數(shù),控制小波函數(shù)的平移;Ψ*為小波母函數(shù)。
其次,提取分解后的各個尺度小波信號的極值(極大值和極小值),利用“極值信號”與對應(yīng)尺度小波信息重構(gòu)原始信號。
(2)
最后,將低頻諧波信息進行能量增強,變換回時間域,并加入到原始數(shù)據(jù)中,這樣即可拓展地震數(shù)據(jù)的頻帶寬度。
(3)
式中:CΨ為小波的可容許條件;Ψ為小波函數(shù)。
圖3給出了地震數(shù)據(jù)拓低頻處理前、后的剖面。由圖3可以看出,經(jīng)過拓低頻處理后,礁灘體內(nèi)部丘狀反射得以體現(xiàn)(藍色箭頭),提高了礁灘體的識別精度。
圖3 地震數(shù)據(jù)拓低頻處理前(a)、后(b)地震剖面
2.1.1 巖石學特征
五百梯地區(qū)長興組巖石類型可劃分為灰?guī)r、白云巖及過渡類型三大類,其中灰?guī)r類包括礁灰?guī)r、顆粒灰?guī)r、泥晶灰?guī)r等,造礁生物以苔蘚蟲、海綿為主,少量珊瑚,附礁生物包括介形蟲、棘皮、雙殼等。云巖類包括礁云巖、顆粒云巖、晶粒云巖、泥晶云巖等。溶蝕孔洞與礁云巖、顆粒云巖、晶粒云巖關(guān)系密切,且多以未充填為主(圖4)。
圖4 五百梯地區(qū)長興組典型巖石類型巖心圖版a 生屑灰?guī)r; b 生屑云巖; c 泥晶灰?guī)r; d 生物礁云巖; e 生物礁灰?guī)r; f 中晶云巖; g 灰質(zhì)粉晶云巖; h 含云砂屑灰?guī)r; i 細晶棘屑云巖
2.1.2 沉積相發(fā)育特征
晚二疊世長興組沉積時期,川東地區(qū)為碳酸鹽巖臺地、海槽、臺地邊緣沉積體系,根據(jù)研究區(qū)巖石學、沉積構(gòu)造、測井相及區(qū)域格局分析,研究區(qū)可以劃分為臺地邊緣、前緣斜坡和開闊臺地3個相,同時受海平面變化的影響,發(fā)育了多個礁、灘相沉積旋回(圖5)。研究區(qū)臺地邊緣及臺內(nèi)高部位為礁灘體的集中發(fā)育區(qū),其中臺緣帶長興組中上部廣泛發(fā)育生物礁,連續(xù)性較好。橫向上,礁相主要發(fā)育在長三段中上部,自西向東表現(xiàn)為臺緣帶橫向連續(xù)性好、分布穩(wěn)定,臺內(nèi)厚度薄、期次多。礁灘體自長二期開始發(fā)育,垂向上發(fā)育多個灘間海-粒屑灘的垂向變淺旋回,開闊臺地內(nèi)灘間低能沉積,礁欠發(fā)育,單個礁灘體規(guī)模小,套數(shù)增多(圖6)。總體看,有利相帶位于上部的長三段,橫向上發(fā)育在臺緣及臺內(nèi)局部高地貌部位,且臺緣外帶優(yōu)于臺緣內(nèi)帶。
圖5 TD11井第1筒心沉積微相垂向發(fā)育序列
2.2.1 單井分析
圖7顯示了長興組典型井合成記錄。由圖7可見,斜坡相帶內(nèi)長興組頂為波峰反射,中強振幅、連續(xù)性好、時間厚度小;臺緣相帶內(nèi)長興組厚度明顯增大,長興組頂為波谷反射,能量明顯變差,連續(xù)性變差,長興組內(nèi)部同相軸增多、存在復(fù)波。
圖7 長興組典型井合成記錄a 斜坡相合成記錄; b 臺緣相合成記錄
將工區(qū)內(nèi)各井合成記錄進行統(tǒng)計并分析后發(fā)現(xiàn):位于臺緣相帶內(nèi)的井,長興組地層厚度為157~270m不等,平均厚度為218m,長興組頂能量弱、存在復(fù)波,部分井存在空白、雜亂反射;位于臺地二排礁的井,地層厚度為160~300m不等,平均厚度為230m,長興組頂能量弱,存在復(fù)波;位于斜坡相帶內(nèi)的井,地層厚度為92~133m不等,平均厚度95m,時間厚度明顯變小,長興組頂能量強,連續(xù)性好??傮w看,無論是臺緣還是臺內(nèi)發(fā)育礁灘的部位,長興組地震反射特征表現(xiàn)為頂部能量明顯變?nèi)酢⒖瞻追瓷?、連續(xù)性差,地層明顯增厚。
2.2.2 剖面分析
過井剖面可印證上述分析的反射特征(圖8)。圖9顯示了長興組臺內(nèi)二排礁地震反射特征及生長演化模式。從圖9a可以清晰地劃分出礁核和礁翼的分布特征。礁核位于古地貌高部位,儲層較發(fā)育,速度較低,振幅變?nèi)?是最有利的勘探區(qū)域,而礁翼分布于高地貌的兩側(cè),礁欠發(fā)育,灘體較發(fā)育,地震反射不明顯。綜上研究,建立了研究區(qū)礁灘體生長演化模式(圖9b),該演化模式顯示造礁生物通常沿砂屑灘、生屑灘等地貌高點建造,并通常經(jīng)歷定殖期、拓殖期、泛殖期等達到繁盛,形成較明顯的正地形。隨著礁體進一步生長及海平面變化,可能因出露海平面之上或水體過深而進入消亡期,生物礁一經(jīng)暴露于海平面之上,常受到大氣水溶蝕或混合水白云化作用而形成儲層。結(jié)合以上地震地質(zhì)認識,識別出礁灘體地震有利相帶50余km2。
圖8 長興組地震反射特征a 東西向地震剖面; b 東西向地震測線位置
圖9 長興組臺內(nèi)二排礁地震反射特征(a)及生長演化模式(b)
古地貌是控制長興組儲層發(fā)育的關(guān)鍵因素。結(jié)合研究區(qū)地質(zhì)特征及前人研究基礎(chǔ),本文選擇殘厚法繪制了研究區(qū)古地貌圖,分別以龍?zhí)督M底和長興組底進行層拉平計算,繪制各時期古地貌圖。圖10給出了研究區(qū)主要沉積期的古地貌。由圖10a可見,長興組沉積前研究區(qū)為緩坡沉積,水深自西向東逐漸變淺。受斷層影響,緩坡內(nèi)可見局部沉積高地貌區(qū)域,沉積地貌高部位主要沿北西—南東向和北東—南西向展布,以前者為主,其中,TD7-TD17-H4一線地貌較為明顯,奠基了長興組海槽的雛形。該期沉積地貌分布與長興組沉積前斷層分布密切相關(guān),也為后期礁灘體生長建造提供了基礎(chǔ)。
圖10 研究區(qū)主要沉積期古地貌a 龍?zhí)督M古地貌; b 長興組古地貌; c 長一段+長二段古地貌; d 長三段古地貌
長興組古地貌(圖10b)反映了長興組沉積模式及礁灘體建造,受礁灘體建造的影響,臺地邊緣沿TD002-1-TD021-5井一線及以東地區(qū)分布,整體呈現(xiàn)出北窄陡、南寬緩的特征,中部地區(qū)臺地邊緣水體淺、能量高,礁灘體發(fā)育,形成礁灘復(fù)合體,連片分布,以垂向加積作用為主。南部區(qū)域水體相較于中部較深,能量低,礁灘體建造相對較弱,礁灘體在長興期具有向內(nèi)帶遷移的特征,礁灘體間廣泛分布低能的灘間海、礁后瀉湖沉積。整體上臺緣外帶比臺緣內(nèi)帶更適于礁灘建造,形成的沉積地貌更高,水體更淺。開闊臺地內(nèi),整體以相對低能的灘間海沉積為主,礁灘復(fù)合體發(fā)育相對富集,形成二排礁分布。二排礁相較于臺緣礁更窄、更陡,二排礁復(fù)合體整體為北西—南東向分布,是在長興組沉積前的高地貌基礎(chǔ)上建造和演化。臺地內(nèi)的其它區(qū)域偶見點礁、馬蹄型礁。長興組由于礁灘體的生長建造,形成了鑲邊明顯的臺地邊緣及臺內(nèi)礁灘復(fù)合體分布帶,但受水體深度、沉積前地形及斷層的影響,不同區(qū)域礁灘體發(fā)育規(guī)模、分布范圍具有較大差異。
進一步研究發(fā)現(xiàn),長興組早中期高地貌格局已經(jīng)形成(圖10c),長三時期生物礁繼續(xù)發(fā)育(圖10d),形成了研究區(qū)現(xiàn)今的長興組生物礁格局。圖10d顯示晚期生長的礁體主要集中在臺緣外帶,臺緣內(nèi)帶和二排礁發(fā)育相對較少,這意味著優(yōu)質(zhì)儲層將會集中在臺緣外帶,與區(qū)內(nèi)高產(chǎn)井分布一致,也為臺內(nèi)礁灘儲層的勘探提供了方向。
隨著海平面下降及礁灘體營建,地貌高部位出露海平面,接受大氣淡水溶蝕及混合水白云化作用,形成礁灘體孔洞型、孔隙型儲層。長興組末期,受區(qū)域性暴露影響,長興組頂部受近地表巖溶影響形成溶洞型、縫洞型儲層,部分被充填。埋藏期,地層水及深部熱液流體沿斷層、裂縫及先期孔洞進入地層,進一步改善了儲集性能(圖11)。顯然古地貌是控制生物礁發(fā)育的最關(guān)鍵因素,控制著長興組有利相帶的分布,古地貌的高部位是勘探的有利區(qū)域,相帶變化的坡折帶是發(fā)育儲層的最佳部位,而且臺緣外帶(靠近斜坡帶)礁灘儲層優(yōu)于臺緣內(nèi)帶(靠近臺地一側(cè)),這為礁灘儲層預(yù)測提供了方向。
圖11 五百梯地區(qū)長興組礁灘體發(fā)育模式
針對礁灘儲層與古地貌的關(guān)系,選擇了工區(qū)內(nèi)已鉆井資料進行分析(表1)。表1中測試產(chǎn)量在20×104m3/d以上的高產(chǎn)井基本分布在臺緣外側(cè),位于古地貌高部位,長興組頂同相軸都存在雜亂、空白反射等特征。而臺緣內(nèi)側(cè)的TD75井儲層不發(fā)育,該井位于古地貌低部位。發(fā)育于臺緣中間的TD2井雖有16.70m的儲層,但位于古地貌低部位,儲層物性差,未能獲得規(guī)模產(chǎn)量。值得一提的是位于臺地內(nèi)二排礁中間的TD002-X18井,雖只有5.30m的儲層,但位于臺內(nèi)古地貌高部位,獲得了16.90×104m3/d的測試產(chǎn)量,預(yù)示著臺內(nèi)二排礁同樣具有勘探潛力,且儲層發(fā)育模式與臺緣相似。
表1 已鉆井特征分析
以上分析結(jié)果進一步佐證了古地貌是控制礁灘儲層發(fā)育的關(guān)鍵因素的認識,古地貌高部位后期容易發(fā)生白云化作用,為礁灘儲層改造提供了有利條件。但白云化作用也存在橫向非均質(zhì)性,從而導(dǎo)致長興組頂呈現(xiàn)出振幅變化大、連續(xù)性差的反射特征??傮w來看,礁灘儲層發(fā)育在臺緣和臺內(nèi)二排礁相帶內(nèi),地震反射特征較為明顯,且外側(cè)優(yōu)于內(nèi)側(cè),古地貌高部位才可發(fā)育較好的儲層。
從以上研究來看,五百梯地區(qū)礁灘儲層橫向上集中在臺緣礁和臺內(nèi)二排礁區(qū)域,縱向上主要分布在長興組頂部的長三段和長二段上部區(qū)域。為了進一步細化儲層的展布特征,我們結(jié)合地震反演和油氣檢測等方法進行了長興組儲層預(yù)測。首先利用測井資料分析出工區(qū)內(nèi)儲層測井響應(yīng)特征:速度5400~6300m/s、波阻抗小于17000m/s·g/cm3、孔隙度大于1.5%。并根據(jù)井數(shù)據(jù)建立了波阻抗、孔隙度和速度之間的擬合關(guān)系,進行長興組儲層厚度預(yù)測。預(yù)測結(jié)果顯示,厚度大的儲層發(fā)育在臺緣帶的中南部,該區(qū)域典型井如TD002-11井獲氣16.90×104m3/d,儲層預(yù)測厚度為51.60m,預(yù)測吻合度高。臺內(nèi)二排礁儲層厚度穩(wěn)定,優(yōu)勢儲層分布在工區(qū)的中南部,為二排礁有利勘探區(qū)域。
綜上分析,依據(jù)“一礁一藏”的勘探思路,劃分了研究區(qū)一類和二類有利區(qū)(圖12中青色和綠色區(qū)域)。表2給出了長興組有利勘探區(qū)的劃分標準。臺緣外側(cè)中部古地貌高部位容易產(chǎn)生白云化作用,高孔隙度、低阻抗分布,為工區(qū)內(nèi)的最有利區(qū)域,臺緣內(nèi)側(cè)的局部古地貌高部位也可劃分為一類有利區(qū)。臺內(nèi)二排礁中南部儲層分布較穩(wěn)定,與臺緣外帶的儲層特征一致,綜合分析該區(qū)同為一類有利區(qū)。
圖12 五百梯區(qū)塊長興組礁灘儲層綜合評價結(jié)果
1) 研究發(fā)現(xiàn)長興組早中期高地貌格局已經(jīng)形成,長三時期礁灘體繼續(xù)發(fā)育,形成了研究區(qū)現(xiàn)今的長興組礁灘發(fā)育格局,晚期生長的礁灘體主要集中在臺緣外帶及臺內(nèi)二排礁的高部位。
2) 長興組沉積末期古地貌整體呈現(xiàn)斜坡相帶低,臺緣及臺內(nèi)二排礁局部高的特征,高地貌是礁灘體發(fā)育的主要區(qū)域,其頂部的礁核發(fā)育區(qū)為儲層發(fā)育的最有利部位。
3) 古地貌控制下的礁灘體才是儲層發(fā)育的主體部位,研究區(qū)最有利的儲層發(fā)育區(qū)位于臺緣外帶及臺內(nèi)二排礁部分區(qū)域。
4) 古地貌是控制礁灘發(fā)育的關(guān)鍵因素,控制著長興組儲層分布。高地貌是勘探的最有利區(qū)域,這一現(xiàn)象在地震、地質(zhì)資料中都有體現(xiàn),為臺內(nèi)礁灘儲層的勘探提供了依據(jù)及借鑒意義,也意味著臺內(nèi)局部高地貌地區(qū)仍然發(fā)育了規(guī)模性礁灘儲層。