谐波干扰压制方法、装置、存储介质及电子设备制造方法及图纸

本申请涉及地震勘探技术领域，具体涉及一种谐波干扰压制方法、装置、存储介质及电子设备，包括：获取目标SEGD道头记录的放炮时间顺序及SPS数据；根据所述目标SEGD道头记录的放炮时间顺序及SPS数据，得到谐波干扰分布的空间范围内和时间影响范围；采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间影响范围内，进行谐波干扰压制。通过本方法可以自动的计算出滑动扫描采集中第n+1炮产生的谐波干扰对第n炮的影响范围，并且精准的在该影响范围内用时空变分频中值方法进行干扰压制，可以大大提高谐波压制的速度和效果。可以大大提高谐波压制的速度和效果。可以大大提高谐波压制的速度和效果。

【技术实现步骤摘要】
谐波干扰压制方法、装置、存储介质及电子设备


[0001]本申请涉及地震勘探
，特别地涉及一种谐波干扰压制方法、装置、存储介质及电子设备。

技术介绍

[0002]随着可控震源地震勘探的出现，人们很早就认识了谐波，由于滑动扫描出现之前的扫描方式中，谐波干扰对数据品质的影响很小，几乎可以忽略，所以谐波压制技术发展缓慢。为了提高生产效率，出现了滑动扫描等高效采集技术，由于人们进一步追究采集速度，滑动时间越来越小，谐波干扰影响也就越来越大，因此压制谐波干扰的技术不断更新，特别是国外，近些年提出了一系列压制谐波干扰的方法。Li等人(1994)基于线性频率扫描方式，提出了纯相移滤波(PPSM)方法来压制谐波干扰，该方法简单高效直观稳定，特别适用于VSP资料，但不适用于滑动扫描记录。对此，黄建平等人(2012)对其进行改进，提出了一种可以适用于滑动扫描技术的相移滤波法，该方法能够有效的压制本炮干扰，但是对于邻炮干扰的压制效果不是很好。
[0003]因此，谐波干扰是目前可控震源高效采集数据中存在的主要干扰。目前针对谐波压制，一般是利用相关前的母记录和力信号来压制谐波干扰的，但在生产中不适合快速应用。根据谐波的能量异常的特点可以应用中值滤波进行压制，但是中值滤波不能识别谐波干扰的位置，往往需要进行整炮的干扰压制，降低了计算效率和去噪的保真度，并且在去噪过程中可能会伤害的有效波。

技术实现思路

[0004]针对上述问题，本申请提供一种谐波干扰压制方法、装置、存储介质及电子设备，解决了相关技术中由于中值滤波需要进行整炮的干扰压制，降低了计算效率和去噪的保真度，并且在去噪过程中可能会伤害的有效波的技术问题。
[0005]第一方面，本申请提供了一种谐波干扰压制方法，所述方法包括：
[0006]获取目标SEGD道头记录的放炮时间顺序及SPS数据；
[0007]将滑动扫描的单炮记录按照所述放炮时间顺序进行排序，得到单炮记录时间序列；
[0008]根据所述SPS数据确认每个单炮记录的炮检点坐标，所述炮检点坐标包括炮点和检波点的坐标信息；
[0009]根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，其中n为大于或等于1的整数；
[0010]在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围；
[0011]根据谐波干扰在时间上的分布规律计算出第n+1炮对第n炮在时间上的影响范围；
[0012]采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间上的影响范围内，进行谐波干扰压制。
[0013]优选的，所述根据所述SPS数据确认每个单炮记录的炮检点坐标，包括：
[0014]根据所述SPS数据定义观测系统，然后根据SPS数据编号与SPS观测系统的对应关系，对所述目标SEGD道头记录进行加载，获得每一个单炮的炮点和检波点坐标信息。
[0015]优选的，所述根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，包括：
[0016]根据计算公式计算第n+1炮在第n炮的位置，所述位置包括谐波干扰所在的道号和谐波干扰所在的排列号，计算公式为：
[0017]谐波干扰所在的道号为：
[0018]Harm_c＝int(x_s
n
‑
x_s
n
_c1)/r_dis+int(x_s
n+1
‑
x_s
n
)/r_dis
[0019]谐波干扰所在的排列号为：
[0020]Harm_rline＝int(y_s
n+1
‑
y_s
n
_c1)/rline_dis+first_rline_no.
[0021]其中，x_s
n
为第n炮的x坐标，x_s
n
_c1为第n炮第一道的x坐标，r_dis为道距离，x_s
n+1
为第n+1炮的x坐标；y_s
n+1
为第n+1炮的y坐标，y_s
n
_c1为第n炮第一道的y坐标，rline_dis为接收线距离，first_rline_no为第一条接收线的线号。
[0022]优选的，所述在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围，包括：
[0023]根据计算公式计算谐波干扰所在排列中每个道的均方根振幅值，计算公式：
[0024][0025]其中，f(x)为时间域函数，N为样点个数；
[0026]计算所述第n+1炮和所述第n炮之间所有道的平均均方根振幅值；
[0027]根据所述每个道的均方根振幅值和所述平均均方根振幅值，计算得到所述谐波干扰分布的空间范围。
[0028]优选的，所述根据所述每个道的均方根振幅值和所述平均均方根振幅值，计算得到所述谐波干扰分布的空间范围，包括：
[0029]分别将所述每个道的均方根振幅值除以所述平均均方根振幅值，得到多个半径值R；
[0030]获取所述多个半径值R中最大的两个R值对应的两个道号，两个道号之间的范围即为所述谐波干扰分布的空间范围。
[0031]优选的，所述根据谐波干扰在时间上的分布规律计算出第n+1炮对第n炮在时间上的影响范围，包括：
[0032]根据如下两个公式计算谐波干扰的时间影响范围：
[0033][0034][0035]其中：T1为谐波干扰的起始时间，T2为谐波干扰的终止时间，T为扫描长度，k为谐波阶次，f
l
为扫描起始频率，f
u
为扫描终止频率，W为扫描频宽。
[0036]优选的，所述采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间上的影响范围内，进行谐波干扰压制，包括：
[0037]根据如下中值滤波计算公式进行谐波干扰压制：
[0038]g(x,y)＝med{f(x
‑
k,y
‑
l),(k,l∈W)}，
[0039]其中，f(x,y)，g(x,y)分别为原始图像和压制后图像，k为谐波阶次，W为扫描频宽。
[0040]第二方面，一种谐波干扰压制装置，所述方法包括：
[0041]获取单元，用于获取目标SEGD道头记录的放炮时间顺序及SPS数据；
[0042]排序单元，用于将滑动扫描的单炮记录按照所述放炮时间顺序进行排序，得到单炮记录时间序列；
[0043]确认单元，用于根据所述SPS数据确认每个单炮记录的炮检点坐标，所述炮检点坐标包括炮点和检波点的坐标信息；
[0044]第一计算单元，用于根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，其中n为大于或等于1的整数；
[0045]分析单元，用于在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围；
[0046]第二计算单元，用于根据本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种谐波干扰压制方法，其特征在于，所述方法包括：获取目标SEGD道头记录的放炮时间顺序及SPS数据；将滑动扫描的单炮记录按照所述放炮时间顺序进行排序，得到单炮记录时间序列；根据所述SPS数据确认每个单炮记录的炮检点坐标，所述炮检点坐标包括炮点和检波点的坐标信息；根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，其中n为大于或等于1的整数；在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围；根据谐波干扰在时间上的分布规律计算出第n+1炮对第n炮在时间上的影响范围；采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间上的影响范围内，进行谐波干扰压制。2.根据权利要求1所述的方法，其特征在于，所述根据所述SPS数据确认每个单炮记录的炮检点坐标，包括：根据所述SPS数据定义观测系统，然后根据SPS数据编号与SPS观测系统的对应关系，对所述目标SEGD道头记录进行加载，获得每一个单炮的炮点和检波点坐标信息。3.根据权利要求1所述的方法，其特征在于，所述根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，包括：根据计算公式计算第n+1炮在第n炮的位置，所述位置包括谐波干扰所在的道号和谐波干扰所在的排列号，计算公式为：谐波干扰所在的道号为：Harm_c＝int(x_s
n
‑
x_s
n
_c1)/r_dis+int(x_s
n+1
‑
x_s
n
)/r_dis谐波干扰所在的排列号为：Harm_rline＝int(y_s
n+1
‑
y_s
n
_c1)/rline_dis+first_rline_no.其中，x_s
n
为第n炮的x坐标，x_s
n
_c1为第n炮第一道的x坐标，r_dis为道距离，x_s
n+1
为第n+1炮的x坐标；y_s
n+1
为第n+1炮的y坐标，y_s
n
_c1为第n炮第一道的y坐标，rline_dis为接收线距离，first_rline_no为第一条接收线的线号。4.根据权利要求3所述的方法，其特征在于，所述在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围，包括：根据计算公式计算谐波干扰所在排列中每个道的均方根振幅值，计算公式：其中，f(x)为时间域函数，N为样点个数；计算所述第n+1炮和所述第n炮之间所有道的平均均方根振幅值；根据所...

【专利技术属性】
技术研发人员：莫延钢，张媛，姜大建，张丰麒，刘来祥，
申请(专利权)人：中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：