氣固兩相流下砂巖co

0含束縛水飽和度對滑脫效應(yīng)的影響隨著低滲氣田的開發(fā)和氣驅(qū)技術(shù)的發(fā)展，氣體泄漏規(guī)律的研究越來越受到重視?；撔?yīng)是氣體滲透流的主要影響因素之一。氣驅(qū)流規(guī)律的研究和氣藏出井的生產(chǎn)能力計算非常重要。氣體滲流性質(zhì)僅與自身的飽和度有關(guān),不受液相性質(zhì)影響[2,3,4,5,6,7,8,9,10,11].因此,需要研究氣體滑脫效應(yīng)與氣相飽和度之間的關(guān)系,目前存在兩種截然相反的觀點:一種觀點認(rèn)為氣體滑脫效應(yīng)隨氣相飽和度的增加而減小.KlinkenbergLJ、萬軍鳳等認(rèn)為,由于氣—液的分子間作用力比氣—固的小,含束縛水時管壁的氣體分子受到的束縛力降低,氣體滑脫效應(yīng)隨束縛水飽和度的增加更加明顯;LiKewen等在20~170℃溫度內(nèi)進行N2—水和蒸汽—水穩(wěn)態(tài)法實驗,得到與克氏理論一致的結(jié)果.另一種觀點認(rèn)為氣體滑脫效應(yīng)隨氣相飽和度的增加而增加.RoseWD通過人造巖心和天然巖心水氣兩相非穩(wěn)態(tài)驅(qū)替實驗,發(fā)現(xiàn)滑脫因子隨氣相飽和度的增加而增加;FultonPF通過在含水飽和度低于30%的條件下滑脫因子測定實驗得到同樣的結(jié)論.RushingJA等利用穩(wěn)態(tài)法和非穩(wěn)態(tài)法,研究不同含水飽和度條件下低滲巖心的滑脫效應(yīng),在滲透率為(0.01~0.10)×10-3μm2、含水飽和度為5%~40%條件下滑脫效應(yīng)隨含水飽和度增加而降低.得到與克氏理論一致結(jié)果的實驗為兩相穩(wěn)態(tài)驅(qū)替實驗,得到相反結(jié)果的實驗處于較高氣相飽和度范圍內(nèi).肖曉春等針對煤層滲流,在滲透率為0.3×10-6~0.2×10-3μm2、含水飽和度為29.5%~51.3%條件下測定滑脫因子,發(fā)現(xiàn)圍壓為4~6MPa時滑脫效應(yīng)隨氣相飽和度的增加而減小;圍壓為8~10MPa時滑脫效應(yīng)隨氣相飽和度的增加而增加.筆者設(shè)計驅(qū)替實驗,在氣相飽和度為24.45%~100%條件下測定滑脫因子,研究氣相飽和度對滑脫效應(yīng)的影響,并通過微觀實驗解釋滑脫效應(yīng)隨氣相飽和度變化的原因.1氣體滑脫效應(yīng)KlinkenbergLJ通過實驗發(fā)現(xiàn)氣體滑脫效應(yīng),提出考慮氣體滑脫效應(yīng)的滲透率表達式:式中:Kg為含有滑脫效應(yīng)的氣測滲透率;Kg0為絕對滲透率;b為氣體滑脫因子;pc為驅(qū)替壓力.其中,氣體滑脫因子b為式中:c為比例因子;λ為氣體分子平均自由程,受溫度和壓力影響;rc為氣體流動通道半徑;Kn為克努森數(shù),其值越大氣體滑脫效應(yīng)越大.2氣體滑動脫離因子的測定2.1高相沖2.1.1氣體滑脫效應(yīng)實驗所用天然巖心和地層水來自延長油田,地層水組成見表1.所用氣體為濃度大于99.9%的CO2.在地層條件下,CO2常處于氣態(tài)(溫度T>31.1℃、壓力p<7.382MPa)或超臨界狀態(tài)(T>31.1℃、p>7.382MPa),在兩種狀態(tài)下CO2存在氣體滑脫效應(yīng).用CO2氣體測定107塊取心巖心的孔隙度、滲透率及滑脫因子,統(tǒng)計滑脫因子與滲透率、孔隙度之間的關(guān)系,擬合結(jié)果為用DO7氣體流量計計量氣體流量;用JYB-K壓變傳感器計量巖心兩端的壓力差.2.1.2壓力差及巖心滲透率壓力變化影響分子自由程和滑脫因子.為了降低壓力變化的影響,在較低壓力下進行實驗:(1)將抽提后的天然巖心烘干、稱重,并用氣態(tài)CO2以50.0、100.0、200.0、300.0、400.0、500.0kPa的壓力差測定巖心滲透率(根據(jù)靖邊采油廠某區(qū)塊地層條件,設(shè)定溫度為45℃、回壓為6MPa);(2)抽真空并飽和地層水后稱重,分別以0.2、0.5、1.0、1.5、2.0、3.0、4.0、5.0mL/min的流速測定水測滲透率,比較相同流速下的氣測和水測滲透率,計算滑脫因子;(3)氣驅(qū)巖心至不再出水并稱重,在殘余水飽和度下測定氣測滲透率(與步驟(1)相同),計算氣相飽和度及滑脫因子;(4)烘干巖心至不同含水飽和度、稱重,測定氣測滲透率(與步驟(1)相同),計算氣相飽和度及氣體滑脫因子.2.1.3在古代環(huán)境下氣體滑脫效應(yīng)針對滲透率為25.7×10-3μm2的延長油田巖心,氣測和水測滲透率結(jié)果見圖1(a),通過式(1)計算滑脫因子b為149.6kPa.在不同氣相飽和度下氣測滲透率部分結(jié)果見圖1(b),通過式(1)的導(dǎo)函數(shù)計算滑脫因子:得到氣相飽和度Sg為89.02%、80.26%、71.29%對應(yīng)的氣體滑脫因子分別為150.3、150.7、161.4kPa.近似認(rèn)為圖1(b)中各壓力曲線較高點的趨勢線與縱坐標(biāo)的交點是滲透率Kg0(即壓力極大時的滲透率).由圖1(b)可以看出,隨著含水飽和度的增加,氣體滑脫效應(yīng)逐漸增加.在含水飽和度下水相為束縛水,占據(jù)孔隙多為無效孔隙,孔隙對滲流的貢獻率很低;對滑脫效應(yīng)的主要影響來自水膜的作用,影響相對較小.2.2低相量2.2.1氣體滑脫因子的測定CO2不僅能降低油相黏度,對糖溶液、甘油等高黏流體也具有降黏作用.因此,選擇改性淀粉溶液調(diào)剖劑測定低氣相飽和度下的氣體滑脫因子.改性淀粉調(diào)剖劑作為強凝膠體系,主要成分包括改性淀粉、丙烯酰胺單體、交聯(lián)劑和引發(fā)劑.該調(diào)剖體系在凝膠前是小分子溶液,注入性好,可以模擬正常的氣液分布;在凝膠后強度大,且CO2對體系的凝膠強度幾乎沒有影響,可以在一定程度上保留氣體竄流通道形態(tài).2.2.2氣竄通道的恢復(fù)當(dāng)氣相飽和度較低時,凝膠后的氣體通過巖心時所受阻力大于凝膠前的;原因是氣體停止驅(qū)替后,較細(xì)的氣竄通道在卡斷效應(yīng)下被截成小段,在凝膠后將油水分布狀態(tài)保留下來,通過注水突破封堵,可以恢復(fù)氣竄通道.實驗步驟:(1)用CO2驅(qū)替飽和過凝膠體系巖心,至不再有液相采出,溫度升至80℃,并候凝8h;(2)用地層水注入巖心,氣驅(qū)至不再出水,稱重并計算氣相飽和度,在殘余水飽和度下測定氣測滲透率,用式(3)計算氣體滑脫因子;(3)將黏度為50.00、100.00、150.00、200.00mPa·s的改性淀粉溶液過濾并注入巖心,用CO2驅(qū)替改性淀粉溶液至不再出液,計量出液量并計算氣相飽和度,用式(3)計算氣體滑脫因子.2.2.3含氣飽和度對滑脫因子的影響?zhàn)ざ燃s為35.30mPa·s的凝膠體系占據(jù)的巖心體積,與地層水條件下的束縛水體積相近,且成膠后的滑脫因子也相似,說明強凝膠體系能夠模擬CO2驅(qū)替過程的氣液兩相分布,并能夠?qū)⑺Ａ粝聛?見圖2).利用高黏流體制造束縛水,使含氣飽和度迅速下降,且滲透率急劇降低.實驗得到的氣相飽和度分別為41.37%、32.64%、28.25%、24.45%,所對應(yīng)的滑脫因子分別為267.2、337.7、321.4、190.4kPa.這說明在低氣相飽和度條件下,滑脫因子隨氣相飽和度的變化規(guī)律與高氣相飽和度條件下的完全不同(見圖2).當(dāng)氣相飽和度低于束縛氣相飽和度時不再發(fā)生氣體滑脫效應(yīng).在低氣相飽和度條件下,隨著氣相飽和度增加,滑脫因子先增加后降低.當(dāng)水相變?yōu)槭`水時,滑脫因子隨著氣相飽和度增加而降低,但氣體滑脫效應(yīng)的變化幅度較小.3微觀驅(qū)替實驗通過微觀驅(qū)替實驗分析氣體滑脫因子隨氣相飽和度變化的原因.3.1模擬蒸餾所需水的組分采用南通華興石油儀器有限公司定制的微觀驅(qū)替模型,可放大1000倍,承壓上限為50MPa,含有保溫裝置,設(shè)備內(nèi)置復(fù)合式光源,其中透射光源為面熒光燈,反射光源為正交照射光源.微觀模型是以延長油田天然巖心的鑄體薄片為依據(jù)刻蝕而成的玻璃模型.實驗所用油、水來自延長油田,地層水組成見表1.原油密度為0.858t/m3、脫氣油黏度為11.54mPa·s,含氣原油黏度為4.87mPa·s,模擬蒸餾得到的脫氣原油碳數(shù)組分見表2.3.2水驅(qū)模型的建立文獻[4-10]利用穩(wěn)態(tài)法和非穩(wěn)態(tài)法研究滑脫效應(yīng)與氣相飽和度的關(guān)系,并得到不同的結(jié)果,文中采用不同驅(qū)替方式的微觀實驗研究滑脫效應(yīng).為了保持CO2的氣體自由程不變,實驗采用恒壓驅(qū)替:(1)用真空泵抽出各組模型空氣,直接飽和延長原油.(2)第一組實驗通過控制氣瓶以10.0kPa的壓力恒壓注入CO2,觀察驅(qū)替過程.(3)第二組實驗首先用微量計量泵以10.0μL/min的流速將水注入水驅(qū)模型,至含水率達到90%;然后通過控制氣瓶以10.0kPa的壓力恒壓注入CO2,觀察驅(qū)替過程.(4)第三組實驗通過控制氣瓶以10.0kPa的壓力恒壓注入CO2,同時用微量計量泵以5.0μL/min的流速將水注入水驅(qū)模型,觀察驅(qū)替過程.3.3結(jié)果與分析3.3.1注氣階段ab氣體滑脫效應(yīng).在CO2驅(qū)替初期,在注氣壓力下油相被啟動并開始運移,注入氣體在毛管中間隨著油相運動.由于氣體的黏滯力小于油相的,注入氣體在卡斷效應(yīng)下被截成小段,部分氣體處于束縛氣狀態(tài)(見圖3(a)).此時氣相不再是連續(xù)相,自由運動的氣體分子受到氣液界面的阻隔,從而嚴(yán)重削弱氣體滑脫效應(yīng).隨著注氣進行,CO2逐漸形成竄流通道,隨著竄流通道數(shù)量增多,氣體滑脫效應(yīng)逐漸增大(見圖3(b));在CO2驅(qū)替中后期,由于CO2竄流通道逐漸擴大,氣體分子自由程幾乎沒變,根據(jù)式(2)克努森數(shù)不斷減小,氣體滑脫效應(yīng)的影響逐漸減小(見圖3(c)).3.3.2氣竄通道與氣體滑脫效應(yīng)在水驅(qū)過程中,由于巖心具有非均質(zhì)性,形成水流優(yōu)勢通道(見圖4(a)).水驅(qū)后開始CO2驅(qū)替,由于水相的流度遠(yuǎn)大于油相的,水相更容易被氣相驅(qū)動,因此CO2優(yōu)先沿著水竄通道前進,迅速取代連續(xù)的水相并形成氣竄通道(見圖4(b)),從而產(chǎn)生氣體滑脫效應(yīng).由于新形成的氣竄通道半徑已經(jīng)很大,因此氣竄通道被拓寬而使氣體滑脫效應(yīng)減小的能力減弱;隨著氣相飽和度的增大,不斷形成新的氣竄通道,使氣體滑脫效應(yīng)逐漸增加.在CO2驅(qū)發(fā)生氣竄前,氣體滑脫效應(yīng)對氣相滲流的影響很小;隨著氣竄通道的形成,氣體滑脫效應(yīng)迅速增大,使氣竄的氣體流量急劇增大;在高氣相飽和度區(qū)域,滑脫效應(yīng)的影響趨于穩(wěn)定.3.3.3中央水相分離,不與油相接觸CO2驅(qū)替實驗和水驅(qū)后CO2驅(qū)替實驗?zāi)M非穩(wěn)態(tài)排驅(qū)過程的氣相分布.在水—CO2同注驅(qū)替時能夠模擬穩(wěn)態(tài)驅(qū)替過程.水被CO2氣體攜帶進入巖心后,CO2與水相可以分離并與油相接觸(見圖5(a)).氣相與水相進入巖心后,氣相段塞的兩側(cè)被水相封閉,形成交替的水氣相小段塞流(見圖5(b)).此時氣體不再是連續(xù)相,因此滑脫效應(yīng)被削弱;在驅(qū)替過程中產(chǎn)生賈敏效應(yīng),使驅(qū)替阻力增加,也影響滑脫效應(yīng)的測定.因此,用穩(wěn)態(tài)法測定氣體滑脫效應(yīng)與飽和度的關(guān)系與非穩(wěn)態(tài)法測定的結(jié)果存在差異.4氣體滑脫因子與相對滲透率利用油藏數(shù)值模擬研究CO2驅(qū),除了考慮CO2的相態(tài)性質(zhì)、CO2在油中的溶解性質(zhì)、溶解造成的油相黏度變化等,還要考慮氣體滑脫效應(yīng)的影響,當(dāng)修正相對滲透率時,一般將滑脫因子作為定值.由于在不同含水飽和度下的氣體滑脫效應(yīng)不同,在相對滲透率曲線測定過程中滑脫因子也不是定值,因此無論是實驗測定數(shù)據(jù)還是現(xiàn)場實測數(shù)據(jù),需要將氣體滑脫效應(yīng)隨氣相飽和度的變化規(guī)律引入相對滲透率的計算.在延長油田地層條件下,修正取心巖心的非穩(wěn)態(tài)排驅(qū)實驗數(shù)據(jù),得到考慮氣體滑脫因子與氣相飽和度關(guān)系的相對滲透率曲線(見圖6).由圖6可以看出,忽略滑脫因子隨氣相飽和度的變化,使用在CO2飽和度為100%條件下測定滑脫因子計算的油氣相對滲透率,其氣相相對滲透率偏大,并且曲線也更直.在氣相飽和度小于20%時,氣相處于束縛狀態(tài)或非連續(xù)相狀態(tài),幾乎不存在氣體滑脫效應(yīng);在驅(qū)替前期滑脫效應(yīng)隨氣相飽和度增加而增大;在氣相飽和度為30%左右時滑脫效應(yīng)達到最大;在驅(qū)替后期氣體竄流通道不斷擴大,由式(2)可知氣體流動通道半徑擴大使氣體滑脫因子減小,氣體滑脫效應(yīng)的影響逐漸被削弱.根據(jù)靖邊采油廠某試驗區(qū)長6層特低滲油藏的測井資料及地質(zhì)資料,同時考慮啟動壓力梯度、非均質(zhì)性和原油性質(zhì)變化等因素,利用Eclipse建立相控地質(zhì)模型及屬性模型;使用模擬蒸餾結(jié)果(見表2)建立組分模型,結(jié)合PVT實驗數(shù)據(jù),利用PVTi模擬油氣藏流體相態(tài)特性及溶解性質(zhì).將氣體滑脫因子與氣相飽和度的關(guān)系引入建立的模型,根據(jù)各網(wǎng)格的氣相飽和度修正CO2滑脫因子和流速,并引入修正的相對滲透率曲線分析每個網(wǎng)格動態(tài)和總體動態(tài).根據(jù)注采壓力、注入量等生產(chǎn)數(shù)據(jù)對CO2驅(qū)油效果進行模擬.在未通過歷史擬合修正相對滲透率時,與定滑脫因子的模擬結(jié)果比較,采用變滑脫因子進行數(shù)值模擬能更好地擬合生產(chǎn)氣油比、換油率等生產(chǎn)數(shù)據(jù),在一定程度上提高歷史擬合的效率,提高預(yù)測提高采收率的效果和氣竄通道形成的準(zhǔn)確性.5氣相飽和度的影響(1)針對CO2非混相驅(qū),研究氣體滑脫效應(yīng)與氣相飽和度的關(guān)系,測定高氣相飽和度條件下的滑脫因子,氣體滑脫效應(yīng)隨著氣相飽和度的降低而降低.