黃鎮(zhèn)逢

（廣東省有色金屬地質(zhì)局水文地質(zhì)隊(duì)，廣東 廣州 510800）

地質(zhì)災(zāi)害會(huì)帶來(lái)地質(zhì)環(huán)境的惡化，并使人類遭受巨大損失[1]。研究地質(zhì)災(zāi)害的形成根源時(shí)，區(qū)域性活動(dòng)斷裂帶是一種必須考慮的因素，這種斷裂帶不僅是地裂縫等重大災(zāi)害的誘發(fā)根源，也是控制常見(jiàn)地質(zhì)災(zāi)害實(shí)際分布情況的重要因素。與區(qū)域性活動(dòng)斷裂帶相關(guān)的三大類地質(zhì)災(zāi)害分別為活動(dòng)斷裂帶次生災(zāi)害、活動(dòng)斷裂帶緩慢活動(dòng)災(zāi)害以及活動(dòng)斷裂帶快速活動(dòng)災(zāi)害。隨著我國(guó)實(shí)施西部大開(kāi)發(fā)戰(zhàn)略，西部建設(shè)了大量高速公路，其規(guī)模可以說(shuō)是前所未有，也使西部城市的經(jīng)濟(jì)進(jìn)一步得到發(fā)展。然而該地區(qū)是頻發(fā)地質(zhì)災(zāi)害的區(qū)域，特別是城市區(qū)域性活動(dòng)斷裂帶十分發(fā)育，為公路的建設(shè)與城市的發(fā)展帶來(lái)極大阻礙與威脅，不僅帶來(lái)了直接損失，也造成了很多間接性損失[2]。與此同時(shí)，我國(guó)其他地區(qū)的城市區(qū)域性活動(dòng)斷裂帶也引發(fā)了很多地質(zhì)災(zāi)害，“5.12”汶川、“4.14”玉樹(shù)地震造成的巨大人員傷亡和財(cái)產(chǎn)損失讓人記憶猶新，而造成如此重大損失的根源就是構(gòu)筑物建在了活動(dòng)斷裂帶上，從而給予我們深刻警示，城市建設(shè)選址首先應(yīng)查明城市擬建區(qū)的地質(zhì)結(jié)構(gòu)，尤其是斷裂構(gòu)造，才能趨利避害，選擇避讓或提高防御級(jí)別，從而更有利于保障城市建設(shè)的安全。因此對(duì)城市區(qū)域性活動(dòng)斷裂帶實(shí)施危險(xiǎn)性評(píng)價(jià)是很有很大必要的。對(duì)于城市區(qū)域性活動(dòng)斷裂帶危險(xiǎn)性評(píng)價(jià)問(wèn)題，很多國(guó)家都進(jìn)行了大量研究，并取得了重要進(jìn)展。在研究中，新方法、新技術(shù)的應(yīng)用越來(lái)越深入和廣泛，包括3S技術(shù)等，這些方法的應(yīng)用，使評(píng)價(jià)結(jié)果變得更加全面。目前其研究方向主要朝著定量化發(fā)展，評(píng)價(jià)可靠性與評(píng)價(jià)精度都有待進(jìn)一步提高。因此綜合以往研究成果，提出一種基于地質(zhì)數(shù)據(jù)庫(kù)的城市區(qū)域性活動(dòng)斷裂帶危險(xiǎn)性評(píng)價(jià)方法。

1 城市區(qū)域性活動(dòng)斷裂帶危險(xiǎn)性評(píng)價(jià)

1.1 地質(zhì)數(shù)據(jù)庫(kù)構(gòu)建

構(gòu)建各地城市區(qū)域性活動(dòng)斷裂帶相關(guān)的城市活動(dòng)斷裂帶地質(zhì)數(shù)據(jù)庫(kù)，在數(shù)據(jù)庫(kù)的構(gòu)建中，首先需要對(duì)存在城市區(qū)域性活動(dòng)斷裂帶的各地的基礎(chǔ)地質(zhì)進(jìn)行調(diào)查，也就是實(shí)施地質(zhì)填圖。地質(zhì)填圖以穿越法為主，填圖路線的布置大致垂直于地層或斷層走向[3]。布設(shè)多條主干路線，調(diào)查人員統(tǒng)一參與主干路線地質(zhì)填圖，統(tǒng)一分層標(biāo)志，統(tǒng)一野外工作方法。對(duì)區(qū)內(nèi)主要地質(zhì)界線、重要標(biāo)志層等實(shí)施追索觀測(cè)。其中線距、點(diǎn)距根據(jù)野外實(shí)際情況布置，一般點(diǎn)距不大于線距。在第四系大面積分布區(qū)，以河流階地兩側(cè)或已開(kāi)挖基坑為重點(diǎn)布置填圖路線。在地質(zhì)填圖中，還要對(duì)人類工程活動(dòng)情況進(jìn)行同步記錄。

接著實(shí)施城市區(qū)域性活動(dòng)斷裂帶的野外調(diào)查。使用的物探方法是淺層地震多次迭加法，也就是采用6次覆蓋觀測(cè)系統(tǒng)進(jìn)行測(cè)量。在測(cè)量中根據(jù)試驗(yàn)成果確定偏移距。以時(shí)間剖面圖上反射波組的實(shí)際反射情況為依據(jù)推測(cè)斷裂構(gòu)造的具體位置，之后布設(shè)物探驗(yàn)證鉆孔進(jìn)行鉆探驗(yàn)證[4]，其中鉆孔孔位及標(biāo)高采用Garmin佳明errex201xGPS進(jìn)行現(xiàn)場(chǎng)測(cè)放，鉆孔采用劈芯法，采取一半巖芯作為樣品，樣長(zhǎng)5~8cm，對(duì)其實(shí)施單軸抗壓強(qiáng)度試驗(yàn)。

將調(diào)查獲取數(shù)據(jù)作為構(gòu)建城市活動(dòng)斷裂帶地質(zhì)數(shù)據(jù)庫(kù)的基礎(chǔ)數(shù)據(jù)，通過(guò)ArcSDE數(shù)據(jù)庫(kù)引擎將以上基礎(chǔ)數(shù)據(jù)以統(tǒng)一的數(shù)據(jù)模型形式存儲(chǔ)至PostgreSQL 關(guān)系數(shù)據(jù)庫(kù)中，實(shí)現(xiàn)數(shù)據(jù)的統(tǒng)一管理，完成城市活動(dòng)斷裂帶地質(zhì)數(shù)據(jù)庫(kù)的搭建。搭建的數(shù)據(jù)庫(kù)構(gòu)造具體如圖1所示。

其中數(shù)據(jù)預(yù)處理指的是數(shù)據(jù)質(zhì)量檢查，主要是對(duì)數(shù)據(jù)情況、時(shí)間精度以及數(shù)據(jù)完整性進(jìn)行檢查，經(jīng)過(guò)數(shù)據(jù)質(zhì)量檢查后再將數(shù)據(jù)存儲(chǔ)至PostgreSQL 關(guān)系數(shù)據(jù)庫(kù)中。

1.2 構(gòu)建評(píng)價(jià)模型

根據(jù)構(gòu)建的地質(zhì)數(shù)據(jù)庫(kù)對(duì)一級(jí)評(píng)價(jià)指標(biāo)進(jìn)行選取，包括已發(fā)生的地質(zhì)災(zāi)害與潛在地質(zhì)災(zāi)害情況、人類工程活動(dòng)情況、工程與水文地質(zhì)條件、巖性與地層構(gòu)造情況、地形地貌情況，并根據(jù)地質(zhì)數(shù)據(jù)庫(kù)數(shù)據(jù)對(duì)各一級(jí)評(píng)價(jià)指標(biāo)進(jìn)行詳細(xì)劃分，確定二級(jí)評(píng)價(jià)指標(biāo)，完成城市區(qū)域性活動(dòng)斷裂帶危險(xiǎn)性評(píng)價(jià)模型的構(gòu)建。構(gòu)建的評(píng)價(jià)模型具體如表1所示。

表1 構(gòu)建的評(píng)價(jià)模型

完成評(píng)價(jià)模型的構(gòu)建后，利用模型實(shí)施城市區(qū)域性活動(dòng)斷裂帶危險(xiǎn)性評(píng)價(jià)。

1.3 指標(biāo)權(quán)重的確定與評(píng)價(jià)的定量描述

計(jì)算各指標(biāo)權(quán)重，通過(guò)可拓學(xué)將城市區(qū)域性活動(dòng)斷裂帶危險(xiǎn)性評(píng)價(jià)的對(duì)象、指標(biāo)及指標(biāo)權(quán)值構(gòu)成一個(gè)整體，利用可拓集合的關(guān)聯(lián)函數(shù)值實(shí)現(xiàn)評(píng)價(jià)結(jié)果的定量描述[5]。

在各指標(biāo)權(quán)重的計(jì)算中，使用的方法為因子分析法。具體步驟如下：

（1）查閱文獻(xiàn)，對(duì)指標(biāo)相關(guān)數(shù)據(jù)進(jìn)行收集。

（2）對(duì)收集數(shù)據(jù)實(shí)施預(yù)處理，具體包括坐標(biāo)統(tǒng)一處理、缺失值處理、標(biāo)準(zhǔn)化處理。

在相關(guān)數(shù)據(jù)的坐標(biāo)統(tǒng)一處理中，使用的坐標(biāo)平移方法具體如下式：

式中：A——平移的直線距離；

——平移的坐標(biāo)距離；

X——平移的橫坐標(biāo)距離；

Y——平移的縱坐標(biāo)距離；

——該空間實(shí)體的原坐標(biāo)位置；

——該空間實(shí)體平移后的坐標(biāo)位置；

X2——平移中第二個(gè)坐標(biāo)中的橫坐標(biāo)；

Y2——平移中第二個(gè)坐標(biāo)中的縱坐標(biāo)；

X1——平移中第一個(gè)坐標(biāo)中的橫坐標(biāo)；

Y1——平移中第一個(gè)坐標(biāo)中的縱坐標(biāo)。

在缺失值處理中，使用高維映射法將缺失值補(bǔ)全。在標(biāo)準(zhǔn)化處理中，使用z-score 標(biāo)準(zhǔn)化處理方法將量綱消除。具體步驟如下：

①計(jì)算各指標(biāo)的標(biāo)準(zhǔn)差，用si來(lái)表示；

②計(jì)算各指標(biāo)的數(shù)學(xué)期望值，用xi來(lái)表示；

③實(shí)施標(biāo)準(zhǔn)化處理，具體如下式：

式中：xij——實(shí)際變量值；

zij——標(biāo)準(zhǔn)化處理以后的變量值。

④對(duì)調(diào)逆指標(biāo)前的正負(fù)號(hào)。

（3）對(duì)數(shù)據(jù)指標(biāo)之間的相關(guān)系數(shù)進(jìn)行計(jì)算，對(duì)選取的指標(biāo)是否能夠?qū)嵤┮蜃臃治鲞M(jìn)行判斷，能夠?qū)嵤┮蜃臃治龅那疤釛l件是指標(biāo)間的關(guān)聯(lián)性較強(qiáng)[6]。

（4）在關(guān)聯(lián)性較強(qiáng)的情況下，對(duì)公共因子進(jìn)行選取。在選取中首先對(duì)各因子的初始載荷矩陣進(jìn)行計(jì)算，并對(duì)相關(guān)系數(shù)構(gòu)成的矩陣的特征向量與特征根進(jìn)行計(jì)算，綜合兩種計(jì)算結(jié)果獲得因子的累積貢獻(xiàn)率。根據(jù)累積貢獻(xiàn)率對(duì)公共因子進(jìn)行選取，確定因子個(gè)數(shù)。

（5）計(jì)算公共因子得分，根據(jù)計(jì)算結(jié)果計(jì)算各因子的對(duì)應(yīng)方差貢獻(xiàn)率值，并將該值作為評(píng)價(jià)中的權(quán)重值。

接著利用可拓學(xué)思想定義評(píng)價(jià)中的待評(píng)物元、經(jīng)典域、節(jié)域。其中經(jīng)典域用Foij來(lái)表示，是指危險(xiǎn)性評(píng)價(jià)等級(jí)關(guān)于評(píng)價(jià)指標(biāo)的規(guī)定量化范圍。節(jié)域是指危險(xiǎn)性全體評(píng)價(jià)等級(jí)關(guān)于評(píng)價(jià)指標(biāo)的量值隸屬度化范圍，用Dqj來(lái)表示。將城市區(qū)域性活動(dòng)斷裂帶危險(xiǎn)性用物元O來(lái)表示，稱為待評(píng)物元。

通過(guò)下式計(jì)算城市區(qū)域性活動(dòng)斷裂帶危險(xiǎn)性這一評(píng)價(jià)單元q與危險(xiǎn)度等級(jí)之間的關(guān)聯(lián)函數(shù)值：

式中：cj——評(píng)價(jià)單元q關(guān)于區(qū)劃指標(biāo)的量值；

γ(.)——關(guān)聯(lián)閾值。

根據(jù)關(guān)聯(lián)函數(shù)值與指標(biāo)權(quán)值計(jì)算評(píng)價(jià)單元q關(guān)于危險(xiǎn)性等級(jí)i的實(shí)際關(guān)聯(lián)度值。根據(jù)計(jì)算的關(guān)聯(lián)度值對(duì)危險(xiǎn)性等級(jí)進(jìn)行計(jì)算，完成區(qū)域性活動(dòng)斷裂帶的危險(xiǎn)性評(píng)價(jià)。

2 評(píng)價(jià)案例測(cè)試

2.1 實(shí)驗(yàn)過(guò)程

對(duì)于設(shè)計(jì)的基于地質(zhì)數(shù)據(jù)庫(kù)的城市區(qū)域性活動(dòng)斷裂帶危險(xiǎn)性評(píng)價(jià)方法，通過(guò)其應(yīng)用對(duì)其評(píng)價(jià)性能進(jìn)行測(cè)試。將其應(yīng)用于城市區(qū)域性活動(dòng)斷裂帶比較發(fā)育的某地區(qū)中，利用該方法對(duì)實(shí)驗(yàn)地區(qū)的活動(dòng)斷裂帶進(jìn)行危險(xiǎn)性評(píng)價(jià)。實(shí)驗(yàn)地區(qū)的崩塌、滑坡主要分布于工作區(qū)北東側(cè)的低山、丘陵區(qū)，因此可以判定其區(qū)域性活動(dòng)斷裂帶主要發(fā)育于北東側(cè)。由于該地區(qū)已經(jīng)遭受過(guò)多次地質(zhì)災(zāi)害，因此對(duì)其區(qū)域性活動(dòng)斷裂帶實(shí)施危險(xiǎn)性評(píng)價(jià)是很有必要的。

首先構(gòu)建該地區(qū)的城市活動(dòng)斷裂帶地質(zhì)數(shù)據(jù)庫(kù)，分別實(shí)施該地區(qū)的基礎(chǔ)地質(zhì)調(diào)查與區(qū)域性活動(dòng)斷裂帶的野外調(diào)查，將調(diào)查數(shù)據(jù)作為基礎(chǔ)數(shù)據(jù)完成地質(zhì)數(shù)據(jù)庫(kù)的構(gòu)建。接著構(gòu)建評(píng)價(jià)模型，通過(guò)因子分析法計(jì)算實(shí)驗(yàn)地區(qū)的指標(biāo)權(quán)重，計(jì)算結(jié)果如表2所示。

表2 實(shí)驗(yàn)地區(qū)的指標(biāo)權(quán)重計(jì)算結(jié)果

然后結(jié)合表2 結(jié)果，依據(jù)可拓學(xué)實(shí)施危險(xiǎn)性評(píng)價(jià)，評(píng)價(jià)結(jié)果如表3所示。

表3 危險(xiǎn)性評(píng)價(jià)結(jié)果

表3數(shù)據(jù)說(shuō)明，通過(guò)設(shè)計(jì)方法能夠?qū)崿F(xiàn)城市區(qū)域性活動(dòng)斷裂帶的危險(xiǎn)性評(píng)價(jià)。完成危險(xiǎn)性評(píng)價(jià)后，對(duì)評(píng)價(jià)結(jié)果的準(zhǔn)確性進(jìn)行測(cè)試。

2.2 評(píng)價(jià)準(zhǔn)確性測(cè)試結(jié)果

根據(jù)專家意見(jiàn)對(duì)評(píng)價(jià)結(jié)果的準(zhǔn)確性進(jìn)行測(cè)試，測(cè)試結(jié)果如圖2所示。

圖2 評(píng)價(jià)準(zhǔn)確性

根據(jù)圖2測(cè)試結(jié)果，設(shè)計(jì)方法對(duì)于六個(gè)區(qū)域性活動(dòng)斷裂帶的危險(xiǎn)性評(píng)價(jià)準(zhǔn)確性均較高，六種評(píng)價(jià)結(jié)果的準(zhǔn)確性均高于96%，證明了設(shè)計(jì)方法有著出色的評(píng)價(jià)性能。

3 結(jié)束語(yǔ)

城市區(qū)域性活動(dòng)斷裂帶是一種需要重視的構(gòu)造因素，很容易引發(fā)各種地質(zhì)災(zāi)害，從而造成極大損失。設(shè)計(jì)了一種基于地質(zhì)數(shù)據(jù)庫(kù)的城市區(qū)域性活動(dòng)斷裂帶危險(xiǎn)性評(píng)價(jià)方法，既實(shí)現(xiàn)了地區(qū)城市活動(dòng)斷裂帶地質(zhì)數(shù)據(jù)庫(kù)的構(gòu)建，又實(shí)現(xiàn)了活動(dòng)斷裂帶的準(zhǔn)確危險(xiǎn)性評(píng)價(jià)，取得了一定研究成果。今后將會(huì)對(duì)該方法進(jìn)行進(jìn)一步實(shí)踐與應(yīng)用。