【学习课件】第四章-降雨径流分析汇流计算

第四章降雨径流分析——汇流计算流域水文学核心课程·从产流到汇流的完整分析框架Contents本章内容概览第四章降雨径流分析——汇流计算01流域汇流基本概念与原理02等流时线法原理与应用03单位线法核心理论与方法04瞬时单位线理论05地下径流汇流与综合案例CHAPTER01流域汇流基本概念与原理从净雨到出口断面流量的物理过程与数学表达Confluence汇流的定义与基本分类汇流是净雨从流域各处出发、经坡面和河网汇集至出口断面的完整物理过程。按空间路径分为坡面汇流与河网汇流两阶段，按物理途径分为地面汇流和地下汇流。流域各点净雨经坡面汇入河网，再由河网流达出口断面，是连接产流与洪水过程线的关键环节净雨在重力作用下沿坡面向最近河槽运动，流速受坡度、糙率和地表覆盖影响显著水流进入河道后沿河网向下游传播，受河道形态、纵比降和河床糙率等水力因素控制经由地表和土壤表层汇入河网，流速快、历时短，是洪峰形成的主要贡献经由地下途径缓慢注入河网，历时长达数天甚至数十天，主要影响退水段和基流补给自然流域中支流汇入干流的航拍场景HYDROLOGY·CHAPTER04流域汇流时间与调蓄作用流域汇流时间的空间分布不均匀性是流域调蓄作用的物理根源，这种天然"水库效应"使出口断面洪水过程线较降雨过程更为平缓展宽。山区流域河道洪水过程—流域调蓄作用下的水位变化01汇流时间定义—雨滴从起始位置汇集至出口断面所经历的时间，受流速和路径长度共同决定02空间差异—最远点汇流可达数小时，出口附近仅数分钟，构成连续的时间分布谱03调蓄本质—降水非单点注入且流速分布不均，入流和出流在时间上存在错位04涨洪阶段—i(t)dt>Q(t)dt，蓄水量持续增加，洪水上涨05落洪阶段—i(t)dt<Q(t)dt，蓄水量逐步释放，洪水消退CONVOLUTIONTHEORY流域出口断面流量的数学表达出口断面流量Q(t)可由净雨强度与汇流曲线的卷积积分精确表达：Q(t)=∫u(τ)·i(t-τ)dτ。这一公式揭示了汇流计算的本质——求解流域的时间响应函数u(t)，所有汇流方法都是对这一核心公式的不同离散化或参数化求解策略。01卷积积分公式Q(t)=∫₀ᵗu(τ)·i(t-τ)dτ，其中u(τ)为汇流曲线，i(t-τ)为净雨强度，物理含义为各时段净雨贡献的线性叠加。该积分形式将连续的降雨输入转化为连续的流量输出，是水文汇流理论的数学基石。02汇流曲线的物理意义u(t)表示单位时间内单位净雨在出口断面形成的流量过程，是流域汇流特性的完整数学刻画。其形状综合反映了流域地形、河网形态及调蓄能力对径流过程的塑造作用。03线性时不变假设经典汇流计算假设流域为线性时不变系统，汇流曲线不随净雨强度和起始时间变化，简化了实际复杂的非线性问题。这一假设使得叠加原理适用，为工程水文计算提供了可行路径。04离散化表达实际计算中将连续积分离散为Q(n)=Σu(k)·i(n-k)，便于利用实测降雨和流量数据进行数值计算。离散形式是计算机编程实现的基础，也是单位线法等实用方法的理论来源。ChapterIV·汇流计算汇流计算方法体系总览流域汇流计算形成了以等流时线法、单位线法和瞬时单位线法为核心的方法体系。三类方法均以卷积公式为理论基础，但在汇流曲线u(t)的求解策略上各有侧重：等流时线法从几何特征出发，单位线法基于实测资料反推，瞬时单位线法则以参数化数学模型替代离散曲线，三者构成了从简化到精确的递进关系。等流时线法从流域地形和几何特征出发，在流域上勾绘等流时线划分汇流面积，概念直观、物理意义清晰。适用于缺乏实测流量资料的小流域，可通过地形图或DEM提取等流时线进行汇流计算。地形驱动单位线法工程实践中应用最广泛的汇流计算方法，通过实测降雨-洪水资料反推单位线，计算流程成熟可靠。包含经验单位线、时段单位线、S曲线时段转换等多种变体，可适应不同时段和精度需求。实测反推瞬时单位线法由Nash于1957年提出，用连续数学函数描述汇流曲线，克服了离散单位线的时段依赖问题。仅需n、K两个参数即可完整描述流域汇流特性，便于参数的地区综合和无资料流域的推广应用。Nash,1957CHAPTER02等流时线法原理与应用从流域几何特征出发刻画汇流时间与面积的空间分布ISOCHRONE·汇流分析基础等流时线的基本概念与勾绘方法等流时线是流域上连接"水质点同时到达出口断面"各点的等值线，相邻等流时线间的面积构成等流时面积ωᵢ，按汇流时间由短到长排列，构成时间-面积分布图。流域地形等高线——等流时线的勾绘基础等流时线定义：在流域地形图上勾绘的一组等值线，线上各点的水质点具有相同的汇流时间，将同时到达出口断面。等流时面积划分：相邻两条等流时线之间的面积称为等流时面积，按汇流时间由短到长依次标记为ω₁(Δt)、ω₂(2Δt)、ω₃(3Δt)……汇流时间确定：根据流域最远点到出口的汇流路径长度和平均流速确定总汇流时间τ，再按Δt等间隔划分n=τ/Δt个时段。勾绘方法：传统方法根据河网密度、坡度变化和分水岭位置手工勾绘；现代方法可利用DEM通过GIS水文分析工具自动提取。ISOCHRONEMETHOD等流时线法的流量计算流量计算本质为净雨与等流时面积的卷积叠加，物理概念清晰但未充分考虑流域调蓄，过程线偏尖瘦，工程常需调蓄校正。流量叠加原理出口断面第i时段流量Qᵢ由各等流时面积上不同时段净雨同时到达出口的贡献叠加组成，体现汇流时间差异对洪峰形成的影响。Qᵢ=Σr·ω卷积表达等价性等流时面积序列可视为离散化的汇流曲线u(t)，流量计算等价于净雨与汇流曲线的离散卷积，数学处理简洁高效。u(t)Discrete方法局限性仅考虑汇流时间差异，未体现流域对洪水的调蓄展平作用，计算结果偏尖偏瘦、洪峰偏高，与实测过程存在偏差。⚠洪峰偏高·过程线偏瘦工程改进措施可引入线性水库等调蓄模型对计算结果进行二次校正，通过调蓄系数优化，有效提高洪水过程线的模拟精度。✓二次校正·精度提升第四章·汇流计算等流时线法计算示例通过等流时面积与分时段净雨的斜线叠加计算，可逐步推求出口断面的完整出流过程线。等流时线法流量计算表（Δt=1h）时段irᵢ×ω₁rᵢ₋₁×ω₂rᵢ₋₂×ω₃rᵢ₋₃×ω₄rᵢ₋₄×ω₅Qᵢ(m³/s)115×20=300----83.3225×20=50015×40=600---305.6310×20=20025×40=100015×60=900--583.340×20=010×40=40025×60=150015×80=1200-861.1500×40=010×60=60025×80=200015×160=24001388.9600010×80=80025×160=40001333.37000010×160=1600444.45块等流时面积（20/40/60/80/160km²）配合3时段净雨（15/25/10mm），洪峰出现在第5时段Q=1388.9m³/sCHAPTER03单位线法核心理论与方法从Sherman经验单位线到时段单位线的完整理论与工程应用Chapter4·RunoffRouting单位线的定义与基本假定Sherman单位线定义为流域上单位时段、单位净雨深（通常10mm）在出口断面形成的地面径流过程线。其理论基础是"倍比假定"与"叠加假定"，二者共同构成线性时不变系统假设。01单位线定义Δt时段内均匀降落10mm净雨所形成的出口断面地面径流过程线q(t)，也称Δt时段单位线。02倍比假定净雨深为单位深度的k倍时，出流过程线各时刻流量均为单位线对应值的k倍，即Q(t)=k·q(t)。03叠加假定连续多时段净雨各自产生的出流过程线，可按时间先后错开后线性叠加，得到总出口流量过程。04线性时不变汇流曲线不随净雨强度和时间变化，简化了实际的物理非线性过程。水文测验站流量观测场景汇流计算·CONFLUENCE单位线法推流计算公式单位线推流公式Qᵢ=Σ(hⱼ/10)·qᵢ₋ⱼ₊₁将多时段净雨的汇流计算转化为单位线的比例缩放与时序错开叠加，是洪水预报和设计洪水计算的核心工具。FORMULAQi=Σj=1m(hj/10)·qi−j+1hj为第j时段净雨深(mm)，q为Δt时段单位线纵坐标(m³/s)推流公式多时段净雨出流等于各时段净雨深按比例缩放单位线纵坐标后的叠加求和，实现净雨到出流的线性转换。h/10·q计算步骤确定单位线与净雨，按比例缩放得分量过程线，按时间列纵向求和得总出流，步骤清晰可操作。三步完成错开规则第j时段净雨产生的出流从该时段才开始，分量过程线需向后错开(j−1)个时段，体现汇流滞后效应。j−1offset适用条件需扣除地下径流取地面净雨，单位线须匹配净雨时段Δt，流域面积一般<500km²，保证计算精度。<500km²汇流计算单位线的推求方法单位线推求是从实测暴雨洪水资料中反推流域汇流特性的逆过程。核心步骤包括选样（单峰孤立洪水）、基流分离、净雨计算和单位线反推。暴雨后河流水位上涨—洪水分析的实际观测对象01洪水选样要求：选择雨型均匀、单峰孤立、洪量集中、历时较短的典型洪水过程，确保净雨时段少、反推结果稳定可靠单峰孤立02基流分离与净雨计算：用斜线割除基流得到地面径流过程线，结合降雨资料和产流模型（蓄满或超渗）计算分时段地面净雨深斜线割除03单时段净雨反推：当净雨仅1个时段时，单位线q(t)=[Q(t)/h]×10，直接由出流过程线按净雨深等比例缩放即可等比缩放04多时段净雨反推：采用最小二乘法使推流误差平方和最小，或用"正向逐时段扣除法"逐段解出单位线各纵坐标值最小二乘HYDROLOGY·UNITHYDROGRAPHS曲线与单位线时段转换S曲线通过时段错开求差，实现任意时段单位线转换，突破对特定时段的依赖。01S曲线定义每隔Δt连续降10mm净雨，在出口断面形成的累积出流过程线S(t)，是无穷多条单位线按时序叠加的结果。时序叠加原理02稳定值Qₛ当t≥流域汇流时间τ后，S(t)趋于常数Qₛ=10F/(3.6Δt)m³/s，物理含义为入流等于出流的平衡状态。水量平衡条件03时段转换公式q'(t)=[S(t)−S(t−Δt')]×Δt/Δt'，本质是用两条错开Δt'的S曲线相减，截取该时段的累积增量。错开求差计算04工程意义将实测Δt时段单位线转换为任意所需时段Δt'的单位线，解决设计洪水计算中的时段匹配问题。灵活时段适配Hydrology·Chapter04单位线推流案例分析多时段净雨叠加使洪峰较单位线峰值滞后约1–2个时段，退水段延长，验证流域调蓄对洪水波形的展平效应。基于Δt=3h时段单位线，配合多时段净雨序列推求洪水过程线，验证流域调蓄展平效应。单位线峰值：第3时段达82m³/s，形态紧凑，反映流域快速集流响应特征。洪峰滞后展平：推求洪峰154.4m³/s滞后至第6时段，滞后约3个时段，退水段明显延长。总洪量估算：总洪量约2050万m³，体现流域调蓄对洪水波形的展平效应。单位线及推求洪水过程线（Δt=3h）单位线峰值第3时段·洪峰第6时段·Δt=3hLIMITATIONS&IMPROVEMENTS单位线法的局限性与改进方向经典单位线法面临三大局限：倍比和叠加假定与实际非线性汇流存在偏差，时段依赖性限制了通用性，大面积流域降雨不均匀性降低了代表性。ISSUE01非线性偏差大洪水流速快、汇流时间短，单位线尖瘦且峰值高；小洪水反之，单位线矮胖——倍比假定无法刻画这种雨强依赖性改进方向：可按洪量大小分组推求多条单位线，或引入非线性水库模型进行校正实际应用中需根据流域特性选择合适方法，综合考虑计算精度与操作复杂度非线性水库校正ISSUE02时段依赖与空间局限不同Δt推求的单位线形状各异，S曲线可转换但增加操作复杂度；Δt选取不当会导致计算误差累积大面积流域（>500km²）降雨时空分布不均匀，单一单位线难以代表全流域汇流特性，需分单元处理分单元汇流计算需合理划分计算单元，并考虑单元间的洪水传播与调蓄作用>500km²分单元CHAPTER04瞬时单位线理论从离散经验曲线到连续参数化模型的汇流理论突破Chapter4·汇流计算Nash瞬时单位线的概念与数学模型Nash瞬时单位线将流域汇流概化为n个串联线性水库的逐级调蓄过程，参数n与K共同描述流域汇流特性，且与时段Δt无关。梯级水库串联示意—Nash模型将流域概化为n个相同线性水库01瞬时单位线定义：Δt→0时单位线的极限形式u(0,t)，表示流域上瞬时降落单位净雨在出口断面形成的流量过程02Nash概化模型：将流域汇流概化为n个相同的串联线性水库（S=KQ），每级水库对入流进行一次调蓄，n级串联后的出流即为瞬时单位线03解析公式参数：n为线性水库个数（正整数或正实数），K为单库蓄泄系数(h)，Γ(n)为Gamma函数值04参数物理意义：n增大→调蓄次数增多→曲线更对称（趋近正态分布）；K增大→单库调蓄能力增强→曲线更宽平、峰值降低参数估计·PARAMETERESTIMATION瞬时单位线参数n、K的确定方法矩法通过实测地面径流过程线的统计矩反推Nash参数：一阶原点矩M₁=nK，二阶中心矩M₂=nK²，联立即可求解。多场洪水分析后可建立n、K与流域特征的地区综合关系式，实现无资料流域的参数估算。01矩法原理利用实测地面径流过程线Q(t)计算一阶原点矩M₁与二阶中心矩M₂，通过统计矩描述汇流过程特征。M₁=∫tQ(t)dt/∫Q(t)dt02参数求解公式由Nash模型理论矩关系M₁=nK和M₂=nK²，联立方程组直接解得两个汇流参数。K=M₂/M₁|n=M₁²/M₂03多场洪水综合分析对多场典型洪水分别求解n和K，取平均值或建立与洪量、雨强的相关关系以提高参数稳定性。多场平均04地区综合移用建立n、K与流域面积F、主河道长L、坡度J的经验公式，使无资料流域可通过地形特征估算汇流参数。n=f(F,L,J)ConversionMethod从瞬时单位线到时段的转换方法瞬时单位线通过积分可得到S曲线，再利用S曲线时段转换公式即可获得任意时段的离散单位线。这一转换流程将参数化的连续函数与工程计算所需的离散时段完美衔接。01S曲线求取将瞬时单位线u(0,t)从0到t积分得S(t)=∫₀ᵗu(0,τ)dτ，数学上等价于Gamma分布的不完全函数，可查表或数值积分实现快速计算。02时段单位线转换qΔt(t)=[S(t)−S(t−Δt)]×(10/Δt)，用两条错开Δt的S曲线之差截取该时段净雨的出流贡献，实现连续到离散的转换。03实用查算工具我国各省水文手册编制了以n、K为参数的S(t)查算表和时段单位线查算表，避免了复杂的Gamma函数积分计算，大幅提升工程应用效率。04完整计算流程确定n和K→查算S(t)→转换为Δt时段单位线q(t)→配合净雨过程推求洪水过程线Q(t)，形成从参数到成果的完整技术链条。HYDROLOGY·UNITHYDROGRAPH经验单位线与瞬时单位线对比分析经验单位线与瞬时单位线各有优劣：前者直接从实测资料反推、无需函数假设但依赖充足资料且受时段限制；后者参数化程度高、便于地区综合但依赖Nash模型概化精度。工程实践中二者通常配合使用。经验单位线EMPIRICALUNITHYDROGRAPH核心优势直接从实测暴雨洪水资料反推，无需预设数学函数形式，能较好地反映流域的真实汇流特性。基于实际观测数据，物理意义明确，结果可信度高。应用局限对资料质量要求高，需选多场典型洪水推求并综合，时段依赖性需通过S曲线转换解决。资料不足地区难以直接应用。瞬时单位线INSTANTANEOUSUNITHYDROGRAPH核心优势仅需n、K两个参数即可完整描述汇流特性，便于建立地区综合公式实现无资料流域的参数移用。数学表达简洁，计算效率高。应用局限Nash串联线性水库为物理概化模型，对河道调蓄显著或形态特殊的流域可能存在系统偏差。参数物理意义相对抽象。Chapter05地下径流汇流与综合案例地面与地下汇流的综合计算及典型洪水推求完整流程HYDROLOGY·SUBSURFACEFLOW地下径流汇流的特点与计算方法地下径流汇流历时长（天至周级）、流速慢，可用线性水库模型精确描述，将总净雨分为地面和地下两部分分别推流后叠加，即可获得综合洪水过程线。01地下汇流特征：经由土壤和含水层缓慢渗入河网，流速极低（cm/h量级），汇流历时可达数天至数十天，主要影响退水段形态02线性水库模型：假设S=Kg·Q，结合水量平衡方程dS/dt=Ig−Q，递推得出流过程03Kg参数确定：利用实测退水段Q(t)=Q₀·e−t/Kg，半对数坐标上近似直线，斜率倒数即为Kg04综合汇流流程：总净雨→产流分割（地面+地下）→分别用单位线和线性水库推流→叠加得总Q(t)地下水渗出汇入溪流，展示地下径流缓慢汇入河网的自然过程BASEFLOWSEPARATION基流分割方法与前处理基流分割是将实测洪水过程线中的地下径流分离、仅保留地面径流的前处理步骤，是单位线法推流和推求的前提。01水平线分割法从起涨点A画水平线至退水段交于B点，AB以下视为基流。该方法操作简便、易于实现，但由于假设基流在洪水期间保持恒定，未考虑地下水对洪水的补给动态变化，因此可能低估实际的基流量，适用于快速估算或资料缺乏的情况。02斜线分割法从起涨点A连斜线至退水段拐点N（地面径流终止点），更符合基流随洪水上涨而逐渐增大的物理实际。该方法考虑了地下水与地表水的交互作用，分割结果较水平线法更为合理，是工程实践中常用的基流分割方法。03退水曲线法利用退水段指数衰减特征Q(t)=Q₀·e−t/Kg，将基流退水曲线向上外推至起涨时刻。该方法基于地下水动力学原理，通过分析历史退水过程确定衰减参数，能够较为准确地刻画基流的变化规律。04水量平衡校核地面径流量Wg应与hg×F相等，偏差过大需调整分割方案或检查产流计算。通过水量平衡原理对分割结果进行验证，是确保基流分割合理性的重要步骤，可有效避免因分割不当导致的洪水预报误差。CASESTUDY综合案例：设计洪水过程线推求以F=200km²流域为例，经产流计算得4时段地面净雨（15/35/20/5mm），配合Δt=3h单位线推求地面径流过程线。推流结果显示洪峰出现在第5时段约286m³/s，峰现时间较净雨峰值滞后约6小时，总洪量约1620万m³，叠加地下径流后构成完整设计洪水过程线。设计洪水过程线推求（F=200km²，Δt=3h）4时段净雨推流：洪峰Q=286m³/s出现在第5时段，较净雨峰值滞后约6h流域参数：F=200km²，Δt=3h，4时段地面净雨序列为15/35/20/5mm洪峰特征：Qp=286m³/s，出现在第5时段，峰现时间较净雨峰值滞后约6小时洪水总量：总洪量约1620万m³，叠加地下径流后构成完整设计洪水过程线Chapter4·降雨径流分析汇流计算常见问题与操作要点汇流计算的精度受净雨时段匹配、洪水选样代表性、基流分割方案和流域空间概化等多个环节影响。关键操作要点包括：净雨与单位线时段必须一致，多场洪水综合分析提高稳定性，基流分割后做水量平衡校核，大流域需分单元处理，单位线总体积必须等于10mm×F的约束条件。时段一致性要求净雨时段划分必须与单位线时段Δt严格一致，不匹配时需通过S曲线转换单位线，切忌混用不同时段的单位线。洪水选样代表性单场洪水推求的单位线受雨型和前期条件影响较大，应选取5场以上典型洪水分别推求后综合分析取平均值。单位线体积校核单位线各纵坐标之和乘以Δt应等于10mm×F/3.6（F为流域面积km²），偏差过大说明基流分割或净雨计算有误。大流域分单元处理面积超过500km²时，需将流域划分为若干子流域，分别推求单位线并逐级河道演进汇合，避免降雨不均匀性的影响。GIUH·UngaugedBasin拓展：地貌单位线与无资料流域应用地貌单位线（GIUH）直接从流域地形地貌特征推导汇流曲线，无需实测洪水资料，结合DEM与GIS技术在无资料流域洪水估算中具有重要应用前景。地貌瞬时单位线（GIUH）由Rodriguez-Iturbe等1979年提出，将流域汇流时间分布与河网拓扑结构（Horton定律）建立定量关系，是水文地貌学的重要理论基础。核心思路利用流域DEM提取河网分级、河长分布和坡度信息，推求水滴从各点到达出口的概率密度函数即为地貌单位线。无资料流域优势仅需地形数据即可估算汇流参数，避免了传统单位线对实测洪水资料的依赖，特别适用于测站稀少或资料缺乏的偏远地区。现代应用发展结合高分辨率DEM、遥感数据和分布式水文模型，在小流域山洪预警、城市雨洪管理及气候变化影响评估中应用日益广泛。CHAPTERREVIEW本章知识体系总结本章围绕"从净雨推求出口断面洪水过程"这一核心问题，构建了等流时线法（几何路径）、经验单位线法（资料驱动）和瞬时单位线法（参数模型）三大方法体系。学习目标是掌握从基流分割、净雨计算到单位线推求、时段转换、推流计算的完整汇流计算链条，并能灵活选择方法解决工程实际问题。核心方法体系等流时线法几何路径从流域几何特征出发，概念直观但未充分考虑调蓄，适用于小流域初步估算经验单位线法资料驱动基于实测资料反推，工程应用最广泛，需掌握推求、验证和S曲线转换技术瞬时单位线法参数模型参数化连续模型，便于地区综合和无资料流域应用，需掌握n、K参数的矩法确定关键计算技能前处理基流·净雨基流分割（水平线/斜线法）与地面净雨计算，确保水量平衡校核通过