【技术实现步骤摘要】
本专利技术涉及地震勘探领域，特别涉及一种地震资料被动源干扰自适应压制方法。
技术介绍
1、海量数据存储及gpu计算能力的提高，推动了采集技术的进步，地震勘探已经从常规的窄方位采集、发展到宽方位、宽频带、高密度采集(“两宽一高”)，并正在向高密度全方位采集方向发展。“两宽一高”地震技术优势明显，既是解决当前地震勘探所面临的复杂问题的系统对策，也是今后地震勘探的发展趋势。高效的采集方式，实现了地震信号的空间充分采样，全面记录有效信号及噪声的波场细节，为后续的研究提供了丰富信息。但由于单点接收、小道距、高道数、高覆盖、宽频大动态范围的采集方式，使得原始记录中干扰波异常发育，这也对后续的去噪效果提出了挑战。地震记录中接收到的干扰波主要分为两类：一类是与震源有关的面波、折射波、异常振幅等；另一类是与激发震源无关的各种地面震动源产生的干扰波，称为“被动源干扰”。
2、被动源干扰的发育程度与干扰源的数目、特性及其与地震测线的相对位置关系有关，其能量强、频带宽，严重影响整体地震资料的品质，因此，被动源干扰的压制是后续去噪处理的重点和难点，必须寻找有效的精细保幅噪声压制方法。
3、近年来，国内外有很多学者开始对被动源干扰波的压制方法进行关注和研究。burch(burch d,kappius r,tessman j.successful noise attenuation using pointreceivers in the presence of aliased noise.2004，seg technical prog
技术实现思路
1、本专利技术所要解决的技术问题是克服
技术介绍
中存在的现有处理方法会因滤波窗口选取不当造成低频有效信号损失或噪声压制不完全有残余的问题，而提供一种地震资料被动源干扰自适应压制方法，该地震资料被动源干扰自适应压制方法，为海量采集、被动源干扰波严重的地震资料提供一种高效、保幅的压制方法，有效提高地震资料的信噪比和保真度。
2、本专利技术解决其问题可通过如下技术方案来达到：该地震资料被动源干扰自适应压制方法，包括以下步骤：
3、s1.在时间域从原始数据中分选出受被动源干扰影响的地震数据；
4、s2.将分选出的受被动源干扰影响的地震数据转换到频率域；
5、s3.将s2中得到的频率域地震数据，在频率域分析识别干扰波的频带及视速度范围，并拟合算子；
6、s4.将得到的频率域地震数据与拟合算子进行褶积，并将结果进行反傅立叶变换到时间域，得到拟合后的时间域线性干扰数据；
7、s5.将步骤s4得到的时间域线性干扰数据在时间域与分选出的受被动源干扰影响的地震数据进行自适应匹配相减，从而得到被动源干扰压制后的地震道；
8、s6.将被动源干扰波压制后的地震道合并到原始地震数据当中，完成数据重排列处理；
9、s7.将原始地震数据与重排列处理后的地震道进行质控对比分析，确定被动源干扰是否被压制彻底。
10、进一步的，所述步骤s1在时间域从原始数据中分选出受被动源干扰影响的地震数据的方法，包括：
11、输入数据集应为时间域炮集数据，进行数据处理前，需在原始单炮记录中，按照初至波前的能量，统计出能量超过0的数据，自动分选出能量超过0的数据；
12、将自动分选出能量超过0的这些数据单独分选出来作为受被动源干扰影响的地震道；
13、进一步的，所述输入数据为置道头处理后未经任何噪声压制处理的原始数据。
14、进一步的，所述步骤s2将时间域分选出来受被动源干扰影响的地震数据变换到频率域的方法，采用以下的公式：
15、
16、式中：d(t,x)为时间域地震道；d(f,x)为频率域地震数据；
17、f为频率；i为采样点数。
18、进一步的，将s2中得到的频率域地震数据，在频率域进行分析识别干扰波的频带及视速度范围的方法，包括：
19、将s2中得到的频率域的数据，进行干扰波的频带及视速度分析，在频率域分别统计被动源干扰波的最大、最小频率及视速度的范围；
20、分析时，应尽量选择能量较强，干扰波特征明显的数据进行统计；
21、同时要尽可能的多分析一些干扰波，以便能覆盖到全区被动源干扰的特征，保证预测的算子更合理、准确。
22、进一步的，s3计算拟合算子的方法：
23、通过求解下列方程计算拟合算子：
24、
25、式中：“*”代表复共轭；a(m)为拟合算子；r(m)＝ei2πkmδx，m＝1,2…,p-1；p为算子长度；k为波数，即频率与确定的视速度比值；δ2为白噪因子；波数陷波宽度由算子长度p控制，p越大，陷波宽度越窄。
26、进一步的，所述步骤s4将得到的频率域地震数据与拟合算子通过公式(3)进行褶积处理：
27、
28、式中：d(f,x-m)为频率域地震数据，a(m)为拟合算子，d'(f,x)为褶积处理结果，f为频率，其中x＝1,2,…,n。
29、进一步的，所述s4将褶积结果进行反傅立叶变换到时间域，得到拟合后的线性干扰的方法，包括：
30、将褶积结果d'(f,x)进行反傅立叶变换采用的变换公式(4)，将数据由频率域变换到时间域，得到拟合后的线性干扰；
31、
32、式中：d'(t,x)为时间域线性干扰数据；d'(f,x)为频率域线性干扰数据；i＝1,2,…n；w为圆频率；t为时间。
33、进一步的，所述步骤s7进行质控对比分析的方法，包括：
34、将压制被动源干扰后的数据与未进行被动源干扰压制的原始地震数据进行差剖面对比，分析被动源干扰波是否被压制彻底；
35、同时对减出的差剖面分析是否携带有效信号。
36、进一步的，所述方法适用于宽方位、宽频带、高密度地震资料被动源干扰自适应压制。
37、本专利技术与上述
技术介绍
相比较可具有如下有益效果：
38、本专利技术地震资料被动源干扰自适应压制方法，本文档来自技高网...
【技术保护点】
1.一种地震资料被动源干扰自适应压制方法，其特征在于：包括以下步骤:
2.根据权利要求1所述的地震资料被动源干扰自适应压制方法，其特征在于：步骤S1在时间域从原始数据中分选出受被动源干扰影响的地震数据的方法，包括：
3.根据权利要求1所述的地震资料被动源干扰自适应压制方法，其特征在于：所述输入数据为置道头处理后未经任何噪声压制处理的原始数据。
4.根据权利要求1所述的地震资料被动源干扰自适应压制方法，其特征在于：S2将时间域分选出来受被动源干扰影响的地震数据变换到频率域的方法，采用以下的公式：
5.根据权利要求1所述的地震资料被动源干扰自适应压制方法，其特征在于：将S2中得到的频率域地震数据，在频率域进行分析识别干扰波的频带及视速度范围的方法，包括：
6.根据权利要求1所述的地震资料被动源干扰自适应压制方法，其特征在于：S3计算拟合算子的方法：
7.根据权利要求1所述的地震资料被动源干扰自适应压制方法，其特征在于：所述步骤S4将得到的频率域地震数据与拟合算子通过公式(3)进行褶积处理：
8.根据权
9.根据权利要求1所述的地震资料被动源干扰自适应压制方法，其特征在于：所述步骤S7进行质控对比分析的方法，包括：
10.根据权利要求1所述的地震资料被动源干扰自适应压制方法，其特征在于：所述方法适用于宽方位、宽频带、高密度地震资料被动源干扰自适应压制。
...【技术特征摘要】
1.一种地震资料被动源干扰自适应压制方法，其特征在于：包括以下步骤:
2.根据权利要求1所述的地震资料被动源干扰自适应压制方法，其特征在于：步骤s1在时间域从原始数据中分选出受被动源干扰影响的地震数据的方法，包括：
3.根据权利要求1所述的地震资料被动源干扰自适应压制方法，其特征在于：所述输入数据为置道头处理后未经任何噪声压制处理的原始数据。
4.根据权利要求1所述的地震资料被动源干扰自适应压制方法，其特征在于：s2将时间域分选出来受被动源干扰影响的地震数据变换到频率域的方法，采用以下的公式：
5.根据权利要求1所述的地震资料被动源干扰自适应压制方法，其特征在于：将s2中得到的频率域地震数据，在频率域进行分析识别干扰波的频带及视速度范围的方法，包括...
【专利技术属性】
技术研发人员：赵忠华，初海红，包燚，王成，晏丰，张芝铭，
申请(专利权)人：大庆油田有限责任公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。