[發(fā)明專利]基于Lp范數(shù)約束的相位識別方法及系統(tǒng)

在審

• 申請?zhí)枺?/td>: 201910983984.2
• 申請日: 2019-10-16
• 公開（公告）號: CN112666603A
• 公開（公告）日: 2021-04-16
• 發(fā)明（設(shè)計）人: 凡睿;劉勇;謝紅斌;王威;潘蓓;蘇克露;季春輝
• 申請（專利權(quán)）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勘探分公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京思創(chuàng)大成知識產(chǎn)權(quán)代理有限公司 11614
• 代理人: 高爽

1.一種基于Lp范數(shù)約束的相位識別方法，其特征在于，包括：

根據(jù)地震記錄褶積模型，構(gòu)建非平穩(wěn)地震褶積模型；

根據(jù)所述非平穩(wěn)地震褶積模型與Lp范數(shù)，建立目標函數(shù)；

通過迭代軟閾值法求解所述目標函數(shù)的復(fù)數(shù)解；

根據(jù)所述目標函數(shù)的復(fù)數(shù)解，識別時頻相位譜。

2.根據(jù)權(quán)利要求1所述的基于Lp范數(shù)約束的相位識別方法，其中，所述非平穩(wěn)地震褶積模型為：

其中，s(t)為地震數(shù)據(jù)，t為時間，N為子波的個數(shù)，f_i為第i個子波的主頻，[W(t,f_i)]為子波褶積矩陣，[r(t,f_i)]為偽反射系數(shù)矩陣，G為子波褶積矩陣[W(t,f_i)]，m為偽反射系數(shù)矩陣[r(t,f_i)]。

3.根據(jù)權(quán)利要求1所述的基于Lp范數(shù)約束的相位識別方法，其中，所述Lp范數(shù)為：

其中，x_i∈(x₁,x₂,x₃…x_n)。

4.根據(jù)權(quán)利要求1所述的基于Lp范數(shù)約束的相位識別方法，其中，所述目標函數(shù)為：

其中，表示L2范數(shù)的平方，||·||_p表示Lp范數(shù)，μ為正則化參數(shù)。

5.根據(jù)權(quán)利要求1所述的基于Lp范數(shù)約束的相位識別方法，其中，所述迭代軟閾值法的最終迭代方程為：

其中，η的值為大于G^HG的最大特征值，m_k為第k次迭代的計算結(jié)果，soft(·)為軟閾值迭代法。

6.一種基于Lp范數(shù)約束的相位識別系統(tǒng)，其特征在于，該系統(tǒng)包括：

存儲器，存儲有計算機可執(zhí)行指令；

處理器，所述處理器運行所述存儲器中的計算機可執(zhí)行指令，執(zhí)行以下步驟：

根據(jù)地震記錄褶積模型，構(gòu)建非平穩(wěn)地震褶積模型；

根據(jù)所述非平穩(wěn)地震褶積模型與Lp范數(shù)，建立目標函數(shù)；

通過迭代軟閾值法求解所述目標函數(shù)的復(fù)數(shù)解；

根據(jù)所述目標函數(shù)的復(fù)數(shù)解，識別時頻相位譜。

7.根據(jù)權(quán)利要求6所述的基于Lp范數(shù)約束的相位識別系統(tǒng)，其中，所述非平穩(wěn)地震褶積模型為：

其中，s(t)為地震數(shù)據(jù)，t為時間，N為子波的個數(shù)，f_i為第i個子波的主頻，[W(t,f_i)]為子波褶積矩陣，[r(t,f_i)]為偽反射系數(shù)矩陣，G為子波褶積矩陣[W(t,f_i)]，m為偽反射系數(shù)矩陣[r(t,f_i)]。

8.根據(jù)權(quán)利要求6所述的基于Lp范數(shù)約束的相位識別系統(tǒng)，其中，所述Lp范數(shù)為：

其中，x_i∈(x₁,x₂,x₃…x_n)。

9.根據(jù)權(quán)利要求6所述的基于Lp范數(shù)約束的相位識別系統(tǒng)，其中，所述目標函數(shù)為：

其中，表示L2范數(shù)的平方，||·||_p表示Lp范數(shù)，μ為正則化參數(shù)。

10.根據(jù)權(quán)利要求6所述的基于Lp范數(shù)約束的相位識別系統(tǒng)，其中，所述迭代軟閾值法的最終迭代方程為：

其中，η的值為大于G^HG的最大特征值，m_k為第k次迭代的計算結(jié)果，soft(·)為軟閾值迭代法。