南京长江4桥北锚碇沉井的排水下沉分析
2022-08-11
来源:欧得旅游网
第32卷第5期 2O10年1O月 土木建筑与环境工程 Journal of Civil.Architectural&Environmental Engineering Vo1.32 NO.5 0ct.2O10 南京长江4桥北锚碇沉井的排水下沉分析 维明 ,牛亚洲 穆保 岗 ,朱建民 , 龚 交第二公路工程局有限公司,西安710075) (1.东南大学土木工程学院,南京210096;2.中 摘要:结合南京长江4桥北锚碇沉井的施工特点,在现场抽水试验基础上,综合确定了沉井排水 下沉期间的渗透系数,由此理论计算单井出水量、总涌水量等并进行排水设计,实际涌水量略大于 理论值;根据有效应力原理,认为沉井降水对潜水含水层的影响很小,沉降主要由承压水含水层压 缩引起,由此得到的沉降曲线与实测结果总体吻合较好。 关键词:沉井;渗透系数;降排水下沉;涌水量;沉降 中图分类号:TU473.2 文献标志码:A 文章编号:1674—4764(2010)05—0135—07 Dewater Sinking Analysis of North Anchorage Caisson for the Fourth Nanj ing Yangtze River Bridge MU Bao—gang ,ZHU Jian—m『n ,GONG Wei—ming ,NIU ya—zhou。 (1.School of Civil Engineering,Southeast University,Nanjing 210096,P.R.China; 2.CCCC Second Highway Engineering Co.,Ltd,Xi’an 7 10075,P.R.China) Abstract:Combined with construction of north anchorage caisson,field pumping tests were carried out. Permeability coefficient was acquired for dewater sinking.Water yield of single well and total water inflow of caisson excavation were calculated.It is found that the measured amount of total water inflow is slightly larger than the theoretical value.Based on effective stress principle,the compression of phreatic water aquifer caused by dewater is almost negligible,and the ground settlement is mainly composed of confined water aquifer compression.The calculated settlement curve along the levee was in good agreement with the measured data. Key words:caisson;permeability coefficient;dewater sinking;water yield;settlement 南京长江第4大桥为双塔3跨悬索桥,主桥跨 径(166+410.2)+1 418+(363.4+118.4)一2 476 m,北锚碇为重力式锚碇,基础采用图1所示的沉井 结构,其顺桥向长69 m,横桥向宽58 m,高52.8 m, 平面尺寸为世界第一。 沉井可以分为排水下沉和不排水下沉。排水下 沉作业面可见,利于对沉井姿态控制,且费用相对低 廉。降排水下沉的综合效率优势明显,故北锚碇沉 井拟采用前期降排水下沉,后期不排水下沉的施工 方案,并力争在汛期到来之前完成降排水下沉。 图1沉井半结构图 收稿日期:2010—03—15 作者简介:穆保岗(1974一),男,博士,副教授,主要从事岩土工程和地下结构方面的研究, (E—mail)mubaogang@seu.edu.cn。 136 土木建筑与环境工程 第32卷 该沉井平面上共分2O个井孔,标准节段井孔顺 桥向长14.2 m,横桥向宽8.68 m;沉井竖向共划分 11节,首节为钢壳沉井,高6 m,其余节段为混凝土 沉井,第2~10节高5 m,第11节高1.8 m。 州大桥南、北锚碇沉井l_5蜘共3座,目前正准备修建 的马鞍山大桥亦采用大型沉井。由于工程地质和水 文地质条件不同,对周边环境影响的计算,尤其是附 近的长江大堤沉降,均建立在现场监测基础上_7J,对 拟采用沉井接高和下沉方案为:地面接高前4 于缺乏沉降观测资料的场地区域,抽水引起地面沉 降量的估算方法和模型较为繁杂,尚无一套简明有 节(21 m),降排水下沉一接高5、6节(10 m),一接高 7、8节(10 m),不排水下沉一接高9~11节(1】.8 m),不排水下沉。 由于沉井的平面尺寸大,下沉深度大,降水深度 大,地质条件复杂,加之沉井紧临长江,地下水与长 江水互相补给,其降排水设计与施工都具有一定的 难度;沉井南侧120 m即为长江大堤,如何在下沉期 间保证大堤的安全也是沉井降排水设计和施工的重 要考虑因素。 表1为北锚锭的地质情况表,北锚碇区地层大 体可分为上部软弱土层和下部砂土层,并以第9层 密实卵砾石层作为沉井的持力层。 地下水分为潜水和承压水2类,其中潜水水位 埋深约1.8~3.3 m,接受大气降水与田间沟塘水补 给,与长江等地表水体有一定水力联系;承压水水位 埋深约12 m,与长江水力联系密切。 表1北锚碇区地质情况表 粉砂:灰色,饱和,中密,局部松散, …‘夹亚粘土 4.5 细砂:灰色,饱和,中密 6.2 中砂:灰色,饱和,密实 17,1细砂;灰色,饱和,密实,夹亚粘土 承压水 2.7 卵砾石、砾砂:黄灰色,饱和,密实 4.3 细砂:灰色,饱和,密实 5.2 卵砾石、砾砂:黄灰色,饱和,密实 2.6 强风化粉砂岩 国外建造的大型沉井较少,而偏重于室内模型 试验 和小型沉井 ,大型沉井主要集中在中国深 水桥梁基础上。中国建造的可供借鉴的同类规模的 陆地沉井,到目前为止有江阴大桥北锚碇沉井 ,泰 效的通用手段而采用个案试验实测的方法的 。 结合南京长江第4桥北锚碇沉井的施工,对降排水 设计和长江大堤沉降计算2个问题进行了分析。 1现场抽水试验 1.1 试验概况 对北锚碇区水文地质条件和工程特点的分析可 知,地下水对沉井下沉影响的主要为上部承压水,为 了获得准确的渗透系数和确定抽水工艺,针对上部 承压水设计了抽水试验。设抽水井1个(JS06),观 测井2个G1、G2,G3为备用井,具体布置见图2。 图2试验井布置图 主井和观测井分别采用 ̄p600 mm、 35O mm 3 翼螺旋合金钻头施工,一次性成孔。循环液采用自 造浆,泥浆较稀,利于洗井,试验井参数见表2。 表2抽水试验井参数表 抽水试验共进行3个降深,试验时间为2008年 1O月15日至1O月19日。各次降深分别为2.47、 8.59和14.78 1TI,出水量分别为20 rn。/h、80 m。/h 和125 m。/h,稳定时间分别8 h、8 h、12 h。 抽水历时曲线见下图3、4、5,抽水试验数据见 第5期 穆保岗,等:南京长江4桥北锚碇沉井的排水下沉分析 137 表3。 图3 20 m /h抽水历时曲线 图4 80 m /h抽水历时曲线 图5 25 m /h抽水历时曲线 表3主井和观测井抽水试验结果 流量Q和降深S关系见图6。 上述曲线表明,试验井成井质量较好,主井出水 量大,观测井反映灵敏。 1.2渗透系数(K)的计算 1.2.1 采用稳定流求参承压水非完整井(含水层 上部)单井渗透系数计算公式[9]为 流量/m ,}l 图6 s曲线 K一 g a 6 (1) 承压水非完整井(含水层上部)多孔渗透系数计 算公式为: K一 (arsh,L .-arsh / ) (2) 其中:K为渗透系数(m/d);Q为涌水量(rn。/d) 、 r:为抽水孔至观测孔之问的距离(m);S 、S 为观测 孔内的水位降深(m);z为过滤器有效长度(m)。 代人数据,计算结果见表4—6。 表4按照承压非完整井单井计算结果 表5承压非完整并带一个观测井计算结果 表6 承压非完整井带2个观测井计算结果 可见,采用稳定流求参的方法,前2种计算结果 较为接近,而采用2个观测井的计算结果明显偏大。 138 土木建筑与环境工程 第32卷 1.2.2采用非稳定流求参 1)直线解析法 水层可视为一个各向同性的承压含水层;在该含水 层中由于地层岩性的分布不均匀,沉积厚度的不同, 又造成含水层在不同地段含水层渗透系数相差较 大;因此,该含水层应视为一个非均质各向同性的承 在抽水流量为定值时,且观测井距抽水井的距 离r已测定后,降深和时间是对数函数关系,绘在单 对数纸上应是一条直线_】 。设在t 、t:时测得观测井 的水位降深为s 、S:,则有: s 一 lg ㈦ s:一2 .3Qlg (4) 解上2式得: s2一S 一 1g tL (5) ,当取t2=10t 时,则有: T一 △S AS—S2一S (6) ~ 再根据公式K—T/M求得承压水含水层的渗 透系数,结果详见表7。 表7直线解析法计算结果 孔号 O M tl £2 S1 S2 K/(m・d一 )平均 Gl 31.9 20 0.31 0.12 G2 114.9 20 0.38 0.26 Gl 195.6 20 0.42 0.34 G2 31.9 20 1.04 0.91 G1 114.9 20 0.58 0.44 G2 195.6 20 1.52 1.40 2)水位恢复法(两点法) 设在水位开始恢复后t 、t 时测得观测井的水 位降深为s 、sz,当“<0.1时,则有 T一 ・g ㈩ 由K—T/M可得到所求的含水层渗透系数K, 详见表8。 表8两点法计算结果 孔号 Q M fl 5I S2 K/(m・d一1)平均 G1 31.9 55 G2 114.9 55 G1 195.6 55 G2 31.9 55 G1 114.9 55 G2 l95.6 55 1.3计算结果讨论 根据南京长江沿岸水文地质勘察资料,平行于 长江方向与垂直于长江方向的含水层渗透系数变化 较小,其中平行于长江方向的含水层渗透系数略小 于垂直长江方向的含水层渗透系数,从整体上该含 压含水系统。 分析表明:针对拟建锚碇区的水文地质条件,选 择不同的数学计算公式,所计算的结果存在一定的 差异,渗透系数K在13.42~36.99 m/d范围内,用 抽水井参与计算的数值均偏小,考虑到抽水井近距 离周边流速快,呈紊流状态,形成的水跃值较大,因 此计算结果比实际偏小将以上各种算法累加求和, 取其平均值为19.9 m/d,接近采用非稳定流解得的 平均值。 2 沉井降水计算 2.1沉井涌水量 沉井总涌水量由潜水含水层和承压水含水层2 部分组成,各含水层的涌水量可按“大井法”计算。 沉井涌水量计算等效半径 ro一0.29(口+b)一36.83 m (8) 1)潜水含水层涌水量 本含水层渗透系数取K=1.84 m/d,含水层厚 度H一6.6 m,水位降深S一6.1 m,降水影响半径: R一2S ̄/HK一42.51 m (9) 涌水量可按均质含水层潜水完整井,降水远离 边界的情况计算El2]: Q一1.366 K 一326.6 m。/d(1o) lgf\ 1+ 1ro/ 可见其涌水量不大,可不作为考虑重点。 2)承压水含水层涌水量 含水层渗透系数取K:19.9O m/d,含水层厚度 M一55 1TI,承压水静止水位埋深1.1 ITI,考虑开挖深 度20 m,降水面在开挖面以下1.0 m,则沉井中心水 位降深s一19.9 m,降水影响半径 R一10S、/ 一887.7 133. (11) 按照基坑规范验算,涌水量可按均质含水层承 压一潜水非完整井计算I】 : Q一1.366K—(2HM) M --h2 ——一lgf\. _1+ 1r0, 一39 132.5 1TI。/d (12) 按照基坑手册验算,涌水量可按均质含水层承 压一潜水非完整井,岸边降水模型计算Il ,近似取 井中心至直线补给边界距离为120+36.83— 156.83 m,则有: Q一—L—_ l(2H"一he)M"— 一(寺)I + ln—— 第5期 2 ̄rKTS穆保岗,等:南京长江4桥北锚碇沉井的排水下沉分析 2 叭 ——:・—一一 37 909.4 s/d. 4 m/d .一in i20+ (13)Ll J 并参考以往工程经验[14 16],降水井分2排布设于沉 m 井环道外侧,距沉井侧边的距离分别为17 m和22 In,见图7。 L 2 可见承压水含水层涌水量很大,且按2种方法 的计算结果较为接近,取其中较大者控制降排水设 计。 2.2计算参数讨论 需要说明的是,含水层影响半径R并不是固定 不变的,随管井的水位降深的增大,含水层的影响半 径也随之增大。 有2个观测孔计算承压含水层影响半径的公 式: lgR一 …) 式中,R为影响半径,m;S 、S 为观测孔内水位降 深,m;r 、r 为抽水孔至观测孔之问的距离,rn。 按照上述实验结果计算如表9。 表9 2个观测孔计算影响半径结果 S1 S2 r1 r2 Q R 由于抽水试验工程紧邻长江大堤,抽水含水层 顶板埋深仅有12.0 m,长江河床的厚度远超过了该 埋置深度,长江河床穿过了该含水层的顶板,形成的 地下水与长江水之间的水力联系较为密切,即一旦 抽取该地段的地下水,长江水就构成了定水头补给 边界。可以判定,在抽水试验该地区,其地下水的影 响半径是不对称的,因此,随着抽水水位下降,其影 响半径是不断的向外扩展,而南部只能扩展到长江 水域。 2.3降水井布设 1)管井出水量 按沉井外侧10 m左右布井,取1:10的水力坡 降,降水井处水位埋深约25 m;降水井过滤器置于 承压水含水层,深度为13~34 m,管径r 一325 mm, 暂不考虑井损,则有效过滤器长z===9 In,得到单井出 水量 q一120nr z ===1 493 ma/d (15) 2)降水井数量 降水井的井点数 一1.1Q/q一29 (16) 考虑到计算参数和计算方法的可能偏差,总井 数最终取为32。 平面布井时应做到近大堤侧稀,远大堤侧密,以 便在能够满足沉井降水要求的前提下,使降水对长 江大堤产生的影响最小。最终根据现场条件情况, 加 弘 驭 。降水井 虬 ・沉降观测点 9 Ⅲ 9 6 图7大堤位置 姐 3长江大堤沉降计算 当降水引起地下水位下降时,土中的孔隙水应 力减小,而总应力基本维持不变,由此将导致土体中 的有效应力增加,土体骨架压缩引起地面沉降。目 前普遍采用Terzaghi(1925年)提出的有效应力原 理,认为地下水压力的降低导致有效应力的增大是 含水层压密的原因。 对于该场地而言,承压含水层基本为砂土层,含 水砂层的变形与水位变化相关性较好,一般认为其 峰值能互对应 。 。含水砂层的压缩量较小,满足 弹性变形规律,可用简约的分层总和法求解。 3.1 潜水含水层压缩量 对潜水含水层而言,沉井的降水影响半径仅为 42.51 m,远小于沉井到长江大堤的距离120 m,即 该处潜水水位不会因沉井降水而改变;若隔水层效 果好,不存在潜水和承压水的互相补给,则承压含水 层降水也不会引起潜水水位的变化;这2种情况都 不会使潜水含水层的有效应力改变,从而也不会产 生压缩沉降。 较为不利的情况则是隔水层效果不好,潜水向 下补给承压水造成其水位下降,但因为上部软土固 结的竖向渗透系数很小,这一过程通常发生的很缓 慢,而沉井的降水持续时问只有半个月左右,而且是 分阶段降低地下水位,即使土层产生压缩固结,其值 也会很小。 综上所述,可忽略潜水含水层在降水期间的沉 降。 3.2承压水含水层压缩量 在沉井降水过程中,长江大堤附近水位降深仍 将处于承压顶板之上(即水压下降,但静止水位仍高 于隔水层顶板),若水头下降了 ,则整个承压含水 层孔隙水应力均降低r h,其有效应力均增加 r ̄21h。 由于整个承压水含水层以砂土为主,产生的压 14O 土木建筑与环境工程 第32卷 缩量为土骨架的压缩,并且会很快完成,因此当作为 沉降计算的重点。 1)水位降深计算 随着下沉深度的不断增加再逐步加大降水力度,沉 井下沉到位时最终开启了25口管井降水,根据各井 的工作状况统计最大日抽水量约42 000 m。。 对于本工程的降水特点而言(承压~潜水非完 整井),现有理论尚无法准确处理水位降深问题,考 虑到沉井离大堤尚有一定距离,故可以近似取完整 从降水方案和实际情况的比较来看,设计方案 总体上是成功的,它可为后续的类似工程提供指导 和借鉴;沉井总涌水量计算结果比实际情况偏小 8 左右,单井出水能力计算公式的结果比真实值偏 小约12 ,由于布井时的富余量,仍保证了降水目 井计算 ] —S一 2 ̄KM骞Q ・n(鲁)) (17) 标的顺利实现。 式中Q 为单井涌水量,按总涌水量取平均.rJ为 降水井到计算点距离;R 为降水井影响半径。 2)压缩量计算 按(17)式计算出大堤沿线点水位降深△ 后,即 得该处有效应力增量r Ah,进而可根据土体的压缩 模量求解总压缩量。 AS=奎r ̄,Ah (18) i=1 考虑到计算点埋深、土体原位应力状态以及原 状土的结构强度等因素时,现场压缩模量E 值要大 于室内试验的结果。现把承压水含水层做为整体计 算,根据北锚碇区地质土层的标贯击数N与压缩模 量问的验公式,可以取整体加权平均的砂层压缩模 量 ,取E 一1OO MPa,得到图8中所示的理论沉 降曲线。 3.3 长江大堤沉降控制 目前中国堤防运行管理期变形控制尚未形成统 一标准,按照相关经验和规定,对重要堤防控制沉降 差小于1 ,可以避免不均匀沉降引起的堤防开裂, 因此大堤沉降控制标准取为堤身高度的1 ,即45 mm作为警戒值。 为保证长江堤防的安全,一方面施工现场应以 观测为主,加强降水过程中的大堤监控,确保沉降不 超标;另一方面,若长江大堤最大沉降接近警戒值, 而沉井尚未下沉到位,则可采用回灌井等应急措 施 。 4降水和沉降监控结果 4.1降水监控情况 沉井前4节21 m接高完毕后,从2009年5月 18日开始首次降排水下沉,到6月2日下沉结束, 历时15 d,沉井下沉深度19.6 m,下沉量15.6 m。 最大日下沉量1.95 m,最小日下沉量0.06 m,平均 日下沉量1.04 m。 下沉到位时的沉井顺桥向倾斜度]/335,横桥 向倾斜度1/707,沉井扭转角3"6”,沉井顺桥向最大 偏位53 mm,横桥向最大偏位91 mm,几何姿态正 常。 沉井降排水下沉时,首先开启周边8口井降水, 造成涌水量计算结果偏差的原因可能在于K 的试验取值方法,也可能在于计算公式本身的误差, 以及群井间的互相干挠等,可积累更多资料后再详 细研究。 4.2长江大堤沉降监控结果 沉井下沉过程中每日两次对长江大堤沉降进行 了监控测量,2009年6月2日沉井下沉到位时,沉 降也发展到了最大程度,图8为危险地段500 ITI范 围内的堤身沿线沉降曲线对照图,虚线圈出范围内 的测点因受施工便道和往来车辆碾压影响,其降水 引发的沉降当比给出的测量值要小。从沉降预测和 图8堤身沿线沉降比较 实际情况的比较来看,事先根据有效应力原理计算 的沉降趋势和数值大小还是较为吻合的,计算思路 和计算方法可为后续的类似工程提供指导和借鉴; 预测的沉降最大值约30 mm,位置在D9点附近,监 测结果表明最大沉降为D9点的37 mm,峰值误差 为19 ,由于实测曲线峰值点附近受施工运输车辆 影响较大,但较难扣除这部分的影响量;从两条曲线 形状来看,非峰值部位预测结果略偏大,总体上二者 的吻合程度还是较好的。 5 结论 南京长江4桥北锚碇沉井下沉难度大,对沉井 本身安全控制和对周边大堤影响显得尤为重要。在 进行理论分析计算和现场实测情况对比后,得到如 下结论: 1)深层降水的渗透系数对涌水量、影响半径、水 位降深、大堤沉降等均有直接影响,但其大小受计算 方法影响较大,必须结合本地区的经验数值确定。 就该工程而言,实际涌水量和理论涌水量的差别较 小,可以验证采用非稳定渗流解的平均值较为接近 第5期 穆保岗,等:南京长江4桥北锚碇沉井的排水下沉分析 141 实际情况。 2)北锚碇沉井降排水设计计算表明,按基坑规 范规定计算方法得到的总涌水量比现场量测结果偏 小8 左右,按经验公式计算的单井出水量比现场 值偏小12 左右。 3)从有效应力原理出发计算沉井附近长江大堤 堤身沉降。根据降水计算结果,认为潜水含水层压 缩量很小,大堤沉降主要由承压水含水层压缩引起, 由此得到的沉降分布曲线总体上与实测结果吻合程 度较高,对短期的深层降水引起的地面沉降可以采 用分层综合法进行估算。 参考文献: [1]REESE.Analysis and Design of Shallow and Deep Foundations[M].New Jersey John Wiley&Sons,Inc. 2006. r 2]CHEN W,ZH0U H,RANDOLPH M F.Effect of installation method on external shaft friction of caissons in soft clay[J].The Journal of Geotechnical and Geoenvironmental Engineering,2009,135(5):605—615. r 3]ALLENBY D,WALEY G,KILBURN D.Examples of open caisson sinking in Scotland[C]//Geotechnical Engineering 162 issue GEI,2009. [4]陈光福.江阴长江公路大桥特大型沉井施工述评[J].土 工基础,1999,13(3):24—28,13. CHEN GUANG—FU.The construction of the large caisson in Jiangyin Highway Bridge[J]. Soil Engineering and Foundation,1999,3:24—28,13. [5]陶建山.泰州大桥南锚碇巨型沉井排水下沉施工技术 [J].铁道工程学报,2009(1):63—66. TA0 JIAN—SHAN. Construction technology for draining—sinkage for south caisson anchorage to Taizhou Yangtze River Highway Bridge with large-size sunk well[J].Journal of Railway Engineering Society,2009 (1):63 66. [6]王卫忠.沉井下沉技术在泰州长江公路大桥北锚碇中的 应用[J].交通科技,2009,2(233):21—24. WANG WEI—ZHONG.Application of caisson sinking technology at north anchorage to Taizhou Yangtze River Highway Bridge[J]. Transportation Science& Technology,2009,2(233):21-24. [7]陈缪芬.沉井施工降水时地下水动态分析及对江堤影响 的研究[D].扬州大学,2008. [8]黄雅虹.吕悦军.北京亦庄轻轨工程场地水位下降引起 地面沉降量的评估方法探讨EJ]。岩土力学,2009,8 (3O):2457—2462. HUANG YA—H0NG,LU YUE—TUN.A method for estimating land subsidence induced by groundwater extraction and its application to site evaluation of Yizhuang light railway in Beijing[J]Rock and Soil Mechanics,2009,8(30):2457-2462. [9]王建秀,吴林高,等.地铁车站深基坑降水诱发沉降机制 及计算方法[J].岩石力学与工程学报,2009,28(5): 101O—l019. WANG JIAN—XIU,WU LIN—GA0,et a1.Mechanism of dewatering induced ground subsidence in deep subway station pit and calculation method[J]Chinese Journal of Rock Mechanics and Engineering 2009,28 (5):1010—1019. [1O]唐益群,栾长青,等.上海某地铁站试降水对周边环境的 影响分析[J]武汉理工大学学报,2008,30(8):147—151. TANG YI—QUN,LUAN CHAN(、rQING, et a1. Analysis of the effects of environments for dewatering in a metro station in Shanghai[J].Journal of Wuhan University of Technology 2008,30(8):147—151. [11]姚天强,石振华.基坑降水手册[M].北京:中国建筑工 业出版社,2006. [12]ALLAN R,Cherry A JOHN Groundwater[M].New Jersey,Prentice—Hall lnc,1979. [13]JGJ 120—99建筑基坑支护技术规程Is]. [14]刘建航,候学渊.基坑工程手册[M].北京:中国建筑工 业出版社,1997. [15]夏明耀,曾进伦.地下工程设计施工手册[M].北京:中 国建筑工业出版社,2002. [16]汪jE华.润扬长江大桥南锚深基坑降水及周边沉降控制 研究[D].南京:河海大学,2004 [17]LEUNG C F,LEE F H,KHOO E.Behavior of gravity cassion on sand[J].The Journal of Geotechnical and Geoenvironmental Engineering,1997,123(3):187—196. rl8]CRAWFORD C B,MCCAMMON N R,BUTLER R C. DEEp-Seated conso1idation settlements in the Fraser River delta[J].Can.Geotech.J.1991,28(2):298—303. [19]王翠英,黄理兴,段卫昌,等.深基坑降水引起周边地面 沉降量值计算修正系数Ms的确定[J].岩土力学, 2006,27(6):1011 1016. WANG CUI-YING,HUANG LI—XING.DUAN WE CHANG,et a1.Determinationof coefficient Ms during calculation of surrounding ground settlement due to foundation pit dewatering[J]. Rock and Soil Mechanics,2006,27(6):1011-1016. [2O]吴林高.工程降水设计施工与基坑渗流理论[M].北京: 人民交通出版社,2003. (编辑 胡 玲)