劉宇恒，鄧 輝，熊倩瑩

(成都理工大學(xué) 地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，成都 610059)

基于層次分析法的茂縣斜坡地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)

劉宇恒，鄧 輝，熊倩瑩

(成都理工大學(xué) 地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，成都 610059)

茂縣是西南深切河谷地區(qū)地質(zhì)災(zāi)害的重災(zāi)縣之一，亦是對(duì)區(qū)域地質(zhì)災(zāi)害進(jìn)行特征分析及易發(fā)性評(píng)價(jià)的良好試驗(yàn)基地。在茂縣斜坡地質(zhì)災(zāi)害詳查及遙感數(shù)據(jù)解譯的基礎(chǔ)上，選取高程、地形起伏度、坡形、坡度、地層巖性、水系因子及地質(zhì)構(gòu)造7個(gè)因素作為評(píng)價(jià)指標(biāo)，采用層次分析法對(duì)研究區(qū)斜坡地質(zhì)災(zāi)害易發(fā)性進(jìn)行評(píng)價(jià)與區(qū)劃。結(jié)果表明：地質(zhì)災(zāi)害高易發(fā)區(qū)占地面積約46.1 km2，占10.02%，主要分布于研究區(qū)的岷江干流沿岸及其部分支流沿岸等；中易發(fā)區(qū)面積約149.59 km2，占32.52%；低易發(fā)區(qū)面積約205.25 km2，占44.62%；非易發(fā)區(qū)面積約59.06 km2，占12.84%。評(píng)價(jià)結(jié)果與野外調(diào)查基本吻合，可為茂縣斜坡地質(zhì)災(zāi)害防治及土體規(guī)劃利用提供依據(jù)。

斜坡地質(zhì)災(zāi)害；易發(fā)性評(píng)價(jià)；四川茂縣；地理信息系統(tǒng)；層次分析法

1 研究背景

四川省茂汶羌族自治縣地處西南深切河谷地區(qū)，境內(nèi)自然地質(zhì)背景環(huán)境特殊，地質(zhì)構(gòu)造復(fù)雜，新構(gòu)造運(yùn)動(dòng)活躍，河谷強(qiáng)烈侵蝕切割，各種滑坡、崩塌等斜坡地質(zhì)災(zāi)害普遍發(fā)育。尤其在2008年5月12日汶川地震之后，斜坡地質(zhì)災(zāi)害的發(fā)生越來(lái)越多，茂縣人民的生命財(cái)產(chǎn)安全受到很大威脅。由此可知，進(jìn)行地質(zhì)災(zāi)害調(diào)查和評(píng)價(jià)工作，科學(xué)分析地質(zhì)災(zāi)害的易發(fā)性，對(duì)茂縣斜坡地質(zhì)災(zāi)害的預(yù)防和減輕具有深遠(yuǎn)意義。

合理評(píng)價(jià)地質(zhì)災(zāi)害易發(fā)性十分重要，近年來(lái)許多專(zhuān)家學(xué)者對(duì)此進(jìn)行了大量研究[1-8]，但在對(duì)評(píng)價(jià)地質(zhì)災(zāi)害易發(fā)性的方法選擇上一直沒(méi)達(dá)成共識(shí)。由于層次分析法簡(jiǎn)易、靈活且邏輯性和實(shí)用性較強(qiáng)，已是一套較為成熟的方法，因此，本文基于層次分析法，對(duì)茂縣區(qū)內(nèi)斜坡地質(zhì)災(zāi)害易發(fā)性進(jìn)行評(píng)價(jià)。

2 研究區(qū)概況

茂汶羌族自治縣位于四川省西北部，是阿壩藏族羌族自治州下轄縣，地理坐標(biāo)在31°24′N(xiāo)─32°17′N(xiāo)， 102°56′E—104°10′E之間，南北寬約94.8 km，東西長(zhǎng)約116.5 km，占地面積約3 885.6 km2。研究區(qū)遙感圖如圖1所示。

圖1 研究區(qū)遙感圖Fig.1 Remotely-sensed image of the research area

茂縣屬于亞熱帶季風(fēng)氣候，多年年平均降水量為484.1 mm，容易出現(xiàn)冬干春旱的現(xiàn)象?？h境內(nèi)地質(zhì)結(jié)構(gòu)較為繁雜，大部分屬于馬爾康地質(zhì)分區(qū)，其東南部的狹窄部分屬于龍門(mén)山地質(zhì)分區(qū)。區(qū)內(nèi)巖土體類(lèi)型主要由夾沙灰、黑色頁(yè)巖、砂板巖及火成巖等組成。土壤主要為暗棕壤、褐土、棕壤。研究區(qū)地處龍門(mén)山地震地段，是中國(guó)地震頻繁的地區(qū)之一。岷江從北往南縱貫全縣，黑水河、赤不蘇河及松坪河分別在大、小兩河口與較場(chǎng)鄉(xiāng)匯進(jìn)岷江，土門(mén)河自西往東貫穿土門(mén)全區(qū)，匯進(jìn)涪江水系。這種奇特的區(qū)域地質(zhì)環(huán)境與頻繁的內(nèi)外動(dòng)力地質(zhì)條件，使得研究區(qū)內(nèi)深切河谷斜坡地質(zhì)災(zāi)害十分發(fā)育，且災(zāi)害類(lèi)型以滑坡、崩塌為主。經(jīng)現(xiàn)場(chǎng)調(diào)查，斜坡地質(zhì)災(zāi)害總共147處，其中滑坡災(zāi)害有98處，約占總數(shù)的67%，其他為崩塌災(zāi)害，共有49處，約占總數(shù)的33%，分布統(tǒng)計(jì)見(jiàn)表1。

表1 茂縣斜坡地質(zhì)災(zāi)害分布統(tǒng)計(jì)Table 1 Statistics of slope geo-hazard distribution in Maoxian county

3 茂縣斜坡地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)

3.1 層次分析法簡(jiǎn)介

層次分析法是一種把定性的問(wèn)題進(jìn)行定量化分析的多重準(zhǔn)則決議方法，其特點(diǎn)是將一個(gè)復(fù)雜的問(wèn)題分解成一系列有序?qū)哟蔚慕M成因素，再將專(zhuān)家的主觀經(jīng)驗(yàn)與研究區(qū)地質(zhì)環(huán)境特征基于一定客觀化的判別矩陣聯(lián)系起來(lái)，使用數(shù)學(xué)運(yùn)算方式定量地描述各層中影響因子的相對(duì)重要程度，并賦予其對(duì)應(yīng)的權(quán)重值，最后根據(jù)全部影響因子相對(duì)總目標(biāo)的相對(duì)權(quán)重值進(jìn)行排序。

本文嘗試采用層次分析法來(lái)建立斜坡地質(zhì)災(zāi)害評(píng)價(jià)因子的遞階層次結(jié)構(gòu)，構(gòu)造兩兩評(píng)價(jià)因子比較的判斷矩陣，給予其相對(duì)應(yīng)的權(quán)重值且進(jìn)行一致性檢驗(yàn)，再計(jì)算各層因素的組合權(quán)重[9]。

3.2 評(píng)價(jià)因子的選取

為突出研究區(qū)斜坡地質(zhì)災(zāi)害發(fā)生的主要影響因素，結(jié)合前人的研究成果，根據(jù)研究區(qū)內(nèi)地質(zhì)環(huán)境特征與斜坡地質(zhì)災(zāi)害的分布規(guī)律及形成條件，選取地形地貌(包括高程、地形起伏度、坡形及坡度)、地層巖性、水系因子及地質(zhì)構(gòu)造7個(gè)因素作為評(píng)價(jià)因子。

其中，在評(píng)價(jià)因子原始數(shù)據(jù)提取時(shí)，基于野外詳查，利用GIS軟件空間分析中的表面分析功能從DEM數(shù)據(jù)中提取出高程；利用ArcGIS軟件的鄰域分析功能，得到一定區(qū)域內(nèi)最低海拔與最高海拔的差值，作為中心柵格的地形起伏度；把剖面曲率值≤-0.5的歸為凹型坡，-0.5<剖面曲率值≤0.5的歸為直線型坡，剖面曲率值>0.5的歸為凸型坡；在GIS系統(tǒng)中利用緩沖分析工具，以河流為中心建立緩沖區(qū)，根據(jù)距離構(gòu)造斷裂帶的遠(yuǎn)近關(guān)系，以斷裂為中心創(chuàng)建緩沖區(qū)。

碳酸鹽礦物溶解過(guò)程中，生成游離態(tài)的Ca2+、Mg2+、CO2和H2O，消耗大部分礦物質(zhì)，且溶液中的Na+、H+隨溶液不斷滲入煤體中，置換出Ca2+、Mg2+，Na+、H+。由于煤樣中大部分碳酸鹽礦物發(fā)生溶解和吸附，溶液中離子不斷發(fā)生擴(kuò)散、遷移等離子交換反應(yīng)，煤樣孔裂隙周?chē)粩嘈纬煽锥矗姑簶涌紫督Y(jié)構(gòu)改變，孔隙率增大，孔隙體積等于反應(yīng)運(yùn)移溶質(zhì)體積。碳酸鹽礦物的溶解量與煤樣孔隙率之間的關(guān)系可以用以下公式表示：

3.3 層次分析法步驟

3.3.1 遞階層次結(jié)構(gòu)分析模型創(chuàng)建

創(chuàng)建遞階層次結(jié)構(gòu)是層次分析法的步驟中最重要的一步，必須要使同一層的各個(gè)單元附屬于其對(duì)應(yīng)的上一個(gè)層次的分類(lèi)，且必須確定與之對(duì)應(yīng)的下一層次中的每一個(gè)單元[10]。將前文中所選取的7個(gè)影響因子分類(lèi)并歸化，便可繪制斜坡地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)的3層遞階層次結(jié)構(gòu)分析模型(如圖2)。其中A層為目標(biāo)層，B層為準(zhǔn)則層，C層則為措施層。

圖2 層次結(jié)構(gòu)分析模型Fig.2 Analytical hierarchy model

3.3.2 判斷矩陣構(gòu)造

在準(zhǔn)則A的支配下按照其從屬元素兩兩之間的相對(duì)重要程度對(duì)B1，B2，B3，…，Bn賦予其相對(duì)應(yīng)的權(quán)重值。本文采用1—9標(biāo)度法[10](如表2)，針對(duì)各層中的要素對(duì)上一層次目標(biāo)的相對(duì)重要性進(jìn)行兩兩比較(表3至表5)，進(jìn)而構(gòu)造判斷矩陣。

表2 判斷矩陣標(biāo)度(重要性指標(biāo))及其含義Table 2 Scales of judgment matrixes and their connotations

表3 A-B層判斷矩陣Table 3 Judgment matrixes in A-B level

表4 B1-C層判斷矩陣Table 4 Judgment matrixes in B1-C level

表5 B2-C層判斷矩陣Table 5 Judgment matrixes in B2-C level

3.3.3 基于判斷矩陣的被比較因素相對(duì)權(quán)重的計(jì)算

運(yùn)用求和法對(duì)每個(gè)兩兩成對(duì)比所形成的判斷矩陣計(jì)算，從而可以得到與其相對(duì)應(yīng)的特征向量和最大特征根近似值λmax：

(1) 對(duì)兩兩相比較的矩陣A中的每一列向量進(jìn)行歸一化處理，得到

(1)

將矩陣進(jìn)行一致性檢驗(yàn)，計(jì)算公式為

(2)

其中，

(3)

式中：CR為一致性比例，只有當(dāng)CR<0.1時(shí)，矩陣的一致性檢驗(yàn)才算合格；n為兩兩比較的元素的個(gè)數(shù)；CI為一致性指標(biāo)；RI為隨機(jī)一致性指標(biāo)，其值由表6確定。

表6 隨機(jī)一致性指標(biāo)RI值Table 6 Numerical values of random consistent index RI

經(jīng)過(guò)計(jì)算，A-B層判斷矩陣的最大特征根近似值λmax=2，把最大特征向量歸一化得到B1，B2評(píng)價(jià)指標(biāo)的權(quán)重Wi，均為0.5。B1-C層判斷矩陣λmax=5.068 3，一致性指標(biāo)CI=0.017，隨機(jī)一致性指標(biāo)RI=1.12，則一致性比例CR=0.015 25<0.1，一致性檢驗(yàn)合格，該矩陣具有較好的判斷一致性。把最大特征向量歸一化得到C1，C2，C3，C4，C5這5個(gè)影響因素的權(quán)重值Wi，分別為0.062 4，0.098 6，0.261 8，0.416 2，0.161 0。而B(niǎo)2-C層判斷矩陣的最大特征根為λmax=2，將最大特征向量歸一化得到C6，C7評(píng)價(jià)指標(biāo)的權(quán)重Wi，分別為0.142 9，0.857 1。

3.3.4 組合權(quán)重計(jì)算

若要獲取每一個(gè)層次中全部因素相對(duì)總目標(biāo)的相對(duì)權(quán)重值，則需要計(jì)算最下層對(duì)目標(biāo)的組合權(quán)向量，并根據(jù)公式做組合一致性檢驗(yàn)，若檢驗(yàn)通過(guò)，則可按照組合權(quán)向量表示的結(jié)果進(jìn)行決策，否則需要重新考慮模型或重新構(gòu)造那些一致性比率較大的成對(duì)比較陣[10]。將B1，B2層歸一化處理后的權(quán)重值在A層次下再次歸一化處理即可得到7項(xiàng)評(píng)價(jià)要素的權(quán)重總排序，如表7所示。

表7 組合權(quán)重結(jié)果Table 7 Result of combinatorial weights

3.3.5 層次分析法定量評(píng)分

地質(zhì)災(zāi)害的易發(fā)性評(píng)價(jià)可以根據(jù)7個(gè)影響因子的權(quán)重值和打分完成，在分析評(píng)價(jià)的過(guò)程中，將每一個(gè)因子根據(jù)一定的準(zhǔn)則進(jìn)行打分，再將每一個(gè)指標(biāo)因子的賦值與其相對(duì)應(yīng)的權(quán)重相乘就可以求出地質(zhì)災(zāi)害的易發(fā)性值R，即

(4)

式中Bi為各評(píng)價(jià)因子的定量打分。

對(duì)指標(biāo)因素定量打分的準(zhǔn)則為將7個(gè)評(píng)價(jià)因素的權(quán)重值排序。將權(quán)重因子最小的因子開(kāi)始賦值為1分，權(quán)重倒數(shù)第二小的因素賦值為2分，以此類(lèi)推，依次從1到7分。按照每個(gè)因素的分級(jí)準(zhǔn)則，如果因子的測(cè)量值是最小的分級(jí)，就統(tǒng)一都賦值為1分，假如是屬于上一個(gè)分級(jí)，就按照它所在的7個(gè)因素中的排序位置加之其相應(yīng)的分值(比如，其權(quán)重值在7個(gè)因素中排在了第5位，它的評(píng)分值就賦予5分，如果測(cè)量值是屬于上一個(gè)分級(jí)，就在原來(lái)的基礎(chǔ)上加上5分)，之后向上按等差級(jí)數(shù)遞增[11]。但是，因?yàn)樗狄蜃訉?duì)地質(zhì)災(zāi)害發(fā)生的敏感性最高，且離水系越近，地質(zhì)災(zāi)害越易發(fā)生，因此，水系因子的分?jǐn)?shù)從35分倒次打到1分。層次分析法定量評(píng)分結(jié)果見(jiàn)表8。

表8 層次分析法定量評(píng)分結(jié)果Table 8 Result of quantitative rating by AHP

為了更為簡(jiǎn)單地進(jìn)行對(duì)比分析，將易發(fā)性指數(shù)歸一處理化，從而得到茂汶羌族自治縣幅區(qū)域地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)圖(圖3)。

圖3 層次分析法的易發(fā)性評(píng)價(jià)結(jié)果Fig.3 Result of susceptibility evaluation by AHP

3.4 斜坡地質(zhì)災(zāi)害易發(fā)性區(qū)劃

將層次分析法得出的區(qū)域地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)結(jié)果運(yùn)用等間距分級(jí)方法進(jìn)行分級(jí)，得到4個(gè)易發(fā)性評(píng)價(jià)等級(jí)：高易發(fā)區(qū)(易發(fā)性指數(shù)為(0.75,1.00])、中易發(fā)區(qū)(易發(fā)性指數(shù)為(0.50,0.75])、低易發(fā)區(qū)(易發(fā)性指數(shù)為(0.25,0.50])、非易發(fā)區(qū)(易發(fā)性指數(shù)為(0,0.25])。在此基礎(chǔ)上，結(jié)合茂縣的自然地質(zhì)環(huán)境如地形地貌、地質(zhì)構(gòu)造、巖石類(lèi)型、河流流域等，考慮環(huán)境保護(hù)及治理工作規(guī)劃，對(duì)上述4個(gè)易發(fā)性等級(jí)進(jìn)行輕微調(diào)整，去除異常點(diǎn)，得出茂縣區(qū)域斜坡地質(zhì)災(zāi)害易發(fā)性程度分區(qū)：高易發(fā)區(qū)、中易發(fā)區(qū)、低易發(fā)區(qū)、非易發(fā)區(qū)(見(jiàn)圖4)。

圖4 茂縣幅斜坡地質(zhì)災(zāi)害易發(fā)性區(qū)劃Fig.4 Zoning of slope geo-hazard susceptibility in Maoxian County

由圖4可知，高易發(fā)區(qū)占地面積約46.1 km2，占研究區(qū)總占地面積的10.02%，主要分布于研究區(qū)的岷江干流沿岸及其部分支流沿岸等；中易發(fā)區(qū)面積約149.59 km2。占研究區(qū)總占地面積的32.52%,主要緊鄰地質(zhì)災(zāi)害高易發(fā)區(qū)域分布；低易發(fā)區(qū)分布密度相對(duì)較低，面積約205.25 km2，占研究區(qū)總面積的44.62%；非易發(fā)區(qū)面積約59.06 km2，占研究區(qū)總面積的12.84%。

圖4斜坡地質(zhì)災(zāi)害主要集中分布在河流沿岸、人類(lèi)活動(dòng)強(qiáng)烈的地區(qū)，與現(xiàn)場(chǎng)調(diào)查成果擬合度較高，可見(jiàn)運(yùn)用層次分析法對(duì)茂縣斜坡地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)的結(jié)果較為理想。

4 結(jié) 論

(1) 本文基于層次分析法，建立以斜坡地質(zhì)災(zāi)害易發(fā)性綜合評(píng)價(jià)指數(shù)為目標(biāo)層，以高程、地形起伏度、坡形、坡度、地層巖性、水系因子及地質(zhì)構(gòu)造7個(gè)因素作為指標(biāo)層的茂縣斜坡地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)的層次分析模型，并確定了各基礎(chǔ)指標(biāo)層的因子對(duì)目標(biāo)層的權(quán)重。其中高程、地形起伏度、坡度、坡形、地層巖性、地質(zhì)構(gòu)造及水系因子所占權(quán)重分別為0.031 2，0.049 3，0.130 9，0.208 1，0.080 5，0.071 4，0.428 6。

(2) 對(duì)茂縣斜坡地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)及區(qū)劃表明，高易發(fā)性區(qū)占地面積約46.1 km2，占研究區(qū)總占地面積的10.02%，主要分布于研究區(qū)的岷江干流沿岸及其部分支流沿岸等；中易發(fā)區(qū)面積約149.59 km2，占研究區(qū)總占地面積的32.52%；低易發(fā)區(qū)面積約205.25 km2，占研究區(qū)總面積的44.62%；非易發(fā)區(qū)面積約59.06 km2，占研究區(qū)總面積的12.84%。

(3) 茂縣斜坡地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)及區(qū)劃成果與野外調(diào)查的地質(zhì)災(zāi)害隱患點(diǎn)分布基本切合，證明基于層次分析法的地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)結(jié)果是有效合理的，可為茂縣斜坡地質(zhì)災(zāi)害防治及土體規(guī)劃利用提供依據(jù)。

[1] 柴宗新. 山地災(zāi)害概念之我見(jiàn)[J]. 山地學(xué)報(bào), 1999, 17(1):91-94.

[2] 孟慶華, 孫煒?shù)h, 王 濤. 陜西鳳縣地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)研究[J]. 工程地質(zhì)學(xué)報(bào), 2011, 19(3):388-396.

[3] 王 哲, 易發(fā)成. 基于層次分析法的綿陽(yáng)市地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)[J]. 自然災(zāi)害學(xué)報(bào), 2009, 18(1):14-23.

[4] 唐 川, 許 強(qiáng). 強(qiáng)震區(qū)城市地質(zhì)災(zāi)害風(fēng)險(xiǎn)管理的研究?jī)?nèi)容與方法探討[J]. 工程地質(zhì)學(xué)報(bào), 2009, 17(1):56-61.

[5] 趙 成, 張永軍, 趙玉紅. 層次分析法在甘肅省地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)中的應(yīng)用[J]. 冰川凍土, 2009, 31(1):186-192.

[6]CASCINI L, BONNARD C. Landslide Hazard and Risk Zoning for Urban Planning and Development[J]. Landslide Risk Management, 2005,(7)：199-235．

[7] CARRARA A, GUZZETTI F. Use of GIS Technology in the Prediction and Monitoring of Landslide Hazard[J] . Natural Hazards, 1999,20(2)：117-135．

[8] 唐 川, 朱 靜. 瀾滄江中下游滑坡泥石流分布規(guī)律與危險(xiǎn)區(qū)劃[J]. 地理學(xué)報(bào), 1999，54(6):84-92.

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

[10]谷天峰, 王家鼎. 地理信息系統(tǒng)在黃土工程地質(zhì)中的應(yīng)用[M].北京：地震出版社, 2013.

[11]鐵永波, 唐 川. 層次分析法在單溝泥石流危險(xiǎn)度評(píng)價(jià)中的應(yīng)用[J]. 中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào), 2006, 17(4):79-84.

(編輯：劉運(yùn)飛)

AHP-based Evaluation of Slope Geo-hazard Susceptibility ofMaoxian County, Sichuan, China

LIU Yu-heng, DENG Hui, XIONG Qian-ying

(State Key Laboratory of Geohazard Prevention and Geoenvironment Protection，Chengdu University of Technology，Chengdu 610059，China)

Maoxian county is located in the deep valley region of southwest China with severe slope geological hazards, hence is a good experimental base for characteristic research and susceptibility evaluation of geological disasters. In this research, Analytical Hierarchy Process (AHP) was adopted for the evaluation and zoning of geological disaster susceptibility in the study area based on geo-hazard investigation and remote sensing data interpretation. Seven indices inclusive of height, land relieve, slope shape, gradient, lithology, water, and geological structure were selected as evaluation indicators. Result show that the area of high susceptibility covers 46.1 km2, accounting for 10.02% and distributing mainly along mainstream Minjiang River and some tributaries; the area of medium susceptibility covers 149.59 km2, accounting for 32.52%; and area of low susceptibility is 205.25 km2, taking up 44.62%; and the stable area covers 59.06 km2, comprising 12. 84%．The results are consistent with field investigation results in general, and thus could be taken as basis for land use planning and slope geo-hazard prevention and reduction in Maoxian county.

slope geo-hazard; susceptibility evaluation; Maoxian county; GIS; AHP

2016-02-26；

2016-04-05

中國(guó)地質(zhì)調(diào)查局項(xiàng)目(12120113010500)

劉宇恒(1991-)，男，四川眉山人，碩士研究生，從事巖土工程及地質(zhì)災(zāi)害方面的研究，(電話)18121960708(電子信箱)348196128@qq.com。

10.11988/ckyyb.20160155

2017,34(5):31-35

P694

A

1001-5485(2017)05-0031-05