CN109724820B - 一种水电机组在线监测***的水轮机工作效率计算方法 - Google Patents
一种水电机组在线监测***的水轮机工作效率计算方法 Download PDFInfo
- Publication number
- CN109724820B CN109724820B CN201811627446.1A CN201811627446A CN109724820B CN 109724820 B CN109724820 B CN 109724820B CN 201811627446 A CN201811627446 A CN 201811627446A CN 109724820 B CN109724820 B CN 109724820B
- Authority
- CN
- China
- Prior art keywords
- water turbine
- unit
- efficiency
- working
- characteristic curve
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Control Of Water Turbines (AREA)
Abstract
本发明提出的一种水电机组在线监测***的水轮机工作效率计算方法,其特征在于,包括以下步骤:测量水电站的工作参数;获取水轮机工作区域完整的真机特性曲线;基于真机特性曲线、水轮机工作方程、单位参数方程迭代求解水轮机工作水头和工作效率。本发明不需要知道水轮机的工作水头,只需要水轮机除水头以外的基本工作参数和机组特性曲线,即可计算水轮机效率。采用机组运行实测数据即可计算水轮机效率,计算过程简单,易于在实际运行机组上实现水轮机内部能量变化的实时监测。在计算水轮机效率的同时得到水轮机工作点,能评估当前机组在特性曲线上的运行区域,便于及时控制机组避开不利区域运行。
Description
技术领域
本发明涉及水力发电能效分析研究领域,尤其涉及一种水电机组在线监测***的水轮机工作效率计算方法。
背景技术
水轮机将水流势能转化为旋转机械能,带动发电机旋转将机械能转化为电能,是水电站中的核心设备。水轮发电机组的效率取决于水轮机和发电机的综合效率特性。其中发电机的效率特性比较稳定、随时间变化较小,一般效率都在95%以上。而水轮机的设计效率虽然可达到95%以上,但运行过程中因容积损失、水力损失和机械损失的影响,使其偏离高效率区,同时泥沙的磨损和空蚀会导致水力性能变差,水轮机的设计选型不合理和制造安装中的技术缺陷都会降低水轮机效率。机组长期在低效率区运行,不仅会影响水能的利用,还会加重空蚀导致的破坏,甚至引起机组振动,导致机组发生严重故障,威胁水电站的安全运行。因此需要对水轮发电机组进行能效分析,以及时了解机组运行状态、排除故障和制定检修计划,确保水电机组安全、高效地运行。
在水轮机能效分析过程中水轮机效率按下式计算:
ηsl-sh=Pdq/gHQfdηfd
其中,Pdq为发电机实际功率;Qfd为机组实际发电用水流量,由超声波流量计获取;ηfd为机组在当前水头和出力下的发电机效率,根据发电机效率曲线获取。
H为水轮机工作水头,按下式计算:
其中,Pwk为蜗壳进口压力,kPa;Pws为尾水管出口压力,kPa;Zwk为蜗壳进口压力表高程,m;Zws为尾水管出口压力表高程,m;Dgg为流量计所在位置钢管直径,m;Qfd为机组实际发电流量,m3/s;γ为水的重度,kg/m3。
在测量蜗壳进口压力和尾水管出口压力时,因传感器安装位置不当、水轮机内部水力性能改变使测点处压力分布发生变化或转轮内水流搅动产生的水击现象,会导致压力测量不准确,直接影响能效分析的结果,在工程实际效率监测中甚至会出现水轮机效率大于1的情况。受制于水电站的运行条件,在已投入运行的水电站中不宜改变压力传感器的安装位置或添加新的传感器,因此研究在工作水头测量失误时也能够对水轮机效率进行有效评价的方法很有必要。
发明内容
本发明正是针对现有技术存在的不足,提供了一种无需改变压力传感器安装位置或添加新的传感器,在工作水头测量结果失真时也能够对水轮机效率进行有效评价的水电机组在线监测***的水轮机工作效率计算方法。
为解决上述问题,本发明所采取的技术方案如下:
一种水电机组在线监测***的水轮机工作效率计算方法,包括以下步骤:
S1.测量水电站的工作参数;其中工作参数包括发电机功率Pdq、发电机效率ηdq、水轮机工作转速nr、水轮机发电流量Qfd、转轮进口直径D1;
S2.获取水轮机工作区域完整的真机特性曲线;
S3.基于真机特性曲线、水轮机工作方程、单位参数方程迭代求解水轮机工作水头和工作效率;
其中,S3具体步骤如下:
S31.根据水轮机工作方程P=gQHη、发电机效率方程P=Pdq/ηdq计算ηsHs=Pdq/ηdqgQfd,其中,ηs为水轮机实际工作效率,Hs为水轮机实际工作水头;
S33.假设水轮机当前效率是ηl,带入步骤S31和S32中的方程式,得到水轮机工作水头H、力矩M、单位转速n′1s、单位流量Q′1s和单位力矩M′1s;
S34.将单位转速n′1s、单位力矩M′1s在水轮机单位参数特性曲面上插值,得到与之对应的单位流量Q′1l;
S35.根据水轮机工作方程P=gQHη和单位流量方程在功率P和流量Q一定的情况下,水轮机工作水头H和效率η反相关、单位流量Q′1和水头H的平方根反相关,单位流量Q′1和效率η正相关,根据单位流量Q′1s和Q′1l的差别对效率进行修正,得到
优选地,在S2中,若水轮机运行区域内的真机特性曲线不完整,采用与当前机组对应的模型特性曲线对当前真机特性曲线进行修正,直到得到完整、连续的真机特性曲线。
优选地,在S2中,若水轮机运行区域内的真机特性曲线不完整,采用与当前机组比转速相近、具有相同工作区的机组特性曲线对当前真机特性曲线进行修正,直到得到完整、连续的真机特性曲线。
优选地,S34的具体步骤如下:
对水轮机特性曲线进行描点,将图转换成数据;利用三次B样条函数插值法对三维数据进行均匀采样,形成综合特性B样条曲线:
以单位转速、单位流量、水轮机效率为坐标,导叶开度为参变量,参数方向取u向为沿开度线方向,v向为开度变化方向,参数u为相对累积弦长,参数v为相对开度,二者在[0,1]区间单调递增,将均匀采样后的三组B样条曲线转换成B样条曲面:
其中,m为沿开度线方向的数据点个数,n为开度线条数,对边缘部分数据线性外延,生成单位转速-单位流量-效率曲面;
根据水轮机工作方程P=gQHη、力矩-功率关系方程P=M2πn/60和单位参数方程求得带入单位转速-单位流量-效率曲面中的数据,得到水轮机单位力矩,生成单位转速-单位流量-单位力矩曲面,即水轮机单位参数特性曲面;
将单位转速n′1s、单位力矩M′1s在水轮机单位参数特性曲面上插值,得到与之对应的单位流量Q′1l。
本发明与现有技术相比较,本发明的实施效果如下:
(1)本发明不需要知道水轮机的工作水头,只需要水轮机除水头以外的基本工作参数和机组特性曲线,即可计算水轮机效率,计算量小。
(2)采用机组运行实测数据即可计算水轮机效率,计算过程简单,易于在实际运行机组上实现水轮机内部能量变化的实时监测。
(3)在计算水轮机效率的同时得到水轮机工作点,能评估当前机组在特性曲线上的运行区域,便于及时控制机组避开不利区域运行。
附图说明
图1为本发明的真机特性曲线图。
图2为本发明的单位转速-单位流量-效率曲面图。
图3为本发明的水轮机单位参数特性曲面图。
图4为本发明的计算流程图。
图5为本发明的效率计算过程中单位流量Q′1、工作水头H和效率η变化图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整性地描述。当然,所描述的实施例只是本发明一部分实施例,而不是全部的实施例。
本发明提出的一种水电机组在线监测***的水轮机工作效率计算方法,包括以下步骤:S1.测量水电站的工作参数;S2.获取水轮机工作区域完整的真机特性曲线;S3.基于真机特性曲线、水轮机工作方程、单位参数方程迭代求解水轮机工作水头和工作效率。
下面结合具体步骤对本实施例做进一步说明。
水轮机综合特性曲线处理
结合图1,将真机特性曲线转换成数据,内插采用样条曲线法插值,外延采用线性插值法延拓;将真机特性曲线处理成单位转速-单位流量-效率的曲面。
具体处理过程为:首先对水轮机特性曲线进行描点,将图转换成数据。如果原始特性曲线参数以数据形式提供则省去这一步。然后利用三次B样条函数插值法对三维数据进行均匀采样,形成综合特性B样条曲线:
将均匀采样后的三组B样条曲线转换成B样条曲面。转换过程以单位转速、单位流量、水轮机效率为坐标,导叶开度为参变量,参数方向取u向为沿开度线方向,v向为开度变化方向,参数u为相对累积弦长,参数v为相对开度,二者在[0,1]区间单调递增。
全特性B样条曲面:
其中,m为沿开度线方向的数据点个数,n为开度线条数。对边缘部分数据的外延,可以直接线性延拓。
处理后的(n′1,Q′1,η)曲面如图2所示。
水轮机效率计算
在准确测量发电机功率Pdq、发电机效率ηdq、水轮机工作转速nr、水轮机发电流量Qfd、转轮进口直径D1等参数后,本发明中水轮机效率计算方法按如下过程进行:
步骤一.根据水轮机工作方程P=gQHη、发电机效率方程P=Pdq/ηdq计算ηsHs=Pdq/ηdqgQfd,其中,ηs为水轮机实际工作效率,Hs为水轮机实际工作水头;
步骤三.假设水轮机当前效率是ηl,带入前述方程式确定水轮机工作水头H、力矩M、单位转速n′1s、单位流量Q′1s和单位力矩M′1s;
步骤四.根据单位转速n′1s、单位力矩M′1s在水轮机单位参数特性曲面图3上插值,求得对应的单位流量Q′1l;
步骤五.根据水轮机工作方程P=gQHη和单位流量方程可知,在功率P和流量Q一定的情况下,水轮机工作水头H和效率η反相关、单位流量Q′1和水头H的平方根反相关,因此单位流量Q′1和效率η正相关,根据单位流量Q′1s和Q′1l的差别对效率进行修正,
为了验证本发明的计算效果,下面结合实例进行说明。
某水轮机的基本参数如下:
假设在额定工况下,水轮机效率下降至85%,在额定水头、额定出力和额定转速下,水轮机引用流量为54.43m3/s。
在缺失实测水轮机工作水头的情况下,假设水轮机当前效率为100%,参考该水轮机真机特性曲线,根据前述步骤对水轮机效率进行计算,当迭代精度为10e-4时,迭代5次,效率计算结果为85.79%,误差为0.93%,水轮机工作点(n′1,Q′1,M′1)在图3中的坐标为(83.152,1.0134,99.838)。迭代计算中工作水头H、单位参数Q′1s及效率η变化过程如图5所示。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种水电机组在线监测***的水轮机工作效率计算方法,其特征在于,包括以下步骤:
S1.测量水电站的工作参数;其中工作参数包括发电机功率Pdq、发电机效率ηdq、水轮机工作转速nr、水轮机发电流量Qfd、转轮进口直径D1;
S2.获取水轮机工作区域完整的真机特性曲线;
S3.基于真机特性曲线、水轮机工作方程、单位参数方程迭代求解水轮机工作水头和工作效率;
其中,S3具体步骤如下:
S31.根据水轮机工作方程P=gQHη、发电机效率方程P=Pdq/ηdq计算ηsHs=Pdq/ηdqgQfd,其中,ηs为水轮机实际工作效率,Hs为水轮机实际工作水头;
S33.假设水轮机当前效率是ηl,带入步骤S31和S32中的方程式,得到水轮机工作水头H、力矩M、单位转速n′1s、单位流量Q′1s和单位力矩M′1s;
S34.将单位转速n′1s、单位力矩M′1s在水轮机单位参数特性曲面上插值,得到与之对应的单位流量Q′1l;
S35.根据水轮机工作方程P=gQHη和单位流量方程在功率P和流量Q一定的情况下,水轮机工作水头H和效率η反相关、单位流量Q′1和水头H的平方根反相关,单位流量Q′1和效率η正相关,根据单位流量Q′1s和Q′1l的差别对效率进行修正,得到
2.根据权利要求1所述的一种水电机组在线监测***的水轮机工作效率计算方法,其特征在于,在S2中,若水轮机运行区域内的真机特性曲线不完整,采用与当前机组对应的模型特性曲线对当前真机特性曲线进行修正,直到得到完整、连续的真机特性曲线。
3.根据权利要求1所述的一种水电机组在线监测***的水轮机工作效率计算方法,其特征在于,在S2中,若水轮机运行区域内的真机特性曲线不完整,采用与当前机组比转速相近、具有相同工作区的机组特性曲线对当前真机特性曲线进行修正,直到得到完整、连续的真机特性曲线。
5.根据权利要求1所述的一种水电机组在线监测***的水轮机工作效率计算方法,其特征在于,S34的具体步骤如下:
对水轮机特性曲线进行描点,将图转换成数据;利用三次B样条函数插值法对三维数据进行均匀采样,形成综合特性B样条曲线:
以单位转速、单位流量、水轮机效率为坐标,导叶开度为参变量,参数方向取u向为沿开度线方向,v向为开度变化方向,参数u为相对累积弦长,参数v为相对开度,二者在[0,1]区间单调递增,将均匀采样后的三组B样条曲线转换成B样条曲面:
其中,m为沿开度线方向的数据点个数,n为开度线条数,对边缘部分数据线性外延,生成单位转速-单位流量-效率曲面;
根据水轮机工作方程P=gQHη、力矩-功率关系方程P=M2πn/60和单位参数方程求得带入单位转速-单位流量-效率曲面中的数据,得到水轮机单位力矩,生成单位转速-单位流量-单位力矩曲面,即水轮机单位参数特性曲面;
将单位转速n′1s、单位力矩M′1s在水轮机单位参数特性曲面上插值,得到与之对应的单位流量Q′1l。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811627446.1A CN109724820B (zh) | 2018-12-28 | 2018-12-28 | 一种水电机组在线监测***的水轮机工作效率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811627446.1A CN109724820B (zh) | 2018-12-28 | 2018-12-28 | 一种水电机组在线监测***的水轮机工作效率计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109724820A CN109724820A (zh) | 2019-05-07 |
CN109724820B true CN109724820B (zh) | 2020-09-22 |
Family
ID=66295963
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811627446.1A Active CN109724820B (zh) | 2018-12-28 | 2018-12-28 | 一种水电机组在线监测***的水轮机工作效率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109724820B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112113620A (zh) * | 2020-09-15 | 2020-12-22 | 甘肃盛御水利水电科技有限公司 | 一种水电厂流量在线监测*** |
CN112525567B (zh) * | 2020-09-29 | 2022-11-08 | 华能澜沧江水电股份有限公司 | 一种检测水轮发电机各稳定性区域累计运行时长的方法 |
CN113987802B (zh) * | 2021-10-28 | 2023-04-11 | 华中科技大学 | 一种水泵水轮机特性系数的计算方法及*** |
CN114417516B (zh) * | 2021-12-09 | 2023-02-14 | 中国电建集团华东勘测设计研究院有限公司 | 一种可逆式水轮机组特性曲线的获取方法、装置、设备及存储介质 |
CN115935851A (zh) * | 2022-12-23 | 2023-04-07 | 中国电建集团成都勘测设计研究院有限公司 | 用于水电站复杂管网恒定流仿真的改进全局梯度方法 |
CN116593046B (zh) * | 2023-07-17 | 2023-10-13 | 河海大学 | 水轮机下降效率确定方法、装置及电子设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101672239A (zh) * | 2009-09-24 | 2010-03-17 | 昆明理工大学 | 水轮机内部能量损失的动态监测方法 |
CN103306886A (zh) * | 2013-05-29 | 2013-09-18 | 郑程遥 | 一种水轮发电机组的全参数调节控制方法 |
WO2016171590A1 (ru) * | 2015-04-24 | 2016-10-27 | Анатолий Михайлович КРИШТОП | Регулируемая гидромуфта с блокировкой |
CN107191308A (zh) * | 2017-07-25 | 2017-09-22 | 中国水利水电科学研究院 | 一种混流式水泵水轮机全特性曲线的预测方法 |
CN108153998A (zh) * | 2018-01-25 | 2018-06-12 | 哈尔滨工业大学 | 离心鼓风机叶轮的全三维气动优化设计方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508159B (zh) * | 2011-11-17 | 2013-09-25 | 河海大学 | 大型发电机组自动调节装置通用辅助测试***及方法 |
CN104533701B (zh) * | 2014-12-23 | 2017-03-08 | 华中科技大学 | 一种水轮机调速***控制参数的自动整定方法 |
-
2018
- 2018-12-28 CN CN201811627446.1A patent/CN109724820B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101672239A (zh) * | 2009-09-24 | 2010-03-17 | 昆明理工大学 | 水轮机内部能量损失的动态监测方法 |
CN103306886A (zh) * | 2013-05-29 | 2013-09-18 | 郑程遥 | 一种水轮发电机组的全参数调节控制方法 |
WO2016171590A1 (ru) * | 2015-04-24 | 2016-10-27 | Анатолий Михайлович КРИШТОП | Регулируемая гидромуфта с блокировкой |
CN107191308A (zh) * | 2017-07-25 | 2017-09-22 | 中国水利水电科学研究院 | 一种混流式水泵水轮机全特性曲线的预测方法 |
CN108153998A (zh) * | 2018-01-25 | 2018-06-12 | 哈尔滨工业大学 | 离心鼓风机叶轮的全三维气动优化设计方法 |
Non-Patent Citations (1)
Title |
---|
基于双向流固耦合的贯流式水轮机动力特性分析;郑小波等;《农业工程学报》;20160223;第78-79页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109724820A (zh) | 2019-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109724820B (zh) | 一种水电机组在线监测***的水轮机工作效率计算方法 | |
CN110821746B (zh) | 一种转桨式水轮机变转速运行的桨叶协联方法 | |
CN108301955B (zh) | 轴流式pat发电模式最优效率点参数和性能曲线预估方法 | |
CN111859746B (zh) | 一种基于流场重构的透平机械变工况性能预测方法 | |
Rai et al. | Analyzing hydro abrasive erosion in Kaplan turbine: A case study from India | |
CN103926079B (zh) | 一种混流式水电机组出力异常检测方法 | |
Batista et al. | Self-start evaluation in lift-type vertical axis wind turbines: Methodology and computational tool applied to asymmetrical airfoils | |
CN113446146A (zh) | 一种在线水轮机效率试验方法 | |
Muratoğlu et al. | Numerical analyses of a straight bladed vertical axis Darrieus wind turbine: Verification of DMS algorithm and Qblade code | |
CN110610034B (zh) | 一种混流式水轮机水力特性的生成方法 | |
Thakker et al. | Experimental investigation of CA9 blades on a 0.3 m wells turbine rig | |
CN107515996B (zh) | 一种达里厄水轮机导流罩型线的优化设计方法 | |
KR102343779B1 (ko) | 수력 에너지를 변환하기 위한 수력학적 기계 및 시설의 작동점을 결정하기 위한 방법 | |
Ivetić et al. | Novel discharge measurement system at the turbine intakes of Iron Gate 2 hydropower plant: a system description | |
CN117539168B (zh) | 一种基于半实物仿真的水轮机空化诊断***及方法 | |
Fernando et al. | Torque performance analysis of a hydrokinetic turbine in a tank set-up using CFD | |
CN111307493B (zh) | 基于知识的塔式太阳能熔盐储热***故障诊断方法 | |
CN113153453B (zh) | 汽轮机末级叶片容积流量估计方法、颤振预警方法及***和装置 | |
CN113673062B (zh) | 水斗式水轮机转轮外径、水斗外宽及转轮重量的计算方法 | |
Quishpe et al. | Evaluation of guide vanes effect over runner Francis turbine sediment erosion using a quasi-two dimensional approach | |
Li et al. | Preliminary results of a vortex method for stand-alone vertical axis marine current turbine | |
Wu et al. | A Rapid 3D-Modeling Method for Impeller of Wind Turbine | |
Kumar et al. | Influence of guide vane clearance on internal flow of medium-specific speed Francis turbine | |
CN115559862A (zh) | 一种海上风力监测方法及*** | |
Adamkowski et al. | Performance and diagnostic tests on hydraulic gensets in Polish hydropower plants |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |