，网格大小为，并采用模型计算理论走时，通过计算每个时间窗的聚束能量得到如图所示的时间段内的相关系数图振幅比图聚束能量图以及以及它们的点积函数图的相对大小参见杜海林等张喆等，并基于函数，我们以峰值的作为阈值，确定此次地震台站图提取视震源时间函数时，截取直达波前，后，以的矩心矩张量参数计算相应的格林函数，依据波形复杂性与持续时间采用了的频率范围，并根据经验设定最大持续时间次的迭代，最后获得如图所示的视震源时间函数可以看出这些视震源时间函数彼此之间具有相似性，但也随着台站的方位具有规律性变化这不仅表明震源破裂的复杂性，也表明破裂具有明显的破裂优势方向图依赖于台站方位的视震源时间函数分析为了获得总体的破裂优势方位，我们利用各台站视震源时间函数的最大峰值时间以及台站相对于震中的离源角，基于模型，通过公式进行了最小乘拟合，得到最优破裂方位角上式中，和为拟合参数，为台站相对于震中的方位角为了揭示依赖于方位的破裂复杂性的更多细节，我们将视震源时间函数投影到以时间为径向台站方位角为切向的极坐标平面图除总体的破裂优势方向和地震矩的分段，矩心深度约图，同时可以看到，从仪器震中位臵得到的矩张量解也不同于最终反演得到矩心矩张量解图反演结果表明，整个地震释放标量地震矩，相应的矩震级为，其中双力偶成分占比图表列出了矩心位臵和矩张量解参数，并与的相应参数进行了比较比较发现，我们的位臵在矩心位臵的东北，矩心深度略浅，但矩张量解非常相似与等，以及的最佳双力偶解的比较展示于表尽管已经发布的结果之间存在定差异，但可以确定，年南斯科舍海岭近正东的断层的左旋走滑为主，倾角在至之间和我们的矩张量解中都存在很强的非双力偶成分，这可能是由多个不同震源机制的子事件叠加所致图本研究中使用的地震台站分布及其长周期波形记录图矩心矩张量反演过程与结果利用得到的矩心矩张量解计算的合成波形与观测波形的比较展示于图个台站观测波形和合成波形连接起来进行比较，发现者的均方根多点震源机制反演中年南斯科舍海岭地震的研究分析地质灾害论文这项研究中我们利用震相反演矩心矩张量式中，ξξ，ξ，ξξ，为格林函数场点与源点之间的路径效应，表示与时间无关的矩张量，为狄拉克函数，为归化震源时间函数形状通常设臵为角形可见，地震矩张量与地震位移之间具有线性关系，但是矩心时间和位臵与地震位移之间是非线性关系，因此，需要个迭代过程惯用的做法是，首先在给定矩心位臵的情况下搜索矩心时间，然后在给定矩心时间的情况下搜索矩心的水平位臵，最后再搜索矩心深度迭代过程反复进行，直到观测资料和合成资料之间的差别最小综上所述，单点的矩心矩张量反演可写成如下矢量方程为了获取年南斯科舍海岭地震的矩心矩张量，我们从全球台网以及数字地震台网联合会数据网站下载了震中距范围在的长周期地震记录和仪器响应数据，并根据背景噪声水平，以为方位间隔，挑选了个台图高质量的垂直向长周期记录图，沿断层走向的条直线上，通过破裂速度控制子事件位臵，在给定子事件的震源时间函数的情况下反演了年汶川地震震源机制的时空变化张勇等，但是，无论如何，还没有种通用的方法能够解决多点震源机制反演的问题图年南斯科舍海岭地震的构造与地震活动背景广义台阵反投影技术是近年来研究大地震破裂时空过程的种比较成熟的方法，它可以用来快速确定高频源的时空分布自年苏门答腊特大地震，至今，已衍生诸多新技术王晓欣等刘志鹏等，我们也曾多次在些地震中应用该技术研究其破裂时空特征杜海林杜海林等张旭等，视震源时间函数也常常用来分析大地震破裂的复杂性，它不但可以揭示地震矩或能量释放随时间的起伏变化，还可以揭示破裂传播的方向变化张旭张旭等，本文将介绍种新的多点矩张量反演的方法并应用于年南斯科舍海岭地震的多点矩张量解的反演首他类型这说明斯科舍板块与南极板块的交界处应力环境并不均匀，构造也不简单年发生的以来最大的次事件，其余震的空间分布以及和的震源机制结果表表明，这次地震断层近东西走向，断层向南倾斜，倾角在之间地震破裂长度超过，明显高于同级地震的平均长度其破裂持续时间也远大于基于标量地震矩估计的时间，年南斯科舍海岭地震矩心东侧约处曾于年月日发生过次具有类似震源机制的地震下称年南斯科舍海岭地震有学者使用震相种在上地幔与地表之间传播的类似于声学中回音的长周期波，反演了这两次地震的矩心矩张量，用反投影技术和面波视震源时间函数分析了破裂的传播特征，并利用远场体波面波以及可用的静态位移反演确定了它们的滑动分布，结果表明年南斯科舍海岭地震持续时间约，平均破裂速度，造成了两个滑动凹摘要年月日，在南极南奥克尼群岛北南极板块与斯科舍板块之间发生了次年南斯科舍海岭，我们利用全球分布的长周期和宽频带地震记录反演确定了这次地震随时间和空间变化的震源机制，验证了提出的种多点震源机制反演的新方法首先利用长周期记录的震相反演了这次地震的矩心矩张量解并利用体波提取了视震源时间函数，同时利用台阵反投影技术从宽频带记录中获得了这次地震的高频源的时空分布，然后基于矩心矩张量解视震源时间函数以及高频源的时空分布，实现了采用新方法对年南斯科舍海岭地震的多点震源机制反演矩心矩张量解表明，地震矩心在，矩心深度，半持续时间，释放标量地震矩，发震断层走向，倾角，滑动角视震源时间函数清楚地揭示了地震矩随时间变化的方位依赖性，总体上可以将时间过程分为前和后两个阶段，但前可细分为两次子事件根据台阵反投影结果，这次地震为沿海沟从西到东的单侧破裂，破裂长心矩张量解的反演，是观测波形资料约束下的矩张量解的逐级拆分第步，利用可以将地震整体事件作为点源的低频波形记录反演确定事件的矩心矩张量解矩心位臵，半持续时间和矩张量解，以备为后续多点源震源机制的反演提供模型数据第步，逐步增加波形资料中的高频成分，允许反映震源机制变化的信息包含于观测资料，同时也隐含着待求解的震源机制解不能远离矩心矩张量解的光滑约束引入视震源时间函数和高频源时空分布提供的信息，为待求震源机制解的矩心位臵和半持续时间提供可能的约束，以得到稳定且具有物理意义的反演结果需要强调的是，这种方法可应用于非线型甚至离散破裂的震源的变机制过程反演年南斯科舍海地震多点源震源机制解反演的成功为这种方法的有效性和可行性提供了很好的验证，随时空变化的震源机制结果表明，这次地震所处的构造背景和应力环境都比较复杂参考文献杜海林年苏门答腊安达曼大地震能量辐射源的时间域台阵技术分析硕士论文北京中国地震局地球物理研究所杜，破裂持续时间约高频源的时空变化有学者已经分别利用分布在澳洲的个宽频带台和北美拉丁美洲的个宽频带地震台的记录对这次地震的高频源进行了反投影成像，得到的两个结果均表明这是次向东传播的单侧破裂，破裂速度在之间，然而高频源时间函数明显不同，能量的空间分布也截然不同根据等，的结果，破裂持续时间为，破裂长度约个宽频带台站构成的台阵根据我们的计算分析，这些台站构成的台阵响应优于澳洲和拉丁美洲的台站构成的台阵图事实上，非洲台阵数据的反投影可以直接把高频源确定在余震分布的海沟带，而无需海沟带的强制约束，我们的反投影结果表明，这次地震持续时间，破裂尺度，优势破裂方位，平均破裂速度次子事件，能量峰值点分别出现在和，各子事件的平均破裂速度分为别和震源机制的时空变化截至目前，还没有关于这次地震的震源机制随时间或者空间变化的研究本研究借个矩心时间分别为和时，观测资料和合成资料的相关性达到最高然后，固定这个子事件的矩心时间，分别在时间范围和，内进行维搜索确定半持续时间结果表明，和分别为子事件和的最优半持续时间为了阐释多点源矩张量反演的必要性和合理性，我们利用单点源矩心和多点源矩张量反演结果分别计算的合成地震图与观测地震图进行了对比如图所示，矩心矩张量解能够使和，能够使频带范围内的相关系数和均方根误差达到和矩心矩张量解，这个数依次为和，个子事件的矩张量解能够使频带范围内的合成地震图和观测地震图的整体相关系数和均方根误差达到和图次子事件的震源机制与震源时间函数图观测资料与合成资料的比较讨论单点源震源机制年南斯科舍海岭地震已经有多个震源机制反演结果表，表这些结果表明这个地震的走向在至方位，倾角在与之间，滑动角在到空特征以及平均视震源时间函数时间特征，在反演震源机制的同时，选择性地确定矩心时间矩心位臵以及半持续时间综合考虑平均视震源时间函数与高频源的时空特征，将事件的空间位臵设臵为子事件和的峰值点的中间位臵，矩心时间定为与的均值将事件的位臵设为子事件和的峰值点的平均位臵，矩心时间为为了进步拆分子事件和，将资料的频率上限提高至拆分为子事件和根据高频源的分布将者的矩心位臵设定为与与结果的稳定性，反演过程需要在约束事件的震源机制不变的条件下进行这个约束类似于等，从原始记录中逐渐剥离子事件的技术事件和半持续时间通过维搜索来确定事件的搜索范围为事件的搜索范围为，反演结果如图所示，半持续时间分别为和，标量地震矩分别为和，对应的矩震级为和根据高频源的震源时间函数图以及视震源时多点震源机制反演中年南斯科舍海岭地震的研究分析地质灾害论文海林，许力生，陈运泰年汶川大地震的破裂过程地球物理学报，刘志鹏，宋超，盖增喜年级新西兰地震北京大学学报自然科学版，王晓欣，丁志峰，马延路年月日尼泊尔地震破裂能量中心运动轨迹地球物理学报，许力生，陈运泰年中国西藏玛尼地震的时空破裂过程地震学报，许力生，陈运泰，高孟潭月日印度古杰拉特破裂过程地震学报，张喆，许力生，杜海林高频辐射过程和发生原因地球物理