[發(fā)明專利]一種利用自適應(yīng)分析時窗提取地震屬性的方法、存儲介質(zhì)

有效

• 申請?zhí)枺?/td>: 201910343618.0
• 申請日: 2019-04-26
• 公開（公告）號: CN111856565B
• 公開（公告）日: 2022-04-22
• 發(fā)明（設(shè)計）人: 王猛;王小品;魏三妹;鄭四連
• 申請（專利權(quán)）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術(shù)研究院
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京聿宏知識產(chǎn)權(quán)代理有限公司 11372
• 代理人: 吳大建;張杰

1.一種利用自適應(yīng)分析時窗提取地震屬性的方法，其特征在于，包括：

步驟S100，獲取目的層的地震信號的平均瞬時頻率；

步驟S200，根據(jù)平均瞬時頻率求取分析時窗的長度以及確定分析時窗的頂端和底面的原始層位面；

步驟S300，對分析時窗的頂端和底面的原始層位面分別進(jìn)行平滑處理，獲得頂端和底面的優(yōu)化層位面；其中，所述步驟S300包括以下步驟：

步驟S310，對于分析時窗的頂端/底面的原始層位面中的任意一點(x₀,y₀,t₀)，按照下式計算表征該點的曲面能量的六個分量：

基于上述六個分量形成對稱的二階張量，記為：

令：以建立該點的特征方程組，

其中，λ為求解所述特征方程組獲得的特征值，λ為多解；

取上述特征值中最大的值λ_max作為該點的特征值；

步驟S320，當(dāng)獲得分析時窗的頂端/底面的原始層位面中所有點的特征值后，建立如下趨勢擬合目標(biāo)函數(shù)F：

F＝∑λ_max+∑(H-H′)

其中，H為分析時窗的頂端/底面的原始層位面，H′為頂端/底面的優(yōu)化層位面，

步驟S330，令趨勢擬合目標(biāo)函數(shù)F的值達(dá)到全局最小，由此確定頂端和底面的優(yōu)化層位面H′；

步驟S400，保持分析時窗的長度不變，結(jié)合頂端和底面的優(yōu)化層位面構(gòu)建新的分析時窗；

步驟S500，利用新的分析時窗提取目的層的地震屬性，以進(jìn)行地震屬性分析。

2.根據(jù)權(quán)利要求1所述的利用自適應(yīng)分析時窗提取地震屬性的方法，其特征在于，所述步驟S100包括以下步驟：

步驟S110，將目的層的地震信號轉(zhuǎn)換為離散地震信號；

步驟S120，利用小波變換對目的層的離散地震信號進(jìn)行變換；

步驟S130，基于下式求取變換后的目的層的離散地震信號的平均瞬時頻率：

其中，F(xiàn)T(f,t)為變換后的目的層的離散地震信號，f為瞬時頻率，t為瞬時時間點，f_m為平均瞬時頻率。

3.根據(jù)權(quán)利要求2所述的利用自適應(yīng)分析時窗提取地震屬性的方法，其特征在于，所述步驟S200中，根據(jù)下式確定分析時窗的長度：

其中，n為根據(jù)分析時窗的長度需求而選取的常數(shù)值，T為周期，L為分析時窗的長度。

4.根據(jù)權(quán)利要求3所述的利用自適應(yīng)分析時窗提取地震屬性的方法，其特征在于，n取值范圍為0.1至2。

5.根據(jù)權(quán)利要求3所述的利用自適應(yīng)分析時窗提取地震屬性的方法，其特征在于，所述步驟S200中，分析時窗的頂端和底面的原始層位面分別為：

其中，H_top和H_bot分別為分析時窗的頂端和底面的原始層位面，H₁為目的層的層位面。

6.一種計算機(jī)可讀存儲介質(zhì)，其中存儲有系統(tǒng)程序，所述系統(tǒng)程序在被處理器執(zhí)行時實現(xiàn)如權(quán)利要求1至5中任意一項所述的利用自適應(yīng)分析時窗提取地震屬性的方法。