[發明專利]一種三維地質體空間插值方法及系統

有效

• 申請號: 201810048606.0
• 申請日: 2018-01-18
• 公開（公告）號: CN110058298B
• 公開（公告）日: 2020-08-07
• 發明（設計）人: 孫旭東;何江;孫旭;鄒本友;張好林
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油工程技術研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京聿華聯合知識產權代理有限公司 11611
• 代理人: 張文娟

1.一種三維地質體空間插值方法，其特征在于，所述方法包括如下步驟：

數據預處理步驟，對執行三維地質體空間插值所需的數據進行預處理，所述數據包括層位數據、井數據、測井屬性資料、測區數據和地震數據；

三維地質體空間插值步驟，基于預處理后的層位數據、井數據、測井屬性資料、測區數據和地震數據，利用插值算法完成層控和井控、或者結合地震約束的插值約束得到插值后的屬性體；

其中，在利用插值算法完成層控和井控的插值約束時，所述三維地質體空間插值步驟進一步包括如下步驟：

層位網格化步驟，對預處理后的層位數據進行逐層網格化，得到層位網格面；

層位校正步驟，通過校正算法利用預處理后的井數據中的井分層數據對所述層位網格面進行插值運算，得到校正后的構造層位曲面數據；

屬性插值步驟，對預處理后的測井屬性資料與校正后的構造層位曲面數據進行數據提取，并對提取的數據進行屬性插值得到插值后的屬性曲面，此時的屬性曲面是在構造面的基礎上進行的二次運算，任意空間點具有屬性值，屬性值包括巖石力學參數、地應力、壓力預測和可壓性評價屬性；

屬性校正步驟，利用預處理后的測井屬性資料與所述屬性曲面運算對井數據進行校正，得到校正后的層位屬性面。

2.根據權利要求1所述的方法，其特征在于，在利用插值算法完成層控、井控和地震約束的插值約束時，所述三維地質體空間插值步驟還包括如下步驟：

層控體插值步驟，將校正后的層位屬性面與預處理后的測區數據/地震數據進行運算，實現層控體插值，最終完成三維地質體空間插值。

3.根據權利要求2所述的方法，其特征在于，在所述層控體插值步驟中，通過如下步驟完成層控體插值：

將校正后的層位屬性面與預處理后的測區數據進行坐標匹配運算，將層位坐標與測區的線道號對應；

對于每個線道號上的共深度點CDP位置，按照深度從小到大依次對每個CDP進行排序，在排序后對線道號上所有的CDP進行插值，完成該屬性的線性體插值。

4.根據權利要求2所述的方法，其特征在于，在所述層控體插值步驟中，通過如下步驟完成層控體插值：

將校正后的層位屬性面與預處理后的測區數據進行坐標與線道號匹配運算，對于同一線道號的每個共深度點CDP位置，按照深度從小到大依次進行排序；

依據層位的排列順序，依次按照層位控制，將預處理后的地震數據抽取成分段數據，進行波形相干計算；

對波形相干計算結果進行三維屬性體插值。

5.根據權利要求1～4中任一項所述的方法，其特征在于，在層位網格化步驟中，采用如下任一算法進行逐層網格化：反距離加權法、離散光滑約束法、薄板樣條法和雙線性插值。

6.根據權利要求1～4中任一項所述的方法，其特征在于，在屬性插值步驟中，采用如下任一算法進行插值：反距離加權法、克里金和最小曲率法。

7.一種三維地質體空間插值系統，其特征在于，所述系統包括如下模塊：

數據預處理模塊，其對執行三維地質體空間插值所需的數據進行預處理，所述數據包括層位數據、井數據、測井屬性資料、測區數據和地震數據；

三維地質體空間插值模塊，其基于預處理后的層位數據、井數據、測井屬性資料、測區數據和地震數據，利用插值算法完成層控和井控、或者結合地震約束的插值約束得到插值后的屬性體；

其中，所述三維地質體空間插值模塊包括：

層位網格化子模塊，其對預處理后的層位數據進行逐層網格化，得到層位網格面；

層位校正子模塊，其通過校正算法利用預處理后的井數據中的井分層數據對所述層位網格面進行運算，得到校正后的構造層位曲面數據，此時的屬性曲面是在構造面的基礎上進行的二次運算，任意空間點具有屬性值，屬性值包括巖石力學參數、地應力、壓力預測和可壓性評價屬性；

屬性插值子模塊，其對預處理后的測井屬性資料與校正后的構造層位曲面數據進行數據提取，并對提取的數據進行屬性插值得到插值后的屬性曲面；

屬性校正子模塊，其利用預處理后的測井屬性資料與所述屬性曲面運算進行井校，得到校正后的層位屬性面。