西藏拉魯濕地水質(zhì)綜合狀況時(shí)空變化特征及污染源分析

李梁婷 邵蓓 卓瑪曲西 姜貞貞 鄭琴 德珍 次巴拉松 邱豪

摘要 依據(jù)2019年1月（枯水期）、5月（平水期）和7月（豐水期）西藏拉魯濕地的采樣點(diǎn)數(shù)據(jù)，利用改進(jìn)后的內(nèi)梅羅污染指數(shù)評(píng)價(jià)方法，研究不同時(shí)期拉魯濕地水質(zhì)綜合狀況時(shí)空變化特點(diǎn)，分析大致的污染物來(lái)源。結(jié)果表明，濕地不同水體環(huán)境指標(biāo)濃度均表現(xiàn)為降低趨勢(shì)，其N(xiāo)H3-N、TP、COD和TN的平均值處于隨時(shí)間推移呈現(xiàn)降低的趨勢(shì)，而水體pH和DO的平均值表現(xiàn)為隨時(shí)間變化先升高再降低的趨勢(shì);濕地在平水期和豐水期的水質(zhì)情況比枯水期表現(xiàn)更為良好?？菟?，東北部表現(xiàn)為Ⅱ類(lèi)水質(zhì)，其余地區(qū)皆為Ⅰ類(lèi)水質(zhì)，其水質(zhì)綜合狀況區(qū)域分布特點(diǎn)為由東北部向西部逐漸變好;平水期和豐水期，整體水質(zhì)表現(xiàn)為Ⅰ類(lèi)水質(zhì)?？菟?，TN與NH3-N存在極顯著正相關(guān)（P<0.01）;平水期，TP與NH3-N存在顯著正相關(guān)（P<0.05），TN與NH3-N存在極顯著正相關(guān)（P<0.01）;豐水期，TN與NH3-N存在極顯著正相關(guān)（P<0.01）。從枯水期到豐水期，西藏拉魯濕地水體中DO、NH3-N、TP、COD和TN質(zhì)量濃度分別降低了61.05%、54.55%、50.00%、2.41%和50.00%，并且其水質(zhì)綜合狀況由枯水期到豐水期逐漸變好。拉魯濕地東北部北入水口水體質(zhì)量治理和附近的污染防治是下一步西藏拉魯濕地水體質(zhì)量治理和管控的重點(diǎn)。

關(guān)鍵詞 水質(zhì)綜合狀況;污染等級(jí);內(nèi)梅羅污染指數(shù);時(shí)空變化;污染源;西藏拉魯濕地

中圖分類(lèi)號(hào) X 52? 文獻(xiàn)標(biāo)識(shí)碼 A? 文章編號(hào) 0517-6611（2022）12-0067-08

doi：10.3969/j.issn.0517-6611.2022.12.018

開(kāi)放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Analysis of Spatial-temporal Variation Characteristics and Pollution Source of Water Quality Comprehensive Conditions in Lalu Wetland， Tibet

LI Liang-ting，SHAO Bei，Zhuomaquxi et al

（Central Laboratory of Geology and Mineral Exploration and Development Bureau of Tibet Autonomous Region， Lhasa， Tibet 850000）

Abstract According to the sampling data of Lalu Wetland in January 2019 （dry season）， May 2019 （normal season） and July 2019 （wet season）， the improved Nemerow pollution index evaluation method was used to study the temporal and spatial variation characteristics of the comprehensive water quality of Lalu wetland in different periods， and to analyze the approximate sources of pollutants. The results showed that the concentrations of environmental indicators in different water bodies of wetland showed a decreasing trend， and the average values of NH3-N， TP， COD and TN showed a decreasing trend with time， while the average values of pH and DO in water bodies showed a trend of increasing first and then decreasing with time. The water quality of wetlands in normal and wet periods was better than that in dry periods. During the dry season， the northeast showed Class II water quality， while the rest areas were Class I water quality， and the regional distribution characteristics of its comprehensive water quality were gradually getting better from northeast to west. During the normal season and wet season， the overall water quality was class I.During the dry season， there was a significant positive correlation between TN and NH3-N （P < 0.01）. During the normal season，there was a significant positive correlation between TP and NH3-N （P < 0.05），and there was a very significant positive correlation between TN and NH3-N （P < 0.01）. During the wet season，there was a significant positive correlation between TN and NH3-N during the wet season （P < 0.01）. From dry season to wet season， the mass concentrations of DO， NH3-N， TP， COD and TN decreased by 61.05%， 54.55%， 50.00%， 2.41% and 50.00%， respectively. Moreover， the comprehensive situation of its water quality gradually improved from dry season to wet season. Water quality control and pollution prevention in the north inlet of northeast Lalu Wetland are the key points of water quality control and control in Lalu Wetland of Tibet in the next step.

Key words Comprehensive conditions of water quality;Pollution level;Nemerow pollution index;Spatio-temporal variation;Pollution source;Lalu Wetland in Tibet

濕地是水陸相互作用形成的獨(dú)特的生態(tài)系統(tǒng)[1]，不但具有豐富的資源，還具有巨大的環(huán)境調(diào)節(jié)功能和生態(tài)、社會(huì)、經(jīng)濟(jì)效益[2]。近年來(lái)，隨著經(jīng)濟(jì)社會(huì)的快速發(fā)展和城鎮(zhèn)化進(jìn)程的加快，導(dǎo)致濕地面臨著區(qū)域生態(tài)環(huán)境破壞、水體功能弱化、面積不斷減少等問(wèn)題，從而制約了濕地的進(jìn)一步保護(hù)、開(kāi)發(fā)和利用[3]。因此，許多學(xué)者將濕地的水質(zhì)綜合狀況分析與評(píng)價(jià)列為重點(diǎn)研究課題，陳瑞紅等[4]對(duì)會(huì)仙巖溶濕地水質(zhì)進(jìn)行了綜合狀況監(jiān)測(cè)與分析，王佐霖等[5]對(duì)深圳壩光濕地水資源環(huán)境綜合狀況及其污染程度進(jìn)行了分析與評(píng)價(jià)，楊海江等[6]利用貝葉斯評(píng)價(jià)法對(duì)星海湖濕地水質(zhì)綜合狀況進(jìn)行分析與評(píng)價(jià)。由此可見(jiàn)，對(duì)濕地開(kāi)展水質(zhì)綜合狀況分析與評(píng)價(jià)也成為水資源環(huán)境研究熱點(diǎn)。針對(duì)濕地地表水常用的評(píng)價(jià)方法有單因子評(píng)價(jià)法[7]、綜合污染指數(shù)法[8]、模糊數(shù)學(xué)綜合評(píng)價(jià)法[9]和內(nèi)梅羅指數(shù)法[10]等。其中成熟的內(nèi)梅羅指數(shù)法是目前國(guó)內(nèi)應(yīng)用較為廣泛的水污染評(píng)價(jià)方法，該方法是一種統(tǒng)籌了極值或者突出最大值的計(jì)權(quán)型多因子環(huán)境質(zhì)量指數(shù)[11]，但因?yàn)樵撝笖?shù)法過(guò)于注重最大污染因子對(duì)水體質(zhì)量的影響[12]，進(jìn)而忽視了其他污染因子對(duì)于水體質(zhì)量的貢獻(xiàn)，而改進(jìn)后的內(nèi)梅羅指數(shù)法將污染因子的相關(guān)權(quán)重重新計(jì)入，以此來(lái)權(quán)衡不同水質(zhì)因子對(duì)整體水體質(zhì)量的影響，同時(shí)統(tǒng)籌考慮了不同評(píng)價(jià)指標(biāo)相互間的差異性與相對(duì)重要性[13]，以此使研究結(jié)果更具有客觀性。

西藏拉魯濕地國(guó)家級(jí)自然保護(hù)區(qū)是世界上海拔最高、面積最大的城市天然濕地。與拉魯濕地相關(guān)的河流主要有拉薩河、流沙河和堆龍河，相關(guān)的渠系主要包括從拉薩河引水的北干渠、中干渠和南干渠[14]。拉魯濕地在增加拉薩市空氣中氧氣含量和相對(duì)濕度，調(diào)蓄流沙河洪水和濕地區(qū)域地下水，維持拉薩市生態(tài)平衡等方面發(fā)揮巨大作用，同時(shí)具有科研教育等功能[15]。但隨著拉薩市高原開(kāi)放型城鎮(zhèn)化發(fā)展不斷推進(jìn)，由于拉魯濕地地處拉薩市城區(qū)，人口的增加、農(nóng)業(yè)和生活污水的排放，致使該濕地水環(huán)境質(zhì)量不容樂(lè)觀，并且路飛等[16]的研究也表明其面積日漸萎縮，可能對(duì)拉魯濕地的生態(tài)環(huán)境進(jìn)一步造成影響。由此許多學(xué)者從不同角度對(duì)拉魯濕地開(kāi)展了水體質(zhì)量研究，如張文駒等[17]采用標(biāo)準(zhǔn)指數(shù)法對(duì)拉魯濕地水體環(huán)境質(zhì)量情況進(jìn)行了分析與評(píng)價(jià);白永飛等[18]測(cè)定了拉魯濕地水體環(huán)境中不同重金屬含量并對(duì)其分布和來(lái)源進(jìn)行了分析;巴桑等[19]采用Pantle和Buck污生指數(shù)法對(duì)拉魯濕地水體環(huán)境質(zhì)量情況進(jìn)行了綜合評(píng)價(jià)。然而，以往的研究主要關(guān)注拉魯濕地某幾種指標(biāo)的水體環(huán)境狀況，從時(shí)空變化角度對(duì)拉魯濕地不同水體環(huán)境指標(biāo)開(kāi)展綜合評(píng)價(jià)的研究鮮見(jiàn)報(bào)道。筆者對(duì)2019年1月、5月和7月拉魯濕地水體pH、溶解氧（DO）、氨氮（NH3-N）、總磷（TP）、化學(xué)需氧量（COD）和總氮（TN）這6個(gè)指標(biāo)進(jìn)行分析，采用改進(jìn)后的內(nèi)梅羅指數(shù)法對(duì)拉魯濕地水質(zhì)綜合狀況進(jìn)行評(píng)價(jià)，研究其時(shí)空變化趨勢(shì)，并對(duì)污染源進(jìn)行分析，以期為拉魯濕地管理制定綜合保護(hù)和防治措施提供科學(xué)依據(jù)，同時(shí)將對(duì)西藏地區(qū)乃至整個(gè)青藏高原的濕地水體環(huán)境保護(hù)具有良好的示范作用。

1 資料與方法

1.1 研究區(qū)概況

西藏拉魯濕地國(guó)家級(jí)自然保護(hù)區(qū)（91°03′41″～91°06′48″E，29°39′25″～29°42′08″N）位于拉薩市區(qū)西北角，呈東西帶狀分布，海拔3 650 m，年降水主要集中在6—9月。濕地北面為高山;東北面與發(fā)源于拉薩市西北部山區(qū)的奪底和娘熱2條山溝支流匯合而成的流沙河相接;東面與拉薩市城關(guān)區(qū)拉魯鄉(xiāng)拉魯村居民生活區(qū)接壤;南面緊鄰拉薩市區(qū)，以中干渠和當(dāng)熱路為界;西面以北京西路與北面高山交接處的順通建材城為終點(diǎn)。同時(shí)拉魯濕地被通往西藏自治區(qū)黨校的魯定北路分為東西兩片濕地[18]，兩片濕地之間通過(guò)涵洞相連。拉魯濕地水源補(bǔ)給主要由中干渠與流沙河提供。中干渠補(bǔ)給水引自獻(xiàn)多水電站尾水[20]，引用流量為55.44 m3/s。流沙河補(bǔ)給水具有典型的山區(qū)河流特點(diǎn)，其流量變化存在明顯的枯平豐季現(xiàn)象，冬季基本無(wú)流量處于枯水期，春末處于平水期，夏季進(jìn)入汛期，降水較多處于豐水期，由此可見(jiàn)，流沙河是一條典型的季節(jié)性河流。綜上所述，拉魯濕地受季節(jié)性河流影響，其存在明顯的枯水期、平水期和豐水期。

1.2 數(shù)據(jù)來(lái)源及處理

結(jié)合拉魯濕地地形及水流流向特征，于2019年1月（枯水期）、5月（平水期）和7月（豐水期）在拉魯濕地北入水口（1號(hào)采樣點(diǎn)）、中干渠拉魯濕地進(jìn)水口（2號(hào)采樣點(diǎn)）、拉魯濕地觀景臺(tái)區(qū)域（3號(hào)采樣點(diǎn)）、拉魯濕地魯定北路出水口（4號(hào)采樣點(diǎn)）、拉魯濕地西區(qū)出水口（5號(hào)采樣點(diǎn)）、拉魯濕地西北門(mén)（6號(hào)采樣點(diǎn)）和拉魯濕地魯定北路西門(mén)（7號(hào)采樣點(diǎn)）共設(shè)置7個(gè)采樣點(diǎn)進(jìn)行水質(zhì)采樣（圖1）。利用采水器于各采樣點(diǎn)采集5～15 cm處水質(zhì)樣品1 L，現(xiàn)場(chǎng)測(cè)定pH，所有采樣點(diǎn)都取3組平行樣，并放置于聚乙烯瓶中帶回實(shí)驗(yàn)室，用于測(cè)定DO、NH3-N、TP、COD和TN。采用玻璃電極法測(cè)定水樣pH;采用典量法測(cè)定DO;采用納氏試劑比色法測(cè)定NH3-N;采用鉬酸銨分光光度法測(cè)定TP;采用重鉻酸鹽法測(cè)定COD;采用堿性過(guò)硫酸鉀消解紫外分光光度法測(cè)定TN。同時(shí)利用Excel 2007和SPSS 24對(duì)1月、5月和7月的水質(zhì)指標(biāo)進(jìn)行Pearson相關(guān)系數(shù)分析并用ArcGIS 10.7中克里金插值分析法（空間局部插值法）對(duì)不同單項(xiàng)水質(zhì)指標(biāo)和水質(zhì)綜合狀況評(píng)價(jià)結(jié)果進(jìn)行空間分析。

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

改進(jìn)后的內(nèi)梅羅指數(shù)法是基于傳統(tǒng)內(nèi)梅羅指數(shù)法上統(tǒng)籌考慮最大污染因子、不同因子平均指數(shù)與污染因子權(quán)重，使其水質(zhì)評(píng)價(jià)結(jié)果能夠?qū)⒆畲笪廴疽蜃优c權(quán)重相對(duì)較大的污染因子對(duì)于水質(zhì)的影響更好地反映出來(lái)。因此該研究采用改進(jìn)后的內(nèi)梅羅污染指數(shù)法對(duì)拉魯濕地不同污染因子（DO、NH3-N、TP、COD和TN）的權(quán)重進(jìn)行計(jì)算，同時(shí)綜合考慮不同評(píng)價(jià)指標(biāo)之間的差異性。

1.3.1 單項(xiàng)污染指數(shù)確定。

單項(xiàng)污染指數(shù)的計(jì)算公式如下[21]：

Fi=CiSi（1）

式（1）中，F(xiàn)i表示第i項(xiàng)水質(zhì)指標(biāo)的單項(xiàng)污染指數(shù);Ci表示第i項(xiàng)水質(zhì)指標(biāo)的實(shí)測(cè)值;Si表示第i項(xiàng)水質(zhì)指標(biāo)的標(biāo)準(zhǔn)值，按照西藏拉魯濕地整體水域功能定位要求，結(jié)合其為國(guó)家級(jí)自然保護(hù)區(qū)，此次以《地表水質(zhì)量標(biāo)準(zhǔn)》（GB 3838—2002）Ⅰ類(lèi)為計(jì)算標(biāo)準(zhǔn)。

1.3.2 污染因子權(quán)重計(jì)算。

按照不同水質(zhì)標(biāo)準(zhǔn)[22]，計(jì)算不同污染因子權(quán)重，此次選?、耦?lèi)水環(huán)境為標(biāo)準(zhǔn)，計(jì)算方法如下[11，21]：

Wi=SmaxSini=1SmaxSi（2）

式（2）中，Wi表示第i個(gè)污染因子權(quán)重;n表示污染因子個(gè)數(shù);Si表示第i個(gè)污染因子計(jì)算標(biāo)準(zhǔn)值;Smax表示n個(gè)污染因子計(jì)算標(biāo)準(zhǔn)值最大值。

1.3.3 改進(jìn)后的內(nèi)梅羅污染指數(shù)確定。改進(jìn)后的內(nèi)梅羅污染指數(shù)計(jì)算方法如下[21]：

P*=Fi，max+Fw22+（i）22（3）

式（3）中，P*表示第i個(gè)水質(zhì)樣品改進(jìn)后的內(nèi)梅羅污染指數(shù);Fi，max表示單項(xiàng)污染指數(shù)最大值;i表示單項(xiàng)污染指數(shù)平均值;Fw表示權(quán)重Wi最大的污染因子的F值。

2 結(jié)果與分析

2.1 拉魯濕地枯水期、平水期和豐水期區(qū)域水體監(jiān)測(cè)指標(biāo)變化特征

從拉魯濕地枯水期、平水期和豐水期的水體pH、DO、NH3-N、TP、COD和TN的變化特征（表1）可以看出，該濕地從枯水期到豐水期不同監(jiān)測(cè)指標(biāo)間變化較大，主要表現(xiàn)為pH平均值由8.27降低至8.24，DO平均值由6.11 mg/L降低至2.38 mg/L，NH3-N平均值由0.22 mg/L降低至0.10 mg/L，TP平均值由0.02 mg/L降低至0.01 mg/L，COD平均值由4.97 mg/L降低至4.85 mg/L，TN平均值由0.26 mg/L 降低至0.13 mg/L?？傮w而言，在枯水期到豐水期，拉魯濕地不同水體環(huán)境指標(biāo)濃度均表現(xiàn)為降低趨勢(shì)，其N(xiāo)H3-N、TP、COD和TN的平均值處于隨時(shí)間推移呈現(xiàn)降低的趨勢(shì)，其緣由為豐水期濕地整體水位上升，水體本身稀釋功能得到增強(qiáng)，故而致使污染物濃度降低，同時(shí)由于水生植物本身對(duì)水質(zhì)具有凈化能力，在豐水期水生植物處于大量繁殖階段，通過(guò)水生植物除去污染物的作用得到加強(qiáng)，也使豐水期污染物濃度低于枯水期，這與李紅海等[23]的研究結(jié)果基本一致;而水體pH和DO的平均值表現(xiàn)為隨時(shí)間變化先升高再降低的趨勢(shì)，其中DO降低幅度達(dá)到61.05%，DO降低幅度較大的原因?yàn)樨S水期濕地中的水生動(dòng)植物處于大量繁殖階段，需要損耗較多氧氣，同時(shí)由于7月和8月拉薩市處于主汛期，天氣狀況不是很好，多為陰雨天，水生植物光合作用較少導(dǎo)致產(chǎn)生氧氣不足，致使水體DO濃度降低，該研究結(jié)果與陳瑞紅等[4]的研究結(jié)果類(lèi)似，即豐水期濕地中水生動(dòng)植物處于繁殖階段，需要損耗大量氧氣，同時(shí)由于夏季陰雨天較多，水生植物光合作用較少，致使水體DO濃度降低。

根據(jù)圖2～7可以得出，枯水期濕地東北部水體pH較大，平水期濕地東北部、南部中段和西部水體pH較大，豐水期濕地東南部和北部中段水體pH較大?？菟诤推剿跐竦爻鞑緿O濃度處于較低水平外，其余DO濃度處于較高水平;豐水期濕地東北部、北部中段DO濃度較高。枯水期濕地東北部和西南部NH3-N濃度較高，平水期濕地東北部NH3-N濃度較高，豐水期濕地西南部NH3-N濃度較高。枯水期濕地西南部TP濃度較高，平水期和豐水期濕地東北部TP濃度較高?？菟凇⑵剿诤拓S水期濕地西部COD濃度都處于較高水平?？菟跐竦貣|北部和西南部TN濃度較高，平水期濕地東北部TN濃度較高，豐水期濕地西南部TN濃度較高。整體而言，枯水期濕地污染主要集中在東北部和西南部，污染原因可能由于東北部靠近城市居民生活區(qū)，該區(qū)域?yàn)榱魃澈尤肟?、西南部靠近道路交通十字路口且水位處于較低水平。平水期和豐水期濕地較高濃度污染所占面積較小且分散分布在東部和中部，其緣由為豐水期濕地整體水位上升，水體本身稀釋功能得到增強(qiáng)，故而致使污染物濃度降低，該研究結(jié)果與曾凱等[24]的研究結(jié)果類(lèi)似。從枯水期到豐水期拉魯濕地不同水體環(huán)境指標(biāo)濃度均表現(xiàn)為降低趨勢(shì)。

2.2 拉魯濕地不同時(shí)期區(qū)域水質(zhì)綜合狀況變化特征

西藏拉魯濕地水體質(zhì)量按照 Ⅰ 類(lèi)水環(huán)境為標(biāo)準(zhǔn)不同污染因子權(quán)重見(jiàn)表2，可以得出污染因子TP的權(quán)重最高，由此可以判定TP污染貢獻(xiàn)率最大。

按照改進(jìn)的內(nèi)梅羅污染指數(shù)法進(jìn)行計(jì)算，參照李蘇等[25]和關(guān)云鵬[26]對(duì)各類(lèi)別污染等級(jí)的表達(dá)，可以得出西藏拉魯濕地不同水質(zhì)類(lèi)別臨界值的劃分關(guān)系：Ⅰ類(lèi)水質(zhì)標(biāo)準(zhǔn)，P*≤3.961 3;Ⅱ類(lèi)水質(zhì)標(biāo)準(zhǔn)，3.961 3

15.582 4。

根據(jù)改進(jìn)的內(nèi)梅羅污染指數(shù)法對(duì)每個(gè)采樣點(diǎn)進(jìn)行評(píng)價(jià)，結(jié)果發(fā)現(xiàn)，枯水期，1號(hào)、2號(hào)、3號(hào)、4號(hào)、5號(hào)、6號(hào)、7號(hào)采樣

點(diǎn)的P*分別為3.983、1.563、0.927、0.931、1.431、0.627、0.935;平水期，1號(hào)、2號(hào)、3號(hào)、4號(hào)、5號(hào)、6號(hào)、7號(hào)采樣點(diǎn)的P*分別

為1.613、1.211、0.837、0.646、0.843、0.661、1.226;豐水期，1號(hào)、2號(hào)、3號(hào)、4號(hào)、5號(hào)、6號(hào)、7號(hào)采樣點(diǎn)的P*分別

為0.805、0.728、0.573、0.431、0.781、0.707、0.629。根據(jù)不同水質(zhì)類(lèi)別臨界值的劃分關(guān)系可以得出，枯水期，除1號(hào)采樣點(diǎn)位Ⅱ類(lèi)水質(zhì)標(biāo)準(zhǔn)外，其余采樣點(diǎn)皆為Ⅰ類(lèi)水質(zhì)標(biāo)準(zhǔn);平水期，所有采樣點(diǎn)皆為Ⅰ類(lèi)水質(zhì)標(biāo)準(zhǔn);豐水期，所有采樣點(diǎn)皆為Ⅰ類(lèi)水質(zhì)標(biāo)準(zhǔn)。

根據(jù)圖8可以得到，枯水期，西藏拉魯濕地除東北部表現(xiàn)為Ⅱ類(lèi)水質(zhì)外，其余地區(qū)皆為Ⅰ類(lèi)水質(zhì)，可見(jiàn)在枯水期拉魯濕地水質(zhì)綜合狀況區(qū)域分布特點(diǎn)為由東北部向西部逐漸變好;在平水期和豐水期，西藏拉魯濕地整體水質(zhì)表現(xiàn)為Ⅰ類(lèi)水質(zhì)。由此可見(jiàn)，西藏拉魯濕地在平水期和豐水期的水質(zhì)情況比枯水期更為良好，這進(jìn)一步表明濕地水體質(zhì)量的改變趨向與水位緊密相連。王瑋[27]通過(guò)開(kāi)展水生植物對(duì)污水中的污染物去除凈化能力進(jìn)行研究，發(fā)現(xiàn)水生植物對(duì)N、P的吸收率皆達(dá)到70%以上，表明水生植物能有效去除凈化污水中的N、P和有機(jī)物等污染物;薩茹拉[28]通過(guò)對(duì)媯水河水生植物群落修復(fù)其受損生境研究，結(jié)果表明水生植物對(duì)濕地水體質(zhì)量與生態(tài)環(huán)境具有明顯的凈化效果。因此，拉魯濕地在豐水期水生植物處于大量繁殖階段，通過(guò)水生植物除去污染物的作用得到加強(qiáng)，使得豐水期污染物濃度低于枯水期，同時(shí)豐水期濕地整體水位上升，水體本身稀釋功能得到增強(qiáng)[29]，故而進(jìn)一步使污染物濃度降低，使其水體質(zhì)量好于枯水期。

2.3 拉魯濕地不同時(shí)期水體監(jiān)測(cè)指標(biāo)間相關(guān)性分析

利用Pearson相關(guān)系數(shù)[30]分析各個(gè)監(jiān)測(cè)指標(biāo)之間是否具有同源性是學(xué)術(shù)研究中經(jīng)常用到的相關(guān)性分析方法，不同監(jiān)測(cè)指標(biāo)之間相關(guān)系數(shù)越高則表明其具有相近的污染源或者移動(dòng)特點(diǎn)。從西藏拉魯濕地不同時(shí)期各個(gè)因子之間的相關(guān)系數(shù)（表3）可以看出，枯水期，TN與NH3-N存在極顯著正相關(guān)（P<0.01），其余監(jiān)測(cè)指標(biāo)之間均無(wú)明顯的相關(guān)性，可以推斷出在枯水期TN和NH3-N這2個(gè)指標(biāo)間可能具有相同污染源;平水期，TP與NH3-N存在顯著正相關(guān)（P<0.05），TN與NH3-N存在極顯著正相關(guān)（P<0.01），其余監(jiān)測(cè)指標(biāo)之間均無(wú)明顯的相關(guān)性，可以推斷出在平水期TP、TN和NH3-N這3個(gè)指標(biāo)之間可能具有相似的污染源或存在相互轉(zhuǎn)移變化的作用;豐水期，TN與NH3-N存在極顯著正相關(guān)（P<0.01），其余監(jiān)測(cè)指標(biāo)之間均無(wú)明顯的相關(guān)性，可以推斷出在豐水期TN和NH3-N也許有著相同的污染源或存在相互移動(dòng)變化的作用。同時(shí)從表3可以得到，NH3-N與TN從枯水期到豐水期均存在極顯著正相關(guān)，表明濕地水體中TN主要以NH3-N的形態(tài)存在，該結(jié)果與文澤偉[31]的研究結(jié)果一致。并且從表3也可以得出，西藏拉魯濕地不同時(shí)期各個(gè)因子之間的相關(guān)性主要隨著濕地的整體水位和季節(jié)性水量的變化而變化，這與馬廣文等[29]的研究結(jié)果基本相似。

2.4 拉魯濕地水質(zhì)綜合狀況變化與污染源分析

從枯水期到豐水期，西藏拉魯濕地水體中DO、NH3-N、TP、COD和TN質(zhì)量濃度分別降低了61.05%、54.55%、50.00%、2.41%和50.00%（表1），并且其水質(zhì)綜合狀況由枯水期到豐水期逐漸變好。西藏拉魯濕地東北部（1號(hào)采樣點(diǎn)）在枯水期表現(xiàn)為Ⅱ類(lèi)水質(zhì)，而在平水期和豐水期則表現(xiàn)為Ⅰ類(lèi)水質(zhì)，其原因可能與東北部緊鄰城鄉(xiāng)居民生活區(qū)易產(chǎn)生生活污水，并且由于枯水期拉魯濕地水位處于較低水平，導(dǎo)致其污染物濃度增加，而在平水期和豐水期水位相比于枯水期處于較高水平，通過(guò)水體DO的增加和水量稀釋作用的增強(qiáng)進(jìn)而提高水質(zhì)綜合狀況。這與吳青梅等[32]的研究結(jié)果一致。因此，拉魯濕地北入水口水體質(zhì)量治理和附近的污染防治是下一步西藏拉魯濕地水體質(zhì)量治理和管控的重點(diǎn)。同時(shí)隨著拉薩市高原開(kāi)放型城鎮(zhèn)化發(fā)展不斷推進(jìn)，由于拉魯濕地地處拉薩市城區(qū)，人口的增加、農(nóng)業(yè)和生活污水的排放，以及高原旅游業(yè)加速推進(jìn)帶來(lái)的人類(lèi)活動(dòng)越來(lái)越頻繁等原因，也給拉魯濕地水環(huán)境質(zhì)量帶來(lái)更大挑戰(zhàn)。

3 結(jié)論

（1）西藏拉魯濕地從枯水期到豐水期不同監(jiān)測(cè)指標(biāo)間變化較大，主要表現(xiàn)為pH平均值由8.27降低至8.24，DO平均值由6.11 mg/L降低至2.38 mg/L，NH3-N平均值由0.22 mg/L降低至0.10 mg/L，TP平均值由0.02 mg/L降低至0.01 mg/L，COD平均值由4.97 mg/L降低至4.85 mg/L，TN平均值由0.26 mg/L降低至0.13 mg/L。總體而言，西藏拉魯濕地不同水體環(huán)境指標(biāo)濃度均表現(xiàn)為降低趨勢(shì)，其N(xiāo)H3-N、TP、COD和TN的平均值處于隨時(shí)間推移呈現(xiàn)降低的趨勢(shì)，而水體pH和DO的平均值表現(xiàn)為隨時(shí)間變化先升高再降低的趨勢(shì)。

（2）枯水期，西藏拉魯濕地除東北部表現(xiàn)為Ⅱ類(lèi)水質(zhì)外，其余地區(qū)皆為Ⅰ類(lèi)水質(zhì)，其水質(zhì)綜合狀況區(qū)域分布特點(diǎn)為由東北部向西部逐漸變好。平水期和豐水期，西藏拉魯濕地整體水質(zhì)表現(xiàn)為Ⅰ類(lèi)水質(zhì)。由此可見(jiàn)，西藏拉魯濕地在平水期和豐水期的水質(zhì)情況比枯水期表現(xiàn)更為良好。

（3）枯水期，西藏拉魯濕地TN與NH3-N存在極顯著正相關(guān)（P<0.01）。平水期，TP與NH3-N存在顯著正相關(guān)（P<0.05），TN與NH3-N存在極顯著正相關(guān)（P<0.01）。豐水期，TN與NH3-N存在極顯著正相關(guān)（P<0.01）。濕地水體中TN主要以NH3-N的形態(tài)存在。綜合來(lái)看，西藏拉魯濕地不同時(shí)期各個(gè)因子之間的相關(guān)性主要隨著濕地的整體水位和季節(jié)性水量的變化而變化。

（4）從枯水期到豐水期，西藏拉魯濕地水體中DO、NH3-N、TP、COD和TN質(zhì)量濃度分別降低了61.05%、54.55%、50.00%、2.41%和50.00%，并且其水質(zhì)綜合狀況由枯水期到豐水期逐漸變好。但是拉魯濕地東北部北入水口水體質(zhì)量治理和附近的污染防治是下一步西藏拉魯濕地水體質(zhì)量治理和管控的重點(diǎn)。

參考文獻(xiàn)

[1] 殷書(shū)柏，李冰，沈方.濕地定義研究進(jìn)展[J].濕地科學(xué)，2014，12（4）：504-514.

[2] COSTANZA R，D’ARGE R，DE GROOT R，et al.The value of the wold’s ecosystem services and natural capital[J].Nature，1997，387：253-260.

[3] 張永澤，王烜.自然濕地生態(tài)恢復(fù)研究綜述[J].生態(tài)學(xué)報(bào)，2001，21（2）：309-314.

[4] 陳瑞紅，莫德清，李金城，等.會(huì)仙巖溶濕地水質(zhì)監(jiān)測(cè)及評(píng)價(jià)[J].山東化工，2018，47（6）：156-160.

[5] 王佐霖，馬鵬飛，張衛(wèi)強(qiáng)，等.深圳壩光濕地地表水水環(huán)境質(zhì)量評(píng)價(jià)[J].林業(yè)與環(huán)境科學(xué)，2019，35（4）：9-17.

[6] 楊海江，鐘艷霞，羅玲玲，等.基于貝葉斯的星海湖濕地水質(zhì)評(píng)價(jià)及特征分析[J].節(jié)水灌溉，2018（4）：92-95，104.

[7] 徐好.南四湖水質(zhì)時(shí)空分布及評(píng)價(jià)研究[D].濟(jì)南：濟(jì)南大學(xué)，2019.

[8] 楊梅玲，胡忠軍，劉其根，等.利用綜合營(yíng)養(yǎng)狀態(tài)指數(shù)和修正的營(yíng)養(yǎng)狀態(tài)指數(shù)評(píng)價(jià)千島湖水質(zhì)變化（2007年—2011年）[J].上海海洋大學(xué)學(xué)報(bào)，2013，22（2）：240-245.

[9] 溫曉君，陳輝，白軍紅.基于模糊矩陣的衡水湖水環(huán)境質(zhì)量評(píng)價(jià)及分析[J].水土保持研究，2016，23（2）：292-296.

[10] 寇文杰，林健，陳忠榮，等.內(nèi)梅羅指數(shù)法在水質(zhì)評(píng)價(jià)中存在問(wèn)題及修正[J].南水北調(diào)與水利科技，2012，10（4）：39-41，47.

[11] 馬騰飛，邱一富.改進(jìn)的內(nèi)梅羅污染指數(shù)法在鑒江水質(zhì)評(píng)價(jià)中的應(yīng)用[J].廣東石油化工學(xué)院學(xué)報(bào)，2019，29（3）：91-94.

[12] 韓術(shù)鑫，王利紅，趙長(zhǎng)盛.內(nèi)梅羅指數(shù)法在環(huán)境質(zhì)量評(píng)價(jià)中的適用性與修正原則[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2017，36（10）：2153-2160.

[13] 吳喜軍，董穎，張亞寧.改進(jìn)的內(nèi)梅羅污染指數(shù)法在黃河干流水質(zhì)評(píng)價(jià)中的應(yīng)用[J].節(jié)水灌溉，2018（10）：51-53，58.

[14] 宗浩，王成善，黃川友，等.西藏拉薩市拉魯濕地的生態(tài)特征與退化機(jī)理的探討[J].西南民族大學(xué)學(xué)報(bào)（自然科學(xué)版），2005，31（1）：72-78.

[15] 李春.拉魯濕地生態(tài)環(huán)境及動(dòng)植物物種資源變化特征研究[J].自然資源學(xué)報(bào)，2005，20（1）：145-151.

[16] 路飛，許先鵬，鮮明睿，等.拉薩市城市規(guī)劃區(qū)濕地保護(hù)現(xiàn)狀與對(duì)策[J].林業(yè)調(diào)查規(guī)劃，2019，44（1）：143-147，152.

[17] 張文駒，張瀟瀟.拉魯濕地水質(zhì)評(píng)價(jià)與分析[J].安徽農(nóng)業(yè)科學(xué)，2014，42（5）：1477-1478，1502.

[18] 白永飛，黃偉，羅笑娟，等.拉魯濕地水中主要金屬元素分布分析評(píng)價(jià)[J].環(huán)境科學(xué)與技術(shù)，2016，39（6）：174-178.

[19] 巴桑，黃香，普布，等.拉魯濕地肉鞭蟲(chóng)群落特征及其水環(huán)境評(píng)價(jià)[J].濕地科學(xué)，2014，12（2）：182-191.

[20] 王佳俊，田瀚鑫，周磊，等.拉魯濕地水生植物群落多樣性與水環(huán)境因子的關(guān)系[J].環(huán)境科學(xué)，2020，41（4）：1657-1665.

[21] 周歡歡，劉引鴿，胡浩楠.寶雞市東沙河流域水質(zhì)污染評(píng)價(jià)[J].四川環(huán)境，2019，38（3）：30-35.

[22] 國(guó)家環(huán)境保護(hù)總局.地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)：GB 3838—2002[S].北京：中國(guó)環(huán)境科學(xué)出版社，2002.

[23] 李紅海，夏夢(mèng)雨，馮德金.襄陽(yáng)漢江國(guó)家濕地公園水質(zhì)分析及植物分布與水質(zhì)的關(guān)系[J].江漢大學(xué)學(xué)報(bào)（自然科學(xué)版），2019，47（6）：572-576.

[24] 曾凱，王家生，章運(yùn)超，等.華陽(yáng)河湖群水位變化對(duì)水質(zhì)的影響[J].長(zhǎng)江科學(xué)院院報(bào)，2020，37（8）：49-53，68.

[25] 李蘇，閆志宏，徐丹，等.改進(jìn)的內(nèi)梅羅指數(shù)法在水庫(kù)水質(zhì)評(píng)價(jià)中的應(yīng)用[J].科學(xué)技術(shù)與工程，2020，20（31）：13079-13084.

[26]

關(guān)云鵬.利用內(nèi)梅羅指數(shù)法模型評(píng)價(jià)地下水水質(zhì)的探討[J].山西水利科技，2012（1）：81-84.

[27] 王瑋.水生和陸生植物對(duì)污水中污染物的凈化功能及其機(jī)理[D].南寧：廣西大學(xué)，2019.

[28] 薩茹拉.媯水河受損生境特征分析與水生植物群落修復(fù)生境研究[D].鄭州：華北水利水電大學(xué)，2019.

[29] 馬廣文，王業(yè)耀，香寶，等.阿什河豐水期氮污染特征及其來(lái)源分析[J].環(huán)境科學(xué)與技術(shù)，2014，37（11）：116-120，155.

[30] 王磊，汪文東，劉懂，等.象山港流域入灣河流水體中重金屬風(fēng)險(xiǎn)評(píng)價(jià)及其來(lái)源解析[J].環(huán)境科學(xué)，2020，41（7）：3194-3203.

[31] 文澤偉.典型黑臭風(fēng)水塘水體pH值對(duì)底泥中氨氮、總磷釋放影響的研究[J].廣東化工，2019，46（6）：158-159.

[32] 吳青梅，羅慧東，孫國(guó)萍，等.典型感潮內(nèi)河涌水質(zhì)污染特征調(diào)查研究[J].環(huán)境科學(xué)學(xué)報(bào)，2011，31（10）：2210-2216.