范燦鵬，熊杰

(珠海市西部生態(tài)環(huán)境監(jiān)測中心，廣東 珠海 519040)

金灣區(qū)是珠海市3大行政區(qū)之一，位于珠海市西南部，地處珠江出?？谀サ堕T與崖門之間的南海之濱。2016—2020年，金灣區(qū)直接管轄三灶、紅旗2個鎮(zhèn)，2019年常住人口19.46萬人。金灣區(qū)現(xiàn)有主要集中式飲用水水源地包括三灶鎮(zhèn)木頭沖水庫(以下簡稱“木頭沖水庫”)、黃綠背水庫和紅旗鎮(zhèn)愛國水庫(以下簡稱“愛國水庫”)，集中式飲用水水源地水質(zhì)的優(yōu)劣直接影響著人體健康，而水質(zhì)監(jiān)測是保證飲用水安全的重要手段[1-2]?，F(xiàn)根據(jù)這3個水庫2016—2020年的水質(zhì)監(jiān)測數(shù)據(jù)，對其水質(zhì)狀況及變化趨勢進(jìn)行分析和評價(jià)，為有針對性地保護(hù)飲用水水源，進(jìn)一步改善水質(zhì)，優(yōu)化水源管理提供科學(xué)依據(jù)。

1 研究方法

1.1 研究區(qū)概況

金灣區(qū)中心位置位于東經(jīng)113°05′—113°25′、北緯21°59′—22°25′，毗鄰香港、澳門，南連大西國際水道，北靠珠江三角洲腹地，轄區(qū)總面積447.6 km2。境內(nèi)屬多水地帶，有黃楊河、泥灣門、磨刀門、崖門等水系分布。

1.2 數(shù)據(jù)來源

2016—2020年珠海市金灣區(qū)環(huán)境監(jiān)測站根據(jù)《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002)中31項(xiàng)指標(biāo)的歷史監(jiān)測數(shù)據(jù)[3]。

1.3 監(jiān)測點(diǎn)位、頻次

木頭沖水庫、黃綠背水庫和愛國水庫每個水庫各布設(shè)一個監(jiān)測斷面，每月監(jiān)測1次。

1.4 評價(jià)方法

采用單因子指數(shù)評價(jià)法、綜合污染指數(shù)評價(jià)法和綜合營養(yǎng)狀態(tài)指數(shù)法。

1.4.1 單因子評價(jià)法

即最差因子判別法，根據(jù)實(shí)測得到的數(shù)據(jù)，按《GB 3838—2002》的標(biāo)準(zhǔn)限值進(jìn)行對比分類，根據(jù)評價(jià)時(shí)段內(nèi)該斷面參評的指標(biāo)中類別最高的一項(xiàng)來確定類別，參與評價(jià)的指標(biāo)為該標(biāo)準(zhǔn)中除水溫、總氮、糞大腸菌群以外的 21 項(xiàng)指標(biāo)[4]。

1.4.2 綜合污染指數(shù)評價(jià)法

使用均值綜合污染指數(shù)法[5-6]，對各污染指標(biāo)的相對污染指數(shù)進(jìn)行統(tǒng)計(jì)，得出代表水體綜合污染程度的數(shù)值。評價(jià)指標(biāo)同單因子評價(jià)法，根據(jù)歷年監(jiān)測結(jié)果選取溶解氧、高錳酸鹽指數(shù)、化學(xué)需氧量、五日生化需氧量、氨氮、總磷、銅、鋅、砷、氟化物、揮發(fā)酚、石油類、陰離子表面活性劑共13個主要污染因子作為評價(jià)指標(biāo)，其他項(xiàng)目常年未檢出或偶有檢出濃度很低，故不參與統(tǒng)計(jì)分析，各項(xiàng)目的評價(jià)標(biāo)準(zhǔn)限值為《GB 3838—2002》[3]中Ⅲ類水標(biāo)準(zhǔn)限值。綜合污染指數(shù)公式如下：

(1)

(2)

式中：P——均值綜合污染指數(shù)；pi——單因子污染指數(shù)；ci——第i項(xiàng)指標(biāo)的實(shí)際濃度測量值；c0——第i項(xiàng)指標(biāo)的評價(jià)標(biāo)準(zhǔn)限值；n=13。

1.4.3 綜合營養(yǎng)狀態(tài)指數(shù)法

根據(jù)《地表水環(huán)境質(zhì)量評價(jià)辦法(試行)》規(guī)定，采用綜合營養(yǎng)狀態(tài)指數(shù)法[TLI(∑)]，以葉綠素a(Chl-a)、總氮(TN)、總磷(TP)、透明度(SD)、高錳酸鹽指數(shù)(CODMn)5個項(xiàng)目為指標(biāo)，對3個飲用水水源地富營養(yǎng)化水平進(jìn)行評價(jià)。綜合營養(yǎng)狀態(tài)指數(shù)計(jì)算公式如下[4,7]：

(3)

式中：TLI(∑)——綜合營養(yǎng)狀態(tài)指數(shù)；Wj——第j種參數(shù)的營養(yǎng)狀態(tài)指數(shù)的相關(guān)權(quán)重；TLI(j)——為第j種參數(shù)的營養(yǎng)狀態(tài)指數(shù)。

以Chl-a作為基準(zhǔn)參數(shù)，則第j種參數(shù)的歸一化的相關(guān)權(quán)重計(jì)算公式為：

(4)

式中：rij——第j種參數(shù)與基準(zhǔn)參數(shù)Chl-a的相關(guān)系數(shù)；m——評價(jià)參數(shù)的個數(shù)。

表1 中國湖泊(水庫)部分參數(shù)與Chl-a的相關(guān)關(guān)系和Wj值

各項(xiàng)目營養(yǎng)狀態(tài)指數(shù)計(jì)算公式見(5)—(9)：

TLI(Chl-a)=10[2.5+1.086ln(chl-a)]

(5)

TLI(TP)=10[9.436+1.624ln(TP)]

(6)

TLI(TN)=10[5.453+1.694ln(TN)]

(7)

TLI(SD)=10[5.118-1.94ln(SD)]

(8)

TLI(CODMn)=10[0.109+2.661lnρ(CODMn)]

(9)

式中：Chl-a單位為mg/m3，SD單位為m；其他指標(biāo)單位均為mg/L。采用 0～100 的一系列連續(xù)數(shù)字對湖泊(水庫)營養(yǎng)狀態(tài)進(jìn)行分級：TLI(∑)<30 ，為貧營養(yǎng)；30≤TLI(∑)≤50 ，為中營養(yǎng)；TLI(∑)>50，為富營養(yǎng)；5070，為重度富營養(yǎng)。

1.5 數(shù)據(jù)統(tǒng)計(jì)分析

采用EXCEL和 SPSS 20.0 軟件對數(shù)據(jù)進(jìn)行處理。水庫綜合污染指數(shù)年度變化趨勢分析采用 Spearman 秩相關(guān)系數(shù)法進(jìn)行變化趨勢分析，當(dāng)N=5時(shí)，顯著性水平為 0.05的臨界值Wp= 0.900，統(tǒng)計(jì)檢驗(yàn)用的秩相關(guān)系數(shù)計(jì)算公式如下[4,8]：

(10)

di=Xi-Yi

(11)

式中：di——變量Xi與Yi的差值；Xi——周期I到周期N按濃度值從小到大排列的序號；Yi——按時(shí)間排列的序號。若rs為負(fù)值，則表明在評價(jià)時(shí)段內(nèi)有關(guān)統(tǒng)計(jì)量指標(biāo)變化呈下降趨勢或好轉(zhuǎn)趨勢；若rs為正值，則表明在評價(jià)時(shí)段內(nèi)有關(guān)統(tǒng)計(jì)量指標(biāo)變化呈上升趨勢或加重趨勢；當(dāng)rs≤Wp則表明變化趨勢沒有顯著意義，說明在評價(jià)時(shí)段內(nèi)水質(zhì)變化穩(wěn)定或平穩(wěn)。

2 結(jié)果與討論

2.1 單因子評價(jià)結(jié)果

2016—2020年金灣區(qū)飲用水水源主要指標(biāo)監(jiān)測結(jié)果見表1。表1選取的是影響水質(zhì)類別的主要監(jiān)測指標(biāo)。依據(jù)各監(jiān)測斷面 21 項(xiàng)指標(biāo)年均值統(tǒng)計(jì)數(shù)據(jù)發(fā)現(xiàn)，2016—2020年金灣區(qū)水質(zhì)比較穩(wěn)定。由表1可見，除2016年黃綠背水庫超出地表水Ⅲ類標(biāo)準(zhǔn)外，超標(biāo)因子為總磷，其他年份和其他水庫水質(zhì)類別均維持在Ⅱ—Ⅲ類，達(dá)到或優(yōu)于地表水Ⅲ類斷面占比93.3%。

表2 2016—2020年金灣區(qū)飲用水水源主要指標(biāo)監(jiān)測結(jié)果 mg/L

金灣區(qū)飲用水水源2016—2020年水質(zhì)類別見圖1。由圖1可見，2016年3個水庫超出地表水Ⅲ類標(biāo)準(zhǔn)的比例最高，其中黃綠背水庫超出地表水Ⅲ類標(biāo)準(zhǔn)占50%，水質(zhì)最差。2016年影響黃綠背水庫和愛國水庫的主要污染指標(biāo)為總磷，木頭沖水庫超標(biāo)因子有化學(xué)需氧量、溶解氧和總磷。2018年4月份木頭沖水庫超出Ⅲ類標(biāo)準(zhǔn)，主要污染指標(biāo)為化學(xué)需氧量、五日生化需氧量和總磷。針對2016年水庫達(dá)標(biāo)率低的情況，已從點(diǎn)源、面源和水土流失等方面著手，加強(qiáng)入庫面源污染的控制，通過設(shè)置集中式飲用水水源地保護(hù)區(qū)保護(hù)范圍圍網(wǎng)和界樁、界碑等警示標(biāo)志，植樹綠化生態(tài)護(hù)土，加強(qiáng)人工巡查、蓄水調(diào)節(jié)等工作，加強(qiáng)對飲用水水源保護(hù)區(qū)的監(jiān)管，確保水體質(zhì)量總體穩(wěn)定。從后期實(shí)際監(jiān)測結(jié)果看，水質(zhì)保持穩(wěn)定性較好，總磷超標(biāo)情況有所改觀。

圖1 2016—2020年金灣區(qū)飲用水水源水質(zhì)類別

2.2 水質(zhì)綜合污染指數(shù)評價(jià)結(jié)果

2016—2020年金灣區(qū)3個水庫水質(zhì)綜合污染指數(shù)趨勢變化見圖2。由圖2可見，3個水庫的水質(zhì)綜合指數(shù)均在2016年達(dá)到最高，水質(zhì)最差，這與單因子評價(jià)結(jié)果一致。其中，黃綠背水庫的綜合指數(shù)范圍變化為0.16～0.32，在2017年達(dá)到最低值，水質(zhì)最好；愛國水庫為0.16～0.30，在2017年水質(zhì)最好；木頭沖水庫為0.22～0.32，2017年和2019年水質(zhì)最好。為定量分析3個水庫水污染變化趨勢，對水質(zhì)綜合污染指數(shù)的年際變化趨勢做秩次相關(guān)檢驗(yàn)，結(jié)果表明，黃綠背水庫和愛國水庫秩相關(guān)系數(shù)均為-0.55，木頭沖水庫秩相關(guān)系數(shù)為-0.50，取置信度為95% 的Wp=0.9 作為判斷標(biāo)準(zhǔn)，表明在2016—2020年期間3個水源地水質(zhì)總體呈好轉(zhuǎn)趨勢，但改善趨勢不顯著，水質(zhì)變化平穩(wěn)。

圖2 2016—2020年各水庫水質(zhì)綜合污染指數(shù)變化

2016—2020年飲用水水源污染指數(shù)見表3。由表3可見，影響水質(zhì)的主要因素為溶解氧、總磷、化學(xué)需氧量和五日生化需氧量，這4項(xiàng)指標(biāo)分擔(dān)率之和每年均超過60%；從單因子污染指數(shù)可看出，除了2016年黃綠背水庫總磷的污染指數(shù)>1，超出Ⅲ類水標(biāo)準(zhǔn)外，其他指標(biāo)均達(dá)到Ⅲ類水要求。

2.3 綜合營養(yǎng)狀態(tài)指數(shù)評價(jià)結(jié)果

2016年3個水庫均未對葉綠素a指標(biāo)進(jìn)行監(jiān)測，因此該年度的綜合營養(yǎng)狀態(tài)指數(shù)無法計(jì)算統(tǒng)計(jì)。2017—2020年金灣區(qū)飲用水水源營養(yǎng)狀態(tài)指數(shù)見表4。由表4可見，2017—2020年3個水庫均處于貧營養(yǎng)狀態(tài)，未呈現(xiàn)富營養(yǎng)化。黃綠背水庫TLI(∑)在15～17范圍內(nèi)，無明顯差異；愛國水庫TLI(∑)在10～19范圍內(nèi)，且呈上升趨勢，主要是TLI(SD)和TLI(Chl-a)指數(shù)逐年升高引起的綜合營養(yǎng)狀態(tài)指數(shù)變大。木頭沖水庫TLI(∑)在16～24范圍內(nèi)波動，在2018年值最高。2017年和2018年3個水庫營養(yǎng)狀態(tài)優(yōu)劣排序均為愛國水庫>黃綠背水庫>木頭沖水庫。2019年3個水庫綜合營養(yǎng)狀態(tài)指數(shù)基本持平，2020年?duì)I養(yǎng)狀態(tài)優(yōu)劣排序?yàn)辄S綠背水庫>愛國水庫>木頭沖水庫。

表4 2017—2020年飲用水水源營養(yǎng)狀態(tài)指數(shù)

3 結(jié)論

(1)單因子評價(jià)結(jié)果表明，除2016年黃綠背水庫超出地表水Ⅲ類標(biāo)準(zhǔn)外，超標(biāo)因子為總磷，其他年份和其他水庫水質(zhì)類別均維持在Ⅱ—Ⅲ類，2016年度3個水庫超出地表水Ⅲ類標(biāo)準(zhǔn)的比例最高。

(2)3個水庫的水質(zhì)綜合污染指數(shù)均在2016年達(dá)到最高，水質(zhì)最差，2016—2020年3個水源地水質(zhì)總體呈好轉(zhuǎn)趨勢，但改善趨勢不顯著。影響水質(zhì)的主要因素為溶解氧、總磷、化學(xué)需氧量和五日生化需氧量。

(3)2017—2020年期間3個水庫均處于貧營養(yǎng)狀態(tài)，未呈現(xiàn)富營養(yǎng)化。