交错网格中基于波数域插值的波场分离方法研究（英文）

应用多分量地震资料进行成像时通常需要先做波场分离,然后再对分离的波型进行成像。其中,波场分离可以在空间域或波数域实现。然而,由于用交错网格有限差分进行弹性波场数值模拟时,用来进行波数域波场分离的质点振动速度分量定义在不同网格节点上,本文提出了利用波数域插值方法来估算同一网格节点所需质点振动速度值;进而给出了先进行波数域插值后进行波场分离的波数域保幅波场分离方案。数值实验结果表明波数域插值方法具有较高的插值精度且保幅波场分离方法具有较好的保幅性,将本文方法进一步应用于弹性波逆时偏移可以获得保幅性较好的成像结果且对存在一定程度速度误差情况具有较好的适应性。     

基于Contourlet域自适应Wiener阈值的同时震源波场分离

同时震源数据包含了多炮之间的串扰噪声,不能直接用于常规数据处理流程.因此,需要对混叠的波场进行分离得到常规采集的单炮记录.本文基于稀疏迭代反演分离,提出了一种具有尺度与空间自适应的Wiener阈值选取方法.该阈值选取方法能够根据不同迭代环境计算不同尺度下串扰噪声的方差和不同空间位置有效信号的方差,从而自适应调整阈值大小,最终通过对变换域系数进行收缩来达到去除串扰噪声的目的.理论模型数据和实际数据测试结果表明,本文方法能够快速有效地压制串扰噪声和保护弱有效信号,取得了比Contourlet域子带一致Wiener阈值方法和Curvelet域指数衰减阈值方法更好的分离效果.

     

Practical aspects of wavefield separation of two-component surface seismic data based on polarization and slowness estimates

The processing of multicomponent seismic data, carried out individually on the different wavetypes (P-, S- and converted waves), should result in an improved image of the subsurface. We examine the wavefield‐separation method proposed by Cho and Spencer . We discuss practical aspects related to the separation of interfering waves on two-component surface seismic data and illustrate these using synthetic data. A sliding spatial window is used for analysis. The choice of its width represents a trade-off between stabilizing the method in the presence of random noise and ensuring a good spatial resolution. No a priori knowledge of the subsurface is required, but locally, the characteristic parameters of the waves, i.e. horizontal slowness and polarization, are assumed to be constant inside the analysis window. These parameters are estimated at each frequency, but a statistical analysis provides a more robust estimate, especially in the presence of random noise. This approach also solves the problem of eigenvalue sharing and switching. Additional smoothing of the estimates according to a model may further improve the results. The width of the analysis window may be decreased if the waves inside the data window differ significantly in amplitude. The dominant wave in each case is separated from the lower-amplitude waves and subtracted from the data. This novel iterative approach thereby allows for the isolation of low-amplitude events.     

Surface wave attenuation based polarization attributes in time-frequency domain for multicomponent seismic data

In the paper, we propose a surface wave suppression method in time-frequency domain based on the wavelet transform, considering the characteristic difference of polarization attributes, amplitude energy and apparent velocity between the effective signals and strong surface waves. First, we use the proposed method to obtain time–frequency spectra of seismic signals by using the wavelet transform and calculate the instantaneous polarizability at each point based on instantaneous polarization analysis. Then, we separate the surface wave area from the signal area based on the surface-wave apparent velocity and the average energy of the signal. Finally, we combine the polarizability, energy, and frequency characteristic to identify and suppress the signal noise. Model and field data are used to test the proposed filtering method.     

Localization and characterization of cracks in clay‐rocks using frequency and time‐domain induced polarization

Induced polarization (IP) is a geophysical method that is potentially sensitive to the presence of cracks in porous rocks and therefore to damage. We performed time‐domain and frequency domain IP measurements at the Tournemire Underground Research Laboratory (URL, Aveyron, France) in areas where different types of cracks are observed. These cracks correspond to both tectonic fractures and new cracks associated with stress release and desiccation resulting from the excavation of a gallery. These measurements were performed both in eastern and northern galleries of the test site. The eastern gallery was excavated in 1996 while the northern gallery was excavated recently in 2008. This gives us the opportunity to study the electrical characteristics of the excavation damaged zone surrounding the galleries with respect to the age of the excavation. Longitudinal profiles were performed along the floor of the galleries with 48 Cu/CuSO₄ electrodes separated by a distance of 20 cm. Chargeability and resistivity were inverted using a Gauss‐Newton iterative approach assuming an isotropic heterogeneous clay‐rock material. The resulting IP tomograms show a correlation between high values of chargeability and the presence of calcite‐filled tectonic fractures. X‐ray analysis indicates that the presence of pyrite in these fractures is a potential source of the observed IP signals. The cracks associated with the mechanical damage of the formation exhibit low values of chargeability, on the same order of magnitude than the chargeability of the clay‐rock matrix and are therefore hardly observable. A smaller IP response associated with the presence of these cracks is observed in the older gallery and this observation is qualitatively related to the desaturation process associated with these cracks. In a specific area of one of the galleries, the presence of calcareous nodules is observed to be an important source of anomalous chargeability. This signature seems to be associated with the presence of pyrite.     

Filter formulation and wavefield separation of cross-well seismic data

Multichannel filtering to obtain wavefield separation has been used in seismic processing for decades and has become an essential component in VSP and cross-well reflection imaging. The need for good multichannel wavefield separation filters is acute in borehole seismic imaging techniques such as VSP and cross-well reflection imaging, where strong interfering arrivals such as tube waves, shear conversions, multiples, direct arrivals and guided waves can overlap temporally with desired arrivals. We investigate the effects of preprocessing (alignment and equalization) on the quality of cross-well reflection imaging wavefield separation and we show that the choice of the multichannel filter and filter parameters is critical to the wavefield separation of cross-well data (median filters, f–k pie-slice filters, eigenvector filters). We show that spatial aliasing creates situations where the application of purely spatial filters (median filters) will create notches in the frequency spectrum of the desired reflection arrival. Eigenvector filters allow us to work past the limits of aliasing, but these kinds of filter are strongly dependent on the ratio of undesired to desired signal amplitude. On the basis of these observations, we developed a new type of multichannel filter that combined the best characteristics of spatial filters and eigenvector filters. We call this filter a ‘constrained eigenvector filter’. We use two real data sets of cross-well seismic experiments with small and large well spacing to evaluate the effects of these factors on the quality of cross-well wavefield separation. We apply median filters, f–k pie-slice filters and constrained eigenvector filters in multiple domains available for these data sets (common-source, common-receiver, common-offset and common-midpoint gathers). We show that the results of applying the constrained eigenvector filter to the entire cross-well data set are superior to both the spatial and standard eigenvector filter results.     

Preserved‐amplitude angle domain migration by shot‐receiver wavefield continuation

We present preserved‐amplitude downward continuation migration formulas in the aperture angle domain. Our approach is based on shot‐receiver wavefield continuation. Since source and receiver points are close to the image point, a local homogeneous reference velocity can be approximated after redatuming. We analyse this approach in the framework of linearized inversion of Kirchhoff and Born approximations. From our analysis, preserved‐amplitude Kirchhoff and Born inverse formulas can be derived for the 2D case. They involve slant stacks of filtered subsurface offset domain common image gathers followed by the application of the appropriate weighting factors. For the numerical implementation of these formulas, we develop an algorithm based on the true amplitude version of the one‐way paraxial approximation. Finally, we demonstrate the relevance of our approach with a set of applications on synthetic datasets and compare our results with those obtained on the Marmousi model by multi‐arrival ray‐based preserved‐amplitude migration. While results are similar, we observe that our results are less affected by artefacts.     

交错网格与旋转交错网格对VTI介质波场分离的影响分析

利用交错网格有限差分和旋转交错网格有限差分进行各向异性介质弹性波场数值模拟时, 质点振动速度分量与应力张量的网格节点定义方式均不相同, 从而对各向异性波场分离效果产生不同的影响. 针对这一问题, 本文以具有垂直对称轴的横向各向同性(VTI)介质的波场分离为例, 首先分析了两种网格的参数定义方式以及VTI介质波场的分离过程; 其次, 详细研究和分析了这两种网格的参数定义方式对各向异性介质波场分离的影响, 并依据波前面连续性以及波场分离效果等方面, 通过数值模拟实验对该影响进行分析验证. 结果表明, 旋转交错网格的参数定义方式更有利于进行各向异性介质波场数值模拟和波场分离.      

基于频率域峰值属性的河道砂体定量预测及应用(英文)

河道砂体是陆相含油气盆地最重要的储集类型之一,其边界识别和厚度定量预测是储层预测的热点难题。本文在总结现有方法技术的基础上,提出一种利用频率域峰值属性进行河道砂体边界识别和厚度定量预测的新方法。对典型河道薄砂体地震反射进行了正演模拟,构造了一种新的地震属性——峰值频率-振幅比,研究表明:峰值频率属性对地层厚度变化敏感,振幅属性对地层岩性变化敏感,两者比值突出河道砂体的边界,同时,借助峰值频率与薄层厚度间存在的定量关系进行薄砂体厚度计算。实际数据应用表明,地震峰值频率属性可以较好的刻画河道的平面展布特征;峰值频率-振幅比属性可以提高对河道砂体边界的识别能力;利用频率域地震属性进行砂体边界识别及厚度定量预测是可行的。     

Image‐domain wavefield tomography with extended common‐image‐point gathers

Waveform inversion is a velocity‐model‐building technique based on full waveforms as the input and seismic wavefields as the information carrier. Conventional waveform inversion is implemented in the data domain. However, similar techniques referred to as image‐domain wavefield tomography can be formulated in the image domain and use a seismic image as the input and seismic wavefields as the information carrier. The objective function for the image‐domain approach is designed to optimize the coherency of reflections in extended common‐image gathers. The function applies a penalty operator to the gathers, thus highlighting image inaccuracies arising from the velocity model error. Minimizing the objective function optimizes the model and improves the image quality. The gradient of the objective function is computed using the adjoint state method in a way similar to that in the analogous data‐domain implementation. We propose an image‐domain velocity‐model building method using extended common‐image‐point space‐ and time‐lag gathers constructed sparsely at reflections in the image. The gathers are effective in reconstructing the velocity model in complex geologic environments and can be used as an economical replacement for conventional common‐image gathers in wave‐equation tomography. A test on the Marmousi model illustrates successful updating of the velocity model using common‐image‐point gathers and resulting improved image quality.     

基于解析时间波场外推与波场分解的逆时偏移方法研究

双程波方程逆时深度偏移是复杂介质高精度成像的有效技术, 但其结果中通常包含成像方法引起的噪音和假象, 一般的滤波方法会破坏成像剖面上的振幅, 其中的假象也会给后续地质解释带来困扰.将波场进行方向分解然后实现入射波与反射波的相关成像能够有效地消除这类成像噪音, 并提高逆时偏移成像质量.波传播方向的分解通常在频率波数域实现, 它会占用大量的存储和计算资源, 不便于在沿时间外推的逆时深度偏移中应用.本文提出解析时间波场外推方法, 可以在时间外推的每个时间片上实现波传播方向的显式分解, 逆时深度偏移中利用分解后的炮检波场进行对应的相关运算, 实现成像噪音和成像信号的分离.在模型和实际数据上的测试表明, 相比于常规互相关逆时偏移成像结果, 本文方法能够有效地消除低频成像噪音和特殊地质构造导致的成像假象.     

线性同相轴波场分离的高分辨率τ-p变换法

基于最小二乘τ-p变换和τ-p域模型稀疏分布的假设，本文给出高分辨率τ-p变换的推导及其模型空间域的离散采样公式，同时给出了保振幅线性同相轴波场分离的算法流程.在求解本文给出的高分辨率τ-p正变换时，由于待求解的矩阵不具备最小平方法所具有的Toeplitz结构，故采用Cholesky分解法进行计算.本文模拟了井间地震和阵列声波测井中的Stoneley上下行波的分离算法过程，高分辨率正反τ-p变换且滤波所得结果显示本文算法误差小和保振幅的特点.对于在τ-p域距离很近或时间域同相轴近于水平的线性波场，高分辨率算法的聚焦作用使得所分离波场畸变小，体现本文算法精度高的优点.理论模型试算表明本文给出的高分辨率τ-p变换线性波场分离算法具有稳定性、精度高和保振幅的特点.     

Research Note: Internal/external seismic source wavefield separation and cancellation

A compact representation is obtained for the separation of a scalar wavefield on a closed surface into parts due to internal and external sources. The formula assumes that the total field and its gradients are known on the surface, as is the exact Green function of the medium. The derivation involves four rather straightforward applications of Green’s theorem or the representation theorem, though it is a remarkable result in that waves from either source that traverse the boundary many times are appropriately separated. The intermediate results at the four steps of the derivation also shed light on the possibility of acoustic shielding from unwanted sources without knowledge of the Green function for the medium.     

基于Bregman迭代的复杂地震波场稀疏域插值方法

在地震勘探中,野外施工条件等因素使观测系统很难记录到完整的地震波场,因此,资料处理中的地震数据插值是一个重要的问题。尤其在复杂构造条件下,缺失的叠前地震数据给后续高精度处理带来严重的影响。压缩感知理论源于解决图像采集问题,主要包含信号的稀疏表征以及数学组合优化问题的求解,它为地震数据插值问题的求解提供了有效的解决方案。在应用压缩感知求解复杂地震波场的插值问题中,如何最佳化表征复杂地震波场以及快速准确的迭代算法是该理论应用的关键问题。Seislet变换是一个特殊针对地震波场表征的稀疏多尺度变换,该方法能有效地压缩地震波同相轴。同时,Bregman迭代算法在以稀疏表征为核心的压缩感知理论中,是一种有效的求解算法,通过选取适当的阈值参数,能够开发地震波动力学预测理论、图像处理变换方法和压缩感知反演算法相结合的地震数据插值方法。本文将地震数据插值问题纳入约束最优化问题,选取能够有效压缩复杂地震波场的OC-seislet稀疏变换,应用Bregman迭代方法求解压缩感知理论框架下的混合范数反问题,提出了Bregman迭代方法中固定阈值选取的H曲线方法,实现地震波场的快速、准确重建。理论模型和实际数据的处理结果验证了基于H曲线准则的Bregman迭代稀疏域插值方法可以有效地恢复复杂波场的缺失信息。     

Velocity–depth and time–depth relationships for a decompacted uplifted unit

A linear instantaneous velocity model is used to describe the velocity variations in an uplifted unit that has been partly decompacted as a result of the reduction in overburden that often accompanies uplift. The model results in a series of equations for deriving values for the function parameters in the velocity–depth and the time– depth domains and for carrying out time‐to‐depth conversions. The formulation uses the base of the unit as a reference level to generate the reference datum from a combination of the depth of the base of the unit and a parameter that represents the decompaction factor.     

散射系数矩阵的计算及其应用

利用Gabor-Daubechies（G?-D）小波束域波场分解和传播在空间和方向上的双重局域性,提出了基于G?-D小波束域叠前深度偏移进行角度域成像和计算局部散射系数矩阵的方法. 以简单分层模型为例,对不同探测系统的局部散射系数矩阵分布特征进行分析. 分析结果表明,在一定的探测系统几何布局下,由本文方法得到的局部散射系数矩阵能够较真实地反映局部结构的散射(或反射)特性. 通过局部散射系数矩阵进一步外推具有不同速度反差的水平界面随角度变化的反射系数,并估计界面的空间位置和倾角等说明局部散射系数矩阵的潜在应用.     

Dynamic response analysis in the frequency domain

Analysis in the frequency domain using discrete Fourier transforms is an efficient means of calculating the dynamic response of linear systems. In fact, for systems with frequency dependent parameters and also in those cases where the complex frequency response functions are more easily determined, frequency domain analysis may be the only effective means of determining the dynamic response. The use of discrete transforms along with finite summation requires that the forcing function and the impulse function be converted into periodic forms. This modification may introduce unacceptable errors in the results of analysis, unless appropriate steps are taken to avoid or minimize the effect of aliasing or overlapping. For single-degree-of-freedom systems, procedures that will eliminate the effect of aliasing have already been developed. However, problems related to frequency domain analysis still exist for multi-degree-of-freedom systems with non-proportional damping, in analysis through substructuring and in those cases where a continuum solution is involved. A new procedure which addresses these problems and is applicable to both single- and multi-degree-of-freedom system as well as to analysis through substructuring is presented here.     

基于广义Curvelet变换和相空间特征聚类的VSP波场分离方法

垂直地震剖面(Vertical Seismic Profiling, VSP)蕴含大量的地层地质信息, 是连接地震反射数据与测井资料之间的桥梁.VSP数据分辨率高, 资料丰富, 包括上行波、下行波、转换波等多种类型的波场, 常被用于地层反射界面标定、介质衰减参数反演等.其中, 如何有效地分离波场是将VSP资料应用于油气勘探的关键之一.广义Curvelet变换是一种具有多尺度、多角度的相空间变换, 能够进行相空间域波场特性的提取.本文聚焦于VSP波场的上下行波分离问题, 提出一种基于广义Curvelet变换和相空间特征聚类方法.首先, 采用广义Curvelet变换对VSP波场进行特征提取, 构造融入相空间角度信息的特征数据集; 然后, 利用K-means聚类算法在相空间对波场特征进行分类; 最后, 对分类结果进行反变换, 完成VSP上下行波波场的自适应分离.为了验证方法的有效性, 本文将所提出的方法用于合成数据和实际VSP数据的波场分离, 并与常用的基于F-K变换的波场分离方法进行对比.处理结果表明, 本文方法角度分辨率高、抗噪能力强, 波场分离结果的保真保幅性好, 这为后续的成像与地层的特征分析提供了重要基础资料.

     

地震计阻尼和自振频率的频域测定

利用地震观测系统中的阶跃标定响应记录 ,在频域基于阶跃信号的幅频特性确定地震计的阻尼常数和自振周期。介绍了该方法的原理 ,采用非线性最小二乘拟合方法和具体的实际应用