摘 要：鹽巖蠕變對(duì)鹽穴地下儲(chǔ)氣庫(kù)的溶腔形態(tài)有著較大的影響，易造成儲(chǔ)氣庫(kù)失效、泄漏、破裂等安全問(wèn)題。為此，利用美國(guó)ITASCA咨詢集團(tuán)公司開(kāi)發(fā)的FLAC3D陜速拉格朗日差分軟件，建立了鹽穴儲(chǔ)氣庫(kù)模型，對(duì)儲(chǔ)氣庫(kù)溶腔蠕變規(guī)律進(jìn)行了數(shù)值模擬，分析了在不同蠕變時(shí)間內(nèi)、不同內(nèi)壓下儲(chǔ)氣庫(kù)溶腔在不同直徑處的變形量，得到在鹽巖蠕變下，鹽穴儲(chǔ)氣庫(kù)溶腔形態(tài)的變化規(guī)律，并給出了相應(yīng)操作建議：①鹽穴儲(chǔ)氣庫(kù)溶腔不同直徑處產(chǎn)生的變形量是不同的，溶腔最大直徑處變形量最大，隨著直徑的增大變形量增大，且呈非均勻性增大，因此在建設(shè)鹽穴儲(chǔ)氣庫(kù)時(shí)，應(yīng)對(duì)儲(chǔ)氣庫(kù)溶腔的高徑比進(jìn)行嚴(yán)格地計(jì)算；②溶腔內(nèi)壓對(duì)儲(chǔ)氣庫(kù)溶腔變形量有直接作用，較高壓力下的鹽巖蠕變對(duì)儲(chǔ)氣庫(kù)溶腔形態(tài)變化的影響比低壓時(shí)小，因此建議在鹽穴儲(chǔ)氣庫(kù)實(shí)際運(yùn)行過(guò)程中，盡量避免低壓運(yùn)行，以保證儲(chǔ)氣庫(kù)溶腔形態(tài)的完整性。該成果為鹽穴儲(chǔ)氣庫(kù)的安全運(yùn)行與維護(hù)提供了可行的依據(jù)。
關(guān)鍵詞：鹽穴地下儲(chǔ)氣庫(kù) 鹽巖蠕變 溶腔形態(tài) 變化規(guī)律 數(shù)值模擬 溶腔內(nèi)壓 完整性 操作壓力
Numeilical simuhaion of the creep defonnation of underground salt cavern gas storages
Abstract：With the rapid development of underground salt cavern gas storages in China，how to scientifically analyze the deformation rule of salt cavern gas storages and ensure their shape integrity have become important issues in gas storage construction，operation maintenance and management．In the operation of salt cavern gas storages，due to the role of different internal pressures and external formation pressures，the deformation in different diameters of cavity gypsum layers is different．In this paper,FLAC3D was used to numerically simulate the creep law of salt cavern gas storages and analyze the deformation in different diameters of cavity gypsum layers under different pressures．The results showed that the internal pressure has a direct effect on the deformation of salt cavern gas storages cavity，and the deformation changes with different diameters of cavity which is inhomogeneous with the increase of diameters．Thus，the ratio of height to diameter of cavity must be controlled strictly in the construction of underground salt cavern gas storages and the operation pressures of gas inj ection and mining must be controlled strictly in the running processes of gas storages．
Keywords：underground salt cavern gas storage，salt rock creep，cavity form，deformation rule，numerical simulation，internal pressure，integrity，operation pressure
鹽穴地下儲(chǔ)氣庫(kù)(以下簡(jiǎn)稱鹽穴儲(chǔ)氣庫(kù))通過(guò)常規(guī)鉆井方法鉆穿鹽巖層，注入淡水進(jìn)行沖蝕，抽出鹵水后利用形成的溶腔進(jìn)行氣體儲(chǔ)存。由于鹽巖具有孔隙度低、滲透性小、塑性大等優(yōu)點(diǎn)，運(yùn)用鹽穴儲(chǔ)氣庫(kù)儲(chǔ)存氣體安全可靠、不易滲漏[1]。
鹽巖具有較大塑性，造成了儲(chǔ)氣庫(kù)溶腔注采氣時(shí)的可收縮性。同時(shí)，鹽巖具有蠕變特性，鹽巖蠕變對(duì)鹽穴地下儲(chǔ)氣庫(kù)的溶腔形態(tài)有較大影響，造成儲(chǔ)氣庫(kù)失效、泄漏、破裂等安全問(wèn)題[2]。鑒于此，利用美國(guó)ITASCA咨詢集團(tuán)公司開(kāi)發(fā)的FLAC3D快速拉格朗日差分軟件，建立鹽穴儲(chǔ)氣庫(kù)模型，對(duì)儲(chǔ)氣庫(kù)溶腔蠕變規(guī)律進(jìn)行數(shù)值模擬，分析在不同蠕變時(shí)間內(nèi)、不同內(nèi)壓下儲(chǔ)氣庫(kù)溶腔在不同直徑處的變形量，得到在鹽巖蠕變下，鹽穴儲(chǔ)氣庫(kù)溶腔形態(tài)的變化規(guī)律，為鹽穴儲(chǔ)氣庫(kù)的安全運(yùn)行與維護(hù)提供可行的依據(jù)。
1 建模原理
利用數(shù)值模擬方法對(duì)鹽穴儲(chǔ)氣庫(kù)溶腔蠕變規(guī)律及溶腔變形量開(kāi)展研究，基本原理是采用拉格朗日法研究每個(gè)流體質(zhì)點(diǎn)隨時(shí)間變化的情況，著眼于某個(gè)流體質(zhì)點(diǎn)，研究任一段時(shí)間內(nèi)質(zhì)點(diǎn)運(yùn)動(dòng)的軌跡以及質(zhì)點(diǎn)的速度和壓力。
FLAC3D將求解物體離散為一系列如圖1所示的四面體單元，并采用插值函數(shù)
進(jìn)行求解，函數(shù)中ui、vi分別代表四面體中節(jié)點(diǎn)的位移、速度。FLAC3D數(shù)值模擬基本原理是[3]：①通過(guò)對(duì)三維介質(zhì)的離散，使所有外力與內(nèi)力集中于三維網(wǎng)絡(luò)節(jié)點(diǎn)上，進(jìn)而將連續(xù)介質(zhì)運(yùn)動(dòng)定律轉(zhuǎn)化為離散節(jié)點(diǎn)上的牛頓定律；②時(shí)間與空間的層數(shù)采用沿有限空間與時(shí)間間隔線性變化的有限差分來(lái)近似地表示；③將靜力問(wèn)題當(dāng)做動(dòng)力問(wèn)題來(lái)求解，運(yùn)動(dòng)方程中慣性項(xiàng)用來(lái)作為達(dá)到所求靜力平衡的一種手段。
2 溶腔模型的建立
我國(guó)鹽巖礦藏以多夾層鹽巖為主，由于主要分析鹽穴儲(chǔ)氣庫(kù)溶腔在鹽巖蠕變下的變形規(guī)律，因此在建立儲(chǔ)氣庫(kù)溶腔模型時(shí)將復(fù)雜的鹽巖夾層簡(jiǎn)化為簡(jiǎn)單的單層鹽巖層。借鑒江蘇金壇鹽穴儲(chǔ)氣庫(kù)的建庫(kù)參數(shù)[4]，建立了如圖2所示的腔體幾何模型。建立溶腔及圍巖地層模型后，將地層巖石力學(xué)參數(shù)、鹽層的蠕變關(guān)系、溶腔周圍的地應(yīng)力參數(shù)及溶腔運(yùn)行參數(shù)帶入模型中進(jìn)行模擬計(jì)算。
2．1 網(wǎng)格劃分及參數(shù)設(shè)置
利用ANSYS軟件進(jìn)行模型網(wǎng)格劃分[5]，并接入到FALC3D中，劃分后的網(wǎng)格如圖3所示。
圖3中，1/4圓柱體模擬儲(chǔ)氣庫(kù)溶腔周圍的地質(zhì)環(huán)境，其直徑是儲(chǔ)氣庫(kù)溶腔模型最大直徑的20倍，在網(wǎng)格中間位置可以看到儲(chǔ)氣庫(kù)溶腔輪廓的1/4。在網(wǎng)格劃分時(shí)將儲(chǔ)氣庫(kù)溶腔網(wǎng)格一起劃分，有助于儲(chǔ)氣庫(kù)溶腔的形成。在此網(wǎng)格中，定義了2個(gè)組，分別是group1和group2，前者對(duì)儲(chǔ)氣庫(kù)溶腔進(jìn)行網(wǎng)格劃分，后者對(duì)儲(chǔ)氣庫(kù)溶腔周圍地質(zhì)環(huán)境進(jìn)行網(wǎng)格劃分。選擇采用各向同性彈性模型，其材料特性常數(shù)見(jiàn)表1。
2．2 模型求解及驗(yàn)證
數(shù)值分析采用迭代方法計(jì)算，在迭代過(guò)程中監(jiān)控一些變量或參數(shù)的變化，以判斷分析結(jié)果是否正確，模型是否與實(shí)際相符[6-7]。監(jiān)控模型最大不平衡力，所謂體系最大不平衡力，是指每一個(gè)計(jì)算循環(huán)中，外力通過(guò)網(wǎng)格節(jié)點(diǎn)傳遞分配到體系各節(jié)點(diǎn)時(shí)，所有節(jié)點(diǎn)的外力與內(nèi)力之差中的最大值。如圖4～6監(jiān)控了最大不平衡力，并記錄了位移等值線圖和應(yīng)力等值線圖。
由圖4可以得出，該模型中最大不平衡力逐漸收斂至零，表示所建儲(chǔ)氣庫(kù)溶腔模型具有收斂性。由圖4、5可以看出，在建立的土體模型中，由于重力場(chǎng)的作用，最大位移和最大應(yīng)力分別出現(xiàn)在模型的最底部。
2．3 腔體的形成及導(dǎo)入蠕變模型
利用空模型形成腔體模型，通過(guò)設(shè)置groupl為空模型進(jìn)行開(kāi)挖，開(kāi)挖后的模型如圖7、8所示。
形成腔體模型后，需要對(duì)溶腔施加內(nèi)壓，通過(guò)施加法向應(yīng)力將儲(chǔ)氣庫(kù)內(nèi)壓分別設(shè)置為10MPa、15MPa、17MPa、20MPa、23MPa、25MPa和30MPa，選擇經(jīng)典黏彈性蠕變模型[8-9]計(jì)算在不同內(nèi)壓條件下溶腔不同直徑處的蠕變變形量。經(jīng)典黏彈性蠕變模型材料參數(shù)如表2所示。
3 蠕變計(jì)算及結(jié)果分析
3．1 蠕變計(jì)算
根據(jù)一些學(xué)者的研究[10-12]，進(jìn)行蠕變計(jì)算時(shí)，分別計(jì)算了儲(chǔ)氣庫(kù)溶腔內(nèi)壓為10MPa、15MPa、17MPa、20MPa、23MPa、25MPa、30MPa，蠕變時(shí)間分別為2a、4a、6a、8a、10a、15a、20a時(shí)，儲(chǔ)氣庫(kù)溶腔直徑為20m、30m、40m、50m處X軸向的位移變化量。
3．2 結(jié)果分析
根據(jù)模擬計(jì)算所得結(jié)果，可以得到不同壓力下儲(chǔ)氣庫(kù)溶腔不同直徑處變形量圖(圖9～11)，經(jīng)總結(jié)分析，僅列出具有代表性的分析圖。圖9～11分別表示溶腔內(nèi)壓為10MPa、20MPa、30MPa時(shí)，不同蠕變時(shí)間鹽穴儲(chǔ)氣庫(kù)溶腔不同直徑處的變形量。
對(duì)比分析圖9～11，可以得出：①溶腔不同直徑處變形量隨時(shí)間變化規(guī)律比較一致，溶腔直徑為50m處的變形量最大；②蠕變時(shí)間為6a時(shí)，不同壓力下溶腔不同直徑處變形量曲線均出現(xiàn)一個(gè)明顯的轉(zhuǎn)折點(diǎn)。蠕變時(shí)間小于6a時(shí)，溶腔變形量曲線斜率較大，表明在此時(shí)間內(nèi)溶腔變形速率較大，變形比較顯著；蠕變時(shí)間大于6a時(shí)，溶腔變形量曲線較平緩，溶腔變形量隨時(shí)間變化不大；③內(nèi)壓為10MPa時(shí)，蠕變20a后，溶腔直徑為50m處變形量為360mm。內(nèi)壓為20MPa時(shí)，蠕變20a后，溶腔直徑為50m處變形量為620mm。內(nèi)壓為30MPa時(shí)，蠕變20a后，溶腔直徑為50m處變形量為640mm。數(shù)據(jù)表明，溶腔內(nèi)壓小于20MPa時(shí)，溶腔變形量隨壓力增大顯著增加；溶腔內(nèi)壓大于20MPa時(shí)，溶腔變形量變化不大。
4 結(jié)論
1)建立了鹽穴儲(chǔ)氣庫(kù)溶腔模型，計(jì)算分析了在鹽巖蠕變下鹽穴儲(chǔ)氣庫(kù)溶腔形態(tài)的變形量，得到了鹽穴儲(chǔ)氣庫(kù)在不同壓力下溶腔不同直徑處的變化規(guī)律。由于受鹽巖蠕變的影響，在鹽穴儲(chǔ)氣庫(kù)溶腔不同直徑上產(chǎn)生的變形量是不同的，溶腔最大直徑處變形量最大，隨著直徑的增大變形量增大，且呈非均勻增大，因此在建設(shè)鹽穴儲(chǔ)氣庫(kù)時(shí)，應(yīng)對(duì)儲(chǔ)氣庫(kù)溶腔的高徑比進(jìn)行嚴(yán)格的計(jì)算。
2)鹽穴儲(chǔ)氣庫(kù)在不同的運(yùn)行壓力下變形量是不同的。當(dāng)壓力小于20MPa時(shí)，隨著內(nèi)壓增加，溶腔整體變形量增加；當(dāng)內(nèi)壓大于20MPa時(shí)，隨著內(nèi)壓增加，溶腔整體變形量基本保持不變。說(shuō)明在較高壓力下，鹽巖蠕變對(duì)溶腔形態(tài)變化的影響比低壓時(shí)小。因此，建議在鹽穴儲(chǔ)氣庫(kù)實(shí)際運(yùn)行過(guò)程中，盡量避免低壓運(yùn)行，以保證儲(chǔ)氣庫(kù)溶腔形態(tài)完整性。
參考文獻(xiàn)
[1]丁國(guó)生．中國(guó)地下儲(chǔ)氣庫(kù)現(xiàn)狀與發(fā)展展望[J]．天然氣工業(yè)，2006，26(6)：111-113．
DING Guosheng．Current situation and prospect of Chinese underground natural gas storage[J]．Natural Gas Industry，2006，26(6)：111-113．
[2]王同濤．多夾層鹽巖體中儲(chǔ)氣庫(kù)圍巖變形規(guī)律及安全性研究[D]．北京：中國(guó)石油大學(xué)，2011．
WANG Tongtao．Study on deformation and safety of rack mass surrounding gas storage in bedded rock salt[D]．Beijing：China University of Petroleum，2011．
[3]丁國(guó)生．鹽穴儲(chǔ)氣庫(kù)溶腔收縮規(guī)律分析[J]．天然氣工業(yè)，2007，27(11)：94-96．
DING Guosheng．The creep deformation rule of gas storage salt cavern[J]．Natural Gas Industry，2007，27(11)：94-96．
[4]尹雪英，楊春和，陳劍文．金壇鹽礦老腔儲(chǔ)氣庫(kù)長(zhǎng)期穩(wěn)定性分析數(shù)值模擬[J]．巖土力學(xué)，2006，27(6)：869-874．
YIN Xueying,YANG Chunhe,CHEN Jianwen．Numerical simulation research on long term stability of gas storage in Jintan Salt Mine[J]．Rock and Soil Mechanics，2006，27(6)：869-874．
[5]譚羽非．天然氣地下儲(chǔ)氣庫(kù)技術(shù)及數(shù)值模擬[M]．北京：石油工業(yè)出版社，2007．
TAN Yufei．Underground gas storage technology and numerical simulation[M]．Beijing：Petroleum Industry Press，2007．
[6]王新勝．鹽巖儲(chǔ)氣庫(kù)運(yùn)營(yíng)期穩(wěn)定性研究[D]．重慶：重慶大學(xué)，2009．
WANG Xinsheng．Study on operation period stability of gas storage of salt rock[D]．Chongqin9：Chongqing University，2009．
[7]李后榮．鹽巖地下儲(chǔ)氣庫(kù)穩(wěn)定性分析[J]．貴州大學(xué)學(xué)報(bào)，2011，28(1)：122-125．
LI Hourong．The stability analysis of underground salt rock gas storage[J]．Journal of Guizhou University，2011，28(1)：122-125．
[8]梁衛(wèi)國(guó)，徐素國(guó)．金壇儲(chǔ)氣庫(kù)鹽巖蠕變特性及其實(shí)用本構(gòu)研究[J]．遼寧工程技術(shù)大學(xué)學(xué)報(bào)，2007，26(3)：354-356．
LIANG Weiguo，XU Suguo．Research on creep property and practical constitutive equation of rock salt in Jintan gas storage of ChinaEJ]．Journal of Liaoning Technical University，2007，26(3)：354-356．
[9]楊春和，陳鋒，曾義金．鹽巖蠕變損傷關(guān)系研究[J]．巖石力學(xué)與工程學(xué)報(bào)，2002，14(11)：1602-1604．
YANG Chunhe，CHEN Feng，ZENG Yijin．Investigation on creep damage constitutive theory of salt rock[J]．Chinese Journal of Rock Mechanics and Engineering，2002，14(11)：1602-1604．
[10]KEELEY DEBORAH．Failure rates for underground gas storage[M]．London：Health and Safety Laboratory，2008．
[11]MOGHADAM S N，MIRZABOZORG H．Modeling time dependent behavior of gas caverns in rock salt considering creep，dilatancy and failure[J]．Tunnelling and Underground Space Technology，2013(33)：171-185．
[12]姜德義．巖鹽溶腔穩(wěn)定性及失穩(wěn)控制研究[D]．重慶：重慶大學(xué)，2003．
JIANG Deyi．Research on stability and destabilization control of rock salt cavern[D]．Chongqing：Chongqing University，2003．
本文作者：梁光川 王夢(mèng)秋 彭星煜 張碩
作者單位：西南石油大學(xué)石油與天然氣工程學(xué)院
中國(guó)石油兩氣東輸管道公司蘇北管理處
您可以選擇一種方式贊助本站
支付寶轉(zhuǎn)賬贊助
微信轉(zhuǎn)賬贊助
