Dinghui Yang

Jing Chen, Yifan Chen, Hao Wu et al.

In [Engquist et al., Commun. Math. Sci., 14(2016)], the Wasserstein metric was successfully introduced to the full waveform inversion. We apply this method to the earthquake location problem. For this problem, the seismic stations are far from each other. Thus, the trace by trace comparison [Yang et al., arXiv(2016)] is a natural way to compare the earthquake signals. Under this framework, we have derived a concise analytic expression of the Frèchet gradient of the Wasserstein metric, which leads to a simple and efficient implementation for the adjoint method. We square and normalize the earthquake signals for comparison so that the convexity of the misfit function with respect to earthquake hypocenter and origin time can be observed numerically. To reduce the impact of noise, which can not offset each other after squaring the signals, a new control parameter is introduced. Finally, the LMF (Levenberg-Marquardt-Fletcher) method is applied to solve the resulted optimization problem. According to the numerical experiments, only a few iterations are required to converge to the real earthquake hypocenter and origin time. Even for data with noise, we can obtain reasonable and convergent numerical results.

Hao Wu, Jing Chen, Xueyuan Huang et al.

In this paper, a new earthquake location method based on the waveform inversion is proposed. As is known to all, the waveform misfit function is very sensitive to the phase shift between the synthetic waveform signal and the real waveform signal. Thus, the convergence domain of the conventional waveform based earthquake location methods is very small. In present study, by introducing and solving a simple sub-optimization problem, we greatly expand the convergence domain of the waveform based earthquake location method. According to a large number of numerical experiments, the new method expands the range of convergence by several tens of times. This allows us to locate the earthquake accurately even from some relatively bad initial values.

Hao Wu, Jing Chen, Hao Jing et al.

This paper introduces the auxiliary function method (AFM), a novel, fast and simple approach for waveform based earthquake location. From any initial hypocenter and origin time, we can construct the auxiliary function, whose zero set contains the real earthquake hypocenter and the origin time. In most of situations, there are very few elements in this set. The overall computational cost of the AFM is significantly less than that of the iterative methods. According to our numerical tests, even for large noise, the method can still achieve good location results. These allow us to determine the earthquake hypocenter and the origin time extremely fast and accurate.