王海亮，張　 亮，張祖遠

(1.山東科技大學(xué) 礦山災(zāi)害預(yù)防控制國家重點實驗室培育基地，山東 青島 266590;2.上海建工二建集團有限公司南昌分公司，江西 南昌 330000)

城市隧道爆破震動影響因素的重要度研究

王海亮1，張 亮1，張祖遠2

(1.山東科技大學(xué) 礦山災(zāi)害預(yù)防控制國家重點實驗室培育基地，山東 青島 266590;2.上海建工二建集團有限公司南昌分公司，江西 南昌 330000)

摘要：爆破震動對建(構(gòu))筑物會產(chǎn)生不同程度的影響，有效降低爆破震動的有害效應(yīng)是控制爆破震動的首要問題。本研究對影響爆破震動強度的相關(guān)因素進行分類，采用灰色關(guān)聯(lián)分析方法，以青島市兩地鐵工程部分實測數(shù)據(jù)為樣本，結(jié)合MATLAB程序編程計算，得到應(yīng)優(yōu)先選擇的控制因素按重要度從大到小依次為：自由面數(shù)目、段最小間隔時間、總藥量、地形地質(zhì)條件、雷管段數(shù)、孔網(wǎng)參數(shù)、爆破持續(xù)時間、單段最大起爆藥量、起爆方向。

關(guān)鍵詞：隧道;爆破震動;重要度;灰色關(guān)聯(lián)分析

爆破在城市地下軌道交通工程中仍占據(jù)較為重要的地位。因城市周圍環(huán)境條件限制，爆破震動控制的要求較高[1]，有效降低爆破震動的影響成為城市隧道工程迫切需要解決的問題。由于城市隧道工程爆破實際施工中的有害效應(yīng)較多[2]，影響爆破震動強度的因素呈離散關(guān)系，同時因受各種客觀因素的影響，難以全面統(tǒng)計到各種影響爆破震動強度因素的數(shù)據(jù)。

目前針對離散系統(tǒng)較多的采用BP(back propagation)神經(jīng)網(wǎng)絡(luò)、GRA(grey relational analysis)灰色度關(guān)聯(lián)等研究方法。沈薇等[3]通過對灰色關(guān)聯(lián)度矩陣的分析，指出臺階爆破爆源與待測點的高差是影響爆破震動參數(shù)的準優(yōu)因素；谷任國等[4]分析了水下爆破震動影響因素的主次關(guān)系，指出不同爆破地區(qū)及爆破條件情況下，爆破參數(shù)影響程度不同；趙同彬等[5]根據(jù)現(xiàn)場爆破實測數(shù)據(jù)采用的灰局勢方法優(yōu)選巖巷掏槽爆破方案，為巖巷掏槽爆破方案優(yōu)選提供了參考；黃成林等[6]采用GRA方法分析了蒙庫鐵礦露天邊坡爆破震動的主要影響因素，指出主要因素為炸藥總量、單段藥量等；楊珊等[7]以廈門某隧道工程為背景，通過灰色關(guān)聯(lián)計算得出爆破震動影響因素的主次順序，確定爆心距為準優(yōu)因素，同時根據(jù)可控準優(yōu)參數(shù)對爆破參數(shù)進行了優(yōu)化。

GRA灰色關(guān)聯(lián)分析法作為影響因素重要度研究量化的重要手段，其在分析中要收集的樣本量較少，適合爆破震動影響研究。本研究采用灰色關(guān)聯(lián)法來分析城市隧道工程中影響爆破震動強度可控因素的重要度，為制定科學(xué)合理的爆破震動控制措施、控制爆破震動有害效應(yīng)提供參考。

1隧道爆破震動強度的影響因素

影響隧道爆破震動強度的因素有很多[8]，大體可分為兩類：第一類是爆源因素，包括炸藥爆速、單段最大起爆藥量、總藥量、雷管段數(shù)、段最小間隔時間、孔網(wǎng)參數(shù)、裝藥結(jié)構(gòu)、自由面數(shù)目、爆破方向和爆破持續(xù)時間等因素；第二類為傳播途徑方面，其中包含了高程、爆源與監(jiān)控點距離、工程地質(zhì)等因素。

采取減振措施時不加以區(qū)分，把每一個因素都進行控制不太可行[9]。因此，首先應(yīng)當把影響因素按能否采取控制措施分為可控因素和不可控因素。分析可知，爆源因素全部為可控因素，而傳播途徑因素除地形地質(zhì)條件外，都為不可控因素。針對諸多可控影響因素，需要通過分析計算進行排序。進而制定針對性安全技術(shù)措施，是實現(xiàn)城市地下隧道工程爆破震動控制的有效途徑之一。

2影響因素重要度計算

2.1GRA運算方法

GRA關(guān)聯(lián)度反映不同因素間相互發(fā)展變化的異同度，其數(shù)值表征了異同度的大小。將影響爆破震動效應(yīng)的各個因素間的灰色絕對關(guān)聯(lián)度數(shù)值按遞減排序，來表征各影響因素重要程度的次序。其運算步驟具體為[10-12]：

1)確定系統(tǒng)特征變量數(shù)據(jù)序列和相關(guān)可控因素數(shù)據(jù)序列

假定xj(k)為第j個關(guān)聯(lián)因素變化值在第k次試驗時監(jiān)測的樣本，假定yi(k)為第i個系統(tǒng)特征值在第k次試驗時監(jiān)測的樣本，表示為：

Yi=[yi(1),yi(2),…，yi(k),…yi(n)],i=1,2,3,…m,

(1)

Xj=[xj(1),xj(2),…，xj(k),…xj(n)],j=1,2,3，…l。

(2)

其中：Yi是第i個系統(tǒng)特征因素變量；Xj是第j個相關(guān)因素變量；n值為樣本的容量值；m表示系統(tǒng)特征因素變量的數(shù)目；l為相關(guān)因素變化值的數(shù)目。

2)正相關(guān)因素與負相關(guān)因素轉(zhuǎn)化的轉(zhuǎn)換

(3)

式中，D0是對稱化算子，Xj是對稱化像。

若Xj(k)與Yi(k)為負相關(guān)關(guān)系，就可以通過算子D0的作用，使Xj(k)的對稱化像XjD0與Yi(k)有正相關(guān)關(guān)系。

3)無量綱化轉(zhuǎn)換

主要是借助算子作用，轉(zhuǎn)換為數(shù)量級上近似的無量綱數(shù)值，同時把負相關(guān)因素變成正相關(guān)。普遍采用的方法包含初值化法、區(qū)間值化法、均值法等計算方法。

進行相關(guān)影響因素數(shù)據(jù)Xj無量綱化轉(zhuǎn)換

(4)

式中，D1是初值化算子，Xj是原象，XjD1是Xj在初值化算子下的像，簡稱初值像。

(5)

其中，D2是初值化算子，Xj是原象，XjD2是Xj在初值化算子下的像，簡稱均值像。

(6)

式中，D3是初值化算子，Xj是原象，XjD3是Xj在初值化算子下的像，簡稱區(qū)間值像。

4)差序列、最大差、最小差的計算

在進行無量綱轉(zhuǎn)換的基礎(chǔ)上，系統(tǒng)特征變量數(shù)據(jù)序列表示為Xo(k)，相關(guān)因素數(shù)據(jù)序列表示為xj(k)，那么其差序列表示為：

Δ0i=‖x0(k)-xi(k)‖。

(7)

式中，i=1,2,3，…m,k=1,2,3，…n。則其最大差和最小差分別為：

(8)

5)計算灰色關(guān)聯(lián)系數(shù)

計算灰色關(guān)聯(lián)系數(shù)(x0(k),xj(k))的數(shù)值，其計算公式為：

(9)

6)求得灰色關(guān)聯(lián)度

在計算灰色關(guān)聯(lián)度時，一般采用平權(quán)處理，式(10)正是運用平權(quán)處理得出權(quán)重值，既而得出加權(quán)灰色關(guān)聯(lián)度。

(10)

7)將關(guān)聯(lián)度按照數(shù)值大小進行前后排序。

8)優(yōu)點評價

根據(jù)前七個步驟可以計算得出各因素間的關(guān)聯(lián)度，但考慮相關(guān)參數(shù)數(shù)列不局限于單一情況，其關(guān)聯(lián)度可以組成一個矩陣E。

(11)

若存在h，j∈{1,2,…,n}，滿足：

εih≥εij,i=1,2,3,…,m。

(12)

則稱相關(guān)因素Xh優(yōu)于Xj，即相關(guān)因素Xh的重要度要大于Xj。若滿足：

(13)

則稱相關(guān)因素Xh準優(yōu)于Xj，即相關(guān)因素Xh的重要度要大于Xj。最優(yōu)因素為Xh，其判定條件為j∈{1,2,3,…,n}，使得j=1,2,3,…,n，j≠n，Xh都較Xj優(yōu)。若存在j∈{1,2,3,…,n}滿足對任意的j=1,2,3,…,n，j≠n，Xh較Xj優(yōu)，即Xh是準優(yōu)因素，而且Xh重要度為最高。

圖1　基于GRA的流程圖

2.2程序基本步驟

在對GRA的計算中采用計算機編程的方式，以便快速準確地得出結(jié)果。綜合分析效率、可視化、科學(xué)計算等因素，決定采用MATLAB軟件進行程序語言的編寫，具體參考流程圖1。

3隧道爆破震動影響因素的GRA

在明確影響因素重要度計算方式和流程后，結(jié)合青島市地鐵一期工程2號線某車站1號風(fēng)道側(cè)穿民房段和3號線下穿某醫(yī)院段為實際工程背景，選取部分典型實測數(shù)據(jù)進行灰色關(guān)聯(lián)分析。兩個工程背景的地質(zhì)條件基本相似，掌子面都處于強風(fēng)化及中風(fēng)化花崗巖層。

3.1系統(tǒng)特征變量和相關(guān)因素變量的選擇

目前，在國內(nèi)外研究中，振動影響依據(jù)振動速度、主頻、持續(xù)時間等，根據(jù)我國《爆破安全規(guī)程》(GB6722—2014)中爆破震動對地面建筑物影響的判斷標準主要還是最大振動速度和主振頻率兩大指標[13]。因此采用質(zhì)點峰值振動速度和主振頻率為系統(tǒng)特征變量，記為Yi，(i=1,2)，Y1表示最大振動速度，Y2即爆破震動的主振頻率。

Xj表示作用于爆破最大振動速度和爆破主振頻率的相關(guān)因素變量，(j=1,2,3,…l)。對于隧道爆破，相關(guān)影響因素主要分為爆源因素和傳播途徑因素兩大類。國內(nèi)外較多研究人員把羅列的影響因素都定為相關(guān)因素變量[14-16]。本研究考慮到實際控制時能否采取措施，分為可控影響因素和不可控影響因素。從爆破震動控制的實際角度出發(fā)，只定義可控影響因素為相關(guān)因素變量。可控影響因素包括炸藥爆速、單段最大藥量值、總起爆藥量、雷管段別數(shù)目、雷管段別間隔時間最小值、孔網(wǎng)爆破參數(shù)、裝藥結(jié)構(gòu)、自由面數(shù)目、爆破方向、爆破持續(xù)時間、地形地質(zhì)條件等。由于兩個工程都使用同一公司生產(chǎn)的2號巖石乳化炸藥，炸藥爆速是不變的，因此不予考慮炸藥爆速。裝藥結(jié)構(gòu)也基本不變，不予考慮。其余9個可控因素全部定義為相關(guān)因素變量，即l=9。單段最大起爆藥量為X1，總藥量為X2，雷管段數(shù)為X3，段最小間隔時間為X4，孔網(wǎng)參數(shù)為X5，自由面數(shù)目為X6，爆破方向為X7，爆破持續(xù)時間為X8，地形地質(zhì)條件為X9。

圖2　測點布置示意圖

孔網(wǎng)參數(shù)包含的內(nèi)容較多，包括孔深、孔間排距、布孔形式等，孔深與裝藥量成正比，可以在總裝藥量和單段最大起爆藥量的量化上體現(xiàn)[3]。孔間排距則與鉆眼密度有關(guān)，可用比鉆眼數(shù)來進行量化。布孔形式呈多樣性，難以量化，暫不考慮。因此，本章將孔網(wǎng)參數(shù)近似量化為比鉆眼數(shù)。

量化爆破方向時，將中心定在裝藥炮孔中心處，把炸藥起爆時最小抵抗線方向確定為基準線。從基準線到測點的角度為爆破方向，按順時針旋轉(zhuǎn)確定角度，大小以弧度計。如圖2所示，爆破方向β的取值范圍為0≤β≤6.28。

關(guān)于地形地質(zhì)條件的量化較難實現(xiàn)，嘗試性地根據(jù)地質(zhì)體的破碎程度及是否存在空腔考慮。將地質(zhì)條件作為作用于爆破地震波傳播影響因素變量的文獻并不多見，對地質(zhì)條件采用量化指標，并入相關(guān)因素變量中，來取得較為吻合工程的計算結(jié)論。根據(jù)地質(zhì)條件得到的具體指標0.2～1.0分五個指標，幅度差為0.2。例如在開挖斷面下臺階中，上臺階已超前開挖，地形地質(zhì)條件相當于發(fā)生了改變，量化為0.4。若采用CRD法施工時，CRD法上層開挖完成后，下層開挖時地形地質(zhì)量化為0.6。表1是青島市地鐵兩個工程背景隧道爆破時隨機抽取的部分實測數(shù)據(jù)樣本。

表1　爆破震動實測數(shù)據(jù)樣本

3.2關(guān)聯(lián)度的統(tǒng)計運算

根據(jù)表1青島市兩工程中的實測數(shù)據(jù)樣本，在前述分析編寫的MATLAB程序編輯器中運算。無量綱轉(zhuǎn)換采用均值法，選取0.5的分辨率系數(shù)代入。由于峰值振動速度與最小間隔時間和自由面數(shù)呈負相關(guān)，而主振頻率與總藥量呈負相關(guān)[17-18]，在進行灰色關(guān)聯(lián)度計算前，應(yīng)先采用軸對稱方法將負相關(guān)因素轉(zhuǎn)化為正相關(guān)因素。得到隧道爆破震動影響因素灰色關(guān)聯(lián)度計算結(jié)果如表2所示。

3.3計算結(jié)果分析

由表2的計算結(jié)果，可得到各相關(guān)影響因素重要度順序圖(圖3)。由表2和圖3的結(jié)果對相關(guān)因素進行灰色關(guān)聯(lián)度排序和優(yōu)勢分析，可得到如下結(jié)論：

1)9個可控影響因素對峰值振動速度的影響主次順序為：單段最大藥量值、爆破自由面的數(shù)量、雷管段別間隔時間最小值、總藥量、地形地質(zhì)條件、雷管段數(shù)、孔網(wǎng)參數(shù)、爆破持續(xù)時間、起爆方向。

2)9個可控影響因素對主振頻率的影響程度從大到小依次為：自由面數(shù)目、總藥量、段最小間隔時間、雷管段數(shù)、孔網(wǎng)參數(shù)、地形地質(zhì)條件、爆破持續(xù)時間、起爆方向、單段最大起爆藥量。

表2　隧道爆破震動影響因素GRA結(jié)果

圖3　影響因素重要度順序

3)9個可控影響因素中沒有最優(yōu)因素，準優(yōu)因素為自由面數(shù)目，其次為段最小間隔時間和總藥量。因此在進行隧道爆破施工時，若需要考慮爆破震動對周圍建(構(gòu))筑物的影響，最有效地控制爆破震動有害效應(yīng)的措施是充分利用和創(chuàng)造盡可能多的自由面、并嚴格控制段最小間隔時間，依次按圖3中的重要度順序采取振動控制措施。

隧道爆破施工時，單段最大起爆藥量和自由面數(shù)目對爆破峰值振速的影響最大，自由面數(shù)目和總藥量對爆破震動主頻的影響最大，若同時考慮對峰值振速和振動主頻的影響，則自由面數(shù)目和段最小間隔時差影響最大。實際施工時，可根據(jù)現(xiàn)場實測數(shù)據(jù)采取相應(yīng)的措施，使爆破震動控制的手段更具有針對性。

4結(jié)論

利用灰色關(guān)聯(lián)分析(GRA)方法，得到爆破震動強度可控影響因素的重要度排序，主要結(jié)論為：

1)隧道爆破施工時，可控影響因素對峰值振動速度的影響主次順序為：單段最大藥量值、爆破自由面的數(shù)量、雷管段別間隔時間最小值、總藥量、地形地質(zhì)條件、雷管段數(shù)、孔網(wǎng)參數(shù)、爆破持續(xù)時間、起爆方向。

2)可控影響因素對主振頻率的影響程度從大到小依次為：自由面數(shù)目、總藥量、段最小間隔時間、雷管段數(shù)、孔網(wǎng)參數(shù)、地形地質(zhì)條件、爆破持續(xù)時間、起爆方向、單段最大起爆藥量。

3)綜合考慮隧道爆破峰值振動速度和主振頻率對周邊建(構(gòu))筑物的有害效應(yīng)時，應(yīng)優(yōu)先選擇的控制因素按重要度從大到小依次為：自由面數(shù)目、段最小間隔時間、總藥量、地形地質(zhì)條件、雷管段數(shù)、孔網(wǎng)參數(shù)、爆破持續(xù)時間、單段最大起爆藥量、起爆方向。

本課題僅針對特定爆破工程分析，應(yīng)針對不同條件、根據(jù)具體情況進行深入研究，以確保得到更符合實際工程條件的爆破震動強度優(yōu)勢影響因素。

參考文獻

[1]孫松,王海亮,彭志燕.豎井爆破控制技術(shù)[J].山東科技大學(xué)學(xué)報(自然科學(xué)版),2015,34(1):106-110.

SUN Song,WANG Hailiang,PENG Zhiyan.The control of flyrock in shaft blasting[J].Journal of Shandong University of Science and Technology (Natural Science),2015,34(1):106-110.

[2]王海亮.工程爆破[M].北京:中國鐵道出版社,2008:179-185.

[3]沈蔚,徐全軍,黃文華.用灰關(guān)聯(lián)分析法確定爆破震動參數(shù)的主要影響因素[J].工程爆破,2000,6(4):6-13.

SHEN Wei,XU Quanjun,HUANG Wenhua,et al.Determination of main influencing factors on blasting vibaration parameters by Grey correlation analysis[J].Engineering Blasting,2000,6(4):6-13.

[4]谷任國,周傳波.水下爆破震動影響因素的灰關(guān)聯(lián)分析[J].爆破,2001,18(4):15-17.

GU Renguo,ZHOU Chuanbo.The influential factors of blasting vibration underwater by Grey correlation analysis[J].Blasting,2001,18(4):15-17.

[5]趙同彬,馬植勝,佟養(yǎng)貞,等.優(yōu)選巖巷掏槽爆破方案的灰局勢方法[J].山東科技大學(xué)學(xué)報(自然科學(xué)版),2004,23(4):30-32.

ZHAO Tongbin,MA Zhisheng,TONG Yangzhen,et al.The Grey situation method of optimizing cutting blasting plan of rock roadways[J].Journal of Shandong University of Science and Technology (Natural Science),2004,23(4):30-32.

[6]黃成林,陳建平,羅學(xué)東,等.露天礦邊坡爆破震動影響因素的灰關(guān)聯(lián)分析[J].爆破,2010,27(2):96-98.

HUANG Chenglin,CHEN Jianping,LUO Xuedong,et al.Grey relational analysis on influence factors of blasting vibration in open-pit slope[J].Blasting,2010,27(2):96-98.

[7]楊珊,陳建宏,郭宏斌,等.隧道爆破震動影響因素的灰色關(guān)聯(lián)分析[J].中國安全科學(xué)學(xué)報,21(6):65-71.

YANG Shan,CHEN Jianhong,GUO Hongbin,et al.Grey relational analysis on influence factors of tunnel blasting vibration[J].China Safety Science Journal,21(6):65-71.

[8]周宜,王海亮,張祖遠,等.城市淺埋隧道下穿建筑物的微振動爆破技術(shù)研究[J].隧道建設(shè),2015,35(1):89-92.

ZHOU Yi,WANG Hailiang,ZHANG Zuyuan,et al.Study on micro-vibration blasting technology for shallow tunnels crossing underneath existing buildings in urban area[J].Tunnel Construction,2015,35(1):89-92.

[9]范孝鋒,周傳波,陳國平.爆破震動影響因素的灰關(guān)聯(lián)分析[J].爆破,2005,22(6):100-102.

FAN Xiaofeng,ZHOU Chuanbo,CHEN Guoping.The influentialfactorsofblasting vibration by Grey correlation analysis[J].Blasting,2005,22(6):100-102.

[10]肖新平,宋中民,李峰.灰技術(shù)基礎(chǔ)及其應(yīng)用[M].北京:科學(xué)出版社,2005:26-50.

[11]劉思峰,郭天榜.灰色系統(tǒng)理論及其應(yīng)用[M].4版.北京:科學(xué)出版社,2008:44-71.

[12]劉思峰,郭天榜,黨耀國,等.灰色系統(tǒng)理論及其應(yīng)用[M].2版.北京:科學(xué)出版社,1999:256-261.

[13]中華人民共和國國家標準.GB6722-2014,爆破安全規(guī)程[S].

[14]羅正.城市淺埋大跨隧道爆破震動監(jiān)測與控制技術(shù)研究[D].長沙：中南大學(xué),2011:53-65.

[15]李付勝.城市隧道掘進爆破數(shù)值模擬及震動效應(yīng)控制研究[D].重慶：重慶大學(xué),2008:53-63.

[16]趙新濤.城市巖體開挖爆破震動效應(yīng)及安全控制研究[D].重慶：重慶大學(xué),2010:111-130.

[17]張奇,白春華,劉慶明.爆炸地震波頻譜特性研究[J].北京理工大學(xué)學(xué)報(自然科學(xué)版),1999,19(3):306-308.

ZHANG Qi,BAI Chunhua,LIU Qingming.Frequency spectrum characteristics of explosion earthquake wave[J].Journal of Beijing Institute of Technology (Natural Science Edition),1999,19(3):306-308.

[18]李孝林,穆太升,董鑫,等.頻率在爆破震害中的作用及其影響因素分析[J].工程爆破,2001,7(3):15-18.

LI Xiaolin,MU Taisheng,DONG Xin,et al.Role of frequency in harm of blasting vibration and analysis on its influence factors[J].Engineering Blasting,2001,7(3):15-18.

(責(zé)任編輯：呂海亮)

Importance Degree of Factors Influencing Tunnel Blasting Vibration

WANG Hailiang1,ZHANG Liang1,ZHANG Zuyuan2

(1.State Key Laboratory of Mining Disaster Prevention and Control Co-founded by Shandong Province and the Ministry of Science and Technology,Shandong University of Science and Technology,Qingdao,Shandong 266590,China;2.Nanchang Branch of Shanghai Construction Group Erjian Group Co.Ltd,Nanchang,Jiangxi 330000,China)

Abstract:Since blasting vibration has different influences on buildings,to effectively reduce its harmful effects becomes the primary problem of blasting vibration control.The paper classified the relevant factors influencing blasting vibration strength and obtained the controlling factors that should be selected primarily by using the Grey relational analysis (GRA) method and taking the measured data of Qingdao subway project as samples.After MATLAB programming calculation,the sequence of these factors in decreasing order according to importance degree was determined.They were number of free surfaces,minimum time interval,total explosive charge,topographic and geologic conditions,number of detonators,hole pattern parameters,blasting duration,largest single segment initiation dose,and initiation direction.

Key words:tunnel;blasting vibration;importance degree;Grey relational analysis

收稿日期：2015-11-09

基金項目：國家自然科學(xué)基金項目(10672091)；青島巿建委科技項目(JK2013-4)

作者簡介：王海亮(1963—)，男，河北石家莊人，教授，博士生導(dǎo)師，主要從事工程爆破、地下工程、安全評價理論研究， E-mail:tlgcbp@263.net E-mail：kedazhangliang@126.com

中圖分類號：TB41

文獻標志碼：A

文章編號：1672-3767(2016)02-0073-06

張亮(1990—)，男，山東泰安人，碩士研究生，主要從事隧道工程、災(zāi)害預(yù)測防治等方向研究，本文通信作者.