層狀鹽穴難溶夾層爆破效應(yīng)數(shù)值模擬研究

潘 越， 李二兵， 尤業(yè)超， 濮仕坤， 段建立，譚躍虎

(1.中國人民解放軍陸軍工程大學(xué) 國防工程學(xué)院，南京 210007；2.蘇交科集團(tuán)股份有限公司，南京 210017)

由于鹽巖具有良好的低滲透性、蠕變性、損傷自愈性、易溶性等特性，深部鹽礦水溶開采所形成的鹽穴是國際上公認(rèn)的石油、天然氣等能源地下儲備首選方式，已在美國、德國、法國、加拿大等國家廣泛應(yīng)用[1]，我國目前也在江蘇金壇、江蘇淮安、湖北應(yīng)城、湖北潛江、河南平頂山等地規(guī)劃了大量鹽穴油氣地下儲庫。我國鹽巖礦床大多屬于湖相沉積所形成的層狀結(jié)構(gòu)，硬石膏、泥巖和鈣芒硝等難溶夾層較多[2]。在深部鹽礦水溶建腔過程中，難溶夾層的存在不僅影響溶腔速度，導(dǎo)致腔體形狀不規(guī)則，甚至突然垮塌給水溶建腔造成嚴(yán)重的安全隱患。研究表明，鹽井聚能爆破技術(shù)不失為促使難溶夾層可控垮塌和消除建腔安全隱患的有效辦法[3]。但由于鹽穴地下儲庫需要滿足穩(wěn)定性和密閉性要求，鹽井爆破一方面要促使難溶夾層可控垮塌，另一方面不能因爆破效應(yīng)過大而影響整個儲庫密閉性，因此在開展現(xiàn)場先導(dǎo)性試驗(yàn)之前，進(jìn)行爆破效應(yīng)數(shù)值模擬分析顯得十分必要和重要，更具有直接的工程指導(dǎo)意義。

目前，國內(nèi)外針對巖體爆破數(shù)值模擬開展了廣泛研究，肖正學(xué)等[4]根據(jù)室內(nèi)試驗(yàn)和國內(nèi)外礦山爆破的案例研究了不同初始應(yīng)力對爆破作用的影響；謝源[5]研究了高地應(yīng)力場對爆破裂紋的影響；袁青[6]基于不確定理論研究了花崗巖地下儲氣庫洞室圍巖爆破特征參量，獲取了最佳特征參量預(yù)測模型；Chen等[7]研究了深埋隧道爆破開挖圍巖損傷區(qū)與地應(yīng)力關(guān)系；夏祥等[8]通過數(shù)值模擬研究了巖體爆破裂紋的產(chǎn)生和擴(kuò)展機(jī)制，得到了粉碎區(qū)和裂隙區(qū)范圍；孫金山等[9]在Flac3D中模擬了爆破應(yīng)力波對鄰近隧洞的動力擾動；郭文章等[10]基于DYNA-2D模擬節(jié)理巖體爆破過程，量化了節(jié)理巖體在爆炸荷載作用下的應(yīng)力場；高科[11]進(jìn)行了巖石SHPB試驗(yàn)?zāi)M研究，證明了半正弦波加載方式是用于巖石類脆性材料大直徑SHPB測試的理想形式；Zhu等[12]利用數(shù)值分析方法研究了動靜組合加載SHPB試驗(yàn)過程中巖體的破壞和損傷特性，得到復(fù)合靜力作用下巖石動強(qiáng)度增加機(jī)理；Zhu等[13]基于模擬分析手段探討了軸壓、圍壓等因素對巖石動態(tài)特性的影響；詹金武等[14]通過LS-DYNA模擬了破碎泥巖注漿結(jié)石體SHPB試驗(yàn)，驗(yàn)證了與試驗(yàn)結(jié)果吻合較好的HJC模型及其參數(shù)；李瑩[15]借助LS-DYNA分析不同初始地應(yīng)力場條件對爆破的影響，結(jié)果表明初始地應(yīng)力對爆破作用影響在近爆區(qū)域較小而在遠(yuǎn)爆區(qū)域較大；張鳳鵬等[16]借助有限元軟件研究節(jié)理?xiàng)l件下巖體爆破機(jī)理，發(fā)現(xiàn)節(jié)理面產(chǎn)生的反射拉伸應(yīng)力波加強(qiáng)了破巖效果。陳江海等[17]對洋山港水下爆破震動進(jìn)行測試分析；陳春歌等[18]對水下爆破進(jìn)行數(shù)值模擬，并提出了人工干預(yù)措施；王宏[19]模擬分析了水下鉆孔爆破得到了水擊波的傳播規(guī)律。

上述分析可知，爆破效應(yīng)數(shù)值模擬研究主要集中在隧道與地下工程、礦山邊坡等領(lǐng)域。關(guān)于動載作用下鹽穴穩(wěn)定性研究分析很少[20]，且由于層狀鹽穴地質(zhì)構(gòu)造特殊性，鹽穴內(nèi)部注滿鹵水，使得上述研究不能直接運(yùn)用。為此，本文首先進(jìn)行三軸壓縮狀態(tài)下泥巖的SHPB室內(nèi)試驗(yàn)，獲得難溶夾層塑性隨動硬化材料模型參數(shù)，利用LS-DYNA建立鹽穴實(shí)例模型模擬分析泥巖夾層和鹽穴腔體爆破振動等效應(yīng)，以期為層狀鹽穴難溶夾層爆破方案設(shè)計(jì)和穩(wěn)定性控制提供依據(jù)。

1 數(shù)值模擬模型建立

1.1 鹽穴實(shí)例模擬背景

數(shù)值模擬實(shí)例為湖北應(yīng)城某埋深1 000 m的層狀鹽穴，其存在厚度8 m的難溶夾層(主要成分是泥巖)，夾層上部鹽巖層厚度約為90 m，夾層下部鹽巖層厚度約為20 m。建腔過程中不斷提升導(dǎo)管，使得溶腔由下而上逐層溶蝕，并在完成夾層下部鹽腔溶蝕后即裝藥進(jìn)行爆破(見圖1(a)) 促使難溶夾層可控垮塌。保證形成的夾層通道暢通，加快建腔速度，并最終建成如圖1(b)所示的鹽穴儲庫。為簡化溶腔建模，近似看作以距泥巖夾層底部上方21 m的點(diǎn)為圓心作半徑29 m的球體并由泥巖夾層底面剪裁得到的下半部分。

圖1 層狀鹽穴爆破建腔示意圖Fig.1 Schematic diagram of blastingconstruction of bedded salt cavern

1.2 動態(tài)力學(xué)試驗(yàn)裝置

考慮到層狀鹽穴難溶夾層埋深約1 000 m，所處地應(yīng)力約25 MPa。方秦等[21]通過三向壓力下SHPB試驗(yàn)得出了鹽巖的DIF隨應(yīng)變率變化曲線，得到了25 MPa下鹽巖動態(tài)力學(xué)特性參數(shù)。因此，本文主要針對夾層泥巖進(jìn)行圍壓25 MPa下的SHPB試驗(yàn)。

本次試驗(yàn)采用Φ75 mm帶主動圍壓的SHPB試驗(yàn)裝置(見圖2)。裝置主要由動力系統(tǒng)(提供初速度的氮?dú)馄亢蛷椞?、彈性桿件(子彈、入射桿和透射桿)、傳感器(應(yīng)變片、動態(tài)測速儀等)，入射桿和透射桿組成的壓桿系統(tǒng)，長度分別為300 mm，3 500 mm和2 000 mm，材質(zhì)均為高強(qiáng)度不銹鋼，密度為7 850 kg/m3，彈性模量210 GPa。

1-入射桿端部液壓油缸；2-試件環(huán)向圍壓裝置；3-透射桿端部液壓油缸；4-單作用液壓油缸圖2 帶主動圍壓SHPB試驗(yàn)裝置Fig.2 SHPB experimental device with confininge pressure system

1.3 試樣制作及試驗(yàn)過程

本試驗(yàn)試件均采自湖北應(yīng)城鹽礦泥巖夾層，其主要成分為灰質(zhì)泥巖平均含量為40%～50%，次要成分為CaSO4、Na2SO4和NaCl。

為了降低摩擦效應(yīng)(試件與壓桿之間的界面摩擦力)、慣性效應(yīng)，使試驗(yàn)過程中試件能較好的滿足均勻性假定，試件的長徑比宜選為0.5。本試驗(yàn)選用Φ75 mm的壓桿系統(tǒng)。為了使試件截面完全受壓，在加工試件時，直徑略小于壓桿的直徑(直徑差不超過3 mm，否則導(dǎo)致數(shù)據(jù)失真)，因此，試驗(yàn)前將試件加工成標(biāo)準(zhǔn)尺寸Φ74 mm×37 mm，打磨過程實(shí)時進(jìn)行平整度檢查(見圖3)，完成后，刷一層油漆，以防空氣中的水分影響試件力學(xué)性能。選取5枚理想試樣如圖4所示，試樣各物理特性及巖性特征見表1。

試驗(yàn)步驟：架設(shè)SHPB壓桿系統(tǒng)，在壓桿兩側(cè)架設(shè)提供反力的拉桿，沖擊桿對中→架設(shè)圍壓系統(tǒng)→粘貼應(yīng)變片、連接線路→在接觸斷面涂抹凡士林進(jìn)行空打試驗(yàn)→加軸壓空打試驗(yàn)→試樣準(zhǔn)備→三向壓力下沖擊試驗(yàn)，沖擊氣壓分別設(shè)定為0.7 MPa，0.9 MPa，1.1 MPa→波形信號分析→數(shù)據(jù)分析處理。

(a)泥巖巖芯

(b)取芯及巖樣打磨圖3 泥巖試樣制作Fig.3 Mudstone specimen making

圖4 試樣成品Fig.4 Specimens of the pre-test

編號厚度L/mm直徑d/mm質(zhì)量/g密度/(g·cm-3)試樣巖性18-435.4271.92327.722.276灰質(zhì)泥巖(純泥巖)4-12-336.3073.50373.582.424灰質(zhì)泥巖(含鹽)5-1-134.2473.57357.52.455灰質(zhì)泥巖(純泥巖)137-435.4573.47345.532.298灰質(zhì)泥巖(純泥巖)30-6-535.5873.27351.572.343灰質(zhì)泥巖(含鹽)

1.4 試驗(yàn)結(jié)果分析

將試驗(yàn)獲得的電信號波形圖數(shù)據(jù)導(dǎo)入SHPB數(shù)據(jù)處理軟件，在平衡圖中比較波與透射波近似重疊，說明試驗(yàn)過程滿足SHPB試驗(yàn)應(yīng)力均勻性假設(shè)，結(jié)果具有準(zhǔn)確性。利用三波法處理數(shù)據(jù)并導(dǎo)出動態(tài)應(yīng)力應(yīng)變曲線，如圖5所示，試驗(yàn)處理結(jié)果見表2。從圖5和表2可知，試樣具有明顯的應(yīng)變率效應(yīng)，隨著應(yīng)變率的增加，峰值應(yīng)力也普遍增加。

圖5 泥巖夾層動態(tài)應(yīng)力應(yīng)變曲線Fig.5 Dynamic stress-strain curve of mudstone interlayer

編號圍壓/MPa沖擊氣壓/MPa子彈速度/(m·s-2)應(yīng)變率/s-1峰值應(yīng)力/MPa18-4251.124.02435.3193.734-12-3251.117.97345.32224.275-1-1250.924.7331.87183.69137-4250.721.24219.84168.2830-6-5250.718.05272.28146.06

1.5 模型參數(shù)驗(yàn)證

王健等[22]通過對不同圍壓下泥巖靜態(tài)強(qiáng)度的擬合，獲得25 MPa圍壓下泥巖夾層峰值應(yīng)力為116.12 MPa。結(jié)合表2各應(yīng)變率下峰值強(qiáng)度計(jì)算得到25 MPa圍壓下動力放大系數(shù)DIF(見表3)。

表3 25 MPa圍壓下泥巖夾層DIF

在LS-DYNA中采用非線性塑形材料模型(PLSTIC-KINEMATIC)可以反應(yīng)應(yīng)變率效應(yīng)，并且模型簡單，材料模型參數(shù)較易獲取。巖石類材料的塑形屈服點(diǎn)與破壞點(diǎn)近似，故不考慮硬化，模型簡化為Cowper-Symonds關(guān)系式

(1)

式中：σy為動態(tài)條件下的屈服應(yīng)力(取峰值應(yīng)力)；σ0為靜態(tài)下的初始屈服應(yīng)力；C和p為應(yīng)變率效應(yīng)提高系數(shù)。用Cowper-Symonds關(guān)系式對DIF進(jìn)行擬合得到C和p參數(shù)值分別為254.4和1.7。根據(jù)應(yīng)城泥巖夾層靜態(tài)三軸試驗(yàn)結(jié)果[23]擬合得到25 MPa圍壓下靜態(tài)力學(xué)參數(shù)，綜合得到25 MPa圍壓下泥巖材料參數(shù)如表4所示。

表4 25 MPa圍壓下泥巖材料參數(shù)

在LS-DYNA中建立與SHPB實(shí)體試驗(yàn)相同尺寸的壓桿模型(入射桿、透射桿長度分別為3 500 mm,2 000 mm；直徑均為75 mm)，試件采用量取的尺寸。所有模型均采用Solid164單元。采用映射網(wǎng)格劃分方法將入射桿沿長度方向劃分350段，透射桿沿長度方向劃分200段，沿桿半徑方向劃分32段。試件軸向分成48段，徑向分成48段。整個系統(tǒng)為對稱結(jié)構(gòu)，為減少計(jì)算時間，建立1/4的模型進(jìn)行計(jì)算(見圖6)。為在接觸面單元發(fā)生材料失效時剩余單元保持接觸狀態(tài)，壓桿與試件之間采用面面侵蝕接觸。另外，考慮到透射桿比較短，而波相對較長，在透射桿的端部加入非反射邊界條件。對稱面采用對稱邊界條件，入射桿自由端端部為應(yīng)力邊界(加載面)。

圖6 SHPB模型及網(wǎng)格劃分示意圖Fig.6 SHPB model and grid partition schematic diagram

將試驗(yàn)測得的入射波作為沖擊荷載施加在入射桿的端部可以避免入射波波形對計(jì)算結(jié)果產(chǎn)生影響。完成荷載施加后，設(shè)置接觸、分析時間和能量控制等，通過調(diào)節(jié)材料模型參數(shù)獲得模擬工況下的波形圖，三波法處理得到模擬條件下的應(yīng)力應(yīng)變曲線。將數(shù)值模擬結(jié)果與SHPB實(shí)體試驗(yàn)結(jié)果比對(見圖7)，大部分曲線與試驗(yàn)所獲曲線較為接近，說明塑形隨動強(qiáng)化模型可以較好地?cái)M合泥巖動力特性，材料參數(shù)合理并可為后續(xù)爆破模擬提供參考。

圖7 圍壓25 MPa下模擬應(yīng)力應(yīng)變曲線與試驗(yàn)應(yīng)力應(yīng)變曲線Fig.7 Simulated stress-strain curves and experimental stress-strain curves under confining pressure of 25 MPa

1.6 材料參數(shù)確定

巖體材料模型均采用PLASTIC-KINEMATIC模型。上文對泥巖動力特性進(jìn)行了研究和驗(yàn)證，用25 MPa圍壓下泥巖的塑性隨動硬化模型參數(shù)填充泥巖夾層和下臥層材料；通過方秦等的研究，將25 MPa圍壓下鹽巖的塑形隨動硬化模型參數(shù)填充鹽巖材料，如表5所示。

表5 巖體材料參數(shù)

炸藥采用MATHIGH_EXPLOSIVE_BURN和EOS_JWL來描述高性能炸藥的材料模型和狀態(tài)方程。根據(jù)文獻(xiàn)[24]，獲得某高性能TNT炸藥材料參數(shù)和狀態(tài)方程參數(shù)，如表6所示。

鹵水通常以低濃度注入，為了簡化計(jì)算，鹵水狀態(tài)方程采用純水的參數(shù)，而密度略大于水的密度，取值1.025 g/cm3。

表6 炸藥材料參數(shù)

1.7 計(jì)算模型建立

在LS-DYNA中選用半徑的3.5倍作為模擬的計(jì)算模型范圍(模型半徑90 m，高度100 m，見圖8)，由于模型較大，內(nèi)部結(jié)構(gòu)并非線性尺寸，為了保證網(wǎng)格長短邊尺寸之比滿足計(jì)算要求同時達(dá)到計(jì)算精度要求，劃分80萬個單元，其中最小單元尺寸為0.5 m(指Lagrangian算法網(wǎng)格，ALE網(wǎng)格最小單元尺寸為0.2 m)。模型邊界條件設(shè)為25 MPa應(yīng)力邊界，對稱面設(shè)為對稱邊界條件。采用多物質(zhì)流固耦合法(流體和固體不共節(jié)點(diǎn))進(jìn)行爆炸數(shù)值模擬，避免了Lagrange方法中因水和炸藥產(chǎn)生大變形而導(dǎo)致的計(jì)算終止。

圖8 數(shù)值模擬模型Fig.8 Numerical simulation model

2 數(shù)值模擬結(jié)果分析

2.1 夾層爆破試算分析

為了確定爆破的效果以及為后期裝藥量研究設(shè)置合理的計(jì)算梯度，先進(jìn)行爆破試算。采用柱狀裝藥，假設(shè)炸藥半徑為7.5 cm，高度6 m，重量為43.2 kg(為四分之一藥柱的重量)，考慮到藥柱下部為鹽腔，為了降低端部爆破效應(yīng)，同時充分利用炸藥，炸藥端部離夾層底部約1 m的位置。由于模型中液體壓縮產(chǎn)生的沖擊波會對爆破作用分析產(chǎn)生干擾，待計(jì)算到0.115 s后所有波形均耗散，整個模型應(yīng)力趨于平穩(wěn)。故在LS-DYNA中利用重啟動將炸藥起爆時間更改為0.15 s，研究0.15～0.2 s過程中的爆破波形圖。依據(jù)文獻(xiàn)[25]將振動峰值速度和等效應(yīng)力作為兩項(xiàng)重要的判據(jù)。

2.1.1 爆破作用下腔體圍巖分析

考慮到鹽腔圍巖損傷過大會影響儲庫的密閉性，在鹽腔圍巖內(nèi)壁上均勻選取平均弧長5.23 m的5個點(diǎn)A～E(見圖9)。圖10為各個點(diǎn)的有效應(yīng)力和合速度的時程曲線。從圖中可以看出，鹽腔圍巖各點(diǎn)有效應(yīng)力變化較小、受沖擊波影響小，故鹽腔圍巖的安全主要受振動速度控制。對比兩圖，可以明顯看出E點(diǎn)距離藥柱間距最大，但是最先受到爆破沖擊的影響產(chǎn)生應(yīng)力的變化和振動，而A點(diǎn)距離藥柱端部僅9 m，卻是最后受到爆破沖擊的影響。分析可知，沖擊波傳遞到鹽腔圍巖主要通過兩種路徑，一是通過鉆井中鹵水傳遞，而作用于鹽腔圍巖，另一種則是通過夾層巖體傳遞。應(yīng)力波在巖體沖傳播的速度明顯大于在水中傳播的速度。故出現(xiàn)了A點(diǎn)在爆破作用下的響應(yīng)滯后于E點(diǎn)。從振動的絕對值來看，在該當(dāng)量下，A點(diǎn)的振動峰值最大，為15.31 cm/s，其次是E點(diǎn)，為13.61 cm/s。其余各點(diǎn)受爆破作用影響較小，振動速度均不超過 5 cm/s。

2.1.2 爆破作用下夾層動力分析

為了對爆破作用下夾層動力響應(yīng)有全面認(rèn)知，在夾層下部和距離夾層底部高度4.5 m處分別布置觀察單元F～I(xiàn)和J～M，上下對應(yīng)的各點(diǎn)距藥柱間距分別

圖9 鹽腔圍巖觀察單元位置示意圖Fig.9 Schematic diagram of the location of the observation unit of the surrounding rock of a salt cavity

(a) 有效應(yīng)力時程曲線

(b) 合速度時程曲線圖10 鹽腔圍巖觀察單元時程曲線Fig.10 The time history curve of the observation unit of the surrounding rock of the salt cavity

約為4.6 m，9.4 m，14.1 m和18.9 m。圖11為夾層觀察單元位置示意圖。圖12為各觀察單元的有效應(yīng)力和合速度的時程曲線。從圖中看出隨著與藥柱間距的增大，應(yīng)力波峰值和振動合速度峰值迅速降低。

圖11 夾層中部觀察單元位置示意圖Fig.11 Schematic diagram of the position of the observation unit in the middle of the interlayer

圖12 夾層觀察單元時程曲線Fig.12 Time history curve of interlayer observation unit

2.2 炸藥當(dāng)量對爆破作用的影響

為了探究炸藥當(dāng)量對夾層爆破產(chǎn)生的影響，依次選用4.8 kg，19.2 kg，49.2 kg和79.2 kg 4個當(dāng)量(為四分之一藥柱的重量)的炸藥模擬對破壞效果的影響。根據(jù)3.1節(jié)所選的觀察單元，提取爆破作用下各點(diǎn)的峰值有效應(yīng)力和峰值振動合速度，并分析其與炸藥當(dāng)量的關(guān)系。圖13為各觀察單元峰值應(yīng)力與炸藥當(dāng)量的關(guān)系。

圖13 觀察單元峰值應(yīng)力與炸藥當(dāng)量的關(guān)系Fig.13 The relationship between the peak stress of the observation unit and the equivalent of the explosive

從圖中可以看出，隨著炸藥當(dāng)量的增加，各觀察單元的有效應(yīng)力峰值一般都增加，但由于各測點(diǎn)的位置不同，其變化的規(guī)律也不盡相同。如圖13中的A～E點(diǎn)，隨著炸藥當(dāng)量變化，其有效應(yīng)力變化很小，故而說明，爆破作用下炸藥對鹽腔的直接力學(xué)破壞作用小，可以忽略。另外，從圖中可以看出，在離炸藥較近的觀察單元(如F點(diǎn)、J點(diǎn))，其峰值有效應(yīng)力受到爆破的作用影響較大，受到炸藥當(dāng)量的影響比較明顯，而距離藥柱越遠(yuǎn)，其峰值應(yīng)力受到爆破炸藥作用越小。

圖14是觀察單元振動合速度與炸藥當(dāng)量的關(guān)系，從圖中可以看出，隨著炸藥當(dāng)量的增加，同一個點(diǎn)的振動和速度也增加，距離藥柱較近的位置速度增加的幅度大，較明顯，而距離藥柱遠(yuǎn)的位置，振動速度增加的幅度較小。圖14(a)可以看出，鹽腔圍巖B，C和D振動速度受炸藥當(dāng)量影響較小，而A點(diǎn)和E點(diǎn)受炸藥當(dāng)量影響較大。圖14(b)和圖14(c)中反應(yīng)距離效應(yīng)比較明顯，距離近的觀察點(diǎn)隨藥柱當(dāng)量的曲線在距離遠(yuǎn)的曲線之上，而在圖14(c)中，A點(diǎn)和E點(diǎn)沒有這個規(guī)律，其原因在于，巖石傳播的波速及阻抗與水中傳播的波速及阻抗不同。

圖14 觀察單元振動合速度與炸藥當(dāng)量的關(guān)系Fig.14 Observe the relationship between unit vibration speed and explosive equivalence

2.3 炸藥位置對爆破作用的影響

基于上部分計(jì)算結(jié)果，振動速度是爆破損傷范圍的主要控制指標(biāo)。選用半徑7.5 cm，高度為6 m的炸藥，藥柱底端到夾層底部的距離直接關(guān)系到鹽腔圍巖的振動速度。為此，根據(jù)藥柱下端面與夾層板底部的距離(分別為0 m，0.5 m，1.0 m和1.5 m)設(shè)置為4個工況。分別討論藥柱端部對鹽腔圍巖及夾層振動的影響。

圖15為A～I(xiàn)觀察單元圍巖振動情況，從圖中可以看出，藥柱下端部與夾層的間距對圍巖振動影響較大。隨著藥柱下端部與夾層的間距逐漸增大，A～D點(diǎn)的峰值振動合速度逐漸降低，主要是由于，藥柱距離下端部的距離越小，藥柱端部距離鹽腔越近，沖擊擾動就越大。另外，E點(diǎn)位于鹽腔邊緣部位，其峰值應(yīng)力主要受到巖石中應(yīng)力波傳播的影響，而鹵水傳遞的沖擊波對其影響較小，故而隨著間距的變化，其值變化相對較小。圖15(b)中，隨著距離增大，峰值合速度先增大后減小。主要是當(dāng)間距為0 m時，夾層底部巖體主要受到?jīng)_擊壓作用，所以其振動速度較大。隨著距離的增大，夾層底部在沖擊波的作用下，會發(fā)生反射，從而提高了振動速度。當(dāng)距離較大后，夾層底部受到?jīng)_擊波和反射波都較小，故而其振動速度又發(fā)生降低。從H點(diǎn)和I點(diǎn)可以看出，觀察單元距離藥柱的絕對距離較大，其受到的沖擊波和反射波均較小，藥柱下端部與夾層的間距對其影響就相對較小，反映出來的速度變化就較小。

圖15 藥柱下端部與夾層底部間距對爆破振動的影響Fig.15 The influence of the distance between the bottom of the grain and the bottom of the interlayer on blasting vibration

2.4 爆破效應(yīng)綜合分析

爆破主要受兩個因素控制，一方面，需要防止鹽腔圍巖損傷過大，影響鹽腔的密閉性；另一方面，需要實(shí)現(xiàn)預(yù)期的爆破效果，獲得理想的泥巖夾層破壞半徑。由于巖石類材料抗拉強(qiáng)度遠(yuǎn)低于抗壓強(qiáng)度，以抗拉強(qiáng)度為準(zhǔn)則，根據(jù)簡單的平面波估算鹽腔圍巖的振動速度承受范圍，依此設(shè)定裝藥量和裝藥位置，驗(yàn)算夾層破壞是否符合預(yù)期效果。

(2)

式(2)可改寫為

(3)

式中：V為巖石的振動速度；[V]為臨界振動速度；c為應(yīng)力波在巖石長傳播的速度；[σt]為極限抗拉強(qiáng)度；E為彈性模量。

超聲波在巖石中傳播波速與圍壓的關(guān)系?？紤]到應(yīng)力波在巖石中傳播波速與聲波相近，此處引用超聲波波速代替應(yīng)力波波速，近似取4 350 m/s。

2.4.1 鹽腔破壞控制

根據(jù)李二兵的研究，鹽巖的抗拉強(qiáng)度為1.04 MPa。由于巖石的動態(tài)抗拉強(qiáng)度的應(yīng)變率效應(yīng)小，近似取動態(tài)抗拉強(qiáng)度等于準(zhǔn)靜態(tài)抗拉強(qiáng)度。根據(jù)獲得的鹽巖抗拉強(qiáng)度(直接拉伸強(qiáng)度)和鹽巖的彈性模量，近似獲取臨界振動速度為16 cm/s。選取圖9中各單元為鹽腔破壞控制點(diǎn)。從圖15(a)可以看出，A點(diǎn)受到藥柱端部影響最大，且當(dāng)藥柱底部距離夾層底部距離大于1.0 m時，距離的影響變小。因此，藥柱端部距離夾層底部1.0 m處為合理的安放位置。分析圖14(a)可知，在藥柱下端部距離夾層下部1 m時，裝藥量應(yīng)控制不超過49.2 kg(為四分之一藥柱的重量)。

2.4.2 夾層破壞控制

依據(jù)夾層破壞效應(yīng)，為了充分反映夾層的破壞范圍(包含爆破產(chǎn)生的直接效應(yīng)和間接效應(yīng))，以使夾層產(chǎn)生裂縫為破壞標(biāo)準(zhǔn)，即夾層拉應(yīng)力達(dá)到泥巖的抗拉強(qiáng)度。李二兵的研究給出泥巖抗拉強(qiáng)度2.19 MPa。為了更好的對計(jì)算結(jié)果定量分析，在模型圖中選擇高度4.5 m處水平面為觀察面，如圖16所示。在LS-DYNA中，以受壓為負(fù)，受拉為正，故選擇單元的最大主應(yīng)力為觀察項(xiàng)目，最大主應(yīng)力達(dá)到極限抗拉強(qiáng)度，即認(rèn)為夾層受拉開裂。

圖16 破壞范圍觀察示意圖Fig.16 Observation map of damage range

圖17為破壞效果與裝藥量的關(guān)系。從圖中可以看出，破壞效果與裝藥量有直接的關(guān)系，隨著裝藥量的增加，破壞范圍增大，但隨著裝藥量增加，破壞半徑增大幅度減小。按照四分之一裝藥量考慮，19.2 kg時破壞半徑為8.66 m，49.2 kg時破壞半徑為11.3 m，可見裝量從19.2 kg增加1.56倍到49.2 kg，破壞半徑僅從8.66 m擴(kuò)大了0.3倍到11.3 m。難溶夾層厚度8 m，一般8 m左右的破壞半徑足以滿足難溶夾層可控垮塌要求。因此，難溶夾層爆破四分之一裝藥當(dāng)量可選在19.2 kg左右，整體裝藥量約76.8 kg。

圖17 破壞范圍與裝藥量的關(guān)系Fig.17 The relationship between the range of damage and the amount of charge

3 結(jié) 論

本文通過SHPB試驗(yàn)獲取難溶夾層材料參數(shù)。使用塑性隨動硬化材料模型擬合，并使用LS-DYNA驗(yàn)證材料模型的正確性和模型參數(shù)的準(zhǔn)確性。對鹽穴實(shí)例建立數(shù)值模型進(jìn)行研究分析，具體結(jié)論如下：

(1)在LS-DYNA中采用塑形隨動強(qiáng)化模型可以較好地?cái)M合泥巖動力特性，與25 MPa圍壓下SHPB試驗(yàn)結(jié)果具有較好的一致性。

(2)作用于鹽腔圍巖應(yīng)力波通過鹵水和夾層巖體兩種介質(zhì)傳遞，應(yīng)力波在巖體沖傳播的速度明顯大于在水中傳播的速度。

(3)通過對本鹽穴實(shí)例進(jìn)行數(shù)值模擬分析，為滿足鹽腔穩(wěn)定性和密閉性的要求，若選取藥柱半徑7.5 cm，高度6 m，爆破炸藥當(dāng)量宜選在76.8 kg左右，并布置在藥柱下端面與夾層板底部距離1 m的位置可以達(dá)到促使夾層垮塌的目的，同時保證鹽巖溶腔圍巖的氣密性。