悬挂式止水帷幕深基坑降水方案模拟优化方法与流程

技术特征：
1.一种悬挂式止水帷幕深基坑降水方案模拟优化方法，其特征在于：包括如下步骤：s1.根据研究域的地质、水文地质特征和含水层间的水力联系，建立深基坑降水水文地质概念模型；s2.根据水文地质概念模型，取坐标x、y和z轴的取向与含水层介质各向异性的主方向一致，建立三维非稳定渗流数学模型和含水层土体压缩量的地面沉降数学模型，将三维非稳定渗流数学模型与地面沉降数学模型进行耦合，并通过有限元程序研制、求解，用于通过渗流场、应力场的耦合作用，完整模拟研究域地下水水位和地面沉降的动态趋势；s3.根据非稳定渗流与地面沉降三维耦合数值模型，初始化模型参数，结合已知的地下水位观测值和地面沉降量观测值，通过抽水试验数据进行参数反演，获取满足预测精度和置信度且具有良好收敛性和稳定性的非稳定渗流与地面沉降三维耦合数值模型；s4.进行悬挂式止水帷幕和降水井过滤器不同组合深度的基坑降水与地面沉降模拟计算，结合施工技术难度、成本、环境影响，获取最佳深基坑高效降水技术方案。2.根据权利要求1所述的悬挂式止水帷幕深基坑降水方案模拟优化方法，其特征在于：步骤s2具体包括如下步骤：s21.以研究域中心海平面为坐标原点，根据研究域水文地质特性、研究域维护连续墙埋藏深度及抽水井滤水管设置位置，将研究域进行垂直分层；s22.根据水均衡原理和达西定律，建立三维非稳定渗流数学模型，具体包括用于描述研究域地下水运动规律的偏微分方程和定解条件，并通过有限单元法求解；s23.建立研究域地下水位下降引起的含水层土体压缩量的地面沉降数学模型，具体包括潜水含水层的弹性变形量、非弹性变形量和承压含水层的弹性变形量、非弹性变形量，并由各含水层的压缩量计算总地面沉降量；s24.综合孔隙度、渗透系数和储水率的动态变化因素，将三维非稳定渗流数学模型和地面沉降数学模型通过水头项进行耦合，建立非稳定渗流与地面沉降三维耦合数值模型，并通过有限元程序研制求解。3.根据权利要求2所述的悬挂式止水帷幕深基坑降水方案模拟优化方法，其特征在于：步骤s22中，所述三维非稳定渗流数学模型包括控制方程和定解条件，控制方程如下：其中，(x，y，z)为空间位置坐标；k
xx
、k
yy
、k
zz
为各向异性主方向渗透系数；h为点(x，y，z)在t时刻的水头值；w为源汇项；t为时间；ω为研究域；ss为点(x，y，z)处的贮水率，具体为：s
s
＝ρ
w
g(α nβ
w
)，n为土体的孔隙度，α为多孔介质的体积压缩系数，β
w
为水的体积压缩系数，即为水的密度，g为重力加速度；定解条件包括初始条件和至少一种边界条件，将研究域由中心向外进行水平剖分，确定边界条件，其中，研究域内部边界为降水井点，外部侧向边界为第一类边界条件，底边界为第二类边界条件，顶部为自由面边界，所述初始条件为：h(x,y,z,t)|
t＝0
＝h0(x,y,z,t0)，(x,y,z)∈ω，t＝0；所述边界条件包括第一类边界条件、第二类边界条件和自由面边界条件，其中：第一类边界条件为：
第二类边界条件为：自由面边界条件为：自由面边界条件为：其中，h0(x，y，z，t0)为点(x，y，z)处的初始水头值，h1(x，y，z，t)为点(x，y，z)在t时刻第一类边界上的已知水头值，q(x，y，z，t)为t时刻第二类边界上单位面积的补给量，cos(n，x)、cos(n，y)、cos(n，z)分别为第二类边界外法线方向与坐标轴方向夹角的余弦，μ为饱和差或给水度，表示自由面改变单位高度下从含水层单位截面积上吸收或排出的水量，n
z
为自由面外法线方向向量n＝{n
x
,n
y
,n
z
}的第3个分量，γ1、γ2、γ3分别为第一类边界、第二类边界和自由面边界。4.根据权利要求2所述的悬挂式止水帷幕深基坑降水方案模拟优化方法，其特征在于：步骤s22中，具体包括如下步骤：s221.将研究域离散化，即将研究域剖分成多个单元体，单元体与单元体之间通过单元体顶点一节点相联系；s222.在标准的单元体内建立任意一点的坐标与节点之间相互关系的形函数n
i
，既然满足在节点i，ni＝1，在其他节点ni＝0，在单元体内任意点上的形函数之和等于1；s223.对于不规则的单元体利用等参变换，将单元体坐标的变换与单元体内水头插值函数采用相同数目的节点参数和相同的形函数进行变换，在空间上，使(ε，η，ζ)局部坐标系中形状简单、规则的标准单元体，在(x，y，z)整体坐标系中变换为具有曲线或曲面边界的形状复杂的单元；s224.对三维非稳定渗流数学模型的边值问题的求解转化为对泛函i(h)的极值问题的求解，即泛函i(h)的偏导数为零，其中，泛函i(h)的表达式如下：单元泛函i
e
(h)的表达式如下：s225.对单元泛函i
e
(h)偏导数为零，得到分别对求偏导，并整理得出单元支配表达式为：集成后可得整体有限元支配表达式为：
并转化为时间差分格式，其中，[k]为总体渗透矩阵，由单元渗透矩阵[k]
e
集合而成；[s]为储水矩阵，由各单元储水矩阵[s]
e
迭加而得；[g]为流量补给矩阵，由自由面穿过的各单元流量补给矩阵[g]
e
迭加而得；{f}为节点流量矩阵，由各单元水量矩阵{f}
e
迭加而得；{h}为待求节点水头向量。5.根据权利要求4所述的悬挂式止水帷幕深基坑降水方案模拟优化方法，其特征在于：步骤s225中，若是内部单元，则若为自由面边界单元，则若为有压渗流场，则针对自由面边界，6.根据权利要求4所述的悬挂式止水帷幕深基坑降水方案模拟优化方法，其特征在于：步骤s22中，边界条件与控制方程联立前需对选取的边界调节进行处理，其中，对降水井点采用置大数法进行处理，对自由面边界采用改进复合单元渗透矩阵调整法进行处理以稳定三维非稳定渗流数学模型，其中，将自由面边界的几何区域分为虚单元、实单元和过渡单元并通过迭代求解，其中对于过渡单元采用连续的折现型罚函数法将计算区域拓展至整个单元，使待定自由面边界固化。7.根据权利要求6所述的悬挂式止水帷幕深基坑降水方案模拟优化方法，其特征在于：整体有限元支配表达式的时间差分格式包括显示差分格式、中心差分格式和隐式差分格式，其中，隐式差分格式是无条件稳定的。8.根据权利要求2所述的悬挂式止水帷幕深基坑降水方案模拟优化方法，其特征在于：步骤s23中，潜水含水层的弹性变形量为：δb＝
‑
δh(1
‑
n n
w
)μ
ske
b0＝
‑
δhμ
fe
；潜水含水层的非弹性变形量为：δb
*
＝
‑
δh(1
‑
n n
w
)μ
skv
b0＝
‑
δhμ
fv
；承压含水层的弹性变形量为：δb＝
‑
δhμ
ske
b0＝
‑
δhμ
fe
；承压含水层的非弹性变形量为：δb
*
＝
‑
δhμ
skv
b0＝
‑
δhμ
fv
；其中，μ
ske
为土体骨架成分的弹性储水率，μ
skv
为土体骨架成分的非弹性储水率，μ
fe
为土体骨架成分的弹性储水因子，μ
fv
为土体骨架成分的非弹性储水因子；μ
fe
、μ
fv
为将含水层及包含于其中的弱含水层视为一个整体后，在整个含水层规模上的等效参数；而μ
ske
与μ
skv
的模型为：其中，g为剪切模量，v为泊松比，g为重力常数，c
c
为土体的压缩系数，σ0′
为初始有效应力，e0为初始孔隙比，ρ
w
为水的密度，σ
′
max
为土层的前期固结应力。9.根据权利要求7所述的悬挂式止水帷幕深基坑降水方案模拟优化方法，其特征在于：步骤s24的具体步骤如下：s241.结合孔隙流体在多孔介质中的流动规律及其对多孔介质本身的变形或者强度造成的影响进行土体非线性固结模型改进，即考虑多孔介质内应力场与渗流场之间的相互耦
合作用来处理土体的非线性渗流问题，其中：固结过程中的孔隙度n为：其中，土层厚度为m
i
，由于降水引起的垂直沉降量为s
i
；渗透系数k的变化公式为：其中，n0为初始孔隙度；k0为初始渗透系数；地下水渗流连续性方程的右端项中的储水率μ
s
为：μ
s
＝γ
w
(α nβ
w
)，其中，γ
w
为水的重度；α为土的体积压缩系数；n为孔隙度；β
w
为水的体积压缩系数；根据固结试验曲线e
‑
lgσ
′
的斜率c
c
代入储水率计算公式中，可得储水率μ
s
随e变化的关系为：其中，e0为初始孔隙比；σ
′
为有效应力；s242.将整体有限元支配表达式的时间差分格式中的隐式差分格式中的系数矩阵进行处理，获取对称正定方程组，并结合边界条件，将对称正定方程组的系数矩阵进行预处理，用于减少等价问题的条件数，之后采用共轭梯度法初始化参数后进行迭代求解，获取初始条件下的渗流场水头，并计算任意时刻渗流场水头，用于提高收敛速度和稳定性；s243.采用fortran95语言将上述整个计算过程编制成计算机程序，形成计算机模型，所述计算机模型包括主程序和子程序，所述主程序用于读取文件、输入参数和程序结束时输出计算值，子程序包括形成主元素序号指示矩阵模块、形成整体劲度矩阵模块、形成常数项模块、置大数法修正已知水头模块、方程求解模块、沉降计算模块和参数变化模块。10.根据权利要求9所述的悬挂式止水帷幕深基坑降水方案模拟优化方法，其特征在于：步骤s242中，具体步骤如下：s2421.整体有限元支配表达式的时间差分格式中的隐式差分格式为：对其系数矩阵进行处理，令则隐式差分格式可简化为：[a]{h}
t δt
＝{b}，其中，[a]称为总刚度矩阵；{b}称为常数项；s2422.假设a为对称正定矩阵，则对称正定方程组为ax＝b；s2423.采用对称逐步超松驰迭代法的分裂矩阵作为预处理矩阵m，具体公式为：m＝(2
‑
ω)
‑1(d/ω l)(d/ω)
‑1(d/ω l)
t
，其中，d为a的对角阵；l为a的严格下三角矩阵；ω为松驰因子且0<ω<2；s2424.采用预处理共轭梯度法进行迭代求解，具体迭代格式为：
其中，预处理共轭梯度法每步迭代的主要计算量为(2r
a
6)n次乘法运算，r
a
为对称正定矩阵a的各行非零元素个数的平均值。

技术总结
本发明提供一种悬挂式止水帷幕深基坑降水方案模拟优化方法，包括如下步骤：S1.根据研究域的地质、水文地质特征和含水层间的水力联系，建立深基坑降水水文地质概念模型；S2.根据水文地质概念模型，建立三维非稳定渗流数学模型和含水层土体压缩量的地面沉降数学模型，将三维非稳定渗流数学模型与地面沉降数学模型进行耦合，并通过有限元程序研制、求解；S3.通过抽水试验数据进行参数反演，获取优化的非稳定渗流与地面沉降三维耦合数值模型；S4.进行悬挂式止水帷幕和降水井过滤器不同组合深度的深基坑降水与地面沉降模拟计算，获取最佳深基坑高效降水技术方案。本发明构建的耦合模型满足一定精度和置信度，提高模型的计算速度和稳定性。稳定性。稳定性。


技术研发人员：徐成华 骆祖江 于丹丹 成磊 张也 张兴旺 顾问 周玲玲
受保护的技术使用者：河海大学
技术研发日：2021.07.30
技术公布日：2021/10/23