[發明專利]基于旋轉交錯網格的雙相介質地質數據獲取方法和裝置

在審

• 申請號: 201610041813.4
• 申請日: 2016-01-21
• 公開（公告）號: CN105676280A
• 公開（公告）日: 2016-06-15
• 發明（設計）人: 林朋;彭蘇萍;盧勇旭
• 申請（專利權）人: 中國礦業大學（北京）
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京超凡志成知識產權代理事務所(普通合伙) 11371
• 代理人: 王術蘭

1.基于旋轉交錯網格的雙相介質地質數據獲取方法，其特征在于，包括：

根據雙相介質運動方程，確定目標區域的地質數據計算式；

使用旋轉交錯網格對所述地質數據計算式進行離散化處理；

將所述離散化處理后的地質數據計算式在所述目標區域所對應的模型下進行數值模 擬，以確定所述目標區域的優化地質數據，所述優化地質數據包括優化速度值和優化應力 值。

2.根據權利要求1所述的基于旋轉交錯網格的雙相介質地質數據獲取方法，其特征在 于，所述使用旋轉交錯網格對所述地質數據計算式進行離散化處理包括：

按照將速度和密度置于整網格點處，將彈性系數和應力置于半網格點處的方式，計算 所述地質數據計算式中，初步地質數據；所述初步地質數據包括初步速度值和初步應力值。

3.根據權利要求2所述的基于旋轉交錯網格的雙相介質地質數據獲取方法，其特征在 于，所述按照將速度和密度置于整網格點處，將彈性系數和應力置于半網格點處的方式，計 算所述地質數據計算式中，初步速度值和初步應力值包括：

將速度和密度賦值于整網格點處，將彈性系數和應力賦值于半網格點處；

使用所述整網格點沿旋轉交錯網格的對角線的四個應力分量，通過中心差分計算的方 式計算速度對時間的偏導數，即初步速度值；

使用所述半網格點沿旋轉交錯網格的對角線的四個速度分量，通過中心差分計算的方 式計算應力對時間的偏導數，即初步應力值。

4.根據權利要求1所述的基于旋轉交錯網格的雙相介質地質數據獲取方法，其特征在 于，所述將所述離散化處理后的地質數據計算式在所述目標區域所對應的模型下進行數值 模擬包括：

按照隨地震波傳播時間變化的順序，依次循環迭代調整所述離散化處理后的地質數據 計算式，每次迭代后，均保存迭代計算出的速度值和應力值，以確定優化速度值和優化應力 值。

5.根據權利要求1所述的基于旋轉交錯網格的雙相介質地質數據獲取方法，其特征在 于，所述地質數據計算式為：

ρ11∂ui∂t2+ρ12∂Ui∂t2=σij′j-bij(∂Uj∂t-∂uj∂t)]]>

ρ12∂ui∂t2+ρ22∂Ui∂t2=s′i+bij(∂Uj∂t-∂uj∂t);]]>

其中，i和j分別表示水平方向x,y或垂直方向z中的一個，u_j和U_j分別是固相和流相的位 移在j方向的分量，b_ij為流體相對固體骨架運動時的耗散系數，σ_ij'j為作用在固相應力分量 在j方向上的偏導，s為作用在流體單元側面上的應力，ρ₁₁和ρ₂₂分別表示介質單位體積內固 相和流相部分的有效質量，ρ₁₂為視質量，即流相相對固相運動時的質量耦合系數。