明渠混合污染物側(cè)向和垂向擴(kuò)散系數(shù)的計算方法及其應(yīng)用

武周虎，胡德俊，徐美娥

（青島理工大學(xué)環(huán)境與市政工程學(xué)院，山東青島 266033）

明渠混合污染物側(cè)向和垂向擴(kuò)散系數(shù)的計算方法及其應(yīng)用

武周虎，胡德俊，徐美娥

（青島理工大學(xué)環(huán)境與市政工程學(xué)院，山東青島 266033）

明渠的側(cè)向和垂向擴(kuò)散系數(shù)是衡量岸邊水流對污染物質(zhì)混合輸移能力的重要水質(zhì)參數(shù)之一，其值的準(zhǔn)確與否直接關(guān)系到明渠水質(zhì)預(yù)測預(yù)報成果的可靠性。基于污染混合區(qū)的理論計算方法，推導(dǎo)了污染混合區(qū)外邊界標(biāo)準(zhǔn)曲線和曲面的統(tǒng)一方程，包含最大長度Ls、最大寬度bs和最大深度ds等特征尺度。該曲線形狀近似于半橢圓，曲面形狀為近似橢球體的一部分，表明污染混合區(qū)具有相似性。給出了由岸邊污染混合區(qū)外邊界最大長度、最大寬度或最大深度和平均流速確定側(cè)向或垂向擴(kuò)散系數(shù)的計算公式，提出了采用污染混合區(qū)面積或體積進(jìn)行總體控制的側(cè)向或垂向擴(kuò)散系數(shù)計算方法和采用水面橫向積分濃度確定垂向擴(kuò)散系數(shù)的實用方法。通過現(xiàn)場觀測結(jié)果分析，給出了府河下游河段枯水期的側(cè)向擴(kuò)散系數(shù)0．27 m2／s。

明渠；污染混合區(qū)；標(biāo)準(zhǔn)曲線；側(cè)向擴(kuò)散系數(shù)；垂向擴(kuò)散系數(shù)；計算方法

1 概 述

明渠沿岸地區(qū)工業(yè)或城鎮(zhèn)生活污水經(jīng)處理達(dá)到相應(yīng)的排放標(biāo)準(zhǔn)后，大多數(shù)情況是通過管道或明渠實施岸邊排放。污水首先在排污口近區(qū)稀釋混合，其次在水域的長度與寬度、深度方向逐漸移流擴(kuò)散，在排污口附近水域會形成污染混合區(qū)［1］。對于寬闊河流的污染混合區(qū)分布多屬于二維問題［2，3］，而對于寬闊深水水庫的污染混合區(qū)分布則多屬于三維問題［4，5］。另一方面，明渠岸邊附近水體常常是人們生產(chǎn)生活用水對水質(zhì)要求較高的區(qū)域，對大江大河來說，“總體水質(zhì)”不超標(biāo)，并不意味著“岸邊水質(zhì)”不超標(biāo)?！鞍哆吽|(zhì)”對應(yīng)“岸邊環(huán)境容量”［6］。

明渠的橫向和垂向擴(kuò)散系數(shù)是衡量水流對污染物質(zhì)混合輸移能力的重要水質(zhì)參數(shù)之一，其值的準(zhǔn)確與否直接關(guān)系到明渠水質(zhì)預(yù)測預(yù)報成果的可靠性。目前國內(nèi)外確定河流橫向和垂向擴(kuò)散系數(shù)的主要方法有理論公式、經(jīng)驗公式法和示蹤實驗法，艾爾德采用對數(shù)流速分布函數(shù)，在各向同性紊動條件下得出垂向擴(kuò)散系數(shù)為Ez＝0．067HU （式中H為斷面平均水深，U 為摩阻流速）。由于無量綱橫向擴(kuò)散系數(shù)經(jīng)驗公式αy＝Ey／HU 一般認(rèn)為其值為0．3～0．9［7］，所以橫向擴(kuò)散系數(shù)的計算結(jié)果是一個范圍值，其準(zhǔn)確性難以判定。示蹤實驗法分為現(xiàn)場實驗和室內(nèi)實驗，示蹤實驗完成后可采用矩法［8－14］、直線圖解法［15，16］、線性回歸法［17］、曲線擬合法［18］、遺傳算法［19］、人工神經(jīng)網(wǎng)絡(luò)法［20］等計算橫向擴(kuò)散系數(shù)?，F(xiàn)場實驗法計算結(jié)果比較可靠，但河流實驗受示蹤劑的投放、取樣等條件的限制；室內(nèi)實驗多用于水槽中橫向擴(kuò)散系數(shù)與粗糙形式、寬度、水深、流速、摩阻流速等水力要素的關(guān)系研究。鄭旭榮等［21］借助于拋物線型斷面形態(tài)方程，提出了橫向紊動擴(kuò)散系數(shù)的斷面分布及其平均值表達(dá)式?，F(xiàn)行關(guān)于橫向擴(kuò)散系數(shù)的關(guān)系式和數(shù)值都是針對河流全斷面混合特性的研究成果。實際工程中常見的多為岸邊排放，對于寬闊的明渠（如長江、黃河、三峽水庫等）水流對污染物質(zhì)混合輸移能力主要取決于岸邊混合擴(kuò)散特性，大江大河岸邊污染混合區(qū)和岸邊環(huán)境容量的計算往往只涉及到全斷面寬度的數(shù)十分之一甚至于數(shù)百分之一，選擇采用全斷面平均水深等水力要素確定的橫向擴(kuò)散系數(shù)顯然不盡合理，側(cè)向擴(kuò)散系數(shù)比橫向擴(kuò)散系數(shù)能更好地反映岸邊水流的混合擴(kuò)散特性。如三峽水庫萬州段的平均水深為71m，而岸邊擴(kuò)散區(qū)的水深比較淺，并且往往處于一個由傾斜岸坡形成的角形變化域。

本文以順直寬闊的明渠水流對污染物質(zhì)的混合輸移研究為背景，基于武周虎等［2－5］提出的岸邊污染混合區(qū)理論計算方法，推導(dǎo)污染混合區(qū)外邊界等濃度標(biāo)準(zhǔn)曲線和曲面的統(tǒng)一方程，提出確定側(cè)向擴(kuò)散系數(shù)和垂向擴(kuò)散系數(shù)的計算公式。采用府河岸邊污染混合區(qū)的現(xiàn)場觀測結(jié)果計算其側(cè)向擴(kuò)散系數(shù)，為進(jìn)一步解決大江大河、水庫岸邊污染混合區(qū)和岸邊環(huán)境容量的計算問題，提供一定的參考研究。

2 側(cè)向擴(kuò)散系數(shù)計算方法

2．1 河流二維問題

式中：x為自排污口沿河流主流向的縱坐標(biāo)；y為垂直于x的橫向坐標(biāo)，坐標(biāo)原點取在岸邊排污點；m為單位時間的排污強(qiáng)度；U為岸邊水域平均流速；H為平均水深；Ey為岸邊排污的橫向擴(kuò)散系數(shù)，又稱側(cè)向擴(kuò)散系數(shù)①。

由文獻(xiàn)［2］給出的岸邊污染混合區(qū)最大長度Ls、最大寬度bs和相應(yīng)縱坐標(biāo)Lc以及面積S的理論計算公式分別為：

式中Cd為河流排污引起的允許濃度升高值，該值等于水環(huán)境功能區(qū)所執(zhí)行的濃度標(biāo)準(zhǔn)值Cs減去背景濃度Cb，即Cd＝Cs－Cb，其等濃度線所包圍的區(qū)域為污染混合區(qū)。由式（2）可以得到岸邊污染混合區(qū)最大長度與最大寬度的比值為

其中Pe＝ULs／Ey稱為貝克來數(shù)，表征物質(zhì)的縱向移流通量與側(cè)向擴(kuò)散通量的比值［3］。即岸邊污染混合區(qū)的長寬比與貝克來數(shù)Pe的0．5次方成正比。

將式（2）代入式（1）并令C＝Cd，化簡整理得到河流岸邊污染混合區(qū)外邊界標(biāo)準(zhǔn)曲線（見圖1）方程為

圖1 岸邊污染混合區(qū)外邊界標(biāo)準(zhǔn)曲線Fig．1 The standard curve of pollutantm ixing zone’s outer boundary in the bank side

由圖1可以看出，河流岸邊污染混合區(qū)外邊界標(biāo)準(zhǔn)曲線形狀近似于半橢圓，在靠近排污口一端出現(xiàn)鈍頭，在污染混合區(qū)下游邊界出現(xiàn)稍尖形狀，污染混合區(qū)最大寬度相應(yīng)的縱坐標(biāo)為最大長度的1／e≈0．368，說明污染混合區(qū)具有相似性。

由式（2）可以得到順直寬闊河流側(cè)向擴(kuò)散系數(shù)的計算公式為：

式（5）給出了側(cè)向擴(kuò)散系數(shù)的3個計算公式。表明順直寬闊河流側(cè)向擴(kuò)散系數(shù)與岸邊污染混合區(qū)最大寬度的2次方成正比，與平均流速的1次方成正比，與最大長度的1次方成反比。說明在平均流速相同的條件下，污染混合區(qū)最大長度越大側(cè)向擴(kuò)散系數(shù)越小，最大寬度和面積越大側(cè)向擴(kuò)散系數(shù)越大。在實際應(yīng)用中，可根據(jù)岸邊污染混合區(qū)的現(xiàn)場觀測結(jié)果采用式（5）計算側(cè)向擴(kuò)散系數(shù)。

通常，受天然河道地形的影響，污染混合區(qū)的形狀一般不會像式（4）、圖1所示的那么規(guī)則，污染混合區(qū)最大寬度的測量受多種因素的影響也會存在變數(shù)。為了減小污染混合區(qū)最大寬度測量不精確給側(cè)向擴(kuò)散系數(shù)計算帶來的誤差，可以根據(jù)現(xiàn)場觀測數(shù)據(jù)繪制污染混合區(qū)外邊界等濃度曲線，在污染混合區(qū)形狀大體為上游出現(xiàn)鈍頭和下游出現(xiàn)稍尖的近似橢圓形情況下，即在最大寬度相應(yīng)的縱坐標(biāo)大約為最大長度的1／e時確定其面積（其它形狀污染混合區(qū)的側(cè)向擴(kuò)散系數(shù)較為復(fù)雜，另文研究），然后按照式（5）計算側(cè)向擴(kuò)散系數(shù)，即采用污染混合區(qū)面積進(jìn)行總體控制的側(cè)向擴(kuò)散系數(shù)計算方法，下同。

①M(fèi)ixing Zone Technical Procedures Manual（DRAFT），Idaho Department of Environmental Quality，USA，2008．8

2．2 傾斜岸水庫三維問題

式中：x為沿水庫主流向的縱坐標(biāo)；y，z為垂直于x的橫向和垂向坐標(biāo)，坐標(biāo)原點取在庫岸水面排污點；β為角域映射系數(shù)，當(dāng)岸坡線與水平線之間的夾角為θ時，β＝360／θ，其他符號含義同前。

由文獻(xiàn)［4］可知，傾斜岸水庫岸邊污染混合區(qū)在z＝0水面上最大長度Ls、最大寬度bs、最大深度ds和相應(yīng)縱坐標(biāo)Lc以及面積S的理論計算公式依次為：

傾斜岸水庫岸邊污染混合區(qū)外邊界等濃度標(biāo)準(zhǔn)曲線（寬度）方程與河流二維問題的式（4）、圖1相同，污染混合區(qū)的空間形狀為該等濃度標(biāo)準(zhǔn)曲線在θ角度上的旋轉(zhuǎn)體，是一近似橢球體的扇形截圖。

由式（7）可以得到順直寬闊傾斜岸水庫側(cè)向擴(kuò)散系數(shù)的計算公式為

順直寬闊傾斜岸水庫側(cè)向擴(kuò)散系數(shù)與平均流速和水面上污染混合區(qū)特征尺度的關(guān)系與河流情況相同，只是系數(shù)減半。在實際應(yīng)用中，可根據(jù)傾斜岸水庫岸邊水面上污染混合區(qū)的現(xiàn)場觀測結(jié)果采用式（8）計算側(cè)向擴(kuò)散系數(shù)。

2．3 垂直岸水庫三維問題

式中符號含義同前。

由文獻(xiàn)［5］可知，垂直岸水庫岸邊污染混合區(qū)在z＝0水面上最大長度Ls、最大寬度bs以及面積S和在y＝0垂直面上最大深度ds以及相應(yīng)縱坐標(biāo)Lc的理論計算公式依次為：

垂直岸水庫岸邊污染混合區(qū)外邊界等濃度標(biāo)準(zhǔn)曲線（寬度和深度）方程分別與河流二維問題的式（4）、圖1相同，污染混合區(qū)的空間形狀為1／4近似橢球體的扁蛋形。將式（10）代入式（9）并令C＝Cd，化簡整理得到該曲面方程為

由式（10）可以得到順直寬闊水庫側(cè)向擴(kuò)散系數(shù)的計算公式為：

順直寬闊垂直岸水庫側(cè)向擴(kuò)散系數(shù)與平均流速和水面上污染混合區(qū)特征尺度的關(guān)系與傾斜岸水庫情況相同。在實際應(yīng)用中，可根據(jù)垂直岸水庫岸邊水面上污染混合區(qū)的現(xiàn)場觀測結(jié)果采用式（12）計算側(cè)向擴(kuò)散系數(shù)。

3 垂向擴(kuò)散系數(shù)計算方法

針對垂直岸水庫三維問題，在橫向和垂向擴(kuò)散系數(shù)不相等的情況下，討論垂向擴(kuò)散系數(shù)的計算方法。

3．1 垂直面上污染混合區(qū)法

根據(jù)垂直面上岸邊污染混合區(qū)外邊界標(biāo)準(zhǔn)曲線的特征尺度，按照與側(cè)向擴(kuò)散系數(shù)相同的方法，由式（10）可以得到順直寬闊水庫垂向擴(kuò)散系數(shù)的計算公式為：

式（13）表明順直寬闊垂直岸水庫垂向擴(kuò)散系數(shù)與垂直面上岸邊污染混合區(qū)最大深度的2次方成正比，與平均流速的1次方成正比，與最大長度的1次方成反比。說明在平均流速相同的條件下，污染混合區(qū)最大長度越大垂向擴(kuò)散系數(shù)越小，最大深度和垂直面積Sz越大垂向擴(kuò)散系數(shù)越大。在實際應(yīng)用中，根據(jù)垂直岸水庫岸邊垂直面上污染混合區(qū)的現(xiàn)場觀測結(jié)果采用式（13）計算垂向擴(kuò)散系數(shù)。

3．2 岸邊污染混合區(qū)體積法

根據(jù)岸邊污染混合區(qū)體積計算垂直岸水庫的垂向擴(kuò)散系數(shù)是考慮到受天然水庫地形的影響，岸邊污染混合區(qū)的空間形狀一般不會像式（11）給出的那么規(guī)則，垂直面上岸邊污染混合區(qū)最大深度和面積的測量受多種因素的影響也會存在變數(shù)。為了減小垂直面上岸邊污染混合區(qū)最大深度和面積測量不精確給垂向擴(kuò)散系數(shù)計算帶來的誤差，可以根據(jù)現(xiàn)場觀測數(shù)據(jù)采用數(shù)值積分求和方法計算岸邊污染混合區(qū)（云團(tuán)）的體積，然后根據(jù)岸邊污染混合區(qū)體積進(jìn)行總體控制計算垂向擴(kuò)散系數(shù)。

由文獻(xiàn)［5］可知，垂直岸水庫三維問題岸邊污染混合區(qū)體積的理論計算公式為

式（15）表明順直寬闊垂直岸水庫垂向擴(kuò)散系數(shù)與平均流速和岸邊污染混合區(qū)體積的2次方成正比，與最大長度的4次方成反比，與側(cè)向擴(kuò)散系數(shù)的1次方成反比。說明在平均流速相同的條件下，岸邊污染混合區(qū)最大長度和側(cè)向擴(kuò)散系數(shù)越大垂向擴(kuò)散系數(shù)越小，污染混合區(qū)體積越大垂向擴(kuò)散系數(shù)越大；在側(cè)向擴(kuò)散系數(shù)相同的條件下，平均流速和體積越大垂向擴(kuò)散系數(shù)越大，而岸邊污染混合區(qū)最大長度越大垂向擴(kuò)散系數(shù)越小。

在實際應(yīng)用中，根據(jù)垂直岸水庫岸邊污染混合區(qū)的現(xiàn)場觀測結(jié)果采用式（15）計算垂向擴(kuò)散系數(shù)。觀測數(shù)據(jù)包括岸邊污染混合區(qū)范圍內(nèi)的平均流速、岸邊污染混合區(qū)的最大長度和體積以及由2．3節(jié)中計算方法觀測求得的側(cè)向擴(kuò)散系數(shù)。

3．3 水面橫向濃度積分法

當(dāng)z＝0時，由式（9）得到垂直岸水庫三維問題的岸邊水面等強(qiáng)度時間連續(xù)點源條件下的水面污染物濃度分布為

由式（10）和式（14）可以得到順直寬闊水庫垂向擴(kuò)散系數(shù)的計算公式為

將式（16）對y積分得到水面橫向積分濃度CwI為

由式（18）可以看出，水面橫向積分濃度公式中沒有出現(xiàn)側(cè)向擴(kuò)散系數(shù)Ey。這樣就可以根據(jù)垂直岸水庫岸邊排污口下游距離x0處水面橫向濃度分布的觀測資料求得水面橫向積分濃度CwI，然后利用下式計算垂向擴(kuò)散系數(shù)

式（19）表明順直寬闊垂直岸水庫垂向擴(kuò)散系數(shù)與平均流速和排污口下游距離的1次方成反比，與水面橫向積分濃度的2次方成反比，與排污強(qiáng)度的2次方成正比。說明在排污強(qiáng)度相同的條件下，水面橫向積分濃度越大垂向擴(kuò)散系數(shù)越小。即垂向擴(kuò)散系數(shù)越小，向深水區(qū)擴(kuò)散的污染物就會越少，大量污染物就會聚集在水面附近的淺層水體。對水庫水面中心排放條件式（19）同樣適用，此時水面橫向積分濃度CwI的積分區(qū)間是斷面全寬度。

根據(jù)水面橫向濃度積分法計算垂向擴(kuò)散系數(shù)，使現(xiàn)場觀測工作得到了大大地簡化，無需在深水區(qū)大量采樣觀測其空間濃度分布。觀測水面橫向濃度分布既省時省力，又容易保證觀測精度。建議水面橫向濃度分布的觀測選在污染混合區(qū)最大寬度出現(xiàn)的斷面附近，該斷面污染物濃度相對較高，污染物濃度分布較寬，有利于保證精度。

4 現(xiàn)場觀測結(jié)果分析

氵光府河發(fā)源于泰山山脈，流經(jīng)泰安、濟(jì)寧兩地市入南四湖的南陽湖北端，流域面積1 331 km2。作者于2005年3月枯水期調(diào)研時，在濟(jì)寧市河路大橋附近右岸發(fā)現(xiàn)某淀粉廠向府河大量排放超標(biāo)淀粉廢水，在排污口附近下游形成清晰可辨的白色污染混合區(qū)，當(dāng)即進(jìn)行了現(xiàn)場拍照和測量，圖2給出了該淀粉廠岸邊排污口污染混合區(qū)現(xiàn)場觀測示意圖。經(jīng)現(xiàn)場測量該淀粉廠岸邊高濃度白色污染混合區(qū)外邊界最大長度Ls和最大寬度bs分別為32 m和5 m，最大寬度相應(yīng)的縱坐標(biāo)Lc為12 m，污染混合區(qū)范圍內(nèi)的平均流速U為0．25 m／s，該河段水面寬度B為72 m。

氵光府河順直寬闊，而且縱向移流作用較強(qiáng)，排放岸邊高濃度白色污染混合區(qū)的色度基本不受對岸邊界反射的影響。將現(xiàn)場觀測的最大長度和最大寬度代入式（4）得到該淀粉廠岸邊高濃度白色污染混合區(qū)外邊界（等色度）曲線方程為

由式（20）繪制的污染混合區(qū)理論曲線見圖2。由圖2可以看出，由現(xiàn)場觀測的最大長度和最大寬度確定的污染混合區(qū)理論曲線與拍照的該淀粉廠岸邊高濃度白色污染混合區(qū)外邊界吻合良好，最大寬度相應(yīng)的縱坐標(biāo)理論值Lc＝Ls／e＝11．8 m與現(xiàn)場觀測值12 m非常接近，說明河流保守物質(zhì)平面二維移流擴(kuò)散簡化方程以及由此得到的理論公式適用于府河污染混合區(qū)的計算。因此，可采用式（5）計算府河側(cè)向擴(kuò)散系數(shù)為

圖2 氵光府河某淀粉廠岸邊排污口污染混合區(qū)現(xiàn)場觀測示意圖Fig．2 The observation of pollutantm ixing zone at some starchy factory’s outfall near Guangfu River bank

順直寬明渠混合污染物側(cè)向擴(kuò)散系數(shù)計算方法涉及的現(xiàn)場觀測數(shù)據(jù)主要是岸邊污染混合區(qū)外邊界最大長度、最大寬度和污染混合區(qū)范圍內(nèi)的平均流速共3個要素，對已有排污口的現(xiàn)場觀測在岸邊污染混合區(qū)外邊界線清晰可辨時更為簡單方便，快捷靈活，利于掌握。一般情況下，需進(jìn)行現(xiàn)場水流水質(zhì)同步監(jiān)測，即在已有排污口鄰近上游和下游污染影響區(qū)布設(shè)若干個斷面和垂線，監(jiān)測流速和水質(zhì)指標(biāo)，同時測定排污強(qiáng)度以及河流水文水力要素。然后，根據(jù)水質(zhì)監(jiān)測結(jié)果繪制一系列等濃度線并計算相應(yīng)的平均流速，在等濃度線與式（4）表示的標(biāo)準(zhǔn)曲線吻合較好時，采用式（5）計算側(cè)向擴(kuò)散系數(shù)，分析側(cè)向擴(kuò)散系數(shù)與河流水文水力要素和污染混合區(qū)尺度的關(guān)系。

5 結(jié) 論

（1）推導(dǎo)了順直寬闊明渠岸邊污染混合區(qū)外邊界標(biāo)準(zhǔn)曲線和曲面的統(tǒng)一方程，包含最大長度Ls、最大寬度bs和最大深度ds等參數(shù)，曲線形狀近似于半橢圓，空間形狀為近似橢球體的一部分，表明污染混合區(qū)具有相似性。

（2）給出了順直寬闊明渠混合污染物側(cè)向和垂向擴(kuò)散系數(shù)的理論計算公式，表明側(cè)向和垂向擴(kuò)散系數(shù)可以根據(jù)污染混合區(qū)特征尺度和平均流速來計算，還提出了采用水面橫向積分濃度確定垂向擴(kuò)散系數(shù)的實用方法。

（3）通過現(xiàn)場觀測結(jié)果分析，給出了氵光府河下游河段枯水期的側(cè)向擴(kuò)散系數(shù)Ey＝0．27 m2／s。

［1］ 張永良，李玉梁．排污混合區(qū)分析計算指南［M］．北京：海洋出版社，1993：3－9．（ZHANG Yong-liang，LI Yu-liang．Companion of Analytical Calculation to the Pol-lutant Mixing Zone in the Outfall’s［M］．Beijing：Ocean Press，1993：3－9．（in Chinese））

［2］ 武周虎，賈洪玉．河流污染混合區(qū)的解析計算方法［J］．水科學(xué)進(jìn)展，2009，20（4）：544－548．（WU Zhou-hu，JIA Hong-yu．Analytic method for pollutant mixing zone in river［J］．Advances in Water Science，2009，20（4）：544－548．（in Chinese））

［3］ 武周虎．明渠移流擴(kuò)散中無量綱數(shù)與相似準(zhǔn)則及其應(yīng)用［J］．青島理工大學(xué)學(xué)報，2008，29（5）：17－22．（WU Zhou-hu．The concept and its application of the dimen-sionless numbers＆conform rules by advection-diffusion in the open channel flow［J］．Journal of Qingdao Tech-nological University，2008，29（5）：17－22．（in Chi-nese））

［4］ 武周虎．傾斜岸水庫坡地形污染混合區(qū)的三維解析計算方法［J］．科技導(dǎo)報，2008，26（18）：30－34．（WU Zhou-hu．3-D Analytic computationalmethod of pollutant mixing zone for reservoir with incline bank［J］．Science＆Technology Review，2008，26（18）：30－34．（in Chi-nese））

［5］ 武周虎．水庫鉛垂岸地形污染混合區(qū)的三維解析計算方法［J］．西安理工大學(xué)學(xué)報，2009，25（4）：436－440．（WU Zhou-hu．3-D analytic computationalmethod of pol-lutantmixing zone for reservoir with vertical bank［J］．Journal of Xi’an University of Technology，2009，25（4）：436－440．（in Chinese））

［6］ 黃真理，李玉粱，李錦秀，等．三峽水庫水環(huán)境容量計算［J］．水利學(xué)報，2004，35（3）：7－16．（HUANG Zhen-li，LIYu-liang，LI Jin-xiu，et al．Water environ-mental capacity for the reservoir of Three Gorges Project［J］．Journal of Hydraulic Engineering，2004，35（3）：7－16．（in Chinese））

［7］ FISCHER H B，IMBERGER J，LIST E J，et al．Mixingin Inland and CoastalWaters［M］．New York：Academic Press，1979，107－112．

［8］ ELDER JW．The diffusion ofmarked fluid in turbulent shear flow［J］．Fluid Mechanics，1959，5（4）：544－560．

［9］ LAU Y L，KRISHNAPPAN B G．Transverse dispersion in rectangular channels［J］．Journal of Hydraulics Divi-sion，ASCE，1977，103（10）：1173－1189．

［10］HOLLEY E R，SIEMOUS J，ABRAHAM G．Some as-pects of analyzing transverse diffusion in rivers［J］．Jour-nal of Hydraulic Research，1972，10（1）：27－57．

［11］DEMETRACOPOULOS A C，STEFAN H G．Transverse mixing in wide and shallow river：case study［J］．Journal of Environmental Engineering，ASCE，1983，109（3）：685－697．

［12］周 云．天然河流橫向混合系數(shù)的研究［J］．水利學(xué)報，1988，19（6）：54－60．（ZHOU Yun．Study on trans-verse diffusion coefficient in natural river［J］．Journal of Hydraulic Engineering，1988，19（6）：54－60．（in Chi-nese））

［13］李 玲，李玉梁，陳嘉范．梯形斷面水槽中橫向擴(kuò)散系數(shù)的實驗［J］．清華大學(xué)學(xué)報（自然科學(xué)版），2003，43（8）：1124－1126，1152．（LI Ling，LI Yu-liang，CHEN Jia-fan．Experimental study on transverse diffusion coefficients in an open trapezoidal channel［J］．Journal of Tsinghua University（Science and Technology），2003，43（8）：1124－1126，1152．（in Chinese））

［14］郭建青，鄭鵬海．確定河流橫向擴(kuò)散系數(shù)的改進(jìn)矩法［J］．水電能源科學(xué)，1998，16（4）：32－35．（GUO Jian-qing，ZHENG Peng-hai．The improved moment method for determining the transverse diffusion coefficient of river stream［J］．International Journal Hydroelectric Energy，1998，16（4）：32－35．（in Chinese））

［15］梁秀娟，劉耀瑩，張文靜，等．室內(nèi)模擬試驗確定橫向擴(kuò)散系數(shù)的研究［J］．吉林大學(xué)學(xué)報（地球科學(xué)版），

2004，34（4）：560－565．（LIANG Xiu-juan，LIU Yao-ying，ZHANGWen-jing，et al．The study of determining the transverse diffusion coefficient of river through the in-door simulation experiments［J］．Journal of Jilin Univer-sity（Earth Seience Edition），2004，34（4）：560－565．（in Chinese））

［16］郭建青，張 勇．確定河流橫向擴(kuò)散系數(shù)的直線圖解法［J］．水利學(xué)報，1997，28（1）：62－67．（GUO Jian-qing，ZHANG Yong．Linear graphicsmethod for determi-ning the transverse diffusion coefficient of river［J］．Jour-nal of Hydraulic Engineering，1997，28（1）：62－67．（in Chinese））

［17］鄭鵬海，郭建青．確定河流橫向擴(kuò)散系數(shù)的線性回歸法［J］．水電能源科學(xué)，1999，17（3）：17－19．（ZHENG Peng-hai，GUO Jian-qing．The linear regression method for determining the transverse diffusion coefficientof river stream［J］．International Journal Hydroelectric Energy，1999，17（3）：17－19．（in Chinese））

［18］龍炳清，趙仕林，周后珍，等．岷江蕨溪段環(huán)境水力學(xué)研究——橫向擴(kuò)散［J］．四川師范大學(xué)學(xué)報（自然科學(xué)版），2001，24（4）：372－374．（LONG Bing-qing，ZHAO Shi-lin，ZHOU Hou-zhen，et al．The study of en-vironmental hydraulics for Juexi section in min river——transverse diffusion［J］．Journal of Sichuan Normal Uni-versity（Natural science），2001，24（4）：372－374．（in Chinese））

［19］金保明，楊曉華，金菊良，等．確定河流橫向擴(kuò)散系數(shù)的實碼遺傳算法［J］．水電能源科學(xué)，2000，18（1）：9－12．（JIN Bao-ming，YANG Xiao-hua，JIN Ju-liang，et al．Real coding genetic algorithm for determining the transverse diffusion coefficient of river［J］．International Journal Hydroelectric Energy，2000，18（1）：9－12．（in Chinese））

［20］龍騰銳，郭勁松，馮裕釗，等．二維水質(zhì)模型橫向擴(kuò)散系數(shù)的人工神經(jīng)網(wǎng)絡(luò)模擬［J］．重慶環(huán)境科學(xué)，2002，24（2）：25－28．（LONG Teng-rui，GUO Jing-song，F(xiàn)ENG Yu-zhao，et al．Modulus of transverse diffuse sim-ulation based on artificial neural network［J］．Environ-mental Science of Chongqing，2002，24（2）：25－28．（in Chinese））

［21］鄭旭榮，鄧志強(qiáng)，申繼紅．順直河流橫向紊動擴(kuò)散系數(shù)［J］．水科學(xué)進(jìn)展，2002，13（6）：670－674．（ZHENG Xu-rong，DENG Zhi-qiang，SHEN Ji-hong．Transverse turbulent dispersion coefficient in straight channel［J］．Advances in Water Science，2002，13（6）：670－674．（in Chinese））

（編輯：周曉雁）

Calculation M ethod of Lateral and Vertical Diffusion Coefficients for M ixed Pollutant in Open Channel Flow and Its App lication

WU Zhou-hu，HU De-jun，XU Mei-e
（Institute of Environmental Engineering，Qingdao Technological University，Qingdao 266033，Shandong）

Lateral and vertical diffusion coefficients of open channel is one of importantwater quality parameters to measure themixedly transporting capacity of pollutants near the bank flow．Whether the value is accurate or not is directly related to the reliability of thewater quality forecast in open channel flow．Based on the theoretical calcula-tionmethod of pollutantmixing zone，the unified equation of the standard curve and curved surface outside the bor-ders of pollution mixing zone is deduced，including the characteristic parameters of themaximum length Ls，maxi-mum width bsandmaximum depth dsetal．The shape of the curve approximates to a semi-ellipse，and the shape ofcurved surface is akin to a partof ellipsoid，which shows the pollution of themixing zone has similarity．The calcu-lation equation for determining the lateral or vertical diffusion coefficient by themaximum length，maximum width ormaximum depth of the outer boundary of pollutantmixing zone and themean velocity is given．The calculation method of the lateral or vertical diffusion coefficient which is overall controlled by the area or volume of pollutant mixing zone，and the practicalmethod of vertical diffusion coefficient which is determined by the concentration of surface transverse integration，are proposed．By analyzing the data of field observation，the lateral diffusion coeffi-cient of the downstream of Guangfu River in low water season is 0．27 m2／s．

open channel；pollutantmixing zone；standard curve；lateral diffusion coefficient；vertical diffusion coefficient； calculation methods

X522，X524

A

1001－5485（2010）010－0023－07

2010-08-23

國家自然科學(xué)基金項目（50979036），國家水體污染控制與治理科技重大專項（2009ZX07210－007－01）

武周虎（1959-），男，陜西岐山人，教授，主要從事環(huán)境水力學(xué)與水污染防治研究，（電話）13061205508（電子信箱）wu＿zh2008＠yahoo．com．cn。