[发明专利]一种基于模糊域的地震信号分数域最优阶计算方法

摘要

一种基于模糊域的地震信号分数域最优阶计算方法，属于地震信号处理领域，旨在利用模糊域快速确定地震信号分数域的最优阶，简化最优阶计算过程。首先读入单道地震信号，并对输入信号做模糊变换；然后计算信号的时宽和带宽，求取信号基于时频带宽积最小的窗函数，并对窗函数做模糊变换；利用窗函数的模糊函数对信号的模糊变换做低通滤波；计算出地震信号自项的模糊函数的水平和竖直两个方向的方向梯度；对方向梯度矩阵进行奇异值分解，计算地震信号自项模糊函数的主方向；最后通过主方向计算出地震信号的分数域最优阶。

主权项

1.一种基于模糊域的地震信号分数域最优阶计算方法，其特征在于，所述方法包括以下步骤：步骤1：输入一单道地震信号s(t)，对信号s(t)做模糊变换，得到信号的模糊函数As；步骤2：计算地震信号s(t)的时宽Ts和带宽Bs，求取地震信号s(t)基于时频带宽积最小的窗函数hTBP(t)，并对窗函数hTBP(t)做模糊变换，得到窗函数hTBP(t)的模糊函数Ag；步骤3：利用窗函数的模糊函数Ag对信号的模糊函数As做低通滤波，得到信号自项的模糊函数A0；步骤4：计算信号自项的模糊函数A0的方向梯度矩阵G；步骤5：对方向梯度矩阵G进行奇异值分解，利用奇异值分解结果计算信号自项的模糊函数A0的主方向θ；具体理论如下：G＝USVT  (1)式中，U是N×2的矩阵，V是2×2的矩阵，S是2×2对角的奇异值矩阵，矩阵V的第一行向量为V1＝[V1，1，V1，2]信号自项的模糊函数A0的主方向θ的计算公式为：步骤6：利用信号自项的模糊函数A0的主方向θ计算信号s(t)的分数域最优阶popt，