三峽回水變動區長壽-洛磧河段水沙數值模擬

畢 方 全

(長江重慶航道局，重慶401147)

三峽回水變動區長壽-洛磧河段水沙數值模擬

畢 方 全

(長江重慶航道局，重慶401147)

三峽庫區泥沙累積性淤積給航道的日常維護管理以及航道整治帶來了新的問題。建立了二維水沙數學模型對三峽庫區回水變動區長壽至洛磧河段的泥沙輸移進行了數值模擬。敏感性分析表明輸沙因子對泥沙淤積量影響最大。水位、流速流向和地形驗證表明：數模計算成果較好地反映了實際水文、泥沙淤積情況。計算表明：受三峽回水影響，長壽至洛磧河段發生累積性的泥沙淤積，淤積區主要位于岸線凹凸不平的彎沱、回流區及河道的寬淺河段，淤積部位主要是深槽和邊灘。淤積時間在每年的汛期6—9月份。受泥沙淤積影響，上、下洛磧和王家灘在低水位時航深不足，有礙航現象，建議加強監測，及時采取整治措施。

航道工程；水沙數學模型；三峽工程；回水變動區；河床演變

三峽水庫正常蓄水運行后，長壽-洛磧河段成為三峽工程變動回水區，是長江上游重點礙航河段之一。據有關學者研究[1-2]，變動回水區河段泥沙淤積較為嚴重，河道兩岸形成淤積邊灘，影響港口作業，部分航道條件惡化，為航道的日常維護管理以及航道整治帶來了新的問題。水沙條件和水庫調度方式是影響三峽工程變動回水區河床演變的重要因素。長江上游大型梯級水庫群及嘉陵江航電樞紐的建設改變了三峽工程變動回水區的來水來沙條件[3-6]，庫區淤積情況與已有的研究成果相比發生了較大的變化。因此對泥沙問題的治理，還需做進一步的包括蓄水后的原型觀測[7]、數值模擬、河工模型試驗等方面的深入研究。筆者建立了二維水沙數學模型，在新的水沙條件下，對長壽-洛磧河段的沖淤趨勢進行分析。

1 河段概況

長壽-洛磧河段位于重慶下游約50 km，長約30.5 km，出口距三峽大壩約530 km，是三峽工程變動回水區中寬淺、多灘的典型河段之一。天然情況下長壽-洛磧河段基本遵循年內洪淤枯沖、年際間沖淤交替且基本平衡的沖淤特性。

三峽水庫蓄水運行當壩前水位低于148 m時，長壽-洛磧河段恢復為天然河段；當壩前水位148～158 m時，長壽-洛磧河段部分受壅水影響；當壩前水位在158～175 m時，整個長壽-洛磧河段處于壅水狀態。

2 二維水沙模型建立及驗證

2.1 二維水沙模型建立

建立的二維水沙模型基本方程包括水流連續方程、動量方程、懸移質不平衡輸沙方程、推移質不平衡輸移方程及河床變形方程，有關模型的基本控制方程、數值解法、關鍵問題環節的處理、初始邊界條件及動邊界技術等詳見文獻[8]。計算河段進口位于南坪壩(上游航道里程609 km)，下游出口位于長壽釣魚嘴(上游航道里程580.5 km)，全長28.5 km。采用三角形非結構網格對計算域進行網格剖分，建筑物附近5 m、灘段10～20 m、其它河段40～50 m，平均20 m，按疏密漸變的方式劃分，共剖分節點82 958個，單元168 046個。

2.2 糙率研究

糙率值是水力計算的關鍵參數之一，它綜合反映了河道對水流的阻力。采用一維數學模型對工程河段的糙率進行率定，長壽-扇沱、扇沱-麻柳嘴和麻柳嘴-太洪崗各河段平均糙率分別為0.036 6，0.033 3和0.039 4。

2.3 敏感性分析

一般影響泥沙淤積量的因素主要有床沙粒徑d50、來沙條件Q沙和輸沙因子(Kb是推移質輸沙因子，Ks是懸移質輸沙因子)。目前常用的輸沙公式有Engelund and Hansen模型、Van-Rijn模型、Engelund and Fredse模型和Meyer-Peter and Müller模型。為校證各輸沙公式計算結果與實測輸沙量之間的誤差，在各輸沙公式中引入輸沙因子，其取值范圍為0.1～2.0。通過對4種泥沙粒徑、兩個來沙條件和6組輸沙因子組合的淤積厚度進行比較。限于篇幅，僅列出不同輸沙因子淤積厚度比較，如圖1。分析認為輸沙因子對淤積強度的影響最為突出。使用Van-Rijn模型，選取不同輸沙因子對2006—2009年河床地形進行驗證計算，確定Kb=Ks=0.15；床徑d50=0.123 mm；來沙條件Q沙采用寸灘水文站實測來沙量。

圖1 不同輸沙因子取值淤積厚度比較(2009年)Fig.1 Deposition comparison of different sediment load factor

2.4 水位驗證

圖2比較了扇沱和麻柳嘴水位站計算水位和實測水位，可見計算值與實測值相位符合良好，水位值基本相當。誤差較大的時段主要出現在汛期，分析原因是汛期水位變幅較大，數模只給定平均糙率不能完全反應不同水位糙率的影響。

圖2 水位驗證(2008年)Fig.2 Verification of water level

2.5 流速流向驗證

圖3比較了計算流速與浮標實測值的大小和方向，可見流速分布及最大值、最小值的位置均與實測資料較為一致。

圖3 流速流向驗證Fig.3 Verification of water level

2.6 地形驗證

采用2006—2009年實測地形資料進行地形驗證。圖4是選取的典型斷面淤積地形比較，總體來看數模計算與實測主要淤積部位基本一致，實測平均淤高1.16 m，數值計算平均淤高1.06m，數模計算淤積量為3 863.28萬m3，實測地形淤積3 846.72萬m3，兩者相差16.57萬m3，相對誤差在0.43%。

圖4 地形驗證Fig.4 Verification of topography

3 河床沖淤規律分析

3.1 邊界條件

為了分析工程河段淤積規律，進行了2009—2011年這3年的模擬計算。圖5顯示了2009—2011年進口流量過程線、進口含沙量過程線和出口水位過程線。

圖5 邊界條件Fig.5 Boundary conditions

3.2 沖淤平面分布

圖6顯示了2009年初至2011年末河床淤積厚度分布。整體來看工程河段河床平均淤厚1.15 m。局部來看，在南坪壩左槽發生淤積，平均淤厚約3～4 m；南坪壩后至金錢罐右槽有淤積帶，平均淤厚2～3 m，過年石與下迎春石之間最大淤厚達6 m；小磧包、小金堆、大金堆和荷葉坪一帶淤厚約 2～3 m；沙公溪、大灘嘴至占燈石原深槽一帶有比較明顯的淤積，平均淤厚在3～5 m，最大淤厚達7 m；渝懷鐵路下游深潭部位出現約7 m的淤積?；Ｆ反a頭前沿有2～3 m的淤積。周家磧右槽深潭淤厚3～5 m；白鷺嘴和重鋼碼頭上游淤積明顯，最大淤厚7～9 m；騎馬橋、大沙壩和螺絲口左槽帶有較大的淤積帶，平均淤厚約5～6 m，騎馬橋處最大淤厚達9 m；忠水磧右槽有1～2 m的淤積；龍舌梁至羊角灘左深槽有約5～6 m的淤積，最大淤厚達9 m?？偟膩砜矗俜e區主要位于岸線凹凸不平的彎沱(扇沱～忠水磧)、回流區(南坪壩后)及河道的寬淺河段，淤積后的河床岸線趨于平順，邊灘淤高長大，河道變窄，主槽位置一般不變。

圖6 工程河段2009—2011年淤積平面分布Fig.6 Siltation plane distribution in works river section from Changshou to Luoqi during 2009 to 2011

3.3 沖淤斷面變化

圖7顯示了8個典型斷面淤積分布?？傮w來看，淤積主要發生在深槽和邊灘。具體來看，SC1斷面(過年石附近)右岸深槽淤積約4 m，右岸邊灘淤積約2 m，在深槽左側淺區有少量沖刷；SC2斷面(金錢罐附近)主要在邊灘淤積，最大淤厚約6 m，斷面形狀向U形發展；SC3斷面(石船梁附近)主要在深槽有少量淤積，淤厚約3 m；SC4斷面(薛家梁附近)在左岸深槽有大量淤積，淤厚近8 m；SC5斷面(渝懷鐵路橋上游600 m)淤積發生在兩岸邊灘，淤厚約4 m；SC6斷面(化?？诖a頭附件)左側深槽和右岸邊灘，最大淤積約4 m；SC7斷面(重鋼碼頭下游400 m)左深槽有8 m淤積；SC8斷面(橫石板)在忠水磧兩側深槽有淤積，淤厚約 5 m。

圖7 工程河段典型斷面淤積分布Fig.7 Typitcal sectional siltation distribution in works river section

3.4 淤積量變化

圖8顯示了2009—2011年底累積淤積量，可見工程河段有比較明顯累積性淤積的趨勢。圖9是平均每月淤積量比較，可見淤積主要發生在每年的6—9月的汛期。分析淤積量隨時間的變化規律為汛期來水來沙量增大，以落淤為主；汛后壩前水位由145 m 逐步抬升至175 m，水位升高，流速減小，挾沙能力減弱，泥沙淤積；每年的汛前壩前水位從175 m降至防洪限制水位145 m，河床主槽及邊灘略有沖刷。因此，淤積量過程每年均有沖淤交替的變化，總體上呈上升趨勢。

圖8 2009—2011年累積淤積量Fig.8 Cumulative siltation volume from 2009 to 2011

圖9 2009—2011年月均淤積量Fig.9 Average monthly siltation volume from 2009 to 2011

4 泥沙淤積對重點淺灘航道的影響

洛磧和王家灘是川江航道上最著名的低水位期彎道淺灘。洛磧河段平面形態為寬窄相間，窄段水深較大，寬段多江心洲和邊灘，在洲尾和邊灘，或者兩個邊灘之間形成過渡段淺區。上洛磧河段右岸南坪壩尾的狹長淤積帶迫使低水位航槽逐漸向河心移動，若遇到大沙小水年份時，在次年消落期可能出現礙航問題，建議加強監測，及時采取整治措施。下洛磧河段左邊灘的淤積使得低水位航槽向右岸位移，航線貼近五金堆、金錢罐等礁石區，對船舶安全航行有較大的危害，需加強關注。王家灘段中低水位時河道被忠水磧分成兩汊，左汊為低水位航道。由于泥沙淤積，左槽航深不足，有礙航現象。

5 結 語

建立了二維水沙數學模型對三峽庫區回水變動區長壽至洛磧河段的泥沙輸移進行了數值模擬，敏感性分析認為輸沙因子取值對泥沙淤積量影響最大。水位、流速流向和地形驗證表明數模計算成果較好反映實際水文、泥沙淤積情況。

對長壽-洛磧河段2009—2011年泥沙淤積的計算表明，該河段發生累積性的泥沙淤積，淤積區主要位于岸線凹凸不平的彎沱(扇沱～忠水磧)、回流區(南坪壩后)及河道的寬淺河段，淤積部位主要是深槽和邊灘。淤積后的河床岸線趨于平順，邊灘淤高長大，河道變窄，主槽位置一般不變。淤積發生時間在每年的汛期6—9月份。受泥沙淤積影響，上、下洛磧和王家灘在低水位航深不足，有礙航現象，建議加強監測，及時采取整治措施。

[1] 袁美琦,張秀芹.三峽工程變動回水區施工期和175m蓄水初期洛磧至王家灘河段演變規律及治理措施分析[C]// 國務院三峽工程建設委員會辦公室泥沙專家組.三峽工程泥沙問題研究“九五”成果匯編第五卷.北京:知識產權出版社,2002. YUAN Meiqi,ZHANG Xiuqin.Analysis of river bed evolution and control measures about Luoqi to Wangjiatan Reach in fluctuating backwater area of Three Gorges Project in period of construction and 175m initial stage[C]// The Sediment Expert Group of the Three Gorges Project Construction Committee of the State Council.The"NineFive"ResearchResultsAssemblyontheSedimentProbleminThreeGorgesProject.Vol.5.Beijing: Intellectual Property Publishing House Co.,Ltd.,2002.

[2] 竇國仁,萬聲淦,陸長石.三峽工程變動回水區泥沙淤積的試驗研究[J].水利水運科學研究,1995(4):327-335. DOU Guoren,WAN Shenggan,LU Changshi.Investigation on sedimentation in varying backwater zone of Three Gorges Project[J].JournalofNanjingHydraulicResearchInstitute,1995(4):327-335.

[3] 李義天,陳建,鄧金運,等.上游建庫對三峽水庫變動回水區航道的影響[J].泥沙研究,2008(5):31-37. LI Yitian,CHEN Jian,DENG Jinyun,et al.Influence of reservoir construction upstream on the navigation channel in the fluctuating backwater area of the Three Gorges Reservoir[J].JournalofSedimentResearch,2008(5): 31-37.

[4] 陳建,李義天,鄧金運,等.水沙條件變化對三峽水庫泥沙淤積的影響[J].水力發電學報,2008,27(2):97-102. CHEN Jian,LI Yitian,DENG Jinyun,et al.Influence on deposition of the Three Gorges Reservoir caused by the change of water silt conditions[J].JournalofHydroelectricEngineering,2008,27(2): 97-102

[5] 陸永軍,左利欽,季榮耀,等.水沙調節后三峽工程變動回水區泥沙沖淤變化[J].水科學進展,2009,20(3):318-324. LU Yongjun,ZUO Liqin,JI Rongyao,et al.Changes of sediment deposition and erosion at Chongqing reach in backwater area of Three Gorges Project after reservoir adjusting of the upstream in the Yangtze River[J].AdvancesinWaterScience,2009,20(3) :318-324.

[6] 張緒進,何進朝,母德偉.上游來水來沙變化及對重慶河段泥沙淤積的影響[J].水利水運工程學報,2010(1):23-29. ZHANG Xujin,HE Jinchao,MU Dewei.Impact of incoming runoff and sediment on sediment deposition along the Chongqing reach[J].Hydro-ScienceandEngineering,2010(1):23-29.

[7] 許全喜.2011年度長江三峽工程水文泥沙原型觀測分析報告[R].武漢:長江水利委員會水文局,2012. XU Quanxi.AReportontheHydrologicalandSedimentPrototypeObservationoftheYangtzeRiverThreeGorgesProjectin2011[R].Wuhan: Hydrology Bureau of Yangtze River Conservancy Committee,2012.

[8] 熊偉.非結構網格下三維潮流泥沙數值模擬[D].大連:大連理工大學,2012. XIONG Wei.3DNumericalSimulationofTidalCurrentsandSedimentUsingUnstructuredMeshModel[D].Dalian: Dalian University of Technology,2012.

Sediment Numerical Simulation of Changshou to Luoqi Reach in FluctuatingBackwater Area of Three Gorges Project

BI Fangquan

(Changjiang Chongqing Waterway Bureau, Chongqing 401147, P.C. China)

The accumulated siltation and sediment in Three Gorges reservoir zone caused new problem to the daily management, maintenance and regulation of waterway. In this paper, the 2D flow-sediment mathematical model was established to simulate the muddy sand transport in the Three Gorges Reservoir regional backwater variation zone in the river reach from Changshou to Luoqi. Sensitivity analysis shows that sediment input factor is key parameter influencing sedimentation volume. Water level, flow velocity and direction and topography verifications show that this mathematical model calculation results well reflect the reality hydrology and siltation. The calculation results show that subject to the influence of backwater of the Three Gorges Reservoir, cumulative deposition happens in the reach from Changshou to Luoqi and siltation zones are located primarily in the river sinuses，the backflow zones and broad-shallow reaches of the river and siltation mainly occurred in deep trench and side beach in flooding season from June to September. Because of the siltation and sediment, water depth at upper and lower reaches of Luoqi and Wangjiatan is insufficient so as to impede navigation. It is suggested to enhance monitoring and take timely rectification and regulating measures.

channel engineering; flow-sediment mathematical model; Three Gorges project; backwater variation zone; river bed evolution

10.3969/j.issn.1674-0696.2016.02.15

2015-02-02；

2015-06-01

畢方全(1962—)，男，重慶人，高級工程師，碩士，主要從事航道、港口工程方面的研究。E-mail：993182321@qq.com。

TV147. 5

A

1674-0696(2015)02-064-06