羅志潔，杜世鵬

（1.浙江同濟(jì)科技職業(yè)學(xué)院，浙江 杭州 310002；2.杭州市蕭山區(qū)浦陽(yáng)江流域管理中心，浙江 杭州 310002）

1 問(wèn)題的提出

杭嘉湖平原歷來(lái)是“魚米之鄉(xiāng)、富饒之地”，是中國(guó)長(zhǎng)三角經(jīng)濟(jì)區(qū)的重要組成部分。隨著越來(lái)越多污水處理廠的興建，其處理尾水的排污口設(shè)置對(duì)上下游、左右岸的影響研究和對(duì)區(qū)域納污能力的影響研究變得十分緊要[1-2]。

目前多采用先進(jìn)的水動(dòng)力數(shù)值模擬軟件，根據(jù)排污口項(xiàng)目的污水排放總量、濃度，建立水動(dòng)力學(xué)模型，模擬污染物在河道中的擴(kuò)散和降解過(guò)程，分析其對(duì)河道水質(zhì)的影響范圍和影響程度[3-4]。

本文以湖州市南潯區(qū)為研究區(qū)域，對(duì)區(qū)域內(nèi)污水處理廠進(jìn)行摸排，采用Mike11 構(gòu)建區(qū)域水動(dòng)力水質(zhì)模型，模擬計(jì)算不同工況下，新建舊館污水處理廠廢水排放對(duì)周邊水環(huán)境的影響，尤其對(duì)下游江蘇、浙江兩省交界斷面的水質(zhì)影響[5-8]。該研究方法對(duì)平原河網(wǎng)地區(qū)開展排污口影響分析工作有一定的借鑒意義。

2 區(qū)域水系概況

南潯區(qū)位于太湖流域水利分區(qū)中的“杭嘉湖區(qū)”，杭嘉湖區(qū)面積7 436 km2，其中浙江省范圍內(nèi)為6 481 km2（稱為杭嘉湖東部平原），其余屬于江蘇省和上海市。

南潯區(qū)屬杭嘉湖平原，北接太湖，東臨江蘇省，南接杭嘉湖平原腹地，域內(nèi)河網(wǎng)縱橫交錯(cuò)、四通八達(dá)，具有江南水鄉(xiāng)特色的網(wǎng)絡(luò)狀結(jié)構(gòu)。湖州市區(qū)水功能區(qū)劃見圖1。根據(jù)《浙江省水功能區(qū)水環(huán)境功能區(qū)劃分方案（2015）》，排污口所在河段水功能區(qū)名稱為頔塘湖州農(nóng)業(yè)、工業(yè)用水區(qū)。近年來(lái)《湖州市水資源公報(bào)》中城市內(nèi)河Ⅲ類水質(zhì)斷面達(dá)標(biāo)率均為100%，但是平原末端水體流動(dòng)性差，水利自凈能力差，加上排污量的提升，整體而言，平原區(qū)水體水環(huán)境容量小，水質(zhì)總體不容樂(lè)觀，見圖2～3。

圖1 湖州市區(qū)水功能區(qū)劃圖

圖2 2014、2017年頔塘南潯國(guó)控?cái)嗝鍯OD濃度圖

圖3 2014、2017年頔塘南潯國(guó)控?cái)嗝鍺H3-N濃度圖

3 區(qū)域排污口調(diào)查

根據(jù)浙江省《關(guān)于做好入河排污口設(shè)置審核登記工作的指導(dǎo)意見（2017）》，南潯區(qū)對(duì)轄區(qū)內(nèi)的企業(yè)入河排污口進(jìn)行了摸排及信息入庫(kù)，轄區(qū)內(nèi)有13 家企業(yè)設(shè)置入河排污口。本次新建的舊館鎮(zhèn)污水處理廠位于頔塘右岸，上游有中環(huán)水務(wù)公司排放口（2 812 m），下游有金潔水務(wù)公司排放口（3 464 m）（見圖4）。南潯區(qū)入河排污口信息見表1。

表1 南潯區(qū)入河排污口信息表

圖4 舊館鎮(zhèn)污水廠附近排污口位置示意圖

新建的舊館污水處理廠工程總規(guī)模為2.5×104m3/d，目前污水處理規(guī)模為1.5×104m3/d。

4 水動(dòng)力模型及水質(zhì)模型

4.1 水動(dòng)力模型

水動(dòng)力模型所用的描述一維非恒定水流運(yùn)動(dòng)規(guī)律的基本方程為圣維南方程組[9-10]，其數(shù)學(xué)表達(dá)式為：

式中：Q為流量，m3/s；A為斷面面積，m2；q為旁側(cè)入流，m3/s；h為水深，m；C為謝才阻力系數(shù)；R為水力半徑，m；g為重力加速度，m/s2；x為水流方向的距離變量，m；t為時(shí)間變量，s；α為動(dòng)量系數(shù)，一般取值為1。

4.2 水質(zhì)模型

描述物質(zhì)在水體中輸運(yùn)的一維非恒定流對(duì)流擴(kuò)散基本方程為：

式中：C為模擬水質(zhì)指標(biāo)濃度，mg/L；D為擴(kuò)散系數(shù)，m2/s；K為綜合衰減系數(shù)，d-1；C2為源匯項(xiàng)濃度，mg/L。

4.3 模型范圍及邊界條件

模型計(jì)算范圍包括南潯區(qū)域內(nèi)主要河道，模型計(jì)算區(qū)域基本覆蓋頔塘及其兩岸區(qū)域。水動(dòng)力上邊界為湖州船閘2014 年實(shí)測(cè)逐日流量過(guò)程，下邊界為菱湖站、烏鎮(zhèn)站、平望站、小梅口2014 年實(shí)測(cè)逐日水位過(guò)程；水質(zhì)上邊界采用三里橋站2014 年實(shí)測(cè)水質(zhì)數(shù)據(jù)，下邊界采用重兆大橋站、練市站、平望站實(shí)測(cè)水質(zhì)數(shù)據(jù)（見圖5）。驗(yàn)證過(guò)程結(jié)果見圖6～8。

圖5 水動(dòng)力水質(zhì)模型概化示意圖

圖6 2014年南潯站水位驗(yàn)證過(guò)程結(jié)果圖

經(jīng)參數(shù)率定，確定模擬河段內(nèi)糙率取值0.025，擴(kuò)散系數(shù)為12 m2/s，綜合衰減系數(shù)為0.072 d-1。由圖7～8 分析可知，模型計(jì)算初期，計(jì)算值小于實(shí)測(cè)值，模型穩(wěn)定后，計(jì)算值與實(shí)測(cè)值吻合較好，說(shuō)明模型參數(shù)設(shè)置基本合理，可以用該模型進(jìn)行預(yù)測(cè)分析。

圖7 2014年舊館站COD驗(yàn)證過(guò)程結(jié)果圖

圖8 2014年舊館站NH3-N驗(yàn)證過(guò)程圖

4.4 模型計(jì)算

在對(duì)區(qū)域水質(zhì)影響進(jìn)行分析之前，先分析頔塘歷年流量、水位、水質(zhì)過(guò)程，其中河道流量采用最不利的情況。

根據(jù)潯溪大橋站2001—2017 年長(zhǎng)系列流量資料，對(duì)歷年最小月平均流量進(jìn)行排頻，并采用P-III型頻率曲線進(jìn)行適線，得到均值為24.6 m3/s，Cv=0.23，Cs/Cv= 2.0（Cv為變差系數(shù)，Cs為偏差系數(shù)），對(duì)應(yīng)頻率90%的年最小月平均流量為17.6 m3/s，頻率曲線見圖9。

圖9 潯溪大橋站歷年月最小流量頻率曲線圖

根據(jù)南潯站、雙林站、烏鎮(zhèn)站、菱湖站（1991—2016 年）實(shí)測(cè)長(zhǎng)系列水位資料，對(duì)歷年最低水位進(jìn)行P-III 型頻率曲線適線，得出頻率90%的年最低水位分別為0.58 m、0.62 m、0.55 m、0.65 m，各站不同頻率統(tǒng)計(jì)成果見表2。

表2 水文要素表

4.5 計(jì)算工況

為研究排污口設(shè)置對(duì)頔塘河段的水質(zhì)影響，主要計(jì)算工況見表3。

表3 計(jì)算工況匯總表

4.6 計(jì)算成果

根據(jù)已建立的南潯舊館鎮(zhèn)污水處理廠局部河網(wǎng)水質(zhì)數(shù)學(xué)模型，計(jì)算得到不同工況下頔塘各斷面污染物濃度，分析排污口建設(shè)對(duì)河道水質(zhì)影響程度及范圍。不同工況下排污口入頔塘后COD 和NH3-N濃度的沿程變化情況見圖10～11，工程建設(shè)前后頔塘COD 和NH3-N 濃度的變化情況見表4～5。

圖10 頔塘COD濃度沿程變化

圖11 頔塘NH3-N濃度沿程變化圖

表4 工程建設(shè)前后頔塘COD濃度變化情況表

表5 工程建設(shè)前后頔塘NH3-N濃度變化情況表

續(xù)表5

分析可知，由于枯水期東部平原水位低于太湖水位，尤其在五大溇港工程實(shí)施后，五大溇港的水流流向總體由北向南，太湖水補(bǔ)充東部平原，頔塘水流的流向由西向東。中環(huán)水務(wù)的排污口和舊館鎮(zhèn)污水處理廠的排污口均位于羅溇與幻溇之間，因此幻溇的來(lái)水在頔塘交匯時(shí)，對(duì)上游頔塘來(lái)水形成頂托，導(dǎo)致該段流速、流量明顯減小。故中環(huán)水務(wù)公司排污口以下至幻溇與頔塘交匯處的水質(zhì)明顯較差，而經(jīng)過(guò)幻溇、濮溇、邢窯塘與頔塘的交匯處后，流速、流量變大，水質(zhì)明顯變好。

在舊館鎮(zhèn)污水處理廠項(xiàng)目完工且正常排放的情況下（按照一級(jí)A 標(biāo)準(zhǔn)排放），頔塘沿線COD濃度增幅較小，排污口上游影響至3 317 m 處，COD 濃度增加0.01 mg/L；增加最大的是排污口下游1 237 m 處，增加2.71 mg/L，增幅14.2%；排污口下游浙江省與江蘇省交界處COD 濃度沒有變化。頔塘沿線NH3-N 濃度增幅較小，排污口上游影響至2 812 m 處，NH3-N 濃度增加0.03 mg/L；增加最大的是排污口下游1 237 m 處，增加0.31 mg/L，增幅21.2%；排污口下游浙江省與江蘇省交界處NH3-N 濃度沒有變化。

本次事故為管道破損或污水處理系統(tǒng)未能正常運(yùn)行導(dǎo)致超標(biāo)排放。事故工況下，污水處理廠的廢水排放濃度相對(duì)正常工況較高。舊館鎮(zhèn)污水處理廠項(xiàng)目完成后事故排放與建設(shè)前正常排放相比，頔塘沿線COD 濃度增幅較大，排污口上游影響至3 317 m 處，COD 濃度增加0.01 mg/L；增加最大的是排污口下游390 m 處，增加23.28 mg/L，增幅104.0%；排污口下游浙江省與江蘇省交界處COD濃度增加0.03 mg/L。頔塘沿線NH3-N 濃度增幅亦較大，排污口上游影響至2 812 m 處，NH3-N 濃度增加0.05 mg/L；增加最大的是排污口下游390 m處，增加1.59 mg/L，增幅84.1%；排污口下游浙江省與江蘇省交界處NH3-N 濃度沒有變化。

因此，在現(xiàn)有舊館鎮(zhèn)污水處理廠未建時(shí)，中環(huán)水務(wù)排污導(dǎo)致羅溇與幻溇之間頔塘的水質(zhì)較差。舊館鎮(zhèn)污水處理廠新建后，正常排放工況下，水質(zhì)沒有明顯變差，對(duì)頔塘沿線水質(zhì)影響有限；但在事故排放工況下，由于排放水質(zhì)明顯低于排放標(biāo)準(zhǔn)，因此對(duì)羅溇與幻溇間的頔塘沿線水質(zhì)影響較大。

5 結(jié) 語(yǔ)

水污染防治是一項(xiàng)復(fù)雜的系統(tǒng)性工作，影響因素多，治理難度大。在加強(qiáng)入河排污口的監(jiān)督管理時(shí)必須充分論證排污口設(shè)置的可行性、合理性和對(duì)論證區(qū)域的水生態(tài)環(huán)境以及第三者的影響，須符合水域管理要求。

平原河網(wǎng)輸運(yùn)污水的速度較慢，工程排放污水會(huì)對(duì)周邊河網(wǎng)產(chǎn)生影響。污水處理廠正常排放時(shí)，對(duì)排污口下游一段區(qū)域的水質(zhì)有一定影響，遇到河道交匯，河流匯入后水質(zhì)緩解較快；但若污水處理廠發(fā)生事故，則對(duì)排污口上下游及其附近河道水質(zhì)影響較大，影響范圍較遠(yuǎn)。因此，在日常排污企業(yè)管理中，應(yīng)加強(qiáng)管理，提高安全生產(chǎn)意識(shí)，定期檢查、維護(hù)，確保處理系統(tǒng)安全、穩(wěn)定、可靠運(yùn)行，避免事故發(fā)生。