油藏數(shù)值模擬基礎(chǔ)培訓(xùn)第六講課件

油藏數(shù)值模擬基礎(chǔ)培訓(xùn)（第六講）油藏數(shù)值模擬基礎(chǔ)培訓(xùn)求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCYREGION/SOLUTION要點分析REGION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCY油藏不同部位可能有不同的流體屬性，比如不同斷塊的油密度，粘度不同，或你的油藏巖性在縱向或垂向有變化，那你就需要在你的模型中設(shè)置流體或巖性分區(qū)。分區(qū)設(shè)置REGION：儲量分區(qū)FIPNUM：如果你想輸出模型不同部位的儲量，你需要設(shè)置儲量分區(qū)

流體分區(qū)PVTNUM：如果你的模型不同部位流體PVT屬性不同，你需要設(shè)置流體分區(qū)。

巖性分區(qū)SATNUM：如果你的模型不同部位巖性不同，需要用不同的相滲曲線和毛管壓力曲線，你需要設(shè)置巖性分區(qū)。

平衡區(qū)分區(qū)EQLNUM：如果你的模型有不同的油水或油氣界面，你需要設(shè)平衡區(qū)分區(qū)。REGION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCY分區(qū)設(shè)置REGION：REGION部分關(guān)鍵字REGION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCY儲量分區(qū)：REGION部分FIPNUM123…/流體分區(qū)：PROPS部分PVTOtable1/table2/table3/PVDGtable1/table2/table3/PVTWtable1/table2/table3/PROPS部分PVTWtable1/table2/table3/DENSITYtable1/table2/table3/ROCKtable1/table2/table3/REGION部分PVTNUM123…/REGION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCY巖性分區(qū)：PROPS部分SWOFtable1/table2/table3/SGOFtable1/table2/table3/REGION部分SATNUM123…/平衡區(qū)分區(qū)：SOLUTION部分EQUILline1/line2/line3/RSVDtable1/table2/table3/REGION部分EQLNUM123…/SOLUTION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCY（1）計算過渡帶高度。由油水界面和油氣界面深度以及相滲曲線提供的最大毛管壓力計算。

（2）計算每一個網(wǎng)格初始的油相，水相，氣相壓力分布。首先將在流體屬性部分提供的油，氣，水地面密度折算為地下密度?；趨⒖键c的深度和對應(yīng)壓力以及油水界面，油氣界面深度，過渡帶高度，結(jié)合油，氣，水地下密度計算其他深度處的油，氣，水相壓力。

（3)由每個網(wǎng)格的油，氣，水壓力計算油水和油氣毛管壓力

（4)計算飽和度分布。這部分計算主要用你提供的相滲曲線端點值。將油水界面以下的含水飽和度設(shè)為你在油水相滲曲線中提供的最大含水飽和度，通常為1。將油氣界面以上的含氣飽和度設(shè)為你提供的油氣相滲曲線的最大值。油氣界面以上的含水飽和度為束縛水飽和度。在油區(qū)的含油飽和度為1減束縛水飽和度。在過渡帶的含油和含水飽和度由你提供的毛管壓力曲線得到。

（5）計算初始溶解油氣比或泡點壓力的分布。模型初始化過程：SOLUTION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCY模型初始化過程：●平衡初始化根據(jù)已知的油藏壓力系數(shù)、油水及油氣界面、界面處的毛管壓力和平衡系統(tǒng)的相對滲透率曲線等，考慮重力及毛管壓力影響，建立油藏原始壓力及飽和度分布。油區(qū)：SWI氣頂：SWI油氣過渡帶：PCOG+SWI油水過渡帶：PCOW+SWI水區(qū)：100%SWSOLUTION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCYSOLUTION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCY各網(wǎng)格的初始壓力和相初始飽和度；各網(wǎng)格塊的初始溶解氣油比和（或）溶解油氣比；油藏流體屬性對深度的依賴關(guān)系，包括API、飽和溶解氣油比、泡點壓力、飽和溶解油氣比OGR和露點壓力隨深度的變化；解析水體的初始狀況。●平衡初始化SOLUTION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCY平衡法指定界面的深度、一個參考深度及此參考深度處的壓力；假設(shè)滿足流體靜力學(xué)平衡；平衡法適用于有水平界面的初始未鉆開油藏；這里可能需要輸入不同的毛管力；用黑油的狀態(tài)方程來計算各處各相的靜壓力；從PROPS部分的飽和度方程中得到各區(qū)域的相飽和度；過渡區(qū)中的Sw和Sg通過反查飽和度函數(shù)中的毛管力曲線得到?！衿胶獬跏蓟疭OLUTION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCY●平衡初始化平衡法計算初始壓力和飽和度是基于飽和度函數(shù)、流體界面、參考深度和參考壓力的。如果已知的是處于油區(qū)的一個參考深度處的壓力，則會循環(huán)求解油相的黑油狀態(tài)方程和油相的靜壓力方程。黑油狀態(tài)方程：油相靜壓力方程：通過以上方程得到各處的油相壓力。其中（r）和（s）分別代表油藏和地面條件，h代表深度。給定油水界面，用同樣的方法求解氣相和水相的狀態(tài)方程，最終可以得到油藏中各個位置處的各相靜壓力。SOLUTION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCY相初始飽和度在水區(qū)，Sw為最大水相飽和度SWU，是輸入的水相飽和度函數(shù)表中的最大含水飽和度值，此值常取值為1。Sg為最小值SGL，是輸入的氣相飽和度函數(shù)表中的最小含氣飽和度值，即原生氣飽和度Sgco，此值常取為0。而油相飽和度為So=1-SWU-SGL。在油區(qū)，含水和含氣飽和度都處于最小值，Sg=SGL=Sgco，Sw=SWL=Swco。則含油飽和度為So=1-SWL-SGL。通常，So=1-Swco。SOLUTION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCY過渡區(qū)內(nèi)的相初始飽和度在過渡區(qū)，相的飽和度是由毛管力控制的。毛管力的定義為：Pcow=Po-Pw和Pcog=Pg-Po。毛管力可以通過計算任意深度處相間的靜壓力差得到。然后，ECLIPSE會反查輸入的飽和度函數(shù)中的毛管壓力曲線來確定過渡區(qū)的含水和含氣飽和度。理論上，在Sw=0時，水相的毛管力曲線Pcow趨于無窮大。而實際上，毛管力曲線應(yīng)該在Swco處終止，以保證不會在飽和度函數(shù)的最低含水飽和度處出現(xiàn)不連續(xù)的情況。讓我們考慮一個穿越油水界面OWC的傾斜網(wǎng)格。通過計算網(wǎng)格的中心深度，及前面的方法所計算出來的初始飽和度Soi、Swi、Sgi，就可以得到此網(wǎng)格塊的飽和度。這種方法提供了一種完全穩(wěn)定的初始化方法，其相飽和度與靜壓力是一致的。這種方法為塊中心平衡。SOLUTION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCYSOLUTION=========================================EQUIL--DepthPiWOCPcow|wocGORPcog|GocRSVDRVVDACC1500150.015800.015000.01*1*10

SOLUTION=========================================EQUIL--DepthPiWOCPcow|wocGORPcog|GocRSVDRVVDACC1500150.015800.00.00.011*10

RSVD2000.096.113000.096.11/SOLUTION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCY重啟計算歷史擬合結(jié)束后需要進行產(chǎn)量預(yù)測，在進行產(chǎn)量預(yù)測計算時，不需要再從歷史擬合開始時進行計算，可以直接從歷史擬合結(jié)束的時間接著往下算。這種應(yīng)用上一次計算的輸出作為下一次計算的初始輸入計算就叫重啟計算。

要進行重啟計算，首先要定義重啟時間步的輸出?？梢杂肦PTRST來定義輸出每時間步，每月，每年或每隔幾月幾年重啟時間步文件。如果采用多文件格式輸出，則文件后綴為：.X0000,.X0001等，如果是單文件輸出，則輸出文件為.UNRST.

重啟文件記錄了每時間步模型壓力分布，飽和度分布，溶解油氣比分布，同時也記錄所有井的井位，射孔位置，產(chǎn)量控制。不過重啟文件沒有記錄垂直管流表（VFP表），所以在應(yīng)用垂直管流表時要記住重啟時需加上垂直管流表。SOLUTION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCY重啟計算

ECLIPSE有兩種重啟計算方法，快速重啟和完全重啟。

快速重啟不需要重新處理RUNSPEC,GRID,EDIT,PROPS和REGIONS部分，如果在歷史擬合計算時設(shè)了SAVE關(guān)鍵字，這些部分將保存在輸出的SAVE文件中，這樣在重啟計算時不用再計算傳導(dǎo)率。

快速重啟步驟:

在歷史擬合部分用SAVE和RPTRST要求輸出SAVE文件和重啟文件。

將歷史擬合文件拷貝為重啟文件。

刪掉所有SUMMARY以前部分。

用LOAD關(guān)鍵字裝載SAVE文件。

用RESTART設(shè)重啟時間。

在SCHEDULE部分用SKIPREST或刪掉所有歷史擬合時間步。

如果有VFP表，要保留VFP表。

增加新時間步進行預(yù)測計算。SOLUTION部分要點分析求實創(chuàng)新超越REALITY，INNOVATION,TRANSCENDENCY重啟計算

ECLIPSE有兩種重啟計算方法，快速重啟和完全重啟