曲波域统计量自适应阈值探地雷达数据去噪技术

非线性、非平稳探地雷达数据常掺杂各种复杂噪声源,其对精确提取弱反射波信号、识别绕射波双曲线同相轴特征具有严重影响,忽略噪声影响给探地雷达探测数据全波形偏移成像及后续解译造成较大误差.采用传统阈值函数的曲波变换去噪需要根据数据噪声水平人为确定合理阈值控制系数.对此,本文开展自适应阈值函数的曲波变换去噪算法研究.引入块状复数域阈值函数算法,分析传统阈值函数曲波变换去噪的效果随阈值控制系数变化的规律;利用高阶统计量理论,对曲波变换系数在尺度、方向上进行相关性叠加,通过相关性统计量自适应确定有效信号在曲波变换系数分布尺度、旋转方向,由此确定清除噪声成分阈值范围,构建统计量自适应阈值函数曲波变换去噪算法.针对包含随机噪声、相关噪声合成探地雷达数据及实测探地雷达数据,采用传统阈值函数曲波变换去噪与本文提出去噪算法处理结果对比分析,检验了本文算法的有效性及可行性.研究成果对复杂探地雷达数据精确推断解译具有指导意义.

Scale and rotation statistic-based self-adaptive function for ground penetrating radar denoising in curvelet domain

Nonlinear and non-stationary ground penetrating radar (GPR) data for geophysics exploration are often mixed with various complex noise sources, such as random and coherent noise. Those complex noise sources are introduced by the acquisition system and other sources of measurement uncertainty. The data sets which are contaminated by the above noises have a serious influence on the accurate extraction of weak reflected signals and the effective identification of diffracted wave hyperbolic coaxial characteristics. The ignorance of the influence of noise will cause great errors in the interpretation of GPR data and the subsequent migration imaging with full waveform. The curvelet transform not only is related to position and frequency as compared with the wavelet transform, but also is controlled by the translation angle. With such a unique advantage, curvelets are used for ground roll whitening, coherent noise denoising and separation of overlapping events. The traditional curvelet transform denoising with a hard threshold function needs to give a reasonable threshold function control coefficient according to the noise level of GPR data. As a result, an appropriate choice of a hard threshold is a basic requirement, and this presents a challenging task in curvelet denoising. To overcome these shortcomings, an self-adaptive threshold function for GPR data denoising with curve transform is proposed in this paper. For detailing the reasonable control coefficient of the threshold function, the complex block threshold function algorithm is used to analyze the distribution of peak-signal-to-noise ratio value of the noisy synthetic GPR data contaminated with random and coherent noise by using the traditional threshold function curvelet transform. Based on the high order statistical theory, the accuracy distribution of the curvelet coefficient for useful signals in scale and rotation direction are correlatively stacked and determined by using the statistical self-adaptive function. And then the threshold range of noise removal components is given by the statistical self-adaptive function. The effectiveness of the proposed denoising algorithm for the noisy synthetic contaminated with different types of noise (i.e., Gaussian random and coherent), and field GPR data is demonstrated by comparing the denoising results via curvelet transform with those from traditional thresholding function. The presented denoising results by the statistical self-adaptive function is helpful and instructive for the accurate inference and interpretation of complex GPR data.