文章快速检索     高级检索
  大地测量与地球动力学  2023, Vol. 43 Issue (3): 221-225  DOI: 10.14075/j.jgg.2023.03.001

引用本文  

颜丙囤, 殷海涛, 冯兵, 等. InSAR数据约束下的2022年泸定MS6.8地震震源参数及滑动分布[J]. 大地测量与地球动力学, 2023, 43(3): 221-225.
YAN Bingdun, YIN Haitao, FENG Bing, et al. Source Parameters and Slip Distribution of the Luding MS6.8 Earthquake in 2022 Constrained by InSAR Data[J]. Journal of Geodesy and Geodynamics, 2023, 43(3): 221-225.

通讯作者

殷海涛,博士,研究员,主要从事GNSS地壳形变分析与地球动力学应用研究,E-mail:yinhaitao121@163.com

Corresponding author

YIN Haitao, PhD, researcher, majors in crustal deformation of GNSS and geodynamics application, E-mail: yinhaitao121@163.com.

第一作者简介

颜丙囤,工程师,主要从事遥感地震监测研究,E-mail:806056236@qq.com

About the first author

YAN Bingdun, engineer, majors in remote sensing earthquake monitoring, E-mail: 806056236@qq.com.

文章历史

收稿日期:2022-09-16
InSAR数据约束下的2022年泸定MS6.8地震震源参数及滑动分布
颜丙囤1     殷海涛2     冯兵3     冯恩国1     刘保华1     李峰1     
1. 山东省地震局聊城地震监测中心站,山东省聊城市中华南路24号,252001;
2. 山东省地震局,济南市文化东路20号,250014;
3. 中国地震局第二监测中心,西安市西影路316号,710054
摘要:基于Sentinel-1A升降轨数据,使用D-InSAR反演泸定MS6.8地震同震形变场。首先基于弹性半空间位错模型反演地震震源参数,然后利用分布滑动模型确定断层面上的滑动分布,最后计算周边断层的静态库仑应力变化。结果表明,此次地震造成的雷达视线向最大形变量约为18 cm,同震位错以左旋走滑为主,主要发生在6~24 km深度,最大滑动量约为2.5 m,位于16 km深度,发震断层主要为鲜水河断裂。静态库仑应力触发关系显示,鲜水河断裂南段和安宁河断裂北端的库仑应力显著增强,强震趋势值得关注。
关键词同震形变Sentinel-1AInSAR泸定地震震源机制库仑应力变化

据中国地震台网测定,2022-09-05四川泸定县发生6.8级地震,震中位置29.59°N、102.08°E,震源深度16 km。此次地震发生在鲜水河断裂带附近,初步震源机制解显示为走滑型破裂。鲜水河断裂带位于青藏高原东部,是一条大型左旋走滑断裂,由甘孜向南经侏倭、炉霍、道孚、康定,最后延伸至泸定的磨西以南部分,呈略向NE向凸出的弧形[1-2]。由于采用的数据和方法不同,各机构和团队给出的断层震源机制解存在差异。张喆等(https://www.cea-igp.ac.cn/kydt/279423.html)利用矩心矩张量反演方法获得的最优模型结果显示,此次地震的标量矩震级为MW6.61,矩心深度为16.0 km;王卫民等(http://www.itpcas.cas.cn/new_kycg/new_kyjz/202209/t20220906_6509485.html)利用地震数据资料研究震源机制和震源破裂过程,得出此次地震类型为高倾角走滑型,地震矩为MW6.65,最大滑动量为192 cm,震源深度约13 km;王洵等(https://www.ief.ac.cn/kydts/info/2022/69528.html)采用IRIS远场波形资料进行研究,认为此次地震矩震级为6.7级,初始破裂点深度约为12.8 km,为高倾角左旋走滑型地震。USGS和GCMT的详细震源参数信息见表 1,区域构造背景见图 1

表 1 泸定MS6.8地震震源参数 Tab. 1 Source parameters of the Luding MS6.8 earthquake

图 1 区域构造背景 Fig. 1 Regional tectonic background
1 InSAR数据处理与同震形变场 1.1 InSAR数据处理

选用Sentinel-1A中C波段SAR数据,极化方式为VV[3],同时选择短时间基线,有助于保持干涉图的相干性,图像参数见表 2。首先采用SARscape软件进行数据处理,然后采用Goldstein滤波方法对去地形后的干涉图进行滤波处理,再采用最小费用流法实现绝对相位恢复,最后将SAR坐标系统结果投影到WGS-84地理坐标系中[3]

表 2 升降轨干涉影像参数 Tab. 2 Interference image parameters of ascending and descending
1.2 InSAR同震形变场

D-InSAR反演的泸定地震升降轨InSAR干涉图和升降轨同震形变场见图 23。由图 2可见,发震区干涉条纹失相干现象严重。由图 3可见,形变场主要分布在震中区域,具有非对称分布特征,影响范围达60 km。其中,升轨干涉图的同震形变场显示,NW方向存在一个长约30 km的近长方形LOS向下沉区域,SE方向存在一个直径约18 km的近圆形上升区域,T26轨道显示最大沉降量约为10 cm,最大上升量约为12 cm。降轨干涉图与升轨干涉图的形变趋势相反,升轨干涉图的下沉区域在降轨干涉图中表现为上升,最大形变量约为9 cm;升轨干涉图的上升区域在降轨干涉图中表现为下沉,最大形变量约为18 cm[4]。根据升降轨成像几何模式可知,形变场的上、下盘呈现出不同的运动状态,说明地震引起的地表形变主要以水平运动为主,符合走滑断层地震的运动特征[5-6]。发震区域InSAR数据失相干严重,上升轨道近场信号噪声较大,升轨干涉对的空间基线为203.7 m,降轨干涉对的空间基线为49.3 m,升轨同震形变场存在形变异常区。为确定该形变异常区是否由相位解缠误差引起,参考韩炳权等[7]的多时相干涉结果发现,干涉对均存在上述形变异常区,说明相位解缠误差并不是引发形变异常区的因素。同时,结合光学影像数据发现,该异常区位于沟谷地带,区域形变信号与地形具有较强的相关性,因此形变异常区极可能是区域地形引起的大气误差(如聚集在沟谷上方的团雾等)所致[7]。基于此,本文最终采用降轨干涉对的形变场作为同震建模的数据源。

图 2 InSAR升降轨干涉条纹图 Fig. 2 InSAR interference fringe pattern of ascending and descending

图 3 InSAR升降轨同震形变场(LOS方向) Fig. 3 InSAR coseismic deformation field of ascending and descending (LOS direction)
2 InSAR震源机制反演和发震构造及区域地震危险性分析 2.1 InSAR震源机制反演

本文在InSAR同震形变分析结果的基础上,首先利用Okada弹性半空间位错模型反演发震断层的几何参数[7-9],然后基于约束条件下的最小二乘原理及最速下降法进一步反演震源断层破裂面上的滑动分布特征。

由于InSAR同震形变场数据量庞大,在空间上高度相关,因此在获取断层几何参数之前,需要对干涉图进行降采样处理,获取合适的InSAR形变数据集[3]。在采样过程中,为了增加震中附近区域的采样点密度、减少远离震中区域的采样点密度,采用均匀矢量网格对降轨形变场进行降采样处理,共得到5 680个降轨形变数据点(图 4)。基于降采样InSAR视线向形变数据的点位形变场结果,结合GCMT推荐的断层几何参数设置搜索间隔,断层走向设置为120°~200°,断层长宽和位置根据InSAR形变场特征确定,断层顶部深度设置为0~20 km,断层滑动角设置为-30°~30°,断层倾角根据区域构造背景设置为30°~90°。将非线性反演得到的均匀滑动断层几何参数引入到线性反演模型中,均匀反演得到的断层参数见表 3。降轨形变场采样及反演残差见图 5。由图可见,模拟结果与反演结果的拟合度较好。

图 4 降轨形变场采样点 Fig. 4 Sampling points of descending deformation field

表 3 均匀滑动反演断层参数 Tab. 3 Fault parameters by uniform slip inversion

图 5 降轨形变场采样及反演残差 Fig. 5 Sampling and inversion residuals of descending deformation field

为获取断层面上的精细滑动分布,将固定滑动角设置为169°,矩形断层面的长度(沿走向)和宽度(沿倾向)分别延长至80 km和50 km,走向和倾角细分为2 km×2 km,阻尼系数设置为0.02[3],所得结果见图 6。由图可见,断层同震错动以左旋走滑为主,同时具有少许逆冲分量。地震宏观震中为102.28°E、29.59°N,矩震级为MW6.9,同震滑动主要集中在6~24 km深度,断层最大滑动量约为2.5 m,位于地下约16 km处。断层在近地表处滑动量较小,说明同震错动未破裂至地表[6]。由表 1可见,张喆等、王卫民等、王洵等的研究成果与本文InSAR反演结果相差不大,证明了本文滑动分布的合理性。

图 6 同震滑动分布 Fig. 6 Coseismic slip distribution
2.2 发震构造与区域地震危险性分析

获得的同震形变场和断层滑动分布表明,青藏高原东南缘的顺时针旋转对区域断层运动方式具有控制作用。从构造单元上看,泸定地震发生在鲜水河断裂带西侧,震中所处的构造力学环境为NWW-SEE向挤压、NNE-SSW向拉张的走滑型应力结构,反演确定的震源机制与鲜水河断层的运动性质一致,初步判定此次泸定地震发震断层为鲜水河断裂。

强震发生后,区域应力场随同震断层错动发生改变,引起的库仑应力变化对余震的时空分布研究具有参考意义。本文基于分层半无限空间粘弹性地壳介质模型[6],运用Coulomb3.3软件包进行库仑应力分析。处理过程中,将泊松比设置为0.25,摩擦系数设置为0.4,分别计算10 km、16 km深度处泸定地震引起的库仑应力变化[5]。静态库仑应力计算结果显示(图 7),鲜水河断裂南段和安宁河断裂北端的库仑应力显著增加,鲜水河断裂带上应力加载最大值达0.5 bar,库仑应力的改变量大于0.1 bar,对余震发生的位置有所影响,库仑应力增加区域的地震危险性更大[6]。因此,鲜水河断裂南段和安宁河断裂北端的强震趋势值得关注。

图 7 泸定地震引起的邻近区域静态库仑应力变化 Fig. 7 Static Coulomb stress changes in neighboring regions caused by the Luding earthquake
3 结语

1) 升轨同震形变场因干涉对失相干严重,存在形变异常区。将降轨形变场作为同震建模的数据源,结果表明,形变场以左旋走滑为主,同时具有少量逆冲分量,形变主要位于6~24 km深度,最大滑动量约为2.5 m,位于16 km深度处,反演的矩震级为MW6.9。其他地震学方法反演得到的震源机制与本文反演结果相近。

2) 从区域应变分配和构造动力背景上看,泸定地震发生在鲜水河断裂带西侧。结合发震断层的几何产状和运动性质可知,反演确定的震源机制与鲜水河断层的运动性质一致,因此认为此次泸定地震的发震断层为鲜水河断裂。

3) 此次泸定地震引起的不同深度上的库仑应力变化结果表明,鲜水河断裂南段和安宁河断裂北端的强震趋势值得关注。

参考文献
[1]
闻学泽. 四川西部鲜水河-安宁河-则木河断裂带的地震破裂分段特征[J]. 地震地质, 2000, 22(3): 239-249 (Wen Xueze. Character of Rupture Segmentation of the Xianshuihe-Anninghe-Zemuhe Fault Zone, Western Sichuan[J]. Seismology and Geology, 2000, 22(3): 239-249) (0)
[2]
徐锡伟, 闻学泽, 郑荣章, 等. 川滇地区活动块体最新构造变动样式及其动力来源[J]. 中国科学: 地球科学, 2003, 33(增1): 151-162 (Xu Xiwei, Wen Xueze, Zheng Rongzhang, et al. The Latest Structural Change Style and Its Power Source of Active Blocks in Sichuan Yunnan Region[J]. Science China: Earth Sciences, 2003, 33(S1): 151-162) (0)
[3]
颜丙囤, 季灵运, 蒋锋云, 等. InSAR数据约束的2022年1月8日青海门源MS6.9地震发震构造研究[J]. 地震工程学报, 2022, 44(2): 450-457 (Yan Bingdun, Ji Lingyun, Jiang Fengyun, et al. Seismogenic Structure of Menyuan, Qinghai MS6.9 Earthquake on January 8, 2022 Constrained by InSAR Data[J]. China Earthquake Engineering Journal, 2022, 44(2): 450-457) (0)
[4]
季灵运, 刘传金, 徐晶, 等. 九寨沟MS7.0地震的InSAR观测及发震构造分析[J]. 地球物理学报, 2017, 60(10): 4 069-4 082 (Ji Lingyun, Liu Chuanjin, Xu Jing, et al. InSAR Observation and Inversion of the Seismogenic Fault for the 2017 Jiuzhaigou MS7.0 Earthquake in China[J]. Chinese Journal of Geophysics, 60(10): 4 069-4 082) (0)
[5]
李振洪, 韩炳权, 刘振江, 等. InSAR数据约束下2016年和2022年青海门源地震震源参数及其滑动分布[J]. 武汉大学学报: 信息科学版, 2022, 47(6): 887-897 (Li Zhenhong, Han Bingquan, Liu Zhenjiang, et al. Source Parameters and Slip Distributions of the 2016 and 2022 Menyuan, Qinghai Earthquakes Constrained by InSAR Observations[J]. Geomatics and Information Science of Wuhan University, 2022, 47(6): 887-897) (0)
[6]
徐晓雪, 季灵运, 朱良玉, 等. 漾濞MS6.4地震同震形变特征及发震构造探讨[J]. 地震地质, 2021, 43(4): 771-789 (Xu Xiaoxue, Ji Lingyun, Zhu Liangyu, et al. The Co-Seismic Deformation Characteristics and Seismogenic Structure of the Yangbi MS6.4 Earthquake[J]. Seismology and Geology, 2021, 43(4): 771-789) (0)
[7]
韩炳权, 刘振江, 陈博, 等. 2022年泸定MW6.6地震InSAR同震形变与滑动分布[J]. 武汉大学学报: 信息科学版, 2023, 48(1): 36-46 (Han Bingquan, Liu Zhenjiang, Chen Bo, et al. Coseismic Deformation and Slip Distribution of the 2022 Luding MW6.6 Earthquake Revealed by InSAR of Observations[J]. Geomatics and Information Science of Wuhan University, 2023, 48(1): 36-46) (0)
[8]
Okada Y. Internal Deformation Due to Shear and Tensile Faults in a Half-Space[J]. Bulletin of the Seismological Society of America, 1992, 82(2): 1 018-1 040 (0)
[9]
Ji L Y, Wang Q L, Xu J, et al. The 1996 MW6.6 Lijiang Earthquake: Application of JERS-1 SAR Interferometry on a Typical Normal-Faulting Event in the Northwestern Yunnan Rift Zone, SW China[J]. Journal of Asian Earth Sciences, 2017, 146(15): 221-232 (0)
Source Parameters and Slip Distribution of the Luding MS6.8 Earthquake in 2022 Constrained by InSAR Data
YAN Bingdun1     YIN Haitao2     FENG Bing3     FENG Enguo1     LIU Baohua1     LI Feng1     
1. Liaocheng Earthquake Monitoring Center of Shandong Earthquake Agency, 24 South-Zhonghua Road, Liaocheng 252001, China;
2. Shandong Earthquake Agency, 20 East-Wenhua Road, Ji'nan 250014, China;
3. The Second Monitoring and Application Center, CEA, 316 Xiying Road, Xi'an 710054, China
Abstract: Based on Sentinel-1A ascending and descending data and D-InSAR, the coseismic deformation field of the Luding MS6.8 earthquake is obtained. The source parameters of the earthquake are inverted based on the elastic half-space dislocation model, the slip distribution on the fault plane is determined by using the distributed slip model, and the static Coulomb stress change of the coseismic dislocation to the surrounding faults is calculated. It is shown that the maximum deformation cause by the earthquake is about 18 cm in the direction of radar line of sight (LOS). The coseismic dislocation is mainly left-lateral strike-slip at the depth of 6-24 km. The maximum slip is about 2.5 m at the depth of 16 km and the main seismogenic fault is the Xianshuihe fault. The static Coulomb stress triggering relationship shows that the Coulomb stress of the southern segment of Xianshuihe fault and the northern end of Anninghe fault increases significantly, and the strong earthquake trend is worthy of attention.
Key words: coseismic deformation; Sentinel-1A; InSAR; Luding earthquake; focal mechanism; Coulomb stress change