水文预报课程设计说明书(自动保存的)3

水文预报课程设计计算说明书学院：三峡大学水利与环境学院专业：水文与水资源工程指导老师：刘冀老师班级：20111502学号：2011150209姓名：杨芝辰2014年01月14日一、设计目的1、流域水文模型的用途：洪水预报方案是现代实时洪水预报调度系统的核心部分，是提高预报精度和增长预见期的关键技术。对水资源可持续利用：流域水文模型是水资源评价、开发、利用和管理的理论基础。对水环境和生态系统保护：流域水文模型是构建面污染模型和生态评价模型的主要平台。流域水文模型还是分析研究气候变化和人类活动对洪水、水资源和水环境影响的有效工具。本次课程设计的目的是通过一个具体的降雨～径流预报方案的制作，使学生了解生产单位对预报任务的要求。2、通过课程设计，要求掌握如下内容：（1）流域综合退水曲线、地下水退水曲线的制作以及次洪分割方法；（2）熟悉降雨～径流相关图编制的完整过程；（3）新安江两水源模型结构及产流参数率定方法；（4）流域经验单位线的推求方法；（5）洪水预报方案精度评定方法；（6）利用预报方案进行实时洪水预报方法；（7）利用马斯京根分段连续演算法进行长河段洪水演进预报。二、设计基本资料某流域集水面积为XXkm2（具体数值见给定数据），干流河长约273km，地貌属于低山丘陵区，以低山分布较广，相对高差多在500km左右，最上游尚有部分原始森林，以松、衫、竹类为主，灌木次之。山涧盆地种植农作物，植被总覆盖面积占全流域的80%以上，土壤多为红色沙壤土，间有亚黏土层。山坡残积坡坡积一般厚度为1~2米，最厚约为4~5米，在山麓坡积层与基岩接触面上，或河流下切至接触部位时，常有泉水流出，涌水量较大每秒可达数升。本流域气候温湿，年降雨量在1700毫米左右，汛期为4~9月，降雨量占全年的73%左右，冬季有降雪，量很少。地下水位较高，且季节变幅小，因此，一般情况下，土壤含水量较大。洪水流量过程线极不对称，涨洪历时仅数小时至十多小时，而退水历时竟达十余日，洪水末端的底水明显抬高，说明洪水期潜水、壤中流补给十分丰富。说明本流域以蓄满产流为主。根据流域的自然地理情况和气候条件，以及洪水流量过程线的分析，可知流域产流规律符合湿润地区的蓄满产流特征。采用降雨径流相关图制作蓄满产流方案，用二层蒸发模型计算蒸发，水源划分暂考虑为两水源。三、课程设计资料1、一场历史洪水的流量过程，相应的流域平均时段雨量（时段长为3h）2、洪水起涨前的若干日降雨量，日蒸散发能力资料（历史洪水与实时预报洪水）3、XX场洪水的退水过程4、部分场次洪水降雨、径流特征值成果表5、干流河段的马斯京根参数及分段数（3段）四、产流方案的编制1、次径流深的计算（1）制作地下水退水曲线并计算地下水消退系数Cg。由给定的多场洪水的退水数据，每场洪水相同流量数值即为退水段，从数据中摘取出来进行绘图，采用最小二乘法公式计算。数据（表一）及退水曲线图（图1）如下：表一地下水退水数据退水过程24.821.117.915.31311.19.486.85.84.94.23.632.62.21.91.61.41.21图1地下水退水过程曲线,1,212,21niiiniiQQCgQ（式一）（2）制作流域平均退水曲线并计算流域平均Cg。多场洪水计算退水平均值即可绘图，如图2，同理应用式一即可算出流域平均Cg。05101520253001234567891011121314151617181920212223流量时段地下水退水过程图2流域平均退水曲线（3）制作蓄泄关系曲线。假定一个始退流量值Q0=500，已知2）中得到的流域平均Cg值，由式二即可求出相邻时段流量，并由式三计算各流量对应的R0值，如表二所示，即可绘出Q0-R0关系图（图3）1ttQQCg（式二）0.020.040.060.080.0100.0120.0140.0160.00123456789101112131415161718192021222324252627282930313233流量时段流域平均退水曲线表二Q047.6257.0668.3781.9298.17117.63140.95168.90202.39242.52290.60348.22417.27500.00R01.902.292.763.323.994.805.776.938.329.9811.9714.3617.2220.6611023.6()/2nniiQQRtQA（式三）（4）次洪分割及次洪径流深的计算。由给定的降雨及流量过程即可绘出降雨流量过程曲线图（图四），已知可采用下述式四确定0.0005.00010.00015.00020.00025.0000.0100.0200.0300.0400.0500.0600.0RtQ0图3Q0—Rt关系曲线Q0—Rt关系曲线地面径流终止点B，式中N为洪峰至B点的时距，单位为日，可以为小数；F为流域面积，km20.20.8NF（式四）得到14N，可以在图中确定B点，径流分割如图4所示。图4降雨流量过程曲线2、流域蓄水量W0的计算。由公式五以及给定的已知的5日前Pa=60计算每日的流域初始蓄水量，表格如下（表三）。,1,()attatPkPP1/kEMWM（式五）其中，Pa为第t日降雨，Pa,t为第t日流域初始蓄水量，EM为流域蒸散发能力，WM为流域土壤含水量。0102030405060708090100020040060080010001200140001234567891011121314151617181920212223降雨流量(m3/s)02004006008001000120001234567891011121314151617181920212223流量Rs终止点B起始点A表三流域初始蓄水量日降雨量（mm）5.1020.36.62.5流域初始蓄水量60567172693建立降雨径流相关图及误差评定1）b值的试算。采用单变量求解的方法，假设一个初值，根据蓄满产流模型计算不同初始含水量W0下的不同降雨条件下的径流深，参考给定的已知的洪水的P-R数据，取最优的参数b值，计算出最终的降雨径流相关数据表（表四），经试算b=0.4，最后绘出降雨径流相关图（图5）。表四降雨径流相关计算表W00102030405060708090100110120pa0.0010.1720.7031.6843.2055.4068.5182.9799.77126.00126.00126.00126.0000.000.000.000.000.000.000.000.000.000.000.000.000.00100.160.500.871.291.762.302.963.815.1110.0010.0010.0010.00200.661.352.122.983.955.096.498.3411.3420.0020.0020.0020.00301.502.583.775.116.638.4410.6813.7620.0030.0030.0030.0030.00402.734.215.867.719.8512.4115.6720.4930.0040.0040.0040.0040.00504.356.278.4110.8413.6817.1321.7330.0040.0050.0050.0050.0050.00606.408.8011.4914.5718.2222.8130.0040.0050.0060.0060.0060.0060.00708.9211.8515.1519.0023.6630.0540.0050.0060.0070.0070.0070.0070.008011.9615.4819.5124.2830.4440.0050.0060.0070.0080.0080.0080.0080.009015.5819.7924.7030.8040.0050.0060.0070.0080.0090.0090.0090.0090.0010019.8824.9331.0640.0050.0060.0070.0080.0090.00100.00100.00100.00100.0011025.0131.2240.0050.0060.0070.0080.0090.00100.00110.00110.00110.00110.0012031.2740.0050.0060.0070.0080.0090.00100.00110.00120.00120.00120.00120.00图5降雨径流相关图（2）径流深预报误差评定。由给定的次洪Pa值，有效降雨PE值，次洪径流深R值，利用（1）中得到的降雨径流相关计算表进行线性插值或者查降雨径流相关图可得到次洪的模拟径流深，与实测值进行比较即可制定误差评定数据表（表五）。根据表格合格率及确定性系数，我们可以知道预报方案为甲等。0204060801001201400.0020.0040.0060.0080.00100.00120.00表五径流深误差评定表Pa(mm)有效降雨PE(mm)次洪径流深(mm)许可误差模拟R预报误差合格否40100.151.910.3850.001.9合格74.796.182.816.5680.82合格62.378479.450.33.3合格84.686.581.916.3881.10.8合格80.964.954.210.8455.81.6合格53.669.828.95.7833.44.5合格684015.63.1219.33.7不合格75.1100.88617.285.90.1合格75.610084.716.9485.60.9合格90101.898.419.68101.83.4合格61.771.144.48.8842.81.6合格73.7103.8851787.52.5合格74.9104.793.118.6289.63.5合格62.8100.976.315.2673.72.6合格86.6150145.520146.61.1合格56.977.547.19.4244.42.7合格8092.883.516.782.80.7合格6367.840.38.0640.80.5合格（续表五）79123.4113.720112.41.3合格62132.810420104.80.8合格76.863.3469.250.14.1合格48.294.154.910.9852.32.6合格55.7139.210620104.91.1合格65.5100.974.514.976.41.9合格75.298.284.316.8683.40.9合格合格率0.96确定性系数Dc0.99五、流域汇流方案的编制1、计算地面径流和地下径流用退水曲线中的洪水，分析地面径流和地下径流的分离点。从起涨点到拐点连直线，查取地面流量Qs。则应用式六计算实测地面径流为：21()*tttRssQsQgFRgsRsRss（式六）式中：Rs为由过程线上计算的径流深，Δt为时段长。Rss为实际的地面径流；Rgs为实际的地下径流。其中Qs的计算过程采用线性插值得到。地面和地下径流推求计算数据如下表表六所示：表六径流计算表Qs1.902.0030.30225.40704.701167.501191.501088.00878.50676.80498.00375.10269.00178.0093.8064.00Qg1.902.0010.1814.3218.4622.6026.7430.8835.0239.1643.3047.4451.5855.7259.8664.00Qs-Qg0.000.0020.12211.08686.241144.901164.761057.12843.48637.64454.70327.66217.42122.2833.940.00总径流深R(mm)58.00Rss(mm)51.77Rgs(mm)6.232、fc的试算。（1）采用试算法推求流域平均稳定下渗率，有时段降雨量及第五日Pa值，计算第一时段的蒸发量、时