[发明专利]岩石声发射的混合矩张量反演计算方法、存储介质和系统

权利要求书

1.一种岩石声发射的混合矩张量反演计算方法，其特征在于，包括：

利用最小二乘法进行绝对矩张量反演计算，得到修正前的震源机制解；

引入迭代参数，根据迭代参数计算权重修正值；

通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解。

2.根据权利要求1所述的岩石声发射的混合矩张量反演计算方法，其特征在于，

利用最小二乘法进行绝对矩张量反演计算，得到修正前的震源机制解的公式为：

M_i＝G_i^-1u_i^obs，其中G_i为格林函数，u_i为接收到的弹性波的波形位移振幅，obs代表观测值，M_i为震源机制解。

3.根据权利要求1所述的岩石声发射的混合矩张量反演计算方法，其特征在于，

引入迭代参数，根据迭代参数计算权重修正值至少包括：

引入第一迭代参数和第二迭代参数，第一迭代参数的计算公式为：

其中为第一迭代系数，median为取中位值，u_ijkl^th为弹性波的波形位移振幅理论值，u_ijkl^obs为弹性波的波形位移振幅实际值，N_eq是指所需计算的方程总数。

4.根据权利要求3所述的岩石声发射的混合矩张量反演计算方法，其特征在于，计算第二迭代参数的公式为：

其中，ω_IterNo为第二迭代参数，IterNo为自然数。

5.根据权利要求4所述的岩石声发射的混合矩张量反演计算方法，其特征在于，根据迭代参数计算权重修正值的公式为：

其中，a_jkl为权重修正值。

6.根据权利要求5所述的岩石声发射的混合矩张量反演计算方法，其特征在于，

通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解包括：先计算修正后的弹性波的波形位移振幅实际值，再计算修正后的震源机制解；

修正后的弹性波的波形位移振幅实际值的计算公式为：

u_ijkl^new＝u_ijkl^old+a_jkl·u_ijkl^old

其中，u_ijkl^old为弹性波的波形位移振幅实际值，u_ijkl^new为修正后的弹性波的波形位移振幅实际值。

7.根据权利要求6所述的岩石声发射的混合矩张量反演计算方法，其特征在于，

计算修正后的震源机制解的计算公式为：

M_i^new＝G_i^-1u_i^new其中，u_i^new为修正后的弹性波的波形位移振幅实际值，M_i^new为修正后的震源机制解。

8.根据权利要求7所述的岩石声发射的混合矩张量反演计算方法，其特征在于，

利用公式u_i^th＝G_iM_i重新计算弹性波的波形位移振幅理论值，重复引入迭代参数，根据迭代参数计算权重修正值，通过权重修正值对修正前的震源机制解进行修正，得到修正后的震源机制解的步骤，通过多次修正得到多次修正后的震源机制解。