基于三角网格的有限差分法叠后逆时偏移   总被引：1，自引：1，他引：0  

郭书娟  李振春  孙小东  叶月明  滕厚华  李芳 《应用地球物理》2008,5(2):115-120

Compared with other migration methods, reverse-time migration is based on a precise wave equation, not an approximation, and performs extrapolation in the depth domain rather than the time domain. It is highly accurate and not affected by strong subsurface structure complexity and horizontal velocity variations. The difference method based on triangular grids maintains the simplicity of the difference method and the precision of the finite element method. It can be used directly for forward modeling on models with complex top surfaces and migration without statics preprocessing. We apply a finite difference method based on triangular grids for post-stack reverse-time migration for the first time. Tests on model data verify that the combination of the two methods can achieve near-perfect results in application.  相似文献   

格子法在起伏地表叠前逆时深度偏移中的应用   总被引：6，自引：2，他引：4       下载免费PDF全文

徐义 《地球物理学进展》2008,23(3):839-845

基于全程波波动方程的逆时偏移(Reverse Time Migration)可以对回转波、多次反射波成像,不受横向速度变化影响,没有倾角限制,随着计算机软硬件技术的进步,再次成为偏移方法研究热点.本文将格子法用于叠前逆时深度偏移成像.格子法作为波场延拓方法,处理起伏地表边界条件容易,可用于含起伏地表边界条件的逆时波场延拓;可利用变尺度非规则对计算域进行离散,因此可根据速度模型调整网格尺度来降低存储量,放大时间步长,降低计算量.采用光滑的曲人工边界,也可避免常规的PML吸收边界存在的角点区域需特别处理的麻烦.本方法通过事先计算和存储边界单元的局部几何参数,与直边界PML方法相比不增加任何计算量.格子法还具有容易实现并行计算的特点,非常适用于叠前逆时偏移.本文给出了二维问题算例.  相似文献   

地震波逆时偏移中两种成像条件应用效果对比     

宋宗平  陈可洋  杨微  李来林  吴清岭  范兴才 《物探与化探》2019,(3):618-625

成像条件是实现地震波场准确成像的关键环节,而优化逆时成像条件在改善逆时成像应用效果中的作用尤为明显。为此,从零井源距和非零井源距脉冲响应出发,系统对比了常规相关法逆时成像条件和基于行波分离法逆时成像条件的逆时偏移脉冲响应波场差异,分析了逆时偏移低波数噪声的成因机理及波场规律,同时以复杂断陷理论模型和SZ工区实际三维地震资料为例,系统对比了两种逆时成像条件的成像应用效果。研究结果表明,采用基于行波分离法逆时成像条件的成像结果能够有效衰减逆时偏移低波数背景噪声,恢复掩盖在这种背景噪声之下的有效地震成像信息,地层细节刻画更加清晰,这可为当前复杂波场和复杂构造高精度地震成像提供方法指导。  相似文献   

变时间步长有限差分法逆时偏移   总被引：3，自引：2，他引：1  

董春晖  王尚旭 《地球物理学进展》2010,25(4)

逆时偏移因其可以利用多次波和回转波等信息,不受倾角限制对构造进行成像而逐渐成为一个研究的热点.但相对极低的计算效率,成为限制其实际应用的瓶颈之一.针对该问题,本文提出了一种时变的计算范围及可变时间步长的思路,即在不同的时间段内利用不同的时间步长计算相应空间范围内的波场,从而在算法上节省大量的计算时间.该思路简单易行,本文将其应用于有限差分法逆时偏移;通过对理论模型数据进行处理,证明该方案不影响波场计算的精度,并且可以在一定程度上提高计算的效率.  相似文献   

多分量联合逆时偏移最佳匹配层吸收边界   总被引：3，自引：2，他引：1  

杜启振  秦童  朱钇同  毕丽飞 《应用地球物理》2010,7(2):166-173

有限空间内的波动方程逆时偏移需要利用有效的边界处理技术用以消除人工截断对偏移结果产生的影响。本文以横向各向同性介质弹性波速度-应力方程为基础,依据传统分裂式最佳匹配层（Perfect Matched Layer,PML）吸收边界技术的思想,推导了应用于逆时偏移的完全匹配层波动方程,并给出了其高阶交错网格有限差分格式。针对由边界处向计算区域内传播的＂反射波＂,以及地震记录排列两端地震同相轴突变对计算区域的影响这两方面问题,本文给出了逆时偏移中吸收层的布设方式。模型和实际资料的弹性波叠前多分量联合逆时深度偏移结果表明本文的边界处理方法取得了较好的吸收效果,获得了好的联合偏移成像结果。  相似文献   

地震叠前逆时偏移算法的CPU/GPU实施对策   总被引：16，自引：8，他引：8       下载免费PDF全文

李博  刘红伟  刘国峰  佟小龙  刘洪  郭建  裴江云 《地球物理学报》2010,53(12):2938-2943

相较于单程波偏移算法而言,逆时偏移成像方法以其物理基础为依托优势,几十年来一直备受国内外地球物理学家的青睐.目前的逆时偏移(RTM)若直接采用双程波动方程进行延拓,尽管可以回避上下行波的分离处理,然就已有算法而言,其计算量和I/O(输入/输出)量却是最大的.针对此问题,本文在分析现行逆时偏移的多种算法基础上,提出利用CPU/GPU(中央处理器/图形处理器)作为数值计算核心,建立随机边界模型,从而克服存储I/O难题和提高计算效率.在实际的数据测试中,本文的方法可以大幅度的提高计算效率和减少存储单元,从而促使其高效地应用于生产实际.  相似文献   

基于阈值混合范数的黏声最小二乘逆时偏移     

宋利伟  侯晓慧  石颖  柯璇  陈树民  李松龄 《地球物理学进展》2024,(1):334-343

最小二乘逆时偏移方法在重建高分辨地下构造图像方面具有理论优势.然而,常规方法忽略了介质对地震波的吸收衰减效应,且受噪声干扰易产生偏移假象,以致成像分辨率较低.为提高最小二乘逆时偏移方法服务于油气勘探开发的能力,本文将常分数阶拉普拉斯算子波动方程作为黏滞性介质中波场的控制方程以降低模拟数据和观测数据匹配的不确定度.此外,利用阈值混合范数构建最小二乘逆时偏移的目标函数以提高成像方法的抗噪性.通过合成和实际数据测试表明所提出的方法能有效地生成高分辨的构造图像.  相似文献   

逆时偏移计算中的边界处理分析及应用   总被引：3，自引：3，他引：0       下载免费PDF全文

胡昊  刘伊克  常旭  王一博  杜向东  杨仁虎 《地球物理学报》2013,56(6):2033-2042

在地震资料的处理中,逆时偏移方法可以实现复杂构造高精度成像,但计算量和存储量两大问题影响了该方法的实际应用.逆时偏移算法中的边界处理方式的优化可以大大减少存储量,本文讨论了记录波场边界信息的吸收边界和随机边界这两种方法的计算效率和成像效果,提出了吸收边界中只记录边界范围内未衰减单层波场的方法并通过数值实验验证了其可行性,选择了不同的边界策略应用于模型数据和实际资料的处理.结果表明:记录单层波场边界信息的吸收边界的成像效果同传统存储波场历史的方式几乎无差别,但要额外存储每个时刻的波场边界信息;随机边界不需要额外存储波场信息,但会带来边界漫反射影响和计算区域的增加;记录单层波场边界可以明显减少存储量,并且不影响成像效果.  相似文献   

各向异性介质qP波传播描述I：伪纯模式波动方程   总被引：6，自引：6，他引：0       下载免费PDF全文

程玖兵  康玮  王腾飞 《地球物理学报》2013,56(10):3474-3486

地球介质相对于地震波波长尺度的定向非均匀性会导致波速的各向异性,进而影响地震波场的运动学与动力学特征．各向异性弹性波动方程是描述该类介质波场传播的基本工具,在正演模拟、偏移成像与参数反演中起着关键作用．为了面向实际应用构建灵活、简便的各向异性波场传播算子,人们一直在寻求简化的各向异性波动方程．本文借鉴各向异性弹性波波型分离思想,通过对平面波形式的弹性波方程(即Christoffel方程)实施一种代表向波矢量方向投影的相似变换,推导出了一种适应任意各向异性介质、运动学上与原始弹性波方程完全等价,在动力学上突出qP波的新方程,即qP波伪纯模式波动方程．文中以横向各向同性(TI)介质为例,给出了相应的qP波伪纯模式波动方程及其声学与各向同性近似,并在此基础上开展了正演模拟和逆时偏移试验,展示了这种描述各向异性波场传播的新方程的特点与优势．  相似文献   

弹性波逆时偏移中的稳定激发振幅成像条件   总被引：6，自引：5，他引：1       下载免费PDF全文

张智  刘有山  徐涛  白志明 《地球物理学报》2013,56(10):3523-3533

本文针对弹性波逆时偏移,提出稳定的激发振幅成像条件.在震源波场的正向传播过程中,计算每个网格点的能量,并保存最大能量密度的时刻和相应的波场值;在检波器波场的逆时传播过程中,在每个网格点提取最大能量密度时刻的检波器波场值,并利用保存的最大能量震源波场做归一化,获得角度依赖的反射系数成像剖面.相比于归一化互相关成像条件,该成像条件在震源波场的正向传播过程中无需存储波场快照,节省大量磁盘空间和I/O吞吐任务,提高了计算效率;相比于弹性波的激发时间成像条件,该成像条件自动校正了水平分量在震源两侧的极性反转,在多炮叠加时避免振幅损失.数值试验表明,与归一化成像条件相比,稳定激发振幅成像条件具有更小的计算量,偏移剖面的低频假象更弱,水平分量的成像能力更优,具有更高的空间分辨率.  相似文献