基于广义S变换的地震资料信噪分离方法

基于S变换具有良好的时频聚焦性和分时分频性,将可灵活选取窗函数的广义S变换引入到地震信号去噪处理中,系统研究了广义S变换在地震资料信噪分离中的应用。首先对地震数据进行广义S变换,然后对含噪频率剖面选取适当阈值压制噪声干扰,提取有效信号,最后重构得到去噪后的记录。经合成记录和实际地震资料处理实验证明,该方法能有效地进行信噪分离,提高地震剖面信噪比和分辨率。     

内蒙古中部地区形变主要干扰的时频响应特征分析

选取内蒙古中部地区定点形变观测中受降雨、大风、气压干扰以及地震波影响的典型数据,采用S变换,对该地区形变观测中4种主要干扰信号的时频响应特征进行分析与讨论。结果表明,降雨干扰信号主要集中在低频区域,其频率先增大后减小;大风干扰表现为高频干扰,时频域能量强度与观测曲线受风扰影响变化幅度成正比;气压干扰信号的优势频率分布在低频段内;地震波影响中VP垂直摆和JCZ地震计记录的地震波较为相似,但VP垂直摆记录的频段要小于JCZ地震计,VP垂直摆主要记录的是低频地震波。同时,4种干扰信号的时频域频率、能量强度随时间变化趋势与数据观测时域变化趋势相一致。     

2013年岷县—漳县6.6级地震前地磁日变化异常及机理分析

利用空间相关法对2013年岷县—漳县6.6级地震前后151个台站地磁垂直分量日变化进行分析,结果显示:(1)2013年7月3日、4日在震中附近分别出现16和17个台站的地磁Z分量日变化明显反相位异常,具有明显的分界线。(2)利用S变换方法对分界线两侧的天水和兰州台的Z分量日变化进行分频率相关系数计算,发现7月3日、4日的地磁Z分量的日变化均在6 h左右的周期出现负相关变化。其他时段未出现类似变化,认为岷县—漳县6.6级地震前地磁Z分量日变化异常产生的反相位异常主要在6 h左右。此外,还对异常产生的机理进行了解释。     

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

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

基于含可变因子广义S变换的瑞雷面波频散曲线提取方法

面波频散曲线的提取是面波资料处理的关键.鉴于S变换和广义S变换不能全时段兼顾频率分辨率和时间分辨率的缺点,提出了含可变因子广义S变换进行瑞雷面波频散曲线的提取方法.可变因子的引入使得高斯窗函数的宽度随频率发生变化时具有目的性,而不是简单地随着频率的增大而变窄.该方法可以有针对性地改善局部频段(特别是低频段和高频段)的频率分辨率及时间分辨率.通过理论模型和实际资料试算表明:含可变因子广义S变换提取面波频散曲线的方法具有可行性和实用性.     

基于改进广义S变换二次时频谱反褶积在南海工区的应用

传统的时频域反褶积模型虽然吻合地震子波在粘弹性介质传播过程中的衰减特征,但受拟合函数和时频分析窗函数的影响,使得其算法缺乏稳定性和准确性。为了避免这一问题,推导出一种基于改进广义S变换二次时频谱反褶积方法。利用改进广义S变换的灵活性来获取非平稳地震记录时频谱,再对其傅里叶变换得到二次时频谱,从而以一种低通滤波的方式提取子波时频谱,实现二次时频谱反褶积。经理论模型论证和实际数据处理结果表明,利用二次时频谱提取子波具有较高的准确性,同时能够有效地分辨薄层,在不用考虑Q值的条件下直接恢复地震波被地层吸收所引起的衰减。     

Reliability Evaluation on Weak Signal Extraction for Airgun Source Surveys

Based on the data synthesis simulation and the actual processing of the airgun seismic source signal, three quantitative indicators of signal-to-noise ratio, waveform correlation coefficient and phase offset, are superimposed. We systematically evaluate the functions of the following three stack methods including linearity, phase weighting and S-transform in the extraction of weak signals under strong background noise and quantitatively estimate the reliability of the stack results. Through the comprehensive discussion of the above three methods of stack results, the preliminary comparative analysis believes that the linear stack signal-to-noise ratio is low, but the waveform distortion is minimal; the phase-weighted superimposed signal-to-noise ratio is high and the phase offset is small, but the results of the waveform quality and linear stack are larger than the deviation; the S-transform stack has a relatively higher signal-to-noise ratio and a small loss of waveform amplitude, but there is a certain phase shift phenomenon. It is therefore suggested that linear stack technology should be used when the requirements of both waveform quality and time precision are high. However, the selection of the stack method when the airgun source excitation is limited should be emphasized. If high fidelity is required, the S-transform stack method should be selected; if the required time is high, accuracy can be selected by phase-weighted stack method to achieve reasonable extraction of weak signals.     

Rayleigh wave separation using high-resolution time-frequency polarization filter

Obtaining geological information from seismic data motivates researchers to innovate and improve seismic wave processing tools. Polarization-based methods have received much attention regarding their ability to discriminate between different phases of the seismic wave based on polarity. Combining the intuitive definition of polarity in the frequency domain (monochromatic waves) with the non-stationary properties provided by time-domain methods, time-frequency approaches are attracting widespread interest because they localize the information extracted from the seismic waves in the joint time and frequency domains. Due to the lack of high-resolution time-frequency maps, the time-frequency polarization approach was not able to resolve specific temporal polarity changes in the seismic signal. The main objective of this study was to devise a robust time-frequency-based polarization filtering method using high-resolution polarization attributes obtained directly from the sparse time-frequency map without using Eigen analysis or analytic signals. The method proposed here utilizes a computationally effective sparsity-based adaptive S-transform to obtain a high-resolution polarization map of an inherently non-stationary seismogram for the entire frequency content of the signal at different times. The superiority of the proposed method over the S-transform method was verified using synthetic and real data sets to calculate the polarization attributes in the time-frequency domain and separate the Rayleigh waves from the seismogram.     

山西运城振动事件S变换时频分析

基于S变换,对2005年以来山西南部运城及其附近地区的振动事件波形进行时频分析,并与山西南部地区的天然地震、爆破和塌陷地震波谱特征进行对比分析。结果表明：天然地震一般S波携带能量较大,衰减较慢,震相高低频成分呈现均匀分布;一般近台记录的人工爆破P波比S波发育,能量衰减较快,震相急促短暂;塌陷地震波列能量随时间和频率的展布相对集中,一般分布在频率较低区域;振动事件震相简单,被不同台站记录的波形能量优势分布及频率分布范围差异较大,衰减特征不明显。据此,基本可以排除山西南部及附近区域的振动事件为天然地震、爆破和塌陷事件的可能。     

巴伦台分量钻孔应变在呼图壁MS6.2地震前的短临异常初探

呼图壁地震前巴伦台分量钻孔应变出现显著压性短临异常变化，排除仪器、气象及环境等干扰因素，认为该异常可能为地震前兆异常。利用超限率和S变换时频分析方法提取该异常，结果显示，信号集中在1 440~10 min频段，NS分量最大信号强度为185×10-10，EW分量最大信号强度为140×10-10。探讨此类异常变化特征，并对其原因进行分析。