成兰铁路老采空区铁路路基稳定性分析.docxVIP

成兰铁路老采空区铁路路基稳定性分析 一、 注浆后地表沉降变形预测 近年来，随着中国经济的发展，交通基础设施建设越来越多。许多地方不可避免地修建了铁路和公路，尤其是在中国东部、华北、东北地区和西南部的老矿区和新矿区的区域。因此采空区问题作为铁路建设的不良地质问题已越来越突出, 如果不对地下采空区进行加固处理就可能引起采空区的地表沉陷并导致路基倾斜, 影响铁路工程建设质量和经济效益。 采空区属于隐蔽工程, 运用一般的连续介质力学难以计算注浆后地表的沉陷变形特征, 因此数值计算成为一种简单易行的注浆后地表沉陷预测方法。闵卫鲸等运用Flac2d对胶济线采空区进行了注浆处理的数值模拟与评价, 验证了路基加宽宽度对地表沉降影响;赵树森等对福建省漳龙高速公路下伏煤矿采空区采用水泥浆、水泥沙浆和水泥玻璃浆液进行充填处理, 较好的控制路基沉降变形;王生俊等运用Flac3d预测某高速公路下伏采空区地表剩余沉降量;谢和平等人运用Flac计算方法在煤矿开采沉陷预测中取得了较好的结果, 弥补一般经典方法不能考虑断层影响的不足。本文以拟建成都-九寨沟线雎水镇韩家巷子矿区老采空区为研究对象, 结合现场工程地质调查, 运用Flac3d对DK73+050~DK73+200路基下覆采空区注浆前拉应力变化特征、注浆前后地表沉陷变形规律进行研究。 二、 采空区的特点 雎水韩家巷子煤矿于1976年建井1978年投产, 开始开采范围为韩家巷子煤矿以东、标高+25m以下煤层, 矿区范围内自上向下主要可采煤层为1#、2#、3#、4#、5#煤层, 均以18o±北东向倾斜。各煤层开采方式一次采全高的多为长臂式开采, 回采率可达80~90%, 少量采用房柱式或刀柱式开采, 回采率可达70~80%。由于开采时间和位置不同, 顶板管理方式亦不同, 根据收集资料、物探、钻探等成果判定, 既有公路两侧各煤层均已大量开采。煤系地层中煤层采空后形成采空区, 其上覆岩层工程地质稳定性变差, 易发生冒落塌陷, 该采空区具有煤层采空区缓倾、多层采空区、重复采动、采空区埋深较浅等特点。 矿区位于雎水煤田的东段, 汉旺向斜的东南翼, 主要地层由老至新有中生界三叠系石灰岩、泥灰岩、泥岩;古生界二叠系茅口组含煤地层;中生界上侏罗系红砂岩及第四系覆土。总体构造形态复杂, 地层产状NW290~310°/NE∠18°, 在龙门山前山断裂的基础上发育有F1及F2两条断层, F1和F2断层由数条小断层组成, 构成密集的断层组, 且具逆断层性质。 三、 数值模型的构建和过程 1. 结构效应的相关研究 Flac3d数值计算结果的可靠程度与岩体力学参数、原岩应力、岩体结构类型、计算域的结构效应及其开挖、加载等约束条件有关。根据采空区工程地质条件及采空区的实际分布情况, 结合研究区域的岩层柱状图、采空区岩体结构面分布状况和岩土体物理力学参数 (如表1) , 综合确定计算模型的基本参数。 2. 破坏函数模型 张鲁渝、郑颖人等人研究表明, 模型边界应该在基准方案基础上扩大1~3倍, 模型边界对模拟结果的影响可以降低到工程可接受的范围内。根据上述原则, 采空区模型边界以铁路线路为中心, X轴与线路走向重合, 长度500m;Y轴垂直于线路走向, 宽度500m;Z轴为地层方向, 考虑模型边界影响, 距地面400m, 模型尺寸为:500×500×400m3(如图1) 。 该模型的力学边界采用两边侧面 (X方向) 、前后缘侧面 (Y方向) 水平约束, 底面 (Z方向) 垂直方向约束, 地表为自由面。整个三维模型共划分11, 332个节点和59, 185个单元, 其中采空部分采用零单元 (null) , 注浆体采用壳体单元 (shell) , 本构模型采用弹塑性本构模型, 破坏准则采用摩尔库仑屈服准则。摩尔-库伦屈服函数破坏包络线可表示为: 式中:σ1、σ3为材料破坏时最大、最小主应力;c为岩土材料的内聚力;Nφ= (1+sinφ) / (1-sinφ) , ?为岩土材料的内摩擦角。 3. 采空区初始充放电 采空区不是随着地质体的形成而固有的, 而是在煤层开采过程中逐步形成的, 即:采空区地表沉陷是一个过程。本文通过FLAC3d中的时步 (step) 来控制模拟采空区地表沉陷变形过程, 每一个过程变形时间差异通过时步的不同来反应。模拟计算首先按计算深度对模型施加自重荷载, 让模型在自重应力作用下稳定, 建立采空区的初始自重应力场, 然后清除历史上自重 (初试应力) 作用造成的位移、速度场, 再按照实际开采的顺序, 即1#煤层、2#煤层、3#煤层、4#煤层和5#煤层每一循环模拟一层开采, 直到5#层开采完成, 预测现今采空区地表剩余沉降变形量和各层采空区顶、底板拉应力的发育规律;然后采用壳体单元 (shell) 来模拟粉煤灰注浆体, 仍然