共查询到19条相似文献，搜索用时 125 毫秒
1.
逆时偏移是识别地下高陡构造的有效方法,通常需利用高精度双程波动方程进行波场模拟是该算法的核心,有限差分是求解波动方程常用方式,但会受到模型速度参数、空间网格间距和时间采样间隔的限制,若采用互相关作为逆时偏移的成像条件,虽能够对复杂构造精确成像,可传统的互相关成像条件会引入低频噪声,影响成像结果信噪比.为提升逆时偏移的成像精度,压制低频噪声干扰,打破稳定性条件对地震波场求解的制约,提出了一种解析波场分离方法,基于一步法进行解析波场模拟,利用Hilbert变换实现解析波场的波场分离,进而形成了基于波场分离的逆时偏移成像方法.结合理论模型和VSP实际地震数据,给出了不同波场成分的逆时偏移成像结果.测试结果表明:提出的方法波场模拟稳定性强,且能够高效的进行上下行波波场分离,所得逆时偏移结果不含低频噪声和虚假反射,具有较高的信噪比. 相似文献
2.
逆时偏移作为一种先进的地震偏移成像方法,其成像结果的好坏取决于很多因素,其中成像奈件是关键的一个因素.成像条件一般采用互相关成像条件,它是将源波场沿时间进行正推,接收波场沿时间进行逆推,然后将源波场和接收波场进行互相关,从而得到偏移成像结果.源波场正推和接收波场逆推都是采用同一个波动方程,这里采用波动方程一阶应力-速度形式,数值计算方法采用交错网格有限差分方法.采用振幅补偿拉普拉斯滤波方法压制逆时偏移成像中的低频噪声.采用三次样条插值方法解决逆时偏移成像中的波形不光滑问题.在接收波场逆推过程中,是将地震记录作为边界条件.而在源波场正推过程中,需要给定震源子波作为初值条件.理论上来讲,震源子波为脉冲时是最佳的,但在实际当中难以实现.在源波场正推时,一般采用子波函数.子波函数形式非常多,最常用的是雷克子波.不同主频震源子波产生的源波场频率也不同,从而导致逆时偏移互相关成像效果也不同.那么震源子波主频对逆时偏移成像影响到底有多大呢?为了简便,本文基于不同主频雷克子波,初步分析了相应的逆时偏移成像效果,同时对速度模型存在不同误差时的情形也进行了试验分析.结果 表明:在不同速度误差情况下,震源子波主频越低,成像同相轴连续性越好,散射点聚焦能量越强;震源子波主频越高,虽然分辨率提高,但成像同相轴连续性变差,散射点聚焦能量变弱,同时背景噪声增强.本文研究结果对逆时偏移成像技术在实践中的应用具有一定的参考价值. 相似文献
3.
基于Poynting矢量的归一化波场分离互相关逆时偏移成像条件(英文) 总被引：3，自引：2，他引：1
叠前逆时偏移是目前精度较高的成像方法,然而严重的低频噪声降低了逆时偏移的构造成像精度,偏移噪声的压制是逆时偏移必需要考虑的问题。分析了低频噪声的产生机理,并根据声波方程Poynting矢量的方向指示地震波场的传播方向的原理,分离出上、下、左、右行波,该方法计算量和存储量都远远小于常用的二维傅里叶变换分离方法,进而提出归一化的波场分离互相关成像条件,以压制逆时偏移低频噪声,提高成像精度。实现了Marmousi模型的试算,表明在波场延拓过程中利用Poynting矢量能够较好的分离上、下、左、右行波,与常规方法、拉普拉斯滤波、二维傅里叶变换波场分离的成像结果对比表明,成像时使用归一化的波场分离互相关成像条件能更好的压制偏移噪声,得到精度更高的逆时偏移成像结果。 相似文献
4.
逆时偏移作为一种先进的地震偏移成像方法,其成像结果的好坏取决于很多因素,其中最关键的是速度.速度越准确,成像效果就越好.但在实践当中,速度往往是未知的,只能通过速度反演等各种手段来估计近似速度.由此导致逆时偏移成像结果存在误差,从而降低后续地震解释的精度.那么速度误差对逆时偏移成像的影响到底有多大呢?实际当中一般估计出来的近似速度可能在模型每一个地方都不同,很难直接进行逆时偏移速度误差系统分析.为了简化,本文采用等效速度误差——速度模型整体平均误差,进行逆时偏移成像分析.逆时偏移成像采用互相关成像条件,低频噪声压制采用振幅补偿拉普拉斯滤波方法,源波场正推和接收波场逆推采用波动方程一阶应力.速度形式及交错网格有限差分方法.首先对比分析不同速度误差对逆时偏移成像结果的影响,并和单程波波动方程偏移进行对比分析;然后分析速度误差对逆时偏移和单程波波动方程偏移成像位置的影响;最后基于偏移速度分析方法估计得到的近似速度模型,进行逆时偏移和单程波波动方程偏移成像试验.结果 表明:在不同速度误差情况下,逆时偏移成像在同相轴连续性和能量聚焦等方面要好于单程波波动方程偏移;速度误差对逆时偏移和单程波波动方程偏移成像位置的影响程度相当.本文研究成果对逆时偏移在实践中的应用具有一定的参考价值. 相似文献
5.
逆时偏移是目前应用非常广泛的成像方法,通常利用互相关成像条件进行成像,然而低频噪声的存在降低了成像的精度.为了改善逆时偏移的成像效果,可以通过对波场进行分离,再选择合适的波场成分进行互相关成像.由于Poynting矢量可以指示地下介质中地震波的传播方向,利用Poynting矢量对波场进行了上、下、左、右行波分离.然后,本文将基于Poynting矢量的全波场分离成像条件应用于VTI介质叠前逆时偏移,该成像条件在保证构造准确成像的前提下有效地压制了低频噪声.在进行正演模拟的过程中,本文通过最小二乘差分系数优化提高了模拟的精度.通过模型试算,验证了最小二乘差分系数优化对提高数值模拟精度的效果,并证明了基于Poynting矢量的全波场分离成像条件的有效性. 相似文献
6.
逆时偏移成像是基于双程波波动方程,通过对波场正向和反向延拓来实现偏移的高精度成像方法.逆时偏移能有效利用全波场信息,且相对于单程波偏移方法而言,无传播方向限制,不需要对波场进行上、下行波的分离;将波场看作矢量波的叠加,避免了倾角限制,并能较好地实现回转波、多次波成像.在给定速度场足够准确的情况下,逆时偏移理论上能对任意复杂介质精准成像,是现在最为精确的偏移方法之一.但是,低频噪声是制约其发展的因素之一.本文在分析各种成像条件优缺点的基础上,选用归一化互相关成像条件,对多个典型模型进行逆时偏移成像测试,以验证程序的正确性;选取波场分离、Poynting矢量、拉普拉斯算子滤波等三种去噪方法,分别对组合模型和Marmousi2模型进行测试,对比三种去噪方法在叠前逆时偏移成像中的去噪效果,为实际地震资料处理提供方法借鉴. 相似文献
7.
双程波逆时偏移方法具有成像精度高、无倾角限制和能够适应任意复杂速度介质等诸多优点,受到人们的广泛关注。常规逆时偏移采用互相关成像条件,其成像结果中包含了严重的低频噪音和成像噪音。现有的波场分离成像条件可以压制这些噪声,但成像结果中依然存在着不同程度的噪音。本文在常规波场分离成像条件的基础上,提出了一种二维的多方向波场分离成像条件,将震源和检波点波场分别分离为八个方向的子波场,并进行选择性的互相关成像,叠加得到最终的成像结果。数值算例结果表明,本文方法能够有效地消除低频噪音和成像噪音,得到高精度的成像剖面。 相似文献
8.
采用曲线网格有限差分法描述复杂起伏地形(或不规则波阻抗界面)时,波场正演中可以避免因阶梯近似导致的虚假散射,进而波场逆时偏移可对起伏地表模型进行准确成像.文中以弹性波逆时偏移理论为基础,求解一阶速度-应力方程,推导出了弹性波正向传播和逆时传播的曲线网格差分格式,使用完全匹配吸收边界压制边界反射,采用互相关成像条件,实现了起伏层状介质中的波场逆时偏移.三层起伏、尖灭模型,以及起伏地表条件下的部分盐丘模型结果表明:曲线网格有限差分法逆时偏移法是一种高效、准确的逆时偏移法. 相似文献
9.
10.
尽管叠前逆时偏移成像精度高,但仅针对单一纵波的成像也可能形成地下介质成像盲区,由于基于弹性波方程的逆时偏移成像可形成多波模式的成像数据,因此弹性波逆时偏移成像可提供更为丰富的地下构造信息.本文依据各向同性介质的一阶速度-应力方程组构建震源和检波点矢量波场,再利用Helmholtz分解提取纯纵波和纯横波波场,使用震源归一化的互相关成像条件获得纯波成像,避免了直接使用坐标分量成像而引起的纵横波串扰问题.针对转换波成像的极性反转问题,文中提出一种共炮域极性校正方法.为有效节约存储成本,也提出一种适用于弹性波逆时偏移的震源波场逆时重建方法,在震源波场正传过程中,仅保存PML边界内若干层的速度分量波场,进而逆时重建出所有分量的震源波场.本文分别对地堑模型和Marmousi2模型进行了弹性波逆时偏移成像测试,结果表明:所提出的共炮域极性校正方法正确有效,基于波场分离的弹性波逆时偏移成像的纯波数据能够对复杂地下构造准确成像. 相似文献
11.
Directional‐oriented wavefield imaging: a new wave‐based subsurface illumination imaging condition for reverse time migration 
下载免费PDF全文
下载免费PDF全文 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. 相似文献
12.
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. 相似文献
13.
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. 相似文献
14.
纵横波波场分离是弹性波偏移方法的必要条件,通过纵横波成像的差异可以获取更多地下介质的信息.目前所用的纵横波波场分离方法多采用Helmholtz分解,这样得到的波场不仅物理意义发生了变化,振幅和相位也会发生改变.本文采用纵横波解耦的弹性波方程,将其应用于三维介质,对比分析了纵横波解耦方法相对传统Helmholtz分解方法在相位、振幅上的优势.将该解耦的波场分离方法应用于弹性波逆时偏移,能得到相位、振幅和物理意义不受改变的偏移结果.但是该解耦方法分离得到的纵横波波场均为矢量场,将该波场分离方法用于弹性波逆时偏移,还需要解决矢量场如何得到标量成像结果的问题.本文引入了Poynting矢量,通过Poynting矢量对矢量波场进行标量化,这样就能得到保振幅、相位,且无极性反转的标量PP和PS成像结果.同时针对S波Poynting矢量求取不准确的问题,采用拟S波应力场和S波速度场得到了更加准确的S波Poynting矢量.理论计算证明了本文采用的3D波场解耦的矢量波场分离方法的正确性和引入Poynting矢量对矢量波场进行标量成像的有效性. 相似文献
15.
叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践. 相似文献
16.
与地面地震资料相比,VSP资料具有分辨率高、环境噪声小及能更好地反映井旁信息等优点.常规VSP偏移主要对上行反射波进行成像,存在照明度低、成像范围受限等问题.为了增加照明度、拓宽成像范围、提高成像精度,本文采用直达波除外的所有声波波场数据(全波),包括一次反射波、多次反射波等进行叠前逆时偏移成像.针对逆时偏移中的四个关键问题,即波场延拓、吸收边界条件、成像条件及低频噪声的压制,本文分别采用自适应变空间差分算子长度的优化有限差分方法(自适应优化有限差分方法)求解二维声波波动方程以实现高精度、高效率的波场延拓,采用混合吸收边界条件压制因计算区域有限所引起的人工边界反射,采用震源归一化零延迟互相关成像条件进行成像,采用拉普拉斯滤波方法压制逆时偏移中产生的低频噪声.本文对VSP模型数据的逆时偏移成像进行了分析,结果表明:自适应优化有限差分方法比传统有限差分方法具有更高的模拟精度与计算效率,适用于VSP逆时偏移成像;全波场VSP逆时偏移成像比上行波VSP逆时偏移的成像范围大、成像效果好;相对于反褶积成像条件,震源归一化零延迟互相关成像条件具有稳定性好、计算效率高等优点.将本文方法应用于某实际VSP资料的逆时偏移成像,进一步验证了本文方法的正确性和有效性. 相似文献
17.
Yanadet Sripanich Sergey Fomel Junzhe Sun Jiubing Cheng 《Geophysical Prospecting》2017,65(5):1231-1245
The goal of wave‐mode separation and wave‐vector decomposition is to separate a full elastic wavefield into three wavefields with each corresponding to a different wave mode. This allows elastic reverse‐time migration to handle each wave mode independently. Several of the previously proposed methods to accomplish this task require the knowledge of the polarisation vectors of all three wave modes in a given anisotropic medium. We propose a wave‐vector decomposition method where the wavefield is decomposed in the wavenumber domain via the analytical decomposition operator with improved computational efficiency using low‐rank approximations. The method is applicable for general heterogeneous anisotropic media. To apply the proposed method in low‐symmetry anisotropic media such as orthorhombic, monoclinic, and triclinic, we define the two S modes by sorting them based on their phase velocities (S1 and S2), which are defined everywhere except at the singularities. The singularities can be located using an analytical condition derived from the exact phase‐velocity expressions for S waves. This condition defines a weight function, which can be applied to attenuate the planar artefacts caused by the local discontinuity of polarisation vectors at the singularities. The amplitude information lost because of weighting can be recovered using the technique of local signal–noise orthogonalisation. Numerical examples show that the proposed approach provides an effective decomposition method for all wave modes in heterogeneous, strongly anisotropic media. 相似文献
18.
叠前逆时偏移是当前最为准确的地震成像方法,由于计算量大、存储量大等原因需要合适的实现策略和高效的计算平台.本文以高阶有限差分逆时偏移为基础,重点讨论了在GPU上实现需要解决的显存不足问题和人工边界问题.利用区域分解技术可以在当前GPU上高效地实现任意生产规模的三维逆时偏移成像,不会受到GPU显存规模的制约.常规最佳匹配层边界条件边界区域控制方程与内部区域差异较大,不适于GPU高速运算.本文在GPU上实现近似最佳匹配层(NPML)边界条件,使得高阶有限差分计算不需要分支判断,边界区域辅助波场的存储量也较低,保证了在GPU上进行波场传播的高效性.三维理论数据和实际资料成像结果表明了本文方法的正确性. 相似文献
19.
角度域共成像点道集是衔接叠前地震数据与储层特征的重要桥梁,对地震偏移成像与储层描述具有重要意义.与克希霍夫偏移和单向波动方程偏移相比,逆时偏移是复杂地区最精确的成像方法.高效稳健地生成逆时偏移角度道集目前仍然是一个挑战.本文主要讨论如何采用光学流方法高效、高质量地提取角度道集.在逆时偏移波场外推过程中,光学流方法可以估计波场传播方向.其中Lucas-Kanade(LK)和Horn-Schunck(HS)方法是光学流方法中两种典型的方法.LK光学流方法是一种局部方法,该方法依赖于局部点的梯度值,但是容易出现奇异现象,HS光学流方法属于全局方法,波场方向估计依赖于整个波场,易受噪声影响,对异常值比较敏感,导致整体波场方向计算精度不高.本文提出采用局部和整体结合(Combining Local and Global,CLG)的光学流方法估计波场传播方向.该方法可以有效地提高波场方向的精度,并且简单高效,便于并行处理.对比HS光学流方法,CLG光学流方法几乎不增加额外的计算量.另外,为了弱化光学流方法无法处理波前重叠问题,本文利用解析波场和方向滤波对波场进行方向分解,仅需波场的空间傅里叶变换即可实现任意波场方向分解,将分解后的波场分别估计波场反向,提取成像结果.进一步地,在估计反射张角和方位角时,本文提出有效的归一化方法和改进的最小二乘除法,提高角度估计的精度和稳定性.最后,理论和实际资料例证了本文提出方法的有效性. 相似文献