岳城水庫(kù)入庫(kù)徑流與觀臺(tái)站實(shí)測(cè)徑流對(duì)比分析

王玉梅

（河北省邯鄲水文勘測(cè)研究中心，河北 邯鄲 056001）

2021 年漳河流域普降暴雨，洪水場(chǎng)次多，持續(xù)時(shí)間長(zhǎng)，出現(xiàn)罕見(jiàn)的夏秋聯(lián)汛。“7·20”暴雨洪水中，觀臺(tái)斷面實(shí)測(cè)洪峰流量3150 m3/s，而下游岳城水庫(kù)反推入庫(kù)洪峰流量4860 m3/s，二者相差超過(guò)54%；10 月初，觀臺(tái)斷面沖刷下切，岳城水庫(kù)超過(guò)歷史最高蓄水位，對(duì)觀臺(tái)斷面造成頂托，影響測(cè)驗(yàn)精度。故本文對(duì)岳城水庫(kù)反推入庫(kù)徑流與觀臺(tái)站實(shí)測(cè)徑流進(jìn)行對(duì)比，分析二者合理性。

1 流域概述

漳河位于海河流域西南部，是漳衛(wèi)南運(yùn)河的重要組成部分之一，發(fā)源于山西省太行山腹地，涉及山西、河北、河南三省，全長(zhǎng)460 km，流域面積18284 km2［1-3］。漳河有濁、清漳兩大支流，太行山主脈南北貫穿流域中部，將流域分東、西兩個(gè)不同自然地理特征的區(qū)域。東區(qū)包括辛安以下濁漳河干流、清漳河大部（除西支）及漳河干流，為石質(zhì)山區(qū)，山高谷深，巖石裸露，坡陡流急，含沙量小；西區(qū)包括濁漳河辛安以上，為山丘和盆地區(qū)，盆地內(nèi)黃土覆蓋較厚，植被較差，水土流失嚴(yán)重，洪水挾帶泥沙較多。清漳河和濁漳河在河北省涉縣合漳匯合成漳河干流，繼續(xù)沿河北、河南兩省的邊界在山谷中蜿蜒穿行，兩岸分別是河南省的林州市、安陽(yáng)縣，河北省的涉縣、磁縣。

觀臺(tái)水文站是漳河干流出山口控制站，也是岳城水庫(kù)入庫(kù)站，1941 年設(shè)站，測(cè)站位于磁縣觀臺(tái)鎮(zhèn)冶子村漳河邊，控制面積17800 km2，距離下游岳城水庫(kù)大壩17 km。

岳城水庫(kù)是漳河干流出山口處的一座大（Ⅰ）型水庫(kù)，水庫(kù)總庫(kù)容13 億m3，控制流域面積18100 km2，占漳河流域面積的99.4%，主要功能為防洪，兼有灌溉、發(fā)電、城市供水等綜合效益，為不完全年調(diào)節(jié)的大型水利樞紐工程。

觀臺(tái)～岳城水庫(kù)區(qū)間徑流面積300 km2，存在都黨溝、黃沙溝等大小不同河溝，區(qū)間產(chǎn)流直接匯入岳城水庫(kù)。

2 洪水分析

2.1 暴雨洪水分析

2.1.1 “7·11”暴雨洪水

隨著副熱帶高壓北抬，雨帶北移，2021 年7 月10—12 日岳城水庫(kù)以上流域發(fā)生強(qiáng)降水過(guò)程，局部暴雨或大暴雨。區(qū)域平均降水量55.7 mm，主要集中在下游，合漳至岳城水庫(kù)壩前部分，最大日雨量為河北省涉縣的匡門(mén)口站182.0 mm，最大1 h 降水量82.0 mm。區(qū)域降水量等值線(xiàn)如圖1。

圖1 “7·11”降水等值線(xiàn)

此次降水形成漳河2021 年1 號(hào)洪水。根據(jù)觀臺(tái)站實(shí)測(cè)徑流，選用7 月11 日23：00 至20 日22：10 徑流過(guò)程，實(shí)測(cè)洪峰流量163 m3/s，徑流量2712 萬(wàn)m3；岳城水庫(kù)對(duì)應(yīng)入庫(kù)時(shí)間段為7 月12 日09：00 至21 日08：00，反 推 入 庫(kù) 徑 流 量3649 萬(wàn)m3，二 者 相 差937 萬(wàn)m3。

根據(jù)岳城水庫(kù)壩上水位～水面面積關(guān)系表，可以推出此段時(shí)間內(nèi)岳城水庫(kù)水面面積約17.44 km2。根據(jù)圖2 可知，“7·11”暴雨過(guò)程中岳城水庫(kù)水面平均降水量為115.9 mm，則計(jì)算可得，岳城水庫(kù)庫(kù)區(qū)水面接納降水總量為202 萬(wàn)m3。

圖2 “7·20”降水等值線(xiàn)

扣除岳城水庫(kù)水面面積17.44 km2，觀臺(tái)～岳城水庫(kù)區(qū)間匯流面積為282.56 km2，區(qū)間平均降水量134.1 mm，則降水徑流系數(shù)為0.19。

2.1.2 “7·20”暴雨洪水

7 月以來(lái)，西北太平洋副熱帶高壓、大陸高壓分別穩(wěn)定維持在日本海和我國(guó)西北地區(qū)，長(zhǎng)時(shí)間阻擋上游系統(tǒng)的移動(dòng)，其次受深厚的東風(fēng)急流河穩(wěn)定低渦切變影響，加上太行山東側(cè)特殊地形抬升作用，使得河南省北部、西北部及河北南部地區(qū)發(fā)生強(qiáng)降水。7 月10—19 日期間，岳城水庫(kù)以上流域，尤其是中下游地區(qū)連續(xù)陰雨天氣，土壤趨于飽和。7 月20—22 日岳城水庫(kù)以上流域普降暴雨，區(qū)域平均降水量142.0 mm，降水主要集中在中下游，尤其是合漳～岳城水庫(kù)區(qū)間；最大日雨量為河北省磁縣的中清流站523.0 mm，最大1 h 降水量63.5 mm，最大24 h 降水量536.5 mm。區(qū)域降水量等值線(xiàn)如圖2。

此次暴雨過(guò)程形成漳河2021 年第2 號(hào)洪水。根據(jù)觀臺(tái)站實(shí)測(cè)徑流，選用7 月20 日22：10～24 日13：00 徑流過(guò)程，實(shí)測(cè)洪峰流量3150m3/s，發(fā)生在7月22 日07：36，徑流量26426 萬(wàn)m3；岳城水庫(kù)對(duì)應(yīng)入庫(kù)徑流過(guò)程為7 月21 日08：00—24 日18：00，反推入庫(kù)最大洪峰4860 m3/s，發(fā)生在7 月22 日08：30，入庫(kù)徑流量為32927 萬(wàn)m3，二者相差6501 萬(wàn)m3。觀臺(tái)站至岳城水庫(kù)壩前洪峰傳播時(shí)間約為1 h。

根據(jù)岳城水庫(kù)壩上水位～水面面積關(guān)系表，可推出此段時(shí)間內(nèi)岳城水庫(kù)水面面積約25.13 km2。據(jù)圖3 可知，“7·20”過(guò)程中岳城水庫(kù)庫(kù)區(qū)水面平均降水量為465.3 mm，則岳城水庫(kù)水面接納降水總量為1378 萬(wàn)m3。

圖3 “10·07”降水等值線(xiàn)

2.1.3 “10·07”暴雨洪水

異常偏強(qiáng)的西太平洋副熱帶高壓先西伸北抬后穩(wěn)定維持在黃淮地區(qū)，與西風(fēng)帶低值系統(tǒng)在山西形成穩(wěn)定的東高西低的環(huán)流形勢(shì)，副熱帶高壓西側(cè)的偏南氣流和低層西南急流將南海和孟加拉灣的水汽向北經(jīng)過(guò)西南地區(qū)源源不斷輸送到山西中南部地區(qū)，提供充沛的水汽來(lái)源，造成漳河上游地區(qū)產(chǎn)生強(qiáng)降水過(guò)程。

此次降水過(guò)程從10 月2 日持續(xù)到6 日，主要發(fā)生在濁漳河、清漳河上游地區(qū)，山西省境內(nèi)，區(qū)域平均降水量99.2 mm；最大日雨量為山西省晉中市榆社縣的社城站84.8 mm，最大1 h 降水量為山西省長(zhǎng)子縣長(zhǎng)子站35.0 mm，最大24 h 降水量97.6 mm（社城站）。區(qū)域降水量等值線(xiàn)如圖3。

此次暴雨過(guò)程形成漳河2021 年第3 號(hào)洪水，主要為上游來(lái)水及區(qū)間基流量。根據(jù)觀臺(tái)站實(shí)測(cè)徑流，選用10 月5 日23：45—15 日11：30 徑流過(guò)程，實(shí)測(cè)洪峰流量2240 m3/s，發(fā)生在7 日11：34，徑流量71452 萬(wàn)m3；岳城水庫(kù)對(duì)應(yīng)入庫(kù)徑流過(guò)程為10 月6日04：00～15 日16：00，反推入庫(kù)最大洪峰2070 m3/s，入庫(kù)徑流量為72455 萬(wàn)m3，二者相差1003 萬(wàn)m3。

2.2 產(chǎn)匯流分析

2.2.1 月年徑流分析

7 月11 日之前觀臺(tái)站處于河干狀態(tài)，故本次僅對(duì)比7—12 月份的徑流量。根據(jù)岳城水庫(kù)來(lái)（蓄）水量月年統(tǒng)計(jì)表可知，岳城水庫(kù)下半年反推入庫(kù)徑流量為28.29 億m3。觀臺(tái)站年徑流量26.31 億m3，如表1。根據(jù)岳城水庫(kù)壩上水位～水面面積關(guān)系表，可推出各月岳城水庫(kù)平均水面面積，從而計(jì)算出7—12 月水庫(kù)水面接納降水總量為3376 萬(wàn)m3，如表1。

表1 岳城水庫(kù)逐月徑流和水面接納降水總量

根據(jù)邯鄲水文中心2021 年暴雨洪水調(diào)查報(bào)告可知，都黨溝、黃沙溝、寺溝均位于觀臺(tái)～岳城水庫(kù)未控區(qū)間，流域面積分別為128.4，20.1，6.26 km2，合計(jì)占觀臺(tái)～岳城水庫(kù)未控區(qū)間面積的56.3%。暴雨洪水調(diào)查面積154.76 km2，產(chǎn)流3726.2 萬(wàn)m3，根據(jù)面積同倍比縮放，觀臺(tái)～岳城水庫(kù)區(qū)間7 月份降水產(chǎn)流6803 萬(wàn)m3。另有水量調(diào)查斷面黃沙，位于觀臺(tái)站下游，黃沙礦場(chǎng)長(zhǎng)期排放疏干水，直接排入岳城水庫(kù)，據(jù)調(diào)查結(jié)果顯示，黃沙斷面7—12 月徑流量為1457 萬(wàn)m3。

岳城水庫(kù)入庫(kù)徑流應(yīng)包含觀臺(tái)站實(shí)測(cè)徑流、水庫(kù)水面接納降水總量、區(qū)間暴雨產(chǎn)流、黃沙疏干水。對(duì)以上各項(xiàng)進(jìn)行逐月加和計(jì)算，可知岳城水庫(kù)反推入庫(kù)徑流量較各項(xiàng)相加值偏大8275 萬(wàn)m3，偏大2.92%，如表2。

表2 岳城水庫(kù)逐月入庫(kù)徑流誤差分析單位：萬(wàn)m3

8 月份，岳城水庫(kù)反推入庫(kù)徑流小于各項(xiàng)相加值，偏小5.21%，其余各月均偏大，尤其是11 月份，岳城水庫(kù)反推入庫(kù)徑流偏大5.28%，其余月份差值均未超過(guò)5%。

通過(guò)圖12對(duì)比發(fā)現(xiàn)，灌漿套筒連接件破壞過(guò)程分為彈性、屈服、強(qiáng)化、頸縮4個(gè)階段，試件H500-14承載力高出H400-14承載力10%左右，表明鋼筋強(qiáng)度的變化對(duì)試件承載力影響明顯，鋼筋強(qiáng)度越大，試件承載力越高。

2.2.2 次洪徑流分析

根據(jù)上游來(lái)水、降水場(chǎng)次及岳城水庫(kù)放水情況，對(duì)岳城水庫(kù)反推入庫(kù)徑流過(guò)程進(jìn)行場(chǎng)次劃分，共劃分10 個(gè)階段，扣除水庫(kù)水面接納降水量、區(qū)間產(chǎn)匯流及黃沙調(diào)查斷面匯入徑流量，與觀臺(tái)站實(shí)測(cè)徑流進(jìn)行對(duì)比分析。對(duì)比情況如表3。

表3 岳城水庫(kù)與觀臺(tái)站分階段徑流對(duì)比分析 單位：萬(wàn)m3

根據(jù)表3 可知，2021 年觀臺(tái)站與岳城水庫(kù)反推入庫(kù)徑流過(guò)程共分為10 個(gè)階段，分別對(duì)其進(jìn)行分析。

第1 階段與第2 階段岳城水庫(kù)反推入庫(kù)徑流分別對(duì)應(yīng)觀臺(tái)站“7·11”“7·20”洪水過(guò)程，岳城水庫(kù)反推入庫(kù)徑流與觀臺(tái)站相差較大，單位差分別為8.31 m3/s、-17.1 m3/s，相差較大。

第3 階段（7 月24 日至8 月2 日）：觀臺(tái)站處于落水階段，7 月23 日17：00 岳城水庫(kù)開(kāi)閘泄洪，最大泄量700 m3/s （27 日9：00），28 日23：00 停止泄洪。岳城水庫(kù)反推入庫(kù)徑流較觀臺(tái)站偏大34.3 m3/s，差值較大。

第4 階段（8 月2 日至9 月19 日）：觀臺(tái)站處于低水期，岳城水庫(kù)未泄洪，岳城水庫(kù)反推入庫(kù)徑流與觀臺(tái)站相差1.96 m3/s。

第5 階段（9 月19 日至9 月28 日）：9 月18 日18 時(shí)岳城水庫(kù)開(kāi)始泄洪，泄洪流量逐漸加大，至27日11:00 漲至502 m3/s。觀臺(tái)站經(jīng)歷“9·20”“9·26”兩次漲水過(guò)程。岳城水庫(kù)反推入庫(kù)徑流較觀臺(tái)站偏大20.0 m3/s，差值較大。

第6 階段（9 月28 日至10 月6 日）：岳城水庫(kù)蓄水位上漲，泄洪流量相對(duì)穩(wěn)定，保持在450～552 m3/s之間，觀臺(tái)站受水庫(kù)蓄水頂托影響明顯，處于緩慢落水期，岳城水庫(kù)反推入庫(kù)徑流較觀臺(tái)站相差無(wú)幾。

第7 階段（10 月6 日至10 月15 日）：觀臺(tái)站迎來(lái) “10·07”漲水過(guò)程，10 月7 日11:34 達(dá)到洪峰2240 m3/s；10 月6 日4:00 岳城水庫(kù)調(diào)整出庫(kù)，逐漸加大，10 月7 日54:31 穩(wěn)定在848 m3/s，至10 月11日12:00 達(dá)到最大泄量905m3/s，隨后緩慢減小泄量。岳城水庫(kù)反推入庫(kù)徑流較觀臺(tái)站偏大9.72 m3/s。

第8 階段（10 月15 日至10 月25 日）：觀臺(tái)站處于緩慢落水期，岳城水庫(kù)泄洪流量從650 m3/s迅速減至150 m3/s，反推入庫(kù)徑流較觀臺(tái)站偏大21.9 m3/s。

第9 階段（10 月25 日至11 月26 日）：觀臺(tái)站處于緩慢落水期，岳城水庫(kù)泄洪流量相對(duì)穩(wěn)定，保持在60.4～150 m3/s 之間，反推入庫(kù)徑流較觀臺(tái)站偏大4.25 m3/s。

第10 階段（11 月26 日至12 月31 日）觀臺(tái)站處于緩慢落水期，岳城水庫(kù)泄洪流量相對(duì)穩(wěn)定，保持在30～70.4 m3/s 之間，反推入庫(kù)徑流較觀臺(tái)站偏小0.290 m3/s。

2.3 誤差分析

2.3.1 暴雨徑流

根據(jù)邯鄲水文中心2021 年暴雨洪水調(diào)查報(bào)告可知，都黨溝、黃沙溝、寺溝均位于觀臺(tái)～岳城水庫(kù)未控區(qū)間，流域面積分別為128.4，20.1，6.26 km2，合計(jì)占觀臺(tái)～岳城水庫(kù)未控區(qū)間面積的56.3%?！?·20”期間，都黨溝、黃沙溝、寺溝流域洪峰流量分別為446，164，64.4 m3/s，降水徑流系數(shù)分別為0.50，0.69，0.72。則推算出的“7·20”觀臺(tái)～岳城水庫(kù)區(qū)間降水徑流系數(shù)0.54，與暴雨洪水調(diào)查報(bào)告中計(jì)算得到的降水徑流系數(shù)較為接近，相對(duì)合理。

2.3.2 次洪

根據(jù)次洪徑流分析可知，岳城水庫(kù)2021 年10個(gè)階段徑流過(guò)程中，第1～3、5、8 階段反推入庫(kù)流量與觀臺(tái)站實(shí)測(cè)徑流相差較大，差值超過(guò)5%，其余階段差值較小，如表4。

表4 岳城水庫(kù)反推入庫(kù)與觀臺(tái)站徑流差值分析

第1 階段岳城水庫(kù)反推入庫(kù)徑流對(duì)應(yīng)觀臺(tái)站“7·11”暴雨洪水過(guò)程。由于沒(méi)有“7·11”暴雨洪水調(diào)查資料，故第1 階段未扣除區(qū)間暴雨洪水產(chǎn)流，造成第1 階段岳城水庫(kù)反推入庫(kù)徑流與觀臺(tái)站實(shí)測(cè)徑流相差較大。

第2 階段岳城水庫(kù)反推入庫(kù)徑流對(duì)應(yīng)觀臺(tái)站“7·20”暴雨洪水過(guò)程。由邯鄲水文中心2021 年暴雨洪水調(diào)查報(bào)告可知，由于觀臺(tái)～岳城水庫(kù)區(qū)間右岸隸屬于河南省，且無(wú)雨量站點(diǎn)，未做調(diào)查，調(diào)查斷面流域面積僅位于漳河左岸，覆蓋觀臺(tái)～岳城水庫(kù)區(qū)間未控面積的56.3%。故第2 階段區(qū)間暴雨洪水產(chǎn)流估算不準(zhǔn)確，造成第2 階段岳城水庫(kù)反推入庫(kù)徑流與觀臺(tái)站實(shí)測(cè)徑流相差較大。

對(duì)10 個(gè)階段差值與岳城水庫(kù)出庫(kù)泄洪流量過(guò)程進(jìn)行對(duì)比，如圖4。

圖4 岳城水庫(kù)反推入庫(kù)差值與泄洪流量過(guò)程對(duì)比

由圖4 可知，第3 階段，岳城水庫(kù)泄洪從400 m3/s逐步加大至700 m3/s，28 日23:00 下降到300 m3/s，隨后閉閘；第5 階段，岳城水庫(kù)開(kāi)始泄洪，泄洪流量從50.2 m3/s 加大至300 m3/s 后，下降至150 m3/s，隨后又加大至502 m3/s，對(duì)應(yīng)觀臺(tái)站“9·20”“9·26”兩次漲水過(guò)程；第8 階段，觀臺(tái)站處于緩慢落水期，岳城水庫(kù)泄洪流量從650 m3/s 持續(xù)減小至150 m3/s。以上3個(gè)階段觀臺(tái)站均處于穩(wěn)定落水階段，岳城水庫(kù)啟閉閘門(mén)較為頻繁，泄洪流量不穩(wěn)定，岳城水庫(kù)反推入庫(kù)徑流與觀臺(tái)站實(shí)測(cè)徑流相差較大。

由圖4 可知，第6 階段，岳城水庫(kù)泄洪流量相對(duì)穩(wěn)定，保持在450～552 m3/s 之間；第7 階段，岳城水庫(kù)調(diào)整出庫(kù)，由450 m3/s 升至905 m3/s，隨后調(diào)小出庫(kù)，降至650 m3/s，對(duì)應(yīng)觀臺(tái)站“10·07”漲水過(guò)程；第9 階段，岳城水庫(kù)泄洪流量相對(duì)穩(wěn)定，保持在80.4～150 m3/s 之間；第10 階段，岳城水庫(kù)泄洪流量相對(duì)穩(wěn)定，保持在30.0～70.4 m3/s 之間。以上4 個(gè)階段，觀臺(tái)站徑流有漲落變化，但岳城水庫(kù)閘門(mén)變動(dòng)不頻繁，時(shí)段內(nèi)泄洪流量相對(duì)穩(wěn)定，岳城水庫(kù)反推入庫(kù)徑流與觀臺(tái)站實(shí)測(cè)徑流相差較小。

2021 年岳城水庫(kù)蓄水位達(dá)152.30 m，蓄水量8.401 億m3，為建庫(kù)以來(lái)歷史最高水位，同時(shí)觀臺(tái)站斷面受采砂影響下切，受水庫(kù)蓄水頂托影響顯著。岳城水庫(kù)的高蓄水與頂托影響均對(duì)觀臺(tái)斷面水流流態(tài)產(chǎn)生影響，尤其是流速垂線(xiàn)分布發(fā)生變化，進(jìn)而影響到實(shí)測(cè)流量，造成一定誤差。

岳城水庫(kù)前期采用泄洪洞泄洪，10 月初采用泄洪洞與溢洪道同時(shí)泄洪，閘門(mén)泄洪曲線(xiàn)啟用于1964年，閘門(mén)全開(kāi)的資料不夠充分，且多年來(lái)未修正，隨著多年上游來(lái)水來(lái)沙及水庫(kù)泄洪洞、溢洪道形態(tài)變化的影響，閘門(mén)泄洪曲線(xiàn)可能存在一定誤差，引起反推入庫(kù)徑流不準(zhǔn)確。

3 結(jié)語(yǔ)

（1）2021 年漳河岳城水庫(kù)以上流域主要發(fā)生“7·11”“7·20”“10·07”3 場(chǎng)暴雨洪水，前兩者暴雨中心位置均位于中下游地區(qū)，岳城水庫(kù)入庫(kù)主要來(lái)自中下游區(qū)域產(chǎn)流；后者暴雨中心位于清漳河、濁漳河上游，岳城水庫(kù)入庫(kù)主要來(lái)自于上游來(lái)水。

（2）根據(jù)暴雨徑流分析可知，“7·20”暴雨洪水，觀臺(tái)～岳城水庫(kù)區(qū)間為暴雨中心，區(qū)間內(nèi)包含都黨溝等多條支流，產(chǎn)流量較大，根據(jù)觀臺(tái)站實(shí)測(cè)徑流與岳城水庫(kù)反推入庫(kù)徑流計(jì)算得到觀臺(tái)～岳城水庫(kù)區(qū)間降水徑流系數(shù)為0.54，與暴雨洪水調(diào)查中都黨溝流域降水徑流系數(shù)較為接近，認(rèn)為合理。

（3）岳城水庫(kù)反推入庫(kù)徑流扣除水庫(kù)水面接納降水總量、區(qū)間暴雨產(chǎn)流及黃沙疏干水入庫(kù)量，與觀臺(tái)站實(shí)測(cè)資料進(jìn)行月年徑流對(duì)比，可知岳城水庫(kù)反推入庫(kù)徑流量偏大8275 萬(wàn)m3，偏大2.92%，符合有關(guān)規(guī)范要求。

（4）根據(jù)次洪徑流誤差分析可知，岳城水庫(kù)2021年10 個(gè)階段徑流過(guò)程中，當(dāng)岳城水庫(kù)閘門(mén)啟閉頻繁，泄洪流量變動(dòng)快速時(shí)，反推入庫(kù)徑流與觀臺(tái)站實(shí)測(cè)徑流相差較大，當(dāng)岳城水庫(kù)泄洪流量相對(duì)穩(wěn)定時(shí)，二者相差較小。蓄水頂托對(duì)觀臺(tái)斷面水流流態(tài)造成一定影響，降低實(shí)測(cè)流量精度，且閘門(mén)泄洪曲線(xiàn)使用時(shí)間過(guò)長(zhǎng)，多年未修訂，隨著多年上游來(lái)水來(lái)沙及水庫(kù)泄洪洞、溢洪道形態(tài)變化的影響，閘門(mén)泄洪曲線(xiàn)可能存在一定誤差，引起反推入庫(kù)徑流不準(zhǔn)確。建議對(duì)岳城水庫(kù)泄洪流量進(jìn)行比測(cè)，重新修訂閘門(mén)泄洪曲線(xiàn)。