We study the uniqueness and accuracy of the numerical solution of the problem of reconstruction of the shape and trajectory of a reflecting obstacle moving in an inhomogeneous medium from travel times, start and end points, and initial angles of ultrasonic rays reflecting at the obstacle. The speed of sound in the domain when there is no obstacle present is known and provided as an input parameter which together with the other initial data enables the algorithm to trace ray paths and find their reflection points. The reflection points determine with high-resolution the shape and trajectory of the obstacle. The method has predictable computational complexity and performance and is very efficient when it is parallelized and optimized because only a small portion of the domain is reconstructed.
[1]
G. Eskin.
Lectures on Linear Partial Differential Equations
,
2011
.
[2]
Vasili M. Babič,et al.
Short-Wavelength Diffraction Theory
,
1991
.
[3]
V. Červený,et al.
Seismic Ray Theory
,
2001,
Encyclopedia of Solid Earth Geophysics.
[4]
O. Basset,et al.
Second harmonic inversion for ultrasound contrast harmonic imaging
,
2011,
Physics in medicine and biology.
[5]
Jack K. Cohen,et al.
Mathematics of Multidimensional Seismic Imaging, Migration, and Inversion
,
2001
.
[6]
Malcolm Sambridge,et al.
Boundary value ray tracing in a heterogeneous medium: a simple and versatile algorithm
,
1990
.
[7]
David Gubbins,et al.
Three-dimensional seismic ray tracing
,
1977
.