濟南淺層地基微觀結(jié)構(gòu)及分形特征研究

王海林 路林海 王國富 孫紅

摘要：天然沉積土的宏觀力學(xué)特性與微觀結(jié)構(gòu)緊密相關(guān)，研究土體微觀結(jié)構(gòu)對深入揭示土的宏觀力學(xué)特性具有重要意義。首先研究濟南地區(qū)深度在50 m以內(nèi)黃河沖積地層天然沉積土的物理力學(xué)參數(shù)分布特征，再采用場發(fā)射掃描電鏡試驗，分析各土層不同方向土樣的微觀結(jié)構(gòu)特征，最后應(yīng)用分形理論，研究各層土分形維數(shù)的演化規(guī)律及其與物理特性參數(shù)間的相關(guān)性。結(jié)果表明：第⑦層黃土是以粗粉粒為主體骨架的蜂窩狀大孔隙結(jié)構(gòu)，第⑧層粉質(zhì)黏土是碎屑粉粒為骨架的蜂窩結(jié)構(gòu)，第⑨、⑩和⑩層粉質(zhì)黏土是以黏粒集合體為基本單元的薄層狀絮凝結(jié)構(gòu)。不同土層土體土粒和孔隙的分形行為明顯，隨著地基深度增大，土粒和孔隙的分布分維數(shù)均逐漸加大，形態(tài)分維數(shù)均逐漸減小，土粒的粒徑變化比孔隙的孔徑變化明顯，而孔隙的形態(tài)變化比土粒復(fù)雜。土體的分維數(shù)與含水率、孔隙比、液塑限、黏聚力和內(nèi)摩擦角之間存在較強的相關(guān)性。

關(guān)鍵詞：天然沉積土;微觀結(jié)構(gòu);分形維數(shù);物理特性;相關(guān)性

中圖分類號：TU411

文獻(xiàn)標(biāo)志碼：A

doi：10.3969/j .issn. 1000- 1379.2019.01. 026

濟南市地處黃河沖積平原和山前平原的交匯處，具有典型的黃河沖積地層地質(zhì)條件，該地區(qū)地下結(jié)構(gòu)所處的土層一般以黏土、粉土和粉質(zhì)黏土為主[1]。地層所在區(qū)域土體強度較為軟弱，土體穩(wěn)定性差，部分地區(qū)地下水位高，施工風(fēng)險大，易出現(xiàn)工程事故，因此有必要對濟南地區(qū)土體的工程性質(zhì)進行研究。

土體的工程性質(zhì)是結(jié)構(gòu)單元體性質(zhì)的綜合表現(xiàn)，結(jié)構(gòu)單元體的性質(zhì)受到土粒集合體甚至更小的單粒礦物性質(zhì)的影響[2]。大量研究已經(jīng)證實，土體的一些宏觀特性與微觀結(jié)構(gòu)特征緊密相關(guān)，微觀結(jié)構(gòu)的表現(xiàn)會對土體的宏觀性質(zhì)產(chǎn)生一定影響[3]。土體強度[4-5]、熱傳導(dǎo)特性[6]、地基承載力[7]等方面的研究表明，土體微觀結(jié)構(gòu)的特征及分布是影響土體相關(guān)特性的重要因素。因此，除了研究土體的宏觀特性，還勢必要對土體的微觀結(jié)構(gòu)進行分析。

土體微觀結(jié)構(gòu)的研究目前主要側(cè)重于顆粒和孔隙兩個方面，通過SEM和MIP等技術(shù)，許多學(xué)者[8-10]已經(jīng)對顆粒的表面積、體積和孔隙率、分布特征進行了研究，且有了深入的認(rèn)識。由于土體微觀結(jié)構(gòu)存在自相似性[2]，因此可以通過分形理論來對土體微觀結(jié)構(gòu)進行探究，從而揭示其內(nèi)在特性。目前，結(jié)合分形理論，在土體及微觀結(jié)構(gòu)方面已開展了大量的研究，并在應(yīng)力路徑影響[10]、膨脹土分類[11]、孔隙分維變化與土體演化[12]、土體宏觀強度的相關(guān)性[13]、凍融循環(huán)[14]、非飽和土有效應(yīng)力[15]、滲透系數(shù)、土水特征曲線[16]等方面取得了眾多成果，反映了分形理論在土體微觀結(jié)構(gòu)研究方面的巨大作用。

為分析與評價濟南地區(qū)不同土體的工程特征，筆者對濟南地區(qū)淺層地基進行一系列宏微觀試驗，研究濟南地區(qū)不同土層土體微觀結(jié)構(gòu)特征，探究土體孔隙和土粒的分形演化機理。

1 試驗方案

土樣取自濟南某工地，鉆孔后利用薄壁取土器取得第⑦（黃土）、⑧～11（粉質(zhì)黏土）層土樣，深度為2.1～50.0 m.其中：⑦～⑨層屬于第四系全新統(tǒng)沖洪積層，⑩、11層屬于第四系上更新統(tǒng)沖洪積層。

根據(jù)土工試驗方法標(biāo)準(zhǔn)[17]進行試驗，獲得50 m以內(nèi)土層的基本物理參數(shù)，然后從原狀土樣的橫斷面和縱斷面上選取試樣進行場發(fā)射掃描電鏡試驗，研究各土層不同方向土樣的微觀結(jié)構(gòu)。微觀掃描試驗采用TESCAN拉曼一電鏡一體化場發(fā)射掃描電鏡。采用液氮冷凍升華干燥法[8]制取電鏡掃描所用樣品，然后用小刀將干燥后的土樣輕輕切開，獲得用于觀察的新鮮斷面，用導(dǎo)電膠固定后進行噴金處理，最后進行電鏡掃描試驗。

為了獲得能夠?qū)ν馏w微觀結(jié)構(gòu)進行定量分析的黑白二值化圖像，借助IPP圖像處理軟件對掃描電鏡試驗獲得的灰度圖片進行灰度閾值處理，采用目視分割法[18]進行多次試驗，得出電鏡圖片的最佳分割閾值，而后對土樣的二值化圖片進行定量分析。

2 宏觀特性分析

淺層地基的土層剖面及基本物理參數(shù)見圖1.其中：ρ為密度;e為孔隙比;w為含水率;w_L為液限;w為塑限;K為滲透系數(shù);K為垂直滲透系數(shù);c為黏聚力;φ為內(nèi)摩擦角。沿著深度方向，土層呈現(xiàn)多個明顯的分層，隨著深度的增大，孔隙比、含水率、液塑限和黏聚力均呈現(xiàn)出一定程度的增大：對于滲透系數(shù)，隨著深度加深，橫斷面和縱斷面的滲透系數(shù)均呈現(xiàn)出先增大而后逐漸減小的趨勢。其中，第⑨層土的滲透系數(shù)最大，原因是第⑨層土含有較多粒狀的鈣質(zhì)結(jié)核，加大了土體的滲透能力。對比不同斷面的滲透系數(shù)可以看出，縱斷面的滲透系數(shù)均略大于橫斷面的滲透系數(shù)，說明土體的滲透系數(shù)具有各向異性。

3 微觀結(jié)構(gòu)特征分析

以橫斷面土體為例，定性分析不同土層土體的微觀結(jié)構(gòu)特征。圖2為不同土層的土樣橫斷面在放大倍數(shù)為2 000倍和10 000倍時掃描電鏡圖片。

第⑦層土（黃土）呈現(xiàn)出團聚體狀，土粒表面上間或附有黏粒狀礦物，土粒之間的接觸以點一點接觸為主，相互間的接觸面積較?。嚎紫兑约芸諣顟B(tài)的孔隙分布為主，孔隙區(qū)域明確，發(fā)育完全，層次清晰，偶有鑲嵌狀態(tài)的孔隙。總的來說，第⑦層黃土是以粗粉粒為主體骨架的蜂窩狀大孔隙結(jié)構(gòu)。第⑧、⑨、⑩和11層土均為粉質(zhì)黏土。第⑧層土土粒較小，非均質(zhì)排列的粉粒和片狀的黏性礦物（圖2（d））聚集或黏附在土粒表面或者接觸點上，土粒之間的接觸關(guān)系為鑲嵌與分散交錯分布，孔隙類型以開放性孔隙為主，與第⑦層土對比，大孔隙數(shù)量開始減少，小孔隙數(shù)量明顯增加，孔隙開始逐漸被填充，土體的整體微觀結(jié)構(gòu)類型屬于碎屑粉粒骨架狀蜂窩結(jié)構(gòu)。第⑨層土除了個別土粒外，土粒整體上比較模糊，相互間接觸面上的碎屑和片狀黏性礦物較多，土粒呈現(xiàn)浸潤狀態(tài)，孔隙大范圍減少。觀察第⑩層土可以發(fā)現(xiàn)，土體結(jié)構(gòu)趨于致密，孔隙主要以鑲嵌孔隙和粒內(nèi)孔隙的形式存在，土粒以致密狀態(tài)的片狀和書卷狀礦物為主，原因是地質(zhì)歷史演變的過程中，在長期的地質(zhì)沉積和外荷載作用下，土粒之間相互作用形成團粒，同時在逐步加壓的狀態(tài)下，土體長期固結(jié)，進而形成扁平狀的聚集體，土體土粒間的接觸以面一面、面一邊形式接觸為主，土體結(jié)構(gòu)類型屬于以黏粒集合體為基本單元的薄層狀絮凝結(jié)構(gòu)。這種結(jié)構(gòu)在第11層土中尤其明顯。這種土一般具有較高的強度。

總的說來，隨著深度的增大，土粒間的大孔隙數(shù)量減少，小孔隙數(shù)量增多：黃土骨架間接觸形式基本以點一點接觸為主，而粉質(zhì)黏土土粒間接觸形式基本以面一面、面一邊接觸為主，深部土體趨于致密。

4 微觀結(jié)構(gòu)的分形特征

4.1 分形維數(shù)的計算

土體的微觀結(jié)構(gòu)整體上規(guī)則性差，復(fù)雜性和無序性高，但其土粒與孔隙均具有一定的自相似性，因此可以采用分形理論對其進行分析。這里選用分布分形模型和形態(tài)分形模型來研究不同土層不同斷面土體微觀結(jié)構(gòu)中孔隙和土粒的分形特征。

分布分形模型[19-20]是將測量圖像中所研究物體的直徑和數(shù)量聯(lián)系起來：

4.2 土體孔隙的分形特征

圖3為各土層的孔隙孔徑分布分維數(shù)圖。所有土層的數(shù)據(jù)都基本成線性關(guān)系，相關(guān)系數(shù)為0.93～ 0.97。雖然試樣土層不同，但其孑L隙的孑L徑分布仍具有明顯的分形特征。

繪制不同斷面孔隙的孔徑分布分維數(shù)和深度的關(guān)系，見圖4。可見，孔徑分布分維數(shù)集中在1. 21～1.45范圍內(nèi)，波動較小，平均值約為1.36。隨著深度的增大，土體的橫斷面和縱斷面孔徑分布分維數(shù)均增大。這主要是因為隨著土層深度的加深，大孔隙數(shù)量減少，而小孔隙數(shù)量逐漸增多，孔隙直徑逐漸趨于細(xì)化，所以孔徑分布分維數(shù)逐漸增大。

根據(jù)形態(tài)分維數(shù)的定義和計算方法，繪制不同土層不同斷面孔隙的周長P和面積A在雙對數(shù)坐標(biāo)系中的散點分布圖（見圖5）。以⑧層粉質(zhì)黏土橫斷面為例，很明顯圖5中土體孔隙lg P-lg A的線性關(guān)系非常明顯，確定系數(shù)R在0.96以上，說明土體孔隙的形態(tài)分布具有較強的分形特征。

繪制不同土層土的孔隙形態(tài)分維數(shù)與深度的關(guān)系，見圖6。可以發(fā)現(xiàn)，整個淺層地基的孔隙形態(tài)分維數(shù)為1.34～1.52，波動范圍比較小，平均值為1.43。隨著土層深度的增大，形態(tài)分維數(shù)D逐漸減小，說明孔隙的復(fù)雜程度在降低。原因是在漫長的地質(zhì)歷史過程中，隨著土體固結(jié)壓密，土粒之間的接觸程度加深，排列和分布趨于緊密，土體中的孔隙受土粒不斷聚集和定向排列的影響，逐漸變化為扁圓形的狀態(tài)，導(dǎo)致形態(tài)的復(fù)雜程度減弱。

從微觀角度看，由于土層深度增大，孔隙細(xì)化，孔隙的比表面積加大，土粒表面容易形成一定厚度的結(jié)合水膜，同時土粒排列緊密，土粒間結(jié)合水聯(lián)結(jié)程度增強，減弱了自由水的運動，因此引起宏觀上不同土層土體的含水率、液塑限值隨著深度增大而逐漸增大。由圖6還可以發(fā)現(xiàn)，橫斷面孔隙的孔徑分布分維數(shù)和形態(tài)分維數(shù)與縱斷面的數(shù)值間存在一定差異，原因是土體固結(jié)過程中土粒定向排列，土體豎直方向（橫斷面）受到的應(yīng)力一般大于水平方向（縱斷面）受到的應(yīng)力，導(dǎo)致橫、縱斷面上孔隙分布不同，進而引起孔隙分形維數(shù)不同。這反映出土層在水平和垂直方向上的固結(jié)程度存在差異，土體呈現(xiàn)出各向異性的特點，進一步說明了土體的滲透系數(shù)具有各向異性。

4.3 土粒的分形特征

不同深度土粒粒徑分布分維數(shù)見圖7。整個淺層地基的土粒粒徑分布分維數(shù)為1.62～1.86.波動范圍比較小，平均值為1.74。第⑦層黃土由于以大土粒為主，土粒粒徑間差異小，土粒粒徑整體趨于均勻，分選性較好，因此其粒徑分布分維數(shù)相對較小，為1.63左右;而其下的⑧、⑨、⑩和11層粉質(zhì)黏土顆粒大小存在差異，且土粒間的孔隙被細(xì)土粒填充，整體上土粒粒徑分化顯著，粒徑分布分維數(shù)隨著深度加大逐漸增大，第11層土粒徑分布分維數(shù)最大。同時，橫斷面的粒徑分布分維數(shù)均略大于縱斷面的，說明土體的粒徑分布具有一定的各向異性。

不同土層土粒的形態(tài)分維數(shù)見圖8。整個淺層地基的土粒形態(tài)分維數(shù)為1.21～ 1.29，波動范圍非常小，平均值為1.25。從電鏡圖片可以看出，第⑦層黃土土粒較大，形態(tài)復(fù)雜，因此其形態(tài)分形維數(shù)較大，為1.28左右。隨著地基深度的增大，土粒的形態(tài)分形維數(shù)逐漸減小，主要原因是處于較深處的土體小粒徑含量一般比較高，而小粒徑一般具有較大的比表面積，在相互間的摩擦、破壞作用下其表面起伏減小，更容易趨向于圓球形，扁圓度升高，土粒形態(tài)的復(fù)雜程度降低，形態(tài)分形維數(shù)減小，同時土體變得更加密實，宏觀上土體強度提高。

4.4 土粒和孔隙的分形維數(shù)對比分析

結(jié)合圖4和圖7可以發(fā)現(xiàn)，深度方向上，土粒和孔隙的分布分維數(shù)均隨深度增大緩慢增大，這反映了土粒和孔隙兩者之間在直徑變化上具有一定的協(xié)調(diào)性：從第⑩層土到第11層土，土粒粒徑分布分維數(shù)增大變緩，說明粒徑的復(fù)雜化程度有減緩的趨勢，而孔徑分布分維數(shù)變化趨勢則基本不變。從整個淺層地基的角度來看，對于任一層土，土粒的粒徑分布分維數(shù)均大于同層土體孔隙的孔徑分布分維數(shù)，這說明土粒粒徑的復(fù)雜程度比孔隙孔徑的復(fù)雜程度要高，土粒細(xì)化程度比孔隙明顯。

結(jié)合圖6和圖8可見，土粒和孔隙的形態(tài)分形維數(shù)變化有相似的趨勢，這反映了土粒和孔隙之間在形態(tài)上具有一定的協(xié)調(diào)性：對于整個淺層地基，土粒的形態(tài)分維數(shù)變化范圍遠(yuǎn)小于孔隙的形態(tài)分維數(shù)變化范圍;對于同一層土，土粒的形態(tài)分維數(shù)小于孔隙的形態(tài)分維數(shù)，表明孔隙形態(tài)的復(fù)雜程度比土粒的高，一定程度上說明土粒的形態(tài)與孔隙相比更具有穩(wěn)定性。

4.5 分形維數(shù)與基本物理特性的關(guān)系

為研究土體基本物理特性指標(biāo)與相應(yīng)的分維數(shù)之間的關(guān)系，進行了多元回歸分析。

定義孔隙的孔徑分布分維數(shù)和形態(tài)分維數(shù)分別為D1和D2，土粒的粒徑分布分維數(shù)和形態(tài)分維數(shù)分別為D3和D4。根據(jù)統(tǒng)計學(xué)原理，考慮變量最高次數(shù)為二次的情況，將某項基本物理特性指標(biāo)F表示為

分析結(jié)果顯示，密度和滲透系數(shù)這兩個指標(biāo)的回歸模型關(guān)于D1、D2、D3和D4的所有組合均不滿足模型檢驗的顯著性要求，因此這兩個基本物理指標(biāo)無法用式（3）來表示，剩余6個指標(biāo)的回歸分析結(jié)果見表1。

由表1可知，土體的含水率、液塑限、孔隙比、黏聚力和內(nèi)摩擦角與相關(guān)的分維數(shù)間存在一定的相關(guān)性。其中：含水率、孔隙比分別與孔隙的形態(tài)分維數(shù)和土粒的形態(tài)分維數(shù)成線性負(fù)相關(guān)，說明孔隙的形態(tài)分維數(shù)是土體含水率高低程度的表現(xiàn)，孔隙的形態(tài)分布越整齊，土體的含水率越高;土粒的形態(tài)分維數(shù)是土體孔隙比變化的重要反映。黏聚力、內(nèi)摩擦角與粒徑的分布分維數(shù)和形態(tài)分維數(shù)相關(guān)，土體粒徑細(xì)化程度高，形態(tài)復(fù)雜，土粒表面起伏平緩，土粒間相互鑲嵌緊密，相應(yīng)的土體密實度高，土體強度增強。

液限取決于孔隙的孔徑分布分維數(shù)與土粒的形態(tài)分維數(shù)，而塑限與土粒的分維數(shù)相關(guān)，說明土粒的粒徑和形態(tài)的分布對土體塑限起了重要作用。原因是對于粒徑分布分維數(shù)大、形狀復(fù)雜的土體，其粒徑細(xì)化程度高，土粒的比表面積大，對水分子的吸附能力更強，需要吸附更多的水分子才能形成吸附水膜，進而提高了土體的塑限。當(dāng)含水率到達(dá)液限時，組成吸附水膜的強結(jié)合水黏滯性大，沒有流動性[21]，影響土體含水率變化的只有土體內(nèi)部自由水的含量，而自由水的含量跟孔隙的孔徑分布緊密相關(guān)，孔隙孔徑分布維數(shù)高的土體孔徑分布復(fù)雜，擁有較多的自由水填充通道，自由水含量高，故其液限相對較高。這就從微觀角度解釋了土體土粒和孔隙的分形特征與液塑限之間的相關(guān)性。

5 結(jié)論

對濟南地區(qū)不同土層土樣進行了一系列系統(tǒng)的試驗，并結(jié)合場發(fā)射掃描電鏡試驗，分析了不同土層土體微觀結(jié)構(gòu)特征，同時基于分形理論，探討了土體孑L隙和土粒的分形特征及其演化規(guī)律，主要結(jié)論如下：

（1）濟南地區(qū)淺層地基各土層的微觀結(jié)構(gòu)差異明顯。第⑦層黃土是以粗粉粒為主體骨架的蜂窩狀大孔隙結(jié)構(gòu)，第⑧層粉質(zhì)黏土是以碎屑粉粒為骨架的蜂窩狀結(jié)構(gòu)，第⑨、⑩和11層土是以黏粒集合體為基本單元的薄層狀絮凝結(jié)構(gòu)。

（2）各土層土體的微觀結(jié)構(gòu)具有明顯的分形特征，土體呈現(xiàn)明顯的各向異性，分形維數(shù)與土層埋深存在明顯相關(guān)性。隨深度加深，孔隙比、含水率、液塑限和黏聚力均一定程度增大，粒徑和孔徑的分布分維數(shù)均逐漸變大，形態(tài)分維數(shù)均逐漸變小，粒徑變化比孔徑變化明顯，而孔隙的形態(tài)變化比土粒的復(fù)雜。

（3）物理參數(shù)與分形參數(shù)具有相關(guān)性。含水率、孔隙比分別與孔隙的形態(tài)分維數(shù)和土粒的形態(tài)分維數(shù)成線性負(fù)相關(guān)，液限取決于孔隙的孔徑分布分維數(shù)與土粒的形態(tài)分維數(shù)，黏聚力、內(nèi)摩擦角與粒徑的分布分維數(shù)與形態(tài)分維數(shù)相關(guān)，塑限只與土粒的兩個分維數(shù)相關(guān)，以分形維數(shù)為中介可以較好地分析土體宏觀物理特性和微觀結(jié)構(gòu)間的聯(lián)系。

參考文獻(xiàn)：

[1] 王渭明，嚴(yán)紀(jì)興，趙增輝.黃河沖積地層深基坑支護優(yōu)化及應(yīng)用研究[J].地下空間與工程學(xué)報，2016，12（1）：176-180.

[2]劉宏偉，劉文白，崔勇濤，等，固化吹填土的滲透性及其分形維數(shù)的定量研究[J].科學(xué)技術(shù)與工程，2017，17（13）： 60-65.

[3] 周萃英，土體微觀結(jié)構(gòu)研究與土力學(xué)的發(fā)展方向：若干進展與思考[J].地球科學(xué)，2000（2）：215-220.

[4] AHMED A.Compressive Strength and Microstructure of SoftClay Soil Stabilized with Recycled Bassanite[J].AppliedClay Science， 2015， 104： 27-35.

[5] WU Zhilong， DENG Yongfeng， LIU Songyu， et al.Strengthand Micro-Structure Evolution of Compacted Soils Modifiedby Admixtures of Cement and Metakaolin[J].Applied ClayScience， 2016， 127-128：44-51.

[6] 陳琳，喻文兵，楊成松，等，基于微觀結(jié)構(gòu)的青藏高原風(fēng)積沙導(dǎo)熱系數(shù)變化機理研究[J].冰川凍土，2014，36（5）：1220-1226.

[7]石剛，支喜蘭，馬建秦，等，黃土的組分及微觀結(jié)構(gòu)特征與地基承載力[J].公路交通科技，2005，22（12）：46-49.

[8] 張先偉，孔令偉，王靜，針對黏性土膠質(zhì)聯(lián)結(jié)特征的SEM-EDS試驗研究[J].巖土力學(xué)，2013，34（增刊2）：195-203.

[9]徐日慶，鄧袆文，徐波，等，基于SEM圖像的軟土三維孔隙率計算及影響因素分析[J].巖石力學(xué)與工程學(xué)報，2015 ，34（7）：1497-1502.

[10] 胡海軍，蔣明鏡，彭建兵，等，應(yīng)力路徑試驗前后不同黃土的孔隙分形特征[J].巖土力學(xué)，2014，35（9）：2479-2485.

[11]趙月平，膨脹土分類的分形插值評價方法[J].人民黃河，2015，37（6）：91- 94.

[12]劉松玉，張繼文，土中孑L隙分布的分形特征研究[J].東南大學(xué)學(xué)報，1997（3）：129-132.

[13] 呂海波，趙艷林，孔令偉，等，海相軟土孑L隙分布的分形特征及應(yīng)用[J].桂林工學(xué)院學(xué)報，2006（1）：50-53.

[14] LU Yang， LIU Sihong， WENC Liping， et al.Fractal Anal-ysis of Cracking in a Clayey Soil Under Freeze - Thaw Cycles[J].Engineering Geology， 2016， 208：93-99.

[15] XU Yongfu， CAO Ling. Fractal Representation for EffectiveStress of Unsaturated Soils[J].Intemational Joumal of Ge-omechanics， 2014， 15（6）：98.

[16] ALFARO SOTO M A， CHANC H K， VAN GENUCHTENM T. Fractal-Based Models for the Unsaturated Soil Hy-draulic Functions[J].Ceoderma， 2017， 306： 144- 151.

[17] 中華人民共和國水利部，土工試驗方法標(biāo)準(zhǔn)：GB/T50123-1999[S].北京：中國標(biāo)準(zhǔn)出版社，1999：50-52.

[18] 張季如，祝杰，黃麗，等，土壤微觀結(jié)構(gòu)定量分析的IPP圖像技術(shù)研究[J].武漢理工大學(xué)學(xué)報，2008，30（4）：80-83.

[19] 陶高梁，張季如，表征孑L隙及顆粒體積與尺度分布的兩類巖土體分形模型[J].科學(xué)通報，2009，54（6）：838-846.

[20] MIAO Tongjun， YANC Shanshan， LONG Zhangcai， et al.Fractal Analysis of Permeability of Dual-Porosity MediaEmbedded with Random Fractures [J]. IntemationalJoumal of Heat and Mass Transfer， 2015， 88：814-821.

[21] 黨發(fā)寧，劉海偉，王學(xué)武，等，基于有效孔隙比的黏性土滲透系數(shù)經(jīng)驗公式研究[J].巖石力學(xué)與工程學(xué)報， 2015 ，34（9）：1909-1917.