用CAP方法研究久治-玛曲MS4.1级地震震源机制
CAP方法反演震源机制的误差分析:以胶东半岛两次显著中等地震为例
CAP方法反演震源机制的误差分析：以胶东半岛两次显著中等地震为例郑建常;林眉;王鹏;徐长朋【摘要】利用区域波形数据使用CAP方法反演中强地震的震源机制正逐渐得到广泛应用.本文以胶东半岛近期发生的两次显著中等地震为例,讨论了使用CAP方法反演震源机制时的误差估计,展示了反演结果的不确定性分析过程.2013年l1月23日和2014年1月7日在山东莱州和乳山分别发生了M4.6和M4.3级中等地震,两次事件均造成了较大影响.我们基于CAP方法,使用自助抽样(bootstrap)技术多次重复反演过程,得到大样本量的震源机制解数据;基于这些数据,使用粒子群算法和聚类分析技术给出了优化解,估计了震源机制解的误差范围,并利用震源机制解的P、T轴给出了震源球上的概率密度分布.【期刊名称】《地球物理学报》【年(卷),期】2015(058)002【总页数】10页(P453-462)【关键词】莱州地震;乳山地震;波形反演;聚类分析;不确定性【作者】郑建常;林眉;王鹏;徐长朋【作者单位】山东省地震局,济南 250014;山东省地震局,济南 250014;山东省地震局,济南 250014;山东省地震局,济南 250014【正文语种】中文【中图分类】P3151 引言据山东台网测定，2013年11月23日13时44分在山东省莱州市（37.10°N，120.02°E）发生 M4.6级地震，这次地震是山东陆地地区自1995年苍山5.2级地震后发生的最大地震，影响范围广，山东东部市地普遍有感；2014年1月7日22时24分在山东乳山（36.80°N，121.70°E）发生 M4.3级地震，这次地震也造成胶东地区大面积有感.这两次事件是1970年以来胶东半岛陆地及近海地区发生的最强烈的地震活动，其中莱州地震震中区在1970年以来的小震目录上属于典型的少震、弱震区，活动水平不高，很少有ML≥3.0地震发生，仅在1991年2月以及2012年7月分别发生最大ML3.8级小震序列各一次；乳山地震震中区历史上曾发生公元1046年岠嵎山51／2级和1939年乳山下初51／2级地震.虽然这两次地震的震中区历史上没有强烈地震活动，但胶东半岛北部近海曾发生多次6、7级强震，如1548年渤海海峡7.0、1948年威海近海6.0以及1969年渤海7.4级等.因此确定这两次显著中等地震的震源机制对于研究区域地质构造的活动特征，以及研判该地区的地震危险性等具有重要的科学价值.我们使用近年来在国内得到广泛使用的CAP（Cut and Paste）方法反演这两次地震的震源机制.为了得到更准确的解，并且合理地估计反演结果的不确定性，我们使用自助抽样法（bootstrap）对反演过程随机重复，在大样本量的反演结果基础上，使用粒子群算法搜索优化解，利用动态聚类技术对结果进行聚类分析，从而得到了更加稳定可靠的断层面解，给出了可能的误差范围，并进一步给出了震源球上P、T轴的概率密度分布.2 理论与方法2.1 CAP方法反演震源机制震源机制和传播效应决定了观测波形的变化.如果地壳模型已知，可以准确地计算波形传播过程中的效应，因此我们可以通过理论波形s（t）和观测波形u（t）的拟合来估计震源的断层面参数.双力偶震源产生的理论位移s（t）可以表示为（Zhu and Helmberger，1996）：其中，i＝1，2，3对应三种基本断层响应，即：垂直走滑、垂直倾滑以及倾角为45°的倾滑；Gi为格林函数，Ai是辐射系数，φ是台站方位角，M0为标量地震矩.θ，δ，λ分别为断层的走向、倾角、滑动角.系数Ai由6个矩张量分量和台站方位角表示如下（Jost and Herrmann，1989）：走向θ、倾角δ、滑动角λ，以及标量地震矩M0等可以通过求解以下方程进行估计：波形反演可以使用全波形数据，也可以单独使用体波或面波震相进行拟合.CAP方法是一种联合使用体波和面波进行反演的方法，近年来在国内得到了广泛的应用（吕坚等，2008；黄建平等，2009；郑勇等，2009；龙锋等，2010；韩立波等，2012），由于该方法分别截取波形的Pnl部分和面波部分分别拟合（Zhao and Helmberger，1994；Zhu and Helmberger，1996），并在反演的过程中允许它们在适当的时间变化范围内相对移动，在一定程度上避免了因为地壳模型不准确而引起的震相到时的误差因素，对速度模型和地壳横向变化的依赖性较小，因此在实际的区域地震震源机制求解中有明显的优势.CAP方法使用频率F－波数K法（Zhu and Rivera，2002）计算格林函数，使用网格搜索方法搜寻最优震源机制参数和震源深度.考虑到波形随震中距的衰减，方法定义误差函数如下：式中，r为台站震中距，r0为选定的参考震中距，p为指数因子.参考有关研究，对体波p＝1，面波p＝0.5（韩立波和蒋长胜，2012）.2.2 CAP方法的优化解及其不确定性估计地球内部的任意震源可以表示为6个独立分量的矩张量，由于CAP方法限制震源为双力偶模型，并且无需发震时刻的对齐，因此只需对震源模型的三个角度，即走向θ、倾角δ、滑动角λ，以及标量地震矩M0进行搜索，理论上而言，仅需要2个台站的三分向波形就可以求解；虽然研究显示，对于大多数3个三分向台的组合，使用波形反演就可以得到相对准确的震源机制，但实际情况也显示，不同的台站组合波形反演得到的解之间仍然存在一定的差异（Godano et al.，2009；郑建常和陈运泰，2012）.目前国内台网密度已经达到相当水平，在东部地区，一个中等地震通常有数十甚至上百个宽频带台能够记录到清晰的波形，以此次莱州地震为例，通过对原始波形进行去均值、去趋势、积分等简单变换后，根据直观的观察，震中距300km范围内，采样率100Hz的宽频带三分向波形有近40个台站的资料可用.在使用CAP方法求解震源机制时，一般的做法是选择部分波形拟合较好的台进行反演，删去拟合不好的台或者震相；有些情况下，甚至仅使用面波部分而删除体波震相，需知面波尤其是径向和切向分量，很容易受到台站下方浅层地壳结构的影响.由于CAP方法是采用网格搜索的方法，因而这种人为的选择，必然会为反演结果增加主观的不确定性因素.我们无法令人信服地说明，拟合不好的波形究竟是数据本身确实存在干扰，还是说搜索到的解无法满足该条数据.另外，CAP方法虽然可以在最后的输出结果中给出断层面参数的不确定性，但该估计值只是面向所使用的台站数据的结果，在上述的人为选择下，该不确定性估计能够在多大程度上客观地反映最终解的整体不确定性，是无法说明的.为了求得更加稳定可靠的解并且合理客观地给出解的误差估计，我们在相对丰富的观测数据基础上，采用自助抽样统计方法进行分析.具体方法是在可用的观测台站中可重复地随机抽取一定数量的台站组成新的台站组合，使用该台站组合的观测数据重复反演过程.在大量的重复计算后（例如，超过1000次），可以有效地排除观测质量不高或存在较大干扰误差的数据的影响，从而得到更加接近真实解的结果，并且可以有效地给出解的不确定性.另外，由于CAP方法在搜索断层面解时采用的是网格搜索的方法，然后通过插值计算误差函数e的最小值，并且由于固定步长的走向、倾角、滑动角的尝试位置在震源球上的分布是不均匀的（许向彤等，1995），因此在最终解中可能会有空缺（gap）的存在.为了求解优化解，我们进一步使用Kagan（1991）定义的双力偶模型最小空间旋转角，对上面自助抽样得到的大量满足条件的震源机制结果进行分析，定义与所有解的空间偏转角度和为目标函数，使用粒子群非线性优化方法搜索该目标函数最小的结果，视为最优解.2.3 聚类分析在震源机制求解中，常见的情况是在震源球上存在几簇相对集中分布的解，对这些可能的解直接取数学平均是不甚合理的，并且在数据存在较大误差或干扰的情况下，满足条件的可能解的分布范围也许会相当大.因此针对这一现象，刁桂苓等（1992）、俞春泉等（2009）分别使用系统聚类和动态聚类技术，对所有的可能解进行聚类分析，求取聚类中心作为反演的优化解，数值试验和实际应用都有很好的效果.聚类分析可以很好地排除孤立解和错误解，从而在大量的数据中获取更加接近真实解的结果.本文在使用不同台站组合重复进行波形反演后，同样得到了大量的震源机制解数据，受台站布局和数据误差的影响，这些解或多或少存在差别，因此对这些结果进行聚类分析是很有必要的.3 数据与资料本文使用了山东台网提供的波形资料，其中还包括了邻省如辽宁、河北、江苏等省交换资料的部分台站.图1给出了本项研究使用的台站分布，其中个别台如JIM、ZSL、HUD等为短周期台，在波形反演中没有使用.本文研究中，首先由观测记录直接读取初动符号，用于约束波形反演；然后将观测数据扣除仪器响应，经过去均值、零漂等预处理后积分至位移记录，旋转到Z－R－T坐标系，对观测波形和理论波形同样进行带通滤波，然后用于反演.图1 本文研究的两次地震震中及山东台网台站分布图Fig.1 Map of stations in Shandong Network and two earthquakes studied in this paper.Red circles denote epicenter，triangles are stations，and solid black lines are faults.使用Chang等（2006）给出的朝鲜半岛南部至黄海地区的中上地壳速度结构模型用于本文的震源机制反演.相关地质资料显示，胶东半岛、南黄海以及朝鲜半岛南部在大地构造分区上都属于下扬子地块，地质构造属性相对较为一致（Ree et al.，1996）.4 结果与分析4.1 乳山M4.3震源机制选择震中距在250km以内的15个台站的宽频带波形记录进行反演，Pnl和面波的反演波段分别选择0.05～0.15Hz和0.033～0.067Hz频段.图2给出了不同深度的最佳双力偶解，及拟合误差随不同深度变化的关系，由图可见，震源深度在4km时观测波形和理论波形的错配值最小，说明事件的震源深度较浅.由CAP方法反演得到的最佳震源机制：节面A的参数为：走向202°、倾角75°、滑动角153°；节面B的参数为：走向299.5°、倾角64°、滑动角16.7°；参考乳山序列的双差定位结果（李冬梅和郑建常，2014）分析认为，节面B可能是乳山地震的发震断层；震源机制显示为左旋走滑型，反演得到此次地震的矩震级MW＝4.2.图2 2014年1月7日乳山M4.3级地震不同震源深度的波形拟合误差及最佳震源机制解Fig.2 Waveform fit errors and best focal mechanisms as function of depth for Jan.7，2014Rushan M4.3event图3给出了对应最优解的理论波形和观测波形的拟合情况.15个台一共75个震相，其中理论波形与观测波形相关系数大于0.9的有39个，超过50%；相关系数大于0.6（相关性较好）的有67个，约占89.3%；最佳解的方差减少（variance reduction）为70.3%，说明理论波形很好地拟合了观测波形，反演结果是可靠的.个别台（如WEH）平均相关系数较差，可能与台站位于震源机制解的节面线附近，振幅相对较小所致；另外如CHD台的拟合程度不好，可能与该台处于海域、噪声干扰较大有关.使用自助抽样的统计方法，对乳山地震震源机制解的不确定性进行估计.选用震中距在300km以内的采样率为100Hz的22个三分向宽频带台的观测波形组成原始数据集，为了保证用于反演的数据的样本量，设用于反演的台站数为20个，对原始数据集进行每个台站等概率、可重复地随机抽取，抽取出的台站波形组成新的数据集，然后用于CAP方法的波形反演.对上述的抽取台站反演过程重复1000次，将反演得到的震源机制的断层节面解和P、T轴绘制在一个震源球上，见图4.可以看出，反演中除去个别反演过程的断层面解出现一定程度的偏离外，其余结果集中分布，均显示为近走滑的机制；图4中的P、T轴位置和断层节面线集中成丛，大致显示出断层面解的误差范围.使用粒子群非线性优化方法，以与自助抽样给出的1000个机制解（图4）的Kagan角之和为目标函数，搜索最优解.结果显示最优解为，节面A：走向208.4°，倾角89.7°，滑动角154.3°；节面B：走向298.5°，倾角64.3°，滑动角0.3°；最优解与图4所示1000个解的平均夹角4.37°，以其与所有解Kagan角的2倍标准差为震源机制解的误差范围，结果显示不确定性为6.44°（图5）.表1 波形反演乳山M4.3地震震源机制解结果Table 1 Parameters of focal mechanism results from waveform inversion for Jan.7，2014RushanM4.3event节面A方法节面B走向（°）倾角（°）滑动角（°）走向（°）倾角（°）滑动角（°）CAP 202 75 153 299.5 64.0 16.7粒子群优化 208.4 89.7 154.3 298.5 64.3 0.3聚类分析208.0 89.3 154.1 298.4 64.1 1.8对自助抽样结果进行动态聚类分析，结果显示最优解为：节面 A：走向208.0°，倾角89.3°，滑动角154.1°，节面B：走向298.4°，倾角64.1°，滑动角1.8°，与粒子群优化解非常一致（见表1）.将自助抽样结果中的P、T轴投影到震源球上（图4），对其进行概率密度统计分析，结果见图6.4.2 莱州M4.6震源机制使用CAP方法对2013年11月23日莱州M4.6地震进行反演（郑建常等，2015），同样进行CAP反演情况的自助抽样统计分析.选用震中距在270km以内的采样率为100Hz的22个三分向宽频带台的观测波形组成原始数据集，采用全样本随机抽取方法，自助抽样反演1000次，图7给出了反演得到的震源机制的断层节面解和P、T轴在震源球上的分布情况.结果显示，莱州地震的自助抽样结果同样很好地显示出了反演得到震源机制解的误差范围，相对于乳山地震，出现了极个别反演过程的结果偏离较大的情况.使用粒子群非线性优化方法，搜索与自助抽样结果的旋转角最小的解.结果显示最优解为，节面A：走向236.9°，倾角76.2°，滑动角－169.3°；节面B：走向144.3°，倾角79.6°，滑动角－14.0°；与所有自助抽样解的平均偏转角17.4°（图8），以其与所有解Kagan角的2倍标准差（图8红色虚线所示）为误差范围，结果显示震源机制解的不确定性为23.7°.从自助抽样得到的所有机制解在震源球上的分布情况（图7）可以直观地看出，断层节面线尤其是北西向节面呈现出两组集中.由于我们定义的粒子群优化的目标函数是搜索与所有自助抽样解的空间旋转最小，因此从图7可以看出，最优解的节面位置处于其中一组的边缘位置，在此情况下，对自助抽样结果进行聚类分析是有意义的.图3 2014年1月7日乳山M.3地震最优解的理论波形（红）与观测黑波形图下方第一行数字为各段理论地震波形相对实际观测波形的移动时间，正值表示理论波形相对观测波形超前.第二行数字为理论波形与观测波形的相关系数（百分比）.波形图左侧字母为台站，其下数字分别为台站震中距（km）和方位角（°）.图左侧的震源球上红色区域代表压缩区，白色代表拉张区，震源球采用下半球投影.震源球上标注的“＋”和“－”表示反演使用台站的P波初动.Fig.3 Comparison between synthetics（red）and observed（black）seismograms of Jan.7，2014Rushan M4.3event The numbers on the lower left side of each seismogram are the time shifts（upper）and cross－correlation coefficient in percent（lower）.Positive time shifts mean that the observed data have been delayed.The letters on the left side are stations，the numbers below it are epicentral distance（in km）and azimuth（degree）.The red color in beach－ball denotes compression area，while white isextension.The‘＋’and‘－’signs on beach－ball indicate polarities on inversion used stations.Lower hemisphere projection is used.图4 自助抽样得到的1000次乳山地震震源机制解及粒子群最优解（下半球极射投影）震源球上黑色细线条表示自助抽样结果的断层节面线，红色线条表示粒子群优化解的节面线.Fig.4 1000focal mechanisms of Rushan M4.3event retrieved by a bootstrapping process，all nodal lines（black）and P，Taxes（blue and red points，respectively）are plotted on one beach－ball.The red lines on the beach－ball show the optimized solution given by a Particle Swarm Optimization method.Its corresponding P，Taxes are also displayed on the beach－ball （yellow and green point，respectively）.Lower hemisphereprojection is used.图5 乳山地震自助抽样结果与粒子群优化解的Kagan角分布Fig.5 Kagan angles of bootstrap results to the PSO solution for Rushan event以两个震源机制解之间的最小空间旋转角（即Kagan角）为距离的定义，对1000次自助抽样结果进行聚类分析，图9给出了聚类谱系图（由于完整的聚类树过于密集和庞大，因此我们只显示了Kagan＞7°的部分），以50°为阈值，可以将结果分为5类.图10给出了聚类分析的结果，属于Ⅰ类的数据占32.1%，Ⅱ类66.6%，其余三类数据合计仅有1.3%.由图10可以看出，其余三类的断层节面线和P、T轴位置明显偏离集中区且机制解类型与绝大部分结果（走滑型）不一致，是典型的孤立解.孤立解（或错误解）的出现，可能说明我们使用的数据中个别台站（或分向）存在较大干扰.使用俞春泉等（2009）的方法求取了四类解的聚类中心，其中I类解的聚类中心的断层面参数（设北东向节面为发震断层面）为：走向231.6°，倾角88.7°，滑动角－168.2°；Ⅱ类解的聚类中心为：走向238.5°，倾角74.1°，滑动角－164.8°.图6 乳山地震自助抽样结果的震源球概率密度分布俯视图（未进行极射投影）色标中正值表示T轴的概率密度分布，负值表示P轴的概率密度分布.Fig.6 Probability density distribution of solutions on beach－ball（top view of lower hemisphere，without projection）Positive values on the color scale （corresponding to the red area on beach－ball）indicate probability of T axis，while negative values（corresponding to blue area）mean probability of Paxis.图7 自助抽样得到的莱州地震震源机制解及粒子群最优解（下半球极射投影）黑色节面线为自助抽样得到的震源机制解；红色节面线为粒子群最优解.Fig.7 1000Focal mechanisms of Laizhou M4.6event retrieved by a bootstrappingprocess，all nodal lines（black）and P，Taxes（blue and red points，respectively）are plotted on one beach－ball.The red lines on the beachball show the optimized solution given by a PSO method.Its corresponding P，Taxes are also displayed on the beach－ball（cyan and yellow point，respectively）.Lower hemisphere projection is used.图8 莱州地震自助抽样结果与粒子群优化解的Kagan角分布Fig.8 Kagan angles of bootstrap results to the PSO solution for Laizhou event图9 莱州地震自助抽样结果的聚类谱系图Fig.9 Dendrogram plot of the hierarchical binary cluster tree for Laizhou event图10 莱州地震自助抽样结果的聚类分析（a）断层节面线的分类显示；（b）机制解P（＋）、T（⊙）轴位置的分类显示.断层节面线和P、T轴颜色表示分类，与图9分类颜色一致.Fig.10 Clustering results of focal mechanisms from a bootstrap process for Laizhou event The different colors of nodal lines and P （“＋”sign in right panel），T （“⊙”in right panel）axes denote different classes，which are corresponding to colors shown in Fig.9将自助抽样结果中的P、T轴投影到震源球上，对其进行概率密度统计分析，结果见图11.可见，同聚类分析的结果一致，P轴位置的概率密度在震源球上出现了两个极值区，分别对应Ⅰ类和Ⅱ类两个聚类中心.5 讨论与结论基于山东省宽频带数字地震波形资料，本文首先使用CAP方法反演了近期胶东半岛地区发生的两次显著中等地震活动的震源机制，讨论了如何合理地估计CAP方法反演震源机制的误差范围以及如何确定优化解的问题.我们首先使用自助抽样方法，对原始数据进行等概率随机抽样，多次重复波形反演过程，排除了人为选择数据的干扰，得到大样本量的震源机制解数据；在此基础上，我们（1）使用了粒子群优化算法从中搜索与这些机制解空间偏转角最小的解当作优化解，以Kagan角的二倍标准差作为反演结果的不确定性范围，结果显示：乳山地震的粒子群优化解为：走向298.5°，倾角64.3°，滑动角0.3°，不确定性为±6.4°；莱州地震的优化解为：走向236.9°、倾角76.2°、滑动角－169.3°，不确定性为±23.7°.（2）对自助抽样结果进行聚类分析，其中：乳山地震结果的聚类中心与粒子群优化解基本一致；莱州地震结果存在多个聚类，排除孤立解后，有两个聚类中心，其对应两类数据合计占结果的98.7%，说明此次地震的真实解在这两类数据范围内. 图11 莱州地震自助抽样结果的震源球概率密度分布（a）下半球俯视图（未进行极射投影）；（b）东南方向，45°三维侧视图.色标说明同图6.Fig.11 Probability density distribution of solutions on beach－ball（a）Top view of lower hemisphere，without any projection；（b）Side view of whole beach－ball，from south－east 45°po sition.Positive values on the color scale （corresponding to the red area on beach－ball）indicate probability of Taxis，while negative values（corresponding to blue area）mean probability of Paxis.（3）将自助抽样结果中的P、T轴投影到震源球，对其进行概率密度统计，给出了机制解在震源球上的概率密度分布图.本文方法不单可以得到更准确的震源机制优化解、给出科学合理的误差估计，而且可以有效地排除孤立解和错误解，克服存在较大干扰台站的数据的影响，因此在震后应急的震源机制自动化求解中可以发挥作用.在强大计算能力的支持下，无须人工干预即可得到准确可靠的震源机制结果，从而为震害评估、趋势分析等提供重要的科学依据.自助抽样结果显示乳山地震震源机制解的不确定性要小于莱州地震，笔者推测可能有莱州地震使用的台站中个别台的干扰较大的原因，另外也无法排除莱州地震的震源破裂过程可能更加复杂的可能.CAP方法中用于计算格林函数的F－K方法使用狄拉克－Delta函数作为震源时间函数（Zhu and Rivera，2002），对于小震级的事件该简化方案更为适用，莱州地震（M4.6）与乳山地震（M4.3）震级相差不大，但莱州地震的自助抽样结果出现了两个概率较高的聚类中心，在使用大部分相同台站的情况下，这可能意味着莱州地震的震源破裂随时间的变化可能与狄拉克－Delta函数存在一定的偏离.Rodríguez－Lozoya等（2008）的研究显示，区域中等地震也可能有复杂的震源破裂过程，在该问题上的深入研究需要更进一步的工作. 致谢感谢两位匿名审稿专家提出的宝贵意见.聚类分析中使用了俞春泉等（2009）提供的部分开放代码，粒子群搜索使用了S Chen给出的Matlab软件包，在此一并表示感谢！ReferencesChang S J，Baag C E.2006.Crustal structure in Southern Korea from joint analysis of regional broadband waveforms and traveltimes.Bull.Seism.Soc.Amer.，96（3）：856－870.Diao G L，Yu L M，Li Q Z.1992.Hierarchical clustering analysis of the focal mechanism solution－Taking the Haicheng Earthquake Sequences for example.Earthquake Research in China （in Chinese），8（3）：86－92. Godano M，Regnier M，Deschamps A，et al.2009.Focal mechanisms from sparse observations by nonlinear inversion of amplitudes：method and tests on synthetic and real data.Bull.Seism.Soc.Amer.，99（4）：2243－2264.Han L B，Jiang C S，Bao F.2010.Source parameter determination of2010Taikang MS4.6earthquake sequences.Chinese J.Geophys.（inChinese），55（9）：2973－2981，doi：10.6038／j.issn.0001－5733.2012.09.016.Han L B，Jiang C S.2012.Focal mechanism inversion of 8Jun 2011 Toksun MS5.3earthquake.Acta Seismologica Sinica （in Chinese），34（3）：415－422.Huang J P，Ni S D，Fu R S，et al.2009.Source mechanism of the 2006MW5.1Wen′an earthquake determined from a joint inversion of local and teleseismic broadband waveform data.Chinese J.Geophys.（in Chinese），52（1）：120－130.Jost M L，Hermann R B.1989.A student′s guide to and review of moment tensors.Seism.Res.Lett.，60（2）：37－57.Kagan Y Y.1991.3－D rotation of double－couple earthquakesources.Geophys.J.Int.，106（3）：709－716.Li D M，Zheng J C.2014.Relocation analysis of Oct.1，2013 Rushan swarms.Qilu Earthquake Sciences （in Chinese），10（1）：26－28.Long F，Zhang Y J，Wen X Z，et al.2010.Focal mechanism solutions ofML≥4.0events in the MS6.1Panzhihua－Huili earthquake sequence of Aug 30，2008.Chinese J.Geophys.（in Chinese），53（12）：2852－2860，doi：10.3969／j.issn.0001－5733.2010.12.008.LüJ，Zheng Y，Ni S D，et al.2008.Focal mechanisms and seismogenic structures of the MS5.7and MS4.8Jiujiang－Ruichang earthquakes of Nov.26，2005.Chinese J.Geophys.（in Chinese），51（1）：158－164，doi：10.3321／j.issn：0001－5733.2008.01.020.Ree J－H，Cho M，Kwon S－T，et al.1996.Possible eastward extension ofChinese collision belt in South Korea：The Imjingang belt.Geology，24（12）：1071－1074.Rodríguez－Lozoya H E，Quintanar L，Ortega R，et al.2008.Rupture process of four medium－sized earthquakes that occurred in the Gulf of California.J.Geophys.Res.，113（B10301），doi：10.1029／2007JB005323. Xu X T，Xu Z H，Zhang D N.1995.A probabilistic grid test method of determining earthquake focal mechanisms using P－wave onset polarity data.Seismological and Geomagnetic Observation and Research （in Chinese），16（4）：34－42.Yu C Q，Tao K，Cui X F，et al.2009.P－wave first－motion focal mechanism solutions and their quality evaluation.Chinese J.Geophys.（in Chinese），52（5）：1402－1411，doi：10.3969／j.issn.0001－5733.2009.05.030.Zhao L S，Helmberger D V.1994.Source estimation from broadband regional seismograms.Bull.Seism.Soc.Amer.，84（1）：91－104.Zheng J C，Li D M，Wang P，et al.2015.Focal mechanisms and seismic tectonic features of the 2013Laizhou M4.6earthquakesequence.Seismology and Geology （in Chinese），in press.Zheng J C，Chen Y T.2012.Stability of sparse station data inversion for deviatoric moment tensor solution of regional earthquakes.Acta Seismologica Sinica （in Chinese），34（1）：31－43.Zheng Y，Ma H S，LüJ，et al.2009.Source mechanism of strong aftershocks（Ms≥5.6）of the 2008／05／12Wenchuan earthquake and the implication for seismotectonics.ScienceinChinaSeriesD：Earth Sciences，52（6）：739－753.Zhu L P，Helmberger D V.1996.Advancement in source estimation techniques using broadband regional seismograms.Bull.Seism.Soc.Amer.，86（5）：1634－1641.Zhu L P，Rivera L A.2002.A note on the dynamic and static displacements from a point source in multilayered media.Geophys.J.Int.，148（3）：619－627.附中文参考文献刁桂苓，于利民，李钦祖.1992.震源机制解的系统聚类分析——以海城地震序列为例.中国地震，8（3）：86－92.韩立波，蒋长胜，包丰.2012.2010年河南太康MS4.6地震序列震源参数的精确确定.地球物理学报，55（9）：2973－2981，doi：10.6038／j.issn.0001－5733.2012.09.016.韩立波，蒋长胜.2012.2011年6月8日新疆托克逊MS5.3地震震源机制解反演.地震学报，34（3）：415－422.黄建平，倪四道，傅容珊等.2009.综合近震及远震波形反演2006文安地震（MW5.1）的震源机制解.地球物理学报，52（1）：120－130.李冬梅，郑建常.2014.2013年10月1日乳山震群重新定位结果分析.齐鲁地震科学，10（1）：26－28.龙锋，张永久，闻学泽等.2010.2008年8月30日攀枝花—会理6.1级地震序列ML≥4.0事件的震源机制解.地球物理学报，53（12）：2852－2860，doi：10.3969／j.issn.0001－5733.2010.12.008.吕坚，郑勇，倪四道等.2008.2005年11月26日九江－瑞昌MS5.7，MS4.8地震的震源机制解与发震构造研究.地球物理学报，51（1）：158－164，doi：。
2017年6月16日湖北秭归M_S4.3地震成因初探
为更好地认识本次地震的发震 背 景，本 文 对 其 震 源 机 制 解 及 其 震 源 深 度 进 行 了 研 究。 首先获取了近震宽频带固定台站的三分量数据 ，基于两种１－Ｄ 速度模型计算不同深度下的 格林函数，然后使用 ＣＡＰ（Ｃｕｔ Ａｎｄ Ｐａｓｔｅ）反演得到本次地震的震源机制和震源深度 。采用 基阶 Ｒａｙｌｅｉｇｈ波振幅谱进一步约束质心深度，最后将本次地震与巴东地震、２０１４年秭归地 震 以 及 ２０１１ 年 阳 新 地 震 进 行 对 比 。结 果 表 明 ，本 次 地 震 可 能 是 水 库 蓄 水 诱 发 的 地 震 ，但 要 确认其成因还需进一步研究。
＊ 收稿日期：２０１８－０７－２３；修改回日期：２０１９－０１－１６ 基 金 项 目 ：国 家 自 然 科 学 基 金 （４１６６１１６４０３５）资 助 作 者 简 介 ：李 伟 （１９９４－），男 ，四 川 宜 宾 人 ，在 读 硕 士 研 究 生 ，主 要 从 事 地 震 学 研 究 。 通 讯 作 者 ：储 日 升 ，研 究 员 。Ｅ－ｍａｉｌ：ｃｈｕｒ＠ａｓｃｈ．ｗｈｉｇｇ．ａｃ．ｃｎ
阿克塞MS4.5地震震源机制及其指示意义
Ⅱ-91阿克塞M S 4.5地震震源机制及其指示意义徐 溶※ 王朋涛 王丽霞（甘肃省地震局，兰州 730000）中图分类号：P315.3+3 文献标识码： A doi ：10.3969/j.issn.0253-4975.2018.08.0552018年6月17日11时12分，甘肃酒泉市阿克塞县发生M 4.5地震，据中国地震台网中心测定的结果，震中位置为（38.87°N ，94.93°E ），震源深度15 km 。
阿克塞M 4.5地震是一次发生在祁连山地震带内的显著地震，同时，也是在祁连山西段M L 4.0地震长期平静背景下发生的显著地震事件。
有研究表明祁连山地震带显著性地震事件的出现与区域地震活动增强和成组前兆同步可以作为该区中强地震的预测指标。
基于甘肃及邻区地震台网记录的波形资料，采用CAP （Cut and Paste ）方法反演阿克塞4.5级地震震源机制解。
受本次地震发生所在区域监测条件的限制，其台站分布较少且方位角不均匀。
而CAP 方法的优势在于综合地震波形反演对震源机制结果的多方面影响因素，即使在台站分布不均匀的情况下也可得到稳定的可靠结果。
地震后共挑选了信噪比较高的7个台站得到此次地震震源机制，获取震源矩心深度为6 km ，最佳双力偶解：节面Ⅰ走向332°、倾角19°、滑动角97°，节面Ⅱ走向145°、倾角71°、滑动角88°，矩震级M W =4.52。
节面Ⅰ走向为NNW —SSE 向及逆冲型的破裂类型，与其附近的党河南山活动断裂的性质基本一致，分析认为，本次地震的发生与该断裂的活动有关，判定节面Ⅰ代表了本次的发震断层面。
祁连山地震带自2010年以来超过300天以上的4.0级平静共有2次，第1次363天的地震平静被2011年12月11日甘肃酒泉阿克塞M L 4.2地震打破，4个月后2012年5月3日发生了金塔5.4级地震；第2次301天的地震平静被2015年8月2日甘肃林泽M L 4.1地震打破，3个月后2015年11月23日发生青海乌兰5.2级地震。
2017年8月8日四川九寨沟M7.0地震及其余震序列的震源参数
①ｈｔｔｐ：／／ｗｗｗ．ｃｅｎｃ．ａｃ．ｃｎ／ｃｅｎｃ／ｒｅｓｏｕｒｃｅ／ｃｍｓ／２０１８／０６／２０１８０６０８１１４２１５４５６４１．ｐｄｆ。
②ｈｔｔｐ：／／ｗｗｗ．ｃｅａ．ｇｏｖ．ｃｎ／ｃｅａ／ｄｚｐｄ／ｄｚｚｔ／３６９８６１／３６９８６２／３５８３３２０／ｉｎｄｅｘ．ｈｔｍｌ。
第４２卷 第２期２０２０年４月地 震 地 质ＳＥＩＳＭＯＬＯＧＹＡＮＤＧＥＯＬＯＧＹＶｏｌ．４２，Ｎｏ．２Ａｐｒ．，２０２０ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．０２５３－４９６７．２０２０．０２．０１５２０１７年８月８日四川九寨沟Ｍ７ ０地震及其余震序列的震源参数吴微微 魏娅玲 龙 锋 梁明剑陈学芬 孙 玮 赵 晶（四川省地震局，成都 ６１００４１）摘 要 ２０１７年８月８日２１时１９分４６秒，四川北部阿坝州九寨沟县发生Ｍ７ ０地震。
该地震发生在东昆仑断裂带东段南端的一条分支断裂上，地处玛曲—玛沁地震空区和松潘—平武地震破裂区的过渡部位。
为了解这次地震的震源参数特征以及震区构造环境，文中基于四川数字波形资料及中国地震台网中心的震相数据，采用Ｌｏｃ３Ｄ定位方法获取初始地震目录参数，并利用双差定位方法重新定位地震序列中的较小地震，选用ＣＡＰ（Ｃｕｔ ａｎｄ Ｐａｓｔｅ）方法分析２９次中强地震的震源机制解与震源矩心深度，同时对１８６次地震（２ ０≤ＭＬ≤５ ５）进行了震源谱恢复，绘制出余震序列的震源应力降空间分布图像。
重新定位结果显示，该地震序列以Ｍ７ ０主震震中为界表现出明显的空间分段差异：余震区的ＮＷ段地震数量较少、震源应力降较大、震源机制解类型复杂；ＳＥ段分布着多数４ ０级以上中强地震，且几乎全部为纯走滑型地震，余震序列的应力降随时间逐渐减小。
结合震源机制解的研究结果分析，发震构造可能为ＮＷ－ＳＥ向、断续的隐伏断裂或虎牙断裂的北延推测段，断层面在主震南、北两侧发生明显转折，同时由挤压逆冲为主逐渐转换为以走滑运动为主。
祁连山中东段地区中小地震震源机制和区域应力场
祁连山中东段地区中小地震震源机制和区域应力场荣代潞【期刊名称】《地震地磁观测与研究》【年(卷),期】2014(35)3【摘要】利用哥伦比亚大学 GCMT 目录给出的祁连山中东段地区中强地震震源机制资料，研究较大区域（34°-41°N，100°-106°E）的应力场；利用该地区布设的中法微震数字监测台网多年监测资料和甘肃数字监测台网资料，使用 P 波和 S 波初动及振幅比联合反演方法，反演中小地震震源机制解和发震应力场。
结果表明，地区构造应力大致为北东40°-45°水平向压应力；景泰地区主压应力方向约北东45°，绝大多数地震为走滑型。
天祝-古浪地区有相当部分的逆断层地震分布，主压应力方向约60°，P 轴仰角在10°左右优势分布，大致为水平应力场。
这与大区域构造应力场和断层实际分布基本一致。
%The stress field in larger area (from 34°N,100°E to 41°N,106°E)of Qilian Mountain area have been studied using the data of mechanisms of moderate-strong earthquakes in the area from the GCMT catalogs of Columbia University.The focal mechanism solutions and the seismogenic stress field have been inverted for the small-moderate events based on the data from the networks set by China and France and the Gansu Digital Monitoring Network using the joint inversion method of first motion and amplitude ratio of P and S wave. The results shoe that the tectonic stress of this area is a horizontal compressive stress with north-east di-rection of 40°-45°. The direction of main pressure stress in the area of Jingtai is about north-east with 45°,the type of mostearthquakes are strike. There are quite a part earth-quakes of reverse fault distributed in Tianzhu-Gulang area with horizontal stress of main compressive with about 60°and advantage elevation angle about 10°. It is accordance with the tectonic stress and fault distribution in larger area.【总页数】7页(P1-7)【作者】荣代潞【作者单位】中国甘肃 730000 中国地震局兰州地震研究所【正文语种】中文【相关文献】1.北天山中东段中小地震震源机制解及应力场反演 [J], 龙海英;高国英;聂晓红;李莹甄;沈军2.祁连山中东段基于CAP方法的中小地震震源机制研究 [J], 张辉;徐辉;姚军;代炜3.川滇地区中小地震震源机制解及构造应力场的研究 [J], 郭祥云;陈学忠;王生文;王恒信4.川滇地区中小地震震源机制解及构造应力场的研究 [J], 郭祥云;陈学忠;王生文;王恒信;5.张家口地区近年中小地震震源机制解及应力场反演 [J], 王曰风;张秀萍;李峰;马利军;宋晓煜;刘燕翔;张珊珊因版权原因，仅展示原文概要，查看原文内容请购买。
2018年佳木斯汤原ML 4.4地震震源机制及地震趋势分析
2018年佳木斯汤原ML 4.4地震震源机制及地震趋势分析Chang Jinlong;Gan Weijun;Jiang Bo;Hu Baohui
【期刊名称】《地震地磁观测与研究》
【年(卷),期】2018(039)006
【摘要】2018年1月15日佳木斯汤原县发生ML4.4地震,采用ISOLA方法反演矩张量解,结果显示为逆冲型地震,兼少量走滑成分.结合震源区地质构造背景,推断依舒断裂北段通河—汤原段西支断裂为发震构造,主压应力轴与背景应力场方向不一致.分析认为,黑龙江亚板块由东南转向东北运动,造成断层两侧应力产生差异,区域应力场调整诱发汤原ML 4.4地震.
【总页数】8页(P70-77)
【作者】Chang Jinlong;Gan Weijun;Jiang Bo;Hu Baohui
【作者单位】
【正文语种】中文
【相关文献】
1.2016年10月20日射阳Ms4.4地震的震源机制与地震序列的时空分布特征 [J], 李锋;李迎春;李强;冯志生;王俊菲
2.采用CAP方法反演2010年贵州罗甸4.4级地震震源机制解 [J], 郭培兰;谢夜玉;龙政强;向巍
3.龙羊峡震群及
4.4级主震格尔木东南
5.4级地震的震源机制 [J], 都昌庭
4.玛曲MS4.4地震震源机制解及其构造意义探讨 [J], 陈晓龙;蒲举;尹欣欣;陈继锋
5.玛曲M_(S)4.4地震震源机制解及其构造意义探讨 [J], 陈晓龙;蒲举;尹欣欣;陈继锋
因版权原因，仅展示原文概要，查看原文内容请购买。
文秘知识-2019年洱源Ms5.5地震序列ML≥3.4地震的震源机制解研究 精品
2.1 CAP方法基本原理
近年来，CAP方法在求解震源机制解的方法中应用较多，CAP方法（Cut and Paste）相比以往的P波初动、体波反演或是面波反演，算得上是一种全波形反演方法，它将宽频带地震记录分成Pnl和面波两个部分进行反演并允许他们相对浮动，在适当的时间变化范围内，搜索出合成地震图和观测地震图全局差异最小的震源机制解。CAP方法具有对参与反演的台站数要求不多、方位角分布不需太均匀、反演结果对速度模型和地壳结构横向不均匀性依赖较小的特点，日前国内的一些研究结果充分证明了CAP方法在震源机制解与地震矩心深度研究方面的有效性与可靠性（Zhao，Helmgerger，1994；Zhu，Helmgerger，1996；龙锋等，2019；吕坚等，2019）。CAP方法的主要原理为：任意一个双力偶震源的理论合成位移s（i）可表示为
洱源Ms5.5地震发生于滇西北裂陷带南段，位于印支板块和川滇菱形块体的交汇区域，该区域地质构造复杂，活动断裂较多，主要分布有NW向红河断裂、NE向程海断裂、北北西的维西一巍山断裂及走向基本与其一致的其它次级断裂。
可靠的震源机制解结果可反映震源区的受力状态，发震断层及其运行方式，进而可分析地震的成因，对于了解地震的孕育、发生机理及区域构造动力学环境等具有十分重要的意义（郑勇等，2019；龙锋等，2019；吕坚等，2019），单个地震的震源机制解可反映出在构造应力场作用下地震断层活动的性质，多个地震的震源机制解可反映出较人区域内的平均构造应力状况（龙海英等，2019；张建国等，2019；曹颖等，2019；虎雄林等，2019）。本文采用CAP方法，利用云南数字台网宽频带波形资料，反演计算了洱源Ms5.5地震序列的主震及其发生前后ML≥3.4地震的震源机制解，并对比了研究区域历史中强地震震源机制解结果，同时结合主震的余震分布情况及区域地质构造，探讨了该地震序列的发震断层面。
祁连山中东段基于CAP方法的中小地震震源机制研究
祁连山中东段基于CAP方法的中小地震震源机制研究
张辉;徐辉;姚军;代炜
【期刊名称】《山西地震》
【年(卷),期】2010(000)001
【摘要】利用甘肃省数字地震台网记录的区域波形资料,采用CAP方法计算了2001年至2009年发生在祁连山地震带中东段的37次中小地震震源机制解;采用格点搜索法对这37次地震震源深度进行了重新计算,获得了较理想的结果.结果表明,震源机制解具有分区特征,N38°以北地震以走滑型地震为主,N38°以南的地震以走滑型及走滑兼逆冲型为主.主压应力主要有两个优势方向NE和NEE向,仰角大多数小于30°,说明作用力以水平作用为主.通过与哈佛大学中强地震震源机制结果对比发现,该地区的中小地震震源机制的反演结果与强震震源机制结果有较好的一致性.
【总页数】5页(P25-29)
【作者】张辉;徐辉;姚军;代炜
【作者单位】甘肃省地震局,兰州,730000;甘肃省地震局,兰州,730000;甘肃省地震局,兰州,730000;甘肃省地震局,兰州,730000
【正文语种】中文
【中图分类】P35.63
【相关文献】
1.基于CAP方法的震源机制研究 [J], 赵凌云;邓津;陈俊华;申学林;戴苗;魏贵春
2.北天山中东段中小地震震源机制解及应力场反演 [J], 龙海英;高国英;聂晓红;李莹甄;沈军
3.利用CAP方法快速计算云南地区中小地震震源机制解 [J], 杨军;苏有锦;陈佳;叶泵;李孝宾;金明培;王宝善
4.祁连山中东段地区中小地震震源机制和区域应力场 [J], 荣代潞
5.祁连山中东段地区非弹性衰减系数、震源参数和场地响应研究 [J], 郭晓;张元生;莘海亮;张辉
因版权原因，仅展示原文概要，查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
去除归一化振幅并使用距离比例系数使得反演结 果更加稳定 。Ｃ Ａ Ｐ方法 【 ９ ＿ ２ Ｉ 】 目前 已得 到广泛的
应用 。
个有效途径。因此 ， 确 定 地 震 的 震 源 机 制
对 于地 震 本 身 的研 究 、 地 震 孕 震 机理 的 解 释 及 震后 应 力 的 分 布 等 具 有 十 分 重 要 的意 义 。 Ｐ波 初 动和 波形 反 演方 法 是 目前研 究 地 震 震 源 机制 比较 常 用 的 手 段 ， 波 形 反演 的 方 法 有 多种 ， 最典 型 的 是 体 波 反 演¨ 或 面 波 反 演 法ｌ ３ 】 ， 国 内 外 地 震 学 者 对 此 开 展 了 大量 的研 究工作 。 ’ 】 。２ ０世纪 ９ Ｏ年代 以前 ， 主要 采用模 拟记录资料 中的 Ｐ波初动及振幅 比信息研究地
合和反演确定地震矩张量的方法。陈运泰等 利用近震源宽频带资料进行了矩 张量 反演研究。
近年来 ， 区域数字地震 台网布设 了大量的宽频带
地震台站 ， 获得了高质量的地震数据 ， 利用近震地 震波形反演震源机制成为可能。Ｚ ｈ ａ ｏ 等 提 出
了通过 分 割波形 记 录为 Ｐ ｎ ｌ和 Ｓ ｎ ｌ部 分 ， 赋 予 不
姚军 ， 徐辉’
（ 甘肃省地震局 ， 甘 肃 兰州 ７ ３ ０ ０ ０ ０ ）
摘要： 选 用区域数 字地震 台网的 波形 资料 ， 利用 Ｃ Ａ Ｐ方 法反 演 了 ２ ０ １ １ 年 ３月 ３日久 治一 玛 曲 ４ ． １ 级 地震 震 源机 制解 。反 演结 果 显 示 ， 本 次地 震 的 最佳 双 力偶 解 为 节 面 Ｉ走 向 ３ ０ ６ 。 、 倾角８ ２ 。 、 滑动 角
一
５ 。 ， 节面 Ⅱ走向 ３ ７ 。 、 倾角 ８ ５ 。 、 滑动 角 一 １ ７ ２ 。 ， 地震类型 为走滑型 ， 震源深度为 ５ ． ５ ｋ ｍ， 属于汶川地震
关键词 ： 数字地震台网； Ｃ Ａ Ｐ方法； 震源机制
中图分类号 ： Ｐ ３ １ ５ ． ３ 文献标 识码 ： Ａ 文章编号 ： １ ０ ０ ５— ５ ８ ６ Ｘ｛ ２ ０ １ ３ ） ０ ２— ０ ０ ０ ８—０ ５
震 的震 源 机 制 。Ｈ ｅ ｌ ｍ ｂ ｅ ｒ ｇ ｅ ｒ等 利 用 地 震 Ｐ波 （ Ｐ 和Ｐ ｎ １ ） 模拟， 通 过 理 论 和 实 际 波 形 的对 比 ， 采 用 格 点搜 索 的方 法 ， 得 到 了地 震 的震 源 机 制 Ｊ 。随 着 矩张 量方 法 的 出现 ， 利 用 长 周 期 Ｐ波 资 料 ， 根 据 广义 射 线 理 论 】 ， 姚振兴 等ｌ １ 。 。 “ 通 过 波 形 拟
２ ０ １ １ 年 ３月 ３日在青海省果洛州久治县与 甘肃省甘南州玛曲县交界发生了 ９ ４ ． １ 级地震 ，
距玛 曲县城约 ９ ３ ｋ ｍ。甘肃省地震台网初步测定 该次 地震 的震 源 深 度 为 ５ ｋ ｍ。本 文 利 用 甘 肃 省
“ 十五” 宽频带数字化地震记录 ， 采用 Ｃ Ａ Ｐ 方法反 演本次地震的震源机制解 ， 通过空 间搜索确定本
第２ ５卷 第 ２期 ２ ０ １ ３年 ６月
高 原 地 震
Ｐ Ｌ Ａ Ｔ Ｅ Ａ Ｕ Ｅ Ａ Ｒ Ｔ Ｈ Ｑ Ｕ Ａ Ｋ Ｅ Ｒ Ｅ Ｓ Ｅ Ａ Ｒ Ｃ Ｈ
Hale Waihona Puke Ｖｏ １ ． ２ ５ Ｎ ｏ ． ２
Ｊ ｕ ｎ ． ２ ０ １ ３
用 Ｃ Ａ Ｐ方 法 研 究 久 治 一 玛 曲 Ｓ ４ ． １级 地 震 震 源 机 制
・ 通讯作者 ： 徐辉 （ １ ９ ７ Ｏ 一
） ， 女， 陕西汉中人 ， 工程师 ， 主要从事地震分析预测研究 。
第 ２期
姚军， 等： 用Ｃ Ａ Ｐ方 法研 究 久治一 玛 曲 ｓ ４ ． １ 级 地震 震 源机制
９
ｓ （ ｔ ）＝ 眠
．
Ａ （ 一 ， ， Ａ ） Ｇ （ ｔ ）
次地 震 的震源 深度 。
１ 方法与原理
Ｃ Ａ Ｐ （ Ｃ ｕ ｔ ａ ｎ ｄ Ｐ ａ ｓ ｔ ｅ ） 方 法是 一 种全 波 形 反演 的方法 ， 具 有 所需 台站 少 、 反 演结果 对速 度模 型和
地壳结构横向变化 的依赖性相对较小等优点 。该 方 法将 宽频带 地震 记 录分成 Ｐ ｎ ｌ 和Ｓ ｎ ｌ 两个 部分 进行反演并允许他们相对 浮动， 在适 当的时间变 化范围内 ， 搜索出合成地震 图和观测地震图全局 差异最小 的震源机制解 ， 提高了Ｐ ｎ ｌ 的权重 ， 且对 地震深度有 比较好 的约束。 基本 原理 ： 任 意 一个 双力 偶 震 源 的理 论 合 成 位移 ｓ （ ｔ ） 表示为 ：
收稿 日期 ： ２ ０ １ ３一ｏ ４－２ ７
基金项 目： 国家 自然基金（ ４ １ １ ７ ４ ０ ５ ９ ） 和中国地震局兰州地震研究所地震科技发展基金 （ ２ ０ １ ２ Ｍ ０ １ ） 共 同资助 。
作者简 介： 姚军（ １ ９ ６ ｏ 一 ） ， 女， 山西太原人 ， 工程师 ， 主要 从事地震分析预测研究。
（ １ ）
２ 数据分析 与反演结果
２ ． １ 数 据基 本情 况
式 中， ＝ｌ ， ２ ， ３时 ， 分别 对应 垂直 走 向滑 动 断 层、 垂 直倾 向滑 动 断层 和 ４ ５ 。 倾 向 滑 动 断 层 三 种 最基本 的 断层 类 型 。 为标 量地震矩 ， Ｇ ｉ 为 格
后 区域 应 力场调 整 引发 的一次地 震 事件 。
０ 引 言
地 震 震 源 机 制 从 一 定 角 度 描 述 了 震 源 的 性 质及 其 破 裂 过 程 ， 是 了 解 地 震 应 力 状 态 的
一
同的权重 ， 通过格点搜索的方法 ， 进行近震震源机 制的反演。Ｚ ｈ ｕ 等【 １ 进一步改进了该方法 ， 通过