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

.专业资料水文分析计算课程设计报告书学院：水文水资源专业：水文与水资源工程学号：姓名：指导老师：梁忠民、国芳2015年06月12日南京目录HYPERLINK1、设计任务 1HYPERLINK2、流域概况 13、HYPERLINK资料情况及计算方案拟定 1HYPERLINK4、计算步骤及主要成果 2HYPERLINK4．1设计暴雨Xp(t)计算 2HYPERLINK4．1．1区域降雨资料检验 2HYPERLINK4．1．2频率分析与设计雨量计算 3HYPERLINK4．2计算各种历时同频率雨量Xt,P 9HYPERLINK4．3选典型放大推求XP(t) 9HYPERLINK4．4产汇流计算 9HYPERLINK4．4．1径流划分及稳渗μ值率定 12HYPERLINK4．4．2地表汇流 17HYPERLINK4．5由设计暴雨XP(t)推求QP(t) 18HYPERLINK4．5．1产流计算 18HYPERLINK4．5．2地面汇流 18HYPERLINK4．5．3地下汇流计算 19HYPERLINK4．5．4设计洪水过程线 20HYPERLINK5、心得体会 22专业资料1、设计任务推求良田站设计洪水过程线,本次要求做P校,即推求Q0.01%(t)。2、流域基本概况良田是赣江的支流站。良田站以上控制的流域面积仅为44.5km2，属于小流域，如右图所示。年降水均值在1500～1600mm之，变差系数Cv为0.2，即该地区降雨充沛，年际变化小，地处湿润地区。暴雨集中。暴雨多为气旋雨、台风雨，季节为3～8月，暴雨历时为2～3日。3、资料情况及计算方案拟定3.1资料情况设计站（良田）流量资料缺乏，邻近站雨量资料相对充分，具体如表3-1：表3-1良田站及邻近地区的实测暴雨系列、历时洪水、特大暴雨资料站名实测暴雨流量系列特大暴雨、历史洪水良田75～78（4年）Q=216m3/s，N=80（转化成X1日，移置峡江站）峡江53～80（28年）36～80（45年）桑庄57～80（24年）X1日=416mm，N=100～150（74.8.11）寨头57～80（24年）沙港特大暴雨X1日=396mm，N=100～150（69.6.30）（移置到寨头站）（设计站(良田)流量资料缺乏，邻近站雨量资料相对充分。）3.2方案拟定本次课设采用间接法推求设计洪水，即是由推求的设计暴雨，经过产汇流计算得到设计洪水。示意图如下：4、设计暴雨XP(t)的计算4.1设计暴雨Xp(t)计算4.1.1区域降雨资料检验为推求该区域设计面降雨量，选取、桑庄、寨头与峡江四站降雨检验该区降雨是否选同一总体。选择四站1981～2013年数据（74年出现极值暴雨，不参加检验），对各站数据取自然对数，对转换后数据进行均值与方差检验，各站转换后系列的均值及方差见表4-1。表4-1 、桑庄、寨头与峡江站最大一日降雨资料均值、方差情况均值检验选取均值差异最大的站（X1）和峡江站（X2）两站进行检验。假设H:X1=X2构造统计变量：=0.78取α=0.10，查得|tα/2|=1.68>|t|，接受假设H，即可认为、桑庄、寨头与峡江站均值相等。方差检验选取方差差异最大的桑庄站（S1）和寨头站（S2）两站进行检验。假设H:S1=S2构造统计变量：=1.19查得Fα/2(55,55)=1.567,QUOTEF1-α/2(55,55)=0.638,所以接受域为（0.638,1.567）,则可以认为统计量F满足接受域,即认为四个站的降雨量数据满足方差相等。综上所述，可认为区域降雨资料来自同一总体，可以进行综合。4.1.2频率分析与设计雨量计算（1）特大值处理 峡江站特大暴雨由良田站历史洪水转换而来，则良田站（峡江站）的X1日=293.9mm。公式如下（其中，Q为地表净峰流量（m3/s），m为汇流参数，取0.7，F为流域面积（km2），L为出口断面沿主河道至分水岭的最长距离（km），J为沿L的坡面和河道平均比降，tc为净雨历时（h），QUOTE为汇流历时（h），R为地表径流深（mm），n为暴雨参数，取0.6，QUOTE为稳渗率（mm/h），取4.5，QUOTE用良田站计算）：沙港站特大暴雨取重现期N=150年，放置寨头站进行频率计算。桑庄站最大一日暴雨取重现期N=150年。（2）排位分析及频率计算先对、峡江、桑庄、寨头四个站进行频率计算，要考虑各站可能存在的特大暴雨系列值。各站的频率计算见表4-2表4-2各站频率计算序号峡江桑庄寨头雨量值频率雨量值频率雨量值频率雨量值频率mm％mm％mm％mm％1177.51.3293.91.2415.51.03961.02169.82.6272.92.8227.24.4173.74.431613.8192.94.4168.26.1170.86.14159.55.1185.76.0158.27.8163.07.85154.66.4183.37.6150.29.5155.19.56145.27.7178.99.2142.711.2148.611.27143.59.0169.210.8142.612.9146.212.9812810.3156.212.4133.214.6145.114.69126.711.5152.714.0130.616.4140.716.410125.612.815.6129.518.1135.418.111123.214.1135.917.2124.819.8133.319.812120.315.412918.8123.821.5126.621.51312016.7127.820.412123.2122.223.214115.817.912521.911524.9116.924.915115.619.2121.823.5112.726.611426.616115.320.5115.125.1111.428.3113.328.317114.521.8113.826.7111.230.0110.230.018114.323.1112.928.3109.831.7110.031.71911424.411229.9109.433.4107.133.420113.225.6111.631.5108.735.1106.535.121112.926.9110.633.1103.836.8106.036.822109.828.2110.434.7101.238.5104.738.523109.529.5106.236.399.940.3104.540.324107.530.8105.537.997.242.09942.025107.432.1105.239.597.143.797.643.726106.633.3103.541.195.445.497.545.42710534.610342.794.347.197.247.128102.335.9101.544.29448.896.748.829101.537.299.345.892.550.595.450.530101.438.59947.490.952.295.452.23199.539.795.349.090.753.990.253.93299.141.090.950.688.355.689.555.63397.642.390.652.285.157.388.157.3349543.687.553.885.159.084.859.03588.544.987.555.483.260.784.760.73688.446.286.857.082.562.483.762.43787.747.48658.682.364.283.464.23887.548.785.860.279.565.982.865.9398650.085.661.87867.682.267.64084.451.384.563.477.669.379.569.34184.352.683.765.076.771.079.171.04282.253.882.166.574.572.778.572.74381.955.180.968.172.474.476.574.44481.356.480.469.772.476.175.676.14580.857.77871.370.877.874.177.84680.359.077.472.970.779.570.679.54779.660.377.274.570.281.270.581.2487961.576.976.17082.969.782.94976.662.87577.76984.669.584.6507664.17479.364.486.368.286.35174.265.471.880.963.788.163.588.15273.966.770.982.560.889.862.089.85373.167.970.484.160.691.559.791.55472.669.269.985.759.493.259.093.25571.770.569.487.357.894.958.494.95671.571.861.588.854.196.657.796.65770.973.160.690.441.398.350.498.35870.574.460.392.0597075.658.893.66069.976.950.795.26169.878.245.296.8626979.544.598.46368.780.86467.982.16564.383.36663.684.6676385.96862.287.26960.988.57059.589.77158.491.07257.992.37357.293.67456.594.9755696.27655.697.47744.398.7（3）四站适线结果图站：峡江站：桑庄站：寨头站：四站均值计算结果及其适线图图4-25个雨量站年最大雨量频率曲线（5）推求X1日，P=0.01%由上表查处X（1日，P=0.01%）=461.9mm，则X（24h，P=0.01%）=1.1*X（1日，P=0.01%）=508.09mm。4．2计算各种历时同频率雨量Xt,P由暴雨公式QUOTE推算t=3h，6h，9h，12h，15h，18h，21h，24h的设计暴雨值，结果如表4-3。表4-3各短历时设计暴雨4．3选典型放大推求XP(t)4．4产汇流计算在设计暴雨中，由于稀遇频率的设计暴雨量很大，损失相对较小，因此，一般采用简化模型。即前段降雨尽量满足土壤蓄水量，即初损，而后假定稳渗率QUOTE，算定地面径流深R上和地下径流深R下，再列表求出QUOTE，如与假定的QUOTE相符，则假定的QUOTE即为所求值。径流分割点绘良田站76.6、毛背站75.5、76.7以及77.6的流量和雨量过程。利用平割法计算Rt，利用斜线分割法进行水源划分。先寻找洪水过程的直接径流终止点B，然后用斜线连接起涨点A与终止点B，将实测流量过程线分为两部分，斜线AB上部分为直接径流RS，下部分为地下径流RG。本次课设先采用梯形面积法求得RG，再用Rt-RG求得RS。计算初损QUOTE，其中,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下，则μ值即为良田站76.6所确定的μ值；否则重新假定μ，转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地表汇流地面汇流的计算方法有经验公式法（如单位线、经验公式等）和推理公式法（如等流时线法、水科院推理公式法、推理△过程线法、汇流系数法等）。本次采用八省一院公式。QUOTE（1）QUOTE（2）其中，Q为地表洪峰流量（m3/s），m为汇流参数，取0.7，F为流域面积（km2），L为出口断面沿主河道至分水岭的最长距离（km），J为沿L的坡面和河道平均比降，tc为净雨历时（h），QUOTE为汇流历时（h），R为地表径流深（mm），n为暴雨参数，取0.6，QUOTE为稳渗率（mm/h），取4.5。对于良田流域，QUOTE，因此采用式（2）。（1）m初值的确定表4-9m初值确定故四站综合,=(0.31+0.6+0.8+0.75)/4=0.62（2）m值的检验（以毛背站76.7为例）计算步骤（1）根据表18中计算的初值m、八省一院公式，对该次的降雨过程的每个QUOTE推求QUOTE，QUOTE（2）对每个Qi，假定过程线为三角形，底宽为QUOTEQi出现在QUOTE处。（3）将各时段的三角形过程进行叠加，与实测洪水（扣除地下径流）对比。如相差太大重新假定m，重新计算。表4-10m=0.62时的地表流量Qs及对应底宽T表4-11m=0.62时的地表流量演算毛背站1976年7月9日地表径流过程线（m=0.62）表4-12m=0.56时的地表流量Qs及对应底宽表4-13m=0.56时的地表流量演算毛背站1976年7月9日地表径流过程线（m=0.56）可以看出m=0.56时实测与计算径流量线拟合的更好。（3）m值综合（以毛背站76.7为例）点绘各次洪水的Q/F～m图，取上端趋于稳定的m值，为设计暴雨之m值。表4-14各场次洪水Qs/F～m关系统计由各次洪水的Q/F～m图，取上端趋于稳定的m值（为0.8），设计暴雨之m值（实际运用时，选用通过检验后的四站m值进行