基于改进的散射积分算法的初至波走时层析.pdfVIP

第 59卷 第 1O期 地 球 物 理 学 报 Vo1．59。NO．10 2016年 1O月 CHINE『SE JOURNAL OF GEOPHYSICS Oct．，2O16 李勇德 ，董 良国，刘玉柱．2016．基于改进 的散射积分算法 的初至波走时层析．地球物理学报 ，59(10)：3820～3828，doi：10 6038／cj LiY D。DongL G．LiuY Z．2016．Firstarrivaltraveltimetomographybasedonanimprovedscattering—integralalgorithm． ChineseJ．Geophys．(inChinese)，59(10)：3820—3828，doi：10．6O38／cjg2O161026． 基于改进的散射积分算法的初至波走时层析 李勇德，董 良国，刘玉柱 同济大学海洋地质 国家重点实验室，上海 200092 摘要 初至波走时层析是获取近地表速度结构的一种常用方法．随着采集技术 的不断发展 ，可使用的数据量迅速 增多，传统的基于射线追踪和解方程组的地震走时层析成像方法面临着内存 占用大、方程求解不稳定等 问题 ．为 了 解决这些问题 ，本文基于前人在波形反演研究中提出的一种改进 的散射积分算法 ，提 出了一种预条件最速下降法 初至波走时层析．该方法无需存储核函数矩阵与Hessian矩阵即可方便地实现 目标函数梯度的计算与预条件，且该 方法计算效率高、求解稳定、易于并行．数值实验结果表明，该方法可以获得与传统方法精度相当的反演结果 ，但所 占用的内存大幅减小． 关键词 地震层析成像 ；散射积分法 ；初至波走时；近地表速度 ；海森矩阵；预条件 ；最速下降法 doi：10．6O38／cjg2O161O26 中图分类号 P631 收稿 日期 2015—1205，2016—08—26收修定稿 Firstarrivaltraveltimetomographybased onan improvedscattering-integralalgorithm LIYong—De，DONG Liang—Guo，LIU Yu—Zhu StateKeyLaboratoryofMarineGeology，TongjiUniversity，Shanghai200092．China Abstract Firstarrivaltravehimetomographyisacommonlyusedmethodatpresenttoobtainthe near—surfacevelocitystructure．W iththedevelopmentofacquisitiontechnology，thedataquantity isgrowingrapidly．Thus。traditionalray—based equations—solvingmethodfacesthechallengesof memoryconsuming and instability．To solve theseproblems，an improved scattering—integral (SI)algorithm proposedinaformerstudyonfull—waveform inversion (FW I)isemployedhereto developanewpreconditionedsteepest—descentray—basedfirstarrivaltraveltimetomography ． This methodcan computethegradientoftheobjectivefu