技術(shù)特征：

1.一種黑油油藏模擬方法，其特征在于，包括：

基于多個(gè)組分滲流參數(shù)建立描述黑油油藏流體滲流規(guī)律的組分模型；

根據(jù)基于黑油模型的黑油模擬所需數(shù)據(jù)對(duì)所述組分模型進(jìn)行數(shù)值離散，得到基于所述黑油油藏的地質(zhì)模型離散網(wǎng)格的組分?jǐn)?shù)值模型；

利用所述黑油模擬所需數(shù)據(jù)中的流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表；

根據(jù)所述組分相平衡參數(shù)表和所述組分滲流參數(shù)的估計(jì)值進(jìn)行閃蒸計(jì)算，得到所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)；

根據(jù)所述流體相態(tài)和相應(yīng)物性參數(shù)調(diào)整所述組分?jǐn)?shù)值模型，以簡(jiǎn)化組分?jǐn)?shù)值模型；

利用所述組分滲流參數(shù)的估計(jì)值和調(diào)整后的所述組分?jǐn)?shù)值模型，計(jì)算得到油藏壓力變化值及各流體相飽和度變化值；

判斷所述的油藏壓力變化值及各流體相飽和度變化值是否收斂，若是，根據(jù)所述油藏壓力變化值和相應(yīng)的所述估計(jì)值求和得到所述離散網(wǎng)格的油藏壓力，根據(jù)所述各流體相飽和度變化值和相應(yīng)的所述估計(jì)值求和得到所述離散網(wǎng)格的各流體相飽和度。

2.如權(quán)利要求1所述的黑油油藏模擬方法，其特征在于，所述多個(gè)組分滲流參數(shù)包括油藏壓力、氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度及油組分在油相中摩爾濃度；

基于多個(gè)組分滲流參數(shù)建立描述黑油油藏流體滲流規(guī)律的組分模型，包括：

基于油藏壓力、氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度及油組分在油相中摩爾濃度中的一個(gè)或多個(gè)，建立氣組分守恒方程、油組分守恒方程、水組分守恒方程、氣組分相平衡方程、油組分相平衡方程、氣相內(nèi)部組分約束方程及油相內(nèi)部組分約束方程，作為所述組分模型。

3.如權(quán)利要求2所述的黑油油藏模擬方法，其特征在于，根據(jù)基于黑油模型的黑油模擬所需數(shù)據(jù)對(duì)所述組分模型進(jìn)行數(shù)值離散，得到基于所述黑油油藏的地質(zhì)模型離散網(wǎng)格的組分?jǐn)?shù)值模型，包括：

根據(jù)所述黑油模擬所需數(shù)據(jù)，基于所述黑油油藏的地質(zhì)模型離散網(wǎng)格，對(duì)偏微分方程形式的所述氣組分守恒方程、所述油組分守恒方程及所述水組分守恒方程進(jìn)行空間數(shù)值離散，得到相應(yīng)的非線(xiàn)性代數(shù)方程；

根據(jù)所述氣組分相平衡方程、所述油組分相平衡方程、所述氣相內(nèi)部組分約束方程、所述油相內(nèi)部組分約束方程及所有所述非線(xiàn)性代數(shù)方程，按牛頓迭代法計(jì)算得到包含雅可比矩陣單元的線(xiàn)性代數(shù)方程組，作為所述離散網(wǎng)格的組分?jǐn)?shù)值模型。

4.如權(quán)利要求1所述的黑油油藏模擬方法，其特征在于，利用所述黑油模擬所需數(shù)據(jù)中的流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表，包括：

根據(jù)所述流體PVT數(shù)據(jù)中的油相密度、油相中溶解氣油比、氣相密度及氣相中揮發(fā)油氣比，確定油相中溶解氣組分與油組分的摩爾比和氣相中揮發(fā)油組分與氣組分的摩爾比；

根據(jù)所述的油相中溶解氣組分與油組分的摩爾比和所述的氣相中揮發(fā)油組分與氣組分的摩爾比，計(jì)算得到氣組分平衡常數(shù)和油組分平衡常數(shù)；

根據(jù)所述的氣組分平衡常數(shù)和油組分平衡常數(shù)，計(jì)算得到氣組分總摩爾分?jǐn)?shù)和油組分總摩爾分?jǐn)?shù)，所述組分相平衡參數(shù)表包括所述氣組分平衡常數(shù)、所述油組分平衡常數(shù)、所述氣組分總摩爾分?jǐn)?shù)、所述油組分總摩爾分?jǐn)?shù)及所述流體PVT數(shù)據(jù)中相應(yīng)的飽和壓力。

5.如權(quán)利要求1所述的黑油油藏模擬方法，其特征在于，根據(jù)所述組分相平衡參數(shù)表和所述組分滲流參數(shù)的估計(jì)值進(jìn)行閃蒸計(jì)算，得到所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)之前，還包括：

根據(jù)所述黑油模擬所需數(shù)據(jù)中的所述黑油油藏的地質(zhì)模型確定所述離散網(wǎng)格的各流體相飽和度初始值，利用靜水柱壓力公式計(jì)算所述離散網(wǎng)格的油藏壓力初始值，通過(guò)查詢(xún)所述組分相平衡參數(shù)表得到所述離散網(wǎng)格的各組分總摩爾分?jǐn)?shù)初始值，各流體相飽和度初始值、油藏壓力初始值及各組分總摩爾分?jǐn)?shù)初始值一同作為初始的所述多個(gè)組分滲流參數(shù)的估計(jì)值。

6.如權(quán)利要求1所述的黑油油藏模擬方法，其特征在于，根據(jù)所述組分相平衡參數(shù)表和所述組分滲流參數(shù)的估計(jì)值進(jìn)行閃蒸計(jì)算，得到所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)，包括：

根據(jù)所述離散網(wǎng)格的所述組分滲流參數(shù)的當(dāng)前估計(jì)值計(jì)算得到所述離散網(wǎng)格的組分總摩爾分?jǐn)?shù)；

根據(jù)所述組分總摩爾分?jǐn)?shù)查詢(xún)所述組分相平衡參數(shù)表，得到所述離散網(wǎng)格中組分構(gòu)成所對(duì)應(yīng)的飽和壓力；

通過(guò)對(duì)比所述離散網(wǎng)格的飽和壓力和所述組分滲流參數(shù)中的油藏壓力的當(dāng)前估計(jì)值，確定所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)。

7.如權(quán)利要求6所述的黑油油藏模擬方法，其特征在于，通過(guò)對(duì)比所述離散網(wǎng)格的飽和壓力和所述組分滲流參數(shù)中的油藏壓力的當(dāng)前估計(jì)值，確定所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)，包括：

根據(jù)所述離散網(wǎng)格的油藏壓力的當(dāng)前估計(jì)值大于所述飽和壓力判定所述離散網(wǎng)格的流體相態(tài)為油相；

當(dāng)所述離散網(wǎng)格的油藏壓力的當(dāng)前估計(jì)值小于所述飽和壓力時(shí)，根據(jù)所述離散網(wǎng)格的油藏壓力的當(dāng)前估計(jì)值查詢(xún)所述組分相平衡參數(shù)表，得到所述離散網(wǎng)格的氣組分平衡常數(shù)和油組分平衡常數(shù)，并根據(jù)所述的氣組分平衡常數(shù)和油組分平衡常數(shù)計(jì)算所述離散網(wǎng)格的液相比例；

根據(jù)所述液相比例大于零判定所述離散網(wǎng)格的流體相態(tài)為油氣兩相共存，根據(jù)所述液相比例等于或小于零判定所述離散網(wǎng)格的流體相態(tài)為氣相；

根據(jù)所述離散網(wǎng)格的流體相態(tài)的判定結(jié)果和油藏壓力的當(dāng)前估計(jì)值，查詢(xún)所述流體PVT數(shù)據(jù)得到相應(yīng)物性參數(shù)。

8.如權(quán)利要求1所述的黑油油藏模擬方法，其特征在于，還包括：

若所述的油藏壓力變化值及各流體相飽和度變化值不收斂，則用所述油藏壓力變化值和所述各流體相飽和度變化值作為相應(yīng)的所述估計(jì)值的增量更新所述估計(jì)值，根據(jù)所述組分相平衡參數(shù)表和更新后的估計(jì)值重新進(jìn)行閃蒸計(jì)算得到所述離散網(wǎng)格的流體相態(tài)及相應(yīng)物性參數(shù)，根據(jù)重新計(jì)算的流體相態(tài)及相應(yīng)物性參數(shù)重新調(diào)整組分?jǐn)?shù)值模型，利用更新后的估計(jì)值和重新調(diào)整后的組分?jǐn)?shù)值模型重新計(jì)算油藏壓力變化值及各流體相飽和度變化值，依次迭代計(jì)算，直到重新計(jì)算的油藏壓力變化值及各流體相飽和度變化值收斂。

9.如權(quán)利要求8所述的黑油油藏模擬方法，其特征在于，判斷所述的油藏壓力變化值及各流體相飽和度變化值是否收斂，包括：

根據(jù)變化值δX^l和相應(yīng)估計(jì)值X^l是否滿(mǎn)足判定所述的油藏壓力變化值及各流體相飽和度變化值收斂，其中，變化值δX^l表示第l次迭代后油藏壓力變化值或各流體相飽和度變化值，相應(yīng)估計(jì)值X^l表示前次迭代后油藏壓力的估計(jì)值或各流體相飽和度的估計(jì)值，迭代次數(shù)l≥1，X¹為初始的估計(jì)值，a為設(shè)定值。

10.如權(quán)利要求2所述的黑油油藏模擬方法，其特征在于，

所述氣組分守恒方程為：

所述油組分守恒方程為：

所述水組分守恒方程為：

所述氣組分相平衡方程為：

$F_f^g(X_{gg},X_{go})=X_{gg}{P}_g-K_g{X}_{go}{P}_g=X_{gg}-K_g{X}_{go}=0,$

所述油組分相平衡方程為：

$F_f^o(X_{og},X_{oo})=X_{og}{P}_g-K_o{X}_{oo}{P}_g=X_{og}-K_o{X}_{oo}=0,$

所述氣相內(nèi)部組分約束方程為：

$F_c^g(X_{og},X_{gg})=X_{og}+X_{gg}-1=0,$

所述油相內(nèi)部組分約束方程為：

$F_c^o(X_{oo},X_{go})=X_{oo}+X_{go}-1=0,$

其中，F(xiàn)_g表示描述氣組分守恒關(guān)系的函數(shù)，F(xiàn)_o表示描述油組分守恒關(guān)系的函數(shù)，F(xiàn)_w表示描述水組分守恒關(guān)系的函數(shù)，描述氣組分相平衡關(guān)系的函數(shù)，表示描述油組分相平衡關(guān)系的函數(shù)，表示描述氣相內(nèi)部油氣組分約束關(guān)系的函數(shù)，表示描述油相內(nèi)部油氣組分約束關(guān)系的函數(shù)；P_g為氣相壓力，S_g為氣相飽和度，S_o為油相飽和度，X_gg為氣組分在氣相中摩爾濃度，X_go為氣組分在油相中摩爾濃度，X_og為油組分在氣相中摩爾濃度，X_oo為油組分在油相中的摩爾濃度；t為時(shí)間，為孔隙度，ρ_g為氣相密度，ρ_o為油相密度，K為滲透率，K_rg為氣相相對(duì)滲透率，μ_g為氣相粘度，g為重力加速度，D為油藏深度，K_ro為油相相對(duì)滲透率，μ_o為油相粘度，P_cgo為氣油兩相之間的毛管力，q_g為氣相產(chǎn)量或注入量，q_o為油相產(chǎn)量或注入量，ρ_w為水相密度，K_rw為水相相對(duì)滲透率，μ_w為水相粘度，P_cow為油水兩相之間的毛管力，q_w為水相產(chǎn)量或注入量。

11.如權(quán)利要求3所述的黑油油藏模擬方法，其特征在于，

所述氣組分守恒方程所對(duì)應(yīng)的非線(xiàn)性代數(shù)方程為：

所述油組分守恒方程所對(duì)應(yīng)的非線(xiàn)性代數(shù)方程為：

所述水組分守恒方程所對(duì)應(yīng)的非線(xiàn)性代數(shù)方程為：

其中，F(xiàn)_g表示描述氣組分守恒關(guān)系的函數(shù)，F(xiàn)_o表示描述油組分守恒關(guān)系的函數(shù)，F(xiàn)_w表示描述水組分守恒關(guān)系的函數(shù)，P_g為氣相壓力，S_g為氣相飽和度，S_o為油相飽和度，X_gg為氣組分在氣相中摩爾濃度，X_go為氣組分在油相中摩爾濃度，X_og為油組分在氣相中摩爾濃度，X_oo為油組分在油相中的摩爾濃度；Ω為空間積分變量，為孔隙度，ρ_g為氣相密度，ρ_o為油相密度，K為滲透率，K_rg為氣相相對(duì)滲透率，μ_g為氣相粘度，g為重力加速度，D為油藏深度，K_ro為油相相對(duì)滲透率，μ_o為油相粘度，P_cgo為氣油兩相之間的毛管力，q_g為氣相產(chǎn)量或注入量，q_o為油相產(chǎn)量或注入量，ρ_w為水相密度，K_rw為水相相對(duì)滲透率，μ_w為水相粘度，P_cow為油水兩相之間的毛管力，q_w為水相產(chǎn)量或注入量，∑_l′表示對(duì)離散網(wǎng)格所有表面l′的通量求和，A_l′表示離散網(wǎng)格表面l′的面積，V表示離散網(wǎng)格的體積，Q_g＝V×q_g表示離散網(wǎng)格的氣相注入或產(chǎn)出率，Q_o＝V×q_o表示離散網(wǎng)格的油相注入或產(chǎn)出率，Q_w＝V×q_w表示離散網(wǎng)格的水相注入或產(chǎn)出率。

12.如權(quán)利要求3所述的黑油油藏模擬方法，其特征在于，所述線(xiàn)性代數(shù)方程組為：

${\left(\frac{\∂F}{\∂X}\right)}^l\mathrm{\δ}X=-F^l,$

其中，F(xiàn)表示函數(shù)，X表示變量，δX為變化量，表示雅可比矩陣；

其中，F(xiàn)_g表示描述氣組分守恒關(guān)系的函數(shù)，F(xiàn)_o表示描述油組分守恒關(guān)系的函數(shù)，F(xiàn)_w表示描述水組分守恒關(guān)系的函數(shù)，描述氣組分相平衡關(guān)系的函數(shù)，表示描述油組分相平衡關(guān)系的函數(shù)，表示描述氣相內(nèi)部油氣組分約束關(guān)系的函數(shù)，表示描述油相內(nèi)部油氣組分約束關(guān)系的函數(shù)，P_g為氣相壓力，S_g為氣相飽和度，S_o為油相飽和度，X_gg為氣組分在氣相中摩爾濃度，X_go為氣組分在油相中摩爾濃度，X_og為油組分在氣相中摩爾濃度，X_oo為油組分在油相中的摩爾濃度，l用于表示求解線(xiàn)性代數(shù)方程組的不同次數(shù)。

13.如權(quán)利要求12所述的黑油油藏模擬方法，其特征在于，在調(diào)整后的所述組分?jǐn)?shù)值模型中：

當(dāng)所述離散網(wǎng)格的流體相態(tài)為油氣兩相共存時(shí)，變量所述雅可比矩陣單元為：

當(dāng)所述離散網(wǎng)格的流體相態(tài)僅為油相時(shí)，變量所述雅可比矩陣單元為：

當(dāng)所述離散網(wǎng)格的流體相態(tài)僅為氣相時(shí)，變量所述雅可比矩陣單元為：

14.如權(quán)利要求6所述的黑油油藏模擬方法，其特征在于，所述離散網(wǎng)格的組分包括油組分和氣組分，

油組分的總摩爾分?jǐn)?shù)為：

$Z_o=\frac{(S_g{\mathrm{\ρ}}_g{X}_{og}+S_o{\mathrm{\ρ}}_o{X}_{oo})}{(S_g{\mathrm{\ρ}}_g{X}_{og}+S_o{\mathrm{\ρ}}_o{X}_{oo}+S_g{\mathrm{\ρ}}_g{X}_{gg}+S_o{\mathrm{\ρ}}_o{X}_{go})},$

氣組分的總摩爾分?jǐn)?shù)為：

$Z_g=\frac{(S_g{\mathrm{\ρ}}_g{X}_{gg}+S_o{\mathrm{\ρ}}_o{X}_{go})}{(S_g{\mathrm{\ρ}}_g{X}_{og}+S_o{\mathrm{\ρ}}_o{X}_{oo}+S_g{\mathrm{\ρ}}_g{X}_{gg}+S_o{\mathrm{\ρ}}_o{X}_{go})}=1-Z_o,$

其中，Z_o為油組分的總摩爾分?jǐn)?shù)，Z_g為氣組分的總摩爾分?jǐn)?shù)，S_g為氣相飽和度，S_o為油相飽和度，ρ_g為氣密度，ρ_o為油密度，X_og為油組分在氣相中摩爾濃度，X_oo為油組分在油相中的摩爾濃度，X_gg為氣組分在氣相中摩爾濃度，X_go為氣組分在油相中摩爾濃度。

15.如權(quán)利要求7所述的黑油油藏模擬方法，其特征在于，所述液相比例為：

$V^{\′}=\frac{Z_g}{1-K_o}+\frac{Z_o}{1-K_g},$

其中，V'表示液相比例，Z_g為氣組分的總摩爾分?jǐn)?shù)，Z_o為油組分的總摩爾分?jǐn)?shù)，K_o油組分平衡常數(shù)，K_g為氣組分平衡常數(shù)。

16.一種黑油油藏模擬裝置，其特征在于，包括：

組分模型建立單元，用于執(zhí)行：基于多個(gè)組分滲流參數(shù)建立描述黑油油藏流體滲流規(guī)律的組分模型；

組分?jǐn)?shù)值模型生成單元，用于執(zhí)行：根據(jù)基于黑油模型的黑油模擬所需數(shù)據(jù)對(duì)所述組分模型進(jìn)行數(shù)值離散，得到基于所述黑油油藏的地質(zhì)模型離散網(wǎng)格的組分?jǐn)?shù)值模型；

組分相平衡參數(shù)表建立單元，用于執(zhí)行：利用所述黑油模擬所需數(shù)據(jù)中的流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表；

流體相態(tài)確定單元，用于執(zhí)行：根據(jù)所述組分相平衡參數(shù)表和所述組分滲流參數(shù)的估計(jì)值進(jìn)行閃蒸計(jì)算，得到所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)；

組分?jǐn)?shù)值模型調(diào)整單元，用于執(zhí)行：根據(jù)所述流體相態(tài)和相應(yīng)物性參數(shù)調(diào)整所述組分?jǐn)?shù)值模型，以簡(jiǎn)化組分?jǐn)?shù)值模型；

油藏參數(shù)變化值生成單元，用于執(zhí)行：利用所述組分滲流參數(shù)的估計(jì)值和調(diào)整后的所述組分?jǐn)?shù)值模型，計(jì)算得到油藏壓力變化值及各流體相飽和度變化值；

油藏參數(shù)確定單元，用于執(zhí)行：判斷所述的油藏壓力變化值及各流體相飽和度變化值是否收斂，若是，根據(jù)所述油藏壓力變化值和相應(yīng)的所述估計(jì)值求和得到所述離散網(wǎng)格的油藏壓力，根據(jù)所述各流體相飽和度變化值和相應(yīng)的所述估計(jì)值求和得到所述離散網(wǎng)格的各流體相飽和度。

17.如權(quán)利要求16所述的黑油油藏模擬裝置，其特征在于，所述多個(gè)組分滲流參數(shù)包括油藏壓力、氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度及油組分在油相中摩爾濃度，所述組分模型建立單元，包括：

組分模型建立模塊，用于執(zhí)行：基于油藏壓力、氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度及油組分在油相中摩爾濃度中的一個(gè)或多個(gè)，建立氣組分守恒方程、油組分守恒方程、水組分守恒方程、氣組分相平衡方程、油組分相平衡方程、氣相內(nèi)部組分約束方程及油相內(nèi)部組分約束方程，作為所述組分模型。

18.如權(quán)利要求17所述的黑油油藏模擬裝置，其特征在于，所述組分?jǐn)?shù)值模型生成單元，包括：

組分模型離散模塊，用于執(zhí)行：根據(jù)所述黑油模擬所需數(shù)據(jù)，基于所述黑油油藏的地質(zhì)模型離散網(wǎng)格，對(duì)偏微分方程形式的所述氣組分守恒方程、所述油組分守恒方程及所述水組分守恒方程進(jìn)行空間數(shù)值離散，得到相應(yīng)的非線(xiàn)性代數(shù)方程；

組分?jǐn)?shù)值模型生成模塊，用于執(zhí)行：根據(jù)所述氣組分相平衡方程、所述油組分相平衡方程、所述氣相內(nèi)部組分約束方程、所述油相內(nèi)部組分約束方程及所有所述非線(xiàn)性代數(shù)方程，按牛頓迭代法計(jì)算得到包含雅可比矩陣單元的線(xiàn)性代數(shù)方程組，作為所述離散網(wǎng)格的組分?jǐn)?shù)值模型。

19.如權(quán)利要求16所述的黑油油藏模擬裝置，其特征在于，所述組分相平衡參數(shù)表建立單元，包括：

組分摩爾比生成模塊，用于執(zhí)行：根據(jù)所述流體PVT數(shù)據(jù)中的油相密度、油相中溶解氣油比、氣相密度及氣相中揮發(fā)油氣比，確定油相中溶解氣組分與油組分的摩爾比和氣相中揮發(fā)油組分與氣組分的摩爾比；

組分平衡常數(shù)生成模塊，用于執(zhí)行：根據(jù)所述的油相中溶解氣組分與油組分的摩爾比和所述的氣相中揮發(fā)油組分與氣組分的摩爾比，計(jì)算得到氣組分平衡常數(shù)和油組分平衡常數(shù)；

組分相平衡參數(shù)表建立模塊，用于執(zhí)行：根據(jù)所述的氣組分平衡常數(shù)和油組分平衡常數(shù)，計(jì)算得到氣組分總摩爾分?jǐn)?shù)和油組分總摩爾分?jǐn)?shù)，所述組分相平衡參數(shù)表包括所述氣組分平衡常數(shù)、所述油組分平衡常數(shù)、所述氣組分總摩爾分?jǐn)?shù)、所述油組分總摩爾分?jǐn)?shù)及所述流體PVT數(shù)據(jù)中相應(yīng)的飽和壓力。

20.如權(quán)利要求16所述的黑油油藏模擬裝置，其特征在于，還包括：

初始估計(jì)值生成單元，用于執(zhí)行：根據(jù)所述黑油模擬所需數(shù)據(jù)中的所述黑油油藏的地質(zhì)模型確定所述離散網(wǎng)格的各流體相飽和度初始值，利用靜水柱壓力公式計(jì)算所述離散網(wǎng)格的油藏壓力初始值，通過(guò)查詢(xún)所述組分相平衡參數(shù)表得到所述離散網(wǎng)格的各組分總摩爾分?jǐn)?shù)初始值，各流體相飽和度初始值、油藏壓力初始值及各組分總摩爾分?jǐn)?shù)初始值一同作為初始的所述多個(gè)組分滲流參數(shù)的估計(jì)值。

21.如權(quán)利要求16所述的黑油油藏模擬裝置，其特征在于，所述流體相態(tài)確定單元，包括：

網(wǎng)格摩爾分?jǐn)?shù)生成模塊，用于執(zhí)行：根據(jù)所述離散網(wǎng)格的所述組分滲流參數(shù)的當(dāng)前估計(jì)值計(jì)算得到所述離散網(wǎng)格的組分總摩爾分?jǐn)?shù)；

網(wǎng)格飽和壓力生成模塊，用于執(zhí)行：根據(jù)所述組分總摩爾分?jǐn)?shù)查詢(xún)所述組分相平衡參數(shù)表，得到所述離散網(wǎng)格中組分構(gòu)成所對(duì)應(yīng)的飽和壓力；

網(wǎng)格流體相態(tài)確定模塊，用于執(zhí)行：通過(guò)對(duì)比所述離散網(wǎng)格的飽和壓力和所述組分滲流參數(shù)中的油藏壓力的當(dāng)前估計(jì)值，確定所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)。

22.如權(quán)利要求21所述的黑油油藏模擬裝置，其特征在于，所述網(wǎng)格流體相態(tài)確定模塊，包括：

油相相態(tài)判斷模塊，用于執(zhí)行：根據(jù)所述離散網(wǎng)格的油藏壓力的當(dāng)前估計(jì)值大于所述飽和壓力判定所述離散網(wǎng)格的流體相態(tài)為油相；

液相比例生成模塊，用于執(zhí)行：當(dāng)所述離散網(wǎng)格的油藏壓力的當(dāng)前估計(jì)值小于所述飽和壓力時(shí)，根據(jù)所述離散網(wǎng)格的油藏壓力的當(dāng)前估計(jì)值查詢(xún)所述組分相平衡參數(shù)表，得到所述離散網(wǎng)格的氣組分平衡常數(shù)和油組分平衡常數(shù)，并根據(jù)所述的氣組分平衡常數(shù)和油組分平衡常數(shù)計(jì)算所述離散網(wǎng)格的液相比例；

油氣兩相及氣相判斷模塊，用于執(zhí)行：根據(jù)所述液相比例大于零判定所述離散網(wǎng)格的流體相態(tài)為油氣兩相共存，根據(jù)所述液相比例等于或小于零判定所述離散網(wǎng)格的流體相態(tài)為氣相；

物性參數(shù)確定模塊，用于執(zhí)行：根據(jù)所述離散網(wǎng)格的流體相態(tài)的判定結(jié)果和油藏壓力的當(dāng)前估計(jì)值，查詢(xún)所述流體PVT數(shù)據(jù)得到相應(yīng)物性參數(shù)。

23.如權(quán)利要求16所述的黑油油藏模擬裝置，其特征在于，還包括：

迭代計(jì)算單元，用于執(zhí)行：若所述的油藏壓力變化值及各流體相飽和度變化值不收斂，則用所述油藏壓力變化值和所述各流體相飽和度變化值作為相應(yīng)的所述估計(jì)值的增量更新所述估計(jì)值，根據(jù)所述組分相平衡參數(shù)表和更新后的估計(jì)值重新進(jìn)行閃蒸計(jì)算得到所述離散網(wǎng)格的流體相態(tài)及相應(yīng)物性參數(shù)，根據(jù)重新計(jì)算的流體相態(tài)及相應(yīng)物性參數(shù)重新調(diào)整組分?jǐn)?shù)值模型，利用更新后的估計(jì)值和重新調(diào)整后的組分?jǐn)?shù)值模型重新計(jì)算油藏壓力變化值及各流體相飽和度變化值，依次迭代計(jì)算，直到重新計(jì)算的油藏壓力變化值及各流體相飽和度變化值收斂。

24.如權(quán)利要求23所述的黑油油藏模擬裝置，其特征在于，所述油藏參數(shù)確定單元，包括：

收斂判斷模塊，用于執(zhí)行：根據(jù)變化值δX^l和相應(yīng)估計(jì)值X^l是否滿(mǎn)足判定所述的油藏壓力變化值及各流體相飽和度變化值收斂，其中，變化值δX^l表示第l次迭代后油藏壓力變化值或各流體相飽和度變化值，相應(yīng)估計(jì)值X^l表示前次迭代后油藏壓力的估計(jì)值或各流體相飽和度的估計(jì)值，迭代次數(shù)l≥1，X¹為初始的估計(jì)值，a為設(shè)定值。