[發(fā)明專利]一種基于模糊域的地震信號(hào)分?jǐn)?shù)域最優(yōu)階計(jì)算方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201710071409.6
• 申請(qǐng)日: 2017-02-09
• 公開（公告）號(hào): CN106842305B
• 公開（公告）日: 2018-12-18
• 發(fā)明（設(shè)計(jì)）人: 田琳;李禎;王慧玲
• 申請(qǐng)（專利權(quán)）人: 伊犁師范學(xué)院
• 主分類號(hào): G01V1/28
• 分類號(hào): G01V1/28
• 代理公司: 成都弘毅天承知識(shí)產(chǎn)權(quán)代理有限公司 51230
• 代理人: 王正楠

1.一種基于模糊域的地震信號(hào)分?jǐn)?shù)域最優(yōu)階計(jì)算方法，其特征在于，所述方法包括以下步驟：

步驟1：輸入一單道地震信號(hào)s(t)，對(duì)信號(hào)s(t)做模糊變換，得到信號(hào)的模糊函數(shù)A_s；

步驟2：計(jì)算地震信號(hào)s(t)的時(shí)寬T_s和帶寬B_s，求取地震信號(hào)s(t)基于時(shí)頻帶寬積最小的窗函數(shù)h_TBP(t)，并對(duì)窗函數(shù)h_TBP(t)做模糊變換，得到窗函數(shù)h_TBP(t)的模糊函數(shù)A_g；

步驟3：利用窗函數(shù)的模糊函數(shù)A_g對(duì)信號(hào)的模糊函數(shù)A_s做低通濾波，得到信號(hào)自項(xiàng)的模糊函數(shù)A₀；

步驟4：計(jì)算信號(hào)自項(xiàng)的模糊函數(shù)A₀的方向梯度矩陣G；

步驟5：對(duì)方向梯度矩陣G進(jìn)行奇異值分解，利用奇異值分解結(jié)果計(jì)算信號(hào)自項(xiàng)的模糊函數(shù)A₀的主方向θ；

具體理論如下：

G＝USV^T (1)

式中，U是N×2的矩陣，V是2×2的矩陣，S是2×2對(duì)角的奇異值矩陣，矩陣V的第一行向量為V₁＝[V_1，1，V_1，2]

信號(hào)自項(xiàng)的模糊函數(shù)A₀的主方向θ的計(jì)算公式為：

步驟6：利用信號(hào)自項(xiàng)的模糊函數(shù)A₀的主方向θ計(jì)算信號(hào)s(t)的分?jǐn)?shù)域最優(yōu)階p_opt，

2.根據(jù)權(quán)利要求1所述的一種基于模糊域的地震信號(hào)分?jǐn)?shù)域最優(yōu)階計(jì)算方法，其特征在于，所述步驟1中對(duì)信號(hào)s(t)做的模糊變換所采用的變換公式如下：

式中，τ是時(shí)延，θ是頻移，j是虛數(shù)單位，e是自然常數(shù)，t是時(shí)間。

3.根據(jù)權(quán)利要求1所述的一種基于模糊域的地震信號(hào)分?jǐn)?shù)域最優(yōu)階計(jì)算方法，其特征在于，所述步驟2中

對(duì)窗函數(shù)h_TBP(t)做模糊變換，得到窗函數(shù)h_TBP(t)的模糊函數(shù)A_g，變換公式如下：

4.根據(jù)權(quán)利要求1所述的一種基于模糊域的地震信號(hào)分?jǐn)?shù)域最優(yōu)階計(jì)算方法，其特征在于，所述步驟3中的濾波方法，計(jì)算方法如下：

A_o＝A_sA_g (7)。

5.根據(jù)權(quán)利要求1所述的一種基于模糊域的地震信號(hào)分?jǐn)?shù)域最優(yōu)階計(jì)算方法，其特征在于，所述步驟4中信號(hào)自項(xiàng)的模糊函數(shù)A₀是N×N的矩陣信號(hào)，其方向梯度為：

其中g(shù)_x(k)表示矩陣A₀在點(diǎn)(x_k，y_k)水平方向的梯度，g_y(k)表示矩陣A₀在點(diǎn)(x_k，y_k)豎直方向的梯度。