基于聲波時差法恢復(fù)延安地區(qū)白堊紀(jì)末期剝蝕量

謝 沖，張鳳奇，李艷霞，謝紅兵

（1.西安石油大學(xué)地球科學(xué)與工程學(xué)院，陜西西安 710065；2.陜西省油氣成藏地質(zhì)學(xué)重點(diǎn)實(shí)驗(yàn)室，陜西西安 710065；3.北京中科鑫宇科技發(fā)展有限公司，北京 100083）

地層剝蝕量恢復(fù)是研究含油氣盆地埋藏史、熱史和油氣成藏史的基礎(chǔ)工作，長期以來都是含油氣盆地分析的熱點(diǎn)和重點(diǎn)問題[1-2]。地層剝蝕現(xiàn)象在沉積盆地中普遍存在，對沉積盆地中油氣的生成、運(yùn)移和保存具有重要影響[3-4]。適宜的構(gòu)造運(yùn)動有利于形成構(gòu)造圈閉，成為油氣運(yùn)移和聚集的有利場所；而強(qiáng)烈的構(gòu)造運(yùn)動則會造成地層嚴(yán)重剝蝕，破壞油氣的保存條件，從而會導(dǎo)致油氣泄漏。因此，恢復(fù)地層剝蝕量對研究盆地構(gòu)造演化史和油氣成藏機(jī)理具有重要意義[2，5]。

延安地區(qū)位于鄂爾多斯盆地陜北斜坡東南緣，該地區(qū)內(nèi)部構(gòu)造簡單，少有斷層或鼻狀構(gòu)造。目前各學(xué)者對延安地區(qū)的研究主要集中于沉積模式[6]、儲層物性特征[7-10]、巖相特征[11]、成藏模式[12]等方面，但對于剝蝕量恢復(fù)的研究較少。根據(jù)已有的研究表明，陜北斜坡自中生代以來，主要發(fā)生了四期不均勻抬升與剝蝕事件，分別是三疊世末期、早侏羅世末期、中侏羅世末期和白堊紀(jì)末期，前三期剝蝕量均較小，后一期剝蝕量較大[13]。前人分別采用聲波時差法、地層對比法[13]和鏡質(zhì)體反射率法[14]對鄂爾多斯盆地四期不均勻抬升事件導(dǎo)致的地層剝蝕量進(jìn)行了計(jì)算，但對延安地區(qū)白堊紀(jì)末期的地層剝蝕量至今仍缺乏精細(xì)研究。為此，本文利用60口井的鉆井、測井、錄井等資料，選取了聲波時差法定量恢復(fù)了延安地區(qū)白堊紀(jì)末期的地層剝蝕量，總結(jié)了延安地區(qū)剝蝕量分布規(guī)律。該研究成果可為該區(qū)油氣勘探及資源評價提供依據(jù)。

1 方法原理及其適用性

1.1 方法的選擇

目前，國內(nèi)外能夠用于恢復(fù)地層剝蝕量的方法眾多[1，15-17]，歸納起來主要分為四類：（1）地球物理學(xué)類，包含聲波時差法[16，18-22]、孔隙度法等[23-24]；（2）地質(zhì)學(xué)類，包含地層對比法[20，25-26]、沉積速率法、沉積波動分析法等[27-28]；（3）地?zé)釋W(xué)類，包含鏡質(zhì)體反射率法[2-3，29-31]、流體包裹體法[32-35]、磷灰石裂變徑跡法等[36-38]；（4）地球化學(xué)類，包含天然氣平衡濃度法、宇宙成因核素法等[39]。

雖然恢復(fù)地層剝蝕量的方法眾多，但是這些恢復(fù)剝蝕量的方法往往有著自身的局限性，需要在一定的地質(zhì)條件下才能夠進(jìn)行操作應(yīng)用?；诼暡〞r差法具有應(yīng)用原理成熟可靠，獲取資料迅速，操作步驟簡單易行等特點(diǎn)，因此本文選擇應(yīng)用聲波時差法對白堊系地層的剝蝕量進(jìn)行估算。

1.2 聲波時差法及其原理

泥巖壓實(shí)外推法是MAGARA[40]在1976 年提出的一種方法，認(rèn)為用泥巖聲波時差法恢復(fù)剝蝕量時聲波時差與地層埋深的關(guān)系可以用下式表征：

式中：Δt0-研究區(qū)的地表聲波時差值，理論值為600～650 μs/m，針對一個地區(qū)，可根據(jù)正常壓實(shí)趨勢線向上外推至地表求得；Δt-深度H 處的泥巖聲波時差值，μs/m；C-正常壓實(shí)曲線斜率；H-泥巖埋深，m。

聲波時差法的原理為：在地層有剝蝕的區(qū)域，當(dāng)剝蝕面以上沉積物的厚度小于地層被剝蝕的厚度時，將剝蝕面以下泥巖的壓實(shí)趨勢線向上外延至古地表（原始地表）Δt0處，古地表與剝蝕面之間的距離為所求的剝蝕量[41]（圖1）。

圖1 聲波時差法推算地層剝蝕量原理示意圖Fig.1 The schematic diagram of the principle of calculating the amount of erosion by interval transit time method

1.3 適用條件

該方法的使用需要滿足一定的適用性要求。一般情況下，使用該方法的前提條件是剝蝕面上部的新地層對其以下殘留地層所施加的壓力應(yīng)小于或等于被剝蝕地層在剝蝕前對殘留地層所施加的壓力，其目的是為了保證下部殘留地層的泥巖壓實(shí)曲線特征沒有受到破壞和改造。當(dāng)剝蝕面上部新沉積的地層厚度大于下部殘留地層的厚度時，這種情況下可以依據(jù)上部新地層的壓實(shí)曲線斜率（C1）與下部殘留地層的壓實(shí)曲線斜率（C2）的關(guān)系來進(jìn)行判定。

（1）當(dāng)C2=C1時，下部地層的壓實(shí)率C2位于上部地層壓實(shí)率C1左側(cè)或在同一條直線上時可以使用聲波時差法恢復(fù)剝蝕量。

（2）當(dāng)C2>C1時，可以使用聲波時差法恢復(fù)剝蝕量。

（3）當(dāng)C2

從研究區(qū)鉆井資料可得，在白堊系上部新沉積的地層厚度在350 m 上下，其沉積厚度遠(yuǎn)遠(yuǎn)小于其已經(jīng)被剝蝕掉的地層厚度。因此可以選用聲波時差法恢復(fù)剝蝕量。鄂爾多斯盆地延安地區(qū)地表聲波時差一般取600 μs/m[13，42]。

2 實(shí)例分析

2.1 鄂爾多斯盆地延安地區(qū)地質(zhì)背景

鄂爾多斯盆地位于華北地臺的西部，面積約為25×104km2。依據(jù)盆地的現(xiàn)今構(gòu)造形態(tài)，可將鄂爾多斯盆地劃分為六個一級構(gòu)造單元，分別為北部的伊盟隆起、南部的渭北隆起、西部的天環(huán)坳陷和西緣逆沖帶、東部的晉西撓褶帶、中部的陜北斜坡（伊陜斜坡）。陜北斜坡作為盆地內(nèi)最大的一個構(gòu)造單元，面積約為9×104km2。延安地區(qū)位于鄂爾多斯盆地陜北斜坡東南緣，地跨陜西省延安市、志丹縣、延長縣、甘泉縣和富縣等地區(qū)，面積約為3.3×104km2。地面的平均海拔在1 000～2 000 m，相對高差在100～500 m。該地區(qū)內(nèi)部構(gòu)造簡單，地層相對平緩，為向西傾伏且傾角不到1°的單斜構(gòu)造。

2.2 聲波時差法恢復(fù)剝蝕量

由于延安地區(qū)位于鄂爾多斯盆地的東部，而對于盆地東部地區(qū)來講，該地區(qū)白堊統(tǒng)的地層厚度較薄甚至已經(jīng)被剝蝕，能夠制作的泥巖壓實(shí)曲線少且短，僅僅使用這些數(shù)據(jù)想要得到比較準(zhǔn)確的白堊紀(jì)末期的剝蝕量是很難的，因?yàn)榘讏紫蹬c它的下部地層（主要指代的是三疊系-侏羅系）之間存在著沉積間斷（不連續(xù)），但通過研究區(qū)的大量實(shí)際資料得到在該間斷期后的下白堊統(tǒng)沉積比較厚的地層區(qū)域，它的下部地層處于一個超補(bǔ)償?shù)臓顟B(tài)。在間斷前下部地層在上部地層白堊系的荷載壓實(shí)作用下，壓實(shí)程度和現(xiàn)今地層目前的實(shí)際壓實(shí)情況已無太大差別，并且在后期也沒有遭到破壞。因此，恢復(fù)延安地區(qū)白堊紀(jì)末期的剝蝕量也可利用下部地層三疊系-侏羅系的泥巖正常壓實(shí)曲線進(jìn)行求出[42]。使用聲波時差法恢復(fù)地層剝蝕量的關(guān)鍵是讀值和泥巖段的準(zhǔn)確選取。結(jié)合鉆井、測井與錄井等資料，選擇延安地區(qū)60 口井的較純泥巖層，采集聲波時差數(shù)據(jù)進(jìn)行制作泥巖壓實(shí)曲線。采集聲波時差數(shù)據(jù)時遵循以下3條原則：（1）讀取數(shù)據(jù)時應(yīng)避免薄層泥巖段，應(yīng)選取厚度大于5 m 的純泥巖段；（2）為了避免誤差，應(yīng)讀取被選泥巖段聲波時差平均值；（3）讀取的泥巖段不能有明顯的井徑擴(kuò)大現(xiàn)象，因?yàn)榫畯綌U(kuò)大會造成聲波時差相應(yīng)增大。

以WY1 井為例對單井恢復(fù)進(jìn)行說明（圖2）。選取白堊系中的正常壓實(shí)段進(jìn)行線性回歸擬合，將其泥巖壓實(shí)趨勢線向上外推與對數(shù)橫坐標(biāo)相交于600 μs/m（600 μs/m 為地表聲波時差），將此點(diǎn)繼續(xù)向上順其延伸即為古地表（原始地表），現(xiàn)今殘留的下白堊統(tǒng)地層頂面到原始地表之間的垂直距離即為所求的地層剝蝕量，可以通過該單井恢復(fù)圖的深度縱坐標(biāo)直接讀出剝蝕量值，為了確保求出的剝蝕量值正確合理，筆者使用了泥巖壓實(shí)曲線的聲波時差（AC）和埋深（H）方程式進(jìn)行了計(jì)算，WY1 井泥巖壓實(shí)曲線的線性公式為y=-1.834 3x+874.89，其中：x-地表聲波時差；y-現(xiàn)今地表至古地表的厚度，將x=600 μs/m 代入y=-1.834 3x+874.89 計(jì)算后得到y(tǒng) 等于226 m，而WY1 井第四系地層厚度為240 m，兩者相加為466 m，即可計(jì)算得到WY1 井地層剝蝕量為466 m。相關(guān)系數(shù)為0.901 5，相關(guān)性較高。采用同樣的方法對研究區(qū)60 口探井的單井剝蝕量進(jìn)行了恢復(fù)。

圖2 延安地區(qū)WY1 井地層剝蝕量估算Fig.2 Estimation of formation erosion of well WY1 in Yan'an area

2.3 連井對比

根據(jù)研究區(qū)部分單井聲波時差-深度回歸方程（表1）和恢復(fù)的剝蝕量圖（圖3），以南西-北東、北西-南東方向兩條剖面為例，南西-北東方向選取的單井FY2井、YY10 井、X57 井、Y311 井、YY3 井，求取的剝蝕量分別為791 m、897 m、1 038 m、1 371 m、1 453 m；北西-南東方向選取的單井YY2 井、YY508 井、YY509井、Y310 井、YY20 井，求取的剝蝕量分別為614 m、985 m、1 062 m、1 349 m、1 624 m；可以看出南西-北東方向的地層剝蝕量從791 m 增大到了1 453 m，北西-南東方向的地層剝蝕量從614 m 增大到了1 624 m，可以得到研究區(qū)剝蝕量自西向東是逐漸增大的，反映出延安地區(qū)剝蝕量呈現(xiàn)自西向東逐漸增大的趨勢。

圖3 聲波時差法恢復(fù)延安地區(qū)地層剝蝕量Fig.3 Recovery of strata erosion in Yan'an area by interval transit time method

表1 延安地區(qū)部分探井聲波時差-深度回歸方程Tab.1 Acoustic transit time-depth regression equation of some exploration wells in Yan'an area

2.4 剝蝕量平面分布特征

本文通過對延安地區(qū)60 口探井的單井剝蝕量估算分析，得到了該地區(qū)白堊紀(jì)末期的地層剝蝕總量。研究區(qū)白堊紀(jì)末期的剝蝕量具有以下規(guī)律：靠研究區(qū)西部邊界，剝蝕量由西向東是逐漸增大的，分布在400～1 800 m。剝蝕量較大的區(qū)域在清澗、延川、延長等東部區(qū)域及宜川的東南部區(qū)域，地層剝蝕量可高達(dá)1 300 m 以上，在志丹以西區(qū)域，地層剝蝕量逐漸降低，在400 m 以下，恢復(fù)結(jié)果反映出延安地區(qū)白堊紀(jì)末期東部遭受剝蝕較強(qiáng)，向西逐漸減弱的總體格局。

2.5 不同方法估算結(jié)果對比

本文估算結(jié)果與已有文獻(xiàn)相比較，延安地區(qū)白堊紀(jì)末期地層剝蝕的分布與強(qiáng)度與前人得出的結(jié)果較一致，即呈現(xiàn)由西向東剝蝕量逐漸增大的趨勢。如翁望飛等[42]利用聲波時差法計(jì)算得出鄂爾多斯盆地東部、東南部剝蝕量最大，為1 400 m 以上，且向西逐漸減小。黃志剛[43]利用聲波時差、EASY%Ro 等方法恢復(fù)了鄂爾多斯盆地東南緣中生代晚期以來的剝蝕量，認(rèn)為其東部剝蝕量較大，最大可達(dá)到2 000 m 以上。陳瑞銀等[13]以地層對比法為主，結(jié)合聲波時差法，認(rèn)為鄂爾多斯盆地地層剝蝕量由西往東是逐漸增大的，在東部地區(qū)剝蝕量達(dá)到1 200 m 以上。高勝利等[14]利用鏡質(zhì)體反射率法恢復(fù)剝蝕量，認(rèn)為鄂爾多斯盆地東部剝蝕量較大、東南部剝蝕量最大，達(dá)到1 400 m 以上。任戰(zhàn)利[44]利用鏡質(zhì)體反射率法、流體包裹體法、磷灰石裂變徑跡法等方法研究了鄂爾多斯盆地的熱演化史，認(rèn)為白堊紀(jì)沉積之后的剝蝕量在陜北斜坡較大，在1 000～2 000 m。從上述前人的研究成果來看，計(jì)算結(jié)果或大或小存在著一定的差異，但整體上都體現(xiàn)出該研究區(qū)在白堊紀(jì)末期的剝蝕量較大，且在該研究區(qū)的整體趨勢是西部剝蝕量較小，東部剝蝕量較大的特征。計(jì)算結(jié)果存在差異，分析其原因主要包括兩個方面。第一個原因是選取的剝蝕量恢復(fù)方法不同，學(xué)者分別采用了聲波時差法、EASY%Ro、地層對比法、鏡質(zhì)體反射率法、流體包裹體法、磷灰石裂變徑跡法等方法對盆地東部的剝蝕量進(jìn)行了恢復(fù)估算，不同方法之間存在優(yōu)劣性，最后得出的計(jì)算結(jié)果也會稍有差異；第二個原因是巖層本身性質(zhì)影響，聲波時差測井曲線幅度值的大小與地層的巖性、孔隙、密度等因素關(guān)系密切[42]，這些因素都會對聲波時差數(shù)據(jù)的準(zhǔn)確性產(chǎn)生影響，從而會對最終的剝蝕量恢復(fù)結(jié)果產(chǎn)生影響。

3 結(jié)論

利用聲波時差法恢復(fù)延安地區(qū)白堊紀(jì)末期的地層剝蝕量，同時利用連井對比驗(yàn)證，恢復(fù)結(jié)果是可信的。結(jié)果表明：鄂爾多斯盆地延安地區(qū)白堊紀(jì)末期的剝蝕強(qiáng)度大致可分為東部的強(qiáng)剝蝕區(qū)和西部的弱剝蝕區(qū)，主要分布在400～1 800 m。延安地區(qū)剝蝕量較大的區(qū)域在清澗、延川、延長等東部區(qū)域，地層剝蝕量在1 300 m 以上，由東往西剝蝕量逐漸降低，降低為400 m以下。反映出延安地區(qū)白堊紀(jì)末期的剝蝕強(qiáng)度呈現(xiàn)東部剝蝕較強(qiáng)，西部剝蝕程度較弱的特征。