陳祚偉

（江西省地質(zhì)局第三地質(zhì)大隊(duì) 江西九江 332000）

泰和縣地處吉泰盆地中部，縣域地勢總體為東部和西部高、中部低，地形以低山、丘陵為主，面積1789.07km2，占全縣總面積的67.12%，縣域內(nèi)地層巖性復(fù)雜，巖土體工程力學(xué)性質(zhì)較差，斷層褶皺等地質(zhì)構(gòu)造較發(fā)育，區(qū)內(nèi)降雨集中，切坡筑路、切坡建房、礦山開采、水利水電建設(shè)等人類工程活動(dòng)強(qiáng)烈［1］。受上述因素的綜合影響，導(dǎo)致區(qū)內(nèi)環(huán)境地質(zhì)問題較嚴(yán)重。地質(zhì)災(zāi)害類型主要為滑坡、崩塌、地面塌陷，存在滑坡隱患和崩塌隱患，已成為贛中地質(zhì)災(zāi)害高易發(fā)縣之一，每年因地質(zhì)災(zāi)害造成的經(jīng)濟(jì)損失較嚴(yán)重。

1 評價(jià)原則與方法

1.1 評價(jià)原則

工作區(qū)內(nèi)自然地質(zhì)環(huán)境條件復(fù)雜，切坡建房、工程建設(shè)等人類活動(dòng)差異大，導(dǎo)致地質(zhì)災(zāi)害造成的影響和危害方式存在一定差異性。為此，在充分認(rèn)識區(qū)內(nèi)孕災(zāi)地質(zhì)環(huán)境條件和地質(zhì)災(zāi)害形成機(jī)理的基礎(chǔ)上，進(jìn)行地質(zhì)災(zāi)害易發(fā)性評價(jià)時(shí)，需首先遵循以反映該區(qū)地質(zhì)災(zāi)害分布規(guī)律與成因機(jī)理為主，并且在突出重要誘發(fā)因素的原則下，圍繞“以人為本”的基本原則，同時(shí)，應(yīng)兼顧為地質(zhì)災(zāi)害防災(zāi)減災(zāi)服務(wù)的原則，以及對全區(qū)進(jìn)行綜合性評價(jià)與分析等5 個(gè)評價(jià)原則［2-3］。通過對區(qū)內(nèi)已查明的災(zāi)害與之災(zāi)害成因進(jìn)行關(guān)聯(lián)性分析，總結(jié)地質(zhì)災(zāi)害發(fā)育特征及分布規(guī)律，選擇適當(dāng)?shù)亩吭袨?zāi)因子指標(biāo)，體現(xiàn)地質(zhì)災(zāi)害孕災(zāi)地質(zhì)環(huán)境條件，再采取分析類比法，定性分析區(qū)內(nèi)地質(zhì)災(zāi)害易發(fā)程度，最終采用定性、定量結(jié)合進(jìn)行地質(zhì)災(zāi)害易發(fā)性評價(jià)。

1.2 評價(jià)方法

由于地質(zhì)災(zāi)害的孕災(zāi)因素很多，既包括孕災(zāi)體的地質(zhì)作用，也包括外部降雨等誘發(fā)作用，還會受到人類工程活動(dòng)的影響，因此地質(zhì)災(zāi)害易發(fā)程度分區(qū)評價(jià)是一個(gè)復(fù)雜的多元評價(jià)系統(tǒng)。本次運(yùn)用地質(zhì)災(zāi)害綜合危險(xiǎn)性指數(shù)法對地質(zhì)災(zāi)害易發(fā)性評價(jià)，評價(jià)流程如圖1所示，總體思路如下。

圖1 地質(zhì)災(zāi)害易發(fā)評價(jià)流程圖

（1）首先，確定控制和影響地質(zhì)災(zāi)害發(fā)生的主要孕災(zāi)因素，以專家定性打分的評價(jià)方法，確定各種孕災(zāi)因素的權(quán)重。

（2）將研究區(qū)整體進(jìn)行單元網(wǎng)格化，然后將每個(gè)單元內(nèi)地質(zhì)災(zāi)害影響因素進(jìn)行分類定性評價(jià)，根據(jù)所等的定性評價(jià)結(jié)果，采用重分類對每個(gè)單元進(jìn)行量化。

（3）將量化后的每個(gè)單元格在GIS 平臺上進(jìn)行各類地質(zhì)災(zāi)害信息量疊加分析。同時(shí)，根據(jù)個(gè)別影響因子進(jìn)行合理的修正，最終得出較為合理的易發(fā)性評價(jià)分區(qū)結(jié)果。

2 評價(jià)模型的建立與影響因子選取

2.1 評價(jià)單元網(wǎng)格劃分

根據(jù)泰和縣的孕災(zāi)地質(zhì)環(huán)境條件與地質(zhì)災(zāi)害分布發(fā)育特征，采用1cm×1cm（實(shí)際大小500m×500m）規(guī)則化單元方格進(jìn)行網(wǎng)格化地質(zhì)災(zāi)害易發(fā)性評價(jià)，全縣共劃分出11 187個(gè)計(jì)算單元格，其中，有災(zāi)害點(diǎn)的單元格數(shù)為1289個(gè)。

2.2 分析計(jì)算方法

地質(zhì)災(zāi)害綜合危險(xiǎn)性指數(shù)的計(jì)算方法：

式中：Z表示地質(zhì)災(zāi)害綜合指數(shù)；Zq表示孕災(zāi)條件強(qiáng)度指數(shù)；r1表示孕災(zāi)條件強(qiáng)度權(quán)值，取值0.65；Zx表示現(xiàn)狀地質(zhì)災(zāi)害發(fā)育強(qiáng)度指數(shù)；r2表示現(xiàn)狀地質(zhì)災(zāi)害發(fā)育強(qiáng)度權(quán)值，取值0.35。

2.3 孕災(zāi)地質(zhì)環(huán)境條件強(qiáng)度指數(shù)計(jì)算

地質(zhì)災(zāi)害孕災(zāi)地質(zhì)環(huán)境條件強(qiáng)度指數(shù)（Zq）按以下公式計(jì)算：

式中，Ti分別為控制評價(jià)單元地質(zhì)災(zāi)害孕災(zāi)地質(zhì)條件（D）、地形地貌條件（X）、人類工程活動(dòng)條件（R）的強(qiáng)度表征值，Ai分別為各孕災(zāi)地質(zhì)環(huán)境條件的權(quán)重值。

2.3.1 孕災(zāi)因素的確定

根據(jù)泰和縣地質(zhì)災(zāi)害的分布特征與發(fā)育規(guī)律，確定其基本孕災(zāi)條件因子有微地貌地形坡度、地貌類型、工程地質(zhì)巖組、地質(zhì)構(gòu)造等；誘發(fā)因素有大氣降雨、地震、孕災(zāi)的人類工程活動(dòng)等［4-5］。因全區(qū)位于地震烈度小于Ⅵ度區(qū)內(nèi)，大氣降雨則具有隨機(jī)性、變動(dòng)性及不確定性等因素，故本次評價(jià)將降雨、地震兩項(xiàng)因素予以剔除。

2.3.2 評價(jià)孕災(zāi)因子權(quán)重的確定

Ai分別為各孕災(zāi)條件的權(quán)重值，孕災(zāi)因子有工程地質(zhì)條件（G）、地形地貌條件（D）和人類工程活動(dòng)條件（R），然后，采用層次分析法，確定各個(gè)孕災(zāi)因子的權(quán)重。

首先，建立準(zhǔn)則層各要素，即工程地質(zhì)條件、地形地貌條件、人類工程活動(dòng)條件對目標(biāo)層的判別矩陣A。

計(jì)算該矩陣中各要素的特征向量并進(jìn)行歸一化處理后，得到各評價(jià)因子對孕災(zāi)地質(zhì)條件強(qiáng)度指數(shù)的權(quán)重向量，通過計(jì)算，得到權(quán)重向量Ab值，與其對應(yīng)的CI、CR值，再通過一致性檢驗(yàn)，得出判斷矩陣A為有效。

在構(gòu)建A矩陣后，再分別進(jìn)行構(gòu)建底層要素矩陣B，各要素要求一一對應(yīng)并歸一化處理后，分別得出底層要素權(quán)重向量BG、D、R，因人類工程活動(dòng)只進(jìn)行范圍影響分析，單獨(dú)為1個(gè)因素。

在得到權(quán)重A、B后，自上而下將單準(zhǔn)則下的權(quán)重進(jìn)行計(jì)算，得到底層中工程地質(zhì)巖組類型、地質(zhì)構(gòu)造等各孕災(zāi)因子對應(yīng)的強(qiáng)度指數(shù)權(quán)重Ai。

2.3.3 致災(zāi)因素的指數(shù)化

對每一個(gè)含有已發(fā)生地質(zhì)災(zāi)害的評價(jià)單元網(wǎng)格內(nèi)的各類致災(zāi)地質(zhì)環(huán)境條件因素，按照各自的影響程度賦以與之相匹配的權(quán)重，并對每類致災(zāi)因素進(jìn)行分級賦相應(yīng)值，將致災(zāi)因素權(quán)重值與相應(yīng)的分級值作乘法計(jì)算，從而得出各孕災(zāi)因素相應(yīng)對應(yīng)的指數(shù)值，再利用ARCGIS，將各因素的指數(shù)值進(jìn)行疊加分析，得出致災(zāi)強(qiáng)度指數(shù)，從而判斷出評價(jià)單元或該區(qū)域地質(zhì)災(zāi)害發(fā)生的易發(fā)性［6］。各孕災(zāi)因子權(quán)重及分級評分值見表1。

表1 各評價(jià)因子權(quán)重與分級賦值一覽表

采用ARCGIS 軟件分區(qū)統(tǒng)計(jì)工具對網(wǎng)格區(qū)的各評價(jià)因子進(jìn)行屬性賦值，然后，通過對屬性要素加權(quán)疊加計(jì)算，得出各網(wǎng)格區(qū)的Zq值，Zq值分布范圍為0.9553～6.3436，Zq值主要集中分布在1.5～4.5。

2.4 現(xiàn)狀地質(zhì)災(zāi)害強(qiáng)度指數(shù)計(jì)算

現(xiàn)狀地質(zhì)災(zāi)害強(qiáng)度指數(shù)（Zx）可以用已發(fā)生地質(zhì)災(zāi)害點(diǎn)歸一化處理后的點(diǎn)密度、面密度及體積密度，利用ARCGIS密度分析求得。

首先，將災(zāi)害點(diǎn)進(jìn)行歸一化處理并得出災(zāi)害點(diǎn)數(shù)密度的系數(shù)，即每個(gè)單元網(wǎng)格內(nèi)地質(zhì)災(zāi)害個(gè)數(shù)與整個(gè)評價(jià)區(qū)內(nèi)的災(zāi)害點(diǎn)個(gè)數(shù)之間的比值，同理求出災(zāi)害點(diǎn)面密度系數(shù)與體積密度系數(shù)。經(jīng)密度分析計(jì)算，得出網(wǎng)格內(nèi)有值的為1289個(gè)網(wǎng)格，占網(wǎng)格總數(shù)的11.52%，Zx值分布在0～180，單元格內(nèi)Zx值主要集中分布在0～20。

2.5 地質(zhì)災(zāi)害綜合危險(xiǎn)性指數(shù)

根據(jù)各單元的工程地質(zhì)條件、地形地貌及人類工程活動(dòng)等孕災(zāi)地質(zhì)環(huán)境條件，利用ARCGIS 中的空間分析功能，取得各評價(jià)單元網(wǎng)格內(nèi)的孕災(zāi)地質(zhì)條件強(qiáng)度指數(shù)與現(xiàn)狀地質(zhì)災(zāi)害致災(zāi)強(qiáng)度指數(shù)［7］，根據(jù)公式Z=Zq·0.65+Zx·0.35，將各強(qiáng)度指數(shù)進(jìn)行疊加換算，從而獲得各評價(jià)單元網(wǎng)格內(nèi)地質(zhì)災(zāi)害綜合危險(xiǎn)性指數(shù)（Z值），Z=0.0004～179.7895（統(tǒng)計(jì)時(shí)排除了Z＞15的值，統(tǒng)計(jì)網(wǎng)格為11 058個(gè)），最終得出Z值集中分布在0.55～3.98 之間，一共有9763 個(gè)單元網(wǎng)格，占統(tǒng)計(jì)網(wǎng)格數(shù)的88.29%。

依據(jù)地質(zhì)災(zāi)害綜合危險(xiǎn)性指數(shù)Z值分布特征及曲率圖，取主要拐點(diǎn)相對于的值，再根據(jù)各取值范圍一一對應(yīng)地質(zhì)災(zāi)害易發(fā)性等級，從而得出各地質(zhì)災(zāi)害易發(fā)性等級區(qū)間值。經(jīng)過曲率分析，主要拐點(diǎn)分別為：＜0.9424、0.9424～1.3451、1.3451～1.8702、＞1.8702，將Z值劃分為Ⅳ級。Z＞1.8702 為地質(zhì)災(zāi)害高易發(fā)區(qū)，有2945 個(gè)單元網(wǎng)格，占整個(gè)單元網(wǎng)格數(shù)的26.33%；Z值為1.3451～1.8702 為地質(zhì)災(zāi)害中易發(fā)區(qū)，有單元網(wǎng)格4379個(gè)，占整個(gè)單元網(wǎng)格數(shù)的39.14%；Z值為0.9424～1.3451 為地質(zhì)災(zāi)害低易發(fā)區(qū)，有單元網(wǎng)格1944 個(gè)，占整個(gè)單元網(wǎng)格數(shù)的17.38%；Z＜0.9424為地質(zhì)災(zāi)害不易發(fā)區(qū)，共有1790 個(gè)單元網(wǎng)格，占單元網(wǎng)格總數(shù)的17.15%，如圖2所示。

圖2 泰和縣易發(fā)性指數(shù)值分布單元網(wǎng)格圖

3 地質(zhì)災(zāi)害易發(fā)性評價(jià)結(jié)果

地質(zhì)災(zāi)害綜合危險(xiǎn)指數(shù)Z值得出的易發(fā)分區(qū)結(jié)果在一定程度上較離散，為了更好地反映泰和縣地質(zhì)災(zāi)害在區(qū)域上的一致性，將全縣劃分為145個(gè)流域，再根據(jù)流域內(nèi)各類方格的占比，結(jié)合定性分析，劃分出地質(zhì)災(zāi)害易發(fā)性分區(qū)。將全縣劃分為20 個(gè)地質(zhì)災(zāi)害易發(fā)性分區(qū)（見圖3），其中，有6個(gè)地質(zhì)災(zāi)害高易發(fā)區(qū)，面積1130.16km2，占全縣面積的42.48%，區(qū)內(nèi)有已發(fā)生地質(zhì)災(zāi)害點(diǎn)487處，占全縣災(zāi)害點(diǎn)總數(shù)的80.23%；中易發(fā)區(qū)6個(gè)，面積452.73km2，占全縣面積的17.02%，區(qū)內(nèi)有已發(fā)生地質(zhì)災(zāi)害點(diǎn)49處，占全縣災(zāi)害點(diǎn)總數(shù)的8.07%；低易發(fā)區(qū)5個(gè)，面積628.84km2，占全縣面積的23.64%，區(qū)內(nèi)有已發(fā)生地質(zhì)災(zāi)害點(diǎn)30處，占全縣災(zāi)害點(diǎn)總數(shù)的4.94%；不易發(fā)區(qū)3 個(gè)，面積448.57km2，占全縣面積的16.86%，區(qū)內(nèi)有已發(fā)生地質(zhì)災(zāi)害點(diǎn)41 處，占全縣災(zāi)害點(diǎn)總數(shù)的6.75%。

圖3 泰和縣地質(zhì)災(zāi)害易發(fā)性評價(jià)結(jié)果圖

4 結(jié)論與建議

（1）本次評價(jià)工作主要利用收集前人資料與縣域1/50 000 地質(zhì)災(zāi)害調(diào)查資料，在分析該區(qū)地質(zhì)災(zāi)害發(fā)育強(qiáng)度與分布規(guī)律的基礎(chǔ)上，選取地形坡度、地形地貌、地質(zhì)構(gòu)造、工程地質(zhì)巖組、人類工程活動(dòng)等5 項(xiàng)孕災(zāi)因子，采取層次分析的綜合指數(shù)法進(jìn)行的地質(zhì)災(zāi)害易發(fā)程度評價(jià)。

（2）本次評價(jià)結(jié)果與該區(qū)地質(zhì)災(zāi)害實(shí)際發(fā)生及分布情況基本吻合，反映出層次分析的綜合指數(shù)法適合本區(qū)地質(zhì)災(zāi)害易發(fā)性評價(jià)，其結(jié)果合理可靠，可作為地質(zhì)災(zāi)害防治與工程建設(shè)的依據(jù)。

（3）建議泰和縣在地質(zhì)災(zāi)害中、高易發(fā)區(qū)進(jìn)行工程建設(shè)時(shí)，必須對場地或線路進(jìn)行更為詳細(xì)的地質(zhì)災(zāi)害危險(xiǎn)性評估工作，以保障人民生命財(cái)產(chǎn)不受地質(zhì)災(zāi)害破壞，避免損失。