基于Landsat 8光学影像的巴基斯坦Awaran M_w7.7地震形变监测及参数反演研究

2013年9月24日巴基斯坦俾路支省(Balochistan)境内的阿瓦兰县(Awaran)发生了Mw7.7级地震.本文利用覆盖该地区的Landsat 8数据,基于影像配准的方法获取了该次地震的同震形变场,并运用地统计的方法对形变结果进行精度评定.针对传统四叉树算法中近场和远场中采样密度的不均匀性,以及噪音区域对数据降采样和反演结果收敛性的影响,本文提出了改进的四叉树算法对点的密度和形变梯度进行合理兼顾.最后利用光学影像获取的形变结果和数据的精度水平,基于Okada弹性半空间形变模型反演了该地震的震源参数和断层滑动分布.结果表明,地震断层北倾47°,滑动以左旋走滑为主,断层的西南部兼具少量的倾滑运动分量,断层滑动主要集中分布在断层面0~15km深度范围,最大滑动量达10m.反演获得的地震标量矩为4.68×1020 N·m,震级约为Mw7.75级.本文的研究结果可以为该地区的地壳应力变化研究和地震灾害评估提供依据,同时为Landsat 8光学影像应用于地震的形变研究提供参考.     

断层滑动分布与震源破裂过程联合反演研究

本论文的主要研究领域集中于地震震源破裂过程的联合反演。具体来说,就是试图通过对发生在中国大陆内部的3个典型震例的地震震源参数、地震断层几何结构及其地表形变特     

谱比法反演震源谱拐角频率的改进研究

由于地震波的信噪比,传播路径效应,震源机制差异和震源过程等因素的影响,传统的谱比法反演的拐角频率结果不精确,离散性大。以谱比法为基础,引入Bootstrap重采样法,KolmogorovSmirnov检验以及谱比趋势检验的筛选方法,对谱比法的反演结果进行精度评定和筛选,提高谱比法反演的拐角频率精度。将此方法应用于汶川地震余震的拐角频率、地震辐射能量、辐射效率以及能矩比等震源参数反演,建立拐角频率与地震矩之间的对应关系,验证谱比法以及三种筛选方法的可靠性。汶川地震余震震源参数反演结果表明,汶川地震余震的拐角频率与地震矩对应关系M₀∝f_c^-3.295±0.25存在与自相似关系M₀∝f_c^-3偏离现象。在地震释放能量方面,存在能矩比随地震矩增大而增大,以及大地震的辐射效率大于小地震的现象。     

基于组稀疏约束的微地震震源参数谱投影梯度反演

震源参数反演是微地震监测中的关键技术,常规走时或逆时定位方法可以快速获取震源的空间位置,但是会忽略震源的时间信息.全波形反演(FWI)是一种有效的工具,利用完整的波形信息,通过选用合适的优化算法对微地震事件震源参数进行迭代反演,虽然存在计算量大的问题,但是反演出的结果信息丰富并且精度较高.本文依据微地震震源的特点,提出...     

利用GPS和InSAR数据反演2011年日本东北MW9.0地震断层的同震滑动分布

本文首先对Envisat/ASAR数据进行干涉处理, 获取2011年日本东北M_W9.0地震的地表InSAR同震形变场; 然后通过对InSAR同震形变数据重采样方法的深入分析, 选择条纹率法结合干涉图的空间相干性对InSAR同震形变数据进行重采样; 最后基于弹性半空间位错模型, 联合InSAR与GPS形变数据, 采用最小二乘法反演发震断层的滑动分布. 研究结果表明: 考虑相干性的条纹率重采样方法, 更适用于形变场中存在除断层外的有限边界、 且形变场范围较大的InSAR数据重采样处理; 断层滑动主要发生在地表以下50 km范围内, 最大滑动量为49.9 m, 矩张量为4.89×1022 N·m, 所对应的矩震级为M_W9.1, 与地震学反演的结果比较吻合.      

Source parameters of the 2009 L'Aquila earthquake,Italy from Envisat and ALOS satellite SAR images

2009年4月6日意大利L’Aquila地区发生了Mw6.3级地震,该地震造成了300余人的人员死亡.本文联合不同波长、不同入射倾角的升降轨Envisat和ALOS卫星的差分干涉数据对该地震进行震源机制解的反演研究.研究首先对卫星雷达影像进行二通差分干涉处理,获取了覆盖L' Aquila地震震区的完整InSAR同震形变场,然后结合四叉树和均匀采样方法对原始观测数据进行降采样.在此基础上,联合GPS形变观测数据,利用弹性半空间矩形和三角位错模型,以及断层自动剖分技术对断层面进行最优离散剖分,反演获取了发震断层的精确几何参数和最优断层滑动分布,结果显示分布式三角位错滑动模型能够很好地解释观测到的地表形变场.反演结果表明发震断层是一个以正倾滑为主兼有少量右旋走滑的盲断层;基于观测数据最优确定的断层单元的最短边长为0.4 km,最长边长为6.3 km;此次地震的滑动分布主要发生在5～14 km深度的范围内,最大滑移量为1.07 m,释放的能量为3.43×1018N·m(Mw6.32),与地震学的研究结果非常一致.     

地震震源的研究:/

简述了1999~2002年中国学者在地震震源的观测、实验和理论方面所取得的一系列成果. 他们对精确定位、地震矩张量反演和震源参数反演等方法进行了改进和发展，并利用这些方法，研究了许多重要地震事件的破裂过程，并结合地面形变测量，研究了地震破裂与区域应力场、地质构造运动的关系. 总之, 在实验和理论研究领域中，中国学者在地震成因、地震发生的应力与构造条件、地震破裂动力学、岩石断裂力学、地震破裂成核等方面取得了一系列进展.      

联合Envisat和ALOS卫星影像确定 L'Aquila地震震源机制

2009年4月6日意大利L'Aquila地区发生了M_w6.3级地震,该地震造成了300余人的人员死亡. 本文联合不同波长、不同入射倾角的升降轨Envisat和ALOS卫星的差分干涉数据对该地震进行震源机制解的反演研究. 研究首先对卫星雷达影像进行二通差分干涉处理,获取了覆盖L'Aquila地震震区的完整InSAR同震形变场,然后结合四叉树和均匀采样方法对原始观测数据进行降采样. 在此基础上,联合GPS形变观测数据,利用弹性半空间矩形和三角位错模型,以及断层自动剖分技术对断层面进行最优离散剖分,反演获取了发震断层的精确几何参数和最优断层滑动分布,结果显示分布式三角位错滑动模型能够很好地解释观测到的地表形变场. 反演结果表明发震断层是一个以正倾滑为主兼有少量右旋走滑的盲断层;基于观测数据最优确定的断层单元的最短边长为0.4 km,最长边长为6.3 km;此次地震的滑动分布主要发生在5～14 km深度的范围内,最大滑移量为1.07 m,释放的能量为3.43×10¹⁸ N·m(M_w6.32),与地震学的研究结果非常一致.     

多孔径InSAR技术电离层校正方法及二维形变场应用研究——以玉树地震为例

差分合成孔径雷达干涉测量技术只能监测地表在雷达视线方向上的一维形变,这极大的限制了其发展和应用.本文详细介绍了一种利用单对干涉影像监测地表二维形变的方法——多孔径InSAR(MAI)技术,研究了基线误差和电离层对该技术的影响,发展了利用顾及方向性的滤波和插值方法改正电离层的影响.最后,以2010年玉树地震为研究对象,利用两景PALSAR影像给出了地表的同震二维形变场.研究结果表明,MAI技术要明显优于传统的灰度偏移量法,在同震形变监测和震源参数反演中具有巨大的潜力.     

2008年仲巴地震形变及其揭示的构造运动

2008年8月25日拉萨块体中部的仲巴县发生了MW6. 7地震,之后发生了MW6. 0余震。然而,针对此次地震,不同机构给出的震源位置和发震断层参数存在一定差异,同时对其所揭示的构造意义也需要更深刻的认识。文中利用不同波长、不同入射倾角的升降轨ENVISAT ASAR和ALOS PALSAR卫星数据,基于In SAR技术获取了该地震的8幅同震形变场图像,得到此次地震引起的地表形变长约50km,且分为南、北2瓣。为了更好地反演主震的破裂滑动分布,在数据降采样过程中,设置近场区域的采样间隔远小于远场区域,以保证主要形变区域的观测数据特征和采样密度,并且扣除了余震形变的影响。基于Okada均匀弹性半空间位错模型和SDM方法反演的结果表明,此次地震的发震断层为E倾的帕龙错断裂,破裂滑动主要为正断倾滑兼部分左旋走滑性质,南段最大滑动量为1. 15m,北段最大滑动量为0. 52m。为了进一步理解此次地震的构造背景,文中利用SBAS-In SAR技术获得了震后仲巴地区2008—2010年的形变时间序列,发现形变主要是由剥蚀导致的沉降区堆积所引起的,而在断裂两侧并不明显。结合1991—2015年的GPS数据分析,认为拉萨块体内部EW向的伸展作用是明显且不均匀的,地堑附近伸展量逐渐减小。正断型地震在高原隆升过程中可能起着重要的调节作用。