顆粒形狀對多孔介質(zhì)孔隙特征和滲流規(guī)律影響研究的探討

張家發(fā)，焦赳赳

顆粒形狀對多孔介質(zhì)孔隙特征和滲流規(guī)律影響研究的探討

張家發(fā)1，焦赳赳2

（1．長江科學(xué)院水利部巖土力學(xué)與工程重點(diǎn)實(shí)驗(yàn)室，武漢 430010；2．香港大學(xué)地球科學(xué)系，香港）

討論了顆粒形狀對于多孔介質(zhì)孔隙系統(tǒng)及其滲流規(guī)律的可能影響及其研究意義；綜述相關(guān)領(lǐng)域的研究成果，探討了顆粒形狀的測量評定方法、顆粒形態(tài)和填筑的模擬方法、孔隙系統(tǒng)的可視化研究方法，以及滲流規(guī)律的試驗(yàn)對比研究方法；分析指出了顆粒形狀對滲流規(guī)律影響的研究在促進(jìn)水處理濾料開發(fā)和碎石料作為大壩填料利用方面的應(yīng)用前景。

顆粒形狀；多孔介質(zhì)；滲流；試驗(yàn)；模擬

1 概 述

J．Bear在其經(jīng)典著作［1］中指出：固結(jié)物質(zhì)的孔隙率主要取決于膠結(jié)程度；而非固結(jié)物質(zhì)的孔隙率則依賴于顆粒的形狀、粒徑分布和顆粒的排列方式。土體是典型的非固結(jié)多孔介質(zhì)。長期以來的滲流研究大多有意無意地忽視了介質(zhì)顆粒形狀的影響，而主要關(guān)注顆粒的大小、級配組成、排列方式的影響。雖然土體顆粒以具有一定的磨圓度為常見，可以忽略其形狀的影響，但是仍然存在著磨圓度差，甚至未經(jīng)磨蝕而棱角分明、形態(tài)各異的土體顆粒，冰磧土、崩塌堆積物以及日益廣泛利用的巖石爆破料就是典型。

鄭穎人指出［2］，巖土類材料是由顆粒材料堆積或膠結(jié)而成，屬摩擦型材料。摩擦材料的特點(diǎn)是抗剪強(qiáng)度中含有摩擦力項(xiàng)，摩擦力項(xiàng)就與顆粒的力學(xué)性質(zhì)、大小、形狀及排列關(guān)系直接相關(guān)。顆粒形狀會影響顆粒之間的接觸關(guān)系和填充關(guān)系，當(dāng)然也就會影響孔隙的形狀和大小，從而影響到其中的水流運(yùn)動(dòng)。

劉杰曾經(jīng)通過比較發(fā)現(xiàn)砂礫石料和人工碎石料在相同級配條件下的壓實(shí)性及滲透性的差異［3］，說明了顆粒形狀會對土體宏觀滲流特性造成影響，但是沒有從機(jī)理上和微細(xì)觀層面進(jìn)行解釋。M．Stew- art等［4］用格子Boltzmann方法模擬扁橢球顆粒填筑體中的水流運(yùn)動(dòng)，發(fā)現(xiàn)顆粒三維尺寸之間的比值對于顆粒的排列結(jié)果及滲透張量的各向異性程度都有顯著影響。這些研究都表明，在一定條件下顆粒形狀對于多孔介質(zhì)滲流規(guī)律的影響可能是顯著的。

本文將討論研究顆粒形狀對多孔介質(zhì)滲流規(guī)律的意義，綜述已有的相關(guān)研究成果，探討顆粒形狀對多孔介質(zhì)滲流規(guī)律影響的研究方法和途徑。

2 顆粒形狀對多孔介質(zhì)滲流影響的研究意義

多孔介質(zhì)滲流是自然界廣泛存在的現(xiàn)象和過程。J．Bear對多孔介質(zhì)定義進(jìn)行了詳細(xì)的討論［1］。本文將多孔介質(zhì)簡單定義為：固體構(gòu)成復(fù)雜孔隙系統(tǒng)的骨架，流體可以在連通的孔隙系統(tǒng)中運(yùn)動(dòng)的物質(zhì)空間。除了環(huán)境因素和邊界條件外，多孔介質(zhì)滲流規(guī)律的主要影響因素是流體和多孔介質(zhì)雙方自身的性質(zhì)。人類生活和工程建設(shè)中經(jīng)常遇到多孔介質(zhì)滲流作用帶來的挑戰(zhàn)，例如大壩及其基礎(chǔ)的滲流；或者有意利用多孔介質(zhì)滲流的作用，例如多孔陶瓷和濾料的運(yùn)用。在地下水和石油資源開采、農(nóng)田水利、地下水環(huán)境、地下空間的利用、自然地質(zhì)災(zāi)害、文物保護(hù)、甚至生物滲流等很多領(lǐng)域也利用多孔介質(zhì)滲流理論來研究解決自己遇到的問題［5］。

隨著人類社會的進(jìn)化和發(fā)展，人類逐漸偏好于在濱海、平原和河流中下游地區(qū)活動(dòng)，所接觸的土體大多是經(jīng)過長距離搬運(yùn)或長期磨蝕后的顆粒沉積物。法國工程師達(dá)西于1856年運(yùn)用石英砂柱開展水力試驗(yàn)［6］，揭示了滲流速度與滲流比降成線性關(guān)系的規(guī)律，這就是如今被廣泛運(yùn)用的達(dá)西定律。達(dá)西試驗(yàn)所采用的石英砂來自圣納河。達(dá)西只給出了其級配組成，沒有描述顆粒的形狀。但是，石英是耐磨的堅(jiān)硬礦物，石英砂是巖石碎屑物經(jīng)歷長距離搬運(yùn)和磨圓作用后在河流中下游形成的典型沉積物。

隨著我國水利水電工程建設(shè)事業(yè)的發(fā)展，高山峽谷地區(qū)的工程建設(shè)項(xiàng)目越來越多。這些地區(qū)都很偏遠(yuǎn)，水泥等工業(yè)材料的大量利用要付出昂貴的運(yùn)費(fèi)，采用當(dāng)?shù)夭牧暇哂忻黠@的成本優(yōu)勢。而在這些地區(qū)，動(dòng)力地質(zhì)作用以剝蝕為主，就地堆積的土體（如崩塌堆積物或冰磧土）沒有經(jīng)過長距離的搬運(yùn)，磨圓度差。這些磨圓度差的堆積物有時(shí)還在性能或數(shù)量上不能滿足工程的要求，于是巖石爆破料等人工碎石料得到越來越廣泛的運(yùn)用。碎石料是典型的異形顆粒，不僅沒有磨圓度，而且棱角分明，形態(tài)各異。

顆粒形狀可以分為球體、圓柱體、橢球體和多面體。用于過濾的化工材料大多制作為球體。自然界的土顆粒大多數(shù)是橢球體，磨圓度越高，越接近球體。完全沒有經(jīng)過磨蝕的土顆?；蛘咚槭瘎t是多面體。

幾何上定義多面體為4個(gè)或4個(gè)以上多邊形所圍成的立體。其中正多面體屬于規(guī)則形狀，其各個(gè)面都是全等的正多邊形，并且各個(gè)多面角都是全等的多面角。正多面體的種數(shù)很少，只有正四面體、正六面體、正八面體、正十二面體、正二十面體5種。不規(guī)則多面體則無數(shù)。將不規(guī)則多面體顆粒稱為異形顆粒。

球形顆粒多孔介質(zhì)孔隙系統(tǒng)的影響因素主要包括顆粒的大小及級配、排列方式和緊密程度。非球形顆粒形成的多孔介質(zhì)孔隙系統(tǒng)則更加復(fù)雜。宏觀上，孔隙性取決于顆粒級配及其排列緊密程度；細(xì)觀和微觀上看，孔隙的形狀和大小與顆粒之間的接觸關(guān)系、填充程度、顆粒形狀，甚至顆粒表面粗糙度等多種復(fù)雜因素有關(guān)。

直觀地理解，異形顆粒更容易架空，其排列關(guān)系和形成的孔隙形態(tài)更具有多樣性，孔隙形態(tài)和大小變異性更大；同時(shí)，異形顆粒的比表面積大于球形顆粒，在其多孔介質(zhì)中固體與流體之間的界面接觸作用更顯著。這些特點(diǎn)毫無疑問會影響流體在孔隙系統(tǒng)中的運(yùn)動(dòng)，從微觀上影響到孔隙系統(tǒng)中流體質(zhì)點(diǎn)運(yùn)動(dòng)的軌跡、能量耗散及固液相的相互作用，從宏觀上影響到多孔介質(zhì)的滲透性、滲透變形特性，以及能夠發(fā)揮的過濾和吸附作用等。這就意味著在某些條件下，利用已經(jīng)掌握的多孔介質(zhì)滲流規(guī)律和理論無法準(zhǔn)確地反映異形顆粒多孔介質(zhì)滲流規(guī)律，甚至?xí)碇囌`。

開展顆粒形狀對多孔介質(zhì)滲流規(guī)律影響研究，將有助于拓展對滲流規(guī)律的認(rèn)識，豐富和完善多孔介質(zhì)滲流理論，也將為堆石壩有關(guān)滲流問題的解決提供更堅(jiān)實(shí)的理論基礎(chǔ)，同時(shí)也可以為工業(yè)和環(huán)境工程中開發(fā)多孔介質(zhì)材料和有效利用多孔介質(zhì)的特性提供指導(dǎo)，因而具有理論意義和實(shí)際應(yīng)用前景。

3 顆粒形狀及異形顆粒多孔介質(zhì)滲流研究方法

3．1 顆粒形狀的測量和評定

異形顆粒的形狀可能很復(fù)雜，顆粒形狀的模擬和量化是研究異形顆粒多孔介質(zhì)孔隙系統(tǒng)及其滲流規(guī)律的基礎(chǔ)。

機(jī)械工程領(lǐng)域常用的圓度誤差評定方法可借以刻畫任意顆粒的現(xiàn)狀。圓度誤差是指回轉(zhuǎn)體的同一正截面上實(shí)際被測輪廓對其理想圓的變動(dòng)量，理想圓的位置應(yīng)符合最小條件［7］。目前常用的有最小區(qū)域圓、最小二乘圓、最小外接圓和最大內(nèi)切圓4種評定方法。在測量數(shù)據(jù)的基礎(chǔ)上用這些方法可以對任意顆粒的復(fù)雜幾何形態(tài)進(jìn)行量化。4種方法都有其適用條件，相同條件下評定結(jié)果也可能有一定差異［8］，要根據(jù)具體情況研究選用。

在相同的顆粒密度和物質(zhì)空間基礎(chǔ)上，不同形狀顆粒多孔介質(zhì)孔隙的對比可以建立在相同的固體體積或者質(zhì)量基礎(chǔ)上，也就是在相同孔隙率基礎(chǔ)上開展孔隙特征的對比。

利用規(guī)則顆粒多孔介質(zhì)與球形顆粒多孔介質(zhì)開展對比研究，可以初步揭示顆粒形狀對滲流規(guī)律的影響，同時(shí)也可以用于指導(dǎo)人工濾料的研發(fā)，從而兼具理論和現(xiàn)實(shí)意義。在相同體積條件下，正四面體、立方體、正八面體、正十二面體、正二十面體的比表面積與球體比表面積的比例分別為1．49，1．24，1．18，1．10，1．06。所以如果顆粒之間表面不契合，多面體顆粒多孔介質(zhì)比球形顆粒多孔介質(zhì)與水流接觸更充分，界面作用更加顯著。

由于粘性土本身具有結(jié)構(gòu)性，討論顆粒形狀影響的意義不大。對于無粘性土，土力學(xué)中常常是用粒徑來表達(dá)顆粒大小，通過篩分試驗(yàn)，用顆粒粒徑分布曲線來表達(dá)土體的顆粒組成，也稱為級配。但是非球形顆粒在篩分過程中是由其最小截面的尺寸決定其所在粒徑組的。正是由于粒徑組的劃分，篩分試驗(yàn)結(jié)果只具有較粗略的統(tǒng)計(jì)意義。即使是在顆粒形狀相同的條件下，相同級配線的寬級配土料之間，顆粒大小也不具有嚴(yán)格的對等關(guān)系，這是干擾土的力學(xué)和滲流試驗(yàn)可重復(fù)性的重要因素之一。對于寬級配天然土料來說，篩分試驗(yàn)造成的顆粒粒徑上的誤差會在很大程度上掩蓋或者干擾顆粒形狀差異帶來的影響。另一方面，大量較細(xì)土顆粒形狀的測量和評定工作量太大，以至于缺乏可操作性。所以，宜從窄級配粗粒料著手，從上述圓度誤差評定方法中選取合適的方法，研究土顆粒形狀，以便進(jìn)一步比較其對滲流規(guī)律的影響。

3．2 顆粒形態(tài)及填筑的模擬研究

異形顆粒形態(tài)的模擬是實(shí)現(xiàn)異形顆粒多孔介質(zhì)孔隙系統(tǒng)模擬的基礎(chǔ)。史旦達(dá)等［9］采用角粒團(tuán)顆粒方法在顆粒流程序（PFC2D／3D）基礎(chǔ)上開發(fā)了非純圓顆粒形狀的模擬功能。Jia X．等提出用顆粒成像的像素集來表達(dá)任意形狀的顆粒［10］。張?jiān)戚x等［11］也提出利用物體的數(shù)字圖像進(jìn)行邊緣檢測，取得圓度測量數(shù)據(jù)。這些方法對顆粒形狀的模擬具有重要借鑒意義和參考價(jià)值。

顆粒填筑（grain packing）是計(jì)算機(jī)模擬領(lǐng)域的熱點(diǎn)研究問題，主要是研究物體在有限空間內(nèi)互不重疊的優(yōu)良布局［12］。已經(jīng)有研究人員將其應(yīng)用于巖土多孔材料的研究。例如，D．Coelho等［13］模擬了球體、橢球體、圓柱體和扁橢球體的填筑，分析了各自的不同幾何尺寸比例與孔隙率、比表面積、排列方向性的關(guān)系和相關(guān)函數(shù)，得到的認(rèn)識包括：填筑過程中顆粒的位移和旋轉(zhuǎn)程度大就會使得孔隙率更大；幾何尺寸比例為1的條件下，填筑后得到的孔隙率最小；扁橢球體的定向排列對結(jié)果有明顯影響等。這些認(rèn)識都有一定的普遍意義和參考價(jià)值。但是由于其研究對象限定為幾種規(guī)則形狀的顆粒，研究手段限于數(shù)值模擬，而且填筑體的模擬以在有限空間內(nèi)獲得最大填筑密度為目標(biāo)，這些與巖土工程中關(guān)心的問題和背景都有很大區(qū)別：顆粒形狀實(shí)際上更復(fù)雜，顆粒大小不均一，土體的密度取決于土的擊實(shí)或固結(jié)程度等，所以，他們關(guān)于熱傳導(dǎo)性、滲透性和彌散系數(shù)與顆粒形狀無關(guān)的認(rèn)識就與客觀規(guī)律相差甚遠(yuǎn)。

R．Williams等［14］提出通過對顆粒位圖的數(shù)字化來表達(dá)顆粒的任意形狀，認(rèn)為填筑體的模擬結(jié)果可以用于預(yù)測微結(jié)構(gòu)的物理性質(zhì)，如熱傳導(dǎo)率、氣體和液體的滲透性等。M．Stewart等［4］用格子Boltz-mann方法模擬扁橢球體顆粒填筑體中的水流運(yùn)動(dòng)，企圖研究孔隙尺度上由于顆粒排列引起的各向異性。其結(jié)果表明，顆粒三維尺寸之間的比值對于顆粒的排列結(jié)果及滲透張量的各向異性程度都有顯著影響，雖然同樣僅限于對規(guī)則形狀及均一顆粒填筑體的數(shù)值模擬，但是取得的認(rèn)識和結(jié)論比起 D．Coelho等的研究［13］顯得更合理。

土體顆粒之間的接觸關(guān)系既取決于顆粒形狀、顆粒大小和級配，更取決于密實(shí)程度，所以實(shí)際問題中土顆粒填筑的模擬是非常復(fù)雜的，填筑體中固相物質(zhì)之間孔隙系統(tǒng)的模擬則更加復(fù)雜。短期內(nèi)，除了開展較規(guī)則顆粒、簡單條件下的顆粒填筑和孔隙系統(tǒng)模擬方法研究外，采用可視化方法直接研究孔隙系統(tǒng)更具有可操作性。

3．3 孔隙系統(tǒng)的可視化研究

研究顆粒形狀對孔隙系統(tǒng)的影響，借助微細(xì)觀試驗(yàn)研究手段獲取多孔介質(zhì)的孔隙系統(tǒng)圖像不僅是更直觀有效的途徑，還可以驗(yàn)證數(shù)值模擬結(jié)果。

土的微細(xì)觀研究是一直是土質(zhì)學(xué)和土力學(xué)研究的重要方面之一［15］。施斌等對掃描電鏡分析的有關(guān)技術(shù)開展了大量研究工作［16］。一些學(xué)者采用掃描電鏡、自x－光衍射技術(shù)等研究粘性土的成果在揭示土的物理力學(xué)特性及其內(nèi)在機(jī)理方面起到了很大的作用［17－19］。然而，掃描電鏡分析方法也有很多技術(shù)上的限制［16］。周健等［20－22］采用立體顯微鏡數(shù)字采集系統(tǒng)結(jié)合攝像／照相機(jī)觀察和記錄飽和砂土液化和砂土管涌試驗(yàn)中的細(xì)觀現(xiàn)象，并采用顆粒流方法進(jìn)行數(shù)值模擬，得到了一些有益的認(rèn)識，也說明了細(xì)觀研究的意義，但是這種方法只能采集模型表面的信息。采用計(jì)算機(jī)控制X射線斷層分析技術(shù)（computerized tomography，簡稱CT）從細(xì)觀上研究土的結(jié)構(gòu)性，尤其是研究粗粒土的組構(gòu)問題，不僅可以做到無損檢測，還可以實(shí)現(xiàn)多層面掃描。

程展林等［23］認(rèn)為，當(dāng)前土體的結(jié)構(gòu)性研究處于停滯狀態(tài)，其根本的原因在于試驗(yàn)過程中對土的結(jié)構(gòu)變化難以動(dòng)態(tài)定量觀測。葛修潤院士等［24］詳細(xì)介紹了CT技術(shù)用于巖土問題細(xì)觀研究的原理和優(yōu)勢。A．Kilfeather等［25］利用掃描電鏡結(jié)合μCT技術(shù)研究了冰磧土的孔隙大小、形狀和連通性，以及它們與變形過程的關(guān)系，得出的主要結(jié)論包括：由顆粒分析結(jié)果不能可靠地預(yù)測孔隙率，與顆粒填筑理論研究結(jié)果相比，碎屑物含量高的冰磧土具有更高的孔隙性；對于大多數(shù)冰磧土來說，微觀結(jié)構(gòu)是控制孔隙性的極重要因素。吳恒等［26］采用CT技術(shù)研究溶蝕作用下土體孔隙的變化，是CT技術(shù)用于滲流問題和水巖相互作用研究的有益嘗試。

長江科學(xué)院新近購置了高分辨率的西門子40層CT機(jī)，正在探討用以研究粗粒土組構(gòu)、膨脹土的損傷機(jī)制、土的滲透變形演化過程、土與結(jié)構(gòu)物結(jié)合面的相互作用機(jī)理以及巖石損傷破壞機(jī)理、疲勞損傷斷裂機(jī)理等。考慮到粗粒土的顆粒尺度相對較大，結(jié)構(gòu)特征相對簡單，主要體現(xiàn)在顆粒本身及顆粒間幾何排列方面，程展林等［23］建議土體結(jié)構(gòu)力學(xué)研究從粗粒土研究入手，利用CT技術(shù)研究粗粒土試樣受力變形過程中內(nèi)部結(jié)構(gòu)的動(dòng)態(tài)變化。為此，長江科學(xué)院已經(jīng)開發(fā)了粗粒土CT三軸試驗(yàn)儀和計(jì)算機(jī)圖像測量系統(tǒng)，初步解決了粗粒土組構(gòu)的量化問題［27］。已經(jīng)取得的粗粒土CT圖像清晰，能夠反映各顆粒的位置和形態(tài)，也為孔隙系統(tǒng)及滲流問題的研究提供了重要手段。

3．4 顆粒形狀對多孔介質(zhì)滲流規(guī)律的影響研究

多孔介質(zhì)本身就是應(yīng)宏觀研究需要而提出的概念。J．Bear甚至指出，盡管固體表面對流動(dòng)起邊界作用，但任何試圖以精確的方式對其幾何形狀進(jìn)行描述的想法都是徒勞的，研究孔隙中的流體并描述有關(guān)的運(yùn)動(dòng)、質(zhì)量運(yùn)移等現(xiàn)象也會遇到同樣困難［1］。然而，只要定性地分析一下非球形顆粒、尤其是多面體顆粒與球形顆粒對水流質(zhì)點(diǎn)的作用存在的差異，就可以理解在某些條件下不應(yīng)忽視顆粒形狀的影響。

當(dāng)水流質(zhì)點(diǎn)與顆粒表面接觸時(shí)，球形顆粒多孔介質(zhì)孔隙中水流質(zhì)點(diǎn)主要是沿切向運(yùn)動(dòng)；而在多面體顆粒多孔介質(zhì)孔隙系統(tǒng)內(nèi)，水流質(zhì)點(diǎn)可能一時(shí)平行于顆粒表面運(yùn)動(dòng)，從而損失能量很小，而另一時(shí)又正面或斜向撞向固體表面后沿折線運(yùn)動(dòng)，從而損失更大的能量。如果多面體顆粒表面粗糙，或有意制造為凹面，那么水流質(zhì)點(diǎn)在界面作用下將主要沿復(fù)雜折線運(yùn)動(dòng)，使得能量損失甚于球形顆粒多孔介質(zhì)滲流。

至于顆粒形狀對多孔介質(zhì)宏觀滲流特性的影響還是可以采用已有的試驗(yàn)方法。在保證合理的試驗(yàn)精度條件下，通過試驗(yàn)對比研究相同級配、不同形狀顆粒多孔介質(zhì)的滲流規(guī)律，包括滲透性、滲透變形特性、反濾保護(hù)作用效果、過濾除污作用等，可以揭示顆粒形狀對滲流規(guī)律的影響，豐富多孔介質(zhì)滲流規(guī)律的認(rèn)識。

提高濾層的有效性是水處理技術(shù)研究的重要任務(wù)之一，其途徑很多，包括濾料本身的材料研制，濾料顆粒粒徑和級配設(shè)計(jì)，以及濾層厚度設(shè)計(jì)等。田家宇等通過對比試驗(yàn)說明，粒徑較小的濾料形成的濾柱，由于具有更大的表面積，因而具有更高的吸附容量和微生物生長空間，所以對污水具有更優(yōu)越的除污能力［28］。湯利華分析了球形顆粒濾料粒徑對濾層孔隙率、孔隙尺寸、濾料比表面積的影響，認(rèn)為顆粒粒徑對孔隙率沒有影響，與濾料當(dāng)量孔隙直徑正相關(guān)，與濾料的總比表面積負(fù)相關(guān)，而且粒徑越小，水中懸浮物與濾料接觸次數(shù)越多，所以較小粒徑的濾料對保證濾后水質(zhì)有利［29］，但是隨著顆粒粒徑的減小，濾層的滲透系數(shù)降低，會降低清潔水的生產(chǎn)效率。這些研究都是以球形顆粒為研究對象。多面體顆粒多孔介質(zhì)不僅具有上面的滲流規(guī)律，而且顆粒具有更大的比表面積，利用多面體多孔介質(zhì)的這些特性開發(fā)濾料，可以開辟提高濾層效率的新途徑。

除了顆粒本身的力學(xué)性質(zhì)外，碎石料與砂礫石料的顯著差別就是顆粒形狀，其對土體的壓實(shí)性質(zhì)和滲透性的影響已經(jīng)引起了關(guān)注。劉杰通過相同級配砂礫石料和人工碎石料的對比試驗(yàn)和成果分析［3］，認(rèn)為碎石料干密度和滲透系數(shù)的變化范圍都大于砂礫石料，碎石料的滲透性不象一般土料一樣明顯地受粒徑小于5 mm顆粒的含量控制，用砂礫石料的判別準(zhǔn)則來判別碎石料的滲透破壞形式也容易與滲透變形試驗(yàn)結(jié)果發(fā)生偏差。作者在研究水布埡大壩填料的滲透變形特性時(shí)，也注意到同一級配的試樣在相同密度下的滲透系數(shù)甚至?xí)辛考壣系淖兓@應(yīng)該與不同試驗(yàn)的試樣之間顆粒粒徑組成不完全一致有關(guān)，也與顆粒形狀的差異脫不了干系［30］。根據(jù)已有的反濾設(shè)計(jì)準(zhǔn)則，爆破碎石料可能由于內(nèi)部結(jié)構(gòu)欠穩(wěn)定而不滿足反濾料的要求，但是水布埡面板壩墊層料與過渡料的反濾試驗(yàn)表明，雖然在低比降條件下就有細(xì)料流失，但是損失量很小，在實(shí)際工程中可能經(jīng)歷的比降條件下，過渡料仍然可以對墊層料起到良好的保護(hù)作用，大壩滲透穩(wěn)定性可以滿足要求，這就可能與爆破碎石料的孔隙特征和滲透變形特性有關(guān)［31］。進(jìn)一步研究碎石料與一般土料滲流特性的差異，將有助于完善現(xiàn)有的規(guī)律性認(rèn)識，甚至可能促進(jìn)土體滲透變形特征和反濾效果判斷準(zhǔn)則的改進(jìn)。

4 結(jié) 語

顆粒形狀對滲流規(guī)律的影響應(yīng)該是多孔介質(zhì)滲流研究的一個(gè)重要方面，而且已經(jīng)開展的工作還很不夠。本文初步探討了研究思路和方法，目前相關(guān)研究條件已經(jīng)初步具備，可以先由簡單問題開始，逐步向復(fù)雜問題和工程應(yīng)用推進(jìn)。作者正在開展規(guī)則多面體多孔介質(zhì)孔隙特征和滲流規(guī)律的研究，相關(guān)成果見后續(xù)報(bào)道。

［1］ J．貝爾．多孔介質(zhì)流體動(dòng)力學(xué)［M］．李競生，陳崇希，譯．北京：中國建筑工業(yè)出版社．1983．（Bear J．Dy-namics of Fluids in Porous Media［M］．Translated by LI Jing-sheng，CHEN Chong-xi．Beijing：China Architecture＆Building Press，1983．（in Chinese））

［2］ 鄭穎人．巖土塑性力學(xué)的新進(jìn)展——廣義塑性力學(xué)［J］．巖土工程學(xué)報(bào)，2003，（1）：1－10．（ZHENG Ying-ren．New Development of Geotechnical Plastic Mechanics Generalized Plastic Mechanics［J］．Chinese Journal of Geotechnical Engineering，2003，（1）：1－10．（in Chi-nese））

［3］ 劉 杰．混凝土面板壩碎石墊層料最佳級配試驗(yàn)研究［J］．水利水運(yùn)工程學(xué)報(bào)，2001，（4）：1－7．（LIU Jie．Experimental Study of Optimum Gradation of Crushed Stone Cushion for Concrete-face Rockfill Dam［J］．Hydro-science and Engineering，2001，（4）：1－7．（in Chinese））

［4］ STEWARTM，WARD A，RECTOR D．A Study of Pore Geometry Effects on Anisotropy in Hydraulic Permeability Using the Lattice-Boltzmann Method［J］．Advances in Water Resources，2006，（29）：1328－1340．

［5］ 孔祥言，盧德唐．滲流力學(xué)的理論應(yīng)用及其前沿研究［J］．中國科學(xué)技術(shù)大學(xué)學(xué)報(bào)，2007，（10）：1262－1266．（KONG Xiang-yan，LU De-tang．Application and Development of Fluid in Multipore Media Flowing［J］．Journal of Science and Technology of China，2007，（10）：1262－1266．（in Chinese））

［6］ H．達(dá)西．水穿越砂流動(dòng)定律的確定［C］∥地下水水力學(xué)的發(fā)展．張宏仁，編譯．北京：地質(zhì)出版社，1992．（DARCY H．Les Fontaines Publiques de la Ville de Di-jon［C］∥Development of Groundwater Hydraulics．Translated by ZHANG Hong-ren．Beijing：Geological Publishing House，1992．（in Chinese））

［7］ 范淑果，郝宏偉，杜 皓．選擇圓度誤差評定方法的體會［J］．計(jì)量與測試技術(shù)，2009，（10）：28－30．（FAN Shu-guo，HAO Hong-wei，DU Hao．Impressions on Choosing Evaluating Method of Roundness Error［J］．Metrology and Measurement Technique，2009，（10）：28－30．（in Chinese））

［8］ 卓興仁，劉煥文，王永波，等．圓度誤差4種評定方法評定結(jié)果關(guān)系的研究［J］．北京工業(yè)大學(xué)學(xué)報(bào)，1996，（2）：111－114．（ZHUO Xing-ren，LIU Huan-wen，WANG Yong-bo，etal．The Relation of the Different Ver-ifying Results of Circularity Deviation through Four Differ-entWays［J］．Journal of Beijing Polytechnic University，1996，（2）：111－114．（in Chinese））

［9］ 史旦達(dá)，周 健，劉文白，等．循環(huán)荷載作用下砂土液化特性的非圓顆粒數(shù)值模擬［J］．水利學(xué)報(bào)，2009，（9）：1074－1082．（SHIDan-da，ZHOU Jian，LIU Wen-bai，et al．Numerical Simulation of Liquefaction of Sands with Non-circular Particles under Cyclic Loading［J］．Journal of Hydraulic Engineering，2009，（9）：1074－1082．（in Chinese））

［10］JIA X，WILLIAMSR．A Packing Algorithm for Particles of Arbitrary Shapes［J］．Powder Technology，2001，（120）：175－186．

［11］張?jiān)戚x，譚慶昌，王樹山，等．一種基于數(shù)字圖像的圓度測量方法［J］．工程與試驗(yàn)，2008，（3）：10－13．（ZHANG Yun-hui，TAN Qing-chang，WANG Shu-shan，et al．AMethod of Measuring Roundness Based on Digital Image Process［J］．Test Technology and Testing Ma-chine，2008，（3）：10－13．（in Chinese））

［12］何 琨，黃文奇．求解長方體Packing問題的高效算法［J］．計(jì)算機(jī)研究與發(fā)展，2008，45（增刊1）：7－10．

（HE Kun，HUANG Wen-qi．An Efficient Algorithm for Solving Cuboids Packing Problem［J］．Journal of Com-puter Research and Development，2008，45（Sup．l）：7－10．（in Chinese））

［13］COELHO D，THOVERT JF，ADLER PM．Geometrical and Transport Properties of Random Packings of Spheres and Aspherical Particles［J］．Physical Review，1997，55（2）：1959－1978．

［14］WILLIAMS R，JIA X．A New Method for Prediction of Bulk Particle Packing Behavior for Arbitrary-shaped Parti-cles in Containers of Any Shape［J］．Particulate Science and Technology，2003，21：195－205．

［15］施 斌．粘性土微觀結(jié)構(gòu)研究回顧與展望［J］．工程地質(zhì)學(xué)報(bào)，1996，（1）：39－44．（SHI Bin．Review and Prospect on the Microstructure of Clayey Soil［J］．Jour-nal of Engineering Geology，1996，（1）：39－44．（in Chi-nese））

［16］施 斌，姜洪濤．粘性土的微觀的結(jié)構(gòu)分析技術(shù)研究［J］．巖石力學(xué)與工程學(xué)報(bào)，2001，（6）：864－870．（SHI Bin，JIANG Hong-tao．Research on the Analysis Tech-niques for Clay Soil Microstructure［J］．Chinese Journal Rock Mechanics and Engineering，2001，（6）：864－870．（in Chinese））

［17］蔣明鏡，沈珠江，邢素英，等．結(jié)構(gòu)性粘土研究綜述［J］．水利水電科技進(jìn)展，1999，（1）：26－30．（JIANG Ming-jing，SHEN Zhu-jiang，XING Su-ying，et al．Re-view of Study on Structural Clay［J］．Advances in Sci-ence and Technology ofWater Resources，1999，（1）：26－30．（in Chinese））

［18］孔令偉，郭愛國，呂海波．典型紅粘土的基本特性與微觀結(jié)構(gòu)特征［J］．巖石力學(xué)與工程學(xué)報(bào)，2001，（增刊1）：973－977．（KONG Ling-wei，GUO Ai-guo，LU Hai-bo，et al．Fundamental Properties and Microstructure Characteristics of Typical Red Clay［J］．Chinese Journal Rock Mechanics and Engineering，2001，（Sup．1）：973－977．（in Chinese））

［19］姚海林，劉少軍，程昌炳．一種天然膠結(jié)土粘聚力的微觀本質(zhì)［J］．巖石力學(xué)與工程學(xué)報(bào)，2001，（6）：871－874．（YAO Hai-lin，LIU Shao-jun，CHENG Chang-bing．Scopic Essence of Cohesion of A Natural Cemented Soil［J］．Chinese Journal Rock Mechanics and Engineering，2001，（6）：871－874．（in Chinese））

［20］周 健，楊永香，劉 洋．飽和砂土液化過程中細(xì)觀組構(gòu)的模型試驗(yàn)研究［J］．同濟(jì)大學(xué)學(xué)報(bào)（自然科學(xué)版），2009，（4）：466－470．（ZHOU Jian，YANG Yong-xiang，LIU Yang．Model Testing of Meso-fabric of Satu-rated Sand Liquefaction［J］．Journal of Tongji University（Natural Science），2009，（4）：466－470．（in Chi-nese））

［21］周 健，姚志雄，白彥峰，等．砂土管涌的細(xì)觀機(jī)理研究［J］．同濟(jì)大學(xué)學(xué)報(bào)（自然科學(xué)版），2008，（6）：733－737．（ZHOU Jian，YAO Zhi-xiong，BAIYan-feng，et al．Meso-mechanical Study on Piping in Sandy Soils［J］．Journal of Tongji University（Natural Science），2008，（6）：733－737．（in Chinese））

［22］周 健，姚志雄，張 剛．砂土滲流過程的細(xì)觀數(shù)值模擬［J］．巖土工程學(xué)報(bào)，2007，（7）：977－981．（ZHOU Jian，YAO Zhi-xiong，ZHANG Gang．Mesomechanical Simulation of Seepage Flow in Sandy Soil［J］．Chinese Journal of Geotechnical Engineering，2007，（7）：977－981．（in Chinese））

［23］程展林，丁紅順，吳良平．粗粒土試驗(yàn)研究［J］．巖土工程學(xué)報(bào)，2007，（8）：1151－1158．（CHENG Zhan-lin，DING Hong-shun，WU Liang-ping．Experimental Study on Mechanical Behavior of Granular Material［J］．Chi-nese Journal of Geotechnical Engineering，2007，（8）：1151－1158．（in Chinese））

［24］葛修潤，任建喜，蒲毅彬，等．巖土損傷力學(xué)宏細(xì)觀試驗(yàn)研究［M］．北京：科學(xué)出版社，2004．

［25］KILFEATHERA A，J．Van Der Meer．Pore Size，Shape and Connectivity in Tills and Their Relationship to De-formation Processes［J］．Quaternary Science Reviews，2008，（27）：250－266．

［26］吳 恒，張信貴，易念平．水土作用與土體細(xì)觀結(jié)構(gòu)研究［J］．巖石力學(xué)與工程學(xué)報(bào)，2000，（2）：199－204．（WU Heng，ZHANG Xin-gui，YI Nian-ping．Research on Water-soil Interaction and Macrostructure of Soil Mass［J］．Chinese Journal Rock Mechanics and Engineering，2000，（2）：199－204．（in Chinese））

［27］程展林，吳良平，丁紅順．粗粒土組構(gòu)之顆粒運(yùn)動(dòng)研究［J］．巖土力學(xué)，2007，28（增刊）：29－33．（CHENG Zhan-lin，WU Liang-ping，DING Hong-shun．Research on Movement of Particle of Fabric of Granular［J］．Rock and Soil Mechanics，2007，28（Sup．）：1151－1158．（in Chinese））

［28］田家宇，陳偉雄，王 威，等．濾料粒徑對生物活性濾池除污染效能的影響［J］．工業(yè)水處理，2007，（11）：15－18．（TIAN Jia-yu，CHENWei-xiong，WANGWei，et al．Effects of Medium Sizes on Pollutant Removal Rate of Biological Activated Filters［J］．Industrial Water Treatment，2007，（11）：15－18．（in Chinese））

［29］湯利華．濾料粒徑對過濾的影響［J］．安徽建筑工業(yè)學(xué)院學(xué)報(bào)（自然科學(xué)版），2000，（3）：45－48．（TANG Li-hua．Effects of Particle Size of Filter Media on Filtration［J］．Journal of Anhui Institute of Architecture（Natural Science），2000，（3）：45－48．（in Chinese））

［30］張家發(fā)，定培中，張 偉，等．水布埡面板堆石壩墊層料滲透與滲透變形特性試驗(yàn)研究［J］．巖土力學(xué)，2009，（10）：3145－3150．（ZHANG Jia-fa，DING Pei-zhong，ZHANGWei，et al．Studies of Permeability and Seepage Deformation Characteristics of Cushion Material for Shuibuya Concrete Faced Rockfill Dam［J］．Rock and Soil Mechanics，2009，（10）：3145－3150．（in Chi-nese））

［31］定培中，張 偉，張家發(fā)，等．水布埡混凝土面板堆石壩壩體填料滲透穩(wěn)定性試驗(yàn)研究［R］．長江科學(xué)院，2005．（DING Pei-zhong，ZHANG Wei，ZHANG Jia-fa，et al．Experimental Studies on the Seepage Stability of Rockfill Materials for Shuibuya Concrete Face Rockfill Dam［R］．Yangtze River Scientific Research Institute， 2005．（in Chinese） ）

（編輯：曾小漢）

Influence of Grain Shape on Characteristics of Pores and Seepage in Porous M edia

ZHANG Jia-fa1，JIAO Jiu-jiu2
（1．Key Laboratory of Geotechnical Mechanics and Engineering of the Ministry ofWater Resources，Yangtze River Scientific Research Institute，Wuhan 430010，China；2．Department of Earth Sciences，The University of Hong Kong，Hong Kong，China）

This paper discusses the possible influence of grain shapes on the pore system of porousmedia and its seepage rules，aswell as its research significance．By reviewing research results in related areas，discussionswere made on themeasurementmethod of grain shape，simulationmethod of grain shape and grain packing，visualization method of pore system，and experimental comparison method of seepage rules．The research result helps demon-strate the prospect of developing filtermaterials for water treatment and applying aggregate for dam construction．

grain shape；porousmedia；seepage；test；simulation

TU411．4

A

1001－5485（2011）03－0039－06

2010-10-12

國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃（2007CB714106）；國家“十一五”科技支撐計(jì)劃（2008BAB29B02）；中央級公益性科研院所基本科研業(yè)務(wù)費(fèi)（YWF0909）

張家發(fā)（1960-），男，安徽安慶人，教授級高級工程師，主要從事巖土工程和水工滲流研究，（電話）027-82820029（電子信箱）seep＠public．wh．hb．cn。