高潛水位采煤沉陷區(qū)水質評價與污染因子識別

胡 林，陳永春，徐燕飛，李 兵，王 錦，安士凱2,，陳 晨，孫宏杰，苗 偉

(1.煤礦生態(tài)環(huán)境保護國家工程實驗室，安徽 淮南 232001；2.安徽省煤礦綠色低碳發(fā)展工程研究中心，安徽 淮南 232001；3.煤炭開采國家工程技術研究院，安徽 淮南 232001)

安徽淮南礦區(qū)是我國14 個億噸煤炭基地、6 大煤電基地之一，同時又是我國東部重要的糧油生產基地，屬于典型的煤糧復合區(qū)。由于煤層厚度較大、地下水位埋藏深度淺、區(qū)域降水豐沛、地表水系發(fā)達等原因，煤炭地下開采易造成地表沉陷，形成相對獨立或有連通的沉陷積水區(qū)域。淮南礦區(qū)沉陷積水面積大、范圍廣，水環(huán)境問題突出。沉陷區(qū)水環(huán)境質量監(jiān)測評價及污染因子識別是解決水環(huán)境問題的基礎，可為采煤沉陷區(qū)污染水體的防治和管控提供有效和針對性的依據(jù)，對正確分析當前的礦區(qū)水環(huán)境形勢和問題，開展沉陷區(qū)水資源的分類以及水資源保護具有重要作用，同時也可為地方政府部門制定環(huán)保決策提供依據(jù)。

隨著社會對采煤沉陷積水區(qū)水體環(huán)境的重視，針對采煤沉陷區(qū)水環(huán)境問題的研究與報告也越來越多。一段時期以來，學界對沉陷區(qū)水域水質指標的研究集中于DO、NH3-N、CODcr、TP、TN 等常規(guī)理化指標[1-5]。近年來，針對危害人類健康的As、Hg、Cu、Cr、F-等水體重金屬及無機陰離子指標成為研究焦點[6-9]。針對沉陷區(qū)水域的水質評價，研究學者開展了大量基礎性研究及應用，目前常用的評價方法有單因子評價法、綜合污染指數(shù)法和內梅羅污染指數(shù)法等[10-13]。針對沉陷區(qū)水域的污染因子識別也得到了廣泛的研究，例如：采用因子分析和相關性分析對宿州煤礦塌陷區(qū)進行污染來源分析[1]。

雖然許多學者針對采煤沉陷積水區(qū)水質做了大量研究，也取得了豐碩的成果，但大多數(shù)研究區(qū)域多集中在少數(shù)幾個礦區(qū)，或研究的水質指標只集中在少數(shù)幾個參數(shù)，缺少大范圍、長系列、多參數(shù)的沉陷區(qū)整體水質現(xiàn)狀研究，尤其對高潛水位礦區(qū)沉陷區(qū)水域典型污染因子的研究與識別較少?；诖耍P者選取淮南礦區(qū)為研究對象，采集15 對礦井、102 片沉陷積水區(qū)、175 個點位的水樣，開展常規(guī)、重金屬、無機陰離子、有機物4 類13 項指標的測試分析，分析沉陷區(qū)水質特征，總結水質分布規(guī)律，識別污染因子，開展水質評價及水體富營養(yǎng)化評價，以期為地方政府制定生態(tài)環(huán)保決策提供數(shù)據(jù)支撐，也為國內其他高潛水位采煤沉陷區(qū)的水環(huán)境治理工作提供技術借鑒。

1 材料與方法

1.1 研究區(qū)概況

安徽淮南礦區(qū)內煤層厚且發(fā)育穩(wěn)定，產狀平緩，受地下潛水位較高、多煤層重復采動、煤層上覆第四系巨厚松散層、地表河網(wǎng)密布等因素影響，采動后的地表形成了大面積采煤沉陷積水區(qū)。近10 年來，隨著煤礦開采活動的進行，沉陷區(qū)水域面積不斷擴大，水環(huán)境問題得到廣泛關注。截至目前，淮南礦區(qū)范圍內采煤沉陷區(qū)面積共計約253.33 km2，常年積水面積約占一半。采煤沉陷區(qū)在各礦井內均有分布，在八公山一帶、潘謝礦區(qū)潘集和新集、謝橋至劉莊形成3 個呈北西向延伸的、連片分布的較大沉陷區(qū)，這3 個沉陷區(qū)均是煤炭開采程度較高的礦井所在區(qū)域。以淮河為界，淮河以南主要包括李一礦、李二礦、新莊孜礦、李嘴孜礦、謝一礦等關閉礦井；淮河以北分布了謝橋礦、張集礦、顧北礦、顧橋礦、丁集礦、朱集東礦、潘三礦、潘一礦、潘二礦等礦井，除潘一礦外均為在產礦井。研究區(qū)溝渠交錯、河系水網(wǎng)發(fā)達，按照沉陷水體與水系是否溝通的原則，將沉陷區(qū)分為封閉型采煤沉陷區(qū)和開放型采煤沉陷區(qū)2 種類型。

1.2 樣品采集與測試

本次野外采樣范圍為淮南礦區(qū)采煤形成的沉陷積水區(qū)，東西長約70 km，南北傾斜寬約25 km。經(jīng)全面梳理淮南礦區(qū)采煤沉陷區(qū)范圍、分布特征，依據(jù)《水質采樣技術指導》要求，按照“以礦劃分、分片采樣”的思想將研究區(qū)采煤沉陷區(qū)劃分為15 片沉陷區(qū)水域，共布設175 個采樣點，如圖1 所示，其中沉陷區(qū)采樣點145 個、河流采樣點30 個。采樣點布設原則如下：對于封閉型采煤沉陷區(qū)，在沉陷積水區(qū)上中下位置共布設采樣點1～3 個；對于開放型采煤沉陷區(qū)，在來水源頭(河流、湖泊)的上下游布置采樣點1～2 個，在沉陷積水區(qū)布設采樣點1～3 個。

圖1 研究區(qū)采樣點分布Fig.1 Distribution of sampling points in the study area

為了提高采樣效率和精度，本次采用“大疆M300”型無人機進行采樣作業(yè)。由采樣人員操控無人機飛至指定位置，采樣深度控制在水面以下20～50 cm。與傳統(tǒng)人工采樣方法相比，運用無人機采樣省時省力，操作簡單，避免了地形和水域環(huán)境的干擾，采樣的效率和科學性大大提高。在每個采樣點采集4 L 水樣，NH3-N指標測試樣品加硫酸酸化至pH<2，重金屬指標測試樣品中每升水加入10 mL 硝酸進行酸化處理，葉綠素測試樣品利用錫紙膜做好避光處理。歷時20 d，完成175 個點位的水樣采集工作，采樣結束后將樣品立即運至實驗室進行相關水質指標檢測，測試數(shù)據(jù)近3 000 個。

1.3 研究方法

依據(jù)《地表水質量標準》和《漁業(yè)水質標準》，將監(jiān)測的水質指標分為4 類13 項。(1)常規(guī)指標：溶解氧(DO)、氨氮(NH3-N)、總磷(TP)、總氮(TN)、化學需氧量(COD)；(2) 重金屬指標：砷(As)、汞(Hg)、銅(Cu)、鉻(Cr)；(3) 無機陰離子指標：氟離子(F-)；(4) 有機類指標：甲基對硫磷、馬拉硫磷、樂果。指標測試過程主要使用紫外可見分光光度計(氨氮、總磷等常規(guī)指標)、原子熒光儀(砷、汞等重金屬)、原子吸收儀(銅等重金屬)、離子色譜儀(無機陰離子)及氣相色譜儀(甲基對硫磷等農藥殘留)等儀器。

2 沉陷區(qū)水質及水體富營養(yǎng)化評價

2.1 沉陷區(qū)水質評價

2.1.1 方法優(yōu)選

查閱文獻發(fā)現(xiàn)，單因子評價法[14-15]存在一定的片面性，不足以完整反映區(qū)域性沉陷區(qū)的水質現(xiàn)狀；綜合污染指數(shù)法適合評價指標污染程度相當?shù)乃w，但不能直觀判斷水質類別；內梅羅污染指數(shù)法[16-18]在一定程度兼顧了其他水質較好的參評因子對總體結果的貢獻，能判斷綜合水質是否達到功能區(qū)目標，故本次采用內梅羅污染指數(shù)法進行研究區(qū)水質評價，評價標準見表1。計算公式如下：

表1 內梅羅污染指數(shù)法水質評價標準Table 1 Water quality evaluation standard of Nemerow Pollution Index method

式中：Fi為第i種污染因子實測濃度與其在某環(huán)境質量標準第j種標準值的比值；Pj為按某環(huán)境要素質量標準中第j種標準計算得出的評分值，內梅羅污染指數(shù)法水質評價標準見表1；Fmax為在某環(huán)境要素質量標準下Fi的最大值；Fave為在某環(huán)境要素質量標準下Fi的算術平均值；i為參評污染因子數(shù)量，i=1,2,···,n；j為某環(huán)境要素質量標準選取的污染因子所對應的質量類別，j=1,2,···,m；Ci為第i種污染因子的實測濃度；Sij為第i種污染因子在某環(huán)境質量標準中第j種標準下的標準值。

2.1.2 整體水質評價

考慮到采樣區(qū)域分為封閉型采煤沉陷區(qū)及開放型采煤沉陷區(qū)2 大類，為了提高水質評價的科學性和準確性，按照沉陷區(qū)類型分別進行水質評價。對于封閉型采煤沉陷區(qū)，考慮到沉陷積水區(qū)的功能、用途不盡相同，分別對標地表水Ⅲ類、Ⅴ類標準進行評價。對于開放型采煤沉陷區(qū)，參照沉陷區(qū)所連通河流的現(xiàn)狀水質狀況進行評價。

1) 開放型

依據(jù)淮南市生態(tài)環(huán)境局發(fā)布的《2021 年淮南市環(huán)境質量狀況公報》[19]，研究區(qū)開放型采煤沉陷區(qū)所涉及的淮河、泥河、架河、西淝河、濟河等水體水質均為Ⅲ類水，故水質評價參考標準取地表水Ⅲ類。評價結果見表2，水質評價等級為優(yōu)良的指標有3 項，分別為Cu、As 和Cr；評價等級為良好的指標有5 項，分別為DO、COD、NH3-N、Hg 和F-；評價等級為較好的指標有2 項，分別為TN、TP；無水質評價等級為較差和差的指標。

表2 采煤沉陷區(qū)整體水質評價結果Table 2 Overall water quality evaluation results of coal mining subsidence area

2) 封閉型

考慮到封閉型沉陷積水區(qū)的功能、用途不盡相同，本著科學客觀評價的態(tài)度，分別對標地表水Ⅲ類、Ⅴ類標準進行封閉型采煤沉陷區(qū)的水質評價。對于淮南礦區(qū)封閉型采煤沉陷區(qū)：對標地表水Ⅲ類標準，水質評價等級為優(yōu)良的指標有2 項，分別為Cu 和Cr；等級為良好的指標有4 項，分別為DO、NH3-N、As、Hg；等級為較好的指標有2 項，分別為COD 和TN；等級為較差的指標有2 項，分別為TP、F-。對標地表水Ⅴ類標準，水質評價等級為優(yōu)良的指標有5 項，分別為NH3-N、As、Hg、Cu、Cr；等級為良好的指標有4 項，分別為DO、COD、TN、TP；等級為較好的指標有1 項，為F-。

2.1.3 分區(qū)水質評價

考慮到研究區(qū)面積較大(1 591 km2)，且各礦開采方式、開采年限、與地表水系連通等情況不盡相同，為進一步探究不同類型礦區(qū)的水質狀況，根據(jù)礦井生產實際、地理位置和后期水質評價需求，將研究區(qū)劃分為老礦區(qū)、潘集片區(qū)、鳳臺潁上片區(qū)3 個區(qū)域，劃分情況如圖2 所示。其中：老礦區(qū)地理處在淮河以南，均為關閉礦井，沉陷區(qū)水域形成時間長，水體環(huán)境較為穩(wěn)定，含李一、李二、謝一、李嘴孜、新莊孜5 對礦井；潘集片區(qū)地理處在淮河以北，開礦時間較短，基本為生產礦井，沉陷區(qū)水域被泥河貫通，含潘一、潘二、潘四東、朱集東、潘三5 對礦井；鳳臺潁上片區(qū)開礦時間最短，開放型沉陷區(qū)水域數(shù)量較多，且多數(shù)沉陷水體被濟河與西淝河貫通，含丁集、顧橋、顧北、張集、謝橋5 對礦井。采煤沉陷區(qū)分區(qū)水質評價結果見表3。

表3 采煤沉陷區(qū)分區(qū)水質評價結果Table 3 Evaluation results of water quality in mining subsidence area

圖2 不同類型礦區(qū)劃分Fig.2 Division map of different types of mining area

1) 開放型

老礦區(qū)：水質評價等級為優(yōu)良的指標有5 項，分別為TP、Cu、As、F-和Cr；等級為良好的指標有5 項，分別為DO、COD、TN、NH3-N、Hg。

潘集片區(qū)：水質評價等級為優(yōu)良的指標有3 項，分別為Cu、As 和Cr；等級為良好的指標有7 項，分別為DO、COD、TN、NH3-N、TP、Hg 和F-。

鳳臺潁上片區(qū)：水質評價等級為優(yōu)良的指標有3 項，分別為Cu、As 和Cr；等級為良好的指標有5 項，分別為DO、COD、NH3-N、Hg 和F-；等級為較好的指標有2 項，分別為TN、TP。

2) 封閉型

(1) 對標地表水Ⅲ類標準

老礦區(qū)：水質評價等級為優(yōu)良的指標有5 項，分別為NH3-N、As、Cu、Cr 和F-；等級為良好的指標有4 項，分別為DO、COD、TN 和Hg；等級為較好的指標有1 項，分別為TN。

潘集片區(qū)：水質評價等級為優(yōu)良的指標有3 項，分別為As、Cu 和Cr；等級為良好的指標有3 項，分別為DO、NH3-N 和Hg；等級為較好的指標有3 項，分別為COD、TN 和TP；等級為較差的指標有1 項，為F-。

鳳臺潁上片區(qū)：水質評價等級為優(yōu)良的指標有3 項，分別為As、Cu 和Cr；等級為良好的指標有3 項，分別為TP、NH3-N 和Hg；等級為較好的指標有3 項，分別為DO、COD 和F-；等級為較差的指標有1 項，為TP。

(2) 對標地表水Ⅴ類標準

老礦區(qū)：水質評價等級為優(yōu)良的指標有8 項，分別為COD、TP、NH3-N、Hg、As、Cu、Cr 和F-；等級為良好的指標有1 項，為TP；等級為較好的指標有1 項，為DO。

潘集片區(qū)：水質評價等級為優(yōu)良的指標有5 項，分別為NH3-N、Hg、As、Cu、Cr；等級為良好的指標有3 項，分別為COD、TP、TN；等級為較好的指標有2 項，為DO 和F-。

鳳臺潁上片區(qū)：水質評價等級為優(yōu)良的指標有6 項，分別為DO、NH3-N、Hg、As、Cu、Cr；等級為良好的指標有4 項，分別為COD、TP、TN 和F-。

2.2 沉陷區(qū)富營養(yǎng)化評價

根據(jù)研究，淮南礦區(qū)沉陷水體主要污染物為營養(yǎng)元素，故采用綜合營養(yǎng)狀態(tài)指數(shù)法進行研究區(qū)水體的富營養(yǎng)化狀態(tài)評價。營養(yǎng)狀態(tài)參數(shù)的計算采用均值方法，選取COD、TN、TP 和Chl(葉綠素)共4 個指標，用綜合營養(yǎng)狀態(tài)指數(shù)進行營養(yǎng)狀態(tài)評價，計算公式如下：

式中：TLI(∑)為綜合營養(yǎng)狀態(tài)指數(shù)；wb為b種參數(shù)的營養(yǎng)狀態(tài)指數(shù)的相關權重；rb為第b種參數(shù)與基準參數(shù)Chl 的相關系數(shù)；c為評價參數(shù)的個數(shù)；TLI(b)為第b種參數(shù)的營養(yǎng)狀態(tài)指數(shù)，分別由以下公式計算：

根據(jù)水體營養(yǎng)物質的污染程度，分成貧營養(yǎng)、中營養(yǎng)、輕度富營養(yǎng)、中度富營養(yǎng)及重度富營養(yǎng)5 種水平[20-23]，沉陷水體富營養(yǎng)化評價結果見表4。計算結果顯示：175 個采樣點的綜合營養(yǎng)指數(shù)范圍為40～90。其中：15 個采樣點評價等級為輕度富營養(yǎng)；90 個采樣點為中度富營養(yǎng)；4 個采樣點為中營養(yǎng)；66 個采樣點為重度富營養(yǎng)?；茨系V區(qū)沉陷區(qū)水域平均營養(yǎng)指數(shù)大小為68.61，屬于中度富營養(yǎng)，說明沉陷區(qū)水體營養(yǎng)水平適中。

表4 水體富營養(yǎng)化程度評價結果Table 4 Water eutrophication assessment results

3 沉陷區(qū)污染因子識別

3.1 內梅羅污染指數(shù)法

根據(jù)內梅羅污染指數(shù)法對研究區(qū)參評因子的計算結果，初步識別了沉陷區(qū)污染因子有2 項，分別為：TP和F-，污染狀況分布如圖3 所示。對于開放型水域，TP、F-指標較差的區(qū)域集中在張集、潘一沉陷區(qū)；對于封閉型水域，TP 指標較差的區(qū)域集中在潘一、潘三、顧橋、顧北、張集沉陷區(qū)，F(xiàn)-污染分布呈現(xiàn)地域特征，淮河以北沉陷區(qū)指標普遍較差，淮河以南沉陷區(qū)指標普遍優(yōu)良。

圖3 TP、F-污染狀況分布Fig.3 TP,F- pollution distribution

將老礦區(qū)、潘集片區(qū)、鳳臺潁上片區(qū)的開放型、封閉型沉陷區(qū)污染指數(shù)進行對比，分析結果如圖4 所示：研究區(qū)所有水域的3 個片區(qū)中，Cr、Cu、Hg、As、NH3-N、TN、COD、DO 指標均較好，整體處于“優(yōu)良-較好”區(qū)間內；開放型水域3 個片區(qū)F-、TP 指標較好，處于“優(yōu)良-較好”區(qū)間內；封閉型水域潘集片區(qū)F-處于“較差”區(qū)間內，鳳臺潁上片區(qū)TP 處于“較差”區(qū)間內；分析3 個片區(qū)污染指數(shù)分布情況，區(qū)域污染程度表現(xiàn)為：鳳臺潁上片區(qū)>潘集片區(qū)>老礦區(qū)。

圖4 開放型與封閉型沉陷區(qū)污染指數(shù)對比Fig.4 Contrast map of pollution index in flood period

3.2 單因子評價法

對標《地表水環(huán)境質量標準》，分析結果見表5：(1) 常規(guī)指標中，DO Ⅰ—Ⅲ類占比66.86%，超過Ⅴ類占比5.71%；NH3-N Ⅰ—Ⅲ類占比89.14%，超過Ⅴ類占比0.57%。COD Ⅰ—Ⅲ類占比39.43%，超過Ⅴ類占比21.71%；TP Ⅰ—Ⅲ類占比71.43%，超過Ⅴ類占比21.13%；TN Ⅰ—Ⅲ類占比12.57%，超過Ⅴ類占比53.14%。(2) 重金屬指標中，As Ⅰ類占比100%；Cu、Cr Ⅰ—Ⅱ類占比100%；Hg Ⅰ—Ⅳ類占比100%。(3) 無機陰離子指標F-Ⅰ—Ⅲ類占比58.86%，超過Ⅴ類占比16.25%。

表5 沉陷區(qū)水質指標單因子法評價結果Table 5 Evaluation results of water quality index in subsidence area by single factor method

單因子法評價結果顯示，常規(guī)指標中DO、NH3-N超過Ⅴ類占比極小，而COD、TP、TN 超過Ⅴ類占比較大；研究區(qū)所有水域均無重金屬超標現(xiàn)象，測試結果極佳；無機陰離子指標F-超過Ⅴ類占比較大。對比發(fā)現(xiàn)，開放型水域有3 個沉陷區(qū)存在超地表水Ⅴ類標準的水樣；封閉型水域有7 個沉陷區(qū)存在超Ⅴ類水樣。開放型采煤沉陷區(qū)水質狀況整體優(yōu)于封閉型采煤沉陷區(qū)。

甲基對硫磷、馬拉硫磷、樂果3 項農藥指標，可以通過食物鏈傷害人體。有機類指標測試結果顯示，所有水樣中均未檢出上述3 項農藥指標，初步說明研究區(qū)域水體中不含上述3 項農藥殘留。

基于單因子評價法分析13 項指標的測試結果可以發(fā)現(xiàn)，淮南礦區(qū)沉陷區(qū)主要污染因子有COD、TP、TN、F-共4 項，各污染因子超地表水Ⅴ類標準情況如圖5 所示。李一、李二、謝一、新莊孜、李嘴孜、潘二、朱集東沉陷區(qū)無污染因子，水質狀況最佳；污染因子數(shù)最多的為潘一(4 個)、潘三(3 個)、顧 橋(3 個)、張集(3 個)沉陷區(qū)，污染因子分布情況見表6。

表6 沉陷區(qū)污染因子分布情況Table 6 Distribution of pollution factors in subsidence area

圖5 各污染因子超地表水Ⅴ類標準所占比例Fig.5 The proportion of each pollution factor exceeding the surface water class V standard

3.3 識別結果對比分析

對比2 種沉陷區(qū)污染因子識別結果發(fā)現(xiàn)：內梅羅污染指數(shù)法識別污染因子2 項，分別為TP 和F-；單因子評價法識別污染因子4 項，分別為COD、TP、TN和F-，TP 和F-為共同識別的污染因子。不難發(fā)現(xiàn)，單因子評價法選擇水質最差的單項指標所屬類別來確定所屬水域的綜合水質類別，適用于評價污染物較為單一或受某一污染物影響較為嚴重的水體；而內梅羅污染指數(shù)法在一定程度兼顧了其他水質較好的參評因子對總體結果的貢獻，適合大范圍、大尺度的區(qū)域性污染因子識別及水質狀況評價。從保護煤礦區(qū)生態(tài)環(huán)境的角度出發(fā)，綜合考慮2 種評價結果的特性，最終識別出污染因子4 項，分別為COD、TP、TN 和F-。

3.4 沉陷區(qū)污染源淺析

前述研究表明，研究區(qū)主要污染因子有4 項，分別為COD、TP、TN、F-。前3 項指標可以表征水體的富營養(yǎng)化程度，值越高代表水體富營養(yǎng)化越嚴重；第4 項指標主要是以氟化氫及其他氟化物的形式出現(xiàn)，水體中氟化物濃度高會使人體骨骼受害。

針對超標因子，結合現(xiàn)場調查發(fā)現(xiàn)，沉陷區(qū)范圍內無明顯排污口，周邊煤礦均設有獨立的污水處理設施，推測受采礦等工業(yè)活動污染的可能性極??；超標點位附近散布耕地和居民生活用地，容易受到農業(yè)生產面源和居民生活點源的污染；其次，也存在受污染因子本底值較高等環(huán)境因素影響的可能性。下一步需采集不同水期的樣品，深入分析研究區(qū)水質狀況及污染因子是否受到水期的影響；同時采集周邊農田、溝渠、居民點水樣以及沉陷區(qū)底泥等樣品，具體分析環(huán)境因素對沉陷區(qū)水質的影響程度，溯清污染源頭。調查同時發(fā)現(xiàn)水體富營養(yǎng)化程度較高的沉陷水體大多被私人承包，存在魚類、貝類養(yǎng)殖的現(xiàn)象，推測這是導致水體富營養(yǎng)化的主要因素。

針對已污染沉陷區(qū)的治理問題，筆者提出以下建議：一是盡快溯源。將現(xiàn)場調查與采樣分析相結合，查明污染因子的超標原因，控制住污染源頭。二是強化管理。建議政府部門加強對沉陷區(qū)的屬地管理，從政策上引導并規(guī)范沉陷區(qū)水域承包商的養(yǎng)殖行為，最大限度減輕水體的富營養(yǎng)化程度。三是加強治理。采用物理法、化學法、生物法等多種手段，進行污染水域的治理工作，使其恢復正常水體功能，滿足周邊居民的生產生活需求。

4 結論

a.采用了無人機技術，采集安徽淮南礦區(qū)15 對礦井、102 片沉陷積水區(qū)、175 個點位的水樣，開展水質測試分析，無人機技術提供了大規(guī)模數(shù)據(jù)量。

b.優(yōu)選了內梅羅污染指數(shù)法，從整體、分區(qū)2 個維度分別對開放型及封閉型沉陷區(qū)進行水質評價，得出了相關離子在不同地表水環(huán)境質量標準下的污染級別，結果顯示開放型沉陷區(qū)水域水質狀況整體優(yōu)于封閉型。

c.綜合營養(yǎng)狀態(tài)指數(shù)法水體富營養(yǎng)化狀態(tài)評價結果表明：研究區(qū)綜合營養(yǎng)指數(shù)范圍為40～90，平均68.61，屬于中度富營養(yǎng)，沉陷區(qū)水體營養(yǎng)水平適中。

d.采用內梅羅污染指數(shù)法、單因子評價法2 種方法識別污染因子，從保護煤礦區(qū)生態(tài)環(huán)境的角度出發(fā)，最終確定研究區(qū)污染因子有COD、TP、TN 和F-。

e.對于污染因子溯源問題，筆者團隊正在圍繞上游河流污染、沉陷區(qū)底泥污染、水期影響3 個方面開展研究，將連續(xù)3 年分豐、枯2 期對研究區(qū)開展大尺度、長系列采樣測試，調查沉陷區(qū)水域周邊潛在污染源，以期盡快溯清污染源頭。