数值模拟技术动态控制帷幕注浆.docVIP

数值模拟技术动态控制帷幕注浆 容玲聪 （长沙矿山研究院 长沙市 410012） 摘 要：本文介绍利用数值模拟技术优化矿山帷幕线位置、帷幕参数、帷幕注浆孔位置、动态指导施工，达到减少注浆钻探工作量和工程投资，提高帷幕堵水效果，最大限度的控制矿山水害的目的。 关键词：矿山帷幕注浆、数值模拟、动态控制水害、专家鉴定、国际领先水平 数值模拟动态控制帷幕注浆的应用概况 新桥矿业有限公司是一个坑露联合开采的大型多金属硫化矿床，矿区被分为东西两翼，西翼采用地下开采，东翼目前是露天开采，往后则准备转入井下开采。矿区东翼水文地质条件复杂，东翼矿体的顶板为栖霞灰岩含水层，灰岩岩溶发育，该层渗透系数一般为2.4m／d，矿坑涌水量预计达50000m3/d以上，属岩溶大水矿床。在基建和开采过程中，多次发生大的突水事故，突水的同时伴随着地表发生地面塌陷和开裂，致使河流断流、河水倒灌矿坑，矿区公路、铁路、农舍、地质环境被毁。矿农矛盾日益加剧，矿山安全生产受到严重威胁。为解决东翼地下水害问题，矿区东翼采用以矿区注浆帷幕截流为主的防治水方案，并开展了帷幕注浆的可行性试验研究。由于工程规模大，投资高，为提高注浆堵水效果，在设计和施工过程中，我们采用数值模拟技术优化帷幕线和注浆孔位置、动态指导和控制施工，以最大限度的降低工程造价，提高堵水效果。目前的堵水率已高达77.96%，远超过65%的原设计标准。 设计时数值模拟优化帷幕线位置 在岩溶大水矿区建造帷幕，必须具有如下条件： 较集中的过水通道，相对较短的过水断面； 可靠的隔水边界； 可靠的隔水底板。 以前，一般根据勘探孔资料确定帷幕线位置，再在帷幕线上均匀布置注浆孔。但由于岩溶发育的不均匀性，勘探孔的数量又有限，因而，选择的帷幕线和注浆孔布置往往不是投资最少、堵水率最高的最佳位置，许多注浆孔起不到预期的作用，达不到预期的效果。从而造成了投资的浪费。为最大限度的节省投资，达到堵水率最高的目的，我们在设计帷幕线位置时，建立了水文地质模型，通过数值模拟确定帷幕线最佳位置。 2.1 数值模拟优化帷幕线位置的理论依据 2.1.1 建立东翼水文地质模型 （1）边界条件概化 根据新桥矿区地质、水文地质条件，东翼水文地质模型的边界条件概化如下：东南部的高骊山组（C1g）及五通组（D3w）砂页岩为隔水边界， 19线闪长岩株为西南侧隔水边界，北北西方向的大隆组（P2d）砂页岩为北北西方向的隔水边界；东北方向由F1断裂组成弱透水边界，北西及北北东方向由茅口灰岩(P1m)为无限补给边界。区内的新西河、原圣冲河及岩溶塌陷构成东翼的局部越流补给边界。 2.1.2含水体的概化 （1）含水体的空间形态分析 由上述边界围成的区域内，含水岩体大部分被第四系覆盖，露天坑局部裸露，含水体的顶面以第四系底面为界（在5米左右），变化不大。矿体一带较高，取值时按剖面图实际情况取值。区内有Fc张扭性断层，它已被注浆试验孔揭露，岩溶裂隙率高达50.43％，是东翼地下水主要迳流通道，也是帷幕的主要封堵地段。含水岩体底板标高变化很大，在50～-150之间,已揭露的按实际含水层顶底板标高取值，未揭露区域的高低变化与剖面图基本一致。 （2）含水岩体富水性分区 由新桥矿历年的勘探资料可知，矿区灰岩地区渗透系数具有明显的各向异性。 根据D2孔抽水试验及措施井排水情况，并结合各时期的勘探资料、92年帷幕试验、2001年优化试验实际揭露的地质情况、渗透系数、帷幕线位置等综合考虑，将整个计算区域划分为21个区域（见模拟分区示意图）。并确定正南方向为主渗流方向。 图1： 模拟分区图 2.1.3模拟区地下水的排泄 模拟计算时，东翼主要的出水点有四处，一是露天坑13线左右，-86m标高，平均日出水量4800m3/d左右，二是露天坑11线，-85m标高出水120m3/d左右，三是1号供水井的抽水量大约为2286m3/d，3号供水井的抽水量大约703m3/d。 2.1.4模拟区地下水特征 根据几十年的长期观测，在矿山开采期间，由于有部分地下水储量的释放，矿坑涌水量和地下水都处于非稳定状态。加之，随着开采深度的不断加深，降落漏斗范围不断扩大。地下水渗流场应属于三维流动。但考虑到东翼地下水主要沿着Fc运动，而且绝大部分都在露天坑内排出。现有资料也有限，所以近似作为二维处理。 2.1.5 数值模拟 取主渗流方向与X轴方向平行，在上述计算区域内地下水流偏微分方程为： ① 上式中X、Y代表平面坐标值，且(X，Y)∈Ω，Ω为计算渗流区范围，H为含水层水位高度，KX，KY为含水层的主、次渗透系数。Kzi为越流含水层渗透系数。Hzi为越流含水层的水头，Mzi为弱透水层垂直渗透距离。Qi为第i个东翼涌水点(包括供水井)或地表塌陷下灌点的流量(涌水为负，下灌为正)。 在东、东南、南、西南、西北隔水边界