基于CT掃描的斜坡非飽和帶土體大孔隙定量化研究和三維重建

徐宗恒， 徐則民， 李凌旭

(1.昆明理工大學(xué) 建筑工程學(xué)院， 云南 昆明 650500; 2.云南師范大學(xué) 旅游與

地理科學(xué)學(xué)院， 云南 昆明 650500； 3 云南大學(xué) 滇池學(xué)院， 云南 昆明 650228)

基于CT掃描的斜坡非飽和帶土體大孔隙定量化研究和三維重建

徐宗恒1,2， 徐則民1， 李凌旭3

(1.昆明理工大學(xué) 建筑工程學(xué)院， 云南 昆明 650500; 2.云南師范大學(xué) 旅游與

地理科學(xué)學(xué)院， 云南 昆明 650500； 3 云南大學(xué) 滇池學(xué)院， 云南 昆明 650228)

摘要：[目的] 系統(tǒng)開展植被發(fā)育斜坡非飽和帶土體中大孔隙定量化和三維重建技術(shù)研究，為水分沿大孔隙迅速運(yùn)移研究提供可靠理論依據(jù)，對于降雨型滑坡研究具有基礎(chǔ)性意義。[方法] 采用CT技術(shù)對云南省昭通市盤河鄉(xiāng)頭寨滑坡滑源區(qū)左側(cè)斜坡區(qū)兩個尺寸為25 cm×25 cm×50 cm的大尺度土柱進(jìn)行掃描，得到一系列平面和縱向CT切片圖像。圖像處理并計算得到了土柱中大孔隙隨深度的變化情況以及3種物質(zhì)(礫石、土體基質(zhì)和大孔隙)的CT值范圍；借助VolView 3.4體可視化交互系統(tǒng)實(shí)現(xiàn)了三維大孔隙通道系統(tǒng)的重組，研究其在空間內(nèi)的連通性、分支性和復(fù)雜性。[結(jié)果] 包括根系通道、動物通道、干縮裂縫及團(tuán)聚體間的結(jié)構(gòu)性孔隙廣泛發(fā)育于非飽和帶土體中，且大孔隙分布異質(zhì)性明顯，波動變異較大，隨著深度的增加大孔隙呈逐漸減小的趨勢。三維重組清晰可見土柱內(nèi)含有較多獨(dú)立分布的大孔隙通道，它們具有較好的連通性。[結(jié)論] CT掃描技術(shù)作為非破壞性獲得技術(shù)在斜坡非飽和土體大孔隙定量化和三維重建研究中應(yīng)用是可行的。

關(guān)鍵詞：植被發(fā)育斜坡； CT掃描； 土體大孔隙； 定量化； 三維重建

植被發(fā)育斜坡區(qū)非飽和帶土體中，由植被根系腐爛后形成的通道、根系—土體和礫石—土體接觸帶空隙、動物通道、團(tuán)聚體間結(jié)構(gòu)性大孔隙、土體干縮以及氣候凍融交替產(chǎn)生的裂隙組成的大孔隙系統(tǒng)是水分滲透的重要通道。在降雨特別是暴雨條件下，這些大孔隙通道使得斜坡能同步吸收來自降雨過程的全部或大部分水量而基本不產(chǎn)生地表徑流[1]，導(dǎo)致斜坡臨時性飽和條件的發(fā)育，地下水位的持續(xù)變動以及加快加強(qiáng)了斜坡的水—巖相互作用，這對于植被發(fā)育斜坡的失穩(wěn)起到了巨大的推動作用。所以水分沿斜坡土體大孔隙通道迅速下滲是降雨型滑坡孕育的一個基礎(chǔ)性環(huán)節(jié)，而開展大孔隙通道定量化研究和三維重建能為水分沿其快速下滲研究提供可靠理論依據(jù)，這對于降雨型植被發(fā)育斜坡失穩(wěn)的研究、揭示發(fā)生機(jī)理具有重要的意義。

土體大孔隙的研究方法有多種，其中CT掃描作為一種無損檢測技術(shù)，非破壞性、精確性、全面性、快速性等優(yōu)點(diǎn)為研究突破了許多瓶頸問題。例如Warner等[2]用CT技術(shù)掃描確定了發(fā)育有肯塔基藍(lán)草(Poapratensis)，20.3 cm直徑的圓柱形土柱內(nèi)大孔隙的數(shù)目、大小等，實(shí)現(xiàn)了大孔隙的可視化，并用染色示蹤試驗揭示了大孔隙在縱剖面上的連續(xù)性。Anderson等[3]用同樣方法對取自耕田和森林土體截面面積為76.4 mm×76.4 mm的樣品進(jìn)行掃描，得到大孔隙分布情況并與含有已知尺寸人造大孔隙掃描結(jié)果進(jìn)行對比，結(jié)果表明，該方法對土壤結(jié)構(gòu)的評價有很好的運(yùn)用；馮杰等[4]對取自于南京郊區(qū)農(nóng)場中內(nèi)徑為13.1 cm的原狀土柱和含有已知人造直徑大孔隙的填充土柱進(jìn)行了CT掃描實(shí)驗，得到大孔隙在橫截面上的大小、數(shù)目、形狀以及分析了縱斷面上的連通性。吳華山等[5]采集了太湖地區(qū)3種主要水稻土(白土、黃泥土和烏柵土)，利用該方法得出了不同深度土壤剖面上各種大孔隙的數(shù)量、面積以及分布狀況，所供試的土樣直徑為12.8 cm。以上有關(guān)研究所針對對象多集中于農(nóng)業(yè)土壤，植被發(fā)育斜坡土體涉及較少，然而，農(nóng)業(yè)土壤中大孔隙與植被發(fā)育斜坡土體中大孔隙的形成和發(fā)育背景有著明顯的差異，因為農(nóng)業(yè)土壤大孔隙受耕作模式控制，大孔隙發(fā)育受限，再者，農(nóng)業(yè)土壤中植物根系發(fā)育規(guī)模和側(cè)向展布程度較斜坡土體要低，與根系相關(guān)形成的大孔隙數(shù)量也會有所降低，所以針對植被發(fā)育斜坡土體大孔隙進(jìn)行研究具有較好的針對性。另外，既有研究多集中于依據(jù)土柱平面CT切片來分析土體中大孔隙的分布和不同位置處的差異；研究對象尺度較小，缺乏代表性，實(shí)際上，依據(jù)平面CT切片僅能反映出土柱內(nèi)大孔隙的平面分布，而不能反映大孔隙通道在空間的連通性、分支性和復(fù)雜性，所以需要對大孔隙通道進(jìn)行三維可視化研究來實(shí)現(xiàn)，這也是大孔隙定量化研究的一個重要方面。本文對兩個取自天然斜坡的大尺度(25 cm×25 cm×50 cm)方形原狀土柱進(jìn)行掃描，對得到的一系列水平和縱剖面CT切片進(jìn)行圖像處理，分析大孔隙隨深度變化情況以及縱剖面上的連通性，然后將數(shù)據(jù)導(dǎo)入VolView 3.4體可視化交互系統(tǒng)對CT切片圖像進(jìn)行三維重組研究大孔隙在空間的分布，以期為植被發(fā)育斜坡土體滲流等物理特性研究提供理論依據(jù)。

1取樣地點(diǎn)和樣品

取樣點(diǎn)位于云南省昭通市昭陽區(qū)盤河鄉(xiāng)頭寨滑坡滑源區(qū)左側(cè)斜坡區(qū)，該地曾于1991年9月23日發(fā)生過造成216人死亡的大型高速遠(yuǎn)程滑坡，是當(dāng)時20世紀(jì)90年代中國發(fā)生的較為慘重的重大滑坡災(zāi)害事件[6-7]。取樣區(qū)內(nèi)發(fā)育有長勢較好的天然森林植被，通過對林區(qū)進(jìn)行植被調(diào)查，試驗區(qū)植被覆蓋度95%以上，主要植被類型共計31種，有滇青岡(Cyclobalanopsisglaucoides)、云南柳(Salixcavaleriei)、玉山竹屬(Yushaniasp.)、滇白楊(Populusrotundifoliavar.bonati)、牛奶子(Elaeagnusumbellata)、綿毛杜鵑(Rhododendronfloccigerum)等。取樣點(diǎn)地理位置為東經(jīng)27°34′2″，北緯103°51′34″。

既有土柱的取樣方法主要為采用一端帶刃口的PVC管在重力擊打(或液壓系統(tǒng))作用下敲入土壤中，采取原狀樣[4-5,8]，但是該方法在植被發(fā)育斜坡區(qū)并不適用，因為土體中存在的縱橫交錯根系以及礫石很容易阻礙取樣器進(jìn)入土體，多次敲擊會造成土體的擾動。本文取樣的方法為：首先去除土體表層范圍內(nèi)的雜物和土體表層疏松腐殖層，然后從外圍100 cm×100 cm邊界逐漸往選定的25 cm×25 cm的核心區(qū)域緩慢開挖，當(dāng)土柱成型以后用毛刷將剖面上附著土顆粒及雜物清理干凈，使四周平面平整，而后將現(xiàn)場融化好的石蠟澆筑于土柱表面，再將多層保鮮膜、棉布包裹于土柱外表面；最后用特制的木框套取土體并將其送至云南省第三人民醫(yī)院進(jìn)行CT掃描。為了使得研究結(jié)果具有代表性，在距離大約為10 m的區(qū)域內(nèi)(定為TZ1和TZ2)共取得2個樣品(尺寸為25 cm×25 cm×50 cm)。

2CT掃描及原理

云南省第三人民醫(yī)院用于土樣掃描的CT儀為SIMENS Sensation Open型，進(jìn)行土樣掃描時設(shè)置掃描電流為176 mA，電壓為120 kV，每個體素的尺寸為1.0 mm×1.0 mm×1.5 mm，總掃描得到的體素為250×250×334(voxels)，即每一CT切片的層間距為1.5 mm。設(shè)置完參數(shù)以后，將土樣頂部朝向CT掃描架的方向推進(jìn)X—射線管和探測器之間，啟動掃描，放射源發(fā)出的X—射線穿透土樣以后，會引起射線強(qiáng)度、速度和頻率的改變，不同物質(zhì)對X—射線的吸收不同[9-10]，導(dǎo)致投影出來的亮度也不同。在CT切片中可以根據(jù)亮度(CT值)將各物質(zhì)區(qū)分開來，物質(zhì)的CT值越高，表征物質(zhì)密度越高，特殊的，水的CT值為0 HU，空氣的CT值為-1 000 HU，該次掃描土樣的CT值范圍為-1 024~3 071 HU。X射線數(shù)據(jù)的變化通過探測器檢測得到并形成投影，通過旋轉(zhuǎn)射線源，就可以得到一系列不同深度處的土樣CT切片，本次掃描每個土柱共得到334個切片圖像。將得到的切片圖像數(shù)據(jù)導(dǎo)入三維空間內(nèi)，在空間里將每一個坐標(biāo)點(diǎn)處用CT掃描所記錄到的數(shù)據(jù)進(jìn)行填充，賦予顏色、透明度等信息就可以將土體三維重組而成。

3圖像處理和計算結(jié)果分析

3.1　圖像處理

土體中大孔隙定量化研究關(guān)鍵環(huán)節(jié)之一為圖像處理中將大孔隙與土體基質(zhì)相分離，即將圖像處理為二值化圖像，用黑色(或白色)代表大孔隙，白色(或黑色)代表土體基質(zhì)，然后進(jìn)行定量化計算。馮杰等[4]采用的方法為基于CT切片對平面圖像假定一閾值，將大孔隙和基質(zhì)分離開來，與已知大孔隙大小進(jìn)行比較，如果相差太大，再重新假定一閾值進(jìn)行計算，直到相差值與已知值的比值很小為止；張季如和陶高梁等[11-12]采用的方法是對于同一圖像，經(jīng)多人處理后求平均值，其誤差控制在10%以內(nèi)即可。筆者認(rèn)為這些方法都是很好的大孔隙分離方法，但精確度有待提高，因為閾值的選擇是土體內(nèi)大孔隙分布評價的關(guān)鍵環(huán)節(jié)之一，如果閾值的選擇不同會導(dǎo)致出現(xiàn)不同的處理結(jié)果，有時會相差甚多。為了消除選擇閾值帶來的麻煩，用物質(zhì)的CT值來分離大孔隙可達(dá)到很好的分離效果，因為物質(zhì)的CT值代表著物質(zhì)的密度，不同物質(zhì)的CT值不同。為了得到土體內(nèi)各物質(zhì)的CT值范圍，本文對根系、土體基質(zhì)以及礫石單獨(dú)進(jìn)行CT掃描，并對每個土柱的334個切片進(jìn)行了詳細(xì)分析，得到該次土柱掃描CT值范圍為-1 024~3 071 HU，其中土體基質(zhì)CT值范圍為90~1 416 HU，礫石CT值范圍為1 417~3 071 HU，根系和土體中大孔隙CT值范圍為-1 024~89 HU，需特別說明的是，所取土柱內(nèi)含有大量的根系，根系的存在對大孔隙的產(chǎn)生和存在起到積極的作用，且根系通道也是土體內(nèi)大孔隙通道的一部分[13-15]，所以為了研究方便和提高精確度，本文未將根系獨(dú)立開來。

圖1所示為TZ2土柱某深度處的CT切片，定義土柱頂面Z=0，雖然在采集及運(yùn)送樣品過程中控制了對土體的擾動，但還是會對樣品邊角處造成少許損傷，所以去除土柱外圍一定區(qū)域，選定20 cm×20 cm×40 cm(深度)的立方體土柱為研究對象。

圖1　土柱CT切片示例

為了得到大孔隙隨著深度增加的變化情況，將得到的CT切片導(dǎo)入Adobe Photoshop CS4中進(jìn)行裁剪，得到選定的20 cm×20 cm為研究區(qū)域(圖1)，對裁剪以后的圖像此處只需簡單選擇一個閾值就能將圖片處理為二值圖像，使切片中土體基質(zhì)和礫石表現(xiàn)為白色，大孔隙顯示為黑色。進(jìn)行上述步驟處理以后，圖像中可能會出現(xiàn)噪點(diǎn)，即與實(shí)際情況不一致的獨(dú)立的像素點(diǎn)，這些像素點(diǎn)的存在會影響解析的精度，所以使用Image-Pro Plus 6.0軟件圖像處理中的形態(tài)學(xué)開運(yùn)算濾鏡功能，進(jìn)行降噪處理，除去獨(dú)立的噪點(diǎn)，再采用膨脹運(yùn)算(2×2范圍，重復(fù)1次)，將大孔隙邊界復(fù)原[16]，將得到的二值化圖片用matlab計算，可得到CT切片的大孔隙度(定義為圖像中大孔隙面積占整個CT切片圖像面積比例[11])，同時每隔7.5 mm深度計算CT切片上各物質(zhì)CT值的平均值、標(biāo)準(zhǔn)差、最大值和最小值。除了對土柱橫斷面進(jìn)行掃描計算外，還對土柱的2個縱向進(jìn)行了掃描，依此來揭示縱剖面上大孔隙的連通性、復(fù)雜性和分支性[4]。

3.2　計算結(jié)果分析

對于大孔隙的定義，相關(guān)研究并未達(dá)成一個共識，并且一直是一個有爭議問題，Warner等[2]利用CT掃描方法觀察取自肯塔基藍(lán)草地(Kentucky bluegrass)，直徑為20.3 cm的土柱中大孔隙時，所研究對象為孔徑大于1 mm的大孔隙；余長洪等[17]用浸漬劑固化再切片方法，研究了黏土和砂壤土中大于85 m的大孔隙；吳華山等[5]用CT掃描研究的是白土、黃泥土和烏柵土3種土樣中大于0.3 mm孔徑的大孔隙，而本文由于CT掃描儀精度的限制，通過計算，能識別的最小孔隙為孔徑是0.44 mm的大孔隙。由圖1可以看出，平面零星的分布著幾塊礫石，偶見體積較大的礫石；根系在空間內(nèi)穿叉而過的跡象明顯，在水平切片上成圓率較高、獨(dú)立顯色面積寬或者縱截面上呈條狀分布的位置可能為根系分布，水平截面上分布有多條寬而窄，形狀各異，成圓率并不高的大孔隙通道，這些通道在各水平切片上均會發(fā)生變化，并且分布具有聚集性的特點(diǎn)，這些大孔隙在空間上相互連通，形成連通性較好的大孔隙通道。從CT切片的縱截面上孔隙分布的特點(diǎn)也可以看出，大孔隙(特別是根系形成的大孔隙)在深度方向上并不是垂直向下的，而是彎延向下，在某深度處出現(xiàn)又在此深度處消失，隨之轉(zhuǎn)至另一個縱剖面位置處；宏觀上，表層20 cm范圍內(nèi)未見單體體積較大的礫石，隨著深度的增加，礫石的數(shù)量和體積均有增加和增大的趨勢。

CT值代表著物質(zhì)的密度，其變化可以一定程度上反映土柱內(nèi)物質(zhì)組成(包括大孔隙)的變化。在所得的橫截面CT切片中，每隔7.5 mm提取一張圖像進(jìn)行計算，計算得到的的各參數(shù)結(jié)果隨深度的變化情況如圖2所示，圖2中“AV”表示平均CT值，“STD”表示標(biāo)準(zhǔn)差，“MAX”表示最大CT值，“MIN”表示最小CT值，由圖2可以看出，最小CT值沿深度波動較小，介于-1 000~-1 024 HU之間，最大值則波動較大，TZ1土柱為2 153~3 071 HU，TZ2土柱為2 189~3 071 HU，這是由于頭寨試驗區(qū)斜坡非飽和

帶長期處于潮濕環(huán)境條件中，在不同季節(jié)氣態(tài)水分子循環(huán)入侵礫石結(jié)構(gòu)體，導(dǎo)致其內(nèi)部造巖礦物溶解，形成一定厚度的腐巖殼，有的甚至形成了腐巖結(jié)構(gòu)體，腐巖的密度較新鮮巖石要低，它的存在致使CT最大值產(chǎn)生了一定的波動；另外，平均CT值的變化可以說明土柱內(nèi)物質(zhì)隨深度的變異程度，同時在一定程度上可以反映大孔隙隨深度的變化趨勢。為了表現(xiàn)出平均CT值與大孔隙度的變化關(guān)系，將兩者單獨(dú)列于圖3中，從圖3中可以明顯看出，兩者隨著深度的增加均表現(xiàn)出來很強(qiáng)的波動不穩(wěn)定性，TZ2的平均CT值在0~200 mm內(nèi)變化波動很小，而在200 mm以下跳躍著升高，在深度375 mm處達(dá)到最高值396.7 HU；TZ1平均CT值多集中在70~180 HU之間。對于大孔隙度，TZ1分別在深度37.5 mm和底部位置為最大值(37.95%)和最小值(7.04%)，TZ2在深度7.5 mm和底面處取得最大值(32.06%)和最小值(11.37%)，兩個點(diǎn)處大孔隙度均隨著深度的增加呈鋸齒狀整體減小趨勢，主要是由于隨著深度增加，制約著大孔隙形成與發(fā)展的因素(植物根系密度、有機(jī)質(zhì)含量、動物群體數(shù)量和個體大小)也隨之減少?？傮w上，對于同一土柱，隨深度變化，大孔隙度和平均CT值的變化趨勢正好相反，這是由于在不同的深度處，若大孔隙的數(shù)量較多，則-1 024~89 HU之間的CT值數(shù)量就較多，導(dǎo)致平均CT值降低，但也有特殊的情況存在，例如在TZ1深度為190~210 mm處，兩者均呈遞減的變化趨勢，這可能是由于該位置處大孔隙的數(shù)量與CT值較大的物質(zhì)(如礫石)的含量呈正相關(guān)關(guān)系。圖2中所示的“AV+STD”和“AV-STD”表示各個剖面間平均CT值隨著深度增加而表現(xiàn)出來的波動變異程度。由圖2—3可以看出，大孔隙分布在同一深度不同土柱間以及不同深度處同一土柱內(nèi)均表現(xiàn)出明顯的異質(zhì)性，且波動變異均較大。

圖2　土柱物質(zhì)CT值隨深度的變化

圖3　土柱平均CT值和大孔隙度隨深度的變化

4大孔隙結(jié)構(gòu)的三維重組

大孔隙結(jié)構(gòu)的三維重組借助VolView 3.4軟件來實(shí)現(xiàn)，VolView是美國Kitware公司開發(fā)的，是一款開放源代碼、直觀的體可視化交互系統(tǒng)，軟件目的在于讓研究者在Windows, Mac或Linux下迅速探究和分析復(fù)雜三維科學(xué)數(shù)據(jù)[18]，VolView是基于Kitware公司開發(fā)的兩個非常著名的醫(yī)學(xué)影像處理開發(fā)包VTK和ITK上開發(fā)出來的，并且提供免費(fèi)的使用下載[19]，VolView軟件平臺提供了DICOM格式數(shù)據(jù)訪問和處理窗口，這些數(shù)據(jù)多來自CT掃描數(shù)據(jù)，本文正是依此來實(shí)現(xiàn)對三維掃描數(shù)據(jù)進(jìn)行導(dǎo)入、分析和處理，其原理是利用3D工具里的體繪制(volume rendering)功能。

導(dǎo)入CT掃描數(shù)據(jù)以后的三維重組過程其實(shí)就是二維CT數(shù)據(jù)按照位置和角度信息進(jìn)行規(guī)則化處理，在三維空間里將每一個坐標(biāo)點(diǎn)處用CT掃描所記錄到的數(shù)據(jù)進(jìn)行填充，然后賦予顏色，透明度和CT值的過程。VolView的好處在于它提供了很好的人機(jī)交互畫面來進(jìn)行參數(shù)調(diào)節(jié)[20]，可以將其中CT值對應(yīng)的物質(zhì)賦予不同的透明度(取值0~1)，以及不同的顏色。為了研究2個土柱內(nèi)大孔隙的三維展布情況，只需將CT值為-1 024~89之間物質(zhì)透明度設(shè)置為1即可。本文中由于土柱內(nèi)局部多個剖面大孔隙度高達(dá)30%以上，將全部大孔隙顯示出來將不利于大孔隙觀察，所以為了便于觀察和研究，此處僅將CT值為-1 024~-684的物質(zhì)顯示出來。圖4(a—b)和4(c—d)為TZ1和TZ2從不同方位視角觀察部分大孔隙的示意圖，圖4中方位標(biāo)示塊中“A”和 “P”分別表示斜坡下坡、上坡方向的縱剖面；“R”和 “L”分別表示土柱右側(cè)、左側(cè)縱剖面；“S”表示的為土柱的頂部方向。

TZ1土柱　　　　　　　　　　　　　　　　　　　　　　TZ2土柱

圖4不同方位視角下土柱內(nèi)大孔隙三維分布

圖4很好地表現(xiàn)了土柱內(nèi)部分大孔隙通道在土柱內(nèi)的三維展布情況，從圖4(a—b)中可以看出，TZ1土柱內(nèi)部含有較多的錯綜復(fù)雜，有的彎曲往下，有的近水平分布的大孔隙通道，其直徑以及數(shù)量從土柱頂部至底部慢慢變大增加而后變小減少。相比而言，圖4(c—d)中TZ2內(nèi)部大孔隙則較多分布在250 mm深度之上，偶見直徑較大的通道，數(shù)量總體上沒有TZ1土柱多。TZ1和TZ2土柱內(nèi)部均含有較多獨(dú)立的大孔隙，這些大孔隙的成圓率不高，在有些位置能形成連通性好的大孔隙，但是這些孔隙并不是垂直向下，而是彎延向下或水平延伸。從宏觀上觀察，TZ1土柱中大孔隙在較多深度處均比TZ2的多，這與圖3中的結(jié)果一致。CT切片圖像分析與三維觀察得到的結(jié)果相一致驗證了該方法在土體大孔隙定量分析和三維重組研究中的可行性。

5結(jié) 論

(1) CT掃描技術(shù)作為非破壞性獲得技術(shù)在土壤科學(xué)研究中得到廣泛運(yùn)用，其精確性、全面性、快速性等優(yōu)點(diǎn)是顯而易見的，本文成功將此技術(shù)運(yùn)用到大尺度土柱內(nèi)部大孔隙定量化研究和三維重建中，表明CT掃描技術(shù)在斜坡非飽和帶土體大孔隙的研究中是可行的。

(2) 對取自于云南省昭通市頭寨滑坡滑源區(qū)2個大尺度土柱進(jìn)行CT掃描，對獲得的平面CT切片進(jìn)行分析處理和分析，得到土柱內(nèi)土體基質(zhì)CT值為90~1 416 HU，礫石為1 417 ~3 071 HU，根系和土體中大孔隙為-1 024~89 HU；基于CT值范圍用圖像處理方法對土體大孔隙進(jìn)行分離，得到土體內(nèi)大孔隙度、平均CT值、最大CT值、最小CT值以及標(biāo)準(zhǔn)差隨深度增加而變化的情況，結(jié)果表明，大孔隙分布異質(zhì)性明顯，波動變異較大，隨著深度的增加大孔隙度呈逐漸降低的趨勢，TZ1土柱在深度37.5 mm和底部位置處取得最大值(37.95%)和最小值(7.04%)，TZ2樣在深度7.5 mm和底面處取得最大值(32.06%)和最小值(11.37%)。

(3) 借助VolView 3.4 體可視化交互系統(tǒng)實(shí)現(xiàn)了2個土柱內(nèi)部部分大孔隙的三維可視化，展示了土柱內(nèi)部分大孔隙通道空間內(nèi)的連通性、分支性和復(fù)雜性，三維重組清晰可見土柱內(nèi)含有較多獨(dú)立分布的大孔隙通道，它們具有較好的連通性，絕大部分呈彎延向下或水平延伸的狀態(tài)。三維重組定性的表現(xiàn)結(jié)果與二維切片計算結(jié)果一致，隨著深度的增加，兩個樣品土柱內(nèi)大孔隙呈降低趨勢，且在較多深度處，TZ1土樣的大孔隙分布均比TZ2的高。研究結(jié)果能為植被發(fā)育斜坡區(qū)土體大孔隙滲流和物理力學(xué)特性研究提供理論依據(jù)。

[參考文獻(xiàn)]

[1]徐則民,黃潤秋.山區(qū)流域高蓋度斜坡對極端降雨事件的地下水響應(yīng)[J].地球科學(xué)展,2011,26(6):598-607.

[2]Warner G S, Nieber J L, Moore I D, et al. Characterizing macropores in soil by computed tomography[J]. Soil Science Society of America Journal, 1989,53(3):653-660.

[3]Anderson S H, Peyton R L, Gantzer C J. Evaluation of constructed and natural soil macropores using X-ray computed tomography[J]. Geoderma, 1990,46(1):13-29.

[4]馮杰,郝振純.CT掃描確定土壤大孔隙[J].水科學(xué)進(jìn)展,2002,13(5):611-617.

[5]吳華山,陳效民,陳粲.利用CT掃描技術(shù)對太湖地區(qū)主要水稻土中大孔隙的研究[J].水土保持學(xué)報,2007,21(2):175-178.

[6]Zhu jiawen, Xu weiya, Yang xingguo, et al. The 28 October 1996 land slide and analysis of the stability of thecurrent Huashiban slope at the Liangjiaren Hydropower Station, Southwest China[J].Engineering Geology, 2010,114(1):45-56.

[7]徐則民,黃潤秋,唐正光.頭寨滑坡的工程地質(zhì)特征及其發(fā)生機(jī)制[J].地質(zhì)論評,2007,53(5):691-698.

[8]Grevers M C J, Jong E D E, Arnaud R J S T. The characterization of soil macroporosity with CT scanning[J]. Canadian Journal of Soil Science, 1989,69(3):629-637.

[9]Peth S. Applications of microtomography in soils and sediments[J]. Developments in Soil Science, 2010,34:73-101.

[10]王路君，左永振，孔憲勇，等.CT技術(shù)在巖土工程研究中的應(yīng)用[J].地下空間與工程學(xué)報,2009,5(S2):1754-1775.

[11]張季如,祝杰,黃麗，等.固結(jié)條件下軟黏土微觀孔隙結(jié)構(gòu)的演化及其分形描述[J].水力學(xué)報,2008,39(4):394-400.

[12]陶高梁,張季如.表征孔隙及顆粒體積與尺度分布的兩類巖土體分形模型[J].科學(xué)通報,2009,54(6):838-846.

[13]Heeraman D A, Hopmans J W, Clausnitzer V. Three dimensional imaging of plant roots in situ with X-ray computed tomography[J]. Plant and Soil, 1997,189(2):167-179.

[14]Gantzer C J, Anderson S H. Computed tomographic measurement of macroporosity in chisel-disk and no-tillage seedbeds[J]. Soil and Tillage Research, 2002,64(1):101-111.

[15]Luo Lifang, Lin Henry, Li Shangcai. Quantification of 3-D soil macropore networks in different soil types and land uses using computed tomography[J]. Journal of Hydrology, 2010,393(1):53-64.

[16]王偉. 三峽庫區(qū)紫色砂巖林地土壤優(yōu)先流特征及其形成機(jī)理[D]. 北京: 北京林業(yè)大學(xué), 2011.

[17]余長洪,周明耀.兩種土壤大孔隙結(jié)構(gòu)差異的數(shù)字圖像分析[J].廣東水利水電,2007(1):27-35.

[18]Kitware homepage.volview 3.4 now available[EB/OL].[2015-03-05]http:∥www.kitware.com/opensource/volview.html.

[19]李恩中.醫(yī)學(xué)圖像處理與分析軟件平臺綜述[J].計算機(jī)科學(xué)與探索,2008,2(5):467-477.

[20]黃易.三維數(shù)據(jù)場可視化關(guān)鍵技術(shù)研究[D].哈爾濱:哈爾濱工程大學(xué),2006.

Soil Macropores Quantification Study and 3D Reconstruction in Vadose Zones of Hillslope Based on X-Ray Computed Tomography

XU Zongheng1,2， XU Zemin1， LI Lingxu3

(1.DepartmentofCivilEngineering,KunmingUniversityofScienceandTechnology,Kunming,

Yunnan650500,China; 2.SchoolofTourismandGeographicScience,YunnanNormalUniversity,

Kunming,Yunnan650500,China; 3.DianchiCollege,YunnanUniversity,Kunming,Yunnan650228,China)

Abstract：[Objective] Study on soil macropores quantification and 3D reconstruction is particularly important as it can provide reliable theoretical data for preferential flow study. It is also of great importance in precipitation induced landslide studies.[Methods] We took two 25 cm×25 cm×50 cm soil cores from vegetated slope in Panhe Town, Yunnan Province. Base on X-ray computed tomography, a series of plane and longitudinal slice images were collected. By images analysis, variations of soil macropores with the increase of soil depth and CT value of three substances(gravel, matrix, macropore) were studied. Soil macropore 3D reconstruction was achieved by visual interactive system VolView 3.4 to investigate the connectivity, branchedness and complexity of macropore. [Results] The macropores consisting of root channels, faunal tunnels, shrinkage cracks and inter-aggregate porosity were broadly found in vadose zones of vegetated slope soils. Heterogeneity in distribution of soil macropores was significant, and macroporosity decreased with the increasing soil depth. As clearly demonstrated by 3D reconstruction, macropores were well distributed and connected in the soil cores. [Conclusion] We highlight the importance of CT technique in soil macropore quantification and 3D reconstruction.

Keywords:vegetated slope; X-ray computed tomography; soil macropores; quantification; 3D reconstruction

文獻(xiàn)標(biāo)識碼：B

文章編號：1000-288X(2015)01-0133-06

中圖分類號：P642.22

收稿日期：2014-01-06修回日期：2014-01-24

資助項目：國家自然科學(xué)基金云南聯(lián)合基金重點(diǎn)項目“植被發(fā)育斜坡對強(qiáng)降雨過程的地下水響應(yīng)機(jī)制及其災(zāi)害效應(yīng)”(U1033601)； 高等學(xué)校博士學(xué)科點(diǎn)專項科研基金項目“峨眉山玄武巖巖崩災(zāi)害孕育機(jī)理”(20135314110005)

第一作者：徐宗恒(1987—)，男(漢族)，云南省麗江市人，博士，講師，主要從事斜坡地質(zhì)災(zāi)害研究。E-mail：xuzh208@163.com。

通信簡介：徐則民(1963—)，男(漢族)，河北省承德市人，博士，教授，博士生導(dǎo)師，主要從事斜坡穩(wěn)定性研究。E-mail：abc5100@188.com。