高速公路下伏采空区剩余沉降量FLAC_3D_计算方法
采空区剩余沉降量FLAC_3D_数值模拟计算方法
陈朝清
2
李森林
3
2、 福建永福工程顾问有限公司, 福建 福州 350003
(2) 式中: σ1S 是在岩石峰值强度时的最大主应力; σ3 是最小主应力; m、 s 是材料常数, 取决于岩石性 质和原始破裂状况; σc 是岩石单轴抗压强度。 此外, 当拉应力超过材料的抗拉强度时, 材料将发生拉裂 破坏。表 2 给出了各岩层和煤层的物理力学参数。 表 2 各岩土层和煤层物理力学参数
强度的水泥土桩体,将严重影响粉喷桩的强度, 在 这种情况下必须增加粉体的掺入量和采用复搅的 施工工艺。 高含水量、 大孔隙比和粘粒含量多时, 土 周边的束缚力极低, 当钻头反转提升喷灰时, 产生 一个垂直向下挤压力和一个径向水平推力, 由于土 呈流塑状, 束缚力极低, 桩体在成形过程中向下及 向四周水平向排水, 影响形成竖向桩体, 通常形成 所谓 “掉桩” 或 “下沉” , 常为地表下 1~2m。当发生 “掉桩” 或 “下沉” 时, 只要当时采取立即回填土并重 新复喷复搅, 就能克服这种现象, 如果不作处理将 会造成过大的路基沉降。 3 粉喷桩复合地基承载力和粉喷桩单桩承载 力的关系 粉喷桩复合地基的平均允许承载力公式为: =a× [σ 桩] +(1-a) [σ 土] [σ 复合] 式中: σ 复合—— —复合地基的平均允许承载力; a—— —置换率; σ 桩—— —搅拌桩的允许承力力; σ 土—— —天然地基土的允许承载力。 无疑,桩的强度将直接影响复合地基的强度, 假设桩的强度不断增加而土的强度依然不变, 按照 公式的计算, 复合地基的强度也会不断增加, 然而 因为粉喷桩从本质上来讲是属 实际情况并不如此。 摩擦桩类, 当土的强度不变, 而且饱和软粘土的强 度很低时, 在这种情况下, 即使不断增加桩的强度, 这如同一根 但总的复合地基强度也不会随之增加。 筷子在浆糊里和一个钢筋在浆糊里的情况是类同 整个复合地基 的。只有当天然土的强度也增强时, 的强度才会增加。所以桩的强度应适度, 要和天然 这当然 土的强度相互匹配。反过来桩的强度很低, 也是不行的。 4 桩土置换率及粉体掺入量对复合地基强度 的影响 在实际设计运算时往往提供所要求达到的复 拟定的粉喷桩强度及天然地基土的强 合地基强度、 度来求桩土的置换率, 计算出桩数, 然后布置桩位,
浅谈高速铁路路基工后沉降及其计算方法
中圈 分 类号 :U 2 3 8 文 献标 识 号 :A 文 章编 号 :2 3 0 6 - 1 4 9 9( 2 0 1 3 )0 9 - 0 0 7 4 - 2
1 . 路基工后沉降的定义和组成 路基 在填筑过 程 中 ( 至铺 轨前 )所产 生的沉 降称为施 工沉 降,这 部 分沉 降可 以采 用填补加 高来解 决。路基在铺 轨完成 后所产 生的沉 降称 为 工后沉 降 ,这部 分沉 降只能 以抬 道补碴来 调整 ,它 直接影 响到线路养 护 维修工作量和高速铁路 的运营 能力 。路基工后沉降 由路基填土压密下沉、 行车引起 的基床累积变形和软 土地基产生 的工后沉 降三 部分组成 。 1 . 1路基填土压密下沉 路基 填土压密 下沉,是 由填土 的 自重引起 的,它发 生在两个 阶段: 是施 工阶段 的下沉 ,不计入工 后沉 降;二 是施工 完成后对 后期运营 有 影 响的工后沉 降。 由散体 材料填筑 而成 的路基本体 产生一 定的压密 下沉 是 正常的 ,其 大小取 决于填料和 施工质 量。如果 下沉量较 大,说 明填 土 的压实度不足 、强度低 ,容 易造成 不均匀变 形。 目前世界各 国关于路堤 填土 的压 密下沉通 常是通过压 实密度 予以保 证 的。例如其 中较 具 代表 性的 日本 对填土 的压实质量 采用值 作指标 ,为 了保证填土 具有
填料 缀配砂 砾石
釜 - 、 l
一 一 、 一
实压标 准
孔隙 率Ⅱ
l 地基蔡 数瑚 《 鲴 P | 触) l
0 : 7 0
I
>1  ̄9 0
I
( 1 8 %
表 卜2 基床底层填料及压实标准
地基系数 K3 O 《 a 舨) ≥1 1 0 ≥l 3 O ≥l 5 o
高速公路下伏多层采空区地表沉陷的预计评价方法
( n t t f e tc nc l n ie r g o t a tU e y a j g2 0 9 Isi e o o h l gn e i ,S uh s mv  ̄t ,N n n 1 0 6,C ia t u G e aE n e i hn )
Ab t ac : sr t
远 超 过 允 许 值 , 采 空 区必 须 进 行 处 治 . 该 关键 词 :采 空 区;地 表 沉 陷 ;预 计 评 价 ;概 率 积 分 法 中图 分 类 号 :T 4 3 U 3 文 献 标 识 码 :A 文 章 பைடு நூலகம் 号 :10 —0 0 ( 0 2 0 .6 80 0 1 5 5 2 0 )40 4 .5
to u e .Usn he i tg a d lo ls e tng l ra lc n a tt n,h e o a e a a e o rd c d i g t n e rl mo e fc o e r c a ua nd b o kig p rii o t e m t d c n b d ptd t h h h nc fmi i g d pt te c a e o n n e h , m n n ik e s, c a e ic nd s ap ,L l a o l x s v rl srk i igt c n s h o lb d p th a h r £ wel s c mp e e e a ti e S g a , a d t e e oe t e p g mn v ra ii d t e p e c in p e iin C l b m p v . Gru d s sd n e o f n r fr m r e s tlt a h r dito r c so al e i r e h h yn o d on u i e c b o n ln d c a e a a iy b r ditd b o r i ae c n e so d o v n in as d dih a fi c i e o lb d C n e l s e p e c e y c o d n ts o v r in a c n e to a r ie a p e d n l n me o r d fe h t d a e mo i d. i K e o d g a ; g o d s sd n e; p dit v u t ; ob bii n e r to t o y w r s: o l r u u ie c n b e r c e e a a e pr a lt i t g ai n me d d l y h
基于FLAC3D数值模拟的老采空区剩余沉降量分析计算-以山东枣庄安博
区剩余 沉 降量 是采 空 区顶板 覆岩 塌 陷形成 破碎 带后
1 工 程 区 概 况
1 . 1 区域 采 矿 条 件
继续 下沉 压实 所 引起 的地 表 沉 降值 , 其 中也 包 括 破
覆 岩 已塌 陷形成 冒落 带 、 裂 隙带 和弯 曲带 , 在 经过 一
段 时 间 的 自然 沉 压 后 采 空 区 逐 渐 趋 于 稳 定 I 1 ] , 而 使
在误 差 。近 年来 , F L AC数 值 模 拟 在 岩 土 地 下 工 程
设 计研 究 中应用 越 来 越 广泛 , 其 基 于代 数 表 达 式 代 替 空 间离散 点处 不 确 定 的场 变 量 , 利用 有 限 的差 分 法 求解 代数 方 程 求 解 。F L AC 。 。 适 合 非 线 性 材 料 在 外 力下 的屈 服流 动 等 大 变形 模 拟 , 而越 来 越 多 的应 用 研究 到 下伏采 空 区稳定 性模 拟 l _ 6 书 ] , 采空 区注 浆加 固模拟 等 _ 9 。下 面 以山东 枣庄 市 安博 化 工 厂 区建 设 项 目为 例 , 对 老采 空 区剩余 沉 降模拟 分析计 算 。
t i o n a n d An a l y s i s o n R e s i d u a l S u b s i d e n c e A mo u n t i n 0l d —G o a l B a s e d o n F I AC N u me r i c a l S i mu l a t i o n E J ] . S h a n —
高速公路下伏采空区稳定性评价方法探讨
公路 中线 斜 交 , 即有一段 路基 处在 了巷道 的正上 也
方 , 处巷 道 埋深 约 10 巷 道 的宽度 约 为 1 m。 此 0 m, 0 自地 表 向下 , 矿采 空 区大体 可分 为 如下 岩性 组 : 煤 地 表 为粘 土 、 粉土 、 卵石 层 , 度 约为 7 m, 砂 厚 0 地表 州市 , 连接郑 州
对路 线有 一定 的影 响 。
给 公路 路基 、 路面 的稳定 带来 极大 隐患 。因此评价 采 空 区地 基 的稳定 性 ,采用 合理 的工程 技术 措施 ,
是 实 际工作 中必须 解 决 的具体 问题 。 由于地 质条 件 复杂性 、 开采方 式 多样 性 以及 勘察过 程 中搜 集 资
宝旗煤矿工程地 质特征
该矿 于 19 95年建 井投 产 , 04年 闭井停 产 。 20
该区煤柱分层 分布, 煤线 的走 向大 致呈南北 向展布 ,
料局 限性 , 以对 采空区 的稳定性进行准确 的评价 , 难
高速 公 路 穿越煤 矿 采 空区 的 问题也 成 了近 年来 公
埋 深在 7m以下 , 0 煤层 的倾 角 陡直 , 大约在 8 0。左 右。 该矿 井 有一条 开采 巷 道 , 巷道 的 宽度和 高度 为 3 巷 道埋 深约 10 m, 0 m,走 向近 南北 ,向北 延伸 至
【 关键 词】 高速公路 采空区 稳定性 评价
西南绕 城 高速公 路 , 途径 郑 州市 、 昌市 和平顶 山 许
l 引言
近年来 , 随着 高速公 路建设 范 围不 断扩大 以及 人 为坑 洞不 断增加 , 路线经 过采 空区 的情 况也会 经
常遇 到 。高速 公路 下伏采 空 区的存在 , 使上 覆岩层 冒落塌陷 , 地表将产生不 连续沉 降或连续沉 降盆 地,
高速公路路基沉降计算方法
高速公路路基沉降计算方法“高速公路路基沉降计算方法”是一种用于计算高速公路路基可能出现的沉降情况的方法。
该方法首先对道路路面上的沉降和抬升异常进行检测,然后通过分析历史数据、地形因素、材料特性以及施工工艺等因素,进行沉降风险分析和预测,最后给出防止沉降发生的具体措施和建议,以便采取有效的技术措施来减少或避免沉降的发生及其后果。
高速公路路基沉降计算方法主要包括以下几个步骤:1、观测路面沉降和抬升异常首先,根据实际情况选择合适的监测技术,定期观测路面沉降和抬升异常,并以图表形式记录和统计数据。
2、分析历史数据其次,根据历史监测数据,对路基沉降现象进行分析,如果发现路基沉降的规律性,应当加以关注,并对历史数据进行详细的分析,以便进一步了解路基沉降的原因。
3、考虑地形因素在考虑路基沉降计算时,还要考虑地形因素,如:路基处于低洼地带,或者路基处于水库、河流、湖泊等易受水位影响的地带,都会对路基沉降产生影响。
4、分析材料特性路基沉降计算还要考虑材料特性,如:路基所用的基础材料的性质,以及混凝土配合比等,都会影响路基的沉降程度。
5、分析施工工艺此外,路基沉降计算还要考虑施工工艺,如:施工过程中使用的机械设备,是否采用了可以有效阻止沉降的技术措施,以及压实温度是否适宜等,都会影响路基的沉降情况。
6、给出技术措施和建议最后,根据以上数据和分析,给出可以预防路基沉降发生的具体措施和建议,以保证路基的安全及质量。
总之,高速公路路基沉降计算方法是一种用于计算高速公路路基可能出现的沉降情况的方法,包括观测路面沉降和抬升异常、分析历史数据、考虑地形因素、分析材料特性和施工工艺,最后给出防止沉降发生的具体措施和建议,以便采取有效的技术措施来减少或避免沉降的发生及其后果。
高速路基下伏采空区治理的关键技术
西部探矿工程
13 6
高速 路 基 下伏 采 空 区治 理 的关 键 技 术
谷 水 兵
( 南省 岩土工 程 有 限公 司, 南 洛 阳 4 12 ) 河 河 7 0 3
摘
要: 随着高速公路 的建设和完善 , 采空区治理技术在 高速公路建设 中的作用越显重要 , 根据 实例
上架站对所定孔位进行复测 , 如复测 的误差在 限差 范 围之
晋城环城高速公路西北段采空 区治理工 程 合
同段位于甘润村的隆欣煤矿上部, 该矿 自 19 年开采 98 至今 , 生产能力 3 ×1 / , 0 0ta短臂式开采, 开采 3 煤层 , # 采 厚 3 回采率 5 % 。采 空 区治理深 26 3 m。 m, 0 2  ̄28 该工 程路 段 地 表 已出 现 塌 陷 现 象 , 综 合 分 析 计 经
向中部展开施工 , 钻孔分二序次进行。
() 8 钻孔施工结束后 , 应带钻孔开孔通知书、 工程报
作者简介 : 谷水兵 (9 2) 男( 18 一 , 汉族)河 南扶 沟人, , 助理工程师 , 现从事岩土工程及地质灾害治理工程技 术工作 。
14 6
西 部探矿 工 程
21 0 2年第 1 期
() 1为了避免塌孔, 钻孔成孔后应及 时进行浇注孔 口管 。 () 2将一端带有j 10 m ̄j 10 m 法兰托盘 的 2 2m 『 2 3m 『 5rm注浆 管 下 入 孔 内 变径 处 , 内投 入 2 c 厚 砾 0 a 孔 0m 石, 以堵 塞 大 的缝 隙 , 再投 入 3c 厚 的 粘 土 , 0m 防止 浆 液 渗 出 , 后灌 入 水 灰 比为 ( 然 1:1 5 ~ ( . ) 1:2 的 稠 水 泥 ) 浆 , 注长 度为 4 6 浇 ~ m。浇注 要求 达到 注浆 过程 中浆 液 不会从注浆孔孔壁 四周溢 出。水泥浆 液中应加入水泥 重量 2 的速凝剂 , 速将注浆管与孑 壁 固结 , 快 L 起到止 浆作 用 。 () 注完 毕 , 3浇 当浇 注 的水 泥 浆充分 凝 固后 , 技术 经 人员 和监理 工 程师验 收合 格 达到 止浆效 果 后 , 开始 注浆 施工 。 3 2 1 采 空 区注浆 材料 ..
高速公路路基沉降的计算方法浅析
高速公路路基沉降的计算方法浅析高速公路是现代化交通运输中的重要组成部分,大大缩短了人们的出行时间。
然而,长期使用和自然因素的影响都会对高速公路产生一定的影响。
其中最常见的问题莫过于路基沉降。
路基沉降不仅会造成驾驶安全问题,还会对道路寿命产生负面影响。
因此,对高速公路路基沉降的计算方法进行深入的分析和研究显得尤为重要。
首先,路基沉降主要是由于土体固结引起的。
固结是指岩土体在受荷后发生变形,使空隙度减小。
计算路基沉降时,通常采用垂直压缩度和剪切变形进行计算。
垂直压缩度指在一定应力状态下,土体垂直于压应力方向的应变与垂直应力之比。
而剪切变形则表示岩土体受到剪应力而产生的变形。
这两个参数的测定和计算都需要运用到一系列实验装置和计算公式。
其次,高速公路路基沉降的计算还需要考虑到不同使用条件的影响。
为了准确计算路基沉降情况,在考虑路基沉降的时候,还应该考虑到高速公路的使用率、使用时间、车流量、气候等因素的综合作用。
由于实际使用情况的复杂性,计算公式和实验方法的应用十分重要。
最后,计算高速公路路基沉降还需要依据一定的理论和实践经验。
岩土力学和土工材料力学等理论在计算中发挥着十分重要的作用。
同时,根据实际情况,针对不同的地质和气候条件,选择恰当的计算方法和技术手段,才能保证计算的准确性和实用性。
总之,高速公路路基沉降的计算方法是一个综合性问题,在计算中需要涉及岩土力学、土工材料力学等专业领域,结合使用情况和地质气候等多个因素进行综合计算。
随着科学技术的不断发展和理论的不断完善,我们相信,在不久的将来,高速公路路基沉降的计算方法将会更加成熟、准确和实用。
FLAC 3D在工程降水引起地面沉降分析中的应用
・
隧 道/ 下 工 程 ・ 地
F A 3在工程降水 引起地 面沉降 分析中的应用 L CD
李 勇 , 王 文喜
( 阳 中铁 交 通 设 计 咨 询 有 限公 司 ,沈 阳 1 0 1 ) 沈 10 3
摘
要 : 沈 阳地 铁 1号 线 青 年 大 街 站 为 工 程 背 案 , 助 三 维 以 借
不 连续分 布 的地 层忽 略 , 将地层 概化 , 概化后 的模 型含 5层土 , 场地 主要地层 特性 如下 。 ①杂 填土 : 分布 连续 , 厚 0 3~ . l层 底埋 深 层 . 4 0I, T
0. ~4. 。 3值 与 实 际观 测 值 进 行 比 将 LC
沈 阳地 铁 1 线青 年大街 站是与 规划地铁 2号 线 号
的换 乘 站 , 站 总 长 10 0m, 宽 2 . 车 9. 总 6 8m。2号 线 车
站 总长 10 0m, 4 . 总宽 2 . 18m。该 段 地势平 坦 , 面高 地
程 为 4 . 。该 站处于 本市 繁华 地 带 , 4 8I n 工程 场 地范 围
砾 石 风 化 严 重 , 胶 结 性 , 土 量 较 大 。 本 次 勘 察 未 穿 具 含 透 该 层 , 大 揭 露 厚 度 1. 层 顶 埋 深 4 . ~ 最 3 0 m, 06
4 0I。 3. l l
研究 地面沉 降 的非 常实用 的工具 。应 用其流 固耦合 力
学计 算模式 , 能够模 拟 基坑 降水 过 程及 基 坑 降水 引 发 的地 层 的三 维 固结 沉 降 。 因此 , 择 F A 软件 , 选 LC 以 沈 阳地铁 1号线青 年 大街 站 为 背景 , 降 水 引起 的地 对
FLAC_3D_在采空区稳定性分析中的应用
序号 标识
岩石 名称
1 t1 土 层 2 t2 细 砂 3 t3 中 砂 4 t4 2 煤层 5 t5 粗 砂 6 t6 砂 页 7 t7 粉 砂 8 t8 4 煤层 9 t9 砂 页
密度/ ( kg·m - 3)
1 960 2 540 2 580 1 400 2 560 2 600 2 630 1 430 2 660
(1) 根据南马圈采空区物探勘查成果图 , 采空区分布如图 2 所示. (2) 为分析采空区对上覆岩层的扰动影响 , 将南马圈矿和村民自采矿进行三维位移场 、应力场和 破坏场对比分析. (3) 结合地表特征点的高程测量结果 , 在三维数值计算模型的地表分别设置了 4 个沉陷监测点 (见图 3) , 对采动沉陷的历史过程进行监测 , 同时确定南马圈矿和村民自采矿的塌陷中心位置.
2 工程区概况
场内地层上覆地层为第四系表土 , 主要以砂 、砾卵石为主 ; 下伏基岩为二叠系至石炭系的煤系地 层 (二叠系地层主要为茂山组和荒山组) . 物探结果表明 , 场地深部存在着 1 条拟断层 , 北东 55°, 倾 向东南 45°左右 ; 岩体普遍发育近似正交的两组节理 , 间距 0. 2 m , 闭合性较好 , 可见钙质填充物. 场地地下水类型为第四系孔隙潜水 , 主要赋存于表层砂砾石层中 , 稳定水位埋深 4 m 左右. 根据钻探 取样和岩石力学实验结果 , 考虑岩石尺寸效应 , 岩石力学参数见表 1.
第 26 卷第 2 期
河南理工大学学报 (自然科学版)
Vol. 26 No. 2
2007 年 4 月 J OU RNAL OF HENAN POL YTECHN IC UN IV ERSIT Y(NA TU RAL SCIENCE) Apr. 2007
下伏采空区拟建水库的FLAC_3D_稳定性分析
收稿日期2014-03-04作者简介吴冰(1989—),男,硕士研究生。
通讯作者何文(1981—),男,讲师,博士。
下伏采空区拟建水库的FLAC 3D稳定性分析吴冰1何文1高忠1王明1柯新华2吴贵生3(1.江西理工大学资源与环境工程学院,江西赣州341000;2.江西铜业集团公司德兴铜矿,江西上饶334224;3.湖南有色金属股份有限公司黄沙坪矿业分公司,湖南郴州424421)摘要某矿山采空区上方拟建一酸性水库,为了论证其可行性,根据采场中段平面图和上部地形图,采用SUR-PAC 三维矿业软件建立采空区及地表三维模型,划分成块体后,利用SURPAC -FLAC 3D 接口程序导入FLAC 3D生成计算模型后进行计算分析,并对库区稳定性进行评价。
通过数值模拟分析可知,在现有情况下,采空区上部地表最大沉降量约为9.5mm ,在拟建水库蓄水后,地表新增沉降量约7.1mm 。
模拟结果表明:现有采空区是稳定的,在拟建水库蓄水后,采空区上部地表沉降量在建筑物允许的沉降范围内,为拟建水库的可行性提供了理论依据。
为分析爆破振动对库区影响,对爆破振动进行现场测试,得到振动最大合成速度均在质点允许振动速度安全范围内,爆破振动对拟建水库安全性并无明显影响;并在此基础上提出了一些建议措施。
关键词水库FLAC 3D数值模拟最大沉降量爆破振动振动速度稳定性中图分类号TD325,TV697文献标志码A 文章编号1001-1250(2014)-06-050-05FLAC 3D Stability Analysis of a Reservoir Planned to be Constructed over Mined-out AreasWu Bing 1He Wen 1Gao Zhong 1Wang Ming 1Ke Xinhua 2Wu Guisheng 3(1.School of Resources &Environment Engineering ,Jiangxi University of Science and Technology ,Ganzhou 341000,China ;2.Dexin Copper Mine ,Jiangxi Copper Corporation ,Shangrao 334224,China ;3.Huangshapin Mining Branch Corporation ,Hunan Nonferrous Metals Holding Group Corporation ,Chenzhou 424421,China )AbstractAn acid reservoir is planned to be constructed in a mined-out area of mine.In order to demonstrate its feasibil-ity ,a three-dimensional model of mined-out area and surface was built with 3D mining software SURPAC according to the inter-mediate entry plane of every stope and the topographic map.After divided into blocks ,the calculation model was generated in FLAC 3D by the interface program of SURPAC-FLAC 3D ,and then the stability of reservoir area is evaluated by calculation and a-nalysis.Through the numerical simulation ,maximum surface settlement on mined-out area is about 9.5mm in the current situa-tion ,and a new increased settlement is about 7.1mm after water storage in reservoir.The simulation results show that the pres-ent mined-out area is stable.The surface settlement on mined-out area is within the allowable range of building settlement after water storage in the reservoir.This provides a theoretical basis for the feasibility of building the reservoir.The field testing of blasting vibration is made to study the effect of blast vibration on stability of the reservoir area.The results show that the maxi-mum vibration velocities are all within allowable safe range ,and blast vibration has no obvious effect on safety of the acid reser-voir.Some suggested measures were proposed on this basis.KeywordsReservoir ,FLAC 3D ,Numerical simulation ,Maximum settlement ,Blast vibration ,Vibration velocity ,Stability多数金属矿山在矿体回采后,会形成采空区,造成不同程度的地面沉陷,但是由于一些矿区的面积或者地形等条件的限制,使得在采空区上方新建建筑物成为一种选择[1-3]。
采空区路基沉降计算公式
采空区路基沉降计算公式引言。
在矿区开采过程中,采空区的形成是不可避免的。
采空区对周围环境和工程设施会产生一定的影响,其中包括路基沉降。
路基沉降是指由于采空区下方地层的变形而导致路面或路基沉降的现象。
为了有效地预测和控制采空区对路基的影响,需要建立相应的计算公式。
本文将探讨采空区路基沉降的计算公式及其应用。
采空区路基沉降计算公式的建立。
采空区对路基的影响主要是通过地下水位变化和地层变形两个方面来实现的。
地下水位变化会导致路基的软化和下沉,而地层变形则会引起路基的沉降和破坏。
因此,建立采空区路基沉降的计算公式需要考虑这两个方面的影响。
首先,我们来看地下水位变化对路基的影响。
地下水位的变化会导致土壤的湿度发生变化,进而影响土壤的强度和稳定性。
根据地下水位变化对路基的影响可以建立如下的计算公式:Δh = K ×ΔG。
其中,Δh表示路基的沉降量,K为地下水位变化系数,ΔG为地下水位的变化量。
地下水位变化系数K可以通过实地观测和试验确定,其数值与地质条件和路基结构有关。
其次,地层变形对路基的影响也是十分重要的。
地层变形会引起路基的沉降和破坏,因此需要建立相应的计算公式来预测路基的沉降量。
根据地层变形对路基的影响可以建立如下的计算公式:Δh = ∑(Δhi)。
其中,Δhi表示地层变形引起的路基沉降量,∑表示对所有地层变形引起的路基沉降量进行累加。
地层变形引起的路基沉降量可以通过地质勘探和数值模拟得到。
综合考虑地下水位变化和地层变形对路基的影响,可以建立如下的采空区路基沉降计算公式:Δh = K ×ΔG + ∑(Δhi)。
该计算公式综合考虑了地下水位变化和地层变形对路基的影响,能够较为准确地预测采空区对路基的影响。
采空区路基沉降计算公式的应用。
采空区路基沉降计算公式的建立为预测和控制采空区对路基的影响提供了有效的工具。
通过对地下水位变化和地层变形的考虑,可以较为准确地预测采空区对路基的影响,并采取相应的措施进行预防和修复。
基于FLAC^(3D)的公路隧道塌方段工程处置及其优化
第30卷第2期 湖南城市学院学报(自然科学版)V ol.30 No.2 2021年3月 Journal of Hunan City University (Natural Science)Mar. 2021基于FLAC3D的公路隧道塌方段工程处置及其优化赵勇1,杨志刚1,周智辉2, 3,凌同华2,张韦华1,何文超2(1. 中铁南方投资集团有限公司,广东深圳 518000;2. 长沙理工大学土木工程学院,长沙 410114;3. 深圳市地铁集团有限公司,广东深圳 518026)摘要:以浙江省武义县白阳山隧道为工程背景,采用FLAC3D有限差分法,模拟现阶段采用的预留核心土-上下台阶法开挖下围岩的变形及应力变化规律,并从优化施工方法的角度,综合分析开挖过程中隧道围岩位移及塑性区分布情况,评价不同施工方案对围岩稳定性的影响﹒结果表明,预留核心土-上下台阶法存在洞身位移变化大的问题;预留核心土-三台阶法控制围岩位移变化效果最好,且围岩塑性区变化范围最小﹒关键词:隧道塌方;有限差分法;开挖方法;围岩稳定性;支护结构设计中图分类号:U45 文献标识码:A doi:10.3969/j.issn.1672-7304.2021.02.0001文章编号:1672–7304(2021)02–0001–07Engineering treatment and optimization of collapse section of highway tunnelbased on FLAC3DZHAO Yong1, YANG Zhigang1, ZHOU Zhihui2, 3,LING Tonghua2, ZHANG Weihua1, HE Wenchao2(1. China Railway South Investment Group Co., Ltd., Shenzhen, Guangdong 518000, China;2. School of Civil Engineering, Changsha University of Science & Technology, Changsha, Hunan 410114, China;3. Shenzhen Metro Group Co., Ltd., Shenzhen, Guangdong 518026, China)Abstract: Based on the engineering background of Baiyangshan Tunnel in Wuyi County, Zhejiang Province, this paper applied the finite difference method FLAC3D to simulate the deformation and stress variation law of surrounding rock caused by the excavation of reserved core soil and upper and lower steps method, and then from the perspective of optimizing the construction method, comprehensively analyzed the displacement and plastic zone distribution of tunnel surrounding rock during excavation and evaluated the influence on surrounding rock stability of different construction schemes. The results show that the method of reserved core soil and upper and lower steps method has the problem that the displacement of tunnel body varies greatly, while the reserved core soil and three steps method has the best effect to control the displacement change, and the plastic zone extends the smallest change range.Key words: tunnel collapse; finite difference method; excavation scheme; surrounding rock stability; supporting structure design由于隧道工程地质条件的复杂性和不可完全预知性,隧道在施工过程中很有可能会遇到塌方、涌水突泥、大变形等工程事故﹒塌方是隧道施工中最常见的安全事故[1],隧道塌方严重影响工程的进度,威胁施工人员的生命安全并带来巨大经济损失[2]﹒随着公路建设的发展,塌方事故在隧道施工过程中时有发生,已成为隧道施工过程中的一个重要安全隐患[3]﹒为了确保隧道工程建设的安全性,开展山岭隧道施工塌方研究及处治具有重要的理论价值与工程实际意义[4-5]﹒根据塌方形态的不同,隧道塌方可划分为局部塌方、拱形塌方、异形塌方、膨胀岩塌方和岩爆,不同类型的塌方对应着不同处治措施[6]﹒其中,拱形塌方按规模可以划分为2大类:在浅埋隧道中,岩体较破碎时,塌方往往延伸到顶部;在深埋隧道中,易出现摩擦桩效应,高度一般在4~20 m,塌方影响范围较大﹒拱形塌方多发生在层状岩体或碎块岩土中﹒预防拱形塌方破坏,一般采用超前锚杆支护、注浆加固等措施﹒塌方处治必须先对塌方原因有正确的认识,收稿日期:2020-08-27基金项目:国家自然科学基金项目(51678071, 51608183, 51878074)第一作者简介:赵勇(1965—),男,河南柘城人,教授级高级工程师,硕士,主要从事隧道与地下工程施工与管理研究﹒E-mail: ********************湖 南 城 市 学 院 学 报(自然科学版)2021年第2期2否则容易造成二次塌方或其他灾害,这不仅将导致经济损失加大,还可能造成安全生产事故[7]﹒塌方处治涉及的施工工艺主要有反压回填、喷浆固结、超前大管棚、超前小导管、锚杆施工和换拱,以上工艺可按工程施工实际情况灵活应用[8-9]﹒王道远等[10]在强震区隧道施工塌方段回填注浆并且设置初支减震缝和二衬减震缝,在抗震和防塌方方面产生了积极作用﹒王彬彬等[11]认为在红黏土隧道塌方事故中,采取快挖快支、及时封闭成环的原则,强化支撑体系,减少土体暴露时间,能高效快速地治理塌方﹒张旭等[12]提出了密排短管棚、三台阶预留核心土+锁脚锚管、注浆回填和洞内临时支撑的防塌方的综合处治技术,在应对卵石地层塌方中卓有成效﹒李志厚等[13]在大风垭口隧道特大塌方事故中,综合分析该隧道突泥富水和涌水险情,提出地表与洞内处治相结合,以及洞内排水的综合处治措施,并对地层应力和变形进行监测,对公路隧道特大塌方提供了治理经验﹒田志宇等[14]详细介绍了隧道塌方体的荷载计算方法、钢支撑对喷射混凝土的加强原理及隧道初期支护与二次衬砌极限承载力的计算方法等,并给出了塌方体支护参数的选取方法﹒由于隧道工程的复杂性,若能在施工时采取稳定性较强的方法,将对提升隧道施工安全性产生重要作用﹒本文从隧道施工方法着手,结合白阳山隧道的地质情况,通过数值模拟和现场监控量测,比较了预留核心土-上下台阶法、预留核心土法、预留核心土-三台阶法3种方法的稳定性﹒1 白阳山隧道工程概况白阳山隧道是浙江省金华市永康至武义公路改建工程(阳丰至湖塘沿段)的重点和难点之一﹒该隧道为分离式隧道,左洞长480 m,右洞长563 m﹒右洞进出口桩号YK15+255~YK15+818,采用上下台阶法开挖﹒本文主要针对白阳山隧道右洞YK15+380~YK15+400塌方段施工开展数值模拟分析,其断面宽为16.4 m,高为9.6 m﹒白阳山隧道穿越垄状丘陵,山体呈东北走向,线路与山脊正交,丘顶浑圆,丘顶高层约235 m﹒进洞口段缓坡堆积含碎石粉质黏土,厚1.0~2.0 m,下伏中风化英安玢岩侵入岩,围岩以碎裂结构-镶嵌碎裂结构为主,洞室自稳能力差,洞身段岩性复杂﹒围岩及支护结构物理力学参数见表1.表1 围岩及支护结构物理力学参数围岩及支护结构弹性模量E/GPa 泊松比μ密度ρ/(kg·m−3) 黏聚力c/MPa 内摩擦角φ/(°) Ⅴ级围岩 1.1 0.39 2000 0.07 23 Ⅳ级围岩 1.8 0.33 2300 0.42 29 锚杆 210.0 0.31 7800 — — 初期支护 25.7 0.20 2500 — —2 隧道塌方段情况及处治措施自2019-01-21白阳山隧道右线掌子面掘进里程至YK15+380时,隧道拱顶部位不断有小块碎石掉落,使钢筋网破损以及钢拱架发生变形﹒从现场监控量测数据可知,YK15+370~YK15+390段围岩变形速率和累计变形量都异常大,同时拱腰部位混凝土表面多处出现开裂剥落现象﹒2019-01-27,施工方已经停止爆破作业,YK15+380拱顶部位出现长6 m,纵深3 m的坍塌洞,为防止塌洞进一步恶化,采用超前小导管注浆填实,封闭掌子面,待围岩释放压力再进行处理;2019-03-18,施工方施作第2次注浆;2019-04-01,业主与设计方协商决定对YK15+380~YK15+400段提高初期支护强度,采用注浆填筑塌腔、变更钢拱架强度与间距、缩短进尺深度等处理措施﹒3 工作面坍塌采用FLAC3D模拟分析本文从沉降、收敛位移、围岩应力、围岩塑性区分布等方面,探讨分离式隧道塌方经处理后各项指标的变化规律﹒根据工程的实际工况,结合地下结构的计算原理,模型边界条件选取范围为:隧道水平方向114 m,轴线方向70 m,隧道仰拱距模型底部25 m,上表面为自由面﹒考虑地表面坡度较平缓,建立模型过程中将地表简化为等坡脚计算﹒在ANSYS中前期建立模型后,将模型导入FLAC3D 中处理计算,建立的三维数值计算模型如图1所示﹒初期支护采用shell单元,锚杆采用cable单元,其他均采用实体单元模拟﹒屈服准则选用赵勇,等:基于FLAC3D的公路隧道塌方段工程处置及其优化第30卷3Mohr-Coulomb准则,将围岩视为各向同性介质,忽略表层风化影响,沿隧道纵向考虑围岩级别变化分别设定力学参数﹒图1 三维数值计算模型3.1 隧道拱顶下沉结果分析图2为预留核心土-上下台阶法的竖向位移云图,图3为预留核心土-上下台阶法开挖施工步与拱顶竖向位移关系﹒图2 预留核心土-上下台阶法竖向位移云图图3 预留核心土-上下台阶法开挖施工步与拱顶竖向位移关系由图2和图3可以看出,在白阳山隧道开挖时,Z方向沉降值主要集中在拱顶区域,左洞的整体沉降值略大于右洞﹒由于对隧道洞身塌方段进行空腔注浆及延缓开挖进尺,这在一定程度上可控制围岩进一步变形破环﹒隧道左右洞底部均有一定隆起但差别不大,且隆起量较拱顶下沉量小,这是因为隧道底部的隆起主要由隧道开挖后应力释放而产生的弹性变形,而拱顶下沉则包含了部分受扰动围岩的塑性沉降﹒3.2 隧道洞周收敛结果分析图4为预留核心土法-上下台阶法开挖水平位移云图,图5为预留核心土法-上下台阶法开挖水平位移变化﹒图4 预留核心土法-上下台阶法开挖水平位移云图图5 预留核心土法-上下台阶法开挖水平位移变化由图4可知,在白阳山隧道开挖过程中,X 方向位移主要集中在拱腰区域,左洞拱腰累计位移值为2.52 mm,右洞拱腰位移值为2.23 mm,左洞的整体收敛值略大于右洞﹒拱顶水平位移变化均较少,不足0.3 mm﹒围岩水平位移变化呈对称分布,左侧拱腰部位水平应力最大,右洞右侧壁边墙部位水平位移值减少明显﹒由图5可知,左洞拱顶水平位移值为0.07 mm,右洞拱顶水平位移值为0.03 mm,两处水平位移值基本保持平稳趋势﹒由此可知,在施工阶段不需要考虑拱顶水平位移的影响﹒拱腰水平位移在下台阶开挖过程中变化明显,左洞拱腰水平位移为1 mm,右洞拱腰水平位移为0.8 mm,而在围岩开挖支护过程中水平位移值都比较少﹒这表明,下台阶开挖后,封闭仰拱有效地阻止了围岩向内挤压产生破坏,并减少了洞室的水平位移﹒3.3 围岩应力状态分析隧道开挖破坏了岩体平衡状态,使围岩发生卸荷回弹和应力重分布,从而产生新的应力平衡状态﹒分析施工过程中围岩应力分布及位移变化,湖 南 城 市 学 院 学 报(自然科学版)2021年第2期4对实际施工中支护结构设计参数选取及围岩稳定性判断具有重大意义﹒图6和图7分别为围岩最大主应力云图﹒由图6与图7可知,隧道开挖过程中,围岩应力分布呈动态变化,左洞拱顶与拱底小范围出现拉应力﹒其中,左洞拱顶部位最大拉应力值为0.16 MPa,拱底出现的最大拉应力值为0.2 MPa,这些应力值均未超过混凝土的抗拉强度;拱腰部位出现的压应力为3.2 MPa;侧墙及边墙部位出现的压应力为2.5 MPa﹒这说明上台阶开挖对围岩应力重分布影响最大﹒图6 围岩最小主应力云图图7 围岩最大主应力云图3.4 围岩塑性区分析隧道开挖后,稳定塑性区为未出现松动破坏的塑性区,这类围岩塑性区比较稳定,被破坏的可能性小,并且能够作为支护结构承担围岩的部分荷载;非稳定塑性区为出现松动破坏后的塑性区,这类塑性区会破坏围岩稳定性,极易造成突发坍塌破坏﹒因此,施工中不允许出现非稳定塑性区﹒在一定范围内的稳定塑性区对隧道施工支护结构有利,它可以分担部分初期支护承担的围岩压力,有利于保持围岩稳定性﹒图8为隧道施工完成后的围岩塑性区分布﹒图8中以不同颜色代表不同的破坏机制:红色和天蓝色表示单元受剪切破坏;粉红色和深蓝色表示单元受拉伸破坏﹒由图8可知,隧道开挖支护完成后,洞室周边仅在4 m范围内塑性区分布明显,这说明施工支护措施有效限制了塑性区的扩展,也表明采用4 m长的砂浆锚杆加固围岩符合工程实际﹒在施工过程中,围岩大部分都发生了剪切破坏,但只在隧道拱顶及仰拱底部1 m范围内外岩体曾经发生过拉伸破坏﹒图8 隧道施工完成后围岩塑性区分布4 施工方案对比及优化现场监控量测数据表明塌方段洞身位移变化较大,这说明现场实际施工方案存在一定的风险﹒隧道塌方处理时,保证隧道围岩稳定性至关重要,为了防止产生二次塌方,应采取合适的施工开挖方法,并控制进尺深度﹒因此,有必要对隧道塌方段施工方案进行优化分析﹒根据国内外隧道施工方法,同时考虑白阳山隧道工程具有大断面及复杂地质情况的特点,提出3种不同的施工方案进行比较﹒①预留核心土-上下台阶法(原施工方法):先行洞采用预留核心土开挖,后行洞采用上下台阶法开挖﹒②预留核心土法:左右洞均采用预留核心土开挖﹒③预留核心土-三台阶法:先行洞采用预留核心土开挖,后行洞采用三台阶法开挖﹒本文着重分析不同施工方案开挖对隧道围岩稳定性的影响,并运用FLAC3D数值模拟,对位移的分布规律、围岩应力和塑性区变化情况进行分析﹒4.1 沉降位移对比分析为了研究不同开挖方式下隧道塌方段围岩的稳定性,对典型断面YK15+380与YK15+390进行了研究,即在模拟断面Z=35与Z=45布设关键点进行动态监测﹒其中,Z=35断面在不同开挖方案下的关键点位移变化曲线如图9所示﹒赵勇,等:基于FLAC3D的公路隧道塌方段工程处置及其优化第30卷5图9 Z=35断面在不同开挖方案下关键点位移变化曲线由图9可知,从拱顶下沉看,预留核心土-上下台阶法的拱顶下沉为−10.23 mm;预留核心土法的拱顶下沉量为−10.15 mm﹒这表明有无核心土开挖对拱顶位移量影响不大,且其初期支护时间基本相当﹒从拱腰竖向位移可以看出,预留核心土-上下台阶法、预留核心土法、预留核心土-三台阶法的拱腰沉降量分别为−2.22,−2.30和−2.02 mm﹒从竖向沉降来看,3种开挖方案中预留核心土-三台阶法最好,收敛量较其他2种方法降低了10%﹒隧道YK15+380~YK15+430段围岩软弱破碎,为研究塌方腔体注浆区域对掌子面开挖的影响,取其2个特殊断面对围岩稳定性进行分析﹒其中,一个断面为注浆段;另一个断面为掌子面开挖断面﹒文中重点分析在不同开挖方案下这2个典型断面的关键点竖向位移最终值,其值如表2所示﹒表2 不同开挖方案下典型断面的竖向位移最终值mm断面位置拱顶竖向位移拱腰竖向位移(a) (b) (c) (a) (b) (c)Z=35(右洞) −10.23 −10.15 −8.51 −2.22 −2.30 −2.02 Z=45(右洞) −11.75 −11.92 −9.12 −2.31 −2.45 −2.17 注:(a)表示预留核心土-上下台阶法;(b)表示预留核心土法;(c)表示预留核心土-三台阶法表2中,Z=35(右洞)为隧道塌腔注浆区域,Z=45(右洞)为掌子面软弱围岩开挖区域﹒由表2可知,在Z=45(右洞)断面,预留核心土法对限制隧道拱顶下沉的效果甚微,而采用预留核心土-三台阶法开挖对拱顶下沉或拱腰竖向位移增幅较小,相对于Z=35(右洞)断面,仅分别增加0.61 mm和0.15 mm,拱顶最终沉降量为9.12 mm,拱腰最终沉降量为2.17 mm,均小于规定的允许沉降值﹒因此,从控制拱顶沉降考虑,应该优先选用预留核心土-三台阶法﹒4.2 洞周收敛对比分析图10为在不同开挖方案下Z=35断面关键点的水平位移变化曲线,表3为在不同开挖方案下典型断面水平位移最终值﹒由表3可知,从拱腰水平位移看,预留核心土-上下台阶法、预留核心土法和预留核心土-三台阶法的隧道拱腰水平位移分别为−2.34,−2.37和−1.76 mm﹒故从控制水平位移角度,应优先选择预留核心土-三台阶法﹒从图10可知,在Z=35(右洞)断面,在3种开挖方案下,预留核心土-三台阶法施工的隧道断面关键点水平位移值最小,因此该法对预防塌方有积极作用﹒图10 在不同开挖方案下Z=35断面关键点的水平位移变化曲线表3 不同开挖方案下典型断面水平位移最终值mm断面位置拱顶水平位移拱腰水平位移(a) (b) (c) (a) (b) (c)Z=35(右洞) −0.03 −0.02 −0.01 −2.34 −2.37 −1.76 Z=45(右洞) −0.05 −0.04 −0.02 −2.61 −2.56 −1.88 注:(a)表示预留核心土-上下台阶法;(b)表示预留核心土法;(c)表示预留核心土-三台阶法4.3 围岩应力对比分析选取Z=35与Z=45断面,分析在不同开挖方案下其隧道围岩最大主应力分布情况,2个断面围岩最大主应力见图11~图13﹒湖 南 城 市 学 院 学 报(自然科学版) 2021年第2期6(a)Z =35断面围岩最大主应力 (b)Z =45断面围岩最大主应力图11 预留核心土-上下台阶法典型断面围岩的主应力云图(c)Z =35围岩断面最大主应力 (d)Z =45围岩断面最大主应力图12 预留核心法典型断面围岩的主应力云图(e)Z =35断面围岩最大主应力 (f)Z =45断面围岩最大主应力图13 预留核心土-三台阶法典型断面围岩的主应力云图由图11~图13可知,预留核心土-三台阶法、预留核心土法和预留核心土-上下台阶法开挖引起的围岩主应力分布规律较为相似,且在隧道开挖支护完成后,围岩均出现了一定程度的应力集中现象﹒4.4 围岩塑性区对比分析在隧道开挖支护过程中,围岩会出现应力集中现象,如当某区域的应力达到屈服应力,或超过围岩的抗剪强度,隧道周边岩体会产生塑性破坏﹒围岩塑性区是判别围岩稳定性的重要指标﹒采用FLAC 3D 软件对不同开挖支护方案下隧道围岩塑性区分布状况进行模拟分析,可判断围岩是否处于稳定状态﹒施工方案不同,对围岩的扰动程度也不一样,(a)预留核心土-上下台阶法 (b)预留核心土法 (c)预留核心土-三台阶法图14 不同开挖方式下隧道围岩塑性区分布赵勇,等:基于FLAC3D的公路隧道塌方段工程处置及其优化第30卷7从而隧道周边岩体塑性区分布情况也会不一样(见图14)﹒由图14可知,3种开挖方式下,预留核心土-上下台阶法与预留核心土法开挖时,隧洞周边拉应力分布更广泛,隧洞围岩更易发生受拉破坏;预留核心土-三台阶法开挖所对应的隧道周边受拉范围较小,即采用该法开挖时,更利于保持围岩稳定性﹒5 结论1)运用有限差分软件FLAC3D对3种施工方案进行模拟分析,3种方案均表现为拱顶沉降最大,仰拱隆起最明显﹒预留核心土-三台阶法控制围岩位移效果最好,其次是预留核心土法,预留核心土-上下台阶法最差﹒2)采用不同施工方案开挖时,围岩塑性区主要分布于右洞拱脚处、拱顶和拱肩区域﹒预留核心土-三台阶法的塑性区变化范围最小,预留核心土法次之,预留核心土-上下台阶法最差﹒3)塌方段开挖方案由预留核心土-上下台阶法改为预留核心土-三台阶法,可以明显控制围岩位移变化,减少塑性区扩展,从而保障隧道施工安全﹒4)根据白阳山隧道右线塌方事故产生的具体原因,确立了处治该塌方事故的合理方案,即建议隧道短进尺施工,掌子面逐榀开挖,及时施作初期支护,并对拱顶进行注浆加固﹒参考文献:[1]马士伟, 韩学诠, 廖凯, 等. 大断面软弱围岩隧道防塌方实时监测预警标准研究[J]. 铁道工程学报, 2014, 31(9):88-92.[2]姜安民, 董彦辰, 张晓波, 等. 基于粗糙集与模糊理论的公路隧道塌方风险评价[J]. 公路与汽运, 2020(3): 156-160. [3]周宗青, 李术才, 李利平, 等. 浅埋隧道塌方地质灾害成因及风险控制[J]. 岩土力学, 2013, 34(5): 1375-1382.[4]曹文贵, 翟友成, 王江营, 等. 山岭隧道塌方风险的集对分析方法[J]. 中国公路学报, 2012, 25(2): 90-99.[5]张胜, 何文超, 黎永索, 等. 小净距隧道洞口段地表下沉分段曲线拟合研究[J]. 湖南城市学院学报(自然科学版), 2019,28(5): 1-5.[6]郑陈旻. 隧道塌方预警预测体系及治理措施研究[D]. 上海:同济大学, 2008.[7]韦盛. 小头岭隧道塌方原因分析和处治措施研究[D]. 西安:长安大学, 2019.[8]刘智勇, 尹登华, 吴学智, 等. 某公路隧道塌方原因分析及处治技术[J]. 路基工程, 2018(1): 235-239.[9]解超平. 全断面穿越昔格达地层高速公路隧道塌方的预防和处治[D]. 成都: 西南交通大学, 2015.[10]王道远, 崔光耀, 袁金秀, 等. 强震区隧道施工塌方段震害机理及处治技术研究[J]. 岩土工程学报, 2018, 40(2):353-359.[11]王彬彬, 赵小凯, 孙熔正. 海德隧道红黏土地段溶洞超前预支护施工技术[J]. 公路, 2018, 63(10): 158-162.[12]张旭, 龚娟, 李博融, 等. 砂卵石地层隧道塌方机理及处治技术研究[J]. 公路, 2018, 63(5): 318-323.[13]李志厚, 杨晓华, 来弘鹏, 等. 公路隧道特大塌方成因分析及综合处治方法研究[J]. 工程地质学报, 2008, 16(6):806-812.[14]田志宇, 邓承波, 林国进. 隧道塌方处治支护参数选取[J].铁道建筑, 2012(6): 54-57.(责任编校:付贵海)。
采空区剩余沉降量FLAC3D数值模拟计算方法
采空区剩余沉降量FLAC3D数值模拟计算方法
陈琦;陈朝清;李森林
【期刊名称】《中国新技术新产品》
【年(卷),期】2009(000)010
【摘要】采空区剩余沉降量计算是确定采空区危害程度的一种定量的方法,其可靠性、准确性和可探作性直接关系到采空区治理质量的好坏.本文通过FLAC3D数值计算方法对西沟新村基建下覆采空区剩余沉降量进行计算,经过与实测资料的对比,结果相一致,说明该计算方法准确可靠,可为基建下覆采空区质量检测以及剩余沉降量的计算方法提供有益的借鉴.
【总页数】2页(P92-93)
【作者】陈琦;陈朝清;李森林
【作者单位】台州市交通勘察设计院,浙江,台州,318020;福建永福工程顾问有限公司,福建,福州,350003;中勘冶金勘察设计研究院有限责任公司,河北,保定,071051【正文语种】中文
【中图分类】TU7
【相关文献】
1.堆载预压工程最终沉降量的计算方法研究 [J], 张立伟;杨帅东;黄志怀
2.海防堤工程沉降量及挡墙稳定计算方法探析 [J], 陈诗
3.长壁老采空区残余沉降计算方法研究 [J], 邓喀中;谭志祥;张宏贞;范洪冬;张立亚
4.基于FLAC3D数值模拟的老采空区剩余沉降量分析计算--以山东枣庄安博化工项
目为例 [J], 臧浩;种衍飞;冯堂武;于树宾;张昌龙
5.基础最终沉降量三维计算方法 [J], 田跃平;高子坤
因版权原因,仅展示原文概要,查看原文内容请购买。
FLAC3D中采场采空区的处理方法探讨
FLAC3D中采场采空区的处理方法探讨
张磊
【期刊名称】《山西焦煤科技》
【年(卷),期】2017(041)011
【摘要】针对FLAC3D在模拟采场采空区时所面临的问题,在对其软件特征分析基础上,得到了对采场采空区的不同方法,本文对FLAC3D处理采空区的不同方法进行了优缺点及适用性分析,并利用其内置FISH语言进行实现,最后以中煤大屯矿姚桥煤矿地质条件为实例,对不同模拟方法效果对比分析,结果表明,分步开挖软化顶板法并不能准确反映实际采场顶板运动规律,软化采场处理方式较为准确易行.
【总页数】5页(P26-30)
【作者】张磊
【作者单位】西山煤电集团公司官地矿,山西太原 030022
【正文语种】中文
【中图分类】TD32
【相关文献】
1.伏萨尔梁法修正及其在露天采场采空区处理中的应用 [J], 余晨晖
2.基于FLAC3D模拟的大采高采场支承压力分布规律研究 [J], 胡国伟;靳钟铭
3.露天采场高陡边坡下深部采空区的爆破处理研究 [J], 解治宇; 李翰林; 房洪亮; 韩连生; 张德辉; 沙成满
4.伏萨尔梁法修正及其在露天采场采空区处理中的应用 [J], 余晨晖
5.FLAC3D在地下矿山采场稳定性分析中的应用 [J], 罗黎明
因版权原因,仅展示原文概要,查看原文内容请购买。
高速公路下伏采空区剩余沉降量FLAC_3D_计算方法_王生俊
2 采空区剩余沉降量 FLAC3D计算方法
FLAC3D 是由美国 Itasca Consulting Group,Inc. 为岩土工程应用而开发的连续介质显式有限差分计 算机软件。该软件主要适用于模拟计算岩土体材料 的力学行为及岩土材料达到屈服极限后产生的塑性 流动,对大变形情况应用效果更好;可有效地模拟 各种开挖工程或施加支护工程等过程。软件自身设 计多种结构元素,可直接模拟这些加固体与岩(土) 体的相互作用。计算所采用的数学模型是根据弹塑 性理论的基本原理(应变定义、运动定律、能量守衡 定律、平衡方程及理想材料的连续性方程等)而建 立的。本文采用莫尔–库仑弹塑性材料本构模型对 郑少高速公路下伏采空区的剩余沉降量进行了计 算,并依据计算结果采取相应的治理方案对下伏采 空区进行了治理。 2.1 模拟方案
(1. 兰州大学 资源环境学院,甘肃 兰州 730000; 2. 石家庄经济学院 地下工程系,河北 石家庄 050031;3. 天津城市建设学院 土木工程系,天津 300381)
摘要:采空区剩余沉降量计算是确定采空区危害程度的一种定量方法,其可靠性、准确性和可操作性直接关系到
采空区治理方案的选择与评价。根据采空区塌陷冒落带的力学特征,利用横波参数或瑞雷波参数与岩体静力学参
的各层岩土介质的参数如表 1 所示。
2.3 计算结果
通过数值模拟计算,得到采空区及其上覆岩土
层的位移场和应力场,结果见图 1,2 和表 2。根据
其分布规律可得到以下结论:
(1) 研究区在不进行治理状态下,最大沉降可
高速公路下伏多层采空区地表沉陷数值模拟及预测研究
高速公路下伏多层采空区地表沉陷数值模拟及预测研究近年来,随着我国国民经济的不断发展,经济建设对交通运输体系要求的不断提高,高速公路和高速铁路建设也以前所未有的速度推进。
由于道路网越来越密集,一些地区的高速公路和铁路不可避免地需要穿越或建立在老采空区地表。
北京至哈尔滨高速公路(简称京哈高速公路)四平至长春段K70+500—K75+000段位于刘房子煤、膨润土矿区北西侧。
矿山开采形成的采空区,导致高速公路局部路面变形,出现裂缝,道路路基下沉。
本文依托“京哈高速公路四平至长春段改扩建工程刘房子路面塌陷段勘查”项目,在阅读和总结大量前人研究文献基础上,以四平至长春段K73+100—K73+500段路面塌陷为研究对象,完成本论文。
首先分析了研究区区域构造地质条件,煤系地层工程地质条件和水文地质条件,煤层厚度、顶底板岩石物理力学条件,煤矿开采历史及现状和采空区范围。
根据野外现场调查、钻探和物理勘探资料,分析研究区地表沉陷特征、分布规律及其形成机理。
然后根据研究区刘房子煤矿调查资料,运用MATLAB软件,以概率积分法为理论基础,编制多层任意形状采空区地表沉陷预计系统,并将系统用于研究区高速公路地表沉陷计算。
利用有限差分法FLAC数值模拟软件,研究汽车动荷载作用下采空区地表变形规律。
采用改进的遗传算法求解和优化Holt-Winters模型的组合参数,建立高速公路沿线监测点地表沉陷的预测模型。
最后对研究进行总结并提出今后研究的建议。
论文包括以下八个部分:绪论、研究区地质环境条件、高速公路地表沉陷特征及其机理分析、多层采空区地表沉陷预计系统设计、研究路段地表沉陷预计分析、汽车动荷载作用下采空区地表变形数值模拟、采空区地表沉陷非线性预测模型研究及结论和展望。
本文建立了基于概率积分法的多层任意形状采空区地表沉陷预计模型,利用MATLAB编制程序并制作界面窗口,形成了基于概率积分法的开采沉陷预计系统。
以研究路段下地质采矿条件为原型设计模型,研究采空区埋藏深度、采空区尺寸、汽车荷载速度对采空区上覆岩层及高速公路地面变形的影响。
乌鲁木齐市过境段公路下伏采空区地表变形计算分析
乌鲁木齐市过境段公路下伏采空区地表变形计算分析发布时间:2022-09-14T09:39:38.208Z 来源:《新型城镇化》2022年18期作者:吴景伟1 王立鹤2 [导读] 国道216线乌鲁木齐市过境段公路工程位于乌鲁木齐市东侧,道路等级为城市快速路,线路全长约43.029km。
1. 中交第二航务工程局有限公司武汉 430040;2. 新疆交通建设管理局项目执行三处乌鲁木齐市 830049摘要:本文针对乌鲁木齐市过境段公路沿线各采空区仍处于不稳定状态,其剩余变形量均较大,地下尚存在许多空洞的现状,通过对相关评价标准的分析,采用经验公式及调查区覆岩结构特点分析预测地表变形,通过计算得出采空区地表最大剩余下沉量在384至887mm之间,大于规范所规定的沉降量,会使公路路基产生裂缝及沉陷,影响着公路运营安全。
关键词:道路工程;采空区;地表变形;计算1 引言国道216线乌鲁木齐市过境段公路工程位于乌鲁木齐市东侧,道路等级为城市快速路,线路全长约43.029km。
经现场勘察,道路经过丰源煤矿的采空区(K11+100-K11+260)、新疆兵团天发能源有限公司第一煤矿的军区后勤煤矿采空区(K11+260~K11+560)和八道湾煤矿采空区(K11+560~K11+860)。
道路工程在采空区路段以路基通过,路基设计宽度为70m,自然地面算起的路基填筑高度为2.4m,边坡坡率为1:1.5。
公路沿线第四系松散层分布于区内,岩性上部是黄土状粉土,其下部是砂砾石,厚度为5~50m,具有湿陷性。
公路沿线煤矿区属于急倾斜、近距离开采,采空区上覆岩层的稳定性差,易发生变形,见图1、2所示。
图1 矿区塌陷坑群图2 路线沿线圆形塌陷坑2 相关评价标准目前,国内《建筑地基基础设计规范》(GB50007-2011)、《岩土工程勘察规范》(GB50021-2001)未对高速公路煤矿采空区上覆特大桥的变形要求有明确的规定;高速公路主要参照煤炭部颁发的《建筑物、水体、铁路及主要井巷道煤柱留设与压煤开采规程》(2016年),确定公路下覆采空区治理范围和保护煤柱的范围。
下伏采空区高速公路隧道变形特征数值分析
下伏采空区高速公路隧道变形特征数值分析1 引言高速公路隧道为线性工程,由于采空区覆岩的不均匀沉降变形,很容易造成隧道路面错落起伏、衬砌结构变形、开裂等破坏现象,对高速公路隧道的正常运营以及破坏后的补救修复,构成了严重的安全隐患和潜在危害[1-2]。
由文献检索知,国内外在高速公路隧道下伏采空区的治理标准、评价依据、理论计算及处治技术等方面尚处于经验摸索阶段、没有系统成熟的理论可供借鉴[3-5]。
目前,基于技术及经济原因,公路下伏采空区的治理,主要以注浆处治技术为主。
由于地下采空区具有隐伏性强、空间分布特征规律性差、采空区顶板冒落塌陷情况难以预测等特点[6],因此,高速公路隧道下伏采空区的潜在危害性评价和处治技术是目前具有挑战性的难题。
山西省煤炭资源在全省大部分地区均有分布,经长期的采掘活动,留下了大范围的采空区,给省内的公路建设增加了额外的难度和费用。
随着高速公路的不断发展和建设,高速公路隧道工程不可避免地要跨越或穿越采空区。
本文结合青岛至银川高速公路中山西柳林段大雨亮隧道下伏采空区治理工程,对穿越采空区的高速公路隧道变形响应特征进行数值模拟分析,探讨高速公路隧道下伏采空区危害性评价及对其进行有效处治的途径。
2 工程概况青岛至银川国道主干线山西省离石至军渡段高速公路简称离军线,位于山西省柳林县境内,全长38.542 km。
线路途经的大雨亮隧道下伏青龙煤矿采空区。
青龙煤矿采空区分布于离石至军渡段高速公路:右线RK17+415~RK17+840(425 m);左线LK17+490~LK17+840(350 m)。
该煤矿为村办煤矿,1975年开采至2004年底停采。
斜井开拓,开采10#煤层,年实际生产能力1万t,采厚1.8 m,煤层埋深40~140 m,煤层顶板为黑灰色泥岩,底板为灰色细砂岩。
路线设计标高以下采空区埋深53~97 m,采深采厚比30~54,岩层倾角近水平,采区少量淋水,煤柱支撑,仓房式开采,回采率40 %。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第24卷 第19期岩石力学与工程学报 V ol.24 No.192005年10月 Chinese Journal of Rock Mechanics and Engineering Oct.,2005收稿日期:2005–02–25;修回日期:2005–05–08作者简介:王生俊(1971–),男,2000年于兰州大学地质工程专业获硕士学位,现为博士研究生,主要从事岩土工程和道路桥梁等检测、养护和加固等方面的研究工作。
E-mail :shengjunw@ 。
高速公路下伏采空区剩余沉降量FLAC 3D 计算方法王生俊1,贾学民2,韩文峰1,3,崔素敏1(1. 兰州大学 资源环境学院,甘肃 兰州 730000;2. 石家庄经济学院 地下工程系,河北 石家庄 050031;3. 天津城市建设学院 土木工程系,天津 300381)摘要:采空区剩余沉降量计算是确定采空区危害程度的一种定量方法,其可靠性、准确性和可操作性直接关系到采空区治理方案的选择与评价。
根据采空区塌陷冒落带的力学特征,利用横波参数或瑞雷波参数与岩体静力学参数相类比,确定工程岩体参数评价指标。
针对采空区冒落带当前特征,使用FLAC 3D 软件对高速公路下伏采空区剩余沉降量进行计算,并与当前几种采空区剩余沉降计算方法相比较,阐述其优越性。
采空区治理后的实测数据表明:该计算方法准确可靠,可为高速公路下伏采空区剩余沉降定量计算提供一条新途径。
关键词:公路工程;采空区;剩余沉降;计算方法;FLAC 3D ;岩体参数中图分类号:U 412.36+6 文献标识码:A 文章编号:1000–6915(2005)19–3545–06A FLAC 3D METHOD FOR CALCULATION OF RESIDUAL SUBSIDENCE INA MINING AREA BENEATH AN EXPRESSWAYWANG Sheng-jun 1,JIA Xue-min 2,HAN Wen-feng 1,3,CUI Su-min 1(1. College of Earth and Environment Sciences ,Lanzhou University ,Lanzhou 730000,China ;2. Department of Underground Engineering ,Shijiazhuang University of Economics ,Shijiazhuang 050031,China ;3. Department of Civil Engineering ,Tianjin Institute of Urban Construction ,Tianjin 300381,China )Abstract :The calculation of residual subsidence in mining area is a quantitative method to evaluate hazards rank derived from the mining area. The reliability ,accuracy and operation of the calculation directly affect the choice and assessment of the scheme for reinforcing mining areas. Based on the mechanical characteristics of collapse ,the parameters of traverse wave or Rayleigh against static rock mass are used to determine assessment indices for the parameters of engineering rock mass. Under the current status of collapse of a mining area ,FLAC 3D software is used to calculate and evaluate the residual subsidence of the mining area beneath an expressway ;and the advantage of the method is discussed as compared to other calculation methods of residual subsidence. Through a case study of a reinforced mining area ,the testing data show that the proposed method is accurate and reliable ,and it can provide a new way to evaluate and calculate the residual subsidence in mining area beneath an expressway.Key words :highway engineering ;mining area ;residual subsidence ;calculation method ;FLAC 3D ;rock mass parameters1 引 言采空区塌陷沉降是危害高速公路建设的重大工程地质灾害问题,但其勘察、评价、治理难度较大。
采空区的塌陷,除了开采过程中引起的采矿塌陷 外,在采矿沉陷稳定后,由于地下尚存在未充填的空间,在某些特定条件和诱发因素下,地面仍然会• 3546 • 岩石力学与工程学报 2005年产生沉降变形或者突然塌陷。
剩余沉降量是自工程建设开始以后地表所产生的沉陷。
剩余沉降量包括部分已采采空区对地表影响还没有基本结束的沉降量和残余沉降量2部分[1]。
残余沉降量指地表下沉基本结束以后,破碎岩块间微小空间经过漫长时间进一步密实后的地表下沉量。
20世纪70~90年代,在调查和经验的基础上,文[2~6]先后研究了采空塌陷、下伏空洞对公路的危害。
文[7]综述了高速公路下伏采空区国内外研究现状及进展。
文[8]将与空洞接触的岩土简化为理想的线弹性介质,研究了在弹性极限条件下岩土的沉陷破坏机理。
文[9,10]采用相似材料模拟试验,研究了老采空区上方兴建建筑物后,地基的沉降规律。
文[11,12]采用数值仿真方法对采空区的开采沉陷进行了数值分析研究。
文[13,14]分别对铁山隧道采空区和木架山采空区治理进行了研究。
简而言之,采空区的剩余沉降量是采矿作用影响区内岩土体继续变形在地面的表象。
显然,正确认识这一力学过程和科学地确定岩土力学评价参数是解决这一问题的关键。
本文根据采空区塌陷冒落带的力学特征,利用横波参数或瑞雷波参数作为工程岩体评价指标,采用FLAC3D软件对高速公路下伏采空区剩余沉降量进行计算和评价,为采空区优化治理提供科学依据。
2 采空区剩余沉降量FLAC3D计算方法FLAC3D是由美国Itasca Consulting Group,Inc.为岩土工程应用而开发的连续介质显式有限差分计算机软件。
该软件主要适用于模拟计算岩土体材料的力学行为及岩土材料达到屈服极限后产生的塑性流动,对大变形情况应用效果更好;可有效地模拟各种开挖工程或施加支护工程等过程。
软件自身设计多种结构元素,可直接模拟这些加固体与岩(土)体的相互作用。
计算所采用的数学模型是根据弹塑性理论的基本原理(应变定义、运动定律、能量守衡定律、平衡方程及理想材料的连续性方程等)而建立的。
本文采用莫尔–库仑弹塑性材料本构模型对郑少高速公路下伏采空区的剩余沉降量进行了计算,并依据计算结果采取相应的治理方案对下伏采空区进行了治理。
2.1模拟方案选择郑少高速公路SK42+230~SK42+580采空区路段来建立计算模型。
计算模型沿路轴线长350 m,宽200 m,高速公路轴线两侧各100 m,采深为130 m。
路宽28 m,边坡坡度1∶1.75。
模型中考虑了高速公路路堤填方后的界限以及采空区地形起伏。
在剖面上,考虑了各地层的分界线和断层的分界线。
计算模型左右边界为单向位移约束边界(=xU0,=yU0),底部边界为双向位移约束边界(=xU0,=yU0,=zU0),地表为自由边界。
由于采空区冒落带和裂隙带位于地下一定深度,对采空区冒落带与充填胶结体进行现场和室内试验都很困难。
因此一般都是用现场横波波速或者面波波速测定来确定模量。
在采空区治理研究中,采用弹性波速测井法检测治理工程质量时,要求钻孔中注浆处理层的平均剪切波(横波)速smV≥160 m/s。
但是笔者在研究唐津高速公路(二期)采空区治理时发现,采空区治理前在冒落带中面波波速的最小值是205 m/s,注浆治理之后最小值为245 m/s[15]。
横波与面波有很好的相关性,且两者近似相等,说明国内过去对采空区冒落带注浆充填固结后强度评价指标具有一定的局限性。
考虑到采空区的长期稳定性,评价指标横波波速值应该适当地提高。
通过对焦晋高速公路下伏采空区治理工程面波与孔中波速测井资料的研究,结合以往的资料,文[16]建议横波波速最小值smV≥200 m/s。
根据以上资料分析可知,考虑到采空区剩余沉降的危害和工程治理特点,将采用如下3种模拟方案:(1) 采空区不进行注浆治理(采空区冒落带剪切波达到100 m/s)+筑路;(2) 采空区注浆治理(采空区冒落带剪切波达到200 m/s)+筑路;(3)采空区注浆治理(采空区冒落带剪切波达到300 m/s)+筑路。
通过3种方案的对比计算,确定采空区的治理方案。
2.2物理力学参数的选择路堤与各岩层的物理力学参数均依据工程岩性室内与现场试验数据,采空区冒落带的物理力学参数根据波速试验数据(主要为剪切波速sV值),经换算后得:(1) 弹性模量)1(22sνρ+=VE(1)(2) 体积模量)21(3)1(2)21(32svVvEK−+=−=νρ(2)第24卷 第19期 王生俊等. 高速公路下伏采空区剩余沉降量FLAC 3D 计算方法 • 3547 •(3) 剪切模量ρ2s )1(2V v EG =+=(3)式(1)~(3)中:E ,K ,G 分别为弹性模量、体积模量和剪切模量,单位均为kPa ;s V 为剪切波速,单位为m/s ;ρ为密度,单位为kg/m 3。
数值模拟采用的各层岩土介质的参数如表1所示。
2.3 计算结果通过数值模拟计算,得到采空区及其上覆岩土层的位移场和应力场,结果见图1,2和表2。