基于UDEC的岩层与地表移动动态模拟研究_刘吉波

1、第 19 卷 第 3 期 ( 总第 118 期)煤 矿 开 采Vol. 19No. 3 ( Series No. 118)2014 年 6 月COAL MINING TECHNOLOGYJune 2014基于 UDEC 的岩层与地表移动动态模拟研究刘吉波，廉旭刚，戴华阳，杨国柱，刘杰( 中国矿业大学 ( 北京)地球科学与测绘工程学院，北京 100083)摘 要 以 UDEC 为研究平台，讨论了数值计算时的开挖模拟和迭代运算问题。通过计算实验，提出并建立了以开采充分性和推进速度为变量的单元开采迭代次数的计算方法，分析了岩层性质及工作面推进速度等因素对岩层及地表移动变形的影响规律，为岩层及地表动态

2、移动变形的模拟计算提供了参考。关键词UDEC; 岩层与地表移动; 数值模拟; 迭代次数; 开挖模拟中图分类号 TD325文献标识码 A文章编号 1006-6225 ( 2014)03-0104-04Dynamical Simulation of ock Strata and Surface Movement Based on UDECLIU Ji-bo，LIAN Xu-gang，DAI Hua-yang，YANG Guo-zhu，LIU Jie( Earth Science Surveying Engineering School，China University of Mining Tech

3、nology ( Beijing) ，Beijing 100083，China)Abstract: Mining simulation and iteration calculation in numerical simulation was discussed based on UDEC software. A calculation method of zone mining iterations was put forward which took mining sufficiency and mining speed as variables. By this method，the i

4、n-fluence rule of rock properties and mining speed on rock and surface movement and deformation were analyzed，which provided refer-ence for simulation and calculation of rock and surface dynamic movement and deformation.Keywords: UDEC; rock and surface movement; numerical simulation; iteration numbe

5、r; excavation simulation数值模拟软件是研究岩层移动规律的重要手段，很多学者应用数值模拟软件在开采沉陷领域进行了研究。杨天鸿等基于岩体渗流耦合及损伤分析系统，模拟了开采过程中覆岩动态破坏过程及应力场和渗流场演化的过程，提出了岩层渗透性演化受到开采引起的应力重新分配及覆岩损伤累积过程中1结构性破坏的分布特征 。唐春安应用自主开发的岩石破坏分析系统 FPA2D ，分析了采空区上覆岩层的移动变形规律，认为 FPA2D 岩层性质、几何条件、移动变形、破坏形态及支撑结构等方面真实地模拟开采过程中岩层发生变形、离层、断裂及2地表下沉的规律 。李新强等分析了采空区上覆岩层内部移动变形规

6、律，提出了动态分析开采沉陷预计模型，揭示了四大控制岩体运动的要素: 即推进速度、覆岩的移动演变、岩体和矸石的流变压，3缩 并据此编制了专用的有限元分析软件 。谢和平院士利用 FLAC 进行了煤矿开采沉陷预计，并将预计结果和概率积分法做了对比分析，揭示了使用 FLAC 可模拟真实开采地质条件，指出时间因素与计算迭代次数有关，即计算步数与对应的时间长度成正比，计算步数越大，模型形变越大45。但目前对开采沉陷引起的动态岩层及地表移动收稿日期 20140328的数值模拟预计方法研究仍不成熟，未提出基于时间影响的具体模型或方法。为此，本文对数值模拟迭代次数与开采充分度及推进速度的关系进行研究，建立计算匀

7、速开采与非匀速开采时在不同开采充分程度下的地表动态预计模型; 分析开采推进速度、岩层性质对岩层及地表移动变形的动态影响过程。1 瞬时开采岩层及地表移动的传递速度煤矿开采是连续的推进过程，覆岩及地表移动随着开采空间和时间变化。瞬时开采是假设在某一时间，采出特定区域的全部矿体，即将开采空间看作一常量，地表及覆岩的移动变形仅随采后时间而变化，不受开采空间的影响。借助数值计算软件可模拟实际生产中不存在的开采状况，研究不同开采充分度时，覆岩及地表移动变形与开采破坏传递速度的关系。假设采动破坏的传递速度与开采的充分程度 ( 宽深比) 有关，认为地下煤层开采引起的移动变形传至地表所需的数值模拟迭代次数为传递

8、速度。为了验证假设的正确性，进行了不同采深不同充分度条件下模型开采引起地表下沉 10mm 的模拟迭代DOI 10. 13532 / j. cnki. cn113677 / td. 2014. 03. 030基金项目 博士学科点专项科研基金项目 ( 20110023110014) ; 中央高校基本科研业务费专项资金资助 ( 2009YD05，2010QD01)作者简介 刘吉波 ( 1975) ，男，山东招远人，博士研究生，从事开采沉陷、变形监测等方面的研究与应用。引用格式 刘吉波，廉旭刚，戴华阳，等 . 基于 UDEC 的岩层与地表移动动态模拟研究 J . 煤矿开采，2014，19 ( 3) :

9、 104107.104刘吉波等: 基于 UDEC 的岩层与地表移动动态模拟研究2014 年第 3 期次数统计实验。通过表 1 分析可知，当开采尺寸一定时，采深对迭代次数具有一定的影响，但与开采充分度相比，其影响较小。表 1不同采深及开采充分度地表下沉 10mm 迭代次数采深D = 2600m迭代D = 500m迭代D = 100m迭代H / mD / H次数D / H次数D / H次数7003. 716300. 718000. 1413506004. 336300. 838100. 1713005005. 205701. 007900. 2013004006. 505501. 257800.

10、2512303008. 675201. 677800. 331200注:D 为开采尺寸。开采充分度的大小与数值模拟中计算步数成反比，即开采越充分，所需计算步数越少。开采充分度在开采沉陷的数值模拟中，起到了重要的作用。分析上述瞬时开采数值模拟结果，同等程度的地表下沉，开采时步与开采充分程度有很大关系。动态预计模型的建立，需要考虑开采时步与开采充分度及推进速度的关系，从而构建时间与开采尺寸影响的开采时步综合计算公式。元，建立单元开采后数值模拟迭代次数与开采充分度、推进速度的关系。在推进过程中，开采充分程度起主要作用，因此建立以开采充分度为主，推进速度为辅的迭代次数计算公式。ni=A( 1 +v0

11、vi)( 1)L / H02v0式中，ni 为开采单元所需的迭代次数; A 为模型平衡所需迭代次数; L / H0 为推进距离与采深比值; v0 为模型无速度影响中值，m / d; vi 为开采单元的推进速度，m / d。当工作面停止后，迭代次数与时间成正比，具体计算公式为:n =At( 2)L / H0式中， n 为工作面停采后 t 时间内迭代次数;t为工作面停采后延续时间，d。2. 2 模型验证分析2. 2. 1 数值模型建立2 动态预计模型建立与验证分析2. 1 动态预计模型建立模型建立的思想: 基于瞬时开采岩层及地表移动的传递速度，以一定尺寸的开采距离为开采单应用通用离散元数值模拟软件

12、 ( UDEC) ，模型建立尺寸为 4200m750m。表 2 为模型的物理力学参数。地质采矿条件: 采深 700m，松散层厚度400m，基岩厚度 300m，采厚 3m，煤层倾角 0。表 2数值模拟物理力学参数序号岩性厚度 / m密度 /块体尺寸弹性模泊松比抗拉强摩擦角黏聚力( kgm3 )a / mb / m量 / GPa度 / MPa/ ( )/ MPa1松散层400180010100. 010. 300. 002150. 302黏土岩、中砂岩64253820167. 900. 131. 70428. 433粉砂岩、细砂岩1202612151015. 90. 204. 534311. 23

13、4中细砂岩1162570107. 259. 820. 142. 35436. 825煤3. 31400103. 31. 900. 111. 20442. 606粉砂岩、泥岩46. 725301011. 67515. 60. 191. 90397. 592. 2. 2 模拟结果与实测对比分析开采单元距离为 20m，工作面推进过程中，每个开采单元的迭代次数采用公式 ( 1) 计算，工作面停采后，仅时间因素对迭代次数有影响，因此利用公式 ( 2) 计算。图 1 为数值模拟预计地表动态下沉曲线，图 2 为地表动态预计水平移动曲线。通过对比分析，数值模拟结果与实测数据比较吻合。3 推进速度及岩性对模拟结

14、果的影响3. 1 推进速度影响分析3. 1. 1 地表动态移动变形规律工作面推进过程中，引起上覆岩层的移动变形随时间变化，推进速度的大小，直接影响岩层移动传递的时间。按照开采单元尺寸 100m 计算，通过分析匀速推进 2m / d，10m / d 地表动态下沉曲线 ( 见图 3 ( a) 和图 3 ( b) ) 可知，随着推进速度推进距离: 1500m; 21000m; 31500m; 42000m; 52600m图 1数值模拟预计下沉与实测下沉对比的增大，相同的开采空间，地表移动变形呈变小趋势，岩层的移动时间较少，地表下沉不充分。开采速度较小时，岩层移动时间较长，其地表下沉较为充分。工作面的

15、推进速度对地表动态水平移动 ( 见105总第 118 期煤矿开采2014 年第 3 期推进距离: 1500m; 21000m; 31500m; 42000m; 52600m图 2数值模拟预计水平移动与实测水平移动对比采空区的移动变形通过岩层传递至地表，而快速推进时，采区空间向上覆岩层传递的时间较短，采空区形成较大空间的趋势加强，容易形成上覆岩层的整体下沉。低速推进覆岩移动变形更为充分。快速推进形成的开采空间上覆岩层的整体性较强，应力分布均匀，覆岩块体跨度较大。低速推进形成的开采空间上覆岩层整体性较弱，应力分布较为集中。3. 2 岩性影响分析3. 2. 1 地表动态移动变形规律通过改变岩石的物理

16、力学参数研究不同软硬程度岩性条件对移动变形的影响规律，如图 4 所示。推进距离: 1300m; 2600m; 31100m; 41700m; 52600m图 3 不同匀速开采速度地表动态移动变形曲线图 3 ( c) 和图 3 ( d) ) 具有与下沉相似的规律，即水平移动量与推进速度成反比。3. 1. 2 岩层动态移动变形规律低速推进时，整个开采过程具有充足的时间将推进距离: 1300m; 2600m; 31100m; 41700m; 52600m图 4 不同岩性地表动态移动变形曲线图 4 ( a) 和图 4 ( b) 分别为坚硬岩性和软弱岩性地表动态下沉曲线。硬岩条件下，非充分阶段下沉较小，

17、充分及超充分阶段形成整体幅度下沉，当岩性偏软时，地表下沉盆地整体较为充分，下沉幅度明显大于岩性偏硬模型。图 4 ( c) 和图 4 ( d) 分别为坚硬岩性和软弱106刘吉波等: 基于 UDEC 的岩层与地表移动动态模拟研究2014 年第 3 期岩性地表动态水平移动曲线，岩性偏硬时，地表水平移动减小，岩性较软时，地表水平移动增加。3. 2. 2 岩层动态移动变形规律软岩条件下，覆岩竖向移动具有与低速推进相似的现象，竖向位移传递速度快; 硬岩条件下，累计开采空间区域较大，造成上覆岩层整体移动。软岩条件下，采空区上方的覆岩垮落较快，应力分布集中，而硬岩条件，采空区易累计形成较大的区域，上覆岩层整体

18、性较强，应力分布均匀。4 结论( 1) 应用数值模拟软件 UDEC，将开采充分性指标和推进速度指标引入单元开采所需计算迭代次数，提出了基于开采充分程度及推进速度的动态模拟算法。建立了迭代次数的计算公式，工作面推进过程中 ni=A( 1 +v0 vi) ，工作面停采后 n =L / H02v0At 。L / H0( 2) 根据数值模拟结果可知，推进速度和岩性对地表移动变形大小及分布形态均具有显著影响。推进速度较小或覆岩岩性较软时，地表容易达到静态移动变形最大值，反之则不易达到。( 3) 数值计算分析表明，推进速度较快或岩性较硬时，采空区上覆岩层块体跨度较大，应力分布均匀，覆岩整体移动性较强; 而

19、当推进速度较慢或岩性较软时，地表下沉整体较为充分，下沉幅度大于岩性偏硬模型。参考文献1 杨天鸿，赵兴东，冷雪峰，等 . 地下开挖引起围岩破坏及其渗透性演化过程仿真 J . 岩石力学与工程学报，2003，22( S1) : 23862389.2 唐春安，徐曾和，徐小荷 . 岩石破裂过程分析 FPA2D 系统在采场上覆岩层移动规律研究中的应用 J . 辽宁工程技术大学学报 ( 自然科学版) ，1999，21 ( 5) : 456458.3 李新强，高延法，张庆松 . 开采沉陷动态数值仿真研究 J . 岩石力学与工程学报，2004，23 ( 1) : 8690.4 谢和平，周宏伟，王金安，等 . 煤

20、矿开采沉陷的岩石力学数值模拟预计方法 A . 中国岩石力学与工程学会第五届学术大会 C . 上海，1998.5 谢和平，周宏伟，王金安，等 . FLAC 在煤矿开采沉陷预测中的应用及对比分析 J . 岩石力学与工程学报，1999，18( 4) : 2933.责任编辑: 王兴库檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸檸( 上接 114 页)4 结论( 1) 矸石充填开采时，按照围岩赋存特征和充填密实度将顶板的运动破断形式主要分为直接顶垮落、直接顶弯曲下沉和坚硬顶板下层位破断。( 2) 坚硬厚顶板下矸石充填开采时，顶板的挠度远小于充填体的欠接顶量，

21、顶板与充填体间不产生力的作用，实测支架工作阻力普遍小于额定工作阻力，此时支架不再需要提供高工作阻力限制顶板下沉。( 3) 坚硬顶板破断前离层主要出现在下层位，离层量发展呈现反复离层压实的特征，且每次离层都较前一次显著增大。参考文献1 李杨，杨宝贵 .我国现代煤矿充填技术发展及其分类 J . 煤矿开采，2011，16 ( 5) : 14.2 张元功，董凤宝 . 城镇建筑群下矸石充填开采新技术的研究与实践 J . 煤矿开采，2008，13 ( 1) : 3133.3 闫少宏，张华兴 . 我国目前煤矿充填开采技术现状 J . 煤矿开采，2008，13 ( 3) : 13.4 张吉雄，缪协兴 . 煤矿矸石井下处理的研究 J . 中国矿业大学学报，2006，35 ( 2) : 197198.5 钱鸣高，许家林 . 煤炭工业发展面临几个问题的讨论 J .采矿与安全工程学报，2006，23 ( 2) : 127132.6 王家臣，杨胜利 . 固体充填开采支架与围岩关系研究 J .煤炭学报，2010，35 ( 11) : 18211826.7 缪协兴 . 综合机械化固体充填采煤矿压控制原理与支架受力分析 J . 中国矿业大学学报，201