海州露天矿闭坑后地下水污染模拟研究.pdfVIP

海州露天矿闭坑后地下水污染模拟研究 陈建平，曹兰柱 辽宁工程技术大学资源与环境工程学院，辽宁阜新（123000 ） E-mail ：chenjianp@ 摘 要：本文针对海州露天矿闭坑后采场积水造成的地下水均衡和地下水化学平衡的破坏， 建立地质模型和地下水运移模型，应用有限元法对地下水流场和地下水污染范围进行了数值 模拟，可对地下水的污染范围进行预测和预报，基本反映了污染物运移的时空变化规律，减 少了圈定污染范围的工作量。 关键词：闭坑；地下水；数值模拟 0 引言 阜新矿区，在采矿过程中由于矿坑疏水，致使矿区周围地下水位下降，形成一个大的漏 斗区域，并且导致一系列的环境问题。海州露天矿破产闭坑后，在城市中形成有一面积达 2 3 6km ，深度达 350m 的大坑以及采矿形成的堆积量 1 亿 m 以上的大型矸石山 3 座，堆积量 0.1 亿 m3 以上的中型矸石山4 座，堆积量 0.01 亿 m3 的小型矸石山 5 座，不具规模的矸石山 340 多座。当露天矿闭坑后，采场内将充水形成人工湖，水的主要来源有地下水的涌入、采 场周围集水区雨水的汇集，其中雨水进入采场时通过地表径流和渗流进入采场。由于采场周 围大量矸石的堆积，在雨水进入采场时会有矸石的淋溶水进入采场，采场积水后尾矿的浸出 液也会进入水中，致使矸石和尾矿中的铝、砷、铅、硫化铁等可溶性盐和一些有机污染物进 入到采场积水中。随着采场内水位的上升，采矿期间形成的地下水漏斗将逐渐消失，地下水 均衡会重新破坏，当采场内水的补给与地下水处于动态平衡时，采场人工湖水与地下水之间 就会有渗流或径流相互补给，采场内受污染的水就会进入地下而污染地下水，破坏了原来的 水化学平衡[1] 。本文研究露天矿闭坑后矿坑不同积深度对地下水污染情况。 1 污染物扩散模型的建立 1.1 污染物迁移转化基本理论 污染物在包气带土层及其在地下水中的迁移转化规律是十分复杂的。其中，控制污染物 质在包气带和含水层多孔介质中迁移转化的因素和过程有：包气带土层及含水层孔隙的结构 和分布等特性参数，边界和初始条件；污染物质的类型，污染源的几何形状及污染物的释放 方式（对流、水力弥散）；地球化学和生物化学反应；以及放射性衰变、生物降解等[2] 。污 染物在包气带土层及其在地下水中的运移速度和浓度的时空分布，在较多的情况下是上述各 因素和过程综合作用的结果。 1.2 地下水在土岩流介质中流动模型 在考虑地下水中污染物及水的流动时，土岩孔隙空间的某些孔洞应相互连通，但水流通 过多孔相互连通的孔隙时，相互连通孔隙的某些部分对流体的流动亦是无效的。如死端孔隙 或只有一条狭窄的小缝和相互连通的孔隙间联系着的盲孔隙。在这样的死端孔隙中，流体几 乎不发生流动。建立基于连续介质的数学模型，假设满足三个条件：假设水流与真实水流的 推动力相等；通过同一断面的流量相等；在孔隙中所受阻力相等。满足上述条件的假想水流 称为渗流[2] 。这种水流可以看作是连续水流，可以应用流体质点的连续方程和运动方程来描 述水在土岩介质中的流动。这样，宏观上用多孔介质特性的最小体积描述地下水运动相关要 - 1 - 素，将这个最小体积定义为模型单元，含水层介质就可以认为是一种假想的连续介质。同时 在模型时要考虑污染物的对流、弥散、与介质的和离子交换作用。 1.3 污染物运移方程 地下水在含水层中运移的控制方程为： ⎡ ∂ ⎛ ∂h ⎞ ∂ ⎛ ∂h ⎞⎤ ∂h ⎢ ⎜Kxx ⎟+ ⎜K yy ⎟⎥ µs (1) ⎣∂x ⎝ ∂x ⎠ ∂y ⎝ ∂y ⎠⎦ ∂t 式中，Kxx ，Kyy