工程地震震源机制解
利用震源机制解走向判定地震影响场长轴方向
第45卷第1期2022年1月地震研究JOURNAL OF SEISMOLOGICAL RESEARCH Vol.45,No.1 Jan.,2022第1期徐志双等:利用震源机制解走向判定地震影响场长轴方向此手段绘制了2014年新疆于田7.3级、云南鲁甸6.5级地震震动图。
前人研究成果和近年来应急经验表明,基于震源机制解走向判别地震影响场长轴方向在一些震例中是切实可行的,即修正后的地震影响场长轴方向与实际调查的极震区长轴方向相符性较好。
然而,这一方法是否适用于所有震例,以及评估偏差是多少,目前尚未有系统性的研究成果。
鉴于此,本文选取19702020年145个M≥5.0破坏性历史地震,研究其震源机制解走向和极震区长轴方向之间的偏差,考察以震源机制解走向作为判别地震影响场长轴方向是否满足震后快速评估的需求。
地震发生后,由于震源机制解产出需要一定的等待时间(约震后30 min),是否可以借鉴活动断层数据库使用方法,在震前准备好历史地震震源机制解走向数据库,震时直接读取震中附近历史强震震源机制解走向作为地震影响场长轴方向?为了检验该方法的可行性,本文对20102020年等震线长轴方向明显的50个强震及距震中最近的历史强震震源机制解走向与极震区长轴方向差值进行了统计研究。
1研究资料选取和震例分析震源机制解一般包含走向(φS)、倾角(δ)和滑动角(λ)3个参数,走向和倾角是断层的几何参数,描述发震断层的产状;滑动角是断层的运动参数,表示断层面上的滑动方向。
根据滑动角的不同,可将断层分为走滑-近走滑型(λ接近0°或±180°)、正断-近正断型(λ接近-90°)、逆断-近逆断型(λ接近90°)和其他类型(λ为其他值)。
本文以19702020年在中国大陆发生的破坏性地震为研究对象,整理了震源机制解和等震线资料齐全、等震线具有明显方向性的145个地震,历史地震空间分布如图1所示,除1980年2月14日新疆叶城5.8级地震震源深度达90km外,其余地震均为浅源地震。
基于震源机制解的分区地震烈度衰减关系研究
基于震源机制解的分区地震烈度衰减关系研究郑韵;姜立新;杨天青;刘杰【摘要】根据中国大陆构造应力场分区,本文拟合不同震源机制解的烈度衰减关系,收集了我国西部198次5.0级以上震例,共419条的等震线记录;东部47次5.0级以上震例,共84条等震线记录,采用椭圆模型进行衰减关系拟合,得到了地震烈度随震级和长短轴长度的衰减关系.文章初步分析了几种震源机制地震烈度衰减规律的差异,以西部地区为例,对比了西部走滑型、逆断型烈度衰减关系和汪素云西部衰减关系与实际值之间的差距.本文研究结果表明,区分发震方式的烈度衰减关系长短轴可以作为应急地震影响场的修正因子,在震后几小时得到震源机制解时,对烈度圈长短轴进行初步修正.【期刊名称】《震灾防御技术》【年(卷),期】2016(011)002【总页数】11页(P349-359)【关键词】烈度衰减关系;震源机制解;构造应力场分区【作者】郑韵;姜立新;杨天青;刘杰【作者单位】福建省地震局,福州350003;中国地震台网中心,北京100045;中国地震台网中心,北京100045;中国地震台网中心,北京100045【正文语种】中文郑韵，姜立新，杨天青，刘杰，2016.基于震源机制解的分区地震烈度衰减关系研究.震灾防御技术，11（2）：349—359. doi：10.11899/zzfy20160218地震烈度衰减关系在抗震设防和震后烈度快速评估等工作中，具有重要的实际意义。
我国幅员辽阔，地质结构复杂，很多研究者已对不同地区、不同地质构造环境的地震烈度衰减关系进行了研究，得出了适用于各个地区的统计关系。
汪素云等（2000）将中国分为东、西部两个区，建立了分区的烈度衰减关系；有些研究人员考虑了地质构造环境和震害分布特点，将川滇地区分为西南地区和盆地地区，分别确定了烈度衰减规律（雷建成等，2007；孙继浩等，2011）；谭明等（2011）分不同时段和区域，得到了符合新疆地区的衰减关系。
震源机制解的基本概念及图形解释
震源机制解图形解释
• 倾向与倾角
节面凸向哪方，即倾向哪方（下半球投影有效） 节面越弯越靠近边缘，倾角越缓（最小0°） 节面越直越靠近中央，倾角越陡（最大90°）
震源机制解图形解释
• 力轴与错动类型
约定：沙滩求白色部分为拉张区，初动向下，为P轴所在区 域；深色部分为压缩区，初动向上，为T轴所在区域 力轴位于各象限的中心部位，并在对称象限成对出现 两个节面相交的位置为B轴 当P轴位于沙滩球中心，错动类型为正断层；当T轴位于中心， 为逆断层；当B轴位于中心，为走滑类型
震源机制解的基本概念 及图形解释
Source Parameters
• 震源 (纬度, 经度, 深度) • 发震时刻 (初始破裂时间) • 震级(地震大小度量) • 断层类型 (震源机制解) • 断层尺度 (长度,宽度 ,位错) • 地震矩 (地震大小度量) • 应力降(剪切应力改变量) • 震源过程 (破裂过程)
Right-lateral strike slip faulting 右旋走滑
走滑断层——九寨沟M7地震
Reverse faulting 逆断层
上盘
下盘
上盘在水平向构造应力的作用下上升，下盘相对下降。
逆断层——汶川地震
Normal faulting 正断层
下盘
上盘
上盘在重力作用下下降，下盘相对上升。
Surface rupture in 1999 Taiwan, Chi-Chi earthquake (Taken by Prof. Abe, the University of Tokyo)
Stritrike slip faulting 左旋走滑
Why source mechanism?
The source mechanism also guide you in the state of the tectonic stress field and location of the week zone (fault zone).
震源机制解
（°）
40.90
40.44 40.40 39.69 40.16 40.83 40.64 39.59 39.96 40.37 39.66 40.10 40.08 40.94 41.19 47.85 45.54 60.94
地心
F — 震源；O — 震中；S1, S2, … , Sn — 观测点
近震射线、远震射线与投影球面的关系
S2 O S3 S4 S1
至…
远震射线
近震射线
F
至 Sn-1
至 S5
至 Sn
近震射线大多穿过上半球面，远震射线多数穿过下半球面
上半球射线的投影方法
N
S1 上半球射线 i1
W S1′ S N S2′
辐射图案与各个方位上地震波波形的对应关系
Hale Waihona Puke 0°y0°30°
60°
90° 90° 120°
270°
x
180°
150° 180°
倾斜节面（断层面）与水平面的交切
节面 a、b 将两者的法 平面分割为四个象限。 P波初动在法平面内呈 对称的四象限分布。
节面 a、b 倾斜正交
a
c b
c
b
节面 a、b 与平面 c 相 交，将平面 c 划分为四 部分，P波初动在水平 面内也呈对称分布。
－Pn
＋Pn －Pn Pn ＋Pg Pg ＋Pg ＋Pn ＋Pn ＋Pg －Pg －Pg －Pg －P －P －P －P ＋P
地震震源机制与构造应力场演化的时空关联研究
地震震源机制与构造应力场演化的时空关联研究地震是地球内部能量释放的一种现象，研究地震的震源机制和构造应力场演化对于理解地球的内部运动和地震灾害预测具有重要意义。
本文将探讨地震震源机制与构造应力场演化的时空关联，分析地震与构造应力场相互作用的关系。
一、地震震源机制地震震源机制是指地震发生时断层的滑动类型和方向。
地震震源机制的研究可以揭示地壳微观应力和变形的特征，从而帮助科学家理解地震的成因和发展趋势。
通过使用地震仪器记录地震波，地震学家可以推测地震发生时断层的滑动情况。
根据震源机制，地震可以分为正断层震和逆断层震，也可以确定断层的滑动面和运动方向。
二、构造应力场演化构造应力场是指地壳中的应力分布情况。
地壳的构造应力来源于板块运动、重力和岩石变形等因素。
构造应力场的研究可以帮助科学家预测地震活动的地点和强度。
构造应力场的演化与地震的产生和发展密切相关。
随着地壳的变形，构造应力场会发生变化，进而影响地震的发生。
三、地震震源机制与构造应力场演化的关系地震震源机制与构造应力场演化之间存在着紧密的关联。
首先，地震的发生会导致应力场的重新分布。
当地震发生时，断层上原先积累的应力得到释放，周围的应力会重新分布。
这种应力分布的再调整可能会引发附近的断层发生滑动，从而引发更多的地震。
其次，构造应力场的演化也会影响地震的发生。
当地壳受到应力的作用时，对断层的应力状态会有所改变。
如果构造应力场的演化给予断层足够的应力，就会促使地震的发生。
四、地震震源机制与构造应力场演化的时空关联地震震源机制与构造应力场演化具有明显的时空关联。
首先，在时间上，地震的发生和构造应力场的演化是连续进行的。
地壳中的构造应力在演化过程中会不断积累，最终导致地震的发生。
其次，在空间上，地震一般发生在构造应力场较大的区域。
这是因为构造应力场较大意味着断层上的应力积累相对较高，容易引发断层滑动和地震。
综上所述，地震震源机制与构造应力场演化之间存在着密切的时空关联。
张家口地区近年中小地震震源机制解及应力场反演
地 震EARTHQUAKE 第39卷第1期2019年1月Vol. 39, No. 1Jan., 2019张家口地区近年中小地震震源机制解及应力场反演**收稿日期：2017-09-28；修改回日期：201740-09基金项目：河北省地縄局地纔科技星火计划项目"晋冀蒙交界地区典型盆地孕震动力学模型构建(DZ20160621076)"作者简介：王曰风(1974-).男.甘肃民勤人.高级工程师.主要从事地震学和地霍分析预测综合研究等工作王曰风，张秀萍，李峰，马利军，宋晓煜，刘燕翔，张珊珊(河北省地震局张家口中心台，河北 张家口 075000)摘要：利用垂直向Pg 和Sg 波最大振幅比方法.计算得到张家口地区2007 -2015年发生的中小地震震源机制解；通过系统聚类分析.研究震源机制分布特征与断裂活动的关系.并采用网格搜索法对张家口地区的现今地壳应力场进行反演.得到了较为精确的张家口地区的构造应力场。
结果表明：主压应力轴的走向范围：211°〜223°.即N E —SW 向；主张应力轴的走向范围：323°〜332°,即NNW -SSE 向，P 轴、T 轴都近于水平；两组节面的走向分别为N18°E 向和近EW 向(N83°E),破裂以水平走滑为主。
关键词：震源机制；系统聚类分析；网格搜索法；构造应力场；张家口地区中图分类号:P315. 7 文献标识码：A 文章编号:1000-3274(2019)01-0011-09引言地震是由于地球内部的岩体受到构造应力的作用，导致岩体突然断裂错动的结果。
目 前研究地壳应力状态较为可行并普遍采用的途径是通过测定大量的中、小地震的震源机制 解来分析研究发震断层的力学性质和动力学特征，大量的中、小地震的震源机制解的优势 结果可以反映区域构造应力场的特征。
20世纪80年代以来，梁尚鸿等⑴提出利用Pg 、Sg 最大振幅比资料反演中小地震的震源机制参数的方法，该方法利用区域台网地震仪垂直向 记录的Pg 、Sg 振幅比资料计算震源机制解.因其具有精度高、需要的记录资料量少等优点 而受到广泛应用.胡新亮等⑷还对本方法计算的震源机制解的可靠性进行全面深入分析. 经过实例对比分析•证实了利用区域地震台网Pg 、Sg 振幅比资料求得震源机制解的可靠 性。
arcgis画震源机制解
arcgis画震源机制解
在实际地震工作中，有时需要绘制震源机制球（沙滩
球,BEACHBALL）。
为此，编制了一个ARCGISPYTHON脚本工具，可用于绘制SHAPEFILE格式的震源机制球。
图为利用该工具绘制的震源机制球。
由于球体为SHAPEFILE格式，不会产生不清晰的问题，并能叠加到任何地理底图上。
该工具需要准备的数据为包含震级，走向，倾角，滑动角（皆为数值属性）的ArcMap地震目录（SHAPEFILE点要素），工具可一次性批量绘出所有震源球。
1.修正了一个浮点截断误差造成ACOS函数超域的Bug；
2.增加了一个光滑度选项；
3.增加震源球与震中不叠压选项。
修正了一个与不同小数点格式相关的Bug.
1.增加了一个新的震源球叠压模式选项
2.提高了震源球分离速度。
v2.1版主要变化
1.提升震源球分离速度。
2.震源球要素坐标范围有效性检查
3.震源机制解参数有效性检查
4.空震源点(点数与表记录数不匹配)检查
5.其它细节变更。
乌鲁木齐地区中小地震震源机制解及构造应力场
关键词 ： 中小 地 震 ； 源机 制解 ； 造 应 力 场 ； 源 断 错 性 质 ；乌鲁 木 齐 地 区 、 震 构 震
中图 分 类 号 ： ３５ ７ Ｐ １ ． 文 献标 识码 ： Ａ 文 章 编 号 ：０ ０３ ７ （０ ７ ０～０ ９０ １ ０—２ ４ ２０ ）３０ ８—８
引 言
＊
收 稿 日期 ： ０ ７０ ～５ 修 改 回 日期 ：２０ —２１ ２ ０ —１１ ； ０７０ —３ 基 金 项 目：乌 鲁 木 齐 城 市 活 断层 探 测 项 目子 专 题 作 者 简 介 ：龙 海 英 （ ９ ３） １ ７一 ，女 ，四 川仪 陇人 ，副 研究 员 ，１９ 年 获 得 硕 士 学 位 ， 要 从事 地震 综 合 分 析 预 报 等研 ９９ 主
（ 新疆 维 吾 尔 自治 区地 震 局 ， 疆 乌 鲁 木齐 新 ８ ０ ｔ） ３ ０ １
摘 要 ： 要 对 乌 鲁 木 齐地 区 中小 地 震 震 源 机 制 解 及 构 造 应 力 场 进 行 分 析 。结 果 表 明 ，乌 鲁 木 齐 主 地 区 中等 地 震 震 源 断 错 性 质 主 要 以 倾 滑 逆 断 层 为 主 ， 震 破 裂 面 与 其 附 近 构 造 走 向基 本 一 致 ， 地 主 压 应力 Ｐ 轴 方 向 为 Ｎ 向或 近 ＮＳ向 ； 鲁 木 齐 地 区 小 地 震 的 震 源 断错 性 质 表现 出一 定 的 区 Ｅ 乌 域 性 ，破裂 面 以 近 Ｅ 向为 主 ，主压 应 力 Ｐ 轴 方 向 近 ＮＳ向 ， 有 部 分 Ｎ 向 或 者 Ｎ 向 。显 Ｗ 兼 Ｅ Ｗ 示 出 小 地 震 的 发 生 既 受 局 部 地质 构造 的影 响 ， 受 区域 构 造 应 力 场 的 影 响 。 也
利用震源机制解走向判定地震影响场长轴方向
利用震源机制解走向判定地震影响场长轴方向作者：徐志双任静谭专条高小跃陈雅慧杨志高李志强来源：《地震研究》2022年第01期摘要：為了评估以震源机制解走向作为影响场长轴方向的适用性，以1970—2020年145个破坏性地震为研究对象，研究各地震的震源机制解走向和极震区长轴方向之间的偏差，发现两者平均差值为17.0°;走滑型、逆断型和正断型地震分别占研究地震总数的56.6%、19.3%和6.2%，其震源机制解走向与极震区长轴方向的平均差值分别为16.4°、16.2°和20.6°;有82.8%的地震的震源机制解走向与极震区长轴走向差值小于30°，说明该方法是可行的。
通过统计研究2010—2020年等震线长轴方向明显的50个强震震中附近区域历史强震震源机制解走向与极震区长轴方向差值，结果表明：有13个地震震中30 km范围内有历史强震发生，最近历史强震震源机制解走向与极震区长轴方向差值小于30°的有8个（62%）。
如果在距离地震震中30 km以内发生过历史强震，则可以综合考虑距震中最近的历史强震震源机制解走向和活动断裂走向来判定影响场长轴方向。
关键词：地震影响场;长轴方向;震源机制解;等震线图;极震区中图分类号：P315.94 文献标识码：A 文章编号：1000-0666（2022）01-0088-12doi：10.20015/ki.ISSN1000-0666.2022.00100 引言震后的灾害快速评估工作主要是对地震影响范围内的人口、经济、建筑物等进行综合考量，因此绘制地震影响场是震害快速评估的第一步。
地震发生后，地震灾害快速评估系统一般依据震中所在区域的烈度衰减关系（汪素云等，2000;张苏平等，2015）和活动断裂分布情况绘制地震影响场，以作为震害评估的基础和应急救援的依据（刘金龙，林均岐，2012）。
基于地震烈度衰减关系生成的地震影响场一般为一系列同心椭圆，椭圆长轴方向取距震中最近的活动断裂的走向。
云南楚雄M4.7地震震源机制解分析
第１期（总第１７８期）２０２１年３月四川地震 ＥＡＲＴＨＱＵＡＫＥＲＥＳＥＡＲＣＨＩＮＳＩＣＨＵＡＮＮｏ．１Ｍａｒ．２０２１收稿日期：２０２０－０８－１０；修回日期：２０２０－１１－１２作者简介：杨云存（１９８８－），女，工程师，主要从事地震监测与分析工作．Ｅ－ｍａｉｌ：５２６３８１６５９＠ｑｑ．ｃｏｍ．云南楚雄Ｍ４．７地震震源机制解分析杨云存，王华柳，段 毅，孙自刚（云南省地震局，云南昆明６５７１００）摘 要：选用云南省地震台网宽频带地震波形数据，采用ＣＡＰ震源机制解反演软件计算了２０１９年６月２４日云南楚雄Ｍ４．７地震的震源机制解并初步分析了该地震的发震构造特征。
结果显示：地震的矩震级为Ｍｗ４．７６，最优震源矩心深度为１３ ８ｋｍ，表明此次地震发生在上地壳浅部。
反演结果显示节面Ⅰ走向３２９°、倾角７３°、滑动角－１７１°；节面Ⅱ走向２３６ ３°、倾角８１ ４°、滑动角－１７ ２°。
震源机制解揭示此次地震的发震断层呈右旋走滑，与南华—楚雄断裂带走向一致，推测南华—楚雄断裂为此次地震的发震断裂。
关键词：ＣＡＰ方法；震源机制解；发震构造中图分类号：Ｐ３１５．３３ 文献标识码：Ｂ 文章编号：１００１－８１１５（２０２１）０１－００４４－０４ＤＯＩ：１０．１３７１６／ｊ．ｃｎｋｉ．１００１－８１１５．２０２１．０１．００９据中国地震台网中心测定，２０１９年６月２４日２１时２４分２１秒在云南省楚雄彝族自治州楚雄市（１０１ ６７°Ｅ，２４ ９８°Ｎ）发生Ｍ４ ７地震，震源深度１０ｋｍ。
此次地震发生在南华—楚雄断裂附近。
历史资料显示，１５１１年至今该断裂发生过９次５级以上地震，最大地震为１６８０年９月９日发生的Ｍ６ ８地震，最近一次中强地震发生在２００１年，震级为Ｍｓ５ ３。
震源机制解被认为是确定地震发震构造的关键依据，能够描述震源的性质及其破裂过程，并且能为分析孕震机理以及区域构造动力学环境提供参考（祁玉萍等，２０１３）。
基于震源机制解分析邢台地震的发震构造
第４５卷㊀第６期２０２３年１１月地㊀震㊀工㊀程㊀学㊀报C H I N A E A R T H Q U A K EE N G I N E E R I N GJ O U R N A LV o l ．４５㊀N o ．６N o v e m b e r ,２０２３㊀㊀收稿日期:２０２２Ｇ０８Ｇ０６㊀㊀基金项目:中央高校科研业务费专项(Z Y ２０２１５１１７);河北地震科技星火计划项目(D Z ２０２１１２１４００００６,D Z ２０２１１１０５００００２,D Z ２０２３１２０５００００９)㊀㊀第一作者简介:章㊀阳(１９８９－),男,工程师,主要从事数字地震学与地壳应力场方面的研究.E Ｇm a i l :z h a n g y a n g _o u c ＠１６３．c o m .㊀㊀通信作者:万永革(１９６７－),男,研究员,主要从事构造应力场㊁地震应力触发等方面研究工作.E Ｇm a i l :w a n y g ２１７２１７＠v i p．s i n a ．c o m .章阳,万永革,赵志远,等．基于震源机制解分析邢台地震的发震构造[J ]．地震工程学报,２０２３,４５(６):１４３９Ｇ１４４８．D O I :１０．２００００/j．１０００Ｇ０８４４．２０２２０８０６００２Z HA N G Y a n g ,WA N Y o n g g e ,Z HA OZ h i y u a ,e t a l ．D e t e r m i n a t i o no f t h eX i n g t a i e a r t h q u a k e s e i s m o ge n i c s t r u c t u r e b a s e do n t h ef o c a lm e c h a n i s ms o l u t i o n [J ]．C h i n a E a r t h q u a k e E ng i n e e r i n g J o u r n a l ,２０２３,４５(６):１４３９Ｇ１４４８．D O I :１０．２００００/j ．１０００Ｇ０８４４．２０２２０８０６００２基于震源机制解分析邢台地震的发震构造章㊀阳１,２,万永革３,４,赵志远１,２,陈㊀婷１,２(１．河北红山巨厚沉积与地震灾害国家野外科学观测研究站,河北邢台０５４０００;２．河北省地震局,河北石家庄０５００２１;３．防灾科技学院,河北三河０６５２０１;４．河北省地震动力学重点实验室,河北三河０６５２０１)摘要:活动断层几何形状的确定为评估一个地区的地震危险性提供重要的理论依据,邢台老震区内构造背景复杂,前人对邢台地震的发震构造展开大量研究,取得丰富成果.通过震源机制节面聚类得到邢台地震发震断层的走向㊁倾角及其标准差,综合前人研究结果求取发震断层面的平均解;结合区域地壳应力张量,估计该断层的滑动方向及其误差,分析该区域的地震活动危险性.结果表明:邢台地震发震断层的走向为３２．４５ʎ,倾角为７９．４４ʎ,滑动角－１５３．９６ʎ,标准差为４．５５ʎ.该断层为走向N N E ㊁倾向NWW 的高倾角走滑型断层,而非正断性质的新河断裂.区域地壳应力场在断层面上产生的相对剪应力为０．８３,相对正应力为０．５９.该断层不是最大剪应力的断层面形状,但剪应力强度仍较大,表明该断层在地震能量积累和释放过程中仍起到必不可少的作用.关键词:邢台地震;活动断层;震源机制节面聚类;区域应力场;地震危险性中图分类号:P ３１５㊀㊀㊀㊀㊀㊀文献标志码:A㊀㊀㊀文章编号:１０００Ｇ０８４４(２０２３)０６－１４３９－１０D O I :１０．２００００/j．１０００Ｇ０８４４．２０２２０８０６００２D e t e r m i n a t i o no f t h eX i n g t a i e a r t h q u a k e s e i s m o ge n i c s t r u c t u r e b a s e d o n t h ef o c a lm e c h a n i s ms o l u t i o nZ H A N G Y a n g １,２,WA N Y o n g g e ３,４,Z H A OZ h i y u a n １,２,C H E N T i n g１,２(１．H e b e iH o n g s h a nN a t i o n a lO b s e r v a t o r y o nT h i c kS e d i m e n t s a n dS e i s m i cH a z a r d s ,X i n gt a i ０５４０００,H e b e i ,C h i n a ;２．H e b e iE a r t h q u a k eA g e n c y ,S h i j i a z h u a n g ０５００２１,H e b e i ,C h i n a ;３．I n s t i t u t e o f Di s a s t e rP r e v e n t i o n ,S a n h e ０６５２０１,H e b e i ,C h i n a ;４．H e b e iK e y L a b o r a t o r y o f E a r t h q u a k eD yn a m i c s ,S a n h e ０６５２０１,H e b e i ,C h i n a )A b s t r a c t :T h e d e t e r m i n a t i o no f t h e g e o m e t r y o f a na c t i v e f a u l t p r o v i d e s a n i m po r t a n t t h e o r e t i c a l b a s i s f o r a s s e s s i n g t h e s e i s m i c r i s ko f a r e g i o n ．T h e s t r u c t u r a l b a c k g r o u n do f t h eX i n gt a i e a r t h Ｇq u a k e z o n e i sc o m p l i c a t e d ．M a n yp r e v i o u ss t u d i e so nt h es e i s m o g e n i cs t r u c t u r eo f t h eX i n gt a i e a r t h q u a k eh a v eb e e n c o n d u c t e d ,a n d p l e n t y o f r e s u l t s h a v e b e e n r e p o r t e d ．I n t h i s p a p e r ,w e e v a l Ｇu a t e t h e s t r i k e a n dd i p a n g l e s o f t h e s e i s m o g e n i c f a u l t o f t h eX i n g t a i e a r t h qu a k e a sw e l l a s t h e i rs t a n d a r dd e v i a t i o nb y c l u s t e r i n g n o d a l p l a n e s o f f o c a lm e c h a n i s m s．T h e n,t h e a v e r a g e s o l u t i o no f t h e s e i s m o g e n i c f a u l t p l a n e i so b t a i n e db y s y n t h e s i z i n g t h e p r e v i o u sr e s u l t s．F i n a l l y,u s i n g t h e c r u s t a l s t r e s s t e n s o r i nt h e r e g i o nu n d e r i n v e s t i g a t i o n,w ee s t i m a t e t h es l i d i n g d i r e c t i o na n d i t s s t a n d a r dd e v i a t i o na n da n a l y z e t h e s e i s m i c a c t i v i t y r i s k．T h e c a l c u l a t e dv a l u e s o f t h e s t r i k e,d i p, a n d r a k e a n g l e s o f t h e s e i s m o g e n i c f a u l t a r e３２．４５ʎ,７９．４４ʎ,a n d－１５３．９６ʎ,r e s p e c t i v e l y,w i t ha s t a n d a r dd e v i a t i o no f４．５５ʎ．T h e s ev a l u e s i n d i c a t e t h a t t h e f a u l tu n d e r i n v e s t i g a t i o n i s ah i g hＧd i p a n g l e s t r i k eＧs l i p f a u l t e x t e n d i n g f r o m N N Et oNWWr a t h e r t h a nan o r m a lX i n h e f a u l t．T h e e s t iＧm a t e d r e l a t i v es h e a ra n dn o r m a ls t r e s s e sa p p l i e dt ot h ef a u l ta r e０．８３a n d０．５９,r e s p e c t i v e l y．T h e s e v a l u e s a r eb a s e do n t h e r e g i o n a l s t r e s s f i e l d．T h e s h e a r s t r e s s a p p l i e d t o t h e f a u l t i sh i g h, i n d i c a t i n g t h a t t h e f a u l t i s e s s e n t i a l i n t h e a c c u m u l a t i o na n d r e l e a s e o f s e i s m i c e n e r g y．K e y w o r d s:X i n g t a i e a r t h q u a k e;a c t i v e f a u l t;c l u s t e r i n g n o d a l p l a n e s o f f o c a lm e c h a n i s m s;r e g i o n a l s t r e s s f i e l d;s e i s m i c r i s k０㊀引言邢台老震区构造背景复杂,发育N N E向㊁N E 向㊁NWW向和近E W向等多组断裂.１９６６年邢台大地震后,众多学者对其发震构造展开大量研究.最初推测邢台地震主震震源是N E向的右旋断层,与束鹿断陷内的断裂相对应,但实际观察中未发现这一大型走滑断层.陈运泰等[１]根据地形变观测资料认为邢台地震的震源过程是一个复合断层模式,而非简单的单一断层模式;徐杰等[２]根据石油地质勘探资料和人工地震剖面研究认为邢台７．２级地震的发震断裂不是单一形式的新河断裂或其下方的高角度断裂,而是二者的组合,高角度断裂是发震断裂的主要部分;赵成斌等[３]根据深地震反射资料认为邢台地震的可能发震断裂是震源下方的超壳深高倾角走滑断裂;徐锡伟等[４]分析认为邢台地震群是在最新构造应力场作用下,受北西向断层或横向阻碍体作用形成 新生断层 的产物;邵志刚等[５]依据同震观测位移及数值模拟分析认为,邢台地震群发震断层浅层为倾向S E的N E向牛家桥断层和永福庄断层,而深部主断层为倾向NW的N E向断层;王椿镛等[６Ｇ７]研究认为邢台地震发震构造为中下地壳的深部断裂.综合前人研究成果发现邢台地震发震构造是深部断裂,而不是新河断裂.为了定量研究该活动断层的性质,需求取该断层的几何参数.前人提出多种精确测定活动断层几何参数的方法,并取得丰富研究成果.对于出露地表的断层,可以通过实地踏勘测量出露断层面的几何形状[８].对于未出露地表的隐伏断层,确定断层几何形状主要有挖探槽法[９Ｇ１０]和地震反射探测法[１１Ｇ１３].近年来地震学发展了一种利用余震震源位置的空间分布精确勾画断层面形状和位置的方法[１４],被学者广泛应用于地震发震断层面参数的确定[１５Ｇ１９],该方法虽然优点明显,但是要求余震的丛集展布,如果余震未出现丛集,则该方法不适用,而且这种方法对震源定位精度要求较高.万永革[２０]提出一种独立于地质㊁大地测量数据和小震定位结果的确定活动断层几何形状的方法,对发生在断层上的中小地震震源机制节面进行聚类分析,获取断层面走向和倾角,并给出其置信区间.同时将该方法应用到漾濞地震序列中,推测该地震序列的发震断层,与其他学者研究结果基本保持一致,进而验证该方法的有效性.本文通过收集１９６６年邢台地震至今的震源机制解资料,对震源机制的节面进行聚类分析,推测丛集地震序列发震断层的几何形状.与前人研究结果比较,进而求取发震断层面的平均解.结合区域地壳应力张量,估计该断层的滑动方向及其误差.最后通过计算区域应力场在断层面上的相对剪应力和相对正应力分析该地区的地震活动危险性.１㊀震源机制节面聚类１．１㊀数据依据前人研究成果[２１Ｇ２５],收集整理出邢台老震区１２８个地震的震源机制解(表１),地震序列基本沿着N N E向新河断裂展布.邢台及周边区域地质构造背景及震源机制解具体分布见图１.将收集的１２８个地震震源机制解按照震源机制水平应变花面应变[２６]分类,结果显示,正断型地震(N)有１６个,占总体数量的１２．５％;正走滑型地震(N S)有１０个,占总体数量的７．８１％;走滑型地震(S S)有９０个,占总体数量的７０．３１％;逆走滑型地震(R S)有８个,占总体数量的６．２５％;逆断型地震(R)有０４４１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀地㊀震㊀工㊀程㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０２３年表１㊀邢台地震序列震源机制解T a b l e１㊀F o c a lm e c h a n i s ms o l u t i o n s o fX i n g t a i e a r t h q u a k e s e q u e n c e发震时间纬度/(ʎ)经度/(ʎ)震级走向/(ʎ)倾角/(ʎ)滑动角/(ʎ)１９６６Ｇ０３Ｇ０６３７．４６７１１５．０３３M S５．２１９８７９１７５１９６６Ｇ０３Ｇ０８３７．３９６１１４．９５７M S６．８１９８１－１７７１９６６Ｇ０３Ｇ２０３７．２６７１１４．９６７M S５．６１９３８０－１７４１９６６Ｇ０３Ｇ２２３７．５００１１５．０８３M S６．７１８６８３－１７８１９６６Ｇ０３Ｇ２２３７．５６６１１５．１７８M S７．２１９９８９－１７９１９６６Ｇ０３Ｇ２２３７．５１７１１５．０３３M S５．６１８６８０１６５１９６６Ｇ０３Ｇ２２３７．５６７１１４．９８３M S５．６１８３７５－１６９１９６６Ｇ０３Ｇ２４３７．４１７１１４．９００M S５．２１６７７－１５８１９６６Ｇ０３Ｇ２６３７．７３５１１５．１７２M S６．２３０７４－１５７１９６６Ｇ０３Ｇ２６３７．７００１１５．３００M S５．２２０５７０－１７１９６６Ｇ０３Ｇ２７３７．６００１１５．３００M S５．８１９３８２－１７８１９６６Ｇ０３Ｇ２９３７．４８６１１４．９９４M S６．０４８７２１４１１９６６Ｇ０４Ｇ２０３７．２５０１１４．７８３M S５．３３０８６０２３１９６６Ｇ０７Ｇ１９３７．２８３１１４．８１７M S５．１３８６１６９１９６７Ｇ１２Ｇ０３３７．７１５１１５．１５０M S５．７１５８７５１５７１９６８Ｇ０１Ｇ０７３７．６００１１５．１００M S４．７２２８７０－６１９６８Ｇ０１Ｇ１６３７．７７９１１５．５１７M S５．３１９８７２－６０１９６８Ｇ０５Ｇ１５３７．２００１１４．８００M S４．７３１７７－１６３１９６８Ｇ０５Ｇ１５３７．３００１１４．８００M S４．７２０６６２－１２４１９７２Ｇ０９Ｇ０７３７．３００１１４．８１７M S４．６１９９７６１４６１９７４Ｇ０６Ｇ０６３７．６００１１５．１００M S４．９２２０８３１７０１９８１Ｇ１１Ｇ０９３７．３７９１１５．０３６M S５．８３１７２１６２１９８２Ｇ０１Ｇ２６３７．４００１１４．８６７M S４．４２１０８０－１５３１９８５Ｇ１１Ｇ３０３７．２３３１１４．８１７M S５．３２３６８９１６１１９８６Ｇ１２Ｇ１６３７．４００１１５．５００M S３．５８４３９０２００２Ｇ０４Ｇ１６３７．２５０１１４．７８０M L２．８５４６０１７５２００２Ｇ０４Ｇ２２３７．２７０１１４．４２０M L４．７１５６２１－８２２００２Ｇ０４Ｇ２２３７．４３０１１４．６３０M L２．６１４８４０－２６２００２Ｇ０５Ｇ０５３７．５５０１１５．０００M L３．３９７８９－１６２００２Ｇ０８Ｇ１０３７．５２０１１５．０７０M L３．５２１７８８２４２００２Ｇ０９Ｇ０２３７．０７０１１４．６３０M L２．７２６２５１０２００２Ｇ１１Ｇ１０３７．３２０１１４．８７０M L３．７１１９３０－１２１２００２Ｇ１２Ｇ０２３７．４７０１１４．９００M L２．８１０４７５１６２００２Ｇ１２Ｇ２５３７．６００１１５．１００M L４．２１５４８－１０７２００２Ｇ１２Ｇ２７３７．６００１１５．０７０M L３．１２８３７６１６２００２Ｇ１２Ｇ２７３７．６２０１１５．０７０M L２．９６２３１－２５２００２Ｇ１２Ｇ３１３７．５８０１１５．１３０M L２．９１２０８５３９２００３Ｇ０１Ｇ１９３７．５８０１１５．１３０M L３．９３３８１１６０２００３Ｇ０１Ｇ１９３７．６００１１５．０７０M L３．０１０２７３１１２００３Ｇ０２Ｇ２０３７．２３０１１４．７５０M L３．０３１４４８－３８２００３Ｇ０８Ｇ２２３７．３８０１１４．９００M L２．７４１８１１６３２００３Ｇ０９Ｇ０６３７．３３０１１４．９２０M L３．６２９６７３－１４１２００３Ｇ１０Ｇ０２３７．４７０１１４．９５０M L３．７１１１６４－７９２００３Ｇ１０Ｇ０３３７．４７０１１４．９３０M L２．７２４６３１５９２００３Ｇ１１Ｇ２７３７．０５０１１４．３００M L２．５１８０３６９２２００３Ｇ１２Ｇ２８３７．４７０１１５．０００M L２．７６６３５２９２００４Ｇ０２Ｇ１０３６．９８０１１４．４００M L２．６２９２１９－１４２２００４Ｇ０３Ｇ１４３７．４００１１４．３２０M L３．６１９１６５－１５５２００４Ｇ０４Ｇ２１３７．５２０１１４．２８０M L３．２２６２３１１６５２００４Ｇ０５Ｇ０２３７．１３０１１４．４００M L２．５６１２０－３６２００４Ｇ０５Ｇ２４３７．４００１１４．５７０M L３．５２３６３４－８０２００４Ｇ０６Ｇ１７３７．５８０１１５．３５０M L２．８２０２７７－１３５２００４Ｇ０７Ｇ２３３７．４７０１１５．０００M L３．６２５０４５５２００４Ｇ０９Ｇ０３３７．４７０１１４．９２０M L２．８２９３５５７８１４４１第４５卷第６期㊀㊀㊀㊀㊀㊀㊀㊀㊀章㊀阳,等:基于震源机制解分析邢台地震的发震构造㊀㊀㊀㊀㊀㊀㊀㊀㊀２４４１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀地㊀震㊀工㊀程㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０２３年续表１发震时间纬度/(ʎ)经度/(ʎ)震级走向/(ʎ)倾角/(ʎ)滑动角/(ʎ)２００４Ｇ１２Ｇ１９３７．６２０１１５．１３０M L３．９３２６３１６８２００５Ｇ０１Ｇ０７３７．５８０１１５．１００M L２．６３４３１７－７９２００５Ｇ０２Ｇ２３３７．５５０１１５．０２０M L３．１１１５７６２４２００５Ｇ０３Ｇ０４３７．５５０１１５．０２０M L３．０２８５７２－２３２００５Ｇ０３Ｇ０４３７．５５０１１５．０２０M L３．０２５８６６－５２００５Ｇ０３Ｇ１８３７．４２０１１４．９３０M L２．５１８１８５－１３５２００５Ｇ０４Ｇ２７３７．４５０１１４．９８０M L３．２４１８７１５３２００５Ｇ０５Ｇ２０３６．９５０１１４．５３０M L２．７２５３４２１７７２００５Ｇ０６Ｇ２５３７．４５０１１４．３００M L２．６６９４２１６１２００５Ｇ０８Ｇ０５３７．５２０１１５．０３０M L３．０５８４５－４２００６Ｇ０４Ｇ２７３７．５２０１１５．０２０M L２．５２９３７９－４３２００６Ｇ０５Ｇ０６３７．５００１１４．９７０M L３．４１２７９１６０２００６Ｇ１１Ｇ２２３７．３５０１１４．６８０M L２．８６６５４０２００６Ｇ１２Ｇ２８３７．６８０１１５．２７０M L３．２２９０７７１４９２００７Ｇ０１Ｇ０５３７．６３０１１４．２２０M L２．９３４３５５－７２００７Ｇ０１Ｇ２０３７．２５０１１４．８３０M L２．５２３５３２－３２２００７Ｇ０２Ｇ０７３７．３２０１１４．８２０M L２．５８０４１１１２００７Ｇ０２Ｇ２５３７．５２０１１４．９７０M L２．５４２８８－２６２００７Ｇ０３Ｇ０９３７．５２０１１５．０３０M L２．５２９５８８１３７２００７Ｇ０８Ｇ０４３７．５５０１１５．０２０M L２．７１２８９１７９２００７Ｇ０９Ｇ２０３７．３００１１４．８００M L２．８６６２８６２００７Ｇ０９Ｇ２９３７．１５０１１４．７２０M L３．３７１７２－２２００７Ｇ０９Ｇ３０３７．１３０１１４．６３０M L２．９２５７７３－６２００７Ｇ１０Ｇ１３３７．４７０１１５．０８０M L２．５１３１３７－６２２００７Ｇ１０Ｇ１８３７．１２０１１４．６５０M L３．０７４５６０２００８Ｇ０５Ｇ３１３７．５００１１４．９８０M L２．６２０１７８－２０２００８Ｇ０８Ｇ０９３７．３８０１１４．８７０M L２．７５９５６－１７３２００８Ｇ１２Ｇ３１３７．５００１１４．９５０M L３．１３５８０１６２２００９Ｇ０１Ｇ０３３７．５８０１１５．０３０M L２．６２６２５８７２００９Ｇ０１Ｇ１４３７．５２０１１５．０２０M L２．５３１６３３－１７３２００９Ｇ０１Ｇ１５３７．３３０１１５．０３０M L２．７１９３９０１７９２００９Ｇ０２Ｇ１６３７．５３０１１５．１５０M L２．９１７７５７３０２００９Ｇ０５Ｇ１０３７．３２０１１４．８８０M L２．７８０６０－１６９２００９Ｇ０５Ｇ３１３７．２２０１１５．３２０M L２．５３１１６０１６６２０１０Ｇ０３Ｇ０１３６．９８０１１４．４２０M L３．３２４０７３－１７９２０１０Ｇ０３Ｇ３１３７．４００１１４．９５０M L２．６２６５７９６２０１０Ｇ０４Ｇ３０３７．５０６１１５．０３７M L３．９１９０８６１７４２０１０Ｇ０８Ｇ０５３７．１３０１１４．５００M L２．７７０７０－９０２０１０Ｇ１１Ｇ３０３７．４８０１１４．７２０M L２．６２８３８０１５２０１１Ｇ０１Ｇ１４３７．５３０１１５．０３０M L２．８３１８７２３４２０１１Ｇ０３Ｇ１９３７．３００１１４．８５０M L２．９１１４６８－３０２０１１Ｇ０４Ｇ１６３７．７２０１１３．７８０M L２．５１８５８７－１１７２０１１Ｇ０５Ｇ０４３７．５０６１１５．０３３M L２．８１７７８６－１５７２０１１Ｇ０９Ｇ０１３７．７３３１１５．２８３M L３．７１２０８８－２５２０１１Ｇ１２Ｇ２８３８．０２０１１５．５８０M L２．６１８８８８－３６２０１２Ｇ０３Ｇ０９３７．６９９１１２．５２０M L３．２９４５６－８２２０１２Ｇ０３Ｇ２４３７．４７０１１４．６４０M L２．８２９４８７－９２２０１２Ｇ０４Ｇ２８３７．４００１１４．２９５M L２．８９７７８－３０２０１２Ｇ０５Ｇ１２３７．４００１１４．３０４M L２．６６８１－１７０２０１２Ｇ０５Ｇ２０３７．４０４１１４．９５０M L２．５３４０７０９０２０１２Ｇ０９Ｇ０８３７．５１１１１５．０２３M L２．６１３１８９２２０１２Ｇ０９Ｇ１１３７．５６５１１５．３６７M L２．７３４３８７２７２０１２Ｇ１１Ｇ０５３７．０４４１１４．８５１M L２．９５４５１－１１８２０１３Ｇ０１Ｇ０２３７．８８２１１３．０７１M L２．６２６２８０１６２０１３Ｇ０１Ｇ１１３７．３９４１１４．９２６M L２．５２２９６２－１１３续表２发震时间纬度/(ʎ)经度/(ʎ)震级走向/(ʎ)倾角/(ʎ)滑动角/(ʎ)２０１３Ｇ０５Ｇ２５３７．６３２１１４．２６３M L ２．６３０５８９－３８２０１３Ｇ０６Ｇ０３３７．６２９１１４．２６７M L ３．２７４９－１７３２０１３Ｇ１０Ｇ１７３７．５１１１１５．０６６M L ２．６６８５－２８２０１３Ｇ１２Ｇ２６３７．８１３１１４．０７７M L ２．６２６９６７－５０２０１４Ｇ０６Ｇ１０３７．９９２１１５．７６６M L ２．８２９４８９－２４２０１４Ｇ１０Ｇ１１３７．４５４１１３．６９３M L ２．５１４８８０－１００２０１４Ｇ１１Ｇ２７３８．０３１１１５．５０４M L ３．１２８３８２３５２０１５Ｇ０１Ｇ２６３７．４４７１１４．９６６M L ２．５１３０９０－１７２０１５Ｇ１０Ｇ０６３８．０２７１１５．５５１M L ２．９１９６８０－２１２０１５Ｇ１１Ｇ０２３７．４７７１１５．０５２M L ３．２３５８４２１５５２０１６Ｇ０５Ｇ０９３７．６６３１１３．６６８M L ２．７１９２８５－１７６２０１６Ｇ０７Ｇ０９３７．５６２１１４．１３１M L ２．８１８５８９－１７６２０１６Ｇ０９Ｇ２２３７．４４９１１５．３８６M L ２．７２２４４１３６２０１７Ｇ０９Ｇ０４３７．５００１１４．３５０M L ４．４９２８５－１２２０１７Ｇ０９Ｇ１５３７．４６０１１５．０１０M L ３．５４１８３－１７８２０１８Ｇ１０Ｇ３０３７．２１０１１４．６６M L ３．５１２７８９－５９２０１９Ｇ０１Ｇ１５３７．６８０１１４．５４M L ３．５２０７８３１７８２０２０Ｇ０２Ｇ２３３８．２６０１１３．９７M L ３．０８７６－１７５２０２１Ｇ０７Ｇ１５３７．５８０１１５．１M L ４．３１７１８０１３５图１㊀邢台地震序列地质构造背景及震源机制解分布F i g ．１㊀G e o l o g i c a l s t r u c t u r e b a c k gr o u n d a n dd i s t r i b u t i o n o f f o c a l m e c h a n i s ms o l u t i o n s o fX i n g t a i e a r t h q u a k e s e qu e n c e ４个,占总体数量的３．１３％[图２(a )].将地震序列震源机制解进行总体分析,得到邢台老震区地震序列的总体震源机制[图２(a )中 Δ 处]节面１的走向㊁倾角和滑动角分别为１９０．３４ʎ,８６．８８ʎ和－１７９．９６ʎ,节面２的走向㊁倾角和滑动角分别为１００．３４ʎ,８９．９６ʎ和－３．１２ʎ,属于走滑型地震,具体见图２(b).１．２㊀震源机制节面聚类方法震源机制节面聚类方法是一种独立于地质资料㊁大地测量数据和小震定位结果的空间聚类方法.相比较挖探槽法,该方法不存在深部断裂和浅部断图２㊀邢台地震序列的震源机制分类及总体特征F i g．２㊀C l a s s i f i c a t i o na n do v e r a l l c h a r a c t e r i s t i c s o f f o c a l m e c h a n i s m s o fX i n g t a i e a r t h q u a k e s e qu e n c e 裂的差异;相比较地震反射探测法,该方法无需耗费大量人力物力,且不存在地球物理资料反演的多解性问题;相比较余震震源位置确定断层面几何参数３４４１第４５卷第６期㊀㊀㊀㊀㊀㊀㊀㊀㊀章㊀阳,等:基于震源机制解分析邢台地震的发震构造㊀㊀㊀㊀㊀㊀㊀㊀㊀方法,该方法不依赖余震的丛集展布,且对地震定位精度要求不高.万永革[２０]将震源机制节面聚类方法应用于２０２１年云南漾濞地震序列中,得到了与该地震序列的地震分布走向大体一致的两个分支断裂带,验证了该种方法的有效性,同时给出结果的置信区间.随着数字地震台网的不断加密,更多中小地震震源机制资料的获取可以进一步丰富该方法的扩展应用.根据震源机制两个节面的法向矢量,求解多个震源机制法向之间的夹角,对多个震源机制的多个节面按照其法向夹角作为距离进行聚类分析,聚类完成后求得与活动断层走向较为一致的类,其中心就包含活动断层面的几何参数,通过计算节面与类中心夹角的标准差,估计活动断层几何参数的标准差.具体计算方法参考万永革[２０]一文,对震源机制的节面进行聚类分析,提取断层几何信息,得到断层面的走向和倾角,并给出其置信区间.１．３㊀活动断层几何参数的确定采用节面聚类方法对收集到的地震序列震源机制(１２８个震源机制,２５６个节面)进行聚类,可以得到两簇聚类中心.第一簇节面数为２０,其中心节面法向的走向为１０５．９０ʎ,倾伏角为１３．９７ʎ,标准差为１１．１２ʎ,其中心节面的断层面走向为１９５．９０ʎ,置信区间为１９０．９１ʎ~２００．８９ʎ,倾角为７６．０３ʎ,置信区间为７１．１６ʎ~８０．９１ʎ;第二簇节面数为３５,其断层面法向的走向为２００．４２ʎ,倾伏角为３．１３ʎ,标准差为１８．８１ʎ,其中心节面的断层面走向为２９０．４２ʎ,置信区间为２８３．２９ʎ~２９７．５５ʎ,倾角为８６．８７ʎ,置信区间为７７．６１ʎ~８７．１３ʎ;第一类和第二类的类中心节面的夹角为８６．３８ʎ,接近垂直.噪声节面数据个数为７３个,并且在全空间分布,较为随机.具体结果见图３.震源机制节面聚类方法得到两个聚类面,断层面１走向１９５．９０ʎ,倾角７６．０３ʎ,断层面２走向２９０．４２ʎ,倾角８６．８７ʎ,结合震中附近地质构造特征,判断断层面１为发震破裂面,断层面２为破裂辅助面.震源机制节面聚类方法独立于地质资料㊁大地测量数据和小震定位结果,并且计算结果与实际地质构造相符,与前人关于邢台地震发震构造的研究结果也保持一致,进而显示出该方法的有效性和优越性.注:震源机制节面采用绿色弧线表示,聚类中心节面采用红色弧线表示.黑点表示震源机制节面的极点位置,红点表示聚类中心的极点位置.聚类中心极点周围的蓝绿色椭圆为聚类中心的误差范围图３㊀邢台地震序列震源机制聚类结果F i g．３㊀C l u s t e r i n g r e s u l t s o f f o c a lm e c h a n i s m s o fX i n g t a i e a r t h q u a k e s e q u e n c e２㊀邢台地震断层面平均解将邢台地震序列震源机制总体特征㊁震源机制节面聚类结果㊁前人通过不同资料和方法得到邢台地震断层面结果以及７．２级地震断层面解进行比较(表２).发现断层面走向有一定的差异性,但整体显示出较好的一致性,倾角基本一致.为综合分析邢台地震的发震构造,根据表２的１２组断层面结果计算断层面平均解,具体结果见图４.最终得到邢台地震断层面走向为３２．４５ʎ,误差范围３５７．１５ʎ~６７．７６ʎ,倾角为７９．４４ʎ,误差范围６６．３３ʎ~８７．４５ʎ.该断层为N N E走向的高倾角型断层.断层面平均解与震源机制节面聚类结果有一定的差异,造成差异的可能原因是地震震源机制解数据的不确定性,会导致聚类结果的不稳定,早期的地震序列由于观测资料和研究方法的局限,震源机制解结果的精度有待商榷,而且大地震破裂后会导致周边应力场发生较大的变化,使得震源机制辅助面有较大差异.３㊀邢台地震断层面滑动性质为判断断层的滑动性质,我们结合区域构造应力场参数求解该断层的滑动角.首先根据收集的震４４４１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀地㊀震㊀工㊀程㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０２３年表２㊀邢台地震的断层面几何参数T a b l e２㊀G e o m e t r i c p a r a m e t e r s o f t h e f a u l t p l a n e f o rX i n g t a i e a r t h q u a k e序号结果来源走向/(ʎ)倾角/(ʎ)１陈运泰等断层模式断层面(出露到地面)[１]３５４５２陈运泰等断层模式断层面(地面以下)[１]３５８２３李钦祖等综合断层面解节面[２７]１８７８４７．２级地震断层面解(中国震源机制研究小组[２８])９３６６５７．２级地震断层面解(李钦祖等[２１])２７９８７６７．２级地震断层面解(鄢家全等[２９])１９９０７７．２级地震断层面解(张诚等[３０])２９７８８７．２级地震断层面解(叶洪等[３１])１１８５９７．２级地震断层面解(张之立等[３２])２５７９１０７．２级地震断层面解(张诚等[３０])２６７３１１本文震源机制节面聚类结果１９５．９０７６．０３１２地震序列总体特征１９０．３４８６．８８１３断层面平均解３２．４５７９．４４注:黑色虚线为表２断层面解,红色实线为断层面平均解.红色五角星为断层面平均解的中心点图４㊀邢台地震发震断层面平均解F i g．４㊀A v e r a g e s o l u t i o no f t h e f a u l t p l a n e o fX i n g t a i e a r t h q u a k e源机制解反演区域构造应力场[３３Ｇ３４],按照应力张量在断层面上的剪应力方向与断层滑动方向最为一致准则,联合采用擦痕数据和定性断层滑动数据,读取搜索空间间隔和置信度,采用网格搜索法求解应力场,并给出反演参数的一定置信度下的置信区间.求解得到最优应力状态的压轴走向６０．０４ʎ,倾伏角２４．４９ʎ,９５％置信度置信范围分别为５９．３２ʎ~６５．０４ʎ和１６．９９ʎ~２８．２１ʎ;张轴走向１５３．８０ʎ,倾伏角８．２１ʎ,９５％置信度置信范围分别为１５３．３０ʎ~１５８．８０ʎ和４．２１ʎ~１６．２１ʎ;中间轴走向２６１．００ʎ,倾伏角６４．００ʎ,９５％置信度置信范围分别为２６０．５０ʎ~２６６．００ʎ和６０．００ʎ~７２．００ʎ,应力形因子R为０．２,主压应力轴方向为N E E向,主张应力轴方向为NWW向.将地壳应力张量投影到断层面上,最终求得断层面的滑动角为－１５３．９６ʎ,标准差为４．５５ʎ.该断层为走向N N E㊁倾向NWW的高倾角走滑型断层,这与收集的邢台老震区地震序列的总体震源机制一致.本研究得到的这条走滑型断层可以为该地区的地质构造和地震孕育环境分析提供更丰富的资料.４㊀讨论前人利用深地震反射剖面揭示邢台老震区地壳内存在一个从高到低角度的铲状断裂,上地壳到中㊁下地壳存在一个沿NWWＧS E E方向伸展的滑脱面.同时根据反射界面的连续性和倾向,推断存在一条高倾角的走滑性质的超壳断裂,认为这可能是邢台地震的发震断层.本文结合多种数据综合分析,最后确定邢台地震的发震构造为走向N N E㊁倾向NWW方向的高倾角走滑型断层.把浅层的正断性质的新河断裂作为发震断层是不合适的,这也进一步论证前人研究结果.为进一步分析该断层的地震活动危险性,我们按照万永革[３５]提出的震源机制与应力体系关系模拟方法计算该应力张量下的各种形状的断层面上的相对剪应力和相对正应力分布(图５).将震源机制节面的剪应力和正应力用最大剪应力归一化即得到相对剪应力和相对正应力,图５(a)为相对剪应力分布,最大剪应力在走向１００ʎ㊁倾角７０ʎ和走向２００ʎ㊁倾角７０ʎ的断层面上,最小剪应力在走向１５０ʎ㊁倾角７０ʎ和走向３５０ʎ㊁倾角３０ʎ的断层面上.图５(b)为相对正应力分布,最大正应力在走向６０ʎ,倾角９０ʎ和走向２４０ʎ,倾角８０ʎ的断层面上,最小正应力在走向１５０ʎ㊁倾角７０ʎ和走向３３０ʎ㊁倾角９０ʎ的断层面上.区域地壳应力场在该断层面上的相对剪应力为０．８３,相对正应力为０．５９,该断层不是最大剪应力的断层面形状.虽然这条断层与反演的现今应力场的最大释放应力方向不一致,但是剪应力强度仍然比较大,该断层在地震能量积累和释放过程中仍起到必不可少的作用.后续这条断层还将会是该区域的一条主要发震断层,对于研判区域地震活动危险性具有一定的指示意义.５㊀结论本文根据独立于地震波资料㊁大地测量资料和地质资料等之外的确定活动断层几何形状的震源机制节面聚类方法,收集邢台及周边区域的震源机制解资料,确定活动断层的走向和倾角,并给出其标准５４４１第４５卷第６期㊀㊀㊀㊀㊀㊀㊀㊀㊀章㊀阳,等:基于震源机制解分析邢台地震的发震构造㊀㊀㊀㊀㊀㊀㊀㊀㊀差.与前人研究结果对比,综合求取发震断层面的平均解,结合区域地壳应力张量,估计该断层的滑动方向及其误差,进而分析该区域的地震活动危险性.主要获得以下结论:注:N S 为正走滑型,S S 为走滑型,N 为正断型,R S 为逆走滑型,R 为逆断型图５㊀区域应力场在各种形状断层面上的相对剪应力和相对正应力F i g ．５㊀R e l a t i v e s h e a r s t r e s s a n d r e l a t i v en o r m a l s t r e s s o f t h e r e g i o n a l s t r e s s f i e l do nv a r i o u s s h a pe s of f a u l t p l a n e ㊀㊀(１)邢台地震发震断层的走向３２．４５ʎ,倾角７９．４４ʎ,滑动角－１５３．９６ʎ,标准差为４．５５ʎ.该断层为走向N N E ㊁倾向NWW 的高倾角走滑型断层,而非正断性质的新河断裂.该断层几何形状的确定对于后续该区域的地震孕育环境及地震发震构造具有很好的参考价值.(２)计算区域地壳应力场在断层面上产生的相对剪应力为０．８３,相对正应力为０．５９.该断层不是最大剪应力的断层面形状,但应力场在这条断层上的剪应力强度仍较大,该断层在地震能量积累和释放过程中仍起到必不可少的作用.后续这条断层还将会是该区域的一条主要发震断层,对于研判区域地震活动危险性具有一定的指示意义.参考文献(R e f e r e n c e s)[１]㊀陈运泰,林邦慧,林中洋,等．根据地面形变的观测研究１９６６年邢台地震的震源过程[J ]．地球物理学报,１９７５,１８(３):１６４Ｇ１８２．C H E N Y u n t a i ,L I N B a n g h u i ,L I N Z h o n g y a n g ,e ta l ．T h e f o c a l m e c h a n i s mo f t h e１９６６H s i n g t a i e a r t h qu a k ea s i n f e r r e df r o m t h e g r o u n dd e f o r m a t i o no b s e r v a t i o n s [J ]．C h i n e s eJ o u r n a lo f S i n i c a ,１９７５,１８(３):１６４Ｇ１８２．[２]㊀徐杰,方仲景,杨理华．１９６６年邢台７．２级地震的构造背景和发震构造[J ]．地震地质,１９８８,１０(４):５１Ｇ５９．６４４１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀地㊀震㊀工㊀程㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０２３年X UJ i e,F A N G Z h o n g j i n g,Y A N G L i h u a．T e c t o n i cb a c k g r o u n da n d c a u s a t i v e f a u l t o f１９６６X i n g t a i M S７．２e a r t h q u a k e[J]．S e i sＧm o l o g y a n dG e o l o g y,１９８８,１０(４):５１Ｇ５９．[３]㊀赵成斌,孙振国,刘保金,等．邢台地震区浅部构造特征及其与深部构造的耦合关系[J]．地震地质,１９９９,２１(４):４１７Ｇ４２４．Z H A OC h e n g b i n,S U NZ h e n g u o,L I U B a o j i n,e t a l．As t u d y o f t h e s h a l l o ws t r u c t u r a l c h a r a c t e r i s t i c s a n d t h e c o u p l i n g r e l a t i o nb e t w e e nd e e p a n d s h a l l o ws t r uc t u r e s i nX i n g t a i e a r t h q u a k e a rＧe a[J]．S e i s m o l o g y a n dG e o l o g y,１９９９,２１(４):４１７Ｇ４２４．[４]㊀徐锡伟,于贵华,王峰,等．１９６６年邢台地震群的发震构造模型:新生断层形成?先存活断层摩擦粘滑?[J]．中国地震,２０００,１６(４):３６４Ｇ３７８X U X i w e i,Y U G u i h u a,WA N G F e n g,e t a l．S e i s m o g e n i cm o d e lf o r t h e１９６６X i ng t a i e a r th q u a k e s:n u c l e a ti o no f n e wＧb o r n f a u l to r s t r i c kＧs l i p o f p r eＧe x i s t i n g f a u l t[J]．E a r t h q u a k eR e s e a r c h i nC h i n a,２０００,１６(４):３６４Ｇ３７８．[５]㊀邵志刚,李旖雯,王芃．１９６６年邢台强震群震前㊁同震和震后断层运动特征分析[J]．地震,２０１５,３５(３):１Ｇ９．S H A OZ h i g a n g,L IY i w e n,WA N G P e n g．P r eＧ,c oＧ,a n d p o s tＧs e i s m i c f a u l t m o v e m e n tf e a t u r e so ft h e１９６６X i n g t a is t r o n ge a r t h q u a k e s w a r m[J]．E a r t h q u a k e,２０１５,３５(３):１Ｇ９．[６]㊀王椿镛,段永红,吴庆举,等．华北强烈地震深部构造环境的探测与研究[J]．地震学报,２０１６,３８(４):５１１Ｇ５４９．WA N G C h u n y o n g,D U A N Y o n g h o n g,WU Q i n g j u,e ta l．E xＧp l o r a t i o no nt h ed e e p t e c t o n i ce n v i r o n m e n to fs t r o n g e a r t hＧq u a k e s i nN o r t hC h i n a a n dr e l e v a n t r e s e a r c h f i n d i n g s[J]．A c t a S e i s m o l o g i c aS i n i c a,２０１６,３８(４):５１１Ｇ５４９．[７]㊀王椿镛,吴庆举,段永红,等．华北地壳上地幔结构及其大地震深部构造成因[J]．中国科学(地球科学),２０１７,４７(６):６８４Ｇ７１９．WA N GC h u n y o n g,WU Q i n g j u,D U A N Y o n g h o n g,e t a l．C r u sＧt a l a n d u p p e rm a n t l e s t r u c t u r e i nN o r t hC h i n a a n d i t s d e e p t e cＧt o n i c g e n e s i s i n g r e a t e a r t h q u a k e[J]．S c i e n t i aS i n i c a(T e r r a e),２０１７,４７(６):６８４Ｇ７１９．[８]㊀X U X W,W E NXZ,Y U G H,e t a l．C o s e i s m i c r e v e r s eＧa n d o bＧl i q u eＧs l i p s u r f a c e f a u l t i n gg e n e r a t e db y t h e２００８M W７．９W e nＧc h u a ne a r t h q u a k e,C h i n a[J]．G e o l o g y,２００９,３７(６):５１５Ｇ５１８．[９]㊀张培震,闵伟,邓起东,等．海原活动断裂带的古地震与强震复发规律[J]．中国科学(D辑:地球科学),２００３,３３(８):７０５Ｇ７１３．Z H A N GP e i z h e n,M I N W e i,D E N G Q i d o n g,e ta l．P a l e o e a r t hＧq u a k e sa n dr e c u r r e n c e l a w o f s t r o n g e a r t h q u a k e s i n H a i y u a na c t i v e f a u l t z o n e[J]．S c i e n t i a S i n i c a(T e r r a e),２００３,３３(８):７０５Ｇ７１３．[１０]㊀冉勇康,王虎,陈立春,等．龙门山断裂带晚第四纪的大地震活动:来自古地震研究的资料[J]．地球物理学报,２０１８,６１(５):１９３８Ｇ１９４８．R A N Y o n g k a n g,WA N G H u,C H E N L i c h u n,e t a l．L a t eＧQ u aＧt e r n a r y f a u l ta c t i v i t y o ft h eL o n g m e nS h a nf a u l tz o n e:e v iＧd e n c e f r o m p a l e o s e i s m i c t r e n c h i n g[J]．C h i n e s e J o u r n a l o f G e oＧp h y s i c s,２０１８,６１(５):１９３８Ｇ１９４８．[１１]㊀刘保金,赵成彬,酆少英,等．应用三分量浅层地震反射方法探测隐伏活动断裂[J]．地球物理学报,２０１２,５５(８):２６７６Ｇ２６８６．L I U B a o j i n,Z HA OC h e n g b i n,F E N GS h a o y i n g,e t a l．A p p l i c aＧt i o n o f t h e t h r e eＧc o m p o n e n t s h a l l o ws e i s m i c r e f l e c t i o nm e t h o dt o p r o b i n g b u r i e da c t i v ef a u l t s[J]．C h i n e s eJ o u r n a lo fG e oＧp h y s i c s,２０１２,５５(８):２６７６Ｇ２６８６．[１２]㊀李燕,刘保金,酆少英,等．利用地震折射和反射波资料研究银川盆地浅部结构和隐伏断裂[J]．地球物理学报,２０１７,６０(８):３０９６Ｇ３１０９．L IY a n,L I U B a o j i n,F E N G S h a o y i n g,e ta l．E x p l o r a t i o n o fs h a l l o w s t r u c t u r ea n d b u r i e df a u l t si nt h e Y i n c h u a n B a s i nu s i n g s e i s m i c r e f r a c t i o na n dr e f l e c t i o nd a t a[J]．C h i n e s eJ o u rＧn a l o fG e o p h y s i c s,２０１７,６０(８):３０９６Ｇ３１０９．[１３]㊀徐增波,刘保金,姬计法,等．太行山南端浅层速度结构成像和隐伏断裂探测[J]．大地测量与地球动力学,２０１９,３９(１):８８Ｇ９２．X UZ e n g b o,L I U B a o j i n,J I J i f a,e t a l．I m a g i n g f o r s h a l l o wv eＧl o c i t y s t r u c t u r ea n db u r i e df a u l t s i nt h es o u t h e r n m a r g i no fT a i h a n g M o u n t a i n s[J]．J o u r n a l o fG e o d e s y a n dG e o d y n a m i c s,２０１９,３９(１):８８Ｇ９２．[１４]㊀万永革,沈正康,刁桂苓,等．利用小震分布和区域应力场确定大震断层面参数方法及其在唐山地震序列中的应用[J]．地球物理学报,２００８,５１(３):７９３Ｇ８０４．WA N Y o n g g e,S H E NZ h e n g k a n g,D I A O G u i l i n g,e t a l．A n a lＧg o r i t h mo f f a u l t p a r a m e t e r d e t e r m i n a t i o n u s i n g d i s t r i b u t i o n o fs m a l l e a r t h q u a k e s a n d p a r a m e t e r s o f r e g i o n a l s t r e s s f i e l da n di t s a p p l i c a t i o n t oT a n g s h a ne a r t h q u a k es e q u e n c e[J]．C h i n e s eJ o u r n a l o fG e o p h y s i c s,２００８,５１(３):７９３Ｇ８０４．[１５]㊀杨超群,孟凡顺,万永革．采用精确定位小震资料反演伽师地震断层面[J]．地球物理学进展,２０１３,２８(６):２８６５Ｇ２８７１．Y A N GC h a o q u n,M E N GF a n s h u n,WA N Y o n g g e．C a l c u l a t i n gt h e p a r a m e t e r s o f f a u l t p l a n e f o r J i a s h i e a r t h q u a k e b y t h e d a t ao f s m a l l e a r t h q u a k e sa c c u r a t e l y l o c a t e d[J]．P r o g r e s s i nG e oＧp h y s i c s,２０１３,２８(６):２８６５Ｇ２８７１．[１６]㊀盛书中,万永革,王未来,等．２０１０年玉树M S７．１地震发震断层面参数的确定[J]．地球物理学进展,２０１４,２９(４):１５５５Ｇ１５６２．S H E N GS h u z h o n g,WA N Y o n g g e,WA N G W e i l a i,e ta l．T h ef a u l t p l a n e p a r a m e t e r d e t e r m i n a t i o n o f t h e２０１０Y u s h u M S７．１e a r t h q u a k e[J]．P r o g r e s si n G e o p h y s i c s,２０１４,２９(４):１５５５Ｇ１５６２．[１７]㊀潘睿,盛书中,万永革．２００９年姚安地震断层面参数的确定[J]．华北地震科学,２０１５,３３(２):６３Ｇ６８．P A N R u i,S H E N G S h u z h o n g,WA N Y o n g g e．S e i s m i cf a u l tp l a n e p a r a m e t e r s e s t i m a t i n g o f２００９Y a o＇a ne a r t h q u a k e[J]．N o r t hC h i n aE a r t h q u a k eS c i e n c e s,２０１５,３３(２):６３Ｇ６８．[１８]㊀胡晓辉,盛书中,万永革,等．由余震分布确定２０１７年九寨沟M S７．０地震发震断层面参数[J]．地震研究,２０１９,４２(３):３６６Ｇ３７１,４５５．HU X i a o h u i,S H E N G S h u z h o n g,WA N Y o n g g e,e ta l．F a u l tp l a n e p a r a m e t e r s o f２０１７J i u z h a i g o u M S７．０e a r t h q u a k e d e t e rＧm i n e db y a f t e r s h o c kd i s t r i b u t i o n[J]．J o u r n a l o fS e i s m o l o g i c a lR e s e a r c h,２０１９,４２(３):３６６Ｇ３７１,４５５．[１９]㊀章阳,万永革,陈婷,等．河北临城地震序列揭示的一条隐伏断层[J]．中国地震,２０２１,３７(２):４４２Ｇ４５１．７４４１第４５卷第６期㊀㊀㊀㊀㊀㊀㊀㊀㊀章㊀阳,等:基于震源机制解分析邢台地震的发震构造㊀㊀㊀㊀㊀㊀㊀㊀㊀Z HA N G Y a n g,WA N Y o n g g e,C H E N T i n g,e ta l．A b l i n df a u l t r e v e a l e db y t h eL i n c h e ng M L４．４e a r th q u a k e s e q u e n c e o fH e b e i P r o v i n c e[J]．E a r t h q u a k eR e s e a r c hi n C h i n a,２０２１,３７(２):４４２Ｇ４５１．[２０]㊀万永革．断裂带震源机制节面聚类确定断裂带产状方法及在２０２１年漾濞地震序列中的应用[J]．地球物理学报,２０２２,６５(２):６３７Ｇ６４８．WA N Y o n g g e．M e t h o do f a c t i v e f a u l t g e o m e t r y d e t e r m i n a t i o nb yc l u s t e r i n g n od a l p l a ne sof f o c a lm e c h a n i s m so c c u r r e do nt h ef a u l tb e l ta n di t sa p p l i c a t i o nt ot h e２０２１Y a n g b i e a r t hＧq u a k es e q u e n c e[J]．C h i n e s eJ o u r n a lo fG e o p h y s i c s,２０２２,６５(２):６３７Ｇ６４８．[２１]㊀李钦祖,靳雅敏,于新昌．华北地区的震源机制与地壳应力场[J]．地震学报,１９８２,４(１):５５Ｇ６１．L IQ i n z u,J I N Y a m i n,Y U X i n c h a n g．F o c a lm e c h a n i s m sa n dc r u s t a l s t r e s s f i e ld i nN o r t hC h i n a[J]．A c t a Se i s m o l o g i c a S i n iＧc a,１９８２,４(１):５５Ｇ６１．[２２]㊀许忠淮,汪素云,高阿甲,等．我国部分早期震源机制解答的重新测定[J]．地震地磁观测与研究,１９９４,１５(５):１Ｇ９．X UZ h o n g h u a i,WA N GS u y u n,G A O A j i a,e t a l．R e d e t e r m i n aＧt i o no fs o m e e a r l y f o c a l m e c h a n i s m s o l u t i o n s o f C h i n e s ee a r t h q u a k e s[J]．S e i s m o l o g i c a l a n dG e o m a n g n e t i cO b s e r v a t i o na n dR e s e a r c h,１９９４,１５(５):１Ｇ９．[２３]㊀王晓山,陈婷,冯向东,等．邢台地震序列的应力释放调整过程[J]．华北地震科学,２０１７,３５(１):１７Ｇ２３．W A N G X i a o s h a n,C H E NT i n g,F E N GX i a n g d o n g,e t a l．S t r e s s r eＧl e a s i n g a n d a d j u s t i n gp r o g r e s s o fX i n g t a i e a r t h q u a k e s e q u e n c e[J]．N o r t hC h i n aE a r t h q u a k e S c i e n c e s,２０１７,３５(１):１７Ｇ２３．[２４]㊀王晓山,冯向东,赵英萍．京津冀地区地壳应力场特征[J]．地震研究,２０２０,４３(４):６１０Ｇ６１９．WA N G X i a o s h a n,F E N G X i a n g d o n g,Z HA O Y i n g p i n g．C h a rＧa c t e r i s t i c so fc r u s t a ls t r e s sf i e l di n B e i j i n g T i a n j i n H eb e ir e g i o n[J]．J o u r n a lo fS e i s m o l o g i c a lR e s e a r c h,２０２０,４３(４):６１０Ｇ６１９．[２５]㊀陈婷,章阳,王晓山,等．邢台老震区库仑应力演化及地震危险性分析[J]．华北地震科学,２０２１,３９(４):５０Ｇ５６．C H E N T i n g,Z H A N GY a n g,WA N GX i a o s h a n,e t a l．C o u l o m bs t r e s s e v o l u t i o n a n d s e i s m i c h a z a r d a n a l y s i s o n t h eX i n g t a i o l de a r t h q u a k e r e g i o n[J]．N o r t h C h i n a E a r t h q u a k e S c i e n c e s,２０２１,３９(４):５０Ｇ５６．[２６]㊀万永革．震源机制水平应变花面应变的地震震源机制分类方法及序列震源机制总体特征分析[J/O L]．地球科学,２０２２:１Ｇ１６(２０２２Ｇ０７Ｇ１８)[２０２２Ｇ０８Ｇ０６]．h t t p s://k n s．c n k i．n e t/k c m s/d eＧt a i l/４２．１８７４．p．２０２２０７１５．１５３２．０１４．h t m l．W a nY o n g g e．S e i s m i c s o u r c em e c h a n i s mc l a s s i f i c a t i o nm e t h o df o r h o r i z o n t a l s t r a i n f l o w e r s u r f a c e s t r a i n a n d a n a l y s i s o f o v e rＧa l lc h a r a c t e r i s t i c s o fs e q u e n c e s o u r c e m e c h a n i s m[J/O L]．E a r t hS c i e n c e,２０２２:１Ｇ１６(２０２２Ｇ０７Ｇ１８)[２０２２Ｇ０８Ｇ０６]．h t t p s://k n s．c n k i．n e t/k c m s/d e t a i l/４２．１８７４．p．２０２２０７１５．１５３２．０１４．h t m l．[２７]㊀李钦祖,王泽皋,贾云年,等．由单台小地震资料所得两个区域的应力场[J]．地球物理学报,１９７３,１６(１):４９Ｇ６１．L IQ i n z u,WA N GZ e g a o,J I A Y u n n i a n,e t a l．S t r e s s f i e l do bＧt a i n e d f o r t w o r e g i o n s f r o m w e a k e a r t h q u a k e d a t a r e c o r d e d a ta s i n g l e s e i s m i c s t a t i o n[J]．C h i n e s e J o u r n a l o f S i n i c a,１９７３,１６(１):４９Ｇ６１．[２８]㊀中国震源机制研究小组．中国地震震源机制研究[M]．第一集．北京:地震出版社,１９７３．C h i n e s eF o c a lM e c h a n i s m R e s e a r c hG r o u p．R e s e a r c h t h e f o c a lm e c h a n i s m s o f e a r t h q u a k e s i nC h i n a[M]．V o l u m e１．B e i j i n g:S e i s m o l o g i c a l P r e s s,１９７３．[２９]㊀鄢家全,时振梁,环文林,等．强余震的断层面解特征[J]．地震学报,１９８０,２(４):３９５Ｇ４０３．Y A NJ i a q u a n,S H I Z h e n l i a n g,HU A N W e n l i n,e t a l．T h e c h a rＧa c t e r i s t i c s o f f o c a l p l a n es o l u t i o n so fs t r o n g a f t e r s h o c k s[J]．A c t aS e i c m o l o g i c aS i n i c a,１９８０,２(４):３９５Ｇ４０３．[３０]㊀张诚．中国地震震源机制[M]．北京:学术书刊出版社,１９９０．Z H A N GC h e n g．F o c a lm e c h a n i s m o fC h i n ae a r t h q u a k e[M]．B e i j i n g:A c a d e m i cB o o k sa n dP e r i o d i c a l sP u b l i s h i n g H o u s e,１９９０．[３１]㊀叶洪,张文郁．初论华北板内地震断层特征及其与地壳㊁上地幔动力学过程的关系[J]．地震地质,１９８０,２(１):２７Ｇ３８．Y E H o n g,Z H A N G W e n y u．T h ec h a r a c t e r i s t i c so f i n t r a p l a t ee a r t h q u a k ef a u l t s i nN o r t hC h i n a a n d t h e i r r e l a t i o n s h i p t o t h ed y n a m i c a l p r o ce s s e s i ne a r t h＇sc r u s ta n du p p e r m o s t m a n t l e[J]．S e i s m o l o g y a n dG e o l o g y,１９８０,２(１):２７Ｇ３８．[３２]㊀张之立,李钦祖,谷继成,等．唐山地震的破裂过程及其力学分析[J]．地震学报,１９８０,２(２):１１１Ｇ１２９．Z H A N G Z h i l i,L I Q i n z u,G U J i c h e n g,e ta l．T h ef r a c t u r ep r o c e s s e s o ft h e T a n g s h a n e a r t h q u a k ea n di t s m e c h a n i c a la n a l y s i s[J]．A c t aS e i s m o l o g i c aS i n i c a,１９８０,２(２):１１１Ｇ１２９．[３３]㊀万永革．联合采用定性和定量断层资料的应力张量反演方法及在乌鲁木齐地区的应用[J]．地球物理学报,２０１５,５８(９):３１４４Ｇ３１５６．WA N Y o n g g e．A g r i d s e a r c hm e t h o d f o r d e t e r m i n a t i o no f t e cＧt o n i c s t r e s st e n s o ru s i n gq u a l i t a t i v ea n d q u a n t i t a t i v ed a t ao fa c t i v e f a u l t sa n d i t sa p p l i c a t i o nt ot h eU r u m q i a r e a[J]．C h iＧn e s e J o u r n a l o fG e o p h y s i c s,２０１５,５８(９):３１４４Ｇ３１５６．[３４]㊀WA N Y G,S H E N GSZ,HU A N GJC,e t a l．T h e g r i ds e a r c ha l g o r i t h mo f t e c t o n i c s t r e s s t e n s o rb a s e do n f oc a lm e c h a n i s md a t a a n d i t sa p p l i c a t i o n i nt h eb o u n d a r y z o n eo fC h i n a,V ie tＧn a ma n dL a o s[J]．J o u r n a l o fE a r t hS c i e n c e,２０１６,２７(５):７７７Ｇ７８５．[３５]㊀万永革．震源机制与应力体系关系模拟研究[J]．地球物理学报,２０２０,６３(６):２２８１Ｇ２２９６．WA N Y o n g g e．S i m u l a t i o no nr e l a t i o n s h i p b e t w e e ns t r e s sr eＧg i m e s a n d f o c a lm e c h a n i s m s o f e a r t h q u a k e s[J]．C h i n e s e J o u rＧn a l o fG e o p h y s i c s,２０２０,６３(６):２２８１Ｇ２２９６．(本文编辑:任㊀栋)８４４１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀地㊀震㊀工㊀程㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０２３年。
震源机制解
三、P波的四象限分布 在地震学的早期研究中，人们就 已注意到P波到达时地面的初始振动有 时是向上的,有时是向下的。20世纪的 10～20年代，许多地震学者在日本和 欧洲的部分地区几 乎同时发现，同一 次地震在不同地点 的台站记录,所得的 P波初动方向具有四 象限分布特征。
地表垂直向地震仪记录P震相的初始 振动方向:向上的,记为正号；向下的，记 为负号。正号P波是压缩波,因为这种波的 到达使台站受到来自地下的一个突然挤压， 台基介质体积发
十一、求解方法
1、P波初动法 2、P波S波振幅比方法 （根据震源理论，在均匀各向同性介质 中，双力偶震源幅射的远场地震波位移 分量表达式，利用振幅的相对大小，可 算出P、T、B轴方位）3、其它方法（例 如在频率域分析P波或S的振幅谱低频幅 值等）
M
上式表述的是双力偶震源辐射的远场地震波位移在观测点 P(r ,θ,Φ)处的分量表达式,ρ为岩石密度,r表达位移的点至震源 的距离, 是双力偶中一个力偶强度随时间的微商 M P波与S波振幅比反演震源机制的方法 振幅比的幅射花样随空间方位的变化比单种波的辐射花样要强 烈得多.从此意义来说,只要有正确的观测振幅比,且观测值归算 至震源球面上后的位置正确的话,振幅比观测对震源机制参数有 较强的约束力.但用振幅比测定震源机制解时,只能求出两节面 的空间位置,不能确定可能断层面的运动特性,还须借助至少一 个P波初动.
近震震源球示意图
远震震源球示意图
只要记录足够多，且台站对应点在震源 球面上的分布范围足够广，则总可找到两个 互相垂直的大圆面将震源球面上的正、负号 分成四个部分,即四象限, 这两个互相垂直的 大圆面称为 P波初动的节面，节面与地面的 交线称为节线，节面上P 波初动位移为零。二节面 之一 (AA′)与地震的断 层面一致，而另一个面
阿拉善左旗5.0级地震震源机制解与发震构造判定
II-40阿拉善左旗5.0级地震震源机制解与发震构造判定李 娟※ 张 帆 韩晓明（内蒙古自治区地震局，呼和浩特 010000）中图分类号：P315.3+3 文献标识码： A doi ：10.3969/j.issn.0253-4975.2018.08.0432017年6月3日18时11分00秒，内蒙古自治区阿拉善左旗（37.99°N ，103.56°E ）发生5.0级地震，震源深度9 km ，宏观震中位于阿拉善左旗额尔克哈什哈苏木乌尼格图嘎查，未造成明显地表破坏。
根据内蒙古地震台网测定，截止2017年6月20日，阿拉善左旗5.0级地震共发生余震92次，最大余震为6月3日18时54分M L 3.1地震。
余震的震级主要集中分布在M L 0.0—1.0之间，时间上主要集中在6月3—5日。
本次5.0地震发生在Ⅰ级构造单元西域地块的次级地块阿拉善块体西南边界。
该区域附近的主要活动断裂有河西堡—四道山断裂、海原断裂、天桥沟—黄羊川断裂、老虎山断裂等。
为了分析本次地震的震源机制解从而判定发震构造，基于内蒙古测震台网数据，利用波形拟合的矩张量反演方法（CAP 法）和HypoDD 双差定位法，并结合强震动场及烈度区划图综合进行判定发震构造，以期对该区地震研究提供依据。
选取震中距在50—350 km 范围的15个台站的波形记录，测定了阿拉善左旗5.0级地震的震源机制解，结果显示，5.0级地震的最佳双力偶震源机制解参数分别为：节面I 走向89°、倾角82°、滑动角8°；节面Ⅱ走向358°、倾角82°、滑动角172°；P 轴方位43°、仰角0°；T 轴方位313°、仰角11°。
矩心深度在8.6 km 处残差最小，即最佳反演深度为8.6 km 。
GCMT 、中国地震局地球物理研究所、台网中心均发布了此次阿拉善5.0级地震的震源机制解，其中GCMT 与中国地震局地球物理研究所的结果与本文基本一致，震源机制表现为左旋走滑特征，这与阿拉善块体应力场背景相一致。
青海门源Ms6.4地震震源机制解与震源深度研究
青海门源Ms6.4地震震源机制解与震源深度研究尹欣欣;赵林林;杨立明;陈继锋;左可桢;蒲举【期刊名称】《大地测量与地球动力学》【年(卷),期】2018(038)006【摘要】Based on waveform data recorded by 19 broadband digital stations of Gansu,Qinghai and Ningxia regional seismic network,we use the CAP method to calculate the focal mechanism solution and the optimal focal depth of Menyuan Ms6.4 earthquake in January 21,2016.The results show that the double-couple solution is:nodal Ⅰ:strike 143°,dip angle 40°,slip angle 71°;nodal Ⅱ:strike 347.2°,dip angle 52.6°,slip angle105.3°.The moment magnitude is Mw5.9 and the optimal depth is7.7km,which is consistent with the results of others (such as CENC,IGP-CEA and Harvard).Using the arrival time difference between sPn phase extracted by sliding window correlation method,we obtain the focal depth of 8.5 km,which is coincident with the result of CAP inversion,and we validate the feasibility of the method.%选取甘肃、青海和宁夏区域测震台网19个宽频带数字台站的地震波形,采用CAP震源机制解方法,研究2016-01-21青海省门源县Ms6.4地震,得到其震源机制解和最佳震源深度.反演结果显示,最佳双力偶解为,节面I:走向143°,倾角40°,滑动角71°,节面Ⅱ:走向347.2°,倾角52.6°,滑动角105.3°,地震矩震级为Mw5.9,最优深度解为7.7 km,与其他结果(CENC、IGP-CEA、Harvard)基本一致.利用滑动时窗相关法提取Pn、sPn震相,再利用其到时差测定震源深度为8.5 km,与CAP结果基本一致,验证了该方法的可行性.【总页数】6页(P624-628,638)【作者】尹欣欣;赵林林;杨立明;陈继锋;左可桢;蒲举【作者单位】中国地震局兰州地震研究所,兰州市东岗西路450号,730000;甘肃省地震局,兰州市东岗西路450号,730000;中国地震局兰州地震研究所,兰州市东岗西路450号,730000;中国地震局兰州地震研究所,兰州市东岗西路450号,730000;中国地震局兰州地震研究所,兰州市东岗西路450号,730000;中国地震局兰州地震研究所,兰州市东岗西路450号,730000;中国地震局兰州地震研究所,兰州市东岗西路450号,730000【正文语种】中文【中图分类】P315【相关文献】1.2016年1月21日青海门源Ms6.4余震序列重定位和主震震源机制解 [J], 梁姗姗;雷建设;徐志国;邹立晔;刘敬光2.2016年1月21日青海门源Ms6.4地震震源机制解及发震构造初步探讨 [J], 李晓峰;3.2016年1月21日青海门源MS6.4地震震源机制解及发震构造初步探讨 [J], 李晓峰4.2016年青海门源Ms6.4级地震多个震源机制中心解的确定 [J], 高熹微; 张骞; 蒋延林; 王福才; 赵卫红; 郭宇鑫; 李玉5.2016年青海门源M_S6.4地震震源机制与震源深度计算 [J], 李启雷;李玉丽;马丽;罗自浩因版权原因，仅展示原文概要，查看原文内容请购买。
震源机制解
震源机制解
震源机制解是地震学中的一个概念，它指的是研究地震发生的过程，以及地震的能量释放的机制。
在地壳内部，由于岩石层的构造变化和温度压力变化，岩石层会发生变形和断裂。
当岩石层发生断裂时，会产生巨大的能量，这就是地震的能量源。
地震能量释放的过程主要有两种机制：一种是滑动机制，即岩石层发生滑动，释放能量；另一种是压缩机制，即岩石层发生压缩，释放能量。
地震学家通过对震源机制解的研究，可以了解地震的发生原因，并预测地震的发生时间、地点和强度，为人类提供有效的防震减灾措施。
2023年11月—2024年2月中国大陆地区犕≥4．0地震震源机制解测定
梁姗姗，邹立晔，刘艳琼，等. 2023年11月—2024年2月中国大陆地区M ≥4.0地震震源机制解测定[J]. 地震科学进展,2024, 54(3): 229-236. doi:10.19987/j.dzkxjz.2024-036Liang S S, Zou L Y, Liu Y Q, et al. Determination of focal mechanism solutions of the earthquakes with M ≥4.0 occurred in the mainland of China during November 2023 to February 2024[J]. Progress in Earthquake Sciences, 2024, 54(3): 229-236. doi:10.19987/j.dzkxjz.2024-036地震科学数据应用2023年11月—2024年2月中国大陆地区M ≥4.0地震震源机制解测定梁姗姗※ 邹立晔 刘艳琼 任 枭（中国地震台网中心，北京 100045）摘要 本文利用中国地震台网记录的宽频带波形资料，采用近震全波形反演方法得到2023年11月1日—2024年2月29日发生在中国大陆地区的M ≥4.0共62次地震震源机制解。
结果显示逆断型45次，走滑型13次，正断型3次，未知型1次。
关键词 震源机制；震源参数；乌什M S 7.1地震中图分类号：P315.3+3 文献标识码： A 文章编号: 2096-7780(2024)03-0229-08doi ：10.19987/j.dzkxjz.2024-036Determination of focal mechanism solutions of the earthquakes with M ≥4.0occurred in the mainland of China during November 2023 to February 2024Liang Shanshan, Zou Liye, Liu Yanqiong, Ren Xiao (China Earthquake Networks Center, Beijing 100045, China)Abstract In this paper, the regional full waveform inversion using the broadband waveforms recorded by China Seismic Network were conducted, and the focal mechanism solutions of the 62 earthquakes with M ≥4.0 occurred in the mainland of China during November 2023 to February 2024 were obtained. The types of these focal mechanism solutions show 45 reverse faulting, 13 strike-slip faulting, 3 normal faulting and 1 odd earthquakes.Keywords focal mechanism; source parameters; Wushi M S 7.1 earthquake0 引言据中国地震台网测定，北京时间2023年11月1日0时—2024年2月29日24时，中国大陆地区共发生M ≥4.0地震71次（表1）。
2017年九寨沟70级地震序列震源机制解和构造应力场特征
摘 要: 九寨沟余震序列的震源机制和构造应力场有助于认识本次地震的发震构造和孕震机理。本文基于四 川区域地震台网的波形资料, 采用波形拟合(CAP)方法和 P 波初动+振幅比(HASH)方法反演得到 2017 年 8 月 8 日九寨沟 7.0 级地震序列中 59 次 ML≥3.0 地震的震源机制解, 并基于该结果采用阻尼线性逆推法(DRSSI), 计算研究区域的平均构造应力场, 给出该区域的应力场特征。结果显示, 利用 CAP 方法反演得到的本次主震 的最佳双力偶机制解节面 I: 走向 248°/倾角 86°/滑动角–169°, 节面 II: 走向 157°/倾角 79°/滑动角–4°, 矩震 级为 Mw6.31, 矩心深度 5 km, 属走滑型地震事件; 大部分余震的震源机制解错动类型与主震一致, 矩心深度 集中在 3~10 km; 应力场反演结果显示, 该区域周边的应力性质为走滑型, 最大主应力方向呈 NWW–SEE 向, 与该区域的应力场方向一致, 表明本次地震主要受区域应力的控制。结合该区域的地震地质构造等已有研究 成果, 分析认为此次地震的发震断层为走向 NW–SE、倾向 SW 的左旋走滑断裂——树正断裂, 巴颜喀拉块体 向 E-SE 向的水平运动受到华南块体的强烈阻挡导致此次地震的发生, 汶川地震的发生对本次地震具有一定 的促进作用。 关键词: 九寨沟地震; 震源机制解; 构造应力场; 发震构造 中图分类号: P315.1; P315.2 文献标志码: A doi: 10.3975/cagsb.2018.061901
本文由中国地震局监测预报司震情跟踪青年课题(编号: 2018010122)、四川省地震局地震科技专项(编号: LY1802; LY1817)和 2018 年 地震大形势跟踪任务联合资助。
震源机制解综述
震源机制解综述1、引言地震学是一门以观测资料为基础的研究地震的成因及其规律已成为地震预报的一种重要手段，它的发展奠定了地震预报的物理基础。
地震震源和地震波传播介质的各种参数在强震前的变化早就被当作地震预测的地震学前兆指标，随着地震预测的深入研究，以及我国“十五”台站数字化改造的完成，我们在进一步研究地震时空强分布特征的同时,加强对地震波的运动学和动力学特征的研究,从中提取震源，我们意识到加强对地震波的运动学和动力学的研究，从中提取震源信息，对增强地震预测的物理基础，提高地震预测的水平是十分必要的。
地震是地球内部物质运动的结果，这种运动反映在地壳上，使得地壳产生破裂，促成了断层的生成、发育和活动。
地震前后的地形变测量和地震波的观测研究等结果确认，天然构造地震是地下岩层的突然错动引起的。
发生错动的岩层可称为地震断层。
断层活动诱发了地震，地震发生又促成了断层的生成与发育，因此地震与断层有密切联系。
地壳中的断层密如织网。
实际地震断层的几何形状可能很复杂，但对多数地震，特别是小地震，作为初级近似，总体上可将地震看成是沿一个平面断层发生的突然错动引起的。
2、前人对震源机制解的研究历程地震震源处地球介质的运动方式。
通常所说的震源机制是狭义的，即专指研究构造地震的机制而言。
构造地震的机制是震源处介质的破裂和错动。
震源机制研究的内容包括，确定地震断层面的方位和岩体的错动方向，研究震源处岩体的破裂和运动特征，以及这些特征和震源所辐射的地震波之间的关系。
对地震震源的研究开始于20世纪初叶。
1910年提出的弹性回跳理论，首次明确表述了地震断层成因的概念。
在地震学的早期研究中，人们就已注意到P波到达时地面的初始振动有时是向上的,有时是向下的。
20世纪的10～20年代，许多地震学者在日本和欧洲的部分地区几乎同时发现，同一次地震在不同地点的台站记录,所得的P波初动方向具有四象限分布。
日本的中野广最早提出了震源的单力偶力系，第一次把断层的弹性回跳理论和P波初动的四象限分布联系起来。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
c b
c
b
平面 c 水平
a
节面 a、b 与平面 c 相 交，将平面 c 划分为四 部分，P波初动在水平 面内也呈对称分布。
工程地震导论
震源机制解
–震源机制的涵义 –断层破裂震源参数 –震源模型 –震源球面 – 赤平投影原理 – P波初动、震源节面和力轴的赤平投影
震源球面
O
地面
＋
－
F － 震源球
Sn
地心
M—莫氏面 C—康氏面
F — 震源；O — 震中；S1, S2, … , Sn — 观测点
震源球面与地震射线的关系
近震射线、远震射线与投影球面的关系
S2
S1
O
S3
S4
至…
远震射线
至 Sn-1
近震射线
F
至 S5
至 Sn
近震射线大多穿过上半球面，远震射线多数穿过下半球面
工程地震导论
震源机制解
– 震源应力场
力轴（压力轴 P、拉力轴 T、中性轴 N）产状 震前应力 0、震后应力 1 应力降 = (0 1)/2 平均应力 m = (0 +1)/2 视应力 s = m
– 地震矩 M0 = A Dm （ ：断裂岩体剪切模量）
工程地震导论
震源机制解
–震源机制的涵义 –断层破裂震源参数 –震源模型 –震源球面 – 赤平投影原理 – P波初动、震源节面和力轴的赤平投影
d
赤平面
180°
走向
S
— 平面的持平投影：大圆弧
l — 直线的持平投影：点
90° E
di r cos
倾向
di
d 倾向地震导论
震源机制解
–震源机制的涵义 –断层破裂震源参数 –震源模型 –震源球面 – 赤平投影原理 – P波初动、震源节面和力轴的赤平投影
P波初动、震源节面和力轴的赤平投影
–震源机制的涵义 –断层破裂震源参数 –震源模型 –震源球面 – 赤平投影原理 – P波初动、震源节面和力轴的赤平投影
赤平投影原理
U 上极点
N n 上半球
赤平面
W
赤道
下半球
S
投影球
下极点 D
下半球投影
E
直线的投影
l
平面的投影
直线和平面的赤平投影
赤平投影原理
投影结果
N
走向
0°
的法线
l W 270°
断层破裂震源参数
几何参数
– 节面（断面）产状 – 断面尺度（断裂长度L、宽度 W 或半径 r）
运动学参数
– 破裂传播方向 – 破裂传播速度 vf – 断层错动历时 – 错距 D（最大错距 Dmax和平均错距 Dm ）
断层破裂震源参数
力学参数
– 地震能量
震源释放能量 E、地震波能量 EW 地震效率 = EW / E
y
P轴
0°
T轴
270°
90°
x
T轴
180°
P轴
双力偶源P波辐射图案 Ur = Asin2
y
P轴
0°
T轴
90°
270°
x
T轴
180°
P轴
双力偶源S波辐射图案 U = Bcos2
倾斜节面（断层面）与水平面的交切
节面 a、b 倾斜正交
a
节面 a、b 将两者的法 平面分割为四个象限。 P波初动在法平面内呈 对称的四象限分布。
– In the case of a fault-related event it refers to the orientation of the fault plane that slipped and the slip vector and is also known as a fault-plane solution.
震源机制的涵义
震源机制
– 震源区在地震发生时的力学过程。 – 震源断层面的方位和滑动矢量，亦称“断层解”。
Focal mechanism
– the inelastic deformation in the source region that generates the seismic waves.
震源模型 — 集 中 力 系 模 型
单力源 z y x y
x
单力偶源 z y x y
x
双力偶源 z y
x y
x
震源模型 — 集 中 力 系 模 型
270°
y
0°
90°
x
180°
单力源P波辐射图案 Ur = Acos
y
0°
90°
270°
x
180°
单力源S波辐射图案 U = Bsin
辐射图案与各个方位上地震波波形的对应关系
y
0°
270°
180°
0°
30°
60°
90°
x
90°
120°
150°
180°
震源模型 — 集 中 力 系 模 型
y
P轴
0°
T轴
270°
90°
x
T轴
180°
P轴
单力偶源P波辐射图案 Ur = Asin2
y
P轴
0°
T轴
90°
270°
x
T轴
180°
P轴
单力偶源S波辐射图案 U = Bsin2
震源模型 — 集 中 力 系 模 型
＋
F — 震源；O — 震中
剖面上断层错动与P波初动方向的分布
震源球面
P
S2
震源球 PO P
F
地面
P — 直达波
Pg — 康氏面首波
Pn — 莫氏面首 波 Pg
Pn
S1
震源到每一个观测点都有一条射线相连。 每条射线与震源球面都有一个交点。 震源球可以任意小到震源球包围的所有 射线均可视为直线。
Pn
From Wikipedia, the free encyclopedia last modified on 6 November 2010 at 03:26.
工程地震导论
震源机制解
–震源机制的涵义 –断层破裂震源参数 –震源模型 –震源球面 – 赤平投影原理 – P波初动、震源节面和力轴的赤平投影
工程地震导论
中国地质大学（北京） 孙进忠
工程地震导论
绪论 地震和地震现象 地震活动性 震源机制解 场地地震效应 地震影响评价 场地工程地震勘察与评价
工程地震导论
震源机制解
工程地震导论
震源机制解
–震源机制的涵义 –断层破裂震源参数 –震源模型 –震源球面 – 赤平投影原理 – P波初动、震源节面和力轴的赤平投影
O
P
iP
近震直达波
iP
的离源角
震源球
F
Ah Av i P P
射
P
i P Av
线
Ah
P波初动位移图示
震源球
PO P
F
iPn iPn
地面
P iPg
g
iPg
首波的离源角
iPn
Pn
iPn
地心
M—莫氏面 C—康氏面
iPn
P波初动的赤平投影
将震源球面作为投影球，保留下半球投影。近震 射线反向延伸至下半球面。
将各个地震台站（观测点）与震源对应的射线按 直线投影到赤平面上。
台站投影参数
台站方位角： z
射线离源角： i （与射线的倾角互余） P 波初动方向： 向上 ，向下
P波初动、震源节面和力轴的赤平投影