【技术实现步骤摘要】
谐波干扰压制方法、装置、存储介质及电子设备
[0001]本申请涉及地震勘探
,特别地涉及一种谐波干扰压制方法、装置、存储介质及电子设备。
技术介绍
[0002]随着可控震源地震勘探的出现,人们很早就认识了谐波,由于滑动扫描出现之前的扫描方式中,谐波干扰对数据品质的影响很小,几乎可以忽略,所以谐波压制技术发展缓慢。为了提高生产效率,出现了滑动扫描等高效采集技术,由于人们进一步追究采集速度,滑动时间越来越小,谐波干扰影响也就越来越大,因此压制谐波干扰的技术不断更新,特别是国外,近些年提出了一系列压制谐波干扰的方法。Li等人(1994)基于线性频率扫描方式,提出了纯相移滤波(PPSM)方法来压制谐波干扰,该方法简单高效直观稳定,特别适用于VSP资料,但不适用于滑动扫描记录。对此,黄建平等人(2012)对其进行改进,提出了一种可以适用于滑动扫描技术的相移滤波法,该方法能够有效的压制本炮干扰,但是对于邻炮干扰的压制效果不是很好。
[0003]因此,谐波干扰是目前可控震源高效采集数据中存在的主要干扰。目前针对谐波压制,一般是利用相关前的母记录和力信号来压制谐波干扰的,但在生产中不适合快速应用。根据谐波的能量异常的特点可以应用中值滤波进行压制,但是中值滤波不能识别谐波干扰的位置,往往需要进行整炮的干扰压制,降低了计算效率和去噪的保真度,并且在去噪过程中可能会伤害的有效波。
技术实现思路
[0004]针对上述问题,本申请提供一种谐波干扰压制方法、装置、存储介质及电子设备,解决了相关技术中由于中值滤波需要 ...
【技术保护点】
【技术特征摘要】
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炮之间所有道的平均均方根振幅值;根据所...
【专利技术属性】
技术研发人员:莫延钢,张媛,姜大建,张丰麒,刘来祥,
申请(专利权)人:中国石油化工股份有限公司石油勘探开发研究院,
类型:发明
国别省市:
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。