基于DEM的河網(wǎng)結構與地貌特征相關性研究

摘 要：以數(shù)字高程模型數(shù)據(jù)為基礎，運用ArcGIS10.2軟件，通過水文分析模塊提取河網(wǎng)信息，研究河網(wǎng)結構與區(qū)域構造、地形地貌的響應關系。結果表明，河流總體流向與斷層走向相同，河網(wǎng)分布、干流流向與區(qū)域斷層走向具有一致性；實驗得到5級河流，一級河流最多，幼年期、壯年期和老年期流域面積分別約占總面積的15.29%、74.85%、9.86%，以壯年期為主，坡度大于15°的流域面積占總面積的77%；干流區(qū)地形起伏度較大，支流集中分布在平地地區(qū)。

關鍵詞：數(shù)字高程模型；河網(wǎng)；構造；地貌特征

河網(wǎng)分布特征與礦產(chǎn)資源、地質災害等資源環(huán)境研究具有重要意義。數(shù)字高程模型（DEM）作為地形地貌數(shù)字化表達方式，成為與空間地理信息相關研究的必然基礎。同時數(shù)字高程模型數(shù)據(jù)在流域方向的研究越來越受到重視[1]。基于目前對大區(qū)域內河網(wǎng)空間分布特征及其與構造活動、地貌特征的響應關系的研究較少。文章以SRTM-DEM（Shuttle Radar Topography Mission-DEM）為數(shù)據(jù)源，利用ArcGIS10.2軟件的水文分析功能提取研究區(qū)內河網(wǎng)結構，結合流域內構造地形地貌特征，定量地分析西藏境內瀾滄江和怒江流域河流分布與地貌、構造運動之間的響應關系。

1 區(qū)域概況

西藏東部地區(qū)平均海拔3500m以上，以瀾滄江和怒江流經(jīng)的昌都縣、類烏齊縣、八宿縣、察雅縣、左貢縣和芒康縣為研究區(qū)，區(qū)內日照時間長，輻射強、晝夜溫差大；干濕分明，夏季多夜雨，冬春多風，年平均氣溫7.6℃，年降雨量400～600mm。由于水系的切割作用，形成了非常顯著的垂直地帶性特征，有著復雜的地貌結構、氣候類型及植被狀況。

2 研究方法

2.1 水系結構提取與分級

利用水文分析模塊進行洼地填充、水流方向、設置不同匯流閾值，提取研究區(qū)內部的地表水系網(wǎng)絡。水系級別的定義采用現(xiàn)今最為廣泛流行的Strahler分類系統(tǒng)，將所有河網(wǎng)弧段中沒有支流的河網(wǎng)弧段定為第1級，兩個1級河網(wǎng)弧段合成的河網(wǎng)弧段為第2級，如此下去為第3級，第4級，……，一直到河網(wǎng)出水口[2]。

2.2 面積高程積分

面積高程積分（Hypsometric Integral，HI） 是通過統(tǒng)計流域地表的高程組合信息，從而揭示流域地貌形態(tài)與發(fā)育特征的重要指標。Pike于1971年提出以流域高程起伏比作為面積高程積分的簡易算法[3，4]。計算公式為：

HI=（Hmean-Hmin）/（Hmax-Hmin） （1）

式中，Hmean、Hmin、Hmax分別是區(qū)域高程的均值、最小值、最大值。

按照戴維斯地貌旋回理論，當流域高程-面積積分值大于0.60 時，地貌特征變化迅速，水系不斷擴展分支，流域侵蝕劇烈，是地貌發(fā)育的不均衡階段，稱為幼年期。當高程-面積積分值小于0.60大于0.35時，是地貌發(fā)育的壯年期；高程-面積積分值小于0.35時，為老年期[5，6]。

3 結果與分析

3.1 河網(wǎng)提取

河流閾值的確定直接影響河流數(shù)量的提取和分析，通過ArcGIS10.2水文分析模塊進行洼地填充、水流方向、匯流累積量的求取之后，漸變設置匯流閾值，計算Stream order、Stream Link及Watershed，以1000為基礎、步長200累計統(tǒng)計不同匯流閾值下的溝谷段數(shù)、溝谷總長，并計算流域內的溝壑密度等水文信息。

由統(tǒng)計信息可知，隨著匯流閾值越大，得到研究區(qū)溝谷段數(shù)、溝谷總長和溝壑密度值越小。當匯流閾值小于3000時，流域內溝壑密度變化較大，當閾值大于3000時，溝壑密度變動較小，并趨于平穩(wěn)；當匯流閾值為1000時，溝谷段數(shù)為3480，約為閾值為5600時對應溝谷段數(shù)的5.6倍，而溝壑總長是其2.3倍，可知流域閾值的選取對提取的流域溝谷段數(shù)影響較大。選擇研究區(qū)合適的匯流閾值，是流域河網(wǎng)提取與特征分析的基礎，當匯流閾值為5600時，溝谷段數(shù)為626，溝谷總長為5305.85km，溝壑密度為0.23km/km2，當匯流閾值增加到6000時，溝谷段數(shù)和總長變動較小并趨于平穩(wěn)，溝壑密度保持不變，數(shù)據(jù)顯示，確定最佳匯流閾值為5600。

由最佳閾值提取研究區(qū)河網(wǎng)數(shù)據(jù)，得到共5級水系結構，612條河流。其中一級河流307條，二級河流108條，三級河流70條，四級河流91條，五級河流36條。一級河流最多，基本覆蓋整個流域范圍，一二級河流占研究區(qū)的68%。

3.2 河網(wǎng)分布與構造關系

水系發(fā)育程度和流向受區(qū)域地形地貌、構造特征和地表巖性的影響。探究研究區(qū)河網(wǎng)分布特征與區(qū)域斷裂構造的相關性，對河流和斷層線進行玫瑰圖式分析，以15°為一個間隔，統(tǒng)計各個方向上的河數(shù)量和斷裂構造數(shù)量。斷層線性構造走向方位統(tǒng)計圖，揭示斷層主要沿北西向，315°方向上斷層較多。河流主要流向也是北西向，在315°和330°方向上發(fā)育的河流較多，與斷層的主要走向一致，研究區(qū)在北東方向上的流河發(fā)育也較為明顯，這與區(qū)域構造45°的斷層相關。對比分析河網(wǎng)數(shù)據(jù)與區(qū)域構造特征，得到河流流向與區(qū)域構造具有明顯的相關性，河網(wǎng)分布受構造活動的強烈影響，河流流向和支流的發(fā)育與區(qū)域斷層構造方向具有一致性。

3.3 流域內地形地貌特征

地貌特征是水系結構的基礎，研究區(qū)最低高程為1666米，最大高程為6468米，地勢起伏較大。1666～3000m區(qū)域占總面積的22.67%；3000～5000m的區(qū)域最大，占總面積的64.63%；5000～6468m的區(qū)域最小，僅占總面積的12.7%。

地面坡度的級別不同，對于坡地系統(tǒng)的發(fā)育影響不同，對于水系結構的影響也不一樣。研究區(qū)坡度范圍在0～75°之間，可根據(jù)其對水土流失的影響和所發(fā)揮的水土保持功能等，將研究區(qū)的坡度劃分為4個范圍等級：0～8°、8～15°、15～25°，以及>25°。其中坡度小于8°區(qū)域最少，占總面積的8.09%；坡度在8～15°之間占15.19%；坡度在15～25°之間的占30.68%；坡度>25°的占46.04%。可見坡度主要集中在15°以上。

研究區(qū)地面起伏度最大值為837米，將研究區(qū)地面起伏度<200米的確定為平地；地面起伏度在200～500米之間的確定為丘陵；>500米的確定為山地。研究區(qū)主要以平地地形為主，占總面積的74.67%。

研究區(qū)四級河流和五級河流流域地表起伏度較大，河流交匯的地區(qū)地形起伏度較大，以丘陵地形為主；一級河流和二級河流主要分布在平地地區(qū)，地表起伏度較小；研究區(qū)內怒江流域的地表起伏度比瀾滄江流域地表起伏度大，由于受地形影響，瀾滄江干流右側河流發(fā)育較多，左側支流發(fā)育較少，以一級河流為主。

匯流累積量對面積高程積分值的影響不大，基于最佳閾值提取區(qū)域內河流的面積高程積分值。得到研究區(qū)處于幼年期地貌河流86條，主要集中分布在瀾滄江和怒江干流左側1級流域，占流域面積的15.29%；處于壯年期地貌河流416條，主要分布在瀾滄江和怒江干流外側，各級河流都有分布，2級河流較多，占流域面積的74.85；處于老年期地貌特征河流110條，主要分布在干流流域范圍內，占流域面積的9.86%；研究區(qū)以壯年期流域地貌為主。瀾滄江與怒江干流呈平行分布。

4 結束語

（1）基于ArcGIS10.2水文分析，提取研究區(qū)流域內河網(wǎng)分布特征并進行面積高程積分分析。匯流閾值越大，河道數(shù)目越少，河道總長和河網(wǎng)密度均呈現(xiàn)持續(xù)下降趨勢；當匯流閾值為5600時，河網(wǎng)密度為0.23km/km2，趨于平穩(wěn)；不同匯流閾值產(chǎn)生的河網(wǎng)形態(tài)各有差異，閾值越大，生成河網(wǎng)越稀疏，反之越密集。水系劃分為5級結構，河網(wǎng)總數(shù)目626條，總長度為5305.85千米。面積-高程積分結果表明，研究區(qū)內大多數(shù)河流處于發(fā)育壯年期，幼年期、壯年期和老年期流域面積分別約占總面積的15.29%、74.85%、9.86%。

（2）研究區(qū)相對高程較大，高程值在3000～5000m之間的分布較廣，坡度主要集中在15°以上，河流交匯地區(qū)以丘陵地形為主，地形起伏度較大，支流主要分布在平地地區(qū)；干流流域內地形起伏度較大，河流發(fā)育受區(qū)域地貌特征影響較大，研究區(qū)內兩條干流呈平行狀發(fā)展。區(qū)域內河流發(fā)育與斷層等構造活動方向一致，河流侵蝕方向與區(qū)域構造、地貌特征具有明顯相關性。

（3）影響流域發(fā)展的因素是多方面的，流域地貌發(fā)育和河網(wǎng)分布與地表基巖巖性有關，在后續(xù)實驗中將進一步完善。

參考文獻

[1]劉新華，楊勤科，湯國安.中國地形起伏度的提取及在水土流失定量評價中的應用[J].水土保持通報，2001，21（1）：57-59.

[2]劉興旺，袁道陽，史志剛，等.六盤山斷裂帶構造活動特征及流域盆地地貌響應[J].地震工程學報，2015，37（1）：168-174.

[3]劉飛，范建容.地貌特征與土壤侵蝕空間分布直方圖相似度研究[J].水土保持研究，2011，18（6）：18-21.

[4]程三友，李英杰.撫仙湖流域地貌特征及其構造指示意義[J].地質力學學報，2010，16（4）：383-392.

[5]孫希華，姚孝友，周虹，等.基于DEM的山東沂沭泗河流域地貌演化與水土流失研究[J].水土保持通報，2005，25（4）：24-28.

[6]祝士杰，湯國安，李發(fā)源，等.基于DEM的黃土高原面積高程積分研究[J].地理學報，2013，68（7）：921-932.

作者簡介：陳娟（1992-），女，重慶市，成都理工大學，碩士研究生，主要從事3S應用研究。