基于时空双变网格的起伏地表变坐标系正演模拟方法（英文）

坐标变换法通过将物理空间的曲网格映射为计算空间的矩形网格,将起伏地表转化为水平地表,同时将物理空间的波动方程转化为计算空间的波动方程,在计算空间完成数值模拟,坐标变换的方法对处理起伏自由边界具有较好的适应性和应用效果。本文在传统坐标变换方法的基础上,根据计算区域速度差异采用不同的网格大小和采样时间步长,提出了一种基于时空双变网格的起伏地表坐标变换正演模拟方法。在编程实现算法的基础上,通过典型模型波场模拟试算结果分析可知:(1)变网格方法与常规方法波场模拟误差在0.5%左右;(2)变网格方法计算效率视不同的变网格区域面积及变网格大小可提高几倍量级,在本文模型和计算参数下提高约5倍。(3)在满足模拟精度及频散条件要求下,变网格方法较全局细网格算法能显著节约计算内存。为此,针对起伏地表数值模拟,本文方法具有较高的模拟计算精度和一定的适应性。     

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

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

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

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

Lebedev网格改进差分系数TTI介质正演模拟方法研究

本文采用一种新的交错网格-Lebedev网格(LG)进行TTI介质的正演模拟研究,避免了Virieux标准交错网格(SSG)算法在处理TTI、单斜等各向异性介质时波场插值引入的数值误差,提高了模拟精度.在方法实现过程中,本文针对有限差分正演模拟面临的网格频散与边界反射两个关键性问题分别做了优化,并通过模型试算验证了它们的有效性与可行性:(1)结合最小二乘思想推导出新的频散改进差分系数(DIC),该系数比Taylor系数更能有效地压制粗网格引起的数值频散,可以节约内存,提高计算效率;(2)将分裂的多轴完全匹配层(M-PML)吸收边界条件引入到LG算法中,解决了传统PML边界条件在某些各向异性介质中的不稳定现象并且具有较好的边界吸收效果.     

裂缝介质横波分裂时差影响因素定量研究

裂缝介质数值模拟是研究裂缝型油藏各向异性及裂缝检测的关键.横波分裂现象是波在各向异性介质中传播的一个显著特征,横波分裂时差能够反映地下裂隙密度和介质物性等信息.本文给出二维均匀TTI介质在不同极化角度下的横波群速度公式,进而得到偏移距表征的横波分裂时差表达式;同时采用二维三分量交错网格高阶有限差分方法模拟不同极化角度下的三分量炮记录,与理论值对比分析,定量研究分裂时差与极化角、偏移距、各向异性参数的关系,得到如下主要结论:1)在不同各向异性强度下,横波分裂时差与极化角呈现不同的变化规律,而快慢横波振幅比与极化角的关系不受各向异性强度的影响;2)从能量角度,在45°～55°极化角内,横波传播的能量对比最强,是观测横波分裂的最佳角度范围;3)表征各向异性强度的三个参数中γ对横波分裂时差的影响最大,其从0变化到0.3,分裂时差增幅达432.4%.     

基于一阶速度-应力方程的多震源最小二乘逆时偏移

最小二乘逆时偏移(Least-Square Reverse Time Migration,LSRTM)相比于常规偏移具有更高的成像分辨率、振幅保幅性及均衡性等优势,是当前研究的热点之一.然而,目前LSRTM算法大多是基于二阶常密度标量声波方程建立的,忽略了密度变化对振幅的影响,因而基于振幅匹配策略的常规LSRTM很难在变密度介质下取得保真的成像结果.一阶速度-应力方程能够很好地处理变密度介质,但简单地将一阶速度-应力方程应用到LSRTM中缺乏理论基础.为此,本文从LSRTM的正问题入手,提出了基于交错网格的一阶速度-应力方程LSRTM理论方法.首先将一阶波动方程线性化,建立了一阶方程LSRTM的目标泛函,随后推导其伴随方程,并借助伴随状态法给出了迭代更新流程,最终建立了基于一阶速度-应力方程LSRTM的理论框架.进一步,通过在相位编码LSRTM中引入随机最优化思想,极大地减小了计算量、提高了计算效率.最后,通过模型试算验证了本算法的正确性和有效性.     

基于贴体全交错网格的起伏地表正演模拟影响因素

贴体网格有限差分正演模拟算法不仅能够精确模拟任意起伏地形下的波场特征,且计算效率较高,是一种很有应用前景的处理西部复杂地表问题的方法;然而,目前求解波动方程时常用的同位网格和标准交错网格,在处理贴体网格起伏地表正演模拟时存在诸多问题。为此,将全交错网格引入到曲线坐标系下,避免了标准交错网格的插值误差和同位网格中奇偶失联引起的高频振荡现象,提高了模拟精度,减小了算法实现的复杂度。在自由边界条件实施时,采用牵引力镜像法计算速度分量,速度自由边界条件配合紧致交错差分格式更新应力分量,得到了较好的效果。随后,重点研究了贴体全交错网格正演模拟算法的影响因素,考虑了网格正交性、网格间距和网格拼接等的影响,并取得了如下认识：算法对网格的正交性没有过分要求;网格间距的突变会引起虚假反射的产生;不同类型的网格拼接对模拟结果不会造成明显的影响。     

一种基于平面波静态编码的最小二乘逆时偏移方法

平面波偏移是一种面炮偏移方法,相对于常规逐炮偏移,其具有较高的计算效率.然而常规平面波偏移方法成像精度低,且成像时会产生串扰噪音.为此,本文在实现常规平面波偏移算法基础上,引入反演思想实现了基于静态平面波编码的最小二乘偏移理论方法及处理流程,在优化算法基础上对平层模型和复杂砂砾断块模型进行了成像测试并与其他成像策略进行对比.研究结果表明:基于时移编码的平面波最小二乘偏移能有效抑制低频成像噪音和串扰噪音,补偿中深部成像能量,是一种较为有效的保幅成像策略.     

Lebedev网格改进差分系数TTI介质正演模拟方法研究

本文采用一种新的交错网格-Lebedev网格（LG）进行TTI介质的正演模拟研究,避免了Virieux标准交错网格（SSG）算法在处理TTI、单斜等各向异性介质时波场插值引入的数值误差,提高了模拟精度.在方法实现过程中,本文针对有限差分正演模拟面临的网格频散与边界反射两个关键性问题分别做了优化,并通过模型试算验证了它们的有效性与可行性：（1）结合最小二乘思想推导出新的频散改进差分系数（DIC）,该系数比Taylor系数更能有效地压制粗网格引起的数值频散,可以节约内存,提高计算效率;（2）将分裂的多轴完全匹配层（M-PML）吸收边界条件引入到LG算法中,解决了传统PML边界条件在某些各向异性介质中的不稳定现象并且具有较好的边界吸收效果.     

基于Student′s t分布的不依赖子波最小二乘逆时偏移

最小二乘逆时偏移(Least-Squares Reverse Time Migration,LSRTM)与常规偏移相比具有更高的成像分辨率、振幅保真性及均衡性等优势,是当前研究的热点之一.震源子波的估计直接影响LSRTM结果的好坏,在实际情况下考虑到震源子波的空变特性,其估计十分困难.为了消除子波对LSRTM结果的影响,本文发展了基于卷积目标泛函的不依赖子波LSRTM算法.目标泛函由观测记录卷积模拟记录的参考道以及模拟记录卷积观测记录的参考道组成,由于观测子波和模拟子波在目标泛函的两项中同时存在,从而消除了子波的影响.此外,常用的基于L2范数拟合的LSRTM算法对噪声非常敏感,尤其是当地震数据中含有异常值时,常规LSRTM无法得到满意的结果.Student′s t分布相比L2范数具有更好的稳健性,本文将其推广到不依赖子波LSRTM中,提升了算法的稳健性,最后通过理论模型及实际资料试算验证了算法的有效性和对复杂模型的适应性.