頁巖氣滑脫、擴散傳輸機理耦合新方法?

李亞雄1)2)3)? 劉先貴2)3) 胡志明2)3) 高樹生2)3) 端祥剛2)3) 常進2)3)

1)(中國科學院大學,北京100049)

2)(中國科學院滲流流體力學研究所,廊坊 065007)

3)(中國石油勘探開發研究院廊坊分院,廊坊065007)

頁巖氣滑脫、擴散傳輸機理耦合新方法?

李亞雄1)2)3)? 劉先貴2)3) 胡志明2)3) 高樹生2)3) 端祥剛2)3) 常進2)3)

1)(中國科學院大學,北京100049)

2)(中國科學院滲流流體力學研究所,廊坊 065007)

3)(中國石油勘探開發研究院廊坊分院,廊坊065007)

(2016年10月17日收到;2017年3月31日收到修改稿)

針對頁巖氣流動計算中所用耦合機理不同的現狀,且為了厘清滑脫和各種擴散之間的關系,首先采用理論分析和數學模型的方法,根據定義和微觀運動機制對滑脫和各種擴散進行了分析,然后在考慮吸附層頁巖氣分子所占空間對氣體流動影響的情況下,提出了“壁聯擴散”的概念來表征克努森擴散和表面擴散的總效應,并指出壁聯擴散和滑脫效應等同,由此提出了壁聯擴散和滑脫效應在流動計算中可互換而不重復疊加的耦合新方法.實例驗證表明,當毛細管半徑從5 nm增大到2000 nm,壁聯擴散和滑脫效應的質量通量相對誤差較小,在絕大部分范圍內都小于10%,且在整個孔徑范圍內兩者平均值相差1.4×10?6kg·m?2·s?1,即平均值的相對誤差僅為5.8%,該方法可以滿足工程計算的需要.考慮到參數選取、機理數學模型有待完善等方面的影響,新方法的論證存在進一步提升的空間.壁聯擴散的提出具有實際開發意義和多重研究意義,耦合新方法的提出闡明了滑脫和各種擴散之間的關系,防止了頁巖納米級孔隙中流動機理的重復疊加,能較好改變頁巖氣流動計算耦合方法不一致的現狀,為頁巖氣開發定量計算指明新方向.

頁巖氣,擴散,壁聯擴散,滑脫

1 引言

頁巖氣是指主要以游離或吸附狀態聚集在富含有機質的暗色泥頁巖或高碳泥頁巖中的自生自儲天然氣[1,2].全球頁巖氣資源量為456.24×1012m3[3],其中我國的頁巖氣資源非常豐富,大約為26×1012m3[3].頁巖氣藏具有超低滲透率,其滲透率的范圍為10?9—10?3μm2[4,5],對頁巖氣藏巖性特性分析表明,主要的納米孔徑范圍為5—200 nm[4,6],因此具有低孔滲性.常規油氣藏中,可壓縮流體在孔隙介質中常用的流動方程適用于幾十至幾百微米的孔隙介質,但是對于頁巖氣藏,其納米孔隙發育,Darcy滲透率不再適用[7].且部分氣體吸附于儲層巖石顆粒表面以吸附態存在于儲層中[8],在納米和分子量級,存在以下幾種氣體形式:孔隙中的自由氣、孔壁的吸附氣和在干酪根中的溶解氣[9,10],如圖1[11]所示.頁巖氣復雜的賦存形式更加增大了其研究難度.事實上,水力壓裂裂縫的滲透率比頁巖基質高幾個數量級,基質中的氣體流動才是長期生產速度的制約因素[9].因此,為了更加高效地進行頁巖氣的開發,需要重點對頁巖基質中的流動規律進行研究.不同于常規氣藏,除了達西定律,頁巖基質納米級孔隙中還存在滑脫、菲克擴散、克努森擴散等流動機制,對于這幾種機理如何進行耦合還沒有形成統一的觀點.

圖1 頁巖氣分子納米孔賦存示意圖[11]Fig.1.Diagrammatic sketch of shale gas molecules in a nanopore[11].

目前的研究中,頁巖氣各種傳輸機理的耦合方法各有差異.Javadpour[11]提出了考慮克努森擴散和滑脫雙重作用的表觀滲透率計算公式.Shabro等[12]建立了同時考慮干酪根中溶解氣擴散、滑脫、克努森擴散和Langmuir解吸的預測頁巖氣產量的算法.Wu等[13]、吳克柳等[14]、Haghshenas等[15]也考慮了克努森擴散和滑脫的耦合.宋洪慶等[16]在納米孔隙中考慮了達西滲流和克努森擴散,將二者線性疊加得到納米孔中的流動通量.Mi等[17]、糜利棟等[18]在不同克努森數Kn的條件下分別考慮了不同的擴散形式,有克努森擴散、過渡擴散和菲克擴散三種形式,采用氣體擴散系數進行數學表征,而氣體滑脫效應采用氣體滑脫修正系數進行數學表征,并將擴散和滑脫進行線性耦合得出基質滲透率的數學表達式.艾爽等[19]在頁巖基質的數學模型中,將達西速度和滑脫速度線性疊加,并引入表觀滲透率對此進行描述.Li等[20]將連續流、滑移流、過渡流和自由分子流進行線性耦合.Jia等[21]在基質流動中考慮了吸附解吸、達西流和菲克擴散作用,其中將達西流和菲克擴散線性疊加得出表觀滲透率的表達式.Yin等[3]建立了雙重介質數學模型,除了考慮等溫吸附解吸、滑脫效應和菲克擴散,還考慮了啟動壓力梯度,建立基質和裂縫的質量守恒方程,用有限差分解出的產量和某一生產實際相比相對誤差小于5%.

由此可見,目前的研究中頁巖氣在納米孔中的傳輸規律具有不同的耦合形式,不同的研究考慮的流動機理并不相同.其中比較突出的方面是頁巖氣的各種流動機理,如滑脫、菲克擴散、克努森擴散等之間存在怎樣的關系,以及在具體的流動計算中上述各種機理是否存在重復疊加的問題和如何處理它們之間的關系等,關系到頁巖基質中流動能力及可動性的正確估算,從而成為頁巖氣開發過程中的關鍵之一.對于以上問題,目前的文獻還沒有給出答案.因此,本文首先根據微觀運動機制和數學模型對滑脫和各種擴散進行分析,并推導出有效半徑的數學表達式,然后提出“壁聯擴散”的概念和滑脫、擴散傳輸機理耦合新方法,最后用實際頁巖氣藏的數據及相關數學模型對上述方法進行分析.

2 理論基礎與新方法的提出

2.1 滑脫傳輸表征

文獻[22,23]提出根據克努森數Kn的取值范圍可以把流體的流態分為四類:Kn<10?3為連續流態;10?310為自由分子流態.頁巖氣納米孔氣體傳輸機理包括納米孔體相氣體傳輸和吸附氣表面擴散,而納米孔體相氣體傳輸機理包括連續流動、滑脫流動、過渡流和克努森擴散[14,22,24].

當用氣體測量巖石的滲透率時,測量結果不僅比液測值高,而且還出現了較強的壓力依賴性.1941年Klinkenberg[25]最早發現了該現象,并把該現象歸因于氣體在巖石孔隙中的滑脫行為所致,此即所謂的滑脫效應或Klinkenberg效應[26].所以“滑脫效應”使得巖心樣品的氣測滲透率大于液測滲透率.后來很多學者通過實驗對Klinkenberg滑脫因子進行修正,或者引入克努森參數校正[22].根據錢學森對不同Kn數的氣體流動所做的特性劃分[27],朱光亞等[28]認為Klinkenberg方程的理論基礎僅適用于Kn=0.01—0.1的流動范圍,在這一流動區域中,控制流動的方程仍然是經典流體動力學的Navier-Stokes(N-S)方程,在物面上必須應用速度滑移和非等溫條件下溫度跳躍的邊界條件.而當Kn>0.1時,在這一流動區域中,流體分子平均自由程與連續介質的特征長度屬同一量階,氣體分子間的碰撞和氣體分子與物面間的碰撞對氣體運動的影響具有同等重要的意義,連續介質假設[29]已不再成立.

從具體含義上來看,氣體滑脫指氣體在介質孔道流動中出現的近孔道壁面氣體分子對壁面發生相對運動的現象[30].氣體滑脫流動的本質是由于氣體分子與孔道固壁的作用使得氣體在孔道固壁附近的各個氣體分子都處于運動狀態,且貢獻一個附加通量,從而在宏觀上表現為氣體在孔道壁面上具有非零速度,產生滑脫流量[28,31].氣體滑脫效應可以表示為圖2,一部分分子以一個大于零的速度沿著壁面運動,另一部分頁巖氣分子與孔壁發生碰撞后,以某一速度反彈,在壁面附近速度方向與壁面相切或與壁面成某一非零角度,這兩部分分子共同構成了壁面附近由于滑脫作用而處于運動狀態的氣體分子.

圖2 滑脫效應示意圖Fig.2.Diagrammatic sketch of slippage e ff ect.

圓管中理想氣體層流的質量通量為[11]

式中Ja為圓管中理想氣體層流的質量通量(kg·m?2·s?1);r為流管半徑(m);ρavg為平均密度(kg·m?3);μ為氣體黏度(Pa·s);p為壓力(Pa).

滑脫流動使得孔隙中的天然氣滲流更易進行.Brown等[32]引入一個F因子,校正流管中滑移速度,即

式中F為滑移速度校正因子;R為通用氣體常數(J·mol?1·K?1);T為絕對溫度(K);M為分子量(kg·mol?1);pave為孔隙進出口壓力的平均值(Pa);α為角動量調節系數,頁巖中需要相應的測試實驗來確定.

所以滑脫流態下的總質量通量為

式中Jt為滑脫流態下的總質量通量(kg·m?2·s?1).

因此,相較于達西連續流動,由于滑脫效應而導致的流通量為

式中Jh為因滑脫效應而導致的流通量(kg·m?2·s?1).

文獻[10]認為多孔介質由一定百分比的開放孔隙(孔隙度)和導致流通路徑長于直線路徑的互連網絡(迂曲度)所組成,則Jh的修正形式如下:

式中Jhm為Jh的修正形式(kg·m?2·s?1);?為頁巖孔隙度;τ為迂曲度;?/τ為孔隙度-迂曲度因子.

2.2 擴散傳輸表征

2.2.1 克努森擴散

當Kn>10時,分子的平均自由程大于孔隙直徑,此時頁巖氣分子和孔隙壁之間的碰撞占主導地位,而分子之間的碰撞退居次要地位,此即克努森擴散[17,18,33].克努森擴散示意圖如圖3.

圖3 克努森擴散示意圖Fig.3.Diagrammatic sketch of Knudsen di ff usion.

納米孔隙中克努森擴散可以表達成壓力梯度的函數,忽略黏性作用,頁巖氣質量通量為[10]

式中JD為克努森擴散質量通量(kg·m?2·s?1);DK為克努森擴散系數(m2·s?1).

克努森擴散系數DK定義式為[10,34?40]

同理,用?/τ對克努森擴散通量進行修正:

式中JDm為JD的修正形式(kg·m?2·s?1).

2.2.2 菲克擴散

頁巖氣分子在濃度差的作用下,游離相頁巖氣從高濃度區向低濃度區運動,稱為菲克擴散[41].當Kn 6 0.1時,孔隙直徑遠大于甲烷氣體分子的平均自由程,這時甲烷氣體分子的碰撞主要發生在自由甲烷氣體分子之間,而分子和毛細管壁的碰撞機會相對較少,此類擴散遵循菲克定理,稱為菲克型擴散[17,18,33].菲克擴散示意圖如圖4.

圖4 菲克擴散示意圖Fig.4.Diagrammatic sketch of Fick’s di ff usion.

菲克第一定律表示每個時刻基質塊中各處氣體濃度相等,氣體擴散速度與內外濃度差呈線性關系.擴散方程為[1,42]

式中V為氣體濃度(m3·m?3);t為時間(s);Fs為基質塊的形狀因子(m?2);D為菲克擴散系數(m2·s?1);VE為裂縫壁面氣體濃度(m3·m?3).

菲克第二定律表示基質塊中氣體濃度分布為時間和位置的函數,擴散方程為[42]

式中?為基質孔隙度.

2.2.3 過渡擴散

當0.1

Mi等[17]、糜利棟等[18]、陳強等[43]由菲克擴散和克努森擴散擴散系數的調和平均得到過渡擴散的擴散系數.而文獻[40]運用通過實驗研究得出的權重系數將黏性流和自由分子流進行耦合得出過渡流的質量流量.關于過渡流的模型正處于研究之中,目前還沒有一種被廣泛采用的過渡流表達式.

2.2.4 表面擴散

吳克柳等[44]認為不論是氣-固吸附系統,還是液-固吸附系統,表面擴散對多孔介質中流體的質量傳輸均具有重要作用.盛茂等[45]認為由于活性炭納米管與干酪根微-納米孔隙在物質組成(碳質)和結構(微-納米孔)方面類似,且氣體在活性炭納米管中的流動存在表面擴散的現象已被實驗證實[46],由此可推測表面擴散可能是頁巖氣傳質的重要方式之一,并利用Langmuir滑脫邊界條件推導了考慮頁巖吸附層對氣體滑脫影響的滑脫流動質量通量計算式.眾多學者研究表明,頁巖吸附氣量用單層吸附的Langmuir等溫式計算更加合理[1,47?50].因此,跳躍模型更加適合頁巖納米孔吸附氣表面擴散的研究[44].跳躍模型假設吸附氣分子是從孔隙表面一吸附位跳躍到相鄰吸附位,這也視為吸附氣分子的活化過程.如果吸附氣分子獲得足夠的能量,并越過兩個相鄰吸附位之間的能量障礙,則活化過程發生,吸附氣發生表面擴散.表面擴散是擴散粒子在吸附位之間隨機跳躍的連續過程,其中每一個跳躍都需要最小的活化能并經歷活化過渡態,活化能與擴散粒子和孔隙表面的吸附能有關,常呈正比關系[44].圖5為表面擴散的示意圖.

圖5 表面擴散示意圖Fig.5.Diagrammatic sketch of surface di ff usion.

依據文獻[45],表面擴散質量通量可寫為

其中Js為表面擴散質量通量(kg·m?2·s?1);Ds為零負荷下的表面擴散系數(m2·s?1);Cμ為頁巖吸附質濃度(kg·m?3).

Sladek等[51]根據30組氣-固系統的數據擬合出的表面擴散系數公式為

?H為等量吸附熱(J·mol?1);m為根據吸附鍵的種類而定的常數,由文獻[51]對吸附鍵的分類,選擇吸附質為非極性、固體吸附劑為導體的組合,所以此處m=1.

頁巖吸附質濃度可用Langmuir吸附方程表示如下:

其中ρs為頁巖基質密度(kg·m?3);Vstd為標準狀態下的氣體摩爾體積(m3·mol?1);qL為Langmuir體積(m3·kg?1);pL為Langmuir壓力(Pa).

所以,

2.3 滑脫、擴散機理耦合新方法

在頁巖納米級小孔隙中,當考慮吸附層頁巖氣分子的大小對氣體流動的影響,由Langmuir等溫吸附方程可得

其中Va為吸附體積(m3);VL為Langmuir體積(m3).

假設頁巖氣全部由甲烷組成,對于平行毛管束所組成的多孔介質,流動發生在圓柱狀毛細管中,Langmuir等溫吸附理論為單分子層吸附,當吸附體積達到最大值時,假定此時相鄰分子之間的間隙可以忽略.而當某一壓力下氣體分子吸附在頁巖壁面時,設此時氣體分子的存在對實際流動區域的影響和與該壓力下吸附氣體總體積相同的一層等厚無間隙且鋪滿孔壁的等效吸附層相同,則得

其中ri為第i個毛管的半徑(m);dm為甲烷的分子直徑(m);dme為壓力為p時的等效吸附層厚度(m);L為平行毛管束的長度(m);n為毛管數.

由(15),(16)和(17)式得

故孔隙的有效半徑為

其中re為有效半徑(m).

將(20)式代入(5)和(8)式,得到考慮吸附層厚度影響的滑脫效應和克努森擴散通量表達式:

克努森擴散和表面擴散都與氣-固相之間在界面處的作用有關,其中克努森擴散是氣體分子和壁面產生碰撞以后所發生的擴散作用,表面擴散是氣體分子在孔隙表面吸附位之間連續跳躍的過程,兩者的作用都是導致氣體分子在壁面周圍的運移速度不為零,這與滑脫現象的本質是一致的,如果提出新名詞——“壁聯擴散”來表征表面擴散和克努森擴散的綜合效應,則可從以上微觀運動機制的分析得出壁聯擴散與滑脫效應等同的結論,并提出一種滑脫和擴散耦合方法的新認識,即在頁巖氣納米孔中壁聯擴散和滑脫效應對總流動的貢獻相同,亦即在流動計算中兩種機制互不疊加且可相互替換,如圖6所示.所以,在目前研究中出現的將克努森擴散和滑脫進行疊加耦合,或將滑脫和表面擴散進行疊加,抑或將三者都考慮在傳輸機制中的計算方法,都將導致重復疊加的問題.出現這種現象的原因歸根究底源于人們對滑脫和幾種擴散之間相互聯系認識的缺失,從而因為各種機理名稱的不同各持己見所致.壁聯擴散的概念以及耦合新方法的提出能夠厘清滑脫和各種擴散之間的關系,在頁巖氣計算中統一所需考慮的流動機制,避免納米級孔隙中流動機理的重復疊加,從而較好地改變頁巖氣流動耦合方法不一致的現狀,為頁巖氣開發工作中的定量計算指明新方向.

因此由上述及(14),(21)和(22)式可得:

圖6 壁聯擴散和滑脫關系圖Fig.6.Relationship between wall-associated di ff usion and slippage e ff ect.

其中Jb為壁聯擴散質量通量(kg·m?2·s?1).

3 耦合方法實例分析

依據文獻[45,52—55],參考典型頁巖氣藏Barnett頁巖參數作為基礎數據:毛細管半徑為5×10?9—2×10?6m;原始地層壓力為2×107Pa;頁巖孔隙度為5%;甲烷氣的動力黏度為1.8×10?5Pa·s;地層溫度為366.5 K;頁巖密度為2500 kg·m?3;Langmuir體積0.0027 m3·kg?1;Langmuir壓力為2.46×106Pa;頁巖等量吸附熱為8000 J·mol?1;τ=1;α=0.8;M=0.016 kg·mol?1;dm=0.4×10?9m;L=50 m;R=8.3144621 J·mol?1·K?1;由溫壓條件根據REFPROP軟件得ρavg=84.176 kg·m?3.分別得出質量通量隨孔徑的變化曲線、滑脫和壁聯擴散質量通量平均值對比圖和質量通量相對誤差隨孔徑的變化曲線如圖7—圖9所示.

3.1結果與討論

從圖7可知,在毛管半徑從5 nm到2000 nm逐漸增大的過程中,壁聯擴散和滑脫效應的變化范圍分別為1.39×10?7—4.88×10?5kg·m?2·s?1和1.07×10?7—4.60×10?5kg·m?2·s?1,而相應條件下滑脫效應和壁聯擴散的擴散通量相差3.21×10?8—2.77×10?6kg·m?2·s?1.從圖9可見,相對誤差隨著毛管半徑的增大而逐漸減小.半徑為2000 nm時,最小相對誤差為5.7%.半徑為50 nm處,相對誤差曲線存在拐點,即拐點之前曲線迅速下降,拐點之后曲線降低速度急速減緩,相對誤差趨于平穩,其中拐點處的相對誤差為7.6%.在半徑為18—2000 nm范圍內,相對誤差都在10%以下.當半徑為5 nm時,最大相對誤差達到23.1%.半徑小于10 nm的范圍內,相對誤差明顯大于其他范圍的數值,原因可能是微尺度流動中,由于不同壁面力場的疊加[56],甲烷在微孔中的吸附不是簡單的單層吸附,而是以填充方式進行[57],孔隙微尺度效應顯著增強,導致發生表面擴散和克努森擴散的氣體不再受到通常情況下單一壁面作用力的影響,從而滑脫和壁聯擴散公式對氣體流動的描述與實際產生更大偏差,由于兩公式從不同的角度建立,相同微尺度條件下產生的偏差也會不同,這很可能是導致此時滑脫效應和壁聯擴散的數學模型一致性變差的主要原因.分析圖中的曲線趨勢可知,在半徑小于5 nm的孔道中,滑脫效應和壁聯擴散的相對誤差會進一步加大,這也與微尺度效應的進一步增強有關.但同時應看到,由于此時力場疊加作用使數學模型對實際流動的描述能力也在下降,因此耦合方法在小于5 nm孔隙中的適用性還不能做出定論.由于在上述所選取參數的計算中,表面擴散通量與孔徑無關,因此隨著孔徑逐漸變大,克努森擴散通量的增量大于滑脫效應通量的增量,這說明在所研究的條件下,相較于滑脫效應,克努森擴散對孔徑的變化更加敏感,所以在研究孔徑變化對流量的影響時,應更加注意克努森擴散的相應變化對流量的影響.在質量通量隨著孔徑變化的過程中,圖7中的曲線為近似過原點的直線,因此在所研究的大部分孔徑范圍內,頁巖氣吸附層分子所占的空間對于流動的影響較小,只有在孔徑與頁巖氣分子直徑量級相當時,吸附層才會對頁巖氣的流動傳輸產生顯著的影響.這說明在實際頁巖氣藏的計算中,需要評估頁巖孔徑分布情況,如果某一開發區域頁巖在幾納米量級或更小范圍內的孔徑普遍存在且體積百分數可觀,則應該使用有效半徑的計算方法,使計算結果更加可信,反之則無需考慮有效半徑的概念.圖7中的兩條曲線整體上比較接近,從圖8和圖9可以看出,隨著孔徑變化,滑脫和壁聯擴散的相對誤差總體上維持在很低的水平,同樣在整個孔徑范圍內所得到的滑脫效應和壁聯擴散通量的平均值也很接近,僅相差1.4×10?6kg·m?2·s?1,即平均值的相對誤差只有5.8%.

圖7 質量通量隨孔徑的變化曲線Fig.7.Graph of mass fl ux vs aperture.

本文基礎數據來源于實際頁巖氣藏,上述流通量的結果證明,在頁巖微孔的耦合流動計算中,將滑脫效應和壁聯擴散視為具有同等效應的流動機理可以滿足工程計算的需要.而兩者在所研究的孔徑范圍內仍然有一定的相對誤差,分析有如下幾點原因:1)目前頁巖氣流動計算相關參數的缺乏,例如對于角動量調節系數α,目前未見專門針對頁巖的角動量調節系數測試實驗,α與頁巖氣在流管管壁上擴散反射與鏡面反射有關,其變化范圍為0—1,只能借助已有類似實驗結果,例如N2,Ar,CO2在硅微管中的結果[55]來類比甲烷與實際頁巖儲層的作用,這顯然會導致結果的偏差,因此今后隨著頁巖氣角動量調節系數測試實驗或分子模擬等研究的開展,可以為結論提供更可靠的依據;目前關于頁巖氣表面擴散系數詳細數據的報道也很少見,因此本文借用Sladek等[51]根據30組氣-固系統的數據擬合出的表面擴散系數公式套用到頁巖氣流動計算中,這也會導致驗證結果的偏離;所以從目前頁巖氣微尺度流動相關領域的研究進展來看,很多基礎數據仍然匱乏,關于頁巖油氣擴散的分子模擬研究工作正在持續開展[58,59],相關實驗或模擬的跟進可以更好地促進對頁巖氣流動機理的認識;2)數學模型有待改進,氣體相對于液、固等狀態的物質而言,最大的特點是具有壓縮性,溫壓等因素的變化對其影響不可忽略,頁巖氣高溫高壓的賦存環境與常見氣體的溫壓狀態有很大區別,而目前頁巖氣開發工作者所運用的計算公式基本上都不是針對頁巖氣的特點推導而來的,例如相關理論的建立是用來解釋極低壓下稀薄氣體的,相關模型也是在亞大氣壓(接近真空)條件下進行驗證或發展而來的[60],故在使用上必將產生偏差,因此有必要建立針對頁巖氣流動特點的獨立的一套數學模型體系.可以預見隨著頁巖氣開發研究的進一步深入,利用更加嚴謹的參數來驗證滑脫和各種擴散耦合新方法會產生更加滿意的結果.

圖8 滑脫和壁聯擴散質量通量平均值對比Fig.8.Average values of mass fl ux of slippage and wall-associated di ff usion.

圖9 質量通量相對誤差隨孔徑的變化Fig.9.Graph of relative error of mass fl ux vs aperture.

綜上理論和實例分析都表明,壁聯擴散和滑脫效應在流動計算中可互換而不重復疊加的耦合新方法是合理的,在頁巖氣藏流動分析中,視滑脫效應和壁聯擴散在數量上等同可以滿足工程計算的需要.

從上述結果可以看到,壁聯擴散描述了頁巖氣分子與壁面作用所產生的擴散機制,它成為厘清滑脫和各種擴散關系的橋梁,從而糾正和改進了目前頁巖基質中氣體流動能力的估算方法,這是它的實際開發意義所在.從研究意義上來看,壁聯擴散是滑脫的微觀表述形式,它從微觀運動機制上將近孔道壁面氣體分子對壁面發生相對運動的現象分為兩大部分,這兩部分的受力機制和運動形態有顯著差異,因此壁聯擴散不但具有形態描述功能,還為從不同的力學機制來研究滑脫現象提供了可能.另一方面的研究意義在于壁聯擴散突破了滑脫流在高克努森數下不再適用的局限性,因此在極小的納米孔隙中,當滑移流態不再適用的情況下,從物理形態得出的壁聯擴散能夠為探索更高克努森數范圍下的壁面周圍氣體流動形式提供思路,即當N-S方程失效而Boltzmann方程和Burnett方程的高階修正求解困難時,可以從壁聯擴散入手,嘗試運用分子模擬等方法探究極小納米孔隙中的所謂“滑移現象”.

4 結論

本文針對頁巖氣開發工作中的難點——基質流動中滑脫、擴散傳輸機理耦合的問題進行了研究.提出了有效半徑的概念,考慮了吸附層頁巖氣分子直徑對氣體流動的影響.結果表明,在實際頁巖氣藏的計算中,只有針對儲層巖石孔徑普遍與氣體分子直徑量級相當時,才需要運用有效孔徑的概念,而在其他的較大孔徑范圍內用孔徑代入公式計算即可,具體操作的選取需要視頁巖孔徑分布情況而定.結果還顯示,隨著孔徑的變化,在所研究的條件下克努森擴散對孔徑的變化更加敏感,因此當孔徑發生變化時,應更加注意克努森擴散的相應變化對流量的影響.而在頁巖氣各種傳輸機理的耦合方法上,還沒有形成統一的觀點,因此需要提出新觀點來改善現狀.故由微觀運動機制提出了壁聯擴散的概念,來表征表面擴散和克努森擴散效應的總和,并提出壁聯擴散和滑脫效應在流動計算中可互換而不重復疊加的耦合新方法.通過用典型頁巖氣藏的基礎數據進行驗證,結果顯示將滑脫效應和壁聯擴散在數量上視為等同可以滿足工程計算的需要.壁聯擴散的提出具有實際開發意義和多重研究意義,新方法摒棄了過去由于流動機理混淆而造成的頁巖納米級孔隙中流動機理的重復疊加,厘清了滑脫和各種擴散之間的關系,能較好地改變基質中流動計算耦合方法不一致的現狀,因此避免了頁巖基質中氣體流動能力的不當估算,由于頁巖基質中的氣體流動是長期生產速度的制約因素,耦合方法的提出對頁巖氣實際開發具有重要的指導意義.

[1]Zhang L H,Guo J J,Tang H M 2014 Development Foundation of Shale Gas Reservoir(1st Ed.)(Beijing:Petroleum Industry Press)p1,pp61–63,p84(in Chinese)[張烈輝,郭晶晶,唐洪明2014頁巖氣藏開發基礎(第一版)(北京:石油工業出版社)第1頁,第61—63頁,第84頁]

[2]Deng J,Zhu W Y,Liu J X,Zhang Z,Ma Q,Zhang M,Deng K,Ma L 2013 Nat.Gas Geosci.24 456(in Chinese)[鄧佳,朱維耀,劉錦霞,張貞,馬千,張萌,鄧凱,馬麗2013天然氣地球科學24 456]

[3]Yin D Y,Wang D Q,Zhang C L,Duan Y J 2015 Open Pet.Eng.J.8 203

[4]Deng J,Zhu W,Ma Q 2014 Fuel 124 232

[5]Zhu W Y,Qi Q 2016 Sci.Sin.Technol.46 111(in Chinese)[朱維耀,亓倩2016中國科學:技術科學46 111]

[6]Lin M,Jiang W B,Li Y,Yi Z X,Zhang Z B 2015 Bull.Mineral.Petrol.Geochem.34 18(in Chinese)[林緬,江文濱,李勇,易智星,張召彬2015礦物巖石地球化學通報34 18]

[7]Yao T Y,Huang Y Z,Li J S 2012 Chin.J.Theor.Appl.Mech.44 990(in Chinese)[姚同玉,黃延章,李繼山2012力學學報44 990]

[8]Jiang B Y,Li Z P,Diwu P X,Liu G,Wang J N,Gan H H 2014 Sci.Technol.Eng.14 58(in Chinese)[姜寶益,李治平,第五鵬祥,劉剛,王建寧,甘火華2014科學技術與工程14 58]

[9]Swami V,Settari A 2012 SPE Americas Unconventional Resources Conference Pittsburgh,USA,June 5–7,2012 p1

[10]Javadpour F,Fisher D,Unsworth M 2007 J.Can.Petrol.Technol.46 55

[11]Javadpour F 2009 J.Can.Petrol.Technol.48 16

[12]Shabro V,Torres-Verdin C,Sepehrnoori K 2012 SPE Annual Technical Conference and Exhibition San Antonio,USA,October 8–10,2012 p1

[13]Wu K,Li X,Wang C C,Chen Z X,Yu W 2015 AlChE J.61 2079

[14]Wu K L,Li X F,Chen Z X 2015 Acta Petrol.Sin.36 837(in Chinese)[吳克柳,李相方,陳掌星2015石油學報36 837]

[15]Haghshenas B,Clarkson C R,Chen S 2013 SPE Unconventional Resources Conference Calgary,Canada,November 5–7,2013 p1

[16]Song H Q,Liu Q P,Yu M X,Wu P,Zhang Y 2014 J.Univ.Sci.Technol.Beijing 36 139(in Chinese)[宋洪慶,劉啟鵬,于明旭,吳鵬,張雨2014北京科技大學學報36 139]

[17]Mi L D,Jiang H Q,Li J J 2014 J.Natur.Gas Sci.Eng.20 74

[18]Mi L D,Jiang H Q,Li J J,Tian Y 2014 Acta Petrol.Sin.35 928(in Chinese)[糜利棟,姜漢橋,李俊鍵,田野2014石油學報35 928]

[19]Ai S,Cheng L S,Huang S J,Fu L B,Wang T,Du B J,Zhang J,Liu H J 2014 Nat.Gas Geosci.25 166(in Chinese)[艾爽,程林松,黃世軍,傅禮兵,王濤,杜保健,張晉,劉紅君2014天然氣地球科學25 166]

[20]Li Y,Li X,Shi J,Wang H,Wu L 2014 SPE Biennial Energy Resources Conf.Port of Spain,Trinidad,June 9–11,2014 p1

[21]Jia P,Cheng L S,Huang S J,Xue Y C 2015 SPE Asia Paci fi c Unconventional Resources Conference and Exhibition Brisbane,Australia,November 9–11 2015 p1

[22]Zhou X,Zhang S C,Ma X F,Zhang Y 2015 Xinjiang Petrol.Geol.36 612(in Chinese)[周祥,張士誠,馬新仿,張燁2015新疆石油地質36 612]

[23]Roy S,Raju R,Chuang H F,Cruden B A,Meyyappan M 2003 J.Appl.Phys.93 4870

[24]Huang T,Guo X,Wang K 2015 J.Chem.2015 1

[25]Klinkenberg L J 1941 Drilling and Production Practice(New York:American Petroleum Institute)p200

[26]Li C L 2007 Natur.Gas Ind.27 85(in Chinese)[李傳亮2007天然氣工業27 85]

[27]Wu Q F,Chen W F,Huang L,Shi Y Z 2004 Rare fi ed Gas Dynamics(1st Ed.)(Changsha:National University of Defense Technology Press)pp96–97(in Chinese)[吳其芬,陳偉芳,黃琳,石于中2004稀薄氣體動力學(第一版)(長沙:國防科技大學出版社)第96—97頁]

[28]Zhu G Y,Liu X G,Li S T,Huang Y Z,Hao M Y 2007 Nat.Gas Ind.27 44(in Chinese)[朱光亞,劉先貴,李樹鐵,黃延章,郝明強2007天然氣工業27 44]

[29]Bear J(translated by Li J S,Chen C X)1983 Dynamics of Fluids in Porous Media(Beijing:China Architecture&Building Press)pp98–99(in Chinese)[貝爾J著(李競生,陳崇希譯)1983多孔介質流體動力學(北京:中國建筑工業出版社)第98—99頁]

[30]Ge H K,Shen Y H,Song Y,Wang X Q,Jiang C F,Shi P,Wang H,Yang L 2014 Nat.Gas Ind.34 46(in Chinese)[葛洪魁,申潁浩,宋巖,王小瓊,姜呈馥,史鵬,王暉,楊柳2014天然氣工業34 46]

[31]Chen D X 2002 Chin.J.Theor.Appl.Mech.34 96(in Chinese)[陳代珣2002力學學報34 96]

[32]Brown G P,DiNardo A,Cheng G K,Sherwood T K 1946 J.Appl.Phys.17 802

[33]Nie B S,Zhang L,Ma W F 2000 Coal Geol.Explor.28 20(in Chinese)[聶百勝,張力,馬文芳2000煤田地質與勘探28 20]

[34]Ziarani A S,Aguilera R 2012 Transp.Porous Med.91 239

[35]Florence F A,Rushing J A,Newsham K E,Blasingame T A 2007 Rocky Mountain Oil&Gas Technology Symposium(Denver:Society of Petroleum Engineers)p1

[36]Tian L,Xiao C,Liu M J,Gu D H 2014 J.Northeast Petrol.Univ.38 93(in Chinese)[田冷,肖聰,劉明進,顧岱紅2014東北石油大學學報38 93]

[37]Ren F,Wang X H,Ren K,Li Q Q,Wang S S,Zhang L,Wang Y 2013 Fault-Block Oil Gas Field 20 649(in Chinese)[任飛,王新海,任凱,李清泉,王珊珊,張磊,王宇2013斷塊油氣田20 649]

[38]Loucks R G,Reed R M,Ruppel S C,Jarvie D M 2009 J.Sediment.Res.79 848

[39]Civan F 2010 Porous Media and Its Applications in Science,Engineering,and Industry:3rd International Conference(Montecatini:AIP Publishing)p53

[40]Ren N,Shu L,Hu Y Q,Zhao J Z 2014 J.Southwest Petrol.Univ.(Sci.Technol Ed.)36 111(in Chinese)[任嵐,舒亮,胡永全,趙金洲2014西南石油大學學報(自然科學版)36 111]

[41]Song F Q,Zhang X,Huang X H,Long Y Q 2016 Sci.Sin.Technol.46 120(in Chinese)[宋付權,張翔,黃小荷,龍運前2016中國科學:技術科學46 120]

[42]Li Q Q,Wang X H,Yin H,Ren F,Wang S S,Zhang L 2013 J.Northeast Petrol.Univ.37 91(in Chinese)[李清泉,王新海,尹虎,任飛,王珊珊,張磊2013東北石油大學學報37 91]

[43]Chen Q,Kang Y L,You L J,Yu Y F,Liu H L 2013 Nat.Gas Geosci.24 1298(in Chinese)[陳強,康毅力,游利軍,俞楊烽,劉洪林2013天然氣地球科學24 1298]

[44]Wu K L,Li X F,Chen Z X 2015 Sci.Sin.Technol.45 525(in Chinese)[吳克柳,李相方,陳掌星2015中國科學:技術科學45 525]

[45]Sheng M,Li G S,Huang Z W,Tian S C 2014 Acta Petrol.Sin.35 347(in Chinese)[盛茂,李根生,黃中偉,田守嶒2014石油學報35 347]

[46]Do H D,Do D D,Prasetyo I 2001 AlChE J.47 2515

[47]Clarkson C R,Haghshenas B 2013 SPE Unconventional Resources Conference Woodlands,USA,April 10–12,2013 p1

[48]Civan F,Rai C S,Sondergeld C H 2011 Transport Porous Med.86 925

[49]Wang H Y,Ajao O,Economides M J 2014 J.Nat.Gas Sci.Eng.21 874

[50]Guo W,Xiong W,Gao S S,Hu Z M 2013 J.Cent.South Univ.(Sci.Technol.)44 2836(in Chinese)[郭為,熊偉,高樹生,胡志明2013中南大學學報(自然科學版)44 2836]

[51]Sladek K J,Gilliland E R,Baddour R F 1974 Ind.Eng.Chem.Fundam.13 100

[52]Curtis J B 2002 AAPG Bull.86 1921

[53]Kuuskraa V A,Koperna G,Schmoker J W,Quinn J C 1998 Oil Gas J.96 67

[54]Wang F P,Reed R M 2009 SPE Annual Technical Conference and Exhibition(New Orleans:Society of Petroleum Engineers)p1

[55]Arkilic E B 1997 Ph.D.Dissertation(Cambridge:Massachusetts Institute of Technology)

[56]Mosher K,He J,Liu Y,Rupp E,Wilcox J 2013 Int.J.Coal Geol.109 36

[57]Dubinin M M,Astakhov V A 1971 Russ.Chem.Bull.20 8

[58]Zhai Z Q,Wang X Q,Jin X,Sun L,Li J M,Cao D P 2014 Energ.Fuel.28 7467

[59]Wang H,Wang X Q,Jin X,Cao D P 2016 J.Phys.Chem.C 120 8986

[60]Kuila U,Prasad M,Kazemi H 2012 Proceedings of 9th Biennial InternationalConference and Exposition on Petroleum Geophysics Hyderabad,India,Febrary 16–18,2012 p182

PACS:47.61.–k,47.57.eb,47.45.Gx,51.10.+yDOI:10.7498/aps.66.114702

A new method for the transport mechanism coupling of shale gas slippage and di ff usion?

Li Ya-Xiong1)2)3)?Liu Xian-Gui2)3)Hu Zhi-Ming2)3)Gao Shu-Sheng2)3)Duan Xiang-Gang2)3)Chang Jin2)3)
1)(University of Chinese Academy of Sciences,Beijing 100049,China)
2)(Institute of Porous Flow and Fluid Mechanics,Chinese Academy of Sciences,Langfang 065007,China)
3)(Langfang Branch,PetroChina Research Institute of Petroleum Exploration and Development,Langfang 065007,China)

17 October 2016;revised manuscript

31 March 2017)

In view of the current status that di ff erent literature applies di ff erent coupling methods to the calculation of shale gas fl ow,and in order to clarify the relation between slippage and several di ff usions,in this paper the slippage e ff ect and various di ff usions are analyzed fi rst by theoretical analysis and mathematical models according to the de fi nitions and the mechanisms of microscopic motions.Afterwards,allowing for the spatial e ff ect of the adsorbed molecules on gas fl ow,the concept“wall-associated di ff usion”is proposed for the fi rst time to represent the gross e ff ects of Knudsen di ff usion and surface di ff usion,and it is pointed out that wall-associated di ff usion is equivalent to slippage e ff ect.Therefore a new coupling way where wall-associated di ff usion and slippage e ff ect are replaceable and no superposition of them is needed in fl ow calculation,is proposed.The case study shows that when the capillary radius ranges from 5 nm to 2000 nm,the relative error between wall-associated di ff usion and slippage e ff ect mass fl ux is fairly small,namely less than 10%in the vast majority of the range.The di ff erence between mean values of wall-associated di ff usion and slippage e ff ect mass fl ux in the whole aperture range is 1.4×10?6kg·m?2·s?1.That is,the relative error between the mean values is only 5.8%.Therefore,the new method satis fi es the requirements for engineering calculations.Taking parameter selection,un fi nished improvements in mathematical models of relevant mechanisms and other factors into account,there is some room for further promoting the veri fi cation of the proposed method.The development of wall-associated di ff usion has practical signi fi cance and multiple research signi fi cance.And the new coupling way reveals the relation between slippage and di ff usions,which prevents reduplicated superposition of shale gas fl ow mechanisms in nano-scale pores and can well change the status where the current coupling methods for shale gas fl ow are not consistent,thus specifying a new direction in the quantitative calculations for shale gas development.

shale gas,di ff usion,wall-associated di ff usion,slippage

10.7498/aps.66.114702

?國家重點基礎研究發展計劃(批準號:2013CB228000)資助的課題.

?通信作者.E-mail:18510284051@163.com

?2017中國物理學會Chinese Physical Society

http://wulixb.iphy.ac.cn

*Project supported by the National Basic Research Program of China(Grant No.2013CB228000).

?Corresponding author.E-mail:18510284051@163.com