基于DEM流域河網水系提取和面雨量計算研究

曾廣旭 曾月明 夏侯斌 謝葦華

摘要 基于DEM高程數據和ArcGIS 10.2工具軟件，分析了贛江峽江樞紐流域水系的數據預處理、水系提取和面雨量計算的原理、過程和方法，通過水系的提取確定流域控制范圍，為流域降水監測、面雨量計算和預報預警提供科學依據。

關鍵詞 DEM；ArcGIS；水系提取；面雨量；峽江樞紐

中圖分類號：TV125 文獻標識碼：B 文章編號：2095–3305（2023）09–0-03

在精細化的流域面雨量實況監測和預報預測預警氣象服務中，對流域河網的水系特征、流域水系、流域邊界和地形地貌特征等信息進行分析十分重要，流域河網的提取是流域面雨量計算、預報和水文模型建立的重要基礎數據[1-3]。基于ASTER GDEM 30 m分辨率高程DEM數據，利用ArcGIS 10.2工具軟件實現峽江水利樞紐贛江流域的河網提取。

1 流域水系提取研究

1.1 研究區概況

峽江水利樞紐位于吉安市峽江縣，處于贛江流域中游，水域從上游而下流經泰和縣、吉安縣、吉州區、青原區等縣（區），有烏江、孤江等支流匯流至峽江水利樞紐流域，是集船舶運輸、洪水防護、水利發電等于一體的梯級控制性水利樞紐。流域集雨區地形以丘陵為主，南高北低，東西兩邊地勢高。該流域屬亞熱帶丘陵山區濕潤季風氣候區，全年雨量豐沛，近30年年均降雨量1 570.0 mm，但降水量的時空分布不均勻，年降水量的72%集中在3—8月，其中，汛期4—6月雨量占年降雨量的43%。

1.2 數據源和預處理

使用ASTER GDEM 30 m分辨率數字高程數據，投影為GCS_WGS_1984，數據范圍為113.99°～116.00°E，25.99°～28.00°N，數據值域范圍為-105～2 128 m，比例尺為1∶25萬，其垂直精度和水平精度分別為30、20 m。

DEM數據因分辨率精度等會形成偽洼地現象，偽洼地會對準確計算水流方向產生一定的影響[4-5]。因此，在計算確定水流方向前，先對偽洼地的數據進行處理，將中心點逐個與緊鄰周邊點的高程值進行比較，當中心點高程值比緊鄰周邊高程值小時，則將中心點的值修改為緊鄰周邊點的最小值，將洼地柵格填平，保證所有柵格的水流方向通暢，重新生成具有一定流向無洼地的DEM高程數據（圖1）。洼地填充通過ArcGIS水文分析中的填洼工具實現。

1.3 水系的提取

1.3.1 確定水流方向 流域河網水系的提取首先應獲取水流方向，運用D8算法計算（圖2）得到水流方向，在3×3九方柵格中，用中心柵格與其緊相鄰的8個方向柵格的高程值進行計算，得出8個方向的距離權重坡度。如果是上下左右4個緊相鄰的柵格，其距離權落差坡度用計算到的高程差數值除以1.0；如果是4個斜角上相鄰的柵格，其距離權落差坡度用計算到的高程差數值除以1.414，用計算出的距離權落差逐個進行比較，判斷出水流的方向。如果距離權落差的值＞0，則該柵格單元為流出，如果＜0，則柵格單元為流入。

1.3.2 計算匯流累積量 以流向數據柵格為輸入，計算水域內匯流的累積柵格圖。其方法是先掃描每個DEM柵格單元，在流向矩陣的指導下追蹤流向，直至流向到達DEM區域的邊界，從而得出整個DEM區域的徑流累積量。通過ArcGIS水文分析中的流量工具。計算生成匯流累積柵格圖后，即可觀察所形成的河流網絡。

1.3.3 提取流域柵格 利用ArcGIS中的柵格計算器工具，將匯流累積柵格圖作為輸入數據，實現對數據的計算和處理[6]。在柵格計算器工具使用條件分析SetNull（參數1，參數2）語句，將柵格圖中滿足條件的柵格值設置為NoData，輸出所需徑流累計量的河網柵格。SetNull（參數1，參數2）有2個參數，參數1表示需要設置為NoData的條件，條件為true時返回NoData，條件為false時返回參數2指定的柵格。通過SetNul（“Flow_Accumulation” ＜ 閾值，”Flow_Accumulation”），分別設置柵格匯流閾值為1 000、3 000、5 000，輸出流域河網與地形圖疊加，并與實際進行比較，最終確定閾值為3 000與實際最為接近（圖3）。

1.3.4 河網柵格矢量化 以生成的閾值為3000的流量柵格數據作為輸入，運用ArcGIS水文分析中的河流鏈接工具，首先獲取到每個河網弧段的始末點，并計算出匯水區的出流點。其次，使用斯特拉勒方法對河網進行分級，為每個河流分配1個級別值。最后，將流量柵格圖矢量化，使用gird_code字段中的唯一值繪制直觀顯示的河網（表1、圖4）。

1.3.5 提取流域 使用編輯工具標注河網的匯流點（傾瀉點），采用捕捉傾瀉點的工具對河網流域進行分析，捕捉河網流域內累積流量最大的柵格點確定傾瀉點，再以傾瀉點為中心，一定距離范圍內搜索累積流量最大的柵格點，并將傾瀉點移動至此處。使用分水嶺工具進一步處理流域，提取出流域的分水嶺，得到贛江峽江樞紐流域圖（峽江樞紐至石虎塘航電樞紐段）（圖5）。通過ArcGIS計算出的峽江樞紐至石虎塘航電樞紐段控制的水系流域面積為1.63萬 km2。

2 流域面雨量計算分析

提取出的峽江樞紐至石虎塘航電樞紐段的水系流域包括永豐、吉水、吉州區、青原區、吉安縣、安福、永新和泰和北部部分區域。將在峽江樞紐水系流域內的174個雨量站點的經緯度信息導入ArcGIS，采用泰森多邊形法計算各站點所影響的泰森多邊形區域的面積，用多邊形包含的雨量站的降雨強度表示多邊形區域內的降雨強度，按照174個雨量站點的泰森多邊形區域面積占比和實測雨量值計算流域內面雨量（表2、圖6、圖7）。

面雨量=

式中，m表示雨量站的個數；Ri表示單個站點的降雨量；Ai表示單個站點泰森多邊形的面積；A表示流域的總集水區面積。

3 結束語

利用DEM數據和ArsGIS提取的峽

江樞紐至石虎塘航電樞紐段的水系網與現實中的流域水系和集水區面積進行驗證比較。結果表明：制作提取的水系網和集水區面積具有較強的一致性，但也有少量的偏差，存在偏差的原因主要如下：一是DEM垂直和水平數據的精度還不夠高，影響了數據的處理結果；二是ArcGIS軟件使用和計算中，其參數的設定有一定的影響，如匯流累積閾值的設定對水系提取的密度和精度均有影響。

參考文獻

[1] 湯國安，李發源，劉學軍.數字高程模型教程[M].北京：科學出版社，2010.

[2] 李翀，楊大文.基于柵格數字高程模型DEM的河網提取及實現[J].中國水利水電科學研究院學報，2004，2（3）：50-56.

[3] 曾向輝，曲樹國，馮杰.一種改進的DEM數字水系簡易提取方法研究[J].水利水電技術，2015，46（3）：17-21.

[4] 葉愛中，夏軍，王綱勝，等.基于數字高程模型的河網提取及子流域生成[J].水利學報，2005，50（5）：531-537.

[5] 楊華容，文路軍，彭文甫，等.基于DEM和GIS的流域水文信息提取：以巴中市為例[J].人民長江，2016，47（8）：34-38.

[6] 王玉德.基于ArcGIS的泰森多邊形法計算區域平均雨量[J].吉林水利，2014， 34（6）：58-60，63.

Research on Water System Extraction and Surface Rainfall Calculation of River Network Based on DEM Basin

—Taking XiaJiang Hub Basin of Ganjiang River as An Example

Zeng Guang-xu et al（Ji’an Meteorological Bureau， Ji’an， Jiangxi 343000）

Abstract Based on DEM elevation data and ArcGIS 10.2 tool software， the principles， processes， and methods of data preprocessing， water system extraction， and surface rainfall calculation in the Ganjiang Xiajiang River Hub watershed were analyzed. By extracting the water system， the control range of the watershed was determined， providing a scientific basis for accurate and precise precipitation monitoring， surface rainfall calculation， and prediction and warning in the watershed.

Key words DEM; ArcGIS; Water system extraction; Surface rainfall; Xiajiang hub