基于面通量的儲層時變數(shù)值模擬研究

姜瑞忠，喬 欣，滕文超，徐攀登，劉子祎

(中國石油大學(xué)，山東 青島 266580)

?

基于面通量的儲層時變數(shù)值模擬研究

姜瑞忠，喬 欣，滕文超，徐攀登，劉子祎

(中國石油大學(xué)，山東 青島 266580)

針對油田注水開發(fā)過程中，長期水驅(qū)沖刷使儲層物性發(fā)生改變，影響油田開發(fā)效果和剩余油分布的問題，開展室內(nèi)實驗研究、傳統(tǒng)黑油模型改造、軟件編制等工作，建立了物性參數(shù)隨面通量連續(xù)性變化的定量表征方法，形成了能夠描述儲層物性時變的數(shù)值模擬技術(shù)。通過建立概念模型對新模擬器進(jìn)行應(yīng)用，結(jié)果表明：考慮儲層物性變化后，主流線區(qū)水淹程度增大，油藏采出程度提高6.4%；模擬器計算結(jié)果基本不受網(wǎng)格尺寸影響，2種網(wǎng)格模型采出程度和含水率的計算誤差僅為0.079%和0.157%。該技術(shù)解決了現(xiàn)有數(shù)值模擬方法在連續(xù)性表征、方向性表征、計算結(jié)果穩(wěn)定性方面存在的問題，對水驅(qū)油藏開發(fā)指標(biāo)預(yù)測和剩余油分布研究具有重要的應(yīng)用價值。

注水開發(fā)；物性時變；數(shù)值模擬；面通量；剩余油

0 引 言

中國大部分油田采用注水開發(fā)，注入水長期沖刷對油藏產(chǎn)生改造作用，使儲層物性發(fā)生改變，這些改變持續(xù)緩慢進(jìn)行，到中高含水期，物性變化非常明顯[1-2]。儲層物性變化會影響油水運動規(guī)律，對油田開發(fā)效果和剩余油分布產(chǎn)生影響[3-4]，而目前常用的商業(yè)化數(shù)值模擬軟件均未能考慮儲層物性時變現(xiàn)象。據(jù)文獻(xiàn)調(diào)研，目前考慮儲層物性時變的數(shù)值模擬方法可分為4類[5-11]，但在連續(xù)性表征、方向性表征以及計算結(jié)果穩(wěn)定性方面存在某些缺陷，為解決這一問題，需要建立基于儲層參數(shù)隨面通量變化的數(shù)值模擬技術(shù)。

1 儲層物性時變表征方法

由于數(shù)學(xué)模型以及數(shù)值模擬不直接考慮微觀參數(shù)，因此，常用滲透率、孔隙度、相滲曲線等宏觀參數(shù)的變化反映儲層物性時變現(xiàn)象。杜玉洪[12]、王昕立[13]等通過測井解釋以及室內(nèi)巖心實驗等手段，發(fā)現(xiàn)水驅(qū)沖刷后孔隙度的變化幅度較小，而滲透率和相滲曲線的變化比較明顯。此外，考慮到現(xiàn)有數(shù)值模擬方法存在的缺陷，文中提出使用滲透率和相滲曲線隨面通量的變化進(jìn)行定量表征儲層時變現(xiàn)象。

1.1 面通量

面通量定義為累計通過單位面積的水相體積，即：

(1)

式中：M為面通量，m3/m2；Q為通過巖心的總水量體積，m3；S為巖心的橫截面積，m2。

對于任意一個三維空間網(wǎng)格，液體存在x、y、z 3個方向的流動，各個面均存在流體的流入或流出。通過網(wǎng)格的總面通量等于x、y、z 3個方向流出水的面通量之和。

x，y，z方向面通量以及總面通量分別為：

(2)

(3)

式中：Qx、Qy、Qz為3個方向上累計流出的水量體積，m3；Dx、Dy、Dz為3個方向上的網(wǎng)格步長，m。

1.2 滲透率變化規(guī)律

利用勝利油田某中高滲儲層5塊巖心樣品進(jìn)行注水沖刷物理模擬實驗，研究滲透率隨面通量的變化規(guī)律(圖1)。由圖1可知，隨著面通量增大，巖心滲透率逐漸增大，但增大速率呈減小趨勢。面通量小于5m3/m2時，滲透率變化較快，面通量大于5m3/m2時，滲透率趨于平穩(wěn)，變化緩慢。

圖1 滲透率隨面通量的變化關(guān)系

1.3 相對滲透率曲線變化規(guī)律

注入水沖刷會導(dǎo)致油、水相對滲透率變化。圖2為3號巖心樣品的相滲曲線隨面通量的變化關(guān)系(圖中Kro為油相相對滲透率，Krw為水相相對滲透率，M單位為m3/m2)。由圖2可知，隨著面通量增大，束縛水飽和度增大，殘余油飽和度減小，殘余油飽和度對應(yīng)的水相滲透率減小。

2 儲層物性時變數(shù)值模擬

2.1 數(shù)學(xué)模型

由于傳統(tǒng)黑油模型無法考慮儲層物性時變現(xiàn)象，需要對其進(jìn)行改造，建立新的數(shù)學(xué)模型，使之能夠反映滲透率及相對滲透率曲線隨面通量的變化。新模型的輔助方程、初始條件、內(nèi)邊界條件及外邊界條件與傳統(tǒng)黑油模型相同，而油、水、氣三相的連續(xù)性方程發(fā)生變化，見式(4)～(6)。

(4)

(5)

(6)

式中：K為絕對滲透率，μm2；Kro、Krg、Krw分別為油、氣、水的相對滲透率；Bo、Bg、Bw分別為油、氣、水的體積系數(shù)；μo、μg、μw分別為油、氣、水黏度，mPa·s；po、pg、pw分別為油、氣、水的壓力，MPa；ρo、ρg、ρw分別為油、氣、水的密度，kg/m3；So、Sg、Sw分別為油、氣、水的飽和度；qvo、qvg、qvw分別為標(biāo)準(zhǔn)狀況下單位時間產(chǎn)出或注入的油、氣、水的體積，m3/s；Rso、Rsw分別為溶解氣油比、溶解氣水比；g為重力加速度，m/s2；D為從某一基準(zhǔn)面算起的深度，向下為正，m；t為時間，s；φ為孔隙度，%；

傳統(tǒng)黑油模型中絕對滲透率為常數(shù)，相對滲透率由原始的相對滲透率曲線決定，式(4)～(6)中將絕對滲透率和相對滲透率(氣相除外)改造為與面通量有關(guān)的函數(shù)。

2.2 模型求解及模擬器開發(fā)

以儲層物性參數(shù)時變數(shù)學(xué)模型為基礎(chǔ)，開發(fā)了具有自主知識產(chǎn)權(quán)的數(shù)值模擬器。該模擬器利用有限差分方法對微分方程進(jìn)行離散，采用全隱式算法求解壓力和飽和度，具有普通黑油模型的全部功能，且能實現(xiàn)儲層物性時變規(guī)律的有效描述。

與傳統(tǒng)黑油模型相比，該模擬器新增模塊每一時間步的計算過程為：

(1) 首先計算每個網(wǎng)格的壓力和飽和度數(shù)據(jù)場，然后計算每個網(wǎng)格各方向的水相流量，最后計算方向面通量和總面通量。

(2) 根據(jù)滲透率隨方向面通量的變化規(guī)律，計算更新后的滲透率和傳導(dǎo)率數(shù)據(jù)場。

(3) 根據(jù)相對滲透率曲線端點值隨總面通量的變化規(guī)律，計算更新后的相對滲透率曲線。

3 模擬器驗證

為驗證模擬器的準(zhǔn)確性，建立概念模型，x、y方向25個網(wǎng)格，步長為10m，z方向為5個網(wǎng)格，步長為2m，孔隙度為0.25，平面方向滲透率為100×10-3μm2，垂直方向滲透率為10×10-3μm2，初始含油飽和度為0.7，水的密度為1.000g/cm3，油的密度為0.861g/cm3，水的壓縮系數(shù)為4.2×10-4MPa-1，巖石的壓縮系數(shù)為1.5×10-4MPa-1，注采單元為一注一采(I1為注入井，P1為生產(chǎn)井)，定壓生產(chǎn)10a。

在不考慮儲層物性時變的情況下，使用文中軟件和ECLIPSE軟件分別對該模型進(jìn)行模擬計算(圖3)。文中軟件與ECLIPSE軟件的模擬結(jié)果基本一致，采出程度和含水率的誤差僅為0.703%、0.058%，進(jìn)而證明文中軟件的準(zhǔn)確性很高。

圖3 文中軟件與ECLIPSE軟件模擬結(jié)果對比

4 模擬器應(yīng)用

4.1 物性時變影響分析

為研究儲層物性時變對開發(fā)效果以及剩余油分布的影響，將物性參數(shù)隨面通量的變化數(shù)據(jù)加入文中軟件，對上述模型進(jìn)行模擬。定壓生產(chǎn)10a后，考慮時變和不考慮時變的油、水分布模擬結(jié)果見圖4，考慮物性參數(shù)時變后，靠近生產(chǎn)井的主流線區(qū)域水淹程度變大，這是因為注入水的長期沖刷使主流線區(qū)物性變好，流體滲流阻力減小，流量增大。此外，圖5中的計算結(jié)果顯示，考慮時變的采出程度比不考慮時變的采出程度高6.4%，一方面因為儲層物性變好，流體流量增大，另一方面是因為儲層中殘余油飽和度降低使驅(qū)油效率提高。

圖4 不考慮時變和考慮時變剩余油分布對比

圖5 考慮時變與不考慮時變開發(fā)指標(biāo)對比

4.2 計算穩(wěn)定性分析

采用注水沖刷倍數(shù)表征物性參數(shù)的變化存在函數(shù)關(guān)系不確定的缺陷，在數(shù)模中表現(xiàn)為不同網(wǎng)格大小對計算結(jié)果的影響很大，而文中所提方法可解決該問題。例如，某巖心橫截面積為S，孔隙體積為V，水的流量為Q，若將其看作一個網(wǎng)格，則沖刷倍數(shù)為Q/V，面通量為Q/S；若將該巖心劃分為n個網(wǎng)格，對于每一網(wǎng)格，孔隙體積為V/n，水的流量仍為Q，此時，沖刷倍數(shù)為nQ/V，面通量仍為Q/S，由此可見，采用面通量表征物性參數(shù)的變化比采用注水沖刷倍數(shù)表征物性參數(shù)的變化具有穩(wěn)定性好的優(yōu)勢。

為驗證考慮儲層時變情況下網(wǎng)格尺寸對文中軟件計算結(jié)果穩(wěn)定性的影響，建立另一概念模型，x、y方向各20個網(wǎng)格，平面網(wǎng)格步長12.5m，其他參數(shù)不變，與前文x、y方向各25個網(wǎng)格的模型進(jìn)行對比。2種模型的計算結(jié)果見圖6，采出程度和含水率的誤差分別為0.079%和0.157%，可見網(wǎng)格尺寸對計算結(jié)果的影響很小，基本可以忽略，說明基于面通量的儲層時變數(shù)值模擬技術(shù)計算結(jié)果非常穩(wěn)定。

圖6 不同網(wǎng)格大小開發(fā)指標(biāo)對比

5 結(jié) 論

(1) 注入水的長期沖刷會使儲層物性發(fā)生改變，對油田開發(fā)效果和剩余油分布產(chǎn)生影響，分階段地質(zhì)建模和數(shù)值模擬以及基于儲層參數(shù)隨含水飽和度、含水率、注水沖刷倍數(shù)變化的數(shù)值模擬技術(shù)，都具有一定缺陷，無法對儲層物性時變現(xiàn)象進(jìn)行準(zhǔn)確模擬。

(2) 建立了物性參數(shù)隨面通量連續(xù)性變化的定量表征方法，并開發(fā)了基于該方法的數(shù)值模擬技術(shù)，實現(xiàn)了連續(xù)性表征、方向性表征，且計算結(jié)果穩(wěn)定，解決了現(xiàn)有儲層時變數(shù)值模擬技術(shù)存在的問題。

[1] 徐春梅，張榮，馬麗萍，等.注水開發(fā)儲層的動態(tài)變化特征及影響因素分析[J].巖性油氣藏，2010，22(增刊1)：89-92.

[2] 李浩，王香文，劉雙蓮.老油田儲層物性參數(shù)變化規(guī)律研究[J].西南石油大學(xué)學(xué)報(自然科學(xué)版)，2009，31(2)：85-89，184-185.

[3] 陳小凡，劉峰，王博，等.長期沖刷條件下的滲透率變化和微觀模型研究[J].特種油氣藏，2012，19(1)：66-69，138.

[4] 尚朝輝.大孔道非均質(zhì)油藏調(diào)堵封竄參數(shù)優(yōu)化模擬研究[J].油氣地質(zhì)與采收率，2014，21(3)：62-65.

[5] 高博禹，彭仕宓，黃述旺.勝坨油田二區(qū)沙二段3砂層組分階段油藏數(shù)值模擬[J].石油勘探與開發(fā)，2004，31(6)：82-84.

[6] 黃有泉，周志軍，劉志軍，等.葡南三斷塊特高含水期油藏數(shù)值模擬精度提高方法[J].油氣地質(zhì)與采收率，2014，21(5)：65-68.

[7] 姜瑞忠，陳月明，鄧玉珍，等.勝坨油田二區(qū)油藏物理特征參數(shù)變化數(shù)值模擬研究[J].油氣采收率技術(shù)，1996，3(2)：50-56.

[8] 吳素英.長期注水沖刷儲層參數(shù)變化規(guī)律及對開發(fā)效果的影響[J].大慶石油地質(zhì)與開發(fā)，2006，25(4)：35-37.

[9] 崔傳智，耿正玲，王延忠，等.水驅(qū)油藏高含水期滲透率的動態(tài)分布計算模型及應(yīng)用[J].中國石油大學(xué)學(xué)報(自然科學(xué)版)，2012，54(4)：118-122.

[10] 劉顯太.中高滲透砂巖油藏儲層物性時變數(shù)值模擬技術(shù)[J].油氣地質(zhì)與采收率， 2011，18(5)：58-62.

[11] 張喬良，姜瑞忠，姜平，等.油藏流場評價體系的建立及應(yīng)用[J].大慶石油地質(zhì)與開發(fā)，2014，33(3)：86-89.

[12] 杜玉洪，張繼春，侯翠芬.儲層宏觀物性參數(shù)隨注水開發(fā)動態(tài)演化模式研究[J].特種油氣藏，2004，11(5)：52-55，122-123.

[13] 王昕立，姚遠(yuǎn)，劉峰，等.長期沖刷條件下的儲層物性參數(shù)變化規(guī)律研究[J].重慶科技學(xué)院學(xué)報(自然科學(xué)版)，2011，6(3)：8-9，20.

編輯 張耀星

20151116；改回日期：20151225

國家科技重大專項“大型油氣田及煤層氣開發(fā)”子課題“特低滲油藏有效開發(fā)技術(shù)”(2011ZX0513-006)；國家自然基金“基于新模型的低滲透油藏非線性滲流理論”(51174223)

姜瑞忠(1964-)，男，教授，博士生導(dǎo)師，1987年畢業(yè)于西南石油學(xué)院油藏工程專業(yè)，2002年畢業(yè)于西南石油大學(xué)油氣田開發(fā)工程專業(yè)，獲博士學(xué)位，現(xiàn)主要從事油氣田開發(fā)方面的教學(xué)和科研工作。

10.3969/j.issn.1006-6535.2016.02.016

TE

A