2010年高雄地震震源参数的近远震波形联合反演

本文改进了传统基于近震波形数据的点源震源参数反演的Cut And Paste(CAP)方法,实现了近震Pnl波、面波和远震P波、SH波的联合反演的CAPjoint算法.对2010年3月高雄地震,分别进行单独反演以及联合反演,获得各自的震源机制解及深度,其中联合反演所得的最佳双力偶机制解参数为,节面1:走向317°,倾角36°,滑移角52°,节面2:走向181°,倾角62°,滑移角114°,深度为21 km.并对不同震中距波形对本次地震以及几种典型机制解断层几何参数的敏感性进行测试.为验证联合反演方法的可靠性,本文采用重抽样思想发展而来的Bootstrap方法,对近震数据的子集及其与远震数据的联合反演所得的参数进行统计,验证了在稀疏近台条件下联合反演中添加远震数据对地震震源参数约束的作用.     

用转换函数方法研究喜马拉雅地区速度结构

利用流动数字地震台网提供的三分量地震波形记录,应用转换函数及快速模拟退火算法对喜马拉雅山脉地区46个地震站下的地壳横波速度结构进行了反演,为进一步揭示青藏高原喜马拉雅山脉地区的动力学演化过程提供了可靠的地球物理证据.根据本文结果可清晰看到,喜马拉雅山脉地区作为当今地壳活动最活跃的地区,物质交换非常活跃,地下结构远远未达到平衡,地壳速度有很大差异,在板块边界处莫霍界面速度间断不是非常明显,自喜马拉雅南坡向高原腹地,地壳厚度大致从55 km增加到80 km;沿经度方向,莫霍面也有一定的起伏.通过研究得到另外一个证据是,在喜马拉雅的主中央逆冲断裂,由大陆碰撞产生的主要构造,其深度可能要大于80 km.     

智利M_W8.3地震与M_W8.8地震的震源过程及其引起的库仑应力特征

2015年9月17日6时54分32秒(北京时间)智利中部伊拉佩尔附近(震中31.57°S,71.67°W)发生了一次M_w8.3大地震,在此次地震震中以南约500 km处的马乌莱地区曾于2010年2月27日14时34分11秒发生过一次M_w8.8强震(震中36.12°S,72.90°W),两次地震余震分布区之间有约75 km的地震空区.本文利用远场体波与面波波形,基于有限断层模型,反演了这两次地震的震源破裂过程.结果显示这两次地震均为逆冲型大地震,2015年伊拉佩尔M_w8.3地震的平均滑动角度为107°,平均滑动量为2.43 m,平均破裂速度为1.82 km·s~(-1),标量地震矩为3.28×10~(21)Nm,95%的标量地震矩在104 s内得到了释放.最大滑动量约8 m,位于沿走向75 km,深度8 km处.2010年马乌莱M_w8.8地震的平均滑动角度为109°,平均滑动量为4.95 m,平均破裂速度1.90 km·s~(-1),标量地震矩为1.86×10~(22)Nm,95%的标量地震矩在121 s内得到了释放.最大滑动量约12.5 m,位于沿走向100 km,深度21 km处.2015年伊拉佩尔M_w8.3地震浅部更大的滑动量应该是其引起了较大海啸的一个原因.基于破裂滑动分布,我们计算了这两次地震引起的周边俯冲带上静态库仑应力变化,结果显示两次地震均显著增加了周边俯冲带上的库仑应力,2010年马乌莱地震使得2015.年伊拉佩尔地震震源区附近的库仑应力增加了(0.01~0.15)×10~5Pa,从应力积累的角度看,2010年马乌莱地震有利于2015年伊拉佩尔地震的发生,对后者的发生起到了促进作用.     

中国台湾以东地区地震矩张量研究及其构造意义

本文利用长周期Ｐ波波形反演方法,推断了１９８６年１０月至１９９１年期间发生在中国台湾岛东部及其附近海域中８次中强地震的矩张量．矩张量解揭示了发生在这一区域边界处的地震震源性质的规律性变化：沿西部边界的主旋错动,沿北部边界的右旋错动,西北角花莲地区的逆冲断错及显著的体积收缩分量．这一结果表明,台湾地区是菲律宾海板块与欧洲大陆板块边界中的一特殊地段,海洋板块的俯冲在这里受阻．随着菲律宾海板块向西北方向的相对运动,在花莲地区发生了强烈碰撞．     

龙门山断裂带周缘重复地震识别及其在台网定位评价中的应用

利用国家测震台网数据备份中心时间跨度超过7年的连续波形资料, 设定至少三个台站记录的垂直分量波形互相关系数大于0.8的地震对为一组重复地震, 通过波形互相关分析, 识别出龙门山断裂带周缘波形相关意义上的“重复地震”2790次, 构成2907组重复地震。 沿用Schaff和Richards的研究结果, 假定“重复地震”间距很小(小于1 km), 地震目录记录的重复地震对位置之差主要为地震台网的定位误差所致, 基于此误差给出了龙门山断裂带周缘地震台网的定位精度估计: 台网的水平定位精度较高, 水平定位误差约为2.8 km(2倍标准差), 且西南段台网的水平定位精度优于中北段; 垂直定位精度较差, 垂直定位误差约为10 km(2倍标准差), 现有地震定位方法对震源深度的测定有待改善。     

庐枞矿集区大地电磁测深强噪声的影响规律

天然大地电磁场信号微弱,极化方向随机,极易受电磁噪声污染.张量阻抗分析、远参考技术、Robust估计等对随机噪声和不相关噪声有比较好的压制效果,但对强的相关噪声目前还没有有效的压制方法.庐枞矿集区人烟稠密、工业发达,是我国著名的铁、铜、硫等矿产基地之一,区内强烈的工业、通讯、矿山、民用等电磁干扰严重污染了大地电磁测深数据.本文首先根据实测的电磁场时间域波形和卡尼亚电阻率测深曲线形态,挑选出基本未受噪声污染的测点(Y1650).然后利用数学形态学从受严重污染的电磁场时间序列中提取出类方波、三角波、阶跃、脉冲和充放电5种典型噪声的波形,并以不同的方式将这些噪声波形与Y1650的电磁场波形叠加,对比分析加噪后Y1650点电阻率和相位测深曲线的变化,进而归纳出典型强噪声对庐枞大地电磁测深资料的影响规律,为进一步压制强噪声和资料处理提供依据.     

中国及邻近地区地壳结构

利用中国及邻近地区数字地震台网提供的三分量地震波形记录,应用转换函数及快速模拟退火算法对中国及邻区61个地震站下的地壳横波速度结构进行了反演. 结合已发表的人工地震测深资料获得了该地区地壳厚度分布,为进一步揭示中国大陆的动力学演化过程、建立大陆动力学理论提供了可靠的地球物理证据. 结果表明中国及邻区地壳厚度分布的特点为从东到西逐渐增厚,可分为：地壳厚度缓变地区,如蒙古高原,华北,华中,华南,青藏高原内部;地壳厚度递变带,如沿大兴安岭—太行山—秦岭—大巴山—云贵高原一带;地壳厚度陡变带,主要为青藏高原边缘地区. 除青藏高原外,其他地壳厚度缓变区速度结构较简单,壳幔边界明显. 青藏高原及其南缘,台湾—菲律宾一带,台站下方速度结构复杂,反映了板块边界处构造活动、物质交换活跃,表明这些地区还未达到均衡.      

动态波形匹配预测火山岩地层的横向变化

讨论一种根据3D地震数据计算地层横向变化的方法动态波形匹配方法,该方：该方法利用动态规划算法对相邻地震道地震波形进行动态波形匹配并计算波形匹配费用．把波形匹配费用作为相邻地震道地震波形之间差异大小程度的度量．对大庆油田一块800km^2的3D地震数据的计算表明,动态波形匹配费用可以较好地度量相邻地震道地震波形之间的差异,表示地下地质体的横向变化,如断层、裂缝发育带和火山岩岩体等．     

Misfit functionals in Laplace‐Fourier domain waveform inversion,with application to wide‐angle ocean bottom seismograph data

In seismic waveform inversion, non‐linearity and non‐uniqueness require appropriate strategies. We formulate four types of L2 normed misfit functionals for Laplace‐Fourier domain waveform inversion: i) subtraction of complex‐valued observed data from complex‐valued predicted data (the ‘conventional phase‐amplitude’ residual), ii) a ‘conventional phase‐only’ residual in which amplitude variations are normalized, iii) a ‘logarithmic phase‐amplitude’ residual and finally iv) a ‘logarithmic phase‐only’ residual in which the only imaginary part of the logarithmic residual is used. We evaluate these misfit functionals by using a wide‐angle field Ocean Bottom Seismograph (OBS) data set with a maximum offset of 55 km. The conventional phase‐amplitude approach is restricted in illumination and delineates only shallow velocity structures. In contrast, the other three misfit functionals retrieve detailed velocity structures with clear lithological boundaries down to the deeper part of the model. We also test the performance of additional phase‐amplitude inversions starting from the logarithmic phase‐only inversion result. The resulting velocity updates are prominent only in the high‐wavenumber components, sharpening the lithological boundaries. We argue that the discrepancies in the behaviours of the misfit functionals are primarily caused by the sensitivities of the model gradient to strong amplitude variations in the data. As the observed data amplitudes are dominated by the near‐offset traces, the conventional phase‐amplitude inversion primarily updates the shallow structures as a result. In contrast, the other three misfit functionals eliminate the strong dependence on amplitude variation naturally and enhance the depth of illumination. We further suggest that the phase‐only inversions are sufficient to obtain robust and reliable velocity structures and the amplitude information is of secondary importance in constraining subsurface velocity models.     

The April–June 2007 Trichonis Lake earthquake swarm (W. Greece): New implications toward the causative fault zone

On 10 April 2007, three moderate earthquakes with M_w = 4.9–5.1 occurred in the vicinity of Trichonis Lake (W. Greece). A local network composed of 12 three-component digital seismographs was installed in the epicentral area and recorded more than 1600 events. The double-difference algorithm HYPODD, incorporating both catalog and waveform cross-correlation differential travel-time data, was applied for the successful relocation of 1490 earthquakes. The latter led to the distinction of a main NW-SE trending and NE-dipping zone, as well as of three neighboring faults; a conjugate NW-SE striking and SW-dipping marginal fault mapped along the northeastern flanks of the lake; a E-W trending and south-dipping low-angle normal fault, possibly related to the major Agrinio Fault Zone (AFZ), parallel to the northern bank of the lake; a NE-SW striking and NW-dipping normal fault, likely related to a segment of the active Evinos fault, located south of the lake. Calculation of the Coulomb stress induced by the combination of the 1975 M_w = 6.0 event and the three largest events of 10 April 2007 on the inferred structures, reveals that most of the seismicity lies within the “stress-loaded” region, except for the westernmost activity, which probably belongs to the deep part of the AFZ. A total of 178 reliable focal mechanisms were determined by regional and local body-wave modeling (5 largest events) and P-wave first motion polarity data. The types of the obtained focal mechanisms are predominantly normal and strike-slip, however, numerous earthquakes were found to exhibit reverse faulting. Inversion of focal mechanism data showed that the prevailing principal horizontal component σ₃ is quite homogeneous throughout the activated area with a roughly NW-SE trend, parallel to the strike of the Hellenides. On the contrary, the compressional field σ₁ appears in two patterns: NE-SW trending onshore and NW-SE trending beneath the lake. This apparent rotation of σ₁ by 90° reveals a complex system enclosed by the suggested NW-SE trending antithetic faults in depths between 7 and 9 km. The calculated stress ratios beneath the lake imply that vertical forces are close to the overburden pressure. The overall inferred stress pattern is rather linked to topographic variations, locally imposing increase or decrease of the vertical forces. The presence of the water in the lake possibly plays an additional important role, penetrating through the bedrock, reducing the friction coefficient, while the pore pressure and, consequently, the effective stress increase. Thus, shearing along mature fractures is enhanced, likely yielding the observed diversity.