2008年于田MW7.2地震滑动分布反演及三维同震形变场解算

基于ENVISAT ASAR升降轨数据, 利用InSAR获取2008年于田M_W7.2地震同震形变场; 采用SDM反演本次地震断层滑动分布; 使用PSCMP正演获取于田M_W7.2地震南北向模拟形变量, 并结合升降轨同震形变场, 解算三维同震形变场。 同震形变场分析表明, 2008年于田M_W7.2地震以正断为主, 且带有走滑运动特征, 破裂带走向为NNE向。 同震滑动分布反演结果显示, 断层沿走向被分为4段F1、 F2、 F3、 F4, 其滑动分布集中在0~14 km区间, 以F2、 F3段为主, 最大滑动量约5.31 m, 位于F2段深部2.76 km处; 沿破裂带走向, 左旋走滑位移与垂直位移比值有增大的趋势; 反演获得的地震矩M₀=5.58×10¹⁹N·m, 相当于矩震级M_W7.1。 三维同震形变场解算结果显示, 断层上盘整体表现为沉降, 断层下盘整体表现为隆升, 且沉降量明显大于隆升量, 表明地震以正断破裂为主; 除靠近断裂带中上部表现为向东南运动外, 上盘整体上表现为向西南运动; 断层下盘则整体表现为向东北运动, 证明破裂兼有左旋走滑运动。 滑动分布反演、 正演与三维同震形变场解算结果皆表明, 于田M_W7.2地震破裂以正断为主, 且带有一定的左旋走滑。     

2015年尼泊尔地震三维发震构造及地震危险性研究

2015年4月25日尼泊尔M_W7.8地震发生在喜马拉雅造山带上,人们普遍认为该地震不足以释放该造山带上累积的能量,但对该地区后续地震危险性评价多基于二维或是假三维的形变反演计算结果.本研究从2015年尼泊尔地震主震与其最大余震M_W7.3地震之间的关系出发,着重分析讨论了尼泊尔地震的时间和空间的非均匀性,结合震源机制解、地壳速度结构、精定位后的余震分布及InSAR反演结果,建立了三维发震构造模型,利用非线性摩擦有限元方法,对一个地震周期内断层摩擦行为和块体变形进行了模拟,将计算结果和地表同震形变、形变反演的同震破裂、历史地震时空演化进行对比,在确认该三维模型可靠性的基础上,讨论了该地区后续地震的可能位置,认为在1934年比哈—尼泊尔M_W~8.1地震以东区域,还存在发生大地震的可能,在最大余震M_W7.3地震东南部位,还存在发生中大地震的可能.     

2011年日本MW9.1地震前后破裂区内视应力时空变化

2011年3月11日, 一个 M_W9.1地震袭击了日本本州地区, 为了分析这次地震前后主震破裂区内应力时空变化, 我们选取1996年1月~2016年6月期间发生在破裂区内的563个5.0≤M_S≤6.9地震, 研究了视应力随时间的变化和空间分布。 日本M_W9.1地震前从2002年中起视应力开始呈趋势性上升变化, 到2009年初以0.18 MPa/a的速率从0.6 MPa上升到1.76 MPa, 相差约3倍, 直到地震发生前夕一直保持在1.5 MPa之上。 地震发生之后, 直到2016年6月在破裂区内视应力呈缓慢下降变化, 但仍保持在1.5 MPa之上较高水平。 视应力在地面上和断层面上的分布显示, 1996—2005年间破裂区仅存在个别视应力高值, 从2006年到2011年2月, 破裂区大面积出现视应力高值。 在日本M_W9.1地震发生之后的近3个月内, 破裂区视应力整体处于高值水平, 之后在较高的水平上缓慢减弱。 视应力是地震断层面上平均应力的下限, 视应力的高低在一定程度上反映的是震源断层面上平均应力的高低。 在日本M_W9.1地震前, 发生在破裂区内的地震, 其断层面上的平均应力经历了大约8.5年的趋势上升变化过程。 这次大地震前破裂区所在的地壳应力逐渐增加, 最后导致断层面错动发生日本M_W9.1地震。     

利用InSAR技术研究2016年青海门源MW5.9地震同震形变场及断层滑动分布

2016年1月21日, 青海省门源县冷龙岭断裂带附近发生了M_W5.9地震。 基于Sentinel-1A影像, 采用差分干涉雷达测量技术研究了此次地震产生的同震形变场, 结果表明, 门源地震的形变影响范围约20~30 km, 形变态势在升降轨道形变场均显示为隆升, 基本沿冷龙岭断裂呈近似同心圆展布, 推测可能是冷龙岭断裂与民乐—大马营断裂之间的一条逆断层, 沿雷达视线方向最大形变量级约为6 cm。 均匀滑动反演显示门源发震断层长7.3 km, 宽6.2 km, 走向298.6°, 倾角34.5°, 倾向宽度9.5 km, 沿走向滑动量为170 mm, 沿倾向滑动量为460 mm, 矩震级为M_W5.97; 分布式滑动反演显示门源地震以逆冲为主, 兼具少量右旋走滑分量, 滑动量主要集中在沿断层倾向方向, 距离地表5～15 km处, 最大滑动量约0.3 m, 位于断层倾向深度10 km处, 矩震级为M_W5.93。     

2017年伊拉克MW7.3地震的类型界定及其震后趋势分析

2017年伊拉克地震发生在我们划分的巴格达地震区,鉴于不同机构提供的该震震级参数相差较大,本文利用孕震断层多锁固段脆性破裂理论,分情况讨论了该震所属地震类型,并分析了巴格达地震区地震趋势.结果表明：若2017年伊拉克地震为M_W7.3,则该震为第3锁固段向峰值强度点演化过程中发生的1次显著前震,该区未来将发生M_W7.7~8.2（双震型为M_W7.5~8.0）标志性地震,目前已接近临界状态;若2017年伊拉克地震为M_W7.5,除可能为显著前震外,还可能为标志性地震（双震）之一,若如此两年内该区将发生另一次M_W7.5地震;若2017年伊拉克地震为M_S7.8,则该震为第3锁固段峰值强度点对应的标志性地震,与我们对该震的前瞻性中长期预测结果相符.我们判断该震不为主震,预测该区未来还将发生M_W7.8~8.3（双震型为M_W7.6~8.1）标志性地震,目前该区远离临界状态.     

2017年8月8日四川九寨沟7.0级地震InSAR同震形变场及断层滑动分布反演

基于InSAR技术,利用欧空局升降轨Sentinel-1A/IW宽幅数据,获取了2017年8月8日四川九寨沟7.0级地震InSAR同震形变场,并以升降轨InSAR观测结果为约束,反演了断层滑动分布,基于三种不同接收断层计算了同震库仑应力变化.结果表明,同震形变场发生在塔藏断裂、岷江断裂和虎牙断裂交汇的三角地带,升降轨干涉位移均显示本次地震的形变场影响范围约为50 km×50 km,形变场长轴方向为NW向,升降轨观测的形变量相反,反映断层运动性质以走滑运动为主,升降轨数据观测得到的最大LOS （Line of Sight,视线向）形变量分别为~22 cm和~14 cm.非对称形变场反映出断层两侧的运动差异.反演结果显示,最大滑动量约为1 m,平均滑动角为-9°,矩震级为M_W6.5,地震破裂主要集中在地下1~15 km深度范围内,但整体而言本次地震破裂较为充分,基本将该区域1973年及1976年4次 > M_W6.0地震的破裂空区完全破裂.考虑到塔藏断裂和虎牙断裂的运动性质,可初步判定发震断层为虎牙断裂北侧延伸分支.基于三种不同接收断层模型的同震库仑应力变化计算结果反映出该区域以应力释放为主,进一步触发较大走滑型余震的可能性不大.     

2016年台湾美浓MW6.4地震震源参数的InSAR和GPS反演

2016年2月6日台湾西南部高雄市美浓区发生了M_W6.4地震.本文结合ALOS2卫星升降轨、Sentinel-1A升轨SAR数据,采用两轨差分干涉技术获取了该区域的同震形变场,形变结果表明震中西北部以抬升为主,最大视线向形变量约为11.2 cm.基于均匀位错模型和多峰值粒子群（MPSO）算法,利用InSAR和GPS形变数据联合反演了美浓地震的断层几何参数,结果表明震源中心位于22.920°N,120.420°E,深度约12 km,发震断层长度约15 km,走向角307°,倾角16.5°,平均滑动角为51.5°,此次地震是以逆冲倾滑兼左旋走滑的破裂模式.利用格网迭代搜索法得到最优倾角为15.7°,GPS和InSAR最优权比为18：1,最优平滑因子为0.06.基于非均匀位错模型,利用非负最小二乘方法进行线性反演,结果显示最大倾滑和走滑量分别为51.7 cm和55.3 cm,对应矩震级为M_W6.38,略小于GCMT （M_W6.4）的结果.通过与已有文献的比较和对该区域断层构造的分析,发现美浓地震的发震断层为单一断层的解释更为合理,我们推测发震断层是位于左镇、后甲里等断层之间的一条东南-西北走向往东北倾斜的盲断层,并初步推测2010年M_W6.3甲仙地震也同该断层有关.     

2008年10月西藏当雄MW6.3地震震后形变提取与余滑反演

2008年10月6日西藏当雄发生M_W6.3地震.本文利用震后2008年10月26日至2010年8月22日的16期ENVISAT ASAR数据,通过小基线集干涉测量、误差校正与MInTS（Multiscale InSAR Time Series）技术提取高精度的震后形变场,利用SDM（Steepest Descent Method）方法反演断层震后余滑演化过程,并分析震后余滑与同震滑动的关系.结果表明：当雄M_W6.3地震的近场震后形变场主要位于断层西侧,在时间演化上具有明显的对数函数衰减规律;震后余滑主要集中于断层中南段深0~15 km区间,最大的余滑量约0.07 m,位于断层深约9.28 km处,滑动角约-103°;震后余滑引起的地震矩能量M₀与矩震级M_W在时间演化上具有指数函数递增规律;当剪切模量μ=32 GPa,震后665天余滑释放的地震矩能量约为1.92×10¹⁷N·m,约占同震滑动释放地震矩的4.8%,相当于矩震级M_W5.46;虽然震后余滑已经延伸到断层浅部0~5 km区间,但由于余滑量相对较小,没有改变同震滑动在断层浅部区域的滑动亏损现象,这可能是2010年11月30日该区域又发生M_W5.3级余震的主要原因之一.     

2016年新疆阿克陶地震Sentinel-1 InSAR同震形变特征

2016年11月25日新疆克孜勒苏州阿克陶县发生M_W6.6地震。 本文利用合成孔径雷达差分干涉测量技术, 对Sentinel-1卫星获取的升、 降轨雷达数据进行了处理, 提取了该次地震的同震形变场, 并结合形变场特征与震源机制解, 采用梯度下降法反演发震断层的滑动分布。 结果表明, 升、 降轨LOS向同震形变场在发震断层两侧具有明显不同的形变特征, 主要形变区域分布在断层两侧, 升轨LOS向形变量可达-8.2 cm与11.2 cm, 降轨LOS向形变量可达-21.4 cm与13.1 cm; 反演的升、 降轨干涉形变场与InSAR测量值之间的残差得到有效控制, 大部分的残差介于±5 cm之间; 断层滑动分布主要集中于沿断层面深约2~18 km处, 最大滑动量位于沿断层面深约7 km处可达0.96 m; 平均滑动角约182.29°, 最大滑动处的滑动角约197.13°, 两个滑动分布中心的滑动角均接近180°, 表明阿克陶地震为一典型的右旋走滑破裂性事件; 当剪切模量取32 Gpa时, 反演的发震断层地震矩M₀可达9.75×10¹⁸, 相当于矩震级M_W6.60, 与地震波形反演结果一致。     

垂直摆倾斜仪地震效应与倾斜场模拟

在整理垂直摆地震响应的基础上,采用USGS提供的有限断层模型和速度结构模型,运用德国地学研究中心（GFZ）王荣江研究员开发的（PSGRN/PSCMP）计算程序,分析日本M_W 9.0地震和中俄交界M_S 6.6地震倾斜场变化,并与垂直摆记录数据进行对比。通过整理现有地质资料和断层数据,模拟依舒断裂带北段鹤岗周边发生M_S 6.0地震时,该地区同震及震后倾斜场的变化,以了解不同震源性质的地震产生的地倾斜变化及对鹤岗地区产生的影响。     

2011年日本MW9.0地震对沂沭断裂带及其两侧地区地壳运动的同震影响研究

2011年日本M_W9.0地震（简称日本地震）后沂沭断裂带及其两侧地区地震活动显著增强,研究日本地震对该地区地壳运动及地震潜势的影响十分必要.为此,本文通过112个连续GPS观测站获取了研究区高空间分辨率的日本地震同震形变场并得到如下认识：（1）8个定点地球物理观测的同震响应验证了本文同震形变场的可靠性;日本地震的东向拉张使研究区整体上处于张性同震应变状态,但存在局部挤压区域,其中莱州湾至海州湾的挤压条带穿过沂沭断裂带并对断裂带南北两段产生了不同的同震作用,对南段具有拉张作用,对北段产生挤压作用;（2）同震形变场在鲁东隆起和鲁西断块产生了显著的剪应变,地震b值显示上述区域的构造应力在日本地震后增强,因此同震形变场可能改变了这些区域的应力特征;（3）地震矩张量叠加分析显示,同震形变场短期内对鲁西断块、鲁东隆起区和沂沭断裂带南段累积了地震矩,可能有助于上述区域在日本地震以后的地震活动增强;日本地震对沂沭断裂带北段的地震矩具有释放作用,或许是该区域地震活动减弱的原因.     

天山北部强震活动对区域应力扰动的分析以2012年伊犁M6.6地震和2017年精河M6.6地震为例

天山地区为典型地震活跃区, 为定量分析该构造活跃区强震对周边构造变形和地震活动的影响, 本文基于地震位错理论和岩石圈分层模型计算了天山北部近期发生的2012年伊犁和2017年精河两次M6.6地震对周围地壳形变和应力的影响。 计算结果显示伊犁地震和精河M6.6地震引起震中附近地表同震位移达数厘米, 地表同震应变量级约为10^-7; 对比天山北部地区年平均构造形变特征, M6.6强震释放了震中附近近十年的构造主压应变积累; 地震引起震中附近(80 km内)同震库仑应力变化大于1 kPa, 而距震中较远区域活动断层上库仑应力变化微弱。 结合天山北部现今地壳变形特征及区域地震分布, 初步推测两次M6.6地震的发生对震后余震有显著的触发作用, 而对区域后续微震活动的影响微弱。     

利用高频GPS资料研究2016年新西兰凯库拉地震的地表形变及预警震级

2016年11月13日,新西兰凯库拉地区发生了M_W7.8级地震.本文利用1 Hz高频GPS观测数据,基于GAMIT track解算模块,采用主分量分析（PCA）空间滤波方法获取了地震地表形变.同时从trackRTr模拟实时解算的动态位移中,提取P波5秒峰值位移（P_d）和地面峰值位移（PGD）,并根据震级统计回归模型计算预警震级.结果表明：测站的动态形变时间长达2 min,距震源最近的HANM和KAIK站出现二次剧烈形变,震源北部测站的形变幅度大于南部,而高频GPS静态同震形变场表现出先逆冲后走滑的震源机制特征.不同GPS台站的P_d预警震级相差较大,最大震级差为M_W2.5.综合考虑预警震级发布的时效性和可靠性,采用顾及空间分布的四台站PGD联合预警方法,其预警震级在震后23 s达到初始稳定（M_W7.56）,在震后110 s达到最终稳定（M_W7.78）,该震级与USGS矩张量反演震级（M_W7.8）基本一致.     

应用近场速度记录快速估算川滇地区中强地震强度

采用Wu等(2004)提出的地震强度测定方法,应用近场速度波形记录,测定地震参数M_(ew),实现了川滇地区中强地震(M_W≥5.7)强度近实时快速测定。M_(ew)是一个与矩震级M_W对接的地震参数,与传统的应用振幅测定震级的方法相比,不会存在震级饱和问题,与全球矩心矩张量(GCMT)的矩震级M_W具有较好的一致性。该方法解决了应用近场波形资料和矩张量反演方法求取大地震矩震级比较困难的问题,可为区域中强地震速报提供近实时矩震级提供参考,从而提高地震速报的时效性和准确性,同时为强震加速度数据的应用奠定基础。     

地震大小的度量

本文首先对地震大小度量的研究历程作了简要回顾。在此基础上着重指出以下几点: ① 近几十年广泛使用的M_L、m_b(m_B)、M_S震级标度,作为对地震大小的度量,不仅都存在着“以偏概全”和“震级饱和”的问题,而且由于未充分顾及地震波衰减的区域差异,尤其台站场地对地震动放大作用的差异,因此都只是对地震大小的不精确的度量。 ② 地震矩M₀不仅物理含义明确,而且克服了M_L、m_b(m_B)、M_S震级标度所存在的各种问题,是最适合科学度量地震大小的物理量。为了继续延用“震级”这一术语,Kanamori定义了矩震级标度M_W,虽然其假定的前提条件仍有待研究,但也是对地震相对大小的较合理的度量。 ③ 为更加科学地度量地震的大小,必须充分利用现代区域数字地震台网产出的大量波形数据,加强地震波衰减特征、场地效应和震源参数测定及有关定标关系的研究,在完善M_L、m_b(m_B)、M_S震级测定方法的基础上,重点改进、完善M_W标度,逐渐推进以M_W作为统一度量地震相对大小的物理量,为地震科学研究和地震预测研究奠定更加坚实的基础。     

中小地震三维形变场重构方法研究与同震滑动分布反演——以2016年5月22日定日MW5.3地震为例

地震三维形变场对于研究地震发震机制等具有重要意义。 已有的InSAR三维形变场重构研究中, 只有中强地震的实例。 由于形变量级小、 InSAR方位向形变的误差较大, 中小地震的三维形变场重构易受噪声等影响。 本文以2016年5月22日西藏定日M_W5.3地震为例, 开展中小地震三维形变场重构的尝试。 首先基于InSAR技术获取了Sentinel-1升轨和降轨观测模式下的同震形变场, 再结合同震形变场的特点、 区域构造特征等, 添加限定方程(走滑运动为0), 重构了同震三维形变场。 结果显示, 震中附近以下降为主, 幅度达7 cm, 南北方向形变较小(约2 mm), 此区域还伴有2 cm的西向水平运动; 形变中心区域东西两侧部分区域均出现少许东向运动(1.5 cm)。 由同震形变场特征判断此次地震以正断破裂为主。 本文提出了基于连续性分层采样选取样本点方法, 以适应本地震形变场的实际情况。 对所得的LOS向位移场和重构的三维形变场进行降采样, 反演得到了断层面上的滑动分布, 两种数据得到的结果相似, 最优发震断层的走向约181°, 倾角约45°, 断层错动平均滑动角约-87.1°, 平均滑动量约为3.6 cm, 最大滑动量位于深度6.5 km处, 相当于一次M_W5.4的地震。     

2022年芦山MS 6.1、马尔康MS 6.0地震井水位、水温同震响应特征

基于全国地震地下流体台网数据,分析了芦山M_S 6.1、马尔康M_S 6.0地震引起的地下流体井水位、水温同震响应特征。结果表明,对于芦山M_S 6.1地震水位同震响应观测井数量较多,以上升变化为主,同震变化幅度较大;而对于马尔康M_S 6.0地震水位同震响应观测井数量较少,以振荡为主,变化幅度较小;2次地震引起的水位同震响应能力均强于水温测项,水温记录同震响应的前提是同井能记录到水位同震变化;2组地震引起的同震响应特征差异主要与震源机制解、台站分布密度、同震响应机理不同等有关。     

青藏高原东北缘断层系统的大地震迁移概率及断层滑动速度的分段特征

青藏高原东北缘是青藏高原隆升、生长及变形前缘.区域地震活动频繁,且地震在其主要断层带之间时空迁移.为了研究区域大地震在主要断层带之间的迁移规律与概率,以及主要断层带大地震破裂的时空分布特征,本文建立了青藏高原东北缘地区的三维黏弹塑性有限元模型,模拟了区域断层系统的地震循环,得到了人工合成的万年时间尺度的地震目录.根据模拟的地震目录,并结合古地震数据,计算分析了大地震(MW≥7)在研究区各个主要断层带之间的迁移概率,探讨了黏度、高程、统计时间长度等因素对大地震在各主要断层带之间的迁移概率和大地震在各主要断层带上的发生概率的影响,并且初步调查了海原断层带和香山天景山断层带的大地震破裂时空分布特征.研究结果显示:继区域最近两次大地震(1920年海原断层带上的M8.5海原大地震和1927年香山天景山断层带上的M8古浪大地震)之后,下一次大地震(MW≥7)发生在海原断层上的概率最大,约为51%～81%;其次是在香山天景山断层上,概率约为9%～37%.模型结果显示,不同的青藏高原中下地壳上地幔黏度大小,对大地震在各个断层带之间的迁移规律和迁移概率的影响较小;而研究区的高程载荷对地震迁移则有显著的影响:高程载荷易于使得海原断层地震活动减弱及香山天景山断层的地震活动增强.研究结果也显示了青藏高原东北缘地区主要断层带的地震活动与断层滑动速率分布的分段性显著;大地震在断层带上的破裂位置并不固定,呈现不均匀性;并暗示了断层几何形状对地震活动、断层滑动速率分布与大地震破裂位置的控制作用.     

2010智利莫尔M_W8.8地震对其破裂区及附近地区地震活动性的影响

强震的发生与地壳中应力增加有关。库仑破裂应力被用于研究由于某次强震发生所产生的应力变化。2010年智利莫尔M_W8.8地震发生之后,其破裂区以南地区没有发生M≥6.0地震,而在其以北地区于2015年在智利发生了伊亚佩尔M_W8.3级地震。在这两个地区地震活动性悬殊。本文利用库仑破裂应力研究了2010年智利莫尔M_W8.8地震在这两个地区产生的应力变化。结果表明,2010年智利莫尔M_W8.8地震在其破裂区以北地区产生的应力变化为正值,有利于断层错动发震,而在其破裂区以南地区产生的应力变化为负值,具有抑制断层错动发震的作用。