一种矿井水常量离子特征标识及水源判别的方法

技术领域

本发明涉及水文地质及水化学领域，属于矿井水常量离子特征标识及水源判别技术领域。

背景技术

目前，面对煤矿区水文地质条件复杂、矿井水害频发的问题，有必要采用一种简单有效的方法来确定矿井涌水的来源和构成，以用于指导煤矿防治水工作的进行。地下水由于赋存、演化环境的不同，各含水层水物理、化学特征存在差异性，基于此，利用同位素、水温、地下水水位动态、微量元素含量、常量离子浓度等参数，本领域的技术人员做了深入的研究，提出了多种水源确定、判断、识别的方法。其中如何明确各水源含水层水质特征是影响判断方法准确性的基础，为了快速识别矿井突水水样来源构成，有必要提前建立准确的水源含水层水质特征库。

煤矿区地下水常量离子含量相较于其他参数而言易于监测，具有快速、时效性强的特点，是用于判别水样来源的优质选择，但合适的特征离子以及计算方法的选取是影响水源判别准确度的关键。常规方法直接利用水样常量离子浓度进行数学方法计算，忽略了含水层水质可能存在某种离子浓度相似带来的计算误差以及实验室测样带来的测量误差，从而导致计算精度不高。

发明内容

为了克服上述现有技术的不足之处，本发明提供一种矿井水常量离子特征标识及水源判别的方法，公开了一种建立矿井含水层水样特征离子库的方法以及利用合理特征离子定量计算矿井水来源构成比例的方法，能够提高计算矿井水来源构成比例的精度，本发明也是国家重点研发计划项目“煤矿区场地地下水污染防控材料与技术(2019YFC1805400)”的成果之一。

本发明是通过如下技术方案实现的：一种矿井水常量离子特征标识及水源判别的方法，包括：

步骤一、获取矿井主要充水含水层的水质数据，对水样的6大常量离子的浓度值进行统计整合,所述6大常量离子包括K

步骤二、排除各含水层异常水样点，明确各水源常量离子的分界阈值；

对统计的水样数据进行系统聚类分析，将6大常量离子的浓度值作为变量进行聚类；聚类方法采用“系统聚类法”，系统聚类法是先将聚类的样本或变量各自看成一群，然后确定类与类间的相似统计量，并选择最接近的两类或若干个类合并成一个新类，计算新类与其他各类间的相似性统计量，再选择最接近的两群或若干群合并成一个新类，直到所有的样本或变量都合并成一类为止；距离采用平方欧式距离，将距离小于5的水样分为一类，根据聚类结果，排除水样异常点(即水化学特征偏离总体特征的水样点)，明确各水源含水层水质特征界限；

步骤三、确定各水源常量离子的标准浓度，根据数据变异系数，数据变异系数是指原始数据标准差与原始数据平均数的比，筛选出各水源含水层差异性明显的常量离子作为特征离子；

根据步骤二所筛选水样，整理归纳各水源水化学常量离子浓度取值范围，取均值作为其常量离子的标准浓度；依据变异系数这一指标剔除各水源变异系数较大的常量离子，选取变异系数小于10％的常量离子作为用于计算的特征离子；

步骤四、建立二级模糊综合评判数学模型，判别矿井水的水源，并定量计算构成比例。

进一步，所述步骤二中，利用SPSS软件对统计的水样数据进行系统聚类分析。

进一步，所述步骤四中，采用偏标加权法确定代表因素的权重，采用专家推理法确定隶属度，对矿井水源进行模糊综合判别。

进一步的，所述步骤四中，建立二级模糊综合评判数学模型为：

首先，建立特征离子的因素集：

U＝{X，Y，Z，......}；

然后，建立评价集

V＝{水源含水层1、水源含水层2、水源含水层3}

三个含水层分别用I、II、Ⅲ来表示；

再建立权重向量A＝{a

采用偏标加权法确定权重；

式中：a

根据偏标加权法确定的待判定水样各评价因素的权重向量为

A＝{a

隶属度的确定采用专家推理法，采用梯形模糊分布的隶属度函数类型，不受离子种类限制，适用于不同环境下的矿井；当特征离子X在I中浓度最小，在II中浓度中等，在Ⅲ中浓度最高时，X的隶属函数模型如下：

式中：r

同理，分别计算出其余特征离子作为影响因素的评判结果，进而得出模糊评价矩阵R；

利用偏标加权法对矿井水样的评价结果计算为：B＝AR；

根据计算结果，得出待判定水样来源。

本发明的有益效果是：本发明方法中，首先利用聚类分析方法明确各水源含水层常量离子标准浓度，然后通过变异系数这一指标选取差异性明显的特征离子作为计算参数，进一步消除了计算误差及测量误差，最后建立二级模糊综合评判数学模型，定量计算了矿井水来源构成比例。

本发明一方面建立了长久可用的矿区含水层水质特征离子库，另一方面提升了矿井水来源构成比例计算的精度，为矿井突水及涌水量构成识别提供了一种新的实用方法。

附图说明

下面根据附图和实施例对本发明进一步说明。

图1是本发明的建立矿井含水层水样特征离子库及定量计算矿井水来源构成比例方法流程图。

图2是本发明实施例提供的利用SPSS软件所制的水样聚类结果示意图。

具体实施方式

下面将结合说明书附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。以下对至少一个示例性实施例的描述实际上仅仅是说明性的，决不作为对本发明及其应用或使用的任何限制。基于本发明中的实施例，本领域普通技术人员在没有开展创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

对于本领域技术人员已知的技术、方法和设备可能不作详细讨论，但在适当情况下，所述技术、方法和设备应当被视为授权说明书的一部分。

一种矿井水常量离子特征标识及水源判别的方法，包括：

步骤一、获取矿井主要充水含水层的水质数据，对水样的6大常量离子的浓度值进行统计整合,所述6大常量离子包括K

步骤二、排除各含水层异常水样点，明确各水源常量离子的分界阈值；

对统计的水样数据进行系统聚类分析，将6大常量离子的浓度值作为变量进行聚类；聚类方法采用“系统聚类法”(这种方法是先将聚类的样本或变量各自看成一群，然后确定类与类间的相似统计量，并选择最接近的两类或若干个类合并成一个新类，计算新类与其他各类间的相似性统计量，再选择最接近的两群或若干群合并成一个新类，直到所有的样本或变量都合并成一类为止)，距离采用平方欧式距离，将距离小于5的水样分为一类，根据聚类结果，排除水样异常点(即水化学特征偏离总体特征的水样点)，明确各水源含水层水质特征界限；

步骤三、确定各水源常量离子的标准浓度，根据数据变异系数(原始数据标准差与原始数据平均数的比)筛选出各水源含水层差异性明显的常量离子作为特征离子；

根据步骤二所筛选水样，整理归纳各水源水化学常量离子浓度取值范围，取均值作为其常量离子的标准浓度；依据变异系数这一指标剔除各水源变异系数较大的常量离子，选取变异系数小于10％的常量离子作为用于计算的特征离子(此句不清楚是将变异系数小于10％的离子留下还是排除)；

步骤四、建立二级模糊综合评判数学模型，判别矿井水的水源，并定量计算构成比例。

如图1所示，下面将结合实施例以更加清晰、明确表述本发明的技术方案。

步骤一、统计整合矿井主要含水层水样水质资料。

某矿矿井主要充水水源为第四系底砾层水(底含水)、煤层顶底板砂岩含水层水(砂岩水)、煤层底板太原组L

步骤二、利用聚类分析方法，排除各水源含水层异样水样，明确各水源含水层水质特征的明显分界。

利用SPSS软件对统计的水样数据进行系统聚类分析，将6大特征离子浓度值作为变量进行聚类。聚类方法采用“沃德联接方法”，距离采用平方欧式距离，将距离小于5的水样分为一类。根据聚类结果，排除水样异常点，明确各水源含水层水质特征界限。

如图2所示，18组水样被分为了A(1、2、3、4、5、6)、B(7、8、9、10、11、12、13、14)、C(15、16、17、18)三类。A组水样为四灰水，B组为底含水，C组为砂岩水。18组水样聚类结果与含水层分类结果一致，无特殊需要排除水样，说明三种主要含水层水质特征界限具有明显分别，具备用于进一步判别水源构成的条件。将A、B、C三组的离子特征浓度范围进行统合分析。

步骤三、建立各水源含水层水样水质离子标准浓度，根据数据变异系数筛选出各水源含水层差异性明显的特征离子用于计算。

根据步骤二所筛选水样，整理归纳各水源水质特征离子浓度取值范围，取均值作为其水质特征标准值。依据变异系数这一指标剔除各水源含水层变异系数小于10％的常量离子作为用于计算的特征离子。

表1根据本发明实施例建立的水源含水层标准值及待判定水样常量离子浓度表

根据表1所示结果，三类水源含水层水样特征总体差异性较为明显，但部分常量离子浓度仍然较为接近。除此以外，各含水层水样在进行测定的时候可能存在实验室测量误差，为了进一步提升计算精度，消除测量及计算方法误差，本次依据变异系数这一指标剔除各水源含水层差异性小于10％的常量离子作为用于计算的特征离子。最终选取K

步骤四、建立二级模糊综合评判数学模型，定量计算矿井混合水样来源及构成比例。

首先，建立因素集

然后，建立评价集

V＝{第四系底含、顶底板砂岩水、L

三个含水层分别用I、II、Ⅲ来表示；

再建立权重向量A＝{a

采用偏标加权法确定权重。

式中：a

表2根据偏标加权法确定的待判定水样各评价因素的权重表

隶属度的确定采用专家推理法，采用梯形模糊分布的隶属度函数类型，不受离子种类限制，适用于不同环境下的矿井。以K

式中：r

同理，分别计算出Ca

利用偏标加权法对矿井水样的评价结果计算为：B＝AR＝(0.342，0.64，0.018)。

根据计算结果，可以得出本次矿井水的待判定水样来源为34.2％的底含水，64％的砂岩水以及1.8％的四灰水，经查验与矿井生产实际涌水情况相符。

以上所述仅为本发明的示例性实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。