基于矢量预测的SIMO系统混合相位地震子波提取方法研究

针对高阶统计量混合相位地震子波提取方法的局限性, 提出一种基于矢量预测的单输入多输出系统(SIMO)的混合相位地震子波提取方法. 该方法利用二阶循环平稳统计量包含的系统相位信息, 将CMP道集视为一个单输入多输出系统的输出, 利用道集中相邻两道或多道数据通过矢量预测来构建反子波计算的方程式, 进一步进行混合相位子波提取. 利用提取出的子波相位信息对CMP道集进行纯相位滤波, 能够替代相位校正技术, 并且利用获取的反褶积算子对CMP道集进行反褶积处理, 使提高分辨率后的不同道子波振幅、 频率、 波形相一致, 提高叠加的质量. 模型试算和实际资料处理结果表明, 文中方法适用于任意相位的子波提取及反褶积处理, 且处理精度较高, 具有较高的实际应用价值.      

Optimization of sub-band coding method for seismic data compression

Seismic data volumes, which require huge transmission capacities and massive storage media, continue to increase rapidly due to acquisition of 3D and 4D multiple streamer surveys, multicomponent data sets, reprocessing of prestack seismic data, calculation of post‐stack seismic data attributes, etc. We consider lossy compression as an important tool for efficient handling of large seismic data sets. We present a 2D lossy seismic data compression algorithm, based on sub‐band coding, and we focus on adaptation and optimization of the method for common‐offset gathers. The sub‐band coding algorithm consists of five stages: first, a preprocessing phase using an automatic gain control to decrease the non‐stationary behaviour of seismic data; second, a decorrelation stage using a uniform analysis filter bank to concentrate the energy of seismic data into a minimum number of sub‐bands; third, an iterative classification algorithm, based on an estimation of variances of blocks of sub‐band samples, to classify the sub‐band samples into a fixed number of classes with approximately the same statistics; fourth, a quantization step using a uniform scalar quantizer, which gives an approximation of the sub‐band samples to allow for high compression ratios; and fifth, an entropy coding stage using a fixed number of arithmetic encoders matched to the corresponding statistics of the classified and quantized sub‐band samples to achieve compression. Decompression basically performs the opposite operations in reverse order. We compare the proposed algorithm with three other seismic data compression algorithms. The high performance of our optimized sub‐band coding method is supported by objective and subjective results.     

地震数据处理中的相位校正技术综述

提高地震剖面的信噪比和分辨率始终是地震数据处理中研究的重要内容.在动静校正后的同一道集内振幅和相位相同的假设条件下,地震叠加技术能明显提高资料的信噪比.实际资料中相位差问题是影响有效反射信号实现同相叠加的重要因素之一.相位校正技术能够消除子波的相位谱差异,使子波接近或达到零相位,从而达到提高叠加剖面质量的目的.对相位校正的国内外研究现状,影响相位变化的因素,相位校正的方法原理以及判别准则进行了系统的总结和概述.通过分析总结得出:随着勘探形势的发展,地震资料中的相位问题将越来越受到人们的关注和重视,对相位的研究也将成为今后的攻关课题和研究热点.     

A procedure to reduce side lobes of reflection wavelets: A contribution to low frequency information

Side lobes of the wavelets arise from the lack of low frequency content in a reflection wavelet. They tend to increase the time span of an individual reflection event and interfere with the other primary reflections or side lobes. Furthermore, their trace-by-trace consistency may produce pseudo-reflections and may cause misinterpretations of the side lobes as weak reflections.A procedure in order to improve the low frequency content of the seismic traces by suppressing the side lobe amplitudes based on the complex trace envelope is proposed. Using the average energies of the seismic trace and its envelope, the polarity table of the trace is obtained and used to correct the phase of the envelope. The resultant trace is termed “side lobe reduced (SLR) trace”. The method can be applied to the stack or migrated seismic data by a trace-by-trace basis. The only required parameter of the method is the moving average operator length which is used to calculate average energies of the input traces. In general, shorter operator lengths yield better results when the dominant frequency of the input increases.Results from synthetics and real seismic data sets show that the procedure improves the low frequency components of the input trace and side lobes in the output SLR trace are significantly suppressed. The method may be considered as a seismic amplitude attribute, which aids the interpreter to obtain the true seismic signature of the geological formations by removing the side lobes of the wavelet and restoring the low frequency components if the lower frequencies of deeper reflections are of primary concern.     

基于高阶统计的非最小相位地震子波恢复

利用高阶统计包含信号的相位信息特性,并基于信号的四阶累积量及其四阶谱,提出一种地震信号的非最小相位子波的估计方法.在任意高斯噪声环境下,对地震子波进行最小相位和最大相位谱分解,两部分信息完全可以从四阶谱中恢复.计算机数值模拟实验证实了方法的有效性.     

Wavelet correlation filter for wide-angle seismic data

A new filtering technique for single‐fold wide‐angle reflection/refraction seismic data is presented. The technique is based on the wavelet decomposition of a set of adjacent traces followed by coherence analysis. The filtering procedure consists of three steps. In the first, a wavelet decomposition of traces into different detail levels is performed. In the second, the coherence attributes for each level are evaluated by calculating cross‐correlation functions of detail portions contained in a space–time moving window. Finally, the filtered traces are obtained as a weighted reconstruction of the trace details. Each weight is obtained from the coherence‐attributes distribution estimated in a proper interval. A sequence of tests is then conducted in order to select possible optimum or unsuitable wavelet bases. The efficiency of the filter proposed was assessed by calculating some properly designed parameters in order to compare it with other standard de‐noising techniques. The proposed method produced a clear signal enhancement in high‐density wide‐angle seismic data, thus proving that it is a useful processing tool for a reliable correlation of seismic phases.     

Depth imaging in anisotropic media by symmetric non-stationary phase shift

We present a new depth‐imaging method for seismic data in heterogeneous anisotropic media. This recursive explicit method uses a non‐stationary extrapolation operator to allow lateral velocity variation, and it uses the relationship between phase angle and the spectral coordinates of seismic data to allow velocity variation with phase angle. A qualitative comparison of migration impulse responses suggests that, for an equivalent cost, the symmetric non‐stationary phase‐shift (SNPS) operator is superior to the phase‐shift plus interpolation (PSPI) operator, for very large depth intervals. To demonstrate the potential of the new method, seismic data from a physical model acquired over a transversely isotropic medium are imaged using a shot‐record migration based on the SNPS operator.     

The attenuated Ricker wavelet basis for seismic trace decomposition and attenuation analysis

The widely used wavelets in the context of the matching pursuit are mostly focused on the time–frequency attributes of seismic traces. We propose a new type of wavelet basis based on the classic Ricker wavelet, where the quality factor Q is introduced. We develop a new scheme for seismic trace decomposition by applying the multi-channel orthogonal matching pursuit based on the proposed wavelet basis. Compared with the decomposition by the Ricker wavelets, the proposed method could use fewer wavelets to represent the seismic signal with fewer iterations. Besides, the quality factor of the subsurface media could be extracted from the decomposition results, and the seismic attenuation could be compensated expediently. We test the availability of the proposed methods on both synthetic seismic record and field post-stack data.     

陇西及周边地区地震活动与震磁效应的高阶统计量研究

高阶统计量是研究非高斯过程,非最小相位信号和非线性系统的有力工具,其应用领域已涉及通信、地球物理、生物医学、故障诊断等。本文对嘉峪关和乾陵2个地磁台站1998～2002年期间的地磁垂直分量(Z)资料运用高阶统计量方法进行处理,探讨陇西及其周边地区地震活动与地磁场变化之间相互关系,并将其结果与相关分析和线性拟合方法的结果进行对比。结果发现,高阶统计量异常一般早于两台Z分量相关系数和嘉峪关台Z分量线性拟合差异常1～2个月出现,且三阶矩(三阶累积量)的异常变化幅度在5～10之间,四阶矩的异常幅度在50以上,甚至达到150,四阶累积量的异常幅度在10～60之间。这些表明在震磁效应的统计分析研究中引入高阶统计量方法的必要性及其在地震预报中的潜力和良好前景。     

High‐definition frequency decomposition

Spectral decomposition is a widely used technique in analysis and interpretation of seismic data. According to the uncertainty principle, there exists a lower bound for the joint time–frequency resolution of seismic signals. The highest temporal resolution is achieved by a matching pursuit approach which uses waveforms from a dictionary of functions (atoms). This method, in its pure mathematical form can result in atoms whose shape and phase have no relation to the seismic trace. The high‐definition frequency decomposition algorithm presented in this paper interleaves iterations of atom matching and optimization. It divides the seismic trace into independent sections delineated by envelope troughs, and simultaneously matches atoms to all peaks. Co‐optimization of overlapping atoms ensures that the effects of interference between them are minimized. Finally, a second atom matching and optimization phase is performed in order to minimize the difference between the original and the reconstructed trace. The fully reconstructed traces can be used as inputs for a frequency‐based reconstruction and red–green–blue colour blending. Comparison with the results of the original matching pursuit frequency decomposition illustrates that high‐definition frequency decomposition based colour blends provide a very high temporal resolution, even in the low‐energy parts of the seismic data, enabling a precise analysis of geometrical variations of geological features.     

A COMPARISON OF STATISTICAL AND DETERMINISTIC WIENER DECONVOLUTION OF DEEP-TOW SEISMIC DATA*

Wiener ‘spiking’ deconvolution of seismic traces in the absence of a known source wavelet relies upon the use of digital filters, which are optimum in a least-squares error sense only if the wavelet to be deconvolved is minimum phase. In the marine environment in particular this condition is frequently violated, since bubble pulse oscillations result in source signatures which deviate significantly from minimum phase. The degree to which the deconvolution is impaired by such violation is generally difficult to assess, since without a measured source signature there is no optimally deconvolved trace with which the spiked trace may be compared. A recently developed near-bottom seismic profiler used in conjunction with a surface air gun source produces traces which contain the far-field source signature as the first arrival. Knowledge of this characteristic wavelet permits the design of two-sided Wiener spiking and shaping filters which can be used to accurately deconvolve the remainder of the trace. In this paper the performance of such optimum-lag filters is compared with that of the zero-lag (one-sided) operators which can be evaluated from the reflected arrival sequence alone by assuming a minimum phase source wavelet. Results indicate that the use of zero-lag operators on traces containing non-minimum phase wavelets introduces significant quantities of noise energy into the seismic record. Signal to noise ratios may however be preserved or even increased during deconvolution by the use of optimum-lag spiking or shaping filters. A debubbling technique involving matched filtering of the trace with the source wavelet followed by optimum-lag Wiener deconvolution did not give a higher quality result than can be obtained simply by the application of a suitably chosen Wiener shaping filter. However, cross correlation of an optimum-lag spike filtered trace with the known ‘actual output’ of the filter when presented with the source signature is found to enhance signal-to-noise ratio whilst maintaining improved resolution.     

Nonlinear structure‐enhancing filtering using plane‐wave prediction*

Attenuation of random noise and enhancement of structural continuity can significantly improve the quality of seismic interpretation. We present a new technique, which aims at reducing random noise while protecting structural information. The technique is based on combining structure prediction with either similarity‐mean filtering or lower‐upper‐middle filtering. We use structure prediction to form a structural prediction of seismic traces from neighbouring traces. We apply a non‐linear similarity‐mean filter or an lower‐upper‐middle filter to select best samples from different predictions. In comparison with other common filters, such as mean or median, the additional parameters of the non‐linear filters allow us to better control the balance between eliminating random noise and protecting structural information. Numerical tests using synthetic and field data show the effectiveness of the proposed structure‐enhancing filters.     

利用EMD和子波振幅谱与相位谱关系的时变子波提取方法

本文首先分析了地震波在黏弹介质的传播规律,基于黏弹介质地震波动方程总结了时变子波振幅谱和相位谱的关系,从而得出结论,准确估计子波相位谱初值和不同时刻的子波振幅谱是实现时变子波准确提取的必要条件.在此基础上,针对传统方法限制子波振幅谱形态且受限于分段平稳假设的问题,提出了一种利用EMD(Empirical Mode Decomposition)和子波振幅谱与相位谱关系的时变子波提取方法,根据子波对数振幅谱光滑连续而反射系数对数振幅谱振荡剧烈的特点,采用EMD方法将不同时刻地震记录的对数振幅谱分解为一组具有不同振荡尺度的模态分量,通过滤除振荡剧烈分量、重构光滑连续分量提取时变子波振幅谱;再应用子波振幅谱和相位谱的关系提取时变子波相位谱,将分别提取的振幅谱和相位谱逐点进行合成,最终实现时变子波的准确提取.本文方法不需要求取Q值,适用于变Q值的情况,具有良好的抗噪性能.数值仿真和叠后实际资料处理结果表明,相比传统的分段提取方法,利用本文方法提取的时变子波准确度更高,研究成果对提高地震资料分辨率具有重要意义.     

Klauder wavelet removal before vibroseis deconvolution

The spiking deconvolution of a field seismic trace requires that the seismic wavelet on the trace be minimum phase. On a dynamite trace, the component wavelets due to the effects of recording instruments, coupling, attenuation, ghosts, reverberations and other types of multiple reflection are minimum phase. The seismic wavelet is the convolution of the component wavelets. As a result, the seismic wavelet on a dynamite trace is minimum phase and thus can be removed by spiking deconvolution. However, on a correlated vibroseis trace, the seismic wavelet is the convolution of the zero-phase Klauder wavelet with the component minimum-phase wavelets. Thus the seismic wavelet occurring on a correlated vibroseis trace does not meet the minimum-phase requirement necessary for spiking deconvolution, and the final result of deconvolution is less than optimal. Over the years, this problem has been investigated and various methods of correction have been introduced. In essence, the existing methods of vibroseis deconvolution make use of a correction that converts (on the correlated trace) the Klauder wavelet into its minimum-phase counterpart. The seismic wavelet, which is the convolution of the minimum-phase counterpart with the component minimum-phase wavelets, is then removed by spiking deconvolution. This means that spiking deconvolution removes both the constructed minimum-phase Klauder counterpart and the component minimum-phase wavelets. Here, a new method is proposed: instead of being converted to minimum phase, the Klauder wavelet is removed directly. The spiking deconvolution can then proceed unimpeded as in the case of a dynamite record. These results also hold for gap predictive deconvolution because gap deconvolution is a special case of spiking deconvolution in which the deconvolved trace is smoothed by the front part of the minimum-phase wavelet that was removed.     

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.     

Rough seas and time-lapse seismic

Time‐lapse seismic surveying has become an accepted tool for reservoir monitoring applications, thus placing a high premium on data repeatability. One factor affecting data repeatability is the influence of the rough sea‐surface on the ghost reflection and the resulting seismic wavelets of the sources and receivers. During data analysis, the sea‐surface is normally assumed to be stationary and, indeed, to be flat. The non‐flatness of the sea‐surface introduces amplitude and phase perturbations to the source and receiver responses and these can affect the time‐lapse image. We simulated the influence of rough sea‐surfaces on seismic data acquisition. For a typical seismic line with a 48‐fold stack, a 2‐m significant‐wave‐height sea introduces RMS errors of about 5–10% into the stacked data. This level of error is probably not important for structural imaging but could be significant for time‐lapse surveying when the expected difference anomaly is small. The errors are distributed differently for sources and receivers because of the different ways they are towed. Furthermore, the source wavelet is determined by the sea shape at the moment the shot is fired, whereas the receiver wavelet is time‐varying because the sea moves significantly during the seismic record.     

Interferometric interpolation of sparse marine data

We present the theory and numerical results for interferometrically interpolating 2D and 3D marine surface seismic profiles data. For the interpolation of seismic data we use the combination of a recorded Green's function and a model‐based Green's function for a water‐layer model. Synthetic (2D and 3D) and field (2D) results show that the seismic data with sparse receiver intervals can be accurately interpolated to smaller intervals using multiples in the data. An up‐ and downgoing separation of both recorded and model‐based Green's functions can help in minimizing artefacts in a virtual shot gather. If the up‐ and downgoing separation is not possible, noticeable artefacts will be generated in the virtual shot gather. As a partial remedy we iteratively use a non‐stationary 1D multi‐channel matching filter with the interpolated data. Results suggest that a sparse marine seismic survey can yield more information about reflectors if traces are interpolated by interferometry. Comparing our results to those of f‐k interpolation shows that the synthetic example gives comparable results while the field example shows better interpolation quality for the interferometric method.     

Fast performance uncertainty estimation via pushover and approximate IDA

Approximate methods based on the static pushover are introduced to estimate the seismic performance uncertainty of structures having non‐deterministic modeling parameters. At their basis lies the use of static pushover analysis to approximate Incremental Dynamic Analysis (IDA) and estimate the demand and capacity epistemic uncertainty. As a testbed we use a nine‐storey steel frame having beam hinges with uncertain moment–rotation relationships. Their properties are fully described by six, randomly distributed, parameters. Using Monte Carlo simulation with Latin hypercube sampling, a characteristic ensemble of structures is created. The Static Pushover to IDA (SPO2IDA) software is used to approximate the IDA capacity curve from the appropriately post‐processed results of the static pushover. The approximate IDAs allow the evaluation of the seismic demand and capacity for the full range of limit‐states, even close to global dynamic instability. Moment‐estimating techniques such as Rosenblueth's point estimating method and the first‐order, second‐moment (FOSM) method are adopted as simple alternatives to obtain performance statistics with only a few simulations. The pushover is shown to be a tool that combined with SPO2IDA and moment‐estimating techniques can supply the uncertainty in the seismic performance of first‐mode‐dominated buildings for the full range of limit‐states, thus replacing semi‐empirical or code‐tabulated values (e.g. FEMA‐350), often adopted in performance‐based earthquake engineering. Copyright © 2009 John Wiley & Sons, Ltd.     

地震复谱分解技术及其在烃类检测中的应用

谱分解技术在地震解释领域已得到广泛应用,但常用的谱分解方法存在两方面的不足.一是时间分辨率低,难以对薄层进行刻画;二是在烃类检测中多解性强,难以区分流体类型.为了改善该问题,本文提出一种基于地震复谱分解技术的烃类检测方法.复谱分解是指用一个包含多个不同频率Ricker子波的复子波库对地震道进行分解,从而得到时变子波频率和相位信息的过程.借助稀疏反演技术复谱分解可以获得高分辨率的时频能量谱和时频相位谱.本文首先通过拟合算例验证了复谱分解方法刻画薄层的能力以及求取子波频率和相位的准确性.然后利用基于Kelvin-Voigt模型的黏弹波动方程数值模拟对衰减引起子波相位改变的原因进行了分析.最后通过实际资料应用展示了本文方法在储层预测中的高时间分辨率优势,验证了利用子波相位信息识别气藏的有效性.     

Rough seas and statistical deconvolution

The rough‐sea reflection‐response varies (1) along the streamer (2) from shot to shot and (3) with time along the seismic trace. The resulting error in seismic data can be important for time‐lapse imaging. One potential way of reducing the rough‐sea receiver error is to use conventional statistical deconvolution, but special care is needed in the choice of the design and application windows. The well‐known deconvolution problem associated with the non‐whiteness of the reflection series is exacerbated by the requirement of an unusually short design window – a requirement that is imposed by the non‐stationary nature of the rough‐sea receiver wavelet. For a synthetic rough‐sea data set, with a white 1D reflection series, the design window needs to be about 1000 ms long, with an application window about 400 ms long, centred within the design window. Although such a short design window allows the deconvolution operator to follow the time‐variation of the rough‐sea wavelet, it is likely to be too short to prevent the non‐whiteness of the geology from corrupting the operator when it is used on real data. If finely spatial‐sampled traces are available from the streamer, the design window can be extended to neighbouring traces, making use of the spatial correlations of the rough‐sea wavelet. For this ‘wave‐following’ approach to be fruitful, the wind (and hence the dominant wave direction) needs to be roughly along the line of the streamer.