在蒙古国Unkheltseg盆地的地质不连续性的三维数值地下水流模拟

在评估和采矿规划中，从勘探阶段到关闭阶段，高效利用干旱地区的水资源变得非常重要。

数值地下水模型的目标是评估从世界上最干旱地区之一提取地下水以满足矿井用水需求的潜力。

数值地下水模型用于评估地下水资源。

使用MODFLOW-SURFACT代码进行了建模，这是一种基于MODFLOW的先进代码，嵌套在Visual MODFLOW 4.6版本框架中。

建立了一个稳态模型，用于评估区域地下水流模式，并在模型中校准了补给和水力导数参数。

该模型的校准系数为0.997，均方根误差为0.3米。

采用瞬态模拟模型预测了在主要含水层水文地质体系上抽取150万立方米/年水量的影响，包括停泵后的情况。

建模模拟了四种水文地质情景。对“最坏情况”的模型结果表明，Saq砂岩含水层应该能够满足矿井用水需求（150万立方米/年）的10年。

然而，长期的水位下降显示持续减少而未达到稳态条件；因此，大部分水都是从含水层储存中取出的，在长期内将会有来自模型区域北部钻孔场的相互干扰。在该区域，水力梯度相对陡峭，超采已超过28年。

另一种情景显示，在停泵18年后，水位下降的11.6米中将有大约8米的回升，这意味着含水层将受到压力，大部分水将从含水层储存中取出。

为了最小化水文地质影响，建议将生产钻孔分散布置，应将钻孔设在Saq砂岩最深的地方，并以较低的抽水率操作钻孔。

总体而言，该研究提供了有关矿井供水评估的有用数值模型输出，并将有助于在世界上最干旱地区的不同条件下改善水资源管理。

导言

水是矿业发展的关键组成部分，从人类饮用的饮用水到矿产回收的工业用水。需要一种更全面、基于物理的方法，利用场地特定的数据。近年来，在不同环境中使用了各种程度的复杂性的水文模型，以评估在各种气候或土地利用情景下的水预算（例如，Turner等人2015；Whittier等人2010；Yihdego等人2015）。总体而言，水预算模型用于估算补给量，这作为地下水模型的输入。

可持续性指标，如可接受的水头或水质水平，是基于影响水提取和土地利用的气候和管理决策来考察的。由于模型和用于模型校准的数据的限制，可持续性指标的现场监测应是任何可持续性规划工作的重要组成部分。使用数值流模拟的资源评估程序说明了将含水层的结构复杂性（例如，Mair等人2013；Yihdego和Webb2015；Yihdego和Becht2013）纳入考虑的重要性，并在评估地下水产量时，这是使用传统的分析方法无法实现的一种方法（Gingerich2008；Sophocleous2000）。在应用不受现场数据支持的简化假设的情况下，应用分析解决方案会出现困难。数值建模用于评估地下水潜力。在这方面，数值模型作为诊断、管理和预测地下水行为的工具在近年来变得越来越重要。

在研究区域——沙特阿拉伯中部，与许多半干旱地区一样，详细评估地下水资源尤为重要，因为城市和农业发展以及矿业项目对可靠地表水资源的需求不断增加（Al-Salamah等，2010年；Hagedorn等，2011年；Won等，2005年）。

在干旱区域，如研究区域，由于勘探阶段至结束，对水资源的理性和高效利用在评估和矿业规划中变得更加相关。本研究提出了在世界上最干旱地区之一为矿山供水设置数值模型的方法论。数值水文地质模型的目标是评估从跨境地区的Saq砂岩（“目标含水层”）提取地下水的潜力，位于沙特阿拉伯中部Khnaiguiyah附近。为了评估从Saq砂岩向矿山提供每年150万立方米（MCM/年）的水需求在10年内的可行性，开发了一个区域三维数值地下水模型。这项工作包括准备概念性水文地质模型以及数值地下水模型的校准和计算机运行。本文提供的数据和解释来自SMEC（2012年）的顾问报告。

区域背景

研究区域位于沙特阿拉伯中部的阿尔利雅德省，距离首都利雅得约170公里。矿山研究区域位于Al-Quwayiyah省，外部取水区域和输水管道范围延伸到邻近的Ad-Duwadimi省。

该地形多山，海拔从租赁区域下游的800米到流域顶部的900米不等。植被覆盖很少。沙特阿拉伯大部分地区通常属于干旱地区，年平均降雨量在50到100毫米之间变化。沿红海悬崖和利雅得以北地区降雨较多，年平均降雨量分别达到600毫米和200毫米。

在利雅得测得的蒸发量年均为2900毫米，夏季每日峰值值可能超过18.5毫米/天。Abunayyan和BRGM（2008年）使用Penman方程从五个区域气象站估算了蒸发蒸腾量。计算得出的年均蒸发蒸腾量为2376毫米/年。因此，年降水量和蒸发蒸腾量之间存在相当大的差异，约为90.3：2376毫米，比值为1：26.3。

水文地质学

沙特阿拉伯地处阿拉伯板块上。在区域上，该板块可以分为两个地质体系，即阿拉伯盾和阿拉伯平台（图1）。它由晚前寒武纪变质火山岩和火山沉积岩组成，被火山岩侵入。阿拉伯盾基底的倾斜导致这些沉积岩向东倾斜，朝着阿拉伯湾。商业上重要的天然气和油田位于这些沉积岩层在利雅得和阿拉伯湾之间的区域。Saq砂岩是一个地区性广泛分布的地质单元，横跨沙特阿拉伯中部到北部的一半，从约旦延伸到叙利亚，距离约1200公里，表面面积约为6.5万平方公里。Saq砂岩与上覆的二叠纪Khuff组之间有不整合。距离Khnaiguiyah项目最近的主要含水层系统是地区性广泛的Saq砂岩。作为一个地区性广泛的含水层，它储存了大量的地下水（约2800亿立方米）。其最南端距离Khnaiguiyah约10公里，那里它与一处古基底高地相对。Khuff组是一个较不重要的不连续含水层系统。它延伸到矿山研究区域附近。

Saq砂岩水位下降最严重的地区集中在卡西姆和塔布克地区，分别为140、170和210米。水位下降与不透水的基底露头平行。由于抽水和地下水减少，压力面正在缓慢地区域性恢复（图1）。基线文献报告了一个渗透率范围为200到7000平方米/天，渗透率为6到70米/天，孔隙度为20到25%。该含水层在西侧是无限制的，在Khuff组的Ash Shiqqah不透水层下方是有限制的。单个地下水产量通常约为50到100升/秒。

自然地下水流方向是向东北，除非由于主要水位下降导致地下水水力梯度的反转。

历史数据

假定矿山对水的需求为每年1.5百万立方米，为期10年。目前，这样的供水量无法通过废水处理设施、地表水或蓄水坝获得。从Saq砂岩中可能长期获得大量地下水供应，用于农业和城镇供水，包括AlQuwayiyah和Labhka水源地。有限制的Saq砂岩含水层的渗透率在100到200平方米/天之间（1.2到2.3×10^−3 m^2/s），渗透率为1到3米/天，储存系数为10^−4到10^−5。评估的含水层渗透率处于其他地方报告的下限。从个别生产井中可开采的产量在20到50升/秒之间。水盐度范围从1360到17,000毫克/升，后者被认为是来自上覆的Khuff组和排放阿拉伯盾东缘的浅层冲积层的垂直渗漏（AlSalamah等2010年；Abunayyan和BRGM 2008年）。

概念水文地质模型

概念水文地质模型（CHM）是发展数值模型中描述场地水文地质条件的关键先决条件。CHM提供了有关地质、水文地质、地下水的出现和流动、水文化学和当前地下水利用的解释信息。该模型涵盖了约357平方公里的区域（“模型区”）。外部水源区位于模拟区域的西南（图2）。

模型的西、东和南边界是根据地质特征和含水层范围选择的。 Al Quwayiyah水源区和Labkha水源区分别位于北边界的北部和东北部。

地质和含水层

研究区域的地质和水文地质包括第四纪沉积物的地表、风化、冲积和堆积的砾石；波斯湾盖层，在其上面是二叠纪Khuff组，不整合；骆驼带沙砾岩至寒武纪，不整合；和前寒武纪基底阿拉伯盾复合体。

Saq砂岩是矿山供水的主要目标。 Khuff组的下部包含Ash Shiqqah单元，由页岩和块状的砾岩组成。上部包含含有化石的泥灰岩、白云石、砂岩和石膏。Khuff组是一个较不重要的不连续含水层。使用钻井记录、地质/水文地质剖面、已发表的地质图、地球物理数据和地质露头图，使用Hydro GeoBuilder 2009版本的软件包制作了一个三维水文地质模型。位于外部水源区下方的Saq砂岩位于前寒武纪基底高地上。

含水层水力参数

采用的含水层性质基于当前水文地质调查的结果以及该地区以前的研究（Al-Salamah等人，2010年）。

在没有数据的情况下，采用了文献值（Fetter 2001；Freeze和Cherry 1979）和专业判断。采用了采用的数据进行模型校准。校准后的含水层参数总结在表1中。

地下水流系统

地下水流系统存在于所有地形中。

运动发生在以下地区：

& 补给区，水可以进入地下水系统；

& 侧向流动区，可以在大型盆地中非常广泛；和

& 地下水排泄区，地下水通过自然或人工抽水系统重新出现在地表。

流动系统可以是区域性、中等规模或局部规模的，并且在很大程度上受到所考虑地区的地形和地质结构的驱动。Saq砂岩是一个区域性含水层系统，长达1200公里。外部水源区位于该区域性含水层系统的南端。Saq砂岩通常朝东倾斜。含水层的补给发生在露头区域。含水层的排泄主要来自灌溉和市政生产井以及地下径流。Khuff组范围较小，但遍布数值模型的整个区域。

Saq砂岩在外部水源区以西和以南几公里处露头（图2）。模型区域内地下水位和流动模式的潜在主要影响因素有： & 地下水开采，从Saq砂岩中用于灌溉和城镇供水。例如，Al Quwayiyah水源区自1984年以来一直在运营； & 地形，沿Saq砂岩的露头；和 & 降雨，上游集水区和露头处。

Saq砂岩的补给发生在露头区域，通过短暂的地表水系统进行垂直渗透。

从区域尺度上看，地下水提取量的体积与含水层的补给量不相容。因此，在许多地方，含水层正在被“开采”。地下水开采导致区域地下水位稳步下降。

压力面

初始的灌溉前水头是从19个数据点插值得出的，其中的水位是根据美国地质调查局的数字高程模型（DEM）估算的（图3）。从灌溉和市政生产井中提取地下水对2012年的影响（基于69个测量点）如图4所示。在过去的28年中，地下水流动模式和水力梯度发生了显著变化。在Al Quwayiyah水源区观察到最大的水位下降（每年4.2米）。在外部水源区域没有发生大规模的地下水开采。在RL + 718米的位置似乎存在人工地下水分水岭，该地下水分水岭位于模型区域的北部（图4）。自2008年以来，外部水源区Saq砂岩中的地下水位保持相对稳定（水位深度∼距离地表61米），表明这个区域受到从位于至少23公里以北的水源区的生产井大量抽取的影响较小。

地下水流向

在进行大规模灌溉和市政开采之前，地下水流动方向呈东南方向，基本符合地表地形（图3）。模型区域内的地形梯度为0.005；自然地下水水力梯度为0.0009（20米的高程在大约23公里的流动路径距离上）。2012年的地下水等值线图显示了在外部水源区域以北的显著水力梯度反转（图4）。

例如，在Al Quwayiyah水源区，发展出了陡峭的水力梯度为0.02（在10公里的距离上失去了约200米的水头），而平均模型区域的水力梯度为0.0005（在10公里的距离上失去了5米的水头）。

补给和排泄

Saq Sandstone地下水流系统的补给涉及露头上的降雨和地表水流。Abunayyan和BRGM（2008）估计了降雨渗透速率为1到5％。

已定义了补给区域。这些基于地表地质单元，自然补给估算为降水的百分比。

Saq Sandstone被分配为年均降水量的5％。

长期年均降水量假定为110毫米/年。由

于缺乏数据并且为了保守，未向Khuff Formation或通常干燥的地表冲积沉积物分配补给。

对于瞬态状态模型，应用了随时间变化的补给，作为平均每月降水量的百分比（5％）。

来自上游地区的补给贡献是使用每月降雨计算的，降雨补给为0.1％。应用于露头的Saq Sandstone形成的降雨补给详见表2。已评估上游集水区的面积分别为77平方公里（上游补给区域1）、46平方公里（上游补给区域2）和45平方公里（上游补给区域3）。这些显示在图5中，位于模型区域的西部和北部，并表示为0.1％的降雨补给，分别产生22、13和13立方米/天。

地下水产量

从模型域中Saq Sandstone含水层的大部分排泄发生在通过流和水源（用于灌溉和市政提取）。

地下水开采未经测量且文献记录不足。已进行了地下水普查（SMEC 2012）。在外部水源区域以北约20公里处提取了大量地下水（图2）。例如，从Quwayiyah水源区和Labkha水源区分别提取了8500和2000立方米/天。

大多数大型灌溉井每年连续供水长达5个月，每天高达7000立方米。在模型区域北部有109个这样的灌溉井。在外部水源区域附近，从Saq Sandstone含水层提取的灌溉摘要相对较小，估计约为360立方米/天，用于五口小直径抽水井（平均72立方米/天）。这些用于灌溉枣树并用于向当地社区和Quwayiyah城镇销售的卡车装载。

数值建模

水文地质数值模型已经成为与矿业运营中的水有关的预测的经常使用的工具。水文地质数值模型的应用包括矿业供水的地下水流模型。数值模型用于评估拟议矿山及相关工程和设施的地下水可用性。地下水流模型通常用于水资源评估，以评估含水层的长期生产力。本研究介绍了Saq Sandstone及其相关地下水系统的数值模拟的开发和应用。它包括模型校准和对位于外部水源区域下的Saq Sandstone进行长期地下水提取（1.5 MCM/年）的模拟。

建模使用MODFLOWSURFACT代码（HGL 2013），这是一个基于MODFLOW的先进代码（Harbaugh 2000），可以处理格网单元的完全饱和和重新饱和，该代码嵌套在Visual MODFLOW版本4.6的框架中。

模型设计

空间离散化和分层

模型涵盖了一个面积约为357平方公里的区域（模型区域），Al Quwayiyah和Labkha水源区分别位于模型区域的北部和东北部（图2）。该模型旨在稳态和瞬态地模拟地下水流动。

在模型区域采用了大约500×500米的网格，并在提议的水源区域周围对网格进行了20×20米的细化，以适应“外部水源”附近陡峭的水力梯度。在垂直方向上，使用变形的模型层来表示模型中的水文地质框架（图6）。使用变形的模型层进行垂直离散化可在减少单元数的情况下保持水平连续性（Reilly和Harbaugh 2004）。模型网格见图7。

该模型基于一个概念模型，包括对应于表3和图6所示的主要地质单元的三个层。分配了非均匀（100米厚度）的基底底层。这些层是Khuff Formation（第1层）、Saq Sandstone（第2层）和基底（第3层）。

边界和初始条件模型域的范围基于适当的现场特定地质和水文地质边界。模型域的设计旨在反映Saq Sandstone的范围，该范围预计将明显控制主要的地下水流。

模型区域位于一系列无流边界之间， & 西 - 沿着基底的露头（不透水）； & 南和东 - 推断出Saq Sandstone的缺失或减薄；和 & 北 - 与当前地下水分水岭线重合（图4）。

北边界被设置得足够远离提议的水源区（约10公里），以最小化来自北方地下水开采区域的影响。这个边界与海拔718米（asl）等高线以上的当前地下水分水岭形成的诱导水力边界重合（图4）。

补给条件

进入模型区域的上游补给量已转换为网格单元值的补给率，补给区分如图5所示。

建模方法

采用分阶段的方法进行了建模，以改进对CHM中采用参数的校准和细化。该模型用于模拟从生产井处进行四种情景，每年抽水1.5 MCM，持续10年（在外部水源区域），位于Saq Sandstone的相对筛选位置，

- 情景1（包括补给和泄漏 - 最佳情况），
- 情景2（包括补给但不包括泄漏），
- 情景3（无补给和无泄漏 - 最差情况）和
- 情景4 - 与情景3相同，但水位模拟延伸至抽水停止后的18年（恢复期）。

模型校准

关于MODFLOW参数的最小校准是为了测试这些值，特别是水力导数，是否独立于所使用的模型。校准主要与控制稳态和瞬态条件的参数有关。首先开发了一个稳态模型，以评估区域地下水流模式并在模型中校准补给和水力导数参数。

然后开发了一个瞬态模型，以评估抽取/抽水对地下水流模式的潜在影响，使用了抽水试验数据。

然后使用预测性模拟来评估含水层的水力响应。通过应用一组水力参数，边界条件和应力，产生计算机生成的模拟压力头（或水位降落），其匹配实际测量的场地内的水头在可接受的误差范围内。

模型校准通过手动（试错）和自动（反演模型，参数估计（PEST）; Doherty 1994）方法完成。

稳态校准

通过将观察到的地下水水头与模型区域内的九个预测水头进行匹配来校准模型（如图9所示）。为了改进校准，调整了水力导数和补给区域，直到模拟的水头高程能够以可接受的精度匹配观察到的水头高程。

图9呈现了校准结果（即目标点的校准头与观测头的对比），相关系数为0.997（均方根误差为0.3米）。该图显示了观察到的和计算的地下水水头之间的良好相关性。模型校准的残差误差范围在0.009米至0.674米之间。对校准稳态模型输出的分析表明，该模型以可接受的程度模拟了模型区域内的地下水高程和流向（比较图8和图4），鉴于信息有限，模型校准被认为是相当可接受的。这种偏差很可能是由于水力导数值在比模拟中表示的尺度小的尺度上的可变性引起的。然而，据信模型是合理校准的，这一点可由合理的头匹配和与研究区域的各种地质特征相关的水力导数值的空间分布证明。

瞬态状态校准

进行了一个抽水试验以评估含水层的水力导数和储存系数（SMEC 2012）。试验运行了48小时，使用了8.3 L/s的生产井排放（第一小时）和13到15 L/s（39小时）以及一个观测压力计（距抽水井50米）。该模型用于模拟水位降低。模拟时间为90小时，有10个时间步和时间步长乘数为1.2。 SMEC（2012）报告了这次抽水试验中遇到的困难。

使用手动和自动方法（PEST）在模型中调整了含水层参数，以匹配短暂瞬态事件的观测降水。对水力导数和储存系数的调整导致了如图10所示的匹配。曲线匹配在较晚的时间很好，但在较早的时间不太好；后者似乎是储存性能的结果。降低储存系数将对早期降低产生所期望的影响，但也将使其超过后来的观测值。不应期望区域模型能够准确复制由于泵取引起的局部尺度的头部变化（Barnet等，2012; Middlemis2004）。一旦完成了外部水源区域新建成的生产井的抽水试验，模型需要重新校准。

敏感性分析

使用PEST（自动估计方法）进行了敏感性分析，通过改变参数值并检查这些变化对模型结果（头部数据）的影响。绘制了缩放的加权敏感性与参数相比
- 水力导数（k），
- 储存（ss，sy）和
- 补给（par）

以比较它们对模型输出（头数据）的影响。考虑了层1和层2（Khuff Formation和Saq Sandstone）。

为此，从抽水试验期间监测压力计中选择了四个代表性的头部数据（命名为of000029、of000045、of000085、of000100），这些数据是在早期到后期的时间序列中观测到的。

敏感性分析的结果呈现在图11中。对于层1（Khuff Formation），储存（比降），其次是Khuff Formation和Saq Sandstone的水平水力导数（kx_1和kx_2）最为敏感。补给参数的敏感性较小。

模型结果

稳态模拟

表4提供了稳态条件下模拟的水量平衡概要。它显示模型是平衡的，进入Saq Sandstone的水量主要通过降雨补给。模拟表明，流系统保持质量。通过补给进入模型的水量等于通过开采和横向流离开模型的水量。目前从少数现有的生产井中抽取的地下水约为潜在含水层补给的24％。这反映在外部水源区域的饼状地下水面，位于地表以下61米处。

由于假设沿边界的无流条件作为分水岭，因此在高海拔区域的不准确性也是显而易见的。当然，正如MODFLOW的模拟情况一样，额外的校准可能会产生更接近的结果。然而，在没有清晰的参数调整基础的情况下，通过调整参数的空间分布来进行校准只是一个拟合过程。

瞬态状态模拟

进行了瞬态模拟，以估算项目地下水抽取对外部水源区域下方的Saq Sandstone的影响。该模型用于对未来10年进行预测，使用1995年至2005年的月变化补给率（即在露头区和高地区均假定了每月降水的百分比），其涵盖了在校准中获得的水力导数和储存系数的值。

情景包括：
- 情景1：10年内以1.5 MCM/年的速度运行（有降雨补给和来自Khuff Formation的泄漏），
- 情景2：10年内以1.5 MCM/年的速度运行（有降雨补给但无泄漏），
- 情景3：10年内以1.5 MCM/年的速度运行（无降雨补给和无泄漏）以及
- 情景4：在10年的抽水后不再泵取18年（恢复期）。

在模型单元格中包含监测压力计（距泵井50米）的水位降低在图12、13和14中呈现。预测的水位降低在10年连续抽水后的监测压力计模型单元格中分别为7.5、11.5和12米，对应于情景1、2和3。整个模型区域在模拟结束时（10年后）的累积体积预算显示在表5中。

情景1显示，以降雨补给和来自覆盖的Khuff Formation的泄漏的条件，水位降低最小，并且随着时间的推移开始恢复，因为由于泵取而形成的影响半径有助于产生大量的补给。情景1显示，在强加了来自上覆Khuff Formation的降雨补给和泄漏的条件下，水位降低最小，并且随着时间的推移开始恢复。

情景2和3表明，水位下降在抽水间隔内继续，而不会达到平衡/稳态条件。在情景2中，强加的条件下，项目地下水提取（1.5 MCM/年）超过了降雨补给（假设来自Khuff Formation没有垂直泄漏），并且地下水来自含水层存储。情景2和3之间存在轻微差异（分别在图13和14中），表明相对于泄漏而言，补给相对不重要。

在“最差情况”情景3条件下，监测压力计中预测的水位降低不会超过Saq Sandstone中的可用降低，并且含水层保持在受限的水文地质条件下。结果表明，Saq Sandstone应该能够在10年的持续时间内满足矿井的用水需求（1.5 MCM/年）。请注意，模型区域西部和南部的单元格由于抽水而变干（图15）。

MODPATH软件包（Pollock 1994）采用了粒子跟踪方法，使用水流路径来定义对地下水抽提的贡献的影响半径（图15）。流路径的确定有助于理解和量化由于泵井导致的模型区域水源的来源（向回跟踪颗粒到达的补给位置）。在情景3下，Outside Borefield的影响半径在设计后的10年抽水后约为10公里。形状反映了环绕Outside Borefield的多个屏障边界。

潜在影响评估

数值模拟表明，创建在井场周围的降水锥可能会延伸到北部模型区域边界之外，那里受到Quwayiyah和Labkha井场以及模型区域北部的灌溉井的干扰。随着时间的推移，可能会发生累积的水位下降。应定期在外部水源区和北部模型区边界附近进行水位监测，以随时间更新数值模型。

从外部水源区抽取地下水将影响周围井的水位行为。图16显示了在10年内以1.5 MCM/年的速度从外部水源区抽取地下水的模拟影响。相邻井中的水位可能下降约12米（井28）。

生产井的水位估计约为20米（图17）。但是，在外部水源区的井中并未进行多级测试。生产井中的水位降低取决于许多因素，包括生产井的数量、井间距、抽水速率、筛长、井直径和井效率。

情景4（图17）显示，在长期内将有约8米的恢复，超过了由于10年内1.5 MCM/年的抽水而引起的11.6米的水位下降。水的恢复将约为70％，这意味着含水层将受到压力，大部分水将从Saq Sandstone含水层储量中提取。

讨论和结论

已构建了一个三维数值模型，以允许模拟设计矿井水需求的地下水响应。该模型是基于有限的数据和水位测量而构建的。因此，它旨在提供目标含水层水供潜力的初步评估。随着新数据的出现，模型需要进行细化。稳态模型似乎已经取得了很好的校准。需要新的生产井和抽水试验以更好地校准瞬态条件。使用开发的数值模型，对这些单元的管理得出了一些重要的结论。

模型校准是根据水头数据进行的。尽管模型数据存在不确定性，但仍然获得了合理的匹配。稳态模型用于评估区域地下水流模式并校准模型中的补给和水力导数参数。该模型的校准系数为0.997，均方根误差为0.3米。瞬态模型用于评估地下水提取对含水层流动模式的潜在影响，使用抽水试验数据。使用瞬态模拟模型预测了以1.5 MCM/年的速度连续抽取10年对Saq Sandstone水文地质体的影响。

在最坏的情况下，即情景3（即无泄漏和无补给），监测压力计中的模拟水位下降表明，在外部水源区下方的Saq Sandstone含水层应能够在10年内满足矿井的用水需求（1.5 MCM/年）。然而，长期水位下降显示出持续减小而未达到稳态条件；因此，大部分水来自含水层存储；在Saq Sandstone终止的南部和西部边界处，含水层似乎会干涸；Saq Sandstone含水层位于基岩高地上；并且在长期内，将会从模型区域北部的Al Quwayiyah井场产生相互干扰。在该地区，由于超过28年的过度抽水，水力梯度相对较大。

在研究区域的情况下，由于补给减少，地下水储量缺口的增加是可预见的，尽管考虑到未来考虑的管理选项将加速耗竭过程（图14）。因此，在地下水管理规划中，应该考虑压力计水平的普遍下降，即使没有抽水，这样，限制性和不受欢迎的政策也无法在水位上产生积极影响。

与任何建模研究一样，应强调由于模型和数据限制而导致的不确定性。此外，校准是使用水头的一个快照完成的；因此，考虑到建模固有的非唯一性，这种校准应被视为不完整的。

情景4显示，在11.6米的水位下降中将有8米的恢复，这意味着含水层将受到压力，大部分水将从含水层存储中提取。为了更准确地评估目前的地下水状况和抽取对地下水的潜在影响，建议进行以下进一步工作，建造新的生产井，并进行多级和长期泵出试验；在Khuff Formation和Saq Sandstone含水层中建立一个嵌套的压力计，以定义模型区域内含水层之间的相互流系统。目前，尚未检查模型输出的垂直水力梯度的可靠性；并在各种含水层（Khuff Formation和Saq Sandstone）中建立一个地下水监测网络，以便将来更新已校准的模型。

为了最小化水文地质影响，建议将生产井横向分布（至少相距500米）；应将井位于Saq Sandstone最深的东部边界周边；并以低抽水速率操作井。

作为替代方案，可能需要增加井场的规模，即购买Outside Borefield北部更深的地块。

在含水层补给率和含水层补给区域的空间分布方面存在局限性，Khuff Formation和基岩单位的水文地质特性，对地下水流的结构影响（断层，裂缝和节理），浅层含水层的地下水监测以及Saq Sandstone和Khuff Formation之间的水力联系；除了最近的数据外，从北部抽水井中获得的关于地下水抽取的历史信息很少，有些水力参数值是从模型区域外的地点获得的（Al Quwayiyah和Labkha井场）。假定从这些抽水试验中获得的数据是具有代表性的，并且符合模型区域中预期的数据。KMC提供的含水层试验用于瞬态校准目的。当新的生产井建造时，应对模型进行额外的泵测校准；活动模型域北部边界的定义是通过使模型限制与过去28年的泵效果引起的地下水分水岭相一致来完成的。在Outside Borefield系统模拟和相互干扰的流交换之间仍然存在一定程度的不确定性。

本研究的结果为Outside Borefield区域地下水资源的可持续管理提供了指导方针。未来的管理决策应考虑到长时间持续的干旱和气候变化对可持续性指标的严重不利影响。在研究区域的情况下，应更详细地研究在气候变化条件下补给的重要性，因为这是这些系统中最重要的流入。因此，需要对引入的气候变化不确定性进行更全面的分析，以减小现有的不确定性，以改进概念模型、数值模型或参数值，从而获得可靠的支持决策的工具。