基于visu-almodflow的白涧矿区水文地质预测研究.docxVIP

基于visu-almodflow的白涧矿区水文地质预测研究 0 应用市场的心理计算平台 加利福尼亚州的威洛伊木马软件公司（威洛伊木马软件公司）根据美国地质秘书处（usgs）的软件模块进行了开发，并将可视化技术应用于云传播手段中的水循环软件模块。这是一个基于空间数据的模拟地下水流动值的软件。该系统在无缝集成MODFLOW-2000 、WinPEST、 MT3D、 MODPATH等软件的基础上,建立了合理的Windows菜单界面与可视化功能,增强了模型数值模拟能力、简化了三维建模的复杂性。其独特的可视化菜单结构允许操作者简单快捷的定义模拟区范围、选择单位、方便的分配模型的属性和边界条件,不需要对源程序进行任何的修改,Visual-Modflow就可以直接用来解决大多数地下水模拟问题。 河北是钢铁产业大省,邯邢地区则是重要的铁矿资源产地,自1950年以来,经过广大地质工作者的努力,该地区共探明铁矿石资源储量近9亿t,但是浅部矿床已基本采完,正在向水下的深部矿体开采。虽然经过了长期的探索实践,但目前主要的防治水措施还是“疏干排水”和局部注浆的“带压开采”。“大水”矿山的防治水仍然是本地区铁矿开发中的一大关键性难题,严重地制约了本区深部铁矿的开发。白涧铁矿是邯邢地区埋藏较深、规模较大的铁矿床之一,是典型受“大水”威胁的深部矿体,因此采用Visual Modflow预测模型对其进行矿井涌水量预测,对今后该地区深部矿床的勘查开发具有重要的指导意义。 1 含水层对矿床充水性的影响 白涧铁矿矿区主要含水层有第四系和煤系含水层及奥陶系中统灰岩含水层,其中:第四系及煤系含水层地层含水量小且有较厚的隔水层,对矿床充水影响不大。矿床充水岩层主要为奥陶系中统灰岩。 奥陶系中统灰岩含水层为矿床直接顶板,因火成岩侵入和构造的影响,使该含水层空间形态变化很大。造成灰岩含水层透水性、富水性在垂向上和水平上极不均匀。但就整个矿区而言,该奥灰含水层是百泉岩溶水地下水系统的一部分,是一个统一含水层,具有统一的地下水流场。 2 流量边界类型 白涧铁矿矿床赋存于-300~-800m水平,位于当地侵蚀基准面以下,属于埋藏较深的隐伏矿床。石灰岩含水层为矿床顶板,矿区处于百泉泉域南部北洺河—郭二庄—王窑—白涧—中关岩溶地下水强径流带上。模型原型以白涧铁矿矿区为对象,根据区内奥灰含水层流场特征和对地层结构的分析,将矿区侧向边界类型确定为不同性质的流量边界(图1),具体如下: (1)西北部边界:为各灰岩含水层隐伏露头,地表主要为冰川沉积杂色粘土、冰碛砾石层以及黄土等组成,阻隔了大气降水对下部含水层的补给。根据2007年奥灰含水层地下水渗流场特征,西北部位的等水位线与西北边界近于垂直,将其设置为隔水边界。 (2)东北部边界:地下水从王窑矿区经白涧矿区、中关矿区、西郝庄最终流向百泉,因此可以确定东北边界是矿区的排泄边界。 (3)东南部边界:矿区东南有大的侵入岩体,结合地下水渗流场特征,将其设置为隔水边界。 (4)西南边界:区域外有马会河流过,且是北洺河—郭二庄—王窑—白涧—中关岩溶地下水强径流带上的补给区,因此将其设为补给边界。 3 地下水系统的概括和含水层参数的选择 3.1 复配介质的性质 矿区内地下水流整体呈自西南向东北流动,将从以下几个方面进行本区地下水系统概化: (1)含水层水文地质参数随岩性的不同而异,方向上存在差异性,所以将含水介质概化为非均质各向异性。 (2)地下水系统输入、输出随时间变化,故为非稳定流。 (3)矿区水力坡度小,含水层分布广、厚度大,地下水运动符合达西定律。 综上所述,将矿区地下水系统概化为非均质各向异性三维非稳定地下水流系统。不考虑降水入渗、蒸发、河流补给以及在局部地段受断层的影响。 3.2 试验资料 对于含水层参数的选用,主要参考以往的和矿区详查的抽水试验资料。根据白涧矿区奥灰含水层的岩性、富水性、埋深条件、厚度以及水位动态对矿区进行分区(见图2),设置不同的渗透系数和储水率(见表1)。 4 利用观测点资料来反求水文参数与验证边界 模型识别验证亦称“反演”,就是数学运算中的解逆问题,即利用水头函数解算地下水均衡方程。做法上表现为根据观测点的资料,来反求水文地质参数与验证边界,直到误差满足精度为止。 本次选用2007年1月至2010年12月拟合期内奥观12#和3#水源井水位观测孔来进行模型识别与验证,采取计算和地质条件分析相结合来及时指导调参,以取得最佳拟合效果。 (1) 模型拟合曲线 奥观12#、3#水源井水位历时拟合曲线见图3、图4。 通过两井水位拟合曲线(图3和图4)可以看出,模拟的水位动态曲线与实测的水位动态曲线达到了较好的拟合,两者的动态变化过程比较吻合。 (2) 观测孔间距的影响 随机选30天和720天水头梯度场进行拟合,结果如图5所示。