本發(fā)明涉及從油氣藏中采油的方法，更具體地涉及在油氣藏內(nèi)的油藏推掃區(qū)中使用泡沫驅(qū)動增強采油技術(shù)的采油系數(shù)估算方法。

背景技術(shù)：

常規(guī)采油方法不能開采儲油巖中封閉的油附近的任何位置。人們使用很多類型的增強采油方法來提高開采量，包括氣體、液體或蒸汽驅(qū)動油藏以抽取殘余油。在1950年代，提出了泡沫驅(qū)油作為排送系統(tǒng)以提高采油系數(shù)。自從1960年代開始，引入了大量實驗和現(xiàn)場測試，并驗證了泡沫驅(qū)油在提高采油系數(shù)方面的有效性，尤其是對于非均質(zhì)油藏的有效性。

在油藏驅(qū)油操作中，通常使用表面活性劑來幫助在例如含水驅(qū)動溶液和疏水的油之間形成乳液，這會增強油從巖石中的分離。其他氣體或蒸汽驅(qū)動油方法通常也使用泡沫形成表面活性劑，泡沫驅(qū)油是一種提高采油量的方法，在該方法中將泡沫注入油藏中以提高驅(qū)油流體的推掃效率。泡沫能夠在油藏空洞中或注射前的表面處產(chǎn)生。泡沫驅(qū)油減輕了推掃的不均勻性，例如由具有比周圍地層更高的滲透性的層造成的那些或由重力超覆造成的那些不均勻性。

泡沫是在加入至少一種發(fā)泡劑的情況下將氣體(例如空氣、氮氣、co2等)混合到液相中所形成的分散體系。在石油工業(yè)中，用于產(chǎn)生泡沫的氣體選自空氣、氮氣、co2或天然氣，用于產(chǎn)生泡沫的液體主要是水，例如淡水、地層水或鹽水。

泡沫驅(qū)油還作為一種有效開采具有極低滲透率的致密油藏的方法被提出。從這些性質(zhì)的油氣藏中采油涉及的一個問題是對采油系數(shù)的估算。增強建?；蛴嬎銖倪@種驅(qū)動油藏中的采油系數(shù)的方法被認為將肯定會被接受。

當前用于估算從油氣藏中采油的采油系數(shù)的建模方法有兩個關(guān)鍵局限。這些局限包括計算作出的實際速度以及模型的準確性。從準確性的角度來看，準確估算經(jīng)受泡沫驅(qū)動的油氣藏中的采油系數(shù)是困難的，因為需要準確估算通過該油藏的推掃效率。隨著該油藏 的滲透性變低，或者該油藏的地質(zhì)變得更復(fù)雜，更難以找到適合的單一方程來合理準確地估算這些情況中的采油系數(shù)。這樣一來，從更準確建模和使用泡沫驅(qū)動增強采油技術(shù)的角度來看，克服了這些局限中的一些或者能得到比現(xiàn)有技術(shù)的方法更準確的結(jié)果的采油系數(shù)建模方法被認為是更期望的。

當前現(xiàn)有技術(shù)的方法中的其他局限之一是計算推掃效率和估算采油系數(shù)所需的其他變量所需的計算能力。根據(jù)來自油氣藏的地層的大量數(shù)據(jù)點等開發(fā)用于實現(xiàn)采油系數(shù)估算的單一大方程是在規(guī)劃或?qū)嵤┡菽?qū)油工程中有效使用這種建模技術(shù)的另一局限。如果可以找到需要比現(xiàn)有技術(shù)的方法更少的計算能力或更快地獲得結(jié)果的估算泡沫驅(qū)動油氣藏中的采油系數(shù)的改進方法，相信這將肯定會被接受。

技術(shù)實現(xiàn)要素：

本發(fā)明，一種泡沫驅(qū)動油藏中采油系數(shù)的估算方法，表現(xiàn)為對迄今為止的現(xiàn)有方法的改進，因為其提供了更精細的計算方法，導(dǎo)致能比目前使用的方法更精確地估算油氣藏的油藏推掃區(qū)內(nèi)的采油系數(shù)。

本發(fā)明，一種泡沫驅(qū)動油藏中采油系數(shù)的估算方法，還表現(xiàn)為對迄今為止現(xiàn)有的估算油氣藏內(nèi)采油系數(shù)的方法的改進，因為通過基于對該油藏和驅(qū)油參數(shù)的理解將該計算和估算細分為多個推掃子區(qū)使得能夠在驅(qū)油操作中對于每次驅(qū)油循環(huán)預(yù)測更快速地計算該模型的結(jié)果。

本發(fā)明，一種泡沫驅(qū)動油藏中采油系數(shù)的估算方法，達到了其目的，這包括在第一步驟中，將油藏推掃區(qū)定義為油氣藏內(nèi)在至少一個注射井和至少一個生產(chǎn)井之間的區(qū)域，其中在至少一個驅(qū)油循環(huán)中注射到該至少一個注射井中的氣體和液體表面活性劑將會把油通過該油藏推掃區(qū)推掃至該至少一個生產(chǎn)井以便從中采出。如此處所述，油藏推掃區(qū)的定義是油氣地質(zhì)和增強采油驅(qū)油工程規(guī)劃和實施領(lǐng)域的技術(shù)人員能夠理解的。從此處所述的建模方法的角度來看以及在所述驅(qū)油工程的更大規(guī)模規(guī)劃中，油藏推掃區(qū)的定義以及至少一個生產(chǎn)井和至少一個注射井的位置都是重要的初始考慮要素。

本發(fā)明的采油系數(shù)計算方法能夠用于計算驅(qū)油操作內(nèi)各驅(qū)油循環(huán)時油氣藏內(nèi)的采油系數(shù)。各驅(qū)油循環(huán)通常由以下構(gòu)成：將一股氣體和該股表面活性劑流體注入油氣藏的地質(zhì) 或地層中。盡管其計算和估算方法能夠用于估算單一股或驅(qū)油循環(huán)，但特別預(yù)期該方法在計算驅(qū)油操作內(nèi)多個驅(qū)油循環(huán)的油藏采油系數(shù)時將會是最適用的。

在定義了油藏推掃區(qū)之后，針對驅(qū)油操作中的一些或所有驅(qū)油循環(huán)而完成的該方法的下一個步驟包括根據(jù)注入該油藏推掃區(qū)中的氣體和液體表面活性劑的當前狀態(tài)最初測定該油藏推掃區(qū)的油藏和驅(qū)油參數(shù)。該油藏和驅(qū)油參數(shù)能夠尤其根據(jù)操作者、操作或油氣藏的地質(zhì)以及其他因素而變化。如此處所述，在該方法的某些實施方案中以及在本發(fā)明的計算中將會用到的各種油藏和驅(qū)油參數(shù)對于本領(lǐng)域技術(shù)人員而言將是顯而易見的，所有這些都包含在本發(fā)明的范圍內(nèi)。

在本發(fā)明的泡沫驅(qū)動油藏中采油系數(shù)的估算方法中的下一個步驟包括根據(jù)該油藏和驅(qū)油參數(shù)以及該操作或建模驅(qū)油操作的當前狀態(tài)，將該油藏推掃區(qū)細分為三個推掃子區(qū)：泡沫推掃子區(qū)、水推掃子區(qū)和氣體推掃子區(qū)。對于各驅(qū)油循環(huán)，重新定義所述推掃子區(qū)，因為在驅(qū)油操作過程中，各推掃子區(qū)的形狀和大小能夠根據(jù)不同的油藏和驅(qū)油參數(shù)而動態(tài)變化。

在定義了三個推掃子區(qū)之后，提供一個計算步驟，在該步驟中將計算以下各項：

·各推掃子區(qū)內(nèi)的縱向推掃效率和體積推掃效率；

·該油藏推掃區(qū)的總推掃面積，其是在當前驅(qū)油循環(huán)時被泡沫推掃的油藏推掃區(qū)的體積面積；

·泡沫推掃面積，其是在當前驅(qū)油循環(huán)時被泡沫推掃的泡沫推掃子區(qū)的體積面積；

·水推掃面積，其是在當前驅(qū)油循環(huán)時被水推掃的水推掃子區(qū)的體積面積；

·氣體推掃面積，其是在當前驅(qū)油循環(huán)時被氣體推掃的氣體推掃子區(qū)的體積面積；

·泡沫穿透面積，其是在該驅(qū)油循環(huán)內(nèi)被現(xiàn)有的泡沫穿透的泡沫推掃子區(qū)的體積面積，其包括泡沫推掃面積及其周圍由于泡沫與油之間流動性差異造成的泡沫旁路面積；

·水穿透面積，其是在該驅(qū)油循環(huán)內(nèi)被現(xiàn)有的水穿透的水推掃子區(qū)的體積面積，其包括水推掃面積及其周圍由于水與油之間流動性差異造成的水旁路面積；

·氣體穿透面積，其是在該驅(qū)油循環(huán)內(nèi)被現(xiàn)有的氣體穿透的氣體推掃子區(qū)的體積面積，其包括氣體推掃面積及其周圍由于氣體與油之間流動性差異造成的氣體旁路面積；和

·總穿透面積，其是泡沫穿透面積、水穿透面積和氣體穿透面積的總和。

在計算步驟中提供的計算可能也依賴于該油藏和驅(qū)油參數(shù)，且本領(lǐng)域技術(shù)人員能夠理解的是，在不脫離本發(fā)明的預(yù)期范圍的情況下，針對該計算步驟描述的各變量的具體計算能夠包括使用該油藏和驅(qū)油參數(shù)或者能夠以很多方式優(yōu)化以產(chǎn)生可能的最精確的結(jié)果。

在該計算步驟完成之后，此處所述方法中的下一步驟是采油系數(shù)估算步驟，在該步驟中使用計算步驟中確定的變量來估算各推掃子區(qū)中采油系數(shù)的貢獻。特別地，該采油系數(shù)估算步驟首先包括計算泡沫區(qū)采油系數(shù)，其是泡沫推掃子區(qū)貢獻的采油系數(shù)。還計算水區(qū)采油系數(shù)，其是水推掃子區(qū)貢獻的采油系數(shù)；以及氣體區(qū)采油系數(shù)，其是氣體推掃子區(qū)貢獻的采油系數(shù)。

如在計算步驟的情況下那樣，在本發(fā)明的方法的一些迭代中，在采油系數(shù)估算步驟中的一些計算中也可能使用所述油藏和驅(qū)油參數(shù)。在該步驟中使用油藏和驅(qū)油參數(shù)來優(yōu)化或?qū)嵤┰撚嬎愫凸浪銓⒁彩潜绢I(lǐng)域技術(shù)人員所能理解的并且也包含在本發(fā)明的范圍之內(nèi)。

接下來根據(jù)總推掃面積和總穿透面積的比較來選擇總油藏采油系數(shù)的計算方法。根據(jù)這些計算中哪個是最高結(jié)果，選擇該總油藏采油系數(shù)合適的計算方法。在選擇了合適的油藏采油系數(shù)計算方法之后，能夠執(zhí)行該計算。

如果總穿透面積小于總推掃面積，那么選擇的基于子區(qū)域的采油系數(shù)計算方法可能會將泡沫區(qū)采油系數(shù)、水區(qū)采油系數(shù)和氣體區(qū)采油系數(shù)進行加和以得到該油藏采油系數(shù)?？商娲?，如果總穿透面積大于總推掃面積，對于三個推掃子區(qū)域中的各區(qū)域，該計算能夠包括：

1.計算最大面積推掃效率；

2.計算穿透之后的異質(zhì)性推掃效率；

3.計算穿透之后的面積推掃效率和穿透之后的異質(zhì)性推掃效率；

4.計算排送效率；以及

5.使用對應(yīng)于該子區(qū)域的排送效率和體積推掃效率計算所述子區(qū)域的子區(qū)域采油系數(shù)。

然后在該情況中，將會把計算得到的所有三個推掃子區(qū)域的子區(qū)域采油系數(shù)進行加和以得到該油藏采油系數(shù)。

在完成了油藏采油系數(shù)的計算之后，能夠?qū)⒂嬎愕玫降挠筒夭捎拖禂?shù)儲存，在該方法的擴展實施方案中，將會使用該系數(shù)繪制出對于驅(qū)油操作中的所有驅(qū)油循環(huán)計算得到的油藏采油系數(shù)在一個圖表軸上的完整圖。

本發(fā)明的方法能夠用于確定單次驅(qū)油循環(huán)或者由單一驅(qū)油循環(huán)構(gòu)成的驅(qū)油操作的采油系數(shù)，或者該方法能夠用于計算出在包括多個驅(qū)油循環(huán)的驅(qū)油操作中被驅(qū)油的油氣藏的采油系數(shù)。本領(lǐng)域技術(shù)人員能夠理解的是，這種方法能夠涵蓋在上下文以及此處基于上下文所述的權(quán)利要求的范圍內(nèi)。

本發(fā)明的方法能夠針對使用任意數(shù)量的注射井和任意數(shù)量的生產(chǎn)井的油氣藏的驅(qū)油操對采油系數(shù)進行建模。

本發(fā)明的泡沫驅(qū)動油藏采油系數(shù)的估算方法能夠手動執(zhí)行或者更可能地以計算機軟件的形式執(zhí)行。該計算機軟件方法能夠使用為特定目的建立的計算機軟件程序、經(jīng)配置以執(zhí)行所需計算的數(shù)學建模軟件、或者甚至是電子數(shù)據(jù)表或類似方法。軟件設(shè)計領(lǐng)域的技術(shù)人員將認識到任意數(shù)量種類的計算機軟件方法，所有這些方法都包含在本發(fā)明的范圍之內(nèi)。

本發(fā)明，一種用于估算泡沫驅(qū)動油氣藏中采油系數(shù)的非暫時性計算機可讀存儲介質(zhì)，實現(xiàn)了其目的，包括含有在被計算機執(zhí)行時實現(xiàn)以下功能的指令的計算機可讀存儲介質(zhì)：將油藏推掃區(qū)定義為油氣藏內(nèi)在至少一個注射井和至少一個生產(chǎn)井之間的體積區(qū)域，其中在至少一個驅(qū)油循環(huán)中注射到該至少一個注射井中的氣體和液體表面活性劑將會把油通過該油藏推掃區(qū)推送至該至少一個生產(chǎn)井以便從中采出。

在定義了油藏推掃區(qū)之后，該計算機可讀存儲介質(zhì)將會使該計算機針對各驅(qū)油循環(huán)、根據(jù)注入該油藏推掃區(qū)中的氣體和液體表面活性劑的當前狀態(tài)確定該油藏推掃區(qū)的油藏和驅(qū)油參數(shù)，并使用該油藏和驅(qū)油參數(shù)，將該油藏推掃區(qū)細分為三個推掃子區(qū)：泡沫推掃子區(qū)、水推掃子區(qū)和氣體推掃子區(qū)。

然后，該計算機可讀存儲介質(zhì)將會使該計算機在計算步驟中：

6.計算各推掃子區(qū)內(nèi)的縱向推掃效率和體積推掃效率；

7.計算該油藏推掃區(qū)的總推掃面積，其是在當前驅(qū)油循環(huán)時被泡沫推掃的油藏推掃區(qū)的體積面積；

8.計算泡沫推掃面積，其是在當前驅(qū)油循環(huán)時被泡沫推掃的泡沫推掃子區(qū)的體積面積；

9.計算水推掃面積，其是在當前驅(qū)油循環(huán)時被水推掃的水推掃子區(qū)的體積面積；

10.計算氣體推掃面積，其是在當前驅(qū)油循環(huán)時被氣體推掃的氣體推掃子區(qū)的體積面積；

11.計算泡沫穿透面積，其是在該驅(qū)油循環(huán)內(nèi)被現(xiàn)有的泡沫穿透的泡沫推掃子區(qū)的體積面積，其包括泡沫推掃面積及其周圍由于泡沫和油之間流動性差異造成的泡沫旁路面積；

12.計算水穿透面積，其是在該驅(qū)油循環(huán)內(nèi)被現(xiàn)有的水穿透的水推掃子區(qū)的體積面積，其包括水推掃面積及其周圍由于水和油之間流動性差異造成的水旁路面積；和

13.計算氣體穿透面積，其是在該驅(qū)油循環(huán)內(nèi)被現(xiàn)有的氣體穿透的氣體推掃子區(qū)的體積面積，其包括氣體推掃面積及其周圍由于氣體和油之間流動性差異造成的氣體旁路面積；以及計算總穿透面積，其是泡沫穿透面積、水穿透面積和氣體穿透面積的總和。

在完成該計算步驟之后，計算機可讀存儲介質(zhì)隨后將會啟動或指令計算機進行采油系數(shù)估算步驟，具體包括以下步驟：

14.計算泡沫區(qū)采油系數(shù)，其是泡沫推掃子區(qū)貢獻的采油系數(shù)；

15.計算水區(qū)采油系數(shù)，其是水推掃子區(qū)貢獻的采油系數(shù)；

16.計算氣體區(qū)采油系數(shù)，其是氣體推掃子區(qū)貢獻的采油系數(shù)；

17.通過比較總推掃面積和總穿透面積的數(shù)值，選擇將要使用的油藏采油系數(shù)計算方法；和

18.執(zhí)行所選擇的油藏采油系數(shù)計算。

該計算機可讀存儲介質(zhì)執(zhí)行的計算和方法的變化將會是在此處描述的采油系數(shù)估算方法的總體范圍內(nèi)能合理推斷的任何改變。

使用計算機軟件執(zhí)行本發(fā)明的方法的一個重要優(yōu)點是在包括根據(jù)計算的油藏采油系數(shù)繪制一個或多個圖的方法的實施方案中，計算了驅(qū)油操作內(nèi)所需的驅(qū)油循環(huán)的油藏采油系數(shù)之后，根據(jù)特定目的建立的軟件或很多可利用的建模工具可以合理地直接繪制該圖。顯然，根據(jù)本發(fā)明的方法的其余部分使用計算機軟件繪制這些圖也包含在本發(fā)明的范圍之內(nèi)，無論采取的特定方法是包括使用集成軟件工具以及本發(fā)明的計算的剩余部分為用戶界 面或打印報表格式繪制一幅或多幅圖，或者是使用安裝在所存儲的油藏采油數(shù)據(jù)上的單獨的繪圖工具繪制該信息的一幅或多幅圖。

附圖說明

為了容易確定對任何特定要素或行為的討論，將參考編號中最重要的一個或多個數(shù)字表示首先介紹該要素時的附圖標記。

圖1是展示依照本發(fā)明的一種實施方案的縱向推掃效率隨油藏和驅(qū)油參數(shù)的變化的圖表；

圖2顯示了油氣藏中的總推掃面積的一種實施例，本發(fā)明的采油系數(shù)估算模型能夠與其對比描述；

圖3是展示在本發(fā)明的泡沫驅(qū)動油藏的采油系數(shù)估算方法的一種實施方案中的步驟的流程圖；

圖4是展示本發(fā)明的泡沫驅(qū)動油藏的采油系數(shù)估算方法的另一實施方案中的步驟的流程圖，包括使用針對多個驅(qū)油循環(huán)存儲的油藏采油系數(shù)繪圖的步驟；

圖5是展示依照本發(fā)明的一種實施方案的體積推掃效率隨油藏和驅(qū)油參數(shù)的變化的圖表；

圖6是依照本發(fā)明的方法的采油系數(shù)估算的一種繪圖結(jié)果；

圖7是依照本發(fā)明的方法的采油系數(shù)估算的另一繪圖結(jié)果；

圖8是依照本發(fā)明的方法的采油系數(shù)估算的另一繪圖結(jié)果；

圖9是依照本發(fā)明的方法的采油系數(shù)估算的另一繪圖結(jié)果；

圖10是依照本發(fā)明的方法的采油系數(shù)估算的另一繪圖結(jié)果；

圖11是依照本發(fā)明的方法的采油系數(shù)估算的另一繪圖結(jié)果；

圖12是依照本發(fā)明的方法的采油系數(shù)估算的另一繪圖結(jié)果；

圖13是依照本發(fā)明的方法的采油系數(shù)估算的另一繪圖結(jié)果；

圖14是依照本發(fā)明的方法的采油系數(shù)估算的另一繪圖結(jié)果；

圖15是依照本發(fā)明的方法的采油系數(shù)估算的另一繪圖結(jié)果；

圖16是依照本發(fā)明的方法的采油系數(shù)估算的另一繪圖結(jié)果；

圖17是依照本發(fā)明的方法的采油系數(shù)估算的另一繪圖結(jié)果；

圖18是展示由計算機可讀存儲介質(zhì)上的處理器指令執(zhí)行的本發(fā)明的方法的一種實施方案的步驟的流程圖；

圖19是圖18的驅(qū)油循環(huán)回路的步驟的分解圖。

具體實施方式

除非另有定義，此處使用的所有技術(shù)和科學術(shù)語都具有與本發(fā)明所屬領(lǐng)域普通技術(shù)人員通常理解的相同的含義。盡管在本發(fā)明的實施或測試中能夠使用與此處所述類似或等效的任何方法和材料，現(xiàn)在描述的是優(yōu)選的方法和材料。

如此處所述，本發(fā)明包括泡沫驅(qū)動油藏的采油系數(shù)的新估算方法。該方法包括下面進一步詳細描述的多個步驟，其可以準確及時地預(yù)測油氣藏中泡沫驅(qū)油項目中的采油系數(shù)。

泡沫驅(qū)油是用于提高油氣藏中的采油系數(shù)的技術(shù)。泡沫驅(qū)油通常由以下操作構(gòu)成：提供油氣藏中的注射井，能夠向其中注射用于油藏的刺激劑，該刺激劑能夠?qū)⒂筒睾蛶r層中包含的油推向生產(chǎn)井，從生產(chǎn)井中能夠采集推掃的油。通常，驅(qū)油操作由一系列驅(qū)油循環(huán)構(gòu)成，其中將把氣體和/或液體表面活性劑注入注射井中以將油通過地層推掃或推動到生產(chǎn)井中。

將氣體和表面活性劑流體注射到地層中通常會導(dǎo)致在地層中產(chǎn)生泡沫。通過在外部壓差下將注射的液體和氣體混合在一起產(chǎn)生泡沫。在驅(qū)油操作的各驅(qū)油循環(huán)中，注射一股液體表面活性劑(例如水、鹽水等)，然后注射一股游離氣體。由于前者的不穩(wěn)定性，液體表面活性劑中的氣體最終將會分離。將氣體和液體表面活性劑循環(huán)注射到油藏中，由此將在油藏中同時存在泡沫、水和氣體。從對現(xiàn)有技術(shù)以及本發(fā)明的方法的一般理解的角度來看，通過注射井將一股液體表面活性劑和/或一股游離氣體的單一配對注射到地層或油藏中稱作一次驅(qū)油循環(huán)。

在驅(qū)油操作之前或在驅(qū)油操作過程中通常需要計算可能的采油系數(shù)或從油氣藏中沖刷出的油氣或?qū)ζ浣?。來自泡沫?qū)油操作的通過一個或多個生產(chǎn)井從油藏中采油提供的結(jié)果通常比不使用泡沫驅(qū)油或刺激從油藏中采集油氣的結(jié)果更好，此處提供的泡沫驅(qū)動油 藏的采油系數(shù)的估算方法提供了一種根據(jù)使用的泡沫驅(qū)油技術(shù)建模或計算來自油氣藏的可能的采油系數(shù)的方法。

通常，在泡沫驅(qū)油應(yīng)用和其他增強采油工程(例如此處所述的那些)中測量來自油氣藏的產(chǎn)量是通過比較油藏或其一些部分的采油系數(shù)而測定的。采油系數(shù)是用于計量不同的增強采油測量的產(chǎn)量對比或者一般而言該油藏生產(chǎn)率的度量。通過建模或通過某些采油活動來提高采油系數(shù)測量值，一般而言是提高該油氣藏的產(chǎn)量的指標。這樣一來，油氣生產(chǎn)商在這種類型的應(yīng)用中需要找到使得在不同應(yīng)用中采油系數(shù)最大化的方法，為此目的，需要能夠在不同情況中對采油系數(shù)建模。本發(fā)明是泡沫驅(qū)動油藏的采油系數(shù)估算方法，即提供為了掌握和/或最大化特定泡沫驅(qū)油工程的產(chǎn)量，在不同泡沫驅(qū)油條件下并隨著時間對采油系數(shù)進行建模的方法。

為了此處所述的方法的目的，將進行泡沫驅(qū)油和將從中采油的泡沫驅(qū)動油藏的區(qū)域稱作油藏推掃區(qū)。該油藏推掃區(qū)是該油氣藏內(nèi)在注射井和生產(chǎn)井之間的區(qū)域，在驅(qū)油操作的至少一個驅(qū)油循環(huán)中將把注射氣體和液體表面活性劑注射到其中。

一般方法概述：

如上所述，本發(fā)明包括一種泡沫驅(qū)動油藏的采油系數(shù)估算方法，其能夠優(yōu)選地在計算機軟件實施方案中實施以進行計算(盡管該方法本身的手動實施也是新穎的且也將落入本發(fā)明的范圍內(nèi))。參照圖3和圖4，我們現(xiàn)在將展示本發(fā)明的方法的兩種實施方案的一般概述，然后我們將對某些樣品數(shù)據(jù)的建模提供更詳細的描述和分析，展示本發(fā)明的實用性。

圖3是展示本發(fā)明的方法的一種實施方案的步驟的流程圖。如本文其他地方所述，本發(fā)明的方法能夠手動進行或以傳統(tǒng)的基于紙張的形式進行，或者更可能地作為軟件系統(tǒng)進行，特別是在存在由大量驅(qū)油循環(huán)構(gòu)成的驅(qū)油操作的情況下。所有這些方法都在本發(fā)明的范圍內(nèi)。

參照圖3，在步驟3-1首先顯示了油藏推掃區(qū)的定義。如上所述，該油藏推掃區(qū)是油氣藏內(nèi)在至少一個注射井到至少一個生產(chǎn)井之間延伸的區(qū)域，為了從油氣藏中推掃油，將把成股的氣體或液體表面活性劑注射到注射井中，通過推掃和驅(qū)油操作將會從生產(chǎn)井中采集和抽取出泡沫、氣體、水或其他液體表面活性劑和油的組合。將回收的液體分離成油 和其他組分超出了本發(fā)明的方法的范圍之外，但本領(lǐng)域技術(shù)人員知曉如何對通過該區(qū)域驅(qū)油操作從油氣藏內(nèi)采集的流體進行分離。

該驅(qū)油操作本身包括通過至少一個注射井將多股氣體和液體表面活性劑(例如水、鹽水等)注射到油氣藏的地層中。每對氣體和液體股注射稱作一個驅(qū)油循環(huán)。整個驅(qū)油操作通常將由多個驅(qū)油循環(huán)構(gòu)成，盡管應(yīng)該理解的是，本發(fā)明的方法能夠用于計算從少至單一驅(qū)油循環(huán)直到包括多達幾十甚至幾百個驅(qū)油循環(huán)的驅(qū)油操作的任何驅(qū)油操作的油藏推掃區(qū)的采油系數(shù)。再次，對該總體驅(qū)油操作的改變或擴展及其對本發(fā)明的泡沫驅(qū)動油藏的驅(qū)油系數(shù)估算方法的數(shù)學函數(shù)的影響也都包含在本發(fā)明的范圍之內(nèi)。

在步驟3-1顯示的定義油藏推掃區(qū)之后，對于該油藏推掃區(qū)和在驅(qū)油操作中各驅(qū)油循環(huán)的各種油藏和驅(qū)油參數(shù)進行一系列計算。在圖3的步驟3-2到步驟3-10中顯示了一個驅(qū)油循環(huán)回路。這意味著對于各個驅(qū)油循環(huán)，都將進行該流程圖的驅(qū)油循環(huán)回路，然后基于此來進行該采油系數(shù)估算方法。

驅(qū)油循環(huán)回路3-2內(nèi)的第一步是定義油藏和驅(qū)油參數(shù)3-3。如本文其他地方所述，油藏和驅(qū)油參數(shù)可包括油藏孔隙率、有助于該油藏內(nèi)不同相間的穿透的因素和用于驅(qū)油操作的其他參數(shù)(例如包括在特定的股或驅(qū)油循環(huán)等中注射的氣體或液體表面活性劑的體積)。在本發(fā)明的方法的執(zhí)行中可能使用的不同類型的油藏和驅(qū)油參數(shù)將是本領(lǐng)域技術(shù)人員知曉的，且在此處被描述，所有這些都包含在本發(fā)明的范圍之內(nèi)。該油藏和驅(qū)油參數(shù)可能隨驅(qū)油循環(huán)而改變或者也可保持恒定。

在定義了特定驅(qū)油循環(huán)的油藏和驅(qū)油參數(shù)之后，該油藏推掃區(qū)將被細分為多個推掃子區(qū)，即泡沫推掃子區(qū)、氣體推掃子區(qū)和水推掃子區(qū)。這些推掃子區(qū)的定義將根據(jù)知曉的油藏和驅(qū)油參數(shù)以及該油藏內(nèi)的驅(qū)油操作的一般狀態(tài)而計算或確定。有效地，在驅(qū)油操作中，該泡沫推掃子區(qū)、氣體推掃子區(qū)和水推掃子區(qū)將移動通過該油藏推掃區(qū)的地層，從至少一個注射井移動到至少一個生產(chǎn)井。這些推掃子區(qū)的移動速度、前端形狀和大小將都受到各種油藏和驅(qū)油參數(shù)的影響，再次，還在本發(fā)明的范圍內(nèi)提供了這些推掃子區(qū)的建模，無論其采用何種特定方法。

在將油藏推掃區(qū)細分成這三個推掃子區(qū)之后，在步驟3-5所述的該方法的下一步是計算各推掃子區(qū)的縱向推掃效率和體積推掃效率。這些計算的細節(jié)在本文的其他地方描述。除了在步驟3-5所示計算縱向推掃效率和體積推掃效率之外，在步驟3-6還顯示了執(zhí)行其他 計算，這是為了估算在特定驅(qū)油操作或驅(qū)油循環(huán)中各推掃子區(qū)對總油藏采油系數(shù)的采油系數(shù)貢獻所需要的。存在一些特定的模型，其被認為最適用于估算各推掃子區(qū)的采油系數(shù)或其對總油藏的采油系數(shù)的貢獻。這顯示在步驟3-7中，在計算了各推掃子區(qū)的各個采油系數(shù)估算值(泡沫區(qū)采油系數(shù)、氣體區(qū)采油系數(shù)和水區(qū)采油系數(shù))之后，能夠確定總油藏采油系數(shù)。這顯示在步驟3-8。

如說明書中詳細描述的那樣，如步驟3-8所示，需要選擇特定的計算方法來確定該油藏采油系數(shù)。所選擇的計算方法將在很大程度上根據(jù)該油藏推掃區(qū)中各推掃子區(qū)之間的穿透系數(shù)的計算值來確定。在3-8選擇了合適的油藏采油系數(shù)計算方法之后，將執(zhí)行該計算，得到油藏采油系數(shù)的計算值。最后，如步驟3-9所示，可儲存該油藏采油系數(shù)的計算值，接著步驟3-10顯示的是驅(qū)油循環(huán)回路的終點。然后可對該驅(qū)油操作內(nèi)的其他驅(qū)油循環(huán)之行該驅(qū)油循環(huán)回路，直至達到該特定驅(qū)油操作的最大驅(qū)油循環(huán)數(shù)，對于各驅(qū)油循環(huán)，該方法中驅(qū)油循環(huán)的增量(計數(shù))將會增加。

如上所述，本發(fā)明的方法特別地將能夠使用計算機軟件執(zhí)行。存在各種類型的數(shù)學建模軟件能夠定制或編程以用于執(zhí)行本發(fā)明的計算方法，或者該方法能夠使用根據(jù)特定目的建立的計算機軟件而實施，如本文其他地方所述。

現(xiàn)在參照圖4，該圖顯示了對本發(fā)明的計算方法略微更詳細描述的另一實施方案的另一流程圖。在圖4的方法開始時，在步驟4-1顯示了再次定義油藏推掃區(qū)。在定義了油藏推掃區(qū)之后，開始驅(qū)油循環(huán)回路，在本例中其從步驟4-2延伸到步驟4-13，表示在驅(qū)油操作中各選定的驅(qū)油循環(huán)的同組計算的執(zhí)行。驅(qū)油循環(huán)回路4-2中的第一步是定義油藏和驅(qū)油參數(shù)4-3。上面所述的一些油藏和驅(qū)油參數(shù)在整個驅(qū)油操作過程中將保持恒定，其他的將不間斷地改變，且對于每次驅(qū)油循環(huán)都將再次計算。執(zhí)行本發(fā)明的方法所需的油藏和驅(qū)油參數(shù)在此處被專門描述，能夠以這種方式使用的所有油藏和驅(qū)油參數(shù)都包含在本發(fā)明的范圍內(nèi)。

在針對驅(qū)油操作中的特定驅(qū)油循環(huán)捕獲、定義或改進油藏和驅(qū)油參數(shù)之后，再次將該油藏推掃區(qū)分割成3個推掃子區(qū)，即泡沫推掃子區(qū)、水推掃子區(qū)和氣體推掃子區(qū)。如步驟4-5中顯示，計算泡沫推掃子區(qū)、水推掃子區(qū)和氣體推掃子區(qū)的各子區(qū)的縱向推掃效率和體積推掃效率，方法移動通過到額外的面積計算步驟4-6，在步驟4-7顯示了在油藏推掃區(qū)內(nèi)各推掃子區(qū)的采油系數(shù)計算。

步驟4-8再次顯示了對油藏采油系數(shù)計算方法的選擇，更詳細地用決定框4-9顯示，該油藏采油系數(shù)計算方法在一種實施方案中(例如如該圖中所示)是通過比較推掃子區(qū)組合的總穿透面積與油藏推掃區(qū)的總推掃面積而確定的。如果該穿透面積大于總推掃面積，那么步驟4-11顯示了根據(jù)基于各推掃子區(qū)的計算方法選擇和執(zhí)行油藏采油系數(shù)的計算方法?；蛘?，如果總穿透面積不大于總推掃面積，能夠使用一個更基本的集合方程，簡單地將泡沫區(qū)采油系數(shù)、氣體區(qū)采油系數(shù)和水區(qū)采油系數(shù)加在一起。在執(zhí)行了采油系數(shù)計算方法或計算之后，在步驟4-12顯示，能夠存儲計算出的油藏采油系數(shù)。這表示驅(qū)油循環(huán)回路的結(jié)束，對于該驅(qū)油操作中各選定的驅(qū)油循環(huán)進行整個回路。然后能夠?qū)τ谠擈?qū)油操作中的其他驅(qū)油循環(huán)執(zhí)行該驅(qū)油循環(huán)回路，直至對于特定的驅(qū)油操作達到最大的驅(qū)油循環(huán)數(shù)，對于各驅(qū)油循環(huán)，該方法中驅(qū)油循環(huán)的增量(計數(shù))將會增加。

在圖4的方法中還在步驟4-14顯示了根據(jù)存儲的油藏采油系數(shù)的計算值繪制圖表。有效地，可以將各驅(qū)油循環(huán)的油藏采油系數(shù)計算值繪在圖的一個軸上，相對于時間或其他油藏和驅(qū)油參數(shù)等繪圖，以提供不同程度和類型的分析工具用于建模、執(zhí)行和理解特定的驅(qū)油操作。應(yīng)該認識到，步驟4-14顯示的被繪制的圖的特定類型可以發(fā)生變化。任何類型的這種圖都包含在本發(fā)明的范圍內(nèi)，包括圖6-圖17中顯示的那些。

油和泡沫在油藏內(nèi)的推掃效率和行為：

在包括將液體表面活性劑和游離氣體注射到油氣藏中的驅(qū)油操作中，假設(shè)氣體、液體表面活性劑以及由在地層內(nèi)對該氣體和液體表面活性劑施加外部壓力產(chǎn)生的泡沫，所有三種將同時在地層中都存在。由于其流度的顯著差異，假設(shè)在油藏內(nèi)將會有三個注射流體的推掃子區(qū)，包括泡沫推掃子區(qū)、水推掃子區(qū)和氣體推掃子區(qū)。通過假設(shè)在其之間沒有重疊，所有三個推掃子區(qū)的面積之和被認為是該油藏推掃區(qū)的總推掃面積。由于從泡沫、液體到氣體，流度逐漸增大，因此從注射井到生產(chǎn)井這三個推掃子區(qū)的順序是泡沫推掃子區(qū)、水推掃子區(qū)和氣體推掃子區(qū)，氣體前端將要比水或液體表面活性劑(二者是由泡沫的坍塌物產(chǎn)生的)更早到達生產(chǎn)井。如果每次完整的泡沫股(包括一個水股和一個氣體股)被認為是一個驅(qū)油循環(huán)，那么在氣體或水的前端到達生產(chǎn)井時能夠計算驅(qū)油循環(huán)數(shù)。一旦泡沫驅(qū)油的工作計劃決定了，那么能夠根據(jù)水或氣體股大小來估算特定的穿透時間。

本發(fā)明的方法以及建模模擬要進行驅(qū)油操作的油氣藏內(nèi)的采油系數(shù)的關(guān)鍵在于了解地層地質(zhì)內(nèi)的推掃效率的概念。將氣體和流體注入地層中形成泡沫，隨著其通過地層朝向生產(chǎn)位置移動時，使油從該地層被推掃或推向生產(chǎn)位置。

在氣體前端到達生產(chǎn)井之前，假設(shè)累積生產(chǎn)量與累積時間和累積注射體積線性相關(guān)。該時間段內(nèi)的推掃面積計算為與總注射流體體積之比為1:1。在氣體穿透之后，用穿透非均質(zhì)推掃效率的改進方程計算推掃效率，直至水的前端到達生產(chǎn)井。然后在水到達生產(chǎn)井之后，用于液體的另一方程被用于計算推掃效率，引入理論最大推掃效率來設(shè)定曲線在穿透后的限值。

為了確定對推掃效率的影響以及為了確定推掃效率隨著與一個或多個驅(qū)油循環(huán)相關(guān)的時間或注射體積的變化，三種推掃子區(qū)各自的穿透時間和理論最大推掃效率是重要的變量。由于這三種推掃子區(qū)的假設(shè)和泡沫的不穩(wěn)定性，在地層內(nèi)氣體和水的排送前端將會在泡沫之前到達生產(chǎn)井。在效率隨著時間或注射體積而改變的繪圖中兩個更關(guān)鍵的點包括氣體前端的穿透時間和水前端的穿透時間。所有這些特征的繪圖被這兩個關(guān)鍵點分割成三個部分。任何注射流體穿透之前的部分像是斜線。一旦氣體前端到達生產(chǎn)井，其將會轉(zhuǎn)變成曲線。在一段時期之后，一旦水前端到達生產(chǎn)井，曲線的斜率將顯著增大。圖1是依照本發(fā)明的一種實施方案的這些特征的樣品繪圖。

定義油藏中的關(guān)鍵面積：

如上所述，油氣藏的泡沫驅(qū)油在物理上包括產(chǎn)生或使用至少一個注射井和至少一個生產(chǎn)井。該至少一個注射井用于將游離氣體和液體表面活性劑注入油氣藏的地層中，其在地質(zhì)學上位于合適的位置，使得注入的游離氣體和液體表面活性劑將移動通過該油氣藏的地質(zhì)層移動到至少一個生產(chǎn)井，將油從該地層朝向生產(chǎn)井推掃或推動，在該生產(chǎn)井的位置處能夠?qū)⒂筒沙觥４颂幙紤]的大多數(shù)建模和實施方案都基于僅一個注射井和一個生產(chǎn)井的數(shù)學模型，以便對推掃效率和其他變量進行建模。然而，油氣藏建模和增強采油技術(shù)(例如此處所述那些)領(lǐng)域的技術(shù)人員將認識到，在不脫離本發(fā)明的范圍的情況下，依照本發(fā)明的方法的剩余部分也能夠接受和模擬包括多于一個注射井和/或多于一個生產(chǎn)井的其他實施方案。所有這些方法都包含在本發(fā)明的范圍之內(nèi)。本領(lǐng)域技術(shù)人員將能夠理解注射井、生產(chǎn)井在地質(zhì)學定義中的概念以及在泡沫驅(qū)油工程中注射生產(chǎn)位置的設(shè)置。

除了在油氣藏的泡沫驅(qū)油技術(shù)中總體了解注射位置和生產(chǎn)位置的概念之外，通常在例如此處所述的模型和方法中對采油系數(shù)進行建模的本發(fā)明的方法的首要步驟之一是定義油藏推掃區(qū)。油藏推掃區(qū)是至少一個注射井和至少一個生產(chǎn)井之間的油氣藏的體積或面積，其中注射的氣體和液體表面活性劑將在驅(qū)油操作的至少一個驅(qū)油循環(huán)中被注入。有效地，該油藏推掃區(qū)是需要使用泡沫驅(qū)油技術(shù)從其中推掃油的油氣藏的區(qū)域。定義油藏推掃區(qū)的概念將被增強采油技術(shù)領(lǐng)域技術(shù)人員所理解，因為針對驅(qū)油工程定義所需的推掃區(qū)域的這種相同類型的概念是在特定驅(qū)油工程的建模方面將會執(zhí)行的第一框架活動之一。通過為至少一個注射井和至少一個生產(chǎn)井選擇位置，在一定程度上定義該油藏推掃區(qū)，再次，在涉及特定驅(qū)油工程的適合范圍方面，這將被增強采油技術(shù)領(lǐng)域技術(shù)人員所理解。

在油氣藏內(nèi)以及在油藏推掃區(qū)內(nèi)，該油氣藏的油藏推掃區(qū)的采油因素的計算或建模的前期步驟之一是確定將該油藏推掃區(qū)分割成該油藏推掃區(qū)內(nèi)的三個推掃子區(qū)。通過將油藏推掃區(qū)分割成泡沫推掃子區(qū)、水推掃子區(qū)和氣體推掃子區(qū)，在其各個子區(qū)中能夠完成獨立的推掃效率和其他變量的計算，據(jù)信本發(fā)明的方法的精度和速度與現(xiàn)有技術(shù)方法相比有顯著的提高。對油藏推掃區(qū)內(nèi)推掃子區(qū)的定義將會因與該油藏推掃區(qū)相關(guān)的油藏和驅(qū)油參數(shù)的認識而被教導(dǎo)。油氣地質(zhì)領(lǐng)域技術(shù)人員將再次認識到，能夠采用不同的方法具體限定該推掃子區(qū)，盡管通常而言將認識到在油藏推掃區(qū)中對這些區(qū)域的定義將通過對該油藏中氣體、流體和泡沫組分的行為以及將在該特定地質(zhì)內(nèi)在將油推掃到生產(chǎn)井時的行為方式的了解而被教導(dǎo)。

為了提高本方法中建模的準確度，根據(jù)對該油藏推掃區(qū)內(nèi)的泡沫、水和氣體在相鄰的推掃子區(qū)之間的穿透的了解，本發(fā)明的方法的某些實施方案可通過引入泡沫穿透面積、水穿透面積和氣體穿透面積而進一步調(diào)節(jié)地層內(nèi)采油系數(shù)和推掃效率的計算，在泡沫穿透面積、水穿透面積和氣體穿透面積內(nèi)可以進一步改進或改變該行為和局部采油系數(shù)。

圖2顯示了依照本發(fā)明的油藏推掃區(qū)的一種實施方案，且描述了其中的三個推掃子區(qū)的分布。假設(shè)這三種推掃子區(qū)由于其流度的差異而具有相對清晰的界面，圖中的afoam、awater和agas表示泡沫推掃子區(qū)、水推掃子區(qū)和氣體推掃子區(qū)，而周圍的區(qū)域abtfoam、abtwater和abtgas表示泡沫穿透區(qū)域、水穿透區(qū)域和氣體穿透區(qū)域，表示總推掃區(qū)域能夠在特定的驅(qū)油循環(huán)中被泡沫、水和氣體穿透的程度。

盡管不需要假設(shè)區(qū)域abtfoam，但在流度極低的情況下，泡沫可能會穿透某個區(qū)域，此時abtfoam幾乎等于afoam。這種潛在的機理被認為是為了使整個計算過程更均勻。使用圖3中顯示的進一步的簡化來強調(diào)整個過程被分成三個時期的方式。例如，在計算出水推掃的面積的同時，計算從區(qū)域右側(cè)外部abtfoam開始，因為泡沫前側(cè)的水將不會推掃剛被泡沫穿透的區(qū)域。該方法用于計算各區(qū)的推掃面積和穿透面積。

圖2還顯示了注射井和生產(chǎn)井的位置，以提供對此處所述的泡沫驅(qū)動油藏的采油系數(shù)估算方法的油層推掃區(qū)體積面積和總體操作的進一步構(gòu)思。

本說明書其他部分展示了本發(fā)明的采油系數(shù)建模方法在不同參數(shù)下執(zhí)行的各種迭代，對圖2中顯示的推掃子區(qū)和穿透區(qū)域的定義的理解被認為對以下描述是非常重要的。

關(guān)鍵假設(shè)：

通過在外部壓力差下將注射的液體表面活性劑和氣體混合在一起產(chǎn)生泡沫。將氣體和液體表面活性劑各自注入油氣藏中構(gòu)成一個驅(qū)油循環(huán)。對于各驅(qū)油循環(huán)，將一股表面活性劑溶液注入，然后注入一股游離氣體。由于泡沫的不穩(wěn)定性，氣體和液體表面活性劑最終將會分開。將氣體和液體表面活性劑循環(huán)注入該油藏中，造成泡沫、水和氣體同時存在。由于其在流度上的顯著差異，假設(shè)在油氣藏內(nèi)將會有注射流體的三個推掃子區(qū)，包括泡沫推掃子區(qū)、水推掃子區(qū)和氣體推掃子區(qū)。通過假設(shè)在其間沒有重疊，所有三個推掃子區(qū)的總和被認為是油藏推掃區(qū)。

以下是在此處所示的建模實施方案和實施例中用到的關(guān)鍵假設(shè)：

·在注射一股表面活性劑然后注射另一股游離氣體之后，在該位置產(chǎn)生泡沫；

·隨著由半衰期描述的時間的經(jīng)過，泡沫坍塌；

·一旦其坍塌，泡沫不會再生；

·時間是決定剩余泡沫體積的唯一因素。不考慮接觸油而弱化泡沫強度的影響；

·泡沫、游離氣體和水(液體表面活性劑)在流度上具有相對明顯的差異；

·具有較高流度的流體行進至注射流體的前端，其將不會到達已經(jīng)被較低流度的流體推掃過的區(qū)域；

·考慮了兩種穿透情況，其中一種是雙前端情況。由于流度的差異，分別估算氣體和水的穿透時間。第一前端是由氣竄(gaschanneling)造成的。第二前端是由氣體和水 (液體表面活性劑)之間的流度差異造成的，其能夠造成產(chǎn)量的少量提高，因為水的驅(qū)動具有更高的推掃效率。另一種情況是單前端情況，其源自傳統(tǒng)水驅(qū)動的情況。

本領(lǐng)域技術(shù)人員將認識到，本發(fā)明也能夠用對這些假設(shè)的改進來實施，本發(fā)明人認為不脫離預(yù)期范圍的對建模假設(shè)的所有這些改變都在本發(fā)明的范圍之內(nèi)。

術(shù)語命名：

以下術(shù)語命名用于此處所述的公式中，且為了參考的目的所述術(shù)語被提供在示例的實施方案以及此處所述各種計算中呈現(xiàn)的概念和測量值或數(shù)值的展示中：

abtfoam＝泡沫穿透面積：被現(xiàn)有的泡沫穿透的面積，m²

abtgas＝氣體穿透面積：被現(xiàn)有的氣體穿透的面積，m²

abtwater＝水穿透面積：被現(xiàn)有的水或液體表面活性劑穿透的面積，m²

afoam＝泡沫推掃面積，m²

agas＝氣體推掃面積，m²

awater＝水推掃面積，m²

cgasinfoam＝氣體體積因子，在泡沫內(nèi)部捕獲的氣體體積與泡沫外部相同摩爾氣體的體積之比

eaafterbt＝穿透后的面積推掃效率

eabt＝特定流度比的穿透后的面積推掃效率

eafoam＝泡沫的面積推掃效率

eagas＝氣體的面積氣體推掃效率

eamax＝具有特定流度比和滲透率變化的理論最大推掃效率

eawater＝水的面積推掃效率

edfoam＝泡沫的排送效率

edgas＝氣體的排送效率

edwater＝水的排送效率

evfoam＝泡沫的體積推掃效率

evgas＝氣體的體積推掃效率

evwater＝水的體積推掃效率

ezfoam＝泡沫的縱向推掃效率

ezgas＝氣體的縱向推掃效率

ezwater＝水的縱向推掃效率

fw＝含水率

h＝油藏高度,m

m＝流度比

m*＝虛擬流度比

pv＝注射的孔體積

rf＝油藏采油系數(shù)

rffoam＝泡沫區(qū)采油系數(shù)

rfgas＝氣體區(qū)采油系數(shù)

rfwater＝水區(qū)采油系數(shù)

soi＝初始油飽和率

sor＝殘余油飽和率

swc＝原始水飽和率

t＝半衰期，泡沫的存在量降低到一半的時間，h

t＝產(chǎn)生的泡沫在特定循環(huán)之后的經(jīng)過時間，h

ut＝將流體注射到油藏內(nèi)的速度，m/天

ugas＝氣體注射速率(sc)，m³/d

uwater＝液體表面活性劑注射速率(sc)，m³/d

vfoam＝泡沫推掃的孔體積，m³

vgas＝氣體推掃的孔體積，m³

vwater＝水推掃的孔體積，m³

v＝滲透率變化

vi＝一個注射循環(huán)之后新形成的泡沫的體積，m³

vleft＝由注射流體循環(huán)產(chǎn)生的泡沫在時間t之后的剩余體積，m³

vtotalfoam＝剩余泡沫總量，m³

winj＝已經(jīng)注射的流體的總體積(sc)，m³

wibt＝在發(fā)生穿透時已經(jīng)注射的流體的總體積(sc)，m³

μo＝油粘度，cp

μw＝水粘度，cp

φ＝孔隙率。

方法計算詳述：

為了更詳細描述的目的，現(xiàn)在我們將展示在此處所述的泡沫驅(qū)動油藏的采油系數(shù)估算方法的某些實施方案中所用的數(shù)學方法。以下用于展示依照本發(fā)明的總體方法使用的數(shù)學模型的一些實施方案的計算將被理解為能夠被改進以改善或優(yōu)化其結(jié)果，或者通過改變其中的某些參數(shù)或條件以產(chǎn)生在本發(fā)明的總體預(yù)期范圍內(nèi)的特別結(jié)果，本領(lǐng)域技術(shù)人員將認識到的所有這些改變都包含在本發(fā)明的范圍內(nèi)。

總穿透面積：

在驅(qū)油循環(huán)回路的迭代中將要進行的特定驅(qū)油循環(huán)的第一組計算是關(guān)于總推掃面積和總穿透面積的計算，其需要計算水穿透面積、氣體穿透面積和泡沫穿透面積。泡沫、水和氣體在油藏推掃區(qū)中相應(yīng)的推掃面積之間穿透的穿透時間是相關(guān)聯(lián)的計算。氣體與水相比具有更高的流動性，水與油相比具有更高的流動性。氣體推掃區(qū)比水推掃區(qū)距該至少一個注射井更遠，水推掃區(qū)比泡沫推掃區(qū)距該注射井更遠。由于這種流動性比例的特別差異，氣體將在水/表面活性劑和泡沫之前到達該至少一個生產(chǎn)井。下式使用上面所述的術(shù)語命名，其顯示了氣體穿透到至少一個生產(chǎn)井時的情況：

abtfoam+abtwater+abtgas=atotal(1)

atotal是總推掃面積，在計算開始時通過數(shù)據(jù)輸入或其他方式對其進行具體指定。依照圖2，一旦實現(xiàn)該方程，就意味著abtfoam、abtwater和abtgas的總和達到atotal，意味著沒有更多的新區(qū)域供氣體推掃，氣體前端到達生產(chǎn)井。此時，abtfoam、abtwater和abtgas能夠根據(jù)以下方程計算：

在五點工作計劃(five-spotworkingscheme)中，c＝0.718。如果該工作計劃是反九點，則c＝0.525，在反七點工作計劃中，c＝0.743。m*是虛擬流度比，其是非均質(zhì)情況的修正的流度比。該流度比是由流度比(m)和滲透率變化(v)確定的。當v<0.7時，

當v>0.7時，

m表示均質(zhì)情況的流度比，v是滲透率變化。afoam、awater和agas是由方程(11)、(17)和(22)計算的。m*foam是泡沫與油的虛擬流度比；m*water是水與油的虛擬流度比；且m*gas是氣體與油的虛擬流度比。

水推掃子區(qū)中的穿透時間：

水穿透進入氣體推掃子區(qū)中將發(fā)生在氣體穿透之后。特別地，當以下方程滿足時，將在水推掃子區(qū)內(nèi)發(fā)生穿透：

abtfoam+abtwater＝atotal(7)

atotal是總油藏圖案面積。參照圖2，一旦實現(xiàn)該方程，意味著abtfoam和abtwater的總和達到atotal，表明沒有更多的新區(qū)域供水推掃，水前端到達至少一個生產(chǎn)井。

估算穿透之前的采油系數(shù)

如果abtfoam、abtwater和abtgas之和小于atotal，該系統(tǒng)被認為尚未穿透。以下是在采油系數(shù)的評價中用到的接下來的計算。

泡沫推掃面積：

因為成股地將游離氣體和液體表面活性劑注射到油氣藏中，假設(shè)產(chǎn)生的落入油藏推掃區(qū)內(nèi)的體積等于注射氣體和注射表面活性劑液體在油藏壓力時的總和?？紤]到不是所有注射到油藏中的氣體都能與表面活性劑成分形成泡沫，因此引入了一個系數(shù)r來描述氣體過度注射時的情況。r是在同一股中形成泡沫的氣體與注射的氣體總量之比。在計算中，r被認為是常系數(shù)。注射但未產(chǎn)生泡沫的氣體假設(shè)與來自泡沫破裂坍塌產(chǎn)生的氣體混合。二者將結(jié)合在一起并形成氣竄的前端，這也是前面提及的氣體推掃子區(qū)。與相同分子量的游離氣體相比，泡沫內(nèi)的氣體被壓縮。關(guān)于理想氣體：

pbubblev＝nrt。

同時，每個單一氣泡中的氣體壓力與氣泡外的氣體壓力具有如下關(guān)系：

γ是氣泡的表面張力。盡管表面活性劑在開始時降低了表面張力，但表面膨脹最終將降低表面活性劑在表面上的濃度，使γ增大。然后，pbubble增大，氣泡被壓縮。

引入另一個新的系數(shù)，氣泡中的氣體體積系數(shù)cgasinfoam，來描述泡沫內(nèi)捕獲的氣體體積與泡沫外相同摩爾氣體的體積之比，其僅比1小一點。這一系數(shù)被引入用于更精確地描述泡沫的存在。與其他類型的注射流體相比，除了其獨特的剖面控制性能之外，泡沫的不穩(wěn)定性也是泡沫的一個顯著特征。引入泡沫的半衰期來計算油藏中泡沫的剩余體積。

vleft是由一個注射流體循環(huán)產(chǎn)生的泡沫在時間t之后的剩余體積。vi是在一個注射循環(huán)之后新形成的泡沫的體積。t表示半衰期，其是由選定的泡沫或發(fā)泡劑所決定的，t是 在泡沫產(chǎn)生之后所經(jīng)過的時間。該方程用于計算一次循環(huán)注射產(chǎn)生的泡沫的剩余體積。依照循環(huán)注射方法，通過將每次循環(huán)的剩余體積加在一起就能夠計算剩余泡沫的總量。泡沫產(chǎn)生時間假設(shè)為一個完整循環(huán)的兩個股都完成注射時的時間。剩余泡沫的總量是由以下方程計算的：

然后，能夠用以下公式計算泡沫推掃體積vfoam：

swc表示水飽和率。sor是殘余油飽和率，其等于在油的滲透率等于0時在泡沫油相對滲透率曲線中的端點的飽和率。泡沫驅(qū)油區(qū)的推掃面積afoam等于：

泡沫推掃子區(qū)的體積推掃效率：

依照圖2，具有良好剖面控制性能的泡沫推掃區(qū)比油藏推掃區(qū)要小得多。該泡沫推掃子區(qū)的前端將幾乎不能到達該至少一個生產(chǎn)井。假設(shè)注射的流體體積等于推掃的油藏流體的體積，那么該推掃的流體體積等于以下方程，

evfoam＝ezfoamxeafoam

eafoam表示該泡沫推掃子區(qū)內(nèi)的縱向推掃效率，其等于：

ezfoam等于縱向推掃效率，其是依照craig的實驗?zāi)Ｐ?1957)計算的。該模型提供了用不同的流度比和粘性力與重力比計算縱向推掃效率的方法。

下面顯示了粘性力與重力比的定義：

在該式中，ut是油藏內(nèi)注射流體的速率，以m/天計。μo是油的粘度，以pa*s計。x是該至少一個注射井與至少一個生產(chǎn)井之間的距離。k是水平滲透率，以md計。δρ是注射流體和油之間的密度差。從現(xiàn)有方案中所取的顯示該關(guān)系的擬合曲線方程為：

m表示泡沫和油之間的流度比。該方程能夠用于計算泡沫推掃子區(qū)內(nèi)的縱向推掃效率。

泡沫推掃子區(qū)內(nèi)的排送效率：

在輕質(zhì)油的油氣藏中，泡沫具有比鹽水或油高得多的粘度，且用作調(diào)剖劑，為此原因，假設(shè)泡沫驅(qū)油的運動是活塞式的排送。排送效率被認為在該假設(shè)中是恒定的，其等于下面的方程，其中soi是初始油飽和率。sor是殘余油飽和率。它們是從相對滲透率曲線的端點值獲得的：

依照其定義，該泡沫區(qū)采油系數(shù)等于由體積推掃效率和排送效率使用例如以下方程相乘得到的乘積：

rffoam＝edfoam×evfoam(16)

水推掃子區(qū)：

為了實現(xiàn)本發(fā)明的方法的下一個面積和計算涉及水推掃子區(qū)。該水推掃區(qū)來自泡沫的坍塌，位于泡沫推掃子區(qū)的前面，其中λsurfactant表示在每次注射循環(huán)中表面活性劑溶液的體積比。vi是每次注射循環(huán)的總注射流體體積：

該水推掃子區(qū)的體積推掃效率計算為：

evwater＝ezwater×eawater(18)

eawater是該水推掃子區(qū)的體積推掃效率：

ezwater是該水推掃子區(qū)的縱向推掃效率。水推掃子區(qū)的縱向推掃效率的計算與泡沫推掃子區(qū)的縱向推掃效率的計算相同，其也是依照craig的模型用不同的流度比計算的。

水推掃子區(qū)中的排送效率

在穿透之前，所有注射體積都對油氣藏內(nèi)新區(qū)域的采油有貢獻。該前端之前的被排送流體的平均飽和度保持恒定，直至注射流體穿透。為此原因，排送效率也保持恒定。

在該計算中，soi是初始油飽和率，sor是殘余油飽和率。它們是從相對滲透率曲線的端點值獲得的。

水推掃子區(qū)中的采油系數(shù)：

依照其定義，采油系數(shù)等于油的體積或體積推掃效率與排送效率相乘的乘積：

rfwater＝edwater×evwater(21)

氣體推掃子區(qū)：

在氣體推掃子區(qū)中使用同一組的計算。首先，需要計算氣體推掃面積，其來自泡沫的坍塌，并且由于流度差異而位于水推掃子區(qū)之前。

在該計算中，agas是該循環(huán)中氣體的體積與注入的所有流體的體積比。c是考慮反應(yīng)和氣體消耗的系數(shù)。例如，如果選擇空氣泡沫注射，c等于0.9496，表示在低溫氧化工藝過程中產(chǎn)生的co2的吸收。

氣體推掃子區(qū)中的體積推掃效率等于evgas：

evgas＝ezgas×eagas(23)

eagas是氣體推掃子區(qū)中的推掃效率：

ezgas是氣體推掃子區(qū)中的縱向推掃效率。氣體推掃子區(qū)的縱向推掃效率的計算與泡沫推掃子區(qū)中縱向推掃效率的計算相同，也是依照craig的模型用不同的流度比計算的。

在穿透之前，所有注射體積都對油藏內(nèi)新區(qū)域的采油有貢獻。該前端之前被排送的流體的平均飽和度保持恒定，直至注射流體穿透。為此原因，排送效率也保持恒定，其中soi是初始油飽和率，sor是殘余油飽和率：

氣體推掃子區(qū)接下來要計算的是氣體區(qū)采油系數(shù)。依照其定義，采油系數(shù)等于體積推掃效率與排送效率相乘的乘積：

rfgas＝edgas×evgas(26)

依照各區(qū)的采油系數(shù)能夠計算出油藏采油系數(shù)rf：

rf＝rffoam+rfwater+rfgas(27)

如果abtfoam、abtwater和abtgas之和大于atotal，該系統(tǒng)被認為已經(jīng)穿透：

abtfoam+abtwater+abtgas＞atotal

在氣體已經(jīng)穿透但表面活性劑溶液未穿透的情況中，將使用來自均質(zhì)情況的經(jīng)驗方程估算穿透后最大面積推掃效率來計算體積推掃效率。那么，使用圖2將原始方程改變?yōu)榭紤]非均質(zhì)的情況。對于均質(zhì)的油藏，穿透時的最大面積推掃效率能夠計算為：

c是關(guān)于井網(wǎng)分布的系數(shù)。當工作計劃是五點時，c＝0.718。當工作計劃是反九點時，c＝0.525，或者當工作計劃是反七點時，c＝0.743。

對于均質(zhì)的油藏，使用下面的方程計算穿透后的最大面積推掃效率：

其中d是關(guān)于井網(wǎng)分布的系數(shù)。當工作計劃是五點時，d＝0.2749。當工作計劃是反九點時，d＝0.201。

對于非均質(zhì)的油藏，根據(jù)圖案改變結(jié)果。最大面積推掃效率計算的改變聚焦在兩點：穿透時的最大面積推掃效率和特定注射流體的理論最大面積推掃效率。依照非均質(zhì)的油藏的推掃效率計算有關(guān)的參考文獻，選擇滲透率變化(v)作為變量以展現(xiàn)非均質(zhì)的影響。因為現(xiàn)有的穿透時最大面積推掃效率的方程僅具有一個變量(即流度比)，因此得到一個中間變量是關(guān)鍵的。該中間變量能夠?qū)⒘鞫缺群蜐B透性變化與推掃效率聯(lián)系在一起。根據(jù)圖案，隨著滲透率變化增大，最大面積推掃效率降低。同時，隨著流度比的增大，最大面積推掃效率也降低。一旦決定了流度比和滲透率變化，從圖案中得到的最大面積推掃效率等于從均質(zhì)的油藏中的其他更高的流度得到的最大面積推掃效率。這種更高的流度比是前面所述的中間變量，將其稱作虛擬流度比。該虛擬流度比定義為在特定的滲透率變化時能夠用其他流度比得到相同的最大面積推掃效率的均質(zhì)油藏推掃區(qū)中的流度比。

依照該物理解釋，非均質(zhì)油藏推掃區(qū)的滲透率變化總是大于0，因此虛擬流度總是高于原始流度?？紤]了非均質(zhì)之后的穿透時最大面積推掃效率方程示于下面。m*為虛擬流度比。

其中c是關(guān)于井網(wǎng)分布的系數(shù)。當工作計劃是五點時，c＝0.718。當工作計劃是反九點時，c＝0.525，或者當工作計劃是反七點時，c＝0.743。

在穿透后，使用下面的方程計算非均質(zhì)油藏推掃區(qū)中穿透后的最大面積推掃效率：

當工作計劃是五點時，d＝0.2749。當工作計劃是反九點時，d＝0.201。x是設(shè)定在當前的流度比和滲透率變化時趨于最大面積推掃效率eamax的最大面積推掃效率的系數(shù)。

圖2的最大面積推掃效率是當wor等于25或96.15％水百分比時的實驗結(jié)果。假設(shè)此時至少一個生產(chǎn)井停止，考慮最大面積推掃效率。ea(m,v)能夠分成兩部分進行曲線擬合，在較低滲透率變化時的平滑部分l1，隨著非均質(zhì)性的增大，最大面積推掃效率變化不大；和在較高的相對低滲透率時的斜線部分l2，其隨著非均質(zhì)性的較小改變而顯著變化。這兩部分的分割點的位置不是恒定的，該位置是流度比的函數(shù)。

vc＝-0.063×ln(m+0.1)+0.5473(32)

m是初始流度比。vc是光滑部分l1和斜線部分l2之間的分割點。依照曲線擬合結(jié)果，當滲透率變化v小于vc時，該最大面積推掃效率為以下方程：

eamax＝-0.07076v+0.8606-0.0605×ln(m+0.1)(33)

當滲透率變化v大于或等于vc時，該油藏推掃區(qū)的最大面積推掃效率為以下方程：

然后根據(jù)穿透后的非均質(zhì)最大面積推掃效率的方程如下確定x：

然后將x代回到原始方程中，能夠計算出穿透后的非均質(zhì)最大面積推掃效率。

縱向推掃效率

下一組計算涉及縱向推掃效率的計算?？v向推掃效率是用相同的方法計算的：

然后，計算體積推掃效率為ev＝ezxea。

穿透后的排送效率

能夠用以下方程計算排送效率，其中sd是被排送的流體的飽和率，swc表示水飽和率：

能夠用累積注射體積估算累積注射孔體積(pv)，其中bw是鹽水的體積因子，winj是累積注射體積，φ是孔隙率，pv是注射孔體積的倍數(shù)：

依照驅(qū)油前端理論，pv等于水的分流量對水的飽和率的微分的倒數(shù)：

因為油的相對滲透率與水的相對滲透率之比是水飽和率的函數(shù)，因此其能夠?qū)懗梢韵滦问?a、b能夠用線性回歸估算)：

水的分流量能夠用以下方程表示，其中μo和μw是油和鹽水的粘度：

在隨后的下一步中，水的分流量與水的飽和率的微分能夠?qū)懗梢韵滦问剑?/p>

能夠根據(jù)該方程估算特定pv的流度比。能夠用剛計算的流度比估算特定時刻水的分流量。根據(jù)水的擴散方程，能夠用水的分流量計算平均水飽和率：

最后，平均排送效率能夠按如下方式計算：

穿透后的采油系數(shù)：

根據(jù)定義，油藏采油系數(shù)等于體積推掃效率和排送效率相乘的乘積：

rf＝ed×ev(46)

一旦水實現(xiàn)了穿透，根據(jù)流動流體理論，其在進入新區(qū)域之前將會流入被氣體推掃的區(qū)域中。這種機理的作用如同滲透率變化的作用，這意味著在被氣體推掃的區(qū)域中形成優(yōu)勢通道。為了正確表達這一機理，引入了滲透率的阻力因子kx。

kv^*＝kx×v

然后，使用該kv*代替原始的滲透率變化v來計算水滲透后的采油系數(shù)。

圖6-17顯示了依照本發(fā)明的方法用不同的油藏和驅(qū)油參數(shù)產(chǎn)生的驅(qū)油操作的樣品繪圖。這些圖在一般意義上用于展示本發(fā)明的方法的實用性和性能，但不應(yīng)當被認為限制了此處所述的本發(fā)明的范圍，其意在覆蓋權(quán)利要求的整個范圍。

在圖6到圖17中的各圖中，標記為雙前曲線的曲線是通過使用此處所述的泡沫驅(qū)動油藏的采油系數(shù)估算方法形成的，而單前曲線是對該模型的簡化，僅考慮了一種氣體前端穿透。用模擬標記的曲線是通過使用計算機數(shù)值模擬軟件(cmg)給出的。

計算機軟件：

圖18和19是顯示在本發(fā)明的方法的計算機軟件迭代中執(zhí)行的步驟的一種實施方案的流程圖。再次，計算機軟件設(shè)計領(lǐng)域的技術(shù)人員將認識到，能夠使用很多不同類型的軟件，包括為此目的特定建立的軟件，或者使用能夠定制以進行該方法的執(zhí)行所需的計算的現(xiàn)有數(shù)學建模軟件。圖18和19的方法將被理解為總體描述了由計算機執(zhí)行的計算機可讀存儲介質(zhì)內(nèi)的處理器指令將進行的方法，如本文總體描述的，在不脫離其預(yù)期范圍和覆蓋范圍的情況下能夠向該總體方法中引入其他公式和方法。

在圖18的方法開始時，顯示為在步驟18-1處捕捉或數(shù)據(jù)輸入油藏和驅(qū)油參數(shù)。油藏和驅(qū)油參數(shù)的捕獲通常包括捕獲例如流度比和滲透率變化的參數(shù)。這些能夠從傳感器總線、人數(shù)據(jù)輸入終端或顯示器、或磁盤或其他內(nèi)存或存儲介質(zhì)上存儲的數(shù)據(jù)捕獲到。

一旦捕獲到該油藏和驅(qū)油參數(shù)，步驟18-2顯示計算該油藏推掃區(qū)內(nèi)的縱向推掃效率，圖18的流程圖顯示上面所述的方程14是將會用于這些計算中的樣品方程，圖18的流程圖中的各個步驟表示了來自說明書其余部分的數(shù)學模型或方程，其將用于特定步驟的執(zhí)行。再次，毫無疑問，即使用于實現(xiàn)該商業(yè)方法步驟的特定公式發(fā)生改變，使用相同的總體方法仍將會實現(xiàn)，且被認為包含在本發(fā)明的范圍之內(nèi)。

在18-2處計算了縱向推掃效率之后，開始驅(qū)油循環(huán)回路，即在該圖的步驟18-3到18-17之間延續(xù)的驅(qū)油循環(huán)回路中對所選擇的各驅(qū)油循環(huán)都執(zhí)行一次與油藏推掃區(qū)采油系數(shù)相關(guān)的技術(shù)。

圖18中顯示的驅(qū)油循環(huán)回路內(nèi)的第一步是計算各個面積，包括泡沫推掃面積、氣體推掃面積和水推掃面積，以及泡沫穿透面積、氣體穿透面積和水穿透面積。所述穿透面積也更詳細地描述在本文的其他位置。確定穿透的存在或大小對于在該驅(qū)油循環(huán)回路內(nèi)該方法中進一步向下計算影響很大。

在完成了18-4處顯示的面積計算之后，步驟18-5指示計算泡沫區(qū)采油系數(shù)，其是使用方程8-16計算的。

在計算了泡沫區(qū)采油系數(shù)(其是來自泡沫推掃子區(qū)的采油系數(shù)或采油貢獻)之后，步驟18-6顯示了第一邏輯塊或步驟，其用于確定在本發(fā)明的方法中接下來應(yīng)用的計算。將泡沫穿透面積、水穿透面積和氣體穿透面積的總和與油氣藏內(nèi)油藏推掃區(qū)的總推掃面積進行比較。如果三個穿透面積之和大于總推掃面積，那么計算邏輯將會返回步驟18-10，而如果總推掃面積大于這三個穿透面積之和，那么該計算將會移動到步驟18-7。這有效地構(gòu)成了采油系數(shù)計算方法的選擇。

步驟18-7顯示了使用方程17-21計算水區(qū)采油系數(shù)，然后使用方程22-26計算氣體區(qū)采油系數(shù)。計算了泡沫區(qū)采油系數(shù)、氣體區(qū)采油系數(shù)和水區(qū)采油系數(shù)這些各個系數(shù)之后，然后能夠在步驟18-9處使用方程27計算該油藏采油系數(shù)，并在數(shù)據(jù)存儲步驟18-16處將其存儲。

或者參照總穿透面積超過總推掃面積的情況，移動到18-10的第二邏輯步驟，將會把泡沫穿透面積和水穿透面積的總和與總推掃面積進行比較。如果總推掃面積小于泡沫穿透面積和水穿透面積之和，就到步驟18-11使用方程29-46在水相中計算水區(qū)采油系數(shù)，然后通過添加泡沫區(qū)采油系數(shù)和水區(qū)采油系數(shù)來計算油藏采油系數(shù)。在步驟18-10的邏輯塊的“否”分支上，用方程17-21計算，然后到步驟18-14用方程29-46在氣相中計算的氣體區(qū)采油系數(shù)。最后，在步驟18-15處顯示，將會依照方程27計算油藏采油系數(shù)。

在各驅(qū)油循環(huán)完成了油藏采油系數(shù)的計算并在步驟18-16處將其存儲之后，該驅(qū)油循環(huán)回路在步驟18-17處終止。如果該方法選擇的驅(qū)油循環(huán)的最大循環(huán)數(shù)尚未完成，該回路將返回步驟18-3，并通過增加其數(shù)值移動到下一個驅(qū)油循環(huán)的計算。或者，如果已經(jīng)達到了最大循環(huán)數(shù)，該驅(qū)油循環(huán)回路將會終止，在繪圖或其他分析中能夠使用所存儲的計算出的油藏采油系數(shù)。

能夠執(zhí)行數(shù)學模型(例如此處所述的)的計算機程序的總體開發(fā)是軟件設(shè)計領(lǐng)域技術(shù)人員所能夠理解的，且所有用于開發(fā)執(zhí)行此處所述方法的根據(jù)目的建立或訂制的建模軟件的方法全部都將包含在本發(fā)明的范圍之內(nèi)。

總體結(jié)果：

表1和表2顯示了對于不同油藏和驅(qū)油參數(shù)(即不同的滲透率變化和注射速率)，依照常規(guī)cmg方法和本發(fā)明的方法計算出的兩組計算結(jié)果。表1顯示了單前端情況的結(jié) 果，在氣體穿透之后沒有突發(fā)性生產(chǎn)。表2顯示了雙前端情況的結(jié)果，其中考慮了在穿透之后不同流度比的影響。

表中列出了在各種情況下生產(chǎn)5年和生產(chǎn)10年之后的采油系數(shù)。將在先方法cmg的結(jié)果與本發(fā)明的結(jié)果進行比較，在表中對于各種情況還計算了百分比的差異。

依照這些結(jié)果，兩種方法都能用于估算采油系數(shù)，與cmg相比本發(fā)明能提供精確的結(jié)果。根據(jù)與模擬的比較，雙前端方法在較低非均質(zhì)度的油氣藏的情況下或者在油氣藏采油的早期階段中提供了更好的結(jié)果，單前端方法在較高非均質(zhì)度的油氣藏的情況下或者在油氣藏采油的后期階段中提供了更好的結(jié)果。

v是滲透率變化，其用于描述油藏的非均質(zhì)度。ugas是氣體注射速率(sc)，uwater是表面活性劑溶液注射速率(sc)。

表1：單前端方法與模擬之間的結(jié)果對比

表2：雙前端方法與模擬之間的結(jié)果對比

應(yīng)該認識到，此處所述的方法的特定步驟和系統(tǒng)的計算機軟件并不被認為限制了本發(fā)明的范圍。本領(lǐng)域技術(shù)人員將容易認識到，在不需要任何創(chuàng)造性勞動的情況下就能夠選擇滿足本發(fā)明的目的的材料和組分。

對本領(lǐng)域技術(shù)人員而言顯而易見的是，在不脫離此處的發(fā)明概念的情況下可以在已經(jīng)描述的那些之外進行很多其他改進。因此，本發(fā)明的主題僅受到后附權(quán)利要求范圍的限制。而且，在解釋說明書和權(quán)利要求書時，所有術(shù)語都應(yīng)當以與上下文相一致的最寬的可能方式加以解釋。特別地，術(shù)語“包括”應(yīng)當被解釋為表示非窮舉方式的要素、組分或步驟，意指所述要素、組分或步驟可以與未明確述及的其他要素、組分或步驟一起存在或使用或結(jié)合。