頁巖儲層納米級孔隙結(jié)構(gòu)及吸附氣解吸作用研究

頁巖儲層納米級孔隙結(jié)構(gòu)及吸附氣解吸作用研究

近年來，磺氣逐漸被認(rèn)為是一種非常有利的天然氣資源。目前，它已在中國西南部的廣大地區(qū)進(jìn)行了勘探和開發(fā)。頁巖氣藏儲層孔隙結(jié)構(gòu)較為復(fù)雜,孔隙非常小,使用核磁共振、CT掃描,電鏡掃描,高壓壓汞以及比表面儀等研究頁巖的孔隙結(jié)構(gòu),發(fā)現(xiàn)納米級孔隙所占比例很大。Howard發(fā)現(xiàn)Frio頁巖的孔隙半徑主要在5~15nm,Sondergeld等、MilnerM(2010)、ElgmatiM、CurtisME(2011)等人都發(fā)現(xiàn)了頁巖有機(jī)質(zhì)中的納米級孔隙。頁巖中納米級孔隙的存在使得氣體在這些孔隙中的流動方式及控制方程的研究非常重要。有20%~85%的頁巖氣是以吸附氣的狀態(tài)存在,開采后隨著儲層壓力降低氣體逐漸從吸附層中釋放出來并進(jìn)入到納米級孔隙中進(jìn)行擴(kuò)散滲流。頁巖本身兼具煤儲層、致密砂巖儲層和常規(guī)天然氣儲層的特點(diǎn)使其開采過程滲流非常復(fù)雜。JavadpourF(2009)提出了頁巖的表觀滲透率并且對表觀滲透率與達(dá)西滲透率的影響因素進(jìn)行了研究,但是這個(gè)滲透率模型的適用范圍比較窄,只適應(yīng)于氣體的滑脫流且缺乏吸附氣解吸對滲流的影響。筆者研究了頁巖納米級孔隙中氣體流動以及吸附氣解吸過程中滲透率的動態(tài)變化。1氣體擴(kuò)散效應(yīng)頁巖氣在開采過程中,隨儲層壓力的下降滲透率發(fā)生動態(tài)變化?？紫栋创笮】梢苑譃榇罂?直徑大于50nm)、小孔(在2~50nm之間)和微孔(小于2nm)。由于孔隙直徑達(dá)到納米級別,因此除受到吸附氣解吸效應(yīng)影響外還受到納米級孔隙氣體擴(kuò)散效應(yīng)影響。納米級孔隙氣體擴(kuò)散效應(yīng)指孔隙流動通道直徑很小,氣體分子平均自由程與孔隙直徑大小接近時(shí),氣體分子與孔隙壁面分子的碰撞概率大大增加,滲透率變差。吸附氣解吸效應(yīng)指儲層壓力下降到低于氣體臨界解吸壓力后,吸附態(tài)頁巖氣發(fā)生解吸導(dǎo)致頁巖基質(zhì)收縮變形,氣體滲流通道增加,滲透率變好;隨著儲層壓力的降低,頁巖儲層大量的吸附氣開始解吸,頁巖基質(zhì)收縮效應(yīng)逐漸加強(qiáng),使得氣體滲流通道逐漸變寬,滲透率不斷增加;當(dāng)儲層壓力降至更低水平時(shí),低壓條件下氣體擴(kuò)散效應(yīng)加劇,使得滲透率不斷降低。在開采過程中滲透率要受這兩種因素耦合作用影響,孔隙直徑越小耦合作用表現(xiàn)得越明顯。2巖石儲量變化的數(shù)學(xué)模型參數(shù)2.1多孔介質(zhì)滲透率的表征氣體在孔隙內(nèi)流動時(shí),氣體在孔隙內(nèi)滲流時(shí)發(fā)生的相互作用為:氣體分子間的碰撞、氣體分子與孔隙壁面分子的碰撞。兩種碰撞作用的物理機(jī)制不同,表現(xiàn)在滲流規(guī)律上也不同,分別為黏滯流和擴(kuò)散流。氣體分子的自由程與孔隙直徑相比小于1時(shí),主要發(fā)生氣體分子之間的相互碰撞;如果比值大于1,則主要產(chǎn)生氣體分子與孔隙壁面分子之間的碰撞。因此將氣體分子自由程大于孔隙直徑(D)的分子所占總的分子量的比例為α,那么小于D的則占1-α。孔隙內(nèi)符合達(dá)西流動產(chǎn)生的流量為:孔隙內(nèi)由分子擴(kuò)散引起的滑脫流動流量為:故孔隙內(nèi)由兩種流動機(jī)制產(chǎn)生的氣體總流量為:對于真實(shí)氣體,有將式(5)、(6)代入式(4)得:式(7)等式右邊項(xiàng)中ue014(p/Z)可展開成:氣體壓縮系數(shù)(cg)定義如下:將式(8)、(9)代入式(7)后化簡得:根據(jù)達(dá)西公式得到考慮氣體擴(kuò)散后的滲透率的表達(dá)式為:式(11)中K∞與多孔介質(zhì)的結(jié)構(gòu)、孔隙幾何形態(tài)等有關(guān),是多孔介質(zhì)的滲透率。根據(jù)假設(shè)的毛細(xì)管模型的絕對滲透率,通過對泊稷葉定律推導(dǎo)后得到:2.2氣體分子平均自由程達(dá)西流動和分子擴(kuò)散流動是氣體在孔隙通道內(nèi)流動的主要機(jī)制。孔隙直徑不同則兩種流動機(jī)制所發(fā)揮的作用不同。根據(jù)分子運(yùn)動理論,自由程描述了氣體分子在未與其他分子發(fā)生碰撞前經(jīng)過的路程,氣體分子平均自由程的表達(dá)式為:若考慮孔隙直徑為D,假設(shè)氣體分子自由程大于D所占總的氣體量的比例為即分配系數(shù),那么小于D的則占1-α。即那么小于D的分子數(shù)則為1-α=1-e-D/λ。當(dāng)多孔介質(zhì)越致密,或氣體壓力低、氣體分子平均自由程大時(shí),自由程大于孔隙直徑的分子與巖壁碰撞對總流量的貢獻(xiàn)將隨之增大,擴(kuò)散現(xiàn)象愈顯著。2.3儲層變形對儲層滲透率的影響根據(jù)Curtis對Barnett、Antrim等統(tǒng)計(jì),頁巖中含有大量的吸附氣體,占據(jù)到20%~85%。由于未能考慮吸附氣解吸在頁巖開發(fā)中的重要性,因此對頁巖氣的開采預(yù)測產(chǎn)生嚴(yán)重的偏差。隨著儲層壓力降低,頁巖中吸附氣體開始解吸,頁巖基質(zhì)收縮改變滲流通道對滲透率有重要影響。引用Bangham固體變形理論分析壓力下降吸附氣解吸對頁巖氣解吸滲透率影響。儲層巖體形變程度與儲層壓力的關(guān)系式為:假設(shè)儲層吸附氣體為一元?dú)?即甲烷氣),氣體吸附及解吸附Langmuir方程為:將式(16)帶入式(15)中,積分后得到頁巖基質(zhì)收縮程度為:隨著儲層壓力的降低,吸附氣體開始解吸,在表面張力的作用下頁巖開始收縮,同時(shí)裂隙內(nèi)的有效應(yīng)力增加,巖體也產(chǎn)生膨脹變形,則總變形量為:Seidle等模型中研究了煤層中氣體解吸收縮導(dǎo)致裂隙張開,孔隙度變大,得出了孔隙度和儲層形變間的關(guān)系。即將式(18)代入式(20)得到儲層形變與基質(zhì)孔隙度間的關(guān)系:假設(shè)在頁巖儲層中孔隙體積的縮小帶來流體流動通道的成比例變化。根據(jù)毛細(xì)管模型,得到假想巖石孔隙度和孔道半徑間關(guān)系式:因此得到了考慮微觀孔隙氣體擴(kuò)散與吸附氣解吸的頁巖氣基質(zhì)滲透率動態(tài)數(shù)學(xué)模型:3儲層壓力和孔隙直徑對氣體分子擴(kuò)散流分配系數(shù)影響由于頁巖氣主要成分是甲烷,因此模型計(jì)算中所選氣體為甲烷,選取頁巖儲層物性參數(shù)數(shù)據(jù),計(jì)算了頁巖基質(zhì)在不同的初始孔隙直徑下隨著儲層壓力降低其氣體的滲透率動態(tài)。其中的關(guān)鍵參數(shù)是:E=2.68×104MPa,V0=22.4×10-3m3/mol,T=323K,R=0.008314MPa·m3·K/kmol,μ=0.015mPa·s,φi=0.08,b=0.22MPa-1,Vm=18.66m3/t,ρ巖=2.56t/m3,cp=4.35×10-4MPa-1。根據(jù)此數(shù)學(xué)模型及相關(guān)頁巖參數(shù),計(jì)算得到不同壓力和孔隙直徑下甲烷分子進(jìn)行達(dá)西流動和擴(kuò)散流動的分配系數(shù)變化曲線和孔隙的滲透率變化曲線。由圖1可知,在儲層壓力較高時(shí),甲烷氣體分子擴(kuò)散流動分配系數(shù)很低,氣體分子以達(dá)西流動為主,隨儲層壓力下降擴(kuò)散流動分配系數(shù)指數(shù)式增加,孔隙通道越粗?jǐn)U散流分配系數(shù)壓力拐點(diǎn)越低。從圖2可知,孔隙直徑對滲透率變化影響比較大,孔隙直徑10nm時(shí),隨著壓力降低,孔隙通道受基質(zhì)收縮影響不斷增大而滲透率增加,壓力下降后期擴(kuò)散流動作用影響明顯增大,使得滲透率開始減小;滲透率變化存在拐點(diǎn),當(dāng)孔隙直徑從3nm增到10nm時(shí)拐點(diǎn)值從9MPa降到3MPa;孔隙直徑越大,滲透率變化拐點(diǎn)出現(xiàn)得越遲且拐點(diǎn)出現(xiàn)后滲透率下降越快;當(dāng)孔隙直徑較大時(shí),氣體分子的擴(kuò)散流動已對滲透率影響很小,達(dá)西流動占據(jù)主要位置,同時(shí)氣體解吸引起孔隙通道增加滲透率增加;儲層在壓力較低(小于10MPa)的情況下,氣體滲流受分子擴(kuò)散效應(yīng)影響較強(qiáng),而在儲層壓力較高的情況下擴(kuò)散效應(yīng)不明顯,與相關(guān)實(shí)驗(yàn)結(jié)果相同。4儲層滲透率模型1)頁巖氣儲層孔隙結(jié)構(gòu)非常復(fù)雜,孔隙非常小,使用各種測試方法研究頁巖的孔隙結(jié)構(gòu),發(fā)現(xiàn)納米級孔隙比例很大;同時(shí)氣藏吸附態(tài)氣體占據(jù)很大比例。2)根據(jù)頁巖開采滲流特征,建立了頁巖氣在基質(zhì)微觀孔隙滲流過程中滲透率動態(tài)數(shù)學(xué)模型,該模型全面的考慮了頁巖吸附氣解吸和納米孔隙內(nèi)分子擴(kuò)散流動對滲透率的影響。3)模型計(jì)算得出,在孔隙直徑較小條件下滲透率受擴(kuò)散流與解吸雙因素影響隨儲層壓力下降呈現(xiàn)出先增加后減小的趨勢;對于大孔道,氣體流動以達(dá)西流動為主,受氣體解吸影響不斷增加。4)儲層在壓力較低(小于10MPa)的情況下,氣體滲流受分子擴(kuò)散效應(yīng)影響較強(qiáng),而在儲層壓力較高的情況下擴(kuò)散效應(yīng)不明顯,對于微觀孔隙直徑,孔隙直徑越大,滲透率下降拐點(diǎn)壓力越低,且下降速度越快。儲層參數(shù)3q1為達(dá)西流動產(chǎn)生的流量,m3/s;q2為分子擴(kuò)散引起的滑脫流動流量,m3/s;A為流動截面積,m2;K∞為滲透率,D;l為流動長度,m;p為壓力,MPa;ρ為氣體密度,g/m3;D為孔隙直徑,nm;α為擴(kuò)散流分配系數(shù),無量綱;M為摩爾質(zhì)量,g/mol;μ為黏度,mPa·s;Dk為氣體擴(kuò)散系數(shù),m2/s;C為濃度,kmol/m3;Z為氣體壓縮因子,無量綱;cg為氣體壓縮系數(shù),MPa-1;K°為玻爾茲曼氣體常數(shù),1.3806505×10-23J/K;d為分子直徑,nm;T為絕對溫度,K;λ為氣體分子平均自由程,m;φ為孔隙度,%;Δε為頁巖收縮程度,無量綱;Vm為氣體的Langmuir體積,m3/t