[發明專利]疊前彈性參數反演方法、裝置、計算機設備和存儲介質

在審

• 申請號: 202011055756.8
• 申請日: 2020-09-30
• 公開（公告）號: CN114428301A
• 公開（公告）日: 2022-05-03
• 發明（設計）人: 紀永禎;張渝悅;胡華鋒
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京聿宏知識產權代理有限公司 11372
• 代理人: 吳大建;陳敏

1.一種疊前彈性參數反演方法，其特征在于，包括：

獲取疊前地震數據；

獲取彈性參數的反射系數；

基于Fatti近似方程，根據所述疊前地震數據及所述反射系數，構建疊前地震數據的正演矩陣；

根據所述正演矩陣，得到地震數據的似然函數；

獲取自動相關判別的先驗信息；

基于貝葉斯理論，根據所述似然函數及所述先驗信息，得到所述彈性參數反演的目標函數；

基于最大似然估計法，將似然函數轉化成最大化邊緣似然函數；

對所述最大化邊緣似然函數進行求解，得到估計因子；

根據所述估計因子及所述目標函數，反演得到彈性參數。

2.根據權利要求1所述的疊前彈性參數反演方法，其特征在于，在所述根據所述正演矩陣，得到地震數據的似然函數的步驟之前，還包括：

將所述正演矩陣中加入趨勢約束，得到約束后的所述正演矩陣；

所述根據所述正演矩陣，得到地震數據的似然函數的步驟，包括：

根據約束后的所述正演矩陣，得到所述似然函數。

3.根據權利要求2所述的疊前彈性參數反演方法，其特征在于，約束后的所述正演矩陣表達式為：

其中，Ψ＝F^-1ΛFC∈R^K×K,Λ＝diag[H]，Λ為以漢寧窗函數向量作為對角元素的對角矩陣，H為漢寧窗函數向量，C為積分矩陣，LowIP＝F^-1ΛF(log(IP))∈R^K×1為縱波阻抗的趨勢，LowIS＝F^-1ΛF(log(IS))∈R^K×1橫波阻抗的趨勢，λ₁,λ₂分別為趨勢約束的權重參數，F為離散傅里葉變換矩陣，D為差分矩陣，IP和IS分別表示擬反演的縱波阻抗和橫波阻抗參數的對數，

A(θ)＝diag[A(t₁,θ) A(t₂,θ)…A(t_K,θ)]，

B(θ)＝diag[B(t₁,θ) B(t₂,θ)…B(t_K,θ)]，diag[…]為對角矩陣。

4.根據權利要求3所述的疊前彈性參數反演方法，其特征在于，還包括：將約束后的所述正演矩陣表達式簡化成如下形式：

d_(M*N+2*K)＝G_{(M*N+2*K)×(2K+2)}m_2K+2+n

其中d代表G包含傅里葉變換矩陣、角度相關矩陣、阻抗趨勢提取矩陣和差分矩陣，m代表待估計的縱橫波阻抗參數組成的向量，n代表噪聲向量。

5.根據權利要求1所述的疊前彈性參數反演方法，其特征在于，所述獲取彈性參數的反射系數的步驟，包括：

構建縱橫波阻抗的彈性參數模型；

根據所述彈性參數模型，得到所述反射系數。

6.根據權利要求1所述的疊前彈性參數反演方法，其特征在于，所述目標函數的表達式為：

p(m|d,h,σ²,θ)＝C|Σ|^-1/2exp[-1/2(m-μ)^TΣ^-1(m-μ)]

其中，Σ＝(H+σ^-2G^TG)^-1，μ＝σ^-2ΣG^Td，C為一個常數，μ為待估計縱橫波阻抗彈性參數的后驗概率均值，H＝diag(h₁,h₂,…,h_2K)為一個對角矩陣。

7.根據權利要求1所述的疊前彈性參數反演方法，其特征在于，所述先驗信息的表達式為：

其中，h＝[h₁,h₂,…,h_2K]^T包含2K個獨立的參數。