基于matlab的边坡虹吸排水渗流场分析.pdf

第 3 9卷 ， 第 5期 2 0 1 4年 1 0月 公 路 工 程 Hi g h wa y En g i n e e r i n g Vo 1 3 9，No 5 Oc t ， 2 0 1 4 基 于 ma t l a b的边坡虹 吸排水渗流场分析 张世 华 ，孙红 月 ，熊 晓亮 ，张文君 ( 浙江大学 港 口海岸及近海工程研究所 ，浙江 杭州3 1 0 0 5 8 ) 摘要 】地表 降雨人渗 会引起边坡地下水位的上升 ， 进 而降低边坡的稳定性 ， 因此降低地下水 位对增加 边坡 的稳 定性极其重要 。虹吸排水是一种免动 力并且 能排除坡体 深部 地下水的措施 ， 利用 ma t l a b优越 的求 解能力对 虹 吸排水后边坡 内的地下水浸润 曲线分 布进行研 究来 寻求 相对 合理有效的虹吸排水孔位置 ， 研究过程 中对虹 吸排水 采用近似的点汇来处理 ， 比较接 近实际的工程。研究的结果表 明浅部排水孔适合边坡应 急而深部排 水孔可 以有效 提高边坡长期稳定性和缓 冲下一次 降雨 引起的边坡地下水位抬升 ， 可 为实际工程提供参考 。 【 关键词 】虹吸排水 ； m a t l a b ； 浸润 曲线；点汇 中图分类号 U 4 1 6 1 4 文献标识码 A 文章编号 】1 6 7 4 - 0 6 1 0 ( 2 0 1 4 ) O 5 0 0 2 5 - 0 6 An a l y s i s o f S l o p e S e e pa g e u n d e r S i p h o n Dr a i n a g e Ba s e d o n M a t l a b M e t h o d Z HANG S h i h u a ，S U N Ho n g y u e ，XI ONG Xi a o l i a n g ，Z HANG W e n j 1 1 9 ( I n s t i t u t e o f P o r t ，C o a s t a l a n d O f f s h o r e E n g i n e e r i n g ， Z h e j i a n g U n i v e r s i t y ， Ha n g z h o u ， Z h e j i a n g 3 1 0 0 5 8 ， C h i n a ) A b s t r a c t I n f i l t r a t i o n o f r a i n f a l l w i l l c a u s e t h e r i s e o f g r o u n d w a t e r l e v e l s i n s l o p e ， w h i c h c a n r e - d u c e t h e s l o pe s t a b i l i t y， t hu s r e du c i n g t h e g r o u n d wa t e r l e v e l f o r i n c r e a s i n g s l o p e s t a b i l i t y i s e x t r e me l y i m- p o rta n t S i p h o n d r a i n a g e i s a f r e e p o we r me a s u r e s wh i c h c a n g e t r i d o f d e e p g r o u nd wa t e r i n s l o p e， u s i n g ma t l a b s u p e r i o r a b i l i t y f o r s o l v i n g t h e d i s t r i b u t i o n o f i n fil t r a t i o n c u r v e a f t e r s i p h o n d r a i na g e wh a t i s h e l pf u l f o r s e e k i n g a r e a s o na b l e a nd e f f e c t i v e s i p h o n d r a i n a g e l o c a t i o n， t h e r e s e a r c h pr o c e s s f o r s i p ho n d r a i n a g e u s i n g a p p r o x i ma t e p o i n t s i n k t o h a n d l e ， r e l a t i v e l y c l o s e t o t h e a c t u a l p r o j e c t R e s u l t s o f t h e s t u d y s h o w s t h a t t h e s h a l l o w d r a i na g e o f t h e s l o p e i s s u i t a b l e for e me r g e n c y a n d d e e p d r a i na g e o f t he s l o p e c a n e f f e c t i v e l y i mp r o v e t he l o ng - t e r m s t a b i l i t y a n d d e f e r t h e e l e v a t i o n o f wa t e r t a b l e wh e n ne x t r a i n f a l l c o me s ， w h i c h c a n p r o v i d e a r e f e r e n c e for t h e a c t u a l p r o j e c t Ke y w o r d s s i p h o n i c d r a i n a g e ； ma t l a b ；I n f i l t r a t i o n c u r v e ； P o i n t s i n k 0 引言 在边坡 变 形破 坏 的影 响 因素 中地 下水 活 动是 大 量滑坡发生的直接原因也是最复杂和变化最频繁的 影响 因素。地下水通过物理的、 化学 的和力学的 等作用改变边坡的结构 ， 从而给坡体施加静水压力 和动水压力 ， 因此排除地下水 、 降低地下水位对 边坡 的稳定性提高具 有显著 的作 用 。当前 边坡 的主要排水措施有： 地下排水洞、 水平排水孔 、 集 水井抽水、 地表排水 沟和排水隧洞 ， 这些 排水措 施存在的主要问题有 ： 排水措施对坡体 的排水环 境要求较高 ， 许多边坡往往缺乏有利的地貌条件 ， 导 致排水措施的有效性难 以保证 ；地表排水和浅部 排水对提高边坡稳定性的效果差 ；地下排水洞费 用较高 ， 施工工艺复杂 ；集水井抽水需要 动力 和 经常性管理。 虹吸排水是一种免动力 实时排水措施， 能克服 当前一些排水措施 的局限且 可 以实 现边坡深部排 水 ， 其在工程 中也 有一定 的应用。如张永 防， 张朝 林 曾就湘黔线 K 9 3路堑 滑坡 虹吸排水工点进行 了试验研究。陈野鹰 ， 唐红梅 在 山峡水库港 口岸 坡排水方案优化研究中采用 了多级串联虹吸排水方 收稿 日期】2 0 1 3 0 9 1 0 【 基金项 目】“ 十二五” 国家 科技支撑计划( 1 2 2 O B A Kl 0 B O 6 )； 国家 自然科学基 金资助项 目( 4 1 2 7 2 3 3 6 ) 作者简介】张世华 ( 1 9 8 7 一) ， 男 ， 江西 上饶人 ， 硕 士研究生 ， 主要从事岩土工程灾害防治技术研究。 学兔兔 w w w .x u e t u t u .c o m 2 6 公路工程 3 9 卷 案。这些研究和工程实例都表明虹吸排水具有传统 排水方案不具备的优点 ， 能在实际工程 中具有较大 的应用前景。因此对虹吸排水后地下水位分布的研 究就显得十分重要和迫切。 地下水流动的控制方程为偏微分方程 ， 而 ma t - l a b是求解偏微分方程的有力工具， 近些年来 m a t l a b 被一些学者广泛的应用 的求解热传 导方程 中， 如李 灿 ， 高彦栋等 】 。由于热传导方程与地下水流动方 程就要 良好的相似性 ， 故可 以采用 m a t l a b求解虹吸 排水过程 中地下水 的分布情况， 为虹吸排水效果的 研 究 奠定 必不 可少 的理论 基础 。 1基本假设及理论 1 1基本 假 设 Ma t l a b具有强大 的数值计算能力 ， 已被应用于 各个领域工程 实践中。一些学 者将 m a t l a b强大的 计算能力应用于热传导方程求解 。鉴于地下水 渗流方程与热传导方程具有 良好 的相似性 ， 因此考 虑采用 m a t l a b对边坡虹吸排水后 的渗流场进 行分 析。由于渗流控制方程建立和求解需要一些理想的 条件 ， 考虑实际工程的复杂性 ， 做了如下的假设 ： 边坡地下水流符合达西定律。 初始地下水的分布具有潜水分布的特征 ， 由 裘布依公式求解得出。 虹吸排水管的出水假设为一个点汇。 在虹吸排水后考虑 由最高地下水位一侧与 边坡长度形成的矩形 区域 内的地下水头分布 ， 高于 地下水浸润曲线的部分假设其压力水头为负值 。 1 2源 汇 理 论 某一吸收一定强度物理量 的无穷小点称之为汇 点 ， 如果放 出一定强度的物理量 ， 则称之为源点。把 放 出一定强度的物理量看作吸收一定强度的负值物 理量 ， 也可统称之为汇点 。当平面 三有一汇点时， 平 面周围的“ 势” 将随距离汇点的远近 ， 时间的长短而 发生变化。汇点本身吸收物理量的强度 ， 它可 以是 时间的函数， 也直接影响周围“ 势” 的变化 。 潜水抽水井对潜水位的动态变化的影响可以应 用汇点来计算 ， 在平面上任取一个坐标系 ， 令水井位 于这个坐标系的( 。 ， Y 。 ) 点上 ， 见图 1 。水井的抽水 量即为汇点的流量 q 。则通 常用于地下水流动方程 中的汇项 可 以表示 为 ： 图 1 汇点位置示意 图 Fi g u r e 1 Lo c a t i o n o f p o i n t s i n k 2初始浸润 曲线分布 边坡 的地下水位会 随着降雨作用逐渐抬升 ， 本 文研究的初始浸润曲线的两侧定水头为降雨作用下 的最危险地下水位( 见 图2 ) 。在未进行虹吸排水时 仅考虑一维流动 ， 近似满足裘布依假设 。 图 2 初始浸润 曲线示意 图 F i g u r e 2 Di a g r am o f I ni t i a l i n f i l t r a t i o n c u r v e 其数学模型 如下 ： 0 ( K h O hO x ) = 0 J h ； 。=H。 h l ： = H 解上 面 的方 程可 得 ： ： +堕 ( 2 ) 式 中： 日。 为坡 顶的地下水位。 为坡脚 的地 下水 位 。l 为边 坡 的长度 。 3 虹 吸排水后的浸润线分布 3 1 虹 吸排 水 定解 问题 的建 立 计算 中将边坡岩土体视 为多孔介质 ， 边坡 中地 下水的流动满足渗流连续方程 。非均质各向异性 的 岩土体 中含有源项的非稳定渗流的计 算公式 如 下 所示 ， 未 ah + K x O h，+ K xzO h ) + 学兔兔 w w w .x u e t u t u .c o m 第 5期 张世华 ， 等 ： 基 于 m a t l a b的边坡 虹吸排水渗流场分析 2 7 ( 芸 + K aa 歹h + K Oh ) + 0 ( K z O h + K O h + K O h ) + W = X d t ( 3 ) 式中： 为渗流场 中各点的水头值 。 K = K K” Kz K K K K K 为岩土体单元在 整体坐 标系下渗透系数张量。 为单元内的源汇项的值 ， 即该项考虑虹吸排水的排水量。 为重力给水度 。 考虑均质各 向同性 的岩土体 中的渗流连 续方 程 “ 为： K 0 2 h O y 2 a h ( 4 ) 令 = 则上式可以改写为 ： d 口 1 0 2 h 0 2h 】 埘 1 = Oh (5 ) 故虹吸排水后的渗流控制方程可以采用上式表 示 ， 则虹吸排水的定解问题可叙述为 ： 口 O y 2 =警 h l 一 、= h L0 l = t0 h l 一 、= h 2l t 0 I = 一O h l ： 0 一I = U a I y ： o 1 0 l ： 0 一I = U d Y I y ： I o ( ) ， ( ， y ) 其 中I ， ( ， y ) 为水头的初始分布， 由于初始仅考 虑一 维流动 ， 故此处的厂 ( ， Y ) 仅为 的函数。 3 2虹 吸排 水控 制 方程 的差分格 式 由于虹吸排水情况下， 边坡潜水层中的地下水 流为二维非稳定流， 即随着虹吸排水 时间的变化而 变化 ， 故求解此定解 问题采用的差分格式 为交替 隐 格式 ( A D I ) 。 其基本思路为： 分别对 和 Y方向采用不 同的 时间层进行处理。 其处理的过程如下所述 ： 令 ： r ， 口 当 l ， 方向为隐式时， 原偏微分方程可化为如 下的差分方程： 1 4 - 1 J 1十 I ， ，J ，一 h 2 一 2 h + h “ -I ， | l ， j I 一1 ， J y 一 一 ， 一 ! ： ! ! 一 f d ( 6 ) 移项可得 ： 一 r y ( k 一+ 。l J + ki + 。l ) + ( 1 + 2 r ) IIl ： = r ( 。 + 。 ) + ( ) ， W 0 0 0 一 y 1 +2r y ： + + F x ( u 。 + u ) + ( 1 2 ) “ ， r ( “ J 一 。 + n J + 。 ) + ( 一 2 r ) “ r ( “ ： ， + ： J + 。 ) + ( 1 2 r ) u ： ， ! k + l ( M y , + “ 。 ) + r y u 】If y + 1 ， ， + r I 卜 l + “ |】If f + l J + ( 1 2 ) ， 由于上式 只有 3个未知数 - f Ih k 一 + 。 lJ， hk ，； ， ： ， 故可用追赶法解此矩阵方程。 当 X方向为隐式时， 原偏微分方程可化为如 下的差 分 方程 ： h 。一2 + I ， ， +l ，J l ， J一1 戈 一2h +h I4 - 1， J ， J 一1 y ，； 一 I ， ， t 肛d 移项可得 ： + ( 7 ) 一 r r -i , j - 。+ 。 ： 。 ) + ( I + 2 r ) = r ( + k + lJ ) + ( - ) W y o 吖 0 y r y r 2 r ) 吖 吖 0 y 2 r ) ) 0 O 学兔兔 w w w .x u e t u t u .c o m 2 8 公路工程 3 9卷 1 +2r 一 r 0 0 一 r 0 1 +2r 一r 一 r 1 +2r 0 0 0 0 0 一 r 1 +2 r r “ ： + r y ( “ k 一+ 。l ， + “ k i+ 。 l ，：) + ( 1 2 r ) “ ： r y ( u k 一+ 。l ， ， + u k + + 。 l ，) ( 1Y 4 - 一 2 r ) u 3 I u f l ，3 + I + 1 ， 3 J I z r y J u r y ( u k 一+ 。l ， + n k + + l ， ) + ( 1 2Y 1 ， ， ) “ ： I u f l ， 4 + f + 1 ，4 J + 【 一， y J 4 ； 4 - 2 ，k 4 - 1 k+1 、 ， i +l+ r ) ， 【 l ， M + + 1 。 J + ( 1 2 r ) “ 4-2 i ,2 4-2 i , 3 +2 4 ： +2 u i M 由于上式只有 3个未知数 h k 。，； 。 ， h k ，； ， h k 4- +2 。 ， 故可用追赶法解此矩阵方程 。 最后由以上求得的各点的水 头分布 ， 通过寻找 到水头值与 Y 值 相等的点作为浸润曲线上 的点 ， 从 而 得到 虹吸 排水后 地 下水浸 润 曲线 的分 布。 4案例分析 该边坡为一均质含砂黏性土坡 ， 坡 高 2 5 m， 黏 聚力 C= 2 0 k P a ， 土体容重 =1 8 k N m ， 内摩擦角 = 2 5 。 ， 重力 给水度取 为 0 0 5 ， 饱和渗透系数取为 1 2 m d 。边坡两侧边为定水头 ， 地下水位分别为 3 5 n l 和 1 3 n l 。坡体中的初始地下水位为降雨后 的 最高地下水位 ， 此后不考虑降雨对坡体水位的影响。 该坡 体 的模 型示意 图 如图 3所示 。 图 3 边 坡渗流场计算模型示意 图 Fi g ur e 3 S l o p e s e e p a g e c a l c u l a t i o n mo de l da g r a m 为了研究虹吸排水后坡体内地下水位的分布情 况分 为 以下 工况进 行考 虑 ： 工 况 一 ： 在坡 面 上 不 同 的位 置 设 置相 同长 度 和 倾角的排水孑 L ， 研究 同一时间的排水情况及暂时失 效时的降水范围和时间。 工况二 ： 在坡面的同一位置设置长度 和倾角不 同的排水孔 ， 研究相同时间时的浸润曲线分布情况。 工况三 ： 虹吸排水孔的位置一定时 ， 研究不同渗 透系数对虹吸排水后地下水分布的影响。 实际 工程 中边 坡虹 吸排水 孑 L 位 置 的选 择 是不 可 避免的 ， 也是意义重大的， 故本文首先对此进行深入 研究 。 工况 1： 相同的时间下。 在边坡的不 同位置设置 3个虹吸排水孔， 其坐 标分别 为( 1 2 ， 2 5 ) ， ( 1 8 ， 2 0 ) ， ( 2 4 ， 1 5 ) ， 分别标记为 A， B， C点， 对应的排水孑 L 长均为 1 4 1 4 m， 每个孑 L 与 竖线夹角为 4 5 。 ， 故每个孔的虹吸管扬程都为 1 0 m， 其布置如图4所示 。由于边坡虹吸排水后地下水分 布受多种因素影响。故在本工况中应采用控制以下 因素 ： 虹吸排水单位时间的流量均为 5 m。 h ； 虹吸排 水 时间均 为 2 4 h 。 - - 一 图 4虹 吸排 水 孔 的 分 布 图 F i g u r e 4 Di s t r i b ut i o n o f s i ph o n dr a i n a g e ho l e s 从以上虹吸排水后地下水水头分布的云图可以 看出排水位置附近的水头变化的较为剧烈。同时由 A排水 位 置到 C排 水 位 置 ， 地 下水 的影 响 范 围逐 渐 减少 ， 排水孑 L 附近水头降也逐渐减少 。 由以上的水头分布可以求得 4种不同排水孔位 置下的浸润曲线分布图如图 5所示。 2 4 2 8 3 2 3 6 4 0 4 4 4 8 5 2 距 离 IT I 图 5虹吸排水后浸润 曲线分布 图 F i g u r e 5 Di s t r i b u t i o n o f i n f i l t r a t i o n c u r v e a f t e r s i p h on d r a i n a g e 从 图中可以计算出不同排水位置在进行虹吸排 T T ， 如 拍 控加 H g、 饵恒繁 舞 T 上 学兔兔 w w w .x u e t u t u .c o m 第 5期 张世华 ， 等 ： 基 于 m a t l a b的边坡虹吸排水渗流场分析 水后地下水位下降的范围 s即原地下水位与排水后 浸润曲线所夹的面积 。A， ， c的 s分别为 3 7 2 3 0 2 m ， 2 4 2 9 5 2 m ， 1 8 7 2 6 4 m 。从排水孔 的布置 图 看 出 点排水孔离坡脚最远 ， 入水深度也较小 ， 在 2 4 h后排水孔 附近的地下水位 下降 的很 大。在虹 吸排水管未暂时失效时较其他位置较优 。 各虹吸排水孔达到最大降深时( 见图 6 ) 。 图 6最终 浸 润 曲线 分 布 图 Fi g u r e 6 Di s t r i b ut i o n o f t he fin a l i n f i l t r a t i o n c ur v e 由于虹吸排水是一种主动免动力排水 ， 当地下 水位下降到排水孔 的储水管的位置时， 虹吸排水暂 时失去作用 ， 停止排水， 工况 1中 A、 、 c位置 的失 效 时间分别约为 3 6 、1 2 0 、 3 8 4 h ， 最大降水影响范围 分别为 6 6 0 5 7 8 、1 4 2 2 5 7 8 、2 0 5 6 5 7 8 m 。可 以 看 出由坡顶 到坡脚虹吸排水孔的失效时间延长 ， 降 水影响范围加大。由此可知 ： 将虹吸排水孔布置在 距离坡脚一定位置的 c点对 提高边坡 的长期稳定 性 效 果最 好 。 综合可知虹吸排水在提高边坡稳定性方面应考 虑长期排水管布置和应急排水管布置。长期虹吸排 水管宜布置在离坡脚较近的 C点 ， 而应急排水管宜 布置在离坡脚一定距离 的 A点。 虹 吸排水孔 布置的时候还需要考 虑倾 角的问 题 ， 在实际工程中钻孔的倾角及孔 的长度决定 了该 孔的费用， 本文仅从排水的效果在工况 2中对排水 孔倾角进行了研究。 工 况 2： 在边坡的同一位置分别设置了 4个不同倾角下 相 同扬 程 的 排 水 孔 ， ( 2 0 ， 1 3 ) ， ( 2 5 ， 1 3 ) ， ( 3 0 ， l 3 ) ， ( 3 5 ， 1 3 ) ， 其分别标记为 ， N， O， P点 ， 对应 的布置 图如 图 7 ， 图 8所 示 。控 制 每个 虹 吸 排 水 孔 的 出水 流量均为 5 m h 。研究各个 虹吸排水孔降水影响 范 围随时 间的 变化关 系 。 上图可以看出： 各个虹 吸排水孔 降水影响范围 s都随着时间的增长而逐渐增大。同时可 以看出随 着入水深度的增加 ， 在相 同的排水时间下降水影 响 图 7 不 同倾 角下 同一位置排水孔的布置 图 Fi g u r e 7 La y o ut of d r a i n a g e h o l e s u nd e r t h e s a me l o c a t i o n a nd di f f e r e n t i n c l i n a t i o n 0 4 8 1 2 l 6 2 O 2 4 2 8 3 2 3 6 4 0 4 4 4 8 5 2 时 间 T h 图 8降水 影响范围与时间的关系曲线 F i g u r e 8 P r e c i p i t a t i o n a f f e c t e d a r e a v s t i me 范围逐渐减小 ， 与工况一 的情况有较好 的吻合。因 此在选定坡面排水位置的情况下在最高扬程的控制 下可以通过变动倾角来使该虹吸排水孔的入水深度 发生变化 ， 从而决定该排水孔作为应急排水孔还是 长期排水孔。 以上工况都是基于同一渗透系数下对虹吸排水 后地下水的分布情况进行 的研究 ， 为研究不同渗 透 系 数对 虹 吸排水 后地 下水 分布 的影 响还 需要进 行 工 况 3的研 究 。 工 况 3： 为了研究不 同渗透系数对虹吸排水后地下水的 分布的影响( 见 图 9 ) 需要控制 以下 因素 排水孔 的位置为( 2 5 ， 1 3 ) ； 虹 吸排水的出水速率是一个 定值为 6 2 5 m h ； 排水时间为 1 2 0 h 。本文研究 了渗透 系数在 0 6 ， 1 2 ， 2 4 m d共 3种 情况 下虹 0 5 l 0 1 5 2 0 2 5 3 0 3 5 4 0 4 5 5 0 距离 m 图 9地下水浸润 曲线分布图 Fi g u r e 9 Di s t r i bu t i o n of g r o un d wa t e r i n f i l t r a t i o n c ur v e T 下 上 加 加 ：2 0 目、 恒旧 卜簧 弘 勰 加 H g、 住恒*卜霉 学兔兔 w w w .x u e t u t u .c o m 公路工程 3 9卷 吸排水后的地下水分布。 从图中可以计算 得 0 6 ，1 2 ， 2 4 m d这 3种 情况下虹吸排水降水 影响范 围 s分别为 8 0 5 4 7 7 ， 9 6 5 2 8 9 ， 9 8 5 9 7 7 m ， 这表明随着渗透系数的增 大 ， 降水影响范围也逐渐增大， 同时地下水浸润曲线 在虹吸排水管附近降低 的越少， 因此较高的渗透系 数可以保证虹吸排水的效果达到较好的效果。 5 结论 通过采用 ma t l a b编程求解虹吸排水后地 下水 位 ， 起到实时了解边坡地下水的状态， 为边坡的稳定 性分析提供了基础。同时通过算例 中的 3种工况的 分析可以得出以下的一些结论： 边坡虹吸排水方式是一种免动力排水措施 ， 入水深度较浅 的虹吸排水 管能较快的降低地 下水 位 ， 然而由于虹吸排水管的特点， 当排水孔附近 的地 下水位降低到储水管位置时， 该 虹吸排水管暂时失 去排水作用 ， 也 即其影响的范围是有限的。深部排 水管在排水初期降水范围不如浅排水管好 ， 但是其 失去排水作用所持续 的时间较长 ， 可以起到长期排 水的作用 ， 为下一次降雨抬升地下水位提供 较大的 缓冲 ， 在实际工程 中深浅排水孔可以考虑配合使用。 边坡坡顶 、 坡 中以及坡脚排水各有其特点 ， 坡顶虹吸排水管人水深度较小短期排水效果好 ， 且 坡外集水井位置较易布置， 可以形成较大的压力差 ， 提供较大的排水量 ， 但其容易达到极限状态， 提高边 坡稳定性能力有限。而坡脚集水井位置不宜 布置 ， 虹吸管两侧 的压力差也难以保证 ， 但是入水深度较 大。因此综合坡脚及坡顶特点 ， 本文建议虹吸排水 管可以布置在离坡脚一定位置的坡中， 实际工程中 可 以根据初始实测地下水位进行深入的考虑。 虹吸排水措施在渗透系数较大的坡体 中能 较大限度的降低地下水范围 ， 在渗透性较低 的坡体 中其影响的范围有限， 也较快达到极限状态。 本文的研究仅从虹吸排水后对均质边坡 的地下 水位的影响上进行研究 ， 由于实际边坡中各种复杂 的因素的影响还需要进一步的研究。 参考文 献 ： 1 俞伯 汀， 孙红月 ， 尚岳全 含碎石 粘性 土边坡 渗流系 统的物 理 模拟试验 J 岩 土工程学报 ， 2 0 0 6， 2 8 ( 6 ) ： 7 0 5 7 0 8 2 高小育 ， 廖红建 ， 丁春华 渗流对 土质边坡稳 定性的影 响 J 岩土力学 ， 2 0 0 4， 2 5( 1 )6 9 7 2 3 何忠 明 考 虑 降雨入 渗影 响 的土 体 边坡 稳定 性 有 限元 分析 J 矿冶工程 ， 2 0 1 2 ( 1 ) 4 S t a n i 6， B I n fl u e n c e o f D r a i n a g e T r e n c h e s o n S l o p e S t a b i l i t y J J o u r n a l o f G e o t e c h n i e a l E n g i n e e ri n g ， 1 9 8 4 ， 1 1 0 ( 1 1 ) ：1 6 2 4 一 l 6 3 6 5C h o i ， E S e e p a g e A r o u n d H o r i z o n ta l D r a i n s i n Hi l l S l o p e s J J o u r - n a l o f H y d r a u l i c E n g i n e e ri n g ， 1 9 8 3 ， 1 0 9 ( 1 0 ) ： 3 6 3 1 3 6 8 6 孙红月 ， 尚岳全 ， 中永江， 等 破碎 岩质边坡排 水隧洞效果监测 分析 J 岩石力学 与工程学报 ， 2 0 0 8 ， 2 7 ( 1 1 ) ： 2 2 6 72 2 7 1 7 张永防 ， 张朝林 湘黔线 K 9 3路堑滑坡虹吸排 水工点的试验研 究 J 路基 工程 ， 1 9 9 9 ( 4 ) ： 2 6 3 0 8 陈野鹰 ， 唐红梅 ， 周作茂 三峡水 库港 口岸坡排 水方案 优化研 究 J 水运 工程 ， 2 0 0 6 ( 5 ) ： 5 0 5 4 9 李灿 ， 高彦栋 ， 黄素逸 热传 导问题 的 ma t l a b数值 计算 J 华 中科技大学学报(自然科学版 ) ， 2 0 0 2， 3 0 ( 9 ) ： 9 1 9 3 1 O 史策 热传导方程有限差分 法 的 MA T L A B实现 J 咸阳 师 范学院学报 ， 2 0 0 9 ， 2 4 ( 4 ) ： 2 72 9 1 1 张菊明 潜水非稳定流 的汇点 、 汇线 和汇面 及其在水 文地 质 计算 中的应用 J 水文地质工程地质 ， 1 9 8 0 ( 6 ) ： l 2一l 8 1 2 陈崇希， 林 敏 成建梅 地 下水 动力学 ( 第五版 ) 北京 ： 地质 出 版社 ， 2 0 1 1 ， 2 85 5 ( 上接第 2 4页) 一 2 9 5 7 5 S o n me z H，U l u s a y R Mo d i fi c a t i o n s t o t h e g e o l o g i c a l s t r e n g t h i n d e x( G S I )a n d t h e i r a p p l i c a b i l i t y t o s t a b i l i t y o f s l o p e s J I n t e r n a t i o n a l J o u r n a l o f R o c k Me c h a n i c s Mi n i n g S c i e n c e s 。 1 9 9 9， 3 6( 6 )： 7 4 37 6 0 6 S o n m e z H，G o k c e o g l u C，U l u s a y RA n a p p l i c a t i o n o f f u z z y s e t s t o t h e Ge o l o g i c a l S t r e n g t h I n d e x( G S I )s y s t e m u s e d i n r o c k e n g i n e e ri n g J E n g i n e e r i n g A p p l i c a t i o n s A r t i fi c i a l I n t e l l i g e n c e ， 2 0 0 3， 1 6 ( 3 ) ： 2 5 1 2 6 9 7 张永杰 ， 曹文贵 ， 赵明华 。 等 基于地质强度指标 与区间理论 的 岩 体抗剪强度确定方 法 J 岩土 力学 ， 2 0 1 1 ，3 2 ( 8) ： 2 4 4 6 2 45 2 8 林杭 ， 曹平 ， 李江腾 ， 等 基于 H o c k - B ro w n准则的三维边坡变形 稳定性分析 J 岩土力学， 2 0 1 0 ， 3 1 ( 1 1 ) ： 3 6 5 63 6 6 0 9 刘立鹏 ，姚磊华 ， 陈洁， 等 基 于 H o e kB row n准则 的岩质边坡 稳定性分析 J 岩石力学与工程学报 ， 2 0 1 0， 2 9( s 1 ) ： 2 8 7 9 2 8 8 6 1 0 陈昌富， 翁敬 良基于广义 H o e k - B row n准则边坡稳定性分析 强度折减法 J 中 国地 质灾害 与防治学报 ， 20 1 0 ， 2 1 ( 1 ) ： 1 3 一 l 8 1 1 邹 琅 山岭公路隧道围岩力学参数估计 J 公路工程 ， 2 0 1 3 ， 3 8 ( 2 ) ： 1 7 01 7 3 1 2 M C a i ， P KK a i s e r ，Y T a s a k a ， e t a 1 D e t e r mi n a t i o n o f r e s i d u a l s t r e n g t h p a r a me t e r s o f j o i n t e d r o c k ma s s e s u s i n g t h e G S I s y s t e m J I n t e r n a t i o n a l J o u rna l o f R o c k Me c h a n i c s Mi n i n g S c i e n c e s ， 2 0 0 7 ，( 4 4 ) ： 2 4 7 2 6 5 学兔兔 w w w .x u e t u t u .c o m