遼寧省小型河流空間分布及屬性信息研究

董婷婷（遼寧省水利水電科學研究院，遼寧 沈陽 110003）

遼寧省小型河流空間分布及屬性信息研究

董婷婷
（遼寧省水利水電科學研究院，遼寧 沈陽 110003）

河湖是地球表面水循環、碳循環、營養物循環和泥沙循環的載體，也是哺育人類歷史文明的偉大搖籃，在區域生態安全、經濟社會可持續發展、文化傳承等方面具有重要意義。本文選取遼寧省內10～50km2河流作為研究對象，充分利用1∶5萬DEM、水系軸線、遙感影像等資料，采取計算機自動提取與人工目視判讀相結合的方法，實現了10～50km2河流提取，并獲取了河流屬性信息。與現存資料進行對比分析，本文提取的結果與實際情況較為一致，為我省典型河湖健康調查評價提供了數據保障。

河流；信息；提取；研究

河湖是地球表面水循環、碳循環、營養物循環和泥沙循環的載體，也是哺育人類歷史文明的偉大搖籃，在區域生態安全、經濟社會可持續發展、文化傳承等方面具有重要意義。在大規模經濟開發和全球氣候變化雙重因素的交織作用下，河湖生態健康問題十分突出，不同程度地出現了水質污染、形態結構破壞、水文條件惡化、生境退化、生物多樣性受損以及重要或敏感水生生物消亡等問題。因此，開展河湖健康評估工作，具有十分重要的戰略意義。

遼寧省水系發達，河流眾多。第一次水利普查成果顯示，全省流域面積50km2以上河流共有845條，河流總長度28459km。長期以來，受自然環境和人類活動的雙重影響，我省河流健康問題逐漸顯現，主要表現為:河道來水量減少、淤積持續增加、水面不斷萎縮、行洪能力下降、嚴重的甚至出現了中小河流消失的現象；河流污染問題突出，城鄉居民飲水安全受到威脅；河流所處流域生態環境惡化，植被破壞嚴重、水土流失加劇、旱澇災害頻發。針對遼寧省河流特點，遼寧省省政府于2013年10月31號聽取了全省水利工作專題匯報后，指示“水利廳要搞一個調查研究，建國以來有多少河流消失了，消失的原因是什么，采取什么措施能恢復多少河流”。對此，遼寧省水利廳決定對全省10km2以上典型河流開展河湖生態健康狀況調查評價課題，其中10～50km2河流提取是其子題之一。

我省已于第一次全國水利普查期間完成了50km2及以上河流的普查工作，本次河流提取工作在其成果之上展開，充分利用1∶5萬DEM、水系軸線、遙感影像等資料，采取計算機自動提取與人工目視判讀相結合的方法，開展10～50km2河流提取工作。

1 研究區概況

遼寧省位于祖國東北地區南部，地理坐標介于東經118°50′～125°47′，北緯38°43′～43°26′之間。全省水系發達，河流眾多。共有三大流域七大水系。三大流域分別為遼河流域、黑龍江流域和海河流域。七大水系分別為遼河流域的遼河水系、遼東灣西部沿渤海諸河水系、遼東灣東部沿渤海諸河水系、遼東沿黃海諸河水系和鴨綠江水系，黑龍江流域的松花江水系，海河流域的灤河及冀東沿海諸河水系。全省流域面積50km2以上的河流有845條。其中，山地河流787條，平原水網河流57條，混合河流1條。山地及混合河流按流域面積劃分為三類，分別是大型河流、中型河流和小型河流。其中，大型河流（流域面積5000km2以上）有16條，分別是遼河、老哈河、東遼河、清河、柳河、繞陽河、渾河、太子河、大凌河、小凌河、大洋河、鴨綠江、渾江、愛河、輝發河和青龍河；中型河流（流域面積為1000km2～5000km2）有32條；小型河流（流域面積為50km2～1000km2）有740條。在第一次全國水利普查期間，我省已經完成了50km2以上河流的普查工作，掌握了河流的數量、分布及相關屬性信息。此次河流提取工作主要針對10～50km2河流展開。

2 研究內容

河流指陸地表面宣泄水流的通道，是溪、川、江、河等的總稱。此次10～50km2河流提取主要考慮在天然和人工控制條件下常年（經常發生）自然的匯流（排水）功能，并有集水區域的河流。不管河流是否常年有水（如北方地區），不管河流的集水區域邊界能否清晰準確劃定（如平原水網區），也不管河流是天然河道還是人工整治河道。

本研究的主要內容包括兩部分，一是提取河流的空間分布信息，二是填補河流的屬性信息。其中河流的空間分布信息是指10～50km2河流分布圖，河流的屬性信息主要包括:河流名稱、流域面積、長度、級別、岸別、上級河流、河源地址、河口地址、流經、平均比降、河源點高程、河口點高程和河流類型13個字段。

3 研究方法

根據研究內容和我省河流特點，制定了技術路線。對于山地及混合河流采用計算機自動提取方法，對于平原水網河流采用人工目視判讀方法，最后，將兩者提取的結果進行匯總得到全省10～50km2河流，進而利用計算機編程自動獲取河流屬性信息。

3.1 數據源

本次河流提取所用到的數據源主要有以下幾種:

（1）1∶5萬DEM:主要用于山地及混合河流提取，河源點高程和河口點高程信息獲取等；

（2）0.5m分辨率的高清影像和水系軸線:主要用于平原水網河流提取；

（3）50km2以上河流信息，主要用于獲取級別、岸別、上級河流等信息；

（4）50km2河流流域界限分布圖，主要用于判定山地及混合河流和平原水網河流的分布范圍；

（5）灌區分布圖，主要用于輔助提取平原水網河流；

（6）省、市、縣、鄉和村5級行政區劃圖，主要用于獲取河源地址、河口地址、流經等信息；

3.2 山地及混合河流提取方法

山地河流指流域內地形起伏大、流域邊界能清晰界定的河流。混合河流指流域內部分區域流域邊界能清晰界定、部分區域流域邊界無法清晰界定的河流。山地及混合河流提取在50km2以上山地及混合河流所在流域范圍內展開，采用坡面徑流模擬法，在無洼地的DEM上利用最大坡降方法獲取每個柵格單元的流向，按照水流由高向低流動的自然規律，計算每個柵格在水流方向上的匯流累積量，通過設定集水面積閾值提取河網信息，進而對10～50km2的山地及混合河流進行篩選。

山地及混合河流的提取利用Arcgis軟件提供的Hydro1ogy模塊來實現，主要的操作步驟有:無洼地填充、水流流向確定、匯流累積量計算、河網提取和10～50km2山地及混合河流篩選。（1）無洼地填充:在DEM制作過程中，不可避免會產生一些錯誤的洼地和尖峰，對于洼地點只有水流流入沒有流出，對于尖峰點只有水流流出沒有流入，這些錯誤點會影響水文特征提取。因此，本文首先利用fi11功能對1∶5萬DEM數據進行無洼地處理。（2）水流流向確定:水流方向決定了河流走向，是提取流域水文特征的關鍵步驟。本文采用單流向法中的D 8方法，利用最大坡降計算每個柵格單元與其領域內8個柵格的距離權落差，將距離權落差的最大值作為該柵格的流向，此步驟利用f1owdirection功能對無洼地填充的結果進行流向計算。（3）匯流累積量計算:匯流累積量表征了河流的匯集能力，其值越大匯流能力越強，成為河流的可能性越大，是提取河網信息的重要因素。匯流累積量是在水流方向的基礎上實現的，假定每個柵格單元攜帶一個單位的水量，遵循自然水流由上游到下游流動，由支流到干流匯集的規律，根據水流方向數據計算每個柵格單元在流向上的水量累計值，即匯流累積量。本文利用f1owaccumu1ation功能對河流流向結果進行匯流累積量計算。（4）河網提取:當匯流累積量達到一定閾值時，就會產生地表水流，區域內所有滿足條件的水流所形成的河流網絡就是河網。河網的提取以匯流累積量結果為基礎，以集水區域面積大小作為判斷閾值，凡是匯流累積量大于閾值的柵格均被定義為河網內的單元，這些單元彼此連接就形成了河網。本研究集水面積閾值設定為0.2km2，利用 mapa1gebra功能提取河網信息。（5）10～50km2山地及混合河流篩選:流域面積大小是判斷河流類型的關鍵因素。本研究利用watershed功能，根據水流方向信息，將從同一出水口流水的所有柵格單元劃分為同一子流域，通過對流域面積進行篩選，提取出10～50km2山地及混合河流。

3.3 平原水網河流提取方法

平原水網河流指流域內地形起伏小、流域邊界無法清晰界定的河流。平原水網河流的提取在50km2以上平原水網河流所在流域范圍內展開，利用0.5m分辨率的高清影像和水系軸線作為主要數據源，采用目視解譯的方法進行人工判定。

平原水網河流的提取利用arcgis軟件實現，主要的操作步驟有:（1）加載數據源:數據源是進行平原水網河流提取的主要信息來源，這里利用的信息主要有0.5m分辨率的高清影像、水系軸線、灌區分布圖、50km2河流流域界限分布圖等。其中，50km2河流流域界限分布圖主要用于界定平原河網的分布范圍；0.5m分辨率的高清影像、水系軸線和灌區分布圖主要用于判斷河流。（2）矢量化:矢量化是獲取平原水網河流的重要環節。這里首先判斷平原水網河流，對于已經確定的平原水網河流，利用E di t o r工具矢量化河流，并進行保存。（3）10～50km2平原水網河流篩選:在分水線能清晰準確劃分的較大區域邊界內，根據重要性選定常年自然匯流的骨干河流，再由骨干河流逐級向周邊延伸選定河流，直至河流所在區域范圍的面積小于10km2；但該區域邊界內的河流總數不得超過區域邊界總面積與10km2的比值，并適當考慮區域內河流分布的均勻性。工作中發現不少河流、渠系長度很短，這里針對我省河流特點將平原水網的長度標準界定為不低于6km。

3.4 河流屬性信息獲取

河流的屬性信息完善了河流內容，便于日后系統地對河流進行管理。這里采用計算機自動編程的方法獲取河流的屬性信息。（1）流域面積:對于山地及混合河流利用watershed提取的流域分布圖，自動求算河流對應的流域面積。（2）長度:利用arcgis中的ca1cu1ategeometry功能自動計算河流長度。（3）級別:利用50km2和10～50km2河流的提取結果，根據河流編碼，自動提取河流級別。（4）岸別:利用50km2和10～50km2河流的提取結果，計算每條河流與其上級河流的夾角，根據角度大小進而左右岸判斷；（5）上級河流:利用50km2和10～50km2河流的提取結果，根據河流編碼，自動提取上級河流。（6）河源地址和河口地址:利用河流提取結果和行政界分布圖，自動獲取。（7）流經:利用河流提取結果和行政界分布圖，自動提取流經。（8）平均比降:根據河源點、河口點高程和河流長度，自動計算平均比降。（9）河源點高程和河口點高程:利用河流提取結果和DEM數據，自動獲取。（10）河流類型:利用利用河流提取結果和流域界限分布圖，自動獲取河流類型。

4 結論與分析

按照設定的技術路線，首先，分別提取山地及混合河流和平原水網河流；其次，匯總山地及混合河流和平原水網河流，生成全省10～50km2河流；最后，自動獲取河流屬性信息。

此次提取的10～50km2河流共計2581條，其中山地及混合河流2424條，平原水網河流157條。由于屬性信息表的記錄數和字段數都比較多，這里只列舉了1條記錄，如表1所示。

表1 10～50k m2河流屬性信息樣表

與現存資料進行對比分析，本文提取的結果與實際情況較為一致，但對于部分屬性信息，如河流名稱、河口地址、河源地址、流經等尚需進行進一步核對。綜上所述，本文完成了全省10～50km2河流的提取工作，獲取了河流相應的屬性信息，為我省典型河湖健康調查評價提供了數據保障。

［1］遼寧省水利廳.遼寧省水利基礎信息集［M］.北京:中國水利水電出版社，2013.

［2］國務院第一次全國水利普查領導小組辦公室.河湖基本情況普查［M］.北京:中國水利水電出版社，2010.

［3］關穎慧，鄭粉莉，王彬，等.基于DEM的黑龍江賓州河流域水系提取實驗研究［J］.水土保持通報，2012（01）.

［4］孫慶艷，余新曉，胡淑萍，等.基于Arcgis環境下DEM流域特征提取與應用［J］.北京林業大學學報，2008（11）.

［5］徐新良，莊大方，賈紹鳳，等.G I S環境下基于DEM的中國流域自動提取方法［J］.長江流域資源與環境，2004（04）.

［6］王加虎，郝振純，李麗.基于DEM和主干河網信息提取數字水系研究［J］.河海大學學報（自然科學版），2005（02）.

［7］付曉，嚴華，賀新.基于遙感圖像的河流提取方法及應用研究［J］.人民黃河，2014（03）.

［8］羅云祥，蘇保林，李卉，等.平原河網地區 SWAT模型子流域概化方法［J］.北京師范大學學報（自然科學版），2011（08）.

TV 1

：B

：1672-2469（2015）03-0023-03

10.3969/j.issn.1672-2469.2015.03.010

遼寧省政府重點項目:遼寧省河流生命現論與治理保護研究

董婷婷（1980年—），女，高級工程師。