基于降雨-徑流集成過(guò)程的分布式洪水淹沒(méi)模型

張漢濤

(梅州市大埔韓江高陂水利樞紐工程建設(shè)管理處， 廣東 梅州 514000)

1 概述

洪水是最具破壞性的自然災(zāi)害之一[1]。近幾十年來(lái)，由于全球變暖和大規(guī)模的城市化，城市洪水災(zāi)害變得更加頻繁，且難以預(yù)測(cè)和防范，危害性日益加劇[2]。遙感技術(shù)具有快速、大空間范圍、低成本和高效的對(duì)地探測(cè)能力，在洪水災(zāi)害監(jiān)測(cè)中具有巨大的應(yīng)用潛力[3]。隨著GIS和遙感技術(shù)在水資源管理中的應(yīng)用，越來(lái)越多的研究者參與到分布式水文模型的開(kāi)發(fā)中來(lái)[4-5]。水力模型提供了模擬洪水的基礎(chǔ)，特別是淹沒(méi)水動(dòng)力數(shù)值模型在城市地區(qū)的開(kāi)發(fā)利用[6]。

因此，建立一個(gè)適用于任何流域的洪水災(zāi)害綜合評(píng)估模型具有重要意義。本文提出了一個(gè)基于降雨-徑流過(guò)程的分布式洪水淹沒(méi)模擬模型，旨在對(duì)流域范圍內(nèi)的洪水災(zāi)害進(jìn)行評(píng)估分析。

2 研究方法

2.1 模型概況

洪水的時(shí)空信息對(duì)于減輕和評(píng)估洪水災(zāi)害至關(guān)重要[7]。本研究建立了一個(gè)結(jié)合降雨-徑流過(guò)程的洪水淹沒(méi)模擬模型，該模型概化模擬了洪水水深和淹沒(méi)范圍。

源于降水或土壤飽和過(guò)量的地表徑流流入河道。當(dāng)河道內(nèi)水位超過(guò)平灘水位時(shí)，水從河道流出，流入鄰近的地表區(qū)域[8-10]，如圖1所示。因此，這些地表區(qū)域要么充當(dāng)洪水的臨時(shí)蓄水池，要么充當(dāng)輸送水流的路線。在此過(guò)程中，河道水面高程采用一維河道水流模型計(jì)算。整個(gè)流域的地表徑流路徑由一個(gè)二維模型模擬，該模型將每個(gè)網(wǎng)格視為一個(gè)存儲(chǔ)單元，單元體積隨時(shí)間的變化等于降雨過(guò)剩量和進(jìn)出的通量。在河道和地表區(qū)域之間的界面上，根據(jù)相對(duì)水面高度簡(jiǎn)單地模擬了水流交換，也就是說(shuō)，地表水流可以遠(yuǎn)離河流，也可以流入河流，如圖2所示。圖中d為單元中水深(m)，g為地面高度(m)，h為水面的高度(m)。

圖1 洪水演進(jìn)二維示意圖

圖2 網(wǎng)格單元1和2之間的流動(dòng)

地表徑流和河流演進(jìn)的過(guò)程可以用圣維南方程的來(lái)近似模擬[11]。在一些基于運(yùn)動(dòng)波的模擬模型中，網(wǎng)格水流的流向由數(shù)字高程模型數(shù)據(jù)生成。然而，它不能處理回水情況。因此，擴(kuò)散波近似法被用作陸地和河道水流的控制方程。

2.2 計(jì)算方程

水文和水力模型需要清楚地了解集水區(qū)的水文循環(huán)[12-14]。本文提出的模型基于網(wǎng)格數(shù)據(jù)結(jié)構(gòu)，采用網(wǎng)格單元描述分布式流域特征，如高程、土地利用類(lèi)型、土壤類(lèi)型等。其中針對(duì)每個(gè)網(wǎng)格單元模擬了與降雨-徑流過(guò)程相結(jié)合的洪水淹沒(méi)情況。

假設(shè)單元上的降雨過(guò)剩量Re等于降雨量Rg減去攔截量Vi和入滲量Vf，單位均為mm，蒸散發(fā)沒(méi)有考慮在內(nèi)。因此，降雨過(guò)剩量由式(1)計(jì)算。

Re=Rg-Vi-Vf

(1)

在模擬過(guò)程中，如果降雨強(qiáng)度小于攔截強(qiáng)度，則將降雨率設(shè)置為零，并在隨后的時(shí)間中從降雨量中抵扣剩余的攔截量，直到滿足攔截水深。假設(shè)土壤均勻、深厚且在每個(gè)網(wǎng)格單元內(nèi)排水良好，則滲透量是從單層格林-安普特滲透模型計(jì)算得出。

地表徑流由連續(xù)性和動(dòng)量方程控制。采用顯式有限差分法求解地表徑流方程。為了求解與流入或流出網(wǎng)格單元及其體積變化相關(guān)的連續(xù)性方程，根據(jù)曼寧公式計(jì)算每個(gè)方向的動(dòng)量方程，見(jiàn)式(2)和式(3)。

(2)

(3)

式中，hi，j—網(wǎng)格單元的自由水面高度，m；t—時(shí)間，s；△x、△y—網(wǎng)格單元尺寸；n—曼寧摩擦系數(shù)；Qx、Qy—網(wǎng)格單元之間在x和y方向上的流量，m3/s，其符號(hào)取決于流動(dòng)方向；Q上、Q下、Q左、Q右—來(lái)自上、下、左、右相鄰網(wǎng)格單元的流量，m3/s。

(4)

(5)

式中，x—距離，m；A—橫截面積，m2；Q—流量，m3/s；q—橫向流入或流出河道的流量，m3/s；p—渠道的濕周，m；S0—河床坡降。

2.3 模型開(kāi)發(fā)

本研究的主要目的是開(kāi)發(fā)一個(gè)綜合的洪水淹沒(méi)模擬模型，以預(yù)測(cè)和分析洪水災(zāi)害。采用一種松耦合的程序開(kāi)發(fā)方法。在這種情況下，GIS和洪水淹沒(méi)模擬模型保持分離，但通過(guò)數(shù)據(jù)輸入和輸出操作松散地聯(lián)系在一起，如圖3所示。

圖3 模型開(kāi)發(fā)方案

模型程序是在微軟Visual C++和Windows操作系統(tǒng)環(huán)境下，根據(jù)上述方法開(kāi)發(fā)的，其中考慮了降水、攔截和滲透等水文過(guò)程。

模型應(yīng)用于韓江流域。韓江發(fā)源于汕尾陸豐與河源紫金交界的烏突山七星崠，韓江全長(zhǎng)470km，流域面積30112km2，平均坡降為0.39‰。韓江在梅州市境內(nèi)有集雨面積14711km2，河長(zhǎng)343km，總落差164m，水力資源豐富。韓江上游干流為梅江，主要支流汀江，梅江與汀江總集水面積占高陂樞紐集水面積的90%以上[15]。

3 模型應(yīng)用

3.1 數(shù)據(jù)獲取

對(duì)于地面徑流或河道水流的計(jì)算和淹沒(méi)區(qū)域的預(yù)測(cè)，需要模型區(qū)域的高程信息。為了模擬大面積水文過(guò)程，以200m的分辨率對(duì)韓江流域的地形進(jìn)行了提取、投影和重采樣，其中行數(shù)和列數(shù)分別為259、250。SRTM數(shù)字移動(dòng)終端是基于干涉測(cè)量法的，它包括通過(guò)兩個(gè)天線接收反向散射的雷達(dá)信號(hào)。尤其是在洪泛區(qū)，這些信號(hào)包括雷達(dá)散斑或隨機(jī)噪聲，會(huì)降低相對(duì)垂直精度。因此，采用數(shù)據(jù)濾波去除數(shù)據(jù)斑點(diǎn)來(lái)提高SRTM質(zhì)量。

對(duì)于河道模型參數(shù)化，將河道橫截面設(shè)定為梯形。沿河道的每個(gè)網(wǎng)格單元所需的參數(shù)是河岸高度、河道寬度和坡度，這些參數(shù)可以在每個(gè)網(wǎng)格中指定。河道參數(shù)是根據(jù)有關(guān)韓江河流域的文件提供的。河道寬度為20～150m不等，見(jiàn)表1。降雨數(shù)據(jù)如圖4所示，通過(guò)插值法得到每個(gè)網(wǎng)格單元的降雨量。地表網(wǎng)格單元的曼寧系數(shù)值和攔截量是根據(jù)土地利用分類(lèi)規(guī)定的，見(jiàn)表2。

圖4 韓江河流域降雨情況

表1 河道屬性示例

表2 不同土地利用分類(lèi)的曼寧系數(shù)和攔截量

3.2 結(jié)果與分析

綜合模型被用于研究韓江流域的一次洪水案例。模擬結(jié)果顯示了洪水淹沒(méi)的最終空間分布，其中觀測(cè)到的洪水淹沒(méi)面積約為41km2，模擬淹沒(méi)面積為42.49km2，結(jié)果符合良好，表明該模型可以模擬洪水并識(shí)別洪水易發(fā)區(qū)。

根據(jù)模擬的淹沒(méi)范圍和觀測(cè)的淹沒(méi)范圍的疊加分析，模擬淹沒(méi)范圍位于觀測(cè)淹沒(méi)范圍內(nèi)的面積約30.22km2，一致性率約為75%。在模擬中，河道由數(shù)字高程模型生成，與實(shí)際河流略有不同，因此對(duì)洪水淹沒(méi)的空間分布有影響，不容易直接比較洪水淹沒(méi)的空間分布。然而，根據(jù)對(duì)流域洪水泛濫的空間分布的分析，模擬誤差在可接受范圍內(nèi)。

模擬和觀測(cè)得到的河道水深對(duì)比如圖5所示。水深峰值及其變化得到了很好的模擬展現(xiàn)，反映了該模型的有效性。模擬水深下降階段沒(méi)有觀測(cè)的下降速度快，這是因?yàn)槟Ｐ湍壳皼](méi)考慮潰堤等其他因素。

圖5 水深-降雨對(duì)比

4 結(jié)論

洪水災(zāi)害是世界上危害最大的災(zāi)害之一，獲取洪水特性信息對(duì)防災(zāi)減災(zāi)具有重要意義。本文提出了一種基于降雨-徑流過(guò)程的分布式洪水淹沒(méi)模擬模型，對(duì)模型進(jìn)行介紹并應(yīng)用于韓江流域，模型模擬結(jié)果表明模型整體性能良好，為洪水淹沒(méi)模擬提供了一種可行的建模工具，對(duì)洪水災(zāi)害評(píng)估，尤其是對(duì)溢流洪水的評(píng)估具有重要意義。

然而，模型的精度不僅取決于地形圖的精度，還取決于水文模型和水力模型的結(jié)構(gòu)。目前，水文模型只考慮了地表水流，簡(jiǎn)化了水力過(guò)程。為了提高模型對(duì)各種水文環(huán)境的適用性，在進(jìn)一步的研究中需要將地下水流和其他詳細(xì)過(guò)程模型耦合起來(lái)。