- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
MODFLOW在矿坑涌水量预测中应用
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 数值模型
对上述数学模型,采用“Visual Modflow4.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孔拟合曲线图
模型
文档评论(0)