调压室水力计算
(4)--水电站调节保证计算,调压室考核试卷及答案
《水电站》考核(四)学生姓名：班级学号：浙江水利水电学院水电站课程组编制2013年8月(修改)使 用 说 明本考核是《水电站》课程形成性考核的依据，与文字教材配套使用。
考核作业是课程考核的重要组成部分，是强化教学管理，提高教学质量，反馈学习信息，提高学生综合素质和能力的重要保证。
通过形成性考核有助于学生理解和掌握本课程的基本概念、基本理论。
同时，形成性考核对于全面测评学生的学习效果，督促和激励学生完成课程学习，培养学生自主学习和掌握知识的能力也具有重要作用。
全部课程要求完成5次计分考核。
学生应按照教学进度按时完成各次计分考核，教师根据学生完成的情况评定成绩，每次作业以100分计，并按5次考核的平均成绩计算学生的形成性考核成绩。
考核成绩占课程总成绩的20% 。
考核四说明：本部分覆盖引水系统水力计算(水锤和调压室)部分，在学完本单元课程后，先完成与本单元相关的题目， 待学完本模块所有内容后，全部完成此次考核。
一、判断题(20分)1．导叶的关闭时间Ts愈大，水锤压力愈大，机组转速升率愈小。
2．对高水头电站，一般可采用先快后慢的机组关闭规律，以达到降低水锤的目的。
( )3．对低水头电站，一般可采用先慢后快的机组关闭规律，以达到降低水锤的目的。
( )4．水电站甩负荷时，初始开度越大，水锤压力就越大。
起始开度越小，水锤压力越小。
( )5．调压室底部流速对调压室的稳定有利。
( )6．阻抗式调压室的阻抗越大越好。
( )7．调压室离进水口越近，则其水位波动幅值越小，故调压室应当尽量靠近进水口。
( )8．水头愈低，需要的调压室稳定断面越小。
( )9．压力钢管的糙率对调压室的稳定断面没有影响。
( )10．调压室越靠近厂房时，会使波动稳定断面减小。
( )(二)填空题(40分)1．延长机组关闭时间可以使__________减小,但___________将会增大。
2．极限水锤沿管道的分布规律为____ ___第一相水锤为 ____ ____。
调压室水力计算分解
调压室的水力计算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λ即可。
第二节 调压室的工作原理和基本方程
第二节调压室的工作原理和基本方程一、调压室的工作原理水电站在运行时负荷会经常发生变化。
负荷变化时，机组就需要相应地改变引用流量，从而在引水系统中引起非恒定流现象。
压力管道中的非恒定流现象（即水锤现象）在上一章中已经加以讨论。
引用流量的变化，在“引水道-调压室”系统中亦将引起非恒定流现象，这正是本节要加以讨论的。
图13-5为一具有调压室的引水系统。
当水电站以某一固定出力运行时，水轮机引用的流量亦保持不变，因此通过整个引水系统的流量均为，调压室的稳定水位比上游水位低，为通过引水道时所造成的水头损失。
当电站丢弃全负荷时，水轮机的流量由变为零，压力管道中发生水锤现象，压力管道的水流经过一个短暂的时间后就停止流动。
此时，引水道中的水流由于惯性作用仍继续，流向调压室，引起调压室水位升高，使引水道始末两端的水位差随之减小，因而其中的流速也逐渐减慢。
当调压室的水位达到水库水位时，引水道始末两端的水位差等于零，但其中水流由于惯性作用仍继续流向调压室，使调压室水位继续升高直至引水道中的流速等于零为止，此时调压室水位达到最高点。
因为这时调压室的水位高于水库水位，在引水道的始末又形成了新的水位差，所以水又向水库流去，即形成了相反方向的流动，调压室中水位开始下降。
当调压室中水位达到库水位时，引水道始末两端的压力差又等于零，但这时流速不等于零，由于惯性作用，水位继续下降，直至引水道流速减到零为止，此时调压室水位降低到最低点。
此后引水道中的水流又开始流向调压室，调压室水位又开始回升。
这样，引水道和调压室中的水体往复波动。
由于摩阻的存在，运动水体的能量被逐渐消耗，因此，波动逐渐衰减，最后全部能量被消耗掉，调压室水位稳定在水库水位。
调压室水位波动过程见图13-5中右上方的一条水位变化过程线。
当水电站增加负荷时，水轮机引用流量加大，引水道中的水流由于惯性作用，尚不能立即满足负荷变化的需要，调压室需首先放出一部分水量，从而引起调压室水位下降，这样室库间形成新的水位差，使引水道的水流加速流向调压室。
水道水力过渡过程计算大纲范本
FJD34270 FJD水利水电工程技术设计阶段水道水力过渡过程计算大纲范本水利水电勘测设计标准化信息网1998年3月1水电站技术设计阶段水道水力过渡过程计算大纲主编单位:主编单位总工程师:参编单位:主要编写人员:软件开发单位:软件编写人员:勘测设计研究院年月2目次1. 引言 (4)2. 设计依据文件和规范 (4)3.计算基本资料 (4)4.大波动水力过渡过程计算 (7)5.小波动水力过渡过程计算 (15)6.专题研究（必要时） (16)7.应提供的设计成果 (17)3１引言＿＿抽水蓄能电站位于＿＿，在电力系统中的功能是＿＿。
电站总装机容量＿＿MW，单机容量＿＿MW。
机组型号＿＿。
电站开发方式(首部开发、中部开发、尾部开发) ＿＿。
引水系统由＿＿组成。
本工程为＿＿等工程。
可行性研究报告于＿＿年＿＿月审查通过。
2 设计依据文件和规范2.1 有关本工程文件(1) 工程可行性研究报告；(2) 工程可行性研究报告审批文件；(3) 技术设计任务书。
2.2 主要设计规范(1) SD 303—88 水电站进水口设计规范(2) SD 144—85 水电站压力钢管设计规范(试行)(3) DL/T 5058-1996 水电站调压室设计规范(4) SDJ 173—85 水力发电厂机电设计技术规范(5) SD 134—84 水工隧洞设计规范(试行)(6) GB 9652—88 水轮机调速器与油压装置技术条件2.3 参考资料和手册《水电站机电设计手册》(水力机械部分)。
3 计算基本资料3.1 水库(水池)特征水位(1)上库(上水池)水位：正常蓄水位＿ｍ；死水位＿ｍ。
(2) 下库(下水池)水位：正常蓄水位＿ｍ；死水位＿ｍ。
4提示：对于混合式抽水蓄能电站，尚应补充上、下库设计、校核洪水位。
3.2 引水系统布置(1)引水系统平面布置(2)引水系统纵剖面布置(3)引水系统特征参数，见表1表1 引水系统特征参数表管道编号部位直径m面积m2长度m管道末端高程m水头损失系数，×Q2备注局部水头损失沿程水头损失水轮机工况水泵工况最大值平均值最小值1…注：（1）引水系统编号示意图，可表示在上表备注栏中。
2023年注册土木工程师(水利水电)之专业知识提升训练试卷A卷附答案
2023年注册土木工程师（水利水电）之专业知识提升训练试卷A卷附答案单选题（共30题）1、一般情况下，工农业用水中还原计算的主要项目是( )。
A.生活耗水量B.农业灌溉C.工业用水量D.蒸发水量【答案】 B2、每一个需要安置的农业人口的安置补助费标准每亩被征用耕地的安置补助费，最高不得超过被征用前三年平均年产值的( )倍。
A.10B.15C.20D.25【答案】 B3、多泥沙河流上的进水闸(分洪闸、分水闸)闸槛高程可比( )。
A.河(渠)底略低B.河(渠)底高许多C.河(渠)底略高D.A或B【答案】 C4、按山区、丘陵区水利水电工程永久性水工建筑物洪水标准确定水工建筑物洪水标准的条件为( )。
A.挡水高度低于15m，且上下游水头差小于10m时B.挡水高度低于15m，或上下游水头差小于10m时C.挡水高度大于15m，或上下游水头差小于10m时D.挡水高度大于15m，且上下游水头差大于10m时【答案】 D5、导流建筑物级别划分时，当失事后果为淹没基坑，但对总工期及第一台(批)机组发电影响不大，经济损失较小时，导流建筑物的级别为( )级。
A.2B.3C.4D.5【答案】 D6、( )是指单位时段内通过流域出口断面的泥沙总量和相应集水面积的比值。
A.输移比B.容许流失量C.输沙模数D.土壤侵蚀模数【答案】 C7、地下水资源量时，一般将评价范围划分为( )和( )两大类，分别采用( )和( )计算地下水资源量。
A.平原区、山丘区、排泄量法、补给量法B.平原区、高原区、补给量法、排泄量法C.平原区、山丘区、补给量法、排泄量法D.平原区、高原区、排泄量法、补给量法【答案】 C8、溢流坝的反弧段应结合( )选择。
A.直线段B.孔口泄流的抛物线C.幂曲线D.下游消能形式【答案】 D9、《中华人民共和国土地管理法》规定，“征用耕地的土地补偿费，为该耕地被征用前三年平均年产值的六至十倍”。
这里的“该耕地”，是指( )。
水力计算书
水力计算书已知条件：1、居民用户：55户2、总室外楼前架空管长度4米；室外楼前架空管1长度65米； 室外楼前架空管2长度75米；3、中压管道长度150米4、4个单元5层5、立管长度12米6、每户额定耗气量：2.7m ³/h （标）7、中压管道设计压力：0.4MPa8、中压管道运行压力：0.2MPa9、低压管道设计压力：5000Pa 10、低压管道运行压力：2400Pa 一、计算流量Q h1=∑kNQn=0.375*0.7*55+0.75*2*55=96.94m ³/h （标） Q h2=∑kNQn=0.40*0.7*30+0.75*2*30=53.40m ³/h （标） Q h3=∑kNQn=0.43*0.7*25+0.75*2*25=45.03m ³/h （标） Q h4=∑kNQn=0.68*5*0.7+0.92*5*2=11.58m ³/h （标） 由理想气体状态方程可得工况下的流量分别为：()h m T P T P Q Q h I /99.34273200325.101293325.10194.96321121=⨯+⨯⨯==中()h m T P T P Q Q h I /63.1012734.2325.101293325.10194.963211211=⨯+⨯⨯==()h m T P T P Q Q h I /45.552734.2325.101293325.1014.533211222=⨯+⨯⨯==()h m T P T P Q Q h I /76.462734.2325.101293325.10103.453211233=⨯+⨯⨯==()h m T P T P Q Q h I /02.122734.2325.101293325.10158.113211244=⨯+⨯⨯==二、计算管道管径由v D Q 436002π＝工可得vQ D I⨯⨯⨯=π36004（由《化工工艺设计手册》查得中压流速范围为10~20m/s ，本设计取15m/s ；低压流速范围为8~12m/s ，本设计取8m/s ）中压管道内径：mm m vQ D I 29029.01514.3360099.34436004==⨯⨯⨯=⨯⨯⨯=π中查《低压流体输送用焊接钢管》得公称口径DN25的焊接钢管外径为33.7mm ，壁厚为3.2mm ，33.7-3.2=30.5mm ＞29mm 。
1.10室内给水系统的水力计算
谢谢欣赏
THANK YOU FOR WATCHING
（4）有两条或两条以上具有不同最大用水时卫 生器具给水当量平均出流概率给水支管的给水 干管，该管段的最大时卫生器具给水当量平均 出流概率按(1-10)式计算：
1.10.2.2集体宿舍、旅馆、宾馆、疗养院、
幼儿园、养老院、办公楼、商场、客运站、会展 中心、中小学教学楼、公共厕所等建筑
1.10.2.3按卫生器具同时作用系数确 定设计秒流量
工业企业生活间、公共浴室、洗衣房、公共食堂、 实验室、影剧院、体育场等建筑的生活给水管道设计 秒流量，应根据卫生器具给水额定流量、同类型卫生 器具数和卫生器具的同时给水百分数下式计算
1.10.3给水管路的水力计算
1.10.3.1 设计流速 流速不要过大或过小，设计时应综合考虑
以上因素，把流速控制在适当的范围内。生活 或生产给水管道的水流速度宜按表1-20采用。
1.10.3.3管道水头损失的计算
l 给水管道的水头损失包括沿程损失和局部损失 两部分
1、沿程水头损失
ห้องสมุดไป่ตู้
1.10.3.3管道水头损失的计算
2.局部水头损失
1.10.3.3管道水头损失的计算
在实际工作中，给水管网的局部水头损失 一般不作详细计算，而是按给水系统的沿程水 头损失的百分数采用，其值如下： l 生活给水管网——25%～30%； l 生活、消防共用给水管网——20%； l 生产、生活、消防共用给水管网——20%； l 生产、消防共用给水管网——15%； l 消火栓系统消防给水管网——10%； l 自动喷水灭火系统消防给水管网——20%。
（1）计算最大用水时卫生器具给水当量平均出 流概率：
住宅建筑设计秒流量的计算
水电站调压室设计规范
水电站调压室设计规范Specification for design of surge chamber of hydropowerstation中华人民共和国电力行业标准水电站调压室设计规范主编部门：电力工业部华东勘测设计研究院批准部门：中华人民共和国电力工业部中华人民共和国电力工业部关于发布《水电站调压室设计规范》电力行业标准的通知电技［1996］733号各电管局，各省、自治区、直辖市电力局，水电水利规划设计总院，各有关单位：《水电站调压室设计规范》电力行业标准，经审查通过，批准为推荐性标准，现予发布。
其编号为：DL/T5058－1996该标准自1997年5月1日起实施。
请将执行中的问题和意见告水电水利规划设计总院，并抄送部标准化领导小组办公室。
1996年10月31日目次1总则2术语、符号3调压室的设置条件及位置选择4调压室的基本布置方式、基本类型及选择5调压室的水力计算及基本尺寸的确定6抽水蓄能电站调压室的设计7调压室的结构设计、构造、观测及运行要求附录A压力水道水头损失计算公式附录B调压室的涌波计算公式附录C抽水蓄能电站水泵工况断电、导叶拒动时的调压室涌波计算方法本规范用词规定附加说明1总则1.0.1水电站调压室是压力水道系统中一项重要建筑物，为体现国家现行的技术经济政策，积极慎重地采用国内外先进技术和经验，统一调压室设计的标准、要求，特制定本规范。
1.0.2本规范适用于大、中型水利水电枢纽工程中常规水电站和抽水蓄能电站调压室设计，小型水电站的调压室设计可参照执行。
1.0.3水电站调压室设计应根据地形、地质情况、压力水道的布置、机电特性和运行条件等资料，经综合论证，做到因地制宜、经济合理、安全可靠。
1.0.4水电站调压室设计除必须遵守本规范的规定外，还应符合SDJ12—78《水利水电枢纽工程等级划分及设计标准(山区、丘陵区部分)》(试行)及补充规定、SD134—84《水工隧洞设计规范》、SDJ173—85《水力发电厂机电设计技术规范》(试行)、DL/T5057—1996《水工混凝土结构设计规范》、SDJ10—78《水工建筑物抗震设计规范》(试行)等现行的国家、行业有关标准与规定。
- 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 dt =(H R−z−KQ s|Q s|−RQ|Q|)gAL=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+∆t2,Z t+K22,Q t+L22)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）在计算最低涌波水位时。