English   español  
Please use this identifier to cite or link to this item: http://hdl.handle.net/10261/140937
logo share SHARE   Add this article to your Mendeley library MendeleyBASE
Visualizar otros formatos: MARC | Dublin Core | RDF | ORE | MODS | METS | DIDL
Exportar a otros formatos:

Data-driven layer-stripping strategy in 3-D joint refraction and reflection travel-time tomography with TOMO3D

AuthorsMeléndez, Adrià ; Korenaga, Jun; Sallarès, Valentí ; Miniussi, A.; Ranero, César R.
Issue Date14-Apr-2015
PublisherEuropean Geosciences Union
CitationGeophysical Research Abstracts 17: EGU2015-7666 (2015)
AbstractWe present a new 3-D travel-time 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 propagation velocity distribution and the geometry of reflecting boundaries in the subsurface. The combination of refracted and reflected data provides a denser coverage of the study area. Moreover, because refractions only depend on the velocity parameters, they contribute to the mitigation of the negative effect of the ambiguity between layer thickness and propagation velocity that is intrinsic to the reflections that define these boundaries. This code is based on its renowned 2-D version TOMO2D from which it inherited the methods to solve the forward and inverse problems. The forward travel-time calculations are conducted using a hybrid ray-tracing technique combining the graph or shortest path method and the bending method. The LSQR algorithm is used to perform the iterative inversion of travel-time residuals to update the initial velocity and depth models. In order to cope with the increased computational demand due to the incorporation of the third dimension, the forward problem solver, which takes by far most of the run time (~90%), has been parallelised with a combination of MP and MPI standards. This parallelisation distributes the ray-tracing and travel-time calculations among the available computational resources, allowing the user to set the number of nodes, processors and cores to be used. The code¿s performance was evaluated with a complex synthetic case simulating a subduction zone. The objective is to retrieve the velocity distribution of both upper and lower plates and the geometry of the interplate and Moho boundaries. Our tomography method is designed to deal with a single reflector per inversion, and we show that a data-driven layer-stripping strategy allows to successfully recover several reflectors in successive inversions. This strategy consists in building the final velocity model layer by layer, sequentially extending it down with each inversion of a new, deeper reflector. One advantage of layer stripping is that it allows us to introduce and keep strong velocity contrasts associated to geological discontinuities that would otherwise be smoothened. Another advantage is that it poses simpler inverse problems at each step, facilitating the minimisation of travel-time residuals and ensuring a good control on each partial model before adding new data corresponding to deeper layers. Finally, we discuss the parallel performance of the code in this particular synthetic case
DescriptionEuropean Geosciences Union General Assembly 2015 (EGU2015), 12-17 April 2015, Vienna, Austria.-- 1 page
Publisher version (URL)http://meetingorganizer.copernicus.org/EGU2015/orals/17375
Identifiersissn: 1607-7962
Appears in Collections:(ICM) Comunicaciones congresos
Files in This Item:
File Description SizeFormat 
Melendez_et_al_2015.pdf45,04 kBAdobe PDFThumbnail
Show full item record
Review this work

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