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

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

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

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

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

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

1.1 試驗(yàn)裝置和試驗(yàn)條件

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

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

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

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

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

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

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

2 原型觀測(cè)

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

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

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

圖3 河道水體TDG釋放過(guò)程曲線

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

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

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

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

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

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

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

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

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

圖5 下游河道kT沿程變化

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

5 結(jié) 語(yǔ)

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

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

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

、

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

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

[2]譚德彩.三峽工程致氣體過(guò)飽和對(duì)魚(yú)類(lèi)致死效應(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]蔣亮,李嘉,李然,等.紫坪鋪下游過(guò)飽和溶解氣體原型觀測(cè)研究[J].水科學(xué)進(jìn)展,2008,19(3):367-371.

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

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

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

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