密实充填采煤的岩层移动理论研究.pdf
书书书 第4 1卷 第6期 中国矿业大学学报 V o l . 4 1 N o . 6 2 0 1 2年1 1月 J o u r n a l o f C h i n a U n i v e r s i t y o f M i n i n g & T e c h n o l o g y N o v . 2 0 1 2 收稿日期2 0 1 2-0 5-2 8 基金项目高等学校学科创新引智计划项目(B 0 7 0 2 8) ; 国家自然科学基金重点项目( 5 0 8 3 4 0 0 4) 作者简介缪协兴(1 9 5 9-) , 男, 江苏省江阴市人, 教授, 博士生导师, 工学博士,从事采矿工程与技术方面的研究. E-m a i lx x m i a o @c u m t . e d u. c n T e l0 5 1 6-8 3 9 9 5 5 1 8 密实充填采煤的岩层移动理论研究 缪协兴1,黄艳利2,巨 峰1,茅献彪1,郭广礼3,张吉雄2 (1.中国矿业大学 深部岩土力学与地下工程国家重点实验室, 江苏 徐州 2 2 1 0 0 8; 2.中国矿业大学 矿业工程学院,江苏 徐州 2 2 1 1 1 6;3.中国矿业大学 环境与测绘学院,江苏 徐州 2 2 1 1 1 6) 摘要综合机械化固体充填采煤技术能实现密实充填, 其岩层移动和地表沉陷规律与垮落式采 煤显著不同. 在较系统分析充填采煤岩层移动的等价采高理论基础上, 建立了密实充填采煤岩层 移动的结构力学模型, 并进行了连续介质力学分析, 给出了相应的岩层移动和地表沉陷计算公 式.通过实例分析, 提出了在密实充填后, 如还采用概率积分法预计地表沉陷, 则应对传统的地表 下沉预计参数加以修正.可为密实充填采煤覆岩移动和地表沉陷分析提供相关思路以及理论和 方法. 关键词密实充填采煤;岩层移动与地表沉陷;连续介质力学;地表沉陷预计 中图分类号T D 8 2 3文献标识码A 文章编号 1 0 0 0-1 9 6 4(2 0 1 2)0 6-0 8 6 3-0 5 S t r a t a m o v e m e n t t h e o r y o f d e n s e b a c k f i l l m i n i n g M I AO X i e-x i n g 1,HUANG Y a n -l i 2, J U F e n g 1, MAO X i a n-b i a o 1, GUO G u a n g-l i 3, Z HANG J i-x i o n g 2 (1. S t a t e K e y L a b o r a t o r y f o r G e o m e c h a n i c s & D e e p U n d e r g r o u n d E n g i n e e r i n g,C h i n a U n i v e r s i t y o f M i n i n g & T e c h n o l o g y,X u z h o u ,J i a n g s u 2 2 1 0 0 8,C h i n a;2. S c h o o l o f M i n e s,C h i n a U n i v e r s i t y o f M i n i n g & T e c h n o l o g y, X u z h o u,J i a n g s u 2 2 1 1 1 6,C h i n a;3. S c h o o l o f E n v i r o n m e n t S c i e n c e a n d S p a t i a l I n f o r m a t i c s, C h i n a U n i v e r s i t y o f M i n i n g & T e c h n o l o g y, X u z h o u,J i a n g s u 2 2 1 1 1 6,C h i n a) A b s t r a c tS i n c e t h e f u l l y m e c h a n i z e d s o l i d b a c k f i l l a n d m i n i n g t e c h n o l o g y a c t u a l i z e s d e n s e b a c k- f i l l i n g o a f,t h e l a w s o f s t r a t a m o v e m e n t a n d s u r f a c e s u b s i d e n c e i n d e n s e b a c k f i l l m i n i n g m e t h o d a r e d i f f e r e n t f r o m c a v i n g m i n i n g m e t h o d . B a s e d o n t h e s y s t e m s a n a l y s i s f o r e q u i v a l e n t m i n i n g h e i g h t t h e o r y i n b a c k f i l l m i n i n g s t r a t a m o v e m e n t,t h i s p a p e r e s t a b l i s h e s t h e s t r u c t u r a l m e c h a n- i c s m o d e l,a n a l y s e s m e c h a n i c a l p r o p e r t i e s o f c o n t i n u o u s m e d i a,a n d d e d u c e s c o r r e s p o n d i n g f o r- m u l a s o f s t r a t a m o v e m e n t a n d s u r f a c e s u b s i d e n c e . I n a d d i t i o n,t h r o u g h a c a s e a n a l y s i s,t h e p a- p e r p u t s f o r w a r d t h a t t h e t r a d i t i o n a l p r e d i c t i o n p a r a m e t e r s s h o u l d b e a m e n d e d w h e n t h e p r o b a- b i l i t y i n t e g r a l m e t h o d w a s u s e d t o p r e d i c t t h e s u r f a c e s u b s i d e n c e . T h e r e s e a r c h r e s u l t s p r o v i d e i d e a s,t h e o r i e s a n d m e t h o d s f o r s t r a t a m o v e m e n t a n d s u r f a c e s u b s i d e n c e a n a l y s i s i n d e n s e b a c k- f i l l m i n i n g . K e y w o r d sd e n s e b a c k f i l l m i n i n g;s t r a t a m o v e m e n t a n d s u r f a c e s u b s i d e n c e;m e c h a n i c s o f c o n- t i n u o u s m e d i a;s u r f a c e s u b s i d e n c e p r e d i c t i o n 随着综合机械化固体充填采煤技术在我国开 发成功并进入规模化推广应用阶段[ 1-2], 使近百年 来的充填采煤技术[ 3-4]攻关取得了重大突破, 为实 现煤矿安全、 绿色和高效的科学开采提供了一种创 新技术[ 5]. 采用综合机械化固体充填采煤技术能实 现密实充填, 而密实充填后的采动覆岩移动和地表 中 国 矿 业 大 学 学 报 第4 1卷 沉陷规律已显著不同于垮落式开采, 必须深入研究 密实充填采煤的采动覆岩移动和地表沉陷机理, 建 立相应的分析方法, 并给出其计算公式. 本文基于 充填采煤岩层移动分析的等价采高理论, 建立了密 实充填采煤岩层移动的连续介质力学模型, 并进行 了相关力学分析, 给出了相应的计算公式. 1 充填采煤岩层移动分析的等价采高理论 实施采空区充填后, 从岩层移动分析的角度 看, 就是使实际采出的煤层高度降低了. 所谓等价 采高, 是指实际采高减去充入物体的高( 厚) 度. 图 1为充填采煤的等价采高分析模型, 图中M为实 际采高,Me为M减去充入物体厚度后的高度, 即 等价采高. 我们知道, 实际充入物体经采场覆岩下 沉运动后总会有一定的压缩变形量. 以固体充填为 例, 任何固体在破碎状态都有一定的体积松散膨胀 量, 而在不同的压力下都有相应的体积压实收缩 量.假设充填固体经最大( 最终) 覆岩下沉压实后的 高度为M0, 则, 根据定义, 固体充填采煤后岩层移 动分析的等价采高为 Me=M-M0. ( 1) 图1 充填采煤的等价采高分析模型 F i g . 1A n a l y s i s m o d e l o f e q u i v a l e n t m i n i n g h e i g h t i n b a c k f i l l m i n i n g 假设 充填时的采空区已有初始下沉量为Mc; 初始充填体的压实系数为k c; 最终充填体的压实 系数为k 0 . 则 M0为 M0=M-Mc-(kc-k0) (M-Mc)= ( 1+k0-kc) (M-Mc) , ( 2) 式( 2) 代入式(1 ) 得 Me=M-(1+k0-kc) (M-Mc), ( 3) 如果初始充填时已经充分压实, 即k c和k0十分接 近, 则从式( 3) 中可看到 Me=Mc, ( 4) 如果充填时采空区没有初始下沉, 即Mc=0, 则从 式( 3) 中可看到 Me=(kc-k0)M. ( 5) 几点说明 1) 上面提到的kc和k0可以由实验测定, 即在 一定压实力作用下松散固体的残余碎胀系统, 大小 与不同松散体的组合和他们的颗粒级配等有关. 2)如果充填时的顶板下沉量很小, 即Mc很 小; 同时初始压实力( 机械夯实力) 又很大, 即k c和 k0很接近. 这样Me很小, 我们把这种情况称为密 实充填. 3)这里, 把M0与M的比值定义为充实率η, 即 η= M0 M 1 0 0%.(6) 如果η大于8 5%, 则可视为密实充填采煤.目 前, 用我国开发出的综合机械化固体充填采煤技术 就能实现密实充填. 4)所谓等价采高, 在数学上是一个极限的概 念, 是指在覆岩可能形成的最大下沉压力下的最小 k0值, 在实际工程中是岩层移动和地表沉陷分析 的上限值, 其分析结果偏于安全. 2 密实充填采煤岩层移动的连续介质力学 分析 根据大量的现场实测和物理模拟等研究表明, 综合机械化固体密实充填采煤技术可控制到直接 顶仅发育微小裂隙, 老顶岩层仅有弯曲下沉变形而 不发生破断.因此, 可用连续介质力学模型来对老 顶岩层进行岩层移动分析.密实充填采煤岩层移动 的结构模型见图2. 图2 密实充填采煤岩层移动的结构模型 F i g . 2S t r u c t u r e m o d e l o f s t r a t a m o v e m e n t i n d e n s e b a c k f i l l m i n i n g 从图2可以看到, 在充填体和直接顶的支撑作 468 第6期 缪协兴等密实充填采煤的岩层移动理论研究 用下, 老顶岩层仅发生弯曲下沉变形而不发生破 断, 因此, 我们可以将其视为 W i n k l e r弹性地基上 的半无限长岩梁进行力学分析, 其受力简图见图 3. 图3 密实充填采煤老顶岩梁的受力简图 F i g . 3S k e t c h o f m a i n r o o f b e a m i n d e n s e b a c k f i l l m i n i n g 根据 W i n k l e r地基假设[ 6] 地基表面任一点的 沉降量与该点单位面积上所受的压力成正比, 即 p1(x)=k1w(x) (0≤x≤l), ( 7) p2(x)=k2w(x) (l≤x≤ ∞) , ( 8) 式中p 1(x) ,p2(x)分别为老顶岩梁受到的地基 支撑力;w( x) 为老顶岩梁的弯曲下沉挠度;k1,k2 为地基系数,k 1由充填体的变形特性确定,k2由老 顶岩梁下直接顶和煤层等抗变形性能确定. 老顶岩梁的挠曲微分方程为 E Id 4 w1(x) dx 4 +k1w1(x)=q0 ( 0≤ x≤l) ,(9) E Id 4 w2(x) dx 4 +k2w2(x)=q0 (l≤x≤ ∞). ( 1 0) 令 β 1= 4 k1 4槡E I, β 2= 4 k2 4槡E I. 则式( 9) , (1 0 ) 为 d 4 w1(x) dx 4 +4 β 4 1w1(x)= q0 E I ( 0≤ x≤l) ,(1 1) d 4 w2(x) dx 4 +4 β 4 2w2(x)= q0 E I ( l≤x≤ ∞).(1 2) 根据初参数法对式(1 1) 和式(1 2) 求解, 得到 w1(x)=w0φ1( β 1x)+θ0 1 β 1φ 2( β 1x)- M0 1 E I β 2 1 φ3(β 1x)-Q0 1 E I β 3 1 φ4(β 1x)+ q0 k1[ 1-φ1( β 1x) ] ( 0≤ x≤l) , ( 1 3) w2(x-l)=wlφ1[ β 2(x-l) ]+ θl 1 β 2φ 2[ β 2(x-l) ]-Ml 1 E I β 2 2 φ3[β 2(x-l) ]- Ql 1 E I β 3 2 φ4[β 2(x-l) ]-q 0 k2φ 1[ β 2(w-l) ] , ( l≤x≤ ∞) , ( 1 4) 式中w0, θ0,M0和Q0分别为老顶岩梁端点(x= 0)的挠度、 转角、 弯矩和剪力;wl,θl,Ml和Ql分别 为老顶岩梁x=l位置的挠度、 转角、 弯矩和剪力. φ1(β x)=c hβxc o sβx, φ2(β x)= 1 2( c hβxs i nβx+s hβxc o sβx) , φ3(β x)= 1 2 s hβxs i nβx, φ4(β x)= 1 4( c hβxs i nβx-s hβxc o sβx) 烅 烄 烆 . ( 1 5) 端部(x=0) 边界条件 Q1(x)x=0=0, θ1(x)x=0=0 { . ( 1 6) x=l处的连续性条件 w1(x)x=l=w2(x-l)x=l=wl, θ1(x)x=l=θ2(x-l)x=l=θl, M1(x)x=l=M2(x-l)x=l=Ml, Q1(x)x=l=Q2(x-l)x=l=Ql 烅 烄 烆 . ( 1 7) 无穷远处边界条件 w2(x-l)x→∞=0, Q2(x-l)x→∞=0. ( 1 8) 由式(1 6) 可求得Q0=0,θ0=0. 由式( 1 7) , (1 8) 可求得 w0φ1( β 1l)-M0 1 E I β 2 1 φ3(β 1l)+ q 0 k1[ 1-φ1( β 1l) ]=wl, -w04 β 1φ4( β 1l)-M0 1 E I β1 φ2(β 1l)+ q04 β 1 k1φ 4( β 1l)=θl, w04E I β 2 1φ3( β 1l)+M0φ1( β 1l)- q0 1 β 2 1φ 3( β 1l)=Ml, w04E I β 3 1φ2( β 1l)-M04β1φ4( β 1l)- q0 1 β 1φ 2( β 1l)=Ql, wlφ1(∞)+θl 1 β 2φ 2(x)-Ml 1 E I β 2 2 φ3( x)- Ql 1 E I β 3 2 φ4( x)-q 0 k2φ 1(x)|x→∞=0, wl4E I β 2 2φ3(x)+θl4E Iβ2φ4(x)+Mlφ1(x)+ Ql 1 β 2φ 2(x)-q 0 β 2 2φ 3(x)|x→∞=0 烅 烄 烆 . ( 1 9) 求解方程组(1 9) , 可得w0,M0,wl,θl,Ml,Ql 参数的解, 并将其代入式( 1 3) , (1 4) , 即可得老顶岩 梁的挠度w(x)表达式. 由于w0,M0,wl, θl,Ml,Ql 的表达式较复杂, 这里不予以具体给出. 根据某试验区固体充填采煤的开采条件和充 填材料的压实特性测试试验, 可得相关计算参数 纵向长度l=1 0 0m, 老顶岩梁高度h=1 0m, 弹性 模量E=5. 0G P a, 岩梁上作用的载荷q 0=7. 0 568 中 国 矿 业 大 学 学 报 第4 1卷 MP a; 煤层的弹性模量Em=3. 0G P a, 煤层厚度为 3m, 则地基系数可近似取为k2=1. 0GN/m 3 ; 如 充填体的压缩弹性模量分别为E ′=3 0 0,2 0 0,1 0 0 MP a, 则其地基系数k1可分别取为0. 1GN/m 3, 6 6. 7MN/m 3, 3 3. 3MN/m 3 . 由式( 1 3)~(1 9)可 计算得到老顶岩梁的下沉曲线见图4. 图4 老顶岩梁的下沉挠度曲线 F i g . 4D e f l e c t i o n c u r v e o f m a i n r o o f 3 密实充填采煤岩层移动计算公式 从图2密实充填采煤岩层移动的结构模型中 可以看到, 老顶岩层上方覆岩中的任何岩层之间都 不出现离层现象, 仅发生连续微小的弯曲下沉, 且 越向上弯曲挠度越小.因此, 可以假设, 在覆岩下沉 影响角θ之内, 老顶岩梁上方任意一个水平面( y水 平面)的下沉体积总量与老顶岩梁的下沉体积总 量相等,即 ∫ l+hc o tθ 0 w(x)dx= ∫ l+yc o tθ 0 wy(x)dx, ( 2 0) 式中wy( x)为y水平面上岩层的下沉挠度;h为 煤层底板到老顶岩梁顶部的高度. 由于密实充填后老顶岩梁上各岩层的弯曲挠 度很小, 因此, 同一弯曲挠度的曲率变化也很小, 其 挠度曲线可以近似为直线, 如图2所示,也就是可 以把三边形O AD和三边形B C F近似为△O AD和 △B C F,则由式(2 0) 可得 wy(x)=l+y c o tθ-x l+yc o tθ wy(0), ( 2 1) 其中 wy(0)=l+h c o tθ l+yc o tθ w(0). ( 2 2) 4 传统岩层移动与地表沉陷预计的修正 概率积分法是传统岩层移动与地表沉陷预计 的代表性方法[ 7], 其地表下沉系数q是最重要的地 表下沉预计参数, 对近水平煤层在充分采动下q的 表达式为 q=wm a x/M ,( 2 3) 式中 wm a x为在充分采动下地表最大下沉值. 众所周知, 在传统垮落式开采中覆岩运动会形 成“ 三带” , 即垮落带、 破碎带和裂隙带, 因而相对密 实充填开采, 即使垮落式开采采高与等价采高相 等, 传递到地表的下沉变形量相对会小. 这里, 我们 用密实充填开采实测证实的式( 2 3) 计算公式与传 统概率积分法预计的结果加以比较. 目前, 综合机械化固体密实充填采煤技术已经 在我国多个矿区得到了推广应用, 这里, 以济宁花 园矿、 兖州济三矿和平煤十二矿为例. 这3个煤矿 的采矿地质条件见表1. 表1 3个煤矿的实际采矿地质条件情况 T a b l e 1 G e o l o g i c m i n i n g c o n d i t i o n o f t h r e e c o a l m i n e s 充填采煤区 采深/ m 采高/ m 老顶厚/ m 充实率/% 花园煤矿 5 7 5 2 . 5 5 1 8 9 济三煤矿6 8 3 3 . 8 4 2 8 5 十二煤矿3 5 5 3 . 3 3 1 8 8 根据3个煤矿的实际采矿地质条件, 我们做了 2种模拟计算分析.一种是改变采高, 即等价采高 和垮落式开采采高同时改变; 另一种是在采高为 0 . 5m的情况下改变采深, 其两种模拟计算的结果 分别见表2和表3.值得说明的是, 在模拟参数与 实际参数相同情况下的计算结果已得到实践验证. 表2 模拟不同采高条件的下沉系数 T a b l e 2 S i m u l a t i o n o f s u b s i d e n c e f a c t o r v a l u e i n d i f f e r e n t m i n i n g h e i g h t 采高/ m 济宁花园矿 垮落法 密实充填 兖州济三矿 垮落法 密实充填 平煤十二矿 垮落法 密实充填 0. 5 0. 9 2 0. 9 3 0. 8 2 0. 8 3 0. 7 3 0. 7 5 1. 0 0. 9 2 0. 9 3 0. 8 1 0. 8 2 0. 7 2 0. 7 4 1. 5 0. 9 1 0. 9 2 0. 8 0 0. 8 2 0. 7 0 0. 7 3 表3 模拟不同采深条件的下沉系数 T a b l e 3 S i m u l a t i o n o f s u b s i d e n c e f a c t o r v a l u e i n d i f f e r e n t m i n i n g d e p t h 采深/ m 济宁花园矿 垮落法 密实充填 兖州济三矿 垮落法 密实充填 平煤十二矿 垮落法 密实充填 3 0 0 0. 9 5 0. 9 6 0. 8 5 0. 8 6 0. 7 3 0. 7 6 4 0 0 0. 9 4 0. 9 5 0. 8 4 0. 8 5 0. 7 0 0. 7 5 5 0 0 0. 9 3 0. 9 5 0. 8 3 0. 8 5 0. 6 8 0. 7 4 6 0 0 0. 9 2 0. 9 4 0. 8 2 0. 8 4 0. 6 6 0. 7 2 7 0 0 0. 9 0 0. 9 4 0. 8 1 0. 8 3 0. 6 4 0. 7 0 从表2和表3中可以看到, 密实充填开采后, 采场覆岩主要为整体下沉, 因而下沉系数q值一般 大于垮落式开采. 5 结论 实施密实充填采煤后, 可控制到采动覆岩中老 顶岩层不发生破断, 并且其上方岩层中也不发生离 层, 仅为整体连续下沉变形, 因而可在连续介质力 学范畴内研究采动覆岩移动和地表沉陷规律. 我们 668 第6期 缪协兴等密实充填采煤的岩层移动理论研究 把老顶岩梁简化为 W i n k l e r弹性地基梁, 并给出了 挠度的具体表达式. 由于老顶岩层上方的岩层移动 量很小, 基于每一岩层水平上的下沉体积相等和下 沉盆地平面形态近似三角形假设, 得到了任意一个 水平面上的下沉挠曲线计算公式.另外, 如果要采 用概率积分法对密实充填采煤的地表沉陷进行预 计, 则对传统地表下沉预计的参数需加以修正. 参考文献 [1] 缪协兴,张吉雄,郭广礼.综合机械化固体废弃物充 填采煤方法与技术[M].徐州 中国矿业大学出版 社,2 0 1 0. [2] 缪协兴,张吉雄,郭广礼.综合机械化固体充填采煤 方法与技术研究[J].煤炭学报, 2 0 1 0,3 5(1) 1-6. M I AO X i e-x i n g,Z HANG J i-x i o n g,GUO G u a n g-l i . S t u d y o n w a s t e f i l l i n g m e t h o d a n d t e c h n o l o g y i n f u l l y m e c h a n i z e d c o a l m i n i n g[J]. J o u r n a l o f C h i n a C o a l S o- c i e t y,2 0 1 0,3 5( 1) 1-6. [3] 刘同有,蔡嗣经.国内外膏体充填技术的应用与研 究现状[J].中国矿业, 1 9 9 8,7(5) 1-4. L I U T o n g- y o u,C A I S i -j i n g .S t a t u s q u o o f a p p l i c a- t i o n a n d r e s e a r c h o f p a s t e f i l l t e c h n o l o g y i n C h i n a a n d a b r o a d[J]. C h i n a M i n i n g M a g a z i n e,1 9 9 8,7(5) 1- 4. [4] 周华强,侯朝炯,孙希奎,等.固体废物膏体充填不 迁村采煤 [J].中 国 矿 业 大 学 学 报,2 0 0 4,3 3( 2) 1 5 4-1 5 8. Z HOU H u a-q i a n g,HOU C h a o-j i o n g,S UN X i-k u i, e t a l . S o l i d w a s t e p a s t e f i l l i n g f o r n o n e-v i l l a g e-r e l o c a- t i o n c o a l m i n i n g[J]. J o u r n a l o f C h i n a U n i v e r s i t y o f M i n i n g & T e c h n o l o g y,2 0 0 4, 3 3(2) 1 5 4-1 5 8. [5] 缪协兴,钱鸣高.中国煤炭资源绿色开采研究现状 与展望[J].采矿与安全工程学报, 2 0 0 9,2 6(1) 1- 1 4. M I AO X i e-x i n g,Q I AN M i n g-g a o . R e s e a r c h o n g r e e n m i n i n g o f c o a l r e s o u r c e s i n C h i n a c u r r e n t s t a t u s a n d f u t u r e p r o s p e c t s[J]. J o u r n a l o f M i n i n g &S a f e t y E n- g i n e e r i n g,2 0 0 9,2 6(1) 1-1 4. [6] 龙驭球.弹性地基梁的计算[M].北京 人民教育出 版社,1 9 8 2. [7] 邹友峰,邓喀中,马伟民.矿山开采沉陷工程[M]. 徐州 中国矿业大学出版社,2 0 0 3. ( 责任编辑骆振福) 768