本發(fā)明涉及油田開(kāi)發(fā)
技術(shù)領(lǐng)域：
，尤其涉及一種黑油油藏模擬方法及裝置。
背景技術(shù)：
：在油氣田開(kāi)發(fā)過(guò)程中，油藏?cái)?shù)值模擬技術(shù)及裝置主要用于認(rèn)識(shí)和預(yù)測(cè)地下油氣資源的分布和運(yùn)移情況。根據(jù)等溫條件下流體物性不同，油藏?cái)?shù)值模擬可分為黑油模擬和組分模擬。黑油模擬是基于黑油模型實(shí)現(xiàn)，該模型描述油氣水三相三組分滲流問(wèn)題，隨著油藏壓力的變化，油氣組分可單獨(dú)或者共存在于氣/液兩相中；黑油模擬計(jì)算量較少，模擬速度快，但考慮的因素相對(duì)簡(jiǎn)單；對(duì)于混相氣驅(qū)、凝析氣藏等更復(fù)雜的油氣藏開(kāi)發(fā)情況，需要采用組分模擬來(lái)取得更準(zhǔn)確的預(yù)測(cè)結(jié)果。組分模擬基于組分模型來(lái)實(shí)現(xiàn)，該模型考慮了在油藏壓力、組分構(gòu)成變化情況下氣液相平衡關(guān)系及其對(duì)流體粘度、密度等物性參數(shù)的影響，因此，在組分模擬過(guò)程中需要額外通過(guò)閃蒸計(jì)算來(lái)了解流體相態(tài)及物性變化，計(jì)算量大，模擬速度慢?，F(xiàn)有油藏?cái)?shù)值模擬技術(shù)是基于黑油模型開(kāi)發(fā)黑油模擬系統(tǒng)，基于組分模型開(kāi)發(fā)組分模擬系統(tǒng)。這給油藏模擬技術(shù)后續(xù)發(fā)展帶來(lái)了模塊化程度低的弊端。后續(xù)模擬功能升級(jí)都需要分別在黑油模擬裝置和組分模擬裝置中重復(fù)開(kāi)發(fā)實(shí)現(xiàn)，這些功能升級(jí)包括但不局限于：雙重介質(zhì)模擬功能、離散裂縫功能、復(fù)雜結(jié)構(gòu)井模擬、滑脫/擴(kuò)散/吸附/解吸附等因素模擬。這大大提升了新功能重復(fù)開(kāi)發(fā)的工作量，而直接用組分模擬代替黑油模擬，需要額外輸入組分相平衡參數(shù)來(lái)完成復(fù)雜的閃蒸計(jì)算，用戶便利性差，計(jì)算資源消耗大。技術(shù)實(shí)現(xiàn)要素：本發(fā)明提供一種黑油油藏模擬方法及裝置，以提高油藏模擬的可擴(kuò)展性。本發(fā)明提供一種黑油油藏模擬方法，包括：基于多個(gè)組分滲流參數(shù)建立描述黑油油藏流體滲流規(guī)律的組分模型；根據(jù)基于黑油模型的黑油模擬所需數(shù)據(jù)對(duì)所述組分模型進(jìn)行數(shù)值離散，得到基于所述黑油油藏的地質(zhì)模型離散網(wǎng)格的組分?jǐn)?shù)值模型；利用所述黑油模擬所需數(shù)據(jù)中的流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表；根據(jù)所述組分相平衡參數(shù)表和所述離散網(wǎng)格的所述組分滲流參數(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)格的各流體相飽和度。一個(gè)實(shí)施例中，所述多個(gè)組分滲流參數(shù)包括油藏壓力、氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度及油組分在油相中摩爾濃度，基于多個(gè)組分滲流參數(shù)建立描述黑油油藏流體滲流規(guī)律的組分模型，包括：基于油藏壓力、氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度及油組分在油相中摩爾濃度中的一個(gè)或多個(gè)，建立氣組分守恒方程、油組分守恒方程、水組分守恒方程、氣組分相平衡方程、油組分相平衡方程、氣相內(nèi)部組分約束方程及油相內(nèi)部組分約束方程，作為所述組分模型。一個(gè)實(shí)施例中，根據(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)的非線性代數(shù)方程；根據(jù)所述氣組分相平衡方程、所述油組分相平衡方程、所述氣相內(nèi)部組分約束方程、所述油相內(nèi)部組分約束方程及所有所述非線性代數(shù)方程，按牛頓迭代法計(jì)算得到包含雅可比矩陣單元的線性代數(shù)方程組，作為所述離散網(wǎng)格的組分?jǐn)?shù)值模型。一個(gè)實(shí)施例中，利用所述黑油模擬所需數(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)的飽和壓力。一個(gè)實(shí)施例中，根據(jù)所述組分相平衡參數(shù)表和所述組分滲流參數(shù)的估計(jì)值進(jìn)行閃蒸計(jì)算，得到所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)之前，還包括：根據(jù)所述黑油模擬所需數(shù)據(jù)中的所述黑油油藏的地質(zhì)模型確定所述離散網(wǎng)格的各流體相飽和度初始值，利用靜水柱壓力公式計(jì)算所述離散網(wǎng)格的油藏壓力初始值，通過(guò)查詢所述組分相平衡參數(shù)表得到所述離散網(wǎng)格的各組分總摩爾分?jǐn)?shù)初始值，各流體相飽和度初始值、油藏壓力初始值及各組分總摩爾分?jǐn)?shù)初始值一同作為初始的所述多個(gè)組分滲流參數(shù)的估計(jì)值。一個(gè)實(shí)施例中，根據(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ù)查詢所述組分相平衡參數(shù)表，得到所述離散網(wǎng)格中組分構(gòu)成所對(duì)應(yīng)的飽和壓力；通過(guò)對(duì)比所述離散網(wǎng)格的飽和壓力和所述組分滲流參數(shù)中的油藏壓力的當(dāng)前估計(jì)值，確定所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)。一個(gè)實(shí)施例中，通過(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ì)值查詢所述組分相平衡參數(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ì)值，查詢所述流體PVT數(shù)據(jù)得到相應(yīng)物性參數(shù)。一個(gè)實(shí)施例中，還包括：若所述的油藏壓力變化值及各流體相飽和度變化值不收斂，則用所述油藏壓力變化值和所述各流體相飽和度變化值作為相應(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ì)算的油藏壓力變化值及各流體相飽和度變化值收斂。一個(gè)實(shí)施例中，判斷所述的油藏壓力變化值及各流體相飽和度變化值是否收斂，包括：根據(jù)變化值δXl和相應(yīng)估計(jì)值Xl是否滿足判定所述的油藏壓力變化值及各流體相飽和度變化值收斂，其中，變化值δXl表示第l次迭代后油藏壓力變化值或各流體相飽和度變化值，相應(yīng)估計(jì)值Xl表示前次迭代后油藏壓力的估計(jì)值或各流體相飽和度的估計(jì)值，迭代次數(shù)l≥1，X1為初始的估計(jì)值，a為設(shè)定值。一個(gè)實(shí)施例中，所述氣組分守恒方程為：所述油組分守恒方程為：所述水組分守恒方程為：所述氣組分相平衡方程為：Ffg(Xgg,Xgo)=XggPg-KgXgoPg=Xgg-KgXgo=0,]]>所述油組分相平衡方程為：Ffo(Xog,Xoo)=XogPg-KoXooPg=Xog-KoXoo=0,]]>所述氣相內(nèi)部組分約束方程為：Fcg(Xog,Xgg)=Xog+Xgg-1=0,]]>所述油相內(nèi)部組分約束方程為：Fco(Xoo,Xgo)=Xoo+Xgo-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ù)；Pg為氣相壓力，Sg為氣相飽和度，So為油相飽和度，Xgg為氣組分在氣相中摩爾濃度，Xgo為氣組分在油相中摩爾濃度，Xog為油組分在氣相中摩爾濃度，Xoo為油組分在油相中的摩爾濃度；t為時(shí)間，為孔隙度，ρg為氣相密度，ρo為油相密度，K為滲透率，Krg為氣相相對(duì)滲透率，μg為氣相粘度，g為重力加速度，D為油藏深度，Kro為油相相對(duì)滲透率，μo為油相粘度，Pcgo為氣油兩相之間的毛管力，qg為氣相產(chǎn)量或注入量，qo為油相產(chǎn)量或注入量，ρw為水相密度，Krw為水相相對(duì)滲透率，μw為水相粘度，Pcow為油水兩相之間的毛管力，qw為水相產(chǎn)量或注入量。一個(gè)實(shí)施例中，所述氣組分守恒方程所對(duì)應(yīng)的非線性代數(shù)方程為：所述油組分守恒方程所對(duì)應(yīng)的非線性代數(shù)方程為：所述水組分守恒方程所對(duì)應(yīng)的非線性代數(shù)方程為：其中，F(xiàn)g表示描述氣組分守恒關(guān)系的函數(shù)，F(xiàn)o表示描述油組分守恒關(guān)系的函數(shù)，F(xiàn)w表示描述水組分守恒關(guān)系的函數(shù)，Pg為氣相壓力，Sg為氣相飽和度，So為油相飽和度，Xgg為氣組分在氣相中摩爾濃度，Xgo為氣組分在油相中摩爾濃度，Xog為油組分在氣相中摩爾濃度，Xoo為油組分在油相中的摩爾濃度；Ω為空間積分變量，為孔隙度，ρg為氣相密度，ρo為油相密度，K為滲透率，Krg為氣相相對(duì)滲透率，μg為氣相粘度，g為重力加速度，D為油藏深度，Kro為油相相對(duì)滲透率，μo為油相粘度，Pcgo為氣油兩相之間的毛管力，qg為氣相產(chǎn)量或注入量，qo為油相產(chǎn)量或注入量，ρw為水相密度，Krw為水相相對(duì)滲透率，μw為水相粘度，Pcow為油水兩相之間的毛管力，qw為水相產(chǎn)量或注入量，∑l′表示對(duì)離散網(wǎng)格所有表面l′l'的通量求和，l′表示離散網(wǎng)格表面l′的面積，V表示離散網(wǎng)格的體積，Qg＝V×qg表示離散網(wǎng)格的氣相注入或產(chǎn)出率，Qo＝V×qo表示離散網(wǎng)格的油相注入或產(chǎn)出率，Qw＝V×qw表示離散網(wǎng)格的水相注入或產(chǎn)出率。一個(gè)實(shí)施例中，所述線性代數(shù)方程組為：(∂F∂X)lδX=-Fl,]]>其中，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ù)，Pg為氣相壓力，Sg為氣相飽和度，So為油相飽和度，Xgg為氣組分在氣相中摩爾濃度，Xgo為氣組分在油相中摩爾濃度，Xog為油組分在氣相中摩爾濃度，Xoo為油組分在油相中的摩爾濃度，l用于表示求解線性代數(shù)方程組的不同次數(shù)。一個(gè)實(shí)施例中，在調(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í)，變量所述雅可比矩陣單元為：一個(gè)實(shí)施例中，所述離散網(wǎng)格的組分包括油組分和氣組分，油組分的總摩爾分?jǐn)?shù)為：Zo=(SgρgXog+SoρoXoo)(SgρgXog+SoρoXoo+SgρgXgg+SoρoXgo),]]>氣組分的總摩爾分?jǐn)?shù)為：Zg=(SgρgXgg+SoρoXgo)(SgρgXog+SoρoXoo+SgρgXgg+SoρoXgo)=1-Zo,]]>其中，Zo為油組分的總摩爾分?jǐn)?shù)，Zg為氣組分的總摩爾分?jǐn)?shù)，Sg為氣相飽和度，So為油相飽和度，ρg為氣密度，ρo為油密度，Xog為油組分在氣相中摩爾濃度，Xoo為油組分在油相中的摩爾濃度，Xgg為氣組分在氣相中摩爾濃度，Xgo為氣組分在油相中摩爾濃度。一個(gè)實(shí)施例中，所述液相比例為：V′=Zg1-Ko+Zo1-Kg,]]>其中，V'表示液相比例，Zg為氣組分的總摩爾分?jǐn)?shù)，Zo為油組分的總摩爾分?jǐn)?shù)，Ko油組分平衡常數(shù)，Kg為氣組分平衡常數(shù)。本發(fā)明還提供一種黑油油藏模擬裝置，包括：組分模型建立單元，用于執(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)格的各流體相飽和度。一個(gè)實(shí)施例中，所述多個(gè)組分滲流參數(shù)包括油藏壓力、氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度及油組分在油相中摩爾濃度，所述組分模型建立單元，包括：組分模型建立模塊，用于執(zhí)行：基于油藏壓力、氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度及油組分在油相中摩爾濃度中的一個(gè)或多個(gè)，建立氣組分守恒方程、油組分守恒方程、水組分守恒方程、氣組分相平衡方程、油組分相平衡方程、氣相內(nèi)部組分約束方程及油相內(nèi)部組分約束方程，作為所述組分模型。一個(gè)實(shí)施例中，所述組分?jǐn)?shù)值模型生成單元，包括：組分模型離散模塊，用于執(zhí)行：根據(jù)所述黑油模擬所需數(shù)據(jù)，基于所述黑油油藏的地質(zhì)模型離散網(wǎng)格，對(duì)偏微分方程形式的所述氣組分守恒方程、所述油組分守恒方程及所述水組分守恒方程進(jìn)行空間數(shù)值離散，得到相應(yīng)的非線性代數(shù)方程；組分?jǐn)?shù)值模型生成模塊，用于執(zhí)行：根據(jù)所述氣組分相平衡方程、所述油組分相平衡方程、所述氣相內(nèi)部組分約束方程、所述油相內(nèi)部組分約束方程及所有所述非線性代數(shù)方程，按牛頓迭代法計(jì)算得到包含雅可比矩陣單元的線性代數(shù)方程組，作為所述離散網(wǎng)格的組分?jǐn)?shù)值模型。一個(gè)實(shí)施例中，所述組分相平衡參數(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)的飽和壓力。一個(gè)實(shí)施例中，還包括：初始估計(jì)值生成單元，用于執(zhí)行：根據(jù)所述黑油模擬所需數(shù)據(jù)中的所述黑油油藏的地質(zhì)模型確定所述離散網(wǎng)格的各流體相飽和度初始值，利用靜水柱壓力公式計(jì)算所述離散網(wǎng)格的油藏壓力初始值，通過(guò)查詢所述組分相平衡參數(shù)表得到所述離散網(wǎng)格的各組分總摩爾分?jǐn)?shù)初始值，各流體相飽和度初始值、油藏壓力初始值及各組分總摩爾分?jǐn)?shù)初始值一同作為初始的所述多個(gè)組分滲流參數(shù)的估計(jì)值。一個(gè)實(shí)施例中，所述流體相態(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ù)查詢所述組分相平衡參數(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ù)。一個(gè)實(shí)施例中，所述網(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ì)值查詢所述組分相平衡參數(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ì)值，查詢所述流體PVT數(shù)據(jù)得到相應(yīng)物性參數(shù)。一個(gè)實(shí)施例中，還包括：迭代計(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ì)算的油藏壓力變化值及各流體相飽和度變化值收斂。一個(gè)實(shí)施例中，所述油藏參數(shù)確定單元，包括：收斂判斷模塊，用于執(zhí)行：根據(jù)變化值δXl和相應(yīng)估計(jì)值Xl是否滿足判定所述的油藏壓力變化值及各流體相飽和度變化值收斂，其中，變化值δXl表示第l次迭代后油藏壓力變化值或各流體相飽和度變化值，相應(yīng)估計(jì)值Xl表示前次迭代后油藏壓力的估計(jì)值或各流體相飽和度的估計(jì)值，迭代次數(shù)l≥1，X1為初始的估計(jì)值，a為設(shè)定值。本發(fā)明的黑油油藏模擬方法及裝置，獨(dú)創(chuàng)地建立描述黑油油藏流體滲流規(guī)律的組分模型，進(jìn)一步得到離散網(wǎng)格的組分?jǐn)?shù)值模型，奠定了用組分模型進(jìn)行黑油模擬的基礎(chǔ)，再利用流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表，為用組分模型進(jìn)行黑油模擬搭建了橋梁，實(shí)現(xiàn)了組分模擬和現(xiàn)有基于黑油模型黑油模擬的有效結(jié)合。本發(fā)明的黑油油藏模擬方法及裝置能夠在不需要用戶額外輸入組分相平衡參數(shù)和保證計(jì)算效率情況下，在組分模型上實(shí)現(xiàn)黑油模型功能，從而提高油藏模擬的可擴(kuò)展性。附圖說(shuō)明為了更清楚地說(shuō)明本發(fā)明實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作簡(jiǎn)單地介紹，顯而易見(jiàn)地，下面描述中的附圖僅僅是本發(fā)明的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來(lái)講，在不付出創(chuàng)造性勞動(dòng)的前提下，還可以根據(jù)這些附圖獲得其他的附圖。在附圖中：圖1是本發(fā)明一實(shí)施例的黑油油藏模擬方法的流程示意圖；圖2是本發(fā)明一實(shí)施例中對(duì)組分模型進(jìn)行數(shù)值離散得到組分?jǐn)?shù)值模型的方法的流程示意圖；圖3是本發(fā)明一實(shí)施例中的利用流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表的方法示意圖；圖4是本發(fā)明另一實(shí)施例的黑油油藏模擬方法的流程示意圖；圖5是本發(fā)明一實(shí)施例中進(jìn)行閃蒸計(jì)算得到離散網(wǎng)格流體相態(tài)的方法流程示意圖；圖6是本發(fā)明一實(shí)施例中通過(guò)對(duì)比飽和壓力和油藏壓力的估計(jì)值確定離散網(wǎng)格流體相態(tài)的方法流程示意圖；圖7是本發(fā)明又一實(shí)施例的黑油模擬方法的流程示意圖；圖8是本發(fā)明另一實(shí)施例中進(jìn)行閃蒸計(jì)算得到離散網(wǎng)格流體相態(tài)的方法流程示意圖；圖9是利用本發(fā)明實(shí)施例的黑油模擬方法得到的油相飽和度示意圖；圖10是本發(fā)明一實(shí)施例的黑油油藏模擬裝置的結(jié)構(gòu)示意圖；圖11是本發(fā)明一實(shí)施例中組分?jǐn)?shù)值模型生成單元的結(jié)構(gòu)示意圖；圖12是本發(fā)明一實(shí)施例中組分相平衡參數(shù)表建立單元的結(jié)構(gòu)示意圖；圖13是本發(fā)明另一實(shí)施例中黑油油藏模擬裝置的結(jié)構(gòu)示意圖；圖14是本發(fā)明一實(shí)施例中流體相態(tài)確定單元的結(jié)構(gòu)示意圖；圖15是本發(fā)明一實(shí)施例中網(wǎng)格流體相態(tài)確定模塊的結(jié)構(gòu)示意圖。具體實(shí)施方式為使本發(fā)明實(shí)施例的目的、技術(shù)方案和優(yōu)點(diǎn)更加清楚明白，下面結(jié)合附圖對(duì)本發(fā)明實(shí)施例做進(jìn)一步詳細(xì)說(shuō)明。在此，本發(fā)明的示意性實(shí)施例及其說(shuō)明用于解釋本發(fā)明，但并不作為對(duì)本發(fā)明的限定。本發(fā)明的黑油油藏模擬，針對(duì)現(xiàn)有黑油模擬和組分模擬所面臨的擴(kuò)展性差、使用不便等問(wèn)題，創(chuàng)造性地采用組分模型實(shí)現(xiàn)黑油模擬，能夠解決現(xiàn)有技術(shù)的問(wèn)題。圖1是本發(fā)明一實(shí)施例的黑油油藏模擬方法的流程示意圖。如圖1所示，本發(fā)明實(shí)施例的黑油油藏模擬方法，可包括步驟：S110：基于多個(gè)組分滲流參數(shù)建立描述黑油油藏流體滲流規(guī)律的組分模型；S120：根據(jù)基于黑油模型的黑油模擬所需數(shù)據(jù)對(duì)所述組分模型進(jìn)行數(shù)值離散，得到基于所述黑油油藏的地質(zhì)模型離散網(wǎng)格的組分?jǐn)?shù)值模型；S130：利用所述黑油模擬所需數(shù)據(jù)中的流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表；S140：根據(jù)所述組分相平衡參數(shù)表和所述離散網(wǎng)格的所述組分滲流參數(shù)的估計(jì)值進(jìn)行閃蒸計(jì)算，得到所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)；S150：根據(jù)所述流體相態(tài)和相應(yīng)物性參數(shù)調(diào)整所述組分?jǐn)?shù)值模型，以簡(jiǎn)化組分?jǐn)?shù)值模型；S160：利用所述組分滲流參數(shù)的估計(jì)值和調(diào)整后的所述組分?jǐn)?shù)值模型，計(jì)算得到油藏壓力變化值及各流體相飽和度變化值；S170：判斷所述的油藏壓力變化值及各流體相飽和度變化值是否收斂，若是，根據(jù)所述油藏壓力變化值和相應(yīng)的所述估計(jì)值求和得到所述離散網(wǎng)格的油藏壓力，根據(jù)所述各流體相飽和度變化值和相應(yīng)的所述估計(jì)值求和得到所述離散網(wǎng)格的各流體相飽和度。在上述步驟S110中，組分滲流參數(shù)可以是油藏壓力、油氣飽和度和油氣組分摩爾分?jǐn)?shù)等參數(shù)。不同組分可以是氣組分、油組分和水組分中的一個(gè)或多個(gè)。該組分模型可描述油氣水三組分黑油油藏流體滲流規(guī)律。在上述步驟S120中，上述基于黑油模型的黑油模擬是現(xiàn)有的黑油模擬方法。基于黑油模型的黑油模擬所需數(shù)據(jù)可包括：流體PVT數(shù)據(jù)、相滲數(shù)據(jù)、毛管力數(shù)據(jù)、地質(zhì)模型等。其中，流體PVT數(shù)據(jù)可包括：油相流體PVT數(shù)據(jù)、氣相流體PVT數(shù)據(jù)、油藏水相流體PVT數(shù)據(jù)。上述流體PVT數(shù)據(jù)可從實(shí)驗(yàn)室測(cè)量得到。上述相滲數(shù)據(jù)可以是相對(duì)滲透率曲線，上述毛管力數(shù)據(jù)可以是毛管力曲線，上述相滲數(shù)據(jù)和毛管力數(shù)據(jù)可同試驗(yàn)室測(cè)量得到。上述地質(zhì)模型可通過(guò)油藏地質(zhì)建模系統(tǒng)生成。該地質(zhì)模型可以是基于黑油模型的黑油模擬所需數(shù)據(jù)中的黑油油藏地質(zhì)模型。其中，上述油相流體PVT數(shù)據(jù)可包括：標(biāo)準(zhǔn)狀態(tài)下油相密度、不同飽和壓力取值下對(duì)應(yīng)的溶解氣油比、粘度及體積系數(shù)；上述氣相流體PVT數(shù)據(jù)可包括標(biāo)準(zhǔn)狀態(tài)下氣相密度、不同飽和壓力取值下對(duì)應(yīng)的揮發(fā)油氣比、粘度及體積系數(shù)；上述水相流體PVT數(shù)據(jù)可包括標(biāo)準(zhǔn)狀態(tài)下水相密度、不同壓力取值下對(duì)應(yīng)的粘度及體積系數(shù)。上述相對(duì)滲透率曲線可包括：油水兩相通過(guò)油藏巖心的相對(duì)滲透率、油氣兩相通過(guò)油藏巖心的相對(duì)滲透率。上述毛管力曲線可包括：不同含水飽和度下油水兩相毛管力、不同含水飽和度下油氣兩相毛管力。上述地質(zhì)模型可以是以網(wǎng)格單元方式描述油藏幾何形狀、油氣水分布、孔隙度分布情況的數(shù)據(jù)體。一些實(shí)施例中，除了采集上述黑油模擬所需數(shù)據(jù)外，根據(jù)研究需要還可利用其他數(shù)據(jù)對(duì)所述組分模型進(jìn)行數(shù)值離散，例如，還可利用額外采集或者輸入的生產(chǎn)井的信息和模擬控制信息。在上述步驟S120中，可以采用多種方法對(duì)組分模型進(jìn)行數(shù)值離散，數(shù)值離散后的組分模型可為非線性代數(shù)方程，不同的離散網(wǎng)格可對(duì)應(yīng)其各自的非線性代數(shù)方程。在上述步驟S130中，利用所述黑油模擬所需數(shù)據(jù)中的流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表，可以將流體PVT數(shù)據(jù)表轉(zhuǎn)換為組分相平衡參數(shù)表，組分相平衡參數(shù)表中的數(shù)據(jù)可與流體PVT數(shù)據(jù)具有一定的對(duì)應(yīng)關(guān)系，所以依據(jù)該組分相平衡參數(shù)表中的數(shù)據(jù)可以查詢得到流體PVT數(shù)據(jù)中的數(shù)據(jù)。在上述步驟S140中，該物性參數(shù)可以包括流體粘度、密度和體積系數(shù)等參數(shù)。組分滲流參數(shù)的估計(jì)值可以是初始設(shè)定的估計(jì)值，或者可以是根據(jù)黑油模擬計(jì)算結(jié)果計(jì)算得到，用于下次重新進(jìn)行黑油模擬。所述閃蒸計(jì)算是指黑油閃蒸技術(shù)，不等同于現(xiàn)有組分閃蒸計(jì)算，現(xiàn)有組分閃蒸計(jì)算需要求解非線性方程組，耗費(fèi)大量計(jì)算資源。本發(fā)明的黑油閃蒸計(jì)算可主要根據(jù)離散網(wǎng)格中的油氣組分構(gòu)成，通過(guò)查詢組分相平衡參數(shù)表得到流體相態(tài)和流體物性，該方式可大大提高了模擬效率。在上述步驟S160中，黑油模擬的結(jié)果通常包括油藏壓力和流體相飽和度信息，根據(jù)本發(fā)明的方法還可以得到其他多種信息。該各流體相飽和度變化值可以包括氣相飽和度變化值、油相飽和度變化值及水相飽和度變化值。在上述步驟S170中，離散網(wǎng)格的各流體相飽和度可包括氣相飽和度、油相飽和度及水相飽和度，三者之和為1，水相飽和度可根據(jù)氣相飽和度和油相飽和度計(jì)算得到。通過(guò)判斷油藏壓力變化值及各流體相飽和度變化值是否收斂，并根據(jù)收斂的油藏壓力變化值及各流體相飽和度變化值計(jì)算得到黑油模擬結(jié)果，以此可以保證黑油模擬結(jié)果的準(zhǔn)確度?？捎糜?jì)算得到的油藏壓力和流體相飽和度替換之前的相應(yīng)參數(shù)的估計(jì)值進(jìn)行后續(xù)計(jì)算。一個(gè)實(shí)施例中，如圖1所示，上述的黑油油藏模擬方法，還可包括步驟：S180：若所述的油藏壓力變化值及各流體相飽和度變化值不收斂，則用所述油藏壓力變化值和所述各流體相飽和度變化值作為相應(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ì)算的油藏壓力變化值及各流體相飽和度變化值收斂。本實(shí)施例中，在油藏壓力變化值及各流體相飽和度變化值不收斂的情況下，重復(fù)步驟S140～S160迭代計(jì)算，可以進(jìn)一步提高油藏壓力和各流體相飽和度的準(zhǔn)確度。本實(shí)施例的黑油油藏模擬方法，針對(duì)現(xiàn)有組分模擬包含組分相平衡及組分摩爾分?jǐn)?shù)的描述，而黑油模型的黑油模擬無(wú)組分相平衡及組分摩爾分?jǐn)?shù)的描述，導(dǎo)致的組分模擬與基于黑油模型的黑油模擬無(wú)法結(jié)合的問(wèn)題，獨(dú)創(chuàng)地建立描述黑油油藏流體滲流規(guī)律的組分模型，進(jìn)一步得到離散網(wǎng)格的組分?jǐn)?shù)值模型，奠定了用組分模型進(jìn)行黑油模擬的基礎(chǔ)，再利用流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表，為用組分模型進(jìn)行黑油模擬搭建了橋梁，實(shí)現(xiàn)了組分模擬和現(xiàn)有基于黑油模型黑油模擬的有效結(jié)合。一個(gè)實(shí)施例中，所述多個(gè)組分滲流參數(shù)可包括：油藏壓力、氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度及油組分在油相中摩爾濃度中的一個(gè)或多個(gè)。其中，油藏壓力可以用氣相壓力表示。該組分模型可以是包括多個(gè)按組分描述的方程，例如，可以包括守恒方程、烴組分相平衡方程、相內(nèi)部組分構(gòu)成的約束方程等。該組分模型為發(fā)明人通過(guò)創(chuàng)造性勞動(dòng)獨(dú)創(chuàng)得到。一個(gè)實(shí)施例中，該組分模型可描述1個(gè)氣組分、1個(gè)油組分和1個(gè)水組分的三相三組分流體在油藏中的運(yùn)移情況。該組分模型可包括3個(gè)組分守恒方程、2個(gè)烴組分相平衡方程(黑油逸度方程)及相內(nèi)部組分構(gòu)成的約束方程。上述各方程可用包含氣相壓力Pg，氣相飽和度Sg，油相飽和度So，氣組分在氣相中摩爾濃度Xgg，氣組分在油相中摩爾濃度Xgo，油組分在氣相中摩爾濃度Xog，油組分在油相中的摩爾濃度Xoo的7個(gè)變量描述，不同方程變量個(gè)數(shù)可不同。一個(gè)實(shí)施例中，所述多個(gè)組分滲流參數(shù)包括油藏壓力、氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度及油組分在油相中摩爾濃度，在上述步驟S110中，基于多個(gè)組分滲流參數(shù)建立描述黑油油藏流體滲流規(guī)律的組分模型的方法，具體實(shí)施方式可以是：基于油藏壓力、氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度及油組分在油相中摩爾濃度中的一個(gè)或多個(gè)，建立氣組分守恒方程、油組分守恒方程、水組分守恒方程、氣組分相平衡方程、油組分相平衡方程、氣相內(nèi)部組分約束方程及油相內(nèi)部組分約束方程，作為所述組分模型。該組分模型可包括：上述的氣組分守恒方程、油組分守恒方程、水組分守恒方程、氣組分相平衡方程、油組分相平衡方程、氣相內(nèi)部組分約束方程及油相內(nèi)部組分約束方程。一個(gè)實(shí)施例中，3個(gè)組分守恒方程可包括：氣組分守恒方程、油組分守恒方程及水組分守恒方程。所述氣組分守恒方程可為：所述油組分守恒方程可為：所述水組分守恒方程可為：一個(gè)實(shí)施例中，2個(gè)烴組分相平衡方程可包括：氣組分相平衡方程和油組分相平衡方程。在現(xiàn)有黑油模擬中，不存在關(guān)于相平衡的描述，而在現(xiàn)有的組分模擬中需要相平衡描述，所以，本實(shí)施例中通過(guò)獨(dú)創(chuàng)設(shè)計(jì)的相平衡方程能夠有效實(shí)現(xiàn)黑油模擬和組分模擬的結(jié)合。所述氣組分相平衡方程可為：Ffg(Xgg,Xgo)=XggPg-KgXogPg=Xgg-KgXgo=0,---(4)]]>所述油組分相平衡方程可為：Ffo(Xog,Xoo)=XogPg-KoXooPg=Xog-KoXoo=0.---(5)]]>一個(gè)實(shí)施例中，相內(nèi)部組分構(gòu)成的約束方程可包括：氣相內(nèi)部組分約束方程和油相內(nèi)部組分約束方程。本實(shí)施例中，通過(guò)各相內(nèi)部組分構(gòu)成的約束方程可以保證歸一化，從而保證黑油油藏模擬結(jié)果的正確性。所述氣相內(nèi)部組分約束方程可為：Fcg(Xog,Xgg)=Xog+Xgg-1=0,---(6)]]>所述油相內(nèi)部組分約束方程可為：Fco(Xoo,Xgo)=Xoo+Xgo-1=0.---(7)]]>在上述公式(1)～(7)中，其中，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ù)；Pg為氣相壓力(可間接表示油藏壓力)，Sg為氣相飽和度，So為油相飽和度，Xgg為氣組分在氣相中摩爾濃度，Xgo為氣組分在油相中摩爾濃度，Xog為油組分在氣相中摩爾濃度，Xoo為油組分在油相中的摩爾濃度；t為時(shí)間，為孔隙度，ρg為氣相密度，ρo為油相密度，K為滲透率，Krg為氣相相對(duì)滲透率，μg為氣相粘度，g為重力加速度，D為油藏深度(可為油藏地質(zhì)模型中離散網(wǎng)格中心所對(duì)應(yīng)的深度)，Kro為油相相對(duì)滲透率，μo為油相粘度，Pcgo為氣油兩相之間的毛管力，qg為氣相產(chǎn)量或注入量(例如水注入量)，qo為油相產(chǎn)量或注入量(例如水注入量)，ρw為水相密度，Krw為水相相對(duì)滲透率，μw為水相粘度，Pcow為油水兩相之間的毛管力，qw為水相產(chǎn)量或注入量。圖2是本發(fā)明一實(shí)施例中對(duì)組分模型進(jìn)行數(shù)值離散得到組分?jǐn)?shù)值模型的方法的流程示意圖。如圖2所示，在上述步驟S120中，根據(jù)基于黑油模型的黑油模擬所需數(shù)據(jù)對(duì)所述組分模型進(jìn)行數(shù)值離散，得到基于所述黑油油藏的地質(zhì)模型離散網(wǎng)格的組分?jǐn)?shù)值模型的方法，可包括步驟：S121：根據(jù)所述黑油模擬所需數(shù)據(jù)，基于所述黑油油藏的地質(zhì)模型離散網(wǎng)格，對(duì)偏微分方程形式的所述氣組分守恒方程、所述油組分守恒方程及所述水組分守恒方程進(jìn)行空間數(shù)值離散，得到相應(yīng)的非線性代數(shù)方程；S122：根據(jù)所述氣組分相平衡方程、所述油組分相平衡方程、所述氣相內(nèi)部組分約束方程、所述油相內(nèi)部組分約束方程及所有所述非線性代數(shù)方程，按牛頓迭代法計(jì)算得到包含雅可比矩陣單元的線性代數(shù)方程組，作為所述離散網(wǎng)格的組分?jǐn)?shù)值模型。在上述步驟S121中，可以采用多種方法對(duì)偏微分方程形式的所述氣組分守恒方程、所述油組分守恒方程及所述水組分守恒方程進(jìn)行空間數(shù)值離散，例如有限差分法、有限元法等。較佳地，是采用有限體積法對(duì)偏微分方程形式的所述氣組分守恒方程、所述油組分守恒方程及所述水組分守恒方程進(jìn)行空間數(shù)值離散，以此可以提高數(shù)值離散的精度。本實(shí)施例中，該組分?jǐn)?shù)值模型可以是包含雅可比矩陣單元的線性代數(shù)方程組可以簡(jiǎn)化計(jì)算。在其他實(shí)施例中，可以利用其他方法對(duì)所述氣組分守恒方程、所述油組分守恒方程及所述水組分守恒方程進(jìn)行空間數(shù)值離散。在其他實(shí)施例中，可以利用其他方法得到上述線性代數(shù)方程組。一個(gè)實(shí)施例中，通過(guò)為線性代數(shù)方程組選用適當(dāng)?shù)慕M分滲流參數(shù)以及安排多個(gè)組分滲流參數(shù)在線性代數(shù)方程組的變量矩陣中的順序，可使雅可比矩陣單元得到顯著簡(jiǎn)化，從而減少計(jì)算量。一個(gè)實(shí)施例中，可采用有限體積方法對(duì)組分模型中3個(gè)組分守恒方程(偏微分方程)進(jìn)行空間離散。例如，對(duì)上述公式(1)～(3)進(jìn)行空間離散。所述氣組分守恒方程所對(duì)應(yīng)的非線性代數(shù)方程可為：所述油組分守恒方程所對(duì)應(yīng)的非線性代數(shù)方程可為：所述水組分守恒方程所對(duì)應(yīng)的非線性代數(shù)方程可為：在上述公式(8)～(10)中，F(xiàn)g表示描述氣組分守恒關(guān)系的函數(shù)，F(xiàn)o表示描述油組分守恒關(guān)系的函數(shù)，F(xiàn)w表示描述水組分守恒關(guān)系的函數(shù)，Pg為氣相壓力，Sg為氣相飽和度，So為油相飽和度，Xgg為氣組分在氣相中摩爾濃度，Xgo為氣組分在油相中摩爾濃度，Xog為油組分在氣相中摩爾濃度，Xoo為油組分在油相中的摩爾濃度；Ω為空間積分變量，為孔隙度，ρg為氣相密度，ρo為油相密度，K為滲透率，Krg為氣相相對(duì)滲透率，μg為氣相粘度，g為重力加速度，D為油藏深度，Kro為油相相對(duì)滲透率，μo為油相粘度，Pcgo為氣油兩相之間的毛管力，qg為氣相產(chǎn)量或注入量，qo為油相產(chǎn)量或注入量，ρw為水相密度，Krw為水相相對(duì)滲透率，μw為水相粘度，Pcow為油水兩相之間的毛管力，qw為水相產(chǎn)量或注入量，∑l′表示對(duì)離散網(wǎng)格所有表面l′的通量求和，Al′表示網(wǎng)格一個(gè)表面l′的面積，V表示離散網(wǎng)格的體積，Qg＝V×qg表示離散網(wǎng)格的氣相注入或產(chǎn)出率，Qo＝V×qo表示離散網(wǎng)格的油相注入或產(chǎn)出率，Qw＝V×qw表示離散網(wǎng)格的水相注入或產(chǎn)出率。Ω表示對(duì)離散網(wǎng)格進(jìn)行空間積分的變量。其他符號(hào)的可與先前實(shí)施例中同樣符號(hào)具有相同的意義。在上述步驟S122中，所述氣組分相平衡方程、所述油組分相平衡方程、所述氣相內(nèi)部組分約束方程、所述油相內(nèi)部組分約束方程也可為非線性代數(shù)方程。根據(jù)所述氣組分相平衡方程、所述油組分相平衡方程、所述氣相內(nèi)部組分約束方程、所述油相內(nèi)部組分約束方程及所有所述非線性代數(shù)方程，按牛頓迭代法計(jì)算得到包含雅可比矩陣單元的線性代數(shù)方程組，形成離散網(wǎng)格的組分?jǐn)?shù)值模型。在其他實(shí)施例中，可以采用其他方法得到線性代數(shù)方程組。一個(gè)實(shí)施例中，按牛頓迭代法計(jì)算得到線性代數(shù)方程組可為：(∂F∂X)lδX=-Fl,---(11)]]>其中，F(xiàn)表示函數(shù)，X表示變量，一般情況下X根據(jù)組分相態(tài)不同變量取值可發(fā)生變化。δX為變化量，表示雅可比矩陣(方程組F對(duì)基本求解變量X的導(dǎo)數(shù))，l用于表示求解線性代數(shù)方程組的不同次數(shù)，例如l可表示迭代次數(shù)，其他符號(hào)的意義可與其他實(shí)施例中同樣符號(hào)的意義相同。本實(shí)施例中，變量X中選取的參數(shù)及各參數(shù)的排列順序可以大大簡(jiǎn)化雅可比矩陣?yán)?，雅克比矩陣單元的一般形式可表示為：如上述雅可比矩陣所示，采用本?shí)施例所建方程組得到的雅可比矩陣單元可為7行7列矩陣，可共有49個(gè)元素，其中，前17個(gè)元素需要計(jì)算得到，剩余元素中有2個(gè)需要查詢組分相平衡參數(shù)表(即參數(shù)-Kg和-Ko)，其余30個(gè)元素取值固定為0或者1，矩陣的這一特征有利于提高模擬速度。圖3是本發(fā)明一實(shí)施例中的利用流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表的方法示意圖。如圖3所示，在上述步驟S130中，利用所述黑油模擬所需數(shù)據(jù)中的流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表的方法，可包括步驟：S131：根據(jù)所述流體PVT數(shù)據(jù)中的油相密度、油相中溶解氣油比、氣相密度及氣相中揮發(fā)油氣比，確定油相中溶解氣組分與油組分的摩爾比和氣相中揮發(fā)油組分與氣組分的摩爾比；S132：根據(jù)所述的油相中溶解氣組分與油組分的摩爾比和所述的氣相中揮發(fā)油組分與氣組分的摩爾比，計(jì)算得到氣組分平衡常數(shù)和油組分平衡常數(shù)；S133：根據(jù)所述的氣組分平衡常數(shù)和油組分平衡常數(shù)，計(jì)算得到氣組分總摩爾分?jǐn)?shù)和油組分總摩爾分?jǐn)?shù)，所述組分相平衡參數(shù)表包括所述氣組分平衡常數(shù)、所述油組分平衡常數(shù)、所述氣組分總摩爾分?jǐn)?shù)、所述油組分總摩爾分?jǐn)?shù)及所述流體PVT數(shù)據(jù)中相應(yīng)的飽和壓力。本實(shí)施例中，發(fā)明人經(jīng)過(guò)創(chuàng)造性勞動(dòng)，通過(guò)所述流體PVT數(shù)據(jù)的油相中溶解氣油比和氣相中揮發(fā)油氣比計(jì)算得到溶解氣組分與油組分的摩爾比和揮發(fā)油組分與氣組分的摩爾比，進(jìn)而計(jì)算得到氣組分平衡常數(shù)和油組分平衡常數(shù)，再計(jì)算得到氣組分總摩爾分?jǐn)?shù)和油組分總摩爾分?jǐn)?shù)，最終實(shí)現(xiàn)組分相平衡參數(shù)表的建立。利用該組分相平衡參數(shù)表中組分總摩爾分?jǐn)?shù)、組分平衡常數(shù)及飽和壓力的對(duì)應(yīng)關(guān)系，能夠根據(jù)組分模擬中的組分總摩爾分?jǐn)?shù)參數(shù)(表示了組分構(gòu)成)數(shù)據(jù)查找得到黑油模擬中的飽和壓力參數(shù)數(shù)據(jù)和相應(yīng)的物性參數(shù)，從而實(shí)現(xiàn)組分模擬和黑油模擬的結(jié)合，即實(shí)現(xiàn)用組分模型實(shí)現(xiàn)黑油模擬。而且，現(xiàn)有黑油模擬僅能得到油藏壓力和油氣飽和度數(shù)據(jù)，本發(fā)明通過(guò)組分相平衡參數(shù)表還可以得到油氣組分總摩爾分?jǐn)?shù)，能夠滿足更多數(shù)據(jù)需求，具有更好的擴(kuò)展性。一個(gè)實(shí)施例中，在上述步驟S131中，根據(jù)油相中溶解氣油比Rgo(單位：SCF/STB)確定溶解氣組分與油組分的摩爾比rg，計(jì)算公式可如下所示：rg=(Rgo×ρgMrg)/(1×158.987×0.0353×ρoMro),---(13)]]>其中，Mrg為氣組分分子量，例如可為23；Mro為油組分分子量，例如可為190；分子量Mrg、Mro的取值不影響模擬結(jié)果。ρg為氣相質(zhì)量密度(SCF/STB)，ρo為油相質(zhì)量密度(SCF/STB)，其它參數(shù)為單位轉(zhuǎn)換的常數(shù)。在上述步驟S131中，根據(jù)氣相中揮發(fā)油氣比Rog(單位：SCF/STB)換算出揮發(fā)油組分與氣組分的摩爾比ro，計(jì)算公式可如下所示：ro=(Rog×158.987×0.0353ρoMro)/(1×ρgMrg).---(14)]]>一個(gè)實(shí)施例中，在上述步驟S132中，根據(jù)溶解氣組分與油組分的摩爾比rg和揮發(fā)油組分與氣組分的摩爾比ro，計(jì)算出氣組分的平衡常數(shù)Kg和油組分的平衡常數(shù)Ko，計(jì)算公式可如下所示：Kg=1+1/rg1+ro,---(15)]]>Ko=1+rg1+1/ro.---(16)]]>一個(gè)實(shí)施例中，在上述步驟S133中，根據(jù)氣組分的平衡常數(shù)Kg和油組分的平衡常數(shù)Ko，計(jì)算氣組分總的摩爾分?jǐn)?shù)Zg和油組分總的摩爾分?jǐn)?shù)Zo，計(jì)算公式可如下所示：Zg=1-KoKg-Ko,---(17)]]>Zo=Kg-1Kg-Ko.---(18)]]>圖4是本發(fā)明另一實(shí)施例的黑油油藏模擬方法的流程示意圖。如圖4所示，圖1所示的包括步驟S110～S170的黑油油藏模擬方法，在根據(jù)所述組分相平衡參數(shù)表和所述組分滲流參數(shù)的估計(jì)值進(jìn)行閃蒸計(jì)算，得到所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)(步驟S140)之前，還可包括步驟：S190：根據(jù)所述黑油模擬所需數(shù)據(jù)中的所述黑油油藏的地質(zhì)模型確定所述離散網(wǎng)格的各流體相飽和度初始值，利用靜水柱壓力公式計(jì)算所述離散網(wǎng)格的油藏壓力初始值，通過(guò)查詢所述組分相平衡參數(shù)表得到所述離散網(wǎng)格的各組總分摩爾分?jǐn)?shù)初始值，各流體相飽和度初始值、油藏壓力初始值及各組分總摩爾分?jǐn)?shù)初始值一同作為初始的所述多個(gè)組分滲流參數(shù)的估計(jì)值。在上述步驟S190中，上述地質(zhì)模型可以是以網(wǎng)格單元方式描述油藏幾何形狀、油氣水分布、孔隙度分布情況的數(shù)據(jù)體。地質(zhì)模型中的網(wǎng)格單元可以和上述離散網(wǎng)格相對(duì)應(yīng)。離散網(wǎng)格的各流體相飽和度可以包括氣相飽和度、油相飽和度等。油藏飽和度信息可直接利用地質(zhì)模型中飽和度分布信息確定。一個(gè)實(shí)施例中，可依據(jù)流體PVT數(shù)據(jù)中的多種不同數(shù)據(jù)，查詢組分相平衡參數(shù)表得到離散網(wǎng)格的各組分總摩爾分?jǐn)?shù)初始值。在其他實(shí)施例中，上述方法還可包括步驟S180。上述確定的油藏壓力、流體相飽和度和組分總摩爾分?jǐn)?shù)的初始估計(jì)值，在黑油模擬過(guò)程中可只作為迭代計(jì)算的初始值，后期隨著迭代進(jìn)行，油藏壓力、飽和度和組分總摩爾分?jǐn)?shù)的估計(jì)值可不斷調(diào)整。一些實(shí)施例中，油藏飽和度信息可利用前次迭代計(jì)算中的信息得到。一個(gè)實(shí)施例中，油藏壓力信息可利用靜水柱壓力公式計(jì)算得到，其中靜水柱壓力公式可為：P＝ρgh(19)其中，P為氣藏壓力，ρ為氣密度，g為重力加速度，h為高度差或油藏深度(例如距離海平面的高度差)。圖5是本發(fā)明一實(shí)施例中進(jìn)行閃蒸計(jì)算得到離散網(wǎng)格流體相態(tài)的方法流程示意圖。如圖5所示，在上述步驟S140中，根據(jù)所述組分相平衡參數(shù)表和所述組分滲流參數(shù)的估計(jì)值進(jìn)行閃蒸計(jì)算，得到所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)的方法，可包括步驟：S141：根據(jù)所述離散網(wǎng)格的所述組分滲流參數(shù)的當(dāng)前估計(jì)值計(jì)算得到所述離散網(wǎng)格的組分總摩爾分?jǐn)?shù)；S142：根據(jù)所述組分總摩爾分?jǐn)?shù)查詢所述組分相平衡參數(shù)表，得到所述離散網(wǎng)格中組分構(gòu)成所對(duì)應(yīng)的飽和壓力；S143：通過(guò)對(duì)比所述離散網(wǎng)格的飽和壓力和所述組分滲流參數(shù)中的油藏壓力的當(dāng)前估計(jì)值，確定所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)。在上述步驟S141中，在第一次迭代過(guò)程中，組分滲流參數(shù)的當(dāng)前估計(jì)值可以是上述初始的估計(jì)值。上述離散網(wǎng)格的組分的總摩爾分?jǐn)?shù)可以是氣組分的總摩爾分?jǐn)?shù)或油組分的總摩爾分?jǐn)?shù)。在上述步驟S142中，可根據(jù)氣組分的總摩爾分?jǐn)?shù)或油組分的總摩爾分?jǐn)?shù)查詢所述組分相平衡參數(shù)表，得到所述離散網(wǎng)格的組分構(gòu)成所對(duì)應(yīng)的飽和壓力。一個(gè)實(shí)施例中，根據(jù)各組分的總摩爾分?jǐn)?shù)(表示了離散網(wǎng)格的組分構(gòu)成)，可通過(guò)線性插值法查詢組分相平衡參數(shù)表，得到所述離散網(wǎng)格的組分構(gòu)成所對(duì)應(yīng)的飽和壓力，以此可以更準(zhǔn)確地查詢得到離散網(wǎng)格的組分構(gòu)成情況下的飽和壓力。在其他實(shí)施例中，可以采用其他方法查詢組分相平衡參數(shù)表得到離散網(wǎng)格的飽和壓力。一個(gè)實(shí)施例中，所述離散網(wǎng)格的組分可包括油組分和氣組分，所述組分滲流參數(shù)的當(dāng)前估計(jì)值可包括：離散網(wǎng)格當(dāng)前的氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度及油組分在油相中的摩爾濃度?？筛鶕?jù)當(dāng)前網(wǎng)格氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度、油組分在油相中的摩爾濃度，計(jì)算出油組分總的摩爾分?jǐn)?shù)和氣組分總的摩爾分?jǐn)?shù)。一個(gè)實(shí)施例中，計(jì)算油組分總的摩爾分?jǐn)?shù)和氣組分總的摩爾分?jǐn)?shù)公式可如下所示：油組分的總摩爾分?jǐn)?shù)可為：Zo=(SgρgXog+SoρoXoo)(SgρgXog+SoρoXoo+SgρgXgg+SoρoXgo),---(20)]]>氣組分的總摩爾分?jǐn)?shù)可為：Zg=(SgρgXgg+SoρoXgo)(SgρgXog+SoρoXoo+SgρgXgg+SoρoXgo)=1-Zo,---(21)]]>其中，Zo為油組分的總摩爾分?jǐn)?shù)，Zg為氣組分的總摩爾分?jǐn)?shù)，Sg為氣相飽和度，So為油相飽和度，ρg為氣相密度，ρo為油相密度，Xog為油組分在氣相中摩爾濃度，Xoo為油組分在油相中的摩爾濃度，Xgg為氣組分在氣相中摩爾濃度，Xgo為氣組分在油相中摩爾濃度。本實(shí)施例中，獨(dú)創(chuàng)地利用組分滲流參數(shù)的估計(jì)值計(jì)算出離散網(wǎng)格內(nèi)各組分的總摩爾分?jǐn)?shù)，有利于準(zhǔn)確地從組分相平衡參數(shù)表中查找到飽和壓力。圖6是本發(fā)明一實(shí)施例中通過(guò)對(duì)比飽和壓力和油藏壓力的估計(jì)值確定離散網(wǎng)格流體相態(tài)的方法流程示意圖。如圖6所示，在上述步驟S143中，通過(guò)對(duì)比所述離散網(wǎng)格的飽和壓力和所述組分滲流參數(shù)中的油藏壓力的當(dāng)前估計(jì)值，確定所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)的方法，可包括步驟：S1431：根據(jù)所述離散網(wǎng)格的油藏壓力的當(dāng)前估計(jì)值大于所述飽和壓力判定所述離散網(wǎng)格的流體相態(tài)為油相；S1432：當(dāng)所述離散網(wǎng)格的油藏壓力的當(dāng)前估計(jì)值小于所述飽和壓力時(shí)，根據(jù)所述離散網(wǎng)格的油藏壓力的當(dāng)前估計(jì)值查詢所述組分相平衡參數(shù)表，得到所述離散網(wǎng)格的氣組分平衡常數(shù)和油組分平衡常數(shù)，并根據(jù)所述的氣組分平衡常數(shù)和油組分平衡常數(shù)計(jì)算所述離散網(wǎng)格的液相比例；S1433：根據(jù)所述液相比例大于零判定所述離散網(wǎng)格的流體相態(tài)為油氣兩相共存，根據(jù)所述液相比例等于或小于零判定所述離散網(wǎng)格的流體相態(tài)為氣相；S1434：根據(jù)所述離散網(wǎng)格的流體相態(tài)的判定結(jié)果和油藏壓力的當(dāng)前估計(jì)值，查詢所述流體PVT數(shù)據(jù)得到相應(yīng)物性參數(shù)。在上述步驟S1431中，當(dāng)離散網(wǎng)格的油藏壓力當(dāng)前估計(jì)值大于飽和壓力時(shí)，即(Pg-Pcgo＞Pb)，網(wǎng)格內(nèi)部不存在氣體，所有烴類(lèi)組分都存在于油相中，例如所有氣組分都溶解在油相中，可以認(rèn)為離散網(wǎng)格的流體相態(tài)為油相。若判斷已離散網(wǎng)格的油藏壓力的當(dāng)前估計(jì)值大于飽和壓力Pb，則不需作進(jìn)一步判斷。油藏壓力是指離散網(wǎng)格中流體相壓力。當(dāng)離散網(wǎng)格中油氣水三相共存時(shí)，油藏壓力有三個(gè)值，他們之間相差一個(gè)毛管力(只與飽和度有關(guān))，知道一個(gè)流體相壓力和流體飽和度，其它相的壓力可以換算得到。對(duì)一個(gè)流體相而言，油藏壓力和飽和壓力統(tǒng)一，例如對(duì)于油相，提及“當(dāng)油藏壓力大于飽和壓力時(shí)，網(wǎng)格不存在氣體”，這里的油藏壓力就可是油相壓力，可以通過(guò)氣相壓力加上一個(gè)毛管力換算處理，這里的飽和壓力可以是油相PVT表中的飽和壓力。在上述步驟S1432中，當(dāng)離散網(wǎng)格的油藏壓力當(dāng)前估計(jì)值小于飽和壓時(shí)，需要進(jìn)一步確定離散網(wǎng)格的流體相態(tài)是油氣兩相還是氣相，例如可以利用離散網(wǎng)格的液相比例作進(jìn)一步判斷。在上述步驟S1434中，若所述流體相態(tài)的判定結(jié)果為氣相，根據(jù)所述離散網(wǎng)格的流體相態(tài)的判定結(jié)果和油藏壓力的當(dāng)前估計(jì)值，查詢所述流體PVT數(shù)據(jù)中氣相的PVT數(shù)據(jù)，可得到相應(yīng)物性參數(shù)，所述相應(yīng)物性參數(shù)可包括氣相密度和氣體積系數(shù)。若所述流體相態(tài)的判定結(jié)果為油相，根據(jù)所述離散網(wǎng)格的流體相態(tài)的判定結(jié)果和油藏壓力的當(dāng)前估計(jì)值，查詢所述流體PVT數(shù)據(jù)中油相的PVT數(shù)據(jù)可得到相應(yīng)物性參數(shù)，所述相應(yīng)物性參數(shù)可包括油相密度和油相體積系數(shù)。一個(gè)實(shí)施例中，上述步驟S1432中，該液相比例可為：V′=Zg1-Ko+Zo1-Kg,---(22)]]>其中，V'表示液相比例，Zg為氣組分的總摩爾分?jǐn)?shù)，Zo為油組分的總摩爾分?jǐn)?shù)，Ko油組分平衡常數(shù)，Kg為氣組分平衡常數(shù)。一個(gè)實(shí)施例中，如果液相比例V'＞0，則油氣兩相共存，油氣組分按照平衡參數(shù)取值在液氣兩相中分配；如果液相比例V'＝0，則只有氣相存在，所有組分都只存在氣相中。一個(gè)實(shí)施例中，在上述步驟S150中，根據(jù)所述流體相態(tài)和相應(yīng)物性參數(shù)調(diào)整所述組分?jǐn)?shù)值模型。根據(jù)所述流體相態(tài)調(diào)整組分?jǐn)?shù)值模型包括三種情況，在調(diào)整后的所述組分?jǐn)?shù)值模型中：當(dāng)所述離散網(wǎng)格的流體相態(tài)為油氣兩相共存時(shí)，變量可保持上述變量的形式，雅可比矩陣單元可采用上述雅可比矩陣的一般形式(式(12))：當(dāng)所述離散網(wǎng)格的流體相態(tài)僅為油相，即網(wǎng)格內(nèi)僅存在油相時(shí)，上述變量X可調(diào)整為變量所述雅可比矩陣單元可調(diào)整為：當(dāng)所述離散網(wǎng)格的流體相態(tài)僅為氣相，即網(wǎng)格內(nèi)僅存在氣相時(shí)，上述變量X可調(diào)整為變量所述雅可比矩陣單元可調(diào)整為：一個(gè)實(shí)施例中，在確定雅可比系數(shù)矩陣和右端項(xiàng)-Fl的情況下，利用線性求解算法計(jì)算如下方程組，獲得δX的值。(∂F∂X)lδX=-Fl.---(25)]]>一個(gè)實(shí)施例中，根據(jù)流體相態(tài)的不同，變化量δX可分三種情況：情況1：當(dāng)網(wǎng)格里面油氣兩相共存時(shí)，求解變化量取前三項(xiàng)值，即氣相壓力/油藏壓力變化值δPg、氣相飽和度變化值δSg、油相飽和度變化值δSo，分別與油藏壓力估計(jì)值、氣相飽和度估計(jì)值、油相飽和度估計(jì)值相加，即可得到油藏壓力及油、氣飽和度等信息。情況2：當(dāng)網(wǎng)格里面只存在油相時(shí)，求解變化量為取前2項(xiàng)值，即氣相壓力/油藏壓力變化值δPg和油相飽和度變化值δSo，分別與油藏壓力估計(jì)值和油相飽和度估計(jì)值相加，即可得到油藏壓力和油相飽和度信息。情況3：當(dāng)網(wǎng)格里面只存在氣相時(shí)，求解變化量為取前2項(xiàng)值，即氣相壓力/油藏壓力變化值δPg和油相飽和度變化值δSg，分別與油藏壓力估計(jì)值和氣相飽和度估計(jì)值相加，即可得到油藏壓力和氣飽和度信息。本實(shí)施例中，由上述各變化量δX的求解結(jié)果可以看出，調(diào)整變量X后，組分?jǐn)?shù)值模型得到簡(jiǎn)化，計(jì)算過(guò)程也得到簡(jiǎn)化，由此可提高黑油模擬的效率。一個(gè)實(shí)施例中，在上述步驟S170中，判斷所述的油藏壓力變化值及各流體相飽和度變化值是否收斂，具體實(shí)施方式可為：根據(jù)變化值δXl和相應(yīng)估計(jì)值Xl是否滿足判定所述的油藏壓力變化值及各流體相飽和度變化值收斂，其中，變化值δXl表示第l次迭代后油藏壓力變化值或各流體相飽和度變化值，相應(yīng)估計(jì)值Xl表示前次迭代后油藏壓力的估計(jì)值或各流體相飽和度的估計(jì)值，迭代次數(shù)l≥1，X1為初始的估計(jì)值，a為設(shè)定值，例如0＜a＜0.1。較佳實(shí)施例中，設(shè)定值a＝0.01，當(dāng)變化值δXl與估計(jì)值Xl之比小于一設(shè)定值，即判斷收斂，則可進(jìn)行模擬結(jié)果輸出處理；否則認(rèn)為不收斂，更新估計(jì)值，更新后的估計(jì)值Xl+1＝Xl+δXl，重復(fù)步驟S140～S160計(jì)算直到變化值收斂，將最新的變化值和前次計(jì)算的估計(jì)值相加作為最終估計(jì)值。此處，估計(jì)值可以是油藏壓力估計(jì)值、油相飽和度估計(jì)值等。一個(gè)實(shí)施例中，根據(jù)組分?jǐn)?shù)值模型求解得到的油藏壓力和飽和度信息，根據(jù)流體相態(tài)不同分三種情況獲得水相飽和度信息。例如：情況1：油氣兩相共存時(shí)，水相飽和度：Sw＝1-Sg-So情況2：只有油相存在時(shí)，水相飽和度：Sw＝1-So情況3：只有氣相存在時(shí)，水相飽和度：Sw＝1-Sg圖7是本發(fā)明又一實(shí)施例的黑油模擬方法的流程示意圖。如圖7所示，采集黑油模擬所需的流體PVT、相滲、毛管力數(shù)據(jù)、地質(zhì)模型等數(shù)據(jù)(步驟S301)；根據(jù)采集的流體PVT數(shù)據(jù)，確定組分相平衡參數(shù)表(步驟S302)；建立一個(gè)描述黑油油藏流體滲流規(guī)律、包含油藏壓力、油氣飽和度和油氣組分總摩爾分?jǐn)?shù)等信息的組分模型(步驟S303)；根據(jù)采集的數(shù)據(jù)對(duì)組分模型進(jìn)行數(shù)值離散，確定油藏壓力、飽和度和組分總摩爾分?jǐn)?shù)的估計(jì)值，并形成初始的組分?jǐn)?shù)值模型(步驟S304)；根據(jù)組分相平衡參數(shù)表和油藏估計(jì)值進(jìn)行閃蒸計(jì)算，確定流體相態(tài)和參數(shù)(步驟S305)；根據(jù)閃蒸計(jì)算得到的流體相態(tài)，調(diào)整組分?jǐn)?shù)值模型(步驟S306)；根據(jù)組分?jǐn)?shù)值模型，線性求解得到油藏壓力、油氣飽和度信息(步驟S307)；判斷計(jì)算結(jié)果是否收斂(步驟S308)；最后，輸出模擬結(jié)果(步驟S309)。圖8是本發(fā)明另一實(shí)施例中進(jìn)行閃蒸計(jì)算得到離散網(wǎng)格流體相態(tài)的方法流程示意圖。如圖8所示，根據(jù)當(dāng)前網(wǎng)格氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度、油組分在油相中的摩爾濃度等信息，計(jì)算出油組分總的摩爾分?jǐn)?shù)和氣組分總的摩爾分?jǐn)?shù)(步驟S401)；根據(jù)當(dāng)前網(wǎng)格油組分總摩爾分?jǐn)?shù)(或者氣組分總摩爾分?jǐn)?shù))，通過(guò)線性插值方式查詢組分平衡參數(shù)表得到該組分構(gòu)成情況下對(duì)應(yīng)的飽和壓力(步驟S402)。對(duì)比網(wǎng)格壓力和飽和壓力，確定網(wǎng)格內(nèi)部流體相態(tài)(步驟S403)，判斷方法如下：當(dāng)網(wǎng)格壓力大于飽和壓力時(shí)，網(wǎng)格內(nèi)部不存在氣體，所有烴類(lèi)組分都存在于油相中；當(dāng)網(wǎng)格壓力小于飽和壓力時(shí)，進(jìn)入步驟S404計(jì)算液相比例；據(jù)當(dāng)前網(wǎng)格壓力查組分平衡參數(shù)表得到氣組分的平衡常數(shù)和油組分的平衡常數(shù)，計(jì)算液相比例(公式(22))，如果液相比例大于0,則油氣兩相共存，油氣組分按照平衡參數(shù)取值在液氣兩相中分配；如果液相比例等于0，則只有氣相存在，所有組分都只存在氣相中。一個(gè)具體實(shí)施例中，某個(gè)油藏采用注水開(kāi)發(fā)方式，利用本發(fā)明實(shí)施例的采用組分模型實(shí)現(xiàn)黑油模擬的方法預(yù)測(cè)油田開(kāi)發(fā)一定周期后剩余油分布情況。具體步驟如下：步驟1：采集黑油模擬所需的流體PVT、相滲、毛管力數(shù)據(jù)、地質(zhì)模型等數(shù)據(jù)。地質(zhì)模型在X、Y和Z三個(gè)方向的網(wǎng)格數(shù)分別可以是10X10X3，X和Y方向的網(wǎng)格尺寸可均為300米，Z方向三層的網(wǎng)格尺度可分別是6米、9米和15米，孔隙度可為0.28、滲透率按層給出可以分別是500、50和200毫達(dá)西。油藏頂面深度可為2537.5米。油藏油氣界面深度可為2376米，油水接觸面深度可為2756.6米。相滲數(shù)據(jù)描述了不同流體飽和度條件下相對(duì)滲透率變化，由于數(shù)據(jù)量較大，只給出部分?jǐn)?shù)據(jù)。其中，油水兩相的相對(duì)滲透率和毛管力數(shù)據(jù)如表1所示。油氣兩相相對(duì)滲透率和毛管力數(shù)據(jù)如表2所示。氣相流體標(biāo)準(zhǔn)狀況下密度可為0.0702，因數(shù)據(jù)量較大，下面只列出部分PVT表數(shù)據(jù)，如表3所示。油相流體標(biāo)準(zhǔn)狀況下密度可為44.98，因數(shù)據(jù)量較大，下面只列出部分PVT表數(shù)據(jù)，如表4所示。表1油水兩相的相對(duì)滲透率和毛管力數(shù)據(jù)表2油氣兩相相對(duì)滲透率和毛管力數(shù)據(jù)表3氣相流體標(biāo)準(zhǔn)狀況下的PVT數(shù)據(jù)壓力體積系數(shù)粘度油氣比14.71.00001.200.0400.01.01201.17165800.01.02551.14335………4000.01.12000.941500表4油相流體標(biāo)準(zhǔn)狀況下的PVT數(shù)據(jù)油相流體標(biāo)準(zhǔn)狀況下密度可為63.01，考慮到水的性質(zhì)相對(duì)簡(jiǎn)單，且不參與閃蒸計(jì)算，所以可只采集參考?jí)毫ο滦再|(zhì)即可，如表5所示。壓力體積系數(shù)粘度油氣比3600.01.00340.960.0表5水相流體標(biāo)準(zhǔn)狀況下的PVT數(shù)據(jù)一個(gè)實(shí)施例中，針對(duì)油藏模型的第三層進(jìn)行研究，在地質(zhì)模型對(duì)角位置分別布置兩口定壓生產(chǎn)井。其中，一口注水井定壓可為8000Psia，一口生產(chǎn)井定壓可為1000Psia。步驟2：根據(jù)采集的流體PVT數(shù)據(jù)，確定組分相平衡參數(shù)表?？筛鶕?jù)公式(13)～(18)，計(jì)算得到組分相平衡參數(shù)表，可如表6所示：表6組分相平衡參數(shù)表步驟3：建立一個(gè)描述黑油油藏流體滲流規(guī)律、包含油藏壓力、油氣飽和度和油氣組分總摩爾分?jǐn)?shù)等信息的組分模型。所建組分模型如方程可如公式(1)～(7)所示。步驟4：根據(jù)采集的數(shù)據(jù)對(duì)組分模型進(jìn)行數(shù)值離散，確定油藏壓力、飽和度和組分總摩爾分?jǐn)?shù)的估計(jì)值，并形成初始的組分?jǐn)?shù)值模型。可采用有限體積方法離散得到的非線性代數(shù)方程如方程，公式(8)～(10)和公式(4)～(7)所示。油藏壓力可根據(jù)公式(19)計(jì)算得到，從上往下油藏壓力逐步增大，整個(gè)油藏平均壓力值可為4026.3PSia。油藏飽和度根據(jù)油氣、油水接觸面深度給定，所有網(wǎng)格初始含氣飽和度可為0，含水飽和度為束縛水飽和度可為0.18，含油飽和度可為0.88?？筛鶕?jù)網(wǎng)格壓力，查詢組分相平衡參數(shù)表得到組分總摩爾分?jǐn)?shù)的估計(jì)值。例如，網(wǎng)格壓力為4000PSia，查組分相平衡參數(shù)表可得到氣相總摩爾分?jǐn)?shù)為0.775，油相總摩爾分?jǐn)?shù)為0.225。離散后形成的雅可比矩陣及數(shù)值模型數(shù)據(jù)量非常巨大，這里列出其中部分?jǐn)?shù)據(jù)，其中，模型中第一個(gè)網(wǎng)格存在油氣兩相，其雅可比矩陣中數(shù)據(jù)可如下所示：模型中第2個(gè)網(wǎng)格離散后其雅可比矩陣中數(shù)據(jù)如下：步驟5：根據(jù)組分相平衡參數(shù)表和油藏估計(jì)值進(jìn)行閃蒸計(jì)算，確定流體相態(tài)和參數(shù)?？蓪?duì)于模型中第2個(gè)網(wǎng)格，根據(jù)公式(20)和(21)計(jì)算得到氣組分總的摩爾分?jǐn)?shù)為0.70，油組分總的摩爾分?jǐn)?shù)為0.3。查組分相平衡參數(shù)表可知該組分構(gòu)成條件下，飽和壓力3880psia。當(dāng)前網(wǎng)格壓力為4020Psia，大于查表得到的飽和壓力，故該網(wǎng)格相態(tài)為純油相，氣組分和油組分都只存在于油相中。對(duì)于其它網(wǎng)格，可利用同樣的方法判斷其流體相態(tài)。步驟6：根據(jù)閃蒸計(jì)算得到的流體相態(tài)，調(diào)整組分?jǐn)?shù)值模型。因數(shù)據(jù)規(guī)模巨大，按照方法說(shuō)明所述，不失一般性，仍以第2個(gè)網(wǎng)格為例說(shuō)明使用細(xì)節(jié)。對(duì)于第一個(gè)網(wǎng)格，在閃蒸計(jì)算中判斷為油相，求解變量調(diào)整為相應(yīng)的則雅可比矩陣調(diào)整可為：步驟7：根據(jù)組分?jǐn)?shù)值模型，線性求解得到油藏壓力、油氣飽和度信息。采用廣義最小余量方法求解組分?jǐn)?shù)值模型，得到網(wǎng)格2的解壓力相對(duì)變化值(8.58/4020)小于0.01，認(rèn)為收斂。則得到該網(wǎng)格壓力值為：4020-8.58＝4011.42Psia，網(wǎng)格氣相飽和度為：0.0+0.0＝0.0，網(wǎng)格油相飽和度為：0.88-8.4617E(-7)≈0.8799，網(wǎng)格水相飽和度為：1-0.0–0.8799＝0.2201。圖9是利用本發(fā)明實(shí)施例的黑油模擬方法得到的油相飽和度示意圖，如圖9所示，利用本發(fā)明實(shí)施例的方法，可得到某個(gè)時(shí)刻油藏油相飽和度分布情況，模擬方法有效。本發(fā)明實(shí)施例的黑油模擬方法，采集黑油模擬所需的數(shù)據(jù)；根據(jù)采集的流體PVT數(shù)據(jù)，建立組分相平衡參數(shù)表；建立一個(gè)描述黑油油藏流體滲流規(guī)律、包含油藏壓力、油氣飽和度和油氣組分摩爾分?jǐn)?shù)等信息的組分模型；根據(jù)采集的數(shù)據(jù)對(duì)組分模型進(jìn)行數(shù)值離散，確定油藏壓力、飽和度和組分摩爾分?jǐn)?shù)的估計(jì)值，并形成初始的組分?jǐn)?shù)值模型；根據(jù)相平衡參數(shù)和油藏估計(jì)值進(jìn)行閃蒸計(jì)算，確定流體相態(tài)和構(gòu)成；根據(jù)閃蒸計(jì)算得到的流體相態(tài)，調(diào)整組分?jǐn)?shù)值模型；根據(jù)組分?jǐn)?shù)值模型，線性求解得到油藏壓力、油氣飽和度信息；判斷油藏壓力和飽和度信息是否滿足收斂要求，收斂則輸出模擬結(jié)果，不收斂，調(diào)整估計(jì)值，重新計(jì)算。本發(fā)明實(shí)施例的黑油模擬方法，利用組分模型和模擬器架構(gòu)實(shí)現(xiàn)了黑油模擬功能，提高了模擬器集成性和后期擴(kuò)展性；直接利用現(xiàn)有黑油模擬的輸入數(shù)據(jù)，不需要額外給出組分相平衡信息，方便用戶使用；利用組分相平衡參數(shù)表的閃蒸計(jì)算，保證了模擬效率。本實(shí)施例的黑油油藏模擬方法獨(dú)創(chuàng)地建立描述黑油油藏流體滲流規(guī)律的組分模型，進(jìn)一步得到離散網(wǎng)格的組分?jǐn)?shù)值模型，奠定了用組分模型進(jìn)行黑油模擬的基礎(chǔ)，再利用流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表，為用組分模型進(jìn)行黑油模擬搭建了橋梁，實(shí)現(xiàn)了組分模擬和現(xiàn)有基于黑油模型黑油模擬的有效結(jié)合?；谂c圖1所示的油藏?cái)?shù)值模擬方法相同的發(fā)明構(gòu)思，本申請(qǐng)實(shí)施例還提供了一種油藏?cái)?shù)值模擬裝置，如下面實(shí)施例所述。由于該油藏?cái)?shù)值模擬裝置解決問(wèn)題的原理與油藏?cái)?shù)值模擬方法相似，因此該油藏?cái)?shù)值模擬裝置的實(shí)施可以參見(jiàn)油藏?cái)?shù)值模擬方法的實(shí)施，重復(fù)之處不再贅述。圖10是本發(fā)明一實(shí)施例的黑油油藏模擬裝置的結(jié)構(gòu)示意圖。如圖10所示，本發(fā)明實(shí)施例的黑油油藏模擬裝置，可包括：組分模型建立單元210、組分?jǐn)?shù)值模型生成單元220、組分相平衡參數(shù)表建立單元230、流體相態(tài)確定單元240、組分?jǐn)?shù)值模型調(diào)整單元250、油藏參數(shù)變化值生成單元260及油藏參數(shù)確定單元270，上述各模塊順序連接。組分模型建立單元210用于執(zhí)行：基于多個(gè)組分滲流參數(shù)建立描述黑油油藏流體滲流規(guī)律的組分模型。組分?jǐn)?shù)值模型生成單元220用于執(zhí)行：根據(jù)基于黑油模型的黑油模擬所需數(shù)據(jù)對(duì)所述組分模型進(jìn)行數(shù)值離散，得到基于所述黑油油藏的地質(zhì)模型離散網(wǎng)格的組分?jǐn)?shù)值模型。組分相平衡參數(shù)表建立單元230用于執(zhí)行：利用所述黑油模擬所需數(shù)據(jù)中的流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表。流體相態(tài)確定單元240用于執(zhí)行：根據(jù)所述組分相平衡參數(shù)表和所述組分滲流參數(shù)的估計(jì)值進(jìn)行閃蒸計(jì)算，得到所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)。組分?jǐn)?shù)值模型調(diào)整單元250用于執(zhí)行：根據(jù)所述流體相態(tài)和相應(yīng)物性參數(shù)調(diào)整所述組分?jǐn)?shù)值模型，以簡(jiǎn)化組分?jǐn)?shù)值模型。油藏參數(shù)變化值生成單元260用于執(zhí)行：利用所述組分滲流參數(shù)的估計(jì)值和調(diào)整后的所述組分?jǐn)?shù)值模型，計(jì)算得到油藏壓力變化值及各流體相飽和度變化值。油藏參數(shù)確定單元270用于執(zhí)行：判斷所述的油藏壓力變化值及各流體相飽和度變化值是否收斂，若是，根據(jù)所述油藏壓力變化值和相應(yīng)的所述估計(jì)值求和得到所述離散網(wǎng)格的油藏壓力，根據(jù)所述各流體相飽和度變化值和相應(yīng)的所述估計(jì)值求和得到所述離散網(wǎng)格的各流體相飽和度。組分滲流參數(shù)可以是油藏壓力、油氣飽和度和油氣組分摩爾分?jǐn)?shù)等參數(shù)。不同組分可以是氣組分、油組分和水組分中的一個(gè)或多個(gè)。上述基于黑油模型的黑油模擬是現(xiàn)有的黑油模擬方法?；诤谟湍Ｐ偷暮谟湍M所需數(shù)據(jù)可包括：流體PVT數(shù)據(jù)、相滲數(shù)據(jù)、毛管力數(shù)據(jù)、地質(zhì)模型等。其中，流體PVT數(shù)據(jù)可包括：油相流體PVT數(shù)據(jù)、氣相流體PVT數(shù)據(jù)、油藏水相流體PVT數(shù)據(jù)。上述流體PVT數(shù)據(jù)可從實(shí)驗(yàn)室測(cè)量得到。上述相滲數(shù)據(jù)可以是相對(duì)滲透率曲線，上述毛管力數(shù)據(jù)可以是毛管力曲線，上述相滲數(shù)據(jù)和毛管力數(shù)據(jù)可同試驗(yàn)室測(cè)量得到。上述地質(zhì)模型可通過(guò)油藏地質(zhì)建模系統(tǒng)生成。其中，上述油相流體PVT數(shù)據(jù)可包括：標(biāo)準(zhǔn)狀態(tài)下油相密度、不同飽和壓力取值下對(duì)應(yīng)的溶解氣油比、粘度及體積系數(shù)；上述氣相流體PVT數(shù)據(jù)可包括標(biāo)準(zhǔn)狀態(tài)下氣相密度、不同飽和壓力取值下對(duì)應(yīng)的揮發(fā)油氣比、粘度及體積系數(shù)；上述水相流體PVT數(shù)據(jù)可包括標(biāo)準(zhǔn)狀態(tài)下水相密度、不同壓力取值下對(duì)應(yīng)的粘度及體積系數(shù)。上述相對(duì)滲透率曲線可包括：油水兩相通過(guò)油藏巖心的相對(duì)滲透率、油氣兩相通過(guò)油藏巖心的相對(duì)滲透率。上述毛管力曲線可包括：不同含水飽和度下油水兩相毛管力、不同含水飽和度下油氣兩相毛管力。上述地質(zhì)模型可以是以網(wǎng)格單元方式描述油藏幾何形狀、油氣水分布、孔隙度分布情況的數(shù)據(jù)體?？梢圆捎枚喾N方法對(duì)組分模型進(jìn)行數(shù)值離散，數(shù)值離散后的組分模型可為非線性代數(shù)方程，不同的離散網(wǎng)格可對(duì)應(yīng)其各自的非線性代數(shù)方程。利用所述黑油模擬所需數(shù)據(jù)中的流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表，可以將流體PVT數(shù)據(jù)表轉(zhuǎn)換為組分相平衡參數(shù)表，組分相平衡參數(shù)表中的數(shù)據(jù)可與流體PVT數(shù)據(jù)具有一定的對(duì)應(yīng)關(guān)系，所以依據(jù)該組分相平衡參數(shù)表中的數(shù)據(jù)可以查詢得到流體PVT數(shù)據(jù)中的數(shù)據(jù)。該物性參數(shù)可以包括流體粘度等參數(shù)。組分滲流參數(shù)的估計(jì)值可以是初始設(shè)定的估計(jì)值，或者可以是根據(jù)黑油模擬計(jì)算結(jié)果計(jì)算得到，用于下次重新進(jìn)行黑油模擬。黑油模擬的結(jié)果通常包括油藏壓力和飽和度信息，根據(jù)本發(fā)明的方法還可以得到其他多種信息。該各流體相飽和度變化值可以包括氣相飽和度變化值、油相飽和度變化值及水相飽和度變化值。離散網(wǎng)格的各流體相飽和度可包括氣相飽和度、油相飽和度及水相飽和度，三者之和可為1，水相飽和度可根據(jù)氣相飽和度和油相飽和度計(jì)算得到。通過(guò)判斷油藏壓力變化值及各流體相飽和度變化值是否收斂，并根據(jù)收斂的油藏壓力變化值及各流體相飽和度變化值計(jì)算得到黑油模擬結(jié)果，以此可以保證黑油模擬結(jié)果的準(zhǔn)確度。一個(gè)實(shí)施例中，如圖10所示，上述的黑油油藏模擬裝置，還可包括：迭代計(jì)算單元280，連接于流體相態(tài)確定單元240和油藏參數(shù)確定單元270之間。迭代計(jì)算單元280用于執(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ì)算的油藏壓力變化值及各流體相飽和度變化值收斂。本實(shí)施例中，在油藏壓力變化值及各流體相飽和度變化值不收斂的情況下，重復(fù)迭代計(jì)算，可以進(jìn)一步提高油藏壓力和各流體相飽和度的準(zhǔn)確度。本實(shí)施例的黑油油藏模擬裝置，針對(duì)現(xiàn)有組分模擬包含相平衡及摩爾分?jǐn)?shù)的描述，而黑油模型的黑油模擬無(wú)相平衡及摩爾分?jǐn)?shù)的描述，導(dǎo)致的組分模擬與基于黑油模型的黑油模擬無(wú)法結(jié)合的問(wèn)題，獨(dú)創(chuàng)地建立描述黑油油藏流體滲流規(guī)律的組分模型，進(jìn)一步得到離散網(wǎng)格的組分?jǐn)?shù)值模型，奠定了用組分模型進(jìn)行黑油模擬的基礎(chǔ)，再利用流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表，為用組分模型進(jìn)行黑油模擬搭建了橋梁，實(shí)現(xiàn)了組分模擬和現(xiàn)有基于黑油模型黑油模擬的有效結(jié)合。一個(gè)實(shí)施例中，所述多個(gè)組分滲流參數(shù)包括油藏壓力、氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度及油組分在油相中摩爾濃度，所述組分模型建立單元210，可包括：組分模型建立模塊211。組分模型建立模塊211用于執(zhí)行：基于油藏壓力、氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度及油組分在油相中摩爾濃度中的一個(gè)或多個(gè)，建立氣組分守恒方程、油組分守恒方程、水組分守恒方程、氣組分相平衡方程、油組分相平衡方程、氣相內(nèi)部組分約束方程及油相內(nèi)部組分約束方程，作為所述組分模型。一個(gè)實(shí)施例中：所述氣組分守恒方程可為：所述油組分守恒方程可為：所述水組分守恒方程可為：一個(gè)實(shí)施例中，2個(gè)烴組分相平衡方程可包括：氣組分相平衡方程和油組分相平衡方程。在現(xiàn)有黑油模擬中，不存在關(guān)于相平衡的描述，而在現(xiàn)有的組分模擬中需要相平衡描述，所以，本實(shí)施例中通過(guò)獨(dú)創(chuàng)設(shè)計(jì)的相平衡方程能夠有效實(shí)現(xiàn)黑油模擬和組分模擬的結(jié)合。所述氣組分相平衡方程可為：Ffg(Xgg,Xgo)=XggPg-KgXgoPg=Xgg-KgXgo=0,---(4)]]>所述油組分相平衡方程可為：Ffo(Xog,Xoo)=XogPg-KoXooPg=Xog-KoXoo=0.---(5)]]>一個(gè)實(shí)施例中，相內(nèi)部組分構(gòu)成的約束方程可包括：氣相內(nèi)部組分約束方程和油相內(nèi)部組分約束方程。本實(shí)施例中，通過(guò)各相內(nèi)部組分構(gòu)成的約束方程可以保證歸一化，從而保證黑油油藏模擬結(jié)果的正確性。所述氣相內(nèi)部組分約束方程可為：Fcg(Xog,Xgg)=Xog+Xgg-1=0,---(6)]]>所述油相內(nèi)部組分約束方程可為：Fco(Xoo,Xgo)=Xoo+Xgo-1=0.---(7)]]>在上述公式(1)～(7)中，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ù)；Pg為氣相壓力(表示油藏壓力)，Sg為氣相飽和度，So為油相飽和度，Xgg為氣組分在氣相中摩爾濃度，Xgo為氣組分在油相中摩爾濃度，Xog為油組分在氣相中摩爾濃度，Xoo為油組分在油相中的摩爾濃度；t為時(shí)間，為孔隙度，ρg為氣相密度，ρo為油相密度，K為滲透率，Krg為氣相相對(duì)滲透率，μg為氣相粘度，g為重力加速度，D為油藏深度(油藏地質(zhì)模型中離散網(wǎng)格中心所對(duì)應(yīng)的深度)，Kro為油相相對(duì)滲透率，μo為油相粘度，Pcgo為氣油兩相之間的毛管力，qg為氣相產(chǎn)量或注入量，qo為油相產(chǎn)量或注入量，ρw為水相密度，Krw為水相相對(duì)滲透率，μw為水相粘度，Pcow為油水兩相之間的毛管力，qw為水相產(chǎn)量或注入量。圖11是本發(fā)明一實(shí)施例中組分?jǐn)?shù)值模型生成單元的結(jié)構(gòu)示意圖。如圖11所示，所述組分?jǐn)?shù)值模型生成單元220，可包括：組分模型離散模塊221和組分?jǐn)?shù)值模型生成模塊222，二者相互連接。組分模型離散模塊221用于執(zhí)行：根據(jù)所述黑油模擬所需數(shù)據(jù)，基于所述黑油油藏的地質(zhì)模型離散網(wǎng)格，對(duì)偏微分方程形式的所述氣組分守恒方程、所述油組分守恒方程及所述水組分守恒方程進(jìn)行空間數(shù)值離散，得到相應(yīng)的非線性代數(shù)方程。組分?jǐn)?shù)值模型生成模塊222用于執(zhí)行：根據(jù)所述氣組分相平衡方程、所述油組分相平衡方程、所述氣相內(nèi)部組分約束方程、所述油相內(nèi)部組分約束方程及所有所述非線性代數(shù)方程，按牛頓迭代法計(jì)算得到包含雅可比矩陣單元的線性代數(shù)方程組，作為所述離散網(wǎng)格的組分?jǐn)?shù)值模型。在組分模型離散模塊221中，可以采用多種方法對(duì)偏微分方程形式的所述氣組分守恒方程、所述油組分守恒方程及所述水組分守恒方程進(jìn)行空間數(shù)值離散，例如有限差分法、有限元法等。較佳地，是采用有限體積法對(duì)偏微分方程形式的所述氣組分守恒方程、所述油組分守恒方程及所述水組分守恒方程進(jìn)行空間數(shù)值離散，以此可以提高數(shù)值離散的精度。本實(shí)施例中，該組分?jǐn)?shù)值模型可以是包含雅可比矩陣單元的線性代數(shù)方程組可以簡(jiǎn)化計(jì)算。在其他實(shí)施例中，可以利用其他方法對(duì)所述氣組分守恒方程、所述油組分守恒方程及所述水組分守恒方程進(jìn)行空間數(shù)值離散。在其他實(shí)施例中，可以利用其他方法得到上述線性代數(shù)方程組。一個(gè)實(shí)施例中，通過(guò)為線性代數(shù)方程組選用適當(dāng)?shù)慕M分滲流參數(shù)以及安排多個(gè)組分滲流參數(shù)在線性代數(shù)方程組的變量矩陣中的順序，可使雅可比矩陣單元得到顯著簡(jiǎn)化，從而減少計(jì)算量。一個(gè)實(shí)施例中，可采用有限體積方法對(duì)組分模型中3個(gè)組分守恒方程(偏微分方程)進(jìn)行空間離散。例如，對(duì)上述公式(1)～(3)進(jìn)行空間離散。所述氣組分守恒方程所對(duì)應(yīng)的非線性代數(shù)方程可為：所述油組分守恒方程所對(duì)應(yīng)的非線性代數(shù)方程可為：所述水組分守恒方程所對(duì)應(yīng)的非線性代數(shù)方程可為：在上述公式(8)～(10)中，F(xiàn)g表示描述氣組分守恒關(guān)系的函數(shù)，F(xiàn)o表示描述油組分守恒關(guān)系的函數(shù)，F(xiàn)w表示描述水組分守恒關(guān)系的函數(shù)，Pg為氣相壓力，Sg為氣相飽和度，So為油相飽和度，Xgg為氣組分在氣相中摩爾濃度，Xgo為氣組分在油相中摩爾濃度，Xog為油組分在氣相中摩爾濃度，Xoo為油組分在油相中的摩爾濃度；Ω為空間積分變量，為孔隙度，ρg為氣相密度，ρo為油相密度，K為滲透率，Krg為氣相相對(duì)滲透率，μg為氣相粘度，g為重力加速度，D為油藏深度，Kro為油相相對(duì)滲透率，μo為油相粘度，Pcgo為氣油兩相之間的毛管力，qg為氣相產(chǎn)量或注入量，qo為油相產(chǎn)量或注入量，ρw為水相密度，Krw為水相相對(duì)滲透率，μw為水相粘度，Pcow為油水兩相之間的毛管力，qw為水相產(chǎn)量或注入量，∑l′表示對(duì)離散網(wǎng)格所有表面l′的通量求和，Al′表示離散網(wǎng)格表面l′的面積，V表示離散網(wǎng)格的體積，Qg＝V×qg表示離散網(wǎng)格的氣相注入或產(chǎn)出率，Qo＝V×qo表示離散網(wǎng)格的油相注入或產(chǎn)出率，Qw＝V×qw表示離散網(wǎng)格的水相注入或產(chǎn)出率。Ω表示對(duì)一個(gè)離散網(wǎng)格進(jìn)行空間積分的變量。在組分?jǐn)?shù)值模型生成模塊222中，所述氣組分相平衡方程、所述油組分相平衡方程、所述氣相內(nèi)部組分約束方程、所述油相內(nèi)部組分約束方程也可為非線性代數(shù)方程。根據(jù)所述氣組分相平衡方程、所述油組分相平衡方程、所述氣相內(nèi)部組分約束方程、所述油相內(nèi)部組分約束方程及所有所述非線性代數(shù)方程，按牛頓迭代法計(jì)算得到包含雅可比矩陣單元的線性代數(shù)方程組，形成離散網(wǎng)格的組分?jǐn)?shù)值模型。在其他實(shí)施例中，可以采用其他方法得到線性代數(shù)方程組。一個(gè)實(shí)施例中，按牛頓迭代法計(jì)算得到線性代數(shù)方程組可為：(∂F∂X)lδX=-Fl,---(11)]]>其中，F(xiàn)表示函數(shù)，X表示變量，一般情況下，δX為變化量，表示雅可比矩陣(方程組F對(duì)基本求解變量X的導(dǎo)數(shù))，不同的l值用于區(qū)別不同的線性代數(shù)方程組，l可表示迭代次數(shù)，其他符號(hào)的意義可與其他實(shí)施例中同樣符號(hào)的意義相同。本實(shí)施例中，變量X中選取的參數(shù)及各參數(shù)的排列順序可以大大簡(jiǎn)化雅可比矩陣?yán)?，雅克比矩陣單元的一般形式可表示為：如上述雅可比矩陣所示，采用本?shí)施例所建方程組得到的雅可比矩陣單元可為7行7列矩陣，可共有49個(gè)元素，其中，前17個(gè)元素需要計(jì)算得到，剩余元素中有2個(gè)需要查詢組分相平衡參數(shù)表(即參數(shù)-Kg和-Ko)，其余30個(gè)元素取值固定為0或者1，矩陣的這一特征有利于提高模擬速度。圖12是本發(fā)明一實(shí)施例中組分相平衡參數(shù)表建立單元的結(jié)構(gòu)示意圖。如圖12所示，所述組分相平衡參數(shù)表建立單元230，可包括：組分摩爾比生成模塊231、組分平衡常數(shù)生成模塊232及組分相平衡參數(shù)表建立模塊233，上述各模塊順序連接。組分摩爾比生成模塊231用于執(zhí)行：根據(jù)所述流體PVT數(shù)據(jù)中的油相密度、油相中溶解氣油比、氣相密度及氣相中揮發(fā)油氣比，確定油相中溶解氣組分與油組分的摩爾比和氣相中揮發(fā)油組分與氣組分的摩爾比。組分平衡常數(shù)生成模塊232用于執(zhí)行：根據(jù)所述的油相中溶解氣組分與油組分的摩爾比和所述的氣相中揮發(fā)油組分與氣組分的摩爾比，計(jì)算得到氣組分平衡常數(shù)和油組分平衡常數(shù)。組分相平衡參數(shù)表建立模塊233用于執(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)的飽和壓力。本實(shí)施例中，發(fā)明人經(jīng)過(guò)創(chuàng)造性勞動(dòng)，想到通過(guò)述流體PVT數(shù)據(jù)的油相中溶解氣油比和氣相中揮發(fā)氣油比計(jì)算得到溶解油組分與溶解氣組分的摩爾比和揮發(fā)油組分與揮發(fā)氣組分的摩爾比，進(jìn)而計(jì)算得到氣組分平衡常數(shù)和油組分平衡常數(shù)，再計(jì)算得到氣組分摩爾分?jǐn)?shù)和油組分摩爾分?jǐn)?shù)，最終實(shí)現(xiàn)組分相平衡參數(shù)表的建立。利用該組分相平衡參數(shù)表中組分摩爾分?jǐn)?shù)、組分平衡常數(shù)及飽和壓力的對(duì)應(yīng)關(guān)系，能夠根據(jù)組分模擬中的組分摩爾分?jǐn)?shù)參數(shù)數(shù)據(jù)查找得到黑油模擬中的飽和壓力參數(shù)數(shù)據(jù)，從而實(shí)現(xiàn)組分模擬和黑油模擬的結(jié)合，即實(shí)現(xiàn)用組分模型實(shí)現(xiàn)黑油模擬。而且，現(xiàn)有黑油模擬僅能得到油藏壓力和油氣飽和度數(shù)據(jù)，本發(fā)明通過(guò)組分相平衡參數(shù)表還可以得到油氣組分摩爾分?jǐn)?shù)以及流體PVT中的多種數(shù)據(jù)，能夠滿足更多數(shù)據(jù)需求，具有更好的擴(kuò)展性。一個(gè)實(shí)施例中，在組分摩爾比生成模塊231中，根據(jù)油相中溶解氣油比Rgo(單位：SCF/STB)確定溶解氣組分與油組分的摩爾比rg，計(jì)算公式可如下所示：rg=(Rgo×ρgMrg)/(1×158.987×0.0353×ρoMro),---(13)]]>其中，Mrg為氣組分分子量，例如可為23；Mro為油組分分子量，例如可為190；分子量Mrg、Mro的取值不影響模擬結(jié)果。ρg為氣相質(zhì)量密度(SCF/STB)，ρo為油相質(zhì)量密度(SCF/STB)，其它參數(shù)為單位轉(zhuǎn)換的常數(shù)。在摩爾比生成模塊231中，根據(jù)氣相中揮發(fā)油氣比Rog(單位：SCF/STB)換算出揮發(fā)油組分與氣組分的摩爾比ro，計(jì)算公式可如下所示：ro=(Rog×158.987×0.0353ρoMro)/(1×ρgMrg).---(14)]]>一個(gè)實(shí)施例中，在平衡常數(shù)生成模塊232中，根據(jù)溶解氣組分與油組分的摩爾比rg和揮發(fā)油組分與氣組分的摩爾比ro，計(jì)算出氣組分的平衡常數(shù)Kg和油組分的平衡常數(shù)Ko，計(jì)算公式可如下所示：Kg=1+1/rg1+ro,---(15)]]>Ko=1+rg1+1/ro.---(16)]]>一個(gè)實(shí)施例中，在組分相平衡參數(shù)表建立模塊233中，根據(jù)氣組分的平衡常數(shù)Kg和油組分的平衡常數(shù)Ko，計(jì)算氣組分總的摩爾分?jǐn)?shù)Zg和油組分總的摩爾分?jǐn)?shù)Zo，計(jì)算公式可如下所示：Zg=1-KoKg-Ko,---(17)]]>Zo=Kg-1Kg-Ko.---(18)]]>圖13是本發(fā)明另一實(shí)施例中黑油油藏模擬裝置的結(jié)構(gòu)示意圖。如圖13所示，圖10所示的黑油油藏模擬裝置(包括單元210～270)，還可包括：初始估計(jì)值生成單元290，與流體相態(tài)確定單元240連接。其他實(shí)施例中，還可包括單元280。初始估計(jì)值生成單元290用于執(zhí)行：根據(jù)所述黑油模擬所需數(shù)據(jù)中的所述黑油油藏的地質(zhì)模型確定所述離散網(wǎng)格的各流體相飽和度初始值，利用靜水柱壓力公式計(jì)算所述離散網(wǎng)格的油藏壓力初始值，通過(guò)查詢所述組分相平衡參數(shù)表得到所述離散網(wǎng)格的各組分總摩爾分?jǐn)?shù)初始值，各流體相飽和度初始值、油藏壓力初始值及各組分總摩爾分?jǐn)?shù)初始值一同作為初始的所述多個(gè)組分滲流參數(shù)的估計(jì)值。在初始估計(jì)值生成單元290中，上述地質(zhì)模型可以是以網(wǎng)格單元方式描述油藏幾何形狀、油氣水分布、孔隙度分布情況的數(shù)據(jù)體。地質(zhì)模型中的網(wǎng)格單元可以和上述離散網(wǎng)格相對(duì)應(yīng)。離散網(wǎng)格的各流體相飽和度可以包括氣相飽和度、油飽和度等。油藏飽和度信息可直接利用地質(zhì)模型中飽和度分布信息確定。一個(gè)實(shí)施例中，可依據(jù)流體PVT數(shù)據(jù)中的多種不同數(shù)據(jù)，查詢組分相平衡參數(shù)表得到離散網(wǎng)格的各組分總摩爾分?jǐn)?shù)初始值。上述確定的油藏壓力、飽和度和組分總摩爾分?jǐn)?shù)的初始估計(jì)值，在黑油模擬過(guò)程中可只作為迭代計(jì)算的初始值，后期隨著迭代進(jìn)行，油藏壓力、飽和度和組分總摩爾分?jǐn)?shù)的估計(jì)值可不斷調(diào)整。一些實(shí)施例中，油藏飽和度信息可利用前次迭代計(jì)算中的信息得到。圖14是本發(fā)明一實(shí)施例中流體相態(tài)確定單元的結(jié)構(gòu)示意圖。如圖14所示，所述流體相態(tài)確定單元240，可包括：網(wǎng)格摩爾分?jǐn)?shù)生成模塊241、網(wǎng)格飽和壓力生成模塊242及網(wǎng)格流體相態(tài)確定模塊243，上述各模塊順序連接。網(wǎng)格摩爾分?jǐn)?shù)生成模塊241用于執(zhí)行：根據(jù)所述離散網(wǎng)格的所述組分滲流參數(shù)的當(dāng)前估計(jì)值計(jì)算得到所述離散網(wǎng)格的組分總摩爾分?jǐn)?shù)。網(wǎng)格飽和壓力生成模塊242用于執(zhí)行：根據(jù)所述組分總摩爾分?jǐn)?shù)查詢所述組分相平衡參數(shù)表，得到所述離散網(wǎng)格中組分構(gòu)成所對(duì)應(yīng)的飽和壓力。網(wǎng)格流體相態(tài)確定模塊243用于執(zhí)行：通過(guò)對(duì)比所述離散網(wǎng)格的飽和壓力和所述組分滲流參數(shù)中的油藏壓力的當(dāng)前估計(jì)值，確定所述離散網(wǎng)格的流體相態(tài)和相應(yīng)物性參數(shù)。在網(wǎng)格摩爾分?jǐn)?shù)生成模塊241中，在第一次迭代過(guò)程中，組分滲流參數(shù)的當(dāng)前估計(jì)值可以是上述初始的估計(jì)值。上述離散網(wǎng)格的組分的總摩爾分?jǐn)?shù)可以是氣組分的總摩爾分?jǐn)?shù)或油組分的總摩爾分?jǐn)?shù)。在網(wǎng)格飽和壓力生成模塊242中，可根據(jù)氣組分的總摩爾分?jǐn)?shù)或油組分的總摩爾分?jǐn)?shù)查詢所述組分相平衡參數(shù)表，得到所述離散網(wǎng)格的組分構(gòu)成所對(duì)應(yīng)的飽和壓力。一個(gè)實(shí)施例中，根據(jù)組分的總摩爾分?jǐn)?shù)，可通過(guò)線性插值法查詢組分相平衡參數(shù)表，得到所述離散網(wǎng)格的組分構(gòu)成所對(duì)應(yīng)的飽和壓力，以此可以更準(zhǔn)確地查詢得到離散網(wǎng)格的組分構(gòu)成情況下的飽和壓力。在其他實(shí)施例中，可以采用其他方法查詢組分相平衡參數(shù)表得到離散網(wǎng)格的飽和壓力。一個(gè)實(shí)施例中，所述離散網(wǎng)格的組分可包括油組分和氣組分，所述組分滲流參數(shù)的當(dāng)前估計(jì)值可包括：離散網(wǎng)格當(dāng)前的氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度及油組分在油相中的摩爾濃度?？筛鶕?jù)當(dāng)前網(wǎng)格氣相飽和度、油相飽和度、氣組分在氣相中摩爾濃度、氣組分在油相中摩爾濃度、油組分在氣相中摩爾濃度、油組分在油相中的摩爾濃度，計(jì)算出油組分總的摩爾分?jǐn)?shù)和氣組分總的摩爾分?jǐn)?shù)。一個(gè)實(shí)施例中，計(jì)算油組分總的摩爾分?jǐn)?shù)和氣組分總的摩爾分?jǐn)?shù)公式可如下所示：油組分的總摩爾分?jǐn)?shù)可為：Zo=(SgρgXog+SoρoXoo)(SgρgXog+SoρoXoo+SgρgXgg+SoρoXgo),---(20)]]>氣組分的總摩爾分?jǐn)?shù)可為：Zg=(SgρgXgg+SoρoXgo)(SgρgXog+SoρoXoo+SgρgXgg+SoρoXgo)=1-Zo,---(21)]]>其中，Zo為油組分的總摩爾分?jǐn)?shù)，Zg為氣組分的總摩爾分?jǐn)?shù)，Sg為氣相飽和度，So為油相飽和度，ρg為氣相密度，ρo為油相密度，Xog為油組分在氣相中摩爾濃度，Xoo為油組分在油相中的摩爾濃度，Xgg為氣組分在氣相中摩爾濃度，Xgo為氣組分在油相中摩爾濃度。本實(shí)施例中，獨(dú)創(chuàng)地利用組分滲流參數(shù)的估計(jì)值計(jì)算出離散網(wǎng)格內(nèi)各組分的總摩爾分?jǐn)?shù)，有利于準(zhǔn)確地從組分相平衡參數(shù)表中查找到飽和壓力。圖15是本發(fā)明一實(shí)施例中網(wǎng)格流體相態(tài)確定模塊的結(jié)構(gòu)示意圖。如圖15所示，所述網(wǎng)格流體相態(tài)確定模塊243可包括：油相相態(tài)判斷模塊2431、液相比例生成模塊2432及油氣兩相、氣相判斷模塊2433及物性參數(shù)確定模塊2434，上述各模塊順序連接。油相相態(tài)判斷模塊2431用于執(zhí)行：根據(jù)所述離散網(wǎng)格的油藏壓力的當(dāng)前估計(jì)值大于所述飽和壓力判定所述離散網(wǎng)格的流體相態(tài)為油相。液相比例生成模塊2432用于執(zhí)行：當(dāng)所述離散網(wǎng)格的油藏壓力的當(dāng)前估計(jì)值小于所述飽和壓力時(shí)，根據(jù)所述離散網(wǎng)格的油藏壓力的當(dāng)前估計(jì)值查詢所述組分相平衡參數(shù)表，得到所述離散網(wǎng)格的氣組分平衡常數(shù)和油組分平衡常數(shù)，并根據(jù)所述的氣組分平衡常數(shù)和油組分平衡常數(shù)計(jì)算所述離散網(wǎng)格的液相比例。油氣兩相及氣相判斷模塊2433用于執(zhí)行：根據(jù)所述液相比例大于零判定所述離散網(wǎng)格的流體相態(tài)為油氣兩相共存，根據(jù)所述液相比例等于或小于零判定所述離散網(wǎng)格的流體相態(tài)為氣相。物性參數(shù)確定模塊2434用于執(zhí)行：根據(jù)所述離散網(wǎng)格的流體相態(tài)的判定結(jié)果和油藏壓力的當(dāng)前估計(jì)值，查詢所述流體PVT數(shù)據(jù)得到相應(yīng)物性參數(shù)。在上油相相態(tài)判斷模塊2431中，當(dāng)離散網(wǎng)格的油藏壓力Pg的當(dāng)前估計(jì)值大于飽和壓力Pb，即Pg＞Pb時(shí)，網(wǎng)格內(nèi)部不存在氣體，所有烴類(lèi)組分都存在于油相中，可以認(rèn)為離散網(wǎng)格的流體相態(tài)為油相。若判斷已離散網(wǎng)格的油藏壓力Pg的當(dāng)前估計(jì)值大于飽和壓力Pb，則不需作進(jìn)一步判斷。在液相比例生成模塊2432中，當(dāng)離散網(wǎng)格的油藏壓力Pg的當(dāng)前估計(jì)值小于飽和壓力Pb，即Pg＜Pb時(shí)，需要進(jìn)一步確定離散網(wǎng)格的流體相態(tài)是油氣兩相還是氣相，例如可以利用離散網(wǎng)格的液相比例作進(jìn)一步判斷。一個(gè)實(shí)施例中，上述液相比例生成模塊2432中，該液相比例可為：V′=Zg1-Ko+Zo1-Kg,---(22)]]>其中，V'表示液相比例，Zg為氣組分的總摩爾分?jǐn)?shù)，Zo為油組分的總摩爾分?jǐn)?shù)，Ko油組分平衡常數(shù)，Kg為氣組分平衡常數(shù)。一個(gè)實(shí)施例中，如果液相比例V'＞0，則油氣兩相共存，油氣組分按照平衡參數(shù)取值在液氣兩相中分配；如果液相比例V'＝0，則只有氣相存在，所有組分都只存在氣相中。一個(gè)實(shí)施例中，在組分?jǐn)?shù)值模型調(diào)整單元250中，根據(jù)所述流體相態(tài)和相應(yīng)物性參數(shù)調(diào)整所述組分?jǐn)?shù)值模型。根據(jù)所述流體相態(tài)調(diào)整組分?jǐn)?shù)值模型包括三種情況，在調(diào)整后的所述組分?jǐn)?shù)值模型中：當(dāng)所述離散網(wǎng)格的流體相態(tài)為油氣兩相共存時(shí)，變量可保持上述變量的形式，雅可比矩陣單元可采用上述雅可比矩陣的一般形式(式(12))：當(dāng)所述離散網(wǎng)格的流體相態(tài)僅為油相，即網(wǎng)格內(nèi)僅存在油相時(shí)，上述變量X可調(diào)整為變量所述雅可比矩陣單元可調(diào)整為：當(dāng)所述離散網(wǎng)格的流體相態(tài)僅為氣相，即網(wǎng)格內(nèi)僅存在氣相時(shí)，上述變量X可調(diào)整為變量所述雅可比矩陣單元可調(diào)整為：一個(gè)實(shí)施例中，在確定雅可比系數(shù)矩陣和右端項(xiàng)-Fl的情況下，利用線性求解算法計(jì)算如下方程組，獲得δX的值。(∂F∂X)lδX=-Fl.---(25)]]>本實(shí)施例中，由上述各變化量δX的求解結(jié)果可以看出，調(diào)整變量X后，組分?jǐn)?shù)值模型得到簡(jiǎn)化，計(jì)算過(guò)程也得到簡(jiǎn)化，由此可提高黑油模擬的效率。一個(gè)實(shí)施例中，所述油藏參數(shù)確定單元270，可包括：收斂判斷模塊271。收斂判斷模塊271用于執(zhí)行：根據(jù)變化值δXl和相應(yīng)估計(jì)值Xl是否滿足判定所述的油藏壓力變化值及各流體相飽和度變化值收斂，其中，變化值δXl表示第l次迭代后油藏壓力變化值或各流體相飽和度變化值，相應(yīng)估計(jì)值Xl表示前次迭代后油藏壓力的估計(jì)值或各流體相飽和度的估計(jì)值，迭代次數(shù)l≥1，X1為初始的估計(jì)值，a為設(shè)定值，例如0＜a＜0.1。較佳實(shí)施例中，設(shè)定值a＝0.01，當(dāng)變化值與估計(jì)值之比小于一設(shè)定值，即判斷收斂，則可進(jìn)行模擬結(jié)果輸出處理；否則認(rèn)為不收斂，更新估計(jì)值，更新后的估計(jì)值Xl+1＝Xl+δXl，重新計(jì)算直到變化值收斂，將最新的變化值和前次計(jì)算的估計(jì)值相加作為最終估計(jì)值。此處，估計(jì)值可以是油藏壓力估計(jì)值、油相飽和度估計(jì)值等。本發(fā)明實(shí)施例的黑油模擬裝置，采集黑油模擬所需的數(shù)據(jù)；根據(jù)采集的流體PVT數(shù)據(jù)，建立組分相平衡參數(shù)表；建立一個(gè)描述黑油油藏流體滲流規(guī)律、包含油藏壓力、油氣飽和度和油氣組分摩爾分?jǐn)?shù)等信息的組分模型；根據(jù)采集的數(shù)據(jù)對(duì)組分模型進(jìn)行數(shù)值離散，確定油藏壓力、飽和度和組分摩爾分?jǐn)?shù)的估計(jì)值，并形成初始的組分?jǐn)?shù)值模型；根據(jù)相平衡參數(shù)和油藏估計(jì)值進(jìn)行閃蒸計(jì)算，確定流體相態(tài)和構(gòu)成；根據(jù)閃蒸計(jì)算得到的流體相態(tài)，調(diào)整組分?jǐn)?shù)值模型；根據(jù)組分?jǐn)?shù)值模型，線性求解得到油藏壓力、油氣飽和度信息；判斷油藏壓力和飽和度信息是否滿足收斂要求，收斂則輸出模擬結(jié)果，不收斂，調(diào)整估計(jì)值，重新計(jì)算。本發(fā)明實(shí)施例的黑油模擬裝置，利用組分模型和模擬器架構(gòu)實(shí)現(xiàn)了黑油模擬功能，提高了模擬器集成性和后期擴(kuò)展性；直接利用現(xiàn)有黑油模擬的輸入數(shù)據(jù)，不需要額外給出組分相平衡信息，方便用戶使用；利用組分相平衡參數(shù)表的閃蒸計(jì)算，保證了模擬效率。本實(shí)施例的黑油油藏模擬裝置，獨(dú)創(chuàng)地建立描述黑油油藏流體滲流規(guī)律的組分模型，進(jìn)一步得到離散網(wǎng)格的組分?jǐn)?shù)值模型，奠定了用組分模型進(jìn)行黑油模擬的基礎(chǔ)，再利用流體PVT數(shù)據(jù)建立組分相平衡參數(shù)表，為用組分模型進(jìn)行黑油模擬搭建了橋梁，實(shí)現(xiàn)了組分模擬和現(xiàn)有基于黑油模型黑油模擬的有效結(jié)合。在本說(shuō)明書(shū)的描述中，參考術(shù)語(yǔ)“一個(gè)實(shí)施例”、“一個(gè)具體實(shí)施例”、“一些實(shí)施例”、“例如”、“示例”、“具體示例”、或“一些示例”等的描述意指結(jié)合該實(shí)施例或示例描述的具體特征、結(jié)構(gòu)、材料或者特點(diǎn)包含于本發(fā)明的至少一個(gè)實(shí)施例或示例中。在本說(shuō)明書(shū)中，對(duì)上述術(shù)語(yǔ)的示意性表述不一定指的是相同的實(shí)施例或示例。而且，描述的具體特征、結(jié)構(gòu)、材料或者特點(diǎn)可以在任何的一個(gè)或多個(gè)實(shí)施例或示例中以合適的方式結(jié)合。各實(shí)施例中涉及的步驟順序用于示意性說(shuō)明本發(fā)明的實(shí)施，其中的步驟順序不作限定，可根據(jù)需要作適當(dāng)調(diào)整。本領(lǐng)域內(nèi)的技術(shù)人員應(yīng)明白，本發(fā)明的實(shí)施例可提供為方法、系統(tǒng)、或計(jì)算機(jī)程序產(chǎn)品。因此，本發(fā)明可采用完全硬件實(shí)施例、完全軟件實(shí)施例、或結(jié)合軟件和硬件方面的實(shí)施例的形式。而且，本發(fā)明可采用在一個(gè)或多個(gè)其中包含有計(jì)算機(jī)可用程序代碼的計(jì)算機(jī)可用存儲(chǔ)介質(zhì)(包括但不限于磁盤(pán)存儲(chǔ)器、CD-ROM、光學(xué)存儲(chǔ)器等)上實(shí)施的計(jì)算機(jī)程序產(chǎn)品的形式。本發(fā)明是參照根據(jù)本發(fā)明實(shí)施例的方法、設(shè)備(系統(tǒng))、和計(jì)算機(jī)程序產(chǎn)品的流程圖和/或方框圖來(lái)描述的。應(yīng)理解可由計(jì)算機(jī)程序指令實(shí)現(xiàn)流程圖和/或方框圖中的每一流程和/或方框、以及流程圖和/或方框圖中的流程和/或方框的結(jié)合?？商峁┻@些計(jì)算機(jī)程序指令到通用計(jì)算機(jī)、專(zhuān)用計(jì)算機(jī)、嵌入式處理機(jī)或其他可編程數(shù)據(jù)處理設(shè)備的處理器以產(chǎn)生一個(gè)機(jī)器，使得通過(guò)計(jì)算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備的處理器執(zhí)行的指令產(chǎn)生用于實(shí)現(xiàn)在流程圖一個(gè)流程或多個(gè)流程和/或方框圖一個(gè)方框或多個(gè)方框中指定的功能的裝置。這些計(jì)算機(jī)程序指令也可存儲(chǔ)在能引導(dǎo)計(jì)算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備以特定方式工作的計(jì)算機(jī)可讀存儲(chǔ)器中，使得存儲(chǔ)在該計(jì)算機(jī)可讀存儲(chǔ)器中的指令產(chǎn)生包括指令裝置的制造品，該指令裝置實(shí)現(xiàn)在流程圖一個(gè)流程或多個(gè)流程和/或方框圖一個(gè)方框或多個(gè)方框中指定的功能。這些計(jì)算機(jī)程序指令也可裝載到計(jì)算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備上，使得在計(jì)算機(jī)或其他可編程設(shè)備上執(zhí)行一系列操作步驟以產(chǎn)生計(jì)算機(jī)實(shí)現(xiàn)的處理，從而在計(jì)算機(jī)或其他可編程設(shè)備上執(zhí)行的指令提供用于實(shí)現(xiàn)在流程圖一個(gè)流程或多個(gè)流程和/或方框圖一個(gè)方框或多個(gè)方框中指定的功能的步驟。以上所述的具體實(shí)施例，對(duì)本發(fā)明的目的、技術(shù)方案和有益效果進(jìn)行了進(jìn)一步詳細(xì)說(shuō)明，所應(yīng)理解的是，以上所述僅為本發(fā)明的具體實(shí)施例而已，并不用于限定本發(fā)明的保護(hù)范圍，凡在本發(fā)明的精神和原則之內(nèi)，所做的任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本發(fā)明的保護(hù)范圍之內(nèi)。當(dāng)前第1頁(yè)1 2 3