6:15 PM - 7:30 PM
[SSS26-P01] Waveform inversion converged towards the grand minimum: A Wiener-filter approach
Keywords:inverse problem, seismic waveform, Wiener filter, grand minimum, phase unwrapping
We propose a new inversion approach in order to make waveform inversions objective and stable. This consists of two stages. At first, the difference between each observed waveform d(t) and synthetic one p(t) (here we shall assume that all the waveforms are expressed by discrete time series) by a Wiener filter w(t). The Wiener filter is defined by the filter that minimize the sum of squares of p(t)*w(t)-d(t), where * represents the convolution. Coefficient of w(t) can be obtained efficiently by a recursive algorithm (Levinson, 1949). If p(t) and d(t) were same, the Wiener filter w(t) would be a unit filter, that is, 1 for t=0 and zero otherwise (corresponding to a delta function in a continuous case). In the second step, therefore, we introduce a new criterion to minimize the sum of squares of w(t) multiplied by a time lag for all the records together: t•w(t). In practice, we need to normalize the above criterion divided by the square of w(t).
With this new criterion, the time lag of w(t) is only weighted without its sign (i.e., positive peak versus negative trough), there are no ambiguities in phase, unlike any conventional waveform inversions. If we iterate the above inversion procedure by modifying a given model in each step, our model will be converged towards a optimal one without falling into any local minimums.
For large-scale inversion problems, one uses an adjoint matrix or operator rather than the formal inverse matrix of a given forward problem. The modification of a model in each iteration step is obtained by the data residual, p(t)–d(t), multiplied by the adjoint matrix. The adjoint matrix is derived by the partial derivative of by the minimizing criterion (the sum of squares of t•w(t) in this case) partially differentiated by the synthetic p(t). Its result is expressed by - t2•w(t)2* p(t)-1. That is , we only need to deconvolve the sum of squares of t•w(t) by p(t). Other procedures are similar to those of conventional waveform inversions.
We conducted numerical tests in order to confirm the efficiency of our new inversion approach, using a medium with velocity fluctuations. Even without a dense coverage of paths connecting sources and receivers, our inversion leads to a final model close to the true model. In contrast, an initial model very close to the true one is required if we obtain a satisfactory model with conventional least-squares criteria in inversion. That is, the present inversion approach can lead our model search into the true or grand minimum without falling into or being stuck by local minimums around it.