地理信息系統(tǒng)在水文分析中的應(yīng)用探討

王文民

（泰安市水文局，山東 泰安271000）

1 研究區(qū)概況

大汶河古稱汶水，發(fā)源于旋崮山北麓沂源縣境內(nèi)。流域內(nèi)地勢東高西低，自東向西流經(jīng)萊蕪、新泰、肥城、寧陽、汶上、東平等縣（市），匯入東平湖。干流河道長239 km，流域面積9 098 km2。

2 行程時間方程

行程時間方程用于估計在分水嶺內(nèi)的特定點的行程時間。行程時間的影響因素包括流動路徑的長度和坡度、地面粗糙度、降雨強度等。由于水文模型的發(fā)展，出現(xiàn)了兩種不同的行程時間分類，分別為集中時間和滯后時間。

式中：TLAG為滯后時間；TC為集中時間。

2.1 滯后時間方程

滯后時間被定義為從降雨過剩的質(zhì)心到流域峰值徑流的時間。用于從地理數(shù)據(jù)中估計滯后時間的經(jīng)驗公式為：

式中，C為計算流域坡度和蓄水量的差異（通常由不同類別流域的實驗數(shù)據(jù)確定）；L為從參考點到流域上游邊界沿主河道的流域長度；Lca為從參考點到質(zhì)心對面點沿主河道的距離；S為坡度；m為滯后指數(shù)。

徑流參數(shù)被“集中”到一個方程中，以概括整個流域的流動路徑和徑流，大部分參數(shù)可變性由系數(shù)和指數(shù)表示。實際上，流動路徑將穿過不同坡度、不同水力條件的區(qū)域。該滯后時間方程對于同質(zhì)盆地可能是準確的，但對于異質(zhì)盆地則需要劃分區(qū)域或使用不同的方程組。

2.2 集中時間方程

集中時間的定義是水的微粒到達分水嶺中的排放點所需的最長行進時間。常用的一組方程式是TR-55方程組，分為3類流動，分別為片流、淺層集中流和渠道流。

片狀流是河流源頭的流動，該段流路應(yīng)小于300 m：

式中Tt為以小時為單位的流動時間；n為曼寧粗糙度系數(shù)；L為流路長度；P2為日降雨量；S為水力坡度線的坡度。

在初始片流之后，使用描述淺集中流的方程，直到流徑可以表示為明渠流：

式中V為平均速度。

淺集中流的速度：

式中C為常數(shù)；S為平均坡度。

最后，一旦流動路徑到達通道、路緣和排水溝，或其他可以使用明渠流動方程估算行程時間的水力條件，可用曼寧方程表示計算方程（4）使用的速度：

式中k為單位因子；n為曼寧粗糙度；R為水力半徑；S為坡度。

理論上，用于計算集中時間方程的實際流路可以分成任意數(shù)量的流路段，每個流路段都有自己的方程，以便估計總的行進時間。

3 地理信息系統(tǒng)矢量數(shù)據(jù)線的應(yīng)用

地理信息系統(tǒng)矢量數(shù)據(jù)是地理標志及其相關(guān)數(shù)據(jù)庫屬性的組合。由于一條給定的線可能由幾個線段組成，斜率可通過取兩個頂點的高差，除以該線段的長度，對所有線段斜率求平均值而得到。使用地理信息系統(tǒng)軟件可以很容易地計算出來，且線段、長度、坡度、粗糙度等關(guān)鍵參數(shù)可以在地理信息系統(tǒng)數(shù)據(jù)庫中進行維護，如圖1所示。用于計算滯后時間或集中時間的不同方程中所使用的參數(shù)可存儲為線段后存入數(shù)據(jù)庫。然后，可以開發(fā)計算機應(yīng)用程序或腳本，以存儲的數(shù)據(jù)以及作為附加屬性存儲的線段為基礎(chǔ)求解行進時間。

圖1 Arc GIS中的流程線定義和屬性

4 流域建模系統(tǒng)的實現(xiàn)

流域建模系統(tǒng)（WMS）是一個綜合工具，用于組織空間和水文數(shù)據(jù)，可以使用數(shù)字圖像、矢量地理信息矢量數(shù)據(jù)和網(wǎng)格數(shù)據(jù)來圈定盆地，還支持導(dǎo)入、導(dǎo)出和創(chuàng)建地理信息矢量數(shù)據(jù)。水流路徑從最陡下降路徑推導(dǎo)出來，可以通過每個子流域自動生成，然后可以對這些導(dǎo)出的流動路徑段確定長度和斜率。

1）創(chuàng)建流動路徑。在WMS中，河流可以被定義為一系列的三角形邊緣或一個具有相對較大上游貢獻面積的數(shù)字高程模型單元的集合。將高程數(shù)據(jù)導(dǎo)入WMS，可以用于沿著最陡下降路徑生成流路線。

2）確定流路線節(jié)點和頂點的高程。根據(jù)節(jié)點處的高差，并除以流路線的投影長度計算得出坡度。當不使用背景數(shù)字地形模型來提供高程數(shù)據(jù)而創(chuàng)建流動路徑時，必須從地形圖中單獨確定坡度，然后手動輸入。長度是從流路線本身導(dǎo)出的，可以手動更改或編輯。

3）確定旅行時間方程。將用于計算沿流路線行駛時間的方程式分配給每個路段盆地或分水嶺。在WMS中，從數(shù)字地形模型計算的或在其他方程中可用的任何變量都可以作為變量分配給用戶定義的方程作為等式的一部分，其他變量必須作為常數(shù)輸入。

4）輸入其他方程參數(shù)。定義流路線并指定方程式后，必須定義與所選方程式相關(guān)的參數(shù)。WMS將自動輸入根據(jù)線和數(shù)字地形數(shù)據(jù)計算的坡度和長度值。其他參數(shù)單獨輸入，或根據(jù)包含此信息的地理信息系統(tǒng)數(shù)據(jù)層計算得出。

5）計算流域或子流域的集中時間。在方程和參數(shù)被定義之后，每個弧的行進時間可以通過計算得出。滯后時間和集中時間在一個盆地內(nèi)通過對從不同弧線計算的行進時間求和來確定。由于WMS可用于確定流域邊界和開發(fā)水文模型，每個流徑線可根據(jù)其地理位置分配給一個流域。

5 結(jié)論

本文提出了一種利用地理信息系統(tǒng)進行行程時間方程輸入的新方法，并以大汶河支流為例證明了這種確定方程參數(shù)的方法是有效的。將行進時間線進行分組，獲得流域的單個行進時間值，然后在數(shù)值模型中估計滯后時間或集中時間，通過合成單位線方法計算流域徑流的方法，可以在流域建模系統(tǒng)中進行推廣應(yīng)用。