基于CT掃描技術對煤巖的孔隙結構的提取和研究

邢 敏，吳金隨，張辭源，李更川，高 嵩

(1.華北科技學院 安全工程學院，北京 東燕郊 065201；2.準格爾旗云飛礦業有限責任公司 串草圪旦煤礦，內蒙古 鄂爾多斯 010400)

0 引言

煤炭資源是我國現階段最主要的原料和能源，是我國現階段必不可少的經濟來源，促進著我國經濟的快速增長[1]。盡管現在有很多新能源的利用，但煤炭資源在未來的很長時間內在供給人類生活所需和促進經濟發展中仍然處于主要地位，因此煤礦企業是促進國民經濟的重要的力量之一，在經濟發展中具有重要的戰略作用。要想合理的使用煤炭資源，就需要在對煤礦中的煤層進行深入的研究。

煤巖作為煤層中的重要組成部分，其結構復雜多變，類型多種多樣。孔隙是煤巖中的重要結構，通過研究煤的孔隙特征，可以了解煤的很多特性。為了更清楚地了解孔隙的相關特性，現階段可以采用壓汞法、核磁共振法等方法直接測出樣品參數[2-5]，除此之外，為了直觀地看到孔隙模型，對孔隙進行可視化分析，還可以采用三維數字巖心技術利用CT提取出樣品數據，運用計算機完成重構，將重構數據與可視化軟件結合，建立孔隙模型，測量出孔隙參數，然后就可以利用相關參數研究孔隙的特性[6-9]。鄧遠剛等運用CT掃描技術對典型脆性巖樣單軸壓縮破壞后破裂面及全斷面進行掃描，獲取全斷面的細觀形態，分析和總結了典型脆性巖樣破壞細觀裂紋擴展的發展規律[10]；查明等利用CT掃描技術從三維角度表征儲層微觀孔隙結構、裂縫分布、微觀流體特征等，定量分析巖石物性參數，動態地監測微觀流體驅替、巖石裂縫的延伸規律[11]；黃宛瑩等利用三維動態CT圖像數據，通過圖像切割、圖像分割等處理步驟，以團簇為基本表征單位，表征不同形狀、大小孔隙的動態變化過程，分析對比全部孔隙、大團簇、小團簇和裂隙型小團簇的動態圖像[12]。

本次實驗首先使用CT掃描儀對采集的煤樣樣品進行掃描，經重構后得到相關數據；然后將數據導入三維可視化軟件Avizo中，在Avizo中選取合適的交互式閾值，閾值分割后使用Axis Connectivity、Arithmetic等運算模塊提取孔隙結構模型和球棍模型；最后，在各個模型的基礎上計算出孔隙度、連通孔隙度、孤立孔隙度、配位數等參數，由此展開對所采集的煤樣孔隙的分析。

1 煤樣數據采集

本文選擇在準格爾旗云飛礦業有限公司串草圪旦煤礦不同區域中采集的4個煤樣作為本次實驗樣品，為便于分析，分別標號為coal-1、coal-2、coal-3、coal-4。

煤樣采集完成后，采用天津三英科技公司生產的Nano Voxel 3000高分辨率臺式CT掃描儀對所要分析的煤樣進行掃描，按照操作步驟掃描完成后，進行數據重構，重構后得到相關數據和圖像。本文以coal-1為例展示孔隙結構的圖像，重構后得到的煤樣切片和煤樣三維渲染圖像如圖1所示。

圖1 煤樣切片及三維渲染圖像

2 煤樣孔隙結構的提取

為了提取煤樣中的有效孔隙結構，在將重構完成的煤樣數據導入Avizo之后，首先要對樣品進行剪裁，一方面，由于樣品數據比較大，為了便于Avizo軟件的運行與對煤樣孔隙結構的分析，需要對其進行剪裁；另一方面，通過剪裁可以去掉樣品周圍的空氣部分并剪切成規則的形狀，從而確保不會因含有無效空隙而影響煤樣的孔隙結構提取，致使樣品的孔隙度受到影響。

在使用Avizo進行孔隙參數的提取前，首先要分離出孔隙。數據體剪裁后選擇中值濾波進行過濾，然后采用交互式閾值分割分離出孔隙。交互式閾值是整個分析過程中最重要的步驟之一，調整好交互式閾值是后期分析孔隙結構的各種參數的基礎。每個煤樣尺寸不同，切片渲染效果也不完全相同，因此對于每個樣品都需要通過對比其切片圖像調整交互式閾值，對于coal-1、coal-2、coal-3、coal-4等四個樣品分別將閾值調整為12411、21466、2560、19600進行分割。如圖2是對coal-1的煤樣切片及交互式閾值分割圖像。

圖2 煤樣切片及交互式閾值分割圖像

在對煤樣進行交互式閾值分割后，分離出本次研究所需要的孔隙，方便下一步從孔隙結構中分離出連通孔隙與孤立孔隙。本次實驗分別從兩個視角進行展示，渲染后的孔隙結構如圖3。

在分離出孔隙結構的基礎上，利用Avizo中的Axis Connectivity運算模塊提取出連通孔隙，渲染后的連通孔隙結構如圖4。連通孔隙是樣品中的是氣和水流通的孔隙結構部分，由孔喉結構組成。提取連通孔隙是進一步分析孔喉結構特性和煤樣滲透率的關鍵步驟。研究連通孔隙不論是在通過研究滲透率來幫助分析煤頂板突水的情況，還是對研究礦井內瓦斯流動從而進一步防治礦井瓦斯方面都有很重要的作用，因此研究連通孔隙的特性具有重要的意義。

下一步在連通孔隙的基礎上運用Avizo里的運算模塊計算出煤樣的球棍模型，通過調整孔喉的參數來展示球棍模型，可供選擇調整的參數有體積、面積和等效半徑。本次實驗選擇調整孔喉的等效半徑，調整前后的球棍模型如圖5。

孤立孔隙不與外界連通，利用Avizo中的Arithmetic數據運算模塊分離出孤立孔隙，渲染后的孤立孔隙結構如圖6。

圖3 煤樣孔隙結構

圖4 煤樣連通孔隙結構

圖5 球棍模型

圖6 煤樣孤立孔隙結構

為了更直觀的看出煤樣中連通孔隙和孤立孔隙各自的分布情況，再運用Avizo軟件中的數據運算模塊進行連通孔隙和孤立孔隙的渲染，將其放到一起展示，并且為了更好地觀察連通孔隙和孤立孔隙的結構分布，增加x-z、y-z視角的渲染圖像如圖7所示，其中紅色孔隙結構部分為連通孔隙，綠色孔隙結構部分為孤立孔隙。

圖7 煤樣孔隙多視角渲染圖像

3 煤樣孔隙參數的提取與分析

在運用Avizo提取出孔隙結構后，下一步進行孔隙參數的提取。

(1) 煤樣孔隙度及所占比例

在使用Avizo提取出孔隙、連通孔隙和孤立孔隙后，即可提取出與之對應的孔隙度以及連通孔隙和孤立孔隙在總孔隙中的占比情況，如表1是本次研究的四個煤樣的孔隙度及占比情況。

表1 煤樣孔隙度匯總表

通過以上分析所得的孔隙度發現，coal-1的連通孔隙度為11.49%，在四個煤樣中最大，且該煤樣的連通孔隙在總孔隙中占比也較高，這表明該煤樣的連通性較好；此外，coal-3的連通孔隙度為4.90%，在四個煤樣中最小，且該煤樣的連通孔隙在總孔隙中占比較低，因此該煤樣的連通性較差。此外，連通性越好的煤樣一般滲透性能也越好，因此，coal-1的滲透性也越好。

(2) 煤樣孔喉結構的參數提取與分析

提取孔喉結構可以分別得到孔喉的參數以及配位數，其中配位數即一個孔所連接的喉的個數。通過Avizo對四個煤樣提取孔喉結構，可得coal-1的連通孔隙中有12個孔和21個喉道，coal-2的連通孔隙中有5個孔和5個喉道，coal-3的連通孔隙中有11個孔和23個喉道，coal-4的連通孔隙中有6個孔和8個喉道。本次實驗所分析的四個煤樣的連通孔隙中孔的參數和配位數見表2～表5。

另外，為了更直觀地觀察連通孔隙中的孔隙大小與配位數的變化關系，繪制出每個煤樣的連通孔隙體積與對應的配位數變化圖像，如圖8～圖11。

表2 coal-1連通孔隙中孔的參數和配位數

表3 coal-2連通孔隙中孔的參數和配位數

表4 coal-3連通孔隙中孔的參數和配位數

表5 coal-4連通孔隙中孔的參數和配位數

圖8 coal-1連通孔隙體積變化及配位數變化圖

圖9 coal-2連通孔隙體積變化及配位數變化圖

圖10 coal-3連通孔隙體積變化及配位數變化圖

圖11 coal-4連通孔隙體積變化及配位數變化圖

從以上所得參數以及圖表中可以分析出：

四個煤樣中的孔和喉存在數量差異有兩部分原因，一部分原因可能是由于在分析前剪裁的樣品的體積不是完全相同的，致使其中連通孔隙的數量存在差異；另一部分原因由于每個樣品的孔隙特性不同導致孔喉個數的差異。

觀察以上四個煤樣的連通孔隙體積變化圖和對應配位數變化圖像，發現二者變化的趨勢大致相同。通過對比二者圖像發現孔的體積越大，對應的配位數一般越大，連接的喉越多，即體積越大的孔隙能夠連接到多個喉道的概率越大。在孔喉結構提取的過程中，對每個孔進行編號的作用是可以將煤樣中每個孔的體積與配位數一一對應，方便觀察孔的體積與配位數的變化關系。

4 結論

(1) 在使用Avizo提取孔隙結構時，選擇合適的交互式閾值分割，交互式閾值分割的界面有x-y、x-z、y-z三個界面的選擇，通過選擇發現，由于Avizo在分析導入的樣品數據時把其當作一個整體來看，所以不論選擇哪個界面作為主界面，都不影響提取出來的孔隙度、連通孔隙度和孤立孔隙度。只有能提取出連通孔隙的樣品才能提取出球棍模型及孔喉參數。

(2) 通過使用Avizo提取煤樣中的孔隙結構，再運用數據運算模塊提取孔隙度、連通孔隙度和孤立孔隙度，觀察連通孔隙在總孔隙中的占比情況，可以知道：在四個煤樣中coal-1的連通孔隙度最大，為11.49%，連通孔隙占比較高，連通性最好，因此滲透性也更好；coal-3的連通孔隙度最小，為4.90%，連通性最差。

(3) 在提取出孔喉參數后，發現四個樣品的孔喉數量存在差異的原因有兩部分，一是每個煤樣本身的孔隙特性不同，二是每個煤樣剪切后的體積不同。

(4) 對比每個煤樣的連通孔隙體積變化圖和對應的配位數變化圖像發現，煤樣中孔的體積越大，與之對應的配位數一般越大，能夠連接到多個喉道的概率就越大。

(5) 使用CT掃描技術和Avizo提取煤樣的孔隙結構，可以得到直觀的、立體的孔隙渲染圖像和球棍模型，還可以通過相應的運算模塊計算出孔隙的相關參數，有助于從微觀層面上分析煤巖的孔隙特性，幫助分析煤層的組成成分，為防治礦井災害提供理論支持，進一步提高煤礦領域的生產作業的安全性。