周文龍，趙衛(wèi)權(quán)，楊家芳，李 威，李一兵，張 凡，蘇維詞

（貴州省山地資源研究所，貴州 貴陽550001）

1 引言

數(shù)字高程模型 （Digital Elevation Model，簡 稱DEM）已經(jīng)被廣泛地應用到諸如地形特征提取、流域水系分析等眾多數(shù)字地形分析領(lǐng)域［1～5］。以DEM為基礎可提取多種地形特征，如坡度、坡向、坡面形態(tài)、流域邊界、水流路徑、溝壑密度［6～9］等，這些特征在地理信息系統(tǒng)的支持下均可用圖形和屬性數(shù)據(jù)來表示。但低精度的DEM將導致水文地形等信息受損，這將嚴重影響流域水文模型參數(shù)的確立及水文過程模擬分析的精度［10］。因此目前針對DEM的誤差來源［11］、誤差檢測［12］、誤差可視化 以及誤差校正［13］、DEM 及數(shù)字地形分析中的尺度問題［14～16］等方面有大量的研究成果。

水文水系數(shù)據(jù)是地理科學數(shù)據(jù)的核心之一，作為地學領(lǐng)域數(shù)據(jù)的一個重要組成部分，是支撐相關(guān)地學研究的基礎性數(shù)據(jù)。而基于DEM的溝谷網(wǎng)絡提取由于受DEM精度、地形復雜度、溝谷網(wǎng)絡的結(jié)構(gòu)特征及模擬方法的制約［17］，存在很大的隨意性和主觀性［18］。Martz和 Garbrecht認為［19］，集水面積閾值具有一定的地貌局限性，錯誤地使用其它地貌類型適用的閾值，將會得到大量的偽水道，直接影響到生成的水系、溝壑的詳細程度以及子流域的劃分數(shù)目［20］。Tribe指出解決這一問題的方法是考慮地貌參數(shù)的空間變化［21］，不同類型的區(qū)域設置不同的閾值，或者考慮反映地貌空間差異性的其它參數(shù)。

喀斯特地貌具有獨特的地貌空間分異特征，集中體現(xiàn)在其內(nèi)部地貌形態(tài)的有序組合。根據(jù)對喀斯特地區(qū)水系結(jié)構(gòu)的研究，發(fā)現(xiàn)其地表、地下水系結(jié)構(gòu)均符合Horton流水網(wǎng)結(jié)構(gòu)［22～25］，這說明喀斯特地區(qū)兩套水系也和世界上大多數(shù)水系一樣存在自相似性。喀斯特區(qū)地表水系的分維不僅是反映水系形態(tài)特征的參數(shù)，也是表征流域內(nèi)喀斯特發(fā)育程度的量化指標［26］。同時有研究表明：在流域地貌演化的不同階段，切割深度與溝壑密度的關(guān)系不同。故溝壑是一個綜合性很強的地貌指標，既能反映地面破碎程度與土壤侵蝕強度，又能反映地貌的演化階段［27］。因此，本文擬結(jié)合貴州喀斯特高原不同地貌類型區(qū)，研究集水面積閾值隨喀斯特地貌復雜度變化而變化的規(guī)律，并研究匯流閾值與溝谷網(wǎng)絡總長度的量化關(guān)系，嘗試利用DEM提取溝壑特征來分析貴州水系、地貌的空間分異規(guī)律，協(xié)助推斷貴州喀斯特流域地貌的整體演化階段。研究結(jié)果將對于完善喀斯特高原溝谷網(wǎng)絡自動提取的理論與方法，對于深化人們對喀斯特高原溝谷網(wǎng)絡空間分異的認識，以及對水文過程分析、土壤侵蝕計算以及土地利用規(guī)劃都有非常重要的現(xiàn)實意義。

2 材料與方法

2.1 數(shù)據(jù)來源與軟件平臺

貴州省DEM數(shù)據(jù)采用美國地質(zhì)調(diào)查局（The U.S.Geological Survey，簡稱 USGS）發(fā)布的7.5′分帶1∶24萬數(shù)字高程模型（地面分辨率90m），借助ArcGIS 9.3軟件平臺中的水文分析工具（Hydrology）對貴州省境內(nèi)溝壑信息進行智能提取，并基于貴州省八大流域邊界及行政區(qū)劃界線運用疊加分析對溝壑信息進行分區(qū)計算。同時參照貴州省1∶25萬基礎地理信息資料中的水文數(shù)據(jù)對計算結(jié)果進行對比驗證，并結(jié)合2005年貴州省石漠化空間分布數(shù)據(jù)來研究喀斯特不同地貌類型與水系發(fā)育之間的關(guān)系。為便于疊加分析，將上述三套數(shù)據(jù)投影方式均統(tǒng)一為：Albers等積圓錐投影，中央經(jīng)線為105°E，雙標準緯線分別為25°N和47°N，采用Krasovsky橢球體。

2.2 計算方法與流程

常規(guī)溝壑密度的制作流程主要包括水流方向的確定、洼地填平、匯流累積矩陣的生成以及溝谷網(wǎng)絡的生成，最后計算溝谷密度［28］，具體操作及計算結(jié)果如下（圖1）。

（1）在ArcMap中加載原始DEM數(shù)據(jù)。

（2）啟動 ArcToolbox，展開 Spatial Analyst Tools工具箱，打開Hydrology工具集。

（3）洼地填充：雙擊 Hydrology工具集中的Fill工具，進行原始DEM的洼地填充。

（4）無洼地的水流方向的計算：雙擊Hydrology工具集中的Flow Direction工具，對步驟（3）中的填充過的無洼地DEM數(shù)據(jù)進行水流方向計算。

（5）匯流累積量的計算：雙擊Hydrology工具集中的Flow Accumulation工具，選擇輸入步驟（4）中水流方向計算結(jié)果數(shù)據(jù)進行匯流量的計算。

（6）柵格河網(wǎng)的生成。在柵格河網(wǎng)的生成中，需設置一個匯流累計閾值（本研究分別用500和1000提取）。雙擊Spatial Analyst Tools→Map Algebra→Multi Output Map工具。在文本框中鍵入如：F：＼溝壑密度＼Gansu＼str1000＝con（F：＼溝壑密度＼Gansu＼flowacc＞1000，1），得到柵格河網(wǎng)數(shù)據(jù)str1000。

（7）柵格河網(wǎng)矢量化：雙擊Hydrology工具集中的Stream To Feature工具。選擇步驟（6）中的str1000作為河網(wǎng)輸入數(shù)據(jù)，將步驟（4）中計算結(jié)果數(shù)據(jù)作為水流方向輸入數(shù)據(jù)。

（8）溝壑密度計算。在步驟（7）中所得的矢量化數(shù)據(jù)表中添加字段LENGTH，右鍵此字段選擇Calculate Geometry進行溝壑長度計算，任務完成后再右鍵此字段選擇Statistics進行溝壑長度統(tǒng)計求和。這樣樣區(qū)的溝壑密度就是：Ds＝ ∑L/A。

3 結(jié)果與討論

3.1 貴州省石漠化空間分布與河網(wǎng)水系的關(guān)系分析

喀斯特地區(qū)獨特的地表地下空間二元結(jié)構(gòu)導致“地表水貴如油，地下水滾滾流”的工程性缺水現(xiàn)象明顯，而貴州省喀斯特石漠化發(fā)育，是喀斯特地貌類型演化的直接表征，地表水資源的科學配置成為石漠化綜合治理的關(guān)鍵。因此，對貴州省的石漠化分布與河網(wǎng)水系之間的關(guān)系研究顯得尤為重要（圖2）。本研究結(jié)合貴州省2005年石漠化空間分布數(shù)據(jù)，運用ArcMap進行空間疊置分析，分別統(tǒng)計不同等級石漠化面積中的河網(wǎng)長度數(shù)據(jù)（表1），再進行相關(guān)分析得出：貴州省不同等級石漠化面積與貴州省1∶25萬河網(wǎng)水系長度呈顯著正相關(guān)關(guān)系，關(guān)系式為y＝0.5487x－2634.9（R2＝0.9223，n＝6）。由此可以看出喀斯特石漠化程度是影響喀斯特地區(qū)地表水系發(fā)育的重要因子，石漠化嚴重地區(qū)喀斯特漏斗、地下河等發(fā)育，地表河流大部分轉(zhuǎn)由地下流走，導致地表河網(wǎng)密度低。因此，要真實反映喀斯特地區(qū)的地表起伏切割狀況，除了考慮地表顯現(xiàn)的河流，還必須要考慮溝壑狀況，即采用溝壑密度更為合理。

表1 貴州省不同等級石漠化面積與河網(wǎng)長度統(tǒng)計

3.2 貴州省八大流域河網(wǎng)密度及溝壑密度的一般統(tǒng)計特征

無論是河網(wǎng)密度還是溝壑密度均與流域劃分有關(guān)，即反映的均是集雨面的匯流狀況。本研究首先采用貴州省劃分的八大流域為界線，分析自然環(huán)境分區(qū)下不同二級流域河網(wǎng)密度與溝壑密度之間的關(guān)系。采用ArcMap中的空間疊置分析工具分別對河網(wǎng)水系及溝壑分布數(shù)據(jù)進行統(tǒng)計分析，并對計算出的河網(wǎng)密度和溝壑密度（閾值500和閾值1000）進行相關(guān)分析（圖3）。

經(jīng)相關(guān)性檢驗分析發(fā)現(xiàn)，采用二級流域為單元的貴州省1∶25萬水系圖中的河網(wǎng)密度與USGS數(shù)據(jù)智能提取的溝壑密度不存在相關(guān)性，閾值500與閾值1 000分別提取的溝壑密度值之間呈顯著相關(guān)，相關(guān)關(guān)系為y＝0.9685x＋0.1899（R2＝0.8445，n＝8）（表2）。智能提取的溝壑特征為均勻分布，干溝信息也能較為充分體現(xiàn)出來。通過對數(shù)據(jù)的比對，當閾值設定為1 000時，溝壑密度與河網(wǎng)密度值相當，但并不能很好地體現(xiàn)溝壑情況；當閾值設定為500時，可基本反映地表溝壑情況（圖4）。

表2 基于貴州省八大流域分區(qū)的河網(wǎng)密度和溝壑密度特征統(tǒng)計

3.3 貴州省行政區(qū)劃河網(wǎng)密度及溝壑密度的一般統(tǒng)計特征

同上，以貴州省各地、州、市等行政界線為劃分單元，分析不同行政單元河網(wǎng)密度與溝壑密度之間的關(guān)系。經(jīng)相關(guān)性檢驗分析發(fā)現(xiàn)，采用行政單元的貴州省1∶25萬水系圖中的河網(wǎng)密度與USGS數(shù)據(jù)智能提取的溝壑密度同樣不存在相關(guān)性，閾值500與閾值1 000分別提取的溝壑密度值之間呈顯著相關(guān)，相關(guān)關(guān)系為y＝1.0658x＋0.1494（R2＝0.9005，n＝9）（表3）。

表3 基于貴州省行政區(qū)劃的河網(wǎng)密度和溝壑密度特征統(tǒng)計

3.4 貴州省不同分區(qū)河網(wǎng)密度和溝壑密度綜合比較分析

通過綜合分析上述兩種分區(qū)劃分結(jié)果可以發(fā)現(xiàn)：貴州省1∶25萬水系圖中的河網(wǎng)密度與USGS數(shù)據(jù)智能提取的溝壑密度不存在相關(guān)性，閾值500與閾值1 000分別提取的溝壑密度值之間呈顯著相關(guān)，相關(guān)關(guān)系為y＝1.0787x＋0.892（R2＝0.892，n＝17）。貴州省1∶25萬水系圖能很好地表達地表河流情況，即與貴州的水系順應地形大勢由西、中部呈帚狀向北、東、南三面分流，八大水系侵蝕切割著高原主體，區(qū)域地貌及地貌類型組合也隨之不同，表現(xiàn)出有規(guī)律分布的省情相符。但河網(wǎng)密度未能很好地表現(xiàn)干溝以及喀斯特地區(qū)地下河信息，溝壑密度則能在一定程度上彌補這一缺陷，溝壑密度可以作為單因子來刻畫地表起伏狀況。實際水系圖中因喀斯特地貌因素表現(xiàn)出中、強度石漠化區(qū)地表河流稀疏，非喀斯特區(qū)則和溝壑密度有較高的耦合度。如南盤江和紅水河水系所屬的黔西南州，是石漠化發(fā)育最嚴重的地區(qū)之一；沅江和柳江水系值明顯較高，正是黔東南州所在的非喀斯特區(qū)，河網(wǎng)密度高（圖5）。

4 結(jié)論

（1）貴州省不同等級石漠化空間分布與貴州省1∶25萬河網(wǎng)水系分布呈顯著正相關(guān)關(guān)系，喀斯特石漠化程度是影響貴州省地表河網(wǎng)水系發(fā)育的重要因子，因此，地表水資源的科學配置成為石漠化綜合治理的關(guān)鍵。

（2）貴州省1∶25萬水系圖中的河網(wǎng)密度與USGS數(shù)據(jù)智能提取的溝壑密度不存在相關(guān)性，因喀斯特地區(qū)獨特的地表地下二元空間結(jié)構(gòu)以及石漠化廣布，實際水系圖中表現(xiàn)出中、強度石漠化區(qū)地表河流稀疏，而溝壑密度較高的不對稱現(xiàn)象；非喀斯特區(qū)則和溝壑密度存在較高的耦合度。

（3）溝壑密度能在一定程度上彌補河網(wǎng)密度未能體現(xiàn)干溝信息等的缺陷，它可以作為單因子來刻畫地表起伏狀況。

（4）在計算研究貴州省溝壑密度匯流閾值與溝谷網(wǎng)絡總長度的量化關(guān)系時，分別設定閾值為1 000和500，通過數(shù)據(jù)比對發(fā)現(xiàn)：當閾值設定為1 000時，溝壑密度與貴州省1∶25萬河網(wǎng)密度值相當，但并不能很好地體現(xiàn)溝壑情況；當閾值設定為500時，可基本反映地表溝壑情況。但最佳閾值還有待進一步研究確定。

（5）可進一步細化研究單元，可從縣域或小流域尺度來定量研究河網(wǎng)密度與溝壑密度之間的關(guān)系，并提高DEM數(shù)據(jù)精度，從不同比例尺數(shù)據(jù)來對比驗證最佳匯流閾值。另外，喀斯特石漠化綜合治理采取以小流域為治理單元，可利用更高精度數(shù)據(jù)對貴州省小流域進行劃分，為石漠化綜合治理小流域確定提供科學依據(jù)。

［1］周啟鳴，劉學軍.數(shù)字地形分析［M］.北京：科學出版社，2006.

［2］李 麗，郝振純.基于DEM的流域特征提取綜述［J］.地球信息進展，2003，18（2）：251～256.

［3］謝順平，都金康，羅維佳，等.基于DEM的復雜地形流域特征提?。跩］.地理研究，2006，25（1）：96～102.

［4］李 碩，曾志遠，張運生.數(shù)字地形分析技術(shù)在分布式水文建模中的應用［J］.地球科學進展，2002，17（5）：769～774.

［5］畢華興，譚秀英，李笑吟.基于DEM的數(shù)字地形分析［J］.北京林業(yè)大學學報，2005，27（2）：49～53.

［6］孫崇亮，王卷樂.基于DEM的水系自動提取與分級研究進展［J］.地理科學進展，2008，27（1）：118～124.

［7］朱紅春，劉海英，張繼賢，等.基于DEM的流域地形因子提取與量化關(guān)系研究——以陜北黃土高原的實驗為例［J］.測繪科學，2007，32（2）：138～140.

［8］李 俊，湯國安，張 婷，等.利用DEM提取陜北黃土高原溝谷網(wǎng)絡的匯流閾值研究［J］.水土保持通報，2007，27（4）：153～157.

［9］唐從國，劉叢強.基于 Arc Hydro Tools的流域特征自動提取——以貴州省內(nèi)烏江流域為例［J］.地球與環(huán)境，2006，34（3）：30～36.

［10］湯國安，趙牡丹，曹 菡.DEM地形描述誤差空間結(jié)構(gòu)分析［J］.西北大學學報：自然科學版，2000，30（4）：349～352.

［11］Felicisirno A.M.A parametric statistical method for error detection in digit al elevation models［J］.ISPRS Journal of Photogrammetry and Remote Sensing，1994，49（4）：29～33.

［12］Hannah M.J.Error detection and correction in digital terrain models［J］.Photogrammetric Engineering and Remote Sensing，1981，47（1）：63～69.

［13］劉學軍，盧華興，仁 政，等.論DEM地形分析中的尺度問題［J］.地理研究，2007，26（3）：433～442.

［14］湯國安，劉學軍，房 亮，等.DEM及數(shù)字地形分析中尺度問題研究綜述［J］.武漢大學學報：信息科學版，2006，31（12）：1 059～1 066.

［15］湯國安，楊勤科，張 勇，等.不同比例尺DEM提取地面坡度的精度研究——以在黃土丘陵溝壑區(qū)的試驗為例［J］.水土保持通報，2001，21（1）：53～56.

［16］郎玲玲，程維明，朱啟疆，等.多尺度DEM提取地勢起伏度的對比分析——以福建低山丘陵區(qū)為例［J］.地球信息科學，2007，9（6）：1～6.

［17］劉 光，李樹德，張 亮.基于DEM的溝谷提取算法綜述［J］.地理與地理信息科學，2003，19（5）：11～15.

［18］楊錦玲.基于分形的數(shù)字水系集水面積閾值確定研究［J］.測繪科學，2011，36（4）：33～34.

［19］Martz W，Garbrecht J.Short communication，automated recog－nition of valley lines and drainage networks from grid digital elevation models：a review and a new method Comment［J］.Journal of Hydrology，1995（167）：393～396.

［20］秦福來，王曉燕，王麗華，等.DEM流域特征提取及其在非點源污染模擬中的應用［J］.地理空間信息，2006，4（2）：46～49.

［21］Tribe A.Automated recognition of valley lines and drainage networks from grid digital elevation models：a review and a new method Comment［J］.Journal of Hydrology，1992（139）：263～293.

［22］楊明德.論貴州巖溶水賦存的地貌規(guī)律［J］.中國巖溶，1982（2）：81～91.

［23］V.R.Baker.Geomorphology and Hydrology of Karst 15Drain－age Basins and Cave Channel Networks in central New York［J］.Water Resources Res.1973（9）：367～384.

［24］梁 虹，盧 娟.喀斯特流域水系分形、熵及其地貌意義［J］.地理科學，1997，17（4）：310～315.

［25］梁 虹，楊明德.喀斯特流域水文地貌系統(tǒng)及其識別方法初探［J］.中國巖溶，1994，13（1）：1～9.

［26］許 模，王 迪，漆繼紅，等.基于分形理論的喀斯特地貌形態(tài)分析［J］.成都理工大學學報：自然科學版，2011，38（3）：328～333.

［27］張麗萍，馬志正.流域地貌演化的不同階段溝壑密度與切割深度關(guān)系研究［J］.地球信息進展，1998，17（3）：273～278.

［28］湯國安，楊 昕.ArcGIS地理信息系統(tǒng)空間分析實驗教［M］.北京：科學出版社，2006.