水锤计算方法.docx

水锤计算方法.docx

《水锤计算方法.docx》由会员分享，可在线阅读，更多相关《水锤计算方法.docx（62页珍藏版）》请在冰豆网上搜索。

水锤计算方法

第一节概述

一、水电站的不稳定工况机组在稳定运行时，水轮机的出力与负荷相互平衡，这时机组转速不变，水电站有压引水系统（压力隧洞、压力管道、蜗壳及尾水管）中水流处于恒定流状态。

在实际运行过程中，电力系统的负荷有时会发生突然变化（如因事故突然丢弃负

荷，或在较短的时间内启动机组或增加负荷），破坏了水轮机与发电机负荷之间的平衡，机组转速就会发生变化。

此时水电站的自动调速器迅速调节导叶开度，改变水轮机的引用流量，使水轮机的出力与发电机负荷达到新的平衡，机组转速恢复到原来的额定转速。

由于负荷的变化而引起导水叶开度、水轮机流量、水电站水头、机组转速的变

化，称为水电站的不稳定工况。

其主要表现为：

（1）引起机组转速的较大变化由于发电机负荷的变化是瞬时发生的，而导叶的启闭需要一定时间，水轮机出力不能及时地发生相应变化，因而破坏了水轮机出力和发电机负荷之间的平衡，导致了机组转速的变化。

丢弃负荷时，水轮机在导叶关闭过程中产生的剩余能量将转化为机组转动部分的动能，从而使机组转速升高。

反之增加负荷时机组转速降低。

（2）在有压引水管道中发生“水锤”现象当水轮机流量发生变化时，管道中的流量和流速也要发生急剧变化，由于水流惯性的影响，流速的突然变化使压力水管、蜗壳及尾水管中的压力随之变化，即产生水锤。

导叶关闭时，在压力管道和蜗壳中将引起压力上升，尾水管中则造成压力下降。

反之导叶开启时，在压力管道和蜗壳内引起压力下降，而在尾水管中引起压力上升。

（3）在无压引水系统（渠道、压力前池）中产生水位波动现象。

无压引水系统中产生的水位波动计算在第八章已介绍。

二、调节保证计算的任务水锤压力和机组转速变化的计算，一般称为调节保证计算。

调节保证计算的任务及目的是：

（1）计算有压引水系统的最大和最小内水压力。

最大内水压力作为设计或校核压力管道、蜗壳和水轮机强度的依据之一；最小内水压力作为压力管道线路布置、防止压力管道中产生负压和校核尾水管内真空度的依据。

（2）计算丢弃负荷和增加负荷时的机组转速变化率，并检验其是否在允许范围内。

（3）选择水轮机调速器合理的调节时间和调节规律，保证压力和转速变化不超过规

定的允许值。

（4）研究减小水锤压力及机组转速变化率的措施。

第二节水锤现象及其传播速度

、水锤现象

在水电站运行过程中，为了适应负荷变化或由于事故原因，而突然启闭水轮机导叶时，由于水流具有较大惯性，进入水轮机的流量迅速改变，流速的突然变化使压力水管、蜗壳及尾水管中的压力随之变化，这种变化是交替升降的一种波动，这种现象称为水锤。

时段，在dx2、dx3、••管段中流速、压力将相继发生同样的变化，见图9-i（a）。

这样，一段接一段地将阀门关闭的影响向上游传播，压力增加如同波一样自阀门A处沿管道逐渐向上游传播，这就是水锤波，其传播速度称之为水锤波速a，变化的压

力厶H称为水锤压力。

使压力增加的波为增压波，使压力降低的波叫降压波。

经过L/a

时间，水锤波达到管道进口处，此时，整个水管内的流速Vo降为零，压力上升为Ho+

△Ho

（2）t=L/a~2L/ao当t=L/a时，水锤波将传至水库点D处，由于D点右端管道内压

力为Ho+厶H,而左端水库保持不变为H,因此“边界”处的水体不能保持平衡，管道

中的水体在△H压差作用下将逆流向水库。

在t=L/a后的dti时段内，首先是紧靠水库

dxn管段内发生变化，流速将由0变为-Vo,压力由Ho+^H变为Ho；管壁及水体随着水

锤压力的消失恢复至原状。

同理接再经dt2、dt3、…时段，在相应dxn-1、dXn-2、dXn-3…

管段中将发生同样的变化，如图9-1（b）。

直到t=2L/a时刻，整个管道中的压力、流速、

管径及水的密度均恢复到初始状态。

这说明，水锤波在水库处要发生反射，反射特点是“等值异号”反射，即反向波与入射波的数值相同，均为H，但符号相反，升压

波反射为降压波。

（3）t=2L/a~3L/a。

当t=2L/a时，水锤波传播到阀门处A点，由于阀门已关闭，力口

之水流的惯性作用，管道中的水继续流向水库。

在t=2L/a~3L/a时段内，首先是紧近阀

门dxi管段内发生变化，依次传到dx2、dx3・・管段，到3L/a时刻，流速将由-Vo变为

O,压力由Ho变为Ho—AH,管径为D—AD,水的密度变为p-△p。

当阀门全关闭时，水锤波在阀门处的反射特点是“等值同号”反射，即反向波与入射波的数值和符

号不变，从水库传来降压波仍反射为降压波。

（4）t=3L/a~4L/a。

当t=3L/a时，水锤波又回到水库处D点，由于管道压力比水库低

△H,则D点压力不能维持平衡，因此水库的水又向阀方向流动，这时水库将阀门反射

回来的降压波又反射为升压波，到t=4L/a时，管道流速将由O变为Vo,压力由Ho—A

Ho变为Ho，管径、水密度都恢复到初始状态。

T=4L/a称为水锤波的“周期”。

每经一个周期，水锤现象就重复一次上述过程。

水锤波在管中传播一个来回的时间tr=2L/a，称之为“相”，两个相为一个周期T=2tr。

阀门突然开启时，水锤现象与上述情况相反。

如果不存在水力摩阻，则上述的水锤过程将无休止地反复下去，但由于水力摩阻的存在，水锤过程不可能无休止地振荡下去，压力波因摩擦损失而逐渐衰减，在一定时段内逐渐消失。

综上所述，我们可以初步得出以下几点结论：

（1）水锤压力实际上是由于水流速度变化而产生的惯性力。

当突然启闭阀门时，由

于启闭时间短、流量变化快，因而水锤压力往往较大，而且整个变化过程是较快的。

（2）由于管壁具有弹性和水体的压缩性，水锤压力将以弹性波的形式沿管道传播。

（3）水锤波同其它弹性波一样，在波的传播过程中，外部条件发生变化处（即边界

处）均要发生波的反射。

其反射特性（指反射波的数值及方向）决定于边界处的物理特性。

、水锤波的传播速度

在水锤过程的分析与计算中，波速是一个重要的参数。

它的大小与管壁材料、厚度、管径、管道的支承方式以及水的弹性模量等有关。

由水流的连续方程并考虑水体和管壁的弹性后，可导出水锤波的传播速度为

a喜

1DK

\E

1435m/s。

1OOOm/s，埋藏式钢管

900m/s~1200m/s。

在缺乏资料的情况下，露天钢管的水锤波速可近似地取为可近似取为1200m/s，钢筋混凝土管可取

第三节水锤基本方程及边界条件

为求解水锤压力升高问题，需要建立基本方程。

基本方程与相应的边界条件联立,

用解析方法或数值计算方法求解水锤值及其变化过程。

、水锤基本方程

（一卜基本方程

对有压管道而言，不论在何种情况下都应满足水流的运动方程及连续方程。

当水管材料、厚度及直径沿管长不变时，其运动方程为：

HVVf

gVVV0（9-2）

xtx2D

将管道材料及水体当作弹性体考虑，其连续方程为：

（9-3）

卫？

■丄V上0tgxx

式中H――压力水头；

V――管道中的流速，向下游为正；

a——水锤波传播速度；

f——水流摩擦阻力系数；

D――管道直径；

x――距离，其正方向与流速取为一致;

t――时间。

图9-2水击计算示意图

式中Ho和Vo为初始水头和流速；为两个波函数，其量纲与水头H相同，故可视为压力波。

F（t-x/a）表示以波速a沿x轴负方向传播的压力波，即逆水流方向移动的压

力波，称为逆流波；f（t+x/a）表示以波速a沿x轴正方向传播的压力波，即顺水流方向移动的压力波，称为顺流波。

任何断面任何时刻的水锤压力值等于

两个方向相反的压力波之和，而流速值为两个压力波之差再乘以-g/a。

如果知道了t时刻在x位置处的水锤波函数F（t-x/a），则当时间变为ti=t+At，研究

xi=x+aAt处的逆流波函数

二、水锤的边界条件

应用水锤基本方程计算压力管道中水锤时，首先要确定其起始条件和边界条件。

（1）、初始条件

当管道中水流由恒定流变为非恒定流时，把恒定流的终了时刻看作为非恒定流的

开始时刻。

即当t=0时，管道中任何断面的流速V=V°；如不计水头损失，水头H=H0。

（2）、边界条件

1.管道进口。

管道进口处一般指水库或压力前池。

水库水位变化比较慢，在水锤计算中不计风浪的影响，认为水库水位为不变的常数是足够精确的。

压力前池的水位变化情况与渠道的调节类型有关。

自动调节渠道的前池水位变化虽大，但与管道中水锤计算时间相比，变化还是缓慢的。

非自动调节渠道，水位变化较小，一般只有几米，在水锤计算中也认为前池水位不变。

所以管道进口边界条件为：

Hp=Ho

2.分岔管。

分岔管的水头应该相同，即

Hpi=Hp2=Hp3=・・・=Hp

分岔处的流量应符合连续条件，即

工Q=0

3•分岔管的封闭端。

在不稳定流的过程中，当某一机组的导叶全部关闭，或某一机组尚未装机，而岔管端部用闷头封死，其边界条件为：

Qp=O

4.

:

调压室内有自由水

调压室。

把调压室作为断面较大的分岔管，其边界条件为面，而隧洞、调压室与压力管道的交点和分岔管相同。

5、水轮机。

水电站压力管道出口边界为水轮机，水轮机分冲击式和反击式，两种型式的水轮机对水锤的影响不同。

（1）冲击式水轮机。

冲击式水轮机的喷嘴是一个带针阀的孔口。

水轮机转速变化对

孔口出流没有影响，对冲击式水轮机，喷嘴全开时断面积为3max，流量系数为00,根

据《水力学》的孔口出流规律，过流量为：

Qmax0max■.2gH。

当孔口关至

i2g（HoH）

（2）反击式水轮机。

反击式水轮机有如下特点：

（i）反击式水轮机有蜗壳、尾水管及导水叶，过流特性与孔口出流不完全相同。

（ii）反击式水轮机的转速与水轮机的流量互相影响。

（iii）流量突然改变时，不仅在压力管道中，而且在蜗壳、尾水管中也发生水锤。

尾水管中发生的水锤现象与蜗壳相反，即导水叶关闭时发生负水锤，开启时发生正水锤。

蜗壳、尾水管中的水锤影响水轮机的流量，继而又对水锤产生影响。

由此可见，反击式水轮机的过水能力与水头H、导叶开度a和转速n有关。

即

Q=Q（H,a,n），需要综合运用管道水锤计算方程、水轮机运转特性曲线、水轮机组转速方程等进行求解，因此增加了问题的复杂性。

为了简化计算，常假定压力管道出口边界条件为冲击式水轮机，然后再加以修正。

（3）

图9-3水轮机开度变化规律

开度按直线规律变化

水轮机导叶和阀门的关闭规律与调速系统的特性有关，实际的关闭规律如图9-3所示。

从全开（T0=到全关（T=0）的全部历时为Tz,曲线开始一段接近水平，关闭的速度极慢，这是由于调节机构的惯性所决定的，在这段过程中，弓I起的水锤压力很小，对水锤计算没有多大实际意义。

在接近关闭终了时，阀门的关闭速度又逐渐减慢，曲线向后延伸，这种现象只对阀门关闭接近终了时的水锤压力有影响。

因此为了简化计算，常取阀门关闭过程的直线段加以适当延长，得到Ts,Ts称为有效关闭时间。

在缺乏资

料的情况下，可近似取Ts=~Tz。

2L

2Lt

直线规律关闭时，一个相长tr，一个相的开度变化—

aTs

负号表示阀门关闭；正号表示阀门开启。

第四节简单管水锤的解析计算

简单管是指压力管道的管径、管壁材料和厚度沿管长不变。

解析法的要点是采用数学解析的方法，引入一些符合实际的假定，直接建立最大

水锤压力的计算公式。

简单易行，物理概念清楚，可直接得出结果。

、直接水锤和间接水锤

水锤有两种类型：

直接水锤和间接水锤。

（1）直接水锤

水锤波在管道中传播一个来回的时间为2L/a,称为“相”。

当水轮机开度的调节时

间Ts<2L/a时，由水库处异号反射回来的水锤波尚未到达阀门之前，阀门开度变化已经终止，水管末端的水锤压力只受开度变化直接引起的水锤波的影响，这种水锤称为直接水锤。

由于水管末端未受水库反射波的影响，因此基本方程（9-6）和（9-7）中的波函数

f（t+x/c）=O,然后从二式中消去F（t-x/c）得直接水锤公式

HHH。

a（VVo）（9-13）

g

公式（9-13）只适用于Ts<2L/a的情况，由此式可得出如下结论：

（1）当阀门关闭时，管内流速减小，V-Vo<0为负值，△H为正，产生正水锤；反之

当开启阀门时，即V-Vo>O,AH为负，产生负水锤。

（2）直接水锤压力值的大小只与流速变化（V-V0）的绝对值和水管的水锤波速a有

关，而与开度变化的速度、变化规律和水管长度无关。

当管道中起始流速Vo=4m/s,a=1000m/s，终了流速V=0时，压力升高值为：

H（VV0）a/g=-1000（0-4）/=407.7m，因此在水电站中应当避免直接水锤。

（2）间接水锤

若水轮机开度的调节时间Ts>2L/a,当阀门关闭过程结束前，水库异号反射回来

的降压波已经到达阀门处，因此水管末端的水锤压力是由向上游传播的水锤波F和反

射回来的水锤波f叠加的结果，这种水锤称为间接水锤。

降压波对阀门处产生的升压波起着抵消作用，使此处的水锤值小于直接水锤值。

发生间接水锤时，水锤压力波的消减、增加过程是十分复杂的。

间接水锤是水电站中经常发生的水锤现象，也是要研究的主要对象。

工程中最关心的是最大水锤压力。

由于水锤压力产生于阀门处，从上游反射回来的降压波也是最后才达到阀门，因此最大水锤压力总是发生在紧邻阀门的断面上。

下面应用前面的水锤连锁方程（9-10）和（9-11）及管道边界条件，推求阀门处各相水锤压力

的计算公式。

二、计算水管末端各相水锤压力的公式

（一）第一相末的水锤压力

（1）设阀门为A点，水库为

A

0

t=0时，

边界和初始条件:

所以:

B点，水锤波从A到B点的连锁方程为:

tB2（V（AVtB）

0A0;在水库进水口

AB、

（V。

Vt）

A

V0

B点，tB

将A点边界条件V（A

0.10A代入上式:

A

V0

（2）水锤波从B到A的连锁方程:

B

t

A

2t

因tB0,VtB

o和A点边界条件

A

V2t

A

2t

0

B

Vt

B

Vt0

2（VtB

A、

V2t）

2t,12A，上式变为

02t■1

2A）

因2t=tr=2L/a为一个相长，用1表示第一相末，

得到:

（二）第二相末的水锤压力

（9-14）

（l）写出水锤波从AtB的连锁方程式:

AB

2t3t

B

V3t

由b点的边界条件得

=0，上式可改写成:

2A/2

所以:

B

V3t

2t1

A

2t

2t一1

A

2t

2

A

2t

B

V3t

（2）写出水锤波从BtA的连锁方程式：

BA

3t4t

A

2t如A

和V4t

B

V3t

A

V4t

把30,VBt

A

2t

4t.14A代人上式,

并用2代替4t表

示第二相末，得:

（三）第n相末的水锤压力

用同样原理可以得出以后任意n相末的水锤压力计算公式，其一般公式为:

（9-16）

利用式（9-14）~（9-16），可以依次求出各相末阀门处的水锤压力，得出水锤压力随时间的变化关系。

上面是阀门关闭情况，当阀门或导叶开启时，管道中的流速增加，压力降低，产生负水锤，其相对值用y表示，用同样的方法可求出各相末计算公式。

此时VtAtA,1%,求出的y本身为负值。

1...1y10

2

（9-17）

1i

n1

yn

（9-18）

nJyn0

i

yi

i1

2

上述水锤压力计算公式的条件：

（1）

没有考虑管道摩阻的影响，

因此只适用于不计

摩阻（如水头较高、管道较短等）的情况；

（2）采用了孔口出流的过流特性，只适用于冲击式水轮机，对反击式水轮机必须另作修改；（3）这些公式在任意开关规律下都是正确的，可以用来分析非直线开关规律对水锤压力的影响。

、水锤波在水管特性变化处的反射

水锤发生后，水锤波在水管末端和水管特性变化处（水管进口、分岔、变径段、

阀门等）都要发生反射。

当入射波到达水管特性变化处之后，一部分以反射波的形式折回，一部分以透射波的形式继续向前传播。

反射波与入射波的比值称反射系数，以r表示。

透射波与入射波的比值称透射系数，以s表示，两者的关系为

sr1（9-19）

（一）水锤波在水管末端的反射

水锤波在水管末端的反射特性取决于水管末端的出流规律。

对于水斗式水轮机,

其喷嘴的出流规律为v1—，当0.5时，可近似地取为v（1/2）。

在

入射波未达到的时刻，00,v0。

（9-20）

根据水锤常数和任意时刻的开度，可利用式（9-20）确定阀门在任意时刻的反射系

数。

当阀门完全关闭时，0,r1，阀门处发生同号等值反射。

（二）水锤波在管径变化处的反射

如图9-4所示的变径管，入射波F1从1管传来，在变径处发生反射。

反射波为f1,

透射波为F2,由方程（9-6）和（9-7）及水流在变径处的连续性，可推导出反射系数

（9-21）

（三）水锤波在分岔处的反射

如图9-5所示，入射波F1从1管传来，在分岔处发生反射，反射波为f1,透射波

为F2和F3,根据基本方程（9-6）和（9-7）及此处水流的连续性，导出反射系数为

（9-22）

122331

式中iaQ,Q为总管流量，A为水

2gH°Ai

管断面积。

四、开度依直线变化的水锤

进行水锤计算，最重要的是求出其最大

值。

在开度依直线规律变化情况下，不必用连锁方程求出各相末水锤，再从中找出最大值，可用简化方法直接求出。

（一）开度依直线变化的水锤类型

当阀门开度依直线规律变化时，根据最大压力出现的时间可归纳为两种类型：

Aa

（1）最大水锤压力出现在第一相末，max1，如图9-6（a），称为第一相水锤。

（2）最大水锤压力出现在第一相以后的某一相，其特点是最大水锤压力接近极限值

m，即m>1，如图9-6（b），称为极限水锤。

产生这两种水锤现象的原因是由于阀门的反射特性不同造成的，阀门处的反射特性可由其反射系数确定。

1•第一相水锤

根据式（9-20），当。

<1时，r为正，水锤波在阀门处的反射为同号。

在阀门关闭过程中，阀门处任意时刻的水锤压力由三部分组成：

阀门不断关闭所产生的升压波、经水库反射回来压力波、经阀门反射向上游的压力波。

（1）第一相中，根据水库异号反射的特性，升压波到达水库后反射回的降压波还未到达阀门处，因此该处水锤压力即是阀门关闭所产生的升压波，在第一相末达到：

。

（2）第二相末，水库传来的降压波到达阀门处，如果此时阀门处具有同号反射的特性，则在该处反射仍为降压波，两个降压波之和将超过第二相中由于阀门关闭所产生的升压波，因而第二相末的水锤压力AA

丿力2<1。

（3）

图9-6开度为直线关闭时的水击类型

第三相末，由于第二相中阀门同

号反射回去的降压波，经水库异号反射为升压波，这两个升压波共同作用，又使阀门

处的水锤压力开始升高，

AA

3>2。

根据阀门同号反射的规律，水锤压力将环绕某一mA值上下波动，最后趋于

由于最大水锤压力出现在第一相末，1a>mA,故称为第一相水锤。

2.极限水锤

根据式（9-20），当。

>1时，r为负，水锤波在阀门处的反射为异号。

在阀门关

闭过程中，阀门处任意时刻的水锤压力仍由上述三部分组成。

第一相末，水库反射回的降压波还未到达阀门处，该处水锤压力只是阀门关闭所产生的升压波，即。

第二相末，水库传来的降压波到达阀门处，因阀门处为异号反射，则在该处反射为升压波，它和在第二相中阀门继续关闭产生的升压波共同作用，使第二相中阀门处的水锤压力继续升高，使2A>iA。

在以后各相，阀门处水锤压力逐渐增加，趋近某一

极限值m。

由于最大水锤压力为mm，故称为极限水锤。

（二）开度依直线变化时水锤的简化计算

当调节阀门按直线规律启闭，

t与

0的关系为：

当阀门关闭时

t0

tr

2L

0

（9-23）

Ts

aTs

当阀门开启时

t0

tr

2L

0

（9-24）

Ts

aTs

2L/aT；关闭情况

2L/aTs开启情况

1.第一相水锤计算的简化公式

当彳<时，,111-，则（9-14）可简化为：

2

1（1-）0—

22

^_^Vmax

令=,称为水锤特性常数，关闭时用正值，开启时

gHoTs

为负值。

考虑到1和0的关系，代入上式可解得第一相末水锤压力值为：

（9-25）

关闭阀门时

a2

y11-

I0

0=1，有

故在丢弃负荷的情

发生第一相水锤的条件是0<1,对于丢弃负荷情况，aVmax/2gH01。

若a=1000m/s,Vmax=5m/s，贝VHo>25Om,况下，只有高水头电站才有可能出现第一相水锤。

（2）极限水锤计算简化公式

根据（9-15），第n相和第n+1相末的水锤压力计算公式为：

1

1.1

上二式相减，得:

1..1

代入上式，则

代人上式得

a2

m

2——

0

将上式代入式1mAmA中，解得值为

（9-30）

9-7所示。

4。

（1°）

12°

如果公式（9-30）满足，则第1A。

公式（9-30）代表一根曲线，如图

图中同时绘出了°的直线。

曲线表示极限水锤和第一相水锤的分界线，直

线0表示第一相水锤和直接水锤的分界线。

共有五个分区：

I区为极限正水锤；II区为第一相正水锤；III区为直接水锤；

区为极限负水锤；V区为第一相负水锤。

简单判别方法：

0＜时，常发生第一相水锤；

0＞时，常发生极限水锤；

＜0＜时，则随值的不同而发生第一相或极限水锤，个别情况下发生直接水锤。

此时按图9-7判别。

最后，为了方便水锤压力的计算，将计算公式汇总于表10-1。

五、起始开度对水锤的影响

水电站可能在各种不同的负荷情况下运行，当机组满负荷运行时，起始开度0=1；

当机组只担任部分负荷运行时，0v|。

因此机组由于事故丢弃负荷时的起始开度0可

能有各种数值。

从前面的水锤压力计算公式可以绘制出图9-8。

图中的曲线和分界点说

明了起始开度对水锤压力的影响。

表10-1水锤压力计算公式汇总表

O

2

2

ym

y1

计

直接水击

直接水击

间接水击

间接水击

了终

水平线。

对第一相水锤

，随着0的减小而增大,

所以在图中表示为一根

曲线。

对直接水锤,

0，为一通过坐标轴原点的直线，其斜率为

2。

图中三

条曲线的交点为:

（1）直接水锤和第一相水锤:

令dA2o和/

相等，可以解出:

（2）第一相水锤和末相水锤

令A2

11

0

相等，可以解出:

1/

因此可得出以下结论:

（