古构造应力场数值模拟方法探讨——以山西省晋城成庄煤矿区为例.pdfVIP

第 15卷第 3期 地 质 力 学 学 报 V01．15No．3 2009年 9月 JOURNALOFGE0MECHANICS Sep．2009 文章编号 ：1006．6616 (20O9)03—0270．11 古构造应力场数值模拟方法探讨 — — 以山西省晋城成庄煤矿 区为例 周春梅 ，章泽军 ，李先福 (1．武汉工程大学环境与城市建设学院，湖北 武汉 430073； 2．中国地质大学地球科学学院，湖北 武汉 430074) 摘 要：从构造模拟 目标、构造模型建立、构造模拟三个方面探讨 了古构造应力场 数值模拟的方法。以山西省晋城成庄煤矿 区为例，根据野外调查和微观构造形迹， 建立平面应力模型和三维地质模型，确定岩体力学参数、边界条件及岩石破裂准 则，采用大型数值模拟软件 ANSYS和FLAC计算，确定外力的大小及加载顺序。通 过分析研究，揭示了山西晋城成庄煤矿区历史上 曾发生先东西向、后南北 向的两期 大的构造运动，加载外力为先东西向加载60MPa，再南北向加载 110～180MPa，这 个结果能解释 多数地表剪节理和 3 煤中剪节理 同时发生。该分析对矿 区的危险性 分区及瓦斯突出区的预测具有重要的意义。 关键词 ：古构造应力；构造模拟；平面应变；剪节理 中图分类号：TU452 文献标识码：A 0 引言 古构造应力场模拟是根据地质历史时期构造作用所 留下的构造形迹 ，恢复或反演古构造 应力作用方式、方向、大小、边界条件及其介质所处环境 、力学性质在内的构造作用过程。 在对古构造应力场数值模拟分析中，关键问题是地质建模时如何考虑边界条件的应力量值大 小和应力加载顺序，同时必须保证数值模拟体与地质原型问满足相似性原理 。。 构造模拟大致可分为确立构造模拟 目标、抽象或建立构造模型、构造模拟三部分。 构造模拟 目标可以定义为判别构造模拟成功与否的重要标志，通常以野外地质调研为主 并综合其他地质信息确定包括构造变形序列与递次变形过程在内的最终变形结果或样式即为 构造模拟 目标。需要说明的是实现最终变形结果或样式的途经不是唯一的，只能说它是一种 可能的地质作用过程 ，是一种合理的解释。 建立构造模型以构造模拟 目标为依据，抽象力学模型，包括介质 (地质体)的形态、力 学参数的取舍、不同介质体之间的关系与边界条件以及加载外力的作用方式的确定 ，对于多 期构造变形问题需要考虑分阶段的边界条件与外力的作用方式 ，构造模型能否合理地建立是 收稿 日期 ：2009．04．02 基金项 目：国家 自然科学基金项 目 40672187)；武汉工程大学青年基金项 目 (Q2oo8o5)。 作者简介 ：周春梅 (1979．)，女，博士，讲师，主要从事矿山构造应力研究。E-mail：zhouehunmei@163．tom 第 3期 周春梅等 ：古构造应力场数值模拟方法探讨 271 构造模拟的关键。 构造模拟是 “以匀代不匀”的基本思想 ，用有限元法、有限差分法等方法，以海量、大 型计算机为工具实施构造模型中不同阶段或过程的运算 ，提供不同阶段的定量模拟数据。需 要指出的是地质历史时期的构造变形具有两大特点：其一是构造变形速率极为缓慢，是一个 漫长的时间过程 ，在这一漫长的变形过程中，变形状态由弹性转化为与时间有关 的流变状 态；另一特点是构造变形多数或严格地说无一例外的均属大变形问题 J，这种与时间有关的 流变以及大变形 问题是数值模拟中迄今为止尚待探索的难点。 1 构造模拟 目标 山西晋城成庄矿区位于华北断块 中的吕梁一太行断块沁水块坳东部次级构造单元沾尚一 武乡一 阳城 NNE向褶皱带南段。矿区地表仅 出露上古生界下石盒子组 (P，X)石英砂岩至石 千峰组 (P2sh)粉砂岩与细砂互层 ，受区域性构造的制约 ，总体为一走向北北东 (北部)逐 渐转折为北东向 (南部)、地层平缓 (倾角 3。一