Calculation of differential seismograms using analytic partial derivatives—II Regional waveforms. Issue 1 (10th October 2017)

• Record Type: Journal Article
• Authors: Hu, Shaoqian; Zhu, Lupei
• Is Part Of: Geophysical journal international. Volume 212:Issue 1(2018:Jan.)
• Page Start: 390
• Page End: 399
• Publication Date: 2017-10-10
• Subjects: Waveform inversion -- Theoretical seismology -- Wave propagation; Geophysics -- Periodicals; 550
• DOI: 10.1093/gji/ggx433
• Languages: English
• ISSNs: 0956-540X

Abstract

SUMMARY: Solving seismic waveform inverse problem by the linearization approach requires calculating partial derivatives of seismograms (differential seismograms) with respect to velocity model parameters. Here we extend our previous work on calculation of differential teleseismic receiver functions using analytic partial derivatives to regional waveforms. We derived analytic partial derivatives of displacement kernels for both P – SV and SH waves generated by a point source in a multilayered half-space. They are expressed explicitly in terms of the analytic partial derivatives of the Haskell propagator matrix and the source vector. Differential seismograms are then computed using the wavenumber–frequency double integration method. We implemented an efficient algorithm to speed up the computation by storing intermediate matrix products. The total computation time for computing differential seismograms with respect to all the layers is only about one extra forward-problem computation time. Numerical tests show that our derivations are correct and the implementation is efficient. Synthetic and real data waveform inversions using more accurate analytic differential seismograms perform better than using finite-difference differential seismograms.