三维地质模型支持下的滑坡体稳定性分析.pdf
《三维地质模型支持下的滑坡体稳定性分析.pdf》由会员分享,可在线阅读,更多相关《三维地质模型支持下的滑坡体稳定性分析.pdf(11页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第1 5 卷1 期2 0 0 7 年3 月应用基础与工程科学学报J O U R N A LO FB A S I CS C I E N C EA N DE N G I N E E R I N GV 0 1 1 5 N o 1M a r c h 2 0 0 7文章编号:1 0 0 5-0 9 3 0(2 0 0 7)0 1-0 0 6 5-0 9中图分类号:1 7 1 4 5 7;P 6 4 2文献标识码:A三维地质模型支持下的滑坡体稳定性分析钟登华1,秦朝霞1,李明超1,胡兴娥2(1 天津大学建筑工程学院,天津3 0 0 0 7 2;2 中国长江三峡工程开发总公司,湖北宜昌4 4 3 0 0 2
2、)摘要:针对具有复杂几何形状和地质条件的滑坡稳定性分析问题,提出实现基于滑坡地质体三维建模的稳定性三维分析方法与系统根据地质勘探和二维解译资料,运用N U R B S B R e p 混合数据结构建模方法快速、准确地重构滑坡体三维地质模型该模型不仅包含直观丰富的地质信息和反映空间变化的复杂滑动面,完成了滑坡的可视化表达;而且为滑坡稳定性计算提供滑坡体三维剖分几何模型和大量的基础数据,极大地方便了三维稳定计算和敏感性分析的实现应用表明,该方法实现的系统功能强大,能获得大量的计算分析成果,满足实际分析需要,提高了滑坡稳定性三维分析的数字化与可视化水平关键词:滑坡稳定性;三维地质模型;三维分析;敏感
3、性分析滑坡研究的目的是认识滑坡的发育机理和运动破坏方式,进而确定滑坡的防治措施,其可分为三个方面:滑坡表达、滑坡稳定性分析、滑坡工程治理,目前在工程实践中大都仍将其作为二维问题处理,其中滑坡表达是滑坡研究的基础,目的在于描述滑坡体的空间分布特征、滑动面的位置和形态、滑坡体的物质组成和地质结构特征,为正确表达滑坡的时空分布规律和内部结构特征,需进行滑坡体三维地质建模研究滑坡稳定性分析一直就是人们研究的难点和焦点,目前人们广泛接受并实践应用的仍是极限平衡分析方法,而且更多侧重于三维研究,如冯树仁等H 1 提出了简单实用的三维条分法,H u a n g 等悼1 采用两个方向的力和力矩平衡推导出了通用
4、的三维稳定分析方法,C h e n 等一。则将二维S p e n c e r法进行了三维扩展,X i e 等H 1 结合G I S 空间分析功能进行滑坡三维稳定性分析,姜清辉等障1 有机地结合了三维边坡地层信息和极限平衡分析方法研制开发了S L O P E 3 D 系统,等上述研究成果推动了滑坡稳定性三维分析理论方法与应用的发展,但大多数程序或系统在三维数据处理、图形分析等方面难以提供便于实际工程应用的条件1,需要进一步的深入研究本文主要研究滑坡表达和稳定性分析两部分内容,通过滑坡地质体三维建模不仅实现滑坡的可视化表达,而且为滑坡三维稳定分析提供基础几何模型,进而采用三维极限平衡法快速、准确地
5、计算滑坡稳定系数,给出综合的评价结果因此,笔者研制开发了滑坡地质体三维可视化与稳定分析系统V i s u a l S l o p e 来实现上述内容收稿日期:2 0 0 6-0 3-2 0;修订日期:2 0 0 6 1 0-2 3基金项目:国家杰出青年科学基金(5 0 5 2 5 9 2 7,5 0 5 3 9 1 2 0)和国家自然科学基金(5 0 4 7 9 0 4 8,5 0 5 3 9 1 2 0)资助项目作者简介:钟登华(1 9 6 3 一),男,博士,教授 万方数据应用基础与工程科学学报1V i s u a l S l o p e 系统总体结构V i s u a l S l o p
6、 e 系统总体的模块框架结构如图1 所示,由滑坡地质体三维建模、滑坡稳定总控界面滑坡地形模型滑坡地质体三维建模滑坡稳定性分析滑坡数据库管理图1V i s u a l S l o p e 系统总体结构图F i g 1O v e r a l ls t r u c t u r eo fV i s u a I S l o p es y s t e m段整个系统结构规范且具有良好的维护性2 滑坡地质体三维建模性分析和数据库管理三个模块组成三个模块不仅分工明确,易于程序实现,而且相互关联紧密滑坡地质体三维建模模块能够提供直观准确的三维滑坡地质体模型和反映空间变化的复杂滑动面,并提供滑体三维剖分功能;滑坡稳
7、定性分析模块是核心,在滑坡地质体三维建模与剖分的基础上根据实际需要对滑坡稳定性进行三维稳定计算和敏感性分析等多方面的综合分析评价;滑坡数据库管理模块为滑坡地质体三维建模、稳定性分析与表达等提供对象属性管理和各种输入输出手通过地质勘探获得足够的地质信息,是分析滑坡稳定性的基本前提;同时,对这些地质信息的认识和理解程度是能否客观准确地进行滑坡稳定性计算分析的关键实践证明,任何计算方法的成功都必须建立在深入查明滑坡地质特征和作出符合实际情况的机制分析的基础上 1 因此,利用地质勘探和实验分析所得到的一系列空间分布不均的离散数据来进行三维地质建模与可视化分析,描述地质构造的空间分布情况,并实现人机可视
8、化的交互分析与操作,为正确判断、分析滑坡体工程地质问题提供可视化的综合信息,是非常必要的三维地质建模是目前国内外研究具有挑战性的一个热点问题,在分析这方面研究成果的基础上,针对中国工程地质勘探的特点和实际需求情况,笔者提出以非均匀有理B样条(N o n-U n i f o r mR a t i o n a lB S p l i n e s,N U R B S)结构为主、结合边界表示(B o u n d a r yR e p r e s e n t a t i o n,B R e p)结构的混合数据,进行三维地质结构体建模;其中N U R B S 结构用来拟合构造复杂的地质曲面,B R e p
9、结构则组织N U R B S 地质曲面的拓扑关系,构造相应的地质体地质结构面实际上是利用N U R B S 技术坤1 由各层面上的离散点和线数据构建而成的,而地质结构体又是基于不同结构面而组成的一J 若研究区域为n,则基于混合结构的整体地质模型定义为:M n2U M c im iM c i2S i luS 忍u(k U _ lS l 斌),i=1,2,几(1)S d=5(P d),i=1,2,I I,;J=1,2S l 诜=s(膻),钉矾O S i lu0 S i 2,i=1,2,n;南=1,2,m i一一|一一剖曼 万方数据N o 1钟登华等:三维地质模型支持下的滑坡体稳定性分析6 7式中:
10、M 仃表示研究区域力的整体地质模型;,为力中包含的地质结构实体总数;M“表示n 中的第i 个地质结构实体模型;s s。:分别是构成丝;的上、下(或左、右)主结构面,它们由其层面上的点集P P 也通过N U R B S 技术拟合构造而成;O S i 是结构曲面s。上所有边界顶点的集合;S l 谴则表示连接s n、s 记形成闭合实体M。的第k 个边界面,它是由边界顶点集合 秽捕 形成的简单N U R B S 曲面;m i 为第i 个单元中连接边界面的数目基于式(1),采用面向对象的地质分类思想、地质曲面的拟合技术、改进的地质趋势面分析技术和三维几何对象的任意布尔算法等关键技术0|,可以完整、快速、
11、客观地重构滑坡地质体三维模型,并完成对滑坡体的三维自动剖分3 滑坡稳定三维极限平衡分析法3 1 计算公式本文采用的三维极限平衡法仍然继承了二维条分法的思想,根据力的平衡条件进行分析,主要基于以下三点:(1)分别采用两组正交的垂直平面(一组沿滑动方向,另一组垂直于滑动方向)将滑坡体分为若干条柱体;(2)分块柱体界面问作用力仅考虑水平向的侧面压力和剪力,不考虑铅垂向的剪力;(3)根据铅垂向力的平衡条件和底滑面沿滑动方向力的平衡条件,进行公式推导图2 为一滑体条柱受力分析示意图,其选取的坐标系统使O x y 在水平面内,戈轴与滑动方向一致,并与y、名轴构成右手系,f f 轴垂直向上;条柱滑面与O x
12、 z 平面的交线为图2 三维滑体条柱受力分析F i g 2F o r c e sa c t i n go na3 Ds l o p ec o l u m n仇,与O y z 平面的交线为O y7;角_ x O x=似,Z _ y O y=掣,滑面法向与铅垂方向夹角为y 幽条块是由平行于O x z 平面和平行于o y z 平面的两组等间距平面相交组成,条块顶面为地表面(叫7 曰c),底面为滑带面(O A B C);侧面为铅垂面假设条柱编号为(i,J),作用其上的力主要包括:(1)侧面(忽略垂直方向剪力):作用在A7 A B B7 面、0 O C C 面和B7 B C C7 面、A7 A 0 07
13、 面上的法向力只。,P 扪P。,P 小作用在A A B B 面、0 O C C 面和B B C C7 面、A A 0 0 面上的水平方向剪力咒。,t。,。,t。;(2)底滑面:抗滑力S i;作用在底面的水压力U i;作用在底面的反力N o;(3)地表面:表面水压力(垂直于坡面)F。可分解为:O x 方向的投影F 埘(与滑动方向相反时为正)和竖直方向的投影F y i j(向下为正);(4)其它受力:条块重力形,;作用在A A B B 面和0 O C C 面上的水压力代数和玑;,(沿石轴方向为正);水平方向地震力H i i(沿石轴方向为正);记Q;f=u x;+峨为总的水平外力根据铅垂方向力的平衡
14、条件和底滑面沿滑动方向(仇方向)力的平衡条件,三维稳 万方数据应用基础与工程科学学报定系数F s 取“材料强度储备系数”意义,推导得出其计算公式为:几2 荟善 C o A u c o s y 幽“+(+)一U u c o s y 幽n 珥 S e _ C O 万l*i j 五ic o s y 幽i 十万一 Q F F;F+(+)t a n a。d(2)式中:c i 和A#分别为条柱(i,)的粘聚力、摩擦系数和底面面积数据预处理输入读取数据0I 臁警mI 雨卢o量一i=i+I计算中间变量t 器N 上计算F s*t1:N输出结果西图3 三维稳定系数计算流程图F i g 3F l o wc h a
15、r to fc o m p u t i n g3 Ds t a b i l i t yc o e f f i c i e n t三维稳定系数风值通过试算求解,由于该公式未考虑力矩平衡,所以适用于任意形状滑面3 2 程序实现根据式(2)求解风值的计算工作量大,尤其是剖分条柱数较多时,因此应通过计算机编程实现其流程如图3 所示,主要步骤为:(1)数据预处理处理所建立的滑坡三维模型几何数据要输入规范,包括滑块编号、底滑面数据从几何模型中提取每一滑块的底滑面中心点高程、面积、角度和y:如)、地表面数据(每一滑块的顶面中心点高程、角度7。)、地下水位面中心点高程、体积数据(每一滑块的总体积、地下水位以上
16、以下的体积);(2)数据输入分为基本参数(包括滑坡基本信息、范围、剖分信息等)、力学参数(包括不同介质的粘聚力、摩擦角、天然饱和容重等)、受力参数(包括水、地震力的参数信息)和几何参数(每一条柱的几何信息)四大部分进行数据输入读取,其中几何参数直接读取(1)步处理好的数据获得;(3)计算F 流程图中以表示剖分滑块总数,所计算的中间变量包括每一滑块的重力、粘聚力、底面水压力、水平向水压力、表面水压力和水平地震力等;(4)结果输出通过迭代计算,所求解满足了一定的精度要求(如1 0。4),即可输出F 值4 实例研究某崩塌堆积滑坡体位于水库库区一处斜坡地段,总体地形特征为前缘较陡,中部较缓,后部比中部
17、稍陡滑坡体的底面形态复杂,厚度大,一般厚6 0 1 0 0 m,最厚达1 3 0 m;滑坡体前缘宽约7 0 0 m,滑床出露线最低高程2 5 m 左右,位于河床中心附近;后缘高程4 3 0 一 万方数据钟登华等:三维地质模型支持下的滑坡体稳定性分析6 94 5 0 m,宽1 7 0 m;纵向长1 1 0 0 m;滑坡体前缘宽,后缘窄,总体呈一不规则斧形4 1滑坡体三维地质模型据实地调绘及钻孔揭露,该滑坡体所处区域地层结构较为简单,主要由3 组基岩和第四系覆盖层构成区内出露基岩主要为震旦系上统灯影组(z)(分布于滑坡下游山坡一带)、寒武系下统水井沱组(。正)(分布于滑坡后缘山坡及滑坡下伏基岩)和
18、石牌组()地层(分布于滑坡西侧陡崖的下部)根据滑坡所处范围、现有的地质资料及分析需要,所选择研究区域为:沿河流方向长约1 0 0 0 0 I n,垂直河床方向宽约1 7 0 0 0 m,底一高程为一5 0 O m,最高高程约1 1 9 0 O m 基于该区域内的地形等高线测绘、钻孔平硐勘探资料、滑面等高线和二维解译剖面等基础资图4 滑坡体三维地质模型料,采用第2 节的方法建立该滑坡体的三维n g 43 Dg e o l o g i e a lm o d e lo ft h es l o p e地质模型如图4 所示该模型不仅可以全方位、动态地显示(旋转、平移、放大、缩小等),而且能进行任意剖切,
19、如图5 为沿滑坡主滑动方向的三维剖切图;为滑坡分析提供直观、准确的地质信息4 2 滑坡三维稳定系数计算(1)稳定系数评价标准滑坡稳定系数是判断滑坡稳定性、确定图5 沿主滑动方向的三维地质剖切图滑坡处理工程量的一项重要指标,直接关系F i g 53 Dg e o l o g i。a lp r o f i l ea l o n gt h em a i n着工程的安全性、经济性与合理性然而,目s l i d i n gd i r e c t i。n前滑坡稳定系数的标准并无规范可循,需要根据滑坡的自身情况,参考其它滑坡实例,确定滑坡体三维稳定系数n 评价标准为:F s 1 1 5 为稳定;F sE(1
20、 1 0,1 1 5)为基本稳定;几(1 0 5 1 1 0)为临界稳定;F s 1 0 5 为不稳;(2)计算参数该滑坡为典型的块石堆积体,选取图5 剖面方向为主滑动方向,平硐揭露到滑坡底部滑带土物质,是一个典型的滑动面,其边界条件主要考虑整体失稳模式、相关的计算参数为:地下水:根据观测资料,以钻孔常水位为滑坡体地下水位,受降水影响变化不大;库水位:主要取初期蓄水位(1 3 5 m)、防洪限制水位(1 4 5 m)和正常蓄水位(1 7 5 m);强度参数C、妒值:综合考虑反演分析成果和试验成果,选用粘聚力C=2 2 5k P a,摩擦角妒=2 2 0。;根据现场试验,滑坡体总体天然容重为2
21、2 7 t m 3,饱和容重为2 3 9 t m 3;万方数据7 0应用基础与工程科学学报V 0 1 1 5图6 滑坡体三维剖分模型图F i g 63 Dd i v i d e dm o d e lo ft h es l o p e(3)几何模型基于三维地质模型的滑坡体三维剖分模型图如图6 所示,沿主滑动方向最大长度为1 1 9 6 0 7 m,垂直方向最大宽度为5 6 7 2 3 m;沿滑动方向和垂直滑动方向进行滑坡体剖分的边长均为2 0 0 m,剖分滑块总数1 1 6 7 块;每个滑块的几何参数包括滑块编号、总体积、地下水以上体积、地下水以下体积、滑面面积、O t。7:幽、滑面前后缘高程、
22、表面前后缘高程及地下水前后缘高程等共1 3 组数据信息基于上述条件,利用三维极限平衡法计算3 个典型库水位的正常情况下的稳定系数分别为1 1 5 2、1 1 5 6 和1 1 7 7;均大于1 1 5,因此,判断在此种情况下滑坡体是稳定的4 3滑坡稳定系数敏感性分析影响滑坡体稳定系数的因素有多种,对于该滑坡可概括为:底滑面强度参数C、妒值;水库水位变化;降雨作用;地震作用;后缘加载;库水对松散层的改造作用等6 个方面的影响因素根据滑坡体的地质背景及周围环境,其后缘及侧缘陡崖岩体滑坡已趋于稳定,滑坡体基本无加载量,后缘加载因素可不用考虑;而库水对松散层的改造作用机理复杂,难以定量化,暂不予分析因
23、此主要针对前四个影响因素对滑坡稳定系数的影响度进行敏感性分析(1)强度参数C、妒值凡与滑面强度参数C、妒值关系最为密切,对其进行敏感性分析很有必要这里仅讨论水库正常蓄水位1 7 5 0 m 时滑面c、妒值对凡的影响:C 值分析范围为1 5 4 0 k P a,驴值分析值范围为1 5 _ 3 0。经计算,凡值与C、妒值的关系为:妒值每增加2。则凡值增大约0 1 1;C 值每增加5 K P a 则凡值增大约0 0 2 可知稳定系数对妒值的敏感度要远远大于对c 值的敏感度,而且随着妒值的增大,稳定系数的增大程度越大,因此妒取值的准确性是分析滑坡体稳定性的关键;(2)水库水位变化主要考虑库水位从1 7
24、 5 0 m 下降至1 4 5 0 m 的条件下对滑坡稳定性的影响,下降时问为5 7 d、l O d、1 5 d 和2 0 d 等典型情况实践表明,库区水位骤降是促使库岸滑坡发生失稳的一个重要原因,主要是由于骤降时滑坡体内地下水位下降相对滞后,浸泡滑坡前缘的库水在滑坡内将产生不利于其稳定的动水压力作用这是一个不定常的渗流问题,由于边界条件及渗流介质的复杂性,要精确确定其大小比较困难这里采用的方法为:综合考虑与问题相关的诸多因素,如库水的下降速度、滑体的渗透系数和孔隙率等,根据布西涅斯克(B o u s s i n e s g)非稳定渗流基本微分方程和边界条件,得到库水位等速下降时滑坡体内浸润线
25、的简化计算公式1 I,来进行库水位骤降的敏感性分析,计算结果如表l 所示 万方数据N o 1钟登华等:三维地质模型支持下的滑坡体稳定性分析71表1 库水位下降时的敏感性分析T a b l e1S e n s i t i v i t ya n a l y s i su n d e rt h ed r a w d o w nc o n d i t i o n由结果可知,在5 _ 7 d 内从1 7 5 0 m 骤降至1 4 5 0 m,滑坡稳定系数均在临界范围附近,而随着下降时间延长至1 0 2 0 d 内,其稳定性有较大提高,处于基本稳定范围内;(3)降雨作用该滑坡体区地处多降雨区,多年平均降雨
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 三维 地质模型 支持 滑坡 稳定性 分析
限制150内