实现复频移完全匹配层吸收边界条件的递推卷积方法

完全匹配层吸收边界（PML）已经被证明是非常有效的边界吸收技术,对体波和面波的吸收都具有非常好的效果,已经被广泛应用于弹性波的数值模拟中。但是在某些情况下传统的PML技术还是存在一定的问题,比如对掠射情况下的体波和窄区域自由表面条件下的面波的吸收等等。在坐标变换中采用复频移拉伸函数的复频移PML可以有效地改善PML边界条件的吸收性能。基于弹性波一阶速度-应力方程,推导了复频移PML的递推卷积实现方法,并采用交错网格高阶有限差分法对其进行了数值模拟,与传统的PML进行了对比。结果表明：传统的PML对掠射情况下的体波和窄区域自由表面条件下的面波吸收不足,会产生虚假反射,影响真实波场;而基于递推卷积的复频移PML算法能够有效地改善困难情况下的吸收效果,并且在实现过程中不用分裂变量,应用更加方便简单。计算卷积时采用递推的形式,推导过程直观易懂,易于编程,而且不会增加计算量,存储量也没有太大的变化。     

地震资料多次波压制方法研究进展及展望

在反射波地震勘探中,发育的多次波会严重降低地震资料处理的精度,进而误导解释人员给出错误的解释结果.因此,多次波压制方法研究一直以来都是地震资料处理中不可忽视的重要环节.目前,将多次波压制方法分为基于特征差异的滤波方法和基于波动理论的预测相减法是一种较为流行的分类方法.前者计算效率较高,但是多次波压制精度有限.后者理论先进,不需要或者极少需要地层先验信息,多次波压制精度较高,但是往往伴随高计算成本.通常,采用单一方法往往难以满足复杂实际资料处理的需求,汇二者所长的综合法则可以针对实际地震资料的特点,通过不同方法组合解决实际应用问题.据上述分类方法,首先详细论述了每一类方法中经典方法的原理、研究进展及发展方向等;然后介绍了目前实际应用中的综合法多次波压制实例;最后,在回顾传统多次波压制方法的研究进程基础上重点介绍了近年来出现的几种新方法,并展望多次波压制的发展前景.     

复杂地表条件下保幅高斯束偏移

高斯束偏移是一种准确、灵活、高效的深度域成像方法,其不但具有接近于波动方程偏移的成像精度,还保留了Kirchhoff偏移灵活、高效的特点以及对复杂地表条件良好的适应性.本文提出了一种适用于复杂地表条件的且具有相对振幅保持特点的高斯束偏移方法.通过考虑地表高程、倾角以及实际的道间距等信息,推导了基于高斯束表示的波场反向延拓公式,并结合反褶积成像条件,得到了复杂地表条件下的共炮域保幅高斯束偏移公式.同原有方法相比,本文方法不但可以直接在起伏的地表面进行局部平面波的分解,具有更高的成像精度,而且可以得到反映地下随角度变化反射系数的成像结果.数值模型的试算验证了上述结论.     

几种自由边界实施方法在完全 匹配层条件下的对比研究

传统的完全匹配层技术是一种能够较为有效地消除边界反射的边界条件,但是当表层为泊松比较高的自由表面时,该技术可能会产生不稳定的现象.针对传统的完全匹配层技术固有的不稳定和掠射情况下吸收效果不佳等缺陷,发展了多轴完全匹配层、卷积完全匹配层以及将两者结合的多轴卷积完全匹配层等3种边界条件.本文介绍了水平自由表面的不同处理方法以及传统、多轴、卷积和多轴卷积等4种完全匹配层条件的原理,通过二维半无限空间模型的交错网格有限差分正演模拟对比,分析了几种自由边界实施方法在这几种完全匹配层条件下的稳定性,并通过提取单道波形与解析解进行对比,定性分析了水平自由表面几种不同处理方法的准确性以及各自的适用条件. 结果表明,泊松比和水平自由表面实施方法对波场模拟效果及其稳定性有重要影响.     

基于稳定逆时传播算子的黏声介质最小二乘逆时偏移

基于GSLS模型黏声介质二阶拟微分方程,采用伪谱法进行数值模拟。针对黏声介质逆时传播过程中产生的高频不稳定问题,提出加入规则化算子对其进行消除的方法,构建了稳定的逆时传播算子。在最小二乘反演的基础上,将黏声介质逆时偏移与最小二乘思路相结合,发展了带有振幅补偿的黏声介质最小二乘逆时偏移( LSRTM)。 Marmousi模型结果表明:相对于常规最小二乘逆时偏移,黏声介质最小二乘逆时偏移校正了地层的黏滞性,得到了更加精确可靠的保幅成像剖面。     

基于黏声介质的角度域高斯束逆时偏移方法研究

高斯束逆时偏移是一种兼具计算效率和成像精度的深度域成像方法,能够面向目标成像.地下介质中黏滞性普遍存在,利用传统各向同性或完全弹性的成像方法处理黏滞性探区的数据会降低分辨率,并导致成像位置不准确和振幅欠估计等问题.本文在高斯束逆时偏移的基础上,通过对震源点和检波点处的波场进行衰减补偿,并结合高斯束求解时的角度信息,实现了黏声介质角度域高斯束逆时偏移方法.最后通过模型和实际资料试算对本文方法的正确性和适用性进行了验证.试算结果表明:相比于声波高斯束逆时偏移,本文方法能够对黏滞性引起的吸收衰减进行有效补偿,同时提取的角度域共成像点道集(ADCIGs)不仅可以用于分角度叠加成像压制成像噪声,而且能够为后续的偏移速度分析提供支撑.     

基于单程波方程的角度域保幅偏移(英文)

传统叠前深度偏移只能够提供地下的构造信息,但工业界在需要构造信息的同时还要与地下界面反射系数成比例的振幅信息。最近几年,基于单程波方程的保幅叠前深度偏移算法有了一定的发展,但是,基于炮域、单程波的保幅型叠前深度偏移必须应用反褶积型的成像条件,这种成像条件在构造复杂、速度变化剧烈的地区会出现不稳定现象。基于角度域的保幅深度偏移克服了这一不稳定性缺点的同时,还域的保幅深度偏移,模型和实际资料的试算分析验证该思路方法的正确性和有效性。     

一种时空域优化的高精度交错网格差分算子与正演模拟

如何有效压制数值频散是有限差分正演模拟研究中的关键问题之一.近年来,许多学者对二阶声波方程的差分算子开展了大量的优化工作,在压制频散方面取得不错的效果.一阶压强-速度方程广泛用于研究地震波在地下变密度模型中传播规律,目前针对一阶方程的优化工作大多只是在空间差分算子上展开.本文在前人研究的基础上,推导出一阶声波方程中压强场与偏振速度场之间的解析关系,据此在传统交错网格基础上给出一种高精度的显式时间递推格式,该递推格式将时间差分与空间差分算子结合在一起,并采用共轭梯度法得到精确时间递推匹配系数,实现时空差分算子的同时优化.在编程实现算法的基础上,通过频散分析与三个典型模型测试表明：本文方法能够较为有效地压制时间频散与空间频散,提高数值计算精度;同时对复杂模型也有很好适用性.     

弹性矢量波成像域联合层析反演

在地震弹性矢量波场框架下,推导了多波联合层析速度反演方程以及走时残差与角道集剩余曲率的转换关系式,提出了一种利用成像域角道集更新P波、S波速度的走时层析反演方法.其实现过程可以概括为:将弹性波多分量数据作为输入,基于高斯束实现矢量波场成像并提取角道集,利用层析反演方程求解慢度更新量,最终获得多波联合反演结果.模型试算和实际资料处理验证了该方法的反演效果,能够为弹性矢量波联合深度偏移提供高质量的叠前速度场.     

地震照明叠前深度偏移方法综述

对地震照明叠前深度偏移的基本概念、实现方法进行了分类和阐述.地震照明叠前深度偏移是通过对震源和炮记录进行合理的选择和合成,从而进行地震照明成像的一种有效方法.其可以分为平面波偏移和小束波偏移.理论模型的处理效果证明地震照明叠前深度偏移成像技术有很高的计算效率,并且还可以提高地下特定目标的成像质量.