王軍保 ，劉新榮 ，楊欣，黃明

(1. 西安建筑科技大學 土木工程學院，陜西 西安，710055；2. 重慶大學 土木工程學院，重慶，400045；3. 福州大學 土木工程學院，福建 福州，350108)

鹽巖是一種特殊的軟巖，具有結構致密、孔隙度低、滲透性低、力學性質穩(wěn)定、蠕變性良好和損傷自恢復能力強等諸多優(yōu)點，被公認為是能源(石油、天然氣)地下儲存和放射性廢棄物地下處置的理想介質[1]。蠕變特性作為鹽巖最為重要的力學性質之一，是影響鹽巖地下儲存庫長期穩(wěn)定性、安全性以及可用性的關鍵因素：因此，對鹽巖蠕變特性進行研究尤為重要。由于鹽巖在能源地下儲存方面的重要地位，近年來國內外學者對鹽巖蠕變本構關系進行了較多研究，提出了多種蠕變模型[1-9]。但總體上，這些模型可以分為以下三大類[10]：經驗模型、唯象模型和元件組合模型。從整個巖石力學流變領域來看，元件組合模型由于本構方程形式簡單，模型及其參數物理意義明確，且能夠把巖石復雜的力學性質直觀地表現出來，因而應用較為廣泛。但傳統(tǒng)的元件組合模型是一種線性模型，無法描述巖石蠕變的非線性特征。為了克服這一不足，一種有效的方法就是采用非線性流變元件代替常規(guī)線性流變元件，并據此建立巖石蠕變模型。熊良宵等[11]提出了一個可以描述巖石非線性黏彈性流變特征的黏滯體，將該非線性黏滯體替換Burgers 模型中的線性黏滯體，得到了非定常 Burgers 模型。黃明等[12]提出了一種修正Kelvin 模型，并將該模型與黏滯體、彈性體進行串聯(lián)建立了可描述巖石非線性蠕變特性的修正Burgers 模型?？涤绖偟萚13]為了描述巖石蠕變全過程，用非定常黏滯體替換Burgers 模型中的定常黏滯體，給出了一種非定常Burgers 模型。Steipi 等[14]將西原模型中的黏塑性參數描述為黏塑性應變的函數，提出了改進的西原模型。蔣昱州等[15]提出了一個變系數的黏滯體，建立了新的巖石非線性黏彈塑性全過程蠕變模型。楊文東等[16]構建了由彈性體、村山體、非線性黏塑性體串聯(lián)而成的非線性黏彈塑性流變模型，該模型可以描述巖石流變的加速階段。曹平等[17]定義應力與試件長期強度的比值為加速蠕變速率冪級數n，模型發(fā)生加速蠕變時的總蠕變量為蠕變特征長度 εc，進而建立了能夠描述巖石加速蠕變的力學模型。以往對于巖石非線性蠕變模型的研究大多是針對巖石加速蠕變階段展開的，而巖石特別是鹽巖，在前2 個蠕變階段也存在非線性特征，且在圍壓作用下，鹽巖加速蠕變并不容易出現[1]?；诖?，本文作者結合非線性流變力學理論，對常規(guī)黏滯體元件進行改進，提出一種變系數非線性黏滯體，然后通過將該非線性黏滯體替換常規(guī)Burgers 模型中的線性黏滯體，建立一種非線性Burgers 模型即NBurgers 模型，并利用鹽巖三軸壓縮蠕變試驗結果對NBurgers 模型參數進行反演識別，以驗證模型的適用性。

1 非線性黏滯體的提出及其特性分析

傳統(tǒng)的元件組合模型通常無法反映巖石蠕變的非線性特征，這是因為這類模型將黏性元件假定為理想的牛頓流體，認為其黏性系數是固定不變的常數的緣故[15]。近年來，隨著人們對巖石非線性蠕變特征的逐步研究和認識，發(fā)現這一假定不夠準確，巖石材料的黏滯系數在蠕變過程中是不斷發(fā)生變化的。陳文玲等[18]通過對巖石蠕變過程的分析，指出當施加的荷載為中高應力水平但未達到巖石破壞應力時，巖石蠕變僅出現前2 個階段，即衰減蠕變階段和穩(wěn)態(tài)蠕變階段，此時，黏滯系數隨著時間增加而逐漸增大，并最終達到一個穩(wěn)定值后不再變化。

受此啟發(fā)，為了更好地描述鹽巖的非線性蠕變特征，本文提出一個變系數的非線性黏滯體，其本構關系如下：

式中：σ 為應力；dε/dt為應變速率；t 為時間；m 為非線性系數；η0為黏滯系數的最大值。

從式(1)可以看出該黏滯體的黏滯系數隨時間增加不斷發(fā)生變化，且其變化規(guī)律符合下式：

下面對式(2)函數特點進行分析。

當t=0 時，

即 η (t)初始值不為0。

將式(2)對時間求導數并整理，可得

說明η (t)單調遞增。

將式(2)變形可得

由式(5)可知：當t→∞時，η (t)→η0，這說明η (t)存在上限值且上限值為η0。

由以上分析可以看出：當時間t 從0→∞時，黏滯系數η (t)從η0/(1+m)→η0。因此，式(2)所示函數特點與文獻[18]所描述的巖石黏滯系數隨時間的變化規(guī)律較符合。

式(1)變形可得

將式(6)對時間積分，并結合初始條件(t=0，ε=0)，可得該非線性黏滯體的蠕變方程為

1.1 非線性系數對蠕變應變的影響

圖1 所示為σ和η0均相同的情況下，非線性系數m 的變化對蠕變曲線的影響規(guī)律示意圖。由圖1可以看出：當m=0時，此時實質為常規(guī)黏滯體，蠕變曲線為斜率不變的直線，僅存在穩(wěn)態(tài)蠕變階段；當m≠0時，蠕變曲線不再是直線，而是一簇形狀相似的曲線簇，且均存在衰減蠕變和近似穩(wěn)態(tài)蠕變2 個階段。隨著m 增大，衰減蠕變階段越來越明顯，持續(xù)時間越來越長，進入穩(wěn)態(tài)蠕變的時間越來越晚，且相同時刻的變形量也越來越大。這說明非線性黏滯體具有衰減蠕變特性，與常規(guī)黏滯體不同。

圖1 不同非線性系數下的蠕變曲線Fig.1 Creep curves with different nonlinear coefficients

1.2 非線性系數對蠕變速率的影響

式(6)所示為非線性黏滯體蠕變速率的表達式。由式(6)可知：m=0時，此時實質為常規(guī)黏滯體，蠕變速率保持σ/η0不變；m≠0時蠕變速率隨時間增加而單調遞減，且當t=0時，dε/dt=σ(1+m)/η0；當t→∞時，dε/dt=σ/η0。這說明非線性黏滯體的初始蠕變速率為常規(guī)黏滯體蠕變速率的(1 + m)倍，且只有當t→∞時，2 種黏滯體的蠕變速率才相等，其余時刻非線性黏滯體的蠕變速率均大于常規(guī)黏滯體的蠕變速率。

圖2 所示為不同非線性系數下的蠕變速率曲線示意圖。由圖2 可見：當m=0時，蠕變速率最小且保持恒定，蠕變速率曲線為1 條水平線；當m≠0時，蠕變速率在初始時刻最大，隨著時間延長而逐漸衰減，且m 越大衰減過程越顯著。理論上，只有當t→∞時，非線性黏滯體和常規(guī)黏滯體的蠕變速率才會相等；而實際上，由于非線性系數m 是1 個有限值，當時間增加到一定程度后，m/(t+1)已經接近于0，此時兩者的蠕變速率已經非常接近。圖2 所示曲線可以清楚地說明這一點。

圖2 不同非線性系數下的蠕變速率曲線Fig.2 Creep rate curves with different nonlinear coefficients

由以上分析可以看出：非線性黏滯體不僅包含了常規(guī)黏滯體的基本特性，而且還具有常規(guī)黏滯體所不具有的特殊性質。因此，在描述鹽巖蠕變特性時將會具有更好的適應性。

2 鹽巖蠕變試驗及非線性蠕變模型

2.1 鹽巖蠕變試驗

試驗所用鹽巖試樣取自江蘇淮安某鹽礦，主要成分為NaCl，質地較純，呈白色、灰白色，部分略帶灰黑色不溶物雜質。按照巖石力學試驗標準的要求，將巖樣加工成直徑為50 mm、高度為100 mm 的圓柱形標準試件。試驗儀器采用長春朝陽儀器廠生產的RLW-2000 巖石流變試驗機。該設備主要由機架、軸向穩(wěn)壓系統(tǒng)、側向穩(wěn)壓系統(tǒng)、數字控制系統(tǒng)及微機系統(tǒng)等5 部分組成，最大軸向荷載為2 000 kN，最大圍壓為60 MPa。

為了研究該礦鹽巖的蠕變特性，對4 個試件進行不同圍壓下的分級加載蠕變試驗。圖3 所示為圍壓10 MPa 試件的分級加載蠕變試驗曲線，軸向應力分別為20，25 和30 MPa。由圖3 可以看到：鹽巖在各級荷載作用下的變形均由加載過程產生的瞬時應變、衰減蠕變應變以及蠕變速率較為穩(wěn)定的穩(wěn)態(tài)蠕變應變3部分組成。由于試驗所施加的最高應力未達到使巖樣發(fā)生加速蠕變的臨界應力，且受試驗條件限制，加載時間不夠長，鹽巖未進入到加速蠕變階段；在圍壓不變的情況下，隨著軸向應力增大，鹽巖瞬時應變、蠕變應變、穩(wěn)態(tài)蠕變速率等均隨之增大，且軸向應力越大，衰減蠕變階段蠕變曲線的曲率半徑越大，經歷的時間越長，達到穩(wěn)態(tài)蠕變的時間越晚。

圖3 鹽巖分級加載蠕變試驗曲線Fig.3 Creep curve of salt rock under step loading

2.2 非線性蠕變模型的建立

通過前面對鹽巖蠕變曲線特征的分析可以看出，試驗曲線與Bugrers 模型蠕變曲線較為相似，但常規(guī)Bugrers 模型無法反映鹽巖蠕變的非線性特征。為了能夠較好地描述鹽巖的非線性黏彈性蠕變特性，本文將前面提出的非線性黏滯體替換常規(guī)Bugrers 模型中的線性黏滯體，從而構成變系數的非線性Burgers 模型，本文將其命名為NBurgers 模型。NBurgers 模型由彈性體(H)、非線性黏滯體(NN)和非線性黏彈性體(NN/H)共3 部分串聯(lián)組成?？紤]到參數數量過多會給模型參數求解增加難度，本文取NN 體的非線性系數m1和NN/H 體的非線性系數m2為同一值，且均以m 來表示。模型示意圖見圖4。

圖4 NBurgers 模型示意圖Fig.4 Schematic view of NBurgers model

在該模型中，假定彈性體的變形為 εe，非線性黏滯體的變形為 εnv，非線性黏彈性體的變形為 εnve，則模型總變形可表示為

假定在t=0時，施加瞬間應力σ=σ0，則有

(1) 彈性體(H)的本構關系為

式中：E1為彈性模量。

(2) 非線性黏滯體(NN)中，結合式(7)可得其蠕變方程為

(3) 非線性黏彈性體(NN/H)中，結合常規(guī)Kelvin模型，可得其本構方程為

式中：E2為黏彈性模量。

將式(11)變形為

對式(12)解微分方程，可得

式中：C 為積分常數。

對于非線性黏彈性體(NN/H)來說，當t=0時，εnve=0。對式(13)代入初始條件，則可解得

將C 代入式(13)可得非線性黏彈性體的蠕變方程為

由于模型中3 部分串聯(lián)，根據式(8)，NBurgers 模型在一維應力狀態(tài)下的蠕變方程可寫為

在三維應力狀態(tài)下，鹽巖內部應力張量 σij可以分解為應力球張量 σm和應力偏張量 Sij，應變張量 εij分解為應變球張量 εm和應變偏張量 eij。根據廣義虎克定律，則有

其中

式中：K 為體積模量；G 為剪切模量；ν 為泊松比。

假定鹽巖為各向同性材料，其蠕變變形僅由應力偏張量引起，且在蠕變過程中鹽巖泊松比不隨時間變化[16]，則在三維應力狀態(tài)下，考慮σ2=σ3，NBurgers模型的軸向蠕變方程為

在三維應力狀態(tài)下，當σ2=σ3時，Burgers 模型的軸向蠕變方程為[11]

對比式(19)和式(20)可以看出：當非線性系數m=0時，NBurgers 模型退化為常規(guī)Burgers 模型。因此，常規(guī)Burgers 模型是本文非線性模型在m=0時的一種特例。

3 鹽巖非線性蠕變模型參數識別

為了對NBurgers 模型進行參數識別，首先將圖3所示鹽巖分級加載蠕變曲線進行處理，得到了不同軸壓時的單級蠕變曲線(圖5 中散點曲線)。

最小二乘法是確定巖石蠕變參數應用較為廣泛的一種方法。本文根據鹽巖蠕變試驗結果，基于非線性最小二乘法基本原理，利用1stOpt 軟件，采用自定義函數的方法，對式(19)所示NBurgers 模型參數進行反演識別。具體方法及過程如下。

(1) 以待反演的蠕變參數作為設計變量，即

圖5 擬合曲線和試驗曲線對比圖Fig.5 Comparison between fitting curves and test curves

式中：K 和G1可根據式(18)轉變?yōu)榇_定彈性模量E1，鹽巖泊松比取ν=0.3。

(2) 取設計變量如式(21)所示，建立目標函數

式中：N 為試驗數據組數；wi(X,ti)為t 時刻計算變形值；wi為t 時刻試驗實測變形值。

(3) 設定目標函數Q 的控制精度并進行參數迭代求解。若Q 滿足精度要求，則停止迭代，輸出計算結果；若不滿足，則繼續(xù)迭代，直到滿足精度要求為止。

參數反演結果見表1。圖5 給出了NBurgers 模型擬合曲線與試驗曲線的對比情況。圖5 同時給出了采用常規(guī)Burgers 模型對試驗數據的擬合情況。

從表1 和圖5 可以看到：NBurgers 模型擬合效果良好，相對誤差較小，擬合值非常逼近于試驗值，且其與試驗值的吻合程度要明顯比常規(guī)Burgers 模型的高，說明NBurgers 模型能夠更好地描述鹽巖的黏彈性蠕變特性，具有一定的適用性。

表1 NBurgers 模型參數反演結果Table 1 Creep parameters of NBurgers model

4 結論

(1) 基于非線性流變力學理論，對常規(guī)黏滯體元件進行了改進，提出了一種變系數的非線性黏滯體，并對其特性進行了分析。分析表明，由于非線性系數的引入，改進后的非線性黏滯體在描述巖石蠕變特性時具有更好的適應性。

(2) 將非線性黏滯體替換常規(guī)Burgers 模型中的線性黏滯體，建立了非線性Burgers 模型，即NBurgers模型。該模型在一定條件下可退化為常規(guī)Burgers模型。

(3)對鹽巖試件開展了三軸壓縮蠕變試驗，并利用試驗結果對模型參數進行了反演識別。擬合曲線和試驗曲線對比顯示，NBurgers 模型擬合曲線和試驗曲線的吻合程度要明顯高于常規(guī)Burgers 模型，說明該模型能夠更好地描述鹽巖的黏彈性蠕變特性。

[1] 杜超, 楊春和, 馬洪嶺, 等. 深部鹽巖蠕變特性研究[J]. 巖土力學, 2012, 33(8): 2451-2456.DU Chao, YANG Chunhe, MA Hongling, et al. Study of creep characteristics of deep rock salt[J]. Rock and Soil Mechanics,2012, 33(8): 2451-2456.

[2] Chan K S, Brodsky N S, Fossum A F, et al. Damage-induced non-associated inelastic flow in rock salt[J]. International Journal of Plasticity, 1994, 10(6): 623-642.

[3] HOU Zhengmeng. Mechanical and hydraulic behaviour of salt in the excavation disturbedzone around underground facilities[J].International Journal of Rock Mechanic and Mining Science,2003, 40(5): 725-738

[4] Jin J, Cristescu N D. An elastic/viscoplastic model for transient creep of rock salt[J]. International Journal of Plastic, 1998,14(1/2/3): 85-107.

[5] WANG Guijun. A new constitutive creep-damage model for salt rock and its characteristics[J]. International Journal of Rock Mechanic and Mining Science, 2004, 41(Supp 1): 61-67.

[6] ZHOU Hongwei, WANG Chunping, Han B B, et al. A creep constitutive model for salt rock based on fractional derivatives[J].International Journal of Rock Mechanics and Mining Sciences,2011, 48(1): 116-121.

[7] 李萍, 鄧金根, 孫焱, 等. 川東氣田鹽巖、膏鹽巖蠕變特性試驗研究[J]. 巖土力學, 2012, 33(2): 444-448.LI Ping, DENG Jingen, SUN Yan, et al. Experimental study of creep characteristics of rock gypsum-salt and rock salt in Chuandong gas field[J]. Rock and Soil Mechanics, 2012, 33(2):444-448.

[8] 郤保平, 趙陽升, 趙延林, 等. 含高鹽份泥巖夾層的鹽巖蠕變特性研究[J]. 地下空間與工程學報, 2007, 3(1): 23-26.XI Baoping, ZHAO Yangsheng, ZHAO Yanlin, et al. Study on creep property of rock salt with mudstone interlayer[J]. Chinese Journal of Underground Space and Engineering, 2007, 3(1):23-26.

[9] 劉江, 楊春和, 吳文, 等. 鹽巖蠕變特性和本構關系研究[J].巖土力學, 2006, 27(8): 1267-1271.LIU Jiang, YANG Chunhe, WU Wen, et al. Study on creep characteristics and constitutive relation of rock salt[J]. Rock and Soil Mechanics, 2006, 27(8): 1267-1271.

[10] 吳文, 楊春和, 侯正猛. 鹽巖中能源(石油和天然氣)地下儲存力學問題研究現狀及其發(fā)展[J]. 巖石力學與工程學報, 2005,24(增2): 5561-5568.WU Wen, YANG Chunhe, HOU Zhengmeng. Investigation on studied situations associated with mechanical aspects and development for underground storage of petroleum and natural gas in rock salt[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(Supp 2): 5561-5568.

[11] 熊良宵, 楊林德, 張堯. 巖石的非定常Burgers 模型[J]. 中南大學學報(自然科學版), 2010, 41(2): 679-684.XIONG Liangxiao, YANG Linde, ZHANG Yao. Non-stationary Burgers model for rock[J]. Journal of Central South University(Science and Technology), 2010, 41(2): 679-684.

[12] 黃明, 劉新榮. 三峽庫區(qū)泥質粉砂巖非線性蠕變模型辨識[J].巖石力學與工程學報, 2011, 30(增2): 3804-3810.HUANG Ming, LIU Xinrong. Research on nonlinear creep model identification of siltite in three Gorges reservoir area[J].Chinese Journal of Rock Mechanics and Engineering, 2011,30(Supp 2): 3804-3810.

[13] 康永剛, 張秀娥. 基于Burgers 模型的巖石非定常蠕變模型[J].巖土力學, 2011, 32(增1): 424-427.KANG Yonggang, ZHANG Xiue. A non-stationary model for rock creep based on Burgers model[J]. Rock and Soil Mechanics,2011, 32(Supp 1): 424-427.

[14] Steipi D, Gioda G. Visco-plastic behaviour around advancing tunnels in squeezing rock[J]. Rock Mechanics and Rock Engineering, 2009, 42(2): 319-339.

[15] 蔣昱州, 張明鳴, 李良權. 巖石非線性黏彈塑性蠕變模型研究及其參數識別[J]. 巖石力學與工程學報, 2008, 27(4):832-839.JIANG Yuzhou, ZHANG Mingming，LI Liangquan. Study on nonlinear viscoelasto-plastic creep model of rock and its parameter identification[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(4): 832-839.

[16] 楊文東, 張強勇, 陳芳, 等. 輝綠巖非線性流變模型及蠕變加載歷史的處理方法研究[J]. 巖石力學與工程學報, 2011, 30(7):1405-1413.YANG Wendong, ZHANG Qiangyong, CHEN Fang, et al.Research on nonlinear rheological model of diabase and treatment for creep loading history[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(7): 1405-1413.

[17] 曹平, 劉業(yè)科, 蒲成志, 等. 一種改進的巖石黏彈塑性加速蠕變力學模型[J]. 中南大學學報(自然科學版), 2011, 42(1):142-146.CAO Ping, LIU Yeke, PU Chengzhi, et al. An improved accelerated creep mechanical model of viscoelasto-plastic rock[J]. Journal of Central South University (Science and Technology), 2011, 42(1): 142-146.

[18] 陳文玲, 趙法鎖, 弓虎軍. 三軸蠕變試驗中云母石英片巖蠕變參數的研究[J]. 巖石力學與工程學報, 2011, 30(增 1):2810-2816.CHEN Wenling, ZHAO Fasuo, GONG Hujun. Study of creep parameters of mica-quartzoseschist during triaxial creep test[J].Chinese Journal of Rock Mechanics and Engineering, 2011,30(Supp 1): 2810-2816.