于春磊，糜利棟，王川，趙玉云，姜漢橋，田野

（1.中國(guó)石化勝利油田分公司勘探開(kāi)發(fā)研究院，山東東營(yíng)257015；2.中國(guó)石油大學(xué)（北京）石油工程教育部重點(diǎn)實(shí)驗(yàn)室，北京102249）

水驅(qū)油藏特高含水期微觀剩余油滲流特征研究

于春磊1，糜利棟2，王川2，趙玉云2，姜漢橋2，田野2

（1.中國(guó)石化勝利油田分公司勘探開(kāi)發(fā)研究院，山東東營(yíng)257015；2.中國(guó)石油大學(xué)（北京）石油工程教育部重點(diǎn)實(shí)驗(yàn)室，北京102249）

水驅(qū)油藏進(jìn)入特高含水期后仍具有一定挖潛潛力，微觀剩余油流動(dòng)特征及運(yùn)移規(guī)律的研究對(duì)于特高含水期水驅(qū)油藏提高采收率具有重要意義。文中通過(guò)微觀玻璃刻蝕模型實(shí)驗(yàn)和計(jì)算機(jī)圖像識(shí)別處理技術(shù)對(duì)微觀剩余油流動(dòng)形態(tài)進(jìn)行分類(lèi)，分析微觀剩余油流動(dòng)形態(tài)及變化規(guī)律。研究結(jié)果表明：特高含水期，考慮油-水-孔喉接觸關(guān)系，可將剩余油流動(dòng)形態(tài)分為簇狀流、多孔流、柱狀流、膜狀流和滴狀流5類(lèi)，其中簇狀流分布比例最多，是引起相滲曲線下彎的主要因素，并且隨著含水飽和度的上升，簇狀流逐漸轉(zhuǎn)化為多孔流、柱狀流、膜狀流和滴狀流；同時(shí)從微觀上解釋了特高含水期油水微觀流動(dòng)規(guī)律及相滲曲線非線性的原因，為特高含水期油藏提高剩余油動(dòng)用程度和采收率提供指導(dǎo)。

水驅(qū)油藏；特高含水；計(jì)算機(jī)圖像識(shí)別技術(shù)；滲流特征

0　引言

水驅(qū)油藏進(jìn)入特高含水期后仍具有一定挖潛潛力，研究微觀剩余油在孔隙尺度下的流動(dòng)特征及運(yùn)移規(guī)律，對(duì)水驅(qū)油藏特高含水期提高采收率具有重要作用［1-2］。但目前關(guān)于剩余油的研究大多局限于靜態(tài)剩余油［3］，而對(duì)產(chǎn)量具有直接貢獻(xiàn)作用的實(shí)際上是流動(dòng)的剩余油。特別是特高含水期及特高含水后期，油水相對(duì)滲透率之比和飽和度半對(duì)數(shù)曲線不再是線性關(guān)系［4］，動(dòng)態(tài)剩余油的研究顯得更加重要。本文通過(guò)微觀玻璃刻蝕模型實(shí)驗(yàn)和計(jì)算機(jī)圖像識(shí)別處理技術(shù)對(duì)微觀剩余油流動(dòng)形態(tài)進(jìn)行分類(lèi)，為了使成果更具有普遍性，研究了不同孔喉比、孔喉半徑、配位數(shù)、流體性質(zhì)，以及注入條件下微觀剩余油流動(dòng)形態(tài)、變化規(guī)律和不同流動(dòng)形態(tài)剩余油微觀分布特征。

1　物理實(shí)驗(yàn)

玻璃刻蝕模型可視化實(shí)驗(yàn)是實(shí)驗(yàn)室研究微觀滲流的重要手段［5］。為了研究特高含水期微觀剩余油流動(dòng)形態(tài)及變化規(guī)律，設(shè)計(jì)了不同孔喉半徑、孔喉比、配位數(shù)形狀因子等特征參數(shù)的概念模型和均質(zhì)、非均質(zhì)實(shí)際模型，完成了不同原油黏度和驅(qū)替條件下的微觀水驅(qū)油實(shí)驗(yàn)。

1.1模型設(shè)計(jì)

6塊概念模型和3塊實(shí)際模型的具體特征參數(shù)如表1所示。

表1　模型孔喉特征參數(shù)

1.2驅(qū)替實(shí)驗(yàn)

為了研究孔隙特征參數(shù)、流體性質(zhì)以及注入條件等因素對(duì)特高含水后期不同剩余油流動(dòng)形態(tài)的影響，使得結(jié)果更具普遍性和實(shí)用性，在不同孔隙特征參數(shù)模型上完成了不同流體黏度和注入速度的正交實(shí)驗(yàn)。

1.2.1實(shí)驗(yàn)材料與設(shè)備

模擬油由不同比例原油與煤油配制而成，黏度分別為2，4，6 mPa·s；實(shí)驗(yàn)用水為模擬地層水；玻璃刻蝕模型尺寸2.0 cm×1.5 cm；實(shí)驗(yàn)設(shè)備為HAS-200型恒壓恒速泵和微觀驅(qū)替裝置。

1.2.2實(shí)驗(yàn)步驟

1）用注射裝置向模型中緩慢注入模擬地層水，使模型孔隙中充分飽和模擬水。

2）用注射裝置向飽和水的模型中緩慢注入一定黏度的模擬油，使模擬油將模型孔隙中的可流動(dòng)水驅(qū)出，同時(shí)使模擬油占據(jù)模型孔隙，此時(shí)模型處于飽和油和束縛水狀態(tài)。

3）打開(kāi)實(shí)時(shí)截屏錄制軟件，設(shè)置恒壓恒速泵為一定驅(qū)替速度，使用微觀驅(qū)替裝置對(duì)模型進(jìn)行水驅(qū)。

4）觀察玻璃刻蝕模型中的油水分布，孔隙中剩余油不再變化時(shí)停止水驅(qū)，結(jié)束視頻錄制。

5）不同模型用不同原油黏度在不同驅(qū)替速度下，重復(fù)步驟1）—4），完成所有實(shí)驗(yàn)，為下一步通過(guò)計(jì)算機(jī)圖像識(shí)別技術(shù)研究剩余油流動(dòng)形態(tài)分類(lèi)及流動(dòng)規(guī)律做準(zhǔn)備。

2　計(jì)算機(jī)圖像處理

通常情況下，實(shí)驗(yàn)人員通過(guò)經(jīng)驗(yàn)對(duì)微觀玻璃刻蝕模型剩余油進(jìn)行手動(dòng)分類(lèi)標(biāo)記［6-9］，而這種傳統(tǒng)的分類(lèi)標(biāo)記方法一方面不利于圖像批量處理，另一方面由于觀察者經(jīng)驗(yàn)及識(shí)別精度的局限性導(dǎo)致結(jié)果差別很大。通過(guò)計(jì)算機(jī)技術(shù)不僅可以大幅度提高處理速度（圖像批量處理）［10］，而且可以得到統(tǒng)一的處理結(jié)果。計(jì)算機(jī)圖像處理技術(shù)應(yīng)用在剩余油的分類(lèi)識(shí)別時(shí)，首先需要解決的2個(gè)基礎(chǔ)問(wèn)題是巖石、油、水精確分離和孔-喉結(jié)構(gòu)的分離。

2.1巖石、油、水精確分離

在對(duì)每張實(shí)驗(yàn)圖片的剩余油進(jìn)行分類(lèi)識(shí)別統(tǒng)計(jì)之前，需要對(duì)圖片進(jìn)行“三值化”處理，也就是進(jìn)行巖石、油、水的精確分離。首先采用“背景窗口移植法”，通過(guò)對(duì)實(shí)驗(yàn)視頻第1幀（即飽和油狀態(tài)）圖片進(jìn)行精確二值化，得到巖石骨架，實(shí)現(xiàn)油、水整體和巖石骨架之間的精確分類(lèi)；然后借助第1步得到的“骨架”作為背景，對(duì)提取的油、水部分進(jìn)行第2次二值化，實(shí)現(xiàn)水驅(qū)油過(guò)程中不同驅(qū)替階段圖像的三值化處理，從而達(dá)到巖石、油、水精確分離的目的。

圖像處理時(shí)，“分割窗口閾值法”［11］可以有效解決因圖像曝光度不均勻引起的識(shí)別精度不高的的問(wèn)題，而“單點(diǎn)窗口閾值法”［12］可以有效降低“分割窗口閾值法”中邊緣處的識(shí)別誤差，筆者將“分割窗口閾值法”和“單點(diǎn)窗口閾值法”相結(jié)合，使得最終識(shí)別精度達(dá)到95%以上。圖1為水驅(qū)某階段的分離結(jié)果，其中黑色是油（閾值為0），白色是水（閾值為1.0），灰色是巖石骨架（閾值為0.5），對(duì)應(yīng)的灰度值分別為0，255，128。

2.2孔-喉結(jié)構(gòu)分離

基于“邊界收縮法”實(shí)現(xiàn)了孔-喉拓?fù)浣Y(jié)構(gòu)的提取，采用“射線法”和“最大類(lèi)間距法”［13］對(duì)1.1中的6個(gè)概念模型進(jìn)行孔-喉分離處理。結(jié)果表明，該方法分離精度可達(dá)到98%以上，圖2為其中一個(gè)實(shí)際模型的分離結(jié)果。

圖1　巖石、油、水分離結(jié)果

圖2　孔隙-喉道分離結(jié)果

3　油水微觀滲流特征分析

3.1剩余油流動(dòng)形態(tài)劃分及定量表征

分析對(duì)比實(shí)驗(yàn)結(jié)果發(fā)現(xiàn)，特高含水后期考慮油-水-孔喉接觸關(guān)系可將剩余油流動(dòng)狀態(tài)分為膜狀流、滴狀流、柱狀流、多孔流以及簇狀流5種，分別定義為：1）膜狀流，有且只有1個(gè)油水界面，油-水界面約占油相表面一半，厚度小于孔喉直徑的1/3；2）滴狀流，有且只有1個(gè)油水界面，油相不與固相接觸，孔隙和喉道數(shù)不大于1；3）柱狀流，剩余油分布在喉道中，有2個(gè)油水界面，孔隙和喉道數(shù)不大于1；4）多孔流，油相分布于多個(gè)孔喉，孔隙和喉道數(shù)不大于5；5）簇狀流，油相分布于多個(gè)孔喉，孔隙和喉道數(shù)大于5）。這5種類(lèi)型可以包含所有剩余油。

為了方便計(jì)算機(jī)程序的處理，分別引入形狀因子G、接觸比C、歐拉數(shù)EN等特征參數(shù)對(duì)5類(lèi)剩余油進(jìn)行定量表征，并通過(guò)大量數(shù)值實(shí)驗(yàn)，得到不同類(lèi)型剩余油對(duì)應(yīng)不同特征參數(shù)的取值范圍（見(jiàn)表2）。圖3為某水驅(qū)實(shí)驗(yàn)驅(qū)替到某時(shí)刻的分類(lèi)結(jié)果，其中簇狀流和多孔流分布最多，柱狀流次之，膜狀流和滴狀流最少。

表2　不同流動(dòng)形態(tài)剩余定量表征

圖3　剩余油分類(lèi)識(shí)別結(jié)果

3.2不同流動(dòng)形態(tài)變化及相互演化規(guī)律

1）隨著水驅(qū)油過(guò)程中含水飽和度的增加，油相不斷被分割，非連續(xù)程度不斷加劇，表現(xiàn)為數(shù)量多、體積小、分散性強(qiáng)的特點(diǎn)；2）非連續(xù)程度加劇是因?yàn)榇貭盍飨蚨嗫琢鞯确沁B續(xù)流動(dòng)形態(tài)轉(zhuǎn)變，原油的流動(dòng)越來(lái)越分散；3）隨著含水飽和度的增加，非連續(xù)流動(dòng)形態(tài)（滴狀流、膜狀流、柱狀流、多孔流）比例持續(xù)增加，連續(xù)流動(dòng)形態(tài)（簇狀流）比例逐漸減小（見(jiàn)圖4）；4）簇狀流作為連續(xù)相，動(dòng)用能力較強(qiáng)，油相相對(duì)滲透率較大，其他4種形態(tài)的剩余油作為非連續(xù)相，動(dòng)用能力較弱，油相相對(duì)滲透率較小，而相滲曲線拐點(diǎn)的出現(xiàn)時(shí)機(jī)恰為相對(duì)滲透率較小的非連續(xù)相開(kāi)始發(fā)育的時(shí)間，也就是說(shuō)抑制連續(xù)相向非連續(xù)相的轉(zhuǎn)化及非連續(xù)相發(fā)育，把連續(xù)相所占比例維持在一定水平上，可以在某種程度上延遲相滲曲線拐點(diǎn)出現(xiàn)時(shí)機(jī)，有助于提高采收率。

圖4　剩余油數(shù)量和平均體積變化曲線

4　結(jié)論

1）特高含水期，考慮油-水-孔喉接觸關(guān)系可將剩余油流動(dòng)形態(tài)分為簇狀流、多孔流、柱狀流、膜狀流和滴狀流，其中簇狀流分布比例最多，并且隨著含水飽和度的上升，連續(xù)相逐漸轉(zhuǎn)化為非連續(xù)相。

2）抑制連續(xù)相向非連續(xù)相的轉(zhuǎn)化及非連續(xù)相發(fā)育，把連續(xù)相所占比例維持在一定水平上，可以延遲相滲曲線拐點(diǎn)出現(xiàn)時(shí)機(jī)，有助于提高采收率。

3）從微觀上解釋了特高含水期油水微觀流動(dòng)規(guī)律及相滲曲線非線性的原因，為特高含水期油藏提高剩余油動(dòng)用程度和提高采收率提供指導(dǎo)。

［1］武麗麗.特高含水期水驅(qū)油藏層系井網(wǎng)重組優(yōu)化研究［D］.青島：中國(guó)石油大學(xué)（華東），2010.

［2］韓福來(lái).利用大物理模型研究特高含水期水驅(qū)規(guī)律［D］.青島：中國(guó)石油大學(xué)（華東），2012.

［3］姜漢橋，谷建偉，陳月明，等.剩余油分布規(guī)律的精細(xì)數(shù)值模擬［J］.中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2004，23（5）：30-34.

［4］王小林，于立君，唐瑋，等.特高含水期含水率與采出程度關(guān)系式［J］.特種油氣藏，2015，22（5）：104-106.

［5］王志偉，李相方.凝析氣相變微觀孔隙模型實(shí)驗(yàn)研究［J］.工程熱物理學(xué)報(bào)，2006，27（2）：251-254.

［6］曹永.駕駛員眼睛開(kāi)閉狀態(tài)計(jì)算機(jī)圖像識(shí)別技術(shù)開(kāi)發(fā)［D］.青島：青島大學(xué)，2011.

［7］董建國(guó).基于形狀的圖像目標(biāo)識(shí)別方法研究［D］.西安：西安電子科技大學(xué)，2009.

［8］劉金保.形狀匹配算法研究及應(yīng)用［D］.廣州：廣東工業(yè)大學(xué)，2011.

［9］白振強(qiáng)，吳勝和，付志國(guó).大慶油田聚合物驅(qū)后微觀剩余油分布規(guī)律［J］.石油學(xué)報(bào)，2013，34（5）：924-931.

［10］孫先達(dá).儲(chǔ)層微觀剩余油分析技術(shù)開(kāi)發(fā)與應(yīng)用研究［D］.長(zhǎng)春：吉林大學(xué)，2011.

［11］李崢嶸.圖像分割多閾值法在CT圖像重建中的應(yīng)用［D］.合肥：合肥工業(yè)大學(xué)，2007.

［12］SINGH T R，ROY S，SINGH O I，et al.A new local adaptive thresholding technique in binarization［J］.International Journal of Computer Science Issues，2011，8（6）：271-277.

［13］趙秀才.數(shù)字巖心及孔隙網(wǎng)絡(luò)模型重構(gòu)方法研究［D］.青島：中國(guó)石油大學(xué)（華東），2009.

［14］王超，王士同.有局部保持的最大間距準(zhǔn)則特征提取方法［J］.模式識(shí)別與人工智能，2009，22（6）：898-902.

（編輯孫薇）

Percolation characteristics investigation of microscopic remaining oil in water flooding reservoir with ultra-high water cut

YU Chunlei1，MI Lidong2，WANG Chuan2，ZHAO Yuyun2，JIANG Hanqiao2，TIAN Ye2
（1.Research Institute of Exploration and Development，Shengli Oilfield Company，SINOPEC，Dongying 257015，China；2.MOE Key Laboratory of Petroleum Engineering，China University of Petroleum，Beijing 102249，China）

The study of remaining oil percolation characteristics under microscopic pore scale plays an important role in enhanced oil recovery for ultra-high water cut period water flooding reservoir.In this paper，microscopic glass-etching model experiments and computer image processing&recognition techniques are combined to re-classify the microscopic remaining oil flowing patterns to analyze the flowing shape and the connect-relationship among the oil，water and rock.Research result indicates all the remaining oil in ultra-high water cut period can be classified into five categories，clustered stream，multi-porous stream，columnar stream，membranous stream and droplet stream.Among of the five categories，clustered stream possess the largest proportion and is also the main factor of relative permeability curves re-curved.Meanwhile，with water saturation rising，clustered stream gradually transforms into other patterns like multi-porous stream，columnar stream，membranous stream and droplet stream.The flow rules of the remaining oil and the reason of the non-linear percolation curve are explained microscopically，which provides a guiding for the development and recovery of water flooding reservoir with ultra-high water cut.

water-flooding reservoir；ultra-high water cut；computer image processing&recognition technique；percolation characteristics

國(guó)家科技重大專(zhuān)項(xiàng)“勝利油田特高含水期提高采收率技術(shù)”（2011ZX05011）

TE312

A

10.6056/dkyqt201605010

2016-01-23；改回日期：2016-07-19。

于春磊，男，1984年生，高級(jí)工程師，碩士，現(xiàn)從事油藏?cái)?shù)值模擬工作。E-mail：yuchunly@sina.com。

引用格式：于春磊，糜利棟，王川，等.水驅(qū)油藏特高含水期微觀剩余油滲流特征研究［J］.斷塊油氣田，2016，23（5）：592-594，598.

YU Chunlei，MI Lidong，WANG Chuan，et al.Percolation characteristics investigation of microscopic remaining oil in water flooding reservoir with ultra-high water cut［J］.Fault-Block Oil&Gas Field，2016，23（5）：592-594，598.