首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 140 毫秒
1.
Kirchhoff积分时间偏移的两种走时计算及并行算法   总被引:11,自引:7,他引:4       下载免费PDF全文
走时计算的可行性在很大程度上决定了Kirchhoff积分法偏移的效果,而基于PC-Cluster的并行计算是实施叠前偏移数据计算量大的基本工具,本文针对这两个问题的研究,阐明两种旨在提高复杂情况叠前时间偏移聚焦能力的弯曲射线对称走时及非对称走时的计算方法,以及旨在提高叠前偏移计算效率的"数据并行"和"成像点并行"策略,并分别对其结果进行了对比分析.  相似文献   

2.
地震射线走时的求取方法是叠前时间偏移研究的核心问题之一,也是影响计算时间域角道集角度精确性的关键因素之一.本文基于Kirchhoff叠前时间偏移,应用第一类切比雪夫多项式,对弯曲射线对称走时加以改进,引进非对称项,优化后得到切比雪夫非对称走时方程,与高精度走时进行比较和误差分析,再将该走时求取方法应用于时间域角道集的求取中,得到地下较真实的入射角.通过模型计算和实际地震资料处理证明,此种非对称走时及其角道集的求取方法具有精度高、计算量少的优点.  相似文献   

3.
常规长排列非双曲动校正公式是在VTI介质中得到的,它不能满足任意空间取向TI(ATI)条件下的扩展.本文以VTI介质中非双曲动校正公式为基础,基于我们推导得出的ATI介质中精确四次时差系数解析解和NMO速度解析解,给出ATI介质中长排列优化的非双曲动校正公式.通过与各向异性射线追踪方法计算所得出的"精确走时"结果对比,研究表明优化后的非双曲动校正公式能精确地描述任意强弱、ATI介质中随测线方位变化的走时曲线,可以用来替代耗时、多偏移距、多方位的射线追踪方法正演拟合ATI介质中长偏移距反射走时,为利用非双曲时距的各向异性参数反演提供理论基础性认识。  相似文献   

4.
层状各向异性介质转换波克希霍夫叠前时间偏移   总被引:7,自引:2,他引:5       下载免费PDF全文
在克希霍夫叠前时间偏移处理中,地震波走时的计算方法是决定大偏移距地震资料成像品质的重要因素.在常规的三维转换波各向异性叠前时间偏移公式中,走时的计算是基于等效单层各向异性介质的非双曲线方法.用这种方法处理的成像道集,在偏移/深度比超过一定阈值后,成像道集中的反射同相轴将出现过偏现象,这种偏移不平的同相轴将影响偏移叠加的最佳响应,使得偏移成像波组呈低频化特征,最终降低三维转换波偏移成像质量.我们采用层状介质的走时计算方法代替常规算法,并且利用了常规方法的转换波各向异性偏移速度模型.基于层状介质的算法能够提高大偏移距转换波走时计算精度,克服中浅地层大偏移距远道成像道集中反射同相轴逐渐上翘的问题.两个地区的三维转换波资料处理结果证实,基于层状各向异性介质的转换波克希霍夫叠前时间偏移方法,明显改善了反射成像剖面的连续性和分辨率,提高成像剖面构造的可解释性.  相似文献   

5.
转换点位置的计算是转换波资料处理中的一个关键问题. 本文提出了分别基于速度随深度线性变化、速度随垂直走时线性变化、慢度随深度线性变化和慢度随垂直走时线性变化四种等效垂向非均匀介质情况下转换点位置的计算方法. 研究了通过速度拟合、走时近似和相似系数谱三种方式选择合适的等效速度方法. 结合理论模型对非均匀介质转换点计算方法、渐进转换点计算方法、Thomsen近似公式和均匀介质解析计算方法的误差进行了分析,结果表明非均匀介质转换点计算方法能更准确地计算转换点位置.  相似文献   

6.
任意介质中的动态规划法地震波三维走时计算   总被引:1,自引:1,他引:1  
任意介质中的地震波三维走时计算是复杂介质情况下Kirchhoff积分法三维叠前深度偏移及走时层析成像的核心.走时算法的效率及精度决定了成像方法的应用范围及效果,对复杂地质构造区域的地震波成像时需要有稳健的走时计算方法.本文把Schneider等提出的用动态规划法计算二维任意复杂介质中走时的方法推广到三维.此方法的核心是构造从源点到当前计算点的平均慢度,基于Fermat原理,用球面波近似导出走时计算所用的公式,并用动态规划法搜索到达当前计算点的初至走时.它适用于任意复杂的介质情况,对速度差异没有限制,计算过程中考虑到各个可能的方向到达当前计算点的初至波.首波及回转波的初至走时也能正确地计算出来.各种理论速度模型上的走时计算及胜利油田某探区的三维叠前深度偏移的成功实践验证了方法的正确性.  相似文献   

7.
横向均匀介质可由若干水平分层足够好地逼近,本文引入零厚度层来描述成层介质,从而建立了回折波、反射波、首波和直达波统一的时距关系以及地震走时随模型参数变化的公式。本文详细讨论了水平分层介质里地震走时曲线的特征,合成地震走时的试验表明了由地震走时,用阻尼最小二乘法可以成功地反演速度模型参数,所反演的模型含有速度梯度层,也可以含有低速层。  相似文献   

8.
大步长波场深度延拓的理论   总被引:19,自引:4,他引:19       下载免费PDF全文
波场延拓是地震偏移成像的基础. 快速进行目标区波场延拓对石油勘探中急需发展的深部地震勘探和无组合海量地震数据的成像有重要意义. 在目标区成像中,目前已有的波场延拓方法,包括基于走时计算的Dix方法和射线追踪方法,以及基于小步长波场递推的方法,在适应复杂介质、计算精度和计算效率的某一方面还不能完全满足实际需要. 本文提出一种基于“算子相位”李代数积分的快速计算延拓算子的方法,称为大步长波场延拓方法. 在该方法中,指向目标区的波场延拓算子象征的复相位被表示成波数的线性组合. 线性组合的系数是层速度函数及其导数的深度积分,计算和存储较为方便. 波场延拓算子通过相移算子加校正的方法,利用快速Fourier变换在空间域和波数域予以实现. 利用动力学等价关系导出了便于计算的表达式. 本文比较了算子主象征函数用一步法展开和用两步法展开的精度,从而说明大步长方法的精度要高于递推方法. 在横向和纵向线性变化介质中,将大步长方法的脉冲响应与递推法做了比较,说明大步长延拓算子的走时精度主要取决于相移因子中的横向变速校正项;且在各种近似下,大步长算子发生的频散都非常小.  相似文献   

9.
基于绝对走时计算的地面微地震监测地震发射层析成像(Seismic emission tomography,SET)定位算法的定位精度受地层速度模型的准确性影响很大,速度模型中的误差会使定位结果出现偏差.本文提出了一种基于相对走时计算的地面微地震监测SET定位算法,以克服绝对走时计算对速度模型精度的高度依赖.首先,通过分析走时误差推导了相对走时计算公式.然后,利用该公式对射孔初至进行校正,得到射孔点附近其他待监测网格点的走时.通过数值模拟和实际资料测试,结果表明在基于SET定位算法的前提下,本文提出的相对走时计算方法能得到比绝对走时计算方法更准确的定位结果.  相似文献   

10.
从单程波方程出发,推导出角度域叠前时间偏移的走时、入射波与反射波夹角、成像幅值计算方法;构建了可直接生成角度域成像道集的叠前时间偏移方法与偏移流程.文中定量分析了速度梯度对走时、角度、幅值的影响,给出了可更好考虑介质非均匀性的角度、幅值计算方法.理论模型以及实际数据验证了本文方法的有效性.  相似文献   

11.
波动方程偏移速度建模:一种直接反演方法   总被引:1,自引:1,他引:0       下载免费PDF全文
基于炮域波动方程叠前深度偏移所产出的角道集,发展了不需要多次迭代的波动方程偏移速度建模方法.文中分析了炮域波动方程偏移生成角道集的方法,给出了非均匀介质中角道集同相轴曲率与速度误差的定量关系,发展了基于同相轴曲率的速度模型直接更新算法.在初始模型较合理的情况下,应用一次炮域波动方程偏移计算,即可得到较准确的速度模型.这...  相似文献   

12.
The study of wave propagation in media with elliptical velocity anisotropy shows that seismic energy is focused according to the horizontal component of the velocity field while the vertical component controls the time-to-depth relation. This implies that the vertical component cannot be determined from surface seismic velocity analysis but must be obtained using borehole or regional geological information. Both components of the velocity field are required to produce a correctly focused depth image. A paraxial wave equation is developed for elliptical anisotropic wave propagation which can be used for modelling or migration. This equation is then transformed by a change of variable to a second paraxial equation which only depends on one effective velocity field. A complete anisotropic depth migration using this transformed equation involves an imaging step followed by a depth stretching operation. This allows an approximate separation or splitting of the focusing and depth conversion steps of depth migration allowing a different velocity model to be used for each step. This split anisotropic depth migration produces a more accurate result than that obtained by a time migration using the horizontal velocity field followed by an image-ray depth conversion using the vertical velocity field. The results are also more accurate than isotropic depth migration and yield accurate imaging in depth as long as the lateral variations in the anisotropy are slow.  相似文献   

13.
Conventional Kirchhoff prestack time migration based on the hyperbolic moveout can cause ambiguity in laterally inhomogeneous media, because the root mean square velocity corresponds to a one-dimensional model under the horizontal layer assumption; it does not include the lateral variations. The shot/receiver configuration with different offsets and azimuths should adopt different migration velocities as they contribute to a single image point. Therefore, we propose to use an offset-vector to describe the lateral variations through an offset-dependent velocity corresponding to the difference in offset from surface points to the image point. The offset-vector is decomposed into orthogonal directions along the in-line and cross-line directions so that the single velocity can be expressed as a series of actual velocities. We use a simple Snell's law-based ray tracing to calculate the travel time recorded at the image point and convert the travel time to an equivalent velocity corresponding to a pseudo-straight ray. The double-square-root equation using such an equivalent velocity in the offset-vector domain is non-hyperbolic and asymmetrical, which improves the accuracy of the migration. Numerical examples using the Marmousi model and a wide azimuth field data show that the proposed method can achieve reasonable accuracy and significantly enhances the imaging of complex structures.  相似文献   

14.
任意复杂介质中主能量法地震波走时计算   总被引:14,自引:1,他引:13       下载免费PDF全文
积分法叠前深度偏移及层析成像的核心是复杂介质情况下的地震波走时计算. 复杂构造的高精度地震成像需要有稳健的走时计算方法。本文把 Nichols提出的用地震波主能量计算走时的方法由二维推广到三维,并推导出三维波动方程Helmholtz形式在球坐标系下用因式分解法求解的差分表达式.三维SEG/EAGE盐丘模型的理论走时计算和积分法叠前深度偏移的实践都验证了本文方法的正确性.  相似文献   

15.
For pre‐stack phase‐shift migration in homogeneous isotropic media, the offset‐midpoint travel time is represented by the double‐square‐root equation. The travel time as a function of offset and midpoint resembles the shape of Cheops’ pyramid. This is also valid for transversely isotropic media with a vertical symmetry axis. In this study, we extend the offset‐midpoint travel‐time pyramid to the case of 2D transversely isotropic media with a tilted symmetry axis. The P‐wave analytical travel‐time pyramid is derived under the assumption of weak anelliptical property of the tilted transverse isotropy media. The travel‐time equation for the dip‐constrained transversely isotropic model is obtained from the depth‐domain travel‐time pyramid. The potential applications of the derived offset‐midpoint travel‐time equation include pre‐stack Kirchhoff migration, anisotropic parameter estimation, and travel‐time calculation in transversely isotropic media with a tilted symmetry axis.  相似文献   

16.
17.
三维VTI介质中波动方程深度偏移的最优分裂Fourier方法   总被引:17,自引:5,他引:12       下载免费PDF全文
从含Thomsen各向异性参数的qP波相速度表示式出发,建立并求解三维VTI介质中的频散方程,得到三维VTI介质中的相移算子,进而将以相移算子为基础的最优分裂Fourier方法推广到三维VTI介质,发展了一个三维VTI介质的深度偏移方法.文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致,将各向同性、弱各向异性及强各向异性统一在一个模型中.文中提出的偏移算法对相移法引入了高阶校正项来补偿介质横向变化的影响,使该方法可应用于横向非均匀VTI介质的陡角度成像,文中给出的偏移脉冲响应很好地证明了这一点.  相似文献   

18.
韩复兴  孙建国  王坤 《地球物理学报》2018,61(11):4558-4567
本文针对射线类偏移成像当中的速度模型光滑处理问题,借鉴数字图像处理当中的偏微分方程法,基于能量泛函,应用变分方法导出基于速度模型的偏微分方程实现射线类偏移成像当中的速度模型的光滑处理.由于偏微分方程法具有线性叠加特性、模型解的唯一性和局部特征保持性,因此,应用该算法可以实现基于原始速度模型空间结构的模型光滑处理.通过在原始速度模型以及光滑处理后的速度模型上计算速度的空间分布以及地震波走时、射线路径可以得出,偏微分方程法对速度模型的光滑处理能够很好地保持原始模型的空间结构,偏移成像结果也证明了该方法的实用性.  相似文献   

19.
由所建立的三维qP波相速度表示式出发,导出并解析求解各向异性介质中的频散方程,得到三维各向异性介质中的相移算子,进而将以相移算子为基础的对称非平稳相移方法推广到各向异性介质,发展了一个三维各向异性介质的深度偏移方法. 文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致,将各向同性、弱各向异性及强各向异性统一在一个模型中. 所建立的各向异性介质对称非平稳相移波场延拓算子可以同时适应速度及各向异性参数横向变化;文中给出的算例虽然是针对二维VTI介质的,但所提出的算法同样适用于三维TI介质.  相似文献   

20.
双平方根单程波动方程叠前τ偏移方法   总被引:1,自引:0,他引:1       下载免费PDF全文
本文将常规双平方根(DSR)单程波动方程从深度域变换到双程垂直走时(τ)域,由此推导出可从数学上实现“沉降观测”的单程波DSR传播算子. 其递归波场延拓算法包含波数域针对常速背景的相移处理和空间域针对横向速度扰动的相位校正,可以应对上覆地层速度横向变化对构造成像的影响. 结合零炮检距、零时间成像条件,提出了在τ域进行波场延拓与成像的DSR方程叠前偏移新方法. 为了克服其全三维偏移算法在实际应用中可能面临的困难,本文采用稳相近似,在crossline常炮检距偏移理论基础上推导了实用的共方位角叠前τ偏移方法. 数值试验表明,DSR方程叠前τ偏移在强横向非均匀介质中的成像精度与分辨率优于传统的时间域成像技术.  相似文献   

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

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