劉乾甫，杜浩，賴子尼，帥方敏，郭霞

(1.農(nóng)業(yè)農(nóng)村部淡水生物多樣性保護(hù)重點(diǎn)實(shí)驗(yàn)室，中國(guó)水產(chǎn)科學(xué)研究院長(zhǎng)江水產(chǎn)研究所，武漢 430223;2.中國(guó)水產(chǎn)科學(xué)研究院珠江水產(chǎn)研究所，廣州 510380)

珠江是中國(guó)南方最大的水系，干流全長(zhǎng)2 214 km，為中國(guó)第3大河流，其徑流量?jī)H次于長(zhǎng)江，居全國(guó)第2位。珠江流域地勢(shì)大體呈西高東低、北高南低趨勢(shì)，上游處于云貴高原與黔桂地區(qū)的高山峽谷，平均海拔超過(guò)1 000米，中游為山地與丘陵相間，下游為珠江三角洲沖擊平原。西江是珠江的主干流，位于珠江干流的中上游地區(qū)，自西向東流經(jīng)中國(guó)云南、貴州、廣西和廣東四省(自治區(qū))，其主要干支流有南盤江、北盤江、紅水河、柳江、郁江、桂江、黔江、潯江及賀江等[1]。珠江中上游屬于熱帶——亞熱帶季風(fēng)氣候區(qū)，常年溫暖，陽(yáng)光充足，雨量豐沛，加上流域內(nèi)山地較多，河道狹窄，水流湍急，水體自然凈化能力較強(qiáng)，水生態(tài)環(huán)境狀況良好，在歷史上一直是沿江城市居民日常用水的主要來(lái)源。珠江中上游江段優(yōu)質(zhì)的水文地理環(huán)境也孕育了珠江流域眾多珍稀特有經(jīng)濟(jì)魚類資源。早期調(diào)查表明珠江中上游江段產(chǎn)卵場(chǎng)密集，魚類群落多樣性豐富，漁業(yè)資源量較高，其特殊的水文地理、礁群流體環(huán)境為魚類提供了良好的產(chǎn)卵、孵化及育肥等繁育場(chǎng)所[2]。然而，近年來(lái)，隨著河流的過(guò)度開(kāi)發(fā)、珠江沿岸城鎮(zhèn)化的發(fā)展，以及流域內(nèi)農(nóng)業(yè)生產(chǎn)中投放的大量農(nóng)藥化肥原料被雨水帶入河流，西江漁業(yè)資源與環(huán)境受到了巨大的影響。多次調(diào)查數(shù)據(jù)表明，由于受水壩阻隔等多種因素影響，珠江漁業(yè)資源已呈現(xiàn)衰退趨勢(shì)[3-4]。前期關(guān)于珠江水系的環(huán)境調(diào)查研究主要集中于珠江下游的珠三角河網(wǎng)和珠江口區(qū)域[5-6]，對(duì)于珠江中上游較多見(jiàn)的是從魚類和浮游生物群落角度進(jìn)行漁業(yè)生物資源的調(diào)查研究，間接反應(yīng)出環(huán)境的變化及其影響[7-8]，而對(duì)于其水體環(huán)境狀況，尤其是水體主要污染因子及其污染程度均缺乏相應(yīng)的監(jiān)測(cè)分析。

為此，本研究分別于2016年3月(平水期)、6月(豐水期)和9月(枯水前期)對(duì)廣東省肇慶市以上的珠江中上游江段的水體進(jìn)行了調(diào)查取樣與測(cè)定分析，并對(duì)測(cè)定結(jié)果進(jìn)行了水環(huán)境狀況的分析與評(píng)價(jià)，以期全面了解珠江中上游水環(huán)境現(xiàn)狀、水體的污染狀況和主要污染因素，并基于本研究的水體指標(biāo)測(cè)定與分析結(jié)果，為珠江中上游水體環(huán)境保護(hù)與管理、珠江流域漁業(yè)資源保護(hù)與修復(fù)提供數(shù)據(jù)支撐和理論基礎(chǔ)。

1 材料與方法

1.1 采樣站位概況

于2016年3月(平水期)、6月(豐水期)、9月(枯水前期)，分別開(kāi)展對(duì)珠江中上游魯布革、八渡、白層、大化、合山、石龍、武宣、桂平(大騰峽)、石咀(東塔產(chǎn)卵場(chǎng))、平南、藤縣、壩上(梧州長(zhǎng)洲壩上)、梧州(梧州長(zhǎng)洲壩下)、封開(kāi)、郁南、德慶和肇慶等17個(gè)樣點(diǎn)的野外調(diào)查采樣工作。其中，魯布革采樣點(diǎn)處于珠江上游末端魯布革水電站修筑后形成的萬(wàn)豐湖區(qū)域，八渡、白層位于南、北盤江流域，大化、合山、石龍和武宣位于紅水河干流，其他樣點(diǎn)均位于西江干流(圖1)。

圖1 珠江中上游采樣點(diǎn)分布圖LBG, BD, BC, DH, HS, SL, WX, GP, SZ, PN, TX, BS, WZ, FK, YN, DQ及ZQ分別表示采樣點(diǎn)魯布革、八渡、白層、大化、合山、石龍、武宣、桂平、石咀、平南、藤縣、長(zhǎng)洲壩上、梧州、封開(kāi)、郁南、德慶和肇慶。下同。Fig.1 Sampling sites in the upper and middle reaches of the Pearl RiverLBG, BD, BC, DH, HS, SL, WX, GP, SZ, PN, TX, BS, WZ, FK, YN, DQ and ZQ are, represented for sampling stations of Lubuge,Badu, Baiceng, Dahua, Heshan, Shilong, Wuxuan, Guiping, Shizui, Pingnan, Tengxian, Changzhoubashang, Wuzhou, Fengkai,Yunan, Deqing and Zhaoqing, respectively. The same below.

1.2 樣品采集與分析

1.3 水體環(huán)境質(zhì)量分析與評(píng)價(jià)

1.3.1 評(píng)價(jià)指標(biāo)及其限定值

1.3.2 評(píng)價(jià)方法

依據(jù)文獻(xiàn)[16]，采用單項(xiàng)污染指數(shù)和負(fù)荷比對(duì)監(jiān)測(cè)參數(shù)進(jìn)行單項(xiàng)評(píng)價(jià)，采用綜合污染指數(shù)對(duì)水體環(huán)境質(zhì)量進(jìn)行整體評(píng)價(jià)。單項(xiàng)污染指數(shù)計(jì)算公式見(jiàn)式(1)。

Pi=Ci/Ci0

式(1)

式(1)中：Pi為水環(huán)境中污染物i的污染指數(shù)；Ci為水環(huán)境中污染物i的實(shí)測(cè)值；Ci0為水環(huán)境污染物i的限量標(biāo)準(zhǔn)值。Pi≤ 1時(shí)，表示水環(huán)境未受污染，指標(biāo)合格，Pi= 計(jì)算所得值；當(dāng)Pi>1時(shí)，表示水環(huán)境受到污染，指標(biāo)不合格，Pi=1.0 + 5 × lg(計(jì)算所得值)[17]。超標(biāo)率的計(jì)算公式見(jiàn)式(2)。

式(2)

在單項(xiàng)污染指數(shù)評(píng)價(jià)的基礎(chǔ)上，采用兼顧單項(xiàng)污染指數(shù)最大值和平均值的綜合污染指數(shù)Pj進(jìn)行評(píng)價(jià)，其計(jì)算公式見(jiàn)式(3)。

式(3)

式(3)中：Pimax為最大單項(xiàng)污染指數(shù)，Piave為平均單項(xiàng)污染指數(shù)。

污染物負(fù)荷比的計(jì)算公式見(jiàn)式(4)。

式(4)

式(4)中：Qi為污染物i的負(fù)荷比；n為參評(píng)指標(biāo)總數(shù)，Pi為水環(huán)境中污染物i的單項(xiàng)污染指數(shù)。

依據(jù)水體環(huán)境綜合污染指數(shù)，將水質(zhì)狀況分為5個(gè)等級(jí)[18](表1)，并對(duì)所研究的區(qū)域水體的污染程度和污染水平進(jìn)行評(píng)價(jià)。

表1 水質(zhì)狀況分級(jí)

Tab.1 Grades of water environment

等級(jí)Grade綜合污染指數(shù)Comprehensivepollutionindex污染程度Degreeofpollution污染水平Levelofpollution1Pj<0.7清潔清潔20.73.0重污染嚴(yán)重超警戒水平

1.3.3 統(tǒng)計(jì)分析方法

文章數(shù)據(jù)處理與作圖使用Excel 2007軟件，采用MATLAB軟件進(jìn)行聚類分析。

2 結(jié)果與分析

2.1 水體理化環(huán)境狀況

各采樣區(qū)域水體pH變化范圍為7.03～9.04，平水期、豐水期及枯水前期均值分別為8.11、7.83和7.86(圖2a)。水體溶解氧含量變化范圍為5.84～9.57 mg/L，3個(gè)時(shí)期DO均值變化不大，分別為7.48(平水期)、7.58(豐水期)和6.99 mg/L(枯水前期)(圖2b)。3個(gè)時(shí)期水體SD差別較大，豐水期最低，最小值僅為20 cm，均值為45 cm；其次為平水期，SD均值為127 cm；枯水前期SD最高，變化范圍為75～230 cm，均值達(dá)到149 cm(圖2c)。水體總磷含量呈現(xiàn)枯水前期明顯高于平水期和豐水期的分布狀況，均值分別為0.055 7(平水期)、0.095 4(豐水期)及0.173 6 mg/L(枯水前期)，但部分樣點(diǎn)在部分時(shí)期也會(huì)凸顯出較高的總磷水平，如2016年3月(平水期)梧州采樣點(diǎn)TP含量達(dá)到 0.267 3 mg/L(圖2d)。水體總氮含量變化范圍為0.636 0～3.715 3 mg/L，整個(gè)采樣期大部分采樣點(diǎn)總氮含量均高于1.0 mg/L，其中上游魯布革、八渡及白層等采樣區(qū)域的總氮含量高于中下游位點(diǎn)(圖2e)。水體氨氮含量在豐水期較平水期和枯水前期略高，平水期、豐水期及枯水前期均值分別為0.149 7、0.208 9和0.139 1 mg/L，其中豐水期氨氮含量峰值達(dá)到0.507 5 mg/L(桂平大藤峽)(圖2f)。水體Chla含量變化范圍為0.17 ～25.85 μg/L，其中處于魯布革水電站庫(kù)區(qū)內(nèi)水體Chla含量較高，平水期和枯水前期分別達(dá)到25.85 μg/L和17.72 μg/L，呈現(xiàn)富營(yíng)養(yǎng)化狀況，其余樣點(diǎn)在整個(gè)采樣期水體Chla含量大都在10 μg/L以內(nèi)，顯示“貧營(yíng)養(yǎng)—中營(yíng)養(yǎng)”水體特征(圖2g)。

圖2 珠江中上游主要水體理化指標(biāo)變化a:pH; b:溶解氧(DO); c:透明度(SD);d:總磷(TP); e:總氮(TN); f:氨氮 g:葉綠素a(Chla)。Fig.2 Changes of physical-chemical parameters in the upper and middle reaches of the Pearl Rivera:pH; b:dissolved oxygen; c:secchi depth; d:total phosphorus; e:total nitrogen; f:ammonia nitrogen; g:chlorophyll a.

表2 珠江中上游各水情期水體理化因子均值及范圍

Tab.2 Mean and range of physical-chemical parameters of each water regime in the upper and middlereaches of the Pearl Rivern=17

2.2 水質(zhì)狀況分析

表3 珠江中上游水質(zhì)狀況分析結(jié)果

Tab.3 Analyzing results of water quality in the upper and middle reaches of the Pearl Rivern=17

參數(shù)Parameters平水期Levelperiod豐水期Wetseason枯水前期Earlydryseason平均值A(chǔ)veragevalue超標(biāo)率Rateofultrastandard最大超標(biāo)倍數(shù)Multipleofmaxultrastandard平均值A(chǔ)veragevalue超標(biāo)率Rateofultrastandard最大超標(biāo)倍數(shù)Multipleofmaxultrastandard平均值A(chǔ)veragevalue超標(biāo)率Rateofultrastandard最大超標(biāo)倍數(shù)Multipleofmaxultrastandard限定值LimitspH8.11—/7.8311.76%0.067.86—/6.50～8.50DO(mg·L-1)7.48—/7.58—/6.99—/5.00NH3(mg·L-1)0.0117.65%0.440.0323.53%7.970.015.88%0.300.02NH+4-N(mg·L-1)0.15—/0.21—/0.14—/1.00NO-3-N(mg·L-1)1.7294.12%1.461.3994.12%1.861.8594.12%2.131.00NO-2-N(mg·L-1)0.01—/0.035.88%1.140.0811.76%0.640.15TN(mg·L-1)2.07100%1.852.05100%2.692.2694.12%2.721.00TP(mg·L-1)0.065.88%0.340.05—/0.1723.53%0.080.20SD(cm)127.00—/45.0082.35%0.67149.00 —/60.00Chla(μg·L-1)3.605.99%1.592.18—/4.285.88%0.7210.00

注：表中“—”表示沒(méi)有出現(xiàn)超標(biāo)；“/”表示沒(méi)有計(jì)算該數(shù)值。

2.3 水質(zhì)評(píng)價(jià)

從時(shí)間上看，平水期、豐水期和枯水前期珠江中上游水體綜合污染指數(shù)Pj分別為1.82、1.76和1.89，根據(jù)表1的漁業(yè)水質(zhì)分級(jí)，可以看到3個(gè)水情期水質(zhì)均為“輕污染”等級(jí)，污染水平為“輕度超警戒水平”。對(duì)每個(gè)采樣點(diǎn)的綜合污染指數(shù)Pj分析結(jié)果見(jiàn)圖3，由圖可以看到，上游魯布革、八渡、白層采樣站位水體受污染較為嚴(yán)重；中下游區(qū)域的石龍、武宣、藤縣和德慶等樣點(diǎn)在豐水期的水體綜合污染指數(shù)Pj明顯較另外兩個(gè)時(shí)期高。

基于水體綜合污染指數(shù)Pj，對(duì)珠江中上游各采樣點(diǎn)的水體環(huán)境質(zhì)量狀況進(jìn)行SOM(Self-organizing map)聚類分析，結(jié)果見(jiàn)圖4。珠江中上游水體環(huán)境大致可以分為3類：第I類，武宣、桂平、石咀、平南、壩上、封開(kāi)、郁南、德慶和肇慶等樣點(diǎn)，參照?qǐng)D3可知，這些樣點(diǎn)在3個(gè)水情期的水體綜合污染指數(shù)Pj及其均值都處于1～2之間，顯示整個(gè)調(diào)查期間水體環(huán)境質(zhì)量狀況比較穩(wěn)定，處于“輕污染”等級(jí)、“輕度超警戒水平”狀況；第II類，魯布革、八渡、白層、大化、合山和梧州等樣點(diǎn)，這些樣點(diǎn)在3個(gè)水情期的水體綜合污染指數(shù)Pj及其均值大都處于2～3之間，顯示整個(gè)調(diào)查期間水體環(huán)境質(zhì)量狀況比較嚴(yán)重，呈“中污染”等級(jí)、“中度超警戒水平”狀況；第III類，石龍和藤縣兩個(gè)樣點(diǎn)，二者在3個(gè)水情期的水體綜合污染指數(shù)Pj差別較大，處于1～5之間，豐水期水體環(huán)境質(zhì)量狀況比較嚴(yán)重，呈“重污染”等級(jí)、“嚴(yán)重超警戒水平”狀況，而在其他時(shí)期略低，處于1～2之間，“輕污染”等級(jí)、“輕度超警戒水平”狀況。

綜上，所調(diào)查的珠江中上游大部分采樣區(qū)域水體大都呈現(xiàn)“輕污染”、“輕度超警戒水平”狀況；(魯布革、八渡和白層)水質(zhì)受污染較嚴(yán)重，要劣于中下游水質(zhì)；但中下游部分樣點(diǎn)(石龍、藤縣和德慶)豐水期水體綜合污染指數(shù)偏高，水質(zhì)劣于其他時(shí)期(圖4)。

表4 珠江中上游水質(zhì)評(píng)價(jià)及污染等級(jí)

Tab.4 Assessment of major aquaculture pond water quality and pollution level in the upper and middle reachesof the Pearl Rivern=17

參數(shù)Parameters平水期Levelperiod豐水期Wetseason枯水前期Earlydryseason單項(xiàng)污染指數(shù)PiSinglepollutionindex負(fù)荷比/%QDutyration單項(xiàng)污染指數(shù)PiSinglepollutionindex負(fù)荷比/%QDutyration單項(xiàng)污染指數(shù)PiSinglepollutionindex負(fù)荷比/%QDutyrationpH0.587.290.636.580.333.67DO(mg·L-1)0.688.560.666.960.727.95NH3(mg·L-1)0.516.441.1912.510.566.20NH+4-N(mg·L-1)0.151.890.212.190.141.54NO-3-N(mg·L-1)2.1527.021.6417.192.2524.94NO-2-N(mg·L-1)0.070.930.252.580.586.43TN(mg·L-1)2.5732.332.4826.062.6729.66TP(mg·L-1)0.303.720.252.580.879.70SD(cm)0.556.922.0121.070.444.85Chla(μg·L-1)0.394.890.222.290.465.06綜合污染指數(shù)Pj1.821.761.89污染程度輕污染輕污染輕污染污染水平輕度超警戒水平輕度超警戒水平輕度超警戒水平

圖3 珠江中上游各采樣點(diǎn)的水體綜合污染指數(shù)Fig.3 Water comprehensive pollution index variation in the upper and middle reaches of the Pearl River

圖4 基于水體綜合污染指數(shù)的珠江中上游水體環(huán)境質(zhì)量狀況聚類分布圖Fig.4 Clustering graph of water quality based on water comprehensive pollution index in the upper and middle reaches of the Pearl River

3 討論

3.1 珠江中上游水體中的主要污染因子及影響分析

3.2 珠江中上游水體環(huán)境狀況的時(shí)空差異

3.3 珠江中上游水體污染程度與防控措施分析

致謝：本研究的野外調(diào)查取樣中，得到了李躍飛、楊計(jì)平、武智等同事的熱情幫助，并得到了仲愷農(nóng)業(yè)工程學(xué)院本科實(shí)習(xí)生黃志鵬同學(xué)在室內(nèi)樣品處理與分析方面的鼎力相助，在此一并致謝！