基于內(nèi)梅羅指數(shù)與主成分分析的洪湖地區(qū)地下水水質(zhì)評(píng)價(jià)

劉 萌，陳世儉

(1.中國(guó)科學(xué)院 測(cè)量與地球物理研究所， 武漢 430077； 2.中國(guó)科學(xué)院大學(xué)， 北京 100049)

?

基于內(nèi)梅羅指數(shù)與主成分分析的洪湖地區(qū)地下水水質(zhì)評(píng)價(jià)

劉 萌1，2*，陳世儉1

(1.中國(guó)科學(xué)院 測(cè)量與地球物理研究所， 武漢 430077； 2.中國(guó)科學(xué)院大學(xué)， 北京 100049)

地下水污染現(xiàn)象在我國(guó)地下水資源豐富地區(qū)普遍存在.該文以洪湖地區(qū)豐水期、枯水期、平水期3個(gè)時(shí)期的氨氮、硝酸鹽氮、亞硝酸鹽氮、總錳、總鐵、高錳酸鉀指數(shù)和總磷等7個(gè)指標(biāo)為對(duì)象，運(yùn)用內(nèi)梅羅指數(shù)法與主成分分析方法，采用定性分析與定量分析相結(jié)合的方式，對(duì)8個(gè)樣點(diǎn)的的水環(huán)境質(zhì)量進(jìn)行綜合評(píng)價(jià).結(jié)果表明，臨近湖岸居民區(qū)與道路的潛水樣點(diǎn)5，受到來(lái)自水陸兩方面生產(chǎn)生活污染物的影響，水質(zhì)最差；位于臨近長(zhǎng)江的農(nóng)村居民區(qū)內(nèi)的潛水樣點(diǎn)7，主要受到來(lái)自居民生活污水的影響，但長(zhǎng)江切穿含水層頂板，江水與地下水水力聯(lián)系緊密，水體自凈好，水質(zhì)相對(duì)較好；位于農(nóng)田附近的承壓水樣點(diǎn)2，承壓水受外界污染少，在內(nèi)梅羅指數(shù)法中綜合得分最低，水質(zhì)相對(duì)較好.在豐水期高錳酸鉀指數(shù)與氨氮為主要污染物，在枯水期總鐵與總錳為主要污染物，在平水期總錳與硝酸鹽氮為主要污染物.長(zhǎng)江對(duì)含水層的影響、人類地表活動(dòng)、區(qū)域地質(zhì)狀況等都是可能影響洪湖地區(qū)地下水水質(zhì)的原因.

內(nèi)梅羅指數(shù)； 主成分分析； 洪湖地區(qū)； 地下水； 水質(zhì)

地下水資源是人類生產(chǎn)生活的重要資源.據(jù)統(tǒng)計(jì)，全國(guó)20%左右的供水來(lái)自開采地下水，淮河平原和華北平原生活飲用水中地下水占60%～70%[1].江漢平原東部地區(qū)，地下水可采資源量為每年4 000×106m3，每年總用水量約為10 322×106m3，其中，農(nóng)業(yè)用水量為7 520×106m3，工業(yè)用水量為2 231×106m3，生活用水量為571×106m3[2].隨著我國(guó)產(chǎn)業(yè)經(jīng)濟(jì)的發(fā)展，以及水資源利用技術(shù)的深入，我國(guó)水資源的開發(fā)得到了巨大發(fā)展，但是伴隨而來(lái)的問(wèn)題也逐漸凸顯.局部地區(qū)出現(xiàn)地下水超采，地面沉降、沿海地區(qū)海水入侵和地下水污染等問(wèn)題日益突出[3-6]，其中地下水污染問(wèn)題已經(jīng)成為不可忽視的難題.掌握地下水污染成分以及污染程度，對(duì)于我國(guó)地下水資源的長(zhǎng)期可持續(xù)開發(fā)具有重要意義.

水質(zhì)評(píng)價(jià)以定量的方式直觀表征水環(huán)境的質(zhì)量狀況，通過(guò)水質(zhì)評(píng)價(jià)可以了解水環(huán)境質(zhì)量的變化規(guī)律，為水環(huán)境質(zhì)量的管理和規(guī)劃提供依據(jù).目前，應(yīng)用于地下水水質(zhì)的評(píng)價(jià)的方法眾多，大致有：?jiǎn)我蜃釉u(píng)價(jià)法、綜合指數(shù)法、內(nèi)梅羅指數(shù)法、主成分分析法、灰色聚類法和物元可拓法等[7-12].上述評(píng)價(jià)方法各有特點(diǎn)，如單指標(biāo)評(píng)價(jià)方法簡(jiǎn)單、明了，選取各評(píng)價(jià)指標(biāo)中最低的合格率作為最終的評(píng)價(jià)結(jié)果，具有一定的合理性，但不能全面的反映水質(zhì)的整體狀況[7]；綜合指數(shù)法是對(duì)各污染指標(biāo)的相對(duì)污染指數(shù)進(jìn)行統(tǒng)計(jì)，得出代表水體污染程度的數(shù)值結(jié)果，能綜合反映出河流污染污染狀況，但沒有判明水質(zhì)類別[8]；內(nèi)梅羅指數(shù)法的數(shù)學(xué)過(guò)程簡(jiǎn)潔，評(píng)價(jià)時(shí)只需計(jì)算出綜合指數(shù)，再對(duì)照相應(yīng)的分級(jí)標(biāo)準(zhǔn)，便可知該評(píng)價(jià)區(qū)某環(huán)境要素的綜合環(huán)境質(zhì)量狀況，便于決策者做出綜合決策[9]；主成分分析法在確保系統(tǒng)原有數(shù)據(jù)信息量丟失最小的原則下，將多個(gè)變量的信息壓縮為幾個(gè)能反映原問(wèn)題特征的綜合變量指標(biāo)，從而避免了人為地確定各指標(biāo)權(quán)重的主觀隨意性，具有降維、簡(jiǎn)化變量之優(yōu)點(diǎn)，在太湖流域、武漢東湖等多區(qū)域水質(zhì)分析中取得了較好的效果[10，13-14].由于地下水中多種污染因子之間關(guān)系復(fù)雜，運(yùn)用單一方法評(píng)價(jià)難以解釋污染因子與水質(zhì)之間的關(guān)系，所以，本文選擇洪湖地區(qū)地下水為研究對(duì)象，運(yùn)用內(nèi)梅羅指數(shù)法和主成分分析法對(duì)區(qū)域地下水水質(zhì)進(jìn)行評(píng)價(jià)，探討其在豐水期、枯水期、平水期的污染程度，了解以及個(gè)采樣點(diǎn)之間、各水質(zhì)指標(biāo)之間的相互關(guān)系，明確各時(shí)期主要污染物，討論可能造成污染的原因，以便為控制地下水污染提供科學(xué)依據(jù).

1 研究區(qū)域與研究方法

1.1研究區(qū)域

洪湖地區(qū)南臨長(zhǎng)江，北臨漢江支流東荊河，位于地勢(shì)低洼的江漢平原四湖地區(qū)的東南部[15]，境內(nèi)水網(wǎng)密布，為典型的泛濫平原，屬亞熱帶季風(fēng)氣候，年平均降水量達(dá)1 100～1 300 mm[16].區(qū)內(nèi)水資源豐富，水產(chǎn)漁業(yè)成為當(dāng)?shù)亟?jīng)濟(jì)發(fā)展的重要產(chǎn)業(yè)[17]，地下水位高且地下水資源豐富，是區(qū)域農(nóng)業(yè)生產(chǎn)重要水源.

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

為監(jiān)測(cè)區(qū)域水質(zhì)變化，在區(qū)域內(nèi)選擇8個(gè)點(diǎn)位

(如圖1)：樣點(diǎn)1、樣點(diǎn)2、樣點(diǎn)3位于農(nóng)田附近；樣點(diǎn)4位于農(nóng)村居民區(qū)內(nèi)；樣點(diǎn)7、樣點(diǎn)8位于靠近長(zhǎng)江的農(nóng)村居民區(qū)內(nèi)；樣點(diǎn)5位于臨湖的農(nóng)村居民區(qū)內(nèi)；樣點(diǎn)6位于城鎮(zhèn)居民區(qū)內(nèi)；其中樣點(diǎn)5、樣點(diǎn)7為潛水，其余樣點(diǎn)為承壓水.于豐水期(2014年7月)、枯水期(2015年1月)、平水期(2015年5月)各采樣一次進(jìn)行監(jiān)測(cè).用聚乙烯瓶采集水樣并酸化保存，在室內(nèi)通過(guò)實(shí)驗(yàn)分析.選取總磷(TP)、氨氮(NH3-N)、亞硝酸鹽氮(NO2-N)、硝酸鹽氮(NO3-N)、高錳酸鉀指數(shù)(CODMn)、總鐵(Fe2+)和總錳(Mn2+)7項(xiàng)主要水質(zhì)指標(biāo)對(duì)區(qū)域水質(zhì)進(jìn)行分析，其中，總磷采用鉬酸銨分光光度法測(cè)定，氨氮采用納氏試劑法測(cè)定，亞硝酸鹽氮采用N-(1-萘基)-乙二胺光度法測(cè)定，硝酸鹽氮采用紫外分光光度法測(cè)定，高錳酸鹽指數(shù)采用高錳酸鉀酸性滴定法測(cè)定，總鐵采用鄰菲羅啉分光光度法測(cè)定，總錳采用高碘酸鉀分光光度法測(cè)定.

圖1 采樣點(diǎn)分布圖Fig.1 The distribution of sampling points

1.3研究方法

1.3.1內(nèi)梅羅指數(shù)法 內(nèi)梅羅指數(shù)法具備數(shù)學(xué)過(guò)程簡(jiǎn)潔、運(yùn)算方便等優(yōu)點(diǎn)，是以單因子評(píng)價(jià)法為基礎(chǔ)，對(duì)水體中的各種指標(biāo)參數(shù)進(jìn)行水質(zhì)評(píng)價(jià)的方法，評(píng)價(jià)地下水步驟如下[18-19].

1) 首先進(jìn)行水質(zhì)的單項(xiàng)組分評(píng)價(jià)，依據(jù)地下水質(zhì)量標(biāo)準(zhǔn)劃分單項(xiàng)組分所屬類別，對(duì)各類別按表

1分別確定單項(xiàng)組分評(píng)分值.

2) 利用內(nèi)梅羅指數(shù)計(jì)算公式計(jì)算綜合評(píng)分值F：

(1)

式中,F為綜合評(píng)分值；Fmax為單項(xiàng)組分Fi的最大值；Fave為各項(xiàng)組分評(píng)分值F的平均值；n為項(xiàng)數(shù).

表1 地下水質(zhì)量評(píng)分

3)根據(jù)F值按表2規(guī)定劃分地下水質(zhì)量級(jí)別.

表2 地下水質(zhì)量分級(jí)

1.3.2主成分分析法 主成分分析法是在確保不損失原有信息的前提下，將多種影響水質(zhì)的指標(biāo)重新組合成一組新的、相互無(wú)關(guān)的、較少的綜合指標(biāo)來(lái)反映指標(biāo)的信息，以達(dá)到降維、簡(jiǎn)化數(shù)據(jù)和突出關(guān)鍵信息的目的.

主成分分析法在水質(zhì)中的應(yīng)用是建立總評(píng)價(jià)指標(biāo)，評(píng)價(jià)各采樣點(diǎn)間的相對(duì)污染程度，并對(duì)各采樣點(diǎn)的污染程度進(jìn)行分級(jí)[20].假設(shè)有n個(gè)水樣，每個(gè)水樣檢測(cè)p個(gè)指標(biāo)，原始數(shù)據(jù)為n×p階矩陣：

(2)

1)對(duì)原始數(shù)據(jù)矩陣進(jìn)行標(biāo)準(zhǔn)化處理，消除量綱的影響：

(3)

式中，

(4)

計(jì)算采樣點(diǎn)各主成分得分：把各采樣點(diǎn)的標(biāo)準(zhǔn)化數(shù)據(jù)分別代入采樣點(diǎn)各主成分的表達(dá)式中，即可得出采樣點(diǎn)各主成分的得分Fj；以方差貢獻(xiàn)率為系數(shù)求得綜合得分F.得分值越大，表明污染程度越嚴(yán)重.

(5)

(6)

2 結(jié)果與分析

按照樣點(diǎn)的不同，在各樣點(diǎn)均選擇亞硝酸鹽氮、氨氮、總鐵、總錳、硝酸鹽氮、高錳酸鉀指數(shù)、總磷七個(gè)指標(biāo)，分豐水期、枯水期、平水期三個(gè)時(shí)期，對(duì)各樣點(diǎn)的水質(zhì)情況運(yùn)用內(nèi)梅羅指數(shù)法以及主成分分析法進(jìn)行分析.

2.1地下水化學(xué)特征

通過(guò)對(duì)研究區(qū)域不同時(shí)期各樣點(diǎn)地下水有關(guān)水化學(xué)參數(shù)進(jìn)行描述性統(tǒng)計(jì)分析，可以大致了解區(qū)域各污染指標(biāo)濃度狀況.分析計(jì)算結(jié)果見表3.

表3 地下水水化學(xué)參數(shù)描述性統(tǒng)計(jì)特征值

由表3可知，各樣點(diǎn)地下水亞硝酸鹽氮在0～0.418mg/L之間，在豐、枯水期略有超標(biāo)；氨氮含量在0.266～2.931mg/L之間，3個(gè)時(shí)期均有超標(biāo)，但超標(biāo)幅度不大；總鐵含量在0～7.662mg/L之間，在枯水期出現(xiàn)個(gè)別超標(biāo)值，且超標(biāo)幅度大；總錳含量在0.211～4.193mg/L之間，在枯水期出現(xiàn)個(gè)別超標(biāo)值，超標(biāo)幅度較大；硝酸鹽氮含量在0～17.735mg/L之間，在枯水期出現(xiàn)個(gè)別嚴(yán)重超標(biāo)值；高錳酸鉀指數(shù)含量在0.123～6.696mg/L之間，在枯水期與平水期出現(xiàn)超標(biāo)值，枯水期超標(biāo)幅度大；總磷含量在0.08～0.375mg/L之間，未出現(xiàn)超標(biāo).總體上各指標(biāo)超標(biāo)率較低，但對(duì)于區(qū)域水質(zhì)的評(píng)價(jià)還需要進(jìn)一步的分析.

2.2內(nèi)梅羅指數(shù)法

運(yùn)用內(nèi)梅羅指數(shù)法，根據(jù)不同污染物在不同類別水質(zhì)的標(biāo)準(zhǔn)濃度，以及各樣點(diǎn)實(shí)測(cè)污染濃度，計(jì)算各樣點(diǎn)水質(zhì)的內(nèi)梅羅指數(shù)綜合評(píng)分F值，得出各樣點(diǎn)在豐水期、枯水期、平水期的水質(zhì)分級(jí)結(jié)果如表4所示.

表4 3個(gè)時(shí)期地下水水質(zhì)的內(nèi)梅羅指數(shù)綜合評(píng)分F值以及質(zhì)量分級(jí)

由表中的質(zhì)量分級(jí)，可以得出在不同時(shí)期，所有8個(gè)樣點(diǎn)的水質(zhì)都為極差，而依據(jù)內(nèi)梅羅指數(shù)綜合得分F值，得出各樣點(diǎn)水質(zhì)在極差基礎(chǔ)上的個(gè)體差異.其中，樣點(diǎn)5在3個(gè)時(shí)期F值得分均為最高，水質(zhì)最差；樣點(diǎn)2在3個(gè)時(shí)期F值得分均為最低，水質(zhì)相對(duì)較好.依據(jù)內(nèi)梅羅指數(shù)法，我們得出了本區(qū)的大致水質(zhì)狀況，但是我們還需要對(duì)于影響各點(diǎn)的主要污染物進(jìn)行分析，所以我們繼續(xù)運(yùn)用主成分分析法對(duì)區(qū)域地下水水質(zhì)進(jìn)行分析.

2.3主成分分析法

通過(guò)統(tǒng)計(jì)軟件SPSS18.0對(duì)采樣所得到的數(shù)據(jù)進(jìn)行主成分分析，根據(jù)數(shù)據(jù)的特征值與主成分貢獻(xiàn)率，提取主成分并計(jì)算各項(xiàng)指標(biāo)的相關(guān)系數(shù)矩陣和特征值，完成各樣點(diǎn)的水質(zhì)綜合評(píng)價(jià).指標(biāo)與某一主成分的聯(lián)系系數(shù)的絕對(duì)值越大，則該主成分與指標(biāo)之間的聯(lián)系越緊密.

2.3.1豐水期水質(zhì)綜合分析 亞硝酸鹽氮、硝酸鹽氮、氨氮、總鐵、總磷指標(biāo)在第一主成分上的載荷較高，總錳、高錳酸鉀指數(shù)在第二主成分上的載荷較高.根據(jù)各主成分的得分與各主成分的方差貢獻(xiàn)率：

F=0.47902×F1+0.23781×F2，

計(jì)算可得水質(zhì)污染綜合得分F并進(jìn)行排序(表5).由表5可知這一時(shí)段樣點(diǎn)7得分最低、水質(zhì)相對(duì)較好；樣點(diǎn)3得分最高、水質(zhì)最差.

表5 豐水期各主成分得分及綜合得分

F1為第一主成分得分，F(xiàn)2為第二主成分得分，F(xiàn)為綜合得分.

2.3.2枯水期水質(zhì)綜合分析 分析枯水期數(shù)據(jù)，根據(jù)主成分載荷大小選取前三個(gè)主成分，亞硝酸鹽氮、硝酸鹽氮與第一主成分密切相關(guān)，氨氮、總鐵、高錳酸鉀指數(shù)在第二主成分上的載荷較高，總磷、總錳在第三主成分上的載荷較高，根據(jù)主成分特征向量，計(jì)算各主成分得分和水質(zhì)污染綜合得分(表6).

表6 枯水期各主成分得分及綜合得分

F1為第一主成分得分，F(xiàn)2為第二主成分得分，F(xiàn)3為第三主成分，F(xiàn)為綜合得分.

F=0.53323×F1+0.19665×F2+0.15813×F3.

按表6可知這一時(shí)段樣點(diǎn)4水質(zhì)相對(duì)較好，樣點(diǎn)7水質(zhì)最差.

2.3.3平水期水質(zhì)綜合分析 對(duì)2015年5月平水期數(shù)據(jù)進(jìn)行分析得出，氨氮、高錳酸鉀指數(shù)均與第一主成分密切相關(guān)，硝酸鹽氮、亞硝酸鹽氮、總鐵、總錳在第二主成分上的載荷較高，總磷在第三主成分上的載荷較高.根據(jù)各主成分的得分與各主成分的方差貢獻(xiàn)率：

F=0.40038×F1+0.26760×F2+

0.19022×F3，

得出2015年5月各主成分得分和水質(zhì)綜合得分及排序(表7).按照排名可知，樣點(diǎn)8水質(zhì)相對(duì)較好，樣點(diǎn)7水質(zhì)最差.

對(duì)不同時(shí)期各樣點(diǎn)的水質(zhì)進(jìn)行主成分分析，得出不同時(shí)期各樣點(diǎn)水質(zhì)的相對(duì)好壞，以及研究區(qū)域內(nèi)影響水質(zhì)的主要污染物，即第一主成分包含的污染物.但是由于在第一主成分上污染物指標(biāo)眾多，需要對(duì)其進(jìn)行進(jìn)一步的分析來(lái)確定對(duì)水質(zhì)影響最大的污染指標(biāo).

2.3.4主要污染物分析 依照時(shí)序變化，對(duì)各時(shí)期污染物進(jìn)行主成分分析，得到各時(shí)期污染指標(biāo)得分值，明確各時(shí)期最主要污染指標(biāo).

表7 平水期各主成分得分及綜合得分

F1為第一主成分得分，F(xiàn)2為第二主成分得分，F(xiàn)3為第三主成分，F(xiàn)為綜合得分.

總鐵Fe2+ 總錳Mn2+ 亞硝酸鹽氮NO2-N 硝酸鹽氮NO3-N 氨氮NH3-N 總磷TP 高錳酸鉀指數(shù)CODMn圖2 各時(shí)段污染指標(biāo)得分Fig.2 Composite score of pollutant indexes

依據(jù)3個(gè)時(shí)段各污染指標(biāo)的主成分分析綜合得分，結(jié)果如圖2，在豐水期，高錳酸鉀指數(shù)的綜合得分最高，為主要污染指標(biāo)，次要污染指標(biāo)為氨氮；在枯水期總鐵指標(biāo)的綜合得分最高，其次為總錳；在平水期，總錳與硝酸鹽氮得綜合得分排名靠前.綜合來(lái)說(shuō)，在豐水期高錳酸鉀指標(biāo)與氨氮為主要污染物，在枯水期總鐵與總錳為主要污染物，在平水期總錳與硝酸鹽氮為主要污染物.

2.4內(nèi)梅羅指數(shù)法與主成分分析法結(jié)果對(duì)比

依據(jù)采樣點(diǎn)的不同對(duì)3個(gè)時(shí)段的兩種水質(zhì)綜合得分求均值，并對(duì)得到的均值進(jìn)行歸一化處理，使其映射到0～1范圍內(nèi)進(jìn)行比較(圖3).潛水樣點(diǎn)5在兩種評(píng)價(jià)方法中綜合水質(zhì)得分均最高，污染最為嚴(yán)重；潛水樣點(diǎn)7在主成分分析法中水質(zhì)綜合得分最低，水質(zhì)相對(duì)較好；承壓水樣點(diǎn)2在內(nèi)梅羅指數(shù)法中綜合得分最低，水質(zhì)相對(duì)較好.

圖3 主成分分析法與內(nèi)梅羅指數(shù)法綜合得分均值排序Fig.3 Sort of average composite score of principal component analysis and Nemerow index method

兩種評(píng)價(jià)方法各有側(cè)重：內(nèi)梅羅指數(shù)法注重于對(duì)水質(zhì)的總體特征的評(píng)價(jià)，兼顧了污染最重因子和一般因子對(duì)于評(píng)價(jià)結(jié)果的影響，但忽視了各因子的重要程度；主成分分析法突出了水質(zhì)評(píng)價(jià)中的關(guān)鍵因子，強(qiáng)調(diào)主要污染因子對(duì)地下水水質(zhì)的影響，但忽略了次要因子的作用，因此兩種方法所得結(jié)果存在不同.內(nèi)梅羅指數(shù)法對(duì)區(qū)域水質(zhì)進(jìn)行定性分析，首先明確了區(qū)域水質(zhì)的等級(jí)，主成分分析法在此基礎(chǔ)上對(duì)于各樣點(diǎn)、各污染物進(jìn)行比較分析，得出了各時(shí)期各樣點(diǎn)的水質(zhì)相對(duì)好壞情況，以及各時(shí)期主要的污染物，兩種方法相結(jié)合對(duì)研究區(qū)域地下水水質(zhì)做出了較為合理的評(píng)價(jià).

3 討論

對(duì)可能造成本區(qū)污染的原因從以下幾個(gè)方面討論.

1) 本區(qū)域地勢(shì)低、水位較高，地表污染狀況對(duì)地下水的影響較大[16].在區(qū)域內(nèi)，豐水期長(zhǎng)江水位高于地下水水位，長(zhǎng)江補(bǔ)給地下水，且長(zhǎng)江切穿了隔水層頂層，與地下水有直接的水力聯(lián)系[21].如潛水樣點(diǎn)7位于臨近長(zhǎng)江的農(nóng)村居民區(qū)內(nèi)，該點(diǎn)地下水與長(zhǎng)江水交換密切，水資源更新速度快，自凈能力好，水質(zhì)較好.

2) 研究區(qū)域地勢(shì)較低，居民多開挖人工池塘，發(fā)展水產(chǎn)養(yǎng)殖業(yè).水產(chǎn)養(yǎng)殖產(chǎn)生大量殘餌，其中所含的氮磷以及耗氧有機(jī)物等導(dǎo)致水體氮磷污染，消耗水中溶解氧[22]，在溫度較高生物活動(dòng)頻繁的豐水期，高錳酸鉀指數(shù)成為主要污染指標(biāo).本區(qū)湖沖擊淤泥質(zhì)中有機(jī)質(zhì)較多[23]，地下水在厭氧細(xì)菌的作用下形成還原環(huán)境，易發(fā)生反硝化作用，水體中溶解性鐵、錳化合物與有機(jī)質(zhì)絡(luò)合形成穩(wěn)定的膠體體系，溶解度增加，遷移性增強(qiáng)[24]，影響水中鐵錳含量.在農(nóng)業(yè)活動(dòng)減少的枯水期，氮磷污染影響減弱，鐵錳成為主要污染指標(biāo).

3) 種植業(yè)是本區(qū)農(nóng)業(yè)中的支柱產(chǎn)業(yè)，豐水期正是農(nóng)事活動(dòng)最為繁忙的季節(jié).農(nóng)業(yè)上不合理施用化肥是造成地下水氮污染的重要原因[25].施用化肥增加了土壤表層硝酸鹽含量，并在雨季淋濾進(jìn)入地下水[26]，淺層地下水最先受到污染，深層地下水污染相對(duì)滯后[27].如位于農(nóng)田附近的承壓水樣點(diǎn)2，主要受農(nóng)業(yè)生產(chǎn)污染影響，但承壓水受外界影響較小，在內(nèi)梅羅指數(shù)法分析結(jié)果中中水質(zhì)相對(duì)較好.反硝化作用降低了地下水中硝酸鹽濃度，而低溶解氧、富含有機(jī)質(zhì)的的含水層適合于反硝化作用的發(fā)生[28].所以在水位較高降水較多的豐水期，還原狀態(tài)使氨氮成分較高，成為主要污染指標(biāo).而在水位較低的枯水期，或降水相對(duì)較少的平水期，硝態(tài)氮較多，硝酸鹽氮是這一時(shí)段的主要污染物.

4) 本區(qū)域缺乏完整的生活污水收集處理系統(tǒng)，生活污水隨意排放，生活垃圾的收集處理設(shè)施落后，相當(dāng)一部分生活垃圾隨意傾倒堆積，是潛在的主要污染源[29-30].如潛水樣點(diǎn)5位于湖岸居民區(qū)與道路之間，附近有眾多船舶停泊，居民生活生產(chǎn)污水影響地下水水質(zhì)，且附近水域易發(fā)生富營(yíng)養(yǎng)化，導(dǎo)致潛水污染，該點(diǎn)水質(zhì)在兩種評(píng)價(jià)結(jié)果中均為最差.

4 結(jié)論

地下水環(huán)境復(fù)雜多變，地下水水質(zhì)受到多方面因素的影響.本文采用內(nèi)梅羅指數(shù)法與主成分分析法相結(jié)合的分析方法，首先得出區(qū)域各點(diǎn)的污染等級(jí)指數(shù)，對(duì)整體污染情況有了定性的了解，再對(duì)各樣點(diǎn)相對(duì)污染情況和主要污染物進(jìn)行定量分析，較為合理地反映了區(qū)域地下水水質(zhì)狀況.臨近湖岸居民區(qū)與道路的潛水樣點(diǎn)5，受到來(lái)自水陸兩方面生產(chǎn)生活污染物的影響，水質(zhì)最差；位于臨近長(zhǎng)江的農(nóng)村居民區(qū)內(nèi)的潛水樣點(diǎn)7，主要受到來(lái)自居民生活污水的影響，但江水切穿含水層頂板，江水與地下水水力聯(lián)系緊密，水體自凈好，水質(zhì)相對(duì)較好；位于農(nóng)田附近的承壓水樣點(diǎn)2，主要受農(nóng)業(yè)生產(chǎn)污染影響，但承壓水受外界影響較小，在內(nèi)梅羅指數(shù)法中綜合得分最低，水質(zhì)相對(duì)較好.在豐水期高錳酸鉀與氨氮為主要污染物，在枯水期總鐵與總錳為主要污染物，在平水期總錳與硝酸鹽氮為主要污染物.

地質(zhì)條件、農(nóng)業(yè)生產(chǎn)、居民生活等因素綜合作用影響著地下水水質(zhì)，導(dǎo)致了污染物的產(chǎn)生與遷移.本區(qū)域應(yīng)調(diào)整農(nóng)業(yè)發(fā)展模式，加大科技投入，改變粗放式的經(jīng)營(yíng)，使農(nóng)業(yè)經(jīng)濟(jì)向可持續(xù)方向發(fā)展.合理管理農(nóng)業(yè)生產(chǎn)及生活污水的處理，逐步改善地下水水質(zhì)狀況.

[1] 文東光， 林良俊， 孫繼朝， 等. 中國(guó)東部主要平原地下水質(zhì)量與污染評(píng)價(jià)[J].地球科學(xué)——中國(guó)地質(zhì)大學(xué)學(xué)報(bào)， 2012， 37(2):220-228．

[2] 曾昭華. 江漢平原東部地區(qū)地下水資源的開發(fā)利用與保護(hù)[J].長(zhǎng)江流域資源與環(huán)境， 1996， 5(4):375-378．

[3] 石建省， 王 昭， 張兆吉， 等. 華北平原深層地下水超采程度計(jì)算與分析[J].地學(xué)前緣， 2010， 17(6):215-220．

[4] 殷躍平， 張作辰， 張開軍. 我國(guó)地面沉降現(xiàn)狀及防治對(duì)策研究[J].中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào)， 2005， 16(2):1-8．

[5] 章 斌， 宋獻(xiàn)方， 韓冬梅， 等. 運(yùn)用數(shù)理統(tǒng)計(jì)和模糊數(shù)學(xué)評(píng)價(jià)秦皇島洋戴河平原的海水入侵程度[J].地理科學(xué)， 2013， 33(3):342-348．

[6] 楊 彥， 于江云， 王宗慶， 等. 區(qū)域地下水污染風(fēng)險(xiǎn)評(píng)價(jià)方法研究[J].環(huán)境科學(xué)， 2013， 34(2):653-661．

[7] 陳仁杰， 錢海雷， 袁 東， 等. 改良綜合指數(shù)法及其在上海市水源水質(zhì)評(píng)價(jià)中的應(yīng)用[J].環(huán)境科學(xué)學(xué)報(bào)， 2010， 30(2):431-437．

[8] 孫 濤， 張妙仙， 李苗苗， 等. 基于對(duì)應(yīng)分析法和綜合污染指數(shù)法的水質(zhì)評(píng)價(jià)[J].環(huán)境科學(xué)與技術(shù)， 2014， 37(4):185-190．

[9] 張亞松， 張兆吉， 費(fèi)宇紅， 等. 地下水質(zhì)量綜合評(píng)價(jià)方法優(yōu)選與分析——以滹沱河沖洪積扇為例[J].水文地質(zhì)工程地質(zhì)， 2011， 38(1):6-10．

[10] 劉德林， 劉賢趙. 主成分分析在河流水質(zhì)綜合評(píng)價(jià)中的應(yīng)用[J].水土保持研究， 2006， 13(3):124-126．

[11] 王洪梅， 盧文喜， 辛 光， 等. 灰色聚類法在地表水水質(zhì)評(píng)價(jià)中的應(yīng)用[J].節(jié)水灌溉， 2007(5):20-22．

[12] 湯 潔， 李艷梅， 卞建民， 等. 物元可拓法在地下水水質(zhì)評(píng)價(jià)中的應(yīng)用[J].水文地質(zhì)工程地質(zhì)， 2005(5):1-5．

[13] 周雪麗， 丁賢榮， 程立剛. 太湖水質(zhì)站網(wǎng)數(shù)據(jù)的主成分分析應(yīng)用[J].地球信息科學(xué)， 2008， 10(2):142-146．

[14] 姚煥玫， 黃仁濤， 甘復(fù)興， 等. 用改進(jìn)的主成分分析法對(duì)東湖的水質(zhì)污染進(jìn)行評(píng)價(jià)[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版)， 2005， 30(8):732-735．

[15] 胡望斌， 王學(xué)雷， 陳世儉. 洪湖地區(qū)淺層承壓水動(dòng)態(tài)模擬研究[J].華中師范大學(xué)學(xué)報(bào)(自然科學(xué)版)， 2003， 37(1):119-122．

[16] 張 婷， 陳世儉， 傅嬌鳳. 四湖地區(qū)地下水“三氮”含量及時(shí)空分布特征分析[J].長(zhǎng)江流域資源與環(huán)境， 2014， 23(9):1295-1300．

[17] 姜?jiǎng)⒅荆?王學(xué)雷， 厲恩華， 等. 生態(tài)恢復(fù)前后的洪湖水質(zhì)變化特征及驅(qū)動(dòng)因素[J].濕地科學(xué)， 2012， 10(2):188-193．

[18] 張曉葉， 張永祥， 任仲宇， 等. 不同地下水水質(zhì)評(píng)價(jià)方法的比較及實(shí)例應(yīng)用[J].水資源與水工程學(xué)報(bào)， 2014， 25(4):98-101．

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

[20] 伊元榮， 海米提·依米提， 王 濤， 等. 主成分分析法在城市河流水質(zhì)評(píng)價(jià)中的應(yīng)用[J].干旱區(qū)研究， 2008， 25(4):497-501．

[21] 傅嬌鳳， 陳世儉， 張 婷. 四湖地區(qū)水質(zhì)評(píng)價(jià)與空間分析[J].華中師范大學(xué)學(xué)報(bào)(自然科學(xué)版)， 2013， 47(6)：868-872．

[22] 吳 偉， 范立民. 水產(chǎn)養(yǎng)殖環(huán)境的污染及其控制對(duì)策[J].中國(guó)農(nóng)業(yè)科技導(dǎo)報(bào)， 2014， 16(2):26-34．

[23] 鄧青軍， 唐仲華， 吳 琦， 等. 荊州市地下水動(dòng)態(tài)特征及影響因素分析[J].長(zhǎng)江流域資源與環(huán)境， 2014， 23(9):1215-1221．

[24] 陳心鳳， 邵衛(wèi)云， 宋瑞平， 等. 舟山飲用水源中鐵錳污染及其成因分析[J].地球與環(huán)境， 2011， 39(2):181-187．

[25] 劉冬碧， 范先鵬， 楊 利， 等. 江漢平原水稻肥水管理現(xiàn)狀與技術(shù)對(duì)策[J].湖北農(nóng)業(yè)科學(xué)， 2010， 49(8):1831-1835．

[26] 沃 飛， 陳效民， 吳華山， 等. 太湖流域典型地區(qū)農(nóng)村水壞境氮、磷污染狀況的研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)， 2007， 26(3):819-825．

[27] 劉宏斌， 李志宏， 張?jiān)瀑F， 等. 北京平原農(nóng)區(qū)地下水硝態(tài)氮污染狀況及其影響因素研究[J].土壤學(xué)報(bào)， 2006， 43(3):405-412．

[28] 陳新明， 馬 騰， 蔡鶴生， 等. 地下水氮污染的區(qū)域性調(diào)控策略[J].地質(zhì)科技情報(bào)， 2013， 32(6):130-149．

[29] WANG X L， LU Yonglong， HAN Jingyi， et al. Identification of anthropogenic influences on water quality of rivers in Taihu watershed[J].Journal of Environmental Sciences， 2007， 19(4):475-481．

[30] 凌 霄， 楊細(xì)平， 陳 滿， 等. 廣東省農(nóng)村生活污水治理現(xiàn)狀調(diào)查[J].中國(guó)給水排水， 2009， 25(8):8-10．

Groundwater quality assessment of Honghu Area based on the Nemerow Index and principal component analysis method

LIU Meng1,2, CHEN Shijian1

(1. Institute of Geodesy and Geophysics, Chinese Academy of Science, Wuhan 430077；2.University of Chinese Academy of Sciences, Beijing 100049)

Groundwater pollution is a widespread phenomenon in the areas with plenty of water resources in China. In the present work, comprehensive evaluation of 8 sampling points in Honghu Area of high-water, normal-water and low-flow period is taken based on Nemerow Index and principal component analysis. Seven indexes are used, including ammonia nitrogen, nitrate nitrogen, nitrite nitrogen, iron, manganese, potassium permanganate index and total phosphorus. The results indicate that the water quality of phreatic water sample 5 located near residential area and road on the lakeshore is the worst, due to the industrial and domestic pollutants from both water and land. The water quality of phreatic water sample 7, which is near the Yangtze River and mainly affected by domestic wastewater, is relatively better, as the Yangtze River cut through the aquifer top layer and the groundwater and river water are closely linked. The confined water sample 2 also shows relatively better water quality. During the high-water period, main pollutants are potassium permanganate index and ammonia nitrogen, while that of low-flow period changed to total iron and manganese. Additionally，in normal-water period total manganese and nitrate nitrogen are the main pollutants. These suggest that the the ground water quality might be affected by complicated factors, such as, influence of the Yangtze river on aquifer, human activities and regional geological conditions.

Nemerow Index; principal component analysis; Honghu Area; groundwater; water quality

2016-04-15.

國(guó)務(wù)院三峽辦項(xiàng)目(JJ2011-024).

1000-1190(2016)04-0633-08

X523;P33

A

*E-mail: liumenglmlmlm@126.com.