喻 超，黃小蘭，，馬洪嶺

(1.武漢工業(yè)學院土木工程與建筑學院，湖北武漢430023;2.中國科學院武漢巖土力學研究所巖土力學與工程國家重點試驗室，湖北武漢430071)

鹽巖結(jié)構(gòu)致密，孔隙率低，滲透一般小于10-20m2，力學方面具有良好的蠕變性可以保證洞庫的密封性，損傷與自我恢復(fù)的穩(wěn)定力學性能可適應(yīng)儲存壓力的變化，被世界各國公認為能源儲存和高放射核廢料處置的理想介質(zhì)。西方發(fā)達國家都已建起大量鹽巖地下儲氣(油)庫，用于國家戰(zhàn)略能源儲備和商業(yè)油氣儲備。我國地下鹽巖資源十分豐富，分布范圍廣，埋藏于地下數(shù)十米至4000米，具有良好的建設(shè)地下能源儲存或廢棄物處置庫的地質(zhì)條件。但與發(fā)達國家巨鹽丘構(gòu)造的鹽巖礦床相比，我國鹽巖層具有層數(shù)多、雜質(zhì)多的特點，鹽巖為層狀分布，內(nèi)含多個泥巖、泥質(zhì)硬石膏夾層［1］。由于與國外鹽巖礦床地區(qū)條件不一致，因此，在我國特殊地層中建腔進行長期能源儲存面臨更加復(fù)雜的力學和技術(shù)問題。

互層鹽巖體存在層狀不均勻巖性問題和不同巖層之間差異變形問題［2］，鹽巖與夾層間存在的力學特性差異及變形不協(xié)調(diào)，會造成兩種巖性交界處的剪切破壞，影響儲庫安全性［3］。蠕變試驗研究認為，層狀鹽巖蠕變過程中泥巖層對鹽巖層的蠕變有一定抑制作用，在鹽巖層及夾層交界面上將因蠕變率不匹配而產(chǎn)生變形不協(xié)調(diào)現(xiàn)象，從而影響油氣儲庫的致密、長期穩(wěn)定性和安全性［4］。本文在前期學者試驗研究和理論分析基礎(chǔ)上，從數(shù)值模擬角度，對含不同夾層厚度鹽巖儲庫的穩(wěn)定性進行分析，考慮不同恒壓長期蠕變情況下，探究腔體內(nèi)壓及夾層厚度對腔體收縮率、塑性區(qū)發(fā)展等的影響。

1 鹽巖儲庫穩(wěn)定性評判標準

通常地下洞室的穩(wěn)定性準則主要以相對收斂變形、收斂變形速率、塑性區(qū)范圍大小和拱頂是否發(fā)生較大塌落等標準進行判斷。如淺埋大跨度的地下工程一般以塑性區(qū)面積、最大塑性區(qū)半徑、拉裂破壞區(qū)面積以及頂部相對下沉量為標準;埋深大、圍巖完好的隧道一般以最大允許位移為評判標準。

鹽巖作為一種特殊的巖體，其層狀復(fù)合體表現(xiàn)出與普通層狀巖體不一樣的破壞狀態(tài)，且采用水溶開挖形成的洞室不可能采用傳統(tǒng)的支護結(jié)構(gòu)增加洞室的穩(wěn)定性及限制圍巖的位移，其穩(wěn)定性評判準則與其他巖體有一定的差別。目前對鹽巖地下儲庫的穩(wěn)定性評判內(nèi)容包括［5-7］：局部穩(wěn)定性，整體穩(wěn)定性及可使用性。整體穩(wěn)定性主要指儲庫群腔間礦柱的穩(wěn)定性，是保障儲庫正常運營生產(chǎn)的關(guān)鍵。礦柱設(shè)留太大，會浪費大量鹽巖資源，影響經(jīng)濟效果。如果設(shè)留太小，則會突跳失穩(wěn)破壞，從而導(dǎo)致頂板大面積破壞，致使群腔報廢。局部穩(wěn)定性評判標準則包括無拉應(yīng)力準則，以及鹽巖的剪脹破壞準則?？墒褂眯园ㄈ齻€方面：腔體體積收縮率，地面沉降以及儲庫的密封性。運行期間由于鹽巖的蠕變使得體積發(fā)生收斂減少影響其調(diào)峰能力，甚至不能正常使用，因此腔體的體積收縮率要保證小于某個值。儲庫在儲存壓力的變化和上覆巖層的壓力作用下發(fā)生變形，將會引起地面沉降，進而引起地面建筑物的損壞，因此地面沉降亦應(yīng)控制在一定范圍之內(nèi)。密封性是保證運行期間儲庫中的油、氣不泄漏。雖然鹽巖自身具有優(yōu)越的力學性能，但鹽巖與夾層間存在的力學特性差異及變形不協(xié)調(diào)，在長期交替壓力運行條件下會造成兩種巖性交界處的剪切破壞，進而影響油氣儲庫的致密、長期穩(wěn)定性和安全性。

目前雖對鹽巖儲庫穩(wěn)定性評判有了一些研究，但在各國還沒有統(tǒng)一的標準和設(shè)計規(guī)范，多數(shù)做法是根據(jù)具體的儲庫及其巖體力學特性，預(yù)先設(shè)置儲庫穩(wěn)定性的一些標準，然后通過試驗研究獲得所需參數(shù)，最后通過數(shù)值計算確定儲庫的穩(wěn)定性。

2 層狀鹽巖儲庫穩(wěn)定性模擬計算分析

2.1 計算模型及參數(shù)選取

為比較分析不同夾層厚度對鹽巖腔體穩(wěn)定性的影響，筆者建立了圓柱形儲庫的單腔三維模型，模型以夾層的中心為坐標原點(距地表1000 m，即腔體中心埋深)，夾層厚度分別取為2 m，4 m，6 m三種工況，腔體高度為60 m，腔體直徑30 m，腔體上部鹽巖厚度120 m，下部鹽巖厚度120 m，腔外巖體直徑200 m?？紤]到問題的軸對稱性，數(shù)值模擬過程中選取四分之一腔體作為研究對象，側(cè)面約束其相應(yīng)的橫向位移，底面約束其豎向位移，頂面施加均布荷載，荷載的大小等于模型上部巖體的自重。所建模型的計算網(wǎng)格劃分情況如圖1所示。

圖1 計算模型網(wǎng)格劃分

靜力計算選用Mohr-Coulomb模型。取層狀鹽巖體的地應(yīng)力為靜水壓狀態(tài)，在約束相應(yīng)邊界并施加上覆巖層靜壓力之后，求解初始應(yīng)力場。隨后對腔壁及夾層表面施加靜壓力，對腔體部分進行開挖。腔體開挖后，取腔內(nèi)氣體壓力分別為5 MPa，10 MPa，15 MPa，20 MPa四種工況長期蠕變。鹽巖和夾層的蠕變模型，采用冪函數(shù)指數(shù)模型。

A，n—材料特性;

結(jié)合相關(guān)地質(zhì)資料，鹽巖密度取2.30 g/cm3，夾層密度取2.50 g/cm3，重力加速度取9.8 N/kg。鹽巖和泥巖夾層力學參數(shù)及蠕變參數(shù)選取如表1所示。

表1 數(shù)值計算力學參數(shù)［6］

2.2 腔體收縮率計算結(jié)果分析

研究工作表明鹽巖即使在非常小的偏應(yīng)力作用下，也具有流變特性;鹽巖的蠕變速率是偏應(yīng)力和溫度的高階非線性函數(shù)［8-12］。即只要鹽巖儲存庫應(yīng)力存在差值，鹽巖溶腔體積減少就不可避免。圖2給出了不同內(nèi)壓工況條件下腔體含不同夾層厚度時流變20年體積收縮率對比情況。

圖2 不同內(nèi)壓下含不同夾層厚度腔體蠕變收縮率對比

從圖2對比曲線發(fā)現(xiàn)：恒定內(nèi)壓下，隨著夾層厚度的增大，腔體收縮率呈略微上漲，影響并不顯著，且隨內(nèi)壓的增大，影響值愈小。蠕變20年，夾層厚度對腔體收縮率的影響最大值發(fā)生在內(nèi)壓為5 MPa時，腔體收縮率分別為45.57%(夾層厚2 m)，46.27%(夾層厚4 m)，46.96%(夾層厚6 m)，最大差值1.39%;當內(nèi)壓為20 MPa時，腔體收縮率分別為0.412%(夾層厚 2 m)，0.424%(夾層厚 4 m)，0.435%(夾層厚6 m)，最大差值不足0.03%。

考慮到相同內(nèi)壓條件下夾層厚度對腔體收縮率的影響并不太，現(xiàn)選取夾層厚度為4 m情況時，對腔體收縮率隨內(nèi)壓的變化情況進行分析。如圖3及表2所示。

圖3 不同恒定內(nèi)壓下腔體體積收縮率隨時間變化曲線

表2 不同內(nèi)壓不同時間腔體體積收縮率 /%

從圖3及表2可以看出，在恒定內(nèi)壓作用下，溶腔體積收縮率隨時間而增加，但增長率呈下降趨勢;隨著內(nèi)壓增加，相同時間段內(nèi)腔體收縮率顯著減少，當內(nèi)壓為20 MPa時，腔體甚至在前1.5年時間內(nèi)出現(xiàn)了擴張，20年流變腔體收縮率僅為0.42%(當內(nèi)壓為5 MPa時，腔體收縮率高達46.27%)，可見增加內(nèi)壓對腔體的體積減少有明顯的抑制作用。從保持腔體體積考慮，在儲氣庫運行期間，維持較高的運行壓力，盡可能降低低壓下運行時間，對腔體的可用性是有利的。

體積收縮率是保證溶腔儲庫正常運營的關(guān)鍵指標之一。儲庫在運行期間，由于鹽巖的蠕變使得體積發(fā)生收斂，進而影響儲庫的正常使用。目前通常規(guī)定儲庫運行5年后腔體體積收縮率控制在5%范圍之內(nèi)作為儲庫穩(wěn)定性評價標準［5，6］。由表 2可知，當內(nèi)壓大于10 MPa時，從收縮率角度評價，腔體滿足穩(wěn)定性要求。

2.3 塑性區(qū)及位移矢量計算和分析

根據(jù)計算結(jié)果，隨夾層厚度增加，夾層位置處變形呈增大趨勢，但隨內(nèi)壓的增加，這種變形增大趨勢明顯減小，內(nèi)壓達到15 MPa時，影響值已不明顯;圍巖塑性區(qū)對夾層厚度的變化反應(yīng)并不明顯?，F(xiàn)以夾層厚度為4 m情況進行分析，如圖4所示。

圖4 不同內(nèi)壓時20年蠕變圍巖塑性區(qū)及位移矢量分布

由4圖可以發(fā)現(xiàn)，隨著溶腔內(nèi)壓的增加，腔周塑性區(qū)及位移量明顯減小。當腔體內(nèi)壓達到或超過10 MPa時，溶腔圍巖已無明顯破損區(qū)，說明溶腔內(nèi)壓大小對腔體的穩(wěn)定性有較大影響。從位移矢量分布圖可以看出，較大的位移主要集中在腔體開挖面上，向腔內(nèi)方向移動，從開挖面到圍巖內(nèi)部，位移逐漸減小。由于泥巖強度和剛度相對鹽巖均較低，因此在夾層位置變形較大。另外在腔體底面和頂面均出現(xiàn)較大位移，且愈靠近腔體中心位移愈大，這可能導(dǎo)致大幅度的底板隆起和地表沉降。

3 結(jié)論

本文建立了含不同夾層厚度的圓柱形儲庫單腔模型，考慮不同恒壓作用下長期蠕變，分析了運行內(nèi)壓及軟夾層厚度對鹽巖腔體穩(wěn)定性的影響，得到以下基本結(jié)論。

3.1 在恒定內(nèi)壓作用下，溶腔體積收縮率隨時間而增加，增長率呈下降趨勢;隨著內(nèi)壓增加，相同時間段內(nèi)腔體收縮率顯著減少。因此在儲氣庫運行期間，維持較高的運行壓力，盡可能降低低壓下運行時間，對腔體的可用性是有利的。

3.2 隨著溶腔內(nèi)壓的增加，腔周塑性區(qū)明顯減小;當腔體內(nèi)壓達到一定值時(本文為10 MPa)，溶腔圍巖已無明顯破損區(qū)，說明溶腔內(nèi)壓大小對腔體的穩(wěn)定性有較大影響。

3.3 同恒壓作用下，隨著夾層厚度在一定范圍內(nèi)增長(本文夾層厚度最大值與最小值比值為3)，腔體收縮率呈略微上漲，影響并不顯著，且隨內(nèi)壓的增大，影響值愈小。

3.4 隨夾層厚度增加，夾層位置處變形呈增大趨勢，但隨內(nèi)壓的增加，這種變形增大趨勢明顯減小，當內(nèi)壓達到一定值時(本文為15 MPa)，影響值已不明顯;而圍巖塑性區(qū)對夾層厚度的變化反應(yīng)并不明顯。

3.5 腔體較大位移主要集中在開挖面上，向腔內(nèi)方向移動，從開挖面到圍巖內(nèi)部，位移逐漸減小;因夾層力學特性與鹽巖不匹配，夾層位置出現(xiàn)較大變形。

［1］ 劉江.層狀鹽巖力學特性試驗研究及理論分析［D］.武漢：中國科學院武漢巖土力學研究所，2006.

［2］ Guangjie Y，Shen R，Tian Z，et al.Review of underground gas storage in the bedded salt deposit in China［C］//SPE Gas Technology Symposium，Calgary Alberta Canada：Society of Petroleum Engineers，2006.

［3］ 梁衛(wèi)國，楊春和，趙陽升.層狀鹽巖儲氣庫物理力學特性與極限運行壓力［J］.巖石力學與工程學報，2008，27(1)：22-27.

［4］ 王安明.層狀鹽巖變形機理及非線性蠕變本構(gòu)模型［D］.武漢：中國科學院武漢巖土力學研究所，2008.

［5］ 吳文，候正猛，楊春和.鹽巖中能源(石油和天然氣)地下儲存庫穩(wěn)定性評價標準研究［J］.巖石力學與工程學報，2005，24(14)：2497-2505.

［6］ 楊春和.層狀鹽巖力學理論與工程［M］.北京：科學出版社，2009.

［7］ 楊花.壓氣蓄能過程中地下鹽巖儲氣庫穩(wěn)定性研究［D］.武漢：中國科學院武漢巖土力學研究所，2009.

［8］ Wawersik W R.Alternatives to a power-low creep model for rock salt at temperature below 160℃［C］//The Mechanical Behavior of Salt，Proceedings of the Second Conference，Clausthal：Trans.Tech.Publications，1984.103-131.

［9］ Stein E，Wetjen D.Equations for the description studies of thermo-mechanical creep of rock salt and parametric studies by the aid of finite element computations［C］//The Mechanical Behavior of Salt，Proceedings of the Second Conference，Clausthal：Trans.Tech.Publications，1984.445-459.

［10］ YANG Chunhe，Daemen J K，YIN Jianhua.Experimental investigation of creep behavior of salt rock ［J］.International Journal of Rock Mechanics and Mining Sciences，2000，36(2)：336-341.

［11］ 高小平，楊春和，吳文.鹽巖蠕變特性溫度效應(yīng)的實驗研究［J］.巖石力學與工程學報，2005，24(12)：2054-2059.

［12］ 陳鋒，李銀平，楊春和.云應(yīng)鹽礦蠕變特性試驗研究［J］.巖石力學與工程學報，2006，25(S1)：3022-3027.