黄土地区水工结构湿陷沉降的可靠性修正方法研究

《Geotechnics》：Reliability-Based Correction of Collapse Settlement for Hydraulic Structures on Loess

　　

采用单一平均湿陷应变（collapse strain）来描述空间变异性黄土（loess）地基的沉降，可能导致对水工结构（hydraulic structures）设计安全性的误判，并引发湿陷性沉降损害。本研究考察了一种按地基与结构条件分类的经验性修正方法。修正系数的分类依据来自12个场地的工程地质记录，并以7座灌溉结构的沉降监测数据进行验证。所得复合修正系数（composite factors）在指定地基与结构类别间介于1.30至1.44。针对7个现场案例，修正后的预测结果平均绝对百分比误差（mean absolute percentage error, MAPE）为12.7%，而两套传统计算方案分别为35.2%和32.4%。档案物理模型观测结果独立支持了所提出的荷载重分布（load-redistribution）机制。该修正方法为在初步设计及设计阶段的沉降评估中考虑地基变异性与结构响应提供了一种透明的途径，前提是具备代表性的湿陷应变数据及明确的浸湿（wetting）情景。

论文解读文章

一、研究背景与问题

黄土是一种亚稳定、通常处于非饱和状态的风成土，在天然状态下能够承受较大荷载，但当浸水导致吸力支撑及胶结颗粒接触破坏时，会发生快速体积压缩。这种湿陷行为对水利和灌溉工程尤为重要，因为渗漏、首次充水、接缝渗滤及长期渗流会提供持续性的局部或渐进饱和水源。由此产生的差异沉降可能使刚性结构段发生转动、接缝张开、集中渗流加剧，甚至中断工程运行。

当前研究面临的核心问题在于：采用单一平均湿陷应变无法充分反映天然黄土体变形潜力的空间复杂性。已有研究表明，黄土湿陷受初始组构、应力状态、饱和度、吸力及浸湿路径等多因素共同控制，而非由任何单一指标性质决定。同时，室内试验指标与现场实际沉降之间存在显著差异，现场浸水试验表明沉降具有空间非均匀性，且变形范围可延伸至直接浸湿区域之外。

此外，水工结构本身对应力场产生影响，同时控制浸湿的位置和发展过程，结构-地基相互作用使问题更加复杂。尽管可靠性理论为显式表征这些不确定性提供了框架，但在黄土湿陷性地区相对小型水工结构的常规设计中，仍存在从复杂随机场或耦合水-力分析到实际工程应用之间的空白。

二、研究内容与意义

研究人员基于12个黄土地区水利工程场地的档案工程地质资料，开发了一套将设计湿陷应变分解为三个具有物理意义的组成部分的可靠性修正方法：方向各向异性（directional anisotropy）、天然空间变异性（natural spatial variability）和结构-地基相互作用（structure–foundation interaction）。通过7座已建成水工结构的沉降监测数据进行验证，修正后的预测平均绝对百分比误差从传统方法的35.2%和32.4%降至12.7%。该研究为初步设计和设计阶段的沉降评估提供了一种透明的、能够同时考虑地基变异性和结构响应的实用方法。

三、主要技术方法

本研究采用档案数据分类与统计分析方法。研究人员从12个黄土地区水利工程场地的工程地质档案中提取数据，按岩性连续性、排水与地下水条件、湿陷厚度H_se与预期侧向浸湿宽度B之比进行分类。场地变异性系数以0.95置信概率评估，筛选分辨率为0.01。采用经验查找矩阵方法确定复合修正系数C_ij，并通过可分离性检验（interaction ratio）分析矩阵内部结构。验证部分使用7座已建成水工结构（包括陡坡进水口、渠道分水器、出水口、支座、埋设结构及土堤）的沉降观测数据，将修正方法与传统SNiP计算法和VSN水工结构计算法进行对比。此外，重新分析了档案物理模型试验（几何比尺M_L = 0.02，应力比尺M_σ = 1，模型尺寸约80×25×35 cm），作为结构-地基相互作用机制的独立验证。

四、研究结果

4.1 复合可靠性系数

研究人员获得了25种地基-结构组合的档案复合系数C_ij，取值范围为1.30至1.44。该矩阵表现出地基条件与结构类型之间的交互作用，不可分解为独立的地基乘子和结构乘子。通过交互比率检验，异质浅层地下水地基上独立支座组合的I_ij值为0.936，表明矩阵不满足精确可分离性。若忽略交互作用，该组合的预测值将被高估6.9%。

4.2 沉降验证

7个现场案例的验证结果表明，传统方法对大多数现场沉降存在高估。修正方法在案例1-4和7中更接近实测值，但对案例5和6两种刚性结构仍偏保守。修正方法的误差范围为-1.9%至+27.2%，未出现类似SNiP方法在案例1中+62.0%的极端误差。总体性能指标方面，修正方法的MAPE为12.7%，较SNiP方法的35.2%和VSN方法的32.4%分别降低约63.9%和60.8%；均方根误差（RMSE）从30.7 cm和26.8 cm降至10.1 cm，分别降低约67.1%和62.3%。修正方法保留了+12.1%的正向平均偏差，表明其具有适度的保守性。

4.3 物理模型的结构-地基相互作用证据

档案物理模型试验揭示了湿陷机制对局部应力状态的明显依赖。在荷载较重的进水口段下方，浸湿产生局部竖向变形柱，其边界大致为基脚边缘下的竖向变形面，响应类似冲切或极限状态机制。在荷载较轻的消力池段下方，响应更接近无显著侧向鼓出的一维压缩。模型观测中，湿陷过程中水平-竖向应力比实测值为1.50和1.53，较Coulomb–Mohr被动土压力理论估算值2.04低约25-26%。进水口下方土柱总湿陷变形为2.18 cm，其中0.78 cm归因于侧向形状变化，贡献比例为35.8%。总变形与纯压缩变形之比K = 1.55。消力池起端和末端现场沉降分别为76 cm和93 cm，模型换算至原型尺度后分别为90 cm和100 cm，差异约为18.4%和7.5%。进水口段传统计算平均沉降为168.2 cm，物理模型换算后为185.5 cm，差异约9%。

4.4 设计解读

结构类型在各地基行内对系数产生修正作用，但该差异的幅度在地基类别间发生变化，因此单一独立结构乘子无法复现完整矩阵。项目特定选择应同时考虑浸湿条件和结构体系。

4.5 区域变异性和渐进水-力响应

区域资料综合显示，哈萨克斯坦黄土的湿陷相关特性变异系数约为53-125%，远高于30%的强工程地质异质性判据。天然状态平均内摩擦角为22度，饱和后为21度，变化不大；而平均黏聚力从30 kPa降至12 kPa，降低约60%。这表明湿陷机制主要由水敏胶结的破坏和孔隙组构的重排控制。初始湿陷压力P_se按地貌湿陷类别划分：强湿陷组为10-30 kPa，中等组为20-60 kPa，相对弱湿陷组为50-120 kPa。物理模型记录显示，可见浸湿边界处含水量W约为0.27，饱和度S_r约为0.73；湿陷区内升至W = 0.30、S_r = 0.85。湿陷层沉降在浸湿锋到达后约两天趋于稳定。浸湿自两端向陡坡中心推进，中心部分在外侧进水口区域开始湿陷后约四天才开始沉降，导致结构段间产生转动和相对位移。被渐进浸湿的进水口区下方压力从约154 kPa降至110 kPa再降至30 kPa，相邻较干支撑区压力升至约220 kPa，表明结构将荷载从软化支撑区向较干土体转移。现场观测中，进水口部分外缘沉降约162 cm（自重湿陷参考值为114 cm），水平位移约12 cm，结构部分间竖向相对错位约33-50 cm，接缝张开约6-9度。

五、讨论总结

研究表明，将湿陷应变分解为方向各向异性、空间变异性和结构-地基相互作用三个分量，为工程决策提供了有用的层级结构。各向异性系数γ_a = 1.15具有黄土方向性组构的物理基础，但不应视为普适值。空间变异性系数γ_n = 1.12-1.20解决了高度变异土体以单一均值表征时的信息损失问题。结构相互作用系数γ_s = 1.05-1.08虽然数值较小，但决定了异质地基的哪一部分可能控制响应。物理模型证据支持荷载向较干支撑区转移和重载下侧向变形的机制。

研究局限性包括：开发证据由12个场地、7座监测结构和1个物理模型研究组成，源资料未提供可复现系数矩阵的场地级样本和完整标定记录；矩阵存在可测量的地基-结构交互作用，不能推断独立性或精确乘性可分离性；无法从汇总数据计算分量间的协方差或沉降限值超越概率；模型未包含瞬态吸力、水力传导率、渐进饱和及三维应力重分布。

六、结论翻译

• •
  档案地基-结构矩阵提供了经验性沉降修正；其条目无法通过将报告的分量范围相乘来复现。
• •
  在7个监测案例中，表格化修正降低了预测误差，但保留了正向平均偏差。
• •
  模型观测支持荷载向较干支撑区转移以及较重荷载下侧向变形的机制。
• •
  该矩阵反映了地基-结构交互作用，并支持在其历史条件适用范围内的初步筛选；设计标定需要场地特定数据。
• •
  所提出的修正方法适用于常规黄土水工结构。在极端复杂地形、高含水量特殊黄土层及大跨度水工结构条件下的适用性需要进一步验证。