劉家宏，秦大庸，李海紅，楊志勇

（中國水利水電科學(xué)研究院 水資源研究所，北京 100038）

1 研究背景

基于DEM的河網(wǎng)提取是分布式水文模型建立的基礎(chǔ)環(huán)節(jié)，提取河網(wǎng)的精度高低是影響水文模型模擬精度的一個(gè)重要因素。自Peuker［1］開始，河網(wǎng)提取的方法和模型已經(jīng)研究了30多年，主要形成了兩種生成河網(wǎng)水系的方法：一種是在一小窗口內(nèi)對(duì)高程進(jìn)行局部計(jì)算，然后在DEM上移動(dòng)該窗口來確定上凹面或V型面，這些面底部的網(wǎng)格單元被視為水道或河網(wǎng)的組成部份，這種方法的缺陷是生成的水系不連續(xù)，需要重新連接和修整才能生成較合理的河網(wǎng)［2］；另外一種是O’Callaghan和Mark［3］提出的坡面流模擬方法，它是依據(jù)水文學(xué)坡面流概念來判別水流路徑，根據(jù)每一柵格單元與相鄰單元之間的最陡坡度識(shí)別水流路徑，再選擇合適的水道集水面積確定河網(wǎng)，該方法是目前應(yīng)用最廣泛的水系提取方法。O’Callaghan和Mark提出的提取河網(wǎng)的方法，只有在DEM中每一個(gè)單元都確定了排水方向后才能自動(dòng)生成連續(xù)的河網(wǎng)。圍繞如何確定每個(gè)柵格單元的流向，F(xiàn)airfield和Leymarie［4］提出了D8方法、Rho8方法，此外還有Lea［5］方法、DEMON法［6］、D∞法［7］等，無論哪一種方法，DEM中的閉合洼地和平坦區(qū)域都會(huì)出現(xiàn)流路失真的問題，尤其是在強(qiáng)人類活動(dòng)下的平原地區(qū)，往往由于人工河道（例如運(yùn)河，排水溝渠等）的影響，改變了原有天然河道的流向，使得提取的水系與實(shí)際情況大相徑庭，從而使得水循環(huán)模擬結(jié)果精度降低。許多學(xué)者研究了解決這個(gè)問題的辦法，總結(jié)起來主要分為兩種思路：一種是對(duì)河網(wǎng)提取模型進(jìn)行改進(jìn)，例如Turcotte［8］在模型進(jìn)行流向和湖泊的判斷時(shí)，加入數(shù)字化河網(wǎng)和湖泊等輔助信息；另一種是修正輸入，即DEM。后者認(rèn)為之所以出現(xiàn)流路提取誤差，是因?yàn)镈EM高程信息不足造成的，與河網(wǎng)提取模型無關(guān)，即錯(cuò)誤的輸出是由信息不全的輸入造成的，因此后者通過對(duì)DEM進(jìn)行預(yù)處理，增加DEM的高程信息來解決，孔凡哲、芮孝芳［9］提出的閉合洼地和平坦區(qū)域的處理方法以及楊傳國、余鐘波［10］等提出的ZB算法即屬于這一類。本文綜合Turcotte方法和孔凡哲、芮孝芳方法的優(yōu)點(diǎn)，提出一種簡(jiǎn)單易行的平原河網(wǎng)流路強(qiáng)化方法，該方法利用已有的主要干流數(shù)字化的水系藍(lán)圖，將水系藍(lán)圖通過ArcGIS軟件轉(zhuǎn)換為與DEM網(wǎng)格尺寸相等的柵格圖像，將DEM圖像與柵格化的水系圖進(jìn)行疊加分析，“挖深”位于水系河網(wǎng)上的DEM單元的高程值，從而強(qiáng)化河流的主槽，得到正確的河網(wǎng)提取結(jié)果。

2 流路強(qiáng)化方法

基于DEM的河網(wǎng)提取方法需要根據(jù)DEM的高程值來判定每個(gè)單元的流向。在山丘地區(qū)，河谷和兩岸山丘的高差較大，在判定單元流向時(shí)比較容易；但在平原區(qū)，許多DEM并不能反映水下地形，DEM表示的常常是河流水面高程，由于平原河道的河流水面高程與兩岸高程高差并不太大，在多沙河流的下游平原，兩岸的地面高程甚至低于河流水面高程（例如黃河中游開封段就是“地上河”），低精度的DEM又不能反映堤防高程（堤防的寬度一般小于200m，柵格尺寸大于200m的DEM即難以反映河流兩岸的堤防），這就造成平原區(qū)各DEM單元的高程值相差很小，甚至完全一樣的情況，使得單元的流向很難判斷，因而出現(xiàn)河網(wǎng)提取中的“偽河道”和主流偏移等問題。

流路強(qiáng)化處理的基本思想是：人類活動(dòng)主要干擾的是河道中下游，即干流流經(jīng)的地區(qū)，在這些地區(qū)通過已有的較小比例尺的干流河道數(shù)字化地圖，確定平原地區(qū)的河道覆蓋的DEM網(wǎng)格單元集，將這個(gè)集合里的每一個(gè)單元的高程值減去一個(gè)值h，h取該地區(qū)河流的平均水深。通過這種處理，可以去掉河流水面的影響，基本反映河流的水下地形，使得DEM單元流向的判斷不再模糊，直接從兩岸指向河流主槽所在的低處。流路強(qiáng)化處理實(shí)質(zhì)上是在平原地區(qū)本來平坦的DEM上人工開挖出了已知河流電子地圖所反映的干流河槽，降低了平原區(qū)流向判斷的不確定性，使得地表漫流盡數(shù)歸槽。流路強(qiáng)化處理方法僅僅改變了與河流主槽重合的DEM單元的高程值，對(duì)流域整體地貌的影響小，不會(huì)影響其他地貌參數(shù)的正確提取，因而不會(huì)對(duì)流域的坡面產(chǎn)流造成影響，并且可以保證流域干流和主要支流河網(wǎng)與天然河網(wǎng)的一致性。

流路強(qiáng)化的處理方法簡(jiǎn)單易行，完全可以通過常用的ArcGIS軟件快速實(shí)現(xiàn)。處理過程分為如下4個(gè)步驟：

（1）為河網(wǎng)圖層增加表示河流平均深度的屬性字段Depth。在ArcGIS中打開河網(wǎng)電子地圖，編輯圖層屬性表，增加表示河流平均深度的屬性字段Depth，并給該字段賦值h，h取該地區(qū)河流的平均水深，h的值也可以根據(jù)不同的河道級(jí)別分別賦值。

（2）河網(wǎng)電子地圖柵格化。將河網(wǎng)地圖的坐標(biāo)轉(zhuǎn)換為與DEM相同的投影坐標(biāo)系，并將矢量格式的河網(wǎng)電子地圖轉(zhuǎn)化為與DEM格柵大小相同的柵格圖層，轉(zhuǎn)換的字段設(shè)定為Depth，即該區(qū)河道平均水深。轉(zhuǎn)換過程直接利用ArcGIS的投影轉(zhuǎn)換工具和文件格式轉(zhuǎn)換工具。

（3）河網(wǎng)柵格圖層與DEM柵格圖層匹配處理。匹配處理包括兩個(gè)方面：一方面需要將河網(wǎng)柵格圖層的圖幅范圍與DEM柵格圖層的圖幅范圍強(qiáng)制設(shè)定為一致，便于兩圖層的疊加分析，本操作通過ArcGIS的Export Data工具實(shí)現(xiàn)；另一方面需要將河網(wǎng)柵格圖層的NULL值全部轉(zhuǎn)化為0值，這是因?yàn)闁鸥窕暮泳W(wǎng)圖層除河流主流覆蓋的柵格值是h外，其他網(wǎng)格數(shù)值是NULL，與DEM作加減處理時(shí)，得到的值也是NULL，會(huì)造成地形數(shù)據(jù)的大量丟失，因此需要將河網(wǎng)柵格圖層的NULL值轉(zhuǎn)化為0，再作加減計(jì)算，本操作通過ArcGIS的Spatial Analyst中的Raster Calculator工具實(shí)現(xiàn)，計(jì)算公式為“con（isnull［河網(wǎng)柵格圖層］，0，［河網(wǎng)柵格圖層］）”。

（4）DEM柵格圖層流路強(qiáng)化處理。將經(jīng)過匹配處理的DEM柵格圖層與河網(wǎng)柵格圖層相減，得到最終的DEM柵格數(shù)據(jù)，用于河網(wǎng)水系的提取。本操作通過ArcGIS的Spatial Analyst中的Raster Calculator工具實(shí)現(xiàn)，計(jì)算公式為“DEM柵格圖層—河網(wǎng)柵格圖層”。

3 海河流域的應(yīng)用

海河流域位于東經(jīng)112°～120°，北緯35°～43°之間，流域面積31.78萬km2，其中平原區(qū)面積14.96萬km2，占流域總面積的47%。流域的平原區(qū)人口密度大，人類活動(dòng)強(qiáng)烈，對(duì)河流的影響大，還有著名的京杭運(yùn)河縱貫流域東部平原，河網(wǎng)水系極為復(fù)雜，直接通過DEM提取河網(wǎng)水系，通常存在較大的誤差。本文通過流路強(qiáng)化方法來提取河網(wǎng)水系，并給出直接提取和流路強(qiáng)化法取的對(duì)比結(jié)果。

3.1 數(shù)據(jù)準(zhǔn)備 本文采用美國宇航局（NASA）2000年通過航天飛機(jī)雷達(dá)拓?fù)錅y(cè)繪（SRTM）得到的90m×90m精度的DEM數(shù)據(jù)（圖2）。數(shù)字化河網(wǎng)藍(lán)圖采用全國標(biāo)準(zhǔn)的5級(jí)河流電子地圖，在海河流域主要涉及3、4、5級(jí)河流，由于第5級(jí)河流枝杈太密集，本著河流強(qiáng)化主要是干流和主要支流的原則，不考慮第5級(jí)河流，海河流域的3、4級(jí)主要河流見圖3。標(biāo)準(zhǔn)的河網(wǎng)圖中包含了京杭運(yùn)河，由于京杭運(yùn)河分割了海河原有的天然流域，不便于流域劃分，因此在河網(wǎng)處理時(shí)去掉部分京杭運(yùn)河連接線，這樣在流路強(qiáng)化處理時(shí)就不再強(qiáng)化切割天然流域的人工河道，僅強(qiáng)化了原有的干流和主要支流河道。

3.2 河網(wǎng)提取結(jié)果對(duì)比 本文利用TOPAZ［11］模型，分別對(duì)流路強(qiáng)化處理前和處理后的DEM進(jìn)行了河網(wǎng)水系提取操作，提取的主要河流結(jié)果見圖4、圖5。圖中淺灰色的粗線是已知的主要水系數(shù)字化藍(lán)圖，黑色細(xì)線是應(yīng)用TOPAZ模型提取的主要河網(wǎng)。從圖4中可以看出，沒有經(jīng)過流路強(qiáng)化處理的DEM提取的河網(wǎng)與實(shí)際的藍(lán)線河網(wǎng)差別較大，尤其是在人類活動(dòng)劇烈的平原地區(qū)，個(gè)別主要河流的主流發(fā)生了明顯的偏移，例如流域南部的徒駭河和馬頰河，還有潮白河在下游直接襲奪了薊運(yùn)河的水道；一些河流在下游入海口的位置發(fā)生了明顯的偏移，例如灤河下游，入?？谙蚰掀屏思s30km。經(jīng)過流路強(qiáng)化處理的DEM提取的河網(wǎng)結(jié)果見圖5，可以看到提取結(jié)果（黑色細(xì)線）與實(shí)際河網(wǎng)（灰色粗線）基本一致，提取河網(wǎng)由于臨界給水面積（Critical Source Area，簡(jiǎn)稱CSA）是統(tǒng)一設(shè)定的，所以顯示的流路長(zhǎng)短與實(shí)際的數(shù)字化河網(wǎng)有些差異，這主要是顯示，并不影響水系內(nèi)部的匯流結(jié)構(gòu)與天然河網(wǎng)的一致性。

實(shí)際上，通過這種方法不僅可以保證干流和主要支流流路的正確性，而且可以得到比現(xiàn)有數(shù)字化河網(wǎng)更加細(xì)密的河網(wǎng)拓?fù)浣Y(jié)構(gòu)，方法就是在河網(wǎng)提取時(shí)減小CSA的設(shè)定值。基于細(xì)致的河網(wǎng)拓?fù)浣Y(jié)構(gòu)，可以進(jìn)行分布式水文模型水文響應(yīng)單元（Hydrological Response Unit，簡(jiǎn)稱HRU）的劃分，圖6、圖7分別是海河流域加密的河網(wǎng)和相應(yīng)的HRU的劃分結(jié)果，以此為基礎(chǔ)即可進(jìn)行相關(guān)數(shù)據(jù)的空間展布，為流域水循環(huán)的精細(xì)模擬奠定基礎(chǔ)。

4 討論和結(jié)論

平原地區(qū)地勢(shì)平坦，河道的溝谷形態(tài)不明顯，DEM高程值極為接近，再加上平原區(qū)人類活動(dòng)劇烈，各種引水、排水溝渠密布，常規(guī)的河網(wǎng)提取模型往往難以得到正確的河網(wǎng)提取結(jié)果，給分布式水文模型在平原區(qū)的應(yīng)用造成一定困難。而平原地區(qū)往往又是水循環(huán)、水資源研究的重點(diǎn)，因此研究強(qiáng)人類活動(dòng)平原地區(qū)的河網(wǎng)提取方法十分必要。本文提出的流路強(qiáng)化方法汲取了Turcotte方法和孔凡哲、芮孝芳方法的優(yōu)點(diǎn)，用已知的河道干流的數(shù)字化信息，在DEM網(wǎng)格上沿正確的河道線進(jìn)行人工“開挖”，強(qiáng)化流路線，然后用處理后的DEM作為輸入，提取河網(wǎng)水系，能保證提取的河網(wǎng)與已知的數(shù)字化河網(wǎng)一致，從而保證分布式水文模型匯流結(jié)構(gòu)的正確性。流路強(qiáng)化方法操作簡(jiǎn)單，完全可以通過ArcGIS軟件實(shí)現(xiàn)，不需要修改河網(wǎng)提取模型，是一種方便、快捷的處理方法。通過流路強(qiáng)化法，改變的只是河道所在的DEM網(wǎng)格的高程值，在很大程度上是修正了原有DEM不能反映水下地形的缺陷（在河道網(wǎng)格單元上減去h，實(shí)際上是還原河底高程，因?yàn)閔取值為該河段的平均水深），因此不會(huì)對(duì)模型的產(chǎn)匯流參數(shù)造成壞的影響。在還原河底高程時(shí)，由于平均水深h是通過河段屬性表添加設(shè)置的，因此可以根據(jù)不同的河流級(jí)別分別設(shè)置，對(duì)于干流河道，可以設(shè)置得大一些；對(duì)于主要支流，可以適當(dāng)設(shè)置得小一些，避免支流挖得太深，造成匯流坡度變異，影響匯流計(jì)算精度。本文已知的干流河道數(shù)字化信息是指以往通過地理調(diào)查測(cè)繪出來的河流信息，往往只有干流和較大支流的矢量圖，沒有分布式水文模型要求的細(xì)密的河網(wǎng)拓?fù)湫畔ⅲ黄鋵傩孕畔⒂邢?，缺乏河道坡度、匯流面積等信息；這些信息的獲取一般通過基于DEM的河網(wǎng)提取模型提取，為此本研究成果可望對(duì)分布式水文模型在高強(qiáng)度人類活動(dòng)平原區(qū)的推廣應(yīng)用起到積極的推動(dòng)作用。

［1］Andrea T.Automated recognition of valley lines and drainage networks form grid digital elevation models：A review and a new method［J］.Journal of Hydrology，1992，139：263-293.

［2］晉華，楊金海，任煥蓮，等.基于DEM的濁漳河南源水系研究［J］.太原理工大學(xué)學(xué)報(bào)，2006，37（6）：646-648.

［3］O’Callaghan J F，Mark D M.The extraction of drainage networks from digital elevation data［J］.Computer Vision，Graphics，and Image Processing，1984，28：323-344.

［4］Fairfield J，Leymarie P.Drainage networks from grid digital elevation models［J］.Water Resources Research，1991，27（5）：709-717.

［5］Lea N L.An aspect driven kinematic routing algorithm in overland flow：Hydraulics and erosion mechanics［M］.New York：Chapman&Hall，1992.

［6］Costa-Cabral M C，Burges S J.Digital elevation model networks（DEMON）：A model of flow over hillslope for computation of contributing and dispersal areas［J］.Water Resources Research，1994，30（6）：1681-1692.

［7］Tarboton D G.A new method for the determination of flow direction and upslope areas in grid digital elevation models［J］.Water resources research，1997，33（2）：309-319.

［8］Turcotte R，F(xiàn)ortin J P，Rousseau A N，et al.Determination of the drainage structure of a watershed using a digital elevation model and a digital river and lake network［J］.J.Hydrol.，2001，（240）：225-242.

［9］孔凡哲，芮孝芳.處理DEM中閉合洼地和平坦區(qū)域的一種新方法［J］.水科學(xué)進(jìn)展，2003，14（3）：290-294.

［10］楊傳國，余鐘波，林朝輝，等.大尺度分布式水文模型數(shù)字流域提取方法研究［J］.地理科學(xué)進(jìn)展，2007，26（1）：68-76.

［11］Garbrecht J，Martz L W.TOPAZ overview［R］.USA：USDA ARS，Grazinglands Research Laboratory，1999.