矿区岩体各向异性渗透参数获取
主渗透张量的方向只取决于裂隙组的产状。据ЧЕРНЬIЁВ(1979)[5]研究认为,风化作用、开挖卸荷等表生作用并不能对裂隙网络产生根本性的改变,而仅仅是有限范围内的改造,其中较明显的是裂隙面的开度,对结构面的产状、间距和长度的影响不甚明显。因此,对于某一工程地质单元的渗透张量的主方向完全可以通过裂隙样本法来确定,在此基础上可充分利用抽水试验取得的综合渗透系数数值,将两者有效结合起来,为了利用抽水试验取得的平均渗透系数数值来修正裂隙样本法取得的渗透张量,定义修正系数m为m=kmk0(3)根据(1)式和(3)式,得到反映原位地质环境裂隙岩体的修正渗透张量Km为Km=mK0=mkx000mky000mkz(4)依据前面所述方法,对研究区进行裂隙抽样调查,共计调查裂隙229条,通过对优势裂隙进行统计,即裂隙宽度,裂隙间距,倾角和倾向的统计,进行渗透张量的计算,得到Kx=0.010m/d,Ky=0.012m/d,Kz=0.017m/d。根据(2)式,求得综合渗透系数为K0=0.013m/d。对研究区进行钻孔抽水试验资料的数据进行整理,得到综合渗透系数,强风化岩体0.04m/d,微新岩体为0.033m/d。然后根据(3)式和(4)式进行渗透张量的修正。修正结果见表2(表略)。
矿区地下水渗流模拟预测评价
将修正后的渗透主值应用到VisualModflow中,进行各地层的渗透系数赋值,根据地下开采空间结构(见图3),建立各向异性的渗流模型。将采矿区基岩上覆盖层视为均质各向同性地层,渗透系数由双环试坑注水试验得到。模型中各层渗透系数赋值见表3(表略)。根据开采方案,确立地下采矿工程的空间结构,以非稳定流进行运行,预测建设期和运营期各时段的地下水降落漏斗。建设期分为三个时段,地下开采层位从浅至深,以建设期的__阶段和第三阶段的地下水位降深等值线图来表示降落漏斗影响范围的变化趋势,如下图4和图5(图略)。建设期__阶段开挖到海拔-115m,第三阶段开挖到海拔-340m。对比图4和图5可以看出,随着建设向下发展,研究区周边越来越多的地下水汇入开采巷道,地下水降落漏斗的影响范围逐渐扩大,并以近似椭圆形态出现。建设期__阶段,地下水位下降影响到了二郎庙村供水井,降幅为8m,影响面积为2.30×106m2;建设期第三阶段,地下水降落漏斗导致刘庄子村供水井下降7m,影响面积为2.59×106m2。其中建设阶段时间为4年,即降落漏斗扩大速度为7.31×104m2/a。建设期完成,地下开矿进入运营期,地下水继续流入生产巷道,降落漏斗的范围将会进一步发生变化,用运营5a,10a,15a,20a的地下水降深等值线图(图6~9)来显示研究区地下水降落漏斗的变化趋势,可以预测生产对矿区周边居民集中供水点影响程度,进而对可能出现的诸如造成饮水问题的情况进行预防,达到生产生活顺利进行的效果。以上各图进一步表明,运营期间,地下水降落漏斗继续扩大:东至孟官营村,西到蔡园村和张家窑村,北至大杨庄村,南到新庄村,共影响11个敏感点;影响面积从运营5a的5.60×106m2扩大到运营20a的9.00×106m2,扩大速度为1.70×105m2/a;影响幅度较大为二郎庙村,李庄子村和刘庄子村,地下水位降幅为30~75m。并且降落漏斗在x轴,y轴,z轴上不均匀扩大,其影响范围形状继续以椭圆形态出现。因此,地下水在岩体中渗流主要受控于岩体的裂隙,与岩体中裂隙发育形态和方向有关。矿区地下水渗流为在非均质各向异性介质中的非稳定流。
结语
裂隙测量的渗透张量方法,应用统计学方法全面对裂隙岩体做抽样调查,测量裂隙系统的水力参数,能较快获得大量计算渗透张量的基础数据,获得不同位置空间各主渗透方向的渗透系数,构成一个三维渗透系数张量场,准确地刻画裂隙岩体的渗透特性。而本文将裂隙抽样方法得到的渗透张量与抽水试验相结合,进行各主渗透方向的渗透系数的修正,得到反映原位地质环境裂隙岩体的修正渗透张量,从而得到更与实际相符合的各向异性地下水渗流场。将修正了的渗透张量应用于数值模型VisualModflow中,得到较与实际情况相符的各向异性的地下水渗流模型,能够较准确的描述地下采矿工程在生产过程中其周边地下水渗流场的变化趋势,进而较准确的预测出敏感点的地下水水位变化幅度,为工程预防和地下水环境的保护提供较为可靠的依据。
本文做:王清 张强 唐海平 谢怀前 单位:四川地质灾害防治与地质环境保护国家重点实验室 中国电力工程顾问集团西南电力设计院