基于曲面断层结构的2007年所罗门群岛地震同震滑动分布反演

2007年4月1日,南太平洋岛国所罗门群岛的俯冲带地区发生了Mw 8.1级地震.利用日本宇宙航空研究开发机构(Japan Aerospace Exploration Agency,JAXA)的ALOS/PALSAR卫星雷达影像数据,采用二通差分干涉技术获得了覆盖所罗门群岛地震震区的同震地表形变场,根据四叉树方法对原始观...     

基于ALOS PALSAR数据进行雷达差分干涉提取玉树地震同震地表形变场

本文针对2010年4月14日玉树地震引起的地表形变,使用日本ALOS卫星PALSAR L波段雷达影像数据,应用两轨雷达差分干涉(DInSAR)处理得到了以玉树为中心11000km2范围内的同震形变场,空间分辨率为8m,并在此基础上对玉树地震的震源机制和发震机理进行了分析.研究结果表明L波段雷达数据适合在地形起伏较大的地区进行DInSAR形变探测.该同震形变场信息可为玉树地震的同震形变反演提供参考数据.该研究进一步证实DInSAR技术在大规模地表形变探测和地学研究领域具有广阔的应用前景.     

2022年泸定Mw 6.6地震InSAR同震形变与滑动分布

2022-09-05,青藏高原东缘的鲜水河断裂上发生了泸定Mw 6.6地震,该地震是鲜水河断裂上40年来发生的最大地震,研究该地震的运动学和同震破裂模式对理解青藏高原东缘构造形变机制和评估鲜水河断裂以及安宁河断裂的地震危险性具有重要意义。利用Sentinel-1和ALOS-2卫星雷达影像,采用合成孔径雷达干涉技术获取了泸定地震的同震形变场,进而基于弹性半空间的位错模型,确定了本次地震发震断层的几何参数和滑动分布。结果表明,泸定地震是一次典型的左旋走滑事件,发震断层西倾,倾角约为72°,走向沿NNW-SSE方向,约为167°;断层破裂主要集中在0~10 km深度,最大滑动发生在约5.8 km深度,约为2.23 m;同震释放的地震矩约为8.74×1018 N·m,相当于矩震级Mw 6.59。通过对震后光学影像解译,发现此次地震诱发的滑坡多集中分布在发震断层西侧,该现象与余震主要集中在断层西侧的结果相一致,可认为是地震上盘效应的体现。     

大地测量地震断层同震滑动分布反演的两步解法

针对同震滑动分布反演中系数矩阵出现病态的问题,提出两步解法,并在两步解法反演过程中引入拉普拉斯二阶平滑矩阵进行平滑约束。该方法不仅改善了系数矩阵的病态问题,同时也很好地抑制了相邻断层面间出现大的梯度变化。在两步解法反演过程中,用L曲线法确定正则化参数。系统模拟实验表明,对于最大滑动量,该方法的反演结果较一步最小二乘法的反演结果精度提高了3.34%~19%;对于均方根误差,该方法的反演结果较一步最小二乘法减小了3.3%~13.3%。芦山地震反演结果表明,利用两步解法进行滑动分布反演是可行的。     

利用方差分量估计的地震同震滑动分布反演

提出了基于方差分量估计的地震同震滑动分布反演方法,该方法不仅可以确定不同数据集的相对权重,而且还能得到平滑因子的值,突破了方差分量估计方法仅用于各类数据定权的实际应用范围。模拟算例验证了此方法的有效性,从解的概率后验密度分布分析来看,方差分量估计方法得到的结果是线性和非线性混合方法反演结果的一个解。用方差分量估计方法反演了Bam地震同震滑动分布,并与线性非线性混合方法反演的结果进行了对比分析。两种方法得到的Bam地震断层滑动分布之间的空间互相关系数为0.999 9,地震同震的最大滑动量相同,都为3.04 m;其深度仅有微小区别,前者为4.50 km,后者为4.47 km;平均滑动量前者为0.714 m,后者为0.718 m。但前者反演计算用时227 s,后者却达到1.5×10⁶ s,约17 d,前者仅为后者的1/6 608。结果表明,本文方法相对于线性和非线性混合方法具有计算简单、计算量小和计算效率高的优点。     

九寨沟Ms7.0地震的同震形变场误差研究及断层滑动分布反演

北京时间2017年8月8日四川省九寨沟发生Ms7.0级地震。为了提高断层滑动分布反演的可靠性,本文利用升降轨InSAR数据联合反演断层滑动分布,再根据模拟形变值二次反演,对比反演结果来探讨升降轨形变场误差和确定断层滑动分布的结果。结果表明,九寨沟升轨同震形变场的质量较好,降轨形变场受余震、震后粘弹性松弛效应影响明显,利用模拟形变值二次反演确定的断层滑动分布可靠性更高。     

基于多源SAR数据的2022年门源Ms 6.9地震同震破裂模型反演研究

2022年1月8日青海省海北州门源县发生Ms 6.9地震,震中位于青藏高原东北缘祁连-海原断裂中段,属历史地震空区,基于多源合成孔径雷达（synthetic aperture radar,SAR）遥感数据研究该地震的破裂模式对理解青藏高原东北缘构造变形机制、应变释放过程以及地震危险性评估具有重要意义。首先利用Sentinel-1数据和合成孔径雷达差分干涉测量（differential interferometry synthetic aperture radar,D-InSAR）技术获取了门源地震的同震形变场,视线（line of sight,LOS）向形变场显示此次地震造成了约20 km长的地表破裂,最大形变约0.75 m;然后基于Sentinel-2卫星数据,利用光学影像配准和相关技术获取了本次地震的东西向同震形变场,最大同震位移达2.5 m;最后基于均匀弹性半无限位错模型,以LOS向形变场为约束反演了断层的滑动分布模型。结果显示,门源地震是一次典型的左旋走滑型地震,地震破裂主要集中在0~10 km深度范围,最大滑动量3.25 m,滑动角10.44°,对应深度4.89 km;反演给出的矩震量为1.07×1019 N·m,对应矩震级Mw 6.6。结合野外考察和地质资料,初步判定发震断裂为冷龙岭断裂,并引起托莱山断裂发生同震滑动。同震库仑应力结果显示,冷龙岭断裂东段和托莱山断裂西段应力状态为加载,未来具有发生强震的风险。     

地震同震滑动分布反演的总体最小二乘方法

同震滑动分布参数与地表形变间的线性关系依赖于格林函数矩阵的构造,格林函数矩阵元素与破裂面位置、几何参数、破裂方式及位错模型假设等因素有关。本文尝试考虑格林函数矩阵元素的误差来补偿上述原因在一定程度上对反演参数的影响,采用同时顾及系数矩阵(格林函数矩阵)和观测向量两者误差的总体最小二乘方法反演同震滑动分布。首先确定了系数矩阵元素和观测向量的协因数矩阵,考虑到格林函数矩阵的病态性(秩亏),借助拉普拉斯二阶平滑得到正则化矩阵,采用总体最小二乘正则化法反演同震滑动分布。并对2009年意大利中部拉奎拉(L’Aquila)Mw6.3级地震实例进行同震滑动分布反演研究。结果表明,拉奎拉地震的走向为144.37°,倾角为59.06°,滑动分布的最大滑动量为0.95m,平均滑动角为-96.4°,主要滑动深度为4~15km的范围,地震矩为3.63×10~(18)N·m,对应的矩震级为Mw6.34。总体最小二乘与最小二乘法的滑动分布解存在一定差别,但差别的量级在10-4以内。     

地震同震滑动分布反演平滑因子的确定

针对大地测量地震同震滑动分布反演中平滑因子（又称正则化参数）的确定问题,一般采用模型粗糙度与数据拟合残差之间的折中曲线来确定（为便于区分,本文将该方法称为"L曲线"）。本文在L曲线的基础上提出一种确定平滑因子的新方法——折中相交曲线法。模拟试验反演结果表明,利用折中相交曲线法确定的平滑因子反演地震滑动分布各参数精度要优于L曲线法。利用折中相交曲线法确定平滑因子反演拉奎拉与台湾美浓实际地震,并与L曲线法反演结果进行对比分析。分析结果表明,利用折中相交曲线法确定的平滑因子反演拉奎拉与台湾美浓实际同震滑动分布各参数结果均在国内外其他学者研究的范围内,并且利用折中相交曲线法确定平滑因子较L曲线法具有计算效率高、无需依赖数据拟合度、确定平滑因子大小更为合适等优点。     

青海玉树地震差分干涉雷达同震形变测量

2010-04-14青海玉树发生7.1级地震后,作者利用震前和震后获取的日本ALOS卫星PALSAR遥感数据,开展了差分干涉雷达(D-InSAR)地震同震形变测量与分析。结果表明:玉树地震引起较大范围地表变形,地震变形沿玉树—甘孜断裂带向南东东方向扩展,在N33.7°,E96.81°附近达到最大形变量,D-InSAR监测到雷达视向上的最大形变量为35cm。地表形变特征对于评价玉树地震破坏程度、推断断层性质、研究地震形变和地震孕育特征具有重要的参考价值。     

2015年尼泊尔地震同震滑动及震后余滑的三角位错模型反演

2015年尼泊尔Mw 7.9级地震发生在印度板块向欧亚板块低角度俯冲的喜马拉雅断裂带上。对该地震的滑动模型进行精化,对于理解喜马拉雅断裂带的孕震模式具有重要意义。采用三角形位错元构建主喜马拉雅断裂“双断坡”几何模型,联合全球定位系统（Global Positioning System,GPS）和合成孔径雷达干涉测量（interferometric synthetic aperture radar,InSAR）资料反演2015年尼泊尔地震同震滑移及震后余滑。结果表明,尼泊尔地震最大同震滑移达到7.8 m,深度为15 km,位于中地壳断坡和浅层断坪的接触部位。不考虑中地壳断坡结构会使反演的最大滑移量偏低。震后余滑主要分布在同震破裂区北侧,释放的地震矩为1.02×1020 N·m,相当于一次Mw 7.3级地震,约占主震释放地震矩的12%。同震库伦应力变化和震间断层闭锁分布均表明,尼泊尔地震破裂区南部宽约60 km的区域仍具有较高的地震危险性。     

利用GPS资料反演汶川Mw7.9级地震滑动分布

2008年5月12日,青藏高原东缘龙门山断裂带发生Mw 7.9级汶川地震,该地震使得北川-映秀断裂、灌县-江油断裂和汉旺-白鹿断裂带发生了同震破裂。本文主要利用“中国地壳运动观测网络”项目组提供的GPS同震位移资料,采用基于敏感度迭代拟合法(SBIF)及分段发震断层模型,分别应用均匀半空间和分层地壳结构两种方法反演分析了该次地震的精细同震滑动分布。结果表明：（1）SBIF反演方法应用于相对较小的GPS数据集仍然可以获得理想的同震滑动分布;（2）采用考虑分层地壳结构模型的反演结果整体上要优于均匀半空间模型的结果;（3）同震滑动主要发生在10-19千米深度以上,且存在五个滑动峰值区域及北川断裂南西端存在明显的深部滑动区,与震后野外考察结果及余震分布相吻合,同时能够较好地解释GPS观测数据;（4）青川断裂存在明显的右旋走滑分量,其平均滑动量为1.99米,北川断裂主要发生逆冲滑动,平均滑动量为3.35米,灌县断裂的平均滑动量为0.65米;（5）反演得到的地震矩为8.74×1020Nm(Mw 7.90),与地震学结果吻合。     

利用GPS观测反演Tokachi-Oki地震(日本)及强余震的同震滑移分布

根据GPS观测的Tokachi-Oki地震及强余震引起的地表位移,采用基于敏感性探测的方法反演了两次地震分别引起的断层滑移。反演过程中,考虑了更切合实际的分层地球模型和断层面几何形状。主震和强余震滑移量的最大值分别为3.07m和0.33m,相应的矩震级分别为Mw7.94和Mw7.1,与地震波形反演的结果基本一致。另外,从主震和强余震的滑移分布来看,强余震的滑移主要分布在主震破裂较小的区域,一种合理的解释是,主震的破裂过程增加了邻区断层面上的能量积累,导致了强余震的发生。     

利用PALSAR数据获取拉奎拉地震同震形变场

2009年4月6日发生的意大利拉奎拉Mw6.3地震造成了地表的严重破坏,为了解地震引起的地表变化情况,本文利用ALOS/PALSAR影像数据,采用二轨差分干涉测量技术提取了该次地震的地表同震形变场。通过分析可知:地表在西北-东南方向上发生错动,形变主要发生在40×30 km的范围内,西南部出现了沿视线方向的下降,东北部出现了沿视线方向的抬升,最大值分别为0.337 m和0.122 m。结果表明,获得的同震形变场与地震地质调查的结果一致。     

2013年芦山Ms 7.0级地震断层参数模型反演

2013年4月20日四川芦山地区发生了Ms 7.0级地震。利用GPS三维同震形变数据获取地表形变场,基于位错模型反演芦山地震的断层几何参数及滑动分布。首先采用多峰值颗粒群算法（multiple peak particle swarm optimization,MPSO）得到断层几何参数,其中断层走向206.47°,倾角44.11°,长度21.94 km,离地表最浅处为7.66 km,最深处为17.84 km。为了反演断层面的精细滑动分布,分析地震所在的龙门山断裂西南段破裂面具有的铲状型特征,将芦山地震破裂面确立为铲状模型,即将断层的倾角预设为上陡下缓,倾角变化范围为21°~50°。结果显示,断层破裂面在不同深度区域出现了两个滑动峰值,其中最大滑动量为0.68 m,深度位于13 km。地震释放的能量为1.47×1019 N·m,对应的矩震级为Mw 6.74,与地震学的研究结果一致。