駱莉莎，周 昕，林 軍

（江蘇開放大學(xué)，江蘇 南京，210000）

1 引言

在地下空間開發(fā)中，工程降水是必不可少的工程環(huán)節(jié)；在于土石壩工程、邊坡工程、道路工程中，土體滲流特性也是關(guān)系工程穩(wěn)定性和安全性的關(guān)鍵因素。土體滲透系數(shù)是描述地下水在土體孔隙中滲流特性的關(guān)鍵參數(shù)，也是決定工程降水設(shè)計(jì)的關(guān)鍵。目前確定土體滲透系數(shù)的方法包括經(jīng)驗(yàn)法、室內(nèi)試驗(yàn)法，原位抽水試驗(yàn)法。由于影響

土滲透性的因素非常多，經(jīng)驗(yàn)法往往難以準(zhǔn)確獲得滲透系數(shù)；室內(nèi)試驗(yàn)法很難制備與原位場(chǎng)地具有相同微觀結(jié)構(gòu)的土樣，因而試驗(yàn)結(jié)果準(zhǔn)確性不夠，而且對(duì)于粗粒土試樣必須足夠大，常規(guī)試驗(yàn)設(shè)備不滿足要求。原位抽水試驗(yàn)法是確定現(xiàn)場(chǎng)滲透系數(shù)的最可靠方法，然而由于對(duì)地下地層結(jié)構(gòu)掌握不全面，地下水連通性無(wú)法準(zhǔn)確揭示，根據(jù)抽水試驗(yàn)得到的滲透系數(shù)也存在影響因素多、準(zhǔn)確性無(wú)法保證的缺陷。

對(duì)影響土體滲透性的各因素的作用機(jī)理開展理論研究是當(dāng)前研究的重要方面。李云川等[1]基于計(jì)算流體動(dòng)力學(xué)與離散元耦合的數(shù)值模擬方法，研究了土壤滲透系數(shù)隨土層深度的變化規(guī)律。焦浩然等[2]從室內(nèi)試驗(yàn)和數(shù)值模擬兩個(gè)方面探討了細(xì)粒土在考慮應(yīng)力作用下滲透特性研究進(jìn)展。李華等[3]通過(guò)試驗(yàn)發(fā)現(xiàn)黃土中孔徑為1.3μm

上述研究表明，土體的滲透受很多因素復(fù)雜影響，尋找定量的預(yù)測(cè)模型成為理論研究實(shí)用化的關(guān)鍵步驟。例如，包孟碟等[7]研究了級(jí)配與滲透系數(shù)的關(guān)系，建立了考慮級(jí)配曲線面積的滲透系數(shù)經(jīng)驗(yàn)公式；曹志翔等[8]根據(jù)水頭損失細(xì)觀機(jī)制推導(dǎo)了滲透系數(shù)理論計(jì)算公式，公式以孔隙率和等效粒徑為基本參數(shù)，具有明確的物理意義。Liu 和Jeng[9]以及劉一飛等[10]模擬了不同形態(tài)顆粒構(gòu)成的土體的滲透系數(shù)，并修正了Kozeny-Carman 滲透系數(shù)公式；饒?jiān)瓶档萚11]基于近百組粗粒土的滲透系數(shù)試驗(yàn)數(shù)據(jù)，以全級(jí)配和土樣孔隙比作為神經(jīng)網(wǎng)絡(luò)的輸入變量，利用遺傳算法訓(xùn)練得到了粗粒土滲透系數(shù)預(yù)測(cè)模型。

上述擬合公式考慮的因素各有差異，預(yù)測(cè)結(jié)果精度各有不同，這是由于各個(gè)影響因素往往具有耦合影響。整體而言，目前對(duì)土體滲透性與顆粒形狀關(guān)系的系統(tǒng)研究還較少，綜合考慮各因素相互耦合的作用成果也較少。為此，本文采用數(shù)值模擬方法，探討顆粒形狀以及密實(shí)程度對(duì)粗粒土滲透特性的影響規(guī)律，以促進(jìn)對(duì)土體滲透性的理論認(rèn)識(shí)深度。

2 數(shù)值模型與模擬方法

首先建立孔隙比為0.69 的球形顆粒和多面體顆粒集合體（相當(dāng)于數(shù)值土體試樣），如圖1 所示。其中，孔隙比的計(jì)算公式為

式中，V 為顆粒集合體的總體積，Vs為所有顆粒體積之和。

圖1 數(shù)值土體試樣

顆粒集合體試樣中，球形顆粒粒徑在1.8mm到2.2mm 范圍均勻分布。多面體顆粒形狀共有6種，如圖2 所示，這些顆粒的最長(zhǎng)軸與最短軸之比在1 到1.4 之間。多面體顆粒的等效粒徑為

式中，Vp為多面體顆粒的體積。

為與球形顆粒對(duì)比，多面體顆粒集合體中的顆粒等效粒徑也在1.8mm 到2.2mm 范圍均勻分布，6 種形狀顆粒的數(shù)量相同。

圖2 多面體顆粒的形狀

圖1 中同時(shí)給出了試樣在三個(gè)正交方向的切面?？梢钥吹?，在每個(gè)切面上固體部分（顆粒切面）排列非常不規(guī)律，使得流體在流經(jīng)這些切面時(shí)受到不同程度的阻礙。此外，同一個(gè)試樣三個(gè)切面上的顆粒排列不盡相同。

采用Avizo 軟件的滲透特性模擬（Absolute Permeability Experiment Simulation）功能測(cè)試兩個(gè)試樣的沿三個(gè)正交方向的滲透系數(shù)。模擬方法為在兩相對(duì)面上分別施加不同的孔隙水壓力，兩相對(duì)面上孔隙水壓相差30kPa，其余四個(gè)面為不透水邊界。通過(guò)模擬一定時(shí)間內(nèi)由高水壓面向低水壓面流動(dòng)的水流量，即可由下式估算滲透系數(shù)：

式中，Q 為體積流速，?w為水的重度（10kN/m3），L 為流動(dòng)距離，S 為垂直流向的橫截面積， ?p為壓力差（本文中為30kPa）。

孔隙穩(wěn)態(tài)流場(chǎng)通過(guò)有限體積法求解簡(jiǎn)化的NS 方程得到，假定流體不可壓縮，且為牛頓流體，處于層流狀態(tài)，流體和固體界面不可滑動(dòng)。

3 模擬結(jié)果與討論

圖3 為穩(wěn)定流動(dòng)時(shí)，孔隙內(nèi)流體壓力的空間分布。整體上，孔隙流體壓力沿流動(dòng)方向降低，在同一垂直流動(dòng)方向的界面上，孔隙壓力分布并不均勻；孔隙連通性較好、顆粒稀疏的局部區(qū)域，流動(dòng)阻力較小，孔隙壓力較低，是流體滲透的優(yōu)勢(shì)通道。而對(duì)于孔隙連通性較差、顆粒密集的區(qū)域，流動(dòng)阻力則較大，孔隙壓力較大。盡管多面體顆粒的等效粒徑與球形顆粒一致，但從圖1 及圖3 可以看到，多面體顆粒試樣中的孔隙被顆粒分割得更加零碎，流通孔道尺寸明顯小于球形顆粒試樣，定性上預(yù)示著流體更難通過(guò)多面體顆粒試樣，即它的滲透系數(shù)更小。

圖3 孔隙流體壓力分布

數(shù)值模擬得到的三個(gè)方向滲透系數(shù)結(jié)果如表1 所示。Kozeny-Carman 公式可用于估算球形顆粒集合體的滲透系數(shù)：

式中，? 為動(dòng)力粘滯系數(shù)（0.001Pa·s），d 為平均粒徑（本文取2mm），n 為孔隙率，可由孔隙比算得：

表1 數(shù)值模擬結(jié)果與理論解

表1 中數(shù)值模擬的滲透系數(shù)在礫石、粗砂的常見滲透系數(shù)范圍內(nèi)，說(shuō)明了數(shù)值模擬的可靠性。從表1 可以看出，球形顆粒集合體三個(gè)正交方向的滲透系數(shù)相差不大，且與理論解結(jié)果基本吻合。而多面體顆粒集合體三個(gè)正交方向的滲透系數(shù)相差較大，這主要是由于多面體顆粒的排列整體上具有一定各向異性，特別是長(zhǎng)、短軸差別較大的顆粒，其取向?qū)B透系數(shù)影響較為明顯。此外，多面體顆粒集合體的滲透系數(shù)僅為球形顆粒集合體滲透系數(shù)的40%；這是由于多面體顆粒試樣中的孔隙被顆粒分割得更加零碎，流體在孔隙中流動(dòng)受到的阻力因而更大，滲透性因而降低。

圖4 進(jìn)一步給出了本研究模擬的不同孔隙比試樣三向滲透系數(shù)平均值隨孔隙比的變化。如圖所示，滲透系數(shù)隨孔隙比增加而非線性增長(zhǎng)，增加速率也隨孔隙比增大而增大。球形顆粒集合體的滲透系數(shù)始終小于理論解，但二者差別逐漸減小。多面體顆粒集合體的滲透系數(shù)受孔隙形狀的影響非常明顯，即使孔隙比達(dá)到0.85，其滲透系數(shù)也明顯小于球形顆粒集合體，兩者比值在35%到50%之間。

圖4 滲透系數(shù)隨孔隙比的變化

4 結(jié)論

本文采用數(shù)值模擬方法研究了不同顆粒形狀的粗粒土滲透特性，獲得了滲透系數(shù)數(shù)值解。結(jié)合本文模擬結(jié)果和滲透系數(shù)理論公式可知，粗粒土的滲透系數(shù)與流體粘滯性、土體密實(shí)程度（孔隙比）、顆粒粒徑、顆粒形狀有關(guān)。本文模擬的球形顆粒和多面體顆粒分別代表了磨圓度很好的粗粒土和磨圓度很差（棱角分明）的粗粒土兩種特殊情況。結(jié)果表明，顆粒形狀越復(fù)雜，孔隙形狀的復(fù)雜程度越高，流體滲透越困難，滲透系數(shù)因而越低。多面體顆粒集合體的滲透系數(shù)明顯小于球形顆粒集合體，兩者比值在35%到50%之間。