固结沉降量与历时预测中实验室、原位及观测方法的对比分析：北马尔马拉高速公路案例

《International Journal of Geo-Engineering》：A comparative analysis of laboratory, in-situ, and observational methods for predicting the magnitude and duration of consolidation settlement: Northern Marmara highway case history

　　

准确预测层状软土地基中固结沉降的量值与速率仍具挑战性，因其关键输入参数——应力历史（超固结比，OCR）、主固结速率参数（固结系数，c_v）及有效排水距离（d）——往往存在不确定性，且通常分别通过室内试验、原位指标或沉降记录独立获取。本研究以北马尔马拉高速公路（NMH/KMO）路堤的全面现场勘察与沉降监测为依托，对上述方法进行了案例对比分析。先期固结压力及OCR剖面由室内固结曲线获取，并通过静力触探试验（CPT）和标准贯入试验（SPT）相关关系所获得的独立现场证据予以校核；随后将上述OCR趋势应用于一维沉降计算中，包括一个附加的分层场景，即基于CPT土分类将具有排水功能的砂性/粉性夹层从可压缩黏土层厚度中剔除。对于时间效应评估，室内试验与CPT消散试验所得的c_v与基于液限（LL）的NAVFAC经验范围一并进行评估，监测获得的沉降-时间记录则采用多种观测拟合方法进行分析。结果表明：分层沉积物中的沉降预测对OCR解译及排水距离假设高度敏感。当将基于CPT的排水边界纳入原位指标消解读译时，所得沉降-时间预测与实测响应吻合良好。该研究为在层状土体主固结沉降评估中协调实验室、原位及基于监测的输入参数提供了实践依据。

该论文发表于《International Journal of Geo-Engineering》，研究人员针对软土地基上路堤设计的核心难题——即如何可靠预测主固结沉降的终极变形量及其随时间的发展规律——展开了系统的案例研究。在工程实践中，固结沉降-时间响应通常仍采用太沙基（Terzaghi）一维固结理论进行评估，该理论由应力历史（有效先期固结应力σ′_p/超固结比OCR）、固结系数（c_v）及有效排水路径（d）控制。然而，在厚层状剖面中，这些输入参数往往难以精确界定：应力历史随深度变化，排水边界条件可能不明确；尽管设计时常采用单一代表性c_v值，但其选取须谨慎考虑应力与层构依赖性所导致的压缩性与渗透性耦合效应。已有案例表明，采用恒定c_v结合简化排水假设会导致时间速率估计产生偏差，尤其当层状理想化与排水边界不确定时更为显著。室内试验与原位证据在固结沉降评估中常用且不总收敛：室内c_v受试样质量、取样扰动及试验操作影响显著，且实验室与现场推导的c_v可能存在数量级差异；为从原位响应反推c_v，孔压静力触探（CPT）消散试验虽被广泛采用，但其推断结果依赖于所选解译框架及理想化边界条件。在分层土体中，消散试验往往给出深度变化的c_v范围而非唯一常数。此外，基于液限（LL）的经验指导可作为实用筛选参考，但须确保所选趋势与应力历史一致。

本研究以北马尔马拉高速公路（NMH/KMO）路堤工程为案例，针对上述实践空白展开研究。该工程采用建设-运营-移交（BOT）模式实施，研究段位于库尔泰佩-阿基亚哲主路段第6子段，临近穆杜尔努河穿越处的K26桥梁引道区域。路堤顶宽45米，边坡坡度为2H:1V，分级施工历时约5至7天，最大填筑高度5.5米。地基土为第四纪沉积物，主要由黏土、粉土和砂土组成，呈现褐灰-灰绿色、潮湿状态，从中等僵性至僵性，偶含低-高塑性土、细砂夹层、化石贝壳、褐煤及碳酸盐结核。液限呈双峰分布（约40%及60%），室内固结试验获得的c_v介于0.00005至0.00035 cm²/s之间。

研究人员首先建立了应力历史（σ′_p/OCR）的表征方法，采用卡萨格兰德（Casagrande）图解法与?enol应变能法两种方法由室内固结曲线确定先期固结压力，并通过CPT和SPT经验相关关系进行原位验证。随后，将OCR剖面应用于一维固结沉降计算，包括一个"黏土-仅"理想化模型（基于CPT土行为类型SBT分类，将SBT≥6的砂性/粉性排水层从可压缩厚度中剔除）以及完整分层模型。对于固结速率参数c_v，研究人员综合了三类来源：室内固结试验c_v（采用泰勒√t法由t₉₀计算）、CPT消散试验c_v（采用Burns与Mayne的SCE-CSSM解，以U=50%时的t₅₀及T′₅₀=0.030计算）、以及基于NAVFAC LL-c_v相关图的经验范围（Carter与Bentley提出的幂律关系：OC/再压缩区c_v=25000/LL^3.9，NC/正常固结区c_v=80/LL^2.8）。有效排水距离d则基于CPT推断的土类分布确定，将SBT≥6的层位视为潜在排水边界，本案例的代表性排水距离约为3.5-4.0米。

在观测反分析方面，研究人员采用三种方法对监测沉降-时间曲线进行后分析：浅冈（Asaoka）法（基于等时间间隔Δt的S_i-S_i+1回归，S_ult=β₀/(1-β₁)，c_v=-5H_D²lnβ₁/(12Δt)）、郭与朱（Guo and Chu）修正法（S_i^1.667=β₀+β₁S_i+1^1.667）、以及Chung等的最佳拟合法（以泰勒模型替代太沙基解进行全段或分段拟合）。研究结果表明：浅冈法和郭-朱法对采样窗口高度敏感，部分记录出现非物理/非收敛结果（负值或极大值），郭-朱法的非适用率更高；而Chung等的方法在各采样窗口下表现更为稳定，未出现负值或发散情况，但其同时优化S_ult和c_v存在非唯一性问题，需以原位证据（如CPT消散c_v）加以约束。

在c_v估计的对比中（假设d=3.75 m），CPT消散试验的较低值与反算c_v最为接近，高LL（低c_v）趋势亦与反算值吻合较好，表明高LL层位对固结时间行为起控制作用。排水边界假设的验证显示，随机深度CPT消散试验揭示某些层位的c_v至少为相邻层位的5倍，支持了排水通道化的层位解释，反分析确认d≈3.75米的代表性排水距离有效。时间预测对比表明：CPT较低平均值c_v≈0.0040 cm²/s对应U=90%的历时301-393天，与实测窗口（170-400天）基本吻合；而室内c_v≈0.00035 cm²/s预测历时达9.4-12.3年，为实测值的9-12倍，现场代表性极差；NAVFAC的LL≈60%建议值（c_v≈0.0035 cm²/s）给出344-449天，与观测一致，而LL≈40%建议值（c_v≈0.015 cm²/s）则偏短。

沉降量预测方面，将Chung法反算的S_ult与基于OCR剖面的理论计算对比发现：Casagrande法所得OCR导致沉降估计为修正终极值的3-4倍，对曲线形状和人为判断高度敏感；而?enol法OCR结合"黏土-仅"理想化（剔除CPT分类的排水层）使理论沉降与观测值趋于一致（S_?enol,min/S_ult,Chung介于0.78-1.26，偏差约-22%至+26%）。

在讨论部分，研究人员指出层状软土地基中预测可靠性主要取决于模型输入的转换方式而非形式计算方案。观测反分析存在窗口依赖性，短序列或不规则监测条件下，理论约束的最佳拟合法最为一致。沉降量受OCR解译及排水层处理方式控制，排除排水层可改善理论预测。采用"应力历史确定→CPT层构转化为显式排水模型→c_v以消散试验及LL指导约束→监测校准"的集成框架，可获得更可靠的设计基础。

研究结论概括如下：本案例历史中，OCR解译对沉降预测变异性影响显著，回归支持的应变能法（?enol法）提供了更一致的OCR剖面基础，且排除排水层后结论改善；将CPT层构转化为排水模型是关键建模选择，CPT分类推断的排水边界和代表性排水距离为分层沉积物中控制排水路径的确定提供了实用基础；对于推断的代表性排水距离范围（d≈3.5-4.0 m），CPT消散较低平均值（c_v≈0.0040 cm²/s）的U=90%历时约301-393天，与观测窗口（≈170-400天）吻合，而室内值预测9.4-12.3年不匹配；NAVFAC再压缩线指导中，高LL趋势（LL≈60%，c_v≈0.0035 cm²/s）与观测一致，低LL趋势则偏短；浅冈法和郭-朱法对窗口敏感且部分记录非收敛，Chung等的方法更为稳定；设计可靠性在以下集成框架下得以提升：优先确定OCR（?enol法配合原位校核）、将CPT层构转化为显式排水模型确定代表性排水距离、以消散试验和LL指导约束c_v、辅以监测进行校准。