雨水管網是緩解城市內澇、改善水環境質量的重要市政基礎設施，受管材老化腐蝕和施工過程不當等影響，我國城市雨水管網建成后普遍存在地下水入滲和污水混接問題，制約著雨水管網功能高效發揮。排水管網入流入滲診斷技術可為雨水管網混接入滲問題的識別和定位提供支持。綜述了排水管網入流入滲診斷技術（物探檢測技術、流量分析技術、特征因子分析技術、水動力反演模型技術）的基本原理、研究進展、應用特征等，對比分析了不同技術應用于雨水管網混接入滲診斷的適用性、實施要點、診斷水平和發展階段等，并提出未來雨水管網混接入滲診斷技術的發展方向。
結果表明：
1）特征因子分析技術可在不干擾雨水管網正常運行條件下識別和量化混接入滲問題、提供治理優先級，因此經濟和環境效益顯著；
2）光纖分布式測溫（FDTS）技術和水動力模型反演技術在無干擾定位雨水管網混接入滲方面優勢顯著，前者易于實施且定位精準度高，后者具有成本低且可評級混接入滲等級特點，可根據場景特征因地制宜選用；
3）未來雨水管網混接入滲診斷技術將朝著低成本、無干擾、可量化、可定位的方向發展。通過評估不同技術在當前和未來的綜合應用效果，明確了單一診斷技術難以具備期望特征，而特征因子分析技術耦合FDTS技術或水動力模型反演技術的分級診斷體系具有實現未來技術發展目標的優勢潛力。
研究可為雨水管網混接入滲診斷技術的科學應用和優化創新提供參考。
1 前言
雨水管網是城市的生命線，對降低城市內澇風險和改善水環境質量意義重大。近10年來，我國城市化進程發展迅速，雨水管道規模以10.1%年均增長率快速增加，截至2021年底已達37.9萬km，其中約四成服役期超過10年。受施工過程不當、管材老化腐蝕和地面載荷等因素影響，我國城市雨水管網建成后普遍存在污水混接、高地下水位地區地下水入滲問題，其能夠對水環境健康與公共安全造成顯著的風險脅迫。
目前，我國城市的平均污水集中收集率僅為69%，意味著將近31%的污染物未經處理直接排放至河湖水體，從而引發嚴重的水環境污染問題。調查表明，污水通過混接至雨水管道直接排入河湖水體，是導致污水從收集系統流失和集中收集率偏低的主要原因之一。例如，在我國東部城市通過混接至雨水管道流失的污水量可達污水總量的26.0%~80.7%。然而，上述問題在我國分流制排水系統中普遍存在。2018年，上海市在1.9萬km雨污水管道中查出20290處混接點；2021年，武漢東湖高新區在2550km雨污水管道中發現600處混接點。混接污水往往具有高出地表水環境承載力的污染物濃度，嚴重者甚至與市政污水相當。因此，雨水管網中存在污水混接問題是造成城市受納水體返黑返臭、黑臭水體治理成效功虧一簣的主要原因之一。2021年11月，中共中央、國務院印發《關于深入打好污染防治攻堅戰的意見》，強調城市黑臭水體縱深治理要強化污水收集效能，而探明污水混接位置是污水收集處理提質增效的關鍵著力點。
在我國地下水位較高的南方地區，雨水管道普遍存在地下水入滲問題，通過文獻數據調研發現，地下水入滲量可達82.8~524.2 m3/（d·km）。入滲地下水能夠通過擠占雨水管道的收集輸送容量，損害其防洪排澇功能，導致地面積水無法通過雨水管道被及時排干，從而引發城市內澇問題；加劇雨水管道堵塞進而產生壅水問題，造成的管道水壓增大可促使窨井蓋發生移位，嚴重時釀成“井蓋吃人”的悲劇。此外，在地下水入滲雨水管道的過程中還可能將管道周邊的土體顆粒沖刷進入管道內，導致入滲區域上方形成侵蝕坑，這會加速管道失穩以致爆管，嚴重時甚至引發地面塌陷災害。由此可見，雨水管網存在地下水入滲問題已成為危害人民生命財產安全的重大隱患，及時開展整治工作意義重大，而查明入滲區域是整治工作高效開展的前提。
在上述背景下，筆者通過綜述國內外雨水管網混接入滲診斷技術的基本原理、應用現狀等，旨在為該類技術的科學應用和創新發展提供參考，有助于提升決策者對雨水管網混接入滲問題的治理水平，以及防范化解水環境污染和公共社會安全風險的戰略主動。
2 雨水管網混接入滲診斷技術類型、原理和應用現狀
根據原理不同，雨水管網混接入滲診斷技術主要分為物探檢測技術、流量分析技術、特征因子分析技術、水動力反演模型技術，下面分別進行論述。
2.1 雨水管網混接入滲診斷技術類型、原理和應用現狀
物探檢測技術始于20世紀60年代，其基本原理是采用感知設備采集管道內部的物理參數，直接或間接為排水管道入流入滲問題的診斷提供分析依據。目前，物探檢測技術基本形成了以管道閉路電視成像（CCTV）為主的傳統檢測技術，以紅外成像技術和光纖分布式測溫（FDTS）為代表的新型檢測技術。
2.1.1 CCTV技術
現階段，CCTV應用最為廣泛，其主要由主控器、操縱線纜架、搭載攝像鏡頭的行進機器人組成。該技術可通過實地拍攝影像輔助工作人員直觀判斷排水管道存在的入流入滲問題。但為保證機器人在管道中順利行進和鏡頭正常拍攝，CCTV技術應用前需要進行管道斷水、排空、清淤等復雜操作，不但造成檢測費用昂貴、效率低下，而且在面對流量大、水位高的管道時甚至不能實現斷水檢測。為克服上述難題，國內研發了基于漂浮滾筒螺旋式行進的CCTV帶水檢測技術，但該設備易受水流變化影響而較難被操控，導致數據質量無法保證；對此，蔡兆祝通過解決管道復雜環境下機器人運動適應性低、檢測方式單一、成像效果差的問題，研發了能在少水、滿水甚至有水流沖擊環境下進行檢測的兩棲機器人，然而該技術尚未實現工程應用。
2.1.2 紅外成像技術
紅外成像技術是采用紅外成像設備掃描物質表面的溫度場，并根據溫度突變診斷管道的入流入滲問題。例如，Lepot等采用紅外成像設備對實驗水槽的748次旁側入流事件進行了記錄，結果表明，該技術主要適用于探測導致管道自由水面溫度場發生變化的入流入滲事件，且檢出限非常低，能夠識別的最小入流僅為管道水體流量的2.5%，但無法定量入流流量以及識別侵入型、裂縫型、管道底部型入流事件。紅外成像技術能有效實現不斷水檢測，但在實際應用中需將成像設備固定在漂浮裝置上采集數據，導致設備易受水流流速變化影響而無法獲得精準數據。當前該技術尚無應用于實際管道檢測的案例報道。
2.1.3 FDTS技術
FDTS技術能基于沿管道鋪設的光纖高頻率感知空間溫度變化，并根據溫度變化特征診斷管道的入流入滲問題，其原理如圖1所示。近年來，FDTS技術已在國內外形成許多試驗場景。例如，Hoes等在2009年率先采用FDTS技術對荷蘭Korendijk州和Groningen州雨水管網中存在的生活污水間歇性接入事件進行了有效追蹤，并指出該技術主要適用于對大流量入流水體的探測，然而該研究未給出FDTS技術識別混接事件的溫度變化閾值；隨后，Kessili等在實驗室條件下分析了FDTS技術的溫度變化閾值對管道充滿度、混接污水溫度和體積、混接路徑長度的響應關系，發現管道充滿度、混接污水和管道環境之間的溫度差異、混接路徑長度是影響混接事件被成功識別的主要因素，并綜合確定FDTS技術識別污水混接事件的溫度變化閾值為0.15 ℃，然而該研究未分析溫度變化閾值的時變化特征；對此，尹海龍等在2022年以1條50 m長的污水管道為例開展了基于FDTS技術的污水和雨水入流實時監控，并基于統計學理論分析了FDTS技術識別入流事件時在空間和時間上的溫度變化閾值（分別為±0.2和±0.5 ℃）。然而，當入流入滲水體流量小、與管道原生水體之間溫差不明顯時，FDTS技術不適用。
圖1 基于FDTS技術的排水管道入流入滲診斷示意
2.1.4 技術適用條件
物探檢測技術是診斷雨水管網混接入滲的有效手段，其中，CCTV技術主要適用于水深較淺、淤積較輕等無需實施降水、清淤操作的情形，否則將面臨高昂的檢測成本；紅外成像技術和FDTS技術均支持帶水檢測，前者應用時要求混接入滲引起的溫度變化能夠在管道水體自由水面顯現以及管道水流流速平穩，后者主要適用于混接入滲水量大以及與管道原生水體之間溫差顯著的情形，同時還應盡可能保持光纖漂浮于水體自由表面并避免光纖彎曲引發定位誤差和激光信號損失等。然而，上述方法均難以量化混接入滲的水量信息，以至于無法通過評級問題嚴重程度，實現治理效果與經濟效益的統一。
2.2 流量分析技術
流量分析技術是實現排水管道入流入滲定量解析的重要方法，其基本原理：基于排水管網節點的流量監測數據和污染源的污水排放數據，構建封閉區域內排水系統的水量平衡方程，并通過方程求解量化排水管網中不同來源類型水體的水量輸入情況。例如，Almeida等采用37臺管道流量計對Costa do Estoril的污水管網進行了節點流量監測，并結合用水戶的用水量數據、水質指標分析、夜間最小流量法等，分析了污水管網的地下水入滲量、工業廢水排放量、雨水混接量，在此基礎上提出了污水管網的優化改造方向；Xu等綜合運用管道流量監測數據、污染源調查數據、泵站運行數據、地下水入滲模型，構建了上海市某分流制排水系統的水量平衡方程（圖2），通過方程聯立求解，定量解析了旱季期間雨水管網中的污水混接、地下水入滲和河水倒灌情況，從而為雨水管網混接入滲問題治理提供了科學指導。然而，流量分析技術需要通過獲取大量管道流量數據明確水量平衡方程的邊界條件，具有成本高、人員投入大的問題。對此，肖濤提出了基于管道液位演算流量的“軟測量”模型——單井模型和雙井模型。單井模型是以上游檢查井液位監測數據為模型輸入條件演算流量，適用于下游管道順坡、水力狀態良好的檢查井；雙井模型是以上下游檢查井的液位監測數據為模型輸入條件演算管道流量，在順坡、平坡和逆坡管道中均適用。然而，在實際應用中，測量誤差或管道淤積常導致演算流量與實測數據之間存在較大誤差。
圖2 基于水量平衡分析的雨水管網混接入滲量化模型
除上述不足外，流量分析技術還存在無法識別水體來源類型的缺陷，導致在開展排水管網入流入滲定量解析時，易將未監測到的水體錯誤納入已知來源類型水體中，引發解析結果的較大偏差。對此，工程實踐中常引入水質指標進行輔助分析，以提升診斷結果的準確性。例如，在《城鎮排水管道混接調查及治理技術規程》中給出了通過監測雨水管道水體污染物濃度和負荷量的沿程變化診斷污水混接的技術路線。然而，該類方法仍無法精細化定量不同來源類型水體的水量。
綜上可知，流量分析技術是定量解析雨水管網混接入滲的有效手段，主要適用于水量平衡模型邊界條件較少的情形，否則需要安裝大量管道流量計，導致診斷成本顯著增加；此外，該技術對監測節點水體的流速和水深大小具有特定要求，在測量過程中還需避免傳感器被垃圾纏繞或淤泥淤積引發數據測量誤差。與物探檢測技術相比，流量分析技術無需深入管道內部獲取數據，因此操作易于實現，且能夠通過評級問題嚴重等級明確混接入滲區域治理的優先級，提高治理工作的經濟和環境效益。然而該技術不能識別水體來源類型，導致解析結果的不確定性較大，易誤導工程治理方向；另外，其診斷精度主要為區域水平，為實現定位還需在各檢查井開展流量監測，成本較高。
2.3 特征因子分析技術
2.3.1 基本原理
特征因子分析技術是精細化識別和定量排水管網入流入滲的有效手段。該技術是美國國家環境保護局（US EPA）開展雨水管網混接調查的核心方法之一，由于該技術允許帶水操作且能夠定量解析不同來源水體的貢獻比例，近年來不斷受到重視。特征因子分析技術的基本思想：篩選針對不同來源類型水體的特異性和保守性特征因子，在此基礎上建立涉及特征因子指標和貢獻比例的化學質量平衡（CMB）模型，并通過CMB模型求解，定量解析排水管網中不同來源類型水體的接入比例，在此基礎上進一步結合流量監測數據可量化評估入流入滲水平。該方法架構如圖3所示。
圖3 基于特征因子的排水管網入流入滲診斷技術架構
（1）CMB模型原理
CMB模型的數學表達如式（1）~式（2）所示。
（1）
式中：
C1,1,…,C1,k+1為第1個特征因子在排水管網上游監測節點入流水體中、旁側第1,…,k種來源類型水體中的指標值；Ck+1,1,…,Ck+1,k+1為第k+1個特征因子在排水管網上游監測節點入流水體中、旁側第1,…,k種來源類型水體中的指標值；f上游、f1…fk為排水管網上游監測節點入流水體、旁側第1,…,k種來源類型水體對下游監測節點出流水體的貢獻比例，滿足f上游+∑ｋｉ＝１fi=1；C1，…，Ck+1為排水管網下游監測節點出流水體中第1,…,k+1種特征因子的指標值。
不同來源類型水體的入流入滲水量大小為：
Qi=Q下游fi （2）
式中：Qi為排水管網旁側第i種來源類型水體的日均入流入滲量，m3/d；Q下游為排水管網下游監測節點出流水體的日均流量值，m3/d；fi為排水管網旁側第i種來源類型水體對下游監測節點出流水體的貢獻比例，i取1，…，k。
（2）特征因子篩選準則和常用類型
特征因子分析技術的實施關鍵之一是選擇合適的特征因子，理論上特征因子應具有以下基本特征：1）在不同來源類型水體之間具有顯著的指標值差異；2）在排水管道的生物地球化學作用過程中能保持較好的穩定性；3）檢測限、測試精度、安全性和重現性較理想。參考上述原則，國內外學者篩選出能夠表征生活污水、不同行業廢水、地下水的濃度型特征因子（表1），并廣泛用于排水管網入流入滲診斷。例如，徐祖信等以上海市某雨水管網為研究對象，分別采用總氮、氟化物、硬度作為表征生活污水、半導體工業廢水和地下水的特征因子，結合CMB模型耦合不確定性算法量化了3種水體的貢獻比例；王詩婧以巢湖市城區污水管網為研究對象，采用安賽蜜作為區分生活污水和地下水的特征因子，并通過CMB模型求解明確了局部管段的嚴重入滲是導致污水廠進水濃度偏低的主要原因，并提出通過識別和修復局部嚴重破損管段提升污水處理效能的技術思路。
表1 不同來源類型水體的特征因子及其濃度分布范圍
2.3.2 傳統濃度型特征因子的時空非保守性問題
研究表明，同一來源類型水體的濃度型特征因子指標值常表現出顯著的時空波動性，是造成CMB模型解析結果出現較大誤差的主要原因之一。表2統計分析了在同一點位連續監測的濃度型特征因子指標值的變異系數（CV），結果表明，約46%監測數據的CV大于閾值15%，反映了數據離散程度較大，特征因子保守性偏低。同一來源類型水體的濃度型特征因子指標值在空間上表現出的顯著性差異主要與人類活動、區域環境特征等有關，例如生活污水的濃度型特征因子指標值大小與居民生活水平、居住區建成年代、排放污水中黑水和灰水混合比例等因素相關，地下水的濃度型特征因子指標值易受不同區域含水層物理化學性質的影響等。
表2 在同一點位連續監測的濃度型特征因子指標的CV值
2.3.3 新型指紋型特征因子的優勢前景
為克服濃度型特征因子指標值的時空非保守性，國內外學者將穩定同位素、三維熒光光譜、傅里葉變換離子回旋共振質譜等指紋技術引入該領域中。例如，Houhou等采用水體中氫原子和氧原子的穩定同位素（分別表示為δ2H和δ18O）對法國大南錫地區的自來水、市政污水、河水進行了識別，根據δ2H-δ18O指紋圖譜可知，不同來源類型水體之間的δ2H-δ18O指紋圖譜的差異性顯著，同一來源類型水體的δ2H-δ18O指紋圖譜的穩定較高，這為基于δ2H-δ18O指紋圖譜揭示污水管網水量平衡過程創造了有利條件。然而，該研究未從統計學角度給出δ2H和δ18O保守性的解釋；Kracht等在利用δ2H-δ18O指紋圖譜區分生活污水和地下水的研究中發現，δ2H-δ18O的CV分別小于-2.0%和-1.8%，表明其保守性較高，可作為精準量化污水管網中地下水入滲的特征因子。在熒光指紋應用方面，Chen等率先采用三維熒光光譜耦合平行因子分析法表征了雨水管網雨天出流水體、地表徑流和生活污水中的溶解性有機物熒光組分，并通過蛋白質類熒光峰識別了雨水管道中的生活污水接入。然而，該研究未對熒光指紋信息進行數值化表征，因此無法明確其保守性。對此，Liao等構建了三維熒光光譜指紋的衍生指標——降解指數，并采用顯著性檢驗證明了其對地表徑流（0.047）和生活污水（0.46）的特異性指示功能，然而該指標存在保守性不足問題（CV≤32%）。除上述方法外，傅里葉變換離子回旋共振質譜也是識別排水系統入流入滲的潛在指紋手段，由于該技術能從分子式、元素組成和群組成分的角度刻畫水體中溶解性有機物的來源信息，近年來已有效用于區分污水廠出水和天然水體、地表徑流和生活污水。通過將分子指紋信息數值化，該技術可進一步用于定量解析混合體系中不同來源類型水體的組成比例，但相關研究鮮有報道。
2.3.4 技術適用條件
綜上可知，特征因子分析技術能夠實現對雨水管網混接入滲水體的精細化定量解析，但主要適用于不同來源類型水體的特征因子之間存在顯著性差異的場景。由于該技術僅需開展特征因子檢測和部分點位流量監測，因此其診斷成本較低。然而，該技術的診斷精度一般為區域水平，為實現對混接入滲問題的定位還需在各檢查井開展水質監測，成本較高。
2.4 水動力反演模型技術
為提高對排水管網入流入滲問題的定位效率同時降低成本，水動力反演模型技術被提出。根據水動力學理論，水體在管道上游的時空輸入特性和在管道內的匯流傳輸過程決定了未來在下游點位的動力學形態（流量、水位等）。基于該關聯機制，國內率先發展了融合排水管網水動力模型和反問題理論的水動力反演模型技術，用于追蹤入流入滲問題。例如，2020年Zhao等采用圖4中的方法架構開展了巢湖市城區污水管網的地下水入滲定位和定量解析。具體地，首先采用特征因子分析技術識別出存在嚴重地下水入滲問題的污水管網子系統，同時給出該區域的水動力學邊界；其次通過耦合排水管網水動力模型（SWMM）和微生物遺傳算法構建子系統的地下水入滲反演模型，并用于推演地下水入滲的空間分布和水量大小，從而為嚴重入滲管道的優先修復提供靶點。隨后，Xu等在2021年綜合采用特征因子分析技術和“軟測量”方法識別了XiaoHecha雨水管網中存在嚴重混接入滲問題的子系統，并明確了其水動力邊界條件；在此基礎上，進一步通過耦合SWMM和粒子群算法構建了子系統的水動力反演模型，進而為存在嚴重混接入滲問題管道的優先修復提供了診斷技術支持。然而，上述研究均未開展反演模型定位效果對水動力（如曼寧粗糙系數、管道坡度）和算法（迭代次數、種群大小等）參數變化的響應機制分析，導致解析結果的不確定性較大；此外，當排水管網復雜度增加時，水動力模型運算效率以及尋優算法對解空間的全局遍歷性將隨之下降，成為制約該方法推廣應用的重要瓶頸。
圖4 基于特征因子和水動力反演模型的污水管網地下水入滲定位方法架構
綜上可知，水動力模型反演技術是開展雨水管網混接入滲問題診斷的新型手段，主要適用于水動力邊界明確和物理參數完整、準確的雨水管網，其實施關鍵是精細化、精準化構建雨水管網水動力過程模擬系統。與水量分析技術和水質特征因子分析技術相比，該方法無需逐檢查井開展水量或水質監測即可實現對混接入滲問題的定位，因此具有投資少、工作效率高的優點。此外該技術還支持對混接入滲水量定量的功能，因此能夠突破物探檢測技術難以評級問題等級的技術瓶頸。然而，如何解決反演模型的“異參同效”效應是當前該技術主要面臨的難題。
3 雨水管網混接入滲診斷技術對比及發展趨勢
綜合2節所述，總結了不同技術應用于雨水管網混接入滲診斷的適用場景、實施要點、診斷水平和發展階段等（表3），在此基礎上提出未來雨水管網混接入滲診斷技術將朝著低成本、無干擾、可量化、可定位的方向發展。根據表3可知，單一診斷技術均不具備未來預期特征，而通過技術組合并遵循分級診斷路線有望實現目的。
表3 不同技術應用于雨水管網混接入滲診斷的特征、水平和發展階段對比
根據表3可知，流量分析技術和特征因子分析技術是針對大范圍雨水管網分區開展混接入滲診斷的有效手段，可通過快速量化不同分區雨水管網的混接入滲情況、評級嚴重問題區域，為后續采用CCTV技術、紅外成像技術、FDTS技術、水動力反演模型技術開展問題定位提供靶區。然而，構建何種分級診斷技術組合路線是當前面臨的重要難題。綜合考慮技術成本、診斷精度、診斷水平、發展階段4個評價因子，對不同分級診斷技術組合路線在當前和未來的應用特征進行綜合分析，結果如圖5所示。根據圖5可知，針對大范圍雨水管網分區開展混接入滲診斷時，特征因子分析技術整體優于流量分析技術，且前者未來仍具有優化提升空間。通過對比兩組分級診斷技術組合路線可知，不論是當前還是未來預期，組合路線①~④的綜合性能整體優于組合路線1~4，表明采用特征因子分子技術耦合混接入滲定位技術更符合未來發展趨勢。分級診斷技術組合路線①~④在當前和未來的綜合評分分別遵循以下順序：④>③=①>②和④>③=②>①。上述結果表明，分級診斷技術組合路線③和④的綜合性能較高，且符合成本低、無干擾、可量化、可定位的預期優勢特征，因此二者有望成為未來雨水管網混接入滲診斷的核心方法。雖然技術組合路線②與③未來具有相同的綜合性能，但前者達到預期效果所需優化提升的程度較大。此外，現階段當考慮診斷結果的準確性時，分級診斷技術組合路線①將是首選。
圖5 不同分級診斷技術路線的綜合評價
4 結語
雨水管網存在污水混接和地下水入滲問題可對水環境健康和公共安全產生風險脅迫，因此及時掌握雨水管網混接入滲的準確信息將有助于提升防范化解水環境污染和公共社會安全風險的戰略主動。通過綜述國內外雨水管網混接入滲問題診斷技術的基本原理、研究進展、技術特征等，旨在為該類技術的科學應用和創新發展提供借鑒，并得出以下主要結論：
1）特征因子分析技術是開展雨水管網混接入滲診斷的有效手段，其可在不影響雨水管道正常運行情況下分區識別和量化混接入滲問題。通過進一步評估問題等級，該技術可為問題嚴重區域管道的優先修復提供靶區，因此具有實現治理效果與經濟效益統一的優勢特征。
2）FDTS技術和水動力模型反演技術在無干擾定位雨水管網混接入滲問題方面均展現出了優勢前景，其中前者易于實施且診斷結果的不確定性小，后者具有診斷成本低、可量化混接入滲水量的優勢，因此決策者可根據應用場景特征因地制宜選擇診斷方法。
3）研究認為未來雨水管網混接入滲診斷技術將朝著低成本、無干擾、可量化、可定位的方向發展，而單一診斷技術難以滿足要求，采用特征因子分析技術耦合FDTS技術或水動力模型反演技術的分級診斷體系有望實現目的。
