English   español  
Please use this identifier to cite or link to this item: http://hdl.handle.net/10261/124227
logo share SHARE logo core CORE   Add this article to your Mendeley library MendeleyBASE

Visualizar otros formatos: MARC | Dublin Core | RDF | ORE | MODS | METS | DIDL
Exportar a otros formatos:


TOMO3D: 3-D joint refraction and reflection traveltime tomography parallel code for active-source seismic data-synthetic test

AuthorsMeléndez, Adrià ; Korenaga, Jun; Sallarès, Valentí ; Miniussi, A.; Ranero, César R.
KeywordsInverse theory
Controlled source seismology
Seismic tomography
Issue DateOct-2015
PublisherOxford University Press
CitationGeophysical Journal International 203(1): 158-174 (2015)
AbstractPublished by Oxford University Press on behalf of The Royal Astronomical Society. We present a new 3-D traveltime tomography code (TOMO3D) for the modelling of active-source seismic data that uses the arrival times of both refracted and reflected seismic phases to derive the velocity distribution and the geometry of reflecting boundaries in the subsurface. This code is based on its popular 2-D version TOMO2D from which it inherited the methods to solve the forward and inverse problems. The traveltime calculations are done using a hybrid ray-tracing technique combining the graph and bending methods. The LSQR algorithm is used to perform the iterative regularized inversion to improve the initial velocity and depth models. In order to cope with an increased computational demand due to the incorporation of the third dimension, the forward problem solver, which takes most of the run time (~90 per cent in the test presented here), has been parallelized with a combination of multi-processing and message passing interface standards. This parallelization distributes the ray-tracing and traveltime calculations among available computational resources. The code's performance is illustrated with a realistic synthetic example, including a checkerboard anomaly and two reflectors, which simulates the geometry of a subduction zone. The code is designed to invert for a single reflector at a time. A data-driven layer-stripping strategy is proposed for cases involving multiple reflectors, and it is tested for the successive inversion of the two reflectors. Layers are bound by consecutive reflectors, and an initial velocity model for each inversion step incorporates the results from previous steps. This strategy poses simpler inversion problems at each step, allowing the recovery of strong velocity discontinuities that would otherwise be smoothened. © The Authors 2015
Description17 pages, 11 figures, 1 table, supporting Information http://gji.oxfordjournals.org/content/203/1/158/suppl/DC1
Publisher version (URL)http://dx.doi.org/10.1093/gji/ggv292
Identifiersdoi: 10.1093/gji/ggv292
issn: 0956-540X
e-issn: 1365-246X
Appears in Collections:(ICM) Artículos
Files in This Item:
File Description SizeFormat 
Melendez_et_al_2015.pdf14,76 MBAdobe PDFThumbnail
Melendez_et_al_2015_supp_material.pdf3,4 MBAdobe PDFThumbnail
Show full item record
Review this work

Related articles:

WARNING: Items in Digital.CSIC are protected by copyright, with all rights reserved, unless otherwise indicated.