基于流動擴散互竟關系的基質吸附態瓦斯表觀擴散系數實驗室測定準確性分析

2022-04-06 07:50:38李成武劉晨熙魏誠敏袁廷璐

煤炭學報 2022年2期

趙偉，王凱，李成武，鞠楊，劉晨熙，魏誠敏，袁廷璐

(1.中國礦業大學(北京) 共伴生能源精準開采北京市重點實驗室，北京 100083；2. 中國礦業大學(北京) 應急管理與安全工程學院，北京 100083;3. 中國礦業大學(北京) 煤炭資源與安全開采國家重點實驗室，北京 100083)

擴散系數是描述煤中瓦斯運移特性的關鍵參數，其對煤層氣儲量計算、瓦斯災害危險評價以及溫室氣體排放估算等有著重要的意義。煤是由基質和裂隙組成的雙重孔隙介質，在原始應力環境中瓦斯分別在2個系統中以擴散和流動2種形式運移。在基質中，擴散以表面擴散、諾森擴散、過渡擴散和菲克擴散為主，其通常在尺度較小的孔隙中發生，以吸附瓦斯形成的濃度梯度為主要驅動力進行傳質；而裂隙尺度空間較大，傳質行為較擴散更為連續、劇烈，根據雷諾數的大小,流動又可分為層流、過渡流和紊流。為了描述方便，原位煤層裂隙中的瓦斯流動可簡單假設為符合達西方程的滲流(層流)，其以游離瓦斯形成的壓力梯度為主要驅動力進行運移。游離瓦斯和吸附瓦斯時刻處于動態平衡之中，兩者難以區分。

擴散現象本質上是由分子化學勢的不均而引起的，故而引起化學勢分布不均的因素均可以引起擴散行為的變化。一般而言，影響擴散的因素可分為2類：一類決定了瓦斯的賦存狀態，即擴散源濃度的大小，如壓力、含量、溫度和水分等因素；另一類則決定了瓦斯擴散路徑的難易程度，如煤的孔隙結構、粒徑和變質程度等因素。擴散是一種復雜的傳質現象，因此擴散系數也存在較多分類。文獻[1]中對擴散系數的分類有多種方法，以多孔介質復雜孔隙系統與單個均質孔道中瓦斯分子運移的差異性為分類依據，擴散系數可分為本質擴散系數、有效擴散系數和表觀擴散系數3類。本質擴散系數為氣體分子本身的擴散系數，與溫度等影響分子動能的特性有關，而與擴散路徑等外在因素無關。有效擴散系數是氣體分子經過多孔介質蜿蜒曲折的孔道所產生的擴散系數，其與多孔介質的曲折度、孔隙率和孔長等孔隙性質有關。表觀擴散系數則是大量分子在多孔介質中運移所產生的擴散系數，其不僅能反應單個分子的擴散特征，還從整體上考慮了大量分子與孔壁的相互作用以及分子與分子之間的相互作用。其不單單反應了擴散行為，也反應了流動對擴散的邊界影響。一般而言，通過解吸曲線反演得出的擴散系數，其融合了流動和擴散2種典型的傳質形式，其中又可細分為非線性紊流、線性層流(滲流)、菲克擴散、過渡擴散、諾森擴散以及表面擴散等多種傳質形式，是一種表征宏觀傳質特征的表觀擴散系數。

研究不同因素對煤中瓦斯擴散的影響機制時，常常通過改變常規解吸試驗的某一自變量，得到表觀擴散系數與其的關系曲線，進而進行分析。在實驗室吸附瓦斯表觀擴散系數測定方法中，通常以常規的解吸試驗為基礎，通過釋放游離瓦斯，測算解吸瓦斯曲線，進而得到解吸分數曲線，從而推算基質中吸附態瓦斯的表觀擴散系數。然而，在釋放游離瓦斯過程中，由于損失時間難以把握，游離瓦斯的不同排放時間長短對解吸曲線及擴散模型擬合產生很大的影響，常會得到不同的表觀擴散系數。另外，在現場測量鉆屑解吸指標時，由于是煤屑第1 min的瓦斯解吸量，其與表觀擴散系數存在一定的定量關系，也常因為鉆屑收取過程中損失時間的差異，測出的存在較大誤差。因此，有必要探究不同解吸損失時間下表觀擴散系數的變化規律，來修正表觀擴散系數測定的精確度。

筆者將游離瓦斯的排放時間定義為解吸損失時間，通過對游離瓦斯和吸附瓦斯的運移特征差異進行分析，構建引入游離瓦斯組分貢獻的瓦斯全濃度表觀擴散模型，進而分析損失時間對表觀擴散系數測量準確性的影響，并基于自主研發的瓦斯解吸特性自動檢測儀，試驗模擬測算不同解吸損失時間下的瓦斯解吸曲線和擴散分數曲線形態，得到表觀擴散系數有效值與損失時間的定量關系，從而指導精準化瓦斯表觀擴散系數測定。

1 試驗

1.1 試驗樣品

試驗煤樣取自內蒙古鄂爾多斯煤田神東集團柳塔礦2煤層，屬于中低變質程度煙煤，其相關物性參數見表1。

表1 煤樣基本參數

1.2 解吸試驗

本次解吸試驗采用自主研發的瓦斯解吸特性自動監測儀(圖1)。與常規通過量筒人工讀數相比較，該設備具有自動化、連續化和精準化數據讀取的優點。尤其對于解吸初期的解吸數據，避免了常規方法中由于短時間內大量游離瓦斯逸散，導致氣泡大量生成，引起人工讀數誤差的問題(圖2)。解吸試驗在0.5，1，1.5 MPa的不同平衡壓力條件下進行，得到0，1，5，10 min不同損失時間條件下的解吸曲線。具體試驗步驟如下：① 首先將現場取得的煤樣放入專用破碎機中破碎，之后用實驗篩進行篩分。取得煤粒大小在60～80目(0.178～0.250 mm)的煤樣。② 將試驗煤樣裝入清潔的動態測壓儀器中，進行真空脫氣。③ 由于實驗室安全要求，瓦斯以CO代替。脫氣8 h左右后向煤樣罐中充入CO，充至一定壓力后關閉閥門，然后進行吸附平衡試驗。④ 待吸附平衡后，連接解吸量自動監測儀，該儀器可每秒讀取1次解吸數據。采用計算機監測不同解吸時間的解吸體積，得出解吸曲線。重復上述試驗，測定不同壓力下的解吸曲線。⑤ 通過Matlab編程語言，對解吸曲線進行擴散模型擬合并得出表觀擴散系數。

圖1 自動解吸測定裝置

圖2 不同測試方法解吸初期數據對比示意

2 結果與分析

2.1 不同損失時間條件下的瓦斯解吸特性

不同損失時間條件下的瓦斯解吸特性是分析其擴散特性變化的基礎。為了避免逸散時間人工控制的誤差，在試驗時不同損失時間的解吸曲線通過截取特定時間段的解吸曲線形成。從圖3的試驗結果可以看出：和以往文獻中的試驗現象相似，平衡壓力越大，則固定時間段內的瓦斯解吸量就越大，且瓦斯初期解吸速度越快；而隨著時間推移，瓦斯解吸速度減少得也越快。而隨著損失時間的增加，解吸曲線的形態會發生顯著變化：曲線逐步遠離軸，并靠近軸，初期“直線”解吸段的長度越來越短，瓦斯解吸速度減小的幅度越來越小。曲線形態的變化，意味著控制曲線的基本傳質形式有所改變，有必要對其內部的主控傳質形式進行研究。

圖3 不同損失時間下的解吸曲線形態變化

2.2 不同損失時間條件下的瓦斯表觀擴散系數

(1)

式中，為時刻的解吸量，mL；為極限解吸量，mL；/為解吸分數；為項數。

式(1)中等式左邊時刻解吸量與極限解吸量的比值定義為解吸分數。在得到解吸曲線之后，采用漸進曲線法，依照曲線形態作解吸曲線的漸近線，取得極限解吸量，取時刻解吸量與之相除，最終得到該特定曲線的解吸分數。之后通過利用式(1)進行擬合，便可以得到不同損失時間條件下的表觀擴散系數有效值。從圖4和表2的擬合結果可以看出，試驗煤樣的表觀擴散系數有效值在10～10s。而如果假設擴散距離為煤粒本身半徑(0.1 mm左右)，則求得的有效擴散系數在10～10m/s，和以往文獻中試驗測得的有效擴散系數數量級相近。而對于單孔模型的擬合度，會發現采用單孔模型擬合解吸分數曲線時擬合度不一。在損失時間為1 min時的擬合度最高，隨著損失時間向兩邊偏離于1 min時，越偏離擬合度越差。這主要與單孔模型本身濃度恒定的假設有關。在解吸初期游離瓦斯的強烈流動影響了表觀擴散系數有效值的大小，使得表觀濃度波動過大；而在解吸后期，隨著濃度急劇下降，且基質產生解吸收縮，孔隙結構發生了顯著改變，這又與單孔模型的均質假設相悖，故形成了1 min時擬合度最高的結果。

表2 表觀擴散系數有效值擬合結果

圖4 單孔擴散模型擬合曲線

將擬合得到的表觀擴散系數有效值與損失時間對應起來，形成如圖5的關系曲線。煤樣瓦斯表觀擴散系數有效值隨著損失時間的增加而逐漸降低。此種變化規律與以往文獻中瞬態擴散系數的變化規律相似，其本質是由于初期瓦斯解吸時存在游離態瓦斯，游離態瓦斯高速逸散時，產生的隨體流動加大了整體表觀擴散系數有效值，導致試驗瓦斯解吸量比實際煤樣上吸附的瓦斯解吸量更大，形成了更為劇烈的傳質效果。

圖5 不同損失時間下的表觀擴散系數有效值變化

3 討論

3.1 類瞬態擴散衰減模型的不同損失時間下表觀擴散系數衰減模型

參考瞬態擴散系數的建模思路，采用類偽一級動力學模型對不同平衡壓力條件下不同損失時間的表觀擴散系數有效值進行擬合(式(2))。偽一級動力學方程是在反應動力學基礎上建立的，其以指數函數為特征，建立了反應時間和反應質量的關系，適合描述大多數動力衰減過程，包括氣體吸附傳質過程。LOSKUTOV等曾指出自擴散系數與衰減時間存在指數函數的關系，而自擴散系數與表觀擴散系數在低壓段可近似為等價。因此，參照上述思路可建立下述衰減關系：

=exp(-)+

(2)

式中，為不同時刻的表觀擴散系數有效值，s；為最終時刻的表觀擴散系數有效值，s；，為待定系數。

為了更好地表示損失時間對表觀擴散系數有效值的影響，定義損失時間修正因子，其為每個損失時間下的表觀擴散系數有效值與最終時刻的表觀擴散系數有效值的比值，即

(3)

則不同損失時間的修正因子與最終時刻的表觀擴散系數有效值的乘積即對應損失時間的表觀擴散系數有效值，即

(4)

采用式(2)對圖5中的數據進行擬合，并根據式(3)，得出待定系數與修正因子(表3)。需要指出的是，由于選取點數量及測試方法的限制，使得各系數的變化規律不明顯，這對解釋待定系數的物理意義有所干擾，但仍能定性地說明一些問題。當=0時，初始表觀擴散系數有效值為+1，則決定了表觀擴散系數有效值的初始大小；而則決定了表觀擴散系數有效值衰減的快慢。

表3 不同壓力下最終時刻表觀擴散系數有效值擬合結果

3.2 游離瓦斯流動與吸附瓦斯擴散的競爭關系

怎樣厘定游離瓦斯在整個解吸過程中的比例是確定基質吸附態瓦斯表觀擴散系數精準與否的關鍵，而怎樣描述裂隙中的游離瓦斯流動與基質中的孔隙擴散則又是其中的重中之重。為了簡化建模難度，煤中瓦斯解吸的過程可以假設為原有游離瓦斯先行流動解吸，之后吸附態瓦斯逐漸擴散形成新游離態瓦斯流的2個獨立關聯的過程，此種假設與滲透率模型中常見的雙重孔隙介質假設相似。根據滲透率模型中常見的平板模型，可以假設裂隙中的游離瓦斯是在對于截面積為、長度為的平板孔隙中運移(圖6)，則根據零級動力學方程，有

圖6 游離瓦斯運移示意

(5)

式中，為時刻解吸出的瓦斯總質量，g；為零級動力學系數。

當解吸開始時，解吸面將從右向左運移，由于解吸速度快，可假設瓦斯相和非瓦斯相質量濃度在各自區域內均一且與時間無關。則對于任意時刻，擴散面運移至距出口處時，非瓦斯相中逸出的瓦斯總質量為

=(-)

(6)

式中，為瓦斯相中的質量濃度，g/mL。

將式(6)代入式(5)，有

(7)

對于任意的時刻，積分可得

(8)

所以有

(9)

令=，則解吸分數為

(10)

式中，為系數，=()。

上述推導表明，游離態瓦斯初期解吸速度快，其解吸質量可近似看成與時間呈線性關系。這與張向陽等關于煤樣瓦斯解吸量隨解吸時間變化關系的試驗結果相符：不同壓力下鉆屑在開始解吸前的一段時間內解吸速度基本上趨勢一致。此現象的出現可能是由于解吸初期游離瓦斯占主導地位，解吸速率在該時間段內保持不變。因此，筆者嘗試將單孔擴散模型與游離瓦斯線性流動模型有機結合，形成引入游離瓦斯組分貢獻的瓦斯全濃度表觀擴散模型，即

(11)

式中，為游離態瓦斯對宏觀解吸的貢獻；1-為吸附態瓦斯對宏觀解吸的貢獻；,為煤粒時刻解吸出的游離態瓦斯和吸附態瓦斯質量；,分別為煤粒最終解吸出的游離態瓦斯和吸附態瓦斯質量。

應當指出的是，式(11)的轉換關系雖在形式上符合基質裂隙并聯傳質的關系，但RUCKENSTEIN等在采用相似方法推導雙孔擴散模型時曾指出：在小孔擴散系數與大孔擴散系數比值極小的條件下，意味著此時大小孔2種系統的相互關聯很小，則并聯傳質關系與串聯傳質關系的差異不大。因此，式(11)仍可用于描述裂隙(大孔系統)和基質(小孔系統)中瓦斯串聯傳質的關系。將式(1)和式(10)分別代入式(11)，可得到

(12)

3.3 游離瓦斯流動與吸附瓦斯擴散在宏觀解吸中的控制作用演化

式(12)可用于擬合解吸過程中游離瓦斯與吸附瓦斯共存階段的瓦斯解吸特性。將之用于擬合圖3中的解吸數據，可得到表4和圖7的結果。

表4 式(12)對試驗數據的擬合結果

圖7 游離瓦斯與吸附瓦斯對宏觀解吸貢獻度變化

從圖7表4中可以看出，游離瓦斯流動對宏觀解吸的貢獻在20%以下，這恰好與通常意義上游離瓦斯所占全部瓦斯的比例相似。另外，隨著損失時間的逐步推移，游離瓦斯在宏觀解吸中的貢獻比例越來越小，且隨著平衡壓力的增加，降低幅度有所減小。當平衡壓力為0.5 MPa時，貢獻值從0.186降至0.032，降低了83%；當平衡壓力為1 MPa時，貢獻值從0.191降至0.045，降低了76%，降幅基本不變；當平衡壓力為1.5 MPa時，貢獻值從0.178降至0.126，降低了29%，降幅顯著縮減。此種衰減規律與游離瓦斯含量以及吸附瓦斯含量隨壓力增高的變化規律有關。由于游離瓦斯含量與壓力呈線性關系，而吸附瓦斯含量與壓力呈Langmuir關系，因此兩者與壓力的關系存在3個演化階段(圖8)：① 斜率關系突變點之前，吸附瓦斯含量要大于游離瓦斯含量，且吸附瓦斯含量增長速率也要大于游離瓦斯含量增長速率(>且′>′)；② 斜率關系突變點之后、大小關系突變點之前，吸附瓦斯含量要大于游離瓦斯含量，但吸附瓦斯含量增長速率小于游離瓦斯含量增長速率(>且′<′)；③ 大小關系突變點之后，吸附瓦斯含量要小于游離瓦斯含量，此時成為負吸附現象(<)。試驗結果中，0.5 MPa和1 MPa下的游離瓦斯貢獻降幅相似，說明斜率關系突變點在0.5～1.0 MPa附近，而之后在1.5 MPa時游離瓦斯貢獻降幅明顯縮減，說明游離瓦斯含量較多且逸散更為持久，使得相同損失時間條件下逸散出的游離瓦斯量占總游離瓦斯量的比例縮小，正好處于階段②之后。