一种地震数字信号的复数道集二维滤波方法技术

本发明专利技术属于地震勘探资料处理以及电子信号处理等数字信号处理技术领域，尤其涉及一种对数字信号二维滤波的方法。将采集到的地震数字信号的数据道集，对信号数据道集进行希尔伯特变换，得到所述地震数字信号的余弦相位函数道集和瞬时包络道集；后对所述余弦相位函数道集实施二维滤波步骤，经二维滤波后的余弦相位函数道集与瞬时包络道集进行各道对应相乘处理，将消除噪音及干扰信号，同时消除混叠现象后的地震信号数据道集输出，用于对地震数字信号的保真处理目的。

【技术实现步骤摘要】

本专利技术属于地震勘探资料处理以及电子信号处理等数字信号处理
，尤其 涉及一种对数字信号的二维滤波方法。
技术介绍
现有技术中，对地震数字信号的二维滤波方法是基于实数道集的二维滤波方法， 如τ-ρ变换滤波、双曲变换滤波和f-k变换滤波等等。作为公知技术，二维滤波目的是压 制随机噪音和规则干扰信号。在数字信号处理领域，特别在地球物理地震资料处理中，实数 道集的二维滤波方法得到广泛的应用。但是，实数道集的二维滤波方法在提高信噪比的同 时会产生不同程度的混叠现象，信号失真和假同相轴的出现是混叠现象的主要特征，压制 或消除混叠现象是一难题。
技术实现思路
本专利技术为了解决现有技术中存在的信号失真和混叠现象，研发了一种地震数字信 号的复数道集二维滤波方法。本专利技术所提出的一种二维滤波实施方法，能够在压制随机噪 音和规则干扰信号的同时消除混叠现象，达到保真处理的目的。本专利技术为了实现上述专利技术目的，所采用的技术方案如下，所述滤波方法为数字信号的复数道集二维滤波方法；将采集到的地震数字信号的 数据道集，对信号数据道集进行希尔伯特变换，得到所述地震数字信号的余弦相位函数道 集和瞬时包络道集；后对所述余弦相位函数道集实施二维滤波步骤，经二维滤波后的余弦 相位函数道集与瞬时包络道集进行各道对应相乘处理，将消除噪音及干扰信号，同时消除 混叠现象后的地震信号道集输出，用于对地震数字信号的保真处理的目的。所述滤波方法包括如下步骤，①.地震数字信号数据道集输入步骤；②.分解地震数字信号步骤希尔伯特变换分解信号；设二维复函数Ci(t) = r, (t) +Jhi (t)(1)式中，hjt)称为实数道集ri(t)的希尔伯特变换道集，i = 1，2，3...11为道序号， t为时间序列，j为虚系数；则瞬时包络道集化0)狀⑷+ h]{t)(2)(瞬时相位道集为雄、=arccos(3)W) J八，、 r (f),、相位函数道集为cos650) 二“^(4)W)于是巧(t)= COS θ J (t) · Bi (t)(5)3③.对相位函数道集实施二维滤波步骤；相位函数与信号振幅大小无关，余弦相位函数的数值范围在[-1，1]之间；无论是 强信号同相轴还是弱信号同相轴在余弦相位函数道集中的成像效果都是一致的(见图4)。对相位函数道集COS θ i (t)实施二维滤波COS <9, (0 二 TDF{ COS6>,(0 }(6)式中，TDF(Tm)-DFilter)为二维滤波方法，COS θ. (t)% cos θ j (t)经二维滤波后 的道集；④.合成输出信号步骤；上述公式(5)表明实数道集ri(t)被分解为相位函数道集cos θ , (t)和瞬时包络 道集ajt)的积；瞬时包络具有下列特征a、与信号的相位和频率无关；b、幅值大小反映信号的振幅变化；c、宽度等于信号的延续长度。瞬时包络锁定了信号的振幅和延续长度(见图1)。相反，混叠现象破坏了信号的 振幅和延续长度。因此，可用瞬时包络来抵消混叠现象。于是φ) 二 COS 0X0 . Qi (0(7)公式(7)为滤波后的实函数道集_被分解为滤波后的相位函数道集cosOi与瞬时包络道集屮⑴的积；。其中，@是经复数道集二维滤波后的无混叠现 象的最终输出。上述步骤3中的二维滤波方法采用τ -P变换滤波、双曲变换滤波或f_k变换滤波 中的一种。本专利技术中的研发理念，如表1所示。\集 信COS0.(t)COS 明)rXt) = cos θXt) · a^t)有效有有有有有干扰有有有无无混叠无无无有无表1复数道集二维滤波的技术理念本专利技术结果表明，应用效果明显，其输出结果在有效地压制干扰信号的同时基本 消除混叠现象，达到保真处理的目的。与现有技术的实数道集二维滤波方法相比，复数道集 二维滤波方法更具有优越性。附图说明图1本专利技术实信号的瞬时包络图；图2现有技术理论模型的CMP实数道集；图3对图2进行实数道集双曲变换滤波结果；图4对图2进行希尔伯特变换获得的余弦相位函数道集；图5对图2进行希尔伯特变换获得的瞬时包络道集；图6对图4的余弦相位函数道集进行双曲变换滤波结果；图7图6与图5各道对应相乘的最终输出结果；图8本专利技术实施例-某测线的CMP实数道集；图9对图8进行实数道集τ -ρ变换滤波结果；图10采用本专利技术对图8进行复数道集τ -ρ变换滤波结果；图11本专利技术实施例-某测线的CMP实数道集及其速度谱；图12对图11中CMP实数道集进行复数道集双曲变换滤波结果及其速度谱；图13本专利技术二维滤波方法流程图。上述各幅附图结合
技术实现思路
部分和具体实施部分加以说明具体实施例方式图2是某理论模型的CMP (共中心点)实数道集，道集中存在直达波和折射波等规 则线性干扰波，因此，选择双曲变换滤波是合适的。我们首先采用现有技术直接对图2的实 数道集进行双曲变换滤波，滤波结果如图3，显示线性干扰波已被消除，但存在信号失真现象。采用本专利技术技术，对图2进行希尔伯特变换分别获得的余弦相位函数道集如图4， 和瞬时包络道集，如图5，然后对余弦相位函数道集进行双曲变换滤波，最后将滤波后的余 弦相位函数道集，如图6，与图5瞬时包络道集进行各道对应相乘处理，图7为输出结果。表 明已达到了双曲变换滤波效果，并基本消除了信号失真现象。图8是某测线的CMP实数道集，采用现有技术直接对图8的实数道集进行τ -ρ变 换滤波，图9为滤波结果，显示实数道集被基本还原，但出现假同相轴现象。图10是采用本专利技术技术对图8进行复数道集τ -ρ变换滤波的结果，实数道集被 基本还原，但是没有出现假同相轴现象。图11是某测线的CMP实数道集及其速度谱，由于该实数道集的信噪比较低，速度 谱有效波能量团不够清晰，影响到速度分析精度。图12是对图11中CMP实数道集进行复 数道集双曲变换滤波结果及其速度谱，道集的信噪比有较大提高，且基本无混叠现象，速度 谱有效波能量团清晰可见。试验结果表明，本技术应用效果明显，其输出结果在有效地压制干扰信号的同时基本消除混叠现象，达到保真处理的目的。与现有技术的实数道集二维滤波方法相比，复数 道集二维滤波方法更具有优越性。 上述技术方案只是本专利技术的一种实施方式，对于本领域内的技术人员而言，在本 专利技术公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本 专利技术上述具体实施方式所描述的结构，因此前面描述的方式只是优选地，而并不具有限制 性的意义。本文档来自技高网...

【技术保护点】
一种地震数字信号的复数道集二维滤波方法，其特征在于，所述滤波方法为数字信号的复数道集二维滤波方法；将采集到的地震数字信号的数据道集，对信号数据道集进行希尔伯特变换，得到所述地震数字信号的余弦相位函数道集和瞬时包络道集；后对所述余弦相位函数道集实施二维滤波步骤，经二维滤波后的余弦相位函数道集与瞬时包络道集进行各道对应相乘处理，将消除噪音及干扰信号，同时消除混叠现象后的地震信号数据道集输出，用于对地震数字信号的保真处理目的。

【技术特征摘要】
1.一种地震数字信号的复数道集二维滤波方法，其特征在于，所述滤波方法为数字信 号的复数道集二维滤波方法；将采集到的地震数字信号的数据道集，对信号数据道集进行 希尔伯特变换，得到所述地震数字信号的余弦相位函数道集和瞬时包络道集；后对所述余 弦相位函数道集实施二维滤波步骤，经二维滤波后的余弦相位函数道集与瞬时包络道集进 行各道对应相乘处理，将消除噪音及干扰信号，同时消除混叠现象后的地震信号数据道集 输出，用于对地震数字信号的保真处理目的。2.根据权利要求1所述的一种地震数字信号的复数道集二维滤波方法，其特征在于， 所述滤波方法包括如下步骤，①.地震数字信号数据道集输入步骤；②.分解地震数字信号步骤希尔伯特变换分解信号；设二维复函数Ci (t) = r, (t) +Jhi (t) (1)式中，比⑴称为实数道集巧⑴的希尔伯特变换道集，i = 1，2，3...1!为道序号3为 时间序列，j为虚系数；则瞬时包络道集《( 二 ^r2(O+ Ai2 (O、 ， SS'rα: W瞬时相位道集为雄、=ar...

【专利技术属性】
技术研发人员：刘志成，杨勤勇，吴鹏，张兵，段心标，邵志东，傅星菊，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：11[中国|北京]