袁 璞，戚少先，位寧寧

(1. 安徽理工大學土木建筑學院，安徽 淮南 232001；2. 安徽理工大學礦山地下工程教育部工程研究中心，安徽 淮南 232001；3. 安徽理工大學深部煤礦采動響應與災害防控國家重點實驗室，安徽 淮南 232001)

隨著社會經(jīng)濟的發(fā)展，地下工程開發(fā)的深度也不斷加大，如各種深部開采的礦山、核廢料深層處理和深部地下防護工程等。深部地下工程是較為熱門的研究領域，干濕循環(huán)作用為水巖作用常見形式之一，其劣化效應是不可忽視的問題[1-2]。

目前，許多學者對干濕循環(huán)作用下巖石損傷特性進行了研究[3-4]。文獻[5]對干濕循環(huán)后的砂巖進行單軸壓縮實驗，對劣化后的力學參數(shù)進行了分析；文獻[6-7]得出了飽水風干后砂巖的縱波波速與孔隙率變化規(guī)律； 文獻[8]用核磁共振技術結(jié)合掃描電子顯微鏡對干濕循環(huán)作用后的砂巖微觀變化進行了研究； 文獻[9]對干濕循環(huán)和地震共同作用下的庫區(qū)泥巖邊坡穩(wěn)定性進行了分析； 文獻[10]發(fā)現(xiàn)干濕循環(huán)對巖石造成的損傷是漸進的并且是不可逆的。以上研究主要集中在干濕循環(huán)對巖石宏觀力學性能的影響，缺乏從微觀角度分析干濕循環(huán)對深部砂巖力學特性劣化機理的探討。

因此，針對深部地下工程中干濕循環(huán)作用引起的巖石劣化效應，本文采用室內(nèi)實驗和細觀數(shù)值模擬相結(jié)合的方法，從宏觀和細觀上共同分析干濕循環(huán)作用對深部砂巖單軸抗壓強度劣化的影響，以期為深部地下工程巖體穩(wěn)定性分析提供一定的實驗依據(jù)。

1 試驗方案設計

1.1 深部砂巖試件

試驗巖樣取自淮南市某煤礦-965m處的深部砂巖，整體呈淺灰色，表面沒有明顯的裂紋，按照規(guī)范[11]將其加工成直徑為50mm、高度為100mm的圓柱形試件(見圖1)。

圖1 深部砂巖試件

1.2 試驗過程

挑選合格的試件分為6組，分別進行0、1、5、10、20、30次的干濕循環(huán)試驗；以試件在60℃烘箱中干燥12h后晾曬至室溫，再放入水中浸泡12h，為1次干濕循環(huán)，用n表示干濕循環(huán)次數(shù)。

2 實驗結(jié)果分析

2.1 干濕循環(huán)后巖石物理特性變化

定義砂巖經(jīng)歷不同干濕循環(huán)次數(shù)之后的孔隙體積與自然狀態(tài)下的砂巖體積之比為孔隙率，即

(1)

式中：P為孔隙率，%；V0為砂巖在自然狀態(tài)下的體積，cm3；V為砂巖的絕對密實體積，cm3；ρ0為水的密度，g/cm3，msat為飽和狀態(tài)下的質(zhì)量，g；mdry為干燥后的質(zhì)量，g。

由圖2可知， 巖石在經(jīng)歷干濕循環(huán)之后孔隙率呈不斷增加的趨勢， 但增加的速率逐漸減弱。 在第10～20次干濕循環(huán)過程中， 砂巖的孔隙率增加最大， 此時宏觀上表現(xiàn)為單軸抗壓強度的明顯降低。

圖2 干濕循環(huán)作用后砂巖孔隙率與縱波波速變化

縱波波速可反映出巖石內(nèi)部的完整性，可定量分析內(nèi)部的損傷程度。由圖2可知，經(jīng)歷干濕循環(huán)作用后，砂巖的縱波波速不斷降低，且在第10次干濕循環(huán)之后砂巖波速降低值最大，此時砂巖內(nèi)部的孔隙率增加也最大，砂巖完整性較差，導致聲波在孔隙處發(fā)生繞射或被吸收，傳播速度降低。

2.2 應力-應變曲線分析

砂巖在不同干濕循環(huán)次數(shù)下的應力-應變曲線如圖3所示。根據(jù)應力-應變曲線可以得出不同循環(huán)次數(shù)下的試件峰值應力、峰值應變和彈性模量(見表1)。

圖3 不同干濕循環(huán)次數(shù)砂巖應力-應變曲線

表1 不同干濕循環(huán)次數(shù)下砂巖力學參數(shù)

從圖3和表1可以看出，峰值應力和彈性模量均隨干濕循環(huán)次數(shù)的增加呈現(xiàn)不同程度的減?。辉诘?0～20次干濕循環(huán)時，強度降低最大，巖石內(nèi)部裂隙擴展最快；在第20～30次干濕循環(huán)時，應力降幅有所減緩。在干濕循環(huán)達到一定的次數(shù)后，應力-應變曲線表現(xiàn)出一定的延性破壞特征，而未經(jīng)過干濕循環(huán)作用為脆性破壞。

3 離散元數(shù)值分析

巖石的損傷是指巖石內(nèi)部的孔隙和微裂紋受到外力作用時繼續(xù)發(fā)展并擴大的過程，最終導致巖石力學參數(shù)的降低。然而這些微小孔隙和裂紋的變化是室內(nèi)實驗難以發(fā)現(xiàn)的，因此對巖石的細觀分析是研究巖石損傷必不可少的手段。本文基于顆粒流軟件PFC2D并結(jié)合室內(nèi)實驗數(shù)據(jù)模擬干濕循環(huán)作用下深部砂巖細觀損傷的演化。

3.1 模型的選擇與建立

PFC2D中的平行粘結(jié)模型可很好地反映砂巖的力鏈分布和裂紋發(fā)展的過程[12]31。

第一步，生成砂巖顆粒。根據(jù)文獻[13]可知，受力后砂巖脆性剪切破裂帶的間距與其平均粒徑成正比，且約為平均粒徑的5～10倍。真實砂巖粒徑在0.04～2.0mm之間，根據(jù)實際情況，模擬砂巖顆粒半徑在0.40～0.60mm之間隨機生成。PFC2D數(shù)值模型寬50mm、高100mm，包含5 964個顆粒。

第二步，確定細觀參數(shù)。文獻[14]提出的由壓縮實驗確定模型細觀參數(shù)的方法，經(jīng)過大量試算，先確定出對峰值應力和峰值應變影響較大的細觀參數(shù)，如平行粘結(jié)的抗拉強度和粘聚力；再根據(jù)文獻[12]77可知，摩擦角的變化對試件破壞形式有較大的影響，剛度比和摩擦系數(shù)對應力-應變曲線特征的影響不大；最后微調(diào)以符合實際情況。確定的細觀參數(shù)如表2所示。

表2 深部砂巖模型細觀參數(shù)表

第三步，施加應力。利用墻體伺服機制，不斷調(diào)整墻體速度達到期望的應力值。

3.2 干濕循環(huán)深部砂巖細觀數(shù)值分析

對經(jīng)歷不同干濕循環(huán)次數(shù)的深部砂巖進行單軸壓縮模擬，其峰值強度如表3所示。

表3 單軸抗壓強度模擬值與實驗值

表3表明PFC2D模擬得到的峰值應力與實驗峰值應力相差較小，最大誤差僅為3.6%。因此，表2的細觀參數(shù)可以很好地分析深部砂巖細觀劣化的過程。

1) 干濕循環(huán)深部砂巖力鏈分布 在外荷載作用下，砂巖顆粒相互接觸在巖石內(nèi)部形成不均勻分布的力鏈網(wǎng)絡，圖4為不同干濕循環(huán)次數(shù)下的力鏈分布圖。力鏈是由3個或多個剛性顆粒接觸形成的線性顆粒組。

圖4 深部砂巖力鏈分布圖

圖4中黑色表示壓力鏈，紅色表示拉力鏈，力鏈的粗細表示接觸力的大小。由于砂巖處于受壓狀態(tài)，此時壓力鏈的數(shù)量和粗細程度明顯高于拉力鏈，當干濕循環(huán)次數(shù)增多時，拉力鏈數(shù)量逐漸減少，壓力鏈數(shù)量增多。

當砂巖試件處于干燥狀態(tài)時，點處砂巖力鏈分布較細且分散，力的分布比較均勻，所以其受壓承載力較大。隨著干濕循環(huán)次數(shù)的增加，力鏈分布逐漸集中，形成強接觸力鏈。此時顆粒間需要承受和傳遞的力逐漸增強，宏觀上表現(xiàn)為砂巖試件更容易劣化。

2)干濕循環(huán)深部砂巖裂紋分布 不同干濕循環(huán)次數(shù)深部砂巖受壓破壞時裂紋分布及演化曲線如圖5和圖6所示。R代表剪切裂紋占總裂紋的比例。

n=0 R=16.8%

圖6 裂紋演化曲線

由圖5可知，隨干濕循環(huán)次數(shù)增加，試件的剪切裂紋數(shù)量逐漸增加，其中紅色表示拉伸裂紋，藍色表示剪切裂紋。干濕循環(huán)次數(shù)較少時，砂巖破壞主要以拉伸裂紋為主，表現(xiàn)為拉裂破壞。當干濕循環(huán)次數(shù)達到10次時，砂巖破壞時其內(nèi)部剪切裂紋明顯增多，破壞形式開始向剪切破壞發(fā)展，由于此時拉伸裂紋還占有一定的比例，所以最終破壞形式表現(xiàn)為脆性剪切破壞。

由圖6可知，在加載初期，試件并沒有新增裂紋，此時處于壓密階段。隨著軸向變形逐漸增大，新生裂紋逐漸出現(xiàn)并加速發(fā)展，大部分裂紋是在試件破壞前短時間內(nèi)形成的，當裂紋逐漸積累并貫通時，試件產(chǎn)生宏觀裂縫并喪失承載力。

4 干濕循環(huán)砂巖微觀結(jié)構(gòu)分析

為研究干濕循環(huán)次數(shù)對砂巖內(nèi)部結(jié)構(gòu)劣化的影響，從不同干濕循環(huán)次數(shù)的砂巖中取出若干試樣，先用乙醇沖洗，再放入烘箱干燥，最后采用上海百賀儀器科技有限公司的FEI掃描電鏡進行觀測。

由圖7可知，在干燥狀態(tài)下，砂巖顆粒間排列緊密、孔隙少，因此強度較高。干濕循環(huán)作用后，顆粒間逐漸出現(xiàn)細小的孔隙，微裂紋逐漸開始發(fā)育。隨干濕循環(huán)次數(shù)的不斷增加，鈣質(zhì)膠結(jié)物溶解、溶蝕以及水巖相互作用，砂巖的微細觀孔隙、裂隙逐漸發(fā)育、貫通，尺度逐漸增大，顆粒結(jié)構(gòu)逐漸松散[15]，砂巖由較致密的孔隙膠結(jié)結(jié)構(gòu)逐漸向不規(guī)則的蜂窩狀結(jié)構(gòu)發(fā)展，進而導致了砂巖宏觀變形和強度的劣化。

圖7 干濕循環(huán)作用后深部砂巖SEM圖(1 000倍)

5 結(jié)語

隨著干濕循環(huán)次數(shù)的增加，砂巖的孔隙率、縱波波速以及單軸抗壓強度均逐漸降低，砂巖微觀結(jié)構(gòu)從致密狀發(fā)展為蜂窩狀。在干濕循環(huán)作用初期，砂巖內(nèi)部力鏈分布較為分散，沒有形成強接觸力鏈，此時承載能力較高。在干濕循環(huán)10次以后，剪切裂紋逐漸增多，開始向剪切破壞發(fā)展。

在實際工程中，除干濕循環(huán)作用外，深部砂巖還可能受到化學腐蝕與溫度變化等作用，干濕循環(huán)與多種環(huán)境耦合作用下深部砂巖的力學特性還需進一步研究。