摘要
为更好探明潮汐环境下近海基坑围护结构力学响应问题,提出了以实测墙边地下水位数据为基础的模型边界潮汐地下水位振幅反演分析方法,分析了不同振幅潮汐作用下地下水渗流衰减规律,并进一步研究了不同振幅潮汐地下水作用下基坑围护结构的力学响应规律及响应机理。研究表明:基于地下水渗流衰减规律的潮汐地下水位振幅反演分析方法不用对实际近海区域地质构造信息进行详细考虑,即可得到包含基坑的小范围区域分析所需的边界地下水位,在近海工程潮汐影响分析中具有良好的应用价值;靠海侧地连墙与强透水层直接接触区域以及坑中坑钻孔灌注桩更多受余弦型潮波作用,两者水平侧移摆动随时间呈“类正余弦型”变化规律,而靠海侧地连墙其他区域和背海侧地连墙由于自身协同变形和受其他地层不同相位孔压的作用,呈前后不对称的“不规则槽型”变化规律;潮汐作用下,背海侧围护墙体侧移振幅沿深度增大逐渐降低,但靠海侧围护墙体侧移振幅、弯矩增量分布以及内撑轴力与地层透水性密切相关,特别在进行受较大潮汐影响工程的结构设计时,应尽量避免静力工况最大侧移位置和弯矩峰值位置与强透水层位置重合。
近年来,随着填海造陆工程的兴建,近海深大基坑工程数量也逐年增多,不同于传统富水环境中的深基坑工程,近海深大基坑工程不仅会面临高地下水位、强水力补给问题,受潮汐效应影响,场地地下水位还可能存在滞后于潮汐的周期性波动现象,致使近海基坑渗流场表现为非稳态特点,从而造成基坑周围地层有效应力和超静孔压时刻改变,给基坑围护结构和周围地层的稳定性带来不利影响。
在非稳态渗流边界条件基坑问题中,应宏伟及其团队成员针对潮汐影响下基坑抗渗防突涌问题展开了较为丰富的研
而现有的模拟方法普遍存在地质制约条件过多、计算连续性不足、对潮汐水流在地层中的衰减考虑不全面等问题,不能完全反映潮汐作用下相应结构的响应问题,因此,在工程中的实际运用也较为困难。在此背景下,本文以深圳市滨海大道(总部基地段)交通综合改造项目深基坑工程为依托,采用流固耦合方法,将潮汐考虑为随时间波动的动态边界,研究了潮汐地下水在地层中的渗流衰减规律,建立以实测墙边地下水位数据为基准的潮汐边界振幅反演分析方法,并将反演得到的渗流边界代入基坑数值模型,从而最终探明了潮汐作用下基坑围护结构的力学响应规律及作用机理。
潮汐对近海工程的荷载作用可分为两类,①是潮水压力荷载,②是潮汐引发的地下水位升降荷载。潮水压力荷载为工程结构直接受到水压作用(

图1 潮汐压力荷载
Fig. 1 Load of tidal pressure

图2 潮汐引发的地下水位升降荷载
Fig. 2 Tide-induced rise and fall load of groundwater level
针对上述潮汐荷载模拟上,现有文献主要采用方法有二,较早的方法可称为静态简化,即结合潮汐水位升降缓慢特点,将潮汐变化过程近似考虑并分解为多个典型静载工况进行求
与基坑围护结构静力响应类似,对于距海岸一定距离仅受地下水升降荷载影响的基坑工程,潮汐作用下围护结构力学响应仅与墙边较近距离范围内的土体分布和水位变化情况紧密相关,只要保证基坑周边地质情况及地下水位变化情况的准确模拟,就能保证潮汐作用下基坑围护结构响应的准确模拟。因此,可以实测墙边地下水位数据为基础,通过潮汐地下水渗流衰减规律反演分析的方法,在包含基坑的小范围分析区域内对地下水位动态条件进行较为准确模拟。该方法的基本思想为:通过地连墙周边水位观测孔获得大量实测数据后,将分析模型按坑周地质信息进行水平简化并延伸建模后,按经验预设多组边界条件试算工况对模型的地层潮汐衰减规律进行分析拟合,而后再以围护结构边的实测地下水位数据作为评判标准反演数值分析时模型边界应当输入的动态水位振幅大小,即可较为准确地进行潮汐作用下围护结构力学响应的计算分析。
本文分析案例不涉及基坑开挖,在固定渗流边界环境下完成基坑降水开挖后,引入动态渗流边界进行计算。在潮汐动态水位边界的处理上,结合潮汐时变特征可近似用三角函数描述的特点,首先将反演和后续结构响应计算中的动态水位边界简化为周期12 h的余弦函数如
(1) |
式中:h为模型边界处t时刻潮汐水位,m;A为边界潮汐振幅,m;t为计算时间,s;H为场地最高水位,m。

图3 潮汐水位函数模拟方法示意
Fig. 3 Simulation method of tidal water level function
在进行反演中,先拟定几组边界潮汐振幅A工况施加于靠海侧边界进行试算,提取分析围护结构旁各地层中间厚度处孔压测点振幅响应情况,并与模型边界输入水位振幅拟合出关系函数曲线,随后将实测的墙边地下水位振幅代入函数中,求得能够让墙边模拟水位振幅与实测水位振幅取得一致的数值计算中需要的边界潮汐振幅,并以此潮汐边界振幅开展相应力学响应分析。显然,如此即可在包含基坑的小范围分析区域内对地下水位动态条件进行较为准确模拟,从而保证潮汐环境下基坑围护结构响应分析的可靠性。
滨海大道(总部基地段)交通综合改造项目位于深圳市南山区、福田区规划超级总部片区及附近,由现状滨海大道(总部基地段)进行下沉扩容改造,下沉改造段西起深湾一路,东至深湾五路,总长约1.5 km (

图4 滨海大道(总部基地段)下沉隧道卫星与实景图
Fig. 4 Satellite and real map of sinking tunnel on Binhai Avenue ( Headquarters base section )

图5 基坑暗埋段典型断面示意图
Fig. 5 Typical section of buried section of foundation pit

图6 地层与基坑位置关系图
Fig. 6 Relationship between stratum and foundation pit position
基坑平行于海岸线由西向东呈长条状布置,故分别于基坑西、中、东段三处区域靠海侧地连墙外 2 m布设地下水位观测孔,进行自动化监测,由西向东分别为DSW-1~DSW-3,满朝时刻距海岸线距离约54~58 m,观测孔位置及时间信息如

图7 地下水位观测孔信息图
Fig. 7 Information map of underground water level observation holes
观测孔 | 观测时间段 | 观测天数 |
---|---|---|
DSW-1 | 2020/10/15~2020/11/24 | 41 |
DSW-2 | 2020/11/04~2020/11/25 | 21 |
DSW-3 | 2021/05/15~2021/08/30 | 108 |

图8 靠海侧观测孔地下水位数据实测曲线
Fig. 8 Measured curve of groundwater level data from seaside observation hole
三观测孔位置与海岸线距离并无明显区别,但是各孔数据规律却存在较大差异。而解读基坑南辅道地质剖面图(

图9 南辅道坑底线纵断面与基坑平面位置关系图
Fig. 9 Relationship between vertical section and plane position of pit bottom line of south auxiliary road of foundation pit
而滨海基坑总长近1.5 km,且截面形式比较单一,在已探明基坑中、东段基本不受潮汐影响后,进行全长范围内的基坑建模会耗费大量时间且无意义,故本文研究范围仅为基坑西段ZN(S)K1+485~ZN(S)K1+698里程范围。并按西侧地下水位最大差值0.79 m取整设置墙边地下水位差0.4 m为围护结构力学响应分析工况,模型最高水位则按最不利取为历史最高地下水位-2.4 m。
由于滨海基坑形状狭长,截面特征变化不大,选取长度为24 m的部分区段,采用FLAC3D有限差分软件依据实际结构型式进行建模分析,如
(2) |
式中:D为钻孔灌注桩直径; s为桩的净间距;w为地下连续墙的折算厚度。

图10 计算模型示意图
Fig. 10 Calculation model
基坑涉及横撑结构5道,分别为主坑砼撑1~3,坑中坑钢管撑4~5。模型土体服从Mohr—Coulomb屈服准则,地连墙、砼撑、格构柱以及立柱桩均采用弹性实体单元模拟,钢管撑采用beam单元模拟。力学边界条件为水平位移和底部竖向位移约束,模型的水力边界条件按靠海侧边界透水、其余边界不透水进行设置,降水井采用“以线代井法”模
地层名称 | 密度ρ/(kg∙ | 变形/弹性模量 E/MPa | 黏聚力c/kPa | 内摩擦角φ/(°) | 泊松比ν | 垂直渗透系数Kv/(cm∙ | 水平渗透系数Kh/(cm∙ | 孔隙比e |
---|---|---|---|---|---|---|---|---|
压实填土 | 1 880 | 15 | 14 | 18 | 0.26 |
8.45×1 |
5.95×1 | 0.718 |
淤泥 | 1 708 | 3.5 | 3.62 | 18.42 | 0.40 |
6.33×1 |
5.08×1 | 1.390 |
含黏性土砾砂 | 1 948 | 30 | 4.88 | 30.56 | 0.29 |
6.18×1 |
4.49×1 | 0.591 |
砾质黏性土 | 1 856 | 25 | 23.82 | 25.36 | 0.30 |
6.00×1 |
5.71×1 | 0.954 |
全风化花岗岩 | 1 895 | 65 | 4.18 | 30.98 | 0.29 |
6.26×1 |
1.25×1 | 0.806 |
土状花岗岩 | 1 900 | 160 | 40 | 32 | 0.28 |
1.25×1 |
1.63×1 | 0.693 |
钢筋混凝土 | 2 500 | 31 500 | 0.30 | |||||
钢管撑 | 7 800 | 200 000 | 0.20 |
试算反演工况为:动态水位边界振幅2.0 m、2.5 m、3.0 m、3.5 m和4.0 m,欲得到的力学响应分析工况为:墙边水位振幅0.4 m。
以靠海侧地连墙边强透水层(含粘性土砾砂层)中间厚度处孔压取得峰值时刻为起点,绘制各反演工况下靠海侧地连墙边各地层中间厚度处的孔压时程曲线如

图11 各孔压边界条件下的墙边地层中间厚度处测点孔压变化曲线
Fig. 11 Variation of pore pressure at the measuring point at the middle thickness of the wall stratum under the boundary conditions of pore pressure
由3.2节分析可知,在仅含单层强透水层的地质条件下,基坑周边各地层孔压波动主要出现在强透水层,其余地层孔压波动不明显,故可近似以墙边强透水层孔压对观测孔中的水位波动进行反演。
提取各工况下靠海侧强透水含粘性土砾砂层在距计算模型边界不同距离处的孔压振幅,如
(3) |

图12 强透水层孔压振幅与距计算模型边界的关系曲线
Fig. 12 Relationship between pore pressure amplitude and distance from calculation model boundary of strong permeable layer

图13 边界至边墙含水层孔压振幅关系曲线
Fig. 13 Relationship of pore pressure amplitude from boundary to sidewall aquifer
故根据

图14 靠海侧强透水层孔压时程曲面图
Fig. 14 Time-history surface diagram of pore pressure near seaside strong permeable layer

图15 0.4 m振幅潮汐作用下围护结构侧移三维时程曲面图
Fig. 15 3-D time-history curve of lateral displacement of retaining structure at an amplitude tide of 0.4 m

图16 摆动幅度与地连墙深度关系曲线
Fig. 16 Swing amplitude versus depth of diaphragm wall

图17 各围护结构5等距测点侧移时程曲线
Fig. 17 Lateral time history curves of 5 equidistances
进一步观察
图

图18 各潮时围护结构弯矩及弯矩增量曲线图
Fig. 18 Bending moment and bending moment increment of enclosure structure at different tides
虽然墙体弯矩曲线线性主要取决于墙背土压力、内撑和坑底土体的共同作用,但上述分析表明,若地连墙弯矩峰值出现在强透水层区域,潮汐作用下的弯矩增量会相较于其他区域更大,其弯矩不利将进一步加剧。
各道横撑最大轴力随潮汐响应曲线如

图19 各撑轴力最大点轴力随时间变化曲线
Fig. 19 Axial force at maximum points versus time
结合

图20 内撑与强透水层关系示意图
Fig. 20 Relationship between internal bracing and strong permeable layer
本文针对近海基坑围护结构在潮汐荷载作用下的力学响应问题,提出了基于地下水渗流衰减规律的潮汐地下水位振幅反演分析方法,并以此为基础,研究了不同振幅潮汐地下水作用下基坑围护结构的力学响应规律及响应机理。通过研究,得到了如下结论:
(1)近海潮汐作用下,地连墙边孔压发生波动的地层主要为强透水层以及相邻弱透水层,且与强透水层相邻弱透水地层中的孔压波动也主要由强透水层孔压波动向上下相邻地层的传递所造成。强透水层内部测点的孔压振幅与边界孔压呈强线性相关关系,以实测墙边孔压为基础,通过反演分析,可得到包含基坑的小范围数值分析区域所需的合理边界地下水位条件。
(2)靠海侧地连墙与强透水层直接接触区域以及坑中坑钻孔灌注桩更多受余弦型潮波作用,两者水平侧移摆动随时间呈“类正余弦型”变化规律,而靠海侧地连墙其他区域和背海侧地连墙由于自身协同变形和受其他地层不同相位孔压荷载作用,其侧移时程曲线则呈现出前后不对称的“不规则槽型”变化规律。
(3)潮汐作用下,背海侧围护墙体侧移振幅沿深度增大逐渐降低,但靠海侧围护墙体侧移振幅、弯矩增量分布以及内撑轴力与地层透水性密切相关,表现为:与强透水层接触区域围护墙体侧移振幅更大;墙体静力弯矩峰值位于强透水层区域内时,潮汐作用引起弯矩增量相较于其他区域更大;同型内撑轴力振幅随竖向位置远离强透水层逐渐减小,因此,在进行受较大潮汐影响工程结构设计时,应尽量避免静力工况最大侧移位置和弯矩峰值位置与强透水层位置重合。
作者贡献声明
汪 波:提出概念,方法设计,指导数值建模分析。
高筠涵:论文撰写,方法设计,数值建模分析。
马龙祥:指导论文撰写,修改论文框架,指导数值建模及分析。
王文昊:绘制部分图片,论文修改。
姬 繁:绘制图片,现场配合。
参考文献
应宏伟,章丽莎,谢康和,等.坑外地下水位波动引起的基坑水土压力响应[J].浙江大学学报(工学版),2014,48(3):492. [百度学术]
YING Hongwei, ZHANG Lisha, XIE Kanghe, et al. Pore and earth pressure response to groundwater fluctuation out of foundation pit[J]. Journal of Zhejiang University(Engineering Science), 2014,48(3):492. [百度学术]
应宏伟,王迪,许鼎业,等.动态承压水作用下考虑土体非线性的基坑弱透水层出逸比降研究[J].浙江大学学报(工学版),2020,54(12):2356. [百度学术]
YING Hongwei, WANG Di, XU Dingye, et al. Analysis on exit gradient of aquitard at bottom of foundation pitunder dynamic confined water considering nonlinearity of soil[J]. Journal of Zhejiang University(Engineering Science) ,2020,54(12):2356. [百度学术]
应宏伟,许鼎业,王迪,等.波动承压水下基坑底部弱透水层的非Darcy渗流分析[J].上海交通大学学报,2020,54(12):1300. [百度学术]
YING Hongwei, XU Dingye, WANG Di, et al. Analysis of non-darcy flow in aquitard at bottom of foundation pit under fluctuation of confined water[J]. Journal of Shanghai Jiao Tong University, 2020,54(12):1300. [百度学术]
郑刚,曾超峰.基坑开挖前潜水降水引起的地下连续墙侧移研究[J].岩土工程学报,2013,35(12):2153. [百度学术]
ZHENG Gang, ZENG Chaofeng. Lateral displacement of diaphragm wall by dewatering of phreatic water before excavation[J]. Chinese Journal of Geotechnical Engineering, 2013,35(12):2153. [百度学术]
BERESLAVSKII E N. The flow of ground waters around a Zhukovskii sheet pile[J]. Journal of Applied Mathematics and Mechanics, 2011,75(2):210. [百度学术]
孙晓贺,施成华,曹成勇,等.强透水地层水平帷幕作用下基坑渗流特性[J].东南大学学报(自然科学版),2021,51(6):1002. [百度学术]
SUN Xiaohe, SHI Chenghua, CAO Chengyong, et al. Seepage characteristics of foundation pit with horizontal curtain in strong permeable stratum[J]. Journal of Southeast University(Natural Science Edition), 2021,51(6):1002. [百度学术]
陈从睿. 潮汐作用下围堰与支护结构受力特性研究[D]. 成都:西南交通大学,2013. [百度学术]
CHEN Congrui. The study on mechanical characteristics of cofferdam and supporting structure under the effect of tide[D], Chengdu:Southwest Jiaotong University, 2013. [百度学术]
寇强. 潮汐动力作用下深大基坑渗流场演化特征及支护对策研究[D]. 济南:山东大学,2016. [百度学术]
KOU Qiang. The study on seepage field evolution characterists and supporting technology of large and deep foundation pit under the action of tiday dynamic[D]. Jinan:Shandong University, 2016. [百度学术]
李健,范方方.临海环境基坑渗流特性研究[J].水利水电技术,2019,50(7):202. [百度学术]
LI Jian, FAN Fangfang. Study on seepage characteristics of foundation pit in seaward environment [J]. Water Resources and Hydropower Engineering, 2019,50(7):202. [百度学术]
应宏伟,聂文峰,王啟铜,等.水位波动对临海重力式挡墙基坑稳定性的影响[J].海洋工程,2013,31(4):48 [百度学术]
Ying Hongwei, Nie Wenfeng, Wang Qitong, et al. The influence of water level fluctuation on the stability of the foundation pit of the gravity retaining wall near the sea [J]. Ocean Engineering, 2013,31(4):48. [百度学术]
陈昆,闫澍旺,孙立强,等. 开挖卸荷状态下深基坑变形特性研究[J]. 岩土力学, 2016,37(4):1075. [百度学术]
CHEN Kun, YAN Shuwang, SUN Liqiang, et al. Analysis of deformation of deep foundation pit under excavation unloading condition[J]. Rock and Soil Mechanics, 2016,37(4):1075. [百度学术]