SWMMH水文模型核心：子汇水面积划分与下渗方法选择指南
简介：这份资源是《SWMM用户参考手册》第I卷「水文学」的汉化PDF，面向城市雨水管理、环境工程与城市规划领域的从业者及学生，用于理解SWMMH动态降雨-径流模拟模型的水文计算原理。SWMMH为美国EPA SWMM的汉化版本，可处理单一事件或长期连续模拟，其径流部分在子汇水面积上接受降水并产生径流与污染物负荷，演算部分则通过管道、渠道、蓄水处理设施、水泵和调节器输送径流，并跟踪各管渠的流量、水深与水质。手册由Lewis Rossman与Wayne C. Huber编写，同济大学李树平教授于2016年完成翻译，内容涵盖对象模型、过程模型、气象学等章节，并附译者与原著声明。资源包共1个PDF文件，约2.93MB，结构完整、便于检索。目前已有298人学习下载，适合需要系统掌握SWMM水文建模、对照官方手册查漏补缺的读者参考。
1. 为什么一场雨能算准，靠的不是公式多，而是子汇水面积切得对
很多人第一次翻 SWMM 手册，注意力全被 Horton、Green-Ampt 这些下渗公式吸走，觉得模型准不准就看公式选得高不高级。但真正跑过几年城市排水模拟的人会告诉你一个反直觉的结论：决定一场降雨模拟成败的，往往不是下渗公式，而是你把研究区域切成了什么样的子汇水面积。SWMMH 作为 EPA SWMM 的汉化版本，它的水文模型核心思路是把城市地表拆成一组子汇水面积，每个子汇水面积独立接受降水、产生径流和污染物负荷，再通过管道、渠道、蓄水设施、水泵和调节器组成的系统把径流输送出去。这个「先切块、再汇流」的架构，才是整套水文模型能 work 的底层骨架。搞懂它为什么这么切、切完怎么算，比背十个公式都有用。
2. 子汇水面积：SWMM 水文模型的第一性拆分
2.1 为什么不能把整片城区当成一块地来算
城市地表最麻烦的地方在于它极不均匀。同一场雨落在沥青路面、屋顶、草坪和裸土上，产流时间可能差十几分钟，径流系数能差三倍以上。如果直接把整个汇水区当成一个均质斜面，用一套参数去算，模型会强行把快慢不同的产流过程抹平，峰值流量和峰现时间都会失真。
SWMMH 的解法是空间离散：把整个研究区域划分成若干子汇水面积，每个子汇水面积内部假定水文响应相对均一，拥有自己的面积、坡度、宽度、不渗透比例、曼宁粗糙系数和洼地蓄水深度。降水在每个子汇水面积上独立转化为径流，然后通过指定的排放口进入管渠系统。这样做的代价是参数变多，收益是每个子区域可以按自己的地表特征来标定，模拟结果才有可能对上实测流量过程线。
从对象模型的角度看，SWMMH 把整个系统抽象成几个核心对象：子汇水面积、节点、管渠、蓄水/处理设施、水泵和调节器。子汇水面积是水文计算的起点，节点和管渠是水力演算的载体。理解这个对象划分，是理解后面所有公式的前提。
2.2 控制方程组：非线性水库到底在算什么
SWMMH 的地表径流计算用的是非线性水库模型。它的核心思想是：把每个子汇水面积想象成一个有出口的水池，降水往里灌，下渗和洼地蓄水先扣掉，剩下的水从出口往外流。水量平衡写成微分形式：
d = d0 + (i - f - e) * dt - q * dt其中 d 是地表蓄水深度，d0 是上一时刻的蓄水深度，i 是降雨强度，f 是下渗率，e 是蒸发率，q 是出流率，dt 是时间步长。出流率 q 用曼宁公式表示：
q = W * (d - d_s)^(5/3) * S^(1/2) / nW 是子汇水面积宽度，d_s 是洼地蓄水深度，S 是坡度，n 是曼宁粗糙系数。这个方程的关键在于：出流不是线性的，它和蓄水深度的 5/3 次方成正比。这意味着水深稍微增加，出流会迅速增大，这正是城市地表产流「来得快、去得也快」的数学表达。
常见的实现思路是：在每个时间步长内，先根据降雨和下渗更新蓄水深度，再用曼宁公式算出流，然后用数值积分推进到下一时刻。由于方程是非线性的，SWMMH 内部会用迭代或近似方法求解，具体迭代策略在手册的计算方法一节有说明。参数估计时，宽度 W 是最敏感的参数之一，它直接决定出流快慢，手册专门用一节讨论子汇水面积宽度的估算。
2.3 离散化到什么程度才算够
子汇水面积划分不是越细越好。切得太粗，空间差异被抹平；切得太细，参数标定工作量爆炸，而且很多子区域根本没有实测数据支撑。手册里给出的原则是：子汇水面积的划分应该让每个子区域内部的地表特征相对一致，同时子区域的数量要控制在计算资源和标定能力允许的范围内。
一个实用的判断标准是：如果两个相邻地块的坡度、不渗透比例和地表类型差异不大，就没必要分开；如果差异明显，比如一个是停车场、一个是绿化带，就应该分开。子汇水面积的宽度参数尤其要注意，它不是一个可以直接测量的几何量，而是从面积和漫流路径长度推算出来的等效值。手册里对宽度的估算给了专门的方法，因为它是影响径流过程线形状最直接的参数。
离散化之后，每个子汇水面积还要指定一个排放口，径流通过这个排放口进入管渠系统。排放口可以是一个节点，也可以是另一个子汇水面积。这个连接关系决定了径流在空间上的汇流路径，如果连错了，模拟出来的流量过程线会完全对不上。
3. 下渗：四种方法背后的取舍逻辑
3.1 Horton 方法：为什么指数衰减能描述下渗
Horton 方法的核心假设是：下渗能力从初始最大值按指数规律衰减到稳定最小值。控制方程写成：
f = f_min + (f_max - f_min) * exp(-k * t)f 是当前下渗能力，f_max 是初始下渗能力，f_min 是稳定下渗能力，k 是衰减常数，t 是时间。这个公式的物理直觉是：土壤刚开始很干，吸水快；随着土壤饱和，吸水能力逐渐下降，最后稳定在一个由土壤渗透性决定的下限。
Horton 方法的关键细节在于下渗能力的恢复。如果降雨停了，土壤会慢慢恢复吸水能力，SWMMH 用一个恢复系数来描述这个过程。手册里专门有一节讲下渗能力的恢复，因为如果忽略恢复，连续模拟时第二场雨的产流量会明显偏大。
参数估计方面，f_max、f_min 和 k 三个参数需要根据土壤类型和前期条件来定。手册给出了不同土壤类型的参考值范围，但实际标定时往往需要根据实测流量过程线反推。Horton 方法的优点是参数少、计算快，缺点是它没有明确的物理基础，参数和土壤物理性质之间的对应关系比较模糊。
3.2 Green-Ampt 方法：从物理出发但参数更难搞
Green-Ampt 方法基于达西定律，假设土壤中存在一个 sharp 的湿润锋，湿润锋以上是饱和区，以下是初始含水量区。下渗率写成：
f = K_s * (1 + (psi * delta_theta) / F)K_s 是饱和导水率，psi 是湿润锋处的吸力，delta_theta 是饱和含水量和初始含水量之差，F 是累计下渗量。这个公式的物理意义很清晰：下渗率由重力项和毛细吸力项共同决定，随着累计下渗量增加，毛细吸力项逐渐减小，下渗率趋近于 K_s。
Green-Ampt 方法的优点是物理基础扎实，参数有明确的物理意义。缺点是参数多，而且 psi 和 delta_theta 在实际中很难直接测量。手册里对 Green-Ampt 的参数估计给了详细讨论，包括如何从土壤质地数据推算这些参数。对于需要长期连续模拟的项目，Green-Ampt 通常比 Horton 更稳定，因为它对前期条件的描述更物理。
3.3 曲线数方法：经验公式的适用边界
曲线数方法（SCS Curve Number）是一个纯经验方法，它用一个综合参数 CN 来描述下渗和产流的关系。CN 值越大，产流越多。这个方法的优点是参数少、使用简单，很多地区的降雨径流手册里都能查到 CN 值。缺点是它原本是为农业流域设计的，在城市区域使用时需要做修正，而且它不能描述下渗的时间变化过程，只适合估算总径流量。
SWMMH 里同时提供了 Horton、Horton 修正、Green-Ampt 和曲线数四种方法，这不是为了让你随便选，而是因为不同方法适用于不同场景。单一事件模拟、参数资料少的时候，Horton 或曲线数可能够用；长期连续模拟、需要描述土壤水分动态的时候，Green-Ampt 更合适。手册里对每种方法都给了数值示例，这些示例是理解方法行为差异的最好材料。
4. 地下水与融雪：被低估的两个水文过程
4.1 地下水模块的通量项为什么这么设计
SWMMH 的地下水模块把地下系统分成上下两层，上层接受地表下渗，下层接受上层穿透，两层都可以向上蒸发蒸腾，地下水最终通过深层穿透或侧向排放进入管渠系统。控制方程是一个水量平衡方程，描述两层土壤含水量的变化。
通量项包括：地表下渗 fI、上层蒸发蒸腾 fEU、下层蒸发蒸腾 fEL、上层到下层穿透 fU、下层深层穿透 fL、地下水排放 fG。每个通量项都有对应的计算公式，手册里逐项给出了表达式和参数含义。用户还可以自定义通量公式，这给了一些特殊场景下的灵活性。
参数估计的难点在于土壤湿度限值和穿透参数。土壤湿度限值决定了蒸发蒸腾和穿透什么时候发生，穿透参数决定了水从上层到下层、从下层到更深层的速率。手册里对每个参数都给了估算方法和参考范围，但实际标定往往需要结合地下水位实测数据。
4.2 融雪：能量平衡比水量平衡更麻烦
融雪模块比降雨径流复杂得多，因为它涉及能量平衡。雪深的变化不仅取决于降雪量和融雪量，还取决于积雪的热交换过程。SWMMH 的融雪模型考虑了融化方程和积雪热交换，融化方程描述雪在气温和太阳辐射作用下的融化速率，热交换描述积雪与大气之间的热量交换。
面积亏损是融雪模块的一个特殊概念：当积雪开始融化时，不是整个雪面同时融化，而是从部分区域开始，逐渐扩大。这个过程中，积雪覆盖面积逐渐减小，但雪深可能还没降到零。SWMMH 用面积亏损曲线来描述这个过程，它影响净径流的计算。
融雪模块的参数估计需要气温、风速、太阳辐射等气象数据，这些数据在连续模拟中尤其重要。手册里对气象数据的来源和处理给了专门章节，包括用户提供的数据、来自政府机构的数据和降雨接口文件。对于有降雪的城市，融雪模块是全年连续模拟不可缺少的部分。
5. 避坑：跑 SWMM 水文模块时最容易翻车的几个地方
5.1 子汇水面积宽度拍脑袋填，峰值流量对不上
现象：模拟出来的流量过程线峰值明显偏高或偏低，峰现时间也偏早或偏晚，但下渗参数怎么调都改善不大。
原因：子汇水面积宽度 W 是影响出流速度最直接的参数，但很多人直接拿面积除以某个长度随便填一个值。宽度偏大，出流快，峰值偏高、峰现偏早；宽度偏小则相反。
解决：手册里对宽度估算给了专门方法，核心思路是从漫流路径长度反推等效宽度。实际操作中，可以用实测流量过程线的峰值和峰现时间来反推宽度，把它作为标定参数而不是固定输入。
5.2 下渗能力恢复被忽略，连续模拟第二场雨产流偏大
现象：单一事件模拟结果还行，但连续模拟时第二场、第三场雨的径流量明显偏大，好像土壤一直没恢复吸水能力。
原因：Horton 和 Green-Ampt 方法都有下渗能力恢复机制，但如果恢复系数设置不当，或者模拟时间步长太大，恢复过程会被跳过。
解决：检查下渗能力恢复相关参数，确保恢复系数合理。连续模拟时，时间步长要足够小，让恢复过程能被捕捉到。手册里对每种下渗方法的恢复机制都有说明，对照检查。
5.3 子汇水面积排放口连错，汇流路径完全乱套
现象：单个子汇水面积的产流量看起来正常，但系统总出口的流量过程线形状完全不对，或者某些管段流量为零。
原因：子汇水面积的排放口指定错误，径流没有进入预期的管渠系统，或者进入了错误的节点。
解决：检查每个子汇水面积的排放口设置，确保汇流路径和实际管网拓扑一致。可以用手册里的对象模型说明来对照，确认子汇水面积、节点和管渠之间的连接关系。
5.4 蒸发参数照搬默认值，长期模拟水量平衡偏移
现象：长期连续模拟时，总径流量和实测值有系统性偏差，而且偏差随时间累积。
原因：蒸发是长期水量平衡的重要输出项，如果蒸发参数照搬默认值，没有根据当地气象条件调整，累积误差会很明显。
解决：根据当地蒸发数据标定蒸发参数。手册里对蒸发数据的来源和处理有专门讨论，连续模拟时尤其要注意蒸发的时间变化。
5.5 时间步长选太大，非线性水库求解不稳定
现象：模拟结果出现数值震荡，流量过程线有锯齿，或者某些时间步长内蓄水深度出现负值。
原因：非线性水库的出流和蓄水深度是 5/3 次方关系，时间步长太大时，数值积分会不稳定。
解决：减小时间步长，尤其是降雨强度变化剧烈的时间段。手册里对时间步长考虑有专门一节，建议根据子汇水面积的大小和降雨特性来选择。
6. 从手册到模型：怎么验证你真的理解了这套水文机制
理解 SWMMH 水文模型的最好方式，不是把手册从头读到尾，而是拿一个简单的算例，手动走一遍计算流程。手册里每一章都有数值示例，这些示例是验证理解的最好材料。你可以先按手册的步骤手算一个子汇水面积的产流过程，再用模型跑同样的输入，对比两者结果。如果对不上，回去查是哪一步的参数或公式理解错了。
另一个验证方法是做参数敏感性分析。手册里专门有一节讨论参数敏感性，你可以自己动手改一个参数，看输出怎么变。比如把子汇水面积宽度加倍，看峰值流量和峰现时间怎么变；把下渗参数改一改，看产流总量怎么变。这种动手验证比单纯读公式有效得多。
从设计取舍的角度看，SWMMH 的水文模型在物理性和实用性之间做了很多折中。非线性水库模型不是最精确的地表径流模型，但它计算快、参数少、能抓住城市产流的主要特征。四种下渗方法并存，是因为没有一种方法能在所有场景下都最优。子汇水面积离散化是一个工程判断，没有标准答案，需要根据研究目的和数据条件来定。
我自己的习惯是：每次拿到一个新项目，先不急着调参数，而是先把子汇水面积划分和汇流路径检查一遍。血泪经验告诉我，大部分模拟结果对不上的问题，根源都在空间离散和连接关系上，而不是公式选得不对。从那以后我每次建新模型都强制走一遍「划分检查—连接检查—参数标定」的流程，顺序不乱。希望帮到你。
本文还有配套的精品资源，点击获取