基于多元回归与有效随机维度(Erd)方法的土体阿太堡界限(Atterberg Limits)及一致性指数(Consistency Index)概率特征化——应用于罗马尼亚苏恰瓦(Suceava)与多夫塔纳(Doftana)场地

《Geotechnics》：Continuous Predictive Modeling of Geotechnical Parameters: A Probabilistic Alternative to Discrete Characteristic Value Estimation

　　

摘要：研究人员提出了一种概率框架，将多元数学回归与有效随机维度(Effective Random Dimension, ERD)方法相结合，用于将易获取的土工指标数据（颗粒级配、天然含水率w、碳酸钙含量CaCO3）转化为符合Eurocode 7(EN 1997-1)及罗马尼亚规范NP 122:2010要求的特征值(Characteristic Value, Xk)。研究基于罗马尼亚摩尔达维亚台地未扰动中新统沉积（Suceava场地）与喀尔巴阡山前褶皱带构造复杂复理层(Flysch)沉积（Doftana场地）的土工试验数据集，建立液限(Liquid Limit, wL)、塑性指数(Plasticity Index, Ip)及一致性指数(Consistency Index, IC)对深度(z)、粘粒含量(Clay Fraction, C)、粉粒(Silt)、砂粒(Sand)及天然含水率(w)的非线性多元回归模型。通过泰勒级数(Taylor Series)一阶展开处理多变量不确定性传播，采用回归残差标准偏差(Se)结合学生化统计系数(kn)按NP 122:2010计算沿深度变化的动态特征值剖面。残余分析验证误差呈零均值高斯(Gaussian)同方差分布，满足Eurocode 7特征值推导前提。Suceava场地wL模型决定系数R2=0.98，Doftana场地Ip模型R2=0.9475。研究证明该方法可替代传统经验安全系数，为由数据驱动的数学模型向可靠性设计规范过渡提供量化依据。

论文解读——《Geotechnics》刊载：基于多元回归与有效随机维度(ERD)方法的土体阿太堡界限及一致性指数概率特征化研究

研究背景与立项依据

传统岩土工程确定设计参数时常依赖全局安全系数作用于算术平均值，或凭经验选取特征值时忽略土体空间变异性与多变量相关性，导致设计偏于保守或不安全。Eurocode 7(EN 1997-1)要求采用基于可靠度的"谨慎估计值"即特征值(Characteristic Value, X_k)，可按受控破坏模式取均值95%置信区间或分布5%分位数，但缺乏从易测物理指标（颗粒组成、天然含水率w、CaCO₃含量）系统推求阿太堡界限——液限(Liquid Limit, w_L)、塑限(Plastic Limit, w_P)导出的塑性指数(Plasticity Index, I_p=w_L-w_P)及一致性指数(Consistency Index, I_C=(w_L-w)/(w_L-w_P))再经概率框架转为合规特征值的方法。罗马尼亚规范NP 122:2010明确允许回归趋势线作动态均值剖面但要求用残差惩罚。本研究针对此空白，结合Suceava（稳定中新统台地沉积、具垂向分层）与Doftana（构造扰动的Curvature Subcarpathians复理层Flysch沉积、高离散度）两典型地质单元数据集，建立多元回归预测模型并引入有效随机维度(Effective Random Dimension, ERD)近似处理多变量不确定性传播，最终输出符合EC7与NP 122:2010沿深度变化的特征值剖面。

主要技术方法概要

研究人员收集Suceava场地68个、Doftana场地50～66个土层试样的基本物理参数：深度z(m)、粘粒Clay(%)、粉粒Silt(%)、砂粒Sand(%)、天然含水率w(%)、碳酸钙CaCO₃(%)及实验室测定的w_L、w_P（计算I_p）、I_C。剔除异常值与缺失记录后做描述性统计。以多元非线性回归建立因变量w_L、I_p、I_C对自变量z、Clay、Silt、Sand、w、CaCO₃及w_L（I_p模型）的函数关系，用决定系数R²与残差正态性评估拟合优度。通过一阶泰勒级数(Taylor Series Expansion)线性化回归函数以量化自变量协方差传播至目标参数的方差。按NP 122:2010用回归预测值X_pred(z)加减k_n·S_e（S_e为回归残差标准差，k_n为考虑样本量的t分布分位数）计算特征值X_k(z)；剪切破坏相关参数（低I_C不利）取下限减号，膨胀/压缩相关参数（高w_L不利）取上限加号。用Spearman秩相关分析变量关联，用直方图与Q-Q图验证残差接近N(0,σ²)且同方差(Homoscedasticity)，确认适用EC7统计推导前提。

研究结果

3.1. Multivariable Regression Models（多元回归模型）

研究人员对Suceava数据集建立I_p对w_L的线性模型：I_p=1.06627·w_L-17.4372（n=68, R²>0.98），该式跨场地验证于Doftana数据仍达R²=0.9475。对于一致性指数I_C，纳入w、1/z（反映上覆压力非线性固结）、粘粒含量C与log项构建幂-指数混合非线性模型（式13），Suceava场R²=0.9055，负1/z项印证随深度增加土体刚度增大。对于液限w_L，以z、粘粒Clay、砂Sand（负系数证实粗粒稀释双电层降低塑性）、CaCO₃（正系数反映泥灰岩互层效应）、w及log(z)、log(Clay)、Sand/w等比值构建对数-线性混合模型（式14），Suceava场R²=0.8408。各模型残差围绕零随机散布、无异方差，证明无系统性偏差。

3.2. Stratigraphic Validation and Textural Analysis（地层验证与质地分析）

Doftana场地预测曲线有效追踪10～25 m内高频地层波动，在w_L与I_p峰值约14 m与27 m处与实测吻合；I_C剖面稳居0.90～1.0（硬～坚硬状态）。Suceava场地预测vs实验室值散点紧密聚于1:1线（R²=0.98），反映台地沉积矿物均一性降低随机噪声，验证模型跨区域适用性。

3.3. Residual Analysis and Error Distribution（残差分析与误差分布）

残差频率直方图与正态分布重叠良好、均值≈0且方差不随预测值增大而变化（同方差），证实模型无偏且满足Gauss-Markov条件，可用均匀统计惩罚系数k_n计算特征值，回归残差标准差S_e作为EC7框架下可靠性计算必要输入。

3.4. Probabilistic Characteristic Values（概率特征值）

以回归线为动态均值，按式(11) X_k(z)=X_pred(z)±k_n·S_e得沿深度特征值包络。Suceava场地w_L特征值区间X_k,inf=58.66%～X_k,sup=70.15%（宽度11.49个百分点安全廊道）；Doftana因地质随机性强区间更宽。大体积变形控制极限状态取上限（不利高w_L引发膨胀/压缩），局部剪切滑裂取下限（不利低I_C降低抗剪）。

3.5. Engineering Implementation and Limit State Selection（工程实施与极限状态选择）

论证特征值取上限或下限取决于主导极限状态——大面积沉降/膨胀潜力评估用偏高特征值，浅层剪切破坏用偏低特征值。数值区间量化了场地非均质性认知边界，取代主观最不利假设。

讨论与结论翻译总结

讨论指出Suceava高R²源于中新统台地矿物均一性，Doftana虽构造复杂仍保持I_p模型R²=0.9475，说明多变量（含z、CaCO₃、w）可过滤单一颗粒级配相关无法处理的地质噪声。负砂粒系数与双电层理论吻合，1/z项数学复现上覆固结历史。NP 122:2010允许回归趋势作动态均值并用残差惩罚确定特征值，本研究予之实证实现。相较文献中双变量相关R²≤0.70～0.82，本框架R²提升15～28个百分点且跨两类地质单元验证。局限含区域校准需求、依赖基本物性数据质量、微薄夹层或有机质未单独建模、假定残差高斯（已验证本数据集成立）及暂只针对指数性质而非直推c'、φ'。

结论(Conclusions)：

(i) 高预测保真度——Suceava台地w_L模型R²=0.98，Doftana复理层I_p模型R²≈0.90（每场地50～68样），较已有基准解释方差提升15～28个百分点且跨异质地质域验证；(ii) 统计有效性——残差零均值高斯同方差（图8）确认模型无偏，满足EC7与NP 122:2010特征值推导Gauss-Markov条件；(iii) 量化安全廊道——ERD方法给出Suceava I_p/w_L相关参数95%特征值区间58.66～70.15%，以数学锚定设计界限替代主观最不利取值，Doftana更宽区间定量反映更高地质不确定性。未来工作建议耦合机器学习并扩展罗马尼亚多岩性地层数据库，将ERD框架延伸至直推粘聚力(c')与内摩擦角(φ')。