有限平面内基于镜像法海底隧道非线性渗流场分析方法技术

本发明专利技术公开了一种有限平面内基于镜像法海底隧道非线性渗流场分析方法，包括确定海底隧道开挖断面的地质条件，地质条件为海底隧道开挖断面位于均质地层内，将均质地层视作环向分布于海底隧道断面外围，并确定各均质地层至隧道断面中心的内/外径与各均质地层对应的渗透系数；将实际海底隧道与虚拟海底隧道进行形成海底隧道渗流场，对海底隧道渗流场计算获得围岩渗流场参数、注浆圈渗流参数与衬砌渗流参数；根据围岩渗流场参数、注浆圈渗流场参数以及衬砌渗流场参数，对海底隧道渗流场进行分析预测；解决了海底隧道渗流场无法考虑下边界不透水的条件下，同时考虑注浆圈与衬砌均为非线性渗流的情况，对海底隧道非线性渗流场进行预测的问题。测的问题。测的问题。

【技术实现步骤摘要】
有限平面内基于镜像法海底隧道非线性渗流场分析方法


[0001]本专利技术涉及土木工程
，尤其涉及一种有限平面内基于镜像法海底隧道非线性渗流场分析方法。

技术介绍

[0002]近年来，伴随着我国交通强国战略实施的深入，滨海城市隧道飞速建设，数量已经十分庞大。但是大部分的隧道或被设计成排水隧道，或因工艺本身存在渗水通道，或在长期服役过程中产生大量的变形裂缝，这些因素都会改变隧道周围的水土压力分布。并且隧道开挖往往会对周围土体产生影响，影响土体渗透系数。
[0003]目前大多将海底隧道的围岩考虑为半无限空间，即不考虑下边界不透水边界对渗流场的影响，与实际不符，且无法考虑下边界为不透水边界，同时考虑注浆圈与衬砌均为非线性渗流的情况下，对海底隧道非线性渗流场进行精确预测，降低了其在隧道施工前后的安全性和适用性。

技术实现思路

[0004]本专利技术提供一种有限平面内基于镜像法海底隧道非线性渗流场分析方法，以克服目前大多将海底隧道的围岩考虑为半无限空间，即不考虑下边界不透水边界对渗流场的影响，与实际不符，且无法考虑下边界为不透本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种有限平面内基于镜像法海底隧道非线性渗流场分析方法，其特征在于，包括步骤S1:确定海底隧道开挖断面的地质条件，所述地质条件为海底隧道开挖断面位于均质地层区域内，且所述均质地层区域的顶端所在平面界定义为透水边界，所述均质地层区域的底端所在平面定义为不透水边界；所述均质地层包括围岩层、注浆圈层以及衬砌层；步骤S2：将均质地层区域视作环向分布于海底隧道断面外围，并确定各均质地层至隧道断面中心的内/外径与各均质地层对应的渗透系数；所述渗透系数包括围岩渗透系数、注浆圈渗透系数以及衬砌渗透系数；步骤S3：设定实际海底隧道中心处的水平面为位势零面，以水位线为x轴参考基准线，以垂直于实际海底隧道轴线的直线为y轴参考基准线，建立直角坐标系；基于镜像叠加原理，以透水边界与不透水边界分别为镜面对实际海底隧道进行映射获取虚拟海底隧道，并将所述实际海底隧道与虚拟海底隧道进行势的叠加，形成考虑围岩的透水边界与不透水边界影响的海底隧道渗流场模型，所述海底隧道渗流场模型包括围岩渗流场、注浆圈渗流场以及衬砌渗流场；步骤S4：基于达西渗流定律对所述围岩渗流场进行计算，获得围岩渗流场参数，所述围岩渗流场参数包括海底隧道断面的涌水量、围岩区总水头以及海底隧道围岩内渗流量；步骤S5：基于Hansbo非线性渗流模型，对所述注浆圈渗流场以及衬砌渗流场进行计算，获得注浆圈渗流参数与衬砌渗流参数；所述注浆圈渗流参数包括注浆圈渗流量、注浆圈总水头以及注浆圈孔压；所述衬砌渗流参数包括衬砌渗流量、衬砌总水头以及衬砌孔压；步骤S6：根据所述围岩渗流参数、注浆圈渗流场参数以及衬砌渗流场参数，对海底隧道渗流场进行分析预测。2.根据权利要求1所述的一种有限平面内基于镜像法海底隧道非线性渗流场分析方法，其特征在于，步骤S4中基于达西渗流定律对所述围岩渗流场计算获得围岩渗流场参数，具体为：步骤S4.1：将所述实际海底隧道以透水边界为镜面进行等流量异号反映获得虚拟海底隧道并定义为源，且反映后透水边界保持为等势线；将所述实际海底隧道以不透水边界为镜面进行等流量同号反映获得虚拟海底隧道并定义为汇，且反映后不透水边界保持为分流线；步骤S4.2:土体的渗流模式符合达西渗流得：其中，i表示围岩的水力梯度；k表示渗透系数；v表示围岩的渗流速度；则当海底隧道为单孔圆形隧道时，流经隧道每一断面的流量为：Q
s
＝2πrv
ꢀꢀꢀꢀꢀꢀ
(2)式中：Q
s
表示围岩层计算断面的涌水量；r表示计算断面任意点到海底隧道圆心的距离,将式(1)带入到式(2)中，由可得：
式中：表示围岩层内任意一点总水头；Q
s
表示围岩层计算断面的涌水量；k
s
表示围岩渗透系数；n
s
表示土体扰动系数；对公式(3)进行积分得到：式中：表示围岩层内任意一点总水头；Q
s
表示围岩层计算断面的涌水量；k
s
表示围岩渗透系数；r表示计算断面任意点到海底隧道圆心的距离；n
s
表示土体扰动系数；C表示积分常数；步骤S4.3：由镜像法平面内任意一点水势头可根据势的叠加原则来确定，因此：式中：Q
si
表示第i个隧道的渗流量，流进时取正值，流出时取负值；r
i
表示计算点到第i个隧道中心点的距离；C1表示为待定常数值；镜像反映后各类隧道的中心坐标包括点汇坐标与点源坐标；所述点汇坐标为(0,4nh+b),(0,2h+4nh
‑
b)；点源坐标为(0,4nh
‑
b),(0,2h+4nh+b)；其中，h表示给水边界与不透水边界之间的距离；b表示实际隧道的圆心至给水边界的距离；n＝0,
±
1,
±
2,
···
,
±
∞；因此，式中：r1表示有限平面内任意一点M(x,y)到第一类隧道圆心的距离；所述第一类隧道为通过镜像法获取的圆心坐标为(0,4nh+b)的虚拟海底隧道；r2表示有限平面内任意一点M(x,y)到第二类隧道圆心的距离；所述第二类隧道为通过镜像法获取的圆心坐标为(0,2h+4nh
‑
b)的虚拟海底隧道；r3表示有限平面内任意一点M(x,y)到第三类隧道圆心的距离；所述第三类隧道为通过镜像法获取的圆心坐标为(0,4nh
‑
b)的虚拟海底隧道；r4表示有限平面内任意一点M(x,y)到第四类隧道圆心的距离；所述第四类隧道为通过镜像法获取的圆心坐标为(0,2h+4nh+b)的虚拟海底隧道；步骤S4.4：将式(6)带入到式(5)中，得到渗流稳定时土体内的任意一点M的水头值为：式中：表示围岩层内任意一点总水头；Q
s
表示围岩层计算断面的涌水量；k
s
表示围岩渗透系数；n
s
表示土体扰动系数；C2表示待定常数；步骤S4.5:根据贝塞特公式：可将公式(7)简化为：
由给水边界条件y＝0时，可求出常数C2＝b；再将注浆圈外边界条件x＝0，y＝b
‑
r
g
时，r
g
表示注浆圈外半径；h
g
表示注浆圈外边界处的总水头；代入式(9)，可得到海底隧道围岩内渗流量：式中：Q表示海底隧道围岩内渗流量：h
g
表示为注浆圈外边界处的总水头；h表示给水边界与不透水边界之间的距离；b表示实际隧道的圆心至给水边界的距离；r0表示衬砌内半径。3.根据权利要求1所述的一种有限平面内基于镜像法海底隧道非线性渗流场分析方法，其特征在于，步骤S5中所述基于Hansbo非线性渗流模型对所述注浆圈渗流场以及衬砌渗流场...

【专利技术属性】
技术研发人员：刘翔，王东池，姜谙男，路通，张瑞，黎奎辰，纪敏龙，
申请(专利权)人：大连海事大学，
类型：发明
国别省市：