摘要
根据旅行时数据用迭代方法可以确定地层速度变化的上界和下界.实际上,求得的是与数据一致的取值在一定范围内的模型,而不是在某种意义上“最佳拟合”的单一模型。本文所采用的算法是层析成像反演中常用的行运算法的一种改进形式,对于十万到百万量级的旅行时数据集,该算法只需要较小的内存,并已证明这是可行的.任何不等式约束,诸如速度总是正的,都可以加到公式中.数据误差分为局部误差(作为对每一个数据的严格约束)和整体误差(作为对总误差估算的一个约束),这种方法也可以用来导出使残差的l^(?)范数为最小的速度结构。 Grimsel 井间实验的数据可以用来将1521个旅行时残余值的0.02毫秒置信区间映射成地震速度的上下界.这些界限的宽度是地下界面位置的函数,具有0.1～0.9公里/秒较大的变化范围,界限的分布与奇异值分解(SVD)方法得到的参数精度相近,在旅行时资料层析成像反演时糜棱状低速带已充分成像。尽管对地震数据进行了各向异性校正,但是在井间弱约束区域出现的明显的正负交替速度扰动表明残留的各向异性影响仍然存在。
