城市排水管網(wǎng)外水診斷中特征因子的選擇

姚 俊，王 慶，2

（1.中電建生態(tài)環(huán)境集團有限公司，廣東 深圳 518101；2.中電建水環(huán)境科技有限公司，廣東 深圳 518102）

城市排水管網(wǎng)的外水診斷是當前城市排水系統(tǒng)提質(zhì)增效項目中重點工作之一，主要包括污水管網(wǎng)中地下水等入滲和雨水管網(wǎng)混錯接診斷。污水管網(wǎng)中地下水等外水入流入滲會引起污水廠進廠水量增加和水質(zhì)濃度下降，進廠水量超過污水廠處理量時會產(chǎn)生污水溢流污染，進廠水質(zhì)濃度下降會降低污水廠處理效率[1]。雨水管網(wǎng)中混錯接會導致污水直排或溢流進入河道[2]。外水量評估是排水管網(wǎng)外水診斷的核心內(nèi)容之一。

目前，排水管網(wǎng)外水量評估方法主要分為物理探測、水量平衡法、化學質(zhì)量平衡法3 類，其中物理探測主要是利用潛望鏡（QV）、管道內(nèi)窺檢測（CCTV）等檢測設(shè)備定位到管網(wǎng)外水入流入滲點[3-5]，然后進行水質(zhì)水量檢測，該方法成本較高；水量平衡法[6-7]是通過對比上下游流量變化，判斷進入該段管道的外水量，該方法無法判斷外水種類和比例；化學質(zhì)量平衡法則是通過上下游水質(zhì)指標濃度判斷該段管道的外水類型和比例，如特征因子法[8-9]和同位素法[10]等，其檢測成本低，能評估各種外水的來水量。

特征因子法等化學質(zhì)量平衡法應(yīng)用的關(guān)鍵是要合理選擇特征因子。不同外水中水質(zhì)指標差異較大，通常選擇氨氮、表面活性劑、安賽蜜和大腸桿菌等指標作為生活污水的特征因子，如徐祖信等[2]選擇安賽蜜來診斷上海市某雨水管網(wǎng)；尹海龍等[11]選擇一種微生物指示菌作為特征因子；Field 等[12]推薦氨氮、表面活性劑等作為特征因子。地下水通常選擇總硬度作為特征因子[13-14]。由于工廠生產(chǎn)工藝形式多樣，需結(jié)合具體情況選擇特征因子，通常電子企業(yè)的工業(yè)污水選擇氟化物作為特征因子[15-16]，食品加工企業(yè)則選擇鈉、鉀、氯化物作為特征因子[17]。

選擇不同外水的特征因子時，通常是根據(jù)某種水質(zhì)指標在不同外水中的濃度具有顯著性差異來判別的，是否具有顯著性差異往往需靠經(jīng)驗來判斷。本文利用統(tǒng)計方法給出特征因子的定量判別準則，可作為城市排水管網(wǎng)外水診斷時選擇特征因子的依據(jù)。

1 化學質(zhì)量平衡模型

排水管網(wǎng)某檢查井或泵站等節(jié)點中來水滿足以下水質(zhì)水量平衡方程組

式中：qi為第i 種來水的流量，i=1，2…n，n 為來水種類的數(shù)量；cij為第i 種來水的第j 種水質(zhì)指標濃度，j=1，2…m，m 為水質(zhì)指標的數(shù)量；Q 為該節(jié)點的總來水量；Cj為該節(jié)點中第j 種水質(zhì)指標的混合濃度。將式（1）等號兩邊同時除以Q，得到以下化學質(zhì)量平衡方程組

式中：xi為第i 種來水的流量占比。式（2）是一個（m，n）階非齊次線性方程組，附加一個流量平衡的約束條件，為了使方程組有唯一解，通常選擇的化學指標個數(shù)等于來水種類的個數(shù)，即m=n。當cij為隨機數(shù)時，式（2）不能通過矩陣分解、最小二乘法等線性方程組的解法進行求解，需借助蒙特卡羅等隨機方法求解[2]，求解基本步驟如圖1 所示。從節(jié)點檢測數(shù)據(jù)初步判斷上游來水的種類，從已建立的特征因子庫中挑選并隨機生成特征因子，代入式（2）中求解，循環(huán)此過程1 萬次，即可得到每種來水占比的后驗分布。

圖1 化學質(zhì)量平衡方程的蒙特卡羅計算流程

2 特征因子的特性

為了保證式（2）的方程組有唯一解，選擇的水質(zhì)指標需滿足以下條件[18]，方可作為特征因子：①同一種水質(zhì)指標在不同外水中濃度應(yīng)具有顯著性差異；②同一種外水中不同的水質(zhì)指標應(yīng)相互獨立，即線性無關(guān)；③水質(zhì)指標濃度的變化范圍應(yīng)足夠小，即較小的方差；④水質(zhì)指標應(yīng)具有穩(wěn)定的性質(zhì)，在管網(wǎng)中流動時不發(fā)生物理、化學變化。第④點要求是為了防止水質(zhì)指標濃度從上游到下游流動過程中發(fā)生衰減。第①至③點要求本質(zhì)上需要滿足方程組有唯一解。

將式（2）寫成矩陣形式

式中：A 為系數(shù)矩陣；b 為常數(shù)項向量。要保證方程（3）有唯一解，矩陣的秩需滿足以下條件

當要求①不滿足時，將矩陣A 按式（5）進行行變換，變換后矩陣A 的秩小于n，方程組有無窮解，因此特征因子必須滿足要求①。

當要求②不滿足時，將矩陣A 按式（6）進行行變換，變換后矩陣A 的秩小于n，方程組有無窮解，因此特征因子必須滿足要求②。

要求③與蒙特卡羅抽樣誤差有關(guān)，根據(jù)中心極限定理，誤差需滿足下式

式中：x 為水質(zhì)指標濃度樣本均值；μ 為水質(zhì)指標濃度均值；σ 為水質(zhì)指標濃度樣本標準差；Zα為標準正態(tài)分布分位數(shù)；α 為置信水平。當樣本數(shù)一定時，誤差與σ 成正比，因此當水質(zhì)指標濃度方差越小，計算誤差越小。

上述顯著性差異、相互獨立和較小方差的要求都是定性的，無法為特征因子的選擇提供依據(jù)，因此需要對這些要求提供具體的定量判別標準。

3 特征因子的選擇方法

特征因子顯著性差異要求可用Mann-Whitney U檢驗[19]；相互獨立要求可用Pearson 相關(guān)性檢驗[20]；較小方差要求可通過式（7）來確定標準差的置信范圍。

3.1 顯著性檢驗

Mann-Whitney U 檢驗，用來檢驗兩組數(shù)據(jù)是否服從同一分布，其計算步驟是將A、B 兩組數(shù)據(jù)混合后按升序排列，對于每個數(shù)據(jù)確定其在所屬組的序號，作為該數(shù)據(jù)的秩數(shù)，將兩組數(shù)據(jù)各自的秩數(shù)求和分別得到秩和RA和RB，然后按照式（8）計算U 值

式中：nA和nB分別為A、B 兩組數(shù)據(jù)的樣本數(shù)。通過UA和UB的最小值與預(yù)設(shè)值α 對比來進行假設(shè)推斷，如果p 小于α，則拒絕原假設(shè)H0：A、B 服從統(tǒng)一分布。

3.2 相關(guān)性檢驗

同一種外水中不同水質(zhì)指標的相關(guān)性用以下公式計算Pearson 相關(guān)系數(shù)

式中：X、Y 為同一種外水中兩種水質(zhì)指標濃度；μX和μY為X、Y 的均值；σX和σY為X、Y 的標準差。但通過相關(guān)系數(shù)來判斷的方式在實際應(yīng)用時仍然是定性的，因此需引入下式的相關(guān)性顯著性指標

式中：α 為相關(guān)性顯著性標準，通常設(shè)置為0.05；n 為X、Y 的樣本數(shù)量。當相關(guān)性顯著性指標p 小于0.05時，則拒絕原假設(shè)H0：X、Y 沒有線性相關(guān)性。

3.3 標準差檢驗

式（7）中的誤差可變形為下列形式

式中：ε 為絕對誤差，其他同式（7）。將式（12）等號兩邊同時除以均值μ 得到水質(zhì)指標的變異系數(shù)

式中：δ=σ/μ 為變異系數(shù)；εr為相對誤差。水質(zhì)指標的變異系數(shù)與樣本數(shù)量和相對誤差的關(guān)系如圖2 所示，圖2 是雙對數(shù)坐標。當樣本為30~60 個，相對誤差控制在5%~10%時，變異系數(shù)小于0.3 的才能作為特征因子。

圖2 變異系數(shù)與樣本數(shù)量和相對誤差的關(guān)系

4 結(jié)論

1）分析了城市排水管網(wǎng)外水診斷中特征因子需滿足顯著性差異、獨立性和方差較小要求的原因，目的是保證選擇的特征因子能使化學質(zhì)量平衡方程組有唯一解。

2）利用Mann-Whitney U 和Pearson 相關(guān)性檢驗等統(tǒng)計方法，給出了特征因子選擇流程，明確了特征因子顯著性差異和相關(guān)性判斷的定量標準。

3）闡述了水質(zhì)檢測數(shù)據(jù)的變異系數(shù)、樣本數(shù)量和相對誤差之間的關(guān)系，并建議在排水管網(wǎng)外水診斷實際工作中，當樣本為30~60 個，相對誤差控制在5%~10%時，選擇變異系數(shù)小于0.3 的水質(zhì)指標作為特征因子。