第四节简单式和阻抗式调压室的水位波动计算.docVIP

第四节 简单式和阻抗式调压室的水位波动计算 ? ? 调压室水位波动计算常用的方法有解析法和逐步积分法。解析法较简便，可直接求出最高和最低水位，有时精度差，不能求出波动的全过程，常用以初步决定调压室的尺寸。逐步积分法是通过逐步计算以求出最高和最低水位，其最大优点是可以求出波动的全过程和求解复杂的问题。逐步积分法可分为图解法和列表法(数学积分法)，两者原理相同。图解法简便，醒目，列表法较精确。逐步积分法一般用于后期的设计阶段。近年来随着电子计算机的发展，在工程设计中已越来越多地采用电算法，以同时解决调压井涌波、水锤压力及机组速率上升的复杂计算，特别是研究各参数的影响时，电算法更为优越。下面我们主要介绍解析法和图解法。电算法则在第九节中介绍。? ? 一、水位波动计算的解析法? ? （一）丢弃全负荷情况? ? 当丢弃全负荷后，水轮机的流量Q=0,连续性方程式(15-3）变为? ? 在水流进出调压室时，如考虑由于转弯、收缩和扩散引起的阻抗孔口水头损失K,则动力方程式(15-4)变为式中，其中为水头损失系数（为一常数）。,其中和分别为流量流过引水道和进出调压室所引起的水头损失。? ? 令，则，，将以上关系代入式(15-7)，两边除以 ,并令，则得? ? 将 V=y代入式（15-6），并和式（15-8）消去dt，得? ? 再令，即Z=SX，dZ=SdX，代入上式，得? ? 系数S具有长度因次，用以表示“引水道—调压室”系统的特性。X和Xo均为无因次的比值。? ? 式（15-10）为变数X和的一阶线性微分方程式，积分后得? ? 积分常数C可由起始条件决定。波动开始时，t=0， V=，即y=1,Z=,X=Xo,以? ? y=1,X=Xo代入上式，得? ? 故（15-10）的最后解答为? ? 对于调压室的任何水位（用X表示），可用上式算出与之对应的引水道的流速V = y，也可以进行相反的计算，但不能求出流速V与水位X对于时间t的关系，因此，不能求出水位波动过程。? ? 1.最高水位的计算? ? 欲求波动的最高水位，只需求出即可。在水位达到最高时，V=0，即y=0，代人式(15-11)得? ? 两边取对数得式中的符号在静水位以上为负，在静水位以下为正。? ? 式（15-12）适用于阻抗式调压室。对于简单式调压室，附加阻抗可以忽略不计，即η=0，则式（15-12）变为? ? 如流量不是减小至零，则不能应用上列公式，只好应用图解法或数学积分法求解。值亦可由图15-5中曲线A，根据查出，算出。 图 15-5 简单式调压室丢弃负荷最大振幅计算曲线 ? ? 有压引水道的水头损失对调压室的最大涌波值影响较大。例如对设有简单调压室的水电站，当机组丢弃全部负荷时。考虑与不考虑引水道水头损失的调压室水位最大升高值为与 ，其比值与引水道长度的关系曲线见图15-6。随引水道长度的增大，比值开始迅速降低，随后则渐趋平缓：例如引水道长lkm时，比减小13%，当引水道长20km时，则减小25％。 图 15-6 与引水道长度关系曲线 ? ? 2.波动第二振幅的计算? ? 全弃负荷后，调压室水位先升至最高水位 ，然后又下降至最低水位，称为第二振幅。进行第二振幅的计算，是为了保证调压室水位不致降入隧洞中，以免带人空气。这时水由调压室流出，故和K的符号与前相反，用相同的方法加以处理后得? ? 如η=0,则? ? 求出，随即可求出第二振幅的。? ? 在应用式(15-14)和式(15-15)时，要特别注意和的符号，前者为负，后者为正。值也可从图15-5中曲线A、B求得。? ? （二）增加负荷情况? ? 当突然增加负荷时，波动微分方程式不能像丢弃全负荷那样进行积分，只能作某些假定求出近似解。? ? 当水电站的流量由增至时，若阻抗η=0，规范建议按下面近似公式求解最低涌波水位式中，为无因次系数，表示“引水—调压室”系统的特性，与前面的S相比，，式(15-16)可绘出曲。见图15-7。 图 15-7 简单调压室增加负荷最低振幅计算曲线 ? ? 二、水位波动计算的图解法? ? （一）图解法的基本原理? ? 图解法是用作图的办法来积分波动为分方程式，其实质是用有限差的比值和代替基本方程中的导数，然后逐时段地求解最后得出水位、流速和时间的关系。连续方程和动力方程可写为? ? 并进一步改写为 式中，当计算时段t选定后，A、α、β均为常数。? ? 计算的基本假定是：在时段t的过程中调压室的水位Z和引水道的流速V保持不变，而集中变化在时段的末尾，即假定Z和V的变化是“阶梯式”的，这是图解法近似之处，一般说来，t取得愈小，结果愈精确，但t太小，不仅加大计算工作量，而且会因图解不便而使误差增大，