申请/专利权人：北方工业大学

申请日：2022-01-20

公开（公告）日：2024-03-29

公开（公告）号：CN114417612B

主分类号：G06F30/20

分类号：G06F30/20;G06F17/16;G06F119/14

优先权：

专利状态码：有效-授权

法律状态：2024.03.29#授权;2022.05.20#实质审查的生效;2022.04.29#公开

摘要：本发明涉及一种基于矩张量反演的采场微震震源机制求解方法，利用本发明方法进行程序开发，编制震源机制求解器，通过输入P波初动振幅、震源定位点坐标及传感器坐标直接进行矩张量求解，进而计算出震源破裂类型、破裂倾角、裂纹运动方向、破裂走向及相对裂缝体积等震源参数，并绘制震源机制的矩张量沙滩球图，实现结果可视化。本发明方法，在确定震源定位坐标、传感器位置以及P波初动振幅的条件下即可进行计算，能够方便且准确的对震源破裂类型进行判别，方便快捷的进行震源参数求解，绘制震源机制沙滩球，求解结果准确且直观。

主权项：1.基于矩张量反演的采场微震震源机制求解方法，其特征在于，具体步骤如下：1利用震源定位点坐标、传感器探头坐标及P波初动振幅进行矩张量求解：根据弹性波动力学理论，P波远场位移U由格林函数矩阵和矩张量来表示：U＝GM*δt1式中：G为弹性动力学格林函数矩阵，M为震源矩张量，δt为震源时间函数；通过引入反射系数Re，并引入脉冲函数从而忽略时间函数δt，得到在半无限空间条件下，通过SIGMA法简化的P波初动振幅Zx的表达式： 其中，C为传感器校正系数；R为震源距，表示震源和传感器探头之间距离；γ1、γ2、γ3表示震源和传感器探头之间的方向余弦；Mij表示矩张量的各个元素；Re为界面的反射系数，表达式如下： 其中，p＝vpvs，表示纵波vp和横波vs速度之比；e＝r·t，表示震源到传感器探头的方向向量r与传感器探头的感度方向t的乘积；2震源参数求解a.根据矩张量特征值进行分解，得到各部分所占比例值，判断震源破裂类型： 其中，通过对矩张量Mij进行分解，得到M1、M2、M3为矩张量的三个特征值，其中M1为最大特征值，M3为最小矩张量，Q代表剪切部分比例值，T代表张拉偏量部分比例值，F代表静水拉应力部分比例值；b.计算破裂走向、裂纹运动方向及破裂倾角：矩张量的三个特征值M1、M2、M3对应的特征向量m1、m2、m3与裂纹的位错向量v、破裂面的法线向量n有如下关系式： c.计算相对裂纹体积参数Vr：Vr＝M1+M2+M3v63绘制震源机制沙滩球。

全文数据：

权利要求：

百度查询： 北方工业大学 基于矩张量反演的采场微震震源机制求解方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD