共查询到18条相似文献，搜索用时 187 毫秒
1.
叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践. 相似文献
2.
3.
4.
针对地震偏移算法中单程波算子的特征函数近似，采用最优可分表示法将该特征函数展开为空间变量(g)和水平波数(k)的可分表达式，以此可分近似表达式为基础，运用正反Fourier变换重新构造单程波算子. 为了克服特征函数最优可分近似计算在奇点及其领域产生的数值振荡，引入等价黏性技巧以增强算法的数值计算稳定性，并采用分频最优可分及对空间变量(g)线性插值的方法，不仅提高了计算精度，也节约计算机时. 文中具体研究了单程波算子的脉冲响应和二维叠前深度偏移. 结果表明，在不甚大步长情况下，本文构造的算子具有适应横向强变速的能力，运用Marmousi模型验证了本文方法适合复杂构造的成像. 相似文献
5.
研究不同偏移方法的成像机理是实现复杂构造高精度成像的前提,研究了两类共炮点域的相关型叠前深度偏移成像方法：基于单程波动方程的叠前深度成像方法和基于双程波动方程的叠前深度成像方法,同时对比了它们在计算效率、数据存储量、成像精度、成像机理、速度敏感性等方面的差异及其共性。以复杂构造模型为例,采用了傅里叶有限差分法（FFD）和逆时成像法（RTM）,这两种方法实现了121个共炮点道集的叠前深度偏移成像处理。计算结果表明,当速度准确时,两种深度域波动方程成像方法均可以恢复出各个地质反射界面,其中逆时偏移对陡倾角成像效果显著,当速度存在百分比误差和随机扰动情况时,逆时成像结果要差于单程波方法,因此,逆时偏移方法对速度的敏感性较大,且低频噪声较为严重。 相似文献
6.
目前常用叠前深度偏移主要有四种:基于射线理论的Kirchhoff偏移和高斯束偏移以及基于波动方程理论的单程波偏移和逆时偏移.要想在实际应用中取得好的成像效果,必须选取合理、有效的叠前深度偏移算法.本文首先简单阐述了这四种常用叠前深度偏移方法的基本原理,对比分析了四种方法在陡倾角、速度依赖性、起伏地表、非规则数据、低信噪比数据等方面的优劣性以及计算效率和计算精度的差异,然后利用模型试算和实际资料处理分析说明这四种叠前深度偏移方法的成像效果,最后给出不同情况下叠前深度偏移方法的选取建议. 相似文献
7.
任意介质中的动态规划法地震波三维走时计算 总被引：1，自引：1，他引：1
任意介质中的地震波三维走时计算是复杂介质情况下Kirchhoff积分法三维叠前深度偏移及走时层析成像的核心.走时算法的效率及精度决定了成像方法的应用范围及效果,对复杂地质构造区域的地震波成像时需要有稳健的走时计算方法.本文把Schneider等提出的用动态规划法计算二维任意复杂介质中走时的方法推广到三维.此方法的核心是构造从源点到当前计算点的平均慢度,基于Fermat原理,用球面波近似导出走时计算所用的公式,并用动态规划法搜索到达当前计算点的初至走时.它适用于任意复杂的介质情况,对速度差异没有限制,计算过程中考虑到各个可能的方向到达当前计算点的初至波.首波及回转波的初至走时也能正确地计算出来.各种理论速度模型上的走时计算及胜利油田某探区的三维叠前深度偏移的成功实践验证了方法的正确性. 相似文献
8.
本文将常规双平方根（DSR）单程波动方程从深度域变换到双程垂直走时（τ）域，由此推导出可从数学上实现“沉降观测”的单程波DSR传播算子. 其递归波场延拓算法包含波数域针对常速背景的相移处理和空间域针对横向速度扰动的相位校正，可以应对上覆地层速度横向变化对构造成像的影响. 结合零炮检距、零时间成像条件，提出了在τ域进行波场延拓与成像的DSR方程叠前偏移新方法. 为了克服其全三维偏移算法在实际应用中可能面临的困难，本文采用稳相近似，在crossline常炮检距偏移理论基础上推导了实用的共方位角叠前τ偏移方法. 数值试验表明，DSR方程叠前τ偏移在强横向非均匀介质中的成像精度与分辨率优于传统的时间域成像技术. 相似文献
9.
本文将常规双平方根（DSR）单程波动方程从深度域变换到双程垂直走时（τ）域，由此推导出可从数学上实现“沉降观测”的单程波DSR传播算子. 其递归波场延拓算法包含波数域针对常速背景的相移处理和空间域针对横向速度扰动的相位校正，可以应对上覆地层速度横向变化对构造成像的影响. 结合零炮检距、零时间成像条件，提出了在τ域进行波场延拓与成像的DSR方程叠前偏移新方法. 为了克服其全三维偏移算法在实际应用中可能面临的困难，本文采用稳相近似，在crossline常炮检距偏移理论基础上推导了实用的共方位角叠前τ偏移方法. 数值试验表明，DSR方程叠前τ偏移在强横向非均匀介质中的成像精度与分辨率优于传统的时间域成像技术. 相似文献
10.
11.
本文提出了一种新的基于平均导数优化方法(average-derivative optimal method,简称ADM)的二维VTI介质qP波波动方程频率空间域二阶9点格式,这种新算法将二维VTI介质qP波波动方程中中心空间导数项的差分近似表示为正交方向上3个网格点的加权平均形式.通过最小二乘优化方法求取空间导数项和加速度项的加权优化系数从而使数值频散达到极小化,每个波长所需要的网格点数在1%的误差范围内仅为3.57个网格点数,而VTI介质常规9点差分格式在相同的误差范围内则需要约12个网格点数,新方法的计算精度明显提高.复杂BP2007 2D VTI海洋标准模型数值模拟结果也验证了本文VTI介质9点ADM算法的有效性和准确性. 相似文献
12.
发展了波动方程深度延拓的频率相关变步长深度延拓方法和表驱动的单点波场插值技术.前者通过减少深度延拓的次数减少了波动方程深度偏移的计算量,而后者用很少的计算量实现了等间距、理想采样的深度成像.就同一偏移方法,采用频率相关变步长深度延拓加单点插值,其计算量大约是常规的等间距采样延拓方法的三分之一,但两者的成像效果基本相同.文中以最优分裂Fourier方法为例,用二维理论数据（Marmousi模型）和三维实际地震资料验证了这一方法,但这一方法可适用于各类频率域波动方程深度偏移方法. 相似文献
13.
在数值模拟中,隐式有限差分具有较高的精度和稳定性.然而,传统隐式有限差分算法大多由于需要求解大型矩阵方程而存在计算效率偏低的局限性.本文针对一阶速度-应力弹性波方程,构建了一种优化隐式交错网格有限差分格式,然后将改进格式由时间-空间域转换为时间-波数域,利用二范数原理建立目标函数,再利用模拟退火法求取优化系数.通过对均匀模型以及复杂介质模型进行一阶速度-应力弹性波方程数值模拟所得单炮记录、波场快照分析表明:这种优化隐式交错网格差分算法与传统的几种显式和隐式交错网格有限差分算法相比不但降低了计算量,而且能有效的压制网格频散,使弹性波数值模拟的精度得到有效的提高. 相似文献
14.
各向异性介质纯P波方程完全不受横波的干扰,在一定程度上可以减缓由于介质各向异性引起的数值不稳定,本文推导了具有垂直对称轴的横向各向同性(VTI)介质纯P波一阶速度-应力方程.由于纯P波方程存在一个分数形式的伪微分算子,无法直接采用有限差分法求解.针对该问题,本文采用伪谱法和高阶有限差分法联合求解波动方程,重点分析了混合法求解纯P波一阶速度-应力方程的稳定性问题,并给出了混合法求解纯P波方程的稳定性条件.数值模拟结果表明纯P波方程伪谱法和高阶有限差分混合法能够进行复杂介质的正演模拟,在强变速度、变密度的地球介质中仍然具有较好的稳定性. 相似文献
15.
16.
本文旨在构建适用于滞弹性近场时域波动有限元模拟的高精度人工边界条件：完美匹配层（Perfectly Matched Layer：PML）,其中阻尼介质时域本构基于广义标准线性体建立.与以往研究不同,本文采用复坐标延拓技术变换弱形式波动方程构建了可直接用有限元离散的弱形式时域PML,规避以往独立对无限域内波动方程及界面条件进行延拓可导致的PML场方程和界面条件匹配不合理引发数值失稳、计算精度低下等问题.其次,针对PML中多极点有理分式与频域函数乘积的傅里叶反变换难以计算的问题,利用PML精度对复坐标延拓函数中延拓参数微调不敏感这一特点,明确给出了参数微调准则以规避多重极点,进而利用有理分式分解给出了一种普适、简便的计算方法,极大地简化了PML计算.基于该方法可实现任意高阶PML.最后,将本文构建滞弹性PML与高阶勒让德谱元（高精度集中质量有限元）结合得到滞弹性近场波动谱元离散方案.基于算例验证了滞弹性PML的计算效率、精度及新离散方案的长持时稳定特性.新离散方案可应用于计入实际介质阻尼的地震波动正、反问题数值模拟,提高波形模拟的精度以及地下波速结构反演的精度和可靠性. 相似文献
17.
海洋勘探环境可以抽象为下伏固体与上覆流体相互耦合的介质,本文针对流-固边界耦合介质提出了一种高效、稳定的多参数（速度和密度）全波形反演方法.本文采用弹性波一阶位移-应力方程作为过渡层耦合声波压力方程与弹性波位移方程来模拟耦合环境,相比于传统的交错网格建模方法或者构建连续性条件,本文提出的方法在正演精度和稳定性上凸显出很大优势,极大降低了计算内存.反演策略对多参数全波形反演至关重要,由于不同参数之间的相互耦合使得密度在多参数全波形反演中较难获得,因此本文将非均匀流-固边界耦合介质多参数全波形反演分为两个步骤完成：第一步利用变密度声波方程结合推导出的密度梯度算子进行纵波速度和密度的双参数反演;第二步根据链式法则求取横波速度的梯度,结合第一步的反演结果使用流-固边界耦合方程反演横波速度.最后通过与声波动方程数值模拟结果对比证明正演算法的准确性;上覆流体的Marmousi-2模型的数值试验测试说明反演方法的有效性和适应性. 相似文献
18.
应用高精度人工边界条件可有效提升近场波动数值模拟计算效率.完美匹配层是吸收层形式高精度人工边界条件,匹配层内场方程和界面条件通常分别采用复坐标延伸技术变换强形式无限域内波动方程和界面条件得到,亦曾将无限域界面条件当作匹配层界面条件.场方程和界面条件构建过程相互独立,可能出现匹配不合理而引发数值失稳、计算精度低下等问题.本文提出采用复坐标延伸技术变换弱形式无限域波动方程以构建完美匹配层的方法.弱形式波动方程耦合了波动方程及界面条件,进而规避了变换后所得场方程与界面条件之间的匹配不合理问题.新方法可直接建立弱形式匹配层,在此基础上亦可给出强形式匹配层.弱形式便于有限元离散,强形式便于有限差分离散.基于弱形式完美匹配层,结合勒让德谱元建立了弹性介质近场波动谱元模拟方案.利用算例验证了新方案的精度及数值稳定性.本文工作可直接推广至多相耦合介质近场波动数值模拟. 相似文献