联合地震位错模型和InSAR数据构建2017年九寨沟Mw 6.5地震同震三维形变场

利用哨兵（Sentinel）-1A卫星升、降轨影像,在地震位错模型约束下获取了2017年九寨沟Mw 6.5地震的高质量三维形变场。首先,利用合成孔径雷达干涉测量技术（interferometric synthetic aperture radar,InSAR）提取九寨沟地震升、降轨同震形变场;然后,通过“两步法”反演获取该地震发震断层的几何参数和分布式滑动模型,以此为约束,采用方差分量估计算法联合解算九寨沟地震三维形变场。结果表明,九寨沟地震同震三维形变场以水平位移为主,垂向形变较弱;南北向形变呈拉张趋势,断层上盘向南、下盘向北滑动,最大位移分别为－19.81 cm和14.38 cm;东西向形变不对称性明显,断层上盘西北部向东水平运动,最大位移为18.37 cm,下盘东南部向西运动,最大位移不足8 cm。将南北、东西向形变与6个全球导航卫星系统（global navigation satellite system,GNSS）台站观测数据进行比较,两者一致性较好且均方根误差较小,分别为1.44 cm和1.77 cm,表明联合升、降轨InSAR观测和地震位错模型约束构建同震三维形变场方法具有较高可行性,显著降低了大地测量数据不足、InSAR观测对南北向形变不敏感等问题的影响。     

九寨沟Ms7.0级地震三维同震形变场提取

针对D-InSAR地震形变监测中存在LOS(line of sight)视线向模糊的问题,该文构建了一种融合升降轨不同视线向干涉测量数据获取三维同震形变场的方法。以2017年九寨沟M_S7.0级地震为例,基于研究区的Sentine1_1A数据和SRTM(1″)DEM数据,采用二轨差分的D-InSAR技术,融合升降轨LOS向以及自定义视线向的干涉测量数据进行联合解算,获取了研究区三维同震形变场,解决了单一轨道雷达LOS提取结果不能准确反映地表三维形变场的问题。实验表明,此次地震造成了震中附近一定范围地表的沉降和东南向的滑移,并通过震区形变剖面图和形变等值线图,分析了形变场的空间分布以及此次地震所造成地表断裂的位置,认为此次地震与塔藏断裂、虎牙断裂存在关联性,其运动形式为主动盘逆冲的走滑型地震。     

2017年伊朗Mw7.4地震三维同震形变场及滑动分布

针对如何准确分析2017年11月12日的伊朗Mw7.4地震断层属性的问题,该文从三维同震形变场与断层滑动分布角度出发,采用Sentinel-1A升降轨数据获取该地区三维同震形变场,并基于弹性半空间矩形位错模型,结合最速下降法反演得到其发震断层参数、模拟三维形变场。结果显示,地震东西向和南北向最大形变分别为68.8、43.8 cm,竖直方向上最大沉降和抬升分别达到34.4、96.7 cm,均符合逆冲断层的运动特征。地震主要形变特征为抬升,地震矩1.66×10²⁰ N·m,矩震级Mw7.4,断层倾角16°,走向351°。结果表明,Okada模型解算所得地震三维形变场与观测结果总体特征相吻合,并且该地震为西南部板块挤压东北部板块而引发,属于逆冲为主带有右旋走滑的逆冲断层,其与地震学结果一致。     

DInSAR技术监测玉树地震同震形变场的研究

近二十年来,合成孔径雷达干涉测量技术（DInSAR）作为一种监测地表形变的有效技术得到广泛应用,其视线向精度可达厘米级甚至毫米级。该技术尤其是对监测快速、激烈的地表形变更为有效,如地震、火山等。利用DInSAR技术,选用日本ALOS卫星PALSAR数据,获取了2010—04—14玉树Mw6．9地震的同震形变场。结果表明：地表至少发生过3次地表破裂,沿断层走向的形变分布范围远远大于垂直断层分布范围,形变分布特征以左旋走滑为主。     

地震破裂模型约束的中国阿里地震三维形变场

利用环境卫星（Environmental Satellite,Envisat）的升、降轨数据在地震破裂模型约束下获取2007年阿里地震的高质量同震地表三维形变场。首先,对合成孔径雷达（synthetic aperture radar,SAR）影像进行差分干涉处理,得到地震造成的视线向同震地表形变场;然后,以地震破裂模型为约束条件,采用赫尔默特方差分量估计法来解算阿里地震的高质量地表三维同震形变场。结果表明,震中区域最大下沉达约4.7 cm,东西向位移较小,南北向呈挤压趋势。总体上,三维形变的特征表明阿里地震是一个以正倾滑为主兼少量右旋走滑运动的地震事件。     

唐山7.8级地震三维地表形变和断层运动特征

通过系统收集整理1976年唐山7.8级地震前后区域水准和边角网观测资料,在研究发展动态平差程序的基础上,解算得到了唐山地震震前和同震三维形变定量结果与误差分布,并反演得到了断层面同震滑动分布。结果显示,震前阶段丰润-唐山-唐海-昌黎-迁安一带整体隆升,隆升高值区位于唐山附近,震中区的水平形变量不大,唐山菱形地块相对于周围地块逆时针扭转;同震阶段唐山断裂以右旋错动为主,并呈现西升东降的垂向变形特征;反演结果表明多段断层参与了唐山7.8级地震的应变释放,主震发震断层的南段破裂主要发生在深度约6~18 km、长轴约50 km范围,北段破裂主要发生在深度约7~17 km、长轴约30 km范围,滦县地震断层也发生了少量滑动;主震断层最大水平位错为7.82 m、最大垂直位错为2.04 m,等效震级为Mw 7.58。挖掘得到的唐山震前和同震三维形变场定量结果为深入研究该地震的孕育发生机制及地震影响提供了珍贵资料。     

DInSAR技术对地震同震形变场的研究

合成孔径雷达差分干涉测量（DInSAR）可用于监测厘米级或更微小的地表形变,以揭示许多物理现象,如地震形变、火山运动、大气变化、冰川漂移、地面沉降以及山体滑坡等。DInSAR作为一种新型的空间对地观测技术,具有不受时间和空间的限制、对地物具有一定的穿透性等传统测量所不可比拟的优势,已得到较为广泛的应用。简要介绍DInSAR技术的基本原理及其处理流程,以Barn地震为例提取Bam地震的同震形变场。     

InSAR同震地表三维形变反演：一种顾及形变梯度的联合解算方法

同震地表三维形变场是对地震引起的地表形变在真实空间中最直观的描述,在地震形变研究中,同震地表三维形变场的恢复具有重要意义。针对目前地表三维形变解算方法存在的不足和缺陷,提出了一种顾及形变梯度的同震地表三维形变场联合解算方法。该方法利用形变梯度信息将合成孔径雷达干涉（interferometric synthetic aperture radar,InSAR）观测的形变区域分解成形变梯度近似相等的互不相交的子区域,然后根据子区域内观测值数量选择加权最小二乘法（weighted least square,WLS）或基于应力-应变模型（stress-strain model,SM）和方差分量估计（variance components estimation,VCE）的InSAR三维形变解算方法（SM-VCE）获取同震地表三维形变场。模拟实验和实际震例研究结果表明,相比于传统的WLS和SM-VCE方法,通过形变梯度信息能有效恢复地震破裂带的地表三维形变,得到更加完整和可靠的同震地表三维形变场。     

融合多平台DInSAR数据解算拉奎拉地震三维同震形变场

针对单一平台DInSAR技术仅能获取雷达视线方向同震形变场的问题,根据雷达成像的几何条件,融合不同轨道、不同平台的DInSAR数据解算了拉奎拉地震的三维同震形变场。三维形变结果反映的拉奎拉地震发震断层的特征与地质调查的结果较吻合。将得到的三维形变场数据与该地区GPS观测站数据进行比较,结果表明,得到的拉奎拉地震的三维同震形变场比较可靠且精度较高。     

利用InSAR及其差分技术获取仲巴地震的同震形变场

本文简要介绍了InSAR及其差分技术的基本原理,并以仲巴地震为例,采用已知DEM的二路差分法获取了仲巴地震的同震形变场。结果表明:仲巴地震为是一个以倾滑为主的浅源正断层型地震,断层面走向近乎南北方向。     

宜宾地震同震形变场双轨法D-InSAR技术的研究

为了研究2019年6月17日宜宾Ms6.0级地震引起的地表形变量,在欧空局下载覆盖研究范围内的Sen-tinel-1A卫星C波段雷达影像数据,使用SRTM(3″)DEM,采用双轨法差分干涉测量技术进行地表形变信息的提取,对比不同相位解缠方法对研究结果精度的影响,最终得到宜宾地震在视线方向上的形变图、形变量.结果表明:双轨法D-InSAR技术在地表形变信息提取的应用中,使用合适的解缠方法能够提高成果的精度.     

西藏改则地震升降轨同震形变场特征分析与破裂模式

通过升降轨ASAR数据的获取与差分干涉处理,获取了西藏改则地震的双视线向同震形变场信息。对升降轨同震形变场的特征分析表明,改则地震主、余震均表现为典型的正断裂破裂模式,并先后形成了东、西两条规模不同的走向NE、倾向NW的正断层破裂带(可能未出露地表),以及西北盘的东、西两个沉降形变中心。从干涉条纹的切割关系判断,东沉降形变中心可能受2008年1月9日Ms 6.9级主震的控制,西沉降形变中心可能受2008年1月16日Ms 6.0级余震的控制。Ms 6.9级主震形成东南盘隆升与西北盘沉降的总体格局;Ms 6.0级余震使西北盘余震震中及以西部位进一步沉降,造成西沉降中心沉降形变量大于东沉降中心。升降轨同震形变场的获取为进一步的震源参数模拟与三维形变场解算提供了更好的约束条件与信息。     

汶川地震GPS形变约束的破裂分段特征及滑移

为了解2008年汶川地震破裂分段特征及滑移,采用弹性位错模型和模拟退火算法,数值模拟汶川震区密集的GPS同震形变。结果表明,GPS同震形变场至少需要用铲状的映秀-北川五断裂加上灌县－江油断裂来模拟,该模型对GPS数据的符合程度与汶川地震滑动分布模型相当,基本反映了汶川地震破裂特征。映秀－北川断裂总长255~294km。南段以逆冲为主,分上下两层,最深达30km。中段右旋走滑和逆冲都很大,而北段以右旋走滑为主。映秀段地震周期最短为3000年。破裂上层深度14km可能是龙门山中央断裂的闭锁深度。汶川地震可能源于相邻块体的相对运动和挤压、深部滑脱或浅部闭锁,在薄弱构造处首先爆发地震。     

地形变观测揭示的汶川Ms8.0级地震同震响应与影响

利用GPS、跨断层长水准、跨断层短水准与基线等多手段、不同时空尺度地形变观测资料,研究了2008年5月12日四川汶川县Ms8.0级地震同震响应及对周围地区的影响。分析表明了汶川地震发震区呈相当显著的逆冲兼右旋走滑同震变化,大震对其附近的甘川交界东段——甘川陕交界区、鲜水河断裂南段及与龙门山断裂、安宁河断裂的Y字型交汇区及其附近影响相对明显,对其应变积累可能以一定程度增强作用为主。     

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

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

利用InSAR研究西藏拉孜地震同震形变

采用欧空局ERS—1卫星雷达数据获取了1993年Mw 6.2级西藏拉孜地震同震形变场,并采用弹性位错模型反演了该地震的断层参数。结果显示该地震产生的最大地表位移约为12 cm,断层倾角为42°,不支持该断层为低角度正断层的论断。     

利用SAR影像配准偏移量提取地表形变的方法与误差分析

单一的InSAR观测技术可提取地表沿雷达视线方向(LOS)上的一维位移,而利用SAR影像配准过程中的同名像素偏移量可提取地表沿雷达方位向(近南北向)与距离向(近东西向)的二维形变场,与LOS方向的一维形变形成优势互补。本文在分析SAR影像配准偏移量提取地表形变场方法的基础上,推导建立了雷达方位向与距离向形变提取的误差模型,探讨了该方法提取地表形变的主要误差源。以Bam地震ASAR影像和玉树地震PALSAR影像为数据源,分别开展同震形变场的提取与误差分析试验,结果表明,基于SAR像素配准偏移量提取同震形变场的精度主要受匹配窗口尺寸与过采样因子影响,形变提取误差随匹配计算窗口的增大而减小,形变提取精度随过采样因子的增大有适量提高,地形起伏效应在高差较大的SAR影像区域中表现较为显著。