一种自动衡量地震数据谐波噪声强弱的振幅谱比值方法与流程

技术特征：

1.一种自动衡量地震数据谐波噪声强弱的振幅谱比值方法，其特征在于，包括以下步骤：

步骤01：读取原始地震记录信号的单道数据；

步骤02：对步骤01读取的单道数据做快速傅里叶变换，从时域转换到频率域；

步骤03：确定高低频频谱分界，计算高频能量占总频带能量的比值；

重复步骤01-03直到所有道数据处理完成。

2.如权利要求1所述一种自动衡量地震数据谐波噪声强弱的振幅谱比值方法，其特征在于，步骤02中对地震资料单道数据做快速傅里叶变换，具体包括：

假设x(t)表示待分析信号，X(ω)为x(t)的频谱，使用快速傅里叶变换求取频谱：

X(ω)＝FFT(x(t))；

假设待分析信号为离散信号，长度为N，则离散的频谱长度也为N；由于快速傅里叶变换得到的频谱具有对称性，因此只考虑前N/2个离散点频谱；

假设第i点的频率值为ω_i，则有下式：

${\mathrm{\ω}}_i=\frac{1}{2\×dt}\frac{i}{N/2},$

其中，dt为采样间隔，则为采样频率；采样频率记为f_N，采样频率一半记为f_N/2。

3.如权利要求1所述一种自动衡量地震数据谐波噪声强弱的振幅谱比值方法，其特征在于，步骤03中高低频频谱分界为40Hz。

4.如权利要求1所述一种自动衡量地震数据谐波噪声强弱的振幅谱比值方法，其特征在于，步骤03中计算高频能量占总频带能量的比值，包括：

选择40Hz作为低频与高频分界，接下来计算高频段40Hz-f_N/2的能量，记为f_high：

$\begin{array}{c}{f}_{high}=\underset{i=i_{40}}{\overset{N/2}{\Σ}}ads\left(X\right({\mathrm{\ω}}_i\left)\right)\\ i_{40}=\frac{40}{f_{N/2}}\×\frac{N}{2}\end{array},$

其中，i₄₀是40Hz频率点的位置；

计算总频段0-f_N/2的能量，记为f_all：

$f_{all}=\underset{i=0}{\overset{N/2}{\Σ}}abs\left(X\right({\mathrm{\ω}}_i\left)\right);$

计算高频段能量与总频段能量的比值，记为α：

α＝f_high/f_all。

5.如权利要求4所述一种自动衡量地震数据谐波噪声强弱的振幅谱比值方法，其特征在于，步骤03中在计算高频段能量时将100Hz以上的频率点舍弃掉；此时，f_high、f_all的改进形式为：

$\begin{array}{c}{f}_{high}=\underset{i=i_{40}}{\overset{i_{i00}}{\Σ}}abs\left(X\right({\mathrm{\ω}}_i\left)\right)\\ f_{all}=\underset{i=0}{\overset{i_{i00}}{\Σ}}abs\left(X\right({\mathrm{\ω}}_i\left)\right)\end{array},$

其中，