基于ArcGIS的鐵路橋涵水文分析

宋玉新

（中鐵一院集團新疆鐵道勘察設計院有限公司，新疆　烏魯木齊830011）

基于ArcGIS的鐵路橋涵水文分析

宋玉新

（中鐵一院集團新疆鐵道勘察設計院有限公司，新疆烏魯木齊830011）

鐵路橋涵水文計算是橋涵勘測設計中的重要工作?；跀底指叱棠Ｐ停―EM）和地理信息系統（ArcGIS）技術，對鐵路沿線流域進行了水文模擬分析。提取了流網、流域界線、流域面積、流域坡度等水文特征，將其與人工在地形圖上圈繪的流網和流域面積比較，發現結果一致。然后將ArcGIS得到的水文特征用大中流域洪峰流量計算模量公式計算流域洪峰流量。

鐵路橋涵水文；地理信息系統（GIS）；洪峰流量。

1　概述

20世紀50年代后期以來，DEM （DigitalElevationModel）在測繪、土木工程、地質、礦山工程、景觀建筑、道路設計、防洪、農業、規劃、軍事工程、飛行器與戰場仿真等領域取得了廣泛應用［1］。DEM（Digitalelevationmodel，數字高程模型）數據中包含了豐富的地形、地貌、水文信息，能夠反映各種分辨率的地形特征，通過DEM可以提取大量的地表形態信息，如流域網格單元的坡向、坡度以及單元格之間的關系等［2］。用DEM進行流域分析分的工具很多，ArcGIS的水文分析模塊（Hydrology model）是美國環境系統研究所（ESRI）為ArcGIS推出的水文分析模塊，主要用于地形和河流網系的提取和分析，實現地形模型可視化。運用地理信息系統（ArcGIS）中Hydrologymodel來提取水系，國內外已開展了大量的工作研究地表水文特性［3－5］。

本文的重點在于基于DEM模型應用ArcGIS的水文分析工具對流域進行降雨流線分析和流域匯水面積的計算，這對于鐵路橋涵選線、水文分析具有非常重要的實際意義和應用價值。

2　數據準備

2．1研究區域概況

研究鐵路項目位于新疆維吾爾自治區博爾塔拉蒙古自治州（簡稱：博州）和新疆生產建設兵團第五師境內。地處準噶爾盆地內，盆地內分布有較多的平緩洼地或湖泊，東部有艾比湖，西部賽里木湖。如圖1所示。

圖1　研究流域位置圖

本線地形總體呈現北高南低、西高東低，除線路起始端約三分之一路段走形于山前戈壁區外，其余路段多數走形于耕地區。本線多數河溝屬于間歇性、季節性流水河溝，平時干枯無流水，只有流域內降暴雨或溫度陡升引起的高山積雪融化才可形成洪峰，洪水出山口后迅速擴散，形成漫流或串流，或消失于耕地，匯流最終注入艾比湖。

2．2數據源

DEM是從GLCF網站上下載的數據（http∶//www. gscloud.cn/），柵格大小為90m×90m。柵格數據形式采用ESRIGRID，橫縱坐標的分辨率均為90m。

3　研究區水文特征分析

3．1洼地填補

DEM表面存在著一些凹陷的區域 （DEM本身是插值計算的，很難與現實情況完全符合），由于這些區域異常低值存在，使得該區域在進行水流流向計算時得到不合理的水流方向。因而，在進行流域地形分析以前，必須對原始的DEM進行洼地的處理。通過對洼地的處理可以生成無洼地DEM，如圖2所示。在無洼地DEM中，自然流水可以暢通無阻地流至區域地形的邊緣。因此，借助無洼地DEM可以對原數字地區進行自然流水模擬分析。洼地填充的基本過程是先利用水流方向數據計算出DEM數據中的洼地區域，并計算出其洼地深度，然后，依據這些洼地深度設定填充閥值進行洼地填充［1］。

圖2　流域DEM圖

3．2水流方向分析

流域內各網格流向的確定是利用DEM提取流域地貌特性的關鍵內容，它決定著地表徑流路徑及網格單元間流量的分配［6］。ArcGIS中的水流方向利用D8算法計算。針對每一個柵格，將其高程與周圍八個柵格進行比較，得到水流方向。

3．3匯流能力分析

匯流累積量由水流方向數據計算而來。每一個柵格匯流累積量的值代表著其上游有多少個柵格的水流方向最終匯流經過該柵格。一般而言，計算出來的匯流累積量的數值越大，代表越有可能是河谷。利用Hydrology下的FlowAccumulation函數，通過確定所有流入本單元格的累積上游單元格數目（NIP）來生成流域匯流能力柵格圖。

3．4水流網絡的提取

目前常用的河網提取方法是采用地表徑流漫流模型計算［7-9］。當匯流累積量達到一定值的時候，就會產生地表水流，那么所有那些匯流累積量大于閾值的柵格就是潛在的水流路徑，由這些水流路徑構成的網絡，就是河網。這個閾值是需要人為設定的，與定義的匯水面積有關。這里，我們設置為500，也就是說，如果在某個柵格的上游有500個柵格的水流會流經這個柵格，則將這個柵格定義為河流。由于本文用的柵格是90m×90m，也就是說，如果某個柵格點，其上游的匯水面積超過4.05km2，則認為這里是河流。圖3就是閾值分別為500個網格單元時所提取到的流域水系河網。

圖3　流域水系河網

3．5流域盆地的確定

流域盆地是由分水嶺分割而成的匯水區域。它通過對水流方向數據的分析確定出所有相互連接并處于同一流域盆地的柵格。流域盆地的確定首先是要確定分析窗口邊緣的出水口的位置，也就是說，在進行流域盆地的劃分中，所有的流域盆地的出水口均處于分析窗口的邊緣。當確定了出水口的位置之后，也就找出了所有流入出水口的上游柵格的位置［10］。

3．6子流域的生成

經過上一步得到的流域盆地是一個比較大的流域盆地，在很多的水文分析中，還需要基于更小的流域單元進行分析，那么就需要進行流域的分割。而流域的分割首先是要確定小級別流域的出水口位置［10］。

1）Streamlink的生成。Stream？link記錄著河網中節點之間的連接信息。Stream？link的每條弧段要么連接著兩個作為出水點或匯合點的結點，要么連接著作為出水點的結點和河網起始點。通過Streamlink的計算，即得到每一個河網弧段的起始點和終止點。同樣，也可以得到該匯水區域（流域）的出水口。

2）集水流域的生成。通過Streamlink作為流域的出水口數據所得到的集水區域是每一條河網弧段集水區域，也就是要研究的最小溝谷的集水區域，它將一個大的流域盆地按照河網弧段分為一個個小的集水盆地。將流域柵格轉換成為矢量圖層，并進行符號設置，得到的結果如圖4所示。

圖4　流域集水盆地

3．7坡度分析

利用ArcGIS軟件的3DAnalyst（三維分析）功能模塊，通過矢量等高線數據生成TIN三角格網，最后轉換成Slop（坡度）規則格網數據，對該數據層進行分析和統計匯總。利用ArcGIS分析功能分析流域坡度組成情況，分析結果精度高。不但可以滿足用戶在任意精度、任意尺度下對成果的要求，而且用戶確定分析精度后，所產生的分析結果具有唯一性，避免分析結果的隨機性、多樣性所產生的誤差［11］。流域的坡度分析云圖如圖5所示。

圖5　流域的坡度分析云圖

4　流量計算

大流域經驗公式的簡介及本線選用情況經驗公式一般采用區域回歸方程的形式，確定了水文子區域后，便可以建立流量與地形和氣候因子的區域回歸方程，繼而可以估計無資料點的流量。常用的經驗公式有［12］∶

式中∶Qp%為設計流量，m3/s；F為流域面積km2；C為隨地區和頻率變化的綜合系數；L為流域長度，km；i為河道干流平均坡度；f為流域形狀系數 （f=F/ L2）；a、b、m、n為經驗指數。

結合水文站及收集到的大中河流資料相對較少，根據搜集到的資料情況及多年水文經驗，采用Qp%=CFn進行大流域模量公式推算。根據在博樂地區收集到的相關水庫資料及沿線大中河流、水文站歷年最大流量系列資料，進行整理分析，并據此確定模量公式Q1%=13·F0.5。具體計算結果見表1。

表1　流域面積與流量

由表1可見，基于GIS確定的流域與人工繪制的流域結果一致，可將其應用于鐵路橋涵水文分析。

5　結論與討論

本研究成功地運用ArcGIS的水文分析模塊進行了流域數字高程模型的處理及水文特征分析。對鐵路沿線流域的地表水文特征分析研究的結果表明，利用DEM模擬的水系的空間分布、數字流域界線與實際分布情況基本相符。為檢驗結果的精確性，將ArcGIS提取的與人工圈繪的水系流域面積進行對比，結果一致。

［1］李志林.數字高程模型［M］.武漢測繪科技大學出版社，2000.

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

［3］JensonSK.Applicationsofhydrologicinformationautomatically extracted from digital elevation models［J］.Hydrologica lProcesses，1991，5（1）∶31-44.

［4］DGTarboton，BrasRL，Rodriguez-IturbeI.Ontheextraction ofchannelnetworksfromdigitalelevationdata［J］.Hydrological Processes，1991，5（1）∶81-100.

［5］陳永良，劉大有，虞強源.從DEM中自動提取自然水系［J］.中國圖象圖形學報，2002，7（1）∶91-96.

［6］趙健，賈忠華，羅紈.ARCGIS環境下基于DEM的流域特征提?。跩］.水資源與水工程學報，2006，17（1）∶74-76.

［7］孫慶艷，余新曉，胡淑萍，等.基于ArcGIS環境下DEM流域特征提取及應用［J］.北京林業大學學報，2008（S2）∶144-147.

［8］顧用紅，舒振文，張火青.DEM在流域水文特征分析中的應用［J］.人民珠江，2001（4）∶5-6.

［9］易紅偉，湯國安，劉詠梅，等.河網徑流節點及其基于DEM的自動提?。跩］.水土保持學報，2003，17（3）∶108-111.

［10］鄧必平，嚴恩萍，洪奕豐，等.基于GIS和DEM的東江湖流域水文特征分析［J］.湖北農業科學，2013，52（15）∶3531-3536.

［11］滕利強，王亮.ArcGIS空間分析功能在流域坡度分析中的應用［J］.中國水土保持，2008（4）∶40-41.

［12］徐輝元.新建鐵路克拉瑪依至塔城線大流域水文計算辦法探討［J］.公路工程，2014（2）∶299-302.

U442.3