煤矿微震震源参数反演及震源破裂机理研究
发布时间：2020-08-01 17:14
【摘要】：我国煤矿煤岩动力灾害非常严重。随着开采深度的不断增加,煤岩动力灾害日趋严重及复杂,对深部资源的安全高效开采造成重大影响。微震监测技术作为一种在煤矿中广泛应用的地球物理监测手段,在矿震及冲击地压的震源定位和能量分析、冲击地压演化过程以及监测预警等方面能发挥重要作用。然而,煤系地层非均质性强且多变,波场结构复杂,微震定位精度和稳定性难以满足现场需求,国内外对煤矿微震震源参数及震源破裂机理缺乏深入系统的研究。本文采用实验室实验、理论分析、数值模拟和现场试验验证等方法,以复杂煤岩介质条件下波场传播特征为切入点,深入研究基于高精度震源定位的震源参数和震源破裂机理,并进行现场验证与应用。论文得到的主要研究成果如下:针对煤系地层复杂波场结构,基于声波方程的空间域4阶、时间域2阶高阶有限差分法,对各种复杂煤岩介质条件下的波场传播特征(波场快照和单炮记录)进行了模拟分析。结果表明:在不同的介质分界层产生的反射波、透射波以及各种转换波的形状都不同;在断层处会形成一个新的震源并产生折射波和绕射波;在双相介质中会产生慢P波和慢S波;在随机离散介质下,会产生复杂的散射波。为引入震源联合定位法和求取震源参数(尤其是品质衰减因子Q)打下基础。提出了单纯形和双差联合微震定位方法,对煤矿现场微震震源位置进行了校正,提高了震源定位的精度。结果表明:单纯形和双差联合定位法给出的爆破事件的震源定位误差控制在20m以内,事件残差控制在15ms以内,验证了联合定位法的精度;通过联合定位法对煤矿微震事件重定位后的均方根残差由1.35s降到0.62s,震源定位精度提高,该方法有效地减少了地壳速度结构模型不精确导致的误差,保证了震源定位和求取震源参数的精度。基于ω~2模型,系统地分析了实验室小尺度煤岩体破裂和煤矿微震的震源参数(拐角频率f_0、品质衰减因子Q、震源能量E、地震矩M_0、震源半径R、应力降△ζ、视应力ζa)和震级,揭示了震源大小和应力状态等震源性质及震源参数随地震矩的变化特征。结果表明:小尺度破裂的震源参数f0和ζa随M_0的增加线性减小,震源参数R和E随M_0的增加线性增加,△ζ整体上围绕着某一个值上下波动;老虎台煤矿微震的震源参数f_0随M_0的增加线性减小,震源参数R、E和△ζ随M_0的增加线性增加,但是△ζ与M_0的线性关系较弱,ζa几乎围绕某一值上下波动,而且老虎台煤矿微震的f_0和△ζ值明显小于天然地震对应震源参数的值。研究了小尺度破裂和煤矿微震的震源破裂机制、破裂尺度和破裂面上的滑动分布,揭示了震源破裂机理。结果表明:通过矩张量反演法反演出的震源沙滩球,清楚地判别了小尺度破裂和煤矿微震震源破裂存在的剪切破坏、剪切张拉破坏和剪切挤压破坏特征。根据建立的震源破裂尺度与矩震级的关系式求取了震源的破裂面积、长度、宽度和滑动分布,对震源的破裂范围进行了系统地评价。基于求取的破裂尺度和k2滑动破裂模型得出了不同震级破裂事件的震源滑动分布,揭示了震源的滑动破裂过程。开发了“微震震源参数反演与震源破裂机理分析软件”,实现了对微震监测波形数据的自动处理、定位计算、震源参数反演及破裂机理分析。揭示了千秋煤矿冲击地压破坏的震源参数特征和机制,得到了4次冲击地压事件的震源破裂尺度和滑动分布数据,结果与现场实际情况基本吻合。揭示了矿震破裂面上滑动大的区域容易诱发冲击地压,因此可基于破裂面上滑动分布,对冲击地压进行有效的分区域防治。研究成果对于进一步提高煤矿微震震源定位的精度、揭示矿震及冲击地压灾害震源破裂机理、准确评价和预测冲击地压、有效防治冲击地压灾害具有重要的理论意义和应用价值。
【学位授予单位】：中国矿业大学
【学位级别】：博士
【学位授予年份】：2019
【分类号】：TD326
【图文】：
以推导处时间域 2 阶空间域 4 阶的波动方程有1, 1, . 1 , 1 2, 2, ,1, ,) (k k k k k k k m n m n m n m n m n m n m k km n m n U U U U U U U+ + + + + + + + ，其它符号意义同前。限差分法的不同煤岩介质波场特征分Analysis of Different Coal and Rockte Difference Method）介质模型波场模拟质模型波场模拟层介质模型，VP=1300m/s，VS=751m/s，ρ=1400 s，网格数为 1600×1600，模拟震源位置为（0,在 t=80ms 时的波场快照如图 2-1 所示。
质模型波场模拟层介质模型，VP=1300m/s，VS=751m/s，ρ=1400s，网格数为 1600×1600，模拟震源位置为（0在 t=80ms 时的波场快照如图 2-1 所示。图 2-1 单层煤介质（X 和 Y 分量）波场快照Wavefield snapshots of single layer coal media (X and Y看出，在单一均匀煤介质模型中，震源在中心周传播，P 波和 S 波的波前是一个同心圆，纵横质模型波场模拟
VP=3500m/s，VS=2021m/s，ρ=2600网格数为 1600×1600，模拟震源位置为（0, t=80ms 时波场快照如图 2-2 所示。，单层岩介质和单层煤介质类似，同样为单播，P 波和 S 波的波前也是一个同心圆，波层介质的速度明显大于煤层介质，所以在相层模型中的传播速度要明显的大于单一均匀波场模拟型波场模拟模型 I模型 I，上层（（-800）~100）为煤层 VP=1~800）为岩层 VP=3500m/s，VS=2021m/s，ρ=2网格数为 1600×1600，模拟震源位置为（0, 子波频率为 30Hz。此模型 I 在 t=80ms 时的照如图 2-4 和 2-5 所示，单炮记录如图 2-6
本文编号：2777758
【学位授予单位】：中国矿业大学
【学位级别】：博士
【学位授予年份】：2019
【分类号】：TD326
【图文】：
以推导处时间域 2 阶空间域 4 阶的波动方程有1, 1, . 1 , 1 2, 2, ,1, ,) (k k k k k k k m n m n m n m n m n m n m k km n m n U U U U U U U+ + + + + + + + ，其它符号意义同前。限差分法的不同煤岩介质波场特征分Analysis of Different Coal and Rockte Difference Method）介质模型波场模拟质模型波场模拟层介质模型，VP=1300m/s，VS=751m/s，ρ=1400 s，网格数为 1600×1600，模拟震源位置为（0,在 t=80ms 时的波场快照如图 2-1 所示。
质模型波场模拟层介质模型，VP=1300m/s，VS=751m/s，ρ=1400s，网格数为 1600×1600，模拟震源位置为（0在 t=80ms 时的波场快照如图 2-1 所示。图 2-1 单层煤介质（X 和 Y 分量）波场快照Wavefield snapshots of single layer coal media (X and Y看出，在单一均匀煤介质模型中，震源在中心周传播，P 波和 S 波的波前是一个同心圆，纵横质模型波场模拟
VP=3500m/s，VS=2021m/s，ρ=2600网格数为 1600×1600，模拟震源位置为（0, t=80ms 时波场快照如图 2-2 所示。，单层岩介质和单层煤介质类似，同样为单播，P 波和 S 波的波前也是一个同心圆，波层介质的速度明显大于煤层介质，所以在相层模型中的传播速度要明显的大于单一均匀波场模拟型波场模拟模型 I模型 I，上层（（-800）~100）为煤层 VP=1~800）为岩层 VP=3500m/s，VS=2021m/s，ρ=2网格数为 1600×1600，模拟震源位置为（0, 子波频率为 30Hz。此模型 I 在 t=80ms 时的照如图 2-4 和 2-5 所示，单炮记录如图 2-6
本文编号：2777758
本文链接：https://www.wllwen.com/kejilunwen/kuangye/2777758.html
