各项岩体参数分析矿区地下水环境推测.docVIP

各项岩体参数分析矿区地下水环境推测 岩体渗透张量的获得主要有以下两种方法:(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．30times;106m2;建设期第三阶段，地下水降落漏斗导致刘庄子村供水井下降7m，影响面积为2．59times;106m2。其中建设阶段时间为4年，即降落漏斗扩大速度为7．31times;104m2/a。建设期完成，地下开矿进入运营期，地下水继续流入生产巷道，降落漏斗的范围将会进一步发生变化，用运营5a，10a，15a，20a的地下水降深等值线图(图6～9)来显示研究区地下水降落漏斗的变化趋势，可以预测生产对矿区周边居民集中供水点影响程度，进