大面积堆载预压处理深厚软基工后沉降预测研究

李波; 程文亮; 项存平; 陈军锋; 余建飞

doi:10.11779/CJGE2021S2039

大面积堆载预压处理深厚软基工后沉降预测研究 English Version

1.
长江科学院水利部岩土力学与工程重点实验室,湖北　武汉　430010
2.
珠海市交通勘察设计院有限公司合肥分公司,安徽　合肥　230000
3.
舟山绿色石化基地管委会,浙江　舟山　316000
4.
浙江石油化工有限公司,浙江　舟山　316000

基金项目:

中央级公益性科研院所基本科研业务费项目 CKSF2017012/YT

国家自然科学基金项目 51308067

详细信息

作者简介:
李波(1982— ),男,山东泰安人,博士,主要从事岩土工程和离心模型试验技术方面的研究。E-mail: libo_auliso@126.com

中图分类号: TU472
计量
- 文章访问数: 227
- HTML全文浏览量: 26
- PDF下载量: 144
出版历程
- 收稿日期: 2021-07-31
- 网络出版日期: 2022-12-05
- 刊出日期: 2021-10-31

Prediction of post-construction settlement of soft foundation after large-area surcharge preloading treatment

1.
Key Laboratory of Geotechnical Mechanics and Engineering of Ministry of Water Resources, Changjiang River Scientific Research Institute, Wuhan 430010, China
2.
Zhuhai Traffic Survey and Design Institute Co., Ltd. Hefei Branch, Hefei 230000, China
3.
Zhoushan Green Petrochemical Base Management Committee, Zhoushan 316000, China
4.
Zhejiang Petrochemical Co., Ltd, Zhoushan 316000, China

摘要

摘要: 沉降预测问题一直是软土工程中的难点,而基于实测数据的沉降预测法最为直观有效,但传统的沉降预测方法具有可靠度和预测精度均较低等局限。提出一种新的基于现场实测数据的工后沉降预测模型,利用最小二乘法对模型参数进行求解。结合某石化基地软基处理的现场沉降监测数据,分析地表沉降规律和沉降速率,结果表明新的沉降预测模型预测的沉降曲线与实测曲线非常接近,验证了其合理性。分析表明：较传统的预测方法更为精确,尤其是中后期的预测结果；将整个软土处理沉降期进行分级叠加预测,相对应每一级模型参数也不同,更符合实际软土沉降规律；对比不同地层厚度的长期沉降预测结果,软土工后沉降量与软土厚度密切相关,软土厚度越大沉降越大。
- 堆载预压 /
- 排水固结 /
- 沉降预测 /
- 深厚软基
Abstract: The problem of prediction of soft soil settlement has always been a difficult problem in engineering, but the settlement prediction method based on the measured data is intuitive and effective. The traditional settlement prediction method has the limitation of low reliability and low accuracy, so a new settlement prediction model is proposed. The characteristics of the new settlement prediction model are analyzed, and the model parameters are solved by using the least square method. Based on the site settlement monitoring data of soft foundation treatment in Zhoushan Green Petrochemical Base, the surface subsidence law and settlement rate of monitoring points are analyzed, and the new settlement prediction model is verified theoretically. Finally, the following three conclusions are drawn. Using this model to predict the settlement of soft soil for a long time is more accurate than using some traditional prediction methods. The prediction of the middle and late periods of this method is relatively accurate, and the predicted results of the middle and late periods of actual engineering are more meaningful. The relative error of this method is smaller than that of the three-point method and the hyperbolic method. It is proposed that the whole settlement period of soft soil treatment should be graded and superimposed, and the corresponding parameters of each model are different, which is more in line with the law of actual soft soil settlement. The larger the thickness of soft soil, the greater the settlement.
- surcharge preloading /
- drainage consolidation /
- settlement prediction /
- deep soft foundation

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. Preloading of plastic drain plate

下载: 全尺寸图片幻灯片

图 2 地表沉降随时间关系曲线

Figure 2. Relationship between surface subsidence and time

下载: 全尺寸图片幻灯片

图 3 沉降速率与时间关系

Figure 3. Relationship between settlement rate and time

下载: 全尺寸图片幻灯片

图 4 本文法预测关系简图

Figure 4. Schematic diagram of predicted relation

下载: 全尺寸图片幻灯片

图 5 拟合曲线与实测曲线对比

Figure 5. Comparison between fitted and measured curves

下载: 全尺寸图片幻灯片

图 6 不同软土厚度长期沉降预测

Figure 6. Prediction of long-term settlement of soft soil with different thicknesses

下载: 全尺寸图片幻灯片

表 1 地层的物理力学参数

Table 1 Physical and mechanical parameters of main soil strata

土层名称	含水率/%	重度(kN·m^-3)	压缩模量/MPa	黏聚力/kPa	内摩擦角/(°)
素填土	27.7	20	50	2	36
冲填土	—	19	20	2	28
粉质黏土	46.4	19.6	9.5	23.5	18.5
含砾粉质黏土	24.0	20	31.1	28.1	19.6
凝灰岩	—	20	20000	—	—

下载: 导出CSV

表 2 A3-1预测模型参数

Table 2 Parameters of prediction model

A3-1	c	a	m
第一级	0.17623	0.36635	1.29968
第二级	0.08836	0.37822	1.35654
第三级	0.07781	0.58262	1.44129

下载: 导出CSV

表 3 A3-1本文方法与其他预测方法对比

Table 3 Comparison between proposed method and other prediction methods

实测值	本文方法		三点法		双曲线法
实测值	计算值	相对误差/%	计算值	相对误差/%	计算值	相对误差/%
836.9	820.6	-1.61	695.5	9.0	691.7	9.6
842.9	826.6	-1.60	695.9	10.3	694.1	10.5
848.9	829.2	-1.20	696.0	11.9	695.0	12.0
859.9	831.6	-0.20	696.2	12.5	695.1	12.7
866.9	833.0	0.45	696.3	13.7	695.7	13.8
867.9	834.0	0.45	696.5	14.6	700.0	14.2

下载: 导出CSV

参考文献(9)

[1]	王安明, 李小根, 李彰明, 等. 软土动力排水固结的室内模型试验研究[J]. 岩土力学, 2009, 30(6): 1643-1648. doi: 10.3969/j.issn.1000-7598.2009.06.020 WANG An-ming, LI Xiao-gen, LI Zhang-ming, et al. Indoor model test of dynamic drainage consolidation method of soft soil[J]. Rock and Soil Mechanics, 2009, 30(6): 1643-1648. (in Chinese) doi: 10.3969/j.issn.1000-7598.2009.06.020
[2]	刘吉福. 排水固结法砂垫层厚度需求[J]. 岩土工程学报, 2008, 30(3): 366-371. doi: 10.3321/j.issn:1000-4548.2008.03.010 LIU Ji-fu. Requirements of thickness of sand mats for consolidation method[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(3): 366-371. (in Chinese) doi: 10.3321/j.issn:1000-4548.2008.03.010
[3]	潘林有, 谢新宇. 用曲线拟合的方法预测软土地基沉降[J]. 岩土力学, 2004, 25(7): 1053-1058. doi: 10.3969/j.issn.1000-7598.2004.07.010 PAN Lin-you, XIE Xin-yu. Observational settlement prediction by curve fitting methods[J]. Rock and Soil Mechanics, 2004, 25(7): 1053-1058. (in Chinese) doi: 10.3969/j.issn.1000-7598.2004.07.010
[4]	黄广军. Asaoka法预测软土地基沉降时存在的问题和对策[J]. 岩土力学, 2016, 37(4): 1061-1065, 1074. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201604019.htm HUANG Guang-jun. Problems and their solutions in predicting soft ground settlement based on Asaoka's method[J]. Rock and Soil Mechanics, 2016, 37(4): 1061-1065, 1074. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201604019.htm
[5]	雷学文, 白世伟, 孟庆山. 灰色预测在软土地基沉降分析中的应用[J]. 岩土力学, 2000, 21(2): 145-147. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200002012.htm LEI Xue-wen, BAI Shi-wei, MENG Qing-shan. The application of grey forecasting to analysing soft foundation settlements[J]. Rock and Soil Mechanics, 2000, 21(2): 145-147. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200002012.htm
[6]	钟国强, 王浩, 李莉, 等. 基于SFLA-GRNN模型的基坑地表最大沉降预测[J]. 岩土力学, 2019, 40(2): 792-798, 808. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201902043.htm ZHONG Guo-qiang, WANG Hao, LI Li, et al. Prediction of maximum settlement of foundation pit based on SFLA-GRNN model[J]. Rock and Soil Mechanics, 2019, 40(2): 792-798, 808. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201902043.htm
[7]	汪建斌, 陈忠平. 一种考虑荷载高度的沉降预测新方法[J]. 中外公路, 2003, 23(5): 72-74. https://www.cnki.com.cn/Article/CJFDTOTAL-GWGL200305022.htm WANG Jian-fu, CHEN Zhong-ping. A new method for settlement prediction considering load height[J]. Chinese and Foreign Highway, 2003, 23(5): 72-74. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GWGL200305022.htm
[8]	张仪萍, 曹国强, 李涛, 等. 分级加载条件下沉降预测方法[J]. 中国公路学报, 2005, 18(1): 30-33. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGL200501007.htm ZHANG Yi-ping, CAO Guo-qiang, LI Tao, et al. Prediction method of settlement under multi-stage loading[J]. China Journal of Highway and Transport, 2005, 18(1): 30-33. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGL200501007.htm
[9]	冯怀平, 耿会岭, 韩博文, 等. 非饱和土地区高速铁路路基沉降预测模型[J]. 岩土工程学报, 2017, 39(6): 1089-1095. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201706019.htm FENG Huai-ping, GENG Hui-ling, HAN Bo-wen, et al. Prediction model for settlement of high-speed railway embankment in unsaturated areas[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(6): 1089-1095. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201706019.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 . 百度学术