孫 剛 隋建波 劉晶晶

(中國(guó)地震應(yīng)急搜救中心，北京 100049)

0 引言

唐山市是我國(guó)地震重點(diǎn)監(jiān)視區(qū)，歷史上曾發(fā)生1976年7.8級(jí)城市直下型地震，導(dǎo)致嚴(yán)重的人員傷亡和巨大的經(jīng)濟(jì)損失，整個(gè)唐山市近乎被毀滅［1］。近年來(lái)唐山市經(jīng)濟(jì)發(fā)展迅速，城市擴(kuò)大快，城市急需開(kāi)展具有防御城市直下型地震災(zāi)害能力的城市規(guī)劃與建設(shè)。

本文利用三維可視化地質(zhì)建模軟件將收集的唐山市鉆孔資料進(jìn)行數(shù)值模擬，給出唐山市第四系地下三維結(jié)構(gòu)模型，分析唐山市可能的地下斷層及活動(dòng)性。

1 鉆孔數(shù)據(jù)選取

本文收集和整理唐山市工程地質(zhì)、水文地質(zhì)、城市建設(shè)以及河北省城市活斷層與地震危險(xiǎn)性評(píng)價(jià)項(xiàng)目(唐山市項(xiàng)目)鉆探的鉆孔資料，共計(jì)117個(gè)場(chǎng)地，1 573個(gè)鉆孔，通過(guò)篩選整理，確定最終用于三維建模的鉆孔共199個(gè)。對(duì)這199個(gè)鉆孔進(jìn)行地層劃分，建立唐山市第四紀(jì)鉆孔數(shù)據(jù)庫(kù)(見(jiàn)表1)。對(duì)各類(lèi)鉆孔進(jìn)行統(tǒng)一數(shù)字化，形成數(shù)字化的鉆孔分布圖。唐山市區(qū)主要發(fā)育三條第四紀(jì)斷裂，分別是陡河斷裂(F1)、巍山—豐南斷裂(F2)、唐山地震斷裂帶(F3)。

表1 鉆孔層序分層數(shù)據(jù)表(部分)

2 三維建模流程

本文建模流程主要包括基礎(chǔ)資料的處理、鉆孔數(shù)據(jù)庫(kù)的建立、三維建模數(shù)據(jù)體的生成、三維建模與成果輸出等。建模采用ArcGIS和MVS軟件完成，主要考慮ArcGIS的地理信息處理與分析功能，MVS軟件的地質(zhì)地球物理數(shù)據(jù)綜合分析能力較強(qiáng)，有利于后續(xù)三維綜合建模和數(shù)據(jù)庫(kù)建設(shè)。

建模流程圖見(jiàn)圖1。

圖1 三維建模流程圖

3 三維建模成果輸出

3.1 三維數(shù)據(jù)體的生成

MVS對(duì)三維地質(zhì)數(shù)據(jù)的格式有一定要求，本文首先是將處理好的層序分層數(shù)據(jù)，通過(guò)程序處理為MVS可以導(dǎo)入的數(shù)據(jù)格式。這一數(shù)據(jù)格式主要包括兩個(gè)部分:頭文件區(qū)和數(shù)據(jù)區(qū)，如圖2所示，頭文件區(qū)給出了每列數(shù)據(jù)的屬性，數(shù)據(jù)區(qū)為鉆孔數(shù)據(jù)內(nèi)容，每行代表一個(gè)鉆孔，每列代表鉆孔的一種屬性，如坐標(biāo)、高程、層底埋深等［2-4］。

圖2 MVS三維建模格式（部分）

3.2 數(shù)據(jù)導(dǎo)入與插值

將上述經(jīng)處理的數(shù)據(jù)體，導(dǎo)入MVS軟件中，形成空間的點(diǎn)，這些點(diǎn)具有高程和各層序界面埋深等屬性。在MVS中可對(duì)這些點(diǎn)的高程或某一數(shù)值屬性進(jìn)行插值，形成地表和各層序界面的空間曲面。

3.3 三維模型的建立

經(jīng)過(guò)曲面擬合(克里金插值)之后，可以以各層面為基礎(chǔ)構(gòu)建三維地質(zhì)體，各幾何體代表了不同層序組的厚度空間變化，形成的地質(zhì)體見(jiàn)圖3。在三維地質(zhì)體的基礎(chǔ)上，可以生成三維剖面，沿南北、東西或其他任意方向，剖面可以使我們?cè)诓煌嵌扔^察各層位之間的關(guān)系，包括空間形態(tài)和厚度變化，更好的輔助沉積和構(gòu)造演化分析，圖4是一組剖面的示意圖。最終形成唐山市區(qū)層序—構(gòu)造格架三維地質(zhì)體、各層序界面的空間起伏與厚度變化圖［5，6］。

4 三維建模成果分析

圖3 建模生成的三維地質(zhì)體

圖4 一組三維剖面

4.1 Qh構(gòu)造空間展布分析

圖5為Qh的側(cè)視圖。

從圖5中可以看出:Qh的底界整體上較為平緩，沉積厚度的變化不大。陡河斷裂(F1)、巍山—豐南斷裂(F2)、唐山地震斷裂帶(F3)對(duì)地層的變化的控制不是很明顯，說(shuō)明三條斷裂在Qh以來(lái)沒(méi)有活動(dòng)的跡象。

圖5 Qh構(gòu)造底板起伏形態(tài)（側(cè)視）

圖6 構(gòu)造底板起伏形態(tài)（側(cè)視）

而唐山地震斷裂帶(F3)的東側(cè)的沉降區(qū)域非常明顯，斷層兩側(cè)沉積基底有較大的變化，說(shuō)明唐山地震斷裂帶(F3)在此時(shí)較為活動(dòng)，對(duì)地層沉積都有一定的控制作用。

圖7 構(gòu)造底板起伏形態(tài)（側(cè)視）

5 結(jié)語(yǔ)

本文收集、整理和分析了唐山市區(qū)199個(gè)典型鉆孔的資料，對(duì)鉆孔巖性分層和層序進(jìn)行分層并建立了數(shù)據(jù)庫(kù)。根據(jù)典型鉆孔數(shù)據(jù)庫(kù)和層序—構(gòu)造格架，基于三維構(gòu)造建模MVS和ArcGIS軟件平臺(tái)，建立了唐山市目標(biāo)區(qū)第四系層序—構(gòu)造格架三維模型，給出包括唐山市目標(biāo)區(qū)層序—構(gòu)造格架三維地質(zhì)體、各層序界面的空間起伏與厚度變化、典型地層剖面等圖件，實(shí)現(xiàn)了目標(biāo)區(qū)三維地質(zhì)建模與可視化。

對(duì)唐山市區(qū)晚第四紀(jì)地層各層段特征、地層—沉積格架和三維建模的綜合分析表明，陡河斷裂(F1)、巍山—豐南斷裂(F2)對(duì)晚第四紀(jì)地層的控制不是很明顯，說(shuō)明它們以來(lái)活動(dòng)較弱。唐山地震斷裂帶(F3)在以前，對(duì)地層的控制較為明顯，斷裂兩側(cè)變化較大，表現(xiàn)為斷裂西側(cè)隆起，東側(cè)凹陷，體現(xiàn)出唐山地震斷裂帶(F3)為以來(lái)斷層活動(dòng)性及其對(duì)沉積的控制作用。

［1］劉 亢.唐山地區(qū)強(qiáng)震復(fù)發(fā)特征與發(fā)震構(gòu)造深淺部耦合關(guān)系綜合研究［D］.北京:中國(guó)地震局地質(zhì)研究所博士學(xué)位論文，2011:5-10.

［2］曹代勇，李青元，朱小弟，等.地質(zhì)構(gòu)造三維可視化模型探討［J］.地質(zhì)與勘探，2001，37(4):60-62.

［3］曹代勇，王占剛.三維地質(zhì)模型可視化中直接交互的實(shí)現(xiàn)［J］.中國(guó)礦業(yè)大學(xué)學(xué)報(bào)，2004，33(4):78-80.

［4］李亦綱，曲國(guó)勝，陳建強(qiáng).城市鉆孔數(shù)據(jù)地下三維地質(zhì)建模軟件實(shí)現(xiàn)［J］.地質(zhì)通報(bào)，2005(5):68-74.

［5］王家華.克里金地質(zhì)繪圖技術(shù)——計(jì)算機(jī)的模型和算法［M］.北京:石油工業(yè)出版社，1999.

［6］武 強(qiáng)，徐 華.三維地質(zhì)建模與可視化方法研究［J］.中國(guó)科學(xué) D 輯，2004，34(1):54-60.

［7］楊承先.唐山、湯陰斷陷地質(zhì)結(jié)構(gòu)及其活動(dòng)性［J］.地震地質(zhì)，1984，6(3):59-66.

［8］閆純有.唐山地區(qū)地貌與構(gòu)造［J］.河北建筑科技學(xué)院學(xué)報(bào)，2000，17(4):63-65.