专利名称：一种对数字信号进行仿真加噪的处理方法
技术领域：
本发明属于地震勘探资料处理等数字信号处理技术领域，尤其是针对在地震勘探 资料处理以及雷达、通讯、音响、图像等数字信号处理领域中，满足信号仿真处理的需要。
背景技术：
在地震勘探资料处理以及雷达、通讯、音响、图像等数字信号处理领域，加噪声是 信号仿真处理的需要。现有加噪方法分为两大类，一类是加白噪声，另一类是加有色噪声。白噪声是在无限频率范围内功率密度为常数的随机噪声信号，具有任意两个样本 都不相关的性质，从某种意义上说它表示了数据的随机性。有色噪声则是随着信号频率的 不同功率密度也变化的随机噪声信号，根据对不同频率段的敏感程度可区分，常见的有粉 红噪声、红噪声、橙色噪声、蓝噪声、紫噪声、灰色噪声、棕色噪声和黑色噪声(静止噪声)。 目前对噪声的研究还是认识的阶段，而对新噪声的合成几乎是一个空白。在地震勘探资料处理中，一般需要压制噪声，提高信噪比。特别是针对规则噪声， 如多次波、散射波、面波等，通常需要利用多维滤波方法进行消除或压制。但是，多维滤波方 法会产生一定的混叠效应，表现之一是输出时间剖面过于呆板。因此，对多维滤波后的信号 道集作适当的仿真加噪处理是必要的。而现有技术中对于加噪的处理往往是加入白噪或有 色噪声，均不能真实体现或还原原有波形制，仿真度低。现有技术中，S' At)为常规加噪方法直接将白噪信号道集加入目标信号道集 (图 3)，一般表达式为S' ,(t) =Si(t) + yNi(t)(1)其中，Sdt)为 需要加噪仿真处理的目标信号道集(图1)，Ni (t)为白噪信号道集(图2)。其中μ为比例 系数，t为时间序列，i为道序列。
发明内容
本发明为了解决现有技术中存在的技术问题，生成了一种新噪声，是一种比较自 然、切合实际的随机噪声。用新噪声加噪后的信号道或信号道集具有极高的仿真度。本发 明的目的是为数字信号处理领域提供一种用于加噪仿真处理的新噪声生成方法，其过程如 下—种对数字信号进行仿真加噪的处理方法，所述方法对目标信号进行仿真加噪处 理首先将目标信号与白噪信号进行褶积运算生成变色信号，后将变色噪声信号加入目标 信号中，并对加噪后的信号进行频谱分析，得到仿真加噪频谱特征；所述的变色噪声信号为所述目标信号和白噪信号进行褶积后的数字信号。所述变色噪声信号为Ni (t)与Si (t)进行褶积后生成一种新随机噪声-变色噪声 N;(t)，如图4，其表达式为<(0 = ^,(0*^(0 (2)； 很显然，(X)并不完全是有色噪声类型，当Si(t)也是白噪信号道时，仍是白噪信号道。A^(X)是在白噪声信号与有色噪声信号之间随Si(t)的种类变化而变化的新一类随机噪声。因此，将其命名为变色噪声。将所述变色噪声信号道集加入目标信号道集，表达式为S: (t) = SiHhjuN: it)⑶其中，$Λ0)为已加噪仿真处理的目标信号道集，如图5。Si (t)为需要加噪仿真处理的目标信号道集，Ni (t)为白噪信号道集；μ为比例系 数，t为时间序列，i为道序列。所述仿真加噪处理方法采用时间域包括以下步骤步骤1，采集目标数字信号或目标数字信号道集；步骤2，生成白噪信号或白噪信号道集；步骤3，生成变色噪声信号或信号道集<(0 = ^(0*^(0 (2)步骤4，将所述步骤3得到的变色噪声信号输入目标数字信号或信号道集，进行加 噪的处理；步骤5，输出加噪后的数字信号。所述仿真加噪处理方法采用频率域包括以下步骤步骤1，采集目标数字信号或目标数字信号道集；步骤2，生成白噪信号或白噪信号道集；步骤3，对步骤1和步骤2中的目标数字信号或信号道集以及白噪信号或道集进行傅氏变换，得到Si (ω) = FFT (Si (t)}(4)Ni (ω) = FFT(Ni (t)}(5)步骤4，生成变色噪声信号或信号道集N: O) = Ni (ω) ■ Si (ω)(6)步骤5，对生成的变色噪声信号或道集进行傅氏反变换N； (O = FFT'1 {N； (ω)}(y)步骤6，将所述步骤5得到的变色噪声信号输入目标数字信号或信号道集， 进行加噪的处理；步骤7，输出加噪后的数字信号。本发明应用在地震勘探领域中，所述目标数字信号道集采集于地震数据，并在地 震勘探资料处理中，多维滤波后的信号道集为所述目标数字信号道集。本发明方法应用于在地震勘探资料处理中，对地震波进行处理过程中，对多维滤 波后的地震数字信号进行仿真加噪处理。通过比对加入本发明的变色噪声后的频谱输出，以及现有技术中的加入白噪和有 色噪声的频谱，分析结果表明，变色噪声是一种自然、切合实际的人工合成随机噪声，用变 色噪声加噪后的目标信号道或信号道集具有极高的仿真度。
图1需要加噪仿真处理的目标信号道集(理论模型)；图2白噪信号道集；图3直接将白噪信号道集加入目标信号道集；图4生成的一种新随机噪声-变色噪声道集；图5变色噪声道集加入目标信号道集；图6原始目标信号道集频谱；图7白噪信号道集频谱；图8直接将白噪信号道集加入目标信号道集的频谱；图9变色噪声道集的频谱；图10将变色噪声道集加入目标信号道集的频谱；图11某地震勘探工区的一组原始CMP道集及其速度谱；图12原始CMP道集频谱；图13消除多次波干扰后的CMP道集及其速度谱；
图14去噪后的CMP道集频谱；图15增加白噪声后的CMP道集及其速度谱；图16增加白噪声后的CMP道集频谱；图17增加白噪声的CMP道集经带通滤波后的频谱；图18增加30%自然噪声后的CMP道集及其速度谱；图19增加30%变色噪声后的CMP道集及其速度谱；图20增加变色噪声后的CMP道集频谱；图21本发明的时间域实施流程图；图22本发明的频率域实施流程图。各幅附图将结合说明书内容加以详细说明
具体实施例方式本发明用Ni (t)与Si (t)进行褶积运算，生成一种新随机噪声-变色噪声(图4)。 变色噪声的表达式为N;(t) = NXt)*SXt)(2)很显然，Λ^Ο)并不完全是有色噪声类型，当Si (t)也是白噪信号道时，A^仍 是白噪信号道。A^O)是在白噪声信号与有色噪声信号之间随Si (t)的种类变化而变化的 新一类随机噪声。将此变色噪声道集加入目标信号道集，一般表达式为=+(3)其中，为已加噪仿真处理的目标信号道集，如图5。图6是SJt)的频谱，图7是Ni(t)的频谱。常规加噪的方法是直接将白噪Ni (t)加入目标信号Sjt)，图3和图8分别是S' dt)的时间剖面和频谱图。在时间剖面中存在 明显的加噪处理现象，而从频谱中也可以发现加入的白噪均勻地分布于整个频率域。因此， 常规加噪处理的仿真度较低。如图4所示，本发明所生成的变色噪声Λ^Ο)仍是随机噪声，在时间剖面中 的能量分布特征符合Si (t)的能量分布特征，的频谱特征(图9)同样符合Si (t)的 频谱特征。因此，无论是的时间剖面中，还是在&频谱(图10)中，都几乎看不 出加噪处理的现象。分析结果表明，变色噪声是一种比较自然、切合实际的人工合成随机噪声，用变色 噪声加噪后的目标信号道或信号道集具有极高的仿真度。在地震勘探资料处理中，一般需要压制噪声，提高信噪比。特别是针对规则噪声， 如多次波、散射波、面波等，通常需要利用多维滤波方法进行消除或压制。但是，多维滤波方 法会产生一定的混叠效应，表现之一是输出时间剖面过于呆板。因此，对多维滤波后的信号 道集作适当的仿真加噪处理是必要的。图11是某地震勘探工区的一组CMP道集及其速度谱，图中可见，在3500ms以下出 现严重的多次波干扰，图12是其CMP道集频谱。图13是利用某种多维滤波方法对图11消 除多次波干扰后的CMP道集及其速度谱，图14是去噪后的CMP道集频谱，观察速度谱可知， 多次波干扰已被消除，但是CMP道集已类似合成模型，感觉很不自然。图15是在图13的基础上，增加白噪声的结果，无论是在时间剖面中还是在频谱 (图16)中都存在明显的加白噪声现象。带通滤波(5，10，60，80Hz)可隐藏时间剖面中的加 噪现象，但输出频谱(图17)仍然显示出加白噪声的痕迹。增加原CMP道集中的自然噪声(有色噪声)是目前最常用的加噪方法。图18是 在图13的基础上，增加30%自然噪声的结果，其中，CMP道集中的多次波干扰有所减弱，但 在速度谱中仍有较严重的多次波能量显示。图19是本发明以图13为目标信号道集，增加30%变色噪声的结果，图20是加噪 后的目标信号道集频谱，其中，CMP道集的时间剖面感觉很自然，其速度谱也无多次波显示， 频谱对比可见，增加变色噪声没有改变其频谱特征。因此，无论是在时间剖面中还是在频谱 中都很难看出有加噪处理的现象。图21是本发明的时间域实施流程图，图22是本发明的频率域实施流程图。两种 实施例的结果是等同的。实际应用结果表明，用变色噪声加噪后的目标信号道集的特点是在时间域中有明 显的噪声显示，而在频率域中没有明显的噪声显示。这对数字信号处理中解决噪声压制和 仿真加噪问题非常有帮助。变色噪声是一种比较自然、切合实际的人工合成随机噪声。用 变色噪声加噪后的目标信号道，无论是在时间剖面中还是在频谱中都很难看出有加噪处理 的现象。因此，变色噪声具有极高的仿真度。
权利要求
1.一种对数字信号进行仿真加噪的处理方法，其特征在于，所述方法对目标信号进行 仿真加噪处理首先将目标信号与白噪信号进行褶积运算生成变色信号，后将变色噪声信 号加入目标信号中，并对加噪后的信号进行频谱分析，得到仿真加噪频谱特征；所述的变色 噪声信号为所述目标信号和白噪信号进行褶积后的数字信号。
2.根据权利要求1所述的一种数字信号进行仿真加噪处理方法，其特征在于， 所述变色噪声信号为队⑴与Si(t)进行褶积后生成一种新随机噪声-变色噪声N;(t),其表达式为将所述变色噪声信号道集加入目标信号道集，表达式为 S;(t) = SXt) + MN;(t)(3)其中，Si(t)为需要加噪仿真处理的目标信号道集，NJt)为白噪信号道集；μ为比例系 数，t为时间序列，i为道序列。
3.根据权利要求1所述的一种数字信号进行仿真加噪处理方法，其特征在于， 所述仿真加噪处理方法采用时间域包括以下步骤步骤1，采集目标数字信号或目标数字信号道集； 步骤2，生成白噪信号或白噪信号道集； 步骤3，生成变色噪声信号或信号道集步骤4，将所述步骤3得到的变色噪声信号输入目标数字信号或信号道集，进行加噪的 处理；步骤5，输出加噪后的数字信号。
4.根据权利要求1所述的一种数字信号进行仿真加噪处理方法，其特征在于， 所述仿真加噪处理方法采用频率域包括以下步骤步骤1，采集目标数字信号或目标数字信号道集； 步骤2，生成白噪信号或白噪信号道集；步骤3，对步骤1和步骤2中的目标数字信号或信号道集以及白噪信号或道集进行傅氏 变换，得到Si(Co) = FFTiSi (t)}(4)Ni(Co) = FFTiNi (t)}(5)其中ω为圆频率；步骤4，生成变色噪声信号或信号道集Ν,Λ (ω) = Νχω)-SXa)(6)步骤5，对生成的变色噪声信号或道集进行傅氏反变换N (t) = FFT1 {N (ω)}(ν)步骤6，将所述步骤5得到的变色噪声信号TV,Λ( )输入目标数字信号或信号道集，进行 加噪的处理；步骤7，输出加噪后的数字信号。
5.根据权利要求1-4之一所述的一种数字信号进行仿真加噪处理方法，其特征在于， 所述目标数字信号道集采集于地震数据，并在地震勘探资料处理中，多维滤波后的信号道集为所述目标数字信号道集。
6.应用如权利要求1-5之一的仿真加噪处理方法，其特征在于，所述方法应用于在地震勘探资料处理中，对地震波进行处理过程中，对多维滤波后的 地震数字信号进行仿真加噪处理。
全文摘要
一种对数字信号进行仿真加噪的处理方法，所述方法对目标信号进行仿真加噪处理首先将目标信号与白噪信号进行褶积运算生成变色信号，后将变色噪声信号加入目标信号中，并对加噪后的信号进行频谱分析，得到仿真加噪频谱特征；所述的变色噪声信号为所述目标信号和白噪信号进行褶积后的数字信号。变色噪声是一种比较自然、切合实际的人工合成随机噪声。用变色噪声加噪后的目标信号道的特点是在时间域中有明显的噪声显示，而在频率域中没有明显的噪声显示，无论是在时间剖面中还是在频谱中都很难看出有加噪处理的现象。因此，变色噪声具有极高的仿真度。
文档编号G01V1/36GK102043168SQ20101050799
公开日2011年5月4日 申请日期2010年10月15日 优先权日2010年10月15日
发明者刘志成, 孙武亮, 杨勤勇, 段心标, 谢金娥 申请人:中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院