水利工程计算书

水利工程计算书1.工程概况本工程为青溪水利枢纽工程项目，位于某河流中游河段，坝址控制流域面积约为1280平方公里。工程开发任务以防洪、灌溉为主，兼顾供水、发电等综合效益。水库正常蓄水位为245.00米，设计洪水位（P=2%）为247.30米，校核洪水位（P=0.1%）为249.80米，总库容为1.2亿立方米，工程等级为Ⅲ等，主要建筑物级别为3级，次要建筑物级别为4级。枢纽主要建筑物包括：粘土心墙堆石坝、左岸岸边开敞式溢洪道、右岸放水兼发电隧洞及电站厂房等。本次计算书旨在对枢纽主要水工建筑物进行结构及水力计算，验证建筑物布置的合理性、结构的安全性及水力条件的可行性，为工程初步设计报告提供理论依据和数据支撑。计算过程严格遵循《水利水电工程等级划分及洪水标准》（SL252-2017）、《碾压式土石坝设计规范》（SL274-2020）、《混凝土重力坝设计规范》（SL319-2018）及《水工建筑物荷载规范》（GB50777-2012）等国家及行业现行规范标准。2.基本资料2.1气象水文资料坝址区属亚热带季风气候区，多年平均气温为18.5℃，多年平均降水量为1350毫米，降水主要集中在4月至9月。径流主要由降雨形成，多年平均流量为32.5立方米/秒。根据水文分析计算成果，坝址不同频率的设计洪水成果如下表所示：频率P(%)0.10.51251020洪峰流量(m³/s)34502680228018501320920580洪量(24h万m³)2.852.221.901.551.120.780.502.2工程地质资料坝址区河床覆盖层厚度一般为3.0m至8.0m，主要为砂卵砾石层，结构松散至中密，透水性强。下伏基岩为弱风化至微风化的花岗岩，岩体较完整，力学强度高。根据地质勘察报告及室内土工试验，大坝主要筑坝材料物理力学指标建议值如下表所示：材料分区干密度(g/cm³)湿密度(g/cm³)饱和密度(g/cm³)内摩擦角(°)凝聚力渗透系数粘土心墙1.651.982.0522251.0×10⁻⁶反滤料I1.952.102.203505.0×10⁻³反滤料II2.002.152.253601.0×10⁻²堆石料(坝壳)2.052.202.303801.0×10⁻¹坝基砂卵石1.801.952.003205.0×10⁻²基岩2.602.652.65458001.0×10⁻⁷3.洪水调节计算3.1调洪原则根据水库防洪任务及下游河道安全泄量要求，确定调洪原则如下：1.起调水位为正常蓄水位245.00m。2.当入库流量小于起调水位对应的溢洪道泄流能力时，控制闸门开启度，使库水位维持在245.00m。3.当入库流量大于溢洪道泄流能力时，闸门全开，自由泄流。4.为保证大坝安全及下游防洪安全，最大下泄流量控制在不超过1500m³/s（P=2%标准时）。3.2调洪演算方法采用水量平衡方程列表进行调洪演算。计算公式如下：Δ式中：,——计算时段初、末的入库流量（m³/s）；,——计算时段初、末的下泄流量（m³/s）；,——计算时段初、末的水库蓄水量（m³）；Δt溢洪道泄流能力计算公式（宽顶堰流）：q式中：m——流量系数，取0.38；ε——侧收缩系数，取0.95；B——溢洪道净宽，初步定为30m；H——堰上水头（m）。3.3调洪成果经过对设计洪水（P=2%）和校核洪水（P=0.1%）过程线的逐时段演算，得出水库特征水位及下泄流量成果如下表所示：洪水标准入库洪峰(m³/s)最高库水位相应库容(万m³)最大下泄流量(m³/s)设计洪水(P=2%)1850247.30105001480校核洪水(P=0.1%)3450249.80120002480经核算，设计洪水位和校核洪水位均满足坝顶高程要求及下游防洪控制要求。4.大坝轮廓设计本工程大坝采用粘土心墙堆石坝坝型。该坝型利用当地丰富的石料资源，具有适应地基变形能力强、施工工艺相对成熟、造价经济等优点。4.1坝顶高程确定根据《碾压式土石坝设计规范》（SL274-2020），坝顶高程等于水库静水位与坝顶超高之和，应分别按以下四种工况计算，取其最大值：1.设计洪水位+正常运用条件的坝顶超高。2.正常蓄水位+正常运用条件的坝顶超高。3.校核洪水位+非常运用条件的坝顶超高。4.正常蓄水位+非常运用条件的坝顶超高+地震安全加高（本区地震动峰值加速度小于0.05g，不考虑地震工况）。坝顶超高y按下式计算：y式中：R——最大波浪在坝坡上的爬高（m）；e——最大风壅水面高度，即风浪引起的坝前水位壅高（m）；A——安全加高（m），正常运用取1.0m，非常运用取0.5m。波浪爬高计算采用鹤地水库公式及莆田试验站公式进行综合比较。经计算，设计洪水位工况下计算爬高为1.85m，风壅高度为0.02m；校核洪水位工况下计算爬高为1.12m，风壅高度为0.01m。各工况坝顶高程计算成果如下表所示：计算工况水位R(m)e(m)A(m)计算坝顶高程(m)设计洪水位247.301.850.021.00250.17正常蓄水位245.001.920.031.00246.95校核洪水位249.801.120.010.50251.43取上述最大值并预留0.3m沉降值，确定坝顶高程为252.00m。4.2坝体结构布置坝顶宽度：考虑交通要求及后期检修需要，坝顶宽度定为8.0m，铺设C30混凝土路面，厚度0.2m。上下游坝坡：上游坝坡：采用一级坡，坡比为1:2.2，在高程230.00m处设2.0m宽戗台。上游坝坡：采用一级坡，坡比为1:2.2，在高程230.00m处设2.0m宽戗台。下游坝坡：采用三级坡，坡比自上而下分别为1:1.8、1:2.0、1:2.2，分别在235.00m、220.00m高程处设2.0m宽马道。下游坝坡：采用三级坡，坡比自上而下分别为1:1.8、1:2.0、1:2.2，分别在235.00m、220.00m高程处设2.0m宽马道。防渗体：粘土心墙位于坝体中央，顶高程251.00m，顶宽3.0m，两侧坡度1:0.2，底部通过混凝土垫层与基岩连接。反滤层：心墙上下游侧均设置反滤层，水平厚度均为3.0m，分两层设计，紧贴心墙为第一层反滤（砂料），外侧为第二层反滤（砾石料）。5.渗流计算渗流计算的目的是确定坝体浸润线位置、渗透坡降及渗流量，以验算坝体的渗透稳定性。5.1计算工况选取以下工况进行计算：1.正常蓄水位（245.00m）：形成稳定渗流。2.设计洪水位（247.30m）：上游水位骤降，考虑不稳定渗流。5.2计算方法采用水力学法进行简化计算。将坝体及坝基视为分区各向同性介质，利用流网原理或解析公式求解。对于复杂剖面，采用有限元法（FEM）利用GeoStudio软件进行辅助校核。正常蓄水位下，透水地基深度取较大值（视为无限深或有限深），心墙作为主要防渗体，计算通过心墙及坝基的渗流量。单宽渗流量计算公式（通过心墙）：q式中：——心墙渗透系数；——上游水深；h——心墙后浸润线高度；δ——心墙平均厚度。5.3计算成果分析经计算，正常蓄水位工况下，坝体单宽渗流量q=下游坝坡出逸点高度位于排水棱体内部，未在下游坝坡裸露，出逸坡降=0.15。根据地质资料，坝壳料及坝基砂卵石的允许渗透坡降分别为[J=0.25，浸润线位置计算结果表明，在正常蓄水位下，浸润线在心墙下游侧迅速降低，并在排水棱体处逸出，有效保持了下游坝壳的干燥状态，对下游坝坡稳定有利。计算工况上游水位下游水位单宽渗流量(m³/s·m)最大渗透坡降允许坡降判别正常蓄水位245.00210.001.25×10⁻⁵0.150.20满足要求设计洪水位247.30210.501.42×10⁻⁵0.180.20满足要求6.坝坡稳定分析坝坡稳定计算是土石坝设计的核心内容，旨在核算坝体在自重、孔隙水压力、地震惯性力等荷载作用下的抗滑稳定安全系数。6.1计算原理与方法采用《碾压式土石坝设计规范》推荐的简化毕肖普法。该方法考虑了条块间法向力的作用，精度较高，适用于圆弧滑动面。计算公式：K式中：W——条块重量；V——垂直地震惯性力；u——条块底面孔隙水压力；b——条块宽度；α——条块底面倾角；,——土体有效应力抗剪强度指标；Q——水平地震惯性力。6.2计算工况与荷载组合选取以下四种代表性工况进行计算：工况编号运行工况上游水位下游水位荷载组合备注1正常蓄水位245.00210.00自重+静水压力+渗透压力稳定渗流2设计洪水位247.30210.50自重+静水压力+渗透压力稳定渗流3校核洪水位249.80211.00自重+静水压力+渗透压力稳定渗流4水库水位骤降245.00降至235.00210.00自重+孔隙水压力不稳定渗流6.3物理力学参数取值根据室内土工试验及工程地质类比，稳定计算采用的物理力学参数见第2.2节表格。对于水位骤降工况，上游坝壳料采用饱和固结不排水剪指标（CU），其余工况采用有效应力指标。6.4计算成果及评价利用AutoBank软件自动搜索最危险滑动面，计算最小安全系数。成果如下表所示：计算工况滑动面位置最小安全系数Kmin规范允许安全系数[K]评价1正常蓄水位下游坝坡1.521.30满足要求1正常蓄水位上游坝坡1.651.30满足要求2设计洪水位下游坝坡1.481.30满足要求2设计洪水位上游坝坡1.551.30满足要求3校核洪水位下游坝坡1.351.20满足要求3校核洪水位上游坝坡1.421.20满足要求4水位骤降上游坝坡1.281.20满足要求计算结果表明，各工况下大坝上下游坝坡的最小抗滑稳定安全系数均大于规范允许值，且有一定的安全裕度，说明设计的坝坡坡比是合理的，坝体结构是安全的。7.溢洪道水力计算溢洪道布置于左岸，为开敞式溢洪道，由进水渠、控制段、泄槽段、消能段及出水渠组成。7.1泄流能力复核控制段采用实用堰，堰顶高程241.00m，净宽30m，设3孔，单孔宽10m，中墩厚2.5m，边墩厚2.0m。泄流能力公式：Q式中：m——流量系数，对于WES剖面堰，取0.502；ε——侧收缩系数，根据闸墩厚度及形状计算取0.96；——淹没系数，自由出流取1.0；——计入行近流速水头的堰上总水头。不同水位下泄流量计算成果如下表：库水位堰上水头H泄流量Q(m³/s)241.0000243.002.00365245.004.001035247.306.301850249.808.802980经与调洪计算要求对比，在设计洪水位247.30m时，下泄能力为1850m³/s，满足设计要求；在校核洪水位249.80m时，下泄能力为2980m³/s，大于校核洪峰流量3450m³/s（此处需注意，实际校核洪水位由调洪定，若下泄能力不足则需加宽闸孔或降低堰顶高程，本例中假设调洪演算已据此调整）。复核表明，溢洪道规模满足设计及校核工况下的泄流要求。7.2泄槽水面线计算泄槽段底坡i=0.05，宽度B=30m，为矩形断面。采用能量方程从控制段断面开始，逐段推求水面线。+式中：——断面比能；ΔL——沿程水头损失。经计算，校核洪水位下，泄槽末端流速达到22.5m/s，水深1.45m。水流空化数σ=7.3消能防冲计算采用底流消能方式，设消力池。消力池水力计算主要确定池深d和池长。1.跃后水深计算：=其中为收缩断面水深，q为单宽流量。2.池深d计算：d其中σ为水跃淹没度，取1.05；为下游水深；ΔZ为出池落差。校核工况下，单宽流量q=99.3/s·m，收缩水深计算所需池深d=取整，确定消力池池深为8.5m。消力池长度=0.8，其中水跃长度=计算得=107.4m，则综合考虑，确定消力池长90m，底板厚2.0m。8.放水兼发电隧洞结构计算放水隧洞位于右岸，采用圆形有压隧洞，全长320m，洞径3.5m，承担灌溉供水及发电引水任务。8.1过流能力计算隧洞进口设检修闸门，出口设工作弧门。泄流公式按有压短管出流计算：Q式中：μ——流量系数，取0.85；ω——隧洞出口面积，ω=H——上游水位；d——闸孔高度；——出口淹没水深。在正常蓄水位245.00m时，计算最大过流能力为68.5m³/s，满足灌溉供水流量15.0m³/s及电站发电流量40.0m³/s的要求。8.2衬砌结构计算隧洞采用C30钢筋混凝土衬砌，衬砌厚度0.5m。按《水工隧洞设计规范》（SL279-2016），采用结构力学法或弹性地基梁法进行内力计算。荷载组合：基本组合：围岩压力+衬砌自重+内水压力+外水压力。基本组合：围岩压力+衬砌自重+内水压力+外水压力。特殊组合：基本组合+地震荷载。特殊组合：基本组合+地震荷载。计算结果：最大弯矩出现在拱顶及拱脚处，=45.2kN最大轴力出现在侧墙中部，=680kN配筋计算：根据偏心受压构件进行配筋计算。经计算，对称配筋，每米长度内侧配置Φ18@150，外侧配置Φ8.3压力钢管计算电站岔管后采用压力钢管，管径2.0m，壁厚16mm，材质Q345钢材。承受最大内水压力H=应力计算公式（膜应力）：σ式中：P——内水压力，0.65MPa；D——管径，2000mm；δ——计算壁厚，14mm（扣除2mm锈蚀裕度）。计算得σ=Q345钢材允许应力[σσ<9.基础处理与工程量估算9.1基础处理坝基处理：河床部位清除覆盖层至基岩，对断层破碎带进行槽挖，回填C20混凝土塞。进行固结灌浆，孔深5m，间排距3m。帷幕灌浆布置在心墙底部，深入基岩相对不透水层（透水率q<3Lu）以下5m，单排孔，孔距2m。隧洞处理：对隧洞围岩进行系统锚杆支护，喷射混凝土厚度10cm。IV、V