高倩倩，楊 賢，吳加祿，楊永欽，傅德瑜

（1.四川省德陽(yáng)生態(tài)環(huán)境監(jiān)測(cè)中心站，四川 德陽(yáng) 618000；2.德陽(yáng)市生態(tài)環(huán)境局，四川 德陽(yáng) 618000）

石亭江是德陽(yáng)市境內(nèi)重要的河流，屬于沱江流域，發(fā)源于四川省什邡市紅白鎮(zhèn)，屬于沱江流域，全長(zhǎng)122 km，流域面積1512 km2，涉及的行政區(qū)域包括綿竹市、什邡市、旌陽(yáng)區(qū)、廣漢市。金輪斷面是石亭江流域的重要斷面，作為石亭江什邡綿竹段控制單元出境控制節(jié)點(diǎn)和石亭江廣漢旌陽(yáng)入境控制節(jié)點(diǎn)，金輪斷面水質(zhì)狀況的控制對(duì)石亭江流域水質(zhì)環(huán)境質(zhì)量的改善起到關(guān)鍵作用。因此，對(duì)石亭江金輪斷面水環(huán)境質(zhì)量現(xiàn)狀進(jìn)行監(jiān)測(cè)和分析，及發(fā)展趨勢(shì)的預(yù)測(cè)，可為改善該流域的水環(huán)境質(zhì)量提供數(shù)據(jù)依據(jù)和支撐。

Spearman秩相關(guān)系數(shù)法作為一種常用污染變化趨勢(shì)的定量分析方法，通過(guò)對(duì)各污染因子監(jiān)測(cè)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，可表明其在評(píng)價(jià)時(shí)段內(nèi)變化趨勢(shì)。環(huán)境監(jiān)測(cè)中，許多因素影響環(huán)境質(zhì)量，有些因素是不明確的?；疑到y(tǒng)模型GM （1，1）趨勢(shì)因子a法建立的定量預(yù)測(cè)模型，可以對(duì)污染因子定性分析，預(yù)測(cè)未來(lái)時(shí)間段內(nèi)污染因子的變化趨勢(shì)，反映出未來(lái)時(shí)段內(nèi)需要重點(diǎn)關(guān)注的污染因子，在水質(zhì)環(huán)境監(jiān)測(cè)中，使用灰色模型來(lái)研究是相適應(yīng)的[1]。

通過(guò)對(duì)金輪斷面近幾年監(jiān)測(cè)數(shù)據(jù)的分析，用Spearman秩相關(guān)系數(shù)法評(píng)價(jià)各污染因子在評(píng)價(jià)時(shí)段內(nèi)的變化趨勢(shì)，用灰色系統(tǒng)模型GM（1，1）趨勢(shì)因子a法預(yù)測(cè)未來(lái)時(shí)段內(nèi)污染指標(biāo)的變化趨勢(shì)。

1 石亭江金輪斷面水環(huán)境質(zhì)量評(píng)價(jià)

1.1 評(píng)價(jià)因子和標(biāo)準(zhǔn)

根據(jù)2016年～2021年市站監(jiān)測(cè)數(shù)據(jù)，選用高錳酸鹽指數(shù)（CODMn）、五日生化需氧量（BOD5）、氨氮（NH3-N）、總磷（TP）進(jìn)行趨勢(shì)預(yù)測(cè)和評(píng)價(jià)分析。根據(jù)德陽(yáng)市地表水環(huán)境質(zhì)量功能區(qū)劃分類，金輪斷面屬于Ⅲ類水域功能區(qū)，應(yīng)執(zhí)行《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》（GB 3838-2002）Ⅲ類水質(zhì)標(biāo)準(zhǔn)。

1.2 水環(huán)境質(zhì)量評(píng)價(jià)方法

金輪斷面水質(zhì)環(huán)境質(zhì)量現(xiàn)狀評(píng)價(jià)采用Spearman秩相關(guān)系數(shù)法和灰色系統(tǒng)模型GM（1，1）趨勢(shì)因子a法進(jìn)行分析評(píng)價(jià)[2]。

1.2.1 Spearman秩相關(guān)系數(shù)法

式中：di為變量Xi與Yi的差值；Xi為周期I-N從小到大按濃度值排序的序號(hào)；Yi為按時(shí)間排序的序號(hào)。

對(duì)將秩相關(guān)系數(shù)Rs的絕對(duì)值同臨界值Wp進(jìn)行比較，當(dāng)|Rs|≥Wp時(shí)，則表明變化趨勢(shì)有顯著意義，如果Rs是負(fù)值，則表明在評(píng)價(jià)時(shí)段內(nèi)有關(guān)統(tǒng)計(jì)量指標(biāo)變化呈下降或好轉(zhuǎn)趨勢(shì)；如果Rs為正值，則表明在評(píng)價(jià)時(shí)段內(nèi)有關(guān)統(tǒng)計(jì)量指標(biāo)變化呈上升或加重趨勢(shì)。

1.2.2 灰色系統(tǒng)模型GM（1，1）趨勢(shì)因子a法

構(gòu)造Yn向量如下：

用最小二乘法求取待辨識(shí)參數(shù)列S：

式中：T表示矩陣的轉(zhuǎn)置，-1表示矩陣的逆。

建立各種污染因子的預(yù)測(cè)模型：

對(duì)（1）式預(yù)測(cè)模型進(jìn)行趨勢(shì)判斷：

對(duì)于灰色系統(tǒng)模型GM（1，1）趨勢(shì)因子α法，當(dāng) 大于0，則e-1e-a小于1，即，表明污染因子的濃度值隨時(shí)間推移而降低，且a值越大，降低趨勢(shì)越明顯。當(dāng) 小于0，則e-a大于1，即，表明污染因子的濃度隨時(shí)間推移而變大，且 值越小，變大趨勢(shì)越明顯。當(dāng) 等于0，則e等于1，即，表明污染因子的濃度不隨時(shí)間變化。

2 水環(huán)境質(zhì)量監(jiān)測(cè)值和結(jié)果評(píng)價(jià)

根據(jù)2016年～2020年金輪斷面監(jiān)測(cè)數(shù)據(jù)，通過(guò)相關(guān)系數(shù)Rs計(jì)算，金輪斷面總磷的Rs為-0.900，總磷濃度呈顯著下降趨勢(shì)，高錳酸鹽指數(shù)和氨氮的Rs分別為-0.100和-0.800，兩者濃度均呈不顯著下降趨勢(shì)，五日生化需氧量的Rs為0.100，五日生化需氧量濃度呈不顯著上升趨勢(shì)[3]。通過(guò)Pi和趨勢(shì)因子 計(jì)算得出，金輪斷面總磷項(xiàng)目的Pi大于1、大于等于0，總磷是現(xiàn)實(shí)污染因子，但污染程度將趨于減輕；高錳酸鹽指數(shù)、五日生化需氧量和氨氮項(xiàng)目Pi小于等于1、大于等于0，高錳酸鹽指數(shù)、五日生化需氧量和氨氮濃度將呈下降趨勢(shì)，是非污染因子，維持現(xiàn)狀即可，見(jiàn)表1。由此可見(jiàn)，2016年～2020年期間該斷面污染治理措施基本有效，但總磷還處于超標(biāo)狀態(tài)，且常年處于劣Ⅴ，亟需加強(qiáng)治理。

表1 金輪斷面水質(zhì)環(huán)境結(jié)果評(píng)價(jià)和趨勢(shì)分析 單位：mg/L

表1中：①Spearman秩相關(guān)系數(shù)臨界值Wp0.05為0.900。②Pi=Ci/Co， Pi大于1，小于0，表明污染已經(jīng)發(fā)生，且程度正在加劇，是主要污染因子，是環(huán)境治理的關(guān)鍵污染因子；Pi大于1，大于等于0，表明污染已經(jīng)發(fā)生，程度趨于減輕，是環(huán)境治理的現(xiàn)實(shí)污染因子；Pi小于等于1，小于0，表明污染程度正增加，未超標(biāo)，是環(huán)境污染的潛在因子；Pi小于等于1，大于等于0，表示該因子未超標(biāo)，且在減少，是非污染因子。 ③當(dāng) 大于0，污染濃度隨時(shí)間推移降低，值越大，降低越明顯；小于0，污染濃度隨時(shí)間推移變大，值越小，變大越明顯；等于0，污染濃度隨時(shí)間推移無(wú)變化。

2021年，金輪斷面高錳酸鹽指數(shù)、五日生化需氧量、氨氮和總磷年均濃度分別為2.1 mg/L、1.8 mg/L、0.29 mg/L和0.16 mg/L，均達(dá)到地表水Ⅲ類水質(zhì)標(biāo)準(zhǔn)，尤其是總磷，水質(zhì)質(zhì)量由2020 年劣Ⅴ變?yōu)棰箢?，但個(gè)別時(shí)段仍存在超標(biāo)現(xiàn)象。

為研究其它監(jiān)測(cè)項(xiàng)目對(duì)總磷的影響，選取2021 年金輪斷面監(jiān)測(cè)項(xiàng)目水溫、pH、溶解氧、流量、高錳酸鹽指數(shù)、氨氮、化學(xué)需氧量、五日生化需氧量、石油類和總磷監(jiān)測(cè)數(shù)據(jù)進(jìn)行變量間相關(guān)性分析，通過(guò)分析兩個(gè)或多個(gè)具備相關(guān)性的變量元素，以此衡量變量因素的相關(guān)密切程度。通過(guò)表2數(shù)據(jù)可以看出，總磷與化學(xué)需氧量、五日生化需氧量有低度正相關(guān)，總磷與高錳酸鹽指數(shù)、氨氮相關(guān)系數(shù)分別為0.617和0.743，具有顯著正相關(guān)，說(shuō)明高錳酸鹽指數(shù)和氨氮濃度的消減對(duì)總磷濃度的降低是有益的?？偭诐舛群土髁恐g的相關(guān)系數(shù)為-0.004，兩者之間不存在相關(guān)，反映出金輪斷面總磷濃度不受流量變化影響。

表2 金輪斷面監(jiān)測(cè)項(xiàng)目相關(guān)系數(shù)矩陣

3 “十三五”金輪斷面水質(zhì)污染的主要因素

3.1 本底值超標(biāo)問(wèn)題不容忽視

德陽(yáng)市龍門(mén)山脈屬磷礦富集地，自20世紀(jì)50 年代末就已開(kāi)始磷化工產(chǎn)業(yè)的生產(chǎn)，早先的粗放式發(fā)展和磷礦的大面積開(kāi)采，對(duì)流域總磷本底值有較高貢獻(xiàn)。

3.2 收納上游工業(yè)廢水

德陽(yáng)市是我國(guó)重要的磷化工基地，石亭江流域沿岸分布著眾多磷化工企業(yè)和工業(yè)園區(qū)，主要集中在什邡市和綿竹市，大量含磷廢水排入石亭江，對(duì)石亭江總磷貢獻(xiàn)頗大。而金輪斷面處于涉磷企業(yè)和工業(yè)園區(qū)的下游，收納了上游含磷工業(yè)廢水。

3.3 沿岸污水處理設(shè)施運(yùn)行效率較低

由于多數(shù)鄉(xiāng)鎮(zhèn)污水處理設(shè)施未正常運(yùn)行，部分污水處理設(shè)施設(shè)計(jì)處理能力不能滿足日益增長(zhǎng)的生活污水處理需求，農(nóng)村生活污水處理效率低，對(duì)石亭江金輪斷面污染日益加重。

4 水污染治理措施

通過(guò)近幾年的監(jiān)測(cè)數(shù)據(jù)分析，總磷作為金輪斷面的現(xiàn)實(shí)污染因子，高錳酸鹽指數(shù)和氨氮對(duì)總磷有顯著正相關(guān)，且總磷濃度不受流量的影響。所以，若想改善金輪斷面總磷污染現(xiàn)象，即要從源頭隔絕總磷的來(lái)源，又要同時(shí)削減高錳酸鹽指數(shù)和氨氮的排放。

4.1 以突出問(wèn)題整改為契機(jī)，提高水環(huán)境質(zhì)量

通過(guò)中央、省級(jí)環(huán)保督察，以及長(zhǎng)江經(jīng)濟(jì)帶生態(tài)環(huán)境問(wèn)題反饋和自查，對(duì)石亭江流域兩岸，尤其是石亭江什邡綿竹段控制單元存在的問(wèn)題，健全責(zé)任體系，加強(qiáng)調(diào)度督辦和部門(mén)聯(lián)動(dòng)協(xié)調(diào)，逐一明確整改任務(wù)、整改目標(biāo)、整改責(zé)任和整改時(shí)限，以整改為契機(jī)，加強(qiáng)流域周邊企業(yè)管理，開(kāi)展“雙隨機(jī)”“測(cè)管協(xié)同”監(jiān)管模式，嚴(yán)厲打擊環(huán)境違法行為，不斷改善石亭江流域水環(huán)境質(zhì)量。

4.2 開(kāi)展“三磷”專項(xiàng)整治

開(kāi)展“三磷”污染專項(xiàng)整治，完成礦井井硐封閉工作，拆除地面生產(chǎn)生活設(shè)施。加強(qiáng)石亭江流域兩岸磷化工企業(yè)廢水整治，抓好冷凝水、初期雨水、生產(chǎn)廢水“3”回用，進(jìn)行有效循環(huán)使用不外排；做好雨污、清污“2”分流；做好水平衡，實(shí)現(xiàn)全廠水平衡。2021 年，所有涉磷企業(yè)廢水不再排入石亭江流域，涉磷企業(yè)“地涌水”全部處理回收再利用，不再外排至石亭江流域。

同時(shí)，加大磷化工企業(yè)磷石膏堆場(chǎng)規(guī)范化整治和綜合利用進(jìn)度，探采礦權(quán)依法有序退出，提高磷石膏利用率，實(shí)現(xiàn)了“消增削存”，很大程度上削減了石亭江兩岸的總磷入河負(fù)荷。

4.3 強(qiáng)化污水處理設(shè)施建設(shè)及運(yùn)維

在全市推行“全域覆蓋、規(guī)劃統(tǒng)籌、投建營(yíng)運(yùn)一體、收集處理一體、政府購(gòu)買(mǎi)服務(wù)”的污水治理新模式。沿岸城市不斷加大污水處理設(shè)施的建設(shè)，完善鄉(xiāng)鎮(zhèn)污水管網(wǎng)建設(shè)，針對(duì)問(wèn)題管網(wǎng)制定整改方案，全面開(kāi)展鄉(xiāng)鎮(zhèn)污水管理排查補(bǔ)漏，提高污水收集率，做到污水不溢流、污水達(dá)標(biāo)排放，2021 年實(shí)施管網(wǎng)修復(fù)改造和新建項(xiàng)目共10 余個(gè)。污水處理效率的提高，大大削減了化學(xué)需氧量和氨氮等污染物的排放，顯著改善石亭江流域水質(zhì)。

4.4 沿岸工業(yè)園區(qū)升級(jí)改造

金輪斷面上游沿岸分布著雙盛和新市兩大工業(yè)園區(qū)，針對(duì)工業(yè)園區(qū)污水排放情況，對(duì)污水處理設(shè)施進(jìn)行升級(jí)改造，安裝除磷設(shè)施、二級(jí)生化池，以及在總排口建設(shè)濕地。同時(shí)，安裝在線監(jiān)測(cè)系統(tǒng)，加強(qiáng)監(jiān)管力度。

5 結(jié)語(yǔ)

本文通過(guò)Spearman秩相關(guān)系數(shù)法和灰色系統(tǒng)模型GM（1，1）趨勢(shì)因子a法，對(duì)石亭江流域金輪斷面的水質(zhì)進(jìn)行的分析和討論，并對(duì)各種污染因子進(jìn)行分類，數(shù)據(jù)分析發(fā)現(xiàn)，近幾年金輪斷面的水質(zhì)質(zhì)量趨于變好，但總磷項(xiàng)目存在偶爾超標(biāo)現(xiàn)象，是環(huán)境治理的現(xiàn)實(shí)污染因子。通過(guò)元素相關(guān)性分析，找出對(duì)總磷有影響的監(jiān)測(cè)項(xiàng)目，為石亭江金輪斷面精準(zhǔn)治污提供了依據(jù)。通過(guò)監(jiān)測(cè)數(shù)據(jù)的分析，結(jié)合石亭江流域周?chē)鷮?shí)際環(huán)境條件，分析了金輪斷面污染的原因，探討了石亭江金輪斷面水質(zhì)改善的治理措施，可為其它類似流域污染治理提供參考。