王雪梅　郭良　翟曉燕

摘 要：山丘區(qū)小流域洪峰模擬結(jié)果與降雨時(shí)空分布密切相關(guān)。以河南省滎陽市王宗店流域“７·２０”山洪災(zāi)害為例，構(gòu)建了１５ 種時(shí)空分布的１０ ｍｉｎ 尺度降雨情景，綜合采用中國(guó)山洪水文模型和方差分析理論評(píng)估降雨時(shí)空情景變化對(duì)洪峰模擬的影響。結(jié)果表明：中國(guó)山洪水文模型可以較好地模擬“７·２０” 暴雨洪水過程。不同時(shí)間情景下１０ ｍｉｎ 降雨量均勻分布模擬得到的洪峰流量較小、峰現(xiàn)時(shí)間較晚；非均勻分布時(shí)最大１０ ｍｉｎ 降雨量發(fā)生時(shí)間越晚，模擬得到的洪峰流量越大、峰現(xiàn)時(shí)間越晚?？臻g情景變化主要影響洪峰流量模擬，采用泰森多邊形法模擬得到的洪峰流量最大。時(shí)間情景變化是洪峰流量和峰現(xiàn)時(shí)間模擬結(jié)果不確定性的最主要來源，空間情景變化和時(shí)空情景耦合變化僅影響洪峰流量模擬。

關(guān)鍵詞：暴雨山洪模擬；降雨時(shí)空不確定性；中國(guó)山洪水文模型；方差分析；王宗店流域

中圖分類號(hào)：Ｐ３３３．２ 文獻(xiàn)標(biāo)志碼：Ａ ｄｏｉ：１０．３９６９／ ｊ．ｉｓｓｎ．１０００－１３７９．２０２４．０４．００８

引用格式：王雪梅，郭良，翟曉燕．降雨時(shí)空不確定性對(duì)小流域洪峰模擬的影響評(píng)估［Ｊ］．人民黃河，２０２４，４６（４）：４９－５４．

水文模型利用概化的數(shù)學(xué)物理方程描述流域雨洪響應(yīng)過程，實(shí)現(xiàn)對(duì)暴雨洪水的模擬和預(yù)警預(yù)報(bào)［ １］ 。受模型輸入、模型結(jié)構(gòu)和模型參數(shù)等因素不確定性的影響，流域徑流模擬結(jié)果的不確定性客觀存在且不可忽視［２］ 。山丘區(qū)暴雨洪水時(shí)空異質(zhì)性強(qiáng)、產(chǎn)匯流過程非線性特征顯著［３］ ，暴雨山洪模擬的難度大、不確定性強(qiáng)［４］ 。厘清導(dǎo)致山洪模擬不確定性的關(guān)鍵要素，控制場(chǎng)次模擬結(jié)果的變化區(qū)間，是提高山洪災(zāi)害預(yù)警預(yù)報(bào)可靠度與準(zhǔn)確性的重要途徑。

模型輸入（降雨、流量等）的不確定性是水文模擬結(jié)果不確定性的主要原因［５］ 。國(guó)內(nèi)外已有不少學(xué)者針對(duì)降雨輸入的不確定性對(duì)洪水模擬效果的影響展開研究。陳華等［６］ 基于新安江模型和ＨＢＶ 模型研究了雨量站密度和空間分布對(duì)湘江流域徑流模擬的影響，發(fā)現(xiàn)增加雨量站密度、優(yōu)化站點(diǎn)空間分布可以提升場(chǎng)次徑流過程模擬精度。Ｓｉｌｖｅｓｔｒｏ 等［７］ 基于概率性洪水預(yù)報(bào)框架評(píng)估了降雨時(shí)空尺度對(duì)意大利Ｂｉｓａｇｎｏ Ｃｒｅｅｋ流域洪水模擬的影響，提出了研究區(qū)次洪模擬的最佳時(shí)空尺度配置方案。Ｌｉｕ 等［８］ 構(gòu)建了２０ 種流域降雨輸入方案以探究時(shí)空尺度對(duì)半干旱地區(qū)洪水模擬的影響，指出降雨時(shí)間尺度對(duì)次洪模擬的影響顯著。李大洋等［９］ 評(píng)估了ＭＩＫＥ ＳＨＥ 在淮河黃泥莊流域水文模擬的尺度適應(yīng)性，發(fā)現(xiàn)洪水模擬結(jié)果對(duì)網(wǎng)格尺寸與時(shí)間步長(zhǎng)的敏感性主要體現(xiàn)在洪峰和洪量上。已有研究多關(guān)注不同降雨時(shí)空特征的次洪模擬結(jié)果變化規(guī)律，各種不確定性對(duì)洪水模擬結(jié)果變化的貢獻(xiàn)有待進(jìn)一步量化。

相較于大流域雨洪模擬，山丘區(qū)小流域更側(cè)重于對(duì)洪水洪峰的模擬預(yù)報(bào)［１０］ 。本文以２０２１ 年７ 月２０日發(fā)生于河南省滎陽市王宗店流域的山洪災(zāi)害（“７·２０”山洪災(zāi)害）為例，基于累積降雨量—?dú)v時(shí)之間的回歸關(guān)系和設(shè)計(jì)條件下的雨峰特性對(duì)實(shí)測(cè)小時(shí)降雨量進(jìn)行時(shí)間上的降尺度插值，基于測(cè)站降雨量采用泰森多邊形法、反距離權(quán)重法和等雨量線法進(jìn)行空間上的流域降雨量插值，綜合利用中國(guó)山洪水文模型和方差分析理論定量評(píng)估降雨時(shí)空分布變化對(duì)山丘區(qū)洪水洪峰模擬的影響，以期為山丘區(qū)小流域洪水模擬的不確定性評(píng)估提供指導(dǎo)與幫助。

１ 研究區(qū)概況及數(shù)據(jù)源

王宗店流域位于河南省北部山丘區(qū)（ 東經(jīng)１１３°１７′—１１３°２２′ ，北緯３４°３６′ —３４°３９′ ），見圖１，集水面積２１．７ ｋｍ２。研究區(qū)屬亞熱帶濕潤(rùn)性季風(fēng)氣候區(qū)，多年平均氣溫１３．５ ℃，多年平均降水量７８５ ｍｍ。地勢(shì)南高北低，高程２４８ ～ ６５０ ｍ，最長(zhǎng)匯流路徑９．４ｋｍ，平均坡度１６．２°，平均河段比降２４．７‰。土壤類型為砂黏土（占８６％）和砂壤土（占１４％），土地利用類型為有林地（占７９％）和耕地（占２１％）。流域內(nèi)建有３２處小型塘堰路壩，總蓄水量超過１０ 萬ｍ３。王宗店村位于流域出口，村落現(xiàn)狀防洪能力不足１０ ａ 一遇。２０２１ 年７ 月２０ 日流域突發(fā)暴雨山洪，專家調(diào)查洪水洪峰流量為７６８ ｍ３ ／ ｓ，重現(xiàn)期為２００～３００ ａ［１１］ ，水流平均流速５ ｍ／ ｓ。洪水導(dǎo)致王宗店村２３ 人失蹤或遇難。

收集到的基礎(chǔ)數(shù)據(jù)主要有小流域數(shù)據(jù)、場(chǎng)次雨洪數(shù)據(jù)和《河南省中小流域設(shè)計(jì)暴雨洪水圖集》。小流域數(shù)據(jù)包括基于１ ∶ ５０ 萬ＤＥＭ、ＤＬＧ 和２．５ ｍ 遙感影像分析提取的小流域劃分和屬性信息。場(chǎng)次雨洪數(shù)據(jù)包括研究區(qū)周邊６ 個(gè)雨量站小時(shí)尺度的實(shí)測(cè)降雨數(shù)據(jù)（２０２１－０７－１８Ｔ８：００—０７－２２Ｔ１２：００）和王宗店水位站１５ ｍｉｎ 尺度的實(shí)測(cè)水位數(shù)據(jù)（２０２１－０７－１９Ｔ０：００—０７－２０Ｔ１７：００），均來源于河南省水文局，降雨集中時(shí)段（２０２１－０７－２０Ｔ６：００—１７：００）實(shí)測(cè)雨洪數(shù)據(jù)見表１，２０日５：００ 前各測(cè)站實(shí)測(cè)累積降雨量為１２９．６ ～ ３３１．０ｍｍ?！逗幽鲜≈行×饔蛟O(shè)計(jì)暴雨洪水圖集》包括河南省暴雨洪水設(shè)計(jì)成果，來源于河南省水文局。

２ 研究方法

２．１ 流域水文模擬

２．１．１ 中國(guó)山洪水文模型

中國(guó)山洪水文模型（Ｃｈｉｎａ Ｆｌａｓｈ Ｆｌｏｏｄ ｈｙｄｒｏｌｏｇｉｃａｌｍｏｄｅｌ， ＣＮＦＦ）是由中國(guó)水利水電科學(xué)研究院研發(fā)的分布式水文模型，已廣泛應(yīng)用于全國(guó)及各省（區(qū)）山洪模擬與預(yù)警預(yù)報(bào)研究［１２－１３］ 。本研究采用初損后損法進(jìn)行產(chǎn)流計(jì)算，采用分布式單位線法進(jìn)行坡面匯流計(jì)算，采用動(dòng)態(tài)馬斯京根法進(jìn)行河道洪水演算。

分布式單位線法應(yīng)用１ ∶ ５ 萬ＤＥＭ（數(shù)據(jù)網(wǎng)格為２５ ｍ×２５ ｍ）和２．５ ｍ 土地利用和植被信息等高精度地形地貌數(shù)據(jù)，通過式（１）計(jì)算流域內(nèi)各網(wǎng)格的坡面流流速和匯流時(shí)間，經(jīng)過一次線性水庫(kù)調(diào)蓄得到不同時(shí)段的小流域分布式單位線組（１０ ｍｉｎ、３０ ｍｉｎ、１ ｈ）。

式中：Ｖ 為水流速度；Ｋ 為坡面綜合流速系數(shù)，根據(jù)土地利用類型確定［１４］ ；Ｓ 為水力坡降；ｉ 為無因次雨強(qiáng)；Ｌｍ為第ｍ 個(gè)網(wǎng)格的水流流路長(zhǎng)度；Ｔｊ 為第ｊ 個(gè)網(wǎng)格的水流匯集時(shí)間；Ｖｍ 為第ｍ 個(gè)網(wǎng)格的流速；ｃ ＝ １ 或 ２；Ｍｊ為第ｊ 個(gè)網(wǎng)格水流流路上網(wǎng)格的數(shù)量。

２．１．２ 模型構(gòu)建及評(píng)估

基于流域ＤＥＭ 及山洪溝分布，研究區(qū)共劃分為１０ 個(gè)小流域單元（見圖１），各小流域單元及單元內(nèi)河道的基本屬性見表２。模型參數(shù)主要包括產(chǎn)流參數(shù)初損Ｉ０、穩(wěn)定下滲率ｆｃ，演進(jìn)參數(shù)糙率ｎ，馬斯京根法參數(shù)ｋ、ｘ。Ｉ０、ｆｃ和ｎ 基于場(chǎng)次前期降雨、流域水文下墊面條件等綜合確定，ｋ 和ｘ 采用概化的拋物線斷面和示儲(chǔ)流量動(dòng)態(tài)計(jì)算。

模型評(píng)估關(guān)注對(duì)“７·２０”暴雨洪水洪峰的模擬。基于實(shí)測(cè)水位數(shù)據(jù)和專家調(diào)查洪峰流量，采用洪峰流量相對(duì)誤差（ＥＱ ）和峰現(xiàn)時(shí)間誤差（ＴＰ ）對(duì)洪峰模擬效果進(jìn)行評(píng)估。指標(biāo)計(jì)算公式如下：

式中：Ｑｍｓ 和Ｑｍｅ 分別為模擬和專家調(diào)查洪峰流量，其中Ｑｍｅ ＝ ７６８ ｍ３ ／ ｓ；ＴＱｍｓ 和ＴＺｍｏ 分別為模擬流量過程和實(shí)測(cè)水位過程的峰現(xiàn)時(shí)間。

２．２ 降雨時(shí)空情景構(gòu)建

經(jīng)計(jì)算，“７·２０” 山洪場(chǎng)次匯流時(shí)間約為３１．３３ｍｉｎ。依據(jù)“雨洪模擬計(jì)算步長(zhǎng)不應(yīng)長(zhǎng)于場(chǎng)次匯流時(shí)間的１／３”的原則，選擇以１０ ｍｉｎ 為模型計(jì)算步長(zhǎng)并對(duì)測(cè)站小時(shí)降雨量插值形成１０ ｍｉｎ 降雨量系列。記最大１ ｈ 降雨（Ｍ６０－ｍａｘ）發(fā)生時(shí)段為主雨峰時(shí)段。對(duì)于非主雨峰時(shí)段，自降雨發(fā)生開始逐３ ｈ 滑動(dòng)構(gòu)建累積降雨量—?dú)v時(shí)之間的回歸關(guān)系進(jìn)而實(shí)現(xiàn)１０ ｍｉｎ 降雨量插值。對(duì)于主雨峰時(shí)段，認(rèn)為場(chǎng)次最大１０ ｍｉｎ 降雨量（Ｍ１０－ｍａｘ）發(fā)生在該時(shí)段且與Ｍ６０－ｍａｘ 頻率相同，通過式（４）對(duì)各測(cè)站Ｍ１０－ｍａｘ 進(jìn)行估算。假定Ｍ１０－ｍａｘ 分別發(fā)生在主雨峰時(shí)段的第２～５ 個(gè)１０ ｍｉｎ，結(jié)合主雨峰前后時(shí)段降雨構(gòu)建回歸關(guān)系進(jìn)行插值，形成４ 種主雨峰時(shí)段１０ ｍｉｎ 降雨量分布情景（記為Ｔ２～Ｔ５）。此外，還設(shè)置了小時(shí)降雨平均分配形成１０ ｍｉｎ 降雨量系列的時(shí)間情景（記為ＴＵ），以評(píng)估１０ ｍｉｎ 降雨量是否均勻分布對(duì)洪水模擬的影響。因此，共構(gòu)建了５ 種１０ ｍｉｎ 降雨量時(shí)間分布情景，即Ｔ２～Ｔ５、ＴＵ。

Ｍ１０－ｍａｘ ＝ ｋＭ６０－ｍａｘ （４）

式中：ｋ 為同頻率下Ｍ１０－ｍａｘ與Ｍ６０－ｍａｘ之比。

基于《河南省中小流域設(shè)計(jì)暴雨洪水圖集》計(jì)算ｋ在不同設(shè)計(jì)重現(xiàn)期（５ ～ ２００ ａ） 的取值，并構(gòu)建ｋ 與Ｍ６０－ｍａｘ的回歸關(guān)系（見圖２），即可基于測(cè)站實(shí)測(cè)Ｍ６０－ｍａｘ查得ｋ 相應(yīng)取值。

降雨的空間描述方案是影響分布式水文模型模擬精度的重要因素［１５］ ?；诓逯档玫降臏y(cè)站１０ ｍｉｎ 降雨量，采用常用的泰森多邊形法（Ｓ１）、反距離權(quán)重法（Ｓ２）和等雨量線法（Ｓ３）計(jì)算流域平均降雨量，形成３種降雨量空間分布?？紤]降雨時(shí)空分布的影響，最終形成１５ 種降雨情景進(jìn)行洪水模擬影響評(píng)估。

２．３ 不確定性來源分解

采用方差分析（Ａｎａｌｙｓｉｓ ｏｆ Ｖａｒｉａｎｃｅ，ＡＮＯＶＡ） 解析各降雨不確定性來源對(duì)洪峰模擬的相對(duì)影響［１６］ 。ＡＮＯＶＡ 將不同降雨情景驅(qū)動(dòng)下洪峰要素（洪峰流量、峰現(xiàn)時(shí)間）模擬結(jié)果的不確定性（ＵＦ）分解為源自降雨時(shí)間和空間分布（即主效應(yīng)）及時(shí)空分布耦合作用（即交互項(xiàng)）的不確定性，即

式中：Ｙｉ，ｊ為使用第ｉ 種降雨時(shí)間情景和第ｊ 種降雨空間情景時(shí)洪峰要素的模擬結(jié)果；Ｙ －為所有降雨情景下洪峰要素模擬結(jié)果的均值；Ｙ －ｉ 和Ｙ －ｊ 分別為使用第ｉ 種降雨時(shí)間情景和第ｊ 種降雨空間情景時(shí)洪峰要素模擬結(jié)果的均值；Ｉ 和Ｊ 分別為降雨時(shí)間和空間情景數(shù)，本研究中Ｉ ＝５，Ｊ ＝３。

３ 結(jié)果與分析

３．１ 降雨情景特征分析

由圖２ 可知，ｋ 隨最大１ ｈ 降雨量增大而減小?；跍y(cè)站實(shí)測(cè)降雨量確定各雨量站的ｋ 值為０．３５３（寺溝） ～０．３８１（盆窯），計(jì)算得到雨量站最大１０ ｍｉｎ 降雨量為２０．６（盆窯） ～３３．０ ｍｍ（寺溝）。

鑒于場(chǎng)次匯流時(shí)間不足１ ｈ，對(duì)１５ 種降雨情景下流域降雨特征的分析重點(diǎn)是峰前１ ｈ 時(shí)段，即２０ 日１３：００—１４：００。時(shí)段降雨總量主要受降雨空間情景影響，Ｓ１～ Ｓ３ 情景下時(shí)段雨量分別為９０．６ ～ ９１．４ ｍｍ、８４．６～８５．５ ｍｍ 和８８．０～８８．３ ｍｍ?？偟膩砜?，３ 種空間情景下時(shí)段雨量相差不大，其中Ｓ１ 情景降雨量略大于Ｓ２ 和Ｓ３ 的，其原因是泰森多邊形法插值結(jié)果易受距離較近測(cè)站影響，柿樹灣雨量站距流域最近且主雨峰較大，導(dǎo)致Ｓ１ 情景下時(shí)段流域降雨總量略大。時(shí)段降雨時(shí)程分布主要受降雨時(shí)間情景影響，以采用空間情景Ｓ１ 為例，１３：００—１４：００ 內(nèi)５ 種時(shí)間情景下流域１０ｍｉｎ 降雨量的累積分布見圖３。最大１０ ｍｉｎ 降雨量的量級(jí)和發(fā)生時(shí)段不同，導(dǎo)致不同時(shí)間情景下降雨的時(shí)程分布存在明顯差異。

３．２ 場(chǎng)次模擬效果評(píng)估

不同降雨情景下洪水模擬結(jié)果見表３，“７·２０”山洪模擬過程見圖４。１５ 種降雨情景驅(qū)動(dòng)下ＥＱ 為－３６．９７％～－１４．９６％，ＴＰ為－０．３３～０．１７ ｈ，ＣＮＦＦ 較好地再現(xiàn)了場(chǎng)次洪峰流量和峰現(xiàn)時(shí)間。當(dāng)前模型未考慮流域內(nèi)塘堰路壩潰決的影響，導(dǎo)致模擬得到的洪峰流量略低于專家調(diào)查洪峰流量。

３．３ 降雨時(shí)空情景對(duì)洪峰模擬的影響

“７·２０”暴雨洪水洪峰流量和峰現(xiàn)時(shí)間模擬結(jié)果見表３。１５ 種情景下模擬洪峰流量Ｑｍｓ為４８４．１～６５３．１ ｍ３ ／ ｓ，模擬峰現(xiàn)時(shí)間ＴＱｍｓ 為２０ 日１３：４０～１４：１０。對(duì)于洪峰流量而言，不同時(shí)間情景下１０ ｍｉｎ 降雨量非均勻分布［３ 種空間情景下Ｔ２～Ｔ５ 與ＴＵ 相比Ｑｍｓ分別增加５１．７～１３３．２ ｍ３ ／ ｓ（Ｓ１）、５０．４～１２３．１ ｍ３ ／ ｓ（Ｓ２）和５０．２～１３１．０ ｍ３ ／ ｓ（Ｓ３）］和最大１０ ｍｉｎ 降雨量發(fā)生位置偏后［３ 種空間情景下Ｔ５ 與其余情景相比Ｑｍｓ 分別增加３６．８～１３３．２ ｍ３ ／ ｓ（Ｓ１）、２８．６～１２３．１ ｍ３ ／ ｓ（Ｓ２）和３４．５～１３１．０ ｍ３ ／ ｓ（Ｓ３）］均會(huì)形成較大的洪峰流量，體現(xiàn)了短歷時(shí)峰值降雨的量級(jí)和發(fā)生時(shí)間對(duì)山丘區(qū)小流域洪峰流量模擬的關(guān)鍵性。不同空間情景下Ｑｍｓ 分布規(guī)律與峰前１ ｈ 時(shí)段雨量分布規(guī)律一致，Ｓ１ 情景下Ｑｍｓ較大，其次為Ｓ３、Ｓ２，Ｓ１ 較Ｓ２、Ｓ３ 的Ｑｍｓ分別增加３５．８～４５．９ｍ３ ／ ｓ、１７．６～２４．３ ｍ３ ／ ｓ。

對(duì)于峰現(xiàn)時(shí)間而言，不同時(shí)間情景下１０ ｍｉｎ 降雨量均勻分布和最大１０ ｍｉｎ 降雨量發(fā)生位置偏后均會(huì)導(dǎo)致峰現(xiàn)時(shí)間較晚，５ 種時(shí)間情景下ＴＱｍｓ 分別為１３：４０、１３：５０、１４：００、１４：１０ 和１４：１０?？臻g情景變化不會(huì)引起ＴＱｍｓ 的改變，Ｓ１ ～ Ｓ３ 情景下ＴＱｍｓ 集合均為１３：４０～１４：１０。

３．４ 洪峰模擬不確定性分解

洪峰流量和峰現(xiàn)時(shí)間模擬結(jié)果的不確定來源相對(duì)貢獻(xiàn)率見表４。對(duì)于Ｑｍｓ而言，降雨時(shí)間情景變化是不確定性的最主要來源，相對(duì)貢獻(xiàn)率為８７．３２％；其次為空間情景變化，相對(duì)貢獻(xiàn)率為１２．２８％；時(shí)空情景耦合變化對(duì)模擬不確定性的影響較小，相對(duì)貢獻(xiàn)率僅為０．４０％。對(duì)于ＴＱｍｓ 而言，時(shí)間情景變化是模擬結(jié)果不確定性的唯一來源，空間情景變化不會(huì)影響峰現(xiàn)時(shí)間模擬結(jié)果。

相較于空間情景變化，降雨時(shí)間情景變化是導(dǎo)致王宗店流域“７·２０”事件洪峰模擬不確定性的關(guān)鍵。暴雨山洪產(chǎn)匯流時(shí)間短、形成發(fā)展速度快，場(chǎng)次洪峰的形成發(fā)展與短歷時(shí)降雨特性密切相關(guān)。隨著最大１０ｍｉｎ 降雨量級(jí)和發(fā)生時(shí)間的改變，流域模擬洪峰量級(jí)和峰現(xiàn)時(shí)間均隨之顯著變化，體現(xiàn)了短歷時(shí)雨峰特性對(duì)山丘區(qū)小流域山洪響應(yīng)模擬的顯著影響。王宗店流域面積較小，降雨空間分布的異質(zhì)性較弱［１７］ ，加之流域周邊雨量測(cè)站的數(shù)量及分布均具有較好的代表性，采用不同的空間插值方法所得到的流域降雨在空間分布和平均量級(jí)等方面差別不大，因此不同降雨空間情景下場(chǎng)次洪峰模擬結(jié)果基本一致。

４ 結(jié)論

本文以河南省滎陽市王宗店流域“７·２０”山洪災(zāi)害為例，基于小時(shí)尺度實(shí)測(cè)降雨，綜合利用累積降雨量與歷時(shí)之間的回歸關(guān)系、設(shè)計(jì)條件下的雨峰特性和降雨空間插值方法，考慮降雨時(shí)空分布變化構(gòu)建了１５ 種１０ ｍｉｎ 尺度下的降雨輸入情景，利用中國(guó)山洪水文模型和方差分析理論解析并量化了降雨時(shí)空情景變化對(duì)洪峰模擬不確定性的影響，主要結(jié)論如下：

１）中國(guó)山洪水文模型合理地再現(xiàn)了王宗店流域“７·２０”山洪流域洪峰響應(yīng)過程，１５ 種降雨情景下洪峰流量相對(duì)誤差為－３６．９７％～－１４．９６％，峰現(xiàn)時(shí)間誤差為－０．５０～０．１７ ｈ，因未考慮流域內(nèi)小型塘堰路壩潰決而導(dǎo)致模擬洪峰流量偏低。

２）不同降雨時(shí)間情景下，１０ ｍｉｎ 降雨量均勻分布模擬得到的洪峰流量較小、峰現(xiàn)時(shí)間較晚。對(duì)于非均勻分布的１０ ｍｉｎ 降雨量情景，最大１０ ｍｉｎ 降雨量發(fā)生位置越晚，模擬得到的洪峰流量越大、峰現(xiàn)時(shí)間越晚。不同降雨空間情景下，采用泰森多邊形法插值流域降雨量得到的模擬洪峰流量最大，而模擬峰現(xiàn)時(shí)間受空間情景影響較小。

３）降雨時(shí)間情景變化是洪峰流量和峰現(xiàn)時(shí)間模擬的主要不確定性來源，相對(duì)貢獻(xiàn)率分別為８７．３２％和１００．００％，空間情景變化和時(shí)空情景耦合變化僅影響洪峰流量模擬，相對(duì)貢獻(xiàn)率分別為１２．２８％和０．４０％。研究區(qū)流域面積較小且雨量站空間數(shù)量代表性較好，采用不同空間插值方法得到的流域降雨無顯著差異，因此降雨空間情景變化對(duì)洪峰模擬不確定性影響較小。

今后應(yīng)進(jìn)一步收集水文氣象條件、地形、下墊面條件等差異顯著的山丘區(qū)小流域場(chǎng)次雨洪資料，以探究降雨時(shí)空變化對(duì)次洪模擬不確定性影響的普適規(guī)律。

參考文獻(xiàn)：

［１］ ＷＥＳＴＲＡ Ｓ，ＴＨＹＥＲ Ｍ，ＬＥＯＮＡＲＤ Ｍ，ｅｔ ａｌ．Ａ Ｓｔｒａｔｅｇｙ ｆｏｒ Ｄｉ?ａｇｎｏｓｉｎｇ ａｎｄ Ｉｎｔｅｒｐｒｅｔｉｎｇ Ｈｙｄｒｏｌｏｇｉｃａｌ Ｍｏｄｅｌ Ｎｏｎｓｔａｔｉｏｎａｒｉｔｙ［Ｊ］．Ｗａｔｅｒ Ｒｅｓｏｕｒｃｅｓ Ｒｅｓｅａｒｃｈ，２０１４，５０（６）：５０９０－５１１３．

［２］ 楊小柳，王藝臻．水文模型的不確定性和異因同果性：認(rèn)知的回顧［Ｊ］．水利學(xué)報(bào)，２０１８，４９（１）：１３６－１４３．

［３］ ＺＨＡＩ Ｘ Ｙ，ＧＵＯ Ｌ，ＺＨＡＮＧ Ｙ Ｙ．Ｆｌａｓｈ Ｆｌｏｏｄ Ｔｙｐｅ Ｉｄｅｎｔｉｆｉ?ｃａｔｉｏｎ ａｎｄ Ｓｉｍｕｌａｔｉｏｎ Ｂａｓｅｄ ｏｎ Ｆｌａｓｈ Ｆｌｏｏｄ Ｂｅｈａｖｉｏｒ Ｉｎｄｉｃｅｓｉｎ Ｃｈｉｎａ［Ｊ］．Ｓｃｉｅｎｃｅ Ｃｈｉｎａ Ｅａｒｔｈ Ｓｃｉｅｎｃｅｓ，２０２１，６４（７）：１１４０－１１５４．

［４］ ＢＯＲＧＡ Ｍ，ＡＮＡＧＮＯＳＴＯＵ Ｅ Ｎ，ＢＬ?ＳＣＨＬ Ｇ，ｅｔ ａｌ．ＦｌａｓｈＦｌｏｏｄ Ｆｏｒｅｃａｓｔｉｎｇ， Ｗａｒｎｉｎｇ ａｎｄ Ｒｉｓｋ Ｍａｎａｇｅｍｅｎｔ： ｔｈｅ ＨＹ?ＤＲＡＴＥ Ｐｒｏｊｅｃｔ［Ｊ］．Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｓｃｉｅｎｃｅ ＆ Ｐｏｌｉｃｙ，２０１１，１４（７）：８３４－８４４．

［５］ ＺＨＯＵ Ｘ Ｄ，ＭＡ Ｗ Ｃ，ＥＣＨＩＺＥＮＹＡ Ｗ，ｅｔ ａｌ．Ｔｈｅ Ｕｎｃｅｒｔａｉｎｔｙｏｆ Ｆｌｏｏｄ Ｆｒｅｑｕｅｎｃｙ Ａｎａｌｙｓｅｓ ｉｎ Ｈｙｄｒｏｄｙｎａｍｉｃ Ｍｏｄｅｌ Ｓｉｍｕｌａ?ｔｉｏｎｓ［Ｊ］．Ｎａｔｕｒａｌ Ｈａｚａｒｄｓ ａｎｄ Ｅａｒｔｈ Ｓｙｓｔｅｍ Ｓｃｉｅｎｃｅｓ，２０２１，２１（３）：１０７１－１０８５．

［６］ 陳華，霍苒，曾強(qiáng)，等．雨量站網(wǎng)布設(shè)對(duì)水文模型不確定性影響的比較［Ｊ］．水科學(xué)進(jìn)展，２０１９，３０（１）：３４－４４．

［７］ ＳＩＬＶＥＳＴＲＯ Ｆ，ＲＥＢＯＲＡ Ｎ，ＧＩＡＮＮＯＮＩ Ｆ，ｅｔ ａｌ．Ｔｈｅ ＦｌａｓｈＦｌｏｏｄ ｏｆ ｔｈｅ Ｂｉｓａｇｎｏ Ｃｒｅｅｋ ｏｎ ９ｔｈ Ｏｃｔｏｂｅｒ ２０１４：Ａｎ “Ｕｎｆｏｒ?ｔｕｎａｔｅ” Ｃｏｍｂｉｎａｔｉｏｎ ｏｆ Ｓｐａｔｉａｌ ａｎｄ Ｔｅｍｐｏｒａｌ Ｓｃａｌｅｓ ［ Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ｈｙｄｒｏｌｏｇｙ，２０１６，５４１：５０－６２．

［８］ ＬＩＵ Ｙ Ｈ，ＬＩ Ｚ Ｊ，ＬＩＵ Ｚ Ｙ，ｅｔ ａｌ．Ｉｍｐａｃｔ ｏｆ Ｒａｉｎｆａｌｌ Ｓｐａｔｉｏ?ｔｅｍｐｏｒａｌ Ｖａｒｉａｂｉｌｉｔｙ ａｎｄ Ｍｏｄｅｌ Ｓｔｒｕｃｔｕｒｅｓ ｏｎ Ｆｌｏｏｄ Ｓｉｍｕｌａｔｉｏｎｉｎ Ｓｅｍｉ?Ａｒｉｄ Ｒｅｇｉｏｎｓ［Ｊ］．Ｓｔｏｃｈａｓｔｉｃ Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｒｅｓｅａｒｃｈａｎｄ Ｒｉｓｋ Ａｓｓｅｓｓｍｅｎｔ，２０２２，３６（３）：７８５－８０９．

［９］ 李大洋，梁忠民，周艷．基于ＭＩＫＥ ＳＨＥ 的洪水模擬與尺度效應(yīng)分析［Ｊ］．水力發(fā)電，２０１９，４５（５）：２８－３３，６６．

［１０］ 肖楠，葉磊，吳劍，等．降雨對(duì)山丘區(qū)小流域洪峰模擬不確定性的影響［Ｊ］．中國(guó)農(nóng)村水利水電，２０１８（７）：３５－３８，４３．

［１１］ 吳澤斌，萬海斌．２０２１ 年河南鄭州山區(qū)４ 市“７·２０”特大暴雨災(zāi)害簡(jiǎn)析［Ｊ］．中國(guó)防汛抗旱，２０２２，３２（３）：２７－３１．

［１２］ 翟曉燕，郭良，劉榮華，等．中國(guó)山洪水文模型研制與應(yīng)用：以安徽省中小流域?yàn)槔郏剩荩畱?yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào)，２０２０，２８（５）：１０１８－１０３６．

［１３］ 王雪梅，郭良，翟曉燕，等． 基于中國(guó)山洪水文模型（ＣＮＦＦ）的中小流域洪水頻率分析［Ｊ］．中國(guó)防汛抗旱，２０２２，３２（６）：３１－３５．

［１４］ ＵＳＤＡ－ＳＣＳ．Ｎａｔｉｏｎａｌ Ｅｎｇｉｎｅｅｒｉｎｇ Ｈａｎｄｂｏｏｋ，Ｓｅｃｔｉｏｎ ４，Ｈｙ?ｄｒｏｌｏｇｙ［Ｍ］．Ｗａｓｈｉｎｇｔｏｎ Ｄ Ｃ：Ｓｏｉｌ Ｃｏｎｓｅｒｖａｔｉｏｎ Ｓｅｒｖｉｃｅｓ，ＵＳ Ｄｅｐａｒｔｍｅｎｔ ｏｆ Ａｇｒｉｃｕｌｔｕｒｅ，１９８５：２５－２８．

［１５］ 張雪松，郝芳華，張建永．降雨空間分布不均勻性對(duì)流域徑流和泥沙模擬影響研究［Ｊ］．水土保持研究，２００４，１１（１）：９－１２．

［１６］ ＶＥＴＴＥＲ Ｔ，ＨＵＡＮＧ Ｓ，ＡＩＣＨ Ｖ，ｅｔ ａｌ．Ｍｕｌｔｉ?ｍｏｄｅｌ ＣｌｉｍａｔｅＩｍｐａｃｔ Ａｓｓｅｓｓｍｅｎｔ ａｎｄ Ｉｎｔｅｒｃｏｍｐａｒｉｓｏｎ ｆｏｒ Ｔｈｒｅｅ Ｌａｒｇｅ?Ｓｃａｌｅ Ｒｉｖｅｒ Ｂａｓｉｎｓ ｏｎ Ｔｈｒｅｅ Ｃｏｎｔｉｎｅｎｔｓ［Ｊ］．Ｅａｒｔｈ ＳｙｓｔｅｍＤｙｎａｍｉｃｓ，２０１５，６（１）：１７－４３．

［１７］ ＨＵＡＮＧ Ｙ Ｃ，Ｂ?ＲＤＯＳＳＹ Ａ，ＺＨＡＮＧ Ｋ．Ｓｅｎｓｉｔｉｖｉｔｙ ｏｆ Ｈｙｄｒｏ?ｌｏｇｉｃａｌ Ｍｏｄｅｌｓ ｔｏ Ｔｅｍｐｏｒａｌ ａｎｄ Ｓｐａｔｉａｌ Ｒｅｓｏｌｕｔｉｏｎｓ ｏｆ ＲａｉｎｆａｌｌＤａｔａ［Ｊ］．Ｈｙｄｒｏｌｏｇｙ ａｎｄ Ｅａｒｔｈ Ｓｙｓｔｅｍ Ｓｃｉｅｎｃｅｓ，２０１９，２３（６）：２６４７－２６６３．

【責(zé)任編輯 許立新】

基金項(xiàng)目：國(guó)家自然科學(xué)基金資助項(xiàng)目（４２１７１０４７）