基于DEM的丹江口水源區治理區小流域劃分研究

郭文慧 ，于 泳 ，李 璐 ，楊 偉 ，袁修猛

(1.湖北省水利水電科學研究院 湖北省水土保持工程技術研究中心，湖北 武漢 430070；2.湖北省水利廳，湖北 武漢 430070)

隨著計算機技術和地理信息系統的發展，數字高程模型DEM(Digital Elevation Model)中蘊涵著大量的地形、地貌信息，可以有效地反映區域的基本地形空間分布規律與地貌特征，通過分析提取河網、流域界線甚至提取流域內的地形屬性都是切實可行的。隨著數字水文的興起和分布式水文模型研究的發展，利用DEM獲取水系和子流域特征，進而進行水文模擬的方法越來越廣泛地為水文學者所采用。

本研究以1∶5萬DEM數字高程模型為基礎，在ArcMap軟件的水文模塊(Hydrology)環境下，探索了流域水系自動提取及流域劃分的方法，并依據《小流域劃分及編碼規范》(SL 653—2013)的相關要求對子流域進行了人工修正、歸并，最后完成了符合規范要求的適宜面積的小流域劃分。

1 研究區概況

以湖北省境內國家水土保持重點治理區劃中的丹江口水源區治理區(以下簡稱“治理區”)為研究對象，包括十堰市的張灣區、茅箭區、鄖西縣、鄖縣、丹江口市、竹山縣、竹溪縣、房縣共8個縣(市、區)。治理區位于秦嶺東西向構造體系的南部邊緣，地形的主要特點是高差大、坡度陡、切割深，總的地勢是西北高、東南低、北陡南緩，漢江沿線地貌峽谷和盆地相間。在我國氣候區劃上位于南北氣候過渡地帶的秦巴山區，北亞熱帶季風氣候，多年平均氣溫為15.9 ℃，年均降水量約900 mm，年均水面蒸發量860 mm左右。土壤以山地黃棕壤和黃褐土為主，山地棕壤土次之，并有紫色土發育。

丹江口水庫是我國水資源調整的重大工程——南水北調中線工程水源地，1973年建成，設計蓄水位157 m，水庫面積745 km2，庫容174.5億m3(其中漢江庫區占53.9%、丹江庫區占46.1%)，漢江、丹江回水長度分別為174和83 km。南水北調中線工程實施后，丹江口大壩加高至176.6 m，設計蓄水位170 m，總庫容290.5億m3，水庫面積1 050 km2，回水長度漢江194 km、丹江93 km。漢江庫區的主要入庫支流有金錢河、將軍河、天河、曲遠河、堵河、神定河、泗河、芝河、浪河等,丹江庫區的主要入庫支流有丹江、淇河、老灌河等。

研究區DEM情況見圖1。

圖1 研究區DEM情況

2 數據來源

(1)DEM數據。研究區1∶5萬DEM數據(柵格大小為25 m×25 m，見圖1)，投影坐標系為WGS84坐標系，UTM投影，6度分帶的49號帶。

(2)水系數據。 研究區5 km以上的河流水系圖。

(3)其他數據。研究區1∶5萬各縣級行政區劃矢量數據。

3 研究方法

由DEM提取流域的數字特征，包括DEM預處理、流向分析、 匯流累積分析及河網提取、流域邊界生成等過程。

(1)DEM預處理。從 DEM中自動提取流域的前提是DEM中不能存在洼地和尖峰，否則在計算水流方向時會出現逆流的現象，給水流線和流域界線的確定帶來困難，為此首先要將所有地形改造成由斜坡構成[1]。

(2)流向分析。在填洼后的DEM中，每個柵格單元都有一個可以定義的水流方向值。水流方向的確定有單流向法和多流向法兩種。單流向法包括D8算法、Rh08算法、Lea算法和D∞算法等，多流向法則包括MFD算法、DAEMOM算法和Ding算法等[2]。本研究采用目前較常用的D8算法,該算法是假定每個單元格的水流方向可能流向8個相鄰的單元格，水流沿最陡坡度的方向流動，計算出中心單元格與周圍8個單元格的坡度即可確定中心單元格的水流方向[3]。

(3)匯流累積分析。從每個柵格單元出發依次掃描流向矩陣，沿水流方向追蹤到DEM邊界，當整個水流方向矩陣掃描完畢，就可以得到流域匯流能力的柵格分布圖，匯流柵格上每個單元的值代表上游匯流區內流入該單元的上游柵格單元格的總數，值較大者可視為河谷，值等于零則是較高的地方，可能為流域的分水嶺。根據匯流累積柵格分布圖可以方便地提取流域的各種特征參數[4]。

(4)河網提取。首先設定最小支流上游集水區的面積閾值，以上游集水區面積大于閾值面積的網格點作為該支流的起始點，流域內集水面積超過該閾值的網格點定義為水道[5]。將匯流累積柵格上所有大于或等于最小集水面積閾值的柵格提取出來，即得到了柵格河網，然后將其轉化為矢量河網。

(5)流域邊界生成。流域又稱集水區域，是指流經其中的水流和其他物質從一個公共的出水口排出，從而形成一個集中的排水區域。集水流域的生成需要先確定該集水區域的最低點，然后結合水流方向數據，分析搜索出該出水點上游所有流過該出水點的柵格，直到所有的該集水區的柵格都確定了位置，也就是搜索到流域的邊界或分水嶺的位置。

(6)小流域歸并。本研究以《小流域劃分及編碼規范》(SL 653—2013)為依據進行小流域歸并工作。小流域面積原則上控制在30～50 km2，特殊情況不宜小于3 km2或大于100 km2。

4 流域特征提取

4.1 無洼地DEM的生成

4.1.1 洼深計算

首先借助Hydrology工具集中的Flow Direction工具獲取水流方向數據。ArcGIS中的水流方向是利用D8算法，即最大距離權落差(最大坡降法)來確定的。

運用Sink工具計算獲得洼地區域，并用Watershed工具生成洼地的貢獻區域，分別用Zonal Statistic和Zonal Fill工具計算得到洼地所形成的貢獻區域最低高程和洼地出水口高程。洼地深度即為出水口高程與貢獻區域最低高程之差。通過計算，研究區DEM最大洼地深度為88 m。

4.1.2 填洼閾值確定

分別將30、50、70、90 m作為填充深度值，在Hydrology工具中集中利用Fill工具將其處理，并對填洼后的DEM進行匯流累積分析。對比分析匯流累積圖水系特征信息與已配準的河流水系圖，確定最佳填洼閾值。

比較發現，30 m填洼結果坡面洼地有所減少，50 m填洼結果較為理想，70 m填洼結果與50 m填洼結果基本相同，90 m填洼結果所生成的匯流累積圖中的水系特征與河流矢量圖差異較大。為避免過度填洼，選擇50 m為填洼閾值。

4.2 河網的自動提取

為使獲取的河網更加準確，即水系的分級和水系的長度更加符合實際，在閾值設定時進行了一系列試驗，分別選取5 km2(8 000個柵格)、10 km2(16 000個柵格)、20 km2(32 000個柵格)作為最小匯流面積閾值，河網提取結果依次見圖2、3、4。

將不同閾值條件下的河網提取結果及所提取河網與水系圖比較發現：5 km2匯流面積條件下提取的河網偏密，且河流長度與實際水系一致或稍短；10 km2匯流面積條件下提取的河網稍密，且河流長度較實際水系稍短；20 km2匯流面積下提取的河網稀疏，且河流長度較實際水系短。

圖2 5 km2匯流面積的河網提取結果

圖3 10 km2匯流面積的河網提取結果

圖4 20 km2匯流面積的河網提取結果

4.3 流域邊界生成

對比分析表明，不同閾值條件下所生成的小流域數量與河道數一致，且一一對應。經統計，最小匯流閾值為5 km2時生成子流域2 923個，最大子流域面積為63.83 km2，50 km2以上子流域9個；最小匯流閾值為10 km2時生成子流域1 261個，最大子流域面積為131.04 km2，50 km2以上子流域49個；最小匯流閾值為20 km2時生成子流域615個，最大子流域面積為194.86 km2，50 km2以上子流域161個。

5 小流域邊界處理與劃分結果

5.1 子流域邊界處理與歸并

5.1.1 子流域邊界整合

根據《小流域劃分及編碼規范》(SL 653—2013)中的相關規定：小流域面積原則上控制在30～50 km2，特殊情況下不宜小于3 km2或大于100 km2。為提高數據處理效率并保證小流域邊界的準確性，本研究選擇子流域面積不大于50 km2，且子流域數量相對較少的劃分結果用于小流域的編碼及命名。為此，選擇最小匯流閾值為10 km2時所提取生成的子流域邊界數據，對于大于50 km2的子流域，用最小匯流閾值為5 km2的提取結果進行數據替換，并將小于0.1 km2的子流域與周邊子流域自動合并，最終得到替換合并更新后的數據(子流域1 742個，大于50 km2的子流域2個)用于小流域的修正、分割及歸并。

5.1.2 子流域邊界修正

由于基于DEM的河網提取是采用最大坡降的方法，因此在平地區域(例如谷底等)水流方向是隨機的，就很容易生成平行狀的河流等錯誤形狀，出現大量復雜鑲嵌的流域格局(如圖5)， 此時需要利用Editor工具欄進行手工編輯修正。

圖5 復雜鑲嵌的流域邊界示例

5.1.3 子流域邊界分割

首先在屬性表中添加“亞單元標識碼”，選出跨縣的小流域，賦值為1，其他未跨縣小流域賦值為0；然后用湖北省分縣行政邊界矢量數據對跨縣子流域進行分割。

5.1.4 子流域邊界歸并

依照《小流域劃分及編碼規范》(SL 653—2013)中的原則對子流域邊界進行歸并，受小流域面積的限制，歸并具有一定的主觀性。現將這次歸并研究過程中所遇到的特殊問題處理如下：

(1)被縣界分割形成面積較小的子流域。對此類子流域研究其與周邊子流域的關系，符合歸并原則的與周邊子流域歸并得到面積適合的小流域；對不符合歸并原則，且面積小于3 km2的子流域，為方便流域管理，考慮與周邊小流域合并；面積大于3 km2的子流域可單獨作為小流域。

(2)丹江口水庫及周邊區域。丹江口水庫周邊區域及庫底地形復雜，基于DEM自動提取的子流域邊界也較為復雜，較多子流域邊界跨越水庫水面，與對岸一起形成區間型小流域，不利于小流域的管理。基于此，考慮將170 m淹沒線以內區域(包括丹江口水庫水面及庫周170 m淹沒線以內區域)劃分為丹江口水庫型流域，將周邊的子流域邊界依據水系匯流關系及行政邊界進行歸并。

(3)盲區。在省界邊緣，會出現少數面積較小的流入其他省份的水系，他們與本省內周邊相鄰子流域的匯流關系較遠。這部分子流域多小于5 km2，分析認為造成這種結果一般有兩種原因：①實際地形地貌復雜；②省級邊界與DEM影像之間在地理位置上存在誤差。由于提取出來的水系與湖北省河流分布圖一致性較好，因此可排除第二種原因。考慮到便于流域管理，依據周邊水系的匯流關系進行歸并，形成類似于坡面小流域的類型(如圖6)。

圖6 特殊情況示例

5.2 小流域劃分結果

本研究在丹江口水源區治理區約2.30萬km2的區域內，提取小流域624個(丹江口水庫型流域除外，該流域面積為389.04 km2)，最小的小流域面積為3.41 km2，最大的小流域面積為59.42 km2。其中：<30 km2的小流域有122個，面積約為0.25萬km2，占該區總面積的10.87%；30～50 km2的小流域有455個，面積約為1.80萬km2，占該區總面積的78.26%；≥50 km2的小流域有47個，面積約為0.25萬km2，占該區總面積的10.87%。研究區小流域劃分結果見圖7。

圖7 研究區小流域劃分結果

6 結論與討論

利用1∶5萬比例尺DEM數據提取的河網與實際河網的比較結果一致性較好；在提取水系圖時，由于水系為連續的匯流最低點相連接而形成的線狀圖，因此水系只能反映河流的走向，而不能反映河流的寬度。自動提取的子流域的邊界與分水嶺基本吻合，流域特征明顯，但由于河底或庫底地形的復雜性，出現了大量復雜鑲嵌的流域格局，存在較多子流域需要人工修正，而人工修正又存在主觀性。

小流域是水土流失治理工作的基本單元。本研究對典型區域的小流域劃分技術進行了研究，為湖北省開展類似區域小流域劃分工作提供了方法。

[參考文獻]

[1] 陳加兵，勵惠國，鄭達賢，等．基于DEM的福建省小流域劃分研究[J]．地球信息科學，2007，9(2):74-77.

[2] GONG J Y，XIE J B.Extraction of drainage networks from large terrain datasets using high through put computing[J]. Comp-uters & Geosciences，2009，35:337-346.

[3] 曹玲玲，張秋文．基于SRTM的數字河網提取及其應用[J]．人民長江，2007，38(8):150-152.

[4] 王俊啓．基于1∶1萬地形圖的數字流域劃分研究——以浙江省安吉縣為例[D].杭州：浙江大學，2013：1-55.

[5] 秦福來，王曉燕，王麗華，等．基于DEM的流域特征提取及應用[J].河北師范大學學報(自然科學版)，2006,30(1):109-112.