基于Chebyshev走时逼近的三维多次反射射线计算
  • 【DOI】

    10.13278/j.cnki.jjuese.20170301

  • 【摘要】

    在利用有限差分等基于网格的数值分析方法解地震波走时所满足的程函方程时,由于速度模型的网格化离散等原因,会使走时在各网格节点之间不具有计算射线路径所要求的光滑性,即走时在邻近网格节点之间不具有连续的一阶导数.因此,直接利用网格节点走时计算射线路径会使最终的射线路径不光滑.为解决这个问题,已有研究者提出了基于B样条插值的逆向梯度方案(法).然而,在速度发生突变时,B样条逆向梯度法所计算出的射线路径会具... 展开>>在利用有限差分等基于网格的数值分析方法解地震波走时所满足的程函方程时,由于速度模型的网格化离散等原因,会使走时在各网格节点之间不具有计算射线路径所要求的光滑性,即走时在邻近网格节点之间不具有连续的一阶导数.因此,直接利用网格节点走时计算射线路径会使最终的射线路径不光滑.为解决这个问题,已有研究者提出了基于B样条插值的逆向梯度方案(法).然而,在速度发生突变时,B样条逆向梯度法所计算出的射线路径会具有较大的误差.针对这个问题,首先采用适合于解最小零偏差逼近及最佳平方逼近问题的Chebyshev多项式取代B样条对来自于分区多级计算方案的网格节点走时进行最佳逼近,得到在最小平方意义下的最优走时公式;然后采用与B样条逆向梯度法类似的计算过程得到光滑的射线路径.数值实验表明,利用Chebyshev多项式逼近走时可以得到具有很高精度的多次反射射线路径,在多次波偏移成像研究中具有潜在的价值. 收起<<

  • 【作者】

    孙建国  苗贺 

  • 【作者单位】

    吉林大学地球探测科学与技术学院

  • 【刊期】

    吉林大学学报(地球科学版) ISTIC EI PKU 2018年3期

  • 【关键词】

    多次反射  射线走时  射线路径  走时逼近  Chebyshev多项式  multiple reflections  ray traveltimes  ray trajectories  traveltimes approximation  Chebyshev polynomials 

  • 【基金项目】

    国家自然科学基金项目(41274120)Supported by National Natural Science Foundation of China