论文部分内容阅读
地震射线层折成像方法是一个非线性的反演过程。从方法原理上讲，它主要由两部分组成，即任意速度分布的射线追踪和求解线性化反演问题。本文中，作者在这两个方面均采用了新的方法。在正演中，（1）应用Moser[1]提出的“最短路径”原理模拟地震射线，用较为简单的方法实现任意复杂介质结构的射线追踪，计算出每一对地点一接收点的最小走时及其对应的射线路径；（2）用线性平面函数去近似矩形单元内的速度分布，从而求出走时相对于每个节点上的速度变化的雅可比矩阵元素的解析式。在反演中，给出两种求解加权最小范数、最小二乘和对偶空间约束反演问题[2]的迭代算法。通过我们的计算机模拟实验表明，这种方法比常规的地震层析成像方法（如ART，LSQR，CGLS）不仅在精度上大大提高，而且有较强的抗噪声能力。在计算方面采用了Scales[3]提出的大型稀疏矩阵的算法，既加快了运算速度，又不需要很大的内存空间。