顾及地壳粘弹性结构的地震断层震后形变反演分析

顾及地壳的粘弹性结构特征,建立了基于粘弹性地球模型的地震断层震后反演模型。利用遗传算法重点对走滑断层和逆断层两种主要地震断层进行了震后形变反演研究,并与球体均匀模型反演进行了比较。研究结果表明:遗传算法能在较大范围内搜索出发震断层的位错参数,反演的结果基本可靠;基于粘弹性地球模型特征建立的地震断层震后反演模型,可较好地反演不同类型的地震断层震源位错参数;如果出现模型参数不清楚的情况,需综合考虑衡量反演结果的准则,找寻其他方法来判断,仅以VTPV最小作为反演准则并不可靠。     

基于GNSS/InSAR/GIS的活动断层地震危险性评估系统

本文从GNSS、InSAR、重力、水准等多种大地测量资料联合确定活动断层地壳形变、应变和应力,以及与GIS相结合的角度出发,讨论了基于库仑破裂准则的活动断层地震危险性评估方法,系统阐述了活动断层地震危险性评估系统的设计。该系统主要利用各种大地测量数据联合确定活动断层区域的位移场,反演计算应变场、应力场、活动断层运动参数以及库仑应力;通过确定活动断层的临界应力值,建立地震危险性评估模型,评估可能发生地震的位置、最大震级、复发时间及其对应的概率;最后借助GIS实现断层区域的地质和构造环境、三维形变数据、应变场、应力场和地震危险性综合评价图的管理与可视化。本文最后讨论了所建系统的具体应用,讨论了汶川地震激发的同震静态库仑应力、应力扰动的持续时间和龙门山地区的地震活动率。     

利用空间观测技术研究青海玛多7.4级地震孕育发生变形时空特征

为揭示2021年青海玛多7.4级地震震前变形与同震破裂特征,本文系统收集与处理了青海省CORS站和中国地壳运动观测网络等GNSS资料、以及Sentinel-1和ALOS-2卫星SAR数据,通过多源数据综合处理获取了玛多地震形变场动态结果.震前GNSS形变场结果显示,玛多震源区处于其北侧剪切形变高值区边缘,震前存在左旋剪...     

InSAR技术及应用研究进展

合成孔径雷达干涉测量(InSAR)是一种快速发展的大地测量技术,能够全天候获取高精度、连续覆盖的地面高程和地表形变信息.InSAR已在地形测绘、全球环境变化、灾害监测评估等相关领域得到了广泛应用并取得了一系列成果.主要综述了近20 a来InSAR技术的发展,以及该技术在数字高程模型提取和地球科学应用中的研究进展,并对该技术的发展趋势及应用研究进行了展望.     

综合InSAR和应变张量估计2016年MW7.0熊本地震同震三维形变场

利用合成孔径雷达(Synthetic Aperture Radar,SAR)影像提取与地质活动相关的三维地表形变场,对深入理解地质灾害的形成机制及其潜在灾害风险评估非常重要.目前,利用SAR影像的同震三维形变场提取主要利用单个像素点的多次观测构建观测方程,然后基于加权最小二乘(Weighted Least Squares,WLS)方法分解从而获得同震三维形变场,因此该方法缺乏对相邻像素点空间相关性的约束.考虑相邻同震位移点的应力连续性,研究学者提出了顾及大地测量应变张量和卫星形变观测的SAR同震三维形变场方法(Extended Simultaneous and Integrated Strain Tensor Estimation from geodetic and satellite deformation Measurements,ESISTEM).本文以2016年MW7.0熊本地震为例,收集了覆盖此次地震的ALOS-2卫星升降轨影像,利用传统差分InSAR(DInSAR)方法和子孔径雷达干涉测量(Multiple Aperture InSAR,MAI)方法分别对升降轨SAR影像对进行处理,得到视线向(LOS)形变和方位向形变,最后利用ESISTEM方法获取此次地震的三维同震形变场.此外,利用GPS和野外考察观测对本文的三维形变场结果进行结果精度分析.研究结果表明,与传统WLS方法相比,ESISTEM方法不仅能有效抑制奇异像素点对形变结果的干扰,同时对近断层的失相干信号能进行较好的恢复,更有助于解释地表破裂区的地震形变特征和掌握地震发生机制.本文确定的三维同震形变场结果显示主形变区发生在Futagawa断层中部和Hinagu断层最北端,最大水平位移为2m,抬升为0.55m.断层破裂以NE-SW走向的右旋走滑为主兼有部分正断成分.应变张量分析表明发震断层处受到了明显的收缩力和剪切力的作用.     

川滇地区GPS速度场聚类分析

以川滇地区GPS水平速度场为研究对象,首先采用层次聚类法对其进行直接聚类,然后对在地理空间中出现错分的GPS站点利用F统计检验来纠正分类。在此基础上,考虑DB指标以及层次聚类的不可逆性,对分类结果进行二次聚类,最终将川滇地区划分为6个块体：藏东块体（BYw）、阿坝块体（BYe）、川滇菱形块体（CD）、华南块体（SC）、腾冲-保山块体（DSw）和印支块体（DSe）。其中,CD块体的西、北边界分别是“三江并流”三江断裂区域和巴塘地区的一条速度阶跃线,且巴颜喀拉主峰断裂将巴颜喀拉块体（BY）分为BYw和BYe两个块体。本文给出的川滇地区块体划分模型是数学意义上的一种合理方案,与基于地质资料的结果相吻合,可为该地区的块体划分提供另一维度的解释。     

汶川地震震后形变的GPS反演

通过对汶川地震后震区附近的GPS连续站资料的高精度数据处理,获得各测站的时间序列.利用时间序列分析,求得震后松弛时间约为38 d,并获得各测站的震后位移.基于粘弹性松弛模型,对震后位移进行了模拟,反演出龙门山地区地壳的弹性层厚度的最佳估值为45 km,粘弹性层的粘滞系数的最佳估值为1.8×1019 Pa·s,该结果与地震学以及地球物理学等给出的结果基本一致.     

利用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范数下可以有效抵抗粗差的干扰,充分利用已有先验信息的同时顾及观测、线性化等原因引起的系数矩阵误差的影响。     

利用PSInSAR研究意大利Etna火山的地表形变

利用2003~2010年升轨、降轨ASAR数据对意大利Etna火山的地表形变进行PSInSAR处理,获取了该火山区域的高精度形变速率场,在此基础上,对部分特征区域点位进行时序分析。结果显示,Etna火山的东部区域的形变量变化大,LOS方向形变速率约为7~12 mm/a,西部区域变化较小,LOS方向形变速率仅为1~3 mm/a;火山东南和南侧存在两个较大垂直形变区域,其垂直形变速率分别为-4 mm/a和12mm/a。目前,Etna火山东南侧区域地表形变起伏较大,位于东南区域下的火山源岩浆活动活跃、岩浆体积处于补充期。     

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

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