2017年12月第17卷第4期廊坊师范学院学报(自然科学版）Journal of Langfang Teachers University(Natural Science Edition)Dec.2017Vol.17 No.4地震数据去噪的研究与分析崔少华，方振国，姜恩华(淮北师范大学，安徽淮北235000)【摘要】小波分解利用信号和噪声在时频域内小波系数差异，重构信号从而去除噪声;SV D算法利用信号和噪声奇 异值的差异，选取合适的空间维数重构信号从而去除噪声。
通过模型验证对比了两种算法的去噪效果和优缺点，发现 小波算法对水平同相轴和倾斜同相轴具有相同的去噪效果，但运算时间较长;SV D算法对水平同相轴去噪效果很好，对于倾斜同相轴几乎不能去噪，但运算时间较短。
【关键词】小波分解;SV D算法;水平同相轴;倾斜同相轴Research and Analysis of Seismic Data DenoisingCUI Shao-hua, FANG Zhen-guo, JIANG En-hua(Huaibei Normal University,Huaibei235000 ,China)【A bstract】The wavelet decomposition uses the difference of signal and noise in the time—frequency domain wavelet coefficient，then reconstructs the signal wiping off noise; SVD algorithm uses the difference of the signal and noise singular value, selecting the appropriate dimensions to reconstmct the signal wiping off noise. Through the comparison of two kinds of model validation algorithm for denoising effect, we know that wavelet algorithm for level phase axis and dip phase axis has the same denoising effect, but the operation time is long. The SVD algorithm for level event denoising effect is very good and the operation time is short, but it does not fit for dip event axis denoising.【K eyw ords】wavelet decomposition; SVD algorithm; horizontal event; dip event〔中图分类号〕P315 〔文献标识码〕A〔文章编号〕1674-3229(2017)04-0033-03〇引目实际中采集到的地震数据通常混入大量随机噪 声，几乎不可避免，因此对地震数据进行随机噪声压 制往往是处理数据的首要工作。
小波分解去噪是在 时频域内处理信号，而SVD则是在时域内处理信号。
本文主要针对上述两种去噪方法进行探讨和分析，从而归纳两种算法的优缺点。
1小波分解与重构去噪法1.1 基本原理小波变换显示的是信号的局部特性，通过伸缩和 平移等运算将信号进行多尺度细化分析r a(MultiscaieAnalysis)。
小波分解与重构去噪法是利用小波变换根据实际需要，将信号(含噪声)分解到同一个尺度下，得 到一系列处于不同频带的小波系数（一酿为高频部分对应信号，低频部分应噪声)。
通过将噪声所在频 带置为零，再进行小M构，从而达到去噪的躲假设/4是/(0的采样数据，即/4=c«u(零尺度空 间小波系数),则可将信号/0)的正交小波变换分解 的公式p写为：C i,k= H C j-lnhn-2kndi,k= H dJ-hnS (2)(1)其中=〈外),‘(0〉表示尺度系数，表示小波系数[4W表示对信号进行小波分解的层数，A分别为一对正交镜像滤波器组，信号离散采样点数用斤表示。
小波经上述分解后的[收稿日期]2017-07-16[基金项目]安徽省高等学校自然科学研究项目（KJ2017B008);淮北师范大学质量工程项目(jy2016133)[作者简介]崔少华（1983-)，女，硕士，淮北师范大学物理与电子信息学院讲师，研究方向：电子与通信。
• 33 •2017年12月廊坊师范学院学报（自然科学版)第17卷•第4期逆过程就是小波重构问过程，即：cj-m =Z (cy >V 2n +^,n^-2n)(2)n1.2模型验证构成1个二维地震记录:总共200道，每道501 个采样点，采样间隔1毫秒，其中一条水平同相轴, 一条倾斜同相轴。
地震子波采用峰频为50H z 的Ricker 子波，每道数据加人幅度为0.1的随机噪声， 如图1所示。
对原始数据进行小波分解与重构，去噪 结果如图2所示，噪声如图3所示。
20406080100120140160180200图2小波重构有效信号偏移/m图3小波分解去除噪声由图2可知，小波分解算法可以将混合数据中的大量随机噪声进行压制，对于倾斜同相轴也可以 达到很好的去噪效果。
但由图3可知,去除的噪声中 还含有部分信号的能量，这是由于小波分解时采用 的Symmlet 8小波算子进行8尺度分解，在重构时对 信号有一定损伤。
可见小波分解算法只要选取合适 的小波分解算子,即可达到理想的去噪效果。
2 SVD 分解去噪算法2.1 SV D 原理假设地震记录为mxri 的数据矩阵X ，其中m 表 示地震波通道观察个数，n 表示每个通道上的采样 点数。
将矩阵Z 做分解变换为〜X = U Y .V 1(3)mxnmxmmxnnxn式(3)中[/与F 均为正交归一阵，E 为准对角矩阵，在的情况下，有〇\ 0 …0 0•"0^cr 2 00 〇•"〇m x n： ： *.: : :0 ••• 〇■_.0**•0L m ⑷一般设0^(72称C T ;为数据的X奇异值。
x ^i^iuiv,T (5)/=1式中&和^称为左、右奇异矢量'其维数分别 为 m x l 和 r exl 〇因此，维数为/n X r i 的原始数据就被分解成维 数相同的ni 个子矩阵,这就是奇异值分解'确定有效信号空间维数心进行选择前个本 征图像来进行重构就可以恢复有效信号：FP(s ) = a siuiv，T(6)SVD 分解算法在去噪时的关键是对有效信号 进行重构，所选取的重构信号维数P ,即加权本征 图像的个数，对于去噪结果有很大的影响' 若选用P 的维数小于有效信号的维数，则无法重构信号且 在去除噪声中包含信号能量，即对信号有损伤;若 选用P 的维数大于信号的维数，则会引人过多的随 机噪声，在重构的信号中含有大量噪声未被压制，去噪效果不佳。
因此重构维数P 的选择是SVD 算法 的关键。
• 34•第17卷•第4期崔少华等:地震数据去噪的研究与分析2017年12月2.2模型验证采用图1的含噪原始数据进行SVD分解，由于 数据中同相轴仅含有两种倾角，因此重构维数/>=2, 重构结果如图4所示，噪声如图5所示。
偏移/m图5 SVD去除噪声由图4可知，SVD重构信号中只含有水平同相轴，不含有任何倾斜同相轴，这是由于SVD算法 是在时域进行分解，利用各道数据之间的相关性将 地震记录分解成奇异值并按序排列，信号的有效性 越强，分解出奇异值差异越明显w，越利于提取信 号的奇异值进行重构去噪。
由于倾斜同相轴各道 的相关性减弱(倾角越大，相关性越弱），因此SVD 算法无法对其进行重构。
从图5很明显地看出，倾 斜同相轴完全没有被重构，噪声中却不含有任何水 平同相轴的能量，所以SVD时域分解对倾斜同相轴 有严重的损害。
3小波算法与SVD算法的对比3.1 PSNR的对比采用图像的峰值信噪比（Peak Signal to Noise Ratio,PSNR)对两种算法对比，结果如表1所示。
表1两种算法PSNR对比小波分解算法PSNR=75.7284dBSVD算法PSNR=67.5590dB由表1可知，小波分解算法的PSNR高出SVD 算法的接近8dB，这说明前者还原图像的程度高于 后者，即前者的去噪效果优于后者的去噪效果。
3.2去噪时间的对比在实际中往往处理的是大规模地震记录，因此 X才于处理时间的要求同样严格。
将两种算法的去噪 时间进行对比(基于验证模型），结果如表2所示。
表2两种算法去噪时间对比小波分解算法23.6120sSVD算法8.3830s3.3去噪效果的对比小波算法在分解和重构信号时对有效信号和噪 声采用同一个分解尺度，因此对水平同相轴和倾斜 同相轴去噪效果相同（图2)，同时小波分解算法对 于倾斜同相轴具有良好的去噪效果(图3)〇SVD算法对水平同相轴的去噪效果最佳（图 4)，对倾斜同相轴几乎不能去噪(图5)。
4总结综上所述，小波分解算法和SVD算法各有优缺 点，决定小波分解算法去噪效果的关键因素是小波 算子分解尺度的选取，分解尺度越大，去噪越明显。
但尺度选取过大，在去噪时会引入过多随机噪声，反 而会使去噪效果下降。
决定SVD算法去噪效果的关键因素是重构维 数的选取，且重构信号维数的取值需人工判断，在大 数据量的地震记录中，每个频率切片上对重构维数 均需进行判断和选取。
[参考文献][1]王香云.基于信号去噪的小波算法研究[J].太原师范学院学报（自然科学版）,2015,14(1) :33-37.[2]王超，沈斐敏小波变换在探地雷达弱信号去噪中的研究[J].物探与化探，2015，39(2): 421-424.[3] 崔少华，单巍•基于小波分析的地震资料去噪方法的研究和应用[J].淮北师范大学学报(自然科学版),2016,37⑶:434& [4]蔡正保.一种基于混池加密和小波变换的数字水印技术研究[J]•廊坊师范学院学报(自然科学版),2016,16(2): 17-21.(下转第42瓦）• 35•2017年12月廊坊师范学院学报（自然科学版)第17卷•第4期X 0003.1常开触点恢复高电平，松开+X 按键，G 116.0 常为“1”，再按下RESET 复位按键后，报警解除。
X 0007.3和X 0010.2分别表示+Y 按键和+Z 按键，同 理当出现Y 轴或Z 轴负超程报警时，分析方法类似。