郭 雯 劉永濤 趙俊峰 王秀珍 趙紅格 丁富峰
(①中國(guó)石油勘探開發(fā)研究院,北京100083;②西北大學(xué)大陸動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,陜西西安710069;③西北大學(xué)地質(zhì)學(xué)系,陜西西安710069;④東方地球物理公司研究院長(zhǎng)慶分院,陜西西安710021;⑤中國(guó)石油長(zhǎng)慶油田公司勘探事業(yè)部,陜西西安710018)
地震拉平技術(shù)是一種展平化地震數(shù)據(jù)的分析方法,主要依據(jù)某種原理或算法,對(duì)原始地震數(shù)據(jù)重新移位或采樣,生成全新的數(shù)據(jù)體,以提高地震數(shù)據(jù)處理的質(zhì)量或便于地質(zhì)分析。目前地震拉平技術(shù)主要用于:①在地震資料處理方面,通過改變常規(guī)數(shù)據(jù)疊加方法、Wheeler變換及傾斜地層同相軸校正,提高地震資料信噪比和成像質(zhì)量[1-3];②在地震資料解釋方面,通過數(shù)據(jù)體拉平[4-5]、拾取層位拉平[6-7]或平衡剖面[8]提高地震解釋效率,挖掘資料中潛在的地質(zhì)信息,開辟一條新的解釋途徑。地震層拉平是地震拉平常用的解釋技術(shù)之一,是指將地震數(shù)據(jù)體沿某一拾取層位拉平,在全新的拉平數(shù)據(jù)體上進(jìn)行解釋。自20世紀(jì)80年代層拉平技術(shù)引入解釋工作站以來(lái),地震層拉平技術(shù)廣泛用于儲(chǔ)層預(yù)測(cè)[9-11]、古構(gòu)造恢復(fù)[12-14]、古 地 貌 恢 復(fù)[15-16]、地 質(zhì) 演 化 分 析[17-19]及油藏動(dòng)態(tài)分析[20-21]等油氣勘探領(lǐng)域,呈現(xiàn)出巨大的應(yīng)用潛力。
隨著地震層拉平技術(shù)應(yīng)用領(lǐng)域的拓展,也暴露出一些亟待解決的問題。如:地震層拉平技術(shù)的應(yīng)用條件;地震層拉平在古地表恢復(fù)過程中對(duì)地震解釋精度的影響;地震層拉平與地震切片、地層厚度法、平衡剖面技術(shù)的區(qū)別和聯(lián)系。這些問題的解答,有助于地震層拉平技術(shù)在油氣勘探領(lǐng)域的深化應(yīng)用。本文從技術(shù)原理出發(fā),首先將地震層拉平技術(shù)與其他古地表恢復(fù)技術(shù)進(jìn)行對(duì)比、分析,明確產(chǎn)生恢復(fù)誤差的原因,然后通過增加地質(zhì)約束條件或修正原有技術(shù)方案,降低誤差對(duì)恢復(fù)結(jié)果的影響,從而推動(dòng)地震層拉平技術(shù)在儲(chǔ)層預(yù)測(cè)、古地貌恢復(fù)、古構(gòu)造恢復(fù)、構(gòu)造演化分析等4個(gè)方面的深化應(yīng)用。
從技術(shù)角度講,地震層拉平技術(shù)是以某一拾取層位為固定基準(zhǔn)面,對(duì)地震數(shù)據(jù)垂直移位,生成全新的解釋數(shù)據(jù)。目前,解釋工作站上都有相應(yīng)的操作模塊,操作非常方便。實(shí)際上地震層拉平的過程也是一個(gè)對(duì)數(shù)據(jù)體重新運(yùn)算的過程。設(shè)原始數(shù)據(jù)為輸入函數(shù)f(x,y,t),運(yùn)算過程即利用某種映射算法(函數(shù))t(x,y,T)對(duì)f(x,y,t)重新排序,生成的新數(shù)據(jù)體為

式中:x 為原始時(shí)間域數(shù)據(jù)體水平方向的坐標(biāo);y為原始時(shí)間域數(shù)據(jù)體垂直方向的坐標(biāo);t為原始時(shí)間域數(shù)據(jù)體的雙程旅行時(shí);T 為拉平時(shí)間域數(shù)據(jù)體雙程旅行時(shí)。
由于層拉平只是對(duì)數(shù)據(jù)體的樣點(diǎn)垂直移位,不考慮樣點(diǎn)的真實(shí)歸位,所以
t(x,y,T)=T-t0
式中t0為層拉平選取的基準(zhǔn)面固定值。則式(1)簡(jiǎn)化為

在三維空間內(nèi),要使每個(gè)樣點(diǎn)真實(shí)歸位,必須考慮恢復(fù)地層的真實(shí)厚度,而地震層拉平只是對(duì)數(shù)據(jù)體的樣點(diǎn)垂直移位,勢(shì)必造成一定恢復(fù)誤差。
在地質(zhì)意義上,地震層拉平是一種古地表重建技術(shù),它是將地震數(shù)據(jù)體沿某一特定地質(zhì)界面對(duì)應(yīng)的地震反射層位拉平,恢復(fù)地史時(shí)期的古地貌,進(jìn)而考查下伏地層的古構(gòu)造、古沉積特征及地質(zhì)演化過程,從而在三維空間內(nèi)進(jìn)行儲(chǔ)層預(yù)測(cè)、古地貌恢復(fù)及含油氣預(yù)測(cè)等。需要注意的是,層拉平技術(shù)需要滿足兩個(gè)前提條件:①選取的待拉平地震反射層位的原始產(chǎn)狀必須是水平的或近似水平的,且在地震剖面上可連續(xù)追蹤,否則恢復(fù)的古地貌沒有明確的地質(zhì)涵義;②后期構(gòu)造變形不甚強(qiáng)烈,否則地層變形已發(fā)生構(gòu)造重組,再進(jìn)行地質(zhì)恢復(fù)必然出現(xiàn)假象。當(dāng)滿足上述條件時(shí),還要進(jìn)一步考慮恢復(fù)精度問題。地質(zhì)演化過程一般造成地層長(zhǎng)度(橫向)伸縮和厚度(縱向)變化,而在三維地震數(shù)據(jù)中,利用層拉平技術(shù)恢復(fù)古地貌并不會(huì)改變剖面的長(zhǎng)度和地層的厚度,因此客觀事實(shí)與技術(shù)局限之間的矛盾是造成古地貌恢復(fù)誤差的根本原因。文中以斷陷湖盆的古地貌恢復(fù)為例進(jìn)行說(shuō)明。
圖1展示了層拉平技術(shù)恢復(fù)古地貌原理。由圖可見:①由于構(gòu)造應(yīng)力場(chǎng)的變化,湖盆發(fā)育在基底褶皺之上,隨后進(jìn)入構(gòu)造穩(wěn)定期,斷陷湖盆被“填平補(bǔ)齊”(圖1a);②在后期演化過程中,構(gòu)造應(yīng)力場(chǎng)再次發(fā)生變化,由于受到水平擠壓和垂向壓實(shí)作用,近似水平的層狀地層發(fā)生變形、彎曲,地層厚度發(fā)生變化,而且湖盆基底也被斷層復(fù)雜化,發(fā)生構(gòu)造重組,基底褶皺高點(diǎn)位置已由原來(lái)的a 點(diǎn)遷移到b 點(diǎn)(圖1b);③選取湖盆頂部近似代表古水平面的地震反射層位作為拉平層位,利用地震層拉平技術(shù)恢復(fù)古地貌,得到下伏湖盆充填地層的形態(tài)及基底構(gòu)造樣式(圖1c),大致恢復(fù)了湖盆充填地層的形態(tài),但湖盆的長(zhǎng)度、充填地層的厚度與地史時(shí)期的原始狀態(tài)(圖1a)還存在較大差異,而且基底褶皺已發(fā)生構(gòu)造重組,高點(diǎn)c已無(wú)法恢復(fù)到位置a。

圖1 層拉平技術(shù)恢復(fù)古地貌原理
利用地震層拉平技術(shù)恢復(fù)古地貌,造成湖盆充填地層的恢復(fù)厚度與原始沉積厚度差異的原因?yàn)椋阂皇怯捎趯有蛏稀⑾碌貙觾A角的差異及地層沉積演化過程中的脫水、壓實(shí)作用導(dǎo)致地層厚度變化,這是主因;二是由于構(gòu)造演化過程中發(fā)生的地層橫向伸縮作用,而地震層拉平技術(shù)并不改變地層的長(zhǎng)度,因此在構(gòu)造恢復(fù)中也造成一定誤差。針對(duì)地震層拉平技術(shù)恢復(fù)古地貌產(chǎn)生的誤差,結(jié)合具體實(shí)例,分析這些誤差對(duì)恢復(fù)結(jié)果的影響程度,并通過增加地質(zhì)條件約束或改進(jìn)技術(shù)方案,最大程度地降低誤差對(duì)恢復(fù)結(jié)果的影響,使地震層拉平技術(shù)能更好地應(yīng)用于地震資料解釋的各個(gè)環(huán)節(jié)。
地震切片分為水平(或時(shí)間)切片、沿層切片和地層切片,其中沿層切片又分為沿頂切片和沿底切片。沿層切片定義為沿某一層序界面進(jìn)行層拉平,然后向上、向下水平切割拉平數(shù)據(jù)體而形成的切片。地層切片定義為在層序頂、底界面間按等比例進(jìn)行插值而形成的切片。一般情況下,由于不同切片所切割層序的位置相近,沿層切片和地層切片會(huì)取得相似的儲(chǔ)層預(yù)測(cè)效果[22]。實(shí)際上沿層切片和地層切片度量地層厚度的標(biāo)準(zhǔn)不同。下面具體分析兩種切片技術(shù)的近似定量約束條件。
為了更好地說(shuō)明這個(gè)問題,這里引入Lee[6]的三種厚度概念,即垂直厚度(vertical thickness,VT)、間隔厚度(interval thickness,IT)、恢復(fù)厚度(restoration thickness,RT)。VT即視厚度,IT 即真厚度,RT即垂直于拉平層位的地層厚度。在水平層狀地層模型情況下,三種厚度相等(圖2a);在楔狀傾斜地層模型情況下,三種厚度存在明顯差異(圖2b)。由切片的制作方法可知,水平切片以VT 為度量標(biāo)準(zhǔn),地層切片以IT 為度量標(biāo)準(zhǔn),沿層切片以RT 為度量標(biāo)準(zhǔn)。由幾何關(guān)系可知,IT 和RT 之間的換算公式為

式中θB、θA分別為層序底界面、頂界面與水平面的夾角,則θB-θA為層序頂、底界面的夾角。由式(3)可知,RT 與IT 之間的差異是由IT 和θB-θA決定的,具體為:
(1)當(dāng)θB-θA=0時(shí),即在水平層狀地層模型情況下,RT=IT,這時(shí)沿層切片和地層切片所處位置相同,儲(chǔ)層預(yù)測(cè)結(jié)果一致。
(2)當(dāng)θB-θA≠0時(shí),分兩種情況:①當(dāng)楔狀層序頂、底界面的IT 相等時(shí),隨著θB-θA變大,RT 與IT 的差值逐漸變大,即所處位置的沿層切片與地層切片厚度差值ΔH 越大(圖3),可能導(dǎo)致儲(chǔ)層預(yù)測(cè)結(jié)果存在差異;②當(dāng)θB-θA保持不變時(shí),隨著IT逐漸增大,ΔH 也逐漸增大(圖4),可能導(dǎo)致儲(chǔ)層預(yù)測(cè)結(jié)果存在差異。

圖2 水平層序模型(a)和楔狀層序模型(b)的三種厚度示意圖

圖3 楔狀層序θB-θA 與ΔH 的關(guān)系

圖4 楔狀層序IT 與ΔH 的關(guān)系
在實(shí)際的儲(chǔ)層預(yù)測(cè)中,RT 與IT 的差值主要體現(xiàn)在地震垂向分辨率上。當(dāng)?shù)卣鸫瓜蚍直媛屎芨邥r(shí),RT與IT之間的細(xì)微差別可能被地震切片反映;當(dāng)?shù)卣鸫瓜蚍直媛瘦^低時(shí),RT與IT之間的細(xì)微差別不能被分辨。因此,只有通過限定地震的分辨極限,才能在實(shí)際應(yīng)用中進(jìn)一步討論RT 與IT 的差別。這里取地震的垂向分辨極限為10m,當(dāng)ΔH>10m 時(shí),沿層切片和地層切片的儲(chǔ)層預(yù)測(cè)結(jié)果差異較大;當(dāng)ΔH<10m 時(shí),沿層切片與地層切片的儲(chǔ)層預(yù)測(cè)結(jié)果基本相似。以此為前提條件,分別取θBθA=5°~30°,IT=50~2000m,然后利用式(3)計(jì)算ΔH 和相對(duì)誤差(表1)。可見,當(dāng)薄儲(chǔ)層的最小厚度大于10m、楔狀層序的IT<1000m、θB-θA<15°時(shí),沿層切片與地層切片的儲(chǔ)層預(yù)測(cè)效果相似,兩者可以交替使用。在地震資料解釋中,可以根據(jù)實(shí)際地質(zhì)條件,利用以上方法快速判斷沿層切片和地層切片的儲(chǔ)層預(yù)測(cè)效果,為選取合適的儲(chǔ)層預(yù)測(cè)方法提供依據(jù)。另外,厚度相對(duì)誤差計(jì)算結(jié)果(表1)可為古構(gòu)造恢復(fù)提供定量參考。
鄂爾多斯盆地延長(zhǎng)組長(zhǎng)8—長(zhǎng)7段為一套完整的三級(jí)層序,IT<80m,一般θB-θA<5°,長(zhǎng)8段的水下分流河道砂體厚度一般大于15m。按照上述認(rèn)識(shí),利用沿層切片和地層切片均可預(yù)測(cè)河道砂體厚度。利用三維資料預(yù)測(cè)長(zhǎng)8 段底部的河道砂體厚度,地層切片(圖5a上)和沿層切片(圖5b上)反映了大致相同的砂體平面展布特征,但在局部細(xì)節(jié)上有差異,這些差異不僅與切片制作方法有關(guān),也與層位解釋精度、提取時(shí)窗大小有關(guān),是一種綜合響應(yīng)的結(jié)果。因此,由于地層切片與沿層切片切片制作方法不同造成的差異被其他因素的影響掩蓋,在一定程度上兩者可以相互替代。
古地貌恢復(fù)的常規(guī)方法有沉積學(xué)法、印模法、殘余厚度法、層拉平法、高分辨率層序地層學(xué)法[23]。其中印模法、殘余厚度法、層拉平法都屬于地層厚度法,從地震資料解釋的角度講,地層厚度法的本質(zhì)與地震層拉平技術(shù)存在相通之處,即相當(dāng)于拉平某一沉積標(biāo)志層對(duì)應(yīng)的地震反射層位,利用下伏地層考查地層的古構(gòu)造形態(tài)以反映某一時(shí)期的古地貌。印模法利用古地貌之上充填地層的厚度鏡像地反映古地貌,殘余厚度法利用古地貌之下殘余的地層厚度反映古地貌。當(dāng)古地貌形成階段剝蝕地層的厚度較小時(shí),殘余厚度法與印模法恢復(fù)的古地貌相似,兩者都可用于古地貌恢復(fù);當(dāng)古地貌形成階段地層剝蝕的厚度較大時(shí),兩種方法的恢復(fù)結(jié)果差異很大。

表1 不同層序模型IT與RT數(shù)據(jù)分析統(tǒng)計(jì)表

圖5 鄂爾多斯盆地A 三維區(qū)長(zhǎng)8段底部河道砂巖地層切片(a上)與沿層切片(b上)及相關(guān)地震剖面(a下、b下)
圖6為利用殘余厚度法和印模法得到的鄂爾多斯盆地A 三維區(qū)地層厚度圖及相關(guān)地震剖面。由圖可見,兩種方法恢復(fù)的古地貌形態(tài)差異較大,尤其在工區(qū)西部,殘余厚度法恢復(fù)的古地貌單元為古洼地(圖6a上),印模法恢復(fù)的古地貌單元為古高地(圖6b上),兩種方法的恢復(fù)結(jié)果可謂天壤之別。綜合分析認(rèn)為,印模法恢復(fù)的古地貌更合理,而殘余厚度法由于沒有考慮剝蝕厚度的影響,導(dǎo)致古地貌恢復(fù)誤差大。通過修正殘余厚度法,在恢復(fù)剝蝕地層厚度的基礎(chǔ)上提高古地貌恢復(fù)精度。
這里需要滿足兩個(gè)前提條件:①不整合面上、下地層壓實(shí)差異作用較小,即不整合面之上剝蝕的地層厚度較小或橫向上剝蝕厚度近似相等;②古地貌形成階段沒有發(fā)生過強(qiáng)烈的構(gòu)造變形。如果構(gòu)造變形強(qiáng)烈,則地層的橫向伸縮量在古地貌恢復(fù)中不容忽視,而地震層拉平技術(shù)并不能有效恢復(fù)橫向伸縮量。在滿足上述前提條件下,恢復(fù)古地貌分兩步。首先,恢復(fù)古地貌形成階段地層的剝蝕厚度:①在古地貌上、下位置各選取一個(gè)等時(shí)界面;②選取古地貌上、下位置地層保存最完整的1 口鉆井作為標(biāo)準(zhǔn)井,并將兩個(gè)等時(shí)界面之間鉆井的垂直厚度作為標(biāo)準(zhǔn)厚度;③利用標(biāo)準(zhǔn)厚度減去兩個(gè)等時(shí)界面之間的現(xiàn)存地層厚度即為剝蝕地層厚度。然后,在三維區(qū)內(nèi)恢復(fù)古地貌,即計(jì)算古地貌之下等時(shí)界面與不整合面之間的地層厚度,然后將該厚度與剝蝕地層厚度相加,即代表了剝蝕作用發(fā)生前的古地貌。
鄂爾多斯盆地前侏羅紀(jì)古地貌形成之前,并未發(fā)生強(qiáng)烈的構(gòu)造變形,構(gòu)造運(yùn)動(dòng)以垂直抬升為主。在風(fēng)化剝蝕和河流侵蝕的雙重作用下,形成了千溝萬(wàn)壑的古地貌景觀,隨后逐漸被延安組煤系地層“填平補(bǔ)齊”。在古地貌上、下分別選取延安組延9段煤層底界(TJ9反射層)和延長(zhǎng)組長(zhǎng)3段泥巖底界(TT3反射層)作為兩個(gè)等時(shí)界面。然后利用兩個(gè)等時(shí)界面之間的最大鉆井地層厚度(220m)減去兩個(gè)等時(shí)界面之間的現(xiàn)存地層厚度,得到剝蝕地層厚度(圖7a上)。該計(jì)算過程相當(dāng)于同時(shí)拉平兩個(gè)等時(shí)界面,中間缺失的部分即為剝蝕地層的厚度。將長(zhǎng)3段泥巖底界(TT3反射層)至不整合面之間的地層厚度與剝蝕地層厚度相加,得到前侏羅紀(jì)地層厚度圖(圖7b上)。對(duì)比圖6與圖7可見,修正殘余厚度法(圖7b上)與印模法(圖6b上)恢復(fù)的古地貌非常相似。這是由于以兩個(gè)等時(shí)界面為基準(zhǔn)面,兩種方法恢復(fù)的地層厚度是互補(bǔ)的。如果后期差異壓實(shí)作用強(qiáng)烈,則必須進(jìn)行壓實(shí)校正,還要考慮剝蝕地層對(duì)不整合面之下的殘余地層厚度的影響。因此,壓實(shí)作用對(duì)古地貌恢復(fù)精度的影響要引起充分重視。綜上所述,印模法與殘余厚度法恢復(fù)古地貌的原理及結(jié)果均存在明顯差異。

圖6 利用殘余厚度法(a上)和印模法(b上)得到的鄂爾多斯盆地A 三維區(qū)地層厚度圖及相關(guān)地震剖面(a下、b下)

圖7 利用殘余厚度法(a上)和修正殘余厚度法(b上)得到的鄂爾多斯盆地A 三維區(qū)地層厚度圖及相關(guān)地震剖面(a下、b下)
古構(gòu)造恢復(fù)常用的方法有寶塔圖法、地層厚度法和平衡剖面[24],其中平衡剖面主要針對(duì)構(gòu)造變形強(qiáng)烈的地區(qū)。這里著重討論在構(gòu)造變形不甚強(qiáng)烈的地區(qū),寶塔圖法、地層厚度法與層拉平法的區(qū)別和聯(lián)系。
寶塔圖法是將同一區(qū)塊不同深度的構(gòu)造圖依次相減,利用地層厚度反映古構(gòu)造形態(tài),其原理與地層厚度法相同。從地震資料解釋的角度看,地層厚度法相當(dāng)于“拉平”地層的頂界面,底界面的起伏即直觀地反映了古構(gòu)造形態(tài)。因此,針對(duì)構(gòu)造變形不甚強(qiáng)烈的地區(qū),寶塔圖法、地層厚度法與層拉平法的古構(gòu)造恢復(fù)原理相同,都是利用地層厚度反映古構(gòu)造形態(tài),只是層拉平技術(shù)在表現(xiàn)上更為直觀。
對(duì)于一套連續(xù)沉積的地層而言,如果利用地層厚度恢復(fù)古構(gòu)造,必須考慮壓實(shí)作用對(duì)地層厚度的影響,而壓實(shí)作用校正又涉及剝蝕地層厚度的恢復(fù)。其實(shí)這是一項(xiàng)系統(tǒng)工作,工作量大,而且去壓實(shí)作用的結(jié)果或然性強(qiáng),去壓實(shí)作用校正本身的誤差將會(huì)進(jìn)一步增大古構(gòu)造恢復(fù)誤差。由于現(xiàn)今地震剖面反映的地層厚度代表壓實(shí)作用最強(qiáng)階段的地層厚度,因此如果恢復(fù)該階段的古構(gòu)造,則不用再進(jìn)行去壓實(shí)作用校正的相關(guān)工作,在減少工作量的同時(shí)也提高了古構(gòu)造恢復(fù)精度[25]。
一般情況下,壓實(shí)作用最強(qiáng)的時(shí)期也是烴源巖大量生烴、排烴、油氣大規(guī)模運(yùn)聚的主要時(shí)期[25-26],該時(shí)期的古構(gòu)造形態(tài)與油氣成藏關(guān)系最密切,可作為油氣綜合評(píng)價(jià)的重要依據(jù)。鄂爾多斯盆地在早白堊世中晚期,地層壓實(shí)作用最強(qiáng),長(zhǎng)7段源巖埋深最大,有機(jī)質(zhì)開始大量生烴、排烴及運(yùn)移,形成原生油氣藏[27]。后期雖然經(jīng)歷了一定的構(gòu)造變動(dòng),但改造作用不強(qiáng),次生油氣藏規(guī)模有限。現(xiàn)今地震剖面的地層厚度大致代表生烴、排烴期的地層厚度,恢復(fù)該時(shí)期的古構(gòu)造形態(tài),無(wú)需進(jìn)行去壓實(shí)作用校正,在提高恢復(fù)精度的同時(shí),也增加了古構(gòu)造恢復(fù)的油氣地質(zhì)意義。
在恢復(fù)侏羅系直羅組沉積期的古構(gòu)造形態(tài)時(shí),地層厚度法的含義相當(dāng)于拉平上覆代表該時(shí)期的地震等時(shí)界面,考查下伏延9段底界的構(gòu)造形態(tài),并計(jì)算這兩個(gè)地震反射層位之間的地層厚度,利用地層厚度圖反映古構(gòu)造形態(tài)(圖8a上)。但由于侏羅系直羅組沉積期的古構(gòu)造形態(tài)與成藏關(guān)系不密切,也沒有進(jìn)行去壓實(shí)作用校正,因此恢復(fù)的古構(gòu)造形態(tài)與鉆井?dāng)?shù)據(jù)的吻合程度不高。
通過恢復(fù)早白堊世中期的古構(gòu)造形態(tài)可以增強(qiáng)古構(gòu)造恢復(fù)的油氣地質(zhì)意義。在地震剖面上選取代表該時(shí)期的地震標(biāo)志層,進(jìn)行層拉平后考查下伏延9段底界反射層的構(gòu)造形態(tài),并計(jì)算這兩個(gè)地震層位的地層厚度,得到生、排烴期的地層厚度圖(圖8b上)。該時(shí)期的古構(gòu)造形態(tài)與鉆井?dāng)?shù)據(jù)吻合度高,油流井幾乎都位于古構(gòu)造高部位,而且無(wú)需進(jìn)行去壓實(shí)校正,恢復(fù)精度高。該時(shí)期的古厚度圖可作為井位部署的主要依據(jù)。
平衡剖面是分析復(fù)雜構(gòu)造變形區(qū)構(gòu)造演化的有效方法之一,利用構(gòu)造變形過程中巖石體積守恒的原理恢復(fù)不同時(shí)期的古構(gòu)造,然后明確構(gòu)造演化過程。平衡剖面能恢復(fù)地層的原始長(zhǎng)度,進(jìn)而計(jì)算構(gòu)造伸縮量,而地震層拉平技術(shù)在構(gòu)造恢復(fù)中地層長(zhǎng)度保持不變,這是兩者的最大區(qū)別。是否可用地震層拉平技術(shù)代替平衡剖面分析構(gòu)造演化呢?
在構(gòu)造變形不甚強(qiáng)烈的地區(qū),由于斷層數(shù)量較少、斷距較小、地層的橫向伸縮量不大,因此由平衡剖面的“去斷層”及“去褶皺”操作得到的構(gòu)造伸縮量很小,不會(huì)引起構(gòu)造形態(tài)發(fā)生根本變化,而且“去壓實(shí)”過程與地震層拉平?jīng)]有本質(zhì)區(qū)別。另外,構(gòu)造演化分析的地質(zhì)目標(biāo)尺度較大,橫向伸縮量和壓實(shí)作用對(duì)恢復(fù)結(jié)果的影響一般可以忽略。因此在構(gòu)造變形不甚強(qiáng)烈的地區(qū),可以近似地利用地震層拉平技術(shù)代替平衡剖面分析構(gòu)造演化。
以鄂爾多斯盆地天環(huán)坳陷構(gòu)造演化分析為例。首先,選取侏羅系底部煤層強(qiáng)反射(TJ9)、白堊系底部礫巖強(qiáng)反射(TK)兩個(gè)區(qū)域等時(shí)界面作為地震拉平層位。然后,自下而上進(jìn)行拉平:①拉平TJ9反射層,即恢復(fù)侏羅系延9 段的古地貌,考查拉平層至長(zhǎng)7段泥巖頂界(TT7)的地層厚度變化。恢復(fù)結(jié)果表明,當(dāng)時(shí)地層最厚的位置處于剖面東部,即延長(zhǎng)組湖盆的沉積中心,剖面兩側(cè)的前積反射可作為佐證(圖9a),這表明該時(shí)期天環(huán)坳陷還沒有形成。②拉平TK反射層,即恢復(fù)早白堊世的古地貌,考查拉平層至侏羅系底界(TJ9)的地層厚度變化。恢復(fù)結(jié)果表明,當(dāng)時(shí)地層最厚的位置處于剖面西部,表明天環(huán)坳陷已經(jīng)形成,最低位置即為天環(huán)坳陷的軸部(圖9b)。③在地震剖面上,考查古近系底界(TE)的構(gòu)造形態(tài),即近似代表現(xiàn)今天環(huán)坳陷的形態(tài)。可以看出,天環(huán)坳陷的軸部位于地震剖面的中部,表明其整體已經(jīng)向東發(fā)生遷移(圖9c)。經(jīng)測(cè)量可知,自中生代以來(lái),受西緣逆沖作用的影響,天環(huán)坳陷軸部向盆地內(nèi)(東部)遷移了約34km,與前人計(jì)算結(jié)果(約為30km)相近[28]。

圖8 由地層厚度法得到的鄂爾多斯盆地B二維區(qū)TJZ-TJ9(a上)、TK-TJ9(b上)厚度圖及相關(guān)地震剖面(a下、b下)

圖9 鄂爾多斯盆地天環(huán)坳陷構(gòu)造演化剖面
目前,地震層拉平技術(shù)已作為一種常規(guī)的主流解釋技術(shù)被用于地震解釋的各個(gè)環(huán)節(jié),彰顯出巨大的應(yīng)用價(jià)值。但由于技術(shù)本身的局限性,或?qū)?yīng)用的前提條件認(rèn)識(shí)模糊、把握不準(zhǔn),在實(shí)際應(yīng)用中還存在一些問題。
本文將地震層拉平技術(shù)與地震切片技術(shù)、印模法、殘余厚度法、寶塔圖法及平衡剖面技術(shù)進(jìn)行對(duì)比、分析,找出技術(shù)原理上的異同點(diǎn),然后通過增加地質(zhì)約束條件或修正原有技術(shù)方案,合理而科學(xué)地將地震層拉平技術(shù)應(yīng)用于儲(chǔ)層預(yù)測(cè)、古地貌恢復(fù)、古構(gòu)造恢復(fù)及構(gòu)造演化分析等四個(gè)方面。以鄂爾多斯盆地二維、三維地震區(qū)塊為例進(jìn)行分析、驗(yàn)證,證實(shí)地震層拉平技術(shù)可以在實(shí)際應(yīng)用中取得預(yù)期的效果。
地震層拉平作為一種古地表重建技術(shù),仍有許多潛力有待挖掘。處理與解釋人員緊密結(jié)合,在地震層拉平的基礎(chǔ)上,對(duì)整個(gè)三維數(shù)據(jù)體實(shí)現(xiàn)體拉平,從而更便捷地分析構(gòu)造、沉積演化及油氣成藏史,這將是下一步的努力方向。