基于基本物理指标计算黄土湿陷势的方法

《iScience》：A method for calculating wetting collapse potential of loess based on basic physical indices

　　

黄土广泛分布于中国西北，是黄河流域中游及一带一路廊道基础设施的重要地基土。因其高孔隙、弱胶结、非饱和等结构特征，受降雨入渗、灌溉、地下水位上升及工程供水等影响浸水后可发生显著结构劣化与突缩变形，即湿陷（collapse）。传统评价依赖室内固结仪（odometer）压缩试验，但现场浸水多为渐进、非全饱和过程，直接采用全淹没条件得到的湿陷参数常高估实际变形，导致地基处理偏于保守、成本增加。研究人员以靖边、延安、西安的砂黄土、粉黄土、黏黄土为对象，开展不同初始含水率下的固结仪试验，建立融合含水率、阿太堡界限（Atterberg limits）与干密度的综合物理指标 Z，发现屈服应力（yield stress）PSC 随 Z 增大呈指数衰减，且不同地区黄土经屈服应力与初始孔隙比（void ratio）e0 归一化后压缩曲线趋同。在此基础上，研究人员提出以 Z 及增量 ΔZ 表征初始与浸水状态，给出湿陷起始压力 Pw 经验式，并重建归一化压缩关系，从而实现仅由常规物理指标与施加应力计算湿化塌陷势（wetting collapse potential）Icw。该框架将易测土性参数与变形行为联系起来，可更高效评价黄土湿陷敏感性，支撑黄土场地工程评估。

研究背景方面，黄土（loess）作为风成沉积土，以高孔隙、弱颗粒胶结与部分饱和状态广泛分布于中国黄土高原，在降雨入渗、灌溉渗漏、管道漏水及地下水位缓慢上升等现实浸水条件下易发生结构退化与湿化塌陷（wetting collapse），进而引发地基差异沉降、边坡失稳和道路桥梁等交通基础设施变形。既有规范与工程实践主要依赖室内固结仪（odometer）全淹没浸水压缩试验获取湿陷系数与湿陷势（collapse potential），但现场浸水往往为渐进、非全饱和过程，全淹没试验所得参数会系统性高估实际变形，造成地基处理深度与范围过大、资源利用效率降低。同时，已有基于单参数或少数参数的经验模型难以兼顾含水率、干密度、粒组与阿太堡界限（Atterberg limits，含液限 w_L、塑限 w_P）的耦合作用，且对中国黄土高原不同区域砂黄土、粉黄土、黏黄土的迁移性差。因此，Wang Haiman 等研究人员在《iScience》发表研究，旨在建立由常规物理指标直接估算黄土湿化塌陷势的统一、简化且具跨区适用性的工程评价框架。

关键技术方法方面，研究人员采集陕西靖边（Jingbian）砂黄土、延安（Yan’an）粉黄土、西安（Xi’an）黏黄土原状样，制备不同初始含水率试件，采用九通道并行固结仪施加工况应力序列 12.5、25、50、100、200、400、800、1600、3200 kPa，以 Casagrande 法从 e–log P 曲线判读屈服应力 P_SC，构造综合物理指标 Z=(w?w_P)·ρ_d²/[(w_L?w_P)·ρ_s·ρ_w]，并以 Z₀、Z₁ 及 ΔZ=Z₁?Z₀ 刻画初始与浸水后状态，再用归一化压缩模型反演孔隙比并计算湿化塌陷势 I_cw。

研究结果部分，压缩性与黄土试件小节中，研究人员分析不同地区 e–log P 曲线，发现压缩分缓压缩与快压缩两阶段，转折点为屈服应力 P_SC；同一地区含水率升高则 P_SC 降低，靖边因干密度高而 P_SC 大，延安与西安干密度相近但黏粒含量差异使 P_SC 不同，由此建立 Z 与 P_SC 的指数关系 P_SC=521.68·e^(?1.15·Z)?73.13。黄土湿化塌陷势小节中，研究人员以延安黄土为例，表明相同应力下 I_cw 随浸水含水率增大而增大、随初始含水率增大而减小；低初始含水率时高应力敏感，高初始含水率时敏感区间向低应力转移，且干密度高、砂粒多的靖边黄土塌陷势最低，黏粒多的西安黄土因黏土吸水弱化接触胶结而塌陷势更高。湿化塌陷起始条件小节中，研究人员定义湿化塌陷起始压力 P_w 为 I_cw=0.015 对应的最小竖向应力，发现 P_w 随浸水后 Z₁ 增大而降低；为消除初始状态影响，采用 ΔZ 得到跨区统一经验式 P_w=843.97·e^{[?5.07·(Z₁?Z₀)]}+130.36，并据此划分发生与不发生湿陷的状态空间。湿化塌陷势计算方法小节中，研究人员以 P_SC 归一化应力、以 e₀ 归一化孔隙比，使多地黄土落入唯一带，给出统一压缩方程 e/e₀=exp[?0.11·(P/P_SC)^0.56]，并形成由基本物理指标求 Z₀、Z₁、P_w、P_SC、e_w0、e_w1 再到 I_cw 的计算流程；方法验证小节中，研究人员用文献中 Chen 等、Hu 等独立样本校验，预测与实测湿化塌陷势围绕 1:1 线分布，证明框架具有跨区鲁棒性。

讨论与结论部分，研究人员指出该框架仅需液限、塑限、含水率、干密度、孔隙比等常规参数即可替代大量固结与湿陷试验，适用于黄土场地快速评价，但数据仍限于黄土高原三类典型黄土，未显式纳入应力路径、各向异性、微观结构演化与三维变形，且独立验证样本较少。结论中，研究人员总结：一、黄土压缩具两阶段特征，综合指标 Z 可跨区统一估算屈服应力 P_SC；二、湿化塌陷势受应力水平与初始含水率共同控制，低初始含水率时塌陷集中于高应力，初始含水率升高后敏感区间下移至低应力；三、湿陷起始压力 P_w 随 ΔZ 指数下降，由此可得仅依基本物性判定湿陷是否发生的统一准则；四、经屈服应力与初始孔隙比双重归一化后压缩曲线唯一化，并可由基本物理指标直接算得湿化塌陷势 I_cw，饱和浸水时即传统湿陷势 I_c。该研究为黄土高原基础设施黄土地基评价提供了简化、可迁移的分析工具。