缺資料地區(qū)小型水庫(kù)納雨能力計(jì)算模型構(gòu)建及應(yīng)用

摘要：針對(duì)缺資料地區(qū)小型水庫(kù)預(yù)警難題，構(gòu)建了考慮不同小水庫(kù)運(yùn)行方式的納雨能力計(jì)算模型，并提出適用于缺資料地區(qū)的模型參數(shù)方案，開(kāi)發(fā)考慮未來(lái)降雨的納雨能力迭代試算方法。以2023年8月21日為起始時(shí)刻，計(jì)算湖南省缺資料地區(qū)32座試點(diǎn)小型水庫(kù)的納雨能力，并以和平水庫(kù)為例開(kāi)展納雨能力滾動(dòng)計(jì)算。結(jié)果表明：計(jì)算得到的試點(diǎn)小型水庫(kù)防洪風(fēng)險(xiǎn)與降雨分布合理相關(guān)，能夠直觀展現(xiàn)降雨主落區(qū)內(nèi)的高風(fēng)險(xiǎn)水庫(kù)分布，實(shí)現(xiàn)了庫(kù)水位和前期降雨對(duì)于納雨能力影響的合理量化，并且能夠?qū)崟r(shí)反映不同時(shí)期影響納雨能力的主要因素的動(dòng)態(tài)變化。

關(guān) 鍵 詞：小型水庫(kù)；納雨能力；前期影響雨量；迭代試算；缺資料地區(qū)；湖南省

中圖法分類(lèi)號(hào)：P339

文獻(xiàn)標(biāo)志碼：ADOI：10.16232/j.cnki.1001-4179.2024.08.013

0 引 言

小水庫(kù)量大面廣、建設(shè)標(biāo)準(zhǔn)偏低、運(yùn)行維護(hù)水平較差，成為水庫(kù)安全度汛工作的關(guān)鍵要害。根據(jù)《2023中國(guó)水利發(fā)展報(bào)告》^［1］統(tǒng)計(jì)，截至2022年底，中國(guó)共有小型水庫(kù)93 390座，小水庫(kù)占全部水庫(kù)數(shù)量的94.7%。納雨能力^［2］指考慮水庫(kù)工程參數(shù)、泄水設(shè)施運(yùn)行調(diào)度和產(chǎn)匯流特性等因素，在防汛安全情況下預(yù)見(jiàn)期內(nèi)水庫(kù)允許容納的最大降雨量，分析小型水庫(kù)的納雨能力對(duì)水庫(kù)預(yù)報(bào)預(yù)警及安全度汛具有重要意義。

已有學(xué)者開(kāi)展了小型水庫(kù)納雨能力計(jì)算相關(guān)方法研究，如國(guó)洪琴^［3］對(duì)流域初始土壤含水狀態(tài)作一定的人為預(yù)設(shè)，研究得到了某種理想狀態(tài)下的納雨能力。侯?lèi)?ài)中等^［4］考慮流域當(dāng)前下墊面情況、水庫(kù)當(dāng)前水位、水庫(kù)調(diào)度方式等因素，反查曲線求得水庫(kù)納雨能力。章四龍等^［5］集成定量降雨預(yù)報(bào)、水庫(kù)預(yù)報(bào)調(diào)度及納雨能力計(jì)算，較為全面地論述了小水庫(kù)納雨能力預(yù)報(bào)方法。然而，針對(duì)缺資料地區(qū)的小型水庫(kù)批量納雨能力計(jì)算，鮮有較為系統(tǒng)深入的研究^［6］。因此，本文構(gòu)建適用于不同小水庫(kù)運(yùn)行方式的納雨能力計(jì)算模型，綜合量化降雨以及流域下墊面特征對(duì)于納雨能力模型參數(shù)的影響，開(kāi)發(fā)考慮未來(lái)降雨的納雨能力迭代試算方法，以湖南省為例開(kāi)展缺資料地區(qū)32座試點(diǎn)小型水庫(kù)的納雨能力快速批量計(jì)算，驗(yàn)證納雨能力模型的合理性和適用性。

1 納雨能力計(jì)算模型構(gòu)建

1.1 適用于不同水庫(kù)運(yùn)行階段的納雨能力計(jì)算模型構(gòu)建

水庫(kù)可納雨庫(kù)容指水庫(kù)當(dāng)前水位距某一特征水位之間的庫(kù)容^［7］，本次研究首先計(jì)算當(dāng)前水位至汛限水位或溢洪道底坎高程之間的納雨能力，然后根據(jù)水庫(kù)當(dāng)前水位以及預(yù)見(jiàn)期內(nèi)的降雨判斷水庫(kù)是否考慮泄流，進(jìn)而計(jì)算水位到達(dá)至設(shè)計(jì)洪水位的納雨能力，并考慮降雨對(duì)徑流系數(shù)的影響，實(shí)現(xiàn)適用于不同水庫(kù)運(yùn)行階段的納雨能力計(jì)算。

1.1.1 不考慮泄流的納雨能力計(jì)算方法

當(dāng)庫(kù)水位不超過(guò)汛限水位或溢洪道底坎高程時(shí)，采用不考慮泄流的方法，計(jì)算當(dāng)前庫(kù)水位相對(duì)于汛限水位或溢洪道底坎高程所能容納的空余庫(kù)容，將空余庫(kù)容折算成降雨量，實(shí)現(xiàn)不考慮泄流的納雨能力計(jì)算，如公式（1）所示。

式中：P_c為不超過(guò)汛限水位或者溢洪道底坎高程的納雨能力；W_z1為汛限水位或者溢洪道底坎高程對(duì)應(yīng)的庫(kù)容；W_z0為當(dāng)前庫(kù)水位Z₀相應(yīng)的庫(kù)容；F為水庫(kù)以上流域的面積；α為徑流系數(shù)；k_c為單位換算系數(shù)。

1.1.2 考慮泄流的納雨能力計(jì)算方法

在汛期實(shí)際防洪過(guò)程中，當(dāng)水位超過(guò)汛限水位或者泄洪道底坎高程時(shí)，需要考慮水庫(kù)泄流對(duì)納雨能力計(jì)算的影響^［8-9］。對(duì)于該類(lèi)情況，基于泄流曲線查算庫(kù)水位相應(yīng)的泄流量，同時(shí)估算泄流時(shí)長(zhǎng)，計(jì)算得到預(yù)見(jiàn)期內(nèi)泄流總量，將泄流總量加上空余庫(kù)容后進(jìn)行降雨量折算，則為考慮泄流情況下的納雨能力。

若當(dāng)前水位低于汛限水位或者泄洪道底坎高程時(shí)，計(jì)算庫(kù)水位不超設(shè)計(jì)水位的納雨能力見(jiàn)公式（2）。

式中：P′_c為不超過(guò)設(shè)計(jì)洪水位時(shí)的納雨能力；W_z2為設(shè)計(jì)洪水位Z₂相應(yīng)的庫(kù)容；Q_ave為水庫(kù)平均泄流量；Δt為水庫(kù)泄流時(shí)長(zhǎng)；α₁為水位由當(dāng)前水位漲至汛限水位或者溢洪道底坎高程的平均徑流系數(shù)；α₂為水位由汛限水位或者溢洪道底坎高程漲至設(shè)計(jì)水位的平均徑流系數(shù)；m為單位換算系數(shù)。其中，水庫(kù)泄流時(shí)長(zhǎng)Δt 根據(jù)使水位漲至Z₁的降雨量P₁和預(yù)見(jiàn)期內(nèi)降雨總量P_sum估算得到；Q_ave可由水位達(dá)到設(shè)計(jì)水位時(shí)的泄流量估算得到，見(jiàn)公式（3）～（4）。

式中：T_Len為預(yù)見(jiàn)期時(shí)長(zhǎng)；Q₂為水位達(dá)到設(shè)計(jì)水位時(shí)的泄流量。

若當(dāng)前水位高于汛限水位或者泄洪道底坎高程時(shí)，庫(kù)水位不超設(shè)計(jì)水位的納雨能力計(jì)算公式見(jiàn)式（5）。

式中：水庫(kù)泄流時(shí)長(zhǎng)Δt等于預(yù)見(jiàn)期；Q_ave=0.5×Q₀+Q₂，Q₀為當(dāng)前水位對(duì)應(yīng)的水庫(kù)泄流量。

1.2 考慮未來(lái)降雨影響的納雨能力迭代試算求解

徑流系數(shù)受降雨影響實(shí)時(shí)發(fā)生變化，導(dǎo)致在水庫(kù)納雨能力計(jì)算時(shí)難以直接確定預(yù)見(jiàn)期內(nèi)的平均徑流系數(shù)，不利于量化未來(lái)降雨對(duì)水庫(kù)納雨能力的影響^［10-12］。因此本次研究設(shè)計(jì)納雨能力迭代試算求解方法，先根據(jù)流域下墊面土壤含水量估算納雨能力初始值，然后反查徑流系數(shù)，開(kāi)展納雨能力迭代試算求解，具體步驟包括：

（1）計(jì)算前期影響雨量。對(duì)于以蓄滿(mǎn)產(chǎn)流為主的區(qū)域，影響產(chǎn)流量大小的主要因素為降雨量和初始土壤含水量，就整個(gè)流域而言，初始土壤含水量很難用實(shí)測(cè)土壤含水量值來(lái)表示，因此常用前期影響雨量P_a量化土壤墑情，其計(jì)算方法見(jiàn)式（6）。

P_a，i+1=K_a×P_a，i+P_i（6）

式中：P_a，i，P_a，i+1分別為第 i 天和第 i+1天的前期影響雨量；P_i為第i天的流域降雨量；K_a 為土壤含水量的日消退系數(shù)，每個(gè)月可近似取一個(gè)平均值。

（2）結(jié)合最大蓄水容量I_m計(jì)算得到納雨能力試算開(kāi)始的初始值P_c，init。

（3）根據(jù)初始值P_c，init，結(jié)合小型水庫(kù)P-P_a-α曲線查算對(duì)應(yīng)的徑流系數(shù)α_c，其中查算徑流系數(shù)α?xí)r對(duì)未來(lái)降雨量P進(jìn)行設(shè)置，從而考慮未來(lái)降雨過(guò)程中徑流系數(shù)的變化，進(jìn)而得到納雨能力計(jì)算值P_c，cal。

（4）統(tǒng)計(jì)納雨能力計(jì)算值P_c，cal和初始值P_c，init之間的差距ΔP_c，并根據(jù)ΔP_c迭代試算納雨能力。

ΔP_c=|P_c，cal-P_c，init|（10）

若ΔP_c的值大于試算閾值，則更新納雨能力初始值并再次迭代試算，公式如下：

式中：n為迭代試算的次數(shù)。

若ΔP_n的值小于試算閾值則停止試算，將此時(shí)的P_c，cal作為試算的最終結(jié)果，一般將試算閾值設(shè)置為3 mm。

2 研究區(qū)域概況及數(shù)據(jù)資料

2.1 研究區(qū)域概況

湖南省位于長(zhǎng)江中游，省內(nèi)河網(wǎng)密布，流長(zhǎng)5 km以上的河流5 341條，總長(zhǎng)度9萬(wàn) km，其中流域面積在5 000 km²以上的大河17條。省內(nèi)主要有湘江、資江、沅江、澧水四水及其支流，順著地勢(shì)由南向北匯入洞庭湖、長(zhǎng)江，形成一個(gè)比較完整的洞庭湖水系。計(jì)算涉及32座試點(diǎn)水庫(kù)主要分布在石門(mén)縣、湘鄉(xiāng)市和瀏陽(yáng)市（圖1），均屬于小型水庫(kù)，具有集水面積以及庫(kù)容較小、匯流速度快等特點(diǎn)。

2.2 水庫(kù)及資料概況

本文主要收集了湖南省內(nèi)32個(gè)試點(diǎn)水庫(kù)的死水位、設(shè)計(jì)高水位等特征水位、集水面積、壩高以及溢洪道情況等基本信息。由于試點(diǎn)小水庫(kù)監(jiān)測(cè)設(shè)備布設(shè)時(shí)間不久、監(jiān)測(cè)數(shù)據(jù)較少，本次研究選用的數(shù)據(jù)主要包括2023年6月以來(lái)的逐小時(shí)水位監(jiān)測(cè)記錄、實(shí)測(cè)降雨數(shù)據(jù)（由地面觀測(cè)站點(diǎn)獲得）、數(shù)值降雨預(yù)報(bào)數(shù)據(jù)（通過(guò)歐洲中期天氣預(yù)報(bào)中心得到）。

3 缺資料地區(qū)納雨能力模型參數(shù)方案構(gòu)建

納雨能力模型的主要參數(shù)見(jiàn)表1。

徑流系數(shù)受到流域自然地理?xiàng)l件、前期雨量的影響，較為復(fù)雜多變，是計(jì)算水庫(kù)納雨能力的核心參數(shù)。根據(jù)水利部門(mén)主編的GB/T 50095-2014《水文基本術(shù)語(yǔ)和符號(hào)標(biāo)準(zhǔn)》^［13］，定義徑流系數(shù)為時(shí)段或次流域平均徑流深與相應(yīng)降雨量的比值，是影響產(chǎn)流的重要參數(shù)之一。

基于湖南省水利廳2015年出版的《湖南省暴雨洪水查算手冊(cè)》^［14］進(jìn)行水文分區(qū)，在各個(gè)分區(qū)內(nèi)統(tǒng)計(jì)最大蓄水量、前期土壤含水量以及最大損失量，計(jì)算得到一定前期影響雨量P_a條件下土壤飽和后的徑流系數(shù)α與降雨量P之間的穩(wěn)定關(guān)系，然后開(kāi)展鄰近流域P-P_a-R曲線的移用，在曲線移用的基礎(chǔ)上計(jì)算得到P-P_a-α曲線，最后根據(jù)徑流系數(shù)α與降雨量P之間的穩(wěn)定關(guān)系對(duì)移用后的P-P_a-α曲線進(jìn)行優(yōu)化，用以滿(mǎn)足32座典型小水庫(kù)的納雨能力參數(shù)方案編制的需求。

3.1 湖南省水文產(chǎn)流分區(qū)

在查算手冊(cè)的基礎(chǔ)上，分析湖南省不同區(qū)域降雨量與徑流系數(shù)的關(guān)系。手冊(cè)將全省共劃分為5個(gè)水文分區(qū)（圖2），表2中分別列出了各分區(qū)的參數(shù)最大蓄水量（I_m）及最大損失量（I_o）。

通過(guò)最大損失量I_o可比較出各分區(qū)流域的降雨損失情況。地勢(shì)較為低平的Ⅱ區(qū)和Ⅴ區(qū)降雨損失量較大，地形起伏的Ⅰ區(qū)、Ⅲ區(qū)和Ⅳ區(qū)降雨損失量較小，更易于產(chǎn)流。

3.2 徑流系數(shù)與降雨量間的穩(wěn)定關(guān)系

根據(jù)查算手冊(cè)附表可繪制出全省5個(gè)分區(qū)內(nèi)降雨量與徑流系數(shù)關(guān)系（圖2），各分區(qū)確定P_a后，則可建立P=F（R）曲線。當(dāng)降雨達(dá)到一定量時(shí)，降雨使包氣帶達(dá)到飽和，流域損失量趨于穩(wěn)定，此后徑流系數(shù)與累積降雨量形成穩(wěn)定的線性關(guān)系（表2）。

3.3 鄰近流域的P-P_a-R曲線移用

在數(shù)據(jù)資料較為缺乏的條件下，借用鄰近流域的P-P_a-R曲線開(kāi)展小水庫(kù)參數(shù)方案優(yōu)化。選用湘贛邊界的高安站、萬(wàn)家埠站作為瀏陽(yáng)市小水庫(kù)的移用站點(diǎn)，選用清江流域的隔河巖站、高壩洲站作為石門(mén)縣小水庫(kù)的移用站點(diǎn)（圖3）。計(jì)算高安站、萬(wàn)家埠站、隔河巖站、高壩洲站的P-P_a-R關(guān)系曲線的斜率，即查詢(xún)一定降雨量P、前期影響雨量P_a條件下的徑流R，然后計(jì)算R與P的比值，得到徑流系數(shù)α，進(jìn)而將P-P_a-R曲線轉(zhuǎn)換為P-P_a-α曲線（圖4）。

3.4 缺資料地區(qū)試點(diǎn)水庫(kù)參數(shù)方案構(gòu)建

分別綜合高安站、萬(wàn)家埠站的P-P_a-R曲線，隔河巖站、高壩洲站的P-P_a-R曲線，并結(jié)合《湖南省暴雨洪水查算手冊(cè)》對(duì)小水庫(kù)的P-P_a-R曲線進(jìn)行優(yōu)化。

修編時(shí)首先設(shè)置一系列降雨量，通過(guò)表2中徑流系數(shù)與累積降雨量形成的穩(wěn)定線性關(guān)系，計(jì)算得到降雨量和徑流系數(shù)的組合［P_n，α_s，n］，將其作為計(jì)算參考點(diǎn)；在移用站點(diǎn)綜合P-P_a-R曲線上查算與參考點(diǎn)相同降雨P(guān)_n以及相同前期影響雨量P_a條件下的徑流系數(shù)α_c，n；然后計(jì)算α_s，n與α_c，n之間的比值，并計(jì)算得到其均值β，將綜合P-P_a-R曲線的各個(gè)節(jié)點(diǎn)乘以比例β；同時(shí)調(diào)整P-P_a-R曲線使得降雨量P超過(guò)流域最大蓄水量I_m時(shí)的P-α關(guān)系線與圖5中的黑色虛線基本平行，并且將P_a最大值外延插補(bǔ)至I_m值，從而完成P-P_a-R曲線優(yōu)化（圖5），使其適用于小水庫(kù)所在的水文分區(qū)，其中圖5黑色虛線即表示徑流系數(shù)與累積降雨量形成的穩(wěn)定線性關(guān)系，紅色點(diǎn)即為參考點(diǎn)。

4 模擬結(jié)果及分析

4.1 基于納雨能力分析的小型水庫(kù)防洪風(fēng)險(xiǎn)分布

設(shè)置納雨能力預(yù)見(jiàn)期為未來(lái)24 h，將前15 d作為其計(jì)算前期影響雨量P_a的預(yù)熱期。結(jié)合起始時(shí)刻各小型水庫(kù)的實(shí)時(shí)水位Z_c、前期影響雨量P_a，計(jì)算得到預(yù)見(jiàn)期內(nèi)分別以汛限水位和設(shè)計(jì)水位為特征水位的納雨能力值，并與預(yù)見(jiàn)期內(nèi)ECMWF降雨預(yù)報(bào)值進(jìn)行對(duì)比（圖6）。

由圖6可以看出，以2023年8月21日為起始時(shí)刻，計(jì)算得到的至水庫(kù)泄流底坎高程的納雨能力P_cx均大于預(yù)見(jiàn)期內(nèi)的降雨量P_f，即反映未來(lái)24 h內(nèi)各水庫(kù)的水位超過(guò)泄流底坎高程的可能性較??；至水庫(kù)設(shè)計(jì)水位的納雨能力P_cx遠(yuǎn)大于P_cx和P_f，這主要是由于當(dāng)水位超過(guò)泄流底坎后水庫(kù)將開(kāi)始泄流，導(dǎo)致P_cx的值較大。

計(jì)算納雨能力P_cx與預(yù)見(jiàn)期降雨P(guān)_f之間的差值ΔP_x，用以反映水庫(kù)防洪風(fēng)險(xiǎn)，通常ΔP_x的值越大，防洪風(fēng)險(xiǎn)越低，反之風(fēng)險(xiǎn)越高，試點(diǎn)水庫(kù)的防洪風(fēng)險(xiǎn)分布見(jiàn)圖7。可以明顯看出，2023年8月21～22日期間，預(yù)報(bào)的降雨落區(qū)主要分布在長(zhǎng)沙至瀏陽(yáng)一帶，該區(qū)域內(nèi)的小水庫(kù)防洪風(fēng)險(xiǎn)較高，其中廖家沖、大山、豎溪、金山以及金盆水庫(kù)因其當(dāng)前水位相對(duì)較高、未來(lái)降雨量較大，導(dǎo)致ΔP_x值均在30 mm以下，屬于該輪降雨過(guò)程中需要重點(diǎn)關(guān)注的小型水庫(kù)。相對(duì)而言，位于石門(mén)縣內(nèi)的和平、白果等水庫(kù)未來(lái)降雨較小，ΔP_x值均在60 mm以上，其防洪風(fēng)險(xiǎn)較低。由此可以看出，本研究提出的方法可以較直觀地反映小型水庫(kù)的防洪風(fēng)險(xiǎn)分布，能夠有力支持小型水庫(kù)防洪決策以及風(fēng)險(xiǎn)動(dòng)態(tài)控制。

4.2 典型水庫(kù)納雨能力計(jì)算結(jié)果

除了小型水庫(kù)防洪風(fēng)險(xiǎn)整體把握和分布直觀展現(xiàn)外，根據(jù)實(shí)時(shí)變化的降雨以及水庫(kù)水位滾動(dòng)追蹤納雨能力動(dòng)態(tài)變化，也是開(kāi)展小型水庫(kù)洪水預(yù)警的重要環(huán)節(jié)。以位于石門(mén)縣的和平水庫(kù)為例，逐日滾動(dòng)計(jì)算2023年8月1日至9月30日的納雨能力，并分析納雨能力與水位、前期影響雨量的相關(guān)性（圖8～9）。

圖8和圖9中A、B、C以及D表示的時(shí)刻及其對(duì)應(yīng)的納雨能力、水位、前期影響雨量見(jiàn)表3。

4.3 計(jì)算結(jié)果分析

為分析納雨能力受水庫(kù)當(dāng)前水位以及前期影響雨量的影響，查算各水庫(kù)當(dāng)前水位至泄流底坎之間的庫(kù)容ΔV_x，并將P_cx、ΔV_x以及P_a各變量折算為流域面平均雨量，采用同一縱軸進(jìn)行量化表示，橫軸則為各水庫(kù)（圖10）。

分析得到水庫(kù)P_cx與ΔV_x之間、P_cx與P_a之間的皮爾遜相關(guān)系數(shù)（Pearson correlation coefficient）值分別為0.36和-0.64，可以反映出試點(diǎn)水庫(kù)的水位距泄流底坎高程較近、前期影響雨量較大時(shí)，計(jì)算得到的水庫(kù)納雨能力較小，這與小型水庫(kù)預(yù)警實(shí)際情況相符，反映了本文計(jì)算方法的合理性。

根據(jù)典型小水庫(kù)納雨能力計(jì)算結(jié)果，在2023年8月18日前，和平水庫(kù)所在區(qū)域沒(méi)有降雨，前期影響雨量不斷減少，土壤含水量下降，且?guī)焖怀掷m(xù)消退，計(jì)算得到的納雨能力值不斷增加，至8月18日水庫(kù)納雨能力達(dá)到高點(diǎn)123.6 mm；8月18～22日，出現(xiàn)一次10.6 mm左右的降雨，雖然前期影響雨量?jī)H小幅增加約13 mm，但是庫(kù)水位迅速增加，其漲幅達(dá)到1.3 m，導(dǎo)致水庫(kù)納雨能力快速下降至89.4 mm；8月22～27日，雖然出現(xiàn)降雨過(guò)程，但是由于人工干預(yù)使庫(kù)水位持續(xù)下降，并且受到前期影響雨量增大至高點(diǎn)102 mm的影響，計(jì)算得到的水庫(kù)納雨能力降至72.4 mm；8月27至9月21日，連續(xù)出現(xiàn)2次降雨過(guò)程，前期影響雨量波動(dòng)攀升至116.8 mm，但由于庫(kù)水位波動(dòng)下降，9月21日的納雨能力為68 mm，較8月27日變幅較小。

因此可以看出，水庫(kù)納雨能力受到庫(kù)水位和前期降雨過(guò)程的綜合影響，不僅可以分布展現(xiàn)基于納雨能力的小型水庫(kù)防洪風(fēng)險(xiǎn)，而且可以考慮不同時(shí)期影響納雨能力主要因素的動(dòng)態(tài)變化。本研究提出的方法能及時(shí)追蹤水庫(kù)納雨能力的變化，實(shí)時(shí)滾動(dòng)分析影響納雨能力的因子，實(shí)現(xiàn)了納雨能力動(dòng)態(tài)合理計(jì)算，為小型水庫(kù)滾動(dòng)持續(xù)預(yù)警提供了有力支撐。

5 結(jié) 論

針對(duì)缺資料地區(qū)小型水庫(kù)預(yù)警難題，構(gòu)建適用于不同小水庫(kù)運(yùn)行方式的納雨能力計(jì)算模型，綜合量化降雨以及流域下墊面特征對(duì)于納雨能力模型參數(shù)的影響，開(kāi)發(fā)考慮未來(lái)降雨的納雨能力迭代試算方法，以2023年8月21日為起始時(shí)刻，計(jì)算湖南省缺資料地區(qū)32座試點(diǎn)小型水庫(kù)納雨能力，并以和平水庫(kù)為例開(kāi)展納雨能力滾動(dòng)計(jì)算。結(jié)果表明：水位距泄流底坎高程較近、前期影響雨量較大時(shí)，計(jì)算得到的水庫(kù)納雨能力較小；所建模型合理量化了降雨對(duì)于小型水庫(kù)防洪風(fēng)險(xiǎn)的影響，展現(xiàn)了基于納雨能力的小型水庫(kù)防洪風(fēng)險(xiǎn)空間分布，同時(shí)可以考慮不同時(shí)期影響納雨能力主要因素的動(dòng)態(tài)變化，從而有力提升湖南省小水庫(kù)風(fēng)險(xiǎn)預(yù)警能力，為小水庫(kù)安全度汛提供了有效支撐。

參考文獻(xiàn)：

［1］中華人民共和國(guó)水利部.2023中國(guó)水利發(fā)展報(bào)告［M］.北京：中國(guó)水利水電出版社，2023.

［2］劉林，李國(guó)文，凍芳芳，等.中小型水庫(kù)納雨能力計(jì)算方法研究［J］.中國(guó)防汛抗旱，2022，32（8）：66-71.

［3］國(guó)洪琴.基于不同預(yù)設(shè)條件的水庫(kù)抗暴雨能力淺析［J］.地下水，2016，38（2）：124-125.

［4］侯?lèi)?ài)中，胡智丹，朱冰，等.中小型水庫(kù)抗暴雨能力的概念與計(jì)算方法［J］.水文，2017，37（6）：35-38.

［5］章四龍，侯?lèi)?ài)中，吳志勇，等.中小型水庫(kù)抗暴雨能力關(guān)鍵技術(shù)研究［M］.北京：中國(guó)水利水電出版社，2019.

［6］王漢東，黃瓅瑤，王燕，等.基于預(yù)報(bào)調(diào)度一體化計(jì)算的水庫(kù)納雨能力分析［J］.水利水電快報(bào)，2022，43（12）：31-35.

［7］王剛，任明磊，李京兵，等.小型水庫(kù)防洪預(yù)警指標(biāo)體系設(shè)計(jì)及應(yīng)用研究［J］.中國(guó)防汛抗旱，2023，33（3）：51-55.

［8］王偉，王樂(lè)，田逸飛，等.丹江口水庫(kù)入庫(kù)洪水概率預(yù)報(bào)及調(diào)度風(fēng)險(xiǎn)分析［J］.人民長(zhǎng)江，2023，54（11）：60-67.

［9］林璐.不同納雨能力計(jì)算方法的應(yīng)用對(duì)比研究：以江門(mén)市大型水庫(kù)為例［J］.廣東水利水電，2021（8）：14-19.

［10］邱輝，范維，王樂(lè).漢江流域多模式面雨量預(yù)報(bào)效果評(píng)估及集成應(yīng)用［J］.人民長(zhǎng)江，2023，54（11）：51-59.

［11］楊欣磊，楊鵬.海河流域河北省中型水庫(kù)納雨能力計(jì)算分析［J］.海河水利，2023（1）：63-65.

［12］王海芝，許冰，翟淑花，等.北京市暴雨-突發(fā)地質(zhì)災(zāi)害量化關(guān)系分析［J］.人民長(zhǎng)江，2023，54（10）：81-89.

［13］中華人民共和國(guó)水利部.水文基本術(shù)語(yǔ)和符號(hào)標(biāo)準(zhǔn)：GB/T 50095-2014［S］.北京：中國(guó)計(jì)劃出版社，2014.

［14］湖南省水利廳.湖南省暴雨洪水查算手冊(cè)［Z］.長(zhǎng)沙：湖南省水利廳，2015.

（編輯：謝玲嫻）

Construction and application of rainfall retaining capacity calculation model

of small reservoirs in ungauged region

TONG Bingxing^1，2，3，ZHANG Xiao^1，2，3，ZHANG Jun^{1，2，3}，LI Jie^{1，2，3}，JIANG Dongqing⁴，LYU Kun⁴

（1.Bureau of Hydrology，Changjiang Water Resources Commission，Wuhan 430010，China; 2.Innovation Team for Flood and Drought Disaster Prevention of Changjiang Water Resources Commission，Wuhan 430010，China; 3.Center of Technology Innovation for Intelligent Water Simulation，F(xiàn)orecasting and Operation of the River Basin，Wuhan 430010，China; 4.Hu′nan Provincial Hydrology and Water Resources Survey Center，Changsha 410000，China）

Abstract：To realize the flood risk warning of small reservoirs in ungauged regions，a rainfall retaining capacity calculation model considering different operation modes of small reservoirs and the model parameter schemes have been developed in this study.An iterative trial solution method considering potential rainfall for the constructed model has been proposed for efficient computation.The rainfall retaining capacity of 32 pilot small reservoirs located in the ungauged region of Hu′nan Province was calculated starting from August 21，2023.In addition，a continuous calculation was carried out from August 1 to September 30，2023，using the proposed method for Heping Reservoir.The results show that the flood risk of the pilot small reservoirs calculated is reasonably related to the rainfall distribution，supporting flood risk visualization for small reservoirs within the rain-affected region.Furthermore，the impact of reservoir water level and antecedent precipitation on rainfall retaining capacity can also be dynamically quantified in real-time，which provides a reference for the continuous monitoring of flood risks of small reservoirs.

Key words：small reservoirs; rainfall retaining capacity; antecedent precipitation; iterative trial calculations; ungauged region; Hunan Province