排水管網(wǎng)檢查井死水區(qū)對(duì)溶質(zhì)縱向離散影響綜述

摘要：在排水管網(wǎng)中，下沉式檢查井的底部區(qū)域或者出現(xiàn)超載情況的平流式檢查井，均可為可溶性物質(zhì)提供滯留死水區(qū)，顯著改變?nèi)苜|(zhì)的輸運(yùn)及混合過(guò)程，對(duì)其縱向離散有重大影響。綜述了排水管網(wǎng)檢查井死水區(qū)內(nèi)溶質(zhì)混合擴(kuò)散過(guò)程的相關(guān)理論及模型研究進(jìn)展，并梳理了相關(guān)學(xué)者的研究成果，揭示了檢查井死水區(qū)對(duì)溶質(zhì)縱向離散的影響，詳細(xì)刻畫(huà)了超載檢查井內(nèi)死水區(qū)對(duì)溶質(zhì)混合擴(kuò)散的滯留效應(yīng)，量化了檢查井超載深度與直徑的比值，并根據(jù)不同的水文條件給出了兩個(gè)水力臨界值。針對(duì)當(dāng)前水質(zhì)模型忽略死水區(qū)滯留效應(yīng)而導(dǎo)致水質(zhì)預(yù)測(cè)存在誤差的現(xiàn)狀，提出了排水系統(tǒng)內(nèi)溶質(zhì)輸運(yùn)離散過(guò)程精確預(yù)測(cè)模型進(jìn)一步研究的重點(diǎn)方向。

關(guān) 鍵 詞：排水管網(wǎng)； 縱向離散； 溶質(zhì)； 死水區(qū)； 檢查井

中圖法分類號(hào)： TV134.2

文獻(xiàn)標(biāo)志碼： A

DOI：10.16232/j.cnki.1001-4179.2024.11.028

0 引 言

流動(dòng)中與主流相比速度為零或相對(duì)較低的區(qū)域定義為死水區(qū)，通常出現(xiàn)在積水區(qū)或再循環(huán)區(qū)。在天然河流中，由于幾何不規(guī)則性，死水區(qū)是一個(gè)相當(dāng)常見(jiàn)的特征^［1-3］。排水系統(tǒng)中同樣也存在死水區(qū)，例如檢查井內(nèi)幾何形狀的顯著變化為雨污水的動(dòng)態(tài)存儲(chǔ)提供了額外的體積^［4］，該區(qū)域可以稱為檢查井的死水區(qū)。研究表明，檢查井內(nèi)的溶質(zhì)混合是一個(gè)非高斯過(guò)程，排水系統(tǒng)中的溶質(zhì)隨管道向下游擴(kuò)散，經(jīng)過(guò)檢查井時(shí)死水區(qū)能把一部分溶質(zhì)滯留下來(lái)，當(dāng)溶質(zhì)主體隨水流下移后，死水區(qū)中的溶質(zhì)又會(huì)慢慢釋放出來(lái)。檢查井主要設(shè)置在流向、管道坡度、直徑和標(biāo)高等發(fā)生變化的位置，主要功能是方便進(jìn)入管網(wǎng)進(jìn)行檢查、清潔、維護(hù)和維修，其死水區(qū)對(duì)溶質(zhì)混合、輸運(yùn)綜合影響較為顯著，對(duì)于排水系統(tǒng)內(nèi)物理過(guò)程的精確預(yù)測(cè)、溯源反演也至關(guān)重要。Valentine等^［5］指出有兩個(gè)變量控制著死水區(qū)在離散過(guò)程中的作用，即流量、檢查井直徑與進(jìn)水管直徑的比率，并認(rèn)為為了更好地描述排水系統(tǒng)的水質(zhì)變化，有必要了解檢查井內(nèi)死水區(qū)對(duì)縱向離散系數(shù)的影響。

國(guó)外學(xué)者最先探究天然河流中死水區(qū)對(duì)溶質(zhì)擴(kuò)散的影響機(jī)理。Guymer等^［6］對(duì)檢查井中的溶質(zhì)對(duì)流擴(kuò)散過(guò)程進(jìn)行了較為詳細(xì)的研究，此后陸續(xù)有學(xué)者開(kāi)始探究排水系統(tǒng)檢查井內(nèi)死水區(qū)對(duì)溶質(zhì)擴(kuò)散的影響^［7-8］。國(guó)內(nèi)關(guān)于死水區(qū)對(duì)溶質(zhì)的存儲(chǔ)、滯留特征方面影響的研究還較為鮮見(jiàn)，多數(shù)以河流中具有暫時(shí)存儲(chǔ)區(qū)的潛流為研究對(duì)象進(jìn)行分析^［9-11］。尚無(wú)研究從檢查井內(nèi)死水區(qū)對(duì)縱向離散系數(shù)影響的角度揭示內(nèi)在作用機(jī)制。為此，本文就這一研究領(lǐng)域進(jìn)行綜述，梳理總結(jié)相關(guān)挑戰(zhàn)性問(wèn)題并指出下一步研究方向。

1 死水區(qū)對(duì)溶質(zhì)縱向離散影響的研究進(jìn)展

1.1 死水區(qū)對(duì)示蹤濃度時(shí)間分布的影響

傳統(tǒng)上，溶質(zhì)縱向離散濃度數(shù)據(jù)根據(jù)一維對(duì)流擴(kuò)散方程^［12］（Advection Dispersion Equation，ADE）獲得：

Ct+uCx=K2Cx2（1）

式中：C為主流區(qū)溶質(zhì)濃度，t為時(shí)間，x為擴(kuò)散距離，u為斷面平均流速，K為縱向離散系數(shù)。

許多研究人員對(duì)溶質(zhì)的縱向離散進(jìn)行了大量現(xiàn)場(chǎng)示蹤試驗(yàn)，發(fā)現(xiàn)實(shí)際測(cè)量的濃度-時(shí)間分布曲線一般都是偏態(tài)的^［13-14］，這種偏斜式的濃度-時(shí)間分布不能基于傳統(tǒng)ADE方程獲得，如圖1所示。Fischer等^［15］認(rèn)為，現(xiàn)場(chǎng)示蹤試驗(yàn)測(cè)出的濃度隨時(shí)間變化的曲線偏態(tài)分布，是天然河流中死水區(qū)的滯留作用導(dǎo)致的。

河床上的洼區(qū)、河岸的凹凸處以及其他建筑物附近，將河流分為具有速度剖面的流動(dòng)區(qū)和速度沿著底部或側(cè)面的停滯區(qū)，此停滯區(qū)即為死水區(qū)^［16］。Beer等^［17］使用參數(shù)估計(jì)的時(shí)間序列方法模擬天然河流中非保守污染物濃度的擴(kuò)散過(guò)程，并使用來(lái)自各種河流的實(shí)測(cè)數(shù)據(jù)來(lái)證明此模型對(duì)濃度分布偏態(tài)的解釋優(yōu)于ADE方程。但自20世紀(jì)60年代末開(kāi)始，學(xué)者意識(shí)到在許多情況下傳統(tǒng)的ADE方程^［18］無(wú)法描述觀測(cè)到的示蹤濃度時(shí)間分布，于是為了真實(shí)地反映天然河流的縱向離散過(guò)程，提出了許多修正模型，包括瞬態(tài)存儲(chǔ)（Transient Storage，TS）模型以及死水區(qū)模型（Aggregated Dead Zone，ADZ）等，且發(fā)現(xiàn)這些修正模型中主流區(qū)的縱向離散系數(shù)比傳統(tǒng)一維對(duì)流擴(kuò)散方程中的小，并證明了死水區(qū)對(duì)污染物的滯留作用，以及死水區(qū)的存在引起的長(zhǎng)尾效應(yīng)。

1.2 用于死水區(qū)模擬的瞬時(shí)存儲(chǔ)模型和OTIS模型

為了描述包含死水區(qū)的水體中溶質(zhì)混合擴(kuò)散過(guò)程，Thackston等^［19］提出將瞬態(tài)存儲(chǔ)（TS）模型納入ADE模型中，以考慮死水區(qū)的影響。TS這種分布式物理模型可用下面的方程組表示：

Ct+uCx=K2Cx2+KsCd-C（2）

Cdt=KsAAtCd-C（3）

式中：A是主流的平均橫截面積，At是死水區(qū)的橫截面積，Cd是死水區(qū)溶質(zhì)濃度，Ks是存儲(chǔ)交換系數(shù)，K是縱向離散系數(shù)。方程（2）左側(cè)的兩項(xiàng)和右側(cè)第一項(xiàng)形成ADE模型的基礎(chǔ)，方程（2）右側(cè)的第二項(xiàng)和方程（3）是描述死水區(qū)中混合過(guò)程的元素。但TS模型沒(méi)有考慮用于超載檢查井的縱向離散分析，也沒(méi)有將模型中的參數(shù)值與排水系統(tǒng)的物理特性聯(lián)系起來(lái)。

其他類似利用TS模型對(duì)ADE模型進(jìn)行修正的技術(shù)，如Runkel^［20］提出的OTIS（One-dimensional Transport with Inflow and Storage）模型，已被開(kāi)發(fā)用于分析河道中暫態(tài)存儲(chǔ)作用的影響。Becker等^［21］用NaCl作示蹤劑并結(jié)合OTIS模型比較了天然河流與常規(guī)城市河道的潛流交換規(guī)律，發(fā)現(xiàn)天然河流與一般城市河道相比具有較強(qiáng)的潛流交換能力和瞬時(shí)存儲(chǔ)能力。近些年中國(guó)也有學(xué)者利用OTIS模型探究天然河流中的死水區(qū)對(duì)溶質(zhì)輸移擴(kuò)散的滯留效應(yīng)^［22-23］。

1.3 死水區(qū)溶質(zhì)擴(kuò)散模型

假設(shè)死水區(qū)中溶質(zhì)是完全混合的，天然河道的主流區(qū)與死水區(qū)存在溶質(zhì)交換，Pederson等^［24］假設(shè)死水區(qū)中的溶質(zhì)濃度（用Cd表示）與Cd和主流溶質(zhì)濃度C之間的差值成比例的速率變化：

Cdt=T^-1C-Cd（4）

因此，死水區(qū)模型（ADZ）可用下面方程表示：

Ct+uCx-K2Cx2=εT^-1Cd-C（5）

式中：ε為單位長(zhǎng)度死水區(qū)體積與單位長(zhǎng)度主流區(qū)體積之比，C為主流區(qū)的溶質(zhì)濃度，T為溶質(zhì)在死水區(qū)中滯留的時(shí)間。Davis等^［25］利用死水區(qū)（ADZ）模型，證明了示蹤劑在某條天然河流中的擴(kuò)散最初是受剪切流影響，但在示蹤云團(tuán)的速度減慢至與主流區(qū)域平均斷面流速相同的一段時(shí)間后，死水區(qū)對(duì)示蹤云團(tuán)縱向離散的影響會(huì)變得非常明顯。

2 檢查井死水區(qū)對(duì)溶質(zhì)縱向離散的影響

由于排水系統(tǒng)和天然河流類似，故一直以來(lái)在模擬排水管網(wǎng)水質(zhì)變化時(shí)多數(shù)借鑒河流水質(zhì)模型。在下水道系統(tǒng)的檢查井中同樣存在具有滯留效應(yīng)的死水區(qū)，為了更加準(zhǔn)確地描述系統(tǒng)中溶質(zhì)的對(duì)流擴(kuò)散過(guò)程，許多學(xué)者借鑒并利用河流死水區(qū)的研究成果對(duì)檢查井內(nèi)死水區(qū)開(kāi)展了大量試驗(yàn)研究。在天然河流中，死水區(qū)或儲(chǔ)存區(qū)的邊界形狀是無(wú)規(guī)則的，但與天然河流截然相反的是，檢查井邊界規(guī)則且復(fù)雜，同時(shí)具有明顯的橫向限制。

2.1 檢查井死水區(qū)中的溶質(zhì)擴(kuò)散

2.1.1 檢查井超載現(xiàn)象及影響分析

檢查井超載是指當(dāng)檢查井所連接的管段都是滿流狀態(tài)，或者節(jié)點(diǎn)水位超過(guò)了所有與之相連管道的最高內(nèi)頂時(shí)所出現(xiàn)的現(xiàn)象^［26］。排水系統(tǒng)在超載條件下管道滿流相對(duì)頻繁，這是因?yàn)榻涤臧l(fā)生后，地表徑流會(huì)通過(guò)井口進(jìn)入排水管網(wǎng)中，但是整個(gè)區(qū)域排水管網(wǎng)十分有限，當(dāng)排入的雨水超過(guò)排水管網(wǎng)所能承受的最大負(fù)荷容量時(shí)，雨水就會(huì)從檢查井溢出。排水系統(tǒng)出現(xiàn)超載的原因有兩個(gè)：① 當(dāng)管道沿線的排水量大于管道坡度下的剛好滿容量時(shí)，就會(huì)發(fā)生超載，如圖2（a）所示；② 當(dāng)下游排水組件產(chǎn)生回水效應(yīng)時(shí)進(jìn)而產(chǎn)生一定范圍的內(nèi)澇積水問(wèn)題，如圖2（b）所示。在滿管流的情況下，檢查井中的水位上升到入口管道的拱腹水位以上，可以看到一個(gè)充滿水的檢查井包括兩個(gè)部分：① 允許水流過(guò)的部分；② 位于拱腹上方的儲(chǔ)藏室。儲(chǔ)藏室通過(guò)提供儲(chǔ)物空間為排水系統(tǒng)提供了額外的容量^［18］，也就是死水區(qū)。

2.1.2 下沉式檢查井死水區(qū)對(duì)污染擴(kuò)散影響

下沉式檢查井是雨水排放系統(tǒng)的重要組成部分，當(dāng)污染物隨雨水從管段進(jìn)入檢查井后，會(huì)在井中形成一股旋流，一部分污染物隨主流向下游擴(kuò)散，一部分則在水深較深的區(qū)域充分混合，此時(shí)在檢查井靠近墻壁的部分會(huì)形成一個(gè)流速很低的區(qū)域，被稱為死水區(qū)^［27］，如圖3所示。由于集水井的深度通常大于其直徑，入流口和出流口對(duì)齊，以及集水井中具有滯留效應(yīng)的死水區(qū)，這都給去除污染物提供了機(jī)會(huì)^［28］。但學(xué)者一般更加關(guān)心檢查井清除雨水中顆粒與沉積物的有效性^［29］，而忽略了對(duì)溶質(zhì)縱向離散的具體影響。

2.2 超載檢查井內(nèi)死水區(qū)的滯留效應(yīng)

一些水質(zhì)模擬模型，如SWMM和InfoWorks，僅采用傳統(tǒng)的一維對(duì)流擴(kuò)散方程模擬溶質(zhì)在排水系統(tǒng)中的濃度變化，它們都假設(shè)溶質(zhì)在檢查井內(nèi)是完全混合的，而忽略了排水系統(tǒng)中檢查井造成的滯留效應(yīng)。Stovin等^［30］根據(jù)檢查井內(nèi)的超載深度和檢查井幾何形狀進(jìn)行示蹤試驗(yàn)，發(fā)現(xiàn)了復(fù)雜的混合條件，在這種情況下，完全混合的假設(shè)是不正確的，這將導(dǎo)致對(duì)排水系統(tǒng)污染物濃度的錯(cuò)誤估計(jì)。與無(wú)死水區(qū)溶質(zhì)的一維對(duì)流擴(kuò)散相反，死水區(qū)的滯留作用可以立即使通過(guò)死水區(qū)的溶質(zhì)峰值濃度顯著衰減。圖4是一個(gè)充滿水的檢查井內(nèi)死水區(qū)的滯留效應(yīng)示例，在時(shí)間t=t0時(shí)，示蹤劑由于平流而縱向離散；在時(shí)間t=t1時(shí)，一些示蹤劑被滯留在檢查井的死水區(qū)中，并且在儲(chǔ)存區(qū)域發(fā)生某種形式的混合；隨著大部分示蹤劑通過(guò)檢查井，在時(shí)間t=t2時(shí)，被困在儲(chǔ)存空間中的示蹤劑逐漸釋放，由于釋放速率較慢，溶質(zhì)的擴(kuò)散明顯延長(zhǎng)^［31］。

Guymer等^［32］研究發(fā)現(xiàn)，與沒(méi)有檢查井的直管相比，沿管道增加檢查井使污染物平均經(jīng)過(guò)時(shí)間增加了約2倍，超載檢查井的存在使縱向離散系數(shù)增大了5倍。Lau等^［33］通過(guò)室內(nèi)實(shí)驗(yàn)測(cè)量水頭損失和溶質(zhì)濃度分布，以量化檢查井井口直徑和超載水平變化對(duì)溶質(zhì)縱向分散的影響，結(jié)果表明縱向離散系數(shù)會(huì)隨著流量的增加而增加，且隨著檢查井井口直徑的增大普遍呈增大趨勢(shì)。

2.3 死水區(qū)模型在檢查井中溶質(zhì)分布預(yù)測(cè)的應(yīng)用

Dennis^［34］將死水區(qū)（ADZ）模型應(yīng)用在超載檢查井中，證明了ADZ模型在預(yù)測(cè)下游濃度分布曲線方面明顯優(yōu)于傳統(tǒng)一維對(duì)流擴(kuò)散模型（圖5），但同時(shí)發(fā)現(xiàn)ADZ模型仍然無(wú)法準(zhǔn)確預(yù)測(cè)濃度值的峰值。研究人員在研究中還發(fā)現(xiàn)，溶質(zhì)進(jìn)入檢查井結(jié)構(gòu)內(nèi)的儲(chǔ)存區(qū)域并旋轉(zhuǎn)，隨后滯留的溶質(zhì)隨主流釋放出來(lái)，流態(tài)形成長(zhǎng)尾。

Boxall等^［35］開(kāi)發(fā)了一種模擬下水道中水質(zhì)變化的模型，該模型利用ADZ模型來(lái)描述下水道系統(tǒng)中溶質(zhì)的輸運(yùn)擴(kuò)散過(guò)程，探索了檢查井內(nèi)的溶質(zhì)混合離散的內(nèi)在機(jī)理。Young等^［36］在實(shí)驗(yàn)中發(fā)現(xiàn)，ADZ模型預(yù)測(cè)的峰值濃度與實(shí)際測(cè)量的峰值濃度相似，但到達(dá)時(shí)間有一些延遲，并且由于ADZ模型需要確定的參數(shù)較多，在實(shí)際應(yīng)用中受到了一定限制。

2.4 影響溶質(zhì)在檢查井中縱向離散的因素

Scott^［37］提出通過(guò)記錄溶質(zhì)滯留時(shí)間分布（Retention Time Distribution，RTD）和累計(jì)滯留時(shí)間分布（Cumulative Retention Time Distribution，CRTD）來(lái)揭示流場(chǎng)的基本混合特性。RTD及CRTD可以通過(guò)物理實(shí)驗(yàn)或流體動(dòng)力學(xué)模擬計(jì)算來(lái)獲得，CRTD表示RTD的積分形式，通常需將其歸一化處理。Lau等^［33］將3種模型（ADE、ADZ、CRTD）與實(shí)驗(yàn)室測(cè)量結(jié)果進(jìn)行比較，發(fā)現(xiàn)與ADE和ADZ模擬結(jié)果相比，CRTD模擬的結(jié)果與實(shí)驗(yàn)室結(jié)果更為吻合，可以看出在特定工況下CRTD模型優(yōu)于ADE和ADZ模型。郭長(zhǎng)強(qiáng)等^［38］通過(guò)示蹤試驗(yàn)得到了水力CRTD曲線，計(jì)算相關(guān)濕地水力性能指標(biāo)，分析了濕地水深對(duì)濕地水力性能的影響，指出濕地水力性能越低，越容易造成水流短路和濕地出現(xiàn)死水區(qū)的情況。

利用兩個(gè)歸一化的CRTD模型研究超載檢查井中溶質(zhì)的離散特性，發(fā)現(xiàn)它們均受到檢查井幾何形態(tài)的影響。此外，這些特性可通過(guò)檢查井內(nèi)部超載深度與直徑的比值來(lái)量化描述^［39-40］。Sonnenwald等^［41］研究發(fā)現(xiàn)，對(duì)于D_m/D_pgt;4.4（Dm和Dp分別為檢查井及進(jìn)水管的直徑）的檢查井，存在臨界超載深度。Lau等^［33］通過(guò)實(shí)驗(yàn)發(fā)現(xiàn)，平流檢查井的CRTD曲線斜率在超載深度閾值后出現(xiàn)了急劇轉(zhuǎn)變，當(dāng)檢查井超載深度S≥（2～2.5）D_m時(shí)，水流不再與檢查井下方的區(qū)域發(fā)生較劇烈的溶質(zhì)交換，而且從出水口直接通過(guò)，流態(tài)接近于平推流，此時(shí)的CRTD曲線急速上升，檢查井的上部分出現(xiàn)死水區(qū)溶質(zhì)擴(kuò)散短路現(xiàn)象，如圖6所示。當(dāng)流量達(dá)到一定程度，檢查井內(nèi)這種流場(chǎng)變化已經(jīng)形成，此時(shí)流量即使繼續(xù)增大，整個(gè)流場(chǎng)流態(tài)也不會(huì)發(fā)生顯著變化。

Stovin等^［42］的研究認(rèn)為，完全混合流動(dòng)的臨界超載深度S=0.258D_m，研究中發(fā)現(xiàn)當(dāng)檢查井的水深在臨界超載深度以下時(shí)，CRTD曲線趨于水平，此時(shí)檢查井內(nèi)部的流場(chǎng)為完全混合狀態(tài)。而當(dāng)水深在水力臨界值以上且增大時(shí)，得到的CRTD曲線斜率將逐漸上升（圖7），說(shuō)明檢查井內(nèi)部出現(xiàn)了返混程度相較于完全混合狀態(tài)減輕的趨勢(shì)，一部分溶質(zhì)會(huì)隨著主流的沖擊帶到下游管道進(jìn)水口；剩下一部分則被滯留在檢查井死水區(qū)內(nèi)，經(jīng)返混后流出檢查井。

3 尚待解決的問(wèn)題

大量研究表明，檢查井的幾何特性影響死水區(qū)的滯留效應(yīng)，而滯留效應(yīng)對(duì)溶質(zhì)在排水系統(tǒng)中的輸移擴(kuò)散產(chǎn)生顯著影響。目前尚未解決的問(wèn)題包括以下幾個(gè)方面：

（1） 當(dāng)前提出的修正模型需要確定的參數(shù)較多，沒(méi)有傳統(tǒng)一維對(duì)流擴(kuò)散方程簡(jiǎn)便，以至于在實(shí)際應(yīng)用中受到諸多限制，亟需構(gòu)建一個(gè)既簡(jiǎn)便又能準(zhǔn)確考慮死水區(qū)滯留效應(yīng)的模型。

（2） 之前的物理實(shí)驗(yàn)和數(shù)值模擬研究中一般僅考慮了穩(wěn)態(tài)條件，應(yīng)進(jìn)一步研究非定常流時(shí)檢查井死水區(qū)對(duì)溶質(zhì)縱向離散的影響。

（3） 目前研究成果僅僅應(yīng)用在了簡(jiǎn)單的圓形平流式檢查井中，應(yīng)繼續(xù)探索在其他幾何形狀的檢查井中死水區(qū)如何影響溶質(zhì)的離散，例如下沉式、矩形、“T”形等檢查井。

（4） SWMM及InfoWorks等現(xiàn)有水質(zhì)模型完全忽略了檢查井對(duì)溶質(zhì)擴(kuò)散的影響，因此利用該水質(zhì)模型獲得的污染物濃度數(shù)據(jù)難以準(zhǔn)確評(píng)估水質(zhì)真實(shí)變化。因此，亟需進(jìn)一步開(kāi)發(fā)針對(duì)檢查井內(nèi)水質(zhì)變化的模塊，需要更加精確地復(fù)現(xiàn)污染物在排水系統(tǒng)中的輸運(yùn)過(guò)程，提高污染溯源的可靠性，并根據(jù)不同的水文條件在排水系統(tǒng)內(nèi)分段建立水質(zhì)模型。

4 結(jié) 語(yǔ)

大量針對(duì)檢查井死水區(qū)的研究表明其存在滯留效應(yīng)，使流態(tài)形成長(zhǎng)尾，增大了離散系數(shù)。針對(duì)傳統(tǒng)一維對(duì)流擴(kuò)散方程的自身局限性，學(xué)者們提出多種修正模型，包括瞬時(shí)存儲(chǔ)模型、死水區(qū)模型、累計(jì)滯留時(shí)間分布模型等，揭示了檢查井死水區(qū)內(nèi)溶質(zhì)混合擴(kuò)散的流動(dòng)特性。另外還發(fā)現(xiàn)，檢查井的水深超載深度S與檢查井直徑D_m的比值影響溶質(zhì)在檢查井內(nèi)的縱向離散，導(dǎo)致排水管網(wǎng)系統(tǒng)內(nèi)出現(xiàn)多種水文環(huán)境。為了更加精確地復(fù)現(xiàn)污染物在排水系統(tǒng)中的離散過(guò)程，提高污染溯源的可靠性，針對(duì)不同幾何特征和水文條件的檢查井，建立更加簡(jiǎn)便且精確的死水區(qū)模型是未來(lái)的一個(gè)重點(diǎn)發(fā)展方向。

參考文獻(xiàn)：

［1］ 余常昭.環(huán)境流體力學(xué)導(dǎo)論［M］.北京：清華大學(xué)出版社，1992：176-186.

［2］ LIN Y C，CHANG M S，MEDINA J M A.A methodology for solute transport in unsteady，nonuniform streamflow with subsurface interaction［J］.Advances in Water Resources，2005，28（8）：871-883.

［3］ 薛宇宙，丁于墨，儲(chǔ)家越.黑臭水體生物與化學(xué)治理技術(shù)研究進(jìn)展［J］.水利水電快報(bào)，2023，44（12）：94-100.

［4］ LAU S，STOVIN V，GUYMER I.The prediction of solute transport in surcharged manholes using CFD［J］.Water Science amp; Technology，2007，55（4）：57-64.

［5］ VALENTINE E M，WOOD I R.Longitudinal dispersion with dead zones［J］.Journal of the Hydraulics Division，1977，103（9）：975-990.

［6］ GUYMER I，OBRIEN R.The effects of surcharged manholes on the travel time and dispersion of solutes in sewer systems［J］.Water Science amp; Technology，1995，31（7）：51-59.

［7］ GUYMER I，DENNIS P，O’BRIEN R，et al.Diameter and surcharge effects on solute transport across surcharged manholes［J］.Journal of Hydraulic Engineering，2005，131（4）：312-321.

［8］ HOWARTH D A.The hydraulic performance of scale model storm sewer junctions［D］.Manchester：University of Manchester，1985.

［9］ 于靖，張華.城市小型河流水動(dòng)力彌散和潛流交換過(guò)程［J］.水科學(xué)進(jìn)展，2015，26（5）：714-721.

［10］金光球，李凌.河流中潛流交換研究進(jìn)展［J］.水科學(xué)進(jìn)展，2008，19（2）：285-293.

［11］PURNAMA A.The effect of dead zones on longitudinal dispersion in streams［J］.Journal of Fluid Mechanics，1988，186：351-377.

［12］TAYLOR G.The dispersion of matter in turbulent flow through a pipe［J］.Proceedings of the Royal Society of London，1954，223（1155）：446-468.

［13］顧莉，華祖林.天然河流縱向離散系數(shù)確定方法的研究進(jìn)展［J］.水利水電科技進(jìn)展，2007，27（2）：85-89.

［14］BENCALA K E.Simulation of solute transport in a mountain pool-and-riffle stream with a kinetic mass transfer model for sorption［J］.Water Resources Research，1983，19（3）：732-738.

［15］FISCHER H B，LIST E J，KOH R C Y，et al.Mixing in inland and coastal waters［M］.Amsterdam：Academic Press，1979.

［16］CZERNUSZENKO W，ROWISKI P M，SUKHODOLOV A.Experimental and numerical validation of the dead-zone model for longitudinal dispersion in rivers［J］.Journal of Hydraulic Research，1998，36（2）：269-280.

［17］BEER T，YOUNG P C.Longitudinal dispersion in natural streams［J］.Journal of Environmental Engineering，1983，109（5）：1049-1067.

［18］GUYMER I，O′BRIEN R，MARK O，et al.An investigation of fine sediment mixing within free-flowing and surcharged manholes［J］.Water Science amp; Technology，1998，37（1）：215-222.

［19］THACKSTON E L，SCHNELLE J K B.Predicting effects of dead zones on stream mixing［J］.Journal of the Sanitary Engineering Division，1970，96（2）：319-331.

［20］RUNKEL R L.One-dimensional transport with inflow and storage （OTIS）：a solute transport model for streams and rivers［M］.Virginia：US Geological Survey，1998.

［21］BECKER J F，ENDRENY T A，ROBINSON J D.Natural channel design impacts on reach-scale transient storage［J］.Ecological Engineering，2013，57：380-392.

［22］李如忠，楊源，丁貴珍，等.基于OTIS模型的巢湖十五里河源頭段參數(shù)靈敏性分析［J］.環(huán)境科學(xué)研究，2014，27（11）：1265-1271.

［23］李如忠，丁貴珍，徐晶晶，等.巢湖十五里河源頭段暫態(tài)儲(chǔ)存特征分析［J］.水利學(xué)報(bào)，2014，45（6）：631-640.

［24］PEDERSON F.Prediction of longitudinal dispersion in natural streams［R］.Lyngby：Technical University of Denmark，1977.

［25］DAVIS P，ATKINSON T.Longitudinal dispersion in natural channels：3.An aggregated dead zone model applied to the River Severn，UK［J］.Hydrology and Earth System Sciences，2000，4（3）：373-381.

［26］BEG M N A，CARVALHO R F，LEANDRO J.Effect of surcharge on gully-manhole flow［J］.Journal of Hydro-environment Research，2018，32：224-236.

［27］楊浩銘.排水系統(tǒng)中流動(dòng)特性及其與沉積物相互作用的數(shù)值模擬［D］.杭州：浙江大學(xué)，2018.

［28］HOWARD A K，MOHSENI O，GULLIVER J S，et al.Hydraulic analysis of suspended sediment removal from storm water in a standard sump［J］.Journal of Hydraulic Engineering，2012，138（6）：491-502.

［29］YANG Y，HUANG B，ZHU D Z.Experimental study of sediment washout from stormwater sumps［J］.Water Science amp; Technology，2022，86（9）：2454-2464.

［30］STOVIN V，BENNETT P，GUYMER I.Absence of a hydraulic threshold in small-diameter surcharged manholes［J］.Journal of Hydraulic Engineering，2013，139（9）：984-994.

［31］LAU S，DOUGLAS S T.Scaling dispersion processes in surcharged manholes［D］.Sheffield：University of Sheffield，2007.

［32］GUYMER I，O′BRIEN R.Longitudinal dispersion due to surcharged manhole［J］.Journal of Hydraulic Engineering，2000，126（2）：137-149.

［33］LAU S，STOVIN V，GUYMER I.Scaling the solute transport characteristics of a surcharged manhole［J］.Urban Water Journal，2008，5（1）：33-42.

［34］DENNIS P.Longitudinal dispersion due to surcharged manholes［D］.Sheffield：University of Sheffield，2000.

［35］BOXALL J，SHEPHERD W，GUYMER I，et al.Changes in water quality parameters due to in-sewer processes［J］.Water Science amp; Technology，2003，47（7/8）：343-350.

［36］YOUNG P，JAKEMAN A，MCMURTRIE R.An instrumental variable method for model order identification［J］.Automatica，1980，16（3）：281-294.

［37］SCOTT F H.Elements of chemical reaction engineering［M］.New York：Pearson Educacion，1999.

［38］郭長(zhǎng)強(qiáng)，董斌，劉俊杰，等.水深對(duì)塘堰濕地水力性能的影響［J］.應(yīng)用生態(tài)學(xué)報(bào)，2014，25（11）：3287-3295.

［39］STOVIN V R，GUYMER I，CHAPPELL M J，et al.The use of deconvolution techniques to identify the fundamental mixing characteristics of urban drainage structures［J］.Water Science amp; Technology，2010，61（8）：2075-2081.

［40］GUYMER I，STOVIN V.One-dimensional mixing model for surcharged manholes［J］.Journal of Hydraulic Engineering，2011，137（10）：1160-1172.

［41］SONNENWALD F，MARK O，STOVIN V，et al.Predicting manhole mixing using a compartmental model［J］.Journal of Hydraulic Engineering，2021，147（12）：04021046.

［42］STOVIN V，GUYMER I，LAU S T D.Dimensionless method to characterize the mixing effects of surcharged manholes［J］.Journal of Hydraulic Engineering，2010，136（5）：318-327.

（編輯：胡旭東）

Review on influence of dead water zone in inspection shaft of drainage network on pollutants longitudinal dispersion

HUANG Biao¹，YANG Liyuan¹，LIU Jiachun¹，YANG Yunhua²，CAI Haibo²

（1.School of Civil amp; Environmental Engineering and Geography Science，Ningbo University，Ningbo 315211，China； 2.Ningbo Municipal Engineering Construction Group Co.，Ltd.，Ningbo 315000，China）

Abstract：

In drainage networks，the bottom area of the sinking inspection shaft or the horizontal flow type inspection shaft under overload condition，can provide a stagnant dead water zone for the soluble substance，which significantly changes the transporting and mixing process of the pollutants，and has a significant impact on its longitudinal dispersion.In this paper，we reviewed the related theories and model research progress on pollutants mixing and diffusion process in the dead water area of inspection shafts in drainage networks，and the research results of relevant scholars were sorted out.The influence of the dead water area of an inspection shaft on longitudinal dispersion of the pollutants was revealed.The retention effect of the dead water zone in the overload inspection shaft on the pollutants mixing and diffusion was described in detail.The ratio of depth to diameter of an inspection well was quantified，and two hydraulic thresholds were given according to different hydrological conditions.In view of the current situation that the water quality models ignore the retention effect of dead water zone，which lead to water quality prediction error，key research directions on accurate prediction models for simulating pollutants transport dispersion process in drainage network are put forward.

Key words：

drainage network； longitudinal dispersion； pollutants； dead water area； inspection shaft