劉錦陽，李喜安,2，簡 濤，郭澤澤

(1.長安大學(xué)地質(zhì)工程與測繪學(xué)院，陜西 西安 710054；2.國土資源部巖土工程開放研究實驗室，陜西 西安 710054 )

馬蘭黃土滲氣率與飽和滲透系數(shù)的關(guān)系研究

劉錦陽1，李喜安1,2，簡 濤1，郭澤澤1

(1.長安大學(xué)地質(zhì)工程與測繪學(xué)院，陜西 西安 710054；2.國土資源部巖土工程開放研究實驗室，陜西 西安 710054 )

利用改進(jìn)的ZC-2015型滲氣儀和TST-55型滲透儀進(jìn)行馬蘭黃土滲氣和飽和滲透試驗，確定穩(wěn)定滲氣率對應(yīng)的穩(wěn)定滲徑，進(jìn)而探討滲氣率ka和飽和滲透系數(shù)Kw的關(guān)系。研究結(jié)果表明：ka不僅能描述非飽和土孔隙和結(jié)構(gòu)特征，也能用于預(yù)測Kw。氣體滲氣率隨著滲徑的增大而減小，最終趨于穩(wěn)定，原狀和重塑黃土的穩(wěn)定滲徑均約8 cm。ka和Kw間呈明顯的雙對數(shù)線性關(guān)系，不同深度的原狀風(fēng)干黃土和不同粒組的重塑土的lgka和lgKw間的相關(guān)性都比較明顯，但原狀風(fēng)干黃土和重塑土之間的擬合公式有很大差異。另外，黏粒含量較高時，重塑土擬合直線的斜率β和截距α明顯增大。不同含水率下，重塑土擬合公式有很大的差異，當(dāng)含水率較大時，隨著干密度的增大，ka的變化程度比Kw大，擬合直線的斜率β和截距α都有明顯的減小。

滲氣率；穩(wěn)定滲徑；飽和滲透系數(shù)；雙對數(shù)線性關(guān)系; 馬蘭黃土; 黏粒含量

馬蘭黃土(Q3)作為一種分布廣且厚度穩(wěn)定(一般厚度10～30 m)的多孔材料[1]，其疏松多孔的結(jié)構(gòu)特性為水、氣在土體中的運移、儲存提供了天然通道和空間。由于其位于更新統(tǒng)頂層且具有很好的水敏性，所以工程建設(shè)中經(jīng)常會遇到水的滲透問題，目前常用單環(huán)法和雙環(huán)法測量原位土的滲透系數(shù)，室內(nèi)滲透試驗常用TST-55型滲透儀法，但這兩種方法常因耗時費力不能較好地滿足實際工程要求。滲氣率ka的測量相對簡單、快速，因此，很多學(xué)者曾探討過利用滲氣率ka表征滲透系數(shù)Kw。

目前沒有專門的儀器測量馬蘭黃土的滲氣率，只有一些改進(jìn)的裝置用于測量土壤、砂土及混合材料的滲氣性。有些學(xué)者利用改進(jìn)的三軸滲氣儀研究了含水率、干密度、體積含氣率、應(yīng)力、材料的各向異性、混合材料的摻料率對滲氣系數(shù)的影響[2～5]。劉奉銀等[6]利用改進(jìn)的水氣運動聯(lián)合測定儀得到了溫度和密度雙變化條件下非飽和黃土的滲水和滲氣函數(shù)。施建勇等[7]利用研制的非飽和滲透試驗儀探究了垃圾土氣體滲透系數(shù)隨著有機質(zhì)含量、飽和度、孔隙比、氣壓梯度的變化規(guī)律。楊益彪[8]用自制的測量裝置測試黃土的氣體滲透系數(shù)，試驗結(jié)果表明高壓實度黃土的氣體滲透系數(shù)隨服役含水率的變化明顯。彭爾興等[9]利用改進(jìn)的柔性壁滲透儀研究了高滲透率含砂細(xì)粒土的氣體滲透規(guī)律。這些研究所涉及的試驗裝置只能滿足室內(nèi)滲透試驗，且只是對不同材料的滲氣規(guī)律和滲氣系數(shù)的影響因素進(jìn)行過一些探討，但都沒有研究滲氣率和飽和滲透系數(shù)之間的關(guān)系。

國外學(xué)者對于滲透率ka和飽和滲透系數(shù)Kw之間的關(guān)系有一定的研究，研究者們利用回歸分析模型對ka和Kw的關(guān)系進(jìn)行探討，發(fā)現(xiàn)ka和Kw之間為雙對數(shù)線性關(guān)系[10～14]，但是研究對象為土壤、砂、黏土等材料，對典型馬蘭黃土ka和Kw之間關(guān)系的研究較少。為了研究典型馬蘭黃土的氣體滲透特性，本文主要進(jìn)行馬蘭黃土的室內(nèi)滲氣試驗和飽和滲透試驗，探究馬蘭黃土的滲氣和滲透之間的關(guān)系，本文得出的滲透和滲氣參數(shù)及相關(guān)規(guī)律，對研究馬蘭黃土的水、氣滲透問題以及與滲透問題有關(guān)的工程建設(shè)具有重要的參考價值。

1 試驗材料和方法

1.1試驗材料

本次土樣主要取自粉黃土和黏黃土帶中的甘肅黑方臺、陜西延安、陜西浐河附近的典型風(fēng)積黃土剖面，3個取樣點位于圖1中黃土高原的粉黃土帶和黏黃土帶[15]。

通過試驗得到三種土的基本物理力學(xué)指標(biāo)，見表1，三種土的顆粒組成差異較大，雖然三種黃土都主要以粉粒為主，但黑方臺黃土的砂粒含量明顯高于其它兩種土，三種土中黏粒含量最高的是浐河剖面黃土。

圖1 馬蘭黃土粒度水平分帶及取樣點圖Fig.1 Granularity level zonation of the Malan loess and sampling location

1.2試驗方法及方案

本次試驗采用的滲氣儀是在ZC-2015型滲氣儀的基礎(chǔ)上改進(jìn)的。ZC-2015型滲氣儀主要用于瀝青混凝土防滲測試[16～17]，由于其實際測量范圍有限，故需要對其進(jìn)行改進(jìn)，才可能滿足非飽和黃土的滲氣試驗。本次儀器的改進(jìn)主要為以下兩個方面。

(1) 在原有的小氣室基礎(chǔ)上增加一個10 000 mL的大氣室，擴(kuò)大儀器的測試范圍，氣室1可以在干密度較大的情況下測量黃土的滲氣性，氣室2能夠滿足一般工況下非飽和黃土的滲氣試驗。

(2) 根據(jù)大量黃土滲氣試驗結(jié)果，對儀器重新率定，并改進(jìn)原有的內(nèi)置程序。改進(jìn)后的滲氣儀能夠方便、簡捷的進(jìn)行非飽和黃土的滲氣測量。改進(jìn)的滲氣儀原理圖見圖2。

表1 埋深4 m的3種土的基本物理指標(biāo)

圖2 改進(jìn)的ZC-2015型滲氣儀構(gòu)造簡圖Fig.2 Structure diagram of the improved ZC-2015 air Permeameter1—緩沖腔；2—汞入口；3—標(biāo)尺；4—U形管；5—紅外光電組合模塊；6—智能測控系統(tǒng)；7—1 000 mL氣室1；8—閥組(閥A，閥B)；9—閥C；10—真空泵；11—閥D；12—10 000 mL氣室2；13—閥E；14—蓄電池

改進(jìn)滲氣儀的主要測試原理為：關(guān)閉閥E并打開閥B和真空泵開關(guān)閥C開始抽真空，將U形管中的汞柱抽到預(yù)定高度，氣室和與連接管組成的系統(tǒng)(下文簡稱密閉系統(tǒng))呈現(xiàn)負(fù)壓，與外界大氣壓形成一定壓差。試驗時打開閥E，空氣直接通過土樣進(jìn)入密閉系統(tǒng)，汞柱在壓差的作用下開始下降，汞柱先后會經(jīng)過標(biāo)尺上標(biāo)定好的P0～P3四個點，P0到P3點壓差恒定，進(jìn)入密閉系統(tǒng)的空氣量恒定，汞柱到達(dá)P3試驗結(jié)束，儀器會記錄汞柱從P0到P3點的時間，不同原狀風(fēng)干黃土和重塑土在恒定進(jìn)氣量下汞柱從P0到P3點的時間不同，通過推導(dǎo)的公式計算出土體滲氣率。

1.3滲氣率的計算公式

抽真空后，打開閥E開始試驗，設(shè)任一時刻t，水銀柱高度為h，室溫T=20 ℃，經(jīng)過Δt時間，U形管內(nèi)水銀柱的高度變化量Δh，則儀器內(nèi)部氣壓為

式中：P——t時刻儀器內(nèi)部氣壓；

Patm——標(biāo)準(zhǔn)大氣壓；

ρHg——汞的密度。

由理想氣體的狀態(tài)方程得Δt時間后密閉系統(tǒng)內(nèi)氣體物質(zhì)的量為

式中：VS——密閉系統(tǒng)的體積；

n——t時刻氣體的物質(zhì)的量。

聯(lián)立式(1)和(2)得儀器內(nèi)部壓力變化量和氣體物質(zhì)的量的變化量的微分表達(dá)式分別為：

外界氣體通過測試土樣進(jìn)入密閉系統(tǒng)的氣體物質(zhì)的量的微分表達(dá)式為：

由密閉系統(tǒng)內(nèi)、外氣體的物質(zhì)的量守恒可得氣體的進(jìn)氣量dVg為：

根據(jù)氣體達(dá)西定律，多孔介質(zhì)一維體積流速為：

式中：k——氣體滲透率；

E——空氣粘滯系數(shù)；

L——氣體滲徑；

A——氣體滲透面積。

對(7)式積分化簡得滲氣率為：

式中：h0——P0點汞柱在標(biāo)尺上的刻度；

lnho/h——通過h-t的率定實測曲線得出。

2 試驗結(jié)果與討論

2.1滲氣率與土樣滲徑的關(guān)系

滲氣率屬于土樣的一種基本物理屬性，在試驗過程中滲徑的合理選擇對于滲氣率的確定至關(guān)重要，滲徑選擇過大會造成制樣和測試的極不方便，而滲徑選擇過小造成滲氣率的測量值比實際值偏大，不同滲徑下的滲氣率只是滲氣率的測量值?？紤]到滲氣率測量值同時受土樣結(jié)構(gòu)、干密度、含水率等因素的影響，設(shè)計了相關(guān)的試驗工況，以確定既便于制樣，又能滿足測量精度的合理滲徑。試驗主要通過制備不同高度的土樣控制滲徑，對于原狀風(fēng)干黃土，利用自制環(huán)刀(φ=8 cm)削樣得到不同高度的原狀樣；對于重塑土，利用壓樣器壓實得到不同滲徑的重塑樣。在固相邊界，不存在完全的氣體零流速“固壁”情況，會產(chǎn)生“氣體滑脫”現(xiàn)象[19]。在制樣時，考慮氣體易從環(huán)刀壁繞流和土樣邊緣泄露，因此在環(huán)刀周圍涂抹一層凡士林，主要起到密封作用。

浐河剖面的原狀風(fēng)干黃土(w=1.37%)和重塑土(w=9.6%)的滲氣率與滲徑的關(guān)系見圖3，風(fēng)干原狀黃土和重塑土的滲氣率測量值都是隨著滲徑的增大而減小，最終趨于穩(wěn)定。原狀風(fēng)干黃土和重塑黃土的穩(wěn)定滲徑在8 cm左右，由于原狀風(fēng)干黃土自身的結(jié)構(gòu)性差異，滲徑較大時結(jié)構(gòu)性差異有一定程度的減小，滲氣率測量值隨著滲徑的增大波動性下降，而重塑土的結(jié)構(gòu)和膠結(jié)差，其滲氣率測量值隨著滲徑平穩(wěn)變化。

圖3 浐河剖面黃土滲氣率測量值與滲徑的關(guān)系曲線Fig.3 Relationship of air permeability measurements and permeability diameter of loess in the Chanhe river profile

圖4所示為不同滲徑的重塑土樣。黑方臺和延安重塑黃土的滲氣率測量值和土樣的滲徑關(guān)系曲線如圖5(a)～5(d)所示，隨著滲徑的增大，土樣的滲氣率測量值越小，含水率越高，滲氣率測量值和滲徑的變化規(guī)律越明顯。不同干密度和含水率下，土樣的穩(wěn)定滲徑基本為8 cm，滲徑越小時，氣體通過土樣的時間越短，滲氣率測量值的穩(wěn)定性越差。

圖4 不同滲徑的土樣Fig.4 Soil samples with different permeability diameters

圖5(e)為浐河重塑黃土滲氣率測量值和土樣滲徑的關(guān)系曲線，含水率在7.4%時，滲氣率測量值與滲徑的關(guān)系不明顯，曲線為波動狀。重塑土樣滲徑從4 cm變化到10 cm時，土樣的滲氣率測量值變化范圍較小，基本在同一個量級內(nèi)變化，不同含水率下，土樣穩(wěn)定的滲徑基本為8 cm。含水率越大，隨著滲徑的增大，滲氣率測量值減小的趨勢越明顯。因為基質(zhì)吸力為孔隙氣壓和孔隙水壓的差值，方向與孔隙氣壓相反，在一定程度上對土體滲氣有阻滯作用，而重塑土的基質(zhì)吸力隨著土樣含水率的增大而減小，所以，隨著含水率的增大，阻滯作用會逐漸減小，導(dǎo)致隨著含水率的增大滲氣率測量值增大，且滲氣率測量值和滲徑的趨勢隨含水率的增大而更明顯。

圖5 重塑土滲氣率測量值與滲徑的關(guān)系曲線Fig.5 Relationship between air permeability measurements and sample permeability diameters of the remold loess

圖6為三種重塑土在相同工況下滲氣率測量值與滲徑的關(guān)系曲線，浐河重塑黃土的滲氣率測量值與其它兩種土的滲氣率測量值相差一個量級，主要是因為浐河黃土中的黏粒含量較高，土中主要以小孔隙為主，土中有效滲氣的大孔隙較少，滲氣率測量值較低。原狀風(fēng)干黃土和重塑土的滲徑大于8 cm時，滲氣率測量值基本不變，則滲氣率對應(yīng)土樣的穩(wěn)定滲徑為8 cm，所以不同工況下重塑土樣的滲徑都為8 cm，這樣既使制樣方便，又能降低測量誤差。

圖6 三種重塑土滲氣率測量值與滲徑關(guān)系對比曲線Fig.6 Comparison between air permeability measurements and sample permeability diameters of three kinds of loess

2.2滲氣率和飽和滲透系數(shù)的關(guān)系

土的滲氣性是氣體通過多孔介質(zhì)孔隙的性質(zhì)，而土的水力傳導(dǎo)性是水通過相同多孔介質(zhì)孔隙的性質(zhì)，因此，土的滲氣率一定程度上可以指示土的水力傳導(dǎo)性[20]。Loll[21]利用9種不同地區(qū)的土， 這9種土都以砂粒為主，粉粒和黏粒含量有一定的差異，在95%的置信區(qū)間下得出9種土的滲氣率和飽和滲透系數(shù)(下文簡稱滲透系數(shù))之間通用關(guān)系式。為了探究滲透系數(shù)Kw和滲氣率ka之間的關(guān)系，原狀黃土和重塑黃土的滲氣試驗使用環(huán)刀(φ=8 cm，h=8 cm)制樣,滲透試驗使用環(huán)刀(φ=6 cm,h=4 cm)制樣。

由于ka和Kw之間存在較大的量級差，為了消除這種量級差且更好的描述兩者之間的關(guān)系，所以先對滲透系數(shù)進(jìn)行單位換算，再對兩種滲透參數(shù)進(jìn)行對數(shù)化，使對數(shù)化后的值在同一個量級內(nèi)。描述非飽和土lgka和lgKw之間的函數(shù)關(guān)系一般用簡單的回歸模型——一元線性回歸模型，該模型不考慮時間和空間因素，適合于一般的室內(nèi)試驗情況，且主要從數(shù)據(jù)出發(fā)利用最小二乘法表現(xiàn)lgka和lgKw之間的函數(shù)關(guān)系：

式中：α——直線截距；

β——直線斜率。

2.2.1原狀黃土的滲氣率和飽和滲透系數(shù)的關(guān)系

圖7為延安和黑方臺原狀風(fēng)干黃土在2～10 m 5個深度下，滲氣率和滲透系數(shù)隨深度的變化曲線，可以看出滲氣率和滲透系數(shù)都是隨著深度的增大而減小。因為黃土是風(fēng)積作用形成的，在深度上其粒徑分布為上粗下細(xì)，馬蘭黃土粒級在垂向剖面上也是從上到下有逐漸變細(xì)的趨勢，延安黃土中砂粒含量隨深度減小，黏粒含量隨深度增加，砂粒和黏粒粒組在垂向上相互消長，粉粒顆粒含量保持相對穩(wěn)定，其含量在70%左右(表2)。圖7(a)中延安4 m的原狀黃土相比其它深度的黃土，肉眼可見的大孔和利于水和氣滲透的有效小孔隙較少，土中裂隙較少，結(jié)構(gòu)相對致密，其滲氣率和滲透系數(shù)相對較低，4 m處的曲線呈下凹狀。另外，兩種土的滲氣率和滲透系數(shù)隨深度變化趨勢基本一致，可見這兩種參數(shù)之間存在明顯的相關(guān)性，即一種對應(yīng)的函數(shù)關(guān)系。

圖7 原狀風(fēng)干黃土的滲透參數(shù)與深度的關(guān)系曲線Fig.7 Relationship between permeability parameters and depth of the undisturbed dry soils

圖8為黑方臺和延安不同深度的原狀風(fēng)干黃土lgka和lgKw間的線性關(guān)系曲線，黑方臺黃土的數(shù)據(jù)整體比較分散，數(shù)據(jù)的均勻性較差。由于原狀風(fēng)干黃土的結(jié)構(gòu)性和顆粒組成之間的差異，兩種土線性公式的斜率和截距有一定的差異。黑方臺黃土兩種滲透參數(shù)間的擬合度更高，且擬合公式更接近Loll的通用擬合公式y(tǒng)=1.27x+14.11，因為黑方臺黃土和Loll試驗用的土都是砂土，只是不同粒組的顆粒含量不同。

表2 延安黃土基本物理指標(biāo)

2.2.2重塑黃土的滲氣率和飽和滲透系數(shù)的關(guān)系

延安10 m重塑土的滲氣率和飽和滲透系數(shù)與干密度的變化趨勢如圖9所示，ka和Kw變化趨勢基本一

圖8 原狀風(fēng)干黃土滲氣率和飽和飽和滲透系數(shù)關(guān)系曲線Fig.8 Relationship between air permeability and saturated permeability coefficient of the undisturbed dry loess

圖9 重塑黃土滲透參數(shù)與干密度的關(guān)系曲線Fig.9 Relationship of permeability parameters and dry density of the remolded loess

致，整體變化趨勢都是隨著干密度的增大而減小，干密度越大，重塑土越密實，單位體積內(nèi)的土骨架顆粒越多，有效滲氣孔隙越少，ka和Kw越小。含水率較大時，變化趨勢越明顯。具體的滲透參數(shù)如表3所示。綜上可知重塑土ka和Kw之間也存在著一定的可擬合的關(guān)系。

表3 延安10 m重塑馬蘭黃土的滲透參數(shù)表Table 3 Coefficient of permeability of the Yan’an 10 m remodeling Malan loess

由圖10可知，不同含水率下，擬合公式有很大的差異，隨著含水率的增大，斜率和截距都有一定程度的減小。w大于11.8%時，斜率值都小于1，斜率明顯降低，因為w較大時，ka隨著干密度的變化程度比Kw大，擬合直線比低含水率下的直線平緩。

圖10 延安10 m重塑土滲氣率和飽和滲透系數(shù)關(guān)系曲線Fig.10 Relationship between air permeability and saturated coefficient of permeability of the Yan’an 10 m undisturbed loess

利用黑方臺、延安、浐河剖面三種地區(qū)的4 m黃土，對比砂粒、粉粒、黏粒含量不同的黃土在含水率一定(w=9.6%)干密度變化時，不同重塑的lgka和lgKw的關(guān)系，見圖11，三種土的lgka和lgKw間都呈現(xiàn)線性關(guān)系。如浐河黃土的擬合曲線所示，黏粒含量越高時，直線的斜率β截距α明顯增大，因為黏粒含量較高時，重塑土的ka隨著干密度的變化不明顯，而Kw變化明顯。對比可知，重塑土數(shù)據(jù)的均勻性比原狀風(fēng)干黃土好，說明重塑樣均勻性較好，在滲透和滲氣的測量時誤差較小。的斜率范圍較大。

圖11 w=9.6%時三種重塑黃土滲氣率和飽和滲透系數(shù)關(guān)系曲線Fig.11 Relationships between air permeability and saturated coefficient of permeability of three remold loesses, w=9.6%

3 結(jié)論

(1)試驗發(fā)現(xiàn)，原狀黃土和重塑黃土的滲氣率均隨著滲徑的增大而減小最終趨于穩(wěn)定。土樣含水率越大時，這種變化趨勢越明顯；顆粒組成不同的土樣，黏粒含量越高，變化趨勢越明顯。原狀黃土和重塑黃土滲氣率穩(wěn)定時對應(yīng)的土樣滲徑都為8 cm。

(2)黑方臺和延安不同深度的原狀風(fēng)干黃土lgka和lgKw間的呈現(xiàn)線性關(guān)系，由于兩種黃土的結(jié)構(gòu)性和顆粒組成之間的差異，使兩種土lgka和lgKw線性公式的斜率和截距有一定的差異。黑方臺黃土兩種滲透參數(shù)擬合公式的擬合度更高，且擬合公式最接近Loll的通用擬合公式y(tǒng)=1.27x+14.11。

(3)不同含水率下，延安10 m重塑土的lgka和lgKw擬合公式有很大的差異，隨著含水率的增大，斜率β和截距α都有一定程度的減小。含水率大于11.8%時，斜率明顯降低，含水率w較大時，隨著干密度增大ka的變化程度比Kw大，擬合直線的斜率比低含水率下的直線斜率低。

(4)不同粒組重塑土的lgka和lgKw間都呈現(xiàn)線性關(guān)系。黏粒含量較高時，直線的斜率β截距α明顯增大，黏粒含量較高時，重塑土的Kw隨著干密度的變化比ka明顯，重塑土數(shù)據(jù)的均勻性比風(fēng)干原狀黃土好，重塑樣均勻性較好，對比發(fā)現(xiàn)原狀風(fēng)干黃土和重塑土之間的擬合公式有很大的差異。

[1] 孫建中. 黃土名詞問題[J]. 地球科學(xué)與環(huán)境學(xué)報，1985(2)：42-46. [SUN J Z. Noun problem of loess[J]. Journal of Earth Science and the Environment，1985(2)：42-46.(in Chinese)]

[2] 苗強強，陳正漢，張磊，等. 非飽和黏土質(zhì)砂的滲氣規(guī)律試驗研究[J]. 巖土力學(xué)，2010，31(12)：3746-3750. [MIAO Q Q，CHEN Z H，ZHANG L，etal. Experimental study of gas permeability of unsaturated clayey sand[J]. Rock and Soil Mechanics，2010，31(12)：3746-3757.(in Chinese)]

[3] 姚志華，陳正漢，黃雪峰，等. 非飽和Q3黃土滲氣特性試驗研[J].巖石力學(xué)與工程學(xué)報，2012,31(6)：1264-1273. [YAO Z H，CHEN Z H，HUANG X F，etal. Experimental research on gas permeability of unsaturated Q3loess[J]. Chinese Journal of Rock Mechanics and Engineering，2012，31(6)：1264-1273.(in Chinese)]

[4] 汪龍，方祥位，申春妮，等. 膨潤土-砂混合型緩沖/回填材料滲氣規(guī)律試驗研究[J]. 巖石力學(xué)與工程學(xué)報，2015，34(增刊1)：3381-3388. [WANG L，F(xiàn)ANG X W，SHEN C N，etal. Experimental study of gas permeability of bentonite-sand mixtures as buffer/back filling materials[J]. Chinese Journal of Rock Mechanics and Engineering，2015，34(Sup1)：3381-3388.(in Chinese)]

[5] 陳存禮，張登飛，張潔，等. 等向應(yīng)力條件下原狀Q3黃土的滲氣特性研究[J]. 巖土工程學(xué)報，2017，39(2)：287-294. [CHEN C L，ZHANG D F，ZHANG J，etal. Gas permeability of intact Q3loess under isotropic stress[J]. Chinese Journal of Geotechnical Engineering，2017，39(2)：287-294.(in Chinese)]

[6] 劉奉銀，張昭，周冬. 濕度和密度雙變化條件下的非飽和黃土滲氣滲水函數(shù)[J]. 巖石力學(xué)與工程學(xué)報，2010，29(9)：1907-1914. [LIU FY，ZHANG Z，ZHOU D. Density-saturation-dependent air-water permeability function of unsaturated loess[J]. Chinese Journal of Rock Mechanics and Engineering，2010，29(9)：1907-1914.(in Chinese)]

[7] 施建勇，趙義. 氣體壓力和孔隙對垃圾土體氣體滲透系數(shù)影響的研究[J]. 巖土工程學(xué)報，2015，37(4)：586-593. [SHI J Y，ZHAO Y. Influence of air pressure And void on permeability coefficient of air in municipal solid waste(MSW)[J]. Chinese Journal of Geotechnical Engineering，2015，37(4)：586-593.(in Chinese)]

[8] 楊益彪，詹良通，王順玉，等. 封頂覆蓋黃土氣體滲透特性測試及填埋氣一維運移分析[J]. 巖土力學(xué)，2015，36(7)：1973-1980. [YANGY B，ZHANG L T，WANG S Y，etal. Measurement of gas permeability of cover loess and one-dimensional analysis of land fill gas emission[J].Rock and Soil Mechanics，2015，36(7)：1973-1980. (in Chinese)]

[9] 彭爾興，孫文博，章定文，等. 非飽和含砂細(xì)粒土的氣體滲透特性研究[J]. 巖土力學(xué)，2016，37(5)：1301-1306. [PENG E X，SUN W B，ZHANG D W，etal. Air permeability of unsaturated fine sandy soil[J].Rock and Soil Mechanics，2016，37(5)：1301-1306.(in Chinese)]

[10] Iversen B V，Moldrup P，Schj?nning P，etal. Air and water permeability in differently textured soils at two measurement scales[J]. Soil Science，2001，166：643-659.

[11] Chief K, Ferre T, Nijssen B. Correlation between air permeability and saturated hydraulic conductivity: Unburned and burned soils[J]. Soil Science Society of America Journal, 2008,72(6):1501-1509.

[12] Mohammad R N，Mehdi R，Seyed A R A，etal. Prediction of unsaturated soil hydraulic conductivity using air permeability: Regression approach[J]. Indian J Agric Rec，2015，49(6)：528-533.

[13] Huang M B，Rodger H，Barbour S L. An evaluation of air permeability measurements to characterize the saturated hydraulic conductivity of soil reclamation covers[J]. Can J Soil Sci，2015，95：15-26.

[14] HUANG M B，Zettl J D，Barbour S L，etal. Characterizing the spatial variability of the hydraulic conductivity of reclamation soils using air permeability[J]. Geoderma，2016，262：285-293.

[15] 劉東生. 黃土的物質(zhì)成分與結(jié)構(gòu)[M]. 北京：科學(xué)出版社，1966. [LIU D S. The material composition and structure of loess[M].Beijing：Science Press，1966.(in Chinese)]

[16] 丁樸榮. ZC-5型滲氣儀的研制與滲氣測試方法的研究[J]. 西安理工大學(xué)學(xué)報，1986(4)：37-52. [DING P R. The research on ZC-5 air permeameter development and Air permeability test method[J]. Journal of Xi’an University of Technology，1986(4)：37-52.(in Chinese)]

[17] 丁樸榮. 瀝青混凝土斜墻滲氣檢驗的理論與實踐[J]. 水力發(fā)電學(xué)報，1989(4)：86-97. [DING P R. The test of asphalt concrete inclined wall air permeability theory and practice[J]. Journal of hydroelectric power，1986(4)：86-97.(in Chinese)]

[18] Marc J，Jacob H D. A handheld device for instrusive and nonintrusive field measurements of air permeability[J].Vadose Zone Journal，2003(2)：611-617.

[19] Klinkenberg L J. The permeability of porous media to liquids and gases[C]//Drilling and Production Practice. Washington D C：[s.n.], 1941：200-213.

[20] Blackwell P S，Ringrose A J，Jayawardan N S，etal. The use of air-filled porosity and intrinsic permeability to air to characterize structure of macro pore space and saturated hydraulic conductivity of clay soils[J]. Soil Science，1990,41：215-228.

[21] Loll P，Moldrup P，Schj?nning P，etal. Predicting saturated hydraulic conductivity from air permeability：applcation in stochastic water infiltration modeling[J]. Water Resour Res，1999，35：2387-2400.

責(zé)任編輯

：張若琳

AstudyoftherelationshipbetweengaspermeabilityandsaturatedcoefficientofpermeabilityoftheMalanloess

LIU Jinyang1, LI Xi’an1,2, JIAN Tao1, GUO Zeze1

(1.SchoolofGeologicalEngineeringandGeomatics,Chang’anUniversity,Xi’an,Shaanxi710054,China; 2.OpenResearchLaboratoryofGeotechnicalEngineering,MinistryofLandandResources,Xi’an,Shaanxi710054,China)

A series of gas and saturated penetration tests for the Malan loess are conducted by using the improved ZC-2015 air permeameter and TST-55 permeameter. We firstly determine the stable permeability diameter of the correspondingka, and then explore the relationship between air permeabilitykaand saturation coefficient of permeabilityKw. The research results show thatkanot only describes the pore and structural characteristics of unsaturated soils, but also can be used to predictKw.kadecreases with the increasing permeability diameter, and ultimately this trend becomes stabilized. The stable permeability of the undisturbed and remolded loess is 8 cm. There is a significant double-logarithmic linear relationship betweenkaandKw, The correlation between lgkaand lgKwof the remolded loess with different grain compositions and undisturbed loess in different depths are obvious, but the fitting formula between the undisturbed soil and the remolded soil is very different. In addition, when the content of clay particles is high, the slopeβand interceptαof the remolded sample obviously increase. Under the different water content, the formula of the remolded soil is very different. when the water content, the change inkais more thanKwwith the increase of the dry density, and the slopeβand the interceptαof fitting line are obviously reduce.

gas permeability; stable permeability diameter; saturated permeability coefficient; double-logarithmic linear relationship; Malan loess; content of day particles

10.16030/j.cnki.issn.1000-3665.2017.06.23

P642.13+1；TU 411

A

1000-3665(2017)06-0154-09

2017-00-00;

2017-00-00

國家自然科學(xué)基金資助項目(41572264，41172255，41440044)

劉錦陽(1991-)，男，碩士，主要從事黃土工程地質(zhì)及地質(zhì)災(zāi)害防治。E-mail：2296454851@qq.com

李喜安(1968-)，男，教授，博士生導(dǎo)師，主要從事黃土地質(zhì)災(zāi)害方面的教學(xué)與科研工作。E-mail：dclixa@chd.edu.cn