唐 然,范宣梅
(1.成都大學建筑與土木工程學院,四川 成都 610106;2.非飽和土力學特性及工程技術(shù) 四川省高校工程研究中心(成都大學),四川 成都 610106; 3.成都理工大學地質(zhì)災害防治與地質(zhì)環(huán)境保護國家重點實驗室,四川 成都 610059)
平推式滑坡在四川盆地紅層地區(qū)廣泛分布,主要指滑面傾角小于10°的平面滑動巖質(zhì)滑坡。張倬元等[1]首先提出巖體陡傾結(jié)構(gòu)面中充水的靜水壓力和滑面揚壓力是此類滑坡產(chǎn)生的主要驅(qū)動力。此后眾多學者對不同的滑坡案例從變形破壞機制、成因機理、關(guān)鍵層位的物理力學性質(zhì)和水理性質(zhì)以及地貌演化等方面對平推式滑坡開展了研究[2-18],取得了大量有益的成果。平推式滑坡的成因機理可歸納為:①一般具有巖層傾角近水平的軟硬互層的巖體結(jié)構(gòu);②坡體內(nèi)部含陡傾張性節(jié)理裂隙易于滲水;③下臥層中含泥質(zhì)類軟巖、泥化夾層或?qū)娱g剪切帶等軟弱夾層作為相對隔水層,在與地下水長期作用下形成滑帶;④具有良好的臨空條件。
垮梁子滑坡是平推式滑坡中的典型案例,是川中地區(qū)唯一的特大型巖質(zhì)平推式滑坡,引起了許多學者的關(guān)注。如吉隨旺等[19]研究認為軟硬巖層的差異應變是垮梁子滑坡形成的關(guān)鍵因素;翟國軍[20]詳細查明了垮梁子滑坡的發(fā)育特征,發(fā)現(xiàn)了降雨和地下水對該滑坡日常蠕變特性的影響規(guī)律;郭曉光等[21]將多級平推式滑坡變形模式分為牽引后退式和推動前進式,認為垮梁子滑坡屬于后者;陳思嬌[22]通過滑帶土環(huán)剪蠕變試驗對垮梁子滑坡日常蠕滑速率和變形機理進行了研究;張玲玲等[23]通過對孔內(nèi)加水加鹽前后相應的電位檢測,分析計算得到垮梁子滑坡坡體內(nèi)地下水的優(yōu)勢流向;呂紅賓[24]認為垮梁子滑坡坡體內(nèi)地下水流主要受控于區(qū)內(nèi)構(gòu)造裂隙發(fā)育而成的近于豎向的優(yōu)勢流動通道(集中流模式);冉佳鑫[25]通過分析垮梁子滑坡多年的降雨、位移、地下水水位監(jiān)測數(shù)據(jù),建立了GMD預警預報模型;Xu等[26]通過長期多源監(jiān)測數(shù)據(jù)分析了垮梁子滑坡地下水水位對降雨時間的響應和滑坡穩(wěn)定性變化;Ling等[27]利用高密度電阻率成像法結(jié)合鉆孔、監(jiān)測資料,查明了垮梁子滑坡地下水循環(huán)運移的特征。
以上研究成果驗證了降雨和地下水不僅是平推式滑坡的誘發(fā)因素,而且對滑坡尤其是滑帶形成演化的影響是必不可少的。因此,平推式滑坡大規(guī)?;瑒忧氨囟ň哂辛己猛杆缘钠麦w結(jié)構(gòu),某些結(jié)構(gòu)面是控制滑坡形成演化的核心,故分析平移式滑坡坡體的結(jié)構(gòu)特征并探究其成因就顯得尤為重要。
垮梁子滑坡位于四川省德陽市中江縣馮店鎮(zhèn),距離成都約80 km,處于川中地區(qū)倉山背斜北西翼,在滑坡周圍30 km范圍內(nèi)沒有大斷裂分布。該滑坡發(fā)育的地層為侏羅系蓬萊鎮(zhèn)組上段(J3p2),巖層產(chǎn)狀為340°~330°∠2°~5°,屬于典型的近水平巖層,巖性為泥巖與砂巖、粉砂巖互層,其中厚層和巨厚層砂巖占63%以上。
垮梁子滑坡前、后緣高差約為110 m,滑坡寬約為1 100 m,縱長為360~390 m,滑體平均厚度為50 m,最大厚度大于90 m,面積約為0.51 km2,體積約為2 550×104m3??辶鹤踊率状伟l(fā)生大規(guī)?;略?949年汛期,變形主要發(fā)生在滑坡Ⅰ區(qū),滑坡前緣形成帶狀分布的多個溜滑體、小型滑塌、濕地及泉點,后部形成3條拉陷槽地貌。1981年汛期垮梁子滑坡再次發(fā)生平推式滑動,滑坡I區(qū)北部變形最為劇烈,形成了滑坡Ⅱ區(qū)和滑坡Ⅲ區(qū),滑坡前緣分布帶狀隆起和多處鼓丘??辶鹤踊掳l(fā)育特征平面和工程地質(zhì)剖面圖,見圖1。
圖1 垮梁子滑坡發(fā)育特征平面圖和工程地質(zhì)剖面圖Fig.1 Development characteristics plan and engineering geology section of Kualiangzi landslide
通過調(diào)查統(tǒng)計垮梁子滑坡及周邊區(qū)域巖體內(nèi)節(jié)理的發(fā)育特征,其統(tǒng)計結(jié)果見表1、圖2和圖3。其中,圖2(a)~(f)為陡傾節(jié)理發(fā)育特征;圖2(g)為剖面節(jié)理發(fā)育特征。
由圖2可知:垮梁子滑坡及周邊區(qū)域巖體內(nèi)主要發(fā)育2組優(yōu)勢共軛節(jié)理,對應兩期不同方向的構(gòu)造擠壓作用(見表1),其中第一組共軛節(jié)理為NNW-NS走向(J1)及NWW走向(J2),大量剖面露頭顯示兩者相互交切形成共軛節(jié)理組[見圖3(a)、(c)];第二組共軛節(jié)理為NE走向(J3)及NW走向(J4)[見圖3(b)],低角度剖面節(jié)理傾向均為NNW-SSE向[見圖3(e)~(g)]。
表1 垮梁子滑坡及周邊區(qū)域巖體內(nèi)發(fā)育的共軛剪節(jié)理產(chǎn)狀特征
圖2 垮梁子滑坡及周邊區(qū)域巖體內(nèi)節(jié)理發(fā)育特征統(tǒng)計Fig.2 Statistical analysis of joint characteristics in rock mass of Kualiangzi landslide and its surrounding area(a)節(jié)理調(diào)查點D1~D4;(b)節(jié)理調(diào)查點D5~D6;(c)節(jié)理調(diào)查點D7;(d)節(jié)理調(diào)查點D8~D9;(e)節(jié)理調(diào)查點D10~D17;(f)節(jié)理調(diào)查點D1~D17;(g)剖面節(jié)理統(tǒng)計;(h)節(jié)理調(diào)查點分布圖
圖3 垮梁子滑坡及周邊區(qū)域巖體內(nèi)節(jié)理發(fā)育特征Fig.3 Joints photos of rock mass in Kualiangzi landslide and its surrounding area(a)第一組共軛節(jié)理;(b)第二組共軛節(jié)理;(c)張節(jié)理;(d)穿層張節(jié)理;(e)厚層砂巖內(nèi)低角度剖面節(jié)理;(f)、(g)厚層泥巖內(nèi)共軛剖面節(jié)理
其中,第一組節(jié)理中的J1節(jié)理普遍具正斷張性,一般開度為5~15 cm,沿走向延伸長度普遍大于5 m,沿傾向穿層特征較明顯[見圖3(d)],其形態(tài)多為樹枝狀或多級分叉狀,節(jié)理面粗糙不平,基本都具有劈理、片理的特征,且夾薄片、透鏡體,充填母巖礫石顆粒夾松散—中密黏粒,地下水活動時間越長黏粒含量越高,其特征詳見圖4。
圖4 垮梁子滑坡及周邊區(qū)域巖體內(nèi)J1張節(jié)理發(fā)育特征Fig.4 J1 open joints characteristics of rock mass in Kualiangzi landslide and its surrounding area(b)、(c)為某人工開挖洞室頂面張節(jié)理
川中地區(qū)屬于剛性基底,盆緣構(gòu)造帶對其產(chǎn)生的影響相對川西和川中地區(qū)更小,構(gòu)造變形也較為輕微,巖層產(chǎn)狀多為0°~5°,調(diào)查中發(fā)現(xiàn)的構(gòu)造變形破裂現(xiàn)象對古應力場的恢復有一定的指示作用。
2.2.1 第一期構(gòu)造變形特征
D7節(jié)理調(diào)查點處發(fā)現(xiàn)了一處斷裂(見圖5),附近出露的厚層砂巖內(nèi)節(jié)理特征與第一組共軛節(jié)理相符。從上盤砂巖斷層面具有的反階步特征以及砂巖層內(nèi)鄰近斷層面的裂隙與斷層面相交銳角的指示方向判斷該斷層面為逆沖斷裂。該處斷裂面產(chǎn)狀為315°~325°∠20°~22°,指示第一期構(gòu)造作用為NW-SE向擠壓,主應力方向與第一組共軛節(jié)理走向計算最大主應力方向基本一致。破裂帶上盤為巨厚層粉砂巖,下盤為紫紅色鈣質(zhì)泥巖,破碎帶可見紫紅色泥巖角礫、碎石與灰白色粉砂巖巖屑混合物膠結(jié)而成角礫巖[見圖5(b)、(c)],判斷其為第一期構(gòu)造運動形成。
圖5 逆沖斷裂指示NW-SE向擠壓Fig.5 Thrust fault indicating NW-SE tectonic compression
2.2.2 第二期構(gòu)造變形特征
多處構(gòu)造變形(見圖6)顯示為同期構(gòu)造成因,指示第二期構(gòu)造作用為NNW-SSE向擠壓,主應力方向與第二組共軛節(jié)理走向計算最大主應力方向基本一致。
圖6 指示NNW-SSE向構(gòu)造擠壓作用的構(gòu)造變形Fig.6 Tectonic deformation indicating NNW-SSE tectonic compression
D1節(jié)理調(diào)查點逆沖斷裂如圖6(a)所示,可見明顯的巖層沿破裂面的錯動現(xiàn)象,主破裂面產(chǎn)狀為351°∠67°,上盤可見強風化的粉砂巖與強風化泥巖的分界面,從下盤附近出露的泥巖判斷破裂面錯動距離約為2 m,破裂面內(nèi)為泥粉質(zhì)充填,膠結(jié)程度較低,形成時間應晚于D7節(jié)理調(diào)查點逆沖斷裂。該破裂面附近出露的第二組共軛節(jié)理[見圖3(b)]顯示兩者為同期構(gòu)造形成。D6節(jié)理調(diào)查點出露小向斜構(gòu)造[見圖6(b)],伴有層間剪切破碎變形帶,北側(cè)可測巖層產(chǎn)狀偏轉(zhuǎn)為173°∠21°,南側(cè)巖層產(chǎn)狀為355°∠19°。D13節(jié)理調(diào)查點出露朝SSE向逆沖的階梯狀破裂面[見圖6(c)],下盤為中層粉砂巖及薄層泥巖互層,上盤為厚層鈣質(zhì)砂巖,破裂面產(chǎn)狀為352°∠33°。
從節(jié)理特征及地表構(gòu)造形跡初步判斷,垮梁子滑坡及周邊區(qū)域巖體內(nèi)前后兩期節(jié)理分期與配套,見表2。
表2 垮梁子滑坡及周邊區(qū)域巖體內(nèi)節(jié)理分期與配套
四川盆地是陸殼內(nèi)擠壓應力作用下形成的沉積盆地,自印支期以來,盆地持續(xù)受到邊緣構(gòu)造帶階段性地向盆地內(nèi)部的擠壓作用,盆地的主要構(gòu)造形跡加強或定型于喜馬拉雅期,因此分析喜馬拉雅期構(gòu)造應力場和構(gòu)造形跡的演化可進一步論證構(gòu)造期次。
研究區(qū)地處川西及川中構(gòu)造區(qū)交界附近,區(qū)域構(gòu)造應力場及構(gòu)造形跡演化如圖7所示。區(qū)域內(nèi)主要的褶皺分屬不同的構(gòu)造系,形成時期有所區(qū)別,其中背斜A1~A8和向斜S1~S7屬于川中EW向構(gòu)造系,背斜A9~A12和向斜S8~S11屬于川西NE向構(gòu)造系。
相關(guān)研究發(fā)現(xiàn),喜馬拉雅中期(漸新世—中新世)四川盆地受近NW-SE向的壓應力作用[28-29],NE向構(gòu)造系應在此期間形成[見圖7(a)]。垮梁子滑坡區(qū)歸屬的倉山背斜[見圖7(a)中A6]西段軸部走向為NE向,推測形成于本階段。喜馬拉雅中期應為研究區(qū)經(jīng)歷的第一期構(gòu)造。
喜馬拉雅中晚期(上新世期間)川中地區(qū)以NS向構(gòu)造擠壓為主[28-29],形成川中EW走向構(gòu)造系,如圖7(b)所示??辶鹤踊聟^(qū)歸屬的倉山背斜東段軸部走向為近EW向,推測形成于本階段,最終其軸跡為向NNW凸出的復合構(gòu)造。垮梁子滑坡區(qū)處在兩期構(gòu)造應力場疊加部位,主應力方向為NNW,第二組共軛節(jié)理(J3、J4)在這個時期形成,同時將J1剪性結(jié)構(gòu)面改造為張性結(jié)構(gòu)面[見圖7(d)]。
圖7 研究區(qū)域構(gòu)造應力場及構(gòu)造形跡演化示意圖Fig.7 Regional tectonic stress field and tectonic trace evolution of the study area(a)喜馬拉雅中期(一期)區(qū)域構(gòu)造應力場;(b)喜馬拉雅中晚期(二期)區(qū)域構(gòu)造應力場;(c)四川盆地邊界及分區(qū);(d)共軛節(jié)理的形成及改造1.逆沖斷層;2.隱伏斷層;3.背斜;4.背斜編號;5.向斜;6.向斜編號;7.未形成的構(gòu)造形跡;8.滑坡點位;9.一期構(gòu)造主應力方向;10.二期構(gòu)造主應力方向;11.構(gòu)造分區(qū)界線;A1.老君庵背斜;A2.八角背斜;A3.金華鎮(zhèn)背斜;A4.南充背斜;A5.建中背斜;A6.倉山背斜;A7.蓬萊鎮(zhèn)背斜;A8.攔江背斜;A9.簡陽背斜;A10.龍泉寺背斜;A11.陽化場背斜;A12.人和場背斜;S1.玉皇廟向斜;S2.金孔向斜;S3.金家場向斜;S4.西山向斜;S5.元興場向斜;S6.河邊場向斜;S7.蟠龍河向斜;S8.胡家場向斜;S9.飛龍寺向斜;S10.中天場向斜;S11.中和場向斜
第二期構(gòu)造后區(qū)域上經(jīng)歷了新構(gòu)造運動,大量學者通過研究認為,新構(gòu)造運動總體上控制和影響了現(xiàn)今地形、地貌的形成和發(fā)展,地表水系展布特征能反映出區(qū)域新構(gòu)造應力場[30],有學者利用水系展布反演新構(gòu)造應力場[31-32],取得了良好的效果。本文采用基于GIS技術(shù)改進的反演法[33],選取垮梁子滑坡及周邊區(qū)域水系[見圖8(a)],采用GIS技術(shù)擬合得到每條河流展布方向讀取實際長度[見圖8(b)],并統(tǒng)計不同角度范圍內(nèi)的累計長度得到水系展布玫瑰花圖,進而得出兩個優(yōu)勢方向為350°~360°和310°~320°[見圖8(c)],兩者銳角平分線方向?qū)鲏簯Ψ较?,因此垮梁子滑坡及周邊區(qū)域新構(gòu)造應力場主壓應力方向為330°~340°,與第二期構(gòu)造主應力方向相比產(chǎn)生了少量逆時針轉(zhuǎn)動,仍然與J1節(jié)理小角度相交,保持了原有節(jié)理的力學性質(zhì)。
圖8 垮梁子滑坡及周邊區(qū)域新構(gòu)造應力場反演Fig.8 Inversion of neotectonic stress field in Kualiangzi landslide and its surrounding area(a)垮梁子滑坡及周邊區(qū)域水系展布;(b)擬合以后水系展布;(c)水系展布玫瑰花圖及主應力方向
調(diào)查和鉆探結(jié)果顯示,J1節(jié)理的張性特點在斜坡內(nèi)部也普遍存在,結(jié)構(gòu)面上附著的鈣膜[見圖9(a)~(c)、(f)~(i)]和深部膠結(jié)物[見圖9(d)、(e)]說明地下水開始活躍的時間遠遠早于垮梁子滑坡大規(guī)模滑動的時間。
圖9 垮梁子滑坡及周邊區(qū)域巖體內(nèi)張節(jié)理地下水活動痕跡Fig.9 Groundwater activity traces of the open joint of rock mass in Kualiangzi landslide and its surrounding area(a)、(b)拉陷槽內(nèi)崩落塊石上附著鈣膜;(c) 鉆孔揭露滑體深部結(jié)構(gòu)面附著鈣膜;(d)、(e)鉆孔揭露滑體深部膠結(jié)物;(f)、(g)滑坡后山結(jié)構(gòu)面內(nèi)充填物含鈣質(zhì)結(jié)核;(h)、(i)滑坡后山結(jié)構(gòu)面附著鈣膜
J1節(jié)理走向與垮梁子滑坡拉陷槽總體走向相符,同時拉陷槽后壁的主要結(jié)構(gòu)面均為J1組節(jié)理。因此,暴雨作用下滑坡產(chǎn)生平推式滑動時,這組導水張節(jié)理又成為產(chǎn)生靜水壓力的儲水裂縫。J1節(jié)理是在歷史地質(zhì)構(gòu)造作用下奠定了張性特征,在后期斜坡的形成過程中,隨著應力的釋放,J1節(jié)理還會受到一定的改造作用。
為了模擬新構(gòu)造運動中內(nèi)、外動力地質(zhì)作用下斜坡形成卸荷過程對節(jié)理的改造作用,以垮梁子滑坡所在的條形山脊作為模擬對象,采用三維離散元方法建立了數(shù)值模型,如圖10所示。從地質(zhì)體原始形態(tài)開始分4步下切,直至初步形成現(xiàn)今地貌形態(tài),整個下切過程對模型南北邊界施加NNW方向水平構(gòu)造應力,東西邊界設為變形約束邊界。數(shù)值模型模擬的條形山脊NS方向長為2.5 km,EW向長為1.6 km,構(gòu)造應力設為水平,方向為340°,構(gòu)造應力大小參考相關(guān)文獻[34]對威遠、安岳地區(qū)構(gòu)造應力場的數(shù)值模擬結(jié)果取10 MPa。數(shù)值模型中巖土材料的物理力學參數(shù)取值,見表3[22]。
圖10 垮梁子滑坡數(shù)值模型建立Fig.10 Numerical modeling process for Kualiangzi landslide(a)模型邊界與橫剖面;(b)數(shù)值模型初始形態(tài)與斜坡形態(tài)
表3 數(shù)值模型中巖土材料的物理力學參數(shù)取值[22]
利用數(shù)值模型模擬得到斜坡形成后的變形破壞特征,見圖11。
圖11 數(shù)值模擬斜坡形成后的變形破壞特征Fig.11 Numerical simulation of deformation and failure characteristics after slope formation(a)、(b)、(c)斜坡山脊表面裂縫;(d)、(e)斜坡內(nèi)部裂縫
由圖11可知:斜坡形成后在條形山脊頂部出現(xiàn)大量追蹤NNW走向的張裂縫,裂縫圈閉形態(tài)大致為向東突出的弧形,對應滑坡拉陷槽Ⅰ后側(cè)壁;南側(cè)靠沖溝部位出現(xiàn)剪切裂縫,裂縫兩側(cè)塊體可見明顯的剪切位錯變形,對應滑坡左側(cè)邊界及受剪切作用產(chǎn)生的滑塌區(qū)[見圖11(a)~(c)],分布范圍與垮梁子滑坡邊界大致相同;地表變形破壞多集中在垮梁子滑坡Ⅰ區(qū),說明滑坡Ⅰ區(qū)比滑坡Ⅱ區(qū)具備更良好的垂直優(yōu)勢滲透通道,印證了現(xiàn)實中滑坡Ⅰ區(qū)比滑坡Ⅱ區(qū)更早產(chǎn)生平推式滑動的緣由;從模型1-1′剖面[見圖11(d)、(e)]來看,在山脊深部形成了一條垂向的裂隙帶,組成垂向滲透網(wǎng)絡通道,說明斜坡在形成初期其結(jié)構(gòu)就已具備一定的滲水條件。
(1) 垮梁子滑坡及周邊區(qū)域巖體內(nèi)共發(fā)育2組共軛節(jié)理,分別在一期構(gòu)造NW-SE向擠壓和二期構(gòu)造NNW-SSE向擠壓作用下形成。其中,NNW-NS走向的J1組節(jié)理呈顯著的張性特點,該組張節(jié)理是坡體結(jié)構(gòu)中主要的導水結(jié)構(gòu)面,在坡體內(nèi)部廣泛分布,具有地下水長期活動的痕跡,且該組結(jié)構(gòu)面在滑坡形成演化過程中起到了重要作用。
(2) 通過分析區(qū)域構(gòu)造應力場演化后認為,是歷史地質(zhì)構(gòu)造作用奠定了J1組節(jié)理張性的特征。J1組節(jié)理形成于喜馬拉雅中期(一期構(gòu)造),在喜馬拉雅中晚期(二期構(gòu)造)被改造成為了張節(jié)理,新構(gòu)造運動期間保持了其張性節(jié)理的性質(zhì)。
(3) 數(shù)值模擬結(jié)果顯示:在新構(gòu)造應力場和河谷下切的內(nèi)、外動力聯(lián)合作用下,斜坡表面形成大量的裂縫,其分布范圍與滑坡邊界大致相近;斜坡山脊部位產(chǎn)生大量追蹤NNW走向的張裂縫,并向下延伸直至下伏的軟弱泥巖層,組成一條垂向的地下水滲透網(wǎng)絡,說明拉陷槽發(fā)育在山脊部位并形成如此大規(guī)模的滑坡具有一定的必然性,且斜坡在形成初期,坡體結(jié)構(gòu)就已具有較好的滲透性。