[發明專利]基于網格自適應加密與粗化技術的數值模擬方法和裝置

有效

• 申請號: 201810419159.5
• 申請日: 2018-05-04
• 公開（公告）號: CN108710734B
• 公開（公告）日: 2022-12-20
• 發明（設計）人: 龔斌;李俊超;劉玄;蘭正凱
• 申請（專利權）人: 特雷西能源科技（杭州）有限公司
• 主分類號: G06F30/20
• 分類號: G06F30/20;G06F111/10;G06F113/08
• 代理公司: 北京中和立達知識產權代理有限公司 11756
• 代理人: 祝妍

1.一種基于網格自適應加密與粗化技術的數值模擬方法，其特征在于，包括如下步驟：

S100：加載原始網格模型，所述原始網格模型包括網格孔隙度數據、網格深度數據、網格體積數據和網格間傳導率數據；

S200：構建多級嵌套網格系統，所述多級嵌套網格系統包括順序設置的原始網格模型和將所述原始網格模型進行集合劃分，依次生成的若干級嵌套的粗化網格模型，其中后一級網格模型是對前一級網格模型的粗化；

S300：基于所述多級嵌套網格系統的第一級網格模型啟動數值模擬運算；

S400：在數值模擬運算過程中，當前時間步狀態更新完畢之后，將相鄰網格單元的狀態信息差值與基準參數進行比較，根據比較結果將該網格單元加密為前一級網格模型的網格單元，粗化為后一級網格模型的網格單元或不變，得到新的多級嵌套網格系統；

S500：計算新的多級嵌套網格系統的網格孔隙度數據、網格深度數據、網格體積數據和網格間傳導率數據；

S600：在下一時間步，對新的多級嵌套網格系統執行S300～S500，直至數值模擬運算結束；

其中，所述S400包括：

S410：獲取當前時間步所有網格單元的狀態信息，所述狀態信息包括飽和度、壓力和/或組分濃度；

S420：計算每個網格單元狀態信息與其所有相鄰網格單元的差值，選取絕對值最大的差值作為狀態信息差值；

S430：將所述狀態信息差值與基準參數進行比較，若所有狀態信息差值均大于基準參數，則該網格單元代表油藏劇烈變化的復雜流動區域，包括水驅前沿或近井地帶，將該網格單元加密為前一級網格模型的網格單元，若所有狀態信息差值均小于基準參數，則將該網格單元粗化為后一級網格模型的網格單元；其中，只允許網格單元在相鄰級別狀態上進行調整，使得網格單元變化平滑；

所述S430包括：

S431：將所述狀態信息差值與基準參數進行比較，若所有狀態信息差值均大于基準參數，則將該網格單元標記為加密，若所有狀態信息差值均小于基準參數，則將該網格單元標記為粗化；

S432：若網格單元被標記為加密，且該網格單元不處在第一級網格模型，則將該網格單元加密為前一級網格模型的網格單元；若網格單元被標記為粗化，且該網格單元不處在最后一級網格模型，且其對應的后一級網格單元包含的所有網格單元均標記為粗化，則將該網格單元粗化為后一級網格模型的網格單元；

所述S420中，相鄰的網格單元通過傳導率聯通表信息查找；所述狀態信息包括飽和度和壓力，所述狀態信息差值包括飽和度差ΔS_i和壓力差ΔP_i，其中：

S_i和P_i分別為網格單元i的飽和度和壓力，S_j和P_i分別為與網格單元i所有相鄰網格單元j的飽和度和壓力，t_ij為傳導率，t_ij≠0代表網格單元i和網格單元j相鄰；

所述基準參數包括ΔS_min、ΔS_max、ΔP_min和ΔP_max，其中：若ΔS_i＞ΔS_max且ΔP_i＞ΔP_max，則將網格單元i標記為加密；若ΔS_i＜ΔS_min且ΔP_i＜ΔP_min，則將網格單元i標記為粗化，其他情況，則網格單元i不做標記。

2.根據權利要求1所述的基于網格自適應加密與粗化技術的數值模擬方法，其特征在于：

所述S300進一步為：調用數值模擬器，基于所述多級嵌套網格系統的第一級網格模型啟動數值模擬運算；

網格模型還包括PVT屬性數據、相滲數據、初始化數據和生產控制數據；

所述S500進一步為：運用基于流動的網格粗化方法，計算新的多級嵌套網格系統的網格孔隙度數據、網格深度數據、網格體積數據、網格間傳導率數據和網格單元的狀態信息。