基于有限差分原理宝山铁矿下告矿段地下水数值模拟研究.docVIP

基于有限差分原理宝山铁矿下告矿段地下水数值模拟研究 　　[摘要]在系统分析紫金县宝山铁矿下告矿段地质及水文地质条件的基础上，建立研究区地下水三维数值模型，利用基于有限差分原理对建立的模型进行求解。运用识别后的模型对矿山不同开采阶段地下水流场进行了预报。 　　[关键词]有限差分原理 采矿 地下水流场 　　[中图分类号] TF521[文献码] B [文章编号] 1000-405X（2015）-6-239-2 　　1研究区概况 　　研究区地形上为一四周山体环绕，中间低洼的山间小盆地，矿体即赋存于盆地内。区内地表水系不发育，仅有下告河及其支流，汇水面积约30km2。研究区多年降雨量在1600～1900mm之间。 　　2水文地质概况 　　根据区域地下水赋存条件，含水层水理性质和水力特征，矿区地下水可分为松散岩类孔隙水、岩溶裂隙水和基岩裂隙水三种类型。 　　2.1松散岩类孔隙水 　　赋存于第四系坡积层及冲洪积层中，前者分布于下告盆地山前平缓地段，后者分布于下告河冲洪积地带，平均厚度18.98m。单位涌水量1.823～4.04L/（s?m），渗透系数为12.35～40.82m/d，富水程度中等～丰富。 　　2.2岩溶裂隙水 　　赋存于本区第四系冲洪积层之下，为中上石炭系隐伏可溶岩含水层（C2+3hta、C2+3htb-1），该含水层为似层状岩溶分布岩带，地下水的赋存空隙为充填型溶洞及溶蚀裂隙。地下水单位涌水量为2.591～4.430L/（s?m），富水程度丰富。含水带渗透系数为3.45～10.04m/d。 　　2.3基岩裂隙水 　　基岩裂隙水分布面积广，可分为层状岩类裂隙水和块状岩类裂隙水二个亚类。主要分布于下告盆地四周，富水性弱。 　　2.4断裂破碎带富水性 　　该矿段内共有5组断裂通过，富水性强弱不均。 　　3地下水数值模拟 　　3.1水文地质概念模型 　　3.1.1模拟范围 　　本次模拟范围为整个下告盆地第四系含水层，另外，考虑到下告盆地四周基岩裂隙水对下告盆地第四系潜水有一定的补给，所以将模型范围下告盆地第四系含水层范围向外延20～50m，型范围面积为0.7464km2。 　　3.1.2含水层概化 　　根据下告盆地水文地质条件可知，下告盆地主要有第四系含水层、岩溶含水层及裂隙含水层。第四系主要岩性为砂质粘土、中砂、粗砂和砾砂等。地层平均厚度为20.46m，最大101.48m。其中，粉细砂、泥质、卵砾石混杂为主，等砂性土富水性辍好，赋存有丰富的地下水，为区内第四系主要含水层；下部岩性为粉质粘土及粉质粘土夹碎石，分布于盆地边缘缓坡及盆底砂卵石层下部，厚0～50.30m不等，富水性和透水性均较差，为相对隔水层。本次仅模拟下告盆地第四系潜水含水层。 　　3.1.3边界条件处理 　　模型的四周边界均概化为流量边界，模型上边界为降水补给、蒸发和井排泄边界，模型下边界为隔水边界。由于矿山对下告河进行了改道，并进行了防渗处理，地表水对地下水的补给可忽略，所以本次模拟不考虑下告河对地下水的影响。 　　3.1.4源汇项的概化 　　区内第四系潜水补给项为降水入渗补给和侧向径流补给，排泄项为侧向径流排泄、蒸发排泄及矿山生产排水。 　　3.2地下水数值模型 　　3.2.1数学模型 　　将模拟区地下水流概化成平面非均质各向同性、空间单层结构、非稳定地下水流系统，可用下列微分方程的定解问题来描述： 　　式中： 　　h―为地下水水位标高（m）； 　　K―为地下水水平渗透系数（m/d）； 　　Kz―为地下水垂向渗透系数（m/d）； 　　P―为源汇项强度（m/d），包括降水入渗补给、蒸发、井的抽水量； 　　h0―为潜水含水层的初始水位标高（m）； 　　μ―为给水度； 　　q―为侧向单宽补排量（m2/d），流入为正，流出为负，隔水边界为0； 　　Γ―为地下水第二类边界。 　　3.2.2定解条件的处理 　　初始条件：采用2010年10月矿山生产时的地下水水位监测数据，利用内插法和外推法获得潜水含水层的初始水位。 　　定水头边界：将模型范围四周边界确定为流量边界，下告盆地上下游进出水口流入边界处理为注水井，其它为零流量边界。 　　3.3模型的识别与检验 　　根据2011年9月的地下水流场对模型进行识别，重点通过各含水层参数的调整来达到流场的拟合。拟合效果见图1，从图可见两个流场的形状和绝对水位都相差很小，拟合效果较为理想。 　　地下水均衡分析： 　　通过模型识别验证，得出了模拟区地下水系统水量均衡结果（表3-1）。由表可以看出，研究区地下水系统总补给量为94.84×104m3。总排泄量为114.39×104m3。均衡差为-19.55×104m3，处于负均衡状态。 　　3.4预测结果