Published as Geophysical Journal International, 214, 2207-2223, (2018)

Plane-wave orthogonal polynomial transform for amplitude-preserving noise attenuation

Yangkang Chen % latex2html id marker 4302
\setcounter{footnote}{1}\fnsymbol{footnote}, Weilin Huang % latex2html id marker 4303
\setcounter{footnote}{2}\fnsymbol{footnote}, Yatong Zhou % latex2html id marker 4304
\setcounter{footnote}{3}\fnsymbol{footnote}, Wei Liu % latex2html id marker 4305
\setcounter{footnote}{2}\fnsymbol{footnote}, Dong Zhang % latex2html id marker 4306
% latex2html id marker 4307
\setcounter{footnote}{1}\fnsymbol{footnote}School of Earth Sciences
Zhejiang University
Hangzhou, Zhejiang Province, China, 310027
% latex2html id marker 4308
\setcounter{footnote}{2}\fnsymbol{footnote}State Key Laboratory of Petroleum Resources and Prospecting
China University of Petroleum-Beijing
Beijing, China, 102249
% latex2html id marker 4309
\setcounter{footnote}{3}\fnsymbol{footnote}Tianjin Key Laboratory of Electronic Materials and Devices
School of Electronics and Information Engineering
Hebei University of Technology
5340 Xiping Road, Beichen District
Tianjin 300401, P.R. China


Amplitude-preserving data processing is an important and challenging topic in many scientific fields. The amplitude-variation details in seismic data are especially important because the amplitude variation is directly related with the subsurface wave impedance and fluid characteristics. We propose a novel seismic noise attenuation approach that is based on local plane-wave assumption of seismic events and the amplitude preserving capability of the orthogonal polynomial transform (OPT). The OPT is a way for representing spatially correlative seismic data as a superposition of polynomial basis functions, by which the random noise is distinguished from the useful energy by the high orthogonal polynomial coefficients. The seismic energy is the most correlative along the structural direction and thus the OPT is optimally performed in a flattened gather. We introduce in detail the flattening operator for creating the flattened dimension, where the OPT can be applied subsequently. The flattening operator is created by deriving a plane-wave trace continuation relation following the plane-wave equation. We demonstrate that both plane-wave trace continuation and OPT can well preserve the strong amplitude variation existing in seismic data. In order to obtain a robust slope estimation performance in the presence of noise, a robust slope estimation approach is introduced to substitute the traditional method. A group of synthetic, pre-stack and post-stack field seismic data are used to demonstrate the potential of the proposed framework in realistic applications.