共查询到19条相似文献，搜索用时 484 毫秒
1.
本文针对弹性波逆时偏移,提出稳定的激发振幅成像条件.在震源波场的正向传播过程中,计算每个网格点的能量,并保存最大能量密度的时刻和相应的波场值;在检波器波场的逆时传播过程中,在每个网格点提取最大能量密度时刻的检波器波场值,并利用保存的最大能量震源波场做归一化,获得角度依赖的反射系数成像剖面.相比于归一化互相关成像条件,该成像条件在震源波场的正向传播过程中无需存储波场快照,节省大量磁盘空间和I/O吞吐任务,提高了计算效率;相比于弹性波的激发时间成像条件,该成像条件自动校正了水平分量在震源两侧的极性反转,在多炮叠加时避免振幅损失.数值试验表明,与归一化成像条件相比,稳定激发振幅成像条件具有更小的计算量,偏移剖面的低频假象更弱,水平分量的成像能力更优,具有更高的空间分辨率. 相似文献
2.
基于Poynting矢量的归一化波场分离互相关逆时偏移成像条件(英文) 总被引：3，自引：2，他引：1
叠前逆时偏移是目前精度较高的成像方法,然而严重的低频噪声降低了逆时偏移的构造成像精度,偏移噪声的压制是逆时偏移必需要考虑的问题。分析了低频噪声的产生机理,并根据声波方程Poynting矢量的方向指示地震波场的传播方向的原理,分离出上、下、左、右行波,该方法计算量和存储量都远远小于常用的二维傅里叶变换分离方法,进而提出归一化的波场分离互相关成像条件,以压制逆时偏移低频噪声,提高成像精度。实现了Marmousi模型的试算,表明在波场延拓过程中利用Poynting矢量能够较好的分离上、下、左、右行波,与常规方法、拉普拉斯滤波、二维傅里叶变换波场分离的成像结果对比表明,成像时使用归一化的波场分离互相关成像条件能更好的压制偏移噪声,得到精度更高的逆时偏移成像结果。 相似文献
3.
4.
逆时偏移作为一种先进的地震偏移成像方法,其成像结果的好坏取决于很多因素,其中成像奈件是关键的一个因素.成像条件一般采用互相关成像条件,它是将源波场沿时间进行正推,接收波场沿时间进行逆推,然后将源波场和接收波场进行互相关,从而得到偏移成像结果.源波场正推和接收波场逆推都是采用同一个波动方程,这里采用波动方程一阶应力-速度形式,数值计算方法采用交错网格有限差分方法.采用振幅补偿拉普拉斯滤波方法压制逆时偏移成像中的低频噪声.采用三次样条插值方法解决逆时偏移成像中的波形不光滑问题.在接收波场逆推过程中,是将地震记录作为边界条件.而在源波场正推过程中,需要给定震源子波作为初值条件.理论上来讲,震源子波为脉冲时是最佳的,但在实际当中难以实现.在源波场正推时,一般采用子波函数.子波函数形式非常多,最常用的是雷克子波.不同主频震源子波产生的源波场频率也不同,从而导致逆时偏移互相关成像效果也不同.那么震源子波主频对逆时偏移成像影响到底有多大呢?为了简便,本文基于不同主频雷克子波,初步分析了相应的逆时偏移成像效果,同时对速度模型存在不同误差时的情形也进行了试验分析.结果 表明:在不同速度误差情况下,震源子波主频越低,成像同相轴连续性越好,散射点聚焦能量越强;震源子波主频越高,虽然分辨率提高,但成像同相轴连续性变差,散射点聚焦能量变弱,同时背景噪声增强.本文研究结果对逆时偏移成像技术在实践中的应用具有一定的参考价值. 相似文献
5.
逆时偏移是识别地下高陡构造的有效方法,通常需利用高精度双程波动方程进行波场模拟是该算法的核心,有限差分是求解波动方程常用方式,但会受到模型速度参数、空间网格间距和时间采样间隔的限制,若采用互相关作为逆时偏移的成像条件,虽能够对复杂构造精确成像,可传统的互相关成像条件会引入低频噪声,影响成像结果信噪比.为提升逆时偏移的成像精度,压制低频噪声干扰,打破稳定性条件对地震波场求解的制约,提出了一种解析波场分离方法,基于一步法进行解析波场模拟,利用Hilbert变换实现解析波场的波场分离,进而形成了基于波场分离的逆时偏移成像方法.结合理论模型和VSP实际地震数据,给出了不同波场成分的逆时偏移成像结果.测试结果表明:提出的方法波场模拟稳定性强,且能够高效的进行上下行波波场分离,所得逆时偏移结果不含低频噪声和虚假反射,具有较高的信噪比. 相似文献
6.
7.
逆时偏移是目前应用非常广泛的成像方法,通常利用互相关成像条件进行成像,然而低频噪声的存在降低了成像的精度.为了改善逆时偏移的成像效果,可以通过对波场进行分离,再选择合适的波场成分进行互相关成像.由于Poynting矢量可以指示地下介质中地震波的传播方向,利用Poynting矢量对波场进行了上、下、左、右行波分离.然后,本文将基于Poynting矢量的全波场分离成像条件应用于VTI介质叠前逆时偏移,该成像条件在保证构造准确成像的前提下有效地压制了低频噪声.在进行正演模拟的过程中,本文通过最小二乘差分系数优化提高了模拟的精度.通过模型试算,验证了最小二乘差分系数优化对提高数值模拟精度的效果,并证明了基于Poynting矢量的全波场分离成像条件的有效性. 相似文献
8.
双程波逆时偏移方法具有成像精度高、无倾角限制和能够适应任意复杂速度介质等诸多优点,受到人们的广泛关注。常规逆时偏移采用互相关成像条件,其成像结果中包含了严重的低频噪音和成像噪音。现有的波场分离成像条件可以压制这些噪声,但成像结果中依然存在着不同程度的噪音。本文在常规波场分离成像条件的基础上,提出了一种二维的多方向波场分离成像条件,将震源和检波点波场分别分离为八个方向的子波场,并进行选择性的互相关成像,叠加得到最终的成像结果。数值算例结果表明,本文方法能够有效地消除低频噪音和成像噪音,得到高精度的成像剖面。 相似文献
9.
与地面地震资料相比,VSP资料具有分辨率高、环境噪声小及能更好地反映井旁信息等优点.常规VSP偏移主要对上行反射波进行成像,存在照明度低、成像范围受限等问题.为了增加照明度、拓宽成像范围、提高成像精度,本文采用直达波除外的所有声波波场数据(全波),包括一次反射波、多次反射波等进行叠前逆时偏移成像.针对逆时偏移中的四个关键问题,即波场延拓、吸收边界条件、成像条件及低频噪声的压制,本文分别采用自适应变空间差分算子长度的优化有限差分方法(自适应优化有限差分方法)求解二维声波波动方程以实现高精度、高效率的波场延拓,采用混合吸收边界条件压制因计算区域有限所引起的人工边界反射,采用震源归一化零延迟互相关成像条件进行成像,采用拉普拉斯滤波方法压制逆时偏移中产生的低频噪声.本文对VSP模型数据的逆时偏移成像进行了分析,结果表明:自适应优化有限差分方法比传统有限差分方法具有更高的模拟精度与计算效率,适用于VSP逆时偏移成像;全波场VSP逆时偏移成像比上行波VSP逆时偏移的成像范围大、成像效果好;相对于反褶积成像条件,震源归一化零延迟互相关成像条件具有稳定性好、计算效率高等优点.将本文方法应用于某实际VSP资料的逆时偏移成像,进一步验证了本文方法的正确性和有效性. 相似文献
10.
逆时偏移方法利用双程波算子模拟波场的正向和反向传播,通常采用互相关成像条件获得偏移剖面,是一种高精度的成像方法.但是传统的互相关成像条件会在偏移结果中产生低频噪声;此外,如果偏移速度中存在剧烈速度变化还可能进一步产生偏移假象.为了提高逆时偏移的成像质量,可在成像过程中先对震源波场和检波点波场分别进行波场分离,然后选择合适的波场成分进行互相关成像.本文基于Hilbert变换,推导了可在偏移过程中进行上下行和左右行波场分离的高效波场分离公式以及相应的成像条件,结合Sigsbee 2B合成数据,给出了不同波场成分的互相关成像结果.数值算例结果表明,采用本文提出的高效波场分离算法以及合理的波场成分互相关成像条件可以获得高信噪比的成像结果. 相似文献
11.
双程波方程逆时深度偏移是复杂介质高精度成像的有效技术,但其结果中通常包含成像方法引起的噪音和假象,一般的滤波方法会破坏成像剖面上的振幅,其中的假象也会给后续地质解释带来困扰.将波场进行方向分解然后实现入射波与反射波的相关成像能够有效地消除这类成像噪音,并提高逆时偏移成像质量.波传播方向的分解通常在频率波数域实现,它会占用大量的存储和计算资源,不便于在沿时间外推的逆时深度偏移中应用.本文提出解析时间波场外推方法,可以在时间外推的每个时间片上实现波传播方向的显式分解,逆时深度偏移中利用分解后的炮检波场进行对应的相关运算,实现成像噪音和成像信号的分离.在模型和实际数据上的测试表明,相比于常规互相关逆时偏移成像结果,本文方法能够有效地消除低频成像噪音和特殊地质构造导致的成像假象. 相似文献
12.
Reverse time migration is an advanced seismic migration imaging method. When the source wavefield and the receiver wavefield are cross-correlated, the cross-correlations of direct arrivals, backscattered waves and overturned waves will produce a lot of low-frequency noise, which will mask the final imaging results. Laplacian filtering, as a common method to suppress low-frequency noise, can adapt to any complex media, just adding a little computational cost. However, simple direct Laplacian filtering will destroy the characteristics of the useful signals. Therefore, the amplitude needs to be compensated before filtering when using the Laplacian filtering method. Zhang and Sun proposed an improved Laplacian filtering method and gave a simple calculation formula and explanation. This method can effectively suppress the low-frequency noise in reverse time migration while retaining the useful signal characteristics, but lacks detailed and strict mathematical derivation. Therefore, this paper gives a detailed and rigorous mathematical derivation of the amplitude-compensated Laplace filtering method from the point of view of amplitude-preserved filtering. The source wavelet is used instead of the source wavefield to compensate amplitude, just adding a little calculation cost. Finally, the amplitude-compensated Laplace filtering method is verified by two theoretical models and compared with the direct Laplacian filtering method. 相似文献
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.
Shaoyong Liu Zhe Yan Hanming Gu Yongjie Tang Chuncheng Liu 《Geophysical Prospecting》2019,67(3):496-507
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. 相似文献
15.
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. 相似文献
16.
相对于传统的逆时偏移,最小二乘逆时偏移具有更高的成像质量,这种改善是通过迭代反演来获得的,另外其精度与效率依赖于求解波动方程算法的精度与效率.本文给出了基于非结构化网格的最小二乘逆时偏移的方法,该方法能够充分结合最小二乘算法与非结构化网格精细刻画地下界面以及随速度自适应剖分的优点;并采用带补偿的拉普拉斯滤波算法,来消除梯度计算中的低频噪声,从而加速目标函数的收敛速度.通过简单倾斜模型以及复杂Marmousi模型测试,显示了该方法的有效性和潜力. 相似文献
17.
常规的单程波波动方程偏移成像方法对大角度的高陡构造偏移成像存在内在的限制.根据波动方程在各个空间方向的数学特性和高陡构造反射地震波的传播特征,通过把地震波分解为垂向的上下行波、水平方向的前后行波和左右行波,提出基于波场垂向外推和水平方向外推相结合的单程波波动方程高陡构造偏移成像方法,即用波场垂向外推的单程波波动方程偏移成像方法解决中低角度平缓构造的偏移成像,用波场水平方向外推的单程波波动方程偏移成像方法解决中高角度陡倾构造的偏移成像.这种基于波场垂向和水平方向外推相结合的高陡构造偏移成像方法是常规单程波波动方程叠前深度偏移成像方法的补充和改进,它相对基于全波方程的逆时偏移具有计算效率上的优势. 相似文献
18.
本文研究分析了双程波波动方程偏移成像中广泛存在的三种主要噪声,特别是针对过去研究中没能很好解决的存在于高速盐丘悬垂边界附近的射线状噪声,提出了基于优化成像条件的有效去噪方法。射线状噪声主要来自于震源一侧波场的下行透射波分量和接收阵列一侧波场的上行散射波互相关成像。这一部分能量具有较强的互相关性,但并不携带真实的反射面信息。它广泛存在叠前偏移成像中,与信号在强度上同量级。多数情况下偏移成像中的相关噪声由方向性传播的波场能量产生。利用波场梯度得到的波场传播角度,可以分离出噪声对应的波场能量,并在成像条件中减去。采用这一方法可以有效地去除多种噪声,包括直达波噪声、散射波噪声和射线状噪声。该去噪方法不依赖波场外推算子,在需要时可以方便地运用到几乎所有的波动方程偏移中去。并且该去噪方法针对噪声的物理根源,对信号的损害很小。对去噪后的偏移成像结果额外地进行波数域滤波处理,可以进一步提高叠加图像的质量。这一去噪方法在超广角单程波偏移成像中取得良好效果,我们同时期待其在其他双程波波动方程偏移特别是逆时偏移(RTM)中的成功运用。 相似文献
19.
针对基于互相关成像条件的探地雷达（GPR）逆时偏移计算效率低、存储量大及易产生低频假象的不足，本文将激发振幅成像条件应用于GPR逆时偏移成像中.通过在源点电磁波场正向传播过程计算每个网格点的能量密度，并保存最大能量密度的时刻和相应的电磁波场值；在接收点电磁波场逆向传播过程提取每个网格点最大能量密度时刻及对应的电磁波场值，并利用保存的最大能量源点电磁波场及走时做归一化，从而获得了依赖反射系数成像剖面，避免了源点正向传播电磁波场的存储和重建.此外，为了提高电磁波场的模拟精度，采用了基于三角形剖分的时间域有限单元法（FETD）计算电磁波正向和逆向传播过程.最后通过模型试算表明：激发振幅成像条件相比于归一化互相关成像条件，成像结果低频噪声更弱，空间分辨率更高，计算效率提高了近2倍. 相似文献