[发明专利]一种次声信号幅值季节性偏差校正方法

权利要求书

1.一种次声信号幅值季节性偏差校正方法，其特征在于：其包括以下步骤：

S1、对目标地区历史爆炸事件进行收集，得到各次声台站记录的次声信号；

S2、将收集的次声信号按照台站进行分类，记录事件发生日期、台站名称、台站至事件的距离、台站至事件的方位角、爆炸事件的当量、监测的次声信号幅值；

S3、对次声信号幅值进行当量归一化处理，并取对数，得到归一化幅值amp_adjust：

S4、将事件发生日期转换为日历天Jdate，即全年第x天；以日历天为横坐标，归一化幅值amp_adjust为纵坐标，对所指定台站记录的所有次声信号进行分析；

S5、根据归一化幅值-日历天所生成的图中所展示的各台站信号幅值变化趋势，将次声台站分为先上升后降低型和先降低后上升型；

S6、按照S5得到的两种类型的台站的地理分布划定分割线，将两种类型的台站分割开来，对分割线两侧台站分别进行幅值校正；对需要幅值校正的台站，将归一化幅值分别进行二阶、三阶、四阶拟合，确定最能反映幅值变化规律的多阶拟合函数；

S7、对于同一类型的台站，将S6得到的某一台站多阶拟合函数在其它台站进行适用性进行分析，确认最能反映幅值变化规律的拟合函数阶数；

S8、将S7确定的某一台站的拟合函数在同一类台站进行适用性分析，最终确定适合全部或大部分台站使用的拟合函数，从而获得普遍适用的通用校正函数；

S9、对于某一个具体事件，将其所发生的日历天代入S8的通用校正函数中，即可获得当日次声信号的校正系数Amp2，得到经校正后的次声信号幅值amp_new为：

amp_new＝amp_adjust-Amp2

S10、根据当量计算公式确定经信号幅值季节性偏差校正后的平流层顶次声信号幅值-当量-传输距离关系。

2.根据权利要求1所述的一种次声信号幅值季节性偏差校正方法，其特征在于：

amp为次声信号幅值，w为爆炸事件的当量。

3.根据权利要求1所述的一种次声信号幅值季节性偏差校正方法，其特征在于：根据台站的实际情况对S8获得的通用校正函数进行微调，包括将日期x左右移动和幅值上下移动，得到微调之后的校正函数；

Amp2＝a1*(x+d)²+a2*(x+d)+a3+B

当d为正数时表示x左移d天，d为负数时表示x右移d天，d为整数；B为正数时幅值向上移动，B为负数时幅值向下移动；a1、a2、a3分别为拟合系数。

4.根据权利要求3所述的一种次声信号幅值季节性偏差校正方法，其特征在于：x的取值范围为Amp2进行次声幅值校正的时间区间。

5.根据权利要求3所述的一种次声信号幅值季节性偏差校正方法，其特征在于：S10中，经信号幅值季节性偏差校正后的平流层顶次声信号幅值-当量-传输距离关系式为：

log(Amp)＝-1.2*log(D)+0.5*log(w)+3.6159+Amp2

其中，Amp为次声信号幅值，Amp2为次声信号幅值校正系数；D为信号传输距离。

6.根据权利要求1所述的一种次声信号幅值季节性偏差校正方法，其特征在于：S6中将归一化幅值分别进行二阶、三阶、四阶拟合，得到反映幅值变化规律的多阶拟合函数分别为：

Amp2＝a1*x²+a2*x+a3

Amp2＝a1*x³+a2*x²+a3*x+a4

Amp2＝a1*x⁴+a2*x³+a3*x²+a4*x+a5

a1、a2、a3、a4、a5分别为拟合系数。

7.根据权利要求1所述的一种次声信号幅值季节性偏差校正方法，其特征在于：所述次声信号幅值季节性偏差是指相同能量水平次声事件由于发生的季节性差异所体现出的幅值偏差。