申懷飛,董 雨,楊 梅,施采琦,楊欣柯

(許昌學(xué)院 城市與環(huán)境學(xué)院,河南 許昌 461000)

滑坡作為六大地質(zhì)災(zāi)害之一，其年發(fā)生量在地質(zhì)災(zāi)害中占比高達(dá)50%～86%，嚴(yán)重威脅人類(lèi)生命財(cái)產(chǎn)安全?；乱装l(fā)性評(píng)價(jià)(又稱(chēng)滑坡敏感性評(píng)價(jià))，是通過(guò)研究某一地區(qū)環(huán)境中的各項(xiàng)影響因子對(duì)滑坡發(fā)生的影響程度，綜合考慮各因子的權(quán)重，對(duì)潛在滑坡的類(lèi)型、空間分布和滑坡概率進(jìn)行定性或定量的評(píng)價(jià)，即解決“在何地容易發(fā)生滑坡災(zāi)害”的問(wèn)題。而滑坡風(fēng)險(xiǎn)性評(píng)估不僅回答了什么地方容易滑坡的問(wèn)題，還進(jìn)一步估算了因滑坡對(duì)個(gè)人、群體、財(cái)產(chǎn)或環(huán)境造成的風(fēng)險(xiǎn)[1-2]?；乱装l(fā)性和風(fēng)險(xiǎn)評(píng)估為當(dāng)前地質(zhì)災(zāi)害研究領(lǐng)域的熱點(diǎn)問(wèn)題。

國(guó)外學(xué)者在滑坡災(zāi)害易發(fā)性和風(fēng)險(xiǎn)評(píng)估方面開(kāi)展了一系列研究，Wicaksono等選取包括坡度曲率、坡度梯度和建筑密度等17個(gè)因子，采用層次分析法(Analytic Hierarchy Process，AHP)創(chuàng)建滑坡易發(fā)性圖[3]。Ellen Felizardo Batista等提出將衡量滑坡成本和評(píng)估脆弱性的方法應(yīng)用于公路周邊區(qū)域滑坡風(fēng)險(xiǎn)評(píng)估[4]。Pratap Ram等采用二元統(tǒng)計(jì)法繪制喜馬拉雅山脈附近小鎮(zhèn)及其周邊地區(qū)的滑坡易發(fā)圖并進(jìn)行了風(fēng)險(xiǎn)評(píng)估[5]。Abhirup Dikshit等綜合研究區(qū)實(shí)際情況和專(zhuān)家判斷確定脆弱性因子和危險(xiǎn)性因子的權(quán)重，繪制滑坡危險(xiǎn)圖和脆弱性圖，并對(duì)由降雨引起的山體滑坡風(fēng)險(xiǎn)進(jìn)行評(píng)估[6]。M.Fressard等以法國(guó)諾曼底的奧吉高原為研究區(qū)，在詳細(xì)尺度上使用多元統(tǒng)計(jì)建模方法對(duì)滑坡敏感性進(jìn)行評(píng)估[7]。Mostafa Mousavi等采用定量風(fēng)險(xiǎn)分析方法對(duì)伊朗北部Emamzadeh Ali地區(qū)由地震誘發(fā)的滑坡進(jìn)行脆弱性評(píng)估[8]。Biplab Mandal和Sujit Mandal應(yīng)用層次分析法，以地理空間工具為基礎(chǔ)，編制了大吉嶺喜馬拉雅東部利什河流域滑坡易感性圖[9]。Kayastha等采用層次分析法對(duì)尼泊爾蒂瑙流域進(jìn)行分析，編制了滑坡易感圖，并利用物理和統(tǒng)計(jì)方法對(duì)滑坡易感圖進(jìn)行了驗(yàn)證[10]。國(guó)內(nèi)研究方面，吳樹(shù)仁等確定了近年來(lái)對(duì)于滑坡災(zāi)害的研究逐步向以定性—定量描述為基礎(chǔ)預(yù)測(cè)災(zāi)害發(fā)生概率的方向發(fā)展[11]。向喜瓊建立了較為完善的區(qū)域滑坡地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)方法體系，完善和修正了區(qū)域滑坡地質(zhì)災(zāi)害風(fēng)險(xiǎn)評(píng)價(jià)和風(fēng)險(xiǎn)管理的基本框架[12]。常建娥等通過(guò)應(yīng)用實(shí)例證明應(yīng)用層次分析法可以確定出權(quán)重，從而為決策提供依據(jù)[13]。朱良峰等采用信息量模型對(duì)中國(guó)范圍內(nèi)的滑坡災(zāi)害進(jìn)行危險(xiǎn)性分析，進(jìn)而進(jìn)行區(qū)域社會(huì)經(jīng)濟(jì)易損性分析，并在此基礎(chǔ)上進(jìn)行最終的滑坡災(zāi)害風(fēng)險(xiǎn)評(píng)估[14]。孫強(qiáng)等基于GIS與層次分析法通過(guò)分析滑坡危險(xiǎn)性評(píng)價(jià)結(jié)果、人口及財(cái)產(chǎn)易損性等，對(duì)龍溪流域進(jìn)行滑坡風(fēng)險(xiǎn)性評(píng)價(jià)[15]。廖祥東等采用定量評(píng)價(jià)方法，開(kāi)展十堰市某單個(gè)滑坡的穩(wěn)定性、危險(xiǎn)性、承災(zāi)體易損性分析評(píng)價(jià)，在此基礎(chǔ)上進(jìn)行滑坡風(fēng)險(xiǎn)分析評(píng)價(jià)[16]。上述研究中的滑坡風(fēng)險(xiǎn)評(píng)估多在中小尺度范圍進(jìn)行，以選取某一種方法或模型對(duì)滑坡風(fēng)險(xiǎn)進(jìn)行評(píng)價(jià)為主流，且評(píng)價(jià)模型理論體系尚需進(jìn)一步完善。

甘肅地處黃土高原、青藏高原和內(nèi)蒙古高原三大高原交匯地帶，地形復(fù)雜，山脈縱橫交錯(cuò)，海拔相差懸殊，有滑坡致災(zāi)的充分條件。本文在綜合考慮各方面因素基礎(chǔ)上，選取地形起伏度、土壤、植被、降水、地震等指標(biāo)[17-25]，結(jié)合層次分析法與信息量模型對(duì)甘肅省境內(nèi)進(jìn)行滑坡易發(fā)性評(píng)估，為研究區(qū)滑坡災(zāi)害防治提供科學(xué)依據(jù)，提高人們對(duì)滑坡災(zāi)害的防范意識(shí)，也豐富了地質(zhì)災(zāi)害特別是滑坡災(zāi)害易發(fā)性評(píng)估的理論和方法[26-28]。

1 研究區(qū)概況

甘肅省位于中國(guó)西部，地處黃河上游，位于東經(jīng)92°20′—108°43′，北緯32°36′—42°48′，東西蜿蜒1 655 km，面積42.59萬(wàn)km2。甘肅省大部分地區(qū)位于中國(guó)地勢(shì)第二階梯上，海拔高程556～5 856 m；土壤類(lèi)型主要為漠土、高山土和初育土，占甘肅省總面積的60.28%；地層巖性主要為新生代的第四系和中古生代的泥盆系。甘肅省東西氣候差異大，西部地區(qū)屬溫帶大陸性氣候，東部屬于溫帶季風(fēng)氣候，降雨在時(shí)空分布上不均衡，西北少東南多，冬春少夏秋多，年降水量36.6～734.9 mm。2017年甘肅省滑坡多集中在東南部地區(qū)，共123處，其中小型10處，中型72處，大型38處，特大型3處。研究區(qū)位置見(jiàn)圖1。

圖1 研究區(qū)位置及滑坡點(diǎn)分布情況

2 研究數(shù)據(jù)和方法

2.1 數(shù)據(jù)來(lái)源及預(yù)處理

進(jìn)行滑坡易發(fā)性評(píng)估，基礎(chǔ)數(shù)據(jù)的來(lái)源及其精確程度將會(huì)直接對(duì)評(píng)價(jià)結(jié)果產(chǎn)生影響，本論文所采用數(shù)據(jù)來(lái)源見(jiàn)表1。本研究所用數(shù)據(jù)(表1)均選用WGS_1984_UTM_Zone_47N投影坐標(biāo)系，避免在不同坐標(biāo)系下各指標(biāo)數(shù)據(jù)的偏移和變形。將統(tǒng)一坐標(biāo)系的甘肅省DEM(Digital Elevation Model)數(shù)據(jù)、年平均降水、地震核密度、植被覆蓋度、土壤類(lèi)型、地質(zhì)類(lèi)型、土地利用類(lèi)型、道路密度及人口密度數(shù)據(jù)，基于ArcGIS軟件分別進(jìn)行剪裁、矢量化、拓?fù)錂z查、插值分析、柵格計(jì)算等處理，以便于對(duì)研究區(qū)滑坡易發(fā)性進(jìn)行評(píng)估。本研究將2017年滑坡點(diǎn)數(shù)據(jù)隨機(jī)分為兩部分，一部分用于計(jì)算不同因子對(duì)滑坡提供的信息量值，另一部分用于滑坡易發(fā)性評(píng)估的檢驗(yàn)。

表1 研究區(qū)基礎(chǔ)數(shù)據(jù)來(lái)源

2.2 評(píng)價(jià)指標(biāo)

對(duì)滑坡災(zāi)害易發(fā)性進(jìn)行評(píng)估，首先確定影響滑坡災(zāi)害的重要因子。本文選取滑坡影響因素中的地形起伏度、降水、坡度、植被覆蓋度、巖性、人口密度等因子(圖2)用于甘肅省滑坡易發(fā)性評(píng)估。

圖2 滑坡易發(fā)性評(píng)估指標(biāo)

2.3 研究方法

國(guó)內(nèi)外專(zhuān)家學(xué)者針對(duì)滑坡易發(fā)性評(píng)估提出多種研究方法。目前常用于滑坡易發(fā)性分析的方法主要為定性的層次分析法和定量的信息量模型等。定性分析法人為主觀因素較強(qiáng)，隨意性大；而定量分析法則忽略了不同因素之間的聯(lián)系。本研究采用定性與定量相結(jié)合的方法對(duì)甘肅省滑坡災(zāi)害進(jìn)行易發(fā)性評(píng)估。

2.3.1 層次分析法

(1)建立層次結(jié)構(gòu)模型。在此步驟中將問(wèn)題所含要素進(jìn)行分組，每組作為一個(gè)層次，并將其按目標(biāo)層—準(zhǔn)則層—措施層的次序排列起來(lái)，見(jiàn)圖3。

圖3 層次結(jié)構(gòu)模型

(2)構(gòu)造判斷矩陣。構(gòu)造判斷矩陣是AHP決策中一個(gè)關(guān)鍵的步驟。從目標(biāo)層開(kāi)始，對(duì)同一層的各個(gè)因素重要程度進(jìn)行兩兩比較，評(píng)定該層次中各有關(guān)元素的相對(duì)重要性(表2)。

表2 判斷矩陣標(biāo)度說(shuō)明

(3)一致性檢驗(yàn)。為了確定矩陣的可靠性需要檢驗(yàn)判斷矩陣的一致性，計(jì)算它的一致性指標(biāo)：

CI=(λmax-n)/(n-1)

(1)

式中：λmax為判斷矩陣的最大特征根；n為判斷矩陣階數(shù)；CI為判斷矩陣的一致性指標(biāo)。當(dāng)CI=0時(shí)，判斷矩陣具有完全一致性；反之，CI越大，就表示判斷矩陣的一致性越差。RI為判斷矩陣的平均隨機(jī)一致性指標(biāo)(表3)。

表3 平均隨機(jī)一致性指標(biāo)

一致性指標(biāo)CI與同階的平均隨機(jī)一致性指標(biāo)RI之比記為CR。

CR=CI/RI

當(dāng)CR<0.1時(shí)，認(rèn)為判斷矩陣有令人滿(mǎn)意的一致性。

計(jì)算得各指標(biāo)的權(quán)重向量ω=(0.0618，0.0848，0.0337，0.2197，0.0783，0.1974，0.1243，0.1187，0.0499，0.0314)

2.3.2 信息量法 本研究選定的10個(gè)影響因子對(duì)滑坡災(zāi)害發(fā)生所提供的信息量各不相同。根據(jù)公式可計(jì)算出各個(gè)因子所提供的信息量值：

(2)

式中：Si0為每個(gè)影響因子X(jué)中各個(gè)類(lèi)別中發(fā)生滑坡災(zāi)害的面積；Si為影響因子X(jué)中每個(gè)類(lèi)別的面積；A0為研究區(qū)內(nèi)發(fā)生滑坡災(zāi)害的總面積；A為研究區(qū)總面積。

計(jì)算得各因子對(duì)滑坡災(zāi)害發(fā)生所提供信息量值，計(jì)算結(jié)果見(jiàn)表4。

表4 各因子信息量值

將信息量值與層次分析法所得權(quán)重值進(jìn)行加權(quán)求和，公式如下：

(3)

式中：ICI為評(píng)價(jià)因子的綜合信息量;wi為層次分析法所得各評(píng)價(jià)因子權(quán)重值;Iij為第i個(gè)評(píng)價(jià)因子中第j個(gè)類(lèi)別的信息量。

3 結(jié)果與分析

3.1 評(píng)價(jià)結(jié)果及精度驗(yàn)證

利用公式(3)計(jì)算得到整個(gè)研究區(qū)綜合信息量值，進(jìn)而得到甘肅省滑坡易發(fā)性區(qū)劃圖(圖4)。從圖中可看出，在高易發(fā)區(qū)和極高易發(fā)區(qū)滑坡點(diǎn)分布集中，隨著易發(fā)性等級(jí)降低，滑坡的數(shù)量和密度分布都隨之降低，說(shuō)明甘肅省實(shí)際發(fā)生滑坡情況與評(píng)估模型基本吻合。

受試者工作特征曲線(receiver operating characteristic curve，ROC)又稱(chēng)為感受性曲線。ROC曲線分析的是二元分類(lèi)模型。在對(duì)滑坡進(jìn)行易發(fā)性評(píng)估時(shí)通常使用ROC曲線分析方法對(duì)評(píng)價(jià)結(jié)果進(jìn)行精度驗(yàn)證。ROC曲線是以真陽(yáng)性率(靈敏度)為縱坐標(biāo)，假陽(yáng)性率(1-特異度)為橫坐標(biāo)繪制的曲線[29]。試驗(yàn)精度與曲線下面積(Area Under Curve，AUC)取值成正比，AUC值越大，試驗(yàn)精度越高，即在ROC曲線圖中，ROC曲線越靠近圖表區(qū)域左上角，模型的預(yù)測(cè)精度越高。AUC處于[0.5,0.7]時(shí)預(yù)測(cè)可信但精度較低，AUC處于[0.7，0.9]時(shí)試驗(yàn)精度為中等水平，0.9

為了對(duì)滑坡預(yù)測(cè)結(jié)果進(jìn)行驗(yàn)證，選取2017年的123個(gè)滑坡點(diǎn)數(shù)據(jù)作為滑坡樣本，賦值為1，創(chuàng)建123個(gè)隨機(jī)樣本點(diǎn)作為非滑坡樣本，賦值為0?；谘芯繀^(qū)滑坡易發(fā)性評(píng)價(jià)結(jié)果利用SPSS繪制ROC曲線，見(jiàn)圖5，AUC值為0.884(即ROC曲線下面積)，面積的標(biāo)準(zhǔn)誤差為0.022，在顯著性水平在0.05的情況下置信區(qū)間為0.841～0.928，由此可知評(píng)估具有較高可信度，可作為甘肅省滑坡易發(fā)性評(píng)價(jià)最終結(jié)果。

圖5 甘肅省易發(fā)性區(qū)劃驗(yàn)證ROC曲線

3.2 滑坡易發(fā)性空間分布特征

由甘肅省滑坡易發(fā)性區(qū)劃圖(圖4)可知，甘肅省滑坡災(zāi)害易發(fā)等級(jí)整體空間分布特征呈現(xiàn)自東南向西北逐漸降低的變化規(guī)律，甘肅省滑坡災(zāi)害多以中低等級(jí)為主，多分布在中部和西部地區(qū)；高易發(fā)區(qū)則主要分布在甘肅省東南部的隴南、天水、甘南、定西、平?jīng)龊蛻c陽(yáng)等地，酒泉、張掖和武威等地也有零星分布；極高易發(fā)區(qū)集中分布在隴南地區(qū)。

圖4 甘肅省滑坡易發(fā)性區(qū)劃

3.3 不同等級(jí)易發(fā)性滑坡分布特征

對(duì)各滑坡易發(fā)區(qū)面積及其間滑坡點(diǎn)分布情況進(jìn)行統(tǒng)計(jì)得到表5。由表5和圖4可知，甘肅省滑坡極低易發(fā)區(qū)的面積為10 km2，約占全區(qū)總面積的0.003%，無(wú)滑坡災(zāi)害點(diǎn)分布。該易發(fā)區(qū)主要分布于甘肅省酒泉市西部及張掖市的肅南裕固族自治縣、高臺(tái)縣、臨澤縣地區(qū)。區(qū)域內(nèi)地形起伏度較小，地質(zhì)環(huán)境好，地層巖性主要為三疊系、泥盆系；土壤主要為鈣層土、高山土；植被覆蓋度高，降雨量少；受人類(lèi)活動(dòng)影響小，交通不發(fā)達(dá)，不易發(fā)生滑坡。低易發(fā)區(qū)的面積為1.35×105km2，約占整個(gè)研究區(qū)總面積的32.875%，滑坡點(diǎn)在此區(qū)間共分布1個(gè)，占總滑坡數(shù)的0.813%，滑坡災(zāi)害密度為0.074個(gè)/萬(wàn)km2。該易發(fā)區(qū)主要分布于酒泉市西部、東南部，武威市東北部及張掖市西北部和中部。區(qū)域內(nèi)地形起伏度較小，地層巖性復(fù)雜，土壤類(lèi)型為漠土、初育土和高山土，雖該部分地區(qū)處于山地，植被覆蓋度低，但人類(lèi)活動(dòng)少，交通不發(fā)達(dá)，因而滑坡易發(fā)性低。中易發(fā)區(qū)的面積為1.96×105km2，約占整個(gè)研究區(qū)面積的47.544%，滑坡點(diǎn)在此區(qū)間共分布40個(gè)，占總滑坡數(shù)的32.520%，滑坡災(zāi)害密度為2.041個(gè)/萬(wàn)km2。該易發(fā)區(qū)分布范圍廣，甘肅省各市均有涉及。區(qū)內(nèi)部分地區(qū)地形起伏度相對(duì)較大，土壤類(lèi)型多樣，包含漠土、初育土、鈣層土、高山土、干旱土、半水成土、鹽堿土，地層巖性復(fù)雜，人類(lèi)活動(dòng)相對(duì)較頻繁，植被覆蓋度相對(duì)較低，白銀市、蘭州市、定西市、慶陽(yáng)市及臨夏回族自治州滑坡災(zāi)害的發(fā)生受降雨量和交通影響較大。高易發(fā)區(qū)的面積為6.28×104km2，約占整個(gè)研究區(qū)總面積的15.243%，滑坡點(diǎn)在此區(qū)間共分布59個(gè)，占總滑坡數(shù)的47.967%，滑坡災(zāi)害密度為9.392個(gè)/萬(wàn)km2。該易發(fā)區(qū)主要分布在甘肅省的東南部和南部的慶陽(yáng)市、平?jīng)鍪小⑻焖?、定西市及甘南藏族自治州等地區(qū)。區(qū)內(nèi)地形起伏度較大，土壤類(lèi)型以半淋溶土、人為土、鈣層土、初育土為主，地層巖性主要為二疊系、三疊系、上三系和白堊系，植被覆蓋度較高,但相對(duì)頻繁的人類(lèi)活動(dòng)、過(guò)高的降雨量以及發(fā)達(dá)的交通對(duì)該區(qū)域較高的滑坡易發(fā)性具有一定的影響。極高易發(fā)區(qū)的面積為1.79×104km2，僅占研究區(qū)總面積的4.336%，滑坡點(diǎn)在此區(qū)間共分布23個(gè)，占總滑坡數(shù)的18.699%，其滑坡災(zāi)害密度達(dá)12.870個(gè)/萬(wàn)km2。該易發(fā)區(qū)主要分布在甘肅省的隴南市及甘南藏族自治州東南部地區(qū)。區(qū)內(nèi)土壤類(lèi)型主要為淋溶土，地層巖性以泥盆系、白堊系和第四系為主，植被覆蓋度高，人類(lèi)活動(dòng)相對(duì)較少，地形起伏度大，交通相對(duì)發(fā)達(dá)，年均降雨量達(dá)518.6～589.8 mm，危險(xiǎn)性高，易發(fā)生滑坡。

表5 甘肅省滑坡易發(fā)性評(píng)估檢驗(yàn)表

甘肅省中北部地區(qū)屬溫帶大陸性氣候，全年降水較少，雨水對(duì)土壤的沖蝕能力較弱，且該地區(qū)地形起伏度及坡度小，人類(lèi)活動(dòng)、交通等對(duì)地表影響較小，滑坡易發(fā)性低；甘肅省南部地區(qū)屬溫帶季風(fēng)氣候，全年降水多集中夏季，短時(shí)降水多，且降水強(qiáng)度大，對(duì)土層沖刷力強(qiáng)，該地區(qū)坡度大，土層不穩(wěn)定，再加之人類(lèi)活動(dòng)密集、交通發(fā)達(dá)，滑坡易發(fā)性較高。據(jù)此可推斷出甘肅省滑坡高發(fā)地區(qū)具有年均降雨量高，地形起伏度大、人類(lèi)活動(dòng)頻繁，交通發(fā)達(dá)等特征。

3.4 各地市滑坡易發(fā)性情況

對(duì)各地市各滑坡災(zāi)害易發(fā)區(qū)面積進(jìn)行統(tǒng)計(jì)，由表6可知，天水市、隴南市和甘南藏族自治區(qū)總體上滑坡易發(fā)程度較高，高易發(fā)和極高易發(fā)區(qū)面積分別占該市總面積的67.474%，81.937%和54.856%；平?jīng)鍪泻投ㄎ魇兄懈咭装l(fā)區(qū)面積占比較高，分別為該市面積的98.376%和93.848%；蘭州市、金昌市、白銀市、張掖市、慶陽(yáng)市、臨夏回族自治區(qū)中易發(fā)區(qū)分布最廣，分別為76.375%，78.076%，78.941%，55.430%，71.887%，65.469%；嘉峪關(guān)市中易發(fā)區(qū)和低易發(fā)區(qū)面積相當(dāng)，兩種易發(fā)區(qū)面積分別占該市面積的50.028%和49.959%；武威市和酒泉市大部分的地區(qū)為中低易發(fā)區(qū)，分別占該市面積的92.617%和98.839%。

表6 各地市滑坡易發(fā)情況

綜上，甘肅省滑坡易發(fā)性總體上多為中度易發(fā)，有少數(shù)地區(qū)高易發(fā)區(qū)或低易發(fā)區(qū)面積占比較大。高易發(fā)地區(qū)主要位于天水市、平?jīng)鍪?、隴南市以及甘南藏族自治區(qū)。這些地區(qū)人口密度相對(duì)較高，地形起伏度及坡度較大，受降水影響較大。

4 結(jié) 論

基于層次分析法和信息量法對(duì)甘肅省滑坡災(zāi)害易發(fā)性進(jìn)行評(píng)估，經(jīng)ROC曲線檢驗(yàn)，本研究的AUC值為0.884，具有較高的可信水平，預(yù)測(cè)結(jié)果具有準(zhǔn)確性。研究結(jié)論如下：

(1)甘肅省滑坡災(zāi)害易發(fā)性呈現(xiàn)集中分布特征，主要集中在東南部地區(qū)，且多位于人口密度較高、交通線較密集地區(qū)。

(2)研究區(qū)極高、高、中、低、極低易發(fā)區(qū)所占面積分別為1.787，6.282，19.595，13.549，0.001萬(wàn)km2，其中極低、低、中易發(fā)區(qū)占研究區(qū)總面積的80.421%。由此可知，甘肅省大部分區(qū)域滑坡易發(fā)性較小。

(3)甘肅省極高、高易發(fā)性滑坡的區(qū)域總面積為8.069萬(wàn)km2，占全省面積的19.579%，且大部分預(yù)測(cè)滑坡點(diǎn)多發(fā)生在甘南地區(qū)，與實(shí)際情況相符。

研究中所用數(shù)據(jù)存在數(shù)據(jù)不完備、匹配度不高、分辨率不同等問(wèn)題,對(duì)研究結(jié)果精度也會(huì)產(chǎn)一定影響;滑坡發(fā)生的動(dòng)力因素及物質(zhì)因素較為匱乏，難以全面地、客觀地進(jìn)行滑坡評(píng)估；影響因子的選取方面可能不夠全面，忽視了其他一些因素的影響作用等，以上問(wèn)題均留待進(jìn)一步研究。