总岗山水库病害整治工程计算书.doc

四川大学本科毕业设计 总岗山水库枢纽病害整治工程设计目录1水库洪水复核11.1流域概况11.1.1流域自然地理11.1.2水文气象11.2设计洪水21.2.1设计洪水标准复核21.2.2基本资料的收集和整理21.3水库调洪演算41.3.1推求泄流曲线41.3.2调洪演算51.4坝顶高程复核91.4.1计算主坝正常运用情况下的坝顶超高：101.4.2计算非常运用情况下的坝顶超高111.4.3主坝坝顶高程复核：121.4.4计算副坝正常运用情况下的坝顶超高121.4.5 计算副坝非常运用情况下的坝顶超高131.4.6副坝坝顶高程复核142渗流安全评价152.1主坝渗流计算分析152.1.1目的152.2主坝渗流运行分析152.2.1渗流分析的方法152.2.2计算过程172.3汉王副坝渗流计算分析262.3.1目的262.3.2渗透系数262.3.3副坝渗流计算分析273结构安全评价363.1主坝坝坡稳定计算分析363.1.1计算分析资料363.1.2安全系数363.1.3土体的抗剪强度指标363.1.4计算工况363.2主坝坝坡稳定运行分析373.2.1 运用stab程序计算坝坡稳定安全系数373.3汉王副坝坝坡稳定计算分析413.3.1计算分析资料413.3.2安全系数413.3.3土体的抗剪强度指标423.4副坝坝坡稳定运行分析423.5溢洪道结构安全分析463.5.1溢洪道结构安全分析463.5.2溢洪道水力计算474.大坝整治方案计算504.1综合评价504.2整治的必要性504.3整治措施及方案对比514.3.1主坝整治方案一：514.3.2主坝整治方案二604.3.3 坝坡整治方案比选69总岗山水库枢纽病害整治工程设计计算书701水库洪水复核1.1流域概况1.1.1流域自然地理总岗山水库位于洪雅县境西部，是洪雅县内的一座中型水库。洪雅县地处四川盆地西南边缘，地理坐标为东经1024910332，北纬29243000全县最高海拔3090m，最低海拔417.5m。县内气候温和湿润，年降雨量1435.5mm，(1971-2000)年日照1006.1小时，年无霜期307天，年平均气温16.7，属中亚热带湿润气候。全县幅员面积1896.49km2，地形由西南向东北高低梯次变化形成高山、中山、深丘、浅丘、台地、河谷、平坝， 地貌以山地丘陵为主，河谷平坝分布在青衣江、花溪河两岸，素有“七山二水一分田”之称。总岗山水库位于洪雅县西部总岗山南麓汉王乡境内，青衣江左岸一级支流香樟河中游，坝址地理坐标为东经10311，北纬2955。西北侧向东南，地形为总岗山主动脉相夹而构成的带状小谷地区，流域内多山地，少平川，山地高者海拔在1000m以上，低者在400m至700m之间。流域内冲沟相对发育，受岩性、地质构造及岩层裂缝的控制，大小支沟的分布在平面上略呈网络状曲折相接。1.1.2水文气象香樟河流域气候类型属于亚热带湿润气候,在气候区划上属于四川盆地亚热带湿润气候区的盆地西部边缘区和盆地西部区。根据洪雅县气象站2007年气象资料统计：多年平均气温16.7，年极端最高气温36.2，年极端最低气温3.3,多年平均相对湿度85；多年平均降水量1435.5mm,多年平均水面蒸发量970.9mm（20cm蒸发皿）；多年平均无霜期302天；多年平均风速1.0m/s。设计区域径流主要由降水补给。多年平均年径流深1123mm左右，变差系数0.25。径流年内分配不均匀，洪枯水变幅大。洪水由暴雨形成。暴雨集中在610月，相应的在这时期形成洪水。实测多年年最大24小时暴雨量约410.8mm，多年平均年24小时平均暴雨量为145mm，变差系数0.53；实测多年年最大6小时暴雨量为313.8mm，占到最大24小时暴雨量的76左右。暴雨集中、强度大，一般暴雨过程12天，因此，洪水过程陡涨陡落，多呈单峰过程，持续时间不长。1.2设计洪水1.2.1设计洪水标准复核总岗山水库是一座以灌溉为主，兼有防洪、养殖和生活供水等综合利用的中型水库。水库主坝为均质土坝，最大坝高34.2。水库正常水位689.70，死水位666.70；有效库容2554万。根据水利部sl252-2000水利水电工程等级划分及洪水标准的规定，水库属等枢纽工程，主要建筑物为3级，采用100年一遇洪水设计、2000年一遇洪水校核。1.2.2基本资料的收集和整理（1）流域特征参数根据110000航测图量算出设计流域面积（f）、河长（l）和沿河长的平均坡度（j），如表1.2.2-1。表1.2.2-1 总岗山水库流域特征值表f（km2）l（km）j（）28205.3（2）库容曲线本次根据洪雅县水利局提供的库区1/2000的实测地形图，对总岗山水库的水位库容曲线进行了复核，总岗山水库的水位库容关系见表1.2.2-2和图1.2.2-1。表1.2.2-2 总岗山水库水位-库容关系表水位（m）680681682683684685686687688689690691692693694库容（万m3）95610861225136615201686186420502240244026602900314033803660图1.2.2-1 总岗山水库水位-库容曲线(3)设计洪水过程线根据总岗山水库大坝安全评价报告，取推理公式法（手册查参）计算成果： 表1.2.2-3 推理公式法推求的设计洪水过程线计算成果（手册查参数）序列p=0.05%p=0.1%p=1%p=2%历时（h）流量历时（h）流量历时（h）流量历时（h）流量103.4603.4603.4603.4620.85632.20.86529.60.916210.93718.231.11360.81.12555.71.1938.41.2183341.5411181.5581081.64873.41.68762.651.9692331.9912122.1061432.15612262.3973482.4233172.5642132.62418172.8254622.8564213.0222833.09324083.0825493.1165003.2973363.37428493.4245773.4625263.6633533.749299104.1955494.2415004.4873364.593284115.1374625.1934215.4942835.624240126.6783486.753177.1422137.311181139.162339.262129.79814310.0291221415.4111815.57810816.48373.416.87162.61522.25860.822.50255.723.80838.424.369331627.39532.227.69429.629.3022129.99318.21733.3883.4633.7523.4635.7123.4636.5543.46洪峰流量（）574522350296洪量（万）1825.81683.6712081051.68 图1.1.2-2推理公式法推求的设计洪水过程线（手册查参数）1.3水库调洪演算1.3.1推求泄流曲线总岗山水库溢洪道为正槽式宽顶堰控制放水，堰顶高程689.7m，堰顶现状净宽b=35m，参照水力计算手册（第二版，李炜主编），其流量系数m取为0.385，侧收缩系数取0.88。 q溢=2gmbh32泄洪建筑物为无闸门表面溢洪道，求水库的泄流曲线以及蓄泄方程曲线： 其中： q溢 溢洪道的下泄流量； b 溢洪道堰顶净宽； h 溢洪道堰上水头； m 流量系数。 q溢=2gmbh32以690m高程为例计算： = 29.80.38535690-689.732=8.627m3s 则各高程的计算结果整理如下：表1.3.1 蓄泄方程计算成果表水位(m)堰上水头(m)泄流量 q库容（万m3）689.70025886900.38.62726476911.377.81729126922.3183.12631396933.3314.72433826944.3468.1253662由上表绘制q泄v曲线如下图 图1.3.1 泄流曲线1.3.2调洪演算拟采用列表试算法，逐时段求解得出水库的下泄流量和蓄水量。其主要步骤如下：1 引用水库的设计洪水过程线。2 根据水库容积曲线v=fz和泄洪建筑物的泄流能力，应用公式求出下泻流量与库容的关系曲线q=f(v)。 m3/s3 根据水库汛期的控制运用方式，确定调洪计算的起始条件，即起调水位和相应的库容、下泄流量。本次计算中，采用溢洪道堰顶高程314.53为起调水位，相应的库容为16.56万m3、下泄流量为04 选取合适的计算时段t(以s为计算单位)，由设计洪水过程线摘录q1，q2。5 确定起始计算时刻的q1，v2值，然后列表计算，计算过程中，对每一计算时段的q1，v2 值都要进行试算。具体试算方法见详细的计算过程。6 将入库洪水 qt 过程和计算所得的泄流 qt 过程绘制在同一张图上，若计算所得的最大下泄流量 qm 正好是两线的交点，说明计算正确。否则，应缩小 qm附近的计算时段t，重新进行计算，直到计算的qm正好是两线的交点为止。7 由qm查q=f(v)关系图，可得最高洪水位的库容由减去起调水位相应库容，即得水库为调节该次入库洪水所需的调洪库容 vm 。再由查水库库容曲线，就可以得到最高洪水位 zm 。具体算法如下：1 列出水库水位库容关系：见表1.2.2-12 计算并绘制水库的 q=fv 关系曲线。（图1.2.2-1）3 确定调洪的起始条件。该水库溢洪道无闸门控制，因此起调水位就是堰顶高程，即689.7m。相应的库容为2588万m3 ，初始下泻流量为0。4 计算平均时段入库流量和时段入库流量。首先计算p=5% 洪水过程线，初选计算时段t=1.125h 填入表格第（1）列。 表中：第（1）列为初选计算时段t。 第（2）列为按计算时段摘录的入库洪水流量q。q=q1+q22=3.46+55.72=29.58(m3s) 第（3）列为计算所得的时段平均入库流量 q“=qt=29.5836001.125=11.98（万m3 第（4）列为时段入库水量5 逐时段试算求泄流过程 qt 。 以第一时段为例，假设q2 =1.53 m3s，q1+q22=0+1.532=0.765m3s则 ， 第（6）列为平均下泻流量 tq1+q22=36001.1250.765=0.31 万m3 第（7）列为时段下泻水量v=tq1+q22-tq1+q2 第（8）列为时段内水库存水量变化 v=11.98-0.312=12.60万m3 第（9）列为水库存水量v。初始存水量为2588万m3，则第一个时段间隔之后的存水量 v2=2588+v=2588+12.60=2600.6万m3 第（10）列为水库水位。已知v1=2588万m3, q1=0, q1=0, q2=55.7m3/s(查设计洪水过程线)假设q2=10m3/s,算得经过第一个时间间隔后水库的存水量为2598万m3，查蓄泄方程关系表，可以得出对应的泄流量q2=1.53m3/s，所以与原假设不符，需要重新假设。 假设q2 =1.53 m3s，按上表算得v2=2599万m3, 查蓄泄方程关系表，可以得出当泄流量q=1.69时，误差0.16小于允许误差，由此可以确定该假设成立。 重复类似的试算过程，可依次得到各个时间间隔的泄流过程qt。以表中第（1）、（2）、（5）列相应的数值，绘制qt与qt，但是，两线的交点并不是正好在qt上，所以需要再在23小时内插入时间间隔，重新试算。表1.3.4-2 p=0.1%调洪计算成果表t入库流量q入库流量m3/s时段入库流量(万m3)下泄流量m3/s平均下泄流量m3/s时段下泄水量(万m3)时段内水库存水量(万m3)水位（m）1.1253.4629.5811.9800.7650.312588689.71.29855.71.53689.90.6933.7186.3587.14020.7659.72600690.20.346500408.5101.9504511.22677690.20.77952651363.97562.57.82767690.6500513143.910087.524.52823691.82942691.3图1.3.4-2 p=0.1%调洪计算成果同理，p=1%的调洪演算过程与上述过程类似，其计算结果如下：图1.3.4-3 p=1%调洪计算成果调洪计算结果如表1.3.4-3表1.3.4-3 总岗山水库调洪计算成果表洪水频率 %起调 水位 (m)溢流 堰宽 (m)最高水位 (m)相应库容 (万m3)滞洪库容(万m3)最大下泄流量（m3/s)0.1689.735692.543267.28675.75243.661689.735691.663058.80467.27139.16由表知，校核洪水位（p=0.1%）为692.54，洪水最大下泄流量为243.66m3/s；设计洪水位（p=1%）为691.66，最大下泄流量139.16m3/s。 1.4坝顶高程复核为了保证库水位不超过或溅过坝顶，坝顶高程在水库正常运用和非常运用的静水位以上应有足够的超高。参考碾压式土石坝设计规范（sl2742001），按下式计算坝顶超高： 式中： y为坝顶超高；r为波浪在坝顶上的设计爬高；为安全超高，按规定采用；e为最大风壅水面高度。 波浪要素涉及以下几个因素：风速，吹程，坝前水深，风向与坝轴线法线方向的夹角，波高，上游面坡角，坝坡护面粗糙系数。坝顶高程要分别计算正常情况和非常情况，得出两个坝顶高程计算值后取大值。设计洪水位加正常运用情况的坝顶超高；校核洪水位加非常运用情况的坝顶超高。1.4.1计算主坝正常运用情况下的坝顶超高：1) 计算波浪爬高r：波浪爬高按蒲田试验站公式计算。先计算平均爬高rm,再计算设计爬高r。平均爬高按下式计算： 式中：k斜坡糙率渗透性系数，其中砌石取.0750.8，这里取0.75kv经验系数，与v/gh的无量纲大小有关，v/gh=10.5/9.824.18=0.682,查表得kv=1m边坡系数 取2.5hm平均波高 =0.104mlm平均波长 lm=2:5hm=250.104=2.6m因此平均波浪爬高为： 由于hm=0.104m，h=34.2m，hm/h=0.00304，三级大坝采用累积频率为1%的爬高值，查表得r/rm=2.23，所以得r=0.323m2) 计算风壅水面高度e： 式中： k为综合摩阻系数,取k=3.610-6；为计算风速，该坝为3级坝，故计算风速取多年平均风速的1.5倍，多年平均风速为1.0m/s，故v0=1.01.5=1.5m/sd为风区长度，取300m；为坝前水域的平均水深，hm=691.66-660=31.66；为计算风向与坝轴线的法线间的夹角, 由于风的方向垂直坝面取0则本设计的最大风壅高度为：3) 安全超高hc：按iii级坝和正常运用情况查表得hc=0.7m。综上所述，正常运用情况下的坝顶超高为 h=0.323+0.7+0.0000039=1.023m1.4.2计算非常运用情况下的坝顶超高非常运用情况下的坝顶超高计算公式与正常情况下一样，所不同的是风速采用多年平均最大风速v=1m/s，坝前水深h=692.54-660=32.54m 1) 计算波浪爬高r：k=0.750.80，采用0.75；取g=9.8m/s2，求得无维量，查表得经验系数kw=1；取风向与坝轴垂线的夹角为0，折减系数kb=1；初拟坝坡m=2.5；吹程d=300m，则平均波高为：平均波长为 将以上各值代入rm式得根据hm=0.0668m、h=32.54m得hm/h=0.0668/32,54=0.00205,查表得r/rm=2.23，则r=0.0928, rm=2.230.0928=0.2069m。2）再按下式计算风壅水高度e：3）安全超高hc：按级坝和非常运用情况查表得hc=0.4m。正常运用情况下的坝顶超高为h=0.2069+0.4+0.0000017=0.6069m1.4.3主坝坝顶高程复核：坝顶高程等于水库静水位与坝顶超高之和，取一下最大值设计洪水位加正常运行条件下坝顶超高： 691.66+1.023=692.680m正常蓄水位加正常运行条件下坝顶超高： 689.70+1.023=690.723m校核洪水位加非常运行条件下坝顶超高： 692.54+0.6069=693.147m正常蓄水位加非常运行条件下坝顶超高： 689.70+0.6069=690.307m所以，坝顶高程取693.147坝实际高程694.20m和副坝实际高程693.66m，主副坝顶高程均符合要求，不需要再进行整治。1.4.4计算副坝正常运用情况下的坝顶超高1) 计算波浪爬高r：波浪爬高按蒲田试验站公式计算。先计算平均爬高rm,再计算设计爬高r。平均爬高按下式计算： 式中k斜坡糙率渗透性系数，其中砌石取.0750.8，这里取0.75kv经验系数，与v/gh的无量纲大小有关，v/gh=10.5/9.824.18=0.682,查表得kv=1m边坡系数 取3.0hm平均波高 =0.104mlm平均波长 lm=2:5hm=250.104=2.6m因此平均波浪爬高为：由于hm=0.104m，h=34.2m，hm/h=0.00304，三级大坝采用累积频率为1%的爬高值，查表得r/rm=2.23，所以得r=0.274m2) 计算风壅水面高度e： 式中： k为综合摩阻系数,取k=3.610-6；为计算风速，该坝为3级坝，故计算风速取多年平均风速的1.5倍，多年平均风速为1.0m/s，故v0=1.01.5=1.5m/sd为风区长度，取400m；为坝前水域的平均水深，hm=691.66-679.72=11.94；为计算风向与坝轴线的法线间的夹角, 由于风的方向垂直坝面取0则本设计的最大风壅高度为：3) 安全超高hc：按iii级坝和正常运用情况查表得hc=0.7m。综上所述，正常运用情况下的坝顶超高为 h=0.274+0.7+0.0000138=0.974m1.4.5 计算副坝非常运用情况下的坝顶超高非常运用情况下的坝顶超高计算公式与正常情况下一样，所不同的是风速采用多年平均最大风速v=1m/s，坝前水深h=692.54-660=32.54m 1) 计算波浪爬高r：k=0.750.80，采用0.75；取g=9.8m/s2，求得无维量，查表得经验系数kw=1；取风向与坝轴垂线的夹角为0，折减系数kb=1；初拟坝坡m=3.0；吹程d=400m，则平均波高为： 平均波长为 将以上各值代入rm式得根据hm=0.0761m、h=32.54m得hm/h=0.0761/32,54=0.00234,查表得r/rm=2.23，则r=0.0902, rm=2.230.0928=0.201m。2）再按下式计算风壅水高度e：3）安全超高hc：按级坝和非常运用情况查表得hc=0.4m。正常运用情况下的坝顶超高为 h=0.2011+0.4+0.00000613=0.6011m1.4.6副坝坝顶高程复核 表1.4.6 副坝坝顶高程计算成果表p(%)z洪(m)e(m)r(m)a(m)y(m)所需（m）现z坝顶（m）5.12地震沉降后高程（m）1.0691.661.41050.2740.70.974692.634693.66 693.410.1692.546.11060.2010.40.601693.141由表1.4.6知，总岗山水库副坝现所需坝顶高程为693.141m，副坝实际坝顶高程为693.41m，满足国家防洪标准（gb50201-94）防洪要求。2渗流安全评价2.1主坝渗流计算分析2.1.1目的（1）确定坝体浸润线在正常水位情况下的位置和逸出点的高度，为校核坝坡稳定计算提供必须的资料。（2）计算坝体的渗流量。（3）计算坝体粘土的渗透坡降j，确定渗流稳定状态。2.1.2渗透系数主坝为均质土坝，坝体土渗透系数采用地勘报告中的建议值，具体数值见表2.1.2。表2.1.2 主坝渗透系数表 单位：cm/s部位渗透系数坝体土2.5310-4坝基土6.610-4抛石（夹粉土）3.510-32.2主坝渗流运行分析2.2.1渗流分析的方法坝体渗流计算1）土石坝渗流计算的目的1、确定坝体浸润线及其下游出逸点；2、确定坝体与坝基的渗流量；3、确定坝坡出逸段与下游坝基表面的出逸坡降；4、确定水库水位降落时上下游坝坡内的浸润线位置或孔隙压力；5、确定岸坡绕渗的渗流量和渗透比降。2）渗流计算的基本假定1、由于地基的渗透系数很小，所以可以视为不透水地基层；2、坝体为粘性土取其渗透系数k2.5310-4cm/s，可假定上游浸润线位置与水库相同，下游浸润线高度为h0。 3、下游设有棱体排水，且下游无水，因此可假设浸润线为以排水起点为焦点且通过水面与y轴交点的抛物线。4、土体中渗流流速不大且处于层流状态，渗流服从达西定律。3）计算工况根据水库的水位特征值和运行情况，对正常水位、设计洪水位及校核洪水位情况下做了稳定渗流计算和非稳定渗流计算，计算工况如下：工况一：正常水位689.7m 形成的稳定渗流期；工况二：上游设计水位691.66m 形成的渗流情况；工况三：上游校核水位692.54m形成的渗流情况；工况四：正常水位689.7m降落至死水位时形成的非稳定渗流情况；工况五：设计水位691.66m降落至死水位时形成的非稳定渗流情况；工况六：校核水位692.54m降落至死水位时形成的非稳定渗流情况。4 ）渗流分析的方法采用水力学法进行土石坝渗流计算，应用达西定律和杜平假设，建立各段的相应的方程式，然后根据水流的连续性求解渗透流量和浸润线等。1、单宽渗流计算公式：式中： q 单位宽度渗流量;坝体渗透系数，0.017m/d；l计算长度；h0排水起点处浸润线的高度；h1坝前水深：正常水位时为689.7660=29.7m；设计水位时为691.66660=31.66m；校核水位时为692.54660=32.54m；m1上游坝坡系数，取3.0；2.2.2计算过程稳定渗流情况：工况一：计算简图如下 由公式可得：l=12.729m 由图可知 ： l=72.425m浸润线起点x坐标 =89.1(m)浸润线终点x坐标 =163.962(m)所以 h0=1.192m q=1.718(m3/d.m)将数据代入浸润线方程按该式绘制浸润线：在96.343163.962之间假定一系列的x值，求得一系列的y值，列于下表，将表中的x、y值绘于断面图上，并连成曲线便为土坝浸润线。表2.2.2-1正常水位时土坝浸润线计算成果表x(m)96.34103.5110.8118.0125.3132.5139.7147.0154.2161.5163.9y(m)20.6619.0417.3415.5513.6511.629.4237.0154.3171.1920浸润线如下图：工况二：计算简图如下由公式可得：l=13.586m 由图可知 ： l=66.425m浸润线起点x坐标 =95.100(m)浸润线终点x坐标 =163.967(m)所以 h0=1.378m q=1.986(m3/d.m)将数据代入浸润线方程按该式绘制浸润线：假定一系列的x值，求得一系列的y值，将表中的x、y值绘于断面图上，并连成曲线便为土坝浸润线。x(m)101.7108.3115.0121.6128.3134.9141.5148.2154.8161.5163.6y(m)21.6519.9718.2116.3514.3812.279.9857.4754.6571.3780表2.2.2-1设计水位时土坝浸润线计算成果表浸润线如下图：工况三：计算简图如下由公式可得：l=13.929m 由图可知 ： l=64.025m浸润线起点x坐标 =97.5(m)浸润线终点x坐标 =163.969(m)所以 h0=1.459m q=2.104(m3/d.m)将数据代入浸润线方程 按该式绘制浸润线：在97.5-163.969之间假定一系列的x值，求得一系列的y值，将表中的x、y值绘于断面图上，并连成曲线便为土坝浸润线。表2.2.2-3 校核水位时土坝浸润线计算成果表x(m)103.9110.3 116.7123.1129.5135.9142.3148.7155.1161.5163.9y(m)22.0220.3218.5416.6614.66 12.52 10.20 7.6564.7941.459 0浸润线如下图：非稳定渗流情况：对于非稳定渗流的情况，可根据k/（v)的大小分情况考虑：若k/（v)0.1，库水位骤降，按原浸润线计算；若0.1k/（v)60，属于极缓慢下降，按最终的浸润线计算。其中：k为大坝的渗透系数；为给水度，=0.117；v为水位下降的速度。非稳定渗流计算时，上下游坡率均采用平均坡率，上游为3.0，下游为2.25，计算结果如下：工况四：查水位库容曲线可得：正常水位689.7m对应的库容为2588万方 死水位666.7m对应的库容为0万方由于溢洪道最大放水流量为240.1m/s，所以水库水位降落所用的时间为：2588*104/240.1*24*60*60=1.248d，降落速度为 : （689.7-666.7）/1.248=18.43m/d, 由于=0.117k1/7=0.094，所以k/（v)=0.219/(18.43*0.094)=0.126,属于缓降，需计算各时段的浸润线。简图如下：水位下降所需时间 = 1.25(天) 时间(天) 浸润线最高点 上游坡出渗点高度he(t)(m) 0.25 28.850 23.621时间(天) 浸润线最高点h0(t)+h0(m) 上游坡出渗点高度he(t)(m) 0.50 27.999 20.930时间(天) 浸润线最高点h0(t)+h0(m) 上游坡出渗点高度he(t)(m) 0.75 27.149 18.938时间(天) 浸润线最高点h0(t)+h0(m) 上游坡出渗点高度he(t)(m) 1.00 26.299 17.323时间(天) 浸润线最高点h0(t)+h0(m) 上游坡出渗点高度he(t)(m)1.25 25.449 15.957各时段浸润线如下：表2.2.2-4 正常水位降至死水位时土坝浸润线计算成果表t1=0.25dt2=0.5dt3=0.75dt4=1dt5=1.25dx(m)y(m)x(m)y(m)x(m)y(m)x(m)y(m)x(m)y(m)70.863 24.195 62.791 20.930 56.815 18.938 51.968 17.323 47.872 15.957 72.687 23.621 65.422 21.741 60.043 27.149 55.681 18.418 51.995 17.145 74.511 24.755 68.053 22.522 63.272 20.841 59.394 19.452 56.118 18.255 76.334 25.303 70.684 23.278 66.500 21.730 77.960 20.434 60.241 19.301 78.158 25.840 73.314 24.009 69.729 22.584 66.821 21.370 64.363 20.294 81.805 26.365 75.945 24.719 89.100 23.407 70.534 22.268 68.486 21.240 83.629 28.850 78.576 25.409 76.186 24.201 74.247 26.299 76.732 22.146 85.453 27.386 81.207 26.081 79.414 24.971 89.100 23.130 80.854 25.449 87.276 27.883 86.469 27.999 82.643 25.718 81.674 24.765 84.977 23.854 89.100 28.370 89.100 27.375 85.871 26.443 85.387 25.544 89.100 24.664 79.982 26.881 83.838 26.736 72.957 19.912 63.107 23.962 72.609 23.016 结果如下图：工况五：查水位库容曲线可得：正常水位689.7m对应的库容为2588万方 设计水位691.66m对应的库容为3050万方死水位666.7m对应的库容为0万方由于溢洪道最大放水流量为240.1m/s，所以水库水位降落所用的时间为：3050*104/240.1*24*60*69=1.471d，降落速度为 : （691.7-666.7）/1.471=16.96m/d, 由于=0.117k1/7=0.094，所以k/（v)=0.219/(16.96*0.094)=0.137,属于缓降，需计算各时段的浸润线。计算简图如下： 水位下降所需时间 = 1.47(天) 时间(天) 浸润线最高点h0(t)+h0(m) 上游坡出渗点高度he(t)(m) 0.29 30.756 24.952 时间(天) 浸润线最高点h0(t)+h0(m) 上游坡出渗点高度he(t)(m) 0.59 29.813 22.008 时间(天) 浸润线最高点h0(t)+h0(m) 上游坡出渗点高度he(t)(m) 0.88 28.869 19.829 时间(天) 浸润线最高点h0(t)+h0(m) 上游坡出渗点高度he(t)(m) 1.18 27.925 18.063 时间(天) 浸润线最高点h0(t)+h0(m) 上游坡出渗点高度he(t)(m) 1.47 26.982 16.571各时段浸润线如下表：表2.2.25 设计水位降至死水位时土坝浸润线计算成果表t1=0.29dt2=0.59dt3=0.88dt4=1.18dt5=1.47dx(m)y(m)x(m)y(m)x(m)y(m)x(m)y(m)x(m)y(m)74.856 24.952 66.023 22.008 59.487 19.829 54.189 18.063 49.714 16.571 76.880 25.592 68.931 22.908 63.049 20.910 58.280 27.925 54.253 17.887 80.929 26.216 71.839 29.813 66.610 21.937 62.371 20.420 81.484 19.113 82.954 26.825 74.746 24.611 70.171 22.919 66.462 21.502 63.330 20.264 84.978 30.756 77.654 25.419 73.732 28.869 95.100 22.532 67.868 21.354 87.002 28.005 95.100 26.202 77.294 24.765 70.553 23.517 72.407 26.982 89.027 28.576 83.469 26.963 80.855 25.638 74.644 24.462 76.946 22.390 91.051 29.137 86.377 27.703 84.416 26.483 78.735 25.372 95.100 23.381 93.076 29.686 89.285 28.424 87.977 27.302 82.827 26.251 86.023 24.331 95.100 30.226 92.192 29.127 91.539 28.096 86.918 27.101 90.561 25.246 78.905 27.422 80.562 23.775 95.100 23.860 91.009 19.278 58.791 26.128 结果如下图：工况六：查水位库容曲线可得：正常水位689.7m对应的库容为2588万方 校核水位692.54m对应的库容为3260万方死水位666.7m对应的库容为0万方由于溢洪道最大放水流量为240.1m/s，所以水库水位降落所用的时间为：3260*104/240.1*24*60*69=1.572d，降落速度为 : （692.5-666.7）/1.471=16.49m/d, 由于=0.117k1/7=0.094，所以k/（v)=0.219/(16.49*0.094)=0.141,属于缓降，需计算各时段的浸润线。计算简图如下：水位下降所需时间 = 1.57(天) 时间(天) 浸润线最高点h0(t)+h0(m) 上游坡出渗点高度he(t)(m) 0.31 31.519 25.484 时间(天) 浸润线最高点h0(t)+h0(m) 上游坡出渗点高度he(t)(m) 0.63 30.537 22.438 时间(天) 浸润线最高点h0(t)+h0(m) 上游坡出渗点高度he(t)(m) 0.94 29.556 20.185时间(天) 浸润线最高点h0(t)+h0(m) 上游坡出渗点高度he(t)(m) 1.26 28.574 18.358时间(天) 浸润线最高点h0(t)+h0(m) 上游坡出渗点高度he(t)(m) 1.57 27.593 16.816 各时段浸润线如下：表2.2.2-6 校核水位降至死水位时土坝浸润线计算成果表t1=0.31dt2=0.63dt3=0.94dt4=1.26dt5=1.57dx(m)y(m)x(m)y(m)x(m)x(m)y(m)x(m)y(m)x(m)76.453 25.484 67.315 22.438 60.555 76.453 25.484 67.315 22.438 60.555 78.557 26.800 70.333 23.375 64.249