首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 109 毫秒
1.
单程波近似实际上是一种多次前向散射和单次后向散射近似.利用单程波近似来描述波传播可以极大地节省地震数值模拟的计算时间和内存,实现地震波长距离传播模拟和三维地震模拟快速计算.本文基于单程波近似和波动积分方程的分离变量逼近,从广义Lippmann-Schwinger波动积分方程推导出耦合反射/透射系数的单程波传播算子.该算子由两部分构成:分离变量Fourier单程波传播算子和薄板间的反射/透射系数表达.前者将常规的Fourier分裂步单程波传播算子(SSF)推广适应横向强速度变化介质和大角度传播波场.后者是利用垂直波数来表示反射/透射系数,自然耦合到波场传播的计算过程中,其为地质界面倾角的隐式表达,精确描述振幅随入射角的变化,能适应任意复杂的模型.通过两个数值算例和一个实际地质模型的计算,本文将该方法和边界元法进行了比较,结果表明:在算例给出的介质横向速度变化情况下,本文提出的方法在相位和振幅方面与全波数值方法基本吻合.  相似文献   

2.
快速Fourier变换波动方程基准面校正方法研究   总被引:4,自引:1,他引:3       下载免费PDF全文
当地表起伏剧烈、近地表速度横向变化较大时,基于地表一致性假设的常规静校正方法存在着较大误差.波动方程基准面静校正方法能很好地解决起伏地表和复杂近地表结构问题,但计算量巨大,特别是三维波动方程基准面校正,适应横向任意速度变化、计算精度较高的有限差分或其混合的方法波动方程基准面校正涉及海量的计算和存储操作.为了提高波动方程基准面校正的计算效率,本文研究一类只用快速Fourier变换(FFT)实施波动方程基准面校正的方法,采用相移(PS)、分裂步(SSF)和一阶退化(DP1)三种具有相同算法结构、但不同计算效率、适应不同地表复杂程度的Fourier变换延拓算子.PS和SSF算子只适应于速度横向变化较弱的起伏地表;DP1通过在两个分裂步之间作波数域线性插值来实现波场延拓,将常规的SSF算法推广适应强速度横向变化介质和大角度传播波场.本文着重比较了基于这三种延拓算子的逐层延拓累加波动方程基准面校正方法对地表起伏和近地表速度横向变化的适应能力和计算效率,给出了一个相对定量的评估,以便针对不同的地表复杂程度合理选择合适的FFT波动方程基准面校正方法,既满足了精度又提高了计算效率.  相似文献   

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

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

5.
共炮检距道集波动方程保幅叠前深度偏移方法   总被引:9,自引:9,他引:0       下载免费PDF全文
本文提出了一种基于双平方根算子的共炮检距道集波动方程保幅叠前深度偏移方法,将振幅误差补偿作为偏移的一部分与“运动学偏移”一起在偏移过程中实现.其基本内容包括:(1)从保幅的单平方根算子方程出发,推导出由双平方根算子定义的保幅单程波方程;(2)根据地震波摄动理论把速度场分裂为层内常速背景和变速扰动,分别在频率-波数域和频率-空间域求得波场深度延拓的偏移时移量及振幅校正系数,从而得到最终的DSR保幅波场延拓算子;(3)在高频假设条件下,把DSR保幅波场延拓公式中的积分运算进行稳相近似,得到保幅波场延拓的相移公式.理论分析和模型数值试验表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析(如AVO/AVA)提供更真实的地震信息.  相似文献   

6.
三维非均匀介质中真振幅地震偏移算子研究   总被引:25,自引:6,他引:25       下载免费PDF全文
利用三维非均匀介质中的波动方程,进行振幅保真波场偏移算子分解,得到用于真振幅偏移的单程波方程. 经过数学推理,导出裂步Fourier法真振幅偏移和Fourier有限差分法真振幅偏移的算子方程,并给出其具体的实现过程.  相似文献   

7.
尤加春  曹俊兴  王俊 《地球物理学报》1954,63(10):3838-3848
叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践.  相似文献   

8.
傅里叶有限差分法保幅叠前深度偏移方法   总被引:23,自引:5,他引:23       下载免费PDF全文
地震数据中饱含有丰富的走时信息和振幅信息. 为解决传统偏移方法中几何扩散和入射角变化引起的振幅误差问题,本文提出了一种实用的波动方程保幅地震偏移方法. 该方法从全声波方程出发进行单程波保幅分解,得到直观、高效率的直接面对地震波传播波场的压力分量进行延拓的保幅偏移单程波方程,进而推导出一个含有6项的傅里叶有限差分法保幅偏移的算子方程;修改边界条件和成像条件,使修改后的边界条件和成像方程中考虑振幅补偿,从而从三方面补偿几何扩散损失和入射角变化对振幅的影响. 脉冲响应测试、单炮记录的数值试验以及SEG/EAGE盐丘模型的叠前偏移结果表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析(如AVO/AVA)提供更真实的地震信息.  相似文献   

9.
单程波算子积分解的象征表示   总被引:7,自引:5,他引:7       下载免费PDF全文
单程波波场延拓算子在地震偏移成像中有重要应用.单程波波场延拓算子按其实现方式可分为Kirchhoff积分、空间隐式有限差分和Fourier变换方法,他们代表了算子的不同表示方法,当截断使用这些方法时会得到不同的精度.象征表示对这些方法的导出和精度分析有重要作用.算子作用于正弦波函数所得函数称为算子的象征.算子的象征是褶积算子Fourier变换的推广.Fourier变换方法则直接用象征函数的可分表示求出.空间隐式有限差分则可以用象征函数的Padè近似或部分分式导出.单程波算子在深度域的积分称为单程波算子积分解.本文推导了单程波算子积分解的象征表达式,给出了算子象征的代数运算的头几阶表达式,这些表达式还未在前人文献中发现.Kirchhoff积分所需格林函数可以通过象征函数和鞍点法导出.基于积分解的象征表达式给出了非对称走时公式,对改善Kirchhoff积分的聚焦性能有重要意义.  相似文献   

10.
基于相位编码的Fourier变换波动方程三维地震照明分析   总被引:1,自引:1,他引:0  
复杂构造勘探是当前地震勘探的主要目标,提高深层信噪比和分辨率十分重要.地震照明分析是研究地震波能量在复杂构造中传播的有效手段,可指导观测系统设计进而提高野外数据采集和最终成像质量.本文发展了基于三维单程波动方程的地震照明方法,采用相位编码合成平面波源,一阶分离变量Fourier算子进行波场传播.阐述了照明和深度偏移之间...  相似文献   

11.
混合域单程波传播算子及其在偏移成像中的应用   总被引:6,自引:4,他引:2  
以地震波的单程波传播方程为基础,利用算子近似展开的方法推导出了当前波动方程叠前深度偏移方法研究中广泛使用的裂步Fourier、Fourier有限差分和广义屏传播算子的一般形式及其近似式.讨论了它们间差异、相互关系以及他们的特点,最后给出了基于裂步Fourier、Fourier有限差分和广义屏传播算子的偏移成像方法时Marmousi模型的偏移成像结果,以说明它们间的优劣与计算效率.  相似文献   

12.
— Dual-domain one-way propagators implement wave propagation in heterogeneous media in mixed domains (space-wavenumber domains). One-way propagators neglect wave reverberations between heterogeneities but correctly handle the forward multiple-scattering including focusing/defocusing, diffraction, refraction and interference of waves. The algorithm shuttles between space-domain and wavenumber-domain using FFT, and the operations in the two domains are self-adaptive to the complexity of the media. The method makes the best use of the operations in each domain, resulting in efficient and accurate propagators. Due to recent progress, new versions of dual-domain methods overcame some limitations of the classical dual-domain methods (phase-screen or split-step Fourier methods) and can propagate large-angle waves quite accurately in media with strong velocity contrasts. These methods can deliver superior image quality (high resolution/high fidelity) for complex subsurface structures. One-way and one-return (De Wolf approximation) propagators can be also applied to wave-field modeling and simulations for some geophysical problems. In the article, a historical review and theoretical analysis of the Born, Rytov, and De Wolf approximations are given. A review on classical phase-screen or split-step Fourier methods is also given, followed by a summary and analysis of the new dual-domain propagators. The applications of the new propagators to seismic imaging and modeling are reviewed with several examples. For seismic imaging, the advantages and limitations of the traditional Kirchhoff migration and time-space domain finite-difference migration, when applied to 3-D complicated structures, are first analyzed. Then the special features, and applications of the new dual-domain methods are presented. Three versions of GSP (generalized screen propagators), the hybrid pseudo-screen, the wide-angle Padé-screen, and the higher-order generalized screen propagators are discussed. Recent progress also makes it possible to use the dual-domain propagators for modeling elastic reflections for complex structures and long-range propagations of crustal guided waves. Examples of 2-D and 3-D imaging and modeling using GSP methods are given.  相似文献   

13.
The improvement in accuracy and efficiency of wave-equation migration techniques is an ongoing topic of research. The main problem is the correct imaging of steeply dipping reflectors in media with strong lateral velocity variations. We propose an improved migration method which is based on cascading phase-shift and finite-difference operators for downward continuation. Due to these cascaded operators we call this method‘Fourier finite-difference migration’(FFD migration). In our approach we try to generalize and improve the split-step Fourier migration method for strong lateral velocity variations using an additional finite-difference correction term. Like most of the current migration methods in use today, our method is based on the one-way wave equation. It is solved by first applying the square-root operator but using a constant velocity at each depth step which has to be the minimum velocity. In a second step, the approximate difference between the correct square-root operator and this constant-velocity squareroot operator (the error made in the first step) is implemented as an implicit FD migration scheme, part of which is the split-step Fourier correction term. Some practical aspects of the new FFD method are discussed. Its performance is compared with that of split-step and standard FD migration schemes. First applications to synthetic and real data sets are presented. They show that the superiority of FFD migration becomes evident by migrating steeply dipping reflectors with complex overburden having strong lateral velocity variations. If velocity is laterally constant, FFD migration has the accuracy of the phase-shift method. The maximum migration angle is velocity adaptive, in contrast to conventional FD migration schemes. It varies laterally depending on the local level of velocity variation. FFD migration is more efficient than higher-order implicit FD schemes. These schemes use two cascaded downward-continuation steps in order to attain comparable migration performance.  相似文献   

14.
高阶方程偏移的分裂算法   总被引:23,自引:6,他引:23       下载免费PDF全文
对于用高阶偏微分方程进行地震偏移,本文提出一种分裂算法。有限差分法目前只能解波动方程的二阶和三阶的近似方程的偏移问题,更高阶的近似方程还没有一种有效的解法。利用高阶方程进行偏移会提高偏移的精度和效果。本文就是为了解决这个问题所进行的探索。 本文阐述了分裂算法的原理,并给出了应用该算法解决偏移问题的实例。  相似文献   

15.
作为高精度波形反演或逆时偏移的重要组成部分,地震波数值模拟对计算速度和效率提出了更高要求.GPU通用计算技术的产生及其内在数据并行性,为高效地震波数值模拟应用和研究得以有效开展奠定了基础.本文借助交错网格的Fourier伪谱微分矩阵算子和GPU上高效矩阵乘法,实现了复杂介质中地震波模拟的高效算法.数值试验表明,优化后的GPU计算相比CPU单核计算在大规模二维地震波场计算中获得至少100x以上的加速比.这对我们快速分析目标反射层在地震剖面中同相轴位置,制定优化采集方案具有重要意义.  相似文献   

16.
Various migration methods have been proposed to image high-angle geological structures and media with strong lateral velocity variations; however, the problems of low precision and high computational cost remain unresolved. To describe the seismic wave propagation in media with lateral velocity variations and to image high-angle structures, we propose the generalized screen propagator based on particle swarm optimization (PSO-GSP), for the precise fitting of the single-square-root operator. We use the 2D SEG/EAGE salt model to test the proposed PSO-GSP migration method to image the faults beneath the salt dome and compare the results to those of the conventional high-order generalized screen propagator (GSP) migration and split-step Fourier (SSF) migration. Moreover, we use 2D marine data from the South China Sea to show that the PSO-GSP migration can better image strong reflectors than conventional imaging methods.  相似文献   

17.
波动方程深度偏移的频率相关变步长延拓方法   总被引:9,自引:1,他引:8       下载免费PDF全文
发展了波动方程深度延拓的频率相关变步长深度延拓方法和表驱动的单点波场插值技术.前者通过减少深度延拓的次数减少了波动方程深度偏移的计算量,而后者用很少的计算量实现了等间距、理想采样的深度成像.就同一偏移方法,采用频率相关变步长深度延拓加单点插值,其计算量大约是常规的等间距采样延拓方法的三分之一,但两者的成像效果基本相同.文中以最优分裂Fourier方法为例,用二维理论数据(Marmousi模型)和三维实际地震资料验证了这一方法,但这一方法可适用于各类频率域波动方程深度偏移方法.  相似文献   

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

19.
库车坳陷复杂高陡构造地震成像研究   总被引:2,自引:1,他引:1       下载免费PDF全文
复杂构造地震成像主要取决于叠前地震数据品质、偏移速度可靠性和偏移算子成像精度. 库车坳陷异常复杂的近地表条件导致极低信噪比的地震采集数据. 该区逆冲推覆高陡构造刺穿盐体大面积分布, 盐层厚度变化大、顶底面形态复杂, 盐下断裂带破碎、小断块发育, 形成异常复杂的地震成像问题. 本文重点研究三个关键环节:(1)精细的叠前地震预处理研究: 根据该区地震地质复杂性和地震资料特征, 采用一些新的方法技术和技术组合从振幅与时移的大、中、小尺度变化三个层次来解决资料信噪比问题, 重建深部反射信号; (2)三级偏移速度分析研究:利用库车坳陷盐刺穿逆冲推覆构造建模理论及变速成图配套技术解决叠前时间偏移速度场时深转换问题,利用井约束低频速度地震迭代反演技术解决连井层速度场与偏移速度场的融合问题,实现从DMO速度分析、叠前时间偏移速度分析到叠前深度偏移速度分析的有机衔接,建立拓扑结构相对保持的叠前深度偏移速度模型;(3)基于退化Fourier偏移算子的半解析波动方程叠前时间和深度偏移研究, 极大地改善了地震偏移过程中高波数波的成像问题. 通过对库车坳陷大北、博孜、却勒、西秋4和西秋10等复杂高陡构造的叠前时间和深度偏移地震成像处理,取得了较好的应用效果.  相似文献   

20.
波动方程偏移成像阴影的照明补偿   总被引:6,自引:6,他引:6       下载免费PDF全文
受地下复杂构造和地震数据采集系统的影响,使地震波对地下目标的照明出现不均匀性,地震采集系统难以有效地获取地下某些目标的反射信息,进而使数据偏移成像在这些目标体上出现成像阴影. 根据波场和Green函数的窗口Fourier框架展开,利用角度域波动方程偏移成像和波动方程照明分析,并结合波动方程反演理论,提出一种角度域波动方程偏移成像阴影照明补偿方法. 这种补偿方法能同时考虑地震数据采集系统和波场传播路径对偏移成像的影响,消除复杂构造区的偏移成像阴影,改进波动方程叠前深度偏移成像在复杂构造区的成像效果.  相似文献   

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

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