基于SEM圖像處理的天然硅藻土分形特征分析

王月香， 石修松， 吉 鋒， 顧歡達， 丁建文

(1. 蘇州科技大學 土木工程學院， 江蘇 蘇州 215011； 2. 南京工業(yè)大學， 江蘇 南京 211816； 3. 香港理工大學 土木與環(huán)境工程學院， 香港； 4. 東南大學 交通學院， 江蘇 南京 210018)

基于SEM圖像處理的天然硅藻土分形特征分析

王月香1, 2， 石修松3， 吉 鋒4， 顧歡達1， 丁建文4

(1. 蘇州科技大學 土木工程學院， 江蘇 蘇州 215011； 2. 南京工業(yè)大學， 江蘇 南京 211816； 3. 香港理工大學 土木與環(huán)境工程學院， 香港； 4. 東南大學 交通學院， 江蘇 南京 210018)

為探討天然沉積硅藻土的微觀結構及變形機制，對硅藻土樣進行了SEM電鏡掃描、壓汞試驗和高壓三軸試驗。采用圖像分割原理及最佳閾值法分析了天然沉積硅藻土在不同固結壓力下的分形特征和孔隙結構特征參數，探討了硅藻土分形維數和各向等壓固結壓力之間的關系。結果表明，硅藻土體多孔隙，具有很強的結構性。土中孔隙孔徑大多在0.1～1.0 μm之間。在固結壓力達到屈服壓力后，孔徑大于1 μm的較大孔隙隨固結壓力的增加而減少，硅藻顆粒破壞，孔隙結構坍塌。獲取了灰度閾值及二值化后的SEM圖像。得到天然沉積硅藻土在不同固結壓力下的分形盒子數和盒維數。結果還表明分形盒維數隨固結壓力增加而增加，天然沉積硅藻土在固結壓力下的微觀結構性及孔隙分布可通過分形盒維數來反映。

SEM圖像； 二值化； 硅藻土； 固結壓力； 分形特征

天然沉積的軟黏土一般都具有一定的結構性[1]，包括沉積過程中形成的初始結構性和沉積完成后，由于機械、化學和生物等各種因素影響而產生的結構性[2-3]。由于結構性的存在，天然沉積土的壓縮性與重塑土有著顯著差別。結構性對天然沉積飽和土力學形狀的影響與應力歷史和初始孔隙率的影響同等重要[1]，大量工程實踐表明，土的結構性與其工程性狀密切相關，近年來天然沉積土的結構性及其對土體力學性質影響的研究受到了廣泛關注。

目前關于結構性土力學特性的研究主要集中在宏觀和微觀兩個方面，其中宏觀的力學性質是微觀結構的表現(xiàn)，如何建立兩者之間的橋梁，是一個重要的科學問題[4-6]。隨著土微觀結構研究的深入，許多學者認為從微觀上認識土的結構是結構性定量化研究的重要途徑，微觀結構性研究方法主要是從微觀角度出發(fā)，通過室內壓汞試驗和掃描電鏡試驗(SEM)等對其微觀結構特征及其演變進行定性或定量的研究。王清等[7]通過SEM圖像處理技術研究了黏性土微結構中結構單元體的形態(tài)、定向性和孔隙特征等定量評價指標。洪振舜等[8-9]通過兩組結構性土的掃描電鏡試驗、壓汞試驗結合室內固結和高壓三軸試驗從宏觀和微觀兩個角度對結構性土的屈服特性進行了深入研究，他認為當固結壓力超過結構性土的結構屈服壓力時，土的結構性完全喪失，表現(xiàn)出與重塑土相同的力學性質。陳嘉歐[10]研究了珠江三角洲黏土的微觀結構特征，從黏性土的微觀結構與工程加固效果進行了初步探討。雷華陽等[11]探討了海積黏土微結構參數與力學強度的關系，隨后又從微觀角度出發(fā)解釋了固結曲線上拐點的意義[12]，進而研究了次固結系數的變化規(guī)律。盧佩霞等[13]用分形理論研究分析了非飽和土的孔隙結構，并建立了其微觀結構和宏觀力學性質之間的聯(lián)系。張先偉等[14]根據SEM圖像實現(xiàn)了軟土微觀結構的三維可視化并得到了真實分形維數。王升福等[15]運用分形理論研究了寧波凍融軟黏土的微觀孔隙特征，并采用容量維數分析了孔體積和孔表面積。以上成果不僅從微觀上描述了結構性土的微觀特征，而且建立了描述土體微觀結構的定量評價指標，這些研究成果對結構性土微宏觀特性的進一步研究和結構性土本構模型的建立具有重要的參考價值。

為了對天然沉積硅藻土在不同應力水平下的微觀結構進行定量評價，本文在SEM圖像的基礎上，結合前人的研究成果，探討了二值化圖像在不同固結壓力下的分形特征，以及分形維數與固結壓力的關系。

1 天然硅藻土的SEM成像特征

本文土樣在我國浙江省嵊州市以西約6 km硅藻土礦坑中取得，原上覆土層厚約1.8～2.4 m，采用刀具切削而成。由于其剛度非常大，結構性很強，因而刀具和挖土設備對其擾動很小。其主要參數見表1。Yoshinaka[16]和Tateishi[17]的研究成果也表明，硅藻土具有很強的結構性，它是一種高含水量的土，但是卻具有很高的強度和剛度。

表1 硅藻土的主要物理力學參數Tab.1 Index properties of diatomaceous earth

對硅藻土采用高壓三軸試驗進行各向等壓固結，然后進行SEM電鏡掃描試驗。在進行SEM試樣制備時，先將各向等壓固結后土樣切削成大小約為1 cm3的小塊，放入液氮中迅速冷凍1 h，然后在零下50 ℃抽真空干燥，并采用掰斷方式形成斷口拍攝面，所選取的圖像具有較為明顯的分形特征，因而具有較強的代表性。圖1為在不同固結壓力下，天然沉積硅藻土不擾動樣的電鏡掃描圖像，其中黑色區(qū)域為孔隙，白色區(qū)域為硅藻顆粒和骨架。由圖1(a)可知，天然硅藻土主要由圓形的硅藻顆粒組成，它們通過外側的環(huán)狀物連接在一起而具有很強的結構性。由于這種結構特征，在硅藻顆粒之間和顆粒內部具有很多大小不等的孔隙。HONG等[9]通過壓汞試驗表明，硅藻土在屈服前孔隙的主要分布區(qū)間為0.007～2.0 μm，對應的分布密度存在兩個峰值，分別為0.08和1.40 μm，其中0.08 μm主要是由存在于硅藻顆粒內部的孔隙形成，1.4 μm主要是由硅藻顆粒之間的大孔隙骨架導致。由圖1可知當固結壓力超過結構屈服壓力時，硅藻顆粒破壞，此后隨著固結壓力增加，大孔隙結構發(fā)生坍塌。

圖1 天然硅藻土在不同固結壓力下的SEM圖像Fig.1 SEM images of natural sedimentary diatomaceous earth under different confining pressures

表2為采用壓汞試驗所得各固結壓力下硅藻土的孔隙結構特征參數。王清等[7]根據分形理論，確定黃土的微、小、大孔隙界限，依據 0.02和0.800 μm兩個孔徑劃分。雷華陽等[7]提出結構性軟土在交通荷載作用下的孔徑劃分界限為 0.01，0.50，2.50 μm。根據前人研究成果，本文結合圖 1，將硅藻土孔隙分為：大孔隙，孔徑d>10.0 μm；中孔隙，1.0 μm

在達到屈服壓力后，孔徑大于1.0 μm的較大孔隙隨固結壓力的增加而逐漸減少。這一點與SEM圖像顯示的情況和前述的分析較一致。

表2 不同固結壓力下硅藻土的孔隙結構特征參數Tab.2 Characteristic parameters of pore structure of diatomaceous earth under different consolidation pressures

2 SEM圖像分維數的估算

2.1 數字圖像二值化

考慮已經存儲在計算機內的圖像，它是由大小為k*的像素點組成。首先把圖像進行二值化處理，使得圖像上的每一個像素點為黑或白兩種顏色。對本文jpg格式的圖片，先用最佳分割閾值的方法對其二值化，得到背景像素為土顆粒(白色)和目標像素為孔隙(黑色)的二值圖像。最佳閾值分割法源自圖像處理中的分形理論，其基本原理[18]是確定某個灰度閾值，然后判斷圖像中像素灰度值與閾值的大小關系進行圖像分割。采用該法，可以清晰界定物質邊界輪廓，將圖像轉化為黑白兩色，有利于后續(xù)分析。

在圖像進行二值化處理時，當閾值過大時，有時候會將部分土顆粒誤認為孔隙，因而須確定最佳閾值，確保大部分土顆粒的輪廓被界定出來。

本文采用如下最佳分割閾值的二值化算法：

(1)求出圖像中最小和最大的灰度值fmin和fmax，將閾值的初始值表示為：

l1=(fmin+fmax)/2

(2)由上一步得出的閾值將數字圖像分成前景和背景兩部分，求出兩部分的平均灰度值f1和f2：

(1)

(2)

式中：f(i,j)表示平面圖像對應的(i,j)點灰度值；w(i,j)表示對應的權重系數，一般取w(i,j)=1。

(3)求出新的閾值ln+1=(fA+fB)/2，若|ln+1-ln|≤0.001則循環(huán)結束，否則，轉步驟(2)重新計算。

由上述算法得到固結壓力為0， 2.0， 4.0，8.0 MPa的灰度閾值分別為96，107，101，112。

首先將掃描電鏡圖像轉換為灰度圖像，其中灰度較暗區(qū)域為孔隙，并分析灰度值，再利用上述算法得到閾值l來提取SEM圖像中微觀結構信息：

(3)

式中：l為最佳分割閾值算法得到的閾值。采用上述算法處理后的二值化圖像如圖2所示。由圖2可知，在固結壓力超過結構屈服壓力(1.3 MPa)后，硅藻顆粒受到破壞，土體骨架結構發(fā)生變化，在外加應力的作用下硅藻顆粒之間的大孔隙逐漸被壓縮，內部的孔隙趨于均勻密實。

在Matlab中編寫程序，將上述二值化圖像轉化成對應的數據文件，其行列數分別與二值圖的行列數相對應，而每一格數據取值1或0，取決于它所對應的像素點的顏色(如圖2所示)，1表示二值化圖像中對應的像素點為白色，0表示為黑色。

圖2 二值化后的SEM圖像Fig.2 Binary images of SEM

2.2 數字圖像分維數的計算

自相似性和分形是許多自然事物和現(xiàn)象的客觀特征，存在于自然界中的物理分形往往表現(xiàn)出某種尺度性和隨機性，即僅在某種特定的尺度范圍內從統(tǒng)計角度上表現(xiàn)出分形特征。如何描述研究對象的分形維數，這取決于研究對象的特征以及研究目的。針對不同的研究對象，可采用不同的描述方法計算其分形維數，如Hausdorff維數和相似維數。然而，對很多分形而言，這兩種維數都是難以計算的，在實際計算中我們一般采用盒維數D，它的一種等價計算公式為：

lgNk=C-Dlgk

(4)

式中：C為常數；k為盒子的尺寸；Nk表示覆蓋所有前景圖像需要尺寸為k的盒子數。

在計算時把數據文件劃分為若干塊，每一塊均為正方形，邊長包含像素的個數為n，把包含0的塊的個數記作Nkn。設一個像素點的尺寸為k*，則行和列都為n的塊邊長kn=nk*。以kn為邊長作為塊劃分，可以得到對應的盒子數。

據上述計算原理，在雙對數坐標平面中，用最小二乘法直線擬合數據點(lg(1/k)，lgNk，所得直線的斜率D就是該圖像的盒維數。將圖2中的圖像帶入編制的計算程序，擬合結果如圖3所示，不同固結壓力下的Nk值見表3。

圖3 最小二乘法擬合雙對數Fig.3 Linear fitting figure of bi-logarithm coordinates

固結壓力/MPa不同盒子邊長的Nk值1361225501002000355637320272990323963122213217621306697321035471303064134453216617058211054113030681614422594481792352561130306

圖4 分形盒維數與固結壓力的關系曲線Fig.4 Relationship between fractal dimension and consolidated pressure

由圖3可知，不同固結壓力下的SEM圖像在二值化后具有較好的分形特征。試驗數據點與擬合直線的相關性較好。因而可以用分形盒維數描述天然沉積硅藻土在破壞過程中的微觀結構變化特征，不同固結壓力下分形盒維數如圖4所示。

由圖4可知，分形盒維數隨固結壓力的增加呈上升的趨勢。結合圖1和圖2的SEM圖像可將分維的變化分為兩個階段，在壓力超過硅藻土的結構屈服壓力后，對應的分形維數出現(xiàn)了陡增段，這是由于硅藻顆粒本身破損造成的。此后隨著固結壓力的增加，分形盒維數也出現(xiàn)了顯著的變化，這是由存在于硅藻顆粒間的孔隙結構坍塌所致。因而，對本文討論的天然沉積硅藻土而言，分形盒維數既可以表征結構性的變化，也可以表征結構性消失后孔隙分布的變化。

3 結 語

(1)根據天然硅藻土的SEM成像特征，可知硅藻顆粒為圓形，由外側環(huán)狀物連接，土體內部有很多孔隙，具有很強的結構性。由壓汞試驗可知硅藻土中孔徑為0.1～1.0 μm的孔隙占大部分。

(2)在固結壓力達到屈服壓力后，孔徑大于1.0 μm的較大孔隙隨固結壓力的增加而逐漸減少，硅藻顆粒破壞，孔隙結構坍塌。并通過由灰度閾值算法得出的二值化圖像進一步得到驗證。

(3)在二值化圖像的基礎上，利用最小二乘法擬合得到了天然沉積硅藻土在不同固結壓力下的分形盒子數和盒維數。結果表明分形盒維數隨固結壓力增加而增加，因此分形特征可反映天然沉積硅藻土在固結壓力下的微觀結構性及變化。

[1] LEROUEIL S, VAULGHALL P R. The general and congruent effects of structure in natural soils and weak rocks[J]. Geotechnique, 1990, 40(3): 467- 488.

[2] CHANDLER R J. Clay sediments in depositional basin: the geotechnical cycle[J]. The Quarterly Journal of Engineering Geology and Hydrology, 2000, 33(1): 7- 39.

[3] COTECCHIA F,CHANDLER R J. A general framework for the mechanical behavior of clays[J]. Geotechnique, 2000, 50(4): 431- 447.

[4] 謝定義, 齊吉琳, 張振中. 考慮土結構性的本構關系[J]. 土木工程學報, 2000, 33(4): 35- 41. (XIE Dingyi, QI Jilin, ZHANG Zhenzhong. A constitutive model laws considering soil structural properties[J]. China Civil Engineering Journal, 2000, 33(4): 35- 41. (in Chinese))

[5] QI Jilin, XIE Dingyi. Soil structure and its quantitative parameter[C]∥Proceedings of the Fifteenth International Conference on Soil Mechanics and Geotechnical Engineering. Lisse: Balkema, 2001: 239- 242.

[6] 王國欣, 黃宏偉, 肖樹芳. 軟土微結構特征的試驗研究[J]. 水利學報, 2005, 36(2): 190- 196. (WANG Guoxin, HUANG Hongwei, XIAO Shufang. Experimental study on micro-structural characteristics of soft soil[J]. Journal of Hydraulic Engineering, 2005, 36(2): 190- 196. (in Chinese))

[7] 王清, 王鳳艷, 肖樹芳. 土微觀結構特征的定量研究及其在工程中的應用[J]. 成都理工大學學報, 2001, 28(2): 148- 153. (WANG Qing, WANG Fengyan, XIAO Shufang. A quantitative study of the micro-structure characteristics of soil and its application to the engineering[J]. Journal of Chengdu University of Technology, 2001, 28(2): 148- 153. (in Chinese))

[8] 洪振舜, 立石義孝, 鄧永鋒. 天然硅藻土的應力水平與孔隙空間分布的關系[J]. 巖土力學, 2004, 25(7): 1023- 1026. (HONG Zhenshun, ITATEISHI Yoshitak, DENG Yongfeng. Relationship between entrance pore distribution and stress level of natural sedimentary diatomaceous soil[J]. Rock and Soil Mechanics， 2004, 25(7): 1023- 1026. (in Chinese))

[9] HONG Zhenshun, TATEISHI Yoshitak, HAN Jie. Experimental study of Macro- and Micro-behavior of natural diatomite[J]. Journal of Geotechnical and Geoenvironmental Engineering, ASCE, 2006, 132(5): 603- 610.

[10] 陳嘉鷗, 葉斌， 郭素杰， 等. 珠江三角洲軟土SEM微結構定量研究[J]. 電子顯微學報, 2001, 20(1): 72- 75. (CHEN Jia’ou, YE Bin, GUO Sujie, et al. Quantitative analysis of SEM image of soft ground microstructure in Zhujiang delta[J]. Journal of Chinese Electron Microscopy Society, 2001, 20(1): 72- 75. (in Chinese))

[11] 雷華陽， 肖樹芳. 天津海積軟土微觀結構與工程性質初探[J]. 地質與勘察, 2002, 38(6): 81- 85. (LEI Huayang, XIAO Shufang. Primary exploration on engineering properties of marine soft soil microstructure in Tianjin region[J]. Geology and Prospecting, 2002, 38(6): 81- 85. (in Chinese))

[12] 雷華陽， 肖樹芳. 天津軟土的次固結變形特性研究[J]. 工程地質學報, 2002, 10(4): 339- 343. (LEI Huayang, XIAO Shufang. Study on secondary-consolidation deformation characteristics of soft soil in Tianjin[J]. Journal of Engineering Geology, 2002, 10(4): 339- 343. (in Chinese))

[13] 盧佩霞， 曹玲， 徐永福. 非飽和土力學的分形理論[J]. 地下空間與工程學報, 2015, 11(2): 375- 381. (LU Peixia, CAO Ling, XU Yongfu. Fractal theory of unsaturated soil mechanics[J]. Chinese Journal of Underground Space and Engineering, 2015, 11(2): 375- 381. (in Chinese))

[14] 張先偉, 王常明, 馬棟和. 軟土微觀結構表面起伏的三維可視化及分形維數的計算[J]. 應用基礎與工程科學學報, 2012, 20(1): 104- 112. (ZHANG Xianwei, WANG Changming, MA Donghe. 3D visualization and fractal dimension of soft clay’s microstructure surface undulation[J]. Journal of Basic Science and Engineering, 2012, 20(1): 104- 112. (in Chinese))

[15] 王升福, 楊平, 劉貫榮, 等. 人工凍融軟黏土微觀孔隙變化及分形特性分析[J]. 巖土工程學報, 2016, 38(7): 1254- 1261. (WANG Shengfu, YANG Ping, LIU Guanrong, et al. Micro pore change and fractal characteristics of artificial freeze thaw soft clay[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(7): 1254- 1261. (in Chinese))

[16] YOSHINAKA R. Mechanical properties for sedimentary soft rocks[M]. Tokyo： Committee of Sedimentary Soft Rocks in Geotechnical Engineering, Japanese Society of Soil Mechanics and Foundation Engineering, 1987: 69- 116. (in Japanese)

[17] TATEISHI Y. Geotechnical properties of diatom earth and stability of surface layer for the cut slope[D]. Saga: Saga University, 1997. (in Japanese)

[18] 方磊, 李春苗， 龔文忠， 等. 閾值分割法處理土微結構圖像質量評價[J]. 地下空間與工程學報, 2013, 9(增刊2): 1045- 1049. (FANG Lei, LI Chunmiao, GONG Wenzong, et al. Quality evaluation of soil micro-structure images processed by thresholding[J]. Chinese Journal of Underground Space and Engineering, 2013, 9(Suppl2): 1045- 1049. (in Chinese))

FractalcharacteristicsofnaturalsedimentarydiatomaceousearthbasedonSEMimages

WANG Yuexiang1, 2， SHI Xiusong3， JI Feng4， GU Huanda1， DING Jianwen4

(1.SchoolofCivilEngineering,SuzhouUniversityofScienceandTechnology,Suzhou215011,China; 2.NanjingTechUniversity,Nanjing211186,China; 3.DepartmentofCivilandEnvironmentalEngineering,TheHongKongPolytechnicUniversity,HongKong，China; 4.CollegeofTransportation,SoutheastUniversity,Nanjing210018,China)

SEM test, mecury injection test and high pressure triaxial test are conducted to obtain the microscopic structure and deformation of natural sedimentary diatomaceous earth. Fractal characteristics of natural sedimentary diatomaceous earth and evolutions of pore structural characteristic parameters under different consolidated pressures are analyzed by image segmentation and optimal dividing threshold algorithm. The relationship between fractal dimension and isotropic consolidation is discussed. The results show that there are many pores in the diatomaceous earth which is strongly structural. The pore diameter distribution is mainly in the range of 0.1～1.0 μm. When consolidated pressure increases to yield pressure, pores whose diameters are larger than 0.1 μm are changing obviously. The number of large pores decreases with the increase of consolidated pressure, and the diatomaceous particles are destroyed and the pore structure is collapsed. The imaging characteristic of diatomaceous earth is analyzed, and then the gray threshold under different consolidated pressures is calculated by optimum dividing threshold algorithm. Based on the corresponding gray threshold, the grey images are then transformed into the binary image. The fractal characteristics of natural sedimentary diatomaceous earth under different consolidated pressures are studied. It is proposed that the binary images have obvious fractal characteristics and the fractal dimension increases with the consolidated pressure. Fractal dimension can not only represent the variation of micro structure but also represent the variation of void distribution of natural sedimentary diatomaceous earth.

SEM images; binary; diatomaceous earth; consolidated pressure; fractal characteristics

TU411.92

A

1009-640X(2017)05-0096-07

10.16198/j.cnki.1009-640X.2017.05.014

王月香， 石修松， 吉鋒， 等. 基于SEM圖像處理的天然硅藻土分形特征分析[J]. 水利水運工程學報, 2017(5): 96-102. (WANG Yuexiang, SHI Xiusong, JI Feng, et al. Fractal characteristics of natural sedimentary diatomaceous earth based on SEM images[J]. Hydro-Science and Engineering, 2017(5): 96-102. (in Chinese))

2016-11-15

國家自然科學基金資助項目(51378327，51378118，51508369)； 江蘇省自然科學基金資助項目(BK20150289)； 江蘇省建設系統(tǒng)項目(2016ZD26)； 2017年蘇州市市級建設科研項目

王月香(1977—)， 女， 江西九江人， 講師， 博士， 主要從事巖土工程、地下工程等領域的教學與科研工作。

E-mail： wyxcjy@sina.cn