張振濤，姜漢橋，李俊鍵，趙　林，裴艷麗，馬　康

(中國石油大學(xué)(北京) 石油工程學(xué)院，北京昌平 102249)

?

動態(tài)毛管力對低滲油藏開發(fā)動態(tài)的影響

張振濤，姜漢橋，李俊鍵，趙林，裴艷麗，馬康

(中國石油大學(xué)(北京) 石油工程學(xué)院，北京昌平 102249)

油水兩相滲流過程中，毛管力不僅是飽和度的函數(shù)，還是飽和度變化率的函數(shù)；動態(tài)毛管力與黏度、滲透率、驅(qū)替速率等因素相關(guān)。通過建立考慮動態(tài)毛管力效應(yīng)的油水兩相滲流數(shù)學(xué)模型，利用數(shù)值模擬軟件研究動態(tài)毛管力對油藏開發(fā)動態(tài)的影響，并考察動態(tài)毛管力與儲層非均質(zhì)性的關(guān)系。結(jié)果表明，動態(tài)毛管力造成水驅(qū)油過程中含水率加速上升，產(chǎn)油量下降，產(chǎn)生阻力效應(yīng)，同時加劇了油藏的非均質(zhì)性。

動態(tài)毛管力；低滲油藏；開發(fā)動態(tài)；數(shù)值模擬；非均質(zhì)性

1　研究背景

毛管力在油水兩相滲流過程中起到十分重要的作用，為了充分理解流體在地層中的流動狀態(tài)與流動特點，應(yīng)對毛管力與流體飽和度之間的關(guān)系進行準(zhǔn)確表征。目前普遍認(rèn)為毛管力是潤濕相流體飽和度的函數(shù)，等于潤濕相流體與非潤濕相流體壓力差[1]。 目前，主要通過室內(nèi)實驗確定毛管力曲線，如半滲隔板法[2]；Topp等人[3]發(fā)現(xiàn)，隨著驅(qū)替速率的增加，毛管力數(shù)值隨之變大，因此在平衡態(tài)得到毛管力曲線不能完全適用于動態(tài)過程；Hassanizadeh和Gray從熱力學(xué)角度得到毛管力本構(gòu)方程[4-5]；Stauffer于1978年提出了動態(tài)毛管力系數(shù)的經(jīng)驗公式，認(rèn)為動態(tài)毛管力系數(shù)與多孔介質(zhì)性質(zhì)及流體性質(zhì)相關(guān)[6]。

目前，動態(tài)毛管力數(shù)值模擬主要在孔隙或巖心尺度進行模擬[7-10]，或者在建模過程中沒有考慮動態(tài)毛管力的影響[11-12]。針對以上問題，本文在油藏尺度進行建模，并探究動態(tài)毛管力對油藏開發(fā)動態(tài)的影響。

2　考慮動態(tài)毛管力油水兩相滲流模型

建立了油水兩相三維滲流數(shù)學(xué)模型，并使用MRST (MATLAB Reservoir Simulation Toolbox)[13]進行隱式求解。

2.1模型假設(shè)

(1)滲流過程在等溫條件下進行。

(2)多孔介質(zhì)中僅存在油水兩相，并且滲流過程符合達西定律。

(3)油藏中油相為潤濕相，水相為非潤濕相。

(4)考慮重力的影響。

2.2數(shù)學(xué)模型的建立

油水兩相滲流連續(xù)性方程如式(1)所示：

(1)

式(1)中兩個方程4個未知量，求解時需輔助方程：

Sw+So=1

(2)

考慮動態(tài)毛細管力：

(3)

(4)

初始條件： p(x,y,z,0)=pini=234atm；

井工作條件：pini=0.5 stb/d，ppro=0.5 stb/d。

計算所用的油藏基本參數(shù)如表1所示，根據(jù)油藏基本參數(shù)，采用一注一采工作制度。

表1　油藏基本參數(shù)

3　動態(tài)毛管力對油藏開發(fā)動態(tài)的影響

根據(jù)建立的數(shù)學(xué)模型，研究動態(tài)毛管力對油藏壓力、產(chǎn)量和含水率等開發(fā)動態(tài)的影響，探究油藏非均質(zhì)性對動態(tài)毛管力的影響。

3.1毛管力動態(tài)效應(yīng)

分別建立考慮動態(tài)毛管力的模型與不考慮動態(tài)毛管力(τ=0)的模型，分析動態(tài)毛管力對油藏開發(fā)動態(tài)的影響。

含水率變化如圖1所示，可以看出，考慮動態(tài)毛管力之后，含水率上升變快。由圖2可以看出，考慮動態(tài)毛管力之后，日產(chǎn)油量降低，累計產(chǎn)油量少于不考慮動態(tài)毛管力的情況。地層壓力場如圖3所示，可見考慮動態(tài)毛管力之后地層壓力較大。這說明在親油油藏中，水驅(qū)油過程毛管力為阻力，考慮動態(tài)毛管力之后，增大了這種阻力效應(yīng)，滲流阻力變大，含水率上升變快，從而造成了產(chǎn)量的下降。

圖1　含水率曲線

圖2　產(chǎn)量曲線

3.2非均質(zhì)性影響

為了分析地層非均質(zhì)性對動態(tài)毛管力的影響，建立兩個非均質(zhì)性不同的模型，如圖4所示，非均質(zhì)模型的主要參數(shù)如表2所示，可以看出模型2非均質(zhì)性更強。

圖4　非均質(zhì)模型示意圖

模擬結(jié)果如圖5、圖6所示，非均質(zhì)模型2的含水率上升速度高于非均質(zhì)模型1?？紤]動態(tài)毛管力之后，非均質(zhì)模型2的累計產(chǎn)油量下降率顯著高于非均質(zhì)模型1。這說明在非均質(zhì)性較強的模型2中，動態(tài)毛管力起到了相對更強的阻力效應(yīng)。

表2　非均質(zhì)模型參數(shù)

圖5　非均質(zhì)模型含水率曲線

圖6　非均質(zhì)模型累計產(chǎn)油量對比

由于儲層非均質(zhì)越強，孔隙結(jié)構(gòu)越復(fù)雜，驅(qū)替前緣越難以達到平衡狀態(tài)，動態(tài)毛管力系數(shù)越大，造成的阻力效應(yīng)越大，相當(dāng)于動態(tài)毛管力加劇了地層的非均質(zhì)性。

4　結(jié)論

(1)通過建立數(shù)值模擬理論模型，研究了動態(tài)毛管力對油藏開發(fā)動態(tài)的影響，結(jié)果表明動態(tài)毛管力造成親油油藏含水率上升變快、產(chǎn)油量下降，起到阻力作用。

(3)建立了兩個非均質(zhì)性不同的理論模型，探究非均質(zhì)性與動態(tài)毛管力的關(guān)系，結(jié)果表明非均質(zhì)性較嚴(yán)重的模型，在考慮動態(tài)毛管力之后產(chǎn)量下降率顯著變大，說明動態(tài)毛管力起到了加劇非均質(zhì)性的作用。

1673-8217(2016)04-0115-04

2016-02-26

張振濤，在讀碩士研究生，1991年生，現(xiàn)從事油氣田開發(fā)理論與提高采收率方面的研究。

國家重點基礎(chǔ)研究發(fā)展計劃(973)項目子課題“致密油高效開發(fā)油藏工程理論與方法研究”(2015250905)。

TE34

A