陳昌昀， 張 鵬， 許國輝??， 劉志欽， 許興北， 任宇鵬(中國海洋大學(xué)1.海洋環(huán)境與生態(tài)教育部重點實驗室; 2.山東省海洋環(huán)境地質(zhì)工程重點實驗室，山東 青島266100)

落球法測定波動液化土粘滯性的試驗研究?

陳昌昀1,2， 張 鵬1,2， 許國輝1,2??， 劉志欽1,2， 許興北1,2， 任宇鵬1,2
(中國海洋大學(xué)1.海洋環(huán)境與生態(tài)教育部重點實驗室; 2.山東省海洋環(huán)境地質(zhì)工程重點實驗室，山東 青島266100)

海底的粉質(zhì)土在波浪作用下發(fā)生液化，液化后的土體產(chǎn)生流體一樣的波動。本文在進行液化粉質(zhì)土波動試驗過程中，利用落球法進行了液化土的黏滯系數(shù)的測試試驗。試驗采用了銅質(zhì)、鐵質(zhì)和氧化鋯等不同密度的光滑球體，通過位移傳感器記錄密度球在液化后土體中的位移和運動時間，換算出密度球在液化土體中沉降的平均速度，利用斯托克斯公式計算給出液化土的黏滯系數(shù)。結(jié)果給出試驗用液化土體的黏滯系數(shù)在10 kPa·s的量級。

土體液化；落球法試驗；黏滯系數(shù)；斯托克斯公式

在波浪和極端天氣如風(fēng)暴潮等的作用下，海床土體會發(fā)生液化破壞，由液化引發(fā)的土體破壞變形是導(dǎo)致工程結(jié)構(gòu)破壞的主要原因之一。因此，國內(nèi)外學(xué)者針對土體液化問題進行了相關(guān)研究。Miyamoto結(jié)合Sassa等通過離心機試驗對波浪荷載作用下海床液化的分析，研究了可液化土層在推進波作用下發(fā)生液化的問題[1-2]。Sumer等在推進波荷載作用下粉質(zhì)海床的波浪響應(yīng)研究中，對海床的液化特點進行了分析[3-4]。李安龍等通過配置不同粒徑的底床，發(fā)現(xiàn)砂土在波浪荷載作用下發(fā)生液化時，若黏粒含量小于10%，則對砂土液化有促進作用[5]。馮秀麗等通過現(xiàn)代黃河三角洲砂質(zhì)粉土和粉質(zhì)粘土觸變性的對比試驗，發(fā)現(xiàn)海底塌陷凹坑是快速沉積的粉砂液化的結(jié)果[6]。

之前的學(xué)者對土體液化的研究大多處于液化前及開始發(fā)生液化的階段, 而對土體液化后的研究相對較少，液化后的土體具有流體的特性[7]，而粘滯系數(shù)可以反映液化后土體的流動性[8]。因此，研究者們通過不同的試驗手段開展了液化土體黏滯性的測試研究。Towhata等利用拖管試驗，通過測定施加在鋼管上的阻力獲得粘性流體的表觀粘度，并發(fā)現(xiàn)當(dāng)砂土完全液化后，拖拽力與拖管的運動速度關(guān)聯(lián)性很大，進一步證實了液化土體的流體性質(zhì)[9]。Kawakami等利用黏度計得到了初始液化后的土體的表觀黏度與剪應(yīng)變率具有一定的對應(yīng)關(guān)系[10]。Hwang等在研究液化地表流對樁基的影響時，通過拖管試驗，發(fā)現(xiàn)液化砂土的黏度隨剪應(yīng)變率升高而降低，在數(shù)值分析中將液化砂土的黏性作為一個重要參數(shù)進行分析[11]。Hadush等在總結(jié)了幾組由不同測量方法獲得的粘度與剪應(yīng)變率關(guān)系后，同樣發(fā)現(xiàn)了液化砂土的這種特性[12]。陳育民等通過拖球試驗測量鋼球所受水平阻力來計算液化及液化后砂土的表觀動力黏度，進一步研究了液化砂土的流動特性[13]。前有研究發(fā)現(xiàn)海床土體在波浪荷載作用下會發(fā)生液化，液化狀態(tài)的土體是可以視為流體的，其特性如同黏性流體[14]，而在對流體特性的研究中黏滯系數(shù)是一個重要的參量。已有的液化土體流動性研究大多是利用振動臺試驗，通過拖球或拖管法針對地震作用下砂土液化后的流動特性開展，而波致液化土體黏滯性測定是需要補充的工作，本文就此開展了相關(guān)研究。海底的粉質(zhì)土在波浪作用下會發(fā)生液化，液化后的粉質(zhì)土產(chǎn)生流體一樣的波動[15]。本文采用測定流體黏滯系數(shù)的經(jīng)典方法——落球法，在液化粉質(zhì)土波動的水槽試驗過程中，利用玻璃、氧化鋁、氧化鋯、鐵、銅等5種不同密度的球體，對液化粉質(zhì)土的黏滯系數(shù)進行測定。

1 試驗方法

試驗用土體是取自黃河三角洲的粉質(zhì)土，其粒徑級配曲線見圖1所示。試驗過程中持續(xù)施加波高15 cm、周期1 s的波浪使土體液化。

1.1 試驗底床及落球裝置設(shè)置

將試驗所用的粉質(zhì)土和水混合攪拌成均勻泥漿，控制其含水率在33%左右。隨后用“順板滑塌法”將混合好的泥漿分批移入水槽鋪設(shè)成尺寸為2 m×0.5 m×0.6 m的底床，將6個直徑為3 cm的表面光滑的不同材質(zhì)的密度球(具體材料及規(guī)格見表1)分成3組，分別放置在底床表層的不同位置處，相互之間距離大于20 cm，避免相互影響(示意圖見圖2)。用0.3 mm的鋼絲線將密度球與拉線式移傳感線相連(由于位移傳感器在拉線抽出時自身存在拉力，試驗前在傳感器下方懸掛一定質(zhì)量的配重物進行拉力平衡)。最后，向水槽內(nèi)加水至床面上40 cm處，靜置固結(jié)10 d。

1.2 落球試驗

靜置固結(jié)后，開始施加波浪作用。首先在底床未擾動的情況下，施加波高為15cm的波浪，持續(xù)加波3h后土體未液化。將底面積10cm×15cm的鋼板錘置于底床上，通過不斷提起放下?lián)舸虻状?，加速土體液化。待底床土體發(fā)生液化時釋放密度球，同時使用奧地利DEWESoft瞬態(tài)數(shù)據(jù)采集器(型號為DEWE-43)，以100Hz的頻率連續(xù)記錄密度球的沉降位移與時間數(shù)據(jù)

直至試驗結(jié)束。密度球沉降結(jié)束后，使用自制土體采樣器插入底床，取不同液化深度的土樣并計算平均密度用于粘滯系數(shù)的計算。設(shè)置的部分密度球隨液化土體運動并出現(xiàn)沉降現(xiàn)象，也有出現(xiàn)停滯的情況，通過再次調(diào)整配重物質(zhì)量來平衡拉線式位移傳感器的拉力，可以使密度球繼續(xù)沉降。

2 試驗結(jié)果

試驗底床受到錘擊擾動發(fā)生液化，并且液化現(xiàn)象在土體的橫向和垂向迅速發(fā)展。通過觀察發(fā)現(xiàn)液化時土顆粒與密度球隨波浪做周期性的橢圓運動，在一個波浪周期內(nèi)密度球向左、向右運動趨勢相同，位移相互抵消，可視為垂直方向的沉降運動。從位移傳感器測得的數(shù)據(jù)中發(fā)現(xiàn)各密度球在土體發(fā)生液化后出現(xiàn)了不同的沉降現(xiàn)象，其中玻璃球和氧化鋁球未發(fā)生沉降，而氧化鋯球、鐵1球、鐵2球和銅球均產(chǎn)生沉降位移，并且各密度球都經(jīng)歷了初期加速下沉，中期勻速沉降，最后沉降停止的過程。在勻速沉降過程中，氧化鋯球、鐵1球、鐵2球和銅球的位移分別是75.7、65.6、71.7、67.0 mm。

氧化鋯球、鐵1球、鐵2球和銅球在液化土中沉降時，受到的阻力逐漸增大到與重力平衡，因此出現(xiàn)了從加速沉降逐漸變?yōu)閯蛩俪两档默F(xiàn)象，而玻璃球和氧化鋁球未發(fā)生沉降的原因是由于兩者的密度較小造成的。選擇密度球沉降均勻的時間段(勻速沉降段)的時間與位移數(shù)據(jù)進行計算分析并繪圖如圖3所示，在該時間段內(nèi)密度球的最小位移是6.56 cm，大于小球直徑2倍以上，圖3中紅線為各密度球沉降位移的標(biāo)準(zhǔn)曲線，發(fā)現(xiàn)氧化鋯球、鐵1球、鐵2球和銅球沉降位移與時間關(guān)系，用線性相關(guān)來擬合，其相關(guān)系數(shù)R2值分別為0.95、0.99、0.95和0.97，可以視為勻速運動。

將氧化鋯球、鐵1球、鐵2球、銅球的下落位移與持續(xù)時間之比作為沉降過程中的平均速度(見表2)，并視為各密度球在液化土體中的沉降速度，則計算得到氧

化鋯球的沉降速度為0.168 mm/s，鐵1球的沉降速度為0.218 mm/s，鐵2球的沉降速度為0.228 mm/s，銅球的沉降速度為0.335 mm/s。幾種密度球在試驗土體中處于靜止?fàn)顟B(tài)時沒有沉降下落，當(dāng)土體發(fā)生液化波動時，密度球的沉降速度區(qū)間為0.1～0.4 mm/s，極為緩慢，并且各密度球與液化土體之間的相對運動發(fā)生在試驗土體的黏滯狀態(tài)下，屬于層流區(qū)的范圍，因此，斯托克斯公式適用于本次試驗。

3 討論

本次試驗中底床土體的空間尺寸(見圖2)相對于密度球而言可視為無限廣延的，球體是光滑且剛性的，測試值是小球在不產(chǎn)生渦流的流體中受到的阻力。試驗中小球在自身重力作用下自由沉降，并逐漸達到勻速沉降階段，此時可視為自然狀態(tài)下的受力平衡，排除了人為或機械原因?qū)π∏蜻\動狀態(tài)造成的影響，符合斯托克斯公式的適用條件。因此，利用本文試驗方法測試得到密度球的沉速后，利用斯托克斯公式計算得到液化土體的黏滯系數(shù)是合理的。

從表3可知，在水槽試驗過程中, 利用落球法測得的波動液化粉質(zhì)土的黏滯系數(shù)在9.67～13.47 kPa·s的范圍。劉濤等在利用拖球試驗對波動液化粉土的流動性研究中，得到密度球在速度為25和33 mm/s時，測得液化粉土的黏滯系數(shù)分別為60.45和65.48kPa·s[16]。本次試驗獲得的粉質(zhì)土黏滯系數(shù)與之相比，量級雖均為10 kPa·s，但數(shù)據(jù)存在一定差異。為分析造成數(shù)據(jù)差異的原因，制備了含水率為35%左右的粉質(zhì)土，并將直徑為30 mm的鐵球埋置其中，在豎直方向上以20 mm·s-1的速度勻速拉動鐵球，發(fā)現(xiàn)在拉力方向的球面處有一定厚度的土體附著集聚，且與鐵球形成一個新的整體共同運動，這可能是導(dǎo)致拖球試驗中粘滯系數(shù)計算結(jié)果偏大的原因。本文落球法為球體依靠重力的自然下沉，沉降速度不到1 mm/s，能夠保證液化土體以繞流形式運動，更加符合斯托克斯公式的應(yīng)用條件。

4 結(jié)論

本文在開展波浪作用下粉質(zhì)土底床的液化試驗過程中，利用落球法開展了液化粉質(zhì)土粘滯性的試驗測試，并用斯托克斯公式計算了液化土體的粘滯系數(shù)，給出如下結(jié)論：

(1)粉土在波浪荷載作用下發(fā)生液化時，可以通過落球法試驗測得土體的黏滯系數(shù)。

(2)試驗獲得落球的沉降速度后，可以通過斯托克斯公式計算液化粉質(zhì)土的黏滯系數(shù)。

(3)波浪水槽試驗中液化粉質(zhì)土的黏滯系數(shù)值在9.67～13.47 kPa·s，量級為10 kPa·s。

[1] Miyamoto J, Sassa S, Sekiguchi H．Progressive solidification of a liquefied sand layer during continued wave loading [J]．Geotechnique， 2004， 54(10)： 617-629．

[2] Sassa S, Sekiguchi H. Wave-induced liquefaction of beds of sand in a centrifuge[J]．Geotechnique， 1999， 49(5)： 621-638．

[3] Sumer B M，F(xiàn)reds?e J，Hatipoglu F，et al．The sequence of sediment behaviour during wave-induced liquefaction[J]．Sedimentology， 2006， 53： 611-629．

[4] Sumer B M，F(xiàn)reds?e J，Christensen S，et al．Sinking / floatation of pipelines and other objects in liquefied soil under waves[J]． Coastal Engineering， 1999， 38(2)： 53-90．

[5] 李安龍，楊榮民，楊少麗，等. 波浪加載下海底土質(zhì)特性變化的研究[J] .青島海洋大學(xué)學(xué)報， 2003， 33(1): 101-106． Li Anlong，Yang Rongmin，Yang Shaoli，et al． Study on the geotechnical behavior seabottom sediments under wave-loading[J]．Journal of Ocean University of Qingdao， 2003， 33(1): 101-106．

[6] 馮秀麗，周松望，林霖．現(xiàn)代黃河三角洲粉土觸變性研究及其應(yīng)用[J]．中國海洋大學(xué)學(xué)報(自然科學(xué)版)， 2004， 34(6)： 1053-1056． Feng Xiuli， Zhou Songwang，Lin lin，et al．The thixotropy of silt in Huanghe delta [J]．Periodical of Oecan Univercity of China， 2004， 34(6)： 1053-1056．

[7] 黃雨，毛無衛(wèi)．液化后土體的流態(tài)化特性研究進展[J]．同濟大學(xué)學(xué)報(自然科學(xué)版)， 2011， 39(4)： 501-506． Huang Yu，Mao Wuwei．State of art of fluidization behavior of post-liquefied soils [J]．Journal of Tongji University (Natural Science)， 2011， 39(4)： 501-506．

[8] 周恩全，王志華，陳國興，等．飽和砂土液化后流體本構(gòu)模型研究[J]．巖土工程學(xué)報， 2015， 37(1)： 112-118． Zhou Enquan，Wang Zhihua，Chen Guoxing，et al. Constitutive model for fluid of post-liquefied sand [J]. Chinese Journal of Geotechnical Engineering， 2015， 37(1)： 112-118．

[9] Towhata I, Vargas—Monge W, Orense R P, et a1.Shaking table tests on subgrade reaction of pipe embeded in sandy liquefied subsoil[J]. Soil Dynamics and Earthquake Engineering，1999，18：347-361．

[10] Kawakami T, Suemasa N, Hamada M, et al. Experimental study on mechanical properties of liquefied sand[C]. // Salt Lake City: Proceedings of the 5th US-Japan Workshop on Earthquake Resistant Design of Lifeline Facilities and Countermeasures Against Soil Liquefaction． 1994：285-299．

[11] Hwang J I，Kim C Y，Chung C K，et al．Viscous fluid characteristics of liquefied soils and behavior of piles subjected to flow of liquefied soils[J]. Soil Dynamics and Earthquake Engineering， 2006， 26(2-4)： 313-323．

[12] Hadush S，Yashima A，Uzuoka R. Importance of viscous fluid characteristics in liquefaction induced lateral spreading analysis [J]．Computers and Geotechnics， 2000， 27(3)： 199-224．

[13] 陳育民，劉漢龍，邵國建，等．砂土液化及液化后流動特性試驗研究[J]． 巖土工程學(xué)報， 2009(9)： 1408-1413． Chen Yumin，Liu Hanlong，Shao Guojian， et al．Laboratory tests on flow characteristics of liquefied and post-liquefied sand [J]．Chinese Journal of Geotechnical Engineering， 2009(9)： 1408-1413．

[14] Miyajima M，Kitaura M，Koike T，et al. Experimental study on characteristics of liquefied ground flow[C]. // A A Balkema: The First International Conference on Earthquake Geotechnical Engineering, Rotterdam, 1995：969-974.

[15] 許國輝，孫永福，于月倩，等．黃河水下三角洲淺表土體的風(fēng)暴液化問題探討[J]．海洋地質(zhì)與第四紀(jì)地質(zhì)， 2011， 31(2)： 37-42． Xu Guohui, Sun Yongfu, Yu Yueqian, et al. Storm-induced liquefaction of the surficial sediments in the Yellow River subaqueous delta [J]．Maine Geology and Quaternary Geology， 2011， 31(2)： 37-42．

[16] 劉濤，崔逢，張美鑫．波浪作用下液化粉土流動特性拖球試驗研究[J]．海洋學(xué)報， 2016， 38(2): 123-130． Liu Tao，Cui Feng，Zhang Meixin. Dragging ball test on flow characteristics of liquefied silt under wave loading [J]. Acta Oceanologica Sinica， 2016， 38(2): 123-130．

責(zé)任編輯 徐 環(huán)

The Experimental Study on Viscosity of the Wave-Induced Liquefaction Silty Soil by Falling Ball Method

CHEN Chang-Yun1,2， ZHANG Peng1,2，XU Guo-Hui1,2，LIU Zhi-Qin1,2，XU Xing-Bei1,2，REN Yu-Peng1,2

(Ocean University of China, 1. Key Laboratory of Marine Environment and Ecology, Ministry of Education; 2. Shandong Provincial Key Laboratory of Marine Environment and Geological Engineering, Qingdao 266100, China)

Liquefaction of submarine slity soil occurs under the wave action, while the liquefied soils would fluctuate like liquids. In this study, we used the falling ball method to measure the viscosity coefficient of liquefied soil during liquefied slity-soil fluctuation tests. By using density-variable smooth balls (made of Copper iron, zirconia), we used a to record the displacement and migration time of the balls in the liquefied soils, and converted the average ball falling speeds in the liquefied soils. Then the viscosity coefficients of liquefied soils were determined from the Stokes formula. The viscosity coefficient of liquefied soils is at the level of 10 kPa·s.

soil liquefaction; falling ball method; coefficient of viscosity; stokes formula

國家自然科學(xué)基金項目(41576039)資助 Supported by National Natural Science Foundation of China(41576039)

2016-01-20；

2016-05-16

陳昌昀(1989-)，男，碩士生，從事海洋工程地質(zhì)方面研究。 E-mail：ccyhuanjing@163.com

?? 通訊作者：E-mail：xuguohui@ouc.edu.cn

P736

A

1672-5174(2017)06-119-05

10.16441/j.cnki.hdxb.20160011

陳昌昀，張鵬，許國輝，等.落球法測定波動液化土粘滯性的試驗研究[J].中國海洋大學(xué)學(xué)報(自然科學(xué)版),2017, 47(6): 119-123.

CHEN Chang-Yun， ZHANG Peng，XU Guo-Hui,et al.The experimental study on viscosity of the wave-induced liquefaction silty soil by falling ball method[J].Periodical of Ocean University of China, 2017, 47(6): 119-123.