单幅干涉图像提取震间形变场模拟研究

针对活动断层震间InSAR形变场提取中存在SAR观测数据不足、轨道误差难以合理去除的问题,该文进行了利用单幅InSAR图像提取震间形变场的模拟研究:考虑到InSAR图像中轨道误差和震间形变场的空间分布特征,利用一次多项式模型和深部滑动模型模拟了100组叠加的InSAR形变场;然后采用直接恢复震间形变场、先分离轨道误差再恢复震间形变场、一步法分离轨道误差与恢复震间形变场3种方法对模拟的震间形变场进行了恢复实验。结果表明,基于单幅InSAR图像,利用一步法能够较为准确地恢复震间形变场。     

2016年日本熊本地震破裂时空过程联合反演

为了深入认识2016年4月15日日本熊本地震破裂的复杂性,利用远场体波资料和同震InSAR资料联合反演了此次地震的震源破裂时空过程. 联合反演结果表明:熊本地震的震源破裂持续时间约为25 s,整个破裂过程释放的总标量矩为6.03×1019 N·m,对应于矩震级M_W7.1;同震滑动主要集中分布于浅部,破裂以右旋走滑为主,但在沿倾向0—5 km范围内,破裂呈较强的正断特征;此次地震破裂的最大同震滑动量约为4.9 m,且最大同震位错区位于背离断层走向上、距离起始破裂点约5—10 km的区域;破裂前期（0—7 s）,在倾向上向浅表发生破裂,在走向上向东北和西南两侧扩展;大约7 s后,破裂背离断层走向主要向东北方向扩展. 根据有限断层联合反演结果推测,此次熊本地震破裂可能出露至地表.      

门源Mw5.9级地震形变InSAR观测及区域断裂带深部几何形态

利用Sentinel-1A卫星升轨、降轨合成孔径雷达影像数据,提取了2016年门源Mw5.9级地震的高精度合成孔径雷达干涉同震形变场,利用单纯形法和非负最小二乘法反演确定了地震断层几何和滑动分布,并构建了区域断裂带的深部几何形态模型。结果表明,门源Mw5.9级地震同震形变以地表抬升为主,沿升轨、降轨视线向的最大值分别为5.3 cm、7.1 cm;地震断层走向、倾角分别为133°、43°;地震滑动以逆冲为主,主要发生在地下6.14~12.28 km处,最大滑动量约0.5 m,平均滑动角为66.85°,地震矩为1.0×1018 N·m（Mw5.94）;形变观测拟合残差均方根为0.36 cm;区域断裂带的深部几何形态以花状构造为特征,整体倾向南西,门源地震发震断裂为花状构造中未出露地表的盲断层。相关成果能够为研究区域地壳运动与变形、活动断裂与地震孕育发生等提供参考。     

2022年青海门源MW6.6地震发震构造——来自InSAR和高分影像约束

2022年1月8日,在青海省门源县海原断裂带冷龙岭断裂西段和托莱山断裂东端发生了M_W6.6地震,发育了明显地震地表破裂带.本文基于欧空局哨兵1号雷达影像,利用InSAR技术获取了门源地震的同震形变场,升降轨InSAR同震视线向位移表现出相反的形变特征,量级达到～60 cm左右,同震形变存在非对称性分布特征,结合高分7号观测解译了近断层地表形变.利用InSAR观测反演获得了发震断层参数及详细滑动分布,计算了同震库仑应力变化,并对发震构造及震中区域未来地震危险性进行了分析讨论.结果表明：门源地震至少有两条断裂发生了破裂,主断层对应冷龙岭断裂西段,InSAR确定的最优断层模型显示主断层东段存在沿走向变化特征,西段则在地质解译断层基础上向西延伸,次断层对应地质解译的托莱山断裂东端,两个断裂组成一个平躺的Y型分布.断层最大滑动量约为～3.7 m,断层浅部存在滑动,表明该地震破裂到了地表,地表破裂长度约19 km.主断层滑动主要集中在0～9 km深度范围;次断层滑动主要集中在0～4 km深度范围,InSAR确定的矩震级为M_W6.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以内。     

利用GPS和海底基准点观测形变反演日本大地震(Mw9.0)同震断层滑动分布

收集了国际上发布的通过日本区域密集的地表观测数据获得的高精度同震形变位移,同时收集了靠近震中的5个海底基准点观测到的同震变形.基于以上结果采用分断层滑动模型,基于分层弹性半空间位错理论,采用二阶拉普拉斯平滑算子方法约束断层面上的滑动量,利用最速下降法求解模型的最优解,以期精化该地震断层的滑动分布,进而研究其地震活动机制.研究结果表明,日本地震断层的平均滑动为6.95m,最大滑动量为58.70m,位于东经143.17°,北纬38.25°,地震所释放的能量为4.48×1022N·rn(等效震级为Mw9.07),地震断裂主要集中于东经142°～143.8°,北纬37.2°～39.5°,海平面以下10～30 km之间区域.     

日本2011 Tohoku-Oki Mw 9.0级地震的同震形变及其滑动分布反演:GPS和InSAR约束

利用覆盖整个日本东北部的3个条带的Envisat/ASAR降轨数据和6个条带的ALOS/PALSAR升轨数据,通过二通差分干涉处理,首先得到了2011年日本Tohoku-Oki Mw 9.0级地震的初步同震形变干涉结果,然后,利用GPS同震观测值对InSAR结果进行校正和基准统一。在此基础上,采用弹性半空间矩形位错模型,联合GPS和InSAR观测结果对同震滑动分布进行反演,获取了发震断层的断层滑动分布。结果显示,此次地震的滑动分布主要发生在40～50km深度范围内,最大滑移量为50.3m,释放的能量为3.20×1022 N·m(相当于Mw 8.94级)。     

基于GPS观测的红河断裂带现今分段及闭锁特征

红河断裂带是川滇地区重要的活动边界断裂之一,了解红河断裂带不同段落的活动特性对理解川滇地区以及青藏高原东南缘的运动特征有着重要意义.本文利用2014-2018年红河断裂带加密的GPS站点,结合1999-2016年中国大陆GPS速度场数据,获取了红河断裂带附近的现今GPS速度场.本文首先采用F检验法将研究区域划分为5个相...     

基于InSAR数据的西藏玛尼Ms7．9级地震的地壳不均匀性研究

基于InSAR观测结果,考虑到发震区域南、北两盘地壳介质的横向不均匀性,结合更为真实的地表断层位置,采用三维半解析粘弹性模型和遗传算法,反演了发震断层滑动分布及南北两盘的地壳介质参数(泊松比).反演结果显示,考虑介质横向不均匀性的模型能较好地拟合地表形变场;南北盘的横向介质差异性也可认为是玛尼地震同震形变的非对称性的影响因素之一.     

利用两种z变换算法的PS-DInSAR相位解缠与等价性证明

在介绍PS-DInSAR相位解缠函数模型的基础上,给出了应用LAMBDA方法求解模糊度和形变参数的过程,并将两种改进的z变换降相关算法——逆整乔列斯基和LLL应用于PS-DInSAR相位解缠。以z变换过程的迭代次数、z变换后的模糊度向量间的平均相关系数和协因数阵的谱条件数为准则,对两种算法进行仿真模拟和分析,结果表明逆整乔列斯基算法和LLL算法等价。最后从理论上对两种降相关算法的一致性进行了解释。