Large Dip Calculation via Robust Nonstationary Plane-Wave Destruction

The omnidirectional plane-wave destruction (OPWD) algorithm can estimate the large dip by using circle-interpolating plane-wave destruction (PWD) filter but at the expense of causing potential instabilities due to the small values of the denominator in the regularized division problem. To mitigate the instability, one needs to use a relatively larger smoothing radius for a stronger regularization of the element-wise division, which however significantly decreases the resolution of dip estimation. We propose a new OPWD method without compromising the dip resolution for the large dip calculation by applying a nonstationary smoothness constraint to the model. We use a larger smoothing radius for areas that tend to cause instabilities and vice versa. The nonstationary smoothing is carried out in a simple and efficient recursion way. The synthetic and real seismic data examples demonstrate the performance of the proposed algorithm.