[发明专利]一种煤岩冲击失稳模式的微震识别方法

摘要

一种煤岩冲击失稳模式的微震识别方法，属于冲击矿压微震特性分析技术领域。其特征在于：包括如下步骤：步骤1001，在同一采场不同位置设置多个微震记录仪对采场进行微震监测，每个微震记录仪形成一个微震信号记录；步骤1002，在微震记录仪监测到的微震记录中，筛选多个微震波形；步骤1003，根据震源远场位移计算步骤，计算每个微震波形的震源远场位移；步骤1004，计算得到震源矩张量；步骤1005，对步骤1004中计算得到的震源矩张量进行分解；步骤1006，对步骤1005中分解后的震源矩张量进行分析得到：震源破裂模式，震源P、T轴以及震源破裂面产状。

主权项

1.一种煤岩冲击失稳模式的微震识别方法，其特征在于：包括如下步骤：步骤1001，在同一采场不同位置设置多个微震记录仪对采场进行微震监测，每个微震记录仪形成一个微震信号记录；步骤1002，在微震记录仪监测到的微震信号记录中，根据波形筛选原则筛选出多个微震波形；步骤1003，根据震源远场位移计算步骤，计算每个微震波形的震源远场位移；步骤1004，计算得到震源矩张量；步骤1005，对步骤1004中计算得到的震源矩张量进行分解；步骤1006，对步骤1005中分解后的震源矩张量进行分析得到：震源破裂模式，震源P、T轴以及震源破裂面产状；所述的震源破裂面产状的求取方法包括如下步骤：步骤b1，震源矩张量可以由破裂面的位置与运动矢量表示为：Mij＝uS[λvknkδij+μ(vinj+vjni)]式中：i，j，k＝1，2，3，为空间坐标系中的3个方向；Mij为地震矩张量形式，根据i，j取值不同，分别代表矩张量不同分量；u为破裂面运动方向的位移量；S为破裂面表面积；λ和μ为拉梅常数；v为运动方向，割据i,j,k取值不同，vi,vj,vk表示破裂面运动矢量在各坐标系方向上的分量；n为破裂面法向方向，其分量取值意义通v；步骤b2，对步骤b1的公式进行本征值化，得到：则有(假设M1>M2>M3)：M1+M3‑2M2＝2μuSnvM1‑M3＝2μuS；步骤b3，根据步骤b2得到的特征矢量与破裂面运动方向和方向的关系：式中，e1⊥e2⊥e3，绝对值符号表示矢量大小；×表示矢量乘法；e1、e2、e3分别为矩张量的最大特征值、中间特征值和最小特征值；步骤b4，设步骤b3中得到的矢量v和矢量n的夹角为β，v和e1的夹角以及n和e1的夹角均为β/2，得到：并得到破裂面运动方向和法向方向与矩张量最大特征值对应特征矢量和最小特征值对应特征矢量的关系：步骤b5，根据破裂面法向方向的空间矢量值到破裂面的几何方程表达式，并确定破裂面的方位角和倾角的计算公式：方位角：式中，n(3)为破裂面法向矢量Z轴分量，|n|表示法向向量的模若dip>90°，则破裂面方位角为取值为(π‑dip)倾角:式中，n(1)为破裂面法向向量X轴分量，n(2)为破裂面法向向量Y轴分量；在执行完所述的步骤1006后，绘制微震的震源球；步骤1003中所述的震源远场位移计算步骤，包括：步骤a1：根据所述的步骤1002中筛选出的微震波形，拾取每一个微震波形的P波初振幅值，并记作U1；步骤a2，在筛选的微震波形中剪切P波的时域波形；步骤a3，根据微震记录仪固有采样频率将步骤a2中剪切的P波时域波形进行快速傅里叶变换，转化为频域波形；步骤a4，对步骤a3中转化得到的频域波形利用如下公式进行衰减修正：其中，A(f)为时域速度谱进行快速傅里叶变换的结果；f为步骤a3中转换得到的频域波形的频率；v为P波速度；Q为衰减因子；步骤a5，在步骤a4中计算得到的衰减修正值得到地动速度功率谱并根据地动速度功率谱得出地动位移功率谱步骤a6，计算得到地动速度功率谱积分SV2和地动位移功率谱积分SD2：其中：为考虑自由面影响乘以1/4的修正速度功率谱；为对应的位移功率谱；步骤a7，计算得到各个微震波形的低频位移U2：