巴姆地震变形场和应力场：Ⅱ．用FEPG有限元方法求解

本文利用 FEPG 有限元程序对2003年12月的巴姆地震变形场和应力场进行了模拟计算。2003年12月26日,巴姆地区发生6.6级地震。Nakamura、Suzuki 等人利用余震资料和地震破坏程度的分布图分析了余震震中的统计分布形式,得到了此次地震的发震断层是巴姆断层以西大约5km 的 Arg-e-Bam 隐伏断层的结论。凌勇等人从差分干涉雷达(D-InSAR)出发,得到了此次地震的同震干涉条纹和雷达视线方向的位移场,并从雷达干涉的相干图上确定了发震断层的痕迹。他们还利用 Okada 程序模拟计算了此次地震的位移场,模拟计算的结果与 D-InSAR 测量的位移场比较符合。Okada 程序从断层的位错为出发点进行计算,从运动学角度对此次地震变形场和应力场进行研究;本文从动力学的角度出发,研究远场区域边界的应力状态以及断层的物性参数的变化对巴姆地震形变场的影响,将各种计算结果进行对比的同时也和前人的结果相比较。本文计算中采用了凌勇等人的断层参数,即 Arg-e-Bam 隐伏断层为近南北走向,向东倾斜,近直立的倾角,右旋走滑,其(Strike,dip,slip)为(175°,80°,173°)。断层未延伸到地表,自地面以下1km 处往下延伸16km 的宽度,断层长度为20km。对断层的杨氏模量、泊松比、断层附近区域的边界受力大小、两个主应力的比例以及主应力的方向的改变对计算结果的影响做了讨论。最后通过上述几组计算结果建立了一个与前人结果较为符合的巴姆地震参数模型。通过计算,一方面揭示了应力状态和物性参数的影响,另一方面也说明 Arg-e-Bam 隐伏断层是此次地震的发震断层这一结论是可信的。     

基于差分干涉雷达的汶川地震同震形变特点

地震同震形变场是认识震源机制、确定发震断裂、分析发震断裂活动方式、评估震害损失及揭示未来发震趋势的重要依据。但在龙门山地区复杂的自然环境中, 现有的地震同震形变场获取方法(GPS测量、构造形迹分析、震源机制解结合数值模拟计算)难以快速、全面的捕捉到这种信息, 这突显出差分干涉雷达技术(D-InSAR, Differential Interferometry SAR)在同震形变场研究中的优势。本文介绍了其观测原理, 在地震同震形变监测中的研究现状, 并重点分析了已获取的D-InSAR汶川地震同震形变场的观测结果。基于差分干涉雷达观测及相关资料, 对龙门山断裂在本次地震中活动性进行了初步分析, 结果表明干涉雷达形变图像与地表破裂范围、逆冲角度的变化、上下盘升降关系及大型余震的展布存在很强的相关性。      

玉树地震地表变形InSAR观测及初步分析

采用玉树M_S7.1级地震前后两期PALSAR雷达数据(震前2010年1月15日, 震后4月17日)进行了“两轨+DEM”的InSAR处理, 获得了高质量的差分干涉雷达条纹图像和同震变形场。参考该区的基本构造格局, 根据干涉图像的变形范围、变形量和变形梯度可以初步判断:(1)玉树地震诱发了总体上NWW走向, 全长约70km地表陡变带, 陡变带南段位错及陡变梯度较大, 会在地表产生地表破裂;而西北部4段位错及陡变梯度较小, 不易在地表诱发破裂, 但可能在地下一定层位产生了隐伏破裂带;(2)陡变带两侧的雷达视线向运动方向预示发震断裂以左旋走滑运动为主;(3)宏观震中位于玉树县城西北约16km的地表陡变带上。D-InSAR解译结果与中国地震台网中心震源机制解、野外发震断裂调查结果及地貌特征吻合较好, 证明了干涉雷达解译成果的可靠性, 可以为准确定位玉树地震发震断裂地表行迹和快速评定震害损失提供有力的技术支持。      

利用雷达差分干涉测量技术获取汶川MS 8.0地震形变场

利用ENVISAT/ASAR和ALOS/PALSAR数据,结合两路差分干涉测量技术获得了汶川地震的同震形变信息和震前、震后地面在卫星视线(light of sight,简称LOS)方向上的形变特征。结果表明,汶川地震引起的地面形变范围广,程度大,地震形变沿断层向北北东向扩展,所形成的差分干涉条纹明显。地面的形变特征对于推断断层的性质,研究地震形变和地震发育特征具有重要的参考价值。另外,将同震雷达差分干涉测量的结果与利用USGS发布的汶川地震有限元断层模型(finite fault model)反演的LOS方向的形变进行了对比验证,发现二者在断裂带的上盘具有较好的一致性,但在下盘却有较大的误差。通过不同传感器干涉结果的对比发现,L波段的雷达干涉结果更能够反映汶川大地震在龙门山断裂带附近的地面形变信息。     

汶川5.12地震发震机理的数值模拟

将围岩看成弹性体, 断层看成具有应变软化的弹塑性体, 断层和围岩组成统一的地质介质系统。在给定的地应力、孔隙压力及边界位移的作用下, 应力逐渐积累, 当达到断层强度时, 断层进入塑性状态, 应变随之软化, 应力突然下降, 能量突然释放, 形成地震。根据这个原理, 利用有限差分FLAC软件, 对汶川5.12地震的机理进行了模拟。模拟结果证实:汶川5.12地震是在青藏高原侧向挤压, 四川盆地阻挡, 使龙门山断裂带受到缓慢增加的挤压应力场作用下形成的。计算结果给出了应力降、能量释放量、断层错动量、地震同震位移、震前位移、地震复发周期等重要参数, 结果与野外调查及其它资料给出的数据具有较好的一致性。      

基于D—InSAR的九寨沟地震同震形变场反演分析

利用覆盖九寨沟地区的RadarSat—2数据与Sentinel—1A数据,采用精轨数据进行定轨,消除轨道误差,并结合合成孔径差分(D-InSAR)方法中的双轨差分技术,获取2017年8月8日Mw7. 0级地震的同震形变场。结果表明,视线方向(LOS)最大沉降量约为20 cm,隆起量达9 cm。基于弹性半空间形变模型反演该地震的断层滑动分布,得出该地震断层滑动以左旋走滑为主,走向为330°,倾角为32°,滑动角为-170°,同震滑动分布主要集中在4～12 km深度处,最大滑动量位于9 km处,约为6. 14 m,平均滑动量为0. 57 m。反演获得的地震标量矩为4. 06E+18N·m,震级Mw约为6. 4,深度为19. 5 km。     

青海玉树M_S7.1地震发震过程的数值模拟

根据玉树地区的地应力场、速度场和断层展布,对青海玉树2010年4月14日MS7.1级地震发震机理进行了数值模拟。将围岩看成弹性体,断层看成具有应变软化的弹塑性体,断层和围岩组成统一的地质介质系统。在地应力、孔隙压力及边界位移的作用下,应力逐渐积累,当达到断层摩擦破坏强度时,断层产生应变软化,断层突然滑动,能量突然释放,应力突然下降,形成地震。模拟结果表明:玉树7.1级地震是在印度板块向北推挤,青藏高原向东南侧向挤压,在玉树地区形成主压应力为北东80°方向的水平应力场,使甘孜-玉树断裂带产生左旋走滑错动形成的。计算结果给出了应力降、能量释放量、断层走滑错动量、地震复发周期、应力积累速度等重要参数,模拟结果与野外调查资料具有较好的一致性。     

2022年青海门源MS6.9地震发震构造和地表破裂初步调查

北京时间2022年1月8日01时45分，青海省海北州门源县发生强烈地震（图1），造成数人受伤，房屋倒塌，部分道路、桥梁、隧道等基础设施被破坏或受损。中国地震台网（CENC）测定该地震的震级为MS 6. 9，震中位于37. 77°N，101. 26°E，震源深度为10 km(https://www.cenc.ac.cn/cenc/dzxx/396391/index.html)。利用欧洲航空局哨兵2号雷达卫星的震前、震后SAR数据进行差分干涉处理，得到同震形变场分布图。限定此次地震以左旋走滑运动为主，断层走向NWW，断层面近直立；主体破裂深度在10 km以上并到达地表，形成长度＞35 km的地表变形带，最大滑动量约2 m。2022年门源MS 6. 9地震发生在青藏高原中北部的祁连- 柴达木次级地块的北部（图1）、托莱山断裂带和冷龙岭断裂带的交会部位，是继1986年和2016年两次门源MS 6. 4地震之后在冷龙岭断裂带上发生的震级最高、地表破裂最长的地震事件。     

日本MW9.0级地震同震位移场和应力场的有限元模拟

根据同震位移GPS观测数据, 利用有限元法反演了2011年3月11日本MW9.0级地震的断层滑移模式。在此基础上, 计算了日本MW9.0级地震引起的同震位移场和应力场, 给出了位移和应力的分布, 分析了他们的变化规律并与实测结果进行了对比。计算结果表明: 日本MW9.0级地震的静态断层滑移量最大可达25 m。地震引起断层上盘向东位移, 最大位移在震中附近, 可达24.25 m, 日本东北地区向东位移最大可达6 m。震后地表隆起, 隆起幅度可达5.6 m, 隆起的最高点也在震中附近。日本东北地区东海岸附近有一下沉带, 下沉量可达0.8 m。同震地表位移的计算值与GPS测量结果基本一致。地震引起应力变化, 导致震后应力下降。应力变化是不均匀的, 在震中附近约为9.9 MPa, 在深处可达32 MPa, 在日本东北地区地表应力变化小于4.4 MPa。地震引起的应力变化主要是水平应力, 垂直应力基本不变。     

青藏高原东缘1933年叠溪Ms7.5级地震发震构造再研究

青藏高原东缘1933年叠溪75级地震的发震构造至今仍然难以琢磨,前人或将其归因于NW向松坪沟断裂的左旋走滑活动、或南北向岷江断裂左旋走滑活动,但地表同震破裂、地震地质、地震等烈度图等调查和研究结果都不支持这种走滑型断层的地震成因。本文基于叠溪地震区构造地貌和湖相地层断层调查,结合古地震和历史地震的研究结果,提出了与2013年四川芦山Ms 70级地震类似的发震构造模型,即隐伏断坡型逆冲断层发震构造模型,认为在叠溪震区10~15km深部隐伏一条西倾的逆冲断坡,其向东逆冲作用导致了叠溪地区频繁的地震活动。这个发震模型有待深部地球物理测深资料和地表大地测量资料的验证。     

Retrieving three-dimensional coseismic displacements of the 2008 Gaize,Tibet earthquake from multi-path interferometric phase analysis

In this paper, synthetic aperture radar (SAR) data from ENVISAT ASAR ascending, descending and ALOS PALSAR ascending orbits are collected to investigate the coseismic displacements of the Mw 6.4 earthquake occurred in Gaize, Tibet on January 9, 2008 and the Mw 5.9 aftershock on January 16, 2008. Two interferometric phase analysis techniques, i.e., D-InSAR and multi-aperture InSAR, are employed to process the SAR data, with which the displacement measurements along three different line-of-sight (LOS) and three different azimuth directions are retrieved, respectively. Complete three-dimensional (3-D) coseismic displacement fields caused by the earthquake are then resolved by integrating the obtained LOS and azimuth displacement measurements with a weighted least squares adjustment, whose distributions are conformed to the two north-northeast trending northwest-dipping normal faults detected in previous studies. Ground subsidence and uplift are observed in the hanging wall and footwall of the main fault, respectively, and the subsidence reaches its maximum in the hanging wall of the second fault as a superimposed result of the Gaize earthquake and its aftershock. Anti-symmetric horizontal movements are also detected during the seismic events, which move inward in the focal region, but outward at the marginal. The left-lateral motions near the main fault indicate a small striking slip component caused by the Gaize earthquake. Finally, we discuss the potential of applying the derived spatially continuous 3-D displacement fields to determine the high-resolution 3-D strain fields of the Gaize earthquake, which provide important knowledge for assessing the source mechanism.     

A Characteristic Analysis of the Dynamic Evolution of Preseismic- Coseismic-Postseismic Interferometric Deformation Fields Associated with the M 7.9 Earthquake of Mani, Tibet in 1997

By using the D-InSAR technique,we have acquired the temporal-spatial evolution images of preseismic-cosesimci-postseismic interferometric deformation fields associated with the M 7.9 earthquake of Mani,Tibet on 8 November 1997.The analysis of these images reveals the relationships between the temporal-spatial evolution features of the interferometric deformation fields and locking, rupturing,and elastic restoring of the source rupture plane,which represent the processes of strain accumulation,strain release,and postseismic restoration.The result shows that 10 months prior to the Mani event,a left-lateral shear trend appeared in the seismic area,which was in accordance with the earthquake fault in nature.The quantity of local deformation on the north wall was slightly larger than that on the south wall,and the deformation distribution area of the north wall was relatively large.With the event impending,the deformation of the south wall varied increasingly,and the deformation center shifted eastward.Two and half monthd before the event,the west side of the fault was still locked while the east side began to slide,implying that the whole fault would rupture at any moment.These features can be regarded as short-term precursors to this earthquake.Within the period from 16 April 1996 to two and half months before the earthquake,the most remarkable deformation zones appeared in the north and south walls,which were parallel to and about 40 km apart from the fault,with accumulated local displacements of 344 mm and 251 mm on the north and south walls,respectively.The south wall was the active one with larger displacements.Five months after the earthquake,the distribution feature of interferometric fringes was just opposite to that prior to the event,expressing evident right-lateral shear.The recovered displacements are～179 mm on the north wall and～79 mm on the south wall,close to the east side of the fault.However,in the area of the south wall far from the fault there still existed a trend of sinistral motion.The deformation of the north wall was small but recovered fast in a larger area,while the active south wall began to recover from the east section of the fault toward the WSW.     

Analysis of the Wenchuan Ms8.0 Earthquake’s Co-seismic Stress and Displacement Change by Using the Finite Element Method

The variation of in situ stress before and after earthquakes is an issue studied by geologists. In this paper, on the basis of the fault slip dislocation model of Wenchuan Ms8.0 earthquake, the changes of co-seismic displacement and the distribution functions of stress tensor around the Longmen Shan fault zone are calculated. The results show that the co-seismic maximum surface displacement is 4.9 m in the horizontal direction and 6.5 m in the vertical direction, which is almost consistent with the on-site survey and GPS observations. The co-seismic maximum horizontal stress in the hanging wall and footwall decreased sharply as the distance from the Longmen Shan fault zone increased. However, the vertical stress and minimum horizontal stress increased in the footwall and in some areas of the hanging wall. The study of the co-seismic displacement and stress was mainly focused on the long and narrow region along the Longmen Shan fault zone, which coincides with the distribution of the earthquake aftershocks. Therefore, the co-seismic stress only affects the aftershocks, and does not affect distant faults and seismic activities. The results are almost consistent with in situ stress measurements at the two sites before and after Wenchuan Ms8.0 earthquake. Along the fault plane, the co-seismic shear stress in the dip direction is larger than that in the strike direction, which indicates that the faulting mechanism of the Longmen Shan fault zone is a dominant thrust with minor strike-slipping. The results can be used as a reference value for future studies of earthquake mechanisms.     

Source characteristics of the Yutian earthquake in 2008 from inversion of the co-seismic deformation field mapped by InSAR

On 21 March 2008, an Ms7.3 earthquake occurred at Yutian County, Xinjiang Uygur Autonomous Region, which is in the same year as 2008 Mw 7.9 Wenchuan earthquake. These two earthquakes both took place in the Bayar Har block, while Yutian earthquake is located in the west edge and Wenchuan earthquake is in the east. The research on source characteristics of Yutian earthquake can serve to better understand Wenchuan earthquake mechanism. We attempt to reveal the features of the causative fault of Yutian shock and its co-seismic deformation field by a sensitivity-based iterative fitting (SBIF) method. Our work is based on analysis and interpretation to high-resolution satellite (Quickbird) images as well as D-InSAR data from the satellite Envisat ASAR, in conjunction with the analysis of seismicity, focal mechanism solutions and active tectonics in this region. The result shows that the 22 km long, nearly NS trending surface rupture zone by this event lies on a range-front alluvial platform in the Qira County. It is characterized by distinct linear traces and a simple structure with 1–3 m-wide individual seams and maximum 6.5 m width of a collapse fracture. Along the rupture zone are seen many secondary fractures and fault-bounded blocks by collapse, exhibiting remarkable extension. The co-seismic deformation affected a big range 100 km × 40 km. D-InSAR analysis indicates that the interferometric deformation field is dominated by extensional faulting with a small strike-slip component. Along the causative fault, the western wall fell down and the eastern wall, that is the active unit, rose up, both with westerly vergence. The maximum subsidence displacement is ～2.6 m in the LOS, and the maximum uplift is 1.2 m. The maximum relative vertical dislocation reaches 4.1 m, which is 10 km distant from the starting rupture point to south. The 42 km-long seismogenic fault in the subsurface extends in NS direction as an arc, and it dipping angle changes from 70° near the surface to 52° at depth ～10 km. The slip on the fault plane is concentrated in the depth range 0–8 km, forming a belt of length 30 km along strike on the fault plane. There are three areas of concentrating slip, in which the largest slip is 10.5 m located at the area 10 km distant from the initial point of the rupture.     

九寨沟7.0级地震的地震断裂及震源破裂的构造动力学机理研究

本文针对九寨沟7.0级地震的发震构造及震源破裂的构造动力学问题,从地质调查、测试与震源构造动力学理论分析相结合的角度阐明:此次地震发生在处于高法向应力左旋剪切受力状态的岷山隆起带北端与西秦岭地槽褶皱带之南缘文县-玛沁断裂相交汇部位;沿"九寨天堂"-震中-五花海以及上四寨村-中查-比芒一线发育的两条具有构造破裂特征的NW向地裂缝带,是"隐伏状"地震断裂错动造成的地表同震破裂,前者是本次地震的控震构造,后者是NW向断层的次级同震复活及扩展;此次地震的震源断层错动过程受NW331°走向的陡倾角（∠87°）断层节面控制的走滑型剪切破裂,断裂两端的剪切破裂表现出持续扩展趋势;震源构造动力学过程实质上是地壳深部沿NW方向左旋剪断"凸出体";地震断裂带两端破裂扩展的持续发展,将导致地壳"凸出体"逐步剪断贯通。     

Relationships between InSAR Seismic Deformation and Fault Motion Sense, Fault Strike, and Ascending/Descending Modes

Based on the working principle of satellite radars, the earthquake deformation field measured by interferometric synthetic aperture (InSAR) is the projection of ground displacement associated with the seismogenic fault in the line of sight (LOS) of the satellite. However, LOS projections are complex, and are not only related to the ascending/descending modes and incidence angles of SAR data, but also related to the strike and motion senses of the fault. Even for the same earthquake, the LOS deformation derived from different ascending/descending data can be almost identical in one case, but quite different in another case, which makes the interpretation of InSAR seismic deformation and its comparison with field observations difficult. In this study, we undertook a quantitative analysis of the relationships between LOS observation sensitivity of InSAR and fault strike, fault motion sense, and ascending/descending modes, as well as 3D deformation fields. We studied the features and differences of the LOS deformation fields in different types of earthquakes using ascending/descending modes, with a particularly detailed analysis of the relations for a strike-slip type of earthquake. We also summarized the characteristics of LOS deformation fields of faults with different strikes and optimal observational data modes. Taking the strike-slip Yushu earthquake and the normal Gaize event as examples, we used SAR data of the ascending/descending modes to verify the results of quantitative calculations. These analyses will not only provide a more reasonable interpretation of InSAR seismic deformation fields and but also help understand the differences of seismic deformation fields revealed by data with different observational modes, therefore promoting the application of InSAR technology in seismology.     

玉树地震地表破裂特征及其破裂方式

北京时间2010年4月14日玉树Ms7.1级地震产生了长约50km的地表破裂带,地表破裂系统非常典型。本文通过对果庆益荣松多村附近河滩与结古镇西南河道两岸2个地表破裂点的野外观测和地表破裂的力学分析,探讨玉树地震地表破裂特征和方式。在果庆益荣松多村附近河滩地表破裂与地震陡坎呈间隔式左行右阶排列,具转换压缩性质;在结古镇西南河道北岸主地表破裂带中地表破裂非常发育,出露一处古断层面。地表破裂的野外调查与应力分析结果表明,玉树地震同震地表破裂带总体走向为125°,属于典型的左行走滑破裂带,最大左行走滑位移量1.75m,主压应力为近东西向,这一结果与美国地质调查局和中国地震局地震地球物理研究所公布的结果基本一致。     

郯庐断裂带安徽段现代地震活动及应力场特征

分析了郯庐断裂带安徽段地震活动特征,利用近震直达波最大振幅比方法反演了99次中-小地震震源机制,并作了统计和聚类分析。现代地震活动整体上继承了历史地震分布格局,未显示增强趋势;断裂带及邻区应力场P轴和T轴优势方位大致呈近东西向和近南北向分布,现今承受的作用力以近水平或斜向为主,但也存在部分其它方向的应力场和近垂直方向的作用力。震源断层的破裂类型,带内以逆冲(或正断),或近逆冲(或正断)型为主。带附近较明显地呈现出走滑型或近走滑型破裂的优势;北、南亚段应力场有明显差异,即北压南张,其分界大约在北纬32°附近,北西西走向的桥头集—东关断裂是两者间的构造变换带。     

关于海城地震震源模式的讨论

地震震源模式是对地震形成条件、孕育和暴发过程的认识,它是与地震预报密切相关的一个重要理论课题。讨论一个地震的震源模式应该解决震源应力场、发震构造和震源物理、化学过程等问题,解决这些问题的根本依据来源于震前、震时和震后所暴露出的各种实际材料。地震地质特征是认识这些问题的背景,模拟实验则为我们认识这些问题提供对比研究的线索。