川西深層碳酸鹽巖巖溶儲層預測

杜勁松, 蘭雪梅, 黃天俊, 張連進, 文 雯, 王天琪, 蘇世為, 溫夢晗

(1.北京陽光杰科科技股份有限公司, 北京 100192； 2.西南油氣田分公司勘探開發(fā)研究院, 成都 610041)

碳酸鹽巖巖溶儲層是中外含油氣盆地中重要的油氣儲層之一,中國大型海相油氣田大都發(fā)現(xiàn)有巖溶型油氣藏,儲量大產(chǎn)量高,成為中國碳酸鹽巖油氣勘探的熱點領域。碳酸鹽巖巖溶儲層的發(fā)育及分布規(guī)律主要受構造、古地貌、巖性和巖相等因素控制,其預測方法基本上以古地貌恢復、地震響應模式、地震屬性等定性分析方法為主,定量預測處于探索階段。

川西地區(qū)深層下二疊統(tǒng)茅口組在多個地質構造中獲得勘探突破,多口井在巖溶體儲層中獲得高產(chǎn)工業(yè)天然氣氣流,揭示了川西茅口組良好的巖溶體儲層天然氣勘探前景。由于茅口組的勘探程度較低,資料有限同時巖溶體儲層埋深大,非均質性強,分布發(fā)育受控影響因素多等原因,增大了巖溶儲層地震預測和評估的工作難度[1-3],目前未能形成一套針對茅口組巖溶儲層進行有效地震預測的方法,制約了下一步川西深層下二疊統(tǒng)油氣勘探的工作。川西九龍山地區(qū)雖然茅口組碳酸鹽巖儲層為復雜非均質性強的巖溶性油氣儲層,但其沉積特征和地球物理響應特征都具有一定的規(guī)律性[4]。在鉆井和地震資料分析的基礎上,綜合運用層序地層學方法對研究區(qū)域進行地層構造精細解釋,開展巖溶體上覆地層反射結構的分析,茅二段模糊波形聚類法地震相劃分,茅口組剝蝕后殘余地層厚度法微古地貌恢復,茅二段巖溶體的地震屬性特征分析,對巖溶體進行綜合定性分析,總結巖溶體的分布規(guī)律。定量預測采用兩步法,首先用伽馬反演方法預測茅二段生屑灰?guī)r,并在此基礎上,進一步利用多屬性儲層孔隙度反演進行巖溶體儲層的精細定量預測,提高了預測精度。通過定性-定量預測相結合的巖溶體儲層預測方法,并結合地震曲率裂縫預測的結果,綜合分析、劃分出巖溶體儲層有利的分布區(qū),預測結果經(jīng)實鉆驗證吻合度高,對川西茅口組巖溶氣藏深化勘探開發(fā)具有重要指導意義。

圖1 龍16井茅二段測井曲線及儲層FMI響應特征Fig.1 Logging and reservoir FMI characteristics of P1m2 of L16

1 茅口組古巖溶識別標志

目前九龍山地區(qū)下二疊統(tǒng)茅口組地層已鉆井9口,獲高產(chǎn)工業(yè)性氣井3口,是天然氣滾動勘探開發(fā)增儲上產(chǎn)的重點層系。川西地區(qū)茅口組碳酸鹽巖沉積期處于開闊臺地海環(huán)境,地層由老到新可分為茅一段至茅四段。早二疊世末期,受東吳構造運動的影響,茅口組碳酸鹽巖地層廣泛抬升,出露的碳酸鹽巖遭受長期的風化剝蝕和大氣淡水的淋濾溶蝕作用,九龍山地區(qū)茅口組地層剝蝕至茅二段,發(fā)育有溶蝕殘丘和巖溶體,其發(fā)育程度、分布規(guī)模和連通性主要受古地貌和裂縫的控制[5]。通常可以利用鉆錄井顯示、測井響應特征、地震反射特征、動態(tài)資料分析等方法對茅口組是否存在巖溶體進行判別。

1.1 鉆錄井響應

茅口組碳酸鹽巖遭受古風化巖溶作用,侵蝕后會產(chǎn)生大量的溶蝕縫洞。鉆井過程中,若鉆遇到縫洞發(fā)育區(qū)時,會直接產(chǎn)生井漏、放空等現(xiàn)象,這是判定巖溶體存在的重要標志之一[6]。據(jù)統(tǒng)計,川西北部地區(qū)鉆井過程中,茅口組累計9口井有直接錄井顯,出現(xiàn)放空、井漏等現(xiàn)象,分布深度為茅二段上部至茅口組頂部,具有典型的巖溶型儲層鉆錄井響應特征。

1.2 測井響應

茅口組巖溶儲層在各種常規(guī)測井曲線上,主要表現(xiàn)為高聲波差和較高中子孔隙度,低自然伽馬、中低電阻和低密度的“兩高三低”測井曲線響應特征[7],儲層電性綜合特征為:聲波時差>49 μs/ft,中子>1.5%,自然伽馬<20 API,密度<2.75 g/cm3,電阻率為高值背景上相對低值,一般100～10 000 Ω·m。

巖溶儲層裂縫發(fā)育段,在常規(guī)曲線上顯示為電阻率降低及雙側向正差異等特征,在成像測井(FMI)上表現(xiàn)呈不規(guī)則正弦曲線狀且具高導特征。龍16井茅二段測井解釋儲層厚9.2 m,平均孔隙度2.1%,儲層段常規(guī)曲線特征為聲波時差為49.7～50.7 μs/ft,自然伽馬為24.8～35 API, 中子為1.4～2.1%,電阻率為高阻背景上相對降低,FMI上顯示有溶蝕孔洞且低角度裂縫發(fā)育,如圖1所示。

1.3 地震響應

圖2 龍004-X1井茅口組巖溶儲層地震響應特征Fig.2 Seismic reflection of Karst reservoir in Maokou Formation of L004-X1

九龍山地區(qū)茅口組溶蝕殘丘在地震剖面上具有較為明顯的響應特征,主要表現(xiàn)為孤立突起,突起部分振幅逐漸變?nèi)?常與層間小斷層伴生,平面上呈散點狀分布。雖然巖溶體縫洞的地震響應特征不如溶蝕殘丘特征明顯,但也表現(xiàn)出一定的規(guī)律性,主要反射特征表現(xiàn)為“低頻弱振幅較雜亂反射”、地層時差較大,平面上呈片狀分布。當巖溶體儲層縫洞發(fā)育時,由于波阻抗降低,與上覆地層的阻抗差下降,導致反射系數(shù)降低,從而形成相對較弱的地震反射特征。巖溶體儲層縫洞中通常充滿流體,特別是含氣后,會使地震的高頻成分被吸收,導致頻率下降,從而形成“低頻弱振幅”的反射特征[8]。通常茅口組巖溶體儲層內(nèi)部的縫洞分布是無序且形態(tài)雜亂的,當縫洞發(fā)育規(guī)模較大時,會出現(xiàn)“較雜亂反射”現(xiàn)象,同相軸表現(xiàn)為扭曲狀。龍004-X1井,在鉆至茅口組時出現(xiàn)井漏現(xiàn)象,測試獲天然氣111.65×104m3/d,巖溶體氣層頂界為弱波峰,底界為弱波谷,同相軸存在扭動現(xiàn)象,如圖2所示。

1.4 動態(tài)資料分析

與龍16和龍004-X1兩井相鄰的龍4井,曾經(jīng)間斷放噴528次,累計放噴量達2 171×104m3,每次放噴后關井壓力恢復迅速而沒有衰減,說明該井茅口組巖溶體連通范圍內(nèi)能量充足,儲層分布范圍較大,從另一角度佐證了九龍山主體構造高部位茅二段存在一個成規(guī)模縫洞系統(tǒng)的可能,茅口組頂部地震裂縫預測平面圖上,龍4井位于團塊狀裂縫發(fā)育區(qū)內(nèi)(圖3)。

圖3 九龍山構造茅口組頂部裂縫預測Fig.3 Prediction of fractures at the top of Maokou Formation in Jiulongshan structure

2 巖溶體上覆地層反射結構

通過精細標定,已鉆井中,茅二段巖溶儲層位于龍?zhí)督M底界波峰同相軸的下沿。聯(lián)井飛一段底界層拉平地震剖面上,如圖4所示,清楚顯示出,當龍16井區(qū)茅口組頂部巖溶儲層發(fā)育時,其上覆龍?zhí)督M地層多出一個波峰同相軸,且其振幅有強弱變化；而在龍17井區(qū),其上覆龍?zhí)督M地層變薄,未出現(xiàn)波峰反射,龍17井鉆井揭示,茅口組頂部巖溶儲層不發(fā)育。巖溶體上覆地層反射結構的不同,也說明這兩個井區(qū)在茅口組上部地層剝蝕后的沉積環(huán)境是不同的,從而間接反映出茅口組巖溶儲層發(fā)育的相帶。

圖4 茅口組巖溶體上覆地層多出一個波峰反射軸(飛一段底界層拉平)Fig.4 An additional seismic event in the overlying strata of the Maokou karst (T1f1 flattened)

圖5 龍?zhí)督M多出波峰同相軸分布Fig.5 Distribution of the additional seismic event of Longtan Formation

計算飛一段底界與龍?zhí)督M底界的時差,如圖5所示,可以清楚地反映出龍?zhí)督M多出的波峰同相軸的分布,進而定性指示茅口組頂部的儲層發(fā)育區(qū)。井標定表明,波峰同相軸尖滅線西部區(qū)域(綠色-紅色-黃色區(qū))上覆地層逐漸變厚,多出一個波峰反射軸,鉆井揭示茅二段巖溶儲層較發(fā)育；尖滅線東部區(qū)域(藍色區(qū))上覆地層厚度薄,無波峰同相軸,茅二段巖溶儲層不發(fā)育。

圖6 茅二段地震相圖Fig.6 Seismic facies of P1m2 member

3 地震相劃分巖溶體分布范圍

利用模糊地震波形分類技術對茅二段地層進行地震相預測,得到該段地震相圖,如圖6所示,巖溶殘丘呈孤立的黃色點狀分布,主要位于研究區(qū)的西半部,巖溶體發(fā)育區(qū)位于研究區(qū)的西部,為黃色北西西向長條狀連片分布。地震相圖中巖溶殘丘和巖溶體的分布位置也位于上覆地層有利反射結構分布區(qū)內(nèi)。

4 微古地貌恢復

茅口組古地貌形態(tài)對于研究古巖溶儲層的分布具有重要指導意義。古地貌形態(tài)恢復的方法較多,最簡便易用的方法是印模法和殘余厚度法,但須根據(jù)研究目標的具體地質特征,采用相應的方法。采用殘余厚度法恢復茅口組地層剝蝕后的古地貌形態(tài)是適合的。殘余厚度法的基本原理,是假定所研究地層剝蝕前的原始沉積厚度變化不大,且該地層之下存在一個等時的基準面,用該基準面與研究層頂界面之間的殘余厚度來表征古地貌的高低變化,殘余厚度越大,古地勢越高[8]。下二疊統(tǒng)棲霞組地層沉積時,整個四川盆地為開闊臺地沉積環(huán)境,沉積厚度相差不大,因此將棲霞組頂界面(即茅口組底界)作為一個相對等時基準面,棲霞組頂界與茅口組頂界之間的厚度差值,用于表征茅口組剝蝕后的古地貌形態(tài)。

實際研究中,在精細層位解釋的基礎上,用棲霞組頂界層位減去茅口組頂界層位,得到茅口組殘余厚度時差,用于茅口組剝蝕后的古地貌恢復,如圖7所示。龍104、龍16、龍4、龍004-X1井在茅口組剝蝕之后,殘余地層厚度大,整體處于古地貌高部位,后期巖溶作用較為劇烈,為巖溶發(fā)育區(qū)。龍17井也位于古地貌高部位,由于茅口組剝蝕之后的沉積環(huán)境與龍16井區(qū)不同,巖溶作用較弱,巖溶儲層并不發(fā)育。

古地貌圖上,巖溶殘丘呈孤立的黃紅色點狀分布,龍16井區(qū)巖溶分布在古地貌高部位,形態(tài)呈片狀,分布范圍與地震相中的黃色條帶一致,但刻畫的較其精細。

5 巖溶體刻畫

已鉆井茅口組頂部巖溶儲層標定在龍?zhí)督M底界波峰反射軸的下沿,高產(chǎn)氣井龍16、龍4、龍004-X1井巖溶儲層發(fā)育,地震剖面上,龍?zhí)督M底界波峰振幅較弱,龍?zhí)?井巖溶儲層不太發(fā)育,為強震幅。

針對巖溶體儲層的地震響應特征,通過大量地震屬性提取和分析,優(yōu)選茅二段弧長屬性對茅口組巖溶體進行刻畫。

5.1 茅二段弧長屬性刻畫巖溶體

弧長(arc length)屬性是指時窗內(nèi)每個地震道的波形伸展后的總長度,它是一種頻率與振幅的混合屬性,可用于區(qū)分具有相同振幅特征,但有高低頻變化的地層情況,對地層非均質性、縫洞發(fā)育等反射特征反映較好。沿龍?zhí)督M底界層位,向上開時窗10 ms,向下開時窗20 ms,提取弧長屬性,結合井點斷溶體的弧長屬性值域范圍,截取低值域800～5 750用于刻畫斷溶體的分布,如圖8所示。茅二段弧長屬性刻畫的斷溶體分布在地震相帶和古地貌指示的斷溶體發(fā)育區(qū)范圍內(nèi),與已鉆井揭示的斷溶體非常吻合,可信度高。

5.2 巖溶體地震屬性特征

已鉆井茅口組頂部巖溶儲層標定在龍?zhí)督M底界波峰反射軸的下沿,高產(chǎn)氣井龍16、龍4、龍004-X1井巖溶體發(fā)育,龍?zhí)督M底界波峰振幅較弱,龍?zhí)?井巖溶體不太發(fā)育,為強振幅。沿龍?zhí)督M底界層位,上下時窗各8 ms,提取最大波峰振幅屬性,如圖9所示。巖溶體地震屬性與茅二段弧長屬性刻畫的巖溶體分布形態(tài)基本一致,巖溶體內(nèi)的地震屬性有強弱變化,可以用于指示巖溶體儲集性能的優(yōu)劣。

巖溶體最大波峰地震屬性圖上,龍16井區(qū)位巖溶體弱振幅分布區(qū)。巖溶體弱振幅區(qū)與茅口組頂部裂縫發(fā)育區(qū)吻合度較高,如圖10所示。

6 巖溶儲層定量預測

在茅口組巖溶體定性分析的基礎上,嘗試定量預測巖溶儲層的厚度。九龍山地區(qū)現(xiàn)有九口井鉆遇茅口組地層。井上揭示茅二段殘留生屑灘體厚度不大,最大厚度約30 m,自然伽馬(GR)一般低于28 API,波阻抗(IMP)在15 000～17 800 m/s*g/cm3,其中巖溶儲層波阻抗15 000～16 600,致密層波阻抗16 600～17 800。茅口組灘體中還夾雜一定厚度的泥灰?guī)r,其自然伽馬大于28 API,波阻抗在13 000～17 500,如圖11所示。巖溶儲層波阻抗與泥灰?guī)r波阻抗完全重疊,無法區(qū)分儲層和非儲層。

研究中,針對茅口組巖溶儲層的定量預測采用兩步法,首先通過伽馬反演預測茅二段灘體的分布,在此基礎上,再利用多屬性孔隙度反演來進行儲層預測。

采用特征反演法對茅二段灘體進行伽馬反演。通過對比,伽馬反演剖面與井曲線吻合度較好,如圖12所示。提取茅二段滿足伽馬值小于28的采樣點數(shù),獲得灘體的時間厚度分布,再與井點灘體厚度回歸,得到茅二段灘體厚度分布。平面上生屑灘體厚度分布(圖13)符合地質規(guī)律。

其次,在伽馬反演的基礎上,利用多屬性多權值孔隙度反演來對巖溶儲層進行預測。原理是假設目標曲線上某采樣點值是由地震屬性道上對應點上下多個采樣點共同作用的結果,通過建立目標曲線與多屬性多樣點的非線性方程來達到預測目的[9]。多屬性多權值孔隙度預測計算過程中,屬性個數(shù)、屬性組和算子長度是關鍵參數(shù)。采用逐步回歸方法把相關系數(shù)最大且誤差最小的屬性形成屬性組各屬性之間是線性獨立的[10]。

按上述方法進行孔隙度反演,所得到的反演剖面與井曲線吻合較好,如圖14所示,通過計算茅二段滿足孔隙度>1.8%的數(shù)據(jù)采樣點的累加時間厚度,與井點儲層厚度進行變差函數(shù)回歸[11],得到茅二段儲層厚度分布(圖15),儲層分布與定性儲層預測結果匹配度較高。

圖11 多井茅二段波阻抗-伽馬交會圖Fig.11 IMP-GR crossplot of P1m2 member of multiple wells

研究區(qū)茅二段儲層厚度可達26 m,溶蝕殘丘儲層呈孤立點狀分布,巖溶體儲層呈北西西向條帶狀展布,主要分布在龍16井區(qū)及其西北部,龍17井區(qū)巖溶儲層不太發(fā)育。

7 有利開發(fā)區(qū)綜合預測

通過對茅口組的構造解釋,古地貌恢復、地震屬性分析、儲層及裂縫預測等定性、定量成果的相互匹配,劃分出開發(fā)有利區(qū)。茅口組巖溶儲層有利開發(fā)區(qū)主要分布在九龍山構造的南翼。龍?zhí)督M底部波峰振幅屬性、裂縫發(fā)育區(qū)、構造位置是部署井位的關鍵參數(shù)。

圖12 聯(lián)井伽馬反演剖面Fig.12 Cross-well Gamma inversion section

圖13 茅二段灘體厚度Fig.13 Biological beach thickness of P1m2 member

圖14 聯(lián)井孔隙度反演剖面Fig.14 Cross-well porosity inversion section

圖15 茅二段儲層厚度Fig.15 Reservoir thickness of P1m2 member

8 結論

(1)九龍山地區(qū)茅口組發(fā)育巖溶儲層,鉆錄井中出現(xiàn)明顯的放空、漏失現(xiàn)象,測井曲線表現(xiàn)為“兩高三低”特征,地震剖面上具有低頻、弱振幅、較雜亂、同相軸扭動、時差變大的反射特征。

(2)茅口組上覆龍?zhí)督M地層多出一個波峰反射軸,指示該區(qū)域茅二段可能有巖溶儲層發(fā)育,未出現(xiàn)波峰反射軸區(qū)域,茅二段巖溶儲層不發(fā)育。

(3)茅口組生屑灘相是巖溶儲層發(fā)育的基礎,受巖溶作用改造后的有利古地貌區(qū)是控制巖溶儲層分布的主要因素。殘余厚度法適于茅口組剝蝕后的古地貌恢復,能很好指示巖溶儲層的分布范圍,巖溶儲層分布在古地貌高部位。

(4)龍?zhí)督M底界最大波峰振幅屬性對茅二段巖溶儲層反應靈敏,古地貌高部位且弱振幅區(qū)是巖溶儲層發(fā)育的有利區(qū)。

(5)采用兩步法進行巖溶儲層定量預測,巖溶儲層分布特征與巖溶體定性預測結果匹配度較高。

(6)定性-定量預測相結合的巖溶儲層預測方法,對茅二段巖溶體儲層的刻畫更可靠、更精細。