应用离散元法分析采矿引起厚松散层变形的特征.pdfVIP

第 27 卷第 3 期 煤 炭 学 报 Vol.27 No.3 2002 年 6 月 JOURNAL OF CHINA COAL SOCIETY June 2002 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 文章编号：0253 - 999（3 2002）03 - 0268 - 05 应用离散元法分析采矿引起厚松散层变形的特征 许延春，张玉卓 （天地科技股份有限公司 开采所事业部，北京 100013） 摘 要：介绍了利用离散元法模拟分析从采矿到基岩垮落至厚松散层变形的过程，分析了厚松散 层的变形特征与开采充分度、采厚、松散层厚度和松散层岩性等的相关关系 . 认识到对永久性非 充分采动，相关建筑物保护煤柱可适当减小 . 厚松散层内部移动边界有明显的曲线特征，应选用 分段或曲线形式设计建筑物保护煤柱 . 厚松散层厚度越大则移动角越小，一次采高大则松散层变 形量大 . 关键词：离散元；厚松散层；变形形态特征 中图分类号：TD325.2 文献标识码：A 厚松散层普遍分布于华东地区的各个矿区，例如淮南、淮北、兖州、徐州和大屯等 . 对华东地区众多 村庄保护煤柱的合理留设有很大影响 . 近年来由于放顶煤开采使一次采高显著增大，采用条带法进行村庄 下采煤情况增多以及该地区出现井壁集中破坏现象等原因，人们认识到研究厚松散层的性质、变形特征是 十分重要的 . 由于松散层内部和地表变形十分复杂，因此，希望通过现场观测的方法全面了解松散层的变 形在技术上难以达到，经济上也是难以承受的；采用有限元方法，一是难以合理模拟采场顶板岩体垮落、 离层等大变形；二是在模拟地表下沉盆地边缘的弯曲和拉伸变形时，会产生较大误差，而该部分正是研究 的重点所在. 本文利用离散元适用于模拟、分析大变形和不连续介质受静、动载荷作用的固体稳定性问题 的特点，建立了从采场到地表的全采场剖面分析模型，从而全过程模拟了采矿对厚松散层变形的影响 . 本 次模拟分析的主要目的是通过模拟采矿引起覆岩的变形，研究不同采矿条件下和不同松散层条件下的松散 层、地表的移动变形特征 . ! 离散元分析的方案和模型 采用正交分析法设计分析方案，共制定 11 个模拟方案（表 1）. 以方案 1 为例说明数学模型（图 1）. 由图 1 可见，模型为全采场剖面平面应变模型 . 即从开采煤层、基岩、松散层到地表的整个采场覆岩结 构 . 地层长度为 600 m，基岩厚为 50 m，松散地层厚度分别为 50，150 和 200 m. 在基岩段设置 45 和 135 两个方向节理面 . 节理和单元的最大边长为 10 m. 在松散层段采用垂直节理，节理单元最大边长为 20 m. 在离散单元的内部由程序自动划分有限单元，并按弹性力学进行力学分析 . 破坏准则采用 Domgrapher，循 环次数为 13 万步 . 对节理面输入法向刚度、切向刚度、黏聚力和内摩擦角 . 根据土工试验获得的性质参 数对松散层输入 . 地层的岩性结构见表 2. 计算时首先以方案 1 模拟类似于兴隆庄矿 2308 - 1 工作面地表 移动观测结果条件，然后进行计算分析 . 深厚松散层变形的形态特征 收稿日期