MODFLOW在矿坑涌水量预测中的应用.docVIP

PAGE PAGE 1 MODFLOW在矿坑涌水量预测中的应用 　　摘要：MODFLOW——模块化三维有限差分地下水流动模型，在地下水模拟研究中应用比较广泛，本文针对南洺河铁矿具体水文地质条件，综合考虑各种影响因素，建立水文地质模型和数学模型，运用MODFLOW软件对矿坑涌水量进行预测。 　　关键词：三维非稳定流，数值模型，矿坑涌水量 　　中图分类号：TD742文献标识码：A文章编号： 　　1矿山概况 　　南洺河铁矿地处太行山东麓中段，地势自西向东倾斜，属邯刑地区岩溶地下水系统。矿床属接触交代型磁铁矿，矿体分布于奥陶系石灰岩含水层与下部闪长岩体之间，奥陶系上覆石炭～二叠系砂岩及薄层石灰岩。由于断裂构造发育，石炭～二叠系的砂岩及薄层石灰岩之间水力联系密切，形成一个统一含水岩组，断裂构造也是地下水自上而下运动的通道。南洺河为季节性河流，河床通过矿区，河床渗漏是矿坑涌水量的来源之一。 　　铁矿体自西向东延长2078.96m，延深达765m，平均宽度351m，赋存标高159.5～-535.6m，矿体走向NE，倾向SE。矿山针对一期开采设计了±0m、-40m和-80m三个标高水平，目前主采区位于±0m水平，设计产量50万吨/年。 　　2模型的建立 　　根据区域边界条件，界定模拟区范围为：西至南丛井、龙雾～活水断层，南至南洺河断层，东至鼓山鼓山～紫山断层，北至沿北洺河一带地下分水岭。 　　2.1概念模型 　　模拟区西部、南、东部均为隔水断层，定为隔水边界；北部与西南部为地下分水岭边界。 　　地下水在接受补给后，水平方向上自西向东～东南径流，上覆石炭～二叠系地下水向下垂向运动，地下水运动具有三维流的特征。由于实际勘探中很难获得各含水层各向异性的参数，在宏观上忽略其各向异性差异，将含水层简化为均质各向同性含水层。因此，本区水文地质概念模型按均质各向同性三维非稳定流模型处理。 　　2.2数学模型 　　在水文地质概念模型认识的基础上，建立三维非稳定流的数学模型，可表示为： 　　式中： 　　X—东西方向，y—南北方向，z—垂直方向； 　　K—含水层的渗透系数（m/d）； 　　H—为任意空间点（x,y,z）的水头随时间t函数（m）； 　　H0—为任意空间点（x,y,z）的初始水头函数（m），即t0时刻的水头空间分布状况； 　　W—源汇项，本模型中即为矿坑涌水量（m3/d）； 　　μe—弹性给水度； 　　B—模拟区边界； 　　Ω—模拟范围。 　　2.3数值模型 　　对上述数学模型，采用“VisualModflow4.2”三维渗流模拟软件来实现三维有限差分方法求解，数值模型能否真正的反映和再现实际水文地质条件，还有待进一步识别。 　　2.3.1渗流区剖分 　　选择对物理模型进行矩形剖分，在地层分界、抽水点、观测孔附近加密剖分处理，共分成113行、111列、35层，总计439005个单元格，如图01。 　　图01水文地质模型剖分图 　　2.3.2数值模型的调试与识别 　　（1）数据初值的确定 　　边界条件 　　东、西、南部为隔水边界，北部与西南部为地下分水岭，均按零流量边界处理。 　　初始流场 　　以2003年1月区域岩溶地下水流场作为模型识别的初始流场。 　　垂向补排强度 　　大气降水补给量：根据灰岩裸露区及第四系薄层覆盖区的地表岩性、地貌等因素，概化为几个降水入渗系数不同的小区，并根据降水入渗系数经验值给出初值，待模型识别时确定。 　　河流入渗补给量：参照前人研究成果、降水量和上游水库的放水情况给出入渗强度的初值，待模型识别时确定。 　　地下水人工开采量：根据实际调查资料，分别加在开采结点上。 　　水文地质参数 　　按含水地质体空间分布特征，对模拟区进行水文地质参数分区，如图02。 　　图02模拟区水文地质参数分区图 　　各分区的参数初值，参照前期研究成果给出，如表01，待模型识别时确定。 　　水文地质参数初始给定值表01 　　（2）模型的识别 　　选取矿区地下水位动态观测资料连续性较好的观测孔，用实测地下水位动曲线与模型求解的水位动态曲线进行拟合，以2003年1月至2011年4月期间作为计算时间段，时间步长为1个月，让模型运行，记录下每个时段各观测孔所在结点的水位。若各种计算初值给的合理，计算的（H—t）曲线应与实测的（H—t）曲线基本拟合，否则要反复调整水文地质参数及各种不确定因素进行试算，一直到曲线拟合程度满意为止。经反复调参计算，结果表明模型是可靠、正确的。 　　部分观测孔水位动态拟合曲线见图03、04。 　　模型识别后各分区水文地质参数如表02。 　　图03W1孔拟合曲线图 　　图04W12孔拟合曲线图 　　模型识别后水文地质分区参数值表02 　　3矿坑涌水量预测 　　3.1涌水量的基本设置 　　根据矿山设计开采条件及开采现状，需预测