调压室水力计算分解
调压室水力计算
调压室的水力计算1. 调压室断面计算当上游死水位，下游为最低水位，最小水位H min=188.9m，三台机满发，引水道糙率取最小值，压力管道糙率取最大值，通过水轮机的流量为57m3s⁄，则此时的引水隧洞水头损失的计算如表格1，压力钢管水头损失的计算如表格2。
引水道应选可能的最小糙率0.012，压力管道应选择可能的最大糙率0.013。
表格1引水隧洞水头损失表表格2压力钢管水头损失表F Tℎ>Lf2αgH1=Lf2αg(H0−ℎw0−3ℎwT)=45.548m2其中H0——最小水头损失，H0=188.9m；ℎw0——引水隧洞损失，ℎw0=17.802+0.296=18.098；ℎwT0——压力管道水头损失，ℎwT=3.110+2.805=5.915m；L——引水隧洞长度，12662m；g——重力加速度，g=9.81m/s2 f——引水隧洞面积，16.62m2。
α——引水道阻力系数v0=Qf=5716.619=3.43m s⁄α=ℎw0v02=18.0983.4302=1.5385为了保证大波动的稳定，一般要求调压室断面大于托马斯断面，初步分析时可取（1.0~1.1）F Tℎ，作为调压室的设计断面。
这里选取D=7.8m，则系数k为：F k=47.784k=F k/F Tℎ=1.052. 最高涌波水位计算按正常蓄水位时共用同一调压室的三台机组全部满载运行瞬时丢弃全部负荷（即流量由Q max=57减至流量Q=0）作为设计工况。
引水隧洞的糙率取尽可能的最小值（能耗少，涌波高）。
n=0.012引水道损失由表格1和表格2得：ℎw0=ℎw0程+ℎw0局=17.802+0.296=18.098mv0为时段开始时管中流速v0=Qf=3.43m s⁄；f为引水隧洞断面面积。
F为调压井断面面积，145.267m2；引水隧洞长L=12662m，g=9.81m s2⁄得引水道—调压室系统的特性系数。
λ=Lfv022gFℎw0=12662×16.62×3.4322×9.81×47.784×18.098=145.89令X0=ℎw0λ=0.124，X=zλ，则要求最高涌波水位z max，只需要求出X max=z maxλ即可。
水厂压力调控算法
水厂压力调控算法
水厂压力调控算法是用于调节和优化水厂供水系统压力的一种算法。
其目的是确保供水系统的稳定运行，满足用户的水量和水压需求，同时避免能源的浪费。
水厂压力调控算法通常基于以下几个关键参数：
1.供水压力：通过传感器实时监测供水系统的压力，确保其在设定的合理范围内。
如果压力过高或
过低，算法会进行相应的调整。
2.流量：监测供水系统的流量，以了解当前的水量需求。
流量数据可以用于预测未来的需求，并据
此调整供水压力。
3.水位：监测储水设施（如蓄水池、水塔等）的水位，以确保其在水位过低时及时补水，避免供水
中断。
4.水泵状态：监测水泵的工作状态，包括转速、功率等，以确保其正常运行，并根据需要调整水泵
的运行策略。
基于这些参数，水厂压力调控算法可以采用以下策略：
1.PID控制：利用比例-积分-微分（PID）控制器对供水压力进行实时调节。
PID控制器可以根据当
前压力与目标压力的差值，计算出一个控制量，用于调整水泵的转速或开度，从而实现对供水压力的精确控制。
2.模糊控制：利用模糊逻辑理论，根据供水系统的实际运行情况，对供水压力进行模糊推理和决策。
模糊控制可以处理供水系统中存在的不确定性和非线性问题，提高压力调控的鲁棒性。
3.优化算法：采用遗传算法、粒子群优化算法等优化算法，对供水系统的运行参数进行优化，以实
现节能降耗、提高供水质量等目标。
在实际应用中，水厂压力调控算法需要根据具体的水厂供水系统和用户需求进行定制和优化。
同时，还需要考虑算法的稳定性和可靠性，确保其在实际运行中的稳定性和可靠性。
空调系统水力平衡分解
三、变流量系统水力平衡策略
分集水器侧的调节与平衡：
（2）电动调节阀控制方式：
利用压差变送器采集分集水器之间的压差，然后输出4-20mA的标准电流 信号到控制器。与控制器上设定压差相比较后，再输出控制信号到电动调节阀， 改变其开度。通过调节电动调节阀改变旁通水量，从而保证分集水器之间的压 差△P恒定。这时，如果分集水器上某一分支回路流量变化，其它回路不再受 其影响，进而系统实现动态水力平衡。这种方法和上一种方法可达到相同的效 果，而且由于设定压差可以根据外部环境加以改变，实现变压差控制，进一步 降低了制冷主机的能耗。
空调系统水力平衡分析与调节
前言
随着时代的进步和科技的发展，在经济社会日益繁荣的今天，人 们对生活的舒适程度要求不断提高。中央空调在我国许多大型商场等 建筑中已经成为了标准化的配置。中央空调系统要保持其稳定高效工 作，一个重要条件就是要保证其水力的平衡调节。在能源情况并不乐 观的今天，保证暖通空调系统使用的节能是暖通专业关注的重要问题。
二、水力平衡的调试
静态平衡阀： 动态流量平衡阀： 动态压差平衡阀：
三、变流量系统水力平衡策略
变流量系统：
两个主要任务：调节和平衡。 目的 调节：使得各分支环路的流量实时地变化，力求准确地与冷热负荷的要求匹 配，达到良好的空调及供暖效果； 平衡：克服凋节过程中出现的水力失调，两者必须有机地结合起来。 对系统判断的标准有以下几点： 一、流量控制精度； 二、灵敏度； 三、稳定性， 四、能够有效地降低能耗。
三、变流量系统水力平衡策略
空气处理机组支路的调节与平衡：
三、变流量系统水力平衡策略
空气处理机组支路的调节与平衡：
方法一、采用压差控制阀+电动调节阀。 在空气处理机组分支环路回水端安装压差控制阀，保持图中A、B两点问的压差 不变，使电动调节阀的流量只受末端负荷变化的影响。通过改变电动调节阀的开度 来调节流量，使之满足负荷变化的要求，如果电动调节阀开度不变，流经此支路上 空气处理机组的流量也基本保持不变。 方法二、直接在空气处理机组支路上安装电动调节动态流量平衡阀。 由于它的流量只是开度的单值函数，只需根据负荷需求调节阀门的开度便可达 到所需流量，一经开度设定，即使其他支路发生变化，其流量维持不变。每个电动 调节动态流量平衡阀产品其流量—开度特性曲线在出厂前已精确确定，流量只和自 身有关，而不受末端设备和其他管路的影响，所以很容易做到流量的精确控制。而 对于第一种方案，尽管压差保持恒定时其流量也只和调节阀的开度有关，但流量一 开度的关系具有不确定性，因为即使同样的调节阀，和不同的末端设备和管道组成 的支路具有不同的阻力特性，其流量不仅和阀门的开度有关，还和阀门的阀权度有 关，因而其流量控制精度要低于采用电动调节动态流量平衡阀的方案。
第二节 调压室的工作原理和基本方程
第二节调压室的工作原理和基本方程一、调压室的工作原理水电站在运行时负荷会经常发生变化。
负荷变化时，机组就需要相应地改变引用流量，从而在引水系统中引起非恒定流现象。
压力管道中的非恒定流现象（即水锤现象）在上一章中已经加以讨论。
引用流量的变化，在“引水道-调压室”系统中亦将引起非恒定流现象，这正是本节要加以讨论的。
图13-5为一具有调压室的引水系统。
当水电站以某一固定出力运行时，水轮机引用的流量亦保持不变，因此通过整个引水系统的流量均为，调压室的稳定水位比上游水位低，为通过引水道时所造成的水头损失。
当电站丢弃全负荷时，水轮机的流量由变为零，压力管道中发生水锤现象，压力管道的水流经过一个短暂的时间后就停止流动。
此时，引水道中的水流由于惯性作用仍继续，流向调压室，引起调压室水位升高，使引水道始末两端的水位差随之减小，因而其中的流速也逐渐减慢。
当调压室的水位达到水库水位时，引水道始末两端的水位差等于零，但其中水流由于惯性作用仍继续流向调压室，使调压室水位继续升高直至引水道中的流速等于零为止，此时调压室水位达到最高点。
因为这时调压室的水位高于水库水位，在引水道的始末又形成了新的水位差，所以水又向水库流去，即形成了相反方向的流动，调压室中水位开始下降。
当调压室中水位达到库水位时，引水道始末两端的压力差又等于零，但这时流速不等于零，由于惯性作用，水位继续下降，直至引水道流速减到零为止，此时调压室水位降低到最低点。
此后引水道中的水流又开始流向调压室，调压室水位又开始回升。
这样，引水道和调压室中的水体往复波动。
由于摩阻的存在，运动水体的能量被逐渐消耗，因此，波动逐渐衰减，最后全部能量被消耗掉，调压室水位稳定在水库水位。
调压室水位波动过程见图13-5中右上方的一条水位变化过程线。
当水电站增加负荷时，水轮机引用流量加大，引水道中的水流由于惯性作用，尚不能立即满足负荷变化的需要，调压室需首先放出一部分水量，从而引起调压室水位下降，这样室库间形成新的水位差，使引水道的水流加速流向调压室。
空调系统水力平衡分解
三、变流量系统水力平衡策略
分集水器侧的调节与平衡：
(1) 压差控制阀控制方式：
压差控制阀控制方式是在分集水器旁通管上设压差控制阀，以此来控制系 统中分集水器之间的压差。当系统中某一支环路流量变化时，由于压差控制阀 的调节作用，使分集水器之间的压差△P保持不变，其余支环路的流量并不随 之发生变化，从而使系统实现动态水力平衡。在这种情况下，水泵的工作点维 持不变，扬程和流量均保持恒定，多余的流量通过压差控制阀旁路流回集水器。 系统在部分负荷工况下运行时，流向末端装置的总流量变小，因水泵流量不变， 旁通的水流量增多，由于这部分冷冻水没有经过换热，使得主机的回水温度降 低，主机部分负荷运行，降低能耗。
二、水力平衡的调试
静态平衡阀： 动态流量平衡阀： 动态压差平衡阀：
三、变流量系统水力平衡策略
变流量系统：
两个主要任务：调节和平衡。 目的 调节：使得各分支环路的流量实时地变化，力求准确地与冷热负荷的要求匹 配，达到良好的空调及供暖效果； 平衡：克服凋节过程中出现的水力失调，两者必须有机地结合起来。 对系统判断的标准有以下几点： 一、流量控制精度； 二、灵敏度； 三、稳定性， 四、能够有效地降低能耗。
其中水力失调的解决就是保证节能的重要措施之一。
CONTENTS
01
水力失调
02
水力平衡的调试
变流量系统水力平衡策略 结语
03
04
一、水力失调
定义： 在空调水系统中，各空调设备的实际流量与设计流量
不一致性称为该设备的水力失调。
水力失衡程度： 分类：
λ =q实际/Q设计
静态水力失调和动态水力失调
一、水力失调
动态水力失调
定义：
动态水力失调是指在中央空调系统运行过程中，由于终端空调设备 数量多，当终端空调设备开关或阀门开度变化时，管路流量、压力产生 波动和变化，引起其他管路空调设备流量、压差波动，偏离设计要求而 产生的水力失调。
水击与调压室计算与演示
水击与调压室计算与演示1、调压室简介为了减小水锤压力，常在有压引水隧洞（或水管）与压力管道衔接处建造调压室，如图1所示。
调压室利用扩大的断面和自由水面反射水锤波，将有压引水系统分成两段：上游段为有压引水隧洞，调压室使隧洞基本避免了水锤压力的影响；下游段为压力管道，由于长度缩短可，从而降低了压力管道中的水锤值，改善了机组的运行条件。
图1 水电站调压室调压室的功用有以下几点：（1）、反射水锤波。
基本上避免（或减小）压力管带中的水锤波进入有压引水道。
（2）、缩短压力管道的长度。
从而减小压力管道及厂房过流部分中的水锤压力。
（3）、改善机组在负荷变化时的运行条件及系统供电质量。
调压室的工作原理是：增大的水面反射水锤波，引水道中水体动能和势能相互转换。
2、水击与调压室计算程序演示2.1 甩负荷工况当水电站丢弃全部负荷时，水轮机的流量由Q0变为零，压力管道中发生水锤现象。
此时，上游调压室水位先上升，下游调压室水位先下降。
引（尾）水道中水流在惯性作用下继续流动，从而引起调压室水位上升（下降），当水位达到极值后，由于调压室和水库的水位差作用，水流开始倒流。
如此往复流动，实现动能和势能的转换，并在阻力消耗下衰减。
用程序演示：起始开度设置为0.8，终了开度为0。
计算过程如下（原始数据见图2，计算结果见图3。
）图2 甩负荷工况参数输入图3 甩负荷工况计算结果和过程线2.2 增负荷工况当水电站增加负荷时，水轮机引用流量加大，压力管道中也出现水锤现象。
此时，上游调压室水位先下降，下游调压室水位先上升。
机组首先开始增大引用流量，水流流出（入）上游（下游）调压室，调压室水位变化，当调压室和水库的水位差达到极值后，水位差的作用使水流开始倒流。
如此往复流动，实现动能和势能的转换，并在阻力作用下快速衰减。
用程序演示：设置起始开度为0.9，终了开度为0。
计算过程如下（原始数据见图4，计算结果见图5。
）图4 增负荷工况参数输入。
有压引水系统水力计算
一、设计课题水电站有压引水系统水力计算。
二、设计资料及要求1、设计资料见《课程设计指导书、任务书》；2、设计要求： （1）、对整个引水系统进行水头损失计算； （2）、进行调压井水力计算球稳定断面； （3）、确定调压井波动振幅，包括最高涌波水位和最低涌波水位； （4）、进行机组调节保证计算，检验正常工作状况下税基压力、转速相对值。
三、调压井水力计算求稳定断面<一>引水道的等效断面积：∑=ii fL Lf ， 引水道有效断面积f 的求解表所以引水道的等效断面积∑=ii fL Lf =511.28/21.475=23.81 m 2 <二>引水道和压力管道的水头损失计算： 引水道的水头损失包括局部水头损失h 局和沿程水头损失h 沿两部分 压力管道的水头损失包括局部水头损失h 局和沿程水头损失h 沿两部分1,22g 2h Qϖξ局局=g ：重力加速度9.81m/s 2 Q ：通过水轮机的流量取102m 3/s ω ：断面面积 m 2ξ：局部水头损失系数局部水头损失h 局计算表从上表中可以看出：引水道的h 局=0..037+0.204+2.202=2.713m 压力管道的h 局=4.464m2,23422n h QRlϖ=沿n ：糙率系数，引水道糙率取最小值0.012；压力管道取最大值0.013 l ：引水道长度 m ω ：断面面积 m 2R ：为水力半径 m Q ：通过水轮机的流量m 3/s沿程水头损失h 程计算表取最大值0.013。
上表中栏号1、2、3、4、5、6中的Q=1023m /s ；栏号7中Q=96.93m /s ；栏号 8中Q=64.63m /s ；栏号9中Q=32.33m /s ； 栏号10中Q=32.33m /s ；栏号11中Q=32.33m /s ；1=h +h f h 局程=（0.307+0.203+2.202）+（0.007+0.011+0.018+0.815+0.011+0.032） =3.606m压力管道沿程水头损失：w h =0.109+0.040+0.004+0.003+0.057=0.213m<三>、调压井稳定断面的计算为使求得的稳定断面满足各种运行工况的要求，上游取死水位，下游取正常尾水位情况计算00013wT w h h H H --=h w0：引水道水头损失，大小为3.606mh wT0：压力管道沿程水头损失，大小为0.213mH 0：静水头，H 0=上游死水位—下游正常尾水位=1082.0－1028.5=53.5m则1H =0H 13f w h h --=53.5-3.606-3⨯0.213=49.255m取K=1.5，D=5.5m ，s m A Q v /284.481.23102===， α=h w0/v 2=0.196当三台机组满出力时，保证波动稳定所需的最小断面：F =k12LfgaH其中k=1.2: 1.5，g=9.812m /s , L=511.28m , f=23.81（2m ）, 1H = 47.255(m)F =k 12LfgaH F=255.47196.081.9281.2328.5115.1⨯⨯⨯⨯⨯=100.492m则调压室断面直径D=πF4=14.349.1004⨯=11.311m.四．水位波动的计算：h w0为引水道的水头损失，包括沿程损失和局部水头损失两部分，沿程水头损失 h 程计算表<一>、最高涌波水位：（1）. 当上游为校核洪水位1097.35m ，下游为相应的尾水位1041.32m ，电站丢弃两台机时，若丢荷幅度为30000——0KW ，则流量为63.6——0m 3/s ，用数解法计算。
第十五章 调压室3
第十五章调压室第三节调压室的基本类型一、调压室的基本布置方式根据水电站不同的条件和要求，调压室可以布置在厂房的上游或下游，在有些情况下，在厂房的上下游都需要设置调压室而成双调压室系统。
调压室在引水系统中的布置有以下四种基本方式。
图15-2 调压室的几种布置方式1、上游调压室（引水调压室）调压室在厂房上游的有压引水道上，如图15-1所示，它适用于厂房上游有压引水道比较长的情况下，这种布置方式应用最广泛，后面我们还要较详细地讨论。
2、下游调压势（尾水调压室）当厂房下游具有较长的有压尾水隧洞时，需要设置下游调压室以减小水锤压力，如图15-2(a)所示，特别是防止丢弃负荷时产生过大的负水锤，因此尾水调压室应尽可能地靠近水轮机。
尾水调压室是随着地下水电站的发展而发展起来的，均在岩石中开挖而成，其结构型式，除了满足运行要求外，常决定于施工条件。
尾水调压室的水位变化过程，正好与引水调压室相反。
当丢弃负荷时，水轮机流量减小，调压室需要向尾水隧洞补充水量，因此水位首先下降，达到最低点后再开始回升；在增加负荷时，尾水调压室水位首先开始上升，达最高点后再开始下降。
在电站正常运行时，调压室的稳定水位高于下游水位，其差值等于尾水隧洞中的水头损失。
尾水调压室的水力计算基本原理及公式与上游调压室相同，应用时要注意符号的方向。
3、上下游双调压室系统在有些地下式水电站中，厂房的上下游都有比较长的有压输水道，为了减小水锤压力，改善电站的运行条件，在厂房的上下游均设置调压室而成双调压室系统，如图15-2 (b)所示。
当负荷变化水轮机的流量随之发生变化时，两个调压室的水位都将发生变化，而任一个调压室的水位的变化，将引起水轮机流量新的改变，从而影响到另一个调压室的水位的变化，因此两个调压室的水位变化是相互制约的，使整个引水系统的水力现象大为复杂，当引水隧洞的特性和尾水隧洞接近时，可能发生共振。
因此设计上下游双调压室时，不能只限于推求波动的第一振幅，而应该求出波动的全过程，研究波动的衰退情况，但在全弃负荷时，上、下游调压室互不影响，可分别求其最高和最低水位。
第五章调压室
第十章调压室第一节调压室的功用、要求及设置条件一、调压室的功用在较长的压力引水系统中，为了降低高压管道的水击压力，满足机组调节保证计算的要求，常在压力引水道与压力管道衔接处建造调压室。
这样，从水库到调压室为纵向坡度较缓的压力隧洞，其内压力较低，而从调压室到厂房为坡度较陡的高压管道。
有时如果尾水隧洞的长度较大，也可设置尾水调压室。
调压室利用扩大断面和自由水面反射水击波，它将有压引水系统分成两段：上游段为压力引水道，下游段为压力管道。
调压室的功用可归纳为以下几点：(1) 反射水击波。
基本上避免(或减小)了压力管道传来的水击波进入压力引水道；(2) 缩短了压力管道的长度，从而减小了压力管道及厂房过水部分的水击压力；(3) 改善机组在负荷变化时的运行条件；(4) 由于从水库到调压室之间引水道的水压力较低，从而降低了其设计标准，节省了建设经费。
二、调压室的基本要求根据其功用，调压室应满足以下基本要求：(1) 调压室应尽量靠近厂房，以缩短压力管道的长度；(2) 调压室内水体应有自由表面和足够的底面积，以保证水击波的充分反射；(3) 调压室的工作必须是稳定的。
在负荷变化时，引水道及调压室水体的波动振幅小并迅速衰减，达到新的稳定状态；(4) 正常运行时，水流经过调压室底部造成的水头损失要小。
为此调压室底部和压力管道连接处应具有较小的断面积；(5) 结构安全可靠，施工简单方便，造价经济合理。
三、调压室的设置条件调压室是改善有压引水系统和水电站运行条件的一种可靠措施。
但调压室一般尺寸较大，投资较高，工期长，特别是对于低水头电站，调压室在整个引水系统造价中可能占相当大的比例。
因此是否设置调压室，应在机组过流系统调节保证计算和机组运行条件分析的基础上，考虑水电站在电力系统中的作用、地形及地质条件、压力管道的布置等因素，进行技术经济比较后加以确定。
1．上游调压室的设置条件初步分析时，可用水流加速时间(也可称为压力引水道的水流惯性时间常数)T w 来判断，设置上游调压室的条件：[]wp ii w T gH V L T ≥=∑ (10-1) 式中：L i ——引水道(包括蜗壳和尾水管)各段长度，m ；V i ——上述各段引水道的流速，m/s ；H p ——水轮机设计水头，m ；[T w ]——T w 的允许值，一般取2~4s 。
调压室水力计算的电算法简介
调压室水力计算的电算法简介电算法与常用的解析法和图解法相比，具有计算理论严密，简化假设少，速度快，精度高，可以计算不同类型的调压室在各种工况下的涌波全过程并可与水锤、机组转速变化联合求解等许多优点。
尤其在研究某参数对调压室水位变化过程的影响时，电算法更为便利。
进行调压室水位波动计算时，以水轮机、阀门的出流方程作为边界条件，从某种已知初始状态开始，采用四阶龙格-库塔数值积分法求解调压室水流连续方程和隧洞水流动力方程。
本节仅介绍阻抗式、简单式调压室的水位波动计算，给出IBM-PC/XT机FORTRAN语言的计算程序及计算实例。
对程序稍加修改，便可用于某些布置方式较特殊的调压室波动计算。
一、调压室水位波动的基本微分方程调压室的基本方程为：.连续方程2.动力方程式中Q——隧洞中的流量；——压力管道中的流量；F——调压室的截面积；Z——调压室水位；——上游水库水位；K——调压室阻抗水头损失系数；——调压室中的流量，以进人调压室时为正；R——隧洞的沿程损失和局部损失系数；g——重力加速度；A——隧洞的截面积；L——隧洞的长度。
如已知出流变化规律，则，可以根据四阶龙格-库塔法来逐步求解式(15-53)和式(15-54）。
二、龙格-库塔法计算公式如已知t时刻的、值，则可以根据以下公式来求t十△t时刻的、之值。
三、程序框图程序框图如图15-25所示。
图15-25 计算框图四、程序中所用符号的意义NS——调压室高程—截面曲线已知点数。
ZA (I, J)——调压室高程—截面曲线。
共I行，每行二列，分别为一个高程值及相应的调压室截面积。
ZA(I，1)为高程，ZA(I，2)为面积。
HR——上游水库水位，m。
TS——水轮机导叶或阀门关闭或开启时间，s。
KI——水流流进调压室时的阻抗系赦，即为中的系数，单位为。
K0——水流流出调压室时的阻抗系数，。
Q0——起始流量，。
QE——终止流量，。
Z——起始时调压室中的水位，m。
L——从水库到调压室处的引水隧洞长度，m。
差动式调压室甩荷时水位波动以及大井面积和有效阻抗孔面积的计算方法
& ’ B
再设
? ?
水 ! 动 ! 力 ! 学 ! 研 ! 究 ! 与 ! 进 ! 展 !!!!!!!!!" . 年第 $ 期
’ S ’S S! S 2" S 2 2 2 $ $ $ S * S * S * S 2 S * ’ ’"" ! 2 #S 2 ’S 2 $ S 2 " S 2
将上式代入式 ! " 得 B S! " ’"" ’" $"" 2 2 ’1 "2 S 2 上式也可以写成
" " " " " #" $ ) ! ( % #" 1! !- " ! " "! ", + 1 &- " * "
1 !-"! &
!
"
! " # .
! " # 上式即为水位第一 次 向 上 波 动 时 调 压 室 最 高 涌 波 水 位的计算式 * 下面推求当调压室水位上升到最高时 # 调压室第二振幅的计算公式 * ! " # #
$ $ D 6!,ED F G > H F < I : 2% * I > G < I : 27 < > L & G !!4 !!4 J J K
" % $ % # M G < I2 I G N & ’ O G P LC & H F I * + * M G < IA # ? @$ ; F G I < =* J =
! " # B
+& #则由上式得 ,* "% + &" " .$ $ " -" , 1" *
有压隧洞的水力计算(自编)
一、有压隧洞的水力计算1、沿程水头损失：h f =Lv²/(C²R)=λLv²/(d2g)=Ln²Q²/(F²R^4/3)R=A/χi上游调压室的设置条件λ=8g/C²C=R^(1/6)/n2、局部水头损失：hj=ζv²/(2g)3、有压隧洞的基本计算公式：①自由出流：Q=μω√(2g(T 0-h p ))式中，Tw—压力水道中水流惯性h p =0.5a+p ′/γLi—压力水道及蜗壳和压②淹没出流：Q=μω√(2g(T 0-h s ))vi—压力水道内各分段流 Hp—水轮机设计水头，m 4、①自由出流：μ=1/(1+∑ζj *(ω/ωj )^2+∑2gl i *(ω/ωi )^2/(C i ²*R i ))^0.5； [Tw]—Tw 的允许值，一般②淹没出流：μ=1/((ω/ω2)^2+∑ζj *(ω/ωj )^2+∑2gli*(ω/ωi)^2/(C i ²*R i ))^0.5,式中：ω2—隧洞出口下游渠道断面面积 ω—隧洞出口断面面积 ζj —几部水头损失系数ωj —与 ζj 相应流速之断面面积L i 、ωi 、R i 、C i —某均匀洞段之长度、面积、水力半径、谢才系数压力钢管经济直径D=1.128(Q/v e )^0.5= 或 压力钢管经济直径D=(5.2*Q max ^[]w w T T ＞iw i pL vT gH =∑二、阻抗式调压室（一）、托马断面计算：A=K*A th =K*L*A 1/(2g*(α+1/(2g))*(H 0-h w0-3*h wm ))式中：A th —托马临界稳定断面面积 L—压力引水道长度 A 1—压力引水道断面面积H 0—发电最小静水头（电站上下游水位差）α—自水库至调压室水头损失系数，α=h w0/v²，（包括局部水头损失与沿程摩擦水头损失），在无连接管 v—压力引水道流速h w0—压力引水道水头损失 h wm —压力管道水头损失K—系数，一般可采用1.0~1.1（二）、最高涌波计算（《水电站调压室设计规范》计算公式）：A=K*A th =K*L*A 1/(2g*(α+1/(2g))*(H 0-h w0-3*1、阻抗孔水头损失计算：h c =（Q/(Ψs)^2)/(2g)式中： h c —通过阻抗孔的水头损失 S—阻抗孔断面面积0.6~0.8之间选用2、丢弃全负荷时的最高涌波计算（《水电站调压室设计规范》计算公式）：λ′=2gA(h c0+h w0）/(LA 1v 0²)(1+λ′Z max )-ln(1+λ′Z max )=(1+λ′h w 0)-ln(1-λ′h c 0)(λ′|Z max -1|)＋ln(λ′|Z max |－1)=ln(λ′h c 0-1)-(λ′h w 0+1)34、增加负荷时的最低涌波计算：1+(((0.5ε-0.275m ′^0.5)^0.5)+0.1/ε-0.9)×(1-m ′)(1-m ′/(0.65ε^0.62))m ′=Q/Q 03、甩负荷时的第二振幅Z2m′=Q/Q0ε=LA1v0²/(gAh w0²)上游调压室的设置条件式中，Tw—压力水道中水流惯性时间常数，s；i—压力水道及蜗壳和压力尾水道各分段长度，m ；i—压力水道内各分段流速，m/s ；Hp—水轮机设计水头，m ；Tw]—Tw 的允许值，一般取2~4s式中： v e —经济流速，明钢管和地下埋管为4~6m ∕s ；管经济直径D=1.128(Q/v e )^0.5= 3.140219≈3.1 钢筋砼管为2~4m/s ；坝内埋管为3~7m/s 压力钢管经济直径D=(5.2*Q max ^3/H)^(1/7)=3.434174≈3.4Q max —管道的最大流量[]w w T T ＞iw i pL vT gH =∑二、阻抗式调压室水力计算程摩擦水头损失），在无连接管时用α代替（α+1/（2g））A1/(2g*(α+1/(2g))*(H0-h w0-3*h wm))141216441618 m′)(1-m′/(0.65ε^0.62))管为4~6m∕s；埋管为3~7m/s。
水力计算及实例讲解
2、管道水力计算中，设备阻力降要考虑，尤其是低压供 气的情况。 日常工作中需要考虑的设备阻力降包括：流量计、过滤 器、电磁阀及其他一些管径或流向发生巨大变化的管道 设备。这部分阻力降在低压管道中所占比例很大，如忽 略，可能会造成用气设备前的压力、流量满足不了设备 正常燃烧的要求。 3、设备选型计算要重视。 设备选型不光要考虑压力、流量、计量精度满足，同时 还要考虑价格因素，有些人不重视设备选型计算，其实 调压计量设备在整个管网系统中是很重要的环节，且对 工程造价有很大的影响。
目
录
水力计算基本步骤 水力计算基本公式 水力计算注意事项
水力计算实例讲解
水力计算基本步骤
1、与用户协商确定用气地点பைடு நூலகம்设备型号、适用气种、负荷、 用气压力等技术参数（需要书面确认，避免基础数据出问 题）； 2、制定供气方案及路线图，并与用户沟通确认（初步进行设 备选型）； 3、根据管线布置图编号； 4、计算各管段计算流量； 5、初选管径，计算各管段实际阻力损失△P； 6、根据允许阻力降△Pd来调整管段管径； 7、设备选型计算。
式中 Re—雷诺数； △P—燃气管道摩擦阻力损失（Pa）； λ—燃气管道摩擦阻力系数； l—燃气管道的计算长度（m）； Q—燃气管道的计算流量（m3/h）； d—管道内径（mm）； ρ—燃气的密度（kg/m3）； T—设计中所采用的燃气温度（K）； T0—标准状态的温度（K）；
ν—标准状态时燃气的运动粘度（m2/s）； K—管壁内表面的当量绝对粗糙度。
管道允许阻力降△Pd＝0.75Pn＋150 Pn－低压灶具的额定用气压力（Pa），要根据不同气种、 不同灶具来确定。 天然气灶具一般为2000Pa，故△Pd＝1650Pa，旧燃规里 根据经验把1650划分成庭院＋户内各分别占多少帕，新 规范里没有明确提出，只是要求分配时要根据情况，经 技术经济比较后确定。 对于高层民用户采用二次调压供气时，应根据低低压调 压器的进口压力范围来确定一二级调压间管道的允许阻 力降。 高差大时，水力计算中应考虑附加压力的影响。
空调水管水力计算
一、空调水系统的设计原则：1、力求水力平衡；2、防止大流量小温差；3、水输送符合规范要求；4、变流量系统宜采用变频调节；5、要处理好水系统的膨胀与排气；6、解决好水处理与水过滤；7、切勿忽视管网的保冷与保温效果。
二、冷冻水、冷却水管的计算1、压力式水管道管径计算D=103πνL4（mm ）公式中 L------水流量（m 3/s ）v-------计算流速（m/s ）一般水管系统的管内水流速可参考表13-12的推荐值取用表13-13选择。
2、直线管段的阻力计算Δh=d l λ×22v ρ=R ×l 式中Δh---长度为l （m ）的直管段的摩擦阻力（Pa ）λ---水与管内壁间的摩擦阻力系数；l----直管段的长度（m ）；d----管内径（m ）；ρ----水的密度（kg/m 3），当4℃时为1000kg/m 3R-----长度为1m 直管段的摩擦阻力（Pa/m ）三、空调设备流量计算由Q=CM ΔT 可得出：M=Q/C*ΔT （Kg/S ）Q-----空调制冷或制热量（Kw ）C-----水的比热容，4.2KJ/Kg*℃ΔT---进出空调设备的供回水温差，ΔT =T G -T H四、风机盘管选择1、计算室内空调冷负荷Q （W ），简单依单位面积指标及经验估算。
2、考虑机组的盘管用后积垢积尘对传热的影响，对空调冷负荷要进行修正，冷负荷应乘以系数a仅冷却使用 a=1.10作为加热、冷却两用 a=1.20仅作为加热用 a=1.153、依据空调冷负荷选择风机盘，一般按中档运行能力选择。
4、校核风量：L=)(3600s n h h Q -ρ L-----风机盘管名义风量（m 3/h ）Q-----室内空调冷负荷（KW）h n-----室内空气计算温度下空气焓值（KJ/Kg）h s------室内空气送风温度下空气焓值（KJ/Kg）ρ-----空气密度，取标态下1.2Kg/m3五、送风温差1、一般舒适性空调送风温差：送风高度≤5m 送风温差Δt s≤10℃送风高度＞5m 送风温差Δt s≤15℃2、工艺性空调的送风温差：六、集水器的选择：1、通常用到集水器及分水器时水系统至少要分为三个子系统以上才会考虑用之！集水器与分水器的管径，接其中水的流速大致控制在通常情况下0.5~0.8m/s，并应大于最大接管开口直径的二倍。
第七节调压室水力计算条件的选择
第七节调压室水力计算条件的选择第七节调压室水力计算条件的选择调压室的基本尺寸是由水力计算来确定的，水力计算主要包括以下三方面的内容：(1）研究“引水道—调压室”系统波动的稳定性，确定所要求的调压室最小断面积。
(2）计算最高涌波水位，确定调压室顶部高程。
(3）计算最低涌波水位，确定调压室底部和压力水管进口的高程。
进行水力计算之前，需先确定水力计算的条件。
调压室的水力计算条件，除去水力条件之外，还应考虑到配电及输电的条件。
在各种情况中，应从安全出发，选择可能出现的最不利的情况作为计算的条件。
现讨论如下。
1.波动的稳动性计算调压室的临界断面，应按水电站在正常运行中可能出现的最小水头计算。
上游的最低水位一般为死水位，但如电站有初期发电和战备发电的任务，这种特殊最低水位也应加以考虑。
引水系统的糙率是无法精确预侧的，只能根据一般的经验选择一个变化范围，根据不同的设计情况，选择偏于安全的数值。
计算调压室的临界断面时，引水道应选用可能的最小糙率，压力管道应选用可能的最大糙率。
流速水头、水轮机的效率和电力系统等因素的影响，一般只有在充分论证的基础上才加以考虑。
2.最高涌波水位的计算上游水库水位应取正常高水位，引水道的糙率应取可能的最小值，负荷的变化情况一般按丢弃全负荷设计。
最高洪水位丢弃全负荷或部分负荷进行校核。
如电站的机组和出线的回路数较多，而且母线分段，经过分析，电站没有丢弃全负荷的可能，也可不按丢弃全负荷计算。
对于丢弃全负荷情况，可假定由最大流量减小至空转流量；为了安全，有人认为应按丢弃至零计算。
3.最低涌波水位的计算上游水库水位应取可能的最低水位，引水道的糙率应取可能的最大糙率。
确定最不利的增荷情况比确定最不利的丢荷情况更加困难。
增加负荷对调压室的工作比丢弃负荷更危险，如计算不正确，可能使引水道和压力管道进入空气，破坏建筑物和机组正常的运行。
在技术设计阶段，增加负荷的条件应根据设计电站在系统中的工作情况，经专门研究确定，并应考虑系统将来的发展。
空调水系统水力计算的基本原理
供 热 工 程
p zh
2
1 2 2 G A G zh zh 2 9002 2 d 4 2
按式 p A zh G 2 编制的水力计算表。
第十三讲 水暖系统水力计算原理
在工程设计中，对常用的垂直单管顺流式 系统，由于整根立管与干管、支管以及支管与 散热器的连接方式，在施工规范中都规定了标 准的连接图式；
设 a=
1 s
=
G
p
，得出
供 热 工 程
ab a1 a2 a3
ab —并联管段的总通导数，（kg/h）/Pa1/2
2 p s G s G s G 由 1 1 2 2 3 3 得 2 2
G1: G2: G3=
1 s1
：s ：
2
1
1 s3
a1：a2：a3
第十三讲 水暖系统水力计算原理
供 热 工 程
上式表明：在串联管路中，管路的总阻抗为 各串联管段阻力数之和。
第十三讲 水暖系统水力计算原理
2、并联管路的总阻抗
供 热 工 程
类似并联电压，并联管路的总流量为各并联 管段流量之和。
G G1 G2 G3
2 p sG 把 带入上式得出：
第十三讲 水暖系统水力计算原理
三、塑料管材的水力计算
供 热 工 程
式中 R—热媒在计算温度和流量下的比摩阻，Pa/m； RO—实验条件下的的数据，Pa/m； 比摩阻的水温修正系数a。 地暖系统要注意。
值的 分户计量热水供暖系统常用塑料管材， 计算公式是由实验得出的，再求出R值编制成表， 计算时可查表。 表中的数值是水温为10℃时的，当水温升高时， 对其要进行修正，可按下式修正： R=RO· a
用三维计算流体力学方法计算调压室阻抗系数
1 数学模型和计算方法
1.1 流动特点 水电站引水系统的水流主要是湍流，其雷诺数通常很高，水力损失可分为沿程和局部两类。沿 程损失主要由边界的滞水作用所引起，在高雷诺数的阻力平方区，主要决定于边壁粗糙度。局部损失是由于流 动边界急剧改变，引起流态发生激烈变化，出现分离、扩散、回流，加剧了湍流掺混，从而引起能量损失的增 加；其损失系数主要取决于流道几何形状和流场结构(如分离、扩散、回流、射流等)。在局部损失中也包含沿 呈损失的部分，只是沿程部分所占比重很小。从计算流体力学的观点来看，水力损失计算的准确与否，关键在 于湍流计算方法的选用和湍流边壁的处理。阻抗式调压室流场有分离、扩散、回流、射流等流动特性，只有选 用适合这些流动结构的计算方法才能得到正确的水力损失。 1.2 湍流模型 目前计算湍流的方法有雷诺平均的NS方程法、大涡模拟法和直接数值模拟法。后两种由于需要 的计算资源庞大而不适用于本问题计算。在雷诺平均NS方程的方法中，雷诺应力模型和
的平均旋转率张量；系数A0=4.04， As =
1 ~ ~ 6 cos(φ ), φ = cos −1 ( 6W ),W = S ij S jk S ki / S , S = S ij S ij 。 3
1.3 边壁处理 边界层内流速梯度大，剪切应力作用强，边壁区的准确模拟对摩阻系数、压力变化和分离点 位置等有明显影响，因而边壁的处理对水力损失系数的计算很重要。雷诺平均湍流模拟方法是用时间平均后的 NS方程模拟平均流动，用湍流模型(例如k-ε方程)描述湍动，用它模拟带边壁的流场时，必须对边壁加以特别 处理。目前边壁处理方法有两种，其中两层法适用于低雷诺数流动，且要求近壁层内网格足够细；壁函数法适
* * +
2
模拟效果验证
典型阻抗调压室由岔管、阻抗孔和突扩段组成。为验证以上计算摸型，了解模拟结果的精度，这里分别模 拟T形岔管、圆管孔口、突扩管的局部损失系数，以及某模型调压室的阻抗系数。 2.1 T形岔管水头损失系数 带连接管的阻抗式调压室的连接管通常与引水道及压力管道形成T形岔管。T 形岔管的水头损失主要取决于连接管、引水道和压力管道的断面积比、水流的分流/汇流比、三管连接处的修 圆半径等。若连接管并不垂直于引水道和压力管道，损失系数还与它们的夹角有关。选用图1所示的等直径 (D1=D2=D3=1.0m)T形岔管，给定不同的分流比Q2/Q1和汇流比Q2/Q3进行模拟，得到相应的损失系数(以总管流速水 头为基准)。图1是分流情况，由图1可见，计算结果和经验值变化趋势一致，数值基本接近。随分流比的增大，
空调管道水力计算解读
的计算阻力相等，称为并联管路阻力平衡。对 一般的通风系统，两支管的计算阻力差应不超 过15%；含尘风管应不超过10%。若过上述规 定，采用下述方法进行阻力平衡。 （1）调整支管管径 这种方法通过改变支管管径来调整支管阻力， 达到阻力平衡。调整后的管径按下式计算：
含有粉尘等，会增加设备和管道 的磨损。反
之，流速低，阻力小，动力消耗少；但是风 管断面大，材料和建造费用大，风管占用的 空间也增大。流速过低会使粉尘沉积而堵塞 管道。因此， 必须通过全面的技术经济比较选定合理的流 速。根据经验总结，风管内的空气流速可按 下表确定。若输送的是含尘气流，流速不应 低于其表所列的值。
计算，获得管网特性曲线。
图
圆形伞形罩 800m3/s
1 L=11m
2 L=6m 3 5 L=3m L=4m 4 L=6m
6 7 L=6m L=8m
除尘器
风机
1500m3/s 4000m3/s
图2-3-2 通风除尘系统的系统图
[ 解 ]：
1.对各管段进行编号，标出管段长度和风点
的排风量。 2.选定最不利环路，本系统选择1-3-5-除尘器 - 6-风机-7为最利环路。 3.根据各管段的风量及选定的流速，确定最 不利环路各管段的断面尺寸和单位长度摩擦 阻力。 根据表2-2-3输送含有轻矿物粉尘的空气时， 风管内最小风速为，垂直风管12m/s、水平 风管14m/s.
考虑到除尘器及风管漏风，取5%的漏风
系数，管段6及7的计算量为
6300×1.05=6615m3/h. 管段1 有水平风管，初定流速为14m/s。根据 Q1=1500m/h(0.42m3/s)、V1=14m/s所选管 径按通风管道 统一规格调整为 D1=200mm：实际流速V1=13.4m3/S;由图2-31查得，Rml=12.5Pa/m。 同理可查得管段3、5、6、7的管径及比摩阻， 具体结果见表2-3-5。 4.确定管段2、4的管径及单位长度摩擦力，
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
调压室的水力计算1. 调压室断面计算当上游死水位，下游为最低水位，最小水位H min=188.9m，三台机满发，引水道糙率取最小值，压力管道糙率取最大值，通过水轮机的流量为57m3s⁄，则此时的引水隧洞水头损失的计算如表格1，压力钢管水头损失的计算如表格2。
引水道应选可能的最小糙率0.012，压力管道应选择可能的最大糙率0.013。
表格1引水隧洞水头损失表表格2压力钢管水头损失表F Tℎ>Lf2αgH1=Lf2αg(H0−ℎw0−3ℎwT)=45.548m2其中H0——最小水头损失，H0=188.9m；ℎw0——引水隧洞损失，ℎw0=17.802+0.296=18.098；ℎwT0——压力管道水头损失，ℎwT=3.110+2.805=5.915m；L——引水隧洞长度，12662m；g——重力加速度，g=9.81m/s2 f——引水隧洞面积，16.62m2。
α——引水道阻力系数v0=Qf=5716.619=3.43m s⁄α=ℎw0v02=18.0983.4302=1.5385为了保证大波动的稳定，一般要求调压室断面大于托马斯断面，初步分析时可取（1.0~1.1）F Tℎ，作为调压室的设计断面。
这里选取D=7.8m，则系数k为：F k=47.784k=F k/F Tℎ=1.052. 最高涌波水位计算按正常蓄水位时共用同一调压室的三台机组全部满载运行瞬时丢弃全部负荷（即流量由Q max=57减至流量Q=0）作为设计工况。
引水隧洞的糙率取尽可能的最小值（能耗少，涌波高）。
n=0.012引水道损失由表格1和表格2得：ℎw0=ℎw0程+ℎw0局=17.802+0.296=18.098mv0为时段开始时管中流速v0=Qf=3.43m s⁄；f为引水隧洞断面面积。
F为调压井断面面积，145.267m2；引水隧洞长L=12662m，g=9.81m s2⁄得引水道—调压室系统的特性系数。
λ=Lfv022gFℎw0=12662×16.62×3.4322×9.81×47.784×18.098=145.89令X0=ℎw0λ=0.124，X=zλ，则要求最高涌波水位z max，只需要求出X max=z maxλ即可。
X max的符号在静水位以上为负，以下为正。
ln（1+X max）−X max=−X0运用牛顿切线法求解方程的根令φ(x)=ln（1+x）−x+X0牛顿迭代公式为：x k+1=x k−φ(x)φ′(x)取迭代初值x0=−0.5，计算结果见下表表格3迭代计算结果由表格3可以看出精确到0.001，X max=z maxλ=−0.419|z max|=−λX max=61.128m根据《水电站调压室设计规范》调压室最高涌波水位以上的安全超高不宜小于1m。
所以调压室的顶高程：Z=1279+61.128+1.5=1341.63m3. 最低涌波水位计算3.1 增加负荷时的最低涌波水位上游水位取死水位，下游取最低水位。
引水隧洞糙率取最大值，取n=0.016。
水轮发电机组由2台增至3台满负荷运行。
由表1和表2引水隧洞水头损失为：ℎw0=ℎw0程+ℎw0局=31.649+0.296=31.945m|z min |ℎw0=1+(√ε−0.275√m +0.05ε⁄−0.9)(1−m)(1−m ε0.62)⁄ 式中ε=2λℎw0⁄=9.134，m = 23 ，解得：z min =50.44m所以当增加负荷时水位降到静水位下23.74m 3.2 丢弃全负荷时产生的第二波动振幅X 2+ln (1−X 2)=X max +ln (1−X max )解之X 2+ln (1−X 2)=−0.069X 2=0.327z min =λX 2=47.70m所以当丢弃全负荷时，水位降到静水位下47.70m 综上所述调压室的最低涌波水位应为静水位下50.44m 根据《水电站调压室设计规范》上游调压室最低涌波水位与调压室处压力引水道顶部之间的安全高度应不小于2～3m 。
所以压力引水道顶部高程为：Z=1274-50.44-2.5=1221.06m4. 四阶龙格库塔法4.1调压室水位的基本微分方程dz dt =(Q−Q m)F=f1(t,z,Q)dQ=(H R−z−KQ s|Q s|−RQ|Q|)gA=f2(t,z,Q)已知压力管道的流量变化规律，则调压室流量等于隧洞流量减压力管道流量，即：Q s=Q−Q m。
上述方程两个未知量，Q和Z，采用四阶龙格库塔法进行逐步求解。
4.2四阶龙格库塔法计算公式如果已知t时刻的Q t，Z t值，采用以下公式求解t+∆t时刻的Q t+∆t，Z t+∆t。
Z t+∆t=Z t+16(K1+2K2+2K3+K4)K1=∆tf1(t,Z t,Q t)K2=∆tf1(t+∆t2,Z t+K12,Q t+L12)K3=∆tf1(t+∆t,Z t+K2,Q t+L2)K4=∆tf1(t+∆t,Z t+K3,Q t+L3)Q t+∆t=Q t+16(L1+2L2+2L3+L4)L1=∆tf2(t,Z t,Q t)L2=∆tf2(t+∆t2,Z t+K12,Q t+L12)L3=∆tf2(t+∆t2,Z t+K22,Q t+L22)L4=∆tf2(t+∆t,Z t+K3,Q t+L3)4.3程序框图4.4 VBA程序代码定义函数f1：Public Function f1(t As Double, z As Double, q As Double)Dim a As Doublea = 5 'a为导叶直线关闭时间If t < a Thenf1 = (q - (57 - 57 * t / a)) / 47.784Elsef1 = q / 47.784End IfEnd Function定义函数f2：Public Function f2(t As Double, z As Double, q As Double)Dim r As Doubler = 0.0055 'r为隧洞的沿程损失和局部损失系数f2 = 9.81 * 16.61 * (1279 - z - r * Abs(q) * q) / 12662End Function调压室甩负荷时水位波动程序Sub tyssw()Dim t As DoubleDim z As DoubleDim q As Double '将t,z,q定义为实数Dim k1, k2, k3, k4 As DoubleDim l1, l2, l3, l4 As Double '将k1, k2, k3, k4,l1, l2, l3, l4定义为实数Dim dt, stoptime As Double 'dt为时间步长，stoptime为计算的调压室水位波动的时间，单位（s）stoptime = 200 '计算从导叶开始关闭到t=200s的水位波动过程i = 2 'i用于计算结果计数dt = 1t = 0q = 57z = 1260.9 't,q,z稳态时的初始值Do While t < stoptime '当时间t小于stoptime时执行以下循环过程k1 = dt * f1(t, z, q)l1 = dt * f2(t, z, q)'k2 = dt * f1(t + dt / 2, z + k1 / 2, q + l1 / 2)l2 = dt * f2(t + dt / 2, z + k1 / 2, q + l1 / 2)'k3 = dt * f1(t + dt / 2, z + k2 / 2, q + l2 / 2)l3 = dt * f2(t + dt / 2, z + k2 / 2, q + l2 / 2)'k4 = dt * f1(t + dt, z + k3, q + l3)l4 = dt * f2(t + dt, z + k3, q + l3)'t = t + dtz = z + (k1 + 2 * k2 + 2 * k3 + k4) / 6q = q + (l1 + 2 * l2 + 2 * l3 + l4) / 6 '如上是四阶龙格库塔法计算过程i = i + 1Cells(i, 1) = tCells(i, 2) = zCells(i, 3) = qLoop '循环结束标志End Sub调压室增加负荷时水位波动程序主程序相同，只需要将初始值改为Q=38，Z=1259.8，其中自定义函数改为如下：Public Function f1(t As Double, z As Double, q As Double)Dim a As Doublea = 5 'a为导叶直线关闭时间If t < a Thenf1 = (q - (38 + 19 * t / a)) / 47.784Elsef1 = (q -57)/ 47.784End IfEnd Function定义函数f2：Public Function f2(t As Double, z As Double, q As Double)Dim r As Doubler = 0.001 'r为隧洞的沿程损失和局部损失系数f2 = 9.81 * 16.61 * (1279 - z - r * Abs(q) * q) / 12662End Function4.5计算结果及调压室水位振荡曲线甩负荷增负荷 57减至0 57减至38 57减至19 19增至3819增至5738增至57最高上升水位 相应水位 61.14 1340.14 52.22 1331.22 57.32 1336.32 最低下降水位 相应水位-47.73 1231.27-46.47 1232.53-47.25 1231.75-22.86 1251.14-48.28 1225.72-30.95 1243.051245125012551260126512701275128012850200400600800100012001400增负荷时调压室水位振荡曲线122012401260128013001320134013600200400600800100012001400甩负荷时调压室水位振荡曲线Z min =1231.27T s =400sZ max =1340.14图 1流量由57减至0调压室水位振荡曲线图 2流量由38增至57调压室水位振荡曲线Z min =1243.05T s =200s5. 解析法与四阶龙格库塔法对比表格4计算结果对比单位（m）在计算最低涌波水位时。