采動影響下皖北恒源煤礦八含水化學(xué)特征研究

施亞麗 吳基文 翟曉榮 王廣濤 洪荒 畢堯山

摘 要：為研究礦井開采活動對含水層的擾動影響，以皖北恒源煤礦八含為研究對象，收集了1988a～2019a共61個八含水樣的水質(zhì)數(shù)據(jù)，運用聚類分析方法將所采水樣分成3個不同開采階段的數(shù)據(jù)組，通過分析不同開采階段的八含水化學(xué)特征，從時間維度上研究了采礦活動對八含水質(zhì)的擾動影響，并通過建立不同開采階段八含水質(zhì)和太灰水質(zhì)的偏最小二乘判別（PLS-DA）模型討論了八含水與太灰水之間的水力聯(lián)系。結(jié)果表明：（1）隨著開采活動的進行，八含水的水化學(xué)特征發(fā)生了明顯變化;（2）對比不同開采階段的八含水質(zhì)與太灰水質(zhì)，顯示開采初期八含水水質(zhì)與太灰水水質(zhì)迥異，但隨著采礦活動的進行，開采Ⅲ階段的八含水水質(zhì)趨近于太灰水水質(zhì);（3）不同開采階段八含水質(zhì)具有明顯的聚類效應(yīng)，部分開采Ⅲ階段的八含水樣水質(zhì)呈現(xiàn)出太灰水與八含水的混合水水質(zhì)特征，驗證了八含水在采動影響下與太灰水發(fā)生水力聯(lián)系的結(jié)果。研究成果可為礦井水害防治以及水資源的開發(fā)利用提供科學(xué)依據(jù)。

關(guān)鍵詞：水化學(xué)特征變化;水力聯(lián)系;聚類分析;PLS-DA;采動影響

中圖分類號： TD745.21? 文獻標志碼：A

文章編號：1672-1098（2021）01-0031-10

收稿日期：2020-09-29

基金項目：國家自然科學(xué)基金資助項目（41272278）;安徽高校科研平臺創(chuàng)新團隊建設(shè)基金資助項目（2016-2018-24）

作者簡介：施亞麗（1995-），女，安徽合肥人，在讀碩士，研究方向：水文地質(zhì)與工程地質(zhì)。

Study on? Hydrochemical Characteristics of the Eighth? Aquifer of? Hengyuan Coal Mine in Northern Anhui Province under the? Influence of Mining

SHI Yali1，WU Jiwen1，ZHAI Xiaorong1，WANG Guangtao1 HONG Huang2，BI Yaoshan1

（1. School of Earth and Environment，Anhui University of Science and Technology，Huainan? Anhui? 232001，China;2. Hengyuan Coal Mine of Anhui Hengyuan Coal Power Co.，Ltd， Huaibei Anhui? 235000， China）

Abstract：In order to study the influence of mining activities on aquifer disturbance， with the eighth aquifer of Hengyuan Coal Mine in Northern Anhui Province as the research object， the water quality data were collected of 61 water samples at the eighth aquifer from 1988 to 2019， which were divided into three data groups of different mining stages with the cluster analysis method. By analyzing the hydrochemical characteristics of the eighth aquifer in three different mining stages， the disturbance effect of mining activities on the water quality was studied in the time dimension. And by establishing the PLS-DA recognition model of water quality of the eighth aquifer and the Taiyuan formation limestone aquifer in different mining stages， the hydraulic connection was discussed between the groundwater of the eighth aquifer and the Taiyuan formation limestone aquifer. The results show the following facts：（1）With the development of mining activities， the hydrochemical characteristics of the eighth aquifer change obviously;（2）Comparing the water quality of the eighth aquifer and the Taiyuan formation limestone aquifer in different mining stages reveals that the water quality of the eighth aquifer water is different from that of the Taiyuan formation limestone aquifer in the early stage of mining. With the mining activities going on， the water quality of the eighth aquifer in the third stage of mining is similar to that of the Taiyuan formation limestone aquifer;（3）There is obvious clustering effect in different mining stages and the water quality of the eighth aquifer water samples in the third stage of mining shows the characteristics of mixed water of the Taiyuan formation limestone aquifer and the eighth aquifer water， which proves the hydraulic connection between the eighth aquifer water and the Taiyuan formation limestone aquifer under the influence of mining. The research results will provide scientific basis for mine water disaster prevention and development and utilization of water resources.

Key words：change of hydrochemical characteristics; hydraulic connection; cluster analysis; PLS-DA; mining influence

礦井水害問題嚴重威脅著煤礦的安全生產(chǎn)[1]。隨著淺部煤炭資源的枯竭，逐漸轉(zhuǎn)為煤炭資源的深部開采階段。由于深部開采導(dǎo)致深層地下水受到一定程度的破壞和擾動，使地下水系統(tǒng)的水動力場和水化學(xué)場發(fā)生改變，水文地質(zhì)條件變得更加復(fù)雜，對礦井的正常生產(chǎn)造成了更大的安全威脅[2-3]。因此，開展采動影響下深部含水層的水化學(xué)特征研究以進一步掌握深部地下水水化學(xué)的演化規(guī)律，對礦井水害防治以及水資源的開發(fā)利用具有重要的意義[4-6]。以往對于充水含水層的水化學(xué)特征的研究大多局限于從含水層的靜態(tài)水化學(xué)場分析[7]363，文獻[8]利用傳統(tǒng)圖示及統(tǒng)計方法分析了太原組灰?guī)r水徑流特征及其水化學(xué)變化;文獻[9]采用多元統(tǒng)計、水化學(xué)分析與地下水系統(tǒng)理論相結(jié)合的方法，對巖溶含水系統(tǒng)的水化學(xué)組分特征及其控制因素進行了探討;文獻[10]等采用因子分析、聚類分析和判別分析等方法對皖北桃源礦深部含水層地下水地球化學(xué)特征進行了深入研究。一些學(xué)者也從充水含水層的水巖作用機制出發(fā)，考慮了采礦活動對含水層的擾動影響，如文獻[7]364等運用主成分分析法對任樓井田的充水含水層的水化學(xué)特征進行了分析，文獻[11]等通過主成分分析研究了采動影響下水化學(xué)演化過程中的水-巖作用機制，但此類研究相對較少，對于分析采礦前后全過程中充水含水層的水化學(xué)演化規(guī)律的研究相對更少。

鑒于此，本文以皖北恒源煤礦深部含水層八含（該礦井主采煤層6煤的頂?shù)装迳皫r裂隙含水層）為研究對象，通過收集開采前后不同時期的水質(zhì)資料，嘗試運用聚類分析方法從時間維度上分析不同開采階段下八含的水化學(xué)特征變化，初步研究采礦活動對八含水質(zhì)的擾動影響，并探討在采動影響下八含可能與太灰含水層發(fā)生的水力聯(lián)系。

1 研究區(qū)概況

安徽恒源煤電股份有限公司煤礦（簡稱恒源煤礦）隸屬于皖北煤電集團公司，位于淮北煤田濉肖礦區(qū)。礦井內(nèi)含煤地層為二疊系下統(tǒng)山西組和下石盒子組，主采煤層為4、6煤層，分3個水平開采：-400m、-600m、-800m。目前礦井生產(chǎn)水平為-400m和-600m。礦井內(nèi)褶皺構(gòu)造比較復(fù)雜，斷層不但數(shù)量多、密度大，而且發(fā)育方向各異，性質(zhì)不同，相互改造和制約。礦井為新生界松散層覆蓋下的裂隙充水礦床。根據(jù)含水層賦存介質(zhì)特征自上而下劃分為新生界松散層孔隙含水層，二疊系煤系砂巖裂隙含水層，太原組石灰?guī)r巖溶裂隙含水層（簡稱“太灰”），奧陶系石灰?guī)r巖溶裂隙含水層（簡稱“奧灰”）。其中，二疊系煤系砂巖裂隙含水層，包括第五含水層（簡稱“五含”）、3煤上第六含水層（簡稱“六含”）、4煤上下第七含水層（簡稱“七含”）、6煤上下第八含水層（簡稱“八含”），七含、八含是主采煤層4、6煤的直接充水含水層。根據(jù)礦井生產(chǎn)資料，太灰水是6煤開采的礦坑充水的隱患之一，奧灰遠離主采煤層，在正常情況下對礦坑充水無明顯影響，井田水文地質(zhì)剖面圖如圖1所示。

2 數(shù)據(jù)來源與處理

2.1 水樣數(shù)據(jù)收集

本次研究收集并整理了恒源煤礦自投產(chǎn)以來的八含水樣資料共計61個，收集的水樣點分布如圖1所示，收集水樣的時間跨度為1988a7月～2019a 4月，收集水樣的時間分布均勻，各水樣的水化學(xué)成分測試數(shù)據(jù)如表1所示?？紤]到測量誤差等因素，需對61個水化學(xué)數(shù)據(jù)的有效性進行驗證，本文采用陰陽離子平衡檢驗（電荷守恒檢驗）方法對水樣數(shù)據(jù)進行驗證，計算公式為

E=∑mc-∑ma∑mc+∑ma×100（1）

式中：E為相對誤差，%;mc、ma分別為陽、陰離子的毫克當量濃度，meq·L-1。

經(jīng)計算，61個水樣數(shù)據(jù)的陰陽離子平衡檢驗的相對誤差的絕對值均小于5%，表明61個水樣均為有效水樣數(shù)據(jù)。

2.2 水樣的系統(tǒng)聚類分析處理

系統(tǒng)聚類分析的基本思想是：首先將每個樣本單獨看成一類，在規(guī)定類間距離的條件下，選擇距離最小的一對合并成一個新類，計算新類和其他類的距離;再將距離最小的兩類合并。因此每次會減少一類，直至所有樣本都成為一類，這樣使得同一類別的數(shù)據(jù)的相似度盡可能地高，不同類別間數(shù)據(jù)的相似度盡可能低，結(jié)果可用1張譜系圖來表示[12]。一般包括R型聚類分析和Q型聚類分析。R型聚類分析與Q型聚類分析的區(qū)別在于分類對象的不同，R型聚類分析為的對象為變量，而Q型聚類分析的對象為樣本[13-14]。本文利用離差平方和法對61個原始水樣進行Q型聚類分析，用平方歐式距離為度量標準，確定各水樣之間的親疏關(guān)系，得到原始水樣的譜系圖，如圖2所示。

一般來說，埋藏條件類似的同一含水層在未受到擾動影響的情況下，在鄰近的短暫時間內(nèi)只通過水-巖作用而發(fā)生的水化學(xué)特征的變化不大。由于圖中水樣均為按時間先后順序編號，根據(jù)圖中所示各水樣水質(zhì)間的親疏關(guān)系，同時考慮各水樣在時間上的鄰近程度，將42號水樣作為異常數(shù)據(jù)予以剔除。重在分析八含水質(zhì)在時間維度上受采礦活動影響前后的變化過程，基于Q型聚類結(jié)果與水樣在時間上的對應(yīng)關(guān)系，將剩下的60個水樣數(shù)據(jù)劃分為3類，根據(jù)研究區(qū)內(nèi)采礦活動的進行時間，將水樣編號為1～29的29個八含水樣歸為開采Ⅰ階段的水樣（時間跨度為1988.07～2002.12），將水樣編號為30～43的13個八含水樣（剔除42號水樣后）歸為開采Ⅱ階段的水樣（時間跨度為2003.10～2009.05），將水樣編號為44～61的18個八含水樣歸為開采Ⅲ階段的水樣（時間跨度為2009.06～2019.04），如表1所示。

3 不同開采階段八含水化學(xué)特征變化

基于聚類分析結(jié)果，可知本次收集的八含水樣隨著開采階段的不同顯示出不同的水化學(xué)特征。以下分別分析不同開采階段八含的水化學(xué)特征。

3.1 開采Ⅰ階段八含水化學(xué)特征

對開采Ⅰ階段的29個水樣數(shù)據(jù)進行統(tǒng)計分析，如表2所示?？?cè)芙庑怨腆w（TDS）是水中溶解組分的總量，是常規(guī)離子在水中積累的綜合反映，一般來說，在地下水循環(huán)系統(tǒng)中，TDS大，水流流速緩慢，水動力條件弱，地下水基本處于停滯狀態(tài)[15]。該階段八含的TDS值大，在3 133.35～3 766.98mg/L，均值為3 363.36mg/L，表明此階段八含處于相對封閉的條件，與其它含水層無水力聯(lián)系。該階段陰離子以SO2-4為主，約在70%～80%之間，陽離子以Na++K+為主，約在80%以上，各陰離子含量為SO2-4>HCO3->Cl-，各陽離子含量為Na++K+>Ca2+>Mg2+，其中該階段Ca2+、Mg2+的變異系數(shù)較大（超過0.5）分別為0.83、0.88，即Ca2+、Mg2+離子的離散程度較大，表明其在水化學(xué)組成中不穩(wěn)定。該階段八含水主要為SO4·Na類型。開采Ⅰ階段為礦井開采初期，主要開采淺部一水平，對八含基本無擾動影響或影響較小。本階段八含的水化學(xué)特征基本顯示為其原生狀態(tài)。

3.2 開采Ⅱ階段八含水化學(xué)特征

對該階段的13個水樣數(shù)據(jù)進行統(tǒng)計分析， 如表2所示。 該階段八含的TDS值在2 803.57～ 3 097.39mg/L，均值為2 923.79mg/L，與開采Ⅰ階段相比，該階段TDS值降低439.57mg/L，表明該階段八含水受到補給較充分，即與其它含水層存在一定的水力聯(lián)系。該階段八含水的陰離子以SO2-4為主，SO2-4離子約60%～80%之間，與開采Ⅰ階段相比，SO2-4離子含量降低，Cl-離子含量略微降低，HCO-3離子含量略微增大;陽離子以Na++K+為主，與開采Ⅰ階段相比，Na++K+、Ca2+、Mg2+離子含量均表現(xiàn)有略微降低;各陰離子含量為SO2-4>HCO-3>Cl-，各陽離子含量為Na++K+>Ca2+>Mg2+，其中該階段Ca2+、Mg2+的變異系數(shù)較大分別為0.59、0.97，與開采Ⅰ階段相比Ca2+的離散程度相對變小，Mg2+的離散程度較大。該階段八含水也主要為SO4·Na類型。與開采Ⅰ階段相比，該階段八含水質(zhì)發(fā)生了一定程度的變化，表明采礦活動對八含造成了擾動影響，但總體影響不大。

3.3 開采Ⅲ階段八含水化學(xué)特征

對該階段的18個水樣數(shù)據(jù)進行統(tǒng)計分析， 如表2所示。該階段八含的TDS值在1 941.37～ 2 806.11mg/L，均值為2 440.74mg/L，在前一階段的基礎(chǔ)上繼續(xù)降低，與開采Ⅰ階段相比，該階段TDS值降低922.62mg/L，表明該階段八含水受到充分補給，與其它含水層存在較強的水力聯(lián)系。該階段八含水的陰離子以SO2-4為主，與開采Ⅰ階段、開采Ⅱ階段相比，SO2-4離子含量進一步降低，Cl-離子含量進一步降低，HCO-3離子含量降低;陽離子以Na++K+為主，與開采Ⅰ階段、開采Ⅱ階段相比，Na++K+離子含量降低明顯，Ca2+、Mg2+離子含量均明顯增大，含量約為開采Ⅱ階段八含水的4倍，也均超過開采Ⅰ階段八含中Ca2+、Mg2+含量，其中Ca2+含量約為開采Ⅰ階段八含水的2倍。該階段各陰離子含量表現(xiàn)為SO2-4>HCO-3>Cl-，各陽離子含量表現(xiàn)為Na++K+>Ca2+>Mg2+，其中該階段Ca2+、Mg2+、HCO-3的變異系數(shù)較大分別為1.48、1.16、0.57，與開采Ⅰ階段、開采Ⅱ階段相比，Ca2+、Mg2+、HCO-3的離散程度明顯變大。該階段八含水主要為SO4·Na、SO4·Ca類型。與開采Ⅰ階段相比，該階段八含水質(zhì)發(fā)生了較大的變化，表明采礦活動對八含造成了明顯的擾動影響。

4 討論

通過上述分析，不同開采階段八含水質(zhì)發(fā)生了明顯的變化，主要表現(xiàn)為開采活動前后八含水的TDS值降低明顯，Ca2+離子含量增加明顯，Na++K+、Cl-、SO2-4含量不斷降低，八含水類型由SO4·Na類轉(zhuǎn)變?yōu)镾O4·Na、SO4·Ca類型，表明隨采礦活動的進行對八含造成了擾動影響，導(dǎo)致八含水與其他含水層產(chǎn)生了水力聯(lián)系。根據(jù)礦井地質(zhì)與水文地質(zhì)條件以及現(xiàn)場采掘資料，富水性較強的太灰含水層可能是補給八含的充水含水層。為進一步探討八含與太灰之間是否存在水力聯(lián)系，本次研究收集了2018a的10個太灰水樣資料，對比兩含水層的水質(zhì)特征。

4.1 太灰水化學(xué)特征分析

對10個太灰水樣數(shù)據(jù)進行統(tǒng)計分析，如表3所示。由表3可知，太灰水的TDS值在1 634.48～2 715.65mg/L，均值為2 163.96mg/L;太灰水的陰離子以SO2-4為主，陽離子以Na++K+為主;各陰離子含量表現(xiàn)為SO2-4>HCO-3>Cl-，各陽離子含量表現(xiàn)為Na++K+>Ca2+>Mg2+。太灰水主要為SO4·Na類型。以不同開采階段八含水和太灰水各離子成分含量的均值為背景值，分析兩含水層的水質(zhì)特征，如表4所示。由表4可知，開采初期八含水水質(zhì)與太灰水水質(zhì)迥異，但隨著采礦活動的進行，八含水水質(zhì)趨近于太灰水水質(zhì)，而且太灰水的TDS值小于八含水，表明太灰水存在對八含水的補給。

4.2 基于PLS-DA模型的含水層水力聯(lián)系驗證

為進一步驗證八含水與太灰水存在水力聯(lián)系，將不同開采階段八含水樣的水質(zhì)與太灰水樣水質(zhì)數(shù)據(jù)進行進行偏最小二乘判別分析（PLS-DA）。PLS-DA是一種與主成分回歸相關(guān)的有監(jiān)督的判別分析統(tǒng)計方法[16]。在分析過程中可以根據(jù)設(shè)定的組別進行計算，可在盡可能多地提取自變量信息的基礎(chǔ)上，保證提取成分與因變量間最大的相關(guān)性。PLS-DA將現(xiàn)代分析技術(shù)得到的大量數(shù)據(jù)與多元統(tǒng)計方法相結(jié)合，具有高效的鑒別能力，目前在樣品分類與識別方面發(fā)揮了重要作用[17]。本文將開采Ⅰ階段八含水樣、開采Ⅱ階段八含水樣、開采Ⅲ階段八含水樣、太灰水樣的水質(zhì)資料依次分為第1、第2、第3、第4類，通過建立不同開采階段八含水質(zhì)和太灰水質(zhì)的PLS-DA識別模型，最大化提取不同類別之間的水質(zhì)信息差異，進而實現(xiàn)不同類別水樣的識別，最后根據(jù)識別結(jié)果驗證八含與太灰之間的水力聯(lián)系。

圖3是本次建模得到的PLS-DA模型得分圖，圖3中t[1]、t[2]分別表示為在主成分PC1、PC2投影上的得分。如果兩個水樣點在得分圖上的距離越近，表明這兩個水樣的水質(zhì)特征越相似。由圖3可知，不同開采階段下的八含水水樣具有明顯的聚類效應(yīng)，開采Ⅰ階段八含水樣、開采Ⅱ階段八含水樣、開采Ⅲ階段八含水樣根據(jù)水質(zhì)的差異各自聚成一類，表明該模型識別效果較好，而部分開采Ⅲ階段的八含水樣與太灰水樣存在混合，無法完全區(qū)分。從八含水水質(zhì)變化的角度考慮，不同開采階段下八含水質(zhì)發(fā)生了明顯變化（如圖3中箭頭所指方向變化），且在采動影響下，八含水受到了太灰水的水源補給，導(dǎo)致開采Ⅲ階段八含水的水質(zhì)呈現(xiàn)出太灰水與八含水的混合水水質(zhì)特征。

由于礦區(qū)地下水化學(xué)環(huán)境復(fù)雜，本文主要研究了采動影響后的地下含水層間地下水的混合作用。根據(jù)文獻[7]365，受礦井地質(zhì)與水文地質(zhì)條件與采礦活動的影響，黃鐵礦氧化，碳酸鹽、硫酸鹽溶解，陽離子交替吸附以及脫硫酸作用等水化學(xué)成分的形成作用對礦井含水層的水質(zhì)演化特征也可能存在一定影響。本研究區(qū)內(nèi)其他水化學(xué)成分的形成作用對八含水的水質(zhì)變化特征的影響，有待后續(xù)進一步研究探討。

5 結(jié)論

本文通過收集恒源煤礦1988～2019a共61個八含水樣的水質(zhì)數(shù)據(jù)，運用聚類分析方法將所收集水樣分成3個不同開采階段的水樣數(shù)據(jù)，分別分析了3個不同開采階段八含水的水質(zhì)特征，并討論了八含水與太灰水之間的水力聯(lián)系，主要得到以下結(jié)論：

（1）基于對61個有效原始水樣進行的Q型聚類分析結(jié)果，剔除了異常水樣數(shù)據(jù)的影響，并從時間維度上將水樣數(shù)據(jù)劃分為3個階段的水樣數(shù)據(jù)，依次為開采Ⅰ階段（1988.07～2002.12）、開采Ⅱ階段（2003.10～2009.05）、開采Ⅲ階段（2009.06～2019.04）;

（2）分別分析了不同開采階段八含水的水化學(xué)特征，結(jié)果表明不同開采階段八含水的水化學(xué)特征發(fā)生了變化。礦井生產(chǎn)早期對八含基本無擾動影響或影響較小;與開采Ⅰ階段相比，開采Ⅱ階段采礦活動對八含造成了擾動影響，但總體影響不大;開采Ⅲ階段八含水質(zhì)發(fā)生了較大的變化，采礦活動對八含造成了明顯的擾動影響;

（3）通過對比不同開采階段的八含水質(zhì)與太灰水質(zhì)，結(jié)果顯示開采初期八含水水質(zhì)與太灰水水質(zhì)迥異，但隨著采礦活動的進行，八含水水質(zhì)趨近于太灰水水質(zhì);

（4）建立不同開采階段下八含水質(zhì)和太灰水質(zhì)的PLS-DA識別模型，結(jié)果顯示不同開采階段八含水質(zhì)具有明顯的聚類效應(yīng)，而部分開采Ⅲ階段的八含水樣與太灰水樣存在混合，其水質(zhì)呈現(xiàn)出太灰水與八含水的混合水水質(zhì)特征，驗證了太灰水與八含水之間存在水力聯(lián)系。

參考文獻：

[1] 陳紅影. 我國礦井水害的類型劃分與水文結(jié)構(gòu)模式研究[D].徐州：中國礦業(yè)大學(xué)，2019.

[2] 畢堯山. 大屯徐莊煤礦7煤頂板涌（突）水特征及充水含水層富水性評價[D].淮南：安徽理工大學(xué)，2019.

[3] 劉鑫，陳陸望，林曼利，等.采動影響下礦井突水水源Fisher判別與地下水補給關(guān)系反演[J].水文地質(zhì)工程地質(zhì)，2013，40（4）：36-43.

[4] 尹子悅，林青，徐紹輝.青島市大沽河流域地下水水化學(xué)時空演化及影響因素分析[J].地質(zhì)論評，2018，64（4）：1 030-1 043.

[5] 柳鳳霞，史紫薇，錢會，等.銀川地區(qū)地下水水化學(xué)特征演化規(guī)律及水質(zhì)評價[J].環(huán)境化學(xué)，2019，38（9）：2 055-2 066.

[6] 許冬清. 宿縣礦區(qū)地下水化學(xué)演化特征與控制因素研究[D].合肥：合肥工業(yè)大學(xué)，2017.

[7] 陳陸望，劉鑫，殷曉曦，等.采動影響下井田主要充水含水層水化學(xué)環(huán)境演化分析[J].煤炭學(xué)報，2012，37（S2）：362-367.

[8] 陳松，桂和榮.宿縣礦區(qū)太原組灰?guī)r水地球化學(xué)特征及地質(zhì)意義[J].水文地質(zhì)工程地質(zhì)，2016，43（5）：33-41.

[9] 楊婷婷，許光泉，余世滔，等.煤層下部太原組巖溶水化學(xué)組分特征及其成因分析[J].水文地質(zhì)工程地質(zhì)，2019，46（2）：100-108.

[10] 孫林華，桂和榮.皖北桃源礦深部含水層地下水地球化學(xué)數(shù)理統(tǒng)計分析[J].煤炭學(xué)報，2013，38（S2）：442-447.

[11] 殷曉曦，陳陸望，謝文蘋，等.采動影響下礦區(qū)地下水主要水-巖作用與水化學(xué)演化規(guī)律[J].水文地質(zhì)工程地質(zhì)，2017，44（5）：33-39.

[12] 張好，姚多喜，魯海峰，等.主成分分析與Bayes判別法在突水水源判別中的應(yīng)用[J].煤田地質(zhì)與勘探，2017，45（5）：87-93.

[13] 曹雪春，錢家忠，孫興平.煤礦地下水系統(tǒng)水質(zhì)分類判別的多元統(tǒng)計組合模型——以顧橋煤礦為例[J].煤炭學(xué)報，2010，35（S1）：141-144.

[14] 耿建軍.基于系統(tǒng)聚類分析的煤礦突水水源識別技術(shù)——以潞安礦區(qū)王莊煤礦為例[J].中國煤炭地質(zhì)，2019，31（8）：54-59.

[15] 桂和榮. 皖北礦區(qū)地下水水文地球化學(xué)特征及判別模式研究[D].合肥：中國科學(xué)技術(shù)大學(xué)，2005.

[16] 牛麗影，李大婧，劉春泉，等.鮮食玉米中游離糖和游離氨基酸含量差異的多元統(tǒng)計分析[J].江蘇農(nóng)業(yè)學(xué)報，2020，36（2）：463-470.

[17] 佘僧，李熠，宋洪波，等.低聚糖和多酚結(jié)合偏最小二乘判別分析鑒別油菜蜜產(chǎn)地[J].食品科學(xué)，2019，40（12）：290-295.

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