【技术实现步骤摘要】
本专利技术涉及油田地震资料处理，特别是涉及到一种迭代法引导滤波压制随机噪音方法。
技术介绍
1、引导滤波是由何凯明等人与2010年提出，此方法有良好的边缘保持特性，且不会出现梯度反转现象，可广泛应用于图像降噪、去雾、高动态范围压缩等领域。
2、在地震数据处理领域，引导滤波应用较少，主要原因是缺乏较高信噪比的数据作为引导图像。
3、在申请号：cn201310522885.7的中国专利申请中，涉及到全波形反演噪音滤波算子的提取及使用其噪音滤波的方法。该专利技术在全波形反演中获取单炮的全波形反演梯度之后，对二维梯度进行二维傅里叶变换将数据从空间域转化到波数域；针对全波形反演速度噪音在波数域中的可存在区间，以及有效全波形反演速度在波数域中的分布，确立波数域中垂直波数与水平波数的关系；采用全波形反演噪音滤波算子避免数据在空间域与波数域之间转换产生的吉普斯现象；最后利用二维傅里叶反变换得到适应于地质特征的具有高信噪比、高保真度的全波形反演结果。该专利技术能够结合地质构造特征确定全波形反演噪声存在区间，保持有效全波形反演速度的同时使反演噪音得以消除，效果稳定可靠。
4、在申请号：cn202011094102.6的中国专利申请中，涉及到一种可控震源黑三角噪音压制方法及系统，属于可控震源激发地震资料数据处理领域。所述可控震源黑三角噪音压制方法首先对地震数据进行傅里叶变换，然后采用相关法在不同的频带范围内求解噪音压制算子矩阵，最后采用稀疏约束反演迭代算法，实现可控震源黑三角噪音的有效压制。该专利技术采用稀疏约束反
5、在申请号：cn201410416113.x的中国专利申请中，涉及到一种用于抑制全波形反演震源编码串扰噪音的方法和系统。所述方法包括：对反演迭代过程中更新的速度模型进行低通滤波；对低通滤波之后的速度模型做主分量分析，确定地层地质构造的主信号方向；沿主信号方向做定向滤波以抑制与主信号方向不同的震源编码串扰噪音。该专利技术利用全波形反演速度模型中信号与噪音的方向差异，在有效信号的方向上滤波，既能保护有效信号的完整性，又能去除噪声。
6、在申请号：cn201210109770.0的中国专利申请中，涉及到一种多次波压制方法，属于油气及煤层气地震勘探与开发领域。本专利技术方法包括：(1)选择地震道；(2)截取数据；(3)确定多次波的周期；(4)定义目标多次波时窗；(5)计算目标多次波得滤波算子；(6)褶积与滤波；(7)自相关检测滤波结果；(8)迭代处理；(9)将原始地震数据中的其它所有地震道的地震道数据依次作为选出的地震道数据x(t)，对每个选出的地震道数据x(t)重复步骤(4)至步骤(8)。该专利技术方法针对每个多次波调节衰减系数的值，获得了满意的滤波效果，有效地压制了多次波；本专利技术方法适用于二维地震数据和三维地震数据。
7、以上现有技术均与本专利技术有较大区别，未能解决我们想要解决的技术问题，为此我们专利技术了一种新的迭代法引导滤波压制随机噪音方法。
技术实现思路
1、本专利技术的目的是提供一种可以通用于叠前、叠后地震数据，提高了地震资料的信噪比的迭代法引导滤波压制随机噪音方法。
2、本专利技术的目的可通过如下技术措施来实现：迭代法引导滤波压制随机噪音方法，该迭代法引导滤波压制随机噪音方法包括：
3、步骤1，输入地震数据作为输入图像；
4、步骤2，应用引导滤波，将地震数据自身作为引导图像，以后均以上次去噪后数据作为引导图像，进行引导滤波，提高资料信噪比；
5、步骤3，将随机噪音进行压制；
6、步骤4，输出随机噪音压制后的地震数据。
7、本专利技术的目的还可通过如下技术措施来实现：
8、在步骤2，应用引导滤波，将地震数据p自身作为引导图像i，以后均以上次去噪后数据作为引导图像i，将地震数据p进行引导滤波，提高资料信噪比。
9、步骤2包括：
10、21)计算引导图像i的均值meani，地震数据p的均值meanp，引导图像i的相关系数corri，引导图像i与地震数据p的相关系数corrip；
11、22)计算引导图像i的方差vari，引导图像i与地震数据p的协方差covip；
12、23)计算线性相关参数a和b；
13、24)计算a的均值meana，b的均值meanb；
14、25)计算引导滤波结果q。
15、在步骤21，利用下列公式，分别计算引导图像i的均值meani，地震数据p的均值meanp，引导图像i的相关系数corri，引导图像i与地震数据p的相关系数corrip；
16、meani＝fmean(i))
17、meanp＝fmean(p)
18、corri＝fmean(i.*i)
19、corrip＝fmean(i.*p)
20、fmean是算数均值.*表示矩阵乘法。
21、步骤22，利用下列公式，分别计算引导图像i的方差vari，引导图像i与地震数据p的协方差covip；
22、vari＝corri-meani.*meani
23、covip＝corrip-meani.*meamp。
24、在步骤23，利用下列公式，分别计算线性相关参数a和b；
25、a＝covip./(vari+∈)
26、b＝meanp-a.*meani
27、a和b为线性相关参数，vari为引导图像i的方差，cvip为引导图像i与地震数据p的协方差，meani为引导图像i的均值，meanp为地震数据p的均值，∈为界定平滑区域和边缘区域的阈值。
28、在步骤24，利用下列公式，分别计算a的均值meana，b的均值meanb；
29、meana＝fmean(a)
30、meanb＝fmean(b)。
31、在步骤25，利用下列公式，计算引导滤波结果q；
32、q＝meana.*i+meanb。
33、在步骤3，以去噪后数据作为输入图像p，重复步骤2，直至迭代次数达到预设值，随机噪音得到很好压制。
34、在步骤3，预设值为地震数据自身。
35、本专利技术中的迭代法引导滤波压制随机噪音方法，采用迭代引导滤波的思路，实现了地震数据随机噪音压制，提高了噪音压制精度，提升了处理效率。本专利技术的方法可以通用于叠前、叠后地震数据，提高了地震资料的信噪比，为后续处理及解释过程提供了有效的技术支撑和保障。与现有的地震资料校正技术相比，本专利技术主要有两点优势：
36、(1)利用引导滤波的边缘保持特性，在去噪的同时，保持断层等特殊地质体的成像效果；
37、(2)在去噪的过程中，引入了本文档来自技高网...
【技术保护点】
1.迭代法引导滤波压制随机噪音方法，其特征在于，该迭代法引导滤波压制随机噪音方法包括：
2.根据权利要求1所述的迭代法引导滤波压制随机噪音方法，其特征在于，在步骤2，应用引导滤波，将地震数据p自身作为引导图像I，以后均以上次去噪后数据作为引导图像I，将地震数据p进行引导滤波，提高资料信噪比。
3.根据权利要求2所述的迭代法引导滤波压制随机噪音方法，其特征在于，步骤2包括：
4.根据权利要求3所述的迭代法引导滤波压制随机噪音方法，其特征在于，在步骤21，利用下列公式，分别计算引导图像I的均值meanI，地震数据p的均值meanp，引导图像I的相关系数corrI，引导图像I与地震数据p的相关系数corrIp；
5.根据权利要求4所述的迭代法引导滤波压制随机噪音方法，其特征在于，在步骤22，利用下列公式，分别计算引导图像I的方差varI，引导图像I与地震数据p的协方差covIp；
6.根据权利要求5所述的迭代法引导滤波压制随机噪音方法，其特征在于，在步骤23，利用下列公式，分别计算线性相关参数a和b；
7.根据权利
8.根据权利要求7所述的迭代法引导滤波压制随机噪音方法，其特征在于，在步骤25，利用下列公式，计算引导滤波结果q；
9.根据权利要求1所述的迭代法引导滤波压制随机噪音方法，其特征在于，在步骤3，以去噪后数据作为输入图像p，重复步骤2，直至迭代次数达到预设值，随机噪音得到很好压制。
10.根据权利要求9所述的迭代法引导滤波压制随机噪音方法，其特征在于，在步骤3，预设值为地震数据自身。
...【技术特征摘要】
1.迭代法引导滤波压制随机噪音方法，其特征在于，该迭代法引导滤波压制随机噪音方法包括：
2.根据权利要求1所述的迭代法引导滤波压制随机噪音方法，其特征在于，在步骤2，应用引导滤波，将地震数据p自身作为引导图像i，以后均以上次去噪后数据作为引导图像i，将地震数据p进行引导滤波，提高资料信噪比。
3.根据权利要求2所述的迭代法引导滤波压制随机噪音方法，其特征在于，步骤2包括：
4.根据权利要求3所述的迭代法引导滤波压制随机噪音方法，其特征在于，在步骤21，利用下列公式，分别计算引导图像i的均值meani，地震数据p的均值meanp，引导图像i的相关系数corri，引导图像i与地震数据p的相关系数corrip；
5.根据权利要求4所述的迭代法引导滤波压制随机噪音方法，其特征在于，在步骤22，利用下列公式，分别计算引导图像i...
【专利技术属性】
技术研发人员：于正军，李继光，赵爱国，滕厚华，刘鹏翔，王胜阁，韩智，
申请(专利权)人：中国石油化工股份有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。