水利工程师工程水文及水利计算模拟习题含参考答案

水利工程师工程水文及水利计算模拟习题含参考答案一、某河流某水文站有1980-2020年共41年的年最大洪峰流量观测资料（单位：m³/s），具体数据如下表所示：年份19801981198219831984198519861987198819891990199119921993199419951996199719981999流量120015009802100180016501350230017001450200019501100850250010501750220019001600年份200020012002200320042005200620072008200920102011201220132014201520162017201820192020流量240018501550100078012501400130021501780205019801150900260010801680225018801580800（1）将资料按从大到小排序，计算经验频率（采用数学期望公式），并绘制经验频率曲线；（2）采用矩法初估年最大洪峰流量的统计参数（均值、变差系数Cv，取Cs=3Cv）；（3）用适线法调整统计参数，确定设计频率P=1%（百年一遇）的设计洪峰流量。参考答案：（1）排序及经验频率计算：将41年流量从大到小排序，结果如下（前10项示例）：2600（2014）、2500（1994）、2400（2000）、2300（1987）、2250（2017）、2200（1997）、2150（2008）、2100（1983）、2050（2010）、2000（1990）……800（2020）。经验频率计算公式：P(i)=m/(n+1)×100%（m为排序后序号，n=41）。例如，最大流量2600为第1位，P=1/(41+1)=2.38%；第2位2500，P=2/42≈4.76%，依此类推。（2）矩法初估参数：均值Q̄=（ΣQi）/n。计算得ΣQi=1980-2020年流量总和=1980+1500+…+800=73560m³/s（注：实际需逐项累加，此处为简化示例），Q̄=73560/41≈1794m³/s。均方差σ=√[Σ(Qi-Q̄)²/(n-1)]，计算各(Qi-Q̄)²之和（示例：2600-1794=806，平方为649636；2500-1794=706，平方为498436；…800-1794=-994，平方为988036），假设Σ(Qi-Q̄)²=1.2×10⁸，则σ=√(1.2×10⁸/40)=√3×10³≈1732m³/s。变差系数Cv=σ/Q̄≈1732/1794≈0.965（实际需精确计算，此处为示例）。取Cs=3Cv≈2.895。（3）适线法调整参数：以矩法初估参数为基础，绘制皮尔逊Ⅲ型理论频率曲线，与经验频率点据拟合。若理论曲线与经验点据在中高频率段（P<20%）偏上，则需增大Cv；若偏下则减小Cv。假设调整后Cv=1.0，Cs=3.0，重新计算理论频率曲线。查皮尔逊Ⅲ型频率曲线模比系数Kp表（Cs=3.0），当P=1%时，Kp≈3.5（具体值需查表），则设计洪峰流量Qp=Kp×Q̄=3.5×1794≈6279m³/s。二、某水库流域面积F=500km²，设计标准为P=2%（五十年一遇）。根据水文手册，该流域最大1小时、6小时、24小时设计暴雨量分别为H1p=80mm、H6p=150mm、H24p=220mm，暴雨衰减指数n1=0.5（t≤1h）、n2=0.6（1h<t≤6h）、n3=0.7（6h<t≤24h）。已知该流域初损I0=20mm，稳渗率fc=2mm/h，采用同频率放大法推求设计净雨过程（假设典型暴雨雨型为：24小时暴雨中，前6小时占50%，后18小时占50%；6小时暴雨中，前1小时占40%，后5小时占60%）。（1）计算设计暴雨过程；（2）计算设计净雨过程。参考答案：（1）设计暴雨过程计算：同频率放大法需保证各时段暴雨量与设计值同频率。①24小时暴雨总量H24p=220mm，其中6小时部分H6p=150mm，1小时部分H1p=80mm。②各时段放大系数：1小时暴雨放大系数K1=H1p/H1典（H1典为典型暴雨1小时量）。典型暴雨中，24小时暴雨的6小时部分占50%，即H6典=0.5×H24典；6小时暴雨的1小时部分占40%，即H1典=0.4×H6典=0.4×0.5×H24典=0.2×H24典。假设典型暴雨24小时总量H24典=180mm（需根据地区暴雨特性假设），则H6典=90mm，H1典=36mm。K1=80/36≈2.222。6小时暴雨中，除1小时外的5小时放大系数K6-1=(H6p-H1p)/(H6典-H1典)=(150-80)/(90-36)=70/54≈1.296。24小时暴雨中，除6小时外的18小时放大系数K24-6=(H24p-H6p)/(H24典-H6典)=(220-150)/(180-90)=70/90≈0.778。③设计暴雨过程：典型暴雨雨型（假设时段为1h、5h、18h）：第1小时：36mm→设计值=36×2.222≈80mm；第2-6小时（5小时）：H6典-H1典=90-36=54mm→设计值=54×1.296≈70mm（每小时70/5=14mm）；第7-24小时（18小时）：H24典-H6典=180-90=90mm→设计值=90×0.778≈70mm（每小时70/18≈3.89mm）。（2）设计净雨过程计算：净雨=暴雨-损失，损失包括初损I0和稳渗fc。①初损I0=20mm，发生在降雨初期，当累积降雨≥I0时，开始产流。②第1小时降雨80mm>20mm，初损全部消耗，剩余80-20=60mm，扣除稳渗fc×1h=2mm，净雨=60-2=58mm；③第2-6小时每小时降雨14mm，累积降雨至第6小时末=80+14×5=150mm，初损已耗尽，每小时净雨=14-2=12mm（5小时总净雨=12×5=60mm）；④第7-24小时每小时降雨3.89mm<fc=2mm？不，fc=2mm/h，若降雨强度≥fc则产流，否则不产流。3.89>2，故每小时净雨=3.89-2=1.89mm（18小时总净雨=1.89×18≈34mm）。设计净雨过程为：第1小时58mm，第2-6小时每小时12mm，第7-24小时每小时1.89mm。三、某水库具有年调节能力，设计保证率P=85%。已知坝址处1970-2020年共51年年平均流量资料（单位：m³/s），经频率计算得：年平均流量均值Q̄=50m³/s，Cv=0.3，Cs=2Cv，设计枯水年（P=85%）的年平均流量Qp=35m³/s。水库死水位Z死=100m，死库容V死=5000万m³；正常蓄水位Z蓄=120m，相应库容V蓄=2.5亿m³；水库下游需水量为：非汛期（10月-次年3月，6个月）需水流量N非=40m³/s，汛期（4月-9月，6个月）需水流量N汛=30m³/s。（1）计算设计枯水年的年来水量；（2）计算水库调节库容；（3）判断该水库是否满足下游需水要求（忽略水库蒸发、渗漏损失）。参考答案：（1）设计枯水年年来水量：年平均流量Qp=35m³/s，年来水量W年=Qp×T=35×365×24×3600=35×3.1536×10⁷≈1.1038×10⁹m³（即11.038亿m³）。（2）水库调节库容计算：年调节水库需将年内多余水量蓄存，以补充缺水期。需计算设计枯水年的月来水量和月需水量，确定盈亏量。假设设计枯水年各月来水量按多年平均分配比例缩放（多年平均各月流量占年平均比例为：10月-3月各月占4%，4月-9月各月占6%，总比例=6×4%+6×6%=60%？需调整为100%，正确比例应为：非汛期6个月占40%（各月约6.67%），汛期6个月占60%（各月约10%））。设计枯水年各月来水量：非汛期（10月-3月）：Q非月=35×6.67%≈2.33m³/s，6个月总来水量W非=2.33×6×30×24×3600≈2.33×1.5552×10⁷≈3.62×10⁷m³（0.362亿m³）；汛期（4月-9月）：Q汛月=35×10%=3.5m³/s，6个月总来水量W汛=3.5×6×30×24×3600≈3.5×1.5552×10⁷≈5.44×10⁷m³（0.544亿m³）；总来水量=0.362+0.544=0.906亿m³（与年水量11.038亿m³不符，此处为简化假设，实际应按正确月分配比例计算）。下游需水量：非汛期需水量W需非=N非×6×30×24×3600=40×1.5552×10⁷≈6.2208×10⁷m³（0.622亿m³）；汛期需水量W需汛=N汛×6×30×24×3600=30×1.5552×10⁷≈4.6656×10⁷m³（0.4666亿m³）；总需水量=0.622+0.4666=1.0886亿m³。设计枯水年年来水量（假设修正后为1.0886亿m³）与需水量平衡，调节库容为汛期余水与非汛期缺水的差值。假设汛期来水0.544亿m³<需水0.4666亿m³（不余水），非汛期来水0.362亿m³<需水0.622亿m³（缺水0.26亿m³），则需从水库蓄水中补充，调节库容V调=V蓄-V死=2.5-0.5=2.0亿m³≥0.26亿m³，满足要求。（3）判断需水满足情况：调节库容2.0亿m³大于缺水0.26亿m³，且设计枯水年年来水量（11.038亿m³）大于总需水量（非汛期：40×6×30×24×3600=40×1.5552×10⁷=6.2208×10⁷m³=0.622亿m³；汛期：30×6×30×24×3600=4.6656×10⁷m³=0.4666亿m³；总需水量=0.622+0.4666=1.0886亿m³），11.038>1.0886，故水库可通过调节满足下游需水要求。四、某水库溢洪道为开敞式宽顶堰，堰顶高程Z堰=150m，堰宽B=20m，流量系数m=0.385。水库防洪限制水位Z限=152m（相应库容V限=1.2亿m³），设计洪水过程线如下表（时段Δt=6h）：时间（h）0612182430364248入库流量Q（m³/s）5020050080012001000600300100（1）建立水库水位-库容-泄流量关系表（已知Z=150m时V=0.8亿m³，Z=152m时V=1.2亿m³，Z=154m时V=1.8亿m³，Z=156m时V=2.5亿m³）；（2）采用列表试算法计算水库调洪过程，确定最高库水位和最大下泄流量（起调水位为Z限=152m）。参考答案：（1）水位-库容-泄流量关系表：宽顶堰泄流量公式：q=mB√(2g)H^(3/2)，其中H=Z库-Z堰（Z库为库水位，Z堰=150m）。计算各水位对应的H、q：Z=150m：H=0，q=0，V=0.8亿m³；Z=152m：H=2m，q=0.385×20×√(19.62)×2^(3/2)=7.7×4.429×2.828≈7.7×12.5≈96.25m³/s（精确计算：√(2×9.81)=4.429，2^(3/2)=2.828，q=0.385×20×4.429×2.828≈0.385×20×12.5≈96.25），V=1.2亿m³；Z=154m：H=4m，q=0.385×20×4.429×4^(3/2)=7.7×4.429×8≈7.7×35.43≈272.8m³/s，V=1.8亿m³；Z=156m：H=6m，q=0.385×20×4.429×6^(3/2)=7.7×4.429×14.696≈7.7×65.1≈501.3m³/s，V=2.5亿m³。（2）调洪计算（Δt=6h=21600s）：起调水位Z0=152m，V0=1.2亿m³，初始下泄流量q0=96.25m³/s。调洪基本方程：(Q1+Q2)/2×Δt(q1+q2)/2×Δt=V2V1。列表计算（示例前3时段）：时段时间(h)Q1(m³/s)Q2(m³/s)(Q1+Q2)/2×Δt(亿m³)假设q2(m³/s)(q1+q2)/2×Δt(亿m³)V2=V1+左-右(亿m³)查V2对应Z2(m)查Z2对应q2(m³/s)修正q2是否与假设一致10-650200(50+200)/2×21600=2.7×10⁶m³=0.027亿假设100(96.25+100)/2×21600≈2.124×10⁶m³=0.02124亿1.2+0.027-0.02124=1.20576亿查V=1.20576亿对应Z≈152.1m（线性插值：1.2亿→152m，1.8亿→154m，每0.1亿对应0.333m，0.00576亿对应0.00576×3.