• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    頁巖氣滑脫、擴(kuò)散傳輸機(jī)理耦合新方法?

    2017-08-09 07:34:24李亞雄1劉先貴2胡志明2高樹生2端祥剛2常進(jìn)2
    物理學(xué)報(bào) 2017年11期
    關(guān)鍵詞:通量孔徑頁巖

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

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

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

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

    頁巖氣滑脫、擴(kuò)散傳輸機(jī)理耦合新方法?

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

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

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

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

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

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

    頁巖氣,擴(kuò)散,壁聯(lián)擴(kuò)散,滑脫

    1 引言

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

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

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

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

    2 理論基礎(chǔ)與新方法的提出

    2.1 滑脫傳輸表征

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

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

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

    圖2 滑脫效應(yīng)示意圖Fig.2.Diagrammatic sketch of slippage e ff ect.

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

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

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

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

    所以滑脫流態(tài)下的總質(zhì)量通量為

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

    因此,相較于達(dá)西連續(xù)流動,由于滑脫效應(yīng)而導(dǎo)致的流通量為

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

    文獻(xiàn)[10]認(rèn)為多孔介質(zhì)由一定百分比的開放孔隙(孔隙度)和導(dǎo)致流通路徑長于直線路徑的互連網(wǎng)絡(luò)(迂曲度)所組成,則Jh的修正形式如下:

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

    2.2 擴(kuò)散傳輸表征

    2.2.1 克努森擴(kuò)散

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

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

    納米孔隙中克努森擴(kuò)散可以表達(dá)成壓力梯度的函數(shù),忽略黏性作用,頁巖氣質(zhì)量通量為[10]

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

    克努森擴(kuò)散系數(shù)DK定義式為[10,34?40]

    同理,用?/τ對克努森擴(kuò)散通量進(jìn)行修正:

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

    2.2.2 菲克擴(kuò)散

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

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

    菲克第一定律表示每個時(shí)刻基質(zhì)塊中各處氣體濃度相等,氣體擴(kuò)散速度與內(nèi)外濃度差呈線性關(guān)系.擴(kuò)散方程為[1,42]

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

    菲克第二定律表示基質(zhì)塊中氣體濃度分布為時(shí)間和位置的函數(shù),擴(kuò)散方程為[42]

    式中?為基質(zhì)孔隙度.

    2.2.3 過渡擴(kuò)散

    當(dāng)0.1

    Mi等[17]、糜利棟等[18]、陳強(qiáng)等[43]由菲克擴(kuò)散和克努森擴(kuò)散擴(kuò)散系數(shù)的調(diào)和平均得到過渡擴(kuò)散的擴(kuò)散系數(shù).而文獻(xiàn)[40]運(yùn)用通過實(shí)驗(yàn)研究得出的權(quán)重系數(shù)將黏性流和自由分子流進(jìn)行耦合得出過渡流的質(zhì)量流量.關(guān)于過渡流的模型正處于研究之中,目前還沒有一種被廣泛采用的過渡流表達(dá)式.

    2.2.4 表面擴(kuò)散

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

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

    依據(jù)文獻(xiàn)[45],表面擴(kuò)散質(zhì)量通量可寫為

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

    Sladek等[51]根據(jù)30組氣-固系統(tǒng)的數(shù)據(jù)擬合出的表面擴(kuò)散系數(shù)公式為

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

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

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

    所以,

    2.3 滑脫、擴(kuò)散機(jī)理耦合新方法

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

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

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

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

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

    故孔隙的有效半徑為

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

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

    克努森擴(kuò)散和表面擴(kuò)散都與氣-固相之間在界面處的作用有關(guān),其中克努森擴(kuò)散是氣體分子和壁面產(chǎn)生碰撞以后所發(fā)生的擴(kuò)散作用,表面擴(kuò)散是氣體分子在孔隙表面吸附位之間連續(xù)跳躍的過程,兩者的作用都是導(dǎo)致氣體分子在壁面周圍的運(yùn)移速度不為零,這與滑脫現(xiàn)象的本質(zhì)是一致的,如果提出新名詞——“壁聯(lián)擴(kuò)散”來表征表面擴(kuò)散和克努森擴(kuò)散的綜合效應(yīng),則可從以上微觀運(yùn)動機(jī)制的分析得出壁聯(lián)擴(kuò)散與滑脫效應(yīng)等同的結(jié)論,并提出一種滑脫和擴(kuò)散耦合方法的新認(rèn)識,即在頁巖氣納米孔中壁聯(lián)擴(kuò)散和滑脫效應(yīng)對總流動的貢獻(xiàn)相同,亦即在流動計(jì)算中兩種機(jī)制互不疊加且可相互替換,如圖6所示.所以,在目前研究中出現(xiàn)的將克努森擴(kuò)散和滑脫進(jìn)行疊加耦合,或?qū)⒒摵捅砻鏀U(kuò)散進(jìn)行疊加,抑或?qū)⑷叨伎紤]在傳輸機(jī)制中的計(jì)算方法,都將導(dǎo)致重復(fù)疊加的問題.出現(xiàn)這種現(xiàn)象的原因歸根究底源于人們對滑脫和幾種擴(kuò)散之間相互聯(lián)系認(rèn)識的缺失,從而因?yàn)楦鞣N機(jī)理名稱的不同各持己見所致.壁聯(lián)擴(kuò)散的概念以及耦合新方法的提出能夠厘清滑脫和各種擴(kuò)散之間的關(guān)系,在頁巖氣計(jì)算中統(tǒng)一所需考慮的流動機(jī)制,避免納米級孔隙中流動機(jī)理的重復(fù)疊加,從而較好地改變頁巖氣流動耦合方法不一致的現(xiàn)狀,為頁巖氣開發(fā)工作中的定量計(jì)算指明新方向.

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

    圖6 壁聯(lián)擴(kuò)散和滑脫關(guān)系圖Fig.6.Relationship between wall-associated di ff usion and slippage e ff ect.

    其中Jb為壁聯(lián)擴(kuò)散質(zhì)量通量(kg·m?2·s?1).

    3 耦合方法實(shí)例分析

    依據(jù)文獻(xiàn)[45,52—55],參考典型頁巖氣藏Barnett頁巖參數(shù)作為基礎(chǔ)數(shù)據(jù):毛細(xì)管半徑為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;由溫壓條件根據(jù)REFPROP軟件得ρa(bǔ)vg=84.176 kg·m?3.分別得出質(zhì)量通量隨孔徑的變化曲線、滑脫和壁聯(lián)擴(kuò)散質(zhì)量通量平均值對比圖和質(zhì)量通量相對誤差隨孔徑的變化曲線如圖7—圖9所示.

    3.1結(jié)果與討論

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

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

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

    圖8 滑脫和壁聯(lián)擴(kuò)散質(zhì)量通量平均值對比Fig.8.Average values of mass fl ux of slippage and wall-associated di ff usion.

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

    綜上理論和實(shí)例分析都表明,壁聯(lián)擴(kuò)散和滑脫效應(yīng)在流動計(jì)算中可互換而不重復(fù)疊加的耦合新方法是合理的,在頁巖氣藏流動分析中,視滑脫效應(yīng)和壁聯(lián)擴(kuò)散在數(shù)量上等同可以滿足工程計(jì)算的需要.

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

    4 結(jié)論

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

    [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頁巖氣藏開發(fā)基礎(chǔ)(第一版)(北京:石油工業(yè)出版社)第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天然氣地球科學(xué)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中國科學(xué):技術(shù)科學(xué)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礦物巖石地球化學(xué)通報(bào)34 18]

    [7]Yao T Y,Huang Y Z,Li J S 2012 Chin.J.Theor.Appl.Mech.44 990(in Chinese)[姚同玉,黃延章,李繼山2012力學(xué)學(xué)報(bào)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科學(xué)技術(shù)與工程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石油學(xué)報(bào)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北京科技大學(xué)學(xué)報(bào)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石油學(xué)報(bào)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天然氣地球科學(xué)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新疆石油地質(zhì)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天然氣工業(yè)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稀薄氣體動力學(xué)(第一版)(長沙:國防科技大學(xué)出版社)第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)[朱光亞,劉先貴,李樹鐵,黃延章,郝明強(qiáng)2007天然氣工業(yè)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多孔介質(zhì)流體動力學(xué)(北京:中國建筑工業(yè)出版社)第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天然氣工業(yè)34 46]

    [31]Chen D X 2002 Chin.J.Theor.Appl.Mech.34 96(in Chinese)[陳代珣2002力學(xué)學(xué)報(bào)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煤田地質(zhì)與勘探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)[田冷,肖聰,劉明進(jìn),顧岱紅2014東北石油大學(xué)學(xué)報(bào)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西南石油大學(xué)學(xué)報(bào)(自然科學(xué)版)36 111]

    [41]Song F Q,Zhang X,Huang X H,Long Y Q 2016 Sci.Sin.Technol.46 120(in Chinese)[宋付權(quán),張翔,黃小荷,龍運(yùn)前2016中國科學(xué):技術(shù)科學(xué)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東北石油大學(xué)學(xué)報(bào)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)[陳強(qiáng),康毅力,游利軍,俞楊烽,劉洪林2013天然氣地球科學(xué)24 1298]

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

    [45]Sheng M,Li G S,Huang Z W,Tian S C 2014 Acta Petrol.Sin.35 347(in Chinese)[盛茂,李根生,黃中偉,田守嶒2014石油學(xué)報(bào)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中南大學(xué)學(xué)報(bào)(自然科學(xué)版)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

    ?國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(批準(zhǔn)號:2013CB228000)資助的課題.

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

    ?2017中國物理學(xué)會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

    猜你喜歡
    通量孔徑頁巖
    冬小麥田N2O通量研究
    不同滲透率巖芯孔徑分布與可動流體研究
    分布式孔徑相參合成雷達(dá)技術(shù)
    頁巖氣開發(fā)降溫
    能源(2016年1期)2016-12-01 05:10:02
    基于子孔徑斜率離散采樣的波前重構(gòu)
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    大孔徑淺臺階控制爆破在重慶地區(qū)的應(yīng)用
    重慶建筑(2014年12期)2014-07-24 14:00:32
    春、夏季長江口及鄰近海域溶解甲烷的分布與釋放通量
    我國頁巖氣可采資源量初步估計(jì)為31萬億m3
    頁巖氣是什么?
    99riav亚洲国产免费| 51午夜福利影视在线观看| 亚洲 国产 在线| 看免费av毛片| 久久九九热精品免费| 久热这里只有精品99| 国产在视频线精品| e午夜精品久久久久久久| 色综合婷婷激情| 久热爱精品视频在线9| 一进一出抽搐动态| 精品卡一卡二卡四卡免费| 亚洲精品中文字幕在线视频| 国产成人影院久久av| 老司机亚洲免费影院| 欧美日韩乱码在线| 美国免费a级毛片| av有码第一页| 亚洲成国产人片在线观看| www日本在线高清视频| 精品国内亚洲2022精品成人 | 国产在线精品亚洲第一网站| 丝瓜视频免费看黄片| 狠狠婷婷综合久久久久久88av| 久久草成人影院| 亚洲专区中文字幕在线| 女警被强在线播放| 一本一本久久a久久精品综合妖精| 女人精品久久久久毛片| 新久久久久国产一级毛片| 亚洲av片天天在线观看| 国产精品偷伦视频观看了| 欧美精品高潮呻吟av久久| 成人三级做爰电影| 日本黄色日本黄色录像| 亚洲欧美精品综合一区二区三区| 99久久精品国产亚洲精品| 亚洲精品美女久久久久99蜜臀| 色婷婷久久久亚洲欧美| 日本一区二区免费在线视频| 久久久国产欧美日韩av| 黑人操中国人逼视频| 久久人人爽av亚洲精品天堂| 日韩有码中文字幕| 精品久久久久久久毛片微露脸| 丁香欧美五月| 久久久久久久国产电影| 亚洲午夜精品一区,二区,三区| 在线av久久热| 亚洲五月天丁香| 日韩熟女老妇一区二区性免费视频| 999精品在线视频| 欧美不卡视频在线免费观看 | 变态另类成人亚洲欧美熟女 | 纯流量卡能插随身wifi吗| 国产成人精品久久二区二区免费| 大片电影免费在线观看免费| 午夜精品久久久久久毛片777| 最近最新中文字幕大全免费视频| 一二三四社区在线视频社区8| 乱人伦中国视频| 国产精品99久久99久久久不卡| 久久人人爽av亚洲精品天堂| 91在线观看av| 成熟少妇高潮喷水视频| 91国产中文字幕| 狠狠狠狠99中文字幕| 美女高潮到喷水免费观看| 欧美在线黄色| 水蜜桃什么品种好| 搡老熟女国产l中国老女人| 国产人伦9x9x在线观看| bbb黄色大片| 一区在线观看完整版| 免费看十八禁软件| 极品教师在线免费播放| 欧美老熟妇乱子伦牲交| 久久人人97超碰香蕉20202| 一级a爱视频在线免费观看| 成人精品一区二区免费| 制服人妻中文乱码| 热99re8久久精品国产| 十八禁人妻一区二区| 亚洲国产精品一区二区三区在线| 国产一卡二卡三卡精品| 国产真人三级小视频在线观看| 日韩免费av在线播放| 不卡一级毛片| 亚洲专区国产一区二区| 亚洲av片天天在线观看| 久久久久国产一级毛片高清牌| 国产精品久久视频播放| 免费观看人在逋| 亚洲中文av在线| 黄片小视频在线播放| 99国产极品粉嫩在线观看| 超碰97精品在线观看| 精品国产超薄肉色丝袜足j| 亚洲中文日韩欧美视频| 亚洲欧洲精品一区二区精品久久久| 国产精品美女特级片免费视频播放器 | 两人在一起打扑克的视频| 99在线人妻在线中文字幕 | 男人的好看免费观看在线视频 | 亚洲欧美激情综合另类| av超薄肉色丝袜交足视频| 亚洲精品粉嫩美女一区| 热99久久久久精品小说推荐| 中文字幕制服av| 精品国产亚洲在线| 亚洲国产欧美网| 日本wwww免费看| 国产欧美日韩一区二区精品| 国产精品久久视频播放| 黄频高清免费视频| 国内毛片毛片毛片毛片毛片| 国产成人精品无人区| 狠狠狠狠99中文字幕| 国产一区二区激情短视频| 国产一卡二卡三卡精品| 在线观看免费视频网站a站| 极品少妇高潮喷水抽搐| 国产欧美日韩一区二区三区在线| 久久精品亚洲精品国产色婷小说| 最近最新中文字幕大全电影3 | 精品人妻1区二区| 欧美精品亚洲一区二区| 午夜福利影视在线免费观看| 久久国产亚洲av麻豆专区| 国产1区2区3区精品| 欧美最黄视频在线播放免费 | videosex国产| 久久中文看片网| a级片在线免费高清观看视频| 亚洲成人免费av在线播放| 久久人人97超碰香蕉20202| 久久久精品区二区三区| 曰老女人黄片| 亚洲午夜理论影院| 日日爽夜夜爽网站| 亚洲色图综合在线观看| 高潮久久久久久久久久久不卡| 久久久久精品国产欧美久久久| 国产97色在线日韩免费| 国产熟女午夜一区二区三区| 亚洲欧美激情综合另类| 又大又爽又粗| 国产免费现黄频在线看| 美女午夜性视频免费| 老熟女久久久| 一进一出抽搐gif免费好疼 | 久久精品亚洲精品国产色婷小说| 精品国内亚洲2022精品成人 | 亚洲一区二区三区不卡视频| 黑人操中国人逼视频| 老司机亚洲免费影院| 精品久久久精品久久久| 亚洲中文日韩欧美视频| 在线观看免费高清a一片| 亚洲av片天天在线观看| 天天躁日日躁夜夜躁夜夜| 日韩免费高清中文字幕av| 亚洲专区字幕在线| 搡老熟女国产l中国老女人| 少妇猛男粗大的猛烈进出视频| 丝袜美足系列| 成人永久免费在线观看视频| 午夜91福利影院| 老司机深夜福利视频在线观看| 黄色毛片三级朝国网站| 成人18禁高潮啪啪吃奶动态图| 老司机福利观看| 国产精品欧美亚洲77777| 日韩欧美在线二视频 | 妹子高潮喷水视频| 精品少妇一区二区三区视频日本电影| 麻豆乱淫一区二区| 在线观看66精品国产| 又黄又爽又免费观看的视频| 在线十欧美十亚洲十日本专区| av电影中文网址| 欧美激情久久久久久爽电影 | 国产av一区二区精品久久| 黄色女人牲交| 国产免费av片在线观看野外av| 午夜福利,免费看| 国产三级黄色录像| 亚洲人成77777在线视频| 人人妻人人澡人人爽人人夜夜| 老司机影院毛片| 精品免费久久久久久久清纯 | 啦啦啦在线免费观看视频4| 国产精品美女特级片免费视频播放器 | 99久久综合精品五月天人人| av福利片在线| 午夜福利乱码中文字幕| 亚洲国产欧美日韩在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 久久精品国产a三级三级三级| 亚洲专区中文字幕在线| 少妇 在线观看| av国产精品久久久久影院| 亚洲精品美女久久久久99蜜臀| 午夜91福利影院| 成人永久免费在线观看视频| 国产深夜福利视频在线观看| 久久香蕉国产精品| 人人妻人人澡人人爽人人夜夜| 亚洲中文av在线| 老汉色av国产亚洲站长工具| 欧美黄色片欧美黄色片| 欧美在线黄色| 午夜福利乱码中文字幕| 国产高清激情床上av| 精品一区二区三卡| 亚洲精品久久成人aⅴ小说| 天天躁狠狠躁夜夜躁狠狠躁| 女人高潮潮喷娇喘18禁视频| 身体一侧抽搐| 国产精品久久久久久精品古装| 一区二区三区激情视频| 中文字幕另类日韩欧美亚洲嫩草| 美女高潮喷水抽搐中文字幕| 欧美激情极品国产一区二区三区| 亚洲中文av在线| 啦啦啦 在线观看视频| 久久久久久久午夜电影 | 天堂俺去俺来也www色官网| 久久青草综合色| 男男h啪啪无遮挡| 两性午夜刺激爽爽歪歪视频在线观看 | 成人特级黄色片久久久久久久| 妹子高潮喷水视频| 在线观看一区二区三区激情| 国产亚洲精品第一综合不卡| 在线视频色国产色| 国产精品免费一区二区三区在线 | 亚洲av第一区精品v没综合| 一边摸一边抽搐一进一小说 | 成人特级黄色片久久久久久久| 久久人人97超碰香蕉20202| 一级a爱片免费观看的视频| a级片在线免费高清观看视频| 老司机影院毛片| 视频在线观看一区二区三区| 日韩有码中文字幕| 人妻久久中文字幕网| 看免费av毛片| 久久精品亚洲熟妇少妇任你| 亚洲欧美激情在线| 日韩欧美国产一区二区入口| 美女午夜性视频免费| 日韩欧美一区视频在线观看| 亚洲熟妇熟女久久| 亚洲精品中文字幕一二三四区| 亚洲av欧美aⅴ国产| bbb黄色大片| 搡老乐熟女国产| av免费在线观看网站| 婷婷精品国产亚洲av在线 | 欧美成人免费av一区二区三区 | 亚洲专区国产一区二区| 大型黄色视频在线免费观看| 一进一出好大好爽视频| 丝袜美腿诱惑在线| 日韩大码丰满熟妇| 成人国语在线视频| 两个人看的免费小视频| 国产成+人综合+亚洲专区| 成人特级黄色片久久久久久久| av线在线观看网站| 国产av一区二区精品久久| 国产精品av久久久久免费| 香蕉国产在线看| 国产亚洲欧美98| 人人妻人人添人人爽欧美一区卜| 精品福利观看| 精品亚洲成国产av| av视频免费观看在线观看| 久久国产精品大桥未久av| 如日韩欧美国产精品一区二区三区| 久久久国产一区二区| 色综合婷婷激情| 国产精品久久久久久精品古装| av欧美777| 亚洲成av片中文字幕在线观看| 国产精品香港三级国产av潘金莲| 国产免费av片在线观看野外av| 国产精品香港三级国产av潘金莲| 国产一区二区三区综合在线观看| 亚洲五月天丁香| 日韩一卡2卡3卡4卡2021年| 精品一区二区三区av网在线观看| 亚洲一区中文字幕在线| 久久久久久久国产电影| 亚洲av美国av| 麻豆av在线久日| 欧美一级毛片孕妇| 99国产精品99久久久久| 黑人巨大精品欧美一区二区mp4| 久久国产乱子伦精品免费另类| 国产精品九九99| 成人手机av| 国内久久婷婷六月综合欲色啪| 久久精品亚洲精品国产色婷小说| 12—13女人毛片做爰片一| 1024香蕉在线观看| 满18在线观看网站| 首页视频小说图片口味搜索| 久久精品国产综合久久久| 日韩免费av在线播放| 人妻一区二区av| 麻豆av在线久日| 成熟少妇高潮喷水视频| 美国免费a级毛片| 婷婷精品国产亚洲av在线 | 法律面前人人平等表现在哪些方面| 成人永久免费在线观看视频| 日韩制服丝袜自拍偷拍| 亚洲精品美女久久av网站| 一级片免费观看大全| 中出人妻视频一区二区| 成人永久免费在线观看视频| 乱人伦中国视频| 成年人免费黄色播放视频| 18在线观看网站| 熟女少妇亚洲综合色aaa.| 欧美 亚洲 国产 日韩一| 亚洲伊人色综图| 亚洲av成人av| 国产成人影院久久av| www.熟女人妻精品国产| 国产精品98久久久久久宅男小说| 老汉色av国产亚洲站长工具| 最近最新中文字幕大全免费视频| 午夜精品在线福利| 热99re8久久精品国产| 久久久久久久久免费视频了| 亚洲自偷自拍图片 自拍| 老司机靠b影院| 日本精品一区二区三区蜜桃| 亚洲五月天丁香| 中出人妻视频一区二区| 国产亚洲精品一区二区www | 久久久国产精品麻豆| 国产成人影院久久av| 成熟少妇高潮喷水视频| 亚洲午夜理论影院| 欧美日韩亚洲高清精品| 亚洲第一青青草原| 久久天躁狠狠躁夜夜2o2o| 少妇的丰满在线观看| 看黄色毛片网站| 一区二区三区激情视频| 国产免费av片在线观看野外av| 日韩欧美免费精品| 最新在线观看一区二区三区| 中文字幕高清在线视频| 嫩草影视91久久| 五月开心婷婷网| 精品无人区乱码1区二区| 热re99久久国产66热| 国产亚洲精品第一综合不卡| 亚洲av电影在线进入| 精品无人区乱码1区二区| 国产成人精品久久二区二区免费| 国产日韩欧美亚洲二区| 叶爱在线成人免费视频播放| 最新在线观看一区二区三区| 老司机亚洲免费影院| 午夜福利一区二区在线看| 超色免费av| 99精品在免费线老司机午夜| 日韩制服丝袜自拍偷拍| 少妇裸体淫交视频免费看高清 | 亚洲av片天天在线观看| 亚洲专区中文字幕在线| 欧美日韩中文字幕国产精品一区二区三区 | 久久精品成人免费网站| 国产精品欧美亚洲77777| 亚洲精品成人av观看孕妇| 91精品三级在线观看| 青草久久国产| 丰满迷人的少妇在线观看| 飞空精品影院首页| 久久久国产精品麻豆| 久久久精品国产亚洲av高清涩受| 天天躁日日躁夜夜躁夜夜| 美女午夜性视频免费| 免费看十八禁软件| 99精品在免费线老司机午夜| 午夜视频精品福利| 亚洲第一青青草原| 99香蕉大伊视频| 成人精品一区二区免费| 黑人巨大精品欧美一区二区蜜桃| 亚洲人成77777在线视频| 99热网站在线观看| 成人精品一区二区免费| www.自偷自拍.com| 狠狠婷婷综合久久久久久88av| 自拍欧美九色日韩亚洲蝌蚪91| 九色亚洲精品在线播放| 精品国产一区二区三区四区第35| 久久久水蜜桃国产精品网| 又大又爽又粗| 最近最新中文字幕大全电影3 | 久久草成人影院| www日本在线高清视频| 免费少妇av软件| 日韩大码丰满熟妇| 亚洲精品久久午夜乱码| 午夜精品国产一区二区电影| 成年人免费黄色播放视频| 欧美黄色淫秽网站| 欧美乱色亚洲激情| 午夜福利视频在线观看免费| 国产免费现黄频在线看| 国产成人av教育| 久久精品亚洲精品国产色婷小说| av中文乱码字幕在线| 美女高潮到喷水免费观看| 久99久视频精品免费| 啦啦啦在线免费观看视频4| 国产无遮挡羞羞视频在线观看| 免费高清在线观看日韩| 国产亚洲欧美98| 韩国av一区二区三区四区| 乱人伦中国视频| 亚洲av欧美aⅴ国产| 黄片播放在线免费| 天天躁夜夜躁狠狠躁躁| 精品免费久久久久久久清纯 | 午夜精品国产一区二区电影| 狠狠婷婷综合久久久久久88av| 最新美女视频免费是黄的| 老司机靠b影院| 国产精品av久久久久免费| 久久精品国产99精品国产亚洲性色 | 人成视频在线观看免费观看| 无人区码免费观看不卡| 色播在线永久视频| 欧美另类亚洲清纯唯美| 脱女人内裤的视频| 男人舔女人的私密视频| 成人亚洲精品一区在线观看| 99久久综合精品五月天人人| 色老头精品视频在线观看| 国产主播在线观看一区二区| 免费观看人在逋| 老鸭窝网址在线观看| 国产一区二区激情短视频| 50天的宝宝边吃奶边哭怎么回事| av天堂在线播放| 亚洲黑人精品在线| 久久久国产成人免费| 黑人操中国人逼视频| 在线永久观看黄色视频| 久久精品成人免费网站| av天堂久久9| ponron亚洲| 色94色欧美一区二区| 精品无人区乱码1区二区| 国产精品乱码一区二三区的特点 | 精品久久久久久电影网| 黄色毛片三级朝国网站| 日本欧美视频一区| 久久精品成人免费网站| 国产精品电影一区二区三区 | 一边摸一边抽搐一进一小说 | 亚洲av日韩精品久久久久久密| 久久久久视频综合| 999精品在线视频| 国产伦人伦偷精品视频| 搡老乐熟女国产| 国产日韩欧美亚洲二区| 欧美成狂野欧美在线观看| 亚洲av片天天在线观看| 国产欧美日韩一区二区三| 国产精品秋霞免费鲁丝片| 精品福利观看| 久久国产精品大桥未久av| 国产日韩一区二区三区精品不卡| 欧美av亚洲av综合av国产av| 精品熟女少妇八av免费久了| 变态另类成人亚洲欧美熟女 | 韩国av一区二区三区四区| 在线免费观看的www视频| 自拍欧美九色日韩亚洲蝌蚪91| 午夜精品国产一区二区电影| 国产精品久久视频播放| 一本大道久久a久久精品| 岛国毛片在线播放| 久久久久精品国产欧美久久久| 亚洲少妇的诱惑av| 欧美人与性动交α欧美精品济南到| 老司机在亚洲福利影院| 黑丝袜美女国产一区| 国产精品电影一区二区三区 | 精品久久蜜臀av无| 三上悠亚av全集在线观看| a在线观看视频网站| 欧美日韩瑟瑟在线播放| 久久久久久久午夜电影 | 亚洲avbb在线观看| 国产在线一区二区三区精| 亚洲人成电影观看| 下体分泌物呈黄色| 欧美日韩亚洲综合一区二区三区_| 亚洲在线自拍视频| 久久精品aⅴ一区二区三区四区| 最新在线观看一区二区三区| 欧美亚洲日本最大视频资源| 亚洲av熟女| 韩国av一区二区三区四区| 捣出白浆h1v1| 国产aⅴ精品一区二区三区波| 国产真人三级小视频在线观看| 午夜精品久久久久久毛片777| 色精品久久人妻99蜜桃| 久久精品国产综合久久久| 久久久久精品人妻al黑| 成人三级做爰电影| videos熟女内射| 久99久视频精品免费| 王馨瑶露胸无遮挡在线观看| 欧美黑人精品巨大| 俄罗斯特黄特色一大片| 人妻丰满熟妇av一区二区三区 | 国产一区二区三区在线臀色熟女 | 人人澡人人妻人| 黄色女人牲交| 777米奇影视久久| 亚洲少妇的诱惑av| 操出白浆在线播放| 欧美日韩黄片免| 超色免费av| 久久国产精品影院| 在线看a的网站| 女人爽到高潮嗷嗷叫在线视频| 国产一区二区激情短视频| 久久久久久久久免费视频了| 少妇被粗大的猛进出69影院| 国产男女内射视频| 亚洲欧美色中文字幕在线| 日本vs欧美在线观看视频| 天天躁夜夜躁狠狠躁躁| 欧美日韩视频精品一区| 亚洲欧美激情在线| 一级,二级,三级黄色视频| 激情在线观看视频在线高清 | 午夜日韩欧美国产| 大香蕉久久网| 一级毛片精品| 亚洲成av片中文字幕在线观看| 中文字幕最新亚洲高清| av视频免费观看在线观看| 亚洲综合色网址| 欧美在线黄色| 精品久久蜜臀av无| 免费黄频网站在线观看国产| 精品一区二区三区四区五区乱码| 久久久久久人人人人人| 久久久国产成人免费| 悠悠久久av| 女警被强在线播放| 18禁观看日本| 人人妻人人澡人人看| 中文字幕制服av| 国产国语露脸激情在线看| 久久久久久久国产电影| 99精品久久久久人妻精品| 老熟妇仑乱视频hdxx| 欧美日韩福利视频一区二区| 久久久久久人人人人人| 日韩有码中文字幕| 日韩制服丝袜自拍偷拍| 女警被强在线播放| 中文字幕制服av| 99热国产这里只有精品6| 天天躁夜夜躁狠狠躁躁| 精品熟女少妇八av免费久了| 欧美日韩亚洲高清精品| 一级作爱视频免费观看| 中文字幕制服av| 三级毛片av免费| 精品亚洲成a人片在线观看| 亚洲人成伊人成综合网2020| 99久久人妻综合| 99久久精品国产亚洲精品| 日韩精品免费视频一区二区三区| 黑人猛操日本美女一级片| 人人妻人人爽人人添夜夜欢视频| 搡老乐熟女国产| 人人妻人人爽人人添夜夜欢视频| 怎么达到女性高潮| 国产午夜精品久久久久久| 丁香六月欧美| 国产精品一区二区在线观看99| netflix在线观看网站| 少妇裸体淫交视频免费看高清 | 波多野结衣av一区二区av| 国产精品亚洲av一区麻豆| 国产一区有黄有色的免费视频| 亚洲综合色网址| 黄色女人牲交| 亚洲情色 制服丝袜| 岛国毛片在线播放| 韩国av一区二区三区四区| 免费观看人在逋| 美女国产高潮福利片在线看| 最近最新免费中文字幕在线| 中文亚洲av片在线观看爽 | 亚洲精品一卡2卡三卡4卡5卡| 18在线观看网站|