程琛 劉佳俊 楊凝

摘要：為給礦區(qū)水資源綜合利用提供信息，以安徽宿州兩煤礦蘆嶺和桃園煤礦塌陷水域?yàn)榈湫脱芯繀^(qū)，對(duì)蘆嶺9個(gè)水樣和桃園17個(gè)水樣進(jìn)行常規(guī)離子（F-、Cl-、SO2-4、Na+、K+、Mg2+、Ca2+、NO-2、NO-3、H2PO-4、HCO-3、CO2-3）和部分重金屬（Cd、Cr、Fe）進(jìn)行檢測(cè)，結(jié)果表明：除了F-離子處于劣Ⅴ類水平，大部分采樣點(diǎn)的Cr和Fe和其他元素均能滿足水質(zhì)標(biāo)準(zhǔn)的Ⅲ類標(biāo)準(zhǔn)。選定7個(gè)評(píng)價(jià)參數(shù)利用綜合污染指數(shù)法評(píng)價(jià)塌陷塘的水環(huán)境質(zhì)量，結(jié)果表明：蘆嶺礦塌陷水域綜合污染指數(shù)是0.886，桃園為0.873，氟（F-）的單因子污染指數(shù)分別是1.355、1.306，F(xiàn)-對(duì)污染貢獻(xiàn)最大，整個(gè)水域?qū)儆谥卸任廴舅健２捎靡蜃臃治龊拖嚓P(guān)性分析對(duì)其進(jìn)行污染來(lái)源分析，結(jié)果表明：主要污染來(lái)自F、Cr、Fe，外在污染源主要是煤矸石淋溶作用和煤礦開(kāi)采；主要來(lái)源分為三類，第一類是基于水力聯(lián)系的常規(guī)離子來(lái)源，第二類是大氣降水對(duì)含氟化合物的溶濾作用，第三類是金屬礦物的溶解和人類活動(dòng)的影響。

關(guān)鍵詞：宿州采煤塌陷水域；水質(zhì)評(píng)價(jià)；污染來(lái)源解析

中圖分類號(hào)： X142 文獻(xiàn)標(biāo)志碼：A文章編號(hào)：1672-1098（2017）04-0024-08

Abstract：To provide valuable information for comprehensive utilization of water resources in mining area， Luling and Taoyuan located in Suzhou city， Anhui province. two typical subsidence water areas of coal mine， have been taken to do the research. What have been detected in 9 water samples from Luning and other 17 water samples from Taoyuan were the conventional ions and some heavy metals， with detected conventional ions as F-， Cl-， SO2-4， Na+， K+， Mg2+， Ca2+， NO-2， NO-3， H2PO-4， HCO-3， CO2-3， and heavy metals as Cd， Cr and Fe. The results show that except F- ，which was at the inferior Class Ⅴ quality， the Cr and Fe in most samples and other elements were all in Class Ⅲ standards， which could totally meet the water quality standard. Also， 7 assessment parameters had been selected and comprehensive pollution index method had been used to evaluate the water environmental quality of subsidence pond， the results show that the comprehensive pollution indexes of Luling and Taoyuan subsidence water area were 0.886 and 0.873， the pollution indexes of single factor F-in Luling and Taoyuan were 1.355 and 1.306， and the conventional ion F- had the largest contribution to the pollution， and the pollution of the water area was at a middle level. And then the factor analysis and the correlation analysis were applied into the source analysis of pollution， the results show that the main pollution conventional ions were F， Cr， Fe， and the external pollution sources mainly came from coal gangue eluviation and coal mining； and the main sources were divided into three categories， the first one was the conventional ion based on hydraulic connection， the second one the atmospheric precipitation leached effects of fluoride compounds， the third one the dissolution of metal mineral and the influence of human activities.

Key words：suzhou mining subsidence water area； water quality evaluation；pollution source apportionment

煤礦采掘后導(dǎo)致的地表沉降是煤礦區(qū)最為重要的地質(zhì)災(zāi)害之一。對(duì)我國(guó)而言，每年開(kāi)采1.2×109t煤就會(huì)造成2.4×102km2的采空塌陷面積，據(jù)不完全統(tǒng)計(jì)，截止2015年，全國(guó)煤礦累計(jì)采空塌陷面積超過(guò)9 000km2，造成的損失超過(guò)650億元，到2020年安徽省兩淮礦區(qū)塌陷面積超1 000 km2 [1-4]。礦區(qū)的采煤塌陷水域是煤炭生產(chǎn)建設(shè)的產(chǎn)物。近年來(lái)，隨著煤炭開(kāi)采規(guī)模擴(kuò)大，采煤塌陷水域也不斷擴(kuò)展，礦區(qū)水環(huán)境不斷惡化。作為一種特殊的地表水體，塌陷水域?qū)ΦV區(qū)生態(tài)環(huán)境有著至關(guān)重要的作用：首先，隨著煤炭開(kāi)采及其規(guī)模的擴(kuò)大，配套選煤廠、電廠增加及其規(guī)模擴(kuò)大，充分并合理利用礦區(qū)水資源便成了解決礦區(qū)生產(chǎn)和生活用水需求量增加的有效途徑之一；其次，作為一個(gè)自然和社會(huì)的復(fù)合生態(tài)系統(tǒng)，塌陷區(qū)水質(zhì)狀況是該區(qū)域環(huán)境質(zhì)量的綜合反映。因此，充分認(rèn)識(shí)這些塌陷區(qū)水體，對(duì)水體進(jìn)行水質(zhì)評(píng)價(jià)和元素來(lái)源解析，對(duì)于合理開(kāi)發(fā)、利用塌陷區(qū)水域并使其生態(tài)環(huán)境和資源向可持續(xù)方向發(fā)展十分有利，從而也為地方政府合理有效治理塌陷區(qū)水域提供基礎(chǔ)資料和科學(xué)依據(jù)。

水質(zhì)指數(shù)法（WQI）是水質(zhì)評(píng)價(jià)常用方法，文獻(xiàn)[4]采用該方法分析河流水質(zhì)，文獻(xiàn)[5]則用于模擬地下水水質(zhì)，文獻(xiàn)[6]則用于評(píng)價(jià)海洋水體污染和富營(yíng)養(yǎng)化情況。塌陷水域作為特殊的地表水體，水質(zhì)污染指數(shù)法亦適用于此類水體水質(zhì)評(píng)價(jià)，但諸多學(xué)者對(duì)塌陷水域的研究多集中在塌陷塘富營(yíng)養(yǎng)化的生物研究方面。文獻(xiàn)[7-9]都將安徽淮南煤礦塌陷塘作為研究對(duì)象，評(píng)價(jià)塌陷水域水體富營(yíng)養(yǎng)化程度并做生態(tài)學(xué)研究。本文采用水質(zhì)指數(shù)法對(duì)塌陷水域的水質(zhì)理化特征進(jìn)行綜合評(píng)價(jià)并結(jié)合水質(zhì)結(jié)果進(jìn)行污染來(lái)源解析。

主成分分析、因子分析和相關(guān)性分析通常用來(lái)進(jìn)行元素來(lái)源解析[10-12]。在水質(zhì)評(píng)價(jià)基礎(chǔ)上，采用多種源解析分析方法相結(jié)合對(duì)主要污染物及其來(lái)源“追根溯源”。本論文中，采集來(lái)自宿州煤礦采煤塌陷區(qū)水域26組水樣，對(duì)其進(jìn)行常規(guī)離子含量測(cè)試和重金屬濃度測(cè)試，在進(jìn)行水質(zhì)評(píng)價(jià)時(shí)參照中國(guó)地表水水質(zhì)標(biāo)準(zhǔn)和WTO國(guó)際標(biāo)準(zhǔn)[13]，在此基礎(chǔ)上討論如何充分利用塌陷水體以達(dá)到礦區(qū)水資源綜合利用。合理選擇檢測(cè)項(xiàng)作為評(píng)價(jià)因子，利用綜合水質(zhì)指數(shù)進(jìn)行水質(zhì)評(píng)價(jià)，并采用旋轉(zhuǎn)因子分析和相關(guān)性分析進(jìn)行源解析。

1研究區(qū)概況

宿州位于安徽省東北部，以煤礦開(kāi)采為主，東西長(zhǎng)195km，南北長(zhǎng)151km，總面積9 787 km2。主要有五大煤礦：朱仙莊和蘆嶺煤礦位于東邊，桃園、祁南和祁東位于南面。蘆嶺煤礦建礦時(shí)間較長(zhǎng)，形成塌陷水域面積近10 km2，積水深度8 m左右，塌陷塘周邊土地主要用以農(nóng)耕種植，水體周邊堆有矸石，淮河主要支流沱河流經(jīng)于此但旱季時(shí)常干涸。桃園煤礦南鄰祁南礦，東鄰蘆嶺礦，與淮河支流澮河、沱河相距不遠(yuǎn)，地下水系活動(dòng)較為頻繁，水文地質(zhì)條件復(fù)雜，塌陷形成的積水面積大、深度大，且所研究區(qū)域周邊存在矸石復(fù)墾造田及水產(chǎn)養(yǎng)殖。蘆嶺和桃園兩個(gè)礦區(qū)塌陷水域較為典型，故以此塌陷塘水體為研究對(duì)象，分別采集這兩個(gè)礦區(qū)周邊塌陷水域水樣進(jìn)行測(cè)試評(píng)價(jià)，水樣采集時(shí)間為2014年5月。圖1（a）為蘆嶺煤礦塌陷區(qū)采樣點(diǎn)布設(shè)圖，圖1（b）為桃園煤礦塌陷區(qū)采樣點(diǎn)分布圖，分別設(shè)9個(gè)和17個(gè)采樣點(diǎn)，采樣點(diǎn)周邊布滿農(nóng)田、村莊、采礦廢石和煤矸石。

2研究方法

2.1綜合污染指數(shù)法與源解析方法

綜合污染指數(shù)法是一種很重要的水環(huán)境質(zhì)量評(píng)價(jià)方法，通過(guò)簡(jiǎn)單的數(shù)值代表水體污染等級(jí)。綜合污染指數(shù)P代表水體質(zhì)量，可以是算術(shù)平均指數(shù)、權(quán)重指數(shù)或者通過(guò)向量模型計(jì)算綜合指數(shù)。本論文采用加權(quán)平均計(jì)算綜合污染指數(shù)，該方法一般分為兩個(gè)步驟，首先根據(jù)每個(gè)參數(shù)選擇對(duì)應(yīng)的評(píng)價(jià)等級(jí)，計(jì)算權(quán)重和單因子指數(shù)，再根據(jù)污染指數(shù)選擇評(píng)價(jià)因子并進(jìn)行綜合污染指數(shù)計(jì)算[14]。由最終得到的綜合污染指數(shù)平均值，依照所對(duì)應(yīng)的水質(zhì)等級(jí)，將評(píng)價(jià)體系分五個(gè)等級(jí)。

式中：ci代表第i 個(gè)參數(shù)的測(cè)定濃度，si 表示第i 個(gè)參數(shù)的允許標(biāo)準(zhǔn)濃度，Wi 是第i 個(gè)參數(shù)的影響權(quán)重，n是所選擇的參數(shù)個(gè)數(shù)。當(dāng)計(jì)算單因子污染指數(shù)時(shí)，允許標(biāo)準(zhǔn)濃度的選擇依據(jù)相對(duì)應(yīng)等級(jí)的上限值確定，同時(shí)由單因子污染指數(shù)與污染指數(shù)總和比值確定參數(shù)影響權(quán)重。

對(duì)水樣數(shù)據(jù)進(jìn)行污染來(lái)源解析采用SPSS軟件進(jìn)行因子分析和相關(guān)性分析。因子分析時(shí)，先降維后采用最大方差的因子旋轉(zhuǎn)法得到各因子荷載；相關(guān)性分析時(shí)，采用基于Spearman秩相關(guān)方法，顯著性檢驗(yàn)選擇雙側(cè)T檢驗(yàn)。

2.2實(shí)驗(yàn)方法

水樣采集 24h 內(nèi)進(jìn)行0.22 μm濾膜抽濾預(yù)處理和2% 硝酸酸化處理。對(duì)水質(zhì)指標(biāo)的檢測(cè)，10種常規(guī)離子F-、 Cl-、 NO-2、NO-3、H2PO-4、SO2-4、Na+、 K+、Mg2+、Ca2+ 的濃度通過(guò)離子色譜儀測(cè)定，HCO-3和CO2-3 離子濃度通過(guò)傳統(tǒng)實(shí)驗(yàn)室方法酸堿滴定測(cè)得，重金屬含量Cd、Cr和Fe采用原子吸收石墨爐法測(cè)定。

水樣中常規(guī)陰陽(yáng)離子測(cè)定使用美國(guó)戴安公司 ICS-900 離子色譜儀：AS19 陰離子交換柱（250mm×4mm）、CS12 陽(yáng)離子交換柱（250mm×4mm）、C18陰陽(yáng)離子保護(hù)柱、MMS300型自動(dòng)循環(huán)再生抑制器、DS5型檢測(cè)器；碳酸鹽采用傳統(tǒng)實(shí)驗(yàn)室方法酸堿滴定；重金屬測(cè)試采用原子吸收分光光度法，實(shí)驗(yàn)室所用儀器型號(hào)為北京普析通用公司TAS-990AFG。

3結(jié)果與分析

3.1離子濃度

26組水樣中L1-L9來(lái)自位于宿州東南蘆嶺塌陷區(qū)，T1-T17來(lái)自宿州南桃園煤礦塌陷區(qū)，測(cè)定所得常規(guī)離子濃度和部分重金屬濃度如表1所示。

由表1可知：對(duì)于NO-2和 NO-3濃度變化范圍，桃園塌陷水域明顯比蘆嶺大；同一塌陷塘的H2PO-4含量變化不大，整體濃度變化桃園高于蘆嶺；HCO-3在所有水樣中濃度含量占絕對(duì)優(yōu)勢(shì)。

通過(guò)Origin軟件繪制主要常規(guī)陰陽(yáng)離子F-、Cl-、SO2-4、Na+、K+、 Mg2+、Ca2+ 和重金屬離子Cd、Cr、Fe濃度變化趨勢(shì)圖，如圖2所示。

由圖2（a）可知：F-濃度相對(duì)穩(wěn)定，兩個(gè)礦區(qū)塌陷塘的F-相差不大，都接近2.0mg/L，但超出水質(zhì)標(biāo)準(zhǔn)V類水質(zhì)限值1.5mg/L；Cl-含量波動(dòng)較小，水質(zhì)較好，濃度范圍20mg/L～40 mg/L；蘆嶺SO2-4含量較桃園高，且兩個(gè)煤礦塌陷水域均有含量起伏較大的情況。由圖2（b）可知：兩個(gè)塌陷塘水樣中陽(yáng)離子Mg2+和 Ca2+變化伴隨SO2-4的變化，且兩種離子濃度變化范圍相近；蘆嶺水樣的 Na++K+含量在兩個(gè)采樣點(diǎn)有較大波動(dòng)，且跟Cl-波動(dòng)點(diǎn)變化一致；由圖2（c）可知：重金屬含量變化較大，尤其是桃園塌陷塘水樣，其中Cr濃度范圍1.742～162.8μg/L，F(xiàn)e濃度范圍84.6～231.5μg/L。同一片塌陷水域部分采樣點(diǎn)水樣中離子變化如此之大可能由于額外的污染源，無(wú)論桃園還是蘆嶺礦區(qū)的塌陷塘，都是周邊農(nóng)田土地灌溉用水，各種農(nóng)藥化肥的使用經(jīng)雨水沖刷都將對(duì)周邊水體造成一定污染，同時(shí)采礦圍巖、廢石、煤矸石長(zhǎng)期堆放并在環(huán)境綜合作用下也會(huì)污染水體。

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

根據(jù)我國(guó)《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》（GB3838-2002），選取F-、Cl-、 NO-3、SO2-4、Cd、Cr 和 Fe 作為地表水水質(zhì)評(píng)價(jià)的污染參數(shù)，標(biāo)準(zhǔn)值選擇如表3所示：對(duì)于F-，所有水樣F-濃度均高于1.5mg/L，選擇V類水質(zhì)作為標(biāo)準(zhǔn)限值；對(duì)于Cl-、NO-3和 SO2-4，在飲用水標(biāo)準(zhǔn)限值基礎(chǔ)上，分別選擇250mg/L、10mg/L和 250mg/L 作為標(biāo)準(zhǔn)值；對(duì)于Cd，濃度在0.1～0.4μg/L范圍內(nèi)，選擇 II 類標(biāo)準(zhǔn)作為相應(yīng)標(biāo)準(zhǔn)，標(biāo)準(zhǔn)限值為1.0μg/L；對(duì)于Cr，根據(jù)濃度范圍的變化，選擇I類、II類、III類和V類四個(gè)等級(jí)對(duì)應(yīng)的限值作為標(biāo)準(zhǔn)值；Fe 標(biāo)準(zhǔn)值的選擇方式同Cr，選擇I類、II類、III類三個(gè)等級(jí)。污染等級(jí)的判定依據(jù)綜合污染指數(shù)即P值，污染等級(jí)分為5個(gè)水平，水體清潔程度如表2所示。

對(duì)蘆嶺和桃園塌陷水域水樣所選定的污染因子進(jìn)行單因子污染指數(shù)及因子權(quán)重計(jì)算，分別將蘆嶺和桃園塌陷水域作為整體，其綜合污染指數(shù)及其計(jì)算過(guò)程如表3和表4所示。由表3～表4可知：蘆嶺水樣污染因子的權(quán)重大小為F-> Cr> Fe> SO2-4> Cd> Cl-> NO-3，桃園為 F-> Cr> Fe> Cd> SO2-4> Cl-> NO-3，兩者權(quán)重大小差異僅存在于SO2-4和Cd 之間，前者SO2-4 影響大，后者Cd影響大，其他離子影響權(quán)重一致。由各離子濃度平均值計(jì)算得到的綜合污染指數(shù)值、兩個(gè)礦區(qū)采樣點(diǎn)綜合污染指數(shù)值和兩個(gè)塌陷水體各自作為整體的綜合污染指數(shù)均在0.7

通過(guò)相關(guān)性分析和因子分析進(jìn)行元素來(lái)源解析。表5是采用斯皮爾曼法的相關(guān)性分析結(jié)果，圖3為因子分析的因子載荷圖，結(jié)果顯示蘆嶺塌陷水域所研究元素有三個(gè)共同來(lái)源，桃園有四個(gè)共同來(lái)源。因子分析時(shí)經(jīng)因子旋轉(zhuǎn)后，元素共同來(lái)源解釋總方差分別為92.0% 和81.5%，表明元素來(lái)源因子總貢獻(xiàn)率高達(dá)92.0% 和81.5%，除此之外可能的所有影響因素所占比例較小。蘆嶺塌陷水域三種主要來(lái)源因子所占比例分別為42.6%、33.7% 和 15.8%，桃園四種來(lái)源因子權(quán)重比例分別為23.2%、20.7%、19.7% 和17.9%。

對(duì)于桃園水樣，盡管因子分析結(jié)果顯示有四個(gè)來(lái)源因子，但除去F-離子的影響，主要貢獻(xiàn)因子為HCO-3、Na+、Ca2+、Mg2+、Cl-和 SO2-4，這與蘆嶺水域是一致的。因此，元素來(lái)源分為三類，一類是常規(guī)離子來(lái)源，一類是外在含氟化合物溶解與富集，另外就是重金屬污染源。對(duì)于F-含量較高，與許光泉等人對(duì)淮北平原淺層地下水研究結(jié)論一致[15]，但對(duì)宿州地區(qū)城市地下水及河流地表水中F-均有測(cè)定[16-17]，含量均在1.0mg/L以下，表明塌陷水體F-受外在來(lái)源影響。

因此，元素可能性來(lái)源可總結(jié)為以下幾個(gè)方面：常規(guī)陰陽(yáng)離子來(lái)源，即來(lái)源因子1可能為地表水與地下水水力聯(lián)系及轉(zhuǎn)化作用，地下水與斜長(zhǎng)石和碳酸鹽巖水巖相互作用導(dǎo)致水體中含有部分常規(guī)陰陽(yáng)離子，經(jīng)水力作用進(jìn)入到地表水層[18]。因子2即氟離子來(lái)源主要為大氣降水對(duì)氟化物的淋溶作用；因子3即重金屬來(lái)源主要考慮為地表徑流、地下水和地表水之間的水力聯(lián)系、地表水體周邊含重金屬礦物巖石溶解和人為活動(dòng)影響[19]，其中后兩者影響比較顯著，這由重金屬Cr 和 Fe變異系數(shù)接近100% 即可證實(shí)。

4結(jié)論

1） 長(zhǎng)期煤礦開(kāi)采導(dǎo)致蘆嶺和桃園礦區(qū)部分地表塌陷形成塌陷水域，采集的26組水樣中常規(guī)離子F-、 Cl-、SO2-4、Na+、K+、Mg2+、Ca2+、NO-2、NO-3、H2PO-4、HCO-3、CO2-3均被檢出，重金屬Cd、Cr和Fe含量也經(jīng)測(cè)出。對(duì)于氟離子，兩個(gè)礦區(qū)塌陷水體離子含量均超標(biāo)，但滿足高氟地區(qū)農(nóng)業(yè)灌溉水水質(zhì)標(biāo)準(zhǔn)；重金屬含量較之常規(guī)離子含量波動(dòng)較大，桃園塌陷區(qū)部分水體存在Cr污染。

2） 選定F-、Cl-、NO-3、SO2-4、Cd、Cr 和 Fe作為評(píng)價(jià)因子，對(duì)水體進(jìn)行綜合評(píng)價(jià)，結(jié)果表明：蘆嶺塌陷水域綜合污染指數(shù)P值為0.886，桃園為0.873，介于第III等級(jí)0.7 和1.0之間，水體存在中度污染情況。其中起主要作用的為F、Cr 和 Fe，權(quán)重較大，污染來(lái)源考慮為煤矸石淋溶和煤炭開(kāi)采活動(dòng)。

3） 結(jié)合相關(guān)性分析、水化學(xué)特征分析和水污染評(píng)價(jià)進(jìn)行元素來(lái)源及污染源分析，因子分析結(jié)果表明：兩個(gè)礦區(qū)塌陷水域具有共同的來(lái)源，可歸結(jié)為常規(guī)離子來(lái)源、氟化物來(lái)源、重金屬來(lái)源，具體表現(xiàn)為各類水體間相互轉(zhuǎn)化的水力聯(lián)系是塌陷水域水體中常規(guī)離子的主要來(lái)源，含氟化合物淋溶作用是氟化物主要來(lái)源，含重金屬礦物巖石的溶解和人類活動(dòng)是重金屬主要來(lái)源。

參考文獻(xiàn)：

[1]LIU JZ， HAN DJ， GU ZK. Key Controlling Regions and Main Measures for Soil and Water Conservation in Coal Mining Area [J]. Mining Research &Development;， 2009， 29（1）：65-67.

[2]徐良驥，嚴(yán)家平，高永梅. 安徽省兩淮礦區(qū)地質(zhì)環(huán)境治理技術(shù)[J]. 煤田地質(zhì)與勘探， 2007， 35（6）：37-40.

[3]王福琴. 安徽省兩淮采煤塌陷區(qū)的現(xiàn)狀、存在問(wèn)題及治理措施建議[J]. 安徽地質(zhì)，2010，20（4）：291-293.

[4]GE WQ. Countermeasures for Coal Mining Geologic Subsidence Disaster of China[J]. Energy of China，2004，26（10）：27-29.

[5]GUPTA RC，GUPTA AK，SHRIVASTAVA RK.Assessment and management of water quality of Kshipra River in Ujjain City （Madhya Pradesh）， India[J]. Journal of Environmental Science & Engineering， 2013， 55 （2）：189-196.

[6]SINGH SK，SRIVASTAVA PK，SINGH D，et al. Modeling groundwater quality over a humid subtropical region using numerical indices， earth observation datasets， and X-ray diffraction technique： a case study of Allahabad district， India[J]. Environmental Geochemistry and Health， 2015，37（1）：157-180.

[7]LIU SG，LOU S，KUANG CP，et al. Water quality assessment by pollution-index method in the coastal waters of Hebei Province in western Bohai Sea， China[J]. Marine Pollution Bulletin，2011，62（10）：2 220-2 229.

[8]桂和榮，宋曉梅，王振紅，等. 煤礦塌陷塘環(huán)境生態(tài)學(xué)研究[M]. 北京：地質(zhì)出版社，2009.

[9]張冰，嚴(yán)家平，范廷玉，等. 采煤塌陷水域富營(yíng)養(yǎng)化評(píng)價(jià)與分析——以淮南張集和顧橋礦為例[J]. 煤炭科學(xué)技術(shù)，2012，32（1）：159-161.

[10]YI QT，WANG XM，WANG TT，et al. Eutrophication and nutrient limitation in the aquatic zones around Huainan coal mine subsidence areas，Anhui，China[J].Water Science & Technology，2014，70（5）： 878-887.

[11]SUN LH，GUI HR，PENG WH.Heavy metals in groundwater from the coal bearing aquifer in northern Anhui Province，China：Concentrations and usability[J]. Journal of Chemical and Pharmaceutical Research，2013，5（12）：1 349-1 353.

[12]SUN LH，GUI HR.Hydro-chemical evolution of groundwater and mixing between aquifers：a statistical approach based on major ions[J].Appl Water Sci，2015 （5）：97-104.

[13]SUN LH，PENG WH.Groundwater in rural area of Suzhou，northern Anhui Province，China：Heavy metal concentrations and statistical analysis[J]. Journal of Chemical and Pharmaceutical Research，2014，6（3）：1 501-1 505.

[14]SUN LH，LIU XH，Cheng chen. Quality of water from subsidence area：a case study in the Luling coal mine， northern Anhui Provice， China[J]. Water Practice & Technology， 2015， 10 （4）： 777-786.

[15]B PRASAD，J M BOSE.Evaluation of the Heavy Metal Pollution Index for Surface and Spring Water Near a Limestone Mining Area of the Lower Himalayas [J]. Environmental Geology， 2001， 41（1-2）：183-188.

[16]許光泉，劉進(jìn)，朱其順，等. 安徽淮北平原淺層地下水中氟的分布特征及影響因素分析[J].水資源與水工程學(xué)報(bào)，2009，20（5）：9-13.

[17]程琛. 礦區(qū)城市地下水水質(zhì)特征與元素來(lái)源分析——以宿州市為例[J]. 宿州學(xué)院報(bào)，2015，30（10）：109-112.

[18]陳松，桂和榮，林曼利，等. 淮北煤田河流水化學(xué)特征及意義[J]. 地球與環(huán)境，2016，44（4）：414-421.

[19]張傳奇，溫小虎，高猛，等. 萊州灣東岸地下水化學(xué)及鹽化特征[J]. 海洋通報(bào)，2014，33（5）：584-591.

[20]徐良驥，嚴(yán)家平，高永梅. 煤礦塌陷水域水環(huán)境現(xiàn)狀分析及綜合利用——以淮南礦區(qū)潘一煤礦塌陷水域?yàn)槔齕J]. 煤炭學(xué)報(bào)，2009，34（7）：933-937.

（責(zé)任編輯：李麗，范君）