低滲油藏活性納米顆粒減阻增注解析模擬方法

劉旭鵬，袁 彬，2，3，戴彩麗，2，3，李 躍，韓明亮

1)中國石油大學(xué)(華東)石油工程學(xué)院，山東青島266580；2)中國石油大學(xué)(華東)非常規(guī)油氣開發(fā)教育部重點(diǎn)實(shí)驗(yàn)室，山東青島266580；3)山東省油田化學(xué)重點(diǎn)實(shí)驗(yàn)室，山東青島266580

納米顆粒指粒徑介于1～100 nm之間的納米粒子材料，具有高表面活性、高化學(xué)穩(wěn)定性以及高比表面積等優(yōu)異性能，在石油行業(yè)得到廣泛的研究與應(yīng)用[1-2]．LI等[3]利用納米二氧化硅制備了親水型納米顆粒，通過室內(nèi)實(shí)驗(yàn)研究了親水納米顆粒在提高采收率方面的效果．ZHAO等[4]通過將納米二氧化硅表面改性，制備了疏水型納米顆粒，并對(duì)其降壓增注效果及機(jī)理進(jìn)行了相應(yīng)的研究．除此之外，MOHAMADIAN等[5]發(fā)現(xiàn)納米材料的加入使鉆井液流變性能和過濾性能得到了顯著的提升．納米技術(shù)對(duì)石油行業(yè)的發(fā)展產(chǎn)生了極大的生產(chǎn)力和推動(dòng)力．

低滲油藏現(xiàn)已成為中國油氣資源的主力軍．但低滲油藏存在孔隙度小、產(chǎn)能低和注水壓力大等開發(fā)難題，必須經(jīng)過相應(yīng)措施提高開發(fā)效率．目前，低滲油藏降壓增注方法主要包括儲(chǔ)層改造和化學(xué)改性．其中，化學(xué)改性法是通過表面修飾材料改變儲(chǔ)層孔隙表面潤濕性，降低流動(dòng)阻力，提高注水和生產(chǎn)能力[6-9]．

針對(duì)納米顆粒降壓增注技術(shù)，國內(nèi)外都進(jìn)行了廣泛實(shí)驗(yàn)研究[10-11]．疏水型納米顆粒能夠在巖石孔隙或吼道表面牢固吸附，形成具有疏水性的顆粒吸附層，取代巖石表面水膜，改變巖石孔隙表面潤濕性，降低油水界面張力，同時(shí)也避免了黏土的水化膨脹，有效解決了低滲油藏高壓欠注難題[12-15]．然而，目前對(duì)活性納米顆粒降壓增注能力的研究僅停留在改變界面張力和潤濕性等宏觀特性，尚缺乏系統(tǒng)表征低滲儲(chǔ)層內(nèi)活性納米顆粒滲流及降壓增注的模擬方法，不能有效指導(dǎo)活性納米顆粒降壓增注方案設(shè)計(jì)，限制了它的推廣應(yīng)用．為此，本研究以毛細(xì)管滲流模型為基礎(chǔ)，在考慮孔隙壁面粗糙條件下，模擬了活性納米顆粒的在滲流過程中所發(fā)生的動(dòng)態(tài)行為，分析了活性納米顆粒所造成的正負(fù)面效應(yīng)，建立了活性納米顆粒降壓增注模型，該模型定義的參數(shù)及參考取值見表1．該模型對(duì)活性納米流體的應(yīng)用有一定的指導(dǎo)意義．

表1 本研究模型定義的參數(shù)及參考取值

1 模型建立

本研究對(duì)經(jīng)典毛細(xì)管束模型進(jìn)行表面“粗糙化”處理，提出了周期性矩形凹槽結(jié)構(gòu)的微管模型，如圖1．其中，x為水流方向的距離；ln為法向力臂；ld為拖曳力臂．以單周期壁面為例，單周期水平長度為M， 將壁面分為4個(gè)區(qū)域： 區(qū)域A和C為水平區(qū)域，長度分別為a和b； 區(qū)域B和D是斜面，水平長度為l且傾角為θ． 模型滿足以下假設(shè)：① 一維多孔介質(zhì)為均質(zhì)，處于熱力學(xué)平衡，初始條件孔隙表面水濕，水和納米顆粒兩相等溫流動(dòng)；② 活性納米流體化學(xué)性質(zhì)穩(wěn)定，不與流體或界面發(fā)生反應(yīng)；③ 流體流速足夠大，可忽略顆粒擴(kuò)散影響；④ 注水井周圍含水孔隙度為(1-Sor)φ； ⑤ 孔隙內(nèi)納米顆粒受到4種力作用：流動(dòng)黏性流體造成的拖曳力Fd和升力Fl、 壁面與納米顆粒間靜電力Fe以及顆粒自身重力Fg[16]．

圖1 周期性矩形槽型結(jié)構(gòu)模擬孔隙粗糙壁面Fig.1 Pore rough wall simulated by periodic rectangular groove structure

(1)

定義脫落因子y為拖曳力與靜電力的比值,即

(2)

其中，ri(i=A、B、C、D)分別為4個(gè)區(qū)域處納米顆粒吸附后的微管半徑；Fdi為4個(gè)區(qū)域處(A、B、C、D)納米顆粒所受拖曳力;r為流體流動(dòng)半徑；μ為流體黏度；v為納米顆粒的注入流速．

結(jié)合式(1)和式(2)可得4個(gè)區(qū)域處微管半徑ri與吸附層厚度rci之間的關(guān)系，

(3)

(4)

其中，N為截面上的毛細(xì)管總數(shù)．

結(jié)合式(3)、式(4)以及巖心孔隙度與毛細(xì)管參數(shù)的關(guān)系式，得到4個(gè)區(qū)域處微管表面納米顆粒的吸附體積為

(5)

最后，單位周期上微管壁面上納米顆粒的吸附體積為VA+VB+VC+VD， 則整個(gè)巖心內(nèi)納米顆粒的吸附體積分?jǐn)?shù)可表示為

(6)

圖2為活性納米顆粒在孔隙表面極限附著體積分?jǐn)?shù)與流速之間的關(guān)系．由圖2可知，活性納米顆粒在巖石顆粒表面的極限附著體積分?jǐn)?shù)隨著流速的增大逐漸變小，即流速越大，流體所帶來的拖曳力越大，使得順時(shí)針方向的脫離力矩增大，納米顆粒從壁面脫落下來，使得吸附層變薄．此外，微管表面的粗糙度也對(duì)納米顆粒的吸附體積分?jǐn)?shù)有較大影響，與表面光滑的微管相比，周期性矩形槽型結(jié)構(gòu)壁面附著的活性納米顆粒要更多．

圖2 活性納米顆粒在孔隙表面極限附著體積分?jǐn)?shù)與流速之間的關(guān)系Fig.2 The relationship between the limit adhesion concentration of active nanoparticles on the pore surface and the flow rate

2 活性納米顆粒傳輸-附著-堵塞耦合過程模擬

特征線法(method of characteristics, MOC)又稱行波法或達(dá)朗貝爾法，是一種解決一階線性和二階雙曲型偏微分方程的常用方法．YUAN等[17-19]利用MOC對(duì)納米流體減緩微米顆粒運(yùn)移方面的問題作了詳細(xì)研究，最終得到了精確的解析解，從微觀層面充分地揭示了納米流體對(duì)減緩儲(chǔ)層中微粒運(yùn)移的過程．因此，本研究采用特征線法對(duì)活性納米顆粒在巖心中3種行為機(jī)制進(jìn)行研究表征．

無因次的活性納米顆粒的質(zhì)量守恒方程為

(7)

由經(jīng)典過濾理論中顆粒捕捉動(dòng)力學(xué)方程，可得到活性納米顆粒在巖心中附著和堵塞的速率表達(dá)式[12]為

(8)

式(7)活性納米顆粒的質(zhì)量守恒方程為典型的一階線性的偏微分方程，為此，在給定初始條件和邊界條件的情況下，采用特征線法可以將定解條件下的偏微分方程的求解變?yōu)槌Ｎ⒎址匠痰那蠼猓?/p>

本研究將活性納米流體在巖心中流動(dòng)時(shí)所產(chǎn)生的3種行為機(jī)制體積分?jǐn)?shù)變化(懸浮、附著和堵塞)分為4個(gè)階段(圖3)：① 活性納米流體注入巖心前(初始條件,φ0為注入體積分?jǐn)?shù))；② 入口處活性納米顆粒附著體積分?jǐn)?shù)達(dá)到極限的過程；③ 巖心內(nèi)部活性納米顆粒附著體積分?jǐn)?shù)達(dá)到極限的過程；④ 整個(gè)巖心內(nèi)活性納米顆粒附著體積分?jǐn)?shù)達(dá)到極限后．

初始條件：

邊界條件：

結(jié)合式(7)和式(8)，得到活性納米顆粒懸浮、附著和堵塞3種行為機(jī)制的質(zhì)量守恒方程為

(9)

以活性納米顆粒懸浮體積分?jǐn)?shù)偏微分方程的求解為例，利用特征線法將式(9)中第1個(gè)式子轉(zhuǎn)化為常微分方程的形式，即

(10)

對(duì)式(10)中的xD進(jìn)行積分，并結(jié)合邊界條件，得到階段2中懸浮納米顆粒的體積分?jǐn)?shù)方程為

φ1=φ0exp[-(λa+λs)LxD]

(11)

利用同樣的方法，得到階段 2中活性納米顆粒附著和堵塞體積分?jǐn)?shù)的表達(dá)式分別為

φ2=λaφ(1-Sor)L(tD-xD)·

exp[-(λa+λs)LxD]

(12)

φ3=λsφ(1-Sor)L(tD-xD)·

exp[-(λa+λs)LxD]

(13)

若巖心入口處的活性納米顆粒附著體積分?jǐn)?shù)在tDC時(shí)刻達(dá)到極限，則tDC的表達(dá)式為

(14)

以侵蝕前緣(xcr,tcr)表示巖心中活性納米顆粒到達(dá)極限附著體積分?jǐn)?shù)的位置和時(shí)刻，則侵蝕前緣前后兩端的活性納米顆粒體積分?jǐn)?shù)由連續(xù)性條件得φ-=φ+=φ0． 因此，結(jié)合式(8)和式(9)，得到納米顆粒附著體積分?jǐn)?shù)距離的微分表達(dá)式為

(15)

對(duì)活性納米顆粒的極限附著體積分?jǐn)?shù)的時(shí)間tD進(jìn)行全微分可得

(16)

結(jié)合式(8)、式(15)和式(16)，得到活性納米顆粒附著體積分?jǐn)?shù)侵蝕前緣時(shí)間與距離的表達(dá)式為

(17)

同樣，利用特征線法，結(jié)合式(9)及此時(shí)的邊界條件得到階段 3中活性納米顆粒懸浮、附著和堵塞的體積分?jǐn)?shù)為

(18)

隨著活性納米流體的進(jìn)一步注入，設(shè)在tcr1(此時(shí)xD=1)時(shí)整個(gè)巖心內(nèi)活性納米顆粒的附著體積分?jǐn)?shù)達(dá)到極限，因此出口處的邊界條件為φ(1,tcr1)=φmax． 階段4中納米顆粒懸浮、附著和堵塞的體積分?jǐn)?shù)為

(19)

以φ0=2%， 穩(wěn)定注入速率v=0.008 m/s為例，分別選取tD為0.5、10.0、20.0和35.0時(shí)，觀察活性納米顆粒的懸浮、附著以及堵塞的體積分?jǐn)?shù)變化．

圖4為特征線圖像，當(dāng)tD=0.5時(shí)，對(duì)應(yīng)圖3和圖4中的階段2，由于巖石孔隙表面上納米顆粒的附著體積分?jǐn)?shù)未達(dá)極限，壁面上納米顆粒的附著數(shù)量會(huì)不斷增加，直至到達(dá)tDC=1.04時(shí)，入口納米顆粒的附著量達(dá)到飽和．結(jié)合圖5(b)可以看出，tD=0.5時(shí)納米顆粒附著體積分?jǐn)?shù)并未到達(dá)峰值， 其中，φmax=6.95%． 隨著納米流體的進(jìn)一步注入，在tD分別為10.0和20.0時(shí)，活性納米顆粒在巖心內(nèi)進(jìn)一步飽和，此時(shí)巖心中經(jīng)納米顆粒侵蝕過的區(qū)域的附著量已達(dá)到飽和，當(dāng)tD分別為10.0和20.0時(shí)，活性納米流體流經(jīng)的區(qū)域都已達(dá)到極限附著體積分?jǐn)?shù)．當(dāng)tD=35.0時(shí)，侵蝕前沿到達(dá)巖心出口處，此時(shí)整個(gè)多孔介質(zhì)納米顆粒的附著體積分?jǐn)?shù)已到達(dá)極限，在圖4中表示為階段4，即整個(gè)巖心不會(huì)再吸附更多的納米顆粒，僅發(fā)生納米顆粒堵塞．

圖4 特征線圖像Fig.4 Characteristic line

圖5 不同時(shí)刻巖心內(nèi)φ1、φ2和φ3的變化Fig.5 The concentration of active nanoparticles suspended, attached and blocked in the core at different times

3 活性納米顆粒降壓增注效果驗(yàn)證

3.1 降壓模型推導(dǎo)

如圖6所示，活性納米流體的注入使得巖心壁面潤濕性由親水轉(zhuǎn)化為疏水或超疏水，產(chǎn)生疏水滑移效應(yīng)，增大了毛細(xì)管的有效半徑．滑移長度的求解借助WU等[20]推導(dǎo)的潤濕性與滑移長度ls之間關(guān)系的表達(dá)式

(20)

其中，re為微管產(chǎn)生疏水滑移效應(yīng)后的半徑，即微管的有效半徑；Cw為所用液體通過實(shí)驗(yàn)或分子動(dòng)力學(xué)(moleculardynamics,MD)模擬所得的液體常數(shù)；θ0為潤濕角．

圖6 納米顆粒吸附層表面疏水滑移效應(yīng)Fig.6 The hydrophobic slip effect on the surface of the nanoparticle adsorption layer

如圖5(b)所示，在活性納米流體注入過程中，整個(gè)巖心會(huì)以侵蝕前緣的位置xcr為分界點(diǎn)，分為滲透率不同的兩個(gè)區(qū)塊．以tD=20.0為例，此時(shí)區(qū)塊1孔隙表面的活性納米顆粒附著后體積分?jǐn)?shù)已達(dá)到飽和，滲透率也由原始滲透率k0轉(zhuǎn)變?yōu)閗1， 滲透率增大；而區(qū)塊2由于未能附著足夠多的活性納米顆粒，其滲透率也未發(fā)生改變，可視為初始滲透率k0．當(dāng)整個(gè)巖心內(nèi)的活性納米顆粒附著體積分?jǐn)?shù)都已到達(dá)極限時(shí)，即注入活性納米流體后巖心兩端垂直表面壓強(qiáng)差Δp1可表示為

(21)

其中，βa和βs分別為活性納米顆粒的附著和堵塞儲(chǔ)層損害系數(shù)．

定義降壓率D為注入活性納米流體前、后巖心兩端垂直表面壓強(qiáng)差(Δp0和Δp1)之差與注入活性納米流體前巖心兩端壓降之比,即

(22)

由式(22)可知，降壓率D越大，活性納米流體降低壓力在巖心內(nèi)的損耗效果越好；反之，效果越差．

圖7 潤濕性的影響Fig.7 Influence of wettability

3.2 活性納米流體選用分析及降壓增注效果評(píng)價(jià)

3.2.1 納米顆?；钚詮?qiáng)弱對(duì)降壓效果影響

如圖7所示，分別設(shè)潤濕角θ0為100°、120°和150°表征選用活性納米顆粒的活性(潤濕性)的強(qiáng)弱．通過與未使用活性納米流體的結(jié)果比較，3種不同活性的納米流體都起到了一定降壓的作用，并且潤濕角越大，降壓效果越好．即所選納米顆?；钚栽綇?qiáng)，降壓效果越明顯．在θ0=120°時(shí)降壓率達(dá)到40%左右．

3.2.2 納米顆粒注入體積分?jǐn)?shù)對(duì)降壓效果影響

通過設(shè)置納米顆粒注入體積分?jǐn)?shù)φ0分別為2%、3%和4%，模擬注入體積分?jǐn)?shù)對(duì)納米顆粒減壓效果的影響(圖8)．在3種不同注入體積分?jǐn)?shù)情況下，注入體積分?jǐn)?shù)為2%和3%都起到了一定的降壓作用．相反，注入體積分?jǐn)?shù)為4%時(shí)壓降大于未使用活性納米顆粒的壓降，但在注入過程中壓降也有下降趨勢．

圖8 注入體積分?jǐn)?shù)的影響 Fig.8 Influence of injection volume fraction

圖9 注入流速的影響Fig.9 Influence of injection flow rate

3.2.3 納米顆粒注入速率對(duì)降壓效果影響

以v=0.004、0.006、0.008和0.010m/s向巖心內(nèi)注入納米流體時(shí)，降壓率分別為37.9%、39.0% 、40.5%和42.5%(圖9)．通過比較發(fā)現(xiàn)，注入速率在0.01m/s時(shí)所起到的降壓效果最好．結(jié)合圖3可知，由于注入速率的減小會(huì)使得納米顆粒極限附著體積分?jǐn)?shù)增加，進(jìn)而使得壁面吸附層厚度增大，影響巖心滲透率和孔隙度的大小；相反，如果注入速率過大，則壁面不能吸附足夠的活性納米顆粒，從而不能達(dá)到出色的降壓效果．

3.2.4 納米顆粒對(duì)注入效果影響

本研究通過比較在納米顆粒注入量相同情況下注入指數(shù)J的大小，來表征活性納米顆粒對(duì)注入效果的影響．注水指數(shù)的表達(dá)式為

(23)

其中，Q為納米流體(或水)注入流量．

以注入流速為0.008m/s時(shí)為例，在相同注入流速的情況下，注入活性納米流體后的巖心的注水指數(shù)要明顯大于未使用納米流體的注水指數(shù)，即單位壓差下的注入流量更大(圖10)，由此可知，活性納米流體的注入顯著增強(qiáng)了巖心的注入能力．

圖10 活性納米流體注入前后注水指數(shù)變化Fig.10 Water injection index before and after active nanofluid injection

4 結(jié) 論

1)活性納米顆粒在巖石孔隙表面的附著是在多種力平衡作用下的結(jié)果，在保證礦化度和活性納米流體類型不變的情況下，活性納米顆粒的極限附著體積分?jǐn)?shù)隨注入速率的增大而減小．

2)利用活性納米流體降壓增注模型對(duì)影響降壓增注效果的納米流體活性、注入體積分?jǐn)?shù)和注入速率3個(gè)因素進(jìn)行計(jì)算分析．結(jié)果顯示，計(jì)算值與實(shí)際巖心測試時(shí)所得到的納米顆粒降壓增注效果具有較高的一致性．

3)活性納米流體降壓增注模型綜合考慮活性納米流體注入巖心時(shí)產(chǎn)生的正、負(fù)面效應(yīng)，從微觀角度得出活性納米降壓增注作用機(jī)理，為實(shí)際生產(chǎn)過程中納米顆粒的優(yōu)選提供理論參考．