（成都理工大學(xué)地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護國家重點實驗室，四川成都610059）

（成都理工大學(xué)地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護國家重點實驗室，四川成都610059）

為了揭示風(fēng)載作用下震裂山體的崩塌機制，建立其穩(wěn)定性評價系統(tǒng)，在巖體震裂特征現(xiàn)場調(diào)查和河谷風(fēng)現(xiàn)場測量的基礎(chǔ)上，以茂縣石大關(guān)鄉(xiāng)崩塌為例，開展了不同條件下的大型風(fēng)洞試驗研究.研究結(jié)果表明：石大關(guān)鄉(xiāng)崩塌為震裂山體在風(fēng)載助推作用下發(fā)生潰崩的典型實例，此類災(zāi)害因震裂山體獨特的結(jié)構(gòu)特征，使得結(jié)構(gòu)面對失穩(wěn)模式的控制作用不再明顯，加之河谷風(fēng)受地形影響，分布規(guī)律復(fù)雜多變，具備與通常重力環(huán)境下斜坡失穩(wěn)機制迥異的特征，整體以潰崩式失穩(wěn)，表現(xiàn)為迅速解體—潰散—垮塌；巖體震裂損傷是崩塌產(chǎn)生的基礎(chǔ)，漸近性風(fēng)化是主要致災(zāi)因素，河谷強風(fēng)的助推是崩塌發(fā)生的誘因之一；巖塊的失穩(wěn)風(fēng)速與風(fēng)的入射夾角、巖塊的高寬比和后緣縫隙的寬度成反比.

地震；河谷風(fēng)；石大關(guān)崩塌；風(fēng)洞試驗；穩(wěn)定性評價

汶川地震誘發(fā)了大量崩塌、滑坡次生地質(zhì)災(zāi)害［1-2］，在此同時，由于強震對一定范圍內(nèi)地質(zhì)環(huán)境的沖擊，還發(fā)育了一類分布范圍廣、規(guī)模巨大、潛在威脅較大的災(zāi)害——震裂山體，其具有大變形、大位移等顯著的宏觀特征［3］.經(jīng)震后震區(qū)地質(zhì)災(zāi)害的排查工作，發(fā)現(xiàn)上千處“裂”而未“滑”、“松”而未“動”的震裂松動山體.震裂縫的大量發(fā)育使得斜坡巖體結(jié)構(gòu)受到了嚴(yán)重破壞，并為雨水下滲提供了良好的通道，坡體穩(wěn)定性在外部環(huán)境的影響下逐年劣化.因此，在強震過后，震區(qū)的滑坡、崩塌、泥石流等次生地質(zhì)災(zāi)害呈現(xiàn)出增強的趨勢［4-5］，黃潤秋將此特征稱為強震地質(zhì)災(zāi)害的“后效應(yīng)”，且震后地質(zhì)災(zāi)害往往具有極強的隱蔽性和突發(fā)性，難以識別和預(yù)警［6］.2009年徹底關(guān)大橋在連日大雨后發(fā)生大面積崩塌［7］、2010年在降雨作用下發(fā)生在雁門鄉(xiāng)索橋村的滑坡及同年舟曲發(fā)生的特大泥石流災(zāi)害都是坡體內(nèi)部形成“內(nèi)傷”的震裂巖體在外因的誘導(dǎo)下觸發(fā)的震后地質(zhì)災(zāi)害［8-9］.

汶川地震震區(qū)多高山峽谷，谷風(fēng)盛行，震后殘留的大量松散碎屑物在強風(fēng)作用下飛沙走石，滾落的碎塊石嚴(yán)重威脅著災(zāi)區(qū)人民的生命財產(chǎn)安全. 2014年7月17日下午2時45分，發(fā)生在四川省阿壩州G213線K774+600 m處茂縣境內(nèi)（石大關(guān)鄉(xiāng)超限站附近）的崩塌就是在河谷風(fēng)的助推作用下觸發(fā)的典型震后地質(zhì)災(zāi)害.塌方量高達4×103m3，飛石砸中13輛車，造成10人死亡，22人受傷.此次血淋淋的教訓(xùn)再次為震裂山體的防治工作敲響了警鐘.

針對震裂山體的研究，雖前人在震裂變形的分布規(guī)律、形成機理等方面取得了豐碩的成果［10］，但風(fēng)載作為震裂山體失穩(wěn)的觸發(fā)因素之一，并未有相關(guān)人員對高山河谷地區(qū)河谷風(fēng)的分布規(guī)律及其對震裂山體穩(wěn)定性的影響進行系統(tǒng)的研究.故本文在汶川地震巖體震裂特征大量現(xiàn)場調(diào)查和河谷風(fēng)現(xiàn)場測量工作的基礎(chǔ)上，以石大關(guān)鄉(xiāng)崩塌為例，論述了震裂山體在河谷風(fēng)助推作用下發(fā)生崩塌的特征及成因機制；結(jié)合風(fēng)洞試驗得出風(fēng)載作用下巖塊的啟動機理及失穩(wěn)風(fēng)速計算公式，建立了風(fēng)載作用下震裂山體穩(wěn)定性評價系統(tǒng)，對震裂山體防治工作具有指導(dǎo)性的意義.

1 汶川地震巖體震裂特征

由于地區(qū)地震動強度的差異性，在汶川地震極震區(qū)，地震動極其強烈，斜坡多失穩(wěn)破壞，崩滑體后壁陡立粗糙，后壁殘留的震裂松弛巖體及坡腳處于臨界狀態(tài)堆積體的穩(wěn)定性問題為震后震區(qū)面臨的主要問題.在汶川地震的強震區(qū)，地震強度較高，斜坡多以變形為主，特別是在斜坡陡緩轉(zhuǎn)折部位，由于其對地震波有明顯的放大效應(yīng)，震裂現(xiàn)象尤為突出，具有大變形、大位移等顯著的宏觀特征，規(guī)模普遍較大.

以震裂巖體結(jié)構(gòu)精測為基礎(chǔ)，得出巖體震裂最典型的宏觀特征為結(jié)構(gòu)面的松動和張開，其發(fā)育模式分為沿結(jié)構(gòu)面的延伸方向張開和沿結(jié)構(gòu)面延伸方向斜交發(fā)育.且由于斜坡越高陡，地震波放大效應(yīng)越明顯，在山脊或臨空條件較好的部位往往震裂縫的張開度較大.故在汶川地震中巖體震裂縫的張開度隨著坡高的降低和埋深的增大逐漸減?。▓D1），就同一條震裂縫而言總體上呈上寬下窄的特征.

圖1 張開度隨坡高區(qū)間分布圖Fig.1 Number of different joint openings with the slope height

震裂張開的層面、節(jié)理面絕大部分無充填或者少量充填，充填物質(zhì)以碎石、巖屑為主.由于震裂巖體的強度是在極短的時間內(nèi)銳減的，而其他風(fēng)化作用下形成的松動巖體的強度多是日積月累慢慢衰減的，故震裂巖體具有發(fā)育裂縫新鮮，裂隙面上有時可見延伸較好的擦痕，且多無充填，充填物的膠結(jié)程度都較低，形成時間短、規(guī)模大等顯著異于其他風(fēng)化、卸荷作用形成的松動巖體的特征.

2 石大關(guān)震裂山體崩塌及其河谷風(fēng)特征

2.1 石大關(guān)崩塌概況

石大關(guān)崩塌地處青藏高原東緣，為典型的深切河谷地貌.從地形上看，斜坡中部發(fā)育陡坎，坡度集中在45°～60°之間，崩塌源處近直立.坡體上中下部分別發(fā)育大光壁，猶如3面隨風(fēng)飄揚的旗幟，故被當(dāng)?shù)卮迕穹Q為“三桿旗”.而本次的崩塌源為坡體上部的第一桿旗，為4面臨空的小山包，地勢高陡突兀，高程大約2 340 m，距國道213的高差720 m左右（圖2）.

圖2 石大關(guān)鄉(xiāng)崩塌剖面圖Fig.2 Cross-section map of collapse at Shidaguan

石大關(guān)鄉(xiāng)斜坡由三疊系中統(tǒng)雜谷腦組（T2Z）灰色中厚層石英絹云母千枚巖夾變質(zhì)石英砂巖組成，發(fā)育3組優(yōu)勢機構(gòu)面：①層面，產(chǎn)狀N60°～70°E/SE/75°～85°，薄～中厚層結(jié)構(gòu)；②結(jié)構(gòu)面J1（底面），產(chǎn)狀N30°～40°E/NW/10°～20°；③結(jié)構(gòu)面J2（切割面），產(chǎn)狀N60°～70°W/SW/60°～70°.坡面坡向270°，坡度45°～60°.

由赤平投影圖（圖3）可以，看出層面和切割結(jié)構(gòu)面的傾向與坡向大角度斜交，底部結(jié)構(gòu)面與坡向相同，且傾角小于坡角，此外坡面與層面、坡面與J2、層面與J1、J1與J2的組合交棱的傾向都與坡向相同，且傾角都小于坡角.切割出的塊體不穩(wěn)定.

研究區(qū)受龍門山斷裂帶影響，地質(zhì)構(gòu)造較復(fù)雜，受多次地震破壞的洗禮，且崩塌距離石大關(guān)斷裂僅1 km［11］.加之崩塌源高陡突兀，臨空條件好，地震放大效應(yīng)明顯［12-13］，強烈地震力作用下，巖體結(jié)構(gòu)遭受了嚴(yán)重破壞，使原本強度不高的巖體更加脆弱.通過現(xiàn)場調(diào)查，發(fā)現(xiàn)張開的裂縫多沿陡傾層面和J2結(jié)構(gòu)面發(fā)育，張開程度普遍較大，呈上寬下窄的形態(tài)，縫面新，多平直，充填物較少，局部充填碎塊石，且無膠結(jié)，巖體結(jié)構(gòu)特征與震裂巖體特征吻合，地震作用下形成的震裂結(jié)構(gòu)為崩塌的發(fā)生提供了物質(zhì)基礎(chǔ).

研究區(qū)地屬亞熱帶季風(fēng)氣候，冬冷夏涼，晝夜溫差大，夏季高溫多雨，年降水量500 mm左右.震裂結(jié)構(gòu)為雨水的下滲提供了良好通道，在反復(fù)凍脹及其他風(fēng)化作用下，震裂縫緩慢擴大，進一步劣化了坡體的穩(wěn)定性，據(jù)當(dāng)?shù)鼐用穹从?，在崩塌發(fā)生前崩塌源最大的張開裂縫約達15 mm.

圖3 結(jié)構(gòu)面赤平投影Fig.3 Stereographic projection of structural plane

同時，經(jīng)走訪調(diào)查得知，當(dāng)?shù)?月初至7月中上旬一直大雨連綿.持續(xù)的降雨造成巖體飽和、自重增加、巖體結(jié)構(gòu)面抗剪強度降低，大大增加了山體崩塌發(fā)生的幾率［14］.根據(jù)茂縣氣象局提供的降雨資料（圖4），7月10日至7月12日累計降雨量達到了60 mm，在崩塌發(fā)生前4天內(nèi)都無明顯降雨現(xiàn)象，說明崩塌發(fā)生時并沒有動水壓力或者靜水壓力的作用.但崩塌發(fā)生當(dāng)日氣溫較高，烈日當(dāng)頭，巖體受熱膨脹，也逐漸弱化了損傷巖體的完整性和穩(wěn)定性.

圖4 7月6日～15日石大關(guān)鄉(xiāng)降雨量曲線Fig.4 Curve of rainfall at Shidaguan from July 6 to July 15

2.2 河谷風(fēng)特征

由于氣流在岷江河谷受到高山峻嶺的阻擋和挾持導(dǎo)向作用而谷風(fēng)盛行.本文選取了茂縣飛虹鄉(xiāng)（N31°9′33.60″，E103°3′17.66″）及汶川地區(qū)多個測量剖面，從谷底1 500～2 500 m高程按需布置多個測量點，運用JL-SXJ風(fēng)速記錄儀對岷江河谷風(fēng)速分布規(guī)律進行研究，得出了以下結(jié)論：

①定時起風(fēng)，午前多靜風(fēng)，平均風(fēng)速＜4 m/s；午后多大風(fēng)，平均風(fēng)速＞4 m/s.一天中最大風(fēng)速多出現(xiàn)在午后16時，達8～10 m/s.一年中2～5月及8月風(fēng)速相對較大，7、10、11月相對較小，其中7月風(fēng)速最小，為3.7 m/s.

②岷江河谷內(nèi)，海拔越高，河谷越開闊，河谷對空氣的壓縮作用越弱，流速因此降低，故風(fēng)速自谷底而上并非呈指數(shù)或線性型增長趨勢，其最大值多出現(xiàn)于山腰位置處.且氣流在沿山脈爬升過程中受復(fù)雜地形的影響，在地形突變位置處可形成較強的局部風(fēng).特別是在下陡上緩山體的轉(zhuǎn)折部位，會出現(xiàn)谷風(fēng)的極大值并伴有渦旋氣流現(xiàn)象.

③當(dāng)氣流方向與河谷走向之間夾角＜60°時，狹管效應(yīng)顯著，即氣流在河谷段有明顯的加速現(xiàn)象，且狹管效應(yīng)的強弱與河谷的寬度呈反比.在底部河谷寬度最小，狹管效應(yīng)最明顯，其次為地形的轉(zhuǎn)折部位，而海拔越高，河谷越開闊，狹管效應(yīng)逐漸減弱.

而在石大關(guān)鄉(xiāng)崩塌發(fā)生當(dāng)日下午14時，茂縣觀測站（石大關(guān)鄉(xiāng)政府處，高程1 742 m）監(jiān)測到了高達8.8 m/s的五級勁風(fēng)，這對于夏季（6～9月）平均風(fēng)速才3.4m/s的地區(qū)而言比較異常［15］. 2015年7、8月份作者們先后兩次帶上專業(yè)的風(fēng)速、風(fēng)向測量儀到石大關(guān)鄉(xiāng)收集河谷風(fēng)資料.調(diào)查得出崩塌源處在午后13時后多大風(fēng)，一直到下午17時河谷風(fēng)多集中在6～9 m/s，最大風(fēng)速可達14 m/s，崩塌源處風(fēng)速與茂縣觀測站處風(fēng)速比值集中在1.4～2.0之間.崩塌源處風(fēng)向集中在290°～310°，與主崩方向近似.

結(jié)合岷江河谷風(fēng)速分布規(guī)律可得：石大關(guān)鄉(xiāng)崩塌所在河谷為典型對稱“V”字形深切河谷地貌，河谷寬度較小，氣流方向與河谷走向小于60°，狹管效應(yīng)明顯，加之崩塌源高陡突兀，位于山腰偏上位置，局部地形為陡緩轉(zhuǎn)折端，屬風(fēng)速較大處，狹管效應(yīng)更為突出，最終在巖壁和臨空面兩側(cè)形成的壓力差作用下，向臨空方向發(fā)生崩塌（圖5、6）.

圖5 石大關(guān)鄉(xiāng)崩塌源崩塌后形態(tài)特征（單位：m）Fig.5 Residual rock mass of collapse source at Shidaguan（unit：m）

圖6 石大關(guān)鄉(xiāng)崩塌全景圖Fig.6 Panoramic photos of Shidaguan collapse

綜上所述，石大關(guān)鄉(xiāng)崩塌加深了對河谷風(fēng)對地質(zhì)災(zāi)害觸發(fā)作用的認識，證實了震裂山體在風(fēng)載助推作用下發(fā)生大型崩塌的可能性.此類崩塌的演化過程可歸納為：早期地震作用形成松動的震裂結(jié)構(gòu)，為雨水的下滲提供了良好通道，在進一步漸進性風(fēng)化作用下，裂縫緩慢擴展，逐漸形成結(jié)構(gòu)十分破碎的山體.且因崩塌源獨特震裂結(jié)構(gòu)特征，使得結(jié)構(gòu)面對其失穩(wěn)模式的控制作用不再明顯，加之河谷風(fēng)受地形影響，方向及大小復(fù)雜多變，在河谷風(fēng)的助推作用下，先發(fā)生表層巖體的小崩小塌，后部巖體因在常年的構(gòu)造和風(fēng)化作用下自身穩(wěn)定性本來就較差，當(dāng)一旦具備臨空條件，便迅速解體—潰散—垮塌，就單個巖塊而言雖還是以傾倒、滑移破壞為主，但整體以潰崩式失穩(wěn)，表現(xiàn)出一跨到底的特征，散落范圍極大，崩落持續(xù)時間長，危害性極大，極具突發(fā)性（圖6）.

3 風(fēng)洞試驗研究

考慮到震裂山體結(jié)構(gòu)特征和風(fēng)載作用下觸使孤石發(fā)生二次啟動的工程經(jīng)驗及試驗可操作性，為提出快速、準(zhǔn)確的穩(wěn)定性定量評價標(biāo)準(zhǔn)，在西南交通大學(xué)風(fēng)工程試驗研究中心［16］（圖7）采用完全脫離母巖的巖塊進行試驗研究，探討了風(fēng)載作用下震裂巖體的啟動機理.

3.1 試驗設(shè)計

震裂巖體因震裂縫的發(fā)育，使巖體的結(jié)構(gòu)受到了嚴(yán)重破壞，巖體被切割成不同形狀、與母巖之間的關(guān)系各不相同的巖塊.

從風(fēng)載的入射角度（初始氣流方向與試樣迎風(fēng)面豎軸線之間的夾角）、巖樣高寬比、巖樣與巖壁之間縫隙的寬度3個方面進行測壓試驗和實心模型啟動過程試驗，其中入射角度α取0°、30°、45°、60°，高寬比α=z/x取1、2、2.5（x恒定為30 cm），縫隙寬度取0、20、50 mm.

圖7 風(fēng)洞試驗場示意圖（單位：m）Fig.7 Schematic of wind tunnel test（unit：m）

在測壓試驗中使用電子掃描閥和測壓管來測量試樣表面壓力的分布情況，使用高速攝像機來記錄試樣失穩(wěn)前的氣動現(xiàn)象，從試樣各面不同壓力分布情況來研究其啟動機理，為后續(xù)的力學(xué)分析提供依據(jù).

風(fēng)的入射角度是通過移測架（圖8）中的底部轉(zhuǎn)盤調(diào)整的，下伏支撐板的傾角設(shè)置為θ=15°，并為了測壓孔安裝便利，試驗中采用空心木塊來模擬巖塊.

由于測壓試驗中空心材料的缺陷性，設(shè)計了實心模型啟動過程試驗來更為真實地模擬巖塊在風(fēng)載作用下的啟動過程及失穩(wěn)模式，確定不同類型試樣的失穩(wěn)風(fēng)速.實心試樣采用石英砂巖塊，其中x恒定為50 mm，為測壓試驗的1/6，狹縫寬度同樣縮小6倍，其他條件與測壓試驗相同.試驗后通過力學(xué)分析，提出極限平衡狀態(tài)方程，一方面用于驗證風(fēng)洞試驗結(jié)果的準(zhǔn)確性，另一方面為定量評價工作提供依據(jù).

圖8 移測架結(jié)構(gòu)圖Fig.8 Test frame structure

3.2 試驗結(jié)果分析

3.2.1 入射角度

③當(dāng)α=45°時，試樣正面和右側(cè)面都對稱分布正風(fēng)壓，最大正風(fēng)=72 Pa出現(xiàn)于右側(cè)面靠前部位.而其余面均為負風(fēng)壓，最小負風(fēng)壓=-70 Pa位于頂面.試樣的失穩(wěn)風(fēng)速進一步降低至10.5 m/s，見圖9（c）.

圖9 不同入射角度塊石臨界失穩(wěn)狀態(tài)表面壓力分布Fig.9 Surface pressure distribution in critical state with different incident angles

（2）啟動過程研究

由表1和記錄影像可知，氣流入射角度的變化主要影響試樣開始旋轉(zhuǎn)滑移的啟動風(fēng)速，角度越大，啟動風(fēng)速越低，由于試樣的不斷旋轉(zhuǎn)，不斷改變著氣流作用于試樣的夾角，當(dāng)試樣的一面與氣流垂直時，便會迅速旋轉(zhuǎn)滑移失穩(wěn)墜落.

（3）力學(xué)分析

試驗結(jié)果表明，正方形試樣在不同入射角度下都以滑移方式失穩(wěn)，角度不同，失穩(wěn)風(fēng)速有所差異，力學(xué)分析如下：

①當(dāng)α=0°時，氣流作用于試樣時，將會在正面（迎風(fēng)面）產(chǎn)生正壓力，而其余各面由于氣流的分離和再附和而產(chǎn)生負吸力.由于氣流自身復(fù)雜性，目前尚未用簡單的計算方程來確定吸力的大小，故在力學(xué)分析中暫忽略吸力的作用，僅考慮在重力、氣流在正面產(chǎn)生的風(fēng)壓和底部摩擦力三者的共同作用下試樣的極限平衡狀態(tài).

②當(dāng)α=60°時，氣流在正面產(chǎn)生的風(fēng)壓、重力下滑分力的合力和底部的最大靜摩擦力在一條直線上.

表1 不同入射角度實心模型啟動過程試驗結(jié)果Tab.1 Initiation process results of solid model with different incident angles m/s

而在低速氣流的條件下，忽略體力作用且無粘性時，任一流線上各點的標(biāo)準(zhǔn)伯努利方程為

式中：v為流線的風(fēng)速，m/s；

w1為單位面積的靜壓力，kN/m2；

V為空氣質(zhì)點的體積，m3；

m為空氣質(zhì)點的質(zhì)量，t，m=ρV；

C為常數(shù).令

式中：w為巖塊迎風(fēng)面與背風(fēng)面之間形成的壓力差；

取標(biāo)準(zhǔn)大氣壓下，重力加速度g為9.8 m/s2；

空氣容重γ1為0.012 018 kN/m3.

由此可得巖塊單位面積上所受的基本風(fēng)壓系數(shù)：

設(shè)正方體巖塊邊長為x，則巖塊正面受正壓力為

巖塊的最大靜摩擦力為

式中：μ0為靜摩擦因數(shù)，取0.5；

Ffom為最大靜摩擦力，單位KN；

γ2為木材容重，值取6.86 kN/m3；

θ為試樣下伏結(jié)構(gòu)面的傾角，值取15°；

d取0.5 cm.

巖塊重力的下滑分量為

根據(jù)極限狀態(tài)力學(xué)平衡條件可得：

由極限平衡方程計算得到α=0°、60°時，失穩(wěn)風(fēng)速分別為11.5、8.5 m/s，結(jié)果均比測壓試驗中小約1 m/s.這是因在極限平衡分析中忽略了氣流在通過試樣時，由于產(chǎn)生分離和再附和而在試樣上產(chǎn)生吸力，減小了作用于試樣底部下伏結(jié)構(gòu)面上的正壓力Fw，降低了最大靜摩擦力Ffom.

3.2.3 高寬比

（1）測壓試驗

由壓力分布圖（圖10）可知試樣正面（迎風(fēng)面）均為正風(fēng)壓，其余幾面均為負風(fēng)壓（吸力）.當(dāng)高寬比b=2.0時，正風(fēng)壓較大值出現(xiàn)于試樣正面中下部，最大正風(fēng)壓P+max=60 Pa.往兩側(cè)，正風(fēng)壓逐漸降低至15 Pa，試樣的失穩(wěn)風(fēng)速為9.25 m/s.當(dāng)b=2.5時，其各面風(fēng)壓的分布情況與b=2.0時相似，只是正風(fēng)壓較大值的位置有所上移，風(fēng)壓的作用點有所抬高，相同風(fēng)壓產(chǎn)生的力矩更大，試樣的失穩(wěn)風(fēng)速也就更小.故力矩為風(fēng)載作用下長方體巖塊失穩(wěn)的主導(dǎo)因素.

圖10 不同高寬比塊石臨界失穩(wěn)狀態(tài)表面壓力分布Fig.10 Surface pressure distribution in critical state with different height-width ratios of mass

（2）啟動過程研究

由表2和記錄影像可知：隨風(fēng)速的增大，長方體試樣都最終斜向下滑移一段距離并順風(fēng)劇烈晃動后順風(fēng)傾倒（側(cè)向傾倒）失穩(wěn).

b=2.5時較b=2.0時旋轉(zhuǎn)及滑移的過程相對較短，啟動過程如圖11.

表2 不同高寬比實心模型啟動過程試驗結(jié)果Tab.2 Initiation process results of solid model with different slenderness ratios m/s

圖11 高寬比（b=2.5）實心模型啟動過程Fig.11 Initiation process of solid model with height-width ratio of 2.5

（3）力學(xué)分析

此處同樣忽略背風(fēng)面負風(fēng)壓的作用，并由于此次試驗中試樣的尺寸較小，故將迎風(fēng)面的風(fēng)壓系數(shù)考慮為常數(shù)，長柱狀巖塊的受力情況如圖12所示.

風(fēng)壓作用下產(chǎn)生的傾倒力矩為

由于試樣在兩個高寬比條件下都以側(cè)向傾倒的形式失穩(wěn)，故Gsin θ產(chǎn)生的傾覆力矩等于0，僅考慮Gcos θ產(chǎn)生的抗傾覆力矩MG為

圖12 長柱狀巖塊側(cè)向受力分析示意圖Fig.12 Lateral force analysis of long-column block

極限平衡狀態(tài)下，傾倒力矩與抗傾覆力矩相等，即

由極限平衡方程計算得到當(dāng)高寬比b=2.0，2.5時，失穩(wěn)風(fēng)速分別為11.2、9.8 m/s，均比測壓試驗失穩(wěn)風(fēng)速小1.5 m/s左右.

這是由于實心模型啟動過程試驗中得出長方體試樣在發(fā)生傾倒前都會旋轉(zhuǎn)一定角度，而測壓試驗中當(dāng)試樣有輕微位移時便記錄為失穩(wěn)風(fēng)速，且力學(xué)計算中忽略了背面吸力產(chǎn)生的傾倒力矩Mx，即導(dǎo)致Mw偏大，故計算出的v0偏大.

3.2.2 后緣裂縫的張開程度

（1）測壓試驗

無論從理論計算分析還是野外調(diào)查發(fā)現(xiàn)，都證實了順向傾倒失穩(wěn)模式的存在.

測壓試驗中，由于由木板制作而成的試樣質(zhì)量偏輕，因左右兩側(cè)面的風(fēng)速差產(chǎn)生的傾倒力矩的效果遠不及順風(fēng)向風(fēng)壓產(chǎn)生的傾倒力矩的效果，故均以側(cè)向傾倒方式失穩(wěn).

從圖13的壓力分布情況可以總結(jié)出：

①由于氣流在通過右側(cè)面緊靠巖壁的狹小裂縫時，受兩壁摩擦作用的影響，風(fēng)速大大降低，而試樣左側(cè)臨空，兩側(cè)面風(fēng)壓不再對稱分布.靠巖壁側(cè)（右側(cè)）風(fēng)出現(xiàn)了正力，而臨空側(cè)全為負力.

②縫隙越小，緊靠巖壁一側(cè)（右側(cè)）產(chǎn)生的傾倒力矩越大.隨著縫隙的增大，該力矩對試樣發(fā)生傾倒失穩(wěn)的貢獻越低.

（2）啟動過程研究

由表3可以看出，當(dāng)裂縫寬度d=0 mm時，試樣在先輕微晃動后產(chǎn)生略微順時針旋轉(zhuǎn)，最終順向傾倒；而裂縫寬度d=3.5，8.0 mm時，試樣最終都以側(cè)向傾倒的形式失穩(wěn).

結(jié)果揭示了順向傾倒只發(fā)生在當(dāng)巖塊與巖壁之間縫隙足夠小時，啟動過程如圖14.

圖13 不同縫隙寬度塊石臨界失穩(wěn)狀態(tài)壓力分布Fig.13 Surface pressure distribution in critical state with different crack widths

（3）力學(xué)分析

首先設(shè)未受到擾動的巖塊臨空側(cè)風(fēng)速為v0，受到擾動時的靠縫隙側(cè)風(fēng)速減小為v1.由伯努利方程可得：

將v0、v1分別代入得：

設(shè)巖塊順風(fēng)方向長為x（m），寬為y（m），高為z（m），則兩側(cè)面風(fēng)壓差所產(chǎn)生的傾覆力矩為

由于巖塊以順向傾倒的形式失穩(wěn)，故順風(fēng)向風(fēng)壓產(chǎn)生的傾覆力矩為0（圖15）.

表3 不同縫隙寬度實心模型啟動試驗結(jié)果Tab.3 Initiation process results of solid model with different crack widths m/s

圖14 縫隙寬度（d=0）實心模型啟動過程Fig.14 Solid model with crack width（d=0）

設(shè)下伏結(jié)構(gòu)面傾角為θ，重力產(chǎn)生的傾覆力矩為

重力產(chǎn)生的抗傾覆力矩為

極限平衡狀態(tài)下，傾倒力矩與抗傾覆力矩相等：

圖15 長柱狀塊體順向傾倒力矩示意圖Fig.15 Overturn moment of long-column block

通過式（15）～（18）可以得出，發(fā)生順向傾倒失穩(wěn)的條件極為苛刻，需保證Gsin θ和Fw2足夠大，即巖塊重量、下伏結(jié)構(gòu)面的傾角、巖塊與巖壁之間縫隙的寬度都是決定其是否可能發(fā)生順向傾倒的關(guān)鍵.

故在相同質(zhì)量條件下，巖塊的失穩(wěn)風(fēng)速與入射角度、巖塊的高寬比、后緣縫隙的寬度都成反比.就正方形和長方形巖塊而言，在不同風(fēng)載入射角度、不同危巖高寬比、不同后緣裂縫張開程度的情況下主要以滑移、側(cè)向傾倒、順向傾倒失穩(wěn)為主（圖16），發(fā)生順向傾倒的條件極為苛刻.

圖16 風(fēng)載作用下正方體、長方體失穩(wěn)模式示意圖Fig.16 Failure modes of cube and cuboid under wind loading

4 風(fēng)載作用下震裂山體穩(wěn)定性評價方法

4.1 綜合指標(biāo)評價法

在文獻［17］中，袁進科等選取了7個因子作為崩塌識別指標(biāo)，采用專家取值法對指標(biāo)進行賦值，建立了震后崩塌的快速識別方法（CDI法）：

式中：

UCDI為崩塌識別的綜合指標(biāo)，值越大，發(fā)生崩塌概率越高；

Ri為第i個識別指標(biāo)因子分級指數(shù)；

Si為權(quán)重值；

λi為降雨修正系數(shù)［17］.

雖以震后震區(qū)100個崩塌災(zāi)害點的穩(wěn)定性評價為依據(jù)，驗證了此方法的合理性，但其并未考慮風(fēng)載對震裂山體穩(wěn)定性的影響.由風(fēng)洞試驗結(jié)果表明巖塊的形狀、后緣裂縫的張開程度決定了巖體的失穩(wěn)情況，而巖體震裂損傷程度是影響巖塊的形狀、后緣裂縫的張開程度的重要因素.

由河谷風(fēng)分布特征可知，山體地形突變處是河谷風(fēng)較大的部位，這些部位由于地震波的放大作用，同樣也是巖體震裂損傷程度相對較大的部位.因此，可用震裂損傷程度來表征風(fēng)載對震裂山體穩(wěn)定性的影響，將風(fēng)載修正系數(shù)βi引入CDI法中，得到風(fēng)載作用下震裂山體穩(wěn)定性評價指數(shù)（wind loading shattered rockmass discrimination index，WDI）.βi取1時，表征風(fēng)載對震裂山體的穩(wěn)定性不產(chǎn)生影響.

式中：

UWDI為風(fēng)載作用下震裂山體穩(wěn)定性評價指數(shù)，其值越大，發(fā)生崩塌的可能性越高；

βi為風(fēng)載修正系數(shù)，修正系數(shù)取值可按表4確定.

風(fēng)載作用下震裂山體穩(wěn)定性評價等級同樣分為3級，具體等級分級情況見表5.

表4 風(fēng)載條件下評價指標(biāo)修正系數(shù)Tab.4 Evaluation index correction under wind loading

根據(jù)上述分析，對石大關(guān)鄉(xiāng)崩塌按照評價指標(biāo)取值標(biāo)準(zhǔn)進行賦值，求得天然狀態(tài)和降雨條件下CDI指數(shù)，以及天然狀態(tài)下WDI指數(shù)，計算結(jié)果見表6.

表5 風(fēng)載作用下震裂山體穩(wěn)定性判別等級表Tab.5 Stability grade of shattering slope under wind loading

評價結(jié)果表明，石大關(guān)鄉(xiāng)斜坡在天然和降雨條件下偶有少量掉塊和落石，而在風(fēng)載作用下局部發(fā)生大規(guī)模的崩塌滾石的可能性較大，評價結(jié)果與真實情況相符.

此外，使用WDI評價指標(biāo)對四川藏區(qū)公路沿線21處震裂山體進行了評價，僅有2處評價結(jié)果與真實情況或定性評價結(jié)果有所區(qū)別，驗證了WDI評價指標(biāo)的合理性.且WDI評價指標(biāo)既考慮了降雨修正系數(shù)又考慮了風(fēng)載修正系數(shù)，使得評價體系相對更為完整.

表6 CDI、WDI評價計算表Tab.6 Evaluation calculation of CDI and WDI

4.2 臨界失穩(wěn)風(fēng)速評價法

針對綜合指標(biāo)評價法中WDI指數(shù)大于90的坡體，可以再以風(fēng)載作用下震裂巖體的啟動機理為基礎(chǔ)，由現(xiàn)場調(diào)查得出震裂縫切割出巖塊的形狀及與后緣巖壁之間的關(guān)系，確定其失穩(wěn)模式，對應(yīng)相應(yīng)失穩(wěn)模式的力學(xué)計算公式得出該巖體的臨界失穩(wěn)風(fēng)速，再通過對研究區(qū)風(fēng)的分布規(guī)律收集，確定出作用于該巖塊的實際風(fēng)速，兩風(fēng)速對比判定失穩(wěn)的可能性，評價流程如圖17所示.

綜合指標(biāo)評價法可以較為快速地對風(fēng)載作用下震裂山體的穩(wěn)定性作出評價，而臨界失穩(wěn)風(fēng)速評價法彌補了綜合指標(biāo)評價法中未考慮實際作用風(fēng)速和臨界失穩(wěn)風(fēng)速的不足.兩個方法共同運用，整體與局部分析相互結(jié)合，優(yōu)劣互補，使風(fēng)載作用下震裂山體穩(wěn)定性評價方法更加快速準(zhǔn)確.

圖17 穩(wěn)定性評價流程圖Fig.17 Flow chart of stability evaluation

5 結(jié) 論

震裂山體獨特的結(jié)構(gòu)特征以及河谷風(fēng)復(fù)雜多變的分布規(guī)律，使得風(fēng)載助推下震裂山體的失穩(wěn)機制迥異于通常重力環(huán)境.本文基于現(xiàn)場震裂山體的大量調(diào)查和河谷風(fēng)測量工作，開展了室內(nèi)風(fēng)洞試驗研究，歸納提出了風(fēng)載作用下震裂山體穩(wěn)定性評價方法，主要得出以下4點結(jié)論：

（1）沿原有結(jié)構(gòu)面拉開或者切割結(jié)構(gòu)面發(fā)育的震裂縫是汶川地震震區(qū)巖體震裂最為直觀的宏觀特征.裂縫普遍較新鮮，多無充填，局部還存在架空結(jié)構(gòu).區(qū)別于其他因素產(chǎn)生的松動巖體，震裂巖體具備形成時間短、規(guī)模大的特征.

（2）石大關(guān)鄉(xiāng)崩塌為風(fēng)載助推作用下震裂山體發(fā)生潰崩的典型實例，地震作用為崩塌的發(fā)生提供了物質(zhì)基礎(chǔ)，后期凍脹等風(fēng)化作用進一步劣化了坡體的穩(wěn)定性，河谷風(fēng)的助推作用為崩塌的發(fā)生提供了動力條件.此類災(zāi)害因震裂山體獨特的結(jié)構(gòu)特征，加上河谷風(fēng)的方向及大小復(fù)雜多變，就單個巖塊而言雖還是以傾倒、滑移破壞為主，但整體以潰崩式失穩(wěn)，表現(xiàn)為迅速解體—潰散—垮塌，散落的范圍極廣，崩落持續(xù)時間長，危害性極大，具極強突發(fā)性.

（3）巖塊的失穩(wěn)風(fēng)速與氣流方向與巖塊面垂線之間的夾角、巖塊的高寬比、后緣縫隙的寬度都成反比.就正方形、長方形巖塊而言，在風(fēng)載作用下以滑移、側(cè)向傾倒、順向傾倒失穩(wěn)模式為主.巖塊的重量、下伏結(jié)構(gòu)面的傾角、巖塊與巖壁之間縫隙的寬度都是決定是否發(fā)生順向傾倒的關(guān)鍵，只有縫隙非常小時，才有發(fā)生順向傾倒的可能.

（4）WDI綜合因子評價法和臨界失穩(wěn)風(fēng)速評價法綜合運用于風(fēng)載作用下震裂山體穩(wěn)定性的定量評價工作中，整體與局部分析相互結(jié)合，優(yōu)劣互補，使評價結(jié)果快速而準(zhǔn)確.

［1］ 許強，黃潤秋.“5·12”汶川大地震誘發(fā)大型崩滑災(zāi)害動力特征初探［J］.工程地質(zhì)學(xué)報，2008，16（6）：721-729.

XU Qiang，HUANG Runqiu.Kinetics charateristics of large landlidestriggered by MAY 12th wenchuan earthquake［J］.Journal of Engineering Geology，2008，16（6）：721-729.

［2］ YIN Y，WANG F，SUN P.Landslide hazards triggered by the 2008 Wenchuan earthquake， Sichuan，China［J］.Landslides，2009，6（2）：139-152.

［3］ 馮文凱，許強，黃潤秋.斜坡震裂變形力學(xué)機制初探［J］.巖石力學(xué)與工程學(xué)報，2009，28（增刊1）：3124-3130.

FENG Wenkai， XU Qiang， HUANG Runqiu. Deformation mechanics of the ruptures by earthquakes［J］.Chinese Journal of Rock Mechanics and Engineering，2009，28（Sup.1）：3124-3130.

［4］ 許強.四川省8·3特大泥石流災(zāi)害特點、成因與啟示［J］.工程地質(zhì)學(xué)報，2010（5）：596-608.

XU Qiang.The 13 August 2010 catastrophic debris flowsin Sichuan province：characteristics，genetic mechanism and suggestions［J］.Journal of Engineering Geology，2010，18（5）：596-608.

［5］ HUANG Runqiu， LI WeiLe. Post-earthquake landsliding and long-term impacts in the Wenchuan earthquake area，China［J］. Engineering Geology，2014，182：111-120.

［6］ 黃潤秋.汶川地震地質(zhì)災(zāi)害后效應(yīng)分析［J］.工程地質(zhì)學(xué)報，2011（2）：145-151.

HUANG Runqiu，After effect of geohazards induced by the wenchuan earthquake［J］.Journal of Engineering Geology，2011（2）：145-151.

［7］ 何思明，莊衛(wèi)林，張雄，等.都汶公路徹底關(guān)大橋橋墩抗?jié)L石沖擊防護研究［J］.巖石力學(xué)與工程學(xué)報，2013（S2）：3421-3427.

HE Siming，ZHUANG Weilin，ZHANG Xiong，et al，Research on rockfall impact prevention of chediguan bridge pier，duwhen road［J］.Chinese Journal of Rock Mechanics and Engineering，2013（S2）：3421-3427.

［8］ TANG C，ZHU J，LI W L，et al.Rainfall-triggered debris flows following the Wenchuan earthquake［J］. Bulletin of Engineering Geology and the Environment，2009，68（2）：187-194..

［9］ 余斌，楊永紅，蘇永超，等.甘肅省舟曲8.7特大泥石流調(diào)查研究［J］.工程地質(zhì)學(xué)報，2010，18（4）：437-444.

YU Bin，YANG Yonghong，SU Yongchao，et al. Research on the giant debris flow hazards in Zhouqu county，Gansu province on august 7，2010［J］.Journal of Engineering Geology.2010，18（4）：437-444.

［10］ 馮文凱，黃潤秋，許強，等.震裂斜坡形成機理及變形破壞模式研究［J］.水文地質(zhì)工程地質(zhì)，2009（6）：42-48.

FENG Wenkai，HUANG Runqiu，XU Qiang，et al. Study on formation mechanism and deformation failure models ofshatter slopes［J］. Hydrogeology and Engineering Geology，2009（6）：42-48.

［11］ 王運生，李渝生.岷江上游馬腦頂—兩河口段滑坡、崩塌形成的控制因素分析［J］.成都理工學(xué)院學(xué)報，2000（S1）：205-208.

WANG Yunsheng，LI Yusheng.Analysis on the controlling factorofformation oflandslides and avalanches in Manaoding-Lianghekou，upper reaches of Minjiang［J］. Journal of Chengdu University of Technology，2000（S1）：205-208.

［12］ 張倬元，王士天，王蘭生.工程地質(zhì)分析原理［M］.北京：地質(zhì)出版社，1994：181-185.

［13］ 羅永紅.地震作用下復(fù)雜斜坡響應(yīng)規(guī)律研究［D］.成都：成都理工大學(xué)，2011.

［14］ 郭沉穩(wěn)，姚令侃，段書蘇，等.汶川、蘆山、尼泊爾地震觸發(fā)崩塌滑坡分布規(guī)律［J］.西南交通大學(xué)學(xué)報，2016，51（1）：71-77.

GUO Chenwen，YAO Lingkan，DUAN Shusu，et al. Distribution regularities oflandslides induced by Wenchuan earthquake，Lushan earthquake and Nepal earthquake［J］. Journal of Southwest Jiaotong University，2016，51（1）：71-77.

［15］ 李德友，郭斌.四川省茂縣風(fēng)能資源初步分析［J］.高原山地氣象研究，2009（4）：70-73.

LI Deyou，GUO Bin.A preliminary analysis of wind energy resources in Mao county， Sichuan province［J］. Plateau and Mountain Meteorology Research，2009（4）：70-73.

［16］ 劉君，廖海黎，萬嘉偉，等.檢修車軌道導(dǎo)流板對流線型箱梁渦振的影響［J］.西南交通大學(xué)學(xué)報，2015，50（5）：789-795.

LIU Jun，LIAO Haili，WAN Jiawei，et al.Effect of guide vane beside maintenance rail on vortex-induced vibration of streamlined box girder［J］.Journal of Southwest Jiaotong University，2015，50（5）：789-795.

［17］ 袁進科，裴向軍.基于震裂損傷程度的公路崩塌識別方法［J］.公路交通科技，2013（3）：45-51.

YUAN Jinke，PEI Xiangjun.Identification of highway collapse based on degree of seismic damage［J］. Journal of Highway and Transportation Research and Development.2013（3）：45-51.

風(fēng)載作用下震裂山體崩塌機制及穩(wěn)定性評價方法

黃潤秋， 裴向軍， 羅 璟

Collapse Mechanism and Stability Evaluation of Shattered Slope under Wind Loading

HUANG Runqiu， PEI Xiangjun， LUO Jing
（State Key Laboratory of Geo-hazard Prevention and Geo-environment Protection，Chengdu University of Technology，Chengdu 610059，China）

To reveal collapse mechanisms of shattered slope under wind loading，the slope stability assessment system was built.Based on the field survey of shattered slopes and valley wind measurement，the large-scale wind tunnel tests in various conditions were conducted by using Shidaguan collapse as an example.The results indicate that Shidaguan collapse is a typical case of shattered slope collapsing in a sudden burst under wind loading.This mode of collapse are hardly controlled by structural plane because of the unique slope structure.Meanwhile，since the distribution of river valley wind is affected by the topographical elements，it becomes complicated and varied.Thus the sudden collapse mode undergoes fast breakdown-crumbling-collapse as a whole，quite distinct from slope failure mechanisms under general gravity force.The loosing and cracking of the rock masses is fundamental to the collapse，while the gradual weathering is the main factor.Besides，the strong rivervalley winds aggravate the collapse.Furthermore，the wind speed leading to rock masse instability is in inverse proportion to wind incident angle，as well as the height-width ratio of mass，and the width of the trailing edge crack.

earthquake；river valley wind；Shidaguan collapse；wind tunnel test；stability evaluation

黃潤秋，裴向軍，羅璟.風(fēng)載作用下震裂山體崩塌機制及穩(wěn)定性評價方法［J］.西南交通大學(xué)學(xué)報，2016，51（5）：958-970.

0258-2724（2016）05-0958-13

10.3969/j.issn.0258-2724.2016.05.020

P642.21

A

2016-06-23

國家自然科學(xué)基金資助項目（41572302）；四川省國土資源廳科學(xué)研究計劃資助項目（KJ-2015-18）

黃潤秋（1963—），男，博士，教授，博士生導(dǎo)師，研究方向為工程地質(zhì)和巖土工程，E-mail：hrq@cdut.edu.cn

（中文編輯：徐 萍 英文編輯：周 堯）