宋瑞明，寧亞偉，李敬茹

（黃河勘測規(guī)劃設(shè)計研究院有限公司，450003，鄭州）

水文工作是支撐流域水利工程規(guī)劃的基礎(chǔ)，對工程規(guī)模的確定具有重要影響。 流域數(shù)字河網(wǎng)是構(gòu)建分布式水文模型的基礎(chǔ)數(shù)據(jù)， 流域面積、河長及比降等流域特征值是水文分析計算中推求設(shè)計洪水的制約性因素。 因此，對于無資料地區(qū)的水文工作，如何快速準確地獲得數(shù)字河網(wǎng)和流域特征值數(shù)據(jù)顯得尤為必要。 隨著地理信息系統(tǒng)（GIS）、數(shù)字化技術(shù)的不斷發(fā)展和推進，運用數(shù)字高程模型（Digital Elevation Model，DEM） 進行數(shù)字河網(wǎng)及水文信息提取已經(jīng)成為一種快速可行的方法。 目前，廣泛應(yīng)用的水文分析模型包括ArcGIS Hydrology、ArcSWAT 和Arc Hydro Tools等， 其中Hydrology 為ArcGIS 內(nèi)置模塊，ArcSWAT 和Arc Hydro Tools 均內(nèi)嵌于ArcGIS 軟件中。 DEM 是數(shù)字流域地形特征呈現(xiàn)的常用表達形式，描述了研究區(qū)域地面高程的分布情況，包含了豐富的水文、地貌信息，是提取流域水文特征進行河道匯流計算的重要基礎(chǔ)數(shù)據(jù)。DEM 數(shù)據(jù)類型包括柵格、矢量等高線、不規(guī)則三角網(wǎng)（TINs）等，其中柵格數(shù)據(jù)在數(shù)字河網(wǎng)提取及子流域劃分中較為常用。 根據(jù)前人學(xué)者研究的成果，在地勢起伏較大的山區(qū)， 河網(wǎng)提取的精度較高，而在人為干擾較大的地方或地勢平坦區(qū)，河網(wǎng)提取適用性較差，但在地勢平坦區(qū)如何提取精度較高的流域特征值研究較少。 本文運用GIS 平臺，通過工程實例對山區(qū)、丘陵平原區(qū)分別運用ArcGIS Hydrology 工具、Arc-SWAT 模型進行提取分析， 對提取效果進行評價， 探索地勢平坦地區(qū)河網(wǎng)提取的可行性方法， 為中小河流水利工程前期規(guī)劃中流域特征值的高精度自動提取提供參考。

一、研究區(qū)概況

桂林市位于廣西東北部和湘桂黔交界地帶，地形復(fù)雜多樣，中山、低山、丘陵、巖溶石山和河谷平原交錯分布。 中山主要分布在桂林市北部、西部和中部海洋山等地，低山主要分布在各大山脈四周，丘陵多分布于中低山與河流谷地之間，巖溶石山和河谷平原主要分布于桂江與湘江兩岸。轄區(qū)內(nèi)山脈縱橫， 河流水系發(fā)達，分布有桂江（漓江段）、古宜河、洛清江、湘江和資江五大河流。桂林市氣候溫和，雨量充沛，多年平均氣溫在17.8～18.8℃。 漓江流域多年平均年降雨量介于1 500～2 600 mm 之間，資江流域、湘江流域、古宜河流域、洛清江流域多年平均年降雨量分別為1 779.0 mm、1 704.5 mm、1 521 mm、2 010 mm。 本次研究對象為中小河流， 主要為洛清江一級支流西河、 相思江以及漓江一級支流小溶江、良豐河等。

西河為洛清江的一級支流，上游稱龍江， 發(fā)源于臨桂區(qū)黃沙鄉(xiāng)圍嶺村，流向自北西向東南，在永?？h城與東河匯合后流入洛清江，干流全長95 km，流域面積1 153 km2，多年平均年徑流量16.14 億m3。 主要支流包括灘頭河、壽城河、鳳凰河、舊縣河、山南河等。

相思江為洛清江的一級支流，發(fā)源于臨桂區(qū)廟頭鎮(zhèn)沙塘村的亂石山中，流向由北向南，在永福縣蘇橋鎮(zhèn)匯入洛清江，干流全長45 km，流域面積574.6 km2，多年平均年徑流量8.53億m3。 流域中上游屬峰林平原地貌，下游屬丘陵剝蝕地貌。 主要支流包括沙塘河、蔡塘河、沉橋河、清水河、睦洞河、會仙河、羅錦河等。

小溶江為漓江的一級支流，發(fā)源于資源縣兩水鄉(xiāng)塘洞村南，越城嶺西南麓，戴云山東南的大坳處，自北向南流，經(jīng)資源縣、興安縣、靈川縣境，匯入漓江。 源頭高程1 656 m，河口高程175 m，主河長49 km，流域集水面積269 km2。 小溶江流域穿行于越城嶺西南之山嶺中， 西北高東南低，屬狹長形流域，地勢陡峻。

良豐河是漓江流經(jīng)市區(qū)南部的一級支流，發(fā)源于臨桂、永福、陽朔三縣（區(qū)）交界的香草巖，先由南向北流， 而后至桂林市郊良豐后轉(zhuǎn)向東北，在柘木圩處匯入漓江。 良豐河干流全長58 km， 流域面積534 km2，多年平均年徑流量4.99 億m3。 在大江水庫壩址以上屬于山區(qū)，河流坡降較陡，以下全是平原和丘陵，河道坡降平緩。 主要支流包括奇峰溪、蘇河、周家河、大埠溪等。

二、數(shù)據(jù)來源與研究方法

1.數(shù)據(jù)來源

研究區(qū)ASTER GDEM 30 m 分辨率數(shù)字高程數(shù)據(jù)集來源于美國國家航空航天局（National Aeronautics and Space Administration，NASA），地理坐標系為WGS_1984，采用墨卡托投影，中央經(jīng)線111°。由于無實測河網(wǎng)資料，用于對比分析的河網(wǎng)數(shù)據(jù)來源于谷歌地球上描繪的數(shù)字化矢量河流數(shù)據(jù)。流域水文信息基礎(chǔ)數(shù)據(jù)來自于已批復(fù)的規(guī)劃報告、設(shè)計報告等相關(guān)資料。

2.研究方法

ArcGIS Hydrology 工具通過模擬地表水形成徑流的過程，并利用這一模擬過程實現(xiàn)河流、出水口以及流域的提取，為流域水文信息提取提供了便利條件。 Hydrology 使用的是D8 單流向算法，該算法假定雨水降落在某一柵格單元上，該柵格單元的水流只能流入與之相鄰的8 個柵格單元之間地形最低的單元中。 如果多個柵格單元的最大下降方向均相同，則擴大相鄰柵格范圍，直到確定最陡下降方向為止。 該工具河網(wǎng)提取流程主要為：DEM 數(shù)據(jù)裁剪獲得研究區(qū)范圍內(nèi)DEM，原始DEM 數(shù)據(jù)填洼處理生成無洼地的修正DEM，水流方向分析確定各柵格單元的水流流向，流量計算，設(shè)置匯流精度進行匯流累積量計算定義最小地表徑流， 運用河流鏈接和河網(wǎng)矢量化工具進行河網(wǎng)提取， 通過捕捉傾瀉點工具對出水口數(shù)據(jù)進行校正，運用分水嶺工具進行流域劃分等。

ArcSWAT 擴展模塊是SWAT 模型在ArcGIS 平臺上的圖形用戶界面，SWAT（Soil and Water Assessment Tool）模型是由美國農(nóng)業(yè)部（USDA）農(nóng)業(yè)研究中心和美國國家自然資源保護委員會聯(lián)合開發(fā)的基于GIS 的分布式流域水文模型， 主要由水文循環(huán)過程、土壤侵蝕以及污染負荷三個子模型組成。 其中水文循環(huán)過程分為兩個部分：一是坡面產(chǎn)流過程，指降雨降落到地表后，經(jīng)過植被截留、下滲、蒸發(fā)后，最終流入河流的過程，該過程主要基于水量平衡原理；二是河道匯流過程， 指各子流域內(nèi)河流匯集的水、泥沙和污染物等流至出水口的遷移過程，河道匯流計算采用馬斯京根法方程。 河網(wǎng)提取和子流域劃分是SWAT 模型構(gòu)建的前提， 其結(jié)果對徑流模擬成果有較大影響，是模型計算過程中必不可少的步驟。 該部分主要在Automatic Watershed Delineation模塊內(nèi)完成， 利用模塊內(nèi)Burn-in 算法引入實際矢量河網(wǎng)圖層，脅迫河網(wǎng)水系的生成， 在此基礎(chǔ)上結(jié)合D8 算法可得到更為真實的河網(wǎng)水系結(jié)果。Burn-in 法由Saunders 提出， 也被稱為Stream Burning 法、 凹陷化算法或河道燒錄法。 該方法原理為通過強制降低河網(wǎng)所在位置的DEM 柵格高程，使河網(wǎng)部分匯流能力提高，從而改善河網(wǎng)提取精度，并提高流域邊界劃分的準確性。ArcSWAT 模型提取水系及流域劃分主要步驟包括： 定義SWAT 工程、加載DEM、加載或手繪Mask 添加掩膜指定研究區(qū)范圍、加載實際河網(wǎng)水系以“引導(dǎo)”河網(wǎng)及子流域的生成、指定最小子流域面積生成河網(wǎng)、編輯流域出口節(jié)點進行子流域劃分、計算子流域參數(shù)等。

Hydrology 工 具 和ArcSWAT 模 型計算流程圖見圖1。

圖1 ArcGIS Hydrology 和ArcSWAT 模型流程圖

三、結(jié)果分析

1.ArcGIS Hydrology 河網(wǎng)提取

基于桂林西中部地區(qū)DEM 高程數(shù)據(jù)， 運用ArcGIS Hydrology 工具對研究區(qū)進行流向計算，然后設(shè)置閾值提取河網(wǎng)并進行流域劃分，流域面積及河長提取結(jié)果見表1。 由于研究區(qū)對象為中小河流，研究區(qū)范圍內(nèi)河流水系等基礎(chǔ)資料匱乏，為驗證河網(wǎng)提取精度，將谷歌地球上描繪的數(shù)字化矢量河網(wǎng)與提取河網(wǎng)結(jié)果進行疊加分析。

針對小溶江流域、 西河流域、相思江流域、良豐河流域河網(wǎng)提取及流域劃分結(jié)果， 根據(jù)地形條件分析可知， 小溶江流域地勢呈西北高東南低，河道坡降大；西河區(qū)域內(nèi)大霧山山脈主峰大霧山海拔高程1 291.6 m，為永?？h和融安縣界山，大霧山山脊線明顯，呈南北向展布，山脊為流域分水嶺； 兩流域均為大起伏山地區(qū)，流域河網(wǎng)提取走向基本與實際河網(wǎng)一致。 相思江流域地面高程大部分在146～152 m 之間，河床與岸坡高差1～3 m，流域地形平坦；良豐河流域兩岸地面高程大部在144～150 m 之間，流域地形平坦，河網(wǎng)偏差均較大。

從流域面積劃分結(jié)果來看，小溶江流域、西河流域提取流域面積與實際值偏差率分別為-3.87%、-1.98%，偏差較小，相思江和良豐河流域提取流域面積與實際值偏差率分別為18.24%和53.62%，偏差大。從干流河長提取結(jié)果分析，小溶江流域、西河流域提取值與實際值偏差率分別為-4.90%、-7.79%， 偏差在可接受范圍內(nèi)，相思江和良豐河流域提取值與實際值偏差率分別為-22.07%和59.33%，偏差大。 究其原因，主要是D8 算法采用8 鄰域柵格單元的遍歷順序，直到找到最陡下降方向， 小溶江流域、西河流域山脊線較為明顯，河道受兩岸山地限制， 柵格單元流向容易確定，河網(wǎng)提取結(jié)果與實際情況較為接近；而相思江和良豐河流域多為丘陵平原地貌，河道坡度平緩，流域分界線不分明， 柵格流向受遍歷順序的影響，平緩區(qū)所得柵格單元流向方向相同，易產(chǎn)生偽河道，且DEM 垂直分辨率相對于地形落差精度偏低，造成相思江和良豐河流域河網(wǎng)提取偏差大、河流流向及流域提取范圍均與實際不符的情況。 因此表明，在地勢起伏較大的山地，ArcGIS Hydrology 提取結(jié)果與實際情況較為一致，但在地勢平坦地區(qū)適用性較差。

2.ArcSWAT模型河網(wǎng)提取

基于桂林西中部地區(qū)DEM 高程數(shù)據(jù)，運用ArcSWAT 模型Burn-in 加載實際河網(wǎng)水系， 脅迫河網(wǎng)水系生成， 并在此基礎(chǔ)上進行子流域的劃分，得到小溶江、西河、相思江、良豐河流域河網(wǎng)水系和流域劃分成果，與谷歌地球上描繪的數(shù)字化矢量河網(wǎng)進行疊加。

從河網(wǎng)提取結(jié)果可以看出，運用ArcSWAT 模型Burn-in 方法修正河網(wǎng)后， 河網(wǎng)生成基本接近實際河網(wǎng)。為進一步驗證結(jié)果精度，采用流域劃分結(jié)果與實際流域特征值進行對比（見表1）?？梢钥闯觯∪芙?、西河、良豐河、相思江流域面積提取值與實際值偏差率分別為1.26%、-1.99%、1.29%、0.12%，河流長度提取值與實際值偏差率分別為6.27%、-6.21%、-0.34%、-7.56%。 可見，相比于流域面積來說， 河長提取值總體偏差率較大， 其原因一是由于匯流累積量計算時設(shè)置匯流精度閾值取值較大， 提取的河網(wǎng)分叉較少；二是DEM 高程數(shù)據(jù)分辨率偏低，模型提取精度受限，因而造成河長提取值偏小。 綜合來看，Burn-in 方法修正河網(wǎng)后，模型提取的不同坡度地區(qū)流域河網(wǎng)的流向、形態(tài)等特征更加趨于真實，流域特征值與實際值符合性較高，能夠較好地解決丘陵平原地區(qū)流域河網(wǎng)提取的問題，分析結(jié)果與其他學(xué)者研究成果相符合。

表1 流域特征值提取值與實際值對比分析

3.對比分析

通過運用ArcGIS Hydrology 與ArcSWAT 模型，利用DEM 數(shù)字高程數(shù)據(jù)，提取流域面積、河長等流域特征值與實際值進行對比分析，可知對于位于山區(qū)的小溶江、 西河流域，流域面積提取偏差率均在5%以內(nèi)，河長提取偏差率在10%以內(nèi)，兩模型提取效果均較好；而對于位于丘陵平原區(qū)的相思江、 良豐河流域，ArcSWAT模型流域面積及河長提取偏差率均在10%以內(nèi)，而利用ArcGIS Hydrology工具相思江水系河長提取偏差率達到59.33%， 提取精度明顯低于Arc-SWAT 模型。 因此，在地形平坦地區(qū)，運用ArcSWAT 模型Burn-in 方法對DEM 修正后，再進行河網(wǎng)提取和流域劃分能夠取得較好的效果，在項目前期工作中是可行的。

四、結(jié) 論

DEM 數(shù)據(jù)與ArcGIS 軟件應(yīng)用平臺為數(shù)字河網(wǎng)及流域水文特征的提取提供了有力支持， 本文通過運用ArcGIS Hydrology 及ArcSWAT 模型兩種流域特征值提取方法，基于30 m 分辨率ASTER GDEM 高程數(shù)據(jù)， 對桂林市西中部地區(qū)山區(qū)和丘陵平原區(qū)流域水系進行子流域劃分和河網(wǎng)提取， 并與實際矢量化河網(wǎng)進行對比，結(jié)果表明：

①基于DEM 數(shù)據(jù)， 運用ArcGIS Hydrology 和ArcSWAT 模型對山區(qū)河網(wǎng)進行提取與子流域劃分，水系提取結(jié)果均與實際河網(wǎng)走向、形態(tài)基本一致，流域面積、河長等流域特征提取值偏差均在10%以內(nèi)。 因此兩種方法均可解決研究區(qū)河流水文信息缺乏的問題，廣泛適用于各個山區(qū)流域的水文分析。

②運用ArcGIS Hydrology 和Arc-SWAT 模型對地勢平坦的丘陵平原地區(qū)河網(wǎng)進行提取與子流域劃分，由于地勢平坦地區(qū)水流無法沿某一方向流動形成徑流，ArcGIS Hydrology 河網(wǎng)提取效果與實際吻合度較差，易產(chǎn)生偽河道；流域面積、河長等流域特征提取值偏差較大，該方法對丘陵平原區(qū)適用性較差； 運用ArcSWAT 模型Burn-in 方法對DEM 修正后，河網(wǎng)提取和子流域劃分精度均較高，可作為提取丘陵平原地區(qū)流域面積等的有效方法，為中小流域水利工程規(guī)劃前期工作提供科學(xué)依據(jù)。 ■