Characteristic wavefield decomposition, imaging and inversion with prestack seismic data
Wang Hua-Zhong,Feng Bo,Liu Shao-Yong,Hu Jiang-Tao,Wang Xiong-Wen,Li Hui
DOI: https://doi.org/10.6038/cjg20150617
2015-01-01
Abstract:The sparse expression of seismic data is important to signal processing, migration and parameter inversion. The conventional linear Radon transform is generalized to invert for both the local plane-wave and shot/receiver ray parameters simultaneously, which can compress the prestack seismic data in the local plane-wave domain. Consequently, it will be convenient to the subsequent ray-based imaging and inversion. Based on the combination of shot/receiver ray parameters, a new plane-wave prediction model is presented which can predict the local plane-wave within the spatial window of local shot/receiver center. The basic assumption is the local linear property of the wave-front. The seismic traces within the spatial window of local shot/receiver center are time-shifted and stacked according to the shot/receiver ray parameters and spatial distance, and the seismic phase with the given shot/receiver ray parameter can be separated naturally. Since the presented time-shift and stack procedure relies on the propagation direction of local plane-wave, we call it characteristic wave decomposition (CWD). Besides, the inversion-based CWD under the framework of compressed sensing (CS) theory can invert for both the local plane-wave and shot/receiver ray parameters simultaneously. Due to the pre-estimated ray-parameters for shot/receiver center, the loop for shot/receiver ray-parameters is avoided. Therefore the efficiency of conventional ray-based migration can be increased remarkably, and the ellipsoid migration response in constant media can be simplified to a point response. The new imaging scheme is called the characteristic wave imaging (CWI). First, the proposed CWD method is tested with the Marmousi synthetic data. Numerical examples show that the seismic phases can be separated completely. Meanwhile, even if the original seismic traces are contaminated by severe random noise, the numerical results demonstrate the robustness of CWD method. Then, the CWI method is tested with the Waxian synthetic data. Due to the high S/N ratio, the seismic data are compressed considerably. Compared with Kirchhoff PSDM method, the imaging efficiency is increased over one order of magnitude. Due to the pre-estimated ray-parameters, the angle domain CIGs can be produced directly. Finally, the CWD method is applied to the separation of first-arrival waveform, serving as the input data for wave-equation based transmission traveltime tomography. The proposed CWD method can compress the prestack seismic data in the local plane-wave domain. Under the framework of compressed sensing theory, both the local plane-wave and shot/receiver ray parameters can be inverted simultaneously. However, the local linear property of wave-front is the basic assumption of CWD. Therefore, it is not suitable for complex topography. Besides, in case of sparse shot/receiver sampling, it is difficult to invert for the shot/receiver ray-parameter simultaneously. The double beam forming procedure in CWD will be degenerated to the classical beam forming method. The presented CWI scheme can map the local plane-wave into the model space with the pre-estimated shot/receiver ray-parameters, so that the imaging efficiency is increased and the migration noise is attenuated. Meanwhile, the angle domain CIGs can be produced easily. Combined with the angle-domain tomography theory, it may be beneficial to the automation of velocity model building procedure.