分形理論在煤巖物理形態(tài)表征中的應(yīng)用

王力力

(開灤(集團(tuán))唐山礦業(yè)分公司，河北省唐山市，063000)

?

分形理論在煤巖物理形態(tài)表征中的應(yīng)用

王力力

(開灤(集團(tuán))唐山礦業(yè)分公司，河北省唐山市，063000)

研究總結(jié)了國(guó)內(nèi)學(xué)者利用分形理論，對(duì)煤巖顆粒粒度、表面結(jié)構(gòu)、孔隙結(jié)構(gòu)所建立的分形模型及對(duì)煤巖形態(tài)的表征。研究結(jié)果表明，粒度分形維數(shù)反映了在煤巖解離中的破碎方式、煤巖組成、能量輸入對(duì)其粒度分布的影響，為研究煤巖顆粒選擇性解離、粒度分布及破磨過程中的能量消耗奠定了相關(guān)理論基礎(chǔ)；同時(shí)煤巖表面結(jié)構(gòu)的分形維數(shù)反映了其表面形態(tài)特征結(jié)構(gòu)的物理特性和表面形態(tài)，可以增強(qiáng)對(duì)煤巖熱解、氧化過程發(fā)生化學(xué)反應(yīng)的了解；孔隙結(jié)構(gòu)的分形維數(shù)反映了煤巖變質(zhì)程度、孔隙發(fā)育度以及環(huán)境因素對(duì)煤巖孔隙分布的影響，有助于更清楚地了解煤巖孔隙對(duì)其破碎、燃燒、熱解、瓦斯突出等方面的影響。

分形維數(shù) 顆粒粒度 表面結(jié)構(gòu) 孔隙結(jié)構(gòu)

分形理論主要是針對(duì)自然和社會(huì)等領(lǐng)域的不規(guī)則性、復(fù)雜性以及相似性進(jìn)行描述和研究，至今已在物質(zhì)結(jié)構(gòu)表征、粒度特征分布等具體領(lǐng)域得到了廣泛應(yīng)用。分形是分形理論中重要的概念，被用來(lái)描述大量不規(guī)則事物中蘊(yùn)含的規(guī)律，認(rèn)為事物組成部分都是以某種方式與整體呈現(xiàn)出自相似性。近年來(lái)，分形理論被引入到物理、化學(xué)、醫(yī)藥、生物、材料以及礦物加工等學(xué)科專業(yè)，成為科學(xué)界熱門的研究課題。

在煤炭應(yīng)用方面，為了促使煤炭資源高效化和潔凈化利用，如煤炭液化和稀缺煤資源解離再選等，都必須使煤巖組分達(dá)到分離和有機(jī)質(zhì)富集，而煤巖的物理形態(tài)是煤炭的最基本性質(zhì)，基于分形理論可以更好地認(rèn)識(shí)煤巖物理形態(tài)的內(nèi)在規(guī)律，為此國(guó)內(nèi)眾多學(xué)者探討了煤巖解離前后的粒度、形貌以及孔隙結(jié)構(gòu)等物理形態(tài)，為煤巖解離提供了有益的借鑒。基于此，為了促進(jìn)分形理論在表征煤巖物理形態(tài)方面的研究及進(jìn)展，本文總結(jié)了分形理論在煤巖顆粒粒度、顆粒表面結(jié)構(gòu)以及孔隙結(jié)構(gòu)表征方面的相關(guān)研究成果。

1 分形理論在粒度方面的研究

1.1 粒度分形模型建立

粒度分布特征是煤巖重要的物理性質(zhì)之一，與其形貌和孔隙結(jié)構(gòu)存在著直接關(guān)系。在描述顆粒粒度分布時(shí)，通常認(rèn)為顆粒粒度是連續(xù)分布的，但是實(shí)際上，顆粒粒度呈現(xiàn)離散型分布，這導(dǎo)致二者之間出現(xiàn)矛盾。而分形理論認(rèn)為，煤巖粒度具有分形特征，其粒度分布可用冪律分布的形式表示，見式(1)：

N(xk)=C0xk-D

(1)

式中：N(xk)——粒度大于xk的每塊數(shù)目，即篩孔為xk時(shí)篩上煤塊數(shù)目，個(gè)；

C0——比例常數(shù)；

xk——篩孔直徑，mm；

D——粒度分形維數(shù)，D值越大，表示碎塊越多、粒度越小、破碎越嚴(yán)重。

根據(jù)分形理論，對(duì)于由簡(jiǎn)單冪律分布表示的Gaudin-Schuhmann分布，設(shè)其直徑小于xk的顆粒累積質(zhì)量產(chǎn)率為F(xk)，見式(2)：

(2)

式中：F(xk)——直徑小于xk的顆粒累積質(zhì)量產(chǎn)率， %；

xmax——粒度最大煤塊粒度，mm；

m——分布模數(shù)，與物料性質(zhì)和設(shè)備性能有關(guān)，m=3-D。

即Gaudin-Schuhmann分布的分形表示式見式(3)：

(3)

分形維數(shù)D的求法是先求出分布模數(shù)m后，由m=3-D求出。為此，對(duì)式(2)兩邊求對(duì)數(shù)見式(4)：

lg(F(xk))=mlg(xk)+lgb

(4)

式中：lgb——常數(shù)。

分別以lg(xk)和lg(F(xk))為橫坐標(biāo)和縱坐標(biāo)，并擬合成直線，由直線斜率可求得分形維數(shù)D。

1.2 煤巖顆粒粒度表征

煤巖顆粒的粒度在礦物加工中有著重要的應(yīng)用，不同的煤加工利用方法對(duì)煤巖粒度有著特定的要求，同時(shí)也是判斷煤中礦物質(zhì)是否充分解離以及脫除有害雜質(zhì)的基礎(chǔ)。因此，對(duì)煤巖顆粒粒度合理的表征不僅有助于煤炭的高效和潔凈利用，也有助于研究煤的破碎機(jī)理和粉碎理論，有助于改善煤的破磨設(shè)備性能，以便降低能耗和更好地控制產(chǎn)品粒度。

在煤巖顆粒粒度表征方面，有專家基于分形理論的粉碎模型研究認(rèn)為，煤中有機(jī)質(zhì)與礦物質(zhì)的解離只與礦物顆粒在煤中的分形維數(shù)有關(guān)；另有專家的研究表明，煤炭經(jīng)過粉碎，粒度分布受到粒度分布維數(shù)和最大粒徑的影響，而分布維數(shù)又受到顆粒的粉碎概率以及粉碎相似比的影響，顆粒分形維數(shù)大，則顆粒群的粒級(jí)偏粗，否則顆粒群的粒級(jí)偏細(xì)。

煤巖的粒度分布也受到破碎方式的重要影響。有專家研究表明，分形模型能夠很好地描述煤塊沖擊破碎后的粒度分布特征，粒度分形維數(shù)可以反映沖擊破碎效果的恰當(dāng)統(tǒng)計(jì)特征量，以山東泰山能源有限責(zé)任公司協(xié)莊煤礦(以下簡(jiǎn)稱協(xié)莊礦)、鞏義大峪溝礦(以下簡(jiǎn)稱大峪溝礦)以及徐州礦務(wù)集團(tuán)有限公司夾河煤礦(以下簡(jiǎn)稱夾河礦)為例，煤塊破碎粒度分布維數(shù)及相關(guān)系數(shù)見表1，在控制破碎產(chǎn)品粒度以及判斷破碎效果等方面可作為重要參數(shù)。

表1 煤塊破碎粒度分布維數(shù)及相關(guān)系數(shù)

還有專家通過研究煤和矸石在沖擊條件下的破碎概率差異的分形行為，揭示了煤和矸石顆粒沖擊破碎概率差異的內(nèi)在機(jī)理是由于顆粒的來(lái)源不同以及二者之間的物理機(jī)械性能的差異；部分專家也研究了沖擊條件下煤的破碎行為，研究表明對(duì)塊度分形維數(shù)的影響由大到小依次為沖擊速度、物料硬度和沖擊次數(shù)，沖擊條件下煤的破碎行為試驗(yàn)結(jié)果分析見表2。

表2 沖擊條件下煤的破碎行為試驗(yàn)結(jié)果分析

由表2可以看出，煤樣的硬度越小，沖擊速度越大，則分形維數(shù)越大。另有專家通過對(duì)煤樣破碎能進(jìn)行實(shí)驗(yàn)研究發(fā)現(xiàn)，矸石顆粒的分形維數(shù)和破碎常數(shù)均小于煤顆粒的分形維數(shù)和破碎常數(shù)；在相同沖擊速度下，與矸石顆粒的沖擊破碎概率相比，煤顆粒的沖擊破碎概率較大，并且存在最優(yōu)沖擊速度使二者的沖擊破碎概率差值最大，對(duì)大峪溝礦和夾河礦煤巖沖擊破碎表明，其最優(yōu)破碎速度分別為11.47m/s和15.80m/s；研究不同截割條件下煤粒度的分布規(guī)律后發(fā)現(xiàn)，分形維數(shù)可以很好地辨別出不同截割條件下煤的破碎程度以及與各個(gè)參數(shù)(抗壓強(qiáng)度、切屑厚度)的關(guān)系。

在破碎過程中，煤巖顆粒粒度分布與其自身礦物質(zhì)組成也有關(guān)系。通過研究煤粉碎過程中顆粒形狀的分形特征發(fā)現(xiàn)，不同煤顯微組分組成的顆粒其顆粒形狀分維是不同的，礦物質(zhì)組成的單顆粒形狀分維最大，其次是殼質(zhì)組組成的顆粒。

能量輸入對(duì)煤巖粒度的分布也有重要的影響。在研磨煤炭過程中，為了達(dá)到動(dòng)態(tài)的預(yù)測(cè)和控制任意時(shí)刻的粒度分布，建立了顆粒粒度動(dòng)態(tài)變化的3種模型，即質(zhì)量分布模型、顆粒數(shù)分布模型、表面積分布模型。其中，在研磨過程中顆粒的堆密度和成漿密度可由質(zhì)量分布模型計(jì)算，粒徑主要存在范圍可由顆粒數(shù)分布模型進(jìn)行直觀描述，能耗變化可由表面積分布模型進(jìn)行表征，并且通過實(shí)際顆粒粒度分布與分形粒度分布模型相對(duì)比，驗(yàn)證了分形粒度分布模型的正確性。

2 分形理論在顆粒表面結(jié)構(gòu)方面的研究

2.1 顆粒表面結(jié)構(gòu)分形模型建立

在構(gòu)建顆粒表面結(jié)構(gòu)分形模型時(shí)，通常利用改變探測(cè)粒子大小、相關(guān)密度函數(shù)、根據(jù)測(cè)度關(guān)系這3種方法求分形維數(shù)。

改變測(cè)定粒子的大小求顆粒表面結(jié)構(gòu)的方法是指利用半徑大小為r的粒子來(lái)覆蓋所考察區(qū)域，所需粒子數(shù)計(jì)做N(r)，然后改變測(cè)定粒子大小r，若存在某一值Ds，滿足式(5)

N(r)∝r-Ds

(5)

式中：N(r)——探測(cè)粒子半徑為r時(shí)覆蓋考察區(qū)域所需粒子數(shù)，個(gè)；

r——探測(cè)粒子的半徑，mm；

Ds——顆粒表面的分維。

則該區(qū)域的表面分維為Ds，2≤Ds≤3，而該區(qū)域的表面積S見式(6)：

S=N0r2-Ds

(6)

式中：S——顆粒的表面積，其大小變化與探測(cè)粒子半徑大小r有關(guān)，mm2；

N0——煤表面大小和探測(cè)粒子大小無(wú)關(guān)的測(cè)度；

若Ds=2則表示為一光滑表面，若Ds=3則表示表面為一盤旋狀，所有的空間都會(huì)被充滿。

煤分形表面的實(shí)驗(yàn)技術(shù)方法有間接法和直接法，間接法包括吸附法(單分子層吸附法和多分子層吸附法)、小角度x射線衍射法、孔隙法；直接法包括掃描電子顯微鏡(SEM)、原子力顯微鏡(AFM)、CCD攝像儀等方法。對(duì)于直接法，通常需要利用適當(dāng)?shù)臄?shù)字圖像處理來(lái)研究煤的分形表面結(jié)構(gòu)。

2.2 煤巖顆粒表面結(jié)構(gòu)表征

煤的表面結(jié)構(gòu)特征不僅影響煤表面物理性質(zhì)和化學(xué)性質(zhì)，也影響煤的轉(zhuǎn)化和利用。因此，對(duì)煤巖顆粒表面結(jié)構(gòu)的表征，在煤燃燒過程模擬、潔凈煤技術(shù)和煤的有效利用方面顯得非常重要。

有研究表明，表面分形維數(shù)可以反映超細(xì)化煤粉表面形態(tài)分形特征，可以作為表征表面形態(tài)特征的重要參數(shù)，同時(shí)，表面分形維數(shù)越高代表煤粉顆粒表面結(jié)構(gòu)越復(fù)雜；相關(guān)專家利用SEM圖像分析了煤顆粒表面結(jié)構(gòu)，研究結(jié)果也表明分維分析可定量分析煤顆粒表面不規(guī)則狀況；另有專家對(duì)5種不同變質(zhì)程度煤樣的研究結(jié)果表明，中等程度變質(zhì)程度煤的分形維數(shù)低于高變質(zhì)程度煤和低變質(zhì)程度煤。

在煤巖顆粒表面與升溫氧化之間的關(guān)系方面，研究了原煤以及升溫氧化條件為50℃、100℃、150℃、200℃下煤表面介觀特征，經(jīng)試驗(yàn)分析證明，分形維數(shù)可以作為升溫氧化過程中，煤表面形態(tài)的特征參數(shù)，反映了煤表面結(jié)構(gòu)以及氧化作用，煤樣氧化溫度越高，分形維數(shù)越大。經(jīng)對(duì)淮南低變質(zhì)煤850℃燃燒過程中各階段煤顆粒SEM圖像分析表明，煤顆粒表面的分形維數(shù)有一個(gè)先變低后變高的過程，煤表面分維數(shù)在2.2～2.5區(qū)間內(nèi)變化。

3 分形理論在孔隙結(jié)構(gòu)方面的應(yīng)用

3.1 煤巖孔隙分形模型建立

孔隙的分維測(cè)量方法有離散法、散射法、密度法、吸附法、電子顯微術(shù)等，其中壓汞法最為常用。

壓汞法是多孔介質(zhì)孔隙測(cè)定的常用方法之一，測(cè)量時(shí)，將水銀注入尺寸為R的孔隙，但由于水銀與孔隙表面存在內(nèi)表面壓力，因此，需要施加壓力P，P與R之間滿足Washburn方程見式(7)：

P=-2σcosα/R

(7)

式中：P——注射壓力，N；

σ——汞的表面張力，N/m；

α——水銀與固體介質(zhì)的接觸角，(°)；

R——孔隙尺寸，m。

在給定壓力下，進(jìn)汞量Vp等于尺寸大于r的孔隙的總體積，因此，推導(dǎo)可得見式(8)：

(8)

式中：VP——顆?？紫扼w積，m3。

由分形理論Menger海綿模型可得剩余骨架體積Vs與孔隙半徑r的關(guān)系為：

式中：Vs——顆粒骨架體積，m3。

由式(7)、式(8)、式(9)和式(10)可推算出式(11)：

∝PD-4

(11)

對(duì)式(11)兩邊取對(duì)數(shù)可得式(12)：

(12)

分形維數(shù)可以通過式(12)關(guān)系求得，分形維數(shù)D=4+K，K是直線的斜率，由式(12)可求得顆?？紫兜姆中尉S數(shù)。

3.2 煤巖孔隙結(jié)構(gòu)表征

在煤巖孔隙與分形特征關(guān)聯(lián)方面，有專家從理論和實(shí)驗(yàn)兩方面證明了煤孔隙具有分形特征，并進(jìn)行了定量表征，研究了原生煤和構(gòu)造煤孔隙結(jié)構(gòu)進(jìn)行了特征分析。研究結(jié)果表明，原生煤和空隙煤的孔隙直徑分別在大于1000nm和205000nm時(shí)具有分形維數(shù)，可以利用分形維數(shù)對(duì)孔隙結(jié)構(gòu)進(jìn)行定量表征。原生煤和構(gòu)造煤孔隙分形維數(shù)情況如圖1所示。

由圖1可以看出，在煤孔隙分布與煤變質(zhì)程度關(guān)系的研究中，煤巖變質(zhì)程度越高，孔隙體積分形維數(shù)越??；同等變質(zhì)程度的煤，煤中吸附孔越多，孔隙分形維數(shù)越大。有研究表明石煤滲透分維數(shù)與變質(zhì)程度呈負(fù)相關(guān)關(guān)系，擴(kuò)散分維數(shù)與變質(zhì)程度呈正相關(guān)關(guān)系。

在煤巖孔隙與環(huán)境因素關(guān)系研究中，煤構(gòu)造變形強(qiáng)度越強(qiáng)烈，分形維數(shù)越高，孔表面積越大，微孔含量越多，孔隙結(jié)構(gòu)非均一性越強(qiáng)。由顆粒粒徑分別為4.00～1.70mm、1.70～0.38mm及0.38～0.18的煤顆粒制作成的型煤，隨著型煤顆粒粒徑逐漸減小，型煤中的孔隙半徑逐漸減小，孔隙總數(shù)逐漸增多，分形維數(shù)值逐漸增大，孔隙發(fā)育程度逐漸增大，孔隙分布均勻程度逐漸增大。0.38～0.18mm粒徑型煤表面孔隙分布分形維數(shù)計(jì)算如圖2所示。

由圖2可以看出，煤巖孔隙的分形特征與溫度也有密切關(guān)系，且隨溫度的升高越來(lái)越明顯，分形特征的孔隙尺寸的下界隨溫度升高而下降；同時(shí)，煤孔隙的分形維數(shù)與程度呈線性關(guān)系，在升溫的過程中，大孔隙擴(kuò)張速度慢于小孔隙擴(kuò)張速度。

4 結(jié)語(yǔ)

煤巖顆粒粒度分形維數(shù)反映了在煤巖解離中破碎方式、煤巖組成、能量輸入對(duì)其粒度分布的影響，為研究煤巖顆粒選擇性解離、粒度分布及破磨過程中能量消耗奠定了相關(guān)理論基礎(chǔ)。對(duì)煤巖表面結(jié)構(gòu)的研究，有助于更清楚了解煤巖表面結(jié)構(gòu)特性，以及其熱解、氧化等化學(xué)反應(yīng)。對(duì)煤巖孔隙結(jié)構(gòu)的分形研究，能夠更清楚了解煤中孔隙的分布和變化，以及建立煤巖的孔隙與破碎、燃燒、熱解、瓦斯突出等相關(guān)數(shù)學(xué)模型，用于指導(dǎo)工程實(shí)踐。

圖1 原生煤和構(gòu)造煤孔隙分形維數(shù)情況

圖2 0.38～0.18 mm粒徑型煤表面孔隙分布分形維數(shù)計(jì)算

[1] 史長(zhǎng)亮, 楊凱, 王海娟等. 細(xì)粒鋼渣干式磁選精礦品位分形特性研究[J]. 河南理工大學(xué)(自然科學(xué)版), 2014(5)

[2] 許滿貴,孟然,韓金子等. 軟煤體孔隙結(jié)構(gòu)影響瓦斯吸附特性的試驗(yàn)研究[J]. 中國(guó)煤炭,2015(10)

[3] 鄒明清. 分形理論的若干應(yīng)用[D]. 華中科技大學(xué), 2007

[4] 逯桂平. 分形理論對(duì)煤粉粒度和孔隙的表征[J]. 潔凈煤技術(shù), 2012(3)

[5] 曾凡桂. 煤的分形表面結(jié)構(gòu)[J]. 煤炭轉(zhuǎn)化, 1995(2)

[6] 王云鶴, 梁棟, 肖淑衡等. 升溫氧化過程煤表面介觀特征和分形參數(shù)研究[J]. 深圳大學(xué)學(xué)報(bào), 2006(4)

[7] 郭品坤,程遠(yuǎn)平,盧守青等. 基于分形維數(shù)的原生煤與構(gòu)造煤孔隙結(jié)構(gòu)特征分析[J]. 中國(guó)煤炭,2013(6)

[8] 林海飛,程博,曾強(qiáng)等.新疆阜康礦區(qū)煤的吸附孔分形特征試驗(yàn)研究[J]. 中國(guó)煤炭,2016(5)

[9] 張玉濤, 王德明, 仲曉星.煤孔隙分形特征及其隨溫度的變化規(guī)律[J]. 煤炭科學(xué)技術(shù), 2007(11)

[10] 王文峰, 徐磊, 傅雪海.應(yīng)用分形理論研究煤孔隙結(jié)構(gòu)[J]. 中國(guó)煤田地質(zhì), 2002(2)

[11] 高榮斌,閆明,張東營(yíng)等.構(gòu)造發(fā)育復(fù)雜性評(píng)價(jià)方法研究及應(yīng)用[J]. 中國(guó)煤炭,2013(6)

(責(zé)任編輯 王雅琴)

Application of fractal theory in the physical form characterization of coal rock

Wang Lili

(Tangshan Mining Branch Company of Kailuan (Group) Limited Liability Corporation, Tangshan, Hebei 063000, China)

The authors summarized the fractal models for coal rock particles size, surface structure and pore structure and their characterizations for coal rock form, which were studied by domestic scholars using fractal theory. The results showed that granularity fractal dimension reflected the effect of crushing modes, petrographic constituent of coal rock and energy input on particle size distribution during coal dissociation, which laid a relevant theoretical foundation for studying selective dissociation of coal rock, particle size distribution and energy consumption during crushing and grinding process; the fractal dimension of coal rock surface structure reflected physical property and surface morphology of coal rock, which was helpful to understand pyrolysis reactions and oxidation reaction of coal rock; the fractal dimension of pore structure reflected the effect of metamorphic grade, porosity development degree of coal rock and environmental factors on pore distribution, which was helpful to know clearly the effect of coal pore on crushing, burning, pyrolysis and gas outburst.

fractal dimension, particle size, surface structure, pore structure

王力力.分形理論在煤巖物理形態(tài)表征中的應(yīng)用[J].中國(guó)煤炭，2017，43(6):102-106.WangLili.Applicationoffractaltheoryinthephysicalformcharacterizationofcoalrock[J].ChinaCoal，2017，43(6):102-106.

TD

A

王力力(1967-)，女，河北唐山人，高級(jí)工程師，現(xiàn)任開灤(集團(tuán))唐山礦業(yè)分公司環(huán)保法律事務(wù)部主任，主要從事煤礦節(jié)能環(huán)保工作以及煤炭潔凈利用與加工方面的研究。