厚风化层覆盖区转换波静校正方法

P-SV转换波处理与传统的P-P波处理有很大的不同,如S波静校正、CCP叠加、P-SV速度分析和偏移等,其中最大的难题就是S波静校正问题.S波速度基本不受潜水面的影响,与纵波静校正没有直接相关性,有时横波静校正量能达到纵波静校正量的十倍,用纵波静校正量乘以比例系数来解决横波静校正问题将导致较大误差.同一接收点X和Z分量存在一定的初至时差,该时差代表了P波和S波在低降速带的走时差,可以利用该时差和近地表纵横波速度比信息去除低降速带对横波的影响,得到准确的静校正量.本文利用多分量初至时差推导了较为精确的横波静校正公式,再结合共检波点叠加求取剩余静校正量的方法,形成了完整的转换波静校正配套方法.利用该方法对苏里格气田二维及三维多波地震资料进行了实际处理,数据处理结果证明了该方法的有效性,该方法尤其适用于其他方法难以奏效的风化层较厚地区的横波静校正量求解,该方法也同时考虑了长波长横波静校正问题.     

一种地震P波和S波初至时间自动拾取的新方法

地震P波、S波初至时间的拾取是地震波分析的一项基础性工作.本文提出了一种新的地震波初至时间自动拾取的方法:首先,把地震波的三分量时程曲线变换为一组空间向的能量变化率时程曲线;然后对能量变化率时程曲线进行STA/LTA(Short Time Average/Long Time Average,短时间的均值/长时间的均值)处理,拾取地震P波和S波的大致初至时间;最后提出采用一种二次方自回归模型对初至附近的能量变化率曲线进行二次方自回归处理,精确拾取出P波和S波的初至时间.本文采用了10组芦山地震的记录数据和150组汶川地震的记录数据对此方法的可靠性进行了检验.以人工拾取结果为参考,此方法具有很高的准确率和稳定性,同时,相比于常用的STA/LTA方法和AIC(Akaike Information Criterion,Akaike信息准则)方法,此方法在计算时间效率方面稍微逊色,但是对S波初至时间的拾取精度和可靠性更高.此方法丰富了地震P波、S波初至时间的自动拾取方法.     

Automatic detection of P- and S-wave arrival times: new strategies based on the modified fractal method and basic matching pursuit

In this work, new strategies for automatic identification of P- and S-wave arrival times from digital recorded local seismograms are proposed and analyzed. The database of arrival times previously identified by a human reader was compared with automatic identification techniques based on the Fourier transformation in reduced time (spectrograms), fractal analysis, and the basic matching pursuit algorithm. The first two techniques were used to identify the P-wave arrival times, while the third was used for the identification of the S-wave. For validation, the results were compared with the short-time average over long-time average (STA/LTA) of Rietbrock et al., Geophys Res Lett 39(8), (2012) for the database of aftershocks of the 2010 Maule M _w = 8.8 earthquake. The identifiers proposed in this work exhibit good results that outperform the STA/LTA identifier in many scenarios. The average difference from the reference picks (times obtained by the human reader) in P- and S-wave arrival times is ～ 1 s.     

STATIC CORRECTIONS FOR SHEAR WAVE SECTIONS*

The computation of static corrections requires information about subsurface velocities. This information can be obtained by different methods: surface wave analysis, short refraction lines, downhole times, uphole times and first arrivals from seismograms. For pure shear waves generated by SH sources the analysis of first arrivals from seismograms combined, if necessary, with short refraction lines has proved to be most accurate and economic. A comparison of first-arrival plots from P- and S-wave surveys of the same line measured in areas of unconsolidated sediments in northern Germany illustrates the characteristic differences between the two velocity models. P-waves show a marked velocity increase at the water table from about 600 to 1800 m/s. S-wave velocities of the same strata increase gradually from about 100 to 400 m/s. As a consequence, S-wave models are vertically and laterally more complex and, in general, show no significant velocity increase at a defined boundary as P-wave models do. Therefore, other suitable correction levels with specific velocities must be chosen. A comparison of “tgd-corrections” (correction time between geophone position and datum level) for P- and S-waves in areas of unconsolidated sediments shows that their ratio is different from the P-/S-velocity ratio for the respective correction level because of the greater depth of the S-wave refractor. Therefore, P- and S-waves are influenced by different near-surface anomalies, and time corrections calculated for both wave types are largely independent.     

基于S变换的高频全球导航卫星系统同震数据的震相识别研究

为研究高频全球导航卫星系统（GNSS）信号的P波和S波到时,本文利用2008年汶川M_S8.0地震震时部分站点记录到的1 Hz高频GNSS数据,采用广义S变换将同震信号进行二维时频域平面分解,以此对P波和S波到时进行识别。本文通过调整广义S变换中的调节因子λ_a和p得出,当λ_a＝1.05,p＝1.05时,变换所得图像中的P波和S波到时均较为明显。结果表明,S变换在高频GNSS震相识别中效果明显,可作为GNSS地震学应用中一项有效的技术手段。      

Characteristics of S-wave Envelope Broadening in the Changbaishan Tianchi Volcano

High-frequency S-wave seismogram envelopes of microearthquakes broaden with increasing travel distance,a phenomenon known as S-wave envelope broadening. Multiple forward scattering and diffraction for the random inhomogeneities along the seismic ray path are the main causes of S-wave envelope broadening,so the phenomenon of S-wave envelope broadening is used to study the inhomogeneity of the medium. The peak delay time of an S-wave,which is defined as the time lag from the direct S-wave onset to the maximum amplitude arrival of its envelope,is accepted to quantify S-wave envelope broadening. 204 small earthquake records in Changbaishan Tianchi volcano were analyzed by the S-wave envelope broadening algorithm. The results show that S-wave envelope broadening in the Changbaishan Tianchi volcano is obvious,and that the peak delay time of S-wave has a positive correlation with the hypocenter distance and frequency of the S-wave. The relationships between the S-wave peak delay time and the hypocenter distance for different frequency bands were obtained using the statistics method. The results are beneficial to the understanding of the S-wave envelope broadening phenomena and the quantitative research on the inhomogeneities of the crust medium in the Changbaishan Tianchi volcano region.     

长白山天池火山区S波包络展宽的特征

高频S波随着传播距离的增大其均方根(RMS)包络逐渐变宽,把这种现象称为S波的包络展宽现象。S波传播路径上随机分布的非均匀体对S波的多次前向散射和绕射作用是导致S波包络展宽现象的主要原因,因此可用S波的包络展宽现象来研究介质的非均匀性。采用S波包络峰值延时对S波包络展宽现象进行量化。S波包络峰值延时定义为S波初至与其RMS包络峰值最大值之间的时间差。文中对长白山天池火山区204个小震记录进行了S波包络展宽分析,研究发现长白山天池火山区小震S波包络展宽现象明显,S波包络峰值延时与震中距、频率呈正相关;运用统计学方法获取了长白山天池火山区不同频带下S波包络峰值延时与震中距之间的关系。研究结果对认识S波包络展宽的特征以及定量研究长白山天池火山区介质的非均匀性具有重要意义     

利用祁连山主动源资料研究2016年门源6.4级地震前后波速变化

采用叠加、互相关、插值拟合等方法对祁连山气枪主动源台网数据进行处理,针对2016年1月21日青海门源6.4级地震前后气枪激发地震P波、S波震相的走时变化特征进行了分析,结果表明,地震前约6个月时,震中附近3个台站的相对走时出现下降变化（走时减少）,至震前约3个月时低值异常恢复正常,之后再次出现走时下降变化,地震即发生于走时变化恢复过程中。S波走时变化最大下降幅度达18ms,震后走时变化逐渐恢复正常,且3个台站变化趋势较为一致,其中,距震中最近的台站的S波走时变化最明显（ZDY38台）,较远台站的走时变化幅度较小,其变化特征与震源区位置有关。走时缩短意味着速度增加,可能与区域应力积累间存在一定的关系。     

One-Dimensional qP-Wave Velocity Model of the Upper Crust for the West Bohemia/Vogtland Earthquake Swarm Region

The western part of the Bohemian Massif (West Bohemia/Vogtland region) is characteristic in the relatively frequent recurrence of intraplate earthquake swarms and in other manifestations of past-to-recent geodynamic activity. In this study we derived 1D anisotropic qP-wave model of the upper crust in the seismogenic West Bohemia/Vogtland region by means of joint inversion of two independent data sets - travel times from controlled shots and arrival times from local earthquakes extracted from the WEBNET seismograms. We derived also simple 1-D P-wave and S-wave isotropic models. Reasons for deriving these models were: (a) only simplified crustal velocity models, homogeneous half-space or 1D isotropic layered models of this region, have been derived up to now and (b) a significant effective anisotropy of the upper crust in the region which was indicated recently by S-wave splitting. Both our anisotropic qP-wave and isotropic P-and S-wave velocity models are constrained by four layers with the constant velocity gradient. Weak anisotropy for P-waves is assumed. The isotropic model is represented by 9 parameters and the anisotropic one is represented by 24 parameters. A new robust and effective optimization algorithm - isometric algorithm - was used for the joint inversion. A two-step inversion algorithm was used. During the first step the isotropic P- and S-wave velocity model was derived. In the second step, it was used as a background model and the parameters of anisotropy were sought. Our 1D models are adequate for the upper crust in the West Bohemia/Vogtland swarm region up to a depth of 15 km. The qP-wave velocity model shows 5% anisotropy, the minimum velocity in the horizontal direction corresponds to an azimuth of 170°. The isotropic model indicates the V_P/V_S ratio variation with depth. The difference between the hypocentre locations based on the derived isotropic and anisotropic models was found to be several hundreds of meters.     

天山-准噶尔地区地震层析成像与壳幔结构

通过对东西天山两条天然地震剖面的远震P波与近震P波走时资料以及2001-2007年新疆地方台网的近震走时资料分析,获得了天山-准噶尔地区的三维地震层析图像.层析反演使用的台站数140个,地震事件1132个,P波走时数24904条.检测板试验结果指出,东西天山剖面、天山轴部和东准噶尔地区P波速度扰动恢复能力相对高.纵波波...     

Seismic imaging of the entire arc of Tohoku and Hokkaido in Japan using P-wave, S-wave and sP depth-phase data

In order to better understand seismic structure and seismotectonics of the entire arc of Tohoku and Hokkaido in Japan, we combined arrival time data from earthquakes beneath Tohoku and Hokkaido land areas, and beneath the Pacific Ocean to determine the three-dimensional (3D) velocity structures (Vp and Vs) under the entire Northeast (NE) Japan-Kuril arc. We adopted 176,431 P-wave and 110,953 S-wave arrival times, from 5123 local earthquakes, and 2843 sP depth-phase data from 385 events that occurred beneath the Pacific Ocean. The 385 suboceanic events were accurately relocated by using P-wave, S-wave and sP depth-phase arrival time data jointly. The obtained results confirmed the major features delineated by previous studies and revealed some new features of the structural heterogeneity beneath NE Japan and the Kuril forearcs. High-velocity anomalies of the cold subducting Pacific slab and low-velocity anomalies in the hot mantle wedge were imaged clearly. Strong lateral heterogeneities were revealed on the upper boundary of the Pacific slab under the forearc region, which showed a good correlation with the spatial distribution of large interplate earthquakes. These results indicated that strong coupling sections (or asperities) and weak-coupled or decoupled patches might exist along the upper boundary of the Pacific slab. Widespread low-velocity anomalies were visible in the forearc mantle above the subducting Pacific slab, which might reflect serpentinization of the forearc mantle associated with the dehydration process of the subducting slab. Our results also showed a general tendency for seismic coupling in the asperities to be located around low-velocity areas on the slab boundary under the suboceanic region.     

鲜水河—安宁河断裂带深部构造特征、固体潮应力变化与地震触发相关性研究

本研究通过反演241,561条纵波和209,363条横波高质量的地震走时数据,获得了鲜水河、安宁河断裂带的P波、S波以及泊松比的三维结构模型,并结合相位分析与Schuster测试讨论了该断裂带上地震事件的触发与其深部构造以及固体潮所产生的应力变化三者之间的关系.结果表明,鲜水河—安宁河断裂带26.5°N、28°N、29.5°N以及31.5°N四个位置附近存在低速、高泊松比异常,可能暗示着深部流体(包括部分熔融物质)的上涌.这些大范围分布的流体提高了孕震层的孔隙流体压力,削弱了断层面之间摩擦力,在鲜水河—安宁河断裂带形成了一些具有高度地震活动性的区域.此外,相位分析与Schuster测试的结果表明,鲜水河—安宁河断裂带附近的地震触发与固体潮剪切应力的变化密切相关.综合分析的结果表明,固体潮与构造应力对地震触发可能存在某种“合”的关系.在分布着大量流体的区域,地震事件可能对固体潮更加敏感.     

中国大陆岩石标本在高压下的弹性波速的初步研究

对10种岩石的20个试件和1个氯化银试件测得了高压下纵波和横波的速度,其中3个试件的纵波测至20,000公斤/厘米~2。发现压力高于4000公斤/厘米~2以后,速度基本上随压力呈直线变化,但当压力高于10,000公斤/厘米~2时直线的斜率逐级降低。在横波速度测定中,为提高脉冲初至到时的测定精度,在发射和接收端分别加入耦合电路。利用S波波列的多点测量外推方法有效地提高了测量精度。最后对结果作了细致的误差分析。     

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

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

利用深反射地震数据构建的多阶面波频散曲线反演近地表横波速度结构——以跨班公湖—怒江缝合带深反射地震资料为例

深反射地震剖面法为了获取深部结构特征常常采取大的偏移距采集数据.目前公开发表的相关资料中,鲜有利用深反射地震炮集数据获取近地表的结构特征.为此,本文通过正演测试了相关数据处理流程,即利用有限差分正演了起伏地表模型的大偏移距地震单炮弹性波场特征,通过共检波点域面波信号F-K频谱叠加构建新方法,从深反射地震数据集中提取了高品质的多阶面波频散曲线,再利用多阶面波联合反演获得了近地表的结构特征.在前述正演流程基础上,利用跨越班公湖—怒江缝合带的SinoProbe深反射地震剖面中的实际炮集数据,求取了基阶和一阶瑞利波频散曲线,联合反演后得到近地表横波速度结构.该结果与初至波走时反演获取的纵波速度结构具有较好的一致性,且在近地表的浅层分辨率较纵波速度结构特征更高,而更与已有地质认识相吻合.本文提供的相关数据处理流程表明利用深反射地震炮集数据,也能够获取近地表浅层的横波速度结构.     

三维地震与OBS联合勘探揭示的神狐海域含水合物地层声波速度特征

以三维高分辨地震与海底高频地震仪（OBS）联合勘探数据为基础,获得海底之下沉积层的地震反射成像剖面及多波信息,并以此确定研究区含天然气水合物沉积层的纵、横波速度的变化特征.根据走时反演获得的横波速度与纵波速度对比分析发现,研究区海底之下500 m深度范围内的某些沉积层具有较高的纵横波速度,这一纵波速度升高区域与水合物稳定带对应,而纵波速度下降并且横波速度变化较小的区域,可能与游离气的存在相关.游离气的可能存在与基于这一区域2007年钻探测井结果的普遍认识不完全相符.     

应用三分量浅层地震反射方法探测隐伏活动断裂

多波多分量地震勘探方法在油气田勘探开发中日益受到重视,但在工程和活动断裂探测中过去还没有开展多波多分量地震勘探的先例.本文简要介绍了在北京黄庄—高丽营断裂上开展的三分量浅层地震反射试验数据采集和资料处理方法,以及取得的初步结果.试验结果表明,采用纵波震源或横波震源激发、三分量接收的浅层地震勘探方法,可获得较高信噪比的浅层地震纵波、横波和转换波剖面,特别是水平分量采集的横波和转换波可提供很有意义的地下结构和构造信息.综合利用纵波、横波和转换波剖面特征进行的地质分层和断层解释结果取得了与跨断层的钻孔联合地质剖面结果较好的一致性.     

河南平顶山平煤矿区天然地震和人工爆破的波形特征及识别判据

基于河南测震台网记录的平顶山平煤矿区及周边发生的天然地震和人工爆破资料,按照直观、可快速识别的要求,选用发震时刻、P波初动方向、振幅比、振幅与尾波持续时间比等方法,对该地区天然地震和人工爆破进行对比分析。结果表明,P波与S波最大振幅之比（Pm/Sm）是识别该地区地震类型的最有效判据;P波初动振幅与S波最大振幅之比（Pc/Sm）和S波最大振幅与尾波持续时间之比（Sm/T）2种方法可以作为辅助判据,以提高识别地震类型的准确性。     

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

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

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

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