朱利晴，俞俊平，吳江鵬，何金寶，孫斌

(1.江西省交通科學(xué)研究院， 江西 南昌 330200；2.江西省港航管理局， 江西 南昌 330006)

三峽水庫(kù)正常蓄水后，長(zhǎng)壽到洛磧河段的回水變動(dòng)區(qū)航道條件變得極為復(fù)雜，雖然國(guó)內(nèi)有許多學(xué)者[1-6]從河床演變理論分析、數(shù)學(xué)模型計(jì)算、實(shí)體模型試驗(yàn)等不同角度對(duì)三峽蓄水后變動(dòng)回水區(qū)的泥沙淤積情況進(jìn)行了研究，得到了較為深入的認(rèn)識(shí)，但以往大多采用19世紀(jì)60年代水沙資料，然而近年來(lái)，長(zhǎng)江上游及其支流建庫(kù)、水土保持等使得水沙條件發(fā)生了明顯變化，該水沙資料已難以體現(xiàn)輸沙量減小的最新趨勢(shì)[7]，原有的研究成果可能與實(shí)際情況存在很大差異，不再適用。本文采用二維水沙長(zhǎng)時(shí)段數(shù)學(xué)模型，選用2002年—2011年+1998年的10+1序列對(duì)三峽庫(kù)區(qū)回水變動(dòng)區(qū)水流與泥沙運(yùn)動(dòng)規(guī)律進(jìn)行了模擬研究，通過(guò)對(duì)回水變動(dòng)區(qū)各時(shí)段水流運(yùn)動(dòng)規(guī)律及泥沙運(yùn)動(dòng)形態(tài)的研究，主要包括懸移質(zhì)、推移質(zhì)的運(yùn)動(dòng)和水流流速、流態(tài)等的運(yùn)動(dòng)規(guī)律，進(jìn)而分析該河段河道沖淤長(zhǎng)時(shí)段變化規(guī)律，從而為航道整治及日常航道維護(hù)管理提供科學(xué)依據(jù)與理論基礎(chǔ)。

1 三峽庫(kù)區(qū)泥沙沖淤特性

長(zhǎng)壽-洛磧河段位于重慶下游約50 km，河段全長(zhǎng)約28.5 km，出口距三峽樞紐約530 km。該河段具有典型的山區(qū)河流特點(diǎn)，包括彎道，邊灘，淺灘、深槽、分叉等多種河型，計(jì)算成果可反映多種河型的泥沙沖淤特點(diǎn)，具有較強(qiáng)的代表性，該河段的河勢(shì)圖見(jiàn)圖1，其中SC1～SC9是選取的9個(gè)典型斷面。

圖1 三峽工程回水變動(dòng)區(qū)長(zhǎng)壽到洛磧河段河勢(shì)圖

1.1 沖淤?gòu)?qiáng)度

相關(guān)資料顯示，三峽水庫(kù)蓄水以來(lái)，2003年3月至2011年11月間庫(kù)區(qū)大壩至銅鑼?shí){干流全河段泥沙淤積總量為12.649億m3，其中大壩至李渡鎮(zhèn)段淤積12.417億m3，李渡鎮(zhèn)至銅鑼?shí){河段淤積2326萬(wàn)m3。通過(guò)三峽庫(kù)區(qū)各河段淤積量和淤積強(qiáng)度對(duì)比可知，泥沙淤積分布和淤積強(qiáng)度與河道形態(tài)存在著密切的關(guān)系，河道寬谷段淤積強(qiáng)度相對(duì)較大，窄深段淤積強(qiáng)度較小甚至局部出現(xiàn)沖刷現(xiàn)象，蓄水后大壩～銅鑼?shí){干流寬谷段總淤積量為11.106億m3[8]，占全河段總淤積量的93.1%，窄深河段總淤積量為0.826億m3，僅占全河段總淤積量的6.9%。

1.2 沖淤形態(tài)

（1）淤積現(xiàn)象普遍。蓄水前，三峽庫(kù)區(qū)縱剖面呈鋸齒狀分布。蓄水后泥沙淤積使縱剖面發(fā)生了一定變化。相較于2006年，2011年深泓線(xiàn)高程大，即工程河段存在普遍游積現(xiàn)象。

（2）垂向沖淤變化顯著。三峽庫(kù)區(qū)兩岸一般為基巖組成，故岸線(xiàn)基本穩(wěn)定，斷面變化主要表現(xiàn)在河床的垂向沖淤變化。自蓄水以來(lái)，三峽庫(kù)區(qū)淤積形態(tài)主要有[9]：

①主槽平淤，此淤積方式分布于庫(kù)區(qū)各河段內(nèi)，如壩前段、臭鹽磧河段、黃花城河段等；

②沿濕周淤積，此淤積方式也分布于庫(kù)區(qū)各河段內(nèi)；

③左岸或右岸主槽淤積，此淤積形態(tài)主要出現(xiàn)在某些河道形態(tài)為彎道處，以土腦子河道為典型；沖刷形態(tài)主要表現(xiàn)為主槽沖刷和沿濕周沖刷，一般出現(xiàn)在河道水面較窄的峽谷段和回水末端位置。

工程所處的洛磧河段位于銅鑼?shí){至涪陵河段中段部位，為川江典型的淺灘河段，三峽水庫(kù)進(jìn)入156 m蓄水期后，該河段成為三峽水庫(kù)變動(dòng)回水區(qū)河段，試驗(yàn)性蓄水后該河段枯水位進(jìn)一步抬高，但由于目前該河段受蓄水影響時(shí)間較短，且處于變動(dòng)回區(qū)上段，汛期基本為天然河道狀態(tài)，該河段還未出現(xiàn)明顯的累積性淤積狀態(tài)。

（3）洲灘深槽淤積穩(wěn)定。根據(jù)相關(guān)實(shí)測(cè)河床地形資料，工程河段整體河勢(shì)穩(wěn)定，沒(méi)有出現(xiàn)大的灘槽移位。相較2006年，2009年河床地形普遍淤積1.5 m～2 m，局部深槽邊灘淤積超過(guò)2 m，如麻柳場(chǎng)淤積約3 m，中擋壩淤積約4 m，雀石子淤積約3 m，板凳角淤積3 m，川江駁船廠(chǎng)碼頭前淤積4 m，重鋼碼頭淤積6 m，王家灘右岸邊灘淤積4 m～6 m，龍舌梁淤積3 m；2010年地形較2009年普遍淤積1 m～2 m，過(guò)年石至麻柳場(chǎng)出現(xiàn)2 m淤積帶，中擋壩淤積2 m～4 m，車(chē)家石梁化?？诖a頭、川江駁船廠(chǎng)和重鋼碼頭等水工建筑物附近均出現(xiàn)了 2 m～5 m的淤積，騎馬橋出現(xiàn)2 m～4 m淤積，龍舌梁至羊角堡左岸邊灘有2 m～5 m淤積。

2 二維水沙模型建立及驗(yàn)證

2.1 二維水沙模型

本文所用二維水沙模型通過(guò)參考熊偉[9]等人的研究?jī)?nèi)容建立，利用 Mike軟件對(duì)所建立的二維水沙模型進(jìn)行相關(guān)數(shù)值模擬，因篇幅所限且此處并非本文研究重點(diǎn)，模型建立的基本方法、邊界條件、應(yīng)注意的核心問(wèn)題及模型求解等詳見(jiàn)參考文獻(xiàn)[10]。本文采用三角形非結(jié)構(gòu)網(wǎng)格對(duì)計(jì)算域進(jìn)行網(wǎng)格剖分，計(jì)算區(qū)域從南坪壩（上游航道里程609 km）開(kāi)始，至長(zhǎng)壽釣魚(yú)嘴（上游航道里程580.5 km）處結(jié)束，全長(zhǎng)共計(jì)28.5 km，建筑物附近5 m、灘段10 m～20 m、其它河段40 m～50 m，平均20 m，按疏密漸變的方式劃分，共剖分節(jié)點(diǎn)82 958個(gè)，單元168 046個(gè)。

2.2 糙率研究

河道對(duì)水流的阻力用糙率值表示，在二維水沙模型的水力計(jì)算中，糙率值是關(guān)鍵參數(shù)，糙率的大小與河床質(zhì)組成、河道坡降、平面地形、斷面形狀及流量、水位等參數(shù)相關(guān)，需要用實(shí)測(cè)資料進(jìn)行率定[11]。參考長(zhǎng)江重慶航運(yùn)工程勘察設(shè)計(jì)院“三峽庫(kù)區(qū)（175 m運(yùn)用初期）設(shè)計(jì)最低通航水位計(jì)算與分析”，采用一維數(shù)學(xué)模型對(duì)工程河段的糙率進(jìn)行率定，長(zhǎng)壽-扇沱、扇沱-麻柳嘴和麻柳嘴-太洪崗各河段平均糙率分別為0.0366、0.0333和0.0394。

2.3 敏感性分析

一般影響泥沙淤積量的因素主要有床沙粒徑d50、來(lái)沙條件Q沙和輸沙因子（kb是推移質(zhì)輸沙因子，ks是懸移質(zhì)輸沙因子）。

首先，為了比較不同床沙粒徑對(duì)淤積的影響，固定來(lái)水來(lái)沙條件和輸沙因子kb=ks=1，選取了 4種床沙粒徑進(jìn)行分析，分別為d50=0.28 mm，0.2 mm，0.1 mm和0.05 mm，進(jìn)行了2009年的淤積計(jì)算。圖2是進(jìn)口流量、含沙量和出口水位的邊界條件；圖3比較了2009年4種粒徑的淤積量，可見(jiàn)床沙粒徑0.28 mm和0.2 mm淤積量基本相同，床沙粒徑0.1 mm和0.05 mm淤厚基本相當(dāng)。隨著粒徑增加，淤積量增加。分析認(rèn)為不同粒徑對(duì)淤積厚度有一定的影響，但敏感性不強(qiáng)。

其次，為了比較不同來(lái)沙條件對(duì)淤積的影響，固定床沙粒徑d50=0.2 mm，kb=ks=1，比較了2種來(lái)沙條件的淤積厚度。來(lái)沙量減小，淤積隨之減小，但影響并不顯著。

圖2 2009年水位、流量和含沙量過(guò)程

圖3 不同床沙粒徑淤厚比較

最后，為了比較輸沙因子對(duì)河床輸沙的影響，分別取值計(jì)算2009年6月1日至6月15日共15 d的淤積情況。各組淤積高程比較可知，不同輸沙因子對(duì)淤厚有明顯影響。

2.4 水位驗(yàn)證

驗(yàn)證資料參考值采用 2008年麻柳嘴和扇沱水位站實(shí)測(cè)水位值。通過(guò)模型計(jì)算值與實(shí)測(cè)水位值的對(duì)比分析可知：計(jì)算值與實(shí)測(cè)值整體契合程度良好，水位值基本一致，不過(guò)由于汛期水位通常變化較大，故汛期部分時(shí)段計(jì)算值與實(shí)測(cè)值稍有出入，但總體而言誤差尚在可接受范圍內(nèi)。

2.5 地形驗(yàn)證

本文采用2006年－2009年三峽水庫(kù)回水變動(dòng)區(qū)實(shí)測(cè)地形資料進(jìn)行地形驗(yàn)證?？傮w來(lái)看，數(shù)模計(jì)算與實(shí)測(cè)主要淤積部位基本一致，實(shí)測(cè)平均淤高1.16 m，數(shù)值計(jì)算平均淤高1.06 m，數(shù)模計(jì)算淤積量為3863.28萬(wàn)m3，實(shí)測(cè)地形淤積3846.72萬(wàn)m3，兩者相差16.57萬(wàn)m3，相對(duì)誤差為0.43%。

3 河床沖淤規(guī)律分析

3.1 邊界條件

為了分析工程河段長(zhǎng)時(shí)段的河床沖淤變化規(guī)律，選用2002年—2011年和1998年的10+1序列重復(fù)一次對(duì)22 a間工程河段的沖淤變化計(jì)算，該序列既考慮輸沙量減小的最新趨勢(shì)，也包含各種特征年，考慮因素比較全面，具有較強(qiáng)的代表性。

3.2 淤積量變化

三峽水庫(kù)運(yùn)行5 a、10 a、15 a和22 a淤積量分別為 3867.878 萬(wàn) m3、9774.527 萬(wàn) m3、15 846.700 萬(wàn)m3和24 144.410萬(wàn)m3，對(duì)應(yīng)的累積平均淤積厚度分別為1.22 m、3.09 m、5.01 m和7.63 m，由于每年的來(lái)沙量不同，每年的淤積強(qiáng)度和淤積厚度有所差異，最大的淤積強(qiáng)度發(fā)生在 2021年，淤積強(qiáng)度為1749.350 m3/a，淤厚為0.55 m/a；多年平均淤積強(qiáng)度為1097.473 m3/a，多年平均淤厚為0.35 m/a。

3.3 沖淤斷面變化

總體來(lái)看，淤積主要發(fā)生在深槽和邊灘，當(dāng)三峽工程運(yùn)行 22 a時(shí)，河段內(nèi)的深槽基本上都被淤平，而淺灘河床淤積則相對(duì)較小，致使斷面形狀發(fā)生很大變化，寬深比明顯增大，逐漸向“U”形發(fā)展，河槽向“單一”發(fā)展。具體來(lái)看，SC1斷面（過(guò)年石附近）屬于回流區(qū)，其右岸深槽最大淤積厚度達(dá)25 m，在深槽左側(cè)淺區(qū)有少量沖刷；SC2斷面（金錢(qián)罐附近）整個(gè)斷面均發(fā)生淤積，最大淤厚約10 m；SC3斷面（石船梁附近）屬于寬淺河段，其深槽部位最大淤厚約 13 m；SC5斷面（渝懷鐵路橋上游600 m）整體斷面淤積，最大淤厚約8.5 m；SC7斷面（重鋼碼頭下游400 m）屬于彎沱區(qū)，其左深槽有26 m淤積，右岸邊灘淤厚約5 m。

3.4 沖淤平面分布

三峽水庫(kù)運(yùn)行5 a、10 a、15 a和22 a后，可以看出整個(gè)河段均出現(xiàn)累積性泥沙淤積，這是由于汛后蓄水，水位大幅度提升，水流不能歸槽，汛期淤積的泥沙難以得到有效沖刷造成的，同時(shí)還可以看出泥沙主要淤積在回流區(qū)、彎沱及河道的寬淺河段，主要的淤積區(qū)域位于南坪壩左槽、南坪壩尾至麻柳場(chǎng)、雀石子至點(diǎn)燈石左槽一帶、渝懷鐵路橋下游至車(chē)家石梁、白鷺嘴、騎馬橋至螺絲口左槽、龍舌梁至羊角堡左槽，這些區(qū)域在三峽水庫(kù)運(yùn)行5 a、10 a、15 a和22 a時(shí)，淤高分別為5 m～7 m、6 m～26 m、10 m～32 m和14 m～34 m。

4 泥沙淤積對(duì)重點(diǎn)河段航道的影響

由于泥沙主要淤積在回流區(qū)、彎沱及河道的寬淺河段，故這些河段航道條件的變化相比其他河段更大。對(duì)于回流區(qū)，受累積性淤積影響較大的是過(guò)年石處河段，由于汛期水流趨直，該處回流區(qū)擴(kuò)大，泥沙淤積，而在枯水期由于水庫(kù)蓄水使得沖沙時(shí)間減少，泥沙不能完全沖走，造成了累積性淤積，但由于泥沙主要淤積于深槽，以及河道水位的整體抬升，其航道條件不會(huì)發(fā)生較大變化；對(duì)于彎沱區(qū)，騎馬橋至螺絲口河段航道條件變化較為顯著，兩岸彎沱會(huì)出現(xiàn)較大強(qiáng)度的淤積，隨著淤積年限的延長(zhǎng)，彎沱逐漸被淤平，主槽將趨于順直，航道寬度會(huì)有所減小，尤其是王家灘段，在中枯水時(shí)，其河道被忠水磧分成兩汊，左汊為枯水航道，而左汊通航水域又主要在深槽，三峽水庫(kù)運(yùn)行5 a時(shí)，河段會(huì)因彎沱及深槽的淤積，出現(xiàn)河深及河寬不足的問(wèn)題，需加強(qiáng)關(guān)注，必要時(shí)采取適當(dāng)?shù)恼未胧?，但是隨著河道水位的整體抬升，運(yùn)行20 a后河段整體水位將抬高1 m～6 m，且兩汊深槽逐漸被淤滿(mǎn)，斷面逐漸向“U”形發(fā)展，成為單汊河段，存在的問(wèn)題將會(huì)得以解決；而對(duì)于寬淺河段，由于過(guò)流面積突然增加，流速迅速減小而發(fā)生累積性泥沙淤積，淤積主要發(fā)生在邊灘和江心洲處，在邊灘處淤積，較為典型的是石船梁處河段，其凸岸進(jìn)一步發(fā)展，河型逐漸向凹岸深槽的彎曲型過(guò)渡，寬度雖有所減少，但仍滿(mǎn)足航運(yùn)寬度要求，而在江心洲處淤積，變化較大的是南坪壩處河段，右汊河槽雖有所淤積，但其水深不會(huì)發(fā)生明顯變化，寬度變化也不大，均滿(mǎn)足航運(yùn)要求，左汊則有逐漸淤?gòu)U的趨勢(shì)。

總的來(lái)說(shuō)，除王家灘段外的其他河段，雖然也有累積性淤積，但航道條件不會(huì)產(chǎn)生明顯變化，而對(duì)王家灘段，短期內(nèi)可能會(huì)出現(xiàn)礙航現(xiàn)象，但長(zhǎng)期來(lái)看，因泥沙落淤部位主要在深槽和邊灘，不會(huì)對(duì)航道水深產(chǎn)生不利影響，且河道整體水位抬升，該段的礙航現(xiàn)象可得到緩解，航道條件也不會(huì)發(fā)生明顯變化。

5 結(jié)論

本文通過(guò)建立二維水沙數(shù)學(xué)模型，在分析其敏感因素及其可靠性的基礎(chǔ)上，對(duì)工程河段進(jìn)行了長(zhǎng)時(shí)段水沙數(shù)值模擬，分析了工程河段河床沖淤變化規(guī)律，預(yù)測(cè)了2012年—2033年間累積淤積情況，表明工程河段有比較明顯的累積性淤積趨勢(shì)，泥沙主要淤積在回流區(qū)、彎沱及河道的寬淺河段，短期內(nèi)這些河段可能會(huì)出現(xiàn)礙航現(xiàn)象，但因泥沙落淤部位主要在深槽和邊灘，且河道整體水位抬升，航道條件不會(huì)產(chǎn)生明顯變化。本研究成果與實(shí)際情況十分吻合，可為計(jì)算河段的航道規(guī)劃、航道維護(hù)、航道整治、碼頭選址以及挖砂采石等提供指導(dǎo)，也可為航道管理提供支撐。