一維水驅(qū)油變壓力驅(qū)替過(guò)程推導(dǎo)

劉海龍，由春梅，姜雪巖，王曙光

（中國(guó)石油大慶油田有限責(zé)任公司勘探開(kāi)發(fā)研究院，黑龍江大慶 163712）

對(duì)于水驅(qū)油藏，Buckley－Leverett方程給出了累積注入量與含水飽和度分布的關(guān)系，通過(guò)這個(gè)關(guān)系，在某一時(shí)刻，在已知相滲曲線、油層孔滲數(shù)據(jù)、注采井距、總注入量、注采壓差的情況下，能夠求得含水飽和度分布、水驅(qū)前緣位置（采出端未見(jiàn)水時(shí)）、注入速度、采出端含水率等一系列動(dòng)態(tài)指標(biāo)［1－5］。以上求解是針對(duì)某一瞬時(shí)，而研究驅(qū)替過(guò)程就需要研究驅(qū)替狀態(tài)隨時(shí)間的變化情況，就需要引入時(shí)間這個(gè)變量。以往的解決方法都是設(shè)定注入速度不變，這樣就把時(shí)間與總注入量通過(guò)恒定的注入速度直接聯(lián)系起來(lái)，從而把時(shí)間引入到驅(qū)替理論體系中［6］。

如果要設(shè)定壓力為已知條件去求解驅(qū)替過(guò)程，由于注入速度不再是恒定不變，而且注入速度隨壓力變化沒(méi)有簡(jiǎn)單明了的直接關(guān)系［7－12］，就無(wú)法利用累積注入量、注入速度、時(shí)間之間的關(guān)系把時(shí)間引入，針對(duì)這一問(wèn)題，作者在另一篇文章中推導(dǎo)了一維水驅(qū)油恒壓驅(qū)替過(guò)程［13］，給出了在恒定壓力情況下水驅(qū)油過(guò)程中各項(xiàng)動(dòng)態(tài)數(shù)據(jù)隨時(shí)間變化的情況。本文以其中的中間結(jié)果為基礎(chǔ)，進(jìn)一步推導(dǎo)了變壓力情況下的水驅(qū)油驅(qū)替過(guò)程。

1 采出端見(jiàn)水之前的驅(qū)替

假設(shè)條件：一維油層，一注一采；不考慮巖石、流體的壓縮性、毛管力的作用；截面積為1 m2。通常，壓力的變化可以認(rèn)為是隨時(shí)間不斷變化，這種模式是相對(duì)比較多見(jiàn)，因此本文以壓力隨時(shí)間變化為設(shè)定條件，某時(shí)刻t，注入端壓力pw（t），采出端壓力po（t），注采壓差p（t）＝po（t）－pw（t）；開(kāi)始注水為初始0時(shí)刻；注入端為x＝0位置；注采兩端距離為L(zhǎng)，水黏度為μw，油黏度為μo，束縛水飽和度為Swc，殘余油飽和度為Sor；掌握比較規(guī)范的相滲數(shù)據(jù)并根據(jù)相滲曲線計(jì)算的分流率f（Sw）、分流率導(dǎo)數(shù)f′（Sw），驅(qū)替過(guò)程中t時(shí)刻含水飽和度的分布見(jiàn)圖1，相對(duì)滲透率、分流率、分流率導(dǎo)數(shù)與含水飽和度的關(guān)系曲線見(jiàn)圖2。

圖1 采出端見(jiàn)水前含水飽和度分布

某一時(shí)刻t，由文獻(xiàn)［13］中采出端見(jiàn)水前的積分式可知：

式中：K——砂巖滲透率，10－3μm2；φ——孔 隙度，%；系數(shù)α參見(jiàn)文獻(xiàn)［13］；Swf——水驅(qū)前緣含水飽和度，%；Q——累積注入量，m3。

整理（1）式，并積分求得某時(shí)刻tn的總注入量：

圖2 相對(duì)滲透率、分流率、分流率導(dǎo)數(shù)關(guān)系曲線

定義Ar是壓力曲線與表示時(shí)間的x軸在坐標(biāo)上圍出的面積，可以稱(chēng)為壓時(shí)面積，即Ar＝。則（2）式可以整理為：

考慮條件t＝0時(shí)，Q＝0，則式（3）的解為

式（4）給出了壓時(shí)面積、總注入量之間的函數(shù)關(guān)系，進(jìn)一步根據(jù)Buckley－Leverett方程可求得給定時(shí)刻的注入速度、飽和度分布、見(jiàn)水時(shí)間等相關(guān)數(shù)據(jù)。并由此可見(jiàn)采出端見(jiàn)水前的驅(qū)替進(jìn)程是由壓時(shí)面積決定的。

（5）式給出了一維水驅(qū)油采出端見(jiàn)水前壓時(shí)面積、位置、含水飽和度之間的關(guān)系。水驅(qū)前緣至采出端的含水飽和度為束縛水飽和度。

2 采出端見(jiàn)水之后的驅(qū)替

隨著水驅(qū)油過(guò)程的不斷進(jìn)行，水驅(qū)前緣不斷向前推進(jìn)，設(shè)采出端見(jiàn)水時(shí)刻為t0，見(jiàn)水后含水飽和度的分布見(jiàn)圖3。

圖3 見(jiàn)水后含水飽和度的分布

某一時(shí)刻t，由文獻(xiàn)［13］中采出端見(jiàn)水后的微分式可知：

式中：Swe——出口端含水飽和度；φ1函數(shù)參見(jiàn)文獻(xiàn)［13］。

t0時(shí)刻出口端含水飽和度為水驅(qū)前緣含水飽和度Swf，見(jiàn)水后某時(shí)刻tn1出口端的含水飽和度為Swen，整理（6）式并積分：

式中：Sw為對(duì)應(yīng)于某一位置 的含水飽和度。

（8）式即為一維水驅(qū)油采出端見(jiàn)水后壓時(shí)面積、位置、含水飽和度之間關(guān)系的計(jì)算公式。可見(jiàn)采出端見(jiàn)水后的驅(qū)替進(jìn)程是由壓時(shí)面積決定的。

3 舉例計(jì)算

選取大慶油田滲透率為0.1μm2標(biāo)準(zhǔn)的相對(duì)滲透率曲線［14］，注入端、采出端壓力差與時(shí)間的關(guān)系：P（t）＝6＋0.001t（壓力單位 MPa，時(shí)間單位千秒），油層滲透率0.1μm2，兩端距離250 m，截面積1 m2，孔隙度0.267。以上述條件計(jì)算驅(qū)替過(guò)程，求解日注入量并與數(shù)值模擬結(jié)果對(duì)比，給出第100天和第400天的含水飽和度分布。

3.1 日注入量求解過(guò)程

（1）由 Buckley－Leverett方 程 L ＝f′（Swf）Qw／φ，計(jì)算見(jiàn)水時(shí)刻的總注入量Qw。

（2）計(jì)算每天的壓時(shí)面積，并由（4）式計(jì)算采出端見(jiàn)水前累積注入量，計(jì)算結(jié)果選取小于Qw的數(shù)值。

（3）當(dāng)由（4）式計(jì)算的累積注入量第一個(gè)大于Qw的時(shí)間點(diǎn)，使用見(jiàn)水后的求解方法：計(jì)算見(jiàn)水后每天的壓時(shí)面積，并由（7）式計(jì)算每天的φ2函數(shù)值。

（4）使用作圖法，根據(jù)φ2函數(shù)值求得f′（Swen），并由（8）式計(jì)算見(jiàn)水后每天的累積注入量。（φ2函數(shù)的內(nèi)容參見(jiàn)文獻(xiàn)［13］）

（5）根據(jù)見(jiàn)水前、見(jiàn)水后相鄰兩天的累積注入量的差值，得到每天的日注入量。

（6）使用數(shù)值模擬方法，建立一維模型，按照給出的驅(qū)替條件模擬計(jì)算，并對(duì)比模擬結(jié)果和理論計(jì)算結(jié)果，二者的對(duì)比曲線如圖4所示。

3.2 含水飽和度分布求解過(guò)程

（1）（6）式Lf取注采井距250 m，求得見(jiàn)水時(shí)刻壓時(shí)面積，根據(jù)壓力與時(shí)間給定的關(guān)系，求得見(jiàn)水時(shí)間為156天，可判斷第100天未見(jiàn)水，并計(jì)算第100天的壓時(shí)面積。

（2）100天的壓時(shí)面積代入（6）式計(jì)算第100天水驅(qū)前緣位置為122.95 m

（3）第100天的壓時(shí)面積代入（5）式計(jì)算注入端至水驅(qū)前緣的含水飽和度分布，即每一個(gè)位置x確定對(duì)應(yīng)的含水飽和度Sw，水驅(qū)前緣至采出端為束縛水飽和度。

（4）第400天大于見(jiàn)水時(shí)間156天，采出端見(jiàn)水，計(jì)算第400天的壓時(shí)面積。

（5）第100天的壓時(shí)面積代入（9）式求得φ2函數(shù)，使用作圖法，根據(jù)φ2函數(shù)值求得函數(shù)自變量值，從而得到位置x與含水飽和度Sw的關(guān)系。

計(jì)算結(jié)果如圖5、圖6所示

圖5 第100天含水飽和度分布

圖6 第400天含水飽和度分布

4 結(jié)論與認(rèn)識(shí)

（1）可以通過(guò)壓時(shí)面積把時(shí)間、隨時(shí)間變化的壓力引入到驅(qū)替理論體系中，建立一維水驅(qū)油變壓力驅(qū)替過(guò)程的求解方法，進(jìn)一步完善一維水驅(qū)油理論體系。

（2）壓時(shí)面積決定了一維水驅(qū)油變壓力驅(qū)替的進(jìn)程。

（3）一維水驅(qū)油變壓力驅(qū)替過(guò)程中，驅(qū)替動(dòng)態(tài)與儲(chǔ)層孔隙度、滲透率、注采井距存在確定的數(shù)學(xué)關(guān)系，這一結(jié)論可用于研究油層靜態(tài)參數(shù)對(duì)開(kāi)發(fā)動(dòng)態(tài)的影響。

［1］葛家理，寧正福，劉月田，等.現(xiàn)代油藏滲流力學(xué)原理［M］.北京：石油工業(yè)出版社，2001：127－135.

［2］李傳亮.油藏工程原理［M］.北京：石油工業(yè)出版社，2005：297－307.

［3］宋考平，王立軍，何鮮，等.單層剩余油分布及動(dòng)態(tài)指標(biāo)預(yù)測(cè)動(dòng)態(tài)劈分法［J］.石油學(xué)報(bào)，2000，21（6）：122－126.

［4］宋考平，吳玉樹(shù)，計(jì)秉玉.水驅(qū)油藏剩余油飽和度分布預(yù)測(cè)的 函數(shù)法［J］.石油學(xué)報(bào)，2006，27（3）：91－95.

［5］張順康，劉炳官，鐘思瑛，等.多層油藏水驅(qū)開(kāi)發(fā)理論模型［J］.新疆石油地質(zhì)，2009，30（6）：734－737.

［6］［美］G.鮑爾.威爾海特.注水［M］.北京：石油工業(yè)出版社，2001：127－135.

［7］王家祿，沈平平，劉玉章，等.油水兩相滲流壓力脈動(dòng)的實(shí)驗(yàn)及數(shù)值分析［J］.實(shí)驗(yàn)力學(xué)，2001，16（3）：331－337.

［8］宋考平，高群峰，甘曉飛，等.喇嘛甸油田特高含水期油井合理地層壓力及流壓界限［J］.大慶石油學(xué)院學(xué)報(bào)，2004，28（2）：98－100.

［9］康凱，馮敏，李彥來(lái)，等.旅大5－2油田合理生產(chǎn)壓差確定［J］.新疆石油地質(zhì)，2011，32（2）：170－172.

［10］張繼成，梁文福，趙玲，等.喇嘛甸油田特高含水期開(kāi)發(fā)形勢(shì)分析［J］.大慶石油學(xué)院學(xué)報(bào)，2005，29（3）：23－25.

［11］張東，孫淑霞，李音，等.對(duì)杏樹(shù)崗油田北部主力油層流動(dòng)單元的劃分［J］.大慶石油學(xué)院學(xué)報(bào)，2002，26（3）：18－21.

［12］謝曉慶，姜漢橋，劉同敬，等.非均質(zhì)油藏分層配水理論［J］.大慶石油地質(zhì)與開(kāi)發(fā)，2008，27（6）：83－86.

［13］劉海龍.一維水驅(qū)油恒壓驅(qū)替滲流過(guò)程推導(dǎo)［J］.大慶石油學(xué)院學(xué)報(bào)，2012，36（3）：90－95.

［14］王曙光，趙國(guó)忠，余碧君.大慶油田油水相對(duì)滲透率統(tǒng)計(jì)規(guī)律及其應(yīng)用［J］.石油學(xué)報(bào)，2005，26（3）：78－85.