大同煤田双系开采强矿压与地质动力环境关系研究.pdf

第 卷第期 年月山西煤炭 文章编号：（）犇 犗 犐：大同煤田双系开采强矿压与地质动力环境关系研究苏经纬，杨东辉，樊聪，李立，吕兆恒，程虹铭（山西大同大学 煤炭工程学院，山西 大同 ）摘要：矿井动力灾害是在区域地质构造运动和开采扰动的共同作用下，局部煤岩体力学系统发生的剧烈动力破坏现象。本文从地质动力环境的角度对大同煤田双系开采条件下区域地质构造活动特征、构造应力场分布、天然地震与矿井强矿压之间的关系进行了分析。研究表明：活动地块边界及其活动断裂是强矿压孕育的动力学条件，区域地质动力环境为强矿压发生提供能量基础；大同煤田的构造作用以向的近水平压应力为主，侏罗系煤层水平地应力介于 ，石炭系煤层地应力介于 ，口泉断裂伸展的正断型应力场控制着矿区“双系两硬”煤层开采的动力灾害；大同煤田东北部地震活动频繁，与复杂的地质构造正相关，且最近 年呈上升趋势；大同煤田区域 年间应变释放率 ，相当于每年释放 次级地震；区域应变能现今仍处于积累阶段，其应变值波动式增加。预测数年内地震活动性会明显增强，地质动力环境趋于活跃，矿井动力灾害将逐渐增加。关键词：大同煤田；地质动力环境；强矿压；天然地震；口泉断裂；应力场中图分类号：文献标识码：开放科学（资源服务）标识码（犗 犛 犐 犇）：犚 犲 犾 犪 狋 犻 狅 狀 狊 犺 犻 狆犫 犲 狋 狑 犲 犲 狀犛 狋 狉 狅 狀 犵犚 狅 犮 犽犘 狉 犲 狊 狊 狌 狉 犲犪 狀 犱犌 犲 狅 犾 狅 犵 犻 犮 犪 犾犇 狔 狀 犪 犿 犻 犮犈 狀 狏 犻 狉 狅 狀 犿 犲 狀 狋 狅 犳犇 狌 犪 犾 狊 狔 狊 狋 犲 犿 犕 犻 狀 犻 狀 犵 犻 狀犇 犪 狋 狅 狀 犵犆 狅 犪 犾 犳 犻 犲 犾 犱犛 犝犑 犻 狀 犵 狑 犲 犻，犢 犃 犖 犌犇 狅 狀 犵 犺 狌 犻，犉 犃 犖犆 狅 狀 犵，犔 犐犔 犻，犔 犢 犝犣 犺 犪 狅 犺 犲 狀 犵，犆 犎 犈 犖 犌犎 狅 狀 犵 犿 犻 狀 犵（犛 犮 犺 狅 狅 犾 狅 犳犆 狅 犪 犾犈 狀 犵 犻 狀 犲 犲 狉 犻 狀 犵，犛 犺 犪 狀 狓 犻犇 犪 狋 狅 狀 犵犝 狀 犻 狏 犲 狉 狊 犻 狋 狔，犇 犪 狋 狅 狀 犵 ，犆 犺 犻 狀 犪）犃 犫 狊 狋 狉 犪 犮 狋：，：收稿日期：基金项目：山西省基础研究计划（自由探索类）项目（）；山西省研究生教育创新计划项目（）；山西大同大学研究生教育创新项目（）；山西大同大学校级科研项目（）作者简介：苏经纬（），男，山西朔州人，硕士在读，主要从事矿山岩石力学领域的研究工作，：引文格式：苏经纬，杨东辉，樊聪，等大同煤田双系开采强矿压与地质动力环境关系研究山西煤炭，（）：，（）：，“”（），犓 犲 狔狑 狅 狉 犱 狊：；矿井动力灾害是一种严重的矿井煤岩系统动力失稳灾害，具有突发性强、破坏性大、机理复杂等特点。实践表明，煤层开采越深，地质构造越复杂，发生矿井动力灾害的频率越大，破坏性也越强。目前普遍认为矿井动力灾害是在区域地质构造运动和开采扰动的共同作用下，局部煤岩体力学系统发生的剧烈动力失稳破坏现象，其动力源为地层动力运动过程中所聚积的弹性能。矿井动力灾害在时间、空间和强度上的差异主要取决于矿井的地质动力环境 。最早从地球动力学的角度对冲击地压等煤岩动力灾害开展研究，并初步建立了地质动力区划分理论。张宏伟等通过分析我国部分冲击地压矿区的地质构造条件、现代地壳运动以及浅源地震特征，提出了冲击地压地质动力条件的概念。韩军等通过分析京西煤田新构造运动、构造应力场、地壳应变能等地质动力环境，利用构造反差强度对其进行定量评估，并探讨了煤岩体的破坏形式和应变能的释放条件。断层、褶皱等地质构造发育的区域往往更容易孕育发生强矿压显现，甚至矿井动力灾害。宋卫华等通过典型岩石的双剪摩擦试验，分析了北票矿区主要断层的滑动摩擦性状，采用多区域耦合有限元法对井田区域构造应力场进行了反演。李志华等采用相似模拟和数值模拟手段系统地研究了采矿活动诱发断层滑移失稳规律，为断层强矿压的现场预测和防治提供一定的指导。许多学者在构造动力灾害发生机理和防治方法等方面取得了一定的研究成果，但由于矿井地质条件的复杂性，以及采掘扰动的不确定性，强矿压发生机理研究多集中在局部简化条件下的相似模拟试验和数值计算方面，较少涉及区域地质条件下的强矿压动力环境分析或现场实测研究。而强矿压动力源范围不仅局限于巷道或回采空间，而是远大于采掘空间尺度，因此需要结合宏观尺度对强矿压的地质动力环境进行研究。基于此，本文以大同煤田双系开采为工程背景，在分析矿井动力灾害特征的基础上，结合活动地块学说，从地块活动断裂、地层动力运动、构造应力场、地震活动规律等角度阐述大同煤田的地质动力环境及其对强矿压的影响。大同煤田区域地质构造大同煤田为一轴向 不对称的开阔向斜构造盆地，长 ，宽 ，面积约 ，其区域构造位于华北断块内二级构造单元吕梁太行断块中的云冈断坳（图）。云冈断坳北部与内蒙古断块相邻，东部及南部与桑干河新裂陷为界，西北部与偏关神池断坪接壤。同忻井田位于大同煤田的东北部，毗邻桑干河新裂陷，口泉断裂是二者间的分界 线。同 忻 井 田 范 围 内 双 系 煤 层 重 叠 率 近 ，同忻矿开采 石炭二叠系煤层，开采上覆侏罗系煤层的矿井自北向南依此为晋华宫矿、忻州窑矿、煤峪口矿、永定庄矿、同家梁矿（图）。口泉断裂是大同煤田的主要边界断裂之一，全长 ，除南端走向近外，总体走向 ，倾向，倾角 ，呈现南北不对称的“”形空间展布特征。大同煤田位于云冈断坳与桑干河新裂陷的交界，是地块构造变形的产物。自晚古生代形成大同煤田以来，在侏罗纪后期，经历了印支和燕山强烈的构造运动，煤田东部地质构造复杂，断裂及岩浆活动较多，地层陡峻；北部、西北部构造则相对简单，以单一向斜为主，断层、褶皱较少，整体属简单型构造煤田。张跃刚等 研究表明，鄂尔多斯断块和山西断陷带在 年间表现为张性，第期苏经纬，等：大同煤田双系开采强矿压与地质动力环境关系研究年间表现为压性，相对运动速率分别为 、，且方向一致。由于运动存在速率差，导致两者之间的口泉断裂边界处于压缩状态，易于聚积能量，为强矿压发生提供了能量基础。大同地震台观测结果表明，口泉断裂在 年间上盘下降的平均速度为 ，至今南东盘仍持续下降，北西盘仍持续上升。应用震源位置和速度结构联合反演方法研究表明，口泉断裂中北段活动强而南北两端弱，其中北起圣水沟南至鹅毛口中段，活动性最强，可聚积较强的孕震或诱冲能量，忻州窑矿、同忻矿即位于口泉断裂的最强活动地段。图大同煤田构造纲要图 图同忻井田双系煤层布置图 忻州窑矿位于大同煤田的东北部，井田中部构造较为复杂，东西两翼相对简单。其中，以上落差断层条，的条，小于的近百条，向斜个，不对称的忻州窑大向斜贯穿整个井田，背斜个。据不完全统计，年间矿井发生强矿压 次，年间 次，年间 次，主要发生在 合并煤层，该层煤分叉合并较强烈，厚度 ，平均 ，普氏系数犳 ；顶板普遍存在有 的厚层状整体砂岩，普氏系数犳；底板存在厚砂岩，属典型“三硬”开采条件。忻州窑矿西二 层 综放工作面发育条断层，其走向与工作面倾向垂直或斜交。年工作面回采过程中发 生 次 强 矿 压，且 集 中 在 距 断 层 范围内，冲击能量 ，大部分为 数量级。可知断层附近易导致局部应力叠加，在采动扰动下诱发坚硬顶板破断，从而产生强矿压。同忻矿也位于大同煤田东北部，属大同向斜的东翼，井田为一走向 、倾向、西低东高的单斜构造形态。井田东南部毗邻口泉断裂，处于压缩状态的口泉断裂两侧聚积大量能量，形成了井田典型的地质动力环境。井田南部边界处发育一逆断层，落差较大。井田开采石炭二叠系的层煤，主采煤层厚 ，平均 ，埋深 ，采用综放开采工艺。复杂的构造背景和坚硬顶板的弹性能量聚积，使同忻矿具备发生强矿压的地质动力条件。工作面回采过程中，回风巷道顶板下沉量达 ，底鼓量达 ，帮鼓量 达 ，导致 巷 道断 面高 ，宽 ，严重影响了巷道的正常使用。强矿压与活动断裂关系前人研究 认为，中国大陆现代构造变形以地块运动为主要特征，主要强震大都发生在活动地块边界带上。有史料以来，中国大陆共记录到级以上（含级）地震 次，其中级以上特大地震 次，级大地震 次，级地震 次，不同级别地震中分布在边界带上的次数各占其总数的 、。且各活动地块边界带运动变形的应变速率与地震活动应变能释放呈线性相关，即地块变动过程中应变和应变能积累的速率越大，地震活动水平就越高，可见，活动地块的运动变形对强震总体活动水平具有控制作用。区域构造运动引起活动地块（煤岩体系统）中应力和能量的重新分布，当能量超过岩体材料破坏极限时，产生新的平衡状态。能量聚积往往容易发生在地壳的非连续变形地块，并在切割地壳浅层的活动断裂区域产生最强烈的非连续构造变形，从而诱发矿震、地震等动力灾害。据不完全统计，我国 个强矿压矿区的活动断裂及其特性见表，活动地块边界及其活动断裂容易造成应力应变集中，导致构造应力方向、应力类型发生变化，成为强震、强矿山西煤炭第 卷压孕育发生的重要地带，区域地质动力环境为强矿压发生提供了所需能量。其中，口泉断裂对大同煤田强矿压的控制作用较为明显。表部分强矿压矿区活动断裂及其特征 矿区采深活动断裂长度 断裂性质走向水平滑动速率（）代表性强震（年代地区震级）参考文献大同 口泉断裂 逆冲断层 山西大同 罗晓华，鹤岗 依兰伊通断裂 正走滑 黑龙江穆棱 武晓军，鸡西 密山敦化断裂 左旋走滑 黑龙江穆棱 王金山，营城 依兰伊通断裂 正走滑 吉林珲春南 武晓军，抚顺 密山敦化断裂 逆冲断层 辽宁海城 王金山，沈阳 营口佟二堡断裂 右旋走滑 辽宁海城 单家增，阜新 闾山松岭断裂 逆冲断层 辽宁北票 孙洪德，北票 北票朝阳断裂 逆冲断层 辽宁北票 王喜龙，开滦 唐山断裂 右旋走滑 河北唐山 周海涛，京西 黄庄高丽营断裂 正断层 北京平谷 周永恒，靖远 海原断裂 左旋走滑 宁夏海原 吴东霖，新汶 聊考断裂 正断层 山东菏泽 郭秀岩，徐州 沂沭断裂 左旋走滑 山东郯城 徐建国，淮南 郯庐断裂 右旋走滑 安徽长丰 陆元超，义马 义马硖石断裂 逆冲断层 陕西华县 曾宪涛，平顶山 秦岭北界断裂 逆冲断层 河南许昌 张建国，大同煤田构造应力场分析 地质模型本文在分析大同区域地质构造的基础上，结合研究区的构造格架和地质体主要物性，采用有限元软件 对大同煤田双系煤层的构造应力场分布进行数值模拟，旨在为强矿压的区域预测提供指导。建立地质模型时，采用二维平面模拟，不考虑垂向载荷，选用 结构单元，假定正常岩体单元和断层带为均质弹性体，遵循摩尔库伦破坏准则。正常岩石力学参数由实验室测定，断层带内岩体力学性质与周围岩体相差较大，其弹性模量一般取正常地区的 ，泊松比一般高出正常地区 。根据大同煤田主要地质构造特征，将侏罗系（埋深 ）和石炭系煤层（埋深 ）分别划分多个构造单元，各构造单元岩体力学参数如表所示。根据李长洪等 对大同煤田地应力的实测结果（表），确定模型的边界条件。由表可知：侏罗系煤层、近水平，方位角为 ，大小为 ，大小为 ；石炭二叠系煤层方位角为 ，大小为 ，大小为 。通过三角变换，侏罗系煤层东西边界施加 的载荷，南北边界施加 的载荷；石炭系煤层东西边 界 施 加 的 载 荷，南 北 边 界 施 加 的载荷。表各构造单元岩体力学参数 构造单元侏罗系煤层石炭系煤层弹性模量 泊松比弹性模量 泊松比煤体区块 岩浆侵入带 断裂带 扩展区块 表大同煤田实测主应力结果 测点埋深最大主应力中间主应力最小主应力数值 方位（）倾角（）数值 方位（）倾角（）数值 方位（）倾角（）忻州窑矿 煤峪口矿 同家梁矿 同忻矿 构造应力场分布规律图显示了大同煤田等效应力（地应力）分布特征，双系煤层主要控煤构造不同，地应力分布差别明显。侏罗系煤层东北部主要受口泉断裂带和青磁窑断裂带的控制，西部主要受西南断层带的影响。石炭系煤层整体上受口泉断裂带影响显著。模拟结果第期苏经纬，等：大同煤田双系开采强矿压与地质动力环境关系研究显示：侏罗系煤层水平地应力主要在 之间，高地应力主要分布在断层边界和顶角（可以在图中标出断层或者在上部加一个模型图），数值在 之间，不作讨论。低地应力集中在西南部和东部断层发育的地区。石炭系煤层水平地应力主要在 之间，高地应力主要分布在口泉断裂中部边界附近，数值在 之间。图大同煤田地应力分布云图 双系煤层地应力分布均呈现四周向内部逐渐增大的趋势，总体上区块边界和断层拐点处应力较大，断裂带应力较小，且构造曲率越大，水平应力越小。说明地应力分布主要受地质构造分布的控制。石炭系煤层地应力整体上高于侏罗系煤层，主要是由于石炭系煤层埋深大于侏罗系煤层，煤层水平应力随埋深增大而有所增加。构造主应力迹线方向图显示了大同煤田双系煤层局部区域最大主应力迹线分布特征，最大主应力以压应力为主，整体向展布，与研究区边界张拉方向存在一定夹角，局部地区差异性明显。侏罗系煤层总体上最大主应力方向以为主，忻州窑矿和口泉断裂中上段最大主应力为和向，与地应力实测结果一致。东部和南部局部区域由于受到口泉断裂带及煤层相变的影响，最大主应力转为 和 向，且分布范围有限，对煤层总体应力分布影响不大。（）侏罗系煤层（）石炭系煤层图主应力迹线模拟结果 石炭系煤层总体上最大主应力方向以 为主，同忻矿和煤层中部最大主应力均为 向，南部区域由于受到口泉断裂带及其它断层带的影响，最大主应力转为向，对煤层东部及南部应力分布影响较大。结合地质动力环境分析，大同煤田水平挤压应力方向基本为，后期向转动，是对燕山运动时期构造应力水平挤压应力的继承和发展。大量震源机制解结果显示 ，大同地区存在走向为 和的两条地震节面，其中 向节面与大同盆地走向一致，另一条与其斜交。大同煤田现代构造应力场主压应力为 向，方位角 ，仰角变化较大，在 之间，主张应力为向，仰角较小。而采用有限元法模拟得出的大同煤田构造主应力迹线方向，不同于华北区域性山西煤炭第 卷 向挤压的走滑型应力场，体现了构造应力场的分区性，但与实测地应力方向具有明显的一致性，且受口泉断裂带的控制和影响较为显著。从煤层赋存上来讲，侏罗系煤层标高在 ，石炭系在 或更深，层间距 ，局部约 ，而口泉断裂东侧地表标高在 。侏罗系煤层由于未受到口泉断裂东侧地层的约束，最大主应力方向发生了较大偏转，聚积的弹性能相对较低；石炭系煤层由于受到口泉断裂东侧地层的约束，聚积大量弹性能，导致水平主应力急剧增大，构造活动性明显增强，发生矿井动力失稳破坏的频率和强度显著增加。由此可见，口泉断裂伸展的正断型应力场控制着矿区“双系两硬”煤层开采的动力灾害。大同煤田区域地震特征 区域地震分布研究认为，地震多发生在构造应力作用复杂、强烈地区和应力区边界，其次在局部应力变化区和活动断裂带上。大同煤田属于山西地震带的一部分，山西地震带是华北地震活动较强的地区之一。历史上大同煤田东侧的口泉断裂曾发生次古地震事件，其中次距今 、，另一次距今 ，平均间隔约 ，反映了口泉断裂可能具备准周期的强震活动。整理国家地震科学数据共享中心（：）大同煤田及其邻域（、）在 年间 级的地震记录（年来源于中国地震台网地震目录，年来源于中国地震台网统一地震目录），总计 次，其 中 级 地 震 次，占 总 数 的 ，级地震次，仅 年便发生次，最大为 月 日大同 级。通过 平台对大同煤田区域地震目录和地貌图进行叠加，展现出地震活动和活动断裂的空间分布关系（图）。六棱山北麓和大同煤田东北部地震分布较为密集，与复杂的地质构造呈现显著的正相关性。大同矿区典型动力灾害矿井忻州窑矿、煤峪口矿、同家梁矿、同忻矿等，均位于地质构造条件复杂的大同煤田东北部。可见地震活动性越强，构造运动就越剧烈。根据中国地震台网统一地震目录，大同煤田及其邻域在 年月日 年 月 日共发生 级地震 次，其中天然地震占 ，塌陷型地震占 ，爆破型地震占 ，可知煤炭开采活动诱发地震占总数的将近一半，其中采空区顶板垮落极易诱发地震活动。图显示了大同煤田区域地震时间分布具有一定的规律性，不同年度地震 频 发 交 替 出 现，主 要 集 中 在 、和 三个时间段，最近 地震活动性呈现上升趋势，地质动力环境活跃程度也明显增强，为矿井动力灾害的发生提供能量基础。历史上山西地震频发，且大部分强震发生在断陷盆地。大同煤田虽未发生过级以上地震，但中小地震不断，且空间分布具有一定的规律性，主要发生在边界断裂和隐伏断裂带上。究其原因，主要是均匀区域构造应力场发生了局部应力变化。我国华北地区构造应力场总体上以 向挤压走滑型为主，而山西断陷带以向拉张为主，大同煤田则以向挤压为主，与华北区域构造应力场相比，应力方向和结构类型发生了明显的变化。图大同煤田及其邻域地震分布 图大同煤田及其邻域年度地震次数统计 第期苏经纬，等：大同煤田双系开采强矿压与地质动力环境关系研究 地震应变积累释放特征大同煤田及其邻域地震资料比较完整，以 年为起点，绘制该区域 年地震应变积累释放曲线及 图（图）。在 总计 大 同 煤 田 区 域 平 均 应 变 释 放 率 ，相当于每年释放 次级地震，与华北地区平均应变释放率 相比，应变积累略低，地震活动性稍弱。图（）显示的应变积累释放曲线大致反应了一个平静活跃周期，从 年 开 始 经过次 级 以 上 地 震 和 年次 级地震，应变能释放基本完成，其应变值达到最低，而后应变能一直处于积累阶段。该区域至今应变值波动式增加，很快达到积累阶段峰值并释放出来，预测数年内地震活动性会明显增强，地 质动 力环 境 趋 于活 跃，矿 井 动 力灾 害逐渐增加。（）年应变释放曲线（）年 图图大同煤田及邻域应变释放及 图 结论）矿井动力灾害的发生是区域地质构造运动和开采扰动共同作用的结果。活动地块边界及其活动断裂是强矿压孕育的动力学条件，区域地质动力环境为强矿压显现提供了所需能量。）大同煤田的构造作用以向的近水平压应力为主，侏罗系煤层水平地应力介于 ，石炭系煤层地应力介于 ，口泉断裂伸展的正断型应力场控制着矿区“双系两硬”煤层开采的动力灾害。）大同煤田地震活动具有明显的时空分布特征。空间上煤田东北部地震活动频繁，与复杂的地质构造呈现显著的正相关性；时间上不同年度地震频发 交 替 出 现，主 要 集 中 在 、和 三个时间段，且最近 地震活动性呈现上升趋势。）大同煤田区域 年间应变释放率 ，相当于每年释放 次级地震；区域应变能现今仍处于积累阶段，其应变值波动式增加，很快达到峰值并释放出来。预测数十年内地震活动性会明显增强，地质动力环境趋于活跃，矿井动力灾害逐渐增加。参考文献：袁亮深部采动响应与灾害防控研究进展煤炭学报，（）：，（）：张宏伟，朱峰，韩军，等冲击地压的地质动力条件与监测预测方法煤炭学报，（）：，（）：姜耀东，潘一山，姜福兴，等我国煤炭开采中的冲击地压机理和防治煤炭学报，（）：，（）：韩军煤矿冲击地压地质动力环境研究煤炭科学技术，（）：，（）：，山西煤炭第 卷 ：，：张宏伟，韩军，宋卫华，等地质动力区划北京：煤炭工业出版社，：韩军，张宏伟，兰天伟，等京西煤田冲击地压的地质动力环境煤炭学报，（）：，（）：宋卫华，张宏伟断层滑动准则和稳定性评价及构造应力场的反演岩石力学与工程学报，（）：，（）：李志华，窦林名，陆振裕，等采动诱发断层滑移失稳的研究采矿与安全工程学报，（）：，（）：岑敏，董树文，施炜，等大同盆地形成机制的构造研究地质论评，（）：，（）：张跃刚，胡新康华北地区块体及其边界的相对运动大地测量与地球动力学，（）：，（）：王霞，宋美琴，王亮，等口泉断裂及其邻近地区的地壳速度结构地震地质，（）：，（）：陈蓥同忻矿地质动力环境分析与矿压显现规律研究阜新：辽宁工程技术大学，张培震中国大陆岩石圈最新构造变动与地震灾害第四纪研究，（）：，（）：张培震，邓起东，张国民，等中国大陆的强震活动与活动地块中国科学：地球科学，（）：，（），（）：张培震，邓起东，张竹琪，等中国大陆的活动断裂、地震灾害及其动力过程中国科学：地球科学，（）：，（），（）：张国民，马宏生，王辉，等中国大陆活动地块边界带与强震活动地球物理学报，（）：，（）：陈蓥同忻矿地质动力环境分析与矿压显现规律研究阜新：辽宁工程技术大学，李长洪，张吉良，蔡美峰，等大同矿区地应力测量及其与地质构造的关系北京科技大学学报，（）：，（）：谢富仁，陈群策，崔效锋，等中国大陆地壳应力环境基础数据库地球物理学进展，（）：，（）：赵德安，陈志敏，蔡小林，等中国地应力场分布规律统计分析岩石力学与工程学报，（）：，（）：谢富仁，张红艳，崔效锋，等中国大陆现代构造应力场与强震活动国际地震动态，（）：，（）：马宏生，张国民，刘杰，等中国大陆及其邻区强震活动与活动地块关系研究地学前缘，（）：，（）：马宏生，刘杰，张国民，等用应变积累释放模型研究中国大陆地块分区地震活动地震学报，（）：，（）：（编辑：单婕）第期苏经纬，等：大同煤田双系开采强矿压与地质动力环境关系研究