羅全勝，謝 龍

(1. 小流域水利河南省高校工程技術(shù)研究中心，河南開封 475003； 2. 黃河水利職業(yè)技術(shù)學(xué)院，河南開封 475003； 3. 江蘇省交通科學(xué)研究院股份有限公司，江蘇南京 210017)

三峽水庫175 m方案蓄水后，常年回水區(qū)的末端位于重慶長壽，變動回水區(qū)末端到達重慶江津[1]，整個重慶主城河段都成為變動回水區(qū)，河道特性在庫區(qū)特性及天然河道情況下轉(zhuǎn)換，水流條件極其復(fù)雜，水流特性發(fā)生顯著變化[2]。同時，龍門浩河段位于重慶主城區(qū)域，河段內(nèi)散布著大片卵石淺灘，常因航深、航寬不足而礙航[3]。據(jù)不完全統(tǒng)計，自1957年來，龍門浩河段整治次數(shù)多達6次，但整治效果不好且不易穩(wěn)定，該河段目前基本每年都要進行挖槽疏浚及相關(guān)航道維護。究其原因，主要是因為水動力條件過于復(fù)雜，同時對卵石成灘機理認識不足。因此，針對三峽水庫175 m蓄水后引起龍門浩河段的水動力條件變化，分析其水流特性，將為該河段河床演變的研究、航道治理提供依據(jù)，此項工作具有十分重要的意義。

圖1 龍門浩河段河勢Fig.1 River regime of Longmenghao reach

1 河段自然條件

龍門浩河段上起黃桷渡，下至玄壇廟，全長4.8 km，河段平面形態(tài)微彎，彎曲半徑約為1 800 m，上接九龍灘2個連續(xù)急彎，下連長江與嘉陵江的交匯口，河段內(nèi)老鸛磧、雞翅膀伏于河心，右岸有豬兒石向航槽凸伸(圖1)，航道條件十分惡劣，常因航深、航寬不足而影響通航，基本每年都需要進行疏浚維護，是重慶河段的主要灘險和重點整治河段[4-5]。

根據(jù)豬兒磧河段上游7.6 km處鵝公巖水位站實測數(shù)據(jù)[6]，成庫前2003—2007年，龍門浩河段日均徑流量為7 057 m3/s，日均輸沙量6.94 t/s，成庫后2007—2010年則分別為6 813 m3/s與5.13 t/s，來水量基本不變，泥沙輸移量則明顯下降，降幅達35.3%。

2 河段二維數(shù)模的建立及驗證

2.1 模型建立

綜合考慮模型計算精度、計算結(jié)果認可度以及計算條件的適用性[7-8]，本文采用Aquaveo.SMS(地表水系統(tǒng)模擬軟件)進行模擬計算。選用2003年及2010年實測地形數(shù)據(jù)，分別建立龍門浩河段成庫前后的平面水流二維數(shù)學(xué)模型。計算模塊采用有限單元法(FESWMS)，全河段采用等邊三角網(wǎng)格，為保證精度，網(wǎng)格間距設(shè)為10 m，全河段共有125 664個節(jié)點及76 808個網(wǎng)格。

2.2 模型驗證

根據(jù)蓄水前后的多次實測資料對各時段、各流量級下模型計算精度進行驗證，從模擬計算結(jié)果來看，水位值誤差在0.1 m以內(nèi)的保證率為93.2%，流速值誤差在10%以內(nèi)的保證率為79.8%，可見本文建立的龍門浩河段二維數(shù)學(xué)模型計算精度較高，適用于該河段水流條件的模擬計算。

3 成庫前后水流特性變化分析

將三峽175 m蓄水以后壩前水位調(diào)度過程點繪于圖2。從圖2可見，三峽水庫年內(nèi)調(diào)度可分為4段：① 6月初至9月末為汛期，上游來流量猛增，該時段內(nèi)按防洪限制水位運行，壩前水位降至全年最低，回水末端位于長壽，龍門浩河段恢復(fù)天然河道特性；② 10月初至10月末為蓄水期，此時上游洪峰已過，來流量逐漸減小，庫區(qū)水位受調(diào)度影響迅速上抬，回水末端向上游移動，龍門浩河段又轉(zhuǎn)變?yōu)閹靺^(qū)特性；③ 11月初至次年1月末為蓄水維持期，此階段上游來流量變化甚微，壩前水位維持在175 m高水位運行，壩前至回水末端河段的水流特性基本維持不變；④ 2月初至5月末為消落期，上游來水降至最小，河段進入枯水期，為保證下游通航及供水需求，壩前水位逐漸消落，直至5月末降至消落低水位。

圖2 175 m試驗性蓄水以來壩前水位調(diào)度過程Fig.2 Scheduling process of upstream water level after 175 m experimental water storage

選擇成庫后年內(nèi)變化各階段內(nèi)的平均流量代表該階段的水流特性，參照長江上游水文局的實測資料，將各特征流量下的對應(yīng)信息列于表1。

表1 成庫前后龍門浩河段計算條件Tab.1 Calculation conditions of Longmenghao reach before and after impoundment

考慮壩前調(diào)度具有周期性，因此在選擇計算條件時，為避免時間因子的影響，在流量誤差相差5%的情況下，盡量保證成庫前與成庫后對應(yīng)年內(nèi)時間一致。

3.1 水位變化

根據(jù)表1條件下的計算成果，將各工況下龍門浩河段成庫前后的水位沿程變化點繪于圖3。

圖3 成庫前后龍門浩河段水位變化Fig.3 Changes in water level along Longmenghao reach before and after impoundment

從圖3可見，成庫后汛期水位與成庫前一致，其他時期則大幅提升，提升幅度從幾米至十幾米不等，蓄水維持期水位增幅最大，消落期增幅最小。越靠近下游，水位抬升越明顯。同時，龍門浩河段同一斷面在6 800 m3/s下蓄水前后的水位差，小于3 680 m3/s、卻大于3 350 m3/s對應(yīng)的值，可見成庫后水位抬升幅度主要取決于蓄水程度，即壩前水位，受流量影響較小。

3.2 橫斷面流速變化

選擇CY20斷面作為代表斷面，將各工況下龍門浩河段橫向流速變化繪于圖4。

由圖4可見，成庫后非汛期內(nèi)，同一流量下龍門浩河段水面線變寬，流速減小，斷面流速分布更加均勻，橫向流速狀況較成庫前發(fā)生大幅變化。水動力軸線(主流區(qū))有向左岸擺動的趨勢。

進一步統(tǒng)計可知，蓄水維持期該斷面的平均流速降幅最大，達到89.3%，遠大于蓄水期的55.4%及消落期的60.1%。流速變化最大的區(qū)域集中在航槽附近(斷面平距340～490 m)，在蓄水維持期，局部位置流速降幅最大可達95.3%。

(a) 汛 期 (b) 蓄水期

(c) 蓄水維持期 (d) 消落期圖4 龍門浩河段成庫前后橫斷面流速變化 Fig.4 Changes in cross-section velocity along Longmenghao reach before and after impoundment

3.3 水面比降變化

根據(jù)表2計算結(jié)果，成庫前后龍門浩河段的水面比降有以下特點：① 成庫前，河段比降與河段來流量呈指數(shù)關(guān)系，流量越大河段比降越??；成庫后，河段比降主要取決于尾水(壩前水位)壅高程度，受來流量影響較小，比降與流量對應(yīng)關(guān)系散亂。② 天然情況下，河段比降較大，基本在2以上，在受回水影響下，比降大幅下降，各時期的平均比降均小于1。在同流量下，成庫前的比降通常要比成庫后大幾倍，甚至十幾倍。③ 各時期內(nèi)平均比降變幅差異較大，蓄水維持期降幅最大，達到93.61%，蓄水期及消落期降幅分別為81.9%與72.4%，汛期則基本無變化，這與水位、流速等變化相對應(yīng)。

3.4 水流功率變化

根據(jù)文獻[1]，長江干流重慶河段水流功率可用下式表示：

(1)

圖5 龍門浩河段成庫前后水流功率變化Fig.5 Changes in stream power of Longmenghao reach before and after impoundment

式中：U為斷面平均流速(m/s)；H為斷面平均水深(m)；J為斷面平均水面比降；g為重力加速度(9.81 m/s2)，γS為推移質(zhì)密度(kg/ m3)；γ為水的密度(kg/ m3)；D為推移質(zhì)中值粒徑(m)。

將各時期水力因子平均量代入式(1)，計算得出水流功率，并點繪于圖5。從圖中可看出，成庫后非汛期內(nèi)，水流功率大幅下降，其中，各時期的變化幅度與水位、流速等因子相對應(yīng)。根據(jù)各時期持續(xù)時間進行加權(quán)估算，成庫后年內(nèi)平均水流功率較成庫前下降47.9%，水流功率的年內(nèi)變化分布極不均勻，這將導(dǎo)致河段輸沙規(guī)律的大幅變化，甚至可能引起新的沖淤變化及改變河床演變規(guī)律。

4 結(jié) 語

本文采用SMS平面水流二維數(shù)學(xué)模型，對三峽大壩175 m方案成庫運行后龍門浩河段水流特性變化進行了分析研究，分析結(jié)果表明：

(1)通過成庫前后年內(nèi)壩前水位變化對比，可將三峽水庫年內(nèi)調(diào)度分為4個時期：汛期(6月初至9月末)、蓄水期(10月初至10月末)、蓄水維持期(11月初至次年1月末)和消落期(2月初至5月末)。

(2)成庫后非汛期內(nèi)，龍門浩河段轉(zhuǎn)變?yōu)閹靺^(qū)特性，尾水位抬升，水深增大，過水?dāng)嗝嬖黾?，水面比降放緩，水流流速、水流功率大幅下降，水流動力條件發(fā)生重大變化。汛期內(nèi)河段恢復(fù)天然河道特性，水流特性與成庫前基本一致。

(3)取各時期的平均流量下水流特性代表各時期水力因子的平均狀況，通過成庫前后各時期相同流量下該河段水流特性模擬對比，蓄水維持期水力因子的變幅最大、蓄水期次之、消落期最小、汛期則不發(fā)生變化。

(4)從各斷面橫向流速變化來看，成庫后整個斷面過流區(qū)域拓寬，流速分布趨于均勻，變化最大的區(qū)域主要分布在航槽附近。同時，斷面的水動力軸線有向左岸擺動的趨勢。

(5)通過加權(quán)估算，成庫后年內(nèi)平均水流功率較成庫前下降47.9%，同時由于水流功率的年內(nèi)變化分布極不均勻，將引起原有的輸沙平衡發(fā)生新的變化。

參 考 文 獻：

[1]謝龍. 三峽變動回水區(qū)末端段復(fù)合水動力條件分析及對泥沙輸移的影響[D]. 重慶： 重慶交通大學(xué)， 2013. (XIE Long. A study of compound hydrodynamic conditions in terminal of fluctuating backwater reach of the Three Gorges Reservoir and its influence on sediment transport[D]. Chongqing： Chongqing Jiaotong University， 2013. (in Chinese))

[2]李雙江， 謝龍. 長江上游塘土壩河段水流特性及河床演變分析[J]. 重慶交通大學(xué)學(xué)報： 自然科學(xué)版， 2013， 32(4)： 673-676. (LI Shuang-jiang， XIE Long. Riverbed evolution and water flow characteristics of Tangtuba reach in the upper stream of the Yangtze River[J]. Journal of Chongqing Jiaotong University(Natural Science)， 2013， 32(4)： 673-676. (in Chinese))

[3]HEY R D. Flow resistance in gravel-bed rivers[J]. Journal of the Hydraulics Division， ASCE， 1979， 105(4)： 365-379.

[4]GOODWIN P. Analytical solutions for estimating effective discharge[J]. Journal of Hydraulic Engineering， 2004， 130(8)： 729-738.

[5]倪晉仁， 惠遇甲， 張國生. 嘉陵江入?yún)R對重慶河段水力特征影響的水力學(xué)分析[J]. 泥沙研究， 1991(2)： 29-37. (NI Jin-ren， HUI Yu-jia， ZHANG Guo-sheng. Hydraulic analyses of the effect of inflow of Jialingjiang River on flow characteristics of Chongqing Reach[J]. Journal of Sediment Research， 1991(2)： 29-37. (in Chinese))

[6]郭繼明， 陳子湘. 三峽水庫變動回水區(qū)銅鑼峽河段演變研究[J]. 長江科學(xué)院院報， 1998(3)： 1-5. (GUO Ji-ming， CHEN Zi-xiang. Evolution research of Tongluoxia river reach in changing backwater region of TGP reservoir[J]. Journal of Yangtze River Scientific Research Institute， 1998(3)： 1-5. (in Chinese))

[7]謝龍， 楊勝發(fā)， 付旭輝， 等. 天白水電站溢洪道彎曲段流態(tài)優(yōu)化試驗[J]. 重慶交通大學(xué)學(xué)報： 自然科學(xué)版， 2013， 32(2)： 310-312. (XIE Long， YANG Sheng-fa， FU Xu-hui， et al. Flow pattern optimization experiment of bending section of Tianbai hydroelectric station spillway[J]. Journal of Chongqing Jiaotong University(Natural Science)， 2013， 32(2)： 310-312. (in Chinese))

[8]ACKERS P， WHITE W R. Sediment transport： new approach and analysis[J]. Journal of the Hydraulics Division， ASCE， 1973， 99(11)： 41-60.