USE OF SEISMIC WAVEFORMS AND INSAR DATA FOR DETERMINATION OF THE SEISMOTECTONICS OF THE MAINLING MS6.9 EARTHQUAKE ON NOV.18, 2017

On November 18, 2017, a M_S6.9 earthquake struck Mainling County, Tibet, with a depth of 10km. The earthquake occurred at the eastern Himalaya syntaxis. The Namche Barwan moved northward relative to the Himalayan terrane and was subducted deeply beneath the Lhasa terrane, forming the eastern syntaxis after the collision of the Indian plate and Asian plates. Firstly, this paper uses the far and near field broadband seismic waveform for joint inversion (CAPJoint method)of the earthquake focal mechanism. Two groups of nodal planes are obtained after 1000 times Bootstrap test. The strike, dip and rake of the best solution are calculated to be 302°, 76° and 84° (the nodal plane Ⅰ)and 138°, 27° and 104° (the nodal plane Ⅱ), respectively. This event was captured by interferometric synthetic aperture radar (InSAR)measurements from the Sentinel-1A radar satellite, which provide the opportunity to determine the fault plane, as well as the co-seismic slip distribution, and assess the seismic hazards. The overall trend of the deformation field revealed by InSAR is consistent with the GPS displacement field released by the Gan Wei-Jun's team. Geodesy (InSAR and GPS)observation of the earthquake deformation field shows the northeastern side of the epicenter uplifting and the southwestern side sinking. According to geodetic measurements and the thrust characteristics of fault deformation field, we speculate that the nodal plane Ⅰ is the true rupture plane. Secondly, based on the focal mechanism, we use InSAR data as the constraint to invert for the fine slip distribution on the fault plane. Our best model suggests that the seismogenic fault is a NW-SE striking thrust fault with a high angle. Combined with the slip distribution and aftershocks, we suggest that the earthquake is a high-angle thrust event, which is caused by the NE-dipping thrust beneath the Namche Barwa syntaxis subducted deeply beneath the Lhasa terrane.     

基于解耦传播的波场分解方法在VTI介质弹性波逆时偏移中的应用

由于弹性波逆时偏移更符合实际情况,而且转换波的成像结果有更高的分辨率,因此弹性波逆时偏移越来越受到重视。然而,弹性波逆时偏移需要多波多分量数据,为了减少成像结果中的串扰假象,在逆时偏移过程中进行P波和S波分解就变得非常必要。结合基于向量的激发振幅成像条件,我们把基于解耦传播的波场分解方法应用到弹性波逆时偏移中,并对比了其在各向同性介质和横向各向同性（VTI）介质中的应用效果。结果说明,该方法可以在各向同性介质中完全分解P波和S波,并保留波场的向量信息。虽然在VTI介质中有较小的分解残余,但是该分解残余不会在逆时偏移结果中产生明显的串扰;因此,这种波场分解方法可以应用于各向同性介质和VTI介质弹性波逆时偏移。该方法是在时间空间域实现的,可以在波场传播过程中直接对P波和S波进行分离,应用方便,计算效率高。与利用Helmholtz分解的弹性波逆时偏移相比,该方法避免了在PS波成像结果中的极性反转问题。在复杂Hess VTI模型的逆时偏移结果中,高速岩体和断层的成像清晰,甚至是两个低速薄夹层也能较好成像;这说明该方法对复杂介质具有较好的适应性。PS波成像结果中的各向异性体成像清晰,说明各向异性介质弹性波逆时偏移可以对传统逆时偏移不能很好成像的构造进行成像。     

用连续GPS与远震体波联合反演2015年尼泊尔中部M_S8.1地震破裂过程

由于印度-欧亚板块碰撞,位于板块边界带的喜马拉雅地区大震频繁,但对其活动性的认识仍十分有限.2015年4月25日尼泊尔中东部地区时隔80年再次发生8级地震,为研究板缘地震提供了一次难得机遇.本文用西藏和尼泊尔的GPS连续观测数据和全球分布的远震地震波记录联合反演此次特大地震的破裂过程,结果显示此次地震发生在印度板块与青藏高原接触边界面——喜马拉雅主滑脱断层上.北倾11°、近东西(295°)走向的断层面破裂约100km长(博卡拉到加德满都),130km宽(从加德满都深入我国西藏吉隆县),破裂以逆冲滑动为主,平均幅度达到2.4m,释放的地震矩高达9.4×1020 N·m.反演结果还显示,震源体主要破裂分布深度范围为5~25km,应无地表破裂,属于一次盲地震.基于GPS资料推测的地壳现今运动速率及1833年地震的震源位置,我们推测地震在此次地震破裂区域复发的周期可能为150~200a,而极震区以南的深部滑脱断层仍保持闭锁,未来仍有导致灾害性大震的可能性.     

抛物Radon变换法近偏移距波场外推

本文给出了抛物Radon变换的基本原理,以及部分动校正后的CMP道集抛物线近似有效性的证明,基于带限正反最小平方抛物Radon变换的Levinson递推算法,对缺失的近偏移距地震波场进行叠前重建和外推.给出了抛物Radon变换法地震道重建外推的基本原理和叠前地震数据规则化的处理流程,另外对于Radon域均匀采样的情形,本文给出了均匀层状介质和Marmousi模型的近偏移距外推结果,计算结果验证了算法的稳定性和适用性.     

矿山微地震全波自动定位的层析成像方法(英文)

针对在矿山微地震监测过程中纵横波及其传播速度难以确定的难点,在球坐标系中以波动方程为基础,推导出了矿山微地震全波自动定位的层析成像的计算公式,并给出了参数选择准则,藉此可以获得微地震成像能量最大值及其对应参数。通过判断能量最大值在风险区内的位置,可以将风险区内外的微地震事件区分出来。对于能量最大值位于风险区域内的微地震事件,其震源位置和震动时间以及传播速度与能量最大值对应的坐标参数相同,能量最大值对应于震源函数相对强度;当能量最大值位于风险区域边界上时,微地震事件不在风险区域范围之内。理论模型试算和矿山实际微地震资料的应用表明:层析成像方法无需纵横波识别和传播速度假定,只要给定事故风险区域范围及其传播速度范围,便可完成微地震事件定位,是一种全波全自动定位方法,具有进一步研究价值。     

Stationary phase point extraction based on high-resolution radon transform and its application to multiple attenuation by wavefield extrapolation

Multiple prediction and subtraction techniques based on wavefield extrapolation are effective for suppressing multiple related to water layers.In the conventional wavefield extrapolation method, the multiples of the seismic data are predicted from the known total wave field by the Green function convoluted with each point of the bottom.However, only the energy near the stationary phase point has an effect on the summation result when the convolutional gathers are added.The research proposed a stationary phase point extraction method based on high-resolution radon transform.In the radon domain, the energy near the stationary phase point is directly added along the convolutional gathers curve, which is a valid solution to the problem of the unstable phase of the events of multiple.The Curvelet matching subtraction technique is used to remove the multiple, which improved the accuracy of the multiple predicted by the wavefield extrapolation and the artifacts appearing around the events of multiple are well eliminated.The validity and feasibility of the proposed method are verified by the theoretical and practical data example.     

基于交错网格Fourier伪谱微分矩阵算子的地震波场模拟GPU加速方案

作为高精度波形反演或逆时偏移的重要组成部分,地震波数值模拟对计算速度和效率提出了更高要求.GPU通用计算技术的产生及其内在数据并行性,为高效地震波数值模拟应用和研究得以有效开展奠定了基础.本文借助交错网格的Fourier伪谱微分矩阵算子和GPU上高效矩阵乘法,实现了复杂介质中地震波模拟的高效算法.数值试验表明,优化后的GPU计算相比CPU单核计算在大规模二维地震波场计算中获得至少100x以上的加速比.这对我们快速分析目标反射层在地震剖面中同相轴位置,制定优化采集方案具有重要意义.     

模拟地震波场的伪谱和高阶有限差分混合方法

伪谱法是一种高效、高精度计算非均匀介质地震波传播的数值方法,但是由于它的微分算子的全局性,使得该方法不适用于分散内存的并行计算.本文将有限差分算子的局部性和伪谱法算子的高效、高精度相结合,发展基于两种方法的伪谱/有限差分混合方法.该方法在一个空间坐标方向上利用交错网格高阶有限差分算子,在另外的空间坐标方向上利用交错网格伪谱法算子,既保留了后者的高效、高精度优势,又便于在PC集群上实现并行计算.对二维模型的计算显示,混合方法能有效处理介质不连续面,在保证伪谱法计算精度的情况下,提供了一种并行计算的可能途径.     

基于源波场重建与波场分解的逆时偏移

逆时偏移以其高精度、适于强变速地质体和不受反射层倾角限制的特点,已经成为处理盐丘、高陡倾角反射层和推覆构造等复杂地质体的标准成像手段。然而,其成像条件需要具有同一时间上的源波场与检波点波场,这在大规模数据上会导致存储的不足与I /O 时间过长。这种成像条件也会在波场传播路径上产生假象。鉴于此,提出了一套可以有效解决这些问题的逆时偏移方案。对于存储与I / O 问题,用源波场重建技术来使源波场与检波点波场在时间上同步,即在源波场正演过程中,仅需存储成像区域边界的波场和最后两层时间的波场作为源波场重建的边界条件和初始条件,从而大幅度节约存储与I /O。对于假象问题,利用上下左右行波分解成像条件与Laplace 滤波技术结合,以达到去除假象的目的。经测试,该方案取得较好的成像效果。其意义在于改进了逆时偏移去假象的方案,使其可以有效地对复杂地质情况进行高质量的成像,并通过源波场重建技术来减少实现逆时偏移所需要的存储与I /O 时间。     

利用远震接收函数反演乌江彭水电站地震台下方地壳厚度

乌江彭水电站目前是重庆市最大的水电站,电站总投资121亿元,总装机容量为175万kw.坝高116.5 m,水库正常蓄水水位293 m,总库容13.35亿m3,为河道型高坝水库.已有研究结果表明,彭水电站坝区是稳定的,但库区及外围地区仍存在5级左右的地震活动背景,在水库蓄水后有可能诱发中,强地震.本文利用彭水台网于200...