陳 鉑,彭夢(mèng)輝,王 華

(1.長沙理工大學(xué)水利工程學(xué)院,湖南 長沙 410004; 2.水沙科學(xué)與水災(zāi)害防治湖南省重點(diǎn)實(shí)驗(yàn)室,湖南 長沙 410004)

經(jīng)濟(jì)高速發(fā)展使得我國環(huán)境污染形勢(shì)嚴(yán)峻,而水污染事件已成為我國主要的突發(fā)性環(huán)境污染事件[1-8]。我國可建設(shè)的大型水電站共有1 360多座,其中很多都是高壩深庫。大型深水庫在蓄水后往往出現(xiàn)水體溫度分層,使水庫有著與河流不同的水溫分層結(jié)構(gòu)[9-11]。水庫建成運(yùn)行后,為處理類似新安江水庫苯酚泄漏、墨西哥灣漏油事件等突發(fā)污染事件造成的環(huán)境問題,需要準(zhǔn)確掌握溫度分層水體中污染物運(yùn)動(dòng)規(guī)律,對(duì)污染情況進(jìn)行預(yù)測(cè)和預(yù)警,為控制污染物擴(kuò)散提供事故處理決策支持。

很多突發(fā)環(huán)境事件中,污染物泄漏是以瞬時(shí)點(diǎn)源羽流(thermal)的形式擴(kuò)散。Scorer[12]研究了靜態(tài)均質(zhì)水體中的點(diǎn)源羽流云團(tuán)運(yùn)動(dòng)規(guī)律,認(rèn)為云團(tuán)的運(yùn)動(dòng)會(huì)在一個(gè)頂角約為30°的圓錐體范圍內(nèi);馬福喜等[13]通過試驗(yàn)觀測(cè)與大渦模擬也得到了相近的結(jié)果。Tarshish等[14]采用直接模擬方法研究了紊動(dòng)熱氣團(tuán)(turbulent thermals)運(yùn)動(dòng)規(guī)律,發(fā)現(xiàn)氣團(tuán)上升時(shí)均為橢球體,而氣團(tuán)加速上升時(shí)所受的浮力與橢球的高寬比相關(guān);林偉波[15]對(duì)密度分層流中的浮射流進(jìn)行了試驗(yàn),并對(duì)垂向密度呈線性變化的垂直浮力射流進(jìn)行了數(shù)值模擬,得到射流的密度場(chǎng)、射流軌跡、稀釋度和射流寬度等的分布規(guī)律;黃海寧等[16]為研究堿渣傾倒入海洋后引起的污染情況,在靜態(tài)溫度分層水體中進(jìn)行了堿渣云團(tuán)和泥漿云團(tuán)沉降試驗(yàn),認(rèn)為密度躍層對(duì)云團(tuán)有阻礙作用，當(dāng)分層強(qiáng)度大時(shí),云團(tuán)會(huì)停滯在密度躍層處振蕩;分層強(qiáng)度弱時(shí),云團(tuán)會(huì)穿過躍層,但部分云團(tuán)內(nèi)的細(xì)顆粒仍會(huì)滯留在躍層。

點(diǎn)源羽流云團(tuán)在溫度分層水體中的沉降主要受浮力影響,水溫結(jié)構(gòu)不同時(shí),浮力變化有很大差異,對(duì)云團(tuán)垂向運(yùn)動(dòng)規(guī)律影響很大。本文通過物理試驗(yàn),對(duì)靜態(tài)溫度分層水體中點(diǎn)源羽流云團(tuán)垂向運(yùn)動(dòng)規(guī)律進(jìn)行研究,為水庫和湖泊中污染物擴(kuò)散預(yù)測(cè)提供依據(jù)。

1 試驗(yàn)裝置與試驗(yàn)方法

試驗(yàn)用水箱如圖1所示,水箱中間位置布置有擋板,擋板下方先注入冷水,擋板閉合后在其上方注入熱水,待水體平靜后(約5 min)緩慢打開,形成溫度分層水體。使用電熱耦合溫度計(jì)測(cè)量水箱中間垂線上溫度隨時(shí)間變化情況,結(jié)果如圖2所示,可以看出不同溫差的分層水體,在50 min內(nèi)仍舊保持穩(wěn)定的分層結(jié)構(gòu)。污染物垂向擴(kuò)散試驗(yàn)通常不到1 min,故在試驗(yàn)中忽略熱量損失。所有試驗(yàn)中水箱中部溫度梯度最大可以達(dá)到0.5℃/mm。

圖1 試驗(yàn)水箱示意圖(單位:mm)

圖2 分層水體溫度場(chǎng)

試驗(yàn)使用氯化鈉和少量染色劑作為模擬污染物的溶質(zhì),與水混合形成模擬污染物溶液,模擬污染物室溫下初始密度為1.028～1.126 g/cm3。試驗(yàn)中假定污染物云團(tuán)的平均密度為云團(tuán)內(nèi)環(huán)境水體與污染物溶質(zhì)的質(zhì)量和云團(tuán)體積的比值。使用3種方式在水體表面投放模擬污染物:①試驗(yàn)組R1使用一無底圓筒(圖3(a)),內(nèi)部填充紗網(wǎng),試驗(yàn)時(shí)將圓筒大部置于水面以下,倒入5 mL污染物后提起,形成污染物云團(tuán);②試驗(yàn)組R2使用一硅膠軟管(圖3(b)),一端儲(chǔ)存污染物(0.4 mL)，另一端密封,試驗(yàn)時(shí)將有污染物一端置于水面以上5 cm,打開另一端釋放污染物,形成點(diǎn)源羽流云團(tuán);③試驗(yàn)組R3釋放方法與R2相同(圖3(c)),但其釋放口略低于水面,釋放污染物0.25 mL。由于投放污染物體積遠(yuǎn)小于環(huán)境水體,故認(rèn)為污染物投放后立即與上層水體混摻,形成的云團(tuán)內(nèi)部溫度與上層水體相同。

圖3 污染物表層投放方式

采用攝像機(jī)拍攝羽流云團(tuán)垂向運(yùn)動(dòng)過程,對(duì)試驗(yàn)圖像進(jìn)行分析并得到所需的各項(xiàng)參數(shù)。羽流云團(tuán)沉降時(shí)其形狀近似于橢球(圖4),與Tarshish等[14]數(shù)值模擬結(jié)果相同,可認(rèn)為云團(tuán)是水平向兩軸長度均為w、垂向高度為h的橢球。根據(jù)Beer-Lambert定律,將云團(tuán)圖像與原有背景之間的亮度差與污染物質(zhì)量相關(guān)聯(lián),定性分析污染物濃度分布[17]。

圖4 羽流云團(tuán)垂向運(yùn)動(dòng)示意圖

將試驗(yàn)照片(圖5)與原背景圖片逐像素比較,得到由染色劑引起的亮度變化,使用云團(tuán)內(nèi)最大亮度變化值的50%作為閾值消去噪點(diǎn)(孤立的亮度變化像素點(diǎn)),再將最大亮度變化值的20%作為閾值勾勒出云團(tuán)的邊界。每一幀圖片處理完后,將各像素點(diǎn)位置濃度分別在垂向和水平向疊加(圖6),并按照時(shí)間順序排列得到歷時(shí)變化過程,如圖7所示。根據(jù)邊界位置計(jì)算云團(tuán)寬度w和高度h,由濃度分布確定云團(tuán)重心位置,并根據(jù)重心位置變化計(jì)算云團(tuán)垂向沉速u。

圖5 原始試驗(yàn)照片

圖6 單幀圖片垂向濃度分布和水平濃度分布疊加示意圖

圖7 污染物云團(tuán)濃度分布疊加歷時(shí)變化

2 試驗(yàn)結(jié)果與分析

圖8 羽流云團(tuán)試驗(yàn)照片

圖9 試驗(yàn)環(huán)境水體溫度場(chǎng)

圖8為不同分層程度典型試驗(yàn)照片,不同水溫結(jié)構(gòu)的分層水體中,羽流云團(tuán)的運(yùn)動(dòng)規(guī)律有著顯著差異。當(dāng)上下層水溫差異較小時(shí),溫度梯度引起的浮力作用較弱,擋板開合產(chǎn)生的擾動(dòng)會(huì)在垂向傳播,最終形成近乎線性分層的水溫結(jié)構(gòu)(圖9中弱分層)。此時(shí),羽流云團(tuán)的運(yùn)動(dòng)(圖8(b)(f)(j))與均勻水體中的羽流云團(tuán)(圖8(a)(e)(i))相近,但垂向速度要小于不分層情況;同時(shí)由于線性分層水體內(nèi)的斜壓作用,橫向擴(kuò)散也受到限制。如果增大上下層水體溫差,此時(shí)水箱內(nèi)可以分為3層:上下層不同溫度的均勻水體和中間溫度線性分布的溫躍層,如圖9中的中度分層。此時(shí)的羽流云團(tuán)在溫躍層會(huì)受到明顯的浮力頂托。由于云團(tuán)內(nèi)部濃度分布不均勻,濃度小的部分(密度小于下層水體)停滯在溫躍層后會(huì)沿橫向繼續(xù)擴(kuò)散,濃度大的部分(密度大于下層水體)會(huì)因Rayleigh-Taylor不穩(wěn)定,產(chǎn)生尺寸較大的指狀液滴,侵入下層水體,如圖8(c)(g)(k)所示。上下層水體溫差很大時(shí),溫躍層厚度會(huì)變得很小,此時(shí)可認(rèn)為溫度場(chǎng)為階梯形分布,如圖9中的強(qiáng)分層。由于較強(qiáng)的浮力作用,羽流云團(tuán)會(huì)被阻滯在溫躍層,由于溫度梯度(上熱下冷)和濃度梯度(上層水含溶質(zhì),下層為清水)的共同作用產(chǎn)生雙擴(kuò)散對(duì)流(double-diffusion convection)[18],云團(tuán)下部會(huì)出現(xiàn)大量細(xì)小的指狀對(duì)流液滴,如圖8(d)(h)(l)所示。Chen等[19]測(cè)量了上述指狀侵入液滴的尺寸與侵入量(中度分層試驗(yàn)與強(qiáng)分層試驗(yàn)),認(rèn)為指狀侵入液滴的量約為總污染物量的5%～20%,且其寬度大于雙擴(kuò)散鹽指(salt-finger)寬度。顯然其測(cè)量的指狀液滴既有因雙擴(kuò)散對(duì)流產(chǎn)生的液滴,也包含了中度分層時(shí)Rayleigh-Taylor不穩(wěn)定產(chǎn)生的液滴。

羽流云團(tuán)垂向運(yùn)動(dòng)主要受浮力影響,其值為

(1)

u=f(B,N,zf)

(2)

使用量綱分析可以得到分層水體中羽流云團(tuán)的垂向沉速為

(3)

式中:c1、c2、c3均為系數(shù)。對(duì)于均勻水體,羽流云團(tuán)的垂向沉速為

(4)

式中:c4、c5為系數(shù)。根據(jù)試驗(yàn)采集數(shù)據(jù),采用文獻(xiàn)[20]中的方法計(jì)算每組試驗(yàn)中的系數(shù),所得結(jié)果如圖10和圖11(圖中白點(diǎn)為中位數(shù),黑色粗條為上、下四分位點(diǎn),外部為對(duì)稱的概率密度曲線)所示,系數(shù)c1、c2、c3、c4、c5中位數(shù)分別為0.44、-0.2、8.18、1.65和-3.46。

試驗(yàn)使用了3種不同的污染物投放方式,且污染物投放量也不相同,同時(shí)羽流云團(tuán)在不同環(huán)境水體中的運(yùn)動(dòng)規(guī)律有很大差異(圖8),為了解不同情況時(shí)羽流云團(tuán)的垂向運(yùn)動(dòng)規(guī)律,需要分別分析沉速公式中的系數(shù)。由圖10(a)和圖11(a)可以看出,投放方式R2試驗(yàn)得到的系數(shù)c1和c4與其他投放方式有明顯差異,顯示R2試驗(yàn)中采取水面以上5 cm投放污染物,污染液滴下落產(chǎn)生的較大初始動(dòng)量對(duì)云團(tuán)垂向運(yùn)動(dòng)有影響。而分層水體羽流云團(tuán)沉速公式(式(3))中的冪指數(shù)c2在不同投放方式和水溫結(jié)構(gòu)中變化不大,且都集中在-0.2附近;相反c3和c5的值雖然受投放方式或水溫分層影響不大,但數(shù)值較為離散。式(3)中的c1和c3以及式(4)中的c4和c5,與具體試驗(yàn)時(shí)的初始動(dòng)量、虛源點(diǎn)O位置(見圖4,綜合反映投放位置和投放量)和試驗(yàn)誤差相關(guān),不同情況會(huì)有一定差異;式(3)中的冪指數(shù)c2反映了各物理量之間的關(guān)系,受初始狀態(tài)影響較小。

圖12 典型試驗(yàn)云團(tuán)垂向沉速關(guān)系

圖10 分層水體沉速公式系數(shù)小提琴圖

圖11 均勻水體沉速公式系數(shù)小提琴圖

將分析得到的c1、c2和c3的中位數(shù)代入式(3),即可以預(yù)測(cè)分層水體中污染物云團(tuán)垂向運(yùn)動(dòng)。圖12為弱分層、投放方式為R1的云團(tuán)垂向沉速預(yù)測(cè)值和實(shí)際值,可見預(yù)測(cè)沉速與實(shí)際沉速趨勢(shì)相同,但由于不同條件下,污染物初始速度、污染物總量和污染物與環(huán)境水體密度差都不相同,導(dǎo)致了預(yù)測(cè)沉速與實(shí)際沉速之間的差異,這也是今后需要繼續(xù)研究的方向。

3 結(jié) 論

a. 不同水溫結(jié)構(gòu)中羽流云團(tuán)的垂向運(yùn)動(dòng)有很大差異,在線性分層水體中云團(tuán)下沉過程中橫向擴(kuò)散受到一定限制。

b. 在有溫躍層的水體中,云團(tuán)會(huì)在溫躍層位置停滯,由于云團(tuán)與環(huán)境水體密度之間的差異,可能出現(xiàn)Rayleigh-Taylor不穩(wěn)定或雙擴(kuò)散,并出現(xiàn)指狀侵入,且其尺寸較鹽指要大。