河海大学水文分析与计算课程设计报告

水文分析计算课程设计报告书学院：水文水资源专业：水文与水资源工程学号：姓名：指导老师：梁忠民、李国芳2015年06月12日南京目录3、9278889021、设计任务推求江西良田站设计洪水过程线,本次要求做P校,即推求%(t)。2、流域基本概况良田是赣江的支流站。良田站以上控制的流域面积仅为，属于小流域，如右图所示。年降水均值在1500～1600mm之内，变差系数Cv为，即该地区降雨充沛，年际变化小，地处湿润地区。暴雨集中。暴雨多为气旋雨、台风雨，季节为3～8月，暴雨历时为2～3日。3、资料情况及计算方案拟定资料情况设计站（良田）流量资料缺乏，邻近站雨量资料相对充分，具体如表3-1：表3-1良田站及邻近地区的实测暴雨系列、历时洪水、特大暴雨资料（设计站(良田)流量资料缺乏，邻近站雨量资料相对充分。）方案拟定本次课设采用间接法推求设计洪水，即是由推求的设计暴雨，经过产汇流计算得到设计洪水。示意图如下：4、设计暴雨XP(t)的计算设计暴雨Xp(t)计算为推求该区域设计面降雨量，选取吉安、桑庄、寨头与峡江四站降雨检验该区降雨是否选同一总体。选择四站1981～2013年数据（74年出现极值暴雨，不参加检验），对各站数据取自然对数，对转换后数据进行均值与方差检验，各站转换后系列的均值及方差见表4-1。表4-1吉安、桑庄、寨头与峡江站最大一日降雨资料均值、方差情况站名实测暴雨流量系列特大暴雨、历史洪水良田75～78（4年）Q=216m3/s，N=80（转化成X1日，移置峡江站）峡江53～80（28年）吉安36～80（45年）桑庄57～80（24年）X1日寨头57～80（24年）沙港特大暴雨X1日（移置到寨头站）1）均值检验选取均值差异最大的吉安站（X1）和峡江站（X2）两站进行检验。假设H:X1=X2构造统计变量：=取α=，查得|tα/2|=|t|，接受假设H，即可认为吉安、桑庄、寨头与峡江站均值相等。2）方差检验选取方差差异最大的桑庄站（S1）和寨头站（S2）两站进行检验。假设H:S1=S2构造统计变量：=查得Fα/2(55,55)=,F1-α/2(55,55)=,所以接受域为（,）,则可以认为统计量F满足接受域,即认为四个站的降雨量数据满足方差相等。综上所述，可认为区域降雨资料来自同一总体，可以进行综合。（1）特大值处理峡江站特大暴雨由良田站历史洪水转换而来，则良田站（峡江站）的X1日=。公式如下（其中，Q为地表净峰流量（m3/s），m为汇流参数，取，F为流域面积（km2），L为出口断面沿主河道至分水岭的最长距离（km），J为沿L的坡面和河道平均比降，tc为净雨历时（h），为汇流历时（h），R为地表径流深（mm），n为暴雨参数，取，为稳渗率（mm/h），取，用良田站计算）：沙港站特大暴雨取重现期N=150年，放置寨头站进行频率计算。桑庄站最大一日暴雨取重现期N=150年。（2）排位分析及频率计算先对吉安、峡江、桑庄、寨头四个站进行频率计算，要考虑各站可能存在的特大暴雨系列值。各站的频率计算见表4-2表4-2各站频率计算序号吉安峡江桑庄寨头雨量值频率雨量值频率雨量值频率雨量值频率mm％mm％mm％mm％1396231611/21/2221211()/()()/(2)XYmntmsnsmn统计量456781289101391112129131201211412511515114161718191141122021222324992526271051032894293099313233349535363786383986784041424344457846（3）四站适线结果图吉安站：峡江站：桑庄站：寨头站：（4）四站均值计算结果及其适线图图4-25个雨量站年最大雨量频率曲线（5）推求X1日，P=%由上表查处X（1日，P=%）=，则X（24h，P=%）=*X（1日，P=%）=。4．2计算各种历时同频率雨量Xt,P4748797049756950767451525354555657585970606162696364656667636869707172737475567677由暴雨公式推算t=3h，6h，9h，12h，15h，18h，21h，24h的设计暴雨值，结果如表4-3。表4-3各短历时设计暴雨4．3选典型放大推求XP(t)4．4产汇流计算在设计暴雨中，由于稀遇频率的设计暴雨量很大，损失相对较小，因此，一般采用简化模型。即前段降雨尽量满足土壤蓄水量，即初损，而后假定稳渗率，算定地面径流深R上和地下径流深R下，再列表求出，如与假定的相符，则假定的即为所求值。径流分割点绘良田站、毛背站、以及的流量和雨量过程。利用平割法计算Rt，利用斜线分割法进行水源划分。先寻找洪水过程的直接径流终止点B，然后用斜线连接起涨点A与终止点B，将实测流量过程线分为两部分，斜线AB上部分为直接径流RS，下部分为地下径流RG。本次课设先采用梯形面积法求得RG，再用Rt-RG求得RS。计算初损，其中,x为总降雨量，R为总径流，R=Rt。率定μ1、假定μ，计算R下和R上。时段取△t=1h。若时段的降雨量累积和小于初损量，则全部雨量补充初损值。当累积降雨量大于初损值时，开始产生径流。当X（t）≤μ时，全部产生底下径流，则R下（t）=μ，R上（t）=0；当X（t）μ时，R下（t）=μ，R上（t）=X（t）-μ。2、判断μ值的正确性。若RS=ΣR上，RG=ΣR下，则μ值即为良田站所确定的μ值；否则重新假定μ，转1。4．4．1径流划分及稳渗μ值率定（1）良田站76年6月17日表4-4良田站1976年6月17日次洪μ值率定（2）毛背站75年5月13日表4-5毛背站1975年5月13日次洪μ值率定（3）毛背站76年7月9日表4-6毛背站1976年7月9日次洪μ值率定（4）毛背站77年6月26日表4-7毛背站1977年6月26日次洪μ值率定（5）各场次洪水的径流分割结果及率定所得μ值表4-8各场次洪水的径流分割结果及率定所得μ值4．4．2地表汇流地面汇流的计算方法有经验公式法（如单位线、经验公式等）和推理公式法（如等流时线法、水科院推理公式法、推理△过程线法、汇流系数法等）。本次采用八省一院公式。（1）（2）其中，Q为地表洪峰流量（m3/s），m为汇流参数，取，F为流域面积（km2），L为出口断面沿主河道至分水岭的最长距离（km），J为沿L的坡面和河道平均比降，tc为净雨历时（h），为汇流历时（h），R为地表径流深（mm），n为暴雨参数，取，为稳渗率（mm/h），取。对于良田流域，，因此采用式（2）。（1）m初值的确定表4-9m初值确定故四站综合,=+++/4=（2）m值的检验（以毛背站为例）计算步骤（1）根据表18中计算的初值m、八省一院公式，对该次的降雨过程的每个推求，（2）对每个Qi，假定过程线为三角形，底宽为Qi出现在处。（3）将各时段的三角形过程进行叠加，与实测洪水（扣除地下径流）对比。如相差太大重新假定m，重新计算。表4-10m=时的地表流量Qs及对应底宽T表4-11m=时的地表流量演算毛背站1976年7月9日地表径流过程线（m=）表4-12m=时的地表流量Qs及对应底宽表4-13m=时的地表流量演算毛背站1976年7月9日地表径流过程线（m=）可以看出m=时实测与计算径流量线拟合的更好。（3）m值综合（以毛背站为例）点绘各次洪水的Q/F～m图，取上端趋于稳定的m值，为设计暴雨之m值。表4-14各场次洪水Qs/F～m关系统计由各次洪水的Q/F～m图，取上端趋于稳定的m值（为），设计暴雨之m值（实际运用时，选用通过检验后的四站m值进行综合，为）。4．5由设计暴雨XP(t)推求QP(t)4．5．1产流计算由可知，mm/h，再按前面过程分水源，（I=0）。计算结果见表20。表4-15产流计算及分水源4．5．2地面汇流表4-16Qs、Ti计算表4-17地表汇流计算4．5．3地下汇流计算Ts（Δt=1h）为37，则Tg（Δt=1h）为74，Rg（mm）为，Qgm(m3/s)为。具体计算表格如表4-17。表4-17地下汇流计算表格4．5．4设计洪水过程线计算结果见表4-18。表4-18设计洪水过程线5、心得体会一周的水文分析计算课设很快结束了。由于刚结束水利计算枯燥的课程设计，本以为水文分析的工作量应该会有所减少，但事实证明其过程并不如预想的美好。由于种种原因，我第一天都没有进行课设，所以直到第二天开始入手后，始终有一种追赶者的感觉。可我并没有敷衍过程，每一步都是经过了自己的认真计算。在具体计算分析过程中，我遇到了以下几点问题或收获：（1)利用斜线法进行径流分割时，更多地是靠自我感觉，而并没有定量地算出应选择某一段，所以每个人会存在选取误差；（2）对m值检验时，由实测值和计算值的图形拟合程度来判断是否正确。但拟合程度也是凭个人主观意识，并没有具体的误差要求。（估计是老师很善良，减少了我们的工作量）；（3）在产汇流计算时，涉及到了大量的插值计算工作。之前我是自己在每一行或列都插入公式，后来运用了trend函数，不过在行列树比较多的情况下也十分繁琐。但如果能运用VBA程序，就能有效解决这个问题。这也是我如后需要努力的地方。1960195919581957195619551954注:括号()中数据为缺失数据,本次采用相应的多年平均值代替。