技術(shù)特征：

1.一種自動(dòng)衡量地震數(shù)據(jù)諧波噪聲強(qiáng)弱的振幅譜比值方法，其特征在于，包括以下步驟：

步驟01：讀取原始地震記錄信號(hào)的單道數(shù)據(jù)；

步驟02：對(duì)步驟01讀取的單道數(shù)據(jù)做快速傅里葉變換，從時(shí)域轉(zhuǎn)換到頻率域；

步驟03：確定高低頻頻譜分界，計(jì)算高頻能量占總頻帶能量的比值；

重復(fù)步驟01-03直到所有道數(shù)據(jù)處理完成。

2.如權(quán)利要求1所述一種自動(dòng)衡量地震數(shù)據(jù)諧波噪聲強(qiáng)弱的振幅譜比值方法，其特征在于，步驟02中對(duì)地震資料單道數(shù)據(jù)做快速傅里葉變換，具體包括：

假設(shè)x(t)表示待分析信號(hào)，X(ω)為x(t)的頻譜，使用快速傅里葉變換求取頻譜：

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

假設(shè)待分析信號(hào)為離散信號(hào)，長度為N，則離散的頻譜長度也為N；由于快速傅里葉變換得到的頻譜具有對(duì)稱性，因此只考慮前N/2個(gè)離散點(diǎn)頻譜；

假設(shè)第i點(diǎn)的頻率值為ω_i，則有下式：

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

其中，dt為采樣間隔，則為采樣頻率；采樣頻率記為f_N，采樣頻率一半記為f_N/2。

3.如權(quán)利要求1所述一種自動(dòng)衡量地震數(shù)據(jù)諧波噪聲強(qiáng)弱的振幅譜比值方法，其特征在于，步驟03中高低頻頻譜分界為40Hz。

4.如權(quán)利要求1所述一種自動(dòng)衡量地震數(shù)據(jù)諧波噪聲強(qiáng)弱的振幅譜比值方法，其特征在于，步驟03中計(jì)算高頻能量占總頻帶能量的比值，包括：

選擇40Hz作為低頻與高頻分界，接下來計(jì)算高頻段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頻率點(diǎn)的位置；

計(jì)算總頻段0-f_N/2的能量，記為f_all：

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

計(jì)算高頻段能量與總頻段能量的比值，記為α：

α＝f_high/f_all。

5.如權(quán)利要求4所述一種自動(dòng)衡量地震數(shù)據(jù)諧波噪聲強(qiáng)弱的振幅譜比值方法，其特征在于，步驟03中在計(jì)算高頻段能量時(shí)將100Hz以上的頻率點(diǎn)舍棄掉；此時(shí)，f_high、f_all的改進(jìn)形式為：

$\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},$

其中，