(4)高应力-强卸荷深部破裂机理
近十年来,在西南地区大型水电工程高边坡的地质勘探过程中,常常揭露边坡具有深部卸荷,并伴随深部张裂的现象,即除了边坡浅表部发育的正常卸荷带(一般0~60m)外,在坡体的深部(一般水平距岸坡120-200m深度)还发育有深部的卸荷带,表现为典型的深部张裂。这种现象最早在雅砻江锦屏水电站普斯罗沟坝址和金沙江向家坝水电站马布坎高边坡被揭露)[15,18,19],并引起重视和研究。类似的深部张裂现象后来在白龙江苗家坝水电站、澜沧江小湾水电站、糯扎渡水电站等又有揭露。
大量的勘探实践和研究资料表明,深部张裂带的发育具有多种复杂的机理,包括本文中所提到的几种模式中实际上也包含了深部张裂带的成分。但作为强裂卸荷原因而引发的这种深部卸荷及其所伴随的深部张裂现象有以下三个强烈的背景:一是高地应力(现实边坡内部的水平或近水平应力通常在15Mpa以上,河谷下切释放前应该更高),这种高地应力是驱动边坡发生强烈回弹变形的内在动力;二是边坡深部存在有利于应力释放的结构面(平行或与边坡小角度斜交的近直立或倾坡外的断层或长大裂隙), 这是深部卸荷和张裂带形成的边坡结构基础;三是河流的快速下切,这是导致边坡内在应力快速释放的外部条件[20]。
典型的高应力-强卸荷深部张裂以雅砻江锦屏水电站普斯罗沟坝址左岸反倾边坡最为典型。主要表现为除浅表部的正常卸荷裂隙外,还发育有一系列深部裂缝,这些深部裂缝的分布特征、成因机制对高边坡和坝肩岩体稳定性有重大影响. 近十年的研究表明:锦屏水电站普斯罗沟左岸高边坡所发育的一套深部裂缝体系,实际上是在坝区特定高地应力环境条件下,伴随河谷快速下切过程中,坡体应力强烈释放,而沿坡体内原有的构造结构面(小断层和长大裂隙)卸荷拉裂的产物。
4 高边坡稳定性评价
在高边坡稳定性评价方面,比较一般的边坡稳稳定性而言,首先是在认识上所取得的进步:
(1)首先,人们认识到了高边坡作为一类特殊的地质体,正确认识其变形破坏机理对稳定性评价的重要意义。岩石高边坡的稳定性不是静止的,而是一个动态演化的地质历史过程,这个过程就是伴随时效变形的发生,边坡潜在滑动面不断的孕育、发展演化,最终进入累进性破坏而贯穿的过程。查明这一过程和所发生的内部作用机理,是岩石高边坡稳定性和适宜性评价的基础[21]。
(2)岩石高边坡稳定性受复杂岩体结构的控制,基于均匀介质的传统土力学稳定性评价方法用于这类边坡的稳定性评价与控制分析是不妥当的。因此,在高边坡稳定性研究过程中,要特别重视对岩体结构特征的研究,重视岩体结构对变形破坏机理模型的控制,重视岩体结构与潜在滑动面关系的分析。
(3)岩石高边坡稳定性的评价不仅是一个强度稳定性问题,也是一个变形稳定性问题,不同的发展演化发展阶段,对应了岩石高边坡所处的不同稳定性状态。在岩石高边坡的稳定性评价中要贯彻变形稳定性评价的思想。
对于特定的边坡地质体,其总存在一个临界破坏的最大变形量,这个变形量代表了边坡最大的承担变形的能力,因此,也是高边坡的安全控制与预警标准制定的基础。随着二维和三维数值模拟技术的发展,人们可以在边坡监测数据的帮助下,通过较为准确的模拟建立边坡变形的数值模型,从而对这个最大变形量做出估计。而监测技术的发展已经可以对边坡的表面和内部变形做出全方位的监测,从而给出不同阶段边坡变形安全度。
(4)岩石高边坡稳定性的控制,关键在于控制变形。从变形破坏演化的历史分析,边坡是通过变形的发展逐渐累计内部的“损伤”,并向潜在滑动面转移,随着“损伤”在潜在滑动面的积累,潜在滑动面逐渐孕育、发展并最终形成边坡的控制性破坏面。因此,在特定阶段,边坡的变形控制住了,不具备进一步发展的条件了,潜在滑动面的演化就会在“孕育”或者“发展”阶段结束,而进入不了最终的累进性破坏阶段。
在具体分析评价方法方面,主要有以下及方面突出的成果:
(1)数值模拟技术在高边坡稳定性评价中得到了广泛深入的应用[21,22,23]。一方面数值模拟技术本身的发展为高边坡变形与强度稳定性分析提供了强有力的工具,三维数值模拟技术得到了广泛的应用。随着计算机技术的发展,以有限单元法为基础的数值模拟技术在90年代中期突破了计算机内存和算法上的限制,使得开展适用于高边坡稳定性分析的大规模科学计算 成为可能。到上个世纪末,在高挡微机上已经可以解决十万个节点以内的大型数值计算和数值模拟问题。目前,随着并行计算技术的发展,正在向数十~百万个节点的计算问题迈进。
另一方面,各种更适合于岩体计算的数值模拟方法自20世纪80年代以来得到极大的发展[24],首先是刚性离散单元法(DEM, Distinct Element Method)的提出,使得节理岩体模拟这种更接近于块体运动的过程模拟成为可能,实现了岩体数值模拟技术一次质的飞跃。紧接着, DDA(Discontinous deformation Analysis, 石根华,Goodman ,1989),NMM (Numerical Manifold Method, 石根华, 1995年)相继提出,它适用于不连续介质的大变形问题分析,从而实现了岩石块体的移动、转动、张开、闭合等全过程的模拟,据此可判断出岩体的破坏程度、破坏范围,从而对岩体的整体和局部的稳定性作出正确的评价。在数值模拟技术方面另一标志性成就是FLAC(Fast Lagrangion Analysis of Continue,P.A.Cundal1,1996)方法的提出。这种方法可以考虑材料的非线性和几何学上的非线性,采用混和离散化法使塑性破坏和塑性流动得到体现;并采用显式时间差分解析法,大大提高了运算速度;适用于求解非线性大变形,但节点的位移连续,本质上仍属于求解连续介质范畴的方法[25]。
对应用数值模拟技术的另一方面的促进是人们在认识岩体变形破坏本质和本构关系描述方面取得了长足的进展。自70~80年代弹塑性、断裂、流变问题提出以来,90年年代关于损伤分析、流变~损伤耦合分析、渗流~损伤耦合分析得到了岩石力学和工程地质学家的重视[26,27,28],并开始应用于诸如三峡船闸高边坡等工程的分析。总之,人们对地质体力学行为的描述与刻画正越来越接近于客观实际,尽管这方面还有很长的路要走。
(2)边坡系统非线性动力学分析。20世纪90年代,随着非线性科学的兴起,这一研究复杂系统的理论在许多领域迅速得到应用,这也促使人们思考边坡系统的非线性行为,并开始了这一新理论在工程地质学领域的应用。这一期间,许多作者在这一方面做了很好的开创性工作[11,12,29,30,31],建立了诸多从描述边坡变形、滑动面发展直到突变失稳的非线性模型。典型的成果包括岩石边坡滑动面累进性破坏过程的自组织临界过程描述,从而从理论上很好的刻画了岩石边坡滑动面的形成和破坏机理,也依此建立了基于自组织临界的滑动面扩展损伤判据和失稳前兆识别判据。依据非线性科学的协同理论、突变理论,建立了边坡失稳预报的协同预测模型和针对不同类型边坡的多种突变理论预报模型。但总的来看,虽然我们在应用非线性科学理论理解边坡的失稳过程和建立预报判据方面取得了显著的进展,但由于客观地质体的复杂性和理论模型的局限性,非线性理论的应用还是一个值得进一步探索的领域。
(3)岩石高边坡的动力响应分析。高边坡的动力响应是高地震烈度区一个必须引起重视的问题。关于这方面最早的研究见于文献[32],他们的后续成果建立了边坡块体运动的动力微分方程,讨论了动力稳定的评价原则[33,34]。文献[35]通过动力有限元发现岩石边坡的地震动力系数并不随坡高增高而单调增大,当坡高约100m时,坡顶动力放大系数达到最大值;坡高超过100m动力系数反而有所降低。1991年长江科学院采用有限单元法研究了三峡船闸高边坡的地震动力稳定性。但总的说来,20世纪国内对边坡,特别是复杂的岩石高边坡动力响应问题的分析理论和方法上仍显不足,是今后有待发展的重点领域。
(4)边坡安全可靠度分析。从20世纪90年代开始,针对边坡稳定性评价中诸多不确定性因素的问题,许多作者开展了边坡安全可靠度分析的研究[36,37],并取得了显著的进展,建立了特定条件下的影响因素的概率模型和边坡可靠度分析模型。但由于地质条件的复杂性,目前边坡可靠度分析还仅局限于针对一些条件相对简单的边坡的分析和应用。