低滲儲(chǔ)層油水兩相非達(dá)西滲流數(shù)值模擬研究.pdf

1、目前一般商用模擬器的基本滲流方程是建立在達(dá)西滲流方程基礎(chǔ)之上，因而其用在具有非達(dá)西滲流特征的低滲油藏上將會(huì)對(duì)計(jì)算預(yù)測(cè)結(jié)果產(chǎn)生較大的誤差。低速非達(dá)西滲流數(shù)值方法的研究工作尚少，更難見(jiàn)于成熟的商業(yè)模擬器，因此本文在總結(jié)分析低滲透儲(chǔ)層非達(dá)西滲流機(jī)理及幾種為國(guó)內(nèi)外學(xué)者所認(rèn)可的非達(dá)西滲流數(shù)學(xué)擬合模型在使用上的優(yōu)缺點(diǎn)的基礎(chǔ)上，推導(dǎo)出了具有變啟動(dòng)壓力梯度的油水兩相三維非達(dá)西滲流模型，用控制體積有限差分方法對(duì)該滲流模型進(jìn)行了的離散化處理，并研究探討了非

2、線性項(xiàng)的全隱式線性化方法和井底流壓隱式處理方法；另外，本文應(yīng)用Newton-Raphson迭代算法原理推導(dǎo)了油水兩相的迭代方程，并給出迭代方程聯(lián)立解法矩陣形式，而且還詳盡的描述了組合ILU(0)預(yù)處理技術(shù)和廣義最小殘量算法(GMRES)線性方程組求解方法，并用VC 6.0編制了包含GMRES、ILU(0)及稀疏矩陣壓縮存儲(chǔ)方法組的線性方程靜態(tài)類庫(kù)GMR；此外，本文還運(yùn)用面向?qū)ο蟮某绦蛟O(shè)計(jì)方法設(shè)計(jì)非達(dá)西滲流數(shù)值計(jì)算程序SimCal中網(wǎng)格、

3、流體、井、計(jì)算控制、求解算法、數(shù)據(jù)庫(kù)訪問(wèn)控制等類模塊。最后本文對(duì)比了Eclipse 100與本文SimCal的計(jì)算結(jié)果，并應(yīng)用SimCal研究了達(dá)西流、非達(dá)西流以及三種僅包含油相恒定啟動(dòng)壓力梯度滲流條件下的油井生產(chǎn)情況，研究結(jié)果表明：(1)在達(dá)西流滲流條件，SimCal計(jì)算結(jié)果是可信賴正確的，與Eclipse的平均偏差約為7.66％；SimCal非達(dá)西滲流計(jì)算結(jié)果是合理的；(2)在僅考慮油相啟動(dòng)壓力梯度情況下，啟動(dòng)壓力梯度越大，生產(chǎn)井的