過飽和總?cè)芙鈿怏w釋放過程預(yù)測

黃奉斌,李 然,鄧 云,曲 璐

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

水電工程在發(fā)揮其巨大經(jīng)濟(jì)和社會效益的同時,不可避免地會對環(huán)境產(chǎn)生不利影響。例如,在大壩泄洪時,伴隨泄水產(chǎn)生的過飽和總?cè)芙鈿怏w(TDG)導(dǎo)致水體中含有過量的氣體,易引發(fā)魚類氣泡病,嚴(yán)重威脅它們的生存和繁殖。為了魚類以及其他水生生物的生存和繁衍,大壩泄水產(chǎn)生的過飽和TDG應(yīng)當(dāng)?shù)玫接行У目刂啤?986年美國環(huán)保局[1]規(guī)定河流允許的TDG過飽和度為110%。我國目前雖然沒有相應(yīng)的環(huán)境控制標(biāo)準(zhǔn),但水體溶解氣體過飽和問題將日漸突出并備受關(guān)注。1987年,我國葛洲壩泄洪時下游水體中大量魚類出現(xiàn)氣泡病癥狀,并導(dǎo)致部分死亡。1994年6月,新安江水庫開閘泄流,導(dǎo)致距電站大壩3km的建德市虹鱒場的虹鱒普遍患?xì)馀莶?因病致死的虹鱒占死魚總量的62.2%[2],這些魚類的死亡均與水中TDG過飽和有關(guān)。

為了定量預(yù)測和研究高壩下游TDG過飽和的影響程度和范圍,需要對過飽和TDG在下游河道的輸移釋放規(guī)律進(jìn)行研究。美國陸軍工程兵團(tuán)[3]在2005年提出了縱向一維過飽和TDG釋放模型。Politano等[4]、Weber等[5]用兩相流模擬了溢洪道下游過飽和氣體的釋放過程。國內(nèi)隨著大量的高壩大庫進(jìn)入設(shè)計和建設(shè)階段,也展開了過飽和TDG的相關(guān)研究,李然等[6]、蔣亮等[7-8]進(jìn)行了過飽和TDG釋放的室內(nèi)試驗以及理論分析研究,對二灘、大渡河等電站進(jìn)行了原型觀測和數(shù)值模擬。程香菊等[9]采用數(shù)值模擬手段對葛洲壩電站下游過飽和TDG的沿程釋放規(guī)律進(jìn)行了初步探討。但總體來說,由于下游河道內(nèi)過飽和TDG原型觀測資料缺乏,對河道內(nèi)過飽和TDG釋放過程的預(yù)測尚存在諸多困難和問題。

本文從室內(nèi)試驗、原型觀測結(jié)果出發(fā),探討得到過飽和TDG沿程釋放規(guī)律,并修正了過飽和TDG的縱向一維釋放模型,該模型可用于定量預(yù)測過飽和TDG的釋放過程。

1 室內(nèi)機(jī)理試驗研究

1.1 試驗裝置和試驗條件

試驗裝置如圖1所示。利用空氣壓縮機(jī)制造溶解氣體飽和度約為116%的過飽和水體,容器尺寸為圓柱形,底面半徑15cm,高30cm,水深 25 cm。水體內(nèi)TDG飽和度測量采用加拿大Point Four Systems Inc.生產(chǎn)的 PT4 Tracker with TGP Smart Probe,將測得的水體中TDG的壓力與當(dāng)?shù)貥?biāo)準(zhǔn)大氣壓進(jìn)行比較,得到水體中TDG的過飽和度,測量精度達(dá)到0.1%。將探頭和攪拌器固定于水面下15 cm處?？刂茢嚢杵鞯牟煌D(zhuǎn)速,在轉(zhuǎn)速 n分別為 200 r/min,400r/min,600r/min時測量TDG飽和度隨時間的變化過程,所得結(jié)果如圖2所示。

圖1 室內(nèi)耗散機(jī)理試驗裝置

圖2 TDG飽和度隨轉(zhuǎn)速及時間變化關(guān)系曲線

由圖2可以看出,在初始時段水體過飽和TDG迅速下降,之后隨著水體內(nèi)TDG過飽和度的降低,耗散速率越來越慢,直至趨近于飽和狀態(tài);轉(zhuǎn)速越大,紊動強(qiáng)度越大,耗散速度越快。轉(zhuǎn)速為600r/min時,水體在 26 min后 TDG降至 102%;轉(zhuǎn)速為400r/min時,水體在41min后TDG 降至102%;轉(zhuǎn)速為200 r/min時,需91minTDG方能降至102%。說明紊動是影響水體過飽和TDG釋放的重要因素之一。

3 過飽和TDG釋放的數(shù)學(xué)模型

關(guān)于過飽和TDG在下游河道的釋放過程,美國陸軍工程兵團(tuán)曾提出過飽和TDG在下游河道的釋放服從于一階動力學(xué)過程。公式表述為式中:G為計算時刻的TDG飽和度,%;G0為TDG初始飽和度,%;Geq為TDG平衡飽和度,%;t為滯留時間,h;kT為釋放系數(shù),h-1。

在Columbia河上,美國陸軍工程兵團(tuán)把kT的估算公式形式寫為

2 原型觀測

四川大學(xué)2007年對黃果樹瀑布下游白水河河道水體、二灘電站下游金沙江河道水體及銅街子電站下游大渡河河道水體開展了過飽和TDG釋放的原型觀測。

圖3為河道水體TDG的沿程釋放過程。圖3(a)表明,觀測期間黃果樹瀑布下游河道內(nèi)TDG飽和度較大,為114.4%。水流經(jīng)過約4 km的河道后,TDG飽和度即衰減至102.7%,平均每千米減小約2.2%,可見該瀑布下游河道過飽和TDG釋放速率相當(dāng)快。分析原因,一是由于河道水淺,水流中的過飽和溶解氣體與大氣充分交換所致;另外,河道內(nèi)分布著多個小落差跌水,跌水有助于水體與大氣的充分接觸和質(zhì)量交換,從而加速了過飽和TDG向大氣中釋放的速度。

圖3(b)表明,二灘電站壩下水流經(jīng)過93km河道后,TDG過飽和度僅降低7%,平均每千米減小約0.075%。圖3(c)表明,銅街子電站下游TDG過飽和度最大值為136.7%,水體經(jīng)過9.7 km河道后TDG飽和度降低至129%,平均每千米減小約0.7%?？梢?相對于金沙江的過飽和氣體沿程釋放速度,銅街子電站下游大渡河干流的過飽和氣體釋放速度較快。分析原因,銅街子電站下游河道的流速比二灘電站下游金沙江河道流速大而水深較小,因此加速了水體中TDG的耗散。由此說明在天然河道中,影響大壩泄洪時產(chǎn)生過飽和TDG耗散的主要因子是水深、紊動強(qiáng)度以及靜水壓強(qiáng)。由于試驗裝置內(nèi)水體深度與河道的深度有一定差別,導(dǎo)致了容器內(nèi)的凈水壓強(qiáng)比河道的靜水壓強(qiáng)小。式中:U為斷面平均流速,m/s;H為斷面平均水深,m;Dm為TDG的分子擴(kuò)散系數(shù),m2/s。

圖3 河道水體TDG釋放過程曲線

kT估算公式(式(2))的局限在于其僅考慮分子擴(kuò)散的作用,忽略了紊動擴(kuò)散的作用,因此應(yīng)用于大江大河時誤差較大。

采用二灘電站下游金沙江干流河段和銅街子電站下游大渡河干流河段原型觀測資料,對上述數(shù)學(xué)模型進(jìn)行參數(shù)率定,得到金沙江干流河段和大渡河干流河段的kT值分別為0.0678h-1和0.2013h-1,相關(guān)系數(shù)R2分別為0.9995和0.9694。率定結(jié)果如圖4所示。

通過對天然河道原型觀測值進(jìn)行擬合,得出金沙江干流河段和大渡河干流河段飽和TDG沿程kT分別為0.0678h-1和0.2013h-1。金沙江干流河段深度比大渡河干流河段大,而流速較小,因此金沙江干流過飽和TDG釋放系數(shù)比大渡河小。

根據(jù)馮鏡潔等[10]對金沙江、瀾滄江、岷江、長江等多條河流過飽和TDG釋放過程的原型觀測成果,考慮分子擴(kuò)散作用以外的紊動作用等在過飽和TDG釋放過程中的作用,對 kT的估算公式進(jìn)行修正,得到

式中:φT為考慮分子擴(kuò)散、紊動擴(kuò)散作用等的綜合系數(shù),量綱為T-1。

采用式(3)計算得到金沙江、大渡河原型觀測河段的 φT分別為4.43×10-10s-1和3.13×10-9s-1。

4 某電站下游過飽和TDG沿程釋放預(yù)測

金沙江中游某待建電站最大壩高138m,泄洪建筑物由左岸溢洪道和廠房左右兩側(cè)的泄洪沖沙底孔組成。根據(jù)相關(guān)研究成果,該電站通過表孔挑流單泄流量2188m3/s時,壩下生成的最大TDG飽和度為122.5%[11]。不考慮下游河道支流匯入情況,采用前述過飽和TDG釋放模型對下游過飽和TDG水流沿程釋放過程進(jìn)行預(yù)測。

預(yù)測模型采用式(1),其中kT根據(jù)式(3)計算得到。φT的取值參照前述金沙江河道原型觀測率定的成果4.43×10-10s-1。沿程斷面平均流速和平均水深采用縱向一維計算軟件HEC-RAS計算得到,由此得到過飽和TDG釋放系數(shù)的沿程變化情況,如圖5所示,預(yù)測得到壩下游TDG飽和度的沿程衰減過程如圖6所示。計算結(jié)果表明:經(jīng)過100km的天然河道后,水體TDG飽和度從122.5%降至115%;經(jīng)過172km后TDG飽和度降至110%?？梢?該工程泄洪產(chǎn)生的TDG過飽和的水體在下游河道的釋放速度較慢,影響范圍較大,需要進(jìn)一步采取措施減緩水體中過飽和TDG對魚類的影響。

圖5 下游河道kT沿程變化

圖6 下游河道TDG飽和度的沿程變化

5 結(jié) 語

電站泄洪造成的下游水體中TDG過飽和問題是水電開發(fā)的環(huán)境影響所面臨的一個重要問題。過飽和TDG沿程釋放規(guī)律研究是客觀公正地預(yù)測和評價水體中過飽和TDG對水生生態(tài)影響的重要基礎(chǔ)。本文在室內(nèi)機(jī)理試驗研究和對二灘電站下游河道、銅街子電站下游河道原型觀測的基礎(chǔ)上,對過飽和TDG沿程釋放的定量預(yù)測進(jìn)行探討,研究發(fā)現(xiàn)水深、靜水壓強(qiáng)以及紊動強(qiáng)度是影響氣體釋放的主要因素。過飽和TDG釋放速率隨水體深度的增加而減小,隨紊動強(qiáng)度的增加而增大。并且過飽和TDG釋放過程隨時間的推移逐漸減慢。

本文預(yù)測研究對定量評價水體中TDG的過飽和程度具有重要意義。由于過飽和TDG釋放研究起步較晚,室內(nèi)試驗和原型觀測資料還需要在今后進(jìn)一步完善和豐富。

本文采用的過飽和TDG釋放預(yù)測模型沒有考慮泄水的非恒定過程。同時過飽和TDG釋放除了受流速的影響外,還要受溫度、含沙量

、

河道地形等條件的影響。為了進(jìn)一步完善縱向一維過飽和TDG釋放預(yù)測模型,需積累更多的過飽和TDG原型觀測資料,并對kT及 φT的取值范圍作進(jìn)一步率定。同時,可以對泄水的非恒定流特征進(jìn)行研究,并開展二維或三維天然河道釋放過程預(yù)測,進(jìn)一步完善過飽和TDG沿程釋放模型。

[1]US Environment Protection Agency.Quality criteria for water[R].Washington,D.C.:Environment Protection Agency,1986.

[2]譚德彩.三峽工程致氣體過飽和對魚類致死效應(yīng)的研究[D].重慶:西南大學(xué),2006.

[3]US Army Corps of Engineers.Technical analysis of TDG processes[R].Jacksonville District:US Army Corps of Engineers-Northwest Division,Environmental Resources and Fish Planning Offices,2005.

[4]POLITANO M,CARRICA P,TUR AN C,et al.Prediction of the total dissolved gas downstream of spillways using a twophase flow model[C].World Water Congress.[s.l.].2004.

[5]WEBER L,HUANG H,LAI Y,et al.Modeling total dissolved gas production and transport downstream of spillway:threedimensional development and application[C].IAHR&INBO[s.l.],2004.

[6]LI Ran,LI Jia,LI Ke-feng,et al.Field observation of total dissolved gas supersaturation of high dams[C].33IAHR Congress.[s.l.],2008.

[7]蔣亮,李嘉,李然,等.紫坪鋪下游過飽和溶解氣體原型觀測研究[J].水科學(xué)進(jìn)展,2008,19(3):367-371.

[8]蔣亮,李然,李嘉,等.高壩下游水體中溶解氣體過飽和問題研究[J].四川大學(xué)學(xué)報:工程科學(xué)版,2008,40(5):69-73.

[9]程香菊,陳永燦.大壩泄洪下游水體溶解氣體超飽和理論分析及應(yīng)用[J].水科學(xué)進(jìn)展,2007,18(3):346-350.

[10]馮鏡潔,李然,李克鋒,等.高壩下游過飽和TDG釋放過程研究[J].水力發(fā)電學(xué)報,2010,29(1):7-12.

[11]李然,李嘉,李克鋒,等,高壩工程總?cè)芙鈿怏w過飽和預(yù)測研究[J].中國科學(xué),2009,39(12):1887-2006.