Progressive Failure Mode and Stability Reliability of Strain-Softening Slope
-
摘要: 针对现有考虑应变软化效应的边坡渐进破坏分析模型中计算假定条件不符合实际情况的不足,以三峡库区赵树岭滑坡为研究对象,考虑地下水位波动和地震力的影响,提出用于分析多场耦合条件下应变软化边坡渐进破坏模式及稳定性可靠度的方法.结果表明,地下水位波动和地震力会不同程度地影响滑坡的渐进破坏模式和破坏概率.滑坡在145 m、175 m水位和库水位从175 m陡降为145 m三种工况下整体稳定,但分别有14.119%、20.266%和21.797%的概率发生局部渐进破坏;在烈度为Ⅶ的地震力作用下同时考虑3种蓄水工况,该滑坡分别有34.067%、38.061%和38.405%的概率发生整体渐进破坏;根据滑坡渐进破坏模式指出最佳加固位置应在沿江大道前沿.该分析方法具备可靠性,可为边坡的渐进破坏模式及稳定性评价研究提供参考.Abstract: Aiming at the shortcomings of the calculation assumptions proposed in the existing slope progressive failure analysis model considering the effect of strain-softening, the Zhaoshuling landslide in the Three Gorges Reservoir area is taken as the research object, and the influence of groundwater level fluctuations and seismic forces is considered to propose the progressive failure mode and stability reliability analysis method of strain-softening slope under multi-field coupling conditions.The results show that the groundwater level fluctuation and seismic force would affect the progressive failure mode and failure probability of the landslide to various degrees.The landslide was stable at 145 m, 175 m water level and the steeply reduced reservoir water level from 175 m to 145 m, but it had a local progressive failure probability with 14.119%, 20.266% and 21.797% respectively. And the landslide had a total progressive failure probability with 34.067%, 38.061% and 38.405% respectively under the action of seismic force with intensity VII, also three kinds of water storage conditions were considered. Moreover, according to the progressive failure mode of landslide, it is pointed out that the best reinforcement position should be at the front of Yanjiang Avenue. The analysis method is reliable and can provide a reference for the progressive failure mode and stability evaluation of slopes.
-
桐柏-大别碰撞造山带介于中朝克拉通与扬子克拉通之间, 向西与秦岭造山带相连, 向东与苏鲁超高压带乃至朝鲜半岛的临津江带相通.至今, 关于大别超高压变质岩形成的构造动力学背景, 已有趋于一致的认识, 即将超高压变质岩理解为是扬子克拉通与中朝克拉通于印支期斜向碰撞的产物[1~13].近年来研究表明, 在桐柏-大别地块内, 现今所观察到的构造几何学样式, 类似于北美西部的变质核杂岩带, 以总体具穹隆形态及多层拆离滑脱带的发育为特征[12, 14].现在所看到的高压、超高压变质岩及非高压岩石的分布是由碰撞期后伸展构造格架控制的, 各构造岩石单元间由多层伸展拆离带分隔[12].
桐柏-大别碰撞造山带具有长期的构造-热历史, 其中, 桐柏、大别核部杂岩内表壳岩系的多期变形, 晋宁期广泛的硅铝壳部分熔融和再造作用, 是印支期陆陆碰撞前的主要地质事件[15, 16].中生代大规模的岩浆岩体就位、陆内变形及伸展塌陷构造等, 则是印支期陆陆碰撞后的热-构造事件[15, 16].研究桐柏-大别碰撞造山带的组成, 不仅不能按传统的地层方法划分, 也不能仅考虑组成物质在二维空间的变化, 而必须结合造山带的结构和演化, 从三维空间予以厘定.研究表明, 桐柏-大别造山带的组成与演化, 与扬子克拉通的关系更为密切一些[12], 不过, 现今桐柏-大别造山带的组成与结构, 主要是印支期碰撞及高压、超高压变质期后伸展构造和中新生代热-构造演化的结果.所以, 除了燕山期及其后的岩浆活动和盆地堆积产物以外, 桐柏-大别碰撞造山带的基本组成主要包括核部杂岩(CC) 单元、超高压(UHP) 单元、高压(HP) 单元、绿帘-蓝片岩(EB) 单元和沉积盖层(SC) 单元等, 此外还有一些镁铁质和超镁铁质岩体, 各构造岩石单元间分别由下伸展拆离带、中伸展拆离带、上伸展拆离带和顶伸展拆离带所分隔(图 1).
1. 基本组成
1.1 核部杂岩单元
核部杂岩单元主要分布于大别造山带的中部和北部及桐柏造山带的核部.它主要由桐柏杂岩和大别杂岩组成, 包括变质表壳岩系、变质镁铁质岩石和变质花岗岩.其中变质表壳岩系和变质镁铁质岩石主要包括斜长角闪岩、黑云斜长片麻岩、变粒岩及磁铁石英岩、夕线榴片麻岩、基性及酸性麻粒岩和大理岩等.它们的原岩多是前寒武纪的表壳岩和侵入其中的基性、超基性岩体, 具有麻粒岩相-高角闪岩相变质作用及多期褶皱变形特征, 并经历了强烈的部分熔融和混合岩化作用.它们多作为大小不一的残块包裹于变质的花岗质岩石之中(图 2), 在大别杂岩中所占的比例很少, 在桐柏杂岩中所占的比例更少.核部杂岩单元中变质花岗质岩石以花岗闪长质片麻岩和花岗质片麻岩为主, 还有少量的英云闪长质片麻岩, 其w (SiO2) 为63.34%~75.36%.研究表明, 桐柏杂岩和大别杂岩中变质花岗质岩石的地球化学特征是一致的或相似的, 具有可对比性[17], 它们主要是古老地壳在晋宁期受到强烈再造和部分熔融的产物[15, 18~20].此外, 核部杂岩中还有大量的燕山期花岗质和镁铁质-超镁铁质岩体就位, 因此, 真正的古老结晶基底变质岩石保留很少, 仅以大别山的木子店和黄土岭等地保存较好(图 2).在核部杂岩分布范围内, 至今尚未发现典型超高压-高压岩石.
1.2 超高压单元
超高压单元主要分布于大别造山带南部的潜山、太湖、岳西和英山一带, 西部的麻城—新县一带以及北部的桐城一带, 位于核部杂岩单元与高压单元之间.主要岩石组合为英云闪长质片麻岩、面理化(含榴) 花岗岩和榴辉岩, 还有少量的大理岩、硬玉石英岩及镁铁质岩石等.超高压榴辉岩以含柯石英为特征, 多以透镜状、扁豆状或团块状产于片麻岩中, 少量产于大理岩中, 其形成温压条件为t=750~800 ℃, p > 2.8 GPa[21].据矿物组合和岩石结构构造特点, 榴辉岩大体可分为两种类型, 即块状榴辉岩和面理化榴辉岩.前者的峰期变质矿物组合主要为石榴石+绿辉石+金红石, 块状或具弱面理; 后者的矿物组合上除了石榴石、绿辉石和金红石外, 一般还含有蓝晶石、多硅白云母、黝帘石或滑石等, 并发育明显的面理和线理组构.含榴辉岩透镜体的大理岩经研究证实曾经历过超高压变质作用[3].榴辉岩的围岩英云闪长质片麻岩主要为黑云斜长片麻岩(俗称“超高压片麻岩”), 含不等量的角闪石、绿帘石和石榴石、在化学成分上主要相当于英云闪长岩[22] (图 3), 与桐柏、大别核部杂岩中的花岗质片麻岩在化学成分上有明显区别, 后者主要相当于花岗闪长岩和花岗岩(图 3).英云闪长质片麻岩中常含有具超高压矿物组合的榴辉岩或退变质榴辉岩的残块及阴影状残余, 且榴辉岩与英云闪长质片麻岩呈渐变关系.超高压单元中面理化(含榴) 花岗岩广泛出露, 在新县田铺、麻城福田河、四道河、英山东冲河、潜山双河、三祖寺、岳西碧溪岭以及桐城桃冲等地都可见及, 在整个单元中占有很大的比例.这些面理化含榴花岗岩常常包容各种英云闪长质片麻岩乃至榴辉岩和退变榴辉岩, 或穿插于它们之中[22].在一些条带状或面理化榴辉岩或退变榴辉岩中, 花岗质细脉或条带穿插其中, 并与其寄主岩一起褶皱; 而在一些块状榴辉岩或榴辉质斜长角闪岩中, 这些脉体或条带常呈枝状贯入, 显示了部分熔融的迹象.面理化(含榴) 花岗岩的主要矿物组合为钠-奥长石+微斜长石+角闪石+黑云母+白云母+石英+绿帘石+石榴石, 在化学成分上相当于奥长花岗岩和花岗岩[18] (图 3).在超高压单元中, 从片麻岩到面理化花岗岩, 呈现出奥长花岗岩的演化趋势[22], 与核部杂岩中花岗质岩石的演化趋势(钙碱性演化趋势) 形成鲜明的对照(图 3).面理化花岗岩中常含有榴辉岩、退变榴辉岩和片麻岩的包体.无论是片麻岩还是面理化花岗岩, 都具有统一的缓倾斜区域性面理, 是伸展流动作用下形成的面理, 表现为由含榴辉岩或退变榴辉岩透镜体的片麻岩成分层堆垛成“假单斜”, 而榴辉岩体的最大扁平面平行区域性缓倾斜面理并在垂向上堆垛.
超高压单元内榴辉岩体常显示不同程度的退变质, 转变为榴辉岩质斜长角闪岩和榴辉岩质片麻岩.在榴辉岩的围岩片麻岩的副矿物中发现绿辉石[23], 以及在锆石中发现柯石英及绿辉石等都是明证[24, 25].在有些较大榴辉岩体产出地可追索出从榴辉岩-角闪石化榴辉岩-榴辉岩质斜长角闪岩-斜长角闪岩(有时可保留角闪石+斜长石的后成合晶)- (含榴) 黑云角闪斜长片麻岩(所谓的超高压片麻岩) 的逐渐过渡, 而且这几种岩石的稀土模式具有一定的相似性和继承性, 暗示了它们之间的演化关系.
除了超高压榴辉岩以外, 超高压单元中还产有少量的超高压镁铁质岩石.它们主要是一些经历了超高压变质的超基性岩, 以变质的方辉橄榄岩和二辉橄榄岩等为主, 如大别北部的饶钹寨和中部的碧溪岭及石马等地所见.这些超高压镁铁质岩石常常与超高压榴辉岩共生, 或将超高压榴辉岩包于其中.此外, 超高压单元中还可见到硬玉石英岩(如在双河地区) 和蓝晶石石英岩(如在新县田铺地区), 它们主要以透镜状产在英云闪长质片麻岩和奥长花岗质片麻岩中, 并与榴辉岩或含榴辉岩的大理岩共生.
超高压榴辉岩及退变质榴辉岩透镜体一般成群分布, 它们在垂向上平行堆垛, 其三维空间上的排列格式及组合特征, 主要取决于与包围它们的围岩间的流变性(或强度) 差及应变体制(纯剪、单剪或一般剪切).流变性差明显的, 如大理岩内的榴辉岩块体, 其形体一般较小, 与大理岩构成残斑-基质流变学格式[12]; 流变性差较小的, 如黑云角闪片麻岩内的榴辉岩或退变质榴辉岩块体, 其形体一般较大, 长轴较长, 与围岩构成布丁-基质流变学型式[12].其中, 榴辉岩及退变质榴辉岩体起着强相的角色.当流变性差很小或应变很大时, 局部会出现围岩与榴辉岩构成假层理或条带的现象.
值得指出的是, 在“北大别”也有超高压-高压变质岩石的分布.魏春景等[26]和徐树桐等[27]在桐城西部大麻岩及百丈岩等地区发现了超高压-高压榴辉岩及其退变质产物.最近, 笔者自岳西县水吼岭开始, 经龙井关、桃冲、牯牛背水库向北追索, 分别在罗家冲、胡湾及舒城、桐城和潜山三县交界地区的百丈岩一带, 观察和发现较多的榴辉岩及分布几百平方千米的超高压单位和/或高压单元岩石组合(图 4).除此之外, 于岳西县北部黄尾地区及霍山县饶钹寨地区, 也观察到变形的方辉橄榄岩与榴辉岩相岩石特有的空间组合关系, 其特征与太湖毛屋及与苏鲁地区的基本一致, 只是后期高温变质叠加更为强烈一些, 保留了麻粒岩相的退变质组合.大别山北部罗家冲及百丈岩地区的榴辉岩主要由石榴子石+绿辉石+金红石组成, 矿物变质反应结构、塑性变形及形组构发育, 代表超高压-高压韧性剪切带的残余形迹.
综上所述, 超高压单元主要由经过超高压变质作用的大陆壳及幔源超镁铁质岩石、退变质的超高压变质岩石及减压退变质和部分熔融作用形成的片麻岩及面理化含榴花岗岩组成, 构成一个8~10 km厚的楔状岩片, 在核部杂岩的南部和北部均有出露.它们与下伏的主要由高温变质的大别杂岩构成的核部杂岩带之间, 以下滑脱带相隔.在有些区段, 因地壳薄化及伸展拆离作用影响, 缺失超高压单元岩石, 以致由高压单元岩石直接覆于核部杂岩单位之上[12].
1.3 高压单元
高压单元夹于超高压单元和蓝闪绿片岩单元之间, 在大别山主要分布于河南罗山、湖北大悟、红安和安徽宿松等地, 大致相当于原来所划的宿松群、红安群和苏家河群浒湾组所在的范围; 在桐柏地区主要分布于桐柏山的两侧.在大别山, 该单元主要由白云钠长片麻岩、钠长绿帘角闪岩及以透镜状产于其中的榴辉岩组成, 还有大量的面理化(含榴) 花岗岩和少量的变质磷块岩和大理岩.高压榴辉岩主要有石榴石和绿辉石组成, 常见蓝闪石和多硅白云母, 其形成温压条件为t=450~680 ℃, p=1.2~1.5 GPa[21, 28].高压榴辉岩也经历了不同程度的退变质作用, 可见到由榴辉岩-榴闪岩石榴角闪岩-绿帘角闪岩、钠长绿帘角闪岩乃至蓝闪绿片岩和绿片岩的连续退变质系列.在高压单元中, 可较清楚地辨认高压榴辉岩与(钠长) 绿帘角闪岩之间的演化关系, 但高压榴辉岩与白云钠长片麻岩间的关系还有待进一步研究.不过, 从高压单元中面理化含榴花岗岩的不断识别和其中榴辉岩残留体的发现(如红安河口地区及苏家河地区), 已为此提供了有意义的启示.
与大别造山带的高压单元类似, 桐柏地区的高压单元也主要由4部分岩石组合构成: (1) 榴辉岩或退变质榴辉岩, 局部伴有变形的蛇纹岩化橄榄岩; (2) 经受高压变质的沉积岩及火山岩, 如白云石大理岩、石墨片岩、石英岩及白云钠长片麻岩等; (3) 由榴辉岩强烈退变质及剪切变形而成的含榴二云石英片岩、含榴绿帘斜长角闪片岩、片麻岩等; (4) 面理化含榴或不含榴的花岗岩.地质部门在新近的1∶5万地质测量过程中, 分别称作为肖家庙岩组、马鞍山岩组、鸿仪河岩组及丘沟岩组等.其实, 它们在岩石组合、变质变形特征方面都是一致的, 在三维空间上都是相互联通的.经过区域上的追索和观察证明, 桐柏地区的高压单元, 向东与大别地区出露的高压变质单位是衔接的, 其构造边界也可逐一连接和对比; 向西被南襄盆地的陆相盆地沉积掩盖, 但据钻孔及地球物理资料, 陆相沉积物基底岩石及主要构造边界, 均可与盆地东西两侧山区的岩石及构造对比, 因而, 桐柏山区的高压单元, 有越过南襄盆地向西延展的趋势.
1.4 绿帘蓝片岩单元
绿帘蓝片岩单元是形成压力最低的一个高压变质带, 分布于桐柏-大别碰撞造山带的南侧, 主要由绿帘蓝闪片岩(变质基性火山岩)、蓝闪白云钠长片岩(变质酸性火山岩)、蓝闪白云石英片岩(变质泥质岩) 和蓝闪大理岩(变质碳酸盐岩) 以及绿片岩、白云钠长片岩及白云石英片岩等组成.大致包括了原来所划分的张八岭群、随县群、耀岭河群及武当山群所在的范围.
从绿帘蓝片岩带的绿帘蓝闪片岩和蓝闪白云钠长片岩等的矿物组合变化及矿物环带等的特征说明, 该带至少经历了3期变质作用演化, 即从低绿片岩相→绿帘蓝闪片岩相的进变质作用→绿片岩相、低绿片岩相的退变质作用过程.其峰期变质作用的温压条件为t=350~450 ℃, p > 0.7 GPa[17].特别应引起注意的是, 在绿帘蓝片岩带中还可见有少量残留的榴辉岩透镜体.同样, 在部分高压榴辉岩中仍可见到绿帘蓝片岩相退变质作用的叠加, 表现为在榴辉岩中有蓝闪石和绿帘石变斑晶的形成, 而且在这两个矿物中含有石榴石和绿辉石的包体.这些都暗示了高压榴辉岩与绿帘蓝片岩间的转化关系.
1.5 沉积盖层
大别碰撞造山带内所保存的盖层岩系(SC), 由于构造揭顶作用及侵蚀破坏, 仅在造山带的南缘有残留露头, 另在上述各单元的顶部也偶见出露.这些沉积盖层由晚震旦纪至三叠纪沉积岩组成.据古地理分析资料, 震旦纪至中三叠世阶段, 整个大别和苏鲁地区, 都曾有沉积作用记录, 该阶段还不存在大面积的古陆[29].沉积岩(局部夹火山岩) 的特征与扬子克拉通盖层有亲缘性.据有限的盖层露头构造研究, 其内部变形比较复杂, 具褶皱逆冲带性质, 以顶拆离带分别与下覆的高压或绿帘蓝片岩单元的岩石接触, 底部可见缺层现象, 压溶现象发育.运动学标志指示沿着顶拆离带, 曾发生过大规模的正向滑脱运动[12].盖层只经历成岩和低级-极低级变质, 未卷入超高压和高压变质作用, 不属于超高压和高压变质带, 但盖层确是超高压和高压峰期变质作用后形成的伸展构造框架的重要要素之一, 对认识超高压和高压变质岩石的形成与折返动力学过程有重要的参考坐标意义.
1.6 镁铁质及超镁铁质岩石
大别碰撞造山带内, 尤其是大别山北部广泛分布大小不一的镁铁质及超镁铁质岩石块体(图 4).依据它们的矿物组合、变形变质特点及与围岩的接触关系, 可分为两大类: 一类是饶钹寨、碧溪岭和石马等地的变形的方辉橄榄岩、纯橄榄岩组合为代表, 如前所述, 它们与榴辉岩相岩石有相同的变形变质及几何学特征.如饶钹寨两个垂向上叠置的方辉橄榄岩扁平透镜体的长轴平行区域拉伸线理, 与区域上榴辉岩透镜体形态及堆垛格式一致.地球化学研究表明[30], 这些超镁铁质岩的稀土模式为LREE富集型, 不同于大洋地幔, 其w (87Rb) /w (86Rb) (0.706 8~0.710 4) 和ε (Nd) 值(+1.2~-7.1) 表明不是来自亏损的地幔源区.另一类镁铁质岩石是辉石岩、角闪辉石岩及辉长岩组合, 多为宏观上未变形的侵入体, 与围岩有清楚的侵入接触关系, 并含有围岩捕虏体, 如岳西小河口岩体及霍山祝家铺岩体等.同位素年代学资料表明后一类镁铁质-超镁铁质岩体是燕山期就位的[31].这些镁铁质及超镁铁质的岩石地球化学特征[32]及野外地质体间几何关系、变形行为, 均不具变质蛇绿混杂岩的特征, 因而, 不能作为三叠纪碰撞时期古缝合线的标志.尽管大别碰撞造山带内不存在变质蛇绿混杂岩带, 但熊店、苏家河及浒湾一带榴辉岩中所显示的加里东期同位素年代学数据及ε (Nd) 值(可达-20) 表明, 在大别山很可能保存有在加里东时期华北与扬子地块对接时被消减的古洋盆的残片.
除了上述主要构造岩石单元外, 桐柏-大别碰撞造山带中还有很多燕山期及其后的岩浆活动和盆地堆积产物, 在此不再赘述.但应提及的是, 中生代燕山期大规模的岩浆岩体就位, 是印支期陆陆碰撞后重要热-构造事件的反映, 它在桐柏-大别碰撞造山带留下了深刻的烙印, 这在大别山东部尤为显著.同位素示踪和地球化学研究表明, 大别山内的燕山期花岗岩, 不管是产出于核部杂岩, 还是产出于超高压或高压单元中, 其源区都来自于核部杂岩[33], 这也从另一侧面证明核部杂岩在空间上是位于高压、超高压单元之下的.据Wang等[34]对深反射地震剖面的解析, 大别超高压岩石主要集中于9 km以上的地壳.
2. 基本结构
总体观之, 桐柏-大别碰撞造山带的基本结构格局类似于一个大型的北美科迪勒拉型变质核杂岩(metamorphic core complex)[21].主要结构要素包括核(杂岩) 带、大型伸展拆离带及夹持的构造岩片和沉积盖层三部分[17].
现今在桐柏-大别碰撞造山带观察到的区域构造几何图像, 主要是由角闪岩相变质条件下形成的透入性成分层(具S-L组构) 区域性几何学控制的.在大别山地区, 成分层及其面理迹线勾绘成以罗田为中心的大型不对称穹窿构造, 成分层及面理向周围缓倾斜, 但拉伸线理则分别向北西及南东方向倾伏.穹窿的核部由大别变质杂岩组成, 顶部及翼部由超高压和高压变质岩组成, 由于中生代岩体的侵入及断裂的切错破坏, 穹窿构造受到肢解, 但基本几何学图像仍很清晰.在超高压、高压及绿帘蓝片岩单元中, 榴辉岩及退变质榴辉岩体、变形超镁铁质-镁铁质岩体最大扁平面以及各种平卧及斜卧无根褶皱轴面、剑鞘褶皱的xy面等, 都是与成分层平行的.利用对称布丁、共轭韧性剪切条带及脉体变形特征等应变及运动学标志估算, 片麻岩及花岗质岩石垂向压扁量多为70%~80%, 榴辉岩及退变质榴辉岩体垂向压扁量仅为50%, 区域上具纯剪(共轴) 及一般剪切应变体制.
在桐柏地区, 区域上形成一NWW-SEE向延长的短轴背形构造, 在背形两翼, 面理分别向北东及南西倾斜, 倾角一般在20°~35°之间.在背形转折端, 如太白顶一带, 面理近水平.背形的东南和北西倾伏端, 面理分别向南东及北西倾斜.矿物及拉伸线理特别是窗棂构造极为发育, 趋向NWW及SEE, 倾伏角多在5°~25°之间.总之, 桐柏地区的区域构造样式与大别地区的区域构造样式非常相似, 都具有一个变质核杂岩的特征.其几何形态为开阔的短轴背形, 结构上其核部单位由桐柏杂岩构成, 沉积盖层只遭受极低级变质作用, 高压变质单元界于下、中复合拆离带与上拆离带之间.总体是在垂向缩短及NWW-SEE方向近水平的伸展体制下形成的.
3. 讨论与结论
现今桐柏-大别碰撞造山带的组成与结构, 主要是印支期碰撞及高压、超高压变质期后伸展构造和中新生代热-构造演化的结果.在组成上, 除了燕山期及其后的岩浆活动和盆地堆积产物以外, 主要包括核部杂岩单元、超高压单元、高压单元、绿帘-蓝片岩单元和沉积盖层单元等, 此外还有一些镁铁质和超镁铁质岩体残留或侵入其中.各构造岩石单元间由多层伸展拆离带分隔.前人将大别碰撞造山带划分为“南大别”和“北大别”, 在组成和结构的内涵上都不够确切.从三维空间看, 经历过超高压和高压变质作用的超高压单元、高压单元和绿帘蓝片岩单元作为不同的岩片夹持于核部杂岩和沉积盖层之中, 即超高压、高压和绿帘蓝片岩单元在空间上都位于核部杂岩之上, 而位于沉积盖层之下, 其分布格局受碰撞期后伸展构造格架所制约.大别碰撞造山带的整体结构样式类似于北美西部的变质核杂岩带, 即以总体具穹隆形态及多层拆离滑脱带的发育为特征, 构成了以罗田和桐柏山为核部的两个穹隆.
至今, 关于三叠纪中-朝与扬子克拉通间碰撞缝合线的位置有几种不同的意见.在大别地区, 有些学者主张水吼-五河断裂带是UHP变质地体的顶部边界, 并以此为构造边界, 划分出北大别弧杂岩及南大别碰撞杂岩等岩石构造单位[35, 36].徐树桐等[37]认为, 大别地块北部的镁铁质和超镁铁质岩带, 是变质的蛇绿混杂岩, 代表消减了的洋壳残余及古缝合线的位置.最近, Hacker等[38]则认为, 中-朝与扬子克拉通间的缝合线, 可能是北淮阳构造带内泥盆纪南湾组的北界接触带.桐柏-大别碰撞造山带的组成和结构的研究表明, 水吼-五河断裂带只是下伸展拆离带的一部分, 而镁铁质-超镁铁质岩带内岩石是由含榴辉岩包体的变形的方辉橄榄岩、纯橄榄岩组合及未变形的辉石岩、角闪辉石岩和辉长岩等年轻侵入体两部分构成的, 地球化学资料及地质体间几何关系、变形行为, 都不具备变质蛇绿混杂岩的典型特征, 因而, 不能看作为三叠纪时期古缝合线.依据超高压和高压变质岩石的空间分布, 构造岩石单位间几何学关系, 以及变质相的构造配置[39]等分析资料, 笔者认为大别地块北缘的八里畈-磨子潭-晓天断裂带以及苏鲁地区的五莲-烟台断裂带, 是被强烈改造了的三叠纪碰撞缝合线.在桐柏山地区, 高压变质岩石分布的北限是娘娘庙剪切带, 它向东与大别地区的八里畈-晓天-磨子潭剪切带相连.位于八里畈-磨子潭-晓天缝合线以北20 km的龟山-梅山断裂带, 是商丹俯冲消减带的东延部分, 主要代表古生代阶段的缝合线.
-
表 1 滑坡物理-力学参数统计值
Table 1. Statistics of landslide physical-mechanical parameters
物理参数 天然密度(kg∙m-3) 饱和密度(kg∙m-3) 2 500 2 600 力学参数 天然内聚力(MPa) 残余内聚力(MPa) 饱和内聚力(MPa) 饱和残余内聚力(MPa) 天然内摩擦系数 残余内摩擦系数 饱和内摩擦系数 饱和残余内摩擦系数 0.020 0.015 0.015 0.010 0.364 0.325 0.325 0.287 变异系数 $ {\delta }_{{\mathrm{c}}_{\mathrm{p}}} $ $ {\delta }_{{\mathrm{c}}_{\mathrm{r}}} $ $ {\delta }_{{\mathrm{c}}_{\mathrm{s}\mathrm{p}}} $ $ {\delta }_{{\mathrm{c}}_{\mathrm{s}\mathrm{r}}} $ $ {\delta }_{{\mathrm{f}}_{\mathrm{p}}} $ $ {\delta }_{{\mathrm{f}}_{\mathrm{r}}} $ $ {\delta }_{{\mathrm{f}}_{\mathrm{s}\mathrm{p}}} $ $ {\delta }_{{\mathrm{f}}_{\mathrm{s}\mathrm{r}}} $ 0.13 0.13 0.18 0.18 0.19 0.19 0.20 0.20 相关系数 $ {\rho }_{{\mathrm{c}}_{\mathrm{p}}, {\mathrm{f}}_{\mathrm{p}}}={\rho }_{{\mathrm{c}}_{\mathrm{r}}, {\mathrm{f}}_{\mathrm{r}}}={\rho }_{{\mathrm{c}}_{\mathrm{s}\mathrm{p}}, {\mathrm{f}}_{\mathrm{s}\mathrm{p}}}={\rho }_{{\mathrm{c}}_{\mathrm{s}\mathrm{t}}, {\mathrm{f}}_{\mathrm{s}\mathrm{r}}}=-0.3 $ 注:变异系数与相关系数是根据现有研究成果进行的合理取值,其实际值应通过统计分析获得. 表 2 滑坡计算工况及滑体条分情况
Table 2. Calculation conditions of landslide and slice details of sliding body
滑坡计算工况 滑体条分情况 145 m水库蓄水 145 m水库蓄水+地震作用 175 m水库蓄水 175 m水库蓄水+地震作用 175 m库水位陡降为145 m库水位 175 m库水位陡降为145 m库水位+地震作用 注:红色实线为地下水位线;库水位陡降时,不考虑坡体内孔隙水压力的消散,即假定坡体内的地下水位线与降落前相同,仅考虑库水位的变动(郑颖人等,2004). 表 3 各工况条件下滑坡峰值稳定性系数及渐进破坏路径
Table 3. Stability coefficient for peak parameters and progressive failure path of landslide under various working conditions
计算工况 峰值强度参数下滑坡稳定性系数 渐进破坏路径 145 m水库蓄水 中前部 中后部 1.217 1.322 145 m水库蓄水+地震作用 中前部 中后部 / 0.840 175 m水库蓄水 中前部 中后部 1.112 1.322 175 m水库蓄水+地震作用 中前部 中后部 / 0.840 175 m库水位陡降为145 m库水位 中前部 中后部 1.044 1.322 175 m库水位陡降为145 m库水位+地震作用 中前部 中后部 / 0.840 注:稳定性系数不唯一且小于0.840. 表 4 滑坡最大概率发生的渐进破坏事件及概率值
Table 4. Progressive failure events with maximum probability and their probability values of landslide
计算工况 渐进破坏事件 发生概率(%) 145 m水库蓄水 整体渐进破坏 2.067×10-13 局部渐进破坏:16#条块发生破坏且破坏停止传递 14.119 145 m水库蓄水+地震作用 整体渐进破坏 34.067 175 m水库蓄水 整体渐进破坏 2.796×10-11 局部渐进破坏:16#条块发生破坏且破坏停止传递 20.266 175 m水库蓄水+地震作用 整体渐进破坏 38.061 175 m库水位陡降为145 m库水位 整体渐进破坏 4.104×10-10 局部渐进破坏:20#条块发生破坏且破坏停止传递 21.797 175 m库水位陡降为145 m库水位+地震作用 整体渐进破坏 38.405 注:计算所得概率值并非真实概率值,但计算结果具备较强的参考性. -
Bishop, A.W., 1967.Progressive Failure with Special Reference to the Mechanism Causing It. In: Proceedings of the Geotechnical Conference on Shear Strength Properties of Natural Soil and Rocks. Norwegian Geotechnical Institute, Oslo, 142-150. Bishop, A. W. , 1971. The Influence of Progressive Failure on the Choice of the Method of Stability Analysis. Géotechnique, 21(2): 168-172. https://doi.org/10.1680/geot.1971.21.2.168 Chen, G. Q. , Huang, R. Q. , Zhou, H. , et al. , 2013. Research on Progressive Failure for Slope Using Dynamic Strength Reduction Method. Rock and Soil Mechanics, 34(4): 1140-1146(in Chinese with English abstract). http://www.cnki.com.cn/Article/CJFDTotal-YTLX201304038.htm Chen, J. H. , Li, J. L. , Xu, X. L. , et al. , 2017. Algorithm for Generation Correlative Variables and Monte Carlo Simulation of Slope Reliability. Rock and Soil Mechanics, 38(11): 3341-3346(in Chinese with English abstract). http://www.researchgate.net/publication/322338618_Algorithm_for_generation_correlative_variables_and_Monte_Carlo_simulation_of_slope_reliability Chen, L. H. , Chen, Z. Y. , Liu, J. M. , 2005. Probability Distribution of Soil Strength. Rock and Soil Mechanics, 26(1): 37-40, 45(in Chinese with English abstract). http://www.researchgate.net/publication/293256313_Probability_distribution_of_soil_strength Chowdhury, R. N. , Tang, W. H. , Sidi, I. , 1987. Reliability Model of Progressive Slope Failure. Géotechnique, 37(4): 467-481. https://doi.org/10.1680/geot.1987.37.4.467 Fan, M. Q. , Sheng, J. B. , 1997. Cross-Correlation of Soil Strength Indicators φ, c. Chinese Journal of Geotechnical Engineering, 19(4): 100-104(in Chinese). http://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CJFD&filename=YTGC704.016 Hu, X. W. , Tang, H. M. , Liu, Y. R. , 2005. Physical Model Studies on Stability of Zhaoshuling Landslide in Area of Three Gorges Reservoir. Chinese Journal of Rock Mechanics and Engineering, 24(12): 2089-2095(in Chinese with English abstract). http://www.oalib.com/paper/1484206 Jiang, S. H. , Yao, C. , Yang, J. H. , et al. , 2018. Model Correction Factor Method Based Approach for Reliability Analysis of Spatially Variable Slopes. Engineering Mechanics, 35(8): 154-161(in Chinese with English abstract). http://www.researchgate.net/publication/328651101_Model_correction_factor_method_based_approach_for_reliability_analysis_of_spatially_variable_slopes Li, D. Q. , Jiang, S. H. , Zhou, C. B. , et al. , 2013. Reliability Analysis of Slopes Considering Spatial Variability of Soil Parameters Using Non-Intrusive Stochastic Finite Element Method. Chinese Journal of Geotechnical Engineering, 35(8): 1413-1422(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-YTGC201308008.htm Li, D. Q. , Xiao, T. , Cao, Z. J. , et al. , 2016a. Slope Risk Assessment Using Efficient Random Finite Element Method. Rock and Soil Mechanics, 37(7): 1994-2003(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-YTLX201607021.htm Li, D. Q. , Xiao, T. , Cao, Z. J. , et al. , 2016b. Auxiliary Slope Reliability Analysis Using Limit Equilibrium Method and Finite Element Method. Chinese Journal of Geotechnical Engineering, 38(6): 1004-1013(in Chinese with English abstract). http://www.researchgate.net/publication/305430876_Auxiliary_slope_reliability_analysis_using_limit_equilibrium_method_and_finite_element_method Li, P. , Bai, J. Z. , Griffiths, D. V. , et al. , 2019. Random Finite Element Analysis for the Reliability of Loess Slopes. Journal of Earth Sciences and Environment, 41(1): 116-126(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-XAGX201901011.htm Li, Q., 2008.Reliability Analysis of Loess High Slope: Take Baojixia Division Project as an Example (Dissertation).Northwest A & F University, Yangling (in Chinese with English abstract). Li, S. L. , Xu, Q. , Tang, M. G. , et al. , 2018. Study on Spatial Distribution and Key Influencing Factors of Landslides in Three Gorges Reservoir Area. Earth Science, 45(1): 341-354(in Chinese with English abstract). Li, Y. C. , Ling, D. S. , Chen, Y. M. , et al. , 2005. Slope Stability Analysis Using Monte Carlo Technique with Fem. Chinese Journal of Rock Mechanics and Engineering, 24(11): 1933-1941(in Chinese with English abstract). http://www.researchgate.net/publication/289702087_Slope_stability_analysis_using_Monte_Carlo_technique_with_FEM Liang, R. Y. , Nusier, O. K. , Malkawi, A. H. , 1999. A Reliability Based Approach for Evaluating the Slope Stability of Embankment Dams. Engineering Geology, 54(3-4): 271-285. https://doi.org/10.1016/s0013-7952(99)00017-4 Lin, S. , Wang, W. , Deng, X. H. , et al. , 2019. Geophysical Observation of Typical Landslides in Three Gorges Reservoir Area and Its Significance: A Case Study of Sifangbei Landslide in Wanzhou District. Earth Science, 44(9): 3135-3146(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-DQKX201909026.htm Potts, D. M. , Zdravković, L. , Addenbrooke, T. I. , et al. , 2001. Finite Element Analysis in Geotechnical Engineering: Theory. Thomas Telford, London. Shen, Z. J. , 2000. Theoretical Soil Mechanical. China Water & Power Press, Beijing (in Chinese). Shi, X. G. , Xu, J. H. , Jiang, H. J. , et al. , 2019. Slope Stability State Monitoring and Updating of the Outang Landslide, Three Gorges Area with Time Series InSAR Analysis. Earth Science, 44(12): 4284-4292(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-DQKX201912038.htm Srbulov, M. M. , 1995. A Simple Method for the Analysis of Stability of Slopes in Brittle Soil. Soils and Foundations, 35(4): 123-127. https://doi.org/10.3208/sandf.35.4_123 Tang, H. M. , Ma, S. Z. , Liu, Y. R. , et al. , 2002. Stability and Control Measures of Zhaoshuling Landslide, Badong County, Three Gorges Reservoir. Earth Science, 27(5): 621-625(in Chinese with English abstract). http://d.old.wanfangdata.com.cn/Periodical/dqkx200205024 Tang, Z. H. , Yu, X. L. , Chai, B. , et al. , 2019. Energetic Criterion of Entering Acceleration in Progressive Failure Process of Bedding Rockslide: A Case Study for Shanshucao Landslide. Earth Science, 1-9(in Chinese with English abstract). https://doi.org/10.3799/dqkx.2019.960 Tu, F. , 2004. Reliability of Progressive Failure of Soil Slope. Rock and Soil Mechanics, 25(1): 87-90(in Chinese with English abstract). http://www.researchgate.net/publication/292877979_Reliability_of_progressive_failure_of_soil_slope Vanmarcke, E. H. , 1977a. Probabilistic Modeling of Soil Profiles. Journal of Geotechnical and Geoenvironmental Engineering, 103(11): 1227-1246. https://doi.org/10.1016/0148-9062(78)90012-8 Vanmarcke, E. H. , 1977b. Reliability of Earth Slopes. Journal of the Geotechnical, 103(11): 1247-1265. https://doi.org/10.1016/0148-9062(78)90159-6 Wang, Y. G. , Ren, W. Z. , Chen, H. , et al. , 2006. First-Order Second-Moment Method and Its Application in Slope Reliability Analysis. Journal of China & Foreign Highway, 26(2): 42-46(in Chinese). Wu, X.M., 2006.Reliability Analysis of Progressive Failure of Homogeneous Earth Slope (Dissertation).Zhejiang University, Hangzhou (in Chinese with English abstract). Xie, X. Y. , Feng, X. , Wu, X. M. , 2015. Reliability Analysis of Progressive Failure of Strain Softening Slope. Rock and Soil Mechanics, 36(S2): 679-684(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-YTLX2015S2099.htm Xu, P., 2011. Reliability Study on Stability of Submerged Slide Mass in Three-Gorges Reservoir Area (Dissertation).Chang'an University, Xi'an (in Chinese with English abstract). Xue, H. B. , Dang, F. N. , Yin, X. T. , et al. , 2016. Stability Analysis Methods for Strain-softening Slopes. Chinese Journal of Geotechnical Engineering, 38(3): 570-576(in Chinese with English abstract). Yang, Q. , Jiao, J. K. , Luan, M. T. , et al. , 2000. Reliability Analysis and Risk Evaluation of the Slope Engineering. Journal of Engineering Geology, 8(1): 86-90(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-GCDZ200001015.htm Yang, Z. Y. , Li, D. Q. , Cao, Z. J. , et al. , 2018. System Reliability of Soil Slope Using Generalized Subset Simulation. Rock and Soil Mechanics, 39(3): 957-966, 984(in Chinese with English abstract). http://www.researchgate.net/publication/327061016_System_reliability_of_soil_slope_using_generalized_subset_simulation Yin, Y. P. , Wang, W. P. , 2014. Researches on Seismic Landslide Stability Analysis. Journal of Engineering Geology, 22(4): 586-600(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-GCDZ201404007.htm Zhang, G. , Zhang, J. M. , 2007. Stability Evaluation of Strain-Softening Slope Based on Swedish Slice Method. Rock and Soil Mechanics, 28(1): 12-16(in Chinese with English abstract). http://www.cqvip.com/QK/94551X/20071/23663872.html Zheng, Y. R. , Shi, W. M. , Kong, W. X. , 2004. Calculation of Seepage Forces and Phreatic Surface under Drawdown Conditions. Chinese Journal of Rock Mechanics and Engineering, 23(18): 3203-3210(in Chinese with English abstract). http://d.wanfangdata.com.cn/Periodical/yslxygcxb200418031 Zhu, B. , Pei, H. F. , Yang, Q. , 2019. Gaussian Process Regression-Based Response Surface Method and Reliability Analysis of Slopes. Chinese Journal of Geotechnical Engineering, 41(S1): 209-212(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-YTGC2019S1054.htm Zhu, Y. X. , 1993. Reliability Analysis of Slope. Metallurgical Industry Press, Beijing (in Chinese). 陈国庆, 黄润秋, 周辉, 等, 2013. 边坡渐进破坏的动态强度折减法研究. 岩土力学, 34(4): 1140-1146. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201304038.htm 陈将宏, 李建林, 许晓亮, 等, 2017. 相关变量生成算法及边坡可靠度Monte Carlo模拟. 岩土力学, 38(11): 3341-3346. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201711034.htm 陈立宏, 陈祖煜, 刘金梅, 2005. 土体抗剪强度指标的概率分布类型研究. 岩土力学, 26(1): 37-40, 45. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200501008.htm 范明桥, 盛金保, 1997. 土强度指标φ, c的互相关性. 岩土工程学报, 19(4): 100-104. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC704.016.htm 胡修文, 唐辉明, 刘佑荣, 2005. 三峡库区赵树岭滑坡稳定性物理模拟试验研究. 岩石力学与工程学报, 24(12): 2089-2095. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200512019.htm 蒋水华, 姚池, 杨建华, 等, 2018. 基于模型修正的空间变异边坡可靠度分析方法. 工程力学, 35(8): 154-161. https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX201808020.htm 李典庆, 蒋水华, 周创兵, 等, 2013. 考虑参数空间变异性的边坡可靠度分析非侵入式随机有限元法. 岩土工程学报, 35(8): 1413-1422. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201308008.htm 李典庆, 肖特, 曹子君, 等, 2016a. 基于高效随机有限元法的边坡风险评估. 岩土力学, 37(7): 1994-2003. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201607021.htm 李典庆, 肖特, 曹子君, 等, 2016b. 基于极限平衡法和有限元法的边坡协同式可靠度分析. 岩土工程学报, 38(6): 1004-1013. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201606005.htm 李萍, 白健忠, Griffiths, D. V. , 等, 2019. 黄土边坡可靠度的随机有限元分析. 地球科学与环境学报, 41(1): 116-126. https://www.cnki.com.cn/Article/CJFDTOTAL-XAGX201901011.htm 李琦, 2008. 黄土高边坡稳定可靠性分析(硕士学位论文). 杨陵: 西北农林科技大学. 李松林, 许强, 汤明高, 等, 2018. 三峡库区滑坡空间发育规律及其关键影响因子. 地球科学, 45(1): 341-354. doi: 10.3799/dqkx.2017.576 李育超, 凌道盛, 陈云敏, 等, 2005. 蒙特卡洛法与有限元相结合分析边坡稳定性. 岩石力学与工程学报, 24(11): 1933-1941. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200511022.htm 林松, 王薇, 邓小虎, 等, 2019. 三峡库区典型滑坡地球物理实测及其意义: 以万州区四方碑滑坡为例. 地球科学, 44(9): 3135-3146. doi: 10.3799/dqkx.2019.074 沈珠江, 2000. 理论土力学. 北京: 中国水利水电出版社. 史绪国, 徐金虎, 蒋厚军, 等, 2019. 时序InSAR技术三峡库区藕塘滑坡稳定性监测与状态更新. 地球科学, 44(12): 4284-4292. doi: 10.3799/dqkx.2019.180 唐辉明, 马淑芝, 刘佑荣, 等, 2002. 三峡工程库区巴东县赵树岭滑坡稳定性与防治对策研究. 地球科学, 27(5): 621-625. http://www.earth-science.net/article/id/1175 唐朝晖, 余小龙, 柴波, 等, 2019. 顺层岩质滑坡渐进破坏进入加速的能量学判据. 地球科学, 1-9. https://doi.org/10.3799/dqkx.2019.960 涂帆, 2004. 土坡渐进破坏的可靠度. 岩土力学, 25(1): 87-90. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200401020.htm 王永刚, 任伟中, 陈浩, 等, 2006. 一次二阶矩法及其在边坡可靠性分析中的应用. 中外公路, 26(2): 42-46. https://www.cnki.com.cn/Article/CJFDTOTAL-GWGL200602012.htm 吴晓明, 2006. 均质土坡渐进破坏可靠性分析(硕士学位论文). 杭州: 浙江大学. 谢新宇, 冯香, 吴晓明, 2015. 应变软化土坡渐进破坏的可靠度分析. 岩土力学, 36(S2): 679-684. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2015S2099.htm 徐平, 2011. 三峡库区涉水滑坡体稳定性的可靠度研究(博士学位论文). 西安: 长安大学. 薛海斌, 党发宁, 尹小涛, 等, 2016. 应变软化边坡稳定性分析方法研究. 岩土工程学报, 38(3): 570-576. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201603026.htm 杨庆, 焦建奎, 栾茂田, 等, 2000. 边坡可靠性与经济风险性分析及其应用. 工程地质学报, 8(1): 86-90. https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ200001015.htm 杨智勇, 李典庆, 曹子君, 等, 2018. 基于广义子集模拟的土坡系统可靠度分析. 岩土力学, 39(3): 957-966, 984. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201803024.htm 殷跃平, 王文沛, 2014. 论滑坡地震力. 工程地质学报, 22(4): 586-600. https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ201404007.htm 张嘎, 张建民, 2007. 基于瑞典条分法的应变软化边坡稳定性评价方法. 岩土力学, 28(1): 12-16. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200701003.htm 郑颖人, 时卫民, 孔位学, 2004. 库水位下降时渗透力及地下水浸润线的计算. 岩石力学与工程学报, 23(18): 3203-3210. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200418038.htm 朱彬, 裴华富, 杨庆, 2019. 基于高斯过程回归的响应面法及边坡可靠度分析. 岩土工程学报, 41(S1): 209-212. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2019S1054.htm 祝玉学, 1993. 边坡可靠性分析. 北京: 冶金工业出版社. -