基于Sentinel-1A的新疆阿克陶MS6.7地震同震形变与滑动分布特征

文中利用Sentinel-1A升、降轨数据获取了2016年新疆阿克陶地震的同震形变场。形变场以水平运动为主,形变主要发生于断层南盘;升、降轨最大形变量分别约为12cm和-21cm;基于雷达影像观测右视成像的特点可知阿克陶地震具有右旋走滑的破裂特征,结合形变场形态特征与余震剖面推断认为,阿克陶地震的发震断层为S倾的木吉断裂。基于均匀弹性半空间模型反演双断层面的静态滑动分布,结果显示:2个断层面上各存在1个椭圆状的滑动集中区,破裂基本到达地表,滑动主要发生在沿断层面走向长约50km、沿断层面倾向方向0～20km的范围内,最大滑动量位于约9km深处,量级约为0.7m;西段主要以走滑为主,东段为走滑兼具少量正断性质,反演得到的地震矩约为8.81×10¹⁸N·m,相当于M_W6.57。综合分析认为,本次地震的震源特征为右旋走滑兼有少量正断分量,是发生在公格尔拉张系内拉张环境下的构造地震事件。     

2016年新疆阿克陶地震Sentinel-1 InSAR同震形变特征

2016年11月25日新疆克孜勒苏州阿克陶县发生M_W6.6地震。 本文利用合成孔径雷达差分干涉测量技术, 对Sentinel-1卫星获取的升、 降轨雷达数据进行了处理, 提取了该次地震的同震形变场, 并结合形变场特征与震源机制解, 采用梯度下降法反演发震断层的滑动分布。 结果表明, 升、 降轨LOS向同震形变场在发震断层两侧具有明显不同的形变特征, 主要形变区域分布在断层两侧, 升轨LOS向形变量可达-8.2 cm与11.2 cm, 降轨LOS向形变量可达-21.4 cm与13.1 cm; 反演的升、 降轨干涉形变场与InSAR测量值之间的残差得到有效控制, 大部分的残差介于±5 cm之间; 断层滑动分布主要集中于沿断层面深约2~18 km处, 最大滑动量位于沿断层面深约7 km处可达0.96 m; 平均滑动角约182.29°, 最大滑动处的滑动角约197.13°, 两个滑动分布中心的滑动角均接近180°, 表明阿克陶地震为一典型的右旋走滑破裂性事件; 当剪切模量取32 Gpa时, 反演的发震断层地震矩M₀可达9.75×10¹⁸, 相当于矩震级M_W6.60, 与地震波形反演结果一致。     

基于InSAR分析2015年尼泊尔MW7.8地震形变场特征及断层滑动分布反演

2015年4月25日尼泊尔发生了M_W7.8地震, 本文基于震前、 震后两景Sentinel-1A雷达影像, 采用D-InSAR两轨差分干涉法提取了此次地震的同震形变场。 结果显示, 同震形变场位于喜马拉雅造山带—主边界逆冲断裂(MBT)和主前锋逆冲断裂(MFT)附近, 形变场整体表现为自西北向往东南方向延伸近150 km的纺锤形包络状, 以大面积隆起抬升形变为主, 视线向最大隆升形变达1.18 m, 抬升区北侧存在一小沉陷区, 以InSAR观测值定位同震最大形变中心。 基于均匀介质弹性半空间模型(Okada模型)与InSAR观测数据反演了断层滑动分布。 反演结果表明该地震属于典型逆冲型地震, 发震断层为主喜马拉雅逆冲断裂(MHT), 同震破裂从主喜马拉雅逆冲断裂(MHT)向上沿着主前锋逆冲断裂(MFT)传递。 基于InSAR同震形变场局部形变细节, 结合震区地质背景、 断裂分布及断层运动特征, 获得了同震破裂拟出露地表迹线。     

基于InSAR技术的伽师MS 6.4地震同震形变与滑动反演研究

<正>1研究背景地震震级较大或震源较浅时,地表会产生一定程度的变形,通常称之为同震形变。及时、准确和全面地获取同震形变场,能够快速认识地震地表形变整体宏观特征,甚至有助于研判发震断层位置及其走向变化等,这对政府应急响应和灾后救援部署等至关重要,具有重大的社会效益和经济效益（温少妍等,2020）。合成孔径雷达干涉测量技术（InSAR）是在干涉雷达基础上发展起来的一种微波遥感技术,具有形变灵敏度高、高空间分辨率高、覆盖范围大、几乎不受云雨天气制约等突出的技术优势,在对同震形变等大范围缓慢形变的测量研究中迅速得到广泛应用。     

基于InSAR形变观测反演2015年皮山MS6.5地震滑动分布

2015年7月3日在新疆皮山县发生了M_S6.5地震, 该地震使当地遭受了巨大的经济损失。 本文利用欧空局提供的Sentinel-1A卫星差分干涉数据对该地震的震源机制情况进行了反演研究, 首先运用两轨法对卫星雷达影像进行差分干涉处理, 获取了覆盖皮山地震震区的同震形变场, 然后利用弹性半空间的均匀滑动模型反演获取了发震断层的几何参数, 并对原始观测数据进行降采样处理, 在此基础上运用分布式滑动模型反演获取了更为精细的断层滑动分布, 结果显示分布式滑动模型与观测结果有很高的拟合度。 反演结果表明发震断层是一个以逆冲为主兼有极少量左旋走滑的盲断层, 此次地震断层面的同震活动分布主要集中在7~15 km深度范围内, 同震的地震矩为6.28×10¹⁸N·m, 矩震级为M_W6.46, 与前人的研究结果非常一致。     

基于Sentinel-1A数据的四川长宁MS6.0地震同震形变场分析及断层滑动分布反演

利用大地测量数据研究2019年6月17日四川长宁M_S6.0地震同震形变场特征和发震断层参数, 基于DInSAR技术处理升降轨Sentinel-1A数据获取干涉相位图, 并考虑大气折射效应和余震形变误差实现同震形变场改正。四叉树采样后的形变数据作为反演数据源, 采用弹性半空间分层模型反演发震断层几何面滑动分布。结果表明本次地震发震机制为兼具逆冲和左旋走滑, 矩震级为M_W5.9, 断层破裂尺度达28 km×20 km, 震源深度约9.4 km。升降轨视线向同震形变场在断层两侧呈现形变特征差异, 最大沉降量分别是8.34 cm(升轨)和4.23 cm(降轨), 最大抬升量分别是5.5 cm(升轨)和7.5 cm(降轨); 发震断层走向为302°, 倾角为43°, 平均滑动角为50°, 断层面最大滑动量达到0.28 m。     

利用InSAR技术研究2016年青海门源MW5.9地震同震形变场及断层滑动分布

2016年1月21日, 青海省门源县冷龙岭断裂带附近发生了M_W5.9地震。 基于Sentinel-1A影像, 采用差分干涉雷达测量技术研究了此次地震产生的同震形变场, 结果表明, 门源地震的形变影响范围约20~30 km, 形变态势在升降轨道形变场均显示为隆升, 基本沿冷龙岭断裂呈近似同心圆展布, 推测可能是冷龙岭断裂与民乐—大马营断裂之间的一条逆断层, 沿雷达视线方向最大形变量级约为6 cm。 均匀滑动反演显示门源发震断层长7.3 km, 宽6.2 km, 走向298.6°, 倾角34.5°, 倾向宽度9.5 km, 沿走向滑动量为170 mm, 沿倾向滑动量为460 mm, 矩震级为M_W5.97; 分布式滑动反演显示门源地震以逆冲为主, 兼具少量右旋走滑分量, 滑动量主要集中在沿断层倾向方向, 距离地表5～15 km处, 最大滑动量约0.3 m, 位于断层倾向深度10 km处, 矩震级为M_W5.93。     

2016年阿克陶MS6.7地震震源复杂性与烈度

2016年11月25日在我国新疆克孜勒苏州阿克陶县发生M_S6.7地震（阿克陶M_S6.7地震）.我们收集国内外地震资料,对主震及4级以上余震进行了重新定位和震源机制反演,对434次余震进行了双差定位,对主震震源过程进行了反演确定和复杂性分析,并基于反演确定的有限动态源模型估计了此次地震的烈度分布.结果表明：这次地震发生在当地一个近乎东西向展布的小型盆地内,很可能由一条新断层或隐伏断层的活动所致.发震断层近乎直立,近东西向展布,总体上表现为右旋走滑.破裂首先向西扩展,紧接着向东,随后向东西两个方向同时扩展,然后西侧破裂首先停止,东侧破裂继续,最后破裂在东侧停止,整个过程持续~20 s,释放地震矩1.08×10¹⁹N·m,相当于M_W6.6.破裂过程最终形成两个位错高值区,分别位于初始破裂点的东西两侧,西侧高值区规模较小,东侧区规模较大.根据烈度估计,烈度椭圆长轴方向与主震破裂方向以及余震展布方向一致,最大烈度约为Ⅸ度,主要集中在震中以东很小的区域,Ⅷ度区呈纺锤形,分布于震中东西两侧,Ⅴ至Ⅶ度区呈椭圆形,总体上东侧烈度大于西侧.

     

2016年阿克陶MW6.6地震的同震断层滑动分布及库伦应力变化

2016年11月25日新疆阿克陶地区发生的M_W6.6地震,对当地人民生命造成一定的危害。基于Sentinel-1数据获取该地震的视线向同震形变场,采用贝叶斯方法反演单一断层走向为106.9°、倾角为73.8°、震源深度为17.35 km,在此断层几何模型基础上,以最速下降梯度法（Steepest Descent Method,SDM）反演滑动分布,结果表明断层面上存在两个滑动峰值,其中位于断层西侧的最大滑动量为0.66 m,深度为11.7 km,位于断层东侧的最大滑动量为0.83 m,深度为7.5 km,根据反演结果模拟LOS形变,其最大残差为~5 cm。构建倾角分别为70.79°和55.33°的双断层几何模型,并根据双断层几何模型反演了滑动分布,结果表明单一断层模型与双断层模型的滑动分布具有一致性,但是最大滑动量值有所不同,相对于单一断层模型的滑动分布而言,双断层模型的滑动量在西侧增大,其值为0.68 m;而在东侧减小,其值为0.77 m;最大残差降低了约2 cm。双断层模型库伦应力增加区域与余震的分布比较吻合。     

新疆阿克陶MS6.7地震引起的新10井水位同震响应研究

本文根据新10井数字化高频采样水位仪记录到的2016年11月25日新疆阿克陶M_S6.7地震所引起的水震波,对比分析了该井水位与地表垂向运动的相关性特征,并对二者与井-含水层系统水文参数的关系进行了深入探讨.分析结果显示:①与地震波信号相似,新疆阿克陶M_S6.7地震引起的新10井水震波存在两个显著的周期,即6—10 s和15—30 s;②新10井水震波响应幅度与地表垂向运动幅度整体呈正相关,且在高频阶段(频率大于0.08 Hz)二者的振幅比随着频率的减小而增大,表明该井水位对周期大于12 s的信号放大效能较高;③利用水震波与地震波的振幅比估算新10井观测含水层渗透系数的量级为10-2 cm/s,且在地震波作用过程中含水层的水文参数也存在波动.本研究表明,井水位的同震响应机理较为复杂,在分析水位同震响应特征时,高采样率的水位数据是获得可靠结果与认识的基础.      

2008年于田MW7.2地震滑动分布反演及三维同震形变场解算

基于ENVISAT ASAR升降轨数据, 利用InSAR获取2008年于田M_W7.2地震同震形变场; 采用SDM反演本次地震断层滑动分布; 使用PSCMP正演获取于田M_W7.2地震南北向模拟形变量, 并结合升降轨同震形变场, 解算三维同震形变场。 同震形变场分析表明, 2008年于田M_W7.2地震以正断为主, 且带有走滑运动特征, 破裂带走向为NNE向。 同震滑动分布反演结果显示, 断层沿走向被分为4段F1、 F2、 F3、 F4, 其滑动分布集中在0~14 km区间, 以F2、 F3段为主, 最大滑动量约5.31 m, 位于F2段深部2.76 km处; 沿破裂带走向, 左旋走滑位移与垂直位移比值有增大的趋势; 反演获得的地震矩M₀=5.58×10¹⁹N·m, 相当于矩震级M_W7.1。 三维同震形变场解算结果显示, 断层上盘整体表现为沉降, 断层下盘整体表现为隆升, 且沉降量明显大于隆升量, 表明地震以正断破裂为主; 除靠近断裂带中上部表现为向东南运动外, 上盘整体上表现为向西南运动; 断层下盘则整体表现为向东北运动, 证明破裂兼有左旋走滑运动。 滑动分布反演、 正演与三维同震形变场解算结果皆表明, 于田M_W7.2地震破裂以正断为主, 且带有一定的左旋走滑。     

2016年11月25日新疆阿克陶6.7级地震序列震源机制特征分析

本研究利用新疆区域数字地震台网的波形资料,采用CAP方法反演了2016年11月25日阿克陶6.7级地震的前震、主震及11次M_S ≥ 3.6余震序列的最佳双力偶震源机制解,得到阿克陶6.7级地震最佳双力偶机制解：节面Ⅰ走向20°/倾角69°/滑动角-10°;节面Ⅱ走向114°/倾角81°/滑动角-159°,表明此次阿克陶6.7级地震为一次走滑型地震事件,结合震源区的地震地质构造及余震序列空间分布等已有研究成果,判定节面Ⅱ代表了主震的发震断层面。主震最大主压力轴方位为339°,与震源区附近历史中强震P轴近NW向的优势方位基本一致。其4.8级前震的震源机制解为走滑型,与主震震源机制解具有较高的一致性。11次余震中有6次为走滑型地震,3次为逆断型地震,1次正断型地震,1次混合型地震,且多数地震具有近NW向的P轴方位。此次6.7级地震序列的震源深度分布于6~16km之间,而大部分地震为9~13km,与本文计算得到的主震的震源深度10km相差不大。此外,初步分析了兴都库什-帕米尔地区强震活动与此次阿克陶6.7级地震的关系。     

Co-seismic Slip Distribution of the 2001 Kokoxili Earthquake Inverted by the GPS Data

INTRODUCTIONThe Kokoxili earthquake,which ruptured the East Kunlun fault in northern Tibet on November14,2001,is the largest earthquake that has occurred on the Chinese continent in the past50yearsever since the Moto,TibetM8.0earthquake of August15,1950.The macro hypocenter of thisearthquake is located in the north of Bukadaban Peak at36.2°N,90.9°E,MS=8.1(ChineseSeismological Network),or the west of Bukadaban Peak at36.01°N,90.49°E,MW=7.9(USGSNEIC),or near the Kusai Lake …     

2016年青海门源MS6.4地震震前应变积累及同震变形特征

2016年1月21日青海省门源县发生了M_S6.4地震, 发震断裂为冷龙岭北侧断裂, 震中位置与1986年门源6.4级地震相同。 本文收集了本次地震震中及其周边区域1999—2015年GPS观测资料, 解算了GPS速度场、 跨断裂连续观测站基线时间序列和应变率场。 结果显示, 祁连山断裂带为一条宽度约60 km的连续变形带。 在断裂带南侧地壳运动以顺时针旋转为主, 运动量值没有显著差异; 跨过断裂带到达其北部之后, 地壳运动量值明显减小, 显示出该断裂带的强烈活动特征。 冷龙岭断裂左旋走滑速率为6.17±0.41 mm/a, 挤压速率为1.83±0.38 mm/a, 断层闭锁深度为22.1±3.1 km。 利用GPS连续观测站数据解算的地震同震位移显示, 震中西南侧26.8 km处的青海门源(QHME)测站记录到了明显同震位移, 其水平运动方向为北东向, 与逆冲为主的震源机制解一致。     

2020年伽师MS6.4地震InSAR同震形变与滑动分布特征分析

利用Sentinel-1A升轨和降轨数据,基于D-InSAR技术,获取2020年1月19日伽师M_S6.4地震同震形变场,并结合其他研究机构给出的震源机制解参数和已有研究成果,反演得到伽师地震的发震断层几何特征和滑动分布。研究结果表明,伽师地震同震形变在地表有明显差异;升轨同震形变在卫星视线方向北侧抬升55 mm,南侧下降42 mm;降轨同震形变在卫星视线方面北侧抬升63 mm,南侧下降23 mm。通过反演得到发震断层走向为275°,倾角为20°,地震滑动主要分布在地下5 km处,最大滑动量约为0.32 m,平均滑动角为89.3°,累积地震矩为1.46×1018 N·m,合矩震级M_W6.1,发震构造为具有少量走滑性质的逆冲断裂。从发震构造特征、同震滑动分布推测,伽师地震发震构造是柯坪塔格褶皱带滑脱面以上沉积盖层内的逆冲断裂,支持了柯坪推覆体的薄皮构造模型观点。     

利用考虑地形起伏的反演方法获取2017年伊拉克MW7.3地震的断层参数及滑动分布

2017年11月13日伊拉克北部地区苏莱曼尼亚省发生了M_W7.3地震, 造成了重大的人员和经济损失。 本文利用升降轨的Sentinel-1和降轨的ALOS-2卫星的SAR数据, 通过差分干涉测量技术获取了该地震的同震形变场, 联合DInSAR和MAI技术, 采用抗差最小二乘法求解该地震的同震三维形变场。 基于改进的考虑地形起伏的均匀位错模型反演确定了发震断层的断层参数, 最后基于非均匀位错模型得到了发震断层的分布式滑动分布模型。 结果显示: ALOS-2卫星降轨轨道观测到的伊拉克地震引起的LOS向地表形变最大为55.8 cm抬升和47.9 cm下沉; Sentinel-1卫星观测到的伊拉克地震引起的LOS向地表形变为: 升轨轨道最大为87.9 cm抬升和17.1 cm下沉; 降轨轨道最大为55.6 cm抬升和38 cm下沉; 相对于前人的研究成果, 本文利用改进的考虑地形起伏的反演方法得到的发震断层几何参数表明发震断层为NNW走向, 倾向角为352°, 同震破裂以逆冲为主, 同时兼有一定的左旋走滑分量。 基于均匀位错模型反演得到的断层滑动分布结果表明, 同震破裂未延伸至地表, 主要滑动量集中在12~18 km, 最大滑动量位于15 km深度, 达到4.3m, 反演得到的矩震级为M_W7.35, 与UGSG、 GCMT等机构给出的结果一致。