[發明專利]基于廣義S變換的譜相關系數的Q值估計方法

在審

• 申請號: 201610289700.6
• 申請日: 2016-05-04
• 公開（公告）號: CN107346034A
• 公開（公告）日: 2017-11-14
• 發明（設計）人: 余青露;居興國;李進;鄒少峰;肖盈;劉思思;高艷霞;祝媛媛
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京思創畢升專利事務所11218
• 代理人: 孫向民,廉莉莉

1.一種基于廣義S變換的譜相關系數的Q值估計方法，包括以下步驟：

步驟1：輸入地震數據；

步驟2：對所述地震數據進行廣義S變換，得到變換后的地震數據；

步驟3：獲得變換后的地震數據的估計振幅譜和實際振幅譜；

步驟4：基于所述估計振幅譜和實際振幅譜，構建譜相關系數；

步驟5：求解所述譜相關系數最小時的Q值。

2.根據權利要求1所述的基于廣義S變換的譜相關系數的Q值估計方法，其中所述地震數據為U(ω,z)，其中ω是地震波的角頻率，z是地震波的傳播距離。

3.根據權利要求1所述的基于廣義S變換的譜相關系數的Q值估計方法，其中根據以下公式(8)得到變換后的地震數據S_N(τ,f)：

SN(τ,f)=|f|p2πλ∫-∞∞exp(-πfτQ+λπ22fp-2Q2)exp[-fp2λ(t-τ+λπfp-1Q)2]dt=exp(-πfτQ+λπ22fp-2Q2)---(8)]]>

其中，f表示頻率，t表示地震波從地面開始傳播到待分析地層的旅行時間，τ表示所述窗函數w(t-τ,f)的時間位置，λ表示用于調節所述窗函數的時間延續長度的參數，p表示用于調節所述窗函數的衰減趨勢的參數，其中所述窗函數w(t-τ,f)表示為以下公式(1)：

w(t-τ,f)=|f|p2πλexp[-fp(t-τ)22λ],λ>0,p>0---(1).]]>