基于K均值聚类地下水水化学特征水文地质单元划分分析.docVIP

基于K均值聚类地下水水化学特征水文地质单元划分分析 　　[摘要]文章结合自身工作经历，首先介绍了传统水文地质单元划分的主要影响因素。接着分析了k均值聚类的基本概念和计算方法及步骤。最后联系平顶山煤田水文地质单元划分的实际情况，探讨基于k均值聚类地下水水化学特征的水文地质单元划分应用效果。以期为水文地质单元划分提供新的理论计算依据。 　　[关键词]k均值聚类 水化学特征 水文地质单元 　　[中图分类号] P332.7 [文献码] B [文章编号] 1000-405X（2015）-3-225-2 　　1引言 　　水文地质单元划分是根据水文地质条件的差异性将所要研究的对象分成若干个小的区域，然后将研究对象进行分类，以期找到研究对象的某些规律的手段。一般而言，水文地质条件的差异性主要体现在研究对象的地形地貌、地质构造、地层条件以及气象和水文等方面。然而，在实际操作中，由于地下条件的复杂性，我们根本无法准确了解整个区域的地质和水文条件。如局部区域可能存在小断层和小溶洞等可以大大改变整个区域地下水的补给、径流以及排泄条件的类外因素。但是从另一个角度来讲，由于每个水文地质单元一般都具有统一的径流、补给、排泄条件。 　　2传统水文地质单元划分方法的主要影响因素 　　2.1地形地貌的影响 　　所要研究区域的地形地貌虽不是影响水文地质单元划分的主要因素。但是由于该区域的湖泊、河流等地表水可能成为地下含水层的补给水源，因此也可以将其作为水文地质单元划分的依据。然而，依靠现有的科技水平，我们还无法将地下水和地表水的沟通情况准确确定，所以区域的地形地貌等只能作为进行水文地质单元的辅助划分因素。 　　2.2地层条件的影响 　　地层条件必须首先考虑地层岩性。由于不同岩性的裂隙发育程度会有明显差异，因此其富水性也会有所不同。所以在岩性变化较大的地区应进行分区处理。但如果是要研究某一含水层的地下水时情况时，由于其岩石的性质较为相似，故此时岩石性质应作为水文地质单元的辅助划分因素。其次，还要考虑不同岩性的连通及渗透情况。一般而言，地下水在其中的运行情况会大大影响地下水水害的程度。最后，我们还要考虑隔水层和含水层的分布、产状及其在地表的出露情况。如图1所示为潜水、承水以及滞水含水层的产状及分布情况，这些都会在一定程度上影响到水文地质单元的划分。但在实际考察中发现，由于隔水层和含水层的产状及分布特性可以很容易地在大范围内被掌握。而局部小范围的缺失或是厚度变薄，则无法精确判断，因此可能会造成含水层的沟通，进而影响到整个区域的局部径流、补水和排泄条件。 此时如果还是按照原来的隔水层及含水层的性质来划分水文地质单元，则会造成较大误差。 　　2.3地质构造的影响 　　地质构造是影响水文地质单元划分的主要因素。我们既要考虑研究区域的断裂带、区域性断层的分布情况，还要掌握该区域的断层破碎带以及节理密集带的显著特点。一般来说，断层分布区的规模越大，该区域表现为岩溶水强径流带的概率就越大。但是，如果该断层是不导水的，那么就可以划分为2个水文地质单元。另外，由于断裂或节理的存在，使得断裂或节理两端的地层径流通道堵塞或畅通，从而形成2个不同的径流、补给和排泄的区域，导致该区域可能成为2个水文地质单元。 　　2.4气象及水文因素的影响 　　若在较大区域范围内进行水文地质单元的划分，那么气象及水文因素的影响较大。此时我们必须将其列为主要依据对象。但是若在特定的规模不大的地区内划分水文地质单元，由于此时气象和水文因素在此范围内相似甚至相同，所以不能将其作为主要的参考因素。 　　3 K均值聚类 　　3.1基本概念 　　k均值聚类是目前最著名的区域划分聚类算法，而且由于该方法简洁、效率高，逐渐成为所有聚类算法中应用最广泛的一个。它的基本原理是首先通过给定数据点集合以及需要的聚类数目k（k值由用户指定），然后再根据某个特定的距离函数反复地把数据分入k个聚类中直至结束。 　　3.2计算方法及步骤 　　首先，随机选取K个对象作为开始的聚类中心。然后，计算每个所要划分的对象与子聚类中心之间的距离。并将其分配给距离自己最近的聚类中心。聚类中心以及分配到其中的对象就代表一个聚类。一旦所有对象都被分配完毕，每个聚类的聚类中心就会被重新计算。该过程将不断重复，直至满足相关终止条件为止。（如图1所示）一般而言，终止条件既可以设定为没有对象被重新分配给不同的聚类，也可以设定为没有聚类中心再发生变化，还可以设定为误差平方和局部最小。 　　4实际案例分析 　　4.1基本资料 　　以中国平顶山煤田寒武系灰岩的水文地质单元划分为例，该区域的自然地理条件、构造条件如下： 　　（1）自然地理条件：东起洛岗一号正断层，南起庚组煤（一煤）层露头，西北至郏县断层，东北至襄郏断层。含煤面