[發明專利]基于時頻聯合域反演的楊氏模量和泊松比預測方法及裝置

在審

• 申請號: 201811108113.8
• 申請日: 2018-09-21
• 公開（公告）號: CN110941009A
• 公開（公告）日: 2020-03-31
• 發明（設計）人: 李京南;馬靈偉;朱童
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京聿宏知識產權代理有限公司 11372
• 代理人: 吳大建;張杰

1.一種基于時頻聯合域反演的楊氏模量和泊松比預測方法，包括：

基于貝葉斯定理，構建基于反演模型參數的后驗概率分布、先驗概率分布和似然函數之間的關系，得到初始關系式；

設置所述反演模型參數符合柯西分布，并根據符合柯西分布的反演模型參數獲取基于時頻聯合域反演的似然函數，根據所述符合柯西分布的反演模型參數、基于時頻聯合域反演的似然函數以及所述初始關系式得到基于所述時頻聯合域反演的后驗概率分布函數；

根據所述后驗概率分布函數構建基于所述反演模型參數的目標泛函；

根據所述反演模型參數的約束條件和所述目標泛函得到基于時頻聯合域的目標函數；

對基于時頻聯合域的目標函數進行求解，得到反演模型參數結果。

2.根據權利要求1所述的方法，其特征在于，基于時頻聯合域的目標函數進行求解，得到反演模型參數結果步驟包括：

獲取待預測工區的原始疊前地震數據和測井數據；

根據所述原始疊前地震數據和測井數據計算基于時頻聯合域的目標函數中的各參數的取值；

根據所述各參數的取值對基于時頻聯合域的目標函數進行求解，得到反演模型參數結果。

3.根據權利要求2所述的方法，其特征在于，根據所述原始疊前地震數據計算基于時頻聯合域的目標函數中的各參數的取值步驟包括：

對待預測工區的原始疊前地震數據進行預處理，并對預處理后的數據進行分角度疊加，得到角度疊加道集d_t；對所述角度疊加道集d_t進行傅里葉變換，得到頻率域地震數據d_f。

4.根據權利要求3所述的方法，其特征在于，對待預測工區的原始疊前地震數據進行預處理包括：

對待預測工區的所述原始疊前地震數據依次進行靜校正處理、球面擴散補償處理、吸收衰減補償處理、去噪處理、反褶積處理以及偏移處理。

5.根據權利要求4所述的方法，其特征在于，根據待預測工區的所述測井數據計算基于時頻聯合域的目標函數中的各參數的取值步驟包括：

根據待預測工區的測井數據分別計算楊氏模量、泊松比以及密度參數，并對所述楊氏模量、泊松比和密度參數進行低通濾波處理，根據低通濾波處理后的所述楊氏模量、泊松比和密度參數構建初始反演模型參數m⁰；

通過對所述待預測工區的測井數據進行井震標定獲取地震子波；根據所述地震子波構建地震子波矩陣，并根據所述地震子波矩陣計算時間域正演算子G_t；計算所述地震子波的頻譜，并根據所述地震子波的頻譜計算頻率域正演算子G_f。

6.根據權利要求5所述的方法，其特征在于，根據所述各參數的取值對基于時頻聯合域的目標函數進行求解，得到反演模型參數結果步驟包括：

基于所述角度疊加道集d_t、頻率域地震數據d_f、初始反演模型參數m⁰、時間域正演算子G_t以及頻率域正演算子G_f，根據迭代重加權最小二乘算法對所述時頻聯合域的目標函數進行迭代求解；當所述迭代次數達到預設最大迭代次數時得到的所述反演模型參數作為反演模型參數結果。

7.根據權利要求6所述的方法，其特征在于，所述預設最大迭代次數為4或5。