李青峰,郭 珊,2,張茹玉,閆曉滿,張學(xué)林
(1.廣州珠科院工程勘察設(shè)計(jì)有限公司,廣東 廣州 510610;2.珠江水利委員會(huì)珠江水利科學(xué)研究院,廣東 廣州 510611)
近年來,位于平原河網(wǎng)的城市經(jīng)濟(jì)發(fā)展迅速,快速的城市化進(jìn)程使得水污染問題日益嚴(yán)重[1-2],對(duì)平原河網(wǎng)水系的影響更為顯著[3-4]。目前,上海[5-6]、杭州[7]、武漢[8]、廣州[9]等城市均提出了一系列水利工程措施以改善水質(zhì),因此,研究平原河網(wǎng)地區(qū)的水質(zhì)改善措施具有重要意義。
水質(zhì)模型可應(yīng)用于量化評(píng)估一些具體工程措施的目標(biāo)可達(dá)性,目前用于水動(dòng)力水質(zhì)模擬較為廣泛[10-11],管儀慶等[12]利用MIKE11建立了臺(tái)州市的一維模型,用于計(jì)算區(qū)域水環(huán)境容量,為平原河網(wǎng)水環(huán)境管理提供依據(jù);崔廣柏等[13]的研究表明平原河網(wǎng)水動(dòng)力的提升對(duì)于河道水質(zhì)有明顯改善作用;高嵩等[14]運(yùn)用MIKE11水質(zhì)模型模擬分析了新建泵站對(duì)片區(qū)水質(zhì)改善的效果;田凱達(dá)等[15]基于一維模型建立了合肥十五里河水質(zhì)模型并應(yīng)用于不同水質(zhì)改善方案的效果評(píng)價(jià)。上述學(xué)者的研究表明MIKE11模型水動(dòng)力水質(zhì)模擬在平原河網(wǎng)應(yīng)用良好,因此可用于同樣為平原河網(wǎng)的珠三角河網(wǎng)區(qū)。
珠三角河網(wǎng)為典型的感潮河網(wǎng)區(qū),水環(huán)境問題隨經(jīng)濟(jì)快速發(fā)展日益突出,以往水質(zhì)改善措施分析大多聚焦于長三角平原河網(wǎng)區(qū)域,本文針對(duì)珠三角感潮河網(wǎng)復(fù)雜多變的特點(diǎn),以珠三角地區(qū)的佛山市南海區(qū)流域綜合治理試點(diǎn)——北村水系為研究區(qū)域,選取COD、NH3-N、TP為特征指標(biāo),利用MIKE 11模型分析不同工程措施實(shí)際應(yīng)用于珠三角河網(wǎng)的水動(dòng)力及水質(zhì)改善效果,為珠三角平原河網(wǎng)地區(qū)水生態(tài)環(huán)境的長效改善提供依據(jù)。
佛山市南海區(qū)北村水系流域地處廣州佛山交界,總面積246.79 km2。北村水系屬珠三角平原河網(wǎng),水系復(fù)雜,流域內(nèi)有雅瑤水道、大欖河、松崗河等32條總長度129.86 km的主干河涌和228條總長度184.94 km的支毛涌。工程區(qū)內(nèi)主要水系見圖1。各河涌水流主要匯集于雅瑤水道,經(jīng)北村水閘和北村泵站排入水口水道,另一部分排入佛山水道和西南涌。
圖1 北村水系流域范圍
北村水系河涌縱橫密布、閘站建設(shè)稍顯滯后,排澇閘站調(diào)度管理體系尚不完善,導(dǎo)致流域內(nèi)水體流動(dòng)性差、排澇能力不足;流域控污截污工作進(jìn)展慢,入河污染負(fù)荷大,水環(huán)境惡化、水生態(tài)退化等問題非常突出,大多數(shù)河涌檢測斷面的部分指標(biāo)不滿足地表水V類標(biāo)準(zhǔn)要求,超標(biāo)河涌的比例高達(dá)95%以上,水生態(tài)環(huán)境持續(xù)惡化問題亟待解決。
本文基于MIKE11中的水動(dòng)力模塊和對(duì)流擴(kuò)散模塊建立北村水系水動(dòng)力水質(zhì)模型,利用內(nèi)置的求解器進(jìn)行模型求解,模擬中小尺度河涌水體的水動(dòng)力水質(zhì)過程,為精準(zhǔn)解決珠三角平原河網(wǎng)區(qū)復(fù)雜的環(huán)境問題提供一些參考。
2.1.1水動(dòng)力控制方程
水動(dòng)力計(jì)算基于一維非恒定流Saint-Venant方程組,包括:
連續(xù)方程:
(1)
動(dòng)量方程:
(2)
式中Z——水位,m;t——時(shí)間,h;x、B——河道縱坐標(biāo)及水面寬度,m;Q——流量,m3/s;q——單寬流量,m2/s;A——過水面積,m2;K——流量模數(shù);g——重力加速度,m/s2。
求解方式為Abbott-Ionescu六點(diǎn)隱式有限差分法,隱式差分具有穩(wěn)定性好、計(jì)算精度高等特點(diǎn)。
2.1.2對(duì)流擴(kuò)散方程
MIKE11AD是根據(jù)水動(dòng)力條件模擬水體中對(duì)流擴(kuò)散過程的工具,基于對(duì)流擴(kuò)散方程計(jì)算。一維河流水質(zhì)模型的基本方程[16]為:
(3)
式中C——模擬物質(zhì)的濃度,mg/L;N——降解系數(shù),d-1;x——空間坐標(biāo);t——時(shí)間坐標(biāo);u——平均流速,m/s;Ex——對(duì)流擴(kuò)散系數(shù),m2/s。
Ex是一個(gè)綜合參數(shù)項(xiàng),對(duì)流擴(kuò)散模型通過經(jīng)驗(yàn)公式來估算對(duì)流擴(kuò)散系數(shù):
Ex=aVb
(4)
式中V——水動(dòng)力計(jì)算的流速,m/s;a、b——自行選取的參數(shù)。
2.2.1河網(wǎng)概化
北村水系河道縱橫交錯(cuò),河網(wǎng)極其復(fù)雜,在收集的水動(dòng)力及水文資料基礎(chǔ)上,根據(jù)河網(wǎng)概化的基本原則,進(jìn)行合理概化。本次模型范圍包含了流域內(nèi)雅瑤水道、大范河、松崗河、機(jī)場涌、豐崗公涌等32條主干河涌,以及部分相關(guān)聯(lián)的支涌,共模擬河涌41條,水閘13個(gè),泵站9個(gè),劃分?jǐn)嗝?78個(gè),斷面為實(shí)測河道斷面,間距100~500 m。對(duì)于流域范圍內(nèi)較小的支涌,將其支渠及排污口概化為點(diǎn)源輸入模型,計(jì)算河網(wǎng)平面布置及概化見圖2。
圖2 一維數(shù)學(xué)模型河網(wǎng)概化
2.2.2邊界條件
根據(jù)GB/T 25173—2010《水域納污能力計(jì)算規(guī)程》,設(shè)計(jì)水文條件應(yīng)采用90%保證率最枯月平均流量或近10年最枯月平均流量作為設(shè)計(jì)流量。本報(bào)告收集到大瀝氣象站2001—2016年降雨資料,經(jīng)分析,2014年1月份為近10年降雨量最小的月份,選取該月作為模型計(jì)算時(shí)段。
a)水位邊界。研究河段共有13條河涌通過閘站與外江相連,每條河涌水位邊界作為模型的開邊界水位邊界采用2014年1月實(shí)測的逐時(shí)潮水位,見圖3。
圖3 北村水閘外江實(shí)測水位值
b)水質(zhì)邊界。根據(jù)長江勘測規(guī)劃設(shè)計(jì)研究有限責(zé)任公司及中設(shè)設(shè)計(jì)集團(tuán)股份有限公司2018年5月編制的《南海區(qū)北村水系流域水環(huán)境綜合治理項(xiàng)目可行性研究報(bào)告》成果,外江水質(zhì)邊界條件見表1。
表1 研究區(qū)域各水體污染物指標(biāo)濃度情況 單位:mg/L
2.2.3初始條件
模型的水質(zhì)初始條件按照模擬期引水河道的水動(dòng)力水質(zhì)狀況確定,見表2。
表2 研究區(qū)域各水體污染物指標(biāo)濃度情況 單位:mg/L
2.2.4模型參數(shù)
a)水動(dòng)力模型參數(shù)。根據(jù)中國水利水電出版社2011年出版的《水工設(shè)計(jì)手冊(cè)》,河道糙率取0.025~0.030。計(jì)算時(shí)間步長為10 s。
a) COD濃度
b)水質(zhì)模型參數(shù)。綜合降解系數(shù)N是反映污染物沿程生物降解、沉降和其他物化等變化的綜合系數(shù)。根據(jù)2.2.5節(jié)中模型率定成果,本次采用的降解系數(shù):COD取0.1/d,NH3-N取0.07/d,TP取0.01/d。
2.2.5模型率定驗(yàn)證結(jié)果
根據(jù)北村水系水環(huán)境治理項(xiàng)目地表水檢測數(shù)據(jù),對(duì)7條河道中間位置進(jìn)行采樣,各個(gè)位置采取3個(gè)樣本,進(jìn)行采樣檢測得到COD、NH3-N、TP數(shù)值進(jìn)行平均。模型率定的主要參數(shù)為污染物降解系數(shù),結(jié)合已有研究成果選用降解系數(shù)見2.2.3節(jié)。通過模型計(jì)算與實(shí)測數(shù)據(jù)進(jìn)行對(duì)比,見表3??芍狢OD模擬除大坑涌相對(duì)誤差較大外,其余數(shù)據(jù)相對(duì)誤差在18%以內(nèi),NH3-N模擬相對(duì)誤差在17%以內(nèi),TP模擬相對(duì)誤差在18%以內(nèi),大部分水質(zhì)模擬值與實(shí)測值誤差能夠滿足模型計(jì)算要求。
表3 研究區(qū)域各河道水質(zhì)監(jiān)測實(shí)測數(shù)據(jù)與模型模擬對(duì)比
以北村水系松崗河片市控考核斷面水質(zhì)為考核目標(biāo),選擇COD、NH3-N、TP 3項(xiàng)指標(biāo)進(jìn)行模擬計(jì)算。分別針對(duì)現(xiàn)狀和2組工程方案進(jìn)行分析,見表4。針對(duì)方案1,截污管網(wǎng)建設(shè)將導(dǎo)致點(diǎn)源污染物濃度降低,在模型中減小邊界條件中點(diǎn)源污染物的濃度,針對(duì)疏浚清淤工程,可導(dǎo)致水流更加通暢,模型中通過修改河道斷面進(jìn)行模擬;針對(duì)方案2,是在方案1的基礎(chǔ)上增加調(diào)蓄池和新建泵閘引水調(diào)水,增加調(diào)蓄池使得面源污染物濃度降低,在模型中減小邊界條件中面源污染物的濃度,而新建泵閘工程使得引水水量加大,促進(jìn)水體交換速率。
表4 計(jì)算方案
為分析各河涌現(xiàn)狀及工程后污染物濃度,分別選取主、支河涌典型斷面的污染物濃度進(jìn)行分析,本文以北村水系中大欖河主河涌、鯉崗尾涌干河涌為例,各河涌采樣斷面布置見圖6。
3.2.1大欖河結(jié)果分析
大欖河起于黎崗水閘,止于松崗南海花卉城,涌面寬約40~60 m。圖7分別為大欖河典型斷面COD、NH3-N及TP濃度隨時(shí)間的變化情況,由圖7可知,通過閘站引水以及各指標(biāo)自身的衰減作用,大欖河各斷面位置處的污染物濃度呈震蕩式衰減。
a)COD濃度
大欖河3種方案下污染物濃度情況見表5,其中污染物濃度下降率為各方案污染物相對(duì)現(xiàn)狀污染物濃度下降的百分比。通過式(5)來估算污染物濃度下降率:
表5 大欖河各方案污染物濃度
(5)
式中S——污染物濃度下降率,%;a——現(xiàn)狀方案的污染物濃度,mg/L;b——各改善方案的污染物濃度,mg/L。
現(xiàn)狀模擬結(jié)果顯示,在現(xiàn)狀閘站調(diào)度引水條件下,大欖河水體交換速度較慢。各斷面COD、NH3-N以及TP濃度最終分別在45.0、3.0、0.8 mg/L左右。
方案1通過新建污水管網(wǎng)截污,降低了入河污染物濃度;同時(shí),對(duì)大欖河及沿河支毛涌進(jìn)行清淤疏浚來提高河道過流能力;綜合上述措施,大欖河可達(dá)到Ⅴ類水標(biāo)準(zhǔn),其中大欖河TP濃度下降率最大,為50%,說明截污及清淤疏浚措施效果良好。
方案2中,水系活水工程的實(shí)施,加快水體交換速度,進(jìn)一步稀釋了污染物濃度。各斷面COD、NH3-N以及TP濃度最終分別在30.0、1.7、0.4 mg/L左右,水質(zhì)較方案1稍有提升。
3.2.2湖馬灣污染物濃度變化分析
湖馬灣起于舊水廠,止于雅瑤水道,全長1.75 km,涌面寬約50 m。涌底淤積嚴(yán)重,汛期排水不暢。圖8分別為湖馬灣典型斷面COD、NH3-N以及TP濃度隨時(shí)間的變化情況,由計(jì)算結(jié)果可知,通過閘站向外江引水以及各指標(biāo)自身的衰減作用,湖馬灣各斷面位置處的污染物濃度呈震蕩式衰減。
湖馬灣3種方案污染物濃度改善情況見表6。與現(xiàn)狀方案對(duì)比,實(shí)施方案1與方案2大欖河內(nèi)污染物濃度均有所下降,且在方案2下污染物濃度進(jìn)一步降低。
表6 湖馬灣各方案污染物濃度
在現(xiàn)狀閘站調(diào)度引水條件下,湖馬灣各斷面COD、NH3-N以及TP濃度分別在73.0、5.2、1.8 mg/L左右。
方案1實(shí)施后,COD、NH3-N以及TP濃度較工程前都有所減小,各斷面COD、NH3-N以及TP濃度最終分別在32.0、2.0、0.7 mg/L左右。然而僅靠截污及清淤工程不能使湖馬灣水質(zhì)達(dá)到Ⅴ類水質(zhì)目標(biāo),但基本滿足消除黑臭水體的要求
方案2中,進(jìn)一步水系活水工程的實(shí)施使得雅瑤水道的水量、水質(zhì)都得到了提升,從而對(duì)湖馬灣水質(zhì)帶來正面影響,雅瑤河道的水回流到湖馬灣,進(jìn)一步稀釋了污染物濃度。COD、NH3-N以及TP濃度較工程前都有所減小,各斷面COD、NH3-N以及TP濃度最終分別在20.0、1.3、0.4 mg/L左右。水質(zhì)基本可以達(dá)到Ⅴ類需求。
a)佛山市南海區(qū)北村水系范圍屬于典型的珠三角平原河網(wǎng)區(qū),河網(wǎng)密度較大,但水流流速較慢,水體交換能力較差,水體污染后,通過自身的流通和凈化能力往往很難得到水質(zhì)改善。本文以北村水系為例,探討了珠三角河網(wǎng)區(qū)進(jìn)行水質(zhì)優(yōu)化的相關(guān)工程措施,定量分析各措施對(duì)水質(zhì)改善的效果,為類似工程提供一定的參考。
b)基于Mike11建立的北村水系一維河網(wǎng)水動(dòng)力水質(zhì)模型,分別對(duì)現(xiàn)狀以及各方案條件下相關(guān)水質(zhì)指標(biāo)的變化過程。分析得到通過現(xiàn)狀泵閘調(diào)控對(duì)水質(zhì)的改善效果甚微,而新建截污管網(wǎng)和疏浚工程即可大幅改善河道水質(zhì),對(duì)水質(zhì)改善的效果最為顯著。新建調(diào)蓄池和活水泵閘工程對(duì)于河段水質(zhì)改善影響程度不同,需綜合考慮其建設(shè)的必要性。
c)污染嚴(yán)重的支河涌,僅靠截污和疏浚無法滿足地表水體Ⅴ類水要求,進(jìn)一步實(shí)施活水工程后,水質(zhì)基本上可以達(dá)到Ⅴ類水要求,但水質(zhì)狀況極不穩(wěn)定,受外界擾動(dòng)影響極大??稍诜桨?的基礎(chǔ)上實(shí)施水系連通和生態(tài)修復(fù)工程,進(jìn)一步改善水環(huán)境質(zhì)量。