Seismic data decomposition into spectral components using regularized nonstationary autoregression |

Prony's method of data analysis was developed originally for
representing a noiseless signal as a sum of exponential components
(Prony, 1795). It was extended later to noisy signals, complex
exponentials, and spectral analysis
(Pisarenko, 1973; Beylkin and Monzón, 2005; Marple, 1987; Bath, 1995). The basic idea follows
from the fundamental property of exponential functions:
. In
signal-processing terms, it implies that a time sequence
(with real or complex ) is predictable
by a two-point prediction-error filter
, or, in the Z-transform notation,

they can be predicted simultaneously by using a convolution of several two-point prediction-error filters:

where . This observation suggests the following three-step algorithm:

- Estimate a prediction-error filter from the data by determining filter coefficients
from the least-squares minimization of

- Writing the filter as a polynomial (equation 8), find its complex roots
. The exponential factors
are determined then as

- Estimate amplitudes of different components in equation 7 by linear least-squares fitting.

Prony's method can be applied in sliding windows, which was a technique developed by Russian geophysicists (Mitrofanov and Priimenko, 2011; Gritsenko et al., 2001) for identifying low-frequency seismic anomalies (Mitrofanov et al., 1998). I propose to extend it to smoothly nonstationary analysis by applying the following modifications:

- Using RNAR, the filter coefficients become smoothly-varying functions of time , which allows the filter to adapt to nonstationary changes in the input data.
- At each instance of time, roots of the corresponding polynomial also become functions of time . I apply a robust, eigenvalue-based algorithm for root finding (Toh and Trefethen, 1994).
The instantaneous frequency of different components is
determined directly from the phase of different roots:

- Finally, using RNR, I estimate smoothly-varying amplitudes of different components .

and the local phase corresponds to time integration of the instantaneous frequency determined in Step 2:

For ease of analysis, real signals can be transformed to the complex domain by using analytical traces (Taner et al., 1979).

Seismic data decomposition into spectral components using regularized nonstationary autoregression |

2013-10-09