调压井课程设计.doc

水电站课程设计计算书水电站课程设计计算书 目录一、设计课题3二、设计资料及要求31、设计资料见课程设计指导书、任务书32、设计要求3三、调压井稳定断面的计算31、引水道的水头损失计算3（1）局部水头损失计算3（2）沿程水头损失计算42、引水道的等效断面面积计算63、调压井稳定断面计算7四、调压井水位波动计算81、最高涌波水位计算81）、当丢弃负荷：300000kw时，采用数解法82）、当丢弃负荷为4500015000时，采用图解法：92、最低涌波水位101) 丢弃负荷度为300000KW时（数解法）102）增加负荷度为30000-45000KW时（两种方法）11五调节保证计算131、检验正常工作情况下的水击压力132、检验相对转速升高是否满足规范要求16六、参考文献16七、附图:17附图1:丢弃负荷时调压井水位波动图17附图2:增加负荷时调压井水位波动图17一、设计课题：水电站有压引水系统水力计算。二、设计资料及要求1、设计资料见课程设计指导书、任务书；2、设计要求：（1）对整个引水系统进行水头损失计算；（2）进行调压井水力计算求稳定断面；（3）确定调压井波动振幅，包括最高涌波水位和最低涌波水位；（4）进行机组调节保证计算，检验正常工作状况下水击压力、转速相对值升高是否满足规范要求。三、调压井稳定断面的计算1、引水道的水头损失计算（1）局部水头损失计算表局部水头损失采用如下公式计算：表1局部水头损失计算表栏号引水建筑物部位及运行工况断面面积（）局部水头损失系数局部水头损失 (m)合计（m）1进水口拦污栅61.280.121.630.3072进口喇叭段29.760.105.763闸门井24.000.2017.724渐变段23.880.054.475隧洞进口平面转弯23.760.076.330.2036末端锥管段19.630.1013.247调压井正常运行19.630.1013.242.2048增一台负荷，从调压井流入管道19.631.50198.619 压力引水管道上水平段平面转弯19.630.045.304.46710下水平段平面转弯9.080.0849.5111斜井顶部立面转弯19.630.0911.9212斜井底部立面转弯9.080.0955.6913锥管9.080.0849.5114三台机满发1叉管19.630.3039.7215三台机满发2叉管19.630.4559.5816一台机满发1叉管19.630.2735.7517一台机满发2叉管19.630.2735.7518蝴蝶阀9.080.1486.64本表计算中Q=102m3/s，g=9.8m/s2（2）沿程水头损失h程计算表沿程水头损失采用如下公式计算： 表2 沿程水头损失计算表栏号引水道部位过水断面面积W（）湿周（m）水力半径R（m）引水道长（m）h程 (m)合计（m）1进水口喇叭口进水段29.7621.981.35406.04522.6590.0072闸门井段24.0020.001.20005.67624.1340.0133渐变段23.8818.641.281110.012603.3900.0224隧洞D=5.5M段23.7617.281.3756469.6543724.60.8155锥形洞段21.6516.491.31255.07423.1540.0116调压井前管段19.6315.711.250010.9821161.590.03771号叉管19.6315.711.250035.7468881.180.10981-2号叉管19.6315.711.250029.2156296.010.04092号叉管19.6315.711.250012.2323570.70.00410锥管段13.8513.191.05003.9719392.710.00311D=3.4段9.0810.680.850021.25320107.10.056其中 栏1、2 、3、4、5、6、7的流量Q为102m3/s，根据压力管道相关参数表得7栏的流量为96.9,；8栏的流量为64.6,； 9、10、11栏流量为32.3 查规范和资料得到糙率n，进水口取0.013，隧洞取最小值0.012，压力管道取最大值0.013 调压井前引水道的水头损失压力管道的水头损失（压力管道长度为113.3m,较长不计局部水头损失）整个引水系统的水头损失2、引水道的等效断面面积计算其中L 为调压井前引水道的长度L=拦污栅长度+喇叭口进口段长度+闸门井段长度+渐变段长度+（D=5.5M洞段长度）+锥形洞段长度+调压井前管段长度=4.1+6.0+5.6+10.0+469.6+5.0+10.98=511.28m 栏号引水道部位过水断面fi（m2）Li(m)(1)拦污栅61.284.10.067(2)喇叭口进水段29.766.00.202(3)闸门井段24.005.60.233(4)渐变段23.8810.00.419(5)D=5.5m23.76469.619.764(6)锥形洞段21.655.00.231(7)调压井前管段19.6310.980.559计算表引水道的等效断面面积： m23、调压井稳定断面计算为使求得的稳定断面满足各种运行工况的要求，上游取死水位，下游取正常尾水位情况计算 净水头H0=上游死水位下游正常尾水位=1082.01028.5=53.5m：引水道水头损失，大小为1.415hwT0：压力管道沿程水头损失，大小为0.212m=53.5-1.415-30.212=51.449m当三台机组满出力时，保证波动稳定所需的最小断面：=k其中K的取值为1.01.1，为引水道总阻力系数D=5.5m =取k=1.0则保证稳定所需要的最小断面为：四、调压井水位波动计算1、最高涌波水位计算（1）当丢弃负荷：300000kw时，采用数解法当上游为校核洪水位1097.35m，下游为相应的尾水位1041.32m，电站丢弃两台机时，若丢荷幅度为300000KW，则流量为63.60m3/s，用数解法计算。L-为引水道的长度为511.28mf-引水道等效断面面积v0-引水道水流流速v0=m/sF-调压井稳定断面为156.46m2-引水道水头损失（=）g-取9.8m/s2(1.63+5.76+17.72+4.47+6.33+13.24)10-663.663.6=0.20m=(4522.659+7624.134+12603.390+543724.6+7423.154+21161.59) 0.0120.01263.663.610-6=0.348m=0.20+0.348=0.548m查书本P150图10-4得，则=0.1051.6+1097.35=1112.51m（2）当丢弃负荷为4500015000kw时，采用图解法：当上游为校核洪水位1097.35m，下游为相应的尾水位1041.32m，电站丢弃两台机组时，若丢荷幅度为4500015000KW，则流量为96.531.0m3/s。利用图解法求解1、以横轴表示引水道流速v，以圆点向左为正（水流向调压室），向右为负；以纵轴表示水位z，以向上为正，向下为负，横轴相当于静水2、作辅助线曲线引水道水头损失曲线：，C=87.882 =0.232 绘制曲线计算时段取值范围为,的取值范围为3.94.6选取=4s=,当丢弃负荷为45000kw15000kw时，流量96.531.0，流速4.051.30m/s。=0.20=0.1524-0.204=0.608-0.8绘制曲线 =0.077（）采用matlab编程计算后画图，源代码如下：v(1)=4.05-0.077*(0.608*4.05-0.8)z(1)=-0.232*4.05*4.05+0.608*4.05-0.8for i=1:29 dz(i)=0.608*v(i)-0.8 dv(i)=0.077*(-1)*z(i)-0.232*v(i)*v(i) v(i+1)=v(i)+dv(i) z(i+1)=z(i)+dz(i)endr(1,:)=zr(2,:)=vr=r其中，v(i)为流速矩阵，z（i）为水位壅高矩阵，dz（i）为水位壅高增量矩阵，dv（i）为流速矩阵增量矩阵,v(1)为第一时段末的水的流速,z(1)为第一时段末调压井内水位的壅高，第二个以后时段的水位及流速如下表所示。 V(m/s)Z(m)VZ3.8850-2.1430-0.10461.56213.7804-0.5809-0.21061.49853.56980.9176-0.29831.37043.27152.2880-0.36741.18912.90413.4771-0.41840.96572.48574.4428-0.45250.71132.03335.1541-0.47070.43621.56255.5903-0.47410.15001.08855.7404由表可知，最大壅高水位在5.59m5.74m之间，线性内插得最大壅高水位为5.67（图纸见附图1）。 Zmax=1097.35+5.67=1103.02m2、最低涌波水位：（1）丢弃负荷度为300000KW时（数解法）当上游为死水位，下游为正常尾水位时，若电站丢弃全负荷时（300000，流量变化为67.50），因调压室水位达到最高水位时，水位开始下降，此时隧洞中的水流朝着水库方向流动，水从调压室流向进水口，因此水头损失应变为负值，水位到达最低值称为第二振幅。，=2.835m/s (1.63+5.76+17.72+4.47+6.33+13.24)10-667.567.5=0.220m=(4522.659+7624.134+12603.390+543724.6+7423.154+21161.59) 0.0120.01267.567.510-6=0.391m=0.220+0.391=0.611m=0.012查书本P150图10-4得=0.08，则=0.08(-52.21)=-4.177mZmin=1082-4.177=1077.823m（2）增加负荷度为30000-45000KW时（两种方法）当上游为死水位，下游为正常尾水位，增荷幅度为3000045000KW，流量变化由68.5102.5，流速2.884.3 m/s。 A、 数解法 m=68.5/102.5=0.668(1.63+5.76+17.72+4.47+6.33+13.24)10-6102.5102.5=0.516m=(4522.659+7624.134+12603.390+543724.6+7423.154+21161.59) 0.0120.012102.5102.510-6=0.902m=0.516+0.902=1.418m =) =3.435 调压井的最低水位为 1082-4.871=1077.129mB、 图解法增加负荷时的图解法与丢弃负荷的图解法类似，同样选择坐标系，绘出、引水道水头损失曲线：，C=87.882 =0.232、绘出曲线；则计算时段取值范围为取，的取值范围为3.94.6选取=4s又=,当增加负荷为3000045000kw时，流量68.5102.5m3/s。=0.655=0.1524-0.6554=0.608-2.62、绘制曲线=0.077（）Matlab编写程序：v(1)=2.88z(1)=-0.232*2.88*2.88for i=1:29 dz(i)=0.608*v(i)-2.62 dv(i)=0.077*(-1)*z(i)-0.232*v(i)*v(i) v(i+1)=v(i)+dv(i) z(i+1)=z(i)+dz(i)endv=vz=zr(:,1)=vr(:,2)=z计算结果如下表V(m/s)Z(m)VZ2.946-2.7930.0599-0.82833.0069-3.62150.11730.79183.1242-4.41340.1655-0.72053.2897-5.13390.2020-0.61993.4917-5.75370.2252-0.49713.7169-6.25080.2345-0.36013.9514-6.61090.2301-0.21754.1815-6.82850.2134-0.07764.3950-6.9061 最大下降水位在-6.83m-6.90m之间,采用线性内插法,得到最大下降水位为6.87m。Zmax=1082-6.87=1075.13m。五调节保证计算1、检验正常工作情况下的水击压力相应的正常工作情况下的正常蓄水位为1092.0米，相应的尾水位为1028.5米，则：，由所给资料可知：正常工作情况下的正常蓄水位为1092.0米，相应的尾水位为1028.5米，三台机满发电，通过水轮机的流量为96.9，于是：起始的流速96.9/23.808=4.07m,水锤波速：上式中为水的体积弹性模量，一般为2.1；为水的容重，取9.81KN/m3，r为管道的半径,压力管道半径为2.5m.蜗壳半径为1.22m.尾水管半径为1.7m，K为管道抗力系数。1.4压力管道=102.32米，蜗壳=20.4米，尾水管=16.2米；，V1=；V2 =；V3= ；V5 = =977.073m/sVm=管道特性系数为当机组满负荷运行时，因为阀门从全开到全关的时间为7s，其中有效关闭时间为,一个相长4.68s,为间接水锤。管道的特性系数： 由以上计算并查图9-4水锤类型判别图可知：该水锤为极限正水锤，由式9-15可得：校核尾水管进口处的真空度尾水管进水口处的水锤压力：；尾水管在进口断面出现时的流速水头如下： 尾水管进口处的真空度在89米之内，满足要求。2、检验相对转速升高是否满足规范要求该水电站的保证出力为1.07104kw，水轮机型为HL211-LJ-225，对于混流式水轮机，f 为水锤影响系数，根据管道特性曲线系数由图9-15得出f=1.12 根据书102页转速变化率计算标准：当机组容量占电力系统总容量的比重较大，且担负调频任务时， max宜小于45；当机组容量占电力系统总容量的比重不大或担负基荷时，max宜小于55；对斗叶式水轮机，max宜小于30,所以，满足要求。六、参考文献水电站建筑物第二版 马善定，汪如泽水工设计手册（水电站建筑物） 华东水利学院水电站建筑物课程设计任务指导书 青海大学水利电力学院水电站调压室设计规范DL/T5058-1996工程制图及CAD七、附图:附图1:丢弃负荷时调压井水位波动图附图2:增加负荷时调压井水位波动图18