基于分频编码的弹性波全波形反演

计算量与计算效率限制了全波形反演方法的应用，震源编码技术可以有效减少全波形反演的计算量，提升反演效率。传统震源编码技术要求各炮具有相同的接收点排列，在观测系统适应性方面存在不足，无法直接应用于滚动排列观测系统。为此，提出一种基于分频编码的弹性波全波形反演方法，即在正演过程中仍可同时对多个炮记录进行正向延拓，给每一炮赋以不同频率的谐波震源，使同时进行正向延拓的各炮波场在频谱上互不重合；在炮点波场正传和残差反传过程中，采用相位灵敏度检测技术提取每炮的单频波场，实现各炮波场的完全分离，用于构建波形反演梯度，实现模型更新。同传统震源编码技术相比，该方法不受观测系统限制，能应用于排列滚动的观测系统，具有更广的应用范围。     

一阶速度—胀缩—旋转弹性波方程交错网格数值模拟

弹性波正演在地震波传播机理研究以及多波地震资料采集、处理、解释和反演中发挥着重要作用。现有的弹性波方程正演模拟常常数值求解一阶速度—应力方程或二阶位移弹性波方程,只能直接得到同时包含纵波和横波的三个质点振动速度分量或位移分量,要想得到更直观的纯纵波和纯横波分量记录,还需要在模拟过程中采用波场解耦算子进行纵、横波分离,因此纵、横波模拟精度同时受制于模拟算法和波场解耦算法的精度。为此,推导了一阶速度—胀缩—旋转弹性波方程在三维交错网格空间中的高阶有限差分格式,并给出了相应的稳定性条件;推导了适应该方程的PML吸收边界条件,实现了一阶速度—胀缩—旋转弹性波方程的正演模拟;分析了模拟结果中各分量的物理意义。由于一阶速度—胀缩—旋转弹性波方程不仅包含了质点的振动速度矢量,而且显式地包含了横波振动速度矢量和纵波振动速度矢量,还包含了一个体应变和一个旋转矢量,因此应用该方程模拟除了能得到三个质点振动速度分量外,还可以直接得到解耦后的纵、横波分量,避免了解耦算法对模拟精度的影响。模型试算证明了该模拟方法的正确性和优越性。