華北地區(qū)分布式蓄超空間動(dòng)態(tài)組合TOKASIDE-D模型研究

劉玉環(huán)，劉志雨,2，李致家，張 珂，劉開(kāi)磊

(1. 河海大學(xué)水文水資源學(xué)院，江蘇 南京 210098；2. 水利部信息中心(水利部水文水資源監(jiān)測(cè)預(yù)報(bào)中心)，北京 100053；3. 淮河水利委員會(huì)水文局(信息中心)，安徽 蚌埠 233000)

洪水是中國(guó)最常見(jiàn)的自然災(zāi)害之一，特別是在中國(guó)華北地區(qū)(半濕潤(rùn))，陡漲陡落的洪水極易造成嚴(yán)重的生命和財(cái)產(chǎn)損失[1-3]。在實(shí)際生產(chǎn)中，洪水防治重要的非工程措施之一就是利用水文模型進(jìn)行洪水預(yù)報(bào)。研究發(fā)現(xiàn)，大部分水文模型在濕潤(rùn)地區(qū)應(yīng)用效果較好，但是在半濕潤(rùn)半干旱地區(qū)模擬效果不佳[4-6]。這是由于這類(lèi)地區(qū)的下墊面和降水時(shí)空分布不均勻，超滲產(chǎn)流和蓄滿(mǎn)產(chǎn)流隨時(shí)空變化的現(xiàn)象尤為明顯[7]。另外，華北地區(qū)秋冬季降雨較少，用水量較多，地下水開(kāi)采量大于補(bǔ)給量，地下水位下降。汛期剛開(kāi)始時(shí)，地下水埋深較深，徑流在匯流過(guò)程中不斷滲漏，用于抬升地下水位[8]。因此，華北地區(qū)的洪水預(yù)報(bào)比濕潤(rùn)地區(qū)更具挑戰(zhàn)性[9]。

目前，經(jīng)過(guò)國(guó)內(nèi)外學(xué)者多年的研究，在通過(guò)大量的模型改進(jìn)后，混合產(chǎn)流問(wèn)題取得了階段性成果，從單一產(chǎn)流模型如CASC2D模型[10]、SCS模型[11]、TANK模型[12]等，到混合產(chǎn)流模型包括垂向混合徑流模型[13-14]、增加超滲的新安江模型[15-16]等，專(zhuān)門(mén)針對(duì)海河流域的新安江-海河模型以及河北雨洪模型[17]等，再到靈活框架模型如FUSE[18]、FARM[19]及SCCM[20]等。但是，這些改進(jìn)后模型的產(chǎn)流模式在計(jì)算時(shí)都是固定不變的，沒(méi)有系統(tǒng)地動(dòng)態(tài)識(shí)別洪水過(guò)程中產(chǎn)流模式隨時(shí)間變化的情況。

李致家等[21]認(rèn)為混合產(chǎn)流模擬的主要困難是識(shí)別隨時(shí)空變化的產(chǎn)流機(jī)制。對(duì)于一個(gè)特定流域的某場(chǎng)洪水來(lái)講，流域地形地貌、土壤屬性、植被等下墊面特征是基本不變的，但降雨和土壤含水量等因子會(huì)隨時(shí)間而動(dòng)態(tài)變化，并且這類(lèi)因子的變化與產(chǎn)流模式的空間分布是相互關(guān)聯(lián)的。因此，針對(duì)這一現(xiàn)象，本文在TOKASIDE模型(topgraphic kinematic approximation and saturation-infiltration double excess grid-based distributed model)的基礎(chǔ)上，通過(guò)降雨和土壤含水量之間的判定關(guān)系，動(dòng)態(tài)識(shí)別蓄滿(mǎn)產(chǎn)流和超滲產(chǎn)流網(wǎng)格，并添加滲漏模塊，構(gòu)建基于物理基礎(chǔ)的分布式蓄超空間動(dòng)態(tài)組合模型(TOKASIDE-D)，并在我國(guó)華北地區(qū)的北辛店流域進(jìn)行驗(yàn)證和分析，探究流域產(chǎn)流區(qū)蓄滿(mǎn)、超滲產(chǎn)流模式空間動(dòng)態(tài)組合方法的應(yīng)用效果。

圖1 TOKASIDE-D模型結(jié)構(gòu)Fig.1 Structure diagram of TOKASIDE-D model

1 模 型 構(gòu) 建

1.1 TOKASIDE 模型

TOKASIDE模型是劉志雨等[22]提出的基于地形與運(yùn)動(dòng)波的蓄滿(mǎn)超滲產(chǎn)流機(jī)制的網(wǎng)格化分布式模型。該模型是在TOPKAPI模型的基礎(chǔ)上發(fā)展起來(lái)的，進(jìn)一步考慮了地下水運(yùn)動(dòng)、水庫(kù)入流與調(diào)蓄計(jì)算等模塊以及超滲產(chǎn)流機(jī)制。TOKASIDE模型在空間上采用非線(xiàn)性運(yùn)動(dòng)波方程法，將每一個(gè)計(jì)算單元中的水文過(guò)程概化為3個(gè)“結(jié)構(gòu)上相似”的非線(xiàn)性水庫(kù)方程，用以描述水流在土壤層、飽和坡面以及河道中運(yùn)動(dòng)波方程的融合。3個(gè)方程分別代表上層土壤中的排水、飽和土壤(或不透水土壤層)中的地面徑流和河道徑流。

1.2 TOKASIDE-D模型

TOKASIDE模型在產(chǎn)流計(jì)算模塊中雖然考慮了流域超滲產(chǎn)流，但是這一考慮是針對(duì)流域內(nèi)所有網(wǎng)格，即扣除蒸發(fā)后的凈雨先完成超滲過(guò)程，后繼續(xù)下滲進(jìn)行蓄滿(mǎn)過(guò)程，這與實(shí)際情況中蓄滿(mǎn)產(chǎn)流與超滲產(chǎn)流交替主導(dǎo)的現(xiàn)象不符。在華北地區(qū)開(kāi)展徑流模擬探究時(shí)，需要對(duì)蓄滿(mǎn)產(chǎn)流和超滲動(dòng)態(tài)交替以及地下水超采的情形進(jìn)行考慮。針對(duì)這一問(wèn)題，對(duì)TOKASIDE模型進(jìn)行改進(jìn)，彌補(bǔ)假設(shè)條件中缺失的情況，動(dòng)態(tài)識(shí)別蓄滿(mǎn)產(chǎn)流和超滲產(chǎn)流網(wǎng)格，以提高模型在華北地區(qū)的模擬精度。

在TOKASIDE模型基礎(chǔ)上，通過(guò)降雨和土壤含水量之間的判定關(guān)系，動(dòng)態(tài)識(shí)別蓄滿(mǎn)產(chǎn)流和超滲產(chǎn)流網(wǎng)格，添加滲漏模塊，構(gòu)建基于物理基礎(chǔ)的TOKASIDE-D模型，模型結(jié)構(gòu)見(jiàn)圖1(圖中p為降雨量，pe為降雨強(qiáng)度，θt為網(wǎng)格土壤含水率，θS為網(wǎng)格田間持水率，f為穩(wěn)定下滲率，Q為模擬流量)。將流域網(wǎng)格初步劃分為“蓄滿(mǎn)網(wǎng)格”和“超滲網(wǎng)格”兩類(lèi)。蓄滿(mǎn)網(wǎng)格和超滲網(wǎng)格分別采用TOKASIDE的土壤水非線(xiàn)性水庫(kù)方程和Green-Ampt下滲公式[23]計(jì)算，另外針對(duì)地下水超采出現(xiàn)的滲漏問(wèn)題，在模型中添加滲漏模塊。

1.2.1 基本假設(shè)

基于TOKASIDE模型的假設(shè)，根據(jù)華北地區(qū)下墊面情況，現(xiàn)添加以下假設(shè)條件：(a)網(wǎng)格內(nèi)的降雨、土壤性質(zhì)、土地利用等特征均勻分布，不存在下滲能力分布曲線(xiàn)。(b)土壤分成2層(飽和層、下滲層)計(jì)算。我國(guó)北方半干旱地區(qū)的特色是土層厚，年降雨量?jī)H在400～800 mm，并且多集中于7—9月，所以其土壤基本不可能達(dá)到飽和狀態(tài)。但其表層土壤還是會(huì)達(dá)到飽和狀態(tài)，所以此處將土壤分成2層進(jìn)行模擬計(jì)算[24]。蓄滿(mǎn)和超滲發(fā)生在飽和層(土壤水模塊)；下滲層則為無(wú)限下滲土層(滲漏模塊)，主要針對(duì)地下水超采形成的滲漏問(wèn)題。(c)蓄超網(wǎng)格劃分：河道附近地勢(shì)平坦，土壤水含量相對(duì)較高，河道網(wǎng)格(三級(jí)以上)全部為蓄滿(mǎn)網(wǎng)格；非河道網(wǎng)格可為蓄滿(mǎn)網(wǎng)格或超滲網(wǎng)格。(d)蓄超網(wǎng)格轉(zhuǎn)化：非河道網(wǎng)格在計(jì)算過(guò)程中，通過(guò)判斷網(wǎng)格土壤含水量是否達(dá)到田間持水量、降雨是否超過(guò)下滲率，動(dòng)態(tài)調(diào)整網(wǎng)格的蓄滿(mǎn)或超滲屬性。

1.2.2 土壤水模塊修改

在原始TOKASIDE模型中，網(wǎng)格計(jì)算分為非河道網(wǎng)格和河道網(wǎng)格。從圖1可以看出，在2種網(wǎng)格內(nèi)，均需要完成土壤水、地表水計(jì)算，河道網(wǎng)格另需進(jìn)行河道匯流演算。因此，產(chǎn)流機(jī)制的轉(zhuǎn)換主要發(fā)生在非河道網(wǎng)格的土壤水模塊。

降雨未發(fā)生時(shí)，在計(jì)算過(guò)程中當(dāng)非河道網(wǎng)格土壤含水率未達(dá)到田間持水率時(shí)，為超滲網(wǎng)格，達(dá)到田間持水率后轉(zhuǎn)為蓄滿(mǎn)網(wǎng)格。降雨發(fā)生時(shí)，降雨強(qiáng)度超過(guò)網(wǎng)格下滲率為超滲網(wǎng)格，若網(wǎng)格為蓄滿(mǎn)則修改為超滲；否則原網(wǎng)格屬性不變。蓄滿(mǎn)產(chǎn)流主要包含飽和地面徑流和壤中流，而超滲產(chǎn)流計(jì)算產(chǎn)生超滲地表徑流。

a. 當(dāng)土壤含水率未達(dá)到田間持水率時(shí)(θ<θs)，計(jì)算網(wǎng)格為超滲網(wǎng)格：

(1)

式中：Vs——網(wǎng)格單元的土壤含水量，m3；X——網(wǎng)格單元尺寸，m；Quo、Qus——從上游區(qū)域進(jìn)入當(dāng)前網(wǎng)格單元i的地表、地下徑流量，m3；hout——時(shí)段末土壤水深度，m；Ks——飽和水力導(dǎo)水率，m/s；θ0——初始土壤含水率，m3/m3；Sf——濕潤(rùn)鋒面處土吸力，m；I——累積下滲量，m。

如果降雨強(qiáng)度pe≤f，超滲地表徑流為hexf=0(hexf為超滲土壤水出流深度，m)，下滲量全部補(bǔ)充土壤含水量：

hout=hin+pet-hexf

(2)

如果pe>f，超滲地表徑流為hexf=pe-f，下滲量全部補(bǔ)充土壤含水量至田間持水量，產(chǎn)生壤中流hsoil：

hsoil=hout-hmax

(3)

式中：t——時(shí)間，s；hin——時(shí)段初土壤水深度，m；hmax——土壤水最大深度，m。

b. 當(dāng)土壤含水率達(dá)到田間持水率時(shí)(θ≥θs)，計(jì)算網(wǎng)格轉(zhuǎn)換為蓄滿(mǎn)網(wǎng)格，其計(jì)算流程如下：

(4)

式中：Cs——當(dāng)?shù)氐膫鲗?dǎo)率，m/s；αs——土壤特性決定的參數(shù)。

時(shí)段初土壤水深度hin=hmax；壤中流hsoil=hout-hmax；飽和地表徑流hexf=hout-hsoil-hmax；更新時(shí)段末土壤水深度hout=hmax。

c. 當(dāng)降雨強(qiáng)度pe>f時(shí)，蓄滿(mǎn)網(wǎng)格修改為超滲網(wǎng)格，則計(jì)算模塊由蓄滿(mǎn)模塊轉(zhuǎn)為超滲模塊，按照式(1)～(3)進(jìn)行計(jì)算。

1.2.3 滲漏模塊

華北地區(qū)地下水超采嚴(yán)重，地下水位埋深大，建模時(shí)考慮垂向土壤深層的滲透。當(dāng)土壤水的含水量超過(guò)蓄水能力時(shí)，穿過(guò)飽和層以上的厚土壤層進(jìn)行垂直向下運(yùn)移。假定上層土壤層的下滲率是土壤水含量的一個(gè)函數(shù)，則下滲水量為

(5)

式中：Pr——下滲水量，m/s；Ksv——土壤垂向飽和導(dǎo)水率，m/s；θsat——土壤飽和含水量,m3/m3；αp——取決于土壤類(lèi)型的指數(shù)(沙石:αp?11；泥土:αp?25)。

1.2.4 地表水模塊

土壤表層飽和后的剩余降水構(gòu)成地表水模塊的輸入，地表坡面徑流計(jì)算與土壤部分的描述相似，公式為

(6)

式中：Vo——地表徑流量，m3；ro——由水量平衡方程求解得到的余下水量(或土壤水回流量)，m/s；Co——與地表水流曼寧公式相關(guān)的系數(shù)；αo——使用曼寧公式得出的指數(shù)，取3/5。

1.2.5 河道演算模塊

河道徑流的運(yùn)動(dòng)波近似方程采用運(yùn)動(dòng)學(xué)方法進(jìn)行描述，在這種方法中動(dòng)量方程通過(guò)曼寧公式近似寫(xiě)成：

(7)

式中：Vc——單元網(wǎng)格的河道水量，m3；rc——旁側(cè)入流，包括地表徑流量和土壤水流量，m2/s；Quc——來(lái)自上游單元河道網(wǎng)格的流量，m3/s；S0——河底坡降；nc——河道曼寧糙率，s/m2/3；Cx——濕周，m。

1.2.6 模型參數(shù)

TOKASIDE-D模型的主要參數(shù)包括計(jì)算單元土壤類(lèi)型對(duì)應(yīng)的土壤厚度L(m)、土壤橫向飽和水力導(dǎo)水率Ksh(m/s)、垂向飽和水力導(dǎo)水率Ksv(m/s)，土地利用類(lèi)型對(duì)應(yīng)的地表曼寧糙率ns，河道分級(jí)對(duì)應(yīng)的河道曼寧糙率nc等。Green-Ampt下滲公式的參數(shù)：飽和水力導(dǎo)水率Ks(m/s)，濕潤(rùn)鋒面處的土壤吸力Sf(m)。在構(gòu)建TOKASIDE-D模型時(shí)，將TOKASIDE模型與Green-Ampt下滲公式中物理含義相同的參數(shù)合并，比如：TOKASIDE模型的垂向飽和水力導(dǎo)水率Ksv和Green-Ampt下滲公式的飽和水力導(dǎo)水率Ks，從而減少模型參數(shù)的數(shù)量。TOKASIDE-D是基于物理基礎(chǔ)的水文模型，在進(jìn)行參數(shù)優(yōu)化率定時(shí)，大部分參數(shù)可根據(jù)土壤和土地利用性質(zhì)得來(lái)，而其余參數(shù)使用人工優(yōu)選法和SCE-UA自動(dòng)優(yōu)選法[25]結(jié)合率定得到。

2 研究流域與數(shù)據(jù)資料

2.1 研究流域

北辛店水文站為海河流域大清河水系清水河的控制站，控制流域面積1 650 km2(圖2)。該站實(shí)測(cè)最高洪水位為19.70 m，最大洪峰流量為710 m3/s，多年平均年徑流量為0.29億m3，多年平均年降水量為457.1 mm。該流域?qū)偬猩缴角皟A斜沖積平原區(qū)，地勢(shì)西高東低，緩緩傾斜，土壤以沙壤土為主。流域西部以山地丘陵為主，幾乎全為次生林或次生溫帶灌草叢，東部為平原區(qū)，地勢(shì)平坦開(kāi)闊，土質(zhì)疏松，土層深厚，已開(kāi)墾為農(nóng)田。

圖2 北辛店流域地形及站點(diǎn)分布Fig.2 Topography and station distribution of Beixindian Watershed

近年來(lái)，受人類(lèi)活動(dòng)影響，海河流域下墊面條件發(fā)生了較大變化，流域產(chǎn)匯流規(guī)律也發(fā)生一定程度的改變。受地下水過(guò)量開(kāi)采的影響，海河流域地下水埋深較深，不易蓄滿(mǎn)，產(chǎn)流機(jī)制呈現(xiàn)出“先超后蓄”的特點(diǎn)，超滲地面徑流產(chǎn)流機(jī)制發(fā)揮了重要作用。

2.2 數(shù)據(jù)資料

本文所需下墊面數(shù)據(jù)為數(shù)字高程地圖(DEM)、土壤類(lèi)型和土地利用類(lèi)型數(shù)據(jù)。DEM從地理空間數(shù)據(jù)云網(wǎng)站下載，選擇GDEMV2 DEM數(shù)字高程數(shù)據(jù)產(chǎn)品，分辨率為30 m。土壤數(shù)據(jù)來(lái)源于聯(lián)合國(guó)糧農(nóng)組織(FAO)，1∶100萬(wàn)土壤柵格數(shù)據(jù)資料。土地利用數(shù)據(jù)可從美國(guó)地質(zhì)勘探局(United States Geological Survey，USGS)下載，數(shù)據(jù)為分辨率為1 km的矢量文件。

研究流域水文數(shù)據(jù)包括降雨及流量數(shù)據(jù)(時(shí)段為1 h)，由河北省水文水資源勘測(cè)局提供。氣溫?cái)?shù)據(jù)來(lái)自中國(guó)國(guó)家氣象數(shù)據(jù)中心。流量數(shù)據(jù)來(lái)自北辛店水文站，降雨資料來(lái)自黃龍寺、順平等6個(gè)雨量站(圖2)。選用北辛店流域1980—2002年的12場(chǎng)次洪水?dāng)?shù)據(jù)，其中1980—1996年的9場(chǎng)為率定期，1997—2002年的3場(chǎng)為驗(yàn)證期。

3 TOKASIDE-D模型的應(yīng)用檢驗(yàn)

將TOKASIDE模型、Green-Ampt模型、TOKASIDE-D模型用于北辛店流域進(jìn)行應(yīng)用檢驗(yàn)。依據(jù)傳統(tǒng)的水文模擬與預(yù)報(bào)精度評(píng)定準(zhǔn)則，并參考《水文情報(bào)預(yù)報(bào)規(guī)范》[26]的規(guī)定，結(jié)合半干旱地區(qū)洪水特征，選擇以下評(píng)價(jià)指標(biāo)：徑流深誤差合格率，該誤差以實(shí)測(cè)徑流深的20%作為許可誤差判定預(yù)報(bào)徑流深是否合格；洪峰誤差合格率，以實(shí)測(cè)洪峰的20%作為許可誤差；峰現(xiàn)時(shí)差合格率，以峰現(xiàn)時(shí)間小于3 h為許可誤差；平均確定性系數(shù)，評(píng)價(jià)洪水實(shí)測(cè)過(guò)程與預(yù)報(bào)過(guò)程之間的擬合程度。北辛店流域3種模型參數(shù)率定的結(jié)果見(jiàn)表1，模擬結(jié)果見(jiàn)表2。

表1 北辛店流域各模型的優(yōu)化參數(shù)

表2 北辛店流域3個(gè)模型的模擬結(jié)果統(tǒng)計(jì)指標(biāo)

由表2可得：徑流深誤差合格率方面，率定期TOKASIDE-D模型最高(56%)，而Green-Ampt模型最低(33%)，驗(yàn)證期3個(gè)模型均只有1場(chǎng)合格，合格率33%；洪峰誤差合格率方面，率定期TOKASIDE-D模型也是最高(78%)，TOKASIDE模型最低(44%)，驗(yàn)證期TOKASIDE與Green-Ampt模型合格率均為33%，而TOKASIDE-D模型高于兩者。峰現(xiàn)時(shí)差方面，TOKASIDE模型和TOKASIDE-D模型的率定期和驗(yàn)證期合格率一致；但均低于Green-Ampt模型的合格率，說(shuō)明TOKASIDE-D模型在模擬洪峰時(shí)間方面與原TOKASIDE模型差異不大，土壤入滲后補(bǔ)充土壤水的過(guò)程延長(zhǎng)了飽和地面徑流的出現(xiàn)時(shí)間，進(jìn)而使得模擬洪峰時(shí)間出現(xiàn)滯后。TOKASIDE-D確定性系數(shù)比原模型有提高，說(shuō)明改進(jìn)模型對(duì)洪水形狀的模擬與觀(guān)測(cè)值更接近。

此外，選取2個(gè)典型的洪水事件1996080420和1980081700進(jìn)行詳細(xì)分析，其洪水實(shí)測(cè)過(guò)程線(xiàn)與模擬過(guò)程線(xiàn)如圖3所示。首先，從2個(gè)洪水事件的實(shí)測(cè)值可以看出，2個(gè)洪水的洪峰值差不多，前者洪水起漲快，退水慢；后者的洪水過(guò)程陡漲陡落，呈尖瘦形態(tài)。從實(shí)測(cè)過(guò)程線(xiàn)可以大致推斷出1996080420洪水的降水量大于1980081700，前者降水更多的是入滲到土壤中，而后者洪水明顯由超滲產(chǎn)流主導(dǎo)，降水入滲量少，基本上發(fā)生了超滲產(chǎn)流。

圖3 北辛店流域部分洪水事件實(shí)測(cè)過(guò)程線(xiàn)與模擬過(guò)程線(xiàn)對(duì)比Fig.3 Comparison of observed and simulated hydrographs of partial flood events in Beixindian Watershed

對(duì)比3個(gè)模型的模擬過(guò)程線(xiàn)，在1996080420洪水事件中(圖3(a))，TOKASIDE-D模型的模擬洪水過(guò)程線(xiàn)在起漲階段以及洪峰值效果最好，模擬洪水的整體形態(tài)與實(shí)測(cè)洪水最接近。TOKASIDE模型和Green-Ampt模型的模擬洪水起漲速度都很快，但模擬洪峰值都比實(shí)測(cè)洪峰大得多，且退水速度快，與實(shí)測(cè)退水過(guò)程偏差較大，說(shuō)明這2個(gè)模型在模擬時(shí)將大部分凈雨攔在了地表，轉(zhuǎn)為超滲地表徑流，入滲補(bǔ)充土壤含水量的部分偏少；而TOKASIDE-D模型在漲水階段的土壤蓄水量在退水時(shí)被“釋放”出來(lái)。但相對(duì)實(shí)測(cè)的退水流量還是偏小，這主要是下滲的水量在滲漏模塊的作用下繼續(xù)下滲補(bǔ)充地下水，從而出現(xiàn)部分水量不足的情況。在1980081700洪水事件中(圖3(b))，3個(gè)模型的模擬洪水過(guò)程線(xiàn)都呈現(xiàn)陡漲陡落的形態(tài)，與實(shí)際洪水過(guò)程線(xiàn)擬合較好。這是因?yàn)?980081700的洪水過(guò)程陡漲陡落，呈尖瘦形態(tài)，是比較典型的超滲產(chǎn)流洪水過(guò)程，觀(guān)察降雨過(guò)程，暴雨中心在流域中下游，因此得以快速匯流至流域出口，并沒(méi)有太多的下滲損失。另外，洪水發(fā)生在1980年，人類(lèi)活動(dòng)和水利工程相對(duì)較少，比較接近天然的降雨徑流過(guò)程，因此3種模型在該場(chǎng)次洪水模擬的效果相對(duì)較好。

與原TOKASIDE模型對(duì)所有網(wǎng)格設(shè)置超滲“門(mén)檻”不同的是，TOKASIDE-D模型通過(guò)對(duì)每時(shí)刻的土壤含水量和降雨等動(dòng)態(tài)因子進(jìn)行“監(jiān)控”，動(dòng)態(tài)調(diào)整網(wǎng)格的產(chǎn)流機(jī)制。以1996080420和1980081700這2個(gè)典型的洪水事件為例，展示洪水計(jì)算過(guò)程中產(chǎn)流機(jī)制的動(dòng)態(tài)變化過(guò)程，如圖4和圖5所示。從圖4可以看出，1996080420洪水從8月4日20:00開(kāi)始降水，一直到5日的2:00，降雨8 h，流域基本還是超滲狀態(tài)，在3:00流域出口附近網(wǎng)格基本蓄滿(mǎn)，在隨后的3 h內(nèi)，網(wǎng)格蓄滿(mǎn)數(shù)量大量增加，直到6日5:00達(dá)到最大，此時(shí)模擬洪水也達(dá)到了峰值，隨后，網(wǎng)格土壤含水量開(kāi)始下降，洪水過(guò)程處于退水階段，蓄滿(mǎn)網(wǎng)格數(shù)量逐漸減少。在這里也可以解釋?zhuān)琓OKASIDE-D模型在1996080420洪水的模擬中，網(wǎng)格蓄水程度高，退水時(shí)網(wǎng)格中有大量的壤中流和地下徑流產(chǎn)生，從而在退水階段模擬優(yōu)于其他2個(gè)模型。從1996080420洪水的實(shí)測(cè)過(guò)程可知其由超滲主導(dǎo)，這一點(diǎn)從圖5也可以看出，整個(gè)洪水過(guò)程蓄滿(mǎn)網(wǎng)格數(shù)量確實(shí)遠(yuǎn)遠(yuǎn)少于1996080420洪水，超滲網(wǎng)格占主導(dǎo)地位，進(jìn)而模擬的洪水過(guò)程線(xiàn)呈現(xiàn)尖瘦形態(tài)。

圖4 1996080420洪水的蓄滿(mǎn)網(wǎng)格與超滲網(wǎng)格動(dòng)態(tài)分布Fig.4 Dynamic distribution of saturation-excess and infiltration-excess grids for flood event 1996080420

圖5 1980081700洪水的蓄滿(mǎn)網(wǎng)格與超滲網(wǎng)格動(dòng)態(tài)分布Fig.5 Dynamic distribution of saturation-excess and infiltration-excess grids for flood event 1980081700

綜上可知，Green-Ampt模型在北辛店模擬效果相比其他2個(gè)模型較差，主要因?yàn)楫a(chǎn)流機(jī)制不符合當(dāng)?shù)氐奶卣鳌１毙恋炅饔蚪涤瓿跗谕寥浪康停菀装l(fā)生超滲產(chǎn)流，而后土壤入滲能力強(qiáng)，土壤蓄水量大，轉(zhuǎn)為蓄滿(mǎn)產(chǎn)流主導(dǎo)。由于Green-Ampt模型為超滲產(chǎn)流模型，凈雨大多發(fā)生產(chǎn)流，轉(zhuǎn)為超滲地表徑流，導(dǎo)致模擬洪水過(guò)程遠(yuǎn)大于實(shí)測(cè)過(guò)程。TOKASIDE模型雖然表現(xiàn)尚可，但是其模擬洪峰與實(shí)測(cè)洪峰相差較大；而TOKASIDE-D模型結(jié)合了蓄滿(mǎn)產(chǎn)流和超滲產(chǎn)流2種產(chǎn)流機(jī)制，并且根據(jù)土壤含水量以及降雨的時(shí)空分布，動(dòng)態(tài)調(diào)整產(chǎn)流機(jī)制，能更加準(zhǔn)確地把握北辛店流域的產(chǎn)流機(jī)制的過(guò)程，適應(yīng)性更強(qiáng)。該地區(qū)的滲漏嚴(yán)重，雖然增加了滲漏模塊，但部分場(chǎng)次洪水模擬結(jié)果依舊較差，在后續(xù)研究中，進(jìn)一步探究該流域的存在的問(wèn)題以及獨(dú)有特征，尋找合適的解決方法。

4 結(jié) 語(yǔ)

在基于物理基礎(chǔ)的分布式模型TOKASIDE模型的基礎(chǔ)上，通過(guò)降雨和土壤含水量之間的判定關(guān)系，動(dòng)態(tài)識(shí)別蓄滿(mǎn)產(chǎn)流和超滲產(chǎn)流網(wǎng)格，并添加滲漏模塊，構(gòu)建了基于物理基礎(chǔ)的分布式蓄超空間動(dòng)態(tài)組合TOKASIDE-D模型，并在北辛店流域進(jìn)行試驗(yàn)應(yīng)用。研究結(jié)果表明，與原始模型(TOKASIDE模型、Green-Ampt模型)相比，TOKASIDE-D蓄超空間動(dòng)態(tài)組合模型在北辛店流域?yàn)榇淼娜A北地區(qū)洪水模擬精度有較為明顯的提升。

TOKASIDE-D模型能夠很好地結(jié)合蓄滿(mǎn)模型和超滲模型的優(yōu)點(diǎn)，并且根據(jù)土壤含水量與降雨等因子的變化，動(dòng)態(tài)調(diào)整產(chǎn)流機(jī)制，并針對(duì)地下水超采考慮地下水的滲漏。隨著模型結(jié)構(gòu)與原理的改進(jìn)完善，TOKASIDE-D模型能夠準(zhǔn)確捕捉降雨、土壤含水量等動(dòng)態(tài)因子對(duì)產(chǎn)流模式的控制效應(yīng)，在我國(guó)北方地區(qū)的中小流域，有著很大潛力的研究?jī)r(jià)值以及較高的應(yīng)用價(jià)值。