火山渣顆粒圖像處理及形狀定量分析

聶志紅 廖靖云 周蘇華 李冰曉 安愛軍

摘 ? 要：火山渣顆粒具有形狀不規則、棱角較多的特點，其形狀特征會影響材料工程性質.為探究火山渣顆粒形狀特性，通過AIMS集料圖像測量系統獲得火山渣顆粒圖像以及形狀指標，采用分形維數和正態分布檢驗方法，討論以顆粒長度作為特征粒徑時粒徑變化對顆粒形狀和分形維數的影響規律，并對比火山渣顆粒與一般巖體顆粒的形狀差異.結果表明：隨著特征粒徑減小，火山渣顆粒長扁度和粗糙度逐漸減小，棱角度、球度和分形維數逐漸增大，越趨近于球形和塊狀，且表面棱角越少.所有采用的形狀指標之間均有較好的相關性，其中棱角度對特征粒徑變化最為敏感，可優先選取棱角度作為描述火山渣形狀的指標.此外，通過與卵石和級配碎石的對比表明，在不同粒徑范圍內，火山渣顆粒長扁度和圓度隨粒徑變化不如卵石和級配碎石顯著，且火山渣粗糙度和棱角度均比卵石和級配碎石大，表明火山渣顆粒棱角更多，表面更粗糙.

關鍵詞：火山渣;分形維數;長扁度;棱角度;球度;粗糙度

中圖分類號：TU411 ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼：A

Image Processing and Shape Quantitative Analysis

of Volcanic Cinder Particles

NIE Zhihong1，LIAO Jingyun1，ZHOU Suhua2？覮，LI Bingxiao3，AN Aijun4

（1. School of Civil Engineering，Central South University，Changsha 410075，China;

2. College of Civil Engineering，Hunan University，Changsha 410082，China;

3. School of Mathematics and Statistics，Central South University，Changsha 410075，China;

4. China Road and Bridge Corporation，Beijing 100000，China）

Abstract：The volcanic cinder has the characteristics of irregular shape and sharper angles，and its shape characteristics affect the material engineering properties. In order to explore the shape of volcanic cinder particles，the image and shape index of volcanic cinder were obtained by AIMS aggregate image measurement system. The fractal dimension and normal distribution method were used to discuss the influence of fractal dimension and shape as the particle size change on the particle shape and fractural dimension. The differences in shape between the volcanic cinder and general rock mass were also investigated. The results show that as the particle size decreases，the elongation and texture decrease，the angularity，sphericity and fractal dimension increase，the volcanic cinder approaches spherical and massive，and the surface angularity is reduced. All the shape indicators have a good correlation，and the angularity is the most sensitive to the change of the particle size. Thus，it can be preferentially selected to describe the shape of the volcanic cinder. In addition，comparison with pebbles and crushed stone shows that the change of volcanic cinder elongation and texture is not as significant as that of the pebbles and crushed stone under different particle size. Moreover，the texture and angularity of volcanic cinder are larger than those of pebbles and crushed stone，indicating that the volcanic cinder has more angular edges and a rougher surface.

Key words：volcanic cinder;fractal dimension;elongation;angularity;sphericity;texture

火山渣是火山噴發中巖漿經過高溫燃燒噴發后冷卻形成的一種礦渣狀多孔輕質的顆粒材料，具有形狀不規則、棱角較多的特點.現有關于火山渣的研究多集中于火山渣填料物理力學性質方面[1-3]，對火山渣顆粒形狀特征研究較少.

許多研究證明顆粒形狀顯著影響材料工程特性，如最大孔隙比和最小孔隙比[4，5]，剪切強度[6，7]，剛度[8]以及剪脹性[9]. 為揭示顆粒形狀對于材料工程性質的影響規律，首先需研究如何對顆粒形狀進行評價.目前顆粒形狀定性評價主要通過目測對比完成[10]，該方法受主觀和經驗性因素影響較大. 近年來，圖像分析也廣泛應用于顆粒形狀測量中[11]. 陸厚根等[12]采用顯微鏡獲取粉煤灰顆粒圖像，并分析圓形度和微觀形態指標隨粒徑的變化;劉清秉等[13]通過顯微鏡獲取砂土顆粒圖像，測算了長扁度、球形度和棱角度等;張家發等[14]通過數字測量儀獲得碎石顆粒的長寬比、扁平度、球度等指標，并取多個截面的形狀指標平均值作為計算指標;曾濤等[15]采用光學顯微鏡獲得冰漬土的二維平面投影，并對顆粒軸向系數和棱角性進行統計分析;秦雪[16]采用相機獲得碎石顆粒圖像，對顆粒扁平度、球形度、棱角度等形狀指標進行分析. 此外，針對表面形狀輪廓復雜的顆粒，研究者也引用分形理論進行分析[17-19]，陳江峰等[20]對碎屑顆粒圓度進行分形描述;陳海洋等[21]應用分形理論證明鈣質砂的顆粒形狀具有分形特性.

現有研究取得了許多有價值的成果，然而，這些研究都屬于對單個圖像截面進行二值化處理，不能獲得顆粒微觀紋理，無法反應顆粒表面粗糙程度，在分析表面棱角較多的火山渣顆粒時具有局限性. 而AIMS集料圖像測量系統可以自動測量顆粒尺寸、棱角、微觀紋理、顆粒形狀等特征. 因此本文通過AIMS集料圖像測量系統獲得火山渣顆粒的二維圖像，計算出形狀指標，并采用分形理論和正態分布檢驗方法，討論粒徑變化對顆粒形狀指標和分形維數的影響規律，并探究火山渣顆粒與一般巖體顆粒的形狀差異.

1 ? 試驗方案

1.1 ? 試驗材料

研究對象取自肯尼亞內馬鐵路DK103處火山渣料場，現場火山渣采用顎式破碎機進行加工.試驗前將火山渣顆粒表面的塵土洗凈后進行烘干篩分.篩分結果如圖1所示. 由圖可知，火山渣試樣最大粒徑不大于75 mm，顆粒主要集中在5～60 mm，小于5 mm顆粒含量僅為0.6%.

1.2 ? 試驗儀器

采用美國PINE公司生產的AFA2A 型號AIMS集料圖像測量系統進行顆粒形狀采集，試驗儀器如圖2.該儀器可通過掃描儀自動采集顆粒的尺寸、棱角、微觀紋理、顆粒形狀等特征，并將其性能通過各種特征指標進行客觀描述.

1.3 ? 顆粒形狀參數的選取

通過AIMS集料圖像測量系統能準確獲取單個顆粒的幾何尺寸測值，如面積、周長、長度、寬度、等效直徑等. 表1給出了各類基本尺寸測值.

2 ? 特征粒徑選取

顆粒特征粒徑的選取是顆粒形狀定量分析方法的重要基礎. 因此，在進行形狀分析前應先確定特征粒徑.形狀分析中常用的特征粒徑有dL、dI、d、R. 將火山渣分別按以上四種粒徑計算試樣在篩孔上的通過率，并與實際篩分試驗所得級配曲線進行對比，如圖3所示. 由圖可知，dL對應的通過率與篩分試驗所得級配曲線更接近.因此，研究采用dL作為顆粒特征粒徑對火山渣顆粒進行粒組分類和形狀分析.

3 ? 單顆粒形狀指標分析

火山渣試樣中dL小于5 mm及大于60 mm的顆粒所占比例較小，故只選取dL范圍在5～10 mm，10～20 mm，20～40 mm和40～60 mm的4個粒徑組進行分析.每個粒徑組分別取100個顆粒，共400個顆粒進行測試和分析. 不同粒徑組下火山渣顆粒形狀如圖4（a）～（d）所示.

3.1 ? 顆粒分形特征分析

當周長的對數與面積的對數呈線性關系時，表明該顆粒的形狀具有良好的分形特征.運用上式（5）～（8）對四個粒徑區間上的100個顆粒進行分形維數曲線擬合，計算結果如下表2和圖5. 結果表明，四條曲線線性擬合較好，相關系數均在0.9以上，因此，火山渣顆粒具有良好的分形特征. 此外，分形維數還反映了顆粒周界的曲折程度，顆粒輪廓越是復雜多變，分形維數愈大. 由表2可知，粒徑越小，火山渣顆粒分形維數越大. 因此，粒徑越小，火山渣顆粒形狀越復雜，分形特征和形狀之間的自相似性更強.

3.2 ? 顆粒形狀統計分析

為探究不同粒徑下各形狀指標的分布規律，采用正態檢驗方法對火山渣顆粒的長扁度、棱角度、球度和粗糙度進行分析.火山渣顆粒長扁度分析結果如下圖6所示，統計結果見表3.結果表明，粒徑越大長扁度越大，且長扁度的值均大于0.7.李霜等[24]的研究表明長扁度等于0.4是作為區分顆粒形狀屬于針狀或塊狀的界限值，當長扁度大于0.4時，其形狀為塊狀，反之則為針狀. 因此火山渣顆粒基本為塊狀.

在正態檢驗方法中，峰度系數越接近3，說明其分布曲線越接近正態曲線;峰度系數越大說明分布曲線相對標準正態曲線越陡峭，形狀指標分布更集中.偏度系數大于0表示均值在大于峰值的一側，表現為右偏態（正偏態），偏度系數小于0表示均值在小于峰值的一側，表現為左偏態（負偏態）. 由表3可知，當粒徑大于10 mm時，偏度系數均大于0，頻數分布的峰值向左側偏移，表現為長尾向右側延伸的分布類型，因此，長扁度在粒徑大于10 mm時均為正偏態分布類型. 峰度系數接近于3，表明長扁度在不同粒徑下的分布均接近于正態分布.

火山渣顆粒的球度、棱角度、粗糙度在各個粒徑區間的正態分布情況如圖7～圖9所示，分析結果見下表4.由表可知，粒徑越大，棱角度和球度越小，粗糙度越大. 其中球度的峰度接近于3，分布曲線與正態分布曲線接近.不同粒徑下的火山渣棱角度分布曲線的偏度均大于0，為正偏態分布. 粒徑大于10 mm時，棱角度的峰度均大于3，說明在此粒徑下，棱角度分布曲線比標準正態分布曲線更陡峭.

變異系數可以反映形狀指標對不同粒組顆粒形狀差異的敏感性. 變異系數值越小，統計數據協調程度越高，均勻性越好. 由表4可知，棱角度的變異系數比球度和粗糙度大，表明棱角度對顆粒輪廓形狀的差異更敏感，更能表征不同粒徑下的顆粒形狀變化.因此在實際運用中，可優先選取棱角度作為描述火山渣顆粒形狀指標.

3.3 ? 形狀指標相關性分析

對火山渣顆粒的長扁度、球度、棱角度、粗糙度以及分形維數進行了相關性分析，結果如下表5所示.由表可知，各形狀指標間相關系數均大于0.8，可知火山渣顆粒形狀指標間的相關性較大.

同時，火山渣顆粒的粗糙度和棱角度、分形維數、球度、長扁度呈負相關.長扁度、棱角度、球度之間相互呈正相關關系. 由此可知，隨著粒徑減小，火山渣顆粒越趨近于球形和塊狀，且表面棱角越少.

3.4 ? 與一般巖體顆粒形狀對比分析

為分析火山渣顆粒與一般巖體顆粒形狀的區別，采集了洗凈后的卵石和級配碎石的圖像;取dL的范圍在5～10 mm，10～20 mm和20～40 mm的3個粒徑組，每個粒徑組分別取100個顆粒與火山渣顆粒進行對比. 三種顆粒的形狀對比圖如圖10所示.

對比不同顆粒在粒徑組中的形狀指標均值，結果如下圖11（a）～（d）所示. 由圖可知，在材料粒徑為5 ～ 10 mm時，卵石的長扁度和圓度最大，級配碎石的長扁度和圓度最小，而火山渣居中;實際表現為碎石更加扁平細長，卵石呈類橢圓形. 在粒徑為20 ～ 40 mm時，火山渣的球度最大，其長扁度僅低于級配碎石.粒徑為5 ～ 40 mm時，火山渣的粗糙度和棱角度均比卵石和級配碎石的大，因此火山渣棱角更多，表面更粗糙. 此外，火山渣的長扁度和圓度隨粒徑變化不如卵石和級配碎石顯著，表明火山渣形狀分布更加均勻.

4 ? 結 論

1）火山渣顆粒具有良好的形狀分形特征，火山渣顆粒粒徑越小，分形維數越大.

2）火山渣顆粒特征粒徑越大，顆粒長扁度和粗糙度越大，而棱角度、球度和分形維數越小.各形狀指標間線性相關性較強，其中粗糙度和球度、棱角度、分形維數呈負相關，其他形狀指標之間均呈正相關;

3）火山渣顆粒在不同粒徑下的棱角度均服從正偏態分布;不同粒徑下棱角度的變化比粗糙度和球度更大.因此，在實際運用中，可優先選取棱角度作為描述火山渣顆粒形狀的指標;

4）當粒徑為5～40 mm時，火山渣顆粒的長扁度和圓度隨粒徑變化不如卵石和級配碎石顯著.

參考文獻

[1] ? ?胡江洋，毛君，張浩，等. 火山渣作為路基填料的試驗研究 [J]. 深圳大學學報（理工版），2015，32（6）：625—631.

HU J Y，MAO J，ZHANG H，et al. Research on test characteristic of cinder and subgrade filling[J]. Journal of Shenzhen University（Science and Engineering），2015，32（6）：625—631. （In Chinese）

[2] ? ?蔣明鏡，鄭敏，王闖，等. 不同顆粒級配的某火山灰的力學性質試驗研究 [J]. 巖土力學，2009，30（s2）：58—61.

JIANG M J，ZHENG M，WANG C，et al. Different grain size distribution of experimental study on the mechanical properties of some ash[J]. Rock and Soil Mechanics，2009，30（S2）：58—61. （In Chinese）

[3] ? ?楊福珍，黃清友，張建棟. 火山渣在路基工程中的應用 [J]. 內蒙古公路與運輸，2006，（4）：24—27.

YANG F Z，HUANG Q Y，ZHANG J D. Application of volcanic cinder in subgrade engineering [J]. Highways & Transportation in Inner Mongolia，2006，（4）：24—27. （In Chinese）

[4] ? ?OH T M，JOO G W，YUN K J，et al. Crush characteristics of abrasive particles during abrasive waterjet rock cutting [M]. Boca Raton：Crc Press-Taylor & Francis Group，2015：1577—1581.

[5] ? ?PAYAN M，KHOSHGHALB A，SENETAKIS K，et al. Effect of particle shape and validity of G（max） models for sand：A critical review and a new expression [J]. Computers and Geotechnics，2016，72：28—41.

[6] ? ?ROUSE P C，FANNIN R J，SHUTTLE D A. Influence of roundness on the void ratio and strength of uniform sand [J]. Geotechnique，2008，58（8）：227—231.

[7] ? ?孔亮，彭仁. 顆粒形狀對類砂土力學性質影響的顆粒流模擬 [J]. 巖石力學與工程學報，2011，30（10）：2112—2119.

KONG L，PENG R. Particle flow simulation of influence of particle shape on mechanical properties of quasi-sands [J]. Chinese Journal of Rock Mechanics and Engineering，2011，30（10）：2112—2119. （In Chinese）

[8] ? ?KANDASAMI R K，MURTHY T G. Effect of particle shape on the mechanical response of a granular ensemble [M]. Boca Raton：Crc Press-Taylor & Francis Group，2015：1093—1098.

[9] ? ?KANDASAMI R K，MURTHY T G. Manifestation of particle morphology on the mechanical behaviour of granular ensembles [J]. Granular Matter，2017，19（2）：1—13.

[10] ?KRUMBEIN W C. Measurement and geological significance of shape and roundness of sedimentary particles [J]. Journal of Sedimentary Research，1941，11（2）：64—72.

[11] ?趙明華，戴杰，張玲，等. 基于分形理論的粉煤灰滲透率研究 [J]. 湖南大學學報（自然科學版），2015，42（1）：75—80.

ZHAO M H，DAI J，ZHANG L，et al. Fractal theory-based study of the permeability of fly ash [J]. Journal of Hunan University （Natural Sciences），2015，42（1）：75—80. （In Chinese）

[12] ?陸厚根，馬魁. 用兩個形狀指數表征粉煤灰顆粒形貌的研究 [J]. 硅酸鹽學報，1992，20（4）：293—301.

LU H G，MA K. Characterization of configuration of flyash particles by use of a pair of shape indexes [J]. Journal of the Chinese Ceramic Society，1992，20（4）：293—301. （In Chinese）

[13] ?劉清秉，項偉，M.BUDHU，等. 砂土顆粒形狀量化及其對力學指標的影響分析 [J]. 巖土力學，2011，32（s1）：190—197.

LIU Q B，XIANG W，BUDHU M，et al. Study of particle shape quantification and effect on mechanical property of sand [J]. Rock and Soil Mechanics，2011，32（s1）：190—197. （In Chinese）

[14] ?張家發，葉加兵，陳勁松，等. 碎石顆粒形狀測量與評定的初步研究 [J]. 巖土力學，2016，37（2）：343—349.

ZHANG J F，YE J B，CHEN J S，et al. A preliminary study of measurement and evaluation of breakstone grain shape [J]. Rock and Soil Mechanics，2016，37（2）：343—349. （In Chinese）

[15] ?曾濤，蔣良濰，羅強，等. 冰磧土和冰水堆積物填料顆粒形狀特征分析 [J]. 鐵道標準設計，2018，62（12）：45—51.

ZENG T，JIANG L W，LUO Q，et al. Particle shape characteristic analysis of glacial till and fluvioglacial [J]. Deposits Railway Standard Design，2018，62（12）：45—51. （In Chinese）

[16] ?秦雪. 基于圖像分析的集料形狀特征與分布特征研究 [D]. 青島：青島理工大學土木工程學院，2016：21—54.

QIN X. Reserch on shape characteristcs and distribution characteristcs of aggregate based on image analysis[D]. Qingdao：College of Civil Engneering，Qingdao University of Technological，2016：21—54. （In Chinese）

[17] ?向貴府，許模，崔杰. 昔格達地層巖土顆粒組成及分形特征研究 [J]. 重慶交通大學學報（自然科學版），2018，37（7）：60—75.

XIANG G F，XU M，CUI J. Study on composition and fractal characteristics of geological particles in xigeda formation [J]. Journal of Chongqing Jiaotong University （Natural Science），2018，37（7）：60—75. （In Chinese）

[18] ?王國梁，周生路，趙其國. 土壤顆粒的體積分形維數及其在土地利用中的應用 [J]. 土壤學報，2005，42（4）：19—24.

WANG G L，ZHOU S L，ZHAO Q G. Volume fractal dimension of soil particles and its applications to land use [J]. acta pedologica sinica，2005，42（4）：19—24. （In Chinese）

[19] ?劉曉明，劉凱，梁新權，等. 基于分形理論的紅層軟巖填料縮尺料制備方法 [J]. 湖南大學學報（自然科學版），2018，45（11）：127—132.

LIU X M，LIU K，LIANG X Q，et al. Scale sample preparation method of red bed weak rock filling material based on fractal theory [J]. Journal of Hunan University （Natural Sciences），2018，45（11）：127—132. （In Chinese）

[20] ?陳江峰，王振芬，閆純忠. 碎屑顆粒圓度的分形描述 [J]. 煤田地質與勘探，2002，30（4）：16—17.

CHEN J F，WANG Z F，YAN C Z. Fractal description of roundness of clastic particles [J]. Coal Geology & Exploration，2002，30（4）：16—17. （In Chinese）

[21] ?陳海洋，汪稔，李建國，等. 鈣質砂顆粒的形狀分析 [J]. 巖土力學，2005，26（9）：1389—1392.

CHEN H Y，WANG N，LI J G，et al. Grain shape analysis of calcareous soil [J]. Rock and Soil Mechanics，2005，26（9）：1389—1392. （In Chinese）

[22] ?COMPANY P I. Aggregate image measurement system operation manual [M]. Grove：COMPANY P I，2014：4—9.

[23] ?MANDELBROT B B. The fractal geometry of nature [M]. New York：WH freeman，1982：172—176.

[24] ?李霜，李煥英，李啟令. 粉煤灰顆粒群幾何特征的圖象分析研究 [J]. 建筑材料學報，2001，4（1）：49—54.

LI S，LI H Y，LI Q L. Research on the image analysis of the geometric characteristic parameters of fly ash[J]. Journal of Building Materials，2001，4（1）：49—54. （In Chinese）