彭欣雨，唐 彥，唐德善，楊慶慶

(1.河海大學(xué)水利水電學(xué)院，南京 210098；2.河海大學(xué)計(jì)算機(jī)與信息學(xué)院，南京 210098； 3.揚(yáng)州市勘測(cè)設(shè)計(jì)研究院有限公司，江蘇 揚(yáng)州 225000)

榆林市工業(yè)的發(fā)展、河川徑流量逐步減少、廢污水堵截、儲(chǔ)蓄利用量不斷增加從而導(dǎo)致流域河道內(nèi)污染物質(zhì)不斷積累。流域內(nèi)超標(biāo)超量排放污水污染物的情況時(shí)有發(fā)生，主要污染物為化學(xué)需氧量、氨氮等。為改善榆溪河流域水質(zhì)污染嚴(yán)重的現(xiàn)狀，急需科學(xué)的水質(zhì)管理方法來(lái)指導(dǎo)水污染防治規(guī)劃研究。

近些年來(lái)，許多專家學(xué)者從不同的角度、運(yùn)用不同的方法對(duì)水污染防治進(jìn)行了研究。夏菁[1]等，采用TMDL計(jì)劃，針對(duì)長(zhǎng)湖進(jìn)行了水污染總量控制；牛麗東[2]等采用WARMF水環(huán)境管理模型，制定了污染物削減方法和分配量；陶亞、夏建新[3,4]等采用TMDL計(jì)劃，計(jì)算分析了深圳灣點(diǎn)源與非點(diǎn)源污染負(fù)荷、各子流域污染貢獻(xiàn)率及BOD5最大日污染負(fù)荷；王生愿[5]采用負(fù)荷歷時(shí)曲線法(LDC)，研究梁子湖流域水體納污能力，確定最大日負(fù)荷總量。以榆溪河流域?yàn)檠芯繉?duì)象，采用數(shù)學(xué)模型模擬的方法，計(jì)算水環(huán)境容量值(TMDL)，分析研究各個(gè)水源地對(duì)榆溪河流域的污染入河量、現(xiàn)狀點(diǎn)源和非點(diǎn)源污染負(fù)荷量，從而對(duì)污染負(fù)荷進(jìn)行分配并提出相關(guān)的污染削減措施。

1 材料與方法

1.1 研究區(qū)域概況及污染源調(diào)查

榆溪河是黃河一級(jí)支流無(wú)定河中游的主要支流，發(fā)源于榆陽(yáng)區(qū)小壕兔鄉(xiāng)，在榆陽(yáng)區(qū)魚(yú)河鎮(zhèn)匯入無(wú)定河，河長(zhǎng)97.8km。榆溪河流域西北接內(nèi)蒙古自治區(qū)，西南與橫山縣相連，東北為禿尾河流域，東南為佳蘆河流域，流域總面積4 938 km2，其中榆陽(yáng)區(qū)境內(nèi)匯流面積3 796 km2，內(nèi)蒙古境內(nèi)匯流面積1 142 km2。

將榆溪河流域分為小壕兔、頭道河則、金雞灘、色草灣、趙家灣、芹河、馬合、沙河、水掌灣及牛家梁-紅石峽10個(gè)水源地，進(jìn)行榆溪河流域污染源研究。10個(gè)水源地合計(jì)面積3 637.2 km2，占榆溪河流域(榆林市境內(nèi))面積3 796 km2的95.8%，見(jiàn)圖1。

圖1 榆溪河流域劃分Fig.1 Division map of Yuxi River Basin

1.2 TMDL計(jì)劃

流域 TMDL 計(jì)劃的關(guān)鍵技術(shù)主要包括流域地表水水質(zhì)受損水體和潛在污染源識(shí)別，選擇與構(gòu)建污染負(fù)荷和水質(zhì)關(guān)聯(lián)技術(shù)模型，計(jì)算最大容許污染負(fù)荷總量，確定點(diǎn)源污染負(fù)荷和非點(diǎn)源污染負(fù)荷分配方案，構(gòu)建公平、合理、可行的點(diǎn)源污染負(fù)荷削減分配模型，優(yōu)化非點(diǎn)源污染控制措施[6]。明確了制定流域TMDL計(jì)劃的步驟[7]: ①分析流域特征；②計(jì)算污染負(fù)荷容量；③確定污染負(fù)荷分配；④提交TMDL報(bào)告。

TMDL計(jì)劃的主要目標(biāo)就是將污染負(fù)荷在各個(gè)污染源之間合理分配，包括點(diǎn)源和非點(diǎn)源。進(jìn)行污染負(fù)荷分配依據(jù)的方法為[8]：

TMDL=∑WLA+∑LA+MOS

(1)

式中：TMDL為流域內(nèi)某特定水質(zhì)指標(biāo)的最大日污染負(fù)荷；WLA為允許的現(xiàn)存和未來(lái)點(diǎn)源的污染負(fù)荷；LA為允許的現(xiàn)存和未來(lái)非點(diǎn)源的污染負(fù)荷；MOS為預(yù)留安全值。

2 流域TMDL計(jì)劃模型構(gòu)建

2.1 選擇水質(zhì)指標(biāo)

水質(zhì)評(píng)價(jià)是將水質(zhì)指標(biāo)的監(jiān)測(cè)結(jié)果與水質(zhì)標(biāo)準(zhǔn)進(jìn)行比較，判定水質(zhì)狀況，水質(zhì)評(píng)價(jià)是工程設(shè)計(jì)水資源評(píng)價(jià)中非常重要的內(nèi)容[9]。利用不同水環(huán)境評(píng)價(jià)模型進(jìn)行污染程度評(píng)價(jià)，如單因子評(píng)價(jià)法、污染指數(shù)法、模糊指數(shù)法等[10,11]。其中單因子評(píng)價(jià)法具有能直接反映出水質(zhì)狀況與評(píng)價(jià)標(biāo)準(zhǔn)之間關(guān)系的優(yōu)點(diǎn)，能快速的反映出各河道污染指數(shù)的變化。

單因子評(píng)價(jià)指標(biāo)法。單項(xiàng)水質(zhì)參數(shù)(污染物指標(biāo))i在e點(diǎn)的評(píng)價(jià)指數(shù)Sie的計(jì)算公式為：

Sie=Cie/Ci

(2)

式中：Cie為污染物指標(biāo)i在監(jiān)測(cè)點(diǎn)e的數(shù)值，mg/L；Ci為污染物指標(biāo)i的環(huán)境標(biāo)準(zhǔn)值，mg/L。

研究表明，有機(jī)物污染指標(biāo)COD、氨氮的評(píng)價(jià)指標(biāo)最大，是造成榆溪河水污染的主要元素，所以本次TMDL的計(jì)算選取的指標(biāo)為污染物指標(biāo)COD、氨氮。

2.2 TMDL值計(jì)算

水環(huán)境容量由稀釋容量和自凈容量構(gòu)成，其中自凈容量一般較小，故忽略其影響并以最大稀釋容量代替水環(huán)境容量[12]。水環(huán)境容量計(jì)算公式為：

Mi=Q(Ci/C0i)/100

(3)

式中：i為污染物指標(biāo)(本文為COD、氨氮)；Mi為污染物指標(biāo)i的環(huán)境容量，萬(wàn)t/a；Q為多年平均地表水資源量，億m3/a；Ci為污染物指標(biāo)i的環(huán)境標(biāo)準(zhǔn)值，mg/L，采用《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838-2002)中Ⅱ類水標(biāo)準(zhǔn)值；C0i為污染物指標(biāo)i的環(huán)境背景值，mg/L，這里忽略本底值的影響，設(shè)C0i=0。

水資源量采用榆林市韓家峁水文站和高家堡水文站為參證站，榆林水文站地處韓家峁、高家堡2站之間，采用30 aP=90%保證率推求榆林水文站多年平均徑流量，相應(yīng)得到榆溪河流域COD、氨氮的水環(huán)境容量。

表1 榆溪河COD、氨氮的水環(huán)境容量 萬(wàn)t/a

2.3 安全臨界值確定

目前國(guó)內(nèi)大多數(shù)依據(jù)經(jīng)驗(yàn)定性確定安全臨界值的范圍，一般取水體最大納污量的5%～10%，取7%。

2.4 污染負(fù)荷分配

計(jì)算污染負(fù)荷分配權(quán)重時(shí)，采用客觀賦權(quán)法中的離差和均方差法求得點(diǎn)源和非點(diǎn)源污染的權(quán)重進(jìn)行負(fù)荷分配，見(jiàn)表2。該方法的計(jì)算步驟如下。

(1)求隨機(jī)變量的均值：

(2)求第j種污染源的均方差：

(3)求第j種污染源的權(quán)重：

(4)進(jìn)行多指標(biāo)決策與排序：

表2 榆溪河流域允許負(fù)荷分配 t/a

2.5 現(xiàn)狀污染負(fù)荷估算

污染物入河量計(jì)算：

Rijk=Pijkβj

(4)

式中：i為污染物指標(biāo)(i=1,2,…,n)；j為流域中污染源的種類(j=1,2,…,m)；k為流域切分的水源地(k=1,2,…,h)；Rijk為第k個(gè)水源地第j種污染源的第i種污染物的入河量，t/a；Pijk為第k個(gè)水源地內(nèi)第j種污染源的第i種污染物的總量，t/a；βj為第j種污染源的入河系數(shù)，農(nóng)村生活污染物入河系數(shù)取值0.3，農(nóng)田污染物的入河系數(shù)取值0.2，工業(yè)污染物入河系數(shù)取值0.8，畜禽養(yǎng)殖入河系數(shù)取值0.2，水產(chǎn)養(yǎng)殖入河系數(shù)取值0.3。

計(jì)算污染物入河量后，相應(yīng)計(jì)算污染物削減量和污染物削減率。

2.6 結(jié)果與討論

(1)現(xiàn)狀點(diǎn)源負(fù)荷估算。根據(jù)《榆林市統(tǒng)計(jì)年鑒》及2014年榆溪河流域的污染物濃度調(diào)查統(tǒng)計(jì)分析，點(diǎn)狀污染源主要包括生活污染源、工業(yè)污染源、規(guī)?；B(yǎng)殖場(chǎng)，見(jiàn)表3。榆溪河流域主要的生活污染源為榆林市污水處理廠和榆林市高新區(qū)閆莊則污水處理廠。工業(yè)污染源分為工業(yè)園區(qū)和經(jīng)濟(jì)開(kāi)發(fā)區(qū)，國(guó)控、省控工業(yè)污染源以及各水源地煤礦礦井水排放產(chǎn)生的污染源。根據(jù)現(xiàn)場(chǎng)調(diào)查及資料統(tǒng)計(jì)，規(guī)劃范圍內(nèi)有規(guī)?；B(yǎng)殖場(chǎng)30個(gè)。

(2)現(xiàn)狀非點(diǎn)源負(fù)荷估算。根據(jù)榆溪河流域的污染源調(diào)查，非點(diǎn)源負(fù)荷包括農(nóng)業(yè)非點(diǎn)狀污染源、分散畜禽養(yǎng)殖、農(nóng)村居民生活污染、漁業(yè)非點(diǎn)狀污染源，見(jiàn)表3。榆溪河流域耕地以種植玉米、高粱、大豆等為主。農(nóng)村以分散家庭養(yǎng)殖為主，家庭養(yǎng)殖業(yè)的條件較差，地面無(wú)防滲排水設(shè)施，動(dòng)物的糞便不經(jīng)任何處理隨意堆放。農(nóng)村居民生活污染包括垃圾、污水及糞便等。漁業(yè)以養(yǎng)殖草魚(yú)為主。

表3 榆溪河流域污染源負(fù)荷估計(jì)統(tǒng)計(jì)匯總 t/a

(3)水源地污染物削減量估算。水源地污染物削減量按照等比例分配法進(jìn)行總量分配，即按照各個(gè)水源地的入河量在總?cè)牒恿恐兴嫉谋壤秊闄?quán)重進(jìn)行，按照小壕兔、頭道河則、金雞灘、色草灣、趙家灣、芹河、馬合、沙河、水掌灣及牛家梁-紅石峽10個(gè)水源地進(jìn)行分配，結(jié)果見(jiàn)表4、表5。

(4)制定污染削減措施。榆溪河流域想要實(shí)現(xiàn)“水清、河暢、岸綠、景美”的目標(biāo)，首先，加快對(duì)污水處理設(shè)施的提標(biāo)改造，完善榆溪河沿河鄉(xiāng)鎮(zhèn)的污水處理廠建設(shè)，開(kāi)展污水處理系統(tǒng)建設(shè)，污水處理后回用。其次，加強(qiáng)生活垃圾收集處理，在流域內(nèi)建立起完善的“戶分類、村收集、鎮(zhèn)中轉(zhuǎn)、區(qū)處理”的城鄉(xiāng)垃圾一體化處理模式。針對(duì)流域內(nèi)易截污位置采用人工方式進(jìn)行生活垃圾、建筑垃圾、漂浮物全面打撈。最后，對(duì)散養(yǎng)的畜禽養(yǎng)殖轉(zhuǎn)變?yōu)橐?guī)?；B(yǎng)殖，對(duì)養(yǎng)殖場(chǎng)實(shí)施污染治理保證治理達(dá)標(biāo)。

3 結(jié) 論

榆溪河流域氨氮的污染較COD更為嚴(yán)重。污染物COD和氨氮，削減量要求最大為非點(diǎn)源污染的分散畜禽養(yǎng)殖，其次點(diǎn)源污染中的工業(yè)污染源榆溪河流域COD和氨氮污染物均來(lái)自非點(diǎn)源污染，其中分散畜禽養(yǎng)殖占比重大；點(diǎn)源污染負(fù)荷中，工業(yè)污染源入河量占比重大。

表4 榆溪河水源地COD污染物削減量匯總 t/a

表5 榆溪河水源地氨氮污染物削減量匯總 t/a

榆溪河流域10個(gè)水源地中，COD入河量的高低次序?yàn)轳R合水源地>小壕兔水源地>金雞灘水源地，COD的削減量高低次序?yàn)轳R合水源地>小壕兔水源地>水掌溝水源地。氨氮入河量的高低次序?yàn)轳R合水源地>小壕兔水源地>金雞灘水源地，氨氮削減總量的高低次序?yàn)樗茰纤吹?馬合水源地>小壕兔水源地。

通過(guò)對(duì)榆溪河流域構(gòu)建TMDL模型，對(duì)該流域的點(diǎn)源、非點(diǎn)源污染負(fù)荷及10個(gè)水源地的污染負(fù)荷進(jìn)行了分配且明確了各指標(biāo)削減量，根據(jù)負(fù)荷分配采取污染控制措施,制定水污染治理方案。引入TMDL模型，結(jié)合客觀賦權(quán)法進(jìn)行權(quán)重分配，更加客觀地反映污染物的削減量及負(fù)荷指標(biāo)。其次，該模型還存在不確定的因素，該方法不能完全真實(shí)且公平地分配污染負(fù)荷，故該模型還有改進(jìn)的空間。針對(duì)榆溪河流域水污染治理方案進(jìn)行分析，明晰現(xiàn)狀的不足，并提出相應(yīng)的整改對(duì)策，這對(duì)榆溪河流域改善水污染現(xiàn)狀具有重要意義。

□