殷 飛，金世佳

(1.吉林農(nóng)業(yè)科技學(xué)院，吉林 吉林 132101；2. 沈陽農(nóng)業(yè)大學(xué)水利學(xué)院，沈陽 110866)

隨著我國經(jīng)濟(jì)社會(huì)的快速發(fā)展和人口的增長，河流以及湖庫的水環(huán)境生態(tài)系統(tǒng)出現(xiàn)了惡化的趨勢，局部地區(qū)的經(jīng)濟(jì)社會(huì)發(fā)展受到阻礙，甚至成為影響和制約區(qū)域可持續(xù)發(fā)展的瓶頸問題之一。國家環(huán)保部在2014年《中國環(huán)境狀況公報(bào)》中指出，長江、黃河等7大流域和浙閩片河流、西北諸河、西南諸河共有968個(gè)國控水質(zhì)監(jiān)測點(diǎn)，其中水質(zhì)處于Ⅳ～Ⅴ類和劣Ⅴ類的比例占27.7%和9.2%[1]，我國地表水體依然為輕度污染。嚴(yán)重污染的水體會(huì)導(dǎo)致其喪失水生態(tài)功能和使用價(jià)值，從而威脅生態(tài)平衡，危害人類健康。因此，合理、正確地監(jiān)測和評(píng)價(jià)地表水水質(zhì)，不但可以為流域水體污染防治提供重要的理論依據(jù)，而且對(duì)于生活、生產(chǎn)用水和水生態(tài)安全極其重要。

對(duì)于水質(zhì)評(píng)價(jià)方法的研究非常多，比如綜合污染指數(shù)評(píng)價(jià)法[2,3]、模糊數(shù)學(xué)綜合評(píng)判[4-6]、人工神經(jīng)網(wǎng)絡(luò)法[7,8]、灰色關(guān)聯(lián)度法[9]、層次分析法[10,11]、綜合水質(zhì)標(biāo)識(shí)指數(shù)法[12,13]等，這些方法都存在著這樣或那樣的缺點(diǎn)。綜合污染指數(shù)法只能進(jìn)行定性評(píng)價(jià)；模糊綜合評(píng)價(jià)法以水質(zhì)樣本的隸屬度最大來對(duì)應(yīng)水質(zhì)，定量計(jì)算相對(duì)復(fù)雜；人工神經(jīng)網(wǎng)絡(luò)法因其樣本的選擇不固定而導(dǎo)致定性評(píng)價(jià)結(jié)果不太直觀；灰色關(guān)聯(lián)度法評(píng)價(jià)有可能掩蓋部分超標(biāo)污染因子對(duì)綜合評(píng)價(jià)結(jié)果的影響；層次分析法因其評(píng)價(jià)指標(biāo)體系不易確定問題可能影響評(píng)價(jià)結(jié)果，計(jì)算也較繁瑣。其中以綜合水質(zhì)標(biāo)識(shí)指數(shù)法在水質(zhì)評(píng)價(jià)中應(yīng)用起來最為方便，并且可以比較簡單、準(zhǔn)確地對(duì)河流的綜合水質(zhì)進(jìn)行定性、定量評(píng)價(jià)[14]，但是在應(yīng)用中因四舍五入時(shí)會(huì)出現(xiàn)小數(shù)點(diǎn)后第1位不為個(gè)數(shù)的現(xiàn)象，因而本研究對(duì)徐祖信的綜合水質(zhì)標(biāo)識(shí)指數(shù)法進(jìn)行改進(jìn)，進(jìn)而對(duì)松花江流域吉林省段 2005-2014 年的水質(zhì)監(jiān)測數(shù)據(jù)進(jìn)行評(píng)價(jià)并分析其水質(zhì)變化的趨勢，以期為吉林省水環(huán)境治理與保護(hù)提供理論依據(jù)。

1 研究區(qū)概況與研究方法

1.1 研究區(qū)概況

松花江流域(N 41°42′～51°38′、E 119°52′～132°31′)面積55.68 萬km2，跨遼、吉、黑、內(nèi)蒙古4省(自治區(qū))，見圖1，包括哈爾濱、長春、吉林等105個(gè)市縣(旗)轄區(qū)，多年平均降水量為500 mm左右，汛期的降水量約占全年降水量的60%～80% 。其中吉林省松花江流域面積 7.80 萬km2，占松花江流域面積的14.01%，流經(jīng)長春市、吉林市等10個(gè)地市43個(gè)縣區(qū)，經(jīng)濟(jì)發(fā)展水平較高，是吉林省主要的商品糧產(chǎn)區(qū)和工業(yè)分布區(qū)，主要有吉林化工和長春汽車等重工業(yè)，因而也是松花江流域的主要污染源地，具體位置見圖2。

圖1 松花江流域位置[15]Fig.1 Location map of Songhuajiang river basin

圖2 研究區(qū)位置示意圖[16]Fig.2 Sketch map of location of the study area

松花江流域(吉林省段)內(nèi)工業(yè)布局密集，人口集中，隨著社會(huì)經(jīng)濟(jì)的發(fā)展，工業(yè)及生活污水排放對(duì)流域內(nèi)的水體造成了嚴(yán)重程度的污染，其中長春、吉林2市對(duì)松花江流域(吉林省段)環(huán)境類型區(qū)的污染效應(yīng)系數(shù)分別達(dá)0.55 和0.22，2市共占 0.77[17]，而水體環(huán)境質(zhì)量的好壞直接影響到吉林省經(jīng)濟(jì)社會(huì)的發(fā)展和水生態(tài)的安全。因此本文選取國控水質(zhì)監(jiān)測斷面吉林省6個(gè)斷面(其中2005-2011年3個(gè)斷面，2012-2014年6個(gè)斷面)，分別為延邊圈河、吉林長春松花江、吉林白城白沙灘、吉林敦化新甸、吉林松原、吉林溪浪口，前3個(gè)斷面為共有斷面。

此次水質(zhì)評(píng)價(jià)指標(biāo)主要為溶解氧、高錳酸鹽指數(shù)、氨氮，所有水質(zhì)監(jiān)測數(shù)據(jù)均來源于環(huán)保部數(shù)據(jù)中心的全國主要流域重點(diǎn)斷面水質(zhì)自動(dòng)監(jiān)測周報(bào)，采用各項(xiàng)評(píng)價(jià)因子的年均值對(duì)松花江流域(吉林省段)的水質(zhì)進(jìn)行評(píng)價(jià)，同時(shí)采用秩相關(guān)系數(shù)法來判定水質(zhì)的變化趨勢。

1.2 水質(zhì)標(biāo)識(shí)指數(shù)法

水質(zhì)評(píng)價(jià)方法已有許多研究成果，本研究采用改良的綜合水質(zhì)標(biāo)識(shí)指數(shù)法，具體的步驟為：

(1)確定各水質(zhì)參數(shù)的單因子水質(zhì)標(biāo)識(shí)指數(shù)P，表示為：Pi=X1·X2X3，總的來講Pi數(shù)值越大，表示第i項(xiàng)監(jiān)測指標(biāo)的水質(zhì)污染越嚴(yán)重。

(2)在單因子水質(zhì)標(biāo)識(shí)指數(shù)P的基礎(chǔ)上，確定綜合水質(zhì)標(biāo)識(shí)指數(shù)I：I=X1·X2X3X4，總的來講I數(shù)值越大，表示綜合水質(zhì)越差。

具體各字母含義見徐祖信等的研究[12,13]。

1.3 水質(zhì)變化趨勢檢驗(yàn)

由于水質(zhì)監(jiān)測值會(huì)受到許多因素的影響，常常表現(xiàn)出波動(dòng)性，時(shí)而趨于上升，時(shí)而趨于下降[18]，因而很難確定水質(zhì)的變化趨勢是上升還是下降，水質(zhì)變化趨勢的分析就必須依賴于數(shù)理統(tǒng)計(jì)學(xué)原理來進(jìn)行趨勢檢驗(yàn)分析[19]。

松花江流域吉林省段2005-2014年水質(zhì)變化趨勢分析采用的是我國環(huán)保部公布的《地表水環(huán)境質(zhì)量評(píng)價(jià)辦法》中推薦的Daniel的趨勢檢驗(yàn)，使用Spearman秩相關(guān)系數(shù)法來判定。它的特點(diǎn)是將各監(jiān)測點(diǎn)水質(zhì)數(shù)據(jù)轉(zhuǎn)化為依次的排名，而不用考慮實(shí)際值。采用此法可以全面了解水質(zhì)的變化趨勢，從而判別水質(zhì)是否好轉(zhuǎn)、加重或水質(zhì)變化平穩(wěn)[20]，從而避免單因子評(píng)價(jià)法造成的一些不合理現(xiàn)象。

具體過程如下：給定時(shí)間Y1,…,YN和相應(yīng)的監(jiān)測值X(即年均值C1，…，CN)，由小到大將順序排好，這種排列的次序就稱為秩。具體按下式計(jì)算：

式中：di表示變量Xi與Yi的差值；Xi表示Y1到Y(jié)N期間內(nèi)水質(zhì)監(jiān)測數(shù)據(jù)由小到大排列的序號(hào)；Yi表示按時(shí)間排列的序號(hào)。

將|rs|同臨界值Wp進(jìn)行比較。rs>Wp表明水質(zhì)變化趨勢具有顯著的意義：如果rs<0，則表明水質(zhì)趨于好轉(zhuǎn)；如果rs>0，則表明水質(zhì)趨于惡化。rs≤Wp表明水質(zhì)的變化不顯著，說明水質(zhì)變化趨于穩(wěn)定。

2 水質(zhì)評(píng)價(jià)及其變化趨勢分析

2.1 改良單因子水質(zhì)標(biāo)識(shí)指數(shù)法評(píng)價(jià)

以延邊圈河溶解氧(6.34 mg/L)為例來說明單因子水質(zhì)標(biāo)識(shí)指數(shù)的計(jì)算，對(duì)照中華人民共和國地表水環(huán)境質(zhì)量Ⅱ類標(biāo)準(zhǔn)，X2為2；X2代表6.34 mg/L在Ⅱ類水中所處的位置，7.50、6.00分別代表溶解氧在Ⅱ類水中的上限值和下限值，則X2=100 (7.50-6.34)/(7.50-6.00)=77；X3表示水質(zhì)類別(Ⅱ)好于功能區(qū)類別(Ⅲ)，則X3=0，即單因子水質(zhì)標(biāo)識(shí)指數(shù)為2.770。以此類推，對(duì)松花江流域吉林省段在2005-2014年各斷面水質(zhì)進(jìn)行評(píng)價(jià)，見表1。結(jié)果表明：參與水質(zhì)評(píng)價(jià)的3項(xiàng)水質(zhì)指標(biāo)在10 a里僅有白城白沙灘和延邊圈河2個(gè)斷面的水質(zhì)指標(biāo)存在不同程度的超標(biāo)，超標(biāo)的水質(zhì)指標(biāo)為高錳酸鹽指數(shù)、溶解氧；其中白城白沙灘斷面僅在2005年溶解氧超標(biāo)，延邊圈河斷面高錳酸鹽指數(shù)在10 a中均超標(biāo)，超標(biāo)頻次為10次，其余4個(gè)斷面都達(dá)標(biāo)(未知數(shù)據(jù)年份除外)。從表1可知，高錳酸鹽指數(shù)在松花江流域吉林省段水質(zhì)中超標(biāo)頻次最高，程度較深，是影響水質(zhì)的主要污染指標(biāo)。

表1 2005-2014年單因子水質(zhì)標(biāo)識(shí)指數(shù)法對(duì)各斷面水質(zhì)的評(píng)價(jià)結(jié)果Tab.1 Assessment of the section water quality based on the single factor water quality identification index in 2005-2014

注:“-”表示該年份無自動(dòng)監(jiān)測的水質(zhì)數(shù)據(jù)。

2.2 改良綜合水質(zhì)標(biāo)識(shí)指數(shù)法評(píng)價(jià)

單因子水質(zhì)標(biāo)識(shí)指數(shù)能夠很好地判別水體中主要的污染物，但是進(jìn)行綜合評(píng)價(jià)時(shí)不應(yīng)以污染最嚴(yán)重的指標(biāo)來判定水質(zhì)的類別，需要綜合考慮各種污染物對(duì)水質(zhì)的影響。大量的研究表明，綜合水質(zhì)標(biāo)識(shí)指數(shù)法的評(píng)價(jià)結(jié)果與灰色關(guān)聯(lián)法等的評(píng)價(jià)結(jié)果擬合程度最高，且運(yùn)算結(jié)果更為穩(wěn)定、計(jì)算過程更為簡單、實(shí)用性更好[21,22]。因而，本文采用綜合水質(zhì)標(biāo)識(shí)指數(shù)法對(duì)松花江流域吉林省段的水質(zhì)進(jìn)行綜合定性的評(píng)價(jià)，評(píng)價(jià)結(jié)果見表2。結(jié)果表明： 2005-2014年，僅有延邊圈河斷面在2006年綜合評(píng)價(jià)時(shí)為Ⅳ類水質(zhì)，未達(dá)到功能區(qū)Ⅲ類水質(zhì)標(biāo)準(zhǔn)，其余各斷面水質(zhì)不僅達(dá)到而且超過功能區(qū)的標(biāo)準(zhǔn)。在2014年，綜合水質(zhì)標(biāo)識(shí)指數(shù)延邊圈河斷面最高(3.001 0)，溪浪口最低(2.280 0)，從優(yōu)到劣排序見圖3。

圖3 2014年各斷面綜合水質(zhì)的變化Fig.3 The comprehensive water quality of each section in 2014

2.3 典型方法的比較分析

單因子評(píng)價(jià)法在我國水質(zhì)監(jiān)測公報(bào)中被廣泛采用, 計(jì)算簡便，缺點(diǎn)是它常常以水質(zhì)最差因子所屬的水質(zhì)類別來判定綜合水質(zhì)；污染指數(shù)評(píng)價(jià)法不僅考慮單項(xiàng)最差的水質(zhì)指標(biāo)，而且考慮其他水質(zhì)評(píng)價(jià)指標(biāo)，通過采用加權(quán)等數(shù)學(xué)方法得到各單項(xiàng)污染指數(shù)的綜合指數(shù)來判定水體的綜合水質(zhì)。具體比較結(jié)果見表3。

表2 2005-2014年綜合水質(zhì)標(biāo)識(shí)指數(shù)法對(duì)各斷面水質(zhì)的評(píng)價(jià)結(jié)果Tab.2 Assessment of the section water quality based on the comprehensive water quality identification index in 2005-2014

表3 評(píng)價(jià)方法結(jié)果比較Tab.3 The comparison of evaluation methods

分析表3可知，單因子評(píng)價(jià)法的評(píng)價(jià)結(jié)果突出了污染最嚴(yán)重因子對(duì)水質(zhì)的影響，掩蓋了其他因子的影響，表現(xiàn)對(duì)水質(zhì)判別過為保護(hù)，因此，比較適合定性來確定水體中主要污染物。等權(quán)重污染指數(shù)法則通過比較計(jì)算結(jié)果是否大于1來判斷達(dá)到相應(yīng)水質(zhì)標(biāo)準(zhǔn),但是在確定綜合水質(zhì)類別時(shí)則顯得無能為力。相比較以上2種評(píng)價(jià)方法，綜合水質(zhì)標(biāo)識(shí)指數(shù)法則能直觀反映綜合水質(zhì)的級(jí)別，是否達(dá)到功能區(qū)目標(biāo)以及劣于功能區(qū)的單項(xiàng)指標(biāo)個(gè)數(shù)，計(jì)算過程簡單，無需復(fù)雜的矩陣運(yùn)算和編程語言，具有很強(qiáng)的操作性和可推廣性，優(yōu)勢明顯，在一定程度上得到了廣泛的應(yīng)用。

2.4 水質(zhì)變化趨勢分析

對(duì)松花江流域吉林省國控?cái)嗝娴乃|(zhì)進(jìn)行常規(guī)監(jiān)測和分析、評(píng)價(jià)非常重要，通過水質(zhì)變化趨勢的分析，不僅能夠判定松花江流域目前實(shí)行的水污染治理措施是否得當(dāng)，而且對(duì)進(jìn)一步采取何種有效的管理與治理措施具有重要的意義。本文采用環(huán)境保護(hù)部常用的Spearman秩相關(guān)系數(shù)法對(duì)2005-2014年的水質(zhì)變化趨勢進(jìn)行分析，前提是使用該方法需要足夠的數(shù)據(jù)，一般至少具備4對(duì)及以上的數(shù)據(jù)[23]，因?yàn)樗稍?、敦化新甸、溪浪口斷面只? a數(shù)據(jù)，所以只能選取其余3個(gè)斷面數(shù)據(jù)進(jìn)行分析，結(jié)果見表4。

表4結(jié)果表明，松花江流域吉林省3個(gè)斷面均為|rs|>Wp0.05，說明這些斷面的水質(zhì)在10 a間呈顯著好轉(zhuǎn)趨勢，從而說明流域治理措施效果明顯。主要措施為吉林省“十一五”、“十二五”期間相繼投資50.39 億元和97.87億元建設(shè)了一批重點(diǎn)項(xiàng)目防治松花江水污染，同時(shí)加以生態(tài)恢復(fù)，確保松花江得以休養(yǎng)生息。

表4 2005-2014年松花江流域吉林省各斷面水質(zhì)變化趨勢分析Tab.4 Water quality trend analysis of song hua river basin on jilin province section in 2005-2014

注：Wp0.05=0.56。

3 結(jié) 語

(1)本文采用改良的單因子和綜合水質(zhì)標(biāo)識(shí)指數(shù)法對(duì)松花江流域吉林省段進(jìn)行評(píng)價(jià)，從而實(shí)現(xiàn)對(duì)河流水質(zhì)進(jìn)行簡單而準(zhǔn)確的定性和定量評(píng)價(jià)。

(2)通過對(duì)單因子評(píng)價(jià)法、污染指數(shù)法以及綜合水質(zhì)標(biāo)識(shí)指數(shù)法比較，表明綜合水質(zhì)標(biāo)識(shí)指數(shù)法優(yōu)勢明顯，比較實(shí)用，易于在其他河流等水環(huán)境評(píng)價(jià)中推廣使用；利用評(píng)價(jià)的水質(zhì)因子太少，不能全面反映真實(shí)的水質(zhì)情況，需要進(jìn)一步收集數(shù)據(jù)。

[1] 中華人民共和國環(huán)境保護(hù)部．2014年中國環(huán)境狀況公報(bào)[Z]．北京：中華人民共和國環(huán)境保護(hù)部，2014．

[2] 劉 碩, 朱建平, 蔣火華. 對(duì)幾種環(huán)境質(zhì)量綜合指數(shù)評(píng)價(jià)方法的探討[J].中國環(huán)境監(jiān)測, 2003,15(5):33-37．

[3] 陳詠淑,吳甫成,呂煥哲,等.近20 a來湘江水質(zhì)變化分析[J].長江流域資源與環(huán)境,2004,13(5):508-512．

[4] 孫靖南, 鄒志紅, 任廣平. 模糊綜合評(píng)價(jià)在天然水體水質(zhì)評(píng)價(jià)中的應(yīng)用研究[J].環(huán)境污染治理技術(shù)與設(shè)備,2005,6(2):45-48．

[5] 傅金祥, 陳 喆,馬興冠, 等. 改良模糊綜合評(píng)價(jià)法在水質(zhì)評(píng)價(jià)中的應(yīng)用[J].環(huán)境工程,2011,29(6):120-127．

[6] 張水珍,劉 玲.模糊綜合評(píng)判法在水質(zhì)評(píng)價(jià)中的應(yīng)用----以松花江流域?yàn)槔齕J].環(huán)境科學(xué)與管理, 2011,36(3):163-165.

[7] 劉國東,黃川友,丁 晶.水質(zhì)綜合評(píng)價(jià)的人工神經(jīng)網(wǎng)絡(luò)模型[J].中國環(huán)境科學(xué),1998,18(6):514-517．

[8] 郭勁松,王海霞,龍騰銳.人工神經(jīng)網(wǎng)絡(luò)在水質(zhì)規(guī)劃和管理中的應(yīng)用[J].重慶環(huán)境科學(xué),2002,(4):69-72．

[9] 秦昌波,鄭炳輝,秦延文,等.渤海灣天津段海岸帶水環(huán)境質(zhì)量灰色關(guān)聯(lián)度評(píng)價(jià)[J].環(huán)境科學(xué)研究,2006,19(6):94-99．

[10] 龐振凌,常紅軍,李玉英,等.層次分析法對(duì)南水北調(diào)中線水源區(qū)的水質(zhì)評(píng)價(jià)[J].生態(tài)學(xué)報(bào),2008,28(4)：1 810-1 819.

[11] 喬 俊,岳芳寧,羅水明,等.應(yīng)用層次分析法綜合評(píng)價(jià)河流沉積物污染狀況[J].環(huán)境科學(xué)研究,2012,25(7):840-844.

[12] 徐祖信．我國河流單因子水質(zhì)標(biāo)識(shí)指數(shù)評(píng)價(jià)方法研究[J].同濟(jì)大學(xué)學(xué)報(bào):自然科學(xué)版, 2005,33(3):321-325．

[13] 徐祖信．我國河流綜合水質(zhì)標(biāo)識(shí)指數(shù)評(píng)價(jià)方法研究[J].同濟(jì)大學(xué)學(xué)報(bào):自然科學(xué)版,2005,33(4):482-488．

[14] 趙海生,喬 俊.丹河晉城段水質(zhì)評(píng)價(jià)及變化趨勢分析[J].環(huán)境科學(xué)與管理,2012,37(7):60-62.

[15] 盧亞靈，蔣洪強(qiáng)，吳文俊，等. 流域經(jīng)濟(jì)環(huán)境綜合形勢研究：以松花江流域?yàn)槔齕J].環(huán)境科學(xué)與技術(shù)，2015，38(10)：175-181.

[16] 楊麗花，佟連軍.吉林省松花江流域經(jīng)濟(jì)發(fā)展與水環(huán)境質(zhì)量的動(dòng)態(tài)耦合及空間格局[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2013，24(2):503-510.

[17] 覃雪波，馬立新，孫海彬．松花江水污染及其防治對(duì)策[J].農(nóng)業(yè)環(huán)境與發(fā)展，2007，25(6):76-79．

[18] 梅永進(jìn)，陸英英.用秩相關(guān)系數(shù)法分析沙溪河沙縣段水質(zhì)變化[J].閩江學(xué)院學(xué)報(bào),2002,(6):112-119.

[19] 吳忠勇.環(huán)境監(jiān)測學(xué)概論[M].長沙：湖南大學(xué)出版社,1987:42-46.

[20] 中華人民共和國環(huán)境保護(hù)部.地表水環(huán)境質(zhì)量評(píng)價(jià)辦法(試行)[Z].北京：中華人民共和國環(huán)境保護(hù)部，2011.

[21] 安樂生,趙全升,劉貫群,等.代表性水質(zhì)評(píng)價(jià)方法的比較研究[J].中國環(huán)境監(jiān)測,2010,26(5):47-51．

[22] 尹海龍，徐祖信．河流綜合水質(zhì)評(píng)價(jià)方法比較研究[J].長江流域資源與環(huán)境,2008,(5):729-733．

[23] 張海嬌，李 峰，劉瑞瓊.曲靖市東部山地巖溶水源地水質(zhì)現(xiàn)狀評(píng)價(jià)及演變趨勢分析[J].資源環(huán)境與工程,2008,22(5):524-528.