2025年工程水文学题目及答案.docxVIP

2025年工程水文学题目及答案

题目一：某山区流域设计洪水计算（25分）

某南方山区流域，集水面积F=320km2，主河道长度L=28.5km，河道比降J=4.2‰。流域多年平均年降水量P=1580mm，多年平均流量Q?=18.6m3/s。根据该流域1975-2023年共49年的实测洪水资料（含1998年、2010年两次历史洪水，重现期分别为80年、120年），已整理得年最大洪峰流量系列（单位：m3/s）：580,620,710,850,930,1020,1150,1280,1350,1420（注：为简化计算，仅列出排序后前10项，实际系列含49项）。要求：

（1）采用矩法初估P-Ⅲ型分布参数（Cs=3Cv）；

（2）通过适线法调整参数，确定设计频率P=0.1%（千年一遇）的设计洪峰流量；

（3）说明历史洪水加入对参数估计的影响机制。

答案一

（1）矩法初估参数

首先，计算年最大洪峰流量的均值、变差系数Cv和偏态系数Cs（Cs=3Cv）。

实测系列长度n=49年，历史洪水个数a=2（1998年、2010年），其中1998年洪水流量Q?=1350m3/s（重现期N?=80年），2010年洪水流量Q?=1420m3/s（重现期N?=120年）。历史洪水的调查考证期N取最大重现期120年（N=120年）。

均值计算：

考虑历史洪水的不连续系列，均值公式为：

\[\bar{Q}=\frac{1}{N}\left(\sum_{i=1}^{a}Q_i+\frac{N-a}{n-a}\sum_{j=a+1}^{n}Q_j\right)\]

假设实测系列中除历史洪水外的47项（n-a=47）的和为S=580+620+…+1280=计算得S=49×18.6×365×86400/(320×10?)（注：此处简化，实际应直接求和），但为简化，假设实测系列（含历史洪水）的算术平均为\(\bar{Q}_s=950m3/s\)（根据给定数据估算），则考虑历史洪水后的均值：

\[\bar{Q}=\frac{1}{120}\left(1350+1420+\frac{120-2}{49-2}\times(49\times950-1350-1420)\right)\]

计算得\(\bar{Q}≈1020m3/s\)（具体数值需精确求和，此处为示例）。

变差系数Cv计算：

不连续系列的均方差公式为：

\[\sigma=\sqrt{\frac{1}{N-1}\left[\sum_{i=1}^{a}(Q_i-\bar{Q})^2+\frac{N-a}{n-a}\sum_{j=a+1}^{n}(Q_j-\bar{Q})^2\right]}\]

假设实测系列中各值与均值的离差平方和为\(\sum(Q_j-\bar{Q})^2=4.2×10^6\)，则：

\[\sigma=\sqrt{\frac{1}{119}\left[(1350-1020)^2+(1420-1020)^2+\frac{118}{47}\times4.2×10^6\right]}≈380m3/s\]

因此，Cv=σ/$\bar{Q}$≈380/1020≈0.37，Cs=3Cv=1.11。

（2）适线法调整参数

以矩法初估参数（$\bar{Q}=1020$，Cv=0.37，Cs=1.11）为初始值，绘制理论频率曲线与经验频率点据对比。经验频率采用数学期望公式计算：

对于历史洪水，经验频率\(P_i=\frac{i}{N+1}\)（i=1,2对应2010年、1998年），得P?=1/(120+1)≈0.83%，P?=2/(120+1)≈1.65%；

对于实测洪水（非历史部分），经验频率\(P_j=\frac{a+j}{N+1}\)（j=1~47），最小流量对应的P≈(2+47)/(120+1)≈40.5%。

对比发现，理论曲线在P1%时偏低于经验点据（如2010年洪水P=0.83%对应流量1420，理论曲线计算值为1020×Kp（Kp为模比系数，查P-Ⅲ表，Cs=1.11，P=0.83%时Kp≈1.85），得理论值≈1020×1.85≈1887，与实际1420偏差大），说明Cs可能偏大。调整Cs=2.5Cv=0.925，重新计算Kp（P=0.1%时，Cs=0.925，查得Kp≈3.2），则设计洪峰流量Qp=1020×3.2≈3264m3/s。反复调整后，最终确定$\bar{Q}=1050$，Cv=0.40，Cs=1.0