宋桃云

（北京中水科水電科技開發(fā)有限公司，北京 100038）

地理信息技術（GIS）成為了地理學研究的熱門語言，地理數據是地理信息系統(tǒng)的重要組成部分。DEM流域提取過程是獲取地理數據的一種方式，可獲得豐富的水文信息，能夠反映地形地貌特征[1]。因此，被廣泛應用于景觀設計、流域仿真、防洪災害、測繪制圖等領域，是測繪學、地理學、地貌學研究的熱點。坡面徑流提取河網運用廣泛，學者對精度提高了要求，通過改進處理方法進一步優(yōu)化流域模擬效果[2]。國內相關研究應用相對較晚，但20世紀90年代以來，DEM自動獲取流域特征越來越高效快捷，進入流域數字化迅速發(fā)展階段。

1 DEM預處理

流域水文分析是建立在無洼地上的DEM數據提取水文特征，但是實際操作中普遍存在小平原和凹陷洼地，洼地的存在會阻礙自然水流流向，導致提取發(fā)生錯誤。為了確保提取的流域水系是連續(xù)的，要在流域提取之前對洼地進行填充預處理。

Wang等[2]通過降低洼地邊緣的柵格單元高程進行改進，提出了洼地裂開法。Lindsay等[3]對洼地填充和洼地裂開法進行比較，得到快速有效的處理方式?？追舱艿萚4]認為由于輸入數據質量低、DEM生成過程中插值誤差引起平坦區(qū)和封閉洼地，通過添加高程信息來消除，使得流場結構匹配精確。為避免傳統(tǒng)填洼方法進行增高處理，李輝等[5]提出一種三方向搜索的洼地處理方法，主要通過對洼地和周圍網格高程值進行增減以達到目的。曾向輝等[6]通過在J＆D和M＆V算法基礎上建立改進的W＆L算法，對填洼后DEM升序排序，逐步對每一柵格進行計算。從上述研究中看出，對洼地和小平地的處理方法不同直接影響到流域提取效果，雖對算法有所改進，但不確定算法的通用性。

2 流向判斷

地表徑流總是從高地勢向低流，最后經流域出口匯流。常用的單流向法是D8法，通過計算中心柵格單元與相鄰8個柵格的最大距離落差值來確定。若差為正值，則為流出；否則為流入。

單流向法不適合模擬水流在坡面上漫流情況，會出現(xiàn)大量平行流現(xiàn)象，接著對多流向算法進行研究，有混合算法、矢量柵格算法等。秦承志等[7]在基于柵格DEM的多流向算法述評中，根據水流分配策略，將多流向算法分為4類，水流分配權重隨局域地形特征變化的多流向法精度更高。他提出的（MFDfg）法[8]，水流分配隨著下坡坡度變化，采用最大下坡坡度的線性函數對水流分配進行加權，結果更合理。王玉著等[9]從河網空間形態(tài)和算法效率方面對比了多流向、R＆N法和D8法，表明多流向具有優(yōu)勢。需要指出的是，地形因素影響流向，各種流向算法的應用要基于一些基本假設，應結合實際情況進行考慮。

3 模型與算法

數字流域提取的軟件工具主要有TOPAZ、RiverTools、ArcGIS。ArcGIS應用廣泛，功能強大，水文模塊是其很小一部分，以子流域作為流域特征值的單位統(tǒng)計，對整體區(qū)域參數的整理重復煩瑣[10]。葉愛中等[11]結合圖論與水文學思想，不再進行填洼處理，從流域出口直接向上搜索，通過圖的遍歷確定流向的自動提取河網與流域劃分的方法AEDNM。張鑫等[12]針對早期DEM數據缺乏、精度低、覆蓋范圍小等問題，建立KH-9數據的山地冰川提取流程。趙冰雪等[13]對山區(qū)水體提取出現(xiàn)的陰影問題和河流斷線問題提出解決思路，基于OLI和DEM數據提取方法簡單，面狀水體齊全，線狀水體連續(xù)。馬永明等[14]基于SWAT模型探討河道剖面和地形特征對水系提取精度的影響。

4 分辨率

通過提取河網，進行不同尺度水文模型計算，存儲地形參數。張宏才[15]認為尺度問題具有明顯的制約作用，不同分辨率DEM提取水系有很大差別?？垫面玫萚16]以1︰50 000和1︰250 000兩種尺度的數字化水系修正DEM進行對比，在精度較高、起伏較大的地區(qū)，后者用于修正效果較好。孫愛立等[17]選取我國50個子流域作為樣本，確定最佳集水面積，采用相關分析和多元性回歸分析，得出閾值與降水、坡度、地表覆蓋的定量關系。石磊等[18]進一步分析高程、坡度、地面粗糙度，來確定分辨率合理范圍。王英等[19]基于WEP-L模型計算坡度，利用分形結合半方差函數分析低分辨率向高分辨率的轉換規(guī)律。李思進等[20]基于點云數據建立不同分辨率的DEM數據集，研究對黃土侵蝕溝的地形形態(tài)影響因素。需要指出的是，不同分辨率對提取結果產生的影響較大，需要找到適宜該地區(qū)的分辨率。

5 地形區(qū)域

流域特征提取研究已經經歷了半個世紀，在各個區(qū)域里都有開展。徐新良等[21]將全國流域劃分為14大流域片分別提取。俞雷等[22]以重慶三峽庫區(qū)為研究對象，應用burn-in算法修正河道及兩側柵格單元，對發(fā)生偏差的平坦區(qū)域疊加到實際DEM中。不僅對流域，對城市集水區(qū)域也適用，如楊華容等[23]對巴中市的研究，左文君等[24]對城市河網水系提取的研究。對丘陵山區(qū)等高程地的研究較多，如陳勇等[25]研究丘陵山區(qū)土地平整土方量的計算；黃晶晶等[26]研究丘陵山區(qū)局地水系提取并討論分形分維特征；吳婕等[27]選擇高山、中山、低山3種地形樣區(qū)，凸顯山頂點和鞍部點一體化提取的優(yōu)點。研究的地形主要以山地為主，逐漸向平緩區(qū)改進，但使用DEM對平緩區(qū)的研究誤差較大。

6 結語

DEM數字流域水文特征提取廣泛運用于分布式水文模型，今后該研究依舊是熱點，已經形成了比較成熟的框架模型和算法，但存在一些實際問題需要考慮：1）流向的判斷中出現(xiàn)平原區(qū)時較難判斷，單流向對DEM原始數據要求較高，多流向則對算法要求較高；2）大量研究表明，DEM分辨率的不同會產生影響，具體怎么確定分辨率的合理值需要每個研究區(qū)視實際情況而定，精度的修正效果沒有參考值判定；3）現(xiàn)在使用的DEM提取數據主要針對具有高程的丘陵山地等地形，對于平原區(qū)提取結果存在較大誤差，人類活動也影響原始數據提取。需要進一步對適用于平原區(qū)、河道及河岸高差不明顯的區(qū)域開展研究，產生適合更多不同區(qū)域地貌的方法，為更廣泛的水文分析提供基礎。