Denoising and Detecting Seismic Weak Signal Based on Curvelet Thresholding Method

LIU Lei,LIU Zhen,ZHANG Jun-hua
DOI: https://doi.org/10.3969/j.issn.1004-2903.2011.04.037
2011-01-01
Abstract:It is very important for old exploitative oil regions to prospect deeper layers and low-amplitude paleostructures.Inevitably,we must face the problem of detecting weak signals mixed with strong noise because of the complexity of those geologic objects in exploration.With multi-scale characteristic,curvelet transform could remove most of the noise and detect weak signals.Denoising methods in curvelet domain are based on the principle that the curvelet transform has the optimal sparse representation for image edges.Those methods such as thresholding try to obtain a sparser curvelet coefficient matrix which stands for signals and they can get remarkable effect.But when we put a high threshold in curvelet-denoising to suppress strong random disturbance,the signals estimated will have dim breakpoints and anamorphic events because a part of small curvelet coefficients for valid signal are removed.To protect weak wavefront,we introduce a process which can decrease the threshold and retains some small effective coefficients.Firstly,we transform a noising image to curvlet domain and then suppress the small coefficients by thresholding which has a low threshold value.Because of the lower threshold value,many noise coefficients are remained as wild values in curvelet domain.In the second step we clear wild value among coefficients away using median filter and curvlet domain is smoothed.At last we transform the coefficient matrix to time-space domain and obtain a denoised image.The reconstituted image we got has a higher SNR and a fewer loss of weak information compared to former thresholding when it was put to model and real seismic data.We find that this program can remove most noise and protect weak signals when the threshold value selected as half of that of hard thresholding. Experiments prove that the denoising method based on curvelet transform can be combined with other techniques and it may give a better result.
What problem does this paper attempt to address?