龍啟建，李克鋒，汪青遼

(四川大學(xué) 水力學(xué)與山區(qū)河流開發(fā)保護國家重點實驗室，四川 成都 610065)

電站的日調(diào)節(jié)泄流和大壩泄洪引起的非恒定流將改變下游河道的天然水流條件，對大壩下游航道、港口、通航設(shè)施的正常使用以及船舶的安全航運造成一定的影響，特別是電站突然開啟或關(guān)閉，大壩開閘泄洪等，這種高瞬變流對山區(qū)河流通航水流條件及船舶通航水深的影響非常明顯［1］.因此研究電站調(diào)節(jié)非恒定流對樞紐間河道的通航水流條件產(chǎn)生的影響，對保障船舶的安全航行很有必要.針對電站調(diào)節(jié)非恒定流對通航水流條件的影響，劉新［2］建立了漢江石泉樞紐下游非恒定流一維數(shù)學(xué)模型，計算滿足通航保證率條件下的石泉樞紐調(diào)節(jié)方式;鄭守仁等［3］在研究三峽電站調(diào)峰產(chǎn)生的非恒定流對三峽至葛洲壩兩壩間水流條件影響時，根據(jù)模型試驗結(jié)果得出在兩電站聯(lián)合運行時能滿足萬噸級船隊航運要求;吳曉黎等［4］研究了三峽電站不同調(diào)峰流量下河道水流的航運滿足情況.目前研究多數(shù)是針對樞紐下游，以單個電站日調(diào)節(jié)對下游航道水流的影響為主，考慮梯級水電站聯(lián)合運行、調(diào)節(jié)對通航水流影響的較少.

本文以嘉陵江干流中游段在建的亭子口、蒼溪和沙溪航電樞紐為研究對象，采用數(shù)值模擬方法，研究梯級水電站聯(lián)合調(diào)節(jié)非恒定流對樞紐間河道通航水流條件造成的影響，以及反調(diào)節(jié)樞紐的興建對樞紐間河道通航水流的改善程度.

1 通航水力指標要求

電站調(diào)節(jié)非恒定流的水流條件對船舶航行有著重要的影響，因此必須對船舶所要求的通航水流條件進行分析.評價非恒定流對航運影響的指標主要有水位(航深)、涌浪、流速、流態(tài)及水面比降、水位日變幅和小時變率.根據(jù)嘉陵江航道規(guī)劃及《內(nèi)河航道標準》［5］，嘉陵江主要通航水力技術(shù)指標有:最小水深為1.9 m，最小通航河寬為50.0 m，最大平均流速不宜超過3.0 m/s，最大水面比降不超過2.0‰，最大小時水位變幅不超過1.0 m.

船舶逆流上行與上灘時推力必須克服水流對船舶所產(chǎn)生的水流阻力和水面比降所產(chǎn)生的比降阻力.所謂船舶的允許上灘能力，是指船舶上灘發(fā)揮最大馬力而獲得的對岸航速(0.3～0.5 m/s)時，此位置的水流表面流速和水面比降亦即船舶上灘時允許的最大流速和水面比降.當(dāng)船型、載量和船舶有效推力一定時，在滿足航行阻力小于有效推力條件下，允許流速和比降將會有無數(shù)個組合.計算時，首先擬定幾個不同的流速值，求出其流速阻力，由已知的船舶推力，反推出相應(yīng)的允許比降值，多組流速及比降構(gòu)成的曲線即構(gòu)成船舶的上灘能力曲線.根據(jù)嘉陵江航道規(guī)劃，得出嘉陵江2×500 t級船舶滿載時的上灘能力曲線如圖1所示.

圖1 2×500 t級船舶滿載時的上灘能力曲線Fig.1 Curve of 2×500 t ship ascending rapids with full load

2 非恒定流數(shù)學(xué)模型及其驗證

與船舶航行有關(guān)的電站調(diào)節(jié)非恒定流的水力學(xué)特征參數(shù)，如斷面平均流速、水深等，可采用非恒定流一維偏微分方程組進行求解.非恒定流一維偏微分方程組由連續(xù)方程和動量方程組成.

式中:Z為水位(m);Q為流量(m/s);A為過水?dāng)嗝婷娣e(m);Lq為單位河長的旁側(cè)入流量(m/s);x為沿水流方向的水平距離(m);t為時間(s);g為重力加速度(m/s2);B為斷面河寬(m);α為流速分布不均勻系數(shù);Sf為摩阻坡度.Sf一般根據(jù)曼寧公式可寫為:

式中:K為流量模數(shù)(m3/s);n為糙率;V為斷面平均流速(m/s);R為水力半徑;Q為流量(m3/s).

現(xiàn)行求解上述偏微分方程組最常用的方法是數(shù)值解法.本文采用Preissmann隱格式差分法求解上述方程組.該方法具有計算速度快、穩(wěn)定性強、收斂性好的特點，特別適用于流動變化不太劇烈的長距離、長時段的洪水演進計算，在非恒定流數(shù)值計算中被廣泛使用.

本文的計算范圍為嘉陵江干流在建的3個連續(xù)開發(fā)的亭子口、蒼溪、沙溪航電樞紐所在河段，其中亭子口樞紐位于河段起點，沙溪樞紐位于末端，均按低水頭水電連續(xù)開發(fā)建設(shè).且3座航電樞紐建成后，沙溪電站水庫回水影響至蒼溪電站、蒼溪電站水庫回水影響至亭子口電站.三樞紐間河道全長39.10 km，沿程共布置53個計算斷面，最大間距為1 418 m，最小為445 m，其中上游亭子口壩址為1#斷面，蒼溪壩址為13#斷面，沙溪壩址為53#段面.嘉陵江干流亭子口至沙溪樞紐間河段的綜合糙率值，采用2001年洪、中、枯3級典型流量下亭子口至沙溪河段沿程各斷面的水面線觀測資料反算求得.

為驗證非恒定流數(shù)學(xué)模型及其求解方法的精確性，將該數(shù)學(xué)模型應(yīng)用于甘肅境內(nèi)黃河四龍到龍灣段進行模型驗證.該河段全長110 km，共125個計算斷面將四龍至龍灣劃分為124個河段，最大間距為2 343 m，最小為235 m.以四龍斷面2002年7月22日0時至2002年7月26日0時實測流量過程作為上邊界條件，下游給定龍灣斷面同時段實測的水位過程為下邊界條件.計算結(jié)果與四龍至龍灣段中三灘斷面實測的水位過程比較見圖2.從圖2可見，數(shù)值模擬計算得到的水位與實測值水位過程基本吻合，兩者的最大相對水深誤差都基本控制在5%以內(nèi)，說明該數(shù)學(xué)模型能較精確地模擬該河段非恒定流的傳播規(guī)律和過程.

圖2 三灘計算與實測水位過程比較Fig.2 Comparison between the measured and calculated water-level

3 航電樞紐聯(lián)合運行非恒定流的數(shù)值模擬

為充分研究梯級水電站聯(lián)合調(diào)節(jié)時反調(diào)節(jié)樞紐建設(shè)對上、下游河道航運水流條件的改善程度，針對以下2種情形進行計算并作比較分析:(1)亭子口與沙溪樞紐建成，中間蒼溪樞紐未建成，沙溪庫區(qū)起反調(diào)節(jié)作用;(2)3個梯級均建成，蒼溪庫區(qū)對亭子口起反調(diào)節(jié)作用，沙溪庫區(qū)對蒼溪樞紐起反調(diào)節(jié)作用.通過對這2種情況下非恒定流水力要素變化規(guī)律的計算、分析和比較，從而確定梯級航電樞紐聯(lián)合調(diào)節(jié)非恒定流對樞紐間河道通航水流條件影響的規(guī)律，以及反調(diào)節(jié)樞紐在梯級航電樞紐調(diào)度中對通航水流條件的改善作用.

根據(jù)亭子口水利樞紐的設(shè)計運行方式，亭子口電站在系統(tǒng)中需要承擔(dān)調(diào)峰任務(wù)，因此在蒼溪反調(diào)節(jié)樞紐興建之前，非恒定流計算以亭子口電站的日調(diào)峰流量下泄過程作為上邊界條件，以沙溪電站壩址水位過程作為下邊界條件.蒼溪反調(diào)節(jié)樞紐興建之后非恒定流計算中，亭子口至蒼溪段以亭子口日調(diào)峰流量下泄過程為上邊界條件，以蒼溪電站壩址水位過程為下邊界條件，蒼溪至沙溪段以聯(lián)合調(diào)節(jié)時蒼溪電站流量下泄過程為上邊界條件，沙溪壩址水位過程為下邊界條件.根據(jù)上游來水情況及亭子口電站日調(diào)節(jié)方式，亭子口下泄流量按以下3種情況加以考慮:

(1)當(dāng)上游日均來流量小于電站滿發(fā)時的最大引用流量1 200 m3/s，來流全部用于發(fā)電，此時主要研究亭子口電站日調(diào)節(jié)非恒定流對下游樞紐間河道航運安全的影響;為方便比較梯級水電站聯(lián)合調(diào)節(jié)非恒定流對樞紐間河道通航水流條件的影響，以及反調(diào)節(jié)樞紐的興建對樞紐間河道通航水流的改善程度，在蒼溪電站興建前后亭子口都采用同樣的日調(diào)節(jié)下泄流量過程.根據(jù)梯級電站優(yōu)化調(diào)度原則擬定亭子口5種日調(diào)節(jié)方案，及同時段蒼溪電站壩址水位流量過程.亭子口下泄流量過程線如圖3，方案1～3為電站連續(xù)啟閉單臺機組情況;方案4和5為電站突然啟閉多臺機組情況.根據(jù)亭子口日調(diào)節(jié)的下泄流量過程和蒼溪壩前的水位流量關(guān)系曲線計算得同時段蒼溪電站在聯(lián)合日調(diào)節(jié)過程中下泄流量過程見圖4.

圖3 亭子口電站的典型日調(diào)節(jié)下泄流量過程Fig.3 Discharge hydrographs of each scheme of typical daily regulation of Tingzikou

圖4 亭子口影響下的蒼溪電站典型日調(diào)節(jié)下泄流量過程Fig.4 Discharge hydrographs of each scheme of typical daily regulation of Cangxi under the influence of Tingzikou

圖5 亭子口影響下蒼溪水電站典型日調(diào)節(jié)泄洪流量過程Fig.5 Flood hydrographs of each scheme of typical daily regulation of Cangxi under the influence of Tingzikou

(2)當(dāng)上游來流量大于電站最大發(fā)電流量而小于10年一遇流量11 800 m3/s時，電站滿發(fā)，泄洪閘泄流調(diào)峰，船閘仍繼續(xù)通航，按對航運影響最不利情況進行考慮，取亭子口水電站上游來流量為最大通航流量11 800 m3/s，泄洪閘起始下泄流量為電站滿發(fā)流量1 200 m3/s，此時主要研究閘門開啟歷時分別擬定為5，10，15，20和30 min共5種方案下亭子口下泄流量按線性變化過程時非恒定流對通航水流條件的影響;根據(jù)亭子口泄洪調(diào)節(jié)的下泄流量過程和蒼溪壩前水位流量關(guān)系曲線計算出同時段蒼溪電站在電站聯(lián)合日調(diào)節(jié)時泄洪流量過程見圖5.

(3)當(dāng)上游來流量大于10年一遇流量11 800 m3/s時，機組停發(fā)，船閘停航，閘門全開泄洪，此時水利樞紐接近天然行洪狀態(tài).本文對此情況未作研究.

3.1 電站聯(lián)合發(fā)電調(diào)節(jié)非恒定流對航運的影響

當(dāng)亭子口上游日均來流量小于電站滿發(fā)最大引用流量1 200 m3/s時，蒼溪電站興建前后亭子口至蒼溪電站樞紐間河道，各斷面最大的平均流速和水面比降組合對船舶上灘能力影響見圖6;各發(fā)電方案下通航水力學(xué)指標值的比較見表1.

圖6 蒼溪電站興建前后最大平均流速和最大水面比降組合Fig.6 The maximum mean velocity and water surface slope before and after construction of Cangxi hydropower station

表1 蒼溪電站興建前后通航水力指標比較Tab.1 Navigable hydraulic indexes before and after construction of Cangxi hydropower station

3.2 電站聯(lián)合泄洪調(diào)節(jié)非恒定流對航運的影響

當(dāng)上游來流量大于電站最大發(fā)電流量而小于10年一遇流量11 800 m3/s時，在亭子口各方案閘門啟閉時間的影響下，蒼溪電站興建前后亭子口至蒼溪電站樞紐間河道，各斷面最大的平均流速和水面比降組合對船舶上灘能力影響見圖7，亭子口至蒼溪段各泄洪方案的通航水力學(xué)指標值的對比統(tǒng)計見表2.由于此時啟門歷時均未超過1 h，因此產(chǎn)生的下泄非恒定流量發(fā)生變化的時間較短，為了準確分析和統(tǒng)計非恒定流引起的水力指標變化情況，故采用1 min作為統(tǒng)計時段.

圖7 蒼溪電站興建前后最大平均流速和最大水面比降組合Fig.7 The maximum mean velocity and water surface slope before and after construction of Cangxi hydropower station

表2 蒼溪電站興建前后通航水力指標比較Tab.2 Navigable hydraulic indexes before and after construction of Cangxi hydropower station

4 梯級電站聯(lián)合調(diào)節(jié)時反調(diào)節(jié)樞紐對河道航運的影響

通過對梯級電站聯(lián)合調(diào)節(jié)時反調(diào)節(jié)樞紐興建前后在亭子口發(fā)電、泄洪產(chǎn)生的各方案非恒定流影響下的樞紐間河道水力要素的比較分析，主要有以下結(jié)論:

(1)建蒼溪反調(diào)節(jié)樞紐后，在亭子口發(fā)電、泄洪各方案非恒定流的影響下，河道沿程的最大平均流速和最大水面比降均出現(xiàn)不同程度的減小:在電站發(fā)電調(diào)峰時，最大平均流速由建壩前方案4的3.18 m/s減小至2.81 m/s，最大水面比降由方案4的2.15‰減小至0.98‰;在電站泄洪調(diào)節(jié)時，最大平均流速由方案1的4.35 m/s減小至3.43 m/s，最大水面比降由方案1的3.24‰減小至2.03‰.由此可見，建壩后河道沿程的最大平均流速和水面比降均出現(xiàn)了一定程度的改善.

(2)建反調(diào)節(jié)樞紐后，在亭子口發(fā)電、泄洪各方案非恒定流的影響下，靠近船舶上灘能力曲線的組合點數(shù)量減少明顯，說明亭子口各方案的下泄非恒定流經(jīng)過蒼溪水庫的反調(diào)節(jié)作用后，其流量有所削弱，通航水流變得更平緩，更加有利于船舶通航.建壩前，在亭子口泄洪產(chǎn)生的非恒定流影響下，各方案均有斷面不滿足2×500 t級船舶的上灘要求，而建壩后除啟門歷時為5 min的方案1存在1個斷面不滿足要求外，其他方案在各斷面均能滿足要求，因而反調(diào)節(jié)樞紐興建后泄洪期亭子口電站的閘門啟閉時間不得小于5 min.

(3)建反調(diào)節(jié)樞紐前后，河道沿程的最小水深值變化較為顯著.在亭子口各方案發(fā)電產(chǎn)生的非恒定流影響下，建反調(diào)節(jié)樞紐前沿程最小平均水深均不滿足最小通航水深1.9 m的要求，建反調(diào)節(jié)樞紐后各方案在全河段均能滿足要求;在亭子口各方案泄洪非恒定流影響下，沿程最小平均水深由3.45 m增至3.83 m.

(4)建反調(diào)節(jié)樞紐前后，在亭子口各方案發(fā)電非恒定流影響下，方案4和5的沿程最大小時水位變幅，分別由1.25 m和1.08 m相應(yīng)減小至0.94 m和0.79 m，建壩后，兩方案在全河段均滿足水位變幅要求;在亭子口各方案泄洪非恒定流影響下，沿程最大分鐘水位變幅由方案1的1.06 m減至0.98 m.

(5)在亭子口發(fā)電各方案非恒定流影響下，方案1～3中電站連續(xù)啟閉單臺機組相比方案4和5中電站突然啟閉多臺機組最大水面比降、最大平均流速、最大分鐘水位變幅都相對更小;在亭子口泄洪各方案非恒定流影響下，各方案中隨著閘門啟閉時間的增大最大水面比降、最大平均流速、最大分鐘水位變幅都隨之減小，且方案1～3中不滿足2×500 t級船舶的上灘要求的斷面明顯比方案4和5少.

5 結(jié)語

(1)梯級電站聯(lián)合日調(diào)節(jié)過程閘門突然啟閉及泄洪期間閘門啟閉歷時越短，下泄流量變化率就越大，下游河道受影響的程度和范圍就越大，導(dǎo)致各斷面的最大水面比降和分鐘水位變幅也最大，更不利于船舶上灘.

(2)比較梯級電站聯(lián)合調(diào)節(jié)時反調(diào)節(jié)樞紐的興建前后發(fā)電調(diào)節(jié)及泄洪調(diào)節(jié)期間水流條件，可知反調(diào)節(jié)樞紐的興建對河道通航水流條件有改善作用，更有利于2×500 t級船舶的安全通航.

［1］王志力，陸永軍.向家壩水利樞紐下泄非恒定流的數(shù)值模擬［J］.水利水電科技進展，2008，28(3):12-15.(WANG Zhili，LU Yong-jun.Numerical simulation of unsteady flow in Xiangjiaba hydropower project［J］.Advances in Science and Technology of Water Resources，2008，28(3):12-15.(in Chinese))

［2］劉新.漢江石泉樞紐下游非恒定流數(shù)學(xué)模型研究［J］.水道港口，2005，28(6):425-429.(LIU Xin.Investigation of downstream unsteady flows of Shiquan hydraulic power station with mathematical model［J］.Journal of Waterway and Harbor，2005，28(6):425-429.(in Chinese))

［3］鄭守仁，孫爾雨，楊文俊.三峽電站調(diào)峰與通航問題研究［J］.水電能源科學(xué)，2002，20(2):7-12.(ZHENG Shou-ren，SUN Er-yu，YANG Wen-jun.Study on the navigation and power station for hump modulation of Three Gorges and Gezhouba project［J］.International Journal of Hydroelectric Energy，2002，20(2):7-12.(in Chinese))

［4］吳曉黎，李承軍，張勇傳，等.三峽電站調(diào)峰流量對航運的影響分析［J］.水利水電科技進展，2003，23(6):7-9.(WU Xiao-li，LI Cheng-jun，ZHANG Yong-chuan，et al.Effects of peak regulation at Three Gorges Power Station on navigation［J］.Advances in Science and Technology of Water Resources，2003，23(6):7-9.(in Chinese))

［5］GB 50139-2004，內(nèi)河通航標準［S］.(GB 50139-2004，Navigation standard of inland waterway［S］.(in Chinese))

［6］曹民雄，姜繼紅，唐存本，等.船舶上灘阻力計算方法對比與分析［J］.水利水運工程學(xué)報，2005(2):41-45.(CAO Min-xiong，JIANG Ji-hong，TANG Cun-ben，et al.Comparison and analyses of calculation methods for resistance of ships ascending rapids［J］.Hydro-Science and Engineering，2005(2):41-45.(in Chinese))

［7］趙克玉.天然河道一維非恒定流數(shù)學(xué)模型［J］.水資源與工程學(xué)報，2004，15(1):38-41.(ZHAO Ke-yu.Numerical model of unsteady flow in natural river［J］.Journal of Water Resources＆ Water Engineering，2004，15(1):38-41.(in Chinese))

［8］張防修，王艷平，韓龍喜，等.復(fù)式河道一維洪水演進數(shù)值模擬［J］.水利水電科技進展，2008，28(5):6-9.(ZHANG Fang-xiu，WANG Yan-ping，HAN Long-xi，et al.Numerical simulation of one-dimensional flood routing for a compound channel［J］.Advances in Science and Technology of Water Resources，2008，28(5):6-9.(in Chinese))

［9］吳持恭.水力學(xué)［M］.3版.北京:高等教育出版社，2003.(WU Chi-gong.Hydraulics［M］.3rd ed.Beijing:Higher Education Press，2003.(in Chinese))

［10］金生，王志力，耿艷芳.安康樞紐下游非恒定流數(shù)學(xué)模型研究［J］.水道港口，2005(增刊1):15-19.(JIN Shen，WANG Zhi-li，GENG Yan-fang.Investigation of downstream unsteady flows of Ankang hydraulic power station with mathematical model［J］.Journal of Waterway and Harbor，2005，26(Suppl1):15-19.(in Chinese))