膨脹巖裂隙發(fā)育的現(xiàn)場觀測及描述方法研究

王 軍，龔壁衛(wèi)，張家俊，周小文，譚峰屹

（1.華南理工大學(xué)亞熱帶建筑國家重點實驗室，廣州 510640；2.長江科學(xué)院水利部巖土力學(xué)與工程重點實驗室，武漢 430010）

膨脹巖裂隙發(fā)育的現(xiàn)場觀測及描述方法研究

王 軍1，2，龔壁衛(wèi)2，張家俊1，2，周小文1，譚峰屹2

（1.華南理工大學(xué)亞熱帶建筑國家重點實驗室，廣州 510640；2.長江科學(xué)院水利部巖土力學(xué)與工程重點實驗室，武漢 430010）

基于圖像灰度化及二值化的數(shù)值處理技術(shù)，提出了膨脹巖裂隙度的概念及定量描述方法。在南水北調(diào)中線新鄉(xiāng)潞王墳試驗段進(jìn)行了膨脹巖裂隙觀測試驗，定量分析了膨脹巖邊坡在干濕循環(huán)中的裂隙的變化規(guī)律及其與含水率的相關(guān)關(guān)系。分析表明，膨脹巖裂隙的發(fā)育與天氣變化密切相關(guān)，隨著干濕循環(huán)的發(fā)展，裂隙的寬度和裂隙度逐漸增加，利用Matlab數(shù)字圖像處理方法得到的裂隙率可以較好地描述裂隙的發(fā)育和發(fā)展情況。

膨脹巖（土）；南水北調(diào)中線工程；Matlab；裂隙度

1 概 述

膨脹巖（土）是自然地質(zhì)過程中形成的多裂隙且具有顯著脹縮性的巖體。由于膨脹巖（土）的黏土礦物成分主要為強(qiáng)親水性的蒙脫石或伊利石，因而具有吸水膨脹、失水收縮和反復(fù)脹縮變形的特征，加之分布于土體中的裂隙和層間結(jié)合面，使其對各種工程結(jié)構(gòu)物具有嚴(yán)重的破壞作用［1］，如地基隆起、路基開裂以及邊坡失穩(wěn)等。研究表明：膨脹巖（土）邊坡中的裂隙及其發(fā)展對邊坡的穩(wěn)定起著關(guān)鍵性的作用［2］。大氣降雨和地表蒸發(fā)，使得膨脹巖（土）濕脹干縮，并循環(huán)往復(fù)，逐步使土體產(chǎn)生裂隙，而裂隙的存在又進(jìn)一步為雨水的入滲和水分的蒸發(fā)提供了通道，加快了土中含水量的波動和脹縮的反復(fù)發(fā)生，導(dǎo)致裂隙擴(kuò)大并向土的深部發(fā)展。裂隙的產(chǎn)生、延伸、擴(kuò)展、貫通，破壞了膨脹巖（土）的結(jié)構(gòu)性和整體強(qiáng)度，使得膨脹巖（土）滑坡常常成為難以治理的“頑疾”。

有關(guān)膨脹巖（土）裂隙的研究工作，目前主要集中在兩個方面，一些主要采用當(dāng)前工程力學(xué)領(lǐng)域的成就（如彈性理論、斷裂力學(xué)）來分析膨脹巖（土）裂隙的分布、深度、開展和延伸度，以及裂隙及其周圍土體的強(qiáng)度和變形特性［3，4］，這種方法包含了一定的假定和簡化條件，成果的精度尚難以驗證。另一些少量的定量研究也只是在室內(nèi)針對重塑土開展的［5－8］。由于室內(nèi)試驗中裂隙的發(fā)育與現(xiàn)場實際情況可能有較大的差異，室內(nèi)試驗未必能夠反映現(xiàn)場的復(fù)雜情況，已有的關(guān)于裂隙的定量描述［2，5，6］僅是假借其它的物理量來衡量裂隙度，如灰度熵［5］，分形維數(shù)［6］等，這使得對于膨脹巖（土）的裂隙特性乃至邊坡失穩(wěn)機(jī)理仍沒有得到令人信服的解釋。因此，繼續(xù)就膨脹巖（土）的裂隙問題開展定量的試驗研究，尤其是采取現(xiàn)場試驗觀測手段取得實際的裂隙數(shù)據(jù)，具有重要的意義。為此，從2007年底開始，長江科學(xué)院等有關(guān)單位在南水北調(diào)中線工程新鄉(xiāng)段開展了膨脹巖邊坡破壞機(jī)理及加固方法的試驗研究，針對膨脹巖的裂隙發(fā)育規(guī)律進(jìn)行了觀察、量測和分析研究。

2 裂隙率的表達(dá)方法

袁俊平［5］等利用遠(yuǎn)距光學(xué)顯微鏡對重塑膨脹土在自然風(fēng)干條件下的裂隙發(fā)展變化進(jìn)行了動態(tài)的、定量的觀測，其試驗結(jié)果表明，裂隙圖像的灰度熵可以很好地表征裂隙的發(fā)育發(fā)展程度，適合作為膨脹土裂隙發(fā)育程度的度量指標(biāo)。但是灰度熵實際上反映的是圖片灰度分布的均勻與離散情況，在表征裂隙發(fā)育程度上具有一定的局限性，灰度熵能否直接用作表征裂隙度的變量有待進(jìn)一步研究，而且其容易受亮度和對比度的影響，精度上取決于圖像資料的后處理精度，之前的研究者對此缺乏有效的處理手段。

易順明、徐永福［6，7］等指出，膨脹土裂隙幾何圖像的不規(guī)則程度高，結(jié)構(gòu)復(fù)雜，定量描述極為困難，建議采用分形理論進(jìn)行統(tǒng)計。對膨脹土不同尺度的裂隙，在二維平面上均構(gòu)成為網(wǎng)絡(luò)系統(tǒng)，其分形是建立在一種統(tǒng)計分布上，可用盒維數(shù)法計算其分維值，分維值能表征裂隙發(fā)育的程度。但分維值是建立在膨脹土的土顆粒大小的分布與孔隙分布的基礎(chǔ)之上的，同樣，表征裂隙發(fā)育的程度不夠直觀和直接。

為此，本文提出一種結(jié)合Matlab數(shù)值圖像處理功能的圖像數(shù)值化分析方法，即將現(xiàn)場拍攝的常規(guī)RGB照片進(jìn)行處理降噪，提高圖像數(shù)值的精度，并將此數(shù)值圖像進(jìn)行裂隙骨骼化，進(jìn)而將目標(biāo)區(qū)域的裂隙數(shù)值化的裂隙統(tǒng)計方法。此方法在復(fù)雜的現(xiàn)場具有較強(qiáng)的適用性，而且Matlab圖像處理技術(shù)保證了計算裂隙度的精度，因此，該處理方法不僅能保證裂隙數(shù)值圖像的精度，而且克服了以往的研究者需要人為調(diào)整外界干擾因素的影響。另外，該處理方法直接將裂隙數(shù)值化，裂隙的發(fā)展變化直接反映在裂隙數(shù)值圖像上，通過對裂隙圖像的數(shù)值處理可以得到裂隙的發(fā)展變化。圖像數(shù)值化的實現(xiàn)過程如下。即原始圖像→灰度圖像→裂隙輪廓圖像→裂隙輪廓二值化圖像。

裂隙面拍照采集到的彩色圖像是RGB（Red／Green／Blue）三基色模型，不僅包含亮度信息，還包含顏色信息（圖1（a）），顏色信息包含著分析裂隙面不需要的信息，甚至是環(huán)境影響的雜質(zhì)信息，需要消除這些影響，圖像就需要進(jìn)一步處理。

利用Mat lab圖像處理工具，將采集到的彩色圖像轉(zhuǎn)化為灰度圖，灰度圖像只表示圖像的亮度信息而沒有顏色信息，因此，灰度圖像的每個像素點上只包含一個量化的灰度級（即灰度值），用來表示該點的亮度水平（圖1（b）），同樣的裂隙面不同的環(huán)境拍出的照片灰度值會不一樣，需要將圖像進(jìn)一步處理。

將灰度圖進(jìn)行邊緣檢測處理，提取出裂隙的邊緣，裂隙邊緣是以圖像的局部特征不連續(xù)的形式出現(xiàn)的，即是指圖像局部亮度變化最顯著的部分，邊緣檢測的基本思想是檢測圖像中的裂隙邊緣點，再按照一定的方式將邊緣點連接成輪廓（圖1（c）），這種處理方法消除了不同環(huán)境造成圖片灰度的不同的影響。

為了得到裂隙率，還需要把處理過的灰度圖像轉(zhuǎn)換為二值圖像，如圖1（d）所示。

二值圖像的像素值只有“0”和“1”兩種，非裂隙點的值為“1”，裂隙點的值為“0”。通過統(tǒng)計二值圖像中像素值為“0”的點的個數(shù)，除以圖像總的像素點數(shù)，即可得到該圖像的裂隙率，作為表征裂隙發(fā)育程度（裂隙度）的物理量。裂隙率的表達(dá)式為

式中：D為裂隙率；Σ0為像素值為“0”的總點數(shù)；Σ1為像素值為“1”的總點數(shù)。

以上圖像處理的方法通過在Mat lab中編制程序得以實現(xiàn)。

圖1 裂隙率在M at lab中的實現(xiàn)過程Fig.1 The processing diagram of calculating the crack rate in M at lab

3 試驗方案及方法

在南水北調(diào)中線工程新鄉(xiāng)段，開展了膨脹巖邊坡的現(xiàn)場試驗。根據(jù)取樣試驗結(jié)果，其基本的物理性質(zhì)如下：天然含水率為21.4%；天然干密度為1.65 g／cm3；粒徑≥0.075 m的砂粒含量為22.3%，0.075～0.005 mm的粉粒含量為41.9%，■0.005 mm的黏粒含量為35.8%；其液限、塑限分別為54.9%和25%。X衍射試驗分析結(jié)果表明，土樣的主要黏土礦物成分為蒙／伊混合層（8%）和伊利石（14%），自由膨脹率60.4。根據(jù)國標(biāo)GB112－87膨脹性分類標(biāo)準(zhǔn)，該土樣屬于弱膨脹土。

現(xiàn)場試驗在渠道邊坡位置開挖出一個3 m×3 m的直立剖面，裸露出邊坡垂直臨空面的裂隙作為觀測斷面，如圖2（a，b）所示。在距離觀測斷面1.5 m處，設(shè)立一固定支架，作為數(shù)碼相機(jī)定點攝錄裂隙照片的平臺（如圖2（c）所示）。同時，在觀測斷面右側(cè)，沿不同深度，每間隔0.2 m埋設(shè)ML2X含水率探頭，每天定時觀測不同深度下土體含水率的變化，如圖3所示。在斷面剛開挖后，有數(shù)條裂隙已經(jīng)完全貫通整個斷面，選取其中的3條裂隙作為主裂隙，選取局部貫穿的2條裂隙作為次裂隙，對這5條裂隙分別編號，量測其長度，并沿裂隙開展方向測量不同位置的裂隙寬度。5條裂隙的形態(tài)如圖4所示。

圖2 觀測平面Fig.2 Observation ichnography

圖3 含水量探頭布置Fig.3 M oisture probe layout

圖4 裂隙分布Fig.4 Cracks distribution diagram

對于裂隙平面，每2 d定點定時拍照，記錄裂隙圖像，從2008年11月23日至2009年4月1日，直至3月份雨季到來，觀測斷面被雨水沖垮，觀測終止，共計130 d，經(jīng)歷2次干濕循環(huán)。

4 試驗成果及討論

4.1 裂隙寬度的變化

1＃主裂隙位于坡面表層，產(chǎn)狀為鉛直向；4＃主裂隙位于坡面表層，產(chǎn)狀為順坡向；3＃主裂隙位于坡面深層，產(chǎn)狀為逆坡向。

2＃，5＃次裂隙均為1＃主裂隙的分支。

采用直尺定點測量裂隙的寬度，其中1＃，3＃，4＃主裂隙的初始寬度分別為10，7.7，10 mm；2＃，5＃次裂隙初始寬度分別為7，3 mm。主裂隙和次裂隙的寬度變化曲線分別見圖5和圖6，相應(yīng)的含水率的變化如圖7所示。

圖5 觀測平面主裂隙寬度變化曲線Fig.5 M ain-cracksw idth change curve in observation plane（unit in mm）

圖6 觀測平面次裂隙寬度變化曲線Fig.6 Sub-cracksw idth change curve in observation p lane（unit in mm）

圖7 觀測平面含水率變化曲線Fig.7 Soilwater content change curve in observation p lane

4.1.1 主裂隙寬度變化

由圖5至圖7可以看出，裂隙寬度的變化實際上是隨著降雨-蒸發(fā)的循環(huán)而呈現(xiàn)收縮與張開的反復(fù)循環(huán)的過程。第一次降雨前，1＃裂隙的寬度由10 mm增大到14.3 mm，降雨后裂隙寬度減小到14 mm，并最終達(dá)到8.7 mm；第二次降雨后，1＃裂隙的寬度由8.7 mm增大到15.6 mm，3＃，4＃裂隙也有類似現(xiàn)象。由此可見，持續(xù)干燥的天氣會導(dǎo)致裂隙張開，寬度增大；降雨會導(dǎo)致裂隙寬度減小甚至閉合。對于貫穿整個裂隙面的主裂隙而言，裂隙寬度的變化相對平緩，3＃裂隙寬度甚至一度減少，這主要是由于3＃裂隙的產(chǎn)狀為逆坡向，降雨發(fā)生后，土體強(qiáng)度降低，后方土體主動位移影響使得裂隙閉合，寬度變小。1＃主裂隙產(chǎn)狀為順坡面，裂隙寬度的變化只受干濕循環(huán)的影響，降雨發(fā)生時，裂隙寬度變小，持續(xù)干旱，裂隙寬度持續(xù)增加。而4＃裂隙雖然與1＃裂隙一樣為順坡向裂隙，但由于其位于渠道坡面表層，降雨過后雨水下滲，導(dǎo)致裂隙閉合，裂隙寬度變小，其影響要大于應(yīng)力的影響，故而也一度呈現(xiàn)寬度下降的趨勢。在本次試驗當(dāng)中，經(jīng)過2次的干濕循環(huán)，1＃主裂隙寬度最大達(dá)到15.6 mm。3月份的雨季到來之后導(dǎo)致雨水大量入滲，隨即裂隙被土充填、愈合，觀測終止，如圖8所示。

圖8 觀測終止時的裂隙平面Fig.8 Crack observation plane graph as the observation finished

4.1.2 次裂隙寬度變化

次裂隙受干濕循環(huán)影響導(dǎo)致的寬度變化更為明顯（圖6）。第一次降雨前后，2＃裂隙的寬度由7 mm增大到11mm，隨后裂隙寬度減小到10mm；5＃裂隙的寬度由3 mm增大到5 mm，隨后減小到4 mm。第二次降雨前后，2＃裂隙的寬度由11 mm增大到14 mm，5＃裂隙的寬度由4 mm增大到7 mm，隨后由于大雨導(dǎo)致裂隙被充填，裂隙閉合。次裂隙沒有貫穿整個豎向的觀測斷面，其寬度變化反映的是局部土體受外界因素影響下的土體狀態(tài)的變化，其變化的狀態(tài)與土體含水量變化基本一致。這表明，干濕循環(huán)是控制土體裂隙的最重要的因素。在本次試驗中，由于次裂隙的持續(xù)開展，最大裂隙寬度也達(dá)到14 mm，多條次裂隙的開展導(dǎo)致次裂隙周圍的土體發(fā)生脫落或崩解，進(jìn)一步提供更多的雨水入滲通道，3月份雨季到來之后，裂隙被土充填、愈合，觀測終止。

4.2 觀測斷面裂隙的定量化描述

根據(jù)上述有關(guān)裂隙率的定義，將觀測斷面裂隙率與時間的關(guān)系描述如圖9。由圖7、圖9可以看出，裂隙率的變化趨勢與含水量變化基本是一致的，含水量的持續(xù)降低（2008年11月至2009年1月）會導(dǎo)致土體裂隙持續(xù)開展和發(fā)育，裂隙率不斷增加。同時，裂隙在初期發(fā)展較明顯，后期有所減緩（2009年1月至2月）。干濕循環(huán)會導(dǎo)致裂隙率呈現(xiàn)階梯性的增加，第一次干濕循環(huán)，裂隙率峰值增加到0.08，其峰值是起初始值的2倍；第二次干濕循環(huán)（2009年2月至4月），裂隙率在此基礎(chǔ)上持續(xù)增加，其峰值達(dá)到其初始值的3倍。直至3月份雨季到來之后，裂隙被填充，裂隙閉合，試驗終止。

圖9 裂隙率隨時間變化曲線Fig.9 Duration curve of crack rate

由上述分析，本文提出的利用Mat lab數(shù)字圖像處理方法得到的裂隙率的變化很好地描述了裂隙的發(fā)展變化情況，一次干濕循環(huán)導(dǎo)致裂隙率增加，在對裂隙的觀測上表現(xiàn)為裂隙條數(shù)的增多和裂隙寬度的增大，反復(fù)干濕循環(huán)導(dǎo)致裂隙率持續(xù)緩慢增加，裂隙的寬度持續(xù)增大。裂隙率實際上反映的是裂隙分割土體強(qiáng)度的一個有效指標(biāo)，裂隙率越大，土體被裂隙分割的強(qiáng)度越高。裂隙率對于描述裂隙的發(fā)展程度是十分有效的，合理地解釋了裂隙隨著氣候變化（干濕循環(huán)）而發(fā)育和發(fā)展。觀測終止時的裂隙面，如圖8所示，土體破碎十分劇烈，已經(jīng)呈現(xiàn)類沙土似的淺層滑坡。

5 結(jié)論及建議

（1）膨脹巖（土）裂隙的發(fā)育程度可以用裂隙度較好地描述，裂隙度可以通過將彩色圖像灰度化和二值化計算得到。

（2）膨脹巖（土）裂隙的發(fā)育與天氣變化密切相關(guān)，隨著干濕循環(huán)的發(fā)展，裂隙的寬度和裂隙度逐漸增加，加劇了雨水入滲和蒸發(fā)的程度，進(jìn)一步促進(jìn)裂隙的發(fā)育和膨脹巖（土）結(jié)構(gòu)破壞。

（3）利用Matlab數(shù)字圖像處理方法得到的裂隙率可以很好地描述裂隙的發(fā)育和發(fā)展情況，裂隙度增加，裂隙分隔土體的程度增加，土體完整性降低，粘聚力減小，土體強(qiáng)度降低，邊坡滑坡呈現(xiàn)類沙土的淺層滑坡。

［1］ 劉特洪.工程建設(shè)中的膨脹土問題［M］.北京：中國建筑工業(yè)出版社，1997.（LIU Tei-hong.Expansive Soil Problem in Construction［M］.Beijing：China Architecture and Building Press，1997.（in Chinese））

［2］ 包承綱.膨脹土裂隙性研究膨脹土處治理論、技術(shù)與實踐［M］.北京：人民交通出版社，2004.（BAO Chengang.Study on Expansive Soil Cracks Treatment of Expansive Soil，Theory，Technology and Practice［M］.Beijing：People＇s Communications Press，2004.（in Chinese））

［3］ KONRAD JM，AYAD R.An idealized gramework for the analysis of cohesive soils undergoing desiccation［J］.Can Geotech，1997，34：477－488.

［4］ 鄭少河，金劍亮，姚海林，等.地表蒸發(fā)條件下的膨脹土初始開裂分析［J］.巖土力學(xué)，2006，27（12）：2229－2233.（ZHENG Shao-he，JIN Jian-liang，YAO Hai-lin，et al.Analysis of initial cracking behavior of expansive soil due to evaporation［J］.Rock and Soil Mechanics，2006，27（12）：2229－2233.（in Chinese））

［5］ 袁俊平，殷宗澤.膨脹土裂隙的量化指標(biāo)與強(qiáng)度性質(zhì)研究［J］.水力學(xué)報，2004，（6）：108－113.（YUAN Jun-ping，YIN Zong-ze.Quantitative index of fissure and strength characteristics of fissured expansive soils［J］.Shuili Xuebao，2004，（6）：108－113.（in Chinese））

［6］ 徐永福，孫婉瑩，吳正根.我國膨脹土的分形結(jié)構(gòu)的研究［J］.河海大學(xué)學(xué)報，1997，（1）：18－23.（XU Yongfu，SUNWan-ying，WU Zheng-gen.On fractal structure of expansive soils in China［J］.Journal of Hohai University，1997，（1）：18－23.（in Chinese）

［7］ 易順民，黎志恒，張延中.膨脹土裂隙結(jié)構(gòu)的分形特征及其意義［J］.巖土工程學(xué)報，1999，（5）：294－298.（YI Shun-min，LI Zhi-heng，ZHANG Yan-zhong.The fractal characteristics of fractures in expansion soil and its significance［J］.Chinese Journal of Geotechnical Engineering，1999，（5）：294－298.（in Chinese））

［8］ 謝和平.分形-巖石力學(xué)導(dǎo)論［M］.北京：科學(xué)出版社，1996.（XIE He-ping.Fractal-Introduction to Rock Mechanics［M］.Beijing：Science Press，1996.（in Chinese）

（編輯：王 慰）

Field Observation and Description M ethod of Cracks Development on Expansive Rock

WANG Jun1，GONG Bi-wei2，ZHANG Jia-jun1，ZHOU Xiao-wen1，TAN Feng-yi2
（1.State Key Laboratory of Subtropical Building Science，South China University of Technology，Guangzhou 510640，China；2.Key Laboratory of Geotechnical Mechanics＆Engineering of MWR，Yangtze River Scientific Research Institute，Wuhan 430010，China）

On the basis of numerical processing technology of image grayscale and binarization，the concepts of crack rate and calculation method on expansive rock are presented.The dry-wet cycle testswere carried out in situ at Luwangfen test section of South-to-North Water Transfer Project Central Route，the tests focused on the observations of development changes in cracks and soilmoisture changes.According to the results in the tests，the variation law of crack width on expansive rock slope and its relationship with moisture content during the dry-wet cycle testwere quantitatively analyzed，also the variation law of crack ratio was investigated.The development of crack on expansion rock is closely related to changes in theweather，and with the developmentofwet and dry cycles，the width of cracks and crack ratio will gradually increase.By using the Mat lab digital image processingmethod，the crack ratio can be used to effectively describe the growth and development of cracks on expansive rock.

expansive rock（soil）；South-to-North Water Transfer Project Central Route；Mat lab；crack ratio

TU455

A

1001－5485（2010）09－0074－05

2009-11-24

水利部巖土重點實驗室開放基金項目（G07－02；G07－15）

王 軍（1983-），男，湖北荊州人，碩士研究生，主要從事特殊土方面的研究，（電話）13476807917（電子信箱）wangjun_ust＠qq.com。