水泥土复合管桩竖向承载特性现场试验

梁善斋

doi:10.11779/CJGE2021S2065

水泥土复合管桩竖向承载特性现场试验 English Version

梁善斋

湖北交投高路公司,湖北　武汉　430010

详细信息

作者简介:
梁善斋(1964— ),男,高级工程师,主要从事高速公路建设技术与质量管理工作。E-mail: 1526640725@qq.com

中图分类号: TU473
计量
- 文章访问数: 193
- HTML全文浏览量: 24
- PDF下载量: 118
出版历程
- 收稿日期: 2021-08-15
- 网络出版日期: 2022-12-05
- 刊出日期: 2021-10-31

Field tests on vertical bearing capacity of pipe piles in cement-improved soil

LIANG Shan-zhai

HBCI Expressway Development Co., Ltd., Wuhan 430010, China

摘要

摘要: 依托某水泥土复合管桩工现场载荷试验,分析了竖向荷载作用下水泥土复合管桩轴力、桩侧摩阻力和桩身压缩量沿桩身的分布特征。测试结果表明,竖向荷载下水泥土复合管桩荷载-沉降曲线呈缓变性,桩身轴力沿桩身近似呈线性递减,桩端荷载约占总荷载的15%～20%,水泥土复合管桩桩侧摩阻力承担80%～85%的竖向荷载；水泥土与管桩的界面摩阻力较天然土与管桩界面摩阻力提升6～9倍；水泥土复合管桩桩身压缩量占总沉降量的为20%～30%。研究成果可为水泥土复合管桩竖向承载设计计算提供依据。
- 水泥土复合管桩 /
- 轴力 /
- 侧摩阻力 /
- 压缩量
Abstract: In order to study the load transfer of the pipe piles in cement-improved soil under vertical loads, the axial load tests are conducted to analyze the distribution characteristics of axial force, side friction and compression along the pile body. The results show that the load-settlement curves of the pipe piles in cement-improved soil change slowly and gradually. The axial force of the piles decreases linearly along the pile depth. The vertical loads at the bottom of the piles range from 15% to 20% of the vertical loads, and the side friction of the piles in cement-improved soil bears 80% ~ 85% of the vertical loads. The side friction of the piles in cement-improved soil is 6~9 times that of the piles in natural surrounding soil. The ratio of the pile compression to the total settlement is 20%~30%. The research results may provide technical support for the vertical bearing design practice of the piles in cement-improved soil.
- cement-improved pipe pile /
- axial force /
- side friction resistance /
- compression

HTML全文

0. 引言

随着工业化和城市化的加速发展，人们对化石能源的消耗和全球变暖趋势加剧，对新能源的开发和利用变得日益重要^[1]。准确测定岩土材料的导热系数是岩土工程温度场分析和能源桩设计的重要内容，土体导热系数是一个复杂的性质，主要取决于土体的类型、饱和度、粒径分布和颗粒堆积密度等因素^[2]。导热系数是土体热力学特性的重要指标之一，能够直接反映热能在土体介质中的传递特性^[3]。

近年来，很多学者对土体热传导特性进行大量研究，并提出大量的预测模型。Kersten等^[4]采用单根线柱热源径向加热的热传导测试方法对19种不同类型土体的导热系数进行测试，分析了含水率和干密度对导热系数的影响，得出了基于含水率和干密度预测土体导热系数的经验公式；Cote等^[5]和Johansen^[6]对650多个土样的导热系数测试结果进行分析，并在Johansen^[4]提出的归一化模型的基础上建立了广义的归一化模型；Erzin等^[7]基于人工神经网络模型建立了预测土体导热系数的单个预测模型和广义预测模型，得出人工神经网络模型在预测土体导热系数方面的有效性。综上所述，预测土体导热系数的主要是通过经验关系模型、理论计算模型和人工神经网络模型，经验关系模型由于考虑的因素较少，计算精度较低。理论计算模型在计算土体导热系数时，所涉及的参数较多，计算过程复杂，难以在实际工程中应用。人工神经网络模型可以考虑的影响因素较多，在预测土体导热系数时精度较高。

为建立预测土体导热系数预测模型，首先对土体导热系数的主要影响因素进行分析，然后利用人工智能算法对土体导热系数进行预测，对预测结果进行检验，并与传统经验关系模型进行对比，验证人工智能算法的有效性和稳健性。

1. 土体导热特性

土体导热系数是表征土体热物理性质的重要指标，与众多因素密切相关，如土颗粒矿物成分、饱和度、粒径分布和孔隙率等，其中矿物成分是决定土颗粒导热性能的关键因素^[8]。导热系数能够直接反映热能在土体介质中的传递特性。因此，为了研究土体的传热机理，本文主要介绍土体饱和度、矿物成分、颗粒分布和温度等参数对热传导的影响。

为研究饱和度对土体导热系数的影响，根据Tarnawski等^[9-11]对40种土体导热系数实测数据，选择黏土、粉土和粉砂的测试结果来研究饱和度对导热系数的影响，如图 1所示。从图 1可以看出，当土体饱和度增加时，导热系数也随之增加，饱和度增加会提高土体中的含水率，含水率的增加会使土颗粒表面形成水膜，颗粒之间形成水桥，而空气的占比会减小，由于水的导热系数明显大于空气的导热系数，热量通过颗粒表面水膜及其水桥传递增加，使土体导热系数增加^[12]。土体饱和度的增加会提高颗粒间的接触面积，热量通过颗粒与颗粒间的传递增加，使土体的导热系数增大。

图 1 饱和度对导热系数的影响

Figure 1. Influences of degree of saturation on thermal conductivity

下载: 全尺寸图片幻灯片

在影响土体导热系数的因素中，固相的矿物组成成分对固体颗粒导热系数的影响非常重要^[13]，土体的固体颗粒主要由原生矿物和次生矿物组成，而这些矿物的导热系数不同，固体颗粒矿物的组成成分不同会影响土体的导热系数。在土体的矿物组成中，石英的导热系数为7.69 (W·m^-1·K^-1)，而其他矿物的导热系数通常在1.25~4.00 (W·m^-1·K^-1)，石英的导热系数明显高于其他矿物，石英含量高的土体其导热系数也较大^[14]。从图 1可以看出，粉砂的导热系数明显高于黏土的导热系数，因为粉砂的石英含量高于黏土的石英含量，所以石英含量高的土体导热系数较大。因此，土体固体颗粒的导热系数取决于矿物的类型和比例，尤其是石英的含量。

对土体导热系数的影响还包括颗粒分布、颗粒粒径大小、颗粒形状、各向异性和温度等因素。颗粒分布会影响土体的级配，级配的变化会影响颗粒与颗粒间的接触面积，从而影响土体的导热系数。在相同干密度下，颗粒粒径越小，颗粒间的接触点越多，而热量通过固体颗粒内部传递热量要快于颗粒间的传播，在相同情况下，颗粒粒径越小，导热系数越低，从图 1可以看出黏土的颗粒粒径最小，在相同情况下其导热系数最小；颗粒形状会影响颗粒与颗粒间接触的面积，从而影响热量的传播。Macaulay等^[15]对澳大利亚粉砂岩导热系数测试，得出导热系数随层理方向角增加而增加，但层理方向角增加到40°后，其导热系随层理方向角的增加而数趋于平稳。温度的变化会影响分子的热运动，Xu等^[16]得出膨润土的导热系数随温度升高而增大，当温度为90℃时，其导热系数为5℃时的（1.2~1.5）倍；当温度为5~60℃时，导热系数随温度的变化相对较小；当温度高于60℃时，导热系数显著增加。

2. 土体导热系数数据库

本文使用的数据来源于40种加拿大田间土^[9-11]。其中22种为细粒土（含砂量 < 40%），18种为粗粒土（含砂量 > 40%）。根据前述土体导热系数影响因素的分析，选择干密度、饱和度、黏粒含量、粉粒含量和砂粒含量作为预测模型的输入参数，土体导热系数作为模型的输出参数。表 1列出了预测模型输入参数和输出参数的边界值。

表 1 预测模型输入参数的输出参数的边界值

Table 1. Boundary values of input parameters and output parameters of prediction model

类型	参数项	最小值	最大值
输入参数	干密度ρ_d/(g·cm^-3)	0.999	1.728
	饱和度S_r/%	0	100
	黏粒含量c/%	0	42
	粉粒含量s_i/%	0	83.4
	砂粒含量s_a/%	0	100
输出参数	导热系数λ/(W·m^-1·K^-1)	0.1	3.2

下载: 导出CSV

| 显示表格

3. 计算模型

随着人工智能算法在岩土工程领域快速发展，使用这些算法对岩土工程参数进行计算，将节约工程造价及时间。基于这一目的，通过人工神经网络模型（ANN模型）、基于自适应神经网络的模糊推理系统（ANFIS模型）和支持向量机模型（SVM模型）对土体的导热系数λ进行计算。

3.1 人工神经网络模型

人工神经网络模型是基于模仿人类大脑的结构和功能的信息处理系统^[17]。目前已经开发出很多神经网络算法，如Levenberg-Marqurdt算法、共轭梯度算法、贝叶斯正则化算法^[17-19]。神经网络模型由输入层、隐藏层和输出层组成，层与层之间通过激活函数 $f(x)$ 进行转换，激活函数 $f(x)$ 通常为Sigmoid函数、双曲线正切函数（Tanh）和校正线性函数（Relu）。图 2是一个典型的神经网络结构，其中N₁，N₂和N_i是输入参数，Y是输出参数。使用下列过程对人工神经网络建模：

图 2 人工神经网络模型结构

Figure 2. Structure of artificial neural network model

下载: 全尺寸图片幻灯片

将输入数据和输出数据在0到1之间进行归一化。归一化计算式为^[7]

${x_{{\text{norm}}}} = \frac{{x - {x_{\min }}}}{{{x_{\max }} - {x_{\min }}}} \text{，}$

(1)

式中，x_norm为标准化值，x为实际值，x_max为最大值，x_min为最小值。输入参数N_i在隐藏层的Sigmoid激活函数的计算式为

${S_{\text{k}}} = {\text{sigmoid}}\left( {\sum\limits_{i = 1}^N {{w_x}{N_i} + {b_x}} } \right) \text{，}$

(2)

${\text{sigmoid}} = \left( {\frac{1}{{1 + {{\text{e}}^{ - x}}}}} \right) \text{，}$

(3)

式中，S_k，w_x和b_x分别为隐藏层节点的计算值、权重参数和偏差因子。输出层节点的计算式为

${Y_n} = \sum\limits_{j = 1}^M {{w_y}{S_{\text{k}}} + {b_y}} \text{，}$

(4)

式中，w_y为权重参数，b_x为输出层偏差因子。输出层数据的反归一化：

$Y = {Y_n}[{x_{\max }} - {x_{\min }}] + {x_{\min }} 。$

(5)

本文使用的神经网络模型是基于贝叶斯正则化算法建立的预测模型，贝叶斯正则化算法使用岭回归将非线性回归转换为统计问题^[19]。这种方法的优点是在验证过程中不需要使用交叉验证技术，这可以简化模型训练，节约模型的计算时间。

3.2 基于自适应神经网络的模糊推理系统

基于自适应神经网络的模糊推理系统（ANFIS）是将模糊系统和神经网络相结合而构成的网络，采用BP算法和最小二乘估计法的混合算法进行学习，考虑了基于一阶Sugeno模型的两种模糊If-then后置规则。

规则一，如果x=A₁和y=B₁，则

${f_1} = {p_1}x + {q_1}y + {r_1} 。$

(6)

规则二，如果x=A₂和y=B₂，则

${f_2} = {p_2}x + {q_2}y + {r_2} 。$

(7)

式中x和y为输入参数；A_i和B_i为模糊集；f_i为模糊规则指定的模糊区域内的输出值；p_i，q_i和r_i为训练过程中确定的参数。本文使用减法聚类算法训练ANFIS模型，这种算法可以自动估计聚类数目及其位置^[20]。

3.3 支持向量机模型

支持向量机模型（SVM）是一种机械学习方法，其理论是建立一个最优分类超平面，将线性不可分问题转换成线性可分问题^[21]。训练样本为{ ${x_i}$ ， ${y_i}$ }，i=1到m，其中 ${x_i}$ ∈R， ${y_i}$ ∈R，R为实数域。以 $\varphi$ ( ${x_i}$ )表示映射到高维空间中的样本，构造得到的函数表达式为

${y_i} = W\varphi ({x_i}) + b \text{，}$

(8)

式中，W为可调权值向量，b为偏置值。

寻找最优的分类超平面亦即寻找最优的W和b。引入松弛变量 $\xi$ 和 $\xi _{}^*$ ，依据结构风险最小化准则，采用 $\varepsilon$ -SVR模型建立带有约束条件的模型优化函数：

${f_{{\text{optimization}}}} = \mathop {\min }\limits_{{\text{W}},{\text{b}},{\xi _i},\xi _i^*} \left[ {\frac{1}{2}{W^{\text{T}}}W + C\sum\limits_{i = 1}^m {({\xi _i} + \xi _i^*)} } \right] \text{，}$

(9)

约束条件为

$\left.\begin{array}{l} y_i-W^{\mathrm{T}} \varphi\left(x_i\right)+b \leqslant \varepsilon+\xi_i^* ,\\ W^{\mathrm{T}} \varphi\left(x_i\right)+b-y_i \leqslant \varepsilon+\xi_i ,\\ \xi_i, \xi_i^* \geqslant 0 \quad(i=1, \cdots, m), \end{array}\right\}$

(10)

式中，C为惩罚函数， $\varepsilon$ 为损失函数。

使用二次规划优化算法SMO（sequential minimal optimization）构造得到关于x的预测函数 $f(x)$ ：

$f(x) = \sum\limits_{i = 1}^m {(a_i^* - {a_i})} \varphi {({x_i})^{\text{T}}}\varphi ({x_i}) + b \text{，}$

(11)

式中， $f(x)$ 为代表土体导热系数， $a_i^*$ 为 ${a_i}$ 的伴随矩阵，m为样本数目。

3.4 蒙特卡罗模拟

在复杂的岩土工程中使用人工智能算法时，蒙特卡罗方法常用于分析计算模型输入参数存在可变性时的鲁棒性^[22]。图 3显示了蒙特卡罗方法的基本原理。它的原理是随机输入参数来多次重现输出结果^[23]。输入参数的随机波动可以在输出解上传播，这意味着输出的概率密度函数是由所有输入参数的统计信息构成的。输出的概率分布定量分析对于表征模型的稳健性是非常有效的。在本文中，使用蒙特卡罗方法分析预测模型的稳健性，并对结果进行统计分析。

图 3 蒙特卡罗方法原理示意图

Figure 3. Schematic diagram of principle of Monte Carlo method

下载: 全尺寸图片幻灯片

4. 预测模型的建立与性能检验

4.1 预测模型的建立

使用了人工智能算法建立预测土体导热系数的预测模型。根据第一节对土体导热系数影响因素的分析，所有计算模型的输入参数为干密度、饱和度、黏粒含量、粉粒含量和砂粒含量，输出参数为土体的导热系数。将数据分为两个部分，训练部分使用70%的数据，测试部分使用30%的数据。人工神经网络模型隐藏层神经元的数量从1开始增加，通过训练得出最佳隐藏层和隐藏层神经元数。SVM模型是由三次多项式核函数构建的，选择盒约束的正则化常数c为0.112，SVM模型的ε参数为0.0111。预测模型（训练集）对土体导热系数λ的计算结果与实测结果如图 4所示。

图 4 预测模型训练数据集预测值与实测值比较

Figure 4. Comparison between predicted and measured values in training data set of prediction model

下载: 全尺寸图片幻灯片

从图 4可以看出，每个智能计算模型在训练数据集中都准确的计算出土体的导热系数，为了得到预测结果与实测结果的相关性，对预测值与实测值进行线性拟合，拟合结果如图 4所示。预测模型训练效果最好的是ANN模型，相关系数R²=0.9784；然后是SVM模型的训练效果较好，相关系数R²=0.9747；训练效果最差的是ANFIS模型，相关系数R²=0.9361。线性拟合结果表明ANN模型和SVM模型的线性拟合最靠近Y=X等值线，也说明ANN模型和SVM模型的预测精度高。

4.2 预测模型的性能检验

为了量化计算模型的质量，采用了相关系数R²、均方根误差RMSE、平均绝对误差MAE和方差比VAF来检验计算模型的性能^{[17, 24]}。均方根误差RMSE、平均绝对误差MAE和方差比VAF的计算公式如下所示：

${\text{RMSE}} = \sqrt {\frac{1}{{N - M}}\sum\limits_1^N {{{({y_0} - {y_{\text{p}}})}^2}} } \text{，}$

(12)

${\text{MAE}} = \frac{1}{N}\sum\limits_{i = 1}^N {\left| {{y_0} - {y_{\text{p}}}} \right|} \text{，}$

(13)

${\text{VAF}} = \left[ {1 - \frac{{{var} ({y_0} - {y_{\text{p}}})}}{{{var} ({y_0})}}} \right] \times 100 \text{，}$

(14)

式中， ${y_0}$ 为实测值， ${y_{\text{p}}}$ 为预测值，N为样本的编号，M为回归过程的评估参数，var为方差。如果VAF为100%且RMSE为0，则该模型被视为优秀。表 2给出了各预测模型的性能指标值。预测模型对土体导热系数λ的预测值与实测值对比如图 5所示。

表 2 预测模型性能指标

Table 2. Performance indexes of prediction models

模型代号	数据集	R²	RMSE /(W·m^-1·K^-1)	MAE /(W·m^-1·K^-1)	VAF /%
ANN	训练数据集	0.9784	0.0874	0.0673	97.84
ANN	测试集	0.9746	0.0985	0.0701	97.45
ANFIS	训练数据集	0.9361	0.1391	0.1096	94.50
ANFIS	测试集	0.9078	0.1958	0.1262	90.17
SVM	训练数据集	0.9747	0.0941	0.0725	97.47
SVM	测试集	0.9649	0.1147	0.0767	96.48

下载: 导出CSV

| 显示表格

图 5 预测模型测试数据集预测值与实测值比较

Figure 5. Comparison between predicted and measured values of test data set of prediction model

下载: 全尺寸图片幻灯片

从和可以看出，每个预测模型对土体导热系数 $\lambda$ 的预测精度都很高。预测模型预测精度最高的是ANN模型，SVM模型的预测精度略低于ANN模型；预测精度最低的是ANFIS模型，从图 5（b）中可以看出，预测值与实测值偏差较大，数据较离散。线性拟合方程结果表明ANN模型和SVM模型的预测精度高。为了分析计算模型预测值与实测值之间的误差，预测模型的误差分布如图 6所示。

图 6 预测值与实测值误差分布

Figure 6. Error distribution between predicted and measured values

下载: 全尺寸图片幻灯片

从图 6可以看出，在训练数据集和测试数据集，每个智能计算模型的误差分布主要集中在-0.3~0.3 (W·m^-1·K^-1)，每个智能计算模型的计算误差都很小；ANN模型误差分布达到最高峰值，并靠近零；其次分别是SVM模型和ANFIS模型。然而，从表 2可以看出，在测试测试数据集，ANN模型在均方根误差RMSE方面表现最好，均方根误差RMSE=0.0985 (W·m^-1·K^-1)，而SVM模型和ANFIS模型的均方根误差RMSE分别为0.1958 (W·m^-1·K^-1)和0.1147 (W·m^-1·K^-1)；上述分析结果表明，ANN模型在预测土体的导热系数 $\lambda$ 方面的性能最好，ANN模型的权重和偏差基于贝叶斯正则化优化使模型的性能更好，ANFIS模型和SVM模型也很好的预测了土体的导热系数 $\lambda$ 。

4.3 预测模型的稳健性分析

为了评估智能计算模型的稳健性，将用于训练的70%的数据进行随机组合，作为新的输入数据集，生成了300个蒙特卡罗模拟。统计每个预测模型计算值与实测值之间的偏差，每个预测模型获得了300个相关系数R²、均方根误差RMSE和绝对平均误差MAE，计算了关系数R²和均方根误差RMSE的概率分布，以得出智能计算模型在可变输入参数下的性能。蒙特卡罗模拟的相关系数R²、均方根误差RMSE和绝对平均误差MAE概率分布如图 7所示。模拟结果统计如表 3所示。

图 7 蒙特卡罗模拟检验参数概率分布

Figure 7. Probability distribution of parameters by Monte Carlo simulation tests

下载: 全尺寸图片幻灯片

表 3 蒙特卡罗模拟结果统计

Table 3. Statistics of simulated results by Monte Carlo method

检验参数	预测模型	D25	D50	D75	平均值
R²	ANN	0.8926	0.9512	0.9723	0.9132
	ANFIS	0.8731	0.8951	0.9056	0.8761
	SVM	0.9362	0.9459	0.9552	0.9449
RMSE	ANN	0.0838	0.1338	0.2518	0.1898
	ANFIS	0.1935	0.2105	0.2405	0.2285
	SVM	0.1093	0.1333	0.1623	0.1373
MAE	ANN	0.0432	0.0952	0.1832	0.1322
	ANFIS	0.1234	0.1394	0.1614	0.1464
	SVM	0.0839	0.0959	0.1109	0.0989

下载: 导出CSV

| 显示表格

从图 7（a）和表 3可以看出，在相关系数R²概率分布图中可以看出，ANN模型的相关系数R²分布最靠近1，其次分别是SVM模型和ANFIS模型，其中SVM模型的相关系数R²主要分布在0.9至0.97范围内。从图 7（b）和表 3可以看出，在均方根误差RMSE的概率分布图中ANN模型的概率分布最靠近0，其次是SVM模型和ANFIS模型，分布结果跟相关系数R²的概率分布结果相同。从图 7（c）和表 3可以看出，在绝对平均误差MAE的概率分布图中SVM模型的概率分布最靠近0，其次是ANN模型和ANFIS模型。ANN模型、ANFIS模型和SVM模型相关系数R²和均方根误差RMSE的分位数值和平均值如表 3所示。在300次蒙特卡罗模拟中，对相关系数R²、均方根误差RMSE和绝对平均误差MAE的统计研究表明，SVM模型的稳健性最好，其次分别是ANN模型和ANFIS模型。

5. 预测模型性能评价

目前，计算土体导热系数主要是理论计算模型和经验关系模型，理论计算模型计算土体导热系数时，涉及的参数及要求较多，难以有效计算土体导热系数；经验关系模型计算土体导热系数时，具有计算过程简单、易操作等优点，在工程设计中普遍使用此方法。为了评价智能计算模型预测土体导热系数的性能，本文选取Johansen^[6]模型、Cote^[5]模型和Lu^[25]模型计算土体的导热系数。

Johansen模型是Johansen^[6]基于归一化热传导系数概念提出的经验关系模型，该模型能对冻结和未冻结土的导热系数进行计算，经验关系式为

$\lambda = ({\lambda _{{\text{sat}}}} - {\lambda _{{\text{dry}}}}){K_{\text{e}}} + {\lambda _{{\text{dry}}}} \text{，}$

(15)

式中， ${\lambda _{{\text{sat}}}}$ 为饱和土体热导率， ${\lambda _{{\text{dry}}}}$ 为干土热导率，K_e为土体的归一化热导率，当 ${S_{\text{r}}}$ > 0.05时，K_e对于粗土的计算式为

${K_{\text{e}}} = 0.7\log {S_{\text{r}}} + 1 \text{，}$

(16)

式中， ${S_{\text{r}}}$ 为饱和度。当 ${S_{\text{r}}}$ > 0.1时，K_e对于粗土的计算式为

${K_{\text{e}}} = \log {S_{\text{r}}} + 1 \text{，}$

(17)

Johansen通过改进De Vries^[26]热导系数计算模型，提出了干土的导热系数计算公式：

${\lambda _{{\text{dry}}}} = \frac{{0.137{\rho _{\text{d}}} + 64.7}}{{2650 - 0.947{\rho _{\text{d}}}}} \text{，}$

(18)

式中， ${\rho _{\text{d}}}$ 为干密度。在饱和土导热系数计算方面，Sass等^[27]提出的公式被普遍运用，其计算式为

${\lambda _{{\text{sat}}}} = \lambda _{{\text{water}}}^n\lambda _{{\text{solid}}}^{1 - n} \text{，}$

(19)

式中， $\lambda _{{\text{water}}}^{}$ ， $\lambda _{{\text{solid}}}^{}$ 分别为水和土固体颗粒的导热系数，n为孔隙率。

Cote模型是Cote等^[5]根据Johansen^[6]提出的归一化热导系数计算模型的基础上，建立了广义的归一化预测模型：

$\lambda = ({\lambda _{{\text{sat}}}} - {\lambda _{{\text{dry}}}}){K_{\text{e}}} + {\lambda _{{\text{dry}}}} \text{，}$

(20)

式中，K_e为土体的归一化热导率。但K_e和 ${\lambda _{{\text{dry}}}}$ 的计算公式已经被修改，

${K_{\text{e}}}{\text{ = }}\frac{{\kappa {S_{\text{r}}}}}{{1 + (\kappa - 1){S_{\text{r}}}}} \text{，}$

(21)

${\lambda _{{\text{dry}}}} = \chi ({10^{ - \eta n}}) \text{，}$

(22)

式中，κ为经验参数，用于表达土体种类和冻结或未冻结对K_e的影响，χ和η为考虑土的种类和颗粒形状对干土热导系数影响的参数。

Lu模型是Lu等^[25]对12种不同土性的天然土在各种含水率条件进行了大量导热系数测试，通过试验数据拟合，得到了更为简单的干土导热系数和孔隙率的线性关系，并在Johansen提出的归一化热导系数计算模型的基础上，提出了如下经验公式：

$\lambda {\text{ = }}\left[ {\lambda _{{\text{water}}}^n\lambda _{{\text{solid}}}^{1 - n} - (b - an)} \right]\exp [\alpha (1 - S_{\text{r}}^{\alpha - 1.33})] + (b - an) \text{，}$

(23)

式中，a和b为确定干土热传导系数的参数，建议分别取a=0.56和b=0.51，α反映土种类对Kersten变量的影响，针对粗粒土和细粒土分别取0.96和0.27。

智能计算模型对土体预测所得结果λ与通过使用Johansen模型、Cote模型和Lu模型预测得到的λ值进行比较如图 8所示，智能计算模型和经验关系模型的性能检测如表 4所示。

图 8 预测模型与经验关系模型对比

Figure 8. Comparison between prediction model and empirical relationship model

下载: 全尺寸图片幻灯片

表 4 预测模型的性能指标对比

Table 4. Comparison among performance indexes of prediction models

预测模型	R²	RMSE /(W·m^-1·K^-1)	MAE /(W·m^-1·K^-1)	VAF /%
ANN	0.9746	0.0985	0.0701	97.45
ANFIS	0.9078	0.1958	0.1262	90.17
SVM	0.9649	0.1147	0.0767	96.48
Johansen model	0.6464	0.4251	0.3631	75.05
Cote model	0.7744	0.3841	0.2805	78.34
Lu model	0.6432	0.4075	0.3335	71.49

下载: 导出CSV

| 显示表格

从图 8可以看出，ANN模型、ANFIS模型和SVM模型都准确的预测了土体的导热系数，其预测值与实测值偏差较小；Johansen模型、Cote模型和Lu模型都明显低估了土体的导热系数，当λ < 1 (W·m^-1·K^-1)时，Cote模型和Lu模型的预测结果与实际值偏差较小，Johansen模型的预测结果比实测值偏小，但偏离值较小；当λ > 1 (W·m^-1·K^-1)时，Johansen模型、Cote模型和Lu模型的预测值明显低于实测值，且预测数据较离散，预测精度较低。

从表 4可以看出，ANN模型、ANFIS模型和SVM模型的相关系数R² > 0.9，RMSE < 0.2 (W·m^-1·K^-1)，MAE < 0.13 (W·m^-1·K^-1)，VAF > 90%；Johansen模型、Cote模型和Lu模型的相关系数R² < 0.78，RMSE > 0.38 (W·m^-1·K^-1)，MAE > 0.28 (W·m^-1·K^-1)，VAF < 79%；经验关系模型的性能明显低于ANN模型、ANFIS模型和SVM模型。

上述结果分析表明，预测模型对土体进行预测时，预测精度最高的是ANN模型和SVM模型，其次是ANFIS模型，预测精度最低的是经验关系模型。对于土体导热系数进行计算，建议选择ANN模型或SVM模型进行估算。

6. 结论

基于ANN模型、ANFIS模型和SVM模型对土体的导热系数进行计算研究，对预测模型的误差和稳健性进行分析，主要得出以下结论：

（1）对土体导热系数的影响因素进行分析，确定干密度、饱和度、黏粒含量、粉粒含量和砂粒含量作为预测模型的输入参数。

（2）建立预测土体导热系数的ANN模型、ANFIS模型和SVM模型，模型的计算结果表明，3个预测模型的预测结果精度较高。其中ANN模型的预测精度最高，SVM模型的预测精度略低于ANN模型，ANFIS模型的预测精度最低。基于蒙特卡洛模拟对预测模型的稳健性进行分析，预测模型稳健性最好的是SVM模型。

（3）将ANN模型、ANFIS模型和SVM模型的预测结果和传统经验关系模型的预测结果进行对比，ANN模型、ANFIS模型和SVM模型的预测精度明显高于传统经验关系模型；传统经验关系模型仅以含水率、干密度和土体类别作为计算依据，其预测结果与实际值偏差较大，难以满足工程设计要求。

图 1 地层分布及应变计布设示意图

Figure 1. Stratigraphic distribution and layout of strain gauges

下载: 全尺寸图片幻灯片

图 2 试桩V-I的Q-s曲线

Figure 2. Q-s curves of pipe pile V-I

下载: 全尺寸图片幻灯片

图 3 试桩V-I的s-lgt曲线

Figure 3. s-lgt curves of pipe pile V-I

下载: 全尺寸图片幻灯片

图 4 试桩V-I的桩身轴力

Figure 4. Axial forces of pipe pile V-I

下载: 全尺寸图片幻灯片

图 5 试桩V-I桩侧摩阻力

Figure 5. Side friction of pipe pile V-I

下载: 全尺寸图片幻灯片

图 6 试桩V-I桩身压缩量

Figure 6. Compressions of pipe pile V-I

下载: 全尺寸图片幻灯片

图 7 试桩V-I桩端阻力与位移关系曲线

Figure 7. Relationship curve between end resistance and displacement of pile V-I

下载: 全尺寸图片幻灯片

表 1 土层物理力学参数

Table 1 Physical and mechanical parameters of soil layer

地层名称	层厚/m	含水率/%	重度/(kN·m^-3)	黏聚力/kPa	内摩擦角/(°)	压缩模量/MPa
粉质黏土	2	29.7	19.03	29.6	16.6	5.10
淤泥质粉质黏土	2	37.2	18.28	14.2	10.4	3.72
粉土	7	28.2	18.88	10.9	20.6	8.86
粉质黏土	6	30.2	18.74	29.1	16.9	6.82
淤泥质粉质黏土	25	39.1	17.97	14.0	10.2	3.16
粉质黏土	6	30.6	18.92	33.9	15.9	6.76
粉土	2	25.1	19.57	10.2	21.1	11.05
粉砂	6	—	—	—	—	20
黏土	3	36.4	18.52	31.8	14.3	5.24

下载: 导出CSV

表 2 各桩型桩侧摩阻力

Table 2 Side frictions of various piles

土层	复合桩/kPa	管桩/kPa	灌注桩/kPa	复合桩/管桩	复合桩/灌注桩
粉质黏土	125	16.00	14.00	7.48	8.55
淤泥质粉质黏土	65	7.00	6.00	8.63	10.08
粉土	120	18.00	16.00	6.20	6.96

下载: 导出CSV

表 3 试桩极限荷载下沉降压缩量

Table 3 Settlement compressions of piles under load of 2200 kN

试桩编号	桩顶荷载/kN	桩顶沉降量/mm	桩身压缩量/mm	桩端沉降/mm	桩身压缩沉降比
V-I	2200	10.23	2.82	7.41	27.57%
V-II	2200	12.65	3.02	9.63	23.87%
V-III	2200	12.10	3.14	8.96	25.95%

下载: 导出CSV

参考文献(12)

[1]	董平, 陈征宙, 秦然. 混凝土芯水泥土搅拌桩在软土地基中的应用[J]. 岩土工程学报, 2002, 24(2): 204-207. doi: 10.3321/j.issn:1000-4548.2002.02.017 DONG Ping, CHEN Zheng-zhou, QIN Ran. Use of concrete-cored DCM pile in soft ground[J]. Chinese Journal of Geotechnical Engineering, 2002, 24(2): 204-207. (in Chinese) doi: 10.3321/j.issn:1000-4548.2002.02.017
[2]	李俊才, 张永刚, 邓亚光, 等. 管桩水泥土复合桩荷载传递规律研究[J]. 岩石力学与工程学报, 2014, 33(增刊1): 3068-3076. doi: 10.13722/j.cnki.jrme.2014.s1.068 LI Jun-cai, ZHANG Yong-gang, DENG Ya-guang, et al. Load transfer mechanism of composite pile composed of jet-mixing cement and phc pile with core concrete[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(S1): 3068-3076. (in Chinese) doi: 10.13722/j.cnki.jrme.2014.s1.068
[3]	吴迈, 窦远明, 王恩远. 水泥土组合桩荷载传递试验研究[J]. 岩土工程学报, 2004, 26(3): 432-434. doi: 10.3321/j.issn:1000-4548.2004.03.029 WU Mai, DOU Yuan-ming, WANG En-yuan. A study on load transfer mechanism of stiffened DCM pile[J]. Chinese Journal of Geotechnical Engineering, 2004, 26(3): 432-434. (in Chinese) doi: 10.3321/j.issn:1000-4548.2004.03.029
[4]	李立业. 水泥土复合管桩承载特性研究[D]. 南京: 东南大学, 2016. LI Li-ye. Study on the Bearing Capacity of Stiffened DCM Pile[D]. Nanjing: Southeast University, 2016. (in Chinese)
[5]	叶观宝, 蔡永生, 张振. 加芯水泥土桩复合地基桩土应力比计算方法研究[J]. 岩土力学, 2016, 37(3): 672-678. doi: 10.16285/j.rsm.2016.03.008 YE Guan-bao, CAI Yong-sheng, ZHANG Zhen. Research on calculation of pile-soil stress ratio for composite foundation reinforced by stiffened deep mixed piles[J]. Rock and Soil Mechanics, 2016, 37(3): 672-678. (in Chinese) doi: 10.16285/j.rsm.2016.03.008
[6]	JAMSAWANG P, BERGADO D T, VOOTTIPRUEX P. Field behaviour of stiffened deep cement mixing piles[J]. Proceedings of the ICE-Ground Improvement, 2010, 164(1): 33-49.
[7]	钱于军, 许智伟, 邓亚光等. 劲性复合管桩的工程应用与试验分析[J]. 岩土工程学报, 2013, 35(增刊2): 998-1001. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2013S2190.htm QIAN Yu-jun, XU Zhi-wei, DENG Ya-guang, et al.Engineering application and test analysis of strength composite piles[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(S2): 998-1001. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2013S2190.htm
[8]	VOOTTIPRUEX P, SUKSAWAT T, BERGADO D T, et al. Numerical simulations and parametric study of SDCM and DCM piles under full scale axial and lateral loads[J]. Computers and Geotechnics, 2011, 38(3): 318-329. doi: 10.1016/j.compgeo.2010.11.006
[9]	WONGLERT A, JONGPRADIST P. Impact of reinforced core on performance and failure behavior of stiffened deep cement mixing piles[J]. Computers and Geotechnics, 2015, 69: 93-104. doi: 10.1016/j.compgeo.2015.05.003
[10]	WANG A H, ZHANG D W, DENG Y G. A simplified approach for axial response of single precast concrete piles in cement-treated soil[J]. International Journal of Civil Engineering, 2018, 16(10): 1491-1501. doi: 10.1007/s40999-018-0341-9
[11]	WANG A H, ZHANG D W, DENG Y G.. Lateral response of single piles in cement-improved soil: numerical and theoretical investigation[J]. Computers and Geotechnics, 2018, 102: 164-178. doi: 10.1016/j.compgeo.2018.06.014
[12]	王安辉, 章定文, 谢京臣. 软黏土中劲性复合桩水平承载特性p-y 曲线研究[J]. 岩土工程学报, 2020, 42(2): 381-389. https://www.cnki.com.cn/Article/CJFDTOTAL-XBKJ201129003.htm WANG An-hui, ZHANG Ding-wen, XIE Jing-chen. p-y curves for lateral bearing behavior of strength composite piles in soft clay[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(2): 381-389. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-XBKJ201129003.htm

施引文献(12)

期刊类型引用(11)

1.	崔纪飞，柏林，饶平平，康陈俊杰，张锟. 基于人工智能算法的氯盐侵蚀混凝土预测模型. 硅酸盐通报. 2024(02): 439-447 . 百度学术
2.	段文魁，王来发，晁华俊，明锋. 冻结过程中土体导热系数预测模型. 中国农村水利水电. 2024(05): 47-52 . 百度学术
3.	唐少容，殷磊，杨强，柯德秀. 微胶囊相变材料改良粉砂土的导热系数及预测模型. 中国粉体技术. 2024(03): 112-123 . 百度学术
4.	姚兆明，王洵，齐健. 土体导热系数智能方法预测及影响因素敏感性分析. 工程热物理学报. 2024(05): 1440-1449 . 百度学术
5.	邓志兴，谢康，李泰灃，王武斌，郝哲睿，李佳珅. 基于粗颗粒嵌锁点高铁级配碎石振动压实质量控制新方法. 岩土力学. 2024(06): 1835-1849 . 百度学术
6.	李林，左林龙，胡涛涛，宋博恺. 基于孔压静力触探试验测试数据的原位固结系数物理信息神经网络反演方法. 岩土力学. 2024(10): 2889-2899 . 百度学术
7.	王红旗，李栋伟，钟石明，贾志文，王泽成，陈鑫，秦子鹏. 石灰改良红黏土导热系数影响因素及模型预测. 科学技术与工程. 2023(05): 2084-2092 . 百度学术
8.	王才进，武猛，蔡国军，赵泽宁，刘松玉. 基于多元分布模型预测土体热阻系数. 岩石力学与工程学报. 2023(S1): 3674-3686 . 百度学术
9.	王健翔，任瑞琪. 电学等效的稳态平板导热系数测试实验装置. 电子制作. 2023(11): 105-109 . 百度学术
10.	王才进，武猛，杨洋，蔡国军，刘松玉，何欢，常建新. 基于生物地理优化的人工神经网络模型预测软土的固结系数. 岩土力学. 2023(10): 3022-3030 . 百度学术
11.	徐明，康雅晶，马斯斯，张鹤. 基于贝叶斯优化的XGBoost模型预测路基回弹模量. 公路交通科技. 2023(11): 51-60 . 百度学术