基于正交設計的泥巖膨脹試驗及其影響因素研究

摘 要：以正交試驗設計方法為基礎，對張家口市京新高速公路路塹膨脹性泥巖開展試驗， 重點研究不同上覆荷載下泥巖試樣的膨脹率與浸水時間、初始含水率、壓實度的規律關系，并采用極差分析法對各種因素進行敏感性分析。研究結果表明，膨脹率隨浸水時間的變化曲線呈現出S形走勢，大致可劃分為初始膨脹階段、加速膨脹階段、緩慢膨脹階段；在相同初始含水率條件下，膨脹率隨著壓實度的增大而增大；在相同的壓實度條件下，膨脹率隨初始含水率增大而減小。通過極差分析法可知，試樣浸水時間為影響膨脹率的主導因素，含水率和壓實度為次要影響因素，且小荷載下含水率對試樣膨脹率的影響比較明顯，荷載繼續增大，壓實度超過含水率對試樣膨脹率的影響。

關鍵詞：泥巖； 膨脹變形； 膨脹率； 正交試驗； 極差分析法

中圖分類號：TU41 文獻標識碼：A DOI：10.7525/j.issn.1006-8023.2024.05.018

Mudstone Swelling Test and its Influencing Factors Based on Orthogonal Design

Abstract： Based on the orthogonal test design method， the expansive mudstone in the cutting of Jingxin expressway in Zhangjiakou city was tested， and the relationship between the expansion rate of mudstone samples under different overlying loads and the soaking time， initial water content and compaction degree was studied， and the sensitivity of various factors was analyzed by the range analysis method. The results showed that the curve of expansion rate with the soaking time presented an S-shaped trend， which can be roughly divided into initial expansion stage， accelerated expansion stage and slow expansion stage. Under the same initial water content， the expansion rate increased with the increase of the compaction degree. Under the same compaction degree， the expansion rate decreased with the increase of initial water content. The range analysis method showed that the soaking time of the sample was the main factor affecting the expansion rate. Water content and compaction degree were the secondary influencing factors， and the influence of water content on the sample expansion rate was obvious under small load， and the load continued to increase， the compaction degree exceeded the influence of water content on the sample expansion rate.

Keywords： Mudstone； expansion deformation； expansion rate； orthogonal test； range analysis method

0 引言

膨脹性泥巖［1］是一類在地質工程中常見的軟巖，主要由蒙脫石和伊利石組成，具有吸水膨脹軟化和失水收縮干裂的特性。泥巖的特殊性質可造成地基破壞，引起建筑物的開裂和上拱、路基、渠道的破壞，以及邊坡產生滑坡失穩等不良現象。然而影響泥巖膨脹變形的因素有很多，主要與泥巖的起始狀態有關，其中膨脹率是衡量泥巖膨脹變形能力的重要指標之一，許多學者做了大量的研究。Chertkov［2］建立膨脹巖膨脹和收縮的物理模型，并推導相應的公式，定量化敘述了膨脹和收縮的變化；屈祥［3］開展不同初始含水率下湘西弱膨脹土的膨脹變形時程特性研究，分析上覆荷載及含水率對湘西膨脹土膨脹變形和膨脹速率的影響；張愛軍等［4］通過對不同初始干重度、初始含水率和上覆壓力的一系列膨脹量的試驗，得到了人工壓實膨脹土的膨脹變形隨3種因素變化的規律；徐永福等［5］根據不同含水量、干密度的膨脹土在不同壓力下的膨脹變形特性，研究了膨脹土的膨脹量和膨脹力與起始含水量、干密度和壓力之間的相關關系；饒錫保等［6］對南水北調中線工程南陽中膨脹土擊實樣進行了不同起始含水率與壓實度的無荷膨脹率和有荷膨脹率試驗，采用正交分析法分析了壓實度和含水率對膨脹率的影響。目前，相關學者對泥巖膨脹變形［7-9］、膨脹時程［10-11］變化規律的研究多是圍繞某單一因素或者多因素定性展開的，而各個因素對泥巖膨脹變形影響程度的研究相對較少，所以研究不同上覆荷載下各因素對膨脹性泥巖的膨脹變形的影響程度具有一定的工程意義。

正交試驗［12］是一種多因素、多水平的試驗設計方法，廣泛應用于工程、科學研究等領域，能夠有效地降低試驗次數，減少試驗誤差，提高試驗效率，并且能夠較好地識別出各個因素的影響程度和因素間的相互作用關系。

基于以上背景，通過室內脹縮試驗，針對河北省張家口市京新高速公路路塹膨脹性泥巖，開展不同含水率、壓實度、上覆荷載條件下的膨脹率試驗，研究膨脹性泥巖浸水膨脹變形特性，分析膨脹速率與試樣浸水時間的關系，以及含水率、壓實度對膨脹率的影響。通過正交試驗設計進行相關試驗，以不同上覆荷載下的膨脹率為判別指標，采用極差分析法對膨脹性泥巖變形特性的影響因素進行敏感性研究，為膨脹性泥巖的工程應用提供科學依據和技術支持。

1 材料與方法

為了解不同壓實度、含水率、浸水時間對膨脹性泥巖的膨脹率的影響，對3種壓實度（85%、90%、95%）、3種含水率（16.2%、18.2%、20.2%）下的膨脹性泥巖進行膨脹率試驗。

1.1 試驗土樣

本次試驗所采用的膨脹性泥巖土樣取自河北省張家口市京新高速公路膠泥灣至西洋河（冀晉界）公路路塹淺層，土樣的顏色呈灰白色，如圖1（a）所示。圖1（b）為環刀試樣。將取回來的膨脹性泥巖試樣碾碎、風干，過2 mm篩，通過輕型擊試驗獲得，泥巖土樣最優含水率為18.2%，最大干密度2.06 g/cm3。測其自由膨脹率為54.5%，屬于弱膨脹巖，泥巖試樣顆粒級配曲線如圖2所示。

1.2 膨脹率試驗

本次試驗主要涉及無荷載膨脹試驗和有荷載膨脹試驗2種。有荷載膨脹試驗是在豎向荷載的作用下，測定土體試樣浸水膨脹后的增加高度與初始高度的比值。通過膨脹率試驗，可以評估膨脹性泥巖在受水浸泡后的膨脹性能，了解其對水分的吸收和膨脹變形的特性。此次試驗操作參照土工試驗方法標準［13］（GB/T 50123—2019）方法，試驗在由南京土壤儀器廠有限公司生產的WG型單杠桿低壓固結儀上進行，如圖3所示。

本次試驗設置的荷載分別為0、12.5、25、50、100 kPa。對于荷載為0 kPa的膨脹性泥巖土樣，裝好試樣后立即向水盒中加水，使水由下向上浸過土樣，使水面超過試樣5 mm，對于荷載不為0 kPa的膨脹性泥巖土樣，先讓其在所加荷載下變形穩定后（每小時記錄一次百分表讀數至變形穩定，即連續2次讀數差值不超過0.01 mm）再向水槽中加水，然后記錄不同時刻膨脹性泥巖發生的膨脹變形，變形穩定的標準為連續2次讀數差值不超過0.01 mm。

膨脹率應按式（1）計算

式中：為某級荷載壓力下的膨脹率，%；為某級荷載壓力下膨脹穩定后的百分表讀數，mm；為試驗開始前百分表初始讀數，mm；為某級荷載壓力下固結儀壓縮變形量，mm；為試樣的初始高度，mm。

1.3 正交試驗安排

本次正交試驗設計的方法以不同荷載作用下的膨脹率為評價指標，3個影響因素為研究對象，考察不同含水率、壓實度及浸水時間對評價指標的影響，每個因素選取3個水平進行正交設計，見表1。假設三因素之間無交互作用， 可選三因素正交表進行安排，根據正交試驗結果進行多因素顯著性分析。本研究使用SPSS數據分析軟件中正交模塊生成L9（34）正交表，見表2。

2 結果與分析

2.1 膨脹性泥巖膨脹率試驗結果分析

根據上述試驗步驟及數據處理方法，對3種不同壓實度、含水率的膨脹性泥巖試樣，分別采用0、12.5、25、50、100 kPa 5級上覆荷載進行膨脹率試驗，探討不同上覆荷載下浸水時間、壓實度、含水率的變化對泥巖土樣的膨脹率的影響。

2.1.1 膨脹性泥巖膨脹變形時程特性分析

考慮到9組試驗結果，觀察到3種不同壓實度、含水率的膨脹性泥巖試樣膨脹率與浸水時間呈現出的規律相似，且由于篇幅限制，以試樣的最優含水率18.2%、最大干密度95%為例展開敘述。不同上覆荷載下的膨脹率與浸水時間的關系曲線如圖4所示。由圖4可知，泥巖試樣的膨脹率隨著上覆荷載的增大而減小，相比0 kPa下的膨脹率，其他4級荷載下試樣的膨脹率大幅度減小，膨脹量減少了65%以上；同時能夠看出，有上覆荷載的泥巖試樣膨脹變形至穩定所需要的浸水時間比無荷載的試樣有所縮短，這是因為上覆荷載能夠抵消土體的部分膨脹潛勢使膨脹穩定時間縮短。

為了更好地描述不同上覆荷載下膨脹性泥巖膨脹率與浸水時間的變化規律，在李志清等［14］研究的對數模型基礎上，使用如式（2）中的三參數Logistic函數來描述膨脹性泥巖膨脹時程曲線有較好的效果，在半對數坐標下不同荷載的泥巖試樣膨脹時程曲線如圖5所示。

式中：S（t）為t時刻泥巖試樣的膨脹率，%；t為浸水時間，min，作為自變量；a、b、p均為相關試驗參數，表征膨脹性泥巖試樣在不同浸水時間情況下吸水膨脹的性能。詳細參數值見表3。

由圖5可知，膨脹性泥巖試樣無論是在無荷載作用下還是有荷載作用下，膨脹率隨浸水時間的變化曲線均呈現S形走勢，大致可劃分為初始膨脹階段、加速膨脹階段、緩慢膨脹階段。這一點與王偉等［15］ 所研究的膨脹土一次浸水后的膨脹時程曲線劃分為3階段規律相似。當浸水時間在0～10 min內，泥巖土樣處于初始膨脹階段，此時水分剛剛浸入試樣土體表層，泥巖試樣開始發生膨脹。相比 0 kPa下泥巖試樣，其他4級荷載下泥巖膨脹潛勢差別不大，但都明顯低于0 kPa下泥巖試樣膨脹潛勢。從浸水時間為10～360 min內，泥巖試樣由于水分慢慢浸入土體內部，開始進入加速膨脹階段，此階段泥巖試樣的膨脹量最大，膨脹速率最快，且在到達第3階段時，泥巖試樣的膨脹量接近達到95%左右。最后到達緩慢膨脹階段（360 min后），這個時候，泥巖試樣已經基本上達到膨脹穩定狀態，膨脹速率趨近于0，理論上講，這種膨脹會隨著浸水時間持續下去，上限是其自由膨脹率。然而對于實際工程而言，緩慢膨脹階段帶來的影響不是很大。

考慮到上述Logistic函數描述膨脹率與浸水時間的關系，對式（2）進行求導，可得

式（3）是關于浸水時間t的函數，是t時刻對應泥巖試樣吸水膨脹過程中膨脹速率［16］，膨脹速率反映了膨脹性泥巖試樣吸水膨脹后膨脹變形的快慢程度，膨脹速率的大小和變化趨勢可以反映出膨脹性泥巖的膨脹特性以及對外界環境的響應能力。分析不同浸水時間下對應的膨脹速率，對研究膨脹性泥巖的吸水膨脹特性具有重要的意義。

圖6給出不同上覆荷載下吸水膨脹過程中，膨脹速率隨浸水時間的演化規律。由圖6可知，從整體上看，膨脹性泥巖試樣在不同荷載下，膨脹速率隨浸水時間先增大后減小慢慢趨近于0。對于無荷載0 kPa下泥巖試樣在整個過程中膨脹速率都明顯高于其他4組有荷載的泥巖試樣，這是因為有上覆荷載的存在，削弱了泥巖試樣的膨脹潛勢，抑制了其吸水膨脹性能。觀察圖6中縮放圖，5種荷載下的泥巖試樣膨脹速率到達峰值點均在膨脹第2階段，0 kPa下的泥巖試樣膨脹速率在60 min后開始減小，反觀其他4組有荷載的泥巖試樣，上覆荷載小的試樣膨脹速率最先達到峰值，荷載大的試樣膨脹速率雖說達到峰值的時間較長，但是數值很小。原因是較小的荷載作用就可以抑制泥巖的大部分膨脹量，最先達到膨脹速率峰值；上覆荷載實際上給泥巖試樣一個壓實效果，荷載越大的試樣壓得越密實，膨脹速率達到峰值的浸水時間越長。因此，可以看出不同上覆荷載下泥巖試樣的浸水時間對試樣的膨脹性能有著重要的影響。

2.1.2 壓實度、含水率對泥巖膨脹的影響

對于膨脹性泥巖來說，壓實度和初始含水率［17］也是非常重要的參數。同一含水率下不同壓實度的膨脹率試驗結果如圖7所示，同一壓實度下不同含水率的膨脹率試驗結果如圖8所示。

由圖7和圖8可知，1）在相同初始含水率條件下，隨著壓實度的增大，有荷膨脹率逐漸增大。這是因為，壓實度大的膨脹性泥巖土樣處于緊密狀態，含有發生膨脹變形的物質基礎多，在相同含水率作用下，預先發生的膨脹變形相差不大，因此含有膨脹物質基礎多的試樣膨脹潛勢更強，發生的膨脹變形更大。2）在相同的壓實度下，隨著含水率的增大，有荷膨脹率整體下降。原因是，膨脹泥巖土樣的壓實度相同，因此具有相同的膨脹物質基礎，而含水率大的土樣在試驗之前比含水率小的土樣發生更多的預膨脹，所以，試驗時含水率大的膨脹性泥巖土樣所能發揮的膨脹變形明顯低于含水率低的土樣。

2.2 正交試驗結果及敏感性分析

2.2.1 試驗結果

按照正交設計方案進行9組試驗，對每組試樣及其平行試樣進行測量，取其平均值作為參考，試驗結果見表4。

2.2.2 敏感性分析

敏感性分析［18］通bRDCJwivpwV0+6p+WNJH7A==常涉及改變每個因素的水平，并觀察結果的變化。通過比較不同水平下的結果，可以確定哪些因素對結果有顯著影響，哪些因素對結果沒有影響。在正交試驗中進行敏感性分析可以幫助確定哪些因素對結果的影響最大，以便在進一步的試驗中重點關注這些因素。這可以幫助節省時間和資源，并提高試驗的效率。

本研究采用極差分析法對正交試驗數據進行分析。極差分析法簡稱R法，簡單易懂、實用性強，應用比較廣泛。極差的大小，反映了因素變化時對試驗指標的影響程度。極差越大，表明該因素的水平改變對試驗指標的影響越大，即該因素的敏感性越大；反之，極差越小，因素的敏感性越小。R法的計算原理與步驟如圖9所示，極差分析結果見表5，影響因素敏感性對比如圖10所示。

通過表5數據以及圖9極差分析可知，在無荷（0 kPa）和12.5 kPa下，3個因素對試樣膨脹率的敏感性大小為C>A>B；對于上覆荷載為25、50、100 kPa而言，3個因素對試樣膨脹率的敏感性大小為C>B>A。其中，在5級荷載下均是浸水時間的極差最大，敏感性最大，為影響膨脹率的主導因素；而次要因素含水率和壓實度之間存在如下關系，當試樣在小的上覆荷載作用下，含水率對試樣膨脹率的影響比較明顯；隨著荷載的增大，壓實度對試樣膨脹率的影響超過了含水率對試樣膨脹率的影響。

通過上述分析，膨脹性泥巖試樣遇水后發生膨脹，膨脹變形影響因素主要是試樣的浸水時間，浸水時間的長短對試樣膨脹變形的潛勢發揮著重要的影響。在實際工程中，可通過加強排水系統的設計和施工，確保排水暢通，縮短土樣的浸水時間，從而減少水分對工程建設的影響，可降低膨脹性泥巖的膨脹和收縮。

3 結論

本研究基于正交試驗理論，用極差分析法對 各影響因素進行敏感性分析，采用張家口市京新高速公路路塹淺層膨脹性泥巖，進行不同初始含水率、壓實度、上覆荷載下的膨脹率試驗，得出以下結論。

1）泥巖試樣的膨脹率隨上覆荷載的增大而減小，對比0 kPa下試樣的膨脹率，其他4級荷載下試樣的膨脹率大幅度減小，膨脹量減少了65%以上。

2）采用三參數Logistic函數來描述膨脹性泥巖膨脹時程曲線具有較好的效果，泥巖試樣無論是在無荷載作用下還是有荷載作用下，膨脹率隨浸水時間的變化曲線均呈現S形走勢，可劃分為初始膨脹階段、加速膨脹階段、緩慢膨脹階段。

3）泥巖試樣在不同荷載下，膨脹速率隨浸水時間表現出先增大后減小慢慢接近于0的趨勢，且膨脹速率到達峰值點均是在膨脹第2階段。0 kPa下泥巖試樣在整個過程中膨脹速率都明顯高于其他4組有荷載的泥巖試樣，原因是有上覆荷載削弱了泥巖試樣的膨脹潛勢。

4）在相同初始含水率條件下，泥巖試樣的膨脹率隨著壓實度的增大而逐漸增大；在相同的壓實度條件下，泥巖試樣的膨脹率隨初始含水率越大而越小。

5）通過正交試驗，采用極差分析法分析各因素對不同上覆荷載下泥巖試樣膨脹率的敏感性大小，分析發現在5級荷載下均是浸水時間的極差最大，敏感性最大，為影響膨脹率的主導因素；含水率和壓實度為影響膨脹率的次要因素，當試樣在小的荷載下，含水率對試樣膨脹率的影響比較明顯；當荷載繼續增大，壓實度對試樣膨脹率的影響超過了含水率對試樣膨脹率的影響。

【參 考 文 獻】

［1］ 陳永貴，劉聰，馬婧，等.鹽堿演化環境下高壓實GMZ膨潤土膨脹變形特性［J］.巖土工程學報，2023，45（4）：690-698

CHEN Y G，LIU C，MA J，et al.Swelling characteristics of compacted GMZ bentonite with saline-alkali evolution solutions［J］.Chinese Journal of Geotechnical Engineering，2023，45（4）：690-698.

［2］ CHERTKOV V Y.Physical modeling of the soil swelling curve vs.the shrinkage curve［J］.Advances in Water Resources，2012，44： 66-84.

［3］ 屈祥.湘西膨脹土膨脹變形及膨脹力特性試驗研究［J］.湖南交通科技，2022，48（4）：40-44.

QU X.Experimental study on swelling deformation and swelling force characteristics of expansive soil in western Hunan［J］.Hunan Communication Science and Technology，2022，48（4）：40-44.

［4］ 張愛軍，哈岸英，駱亞生.壓實膨脹土的膨脹變形規律與計算模式［J］.巖石力學與工程學報，2005（7）：1236-1241.

ZHANG A J，HA A Y，LUO Y S，et al.Swelling deformation and calculation models of compacted expansive soil［J］.Chinese Journal of Rock Mechanics and Engineering，2005（7）：1236-1241.

［5］ 徐永福，吳正根，劉傳新，等.膨脹土的擊實條件與膨脹變形的相關性研究［J］.河海大學學報（自然科學版），1997，25（3）：57-60.

XU Y F，WU Z G，LIU C X，et al.Relativity between compaction conditions and swelling deformationof expansive soils［J］.Journal of Hehai University（Natural Sciences），1997，25（3）：57-60.

［6］ 饒錫保，黃斌，吳云剛，等.膨脹土擊實樣膨脹特性試驗研究［J］.武漢大學學報（工學版），2011，44（2）：211-215.

RAO X B，HUANG B，WU Y G，et al.Experimental study of expansion characteristics of compacted of expansive soil［J］.Journal of Wuhan University （Engineering and Technology），2011，44（2）：211-215.

［7］ ZENG Z，KONG L，WANG M，et al.Assessment of engineering behaviour of an intensely weathered swelling mudstone under full range of seasonal variation and the relationships among measured parameters［J］.Canadian Geotechnical Journal，2018，55（12）：1837-1849.

［8］ 王晨，郭小龍.膨脹性泥巖膨脹性試驗研究［J］.國防交通工程與技術，2023，21（1）：37-40.

WANG C，GUO X L.Experimental study on expansibility of expansive mudstone［J］.Defense Transportation Engineering and Technology，2023，21（1）：37-40.

［9］ 馬福榮.取樣卸荷對膨脹性泥巖強度與變形特性影響的試驗研究［J］.水利水運工程學報，2017（5）：109-116.

MA F R.Experimental study of unloading on impact deformation on deformation behavior and strength of mudstone in sampling state［J］.Hydro-Science and Engineering，2017（5）：109-116.

［10］ CHEN S X.Study of evolution law of cracking for swelling mudstone of Sichuan central Redbeds in hydration process［J］.Japanese Geotechnical Society Special Publication，2019，7（2）：96-101.

［11］ 陳占峰.局部膨脹泥巖套管坍塌理論分析［J］.石油科學與工程學報，2021，203：96-100

CHEN Z F.Theoretical analysis of casing collapse in locally expansive mudstone［J］.Journal of Petroleum Science & Engineering，2021，203：96-100

［12］ 刁心宏，胡樹根，劉宇云，等.膨脹巖相似材料的配比試驗研究［J］.鐵道建筑，2017，57（12）：136-140.

DIAO X H，HU S G，LIU Y Y，et al.Experimental study mix proportion of materials similar to swelling rock［J］.Railway Engineering，2017，57（12）：136-140.

［13］ GB/T 50123—2019土工試驗方法標準［S］.北京：中國計劃出版社，2019.

GB/T 50123—2019 Standard of Geotechnical test method ［S］.Beijing：China Planning Press，2019.

［14］ 李志清，李濤，胡瑞林，等.蒙自重塑膨脹土膨脹變形特性與施工控制研究［J］.巖土工程學報，2008，30（12）：1855-1860.

LI Z Q，LI T，HU R L，et al.Expansion characteristics and construction control of remolded Mengzi expansive soil［J］.Chinese Journal of Geotechnical Engineering，2008，30（12）：1855-1860.

［15］ 王偉，盧廷浩，陶菲菲，等.膨脹土一次性浸水膨脹3參數時程模型［J］.應用基礎與工程科學學報，2008，16（6）：840-845.

WANG W，LU T H，TAO F F，et al.Three-parameter model for expansive deformation process of expansive soil［J］.Journal of Basic Sciences and Engineering，2008，16（6）：840-845.

［16］ 張唐瑜，馬麗娜，張戎令，等.高速鐵路膨脹泥巖路基滲流及膨脹特性試驗研究［J］.水利水運工程學報，2019（2）：48-54.

ZHANG T Y，MA L N，ZHANG R L，et al.Experimental studies on seepage and swell ability of expansive mudstone on high speed railway subgrade［J］.Hydro-Science and Engineering，2019（2）：48-54.

［17］ 李明艷，童軍，徐超.膨脹土的有荷膨脹模型研究及其應用初探［J］.工程地質學報，2023，31（5）：1543-1551.

LI M Y，TONG J，XU C，et al.Preliminary study on loading expansion model of expansive soil and application［J］.Journal of Engineering Geology，2023，31（5）：1543-1551.

［18］ 呂文華.膨脹土邊坡穩定性影響因素的表征參數及其敏感性研究［J］.科技風，2023（25）：96-98.

Lü W H.Study on characterization parameters and sensitivity of influencing factors of expansive soil slope stability［J］.Ke Ji Feng，2023（25）：96-98.