基于AHP和書ArcGlS技術(shù)的巖溶塌陷風(fēng)險(xiǎn)評(píng)價(jià)——以遵義永樂鎮(zhèn)為例

潘宗源，賈 龍，劉寶臣

（1.中國地質(zhì)科學(xué)院巖溶地質(zhì)研究所/巖溶塌陷防治重點(diǎn)實(shí)驗(yàn)室，廣西　　桂林 541004；2.桂林理工大學(xué)，廣西　　桂林 541004）

基于AHP和書ArcGlS技術(shù)的巖溶塌陷風(fēng)險(xiǎn)評(píng)價(jià)——以遵義永樂鎮(zhèn)為例

潘宗源1，賈 龍1，劉寶臣2

（1.中國地質(zhì)科學(xué)院巖溶地質(zhì)研究所/巖溶塌陷防治重點(diǎn)實(shí)驗(yàn)室，廣西桂林 541004；2.桂林理工大學(xué)，廣西桂林 541004）

對(duì)遵義永樂鎮(zhèn)巖溶塌陷區(qū)域風(fēng)險(xiǎn)因素進(jìn)行了識(shí)別和確定，選取巖溶條件、覆蓋層條件、構(gòu)造條件、地下水條件等4個(gè)準(zhǔn)則層，共13個(gè)風(fēng)險(xiǎn)因子來建立層次結(jié)構(gòu)模型。結(jié)合當(dāng)?shù)厮牡刭|(zhì)背景，運(yùn)用層次分析法（AHP）計(jì)算出風(fēng)險(xiǎn)因子權(quán)重；確定各區(qū)域巖溶塌陷風(fēng)險(xiǎn)等級(jí)，借助ArcGIS技術(shù)對(duì)研究區(qū)進(jìn)行巖溶塌陷風(fēng)險(xiǎn)評(píng)價(jià)。結(jié)果表明：塌陷極危險(xiǎn)區(qū)呈北北西向、不規(guī)則條帶狀位于研究區(qū)中部；塌陷危險(xiǎn)區(qū)和次危險(xiǎn)區(qū)依次位于極危險(xiǎn)區(qū)的兩側(cè)；最外側(cè)的為次穩(wěn)定區(qū)；最南邊及東邊局地為穩(wěn)定區(qū)。該結(jié)果與近年來發(fā)生巖溶塌陷的地理位置和塌陷危害程度相吻合。

巖溶塌陷；風(fēng)險(xiǎn)因素；層次結(jié)構(gòu)模型；ArcGIS；風(fēng)險(xiǎn)評(píng)價(jià)；遵義

我國可溶巖覆蓋面廣泛，尤其是西南地區(qū)為主要發(fā)育區(qū)。在可溶巖發(fā)育及一定蓋層厚度的基礎(chǔ)地質(zhì)背景下，附加外動(dòng)力影響條件下（如地下水劇烈波動(dòng)、大幅度抽排地下水、強(qiáng)降雨或工程施工等）較容易發(fā)生巖溶塌陷。近年來，巖溶塌陷在巖溶覆蓋區(qū)時(shí)有發(fā)生，其在時(shí)空領(lǐng)域上的隱蔽性和突發(fā)性，給塌陷區(qū)人民生命財(cái)產(chǎn)帶來嚴(yán)重?fù)p失，巖溶塌陷是我國主要地質(zhì)災(zāi)害之一。

我國學(xué)者自20世紀(jì)60年代就開始對(duì)巖溶塌陷進(jìn)行分析研究，并取得了非常豐碩的成果。在碳酸鹽巖覆蓋區(qū)，文獻(xiàn)［1-6］將地質(zhì)調(diào)查與試驗(yàn)手段相結(jié)合，對(duì)巖溶發(fā)育特征、塌陷形成機(jī)制、塌陷分布規(guī)律、塌陷誘發(fā)因素等進(jìn)行了定量的研究分析；在巖溶塌陷危險(xiǎn)性區(qū)域分區(qū)評(píng)價(jià)和預(yù)測(cè)方面，文獻(xiàn)［7-15］運(yùn)用不同的評(píng)價(jià)體系對(duì)巖溶塌陷地區(qū)進(jìn)行了時(shí)間和空間上的定量化預(yù)測(cè)分析和初步評(píng)價(jià)；在巖溶塌陷危險(xiǎn)區(qū)監(jiān)測(cè)方面，文獻(xiàn)［16-18］通過預(yù)埋光纖的手段對(duì)隱伏巖溶土洞及覆蓋層厚度變化進(jìn)行監(jiān)測(cè)，以達(dá)到對(duì)巖溶塌陷進(jìn)行預(yù)警的目的；在對(duì)巖溶發(fā)育特征勘測(cè)方面，文獻(xiàn)［19-22］合理地運(yùn)用地球物理勘探方法對(duì)地下巖溶發(fā)育特征進(jìn)行探測(cè)分析，證明物理勘探能有效地探明巖溶發(fā)育位置及特征。

本文通過總結(jié)前人的工作成果，以貴州省遵義市永樂鎮(zhèn)為研究對(duì)象，分析影響巖溶塌陷的風(fēng)險(xiǎn)因子，將層次分析法（AHP）與ArcGIS相結(jié)合，對(duì)研究區(qū)進(jìn)行巖溶塌陷風(fēng)險(xiǎn)評(píng)價(jià)。

1　研究方法

本文綜合運(yùn)用AHP和ArcGIS技術(shù)對(duì)巖溶發(fā)育區(qū)域進(jìn)行塌陷風(fēng)險(xiǎn)性評(píng)價(jià)分析。層次分析法（AHP）是將對(duì)總目標(biāo)構(gòu)成影響性的基本因子，根據(jù)因子間關(guān)聯(lián)性和隸屬關(guān)系分為多層次結(jié)構(gòu)，以兩兩因子相比較來確定其重要程度，最終通過層層比較和專家經(jīng)驗(yàn)來確定基本因子的權(quán)重。該方法可以將對(duì)總目標(biāo)（目標(biāo)層）有影響性的不確定性基本因子都囊括在內(nèi)，巖溶塌陷本身受控的因子就是復(fù)雜且多變的，大體上可分為：土層厚度、土層結(jié)構(gòu)、巖溶發(fā)育程度、距斷裂構(gòu)造距離、巖層結(jié)構(gòu)、地下水波動(dòng)情況和人類影響因素等7方面［23］，相互間都有著緊密的聯(lián)系，難以用定量的統(tǒng)計(jì)方法對(duì)巖溶塌陷危險(xiǎn)區(qū)進(jìn)行評(píng)價(jià)。而ArcGIS在地理信息統(tǒng)計(jì)分析和應(yīng)對(duì)圖層拆解、疊加方面的具有優(yōu)勢(shì)，結(jié)合AHP能對(duì)巖溶塌陷風(fēng)險(xiǎn)性作較全面和直觀的預(yù)測(cè)評(píng)價(jià)。

2　研究區(qū)概況

2.1環(huán)境地質(zhì)條件

2.1.1地理氣候條件 研究區(qū)位于貴州省遵義市東部永樂鎮(zhèn)，茅官河西岸斜坡上。其中部一帶地形平緩，地形坡度為5°～20°；北西山脊一帶地形抬升，坡度為20°～35°；南東河岸地帶地形陡峭，多呈陡崖，高20～37 m。該地氣候上屬亞熱帶季風(fēng)氣候，年均降雨量為1 200 mm，年平均相對(duì)濕度在82%左右。

2.1.2地質(zhì)構(gòu)造 區(qū)域構(gòu)造上位于黔北臺(tái)隆遵義斷拱風(fēng)岡北北東向構(gòu)造變形區(qū)，地處川黔南北向構(gòu)造帶與風(fēng)岡北北東向構(gòu)造帶交匯部位，蝦子場復(fù)向斜東翼。區(qū)內(nèi)地層為單斜構(gòu)造，巖層總體向南西傾斜，傾角較平緩（6°～20°）。在研究區(qū)內(nèi)共有斷層9條，其中：北東向和北北東向兩組斷裂構(gòu)造最為發(fā)育（6條），處于大斷裂構(gòu)造南西端尖滅部位；另有近東西向斷層3條。

2.1.3地層分布及特征 區(qū)內(nèi)出露地層為志留系、三疊系和二疊系，以三疊系和二疊系為主，缺失志留系上統(tǒng)、泥盆系與石炭系。上覆蓋土層厚度為0～5 m，耕植土層較薄，主要為黃褐色粘土、粉質(zhì)粘土，中夾雜碎石。覆蓋層下為二疊系長興組（P3c）灰色含燧石生物屑泥晶灰?guī)r，下部有厚10 m左右的含碳泥質(zhì)灰?guī)r夾頁巖，頂部為薄層泥灰?guī)r及鈣質(zhì)粘土巖，厚約38 m。

2.1.4地下水補(bǔ)徑排 大氣降水是地下水的主要補(bǔ)給來源，由于區(qū)域內(nèi)碳酸鹽巖地層巖溶發(fā)育，溶蝕裂隙、落水洞等較多，降水一方面通過裂隙入滲補(bǔ)給地下水，另一方面匯成地表水從落水洞流入直接補(bǔ)給地下水。區(qū)內(nèi)為單斜地層，地形北高南低，茅官河從研究區(qū)北面由西向東，再轉(zhuǎn)向南，呈圓弧形繞研究區(qū)邊緣而過，為該區(qū)地下水的排匯基準(zhǔn)面。由金竹山起，過田月包至老鷹巖山脊一帶為地下水分水嶺，出露二疊系龍?zhí)督M（P3l）薄—中層狀粉砂巖、粉砂質(zhì)粘土巖隔水層，呈東西向延伸，將研究區(qū)分成兩個(gè)水文系統(tǒng)單元：Ⅰ單元，即隔水層以北區(qū)域內(nèi)，地下水由南向北匯流；Ⅱ單元，隔水層以南區(qū)域內(nèi)，地下水呈“扇”形向南東方向匯流，于茅官河岸兩個(gè)泉水點(diǎn)集中排泄。

2.2巖溶發(fā)育特征

研究區(qū)內(nèi)上二疊統(tǒng)長興組（P3c）、下三疊統(tǒng)夜郎組第二段（T1y2）巖溶發(fā)育較為強(qiáng)烈，其巖性分別為薄—中厚層狀灰?guī)r、含燧石結(jié)核灰?guī)r和中厚層含燧石灰?guī)r，地理位置上分布于區(qū)域北面山頂?shù)奶镌掳ⅫS泥頂，區(qū)內(nèi)及其周邊地帶，以及南西白土壩、新田壩、華家灣一帶。研究區(qū)巖溶發(fā)育特征如圖1所示。地表上，溶潭、落水洞（消水洞）發(fā)育：溶潭直徑為29～35 m，深度不詳，常年有積水，消水能力差，大雨后水位上升約5 m，之后水位緩慢下降，但一般不會(huì)干枯；落水洞直徑一般為0.2～0.5 m，多為不規(guī)則形，可見深度0.5～1.0 m；落水洞消水量最大1.5 L/s，主要來源于洼地地表水補(bǔ)給。地下則發(fā)育較多巖溶裂隙和溶洞，溶洞埋深一般在13～20 m，高度為0.5～1.8 m，多為半充填或無充填，充填物為黃色軟塑狀黏土。底部基巖溶蝕發(fā)育，多見溶蝕裂隙和孔洞，裂隙有黃色軟塑狀黏土充填。

圖1　研究區(qū)巖溶發(fā)育示意圖Fig.1 Karst development characteristics in research area

3　巖溶塌陷風(fēng)險(xiǎn)定量評(píng)價(jià)

3.1巖溶塌陷層次結(jié)構(gòu)模型

巖溶塌陷是指在巖溶地區(qū)，下部可溶巖層中的溶蝕洞穴或上覆土層中的土洞，在外在營力或自然因素變化時(shí)發(fā)生地面下陷現(xiàn)象。之前有不少學(xué)者選取不同影響因素，對(duì)巖溶塌陷的必要因素和誘發(fā)因素進(jìn)行了分析［3-8］。通過對(duì)本研究區(qū)水文地質(zhì)條件的分析，總結(jié)區(qū)內(nèi)巖溶發(fā)育情況如下：（1）區(qū)內(nèi)地形較緩，地表巖溶發(fā)育特征明顯，巖溶塌陷、溶潭、落水洞、豎井、天窗等巖溶形態(tài)點(diǎn)較多，是降雨匯入地下河的良好運(yùn)輸通道，使得在遭遇極端氣候條件下地下水位容易劇烈波動(dòng)，對(duì)較為脆弱的土層頂板產(chǎn)生壓-吸效應(yīng)；（2）現(xiàn)場踏勘資料顯示，研究區(qū)內(nèi)共有9條斷層，對(duì)地層切割作用強(qiáng)烈，裂隙和溶蝕管道較發(fā)育；（3）根據(jù)地質(zhì)鉆探資料，該區(qū)覆蓋層厚度較薄，土洞埋深淺，因而容易向上擴(kuò)展最終導(dǎo)致塌陷；（4）區(qū)內(nèi)地下水位為10～17 m，該區(qū)段亦為裂隙和溶蝕洞穴發(fā)育強(qiáng)烈區(qū)間，在大量抽排水和暴雨時(shí)地下水動(dòng)力會(huì)急劇增強(qiáng)，地下水位快速抬升或下降，加大巖溶塌陷的風(fēng)險(xiǎn)。

因此，本文選擇4項(xiàng)共13方面的風(fēng)險(xiǎn)因素構(gòu)建層級(jí)結(jié)構(gòu)模型。其中準(zhǔn)則層包括巖溶條件、覆蓋層條件、構(gòu)造條件、地下水條件，這是巖土體發(fā)生巖溶塌陷需要具備的基本因素；子準(zhǔn)則層則是地層巖性、巖體結(jié)構(gòu)特征、地下巖溶發(fā)育情況、覆蓋層結(jié)構(gòu)、覆蓋層厚度、覆蓋層巖性、斷層位置、斷層特征、斷層附近巖溶發(fā)育情況、地下水位變化幅度、地下水距基巖高度、地下水徑流強(qiáng)度、地下水抽取強(qiáng)度等，這是巖溶塌陷發(fā)生的誘發(fā)因素（圖2）。

圖2　研究區(qū)巖溶塌陷層次結(jié)構(gòu)模型Fig.2 Hierarchical structure of karst collapse in research area

3.2建立風(fēng)險(xiǎn)因子判斷矩陣及確定因子權(quán)重

3.2.1判斷矩陣結(jié)構(gòu) 根據(jù)兩兩因子間的重要程度判據(jù)（圖3），對(duì)風(fēng)險(xiǎn)因子進(jìn)行賦值，并結(jié)合圖2的巖溶塌陷層次結(jié)構(gòu)模型，建立不同層級(jí)間關(guān)聯(lián)判斷矩陣，如表1～表5所示。

表　1 目標(biāo)層　A-準(zhǔn)則層　Bi的關(guān)聯(lián)判斷矩陣Table 1 Related judgment matrix of target level A to criteria level Bi

表2　準(zhǔn)則層　B1-子準(zhǔn)則層　Ci的關(guān)聯(lián)判斷矩陣Table 2 Related judgment matrix of criteria level B1to secondary criteria level Ci

表3　準(zhǔn)則層　B2-子準(zhǔn)則層　Ci的關(guān)聯(lián)判斷矩陣Table 3 Related judgment matrix B2of criteria level B2to secondary criteria level Ci

圖3　風(fēng)險(xiǎn)因子重要性判據(jù)及賦值Fig.3 Importance judgment of risk factors and assessment score

表4　準(zhǔn)則層　B3-子準(zhǔn)則層　Ci的關(guān)聯(lián)判斷矩陣Table 4 Related judgment matrix of criteria level B3to secondary criteria level Ci

表5　準(zhǔn)則層　B4-子準(zhǔn)則層　Ci的關(guān)聯(lián)判斷矩陣Table 5 Related judgment matrix of criteria level B4to secondary criteria level Ci

3.2.2風(fēng)險(xiǎn)因子權(quán)重的確定及排序 ①以目標(biāo)層A-準(zhǔn)則層Bi的關(guān)聯(lián)判斷矩陣為例，求得n階矩陣行向量乘積的n次方根數(shù)fai為

同理可得

為

因此，判斷矩陣A滿足一致性檢驗(yàn)要求。

同理，求得其他判斷矩陣一致性檢驗(yàn)結(jié)果為

子準(zhǔn)則層權(quán)重向量為F=［Fb1，F(xiàn)b2，F(xiàn)b3，F(xiàn)b4］T，即F=［0.493 4，0.195 8，0.310 8，0.209 8，0.549 9，0.240 2，0.163 4，0.297 0，0.539 6，0.199 6，0.453 6，0.091 7，0.255 1］T，

歸一化后為

經(jīng)驗(yàn)算后符合一致性檢驗(yàn)要求，因此各子準(zhǔn)則層因子及對(duì)應(yīng)的風(fēng)險(xiǎn)因子權(quán)重如表6所示。對(duì)風(fēng)險(xiǎn)因子進(jìn)行權(quán)重大小排序?yàn)椋焊采w層厚度、地層巖性、斷層附近巖溶發(fā)育、地下水距基巖高度、地下巖溶發(fā)育、斷層特征、地下水抽取強(qiáng)度、覆蓋層巖性、覆蓋層結(jié)構(gòu)、地下水位變化幅度、巖體結(jié)構(gòu)特征、斷層位置、地下水徑流強(qiáng)度（圖4）。可見，覆蓋層厚度、地層巖性、斷層附近巖溶發(fā)育為塌陷內(nèi)在基礎(chǔ)條件，地下水距基巖高度、地下水抽取強(qiáng)度、地下水位變化幅度、地下水徑流強(qiáng)度為塌陷外在動(dòng)力條件，因而在覆蓋層較薄和巖溶發(fā)育程度比較高的區(qū)域進(jìn)行大型工程建設(shè)需要抽排地下水時(shí)應(yīng)采取有效的防護(hù)措施。

表6　子準(zhǔn)則層因子及對(duì)應(yīng)權(quán)重值Table 6 Related secondary criteria level factors to the weights of elements

圖4　巖溶塌陷風(fēng)險(xiǎn)因子權(quán)重排序Fig.4 Weight function ranking model of karst collapse risk factors

3.3巖溶塌陷風(fēng)險(xiǎn)因子等級(jí)劃分

由于各因素的地質(zhì)條件背景有差異，對(duì)塌陷的影響程度也有所不同，因此建立巖溶塌陷風(fēng)險(xiǎn)因子等級(jí)區(qū)劃，明確各部分評(píng)判標(biāo)準(zhǔn)和劃分范圍（表7）。

3.4巖溶塌陷危險(xiǎn)性區(qū)域評(píng)價(jià)

3.4.1構(gòu)建評(píng)價(jià)體系 根據(jù)當(dāng)?shù)丨h(huán)境地質(zhì)情況，結(jié)合表7對(duì)可能造成巖溶塌陷的危險(xiǎn)因子如地層巖性、巖體結(jié)構(gòu)特征、地下巖溶發(fā)育等均進(jìn)行綜合考慮，選取如下評(píng)價(jià)體系公式來求取巖溶塌陷危險(xiǎn)指數(shù)

式中：n=13；U為巖溶塌陷危險(xiǎn)指數(shù)，其值越大表示潛在塌陷危險(xiǎn)性越大，在評(píng)價(jià)圖上表現(xiàn)的顏色也愈深，相反則表示危險(xiǎn)性愈低；A表示巖溶塌陷危險(xiǎn)因子；Bi表示危險(xiǎn)因子權(quán)重值。

3.4.2評(píng)價(jià)實(shí)例 選取遵義市永樂鎮(zhèn)作為研究對(duì)象，對(duì)集鎮(zhèn)及周邊耕地均納入考慮范圍，采用1 m ×1 m的網(wǎng)格將集鎮(zhèn)劃分為若干區(qū)域，對(duì)該區(qū)域及周邊的踏勘、鉆探和地質(zhì)雷達(dá)等資料進(jìn)行分析后確定塌陷危險(xiǎn)因子，運(yùn)用層次分析法得出各因子權(quán)重；借助ArcGIS在圖形分析處理上區(qū)別于其他地理信息系統(tǒng)軟件的優(yōu)越性，對(duì)各危險(xiǎn)因子圖層屬性表進(jìn)行賦值等方式，結(jié)合式（1）得出巖溶塌陷危險(xiǎn)指數(shù)；最后，利用ArcGIS的空間分析功能，對(duì)數(shù)據(jù)參數(shù)進(jìn)行賦值、重分類以及多層面合并后得出巖溶塌陷危險(xiǎn)性評(píng)價(jià)圖（圖5），評(píng)價(jià)圖中圖形顏色的深淺表示發(fā)生巖溶塌陷可能性的大小，即穩(wěn)定區(qū)（Ⅰ）、次穩(wěn)定區(qū)（Ⅱ）、次危險(xiǎn)區(qū)（Ⅲ）、危險(xiǎn)區(qū)（Ⅳ）、極危險(xiǎn)區(qū)（Ⅴ）等5個(gè)區(qū)域。

3.4.3評(píng)價(jià)結(jié)果分析 結(jié)合當(dāng)?shù)氐匦蔚刭|(zhì)條件，對(duì)5個(gè)塌陷危險(xiǎn)性區(qū)域的分析如下。

（1）巖溶塌陷極危險(xiǎn)區(qū)。主要分布在北面法庭周邊、東面敬老院和南面鎮(zhèn)政府、郵電局及計(jì)生站一帶。該區(qū)域第四系覆蓋層厚度＜3 m，土體結(jié)構(gòu)較為混雜，為粘土、含碎石粉質(zhì)粘土、碎石土等。地層為中厚層含燧石灰?guī)r，下部有厚10 m左右的含碳泥質(zhì)灰?guī)r夾頁巖，頂部為薄層泥灰?guī)r及鈣質(zhì)粘土巖（厚約38 m）。區(qū)域內(nèi)洼地、漏斗和落水洞較多，鎮(zhèn)政府一帶為地下水富集區(qū)域，全年地下水徑流模數(shù)為7～13 L/（s·km2）。地下溶洞和裂隙強(qiáng)烈發(fā)育，溶洞最高1.8 m，充填物多為黏土或松散堆積物。綜合考慮該區(qū)域巖溶水文地質(zhì)條件充分滿足發(fā)生巖溶塌陷地質(zhì)災(zāi)害的最低限度，結(jié)合近年來發(fā)生巖溶塌陷的情況，將該區(qū)域劃入塌陷極危險(xiǎn)區(qū)是符合實(shí)際條件的。

表7　巖溶塌陷區(qū)域風(fēng)險(xiǎn)因素等級(jí)劃分Table 7 Gradation division of karst collapse area

圖5　研究區(qū)巖溶塌陷風(fēng)險(xiǎn)性評(píng)價(jià)圖Fig.5 Risk evaluation map of karst collapse in research area

（2）巖溶塌陷危險(xiǎn)區(qū)。分布于北北西和東南部，呈條帶狀覆蓋至國土所、鎮(zhèn)政府西部，該區(qū)第四系覆蓋層厚度為3～5 m，多為粘土和碎石粉質(zhì)粘土?？拷鼛r溶塌陷極危險(xiǎn)區(qū)，受到地下水波動(dòng)影響也較為強(qiáng)烈，下伏基巖小型溶洞及裂隙發(fā)育，綜合考慮將其納入危險(xiǎn)區(qū)范圍。

（3）巖溶塌陷次危險(xiǎn)區(qū)。分布于西北部，向西至派出所、居委會(huì)一帶，及東南邊雜木林一帶。覆蓋層厚度為5 m左右，為粘土或碎石粉質(zhì)粘土。偏離地下河補(bǔ)排管道及巖溶較發(fā)育地區(qū)，經(jīng)全面分析將其納入次危險(xiǎn)區(qū)。

（4）巖溶塌陷次穩(wěn)定區(qū)。主要分布于該區(qū)西部汽車站、南部磚廠附近一帶，覆蓋層為粘土，厚度為3～5 m?；鶐r溶蝕裂隙微發(fā)育，落水洞、漏斗發(fā)育較少，地下水波動(dòng)不大，因而確定為塌陷次穩(wěn)定區(qū)。

（5）巖溶塌陷穩(wěn)定區(qū)。主要分布于覆蓋型灰?guī)r不發(fā)育地區(qū)，位于南西面附近的黃泥坡、東部農(nóng)田一帶，地層巖性主要以黃綠、灰綠色頁巖及鈣質(zhì)頁巖為主，夾少量泥灰?guī)r，下伏基巖裂隙發(fā)育程度較低，發(fā)生巖溶塌陷的危險(xiǎn)很低。

4　結(jié)束語

（1）綜合覆蓋層厚度、地下水系富集程度和波動(dòng)頻率、下伏基巖溶蝕裂隙和孔洞發(fā)育程度，以及周邊斷層構(gòu)造影響等因素，永樂鎮(zhèn)發(fā)生巖溶塌陷可能性較高的地區(qū)主要位于北部和中部。評(píng)估結(jié)果與近年來集鎮(zhèn)發(fā)生巖溶塌陷的地點(diǎn)相吻合，表明運(yùn)用層次分析法和ArcGIS技術(shù)相結(jié)合的方法對(duì)巖溶塌陷危險(xiǎn)評(píng)估是可行的，其結(jié)果能為地質(zhì)災(zāi)害的防治、重點(diǎn)區(qū)域規(guī)劃以及土地開發(fā)利用等提供一定的科學(xué)依據(jù)和技術(shù)支持。

（2）研究區(qū)覆蓋層條件、可溶巖發(fā)育和地下水條件3個(gè)巖溶塌陷基本因素主要通過鉆孔資料獲得，在一定時(shí)段較好地體現(xiàn)了區(qū)內(nèi)巖溶發(fā)育程度。而地下巖溶發(fā)育有著潛伏性和隱蔽性，巖溶塌陷也有其突發(fā)性和不可預(yù)見性，在遭遇極端天氣情況或者樁基施工等外在誘發(fā)因素時(shí)，容易造成地下水位劇烈變化或土洞頂板突然垮塌。因此，建議在巖溶強(qiáng)烈發(fā)育區(qū)域內(nèi)建立巖溶塌陷監(jiān)測(cè)點(diǎn)，有針對(duì)性地對(duì)地下水氣壓力進(jìn)行監(jiān)測(cè)，在巖溶塌陷危險(xiǎn)區(qū)域?qū)?huì)起到很好的預(yù)警作用。

［1］雷明堂，蔣小珍，李瑜.巖溶塌陷模型試驗(yàn)——以武昌為例 ［J］.地質(zhì)災(zāi)害與環(huán)境保護(hù)，1993，4（2）：39-44.

［2］蔣小珍，雷明堂，管振德.單層土體結(jié)構(gòu)巖溶土洞的形成機(jī)理 ［J］.中國巖溶，2012，31（4）：426-432.

［3］朱壽增，周健紅，陳學(xué)軍.桂林市西城區(qū)巖溶塌陷形成條件及主要影響因素［J］.桂林工學(xué)院學(xué)報(bào)，2000，20（2）：100-105.

［4］王柳寧，高武振.桂林市西城區(qū)地下水活動(dòng)與巖溶塌陷的關(guān)系 ［J］.桂林工學(xué)院學(xué)報(bào)，2000，20（2）：106-111.

［5］劉寶臣，宋慶生，陳學(xué)軍，等.桂林市西城區(qū)巖溶塌陷特征分析 ［J］.地質(zhì)災(zāi)害與環(huán)境保護(hù)，2000，11（3）：200 -203.

［6］馮佐海，李曉峰，劉輝利，等.桂林市斷裂構(gòu)造對(duì)巖溶塌陷的控制作用［J］.桂林工學(xué)院學(xué)報(bào)，2000，20（4）：344-349.

［7］雷明堂，蔣小珍，李瑜，等.湘潭市巖溶塌陷的綜合預(yù)測(cè)與評(píng)價(jià) ［J］.中國地質(zhì)災(zāi)害與防治學(xué)報(bào)，1997，8（S）：69-77.

［8］陳學(xué)軍，陳植華，陳先華，等.桂林市西城區(qū)巖溶塌陷模糊層次綜合預(yù)測(cè)［J］.桂林工學(xué)院學(xué)報(bào)，2000，20（2）：112-116.

［9］包惠明，胡長順.巖溶地面塌陷神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)［J］.工程地質(zhì)學(xué)報(bào)，2002，10（3）：299-304.

［10］胡成，陳植華，丁國平，等.GIS技術(shù)在巖溶塌陷預(yù)測(cè)中的應(yīng)用［J］.桂林工學(xué)院學(xué)報(bào)，2000，20（2）：117-119.

［11］蒙彥，黃健民，雷明堂，等.基于灰色Verhulst模型的巖溶塌陷定量預(yù)報(bào)預(yù)測(cè)方法［J］.中國巖溶，2009，28（1）：17-22.

［12］劉秀敏，陳從新，沈強(qiáng)，等.覆蓋型巖溶塌陷的空間預(yù)測(cè)與評(píng)價(jià) ［J］.巖土力學(xué)，2011，32（9）：2785-2790.

［13］趙增玉，潘懋，梁河.二值證據(jù)權(quán)（Wofe）模型巖溶塌陷區(qū)劃研究［J］.北京大學(xué)學(xué)報(bào)：自然科學(xué)版，2010，46（4）：594-600.

［14］金艷珠，談樹成，虎雄崗，等.基于層次分析法與GIS相結(jié)合的巖溶塌陷地質(zhì)災(zāi)害易發(fā)性分區(qū)評(píng)估——以云南省師宗縣為例［J］.熱帶地理，2012，32（2）：173-178.

［15］潘健，周森，林培源，等.廣州市白云區(qū)巖溶塌陷風(fēng)險(xiǎn)初探［J］.巖土力學(xué)，2013，34（9）：2589-2600.

［16］蔣小珍，雷明堂，戴建玲，等.巖溶（土洞）塌陷TDR監(jiān)測(cè)試驗(yàn)研究［J］.水文地質(zhì)工程地質(zhì)，2011，38（1）：118-122.

［17］徐洪鐘，周元，張丹.基于GIS的巖溶塌陷分布式光纖監(jiān)測(cè)系統(tǒng)的研發(fā)［J］.水文地質(zhì)工程地質(zhì)，2011，38（3）：120-123.

［18］管振德，蔣小珍，高明.巖溶塌陷光纖傳感試驗(yàn)裝置的標(biāo)定試驗(yàn)［J］.中國巖溶，2012，31（2）：173-178.

［19］趙永輝，陳軍，王水強(qiáng)，等.高層建筑石灰?guī)r基底的雷達(dá)探測(cè)研究 ［J］.地下空間與工程學(xué)報(bào)，2006，2（4）：700-704.

［20］黎華清，盧呈杰，韋吉益，等.孔間電磁波CT探測(cè)揭示水庫壩基巖溶形態(tài)特征——以廣西靖西大龍?zhí)端畮灬∧还酀{為例［J］.巖土力學(xué)，2008，29（S）：607-610.

［21］李大虎.探地雷達(dá)在探測(cè)墩下隱伏巖溶中的應(yīng)用 ［J］.重慶工學(xué)院學(xué)報(bào)：自然科學(xué)版，2009，23（2）：50-54.

［22］張華，陳小宏，劉松.地震勘探技術(shù)在公路采空區(qū)調(diào)查中的應(yīng)用 ［J］.地下空間與工程學(xué)報(bào)，2011，7（2）：371-374.

［23］趙正君，蔣小珍，鄭世華.美國佛羅里達(dá)州巖溶塌陷保險(xiǎn)法案及啟示［J］.中國巖溶，2009，28（4）：432-438.

Risk evaluation of karst collapse based on technology of AHP and ArcGIS—A case of Yongle Town in Zunyi City

PAN Zong-yuan1，JIA Long1，LIU Bao-chen2
（1.Institute of Karst Geology and Key Laboratory of Karst Collapse Prevention，Chinese Academy of Geological Sciences，Guilin 541004，China；2.Guilin University of Technology，Guilin 541004，China）

The risk factors causing karst collapse of Yongle Town in Zunyi city are recognized and confirmed.In the hierarchical structure model there are 13 factors belonging to 4 conditions：karst development condition，overburden condition，geotectonic condition，groundwater condition.Combined with local hydrogeology background，weight of risk factors was calculated by Analytic Hierarchy Process（AHP）.The risk gradation division of karst collapse of each zone was confirmed，and risk evaluation of karst collapse in research area was defined by ArcGIS technology.The result shows that the most subsidence area is distributed on the central region of research area，which is presented as NNW irregular stripe.The more likely and likely subsidence areas are successively distributed at the two sides of the most likely subsidence area.The less stable and stable areas are respectively distributed at the outside or around the edge of south and the east.The results are adaptable to the geographical location and hazard degree of collapse in recent years.

karst collapse；risk factors；hierarchical structure model；ArcGIS；risk evaluation；Zunyi

P642；X43

A

1674-9057（2016）03-0464-07

10.3969/j.issn.1674-9057.2016.03.008

2015-10-20

國家自然科學(xué)基金項(xiàng)目（41402284）

潘宗源 （1987—），男，碩士，研究方向：巖溶地質(zhì)災(zāi)害防治及GIS應(yīng)用，zypan@karst.ac.cn。

引文格式：潘宗源，賈龍，劉寶臣.基于AHP和ArcGIS技術(shù)的巖溶塌陷風(fēng)險(xiǎn)評(píng)價(jià)——以遵義永樂鎮(zhèn)為例［J］.桂林理工大學(xué)學(xué)報(bào)，2016，36（3）：464-470.