最小二乘逆时偏移(least-squares reverse time migration,LSRTM)具有更高的成像分辨率、振幅保真性及均衡性等优势。然而,目前的LSRTM算法大多基于水平地表假设,在面对复杂地形时无法很好地适应剧烈的起伏地表。基于二阶常密度声波方程...最小二乘逆时偏移(least-squares reverse time migration,LSRTM)具有更高的成像分辨率、振幅保真性及均衡性等优势。然而,目前的LSRTM算法大多基于水平地表假设,在面对复杂地形时无法很好地适应剧烈的起伏地表。基于二阶常密度声波方程的LSRTM算法忽略了密度变化对振幅的影响,很难在变密度介质中取得保真的成像结果。为此,从一阶速度-应力方程出发,在曲线坐标系下推导相应的扰动方程和伴随方程,并通过伴随状态法给出梯度更新公式,最终实现基于贴体网格的起伏地表LSRTM算法。模型试算验证了算法的有效性和对复杂地表的适应性。结果表明,提出的算法能够消除起伏地表的影响、压制低频噪声、恢复高频成分、均衡成像振幅,实现地下变密度介质的高分辨率和高保真度成像。展开更多
从TTI介质一阶应力—速度方程出发,利用旋转交错网格高阶有限差分方法,将非分裂完全匹配层(Non-spliting Perfect Match Layer,简称NPML)边界吸收条件和自由边界条件相结合形成组合边界条件,进行了二维三分量TTI介质弹性波场数值模拟。...从TTI介质一阶应力—速度方程出发,利用旋转交错网格高阶有限差分方法,将非分裂完全匹配层(Non-spliting Perfect Match Layer,简称NPML)边界吸收条件和自由边界条件相结合形成组合边界条件,进行了二维三分量TTI介质弹性波场数值模拟。波场快照和炮记录表明:①采用非分裂式边界条件能较好地消除近地表大角度入射波和瞬逝波;②组合边界条件与NPML边界吸收条件相比,不仅有效地压制了边界反射,同时实现了对自由地表的模拟,获得了丰富的全波场信息,其中在地表产生的PS转换横波作为一种特殊的横波现象,可为近地表结构调查以及多波波场分析等提供有益信息;③自由地表引起的面波以及多次波对偏移结果有着重要影响,因此在实际地震资料处理中应当充分考虑自由地表条件对波场的影响效应。数值模拟结果证实了组合边界条件下二维三分量TTI介质波场数值模拟方法的可行性和正确性。展开更多
文摘最小二乘逆时偏移(least-squares reverse time migration,LSRTM)具有更高的成像分辨率、振幅保真性及均衡性等优势。然而,目前的LSRTM算法大多基于水平地表假设,在面对复杂地形时无法很好地适应剧烈的起伏地表。基于二阶常密度声波方程的LSRTM算法忽略了密度变化对振幅的影响,很难在变密度介质中取得保真的成像结果。为此,从一阶速度-应力方程出发,在曲线坐标系下推导相应的扰动方程和伴随方程,并通过伴随状态法给出梯度更新公式,最终实现基于贴体网格的起伏地表LSRTM算法。模型试算验证了算法的有效性和对复杂地表的适应性。结果表明,提出的算法能够消除起伏地表的影响、压制低频噪声、恢复高频成分、均衡成像振幅,实现地下变密度介质的高分辨率和高保真度成像。
文摘从TTI介质一阶应力—速度方程出发,利用旋转交错网格高阶有限差分方法,将非分裂完全匹配层(Non-spliting Perfect Match Layer,简称NPML)边界吸收条件和自由边界条件相结合形成组合边界条件,进行了二维三分量TTI介质弹性波场数值模拟。波场快照和炮记录表明:①采用非分裂式边界条件能较好地消除近地表大角度入射波和瞬逝波;②组合边界条件与NPML边界吸收条件相比,不仅有效地压制了边界反射,同时实现了对自由地表的模拟,获得了丰富的全波场信息,其中在地表产生的PS转换横波作为一种特殊的横波现象,可为近地表结构调查以及多波波场分析等提供有益信息;③自由地表引起的面波以及多次波对偏移结果有着重要影响,因此在实际地震资料处理中应当充分考虑自由地表条件对波场的影响效应。数值模拟结果证实了组合边界条件下二维三分量TTI介质波场数值模拟方法的可行性和正确性。
