粘土顆粒等效孔徑的算法研究

耿大新，王迎迎，李宇晗

(1.華東交通大學 土木建筑學院巖土工程研究所，江西 南昌 330013；2.南京水利科學研究院 巖土工程研究所，江蘇 南京210029)

?

粘土顆粒等效孔徑的算法研究

耿大新1，王迎迎1，李宇晗2

(1.華東交通大學 土木建筑學院巖土工程研究所，江西 南昌 330013；2.南京水利科學研究院 巖土工程研究所，江蘇 南京210029)

非飽和土中隨機復雜的孔隙特征，對土體的工程性質產生重要的影響。為研究土體微觀特性，以不同黏粒含量的江西重塑紅粘土為研究對象，進行土-水特征實驗和土樣一維毛細上升實驗。研究得到土樣孔隙分布情況，提出粘土顆粒等效孔徑的算法公式，給出增量步取值方式并詳細介紹了算法過程；驗證了等效孔徑公式的合理性，實現了微觀結構的定量化；同時為最大毛細上升高度的計算提供了一種有效的方法。

非飽和土；等效孔徑；毛細上升

針對非飽和土的土體結構，其粒徑方面的研究較多，但土顆粒間的孔隙大小、形狀及分布亦是土體結構組成的重要因素，對土體的毛細作用等[1-3]工程行為有著重要的影響。Arya等[4]建立經驗模型，將粒度分布轉化為孔徑分布進行分析，通過土體的顆粒級配、容重等參數預測其水分特征。Griffiths等[5]對不同種類的粘土進行分析，研究其固結過程中的孔隙尺寸分布情況。張先偉等[6]通過室內壓縮試驗，研究結構性粘土的孔隙變化規(guī)律，以深入了解土的變形機制。梁建偉等[7]經試驗指出，對于極細粘土，孔隙微尺度效應是影響滲流特性的重要因素，而土顆粒比表面積、孔隙尺度等是重要影響參數。陶高梁等[8]采用三種不同的試驗方法，通過壓縮變形影響下土體結構內部孔隙孔徑分布規(guī)律，研究土-水特征曲線的變化規(guī)律，指出在不同的壓縮變形條件下，土體累積孔隙體積及土-水特征曲線變化規(guī)律相近，整體上呈“掃帚型”分布形式。孫德安等[9]運用三種方法研究原狀土和重塑土的土-水特征和孔隙分布，指出兩種土的孔隙結構依次為單峰和雙峰，且在不同干密度條件下，壓實樣內部粒間孔隙分布幾乎相同。宋盛淵等[10]通過壓泵試驗并結合分形理論劃分吹填土的孔徑分布區(qū)間，研究在不同的地基處理方式下，吹填土的孔隙分布特征。

然而，已有的研究多圍繞土樣的粒度、孔隙結構及孔徑分布等方面展開，并沒有給出孔隙孔徑的定量算法。土體內部孔徑分布的隨機性、復雜性、不均勻性等特點給孔徑的計算增加了難度。本文首先通過理論推導，提出不同黏粒含量土樣等效孔徑的計算公式；再結合室內實驗結果，驗證其合理性與可行性。

1　算法研究

眾所周知，Young-Laplace方程是用于量化非零接觸角的氣-液-固三相交界面上壓降的典型分析模型，該方程可表達為：

(1)

式中：ua-uw為基質吸力；Ts為液相表面張力；α為土顆粒與孔隙水之間的接觸角；r1、r2為理想的土顆粒與氣-水交界處彎液面的兩個曲率半徑；對于直徑為d的理想的圓柱形毛細管，r1=r2=d/2，則式(1)可表達為：

(2)

土體孔隙結構決定了處于平衡狀態(tài)的蒸汽壓大小[11]。于楊-拉普拉斯方程的基礎之上，1871年，William Thomson在不違背土體體系原有特性的前提下，把土體孔隙理想化為一個直徑為d、接觸角為α、表面張力為Ts的毛細管，提出了用以描述氣-水交界曲面上的壓力變化與交界曲面上方蒸汽壓之間的聯系的開爾文公式，整理如下：

(3)

式中：R=8.314 J/(mol·K)為通用氣體常數；T=298 K為熱力學溫度；uv0為溫度T條件下平衡狀態(tài)的自由水飽和蒸汽壓值(kPa)；uv1為平衡狀態(tài)的溶液蒸氣壓值(kPa)；vw=18×10-6m3/mol為水蒸氣的偏摩爾體積。定義uv1/uv0為相對濕度RH，根據式(3)可以建立相對濕度與基質吸力之間的函數關系：

(4)

由以上推導可得，毛細直徑d用相對濕度表示為：

(5)

或表示為基質吸力形式：

(6)

(7)

則每單位質量固體內水的排出量為：

(8)

式中：wi為含水量，可借助土-水特征實驗建立基質吸力ua-uw與含水量之間的對應關系。假設固-液相互作用只對土顆粒表層附近水膜密度產生顯著影響，則毛細吸附區(qū)ρw本質上是常數，即ρw=1 g/cm3。

(9)

其中，相對濕度在第i步增量時所對應的吸附水膜厚度ti用Halsey方程表示如下：

(10)

式中：τ表示吸附水分子的有效直徑。

把計算得到的相鄰兩增量步之間的孔隙半徑取平均值，得平均孔隙半徑：

(11)

然而，以上推導是建立在毛細管假設的基礎之上，已有的研究并沒有給出土體孔隙孔徑的定量計算方法?？紤]到土體結構內部孔隙大小不一、分布不均等特性，以上述理論推導為依托，在總結了大量的實驗結果的基礎上，按照統計理論中加權平均值計算公式，給出粘土顆粒等效毛細直徑的計算公式:

(12)

2　實驗方案

為研究粘土土樣孔徑分布情況，現于江西省南昌市孔目湖地區(qū)取土制樣進行實驗。取土后過0.5 cm的細篩，原狀土定名為CL-61(黏粒含量60.8%)，土粒比重為2.71。摻入砂配置另外兩種土樣，根據黏粒含量定名為CL-50(黏粒含量50.2%)、CL-38(黏粒含量38.4%)，土粒比重分別為2.68、2.65。將試樣置于封閉的空間，實驗期間室內溫度控制在20～25 ℃。

2.1土-水特征實驗

在最優(yōu)含水率條件下對三種土樣標準擊實，采用DIK-4303土壤pF水分特征曲線測定儀進行土-水特征實驗。首先將制好的三種土樣放入真空飽和器中24 h，以保證充分飽和；然后基質吸力分別控制為25 kPa、50 kPa、75 kPa、100 kPa、125 kPa、160 kPa、200 kPa、250 kPa、300 kPa、350 kPa、400 kPa、500 kPa進行加壓，并且在每次加壓前稱取每個試樣的質量。

2.2土樣一維毛細實驗

使用高1.1 m、直徑11 cm的實驗用PVC管制樣，在最優(yōu)含水率條件下每放入300 g后進行標準擊實30次，裝樣高度為100 cm。每組土樣按上述方法制作5根，共15根分為3組進行編號。PVC管底部放兩層濾紙用紗網包裹，再用橡皮筋將其固定。

圖1　土樣一維毛細實驗Fig.1　One-dimensional capillary experiment of samples

將3組土樣靜置于水面高度恒為20cm的蓄水箱中(見圖1)，水箱底部鋪一層透水性良好的細砂以保證PVC管底面與水充分接觸。靜置30d后從下往上依次每隔5cm進行切割，并按照《公路土工試驗規(guī)程》(JTG E40-2007)規(guī)定，測定切割后樣塊的含水率和壓實度。

3　等效孔徑計算

為排除每種土樣人為壓實誤差的影響，將每組去掉一個最大值和一個最小值，得3種土樣在脫濕過程中的土-水特征曲線(見圖2)，CL-61(0.95)表示CL-61土樣壓實度為0.95。當孔徑范圍大約為10-9～10-4m，相應的基質吸力處于0～144 000 kPa，相對濕度為35%～100%時，毛細作用為土中孔隙持水的主要因素，適用于毛細孔徑計算方法。

為得到較大吸力范圍內的土-水特征曲線，使用MATLAB軟件Van Genuchten模型(簡稱VG模型)進行擬合[15]，擬合方程為：

(13)

式中：a、m、n均為相關系數；θ、θs、θr分別為含水率、飽和含水率和殘余含水率。

圖2　土-水特征曲線Fig.2　Soil-water characteristic curves

1996年Clayton指出，當m和n為相互獨立的變量時，所得擬合曲線能更好地與室內試驗數據相吻合[9]，因此本文中土-水特征曲線的擬合將a、m和n均視為相互獨立的參數。擬合參數見表1。

表1　SWCC參數擬合結果

以(ua-uw)i=(ua-uw)i-1+50(i-1)，(i=2，3，4，…)，(ua-uw)1=100 kPa的方式選取基質吸力在0～144 000 kPa之間的數據進行計算，以600 kPa、850 kPa的CL-50(0.95)為例，則體積含水率：

+0.014 0=10.99%

+0.014 0=5.24%

轉化為質量含水量：wi-1=0.041 0 (g/g)，wi=0.019 6 (g/g)。

通過式(4)得到對應的相對濕度：

由式(7)計算單位質量土體被氣體或水填充的孔隙體積，即：

則在第i步時，每單位質量固體內水所占的孔隙體積減小量為：

假設液體水分子所占區(qū)域的橫截面面積約為A=10.8 ?2，Avogadro常數NA=6.02×1023/mol，吸附水分子的有效直徑：

代入到式(10)得吸附水膜的厚度為：

再由式(9)計算實際孔隙半徑：

用(11)式計算平均孔隙半徑：

表2　等效毛細直徑

4　公式驗證

在非飽和土體結構內部，曲折連通的微小孔隙的存在使其具備了毛細上升的前提條件。在毛細管力的作用下，水分沿著相互連通的孔隙通道向上遷移。當毛細管內的水柱達到最大高度時，我們認為重力作用下的水柱質量與沿水-固交界面上的表面張力達到平衡。如圖3所示，通過分析虛線圓圈內毛細上升水柱的力學平衡，可以推算毛細上升的最大高度hmax，見式(14)。

(14)

式中：g表示重力加速度。

圖3　小直徑管內毛細上升的力學平衡Fig.3　Mechanical equilibrium of capillary rise in small diameter tubes

將求得的等效毛細直徑代入式(14)，計算最大毛細上升高度，并與實驗結果對比，如表3所示。由表3可知，計算值與實驗值基本符合，說明該種計算方式具有一定的合理性。

表3　計算與實驗結果對比

5　結　論

本文對重塑紅粘土進行研究，通過理論計算和大量室內實驗得到幾點結論。

1)根據實驗總結并借助于統計理論，推導獲得了粘土等效孔徑的經驗公式，增量步以非常數公差遞增數列的方式取值，詳細介紹了計算過程及公式的適用條件。

2)根據公式法求得的毛細水最大上升高度和一維毛細實驗結果較為接近，基本相符，證明了不同黏粒含量土樣等效孔徑計算公式是合理可行的?？蔀轭愃蒲芯刻峁┫鄳獏⒖迹瑫r也可作為估算路基最大毛細上升高度的一種方法。

3)本文獲得的經驗公式是基于較多假設和簡化而建立的，是對非飽和土孔隙定量研究的初步探討，因此在后續(xù)研究中有必要進行完善和全面的分析。

[1]LI J,ZHAO C,CAI G,et al.A model considering solid-fluid interactions stemming from capillarity and adsorption mechanisms in unsaturated expansive clays[J].Chinese Science Bulletin,2014,59(26):3314-3324.

[2]劉杰,姚海林,盧正，等.非飽和土路基毛細作用的數值與解析方法研究[J].巖土力學,2013,34(2):421-427.

LIU Jie,YAO Hailin,LU Zheng，et al.Study of analytic and numerical methods for capillary action of unsaturated soil subgrade [J].Rock and Soil Mechanics,2013,34(2):421-427.

[3]董斌,張喜發(fā),李欣，等.毛細水上升高度綜合試驗研究[J].巖土工程學報,2008,30(10):1569-1574.

DONG Bin,ZHANG Xifa,LI Xin,et al.Comprehensive tests on rising height of capillary water[J].Chinese Journal of Geotechnical Engineering,2008,30(10):1569-1574.

[4]ARYA L M,PARIS J F.A physicoempirical model to predict the soil moisture characteristic from particle-size distribution and bulk density data[J].Soil Science Society of America Journal,1981,45(6):1023-1030.

[5]GRIFFITHS F J,GOSHI R C.Changes in pore size distribution due to consolidation of clays[J].Ceotechnique，1989,39(1):159-167.

[6]張先偉,孔令偉,郭愛國，等.基于SEM和MIP試驗結構性黏土壓縮過程中微觀孔隙的變化規(guī)律[J].巖石力學與工程學報,2012,31(2):406-412.

ZHANG Xianwei,KONG Lingwei,GUO Aiguo,et al.Evolution of microscopic pore of structured clay in compression procession based on SEM and MIP test[J].Chinese Journal of Geotechnical Engineering,2012,31(2):406-412.

[7]梁健偉,房營光.極細顆粒黏土滲流特性試驗研究[J].巖石力學與工程學報,2010,29(6):1222-1230.

LIANG Jianwei,FANG Yingguang.Experimental study of seepage characteristics of tiny-particle clay[J].Chinese Journal of Geotechnical Engineering,2010,29(6):1222-1230.

[8]陶高梁,張季如,莊心善，等.壓縮變形影響下的土-水特征曲線及其簡化表征方法[J].水利學報,2014,45(10):1239-1246.

TAO Gaoliang,ZHANG Jiru,ZHUANG Xinshan,et al.Influence of compression deformation on the soil-water characteristic curve and its simplified representation method[J].Journal of Hydraulic Engineering,2014,45(10):1239-1246.

[9]孫德安,高游,劉文捷，等.紅黏土的土水特性及其孔隙分布[J].巖土工程學報,2015,37(2):351-356.

SUN Dean,GAO You,LIU Wenjie,et al.Soil-water characteristic and pore-size distribution of lateritic clay[J].Chinese Journal of Geotechnical Engineering,2015,37(2):351-356.

[10]宋盛淵,王清,孫鐵，等.不同地基處理條件下吹填土孔隙分布特征[J].長安大學學報(自然科學版),2015,35(2):46-51.

SONG Shengyuan,WANG Qing,SUN Tie,et al.Pore distribution characteristics of dredger fill consolidated by different methods.[J].Journal of Chang’an Univeersity (Natural Science Edition),2015,35(2):46-51.

[11]盧寧，力科思.非飽和土力學[M].韋昌富,侯龍,簡文星，譯，北京:高等教育出版社,2012.

[12]FAYBISHENKO B A.Hydraulic behavior of quasi-saturated soils in the presence of entrapped air:Laboratory experiments [J].Water Resources Research,1995,31(10):2421-2435．

[13]張華,吳爭光.封閉氣泡對一維積水入滲影響的試驗研究[J].巖土力學,2009,30(增刊2):132-137,148.

ZHANG Hua,WU Zhengguang.Experimental research on the effect of entrapped air on one-dimensional ponded water infiltration[J].Rock and Soil Mechanics,2009,30(S2):132-137,148.

[14]談云志,孔令偉,郭愛國,等.壓實紅黏土水分傳輸的毛細效應與數值模擬[J].巖土力學,2010,31(7):2289-2294.

TAN Yunzhi,KONG Lingwei,GUO Aiguo,et al.Capillary effect of moisture transfer and its numerical simulation of compacted laterite soil[J].Rock and Soil Mechanics,2010,31(7):2289-2294.

[15]薛海斌,黨發(fā)寧,尹小濤,等.土水特征曲線形狀對降雨條件下土質邊坡穩(wěn)定性的影響研究[J].西安理工大學學報,2014,30(3):253-260,250.

XUE Haibin,DANG Faning,YIN Xiaotao,et al.Research on the effect of shape of soil-water characteristic curve on soil slope stability under rainfall condition[J].Journal of Xi’an university of Technology,2014,30(3):253-260,250.

(責任編輯周蓓)

Research on the algorithm of equivalent aperture of clay particles

GENG Daxin1，WANG Yingying1，LI Yuhan2

(1.Geotechnical Research Institute from School of Civil Engineering and Architecture,East China Jiaotong University,Nanchang 330013,China；2.Geotechnical Research Institute,Nanjing Hydraulic Research Institute,Nanjing 210029,China)

The characteristics of randomicity and complexity about pores have a significant influence on the engineering properties of unsaturated soil.In order to analyze soil from microcosm,remolded red clay in Jiangxi Province with different clay contents are taken as research objects.The experiment on soil-water characteristic and experiment on one-dimensional capillary are conducted.The distribution of pores is obtained,and the algorithm of equivalent aperture of clay particles is given.Moreover,the method of incremental step and detailed algorithm progress are provided.The rationality of the formula is confirmed,a realization of the quantification of micro-structure.Meanwhile an effective method is provided for the calculation of the maximum height of capillary rise.

unsaturated soil; the equivalent aperture; capillary rise

10.19322/j.cnki.issn.1006-4710.2016.01.016

2015-09-08

江西省自然科學基金資助項目(20142BAB206002)

耿大新，男，博士，副教授，研究方向為巖土工程、隧道工程等。E-mail：gengdaxin@ecjtu.jx.cn

TU4

A

1006-4710(2016)01-0091-05