徐兆輝 胡素云 王 露 趙文智 曾洪流
1中國石油勘探開發(fā)研究院,北京 100083 2德克薩斯大學(xué)奧斯汀分校經(jīng)濟(jì)地質(zhì)局,美國奧斯汀 78757
地震沉積學(xué)自正式被提出開始(Zengetal., 1998a, 1998b), 經(jīng)過與地震地貌學(xué)相結(jié)合(Posamentier, 2001), 其本質(zhì)內(nèi)涵在2004年被明確下來(Zeng and Hentz, 2004), 之后逐漸引起了國內(nèi)外地質(zhì)學(xué)家和地球物理學(xué)家的關(guān)注(董艷蕾等, 2008; 曾洪流, 2011)。目前廣為接受的定義為: 地震沉積學(xué)是地震巖性學(xué)和地震地貌學(xué)的有機(jī)結(jié)合, 是基于地震資料研究沉積巖和沉積史的一門學(xué)科(林承焰和張憲國, 2006; 林承焰等, 2007; 曾洪流等, 2012; 朱筱敏等, 2019)。近年來, 隨著地震采集處理技術(shù)和計(jì)算機(jī)運(yùn)算能力的提升, 地震沉積學(xué)在石油工業(yè)中的應(yīng)用潛力與日俱增。中國含油氣盆地普遍經(jīng)歷了復(fù)雜的沉積演化, 形成了不同的沉積環(huán)境和巖性組合。尤其是在西部復(fù)合含油氣盆地中, 這些不同的巖性組合蘊(yùn)含了豐富的油氣資源(陳永武和何文淵, 2004; 趙文智等, 2005)。對(duì)于油氣勘探而言, 沉積相研究和儲(chǔ)集層預(yù)測(cè)是2項(xiàng)至關(guān)重要的研究內(nèi)容。在碎屑巖、 碳酸鹽巖和混積巖3種巖性的沉積環(huán)境中, 其沉積相和儲(chǔ)集層特征各異, 導(dǎo)致地震反射千差萬別。但是3種環(huán)境中發(fā)育的沉積層均屬地震解釋意義上的薄層, 這為開展地震沉積學(xué)研究提供了必要性。筆者在研究四川盆地川中龍王廟組和嘉陵江組及準(zhǔn)噶爾盆地瑪湖凹陷克拉瑪依組過程中發(fā)現(xiàn), 采取地震沉積學(xué)相關(guān)技術(shù)的不同組合, 可以有效地預(yù)測(cè)3種環(huán)境的沉積相和儲(chǔ)集層分布。通過3個(gè)實(shí)例的解剖, 詳細(xì)介紹了如何根據(jù)3種沉積體系的具體情況, 因地制宜地利用地震沉積學(xué)充分挖掘地震資料中的地質(zhì)信息, 從而更好地研究沉積相和預(yù)測(cè)儲(chǔ)集層。
準(zhǔn)噶爾盆地瑪湖凹陷克拉瑪依組被作為沖積扇—扇三角洲的典型實(shí)例,沉積了巨厚的砂礫巖層。筆者研究發(fā)現(xiàn)在這種以扇—三角洲沉積為主體的環(huán)境中,局部發(fā)育曲流河沉積。如果能夠解釋這種曲流河與沖積扇“跨相并存”現(xiàn)象,準(zhǔn)確預(yù)測(cè)有利儲(chǔ)集層分布,既有益于完善沉積學(xué)理論,又有利于提高油氣勘探開發(fā)效率。在這一研究內(nèi)容中,主要應(yīng)用了地震沉積學(xué)中的地震地貌學(xué)和地震巖性學(xué)2個(gè)方面的理論技術(shù)開展研究。
準(zhǔn)噶爾盆地面積約13×104ikm2,是中國西部重要的含油氣盆地之一,其周緣被哈拉阿拉特山、青格里底山、克拉美麗山、博格達(dá)山、伊林黑比爾根山和扎伊爾山所圍限(張國俊,1995;陳業(yè)全和王偉鋒,2004),盆地包括6個(gè)二級(jí)構(gòu)造單元,分別為: 烏倫古坳陷、陸梁隆起、中央坳陷、南部斷褶帶、西部隆起和東部隆起(圖 1-a)。近年來,針對(duì)砂礫巖油氣藏的勘探持續(xù)獲得重大突破,掀起了新一輪勘探熱潮(唐勇等,2019;楊帆等,2019)。其中,隸屬于中央坳陷的瑪湖凹陷是近年來準(zhǔn)噶爾盆地勘探重點(diǎn)地區(qū)。在地質(zhì)演化過程中,瑪湖凹陷受到附近沖斷活動(dòng)的影響,經(jīng)歷了3個(gè)主要的構(gòu)造旋回,分別為泥盆紀(jì)至早石炭世、晚石炭世至三疊紀(jì)、侏羅紀(jì)至第四紀(jì)。
前人的研究已經(jīng)指出, 瑪湖凹陷二疊紀(jì)至三疊紀(jì)的沉積相類型包括: 沖積扇、 扇三角洲、 河流三角洲和湖泊(丘東洲, 1994; 何登發(fā)等, 2018)。 在前石炭系結(jié)晶基底之上, 凹陷中沉積了從二疊系佳木河組到第四系巨厚地層, 其中, 三疊系克拉瑪依組上段(T2k2)被劃分為5個(gè)砂層組(S1-S5), 總厚度150~200im, 地層單層厚度較薄, 介于10~60im之間, 上覆地層為白堿組(T3b)厚層泥巖。 本文的研究目的層為最頂部的第1個(gè)砂層組(S1), 鉆井和露頭資料揭示其巖性包括: 粗礫巖、 細(xì)礫巖、 砂礫巖、 砂巖、 粉砂巖和泥巖等。 這些巖性被認(rèn)為形成于辮狀河、 三角洲、 扇—三角洲及湖泊相(劉順生等, 1999; 何苗等, 2017), 未見關(guān)于S1中發(fā)育曲流河的報(bào)道。 作為全區(qū)分布的白堿灘組厚層泥巖之下的第1套砂礫巖粗粒沉積, S1的沉積相和儲(chǔ)集層發(fā)育情況是非常重要的一項(xiàng)研究內(nèi)容。
文中所采用的三維地震覆蓋面積約600ikm2,是面元為12.5im×12.5im的高密度地震資料,整體頻帶寬度為10~70iHz,主頻為30iHz,采樣間隔2ims。重點(diǎn)關(guān)注的S1砂層組所在的T2k2,地震主頻高達(dá)35iHz。鑒于目標(biāo)層平均層速度約為4300im/s,因此該地震資料在目標(biāo)層的縱向分辨率約為30im,品質(zhì)較高。
在地震覆蓋范圍內(nèi),鉆穿目的層S1的井總計(jì)35口,所有井都有常規(guī)測(cè)井曲線數(shù)據(jù)(包括GR、RT、AC和RHOB),其中,18口井有孔隙度測(cè)井曲線,14口井有巖心資料(圖 1-b)。這些井資料可以保證利用測(cè)井曲線和地震資料來識(shí)別和預(yù)測(cè)巖性與物性的準(zhǔn)確性。
圖 1 準(zhǔn)噶爾盆地構(gòu)造單元?jiǎng)澐旨把芯繀^(qū)位置(a)和區(qū)內(nèi)鉆測(cè)井分布情況(b)Fig.1 Location of study area and tectonic units in Junggar Basin(a)and wells distribution(b)
采用的研究方法包括地震地貌學(xué)和地震巖性學(xué)2種,前者主要用于定性恢復(fù)沉積相和砂體展布,后者主要用來定量表征儲(chǔ)集層。具體而言,在沉積相和砂體展布定性預(yù)測(cè)過程中,首先是利用層拉平技術(shù),將白堿灘組底界面(即克拉瑪依組頂界面)拉平,得到的圖件可視為S1沉積古地貌,然后對(duì)相位旋轉(zhuǎn)之后的地震數(shù)據(jù)體進(jìn)行分頻處理,分成3個(gè)分頻數(shù)據(jù)體并分別提取平均振幅屬性,將3幅分頻振幅屬性圖進(jìn)行RGB三色融合,最后結(jié)合古地貌和分頻融合屬性得到沉積相與砂體定性展布。在儲(chǔ)集層定量表征過程中,首先是在-90°相位數(shù)據(jù)體上提取目的層S1的一系列地震屬性,然后優(yōu)選與井點(diǎn)儲(chǔ)集層厚度相關(guān)性較好的地震屬性并對(duì)其進(jìn)行主因子分析,最后將前2個(gè)主因子與井點(diǎn)儲(chǔ)集層厚度擬合,把得到的擬合關(guān)系式應(yīng)用到前2個(gè)主因子之上,即可定量算得儲(chǔ)集層厚度分布。
圖 2 瑪湖凹陷西北部三疊系克拉瑪依組S1層分頻地震剖面特征(a、b、c)與分頻平均振幅屬性平面特征(d、e、f)(AA′位置見圖1-b)Fig.2 Characteristics of frequency-decomposed seismic profiles(a、b、c)and lateral average amplitude maps(d、e、f)of layer S1 of the Trassic Karamay Formation in northwestern Mahu sag (location of AA′ in Fig.1-b)
1.3.1 地震地貌學(xué)恢復(fù)沉積相與砂體展布
根據(jù)鉆井揭示的地層和巖性厚度分布和平均層速度統(tǒng)計(jì),將20im、30im和50im作為薄層、中層和厚層的厚度,據(jù)此計(jì)算其所對(duì)應(yīng)的調(diào)諧頻率為54iHz、36iHz和22iHz,對(duì)相位旋轉(zhuǎn)之后的數(shù)據(jù)體進(jìn)行分頻,得到分頻地震數(shù)據(jù)體(圖 2-a,2-b,2-c),然后提取目的層S1的平均振幅屬性,得到分頻振幅屬性(圖 2-d,2-e,2-f)。圖2中可以看出,隨著頻率的增加,剖面上地震同相軸逐漸變細(xì),目的層S1從低頻剖面上的1個(gè)粗波峰漸變成高頻剖面上的1個(gè)細(xì)波峰。從平均振幅屬性平面圖上看,低頻振幅屬性圖反映厚層大面積地分布在研究區(qū)中部,高頻振幅屬性圖反映了薄層主要散布在西北部、東部和東南部地區(qū)。
圖 3 瑪湖凹陷西北部三疊系克拉瑪依組S1層古地貌恢復(fù)結(jié)果(a)與沉積相平面分布圖(b)Fig.3 Recovered palaeotopography(a)and sedimentary facies distribution(b)of layer S1 of the Trassic Karamay Formation in Mahu sag
利用層拉平技術(shù)恢復(fù)的古地貌揭示,研究區(qū)內(nèi)主要發(fā)育3個(gè)古高地: 東北角、西南部和東南部(圖 3-a)。這3個(gè)地區(qū)的薄層、中層和厚層砂體均不甚發(fā)育,在振幅圖上表現(xiàn)為低振幅值(圖 2-d,2-e,2-f中的深藍(lán)色)。將3個(gè)分頻振幅屬性圖進(jìn)行RGB三色融合之后,生成砂/礫巖巖相定性平面分布圖(圖 3-b)。其中,紅色代表厚層、藍(lán)色代表薄層、綠色代表中層。厚層主要發(fā)育在由紅虛線勾勒出的研究區(qū)中部(區(qū)域C),呈扇形。中厚層則發(fā)育在由黑實(shí)線勾勒出的研究區(qū)西北部(區(qū)域A)和由黑虛線勾勒出的研究區(qū)中部(區(qū)域B)。
對(duì)比古地貌圖和巖相定性分布圖可見,中厚層砂礫巖主要分布在地貌低部位的高可容納空間地區(qū)(圖 3-a中的區(qū)域A和C),體現(xiàn)了典型的“溝谷控砂”規(guī)律。另外,利用地震地貌學(xué)基本原理,從巖相定性分布圖上可以看出,明顯存在2個(gè)沉積體系: 西北部的曲流河沉積體系(區(qū)域A)和中部的沖積扇沉積體系(區(qū)域C)。其中,沖積扇沉積體系的沉積物厚度明顯大于曲流河沉積體系,兩者平面距離最小處不足2ikm。將幾口典型井中S1段的測(cè)井曲線投影到圖中,特別是X723和FN401這2口井的巖性垂向疊置關(guān)系表現(xiàn)出典型的沖積扇反旋回和曲流河二元結(jié)構(gòu)的特點(diǎn),進(jìn)一步證實(shí)了2個(gè)沉積體系的存在。
圖 4 瑪湖凹陷西北部三疊系克拉瑪依組S1層過2個(gè)巖相的任意地震剖面與井標(biāo)定情況(BB′位置見圖1-b)Fig.4 Seismic profile with well calibration across two lithofacies of layer S1 of the Trassic Karamay Formation in Mahu sag (location of BB′ in Fig.1-b)
根據(jù)沃爾索相律(Middleton,1973),研究區(qū)目的層中的曲流河和沖積扇是不應(yīng)該緊鄰發(fā)育的。這一看似矛盾的巖相平面分布特征,可以通過精細(xì)的地震剖面標(biāo)定與觀察得到解釋(圖 4)。如前所述,巖相圖是將提取自S1的平均振幅屬性通過三色融合技術(shù)生成的。該過程中使用的時(shí)窗選取原則是白堿灘組厚層泥巖底界下的第1個(gè)波峰強(qiáng)反射,也就是白堿灘組底界之下約20ims(將第1個(gè)波峰包絡(luò)進(jìn)來)的范圍。油田勘探實(shí)踐發(fā)現(xiàn),該界面之下的第1套砂礫巖粗粒沉積具有很好的成藏條件。然而,從地震剖面可以看出,實(shí)際上白堿灘組底界面之下的第1個(gè)波峰并非嚴(yán)格等時(shí)。位于剖面中間的井(X723)中,巖性組合反映了典型的礫巖—泥巖—粉砂巖混雜的沖積扇扇中沉積。位于剖面兩側(cè)的2口井(FN407和M7)巖性旋回揭示,頂部的巖性旋回(S1a)為砂巖—粉砂巖—泥巖構(gòu)成的河流相特征,而真正的含有礫巖—泥巖—粉砂巖混雜的沖積扇沉積特征的旋回則出現(xiàn)在2口井的較深部(S1b),對(duì)應(yīng)于剖面中頂部波峰之下的第1套黑色波峰反射同相軸。也就是說,在側(cè)翼2口井所在剖面上,真正與中間鉆井沖積扇沉積等時(shí)的是2口井下部的同相軸(剖面上由藍(lán)色虛線勾勒出來的范圍),而非頂部的第1個(gè)波峰。
通過以上分析可得出結(jié)論,圖 3所反映的地震巖相圖實(shí)際上是2個(gè)時(shí)期沉積相的疊置。首先是在研究區(qū)中部(乃至全區(qū)范圍)沉積了沖積扇體系,之后,在其兩翼(西北部和東南部)發(fā)育了河流相沉積。
1.3.2 地震巖性學(xué)定量刻畫儲(chǔ)集層分布
利用主因子分析技術(shù)處理地震屬性系列,得到相應(yīng)的主因子。將前2個(gè)主因子與井點(diǎn)處實(shí)際儲(chǔ)集層厚度擬合,進(jìn)而計(jì)算出儲(chǔ)集層厚度分布圖(圖 5)。
圖 5 瑪湖凹陷西北部三疊系克拉瑪依組S1層儲(chǔ)集層厚度定量分布及單井油氣水產(chǎn)量Fig.5 Quantitative distribution of reservoir and the production of single well of layer S1 of the Trassic Karamay Formation in Mahu sag
圖 6 川中高磨地區(qū)安岳氣田井點(diǎn)位置(a)和綜合巖性柱(b)Fig.6 Location of wells(a)and comprehensive lithologic column(b)in Anyue Gas Field in Gaomo area of Central Sichuan Basin
目的層S1的儲(chǔ)集層厚度介于0~10im之間,在研究區(qū)內(nèi)廣泛分布,主要圍繞兩大區(qū)域: 河流相沉積區(qū)(點(diǎn)壩,區(qū)域A)和沖積扇相沉積區(qū)(扇緣,區(qū)域B)。盡管2個(gè)沉積體系中的儲(chǔ)集層平均厚度接近(6~8im),但是2個(gè)區(qū)域內(nèi)單井油氣產(chǎn)量卻差異明顯。與河流體系相關(guān)的鉆井(黑色箭頭所指)油氣產(chǎn)量明顯高于與沖積扇沉積體系相關(guān)的鉆井(白色箭頭所指)(圖5)。結(jié)合定性巖相恢復(fù)結(jié)果,不難解釋這一現(xiàn)象。區(qū)域A所圍限的曲流河相關(guān)儲(chǔ)集層大多形成于河流點(diǎn)壩沉積,其成分成熟度和結(jié)構(gòu)成熟度均應(yīng)高于區(qū)域B所圍限的扇緣沉積,可能導(dǎo)致同等厚度的河流點(diǎn)壩儲(chǔ)集層具有更好的孔滲條件,因而可以較好地解釋其中所鉆井的產(chǎn)能更高。
碳酸鹽巖沉積體系不存在異源碎屑物質(zhì)注入,川中高磨地區(qū)寒武系龍王廟組是典型代表。作為安岳氣田主力產(chǎn)層,龍王廟組的儲(chǔ)集層主要發(fā)育在白云巖地層中。龍王廟組白云巖是淺水高能開闊臺(tái)地環(huán)境形成的顆?;?guī)r,經(jīng)過多期白云石化作用所形成(張建勇等,2015),單一的白云巖排除了巖性對(duì)于物性預(yù)測(cè)的干擾。這種巖性特征為利用地震沉積學(xué)恢復(fù)古地貌,定性、定量預(yù)測(cè)儲(chǔ)集層分布提供了良好的試驗(yàn)區(qū),因而選取龍王廟組作為實(shí)例,研究地震沉積學(xué)在碳酸鹽巖體系沉積相和儲(chǔ)集層預(yù)測(cè)中的應(yīng)用。
高磨地區(qū)地處川中, 構(gòu)造上隸屬于川中古隆起平緩構(gòu)造區(qū)的威遠(yuǎn)至龍女寺構(gòu)造群。從平面上看, 高磨地區(qū)現(xiàn)今構(gòu)造總體為古隆起背景之上的1個(gè)鼻狀突起, 自西向北東方向傾伏。高磨地區(qū)實(shí)際上是由南北2個(gè)構(gòu)造圈閉構(gòu)成, 南部為高石梯構(gòu)造圈閉, 北部為磨溪構(gòu)造圈閉(圖 6-a)。其中, 后者是安岳氣田龍王廟組產(chǎn)層的主要分布地區(qū)。高磨地區(qū)沉積地層發(fā)育在前震旦系花崗巖基底之上, 經(jīng)過陡山沱組泥頁巖和砂礫巖等沉積之后, 在上震旦統(tǒng)燈影組沉積了藻云巖夾泥巖。下寒武統(tǒng)筇竹寺組和滄浪鋪組主要發(fā)育粉砂巖和頁巖間互沉積, 目的層龍王廟組白云巖地層厚度介于10~80im之間, 被上覆中寒武統(tǒng)高臺(tái)組砂泥巖所覆蓋(圖 6-b), 上寒武統(tǒng)洗象池群則以砂屑白云巖和鮞粒白云巖為特征(汪澤成等, 2016)。龍王廟組的白云巖為灘相白云巖, 儲(chǔ)集層的孔隙類型包括溶孔、 晶間孔和殘余粒間孔, 孔隙的存在降低了白云巖的速度和密度(亦即降低了波阻抗), 因此在地震資料上有所響應(yīng)。
圖 7 川中高磨地區(qū)龍王廟組沉積古地貌與同沉積走滑斷裂分布(a)和圣安德烈斯斷裂體系與底形分布(b)Fig.7 Distribution of topography and syndepositional strike-slip faults in the Longwangmiao Formation in Gaomo area of Central Sichuan Basin (a) and modern San Andreas Fault Zone and related topography(b)
采集于2011年的面積為2500ikm2的三維地震數(shù)據(jù),面元20im×20im,在龍王廟組分布范圍內(nèi)的頻率帶寬為10~60iHz,主頻30iHz,地震垂向分辨率約為50im。由于龍王廟組白云巖上下均被波阻抗較低的碎屑巖所圍限,因此地震資料的信噪比相對(duì)較高。這種得天獨(dú)厚的巖性組合,有效提高了地震資料的探測(cè)能力。
研究區(qū)三維地震測(cè)網(wǎng)范圍內(nèi),鉆穿龍王廟組的共計(jì)35口井(圖 6-a),每口井都有AC、RHOB等常規(guī)測(cè)井曲線數(shù)據(jù),為精細(xì)井震標(biāo)定奠定資料基礎(chǔ)。與此同時(shí),每口井在龍王廟組都有測(cè)井解釋孔隙度資料,為計(jì)算龍王廟組孔隙度并將其與地震資料做關(guān)聯(lián)提供了條件。
針對(duì)龍王廟組白云巖地層,在古地貌恢復(fù)的基礎(chǔ)上,為了實(shí)現(xiàn)利用地震沉積學(xué)方法定性和定量計(jì)算儲(chǔ)集層平面分布的目標(biāo),制定了如下研究流程。首先,在-90°地震數(shù)據(jù)體上精細(xì)解釋了龍王廟組,利用層拉平技術(shù),恢復(fù)了沉積古地貌。第2步,利用分頻技術(shù),將相位旋轉(zhuǎn)之后的地震數(shù)據(jù)體分成高、中、低3個(gè)頻段分頻體,在分頻體上提取了對(duì)儲(chǔ)集層厚度敏感的平均振幅屬性,得到薄、中、厚儲(chǔ)集層的平面振幅屬性分布圖。第3步,利用RGB三色融合技術(shù),將3張分頻振幅屬性分布圖融合成定性儲(chǔ)集層厚度預(yù)測(cè)圖。第4步,將-90°相位地震數(shù)據(jù)體分為7個(gè)分頻體,分別提取7個(gè)分頻體的平均振幅屬性。第5步,利用PCA方法分析7個(gè)分頻振幅屬性,將前2個(gè)主因子與井點(diǎn)處解釋的儲(chǔ)集層厚度擬合,得到擬合關(guān)系式。第6步,將前2個(gè)主因子代入擬合關(guān)系式,計(jì)算獲得儲(chǔ)集層厚度定量預(yù)測(cè)結(jié)果。
2.3.1 古地貌及其與走滑斷裂的關(guān)系
結(jié)合研究區(qū)的構(gòu)造演化史和地層發(fā)育情況,通過層拉平的方式,定性恢復(fù)了龍王廟組沉積古地貌(圖 7-a)。研究區(qū)整體西北高、東南低,表現(xiàn)為寬緩的淺水陸棚特征。磨溪構(gòu)造總體位于高地之上斜坡部位,而高石梯構(gòu)造則是位于上斜坡至下斜坡部位。研究區(qū)被幾條大斷裂和若干伴生小斷裂所切割,形成類似地塹—地壘的地貌特征。其中,近南北向的大斷裂(A-A′)是繼承自震旦紀(jì)的一條深大斷裂,是德陽—安岳裂陷槽東緣斷裂(杜金虎等,2016),不具走滑性質(zhì)。
北東—南西向的另一條大斷裂(B-B′)則是1條寒武紀(jì)持續(xù)活動(dòng)的同沉積走滑斷裂。該斷裂平面呈帚狀展布, 由北東向南西方向散開, 具右旋特征。走滑斷裂體系主斷層(B-B′)及伴生的一系列次級(jí)走滑斷裂, 對(duì)古地貌的調(diào)整作用非常明顯。實(shí)際上, 同沉積走滑斷裂控制地形地貌的現(xiàn)代實(shí)例也很常見, 現(xiàn)今美國圣安德烈斯大斷裂就是典型代表(圖 7-b)。在右旋走滑斷裂體系的控制下, 洛杉磯海岸及淺海陸棚之上的次盆(地塹)和水下隆起(地壘)的分布, 與走滑斷裂具有很好的對(duì)應(yīng)關(guān)系。
圖 8 川中高磨地區(qū)利用測(cè)井?dāng)?shù)據(jù)內(nèi)插儲(chǔ)集層厚度(a)與三色融合地震屬性得到的定性儲(chǔ)集層厚度(b)Fig.8 Reservoir thickness derived by inserting well data(a)and qualitative reservoir thickness derived by RGB blending seismic attributes(b)in Gaomo area of Central Sichuan Basin
2.3.2 白云巖儲(chǔ)集層定性預(yù)測(cè)
已經(jīng)有學(xué)者指出,地震振幅屬性能夠反映儲(chǔ)集層(李操,2014),然而直接應(yīng)用振幅預(yù)測(cè)儲(chǔ)集層物性往往會(huì)受到復(fù)雜巖性的干擾。但研究目的層龍王廟組因?yàn)槠鋷r性單一(白云巖),大大降低了巖性對(duì)于振幅屬性預(yù)測(cè)儲(chǔ)集層的干擾。
根據(jù)實(shí)鉆鉆井?dāng)?shù)據(jù)和孔隙度測(cè)井解釋數(shù)據(jù),統(tǒng)計(jì)出井點(diǎn)處的龍王廟組地層厚度和儲(chǔ)集層(孔隙度大于4%)累計(jì)厚度,進(jìn)而利用內(nèi)插法繪制出基于測(cè)井解釋數(shù)據(jù)的儲(chǔ)集層平面分布圖(圖 8-a)。從圖8中可以看出,研究區(qū)內(nèi)龍王廟組白云巖儲(chǔ)集層厚度介于0~60im之間,平均厚度約30im。儲(chǔ)集層厚度與古地貌(圖 7-a)具有較好的對(duì)應(yīng)關(guān)系,厚儲(chǔ)集層基本位于古地貌高部位周邊的上斜坡位置。同時(shí),儲(chǔ)集層平面分布受控于走滑斷裂體系,而繼承自震旦紀(jì)的近南北向大斷裂對(duì)于儲(chǔ)集層厚度的影響不大。
根據(jù)井點(diǎn)解釋儲(chǔ)集層厚度的分布特點(diǎn),選取70im、45im和33im作為厚、中、薄儲(chǔ)集層的代表厚度值,計(jì)算了其所對(duì)應(yīng)的低、中、高調(diào)諧頻率值分別為18iHz、28iHz和38iHz。在這3個(gè)頻段分頻數(shù)據(jù)體上,提取了平均振幅屬性。然后利用RGB三色融合技術(shù)(紅色代表低頻厚儲(chǔ)集層、綠色代表中頻中等厚度儲(chǔ)集層、藍(lán)色代表高頻薄儲(chǔ)集層),得到了研究區(qū)內(nèi)龍王廟組白云巖儲(chǔ)集層厚度定性分布(圖 8-b)。
與井?dāng)?shù)據(jù)內(nèi)插得到的儲(chǔ)集層厚度相比,三色融合儲(chǔ)集層厚度分布圖反映了大致相同的趨勢(shì)(西北厚、東南薄),古地貌和斷層對(duì)儲(chǔ)集層厚度分布的控制作用同樣明顯。同時(shí),融合儲(chǔ)集層厚度圖上反映了更多細(xì)節(jié),特別是反映了少井或者無井地區(qū)的儲(chǔ)集層分布情況,體現(xiàn)了地震資料預(yù)測(cè)儲(chǔ)集層的優(yōu)勢(shì)。對(duì)比磨溪構(gòu)造和高石梯構(gòu)造儲(chǔ)集層分布特征發(fā)現(xiàn): 磨溪構(gòu)造發(fā)育的儲(chǔ)集層厚度更大、面積更廣,表現(xiàn)為塊狀連片厚儲(chǔ)集層特征;高石梯構(gòu)造發(fā)育的儲(chǔ)集層整體厚度較小、面積局限,平面上表現(xiàn)為環(huán)帶狀分布的特征。這可能是由2個(gè)構(gòu)造圈閉所處的古地貌存在差異所導(dǎo)致: 磨溪構(gòu)造緊鄰古地貌高地的上斜坡部位,而高石梯構(gòu)造則地處遠(yuǎn)離古地貌高部位的中下斜坡位置。
圖 9 川中高磨地區(qū)利用振幅屬性主因子分析結(jié)果(a和b)擬合計(jì)算的儲(chǔ)集層厚度定量分布(c)Fig.9 Quantitatively distribution of reservoir(c)calculated by principle components derived from seismic attributes(a and b)in Gaomo area of Central Sichuan Basin
2.3.3 白云巖儲(chǔ)集層定量預(yù)測(cè)
為了更好地指導(dǎo)龍王廟組氣藏的勘探與開發(fā),除了定性刻畫不同厚度規(guī)模儲(chǔ)集層的平面分布之外,還有必要開展儲(chǔ)集層厚度定量預(yù)測(cè),因此,引入了主因子分析技術(shù)。
如前所述, 將-90°數(shù)據(jù)體由低到高分成7個(gè)分頻體之后, 分別提取了目的層龍王廟組的平均振幅屬性, 繼而進(jìn)行了主因子分析。在得到的7個(gè)主因子中, 前2個(gè)主因子含有絕大部分的原有信息, 可以作為自變量進(jìn)行擬合運(yùn)算。實(shí)際上, 第1個(gè)主因子已經(jīng)能基本反映厚儲(chǔ)集層平面分布情況(圖 9-a)。
求取前2個(gè)主因子在井點(diǎn)處的數(shù)值,將其與井點(diǎn)測(cè)井解釋儲(chǔ)集層累計(jì)厚度作擬合,下列關(guān)系式給出了比較理想的相關(guān)關(guān)系(R=0.78)。利用該關(guān)系式對(duì)前2個(gè)主因子進(jìn)行運(yùn)算,得到了研究區(qū)龍王廟組白云巖儲(chǔ)集層的厚度平面分布圖(圖 9-c)。
其中,z是井點(diǎn)處儲(chǔ)集層厚度,x是第1個(gè)主因子在井點(diǎn)處的數(shù)值,y是第2個(gè)主因子在井點(diǎn)處的數(shù)值, 常數(shù)系數(shù)a至g分別是: a=26.40、 b=-64.79、 c=36.75、 d=40.53、 e=-1.58、 f=3.51、 g=4.82。
從圖9中可以看出,厚度大于30im的中厚儲(chǔ)集層的分布情況與定性儲(chǔ)集層預(yù)測(cè)結(jié)果類似,大部分位于磨溪構(gòu)造,少部分處在高石梯構(gòu)造。厚度大于50im的厚儲(chǔ)集層大部分發(fā)育在磨溪構(gòu)造。雖然與定性儲(chǔ)集層分布圖反映的趨勢(shì)大體一致,但是量化預(yù)測(cè)結(jié)果揭示出更多的細(xì)節(jié)差異。例如,在定性圖上,磨溪地區(qū)厚儲(chǔ)集層呈連片分布,平面均質(zhì)性較好。但是在定量圖上,磨溪構(gòu)造的這些厚儲(chǔ)集層并非連片分布,而是存在較強(qiáng)的平面非均質(zhì)性。這一認(rèn)識(shí),為后續(xù)的滾動(dòng)勘探和氣田開發(fā)提供了更為精細(xì)的重要參考。
圖 10 四川盆地三疊系嘉陵江組巖相古地理(a)與層序地層劃分(b)Fig.10 Lithofacies palaeogeography(a)and sequence stratigraphy(b)of the Triassic Jialingjiang Formation in Sichuan Basin
四川盆地面積約18×104ikm2,是在揚(yáng)子地臺(tái)基礎(chǔ)上發(fā)育起來的大型沉積盆地,其周圍被龍門山、米倉山、大巴山、七曜山、婁山和大涼山所環(huán)繞(黃汲清,1954)。在早三疊世嘉陵江組沉積時(shí)期,四川盆地巖相古地理包括深海、古陸、局限臺(tái)地和開闊臺(tái)地(王宓君等,1989)。海水分別從西北、東北和東南3個(gè)方向侵入四川盆地(田在藝和張慶春,1997),盆地從海相環(huán)境向陸相環(huán)境過渡(圖 10-a)。此時(shí),研究區(qū)位于低緯度區(qū),氣候從半干旱向干旱轉(zhuǎn)變(Weidlich,2007;徐國盛等,2011)。從早三疊世開始,研究區(qū)周邊地區(qū)開始逐漸隆升,導(dǎo)致了相對(duì)局限的沉積環(huán)境。盆地西南部地區(qū),季節(jié)性活躍的河流將碎屑物質(zhì)從康滇古陸搬運(yùn)進(jìn)入盆地發(fā)生沉積(吳應(yīng)林等,1989)。
嘉陵江組(T1j)自下而上可分為嘉1段至嘉5段(T1j1-T1j5)共5段,其中,嘉1段和嘉3段由灰?guī)r構(gòu)成;嘉2段包括膏巖、白云巖、灰?guī)r和碎屑巖;嘉4段和嘉5段則發(fā)育膏巖、白云巖、灰?guī)r和鹽巖(圖 10-b)。嘉陵江組地層可劃分為3個(gè)三級(jí)層序,自下而上分別為: 嘉1段至嘉2段(層序1)、嘉3段至嘉4段(層序2)、嘉5段(層序3)。其中,研究目的層嘉2下亞段(T1j2L)的厚度在平面上較為穩(wěn)定,地層厚度相對(duì)較薄(介于10~70im),平均厚約50im,在-90°相位旋轉(zhuǎn)后的地震剖面上對(duì)應(yīng)1個(gè)波谷(Huetal., 2017)。
T1j2L中的碎屑巖呈現(xiàn)灰綠色, 由泥巖和粉砂巖隔夾層組成, 推測(cè)其形成于陸相沉積環(huán)境(Zhaoetal., 2018)。 泥質(zhì)灰?guī)r呈深灰色薄層狀, 推測(cè)其形成于海陸過渡環(huán)境。 淺灰色白云巖孔隙十分發(fā)育, 表現(xiàn)為針孔狀特征, 根據(jù)薄片測(cè)定數(shù)據(jù), 其孔隙度最高可達(dá)20%, 平均為8%左右。 膏質(zhì)云巖和云質(zhì)膏巖也很常見, 發(fā)育典型的揉皺變形層理。 隨著蒸發(fā)強(qiáng)度的增加, 出現(xiàn)了膏巖(硬石膏)沉積, 為白色或者淺灰色塊狀特征。 這些巖性垂向間互頻繁、 橫向相變迅速, 為混積巖相研究提供了良好的解剖區(qū)。
圖 11 四川盆地高磨地區(qū)位置及鉆穿目的層的井分布Fig.11 Location of Gaomo area in Sichuan Basin and distribution of wells drilling through the target layer
為了確定混積環(huán)境下的巖相分布,制定了針對(duì)性的技術(shù)方案。綜合利用巖心分析數(shù)據(jù)、測(cè)井解釋數(shù)據(jù)和地震屬性分析數(shù)據(jù)等,首先找到不同巖性的各項(xiàng)地球物理參數(shù)特征,據(jù)此計(jì)算不同巖性含量在井點(diǎn)處的數(shù)值并利用地震資料預(yù)測(cè)其平面分布,進(jìn)而根據(jù)平面巖性組合確定巖相分布。
3.3.1 巖心分析與測(cè)井?dāng)?shù)據(jù)分析結(jié)果
利用巖心分析(礦物成分、薄片鑒定等),可以根據(jù)礦物組分和物性情況識(shí)別出5種不同的巖性/巖相,分別是膏巖、碎屑巖、致密白云巖、灰?guī)r和多孔白云巖。結(jié)合相應(yīng)深度的測(cè)井曲線參數(shù),分析了不同巖性的曲線參數(shù)分布特征(圖 12)。從圖中可以看出,碎屑巖的GR值最高(圖 12-a),據(jù)此可將GR值高于45 API的曲線對(duì)應(yīng)的巖性解釋為碎屑巖。從不同巖性的密度分布范圍來看,膏巖的密度最高(圖 12-b),將密度值大于2.89ig/cm3的曲線段解釋為膏巖。識(shí)別出碎屑巖和膏巖之后,剩下的巖性包括致密白云巖、灰?guī)r和多孔白云巖,很難通過GR和RHOB將三者區(qū)分開。同時(shí),三者的AC分布范圍也很相近(圖 12-c),因此將其歸為碳酸鹽巖這一大項(xiàng)。
圖 12 四川盆地高磨地區(qū)基于鉆測(cè)井資料的不同巖性關(guān)鍵參數(shù)Fig.12 Key factors for different lithologies in based on well data in Gaomo area of Sichuan Basin
利用RHOB和AC計(jì)算了不同巖性波阻抗(AI)的分布情況(圖 12-d)??梢钥闯?,不同巖性的AI數(shù)值分布范圍有搭接,但大部分巖性分布的主峰是錯(cuò)開的。其中,致密白云巖的AI值最高,其次是膏巖,碎屑巖的AI值最低,多孔白云巖和灰?guī)r的AI主峰基本重疊。這種AI的分布特征,意味著利用地震數(shù)據(jù)幾乎不可能將灰?guī)r和多孔白云巖完全分開。同時(shí),在地震數(shù)據(jù)上應(yīng)該有不同的響應(yīng)(因?yàn)槿叩腁I主峰有差異),盡管如此,地震數(shù)據(jù)很難明確地將混積環(huán)境下的這幾種巖性一一甄別(因?yàn)锳I分布有重疊)。
圖 13 四川盆地高磨地區(qū)種巖性地震屬性主因子分析結(jié)果特征圖Fig.13 PCA results of seismic attributes generated for three lithologies of in Gaomo area of Sichuan Basin
3.3.2 主因子分析結(jié)果
正是由于AI的上述特征,導(dǎo)致利用單一地震屬性難以區(qū)分不同的巖性。為了解決這一問題,本文引入主因子分析技術(shù),對(duì)地震屬性進(jìn)行重構(gòu)分析,將其中的巖性信息優(yōu)化提取形成線性不相關(guān)的主因子。具體實(shí)施過程中,分別針對(duì)碎屑巖(Sc)、膏巖(An)和碳酸鹽巖(Ca),根據(jù)井點(diǎn)處巖性百分含量與不同地震屬性之間的相關(guān)性,分別優(yōu)選了3組(對(duì)應(yīng)3種巖性)相關(guān)性較好的地震屬性,然后對(duì)優(yōu)選出來地震屬性分別進(jìn)行主因子分析。其中,針對(duì)碎屑巖選擇的地震屬性主要為振幅、弧長和頻率;針對(duì)碳酸鹽巖選擇的地震屬性主要為振幅、頻率和能量;針對(duì)膏巖選擇的地震屬性主要為頻率、振幅、能量和有效帶寬。
圖 14 四川盆地高磨地區(qū)不同巖性(a、b、c)百分含量平面分布情況及RGB巖相特征(d)Fig.14 Distribution maps of lithology content(a、b、c)and the RGB lithofacies map(d)in Gaomo area of Sichuan Basin
主因子分析結(jié)果(圖 13)顯示,對(duì)于3種巖性而言,主因子所含有的原有地震屬性逐漸減少,而前2個(gè)主因子涵蓋了原有屬性中絕大部分信息(75%以上),因此可以利用這2個(gè)主因子代表原有地震屬性進(jìn)行擬合運(yùn)算。這種做法不但減少了參與運(yùn)算的自變量個(gè)數(shù)、降低了擬合難度,而且最大限度地保留了原有地震屬性中蘊(yùn)含的巖性信息,在保證精度的前提下,提高了擬合運(yùn)算效率。
3.3.3 巖性含量擬合計(jì)算結(jié)果
3.3.4 混積巖相分布特征
分別獲取3種巖性含量平面分布之后,將這3種巖性作為3個(gè)端元,利用RGB三色融合技術(shù),將3種巖性含量轉(zhuǎn)化成相應(yīng)的巖相,例如膏巖含量高指示膏巖巖相、碎屑巖含量高指示碎屑巖巖相、碳酸鹽巖含量高指示碳酸鹽巖巖相,據(jù)此原則可以將融合圖解釋為巖相圖(圖 14-d)。其中,紅色代表膏巖、藍(lán)色代表碎屑巖、綠色代表碳酸鹽巖。從巖相圖上可以看出,巖相分帶性比較明顯,共分為5個(gè)區(qū)域(即由紅色虛線分割開的A至E)。其中,區(qū)域A以碎屑巖為主,含有痕量的碳酸鹽巖。區(qū)域B是碎屑巖和碳酸鹽巖共存的地區(qū),表現(xiàn)為海陸過渡沉積特征。區(qū)域C以碳酸鹽巖為主,反映了正常海沉積。區(qū)域D是研究區(qū)內(nèi)面積最大的,其巖性組合為膏巖和碳酸鹽巖,北部局部地區(qū)還含有痕量的碎屑巖,這種巖性組合是一種典型的局限環(huán)境。在位于東北角的區(qū)域E,巖性組合相對(duì)簡單,以碳酸鹽巖為主,膏巖含量明顯降低至痕量,表現(xiàn)出由局限環(huán)境向半局限或開闊環(huán)境的過渡。
圖 15 澳大利亞圣誕島上發(fā)育的碎屑巖—碳酸鹽巖—蒸發(fā)巖混積體系的現(xiàn)代實(shí)例Fig.15 Modern analogy for mixed siliciclastic-carbonate-evaporite systems on the Christmas Island in Australia
除了相位旋轉(zhuǎn)、地層切片和分頻3種傳統(tǒng)技術(shù)之外,地震沉積學(xué)還包括主因子分析和RGB三色融合2項(xiàng)新技術(shù)。其中,前3種傳統(tǒng)技術(shù)已經(jīng)被證實(shí)可以在沉積相研究中發(fā)揮巨大作用(董春梅等,2006;彭才等,2011)。后2種新引入的技術(shù)旨在利用井震資料,定性和定量預(yù)測(cè)儲(chǔ)集層分布,其應(yīng)用尚處于起步階段(姜秀娣等,2013;Chopra and Marfurt,2014;Zeng,2018)。通過3個(gè)研究實(shí)例的分析,可以看出3種沉積體系共同的特征是地層薄,導(dǎo)致常規(guī)地震資料解釋困難。同時(shí),由于沉積環(huán)境、巖性組合和后期成巖演化過程的差異,需要采用不同的技術(shù)組合以實(shí)現(xiàn)沉積相和儲(chǔ)集層預(yù)測(cè)之目的。筆者希望通過討論不同沉積體系的具體差異,探討地震沉積學(xué)針對(duì)這3種體系選擇技術(shù)組合時(shí)的依據(jù)和應(yīng)用效果,以期推進(jìn)地震沉積學(xué)在量化研究儲(chǔ)集層分布方面的發(fā)展,為從事類似研究的學(xué)者提供借鑒。
針對(duì)碎屑巖沉積體系而言,因其地形地貌特征明顯和后期成巖作用序列簡單(劉林玉等,2002),可以重點(diǎn)采納地震地貌學(xué)研究思路,采取“相位旋轉(zhuǎn)、地層切片和RGB三色融合”技術(shù)組合,定性恢復(fù)沉積相。研究結(jié)果揭示了瑪湖凹陷克拉瑪依組沖積扇—河流的微地貌和沉積相特征,為準(zhǔn)確識(shí)別巖性提供有效依據(jù)。因波阻抗相近而導(dǎo)致常規(guī)地球物理方法難以準(zhǔn)確識(shí)別礫巖和砂巖(吳兆徽等,2016),但是在微地貌的約束下,兩者變得涇渭分明,即礫巖通常發(fā)育在沖積扇體系而砂巖通常發(fā)育在河流體系。確定了沉積相和巖性平面分布之后,就可以合理解釋有利儲(chǔ)集層和高產(chǎn)井的分布規(guī)律,為下一步的油氣勘探開發(fā)提供指導(dǎo)。
碳酸鹽巖沉積體系中的灰?guī)r和白云巖雖然存在波阻抗差異, 但是白云巖往往發(fā)育孔隙, 同時(shí), 后期演化情況復(fù)雜, 造成2種巖性的波阻抗有時(shí)相差無幾(劉玲等,2020)。 因此, 利用地震資料預(yù)測(cè)儲(chǔ)集層過程中, 巖性差異直接造成物性預(yù)測(cè)的多解性。 要想準(zhǔn)確預(yù)測(cè)儲(chǔ)集層分布, 首先要盡量減少巖性差異帶來的影響。選擇的川中高磨地區(qū)龍王廟組具有得天獨(dú)厚的優(yōu)勢(shì): 首先, 大量鉆井揭示, 龍王廟組沉積地層已經(jīng)完全白云石化(張建勇等, 2015), 不發(fā)育灰?guī)r, 因此沒有巖性帶來的干擾; 其次, 龍王廟組上、 下被碎屑巖所圍限, 頂?shù)捉缑嫣幍牟ㄗ杩共町惷黠@。在這2個(gè)有利條件的保證下, 加之大量鉆井的標(biāo)定, 重點(diǎn)借鑒地震巖性學(xué)研究思路, 采用“相位旋轉(zhuǎn)、 分頻地震屬性RGB融合、 主因子分析”的技術(shù)組合, 高質(zhì)量地完成了沉積古地貌的定性恢復(fù)和白云巖儲(chǔ)集層的定量預(yù)測(cè)。
混積巖沉積體系通常指碎屑巖和碳酸鹽巖的混合沉積(Guertinetal., 2000;Brandanoetal., 2010),而川中高磨地區(qū)嘉陵江組嘉二段中碎屑巖、碳酸鹽巖和蒸發(fā)巖3種巖性的混積則更為復(fù)雜。在這種垂向和橫向巖相快速變化、非均質(zhì)性(巖性非均質(zhì)性和物性非均質(zhì)性)極強(qiáng)的沉積體系中,有效預(yù)測(cè)巖相是預(yù)測(cè)潛在儲(chǔ)集層發(fā)育區(qū)的前提。預(yù)測(cè)研究區(qū)內(nèi)巖相分布的1條有效途徑是充分利用鉆井信息標(biāo)定地震資料,為了實(shí)現(xiàn)這一目標(biāo),采用了“地震屬性+主因子分析”的技術(shù)組合。這種技術(shù)組合有效地將單井巖相定量識(shí)別結(jié)果外推到地震數(shù)據(jù),得到了巖相平面分布,識(shí)別出有利于儲(chǔ)集層發(fā)育的巖相,指出了有利的油氣勘探區(qū)帶。盡管地震沉積學(xué)已經(jīng)具備了定性恢復(fù)沉積相和定量預(yù)測(cè)儲(chǔ)集層(特別是常規(guī)地震方法無能為力的薄地層中的儲(chǔ)集層)的能力,但還有許多可以提升和改進(jìn)之處,例如分頻時(shí)的頻率是依據(jù)地質(zhì)家對(duì)研究區(qū)地層厚度的經(jīng)驗(yàn)評(píng)估值計(jì)算得到的,必然存在一定偏差。如果能研發(fā)出自動(dòng)掃描地層厚度分布的方法,即可據(jù)此計(jì)算更為精確的分頻頻率,效果必然更佳。另外,不同地區(qū)的沉積背景和研究目的有所差異,因此筆者提供的3種不同的應(yīng)用場景僅供相關(guān)的研究者參考,實(shí)際應(yīng)用時(shí)應(yīng)具體問題具體分析,針對(duì)性地選擇技術(shù)組合和研究流程。筆者期望此文能夠拋磚引玉,為地震沉積學(xué)更廣泛地應(yīng)用提供助力。
油氣勘探實(shí)踐助推地震采集處理技術(shù)不斷提高,計(jì)算機(jī)運(yùn)算能力顯著增強(qiáng),地震沉積學(xué)迎來良好發(fā)展契機(jī)?!?(地震地貌學(xué))+1(地震巖性學(xué))”的理論體系和“3(相位旋轉(zhuǎn)、地層切片、頻譜分解)+2(主因子分析、三色融合)”的技術(shù)系列,大大拓展了地震沉積學(xué)的應(yīng)用范疇。研究表明,地震沉積學(xué)方法和技術(shù)可以在不同環(huán)境下,有效恢復(fù)沉積相,定量計(jì)算儲(chǔ)集層厚度。
1)準(zhǔn)噶爾盆地瑪湖凹陷三疊系克拉瑪依組被作為砂礫巖沉積的典型代表,利用地震地貌學(xué)原理和地層切片融合技術(shù),首次在克拉瑪依組上段S1砂層組發(fā)現(xiàn)了曲流河沉積,該河流體系是在沖積扇體系沉積之后形成的。兩者都在溝谷中沉積,但是扇先河后。利用主因子分析技術(shù)、結(jié)合井震標(biāo)定,發(fā)現(xiàn)研究區(qū)儲(chǔ)集層厚度小于10im,主要集中在河流的點(diǎn)壩沉積和沖積扇的扇緣沉積中,河流相關(guān)儲(chǔ)集層因其較好的孔滲條件而油氣產(chǎn)能更高。
2)川中高磨地區(qū)寒武系龍王廟組為單一的白云巖地層,受同沉積走滑斷裂的調(diào)整,形成了大型淺水緩坡背景下的“壘-塹間互”的構(gòu)造格局。利用地震沉積學(xué)的分頻融合和主因子分析技術(shù),定性定量恢復(fù)了儲(chǔ)集層平面分布。發(fā)現(xiàn)磨溪構(gòu)造是龍王廟組儲(chǔ)集層分布主體區(qū)域,大面積分布在古地貌高地周邊的斜坡區(qū),整體呈塊狀,局部存在非均質(zhì)性;高石梯地區(qū)的儲(chǔ)集層呈環(huán)帶狀局限分布,儲(chǔ)集層厚度基本小于40im。
3)川中嘉陵江組嘉二段發(fā)育碳酸鹽巖—碎屑巖—蒸發(fā)巖混積環(huán)境,利用伽馬識(shí)別碎屑巖,利用密度識(shí)別膏巖。主因子分析技術(shù)突出了地震屬性中的巖性差異,據(jù)此預(yù)測(cè)了3種巖性含量的平面分布。三色融合技術(shù)生成的巖相圖顯示,研究區(qū)從西南向東北,逐漸由陸相經(jīng)海陸過渡相和蒸發(fā)潟湖相,最終變?yōu)殚_闊海相。
致謝感謝3位匿名審稿專家對(duì)文章提出的建設(shè)性建議和具體修改意見!