基于Poynting矢量行波分离的最小二乘逆时偏移

因为在逆时偏移中基于双程波动方程构建震源波场和检波器波场,所以在波场延拓过程中地震波遇到波阻抗界面时,背向发育的反射波会与正常传播的波场互相关产生较强振幅的低频噪声.这一特点使得以逆时偏移为基础的最小二乘偏移方法在梯度计算时同样存在着低频噪声的干扰,从而导致反演收敛的速度减慢.考虑到计算量和存储成本的因素,本文借助Po...     

基于Poynting矢量的归一化波场分离互相关逆时偏移成像条件(英文)

叠前逆时偏移是目前精度较高的成像方法,然而严重的低频噪声降低了逆时偏移的构造成像精度,偏移噪声的压制是逆时偏移必需要考虑的问题。分析了低频噪声的产生机理,并根据声波方程Poynting矢量的方向指示地震波场的传播方向的原理,分离出上、下、左、右行波,该方法计算量和存储量都远远小于常用的二维傅里叶变换分离方法,进而提出归一化的波场分离互相关成像条件,以压制逆时偏移低频噪声,提高成像精度。实现了Marmousi模型的试算,表明在波场延拓过程中利用Poynting矢量能够较好的分离上、下、左、右行波,与常规方法、拉普拉斯滤波、二维傅里叶变换波场分离的成像结果对比表明,成像时使用归一化的波场分离互相关成像条件能更好的压制偏移噪声,得到精度更高的逆时偏移成像结果。     

Reverse-time migration using multidirectional wavefield decomposition method

Reverse-time migration has attracted more and more attention owing to the advantages of high imaging accuracy, no dip restriction, and adaptation to complex velocity models. Cross-correlation imaging method is typically used in conventional reverse-time migration that produces images with strong low-frequency noise. Wavefield decomposition imaging can suppress such noise; however, some residual noise persists in the imaging results. We propose a 2D multidirectional wavefield decomposition method based on the traditional wavefield decomposition method. First, source wavefields and receiver wavefields are separated into eight subwavefields, respectively. Second, cross-correlation imaging is applied to selected subwavefields to produce subimages. Finally, the subimages are stacked to generate the final image. Numerical examples suggest that the proposed method can eliminate the low-frequency noise effectively and produce high-quality imaging profiles.     

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

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

基于自适应优化有限差分方法的全波VSP逆时偏移

与地面地震资料相比,VSP资料具有分辨率高、环境噪声小及能更好地反映井旁信息等优点.常规VSP偏移主要对上行反射波进行成像,存在照明度低、成像范围受限等问题.为了增加照明度、拓宽成像范围、提高成像精度,本文采用直达波除外的所有声波波场数据(全波),包括一次反射波、多次反射波等进行叠前逆时偏移成像.针对逆时偏移中的四个关键问题,即波场延拓、吸收边界条件、成像条件及低频噪声的压制,本文分别采用自适应变空间差分算子长度的优化有限差分方法(自适应优化有限差分方法)求解二维声波波动方程以实现高精度、高效率的波场延拓,采用混合吸收边界条件压制因计算区域有限所引起的人工边界反射,采用震源归一化零延迟互相关成像条件进行成像,采用拉普拉斯滤波方法压制逆时偏移中产生的低频噪声.本文对VSP模型数据的逆时偏移成像进行了分析,结果表明:自适应优化有限差分方法比传统有限差分方法具有更高的模拟精度与计算效率,适用于VSP逆时偏移成像;全波场VSP逆时偏移成像比上行波VSP逆时偏移的成像范围大、成像效果好;相对于反褶积成像条件,震源归一化零延迟互相关成像条件具有稳定性好、计算效率高等优点.将本文方法应用于某实际VSP资料的逆时偏移成像,进一步验证了本文方法的正确性和有效性.     

Reverse time migration of ground-penetrating radar with full wavefield decomposition based on the Hilbert transform

The conventional reverse time migration of ground-penetrating radar data is implemented with the two-way wave equation. The cross-correlation result contains low-frequency noise and false images caused by improper wave paths. To eliminate low-frequency noise and improve the quality of the migration image, we propose to separate the left-up-going, left-down-going, right-up-going and right-down-going wavefield components in the forward- and backward-propagated wavefields based on the Hilbert transform. By applying the reverse time migration of ground-penetrating radar data with full wavefield decomposition based on the Hilbert transform, we obtain the reverse time migration images of different wavefield components and combine correct imaging conditions to generate complete migration images. The proposed method is tested on the synthetic ground-penetrating radar data of a tilt-interface model and a complex model. The migration results show that the imaging condition of different wavefield components can highlight the desired structures. We further discuss the reasons for incomplete images by reverse time migration with partial wavefields. Compared with the conventional reverse time migration methods for ground-penetrating radar data, low-frequency noise can be eliminated in images generated by the reverse time migration method with full wavefield decomposition based on the Hilbert transform.     

相对保幅的角度域VSP逆时偏移(英文)

本文介绍了一种改进的角度域VSP逆时偏移方法。对VSP逆时偏移中的逆推公式进行了改进,为方便数值计算出相对保幅的角度域共成像点道集(ADCIGs)。此外VSP记录到的波场信息丰富,包括上行波场、下行波场和直达波场等,本文分析了这些波场的响应特征,发现直达波和下行波在角度域共成像点道集(ADCIGs)上都产生了成像噪音,直达波产生的噪音尤为严重。把该方法用于我国西部地区实际观测的VSP资料,不仅获得相对保幅角度域共成像点道集(ADCIGs),而且压制了成像噪音。通过数值模型试算,实际资料的应用验证了该方法的实用性与有效性,从而为VSP偏移速度分析、VSP AVA/AVO分析和反演等提供可靠的基础资料。     

一种TI介质纯qP波正演方法及其在逆时偏移中的应用

基于Tsvankin提出的精确频散关系,利用近似展开的方法,推导出解耦合的TTI介质纯qP波近似方程,并将方程中的偏微分算子分解成一个laplace算子和一个标量算子,用于代表qP波的精确传播方向,构建时间域二阶纯qP波方程.此推导过程无需设置横波速度为零,能够更加精确地描述qP波的运动学特征.这个方程相比于求解波数域二阶解耦qP波方程,计算效率高,存储需求小;相比于基于Alkhalifah频散关系推导的时间域二阶纯qP波方程,假象干扰压制好,数值误差小,更具一般性.但此方法求解波矢量时采用波场梯度一阶渐近近似,会造成垂直于对称轴方向的波场振幅不准确.为了较正振幅,将椭圆分解方法应用于此方程中,构建纯qP波椭圆分解方程,使得振幅更加均衡,并与Xu等提出的方程比较分析,应用本文构建的纯qP波椭圆分解方程得到的波场振幅值更加准确.本文首先选取了均匀TI介质模型进行了qP波正演模拟,并抽取波场单道波形进行振幅分析,验证了本文构建的纯qP波方程和纯qP波椭圆分解方程的正确性及有效性;然后选取BP TTI模型进行了qP波正演模拟,将其qP波正演结果和均匀TI介质模型振幅分析结果相结合,突出了本文构建的纯qP波椭圆分解方程的优势及适应性;最后选取逆冲模型和BPTTI模型,应用本文构建的纯qP波椭圆分解方程对其进行逆时偏移成像,验证了本文构建的纯qP波椭圆分解方程在逆时偏移中的可行性和适用性.     

基于波场传播方向的波动方程偏移有效去噪方法(英文)

本文研究分析了双程波波动方程偏移成像中广泛存在的三种主要噪声,特别是针对过去研究中没能很好解决的存在于高速盐丘悬垂边界附近的射线状噪声,提出了基于优化成像条件的有效去噪方法。射线状噪声主要来自于震源一侧波场的下行透射波分量和接收阵列一侧波场的上行散射波互相关成像。这一部分能量具有较强的互相关性,但并不携带真实的反射面信息。它广泛存在叠前偏移成像中,与信号在强度上同量级。多数情况下偏移成像中的相关噪声由方向性传播的波场能量产生。利用波场梯度得到的波场传播角度,可以分离出噪声对应的波场能量,并在成像条件中减去。采用这一方法可以有效地去除多种噪声,包括直达波噪声、散射波噪声和射线状噪声。该去噪方法不依赖波场外推算子,在需要时可以方便地运用到几乎所有的波动方程偏移中去。并且该去噪方法针对噪声的物理根源,对信号的损害很小。对去噪后的偏移成像结果额外地进行波数域滤波处理,可以进一步提高叠加图像的质量。这一去噪方法在超广角单程波偏移成像中取得良好效果,我们同时期待其在其他双程波波动方程偏移特别是逆时偏移(RTM)中的成功运用。     

利用微震数据对震源下方构造直接成像

现有的微震记录直接成像方法是将微震记录既当作入射记录,也当作散射记录,从而实现偏移成像.但此方法并不能突出透射波所携带的来自震源下方的深层散射波信息.本文在假设已知微震位置与子波的前提下,提出了对微震下方构造进行逆时偏移的成像方法.该方法类似于常规的逆时偏移,只是震源位置在地下.这使得在成像时,地下更深部的入射波场相比震源在地表时会更为精确,因此能够获得更加准确的成像结果.该方法会给成像结果带来一种尾波高频干扰：地下的震源发出的上行波与上方介质作用后,所产生的多级散射波会干扰反传波场.对此,在成像过程中,对入射场和散射场都进行左右行波分离,以压制该噪声.而在子波信息未知,无法重构入射场时,使用了激发时间成像条件,也能够实现同等效果的偏移成像,且不会出现尾波高频干扰.利用数值实验验证了本文方法的有效性.     

用最大似然法进行波场分解和震相识别

本文以介质中弹性波场的传播特征和多元统计分析为基础,从数学上进行严格的推导,得到了模型与资料拟合的最大似然判据,以检验地震波列中是否存在某种特殊波型的能量。利用概率滤波可得分解后的地震图,从图中可以直接得到震相到时。用此方法处理P波资料,效果较好。对于剪切波资料,有时处理结果不十分理想。这可能是由于径向分量和横向分量之间的线性相关性造成的。对于因介质各向异性而形成的分裂剪切波,用此方法得到的S波到时为快S波到时。此方法对于单台三分向记录资料的处理是十分有效的。     

基于CLG光学流和波场分解的逆时偏移角度道集提取方法研究

角度域共成像点道集是衔接叠前地震数据与储层特征的重要桥梁,对地震偏移成像与储层描述具有重要意义.与克希霍夫偏移和单向波动方程偏移相比,逆时偏移是复杂地区最精确的成像方法.高效稳健地生成逆时偏移角度道集目前仍然是一个挑战.本文主要讨论如何采用光学流方法高效、高质量地提取角度道集.在逆时偏移波场外推过程中,光学流方法可以估计波场传播方向.其中Lucas-Kanade(LK)和Horn-Schunck(HS)方法是光学流方法中两种典型的方法.LK光学流方法是一种局部方法,该方法依赖于局部点的梯度值,但是容易出现奇异现象,HS光学流方法属于全局方法,波场方向估计依赖于整个波场,易受噪声影响,对异常值比较敏感,导致整体波场方向计算精度不高.本文提出采用局部和整体结合(Combining Local and Global,CLG)的光学流方法估计波场传播方向.该方法可以有效地提高波场方向的精度,并且简单高效,便于并行处理.对比HS光学流方法,CLG光学流方法几乎不增加额外的计算量.另外,为了弱化光学流方法无法处理波前重叠问题,本文利用解析波场和方向滤波对波场进行方向分解,仅需波场的空间傅里叶变换即可实现任意波场方向分解,将分解后的波场分别估计波场反向,提取成像结果.进一步地,在估计反射张角和方位角时,本文提出有效的归一化方法和改进的最小二乘除法,提高角度估计的精度和稳定性.最后,理论和实际资料例证了本文提出方法的有效性.     

一种起伏地表条件下的照明补偿方法

起伏的地表条件限制了采集孔径范围并造成深层地震照明不足,为改善该类地区的成像质量,本文提出了一种起伏地表条件下的照明补偿方法.首先,基于小波束波场延拓算子和逐步累加的外推方法在波场延拓过程中解决起伏地表面的影响,并引入空间滤波函数压制虚拟层内的偏移噪音;其次,利用局部指数标架对上、下行波场分解,得到局部角度域成像和照明补偿因子.再次,利用计算出的成像值和照明补偿因子,在局部倾角域完成照明补偿.SEG起伏地表模型测试证明了本方法的有效性,深层构造照明度明显加强,不同角度成像振幅更加均衡,该技术为提高起伏地表地区的成像品质提供了新的手段.     

Imaging artefacts of artificial diving waves in reverse time migration: cause analysis in the angle domain and an effective removal strategy

In areas with strong velocity gradients, traditional reverse time migration based on cross-correlation imaging condition not only produces low-frequency noise but also generates diving wave artefacts. The artefacts caused by diving waves have no typical low-frequency characteristics and cannot be eliminated by simple high-pass filtering approaches. We apply the wave-field decomposition imaging condition to analyse the causes of false images in reverse time migration by decomposing the full wave-field into up-going and down-going components in the angle domain. We find that artificial diving wave imaging artefacts, which are generated by the cross-correlation between the up-going source and down-going receiver wave-fields in areas with strong velocity gradients, arise at large angles. We propose an efficient strategy by means of the wavelength-dependent smoothing operator to eliminate artefacts from artificial diving waves in reverse time migration. Specifically, the proposed method provides more reasonable down-going wave-fields in areas with sharp velocity constructs by considering the factor of varying seismic wavelengths during wave propagation, and the artificial components of diving waves are eliminated in a straightforward manner. Meanwhile, the other wave-field components that contribute to true subsurface images are minimally affected. Benefiting from a smoothed velocity, the proposed method can be adapted to the traditional reverse time migration imaging frame, which reveals significant implementation potential for the seismic exploration industry. A salt model is designed and included to demonstrate the effectiveness of our approach.     

Ground roll wave suppression based on wavelet frequency division and radial trace transform

Ground roll waves interfere with seismic data. The suppression of ground roll waves based on the division of wavelet frequencies considers the low-frequency characteristics of ground roll waves. However, this method will not be effective when the ground roll wave and the effective signal have the same frequency bands because of overlapping. The radial trace transform (RTT) considers the apparent velocity difference between the effective signal and the ground roll wave to suppress the latter, but affects the low-frequency components of the former. This study proposes a ground roll wave suppression method by combining the wavelet frequency division and the RTT based on the difference between the ground roll wave velocity and the effective signal and their energy difference in the wavelet domain, thus making full use of the advantages of both methods. First, we decompose the seismic data into different frequency bands through wavelet transform. Second, the RTT and low-cut filtering are applied to the low-frequency band, where the ground roll waves are appearing. Third, we reconstruct the seismic record without ground roll waves by using the inverse RTT and the remaining frequency bands. The proposed method not only improves the ground roll wave suppression, but also protects the signal integrity. The numerical simulation and real seismic data processing results suggest that the proposed method has a strong ability to denoise while preserving the amplitude.     

Ore-body imaging by crosswell seismic waveform inversion: A case study from Kambalda,Western Australia

Seismic P waveforms from a hardrock crosswell experiment were inverted for the velocity distribution between boreholes. The inversion was carried out by an iterative source estimation scheme in the frequency domain. Six frequencies in the range 900–2000 Hz were employed to yield a tomogram which closely matches the known geology, and clearly discloses the nickel sulphide mineralization. The waveform inversion is superior to traveltime tomography, which was used to obtain the starting model. Pre-processing, to isolate the direct P arrivals, is critical to the success of the method, because noise in the form of tube wave-to-body-wave conversions in the source hole, and body wave-to-tube wave conversions in the receiver hole complicate the seismograms and cannot be satisfactorily accommodated in the wavefield modelling. Such noise also precludes a successful migration. Even though the very fine detailed velocity structure cannot be fully recovered using just the first arrival full waveform, it remains a viable and useful approach.     

Arbitrary-order Taylor series expansion-based viscoacoustic wavefield simulation in 3D vertical transversely isotropic media

Taking the anisotropy of velocity and attenuation into account, we investigate the wavefield simulation of viscoacoustic waves in 3D vertical transversely isotropic attenuating media. The viscoacoustic wave equations with the decoupled amplitude attenuation and phase dispersion are derived from the fractional Laplacian operator and using the acoustic approximation. With respect to the spatially variable fractional Laplacian operator in the formulation, we develop an effective algorithm to realize the viscoacoustic wavefield extrapolation by using the arbitrary-order Taylor series expansion. Based on the approximation, the mixed-domain fractional Laplacian operators are decoupled from the wavenumbers and fractional orders. Thus, the viscoacoustic wave propagation can be conveniently implemented by using a generalized pseudospectral method. In addition, we perform the accuracy and efficiency analyses among first-, second- and third-order Taylor series expansion pseudospectral methods with different quality factors. Considering both the accuracy and computational cost, the second-order Taylor series expansion pseudospectral method can generally satisfy the requirements for most attenuating media. Numerical modelling examples not only illustrate that our decoupled viscoacoustic wave equations can effectively describe the attenuating property of the medium, but also demonstrate the accuracy and the high robustness of our proposed schemes.     

高保真VSP波场分离方法研究

三维VSP资料中,各种不同类型的波混杂一起形成复杂的波场.因此,波场分离是三维VSP数据处理关键的第一步.从不同波场的偏振方向和传播方向之差异着手,提出了一种高保真的VSP波场分离方法.首先通过射线追踪和偏振滤波的结合,把复杂波场(分解为简单波场;然后根据简单波场中不同波的传播方向截然相反的特点,进行方向滤波,达到波场分离的目的.实际数据处理表明,与常规波场分离方法相比,本方法大大降低了混波作用以及由此而生的波形畸变.     

面炮成像、控制照明与AVA道集

基于波场延拓的叠前深度偏移是实现复杂构造地质体成像的最可靠方法,但存在着计算量大、对观测系统适应性差等缺点.面炮偏移是波动方程实现精确叠前成像的另一类方法,具有较高的计算效率,不存在偏移孔径问题,而且可以通过控制照明方法,解决平面波在目标区域的能量补偿问题.本文采用面炮成像技术进行叠前深度偏移,通过对面炮震源下行波场的质量控制和优选射线个数和范围,以达到最佳的成像效果.采用控制照明技术,较大地提高了目标地层的成像精度.与此同时,得到振幅随入射角变化(AVA)道集,有利于叠前振幅解释和储层岩性预测.数据实验表明面炮成像技术是一种快速有效的方法,其成像精度与单平方根算子的共炮点道集偏移和双平方根算子的共中心点道集偏移相当,但在计算速度上要快得多,而且易于并行计算.