利用STLN和InSAR数据反演2008年青海大柴旦Mw6.3级地震断层参数

利用Envisat/ASAR雷达影像处理后获取的2008年11月10日青海大柴旦Mw6.3级地震InSAR同震形变场数据和LN算法反演观测值含有粗差情况下的地震Okada矩形位错模型断层参数;断层深度为19.691 5 km,倾角为56.892 3°,长度为16.499 3km,宽度为6.854 8km,倾滑量为0.832 6 m,中心经度为95.884 7°,中心纬度为37.529 5°,走向方位角为116.411 5°,地震矩为3.108 4×1018 Nm(Mw6.295 0)。反演结果表明,当观测值含有粗差时,LN算法在1范数下可以有效抵抗粗差的干扰,充分利用已有先验信息的同时顾及观测、线性化等原因引起的系数矩阵误差的影响。     

基于GIS的InSAR结果分析方法及在汶川Mw7.9级地震同震解释中的应用

利用GIS对InSAR结果进行分析可以更全面地解释地表形变机制。介绍了利用GIS对SAR影像进行裁剪、切割和拼接处理、剖面和等值分析的方法。在此基础上,以汶川Mw7.9级地震为例介绍了基于GIS的InSAR结果分析方法在地震形变机制解释方面的应用:①根据InSAR距离向偏移量和同震形变场,利用GIS技术提取了汶川Mw7.9级地震发震断层的地表迹线;②将基于位错模型正演得到的视线向形变场与InSAR同震形变场进行比较,得到了汶川地区InSAR同震形变场的上盘和整体改正数,进一步借助GIS进行修正得到了更为真实的InSAR同震形变场;③利用2D/3DInSAR同震形变场、干涉纹图和剖面分析结果对汶川Mw7.9级地震同震形变特征进行了具体分析和解释。     

联合地震位错模型和ESISTEM方法提取地震同震三维形变场

高精度同震三维形变场对于研究地震变形模式、震源机制等具有重要意义。设计并实现了一种联合地震位错模型和扩展融合大地测量、卫星形变观测、应变张量估计(extended simultaneous and integrated strain tensor estimation from geodetic and satellite deformation measurements, ESISTEM)方法的新方法,以2021年Ms 6.4漾濞地震为例,利用哨兵1号A、B星(Sentinel-1A/B)升、降轨影像获得该地震的合成孔径雷达干涉测量形变场,利用地震位错模型正演得到的南北向形变分量进行约束,成功提取了该地震完整的同震三维形变场及应变场。结果表明,漾濞地震断层西南侧主要向西、向北运动,最大形变分别为4.8 cm、9.5 cm;东北侧主要向东、向南运动,最大形变分别为7.4 cm、4.6 cm;垂直向抬升、沉降的最大值分别为3.6 cm、3.4 cm;发震断层以右旋走滑运动为主,兼有少量正断分量;发震断层区域受到显著的膨胀、剪切和旋转作用。     

活动地块运动和应变模型辨识

推导了块体运动和应变的各种模型，讨论了各种模型的适用条件。利用模拟数据反演计算了块体运动和应变参数。基于模型的无偏性、有效性等概念，对块体的刚性旋转运动模型、旋转应变模型进行了辨识，结果表明：同时考虑块体刚性旋转运动和应变的模型比单独考虑刚性旋转的模型更符合实际，模型的适用范围更广。结合中国地壳运动观测网络及华北地区网实际观测所得的GPS速度场，对华北地区块体运动和应变模型进行了辨识。     

同震变形分析中不同地球结构模型差异性研究

顾及地球结构的非均匀性、分层和球体特性，定量研究了球体分层模型、球体均匀模型、半空间均匀模型分别作用于地震断层（走滑断层，逆冲断层）产生的地表同震变形的差异。计算结果显示，相对实际地球模型，利用半空间均匀位错模型计算的同震变形差异很大，大大超出现有高精度测量技术的精度。如果利用实测资料（例如GPS，VIBL，InSAR等高精度测量数据）进行反演研究，必须根据实际地壳构造情况考虑地壳分层结构的影响。     

用多时相InSAR技术研究廊坊地区地下水体积变化

利用多时相InSAR(multi-temporal InSAR,MTI)技术,以河北廊坊市城区为试验区,利用两个Track覆盖的SAR数据获取了该地区2007～2010年间的高精度地面沉降速率场。结果表明,廊坊市中心城区存在一个非常明显的沉降漏斗,其最大地面沉降速率为76.0mm/a。在此基础上,利用自由应力下的体积应变模型进行地下水体积变化反演,显示该地区的体积变化最大达到-1.15×10-4/a,整个研究区的地下水年抽取体积达554万m3。     

利用PALSAR数据研究长白山火山活动性

利用2006-10~2011-03期间共24景ALOS卫星的PALSAR数据对长白山天池火山区的近年形变进行了研究,基于PSInSAR方法进行了数据处理,获取了长白山近年火山的时序形变场和形变速率结果。结果表明,长白山火山近年火山活动较弱,仅在火山口附近存在10mm/a的形变速率。同时利用Mogi模型对火山区的岩浆变化体积速率进行了反演,显示该区域的岩浆以~1×106 m3/a的速率补充。     

升降轨InSAR数据约束下的2007年阿里地震反演分析

2007年5月5日,西藏阿里地区发生了Mw6.1级地震。本文采用Envisat卫星的升、降轨SAR数据获取地震的精确同震地表位移,然后基于弹性半空间位错模型,分别采用MPSO非线性和最小二乘线性反演算法确定了断层几何参数和滑动分布模型。结果表明,分布式滑动分布模型能较好地解释观测到的InSAR地表形变场。地震发震断层为走向158°、倾角43°的西南倾断层,主要的滑动量集中在7~12km的深度,最大滑动量约0.3m,位于9km的深度。反演给出的地震矩为1.24×1018 Nm,与地震学结果相一致。     

利用双台站背景噪声分析2017年墨西哥MW7.1地震震源区的地震波速变化

背景噪声中包含有丰富的地球内部结构信息,可以用来监测地震波速度变化.基于墨西哥中部G.UNM和MX.TLIG两个地震台站2016年4月1日至2018年4月1日的连续波形资料,本文首先对其进行不同频带下的背景噪声互相关处理,然后采用小波变换对获取到的互相关函数进行滤波来解决台站间距过长所引起的互相关函数不稳定问题,最后分别通过Stretching方法和MWCS方法测量得到的走时偏移来提取地震波速度变化.结果表明：利用小波变换能有效提高互相关信噪比和结果精度;采用Stretching方法发现震源区浅部5~10 km的地震波速在地震中上升约0.07%,而在2~5 km的地震波速则有约0.05%的下降趋势;MWCS方法计算结果稳定性稍差,但变化趋势与Stretching结果相同.造成该地震波速变化的主要原因可能是2017年地震所引起的脆性岩层结构破裂.     

两次大柴旦Mw6.3地震间地表形变的InSAR观测及与同震破裂的联合分析

基于2008-11-10、2009-08-28两次大柴旦Mw6.3地震期间的6景Envisat SAR影像,利用小基线集技术提取相应的InSAR地表形变场。结果显示,震中周缘区域未观测到显著的震后形变空间模式,6个特征点区域的形变序列也未显示由快到慢的时变特征。基于震后形变理论对获取的地表形变进行余滑解释,分析两次地震同震破裂的主要参数,建立震时及期间地下断层滑移过程。两次地震分别在不同深度发生破裂且平均滑移量差异较大,前次地震未产生显著的浅层震后余滑。地震区域岩石圈8.2~23.7km深度内在同震阶段的滑移亏损可能以震间无震滑动的形式予以补偿。