朱 杰,權(quán) 錦,龍 巖,雷曉輝,孔令仲
(1.北京工業(yè)大學(xué),北京 100022;2.中國水利水電科學(xué)研究院,北京 100038)
突發(fā)水污染事件一般發(fā)生突然,污染危害較大,難以控制[1]。研究突發(fā)水污染事件的關(guān)鍵,在于研究事故發(fā)生后,水流的特性以及污染物在水體中的輸移擴散過程[2, 3]?,F(xiàn)階段研究水流特性以及污染物輸移擴散規(guī)律的主要方法有:理論分析、模型試驗和數(shù)值模擬等[4]。理論分析方法主要是建立水流運動方程和污染物輸移擴散方程,從而得到方程解析解[5, 6];模型試驗是一種重要研究手段[7, 8],但是模型比尺效應(yīng)[9]問題至今未能得到很好解決;數(shù)值模擬是一種較高效的研究方法,比物理模型、試驗建模速度快、周期短,計算效率通常取決于計算機的CPU效率,數(shù)值模擬便于修改模擬工況,能展現(xiàn)多種工況下的情景,因此,用數(shù)值模擬來研究污染物在水體中的輸移擴散過程越來越被廣泛應(yīng)用[10, 11]。
高學(xué)平、張晨[11]等運用數(shù)值模型對引黃濟津調(diào)水工程河道水質(zhì)進行了預(yù)測,并根據(jù)污染物峰值濃度的變化情況,分析了污染物輸移擴散及降解轉(zhuǎn)化規(guī)律。房彥梅等[12]利用一維水力學(xué)和水質(zhì)數(shù)值模型方法對南水北調(diào)中線工程典型渠段突發(fā)性水污染事故應(yīng)急控制策略進行了數(shù)值模擬,分析不同控制策略下渠道退水量、穩(wěn)定時間、水位變幅、水質(zhì)等參數(shù)的變化情況,制定保證明渠安全運行的水污染事故應(yīng)急閘控方案,為中線工程應(yīng)急預(yù)案的建立提供參考。宋國棟[13]以南水北調(diào)中線污染物運移特性研究為背景,在現(xiàn)場調(diào)查及模型試驗?zāi)M的基礎(chǔ)上,研究了突發(fā)性水污染事件發(fā)生后污染物的輸移規(guī)律,并給出了突發(fā)性水污染事件發(fā)生后的處置措施。王興偉等[14]應(yīng)用Mike水質(zhì)模型軟件對污染物在干渠內(nèi)的污染擴散情況進行模擬,對突發(fā)事故的影響范圍、程度、時間做出定量預(yù)報,揭示突發(fā)事故發(fā)生后,節(jié)制閘未參與調(diào)控情況下的污染物輸移擴散規(guī)律。張秋燕[15]研究了河流中污染物隨時間和空間輸移擴散的規(guī)律及其影響因素,并建立了一維恒定流水質(zhì)模型,通過隱式有限差分法求解該模型。此外,張秋燕以BOD5為預(yù)測因子,將數(shù)值模擬結(jié)果與QUAL2K綜合水質(zhì)模型軟件的模擬結(jié)果進行比較,數(shù)值模擬與軟件模擬結(jié)果基本吻合,表明所建立的水質(zhì)模型能夠較準確的模擬河流中污染物輸移擴散規(guī)律。
突發(fā)水污染事件[16]具有突發(fā)性、多樣性及不確定性,突發(fā)水污染事件發(fā)生后,需要快速采取有效的響應(yīng)措施[17-19],因而快速預(yù)測[20]出污染團的輸移擴散特征指標,對科學(xué)、高效的制定出響應(yīng)策略非常重要。然而目前污染物輸移擴散研究多數(shù)集中在天然河流或者人工調(diào)水渠道的順坡段,對平坡和逆坡的渠道中污染物輸移擴散研究相對較少。
本文運用HECRAS軟件進行數(shù)值模擬,分析平坡和逆坡輸水渠道的污染物輸移擴散規(guī)律;并依據(jù)現(xiàn)有的長距離輸水渠道污染物輸移擴散快速計算公式(龍巖[21, 22],主要面向順坡渠道突發(fā)水污染計算),分析該快速計算公式在平坡和逆坡輸水渠道中的適用性,研究結(jié)果為長距離輸水工程突發(fā)水污染事件應(yīng)急調(diào)控提供快速決策支持。
南水北調(diào)中線全長1 432 km,總干渠渠道主要為全斷面混凝土襯砌的梯形斷面,渠道底寬由28.5 m遞減為7 m,渠道的邊坡系數(shù)也是在不斷變化的,渠道的邊坡系數(shù)變化范圍是0~3.5,其中邊坡系數(shù)為2占44%;中線全線采用順坡渠道,底坡變化范圍大概是1/30 000~1/10 000,渠道底坡有80%以上的位于區(qū)間1/30 000~1/2 000。南水北調(diào)東線輸水干線總長1 156 km,建設(shè)13個梯級泵站,全線最高處東平湖水位與長江水位差約40 m,渠道底坡為逆坡和平坡。本文以南水北調(diào)中線干線渠道斷面特征為基礎(chǔ),分別采用底坡為平坡和逆坡的渠道作為研究對象,渠道的基本要素和初始條件如表1所示。
表1 渠道的基本要素表Tab.1 Basic elements of the channel
艾恒雨[23]對2001-2011年國內(nèi)發(fā)生的1 176起重大突發(fā)性水污染事件進行統(tǒng)計分析,得出突發(fā)性水污染事件的主要污染物是工業(yè)廢水、綜合廢水和油類,分別占比31%、14%和13%;并且污染水團的排放量一般能達到數(shù)百公斤到數(shù)十噸之間。
基于此研究,本文假定污染情景為可溶性污染物、量級為10 t。在數(shù)值模擬中,渠道的離散系數(shù)取值根據(jù)經(jīng)驗公式[21, 24, 25]計算所得。
渠道內(nèi)突發(fā)可溶性水污染事件的污染水團模擬特征如表2所示。
表2 渠道內(nèi)突發(fā)可溶性水污染事件的污染水團模擬特征Tab.2 Simulation characteristics of polluted water masses in sudden water pollution incidents in channel
采用美國陸軍工程兵團開發(fā)的河道水力分析軟件HEC-RAS,構(gòu)建一維水動力水質(zhì)模型。在軟件中,水動力運動方程和污染物輸移擴散方程為:
(3)
式中:c為污染水團濃度,mg/L;c0為側(cè)向入流濃度,mg/L;g為重力加速,m/g2;q為單位長度渠道上的側(cè)向入流流量,m2/s;t為時間,s;x為空間距離,m;z為水位,m;A為過流面積,m2;C為謝才系數(shù),m1/2/s;DL為縱向離散系數(shù),m2;Q為渠道流量,m3/s;R為水力半徑,m;S為源匯項。
基于正常輸水情況下的污染物輸移擴散規(guī)律,通常用污染物峰值輸移距離、縱向長度和峰值濃度來描述污染物的輸移擴散過程。圖1中給出了污染物峰值輸移距離、縱向長度以及峰值濃度的涵義。
圖1 污染物輸移擴散特征Fig.1 Pollutant transport characteristics
2.2.1 峰值輸移距離快速計算公式
研究發(fā)現(xiàn)[22]峰值輸移距離與渠道長度、底坡、流量以及傳播時間成正比,與渠道邊坡系數(shù)、底寬和水深成反比,但渠道的長度、底坡對峰值輸移距離影響甚微,而流量與渠道流速成正比,邊坡系數(shù)、底寬和水深與渠道流速成反比,因此認為峰值輸移距離與渠道流速、水流的傳播時間成正比,其關(guān)系式[22]為:
D=vT
(4)
式中:D為峰值輸移距離,m;v為渠道內(nèi)水流速度,m/s;T為傳播時間,s。
2.2.2 污染物縱向長度快速計算公式
根據(jù)文獻[26]可知,一維瞬時污染物擴散中,污染物濃度隨流動時間呈正態(tài)分布。由正態(tài)分布曲線間的面積比例可知:在4σ(σ為正態(tài)分布的標準差)范圍內(nèi)包含污染物總量的95.44%,在6σ范圍內(nèi)包含污染物總量的99.74%,定義kσ為污染物離散寬度(k為離散寬度系數(shù))。
T時間內(nèi)污染水團沿渠道縱向長度W可表達[27]為:
(5)
式中:DL為縱向離散系數(shù),m2/s;t為時間,s;a為k/2。
2.2.3 污染物峰值濃度快速計算公式
峰值濃度[27]為:
(6)
假設(shè)初始濃度為C0,某時刻峰值濃度為Cm,M為污染物投放量,則定義P為峰值濃度與初始濃度百分比,即:
(7)
經(jīng)該提出者驗證,該公式對突發(fā)水污染事件下長距離順坡輸水渠道的污染物輸移擴散模擬效果良好。但并未研究逆坡、平坡渠道中的適用性。
為研究平坡渠道中的污染物輸移擴散,本文設(shè)置了兩種平坡渠道方案(渠段長度為25 km和50 km。)和污染情景,平坡渠道參數(shù)如表1所示。
經(jīng)模擬計算及運用快速計算公式計算,污染物峰值輸移距離D模擬值及計算值對比結(jié)果如圖2所示??梢钥闯觯斍赖灼聻槠狡聲r,峰值輸移距離D與污染物擴散時間成正比,與渠道長度成反比; HECRAS軟件模擬值與順坡快速計算公式計算值變化過程基本一致,計算結(jié)果較為接近,偏差值較小。
計算值與模擬值的偏差如圖3所示??梢钥闯?,渠段長度為25 km時,偏差值小于5%;渠段長度為50 km時,模擬值與計算值的偏差小于20%。因此,在渠道底坡為平坡時,污染物峰值輸移距離可以用順坡渠道峰值輸移距離快速計算公式進行預(yù)測,當渠段長度規(guī)模為25 km以下時,預(yù)測效果較好,渠段越長,偏差值越大。
圖2 平坡渠道計算值與模擬值對比Fig.2 Comparison of calculated values of flat slope channels with simulated values
圖3 平坡渠道計算值與模擬值偏差Fig.3 Deviation between calculated value and simulated value of flat slope channel
經(jīng)模擬計算及運用快速計算公式計算,污染物縱向長度W隨傳播時間的變化過程的模擬結(jié)果及計算結(jié)果如圖4所示??梢钥闯?,當渠道底坡為平坡時,污染物縱向長度隨時間成正相關(guān)變化,且隨時間的增加,渠道越長,污染物縱向長度就越大,但不同長度的渠道間差別較小。模擬值與計算值的變化規(guī)律基本相同,數(shù)值較為接近。
對模擬值與計算值的偏差進行比較分析(圖5),可以發(fā)現(xiàn),快速計算公式計算值與模擬值在渠段長度為25 km時的相對偏差在5%左右,在渠段長度為50 km時,偏差值低于10%。因此,順坡渠道污染物縱向長度快速計算公式可用于平坡渠段可溶性污染團的污染物縱向長度指標快速預(yù)測,且效果較好。
圖4 平坡渠道計算值與模擬值對比Fig.4 Comparison of calculated values of flat slope channels with simulated values
圖5 平坡渠道計算值與模擬值偏差Fig.5 Deviation between calculated value and simulated value of flat slope channel
經(jīng)模擬計算及運用快速計算公式計算,污染物峰值濃度C模擬值及計算值對比結(jié)果如圖6所示。
圖6 平坡渠道峰值濃度模擬值與計算值對比Fig.6 Comparison of simulated peak value and calculated value of flat slope channel
從圖中可以看出,相同離散系數(shù)下,污染團峰值濃度與時間成反比,同種底坡情況下,污染團峰值濃度與渠道長度成反比??焖俟接嬎阒蹬c模擬值高度吻合,基本一致。圖7反映了平坡情況下,某一時刻的濃度值與初始濃度的比值P的變化過程,可以看出,濃度比值P隨時間的增加逐漸減小,反映了污染物濃度逐漸降低,這與圖6的結(jié)果相吻合。同時圖7還可以看出,模擬值與快速公式計算得到的P值結(jié)果接近,偏差值較小。因此,在渠道底坡為平坡時,污染物峰值濃度及濃度比值P仍然可以用順坡渠道快速計算公式計算,且效果較好。
圖7 平坡渠道濃度比值P模擬值與計算值Fig.7 Flat slope channel simulation and concentration ratio of P value
對污染團峰值濃度的模擬值與計算值偏差進行比較分析(圖8),可以發(fā)現(xiàn),快速計算公式計算值與模擬值在渠段長度為25 km時的相對偏差低于10%,在渠段長度為50 km時,偏差值在12%左右,兩種方案的偏差都有逐步減小的趨勢。峰值濃度比值P模擬值與計算值的偏差對比如圖9。從圖9中可以看出,渠段長度為25 km時的偏差在10%以下,渠段長度為50 km時偏差在12%左右,且偏差逐漸減小。因此,平坡渠道污染物峰值濃度及濃度比值P的快速預(yù)測可以用順坡渠道快速計算公式,且效果較好。
圖8 平坡渠道峰值濃度模擬值與計算值偏差Fig.8 Simulation flat slope channel peak concentration and calculated deviation
圖9 平坡渠道濃度比值P模擬值與計算值偏差Fig.9 Simulation flat slope channel concentration ratio and the calculated value of the deviation of P
為研究逆坡渠道中的污染物輸移擴散,本文設(shè)置了兩種逆坡渠道方案(渠段底坡為-0.000 01和-0.000 04。)和污染情景,逆坡渠道參數(shù)如表1所示。
經(jīng)模擬計算及運用快速計算公式計算,污染物峰值輸移距離D模擬值及計算值對比結(jié)果如圖10所示。從圖結(jié)果所可以看出,在逆坡渠道的不同方案下,污染團峰值輸移距離D與污染物擴散時間呈線性關(guān)系,并與渠道底坡坡度成反比。
HECRAS軟件模擬值與快速計算公式計算值變化情況一致,變化過程較接近,偏差較小。計算值與模擬值的偏差如圖11所示,渠道底坡為-0.000 04時(方案4),偏差在20%左右,且隨時間的推移,偏差越來越?。划?shù)灼聻?0.000 01時(方案3),偏差在12%左右,且逐漸減小。因此,順坡渠道峰值輸移距離快速計算公式可用于快速預(yù)測逆坡渠道污染物峰值輸移距離,且效果較好。
圖10 逆坡渠道計算值與模擬值對比Fig.10 Comparison of calculated value and simulated value of reverse slope channel
圖11 逆坡渠道計算值與模擬值偏差Fig.11 Deviation between calculated value and simulated value of reverse slope channel
經(jīng)模擬計算及運用快速計算公式計算,污染物縱向長度W模擬值及計算值對比結(jié)果如圖12所示。污染物縱向長度隨時間成正相關(guān)變化,不同底坡的渠道間縱向長度差別較小。模擬值與計算值的變化規(guī)律基本一致,數(shù)值較為接近。
圖13反映了運用順坡渠道污染物快速計算公式計算值與模擬值的相對偏差值,可以發(fā)現(xiàn),渠段底坡為-0.000 01和-0.000 04時,運用快速計算公式得到的結(jié)果與模擬結(jié)果相對偏差均在10%左右,且隨時間增加,偏差值逐漸減小。因此,可以用順坡渠道污染物縱向長度快速計算公式快速預(yù)測逆坡渠道中污染物的縱向長度,且效果較好。
圖12 逆坡渠道計算值與模擬值對比Fig.12 Comparison of calculated value and simulated value of reverse slope channel
圖13 逆坡渠道計算值與模擬值偏差Fig.13 Deviation between calculated value and simulated value of reverse slope channel
經(jīng)模擬計算及運用快速計算公式計算,污染物峰值濃度C模擬值及計算值對比結(jié)果如圖14所示。
從圖中可以看出,相同離散系數(shù)下,污染物峰值濃度與時間成反比??焖儆嬎愎接嬎阒蹬c模擬值變化過程基本一致,峰值濃度隨時間的增加逐漸降低。圖15反映了逆坡情況下,某一時刻的濃度與初始濃度的比值P的變化過程,濃度比值P隨時間的增加逐漸減小,意味著污染物濃度逐漸降低,與圖14所反映出來的結(jié)果一致,圖15中,模擬值與快速公式計算得到的峰值濃度P值結(jié)果接近,偏差值較小。
對污染團峰值濃度的模擬值與計算值偏差進行比較分析(圖16),可以發(fā)現(xiàn),快速計算公式計算值與模擬值在渠道底坡為-0.000 04時的相對偏差低于5%,渠段底坡為-0.000 01時,偏差值在10%左右。峰值濃度比值P模擬值與計算值的偏差對比如圖17。從圖中可以看出,渠段底坡為-0.000 04時的偏差在5%以下,渠段底坡為-0.000 01時偏差在10%左右,且偏差逐漸減小。因此,逆坡渠道污染物峰值濃度及濃度比值P的快速預(yù)測可以運用順坡渠道快速計算公式,且效果較好。
圖14 逆坡渠道峰值濃度模擬值與計算值對比Fig.14 Comparison between simulated value and calculated value of peak concentration in reverse slope channel
圖15 逆坡渠道濃度比值P模擬值與計算值Fig.15 P simulation value and calculation value of the concentration ratio of reverse slope channel
圖16 逆坡渠道峰值濃度模擬值與計算值偏差Fig.16 Deviation between simulated peak value and calculated value of reverse slope channel
圖17 逆坡渠道濃度比值P模擬值與計算值偏差Fig.17 Inverse slope channel concentration ratio P simulation value and calculation deviation
總結(jié)國內(nèi)外對于污染團輸移擴散的研究成果,開展平、逆坡渠道中的污染物輸移擴散研究,基于已有的順坡渠道污染物輸移擴散快速計算公式,研究其在平坡、逆坡渠道情況下的適用性,并將計算結(jié)果與HECRAS軟件模擬結(jié)果進行對比,嘗試尋求平、逆坡渠道中污染物輸移擴散的快速計算方法;利用HECRAS軟件進行一維水動力水質(zhì)數(shù)值模型模擬,研究了平坡和逆坡渠道中污染物峰值輸移距離、縱向長度及峰值濃度變化規(guī)律,并與順坡渠道快速計算公式計算值進行了對比,結(jié)果表明:
(1)渠道為平坡和逆坡時,峰值輸移距離的模擬值與計算值有一定差別,當渠段規(guī)模為25 km以下時,偏差值在5%以下,當渠段規(guī)模為50 km時,偏差值在20%以下,因此,順坡渠道峰值輸移距離快速計算公式適用于平坡或逆坡渠道,并在渠段規(guī)模小于25 km時計算效果較好。
(2)在計算平坡或逆坡渠道的污染物縱向長度時,模擬值與順坡渠道計算公式相吻合,偏差值在10%左右,因此,可以用順坡渠道污染物縱向長度快速計算公式快速預(yù)測平坡和逆坡渠道的污染物輸移擴散指標,且效果較好。
(3)在相同的離散系數(shù)下,平坡或逆坡渠道中,污染物峰值濃度數(shù)值模擬值與運用順坡渠道快速計算公式計算值基本一致,因此,在快速預(yù)測平坡和逆坡渠道中的污染物峰值濃度及濃度比值P值時,可以用順坡渠道快速計算公式計算。
□