A procedure for deconvolution of nonminimum phase non-Guassian time series based on the estimation of higher order (greater than two) spectra is given. This can be applied to the analysis of seismograms. The procedure allows estimation of the wavelet. Knowledge of cumulant spectra of order greater than two allows estimation of the phase of the wavelet. In this way one has access to information not available in the ordinary second-order deconvolution procedures. Computational details of the method for estimating the phase of the wavelet are given. There are simulated illustrative examples. One of the examples is based on an actual reflectivity series from a sonic well log. The method is effective asymptotically in the nonminimum phase non-Gaussian context where the Wiener-Levinson procedure does not apply.
