【技术实现步骤摘要】
一种矩心矩张量反演方法和系统
[0001]本专利技术涉及矩张量反演
，更具体的说是涉及一种矩心矩张量反演方法和系统
。
技术介绍
[0002]在震源的矩张量解反演中，过于简单的源函数与存在偏差的震源位置会对矩张量反演产生难以忽略的影响，因此相对简单的震源机制解与点矩张量模型由于包含信息较少，通常不推荐作为首选；与此同时，矩心矩张量解因其包含更多的时空参数和丰富的信息，以及可以更好的近似震源，而成为主要反演方法之一
。
[0003]目前，传统的矩心矩张量反演中，在对矩心时空参数进行搜索时，往往包含一个共同假设，即震源时间函数为一个等腰三角形的时间序列，三角形的中心对应矩心时刻，底边长对应震源时间函数的总持续时间
。
[0004]该种假设对大多数时空过程简单的地震而言是可行的，但是对包含多次子事件和复杂破裂过程的地震而言是不完备的，可能会导致反演得到的震级偏小
、
也可能会导致矩张量解中出现过多的非双力偶成分
。
[0005]因此，如何克服上述问题是本...
【技术保护点】
【技术特征摘要】
1.
一种矩心矩张量反演方法，其特征在于，包括如下步骤：
S1、
根据地震波形数据确定震源区；
S2、
基于所述震源区构建自适应网格；
S3、
对所述自适应网格进行寻优，得到向量的集合，其中为矩心时间，为震源时间函数的半持续时间，为矩心空间坐标；根据向量，确定矩张量以及对应的矩心空间坐标；
S4、
对所述矩张量进行反演，得到矩张量解和震源时间函数
。2.
根据权利要求1所述的一种矩心矩张量反演方法，其特征在于，
S3
中，根据向量，按如下公式确定矩张量；式中，为自适应网格节点处观测位移第
n
个分量的位移记录，表示处与接受点之间的格林函数，表示矩心空间坐标，表示矩心时间，为震源时间函数的半持续时间，表示矩张量，
S
表示震源时间函数
。3.
根据权利要求1所述的一种矩心矩张量反演方法，其特征在于，
S4
中，根据如下公式对所述矩张量进行反演，直至到达最大迭代次数
N
，；式中，为观测数据矢量，
k
为矩张量解的个数，
n
为迭代次数，为第
k
个矩张量解，
S
为震源时间函数，表示卷积，为第一系数矩阵，
B
为第二系数矩阵，其中，；式中，为第道分量生成的格林函数矢量
。4...
【专利技术属性】
技术研发人员：张喆，
申请(专利权)人：中国地震局地球物理研究所，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。