A new stacking technique in ambient noise tomography

We have developed a new stacking technique in ambient noise tomography to obtain high-quality dispersion curves of Rayleigh waves.This technique is used to stack the vertical components of the Estimated Green Functions(EGFs) obtained respectively from cross correlation of the ambient noise data recorded by a remote seismic station and one of the short distance seismic stations of a seismic array.It is based on a phase-matched filter and is implemented by a four-step iterative process:signal compression,stacking,signal extraction and signal decompression.The iterative process ends and gives the dispersion curve of Rayleigh wave when the predicted one and the processing result converge.We have tested the method using the vertical components of synthetic Rayleigh wave records.Results show that this new stacking method is stable and it can improve the quality of dispersion curves.In addition,we have applied this method to real data.We see that the results given by our new technique are obviously better than the ones employing the traditional method which is a three-step process:signal compression,signal extraction and signal decompression.In conclusion,the new method proposed in this paper can improve the signal to noise ratio of EGFs,and can therefore potentially improve the resolution of ambient noise tomography.     

台站地震资料的时频域自适应极化分析和滤波

提出一种自适应协方差的时频域极化滤波方法。该方法在广义S变换时频方法的基础上,构造时频域自适应协方差矩阵,通过特征分析计算时频域瞬时极化参数,设计极化滤波器,实现多分量地震极化分析和滤波。其优势在于协方差矩阵的分析时窗的长度由多分量地震数据的瞬时频率确定,可以自适应于有效信号的周期,在每个时频点计算极化参数不需要进行插值处理;结合时间频率信息,解决在时间域或频率域波形或频率重叠的信号具有明显的直观性。模型数据及实际三分量台站地震数据处理结果表明,该极化滤波方法在台站地震资料分析和处理方面具有很好的直观性和较高的分辨率。     

Dispersion calculation method based on S-transform and coordinate rotation for Love channel waves with two components

Dispersion analysis is an important part of in-seam seismic data processing, and the calculation accuracy of the dispersion curve directly influences pickup errors of channel wave travel time. To extract an accurate channel wave dispersion curve from in-seam seismic two-component signals, we proposed a time–frequency analysis method based on single-trace signal processing; in addition, we formulated a dispersion calculation equation, based on S-transform, with a freely adjusted filter window width. To unify the azimuth of seismic wave propagation received by a two-component geophone, the original in-seam seismic data undergoes coordinate rotation. The rotation angle can be calculated based on P-wave characteristics, with high energy in the wave propagation direction and weak energy in the vertical direction. With this angle acquisition, a two-component signal can be converted to horizontal and vertical directions. Because Love channel waves have a particle vibration track perpendicular to the wave propagation direction, the signal in the horizontal and vertical directions is mainly Love channel waves. More accurate dispersion characters of Love channel waves can be extracted after the coordinate rotation of two-component signals.     

希尔伯特-黄变换在地震资料处理中的应用

希尔伯特-黄变换是最新发展起来的处理非线性非平稳信号的时频分析方法.本文介绍了该方法的基本原理和几个关键技术,以及其优势所在,并将它应用于含随机噪音的地震信号和实际地震资料的分解.分解后的信号能够重构原始信号,重构后在信噪比和分辨率方面有了一定的提高,对地震资料的去噪有重要意义.     

GHM类正交多小波变换及其在地震资料去噪中的应用

多小波是对小波理论的一个新发展,它可以同时满足正交性、对称性、短支撑等良好的特性要求。文中介绍了多小波基本理论、多小波变换具体过程及预处理方法,提出了基于GHM类多小波变换的地震资料软阈值去噪方法,通过对合成数据和实际资料进行处理分析,表明多小波变换在有效压制随机噪声的同时,能较好地保留原信号的特征信息,是一种行之有效的去噪方法     

利用小波变换识别山西测震台网震相

为有效提高山西测震台网震相识别精度,在论述小波变换原理基础上,结合台网记录的震相特征,针对干扰大、能量弱的震相,采用小波变换进行分析。常见噪声干扰相对地震信号频率较高,小波变换能够有效压制噪声,提高震相识别精度;对于能量较弱的震相,可以利用小波变换分解信号,按频率逐层进行震相识别。分析认为,小波变换是识别山西地区震相的有效方法。     

线性化空间变换叠前与叠后地震记录去噪方法

本文提出了一种适应复杂构造地区地震资料叠前与叠后去噪新方法.该方法基于叠加剖面信息,利用正交多项式逼近技术,计算出描述地震同相轴走势的相关系数.并利用这些系数对叠加剖面或炮记录在横向作线性化变换,对变换后地震数据进行奇异值检测压制和F X域有效信号预测.通过理论数据检验表明：本方法能较好适应复杂构造,具有较强去随机干扰和野值能力,保真性较好等特点.该方法已经在实际资料处理生产中得到了应用,并取得了较好的应用效果.     

FFT和DFT在地震数据采集器幅频特性测试中的应用

在进行地震数据采集器的幅频特性测试中,当标定信号频率稍高时(约为采集器采样率1/5以上),采集器输出的正弦波形发生畸变,无法人工读取正弦波峰值.为了解决该问题,本文提出了基于离散傅里叶变换(DFT)的积分恢复算法,将该方法与快速傅里叶变化(FFT)谱分析方法同时应用于一台高精度地震数据采集器的幅频特性测试数据处理.结果表明这两种算法均能标准化、高精度地完成采集器幅频特性标定数据的计算.其中基于DFT的积分恢复算法的计算程序操作简单,能够高效率地完成数据采集器的幅频特性检定工作;并且该方法还可以应用于地震数据采集器的线性度测试.      

基于改进BISQ机制的双相HTI介质波传播伪谱法模拟与特征分析

改进BISQ(Biot-Squirt)机制在不引入特征喷流长度的情况下,将含流体孔隙介质中Biot流动和喷射流动两种重要的力学机制有机地结合起来,且各相关参数具有明确物理意义和可实现性.本文将改进BISQ机制一维孔隙流体压力公式推广到三维具有水平对称轴横向各向同性介质(HTI介质)情况,结合裂缝各向异性理论,给出了基于改进BISQ机制的双相HTI介质模型及其二维三分量波传播方程,采用伪谱法求解该方程,进行了不同相界、不同频率以及双层地质结构情况下该类介质中波场的数值模拟与特征分析.数值模拟结果表明:伪谱法模拟精度高,压制网格频散效果好,可以得到高精度的波场快照和合成记录;基于改进BISQ机制的双相HTI介质模型兼具裂缝各向异性特征和孔隙弹性特征,其为从双相各向异性理论角度深入研究裂缝性储层的地震响应奠定了理论基础.     

Body-wave dispersion: Measurement and interpretation

Generalizing previous studies on short-period data, it is shown that body-wave dispersion can be measured from broad-band records of earthquakes of moderate magnitude. The method is based on the direct measurement of the arrival time of the frequency components of a seismic wave, and the arrival time is defined by its expectation value. The frequency components of the signal are obtained through a narrow band-pass filtering process. Previous to any interpretation, a correction of the arrival time for instrument response and group delay of the filter is needed. In the first step, body-wave dispersion is related to an absorption band to account for intrinsic attenuation, and thereafter we generalize this interpretation by considering a cascade of filters to account for medium parameters (attenuation and a layered crust) and source parameters (source time function and finiteness of fault). An inversion scheme to obtain the filter parameters can be devised by following, in a formal way, the same procedure as for the case of surface wave dispersion.     

Windowless Q-factor tomography by the instantaneous frequency

The estimation of the Q factor of rocks by seismic surveys is a powerful tool for reservoir characterization, as it helps detecting possible fractures and saturating fluids. Seismic tomography allows building 3D macro-models for the Q factor, using methods as the spectral ratio and the frequency shift. Both these algorithms require windowing the seismic signal accurately in the time domain; however, this process can hardly follow the continuous variations of the wavelet length as a function of offset and propagation effects, and it is biased by the interpreter choice. In this paper, we highlight some drawback of signal windowing in the frequency-shift method, and introduce a tomographic approach to estimate the Q factor using the complex attributes of the seismic trace. We show that such approach is particularly needed when the dispersion is broadening the waveforms of signals with a long wave-path. Our method still requires an interpretative event picking, but no other parameters as the time window length and its possible smoothing options. We validate the new method with synthetic and real data examples, involving the joint tomographic inversion of direct and reflected signals. We show that a calibration of the frequency-shift method is needed to improve the estimation of the absolute Q factor, otherwise only relative contrasts are obtained.     

面波频散谱多模式高分辨率成像的多道信号比较法

高分辨率的面波频散谱成像是浅层地震勘探领域基于频散性质反演横波速度结构中的一个关键步骤.在天然地震探测领域,仅利用两个台站记录的线性信号比较法(LSC),被广泛用来计算面波的频散谱,并用于大尺度的面波层析成像.然而互相关的成像方式会造成频散谱在低频端较低的分辨率.非线性信号比较法(NLSC)利用指数函数克服了这个问题,同时极大地提高了频散谱的成像分辨率.然而,在研究中我们发现,仅利用两个台站的地震记录,并不能将面波的频散特性完整地考虑在内,导致LSC和NLSC方法对高阶模式的成像存在较大的误差.由于主动源面波勘探多道采集的方式,基于信号比较理论的多道信号比较法(MSC)充分利用多道地震信号,可以获得准确的多模式成像,然而该方法需要计算任意两道的频散谱并叠加,存在冗余的计算,导致计算效率较低.因此,本研究对MSC方法进行了相应的改进,通过追踪炮集记录上的面波波组提高了原方法的计算效率,同时,利用理论频散曲线进行叠加分析,验证了改进的MSC方法的正确性和有效性.通过与相移法、LSC和NLSC方法的对比分析,展示了MSC方法是一种准确的、高分辨率的面波多模式频散谱成像方法.实际地震资料的应用,揭示了MSC方法在浅层地震勘探中用于提取面波多模式频散信息的巨大潜力.     

地震数据有损压缩与失真分析

文中提出了应用双正交小波变换有损压缩地震数据的方法,该方法具有实现简单、在高压缩比情况下失真小等特点。以信噪比、剩余能量作为评价指标,分析了压缩引起的信号失真。结果表明压缩比在0~31范围内取值时,信噪比较高,剩余能量在99%以上,原始信号与重构信号有很好的相关度。即使在不太理想的情况下,将重构信号用于地震参数快速测定,对震级的计算误差最大为8.6×10-3,震中距的最大计算误差为3.69 km。     

Applications of Savitzky-Golay Filter for Seismic Random Noise Reduction

This article utilizes Savitzky–Golay (SG) filter to eliminate seismic random noise. This is a novel method for seismic random noise reduction in which SG filter adopts piecewise weighted polynomial via leastsquares estimation. Therefore, effective smoothing is achieved in extracting the original signal from noise environment while retaining the shape of the signal as close as possible to the original one. Although there are lots of classical methods such as Wiener filtering and wavelet denoising applied to eliminate seismic random noise, the SG filter outperforms them in approximating the true signal. SG filter will obtain a good tradeoff in waveform smoothing and valid signal preservation under suitable conditions. These are the appropriate window size and the polynomial degree. Through examples from synthetic seismic signals and field seismic data, we demonstrate the good performance of SG filter by comparing it with the Wiener filtering and wavelet denoising methods.     

基于两级尺度粗化的裂缝-多孔隙介质地震响应分析方法

实际地震勘探中,储层物性参数的差异是导致地震波响应特征发生变化的根本原因,而建立储层物性参数与地震响应特征之间的联系,需要跨越微观孔隙尺度、介观测井尺度以及宏观地震尺度等三个不同尺度空间.本文基于已知井的岩石物理实验数据和测井数据,利用复杂多孔隙介质理论将微观尺度孔隙岩石粗化到介观测井尺度,利用Backus平均理论将介观测井尺度的模型进一步粗化到宏观地震尺度,最终,得到地震尺度裂缝-多孔隙介质模型.其数值计算结果与测井数据和地震数据的对比表明：基于两级尺度粗化算法的裂缝多孔隙介质模型在给定参数下是有效的,且基于该模型的地震响应特征分析方法能够对储层的地震响应特征随物性参数的变化进行分析.     

基式地震仪（SKZ）初动特性分析

本文通过计算基式地震仪（SKZ）对于不同形态和不同初始斜率的14种归一化地震信号的响应,说明地震仪的初动特性与地震信号有很大关系。作者建议用伯拉格（Berlage）函数近似远场地动位移,并由基式地震仪的时间特性和伯拉格函数褶积计算出仪器响应。对于频率不同的伯拉格函数得出了仪器的初动特性。最后指出,用与初动分析相同的方法去分析所谓二动、三动特性,没有普遍使用意义。     

东北地震区地震活动趋势数学建模分析

本文按照活动地块边界带和地震区的划分方案,使用Mapsis软件选取了东北地震区1900年以来Ms4.7级以上地震资料,这样得到的地震资料可能具有一个完整的地质系统性。使用数学中19阶多项式模型,对地震资料进行最佳逼近拟合,并递推未来1～2年时间尺度预测,模型中的参数用最小二乘法求取,引用了一种数学方法来描述东北地震区的地震趋势规律。     

A method for absorption compensation based on adaptive molecular decomposition

In this paper, we present a new method for seismic stratigraphic absorption compensation based on the adaptive molecular decomposition. Using this method, we can remove most of the effects resulting from wavelets truncation and interference which usually exist in the common time-frequency absorption compensation method. Based on the assumption that the amplitude spectrum of the source wavelet is smooth, we first construct a set of adaptive Gabor frames based on the time-variant properties of the seismic signal to transform the signal into the time-frequency domain and then extract the slowly varying component (the wavelet’s time-varying amplitude spectrum) in each window in the time-frequency domain. Then we invert the absorption compensation filter parameters with an objective function defined using the correlation coefficients in each window to get the corresponding compensation filters. Finally, we use these filters to compensate the time-frequency spectrum in each window and then transform the time-frequency spectrum to the time domain to obtain the absorption-compensated signal. By using adaptive molecular decomposition, this method can adapt to isolated and overlapped seismic signals from the complex layers in the inhomogeneous viscoelastic medium. The viability of the method is verified by synthetic and real data sets.     

独立分量分析及其在地震信息处理中应用初探

独立分量分析(ICA)是最近才发展起来的一种统计学方法，旨在寻求对非高斯分布数据进行有效表示，使得各个分量在统计学上独立，或者尽最大可能地独立。许多应用中，这种表示意在获取数据的基本结构，包括特征提取和信号分离。本文给出ICA的基本理论和快速算法，并对FastICA稍作改进，在分析地震信号特点的基础上，对其在地震信息处理中应用进行初步探索，表明ICA在地震信号处理中具有应用前景。     

一种边缘保持的地震数据插值方法

在地震数据处理中,地震数据插值方法常常用来解决地震数据空间采样率低和不规则的问题.本文提出了一种基于边缘保持滤波器的地震数据插值方法.在该方法中,对于一个1D信号,逐点滑动一个处理窗口,将信号分成多个信号片段.对于某一个待恢复的缺失采样点,存在多个包含(或邻近)该采样点的信号片段可以用来估计它.采用多项式来拟合这些信号片段,并选择拟合误差最小的信号片段估计此缺失采样点,达到边缘保持的目的.对于2D地震信号,先沿不同方向扫描抽取1D信号,然后采用上述1D边缘保持插值算法分别进行处理,得到沿不同方向的插值结果.对于任一待插值采样点,选取对应拟合误差最小的方向的插值结果作为最后输出的2D数据的插值结果.理论模型和实际资料的处理结果表明,所提方法具有保边缘、抗假频及能够进行不规则数据重建等特点,既能有效的实现不规则地震数据的重建,又能很好的保持原有数据的边缘特征.