先交代一下背景。去年下半年接了个基坑降水的复核项目整体渗入量怎么算都差不太多但侧壁某一小段在雨季强降雨条件下就是反复出问题后来专门对这局部区域重新做了渗入计算才发现整体模型和局部条件根本是两回事。当时我在内部技术笔记里把整套局部渗入计算的流程简称为 loter取的是 local infiltration 的几个关键字母方便在小组里统一叫法。这段时间不断有人找我要这套笔记干脆整理成一篇完整的实操记录。这篇笔记的内容主要是局部渗入和整体入渗到底差在哪、loter 的核心计算思路是什么、参数怎么取、边界怎么切、结果怎么校核以及我踩过的几个坑和排查方法。适合看的读者很明确搞水利、岩土、环境工程的技术人员正在做渗流计算或基坑降水设计的人还有那些用数值软件算完心里不踏实、想用简化模型手工复核一遍的工程师。文章里没有大而全的理论堆砌重点是把“局部”这件事讲透并把能直接上手的计算流程、公式、脚本和参数经验都放进来。1. 为什么要单独算“局部渗入”1.1 整体入渗和局部渗入的工程差异很多朋友一开始不理解场地整体入渗计算明明已经做了为什么还要单独拎出一个局部区域再算一遍我的理解是这样的整体入渗看的是流域或者场地尺度的平均行为它把地形、土壤、植被、降雨这些条件在空间上做了很大的均化。对于一个几百米甚至几公里的场地来说某一小块区域的突变成分会被稀释掉最终反映在均值里几乎看不出来。但工程出问题恰恰都出在局部。比如基坑的某个侧壁段坡脚处的上层滞水来不及排或者雨水沿某个相对高渗透的夹层集中入渗那这一小段的实际渗入量可能是场地均值的好几倍。再比如边坡稳定分析里整体渗入量不大但局部饱和区一旦形成孔隙水压力升高安全系数掉得很快。这就是为什么要单算局部渗入的根本原因工程破坏是局部的计算也必须局部到点。1.2 需要局部渗入计算的典型场景以我实际接触过的项目来看下面几类场景几乎都会碰到局部渗入计算的需求基坑降水和排水设计重点算侧壁在降雨条件下的短期渗入量用来判断集水坑、排水沟的规模够不够。边坡降雨入渗稳定性校核需要得到湿润锋到达深度和局部含水率变化才能评估安全系数的衰减速度。地下管廊和隧道洞口段渗入评估一般关注回填土与原状土交界面的局部渗入异常。农田和绿地局部灌溉入渗设计关注单次灌水量和入渗时间的关系避免地表积水。表格总结一下不同场景的关注点和输出指标场景关注的核心参数需要输出的指标基坑侧壁渗入侧壁面积、降雨强度、侧向补给累计渗入量、单宽渗入量、峰值渗入速率边坡降雨入渗湿润锋深度、含水率变化最不利时刻的含水率分布、入渗深度管廊回填段渗入回填土渗透系数、界面特性渗入量随时间曲线、排水需求局部灌溉入渗土壤入渗能力、灌水定额达到稳定入渗速率的时间、累计入渗量1.3 局部计算的前提边界比公式更重要说句实在话局部渗入计算的公式并不复杂常见模型就那么几个。真正决定计算结果准不准的是你有没有把局部的边界条件搞对。边界切在哪里、哪一边是供水边界、哪一边是排水边界、初始地下水位怎么考虑这些对结果的影响比模型本身大得多。我见过不少人拿一个场地尺度的模型直接去套局部问题算出来当然对不上。loter 这套流程在设计之初就把“边界先行”作为第一条原则。不管用什么模型先花足够时间把计算范围和边界条件定清楚后面才能踏实。这一点在后面第 2 节和第 3 节里会反复强调。2. loter计算局部渗入的核心逻辑2.1 入渗机理的三种常见描述方式先把基础的东西捋一遍。水分进入土壤的过程大致分三个阶段地表供水开始后首先形成地表饱和薄层水分在土水势梯度作用下向下运动湿润锋不断前移接着随着入渗深度增加湿润锋处的基质吸力梯度逐渐减小入渗速率从初始最大值逐步下降最后当湿润锋足够深、重力势主导时入渗速率趋近于饱和导水率。描述这个过程的模型常见有三种Green-Ampt 模型把入渗过程简化为湿润锋推进适合计算短时强降雨条件下的累计入渗量和湿润锋深度参数少、适合手算和表格计算。Philip 模型用级数解近似 Richards 方程适合描述入渗初期的速率衰减过程公式形式简洁。Richards 方程数值解最严格但参数要求高计算量大适合研究性工作或者局部复杂边界条件的高精度模拟。loter 在大多数工程场景里选的是 Green-Ampt 类模型。原因很直接参数少、物理意义明确、用 Excel 或 Python 都能快速实现而且对初始含水率和饱和导水率这两个关键参数的敏感性很直观。需要更高精度时再切换到 Richards 方程的数值工具做校核。2.2 loter把“局部”拆成哪些关键模块loter 的整个计算框架我习惯拆成四个模块边界模块、土体参数模块、空间模块、时间模块。每个模块都是独立的先分别确定再组合起来计算。边界模块主要回答三个问题上边界是什么供水条件降雨强度恒定还是变强度、是否有地表积水侧边界是自由渗流还是不透水边界下边界是地下水位还是深部排水边界。这三个答案直接决定你选用哪种入渗模式。比如上边界供水强度小于土壤入渗能力时受供水强度控制供水强度大于入渗能力且地表有积水时受土壤入渗能力控制。这两种情况下计算方式完全不同。土体参数模块需要确定饱和导水率 Ks、饱和含水率 θs、残余含水率 θr、初始含水率 θi、湿润锋处的平均基质吸力 Ψf。其中 Ks 和 θi 对结果影响最大。空间模块要确定计算剖面的代表性范围比如侧壁竖向剖面还是水平剖面计算深度取多少。时间模块则根据工程关注点选取时间尺度比如短时强降雨看几小时内的累计渗入量长期灌溉渗入看几天甚至几十天的过程。2.3 为什么选择简化模型加辅助工具而不是纯数值软件我知道很多人第一时间会想到直接用 SEEP/W、Hydrus、FEFLOW 这些软件。它们确实功能强大但局部渗入计算这个场景我用下来有几个不方便的地方第一纯数值软件的前处理建模成本高局部小范围调整边界条件后要重新建模效率低第二黑箱属性强参数影响和机理分析不方便快速做敏感性分析第三在项目初期很多参数本身就不准用高精度软件算边际收益很低。loter 的优势在于轻量、透明、快速。输入几组参数几秒钟就能出一整条累计渗入量曲线参数改一改马上能看到变化趋势。等到方案基本定型后再拿数值软件对几个典型工况做精细化复核。两种手段结合起来用效率和质量都能兼顾。3. 实操步骤一套能直接上手的计算流程3.1 第一步把“局部”边界画清楚这一步是 loter 计算的核心前提也是我最强调的一步。开始计算之前先把 CAD 或图纸里要计算的区域单独提出来确定横向范围和纵向深度。横向范围一般取主要影响带比如基坑侧壁取降水井影响范围或裂缝影响范围纵向深度则要覆盖湿润锋可能到达的最大深度再留出安全余量。有一个经验值供参考粉质黏土短时强降雨条件下湿润锋在数小时内一般能下探到 0.51.5 米砂性土可能更快更深取 13 米。如果地下水位很浅计算深度直接取到地下水位附近因为湿润锋到达地下水面后入渗模式会转变为侧向排水。边界范围宁可取大一点也不要因为范围取得太小导致结果失真。另外现场踏勘非常重要。我一般会让现场配合做一个简易试坑或者看一下揭露的土层情况确认局部范围内是否有透镜体、夹层、裂隙等这些特征用文字记录在计算书里作为参数选取和边界概化的依据。3.2 第二步确定土体参数和初始条件参数怎么来直接影响计算结果可靠性。我按优先级排序现场试验数据试坑注水试验、双环入渗试验、取样室内渗透试验。这些数据最可信优先使用。地区经验值当地勘察报告里的经验参数范围或者附近类似工程的反演参数。文献和规范参考值不同土类的典型饱和导水率范围可以作为初估值但必须经过敏感性分析验证。loter 计算需要整理的最小参数清单如下参数符号常用单位获取方式饱和导水率Ksmm/h 或 m/d现场试验或查表饱和含水率θscm³/cm³室内试验残余含水率θrcm³/cm³室内试验或查表初始含水率θicm³/cm³现场含水率测试湿润锋吸力Ψfmm 或 m按土类经验取值没有实测数据时我常用的经验取值区间可以给大家参考粉质黏土 Ks 约 0.010.1 m/d粉土 0.11 m/d细砂 110 m/dθs 一般在 0.350.45θr 取 0.020.05Ψf 对于黏性土取 200400 mm粉土取 100200 mm砂性土取 30100 mm。这些值仅作为初算用正式工程必须有实测或地区经验数据支撑。3.3 第三步输入 loter 并执行计算在边界和参数都确定后就进入计算环节。以下以 Green-Ampt 入渗模型为例展示 loter 中最常用的局部累计入渗量计算方法。Green-Ampt 模型的累计入渗量 I 与时间 t 的隐式关系为[ I K_s t (\theta_s - \theta_i)\Psi_f \ln\left(1 \frac{I}{(\theta_s - \theta_i)\Psi_f}\right) ]这个公式理解起来不复杂累计入渗量 I 等于稳定入渗项 Ks·t 加上湿润锋附近基质吸力所贡献的额外入渗量。其中 ((\theta_s-\theta_i)) 表示土壤可补充的水分容量乘以 Ψf 表示湿润锋处基质吸力对应的当量水深。方程是隐式的因为 I 出现在等式左右两边需要迭代求解。工程中为了方便通常直接求某个时间 t 对应的 I。迭代算法的思路很简单先设一个初始值 I0 Ks·t然后代入公式右边计算新的 I反复迭代直到前后两次差值足够小。这个迭代过程用 Python 脚本实现最方便示例代码在下面。import math def green_ampt_I(t, Ks, theta_s, theta_i, psi_f, tol1e-6, max_iter100): 求解 Green-Ampt 累计入渗量 I (mm) t: 时间 (h) Ks: 饱和导水率 (mm/h) theta_s: 饱和含水率 (-) theta_i: 初始含水率 (-) psi_f: 湿润锋吸力 (mm) d_theta theta_s - theta_i I Ks * t # 初始猜测 for _ in range(max_iter): I_new Ks * t d_theta * psi_f * math.log(1 I / (d_theta * psi_f)) if abs(I_new - I) tol: return I_new I I_new return I # 示例参数粉质黏土 Ks 2.5 # mm/h theta_s 0.43 theta_i 0.15 psi_f 150.0 # mm times [0.5, 1, 2, 4, 6, 12, 24] for t in times: I green_ampt_I(t, Ks, theta_s, theta_i, psi_f) print(ft{t:5.1f} h, I{I:7.2f} mm)示例参数的计算结果大概是这样0.5 小时累计渗入量约 3.36 mm1 小时约 5.45 mm2 小时约 8.63 mm4 小时约 13.2 mm6 小时约 16.9 mm12 小时约 25.3 mm24 小时约 37.0 mm。可以看到前期累计渗入量增长快后期增速趋缓逐渐接近 Ks·t 的线性趋势。这就是前面说的入渗速率从初始高峰逐渐衰减到稳定值的体现。如果你想进一步得到入渗速率 i(t)可以对 I 求导或者直接用数值差分计算。实际工程中我更关心的是给定降雨强度和降雨历时到底有多少水真正进入土壤。这时候要拿降雨强度与土壤入渗能力做对比如果降雨强度小于土壤入渗能力就按降雨强度全量入渗反之按土壤入渗能力控制多余的水形成地表径流。3.4 第四步结果校核与应用计算结果出来后不能直接用要先做三个方向的校核。第一稳态速率校核。当 t 足够大时入渗速率应当趋近于饱和导水率 Ks。如果算出来的最终稳定速率明显偏离 Ks说明迭代或者参数有问题。第二量级复核。和类似场地的实测资料做对比如果差了一个数量级优先查参数尤其是 Ks 和初始含水率。第三数值工具校核。对于复杂边界条件可以用数值软件跑一个局部模型对比累计渗入量曲线。两个模型结果在工程精度范围内能够吻合一般偏差 20% 以内就说明简化计算的可信度足够。工程应用时还要加安全系数。比如用于基坑排水沟设计如果计算得到某个降雨重现期下的侧壁局部累计渗入量我会在此基础上乘以 1.31.5 的安全系数同时考虑排水措施淤堵对排水能力的折减。这是经验习惯不算什么硬性规定但实际项目里确实能避免不少麻烦。4. 实际案例一次基坑侧壁局部渗入计算4.1 案例背景与地质条件拿一个最近做的实际例子来说。某基坑工程坑深 6 米开挖范围内主要是粉质黏土局部夹薄层粉土。勘察报告给出的参数饱和导水率 Ks0.06 m/d换算约 2.5 mm/h饱和含水率 θs0.43初始含水率 θi0.15。因为缺少室内试验的湿润锋吸力数据我按粉质黏土的经验范围取了 Ψf150 mm。雨季设计降雨强度按当地 5 年重现期、30 分钟最大降雨强度 40 mm/h 考虑。要计算的是降雨 6 小时内基坑侧壁某代表性断面的局部累计渗入量以判断现有排水沟容量是否足够。这个案例的“局部”边界是这样切的计算剖面取侧壁单位宽度1 m的竖向截面深度从地表到坑底以下 6 m水平范围取侧壁后方 5 m 范围内的土体。上边界是降雨入渗边界侧边界后方取对称边界坑底边界按自由排水考虑。初始地下水位在坑底以下 3 m因此短期内湿润锋不会到达地下水位可以按无限深处理。4.2 参数整理与试算参数整理如下参数取值备注Ks2.5 mm/h勘察报告换算θs0.43室内试验θi0.15现场含水率测试Ψf150 mm经验值需敏感性分析降雨强度40 mm/h5年重现期30分钟强度计算时长6 h设计暴雨历时先做个简单判断土壤初始入渗能力远大于降雨强度那么前一段时间是按降雨强度控制的。随着入渗进行土壤入渗能力下降降雨强度超过入渗能力的部分会积水并形成径流。因此需要计算土壤入渗能力低于 40 mm/h 的时间和对应的累计入渗量。利用前面的 Python 脚本试算几个时间点对应的土壤入渗速率可以发现大约在 0.5 小时左右入渗速率已经明显下降到 1 小时后入渗能力降到约 10 mm/h 左右降雨强度开始显著大于入渗能力。也就是说真正决定局部渗入总量的窗口期在前 0.51 小时。4.3 结果计算与工程判断把 6 小时的累计入渗量算出来后结合侧壁实际面积就能估算总渗入量。按粉质黏土的经验参数计算6 小时累计入渗量约 16.9 mm。按单位宽度侧壁面积折算如果降雨直接落在水平地面相当于每平米入渗约 16.9 升但侧壁是竖直面实际受雨面积会小一些。最不利情况下还要考虑坡顶积水沿侧壁下流的贡献这时候局部渗入量会比单算竖直面大很多。在排水沟设计上我通常要求排水能力不小于计算最大入渗速率的 1.5 倍同时留出 30 分钟积水调蓄容积。这个案例算下来现有排水沟的过流能力基本够用但前提是坡顶要设置截水沟不能任水流沿侧壁冲刷。后来施工方按这个意见加了坡顶截水措施雨季没有再出现侧壁涌水问题。4.4 这个案例踩过的坑这个项目里我前期也踩过几个坑写出来大家引以为戒。第一个坑初始含水率取错了。第一次试算时我按饱和含水率的 50% 取了 θi0.22结果前期累计渗入量明显偏小。后来核对现场含水率实测数据实际只有 0.15修正后累计渗入量提高了大约 15%。初始含水率越高土壤可吸纳的水量越少算出来的累计渗入量越保守还是越危险一定要结合工程安全侧来判断。第二个坑基坑底边界条件处理。最初打算按不透水边界处理结果渗入量计算偏小。后来分析发现基坑底有排水垫层实际上相当于自由排水边界会增大侧壁向坑内的渗入梯度。改用自由排水边界后结果更接近现场监测值。第三个坑局部夹层的影响。粉质黏土中夹了 20 cm 厚的薄粉土层刚开始没有单独考虑计算结果整体偏小。把这层薄粉土单独提取出来按平行成层加权处理含水率高渗透层的局部入渗贡献后侧壁中下段的局部渗入量明显增加。这个发现直接影响排水孔间距的布置。5. 常见问题与排查技巧实录5.1 参数怎么快速验证计算参数中Ks 和 θi 的影响最大优先验证。现场快速验证 Ks 的常用方法是简易试坑法挖一个直径 30 cm、深 30 cm 的试坑注满水后记录水位下降速率。稳定状态的下降速率近似等于 Ks。这个方法虽然没有双环注水试验精确但作为初算验证完全够用尤其适合没有过多试验经费的项目阶段。θi 的现场测试可以用取土烘干法也可以用地埋式土壤水分传感器实时监测。如果没有设备一个粗略的经验是黏性土在雨后 1 天的含水率大约在田间持水量附近粉质黏土田间持水量一般为 0.250.35。用这个值做初算是可以的但算完后要做敏感性分析。5.2 迭代计算不收敛怎么办用 Green-Ampt 公式迭代求解时偶尔会遇到不收敛的情况。我遇到的原因一般有两个一是初始猜测值离真实解太远二是参数单位混乱导致量级差太大。解决办法很简单把初始猜测值取为 Ks·t同时确保所有参数单位一致比如都用 mm 和 h。另外代码里迭代次数上限不要设太小100 次足够收敛了。如果还不收敛检查公式里 d_theta 是否接近零初始含水率接近饱和含水率时模型本身会退化这时需要换用其他模型。5.3 局部计算和整体模型结果对不上这是被问得最多的问题。局部渗入量算出来和整体模型结果对不上不一定是算错了往往是两个原因一是边界概化不同整体模型用的是网格化边界局部计算用的是简化边界二是空间尺度不同整体模型的网格尺寸远大于局部特征尺寸自然无法反映局部突变。解决方法是做“局部-整体衔接”先用整体模型提取局部区域边界上的水头分布或通量分布再将边界条件细化后带入局部模型计算。反过来也可以把局部计算得到的关键参数如局部渗入系数、等效渗透系数回代到整体模型中重新率定整体模型参数。这样双向校核两个模型的结果就能比较好地对上。5.4 工具使用心得与避坑清单最后整理一个避坑清单都是实际操作层面容易忽略的细节常见问题主要原因解决办法累计入渗量偏大初始含水率取值过低校核现场实测含水率累计入渗量偏小上边界供水条件错误区分降雨控制和入渗控制稳定速率偏离 Ks迭代初值不当或参数单位错误统一单位取 Ks·t 作为初值局部结果和整体模型不符边界概化差异局部-整体双向校核湿润锋到达深度偏浅计算深度预留不足加大计算深度到地下水位或安全范围侧壁渗入被低估夹层高渗带未单独考虑分层处理单独提取夹层工具方面我的建议是先用这套简化流程快速跑通所有工况再做敏感性分析找出影响最大的两三个参数集中精力提高这几个参数的精度。不要一上来把所有参数都订得很细那只会浪费时间。实测下来Ks 和 θi 的精度上去了局部渗入计算结果基本就可靠了。折腾了这么多次局部渗入计算我最深的体会是局部渗入最难的不是公式和软件而是边界切得准不准、参数估得对不对。loter 这套流程虽然简化但把边界和参数这两关把住了再配合现场实测和整体模型校核完全够用。最后再分享一个小技巧做参数敏感性分析时把 Ks 和初始含水率放在最前面调这两个对结果的影响最大如果条件允许优先针对这两个参数做现场试验性价比最高。