開偉萌，梁秀娟，肖長來，孟凡傲，楊偉飛，王 鴿

（1.吉林大學(xué)地下水資源與環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室，長春130021；2.吉林大學(xué)水資源與水環(huán)境重點(diǎn)實(shí)驗(yàn)室，長春130021；3.吉林大學(xué)油頁巖地下原位轉(zhuǎn)化與鉆采技術(shù)國家地方聯(lián)合工程實(shí)驗(yàn)室，長春130021；4.吉林大學(xué)新能源與環(huán)境學(xué)院，長春130021）

0 引 言

近年來，隨著城市化和工業(yè)化進(jìn)程的加快，我國地下水污染問題日趨嚴(yán)重。其成因可簡要?dú)w納為原生水質(zhì)問題和人類活動(dòng)導(dǎo)致的水質(zhì)問題［1］。對(duì)地下水水質(zhì)進(jìn)行檢測(cè)和評(píng)價(jià)，識(shí)別其主要污染源、空間分布及污染風(fēng)險(xiǎn)，在地下水治理和開發(fā)利用中具有重要意義。

以往的水質(zhì)評(píng)價(jià)方法有單因子指數(shù)法、內(nèi)梅洛指數(shù)法、模糊數(shù)學(xué)評(píng)價(jià)法、灰色聚類法等［2，3］。這些方法可以對(duì)當(dāng)前水質(zhì)進(jìn)行評(píng)價(jià)與分級(jí)，但難以識(shí)別地下水污染源，且忽略了地下水污染的空間特性。前人為確定地下水污染源做出了大量的研究。唐軍等［4］將哈斯圖解法與綜合評(píng)價(jià)法結(jié)合用于地下水優(yōu)控污染物識(shí)別；馮茜等［5］利用因子分析法識(shí)別了德惠市淺層地下水的污染源，并分析了污染源的空間分布特征；郭濤等［6］利用正定矩陣因子分解法識(shí)別了拉林河流域地下水的主要污染源。曹陽等［7］利用因子分析法和系統(tǒng)聚類法對(duì)水質(zhì)影響因素進(jìn)行評(píng)價(jià)。綜合來看，因子分析法，不但可以從宏觀上描述水質(zhì)因子，識(shí)別出污染源，還可以利用因子得分將研究區(qū)各污染因子及綜合污染因子的空間分布情況在圖中直觀表達(dá)，是識(shí)別地下水污染源的理想方法。

本文以松原市為研究對(duì)象，選取2017年11月58 口潛水監(jiān)測(cè)井的水質(zhì)檢測(cè)數(shù)據(jù)進(jìn)行分析研究，首次采用因子分析法確定松原市淺層地下水污染因子，分析污染源分布特征，并結(jié)合USEPA 的非致癌風(fēng)險(xiǎn)模型找出污染因子中的高風(fēng)險(xiǎn)因子，為研究區(qū)有關(guān)部門治理和開發(fā)利用地下水提供依據(jù)。

1 研究區(qū)概況及數(shù)據(jù)來源

1.1 研究區(qū)概況

松原市位于松嫩平原以南，東與長春市、四平市接壤，西與白城市接壤，屬溫帶大陸性氣候。區(qū)內(nèi)東部為由王府-伏龍泉砂礫石臺(tái)地和松拉河間地塊組成的高平原，其第四系孔隙潛水含水層為下更新統(tǒng)砂礫石含水巖組（Q1）和中更新統(tǒng)沖積砂礫石含水巖組（Q2）；西部為低平原，其地勢(shì)低平開闊，向二松河谷方向緩傾，其第四系孔隙潛水含水層為上更新統(tǒng)沖湖積粉細(xì)砂、砂礫石含水巖組（Q3）；東部高平原與西部低平原之間為河谷平原，其第四系孔隙潛水含水層為全新統(tǒng)沖湖積、湖沼積淤泥質(zhì)亞黏土含水巖組（Q4）［8］，如圖1（a）所示。區(qū)內(nèi)淺層地下水的補(bǔ)給來源主要為降水入滲補(bǔ)給，在西部低平原區(qū)，由于徑流較為緩慢，淺層地下水的排泄方式主要以蒸發(fā)和人工開采為主。

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

選用松原市2017年11月58組淺層地下水水質(zhì)檢測(cè)數(shù)據(jù)進(jìn)行分析，取樣分布，如圖1（b）所示。檢測(cè)的潛水含水層厚度在7～40 m 之間，按均勻分布的原則，每個(gè)取樣點(diǎn)取三瓶水樣，取樣完后立即用錫紙密封送往譜尼測(cè)試（吉林分公司）化驗(yàn)，檢測(cè)結(jié)果真實(shí)可靠。利用Piper 圖法確定研究區(qū)地下水化學(xué)類型時(shí)選取鉀、鈣、鈉、鎂、碳酸根、重碳酸根、硫酸根、氯離子8項(xiàng)指標(biāo)，利用因子分析法識(shí)別淺層地下水污染源時(shí)選取pH、鋇、鐵、鉬、鈉、鋅、氟化物、氯化物、硝酸鹽、硫酸鹽、溶解性總固體（TDS）、總硬度、耗氧量、氨氮等18項(xiàng)指標(biāo)。

圖1 研究區(qū)地形地貌及采樣分布圖Fig.1 Topography and sampling distribution map

2 研究方法與步驟

2.1 研究方法

本研究采用spss 軟件對(duì)地下水各離子指標(biāo)進(jìn)行統(tǒng)計(jì)性分析；運(yùn)用Aquachem軟件繪制Piper圖，以確定淺層地下水水化學(xué)類型。采用因子分析法識(shí)別研究區(qū)地下水中的污染源；將因子分析提取的污染因子代入非致癌物風(fēng)險(xiǎn)模型，計(jì)算風(fēng)險(xiǎn)指數(shù)。因子分析數(shù)學(xué)模型與非致癌物風(fēng)險(xiǎn)模型如下。

因子分析數(shù)學(xué)模型［1］：

式中：f1，f2，…，fk為因子變量，a11，a12，…，apk為因子xp與變量fk之間的相關(guān)系數(shù)，也稱為載荷；εi為特殊因子，代表公因子以外的影響因素，實(shí)際分析時(shí)可忽略不計(jì)。

非致癌物風(fēng)險(xiǎn)模型［9］：

式中：RI為非致癌物所致的個(gè)人平均風(fēng)險(xiǎn)，無量綱；DI為經(jīng)飲水途徑所導(dǎo)致的單位日均暴露劑量，mg/（kg·d）；RiDI為經(jīng)飲水途徑的參考劑量，mg/（kg·d）；CI為水中污染物的濃度，mg/L；IR為攝入率，一般取值2.1 L/d；EF為暴露頻率，一般取值365 d/a；ED為暴露持續(xù)時(shí)間，a；BW為平均體重，一般取值70 kg；AT平均壽命，取70 a［9］。

2.2 具體步驟

2.2.1 數(shù)據(jù)標(biāo)準(zhǔn)化

將選取的各項(xiàng)指標(biāo)的數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，其公式如下。

式中：xi為第i個(gè)水樣的原始測(cè)量值；x′i為第i個(gè)水樣去除量綱后的結(jié)果；s為其標(biāo)準(zhǔn)差。

2.2.2 數(shù)據(jù)檢驗(yàn)

為檢驗(yàn)變量之間是否有足夠的相關(guān)性，選用KMO 檢驗(yàn)和Bartlett球型檢測(cè)法，作為檢驗(yàn)方法。檢驗(yàn)結(jié)果如表1所示。

表1 KMO 和巴特利特檢驗(yàn)Tab.1 KMO and Bartlett tests

檢測(cè)結(jié)果顯示，KMO 檢驗(yàn)結(jié)果為0.698，大于0.6，巴特利特檢驗(yàn)結(jié)果為0，小于0.05，檢驗(yàn)結(jié)果適合做因子分析。

2.2.3 因子提取

提取其中幾個(gè)能反映原始變量絕大部分信息的公因子，其特征值應(yīng)大于1，為使公因子具有良好的代表性，其累計(jì)方差值應(yīng)大于70%。

2.2.4 因子旋轉(zhuǎn)

采用最大方差法對(duì)矩陣進(jìn)行旋轉(zhuǎn)，旋轉(zhuǎn)后的載荷矩陣應(yīng)使各個(gè)因子的載荷方差之和最大，以便對(duì)評(píng)價(jià)內(nèi)容進(jìn)行求解與分析。

2.2.5 計(jì)算因子得分

因子得分是指所提取的公共因子f1，f2，…，fk在每個(gè)取樣點(diǎn)的具體得分，每個(gè)因子得分由以下因子得分公式給出。

以上步驟借助spss軟件中的因子分析模塊實(shí)現(xiàn)。

2.2.6 因子風(fēng)險(xiǎn)研究

采用USEPA 的非致癌物風(fēng)險(xiǎn)模型［9］對(duì)因子分析法所提取的污染因子進(jìn)行風(fēng)險(xiǎn)指數(shù)計(jì)算，以EPA 推薦的最大可接受風(fēng)險(xiǎn)度（1.0×10-6）［10］為參照依據(jù)，研究分析各污染因子的風(fēng)險(xiǎn)程度。

3 結(jié)果與分析

3.1 區(qū)域水化學(xué)特征

通過分析研究區(qū)淺層地下水Piper 圖，如圖2 所示，可以看出研究區(qū)水化學(xué)類型主要以Ca-HCO3型水，CaNa-HCO3型水為主，少部分地區(qū)為Ca-Cl-HCO3型水。由表2（研究區(qū)主要水化學(xué)指標(biāo)統(tǒng)計(jì)表）可知，該區(qū)淺層地下水的pH均值為7.64，呈弱堿性。TDS 和總硬度的變異系數(shù)較大，說明在全區(qū)范圍內(nèi)礦化度相差較大，應(yīng)存在地下水富集的高污染地區(qū)，研究區(qū)Cl-、NO2-、NO3-、NH4+、Fe、As、Zn 等離子的標(biāo)準(zhǔn)差和變異系數(shù)較大，說明部分地區(qū)相關(guān)離子的含量偏高，可能是潛在的污染源。

表2 研究區(qū)主要水化學(xué)指標(biāo)統(tǒng)計(jì)表Tab.2 Statistical table of main hydrochemical indexes in the study area

圖2 淺層地下水Piper圖Fig.2 Piper diagram of shallow groundwater

3.2 地下水污染源識(shí)別及空間分布情況

在對(duì)選取的18 項(xiàng)指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理，KMO 檢驗(yàn)和Bartlett 球型檢測(cè)后，提取5 個(gè)公因子，累計(jì)方差貢獻(xiàn)率為75.546%，如表3 所示。為使各公因子的典型代表變量更加突出，將因子荷載矩陣進(jìn)行正交旋轉(zhuǎn)，旋轉(zhuǎn)因子荷載矩陣，見表4。

表3 總方差解釋Tab.3 Total variance interpretation

表4 旋轉(zhuǎn)因子載荷矩陣Tab.4 Rotation factor load matrix

將各公因子得分及綜合因子得分，與對(duì)應(yīng)的坐標(biāo)結(jié)合，運(yùn)用ArcGIS 中的克里金插值法得到淺層地下水各公因子與綜合污染因子得分的空間分布圖，如圖3 所示。圖中某區(qū)域單項(xiàng)因子得分較高，說明該地區(qū)某項(xiàng)離子污染程度較為嚴(yán)重；綜合因子得分較高，說明該地區(qū)整體離子濃度偏高，可能存在多項(xiàng)離子污染超標(biāo)的現(xiàn)象。

f1因子的方差貢獻(xiàn)率為38.151%，是區(qū)域淺層地下水的主要成分，f1因子的主要代表變量有溶解性總固體（TDS）、氯化物、總硬度、硝酸鹽氮、鈉離子、硫酸鹽，將f1命名為“主成分”因子。f1因子的得分在區(qū)域內(nèi)呈西高東低，得分較高的地區(qū)為乾安縣及長嶺縣西部區(qū)域，如圖3（a）所示。f1因子的得分與水文地質(zhì)條件、水動(dòng)力條件以及人類活動(dòng)有關(guān)。松原市西部地區(qū)為低平原區(qū)，干旱半干旱氣候，地勢(shì)低洼，地下水水位較低，周邊的地下水與巖石長期發(fā)生溶濾作用攜帶大量離子向此處運(yùn)移，在長期蒸發(fā)濃縮作用下，地下水中的離子濃度不斷增高，導(dǎo)致該地區(qū)f1因子得分較高。其次，生活污水和工業(yè)污水中的酸和鹽類物質(zhì)排入地表水后，由地表水下滲到地下水中，與巖層發(fā)生陽離子交替作用，使巖層中的鈣、鎂離子析出，這使得地下水的總硬度逐漸增高［11］。

f2因子的方差貢獻(xiàn)率為13.036%，主要代表變量有亞硝酸鹽氮和氨氮。f2因子得分自西向東逐漸降低，其中長嶺縣及乾安縣西部及寧江區(qū)因子得分較高，如圖3（b）所示。從空間分布來看，f2因子得分較高的地區(qū)主要分布在西部低平原區(qū)及人口密集的寧江區(qū)，這是由于亞硝酸鹽氮和氨氮主要和農(nóng)藥化肥的使用及污水的排放有關(guān)。在低平原地區(qū)，由于地勢(shì)平坦開闊，農(nóng)田較為集中，單位面積農(nóng)藥化肥的使用量遠(yuǎn)高于東部地區(qū)，過量的氮肥進(jìn)入土壤后隨雨水進(jìn)入地下水，造成地下水中氮含量逐漸增高。在寧江區(qū)，工廠排放的工業(yè)廢水及集中排放的生活污水存在大量的氮化物，處理未達(dá)標(biāo)的污水被排入第二松花江河道內(nèi)后，在遷移時(shí)發(fā)生下滲和越流補(bǔ)給現(xiàn)象，使得市區(qū)附近的淺層地下水中氮含量逐漸升高。

f3因子的方差貢獻(xiàn)率為10.128%，主要代表變量有鉬、鐵、和有機(jī)物。f3因子主要分布在乾安縣西部地區(qū)，如圖3（c）所示。參照《地下水質(zhì)量標(biāo)準(zhǔn)》（GB/T14848-2017），鉬和三氯甲烷的離子濃度均在III類標(biāo)準(zhǔn)以內(nèi)，屬正常范圍。地下水中鐵元素的含量與TDS 的含量呈正相關(guān)，與Ph 的值呈負(fù)相關(guān)，同時(shí)還與地下水中的氧化還原環(huán)境有直接關(guān)系，當(dāng)?shù)叵滤幱谶€原環(huán)境時(shí)鐵元素更容易從土壤中析出。在西部低平原區(qū)，鐵元素主要來源于土壤及整個(gè)含水層［12］。

f4因子的方差貢獻(xiàn)率為7.724%，主要代表變量是氟化物。f4因子主要分布于長嶺縣和乾安縣兩地，其中長嶺縣因子得分最高，如圖3（d）所示。氟元素含量主要受水文地質(zhì)條件和原生地質(zhì)環(huán)境影響。長嶺、乾安兩地，屬低平原區(qū)，為干旱半干旱氣候，其地形低洼，沉積物以沖湖積、風(fēng)積為主，含水層滲透性差，地下水交替緩慢，蒸發(fā)排泄為其主要排泄方式，易于離子的富集。其次，周邊裸露的氟礦物（方解石、石膏、長石類礦物和螢石等）經(jīng)過長期的風(fēng)化、淋溶、并借助水動(dòng)力條件富集于此。發(fā)源于內(nèi)蒙古大興安嶺德福勒罕山北麓的霍林河，是運(yùn)輸大興安嶺含氟巖石的主要載體，為該地區(qū)源源不斷的提供氟離子［13］。

圖3 各公因子與綜合污染因子得分空間分布圖Fig.3 Spatial distribution of the scores of common factors and comprehensive pollution factors

f5因子的方差貢獻(xiàn)率為6.507%，主要代表變量是砷元素。f5因子分布于松原市絕大部分地區(qū)，如圖3（e）所示。參照《地下水質(zhì)量標(biāo)準(zhǔn)》（GB/T14848-2017），砷的離子濃度在III 類標(biāo)準(zhǔn)以內(nèi)，屬正常范圍。原生地質(zhì)作用是地下水中砷的主要來源，含砷巖石經(jīng)風(fēng)化和雨水沖蝕等過程釋放可溶態(tài)砷吸附于土壤中，又在雨水作用下經(jīng)過長時(shí)間的下滲進(jìn)入地下水中。另外，地下水中的鐵元素會(huì)和固態(tài)砷產(chǎn)生鐵砷循環(huán)耦合作用，地下水中鐵元素含量越高，越容易促進(jìn)砷的釋放［14］。

通過公式f綜合=（38.151f1+13.036f2+10.128f3+7.724f4+6.057f5）/75.546 計(jì)算松原市淺層地下水污染綜合得分f。結(jié)果顯示，長嶺縣西部地區(qū)的太平川鎮(zhèn)、北正鎮(zhèn)，以及乾安縣西部地區(qū)永字鎮(zhèn)、大布蘇鎮(zhèn)污染最為嚴(yán)重，長嶺縣縣區(qū)、扶余縣縣區(qū)、及前郭爾羅斯蒙古族自治縣哈拉毛都鎮(zhèn)、長龍鄉(xiāng)、東三家子鄉(xiāng)水質(zhì)較好，如圖3（f）所示。綜上所述，綜合污染因子得分較高地區(qū)主要為低平原區(qū)，其主要受原生地質(zhì)條件、水文地質(zhì)條件及地下水動(dòng)力條件共同作用；其次，人類生產(chǎn)生活所排放的生活污水、工業(yè)污水及過度施用的農(nóng)藥也是造成該地淺層地下水污染的原因。

3.3 因子風(fēng)險(xiǎn)研究

在污染源識(shí)別的基礎(chǔ)上，選取氟、鐵、硝酸鹽氮、氨氮4個(gè)污染因子進(jìn)行風(fēng)險(xiǎn)指數(shù)計(jì)算。依據(jù)美國環(huán)保署公眾健康評(píng)估手冊(cè)［15］，硝酸鹽氮、氨氮、氟、鐵的參考劑量如表5所示。將指標(biāo)含量及各參數(shù)代入非致癌風(fēng)險(xiǎn)模型計(jì)算，結(jié)果如圖4所示。

圖4 各因子風(fēng)險(xiǎn)指數(shù)Fig.4 Risk index of each factor

表5 非致癌物參考劑量 mg/(kg·d)Tab.5 Reference dose of non-carcinogens

以EPA 推薦的最大可接受風(fēng)險(xiǎn)度（1.0×10-6）［12］為參照依據(jù)，在58 個(gè)點(diǎn)的檢測(cè)樣本中，氟化物有13 個(gè)點(diǎn)的風(fēng)險(xiǎn)指數(shù)高于最大可接受風(fēng)險(xiǎn)度，超標(biāo)率達(dá)22.41%，最大超標(biāo)倍數(shù)為2.48倍；硝酸鹽氮有6 個(gè)點(diǎn)的風(fēng)險(xiǎn)指數(shù)高于最大可接受風(fēng)險(xiǎn)度，超標(biāo)率達(dá)10.34%，最大超標(biāo)倍數(shù)為3.40 倍；鐵元素有1 個(gè)超標(biāo)，超標(biāo)率為1.72%，氨氮無超標(biāo)。結(jié)果表明地下水中氟化物、硝酸鹽氮的風(fēng)險(xiǎn)指數(shù)超標(biāo)率偏高，鐵元素和氨氮的風(fēng)險(xiǎn)指數(shù)超標(biāo)率較低。綜上所述，氟化物為風(fēng)險(xiǎn)最大的污染因子，其次是硝酸鹽氮類，鐵和氨氮因子的風(fēng)險(xiǎn)性較小。

4 結(jié) 論

（1）松原市淺層地下水總體上呈弱堿性，部分地區(qū)的礦化度較高。地下水中主要陽離子有Ca+、Na+、Mg2+，主要陰離子有HCO3-、Cl-。區(qū)內(nèi)地下水水化學(xué)類型主要為：Ca-HCO3、CaNa-HCO3型水。

（2）影響松原市淺層地下水水質(zhì)的主要因子有TDS 和硝酸鹽氮，其方差貢獻(xiàn)率達(dá)38.151%；其次是亞硝酸鹽氮、氨氮和氟化物；其方差貢獻(xiàn)率總和達(dá)20.76%；鐵和砷的方差貢獻(xiàn)率之和在17%左右。

（3）污染因子風(fēng)險(xiǎn)研究表明，區(qū)內(nèi)氟化物是風(fēng)險(xiǎn)最大的污染因子，其次是硝酸鹽氮類，鐵和氨氮因子的風(fēng)險(xiǎn)較小。

（4）區(qū)內(nèi)淺層地下水污染較嚴(yán)重的區(qū)域主要集中在乾安縣與長嶺縣西部地區(qū)，其主要受原生地質(zhì)環(huán)境作用和離子的遷移-富集作用影響。另外污水排放及化肥農(nóng)藥的過量施用也對(duì)淺層地下水水質(zhì)造成污染。 □