一种冰层振源信号到时拾取方法和装置

本发明涉及一种冰层振源信号到时拾取方法和装置，尤其涉及一种基于熵函数分析的到时拾取方法，属于极地声学。

背景技术：

1、信号到时拾取是信号分析处理的基础，其精确程度直接影响了振源定位、介质波速测量的准确性。受冰层波导影响，冰层接收信号呈现明显的非平稳以及频散特性，信号峰值对应的频率、波速波动较大，信号相关法同样很难获取准确到时差。对于非平稳信号，到时拾取通常由具有丰富经验的人进行拾取，但其拾取速度和精度受人为因素影响较大，随着数据量增大，人工拾取已无法满足大量振动信号拾取的实时性要求。

2、目前国内外针对非平稳振动信号到时拾取的方法主要有长短时窗平均能量比法(sta/lta)、赤池信息量准则(aic)算法等。但对于冰层接收信号，其初始到达信号幅值较小，与噪声之间不存在幅值突变，因此基于能量变化的sta/lta方法的拾取精度较差；而aic方法需提前确定信号初始到达的时窗位置和时窗长度，否则到时拾取误差很大。专利申请《一种基于长短时窗和ar模型方差激增效应的p波到时拾取方法》提出使用sta/lta方法初步确定信号到达时窗位置，通过截取信号片段结合aic算法实现精确p波到时拾取。但受极地环境下的复杂噪声干扰影响，初始到达信号信噪比较低，导致sta/lta方法很难通过选取阈值去有效确定信号初步到达时窗，同样降低了到时拾取的准确性。针对信噪比较低的信号，专利申请《一种微震震动波到时拾取方法》以及《一种岩体试件超声波初至拾取方法》中提出的先滤波后拾取到时信息的方法，但由于冰层接收信号与极区噪声信号之间还存在时频域混叠，会导致在滤波时丢失有效初至信号，增加到时拾取误差，因此也不适用。

3、综上所述，现有的到时提取方法很难满足冰层振源信号的精准到时拾取需求。因此，需针对极地冰声信号特点提出一种冰层振源信号到时拾取方法，以解决极地声源定位以及冰层波速测量的需求。

技术实现思路

1、本发明的目的是为了突破现有振源信号到时拾取方法面对冰层振源信号时适用性受限，无法有效提取有效信号的到时问题而提出的一种冰层振源信号到时拾取方法。该方法基于模糊熵对极地冰声信号复杂度进行分析，从而有效区分冰层噪声信号和有效信号，实现极地有效信号的到时拾取。

2、本发明对采集的冰声信号进行分段处理分析，对每段信号依次进行模糊熵分析，截取模糊熵值突增的信号段；针对极地冰层信号，设计频域敏感的aic算法特征函数，并对截取信号进行分析，以准确拾取振源信号到时。

3、本发明提出一种冰层振源信号到时拾取方法，该方法包括以下步骤：

4、s1、对冰面上的接收信号x进行分段，得到m段长度为n的分段信号si,i＝1,2,…,m；

5、s2、计算每个分段信号si的模糊熵值fueni；

6、s3、初至时窗截取

7、对每一分段信号si的模糊熵值fueni做差求得模糊熵变化率vi，通过最大值搜寻得到模糊熵变化率vi的最大值vj，将最大值vj对应的分段信号以及前两段分段信号sj-2,sj-1,sj组合得到信号长度为3n的初至振源信号段sig＝[sj-2,sj-1,sj]；

8、s4、aic函数值计算

9、设计aic特征函数cf：

10、cf(sig(p))＝2sig(p)-(sig(p+1)+sig(p-1))，2＜p＜3n-1

11、计算aic函数值aic：

12、aic(q)＝q*log{var(cf(sig([1,q])))}

13、+(3n-q-1)*log{var(cf(sig([q+1,3n])))}

14、其中q＝2,3,4,…,3n-2，var表示方差；

15、s5、振源信号到时拾取

16、记aic函数值aic最小值对应信号时刻为t，则t为振源信号的初至到达拾取时刻。

17、有利地，s1中，对信号长度为n的接收信号做傅里叶变换，确定有效信号的中心频率f0，结合采样率fs，以10倍中心频率f0的周期长度对信号x进行分段，得到分段信号si。

18、有利地，长度n＝10fs/f0。

19、有利地，s2包括以下步骤：

20、s2.1、相空间重组

21、将每个分段信号si以相空间维数m进行相空间重组得到k个时间序列向量其中k＝1,2,3,…,n-m+1；

22、s2.2、计算最大绝对距离以及重组向量模糊隶属值

23、根据下式计算所有两组时间序列向量和之间的最大绝对距离k,l＝1,2,3，…,n-m+1,k≠l：

24、

25、根据下式计算所有两组时间序列向量和之间的模糊隶属值

26、

27、r为相似度容限；

28、s2.3、模糊度计算

29、针对每个时间序列向量求其模糊隶属值的特定平均值

30、

31、通过下式计算分段信号si在相空间维数m下的模糊度

32、

33、s2.4、计算模糊熵

34、令相空间维数m＝m+1，重复s2.1～s2.3得到分段信号si在相空间维数m+1下的模糊度计算每一分段信号si的模糊熵值：

35、

36、有利地，相空间维数m为根据经验的给定值。

37、有利地，r＝0.1*σ，σ为分段信号si的标准差。

38、有利地，s3中，

39、本发明还提供一种冰层振源信号到时拾取装置，该装置包括检波器以及信号处理装置，通过检波器在冰面上获取接收信号x，通过信号处理装置对接收信号x进行处理和分析，执行如权利要求1-7中任一项所述的冰层振源信号到时拾取方法。

40、有利地，所述信号处理装置包括分段模块、模糊熵值计算模块、初至时窗截取模块、aic值计算模块以及振源信号到时拾取模块，其中

41、所述分段模块对冰面上的接收信号x进行分段，得到m段长度为n的分段信号si,i＝1,2,…,m；

42、所述模糊熵值计算模块计算每个分段信号si的模糊熵值fueni；

43、所述初至时窗截取模块对每一分段信号si的模糊熵值fueni做差求得模糊熵变化率vi，通过最大值搜寻得到模糊熵变化率vi的最大值vj，将最大值vj对应的分段信号以及前两段分段信号sj-2,sj-1,sj组合得到信号长度为3n的初至振源信号段sig＝[sj-2,sj-1,sj]；

44、所述aic值计算模块计算aic特征函数cf：

45、cf(sig(p))＝2sig(p)-(sig(p+1)+sig(p-1))，2＜p＜3n-1

46、计算aic函数值aic：

47、aic(q)＝q*log{var(cf(sig([1,q])))}

48、+(3n-q-1)*log{var(cf(sig([q+1,3n])))}

49、其中q＝2,3,4,…,3n-2，var表示方差；

50、所述振源信号到时拾取模块记aic函数值aic最小值对应信号时刻为t，则t为振源信号的初至到达时刻。

51、有利地，所述模糊熵值计算模块执行以下步骤：

52、s1、相空间重组

53、将每个分段信号si以相空间维数m进行相空间重组得到k个时间序列向量其中k＝1,2,3,…,n-m+1；

54、s2、计算最大绝对距离以及重组向量模糊隶属值

55、根据下式计算所有两组时间序列向量和之间的最大绝对距离k,l＝1,2,3，…,n-m+1,k≠l：

56、

57、根据下式计算所有两组时间序列向量和之间的模糊隶属值

58、

59、r为相似度容限；

60、s3、模糊度计算

61、针对每个时间序列向量求其模糊隶属值的特定平均值

62、

63、通过下式计算分段信号si在相空间维数m下的模糊度

64、

65、s4、计算模糊熵

66、令相空间维数m＝m+1，重复s2.1～s2.3得到分段信号si在相空间维数m+1下的模糊度计算每一分段信号si的模糊熵值：

67、

68、与现有技术相比，本发明的有益效果是：

69、1.本发明基于冰层振源信号频域覆盖范围广的特点，利用模糊熵对信号混乱程度进行评估，在没有丢失有效信号的基础上，更为精准的估计了冰层有效信号的到达临界时窗，能有效提高后续到时拾取以及信号特征分析的准确性。

70、2.本发明针对初至冰层振源信号信噪比低、频域广特点，设计了对频率变化更为敏感、对幅值变化感知较弱的特征函数计算赤池信息量准则，有效提高了冰层振源信号到时拾取的精度。

71、已经讨论的特征、功能和优点可在各种示例中独立实现，或者可以在其他示例中进行组合。