引用格式:.荊江太平口逆流現(xiàn)象成因分析[J].水利水電快報,2025,46(7):11-16.
近年來,荊江三口之一的太平口頻繁出現(xiàn)逆流現(xiàn)象[1]。三峽水庫運行后,荊江三口的分流比減小,虎渡河口門段河道微沖[2-5]。三峽水庫蓄水后虎渡河水量明顯減少[6-8]。事實上,三峽水庫運行后,虎渡河的水沙特性均發(fā)生了變化[9-10],而相關(guān)站點監(jiān)測情況顯示,逆流現(xiàn)象也是在三峽水庫運行之后出現(xiàn)的。目前關(guān)于荊江三口逆流現(xiàn)象的研究大多是對個別年份或場次的分析,為進一步探索太平口逆流現(xiàn)象產(chǎn)生的原因,本文對彌陀寺站所有逆流時段的數(shù)據(jù)進行匯總,采用彌陀寺站首次逆流年份(2015年)至2023年的太平口、彌陀寺等主要控制站點實測水位、流量資料,以及三峽水庫運行期間相關(guān)監(jiān)測站點的水位,對太平口的逆流現(xiàn)象進行分析。
1 研究區(qū)域概況
荊江三口原為荊江四口,由荊江南岸的松滋口、太平口、藕池口、調(diào)弦口組成[11-12],其中調(diào)弦口于 1959年封堵[13],因此現(xiàn)今稱為荊江三口。荊江三口是江湖關(guān)系的紐帶[14],其中太平口原稱“虎渡口\"或“虎口”,為虎渡河分流長江之口門。河口下游約 1.5km 有荊江分洪北閘,當(dāng)荊江分洪工程啟用時,對太平口站水位影響較大。北閘下游約 0.6km 有太平口水位站,下游約 8.6km 有彌陀寺水文站,均為虎渡河控制站。彌陀寺站下游約 41km 為中河口,系虎渡河與松滋河(東支)水系連通的串溝,流向順逆不定,下游約 83km 處有虎渡河節(jié)制閘南閘,再往下游與松滋河及澧水匯合后流入洞庭湖(圖1)。復(fù)雜的河網(wǎng)關(guān)系使得彌陀寺站既受長江干流上游來水影響,又受下游洞庭湖頂托影響。
2逆流現(xiàn)象及成因分析
根據(jù)彌陀寺站歷史實測數(shù)據(jù),太平口逆流首次出現(xiàn)時間為2015年,之后除2016年,每年都會發(fā)生不同時長的逆流現(xiàn)象,其中2023年逆流持續(xù)天數(shù)更是長達(dá)70d 。對彌陀寺站 2015~2023 年水文資料及逆流相關(guān)特征值進行統(tǒng)計,包括逆流持續(xù)天數(shù)、年實測最大逆流流量、最大日平均逆流流量、年徑流量等(表1)。
從持續(xù)時間上看,除2021年與2023年,其余年份逆流累計天數(shù)均小于10d,且2023年逆流累計天數(shù)達(dá)70d 。2018年實測最大逆流流量出現(xiàn)在10月底,即汛期末,其余年份均出現(xiàn)在5\~7月主汛期。結(jié)合年徑流量進行比較發(fā)現(xiàn),2022年與2023年徑流量雖小,但逆流持續(xù)天數(shù)相對較長。綜合分析,逆流持續(xù)時間、逆流流量大小與水文年豐枯沒有明顯關(guān)系。本文將對部分年份出現(xiàn)逆流現(xiàn)象的洪水過程進行分析。
在2015年“東方之星”事件中,三峽水庫壓減下泄流量,長江干流荊江段水位迅速止?jié)q轉(zhuǎn)落,逆流時段持續(xù)了 6d[15] 。逆流期間,三峽水庫入庫流量在10000~16000m3/s 范圍內(nèi)波動,出庫流量從17200m3/s 降至 6800m3/s 。逆流期間三峽水庫調(diào)度過程見圖2。在此期間,太平口水位從 33.53m 降至 31.06m 彌陀寺水位隨之降低,但兩站產(chǎn)生了明顯的負(fù)落差,最大負(fù)落差為 -1.169m 。同時,下游黃山頭(閘上)的水位也與彌陀寺水位變化趨勢一致,三站水位對比如圖3所示。
2017年“長江1號洪水\"期間,為緩解長江中下游防洪壓力,三峽水庫調(diào)度及時控制了長江干流荊江河段及洞庭湖區(qū)水位上漲趨勢[16-17],逆流時段持續(xù)4d。
逆流期間,三峽水庫入庫流量從 26000m3/s 下降至17500m3/s ,維持3d后上升至 30 000m3/s 左右,出庫流量從 27 300m3/s 下降至 8000m3/s 左右,持續(xù)3d 后開始上升。逆流期間三峽水庫調(diào)度過程見圖4。在此期間,太平口水位從 37.08m 降至 33.79m ,彌陀寺水位隨之降低,兩站產(chǎn)生了較小的負(fù)落差,最大負(fù)落差為 -0.246m 。結(jié)合下游黃山頭(閘上)的水位變化趨勢,三站水位對比如圖5所示。
2020年的逆流過程集中于6月 9~16 日,期間三峽水庫開展調(diào)度[18](圖6)。三峽水庫入庫流量先是穩(wěn)定在 12000m3/s 左右,后逐漸增大,最大達(dá)18000m3/s ,出庫流量先由 18000m3/s 下降至 11000m3/s 后面呈周期性波動。相比2015年和2017年,這次的調(diào)度過程具有周期性的特點,因此逆流量相對較小,且持續(xù)時間不長。在此期間,太平口水位從 33.82m 降至 31.97m ,彌陀寺水位隨之降低,兩站產(chǎn)生了負(fù)落差,最大負(fù)落差為 -0.306m ,結(jié)合下游黃山頭(閘上)的水位變化趨勢,三站水位對比如圖7所示。
2021年太平口發(fā)生數(shù)次逆流現(xiàn)象,本文選取5月23~30 日的過程進行分析。5月23~30日三峽水庫調(diào)度過程如圖8所示。可以看出,逆流過程與三峽水庫周期性調(diào)度之間存在很大關(guān)系。三峽水庫入庫流量在 11 000m3/s 左右波動,出庫流量呈周期性波動,平均幅度達(dá) 5000m3/s 左右。在此期間,太平口水位與彌陀寺水位兩站產(chǎn)生的負(fù)落差很小,最大負(fù)落差為-0.026m ,結(jié)合下游黃山頭(閘上)的水位變化趨勢,三站水位對比如圖9所示。
2023年的太平口逆流發(fā)生次數(shù)較多,本文選取6月 10~23 日的過程進行分析。2023年6月10\~23日三峽水庫調(diào)度過程如圖10所示。三峽水庫入庫流量從 10300m3/s 減少至 8000m3/s ,出庫流量呈周期性波動,平均幅度達(dá) 4000m3/s 左右。6月10\~17日,太平口水位從 33.02m 降至 29.84m ,太平口水位與彌陀寺水位兩站產(chǎn)生的負(fù)落差持續(xù)時間長,最大負(fù)落差為 -0.556m ,結(jié)合下游黃山頭(閘上)的水位變化趨勢,三站水位對比如圖11所示。
三峽水庫調(diào)度導(dǎo)致長江干流荊江段水位驟降,虎渡河太平口站與彌陀寺站產(chǎn)生倒比降,為逆流創(chuàng)造了條件。2020年開始,三峽水庫的周期性調(diào)度過程使逆流現(xiàn)象更加頻繁。以黃山頭(閘上)站水位代表下游洞庭湖的頂托作用,長江干流上游來水與下游洞庭湖頂托的影響使彌陀寺站逆流現(xiàn)象更加復(fù)雜。除此之外,松滋河與虎渡河交匯的中河口順逆不定,中河口來水也可能會導(dǎo)致短期的少量逆流產(chǎn)生。
3逆流流量計算
為進一步探索逆流流量大小與已知因素之間的關(guān)系,對 2015~2023 年間彌陀寺站實測水位流量數(shù)據(jù)進行統(tǒng)計,發(fā)現(xiàn)共有51次逆流。根據(jù)水力學(xué)明渠恒定流的計算原理,河道流量大小與比降、河相關(guān)系、河道糙率等因素有關(guān)。其中比降可用太平口與彌陀寺之間的落差來表示,而河相關(guān)系、河道糙率在不同水位級下是不同的,且隨河道沖淤變化也會發(fā)生改變。因此,為了簡化逆流流量的計算,按恒定流考慮,建立彌陀寺站逆流指標(biāo),參照曼寧公式的形式,綜合考慮了水位落差、河相關(guān)系之間的數(shù)值關(guān)系。彌陀寺站逆流的主要原因是太平口站與彌陀寺站之間存在倒比降,根據(jù)太平口至彌陀寺的距離( 6.5km )及兩站之間的水位落差來計算倒比降,再加上彌陀寺站大斷面河相關(guān)系(斷面面積、水面寬、斷面平均水深),將這些參數(shù)整合計算出一個值 F ,為彌陀寺站逆流指標(biāo)。 F 與相關(guān)因素的關(guān)系可用如下公式表示:
式中: A 為彌陀寺站斷面面積, m2;ΔZ 為太平口與彌陀寺的水位落差, m;B 為彌陀寺站水面寬, 為彌陀寺站斷面平均水深,H=A B,m;F為彌陀寺站逆流指標(biāo),m4/3
為探索彌陀寺站逆流流量大小 Q 與彌陀寺站逆流指標(biāo) F 之間的相關(guān)關(guān)系,對 2015~2023 年間彌陀寺站實測數(shù)據(jù)進行統(tǒng)計,共有51次逆流資料。通過對不同水位優(yōu)化分級分別進行擬合,得到 F 與 Q 的相關(guān)關(guān)系(圖12~14)??梢钥闯?,按水位級劃分后,逆流流量隨著參數(shù) F 的增大而增加,水位 (Z) 在 30.22~ 32.00m 時,相關(guān)度達(dá)到了0.9782;水位在 32.00~ 33.00m 時,相關(guān)度為0.9209;水位在 33.00~34.66 m 時,相關(guān)度為0.9200,總體相關(guān)性好。
根據(jù)得到的不同水位級下 F 與 Q 之間的關(guān)系,當(dāng)太平口與彌陀寺實時水位監(jiān)測數(shù)據(jù)出現(xiàn)落差為負(fù)時,則可根據(jù)兩站水位快速計算出相應(yīng)的 F 值,進而根據(jù)F 與 Q 的擬合關(guān)系推算出彌陀寺站此時逆流流量 Q 值。據(jù)統(tǒng)計, 34.66m 是彌陀寺站截至2023年測到的最高逆流水位。實際上,當(dāng)水位較高、來水量較大時,逆流更不容易產(chǎn)生,但具體限制條件還需要后續(xù)研究論證。因此,若后續(xù)還有更高水位對應(yīng)的逆流現(xiàn)象出現(xiàn),可適時再進行分析。
4結(jié)論
(1)太平口逆流產(chǎn)生的原因主要有: ① 三峽水庫調(diào)度導(dǎo)致長江干流荊江段水位驟降; ② 虎渡河受洞庭湖的頂托影響。太平口逆流流量大小與多方面因素有關(guān),在彌陀寺初始水位一致條件下,太平口-彌陀寺的落差越大,逆流流量越大;在太平口-彌陀寺的落差相近時,彌陀寺水位越高,斷面過流面積越大,逆流流量則越大。
(2)根據(jù)統(tǒng)計,目前彌陀寺站出現(xiàn)逆流現(xiàn)象實測最高水位為 34.66m 。當(dāng)長江及虎渡河水位較高時、彌陀寺站逆流不容易產(chǎn)生,可能原因: ① 高水位時,長江水位短時間的波動,傳遞至彌陀寺站時已坦化不明顯; ② 高水位時,長江水位出現(xiàn)持續(xù)驟降的情況少。
若后續(xù)還有更高水位對應(yīng)的逆流現(xiàn)象出現(xiàn),可適時再進行分析。
(3)建立了彌陀寺站不同水位級下彌陀寺站逆流指標(biāo) F 與逆流流量 Q 之間的相關(guān)關(guān)系函數(shù),在太平口與彌陀寺落差為負(fù)的情況下,可初步推算彌陀寺站逆流流量 Q 值,預(yù)判逆流是否發(fā)生以及逆流程度,及時開展水文測驗、提高水文監(jiān)測效率,使水文監(jiān)測方案更加科學(xué)合理準(zhǔn)確,為洞庭湖四口水系綜合整治及區(qū)域水資源調(diào)配提供基礎(chǔ)的數(shù)據(jù)支撐。
參考文獻(xiàn):
[1]馮寶飛,牛文靜,許第兵,等.三峽水庫調(diào)蓄作用下荊江三口逆流現(xiàn)象成因解析[J].人民長江,2019,50(7):16-21.
[2]閆依陽.變化環(huán)境下荊江水文情勢演變特征及其生態(tài)效應(yīng)研究[D].鄭州:華北水利水電大學(xué),2024.
[3]DUAN GL,PENG YB.Flow and sediment measurement forthe Three Gorges Project[J]. River,2024,3(1) :69-77.
[4]魏軒,劉瑜,胡家彬.三峽水庫試驗性蓄水后荊江三口分流變化[J].人民長江,2020,51(8):99-103.
[5]朱玲玲,許全喜,戴明龍.荊江三口分流變化及三峽水庫蓄水影響[J].水科學(xué)進展,2016,27(6):822-831.
[6]蔣婕妤,魯瀚友,趙文剛,等.近70年來虎渡河水沙變化特征與沖淤演變分析[J].泥沙研究,2023,48(5):41-47.
[7]段光磊,彭嚴(yán)波.新水沙條件下荊江河段洲灘演變特征及其影響因素[J].長江科學(xué)院院報,2024,41(7):175-181.
[8]趙秋湘,付湘,孫昭華.三峽水庫運行對荊江三口分流的影響評估[J].長江科學(xué)院院報,2020,37(2):7-14.
[9]方春明,胡春宏,陳緒堅.三峽水庫運用對荊江三口分流及洞庭湖的影響[J].水利學(xué)報,2014,45(1):36-41.
[10]柴澤清,郭小虎,朱勇輝.三峽工程運用后荊江與洞庭湖關(guān)系變化研究進展[J].長江科學(xué)院院報,2023,40(4):17 -23,30.
[11]曹文洪,毛繼新.三峽水庫運用對荊江河道及三口分流影響研究[J].水利水電技術(shù),2015,46(6):67-71,78.
[12]李景保,何霞,楊波,等.長江中游荊南三口斷流時間演變特征及其影響機制[J].自然資源學(xué)報,2016,31(10) :1713-1725.
[13]:秦凱,彭玉明,陳儉煌.荊江三口分流能力變化分析[J].人民長江,2015,46(18):34-38.
[14]許全喜,胡功宇,袁晶.近50 年來荊江三口分流分沙變化研究[J].泥沙研究,2009(5):1-8.
[15]尚全民,褚明華,許弟兵,等.“東方之星”號客輪翻沉事件水利應(yīng)急處置實踐[J].中國防汛抗旱,2016,26(1):34 -37.
[16]龐樹森.三峽水庫對長江中下游防洪補償調(diào)度的作用——以2017年長江第1號洪水為例[J].中國防汛抗
旱,2019,29(8):16-19.
[17] 李強,王琴,黃火林.長江上游梯級水庫群應(yīng)急調(diào)度下荊江河段特殊水情分析[J].水利水電快報,2018,39(10):20-23.
[18]魏林云,李強,魏軒,等.長江中游荊江河段2020 年與1998年洪水對比分析[J].水利水電快報,2023,44(10):8-13.
(編輯:江文)
Causes analysis of backflow phenomenon at Taipingkou of Jingjiang River
YE Jiawen1,HUANG Liemin1,F(xiàn)ANG Zhi1,LYU Shiyang2 (1.JingiangBureauofHdrologyndWaterResocesSureyureufHydroogfangangWateresoucesComisiog zhou43400,China;2.SubordinateBranchofJingzhouChangjiang RiverManagementBureau,Jingzhou4340o,China)
Abstract:Toinvestigatethecausesand influencing factors ofbackflowphenomenonatTaipingkouof Jingjiang River,we used monitoring data from thecontrol stationsof Taipingkou and Mituosi Station along the Hudu River,and considered waterlevel variation processs upstreamand downstream.A functionalfting approach was adopted to calculate the backflow discharge.The results showed that the main causesof the backflow phenomenon were a sudden dropin water level in the Jingjiang section of Changjiang River main stream and the backwater efectof Dongting Lake.The operation of the Three Gorges Reservoir caused a sudden drop in water levels inthe middle reaches of the Jingjiang section.This resulted in areversed waterlevel gradient between the Taipingkou Station andthe Mitaoshi Station of the Hudu River,creating conditions forthebackflow.TherewasacorrlationbetweentheflowrateatMituosi Stationunderdiferent waterlevels and factors such asthe river-stage relationshipatthe station and the waterlevel drop from Taipingkou to Mituosi Station.Thebackflow dischargecan be preliminarilycalculated by backflow index.Theresearch findings can provide fundamental data supportforthecomprehensive improvementof thefour-outlet watersysteminDongting Lakeandregional water resources allocation.
Key words: backflow phenomenon ; Jingjiang River; Taipingkou ; Mituosi Station; Dongting Lake