樊文虎 楊平 王升福

摘 要：為探究軟土融沉變形的細(xì)觀機(jī)理，以上海地區(qū)粉質(zhì)黏土為研究對(duì)象，通過(guò)開展凍融（freeze-thaw，F(xiàn)-T）試驗(yàn)和X射線計(jì)算機(jī)斷層掃描成像技術(shù)（X-ray Computed Tomography，X-CT）試驗(yàn)，并結(jié)合圖像處理技術(shù)和分形理論，分析有無(wú)補(bǔ)水條件飽和粉質(zhì)黏土凍融前后細(xì)觀結(jié)構(gòu)的變化規(guī)律。試驗(yàn)結(jié)果表明，不補(bǔ)水條件下發(fā)生了凍融頸縮現(xiàn)象，而補(bǔ)水條件下未發(fā)生；凍融后CT灰度強(qiáng)度平均值的變化與土體含水率、孔隙率和干密度的變化之間均呈良好的線性關(guān)系；基于提出的土體CT圖像細(xì)觀孔隙識(shí)別方法，發(fā)現(xiàn)凍融顯著改變了土體的橫截面細(xì)觀孔隙率和平均孔徑，補(bǔ)水條件最大變化發(fā)生在未凍區(qū)，而不補(bǔ)水條件則發(fā)生在鄰近最終凍結(jié)鋒面處；分形維數(shù)與橫截面孔隙率和平均孔徑均呈良好的線性關(guān)系。此研究表明，有無(wú)補(bǔ)水條件下凍融后土體細(xì)觀結(jié)構(gòu)的變化具有明顯差異性，補(bǔ)水條件下凍融更顯著改變土體細(xì)觀孔隙結(jié)構(gòu)。

關(guān)鍵詞：凍融試驗(yàn)；X-CT試驗(yàn)；細(xì)觀結(jié)構(gòu)；灰度強(qiáng)度；CT圖像處理；細(xì)觀孔隙

中圖分類號(hào)：TU445 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1006-8023（2023）03-0182-09

Abstract：In order to investigate the mesoscopic mechanism of thaw settlement for soft soil， freeze-thaw （F-T） test and X-ray computed tomography （X-CT）， combined with image processing technology and fractal theory， were adopted to study the changes of meso-structure for silty clay in Shanghai before and after F-T with or without water supply during freezing. Results showed that， the freeze-necking phenomenon was observed in soil specimen without water supply， but it did not happen in soil specimen with water supply. There was a good linear relationship between the change of the average CT gray intensity after F-T and the changes in soil moisture content， porosity， and dry density. Based on the proposed soil mesopore recognition method， it was found that freeze-thaw significantly changed the transverse-sectional porosity and average mesopore size， and the biggest change happened in unfrozen area with water supply， but happened near the final freezing front without water supply. Fractal dimensions had better linear relation with transverse-sectional porosity and average mesopore size， respectively. The study showed that the change in soil meso-structure with water supply during freezing differed significantly from that without water supply during F-T， and F-T altered soil mesopore structure obviously under water-supply condition.

Keywords：Freeze-thaw test; X-ray computed tomography; meso-structure; gray intensity; CT image processing; mesopore

基金項(xiàng)目：國(guó)家自然科學(xué)基金項(xiàng)目（52178337；52108323）；江蘇省自然科學(xué)基金項(xiàng)目（BK20210006）；金陵科技學(xué)院高層次人才科研啟動(dòng)基金（jit-b-202125）

第一作者簡(jiǎn)介：樊文虎，博士，講師。研究方向?yàn)榄h(huán)境巖土與地下工程。E-mail： fwh60131@jit.edu.cn

*通信作者：楊平，碩士，教授。研究方向?yàn)榄h(huán)境巖土與地下工程。E-mail： yangping@njfu.edu.cn

0 引言

人工凍結(jié)法已廣泛應(yīng)用于軟土地區(qū)地鐵隧道的修建，確保盾構(gòu)的安全始發(fā)、接收以及聯(lián)絡(luò)通道的安全貫通[1-2]，也被用于解決復(fù)雜地下工程問(wèn)題[3-4]，但工后地層的融沉變形預(yù)測(cè)和控制是亟須解決的關(guān)鍵問(wèn)題。室內(nèi)試驗(yàn)通常采用封閉凍結(jié)（即不補(bǔ)水條件）和開放凍結(jié)（即補(bǔ)水條件），分別模擬凍結(jié)工程中凍結(jié)管間土體的凍結(jié)和凍結(jié)壁邊緣土體的凍結(jié)，獲得土體的融沉系數(shù)來(lái)估算土層的融沉變形[5-6]。凍融后土體細(xì)微觀結(jié)構(gòu)的變化是人工凍土宏觀融沉變形的根本原因，也是凍土工程研究的熱點(diǎn)問(wèn)題，其中凍融對(duì)土體細(xì)觀結(jié)構(gòu)的影響研究對(duì)揭示其融沉變形的機(jī)理起重要作用。研究有無(wú)補(bǔ)水條件下凍融前后土體細(xì)觀結(jié)構(gòu)的差異，可為實(shí)際工程中不同凍結(jié)區(qū)域的土層融沉變形預(yù)測(cè)和控制提供重要理論依據(jù)。

X射線計(jì)算機(jī)斷層掃描成像技術(shù)（X-ray Computed Tomography， X-CT），因其具有無(wú)損觀測(cè)物體內(nèi)部結(jié)構(gòu)的優(yōu)點(diǎn)，除了用于常溫土和凍土的研究，近年來(lái)越來(lái)越多的學(xué)者運(yùn)用其研究?jī)鋈趯?duì)土體滲透特性以及體積變化的影響。Santa等[7]發(fā)現(xiàn)粉質(zhì)黏土經(jīng)凍融循環(huán)后，其孔徑有顯著提高，進(jìn)而增加了孔隙連通性并改變了孔隙方向。Xu等[8]建立了凍融循環(huán)影響下含鹽黃土孔隙與滲透系數(shù)之間的關(guān)系。Nishimura等[9]采用CT試驗(yàn)揭示了凍融循環(huán)對(duì)黏土體積變化影響的機(jī)理。Wang等[10-11]發(fā)現(xiàn)原狀非飽和軟黏土凍融后（不補(bǔ)水條件下凍結(jié)）試樣發(fā)生了凍融頸縮現(xiàn)象，并定量分析了凍融后體積和CT灰度強(qiáng)度變化。有關(guān)凍融對(duì)土體內(nèi)部孔隙及裂隙的影響研究不多。蔡正銀等[12]研究了濕干及濕干凍融耦合循環(huán)作用下膨脹土內(nèi)部裂隙的演化規(guī)律。Fan等[13]基于CT圖像定量分析了凍融前后橫截面孔隙率、孔徑、裂隙率和縱截面裂隙方向的頻數(shù)分布。

分形理論不僅能深化對(duì)土體孔隙結(jié)構(gòu)的認(rèn)識(shí)，還可為解決孔隙結(jié)構(gòu)中的復(fù)雜問(wèn)題提供新的工具。針對(duì)凍融對(duì)土體宏觀特性的影響，眾多學(xué)者多采用微觀試驗(yàn)（如掃描電子顯微鏡（Scanning Electron Microscope， SEM）和壓汞（Mercury Intrusion Porosimetry， MIP）試驗(yàn)），并結(jié)合分形理論，對(duì)土體孔隙、裂隙和顆粒等形態(tài)特征方面開展了大量研究。Zhou等[14]基于MIP試驗(yàn)結(jié)果，采用熱力學(xué)模型計(jì)算得到分形維數(shù)，分析凍融后土體微觀孔隙結(jié)構(gòu)的變化。王升福等[15]發(fā)現(xiàn)軟黏土凍融及壓縮前后孔隙體積和表面積均存在分形特性。張英等[16]基于SEM和MIP試驗(yàn)定量分析了凍融循環(huán)對(duì)土體單軸抗壓強(qiáng)度影響的微觀機(jī)制。張澤等[17]開展了凍融循環(huán)對(duì)黃土微觀孔隙特征的研究。此外，Lu等[18]通過(guò)對(duì)不同凍融循環(huán)次數(shù)下的土體表面進(jìn)行拍照，發(fā)現(xiàn)土體的表面裂隙率與分形維數(shù)呈對(duì)數(shù)關(guān)系。然而有關(guān)結(jié)合分形理論與X-CT試驗(yàn)結(jié)果研究?jī)鋈趯?duì)土體細(xì)觀結(jié)構(gòu)的影響卻少見。

綜上所述，基于X-CT試驗(yàn)研究土體凍融前后細(xì)觀結(jié)構(gòu)的變化不夠深入，尤其是有無(wú)補(bǔ)水條件的差異性缺乏系統(tǒng)研究。因此，本研究通過(guò)凍融試驗(yàn)和X-CT試驗(yàn)，結(jié)合圖像處理技術(shù)和分形理論，研究有無(wú)補(bǔ)水條件的土體凍融前后細(xì)觀結(jié)構(gòu)變化，有助于進(jìn)一步揭示人工凍土的融沉變形機(jī)理。

1 試驗(yàn)土樣及試驗(yàn)方案

1.1 土樣及試樣制備

試驗(yàn)土樣選用上海地區(qū)粉質(zhì)黏土，其顆粒級(jí)配曲線如圖1所示。通過(guò)制備重塑土試樣開展相關(guān)試驗(yàn)，試樣尺寸為直徑79.8 mm、高度100 mm，具體制樣方法如下：1）按照干密度1.4 g/cm3分層擊實(shí)成直徑79.8 mm、高度130mm的圓柱土樣；2）將土樣放置在真空飽和缸內(nèi)充分飽和；3）在土樣頂部施加200 kPa的豎向壓力進(jìn)行排水固結(jié)，待試樣高度穩(wěn)定（1 h內(nèi)高度變化量≤0.01 mm），將其切削成高度100 mm的土樣，即為重塑土試樣，其基本物理特性指標(biāo)見表1。

1.2 凍融（freeze-thaw，F(xiàn)-T）試驗(yàn)

為模擬土體自下而上單向凍結(jié)，冷端設(shè)置在試樣底部，暖端設(shè)置在試樣頂部，該凍融試驗(yàn)設(shè)備主要包括溫控環(huán)境箱、溫度和位移傳感器、試樣筒、銅制頂板和底板、補(bǔ)水系統(tǒng)等，如圖2所示。本研究采用一種簡(jiǎn)易圓形盤管式補(bǔ)水裝置，與銅制頂板中蓄水槽連接，可實(shí)現(xiàn)試樣頂端補(bǔ)水。溫控環(huán)境箱的控溫范圍為-30～50 ℃，且控溫精度達(dá)0.2 ℃；5個(gè)T型溫度傳感器沿試樣高度布置，監(jiān)測(cè)土體內(nèi)部的溫度變化，銅制頂板和底板中各布置1個(gè)T型溫度傳感器，監(jiān)測(cè)頂板和底板的溫度變化，測(cè)量精度為0.1 ℃；每個(gè)試樣采用2個(gè)位移傳感器監(jiān)測(cè)土樣的豎向位移，測(cè)量精度為0.002 mm。

凍融試驗(yàn)的主要步驟如下：1）恒溫階段，將重塑土試樣放入試樣筒中，組裝相關(guān)部件，并置于環(huán)境箱內(nèi)，環(huán)境箱、頂板和底板溫度均設(shè)為1 ℃，恒溫6 h；2）凍結(jié)階段，調(diào)節(jié)底板溫度至-5 ℃，保持環(huán)境箱和頂板溫度不變，按試驗(yàn)方案要求確定凍結(jié)過(guò)程中是否進(jìn)行補(bǔ)水，凍結(jié)72 h；3）融化階段，關(guān)閉補(bǔ)水系統(tǒng)以及冷浴，調(diào)節(jié)環(huán)境箱溫度至20 ℃，使試樣開始融沉，直至每2 h內(nèi)變形小于0.05 mm，即融沉穩(wěn)定。

1.3 X-CT掃描試驗(yàn)

X-CT掃描設(shè)備主要包括X射線發(fā)射源、平板探測(cè)器和旋轉(zhuǎn)臺(tái)上的樣品（本研究中為裝入試樣筒中的土樣）。在試驗(yàn)過(guò)程中，X射線發(fā)射源會(huì)射出錐形X射線束穿透樣品，由于樣品內(nèi)部物質(zhì)組成和密度的差異，造成X射線發(fā)生不同程度的衰減，平板探測(cè)器可接收樣品不同位置處的射線信號(hào)，并在平面投影上呈現(xiàn)不同的灰度。與此同時(shí)，樣品會(huì)隨旋轉(zhuǎn)臺(tái)以一定的速率旋轉(zhuǎn)一周，進(jìn)而平板探測(cè)器可獲得不同旋轉(zhuǎn)角下樣品的投影，經(jīng)過(guò)復(fù)雜的數(shù)學(xué)運(yùn)算[19]，可重建樣品的三維CT灰度圖像。掃描土樣橫、縱斷面分別產(chǎn)生1 024張圖像，每個(gè)橫截面的尺寸為1 024 px×1 024 px，間隔為0.11 mm，細(xì)節(jié)分辨率可達(dá)0.1 mm。為防止掃描過(guò)程中土樣發(fā)生擾動(dòng)，且避免移除試樣筒對(duì)土樣結(jié)構(gòu)造成的損傷，掃描過(guò)程中試樣筒不移除。

1.4 試驗(yàn)流程

為研究有無(wú)補(bǔ)水條件土體凍融前后細(xì)觀結(jié)構(gòu)的差異，2個(gè)重塑土試樣凍融前分別進(jìn)行X射線掃描，隨后進(jìn)行凍融試驗(yàn)（其中一個(gè)試樣在補(bǔ)水條件下凍結(jié)，另一個(gè)試樣在不補(bǔ)水條件下凍結(jié)），凍融后再分別進(jìn)行X射線掃描。為使土樣凍融前后相同位置處的CT數(shù)據(jù)具有可比性，在旋轉(zhuǎn)臺(tái)與試樣筒之間做標(biāo)記，確保同一試樣凍融前后2次掃描中位于旋轉(zhuǎn)臺(tái)上的相同位置。凍融后將土樣自上而下平均劃分為5層（A、B、C、D、E），對(duì)每層土取樣進(jìn)行含水率試驗(yàn)，并在土層A、C和E中取樣進(jìn)行密度試驗(yàn)。

2 CT圖像數(shù)據(jù)處理

通過(guò)VG Studio Max軟件可觀測(cè)土樣的三維CT圖像，獲取土樣上任意一點(diǎn)的CT灰度強(qiáng)度（Gray intensity，GI）和截取任意切面。為消除溫度器和試樣筒側(cè)壁對(duì)土樣結(jié)構(gòu)的影響，使凍融前后數(shù)據(jù)具有類比性，采用圓柱切面在土樣凍融前后相同位置裁剪圓柱形土樣，其橫截面的尺寸為460 px×460 px，高度與試樣高度一致。另外，由于X射線的散射影響，試樣的縱截面兩端均存在厚度5 mm的端部偽影，會(huì)影響數(shù)據(jù)處理的結(jié)果，因此將端部偽影的范圍截除。沿試樣高度選取18個(gè)典型的橫斷面，提取每個(gè)橫截面的平均灰度強(qiáng)度（GIA），對(duì)比分析相同土樣凍融前后的變化，但無(wú)法表征凍融后土體細(xì)觀孔隙結(jié)構(gòu)的變化，需要對(duì)CT圖像進(jìn)一步處理分析。

圖3為識(shí)別土體橫截面細(xì)觀孔隙的圖像處理過(guò)程。首先采用Brun等[20]提出的算法盡可能消除環(huán)狀偽影，然后采用中值濾波降低圖像噪聲，經(jīng)過(guò)自適應(yīng)直方圖均衡化算法提高圖片中局部的對(duì)比度，使細(xì)觀孔隙顯現(xiàn)更清晰。處理后的圖像可以發(fā)現(xiàn)，試樣邊緣位置偏亮而靠近中間位置偏暗，若采用全局閾值分割的方法會(huì)使細(xì)觀孔隙提取的結(jié)果不準(zhǔn)確，因此采用一種局部自適應(yīng)閾值法將圖像二值化，并基于人眼對(duì)比將細(xì)觀孔隙識(shí)別出來(lái)。此外，單個(gè)像素往往是圖片中的噪聲，需要進(jìn)一步去除，因此細(xì)觀孔隙至少由2個(gè)相連像素組成，具體的處理細(xì)節(jié)詳見文獻(xiàn)[13]?；诒狙芯坎捎肅T設(shè)備的圖像精度以及上述處理方法，識(shí)別的細(xì)觀孔隙孔徑均大于0.1 mm，其相關(guān)的特征信息（如孔隙面積、平均孔徑和分形維數(shù)等）可通過(guò)ImageJ軟件（1.53t版本，美國(guó)）獲取。

3 凍融前后CT灰度圖像分析

3.1 土樣CT灰度圖像特征

CT灰度圖像中，土樣的灰度強(qiáng)度反映土體的密度，越亮的區(qū)域?qū)?yīng)的GI越大，反映該區(qū)域土樣的密度越大，反之，越暗的區(qū)域?qū)?yīng)的GI越小，密度越小。圖4為有無(wú)補(bǔ)水條件下試樣凍融前后典型的縱截面和橫截面，圖中偏亮的區(qū)域?yàn)樵嚇樱渲車h(huán)狀偏暗的區(qū)域?yàn)樵嚇油?，試樣?nèi)部偏暗的零星部分即為細(xì)觀孔隙。還可發(fā)現(xiàn)，縱截面土樣頂部和底部均存在端部偽影，橫截面上有環(huán)狀偽影，這些偽影均會(huì)影響識(shí)別土中細(xì)觀孔隙的準(zhǔn)確性，需盡可能消除。

由圖4（a）可見，補(bǔ)水條件下凍融后試樣高度增加了5 mm，不補(bǔ)水條件下僅增加了2.3 mm，相較于凍融前均發(fā)生了膨脹，說(shuō)明凍融改變了土體內(nèi)部的骨架，使土體難以恢復(fù)至凍融前的狀態(tài)。相較于不補(bǔ)水條件，補(bǔ)水條件下凍融后土樣的豎向變形更大，這是由于補(bǔ)水條件下凍結(jié)過(guò)程中吸入了較多的水分，并凍結(jié)形成更多的冰晶和冰透鏡體，但在自重應(yīng)力作用下融沉，土體內(nèi)部孔隙水難以排出。由圖4（b）可知，不補(bǔ)水條件下凍融后土樣頂部（鄰近暖端）發(fā)生凍融頸縮現(xiàn)象，頸縮長(zhǎng)度為30.71 mm，半徑減小了0.85 mm，但在補(bǔ)水條件下并未發(fā)生。Wang等[10]研究發(fā)現(xiàn)在不補(bǔ)水條件下非飽和黏土在暖端也發(fā)生凍融頸縮現(xiàn)象，而本研究采用的飽和粉質(zhì)黏土，結(jié)果充分說(shuō)明，不補(bǔ)水條件下鄰近暖端土體中的水分受凍結(jié)吸力的作用，不斷向凍結(jié)鋒面處遷移，致使其不斷失水，進(jìn)而發(fā)生凍融頸縮現(xiàn)象，但補(bǔ)水條件下，鄰近暖端土體凍結(jié)過(guò)程中失水的同時(shí)，外界水會(huì)補(bǔ)給，其綜合效應(yīng)不足以發(fā)生凍融頸縮現(xiàn)象。

3.2 土樣灰度強(qiáng)度分布

圖5為凍融前后每個(gè)橫截面的GIA（公式中用GIA表示）和灰度強(qiáng)度變化量（公式中用ΔGIA表示）沿試樣高度的分布，計(jì)算公式見式（1），其中最終凍結(jié)鋒面的高度根據(jù)試樣凍結(jié)完成時(shí)的溫度分布確定，約為76.4 mm。

式中，GIA，b和GIA，a分別為同一試樣同一高度橫截面凍融前后的平均灰度強(qiáng)度。

由圖5（a）和圖5（b）可見，有無(wú)補(bǔ)水條件凍融前重塑土樣沿試樣高度的GIA差值均在20以內(nèi)，且不同試樣GIA差值在10左右（補(bǔ)水條件下約為695，不補(bǔ)水條件下約為704），說(shuō)明重塑土試樣相對(duì)較均勻。經(jīng)凍融后，GIA值沿試樣高度發(fā)生不同程度的變化，補(bǔ)水條件下土層A中鄰近暖端截面的GIA值增大，不補(bǔ)水條件下土層A和B中多數(shù)截面的GIA值增大，而其余土層的GIA值均減小，表明土體凍結(jié)過(guò)程中孔隙水從暖端向冷端遷移，融沉后未凍區(qū)更密實(shí)，而凍結(jié)區(qū)更疏松。從圖5（c）可知，補(bǔ)水條件相較于不補(bǔ)水條件，沿試樣高度灰度強(qiáng)度變化量總體偏小，表明補(bǔ)水條件較不補(bǔ)水條件，凍融后土體更加疏松，尤其在鄰近最終凍結(jié)鋒面處更顯著。

3.3 灰度強(qiáng)度與物理指標(biāo)變化量之間的關(guān)系

類似ΔGIA的定義，圖6（a）—圖6（c）分別給出了有無(wú)補(bǔ)水條件下凍融后含水率、孔隙率和干密度的變化量（Δω、Δn、Δρd）沿試樣高度的分布，其凍融前后具體數(shù)值詳見文獻(xiàn)[21]。圖6（d）為凍融后每層土的平均灰度強(qiáng)度變化量（ΔGL，IA）沿試樣高度的分布。不難發(fā)現(xiàn)，ΔGL，IA與Δω、Δn呈負(fù)相關(guān)（即ΔGL，IA越大，對(duì)應(yīng)的Δω和Δn越?。?，而ΔGL，IA與Δρd呈正相關(guān)（即ΔGL，IA越大，對(duì)應(yīng)的Δρd越大）。圖7建立了凍融前后Δω、Δn、Δρd和ΔGL，IA之間的定量關(guān)系，由圖7可以發(fā)現(xiàn)，Δω、Δn、Δρd和ΔGL，IA之間均存在明顯的線性關(guān)系，決定系數(shù)R2均高于0.8，說(shuō)明線性相關(guān)性較好，表明土體凍融后灰度強(qiáng)度的變化可以反映土體凍融后含水率、孔隙率和干密度等物理指標(biāo)的變化。

4 凍融前后細(xì)觀孔隙結(jié)構(gòu)分析

4.1 凍融前后橫截面細(xì)觀孔隙率分布

為定量分析凍融后土體細(xì)觀孔隙率的變化，定義橫截面細(xì)觀孔隙率（Sp）和其凍融后變化量（ΔSp）進(jìn)行比較，可由下式計(jì)算

式中：Ap為孔隙面積；A為橫截面的面積；Sp，b和Sp，a分別為同一試樣同一高度凍融前后的橫截面細(xì)觀孔隙率。

圖8為有無(wú)補(bǔ)水條件下土體凍融前后典型橫截面的Sp和ΔSp沿試樣高度的分布。由圖8（a）和圖8（b）可知，凍融前Sp沿試樣高度分布較均勻（約16%），既表明凍融前的試樣較均勻，同時(shí)也間接證明了圖像處理方法的可行性。凍融后Sp較凍融前均增加，說(shuō)明經(jīng)過(guò)一次凍融后，土體產(chǎn)生了新的細(xì)觀孔隙。不補(bǔ)水條件下，鄰近最終凍結(jié)鋒面處（位于土層B）的孔隙率增加最多，而補(bǔ)水條件下卻發(fā)生在未凍區(qū)（即土層A），不難理解，不補(bǔ)水條件產(chǎn)生的原因是凍結(jié)過(guò)程土層B中形成了更多的冰晶和冰透鏡體，但對(duì)于補(bǔ)水條件，最可能的原因是凍結(jié)過(guò)程中產(chǎn)生的孔隙水壓力差，使未凍區(qū)發(fā)生滲流甚至管涌。

從圖8（c）可以發(fā)現(xiàn)，在凍結(jié)區(qū)中，距離冷端越遠(yuǎn)，ΔSp大體呈增加趨勢(shì)，這種現(xiàn)象是由于單向凍結(jié)過(guò)程中，距離冷端越遠(yuǎn)，凍結(jié)鋒面的發(fā)展速率會(huì)減小，土中的孔隙水有更充足的時(shí)間向凍結(jié)鋒面遷移，形成冰晶和冰透鏡體。除鄰近最終凍結(jié)鋒面的區(qū)域，補(bǔ)水條件下ΔSp明顯大于不補(bǔ)水條件，凍結(jié)區(qū)中產(chǎn)生的原因是補(bǔ)水條件下凍結(jié)階段土體內(nèi)部形成了更多的冰晶，融化后孔隙又無(wú)法恢復(fù)，凍融后含水率的增加可為此提供證據(jù)（圖6（a）），而未凍區(qū)中產(chǎn)生的原因是補(bǔ)水條件下凍結(jié)階段，外界水在負(fù)孔隙水壓力的作用下流入并通過(guò)未凍區(qū)，排開土顆粒，形成新的滲流通道。在鄰近最終凍結(jié)鋒面處，補(bǔ)水條件凍融后ΔSp相對(duì)較小，可能是因?yàn)槲磧鰠^(qū)和最終凍結(jié)鋒面所處土層之間的孔隙水壓力差引起管涌，進(jìn)而使未凍區(qū)的土顆粒發(fā)生遷移，監(jiān)測(cè)凍融過(guò)程中孔隙水壓力的變化已給出了論證[13]。此外，可發(fā)現(xiàn)凍融后沿試樣高度ΔSp的變化與宏觀孔隙率的變化存在不一致的情形，即凍融后Sp增大孔隙率卻減小，其原因只可能為孔隙分析尺度的差異，宏觀孔隙變化是細(xì)觀孔隙和微觀孔隙變化的綜合體現(xiàn)，細(xì)觀孔隙變化僅僅是宏觀孔隙變化中的一部分，下一階段將開展微觀孔隙變化的研究并為此提供佐證。

4.2 凍融前后平均孔徑分布

為定量表征土體凍融前后平均孔徑的變化，將圖像中細(xì)觀孔隙視為等效橢圓，將等效橢圓短軸的長(zhǎng)度視為細(xì)觀孔隙孔徑，并定義平均孔徑（DA）（即橫截面孔徑的平均值）及其凍融后變化量（ΔDA），比較凍融前后細(xì)觀孔徑沿試樣高度的變化，見式（4）。

式中，DA，b和DA，a同一試樣同一高度橫截面凍融前后的平均細(xì)觀孔徑。

由圖9（a）和圖9（b）可見，凍融前土體的DA沿試樣高度分布較均勻，但經(jīng)凍融后DA均增加，一方面是由于冰晶和冰透鏡體的形成和融化，產(chǎn)生了較大孔徑的細(xì)觀孔隙，另一方面細(xì)小的孔隙經(jīng)凍融后連接形成較大孔徑的細(xì)觀孔隙。與沿試樣高度Sp的分布類似，補(bǔ)水條件下凍融前后DA的最大差值發(fā)生在未凍區(qū)（即土層A），而不補(bǔ)水條件下則發(fā)生在鄰近最終凍結(jié)鋒面處。從圖9（c）可發(fā)現(xiàn)，補(bǔ)水條件下凍融后ΔDA顯著大于不補(bǔ)水條件，其沿試樣高度分布趨勢(shì)與ΔSp大致相同。

4.3 凍融前后分形維數(shù)分布

分形維數(shù)能夠定量表征事物的分形特征，研究表明土體的孔隙結(jié)構(gòu)也具有分形特征，利用分形維數(shù)能定量描述孔隙結(jié)構(gòu)的復(fù)雜程度，盒計(jì)數(shù)維數(shù)是一種最常用的分形維數(shù)，即用不同邊長(zhǎng)為r的正方形網(wǎng)格（盒子）覆蓋孔隙圖像，每次覆蓋孔隙占據(jù)的方格數(shù)計(jì)數(shù)為N（r），可用下式進(jìn)行線性回歸得到分形維數(shù)

式中：DF為分形維數(shù)；C為擬合常數(shù)。

基于以上分形維數(shù)的計(jì)算方法，對(duì)土體凍融前后典型橫截面進(jìn)行處理，圖10給出了有無(wú)補(bǔ)水條件凍融前后沿試樣高度分形維數(shù)的分布。由圖10可知，土體凍融前分形維數(shù)沿試樣高度的分布較均勻（約1.66），經(jīng)凍融后，沿試樣高度分形維數(shù)均發(fā)生不同程度的增加，說(shuō)明凍融使土體內(nèi)部的細(xì)觀孔隙結(jié)構(gòu)越來(lái)越復(fù)雜。此外，可以清晰發(fā)現(xiàn)，補(bǔ)水條件凍融后的分形維數(shù)和變化量明顯大于不補(bǔ)水條件，說(shuō)明凍結(jié)階段補(bǔ)水會(huì)顯著改變土體內(nèi)部的細(xì)觀孔隙結(jié)構(gòu)，進(jìn)而引起凍融后土體細(xì)觀結(jié)構(gòu)更復(fù)雜。

4.4 分形維數(shù)與細(xì)觀孔隙結(jié)構(gòu)參數(shù)之間的關(guān)系

分形維數(shù)可用于表征土體細(xì)觀孔隙結(jié)構(gòu)的復(fù)雜程度，而細(xì)觀孔隙的橫截面細(xì)觀孔隙率和平均孔徑均能體現(xiàn)土體的細(xì)觀孔隙結(jié)構(gòu)特征，結(jié)合圖8—圖10可以發(fā)現(xiàn)，凍融前后Sp和DA沿試樣高度的變化趨勢(shì)均與DF一致。圖11分別給出了DF與Sp、DA之間的關(guān)系，可以發(fā)現(xiàn)Sp、DA和DF之間均存在線性關(guān)系，表明分形維數(shù)可用于表征土體細(xì)觀孔隙的橫截面細(xì)觀孔隙率和平均孔徑，分形維數(shù)越大，土體細(xì)觀孔隙的橫截面孔隙率和平均孔徑越大。

5 結(jié)論

本研究以上海地區(qū)典型粉質(zhì)黏土為研究對(duì)象，結(jié)合凍融試驗(yàn)和X-CT試驗(yàn)，研究有無(wú)補(bǔ)水條件凍融前后飽和土體細(xì)觀結(jié)構(gòu)的變化，得到如下主要結(jié)論。

1）不補(bǔ)水條件下試樣發(fā)生了凍融頸縮現(xiàn)象，但補(bǔ)水條件下卻未發(fā)生，說(shuō)明不補(bǔ)水條件下土體暖端的體積收縮源于凍結(jié)過(guò)程中的水分遷移。

2）凍融后沿試樣高度CT灰度強(qiáng)度平均值的變化與土體物理指標(biāo)（即含水率、孔隙率和干密度）的變化之間存在良好的線性關(guān)系，因此凍融后CT灰度強(qiáng)度平均值的變化可用于反映凍融后土體物理特性的變化。

3）基于提出的土體CT圖像細(xì)觀孔隙識(shí)別方法，發(fā)現(xiàn)凍融后橫截面細(xì)觀孔隙率和平均孔徑沿試樣高度均發(fā)生不同程度的變化，且變化量沿試樣高度的變化趨勢(shì)相近，補(bǔ)水條件下最大變化發(fā)生在未凍區(qū)，而不補(bǔ)水條件下則鄰近最終凍結(jié)鋒面。

4）凍融后土體的分形維數(shù)均增大，且補(bǔ)水條件均大于不補(bǔ)水條件，說(shuō)明凍融使土體的細(xì)觀孔隙結(jié)構(gòu)變得更復(fù)雜，尤其在補(bǔ)水條件下。此外，分形維數(shù)與細(xì)觀孔隙的橫截面孔隙率和平均孔徑均呈良好的線性關(guān)系，分形維數(shù)越大，土體細(xì)觀孔隙率和平均孔徑越大。

【參 考 文 獻(xiàn)】

[1]熊輝，胡俊，周禹暄，等.圓形凍結(jié)板加固盾構(gòu)隧道端頭溫度場(chǎng)數(shù)值分析[J].森林工程，2021，37（6）：109-117.

XIONG H， HU J， ZHOU Y X， et al. Numerical analysis of temperature field at the end of shield tunnel reinforced with circular freezing plate[J]. Forest Engineering， 2021， 37（6）： 109-117.

[2]郜新軍，李銘遠(yuǎn)，張景偉，等.富水粉質(zhì)黏土中地鐵聯(lián)絡(luò)通道凍結(jié)法試驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào)，2021，40（6）：1267-1276.

GAO X J， LI M Y， ZHANG J W， et al. Field research on artificial freezing of subway cross passages in water-rich silty clay layers[J]. Chinese Journal of Rock Mechanics and Engineering， 2021， 40（6）： 1267-1276.

[3]YANG P， ZHAO J L， LI L. An artificial freezing technique to facilitate shield tail brush replacement under high pore-water pressure using liquid nitrogen[J]. KSCE Journal of Civil Engineering， 2021， 25（4）： 1504-1514.

[4]周潔，李澤垚，萬(wàn)鵬，等.組合地層滲流對(duì)人工地層凍結(jié)法及周圍工程環(huán)境效應(yīng)的影響[J].巖土工程學(xué)報(bào)，2021，43（3）：471-480.

ZHOU J， LI Z Y， WAN P， et al. Effects of seepage in clay-sand composite strata on artificial ground freezing and surrounding engineering environment[J]. Chinese Journal of Geotechnical Engineering， 2021， 43（3）： 471-480.

[5]FAN W H， YANG Z H， YANG P. A model for evaluating settlement of clay subjected to freeze-thaw under overburden pressure[J]. Cold Regions Science and Technology， 2020， 173： 102996.

[6]陳勇軍，章渺，代楊.凍融循環(huán)下石灰處治軟土動(dòng)態(tài)回彈模量試驗(yàn)研究[J].公路工程，2022，47（1）：79-84.

CHEN Y J， ZHANG M， DAI Y. Experimental study on dynamic resilient modulus of the lime modified soft soil under freeze-thaw cycle[J]. Highway Engineering， 2022， 47（1）： 79-84.

[7]SANTA G D， COLA S， SECCO M， et al. Multiscale analysis of freeze–thaw effects induced by ground heat exchangers on permeability of silty clays[J]. Geotechnique， 2019， 69（2）： 95-105.

[8]XU J， LI Y F， LAN W， et al. Shear strength and damage mechanism of saline intact loess after freeze-thaw cycling[J]. Cold Regions Science and Technology， 2019， 164： 102779.

[9]NISHIMURA S， OKAJIMA S， JOSHI B R， et al. Volumetric behaviour of clays under freeze–thaw cycles in a mesoscopically uniform element[J]. Geotechnique， 2021， 71（12）： 1150-1164.

[10]WANG S F， YANG Z H， YANG P. Structural change and volumetric shrinkage of clay due to freeze-thaw by 3D X-ray computed tomography[J]. Cold Regions Science and Technology， 2017， 138： 108-116.

[11]WANG S F， YANG P， YANG Z H. Characterization of freeze–thaw effects within clay by 3D X-ray Computed Tomography[J]. Cold Regions Science and Technology， 2018， 148： 13-21.

[12]蔡正銀，朱洵，黃英豪，等.凍融過(guò)程對(duì)膨脹土裂隙演化特征的影響[J].巖土力學(xué)，2019，40（12）：4555-4563.

CAI Z Y， ZHU X， HUANG Y H， et al. Influences of freeze-thaw process on evolution characteristics of fissures in expensive soils[J]. Rock and Soil Mechanics， 2019， 40（12）： 4555-4563.

[13]FAN W H， YANG P， YANG Z H. Freeze-thaw impact on macropore structure of clay by 3D X-ray computed tomography[J]. Engineering Geology， 2021， 280： 105921.

[14]ZHOU J， TANG Y Q. Experimental inference on dual-porosity aggravation of soft clay after freeze-thaw by fractal and probability analysis[J]. Cold Regions Science and Technology， 2018， 153： 181-196.

[15]王升福，楊平，劉貫榮，等.人工凍融軟黏土微觀孔隙變化及分形特性分析[J].巖土工程學(xué)報(bào)，2016，38（7）：1254-1261.

WANG S F， YANG P， LIU G R， 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.

[16]張英，邴慧，楊成松.基于SEM和MIP的凍融循環(huán)對(duì)粉質(zhì)黏土強(qiáng)度影響機(jī)制研究[J].巖石力學(xué)與工程學(xué)報(bào)，2015，34（S1）：3597-3603.

ZHANG Y， BING H， YANG C S. Influences of freeze-thaw cycles on mechanical properties of silty clay based on SEM and MIP test[J]. Chinese Journal of Rock Mechanics and Engineering， 2015， 34（S1）： 3597-3603.

[17]張澤，周泓，秦琦，等.凍融循環(huán)作用下黃土的孔隙特征試驗(yàn)[J].吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），2017，47（3）：839-847.

ZHANG Z， ZHOU H， QIN Q， et al. Experimental study on porosity characteristics of loess under freezing-thawing cycle[J]. Journal of Jilin University （Earth Science Edition）， 2017， 47（3）： 839-847.

[18]LU Y， LIU S H， WENG L P， et al. Fractal analysis of cracking in a clayey soil under freeze–thaw cycles[J]. Engineering Geology， 2016， 208： 93-99.

[19]KETCHAM R A， CARLSON W D. Acquisition， optimization and interpretation of X-ray computed tomographic imagery： applications to the geosciences[J]. Computers & Geosciences， 2001， 27（4）： 381-400.

[20]BRUN F， KOUROUSIAS G， DREOSSI D， et al. An improved method for ring artifacts removing in reconstructed tomographic images[J]. IFMBE Proceedings， 2009， 25： 926-929.

[21]FAN W H， YANG P， YANG Z H. Impact of freeze-thaw on the physical properties and compressibility of saturated clay[J]. Cold Regions Science and Technology， 2019， 168： 102873.