堆積碎石土中細小黏粒的細觀孔隙特征

董 輝，李智飛

全球氣候異常變化使地質(zhì)災(zāi)害的防治面臨著更大的挑戰(zhàn)。近3a的地質(zhì)災(zāi)害統(tǒng)計分析表明，我國的地質(zhì)災(zāi)害類型以滑坡、崩塌和地面塌陷為主，約占總體85.3%。自然因素引發(fā)（以降雨為主）的地質(zhì)災(zāi)害是其中的主體，如2013年1—4月全國552起地質(zhì)災(zāi)害中，降雨引發(fā)460起，占總數(shù)的83.3%。降雨誘發(fā)型滑坡多為淺層堆積滑坡，堆積物一般表現(xiàn)為風化卸荷、殘坡積和沖洪積等復雜成因形成的第四紀土與石塊的二重介質(zhì)混合體。由于土石體物質(zhì)組成的復雜性、結(jié)構(gòu)分布的不規(guī)則性給其工程性質(zhì)與力學性質(zhì)的研究帶來了較大的困難。而降雨誘發(fā)淺層堆積體滑坡，既表現(xiàn)為雨中和雨后短歷時的失穩(wěn)，也存在反復干濕交替作用下的長時累積效應(yīng)失穩(wěn)。針對降雨對堆積坡體的影響，許多學者從宏觀角度研究了降雨作用下碎石土（土石體）滑坡解體破壞機理［1－3］。

碎石土滲透性與其孔隙特征的相關(guān)性顯著，其中碎塊石粒組的含量和以粉粒、黏粒為主的細粒土粒組的含量對碎石土滲透系數(shù)影響最大［3］。研究其孔隙特征（孔隙變異度、孔隙復雜度和孔隙度等）對淺層堆積體雨水入滲有著重要的意義。至1973年Tovey［4］等利用電子顯微試驗結(jié)果對土微觀結(jié)構(gòu)進行定量化分析嘗試以后，國內(nèi)外學者在巖土體微觀孔隙結(jié)構(gòu)方面開展了大量的研究工作，獲得了卓有成效的研究成果?；谏皫r的微觀孔隙結(jié)構(gòu)研究，證明了砂巖孔隙具有分形特性［5］。針對黃土、膨脹土浸水前后微觀結(jié)構(gòu)，分析了濕陷的過程中孔隙的變化，獲得浸水后大孔徑孔隙含量減少，中、小孔徑的孔隙含量增加的結(jié)論［6］。在巖土微觀孔隙結(jié)構(gòu)的分形特征領(lǐng)域，驗證了黏性土微觀結(jié)構(gòu)中土顆粒的分布符合分形特征，并指出分形維數(shù)在1～2之間［7－9］。同時，基于描述顆粒和孔隙分布特性的顆粒分維數(shù)、孔隙分維數(shù)和譜維數(shù)，總結(jié)出多孔介質(zhì)滲透系數(shù)的一種預測公式［10］。隨著研究方法改進，以掃描電鏡（scanning electron microscope，SEM）、X 線 計 算 機 斷 層 攝 影 （computed tomography，CT）等為手段分析了包括重塑黏性土的孔隙結(jié)構(gòu)SEM圖像特征［11］，以及溶蝕作用下土體孔隙變化規(guī)律［12］的微細觀研究成果。

已有研究主要集中在一般土體孔隙的固有特征，而對于包含有土石間的空隙，碎石與土體自身的孔隙，以及在坡積和聚集作用下存在碎石土質(zhì)量密度、碎石含量和土石體壓實度空間變異性的土石混合介質(zhì)，目前在影響滲透性的碎石土中細小黏粒的孔隙特征描述方面還相對缺乏，同時，考慮降雨入滲的影響，研究不同含水率水平下的孔隙特征也有利于從細觀角度支撐已有降雨入滲過程中和入滲后滑坡解體破壞機理。本研究針對連續(xù)強降雨雨后殘坡積（堆積）碎石土滑坡失穩(wěn)破壞的災(zāi)害問題，以湖南省湘潭市昭山區(qū)某失穩(wěn)斜坡為研究對象，試圖采用高倍電子光學顯微鏡及數(shù)字圖像技術(shù)對堆積碎石土中細小黏粒的變異度和復雜度孔隙特征進行細觀定量研究，并分析含水率與孔隙特征的相關(guān)關(guān)系，為類似降雨入滲引起滑坡失穩(wěn)機理研究提供理論支持。

1 材料與方法

1.1 原材料選取與制樣

對已有碎石土研究成果，統(tǒng)計其碎石土的級配特征，將具有較好代表性的湖南省昭山區(qū)某斜坡堆積碎石土作為試驗原材料，其天然狀態(tài)含石量為69.3%，含黏粒量為21.7%，含水率為8.7%，母巖為粉質(zhì)砂巖，黏粒為紅黏土，圖1為取土坡體，材料級配曲線如圖2所示。將原狀土取回實驗室，經(jīng)烘干、篩分后，選取粒徑小于2mm的細小黏粒作為試驗用料，考慮到土體孔隙特征與含水率有關(guān)，試驗中按含水率不同分別配制5%，7.5%，10%，12.5%，15%和20%共6組試樣，每組制備5個，共30個試樣，如圖1所示。試樣均用39.1mm的套筒壓實成形，高度為4cm，土粒相對密度Gs為2.72，干密度ρ為1.75g／cm3。

圖1 采樣坡體和試驗樣品

圖2 碎石顆粒級配曲線

1.2 孔隙細觀形態(tài)的獲取與處理

為了獲得土樣孔隙的細觀特征，采用奧林巴斯高倍光學顯微鏡（Olympus－BX51M）放大200倍后對切片表面進行觀察，拍攝圖像（像元）面積為1 280×1 024。

應(yīng)用圖像處理技術(shù)對拍攝的圖像進行處理，包括去除背景（消除噪聲和部分陰影）、灰度化（去除色彩信息）、高斯模糊（增強圖像中目標信息）、二值化（分割圖象為目標和背景）等，從圖形中提取孔隙的面積、周長、圓度、長軸、角度和孔隙度等基本參數(shù)。面積A表示二值圖像中區(qū)域標記為“l(fā)”的像素個數(shù)；周長等于π乘以廣義方向直徑的平均數(shù)P＝πDm，二值圖像中Dm的計算公式為：

1.3 孔隙特征的研究方法

孔隙定向性是由測量對象定向角（長軸與水平線的夾角）的分布θ1，θ2，…，θn，用一定的角密度 Δθ，一般將180°分成18個定向角區(qū)間，計算定向角落入每個區(qū)間的頻率［13－14］。用F（α）來代表各定向角區(qū)間的結(jié)構(gòu)單元體或孔隙的頻率，計算公式為：

式中：nα，n——定向角落入［θn，θn＋1］的測量對象數(shù)量和總的測量對象數(shù)量?？紫兜亩ㄏ蛐杂绊懼謧?cè)向入滲遷移，定向角在各區(qū)間分布越不均勻，越有利于水分的水平滲透遷移。定向角頻率統(tǒng)計結(jié)果詳見表1，定向頻率可用玫瑰圖如圖3所示。

表1 試樣的定向角頻率統(tǒng)計結(jié)果

孔隙變異度E用土壤切片上、下相對表面的孔隙度P之間的相對差異作為評價指標?？紫蹲儺惗菶表達式為［15］：

式中：P上，P下——各切片的上、下表面的孔隙度。E值越高，說明孔隙變異度越大。

土體實際上是具有統(tǒng)計意義上的自相似的分形結(jié)構(gòu)特征［14］，采用統(tǒng)計自相似的方法來定量的描述復雜的孔隙分布特征，可以揭示土體孔隙特征對雨水入滲的影響?？紫稄碗s度可以用孔隙維數(shù)D來描述，即利用孔隙的周長與面積之間的關(guān)系來定量評價［15］，計算公式為：

式中：D，P，A——孔隙分數(shù)維數(shù)、孔隙周長和孔隙面積（孔隙所占像元數(shù)）。D值越高，說明孔隙的扭曲程度越高，孔隙就越復雜。需要指出的是，在此研究中D值反映的是整幅圖像中孔隙的復雜度，不具體涉及圖像中單個的孔隙。

2 孔隙特征分析

2.1 含水率對孔隙度的影響

土體的孔隙度定義為圖像處理二值化后，黑色面積與圖像總面積的百分比，其數(shù)值大小影響著土體的結(jié)構(gòu)及其密實度、強度等物理力學性質(zhì)。對不同含水率試樣的圖像進行處理，得到各含水率條件下的土體的孔隙度。觀察發(fā)現(xiàn)，當含水率（20%）較高時，試樣表面已形成水膜，很難用顯微鏡準確觀測到切片表面孔隙結(jié)構(gòu)，因此只考慮含水率在15%以下的情況。

表2為不同含水率條件下試樣的孔隙度。分析表2可知：（1）不同含水率試樣切片面的孔隙度比下表面的孔隙度高，表明圓柱土樣由上至下越來越密實，這是分層擊實的結(jié)果；（2）不同含水率下試樣的孔隙度值位于30%左右，含水率不同的土體，盡管初始孔隙比相同，在同等外力作用下所產(chǎn)生的壓縮變形有很大差別，一般在最優(yōu)含水量附近隨著壓實含水量的增大，土體的壓縮變形也增加，孔隙度減??；（3）土體壓縮后的孔隙度與含水率關(guān)系密切，土體的孔隙度隨著含水率的增加呈現(xiàn)先增加后減小的趨勢。為了更好地描述在不同含水率下孔隙度的變化特征，對其數(shù)據(jù)進行了多項式擬合，表達式為：

式中：P——孔隙度；εa——含水率。由擬合公式可以得出在含水率為10.41%時土體孔隙度出現(xiàn)峰值為33.52%，表明在此含水率土體還可以繼續(xù)壓實達到最佳密實狀態(tài)。

表2 試樣的孔隙度統(tǒng)計結(jié)果

圖3 不同含水率的定向性玫瑰圖

2.2 孔隙定向性分析

從表1和圖3可見：（1）不同含水率下，孔隙定向角在各個區(qū)間都有不同幾率的分布，且在10°～20°，50°～60°，100°～110°和140°～150°這4個區(qū)間內(nèi)分布相對密集，其中含水率為15%時最大，定向角頻率達31.7%；（2）隨著含水率的增加，孔隙定向性由發(fā)散趨于定向然后又轉(zhuǎn)向發(fā)散，含水率為7.5%和12.5%時，孔隙定向角在0°～20°，40°～60°兩個區(qū)間內(nèi)的分布相對較多，累計分別達28.8%和33.3%，表現(xiàn)出較明顯的定向性，其他區(qū)間內(nèi)分布比較均勻；當含水率為10%時，孔隙定向角在0°～20°，40°～60°，80°～110°和130°～150°這4個區(qū)間內(nèi)定向性明顯表現(xiàn)突出；當含水率為15%時，孔隙定向角在各個區(qū)間都有不同幾率的分布，分布相對不均勻，輻射狀最為明顯，利于水分的水平滲透遷移。分析其原因，主要是隨著含水率的增加，孔隙度先增加后減少，且大孔隙相對增加，而多孔介質(zhì)的滲透性與孔隙度之間存在非線性關(guān)系，在相同孔隙度的情況下孔徑越大則多孔介質(zhì)的滲透性越好。

2.3 孔隙變異度與復雜度分析

表3和圖4為各含水率下試樣的孔隙數(shù)目、孔隙變異度和復雜度的統(tǒng)計結(jié)果。

（1）土體孔隙數(shù)目隨著含水率的增加而減少，這是由于土體結(jié)構(gòu)是三相體系，相同體積相同質(zhì)量的土顆粒下，隨著含水率的增加土顆粒之間的膠結(jié)作用增強，從而減少了孔隙的數(shù)目。

（2）土體孔隙變異度數(shù)值位于4.19～5.16之間，且隨著含水率的增加呈現(xiàn)出減小的趨勢；含水率在7.5%～12.5%之間，變異度對含水率不敏感，其原因是含水率為10%左右，土體的壓實性好；在含水率12.5%處，孔隙變異度下降曲線出現(xiàn)明顯轉(zhuǎn)折，斜率變小，表明含水率超過12.5%后，土顆?？伤苄栽鰪姡馏w軟化程度加劇，壓縮系數(shù)變大。

（3）土樣切片的孔隙復雜度值位于1.625～1.675之間，變化區(qū)間不大，只有0.05，表明含水率對孔隙復雜度的影響不大；不同含水率下，土樣的孔隙復雜度隨著含水率的增加而先降低后趨于平緩，這是由于隨著含水率的增加，顆粒團?；觿?，孔隙數(shù)目減少，且大孔隙數(shù)量相對較多，微觀孔隙變得單調(diào)、簡單，因而孔隙復雜度減小，即含水率越小，孔隙的扭曲程度越高，孔隙就越復雜，且孔隙數(shù)含量多，越有利于表層雨水下滲。

（4）一般孔隙復雜度愈大，滲透性愈好，這是含水率越小，孔隙的扭曲程度越高，孔隙連通性大，且吸附在土顆粒表面上的束縛水越少，液體流動的阻力也隨之減小的緣故；孔隙的變異度實質(zhì)是孔隙度在垂直空間上的變化程度，一般滲透性隨有效孔隙度的增加而增強，即孔隙度越高，且連通性好，流體滲透能力愈好。

表3 試樣孔隙特征統(tǒng)計數(shù)據(jù)

圖4 孔隙特征與含水率相關(guān)曲線

3 結(jié)果討論

土體的孔隙作為滲流通道，其孔隙的大小及其連通性、曲折性等因素直接影響土體的滲透性，即土體的滲透特性主要取決于土的孔隙度和孔隙尺度分布以及孔隙微細觀特征［16］。細小土顆粒作為堆積碎石土的重要組成部分，其細觀孔隙結(jié)構(gòu)影響著土石體的滲透性能及其雨水的滲透運移。

雨水入滲實質(zhì)上是水分在土體飽氣帶中的運動，是一個涉及兩相流的過程，即水在下滲過程中驅(qū)替空氣的過程。一般情況下飽和度愈低，滲透性愈差，這是因為低飽和土的孔隙中存在較多氣泡會減小過水斷面，甚至堵塞細小孔道?？紫抖仁强紫督Y(jié)構(gòu)的重要特征之一，土壤孔隙度隨著含水率的變化而變化，由不同含水率與孔隙度的二次擬合結(jié)果可以看出，孔隙度隨著含水率的增加先增加后減少?？紫抖仍黾邮雇恋倪^水斷面增加，但土體的滲透性不一定增強，滲透性還與孔徑的大小、孔隙復雜度、孔隙定向性等相關(guān)。當孔隙度相同的條件下，孔徑越大，滲透性越好。由圖3可知，當含水率為15%時，孔隙定向角分布相對不均勻，輻射狀最為顯著，有利于水分的水平滲透遷移?？紫兜淖儺惗葘嶋H是孔隙度在垂直空間上的變化程度，孔隙復雜度與孔隙自身的二維面積和周長相關(guān)，表征了孔隙自身的復雜程度，孔隙的扭曲程度越高，孔隙就越復雜，且孔隙含量多，越有利于表層雨水下滲。由表3可知，不同含水率下，孔隙數(shù)目位于206～456之間，孔隙變異度位于4.19～5.16之間，孔隙復雜度位于1.625～1.675之間，變化區(qū)間不大，只有0.05，說明含水率對孔隙數(shù)目影響較大，對孔隙的變異度和復雜度影響相對不明顯；由圖4可見，孔隙數(shù)目、孔隙復雜度以及孔隙變異度都隨著含水率的增加而減小，表明含水率越小，吸附在土顆粒表面上的束縛水越少，液體流動的阻力愈小，孔隙大小分布分形維數(shù)越大，孔隙占據(jù)的面積越大，孔隙數(shù)目越多，更易入滲。

4 結(jié)論

（1）相同質(zhì)量的細小黏粒土體的孔隙度與含水率存在二次曲線關(guān)系，且得出試驗土體在含水率為10.41%壓實性最好?？紫抖仍黾油馏w的過水斷面，但碎石土中細小黏粒的滲透性增強，還與孔徑的大小、孔隙復雜度、孔隙定向性等相關(guān)。

（2）細小黏粒的微觀結(jié)構(gòu)分形維數(shù)與其結(jié)構(gòu)具有密切的關(guān)系，含水率越高，試樣的團?；礁?，則孔隙數(shù)目及其分形維數(shù)越小，孔隙復雜度越低。

（3）孔隙數(shù)目、孔隙復雜度以及孔隙變異度隨著含水率的增加而減小，表明含水率越小，孔隙的扭曲程度越高，孔隙就越復雜，且孔隙數(shù)含量越多，越有利于表層雨水下滲。

（4）不同含水率下孔隙定向角在10°～20°，50°～60°，100°～110°和140°～150°這4個方位上分布相對密集，且含水率越小，孔隙定向性相對較明顯；統(tǒng)計出孔隙方向分布具有輻射狀，含水率為15%最為明顯，有利于水分的水平滲透遷移。

［1］ 許建聰，尚岳全.降雨作用下碎石土滑坡解體變形破壞機制研究［J］.巖土力學，2008，29（1）：106－112，118.

［2］ 寧萬輝，寧健，俞美華，等.降雨對碎石土邊坡穩(wěn)定性的影響分析［J］.水電能源科學，2011，29（1）：83－85.

［3］ 許建聰，尚岳全.碎石土滲透特性對滑坡穩(wěn)定性的影響［J］.巖石力學與工程學報，2006，25（11）：2264－2271.

［4］ Tovey N K.Quantitative analysis of electron micrographs of soil structure［C］∥Minter national Symposium on Soil Structure.Stockholm：Swedish Geotechnical Institute，1973.

［5］ Katz A J，Thompson A H.Fractal sandstone Pores：Implications for conductivity and Pore formation［J］.Phys.Rev.Lett.，1985，54（12）：1325－1328.

［6］ 方祥位，申春妮，汪龍，等.Q2黃土浸水前后微觀結(jié)構(gòu)變化研究［J］.巖土力學，2013，34（5）：1319－1324.

［7］ 王寶軍，施斌，劉志彬，等.基于GIS的黏性土微觀結(jié)構(gòu)的分形研究［J］.巖土工程學報，2004，26（2）：244－247.

［8］ 黃麗.飽和軟黏土微觀孔隙的定量分析及其分形研究［D］.湖北 武漢：武漢理工大學，2007.

［9］ 許勇，張季超，李伍平.飽和軟土微結(jié)構(gòu)分形特征的試驗研究［J］.巖土力學，2007，28（S）：49－52.

［10］ Shi Minglun，Cheng Yongping.Determination of permeability using fractal method for porous media［J］.Science in China（E）：Techological Sciences，2001，44（6）：625－630.

［11］ 熊承仁，唐輝明，劉寶琛，等.利用SEM照片獲取土的孔隙結(jié)構(gòu)參數(shù)［J］.中國地質(zhì)大學學報，2007，32（3）：415－419.

［12］ 吳恒，張信貴，易念平.水土作用與土體細觀結(jié)構(gòu)研究［J］.巖石力學與工程學報，2000，19（2）：199－204.

［13］ 施斌，李生林.黏性土微觀結(jié)構(gòu)SEM圖像的定量研究［J］.中國科學（A）：數(shù)學，1995，25（6）：666－672.

［14］ 趙安平，王清，李楊，等.長春季凍區(qū)路基土微觀孔隙特征的定量評價［J］.工程地質(zhì)學報，2008，16（02）：233－238.

［15］ 李德成，Velde B，張?zhí)伊?利用土壤切片的數(shù)字圖像定量評價土壤孔隙變異度和復雜度［J］.土壤學報，2003，40（5）：678－682.

［16］ 閆小慶，房營光，張平.膨潤土對土體微觀孔隙結(jié)構(gòu)特征影響的試驗研究［J］.巖土工程學報，2011，33（8）：199－204.