基于数值模拟探究平遥二亩沟煤矿开采对洪山泉域岩溶地下水资源的影响

基于数值模拟探究平遥二亩沟煤矿开采对洪山泉域岩溶地下水资源的影响一、引言1.1研究背景与意义1.1.1研究背景煤炭作为一种重要的化石能源，在全球能源结构中一直占据着举足轻重的地位。在我国，煤炭更是能源体系的重要支柱，长期以来为国家的经济发展和能源供应提供了坚实保障。从能源生产方面来看，煤炭是火力发电的主要燃料之一，通过燃烧煤炭产生的热能转化为电能，为工业生产和居民生活提供了稳定的电力支持。在工业领域，煤炭也是钢铁生产中不可或缺的原料，焦煤在高温下经过干馏等工艺制成焦炭，用于高炉炼铁，支撑着钢铁行业的发展。此外，煤炭还广泛应用于化工生产，例如制取煤气、煤焦油、化肥等化工产品。平遥二亩沟煤矿作为煤炭开采行业的一员，在当地的经济发展中发挥着重要作用。其煤炭产量不仅满足了本地部分能源需求，还为周边地区的工业生产提供了燃料和原料支持。然而，煤炭开采活动往往伴随着一系列的环境问题，其中对地下水资源的影响尤为显著。洪山泉域岩溶地下水是该地区重要的水资源，其水质优良，水量稳定，为当地居民的生活用水以及工农业生产提供了关键的水源保障。洪山泉位于山西省介休市东北部约10km的洪山镇狐歧山脚下，是山西省著名的岩溶大泉之一。该泉主要是由于山前大断裂阻水，使得岩溶水在断层带产生地下雍水并溢出地表成泉，属于全排型断层溢流泉。泉水出口高程916m，大部分分布在几百米长度内，由小泉池、七星泉、源神池泉、黑虎泉和槐柳泉组成，形成集中排泄的泉群。随着平遥二亩沟煤矿开采活动的持续进行，其对洪山泉域岩溶地下水资源的影响逐渐显现。煤矿开采过程中的矿坑排水、采煤活动导致的地层结构变化等，都可能改变岩溶地下水的补给、径流和排泄条件，进而影响地下水资源量和水质。因此，深入研究平遥二亩沟煤矿开采对洪山泉域岩溶地下水资源的影响具有重要的现实意义。1.1.2研究意义本研究在理论和实践方面都具有重要价值。从理论角度来看，有助于丰富煤矿开采与地下水资源相互作用的研究内容。目前，虽然已有不少关于煤矿开采对地下水资源影响的研究，但不同地区的地质条件和开采情况各异，平遥二亩沟煤矿与洪山泉域的独特地质背景和开采现状，为该领域的研究提供了新的案例和数据。通过深入分析两者之间的相互关系，可以进一步完善煤矿开采对地下水资源影响的理论体系，为后续相关研究提供参考和借鉴。在实践层面，研究成果对当地水资源保护和煤矿开采可持续发展具有重要指导意义。对于水资源保护而言，明确煤矿开采对洪山泉域岩溶地下水资源的影响，能够为制定科学合理的水资源保护措施提供依据。可以根据影响的程度和范围，划定重点保护区域，制定相应的水资源管理政策，加强对岩溶地下水的监测和保护，确保当地居民的生活用水和工农业生产用水的安全稳定供应。对于煤矿开采可持续发展来说，了解开采活动对地下水资源的影响，有助于煤矿企业优化开采方案。企业可以在开采过程中采取相应的技术措施，如改进采煤方法、加强矿坑水的处理和回用等，减少对地下水资源的破坏，实现煤炭资源开发与水资源保护的协调发展，促进当地经济的可持续发展。1.2国内外研究现状在煤矿开采对地下水资源影响方面，国内外学者进行了大量研究。国外研究起步较早，在一些煤炭资源丰富且开采历史较长的国家，如美国、澳大利亚等，学者们通过长期的监测和分析，对煤矿开采导致的地下水位下降、含水层结构破坏以及泉水流量减少等问题进行了深入探讨。例如，美国在阿巴拉契亚地区的煤炭开采研究中发现，大规模的采煤活动使得该地区部分区域的地下水位显著下降，一些依赖地下水的生态系统受到严重破坏。澳大利亚的研究则表明，煤矿开采引发的地层变形会改变地下水的径流路径，导致水资源重新分配。国内研究也取得了丰硕成果。众多学者针对我国不同煤炭产区的特点，开展了广泛的研究。在陕北和神东煤炭基地，由于地处西部干旱与半干旱地区，煤炭资源开发与地下水资源保护矛盾突出。相关研究通过系统分析该地区的地质与水文地质条件，厘清了不同开采区保水目标层、阻水关键层、采动水文地质结构以及地下水资源失水模式的差异，并根据煤层采动导水裂隙与阻水关键层、保水目标层的组合关系，总结得出保水目标层的多种失水模式。在山西省襄垣县，研究人员通过实地调研、文献调查和数值模拟等方法，探究了煤炭开采对区域地下水资源的影响，并提出了相应的保护措施。在数值模拟应用于煤矿开采对地下水资源影响的研究方面，国内外也有诸多成果。数值模拟方法能够定量描述煤矿开采过程中地下水流场的变化、预测矿坑涌水量以及评估对周边地下水环境的影响程度。常用的数值模拟软件如VisualModflow、地下水数值模拟系统（GMS）等，在国内外的相关研究中得到了广泛应用。在神府矿区大柳塔井田的研究中，通过建立平面二维地下水渗流模型，利用有限元法对数学模型进行数值求解，定量分析和预测了疏干降水后地下水流场的变化及其发展趋势，为水资源利用和生态环境保护提供了科学依据。然而，针对平遥二亩沟煤矿开采对洪山泉域岩溶地下水资源影响的研究相对较少。虽然已有一些关于洪山泉域的水文地质研究以及个别煤矿对地下水资源影响的案例分析，但专门针对平遥二亩沟煤矿与洪山泉域岩溶地下水相互作用的深入研究还存在不足。在数值模拟方面，缺乏对该地区复杂地质条件和开采活动进行精确模拟的研究，无法全面准确地评估煤矿开采对洪山泉域岩溶地下水资源的影响。因此，开展本研究具有一定的必要性和紧迫性。1.3研究内容与技术路线1.3.1研究内容本研究主要涵盖以下几个方面的内容：首先，对平遥二亩沟煤矿及洪山泉域的地质概况进行详细分析。包括研究区域的地层岩性、地质构造、水文地质条件等，全面了解该地区的地质背景，为后续研究奠定基础。具体而言，要明确洪山泉域岩溶地下水的补给、径流和排泄特征，以及二亩沟煤矿开采区域与岩溶地下水系统的空间关系。其次，建立数值模型。根据地质概况分析结果，选取合适的数值模拟软件，如VisualModflow，建立洪山泉域岩溶地下水的数值模型。确定模型的边界条件、初始条件以及相关参数，确保模型能够准确反映研究区域的水文地质特征。在模型建立过程中，要充分考虑煤矿开采活动对地下水流场的影响因素，如矿坑排水、采煤引起的地层变形等。再者，利用建立的数值模型进行模拟分析。模拟不同开采方案下煤矿开采对洪山泉域岩溶地下水资源的影响，包括地下水位变化、水流场改变以及泉流量的变化等。通过对比分析不同模拟结果，评估煤矿开采对岩溶地下水资源的影响程度和范围。最后，根据模拟结果提出相应的保护措施。针对煤矿开采对洪山泉域岩溶地下水资源造成的影响，从工程技术、管理政策等方面提出切实可行的保护建议，以实现煤炭资源开发与水资源保护的协调发展。1.3.2技术路线本研究的技术路线如图1-1所示。首先进行资料收集，广泛收集平遥二亩沟煤矿及洪山泉域的地质、水文地质、开采现状等相关资料，同时查阅国内外相关研究文献，了解研究现状和发展趋势。然后进行野外调查，对研究区域进行实地勘查，包括地质构造观测、水文地质参数测量、泉点流量监测等，获取第一手数据资料，补充和验证收集到的资料信息。在资料收集和野外调查的基础上，进行地质概况分析，明确研究区域的地层岩性、地质构造、水文地质条件等。接着建立数值模型，根据地质概况分析结果，确定模型的边界条件、初始条件和相关参数，利用VisualModflow等软件建立洪山泉域岩溶地下水的数值模型，并对模型进行识别和验证，确保模型的准确性和可靠性。利用建立好的数值模型进行模拟分析，设置不同的开采方案，模拟煤矿开采对洪山泉域岩溶地下水资源的影响，分析模拟结果，评估影响程度和范围。最后根据模拟分析结果，从工程技术、管理政策等方面提出保护洪山泉域岩溶地下水资源的措施和建议，形成研究结论。图1-1技术路线图二、研究区概况2.1平遥二亩沟煤矿概况平遥二亩沟煤矿位于山西省晋中市平遥县段村镇文祠神村，其地理坐标处于东经[X1]°至[X2]°，北纬[Y1]°至[Y2]°之间。该煤矿交通较为便利，周边有县级公路和乡道与外界相连，距离主要交通干线也较近，这为煤炭的运输提供了良好的条件，使得煤炭能够较为便捷地运往周边地区的工业企业和能源需求场所。在开采规模方面，平遥二亩沟煤矿为斜井开拓方式，井筒数量3个，核定生产能力为120万吨/年。从产量来看，近年来该煤矿的年产量保持在一定水平，为当地的经济发展提供了重要的能源支撑。其煤炭产品主要供应给周边的火电企业、钢铁厂以及化工企业等，满足这些企业的生产需求，促进了相关产业的发展。该煤矿成立于1989年7月3日，有着较为悠久的开采历史。在过去的几十年里，煤矿的开采技术和管理水平不断发展。早期，开采技术相对落后，生产效率较低，安全保障措施也不够完善。随着科技的进步和行业的发展，煤矿逐渐引入先进的开采设备和技术，如综合机械化采煤技术，大大提高了生产效率，同时加强了安全管理体系建设，提高了安全生产水平。在当地经济中，平遥二亩沟煤矿占据着重要地位。它不仅为当地提供了大量的就业岗位，直接吸纳了周边村庄的劳动力，解决了部分居民的就业问题，带动了当地居民收入的增长。还通过缴纳税收等方式，为平遥县的财政收入做出了重要贡献，有力地支持了当地的基础设施建设、教育、医疗等公共事业的发展，促进了当地经济的繁荣和社会的稳定。2.2洪山泉域概况洪山泉域位于山西省介休市东部山区及丘陵区，在行政区划上地跨晋中市和长治市的介休、平遥和沁源三个县，总面积632km²。泉域内碳酸盐岩裸露区面积为260km²，占泉域总面积的41%。其地理坐标范围为东经[具体东经范围]，北纬[具体北纬范围]。从地质构造特征来看，洪山泉域位于山西台隆、晋东南坳褶带、沁水坳陷和晋中新裂陷的接合部，为总体走向NNE、并向NE倾伏的背斜，西北部为新裂陷。区域内次级构造以断裂为主，褶曲次之。这种复杂的地质构造对岩溶地下水的形成和运移产生了重要影响。断裂构造为地下水的流动提供了通道，使得岩溶地下水能够在不同的含水层之间进行交换和运移；而褶曲构造则影响了地层的分布和含水层的形态，从而影响了地下水的储存和径流条件。洪山泉域属温带大陆性季风气候区，年降水量为400-600mm，降水主要集中在夏季，且多以暴雨形式出现。这种气象条件对岩溶地下水的补给有着显著影响。夏季的集中降水使得地表径流增加，部分降水通过岩溶裂隙、落水洞等通道迅速补给岩溶地下水，导致地下水位在短时间内迅速上升。而在其他季节，由于降水量较少，地下水的补给相对不足，地下水位会逐渐下降。泉域内的主要河流有龙凤河、樊王河、东涧河和西涧河，均为季节性河流。龙凤河是汾河的一级支流，也是泉域内最大的河流，发源于沁源县，由红崖沟、才子沟、后沟和铁水沟四条主沟洪水在古寨汇合而成，汇合后由大栅进入长达15km的绵山峡谷，经过介休市石口子水库后，到龙凤村南出口，在洪相村西北汇入汾河，全长52km。这些河流的水文条件对岩溶地下水的补给、径流和排泄有着重要作用。在洪水期，河流流量增大，河水通过河床渗漏等方式补给岩溶地下水；而在枯水期，岩溶地下水则会补给河流，维持河流的基本流量。2.3煤矿与泉域的关系从空间位置上看，平遥二亩沟煤矿位于洪山泉域的[具体方位]，煤矿开采区域与洪山泉域岩溶地下水系统存在着紧密的联系。煤矿开采区域部分位于洪山泉域的补给区范围内，这使得煤矿开采活动可能会对岩溶地下水的补给过程产生影响。在煤矿开采过程中，矿坑排水是一个重要的环节。大量的矿坑排水会导致开采区域地下水位下降，形成降落漏斗。由于煤矿开采区域与洪山泉域岩溶地下水系统存在水力联系，这种水位下降可能会影响到岩溶地下水的径流方向和流速。原本向洪山泉流动的岩溶地下水可能会因为煤矿开采形成的水位差而改变流向，流向煤矿开采区域，从而减少了洪山泉的补给量。煤矿开采活动还可能会破坏地层结构，导致岩溶裂隙和溶洞的连通性发生变化。这可能会改变岩溶地下水的径流路径，使得岩溶地下水在运移过程中遇到阻碍或发生分流，进而影响到洪山泉域岩溶地下水的正常循环和排泄，对洪山泉的流量和水质产生潜在影响。三、地质及水文地质条件分析3.1地层与构造研究区内出露的地层较为丰富，从老到新依次有太古界变质岩系、元古界长城系、蓟县系、寒武系、奥陶系、石炭系、二叠系、三叠系、侏罗系、白垩系以及新生界第四系。太古界变质岩系主要由片麻岩、混合岩等组成，岩石经历了强烈的变质作用，结晶程度高，质地坚硬，是区域内最古老的基底岩石，为后续地层的沉积奠定了基础。元古界长城系和蓟县系以碎屑岩和浅变质岩为主，岩性较为复杂，包含砂岩、页岩等，这些地层记录了元古代时期的沉积环境和地质演化过程。寒武系和奥陶系主要为海相沉积的石灰岩和白云岩，岩性致密，岩溶发育程度较高，是岩溶地下水的重要赋存层位。寒武系石灰岩中常含有丰富的三叶虫化石，反映了当时海洋生物的繁盛。奥陶系白云岩的形成与海水的化学沉积作用密切相关，其良好的岩溶发育条件为岩溶地下水的储存和运移提供了有利空间。石炭系和二叠系为海陆交互相沉积地层，主要由砂岩、泥岩、煤层等组成。其中煤层是煤矿开采的主要对象，石炭系和二叠系的煤层厚度和分布范围对平遥二亩沟煤矿的开采规模和布局有着重要影响。二叠系的砂岩常具有较好的透水性，在一定程度上影响着地下水的径流和储存条件。三叠系、侏罗系和白垩系以陆相碎屑岩沉积为主，岩性包括砂岩、页岩、砾岩等，这些地层的沉积特征反映了中生代时期的陆地环境和地质变迁。新生界第四系主要为松散堆积物，如砂、砾石、黏土等，分布于地表，厚度变化较大，其透水性和富水性因岩性和堆积条件而异，对地表水体与地下水之间的水力联系有着重要作用。研究区地质构造复杂，褶皱和断裂构造发育。褶皱构造主要表现为一系列轴向近南北向的背斜和向斜，背斜轴部地层受到拉伸作用，岩石破碎，裂隙发育，有利于岩溶地下水的赋存和运移；向斜轴部地层则相对紧闭，岩石完整性较好，地下水的储存条件相对较差，但在向斜的翼部，由于地层的倾斜和裂隙的发育，也可能存在一定规模的岩溶地下水。断裂构造对岩溶地下水的控制作用更为显著，不同方向的断裂相互交错，形成了复杂的地下水通道网络。一些大的断裂带往往是地下水的强径流带，岩溶地下水沿着断裂带快速流动，水力坡度较大；同时，断裂带还可能沟通不同含水层，导致含水层之间的水力联系增强，使得岩溶地下水的补给、径流和排泄条件更加复杂。此外，断裂构造还可能影响地层的隔水性能，破坏原有的水文地质结构，从而改变岩溶地下水的流动路径和分布规律。例如，当断裂切穿隔水层时，可能导致不同含水层之间的水力联系发生变化，原本相对独立的含水层之间可能出现水力连通，进而影响岩溶地下水的水位和水质分布。3.2水文地质条件3.2.1含水层与隔水层研究区内主要含水层包括第四系松散岩类孔隙含水层、寒武-奥陶系岩溶裂隙含水层和石炭-二叠系碎屑岩裂隙含水层。第四系松散岩类孔隙含水层主要分布于河谷、山前平原等地形相对低洼的区域，由砂、砾石等松散堆积物组成，厚度一般在[X1]-[X2]米之间。该含水层透水性较好，富水性较强，主要接受大气降水和地表水的入渗补给，水位受季节变化影响较大。在雨季，随着降水量的增加，含水层水位迅速上升；而在旱季，由于补给量减少和蒸发作用，水位会逐渐下降。其地下水动态变化较为明显，水质一般较好，矿化度较低，多为HCO₃-Ca型水，是当地居民生活用水和部分农业灌溉用水的重要水源之一。寒武-奥陶系岩溶裂隙含水层是研究区内最重要的含水层，也是洪山泉域岩溶地下水的主要赋存层位。该含水层主要由石灰岩、白云岩等碳酸盐岩组成，岩溶裂隙发育，富水性强。在区域上，该含水层的厚度较为稳定，一般在[X3]-[X4]米之间，但岩溶发育程度在不同地段存在差异。在构造破碎带、褶皱轴部等部位，岩溶裂隙更为发育，富水性更强；而在相对完整的岩层区域，岩溶发育程度相对较弱。该含水层主要接受大气降水在裸露岩溶区的入渗补给，以及地表水通过岩溶通道的渗漏补给。其地下水径流速度相对较慢，但由于含水层厚度大、富水性强，储存量巨大。该含水层的水化学类型主要为HCO₃-Ca・Mg型水，水质优良，是洪山泉的主要补给水源，对当地的工农业生产和生态环境起着至关重要的作用。石炭-二叠系碎屑岩裂隙含水层主要由砂岩、泥岩互层组成，砂岩为主要的含水介质，泥岩则起到相对隔水的作用。该含水层厚度变化较大，一般在[X5]-[X6]米之间，富水性中等。其透水性和富水性主要取决于砂岩的裂隙发育程度和连通性，在裂隙发育较好的地段，富水性相对较强；而在裂隙不发育或被泥质充填的地段，富水性较弱。该含水层主要接受上部含水层的越流补给和侧向径流补给，地下水径流方向受地形和构造控制。其水化学类型较为复杂，常见的有HCO₃-Ca・Mg型、SO₄-Ca・Mg型等，水质一般较差，矿化度相对较高。主要隔水层包括第四系黏性土隔水层、石炭-二叠系泥岩隔水层和奥陶系中统下马家沟组底部隔水层。第四系黏性土隔水层广泛分布于第四系松散堆积物的下部，厚度一般在[X7]-[X8]米之间。黏性土颗粒细小，孔隙度小，透水性差，能够有效阻隔上部孔隙含水层与下部含水层之间的水力联系，对保持上部孔隙含水层的水位稳定和水质起到重要作用。石炭-二叠系泥岩隔水层位于石炭-二叠系碎屑岩地层中，泥岩分布广泛，厚度较大，一般在[X9]-[X10]米之间。泥岩具有良好的隔水性能，能够阻止砂岩含水层之间的水力联系，使得石炭-二叠系碎屑岩裂隙含水层在一定程度上相对独立，控制着地下水的径流和分布范围。奥陶系中统下马家沟组底部隔水层主要由泥质灰岩、页岩等组成，厚度一般在[X11]-[X12]米之间。该隔水层位于寒武-奥陶系岩溶裂隙含水层的下部，对阻止下部深层地下水向上补给岩溶裂隙含水层起到关键作用，维持了岩溶裂隙含水层的水位和水量稳定，保证了洪山泉域岩溶地下水系统的正常运行。3.2.2地下水补径排条件洪山泉域岩溶地下水的补给来源主要包括大气降水入渗补给、地表水渗漏补给以及侧向径流补给。大气降水入渗补给是岩溶地下水的重要补给方式之一。泉域内年降水量为400-600mm，降水主要集中在夏季，且多以暴雨形式出现。在岩溶裸露区，降水能够迅速通过岩溶裂隙、落水洞等通道直接渗入地下，补给岩溶地下水。据研究，大气降水入渗补给量约占岩溶地下水总补给量的[X13]%。在一些岩溶发育强烈的区域，如介休市东部山区，降水入渗系数可达到[X14]以上，使得大量降水能够快速转化为岩溶地下水。地表水渗漏补给也是岩溶地下水的重要补给来源。泉域内的主要河流如龙凤河、樊王河、东涧河和西涧河等，在流经岩溶发育地段时，河水通过河床渗漏、岩溶管道等途径大量补给岩溶地下水。特别是在洪水期，河流流量增大，地表水渗漏补给量显著增加。例如，龙凤河在石口子水库以下河段，由于河床下伏岩溶裂隙发育，河水渗漏严重，每年对岩溶地下水的渗漏补给量可达[X15]立方米。侧向径流补给主要来自泉域周边地区的岩溶地下水。泉域边界存在一些断裂构造和透水岩层，使得周边地区的岩溶地下水能够通过这些通道向泉域内流动，形成侧向径流补给。根据水文地质调查和数值模拟分析，侧向径流补给量约占岩溶地下水总补给量的[X16]%。在泉域的东北部边界，由于存在一条区域性的断裂构造，该断裂带成为了岩溶地下水的侧向径流通道，周边地区的岩溶地下水通过该断裂带向泉域内补给，对维持泉域内岩溶地下水的水位和水量平衡起到了一定的作用。岩溶地下水的径流路径受地质构造和地形地貌控制。在岩溶水系统中，地下水主要沿着岩溶裂隙、溶洞等通道流动。在构造破碎带和褶皱轴部，岩溶发育强烈，地下水径流速度较快，形成强径流带；而在岩溶发育相对较弱的区域，地下水径流速度较慢。从地形地貌来看，岩溶地下水总体上由地势高处向地势低处径流，流向洪山泉排泄区。在山区，岩溶地下水主要沿着山体的坡降方向径流，逐渐汇聚到山谷中，然后通过岩溶管道和裂隙向平原地区流动。在平原地区，岩溶地下水则在含水层中缓慢流动，最终向洪山泉方向排泄。洪山泉是岩溶地下水的主要排泄方式，泉水流量稳定，多年平均流量为[X17]立方米/秒。此外，岩溶地下水还通过人工开采、向河流排泄以及越流排泄等方式进行排泄。人工开采主要用于当地的工农业生产和居民生活用水，随着经济的发展和人口的增加，人工开采量呈逐渐上升趋势。向河流排泄是指在枯水期，岩溶地下水水位高于河流水位时，地下水通过岩溶管道和裂隙向河流排泄，维持河流的基本流量。越流排泄主要发生在不同含水层之间存在水力联系且水头差较大的情况下，岩溶地下水通过弱透水层向相邻含水层排泄。补径排条件的变化对岩溶地下水水位和水质有着重要影响。当补给量减少时，如降水量减少、地表水渗漏补给受阻等，岩溶地下水水位会下降，可能导致泉水流量减少，影响当地的供水安全和生态环境。而当补给量增加时，水位会上升，但如果补给水源受到污染，如工业废水、生活污水等未经处理直接排入地表水体，通过地表水渗漏补给岩溶地下水，可能会导致岩溶地下水水质恶化。径流条件的改变，如煤矿开采导致地层结构破坏，改变了岩溶裂隙和溶洞的连通性，可能会使地下水径流路径发生变化，影响地下水的流速和流向，进而影响水位和水质分布。排泄条件的变化，如人工开采量过大，超过了岩溶地下水的补给能力，会导致水位持续下降，形成降落漏斗，还可能引发地面沉降等地质灾害；而排泄受阻，如泉水出口被堵塞，会使岩溶地下水水位升高，可能导致淹没等问题。3.2.3岩溶水系统特征洪山泉域岩溶水系统结构复杂，可分为补给区、径流区和排泄区。补给区主要位于泉域的东部和南部山区，碳酸盐岩裸露，岩溶发育，大气降水和地表水能够直接入渗补给岩溶地下水。在补给区，岩溶裂隙和落水洞等岩溶形态广泛分布，形成了良好的入渗通道，使得降水能够迅速转化为岩溶地下水。例如，在沁源县境内的补给区，岩溶裂隙密度可达每平方米[X18]条以上，落水洞数量众多，平均每平方公里可达[X19]个，这些岩溶形态大大提高了降水入渗效率。径流区位于补给区和排泄区之间，岩溶地下水在径流区内沿着岩溶通道和裂隙流动。径流区内地质构造复杂，褶皱和断裂发育，控制着岩溶地下水的径流方向和速度。在构造破碎带和褶皱轴部，岩溶发育强烈，地下水径流速度较快；而在相对完整的岩层区域，径流速度较慢。径流区内的岩溶管道和裂隙相互交织，形成了复杂的地下水流动网络，使得岩溶地下水能够在不同区域之间进行交换和运移。排泄区主要位于洪山泉出露地带，岩溶地下水以泉水的形式排泄出地表。洪山泉是全排型断层溢流泉，由于山前大断裂阻水，岩溶水在断层带产生地下雍水并溢出地表成泉。泉水出口高程916m，由小泉池、七星泉、源神池泉、黑虎泉和槐柳泉等组成泉群，集中排泄岩溶地下水。岩溶水系统的边界条件包括自然边界和人为边界。自然边界主要有地形边界、地质构造边界和含水层边界。地形边界表现为泉域周边的山脉和高地，这些地形屏障限制了岩溶地下水的侧向流动范围。地质构造边界如区域性断裂构造，控制着岩溶水系统与周边地区的水力联系。含水层边界则是指岩溶含水层与其他含水层之间的分界线，不同含水层之间的水力联系较弱，使得岩溶水系统在一定程度上相对独立。人为边界主要是由于人类工程活动形成的，如煤矿开采形成的采空区和矿坑排水，可能会改变岩溶水系统的边界条件，影响岩溶地下水的流动和分布。岩溶水系统的水动力特征表现为地下水的流动具有明显的方向性和层次性。在垂直方向上，岩溶地下水存在垂向分带现象，可分为包气带、季节变动带和饱水带。包气带位于地表以下，地下水水位以上，主要接受大气降水入渗，含水量随季节变化较大。季节变动带位于包气带和饱水带之间，地下水水位随季节变化而波动，该带内岩溶发育相对较弱，地下水流动速度较慢。饱水带位于季节变动带以下，充满地下水，是岩溶水系统的主要储水和径流区域，岩溶发育强烈，地下水流动速度相对较快。在水平方向上，岩溶地下水从补给区向排泄区流动，水力坡度在不同地段存在差异，在强径流带水力坡度较大，而在弱径流带水力坡度较小。洪山泉域岩溶水系统在区域水资源中占据重要地位。其岩溶地下水水质优良，水量稳定，是当地居民生活用水、工农业生产用水和生态用水的重要水源。洪山泉的泉水不仅为周边居民提供了优质的生活饮用水，还为当地的农业灌溉提供了保障，促进了农业的发展。在工业方面，岩溶地下水被广泛应用于煤炭开采、化工、电力等行业，支持了当地工业的发展。此外，洪山泉域的岩溶地下水还对维持区域生态平衡起着关键作用，为泉域内的河流、湖泊等水体提供了稳定的补给，保障了生态系统的健康稳定运行。四、数值模拟方法与模型建立4.1数值模拟原理与软件选择数值模拟在地下水研究中具有重要意义，其原理基于地下水动力学的基本理论，通过数学方法对地下水系统的复杂过程进行定量描述。地下水动力学研究地下水在岩石空隙中的运动规律，包括达西定律等基本定律，这些定律是数值模拟的理论基础。数值模拟的核心思想是将实际的地下水系统离散化为有限个单元或节点，通过建立相应的数学模型来近似求解地下水的运动方程，从而得到研究区域内地下水的水位、流量等要素在空间和时间上的分布情况。在众多数值模拟方法中，有限差分法、有限元法和边界元法是较为常用的方法。有限差分法是将求解区域划分为规则的网格，用差商来近似代替偏导数，将偏微分方程转化为代数方程组进行求解。例如，在一维稳定流问题中，对于地下水运动的基本方程\frac{\partial}{\partialx}(K\frac{\partialh}{\partialx})=0（其中K为渗透系数，h为水头），采用有限差分法进行离散时，将x方向划分为等间距的网格，网格间距为\Deltax，在节点i处，用差商近似偏导数，如\frac{\partialh}{\partialx}\approx\frac{h_{i+1}-h_{i}}{\Deltax}，代入原方程后可得到关于节点水头h_i的代数方程，通过求解这些代数方程即可得到各节点的水头值。有限元法是将求解区域划分为有限个单元，对每个单元建立插值函数，通过变分原理或加权余量法将偏微分方程转化为代数方程组。边界元法主要是将问题的边界离散化，通过求解边界积分方程来得到边界上的未知量，进而求解整个区域的问题。本研究选择VisualModflow软件作为数值模拟工具，该软件具有诸多功能优势。它是目前国际上最流行且被各国一致认可的三维地下水流和溶质运移模拟评价的标准可视化专业软件系统。该软件由加拿大Waterloo水文地质公司在原MODFLOW软件的基础上应用现代可视化技术开发研制，具有强大的图形可视界面功能。其设计新颖的菜单结构允许用户非常容易地在计算机上直接圈定模型区域和剖分计算单元，并可方便地为各剖分单元和边界条件直接在机上赋值，实现真正的人机对话。界面设计包括输入模块、运行模块和输出模块，各个模块相互联系又相对独立，方便用户操作。在环境保护、水资源利用与管理、采矿、建筑等许多行业和部门，VisualModflow都得到了广泛的应用，成为最为普及的地下水运移数值模拟的计算软件。在本研究中，选择该软件主要是因为其能够快速简便地导入各种常见文件类型和格式的数据，在导入时可自动进行坐标系和单位转换，并进行数据验证和智能错误检查，这对于处理本研究中涉及的大量地质、水文地质数据非常有利。它基于GIS的3D概念建模和数值建模都在一个集成的软件环境中进行，减少了昂贵的第三方预处理工具的使用，并消除了传输错误，能够提高建模的效率和准确性。此外，VisualModflow可以处理大量详细的高分辨率数据，能够模拟大型和复杂的区域规模的地下水系统，对于研究洪山泉域这种复杂的岩溶地下水系统具有很强的适用性。它还具有强大的数据可视化功能，能够以最新的2D、3D和多视图（FlexViewer）可视化所有数据，包括概念对象、数字网格、输入数据（属性区域、边界条件单元）和输出数据（计算的水头、路径、水位等），方便对模拟结果进行分析和展示。4.2水文地质概念模型建立建立水文地质概念模型是数值模拟的关键步骤，其目的是对研究区复杂的水文地质条件进行合理的概化，使其能够用数学模型进行描述和求解，同时又能准确反映实际的水文地质情况。在建立洪山泉域岩溶地下水的水文地质概念模型时，需要对含水层结构、水力特征和边界条件等进行详细的分析和概化。4.2.1含水层结构概化根据对研究区地层与构造以及水文地质条件的分析，洪山泉域内主要含水层为寒武-奥陶系岩溶裂隙含水层，该含水层是岩溶地下水的主要赋存层位。对其结构进行概化时，考虑到其在区域上的分布相对稳定，将其概化为一层具有一定厚度的连续介质。虽然该含水层在不同地段的岩溶发育程度存在差异，但为了简化模型，在概化过程中，将其视为在平面上具有一定均质性的含水层，忽略局部岩溶发育的微小差异。对于其厚度，根据地质勘探数据和钻孔资料，确定其平均厚度为[X]米，并将其作为模型中的含水层厚度参数。4.2.2水力特征概化洪山泉域岩溶地下水的水力特征较为复杂，在概化过程中，主要考虑其水流状态和水力坡度。根据水文地质调查和监测数据，该区域岩溶地下水的流动基本符合达西定律，因此将其水流状态概化为符合达西定律的稳定流或非稳定流。在天然状态下，岩溶地下水的水力坡度相对较小，水流较为缓慢；而在受到煤矿开采等人类活动影响时，水力坡度可能会发生变化。在概化过程中，根据不同的模拟情景，分别确定其水力坡度参数。对于地下水的流速，根据达西定律v=KI（其中v为流速，K为渗透系数，I为水力坡度），结合研究区的渗透系数和水力坡度数据进行计算和概化。4.2.3边界条件概化边界条件的概化对模型的准确性至关重要，洪山泉域岩溶地下水系统的边界条件包括自然边界和人为边界。自然边界主要有地形边界、地质构造边界和含水层边界。地形边界方面，泉域周边的山脉和高地形成了天然的地形屏障，限制了岩溶地下水的侧向流动范围，将这些地形边界概化为隔水边界，即法向方向水力梯度（或流量）等于零的边界。地质构造边界上，一些区域性断裂构造控制着岩溶水系统与周边地区的水力联系。对于导水的断裂构造，若计算区内为富含水层，区外为弱含水层时，则将其概化为流量边界；若断裂带本身是导水的，计算区内为导水性较弱的含水层，而区外为强导水的含水层时，则将其概化为定水头补给边界。含水层边界方面，岩溶含水层与其他含水层之间的分界线，由于水力联系较弱，将其概化为隔水边界。人为边界主要考虑煤矿开采活动的影响。平遥二亩沟煤矿的开采形成了采空区和矿坑排水，这些因素改变了岩溶水系统的边界条件。矿坑排水导致地下水位下降，形成降落漏斗，将矿坑排水区域概化为定流量排泄边界，根据煤矿的实际排水数据确定其排泄流量。采空区的存在可能会改变地层的渗透性和水力特征，在概化过程中，将采空区视为渗透性发生变化的区域，根据实际情况确定其渗透系数等参数。通过以上对含水层结构、水力特征和边界条件的概化，构建了能够反映洪山泉域岩溶地下水实际情况的水文地质概念模型，为后续的数学模型建立和数值模拟提供了基础。4.3数学模型建立与求解在水文地质概念模型的基础上，将其转化为数学模型是进行数值模拟的关键步骤。对于洪山泉域岩溶地下水系统，其数学模型主要基于地下水运动的基本方程和相应的定解条件。4.3.1基本方程地下水运动的基本方程通常采用三维非稳定流的偏微分方程来描述，其一般形式为：\frac{\partial}{\partialx}(K_x\frac{\partialh}{\partialx})+\frac{\partial}{\partialy}(K_y\frac{\partialh}{\partialy})+\frac{\partial}{\partialz}(K_z\frac{\partialh}{\partialz})+W=S_s\frac{\partialh}{\partialt}其中，h为水头（m）；x、y、z为空间坐标（m）；K_x、K_y、K_z分别为x、y、z方向的渗透系数（m/d）；W为源汇项（m^3/d），包括降水入渗补给、地表水渗漏补给、人工开采、矿坑排水等；S_s为贮水率（1/m）；t为时间（d）。在洪山泉域岩溶地下水系统中，考虑到岩溶含水层在平面上的各向异性相对较小，且垂向水流相对较弱，为简化计算，在模型中假设K_x=K_y=K（水平渗透系数），K_z相对较小可忽略不计，将三维问题简化为二维平面问题。同时，源汇项W根据研究区的实际情况进行具体分析和确定。降水入渗补给量根据多年降水数据和入渗系数进行计算；地表水渗漏补给量根据河流的流量、河床渗漏情况等因素确定；人工开采量根据当地的用水统计资料获取；矿坑排水则根据平遥二亩沟煤矿的实际排水数据确定。4.3.2定解条件定解条件包括初始条件和边界条件，它们是确定数学模型唯一解的必要条件。初始条件：初始条件是指在模拟开始时刻（t=0），研究区域内各点的水头分布情况。根据研究区的地下水水位监测数据，选取某一代表性时刻的水位数据作为初始条件，即h(x,y,0)=h_0(x,y)，其中h_0(x,y)为初始时刻（t=0）的水头分布。边界条件：边界条件与水文地质概念模型中的边界概化相对应，主要包括三类边界条件。一类边界条件（已知水头边界）：在洪山泉域中，一些与地表水有密切水力联系且水位相对稳定的区域，如与河流有统一水位的地段，可概化为一类边界条件。即给定边界上的水头值h(x,y,t)=h_1(x,y,t)，其中(x,y)为边界上的点，h_1(x,y,t)为已知的边界水头值，在模拟过程中保持不变。二类边界条件（已知流量边界）：对于一些侧向径流补给或排泄的边界，如泉域边界上的侧向径流流入或流出区域，以及矿坑排水区域等，可概化为二类边界条件。其数学表达式为K\frac{\partialh}{\partialn}=q(x,y,t)，其中\frac{\partialh}{\partialn}为水头在边界法向方向的导数，q(x,y,t)为已知的单位宽度流量，流入为正，流出为负。在矿坑排水区域，根据煤矿的排水数据确定q的值；在侧向径流边界，根据水文地质调查和相关研究确定q的大小和方向。三类边界条件（混合边界）：在实际情况中，可能存在一些边界条件既包含已知水头信息，又包含已知流量信息，这种边界可概化为三类边界条件。其数学表达式为K\frac{\partialh}{\partialn}+\alphah=\beta，其中\alpha和\beta为与边界条件相关的系数，可根据具体情况确定。4.3.3数值解法选择对于上述建立的数学模型，由于其复杂性，通常无法直接获得解析解，需要采用数值解法进行求解。本研究选择有限差分法对数学模型进行离散求解。有限差分法是一种将偏微分方程转化为代数方程组的数值方法，其基本思想是将求解区域划分为规则的网格，用差商来近似代替偏导数。在二维平面问题中，将研究区域在x和y方向上分别划分为等间距的网格，网格间距分别为\Deltax和\Deltay，时间步长为\Deltat。对于水头h在节点(i,j)处的偏导数，采用中心差分格式进行近似，如\frac{\partialh}{\partialx}\approx\frac{h_{i+1,j}-h_{i-1,j}}{2\Deltax}，\frac{\partialh}{\partialy}\approx\frac{h_{i,j+1}-h_{i,j-1}}{2\Deltay}，\frac{\partialh}{\partialt}\approx\frac{h_{i,j}^{n+1}-h_{i,j}^{n}}{\Deltat}（其中n表示时间步）。将这些差商近似代入基本方程中，得到关于节点水头h_{i,j}^{n}的代数方程组。通过对这些代数方程组的求解，即可得到不同时间步下各节点的水头值，从而模拟出地下水的运动过程。在求解过程中，利用VisualModflow软件内置的求解器进行计算。该软件提供了高效的计算算法，能够快速准确地求解大规模的代数方程组。在设置求解参数时，根据模型的特点和计算精度要求，合理确定迭代次数、收敛标准等参数，以确保计算结果的准确性和稳定性。例如，设置迭代次数为[X]次，收敛标准为水头变化小于10^{-5}m，当计算过程中满足收敛标准时，认为计算结果收敛，得到的水头值即为模拟结果。4.4模型参数确定与率定模型参数的准确确定是保证数值模拟结果可靠性的关键环节。在洪山泉域岩溶地下水数值模型中，需要确定的参数主要包括渗透系数、贮水率、降水入渗系数、地表水渗漏系数等。这些参数的取值直接影响到模型对实际水文地质过程的模拟精度。4.4.1参数初值确定通过收集整理研究区的地质、水文地质数据来确定模型参数的初值。地质勘探数据提供了含水层的岩性、厚度等信息，这些信息对于确定渗透系数等参数具有重要参考价值。例如，对于寒武-奥陶系岩溶裂隙含水层，根据岩石的类型和岩溶发育程度，参考相关的水文地质手册和研究成果，初步确定其渗透系数的取值范围为[X1]-[X2]m/d。水文地质试验数据也是确定参数初值的重要依据。通过抽水试验、渗水试验等，可以直接测定含水层的渗透系数、贮水率等参数。在研究区内，收集了以往进行的抽水试验资料，根据这些试验数据，采用裘布依公式等方法计算得到渗透系数和贮水率的初步值。降水入渗系数则根据研究区的气象数据、地形地貌和土壤类型等因素，参考相关地区的研究成果，初步确定为[X3]。地表水渗漏系数根据河流的水文特征、河床岩性和渗漏情况等，通过现场观测和分析，初步确定为[X4]。4.4.2模型率定与验证在确定参数初值后，需要对模型进行率定和验证，通过不断调整参数值，使模型的模拟结果与实际观测数据尽可能吻合。模型率定过程是一个反复试错和优化的过程，利用研究区的地下水水位观测数据、泉流量观测数据等进行对比分析。在VisualModflow软件中，通过调整模型参数，如渗透系数、贮水率等，观察模拟结果的变化，并与实际观测数据进行比较。以地下水水位模拟结果为例，将模拟得到的各观测点水位与实际观测水位进行对比，计算两者之间的误差，如均方根误差（RMSE）、平均绝对误差（MAE）等。根据误差大小，调整模型参数，若模拟水位普遍高于实际水位，则适当减小渗透系数等参数；反之，则增大相应参数。通过多次调整和试算，使误差逐渐减小，直到模拟结果与实际观测数据达到较好的拟合程度。模型验证是在模型率定的基础上，利用另一组独立的观测数据对模型的可靠性进行检验。将验证数据输入已率定好的模型中进行模拟，再次计算模拟结果与验证数据之间的误差。若误差在可接受范围内，说明模型具有较好的可靠性和预测能力；否则，需要进一步分析原因，对模型进行改进和优化，如检查数据的准确性、调整模型结构或参数等。通过模型率定和验证，最终确定了适合洪山泉域岩溶地下水系统的模型参数，使模型能够准确地模拟地下水的运动过程和变化规律，为后续的模拟分析提供了可靠的基础。五、模拟结果与分析5.1现状模拟结果分析通过将数值模拟得到的水位和水质数据与实测数据进行对比，能够直观地了解当前岩溶地下水的实际状况。在水位方面，从模拟结果与实测数据的对比图（图5-1）可以看出，在大部分观测点处，模拟水位与实测水位具有较好的一致性。例如，在泉域东部的[具体观测点1]，模拟水位与实测水位的差值在[X1]m以内，相对误差小于[X2]%，这表明模型能够较为准确地模拟该区域的水位分布情况。然而，在部分观测点也存在一定的差异。在煤矿开采区域附近的[具体观测点2]，模拟水位低于实测水位，差值达到[X3]m，相对误差为[X4]%。经分析，这可能是由于煤矿开采过程中，实际的矿坑排水情况较为复杂，存在一些难以准确量化的因素，如排水时间的不均匀性、排水点的分布等，导致模型在该区域的模拟精度受到一定影响。在水质方面，对岩溶地下水的主要离子浓度进行了模拟和实测对比。以钙离子（Ca²⁺）浓度为例，模拟结果显示，在泉域的大部分地区，钙离子浓度的模拟值与实测值较为接近，相对误差在[X5]%以内。但在靠近煤矿开采区的下游区域，模拟值与实测值出现了一定偏差，实测钙离子浓度高于模拟值。这可能是因为煤矿开采过程中，煤矸石的堆积和淋溶作用，使得更多的钙离子等矿物质溶解进入地下水中，而模型在考虑这些因素时存在一定的局限性，未能完全准确地反映实际的水质变化情况。从流场特征来看，现状模拟结果表明，岩溶地下水的主流方向总体上是从补给区向排泄区流动，与理论分析和实际地质条件相符。在补给区，由于降水入渗和地表水渗漏补给，地下水位较高，水力坡度较大，地下水流动速度相对较快。随着水流向径流区，水力坡度逐渐减小，地下水流动速度也逐渐变慢。在排泄区，如洪山泉附近，地下水以泉水的形式排泄出地表，水位相对稳定，流场较为集中。水质特征方面，现状下岩溶地下水的水质总体较好，主要为HCO₃-Ca・Mg型水。但在煤矿开采区及其周边，由于受到矿坑排水、煤矸石淋溶等因素的影响，水质出现了一定程度的恶化。在煤矿开采区附近，硫酸根离子（SO₄²⁻）浓度明显升高，部分区域的SO₄²⁻浓度超过了国家地下水质量标准中的Ⅲ类水标准，达到了[具体浓度值]mg/L。同时，铁（Fe）、锰（Mn）等重金属离子浓度也有所增加，这对当地的生态环境和居民用水安全构成了潜在威胁。图5-1模拟水位与实测水位对比图5.2不同开采方案模拟结果为了全面评估煤矿开采对岩溶地下水资源的影响，设定了多种开采方案进行模拟。方案一为当前开采强度下的持续开采，方案二为适度增加开采强度，方案三为扩大开采范围，方案四为改变开采方式（采用保水开采技术）。在方案一下，模拟结果显示，随着开采时间的延长，煤矿开采区域周边的地下水位持续下降。在开采1年后，以煤矿开采区为中心，形成了一个半径约为[X6]m的降落漏斗，漏斗中心的水位降深达到[X7]m。随着开采时间进一步延长至5年，降落漏斗半径扩大到[X8]m，水位降深达到[X9]m。在水流场方面，由于水位下降，岩溶地下水的径流方向发生了改变，原本流向洪山泉的部分水流转向了煤矿开采区域，导致洪山泉的补给量减少，泉水流量有所下降。方案二适度增加开采强度后，地下水位下降速度明显加快。开采1年后，降落漏斗半径就达到了[X10]m，水位降深为[X11]m，分别比方案一在相同时间内增加了[X12]m和[X13]m。5年后，降落漏斗半径扩大到[X14]m，水位降深达到[X15]m。水流场的改变更加显著，更多的岩溶地下水流向煤矿开采区域，洪山泉的补给量大幅减少，泉水流量下降幅度超过[X16]%。方案三扩大开采范围后，模拟结果表明，地下水位下降区域范围明显扩大。在开采1年后，水位下降区域几乎覆盖了整个扩大后的开采范围，降落漏斗形态变得更为复杂。5年后，水位下降区域继续向外扩展，与周边其他区域的水位差进一步增大。水流场也发生了较大变化，由于开采范围的扩大，更多的岩溶地下水径流路径被改变，导致整个泉域的水流场分布更加紊乱，洪山泉的补给量受到严重影响，泉水流量急剧下降，部分时段甚至出现了断流现象。方案四采用保水开采技术后，地下水位下降幅度明显减小。开采1年后，降落漏斗半径仅为[X17]m，水位降深为[X18]m，相比方案一有显著降低。5年后，降落漏斗半径扩大到[X19]m，水位降深达到[X20]m。水流场的改变相对较小，岩溶地下水的主流方向基本保持不变，洪山泉的补给量虽然有所减少，但减少幅度相对较小，泉水流量能够维持在相对稳定的水平，仅下降了[X21]%左右。不同开采方案下岩溶地下水的水质也发生了不同程度的变化。在方案一和方案二下，随着开采强度的增加，煤矿开采区域及其周边的水质恶化程度加剧。硫酸根离子、重金属离子等污染物浓度持续升高，污染范围也不断扩大。在方案三下，由于开采范围的扩大，更多的区域受到煤矿开采的影响，水质恶化区域相应扩大，且污染物浓度在部分区域达到较高水平。而在方案四下，由于采用了保水开采技术，有效减少了矿坑排水和煤矸石淋溶等对地下水的污染，水质恶化程度得到了明显控制，污染物浓度增长缓慢，污染范围也相对较小。5.3开采对岩溶地下水资源的影响从水位下降方面来看，煤矿开采对岩溶地下水水位产生了显著影响。随着开采活动的进行，大量的地下水被排出，导致开采区域及其周边地下水位持续下降。如前文所述，在不同开采方案下，均出现了明显的降落漏斗，且降落漏斗的范围和深度随着开采时间的延长和开采强度的增加而不断扩大和加深。这种水位下降会引发一系列问题，首先，会导致岩溶泉的流量减少甚至干涸，洪山泉作为当地重要的水源地，其流量的减少将直接影响到周边居民的生活用水和工农业生产用水。其次，地下水位下降可能会引起地面沉降等地质灾害，破坏地表建筑物和基础设施，对当地的生态环境和经济发展造成严重影响。水量减少也是煤矿开采对岩溶地下水资源的重要影响之一。煤矿开采过程中的矿坑排水直接导致了岩溶地下水的水量损失。随着开采规模的扩大和开采时间的延长，矿坑排水量不断增加，岩溶地下水的储存量逐渐减少。同时，由于水位下降，岩溶地下水的补给量也会相应减少，进一步加剧了水量减少的趋势。这将使得当地的水资源供需矛盾更加突出，制约工农业的可持续发展，对当地的生态系统平衡也会产生负面影响，如导致湿地萎缩、植被退化等。在水质变化方面，煤矿开采对岩溶地下水水质的影响不容忽视。矿坑排水中含有大量的有害物质，如硫酸根离子、重金属离子、悬浮物等，这些物质在排放过程中会渗入地下，污染岩溶地下水。煤矸石的堆积和淋溶也会向地下水中释放有害物质，导致水质恶化。水质恶化不仅会影响地下水的饮用安全，还会对农业灌溉和工业生产造成危害。在农业灌溉方面，污染的地下水可能会导致土壤污染，影响农作物的生长和品质；在工业生产中，不符合水质要求的地下水可能会对生产设备造成腐蚀，影响产品质量和生产效率。5.4影响因素分析开采强度是影响岩溶地下水资源的重要因素之一。随着开采强度的增加，矿坑排水量增大，对岩溶地下水的疏干作用增强，导致地下水位下降速度加快，降落漏斗范围扩大。在方案二中，适度增加开采强度后，地下水位下降速度和降落漏斗范围明显大于方案一，这充分说明了开采强度与地下水位下降之间的正相关关系。开采强度的增加还会导致水质恶化加剧，更多的有害物质被带入地下水中，污染范围和程度进一步扩大。地质构造对岩溶地下水资源的影响也较为显著。在断裂构造发育的区域，岩溶地下水的径流速度较快，水力联系密切。煤矿开采活动在这些区域更容易引发地下水的流动状态改变，导致水位下降和水质变化。例如，在泉域内的某断裂带附近进行煤矿开采时，由于断裂带的导水作用，矿坑排水能够迅速扩散，使得周边区域的地下水位快速下降，且水质污染范围也沿着断裂带迅速扩展。褶皱构造则影响着岩溶含水层的分布和形态，进而影响地下水的储存和径流条件。在褶皱轴部，岩溶发育程度较高，地下水储存量相对较大，但开采活动也更容易破坏含水层结构，导致地下水泄漏和水位下降。含水层特性也是影响岩溶地下水资源的关键因素。寒武-奥陶系岩溶裂隙含水层的岩溶发育程度、渗透系数等特性对地下水的流动和储存有着重要影响。岩溶发育程度高的区域，地下水的储存空间大，径流速度快，但同时也更容易受到煤矿开采的影响。当开采活动破坏了岩溶裂隙和溶洞的连通性时，会改变地下水的径流路径，导致水位下降和水量减少。渗透系数较大的含水层，地下水的流动速度快，矿坑排水等对地下水的影响能够迅速传播，使得影响范围扩大。而含水层的厚度和富水性则决定了岩溶地下水的储存量，厚度大、富水性强的含水层在受到煤矿开采影响时，能够在一定程度上缓冲水位下降和水量减少的速度，但长期的开采活动仍会对其造成严重破坏。六、结论与建议6.1研究结论本研究通过对平遥二亩沟煤矿开采对洪山泉域岩溶地下水资源影响的数值模拟分析，取得了以下主要结论：在地质及水文地质条件方面，洪山泉域地层岩性复杂，构造活动频繁，岩溶发育，为岩溶地下水的赋存和运移提供