韩复兴， 孙章庆， 陈祥忠， 高正辉， 王瑞兴， 刘俊成， 王丽丽

1 吉林大学地球探测科学与技术学院， 长春 1300262 北京桔灯地球物理勘探股份有限公司， 北京 1022003 吉林油田地球物理勘探研究院， 松原 138006

0 引言

基于弹性波传播方程高频近似的射线理论类偏移成像技术(Hill, 1990, 2001; 高成等, 2015; Gao et al., 2017；韩建光等，2017)，其在实现偏移成像的过程当中，采用射线理论方法求解运动学和动力学射线追踪系统获得波场传播所需的走时振幅和相位信息，然后进行地震波场的延拓外推从而达到最终成像的目的.其实现过程和实现方式相对于波动方程类偏移成像具有很大的时效性和经济性，因此一直是工业界应用的主流成像方法和技术(岳玉波, 2011；杨珊珊等，2015；袁茂林等，2016).基于体波、面波等地震弹性波和电磁波的有效联合，解译了华南深部结构(Di et al., 2021).那么在应用射线理论实现波场延拓成像的过程当中，难免会遇到依赖速度模型进行的射线追踪系统的求解.在实际应用过程当中，一般给定的速度模型都比较复杂，而且是以离散的形式只给定网格节点上的数值，非网格节点上的值只能通过构建的插值算法来进行插值计算.在应用射线理论进行运动学和动力学追踪计算中，波场随着时间的不断外推增加，每一步计算所应用到的速度场信息不一定都落在给定的原离散网格上，如果计算所需要的速度场不在给定的离散节点之上，就需要对速度模型进行插值计算，从而解决在波前外……