3-D simultaneous inversion for velocity and reflector

1 downloads 0 Views 1018KB Size Report
⋅V xe k. А Б. (1) where xk e and V(xk e) (k=1,2,......,8) are the vector coordi- nates and the ..... but 6 was the best) and the damping factor is 0.01 (soft one). The.
3-D simultaneous inversion for velocity and reflector geometry using multiple classes of arrivals in a spherical coordinate frame Guo-jiao Huang, Chao-ying Bai, Xingwang Li & Stewart Greenhalgh

Journal of Seismology ISSN 1383-4649 Volume 18 Number 1 J Seismol (2014) 18:123-135 DOI 10.1007/s10950-013-9406-z

1 23

Your article is protected by copyright and all rights are held exclusively by Springer Science +Business Media Dordrecht. This e-offprint is for personal use only and shall not be selfarchived in electronic repositories. If you wish to self-archive your article, please use the accepted manuscript version for posting on your own website. You may further deposit the accepted manuscript version in any repository, provided it is only made publicly available 12 months after official publication or later and provided acknowledgement is given to the original source of publication and a link is inserted to the published article on Springer's website. The link must be accompanied by the following text: "The final publication is available at link.springer.com”.

1 23

Author's personal copy J Seismol (2014) 18:123–135 DOI 10.1007/s10950-013-9406-z

ORIGINAL ARTICLE

3-D simultaneous inversion for velocity and reflector geometry using multiple classes of arrivals in a spherical coordinate frame Guo-jiao Huang & Chao-ying Bai & Xing-wang Li & Stewart Greenhalgh

Received: 28 December 2012 / Accepted: 30 October 2013 / Published online: 3 December 2013 # Springer Science+Business Media Dordrecht 2013

Abstract The two key requirements in conducting 3-D simultaneous traveltime tomography on real data at the regional and global scale with multiple classes of arrival time information are (1) it needs an efficient and accurate arrival tracking algorithm for multiply transmitted, reflected (or refracted) and converted waves in a 3-D variable velocity model with embedded velocity discontinuities (or subsurface interfaces), and (2) a subdimensional inversion solver is required which can easily search for different types of model parameters to balance the tradeoff between the different types of model parameter updated in the simultaneous inversion process. For these purposes, we first extend a popular grid/cell-based wavefront expanding ray tracing algorithm (the multistage irregular shortest-path ray tracing method), which previously worked only in Cartesian coordinate at the local scale, to spherical coordinates appropriate to the regional or global scale. We then incorporated a fashionable inversion solver (the subspace method) to formulate a simultaneous inversion algorithm, in which the multiple classes of G.

Suggest Documents