趙建平 宋曉東

?

城市地鐵爆破開(kāi)挖對(duì)淺埋地下大直徑管道的安全影響研究1

趙建平 宋曉東

（中南大學(xué)資源與安全工程學(xué)院，長(zhǎng)沙 410083）

在城市淺埋地鐵爆破開(kāi)挖中，經(jīng)常遇到地下管網(wǎng)、涵洞等構(gòu)筑物，而爆破地震效應(yīng)對(duì)其影響范圍和程度的正確評(píng)價(jià)就顯得尤為重要。本文以長(zhǎng)沙地鐵爆破開(kāi)挖為例，以現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)為基礎(chǔ)，采用有限單元法，對(duì)爆破振動(dòng)下大直徑混凝土污水管道的動(dòng)力響應(yīng)、變形和動(dòng)應(yīng)力等進(jìn)行了計(jì)算，評(píng)價(jià)了爆破地震對(duì)管道的安全影響。研究表明：管道在控制爆破作用下是偏于安全的；爆心距是影響管道受到爆炸作用力影響大小的最主要因素；對(duì)于以實(shí)測(cè)數(shù)據(jù)為基礎(chǔ)，采用加速度激勵(lì)的時(shí)程分析方法以爆破振動(dòng)對(duì)埋地管道的影響進(jìn)行評(píng)價(jià)，是一種可行且較為精確的方法。

爆破振動(dòng) 有限元 時(shí)程分析 動(dòng)力響應(yīng)

引言

隨著爆破技術(shù)在城市地下工程施工中的逐步應(yīng)用，爆破振動(dòng)對(duì)埋設(shè)在其周?chē)慕Y(jié)構(gòu)或構(gòu)件造成的危害也越來(lái)越受到關(guān)注。地下爆破施工產(chǎn)生的應(yīng)力波，可以通過(guò)土層對(duì)其上部埋設(shè)的管線結(jié)構(gòu)造成危害，一旦發(fā)生破壞將造成嚴(yán)重的經(jīng)濟(jì)損失和社會(huì)影響。管道通常埋在地下，爆破動(dòng)力響應(yīng)無(wú)法通過(guò)實(shí)測(cè)數(shù)據(jù)得出，破壞后不易被發(fā)現(xiàn)和修復(fù)，容易造成更嚴(yán)重的次生事故，因此如何對(duì)這種影響進(jìn)行正確的評(píng)價(jià)就顯得尤為重要。《爆破安全規(guī)程（GB 6722-2003）》（中華人民共和國(guó)國(guó)家標(biāo)準(zhǔn)，2003）指出，評(píng)價(jià)爆破振動(dòng)對(duì)結(jié)構(gòu)安全的影響主要依據(jù)爆破地震動(dòng)速度，即通過(guò)實(shí)測(cè)數(shù)據(jù)與規(guī)程允許速度做對(duì)比進(jìn)行評(píng)價(jià)，但是僅僅根據(jù)最大振速來(lái)評(píng)價(jià)對(duì)結(jié)構(gòu)安全的影響顯然是不夠的。宋光明等（2000）和史秀志等（2010）根據(jù)爆破振動(dòng)實(shí)地監(jiān)測(cè)及質(zhì)點(diǎn)振速預(yù)測(cè)方程，結(jié)合現(xiàn)場(chǎng)損傷破壞調(diào)查，提出了一種爆破振動(dòng)危害動(dòng)態(tài)應(yīng)力比評(píng)價(jià)方法，然而這種方法對(duì)于地下埋地管線的研究還是有一定的局限性；婁建武等（2003）采用反應(yīng)譜方法研究爆破振動(dòng)信號(hào)的振動(dòng)速度和頻率特征，提出了采用反應(yīng)譜曲線積分值來(lái)評(píng)估爆破振動(dòng)破壞效應(yīng)，但沒(méi)有完全分析出結(jié)構(gòu)的動(dòng)力響應(yīng)。該類(lèi)問(wèn)題實(shí)際上是結(jié)構(gòu)動(dòng)力響應(yīng)問(wèn)題，要正確分析結(jié)構(gòu)動(dòng)力響應(yīng)最佳的分析方法是依據(jù)爆破地震實(shí)測(cè)數(shù)據(jù)，采用動(dòng)力分析的方法準(zhǔn)確分析結(jié)構(gòu)的動(dòng)力響應(yīng)、變形、動(dòng)應(yīng)力等（錢(qián)勝?lài)?guó)等，2006）。

為了更科學(xué)精確地評(píng)價(jià)地下隧道爆破對(duì)埋地污水管道的影響，本文運(yùn)用有限單元法，基于實(shí)測(cè)爆破地震動(dòng)數(shù)據(jù)，對(duì)埋地管道進(jìn)行了動(dòng)力時(shí)程分析。并以長(zhǎng)沙地鐵2號(hào)線暗挖爆破工程為例，評(píng)價(jià)了控制爆破工程中爆破地震動(dòng)對(duì)管道的安全影響。

1 結(jié)構(gòu)動(dòng)力時(shí)程分析原理

時(shí)程分析法是基于地震動(dòng)反應(yīng)數(shù)值仿真，模擬實(shí)際地震作用的方法，是反映建筑結(jié)構(gòu)地震響應(yīng)的最直接的方法（吳小峰等，2011）。在計(jì)算結(jié)構(gòu)動(dòng)力響應(yīng)時(shí)，從地基處輸入實(shí)測(cè)爆破地震波，從爆破地震動(dòng)初動(dòng)開(kāi)始，逐步對(duì)動(dòng)力平衡方程積分，直至地震動(dòng)波振幅為零的時(shí)刻停止，以求得地震動(dòng)過(guò)程中任一時(shí)刻的結(jié)構(gòu)變形和內(nèi)力，該方法能“再現(xiàn)”結(jié)構(gòu)在爆破地震動(dòng)時(shí)的變形狀況，是非常精確的動(dòng)力分析方法（陳穎，2009）。

把爆破振動(dòng)作用下的土-結(jié)構(gòu)作為相互作用的體系進(jìn)行爆破地震反應(yīng)分析（裴星洙，2006），根據(jù)達(dá)朗貝爾原理可以建立該體系在地震作用下任一時(shí)刻的動(dòng)力平衡方程：

式中，［］為系統(tǒng)的質(zhì)量矩陣；［］為阻尼矩陣；［］為系統(tǒng)的剛度矩陣；、、分別為系統(tǒng)的加速度、速度和位移向量；為輸入地震加速度向量。

對(duì)該運(yùn)動(dòng)方程求解時(shí)，采用時(shí)程逐步積分法，將整個(gè)地震過(guò)程分為若干微小時(shí)間段，并假設(shè)在各微小時(shí)間段內(nèi)體系是線性變化，通過(guò)積分可得到體系在整個(gè)地震過(guò)程中各時(shí)刻的運(yùn)動(dòng)狀態(tài)及其變化情況。同時(shí)還可將式（1）轉(zhuǎn)變成增量方程：

2 工程概況

分析實(shí)例為長(zhǎng)沙軌道交通2號(hào)線西湖公園站—溁灣鎮(zhèn)站區(qū)間隧道暗挖工程。區(qū)間隧道斷面形式采用分離式單洞斷面，區(qū)間暗挖段線間距13m，采用分離式單洞斷面形式，馬蹄形斷面結(jié)構(gòu)，采用上下臺(tái)階法施工，開(kāi)挖斷面寬6.30m，高6.60m，開(kāi)挖面積34.5m2。本文主要分析區(qū)間隧道爆破開(kāi)挖時(shí)，爆破振動(dòng)對(duì)埋地污水管道的影響，故選取位于污水管道上方編號(hào)XYCD-4號(hào)監(jiān)測(cè)點(diǎn)，對(duì)3月25日到4月27日監(jiān)測(cè)到的5次爆破實(shí)測(cè)數(shù)據(jù)進(jìn)行分析（以4月23日數(shù)據(jù)為例附圖）。實(shí)測(cè)數(shù)據(jù)見(jiàn)表1；實(shí)測(cè)振動(dòng)波形如圖1所示。

表1 隧道爆破振動(dòng)四號(hào)測(cè)點(diǎn)監(jiān)測(cè)數(shù)據(jù)表

續(xù)表

時(shí)間監(jiān)測(cè)點(diǎn)編號(hào)最大段量（kg）總藥量（kg）爆心距（m）垂直振動(dòng)速度（cm/s）振動(dòng)持續(xù)時(shí)間（s）振動(dòng)主頻（Hz） 4月12日XYCD-411.83555.90.2871.51820.02 4月22日XYCD-49.661.444.70.2161.92441.62 4月23日XYCD-410.050.820.80.6281.82014.65 4月27日XYCD-49.254.672.60.2131.99218.56

圖1 4月23日實(shí)測(cè)垂直方向振動(dòng)波形

圖2 物理模型圖

3 淺埋污水管動(dòng)力響應(yīng)分析

3.1 有限元模型的建立

爆破地震作用下，埋地管道的動(dòng)力響應(yīng)是土與管道結(jié)構(gòu)共同作用、相互影響的結(jié)果。對(duì)于污水管道，其埋在地下，爆破動(dòng)力響應(yīng)無(wú)法通過(guò)實(shí)測(cè)數(shù)據(jù)得出，而借助數(shù)值模擬仿真可分析地下結(jié)構(gòu)在爆破地震動(dòng)作用下的動(dòng)力響應(yīng)。

實(shí)際場(chǎng)地為三維，為便于分析，截取管道橫斷面作為分析對(duì)象，把三維問(wèn)題簡(jiǎn)化為二維平面應(yīng)變問(wèn)題。設(shè)模型中混凝土污水管道管芯外徑為2.40m，厚320mm，埋深4m，上邊界為自由地面，下邊界距上表面18.4m。結(jié)構(gòu)所處場(chǎng)地上表面距離地下10m為土體，往下8m為基巖，物理模型如圖2所示?；诰鶆蚧碚摷坝邢拊椒ǎㄩZ曉鵬等，2011），在ANSYS建模時(shí)，鋼筒按照與混凝土彈性模量之比，將鋼筒等效為混凝土。場(chǎng)地土、混凝土管芯、基巖均采用平面單元PLANE42模擬，預(yù)應(yīng)力鋼絲采用LINK8模擬，采用初始應(yīng)變法施加預(yù)應(yīng)力。預(yù)應(yīng)力鋼絲與混凝土采用COMBIN14單元連接。根據(jù)工程項(xiàng)目提供的地質(zhì)資料并查閱相關(guān)資料，確定模型中材料特性如表2所示。邊界條件底邊和兩個(gè)側(cè)邊界均采用粘彈性人工邊界條件處理（劉晶波等，2006；2007；劉漢龍等，2001）。根據(jù)經(jīng)驗(yàn)做法，通常取地下結(jié)構(gòu)寬度的5—10倍范圍作為計(jì)算區(qū)域（杜英杰，2008），有限元模型中左右兩側(cè)邊界取離地下管道25m，本文計(jì)算的有限元模型共23087個(gè)單元，23522個(gè)節(jié)點(diǎn)。ANSYS模型及網(wǎng)格劃分如圖3所示。

圖3 有限元模型及網(wǎng)格劃分圖

表2 材料參數(shù)

3.2 動(dòng)力響應(yīng)分析結(jié)果

實(shí)測(cè)數(shù)據(jù)表明，本次爆破以豎向振動(dòng)為主，故輸入地震波采用實(shí)測(cè)豎向地震動(dòng)。本文取離污水管道正上方約5m處四號(hào)監(jiān)測(cè)點(diǎn)實(shí)測(cè)豎直方向波形為例進(jìn)行分析，通過(guò)對(duì)實(shí)測(cè)速度波形進(jìn)行一元微分，得到了加速度時(shí)程曲線，然后采用加速度激勵(lì)的時(shí)程方法進(jìn)行加載求解。計(jì)算結(jié)果見(jiàn)圖4—圖6。

圖4 污水管道變形圖

計(jì)算結(jié)果表明，爆破振動(dòng)對(duì)污水管道產(chǎn)生的最大拉應(yīng)力為370.49Pa，主要集中在管道上下兩側(cè)的內(nèi)側(cè)；爆破振動(dòng)對(duì)于整個(gè)爆破模型產(chǎn)生的整體最大位移為1.476mm；爆破振動(dòng)對(duì)于污水管道產(chǎn)生的最大位移為1.213mm，幾乎可以忽略不計(jì)。另外，取污水管道頂部2號(hào)節(jié)點(diǎn)做時(shí)間歷程后處理，分別得到了該節(jié)點(diǎn)的豎向位移、速度、加速度時(shí)程曲線，具體如圖7—圖10所示。

圖7 2號(hào)節(jié)點(diǎn)位置

圖8 2號(hào)節(jié)點(diǎn)豎直方向位移時(shí)程曲線

圖9 2號(hào)節(jié)點(diǎn)豎直方向速度時(shí)程曲線

圖10 2號(hào)節(jié)點(diǎn)豎直方向加速度時(shí)程曲線

從圖9中可以得到，2號(hào)節(jié)點(diǎn)最大速度為0.625cm/s，這與地表4號(hào)測(cè)點(diǎn)豎直方向?qū)崪y(cè)數(shù)據(jù)0.628cm/s相差不大。同時(shí)，分別輸入其他幾次爆破監(jiān)測(cè)實(shí)測(cè)波形，還可分別得到污水管道垂直方向的動(dòng)力響應(yīng)位移、響應(yīng)應(yīng)力和2號(hào)節(jié)點(diǎn)的響應(yīng)速度，結(jié)果見(jiàn)表3。

表3 時(shí)程分析結(jié)果

3.3 結(jié)果分析

從表3可知：①時(shí)程分析得到的最大振速與實(shí)測(cè)數(shù)據(jù)相差很小，都遠(yuǎn)小于爆破規(guī)程中規(guī)定的4cm/s。時(shí)程分析數(shù)值偏小，是因?yàn)樗奶?hào)測(cè)點(diǎn)比地下污水管道離爆源更近，但與實(shí)測(cè)數(shù)據(jù)很接近，說(shuō)明時(shí)程分析結(jié)果基本可靠。筆者建議對(duì)地下污水管道的監(jiān)測(cè)可以采用管道上方附近地面實(shí)測(cè)數(shù)據(jù)作為重要參考依據(jù)。②爆破振動(dòng)對(duì)污水管道產(chǎn)生的響應(yīng)位移非常小，最小位移為0.179mm，爆心距為20m時(shí)最大動(dòng)力響應(yīng)位移才1.213mm，幾乎可以忽略不計(jì)。③爆破振動(dòng)對(duì)管道產(chǎn)生的等效應(yīng)力最大值為680Pa，遠(yuǎn)小于混凝土管道MPa級(jí)的抗拉強(qiáng)度，也遠(yuǎn)小于按《室外給水排水和煤氣熱力工程抗震設(shè)計(jì)規(guī)范》規(guī)定方法的計(jì)算結(jié)果3MPa（薛景宏等，2005）。可以說(shuō)，地鐵隧道爆破振動(dòng)對(duì)管道影響甚小，管道基本處于安全穩(wěn)定狀態(tài)。

另外從表3還可以得出，并不是響應(yīng)速度越大，響應(yīng)位移就越大，這還與爆破振動(dòng)持續(xù)時(shí)間有關(guān)。爆破振動(dòng)持續(xù)時(shí)間越長(zhǎng)，速度越大，則產(chǎn)生的響應(yīng)位移就越大。

比較裝藥量和爆心距對(duì)管道動(dòng)力響應(yīng)的影響時(shí)，當(dāng)最大段藥量相差不大，隨著爆心距的減小，表3中第一次和第四次監(jiān)測(cè)動(dòng)力響應(yīng)速度和位移增幅較大；在爆心距相差不大時(shí)，隨著藥量的增加，第一次和第五次監(jiān)測(cè)動(dòng)力響應(yīng)速度和位移增幅并不明顯，位移幾乎無(wú)變化。時(shí)程分析表明，管道的受力及運(yùn)動(dòng)情況受爆心距的影響遠(yuǎn)大于裝藥量大小的影響。筆者建議在城市地下爆破開(kāi)挖工程中，如遇到管道的安全問(wèn)題，應(yīng)該首先考慮管道到爆源的距離大小，爆心距是影響管道受到爆炸作用力大小的主要因素。

4 結(jié)論

（1）對(duì)地下污水管道進(jìn)行的動(dòng)力時(shí)程分析表明，數(shù)值計(jì)算的管道最大振速與實(shí)測(cè)數(shù)據(jù)相差不大，且都遠(yuǎn)小于爆破規(guī)程中規(guī)定的4cm/s。爆破對(duì)管道產(chǎn)生的動(dòng)力響應(yīng)位移和應(yīng)力也都非常小，幾乎可以忽略不計(jì)。在地鐵控制爆破施工中，爆破振動(dòng)對(duì)埋地管道影響甚小，管道是安全的。建議對(duì)于地下污水管道的監(jiān)測(cè)可以采用管道上方附近地面實(shí)測(cè)數(shù)據(jù)作為重要參考依據(jù)。

（2）地鐵控制爆破施工中，爆破振動(dòng)持續(xù)時(shí)間也應(yīng)當(dāng)引起足夠重視；此外，管道的受力及運(yùn)動(dòng)情況受爆心距的影響遠(yuǎn)大于受裝藥量大小的影響。建議在進(jìn)行地下爆破時(shí)，對(duì)于管道的安全問(wèn)題應(yīng)該首先考慮管道距爆源的距離大小，爆心距是影響管道受到爆破作用力影響大小的主要因素。

（3）以實(shí)測(cè)數(shù)據(jù)為基礎(chǔ)，采用實(shí)測(cè)數(shù)據(jù)結(jié)合加速度激勵(lì)的時(shí)程分析有限元方法，就爆破地震對(duì)埋地污水管道的影響進(jìn)行了評(píng)價(jià)，計(jì)算出了實(shí)際變形、動(dòng)應(yīng)力等。結(jié)果表明本文的計(jì)算方法是有效的和較為準(zhǔn)確的。

陳穎，2009. 某砌體房屋爆破地震動(dòng)的實(shí)測(cè)及結(jié)構(gòu)動(dòng)力反應(yīng)時(shí)程分析. 工程抗震與加固改造，32（6）：117—120.

杜英杰，2008. 埋地預(yù)應(yīng)力鋼筒混凝土管的地震響應(yīng)分析. 天津：天津大學(xué)，39.

賈要紅，2011. 地鐵換乘車(chē)站多層結(jié)構(gòu)地震響應(yīng)分析. 西安：西安科技大學(xué)，13—14.

婁建武，龍?jiān)矗?003. 基于反應(yīng)譜值分析的爆破震動(dòng)破壞評(píng)估研究. 爆炸與沖擊，23（1）：41—46.

劉晶波，谷音等，2006. 一致粘彈性人工邊界及粘彈性邊界單元. 巖土工程學(xué)報(bào)，28（9）：1070—1075.

劉晶波，杜義欣等，2007. 粘彈性人工邊界及地震動(dòng)輸入在通用有限元軟件中的實(shí)現(xiàn). 防災(zāi)減災(zāi)工程學(xué)報(bào)，27（增刊）：37—42.

劉漢龍，高玉峰等，2001. 邊界條件對(duì)土層粘彈性地震反應(yīng)的影響. 巖土力學(xué)，22（4）：408—412.

裴星洙，2006. 高層建筑結(jié)構(gòu)地震響應(yīng)的時(shí)程分析法. 中國(guó)水利水電出版社，187—196.

錢(qián)勝?lài)?guó)，陳玲玲等，2006. 爆破地震工程結(jié)構(gòu)動(dòng)力分析. 北京：中國(guó)水利水電出版社，36—39.

宋光明，史秀志，陳壽如，2000. 露天礦邊坡爆破振動(dòng)破壞判據(jù)新方法及其應(yīng)用. 中南工業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），31（6）：485—488.

史秀志，周健等，2010. 爆破振動(dòng)對(duì)民房破壞效應(yīng)預(yù)測(cè)的BDA模型及應(yīng)用. 振動(dòng)與沖擊，29（7）：60—65.

吳小峰，孫啟國(guó)等，2011. 抗震分析反應(yīng)譜法和時(shí)程分析法數(shù)值仿真比較. 西北地震學(xué)報(bào)，33（3）：275—279.

薛景宏，逮廣東，2005. 埋地管道動(dòng)力反應(yīng)時(shí)程分析結(jié)果與規(guī)范比較. 油氣田地面工程，24（11）：11—12.

閆曉鵬，牛衛(wèi)晶等，2011. 基于均勻化理論的混凝土等效彈性模量的數(shù)值模擬. 太原理工大學(xué)學(xué)報(bào)，42（2）： 212—214.

Security Impact of Blasting Excavation of Urban Subway on Shallow Buried Large Diameter Pipeline

Zhao Jianping and Song Xiaodong

(School of Resources and Safety Engineering, Central South University, Changsha 410083, China)

In the process of shallow buried underground blasting excavation, there often exists underground pipe network, culverts and other buildings. Therefore, the influence range and degree of accuracy evaluation for the blasting is particularly important. Taking Changsha subway tunnel blasting construction as an example, we used finite element method based on spot-measuring data to calculate the pipeline dynamic response on buried sewage pipeline under tunnel blasting vibration as well as deformation and stress, and to evaluate the safety effect of blasting vibration on pipeline. The results show that the pipelines are basically safe. Blasting distance is the main factor to pipeline affected by explosive force size. Finally, time-history analysis method of acceleration excitation is a feasible and accurate method for evaluating the impact on buried pipeline and blasting vibration.

Blasting vibration; Finite element; Time-history analysis; Dynamic response

國(guó)家自然科學(xué)基金青年基金（51104178）

2014-12-19

趙建平，生于1977年。博士，副教授。研究方向：工程爆破、巖土工程。E-mail：jpzcsu@126.com