纵横波独立成像的弹性波逆时偏移方法

本文采用逆时递推,在空间域完成纵横波分离,同时实现了位移场成像和纵、横波的独立成像。模拟数据与实际资料试算表明,对纵波和横波分别成像不仅有理论意义,而且有明显的应用价值。文中还就偏移速度对成像结果的影响和实际资料的预处理问题进行了讨论。     

纵横波独立成像的弹性波逆时偏移方法

本文采用逆时递推,在空间域完成纵横波分离,同时实现了位移场成像和纵、横波的独立成像。模拟数据与实际资料试算表明,对纵波和横波分别成像不仅有理论意义,而且有明显的应用价值。文中还就偏移速度对成像结果的影响和实际资料的预处理问题进行了讨论。     

基于波场分离的弹性波逆时偏移

尽管叠前逆时偏移成像精度高,但仅针对单一纵波的成像也可能形成地下介质成像盲区,由于基于弹性波方程的逆时偏移成像可形成多波模式的成像数据,因此弹性波逆时偏移成像可提供更为丰富的地下构造信息.本文依据各向同性介质的一阶速度-应力方程组构建震源和检波点矢量波场,再利用Helmholtz分解提取纯纵波和纯横波波场,使用震源归一化的互相关成像条件获得纯波成像,避免了直接使用坐标分量成像而引起的纵横波串扰问题.针对转换波成像的极性反转问题,文中提出一种共炮域极性校正方法.为有效节约存储成本,也提出一种适用于弹性波逆时偏移的震源波场逆时重建方法,在震源波场正传过程中,仅保存PML边界内若干层的速度分量波场,进而逆时重建出所有分量的震源波场.本文分别对地堑模型和Marmousi2模型进行了弹性波逆时偏移成像测试,结果表明：所提出的共炮域极性校正方法正确有效,基于波场分离的弹性波逆时偏移成像的纯波数据能够对复杂地下构造准确成像.

     

纵横波独立成像的弹性波逆时偏移方法

本文采用逆时递推,在空间域完成纵横波分离,同时实现了位移场成像和纵、横波的独立成像。模拟数据与实际资料试算表明,对纵波和横波分别成像不仅有理论意义,而且有明显的应用价值。文中还就偏移速度对成像结果的影响和实际资料的预处理问题进行了讨论。     

声反射成像测井弹性波逆时偏移方法对比研究

声反射成像测井中的声波逆时偏移算法在处理单分量偶极横波反射资料时可获得高精度成像结果, 而面对单极声源中的多波多分量弹性波反射资料时, 由于其包含PP、PS、SP和SS等多种波型分量, 基于声波方程的逆时偏移处理结果中出现单反射界面对应多同相轴的偏移假象.首先在声反射成像测井中实现了纵横波耦合弹性波逆时偏移和基于Helmholtz分解的弹性波逆时偏移, 在此基础上首次提出和实现了基于一阶速度-应力纵横波解耦方程的弹性波逆时偏移, 通过纵横波解耦结果和成像结果对三种声反射成像测井弹性波逆时偏移进行对比研究, 数值算例表明基于纵横波耦合的弹性波逆时偏移无法处理转换波(PS或SP)成像问题; Helmholtz分解虽可实现纵横波解耦但同时引起反射波振幅和相位改变, 进而导致包含S波的成像分量出现极性反转问题.基于一阶速度-应力纵横波解耦方程的弹性波逆时偏移可有效处理多波多分量反射资料偏移成像, 消除转换分量成像噪声, 解决S波分量成像极性反转, 为声反射成像测井弹性波逆时偏移算法奠定了有力的理论基础.

     

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

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

解耦纵横波反射波走时反演

基于常规弹性波动方程的反射波走时反演结合走时和反射波信息可以有效的摄取模型参数中的低波数成分,然而纵横波之间的耦合效应以及纵横波速度对波场的敏感性差异,导致反演的非线性问题增强.为此本文研究了基于解耦波动方程的反射波走时反演,并提出改进的时移互相关目标函数,分别隐式计入射波场快照与反传波场快照的时移量,很大程度的降低了纵波、横波之间的耦合关系,并提高纵横波速度低波数信息的反演质量.最后模型测试证明了本文方法的正确性.

     

三维各向异性TI介质中的P/S波快速解耦技术及在弹性波逆时偏移中的应用

随着多分量采集技术的发展,弹性波逆时偏移技术在三维各向异性介质复杂地质构造成像中得到了广泛的应用.然而耦合的P波场和S波场,会在传播过程中产生串扰噪声,降低弹性波逆时偏移的成像精度.为了解决这一问题,本研究针对具有倾斜各向异性对称轴的三维横向各向同性（Transverse Isotropy,TI）介质,提出了一种矢量弹性波场快速解耦方法,可以有效提高偏移剖面的成像质量.该方法首先通过坐标转换,将观测系统坐标系的垂直轴旋转到TI介质的对称轴方向,在新坐标系下,根据具有垂直对称轴的三维横向各向同性（Vertical Transverse Isotropy,VTI）介质中的分解算子,推导出三维TI介质解耦算子表达式.接着引入一种在空间域快速计算分解波场的方法,来实现空间域矢量P波场和S波场分离,极大地提高了计算效率.最后,通过点积成像条件,将提出的P/S波分解方法引入到三维TI介质弹性波逆时偏移中,得到高精度的PP和PS成像.与以往的波场分解方法相比,本文方法具有数值稳定和计算效率高的特点.数值算例表明,应用上述三维TI分解算子得到的偏移剖面有效压制了噪声,提高了成像质量.

     

基于波场分离理论的逆时偏移成像条件研究及应用

逆时偏移成像方法与其它算法的成像方法相比,由于算法不受地层倾角的限制,因此针对陡倾角构造和复杂地质模型地区的成像有着无可比拟的优越性,特别是在盐丘发育的地区优势更加明显.实际应用过程中,逆时偏移通常是应用互相关成像条件构建成像,然而,这种成像条件会产生低频、强振幅的噪音,如果不能彻底消除这些噪音,那它就会严重影响甚至淹没有效信号,特别是在浅层的强反射界面上方往往会出现大量的低频噪音,从而使得浅层构造基本无法识别.基于传统互相关成像条件的这一局限,本文提供了一个新的成像条件,它能在有效成像的同时消除这些低频、强振幅的干扰噪音;其具体实现思路就是将炮点及检波点波场分离成它们的单程波传播分量,然后采用互相关成像条件对分离后的波场进行成像,从而实现逆时偏移成像.通过对模型数据和实际野外数据的实测,证明了该方法的有效性和适用性.     

基于改进余弦组合窗的解耦方程弹性波逆时偏移

弹性波矢量波场逆时偏移可以综合利用纵横波场信息,对地下空间进行清晰成像,且对于成像介质没有角度限制,可以对复杂构造进行更清晰的成像.而弹性波逆时偏移中最重要的就是求解波动方程的算法,其直接影响成像的精度以及效率.本文引进电力系统谐波分析中常用的余弦组合窗函数,并通过一种新的优化算法得到了改进的余弦组合窗函数从而得到优化后的有限差分算子.并将此算子应用于解耦方程的矢量波场分离算法从而提高了成像精度.数值测试表明基于新算法的逆时偏移的成像精度和清晰度得到了明显的提高.

     

Poynting and polarization vectors based wavefield decomposition and their application on elastic reverse time migration in 2D transversely isotropic media

With the progress in computational power and seismic acquisition, elastic reverse time migration is becoming increasingly feasible and helpful in characterizing the physical properties of subsurface structures. To achieve high-resolution seismic imaging using elastic reverse time migration, it is necessary to separate the compressional (P-wave) and shear (S-wave) waves for both isotropic and anisotropic media. In elastic isotropic media, the conventional method for wave-mode separation is to use the divergence and curl operators. However, in anisotropic media, the polarization direction of P waves is not exactly parallel to the direction of wave propagation. Also, the polarization direction of S-waves is not totally perpendicular to the direction of wave propagation. For this reason, the conventional divergence and curl operators show poor performance in anisotropic media. Moreover, conventional methods only perform well in the space domain of regular grids, and they are not suitable for elastic numerical simulation algorithms based on non-regular grids. Besides, these methods distort the original wavefield by taking spatial derivatives. In this case, a new anisotropic wave-mode separation scheme is developed using Poynting vectors. This scheme can be performed in the angle domain by constructing the relationship between group and polarization angles of different wave modes. Also, it is performed pointwise, independent of adjacent space points, suitable for parallel computing. Moreover, there is no need to correct the changes in phase and amplitude caused by the derivative operators. By using this scheme, the anisotropic elastic reverse time migration is more efficiently performed on the unstructured mesh. The effectiveness of our scheme is verified by several numerical examples.     

2D and 3D amplitude-preserving elastic reverse time migration based on the vector-decomposed P- and S-wave records

The elastic reverse time migration approach based on the vector-wavefield decomposition generally uses the scalar product imaging condition to image the multicomponent seismic data. However, the resulting images contain the crosstalk artefacts and the polarity reversal problems, which are caused by the nonphysical wave modes and the angle-dependent reduction of image amplitudes, respectively. To overcome these two problems, we develop an amplitude-preserving elastic reverse time migration approach based on the vector-decomposed P- and S-wave seismic records. This approach includes two key points. The first is that we employ the vector-decomposed P- and S-wave multicomponent records to independently reconstruct the PP and PS reflection images to mitigate the crosstalk artefacts. The second is that we propose two schemes in addressing the issue of polarity reversal problem in the conventional PP image. One solution is to adopt the angle-dependent equation. Another one is to reconstruct an amplitude-preserving PP image with a separated scalar P-wave particle velocity, which has a clear physical meaning. Numerical examples using two-dimensional and three-dimensional models demonstrate that the proposed elastic reverse time migration approach can provide the images with better amplitude-preserving performance and fewer crosstalk artefacts, compared with the conventional elastic reverse time migration approach based on the scalar product imaging condition.     

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

因为在逆时偏移中基于双程波动方程构建震源波场和检波器波场,所以在波场延拓过程中地震波遇到波阻抗界面时,背向发育的反射波会与正常传播的波场互相关产生较强振幅的低频噪声.这一特点使得以逆时偏移为基础的最小二乘偏移方法在梯度计算时同样存在着低频噪声的干扰,从而导致反演收敛的速度减慢.考虑到计算量和存储成本的因素,本文借助Poynting矢量良好的方向指示性实现波场的上下行波分离,并在早期迭代的梯度计算中只保留震源波场和检波器波场沿不同垂直方向传播的组分之间的互相关,有效避免了成像噪声的干扰,提高了算法收敛的速率.数值算例验证了方案的有效性.

     

二维TI介质非结构网格弹性波矢量逆时偏移

波场延拓得到的多分量波场中既包含纵波信息也包含横波信息,能否在全波场中实现纵横波的分离对各向同性和各向异性逆时偏移都有非常重要的意义.传统的散度旋度分离只适应于各向同性介质而对各向异性介质却无效.在非规则、非结构网格的弹性波数值模拟方法的基础上,发展了一种适应于各向异性介质的波场分离方法.该方法通过求解Christoffel方程,得到相角和极化角的关系,再利用群角和相角的关系,直接得到群角和极化角的关系.该方法与现存的各向异性波场分离相比,获得的计算效率改进更显著,而且存储量小.用简单各向异性模型和SEG各向异性Hess模型进行测试,都得到了较好的效果,证明了本文方法的有效性.

     

起伏地表采集数据的三维直接叠前时间偏移方法

提出一种可对起伏地表采集的三维地震资料直接进行偏移成像的叠前时间偏移方法和流程.它用两个等效速度描述近地表和上覆层对地震波传播的影响,可对炮、检点不在同一水平面的三维地震资料直接进行叠前时间偏移处理.该方法不对近地表地震波传播做垂直出、入射假定,因此可适应高速层出露等不存在明显低、降速带情况.描述近地表和上覆层的两个等效速度参数可依据偏移道集的同相轴是否平直来确定,避免了确定近地表速度的困难;而对已知近地表速度的情况,则可进一步修正近地表速度,获得更好的成像效果.用三维起伏地表的理论数据和中国东部某工区实际数据验证了所发展方法和处理流程的有效性和实用性.     

基于面波频散的三维横波速度方位各向异性层析成像方法

地震各向异性是反映地球内部介质特性的重要指针之一。常用的横波分裂法和二维面波方位各向异性层析成像方法很难准确反映各向异性随深度的变化。将与周期相关的区域化面波方位各向异性转换成与深度相关的一维横波速度方位各向异性可以弥补深度信息不足的缺陷。现有三维横波速度各向异性研究多是通过两步方法来实现的,即逐个周期二维面波方位各向异性层析成像以及逐个格点一维横波速度方位各向异性反演。这种分步反演的方式既不利于三维先验约束的引入,也不利于利用原始观测拟合误差对三维模型进行直接评估。因此本文开发了基于面波频散曲线的三维横波速度方位各向异性层析成像方法,并编制了相关正演和反演程序。为了检测方法和程序的有效性,我们对规律分布的三维检测板模型进行了模拟测试。测试结果显示:该方法可以很好地恢复各向同性波速异常、各向异性相对强度和快波方向等三维结构信息;而且反演模型相对于参考模型明显改善了对观测数据的拟合,降低了对观测数据的均方根误差。但对各向同性理论模型进行各向异性反演时,在波速均匀区可产生小于0.5%的假各向异性幅值,在波速非均匀区该假的各向异性幅值会更大,浅部可达3.5%。因此在实际应用中需要谨慎解释(浅...     

Topographic elastic least-squares reverse time migration based on vector P- and S-wave equations in the curvilinear coordinates

Elastic least-squares reverse time migration has been applied to multi-component seismic data to obtain high-quality images. However, the final images may suffer from artefacts caused by P- and S-wave crosstalk and severe spurious diffractions caused by complex topographic surface conditions. To suppress these crosstalk artefacts and spurious diffractions, we have developed a topographic separated-wavefield elastic least-squares reverse time migration algorithm. In this method, we apply P- and S-wave separated elastic velocity–stress wave equations in the curvilinear coordinates to derive demigration equations and gradient formulas with respect to P- and S-velocity. For the implementation of topographic separated-wavefield elastic least-squares reverse time migration, the wavefields, gradient directions and step lengths are all calculated in the curvilinear coordinates. Numerical experiments conducted with the two-component data synthetized by a three-topographic-layer with anomalies model and the Canadian Foothills model are considered to verify our method. The results reveal that compared with the conventional method, our method promises imaging results with higher resolution and has a faster residual convergence speed. Finally, we carry out numerical examples on noisy data, imperfect migration velocity and inaccurate surface elevation to analyse its sensitivity to noise, migration velocity and surface elevation error. The results prove that our method is less sensitive to noise compared with the conventional elastic least-squares reverse time migration and needs good migration velocities as other least-squares reverse time migration methods. In addition, when implementing the proposed method, an accurate surface elevation should be obtained by global positioning system to yield high-quality images.     

三维偏移距平面波有限差分叠前时间偏移

本文提出了中点-半偏移距域内的三维偏移距平面波(offset plane-wave)方程,并给出了其有限差分解法.偏移距平面波可通过对CMP道集进行平面波分解(倾斜叠加或线性Radon变换)生成,然而这样做会产生严重的噪音干扰.本文提出了局部倾斜叠加方法(local slant-stacking)来消除离散线性Radon变换引入的噪音.针对实际三维数据的不规则性(中点-偏移距域内方位角展布不均匀及偏移距采样不规则),本文还提出了与方位角无关的三维倾斜叠加方法(azimuth-independent 3D slant-stacking),解决了三维平面波分解中存在的问题.使用文中提出的平面波分解方法,可以得到高信噪比的偏移距平面波数据体.同时,三维偏移距平面波偏移可以输出偏移距射线参数域共成像点道集,基于此道集的剩余速度分析方法可以用来更新偏移速度场.偏移距平面波偏移具有很高的计算效率,相较Kirchhoff积分叠前时间偏移有较好的保幅特性,可作为水平地表三维叠前时间偏移的一个很好的解决方案.     

Directional‐oriented wavefield imaging: a new wave‐based subsurface illumination imaging condition for reverse time migration

The key objective of an imaging algorithm is to produce accurate and high‐resolution images of the subsurface geology. However, significant wavefield distortions occur due to wave propagation through complex structures and irregular acquisition geometries causing uneven wavefield illumination at the target. Therefore, conventional imaging conditions are unable to correctly compensate for variable illumination effects. We propose a generalised wave‐based imaging condition, which incorporates a weighting function based on energy illumination at each subsurface reflection and azimuth angles. Our proposed imaging kernel, named as the directional‐oriented wavefield imaging, compensates for illumination effects produced by possible surface obstructions during acquisition, sparse geometries employed in the field, and complex velocity models. An integral part of the directional‐oriented wavefield imaging condition is a methodology for applying down‐going/up‐going wavefield decomposition to both source and receiver extrapolated wavefields. This type of wavefield decomposition eliminates low‐frequency artefacts and scattering noise caused by the two‐way wave equation and can facilitate the robust estimation for energy fluxes of wavefields required for the seismic illumination analysis. Then, based on the estimation of the respective wavefield propagation vectors and associated directions, we evaluate the illumination energy for each subsurface location as a function of image depth point and subsurface azimuth and reflection angles. Thus, the final directional‐oriented wavefield imaging kernel is a cross‐correlation of the decomposed source and receiver wavefields weighted by the illuminated energy estimated at each depth location. The application of the directional‐oriented wavefield imaging condition can be employed during the generation of both depth‐stacked images and azimuth–reflection angle‐domain common image gathers. Numerical examples using synthetic and real data demonstrate that the new imaging condition can properly image complex wave paths and produce high‐fidelity depth sections.