彭庄矿区村庄下条带开采地表的移动规律.pdf
第3 4卷 第8期 2 0 2 0年8月 北京测绘 B e i j i n gS u r v e y i n ga n dM a p p i n g V o l . 3 4 N o . 8 A u g u s t 2 0 2 0 引文格式刘承旭, 崔志浩, 李增鑫.彭庄矿区村庄下条带开采地表的移动规律[J].北京测绘, 2 0 2 0,3 4(8) 1 0 6 6 1 0 7 0. 犇 犗 犐1 0. 1 9 5 8 0/j . c n k i . 1 0 0 7 3 0 0 0. 2 0 2 0. 0 8. 0 1 1 [ 收稿日期] 2 0 1 9 1 2 1 7 [ 作者简介] 刘承旭(1 9 9 5-) , 男, 山东潍坊人, 硕士在读, 研究方向为地表变形监测。 犈 犿 犪 犻 犾2 2 1 0 4 5 0 1 1 1@q q . c o m 彭庄矿区村庄下条带开采地表的移动规律 刘承旭 崔志浩 李增鑫 ( 山东科技大学 测绘科学与工程学院,山东 青岛2 6 6 5 9 0) [ 摘 要] 为了研究条带开采工作面的地表移动规律, 本文以4 3 0 2工作面的已有资料和和该工作面的 实际水准及水平观测数据为基础, 对该区域地质条件下的彭庄煤矿4 3 0 2工作面该条采工况推进过程中的地 表沉陷规律进行了研究, 根据已有的理论基础沪指移动角制图并解算出地表各项移动参数, 在理论分析的基 础结合实际观测情况上运用概率积分法对矿区地表形变进行了预计。预计结果表明, 条带开采对于控制地 表移动变形有较好的效果。 [ 关键词] 地表移动;条带开采;沉陷规律;形变预计 [ 中图分类号] P 2 5 8T D 3 2 5 [ 文献标识码] A [ 文章编号] 1 0 0 7 - 3 0 0 0( 2 0 2 0)0 8 - 1 0 6 6 - 5 0 引言 煤炭是我国占比最大的能源消耗品, 近些年 虽然我国能源结构中的煤炭所占的比重有所下 降, 但仍然在我国的能源体系中处于主体地位, 是国家不可或缺的资源。在采出煤炭的过程中, 破坏了开采区域周围岩体的应力平衡, 应力得到 重新分配, 从而达到了新的平衡。在此过程中, 岩 层和地表产生连续的移动、 变形和非连续的破坏, 严重威胁矿区地表居民的人身和财产安全[ 1]。我 国是煤炭储量大国, 然而由于“ 三下” 压煤的存在煤 炭产出率偏低。条带开采法属于局部开采, 此法可 有效控制覆岩及地表移动变形、 保护地表建( 构) 筑 物, 不仅有利于煤矿安全生产、 提高煤炭采出率, 还 降低了生产成本、 扩大了矿井经济效益, 并为最大 限度地减小开采损害提供技术保障[ 2]。近年来, 高 效减沉的条带开采法正在被广泛应用到三下压煤 开采中来[ 3], 因此研究彭庄煤矿该工作面条带开采 的地表移动规律, 对彭庄煤矿储藏的三下压煤的后 续开采工作具有重要意义。 1 观测站布设 4 3 0 2工作面位于彭庄煤矿西二采区南翼, 地 面地势相对平坦, 排水干渠纵横交错, 一条从西 北至东南方向的向阳河从工作面上覆穿过, 两条 高压线路从东北至西南方向从4 3 0 2工作面西北 侧穿过, 其中一条为5 0 0k V高压线路( 郓泰线) 。 相对地面位于张官屯村南侧, 王集村北侧, 部分 在周垓村下。为了确保观测线能有效监测整个 采煤影响区及采煤区附近塌陷半流域的地表运 动和变形, 并获得采动下地表岩石运动的相关参 数。在一定条件下, 该地区布置了许多观测线, 并进行了大量的长期观测。该工作面一共布设了 一条走向线和三条倾向线, 其中倾向线I线共2 8 个测站, 该观测线主要用来定期监测观测两个条采 工作面南部以及工作面南部上方的周垓村在工作 面开采开始后的地表移动情况。E线为沿向阳河 布设的倾向线共3 7个测站, 用以观测河岸在开采 过程中的地表移动情况, 保证开采的顺利进行。F 线是位于工作面中部的倾向线共3 2个测站, 用以 观测工作面开采对附近地面的影响。H线主要位 于周垓村内是一条倾向线共4 0个测站, 实时监测 工作面开采对村内地表建筑物的影响。 2 研究区实测分析 2. 1 下沉规律分析 图1、 图2分别为走向线I线和倾向线F线 的动态下沉曲线。由于初次进行观测时工作面 第3 4卷 第8期刘承旭, 崔志浩, 李增鑫.彭庄矿区村庄下条带开采地表的移动规律 开采尚未开始, 所以从第二期数据开始统计下沉 值。开采初期各水准观测站的地表下沉值较小, 随着工作面开采工作面的推进和水准测量监测 工作的定期进行, 对周围岩体不断破坏, 地表逐 渐进入移动变形活跃期, 下沉速度越来越大, 下 沉量与此同时也逐渐增大[ 4]。工作面停采后下 沉速度迅速减小, 整个水准监测过程中地表下沉 速度呈现先增大后减小的状态。通过分析水准 观测数据可知, 在4 3 0 2 1工作面和4 3 0 2 2两个 工作面间歇条采的过程中, 由于覆岩关键层中细 砂岩对开采工作面围岩的应力平衡起着决定性 的作用, 所以地表的最大下沉量相对较小, 地表 下沉曲线形态为极不充分采动时的状态, 基本呈 岩梁挠度形态[ 5]。 图1 犉观测线下沉曲线图 由图1可以看出观测线两端距离工作面较 远的点受煤矿开采的影响较小, 下沉量较小。离 工作面越近水准观测点受煤炭开采的影响越大, 下沉量越大。最大下沉值距离第一个观测点F 2 1 的距离为5 0 0m, 基本位于两个工作面中央, 且曲 线大致关于峰值对称。图2中I 1点距离工作面 4 3 0 2 2约2 4 0m受工作面开采影响而引起的沉 陷较小, 随着点距离I 1距离的拉近, 沉陷值不断 增大。I 2 8距离工作面不足1 0 0m所以地表变形 较明显。由图1和图2可以清晰地看出第五期数 据和第六期数据下沉值相差较大。分析工作面 的开采进度以及观测规律可知, 由于工作面推进 速度较快且观测相隔时间相对较长, 所以第五期 与第六期数据下沉量相差较大[ 6]。从第十期数 据开始为停采后的观测结果, 此时两个工作面地 表移动尚处于活跃期, 地表移动还将持续, 随着 下沉速度的减缓地表移动逐渐趋于稳定。由图 1、 图2可知最大下沉值偏向于下山方向, 究其原 因是因为该工作面煤层倾角平均为1 2. 4 , 属于 缓倾斜煤层, 随着工作面的推进最大下沉值会偏 向于下山方向, 符合地表沉陷移动规律[ 7]。分析 获得的水准数据可知彭庄煤矿该条采工况下开 采引起的地表下沉较小, 产生的地表移动变形相 对较小[ 8]。 图2 犐观测线下沉曲线图 2. 2 地表移动变形的参数确定 2. 2. 1 下沉系数 彭庄矿区4 3 0 2工作面为充分采动条件下求 下沉系数。充分下沉条件下测得该工作面最大 下沉值(犠m a x) , 下沉系数狇的计算公式如下 狇=犠m a x/犿 c o s δ ( 1) 已知彭庄煤矿4 3 0 2工作面的煤层采厚为 3 . 0 7m, 平均采深为5 5 0m, 采空区倾斜方向的实 际长度为4 2m, 采空区走向方向的实际长度为 12 3 0m, 煤层倾角为1 2. 4 , 实测最大下沉值为 4 1 4mm。 采动系数为采空区倾斜或者走向方向所开 采的实际长度与工作面充分采动相应的倾斜或 者水平方向的最小长度之比。倾向和走向方向 采动系数的计算公式表式为 狀1=犓1(犇1/犎0) ( 2) 狀2=犓2(犇2/犎0) ( 3) 式中犇1为采空区倾斜方向的实际长度;犇2为 采空区走向方向的实际长度;犎0为平均采深; 犓1、犓2一般取0. 8。 经过计算得 狀1=0. 0 6 1;狀2=1. 7 8。 由采动系数可知彭庄煤矿4 3 0 2工作面沿走 向方向达到充分采动, 而该工作面沿倾向方向 未达到充分采动, 所以计算公式为 狇=犠m a x/犿c o sα 狀1狀槡2=0. 1 9 ( 4) 所以,4 3 0 2工作面的下沉系数为0. 1 9。 2. 2. 2 水平移动系数 由最大水平移动的计算公式犝0=犫犠0可 推得 犫=犝0/犠0 ( 5) 7601 北京测绘第3 4卷 第8期 图3 4 3 0 2工作面求参结果图 式中 犝0最大水平移动值;犠0最大下沉值; 犫水平移动系数。 利用在该工作面长时间定期观测获得的数 据, 分 析 得 到 该 工 作 面 的 最 大 水 平 移 动 值 为 1 2 2 . 3mm, 所有水准观测点中的最大下沉值为 4 1 4mm。水平移动系数计算公式为 犫=犝0/犠0=1 2 2. 3/4 1 0=0. 3 0 ( 6) 所以, 该工作面的水平移动系数为0. 3 0。 2. 2. 3 主要影响半径和主要影响角正切 主要影响半径指的是矿区地表变形的影响 范围, 可根据在彭庄4 3 0 2工作面定期水平以及 水准测得的数据, 计算分析得到倾斜的最大值和 下沉的最大值, 利用已知公式求得主要影响半 径。已知由最大倾斜值和主要影响半径求最大 下沉值的公式为 犠0=犻0狉 ( 7) 因此, 主要影响半径的计算公式为 狉=犠0/犻0 ( 8) 式中 狉为主要影响半径;犠0为最大下沉值;犻0为 最大倾斜值。 由观测数据经过软件处理获得的下沉曲线 得到的 结 果 知,4 3 0 2工 作 面 的 最 大 倾 斜 值 为 2 . 0 1mm/m, 带入公式进行计算, 影响半径为 狉=狑0/犻0=4 1 4/2. 0 1=2 0 6犿 ( 9) 相应主要影响角正切为 T a n β=犎/ 狉=5 3 0/2 0 6=2. 5 7 ( 1 0) 结合测量成果及4条观测线的下沉曲线图, 根据工作面的实际开采情况以及经验公式可得到 岩层的移动预计参数为 彭庄4 3 0 2工作面的地表 下沉系数狇= 0 . 1 9, 主要影响角正切t a n β= 2 . 5 7, 水 平移动系数犫= 0 . 3 0, 主要影响半径狉= 2 0 6m。 将分析研究的4 3 0 2工作面作为计算对象, 利用开采沉陷数据处理与预计软件S O D P 3分析 最优参数, 从而获取地表移动变形参数。进一步 得到概率积分法最优参数, 实质上就是求取使误 差函数值最小的参数值[ 9]。概率积分法预计模 型的误差函数是典型的多峰多态函数, 分析预计 模型最优参数的关键在于寻找全局最优参数。 在软件中建好工程并将工作面数据导入软件。 使用直线段, 标注出该工作面走向方向和倾 向方向, 然后使用多段线逆时针绘制出开采区 域。拾取煤层底板等高线上的点或巷道导线点, 输入各个点的采深值。采深点拾取完毕后, 将拾 取的采深点导入工程。将拾取的工程数据导入 求参项目中。输入实测数据的观测日期、 实测数 据的类型, 并将实测数据复制到文本框中, 完成 求参数据的导入。点击“ 开始解算” , 等待最优参 数解算完毕[ 1 0], 求参结果如图4所示。 预计的参数为 下沉系数狇=0. 1 8, 主要影响 角正切t a n β=2. 4 7, 拐点偏移距犽=0. 0 7, 预计中 误差为2 0 8, 残差平方和为17 3 45 9 3。 将预计分析得到的地表变形参数与实测数 据处理计算获得的地表变形参数进行对比, 两种 方法获得的下沉系数与主要影响角正切基本相 同, 对传统方法获得的参数进一步验证[ 1 1]。 3 采煤沉陷预计软件分析 在彭庄煤矿该条采工况所对应的开采条件, 8601 第3 4卷 第8期刘承旭, 崔志浩, 李增鑫.彭庄矿区村庄下条带开采地表的移动规律 运用基于S u f f e r平台的采煤地表移动变形预计 软件对工作面的地表移动进行预计。 表1 4 3 0 2工作面地表变形预计参数表 变形预计参数数值变形预计参数数值 角点个数/个 4 煤层倾角/ ( )1 2. 4 主要影响角正切 2. 5 7 开采影响传播角/ ( )9 0 下沉系数 0. 1 9 水平移动系数 0. 3 0 煤层采厚/m 3. 0 7 平均采深/m 5 5 0 参照该工作面经过观测处理获得如表1所 示的相关地表变形预计参数, 利用地表变形预计 软件对地表各种移动变形进行预计, 所得部分结 果如图4、 图5所示。 图4 该条采条件下地表下沉等值线分布图 图5 该条采条件下地表水平变形等值线分布图 由图5可以看出条采两个工作面地表下沉 值是在2 1 0mm~2 7 0mm之间, 沿着两个工作面 的布置走向呈现增大趋势。两个工作面之间受 围岩应力平衡被破坏而发生的地表移动最大, 最 大下沉值达到了2 7 0mm。由图6可以看出两个 工作面中间地表变形最大, 最大水平变形值达到 了1. 3mm/m。对比两张图可以看出水平变形范 围略大于地表下沉范围, 其原因主要是由于在覆 岩关键层的控制下坚硬厚岩层的挠曲变形使得 地表水平变形范围较下沉范围大。 表2 4 3 0 2工作面各移动变形值统计表 移动变性参数移动变形值 最大下沉值犠/mm 2 7 0 最大倾斜值犻/ (mm/m) 1. 1 最大水平移动μ/mm 7 0 最大曲率值犓/m 0. 0 0 31 0 -3 最大水平变形值ε/ (mm/m) 1. 3 水平变形值范围/ (mm/m) 0~1. 3 根据 建筑物、 水体、 铁路及主要井巷煤柱留 设与压煤开采规范 规程的评价标准, 预计所得 表2中的各项移动变形值均满足该规程建筑物 Ⅰ级损害的要求, 对地表建筑物基本无破坏 [1 2]。 4 结束语 ( 1) 通过对彭庄煤矿4 3 0 2两个条采工作面 地表移动观测站所取得的数据进行处理分析, 探 求了硬岩条件下条带开采前后地表移动变形的 特点和规律, 取得了地表沉降的些移动和变形参 数, 对类似条件矿区建下采煤及保安煤柱留设具 有一定的借鉴价值。 ( 2) 基岩在煤矿采动及上覆厚褐土层载荷影 响的综合作用下, 导致预计所得的地表水平变形 范围大于沉陷范围, 同时在关键层厚硬砂岩的控 制下, 地表移动下沉盆地未出现异常, 协调平缓, 得到的各项地表移动变形参数均在I级破坏范围 内, 可保证开采安全。 ( 3) 通过实地测量数据的计算分析结合地表 移动变形预计所得各项参数的论证分析可得 在 彭庄矿区地质采矿条件下, 运用条带开采方法进 行开采, 在保护地表建筑物的同时提高了经济效 益, 同时验证了原开采方案设计参数的选取是合 适的。 参考文献 [1]邓喀中.变形监测及沉陷工程学[M].江苏徐州 中国矿业大 学出版社,2 0 1 4. 9601 北京测绘第3 4卷 第8期 [2]余学义, 郭文彬, 赵兵朝, 等.巨厚黄土层宽条带开采地表移 动规律及参数优化[J].煤炭科学技术,2 0 1 6,4 4(4) 6 1 0. [3]刘义新, 戴华阳, 郭文兵.巨厚松散层下深部宽条带开采地表 移动规律[J].采矿与安全工程学报,2 0 0 9,2 6(3) 3 3 6 3 4 0. [4]余学义, 刘传杰.亭南煤矿村庄下宽条带开采地表移动规律 研究[J].煤炭工程,2 0 1 9,5 1(5) 1 1 8 1 2 2. [5]王钰.基于D I n S A R及数值模拟的张双楼矿开采沉陷规律 研究[D].江苏徐州 中国矿业大学,2 0 1 9. [6]余学义, 刘传杰.亭南煤矿村庄下宽条带开采地表移动规律 研究[J].煤炭工程,2 0 1 9,5 1(5) 1 1 8 1 2 2. [7]李家琪.菏泽赵楼矿深部开采地表沉陷变形规律与预计参数 研究[D].辽宁阜新 辽宁工程技术大学,2 0 1 7. [8]刘娜, 栾元重, 刘阳.欢城煤矿二号井村庄下压煤开采设计分 析[J].煤矿开采,2 0 1 1,1 6(1) 4 4 4 7. [9]高永格, 牛矗, 张强, 等.厚松散层下采煤地表沉陷特征研究 [J].煤炭科学技术,2 0 1 9,4 7(6) 1 9 2 1 9 8. [1 0]邹友峰, 陈俊杰, 张子月.矿山开采沉陷类项目基金资助现 状与学科发展方向展望 基于国家自然科学基金资助的 视角[J].河南理工大学学报自然科学版,2 0 1 7,3 6(6) 1 7. [1 1]刘欣, 吴昊, 贾勇帅, 等.基于概率积分法的矿区地表沉陷预 计分析[J].北京测绘,2 0 1 7(1) 1 2 3 1 2 6. [1 2]田华, 李淑梅, 张利群.大采深硬岩层条件下地表移动变形 规律研究[J].北京测绘,2 0 1 6(3) 1 1 1 4. 犛 狌 狉 犳 犪 犮 犲犕 狅 狏 犲 犿 犲 狀 狋犔 犪 狑狅 犳犛 狋 狉 犻 狆犕 犻 狀 犻 狀 犵 犻 狀犞 犻 犾 犾 犪 犵 犲 狊 犻 狀犘 犲 狀 犵 狕 犺 狌 犪 狀 犵犕 犻 狀 犻 狀 犵犃 狉 犲 犪 犔 犐 犝犆 犺 犲 狀犵狓 狌,犆 犝 犐犣 犺 犻 犺 犪 狅,犔 犐犣 犲 狀犵狓 犻 狀 (犆 狅 犾 犾 犲 犵犲狅 犳犌 犲 狅 犿 犪 狋 犻 犮 狊,犛 犺 犪 狀 犱 狅 狀犵犝 狀 犻 狏 犲 狉 狊 犻 狋狔狅 犳犛 犮 犻 犲 狀 犮 犲犪 狀 犱犜 犲 犮 犺 狀 狅 犾 狅犵 狔,犙 犻 狀犵犱 犪 狅犛 犺 犪 狀 犱 狅 狀犵2 6 6 5 9 0,犆 犺 犻 狀 犪) 犃 犫 狊 狋 狉 犪 犮 狋犐 狀狅 狉 犱 犲 狉 狋 狅狊 狋 狌 犱狔狋 犺 犲狊 狌 狉 犳 犪 犮 犲犿 狅 狏 犲 犿 犲 狀 狋 犾 犪 狑狅 犳狋 犺 犲狊 狋 狉 犻狆犿 犻 狀 犻 狀犵狑 狅 狉 犽 犻 狀犵犳 犪 犮 犲,狋 犺 犻 狊犪 狉 狋 犻 犮 犾 犲狑 犪 狊犫 犪 狊 犲 犱狅 狀狋 犺 犲 犲 狓 犻 狊 狋 犻 狀犵犱 犪 狋 犪狅 犳 狋 犺 犲4 3 0 2狑 狅 狉 犽 犻 狀犵犳 犪 犮 犲犪 狀 犱狋 犺 犲犪 犮 狋 狌 犪 犾 犾 犲 狏 犲 犾犪 狀 犱犺 狅 狉 犻 狕 狅 狀 狋 犪 犾狅 犫 狊 犲 狉 狏 犪 狋 犻 狅 狀犱 犪 狋 犪狅 犳狋 犺 犲狑 狅 狉 犽 犻 狀犵犳 犪 犮 犲 .犜 犺 犲 犵狉 狅 狌 狀 犱狊 狌 犫 狊 犻 犱 犲 狀 犮 犲 犾 犪 狑狑 犪 狊 狉 犲 狊 犲 犪 狉 犮 犺 犲 犱犱 狌 狉 犻 狀犵狋 犺 犲犪 犱 狏 犪 狀 犮 犲 犿 犲 狀 狋狅 犳犿 犻 狀 犻 狀犵犮 狅 狀 犱 犻 狋 犻 狅 狀 狊 . 犅 犪 狊 犲 犱狅 狀狋 犺 犲犲 狓 犻 狊 狋 犻 狀犵狋 犺 犲 狅 狉 犲 狋 犻 犮 犪 犾 犫 犪 狊 犻 狊,狋 犺 犲犛 犺 犪 狀犵犺 犪 犻 犻 狀 犱 犲 狓犿 狅 狏 犲 犿 犲 狀 狋 犪 狀犵犾 犲狑 犪 狊犿 犪狆 狆犲 犱犪 狀 犱狋 犺 犲犵狉 狅 狌 狀 犱狊 狌 狉 犳 犪 犮 犲犿 狅 狏 犲 犿 犲 狀 狋狆犪 狉 犪 犿 犲 狋 犲 狉 狊狑 犲 狉 犲犮 犪 犾 犮 狌 犾 犪 狋 犲 犱 . 犅 犪 狊 犲 犱狅 狀 狋 犺 犲 狋 犺 犲 狅 狉 犲 狋 犻 犮 犪 犾 犪 狀 犪 犾狔狊 犻 狊 犮 狅 犿 犫 犻 狀 犲 犱狑 犻 狋 犺 狋 犺 犲犪 犮 狋 狌 犪 犾 狅 犫 狊 犲 狉 狏 犪 狋 犻 狅 狀 狊,狋 犺 犲狆狉 狅 犫 犪 犫 犻 犾 犻 狋狔犻 狀 狋 犲犵狉 犪 狋 犻 狅 狀犿 犲 狋 犺 狅 犱狑 犪 狊 犪狆 狆犾 犻 犲 犱 狋 狅狋 犺 犲犿 犻 狀 犻 狀犵狊 狌 狉 犳 犪 犮 犲犱 犲 犳 狅 狉 犿 犪 狋 犻 狅 狀狑 犺 犻 犮 犺狑 犪 狊 犲 狓狆犲 犮 狋 犲 犱 . 犜 犺 犲犲 狓狆犲 犮 狋 犲 犱狉 犲 狊 狌 犾 狋 狊 狊 犺 狅 狑狋 犺 犪 狋 狊 狋 狉 犻狆犿 犻 狀 犻 狀犵犺 犪 犱犪犫 犲 狋 狋 犲 狉犲 犳 犳 犲 犮 狋 狅 狀犮 狅 狀 狋 狉 狅 犾 犾 犻 狀犵狊 狌 狉 犳 犪 犮 犲犿 狅 狏 犲 犿 犲 狀 狋犪 狀 犱犱 犲 犳 狅 狉 犿 犪 狋 犻 狅 狀 . 犓 犲 狔狑 狅 狉 犱 狊狊 狌 狉 犳 犪 犮 犲犿 狅 狏 犲 犿 犲 狀 狋;狊 狋 狉 犻狆犿 犻 狀 犻 狀犵;狊 狌 犫 狊 犻 犱 犲 狀 犮 犲 犾 犪 狑;犱 犲 犳 狅 狉 犿 犪 狋 犻 狅 狀狆狉 犲 犱 犻 犮 狋 犻 狅 狀 0701