隧道开挖对地下水扰动影响与区域水运移规律分析.docxVIP

1?工程概况

本工程位于邢台市，多数地区为中低山区，地势西高东低、北高南低。区域内山峰群立，层峦叠嶂，最高山峰为天河梁，海拔1?753?m，最低点位于石店村前河滩，海拔62?m。

研究对象路罗隧道、前营隧道等隧道洞口导洞海拔在900~1?100?m，邢台大峡谷景区各景点海拔在700~900m，景区下游各村庄海拔在500~700?m。隧道位于太行山主脉东侧，受构造强烈切割影响，峡谷走向基本以西北东南向为主，山体地层复杂，区内地下水主要以第四系孔隙潜水、基岩裂隙水为主。

由于岩层节理裂隙较发育且线路高程相对较低，隧址区段隧道洞身从王山铺村以东约800?m的峡谷地下穿过，王山铺村村民饮水均取半山腰处流出的构造裂隙水，为隧道穿越的沟谷中流水的上游。而沟谷水流向贺坪峡风景区，风景区位于沟谷的下游。为研究隧道施工对村民取水和风景区水源的影响，结合无人机航拍地形图开展地表水径流分析，并基于地形三维建模研究隧道开挖对下游地下水的影响深度。

2?航测方案

采用B厂家无人机对王山铺、周公山景区入口、邢台大峡谷景区、邢台大峡谷以及前营隧道平导口5处进行重点区域精飞，采用GPS坐标系CGS2000_ALBERS，航向重叠率为85?%，旁向重叠率为80?%，相机倾斜角42°，飞行速度3?m/s，航拍过程中会采集区域内的图像资料以及地形点云信息。在此基础上，采用A厂家无人机对上游山西省一部分区域以及大峡谷下游到路罗站入口的区域开展全局点云信息采集，覆盖隧址区上下游水系，测绘范围覆盖面积为52.24?km2。

结合A厂家和B厂家的无人机地形图测绘资料，将航飞工作区域内的航测数据整合拼接出带有三维坐标信息的地形图，再将无人机采集的点云地形信息导入Rhinoceros7中，生成三维地形，发现隧址区内业有3条主干河道，分别位于前营隧道平道口处、邢台大峡谷流域以及周公山下游。路罗隧道正好横穿了三大主干河道所在的流域，因此，本研究重点针对隧址区内三大河道所在流域展开地下水运移模拟。

3?研究区汇水模型与分析

对研究区域进行汇水分析。汇水方向的确定通常采用D8（DeterministicEighthours）算法进行计算，假定被处理栅格内的水只流向其周围各栅格中的一个，并指向坡降最大的方向。将目标栅格同相邻各栅格进行比较，目标栅格中心与相邻栅格中坡降最大栅格中心的连线方向为目标栅格的水流方向。用128、1、2、4、8、16、32、64这8个特征码表示东北、东、东南、南、西南、西、西北和北8个有效的水流方向。

逆向跟踪汇水方向，统计注入该栅格的周围栅格数量，作为该栅格的汇流特征值。汇流特征值数值越大，表明注入的水流就越多。规定集水区为集水能力强栅格所在的区域，若某区域集水能力特征值为零，则代表该区域地势较高，据此可划分出研究区域分水岭的空间分布和集水区的边界。将栅格集水量为0的区域作为集水区分界线，进而划分区流域范围。图1将流域分为5个级别，但可以清晰发现，根据流域的划分同样可以将隧址区划分为3个主要的流域，划分结果和河网分布的分析结果一致。利用GIS软件将数据库中数据写入MODFLOW格式文件中。其中降雨入渗输入Recharge程序包，对应数据存入“*.rch”文件中；边界侧向补给量、开采量输入到Well程序包，对应数据存入“*.wel”文件。

图1?河网流域分布地形

为探究各水域水量汇集分布占比，需要对河网各级河道的倾泻点进行识别，具体原理即从每条河流的最高点出发，沿着坡降线往下，找到与其他支流汇集的交界处，对交界处进行标注。

图2（a）识别出了各级别河道的倾泻点，基于倾泻点，可计算出隧址区内所有河道汇集到各个倾泻点的流量大小，进而计算出地表水流量的分布及各流域的流量占比。如图2（b）所示，根据倾泻点流量统计，可将图1的3大主要流域细化分为20个级别的小流域，根据图中信息可知，图2（a）中颜色较深的流量大，流量单元最高达到了6?964，颜色靠较浅的流量越小，流量单元最低为1。

（a）

（b）

图2?流量分布地形及流量比例分布地形

（a）示意一;(b)示意二

利用盆域分析计算可知，在三大流域中，位于前营隧道平导口附近的流域流量最大，邢台大峡谷流域流量次之，最后是周公山流域一带。根据图2（b）中小球所示比例，计算得出前营隧道平导口流域集水占比为49.99?%，邢台大峡谷流域占比为33.31?%，周公山流域占比为16.70?%。考虑隧址区面积占比，前营隧道平导口流域面积占隧址区24.93?%，邢台大峡谷流域面积占隧址区53.36?%，周公山流域面积占比为21.71?%。

计算得到前营隧道平导口的集水比为2.005，邢台大峡谷流域的集水比为0.624，周公山流域的集水比为0.769。因此，考虑地表水的汇集能力，前营隧道平导口流域大于周公山流域大于邢台大峡谷流域，具体见