2021-4-9 | 矿业研究论文
岩体渗透张量的获得主要有以下两种方法:(1)现场专门压水试验(Louis三段压水试验，Hsich＆Neuman交叉孔压水试验)，虽然可较好地反映试验区段的地质环境，但代表范围有限而且试验费用较高而难以在实际工程中大量使用。(2)根据裂隙几何参数计算渗透张量，不仅在理论上是完备的，而且在实际应用中按照一定的抽样原则大量测量裂隙几何参数，可得大量的岩体渗透张量数据，是本文讨论的渗透张量的获取方式［2］。迁安铁矿区属温带大陆性气候，冬春两季干旱少雨，蒸发量大于降水量，多年平均降水量648．48mm，雨季6～9月份四个月占全年降水总量的82．1%，多年平均蒸发量1629．4mm(1976～2010年)，常形成旱季。矿区位于燕山山脉南麓，属于低山丘陵地貌。地势总体西北高、东南低，自西北向东南倾斜，并由低山丘陵逐渐向冲洪积平原过渡，海拔70～446m，地形起伏较大。矿区广泛出露的地层有太古界迁西群古老的变质岩系，西部高山地带的中元古界长城系地层，二者呈角度不整合接触，北部及中部有侏罗系后城组砾岩分布，第四系位于平地(平原)及沟谷中。根据钻孔资料和研究矿区水文地质结构特征以及风化程度，垂向上有第四系和太古界迁西群三屯营组(Ars)，按渗透性分为覆盖层，强风化岩体，微风化及新鲜岩体，其中第一层覆盖层据钻孔资料层厚1～24m，第二层强风化岩体层厚11～25m，第二层以下的为微新岩体。矿区南部发育近东西展布的逆断层F1，在上盘靠近下盘的地方出露断层泉，流量为0．014L/s，即断层为阻水断层。矿区岩体结构分层剖面见图1（图略）。本次研究共计在矿区测量裂隙229条，裂隙较为密集，间距均较小，平均0．5～1mm。通过对裂隙参数进行统计分析，矿区主要发育6组裂隙(见图2，图略)。通过图1的统计，结合现场裂隙统计数据，得到控制渗流的岩体裂隙的参数(见表1，表略)。
矿区岩体各向异性渗透参数获取
主渗透张量的方向只取决于裂隙组的产状。据ЧЕРНЬ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中，得到较与实际情况相符的各向异性的地下水渗流模型，能够较准确的描述地下采矿工程在生产过程中其周边地下水渗流场的变化趋势，进而较准确的预测出敏感点的地下水水位变化幅度，为工程预防和地下水环境的保护提供较为可靠的依据。
本文做：王清 张强 唐海平 谢怀前 单位：四川地质灾害防治与地质环境保护国家重点实验室 中国电力工程顾问集团西南电力设计院