張前進(jìn), 武 科, 崔帥帥, 于雅琳, 張樂(lè)文

(1.山東大學(xué) 土建與水利學(xué)院，山東 濟(jì)南 250061; 2.山東大學(xué) 海洋研究院，山東 青島 266237)

隨著社會(huì)的不斷發(fā)展，地下工程已成為城市建設(shè)的重要部分。隧道能夠壓縮城市空間，方便人們的出行，在交通運(yùn)輸中起著重要的作用，被大量運(yùn)用于城市建設(shè)中。但是，地下空間有限，有時(shí)會(huì)遇到隧道下穿情況，下穿隧道的建設(shè)勢(shì)必會(huì)給既有隧道帶來(lái)影響，產(chǎn)生地面沉降、既有隧道變形及土層變形等危害。由于用礦山法開挖隧道的爆破會(huì)產(chǎn)生爆轟波，嚴(yán)重影響了既有結(jié)構(gòu)的穩(wěn)定性[1-7]，因此，在地鐵爆破過(guò)程中，需要研究工程爆破對(duì)既有結(jié)構(gòu)的影響，預(yù)測(cè)結(jié)構(gòu)受力變形。

關(guān)于工程爆破振動(dòng)的研究有很多。周俊汝[8]等人基于粘性巖體爆破振動(dòng)頻譜表達(dá)式,并結(jié)合LS-DYNA有限元軟件,對(duì)球狀藥包爆破振動(dòng)主頻和平均頻率的衰減機(jī)制及其規(guī)律進(jìn)行了分析。張震[9]等人采用現(xiàn)場(chǎng)監(jiān)測(cè)和動(dòng)力有限元數(shù)值模擬相結(jié)合的研究方法,對(duì)超淺埋地鐵站通道爆破開挖鄰近埋地混凝土管道的動(dòng)力響應(yīng)進(jìn)行了研究。鄧鍔[10]等人運(yùn)用數(shù)值模擬方法,分別研究了不同頂板厚度下隧道爆破施工引起的圍巖振速分布特征及其對(duì)圍巖的損傷情況,并對(duì)圍巖穩(wěn)定性進(jìn)行了分析。Jun[11]等人為了研究隧道爆破引起的沖擊波傳播衰減，建立了具有爆破效果的三維模型。

有關(guān)學(xué)者研究振動(dòng)載荷的影響集中于地表的振動(dòng)速度和沖擊波的衰減規(guī)律,而研究下穿既有隧道時(shí)爆破施工對(duì)既有結(jié)構(gòu)的影響及選擇合理的施工方法較少[12-14]。作者擬在數(shù)值模擬軟件Flac3D的基礎(chǔ)上，結(jié)合具體工程實(shí)例，對(duì)爆破載荷作用下新建隧道與既有隧道的應(yīng)力和變形情況進(jìn)行分析，研究地層和拱頂沉降，以期對(duì)相關(guān)施工進(jìn)行指導(dǎo)。

1 工程概況

既有地鐵隧道采用土壓平衡盾構(gòu)建造。該盾構(gòu)隧道外徑為6.0 m，內(nèi)徑為5.4 m。隧道管片由六塊拼裝而成，管片寬度為1.5 m。新建地鐵隧道部分區(qū)間采用礦山法建造，新、舊隧道夾角約為85°。新建地鐵隧道下穿既有隧道，新建隧道埋深18.2 m，新建隧道與既有隧道右線與左線最小凈距分別為7.5和1.8 m。為了確保下穿期間既有隧道的安全，穿越過(guò)程中對(duì)既有隧道變形進(jìn)行了實(shí)時(shí)監(jiān)測(cè)。沿既有隧道與新建隧道的交界處布置了斷面檢測(cè)，斷面的間隔5～6 m。此外，對(duì)隧道附近同一埋深帶點(diǎn)的地表沉降進(jìn)行了檢測(cè)，斷面距隧道邊緣2～3 m。隧道的地層和監(jiān)測(cè)點(diǎn)布置分別如圖1，2所示。

圖1 地層Fig. 1 Stratigraphy

圖2 監(jiān)測(cè)點(diǎn)布置(單位：m)Fig. 2 The layout of the monitoring point(unit:m)

2 數(shù)值計(jì)算模型與分析方法

2.1 數(shù)值近似模型

采用 FLAC 3D有限差分軟件，對(duì)下穿既有隧道的爆破過(guò)程進(jìn)行數(shù)值模擬，研究新建隧道爆破對(duì)既有隧道的動(dòng)力響應(yīng)。該模型為馬蹄形隧道，由既有隧道和土體組成。模型中新建隧道橫向取54 m，豎向取50 m，軸向取70 m。模型土體采用Mohr-Coulomb模型，新建隧道襯砌與既有隧道按照彈性模型進(jìn)行計(jì)算。隧道模型如圖3所示。

圖3 隧道模型Fig. 3 Tunnel model

2.2 巖土介質(zhì)力學(xué)參數(shù)

根據(jù)現(xiàn)場(chǎng)地質(zhì)勘查情況，隧址區(qū)土體由素填土、細(xì)粉砂、礫質(zhì)粘性土、全風(fēng)化花崗巖及強(qiáng)風(fēng)化花崗巖等組成。具體土層和支護(hù)參數(shù)見表1。

2.3 數(shù)值計(jì)算分析方法

2.3.1 數(shù)值計(jì)算步驟

工程爆破數(shù)值模擬計(jì)算分為靜力計(jì)算和動(dòng)力計(jì)算2部分。在靜力模擬的基礎(chǔ)上才能進(jìn)行動(dòng)力模擬。Flac3D數(shù)值模擬步驟為：① 建立模型，進(jìn)行網(wǎng)格劃分并確定動(dòng)力參數(shù)；② 定義靜力條件下的初始邊界條件；③ 進(jìn)行重力平衡計(jì)算，然后進(jìn)行隧道靜力開挖模擬，得到開挖后的結(jié)果；④ 選擇動(dòng)力模式，設(shè)置動(dòng)力邊界條件和阻尼條件；⑤ 施加動(dòng)力載荷，進(jìn)行動(dòng)力計(jì)算，得到模擬結(jié)果。

2.3.2 動(dòng)載荷施加

在動(dòng)力計(jì)算中，邊界條件包括粘性邊界條件和自由邊界條件。動(dòng)力計(jì)算模擬開挖使用粘性邊界。采用粘滯邊界條件就是在模擬行的前、后、左、右邊界以及下邊界設(shè)置粘滯邊界條件，上邊界為自由邊界。力學(xué)阻尼包括瑞利阻尼、局部阻尼和滯后阻尼?，F(xiàn)采用瑞利阻尼，其臨界阻尼比為5%，中心頻率取100 Hz。

為了準(zhǔn)確地模擬隧道施工過(guò)程中的沖擊載荷，采用載荷作用時(shí)間和爆破沖擊載荷經(jīng)驗(yàn)公式來(lái)施加動(dòng)力載荷[12]。

(1)

(2)

(3)

式中：PD(t)為爆破沖擊載荷；B為載荷常量；t為作用時(shí)間；PB為實(shí)測(cè)爆轟速度D和密度ρ0時(shí)的炸藥爆轟壓力；Sge為炸藥密度；ve為爆轟速度；μ為泊松比；m為裝藥質(zhì)量[14]。

施加動(dòng)載荷方式為：將動(dòng)載荷施加在垂直隧道壁新建隧道的周圍。在模擬靜載荷下施工完成后，施加動(dòng)載荷，位置為一次爆破循環(huán)作用的管壁，動(dòng)載荷總體作用時(shí)間為0.2 s，動(dòng)載荷施加的方向?yàn)榇怪庇诠鼙谙蛲猓鐖D4所示。通過(guò)實(shí)際測(cè)量，得到的參數(shù)有：實(shí)測(cè)爆轟速度為3.5 km/s；密度為1 g/cm3；泊松比為0.23；裝藥質(zhì)量為28 kg。

圖4 爆破時(shí)程載荷Fig. 4 Blasting time versus the load

3 結(jié)果分析

3.1 新建地鐵隧道穩(wěn)定性分析

3.1.1 應(yīng)變分析

爆破載荷的作用下，新建隧道和土體的應(yīng)變會(huì)發(fā)生改變。隧道的位移云圖如圖5所示。

從圖5中可以看出，靜載下豎向位移最大為2.57 cm，施加爆破載荷后最大豎向位移為2.83 cm，豎向位移增加了0.26 cm。豎向位移沉降最大處在拱頂部分，最大隆起出現(xiàn)在拱底。靜載下水平位移為0.77 cm，施加爆破載荷后水平位移的最大值為1 cm，增加了2.3 mm。

經(jīng)過(guò)施加動(dòng)載荷，豎向沉降增加了約10%，水平位移增加了約30%。這是由于豎向位移受重力作用影響較大，動(dòng)載荷影響沉降不是主要作用因素，因此豎向位移增加得不多。水平位移受重力作用的影響較小，但施加動(dòng)載荷后其所受影響較大。在爆破開挖的隧道中，不僅要考慮豎向位移，更要加強(qiáng)水平方向的支護(hù)措施來(lái)減小爆破振動(dòng)所產(chǎn)生的影響。

圖5 位移云圖Fig. 5 Displacement cloud

在模擬過(guò)程中，對(duì)斷面MC-C爆破區(qū)間內(nèi)5個(gè)點(diǎn)的豎向位移進(jìn)行了檢測(cè)，檢測(cè)點(diǎn)位置和豎向位移的變化分別如圖6，7所示。

圖6 檢測(cè)點(diǎn)位置Fig. 6 Test point location

圖7 豎向位移與時(shí)間的關(guān)系Fig. 7 Relationship between vertical displacement and time

從圖7中可以看出，隧道的拱頂與拱腰處的監(jiān)測(cè)點(diǎn)位移趨勢(shì)相同。經(jīng)過(guò)施加爆破載荷后，位移先減小后增加，最后趨于穩(wěn)定。檢測(cè)點(diǎn)1的位移由-2.5 cm減小到-1.3 cm。隨著爆破載荷逐漸減小，位移又逐漸增加，最后，穩(wěn)定在-2.8 cm附近。拱底的位移在動(dòng)載荷的作用下由隆起變成沉降，然后，隨著載荷的不斷減小，又變?yōu)槁∑?。檢測(cè)點(diǎn)5的位移由-0.8 cm變?yōu)? cm，然后，又逐漸增加，最后，穩(wěn)定在-1.1 cm附近。

監(jiān)測(cè)點(diǎn)的變化趨勢(shì)與動(dòng)載荷的變化趨勢(shì)有關(guān)。在0.02 s處位移出現(xiàn)明顯轉(zhuǎn)折，0.06 s后動(dòng)載荷消失，位移趨于穩(wěn)定。經(jīng)動(dòng)載荷作用后，隧道上部分為沉降且位移增加，拱底位移略有減小。在爆破施工過(guò)程中，由于爆破作用，隧道的位移會(huì)超過(guò)最終的穩(wěn)定值，并且中間波動(dòng)量較大，因此，對(duì)隧道的防護(hù)要考慮到爆破過(guò)程中位移的變化，對(duì)于隧道穩(wěn)定性的防護(hù)要留有一定的安全空間，防止因爆破波動(dòng)過(guò)大而造成的破壞。

3.1.2 應(yīng)力分析

為了分析襯砌的受力情況，對(duì)隧道襯砌的最大正應(yīng)力與最小正應(yīng)力進(jìn)行了分析，應(yīng)力云圖如圖8所示。從圖8中可以看出，在施加動(dòng)載后，最小正應(yīng)力由3.2 MPa變?yōu)?.7 MPa，增加了16%。在爆破的作用下，最大正應(yīng)力由1.7 MPa增加到2 MPa。動(dòng)載荷對(duì)最小正應(yīng)力和最大正應(yīng)力應(yīng)力的影響都比較大，施加動(dòng)載荷后應(yīng)力整體增加了約16%，因此，要加強(qiáng)材料和結(jié)構(gòu)的強(qiáng)度來(lái)降低爆破振動(dòng)的影響。

圖8 應(yīng)力云圖Fig. 8 Stress cloud

新建隧道最小主應(yīng)力受壓，并分布在拱腰附近，拱頂受力較小。在動(dòng)載荷的作用下，最小主應(yīng)力在拱腰部分的壓應(yīng)力逐漸增加，而既有隧道的壓應(yīng)力有所減小。新建隧道的最大正應(yīng)力產(chǎn)生在拱頂位置，拱腰受壓，既有隧道在拱頂和拱底受壓較大。在動(dòng)載荷的作用下，新建隧道拱腰的壓應(yīng)力作用區(qū)域增加，而壓應(yīng)力的最大值減小，拉應(yīng)力增加了18%。因此，動(dòng)載荷使新建隧道和既有隧道的應(yīng)力都增加，但是爆破作用對(duì)新建隧道的影響要大于既有隧道的。處于隧道爆破位置的主應(yīng)力的變化最為明顯，其原因是振動(dòng)作用距離爆破位置越遠(yuǎn)，振動(dòng)作用越小。要重點(diǎn)加強(qiáng)爆破位置的防護(hù)，防止局部結(jié)構(gòu)因超過(guò)極限條件而遭受破壞。

3.2 既有地鐵隧道穩(wěn)定性分析

3.2.1 既有結(jié)構(gòu)地層沉降

對(duì)MC-A和MC-B點(diǎn)的地層沉降進(jìn)行了檢測(cè)分析，地層沉降曲線如圖9所示。

從圖9中可以看出，MC-A斷面最大沉降的實(shí)測(cè)值為9.5 mm，而最大沉降的模擬值為8.1 mm。MC-B斷面最大沉降的實(shí)測(cè)值為7.8 mm，而最大沉降的模擬值為6.2 mm。模擬值比實(shí)測(cè)值小1 mm左右。實(shí)際測(cè)量值與模擬值曲線趨勢(shì)相同，距離隧道中心處的位移最大。隨著與新建隧道距離的增加，地層沉降不斷減小。因此，要在新建隧道附近區(qū)域進(jìn)行重點(diǎn)防護(hù)。另外，既有隧道二附近斷面MC-A的沉降大于既有隧道一附近的斷面MC-B的沉降。表明：距爆破點(diǎn)較遠(yuǎn)的既有隧道一受到爆破的影響比距爆破點(diǎn)較近的既有隧道二受到爆破的影響要小。用Flac 3D建立的數(shù)值計(jì)算模型偏差較小，其結(jié)果較為準(zhǔn)確。

3.2.2 既有隧道拱頂沉降

根據(jù)圖2監(jiān)測(cè)點(diǎn)的布置，對(duì)隧道拱頂和地層的沉降進(jìn)行監(jiān)測(cè)。各點(diǎn)的沉降如圖10所示。

從圖10中可以看出，既有隧道一在靜載下拱頂?shù)淖畲蟪两禐?.5 mm。施加動(dòng)載后，拱頂?shù)淖畲蟪两禐?.9 mm，增加了0.4 mm，約為11.4%。既有隧道二在靜載下拱頂?shù)淖畲蟪两禐?.47 mm。發(fā)生爆破后，拱頂?shù)淖畲蟪两禐? mm，增加了1.53 mm，約為44%。動(dòng)載荷對(duì)既有隧道二的影響大于對(duì)既有隧道一的，其原因是既有隧道二距新建隧道較近，受爆破振動(dòng)的影響較大。因此，在爆破位置周圍要加強(qiáng)支護(hù)，遠(yuǎn)離新建隧道部分的支護(hù)可在允許范圍內(nèi)適當(dāng)減少。

隨著到既有隧道拱頂距離的增加，地層沉降也會(huì)增加。在動(dòng)載作用下，既有隧道一拱頂最大沉降由3.5 mm增加到4.4 mm，既有隧道二拱頂最大沉降由5 mm增加到6.8 mm。其原因是缺少隧道管片的支護(hù)，使得地層沉降增加。在既有隧道拱頂處，由于襯砌支撐作用，其沉降較隧道兩側(cè)的小。因此，要加強(qiáng)距爆破位置近部分的支護(hù)，防止因結(jié)構(gòu)振動(dòng)造成的破壞，還要加強(qiáng)既有隧道周圍土體的強(qiáng)度，以減少地層沉降。

3.2.3 臺(tái)階法模擬爆破

將全斷面開挖爆破改為臺(tái)階法開挖爆破，會(huì)使爆破對(duì)隧道的影響降低。對(duì)既有隧道一斷面MC-6處的振動(dòng)速度進(jìn)行了檢測(cè)，監(jiān)測(cè)點(diǎn)的位置與MC-C的布置相同。分別采用全斷面爆破和臺(tái)階法爆破的方式進(jìn)行隧道開挖，監(jiān)測(cè)點(diǎn)的振動(dòng)速度如圖11所示。

圖11 監(jiān)測(cè)點(diǎn)振動(dòng)速度Fig. 11 The vibration speed at the monitoring point

從圖11中可以看出，全斷面爆破開挖的振動(dòng)速度最大值為5.2 cm/s，臺(tái)階法爆破開挖的振動(dòng)速度的最大值為3.7 cm/s，臺(tái)階法爆破比全斷面爆破的最大振動(dòng)速度減少了1.5 cm/s。隨著時(shí)間的增加，振動(dòng)呈現(xiàn)先增長(zhǎng)后減小的狀況，在7 ms左右出現(xiàn)峰值，然后不斷波動(dòng)，最后趨向于0 cm/s。臺(tái)階法爆破的曲線整體小于全斷面爆破的，臺(tái)階法爆破產(chǎn)生的影響比全斷面爆破產(chǎn)生的影響降低了28.8%，并且臺(tái)階法爆破產(chǎn)生的振動(dòng)較早地消失，施工更安全。

從圖11中還可以看出，較其他監(jiān)測(cè)點(diǎn)的振幅，位于既有隧道拱底監(jiān)測(cè)點(diǎn)6的振幅最大，監(jiān)測(cè)點(diǎn)1的振幅最小。在臺(tái)階法爆破第一次速度峰值中，監(jiān)測(cè)點(diǎn)6的速度為3.7 cm/s，位于拱頂?shù)谋O(jiān)測(cè)點(diǎn)1的振動(dòng)速度為1.05 cm/s，位于隧道中部其他監(jiān)測(cè)點(diǎn)的振動(dòng)速度處于這2個(gè)速度之間。表明：振動(dòng)的幅度與距爆破點(diǎn)的位置有關(guān)。距爆破點(diǎn)越近，受到爆破的影響越大。要加強(qiáng)爆破點(diǎn)附近的襯砌條件，同時(shí)選用的臺(tái)階法能夠降低爆破振動(dòng)所產(chǎn)生的影響。

4 結(jié)論

利用有限差分軟件Flac 3D，對(duì)隧道爆破開挖進(jìn)行數(shù)值模擬，并與實(shí)際數(shù)據(jù)對(duì)比，得到的結(jié)論為：

1) 新建隧道爆破載荷對(duì)水平位移的影響要大于其對(duì)豎向位移的影響。在爆破開挖的隧道中，要加強(qiáng)水平方向的支護(hù)措施，以減小爆破振動(dòng)所產(chǎn)生的影響。

2) 在爆破過(guò)程中，隧道應(yīng)力和應(yīng)變情況中間波動(dòng)量較大，對(duì)隧道的防護(hù)需留有一定的安全空間。

3) 臺(tái)階法爆破施工在一定程度上能夠降低結(jié)構(gòu)振動(dòng)幅度的28.8%，為爆破施工減振提供了一種方法。