楚萬強，李永麗，謝 龍

(1.黃河水利職業(yè)技術(shù)學院，河南 開封 475003；2.江蘇省交通科學研究院股份有限公司，江蘇 南京 210017)

?

長江干流重慶主城河段造床流量計算

楚萬強1，李永麗1，謝 龍2

(1.黃河水利職業(yè)技術(shù)學院，河南 開封 475003；2.江蘇省交通科學研究院股份有限公司，江蘇 南京 210017)

為反映長江干流重慶主城河段河道特性，確定出造床流量，根據(jù)長江上游水文局1961—2007年的實測資料，利用馬卡維耶夫法和頻率分析法分別對造床流量進行了計算，通過對二者方法及計算結(jié)果的綜合分析，推薦了寸灘河段及玄壇廟河段的造床流量值。通過多元回歸分析，建立了造床流量與水力參數(shù)的經(jīng)驗公式，其結(jié)果可為各河段的航道治理以及同類型河段造床流量的確定提供參考。

航道工程；造床流量；頻率分析

重慶主城河段干流段上起大渡口(距壩里程644.5 km)，下至銅鑼峽(距壩里程598 km)，河段長36 km，流域面積2.85 km2。該河段地處長江上游丘陵寬谷河段,平面形態(tài)呈連續(xù)彎曲狀，航道內(nèi)礁石林立，灘險密布，航道條件極為復(fù)雜[1]。近十年來，隨著西部大開發(fā)戰(zhàn)略的實施，重慶地區(qū)水上運輸需求日益旺盛，為了滿足水運發(fā)展需求，交通運輸部提出了“深下游、暢中游、延上游”的建設(shè)發(fā)展思路[2]，為此需要對目前上游部分灘險河段進行整治，切實提高其航道通過能力。在河道整治規(guī)劃中，造床流量作為主要參數(shù)之一，常用來作為設(shè)計河流常年的斷面和平面形態(tài)的依據(jù)。因此筆者主要研究該河段造床流量的特性及計算方法，為治理此類河段提供理論基礎(chǔ)。

1 造床流量計算方法分析

造床流量Qf是造床強度與多年流量過程綜合造床強度相等的單一流量[3]，它是河道演變中最重要的自變量，決定了河流的平均形態(tài)，對河床的塑造起著重要作用。由于造床強度是一個很模糊的概念，造床流量的計算很多，目前代表性的有以下幾種方法。

1.1 馬卡維耶夫法

前蘇聯(lián)學者馬卡維耶夫認為河床的變化與造床強度及造床歷時有關(guān)[4]。由于床面變化是河段輸沙作用的結(jié)果，因此可用輸沙率來表征造床強度；且造床歷時可用該級強度對應(yīng)的流量所占流量過程的頻率確定[5]，因此將造床作用定義為：

G=QmPJ

(1)

式中：Q為流量；m為變指數(shù)，由Q～Qs的關(guān)系確定(Qs為輸沙率)；J為水面比降；P為本級流量的頻率密度，P=ΔP/ΔQ，ΔQ為流量分級，ΔP為ΔQ對應(yīng)的頻率。

1.2 頻率分析法

錢寧[6]提出：“河床的演變不僅取決于上游來水來沙的絕對數(shù)量，還與它們的過程有關(guān)?！痹趯﹂L江上游順直河段多年數(shù)據(jù)統(tǒng)計、分析的基礎(chǔ)上，錢寧認為造床流量接近于4年(P= 25%) 一遇的洪峰流量。

1.3 平灘流量法

韓其為[4]分析認為：河道水流從枯水流量逐漸增大的過程中，流速和床面切應(yīng)力不斷增大，造床作用增強，而水流漫灘后，水流分散且邊灘阻力較大，主槽造床作用減弱，因此平灘流量法把河道水位與河漫灘齊平時水位對應(yīng)的流量作為造床流量。

由于此法適用于河槽寬淺、河漫灘寬闊(一般洪水漫灘長度可達到主槽的5～10倍)的沖擊河段，而重慶河段屬大比降山區(qū)河流，河槽斷面主要呈“V”字形峽谷，灘唇高程與水位-流量關(guān)系不易準確標定，因此筆者不考慮此計算方法。

2 長江干流重慶主城河段造床流量計算

筆者根據(jù)長江上游水文局對長江干流寸灘、玄壇廟水文站1961—2007年的實測資料，對各方法分別進行計算。

2.1 馬卡維耶夫法

2.1.1 計算步驟

將實測的日均流量匯總，找出流量值的最大值與最小值，以500 m3/s為步長劃分流量區(qū)間，確定每一流量區(qū)間內(nèi)的概率P，并計算每一區(qū)間內(nèi)流量的加權(quán)平均值Q*；建立流量Q與坡降J的曲線關(guān)系，查找每一級平均流量Q*所對應(yīng)的比降；建立水流輸沙率Qs與流量Q的指數(shù)關(guān)系，確定變指數(shù)m的取值；然后根據(jù)式(1)計算每一級流量的造床作用，擬合G-Q的關(guān)系曲線，曲線峰值(即G的最大值)對應(yīng)的流量即為造床流量。

2.1.2Q-J關(guān)系確定

天然河流中，由于水面比降難以施測，可根據(jù)該河段附近的兩個水文站，其實測的水位差值除以距離即為水面坡降。寸灘河段的比降采用寸灘站和銅鑼峽站的實施數(shù)據(jù)，玄壇廟河段的比降采用玄壇廟站和鵝公巖站的實測數(shù)據(jù)。

通過實測數(shù)據(jù)擬合各站流量-水面比降關(guān)系，其曲線關(guān)系如圖1。

圖1 Q-J曲線關(guān)系 Fig.1 Q-J curve relation

2.1.3Qs-Q關(guān)系確定

由于泥沙運動的能量來源于水流的切應(yīng)力[7]，因此馬卡維耶夫認為輸沙率Qs可用流量Q和比降J來表征，通常情況下輸沙率隨水面比降變化較小，因此可忽略水面比降變化對其影響。分別建立各河段Qs-Q的關(guān)系(圖2)，可確定寸灘河段變指數(shù)m=2.316，玄壇廟河段的m=2.612。

馬卡維耶夫推薦沖積平原河流m=2，對比可知山區(qū)大比降河流輸沙率受流量影響要大于沖積平原河流。

圖2 Qs-Q曲線關(guān)系Fig.2 Qs-Q curve relation

2.1.4G-Q曲線擬合

根據(jù)式(1)分別計算各河段每級流量的造床作用G，點繪成曲線，如圖3。從圖中可知寸灘河段的造床流量為16 500～17 000 m3/s，玄壇廟河段的造床流量為14 500～16 000 m3/s。對比可知，在重慶主城河段，支流入?yún)R后(寸灘河段)的造床作用要略大于入?yún)R前(玄壇廟河段)。

圖3 G-Q曲線關(guān)系Fig.3 G-Q curve relation

2.2 頻率分析法

根據(jù)洪水經(jīng)驗頻率計算公式〔式(2)〕[8]，采用皮爾遜Ⅲ型水文頻率分布曲線[9]，繪制各河段的洪峰流量經(jīng)驗頻率曲線(圖4)，對應(yīng)可查出寸灘站的4年(P=25%) 一遇的洪峰流量為17 500 m3/s，即為造床流量，玄壇廟站的造床流量為13 800 m3/s：

(2)

式中：P(i)為樣本中流量≥i的頻率；m為樣本中流量≥i的次數(shù)；n為樣本總?cè)萘俊?/p>

圖4 P-Q曲線關(guān)系Fig.4 P-Q curve relation

3 計算結(jié)果分析

由于頻率分析法是在對長江上游山區(qū)河流的統(tǒng)計基礎(chǔ)上建立的經(jīng)驗公式，與實際值會有一定的誤差，且大比降山區(qū)河流不同河段地形、邊界條件差異較大，該經(jīng)驗公式在應(yīng)用到具體河段時也會有一定的誤差。馬卡維耶夫法則較全面的從造床強度和造床歷時兩方面考慮了水流對河床的切應(yīng)力作用[10]，通過具體河段的實測數(shù)據(jù)進行計算，誤差相對較小。兩種方法的計算結(jié)果如表1。

表1 長江干流重慶主城各河段造床流量計算

由表1數(shù)據(jù)看出，這兩種計算方法所得結(jié)果較為相近，據(jù)分析，綜合來看用馬卡維耶夫法所得結(jié)果誤差較小，因此推薦在參考馬氏法計算結(jié)果基礎(chǔ)上，寸灘河段的造床流量為17 000 m3/s，玄壇廟河段的造床流量為15 000 m3/s，前者(下游)略大于后者(上游)，這主要是因為兩河段之間有嘉陵江支流入?yún)R，流量增大所致。

4 長江干流重慶主城河段造床流量特性分析

河段的沖淤變化是水流與泥沙相互作用的結(jié)果，由河道的來水來沙情況、河床邊界條件、河段輸沙能力等因素決定。對具體河段來說，輸水輸沙是緊密相連的，因此可用某一特征流量來表征河段長期水沙條件變化，該特征流量即為造床流量。另外，水流對河床的造床作用與上游來水情況相關(guān)，反映河段來水情況的水力參數(shù)主要有多年平均流量、多年汛期平均流量等，這些特征參數(shù)與造床流量必定存在一定的聯(lián)系。因此可將造床流量定義為：

(3)

式中：Q年為多年平均流量；Q汛為多年汛期平均流量；a，b，c，d，e均為待定系數(shù)。

根據(jù)實測資料，采用麥夸特迭代法對式(3)進行多元回歸分析，可得寸灘河段的造床流量：

Qf寸灘=0.56Q年0.86+0.33Q汛0.53+172.2

(4)

玄壇廟河段的造床流量：

Qf玄壇廟=0.77Q年0.71+0.41Q汛0.60+23.62

(5)

5 結(jié) 論

1)筆者針對大比降山區(qū)河段，根據(jù)1961—2007年的實測資料，分別利用馬卡維耶夫法和頻率分析法對長江干流寸灘和玄壇廟河段的造床流量進行了計算，通過兩種方法及其計算結(jié)果的對比推薦寸灘河段的造床流量為17 000 m3/s，玄壇廟河段的造床流量為15 000 m3/s。為大比降山區(qū)河流造床流量的確定提供了參考。

2)用造床流量這一特征流量來表征河段長期水沙變化，用水力參數(shù)來表征河段多年來水情況，分析這兩者的相關(guān)性，并用多元線性回歸確定待定系數(shù)，建立長江干流重慶主城河段造床流量的經(jīng)驗公式，可為科學預(yù)測該河段的造床流量提供參考。

[1] 謝龍.三峽變動回水區(qū)末端段復(fù)合水動力條件分析及對泥沙輸移的影響[D].重慶:重慶交通大學,2013. Xie Long.Study on Compound Hydrodynamic Conditions in terminal of Fluctuating Backwater Reach of Three Gorge Reservoir and Its Influence of Sediment Transport [D].Chongqing:Chongqing Jiaotong University,2013.

[2] 謝龍,楊勝發(fā),付旭輝,等.天白水電站溢洪道彎曲段流態(tài)優(yōu)化試驗[J].重慶交通大學學報:自然科學版,2013,32(2):310-312. Xie Long,Yang Shengfa,Fu Xuhui,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.

[3] 趙志舟,呂娜.長江上游急彎分汊河段通航整治汊道選擇[J].水運工程,2009(10):112-117. Zhao Zhizhou,Lv Na.Regulation channel choice for the pebble shoals of branching reach with sharp bend shoal on the upper reach of the Yangtze River [J].Port & Waterway Engineering,2009(10):112-117.

[4] 韓其為.黃河下游輸沙及沖淤的若干規(guī)律[J].泥沙研究,2004(3):1-13. Han Qiwei.Some rules of sediment transportation and deposition-scouring in the lower Yellow River [J].Journal of Sediment Research,2004(3):1-13.

[5] Colosimo C,Copertino V A,Vehri M.Friction factor evaluation on gravel-bed rivers [J].Journal of Hydraulic Engineering,1988,114(8):861-869.

[6] 錢寧.河床演變學[M].北京:科學出版社,1989. Qian Ning.River Bed Evolution Learning [M].Beijing:Science Press,1989.

[7] 張紅武,張俊華,卜海磊,等.試論推移質(zhì)輸沙率公式[J].南水北調(diào)與水利科技,2011,9(6):140-145. Zhang Hongwu,Zhang Junhua,Bu Hailei,et al.Discussion of bed-load transport equations [J].South-to-North Water Diversion and Water Science & Technology,2011,9(6):140-145.

[8] 詹道江,徐向陽,陳元芳.工程水文學[M].4版.北京:中國水利水電出版社,2010. Zhan Daojiang,Xu Xiangyang,Chen Yuanfang.Engineering Hydrology [M].4th ed.Beijing:China WaterPower Press,2010.

[9] Wilcock P R.Critical shear stress of natural sediments [J].Journal of Hydraulic Engineering,1993,119(4):491-505.

[10] 李雙江,謝龍.長江上游塘土壩河段水流特性及河床演變分析[J].重慶交通大學學報:自然科學版,2013,32(4):673-676. Li Shuangjiang,Xie Long.Riverbed evolution and water flow characteristics of Tangtuba Reach in the upper stream of Yangtze Rive [J].Journal of Chongqing Jiaotong University:Natural Science,2013,32(4):673-676.

Calculation of Dominant Discharge of Chongqing Urban Reach in Yangtze River Main Stream

Chu Wanqiang1, Li Yongli1, Xie Long2

(1. Yellow River Conservancy Technical Institute, Kaifeng 475003, Henan, China; 2. Jiangsu Transportation Institute, Nanjing 210017, Jiangsu, China)

To reflect the channel characteristics of Chongqing urban reach of the Yangtze River and to determine the discharge of forming the river bed, the dominant discharge was calculated by Makovieve method and frequency analysis method, which was based on the measured data of the Yangtze River committee hydrographic office from 1961 to 2007. Through the comprehensive analysis on the calculation results of the above two method, the bed-building discharge of Cuntan reach and Xuantanmiao reach was recommended. Finally, the empirical formula between the bed-building discharge and the hydraulic parameters was established through the multiple regression analysis, and its results could offer the reference for the channel governance of reaches and the determination of the dominant discharge of the reaches with the same kind.

waterway engineering; dominant discharge; frequency analysis

10.3969/j.issn.1674-0696.2015.02.15

2013-01-05；

2013-12-05

國家自然科學基金項目(50579041)

楚萬強(1976—)，男，河南洛陽人，講師，碩士，主要從事水利水電工程方面的研究。E-mail: 287109894@qq.com。

TV143

A

1674-0696(2015)02-069-03