首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 140 毫秒
1.
弹性波逆时偏移是一种以矢量波理论为基础的深度域偏移算法,该算法在波场传播过程中能够正确处理波形能量的转换,保持地震波的动力学和运动学特征.然而,传统弹性波逆时偏移由于成像时纵横波未分离,导致串扰噪声的产生,降低成像精度.为压制弹性波逆时偏移成像串扰噪声,本文对基于Helmholtz分解的波场分离与基于解耦方程的波场分离方法进行研究和分析,并测试对比相应的弹性波逆时偏移成像效果.均匀各向同性介质波场分离测试表明,基于Helmholtz分解分离波场的振幅与相位均发生改变,需要额外校正;而基于解耦方程的波场分离方法可以保持原始波场相位与振幅特征不变.Marmousi截断模型测试表明,基于Helmholtz分解分离波场的弹性波逆时偏移成像振幅与相位发生改变,且存在极性反转现象;基于解耦方程波场分离方法的弹性波逆时偏移成像保持原始波场矢量特征,无极性反转现象,成像效果较为理想.  相似文献   

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

3.
欲实现基于弹性波方程的矢量波场逆时偏移纵、横波独立成像,必须在波场延拓过程中实现纵、横波场的分离,散度和旋度算子分离的纵、横波出现振幅与相位的畸变,导致输出成像结果的振幅失真。本文提出一种在弹性波场延拓过程中实现纵、横波保幅分离的方法,在传统的弹性波方程中加入纵波压力、纵波振动速度和横波振动速度方程,实现纵横波的矢量分解,再对分解后的矢量纵波和矢量横波做标量化合成得到保幅分离的纵、横波场,对保幅分离的纵、横波场应用成像条件,然后实现矢量波场逆时偏移的保幅纵横波成像。该方法可以保证分离后纵、横波的振幅与相位不变;同时,分解后的纵波压力和纵波振动速度可用于层间反射噪音压制和横波极性校正,提高多分量地震资料联合逆时偏移的纵、横波成像质量,从而实现保幅弹性波逆时偏移的目的,为叠前深度剖面应用于叠前反演工作奠定基础。  相似文献   

4.
弹性波逆时偏移是当前多分量地震资料相对准确的偏移算法,它能够形成多波模式的成像剖面,从而减少纵波勘探的多解性.本文首先依据各向同性介质中矢量分离纵横波场的速度-应力方程组,利用高阶交错网格有限差分数值方法求解弹性波方程,进而构建矢量的纵波和横波波场,不同于散度和旋度算子分离纵横波场的传统方法,文中提出的矢量分离纵横波场方法保持了原始波场的振幅和相位特征.文中也提出将震源归一化的内积成像条件应用于分离后的纯纵波和横波矢量场,由此得到的转换波成像避免了传统弹性波成像方法中出现的极性反转.水平层状和复杂构造模型测试表明,文中提出的基于矢量分离纵横波场的弹性波逆时偏移方法成像精度高,转换波成像PS和SP极性无反转,所形成的多种模式纯波剖面能够准确地对复杂地下构造成像.  相似文献   

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

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

7.
李闻达  孟小红  刘洪  王建  孙军  桂生 《地球物理学报》2019,62(11):4378-4392
弹性波矢量波场逆时偏移可以综合利用纵横波场信息,对地下空间进行清晰成像,且对于成像介质没有角度限制,可以对复杂构造进行更清晰的成像.而弹性波逆时偏移中最重要的就是求解波动方程的算法,其直接影响成像的精度以及效率.本文引进电力系统谐波分析中常用的余弦组合窗函数,并通过一种新的优化算法得到了改进的余弦组合窗函数从而得到优化后的有限差分算子.并将此算子应用于解耦方程的矢量波场分离算法从而提高了成像精度.数值测试表明基于新算法的逆时偏移的成像精度和清晰度得到了明显的提高.  相似文献   

8.
逆时偏移成像是基于双程波波动方程,通过对波场正向和反向延拓来实现偏移的高精度成像方法.逆时偏移能有效利用全波场信息,且相对于单程波偏移方法而言,无传播方向限制,不需要对波场进行上、下行波的分离;将波场看作矢量波的叠加,避免了倾角限制,并能较好地实现回转波、多次波成像.在给定速度场足够准确的情况下,逆时偏移理论上能对任意复杂介质精准成像,是现在最为精确的偏移方法之一.但是,低频噪声是制约其发展的因素之一.本文在分析各种成像条件优缺点的基础上,选用归一化互相关成像条件,对多个典型模型进行逆时偏移成像测试,以验证程序的正确性;选取波场分离、Poynting矢量、拉普拉斯算子滤波等三种去噪方法,分别对组合模型和Marmousi2模型进行测试,对比三种去噪方法在叠前逆时偏移成像中的去噪效果,为实际地震资料处理提供方法借鉴.  相似文献   

9.
在海洋地震勘探中,海底电缆技术采集的多分量地震数据,涉及到多参数的反演与成像问题,本文针对海底电缆多分量地震数据提出了一种弹性波多参数最小二乘逆时偏移方法。该方法的波场延拓算子为混合方程,即在海水介质中采用声波方程进行波场计算,而海底固体介质的波场由纵横波形分离的矢量弹性波方程得到。在海底界面中采用声弹耦合控制方程将两种类型的方程结合起来。通过推导纵横波形分离的弹性波偏移算子、反偏移算子和梯度公式实现基于纵横波形分离的弹性波最小二乘逆时偏移方法。通过模型试算证明了该方法能够得到高质量的纵波速度和横波速度分量的成像剖面,相比与传统弹性波最小二乘逆时偏移方法,多参数耦合造成的成像串扰噪音得到了很好地压制。  相似文献   

10.
逆时偏移中用Poynting矢量高效地提取角道集   总被引:3,自引:1,他引:2       下载免费PDF全文
逆时偏移在提高复杂介质的成像质量方面表现出了优越的性能,但逆时偏移对速度精度的要求比较高.共成像点道集是一种非常重要的叠前深度偏移输出,它除了能为深度偏移处理提供速度信息外,还能够提供振幅和相位等信息,为后续的属性解释提供依据.本文提出一种在逆时偏移成像过程中提取角度域共成像点道集的方法,该方法采用矢量波动方程进行波场传播,并用能流密度矢量(Poynting vector)计算反射角,最后应用互相关成像条件输出角度域共成像点道集.该方法简单易于实现,且几乎不需要额外的计算量和存储量,非常适合于进行逆时偏移速度分析,同时提出的角道集也能用于进行AVA等分析.最后通过模型算例和实际数据检验了方法的有效性和优越性.  相似文献   

11.
Angle-domain common-image gathers (ADCIGs) transformed from the shotdomain common-offset gathers are input to migration velocity analysis (MVA) and prestack inversion. ADCIGs are non-illusion prestack inversion gathers, and thus, accurate. We studied the extraction of elastic-wave ADCIGs based on amplitude-preserving elastic-wave reversetime migration for calculating the incidence angle of P-and S-waves at each image point and for different source locations. The P-and S-waves share the same incident angle, namely the incident angle of the source P-waves. The angle of incidence of the source P-wavefield was the difference between the source P-wave propagation angle and the reflector dips. The propagation angle of the source P-waves was obtained from the polarization vector of the decomposed P-waves. The reflectors’ normal direction angle was obtained using the complex wavenumber of the stacked reverse-time migration (RTM) images. The ADCIGs of P-and S-waves were obtained by rearranging the common-shot migration gathers based on the incident angle. We used a horizontally layered model, the graben medium model, and part of the Marmousi-II elastic model and field data to test the proposed algorithm. The results suggested that the proposed method can efficiently extract the P-and S-wave ADCIGs of the elastic-wave reverse-time migration, the P-and S-wave incident angle, and the angle-gather amplitude fidelity, and improve the MVA and prestack inversion.  相似文献   

12.
Imaging the PP- and PS-wave for the elastic vector wave reverse-time migration requires separating the P- and S-waves during the wave field extrapolation. The amplitude and phase of the P- and S-waves are distorted when divergence and curl operators are used to separate the P- and S-waves. We present a P- and S-wave amplitude-preserving separation algorithm for the elastic wavefield extrapolation. First, we add the P-wave pressure and P-wave vibration velocity equation to the conventional elastic wave equation to decompose the P- and S-wave vectors. Then, we synthesize the scalar P- and S-wave from the vector Pand S-wave to obtain the scalar P- and S-wave. The amplitude-preserved separated P- and S-waves are imaged based on the vector wave reverse-time migration (RTM). This method ensures that the amplitude and phase of the separated P- and S-wave remain unchanged compared with the divergence and curl operators. In addition, after decomposition, the P-wave pressure and vibration velocity can be used to suppress the interlayer reflection noise and to correct the S-wave polarity. This improves the image quality of P- and S-wave in multicomponent seismic data and the true-amplitude elastic reverse time migration used in prestack inversion.  相似文献   

13.
常规陆上VSP(Vertical Seismic Profiling)勘探普遍采用纵波震源激发,三分量检波器接收,主要利用的是纵波和转换横波信息。已有的研究表明,炸药震源在井下激发、可控震源在地面垂向振动,均会产生较强的纯纵波和一定强度的纯横波;泊松比差别较大的分界面有利于形成较强的透射转换横波。本文通过对激发形成的纯横波和下行转换形成的横波进行对比分析,认为纯横波的主频往往低于纯纵波的主频,而下行转换横波的主频通常接近纵波的主频。本文分别对两个陆上纵波源零偏和非零偏VSP资料进行分析,结果表明这些资料中普遍存在纯横波,只是横波的强弱存在不同程度的变化。利用纵波源零偏VSP资料,可以获得横波速度信。最后对VSP纵波和横波联合应用前景进行了分析,应该充分利用纵波源VSP资料中的横波信息。  相似文献   

14.
以多分量地震观测为基础,联合纵波和转换横波数据能更有效地估计地下介质的弹性和物性参数,提升地质构造成像与油气储层描述的精度.在海底多分量地震数据处理过程中,观测记录的上-下行波分解和P/S波分离可压制水层鸣震以及P与S波之间的串扰,对偏移成像和纵横波速度建模至关重要.但受海底环境、仪器与观测因素共同影响,许多海底多分量地震资料都无法基于现有的海底波场分离方法与流程取得合理的结果.本文以海底声波场与弹性波场分离基本原理为基础,通过对方法流程的修正,摆脱常规流程对中小偏移距直达波信号的依赖性.借助模拟数据实验讨论了波场分离对海底介质参数、噪声的敏感性.结合东海YQ探区海底多分量地震资料上-下行P/S波分离及其叠前深度偏移处理,验证了本文方法流程的可行性.  相似文献   

15.
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.  相似文献   

16.
随着我国勘探开发难度逐步增大,勘探目标开始向裂缝油气藏、岩性油气藏等复杂探区转移,研究高精度、适应性强的多波多分量深度偏移算法在后续的地震解释、属性分析及储层预测中具有重要意义.针对多波多分量地震数据,本文提出了一种二维弹性波时空域高斯束偏移方法.时空域高斯束沿中心射线传播时能够面向成像目标描述局部波场,且对振幅和频率可调制的Gabor基函数有天然的适应性,因而将基于Gabor分解的子波重构方法应用于震源波场构建,从而得到任意点源函数产生的时空域高斯束波场.该方法由于直接在时间域进行计算,可以避开频率域中出现的假频和边缘截断效应等问题.基于各向同性弹性波动方程的Kirchhoff-Helmholtz积分解,利用矢量时空域高斯束传播算子构建格林函数和格林位移张量,并结合上行射线追踪策略,实现了检波点波场的反向延拓.针对矢量波成像问题,本文借鉴弹性波逆时偏移方法从矢量延拓波场中分离出纯纵波分量和纯横波分量,进而采用修改后的内积成像条件产生具有明确物理意义的PP、PS成像结果,避免了转换波成像的极性反转问题.最后利用简单两层模型和不含盐体构造的部分Sigsbee2a模型的成像结果,并将其与应用近似纵横波成像条件、标量和矢量势成像条件的偏移剖面进行对比,验证了本文方法的正确性和有效性.  相似文献   

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

18.
基于弹性波解耦延拓方程的波场分离方法不仅可以得到解耦的纵、横波质点振动速度场,还可以获得纵、横波应力场.针对横波纯应力场在利用单一分量成像时不具有明确物理意义的问题,本文将应力偏张量引入到横波应力场中,基于应力偏量第二不变量构建得到横波应力不变量并将其用于应力场的逆时偏移成像中,获得了可完整表征横波应力场的成像结果.模型试算表明,本文构建的横波应力不变量可以有效利用横波应力张量中的波场信息,并得到准确的弹性逆时偏移成像结果.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号