We develop a new phase space method for reconstructing the index of refraction of a medium from travel time measurements. The method is based on the so-called StefanovUhlmann identity which links two Riemannian metrics with their travel time information. We design a numerical algorithm to solve the resulting inverse problem. The new algorithm is a hybrid approach that combines both Lagrangian and Eulerian formulations. In particular the Lagrangian formulation in phase space can take into account multiple arrival times naturally, while the Eulerian formulation for the index of refraction allows us to compute the solution in physical space. Numerical examples including isotropic metrics and the Marmousi synthetic model are shown to validate the new method.

No abstract uploaded!

No abstract uploaded!

No abstract uploaded!

Let$L$be a formally selfadjoint differential operator and$p$a real-valued function, both on$a$≤$x$<∞. The deficiency indices are the numbers of solutions of$Lu$=λ$pu$for Im λ>0 and for Im λ<0 which have a certain regularity at$x$=∞. (A) If$p(x)$≥0 this regularity means that the integral of$p(x)$│$u$│^{2}converges at infinity. (B) If$p$changes its sign for arbitrarily large values of$x$but$L$has a positive definite Dirichlet integral it is natural to relate the regularity to this integral. Weyl’s classical study of the deficiency indices is reviewed for (A) with the help of elementary theory of quadratic forms. Individual bounds are found for the deficiency indices also when$L$is of odd order. It is then indicated how the method carriers over to (B).