高频GNSS形变波的震相识别：模拟实验与实例分析

高频GNSS震时形变波震相及识别是GNSS地震学的重要内容.在实时数据处理基础上,本文利用振动台的高频GNSS观测实验,并结合近期部分大震的高频GNSS形变波震相特征进行研究.数据处理结果表明,实时处理与事后处理的精度在同一量级,且与采样率无关.通过与同址观测强震仪和地震计记录的对比及特定震相的频谱分析,发现高频GNSS可完整记录P波、S波、Love波及Rayleigh波震相,影响震相记录的主要因素是GNSS定位精度与震级,而仅当震中距很小时,采样率将产生一定影响.研究结果得出：基于地震波传播规律,利用高频GNSS台阵记录的形变波空间分段特征,结合震相运动轨迹及其他地震波记录,可实现实时高频GNSS形变波的震相识别.     

利用高频GPS进行地表同震位移监测及震相识别

结合日本3.11大地震,利用近震区和远震区若干测站的高频GPS观测数据,采用动态精密单点定位技术分析了强震时复杂的地表运动时序,获得了地表同震位移.在此基础上提出利用S变换对高频GPS地震波进行震相识别,有效拾取了P波初至,反演了P波传播速度并与地震仪数据、USGS提供的参考值进行了对比和分析,为后续地震应急、地震预警提供了重要的、科学的决策依据.     

基于小波包和峰度赤池信息量准则的P波震相自动识别方法

基于小波包变换和峰度赤池信息量准则(AIC), 提出了一种新的自动识别P波震相的综合方法, 即小波包-峰度AIC方法． 首先对由加权长短时窗平均比(STA/LTA)法粗略确定的P波到时前后3 s的记录进行小波包三尺度的分解与重构, 分别计算每个尺度重构信号的峰度AIC曲线并将其叠加, 叠加曲线的最小值则为P波震相到时; 然后对原始地震记录进行有限冲激响应自适应滤波以提高信噪比和识别精度; 最后将小波包-峰度AIC方法应用到合成理论地震图及实际地震记录的P波初至自动识别中. 结果表明: 初至清晰度对识别精度的影响比信噪比对其影响更大; 与单独使用加权STA/LTA方法和峰度AIC法相比, 小波包-峰度AIC法具有更强的抗噪能力, 识别精度更高; 当初至清晰时, 小波包-峰度AIC法自动识别与人工识别的P波到时平均绝对差值为(0.077±0.075) s．      

基于支持向量机的地震体波震相自动识别及到时自动拾取

面对海量地震资料,自动准确地拾取震相并确定其到时的需求非常迫切.基于支持向量机技术,本文提出了使用两个分类器SSD和SPS自动识别地震体波震相并自动拾取其到时的方法.相比于传统的自动拾取方法,本文方法能够更准确地识别震相并区分P波和S波.进一步地,我们提出了利用台阵资料辅助识别震相的方案,有效地提高了地震震相拾取的准确率.     

利用标准时频变换方法在强噪声环境下无偏拾取地震P波、S波到时

地震P波、S波到时是精确分析地震水平位置、深度与速度结构等的重要参数,如何准确拾取P波和S波到时是地震学的一项重要的基础工作.大数据量与强噪声环境给地震到时的自动拾取带来了很大挑战.在频率域中可将信号与噪声分离,但会造成震相的偏移.针对上述问题,本文在STA/LTA、AIC方法的基础上,引入了标准时频变换(Normal...     

数字化地震记录震相自动识别的方法研究

针对目前震相自动识别方法不能自动给出震相识别区间,以及不能确定识别出震相名称的问题,运用震相的运动学特征,由平均速度模型和J-B走时表数据,自动计算近震、远震和极远震的震相走时及震中距。对多尺度小波分解进行单支重构作为识别不同震相的分析信号。先求出初至震相和最大面波到时,估算出震中距,然后找到S波或PP波到时,求出准确的震中距,即可自动给出各震相的识别区间,采用线性偏振法在给定区间中识别出震相的初至时刻。由于该方法采用的是先明确要找什么震相,再由该震相的走时确定寻找区间,所以找出的初至就是要寻找的震相,自然解决了识别出震相的名称问题,从而实现了对震相的全程自动识别。     

STA/LTA—AIC算法对地震P波震相拾取稳定性影响

选取区域地震台网记录的地震波形数据,使用STA/LTA算法与STA/LTA—AIC算法,进行地震P波震相初至到时自动拾取,对地方震及震中距较大的震相进行P波震相拾取效果分析,发现:STA/LTA算法对于地方震P波震相识别精度较高,与STA/LTA—AIC算法拾取的P波震相初至到时相差不大;震中距变大后,STA/LTA算法对P波拾取位置相对于最佳位置向后延迟,STA/LTA—AIC算法有效矫正了STA/LTA算法拾取位置的延迟问题,与人工拾取位置差别可忽略不计。     

区域地震信号自动识别方法及应用(英文)

地震信号的实时、自动、准确识别对于地震自动速报和地震预警十分重要。仿真信号试验分析表明,观测数据的四阶统计量函数(BKCF)对信号与噪声在能量和(或)频率方面的微弱差异变化具有较高的分辨能力。以此为基础,本文提出了一种新的自动探测区域地震事件的方法和测定直达波震相到时的BKCF-AIC方法。为了进一步提高波震相到时测定的精度,本文首先对指定时段的P-波记录进行偏振特性分析,其次对含有P波的S波记录进行偏振滤波处理,再次应用上述方法测定震相到时。与传统算法相比,基于山东测震台网记录的区域地震震例分析结果表明,使用本文提出的方法能够大幅度降低地震事件误检、漏检率,进一步提高了震相识别精度。     

近震S波震相实时自动识别方法研究

提出了一种用于地震早期预警的S波震相实时自动识别方法. 该方法不对原始信号进行任何滤波处理, 直接对三分向记录进行计算分析. 首先根据P波前0.5 s数据的卓越频率计算适用于该三分向记录的窗长, 采用由偏斜角和水平能量与总能量比值的平方积作为确定S波识别区间的特征函数, 将特征函数已有数据的5倍均值和5倍方差之和作为识别区间的触发阈值; 然后采用VAR-AIC方法对两个水平分向识别区间的数据分别计算分析, 对两个识别结果进行判断, 最终确定S波初动时刻. 经过对118个三分向记录的实际应用验证, 通过自动识别结果与人机交互震相识别结果相比, 本文方法对于S波相对P波尾波信噪比大于5 dB的地震记录, 其识别误差小于0.1 s的概率高达89.39%.      

广义S变换地震信号时频分析

广义S变换继承和发展了S变换的优点,其引入了λ、p两个参数,对S变换中固定的基本小波进行了推广和扩展,使得其时窗更灵活多变,在一定程度上提高了非平稳信号分析处理的时频分辨率。通过线性调频信号模拟测试对比,证实了广义S变换的信号聚焦能力明显优于S变换,能够更好地刻画信号的局部特性。最后,以某区二维单炮地震数据为例进行了地震信号的广义S变换时频分析,取得了良好效果。     

基于三重震相拟合的华南地区上地幔P波与S波速度结构

利用中国数字测震台网（CDSN）记录到的台湾地区两个地震事件6°~30°震中距范围的三重震相波形资料,基于观测与理论波形拟合法,获到华南地区上地幔P波和S波的最佳波形拟合速度模型及其V_P/V_S比值.与AK135模型相比,华南地区410 km深度上方存在明显低速层：S波低速区厚度约为70 km,速度降为2%~5%;而P波低速区厚度为70~230 km,速度降为5%~6%.另外,410 km间断面整体表现为一个梯度层,厚度约为10~40 km,V_P跃增量为4.0%~5.4%,而V_S跃增量为2.6%~11.7%.研究区内,低速层的V_P和V_S异常值大小和410 km间断面速度跃变量由北向南逐步减小.结合以往的接收函数和地震层析成像结果,华南地区410 km间断面上方的低速区可能与太平洋俯冲板块脱水有关.     

南海中北部陆缘横波速度结构及其构造意义

纵横波联合勘探可以得到更多关于岩石圈层岩性、物性等介质属性方面的信息,有效提高地壳物质组成的约束性.在纵波速度结构模型的基础上,通过射线追踪和走时拟合对OBS2006-3地震剖面径向分量的转换震相进行了横波速度结构模拟.结果表明:沉积层1、沉积层2的横波速度分别为0.7～0.9 km/s和1.6～1.7 km/s,波速...     

地震与爆破的小波包识别判据研究

利用sym5小波包基函数对小震级天然地震和人工爆破波形进行4层小波包分解并绘制了时频谱图.通过时频谱图可直观得出, 爆破频率成分简单, 时频谱聚集性较好. 为寻求定量的识别指标, 综合P波和S波小波包变换结果, 提出并定义了P/S能量比. 分析识别效果较好的 P/S能量比判据得出爆破的P波主频集中在频段3.125—9.375 Hz处, 地震频率成分较复杂, S波在高频12.5—23.4375 Hz处也较发育, 在这些频段上, 爆破的P波与S波差异要大于地震的P波与S波差异. 作为小波包判据研究的补充, 文中也提取分析了P波的能量比与S波的能量比. 能量比判据识别结果表明, 人工爆破与天然地震的频率成分存在差异, 通过小波包变换能够提取有效的识别判据.      

共转换点道集S变换谱相似性的多波地震数据横波波场提取技术

横波速度动校正后的共转换点(CCP)道集内,同时刻的各道横波信号S变换(ST)谱与其叠加道ST谱具有相似关系.因此,可基于这种相似关系设计自适应滤波器来提取多波地震数据中的横波波场.首先对共中心点(CMP)道集应用纵波速度动校正并在各道减去叠加道来去除数据中的纵波波场;然后在CCP道集应用横波速度动校正,将地震道振幅水平调整至叠加道振幅水平并做S变换,以叠加道ST谱为参考对地震道ST谱进行自适应滤波,去除数据中的残余纵波和噪声;最后,将滤波结果的振幅水平恢复至滤波前振幅水平.理论和实际数据试算表明,本文方法可有效提取多波地震数据中的横波波场,为多波多分量横波数据处理提供新思路.     

Relationship between P- and S-wave velocities and geological properties of near-surface sediments of the continental slope of the Barents Sea

Seismic velocities ( V _p and V _s) of compressional (P-) and shear (S-) waves are important parameters for the characterization of marine sediments with respect to their sedimentological and geotechnical properties. P- and S-wave velocity data of near-surface marine sediments (upper 9 m) of the continental slope of the Barents Sea are analysed and correlated to sedimentological and geotechnical properties. The results show that the S-wave velocity is much more sensitive to changes in lithology and mechanical properties than the P-wave velocity, which is characterized by a narrow range of values. The correlation coefficients between S-wave velocity and silt and clay content, wet bulk density, porosity, water content and shear strength are higher than 0.5 while the correlation coefficients of P-wave velocity and the same parameters are always lower than 0.4. Although the relationship between V _s and clay content has been widely described, the data show that V _s is better correlated with silt content than with clay content for the sediments of the area investigated. However, they show different trends. While V _s increases with increasing clay content, it decreases with increasing silt content.     

一种基于贝叶斯理论的P波双参数地震预警方法

为了及时、准确、可靠地发出地震预警信息,本文提出了一种基于贝叶斯理论的地震P波双参数预警方法。选取中国地震台网中心记录的四川地区的355条地震数据,统计P波触发后前3 s的P波平均周期τ_c和位移幅值P_d,结合贝叶斯理论建立震级和峰值加速度的预测模型,并以震级M4.5和峰值加速度为120 cm/s2为预警阈值,建立了地震危害性判别模型。与传统拟合方法进行对比仿真分析,并以汶川M_S8.0地震为震例,进行地震危害性判别实验与分析。实验结果表明:本文提出的基于贝叶斯理论的地震 P 波双参数预警方法比传统拟合方法地震漏报率低 15.15%,可以快速、准确地估计震级与峰值加速度,并有效地评估地震的危害性,能够为地震监测预警提供数据支持和决策依据。      

东北日本地震波速度、VP/VS和各向异性结构:对俯冲带水迁移过程的探讨

本文利用区域地震初至波到时数据,通过地震层析成像研究获得了东北日本俯冲带上地幔（深至约150 km）的P波速度（V_P）、S波速度（V_S）、V_P/V_S和P波各向异性结构.结果表明,低速及高V_P/V_S比异常体主要分布在火山下方的下地壳和地幔楔中,其与低频地震的分布吻合,该区域与俯冲板块脱水所释放的流体及其导致的部分熔融密切相关;俯冲的太平洋板块内可能由于脱水脆化导致的双层地震带区域则没有表现出整体的高V_P/V_S值,其可能与俯冲板块内部含水矿物含量有关;俯冲板块内双重地震带区域及上覆地幔楔薄层主要表现为与海沟平行的方位各向异性和正的径向各向异性,其可能是由于含水矿物的脱水使橄榄石晶格结构发生了从A型到B型的变化所引起的.我们研究表明,结合地震波速度和各向异性结构能够加深对俯冲带内水运移过程的认识.     

Simulation of improved pure P-wave equation in transversely isotropic media with a horizontal symmetry axis

Characterizing the expressions of seismic waves in elastic anisotropic media depends on multiparameters. To reduce the complexity, decomposing the P-mode wave from elastic seismic data is an effective way to describe the considerably accurate kinematics with fewer parameters. The acoustic approximation for transversely isotropic media is widely used to obtain P-mode wave by setting the axial S-wave phase velocity to zero. However, the separated pure P-wave of this approach is coupled with undesired S-wave in anisotropic media called S-wave artefacts. To eliminate the S-wave artefacts in acoustic waves for anisotropic media, we set the vertical S-wave phase velocity as a function related to propagation directions. Then, we derive a pure P-wave equation in transversely isotropic media with a horizontal symmetry axis by introducing the expression of vertical S-wave phase velocity. The differential form of new expression for pure P-wave is reduced to second-order by inserting the expression of S-wave phase velocity as an auxiliary operator. The results of numerical simulation examples by finite difference illustrate the stability and accuracy of the derived pure P-wave equation.     

从阵列声波测井波形处理地层纵、横波时差的新方法

阵列声波测井波形中包含丰富的地层纵、横波速度信息.如何从这些波形中提取出这些速度信息,得到比较准确的地层纵、横波时差是资料应用人员非常关心的问题.目前,最常用的纵、横时差处理方法是STC———慢度—时间相关法.该方法只有在波形的速度频散不严重时比较有效.而井内声波传播的基本特征是具有频散特征.利用Ma-trix pencil方法可以求出阵列声波测井波形的二维谱.借助于波形的二维谱,本文给出一种求地层的纵、横波时差的方法.本方法充分利用井内声波传播的规律,有效地减小了频散对纵、横波时差计算的影响,得到的纵、横波时差值精度比较高.     

利用纵、横波震源谱参数研究波速比——以岷县—漳县MS6.6地震为例

利用甘肃省区域地震台网提供的地震波形,计算得到甘肃东南地区2010年1月至2014年6月183个M_L≥2.5地震纵、横波的拐角频率和零频极限等震源谱参数,分析2013年7月22日岷县—漳县M_S6.6地震前后纵、横波拐角频率比值和零频极限比值的变化特征,间接反映震源处波速比的变化。结果如下:(1)拐角频率比值的平均值约为1.32,而由零频极限比值计算得到的波速比平均值约为1.55,后者高于前者;(2)二者得到波速比的时间变化趋势基本一致:两个波速比在岷县地震前后呈现出"明显升高-下降-震前急剧下降-震后缓慢恢复"的变化过程;(3)二者得到的波速比空间分布基本一致:岷县地震震中附近及合作、舟曲等地区在震前出现低值异常,震中附近地区的波速比在震后有所回升。