李南，田冀，譚先紅，劉新光

（中海油研究總院，北京100028）

合理開發(fā)低滲透儲層，已經(jīng)成為當(dāng)前以及今后緩解石油緊缺矛盾的重要手段，也對我國石油行業(yè)的持續(xù)穩(wěn)定發(fā)展具有重要意義［1-3］。由于針對我國低滲透儲層特征進行的CO2驅(qū)相關(guān)研究還處在初級階段，因此加強低滲透儲層CO2驅(qū)油波及特征的研究，顯得尤為重要［4-7］。低滲透儲層孔隙喉道細(xì)小，采用常規(guī)顯微鏡等識別較為復(fù)雜困難，應(yīng)用有限元算法和兩相滲流理論，結(jié)合低滲透油藏CO2滲流的擴散機理，可以有效地進行孔隙喉道的仿真模擬計算，從而精確得到微觀CO2驅(qū)的波及特征，為認(rèn)識低滲透油藏CO2滲流機理奠定了基礎(chǔ)。

1 考慮擴散微觀兩相滲流模型的建立

在考慮擴散的前提下，將CO2驅(qū)油簡化為油氣兩相流，著重研究不同的孔隙喉道、對流擴散系數(shù)下的微觀驅(qū)油特征及波及體積。

1）連續(xù)性方程：

式中：ρ 為流體密度，kg/m3；v 為流體滲流速度，m/s；p為地層壓力，MPa；g 為重力加速度，N/kg；t 為時間，s；μ流體黏度，mPa·s；下標(biāo)o，g 分別表示油、氣相。

2）擴散方程：

式中：cg為氣相濃度，mol/m3；Dc為氣相的對流擴散系數(shù)，m2/s。

3）黏度、密度矯正方程：

式中：S 為流體飽和度。

4）飽和度約束條件：

2 CO2 微觀波及體積的影響因素

前人在研究CO2驅(qū)時，多是宏觀分析對流擴散對產(chǎn)能的影響［8-18］，但沒有詳細(xì)研究對流擴散對微觀滲流的影響。本文通過微觀有限元模擬，以精確了解對流擴散對CO2驅(qū)微觀波及體積的影響程度。

2.1 對流擴散

2.1.1 簡單孔喉模型

基于有限元算法可以更為精確地計算復(fù)雜邊界條件下流體的運動規(guī)律，因此，建立考慮擴散微觀兩相滲流有限元數(shù)學(xué)模型，采用有限元三角形網(wǎng)格剖分方法，選用大規(guī)模稀疏矩陣PARDISO 進行求解。

CO2驅(qū)微觀流動模型采用考慮擴散的微觀兩相滲流模型，模擬CO2與原油沒有混相時的微觀流動，著重考慮不同對流擴散系數(shù)情況下對CO2波及體積的影響。模型參數(shù)如下：CO2地下黏度0.04 mPa·s，原油地下黏度2.70 mPa·s，地層壓力22 MPa，喉道半徑2 μm，CO2地下密度500 kg/m3，原油密度900 kg/m3，混相流體濃度3.5 mol/m3，模型尺寸20 μm×14 μm。模型示意見圖1。

對流擴散系數(shù)很小時（可以忽略不計），在孔喉的末端存在濃度滯后現(xiàn)象；當(dāng)考慮CO2與原油的對流擴散系數(shù)時，濃度滯后現(xiàn)象幾乎不存在了。這說明CO2混相后擴散現(xiàn)象會使流體濃度分布更均勻，同時能達(dá)到提高CO2波及體積及混相效果的目的（見圖1，色譜柱以下同）。

圖1 不同擴散系數(shù)CO2 濃度分布

2.1.2 復(fù)雜孔喉模型

為了研究在復(fù)雜孔喉中對流擴散作用對CO2微觀滲流的影響，在二維簡單孔喉模型基礎(chǔ)上，建立了CO2驅(qū)微觀復(fù)雜孔喉流動模型。按3 種方案分別進行CO2驅(qū)油的微觀波及特征分析：方案1，均勻孔喉，不考慮對流擴散；方案2，非均勻孔喉，不考慮對流擴散；方案3，非均勻孔喉，考慮對流擴散。網(wǎng)格劃分采用自由剖分的三角形網(wǎng)格，采用瞬態(tài)計算，得到不同時間的各種參數(shù)場圖（見圖2）。

圖2 不同方案參考點CO2 濃度分布

由于氣液的黏度差，在多孔介質(zhì)中存在CO2突進的現(xiàn)象。方案2 與方案1 相比，由于孔隙喉道的非均勻性，CO2的波及效果明顯變差；方案3 中，當(dāng)考慮對流擴散后，從相同時間下CO2濃度的分布可以看出，對流擴散作用明顯提高了CO2的波及體積，即對流擴散作用降低了平面各向異性的影響。

2.2 潤濕性

針對低滲透儲層提高CO2波及體積困難的問題，本文采用微尺度仿真孔喉模型，在上述模擬的基礎(chǔ)上，同時考慮微觀擴散作用和潤濕性，分析微觀波及特征和模式，研究低滲透油藏提高CO2微觀波及體積的方法。模型參數(shù)：擴散系數(shù)1×10-8m2/s，模型尺寸32 μm×22 μm，界面張力1.7 mN/m，接觸角分別為30，90，120°，喉道半徑3 μm，含油飽和度大于0.3。

結(jié)合油藏實際的薄片實驗，在考慮邊界潤濕性的條件下，可以明顯發(fā)現(xiàn)，孔隙喉道內(nèi)部剩余油的微觀分布模式主要分為3 種： 在孔隙喉道較窄的區(qū)域出現(xiàn)柱狀剩余油、 在孔隙喉道突出的區(qū)域出現(xiàn)楔形剩余油和在孔隙喉道邊界層出現(xiàn)膜狀剩余油（見圖3）。

圖3 微觀剩余油模式

在相同的時間、不同邊界潤濕性的條件下，模型中的微觀波及體積明顯不同。通過計算發(fā)現(xiàn)：當(dāng)接觸角為30°時，CO2微觀波及體積62%；當(dāng)接觸角為90°時，CO2微觀波及體積53%；當(dāng)接觸角為120°時，CO2微觀波及體積46%（見圖4）。因此，在低滲透儲層孔喉的末端或狹窄的區(qū)域，可以通過轉(zhuǎn)變儲層的潤濕性提高CO2微觀波及體積。

圖4 微觀CO2 濃度分布

3 結(jié)論

1）分別建立了考慮擴散的低滲透簡單、復(fù)雜的孔喉微觀模型。通過微觀有限元模擬認(rèn)為，對流擴散提高了低滲透油藏CO2驅(qū)的微觀波及體積。

2）通過對仿真孔隙喉道的模擬驅(qū)油發(fā)現(xiàn)，其內(nèi)部剩余油的微觀分布模式主要有柱狀、楔形、膜狀。

3）考慮不同潤濕性邊界時，隨著潤濕角的增加，CO2微觀波及體積明顯降低。

［1］孫長宇，王文強，陳光進，等.注CO2油氣藏流體體系油/水和油/氣界面張力實驗研究［J］.中國石油大學(xué)學(xué)報：自然科學(xué)版，2006，30（5）：109-112.

［2］馬濤，湯達(dá)禎，蔣平，等.注CO2提高采收率技術(shù)現(xiàn)狀［J］.油田化學(xué)，2007，24（4）：379-383.

［3］丁景辰.裂縫性油藏注CO2物質(zhì)平衡方程的建立與應(yīng)用［J］.斷塊油氣田，2014，21（2）：213-216.

［4］鄧波，李鴻，曹建，等.注CO2驅(qū)油藏先導(dǎo)性篩選評價方法［J］.西南石油大學(xué)學(xué)報：自然科學(xué)版，2009，31（2）：105-108.

［5］王進安，袁廣均，張軍，等.長巖心注二氧化碳驅(qū)油物理模擬實驗研究［J］.特種油氣藏，2001，8（2）：75-78.

［6］何應(yīng)付，高慧梅，周錫生，等.改善特低滲透油藏注二氧化碳驅(qū)油效果的方法［J］.斷塊油氣田，2011，18（4）：512-514.

［7］張廣杰，武若霞.改善新民裂縫性砂巖油藏注水波及體積研究［J］.石油勘探與開發(fā)，1998，25（3）：88-91.

［8］廖海嬰.腰英臺DB34 井區(qū)CO2驅(qū)替油藏數(shù)值模擬研究［J］.西安石油大學(xué)學(xué)報：自然科學(xué)版，2010，25（5）：50-53.

［9］劉中云，曾慶輝，唐周懷.潤濕性對采收率及相對滲透率的影響［J］.石油與天然氣地質(zhì)，2000，21（2）：148-150.

［10］莊永濤，劉鵬程，郝明強，等.低滲透油藏CO2驅(qū)井網(wǎng)模式數(shù)值模擬［J］.斷塊油氣田，2013，20（4）：477-480.

［11］王建波，高云叢，王科戰(zhàn).腰英臺特低滲透油藏CO2驅(qū)油井見氣規(guī)律研究［J］.斷塊油氣田，2013，20（1）：118-122.

［12］鐘立國，韓大匡，李莉.人工壓裂低滲透油藏二氧化碳吞吐因素敏感性分析［J］.大慶石油學(xué)院學(xué)報，2009，33（4）：125-129.

［13］吳忠寶，甘俊奇，曾倩.低滲透油藏二氧化碳混相驅(qū)油機理數(shù)值模擬［J］.油氣地質(zhì)與采收率，2012，19（3）：67-70.

［14］郭平，黃宇，李向良.滲透率及壓力對低滲油藏CO2驅(qū)油效率的影響［J］.斷塊油氣田，2013，20（6）：768-771.

［15］湯勇，孫雷，周涌沂，等.注富烴氣凝析/蒸發(fā)混相驅(qū)機理評價［J］.石油勘探與開發(fā)，2005，32（2）：133-136.

［16］章星，楊勝來，文博，等.低滲油藏CO2混相驅(qū)啟動壓力梯度實驗研究［J］.石油實驗地質(zhì)，2013，35（5）：583-586.

［17］仵元兵，胡丹丹，常毓文，等.CO2驅(qū)提高低滲透油藏采收率的應(yīng)用現(xiàn)狀［J］.新疆石油天然氣，2010，6（1）：36-39.

［18］張奉東，王震亮.蘇北盆地草舍油田CO2混相驅(qū)替試驗與效果分析［J］.石油實驗地質(zhì)，2010，32（3）：296-300.