基于现场实测数据的断层应力分析.docxVIP

基于现场实测数据的断层应力分析 0 断层周围地应力场 岩体扰动和不连续面的存在影响了岩体的应力变化和迁移过程。在应力适应过程中，扰动和不连续面附近的局部应力和区域应力之间存在较大差异。 从某种意义上来讲, 断层是区域地壳稳定和工程岩体失稳的主控因素之一。 在实际工程中, 为研究断层对应力场的扰动效应, 很多学者从断层周围地应力测试及数值模拟等方面进行了相关研究。 从目前的研究发展趋势来看, 主要有以下3 种途径: ①把断层抽象为被节理及裂隙分割的块体系统进行力学分析; ②把断层的力学性态作为附加条件, 采用连续介质力学的方法进行模拟; ③采用接触面单元模拟断层与岩体之间的不连续面。 在实际工程中, 断层的几何形体是易于确定的, 而其力学性态往往难以确定, 且力学性态是影响断层及其周边应力场分布的主要因素, 特定的断层力学参数与特定的应力分布状态相匹配。 本文利用断层附近及其周边测孔的测试资料, 结合正交设计理论及人工神经网络的非线性反演方法对断层的力学参数进行反演, 并对断层对应力场的扰动效应进行研究。 1 测试孔的配置和测试结果 1.1 常规水压致裂法地应力测试 某隧道位于福建东部重力梯度带上, 隧道走向为N70°W, 山脉主要走向为北北东~南南西, 场区以构造剥蚀中低山为主, 岩性以灰色、 浅灰色钙质粉砂岩为主。 为研究隧道轴线部位岩体应力分布情况及断层对应力场的扰动效应, 布置了1、 2 号2 个测孔进行常规的水压致裂法地应力测试 (见图1) , 其中2 号测孔位于断层F11和F12中间。 F11断层的产状为135°∠75°, 断层带内发育2 条断面, 使岩层产生一定错动, 带内岩石破碎, 断面呈舒缓波状, 为压性逆断层。 断层上盘影响宽度约25 m, 下盘约10 m, 节理裂隙发育, 较破碎。 F12断层产状为157°∠59°, 断层带内为硅化、 叶腊石化碎裂岩, 见断层泥及磨圆的断层角砾, 断面略呈舒缓波状, 断面上发育阶步及擦痕, 为压性逆断层。 断层上、 下盘影响宽度约20 m, 节理裂隙发育, 较破碎。 1.2 最大水平主应力测试 由图1 可知, 1 号测孔距离断层较远, 其应力分布不受断层的影响, 且测试区域内地势起伏相对较小, 其应力分布情况可视为测试区域内的应力分布; 2 号测孔断层位于断层F11和F12之间, 其应力分布受复合断层的影响较大。 由测试成果可知, 1号测孔测试区域内最大水平主应力方位位于N25°~30°W之间, 而2 号测孔位于N7°~23°E之间, 两者夹角较大。 若将1 号测孔测得的最大水平主应力方位视为区域构造应力方位, 则2 号测孔所在部位的最大水平主应力方位受断层的影响偏离区域构造方位高达约45°。 图2 给出了2 个测孔测得的最大和最小水平主应力随深度的变化情况 (图中, σH为最大水平主应力, σh为最小水平主应力, σz为自重应力) 。 由图2可知, 1 号测孔应力量值整体大于2 号测孔。 线性拟合可知, 1 号测孔最大水平主应力侧压系数 λ (λ=σH/σz) 约为1.6, 而2 号测孔约为1.3。 2 号测孔最大水平主应力方位的偏转及应力量值的弱化可归结为断层的扰动效应。 2 反演误差力学参数 2.1 最大水平地位参数确定 为方便建模, 并消除边界效应的影响, 建模时X轴取为隧道走向方位N70°W, 使断层及2 号测孔处于计算模型的中部, 计算范围取1 200 m×1 200m。 根据1 号孔测试结果, 取最大水平主应力方位为N30°W、 侧压系数为1.56, 最小水平主应力侧压系数为1.04。 由于最大水平主应力方位与计算坐标系中的X轴呈40°夹角, 需将计算模型上的最大及最小水平主应力转化为边界上的应力分量 σx、 σy和τxy, 计算公式如下 式中, β 为最大水平主应力方位与计算坐标X轴的夹角。 2.2 基于正交设计的人工神经网络非线性反演 本文采用连续介质力学的方法, 将断层及其影响带内的岩体参数进行弱化。 在一定的构造应力及自重应力组成的应力场中, 由于断层的扰动效应, 断层及其周边的应力将出现应力分异现象, 但应力分布情况与断层的力学参数相适应, 且这种应力-断层参数的响应十分复杂, 不能用简单的线性方程表示。 因此, 采用基于正交设计的人工神经网络非线性反演模型对断层参数进行反演。 分别取断层F11、 F12的弹性模量E11、 E12为0.5~4.5 GPa, F11、 F12的泊松比 μ11、 μ12分别为0.25~0.41。选择U9*(94) 表进行设计, 每个因素取9 个水平, 各水平组合见表1。 进行有限元正向计算, 将每种工况下的2 号测孔测点对应的应力量值的计算结果作为神经网络的输入, 对应的断层力学参数作为输出进行神经网络训练。 神经网络