石泉縣城區(qū)水源地水質(zhì)調(diào)查與評價

胡 娟

（陜西省石泉縣水利局，陜西 石泉 725200）

地表水資源是區(qū)域水資源的重要組成,特別是秦巴地區(qū)近60%人口的生活用水來自地表徑流。隨著人類活動強度加大，地表徑流的污染呈點、線、面日益蔓延，對居民飲水安全乃至環(huán)境發(fā)展帶來重大影響[1]。做好水資源調(diào)查與質(zhì)量測定的工作，有利于深入分析潛在的污染源，并為水資源管理和保護提供基礎(chǔ)信息。近年來針對水資源評價的方法可以分為多元統(tǒng)計分析、高階數(shù)學(xué)模數(shù)和神經(jīng)網(wǎng)絡(luò)等三種[2]，其中多元統(tǒng)計模型如TOPSIS法、熵權(quán)綜合指數(shù)法、投影尋蹤法等，高階數(shù)學(xué)模型如物元分析、科拓分析、D數(shù)理論等，神經(jīng)網(wǎng)絡(luò)應(yīng)用較多的有BP、RBF、SVM等[3]。其中多元統(tǒng)計分析具有操作簡單、評價過程直觀、方案可移植性強的特點，因而受到廣泛應(yīng)用。此外，GIS技術(shù)具有空間可視化能力也越來越多的應(yīng)用于水質(zhì)評價中。以石泉縣城區(qū)水源地水質(zhì)調(diào)查數(shù)據(jù)為基礎(chǔ)，運用因子分析等方法對水源地水質(zhì)進行綜合評價。

1 材料與方法

1.1 研究區(qū)概況

石泉縣地處秦巴山區(qū)中部漢江谷地中上游區(qū)域面積1525km2。區(qū)域?qū)贊h江流域，漢江自西向東穿流全境，南北兩側(cè)分屬巴山、秦嶺，形成谷地平原、山地、丘陵等復(fù)雜地形區(qū)，海拔介于328 m～2008 m之間。受季風環(huán)流與海陸位置效應(yīng)，形成北亞熱帶季風濕潤氣候，季候分明、溫暖濕潤，年均溫氣溫14.5℃，降水量873.9 mm，雨熱同季。石泉縣水資源豐富，以地表水資源位置，區(qū)域林草覆被程度高，達72%，自然生態(tài)環(huán)境良好。近年來隨著開發(fā)區(qū)、養(yǎng)殖業(yè)、工礦業(yè)的發(fā)展，局部徑流污染加劇。

1.2 水樣調(diào)查與測定分析

石泉縣城區(qū)飲用水水源地分兩處，一處位于城區(qū)西部漢江上游地段，一處位于湘子河水庫。根據(jù)水源地水域面積空間展布，遵循樣帶分布均衡性的要求，對漢江水源處采集36個樣點，對湘子河水源處采集24個樣點，并應(yīng)用GPS儀記錄樣點地理坐標。將采集的樣品放置干燥、潔凈的玻璃瓶中，除了pH值、電導(dǎo)率(EC)、溶解性總固體(TDS)等在實地測試外，其他水質(zhì)指標如總硬度(TH)、HCO3-、Cl-、SO42-、K＋、Ca2＋、Mg2＋等參考常規(guī)方法于實驗室測定。

1.3 數(shù)據(jù)分析

以測定的11個指標數(shù)據(jù)為基礎(chǔ)，先應(yīng)用3＋δ的方法剔除異常值隨后進行主成分分析，其具體流程如下[4]：

①應(yīng)用Zscore函數(shù)進行標準化，在此基礎(chǔ)上計算相關(guān)關(guān)系矩陣 R＝(rjk)：

式中，i為水質(zhì)樣點，j為水質(zhì)指標序列，xij為第i個水質(zhì)點的第j項指標含量值。

②計算矩陣特征值，根據(jù)各分量方差貢獻率大于80%的原則確定成分個數(shù)。若降維分析得到m個特征值，則存在m個維度，每個維度有q個指標，則各維度的重要性權(quán)值表示為：

式中，βij表示各指標因子分子系數(shù)，ei表示方差貢獻率，ωi代表各子系統(tǒng)重要性權(quán)值。

③對各維度進行主成分分析，基于權(quán)值進行加權(quán)求和，得到各水質(zhì)樣點綜合得分值，即為水質(zhì)評價量值[5]。

2 石泉縣城區(qū)水源地水質(zhì)分析

2.1 石泉縣水源地水質(zhì)統(tǒng)計特征

表1 石泉縣城區(qū)水源地水質(zhì)統(tǒng)計特征

應(yīng)用經(jīng)典統(tǒng)計學(xué)描述石泉縣城區(qū)水源地水質(zhì)指標基本特征。pH值是飲用水水質(zhì)重要衡量指標，一般要求其為中至弱堿性水，本區(qū)水體pH值介于7.09～8.32之間，平均值為7.59，變異系數(shù)僅為3.32%，表明區(qū)域pH分布差異性較小，且滿足飲水安全標準。EC是水體中礦物質(zhì)含量溶解度的反映，本區(qū)水體重EC含量在0.53 ms～199.72 ms之間，均值為62.17 ms，變異系數(shù)達112.36%，呈強變異，水體中EC含量差異性較大。水庫水體中水陽離子中 Cl-、SO42-、HCO3-的濃度值分別介于 0.36～2.58、1.48～43.2、30.4～85.9、17.1～41.5 mg/L 之間，平均值依次是1.18、9.82、83.8、26.23 mg/L，從變異系數(shù)來看，均呈中等程度變異。。其中，Ca2＋濃度變異系數(shù)較小，且濃度值較大，是地下水中最主要的陽離子。庫區(qū)水體重水陰離子Cl－、SO42-、HCO3-濃度值分別在 2.58～59.29、1.48～56.43 和 109.73～345.73 mg/L 之間，其平均水平依次為12.99、12.27、201.73 mg/L。TH的含量范圍介于185.98～17.97 mg/L之間，平均值為106.94 mg/L。變異系數(shù)顯示，除了EC、K＋和Cl-的含量呈強變異之外，其他指標均為中等程度變異。從水質(zhì)各項指含量來看，石泉縣城區(qū)水源地水質(zhì)相關(guān)指標均低于WHO(世界衛(wèi)生組織)飲用水標準限制，表明其水質(zhì)良好。

2.2 石泉縣城區(qū)水源地水質(zhì)主成分分析

將原始水質(zhì)數(shù)據(jù)予以標準化處理后，采用主成分分析提取方差貢獻率＞80的成分，共計有4特征成分，依次為5.32、4.14、3.22、1.09，與此同時其方差貢獻率依次達到39.51%、27.17%、14.86%、7.52%，累積方差貢獻率為89.06%，各因子得分系數(shù)旋轉(zhuǎn)矩陣如表2所示。

表2 主成分分析得分系數(shù)矩陣

將該系數(shù)值與方差貢獻率進行加權(quán)求和，即可得到水質(zhì)綜合指數(shù)中各項指標Y1的重要性系數(shù)，結(jié)果如表3所示。統(tǒng)計分析表明，Y1的值域范圍介于-0.76～1.86之間，其綜合得分越高，表明水體質(zhì)量越差，反之說明水質(zhì)越優(yōu)?；谄渲涤蚍秶瑒澐炙|(zhì)等級，若水質(zhì)綜合指數(shù)介于[-0.76，0)，說明水體質(zhì)量良好，其值介于[0,1]，表明水質(zhì)呈中等，當其值為(1,1.86]時，說明水體綜合質(zhì)量較差。

表3 綜合水質(zhì)指數(shù)Y1得分系數(shù)

2.3 石泉縣城區(qū)水質(zhì)空間可視化分析

在前述主成分分析的基礎(chǔ)上計算石泉縣城區(qū)水源地各水樣點的水體質(zhì)量分數(shù)，并根據(jù)其指數(shù)大小采用普通kriging插值進行空間可視化表達，得到結(jié)果如圖1所示?？芍瑵h江水源地水質(zhì)空間分布較為復(fù)雜，總體呈點狀離散分布，于水源地局部地區(qū)水質(zhì)較差，水質(zhì)綜合分數(shù)＞1；大部分地區(qū)水質(zhì)綜合分數(shù)在0～1之間，局域有斑點狀分布的水域其水質(zhì)綜合分數(shù)＜0。統(tǒng)計分析顯示，漢江水源地中40個水質(zhì)樣點中有16個水質(zhì)樣點屬于良好級別，有25個屬于中等，另外9個樣點水質(zhì)較差。從湘子河水質(zhì)分布來看，庫區(qū)北部和南部水質(zhì)較差，達到較差級別，中部及大部分水域?qū)儆诹己眉?。統(tǒng)計表明，湘子河水源地的24個水質(zhì)點中，有6個為較差等級，13個樣點為中等，另外5個為良好水平。從其分布形態(tài)來看，局部存在異常水質(zhì)樣點，結(jié)合實際調(diào)查發(fā)現(xiàn)，湘子河水源地局部水質(zhì)較差的樣點為人為活動強烈和地質(zhì)災(zāi)害頻發(fā)區(qū)；而漢江水源地水域局部水質(zhì)較差則是暗樁引起植物殘體富集所致。

圖1 漢江和湘子河水源地水質(zhì)綜合分數(shù)空間分布

3 結(jié)論

采用主成分分析和GIS技術(shù)相結(jié)合的方法對漢江和湘子河水庫兩處水源地水質(zhì)進行評價。結(jié)果表明，該水源地水體中的陰離子、陽離子、礦質(zhì)成分含量以及pH基本達到因素飲用水安全標準；整體水質(zhì)質(zhì)量處于中度以上較好的水平，但局部水質(zhì)較差，從其分布形態(tài)來看可知區(qū)域水質(zhì)存在部分點源污染，雖然污染程度低，但仍不可忽視，應(yīng)加強水環(huán)境治理。