陳念輝，鄭 碩,2，徐德芳

（1.中國(guó)電建集團(tuán)華東勘測(cè)設(shè)計(jì)研究院有限公司，浙江杭州，311122；2.國(guó)家水電站大壩安全和應(yīng)急工程技術(shù)中心，浙江杭州，311122）

0 引言

爆破被廣泛應(yīng)用于各種巖土工程領(lǐng)域。爆破時(shí)，爆炸能將對(duì)周?chē)慕ㄖa(chǎn)生不利的影響。如何避免這種不利影響是爆破施工的一個(gè)重要問(wèn)題[1-6]。池潭水電站擴(kuò)建工程引水發(fā)電系統(tǒng)緊鄰池潭水電站大壩、發(fā)電廠房及開(kāi)關(guān)站，爆破施工可能對(duì)運(yùn)行電站產(chǎn)生不利影響。采用三維動(dòng)力有限元計(jì)算，研究各開(kāi)挖部位典型爆破振動(dòng)荷載作用下原大壩及大壩進(jìn)水口建筑物、廠房、中控室、開(kāi)關(guān)站、尾水平臺(tái)等保護(hù)對(duì)象的動(dòng)力響應(yīng)；對(duì)擴(kuò)建工程進(jìn)水口基坑道路、廠房及尾水出口開(kāi)挖、引水隧洞等關(guān)鍵爆破部位的爆破振動(dòng)影響做出評(píng)價(jià)，得到關(guān)鍵爆破部位的爆破振動(dòng)衰減特性及關(guān)鍵控制保護(hù)對(duì)象，為全面分析擴(kuò)建工程爆破施工對(duì)既有建筑物的安全影響評(píng)價(jià)提供基礎(chǔ)。

1 工程概況

池潭水電站位于福建省西北部泰寧縣境內(nèi)的金溪上，屬梯級(jí)電站的龍頭水電站，以發(fā)電為主，兼顧防洪等綜合利用。樞紐建筑物主要包括混凝土寬縫重力壩、壩后溢流廠房、壩內(nèi)輸水系統(tǒng)及泄水底孔、開(kāi)關(guān)站等，為二等大（2）型工程。水庫(kù)正常蓄水位高程275.00 m，水庫(kù)調(diào)節(jié)庫(kù)容6.26億m3。池潭水電站擴(kuò)建的水電站位于池潭水庫(kù)的左岸，利用池潭大壩擋水發(fā)電。新建樞紐建筑物主要包括輸水系統(tǒng)、地面廠房及升壓開(kāi)關(guān)站等。廠房位于左岸壩軸線下游約300 m處，引水發(fā)電系統(tǒng)由岸坡豎井式進(jìn)水口、引水隧洞、壓力鋼管等組成，廠區(qū)建筑物主要由主廠房、副廠房、220 kV開(kāi)關(guān)站、進(jìn)廠公路、廠前區(qū)及尾水渠等組成。主副廠房位于左岸且沿江布置。圖1為池潭水電站擴(kuò)建工程總布置。

圖1 池潭水電站擴(kuò)建工程總布置Fig.1 Layout of the expansion of Chitan hydropower station

2 鄰近建筑物的爆破振動(dòng)響應(yīng)三維動(dòng)力有限元數(shù)值模擬

新建的水電站引水發(fā)電系統(tǒng)緊鄰池潭水電站大壩、發(fā)電廠房及開(kāi)關(guān)站，基本將原大壩及廠房左岸側(cè)三面包圍。爆源點(diǎn)距原有建筑物60～200 m不等，擴(kuò)建工程施工將存在大范圍、多次反復(fù)爆破影響問(wèn)題。采用三維動(dòng)力有限元計(jì)算，研究各開(kāi)挖部位典型爆破振動(dòng)荷載作用下原大壩及大壩進(jìn)水口建筑物、廠房、中控室、開(kāi)關(guān)站、尾水平臺(tái)等保護(hù)對(duì)象的動(dòng)力響應(yīng)。

2.1 基本控制方程

確定結(jié)構(gòu)的動(dòng)力響應(yīng)最終歸結(jié)為求解結(jié)構(gòu)的三大基本方程問(wèn)題。結(jié)構(gòu)系統(tǒng)的動(dòng)力方程為：

式中，σij為柯西應(yīng)力；fi為體力密度；x¨為加速度。

其邊界條件如下：

（1）受力邊界（?b1）條件：

（2）位移邊界（?b2）條件：

（3）接觸邊界（?b3）條件：

結(jié)構(gòu)的質(zhì)量守恒方程為：

式中，ρ0為初始密度；v為相對(duì)體積，

結(jié)構(gòu)的能量守恒方程為：

其中，ε˙ij為應(yīng)變率張量。

將動(dòng)量方程和邊界條件寫(xiě)成以下Galerkin弱形式的平衡方程：

其中，δxi滿足位移邊界條件。

應(yīng)用散度定理和分部積分：

則Galerkin弱形式平衡方程（8）可改寫(xiě)成虛功原理的積分形式：

相對(duì)于靜力學(xué)問(wèn)題而言，彈性動(dòng)力學(xué)問(wèn)題由于引入了慣性力和阻尼力，其有限元解法要復(fù)雜許多，彈性動(dòng)力有限元方法的任務(wù)就是對(duì)上述基本方程進(jìn)行數(shù)值求解。其單元?jiǎng)恿ζ胶夥匠淌綖椋?/p>

式中，[Me]、[Ce]、[Ke]、[Fe]分別是單元的質(zhì)量矩陣、阻尼矩陣、剛度矩陣和結(jié)點(diǎn)載荷向量。分別對(duì)四個(gè)單元的矩陣進(jìn)行集成，形成整個(gè)系統(tǒng)的質(zhì)量矩陣[M]、阻尼矩陣[C]、剛度矩陣[K]和結(jié)點(diǎn)荷載向量[F]，相應(yīng)得到整個(gè)系統(tǒng)的動(dòng)力平衡方程為：

2.2 計(jì)算模型及方法

如圖2所示建立池潭水電站擴(kuò)建工程三維模型，模型包括引水隧洞、擴(kuò)機(jī)廠房、擴(kuò)機(jī)電站進(jìn)水口、池潭大壩和開(kāi)關(guān)站等建筑物。計(jì)算模型采用六面體單元?jiǎng)澐?，總單元?shù)與節(jié)點(diǎn)數(shù)分別為199 449和203 119。

圖2 桐子林水庫(kù)蓄水初期不同壩前水位時(shí)二灘尾水水位-流量關(guān)系Fig.2 Fitted tail water levels and flow rates of Ertan hydropower station with different water levels upstream of Tongzilin hydropower station

圖2 池潭水電站擴(kuò)建工程計(jì)算模型Fig.2 Three dimensional numerical model of expansion of Chitan hydropower station

模型四周及底面設(shè)置無(wú)反射邊界。計(jì)算中時(shí)間積分采用中心差分格式，采用分配參數(shù)法求解炸藥與巖體的接觸碰撞問(wèn)題。由于本項(xiàng)目中關(guān)注的保護(hù)對(duì)象均距離爆源很遠(yuǎn)，爆破振動(dòng)對(duì)保護(hù)對(duì)象的動(dòng)力響應(yīng)基本屬于彈性振動(dòng)的范圍，因此考慮計(jì)算的可行性和實(shí)際的合理性，選用彈性材料來(lái)模擬巖石介質(zhì)可以滿足本項(xiàng)目的計(jì)算精度和要求。

2.3 計(jì)算工況及爆破振動(dòng)考察點(diǎn)的選擇

根據(jù)引水隧洞和擴(kuò)建工程廠房的平面總體布置方案，考慮到國(guó)內(nèi)巖體開(kāi)挖爆破技術(shù)的發(fā)展現(xiàn)狀，并參考類似工程的實(shí)踐經(jīng)驗(yàn)，在動(dòng)力有限元模擬中的總體開(kāi)挖方法和方案確定上做了以下假定：進(jìn)水口部位的邊坡及基坑爆破開(kāi)挖和預(yù)留巖坎的爆破拆除主要采用深孔臺(tái)階爆破方式，邊坡及基坑爆破開(kāi)挖孔深約6～8 m，巖坎拆除孔深最大30～40 m。爆破采用排孔起爆，臺(tái)階高度8 m，單孔藥量20 kg，孔間距2.5 m，排距2 m，最大單響藥量為80～120 kg。新建廠房基礎(chǔ)的開(kāi)挖臺(tái)階高度6～8 m，單響藥量為40～80 kg，其他參數(shù)同進(jìn)水口部位。對(duì)于引水隧洞的開(kāi)挖，根據(jù)已有工程經(jīng)驗(yàn)，掏槽孔的爆破振動(dòng)效應(yīng)最為明顯，因此考慮引水隧洞掏槽爆破開(kāi)挖對(duì)工程區(qū)域保護(hù)對(duì)象爆破振動(dòng)的最不利影響。掏槽爆破的孔深3.5 m，4孔一響，最大單響藥量為15～20 kg。

根據(jù)工程總體開(kāi)挖方案和爆破方法，結(jié)合工程區(qū)域內(nèi)不同保護(hù)對(duì)象與爆源的相對(duì)位置，并考慮到不同類型爆源的距離和實(shí)際影響，在三維動(dòng)力有限元模擬中分析了以下3個(gè)部位的爆破施工影響：（1）進(jìn)水口部位的深孔臺(tái)階爆破；（2）廠房部位的深孔臺(tái)階爆破；（3）引水隧洞部位（距大壩左岸壩肩最近處）的淺孔爆破。

各種計(jì)算工況的具體爆破參數(shù)見(jiàn)表1，基巖物理力學(xué)參數(shù)見(jiàn)表2。

表1 各計(jì)算工況爆破參數(shù)Table 1 Blasting parameters under different calculation conditions

表2 巖石主要物理力學(xué)參數(shù)建議值Table 2 Suggested values of rock physical and mechanical parameters

為了準(zhǔn)確細(xì)致地對(duì)大壩等需要保護(hù)建筑物的振動(dòng)響應(yīng)規(guī)律進(jìn)行分析，在壩肩、中控室、開(kāi)關(guān)站等關(guān)鍵部位設(shè)置爆破振動(dòng)影響評(píng)價(jià)考察點(diǎn)，通過(guò)分析不同位置部位的振動(dòng)速度響應(yīng)，分析擴(kuò)建施工對(duì)需要保護(hù)建筑物的影響，考察點(diǎn)的布置如圖3所示。

圖3 監(jiān)測(cè)點(diǎn)布置圖Fig.3 Arrangement of monitoring points

2.4 爆破振動(dòng)荷載及施加方式

計(jì)算時(shí)荷載分別施加在進(jìn)水口、引水隧洞、新建廠房中央部位。為簡(jiǎn)化數(shù)值模擬，可以將粉碎區(qū)和破碎區(qū)視為爆破振動(dòng)源的一部分，將簡(jiǎn)化的爆破荷載曲線施加在破碎區(qū)的外邊界上，即彈性邊界上。爆破荷載的作用形式簡(jiǎn)化為三角形，如圖4（a）所示，其中荷載上升時(shí)間為1 ms，荷載持續(xù)時(shí)間為7 ms。不同部位的爆破荷載垂直施加于爆區(qū)等效彈性邊界巖體表面，如圖4（b）、（c）、（d）所示。

圖4 簡(jiǎn)化的爆破振動(dòng)荷載及爆破荷載的施加方式Fig.4 Simplified blasting vibration load and its loading mode

3 計(jì)算結(jié)果分析

以進(jìn)水口深孔臺(tái)階爆破引起的大壩等的振動(dòng)響應(yīng)為例，表3給出了進(jìn)水口典型爆破振動(dòng)荷載作用下的三維動(dòng)力有限元計(jì)算結(jié)果。在單響藥量60 kg的爆破振動(dòng)荷載作用下，保護(hù)對(duì)象中控室、開(kāi)關(guān)站的質(zhì)點(diǎn)峰值振動(dòng)速度分別為0.30 cm/s和0.39 cm/s，對(duì)應(yīng)的附加最大拉應(yīng)力分別為14.33 kPa和47.80 kPa；防滲帷幕處的質(zhì)點(diǎn)峰值振動(dòng)速度為0.52 cm/s，附加最大拉應(yīng)力為64.43 kPa；進(jìn)水口閘門(mén)和左岸壩肩處的質(zhì)點(diǎn)峰值振動(dòng)速度分別為0.31 cm/s和1.39 cm/s，對(duì)應(yīng)的附加最大拉應(yīng)力分別為4.95 kPa和52.57 kPa。除進(jìn)水口閘門(mén)外，其他關(guān)鍵位置質(zhì)點(diǎn)的峰值振動(dòng)速度都出現(xiàn)在水平徑向。此外，雖然壩肩和防滲帷幕的爆心距均為145 m，但由于所處的方位和結(jié)構(gòu)本身性質(zhì)不同，導(dǎo)致質(zhì)點(diǎn)峰值速度相差一倍還多。計(jì)算沒(méi)有考慮結(jié)構(gòu)初始應(yīng)力作用，因此上述應(yīng)力為爆破振動(dòng)引起的附加動(dòng)應(yīng)力。可以看出爆破振動(dòng)產(chǎn)生附加動(dòng)應(yīng)力較小，基本不會(huì)改變大壩及附屬結(jié)構(gòu)的總應(yīng)力狀態(tài)。爆破振動(dòng)速度沿大壩和廠房結(jié)構(gòu)的高度方向存在一定的放大效應(yīng)：?jiǎn)雾懰幜繛?0 kg時(shí)，壩基面的質(zhì)點(diǎn)振動(dòng)速度為0.41 cm/s，而同一壩段壩頂?shù)馁|(zhì)點(diǎn)振動(dòng)速度為1.54cm/s。

表3 進(jìn)水口典型爆破振動(dòng)荷載作用下各考察點(diǎn)的振速峰值（單響藥量60 kg）Table 3 Peak vibration velocity of monitoring points with water inlet under typical blasting vibration load

表4為進(jìn)水口部位深孔臺(tái)階爆破、新建廠房深孔臺(tái)階爆破和引水隧洞全斷面開(kāi)挖爆破時(shí)各保護(hù)對(duì)象峰值響應(yīng)速度。進(jìn)水口部位深孔臺(tái)階爆破施工的主要保護(hù)對(duì)象為防滲帷幕、左岸壩肩、開(kāi)關(guān)站、中控室，其質(zhì)點(diǎn)峰值振速（單響藥量60 kg）分別為0.52 cm/s、1.39 cm/s、0.39 cm/s和0.30 cm/s；新建廠房部位深孔臺(tái)階爆破施工的主要保護(hù)對(duì)象為開(kāi)關(guān)站、發(fā)電機(jī)層、中控室、尾水平臺(tái)，其質(zhì)點(diǎn)峰值振速（單響藥量80 kg）分別為0.41 cm/s、0.49 cm/s、0.28 cm/s和0.38 cm/s；引水隧洞淺孔掏槽爆破施工的主要保護(hù)對(duì)象為壩肩防滲帷幕和開(kāi)關(guān)站，其質(zhì)點(diǎn)峰值振速（單響藥量10.8 kg）分別為0.23 cm/s和0.06 cm/s。各開(kāi)挖部位的爆破振動(dòng)三維動(dòng)力有限元結(jié)果均表明，振動(dòng)響應(yīng)沿高程有一定的放大效應(yīng)，如進(jìn)水口深孔臺(tái)階爆破時(shí)，壩基面和壩頂?shù)恼駝?dòng)峰值分別為0.41 cm/s和1.54 cm/s；新建廠房深孔臺(tái)階爆破時(shí)，壩基面和壩頂?shù)恼駝?dòng)峰值分別為0.17 cm/s和0.50 cm/s；引水隧洞鉆孔爆破時(shí)，壩基面和壩頂?shù)恼駝?dòng)峰值分別為0.07 cm/s和0.20 cm/s。

表4 各保護(hù)對(duì)象振動(dòng)峰值響應(yīng)Table 4 The maximum vibration response of each protected structure

為了分析不同單響藥量下各主要保護(hù)對(duì)象（包括防滲帷幕、大壩、開(kāi)關(guān)站和廠房）振動(dòng)速度的變化規(guī)律，根據(jù)爆破振動(dòng)傳播規(guī)律，分別得到了進(jìn)水口深孔臺(tái)階爆破、引水隧洞爆破和廠房深孔臺(tái)階爆破時(shí)各主要保護(hù)對(duì)象在不同單響藥量和爆心距條件下的峰值振動(dòng)速度，如圖5所示。

由圖5（a）可知，當(dāng)擴(kuò)建工程進(jìn)水口在距開(kāi)關(guān)站最近處（124 m）爆破施工時(shí)，開(kāi)關(guān)站在單響藥量30 kg條件下的振動(dòng)峰值速度為0.51 cm/s，可見(jiàn)開(kāi)關(guān)站是進(jìn)水口爆破施工的控制保護(hù)對(duì)象。當(dāng)擴(kuò)建工程進(jìn)水口在距防滲帷幕最近處（40 m）爆破施工時(shí)，開(kāi)關(guān)站在單響藥量30 kg條件下的振動(dòng)峰值速度為1.30 cm/s，可見(jiàn)防滲帷幕是進(jìn)水口爆破施工的控制保護(hù)對(duì)象。因此，在距開(kāi)關(guān)站和防滲帷幕最近處的進(jìn)水口爆破施工過(guò)程中，需要針對(duì)爆破振動(dòng)對(duì)開(kāi)關(guān)站的影響開(kāi)展跟蹤監(jiān)測(cè)，依據(jù)實(shí)測(cè)資料最終確定合理的振動(dòng)控制標(biāo)準(zhǔn)；由圖5（c）可知，當(dāng)單響藥量達(dá)到120 kg時(shí)，各保護(hù)對(duì)象的振動(dòng)速度均低于其振動(dòng)安全控制標(biāo)準(zhǔn)；當(dāng)藥量增加到130 kg時(shí)，開(kāi)關(guān)站的振動(dòng)速度超過(guò)了控制標(biāo)準(zhǔn)，所以對(duì)于廠房基礎(chǔ)開(kāi)挖，單響藥量須嚴(yán)格控制在120 kg以內(nèi)。

圖5 保護(hù)對(duì)象在不同爆心距和單響藥量下的峰值振動(dòng)速度Fig.5 Peak vibration velocity of protected structures from different distances and explosive charges

圖5表明，各保護(hù)對(duì)象的最大振動(dòng)速度隨單響藥量的增加而增大；相同藥量時(shí)，最大振動(dòng)速度隨爆心距增加而減小。通過(guò)控制單響藥量，可以防止爆破振動(dòng)對(duì)保護(hù)對(duì)象產(chǎn)生可見(jiàn)的有害影響。

4 結(jié)語(yǔ)

根據(jù)池潭水電站擴(kuò)建工程地質(zhì)條件及設(shè)計(jì)方案，基于進(jìn)水口、廠房基礎(chǔ)的深孔臺(tái)階爆破和引水隧洞全斷面開(kāi)挖爆破對(duì)現(xiàn)有鄰近大壩和發(fā)電設(shè)施設(shè)備的影響分析，提出以下結(jié)論和建議：

（1）引水隧洞全斷面淺孔爆破引起的振動(dòng)對(duì)保護(hù)對(duì)象的影響相對(duì)較小，但在距離大壩最近處爆破施工時(shí)，應(yīng)重點(diǎn)關(guān)注壩肩防滲帷幕的振動(dòng)響應(yīng)；進(jìn)水口開(kāi)挖部位由于距原大壩、廠房、開(kāi)關(guān)站等建（構(gòu)）筑物較近，對(duì)保護(hù)對(duì)象的影響較大，為關(guān)鍵控制爆源；新建廠房主體部位開(kāi)挖時(shí)對(duì)原廠房建筑物及設(shè)備的影響較大，且新建副廠房和安裝間距原大壩更近，爆破施工時(shí)需要嚴(yán)格控制該部位引起的振動(dòng)響應(yīng)。

（2）爆破振動(dòng)影響三維動(dòng)力有限元數(shù)值模擬分析發(fā)現(xiàn)，保護(hù)對(duì)象的最大振動(dòng)速度隨單響藥量的增加而增大；相同藥量時(shí)，最大振動(dòng)速度隨爆心距增加而減小。針對(duì)池潭水電站擴(kuò)建工程，通過(guò)控制最大單響藥量、優(yōu)化爆破參數(shù)及優(yōu)化起爆網(wǎng)絡(luò)等措施，能夠?qū)⒈普駝?dòng)對(duì)鄰近建（構(gòu)）筑物的影響控制在允許范圍。對(duì)進(jìn)水口爆源，在距離壩后開(kāi)關(guān)站130.0～150.0 m以內(nèi)的部位爆破時(shí)，最大單響藥量應(yīng)控制在20.0 kg以內(nèi)，其他區(qū)域控制在40.0～60.0 kg；對(duì)引水隧洞的爆破開(kāi)挖，距離防滲帷幕40.0～60.0 m范圍內(nèi)的部位，最大單響藥量應(yīng)控制在15.0～25.0 kg，其他區(qū)域控制在20.0～40.0 kg；對(duì)廠房基礎(chǔ)部位的爆破開(kāi)挖，最大單響藥量應(yīng)控制在70.0～120.0 kg。