谐波干扰压制方法、装置、存储介质及电子设备与流程

1.本技术涉及地震勘探技术领域，特别地涉及一种谐波干扰压制方法、装置、存储介质及电子设备。


背景技术：

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


技术实现要素：

4.针对上述问题，本技术提供一种谐波干扰压制方法、装置、存储介质及电子设备，解决了相关技术中由于中值滤波需要进行整炮的干扰压制，降低了计算效率和去噪的保真度，并且在去噪过程中可能会伤害的有效波的技术问题。
5.第一方面，本技术提供了一种谐波干扰压制方法，所述方法包括：
6.获取目标segd道头记录的放炮时间顺序及sps数据；
7.将滑动扫描的单炮记录按照所述放炮时间顺序进行排序，得到单炮记录时间序列；
8.根据所述sps数据确认每个单炮记录的炮检点坐标，所述炮检点坐标包括炮点和检波点的坐标信息；
9.根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，其中n为大于或等于1的整数；
10.在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围；
11.根据谐波干扰在时间上的分布规律计算出第n+1炮对第n炮在时间上的影响范围；
12.采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间上的影响范围内，进行谐波干扰压制。
13.优选的，所述根据所述sps数据确认每个单炮记录的炮检点坐标，包括：
14.根据所述sps数据定义观测系统，然后根据sps数据编号与sps观测系统的对应关系，对所述目标segd道头记录进行加载，获得每一个单炮的炮点和检波点坐标信息。
15.优选的，所述根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，包括：
16.根据计算公式计算第n+1炮在第n炮的位置，所述位置包括谐波干扰所在的道号和谐波干扰所在的排列号，计算公式为：
17.谐波干扰所在的道号为：
18.harm_c＝int(x_s
n-x_sn_c1)/r_dis+int(x_s
n+1-x_sn)/r_dis
19.谐波干扰所在的排列号为：
20.harm_rline＝int(y_s
n+1-y_sn_c1)/rline_dis+first_rline_no.
21.其中，x_sn为第n炮的x坐标，x_sn_c1为第n炮第一道的x坐标，r_dis为道距离，x_s
n+1
为第n+1炮的x坐标；y_s
n+1
为第n+1炮的y坐标，y_sn_c1为第n炮第一道的y坐标，rline_dis为接收线距离，first_rline_no为第一条接收线的线号。
22.优选的，所述在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围，包括：
23.根据计算公式计算谐波干扰所在排列中每个道的均方根振幅值，计算公式：
[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
为扫描起始频率，fu为扫描终止频率，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]
第二计算单元，用于根据谐波干扰在时间上的分布规律计算出第n+1炮对第n炮在时间上的影响范围；
[0047]
压制单元，用于采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间上的影响范围内，进行谐波干扰压制。
[0048]
第三方面，一种存储介质，该存储介质存储的计算机程序，可被一个或多个处理器执行，可用来实现如上述第一方面所述的谐波干扰压制方法。
[0049]
第四方面，一种电子设备，包括存储器和处理器，所述存储器上存储有计算机程序，所述存储器和所述处理器之间互相通信连接，该计算机程序被所述处理器执行时，执行如上述第一方面所述的谐波干扰压制方法。
[0050]
本技术提供的一种谐波干扰压制方法、装置、存储介质及电子设备，包括：获取目标segd道头记录的放炮时间顺序及sps数据；将滑动扫描的单炮记录按照所述放炮时间顺序进行排序，得到单炮记录时间序列；根据所述sps数据确认每个单炮记录的炮检点坐标，所述炮检点坐标包括炮点和检波点的坐标信息；根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，其中n为大于或等于1的整数；在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围；根据谐波干扰在时间上的分布规律计算出第n+1炮对第n炮在时间上的影响范围；采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间上的影响范围内，进行谐波干扰压制。通过本方法可以自动的计算出滑动扫描采集中第n+1炮产生的谐波干扰对第n炮的影响范围，并且精准的在该影响范围内用时空变分频中值方法进行干扰压制，可以大大提高谐波压制的速度和效果。
附图说明
[0051]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据
提供的附图获得其他的附图。
[0052]
图1为本技术实施例提供的一种谐波干扰压制方法的流程示意图；
[0053]
图2为本技术实施例提供的谐波干扰在单炮记录图中的分布范围；
[0054]
图3为本技术实施例提供的经过谐波干扰压制后的单炮记录图；
[0055]
图4为本技术实施例提供的谐波干扰压制前后的叠加剖面对比图；
[0056]
图5为本技术实施例提供的一种谐波干扰压制装置的结构示意图；
[0057]
图6为本技术实施例提供的一种电子设备的连接框图。
具体实施方式
[0058]
以下将结合附图及实施例来详细说明本技术的实施方式，借此对本技术如何应用技术手段来解决技术问题，并达到相应技术效果的实现过程能充分理解并据以实施。本技术实施例以及实施例中的各个特征，在不相冲突前提下可以相互结合，所形成的技术方案均在本技术的保护范围之内。
[0059]
由背景技术可知，目前针对谐波压制，一般是利用相关前的母记录和力信号来压制谐波干扰的，但在生产中不适合快速应用。根据谐波的能量异常的特点可以应用中值滤波进行压制，但是中值滤波不能识别谐波干扰的位置，往往需要进行整炮的干扰压制，降低了计算效率和去噪的保真度，并且在去噪过程中可能会伤害的有效波。
[0060]
有鉴于此，本技术提供一种谐波干扰压制方法、装置、存储介质及电子设备，解决了相关技术中由于中值滤波需要进行整炮的干扰压制，降低了计算效率和去噪的保真度，并且在去噪过程中可能会伤害的有效波的技术问题。
[0061]
实施例一
[0062]
图1为本技术实施例提供的一种谐波干扰压制方法的流程示意图，如图1所示，本方法包括：
[0063]
s101、获取目标segd道头记录的放炮时间顺序及sps数据；
[0064]
s102、将滑动扫描的单炮记录按照所述放炮时间顺序进行排序，得到单炮记录时间序列；
[0065]
s103、根据所述sps数据确认每个单炮记录的炮检点坐标，所述炮检点坐标包括炮点和检波点的坐标信息；
[0066]
s104、根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，其中n为大于或等于1的整数；
[0067]
s105、在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围；
[0068]
s106、根据谐波干扰在时间上的分布规律计算出第n+1炮对第n炮在时间上的影响范围；
[0069]
s107、采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间上的影响范围内，进行谐波干扰压制。
[0070]
需要说明的是，本发明的目的是提供了一种自适应可控震源滑动扫描谐波压制方法，可以根据sps的信息和谐波干扰的能量分布特点及时间上的分布规律，自动计算滑动扫描谐波干扰的影响范围，并在谐波干扰的影响范围内应用中值滤波对谐波干扰进行压制。实现自动对谐波干扰的范围进行识别，提高谐波压制的针对性和效果。
[0071]
优选的，所述根据所述sps数据确认每个单炮记录的炮检点坐标，包括：
[0072]
根据所述sps数据定义观测系统，然后根据sps数据编号与sps观测系统的对应关系，对所述目标segd道头记录进行加载，获得每一个单炮的炮点和检波点坐标信息。
[0073]
需要说明的是，具体可以根据野外提供的sps文件定义观测系统，然后根据文件号与sps观测系统的对应关系，对道头进行加载获得每一个单炮的炮点和检波点坐标信息。
[0074]
优选的，所述根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，包括：
[0075]
根据计算公式计算第n+1炮在第n炮的位置，所述位置包括谐波干扰所在的道号和谐波干扰所在的排列号，计算公式为：
[0076]
谐波干扰所在的道号为：
[0077]
harm_c＝int(x_s
n-x_sn_c1)/r_dis+int(x_s
n+1-x_sn)/r_dis
[0078]
谐波干扰所在的排列号为：
[0079]
harm_rline＝int(y_s
n+1-y_sn_c1)/rline_dis+first_rline_no.
[0080]
其中，x_sn为第n炮的x坐标，x_sn_c1为第n炮第一道的x坐标，r_dis为道距离，x_s
n+1
为第n+1炮的x坐标；y_s
n+1
为第n+1炮的y坐标，y_sn_c1为第n炮第一道的y坐标，rline_dis为接收线距离，first_rline_no为第一条接收线的线号。
[0081]
需要说明的是，本技术通过计算谐波干扰所在的位置的道号和排列号来确定所述谐波干扰所在的位置。
[0082]
优选的，所述在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围，包括：
[0083]
根据计算公式计算谐波干扰所在排列中每个道的均方根振幅值，计算公式：
[0084][0085]
其中，f(x)为时间域函数，n为样点个数；
[0086]
计算所述第n+1炮和所述第n炮之间所有道的平均均方根振幅值；
[0087]
根据所述每个道的均方根振幅值和所述平均均方根振幅值，计算得到所述谐波干扰分布的空间范围。
[0088]
优选的，所述根据所述每个道的均方根振幅值和所述平均均方根振幅值，计算得到所述谐波干扰分布的空间范围，包括：
[0089]
分别将所述每个道的均方根振幅值除以所述平均均方根振幅值，得到多个半径值r；
[0090]
获取所述多个半径值r中最大的两个r值对应的两个道号，两个道号之间的范围即为所述谐波干扰分布的空间范围。
[0091]
需要说明的是，采用两个最大的r值对应的道号之间的范围，可以最大范围的包含谐波干扰范围，避免遗漏干扰谐波，导致最终的压制效果不理想。
[0092]
进一步需要说明的是，本技术压制谐波干扰从空间和时间两个方面来进行，这样的方式可以更加精准的进行谐波烦扰压制，更有针对性，尽可能的在去噪过程中不伤害有效波。这一部分为确定空间范围，下一部分为确定时间范围。
[0093]
优选的，所述根据谐波干扰在时间上的分布规律计算出第n+1炮对第n炮在时间上的影响范围，包括：
[0094]
根据如下两个公式计算谐波干扰的时间影响范围：
[0095][0096][0097]
其中：t1为谐波干扰的起始时间，t2为谐波干扰的终止时间，t为扫描长度，k为谐波阶次，f
l
为扫描起始频率，fu为扫描终止频率，w为扫描频宽。
[0098]
需要说明的是，本实施例通过确定谐波干扰的起始时间和终止时间来确定谐波干扰的时间范围，在最终的谐波干扰压制时，只在确定的时间范围和上一部分中确定的空间范围进行谐波干扰压制，更加精准，更有针对性，尽可能的在去噪过程中不伤害有效波。
[0099]
优选的，所述采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间上的影响范围内，进行谐波干扰压制，包括：
[0100]
根据如下中值滤波计算公式进行谐波干扰压制：
[0101]
g(x,y)＝med{f(x-k,y-l),(k,l∈w)}，
[0102]
其中，f(x,y)，g(x,y)分别为原始图像和压制后图像，k为谐波阶次，w为扫描频宽。
[0103]
需要说明的是，通过成像可以看出进行谐波干扰压制前后的区别，如图2所示，为谐波干扰在单炮记录图中的分布范围，可以看到图中左下角有一片黑色的区域，即为谐波干扰，而如图3所示，为经过谐波干扰压制后的单炮记录图，对比图2可以看到，黑色的区域消失了，这就表示成功去除了谐波干扰。
[0104]
综上所述，本技术实施例提供了一种谐波干扰压制方法，包括：获取目标segd道头记录的放炮时间顺序及sps数据；将滑动扫描的单炮记录按照所述放炮时间顺序进行排序，得到单炮记录时间序列；根据所述sps数据确认每个单炮记录的炮检点坐标，所述炮检点坐标包括炮点和检波点的坐标信息；根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，其中n为大于或等于1的整数；在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围；根据谐波干扰在时间上的分布规律计算出第n+1炮对第n炮在时间上的影响范围；采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间上的影响范围内，进行谐波干扰压制。通过本方法可以自动的计算出滑动扫描采集中第n+1炮产生的谐波干扰对第n炮的影响范围，并且精准的在该影响范围内用时空变分频中值方法进行干扰压制，可以大大提高谐波压制的速度和效果。
[0105]
实施例二
[0106]
基于上述本发明实施例公开的谐波干扰压制方法，本技术实施例二通过具体实例应用该谐波干扰压制方法。
[0107]
获取目标segd道头记录的放炮时间顺序及sps数据；
[0108]
将滑动扫描的单炮记录按照所述放炮时间顺序进行排序，得到单炮记录时间序列；
[0109]
根据所述sps数据确认每个单炮记录的炮检点坐标，所述炮检点坐标包括炮点和检波点的坐标信息；
[0110]
根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，其中n为大于或等于1的整数；
[0111]
其中，在本实施例中，谐波干扰所在的位置为道号为526，排列号为1327。
[0112]
在第1327线来统计计算每道的振幅值，然后根据求解半径的计算方法，计算得到的半径为15道的空间范围；
[0113]
然后计算得到的2次谐波干扰的时间范围为2000-18000；
[0114]
在所述半径为15道的空间范围和2000-18000的时间范围内进行谐波干扰压制。
[0115]
如图4所示，为谐波干扰压制前后的叠加剖面对比图。从图4中可以看到压制后的叠加剖面相比于压制前的叠加剖面，更加清晰，最终可以得到更准确的有效值。
[0116]
实施例三
[0117]
基于上述本发明实施例公开的谐波干扰压制方法，图5具体公开了应用该谐波干扰压制方法的谐波干扰压制装置。
[0118]
如图5所示，本发明实施例公开了一种谐波干扰压制装置，该装置包括：
[0119]
获取单元501，用于获取目标segd道头记录的放炮时间顺序及sps数据；
[0120]
排序单元502，用于将滑动扫描的单炮记录按照所述放炮时间顺序进行排序，得到单炮记录时间序列；
[0121]
确认单元503，用于根据所述sps数据确认每个单炮记录的炮检点坐标，所述炮检点坐标包括炮点和检波点的坐标信息；
[0122]
第一计算单元504，用于根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，其中n为大于或等于1的整数；
[0123]
分析单元505，用于在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围；
[0124]
第二计算单元506，用于根据谐波干扰在时间上的分布规律计算出第n+1炮对第n炮在时间上的影响范围；
[0125]
压制单元507，用于采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间上的影响范围内，进行谐波干扰压制。
[0126]
需要说明的是，本发明的目的是提供了一种自适应可控震源滑动扫描谐波压制方法，可以根据sps的信息和谐波干扰的能量分布特点及时间上的分布规律，自动计算滑动扫描谐波干扰的影响范围，并在谐波干扰的影响范围内应用中值滤波对谐波干扰进行压制。实现自动对谐波干扰的范围进行识别，提高谐波压制的针对性和效果。
[0127]
优选的，所述根据所述sps数据确认每个单炮记录的炮检点坐标，包括：
[0128]
根据所述sps数据定义观测系统，然后根据sps数据编号与sps观测系统的对应关系，对所述目标segd道头记录进行加载，获得每一个单炮的炮点和检波点坐标信息。
[0129]
需要说明的是，具体可以根据野外提供的sps文件定义观测系统，然后根据文件号与sps观测系统的对应关系，对道头进行加载获得每一个单炮的炮点和检波点坐标信息。
[0130]
优选的，所述根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，包括：
[0131]
根据计算公式计算第n+1炮在第n炮的位置，所述位置包括谐波干扰所在的道号和谐波干扰所在的排列号，计算公式为：
[0132]
谐波干扰所在的道号为：
[0133]
harm_c＝int(x_s
n-x_sn_c1)/r_dis+int(x_s
n+1-x_sn)/r_dis
[0134]
谐波干扰所在的排列号为：
[0135]
harm_rline＝int(y_s
n+1-y_sn_c1)/rline_dis+first_rline_no.
[0136]
其中，x_sn为第n炮的x坐标，x_sn_c1为第n炮第一道的x坐标，r_dis为道距离，x_s
n+1
为第n+1炮的x坐标；y_s
n+1
为第n+1炮的y坐标，y_sn_c1为第n炮第一道的y坐标，rline_dis为接收线距离，first_rline_no为第一条接收线的线号。
[0137]
需要说明的是，本技术通过计算谐波干扰所在的位置的道号和排列号来确定所述谐波干扰所在的位置。
[0138]
优选的，所述在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围，包括：
[0139]
根据计算公式计算谐波干扰所在排列中每个道的均方根振幅值，计算公式：
[0140][0141]
其中，f(x)为时间域函数，n为样点个数；
[0142]
计算所述第n+1炮和所述第n炮之间所有道的平均均方根振幅值；
[0143]
根据所述每个道的均方根振幅值和所述平均均方根振幅值，计算得到所述谐波干扰分布的空间范围。
[0144]
优选的，所述根据所述每个道的均方根振幅值和所述平均均方根振幅值，计算得到所述谐波干扰分布的空间范围，包括：
[0145]
分别将所述每个道的均方根振幅值除以所述平均均方根振幅值，得到多个半径值r；
[0146]
获取所述多个半径值r中最大的两个r值对应的两个道号，两个道号之间的范围即为所述谐波干扰分布的空间范围。
[0147]
需要说明的是，采用两个最大的r值对应的道号之间的范围，可以最大范围的包含谐波干扰范围，避免遗漏干扰谐波，导致最终的压制效果不理想。
[0148]
进一步需要说明的是，本技术压制谐波干扰从空间和时间两个方面来进行，这样的方式可以更加精准的进行谐波烦扰压制，更有针对性，尽可能的在去噪过程中不伤害有效波。这一部分为确定空间范围，下一部分为确定时间范围。
[0149]
优选的，所述根据谐波干扰在时间上的分布规律计算出第n+1炮对第n炮在时间上的影响范围，包括：
[0150]
根据如下两个公式计算谐波干扰的时间影响范围：
[0151][0152][0153]
其中：t1为谐波干扰的起始时间，t2为谐波干扰的终止时间，t为扫描长度，k为谐波阶次，f
l
为扫描起始频率，fu为扫描终止频率，w为扫描频宽。
[0154]
需要说明的是，本实施例通过确定谐波干扰的起始时间和终止时间来确定谐波干扰的时间范围，在最终的谐波干扰压制时，只在确定的时间范围和上一部分中确定的空间范围进行谐波干扰压制，更加精准，更有针对性，尽可能的在去噪过程中不伤害有效波。
[0155]
优选的，所述采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间上的影响范围内，进行谐波干扰压制，包括：
[0156]
根据如下中值滤波计算公式进行谐波干扰压制：
[0157]
g(x,y)＝med{f(x-k,y-l),(k,l∈w)}，
[0158]
其中，f(x,y)，g(x,y)分别为原始图像和压制后图像，k为谐波阶次，w为扫描频宽。
[0159]
需要说明的是，通过成像可以看出进行谐波干扰压制前后的区别，如图2所示，为谐波干扰在单炮记录图中的分布范围，可以看到图中左下角有一片黑色的区域，即为谐波干扰，而如图3所示，为经过谐波干扰压制后的单炮记录图，对比图2可以看到，黑色的区域消失了，这就表示成功去除了谐波干扰。
[0160]
综上所述，本技术实施例提供了一种谐波干扰压制装置，包括：获取目标segd道头记录的放炮时间顺序及sps数据；将滑动扫描的单炮记录按照所述放炮时间顺序进行排序，得到单炮记录时间序列；根据所述sps数据确认每个单炮记录的炮检点坐标，所述炮检点坐标包括炮点和检波点的坐标信息；根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，其中n为大于或等于1的整数；在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围；根据谐波干扰在时间上的分布规律计算出第n+1炮对第n炮在时间上的影响范围；采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间上的影响范围内，进行谐波干扰压制。通过本方法可以自动的计算出滑动扫描采集中第n+1炮产生的谐波干扰对第n炮的影响范围，并且精准的在该影响范围内用时空变分频中值方法进行干扰压制，可以大大提高谐波压制的速度和效果。
[0161]
实施例四
[0162]
本实施例还提供一种计算机可读存储介质，如闪存、硬盘、多媒体卡、卡型存储器(例如，sd或dx存储器等)、随机访问存储器(ram)、静态随机访问存储器(sram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、可编程只读存储器(prom)、磁性存储器、磁盘、光盘、服务器、app应用商城等等，其上存储有计算机程序，该计算机程序被处理器执行时可以实现如实施例一的谐波干扰压制方法，步骤如下：
[0163]
获取目标segd道头记录的放炮时间顺序及sps数据；
[0164]
将滑动扫描的单炮记录按照所述放炮时间顺序进行排序，得到单炮记录时间序列；
[0165]
根据所述sps数据确认每个单炮记录的炮检点坐标，所述炮检点坐标包括炮点和检波点的坐标信息；
[0166]
根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，其中n为大于或等于1的整数；
[0167]
在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围；
[0168]
根据谐波干扰在时间上的分布规律计算出第n+1炮对第n炮在时间上的影响范围；
[0169]
采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间上的影响范围内，进行谐波干扰压制。
[0170]
综上所述，本技术实施例提供了一种谐波干扰压制方法，包括：获取目标segd道头记录的放炮时间顺序及sps数据；将滑动扫描的单炮记录按照所述放炮时间顺序进行排序，得到单炮记录时间序列；根据所述sps数据确认每个单炮记录的炮检点坐标，所述炮检点坐标包括炮点和检波点的坐标信息；根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，其中n为大于或等于1的整数；在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围；根据谐波干扰在时间上的分布规律计算出第n+1炮对第n炮在时间上的影响范围；采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间上的影响范围内，进行谐波干扰压制。通过本方法可以自动的计算出滑动扫描采集中第n+1炮产生的谐波干扰对第n炮的影响范围，并且精准的在该影响范围内用时空变分频中值方法进行干扰压制，可以大大提高谐波压制的速度和效果。
[0171]
实施例五
[0172]
图6为本技术实施例提供的一种电子设备600的连接框图，如图6所示，该电子设备600可以包括：处理器601，存储器602，多媒体组件603，输入/输出(i/o)接口604，以及通信组件605。
[0173]
其中，处理器601用于执行如实施例一中的谐波干扰压制方法中的全部或部分步骤。存储器602用于存储各种类型的数据，这些数据例如可以包括电子设备中的任何应用程序或方法的指令，以及应用程序相关的数据。
[0174]
处理器601可以是专用集成电路(application specific integrated circuit，简称asic)、数字信号处理器(digital signal processor，简称dsp)、数字信号处理设备(digital signal processing device，简称dspd)、可编程逻辑器件(programmable logic device，简称pld)、现场可编程门阵列(field programmable gate array，简称fpga)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述实施例一中的谐波干扰压制方法。
[0175]
存储器602可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，例如静态随机存取存储器(static random access memory，简称sram)，电可擦除可编程只读存储器(electrically erasable programmable read-only memory，简称eeprom)，可擦除可编程只读存储器(erasable programmable read-only memory，简称eprom)，可编程只读存储器(programmable read-only memory，简称prom)，只读存储器(read-only memory，简称rom)，磁存储器，快闪存储器，磁盘或光盘。
[0176]
多媒体组件603可以包括屏幕和音频组件，该屏幕可以是触摸屏，音频组件用于输出和/或输入音频信号。例如，音频组件可以包括一个麦克风，麦克风用于接收外部音频信号。所接收的音频信号可以被进一步存储在存储器或通过通信组件发送。音频组件还包括至少一个扬声器，用于输出音频信号。
[0177]
i/o接口604为处理器601和其他接口模块之间提供接口，上述其他接口模块可以是键盘，鼠标，按钮等。这些按钮可以是虚拟按钮或者实体按钮。
[0178]
通信组件605用于该电子设备500与其他设备之间进行有线或无线通信。无线通信，例如wi-fi，蓝牙，近场通信(near field communication，简称nfc)，2g、3g或4g，或它们中的一种或几种的组合，因此相应的该通信组件605可以包括：wi-fi模块，蓝牙模块，nfc模块。
[0179]
综上，本技术提供的一种谐波干扰压制方法、装置、存储介质及电子设备，该方法包括：获取目标segd道头记录的放炮时间顺序及sps数据；将滑动扫描的单炮记录按照所述放炮时间顺序进行排序，得到单炮记录时间序列；根据所述sps数据确认每个单炮记录的炮检点坐标，所述炮检点坐标包括炮点和检波点的坐标信息；根据所述每个单炮的炮检点坐标和所述单炮记录时间序列，计算第n+1炮在第n炮的位置，其中n为大于或等于1的整数；在所述第n+1炮在第n炮的位置分析得到谐波干扰分布的空间范围；根据谐波干扰在时间上的分布规律计算出第n+1炮对第n炮在时间上的影响范围；采用中值滤波方法在所述谐波干扰分布的空间范围内和所述时间上的影响范围内，进行谐波干扰压制。通过本方法可以自动的计算出滑动扫描采集中第n+1炮产生的谐波干扰对第n炮的影响范围，并且精准的在该影响范围内用时空变分频中值方法进行干扰压制，可以大大提高谐波压制的速度和效果。
[0180]
在本技术实施例所提供的几个实施例中，应该理解到，所揭露的方法，也可以通过其它的方式实现。以上所描述的方法实施例仅仅是示意性的。
[0181]
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括要素的过程、方法、物品或者设备中还存在另外的相同要素。
[0182]
虽然本技术所揭露的实施方式如上，但上述的内容只是为了便于理解本技术而采用的实施方式，并非用以限定本技术。任何本技术所属技术领域内的技术人员，在不脱离本技术所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本技术的专利保护范围，仍须以所附的权利要求书所界定的范围为准。