蒙古Unkheltseg盆地地质间断的三维地下水流数值模拟_中文版

摘要

地下水模型真实地代表了像Unkheltseg盆地这样复杂系统的水文地质，对蒙古至关重要。

它们为研究、政府、管理战略和国内商业发展带来了源源不断的好处。在与Unkheltseg盆地类似的盆地间断控连接中，已有有限的校准三维数值瞬态模拟案例研究，本文提出的模型旨在解决这一问题。该盆地具有独特的地质控制，是Taikh山谷未来经济发展的关键供水资源。商业勘探项目产生了高质量的地质和水文地质数据，这是在整个盆地范围内成功进行模型模拟所必需的。使用MODFLOW-SURFACT中的“DRAIN Package”和“Fracture-Well Package FW4”，实现了充分表示水平和垂直流动方向所需的空间离散化，从而有效地约束了断层屏障上的补给和排放。该模型是在断控盆地建立长期监测方案的重要工具，可以预测短期和长期的区域影响。

介绍

地下水资源在蒙古经济中占有重要地位，工业用水和生活用水主要依靠地下水资源;大约80%的总消耗量来自地下水(Hasiniaina et al . 2010)。

蒙古的地下水资源建模在其他地方也很常见，但由于信息有限，包括出版的和政府的数据，历史上对简单分析评估的依赖，以及复杂的地质系统，对各种建模技术提出了限制，因此具有挑战性。为了有效地管理地下水资源，建立全面、逼真的模型至关重要。

在过去的几年里，许多从事需要地下水供应的项目(Oyu Tolgi, Energy Resources UHG煤矿，Burun Narran煤炭项目)的国际公司已经开始提供用于资源估算的三维数值模拟，并使用这些模拟计算资源寿命和开采极限，这是蒙古目前使用的分析方法无法实现的。

在过去的几年中，许多讨论文件讨论了模型如何最好地为决策支持过程服务的问题。先前的研究(如Gupta等人2012;Kresic and Mikszweski 2014;Nordstrom 2012;Simmons and Hunt 2012;Voss 2011)已经解决了模型参数应该调整的程度，以允许它复制过去的系统行为，作为其用于管理未来系统行为的前兆，以及在这种能力下使用它需要有多复杂(或其他)。在这种情况下，模型的目的是预测管理制度下系统的行为。在模型所需的预测仅部分受历史数据约束的情况下，选择适当的模型复杂性水平是最困难的(Doherty和Simmons 2013)。

断层含水层是分析和解释水力测试数据最复杂的地质环境之一，因为断层不仅具有高传输带的固有能力，而且还具有水力屏障的作用(Bense等人2013;Bredehoeft 1997;Cello等人2000;Evans等人1997;Folch and mass - pla 2008;Shan et al . 1995)。虽然我们对断层含水层的解释本质上仍然是线性的，但通过数值模拟的参数估计突出了与断层相关的水力障碍的存在。这些障碍在恒放电试验数据中很明显。

断裂带和断层长期以来一直被认为对地下水的流动和输送有重大影响，因为它们有助于改变含水层的有效渗透率。由于这个原因，已经有大量的研究旨在描述与断裂和断层有关的各种各样的现象，在各种尺度上(Allen和Michel 1999;Nordqvist et al . 1992;Yihdego and Becht 2013)。在地质复杂地形中，地下水转移的作用是具有挑战性的(例如Nelson和Mayo 2014;Yihdego and Webb 2014)。在断层含水层中，可能存在额外的复杂性，因为断层经常被观察到作为管道(hunton和Lundy 1979;Pimentel and Hamza 2014)和作为障碍(Ran et al . 2014;Rojstaczer, 1987)。因此，很难对断层的影响作出广泛的概括。

盆地地区的地下水资源估算最初是基于抽水试验的产量(Battumur 2009)，并通过对该方法的Cooper - Jacob简化进行分析解释(例如Freeze和Cherry 1979;羁绊1994)。这一估计表明，可能有足够的地下水供应附近的采矿活动，但需要数值模拟来确认和评估影响。雪山工程公司(SMEC)于2011年受委托进行了预可行性水平研究，并于2012年进行了地下水资源可行性水平评估，重点是将Unkheltseg盆地作为巴彦阿拉格勘探公司的矿山供水水源。

Unkheltseg盆地位于蒙古首都乌兰巴托以西约930公里处，与俄罗斯边境以南250公里处，位于Zavkhan省，Erdenekhairkhan Soum，该地区以陡峭的山地地形和宽阔的山谷为特征。该地区以逆断层为主，与山脉平行，垂直位移明显，部分地区可达100米。这种断层的几何形状经常使冲积系统中的地下水含水层断开，并且断层之间的水力连通性通常有限或根本没有。盆地体系似乎以断块为主，形成了一个侧面陡峭、底部相对平坦的盆地，从概念上看，它很像一个浴缸(水箱水库)。

盆地系统向Tost盆地南部的排水受屏障的渗透性和水力条件的控制，因为地下水的排放发生在一个地质不连续面上，这里假定这是一个断层隆起的基底。因此，该盆地与托斯特盆地南部之间不存在直接的水力联系，这种地下水储放关系无法通过简单的解析方程或一维或二维数值模拟程序轻松处理。

三维数值模拟的挑战是产生一个真实地模拟自然条件的模型。利用MODFLOW-SURFACT中的“DRAIN Package”和“Fracture-Well Package FW4”，成功模拟了观测到的水文地质条件。

在流域模型中，地下饱和域通常由线性或非线性的储流量函数表示(例如Fiorillo 2011;Kampf and Burges 2007;Rupp et al . 2009;Singh and Woolhiser 2002;文平等

2011)。在某些情况下，该函数的推导始于对饱和流的基于物理的描述，但此后做出的假设(例如恒定水头梯度或连续稳态)导致单值储流量函数(例如TOPMODEL;Beven and Kirkby, 1979)。这些储存-排放函数的优势在于它们在计算上非常简单，但这种表示意味着在模型元素(例如子集水区)的尺度上，在任何空间维度(x, y或z)中都没有明确的分布信息，并且模型被认为是零或一维的(例如Kampf和Burges 2007)。

将地下水域离散为二维或三维并对饱和流的控制偏方程进行数值求解的模型(例如Palma and Bentley 2007;Singh and Woolhiser 2002;Vandenbohede等人(2011)确实具有空间显式的优势。利用瞬态饱和流模型给出了受时变应力(如补给或抽水)作用的线性偏微分方程的解析解(例如Bidwell et al . 2008;Pulido-Velazquez等人2005;文平等，2011)。然而，这些解决方案可能仅限于含水层的均匀一维表示，这可能需要对系统进行严重简化(Rupp et al . 2009)。

模型是自然系统按比例缩小的简化表示。因此，模型的开发以系统的知识为前提。如何建立Unkheltseg盆地“储罐-水库”模型的解决方案是利用MODFLOW-SURFACT中的三维瞬态和稳态模型，该模型克服了一维或二维分析方程的缺点，因为空间信息(即水头)保留在一次和二次流动方向上，并且没有做出流量与含水层储水量一对一关系的限制性假设(Sloan 2000)。

图1研究区位置图。

这项工作的目的是为特殊/独特的Unkheltseg盆地选择合适的模型结构和参数化水平，并为蒙古学术界、相关行业和政府机构提供一个盆地尺度的三维地下水资源建模的成功例子。详细和准确的建模对蒙古至关重要，因为该盆地系统的不同寻常之处在于补给有限，该系统受地质控制，并且该系统可能具有最大地下水水位，当达到该水位时，可以将其排放到向下梯度的Tost盆地。Tost盆地的补给主要来自Unkheltseg盆地的排放，这一有趣的联系给该地区的资源管理带来了挑战。位于蒙古北部Taikh山谷的Unkheltseg盆地含水层系统(图1)，利用MODFLOWSURFACT，通过计算机三维数值稳态和瞬态流动模拟，评估了地下水资源为附近矿山供水的潜力以及开采的潜在影响。

本文介绍了概念模型和数值方法，并论证了利用瞬态仿真进行模型标定的重要性。目的是证明瞬态数值模拟作为一种重要工具的有效性，并解释来自复杂含水层的数据，特别是来自盆地尺度上的断层含水层的数据。

背景和数据来源本工作的目的不是复制SMEC报告，而是总结相关信息并介绍地下水模拟和结果。本文中提供的大部分信息以及相关的假设和解释都来自以下顾问的报告:EcoTrade(2007年)、Battumur(2009年)、SMEC(2011年、2012年)，这些报告目前是保密的。巴彦阿拉格勘探公司已同意SMEC(2012)开展的工作，并将在论文中发表相关解释。

这项工作提供了蒙古为数不多的已发表的水文地质数值评估之一，使蒙古科学界受益，并为从事类似盆地系统工作的顾问和学者提供了参考工具。

SMEC(2012)的研究重点是利用稳定和瞬态数值地下水流动模型，定量评估含水层作为满足采矿需求的供水水源的潜力。预计矿井平均需水量为9升/秒，最高需水量为17升/秒，最低需水量为4.2升/秒，拟建矿井寿命为7年(不含建设)。因此，该工程需要供水高达(2.82 9105立方米/年)9l /s，为期7年。此外，本研究还通过模拟长期开采来评估地下水开发的影响。

气候和拓扑

该盆地地区的气候特点是夏季温暖，冬季寒冷，温度范围在-38至30 摄氏度之间。在盆地以南约30公里处的Dorvuljin Soum有降水记录。过去10年的平均降雨量为90毫米，其中约90%的降雨量在夏季的6 - 8月。

轶事证据表明，由于地形的影响，盆地以西的山脉(由石灰岩组成)的降雨量可能更高。该地区为山地，海拔高度在1670 - 2160米之间，盆地底部海拔高度在1720 - 1800米之间，最低点在南部。它西边是石灰岩山，东边是丘陵。该盆地向南流入约30公里外的Zavkhan河。

该盆地没有永久的地表水水体(即地下水)。

泉水、湖泊、恒河)。地表流动只有在强降雨事件超过入渗速率后才会发生。

地质与水文地质

该盆地布满了第四纪松散沉积物，覆盖在由石灰岩(西部和西北部)、tonalite、石英、闪长岩和花岗岩(东部和南部)组成的Riphean古老基岩基底上(图2)。确定了两种主要的含水层类型:无约束至半约束洪积/崩积松散沉积物和约束基岩Riphean Shuvuun组(Ecotrade 2007)。松散含水层分布在山坡附近，穿过沟壑和干涸的河床，由砾石和卵石、小粘土透镜体和小巨石组成。预可行性报告和水文地质评估表明，上部和下部含水层由一个漏水、低渗透率的含水层分隔(SMEC 2011)。Shuvuun组是一个低产量的基岩含水层，相互连接和范围有限(Battumur 2009)。

表1总结了Unkheltseg盆地(Unkheltseg钻孔)和盆地南部(Tost和Oortsog钻孔)的钻孔测试结果。资料显示，Unkheltseg井的抽水时间为72 h，降深2 m, 2 h后达到完全采收率。地下水质量主要由钠、镁阳离子和氯化物、亚硫酸盐和碳酸盐阴离子主导。

盐度在800至1800毫克/升之间变化。环境地球科学(2015)73:419 - 4133碳酸盐含水层中没有一个样品具有“典型”的特征。该地区有很多碳质岩石，但水并没有显示出与这种岩石类型有很多重要的相互作用。水质总体较好，矿化度低，补水合理，停留时间短。

Battumur(2009)也描述了包含巴彦Airag矿床的Riphean Shuvuun组含水层。共钻了5个垂直井眼，其中2个为干井，3个与水相交。表2总结了气升试验的结果。

在4口井中进行了井下地球物理测井，以确定淤泥/粘土和砾石互层。采用的方法有:自然伽马(gamma) (GP)、密度伽马(gamma)、卡尺、单点电阻(PR)、自发电位(SP)、电导率。GeoMaster工程有限责任公司进行了诱导极化(IP)测量，使用100米的偶极子间距测量了5条总长24000米的线。

该数据由QGX蒙古有限责任公司根据现有钻孔进行校准，以提供沿线路到基底的深度。然后QGX使用它来绘制基底轮廓。

Unkheltseg盆地被解释为一个由输运沉积物充填的断块盆地。沉积物是松散的到弱硬化的，可能代表了几个阶段和沉积类型，包括:•洪积冰川冲刷沉积;•风成沙丘;•冲积/崩积沉积物从周围的山丘流入盆地;•洪水期间冲积物冲刷;•低能浅湖沉积在湿润期或在盆地最深处发育。

整个盆地的岩性层不是连续的。特别是粘土或砂质粘土层是不连续的，可能代表在盆地最深处形成的短暂湖泊，并可能根据气候和持续的沉降而移动。断层位置和推断断层分别如图3a和图2所示。截面和解释单元如图3b所示。从诱导极化测量、钻井结果和横截面来看，Unkheltseg盆地很明显是封闭的，可以被认为是一个“充满水和沉积物的浴缸”。Unkheltseg盆地被解释为一个充填断块盆地，因此整个盆地的岩性层不是连续的。地下水梯度为0.0015，非常平坦，仅在北部略有增加，达到0.008左右，这可能是对盆地结构和补给和排泄带的响应。盆地的构造与浴盆并没有太大的不同，盆地的边缘陡峭，底部平缓倾斜，这是通过钻井勘探计划和地球物理测量(SMEC 2012)确定的。非常平坦的地下水梯度要么表明高渗透性，至少在上层含水层内，要么表明低补给。

在监测钻孔中观察到水位下降的响应，并表明:•抽水钻孔的下降远远大于附近的监测钻孔，表明钻孔发育不良;•抽水的影响在上层含水层中传播得很快，在下层含水层中传播得较慢;•上层含水层的水位迅速下降，然后以比抽水孔和下层含水层慢得多的速度继续下降;•下部含水层的水位对抽水有延迟响应，并且水位模拟抽水孔中的响应。这表明下部含水层筛管的合理开发。

水位监测表明，地下水位上、下两层之间不存在压力差。这表明两者是相互关联的。

上部含水层被认为是半承压的，下部含水层被认为是漏承压的。采收率和抽水试验数据表明，上部含水层的渗透率(平均)是下部含水层的几倍。

自2010年以来进行的水位监测表明，随着时间的推移，水位变化很小，这表明在此期间，要么补给有限，要么盆地有一个溢出点，控制着地下水的最高海拔。北部地下水高程的增加反映了渗透率的降低和补给的可能性。

流域内的地下水补给主要来自直接入渗降水和周围丘陵和山脉的坡面径流。

这主要是通过地表径流发生的，地表径流从较高的高架山脉流向盆地边缘，并在那里渗透。地表暴露的基岩也可以通过横向和垂直运动进行补给，然后通过沉积物表面以下基岩的裂缝/结构向盆地沉积物进行补给。此外，补给可能来自落在崩积盆地沉积物上的降水的直接入渗或来自坡面流;这被认为只发生在非常强的降雨事件中。由于地形的影响，降水被认为在盆地周围的山脉较高。这表明大部分补给是通过山前渗透和基岩转移进行的。

地下水盐度低，停留时间短，周转温和。因此，补给被认为直接通过山前或直接入渗进行，蒸发损失小，地下水运移距离短。

盆地南端的地下水位为36米，北端的地下水位超过100米，对于通过毛细力或植物蒸腾造成的蒸散损失来说，地下水位太深了。由于地下水位的显著变化，假定盆地沉积物中的地下水通过基底流通过盆地与南托斯特盆地之间的断层屏障排出。盆地内没有永久或短暂的地表水体(如泉水、湖泊和河流)。在地下水高程和梯度稳定的情况下，根据18个月的监测，假设南侧流域的排水发生在流域沉积物的地下水位处于监测指示的水平时。从托斯特盆地流出的水最终到达位于南面30公里处的扎夫汗河，这条河被认为是地下水的主要接收地。少量排放通过抽取地下水供当地居民储存和家庭使用而发生。

水位监测(SMEC 2012)表明，上下含水层的地下水位没有差异，表明两个含水层是相互联系的。采收率和泵送试验数据(SMEC表2)矿区钻孔气举试验结果来源:Battumur(2009)。井位位置如图1所示钻孔数屈服(L/s)下沉(m)含水带厚度(m)导水率(m/d)透光率(m2 /d) CVZDH 122 1.1 51 31.4 0.38 11.93 CVZDH 124 0.027 42.6 28.5 0.0092 0.26 CVZDH 125 1.2 37.7 57.5 0.0857 4.93算术平均值36.33 0.16 5.71 4124 Environ Earth Sci(2015) 73: 419 - 4133 2012)表明，上部含水层的透水性(平均)是下部含水层的数倍。在SMEC调查期间(2010-2011年)，地下水位变化很小，这表明在此期间补给有限，或者是通过断层屏障的泄漏点控制地下水位。

这个泄漏点控制着最大的地下水高度，就像大坝一样。盆地北部较高的地下水位反映了渗透率下降和可能的补给。

概念模型
模型首先是对地下水系统知识的整合和综合，使人们能够深入了解地下水流系统的功能(例如，Bredehoeft 2005;Vandenbohede et al . 2011)。如图3所示，Unkheltseg盆地系统的概念模型从一个基本的理解，即一个充满沉积物和水的浴缸，发展成为一个概念性的二维地质模型，成为数值模拟的三维地质框架。盆地模型面临的挑战是可用的构造和地质信息有限。将勘探钻孔数据与地球物理极化测量的解释深度到基底相结合，绘制地质和地形图，以确定沉积物的厚度。利用钻孔数据和井内伽马测井来解释关键地质层，并确定上下含水层。

创建地质表面的方法与Danis(2012)的方法类似，即通过钻孔数据、地球物理数据和地表地质图的插值创建层。在复杂的岩性中，例如有许多不连续的沉积层，要创建最接近地质现实的三维表示，需要一定程度的简化以保持计算能力。因此，将盆地模型简化为由关键含水层、含水层和基底组成的地质单元。每一层都参考UTM WGS 84波罗的海高度基准面，并使用SurferTM version 9 (Golden Software Inc .)中的Kriging算法创建，在模型域上间隔50米。表面从Surfer导出到Hydro GeoBuilder(斯伦贝谢水服务2009)，在那里使用Kriging插值从不规则间隔的数据创建表面，然后导入MODFLOW-SURFACT。该模型由表3所示的四层组成。

数值模型

建模使用MODFLOW- surfact代码(hydrogeological Inc 2002)，一种先进的MODFLOW (McDonald and Harbaugh 1988;哈博等人

2000年)，基于hydrogeological Inc.开发的代码，处理网格细胞的完全去饱和和再饱和)，在Visual MODFLOW 4.6版本(斯伦贝谢水务公司2010年)的框架内。MODFLOW- surfact解决了执行无约束模拟的增强方程，以严格模拟含水层的去饱和/再饱和，并克服了以前版本的MODFLOW遇到的数值困难。此外，MODFLOW-SURFACT还提供了一个严格的井筒回收包“fracrewell package FLW4”，可以模拟多层井筒的流动。裂缝管表示允许连接含水层细胞。由代码自动计算与裂缝相关联的每个单独节点的体积通量，以计算井眼的总压降之和。此外，严格纳入了规定的孔眼总抽出率。FLW4的另一个特点是，对于超泵无承压系统，当井内水位达到井眼底部时，总井眼排量会自动调整。

调查的模型域是一个包含盆地的子区域(图1)。模型的实际面积是在模型开发的桌面阶段确定的，包括一个大约6公里乘7公里的矩形区域，总面积约为43平方公里。该模型主要用于地下水流的稳态和瞬态模拟。首先建立了一个稳态模型来评估盆地含水层的区域流动模式，并校准补给和透射率。通过试错法，通过应用一组水力参数完成校准;产生与实际现场测量相关联的计算机生成的压头(或水位下降)的边界条件和应力。然后利用稳态模型层建立了一个瞬态模型，表3模型层层描述厚度1泥沙沉积含水层达167 m 2在盆地东侧遇到的半封闭粘土沉积物达38 m 3在洪积沉积物中的半封闭含水层达88 m 4基岩，基于实测抽水试验数据、观测数据和12个月的监测数据，采用水头输出作为初始水头模型，评估抽取/抽水对地下水流型的潜在影响。

域与离散化

模型域的范围是基于适当的特定地点的地质和水文地质边界。模型域的设计是为了反映预计将显著控制主要地下水流量的含水层的范围。该模型包含119行乘144列，共68,544个单元格。在整个模型区域中应用了大约70 - 9 - 70米的网格，并在建议的泵井周围将网格细化到5 - 9 - 5米，以允许泵井附近更陡峭的水力梯度。在垂直方向上，采用变形模型层表示模型中的水文地质格架(图3c)。使用变形模型层的垂直离散化允许用更少的单元保持水平连续性(Reilly和Harbaugh 2004)。模型网格如图4所示。

模型区域位于一系列无流边界之间。在东部和西部，沿着地表划分，北部是由含水层的缺失或变薄推断的，南部是与基底高度一致的。为了模拟地下水系统的实际目的，采用了一种通用的边界制度，分别基于使用上游和下游的恒定水头和排水高度。在第3层活动模型单元的模型区北边界处，应用1735 m的恒定水头边界条件。排水标高为1,710.8 m，校准电导值为100 m2 /天，代表模型区南部边界的地下室高层的所有层都被分配。最低层(基底层)渗透性很低，被认为是不活跃的细胞。模型的下无流边界为基底底部，假设基底渗透率低，流动最小。

排水方案

DRAIN包(McDonald and Harbaugh 1988)被设计用来模拟农业排水和地质屏障等特征的影响。Drain包用于模拟与磁头相关的通量边界(Harbaugh et al . 2000)。在DRAIN包中，如果单元中的水头低于某一阈值，则从漏极到模型单元的通量降至零。假设排水管只部分满溢，因此排水管内的水头大约等于排水管标高的中位数di, j, k [L]。由单元(i, j, k)的模型计算的水头hi, j, k [L]实际上是单元的平均值，通常假设在距离排水口本身一定距离处普遍存在。排水管头di、j、k仅在排水管内局部存在;它不能作为一个整体来描述细胞的特征。

排水管以含水层水头与排水管标高之间的差成比例的速率从含水层中排出水。只要含水层的水头高于排水高度，排水就会继续，但如果水头低于该高度，排水就会停止。DRAIN包的功能由方程对描述:(McDonald and Harbaugh 1988)

式中QDi, j, k [L3 /T]为地下水从单元(i, j, k)排入排水口的排水量通量;hi, j, k [L]是计算出的cell (i, j, k)中的head;di, j, k [L]为排水标高。系数CDi, j, k [L2 /T]是一个集总电导，描述了漏极和电池(i, j, k)区域之间的所有水头损失，其中假设水头hi, j, k占上风。这种集总排水电导被添加到由块中心流(BCF)包计算的电导中，该包来自排水单元节点中心与相邻节点中心之间的水平和垂直阻力。因此，漏极和BCF电导的总和定义了MODFLOW漏极电池的总流动阻力。排水口与其相邻单元之间的水头损失是由流向排水口的会聚流、流过排水口的回填材料以及流过排水口壁造成的(McDonald and Harbaugh 1988)。

所采用的含水层性质要么基于水文地质调查(抽水试验)的结果，要么基于文平等(2011)的数据。在没有数据的情况下，使用行业公认值和最佳估计值。所采用的含水层性质包括导水率、储水系数和孔隙度如表4所示。在非平衡条件下，采用Cooper-Jacob直线法和距离递减法对含水层参数进行了估计。

它们提供了单个井眼和区域视角的透射率和储存系数的评估。结果总结在表2中。由于每个含水层的流量速率未知，根据2010年抽水试验结果(SMEC 2011)，假设每个含水层的流量为4升/秒，或者下部含水层的流量为1升/秒，上部含水层的流量为7升/秒。使用Aqtsolv Pro 4.0软件进行的评估表明，基于曲线匹配，每个含水层的流量分别接近1 L/s和7 L/s(即使用这些流量的曲线匹配优于使用每个含水层的4 L/s)。表2给出了假定承压和泄漏承压含水层的多重评价结果。

对于所有沉积物，垂直导电性设为水平导电性的1/10，以反映沉积沉积中分层引起的垂直各向异性的影响，以及层状非均质性(Freeze and Cherry 1979;齐米。2012)。各向异性表明，断裂带内的裂缝或角化物质必然存在某种优选方向，从而导致沿断层的优异导电性。在垂直于断层面的方向上，电导率要低得多。从地质学上讲，带断层泥的断层角砾岩可能表现出各向异性，就像断层内部的粘土一样(Lopez和Smith 1996)。对于地下室，考虑垂直方向的地下水流量大于水平方向;因此，垂直导水率设为水平导水率的10倍。通用存储术语已用于所有围合层和半围合层。在崩积层和洪积层的情况下，储存条件是从SMEC(2011年、2012年)抽水试验结果或可比岩性的已发表文献值中获得的。为了模拟基底含水层补给，模型域盆地段西部的灰岩露头渗透率相对较高，为0.06 m/d。

补给、排放和蒸散

根据地表地质特征和地下水深度，假设模型区上层地下水位的净补给量为每年0.9 mm，占降水的1%。采用2010 - 2012年逐月随时间变化的补给率(假设为逐月降水的1%)，对降雨进行暂态运行。考虑到地下水的深度，地表以下15米处的蒸发蒸腾不被认为是活跃的。

稳态标定

通过将观测到的地下水水头与模型域内20个钻孔的预测水头进行匹配，对模型进行了校准(图5)。为了改进校准，调整了水力传导率和补给区，直到模拟的水头高程能够与观测到的水头高程匹配到可接受的精度水平。

校准水头与第1层和第3层的观测水头的相关系数为0.98，均方根误差为0.45 m，残差平均误差为0.11 m，被认为是非常好的。基线梯度和流量沿含水层呈南向和西南向。基于这些结果，认为该模型较好地模拟了整个模型域中地下水的高程和流向。

进入盆地含水层的水量主要通过降雨补给，部分通过基岩流补给。模拟表明，流动系统保存质量(表5)，即通过补给和流动进入模型的水量等于通过排水离开模型的水量，代表地下室的高度。断层屏障处地下水位为海拔1710.8 m，标定电导值为100 m2 /day。

总体而言，基于相关系数(0.98)、质量平衡差异(0.02%)和空间随机残差(Barnett et al . 2012)，模型校准被认为是可以接受的。

瞬态标定与仿真

校准:观测头

瞬态校准的目的是在模型和观测到的水头压降之间实现±0.5 m的差异。陷落的波动主要是由于断层屏障处的补给、排水和应用定水头条件的波动引起的。

校准:泵送试验数据

SMEC(2012)进行了抽水试验，以评估含水层的水力传导性和蓄水系数，用于瞬态模型校准。现场泵送测试进行了10天，其中一个生产井眼的平均流量为8l / s，使用了10个观测压力表。这里的模型模拟了26天的水位下降，有10个时间步长，时间步长乘数为1.2。在模型中手动调整含水层参数，以匹配观察到的短期瞬态事件的下降。模拟的下降与观察到的下降的结果表明，大约10天后有一个很好的匹配，但最初的匹配不是很好，这可能是存储属性的结果。根据Middlemis(2004)的说法，不应期望区域模型准确地再现抽水对水头响应的局部尺度变化。

任何模型的价值都来自于分析参数对变化的敏感性，以观察这种变化对模型结果的影响。对补给和渗透率的敏感性进行了评估，结果如图6和7所示，以米为单位绘制了平均绝对误差。对变化最敏感的是第1层和第2层的水力导电性，如预期的那样，是补给层。

排水电导、比产量和比库存量对其敏感性较低。

模拟

从2000年开始到2012年结束，进行了7年的瞬态模拟，以评估矿井抽取地下水的影响。采用月降水1%的月时变补给率和表4的水力传导系数和蓄水系数值。方案如下:3口井旋转抽油7年，抽油速率从4.2升/秒到17升/秒不等(表6)，以满足矿山需求和抽油后20年的回收期。

通过瞬态运行模拟了盆地上下含水层的地下水抽取，分别预测了第1层和第3层在抽水后7年(图8)和20年(图9)后的地下水位和地表压降。

模拟结果:7年抽水结束时，监测压力表模型单元内的最大水位预测降为72 m左右。经过7年的泵送，泵送钻孔的影响半径约为5公里(图8)。注意，对于第一层，模式区中部和东北部的单元由于泵送而变干。对于第三层，只有东北部的几个单体变干。如图8所示，经过7年的抽水，第1层(崩塌含水层)的最大扬程下降约为20米，而位于盆地南部的Unkheltseg钻孔(当地的一个钻孔)的总扬程下降估计约为13米。

恢复后的等高线(图9)显示，即使在20年后地下水位也不会完全恢复。

对于第1层，在20米的下陷中，采收率约为13米，在20年后采收率为65%。

数值计算结果表明，以矿井平均需水量为9 L/s (2.83 9 105 m3 /年)计算，抽水7年后的水文地质影响将使崩落上含水层和洪积下含水层的地下水位分别下降20 m和72 m。

模型结果表明，抽水导致地下水储水量减少，20年后浅层和深层含水层的长期水位分别恢复65 - 85%。这表明由于渗漏到下面的含水层(第3层)，水从第1层垂直移动;对扎夫汗河基流的影响可能可以忽略不计。

该模型包括基于文献值和建模者经验的假设。

存在以下方面的局限性:•含水层补给速率和含水层补给区的空间分布;

•石灰岩和基底岩单元的水文地质性质;•构造对地下水流动的影响(断层、裂缝和节理)。

目前的模型假设降水补给有限，占降水的1%，流量由地质控制，有一个固定的反转点，下伏和周围基岩的影响有限。位于盆地西北侧的大型石灰岩可以提供大量的水。然而，只有长期抽水和监测才能评估这些因素的影响和所采用模型参数的保守程度。地下水补给很少得到很好的定义(Voss 2011)。初始水头数据与用于校准瞬态模型的观测时间跨度有关。由于没有较长期的数据，因此使用了两个水文年的数据;这可能是不确定性的来源，因为长期波动可能没有被观察到。

The geological data and measurements used in creating this model are spatially limited within the Unkheltseg Basin and this should be taken into account when assessing the results.

用于创建该模型的地质数据和测量在空间上受限于Unkheltseg盆地，在评估结果时应考虑到这一点。

结论在蒙古和中国的大型干旱或半干旱内陆盆地中，逆断层一般与山走向平行，垂直位移显著，使第四纪含水层系统的地下水位在断层两侧相差100 m以上。在断层两侧没有直接的水力连接，地下水在断层上落下。下游含水层系统有足够的能力接受来自上游侧的侧向流量，这意味着下游含水层系统的补给等于上游含水层系统的流量。排水量或回灌量仅由上游含水层系统的渗透率和水力条件决定。上游含水层系统的模拟(以及由断层形成的地下水)已成为区域地下水模拟中的一个有趣的挑战。Unkheltseg盆地是位于蒙古北部Taikh山谷的不寻常的地质限制地下水资源的一个例子。这种资源对该地区的经济发展至关重要，但以前的评估仅限于对抽水试验数据的分析解释。由于两部分盆地系统之间没有直接的水力连接，从Unkheltseg盆地到Tost盆地的地下水补给在地质上是通过断层控制的，只有当Unkheltseg盆地的地下水位达到一定水平时才会发生。

为了在区域范围内有效地管理地下水，需要进行校准的三维数值模拟，以评估资源满足矿山需求的能力和评估这种开采活动的潜在影响。编制一个经过校准和验证的地下水模型可以提高对控制水头和流量的主要物理过程和力的概念模型的信心。模型敏感性分析是确定水文地质系统最关键参数以确保适当的现场监测的有用工具。在这种情况下，该方法确定了需要明确的空间离散化和充分表示Unkheltseg盆地南端的水平和垂直流动方向，作为限制补给面积和评估和量化地下水动态所需的关键因素。改变应力，如增加补给，对模拟地下水位的影响很小，因为控制因素是跨越断层屏障的阈值流量水平。

准确预测区域影响的区域地下水模型是蒙古等高度依赖地下水的国家的重要管理工具。地下水模型是证明对与现有数据一致的水文地质系统的定性和定量理解的重要工具。该模型表明，短期高强度提取会在很长一段时间内影响大面积。虽然巴彦埃拉格矿是目前唯一一个计划从该流域系统取水的矿山，但地下水资源已经有了储备和国内用户。该模型展示了现有用户将如何受到影响，以及未来潜在用户也可能受到影响的地方，强调了适当的长期含水层管理战略的重要性。

References Allen DM, Michel FA (1999) Characterizing a faulted aquifer by field testing and numerical simulation. Groundwater 37(5):718–728 Barnett B, Townley LR, Post V, Evans RE, Hunt RJ, Peeters L, Richardson S, Werner AD, Knapton A, Boronkay A (2012) Australian groundwater modelling guidelines, Waterlines report.

National Water Commission, Canberra Battumur (2009) Assessment of water supply options for proposed Bayan Airag gold–copper mine Bense VF, Gleeson T, Loveless SE, Bour O, Scibek J (2013) Fault zone hydrogeology. Earth Sci Rev 127:171–192 Beven KJ, Kirkby MJ (1979) A physically based variable contributing area model of basin hydrology. Hydrol Sci Bull 24:43–69 Bidwell VJ, Stenger R, Barkle GF (2008) Dynamic analysis of groundwater discharge and partial-area contribution to Pukemanga Stream, New Zealand. Hydrol Earth Syst Sci 12:975–987 Bredehoeft JD (1997) Fault permeability near Yucca Mountain.

Water Resour Res 33(11):2459–2463 Bredehoeft JD (2005) The conceptualization model problem-surprise.

Hydrogeol J 13:37–46 Cello G, Gambini R, Mazzoli S, Read A, Tondi E, Zucconi V (2000) Fault zone characteristics and scaling properties of the Val d’Agri fault system (Southern Apennines, Italy). J Geodyn 29(3–5):293–307 Danis C (2012) Sydney–Gunnedah–Bowen Basin deep 3-D structure.Explor Geophys 43:26–35. doi:10.1071/EG11043 Doherty J, Simmons TC (2013) Conceptual modelling in decision support; reflections on a unified conceptual framework. Hydrogeol J 21:1531–1537 EcoTrade (2007) Bayan Airag Environmental and Social Baseline Report, Zavha Aimag, Southern Mongolia Evans JP, Forster CB, Goddard JV (1997) Permeability of faultrelated rocks, and implications for hydraulic structure of fault zones. J Struct Geol 19:1393–1404 Fetter CW (1994) Applied hydrogeology, 3rd edn. Prentice Hall, Upper Saddle River 

Fiorillo F (2011) Tank–reservoir drainage as a simulation of the recession limb of karst spring hydrographs. Hydrogeol J 19:1009–1019 Folch A, Mas-Pla J (2008) Hydrogeological interactions between fault zones and alluvial aquifers in regional flow systems. Hydrol Process 22:3476–3487 Freeze RA, Cherry JA (1979) Groundwater. Prentice Hall Inc, Upper Saddle River Gupta HV, Clark MP, Vrugt JA, Abramowitz G, Ye M (2012) Towards a comprehensive assessment of model structural adequacy. Water Resour Res. doi:10.1029/2011WR011044 Harbaugh AW, Banta ER, Hill MC, McDonald MG (2000) MODFLOW-2000, the US Geological Survey modular ground-water model—user guide to modularization concepts and the groundwater flow process, Open File Report. OFR00-92, US Geological Survey, Denver, Colorado, pp 121 Hasiniaina F, Zhou J, Guoyi L (2010) Regional assessment of groundwater vulnerability in Tamtsag Basin, Mongolia using drastic model. J Am Sci 6(11):65–78 Huntoon PW, Lundy DA (1979) Fracture-controlled ground-water circulation and well siting in the vicinity of Laramie, Wyoming.

Groundwater 17(5):463–469. doi:10.1111/j.1745-6584.1979.

tb03342.x HydroGeoLogic Inc. (2002) MODFLOW-SURFACT: a comprehensive MODFLOW-based hydrologic modelling system: code documentation and user’s manual. HydroGeoLogic Inc., USA Kampf SK, Burges SJ (2007) A framework for classifying and comparing distributed hillslope and catchment hydrologic models. Water Resour Res 43:W05423. doi:10.1029/2006WR005370 Kazemi GA (2012) Hydrogeology—a global perspective. InTech publishing company, Rijeka. ISBN: 978-953-51-0048-5 Kresic N, Mikszweski A (2014) Hydrogeological conceptual site models data analysis and visualization. Environ Earth Sci 72:623. doi:10.1007/s12665-014-3268-8 Lopez DL, Smith L (1996) Fluid flow in fault zones: influence of hydraulic anisotropy and heterogeneity on the fluid flow and heat transfer regime. Water Resour Res 32:3227–3235 McDonald MG, Harbaugh AW (1988) A modular threedimensional finite-difference ground-water flow model. USGS Techniques of Water-Resources Investigations Book 6, Chapter A1 Middlemis H (2004). Benchmarking best practice for groundwater flow modelling. The Winston Churchill Memorial Trust of Australia Nelson ST, Mayo AL (2014) The role of interbasin groundwater transfers in geologically complex terranes, demonstrated by the Great Basin in the western United States. Hydrogeol J 22:807–828 Nordqvist AW, Tsang YM, Tsang CF, Dverstorp B, Andersson J (1992) A variable aperture fracture network model for flow and transport in fractured rocks. Water Resour Res 28(6):1703–1713 Nordstrom DK (2012) Models, validation and applied geochemistry: issues in science, communication, and philosophy. Appl Geochem 27:1899–1919 Palma HC, Bentley LR (2007) A regional scale groundwater flow for the Leon-Chinandega aquifer, Nicaragua

Pimentel ET, Hamza VM (2014) Use of geothermal methods in outlining deep groundwater flow systems in paleozoic interior of Brazil. Hydrologeol J 22:107–128 Pulido-velazquez MA, Sahquillo-Herraiz A, Camilo OJ, Pulidovelazquez D (2005) Modelling of stream–aquifer interaction: the embedded multi reservoir model. J Hydrol 313:166–181 Ran G, Eyal S, Yoseph Y, Amir S, Noam W (2014) The permeability of fault zones: a case study of the Dead Sea rift (Middle East).

Hydrogeol J 22:425–440. doi:10.1007/s10040-013-1055-3 Reilly TE, Harbaugh AW (2004) Guidelines for evaluating groundwater flow models. US Geological Survey Scientific Investigations Report 2004-5038, 30 p Rojstaczer S (1987) The local effects of groundwater pumpage within a fault-influenced groundwater basin, Ash Meadows, Nye County, Nevada, USA. J Hydrol 91:319–337 Rupp DE, Schmidt J, Woods RA, Bidwell VJ (2009) Analytical assessment and parameter estimation of a low-dimensional groundwater model. J Hydrol 377:143–154 Schlumberger Water Services (2009) Hydro geobuilder. A flexible, simulator-independent, hydrogeological modelling environment.

User’s manual. Waterloo, Ontario, Canada Schlumberger Water Services (2010) Visual modflow. Dynamic groundwater flow and contaminant transport modelling software.

User’s manual. Waterloo, Ontario, Canada Shan S, Javandel I, Witherspoon PA (1995) Characterisation of leaky faults: study of water flow in aquifer-fault-aquifer systems.

Water Resour Res 31(12):2897–2904 Simmons CT, Hunt RJ (2012) Updating the debate on model complexity.

GSA Today 22(8):28–29. doi:10.1130/GSATG150GW.1 Singh VP, Woolhiser DA (2002) Mathematical modelling of watershed hydrology. J Hydrol Eng 7:270–292 Sloan WT (2000) A physics-based function for modelling transient groundwater discharge at the watershed scale. Water Resour Res 36:225–242 SMEC (2011) Bayan Airag Gold Project pre-feasibility level mine water supply assessment. Bayan Airag Exploration LLC SMEC (2012) Bayan Airag Gold Project feasibility level mine water supply assessment. Bayan Airag Exploration LLC Vandenbohede A, Hinsby K, Courtens C, Lebbe L (2011) Flow and transport model of a polder area in the Belgian coastal plain: example of data integration. Hydrogeol J 19:1599–1615 Voss CI (2011) Groundwater modelling fantasies—part 1, adrift in the details. Hydrogeol J 19:1281–1284 Wenping L, Zhenying L, Haipeng G, Ning L, Weidong K (2011) Simulation of groundwater fall caused by geological discontinuities. Hydrogeol J 19:1121–1133 Yihdego Y, Becht R (2013) Simulation of lake–aquifer interaction at Lake Naivasha, Kenya using a three-dimensional flow model with the high conductivity technique and a DEM with bathymetry. J Hydrol 503:111–122 Yihdego Y, Webb AJ (2014) Use of a conceptual hydrogeological model and a time variant water budget analysis to determine controls on salinity in Lake Burrumbeet in southeast Australia.

Environ Earth Sci. http://link.springer.com/article/10.1007/ s12665-014-3509-x Environ Earth Sci (2015) 73:4119–4133 4133