顆粒形態(tài)對粗粒土滲透系數(shù)影響的數(shù)值模擬研究

駱莉莎，周 昕，林 軍

（江蘇開放大學，江蘇 南京，210000）

1 引言

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

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

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

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

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

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

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

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

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

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

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

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

圖2 多面體顆粒的形狀

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

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

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

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

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

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

圖3 孔隙流體壓力分布

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

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

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

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

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

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

4 結(jié)論

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