李 媛 丁文秀 林 松,3

1 中國地震局地震大地測量重點實驗室，武漢市洪山側(cè)路40號，430071 2 地震預(yù)警湖北省重點實驗室，武漢市洪山側(cè)路40號，430071 3 武漢地震工程研究院有限公司，武漢市洪山側(cè)路40號，430071

古地貌恢復(fù)有助于分析物源區(qū)的剝蝕程度、沉積體系、斷裂分布與構(gòu)造格局的空間配置關(guān)系，對地質(zhì)構(gòu)造、資源勘探有重要的指導(dǎo)意義[1-3]，主要方法有印模法、填平補齊法、沉積學(xué)分析法、層序地層學(xué)法等[4-7]。

北部灣盆地福山凹陷是一個中新生代次級裂谷型凹陷，內(nèi)部多期構(gòu)造疊加，具有發(fā)生強震的構(gòu)造背景[8-9]。本文基于鉆井、測井及三維地震數(shù)據(jù)，采用層序地層學(xué)古地貌恢復(fù)方法，對福山凹陷古近系流沙港組3個關(guān)鍵界面沉積前的古地貌進行恢復(fù)和三維建模，運用“將古論今”的地質(zhì)學(xué)原理，以期為該地區(qū)的地震地質(zhì)構(gòu)造、資源勘探、地震預(yù)報及防震減災(zāi)工作提供基礎(chǔ)支撐。

1 區(qū)域地質(zhì)概況

福山凹陷總面積約為2 920 km2，南部為海南凸起，西接臨高凸起，東臨云龍凸起，西緣的臨高斷裂和東緣的長流斷裂為一級控盆大斷裂(圖1)。古近系流沙港組為一套巨厚的湖相三角洲沉積，根據(jù)層序地層基準面原理和層序地層學(xué)分析技術(shù)，結(jié)合鉆井、測井及覆蓋整個研究區(qū)的高分辨率三維地震資料，在古近系流沙港組等時層序地層格架中由新到老識別出3個沉積層序：流沙港組一段(Els1)、流沙港組二段(Els2)和流沙港組三段(Els3)[8]。

圖1 北部灣盆地福山凹陷構(gòu)造格局示意圖 及回剝模擬測線位置Fig.1 Schematic diagram of structural pattern and location of stripping simulation line in Fushan sag, Beibuwan basin

2 古地貌恢復(fù)原理與方法

在建立高精度層序地層格架的基礎(chǔ)上，利用地質(zhì)鉆井和三維地震數(shù)據(jù)，使用定量回剝技術(shù)[10-11]恢復(fù)流沙港組3個沉積時期的原始地貌形態(tài)，并應(yīng)用三維建模軟件對古地貌進行立體展示。

2.1 鉆井與地震數(shù)據(jù)

本文使用的基礎(chǔ)數(shù)據(jù)由海南福山油田勘探開發(fā)有限責(zé)任公司提供，包括230口鉆井?dāng)?shù)據(jù)、綜合錄井圖、測井解釋結(jié)果等資料及基本覆蓋主要研究區(qū)的最新成果數(shù)據(jù)體(二維地震1 260.6 km，三維地震408.625 km2)(圖1)。

在整理數(shù)據(jù)的基礎(chǔ)上，首先對重點井進行聲波合成地震記錄與層位標定，保證凹陷內(nèi)地震反射特征與地質(zhì)層位統(tǒng)一；然后對三維工區(qū)內(nèi)目的層段7個三級層序界面的地震資料進行解釋，解釋密度達到800 m×800 m，成圖面積累計達3 680 km2；最后選取南北向及東西向共10條聯(lián)絡(luò)測線進行精細構(gòu)造分析(圖1、圖2)，以確保地震解釋的精度和覆蓋面。

圖2 福山凹陷骨干剖面(WE-1)層序地層精細解釋圖Fig.2 Detailed interpretation map of sequence stratigraphy of backbone section(WE-1) in Fushan sag

2.2 測線及觀測點布設(shè)

本文采取“點、線、面”結(jié)合的方法，同時體現(xiàn)“體”和“時”的概念。通過選取的7條二維地震測線和56個觀測點(圖3)，應(yīng)用EBM盆地模擬軟件模擬流沙港組地層的沉積演化過程，分析各個觀測點沉降速率的特征。

圖3 回剝模擬測線及觀測點分布Fig.3 Back stripping simulation survey line and observation point distribution

2.3 沉降史模擬

通過沉降史分析，可以在三維沉積體系空間內(nèi)對盆地構(gòu)造特征及其形成過程進行細致的剖析[12]，具體過程見圖4。在對二維測線進行回剝計算的基礎(chǔ)上，確定各個時期的總沉降量，即該時期的實際盆地深度；再根據(jù)選取的界面和時段，將沉降量數(shù)據(jù)轉(zhuǎn)繪到平面底圖上，描繪沉降量平面等直線，從而定量恢復(fù)得到每個時期的原始地貌形態(tài)。

圖4 回剝法工作流程Fig.4 Workflow of stripping back method

3 古地貌恢復(fù)結(jié)果

3.1 沉降史分析

通過沉降史分析得出，福山凹陷古近系流沙港組沉積時期在垂向上具有幕式沉降特征。二維地震測線WE-1的垂向沉降速率演化結(jié)果(圖5)顯示，早期凹陷西部的沉積速率大于東部，而晚期凹陷東部的沉降速率大于西部。

圖5 福山凹陷流沙港組東西向聯(lián)絡(luò)測線 沉降速率分析結(jié)果Fig.5 Analysis of settlement rate of EW connecting survey line of Liushagang formation in Fushan sag

將相應(yīng)時段的沉降量數(shù)據(jù)轉(zhuǎn)繪到平面等值線上，得出以下認識：

1)Els3時期，沉降速率值主體范圍在280 m/Ma以內(nèi)，沉降中心的數(shù)目增多，主要發(fā)育4個沉降中心，水體發(fā)育和斷層活動均呈繼承性擴張的趨勢[12](圖6(a))。

2)Els2時期，沉降速率值主體范圍在250 m/Ma以內(nèi)，福山凹陷向北部遷移，南部斜坡帶范圍明顯擴大，西部沉降中心持續(xù)性下沉，東部沉降中心繼承性發(fā)育，水體深度達到最大，中部低凸起略顯下降趨勢[12](圖6(b))。

3)Els1時期，凹陷萎縮，整體面積減小，但逐漸向北邊遷移，南部斜坡帶遭受的剝蝕面積達到最大，沉降速率仍呈現(xiàn)為西高東低[12](圖6(c))。

圖6 福山凹陷古近系流沙港組時期總沉降量分布Fig.6 Distribution of total settlement of Paleogene Liushagang formation in Fushan sag

3.2 古地貌恢復(fù)

本文在沉降史分析的基礎(chǔ)上應(yīng)用Petrel軟件對古地貌進行了立體展示。整體而言，福山凹陷為一北斷南超的箕狀斷陷，具有東西分帶、南北分塊的構(gòu)造格局，古地貌呈幕式變化特征，溝谷縱橫發(fā)育，凹凸分割性明顯。從Els3時期到Els1時期，南部斜坡帶剝蝕面積逐漸變大，中部低凸起的分割作用由強減弱再變強，北部地區(qū)持續(xù)擴張，推測這些都與凹陷強烈的裂陷活動相關(guān)(圖7)。

圖7 福山凹陷流沙港組關(guān)鍵界面古地貌立體圖Fig.7 3-D map of structure and palaeogeomorphology at the bottom interface of Liushagang formation in Fushan sag

4 結(jié) 語

本文綜合利用地質(zhì)鉆井和三維地震資料，采用層序地層學(xué)古地貌恢復(fù)方法對福山凹陷古近系流沙港組3個關(guān)鍵界面沉積前的古地貌進行恢復(fù)和三維建模，并從垂向和平面2個維度對福山凹陷古近系時期的構(gòu)造格局特征進行分析。結(jié)果表明：

1)通過沉降史分析可知，早期凹陷西部沉積速率大，晚期凹陷東部沉降速率大；古地貌具有東西分帶、南北分塊的構(gòu)造特征。

2)對比沉降速率平面圖與古地貌圖可以看出，沉降速率中心與凹陷中心一致，沉降速率越大，凹陷深度越深；反之，沉降速率較小的地方則為相對凸起的地帶[12]。

3)垂向上，流沙港組古地貌起伏具有幕式沉降特征。