李自紅 劉鴻福 張 敏 曾金艷

1)太原理工大學(xué)，太原 030024

2)山西省地震局，太原 030021

3)太原大陸裂谷動(dòng)力學(xué)國(guó)家野外科學(xué)觀測(cè)研究站，太原 030025

4)北京市地質(zhì)礦產(chǎn)勘查開發(fā)總公司，北京 100195

0 引言

地震活動(dòng)受新構(gòu)造運(yùn)動(dòng)特別是活動(dòng)斷裂的控制，據(jù)統(tǒng)計(jì)，約有85%～95%的6級(jí)以上地震是由活動(dòng)構(gòu)造尤其是活動(dòng)斷裂引起的，活動(dòng)構(gòu)造與地震的關(guān)系不僅在空間分布上具有一致性，而且活動(dòng)斷裂的運(yùn)動(dòng)特征與大地震中沿?cái)嗔旬a(chǎn)生的同震破裂也具同一性(鄧起東等，2003)。地震與活動(dòng)斷裂的關(guān)系是地震地質(zhì)領(lǐng)域內(nèi)最重要的研究?jī)?nèi)容之一。以往這類研究大多停留在定性分析的層面上，其中一個(gè)重要的原因是缺乏有效的空間計(jì)算工具。地理信息系統(tǒng)(GIS)強(qiáng)大的空間分析功能為定量或半定量地研究二者之間的關(guān)系提供了有效途徑(鄔寧芬等，2004)。

應(yīng)用G1S技術(shù)研究地震與活動(dòng)斷裂的關(guān)系尚處于起步階段。于貴華等(1996)利用GIS系統(tǒng)建立了中國(guó)活動(dòng)斷裂信息咨詢分析系統(tǒng)，實(shí)現(xiàn)了利用空間分析功能查詢活動(dòng)斷裂的信息。石樹中等(1998)基于ARC/INFO系統(tǒng)平臺(tái)，開發(fā)建成北天山活動(dòng)斷裂地理信息系統(tǒng)，實(shí)現(xiàn)了活動(dòng)斷裂與地震活動(dòng)的空間分析與統(tǒng)計(jì)分析等功能。屈春燕等(2000)利用GIS的緩沖區(qū)空間分析功能技術(shù)，對(duì)活動(dòng)斷裂和歷史地震及新生代盆地的相關(guān)性進(jìn)行了初步分析。彭自正等(2002)基于GIS緩沖區(qū)分析功能，對(duì)江西活動(dòng)斷裂分布與地震活動(dòng)關(guān)系進(jìn)行了研究。鄔寧芬等(2004)利用GIS對(duì)上海市活動(dòng)斷裂走向與地震的關(guān)系以及活動(dòng)年代與地震的關(guān)系進(jìn)行了分析。劉峰等(2009)利用GIS方法研究了南北地震帶和中央造山帶交會(huì)區(qū)活動(dòng)斷裂與地震的關(guān)系。上述學(xué)者都是利用二維地理信息系統(tǒng)(2D GIS)的緩沖區(qū)空間分析功能研究地震與活動(dòng)斷裂的關(guān)系。

在信息化高度發(fā)達(dá)的今天，科學(xué)計(jì)算可視化技術(shù)正在迅猛發(fā)展，而三維地理信息系統(tǒng)(3D GIS)也是當(dāng)前GIS發(fā)展和研究的一個(gè)熱點(diǎn)問題。與科學(xué)計(jì)算可視化和3D GIS密切相關(guān)的三維地質(zhì)建模技術(shù)已經(jīng)越來越多地受到國(guó)內(nèi)外學(xué)者的重視，人們已經(jīng)迫切希望在三維空間中解決復(fù)雜的地質(zhì)問題(吳立新等，2003;姚鑫等，2006;張敏等，2007)。地震與活動(dòng)斷裂本身處在三維地質(zhì)空間，地震震源在三維空間上是不均勻分布的，這種不均勻分布如果投影到二維空間就會(huì)把同一條活動(dòng)斷裂產(chǎn)生的地震投影到不同的斷裂上，就會(huì)對(duì)地震與活動(dòng)斷裂的關(guān)系研究產(chǎn)生誤導(dǎo)。因此，嘗試在三維空間研究地震與活動(dòng)斷裂的關(guān)系對(duì)于活動(dòng)斷裂研究無(wú)疑具有重要的意義。

本文正是針對(duì)上述研究課題，結(jié)合“太原市活斷層探測(cè)與地震危險(xiǎn)性評(píng)價(jià)”項(xiàng)目，利用Arc-GIS強(qiáng)大的三維顯示功能，設(shè)計(jì)構(gòu)建了三維地質(zhì)模型的算法，利用AOVBA開發(fā)方式進(jìn)行了程序的開發(fā)。利用科學(xué)計(jì)算可視化技術(shù)和3D GIS技術(shù)相結(jié)合的三維地質(zhì)建模技術(shù)，對(duì)太原市區(qū)內(nèi)主要斷裂進(jìn)行了三維地質(zhì)可視化建模，并基于三維斷層模型在ArcGIS環(huán)境下嘗試性構(gòu)建了斷層的三維空間緩沖區(qū)模型?；谠撃Ｐ屠每臻g分析以及統(tǒng)計(jì)學(xué)原理對(duì)太原市區(qū)內(nèi)主要斷裂與地震進(jìn)行了空間統(tǒng)計(jì)分析，進(jìn)而對(duì)地震與活動(dòng)斷裂空間分布關(guān)系進(jìn)行了研究。

1 三維斷層模型的設(shè)計(jì)

三維斷層模型的設(shè)計(jì)是實(shí)現(xiàn)三維斷層空間立體緩沖區(qū)的基礎(chǔ)，也是三維地層模型的延伸(張敏等，2007)。本文建立的三維斷層模型是以淺層地震探測(cè)資料和斷層調(diào)查原始資料為依托，其數(shù)據(jù)結(jié)構(gòu)采用了《中國(guó)地震活動(dòng)斷層探測(cè)技術(shù)系統(tǒng)技術(shù)規(guī)程》(中國(guó)地震局，2005)中提供的有關(guān)數(shù)據(jù)結(jié)構(gòu)。在此基礎(chǔ)上，建立了三維斷層模型的數(shù)據(jù)結(jié)構(gòu)并實(shí)現(xiàn)了三維斷層模型構(gòu)建與顯示的算法設(shè)計(jì)。

在三維斷層模型數(shù)據(jù)結(jié)構(gòu)的設(shè)計(jì)過程中，主要依據(jù)淺層地震探測(cè)資料。從淺層地震探測(cè)深度剖面上可以獲取斷層面原始數(shù)據(jù)，結(jié)合斷層調(diào)查資料，綜合判定研究斷裂的空間展布，包括其傾向和走向等，在此基礎(chǔ)上對(duì)斷層進(jìn)行三維建模。斷層的空間數(shù)據(jù)庫(kù)中包括Feature類型的斷點(diǎn)表(eqexpl)和最后生成斷層后用于存儲(chǔ)斷層剖面數(shù)據(jù)的生成斷面表(scdm)。斷點(diǎn)表中包含Line_ID、Faultname、Fd_id和Kz_id 4個(gè)控制字段和斷點(diǎn)本身的屬性信息。其中，Line_ID用來表示斷層是在哪條測(cè)線上發(fā)現(xiàn)的，F(xiàn)aultname表示此立體面代表的是哪條斷層，F(xiàn)d_id是斷層上斷點(diǎn)的編號(hào)，Kz_id用來控制斷層上斷點(diǎn)連接的順序(中國(guó)地震局，2005)。在實(shí)際工作中，通常會(huì)對(duì)被探測(cè)的斷層布設(shè)多條淺層地震測(cè)線進(jìn)行探測(cè)，因此就需要在幾個(gè)淺層地震剖面上找同一斷層斷點(diǎn)的信息。并對(duì)控制同一條斷裂的幾個(gè)淺層地震剖面做對(duì)比分析，結(jié)合斷層資料給出的斷層傾向、走向等特征來構(gòu)建三維斷層模型。

三維斷層模型實(shí)現(xiàn)過程可以分為2個(gè)步驟:訪問空間數(shù)據(jù)庫(kù)中的斷層數(shù)據(jù)，將用于構(gòu)建三維斷層模型的斷面點(diǎn)從其中篩選出來并進(jìn)行排序，并將二維斷點(diǎn)重構(gòu)為三維斷點(diǎn);依據(jù)目標(biāo)數(shù)據(jù)建立數(shù)據(jù)結(jié)構(gòu)，將目標(biāo)數(shù)據(jù)存儲(chǔ)在已經(jīng)設(shè)計(jì)好的數(shù)據(jù)結(jié)構(gòu)中并據(jù)此構(gòu)建斷面，最后將斷面數(shù)據(jù)添加到Feature中顯示。

1.1 構(gòu)建三維斷點(diǎn)

構(gòu)建三維斷點(diǎn)就是將同一斷層在不同淺層地震探測(cè)剖面上的斷點(diǎn)依次提取，并將提取的二維斷點(diǎn)重構(gòu)為三維斷點(diǎn)。即依據(jù)要研究的斷層編號(hào)查找控制此斷層的所有測(cè)線，將其從數(shù)據(jù)庫(kù)中篩選出來，篩選出的測(cè)線還需要進(jìn)行排序，以保證斷層面的延伸方向。另外，還需要聲明IZAware接口，以使同步好的斷點(diǎn)具有三維集合屬性。提取斷面點(diǎn)并將提取后的二維斷點(diǎn)進(jìn)行三維重構(gòu)的核心代碼如下:

1.2 生成三維斷面

在提取斷面點(diǎn)后，就可以生成三維斷面。在前面對(duì)斷面和剖面的數(shù)據(jù)進(jìn)行分析后發(fā)現(xiàn)，由于在控制同一條斷層的所有剖面上，斷層的位置大致相同，傾向和走向基本一致，但是具體到斷層面的連接工作中，連接一個(gè)斷面的所有點(diǎn)未必在一個(gè)平面上。如果這條斷層由3個(gè)以上的淺層地震探測(cè)剖面控制，則它的走向很可能不是一條直線而是一條折線，這樣控制斷層面的斷層線放在一個(gè)平面上時(shí)就會(huì)有交叉的現(xiàn)象，所以采用將斷層面的四邊形平面劃分成2個(gè)三角形平面來實(shí)現(xiàn)。

在算法設(shè)計(jì)時(shí)首先聲明了用于存放斷層信息的ProfileType類型的數(shù)據(jù)結(jié)構(gòu)，里面存放有斷層號(hào)(dcid)、測(cè)線號(hào)(cxid)、斷點(diǎn)號(hào)(ddid)等信息。同時(shí)定義用于暫時(shí)存放斷層信息的Profile-Type類型的動(dòng)態(tài)數(shù)組結(jié)構(gòu)。實(shí)際繪制斷面時(shí)，將預(yù)先存儲(chǔ)在動(dòng)態(tài)數(shù)組中的數(shù)據(jù)分別賦給點(diǎn)p1、p2、p3、p4，通過對(duì)2個(gè)三角形的分別繪制達(dá)到對(duì)斷面繪制的目的(圖1)。核心代碼如下:

圖1 三維斷層模型示意圖Fig.1 Schematic diagram of 3D fault model.

2 三維空間緩沖區(qū)的建立

在三維斷層模型構(gòu)建的基礎(chǔ)上，設(shè)計(jì)出用于實(shí)現(xiàn)三維空間緩沖區(qū)的數(shù)學(xué)模型。三維空間緩沖區(qū)建立的目的是，在三維空間統(tǒng)計(jì)分析地震位于指定緩沖區(qū)內(nèi)的數(shù)量，從而統(tǒng)計(jì)分析在距離活斷層一定距離的范圍內(nèi)，地震發(fā)生的歷史情況和定性地分析斷層的活動(dòng)情況。

數(shù)學(xué)模型設(shè)計(jì)的大體思路是根據(jù)已經(jīng)構(gòu)建好的斷層模型，對(duì)斷層面上的點(diǎn)進(jìn)行向量擬合，再根據(jù)擬合向量計(jì)算與其垂直的法向量。緩沖區(qū)平面上點(diǎn)的坐標(biāo)依據(jù)歸一化后的法向量和緩沖半徑結(jié)合斷面上的點(diǎn)，依據(jù)法向量移動(dòng)獲得的點(diǎn)坐標(biāo)計(jì)算得到，進(jìn)而構(gòu)建緩沖區(qū)平面。斷層正反2個(gè)方向上的緩沖區(qū)平面都依據(jù)此方法構(gòu)建。

2.1 求解斷面法向量

求解斷面法向量是從訪問數(shù)據(jù)庫(kù)中的斷點(diǎn)表作為切入點(diǎn)，對(duì)斷面上的斷點(diǎn)信息進(jìn)行訪問，篩選出目標(biāo)斷點(diǎn)并獲取斷點(diǎn)的三維空間(x，y，z)坐標(biāo)值。在求解所有斷面法向量的過程中，需要根據(jù)淺層地震探測(cè)剖面上的控制點(diǎn)尋找斷面斷點(diǎn)順序，即將按序連接斷面的上斷點(diǎn)的屬性找出，并依據(jù)此屬性將斷點(diǎn)表中的斷點(diǎn)的信息讀取到動(dòng)態(tài)數(shù)組中。斷面四邊形被視為2個(gè)三角形對(duì)待，2個(gè)三角形上的2個(gè)點(diǎn)通過執(zhí)行程序形成2個(gè)方向向量，2個(gè)方向向量叉乘得到垂直于它們的法向量。在向量叉乘的函數(shù)過程中，用到一個(gè)重要接口Ivector3D，此接口提供向量叉乘的功能。

自定義BufVector型數(shù)據(jù)結(jié)構(gòu)用于存放方向向量的分量值，ProfileType型動(dòng)態(tài)數(shù)組gSufDMALL()用于存放所有斷面數(shù)據(jù)，BufVector型動(dòng)態(tài)數(shù)組gFvALL()用于存放所有矢量，公共變量gDMallNum表示所有斷面數(shù)量，方向向量和垂直于方向向量的法向量定義為BufVector型，求解斷面法向量的核心代碼如下所示:

2.2 構(gòu)建三維緩沖區(qū)

在求解斷面法向量的基礎(chǔ)上構(gòu)建三維緩沖區(qū)。首先對(duì)斷面法向量進(jìn)行歸一化，再依據(jù)緩沖區(qū)半徑將斷面上的點(diǎn)沿垂直于斷面的法向量平移到緩沖區(qū)平面所在的位置，進(jìn)而構(gòu)建三維緩沖區(qū)平面。構(gòu)建時(shí)把四邊形平面分解為2個(gè)三角形平面。構(gòu)建三維緩沖區(qū)核心代碼如下:

'緩存區(qū)面一:

'緩存區(qū)面二:

緩沖區(qū)半徑?jīng)Q定了最后繪制的緩沖區(qū)的范圍，當(dāng)緩沖區(qū)半徑取300m時(shí)，三維空間緩沖區(qū)模型如圖2所示。圖中，紅色的立體面代表斷層，在其兩側(cè)的深灰色立體面是距離斷層面300m的立體緩沖區(qū)。有了基于斷層面的三維空間緩沖區(qū)，依據(jù)斷層和三維緩沖區(qū)模型分析斷層與地震的關(guān)系也就有了基礎(chǔ)。

3 地震與活動(dòng)斷裂三維空間分析

3.1 分析算法設(shè)計(jì)

分析地震與活動(dòng)斷裂的三維空間關(guān)系，在程序設(shè)計(jì)時(shí)可以簡(jiǎn)單的等同于判斷地震與活斷層三維空間緩沖區(qū)的位置關(guān)系，即地震是否在活斷層緩沖區(qū)內(nèi)。本文采用的是向量判別法，即通過判斷向量點(diǎn)乘的結(jié)果進(jìn)而得到三維空間內(nèi)點(diǎn)和平面之間空間關(guān)系的算法。

前面已經(jīng)討論過斷層三維空間緩沖區(qū)的構(gòu)建方法，一個(gè)斷層兩側(cè)的三維空間緩沖區(qū)是由若干對(duì)三角形表面所組成的，那么判斷地震是否在三維空間緩沖區(qū)內(nèi)就是判斷地震是否位于構(gòu)成同一緩沖區(qū)的三角形之間。取其中一對(duì)三角形表面(圖3)，由繪制緩沖區(qū)算法中計(jì)算法向量的算法，可以得到緩沖區(qū)上下表面的法向量V n1和V n2(方向相反，并且朝向緩沖區(qū)外側(cè))。地震震源P z1已知，緩沖區(qū)上下表面任意點(diǎn)到地震震源P z1形成的向量為V p1、V p2，V n1·V p1＞0(V n1和 V p1 夾角 ＜90°)，P z1在緩沖區(qū)上表面法向量同側(cè);V n2·V p2＜0(V n2和 V p2夾角 ＞90°)，P z1在緩沖區(qū)下表面法向量異側(cè)，所以P z1在緩沖區(qū)外側(cè)。反之，V n1·V p1＜0，P z1在緩沖區(qū)上表面法向量異側(cè);V n2·V p2＞0，P z1在緩沖區(qū)下表面法向量同側(cè)，P z1也在緩沖區(qū)外側(cè)。如圖3所示，地震震源P z2已知，緩沖區(qū)上下表面任意點(diǎn)到地震震源P z2形成的向量為 V p3、V p4，當(dāng) V n1·V p3 ＜0 并且 V n2·V p4＜0，P z2位于2個(gè)緩沖區(qū)表面之間，即地震在斷層三維空間緩沖區(qū)內(nèi)。在實(shí)際編程中，上、下表面的任意點(diǎn)，選擇的是三角形的頂點(diǎn)，即斷面控制點(diǎn)平移后形成的點(diǎn)。核心代碼如下:

'緩存區(qū)面上表面 三角形二:

圖2 三維空間緩沖區(qū)模型Fig.2 3D space model of the buffer.

圖3 地震震源與緩沖區(qū)關(guān)系圖Fig.3 Relationship diagram between focuses and the buffer.

'緩存區(qū)面下表面 三角形二:

3.2 實(shí)例分析

三維空間緩沖區(qū)的設(shè)計(jì)和實(shí)現(xiàn)方法具有通用性，基于斷層模型依據(jù)上述方法都可實(shí)現(xiàn)三維空間緩沖區(qū)的建立和基于緩沖區(qū)進(jìn)行定量分析。本文利用“太原市活斷層探測(cè)與地震危險(xiǎn)性評(píng)價(jià)”項(xiàng)目資料，對(duì)太原市區(qū)主要斷裂進(jìn)行了三維空間緩沖區(qū)模型的建立，并基于此模型對(duì)斷裂與地震空間分布的關(guān)系進(jìn)行了初步定量分析研究。

研究區(qū)有儀器記錄以來有震源深度的地震共891個(gè)，震級(jí)為1.0～5.0級(jí)，震源深度為1～30km，分布范圍廣泛，基本遍及研究區(qū)范圍。由于該項(xiàng)目淺層地震探測(cè)深度有限，依據(jù)淺層地震探測(cè)深度剖面構(gòu)建的斷層面深度也有限，不能真實(shí)反映斷層的空間展布和地震的分布關(guān)系，所以在進(jìn)行斷層與地震關(guān)系的空間分析時(shí)，結(jié)合其他成果對(duì)斷層面進(jìn)行了適當(dāng)延伸。分別以200m、500m和1 000m緩沖半徑構(gòu)建了三維空間緩沖區(qū)模型，圖4是以200m、500m和1 000m緩沖半徑構(gòu)建的三維空間緩沖區(qū)模型、三維斷層模型以及地震在三維空間的分布。

圖4 三維空間分析結(jié)果示意圖Fig.4 Schematic diagram of 3D space analysis result.

在研究區(qū)廣泛分布的891個(gè)地震中，分別統(tǒng)計(jì)緩沖半徑取200m、500m和1 000m的情況下，目標(biāo)斷層緩沖區(qū)范圍內(nèi)地震的數(shù)量。以研究區(qū)內(nèi)東山山前斷裂為例，在緩沖半徑為200m的情況下，落入緩沖區(qū)的地震數(shù)為123個(gè)，約為13.8%;在緩沖半徑為500m的情況下，落入緩沖區(qū)的地震數(shù)為183個(gè)，約為20.5%;而在緩沖半徑為1 000m的情況下，891個(gè)地震中有264個(gè)落入緩沖區(qū)，約為29.6%。其余6條斷層也按照相同方式統(tǒng)計(jì)和分析，結(jié)果見表1。

表1 斷層與地震統(tǒng)計(jì)分析表Table 1 Statistical analysis table of faults and earthquakes

4 小結(jié)

在ArcGIS環(huán)境下提出了一種利用淺層地震探測(cè)資料構(gòu)建三維斷層模型的新方法，該方法通過對(duì)淺層地震探測(cè)深度剖面進(jìn)行“微分”，結(jié)合其他斷層資料實(shí)現(xiàn)了三維斷層精確還原。并利用ArcGIS強(qiáng)大的三維顯示功能，通過設(shè)計(jì)算法，采用AOVBA開發(fā)方式實(shí)現(xiàn)了三維斷層模型的構(gòu)建。在實(shí)現(xiàn)三維斷層模型的基礎(chǔ)上，在ArcGIS環(huán)境下嘗試性地提出了三維空間緩沖區(qū)模型的概念，通過設(shè)計(jì)算法和程序，實(shí)現(xiàn)了可按需選取緩沖半徑的斷層三維空間緩沖區(qū)模型。

基于三維空間緩沖區(qū)模型采用向量判別法實(shí)現(xiàn)了三維空間內(nèi)點(diǎn)和平面之間空間關(guān)系的算法，進(jìn)而實(shí)現(xiàn)了地震與活動(dòng)斷裂三維空間關(guān)系的分析方法。基于該算法對(duì)太原市區(qū)內(nèi)主要斷裂與地震的空間分布關(guān)系進(jìn)行了初步定量分析研究，并給出了定量的統(tǒng)計(jì)結(jié)果。由分析的實(shí)例結(jié)果可以看出，三維地質(zhì)模型能更真實(shí)、直觀地反映地震與斷裂的三維空間分布。

需要強(qiáng)調(diào)的是，由于沒有進(jìn)行小震精定位，本文最后的分析結(jié)果還不是特別精確。在此主要目的是提出地震與活動(dòng)斷裂三維空間關(guān)系的分析方法，為定量評(píng)價(jià)和研究活動(dòng)斷裂與地震之間的相關(guān)性提供了一種重要途徑。

鄧起東，張培震，冉勇康，等.2003.中國(guó)活動(dòng)構(gòu)造與地震活動(dòng)［J］.地學(xué)前緣，10(特刊):66—73.

DENG Qi-dong，ZHANG Pei-zhen，RAN Yong-kang，et al.2003.Active tectonics and earthquake activities in China［J］.Earth Science Frontiers，10(suppl):66—73(in Chinese).

劉峰，張家聲，黃雄南，等.2009.利用GIS方法研究南北地震帶和中央造山帶交會(huì)區(qū)活動(dòng)斷裂與地震的關(guān)系［J］.中國(guó)地震，25(4):394—404.

LIU Feng，ZHANG Jia-sheng，HUANG Xiong-nan，et al.2009.Research on the relationship between active faults and earthquakes in the junction area of the China North-South Seismic Belt and Central Orogenic Belt based on GIS［J］.Earthquake Research in China，25(4):394—404(in Chinese).

彭自正，趙愛平，胡翠娥，等.2002.基于GIS的江西活動(dòng)斷裂分布與地震活動(dòng)關(guān)系研究［J］.華南地震，22(4):9—18.

PENG Zi-zheng，ZHAO Ai-ping，HU Cui-e，et al.2002.Research on the relationship based on GIS between the distribution of the active faults and the seismicity in Jiangxi［J］.South China Journal of Seismology，22(4):9—18(in Chinese).

屈春燕，葉洪.2000.利用GIS分析活動(dòng)斷裂與地震的相關(guān)性［J］.地震研究，23(1):72—75.

QU Chun-yan，YE Hong.2000.Analysis of correlation between fault and earthquake using GIS ［J］.Journal of Seismological Research，23(1):72—75(in Chinese).

石樹中，宋立軍，楊紀(jì)林，等.1998.北天山活動(dòng)斷裂地理信息系統(tǒng)的設(shè)計(jì)與應(yīng)用［J］.內(nèi)陸地震，12(4):367—370.

SHI Shu-zhong，SONG Li-jun，YANG Ji-lin，et al.1998.Design and application of geographic information system of active faults in northern Tianshan Mountains［J］.Inland Earthquake，12(4):367—370(in Chinese).

鄔寧芬，周祖翼，勞秋元，等.2004.GIS支持下地震與活動(dòng)斷裂關(guān)系研究:以上海及鄰近海域?yàn)槔跩］.大地構(gòu)造與成礦學(xué)，28(3):248—253.

WU Ning-fen，ZHOU Zu-yi，LAO Qiu-yuan，et al.2004.GIS-based study on the relationship between earthquakes and active faults in Shanghai and its adjacent offshore region［J］.Geotectonica et Metallogenia，28(3):248—253(in Chinese).

吳立新，史文中，Christopher Gold.2003.3D GIS與3D GMS中的空間建模技術(shù)［J］.地理與地理信息科學(xué)，19(1):5—11.

WU Li-xin，SHI Wen-zhong，Christopher Gold.2003.Spatial modeling technologies for 3D GIS and 3D GMS［J］.Geography and Geo-Information Science，19(1):5—11(in Chinese).

姚鑫，王崇倡.2006.基于ArcScene的礦床空間三維模型的建立及可視化研究［J］.工程地質(zhì)計(jì)算機(jī)應(yīng)用，41(1):15—18.

YAO Xin，WANG Chong-chang.2006.Study of 3D model and visibility for mine-bed based on ArcScene［J］.Engineering Geology Computer Application，41(1):15—18(in Chinese).

于貴華，鄧起東，鄔倫，等.1996.利用GIS系統(tǒng)建立中國(guó)活動(dòng)斷裂信息咨詢分析系統(tǒng)［J］.地震地質(zhì)，18(2):156—160.

YU Gui-hua，DENG Qi-deng，WU Lun，et al.1996.Building Chinese active fault information system using GIS system［J］.Seismology and Geology，18(2):156—160(in Chinese).

張敏，李自紅，劉鴻福，等.2007.基于ArcGIS的太原斷陷盆地第四系三維地質(zhì)建模［J］.震災(zāi)防御技術(shù)，2(3):243—248.

ZHANG Min，LI Zi-hong，LIU Hong-fu，et al.2007.An ArcGIS based 3-D model of Taiyuan Graben Basin in Quaternary［J］.Technology for Earthquake Disaster Prevention，2(3):243—248(in Chinese).

中國(guó)地震局.2005.中國(guó)地震活動(dòng)斷層探測(cè)技術(shù)系統(tǒng)技術(shù)規(guī)程［S］.北京:地震出版社.

China Earthquake Administration.2005.Stipulation on Technical System for China Earthquake Active Fault Surveying［S］.Seismological Press，Beijing(in Chinese).