Design and implementation of a 2D circular domain seismic tomography general ray tracer algorithm

• Published in: Ingenieria y Desarrollo Vol. 32; no. 2; pp. 242 - 259
• Main Authors: Imhof, Armando Luis, Calvo, Carlos Adolfo
• Format: Journal Article
• Language: English; Spanish
• Published: 01.07.2014
• Subjects: Algorithms; Arrays; Circularity; Construction; Discretization; Inversions; Pixels; Quadrants
• ISSN: 0122-3461; 2145-9371; 2145-9371
• DOI: 10.14482/inde.32.2.5717

Tomographic matrix inversion techniques in travel time are based on the construction of distance matrices containing information of survey array geometry, the discretized domain shape and the distance traveled by rays in each area elements that compose it. The domain considered is generally rectangular for ease of geometric resolution. In this work is circular due to the need to carry out investigation in columns of this shape. Distance matrices are used to solve the linear system, which leads to the determination of velocity values in each element of the domain. Taking in mind discretization, there are several possibilities for the construction of arrays of pixels: may those be defined in one quadrant and obtaining the rest by symmetry; build arrays with uniform pixels using the previous strategy, and finally to be independent from quadrant and generate any distribution and individual size thereof. Each option will provide access to new research possibilities. In this work an original procedure was developed for the construction of matrices with elements of the last type. The procedure was validated with simulated data using the SVD method, which demonstrated the efficiency of the algorithm.