2　全波形反演

本文建立的全波形反演思路与Singh等［4］相近，是一种水平层状弹性介质的截距时间（τ）-水平慢度（p）域反演方法，分两步。第一步是根据τ-p域走时数据应用非常快速模拟退火算法求得速度结构的长波长部分。第二步，基于第一步的速度模型作为初始模型，利用波形资料用共轭梯度法求得速度的短波长扰动部分。这样，最后反演得到的模型包含了速度结构的长波长与短波长部分。反演结合了全局搜索方法与非线性局部搜索方法。(https://www.daowen.com)

这一反演方法包含了走时反演与波形反演，全局反演与局部优化反演。走时反演用于求取长波长速度信息，波形反演求取短波长速度信息。前者非线性程度较高，应用了全局搜索方法，而后者非线性程度较低，应用了局部搜索方法。这也是当前实现地震波形反演的现实途径。全局反演相当耗时，大量模型参数的全局反演不太现实，因此采取了折衷的方法。正演是用Kennett［14］的反射透射矩阵方法实现的。反演的数据是共中心点道集经过平面波分解得到τ-p域地震数据。并需考虑震源子波估算、地震数据标定与组合校正等各种预处理。反演中应用的快速模拟退火算法与共轭梯度算法如下所述。