李衛(wèi)波，侯可斌

（1.中國運載火箭技術(shù)研究院環(huán)境影響評價與監(jiān)測中心北京航天計量測試技術(shù)研究所，北京100076；2.北京市環(huán)境影響評價評估中心，北京100161）

陜西省重力侵蝕災(zāi)害分布特征規(guī)律研究

李衛(wèi)波1，侯可斌2

（1.中國運載火箭技術(shù)研究院環(huán)境影響評價與監(jiān)測中心北京航天計量測試技術(shù)研究所，北京100076；2.北京市環(huán)境影響評價評估中心，北京100161）

以陜西省重力侵蝕災(zāi)害數(shù)據(jù)為基礎(chǔ)，選取地形起伏度與年均大雨天數(shù)作為地形與降雨因子的參數(shù)，研究了陜西省重力侵蝕分布特征。結(jié)果表明，陜西省重力侵蝕依據(jù)其主導(dǎo)因子的不同可以劃分為三個區(qū)，分別為：重力侵蝕低發(fā)區(qū)（Ⅰ區(qū)）、地形因子主控區(qū)（Ⅱ區(qū)）和重力侵蝕高發(fā)區(qū)（III區(qū)）。將重力侵蝕分區(qū)結(jié)果同地貌分區(qū)對比，發(fā)現(xiàn)兩者存在較好的吻合關(guān)系，其中重力侵蝕低發(fā)區(qū)包括關(guān)中平原、風(fēng)沙區(qū)及黃土高原區(qū)，地形因子主控區(qū)包括秦嶺與安康盆地，重力侵蝕高發(fā)區(qū)則主要在大巴山區(qū)。

重力侵蝕；起伏度；年均大雨天數(shù)；趨勢面；災(zāi)害

受復(fù)雜的自然環(huán)境和人為活動影響，我國水土流失問題嚴(yán)重，給國計民生和生態(tài)環(huán)境造成巨大損失[1]。土壤侵蝕可分為風(fēng)力侵蝕、水力侵蝕、凍融侵蝕和重力侵蝕。其中，重力侵蝕包括崩塌、滑坡、泥石流等，這些侵蝕類型都是典型的地質(zhì)災(zāi)害，造成的危害十分顯著。根據(jù)黃委會調(diào)查數(shù)據(jù)，重力侵蝕在很小的侵蝕面積下，可以產(chǎn)生占整個流域50%以上的侵蝕量[2]。因此，針對重力侵蝕災(zāi)害展開研究具有十分重要的意義[3-5]。重力侵蝕的影響因素可以分為兩大類，一類為準(zhǔn)備性因子，包括地形、地貌、地質(zhì)、植被等，第二類因子為觸發(fā)性因子，包括氣候、地震等[6]。其中，坡度是最為典型的地形因子，其與重力侵蝕間的關(guān)系研究也較多[7-8]，在觸發(fā)性因子中，降雨作為氣候的一項重要表征，是誘發(fā)重力侵蝕的另一個重要影響因子[9]。

重力侵蝕的發(fā)生同時受準(zhǔn)備性因子與觸發(fā)性因子的共同作用，單一重力侵蝕的發(fā)生具有相當(dāng)程度的隨機性，而在區(qū)域上則可表現(xiàn)出一定的規(guī)律性。為此，本研究以陜西省為例，選取準(zhǔn)備性因子中的地形與觸發(fā)性因子中的降雨作為重點，研究其與重力侵蝕分布特征的關(guān)系。

1 研究區(qū)域和數(shù)據(jù)介紹

1.1 陜西省概況

陜西省位于中國西北地區(qū)東部的黃河中游，面積為20.56萬km2。該省受大陸性氣候影響較為嚴(yán)重，干旱少雨，降水和干燥度從南向北減少[10]；地表由西北向東南逐級降低，全省劃分為6個地貌區(qū)：關(guān)中平原、風(fēng)沙區(qū)、黃土高原區(qū)、秦嶺、安康盆地及大巴山區(qū)[11]。

1.2 數(shù)據(jù)來源

在本研究中，收集了3個有關(guān)重力侵蝕的數(shù)據(jù)庫，分別是中國科技部、國家計委和國家經(jīng)貿(mào)委災(zāi)害研究組編的1949—1990年中國泥石流災(zāi)害逐年統(tǒng)計表，中國科技部、國家計委和國家經(jīng)貿(mào)委災(zāi)害研究組編的1990—2000年中國地質(zhì)災(zāi)害數(shù)據(jù)庫，中國地質(zhì)科學(xué)院地質(zhì)力學(xué)所整理的中國及鄰區(qū)地質(zhì)災(zāi)害數(shù)據(jù)庫。在這3個數(shù)據(jù)庫中，災(zāi)害記錄的時間跨度為1949—2000年，主要涉及的重力侵蝕包括滑坡、崩塌和泥石流。

地形數(shù)據(jù)選擇的是90m×90m的DEM，并以此數(shù)據(jù)為基礎(chǔ)，提取相關(guān)地形因子?？紤]到選取的DEM精度無法獲取準(zhǔn)確的坡度因子，故選取與坡度因子密切聯(lián)系的起伏度作為地形特征指標(biāo)。

選擇的降雨數(shù)據(jù)來自于國家氣象局統(tǒng)計整理的194個國際交換基本（準(zhǔn)）氣象站逐日氣候數(shù)據(jù)，選取日降雨≥25 mm大雨日數(shù)作為降雨指標(biāo)。在本研究中，還對植被因素進(jìn)行了考慮，用到的數(shù)據(jù)是黃河流域NDVI數(shù)據(jù)，該數(shù)據(jù)的分辨率為8 km×8 km。

1.3 數(shù)據(jù)處理

首先，從不同災(zāi)害數(shù)據(jù)庫中提取發(fā)生于陜西省的災(zāi)害數(shù)據(jù)，刪除其中重復(fù)部分，對于提取的重力侵蝕數(shù)據(jù)以縣市為單位進(jìn)行統(tǒng)計，共得到140次重力侵蝕事件。為了消除由于各個縣（含縣、縣級市，下同）面積不同造成重力侵蝕數(shù)據(jù)上的差異，將每個縣發(fā)生次數(shù)同面積相除，得到各個縣每千平方公里在1949—2000年間發(fā)生重力侵蝕次數(shù)。

2 結(jié)合地形和降雨因子的趨勢面分析

趨勢面分析方法可以綜合考慮兩個要素變化對于因變量的影響作用，在醫(yī)學(xué)、地貌、地質(zhì)勘探、氣候、測量等很多領(lǐng)域都得到了應(yīng)用[12-16]。在本研究中選取重力侵蝕發(fā)生次數(shù)作為因變量，而將起伏度與年均大雨天數(shù)作為自變量進(jìn)行趨勢面分析。

通過2次趨勢面和3次趨勢面的比較檢驗，發(fā)現(xiàn)當(dāng)次數(shù)提高到3次時，計算得到的檢驗值并不顯著，即較高次多項式對于回歸并無新的貢獻(xiàn)，因此選擇2次趨勢面分析方法，結(jié)果如圖1所示。其中橫坐標(biāo)表示地形因子（起伏度），縱坐標(biāo)表示降雨因子（年均大雨天數(shù)），而等高線表示的是每千平方公里發(fā)生重力侵蝕次數(shù)。

該趨勢面接近于馬鞍面，其鞍部位于圖中等高線最疏的地區(qū)（即圖中兩條虛線交叉的地區(qū)），而在x軸和y軸分別達(dá)到最高值時，形成了鞍面的兩個高值區(qū)。沿著接近對角線的直線變化，鞍面的斜率先變緩，經(jīng)過鞍面后，又逐漸增高。進(jìn)一步分析整個鞍面，可以將趨勢面分為4個區(qū)分別描述。

第Ⅰ區(qū)為重力侵蝕低發(fā)區(qū)，是以XY軸交點為中心的類似矩形的區(qū)域，其范圍為：起伏度0～31.57 m，大雨天數(shù)0～10 d。該區(qū)地形因子和降雨因子值在4個區(qū)內(nèi)是最低的，重力侵蝕發(fā)生次數(shù)在4個區(qū)內(nèi)也是最低的。在該區(qū)內(nèi)，地形因子和降雨因子同時增高，重力侵蝕發(fā)生次數(shù)也隨之增高。該區(qū)域等高線形狀類似一個橢圓形，當(dāng)從一條等高線躍遷至另一條線時，兩個因子同時增大所需因子增高值明顯低于一個因子增大所需值，說明當(dāng)?shù)匦魏徒涤晖瑫r增長對于重力侵蝕的促進(jìn)作用要比僅僅一個因子增長形成的促進(jìn)作用更強。

圖1 起伏度和大雨因子的趨勢面分區(qū)

第Ⅱ區(qū)為地形因子主控區(qū)，其范圍為：起伏度31.57～55.30 m，大雨天數(shù)0～10 d。在該區(qū)內(nèi)，地形因子達(dá)到較高的值，成為控制重力侵蝕的主要因素。在該區(qū)內(nèi)，降雨因子變化范圍較廣。從變化趨勢來看，當(dāng)降雨因子一定，地形因子增大時，即沿著X軸增大時，重力侵蝕發(fā)生次數(shù)是增高的，說明在一定降雨條件下，地形起伏越劇烈，重力侵蝕越容易發(fā)生，這一點比較容易解釋，但當(dāng)?shù)匦我蜃右欢〞r，隨著降雨因子的增大，重力侵蝕發(fā)生次數(shù)反而降低，且當(dāng)大雨天數(shù)達(dá)到10 d左右時，重力侵蝕發(fā)生次數(shù)反而降到了最低點，即到達(dá)了趨勢面的鞍部。一般認(rèn)為，滑坡發(fā)生的概率和滑坡數(shù)量與降雨量的大小成正相關(guān)關(guān)系[17]，所以這種“反?！钡淖兓厔菡f明趨勢面分析的方法由于僅僅考慮了兩個因子而產(chǎn)生的不足，因此，必須在此分析的基礎(chǔ)上引入對其他因子的考慮。

圖2 平均NDVI值隨大雨天數(shù)變化曲線

從圖2可以看出，隨著大雨天數(shù)的增高，平均NDVI值同時增高，即說明大雨天數(shù)高的地區(qū)，其植被蓋度越高。進(jìn)一步分析大雨天數(shù)與降雨量之間的關(guān)系，發(fā)現(xiàn)兩者的Pearson相關(guān)系數(shù)高達(dá)0.989（顯著性水平為0.01），即大雨天數(shù)也可以作為一個地區(qū)降雨狀況的指標(biāo)。即大雨天數(shù)高的地區(qū)，其降雨狀況良好，植被狀況要優(yōu)于干旱少雨地區(qū)。植被對于防治土壤侵蝕有重要作用，通過對降雨截留、緩沖作用以及落葉層對集中水流的分散和加糙作用大大減弱了坡面沖刷強度，減緩30°以上陡坡的切溝發(fā)育過程，從而阻礙重力侵蝕發(fā)生條件的形成；植被根系深入土壤，能夠承擔(dān)一部分抗剪強度，從而增強土體的抗變形能力[18]。因此，可以認(rèn)為：隨著降雨量的增高，降雨對于重力侵蝕的觸發(fā)作用與植被增高形成的抗侵蝕作用相互抵抗，最終形成了地形主控區(qū)這種特殊的重力侵蝕現(xiàn)象。

第III區(qū)為重力侵蝕高發(fā)區(qū)，其范圍為：起伏度31.57～55.30 m，大雨天數(shù)＞10 d。該區(qū)地形因子和降雨因子的值在4個區(qū)中是最大的，重力侵蝕發(fā)生次數(shù)很高。隨著地形因子和降雨因子的增高，重力侵蝕發(fā)生次數(shù)同時增高，其基本變化趨勢類似于第Ⅰ區(qū)，該區(qū)是第Ⅰ區(qū)在經(jīng)歷了趨勢面鞍部低值區(qū)后的延伸。當(dāng)?shù)匦我蜃雍徒涤暌蜃油瑫r達(dá)到很高的值以后，限制重力侵蝕發(fā)生的因素作用變得不顯著，從而沒有出現(xiàn)類似于第Ⅱ區(qū)的反?，F(xiàn)象。在該區(qū)，隨著因子的增高，等值線由疏變密，即說明隨著因子值的增大，對于重力侵蝕的促進(jìn)作用是一個加速的過程，而不是簡單的線性關(guān)系。第IV區(qū)在實際情況中并未出現(xiàn)，因此不予討論。分析陜西省各個縣市所屬類型區(qū)，可以發(fā)現(xiàn)，大部分縣屬于第Ⅰ類型區(qū)（共計66個縣），有24個縣屬于第Ⅱ類型區(qū)，只有3個縣（紫陽縣、鎮(zhèn)巴縣和嵐皋縣）屬于第III類型區(qū)，結(jié)果如圖3所示。

圖3 陜西省各縣市重力侵蝕分區(qū)

從地貌學(xué)的角度，該分區(qū)與陜西省的地貌分區(qū)吻合較好。其中，風(fēng)沙高原區(qū)、黃土高原區(qū)、關(guān)中盆地區(qū)大部分、漢中安康低山丘陵區(qū)部分和秦嶺山地區(qū)小部分屬于第Ⅰ類型區(qū)，秦嶺山地區(qū)大部分、漢中安康低山丘陵區(qū)大部分、大巴山地區(qū)一部分和關(guān)中平原區(qū)一部分屬于第Ⅱ類型區(qū)，大巴山地區(qū)大部分屬于第III類型區(qū)。

利用兩個要素進(jìn)行趨勢面分析的分區(qū)結(jié)果與地貌分區(qū)得到很好的吻合，說明不同的地貌類型區(qū)內(nèi)，控制重力侵蝕發(fā)生的要素是不同的，如秦嶺地區(qū)，其屬于地形因子主控區(qū)，在該區(qū)內(nèi)，地形因子已達(dá)到很高的值，因此，降雨狀況成為重力侵蝕發(fā)生與否的關(guān)鍵。進(jìn)一步說明，在該類區(qū)內(nèi)，需要時刻注意大型的降雨事件，因為該類事件很容易觸發(fā)重力侵蝕。而在屬于第Ⅰ類型區(qū)的黃土高原地區(qū)，監(jiān)測重力侵蝕事件需要同時注意降雨和地形兩個要素，從地形的角度，重力侵蝕主要集中于地形起伏較大的地區(qū)，這為監(jiān)測和預(yù)防提供了便利。

3 總結(jié)

本研究利用趨勢面分析的方法，對兩個要素同時變化下的重力侵蝕情形進(jìn)行了探討，得到相關(guān)結(jié)論如下：

（1）趨勢面分析可以得到四個區(qū)，分別為重力侵蝕低發(fā)區(qū)（Ⅰ區(qū)）、地形因子主控區(qū)（Ⅱ區(qū)）、重力侵蝕高發(fā)區(qū)（III區(qū)）和降雨因子主控區(qū)（IV區(qū)）。

（2）從第Ⅰ和III區(qū)的分析中可以發(fā)現(xiàn)，當(dāng)?shù)匦魏徒涤暌蜃油瑫r增加時，對于重力侵蝕具有顯著的促進(jìn)作用，且這種促進(jìn)作用不是簡單的線性關(guān)系。

（3）在第Ⅱ區(qū)中，當(dāng)?shù)匦我蜃右欢?，降雨因子增高時，出現(xiàn)反常的現(xiàn)象，重力侵蝕發(fā)生次數(shù)會隨著降雨因子的增高而降低。這一現(xiàn)象可能與植被覆蓋度隨降雨增多而增大的影響有關(guān)。

（4）將趨勢面分區(qū)同實際縣市相對應(yīng)，第Ⅰ區(qū)有66個縣，所占面積最大；次之為第Ⅱ區(qū)，有24個縣；而第III區(qū)只有3個縣。趨勢面分區(qū)與地貌分區(qū)吻合較好。

本研究相關(guān)結(jié)論可以為不同地貌類型區(qū)內(nèi)，監(jiān)測和預(yù)防重力侵蝕災(zāi)害事件提供服務(wù)。

[1]彭珂珊.水土流失是生態(tài)環(huán)境惡化的禍根[J].中國環(huán)境管理干部學(xué)院學(xué)報，2000，10（3、4）：90-98.

[2]黃河水利委員會黃河中游治理局.黃河水土保持志[M].鄭州：河南人民出版社，1993：57-61.

[3]呂玲莉.秦皇島市地質(zhì)災(zāi)害風(fēng)險區(qū)劃的思考與探索[J].中國環(huán)境管理干部學(xué)院學(xué)報，2010，20（6）：39-41.

[4]任改娟，孟凡杰，孫娟.老鷹巖滑坡基本特征及危害性分析[J].中國環(huán)境管理干部學(xué)院學(xué)報，2011，21（4）：40-42.

[5]劉亮，朗坤，李曉露.Midas/GTS在滑坡穩(wěn)定性分析及治理中的應(yīng)用[J].中國環(huán)境管理干部學(xué)院學(xué)報，2015，25（2）：44-46，81.

[6]DAI F C，LEE C F，NGAI Y Y.Landslide risk assessment and management：an overview[J].Engineering Geology，2002，64：65-87.

[7]FRANKSC A M.Characteristics of some rainfall-induced landslides on natural slopes，Lantau Island，Hong Kong[J].Quarterly Journal of Engineering Geology，1999，32：247-259.

[8]CHAU K T，SZE Y L，F(xiàn)UNGM K，et al.Landslide hazard analysis for Hong Kong using landslide inventory and GIS[J].Computers&Geosciences，2004，30（4）：429-443.

[9]CHEN J，DAI F C，YAO X.Holocene debris-flow deposits and their implications on the climate in the upper Jinsha River valley，China[J].Geomorphology，2008，93：493-500.

[10]楊武圣.陜西省志——氣象志[M].北京：氣象出版社，2001：90-92.

[11]陜西省地方志編纂委員會編.陜西省志——地理志[M].西安：陜西人民出版社，2000：292-313.

[12]彭建，楊明德，梁虹，等.路南巴江喀斯特流域空間趨勢面分析[J].地理與地理信息科學(xué)研究，2005，21（4）：109-112.

[13]胡光海，周興華.趨勢面分析在水深測量數(shù)據(jù)處理中的應(yīng)用[J].測繪工程，2004，13（3）：23-25.

[14]張瑞林，李東印，李小軍.趨勢面擬合法在煤礦瓦斯地質(zhì)變量預(yù)測中的應(yīng)用研究[J].中國安全科學(xué)學(xué)報，2006，16（12）：25-29.

[15]李隨民，姚書振，韓玉丑.Surfer軟件中利用趨勢面方法圈定化探異常[J].地質(zhì)與勘探，2007，43（2）：72-75.

[16]王曉燕，沈毅，陳坤，等.趨勢面分析法在肺癌死亡率地理分布研究中的應(yīng)用[J].中華流行病學(xué)雜志，2007，28（6）：608-612.

[17]張玉成，楊光華，張玉興.滑坡的發(fā)生與降雨關(guān)系的研究[J].災(zāi)害學(xué)，2007，22（1）：82-85.

[18]薛海，王文成，何沛華.黃河中游地區(qū)重力侵蝕影響因素分析[J].人民黃河，2008，30（5）：60-62.

（編輯：周利海）

Gravitational Erosion Disaster Distribution Characteristics of Shaanxi Province

LiWeibo1,Hou Kebin2
（1.Environmental Impact Assessment and Monitoring Centre for China Academy of Launch Vehicle Technology,Beijing Aerospace Institute for Metrology and Measurement Technology,Beijing 100076,China；2.Beijing Municipal Appraisal Centre of Environmental Impact Assessment,Beijing 100161,China）

The distribution characteristics of gravitational erosion in Shaanxi province was studied based on gravitational erosion hazard data of Shaanxi province.Waviness and average annual heavy rain day(AAHRD)were chosen as landform and rainfall factors.As the results showed,gravitational erosion in Shaanxi could be divided into 3 zones according to dominant factor,and the 3 zoneswere:gravitational erosion low occurrence frequency area(ZoneⅠ),landform factor controlling area(ZoneⅡ),and gravitational erosion high occurrence frequency area(Zone III).Comparison wasmade between gravitational erosion zoning and geomorphic division,and good correlation between the two was found:ZoneⅠincluded Guanzhong Basin,Aeolian Sandy Plateau and the Loess Plateau,ZoneⅡincluded Qinling and Ankang Basin,and ZoneIII concentrated in Daba Mountain Region.

gravitational erosion,waviness,average annual heavy rain day,trend surface,disaster

X43

A

1008-813X（2017）02-0033-04

10.13358 /j.issn.1008-813x.2017.02.09

2017-02-15

李衛(wèi)波（1983-），女，山東濰坊人，畢業(yè)北京大學(xué)環(huán)境科學(xué)專業(yè)，碩士，工程師，主要從事環(huán)境影響評價、環(huán)境監(jiān)測的研究工作。