肖 佳，韓凱東，張澤的

（中南大學(xué) 土木工程學(xué)院，湖南 長(zhǎng)沙 410075）

石灰石粉資源儲(chǔ)量大、分布廣泛，作為水泥混凝土摻和料有環(huán)保、節(jié)能、利廢的重要意義［1］.水泥漿體的絮凝與漿體的流變行為密切相關(guān)［2］，顆粒間絮凝-解絮的可逆行為在短時(shí)間內(nèi)決定了漿體的觸變性，并導(dǎo)致流變參數(shù)隨時(shí)間變化［3］.學(xué)者們通常采用微觀參數(shù)表征或顯微鏡觀察水泥漿體的絮凝狀態(tài)，Rous?sel［4］通過(guò)觸變性模型擬合流變參數(shù)得到絮凝結(jié)構(gòu)因子λ和絮凝速率Athix，Ye等［3］利用絮體的平均粒徑研究了水泥基材料靜態(tài)屈服應(yīng)力與絮凝的關(guān)系.然而，這些參數(shù)都只能反映絮體的平均狀態(tài)或尺寸分布，絮體的空間分布也會(huì)影響其性能.近年來(lái)，研究者采用分形理論來(lái)描述絮體顆粒的幾何特征.在絮體生長(zhǎng)過(guò)程中，絮體的空間和尺寸分布時(shí)刻變化，簡(jiǎn)單分形維數(shù)只能描述絮體生長(zhǎng)的平均性和整體性，不能完全揭示絮體分形變化的動(dòng)力學(xué)過(guò)程［5］.研究表明，水泥漿體的分形維數(shù)與測(cè)度有關(guān)，這也暗示了絮體生長(zhǎng)的多 重分形特性［6］.Valentini等［7］通過(guò)水泥漿體的數(shù)字圖像計(jì)算其多重分形譜，定量描述了水化硅酸鈣（C?S?H）凝聚.Vahedi等［8］證明了多重分形譜可以提供絮體的內(nèi)部結(jié)構(gòu)、空間分布、尺寸分布的信息.Gao等［6］也指出多重分形在無(wú)需任何先驗(yàn)假設(shè)的條件下可以描述水泥漿體孔隙結(jié)構(gòu)的大小分布和空間排列.本文觀測(cè)了水泥-石灰石粉漿體的絮體生長(zhǎng)，定量描述了絮體生長(zhǎng)的分布特征，并研究了其絮體生長(zhǎng)特征與多重分形譜特征參數(shù)的關(guān)系.

1 試驗(yàn)

1.1 原材料及配合比

水泥（C）采用P?I 42.5拉法基瑞安基準(zhǔn)水泥，比表面積為347 m2/kg；石灰石粉（L）產(chǎn)自湖北荊門，CaCO3含量1）文中涉及的含量、水膠比等除特殊說(shuō)明外均為質(zhì)量分?jǐn)?shù)或質(zhì)量比.為99%，比表面積分別為411、608、807、1 007 m2/kg，依次記為L(zhǎng)Ⅰ、LⅡ、LⅢ、LⅣ；拌和水（W）為自來(lái)水.水泥-石灰石粉漿體的水膠比為0.4，其配合比見(jiàn)表1.

水泥-石灰石粉顆?？偙缺砻娣e（SSA）計(jì)算公式見(jiàn)式（1）［9］，計(jì)算結(jié)果見(jiàn)表1.

表1 水泥-石灰石粉漿體配合比與顆?？偙缺砻娣eTable 1 Mix proportions and total specific surface area of cement‐ground limestone pastes

式中：sC、sL分別表示水泥、石灰石粉顆粒的比表面積；φC、φL分別表示水泥、石灰石粉顆粒占固體總顆粒的體積分?jǐn)?shù).

1.2 試驗(yàn)方法

按表1配合比制備水泥-石灰石粉漿體：先將漿體慢攪60 s，靜置30 s，接著快攪120 s，再慢攪60 s，靜置90 s，然后慢攪14 min.從水泥水化的第6 min開(kāi)始，每隔2 min用針狀工具取1個(gè)針尖大小的漿體，放在載玻片上，用滴管取1滴無(wú)水乙醇滴在漿體上，再放上蓋玻片，放在顯微圖像分析儀下進(jìn)行觀察，調(diào)整光學(xué)顯微鏡的粗調(diào)和微調(diào)焦距，直到所觀察到的漿體圖像清晰，再調(diào)整載物臺(tái)的位置，找到3個(gè)合適位置的圖像留存［3，10］.每個(gè)試樣采集100像素×100像素的圖像約30幅（見(jiàn)圖1（a））；對(duì)圖像進(jìn)行二值化處理，并用黑色表示絮體區(qū)域，接著進(jìn)行邊界濾波和面積濾波，去除面積過(guò)小且不影響測(cè)試結(jié)果的黑色部分（見(jiàn)圖1（b））；最后選取合適的絮體圖像進(jìn)行分割處理（見(jiàn)圖1（c））.

圖1 水泥漿體的絮體圖像Fig.1 Floc image of cement paste

2 結(jié)果與討論

2.1 水泥-石灰石粉漿體的絮體生長(zhǎng)特征

表2為水泥-石灰石粉漿體絮體面積占比（ξ）與絮體面積平均占比表征絮體生長(zhǎng)的平均尺寸分布狀態(tài).由表2可見(jiàn)，隨著水化的進(jìn)行，漿體的絮體面積占比逐漸增大，絮體的尺寸增大.

表2 水泥-石灰石粉漿體的絮體面積占比與絮體面積平均占比Table 2ξandξˉof cement?ground limestone pastes

為探究石灰石粉對(duì)水泥漿體絮體生長(zhǎng)空間分布的影響，將絮體圖像用無(wú)窮個(gè)大小相等的小方格覆蓋（即小方格的邊長(zhǎng)要無(wú)限?。?，計(jì)算每1個(gè)小方格中絮體所占像素.若每1個(gè)小方格中絮體的像素值均接近于平均值（標(biāo)準(zhǔn)差越小），絮體分布越均勻，反之，絮體分布越分散.絮體分布標(biāo)準(zhǔn)差σ計(jì)算如式（2）所示.

式中：p為絮體圖像像素，取p=100；n為小方格個(gè)數(shù)，取n=256；ξi為第i個(gè)小方格絮體面積占比.

表3為水泥-石灰石粉漿體絮體分布標(biāo)準(zhǔn)差（）與絮體分布平均標(biāo)準(zhǔn)差（），表征絮體生長(zhǎng)的平均空間分布狀態(tài).由表3可見(jiàn)，隨著水化的進(jìn)行，水泥-石灰石粉漿體絮體分布標(biāo)準(zhǔn)差逐漸減小，表明其空間分布越來(lái)越均勻.

表3 水泥-石灰石粉漿體絮體分布標(biāo)準(zhǔn)差與平均標(biāo)準(zhǔn)差Table 3 and of cement?ground limestone pastes

表3 水泥-石灰石粉漿體絮體分布標(biāo)準(zhǔn)差與平均標(biāo)準(zhǔn)差Table 3 and of cement?ground limestone pastes

Sample C0 CLⅠ20 CLⅠ30 CLⅡ20 CLⅡ30 CLⅢ20 CLⅢ30 CLⅣ20 CLⅣ30 σ 6 min 24.66 18.84 25.37 31.79 34.95 25.37 13.78 17.13 15.25 8 min 15.46 14.13 23.23 15.55 16.30 22.21 20.19 19.43 13.52 10 min 13.76 14.42 16.68 20.57 22.31 15.13 12.64 11.98 11.98 12 min 13.54 12.37 15.12 15.12 21.24 17.06 11.15 12.60 9.45 14 min 13.80 14.07 13.65 14.71 14.49 13.08 13.89 13.37 9.43 16 min 12.17 13.44 14.35 13.17 14.26 12.09 9.29 13.36 9.54 18 min 11.74 10.18 11.07 12.90 16.62 11.41 12.49 11.41 11.04 20 min 10.21 9.21 11.34 11.75 13.54 11.04 7.78 7.54 7.27 σˉ 14.42 14.33 16.35 16.94 19.21 15.92 12.65 13.35 10.93

由表1~3可見(jiàn)：隨著顆?？偙缺砻娣e的增加，絮體面積平均占比先減小后增加，分布平均標(biāo)準(zhǔn)差先增加后減小.由此可見(jiàn)，當(dāng)顆?？偙缺砻娣eSSA增加時(shí)，絮體生長(zhǎng)尺寸減小，分布變分散，而隨著SSA的繼續(xù)增加，絮體將有更好的尺寸和空間分布.

2.2 水泥-石灰石粉漿體絮體生長(zhǎng)的多重分形特征

多重分形包括2種關(guān)系［11］：（1）各盒子（用i編號(hào)）內(nèi)物理量的概率Pi（ε）與測(cè)度ε之間為冪函數(shù)關(guān)系Pi（ε）~εα，α為奇異性指數(shù)；（2）盒子數(shù)N（ε）與測(cè)度ε之間為冪函數(shù)關(guān)系N（ε）~ε-f（α），f（α）為多重分形譜函數(shù).假設(shè)絮體圖像邊長(zhǎng)為1，將圖像用尺寸為ε×ε（ε=1/2、1/4、1/8、1/16、1/32、1/64）的網(wǎng)格覆蓋.q為權(quán)重因子，表示各概率子集在配分函數(shù)中的貢獻(xiàn)，理論上，q（-∞

若對(duì)于不同的q值，lnIq（ε）與lnε都有良好的線性相關(guān)性，說(shuō)明絮體具有標(biāo)度不變性，在此基礎(chǔ)上，若τ（q）?q曲線呈上凸形，說(shuō)明絮體具有多重分形特征.水泥漿體的多重分形特征見(jiàn)圖2.由圖2可見(jiàn)：lnIq（ε）與lnε呈線性關(guān)系，尤其對(duì)于階次較高的部分（q>0），表現(xiàn)出非常嚴(yán)格的線性特征，結(jié)合τ（q）?q曲線呈上凸形，證明了水泥漿體絮體生長(zhǎng)具備多重分形特征.

圖2 水泥漿體的多重分形特征Fig.2 Multifractal characteristics of cement paste

多重分形理論包括3個(gè)重要的組成部分：廣義分形維數(shù)Dq、奇異性指數(shù)α以及多重分形譜函數(shù)f（α），計(jì)算公式見(jiàn)式（6）~（8）.

水泥-石灰石粉漿體絮體生長(zhǎng)的多重分形特征見(jiàn)圖3.由圖3可見(jiàn)：水泥-石灰石粉漿體絮體生長(zhǎng)的廣義分形維數(shù)Dq為反“S”遞減函數(shù)曲線，Dq隨著q增大而減小，曲線在q=0處出現(xiàn)向另一方向彎折的趨勢(shì)；當(dāng)q<0時(shí)，Dq隨顆?？偙缺砻娣e的增加先增后減，而在q>0時(shí)會(huì)出現(xiàn)相反的結(jié)果；絮體生長(zhǎng)的多重分形譜均為左鉤狀，且隨著顆?？偙缺砻娣e的增加，多重分形譜先向右發(fā)展，然后向左發(fā)展.

圖3 水泥-石灰石粉體漿體絮體生長(zhǎng)的多重分形特征Fig.3 Multifractal characteristics of flocs growth in cement?ground limestone pastes

多重分形譜中：譜寬Δα可描述絮體的局部特征，代表著絮體的空間分布情況，Δα越大，絮體分布越不均勻；左譜寬ΔαL反映了絮體中大空隙的數(shù)目；右譜寬ΔαR反映了絮體中小空隙的數(shù)目；譜差Δf反映多重分形譜的不對(duì)稱性，代表著絮體的尺寸分布情況，Δf<0時(shí)，多重分形譜曲線呈現(xiàn)左鉤狀，反之多重分形譜曲線呈現(xiàn)右鉤狀；右譜值f（αmax）反映小概率對(duì)象即大粒徑絮體的數(shù)量（αmax為奇異性指數(shù)的最大值）；左譜值f（αmin）反映大概率對(duì)象即小粒徑絮體的數(shù)量（αmin為奇異性指數(shù)的最小值）.將Δα、ΔαL、ΔαR稱為絮體的空間分布參數(shù)，Δf、f（αmax）、f（αmin）稱為絮體的尺寸分布參數(shù).最小分形維數(shù)Dmin與使f（α）達(dá)到最大的奇異性指數(shù)α0一起反映了最或然子集的性質(zhì)，即反映了絮體的整體幾何特征，稱為絮體的整體狀態(tài)參數(shù).

圖4為水泥-石灰石粉顆粒總比表面積對(duì)絮體生長(zhǎng)的空間分布、尺寸分布和整體狀態(tài)參數(shù)的影響.由圖4可見(jiàn)：當(dāng)顆?？偙缺砻娣eSSA增加時(shí)，Δα、ΔαL、ΔαR增加，表明絮體間的小空隙與大空隙數(shù)目增加，絮體分布變得分散，但隨著SSA的繼續(xù)增加，絮體生長(zhǎng)的空間分布出現(xiàn)相反的結(jié)果；f（αmin）先減小后增大，即小粒徑的絮體數(shù)目先減小后增大，f（αmax）大致保持不變，即大粒徑的絮體數(shù)目基本穩(wěn)定，Δf先增大后減小，即絮體的尺寸先減小后增大；Dmin先減小后增大，α0先增大后減小，延緩絮體生長(zhǎng)，而隨著SSA的繼續(xù)增加將促進(jìn)絮體生長(zhǎng).

2.3 水泥-石灰石粉漿體的絮體生長(zhǎng)特征的多重分形表征

結(jié)合表2、3與圖4可知，絮體面積平均占比與f（αmin）正相關(guān)，與Δf負(fù)相關(guān)，即絮體生長(zhǎng)的尺寸平均分布狀態(tài)可以用f（αmin）與-Δf來(lái)描述；絮體分布平均標(biāo)準(zhǔn)差與Δα、ΔαL、ΔαR正相關(guān)，即絮體生長(zhǎng)的平均空間分布狀態(tài)可以用-Δα、-ΔαL、-ΔαR來(lái)表征.由于絮體面積平均占比和絮體分布平均標(biāo)準(zhǔn)差可以描述絮體生長(zhǎng)特征，稱之為絮體生長(zhǎng)特征參數(shù)；絮體生長(zhǎng)的空間分布參數(shù)、尺寸分布參數(shù)和整體狀態(tài)參數(shù)可以表征絮體生長(zhǎng)的多重分形特征稱之為多重分形譜參數(shù).將絮體生長(zhǎng)特征參數(shù)及多重分形譜參數(shù)進(jìn)行歸一化處理（處理后的數(shù)據(jù)以下標(biāo)n表示），其計(jì)算公式為：

圖4 水泥-石灰石粉顆粒總比表面積對(duì)絮體生長(zhǎng)的空間分布、尺寸分布和整體狀態(tài)參數(shù)的影響Fig.4 Effect of total specific surface area of cement?ground limestone on spatial distribution，size distribution and whole state parameters

式中：xi、xin為第i個(gè)歸一化前、后的數(shù)據(jù)；xmin、xmax為原始的最小值和最大值.

表4為歸一化的絮體生長(zhǎng)特征參數(shù)及多重分形譜參數(shù)，采用BP神經(jīng)網(wǎng)絡(luò)對(duì)表4數(shù)據(jù)進(jìn)行訓(xùn)練，以絮體生長(zhǎng)特征參數(shù)為輸入，多重分形譜參數(shù)為輸出，建立水泥-石灰石粉漿體絮體生長(zhǎng)特征參數(shù)及多重分形譜特征參數(shù)的關(guān)系.通過(guò)計(jì)算歸一化后的絮體生長(zhǎng)特征參數(shù)及多重分形譜特征參數(shù)的標(biāo)準(zhǔn)差σ′，對(duì)其進(jìn)行關(guān)聯(lián)度分析，得到絮體生長(zhǎng)尺寸分布與f（αmin）n、-Δfn的關(guān)系見(jiàn)圖5，絮體生長(zhǎng)空間分布-與-ΔαLn、-ΔαRn、-Δαn的關(guān)系見(jiàn)圖6.由圖5可見(jiàn)，-Δfn的標(biāo)準(zhǔn)差小于f（αmin）n，Δfn可以更好地描述絮體生長(zhǎng)的尺寸分布，且與絮體生長(zhǎng)的尺寸分布負(fù)相關(guān).由圖6可見(jiàn)，-Δαn的標(biāo)準(zhǔn)差最小，-ΔαLn居中，-ΔαRn最大，Δαn可以更好地描述絮體生長(zhǎng)的空間分布，且與絮體生長(zhǎng)空間分布負(fù)相關(guān).

圖5 絮體生長(zhǎng)尺寸分布與f（αmin）n、-Δfn的關(guān)系Fig.5 Relationship between size distribution of floc growth and f（αmin）n and（-Δfn）

圖6 絮體生長(zhǎng)空間分布-與-ΔαLn、-ΔαRn、-Δαn的關(guān)系Fig.6 Relationship between spatial distribution of floc growth- and-ΔαLn、-ΔαRn、-Δαn

表4 歸一化后的絮體生長(zhǎng)特征參數(shù)及多重分形譜參數(shù)Table 4 Characteristic parameters of floc growth and multifractal spectrun after normalization

綜上，多重分形理論能有效、準(zhǔn)確地表征水泥-石灰石粉漿體絮體生長(zhǎng)特征，清晰地分析絮體的分布情況，Δα、ΔαL、ΔαR可以用來(lái)表征絮體生長(zhǎng)的空間分布，且譜寬Δα的準(zhǔn)確性更高，f（αmin）與Δf可用于描述絮體生長(zhǎng)的尺寸分布，且譜差Δf的準(zhǔn)確性更高.

3 結(jié)論

（1）當(dāng)顆?？偙缺砻娣e增加時(shí)，水泥-石灰石粉漿體的絮體面積平均占比先減小后增加，絮體分布平均標(biāo)準(zhǔn)差先增加后減小.

（2）當(dāng)顆?？偙缺砻娣e（SSA）增加時(shí)，絮體間的小空隙與大空隙數(shù)目均增加，絮體分布變得分散，小粒徑的絮體數(shù)目減小，大粒徑的絮體數(shù)目基本穩(wěn)定，絮體的尺寸減小，整體上延緩了絮體生長(zhǎng)，而隨著SSA的繼續(xù)增加，絮體生長(zhǎng)將呈現(xiàn)相反的結(jié)果.

（3）水泥-石灰石粉漿體的絮體生長(zhǎng)具備多重分形特征.多重分形譜的譜寬、左譜寬和右譜寬可以用來(lái)表征絮體生長(zhǎng)的空間分布，且譜寬更為準(zhǔn)確，與絮體生長(zhǎng)的空間分布負(fù)相關(guān).譜差與左譜值可用于描述絮體生長(zhǎng)的尺寸分布，且譜差更為準(zhǔn)確，與絮體生長(zhǎng)的尺寸分布負(fù)相關(guān).