基于蓄滿超滲時空動態(tài)組合的網(wǎng)格新安江模型

杜若愚，姚 成，劉玉環(huán)，李致家，張 珂，朱躍龍

(1.河海大學(xué)水文水資源學(xué)院，江蘇 南京 210098； 2.河海大學(xué)計算機與信息學(xué)院，江蘇 南京 210098)

產(chǎn)流模式主要分為蓄滿產(chǎn)流和超滲產(chǎn)流[1]，我國南方濕潤地區(qū)以蓄滿產(chǎn)流為主，北方干旱地區(qū)以超滲產(chǎn)流為主，半干旱半濕潤地區(qū)的產(chǎn)流模式兩者兼而有之[2-3]。半干旱半濕潤地區(qū)降雨時空分布不均，下墊面條件空間異質(zhì)性大，產(chǎn)流模式混合多變，即蓄滿產(chǎn)流和超滲產(chǎn)流隨時空變化明顯，導(dǎo)致多數(shù)水文模型洪水預(yù)報結(jié)果欠佳[4-5]。

半干旱半濕潤地區(qū)，是濕潤地區(qū)與干旱地區(qū)的過渡地帶，如果降雨落在具有很高下滲能力的土壤，就會像濕潤地區(qū)那樣發(fā)生壤中流和地下徑流，土壤蓄滿后發(fā)生飽和地面徑流，即蓄滿產(chǎn)流；而當(dāng)降雨落在下滲能力小，再加上土壤保水能力差，就會發(fā)生超滲產(chǎn)流，這就在空間尺度上形成了蓄滿產(chǎn)流與超滲產(chǎn)流同時發(fā)生的現(xiàn)象。而半干旱半濕潤地區(qū)洪水預(yù)報的難點在于確定產(chǎn)流模式，李致家等[6-8]認為產(chǎn)流模式是隨時空動態(tài)變化的，并在新安江模型的基礎(chǔ)上添加了超滲產(chǎn)流結(jié)構(gòu)，完善了產(chǎn)流理論，使其能夠適應(yīng)半干旱半濕潤地區(qū)復(fù)雜條件的混合產(chǎn)流洪水預(yù)報；李致家等[9]對基于子流域的蓄滿超滲空間組合進行研究，分別采用新安江模型和河北雨洪模型計算不同產(chǎn)流區(qū)的產(chǎn)流[10]，但這種產(chǎn)流計算為空間靜態(tài)組合，因此，Liu等[11-12]研究了基于網(wǎng)格的蓄滿超滲時空組合，考慮了產(chǎn)流模式隨時空變化的特點，基于TOKASIDE模型構(gòu)建TOKASIDE-D模型，并成功應(yīng)用于我國華北地區(qū)洪水預(yù)報。

TOKASIDE-D模型是基于物理基礎(chǔ)的分布式水文模型，其考慮了地形與運動波的蓄滿超滲產(chǎn)流、地下水運動、水庫入流、調(diào)蓄計算和滲漏等模塊，核心思路是產(chǎn)流模式隨時空的動態(tài)變化。由于網(wǎng)格新安江(Grid-XAJ)模型[13-18]自創(chuàng)建以來在全國多數(shù)流域得到了廣泛應(yīng)用，因此本文借鑒TOKASIDE-D模型中產(chǎn)流模式時空動態(tài)組合的思路，構(gòu)建基于蓄滿超滲時空動態(tài)組合的網(wǎng)格新安江(Grid-XAJ-SIDE)模型，并將Grid-XAJ、Grid-XAJ-SIDE和Grid-GA[19](網(wǎng)格格林-安普特)模型應(yīng)用于半干旱地區(qū)綏德流域進行驗證和分析，探討蓄滿、超滲產(chǎn)流模式時空動態(tài)組合方法對流域洪水模擬結(jié)果的影響。

1 模 型 構(gòu) 建

1.1 Grid-XAJ模型

Grid-XAJ模型以DEM網(wǎng)格為計算單元，主要分為蒸散發(fā)、蓄滿產(chǎn)流、分水源、匯流4個計算模塊[13-18]。模型假設(shè)每個網(wǎng)格的降雨和下墊面條件空間分布均勻，不存在張力水蓄水容量分布曲線和自由水蓄水容量分布曲線，只需考慮各要素在不同網(wǎng)格間的變異性。Grid-XAJ模型是蓄滿產(chǎn)流型分布式水文模型，廣泛應(yīng)用于我國濕潤地區(qū)洪水預(yù)報。

1.2 Grid-GA模型

Grid-GA模型[19]主要分為蒸散發(fā)模塊、產(chǎn)流模塊和匯流模塊3個計算模塊，各個模塊均在以DEM網(wǎng)格為基礎(chǔ)的正交網(wǎng)格內(nèi)進行。模型認為每個網(wǎng)格的降雨和下墊面條件空間分布均勻，不存在下滲能力分布曲線，Grid-GA模型是超滲產(chǎn)流型分布式水文模型，主要用于干旱半干旱地區(qū)洪水預(yù)報。

1.3 Grid-XAJ-SIDE模型構(gòu)建

Grid-XAJ-SIDE模型是在Grid-XAJ模型的基礎(chǔ)上，使用徑流曲線數(shù)-地形指數(shù)法劃分流域初始蓄超產(chǎn)流區(qū)分布，并添加超滲產(chǎn)流模塊，同時在產(chǎn)流計算過程中動態(tài)識別蓄滿網(wǎng)格和超滲網(wǎng)格，分別采用Grid-XAJ模型的產(chǎn)流與分水源和Green-Ampt下滲公式計算產(chǎn)流，最后根據(jù)網(wǎng)格間的匯流演算次序，依次將不同徑流成分演算至流域出口。

Grid-XAJ-SIDE模型主要分為蒸散發(fā)、蓄超網(wǎng)格動態(tài)識別、產(chǎn)流、分水源、匯流5個計算模塊，模型計算流程如圖1所示(圖中Pt為降雨量，Ipt為降雨強度，Wt為網(wǎng)格土壤含水量，WM為網(wǎng)格張力水蓄水容量，ft為下滲能力，t為計算時段，n為總時段數(shù))。

圖1 Grid-XAJ-SIDE模型計算流程Fig.1 Calculation flow-process of Grid-XAJ-SIDE model

Grid-XAJ-SIDE模型基于網(wǎng)格新安江模型假設(shè)，并添加以下假設(shè)條件：①網(wǎng)格內(nèi)降雨、土壤性質(zhì)、土地利用等特征均勻分布，不存在下滲能力分布曲線；②流域內(nèi)網(wǎng)格在計算過程中，依據(jù)土壤含水量是否達到田間持水量，降雨是否超過下滲能力，動態(tài)識別蓄滿網(wǎng)格或超滲網(wǎng)格；③每個網(wǎng)格的土壤垂向分布均勻。

1.3.1 蒸散發(fā)計算

Grid-XAJ-SIDE模型采用3層蒸散發(fā)模型計算流域?qū)嶋H蒸散發(fā)量。3層蒸散發(fā)模型將每個網(wǎng)格的土壤分為上層、下層和深層，對應(yīng)的土壤蓄水容量分別為WUM、WLM和WDM，具體計算公式詳見文獻[1]。

1.3.2 初始蓄超網(wǎng)格判定

Grid-XAJ-SIDE模型基于徑流曲線數(shù)和地形指數(shù)，對流域的蓄滿網(wǎng)格和超滲網(wǎng)格的初始空間分布進行確定，徑流曲線數(shù)源自SCS模型，表示某種土壤水分條件下的曲線數(shù)，主要用來描述地面超滲產(chǎn)流發(fā)生的難易程度，徑流曲線數(shù)越大，發(fā)生超滲徑流的可能性越大。初始蓄超網(wǎng)格判定采用的數(shù)據(jù)主要包含DEM、土壤類型和土地利用空間分布。劃分方法如下：計算每個網(wǎng)格內(nèi)的徑流曲線數(shù)和地形指數(shù)，再根據(jù)徑流曲線數(shù)進行初始分類[20]，當(dāng)徑流曲線數(shù)小于60時，將網(wǎng)格劃分為蓄滿網(wǎng)格，否則劃分為超滲網(wǎng)格。在此基礎(chǔ)上，再加入地形指數(shù)進一步修正，當(dāng)蓄滿網(wǎng)格的地形指數(shù)小于7，將該網(wǎng)格修改為超滲網(wǎng)格；當(dāng)網(wǎng)格的地形指數(shù)大于25，將該網(wǎng)格修改為蓄滿網(wǎng)格，最終得到整個流域的初始蓄超網(wǎng)格空間分布。

1.3.3 蓄超網(wǎng)格動態(tài)識別

與初始蓄超網(wǎng)格判定不同，蓄超網(wǎng)格動態(tài)識別是Grid-XAJ-SIDE模型在產(chǎn)流計算過程中，根據(jù)當(dāng)前時段網(wǎng)格單元的土壤含水量、降雨強度與土壤下滲能力之間的關(guān)系，動態(tài)識別蓄滿網(wǎng)格和超滲網(wǎng)格，是一種動態(tài)過程；而初始蓄超網(wǎng)格判定是在模型開始計算之前，基于徑流曲線數(shù)-地形指數(shù)法，根據(jù)下墊面信息進行判定，是一種靜態(tài)過程。

蓄超網(wǎng)格動態(tài)識別的主要原則：當(dāng)降雨強度大于下滲能力時，該網(wǎng)格為超滲網(wǎng)格；當(dāng)網(wǎng)格的土壤含水量達到田間持水量后，該網(wǎng)格則為蓄滿網(wǎng)格。其中，流域初始土壤含水量的值由日模型計算得到。

1.3.4 產(chǎn)流及分水源計算

1.3.4.1 蓄滿產(chǎn)流及分水源

Grid-XAJ-SIDE模型對于蓄滿網(wǎng)格采用蓄滿產(chǎn)流模式，即在降雨過程中，直到土壤含水量達到田間持水量才能產(chǎn)流，而在此之前，所有來水均用于補充土壤含水量。將計算時段內(nèi)網(wǎng)格單元的實測降雨先扣除相應(yīng)時段的蒸散發(fā)量，即可得到實際用于產(chǎn)流計算的時段雨量Pe，則

(1)

式中：R為時段產(chǎn)流量，mm；W0為上一時段土壤含水量，mm。

在Grid-XAJ-SIDE模型中，每個網(wǎng)格單元的R均被劃分為3種徑流成分即地表徑流Rs、壤中流Ri以及地下徑流Rg。與產(chǎn)流計算一樣，在進行分水源計算時，認為每個網(wǎng)格內(nèi)自由水蓄水容量分布均勻。分水源計算公式為

Ri=KiS

(2)

Rg=KgS

(3)

(4)

式中：Ki為自由水含量對壤中流的出流系數(shù)；Kg為自由水含量對地下水的出流系數(shù)；S為自由水含量，mm；SM為自由水蓄水容量，mm。

1.3.4.2 超滲產(chǎn)流

對于超滲網(wǎng)格，采用Green-Ampt下滲公式計算網(wǎng)格下滲能力：

(5)

式中：f(t)為土壤下滲能力，mm/h；Ks為飽和水力傳導(dǎo)度，mm/h；Ψ為濕潤鋒處的土壤吸力，mm；Δθ為飽和含水率與初始含水率之差；F(t)為累計下滲量，mm。計算時無需考慮下滲能力分布曲線，根據(jù)超滲產(chǎn)流原理，當(dāng)降雨強度大于土壤下滲能力時，以下滲能力下滲，超出下滲能力的部分降雨形成地表徑流，計算步驟如下：① 根據(jù)式(5)計算各網(wǎng)格下滲能力，其中，初始土壤含水率由日模型計算得到，初始土壤F(t=0)賦極小值0.000 01，即認為初始時刻的下滲能力無限大；②計算每個網(wǎng)格的R：

(6)

③計算該網(wǎng)格當(dāng)前時段的下滲量：i=Pe-R。該網(wǎng)格τ時段的累積下滲量為：F(t=τ)=F(t=τ-1)+i(t=τ)，其中i為當(dāng)前時段下滲量，mm。將F(t=τ)代入式(5)，計算下一時段的土壤下滲能力。每個網(wǎng)格存在一個最大累積下滲量Fmax，當(dāng)土壤累積下滲量大于最大累積下滲量時，土壤含水量達到田間持水量，多出的水量形成地下徑流。④重復(fù)步驟①，依次計算流域內(nèi)每個網(wǎng)格在各時段下的徑流量。

1.3.4.3 土壤含水量更新

Grid-XAJ-SIDE模型中，蓄滿網(wǎng)格土壤含水量按照蓄滿產(chǎn)流算法進行更新。超滲網(wǎng)格土壤含水量更新方式：根據(jù)超滲產(chǎn)流原理，當(dāng)降雨強度小于下滲能力時，降雨全部下滲到土壤中，地面不產(chǎn)生積水；當(dāng)降雨強度大于下滲能力時，降雨分為兩部分，一部分按照下滲能力下滲到土壤中，補充土壤含水量；另一部分則留在了地面，成為超滲地表徑流。設(shè)網(wǎng)格當(dāng)前時段下滲量為i，上一時段上層、下層和深層的土壤含水量分別為WU0、WL0和WD0，三層土壤的蓄水容量分別為WUM、WLM和WDM，按照先上層后下層的次序更新土壤含水量。入滲量優(yōu)先補充上層土壤含水量，若上層土壤達到蓄滿，多余的水量繼續(xù)下滲補充下層土壤含水量，如果下層土壤同樣已經(jīng)蓄滿，多余的水量繼續(xù)補充深層土壤含水量，直至整個土層土壤含水量達到田間持水量。

1.3.5 匯流計算

Grid-XAJ-SIDE模型匯流分為逐網(wǎng)格坡面匯流和河道匯流2個階段計算。坡面匯流采用一維擴散波方程計算，使用基于兩步法MacCormack算法[21]的二階顯式有限差分格式進行擴散波方程組的求解；河道匯流采用基于網(wǎng)格的馬斯京根河道匯流演算法[22]進行計算。

1.3.6 模型參數(shù)

Grid-XAJ-SIDE模型以Grid-XAJ模型為基礎(chǔ)，并添加了超滲產(chǎn)流計算模塊，實現(xiàn)了蓄超產(chǎn)流時空動態(tài)組合的定量模擬。由于下墊面的空間分布具有異質(zhì)性，研究流域的網(wǎng)格參數(shù)與下墊面的土壤、土地利用、高程等分布具有緊密相關(guān)性，本文涉及的3個模型均為分布式水文模型，其網(wǎng)格的上層張力蓄水容量WUM、下層張力蓄水容量WLM、張力水蓄水容量WM、自由水蓄水容量SM、地下徑流出流系數(shù)KG、壤中流出流系數(shù)KI、土壤總孔隙度θe、土壤有效孔隙度θs、濕潤鋒處土壤吸力ψ和飽和水力傳導(dǎo)度Ks等參數(shù)的取值也為分布式。其中WUM、WLM、WM、SM、KG和KI均根據(jù)Grid-XAJ參數(shù)及其空間分布估計方法進行計算[23-24]；每個網(wǎng)格的θe、θs、ψ和Ks等超滲產(chǎn)流參數(shù)，均與土壤類型相關(guān)，根據(jù)Rawls等[25-27]的研究確定每個網(wǎng)格的參數(shù)值；蒸發(fā)折算系數(shù)K、深層蒸散發(fā)系數(shù)C、最大累積下滲量IO、地下徑流消退系數(shù)CG、壤中流消退系數(shù)CI、河網(wǎng)蓄水消退系數(shù)CS、滯后時間Lag和河道徑流逐網(wǎng)格馬斯京根法演算參數(shù)Kech、Xech等參數(shù)則基于流域?qū)崪y水文資料，采用人工優(yōu)選法在經(jīng)驗取值范圍內(nèi)進行率定。

2 模 型 驗 證

2.1 研究流域

選擇半干旱地區(qū)的綏德流域作為研究流域，流域高程、水系及站點分布如圖2所示。綏德水文站位于陜西省綏德縣，東經(jīng)110°14′，北緯37°30′，建站于1959年6月，流域控制面積3 893 km2。該站多年平均降雨量為486 mm，多年平均氣溫為9.7 ℃，實測最大流量為2 350 m3/s，相應(yīng)水位為818.59 m。綏德流域地形東西跨度大，主要以山地為主，總體呈現(xiàn)西高東低的趨勢；土壤以壤質(zhì)砂土為主。

圖2 綏德流域地形、水系及站點分布Fig.2 Topography, river and station distribution of Suide Watershed

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

Grid-XAJ-SIDE模型所需下墊面數(shù)據(jù)主要包括DEM、土壤類型和土地利用類型。DEM數(shù)據(jù)選擇由美國太空總署(NASA)與國防部國家測繪局(NIMA)聯(lián)合測量的SRTM數(shù)據(jù)，分辨率為30 m，通過地理空間數(shù)據(jù)云(http://www.gscloud.cn)下載。土壤類型數(shù)據(jù)來源于聯(lián)合國糧農(nóng)組織(FAO)和維也納國際應(yīng)用系統(tǒng)研究所(IIASA)構(gòu)建的世界土壤數(shù)據(jù)庫HWSD (harmonized world soil database version)1∶100萬土壤柵格數(shù)據(jù)。土地利用類型數(shù)據(jù)通過美國地質(zhì)勘探局(United States Geological Survey，USGS)下載，數(shù)據(jù)是分辨率為1 km 的矢量文件。

以1 h為計算步長建立模型，選取2010—2018 年站點資料較為齊全且洪峰較大的15 場洪水資料進行計算，將收集到的實測降雨和流量資料插值為1h步長數(shù)據(jù)，其中前10場洪水資料用于模型參數(shù)率定，后5場洪水用于模型驗證。

2.3 驗證結(jié)果與分析

將Grid-XAJ模型、Grid-GA模型和Grid-XAJ-SIDE模型應(yīng)用于綏德流域進行應(yīng)用檢驗和對比分析時，網(wǎng)格大小均采用1 km×1 km，既可以更好地與土壤、植被等數(shù)據(jù)的網(wǎng)格大小保持一致，也可以在保證應(yīng)用精度的前提下，提高模型運算效率。Grid-GA模型和Grid-XAJ-SIDE模型對于超滲產(chǎn)流的計算，計算時段會影響模擬結(jié)果，但綏德流域原始資料大多為2 h步長，目前受限于流域內(nèi)觀測資料的精度，3個模型計算步長均選擇1 h。

依據(jù)GB/T 22482—2008《水文情報預(yù)報規(guī)范》[28]評價模型模擬結(jié)果，選擇徑流深誤差、洪峰相對誤差和確定性系數(shù)等洪水模擬指標(biāo)進行綜合評價。3個模型人工優(yōu)選參數(shù)見表1，率定和驗證的精度見表2。

表1 綏德流域各模型參數(shù)值

表2 綏德流域各模型模擬結(jié)果對比

對于徑流深的模擬，率定期Grid-XAJ-SIDE模型、Grid-XAJ模型和Grid-GA 模型模擬合格率分別為90.0%、80.0%和100.0%，驗證期合格率分別為80.0%、80.0%和100.0%，模擬結(jié)果均較好，主要由于綏德流域?qū)崪y徑流深普遍偏小，15場洪水中，只有20170723號洪水實測徑流深大于15 mm，且3個模型對于徑流深的模擬值與實測值誤差多數(shù)在3 mm以內(nèi)，依據(jù)GB/T 22482—2008《水文情報預(yù)報規(guī)范》[28]，3個模型的徑流深合格率都比較高。綏德流域地處半干旱區(qū)，洪水歷時短、漲落快、總徑流深小但洪峰較大，因此準確模擬洪峰更為重要。在洪峰合格率方面，Grid-XAJ-SIDE模型率定期和驗證期模擬結(jié)果均最好，分別為70.0%和60.0%，相比于Grid-XAJ模型(率定期60.0%，驗證期40.0%)和Grid-GA模型(率定期30.0%，驗證期60.0%)有明顯提高。3個模型對于峰現(xiàn)時間的模擬結(jié)果相同。3個模型的確定性系數(shù)普遍較差，說明在綏德這樣的半干旱流域，確定性系數(shù)不能用來描述陡漲陡落的洪水過程。

20130726號洪水事件中(圖3(a))，Grid-XAJ-SIDE模型和Grid-XAJ模型模擬結(jié)果接近實測過程，此次洪水的前期土壤蓄水程度較高，缺水量小，是以蓄滿產(chǎn)流為主導(dǎo)的洪水過程，初始土壤飽和度如圖4(a)所示。而Grid-GA模型模擬結(jié)果與實測過程相差較大，這是因為Grid-GA模型認為土壤是半無限土柱，即不會達到“蓄滿”狀態(tài)，入滲到土壤中的降雨不會再產(chǎn)生徑流。在20180710號洪水事件中(圖3(b))，Grid-XAJ-SIDE模型和Grid-GA模型模擬結(jié)果明顯優(yōu)于Grid-XAJ模型，這是因為此次洪水的前期土壤蓄水程度低，缺水量大，不易蓄滿，是以超滲產(chǎn)流為主導(dǎo)的洪水過程，初始土壤飽和度如圖4(b)所示。此次洪水過程中Grid-XAJ-SIDE模型和Grid-GA模型在洪水起漲階段模擬結(jié)果相近，且與實測過程接近，但在退水階段，Grid-GA模型退水過程過快，與實測過程偏差較大。Grid-XAJ-SIDE模型認為，即使是超滲產(chǎn)流主導(dǎo)的網(wǎng)格，當(dāng)下滲量補充土壤含水量達到“蓄滿”狀態(tài)后，下滲的水量會產(chǎn)生壤中流和地下徑流，因而比Grid-GA模型更接近實測退水過程。

圖3 綏德流域部分洪水實測過程線與模擬過程線對比Fig.3 Comparison of observed and simulated hydrographs of partial flood events in Suide Watershed

圖4 綏德流域部分洪水場次初始土壤飽和度分布Fig.4 Initial soil saturation distribution of partial flood events in Suide Watershed

Grid-XAJ-SIDE模型通過對逐時刻的土壤含水量和降雨等因子進行實時“監(jiān)控”，動態(tài)調(diào)整網(wǎng)格的產(chǎn)流方式。以20130726號洪水和20180710號洪水過程為例，展示在模型計算過程中蓄超網(wǎng)格的動態(tài)變化，如圖5所示。20130726號洪水過程的整體蓄水程度較高(圖4(a))，容易蓄滿，因此本場洪水從7月26日16:00開始降雨之后，蓄滿網(wǎng)格數(shù)量便迅速增多(圖5(a))，到7月27日2:00達到最大，此時模擬洪水達到頂峰，之后進入退水階段。從圖4(b)和圖5(b)可以看出，20180710號洪水過程整體蓄水程度不高，且下滲較小，由于降雨強度大，歷時短，隨著降雨的進行，從開始降雨(20180710T10:00)到洪峰時刻(20180711T14:00)蓄滿網(wǎng)格減少，整個洪水過程蓄滿網(wǎng)格數(shù)量少于超滲網(wǎng)格，是以超滲網(wǎng)格為主導(dǎo)的洪水過程，進而模擬的洪水過程也呈尖瘦形態(tài)。

圖5 綏德流域部分洪水蓄滿和超滲網(wǎng)格空間動態(tài)分布Fig.5 Spatial dynamic distribution of saturation-excess and infiltration-excess grids in Suide Watershed

Grid-XAJ-SIDE模型以網(wǎng)格為計算單元進行產(chǎn)匯流計算，并基于下墊面特性的分布信息推求模型參數(shù)的空間分布，因此模型在輸出流域出口斷面流量的同時，可以輸出任意網(wǎng)格單元的流量過程。以綏德流域20180710號洪水過程為例，由圖6可以看出，此次洪水事件前期主要發(fā)生在流域上游，隨著洪水過程的進行逐漸轉(zhuǎn)移至中下游，以洪水過程流量空間動態(tài)分布為依據(jù)，不僅可以反映流域各點的流量過程，還可以推測暴雨中心的轉(zhuǎn)移過程，也從側(cè)面反映出半干旱地區(qū)降雨空間分布的不均勻程度。

圖6 綏德流域20180710號洪水過程流量空間動態(tài)分布Fig.6 Spatial dynamic distribution of discharge for flood event 20180710 in Suide Watershed

通過3個模型的對比驗證與分析，Grid-XAJ-SIDE模型結(jié)構(gòu)完善，結(jié)合了蓄滿型產(chǎn)流和超滲產(chǎn)流模型的優(yōu)點，在下墊面條件比較復(fù)雜的半干旱地區(qū)提升了洪水模擬精度，對流域內(nèi)的洪水發(fā)生過程實現(xiàn)了精細描述。

3 結(jié) 語

在網(wǎng)格新安江模型(Grid-XAJ)的基礎(chǔ)上，基于徑流曲線數(shù)和地形指數(shù)劃分流域初始蓄超分布，構(gòu)建基于蓄滿超滲時空動態(tài)組合的網(wǎng)格新安江(Grid-XAJ-SIDE)模型，并以半干旱地區(qū)綏德流域?qū)崪y數(shù)據(jù)對模型進行了驗證。研究結(jié)果表明，與Grid-XAJ模型和Grid-GA模型相比，Grid-XAJ-SIDE模型根據(jù)土壤含水量和降雨等動態(tài)產(chǎn)流影響因子的實時變化，更加精確地描述了網(wǎng)格內(nèi)的產(chǎn)流模式，更好地結(jié)合了蓄滿產(chǎn)流模型和超滲產(chǎn)流模型的優(yōu)點，使洪水模擬精度有了較為明顯的提升，對半干旱半濕潤地區(qū)洪水預(yù)報具有積極意義。