邱居華，呂康玄

(1.三峽大學(xué)(湖北)設(shè)計咨詢研究院有限責(zé)任公司，湖北 宜昌 443002；2.中交第二航務(wù)工程勘察設(shè)計院有限公司，湖北 武漢 430060)

0 前 言

長江是中國第一大河，水資源極其豐富，是我國重要的生態(tài)安全屏障區(qū)。近年來，長江流域水資源的開發(fā)程度日益增大，水環(huán)境受人類活動干擾程度日趨惡化，水環(huán)境污染問題突出。因此，進(jìn)行長江水質(zhì)污染的評價和確定污染源十分重要。

隨著水環(huán)境污染問題日益嚴(yán)重，水質(zhì)與污染源的合理評價與判斷逐漸成為研究熱點之一。目前，已有諸多學(xué)者對此問題開展了系列研究。Lin tao等[1]使用多元統(tǒng)計技術(shù)和綜合水質(zhì)指數(shù)(IWQI)對2018—2020年中國廣東省茅州河流域的地表水水質(zhì)進(jìn)行了評估。Thotagamuwa H.T.B.N等[2]結(jié)合土地開發(fā)利用與水質(zhì)的關(guān)系對地表水進(jìn)行了評價。Wu Tao[3]等利用水質(zhì)綜合評價模型對2008-2018年的伊龍湖水質(zhì)進(jìn)行了分析。除此之外，朱長軍等[4-9]均從不同角度提出了綜合評價模型。

另一方面，如何進(jìn)一步確定河流污染源，對河流生態(tài)環(huán)境治理與防控有著重要意義。高嵩等[10]以泰州老通揚運河引水工程為例,利用一維河網(wǎng)水動力水質(zhì)模型來分析論證新建泵站的必要性和規(guī)模。王偉等[11]系統(tǒng)介紹了國內(nèi)外廣泛運用的三類水質(zhì)耦合模型類型，并提出了展望。潘祥東等[12]基于一維水質(zhì)模型計算了鹽津河納污能力，為鹽津河當(dāng)?shù)厮h(huán)境治理提供參考和科學(xué)依據(jù)。周睿等[13]提出了將一維水質(zhì)模型應(yīng)用于分離監(jiān)測斷面點源污染負(fù)荷和非點源污染負(fù)荷的方法。張穎等[14]對里下河地區(qū)水質(zhì)模型進(jìn)行了研究。然而，以上方法在對污染源進(jìn)行確定與分析時未給出詳細(xì)的隸屬度函數(shù)及權(quán)重計算方法。

本文對隸屬度函數(shù)及權(quán)重的計算進(jìn)行了深入的分析，對長江水質(zhì)綜合評價方法適用性進(jìn)行了驗證，并基于一維水質(zhì)模型對污染源進(jìn)行了確定與分析，為水環(huán)境治理提供方法依據(jù)。

1 模糊綜合評價法

根據(jù)模糊綜合評價方法，可以得到水質(zhì)綜合評價模型如下：

FCI=B×S

(1)

式中：B為各種等級水質(zhì)出現(xiàn)的可能性；S為描述各類水質(zhì)危害程度的水質(zhì)標(biāo)準(zhǔn)類別向量。可分別利用下式求出：

B=W×R

(2)

S=(S1,S2,S3,S4,S5,S6)

(3)

公式(1)～(3)中：W為水質(zhì)評判因素權(quán)重向量；R為隸屬度組成的模糊關(guān)系矩陣；S這是一個主觀因素很強的量，可以人為確定。已知Ⅰ、Ⅱ、Ⅲ類為可飲用水即品質(zhì)最好危害程度最低，此處根據(jù)日常經(jīng)驗可以取S=(1,2,3,4,5,6)。對于式(2)中的W及R是本文所要考慮的重點，最終構(gòu)造權(quán)重矩陣及隸屬度模糊關(guān)系矩陣如下：

對于權(quán)重矩陣，其超標(biāo)倍數(shù)越多權(quán)重越大，因此構(gòu)造其計算方式如下：

W=(w1,w2,w3,w4)

(4)

其中，

表1 水質(zhì)標(biāo)準(zhǔn)值 /(mg·L-1)

另一方面，對于由隸屬度組成的模糊關(guān)系矩陣可得：

(5)

式中:i對j的隸屬度 表示第 種污染物的環(huán)境質(zhì)量數(shù)值可以被評價為第j類環(huán)境質(zhì)量的可能性。例如：r11表示DO對Ⅰ類水的隸屬程度。進(jìn)而，以表1中污染物的順序可求得其隸屬度的計算公式如下：

對i=2,3,有j=1時為：

(6)

j=2,3,4,5時為：

(7)

j=6時為：

(8)

對i=4,j=1,2,3,4,5時為：

(9)

j=6時為：

(10)

當(dāng)i=1時，與上面類似，只不過它是逆向趨勢，這里不再贅述。這里kij是《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》中4個主要項目標(biāo)準(zhǔn)限值。

就本文而言，其水質(zhì)影響因素集即為U={DO,CODMn,NH3-N,pH}，而水質(zhì)判斷集即為水質(zhì)的不同評價等級組成的集合。V={Ⅰ類，Ⅱ類，Ⅲ類，Ⅳ類，Ⅴ類}。

2 基于一維水質(zhì)模型的污染源確定方法

已知干流上各觀測點每個月份的水流速，但由于水流速度是變化的，它的平均速度很難求出，所以根據(jù)上一個站點的水流速度和下一個站點的水流速度來估計這兩個站點之間的平均速度，即：

(11)

(12)

各觀測點每秒流過長江水中污染物的總量為污染物濃度與單位時間該觀測點流量Qij的乘積，所得表達(dá)式如下：

mij=nij×Qij

(13)

公式(15)中:mij為每秒流過長江水中污染物的總量，mol；nij為污染物濃度，mol/m3；Qij為單位時間該觀測點流量，m3/s。

基于GIS技術(shù)的一維水質(zhì)模型，建立一維河流的污染物濃度與降解系數(shù)的偏微分方程為：

(14)

公式(16)中:t為時間，s;C為污染物濃度， mol/m3;k為降解系數(shù);S為污染物源項;E為擴散系數(shù);x為沿河流方向位置。

當(dāng)流速較大時，污染物的推移作用遠(yuǎn)大于擴散作用，故可忽略擴散系數(shù)不計，同時若假定個污染源排放為連續(xù)穩(wěn)定，且 ，即污染物濃度變化與時間無關(guān)，此時一維模型可簡化為：

(15)

給定邊界條件，則下游斷面處污染物與降解系數(shù)的函數(shù)關(guān)系為：

f(t)=mij×e-αtij

(16)

因為一個觀測站(地區(qū))的水質(zhì)污染主要來自于本地區(qū)的排污和上游的污水，所以求出上一個站點的污染物降解后的剩余污染物后，可以兩個站點之間污染物增加情況，所得表達(dá)式如下：

mi(j+1)-f(t)=ni(j+1)×Qij-mij×e-αtij

(17)

3 長江流域水質(zhì)評價及污染源確定分析

3.1 數(shù)據(jù)來源情況

長江流域水資源數(shù)據(jù)和水污染數(shù)據(jù)均來自2011—2017年《中國統(tǒng)計年鑒》及長江干流7個觀測站點的觀測數(shù)據(jù)。

3.2 綜合評價結(jié)果與分析

利用上述模糊綜合評價模型以綜合評定值排序的結(jié)果得出各地區(qū)水污染程度，最終可以判斷各地區(qū)水質(zhì)的污染狀況，如圖1所示。

圖1中橫坐標(biāo)1-17代表表2對應(yīng)的17個區(qū)域。其中，各個地區(qū)水質(zhì)的污染狀況由小到大的排序如表2所示。

表2 污染狀況由小到大排序

由表2可以看出：各觀測城市所在的江段的水質(zhì)污染的情況，水質(zhì)最差的是江蘇南京林；其次是江西南昌滁槎；第三位的是四川樂山岷江大橋地區(qū)。水質(zhì)最好的是湖南岳陽城陵磯，其次是湖南長沙新港地區(qū)，第三位的是湖南岳陽岳陽樓地區(qū)?？梢姾系貐^(qū)水質(zhì)保護(hù)效果較好，江蘇與江西地區(qū)水質(zhì)保護(hù)效果較差。

3.3 污染源確定結(jié)果與分析

引入一維水質(zhì)模型，上游污染物呈衰減形式對下游影響，用本地區(qū)污染物的量減去上游對下游的影響量及得到本地區(qū)排入污染物的量，依次求得各地區(qū)排入的污水量進(jìn)行比較，兩種污染物排入量如表3及表4所示。

表3 高錳酸鹽指數(shù)總量 /mol

表4 氨氮總量 /mol

平均值可以反映事物的一般規(guī)律,以高錳酸鹽和氨氮的平均值變化觀察污染物主要集中地,如圖2所示。

由圖2可以看出，數(shù)據(jù)值越大說明排放的污染物越多，高錳酸鹽和氨氮的排放量都是在湖北宜昌至湖南岳陽這一區(qū)段數(shù)量最大，江西九江河西水廠至安徽安慶次之，由此可下結(jié)論說長江干流近一年多主要污染物高錳酸鹽指數(shù)和氨氮的污染源主要在湖北宜昌至湖南岳陽、江西九江河西水廠至安徽安慶之間。

通過以上對長江流域污染程度評價及污染源的確定與分析，結(jié)果和后期相關(guān)報道[15-16]中關(guān)于長江流域污染程度及污染物種類高度吻合，由此可見，本文介紹的方法對水質(zhì)污染的評價和污染源確定問題具有較好的工程適用性。

4 結(jié) 論

(1) 本文采用模糊數(shù)學(xué)的概念建立了模糊綜合評價模型，基于一維水質(zhì)模型提出了一種水質(zhì)評價及污染源判別的方法,并分析得出長江干流近一年多主要污染物高錳酸鹽指數(shù)和氨氮的污染源主要在湖北宜昌至湖南岳陽、江西九江河西水廠至安徽安慶之間。

(2) 結(jié)合2011—2017年《中國統(tǒng)計年鑒》及長江干流7個觀測站點的觀測數(shù)據(jù)對本文方法進(jìn)行了驗證，結(jié)果表明該方法可準(zhǔn)確評價水質(zhì)污染程度并有效確定污染源，具有較高的工程適用性，可為長江大保護(hù)的水質(zhì)考核及污染物的確定與治理提供數(shù)據(jù)支撐與理論參考。