CN111365158B - 水轮机转轮实时状态评估和生命周期管理预测系统 - Google Patents

水轮机转轮实时状态评估和生命周期管理预测系统 Download PDF

Info

Publication number
CN111365158B
CN111365158B CN202010135942.6A CN202010135942A CN111365158B CN 111365158 B CN111365158 B CN 111365158B CN 202010135942 A CN202010135942 A CN 202010135942A CN 111365158 B CN111365158 B CN 111365158B
Authority
CN
China
Prior art keywords
runner
data
real
rotating wheel
unit
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
Application number
CN202010135942.6A
Other languages
English (en)
Other versions
CN111365158A (zh
Inventor
张猛
郭靖
文树洁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dongfang Electric Machinery Co Ltd DEC
Original Assignee
Dongfang Electric Machinery Co Ltd DEC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Dongfang Electric Machinery Co Ltd DEC filed Critical Dongfang Electric Machinery Co Ltd DEC
Priority to CN202010135942.6A priority Critical patent/CN111365158B/zh
Publication of CN111365158A publication Critical patent/CN111365158A/zh
Application granted granted Critical
Publication of CN111365158B publication Critical patent/CN111365158B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03BMACHINES OR ENGINES FOR LIQUIDS
    • F03B11/00Parts or details not provided for in, or of interest apart from, the preceding groups, e.g. wear-protection couplings, between turbine and generator
    • F03B11/008Measuring or testing arrangements
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/20Hydro energy
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P70/00Climate change mitigation technologies in the production process for final industrial or consumer products
    • Y02P70/50Manufacturing or production processes characterised by the final manufactured product

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Artificial Intelligence (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Computational Linguistics (AREA)
  • Biomedical Technology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • Hydraulic Turbines (AREA)

Abstract

本发明属于水轮发电机组状态监测技术领域,公开了一种水轮机转轮运行状态实时健康评估和全生命周期管理及预测系统,根据电站机组转轮的固有特性、转轮制造信息、机组安装信息、机组实时运行数据、机组过流水质及泥砂含量等综合信息对转轮的运行健康状态进行实时评估,并在显示台上实时数字化显示转轮运行状态信息、转轮的运行健康度、转轮止漏环实时间隙值、转轮劣化程度与疲劳寿命、转轮预测检修周期、转轮运行评估报告、机组运维优化建议等信息,为水电机组的安全、稳定、经济运行、合理检修和故障诊断提供定量的科学决策支持。

Description

水轮机转轮实时状态评估和生命周期管理预测系统
技术领域
本发明专利属于水轮发电机组状态监测与健康评估技术领域,具体涉及为一种水轮机转轮实时状态评估和生命周期管理预测系统。
背景技术
水轮发电机组的转轮对整个机组的安全经济稳定运行至关重要,目前不少电站均出现转轮磨损及汽蚀严重、转轮叶片裂纹、转轮止漏环碰磨等问题。而且,多数水电机组已接入电网AGC,即机组自动发电控制,AGC的投入使得水电机组运行方式发生较大变化,即多台机组长时间带部分负荷运行、机组负荷调节频繁、负荷变动幅度大且调节速度快、机组频繁跨越振动区和机组开停机频繁,加剧转轮叶片的空蚀程度、恶化转轮叶片受力情况并诱发转轮疲劳裂纹。由于目前转轮在线监测和评估技术的限制,多数机组只能在定期检修时才能发现上述问题,这给机组的安全、经济、稳定运行带来很大的隐患。
目前,行业内针对转轮空蚀和裂纹尚无法直接监测,现有技术多集中于流道压力监测、流道压力脉动监测、超声监测、噪音监测、振动摆度监测以及声发射技术等间接监测手段,上述方法均存在以下主要难题:1)故障样本缺乏,信号噪音复杂,有效特征信号提取困难;2)滞后性,如振动摆度监测、噪音监测和声发射技术,只有当转轮出现裂纹或裂纹发展到一定程度后,才能监测到相关特征信号,无法有效的指导电站优化运行和预知性检修。
目前行业内针对转轮的安全评估,多基于合同约束值,针对机组运行若干特定工况对转轮进行刚度、强度、动力特性和疲劳寿命分析。对其疲劳寿命的评估一般基于合同所约束的且极为简化的疲劳荷载谱进行相关计算,由于无法获取真机实时运行工况数据,其计算结果无法满足转轮真实疲劳损伤程度的评估需求。
如在现有技术中,公开号为CN105678025A,公开日为2016年6月15日,名称为“基于动应力测试和稳定性试验水轮机运行优化方法及系统”的发明专利文献,公开了一种基于动应力测试和稳定性试验水轮机运行优化方法及系统,首先结合水轮机模型分析在各运行区域下的运行特点;获取水轮机真机运行状态实时数据,通过水电机组在线监测系统分析统计运行时的状态量;通过对水轮机转轮的静应力进行有限元分析,核查转轮在各种工况下最大静应力数值和部位;然后对水轮机进行动应力测试,分析和评估动态应力和转轮运行寿命;测试水轮机全水头、全负荷段稳定性试验,最后获取水轮机各重要部位的参数,评估出机组运行的振动区和稳定运行区;最大限度的保证机组运行安全,避免转轮裂纹产生延长检修周期,提高水轮机的安全稳定性,保证机组的安全经济运行。该专利方案存在以下不足之处:1)转轮寿命评估仅考虑实测动应力,实际影响因素较多,如机组全流道水力设计、转轮固有特性、转轮材料、制造工艺和残余应力等;2)影响转轮运行健康状态的指标包括空蚀与磨蚀,该专利未涉及和未引入相关的评估参数,如机组过流水质信息和泥砂含量等;3)未涉及转轮止漏环间隙的安全评估。
水轮机组在实际运行中,其工况多复杂多变,机组压力脉动、汽蚀、泥砂磨损、转轮制造质量、机组安装质量、机组轴系运行状态等因素对转轮的振动、应力和疲劳寿命等均有较大影响,且转轮在线流体仿真和在线结构有限元仿真技术尚未可行,目前尚无相关系统能综合考虑上述因素对转轮的运行状态进行实时评估和生命周期管理及预测。
发明内容
本发明提供了一种根据电站机组转轮的固有特性、转轮制造质量、机组安装质量、机组实时运行数据、泥砂含量等综合信息对转轮的运行状态进行实时自动评估和对转轮疲劳寿命进行定期自动评估,并在显示台上实时显示转轮运行状态的优劣程度、转轮剩余疲劳寿命、推荐检修期等信息,为电站机组的安全运行与合理检修提供定量的科学决策的转轮实时运行状态健康评估和全生命周期管理及预测系统。
本发明的这种水轮机转轮实时状态评估和生命周期管理预测系统,其特征在于:包括数据输入端、数据处理单元和数据输出端;所述数据输入端包括机组实时运行数据输入端和转轮固有特性信息数据输入端;
所述机组实时运行数据输入端连接有若干用于监测水轮机组实施运行数据和水文环境数据的实时数据监测单元,即机组实时运行数据输入端是向系统传输采集到的重要参数对应运转状态和环境变化的实时数据的,而所述转轮固有特性信息数据输入端用于输入对应所述水轮机组实时运行数据参数下的机组机械状态数据,以及对应所述水文环境数据参数下的机械特性受环境影响变化数据,即转轮固有特性信息数据是水轮机组根据设计计算、试验记录、历史分析结果等记录方式得到的机组机械特性在受到重要参数影响的历史曲线数据;
所述数据处理单元接收所述数据输入端的输入数据,并将机组实时运行数据输入端的输入数据与转轮固有特性信息输入端的输入数据进行比对,找出转轮固有特性信息数据中对应水轮机组实时运行数据和机械特性受环境影响变化数据的状态,对包括转轮止漏环间隙、转轮疲劳寿命和转轮运行状态进行评估,即将实时监测到的关键数据去对应寻找其在历史曲线数据中对机组机械特性的影响结果,然后通过所述数据输出端输出;
所述数据输出端包括信息显示单元。
具体的,所述机组实时运行数据输入端连接有用于监测机组运行时的工况类数据、状态类监测数据以及机组过流水质泥砂含量和机组过流水质PH值的实时数据监测模块;其中所述工况类数据包括水头、出力、效率和尾水位数值;所述状态类监测数据包括流道压力、全流道压力脉动、机组振动摆度监测数据、推力以及导轴承瓦温数据。
所述电所述站水质的泥砂含量是根据多年每月份统计值确定的。
进一步的,所述转轮固有特性信息数据输入端输入的机组机械状态数据包括机组机械数据、水轮机转轮运转特性曲线和检修历史数据以及机组安装数据;
所述机组机械数据包括弹性模量、材质成分、屈服强度、抗拉强度和S-N曲线的转轮材质信息,和包括转轮尺寸、不平衡重量、表面粗糙度和残余应力的转轮设计与制造数据,以及电站水质和泥砂对转轮磨损影响的相关系数和理论预测曲线;
所述轮机转轮运转特性曲线包括等效率线、等开度线、吸出高度等值线、等空化系数线、压力脉动等值线、出力限制线、转轮叶片转角等值线、关键位置处的平均应力以及关键位置处的动态应力和变形曲线;
所述检修历史数据包括转轮历次故障和修复信息;
所述机组安装数据包括机组安装时的轴系摆度数据、转轮与顶盖止漏环间隙安装值、转轮与底环止漏环间隙安装值。
对应的,所述机组实时运行数据输入端输入的机组机械状态数据也包括机组机械数据、水轮机转轮运转特性曲线和检修历史数据以及机组安装数据;
且所述机组实时运行数据输入端输入的等效率线、等开度线、吸出高度等值线、等空化系数线、压力脉动等值线、出力限制线、转轮叶片转角等值线,是将若干转轮运行工况Sijk汇同水轮机模型转轮综合特性曲线以及真机稳定性实验结果数据,经过数据读取和比对计算得到的,其中,设H为机组运行水头,N为机组出力,Z为机组尾水位,等效率线ηij=η(Hi,Nj)、等开度线αij=α(Hi,Nj)、吸出高度等值线Hsij=Hs(Hi,Nj)、等空化系数线σij=σ(Hi,Nj)、出力限制线(优选的如混流式机组5%出力限制)P5%-ij=P5%(Hi,Nj)、压力脉动等值线
Figure GDA0002959309210000041
转轮叶片转角等值线θij=θ(Hi,Nj),其中Hi、Nj、Zk可根据机组运行工况按任意规则选取并组合,优选的,
Figure GDA0002959309210000042
Hmax为最大水头、Hmin为最小水头、i=0,1,2,...,n1,n1+1为Hi样本数;
Figure GDA0002959309210000043
Nmax为最大出力、Nmin为最小出力、j=0,1,2,...,n2,n2+1为Nj样本数;
Figure GDA0002959309210000044
Zmax为最大尾水位、Zmin为最小尾水位、k=0,1,2,...,n3,n3+1为Zk样本数。
所述机组实时运行数据输入端输入的轮机转轮运转特性曲线中关键位置处的平均应力,是将若干转轮运行工况Sijk经过全流道CFD稳态分析和转轮结构FEM稳态分析得到转轮流道水压力稳态分布和转轮稳态平均应力分布,然后通过对转轮稳态平均应力分布进行提取和计算,得到转轮关键位置处在工况Sijk下的平均应力σm-ijk,基于机器学习或深度学习的智能算法对平均应力σm-ijk进行拟合得到关键位置处的平均应力数据库;σm为转轮关键位置处的平均应力,σm=f(H,N,Z),其中H为机组运行水头,N为机组出力,Z为机组尾水位;Sijk=g(Hi,Nj,Zk),其中Hi、Nj、Zk可根据机组运行工况按任意规则选取并组合,优选的如上。
通过CFD流体仿真和FEA结构仿真,可得到转轮关键位置处在Sijk工况下的平均应力σm-ijk,然后,基于机器学习或深度学习的智能算法,可实现σm=f(H,N,Z)这一复杂非线性函数较为精确的拟合,优选的,可使用人工神经网络建模技术,模型的精度取决于样本Sijk的分布情况和总数量。人工神经网络是一种应用类似于大脑神经突触联接的结构进行信息处理以实现分类的数学模型,其优势在于可以任意精度逼近任意函数;神经网络方法本身属于非线性建模技术,能够适用于各种复杂的数据关系;神经网络具备很强的学习能力,比很多分类方法更好的适用数据空间的变化。
所述机组实时运行数据输入端输入的轮机转轮运转特性曲线中关键位置处的动态应力和变形曲线,是将若干转轮运行工况Sijk经过转轮流固耦合FSI瞬态分析获得转轮动应力分布和轮转变形分布情况,然后通过转轮动应力试验修正,对成果进行提取得到转轮关键位置处在Sijk工况下的动应力Δσijk和相关位置处的变形rijk,所述相关位置处的变形rijk包括转轮上止漏环处的变形值r1 ijk=F1(Hi,Nj,Zk)和\或下止漏环处的变形值r2 ijk=F2(Hi,Nj,Zk),然后基于机器学习或深度学习的智能算法对动态应力Δσijk和动态变形rijk进行拟合得到关键位置处的动态应力数据库、止漏环处的动态变形数据库;其中,H为机组运行水头,N为机组出力,Z为机组尾水位,Hi、Nj、Zk可根据机组运行工况按任意规则选取并组合,优选如上。
所述数据处理单元对转轮止漏环间隙的运行状态进行评估,是根据所述机组实时运行数据输入端输入的转轮与顶盖止漏环间隙安装实测值
Figure GDA0002959309210000051
转轮与底环止漏环间隙安装实测值
Figure GDA0002959309210000052
若干转轮运行工况Sijk经过流固耦合FSI瞬态分析得到并通过机器学习或深度学习智能建模的变形值数据库,包括:转轮止漏环处的实时变形值数据库r1 ijk=F1(Hi,Nj,Zk)和r2 ijk=F2(Hi,Nj,Zk)、顶盖止漏环处的变形值R1 ijk=F3(Hi,Nj,Zk)、底环漏环处的变形值R2 ijk=F4(Hi,Nj,Zk);其中,H为机组运行水头,N为机组出力,Z为机组尾水位,Hi、Nj、Zk可根据机组运行工况按任意规则选取并组合,优选如上;结合所述转轮固有特性信息数据输入端输入的转轮与顶盖止漏环间隙安装值
Figure GDA0002959309210000053
转轮与底环止漏环间隙安装值
Figure GDA0002959309210000054
以及机组安装时的轴系摆度数据为转轮中心处的轴系轨迹
Figure GDA0002959309210000055
得到转轮运行时的上止漏环实时间隙值
Figure GDA0002959309210000056
和转轮运行时的下止漏环实时间隙值
Figure GDA0002959309210000061
水轮机止漏环结构类型较多,止漏间隙一般为0.6~2.5mm之间,现有技术方案中尚无有效监测技术可对运行中转轮与固定件之间的止漏环间隙进行在线监测,本方案基于流体和结构有限元仿真,并结合机组实时运行数据与振摆监测数据,可实现运行中的转轮止漏环间隙的在线评估,可对转轮碰磨故障进行在线预警和趋势分析,并为因止漏环间隙不均而诱发的机组振动提供定量的科学依据。
所述数据处理单元对转轮疲劳寿命的评估,包括转轮快速裂纹安全评估和转轮疲劳裂纹评估;
所述转轮快速裂纹安全评估是转轮结构与机组激励源共振导致的快速裂纹,尤其是叶片出水边卡门涡频率共振,是通过将激振力理论计算频率和激振力电站实测频率与转轮叶片在空气中的固有频率和转轮叶片在水环境中的固有频率进行比较,判断是否共振;
所述转轮疲劳裂纹评估是转轮关键危险点处长期承受动载荷而导致的高周疲劳裂纹,是由转轮固有特性、机组运行工况数据、机组启停信息、仿真计算和/或转轮动应力试验得到的转轮危险点应力-时间荷载曲线,结合转轮材料S-N曲线和历史检修及修复信息,通过雨流荷载计数和疲劳累积损伤理论,对转轮的疲劳寿命进行预测。
智能水电站的水轮机转轮应是基于数字化设计和数字化制造的,数字化设计和数字化设计的转轮与水电站实际运行中的转轮是“数字孪生”型式的,能够完全真实的反映电站运行中转轮的实际工作状态。与现有技术相比,本发明的技术方案综合考虑了水轮机转轮运行状态和使用寿命的各影响因素,从设计、制造、安装、运行等方面进行了考虑,为电站的长期安全、经济、稳定运行提供了定量的科学决策依据,为电站合理安排检修周期和降低检修工作量,实现水轮机的状态检修提供了重要决策支撑。
附图说明
本发明的前述和下文具体描述在结合以下附图阅读时变得更清楚,其中:
图1为本发明一种优选的系统结构关系框图;
图2为本发明一种优选的转轮实时应力数据库构建框图;
图3为本发明一种优选的转轮实时动应力和变形数据库构建框图;
图4为本发明一种优选的数字化机组综合运转特性数据库构建框图;
图5为本发明一种优选的顶盖和底环止漏环变形数据库构建框图;
图6为本发明一种优选的转轮止漏环实时动态间隙预测框图;
图7为本发明一种优选的转轮疲劳寿命预测框图。
具体实施方式
下面通过几个具体的实施例来进一步说明实现本发明目的技术方案,需要说明的是,本发明要求保护的技术方案包括但不限于以下实施例。
实施例1
作为本发明一种最基本的实施方案,如图1,公开了一种水轮机转轮实时状态评估和生命周期管理预测系统,包括数据输入端、数据处理单元和数据输出端;所述数据输入端包括机组实时运行数据输入端和转轮固有特性信息数据输入端。
其中,所述机组实时运行数据输入端连接有若干用于监测水轮机组实施运行数据和水文环境数据的实时数据监测单元,即机组实时运行数据输入端是向系统传输采集到的重要参数对应运转状态和环境变化的实时数据的,而所述转轮固有特性信息数据输入端用于输入对应所述水轮机组实时运行数据参数下的机组机械状态数据,以及对应所述水文环境数据参数下的机械特性受环境影响变化数据,即转轮固有特性信息数据是水轮机组根据设计计算、试验记录、历史分析结果等记录方式得到的机组机械特性在受到重要参数影响的历史曲线数据。
以及,所述数据处理单元接收所述数据输入端的输入数据,并将机组实时运行数据输入端的输入数据与转轮固有特性信息输入端的输入数据进行比对,找出转轮固有特性信息数据中对应水轮机组实施运行数据和机械特性受环境影响变化数据的状态,对包括转轮止漏环间隙、转轮疲劳寿命和转轮运行状态进行评估,即将实时监测到的关键数据去对应寻找其在历史曲线数据中对机组机械特性的影响结果,然后通过所述数据输出端输出。
而所述数据输出端包括信息显示单元。信息显示单元可以显示的信息包括制造厂家、投运时间等转轮出厂信息,最近一次的转轮止漏环实际安装值的转轮安装信息,包括水头、出力、流量、效率、空蚀系数、转轮应力、止漏环变形、压力脉动、机组振动级别、转轮实时运行健康程度信息,以及转轮累计疲劳损伤值、转轮检修周期预测值等;还可以显示至少包括转轮在各运行区域的时间统计图、转轮关键位置处的应力-时间曲线以及转轮剩余疲劳寿命预测值在内信息的运行评估报告,以及运维优化建议。
实施例2
作为本发明一种优选地实施方案,在上述实施例1的基础上,进一步,具体的,所述机组实时运行数据输入端连接有用于监测机组运行时的工况类数据、状态类监测数据以及机组过流水质泥砂含量和机组过流水质PH值的实时数据监测模块;其中所述工况类数据包括水头、出力、效率和尾水位数值;所述状态类监测数据包括流道压力、全流道压力脉动、机组振动摆度监测数据、推力以及导轴承瓦温数据。而所述电所述站水质的泥砂含量是根据多年每月份统计值确定的。
在进一步的,所述转轮固有特性信息数据输入端输入的机组机械状态数据包括机组机械数据、水轮机转轮运转特性曲线和检修历史数据以及机组安装数据;
所述机组机械数据包括弹性模量、材质成分、屈服强度、抗拉强度和S-N曲线的转轮材质信息,和包括转轮尺寸、不平衡重量、表面粗糙度和残余应力的转轮设计与制造数据,以及电站水质和泥砂对转轮磨损影响的相关系数和理论预测曲线;
所述轮机转轮运转特性曲线包括等效率线、等开度线、吸出高度等值线、等空化系数线、压力脉动等值线、出力限制线、转轮叶片转角等值线、关键位置处的平均应力以及关键位置处的动态应力和变形曲线;
所述检修历史数据包括转轮历次故障和修复信息;
所述机组安装数据包括机组安装时的轴系摆度数据、转轮与顶盖止漏环间隙安装值、转轮与底环止漏环间隙安装值。
对应的,所述机组实时运行数据输入端输入的机组机械状态数据也包括机组机械数据、水轮机转轮运转特性曲线和检修历史数据以及机组安装数据;
且如图4,所述机组实时运行数据输入端输入的等效率线、等开度线、吸出高度等值线、等空化系数线、压力脉动等值线、出力限制线、转轮叶片转角等值线,是将若干转轮运行工况Sijk汇同水轮机综合运转特性曲线以及真机稳定性试验结果数据,经过数据读取和比对计算得到的,其中,设H为机组运行水头,N为机组出力,Z为机组尾水位,等效率线ηij=η(Hi,Nj)、等开度线αij=α(Hi,Nj)、吸出高度等值线Hsij=Hs(Hi,Nj)、等空化系数线σij=σ(Hi,Nj)、出力限制线(优选的如混流式机组5%出力限制)P5%-ij=P5%(Hi,Nj)、压力脉动等值线
Figure GDA0002959309210000091
转轮叶片转角等值线θij=θ(Hi,Nj),其中Hi、Nj、Zk可根据机组运行工况按任意规则选取并组合,优选的,
Figure GDA0002959309210000092
Hmax为最大水头、Hmin为最小水头、i=0,1,2,...,n1,n1+1为样本数;
Figure GDA0002959309210000093
Nmax为最大出力、Nmin为最小出力、j=0,1,2,...,n2,n2+1为样本数;
Figure GDA0002959309210000094
Zmax为最大尾水位、Zmin为最小尾水位、k=0,1,2,...,n3,n3+1为样本数。
并且,如图2,所述机组实时运行数据输入端输入的轮机转轮运转特性曲线中关键位置处的平均应力,是将若干转轮运行工况Sijk经过全流道CFD稳态分析和转轮结构FEM稳态分析得到转轮流道稳态分布和转轮稳态平均应力分布,然后通过对转轮稳态平均应力分布进行计算得到转轮关键位置处在工况Sijk下的平均应力σm-ijk,基于机器学习或深度学习的智能算法对平均应力σm-ijk进行拟合得到关键位置处的平均应力数据库;σm为转轮关键位置处的平均应力,σm=f(H,N,Z),其中H为机组运行水头,N为机组出力,Z为机组尾水位;Sijk=g(Hi,Nj,Zk),其中Hi、Nj、Zk可根据机组运行工况按任意规则选取并组合,优选的如上。
通过CFD流体仿真和FEA结构仿真,可得到转轮关键位置处在Sijk工况下的平均应力σm-ijk,然后,基于机器学习或深度学习的智能算法,可实现σm=f(H,N,Z)这一复杂非线性函数较为精确的拟合,优选的,可使用人工神经网络建模技术,模型的精度取决于样本Sijk的分布情况和总数量。人工神经网络是一种应用类似于大脑神经突触联接的结构进行信息处理以实现分类的数学模型,其优势在于可以任意精度逼近任意函数;神经网络方法本身属于非线性建模技术,能够适用于各种复杂的数据关系;神经网络具备很强的学习能力,比很多分类方案更好的适用数据空间的变化。
再者,如图3,所述机组实时运行数据输入端输入的轮机转轮运转特性曲线中关键位置处的动态应力和变形曲线,是将若干转轮运行工况Sijk经过转轮流固耦合FSI瞬态分析获得转轮动应力分布和轮转变形分布情况,然后通过转轮动应力试验修正,并对计算成果进行提取,得到转轮关键位置处在Sijk工况下的动应力Δσijk和相关位置处的变形rijk
所述关位置处的动应力Δσijk=f1(Hi,Nj,Zk),所述相关位置处的变形rijk包括转轮上止漏环处的变形值r1 ijk=F1(Hi,Nj,Zk)和\或下止漏环处的变形值r2 ijk=F2(Hi,Nj,Zk),然后基于机器学习或深度学习的智能算法对动态应力Δσijk和动态变形rijk进行拟合得到关键位置处的动态应力数据库、动态变形数据库;其中,H为机组运行水头,N为机组出力,Z为机组尾水位,Hi、Nj、Zk可根据机组运行工况按任意规则选取并组合,优选如上。
如图5,所述数据处理单元对转轮止漏环间隙运行状态的评估是在线评估,是根据所述机组实时运行数据输入端输入的转轮与顶盖止漏环间隙安装实测值
Figure GDA0002959309210000101
转轮与底环止漏环间隙安装实测值
Figure GDA0002959309210000102
若干转轮运行工况Sijk经过流固耦合FSI瞬态分析得到并通过机器学习或深度学习智能建模的变形值数据库,包括:转轮止漏环处的实时变形值数据库r1 ijk=F1(Hi,Nj,Zk)和r2 ijk=F2(Hi,Nj,Zk)、顶盖止漏环处的变形值R1 ijk=F3(Hi,Nj,Zk)、底环漏环处的变形值R2 ijk=F4(Hi,Nj,Zk);其中,H为机组运行水头,N为机组出力,Z为机组尾水位,Hi、Nj、Zk可根据机组运行工况按任意规则选取并组合,优选如上。结合所述转轮固有特性信息数据输入端输入的转轮与顶盖止漏环间隙安装值
Figure GDA0002959309210000103
转轮与底环止漏环间隙安装值
Figure GDA0002959309210000104
以及机组安装时的轴系摆度数据为转轮中心处的轴系轨迹
Figure GDA0002959309210000105
得到:
转轮运行时的上止漏环实时间隙值
Figure GDA0002959309210000106
转轮运行时的下止漏环实时间隙值
Figure GDA0002959309210000111
水轮机止漏环结构类型较多,且止漏间隙一般为0.6~2.5mm之间,现有技术方案中尚无有效监测技术可对运行中转轮与固定件之间的止漏环间隙进行在线监测,如图6,本方案基于流体和结构有限元仿真,并结合机组实时运行数据与振摆监测数据,可实现运行中的转轮止漏环间隙的在线评估,可对转轮碰磨故障进行在线预警和趋势分析,并可为因止漏环间隙不均诱发的机组振动提供定量的科学依据。
而所述数据处理单元对转轮疲劳寿命的评估,如图7,包括转轮快速裂纹安全评估和转轮疲劳裂纹评估;
所述转轮快速裂纹安全评估是转轮结构与机组激励源共振导致的快速裂纹,尤其是叶片出水边卡门涡频率共振,是通过将激振力理论计算频率和激振力电站实测频率与转轮叶片在空气中的固有频率和转轮叶片在水环境中的固有频率进行比较,判断是否共振;
所述转轮疲劳裂纹评估是转轮关键危险点处长期承受动载荷而导致的高周疲劳裂纹,是由转轮固有特性、机组运行工况数据、机组启停信息、仿真计算和/或转轮动应力试验得到的转轮危险点应力-时间荷载曲线,结合转轮材料S-N曲线和历史检修及修复信息,通过雨流荷载计数和疲劳累积损伤理论,对转轮的疲劳寿命进行预测。
对于转轮运行状态的健康评估,因为水轮机转轮是制造厂完全定制化设计和制造的,且机组运行状态健康度影响因素除设计、制造之外,还与电站运行方式密切相关,因此,转轮运行状态健康评估是无法事先完全确定的,必须依赖于多方尤其是电站真机运行历史和实时运行工况的,应根据转轮实际情况确定相应的影响因素及其权重。优选的,转轮运行状态健康评估准则如下:
1.机组运行工况偏离最优工况点的远近程度,权重w1
2.转轮关键部位的应力和疲劳损伤程度,权重w2
3.转轮止漏环处安装间隙值与实时动态间隙值,权重w3
4.机组运行历史,包括机组启停次数、甩负荷次数、过速次数、已投运时间、各历史运行工况的时间统计,权重w4
5.流道压力脉动大小,权重w5
6.转轮汽蚀状态,权重w6
7.转轮泥砂磨损状态,权重w7
8.转轮历史故障和修复情况,权重w8
9.机组实时振动状态,权重w9
定义:
Figure GDA0002959309210000121
0≤wi≤1,此处n=9;wi值越小,转轮运行工况越好。
因此,可定义转轮运行健康状态,并按不同颜色进行显示:
最优状态:0≤W≤0.05,绿色显示,
优秀状态:0.05≤W<0.25
良好状态:0.25≤W<0.4
一般状态:0.4≤W<0.6
监视运行状态:0.6≤W<0.75
恶劣状态:0.75≤W≤1,红色显示。

Claims (7)

1.水轮机转轮实时状态评估和生命周期预测系统,其特征在于:包括数据输入端、数据处理单元和数据输出端;所述数据输入端包括水轮机组实时的运行数据输入端和转轮固有特性信息数据输入端;
所述水轮机组实时的运行数据输入端连接有若干用于检测水轮机组实时的运行数据和水文环境数据的实时数据检测单元,所述转轮固有特性信息数据输入端用于输入对应所述水轮机组实时的运行数据参数下的机组机械状态数据,以及对应所述水文环境数据参数下的机械特性受环境影响变化数据;进一步的,所述水轮机组实时的运行数据输入端连接有用于检测机组运行时的工况类数据、状态类监测数据以及机组过流水质泥沙含量和机组过流水质PH值的实时数据监测模块;其中所述工况类数据包括水头、出力、效率和尾水位数值;所述状态类监测数据包括流道压力、全流道压力脉动、机组振动摆度监测数据、推力以及导轴承瓦温数据;
所述转轮固有特性信息数据输入端输入的机组机械状态数据包括机组机械数据、水轮机转轮运转特性曲线和检修历史数据以及机组安装数据;
所述机组机械数据包括弹性模量、材质成分、屈服强度、抗拉强度和S-N曲线的转轮材质信息,和包括转轮尺寸、不平衡重量、表面粗糙度和残余应力的转轮设计与制造数据,以及电站水质和泥沙对转轮磨损影响的相关系数和理论预测曲线;
所述轮机转轮运转特性曲线包括等效率线、等开度线、吸出高度等值线、等空化系数线、压力脉动等值线、出力限制线、转轮叶片转角等值线、关键位置处的静态应力以及关键位置处的动态应力和变形曲线;
所述检修历史数据包括转轮历次故障和修复信息;
所述机组安装数据包括机组安装时的轴系摆度数据、转轮与顶盖止漏环间隙安装值、转轮与底环止漏环间隙安装值;
所述数据处理单元接收所述数据输入端的输入数据,并将水轮机组实时的运行数据输入端的输入数据与转轮固有特性信息输入端的输入数据进行比对,找出转轮固有特性信息数据中与水轮机组实时的运行数据和机械特性受环境影响变化数据对应的数据,对包括转轮止漏环间隙、转轮疲劳寿命和转轮运行状态进行定期评估;
所述数据输出端包括信息显示单元。
2.如权利要求1所述的水轮机转轮实时状态评估和生命周期预测系统,其特征在于:所述机组过流水质泥沙含量是根据多年每月份统计值确定的。
3.如权利要求1所述的水轮机转轮实时状态评估和生命周期预测系统,其特征在于:所述水轮机组实时的运行数据输入端输入的机组机械状态数据也包括实时运行的机组的机械数据、水轮机转轮运转特性曲线和检修历史数据以及机组安装数据;
且所述水轮机组实时的运行数据输入端输入的等效率线、等开度线、吸出高度等值线、等空化系数线、压力脉动等值线、出力限制线、转轮叶片转角等值线,是将若干转轮运行工况Sijk汇同水轮机模型转轮综合特性曲线以及真机稳定性实验结果数据,经过数据读取和比对计算得到的,其中,设H为机组运行水头,N为机组出力,Z为机组尾水位,等效率线ηij=η(Hi,Nj)、等开度线αij=α(Hi,Nj)、吸出高度等值线Hsij=Hs(Hi,Nj)、等空化系数线σij=σ(Hi,Nj)、5%出力限制线P5%-ij=P5%(Hi,Nj)、压力脉动等值线
Figure FDA0002959309200000021
转轮叶片转角等值线θij=θ(Hi,Nj),其中,i=0,1,2,...,n1,n1为样本数,j=0,1,2,...,n2,n2为样本数,k=0,1,2,...,n3,n3为样本数。
4.如权利要求1所述的水轮机转轮实时状态评估和生命周期预测系统,其特征在于:所述水轮机组实时的运行数据输入端输入的水轮机转轮运转特性曲线中关键位置处的静态应力,是将若干转轮运行工况Sijk经过全流道CFD稳态分析和转轮结构FEM稳态分析得到转轮流道稳态分布和转轮稳态平均应力分布,然后通过对转轮稳态平均应力分布进行计算得到转轮关键位置处在工况Sijk下的平均应力σm-ijk,基于机器学习或深度学习的智能算法对平均应力σm-ijk进行拟合得到关键位置处的静态应力数据库;σm为转轮关键位置处的平均应力,σm=f(H,N,Z),其中H为机组运行水头,N为机组出力,Z为机组尾水位;Sijk=g(Hi,Nj,Zk),i=0,1,2,...,n1,n1为样本数,j=0,1,2,...,n2,n2为样本数,k=0,1,2,...,n3,n3为样本数。
5.如权利要求1所述的水轮机转轮实时状态评估和生命周期预测系统,其特征在于:所述水轮机组实时的运行数据输入端输入的轮机转轮运转特性曲线中关键位置处的动态应力和变形曲线,是将若干转轮运行工况Sijk经过转轮流固耦合FSI瞬态分析获得转轮动应力分布和轮转变形分布情况,然后通过转轮动应力实验修正并计算成果提取得到转轮关键位置处在Sijk工况下的动应力Δσijk和相关位置处的变形rijk,所述相关位置处的变形rijk包括转轮上止漏环处的变形值r1 ijk=F1(Hi,Nj,Zk)和\或下止漏环处的变形值r2 ijk=F2(Hi,Nj,Zk),然后基于机器学习或深度学习的智能算法对平均应力rijk进行拟合得到关键位置处的动态应力数据库;其中,H为机组运行水头,N为机组出力,Z为机组尾水位,其中,i=0,1,2,...,n1,n1为样本数,j=0,1,2,...,n2,n2为样本数,k=0,1,2,...,n3,n3为样本数。
6.如权利要求1所述的水轮机转轮实时状态评估和生命周期预测系统,其特征在于:
所述转轮止漏环间隙的定期评估,是根据所述水轮机组实时的运行数据输入端输入的转轮与顶盖止漏环间隙安装实测值、转轮与底环止漏环间隙安装实测值将若干转轮运行工况Sijk经过流固耦合FSI瞬态分析得到顶盖和底环变形分布信息,再计算得到顶盖止漏环处的变形值R1 ijk=F3(Hi,Nj,Zk)和底环漏环处的变形值R2 ijk=F4(Hi,Nj,Zk);其中,H为机组运行水头,N为机组出力,Z为机组尾水位,其中,i=0,1,2,...,n1,n1为样本数,j=0,1,2,...,n2,n2为样本数,k=0,1,2,...,n3,n3为样本数;
结合所述转轮固有特性信息数据输入端输入的转轮与顶盖止漏环间隙安装值
Figure FDA0002959309200000031
转轮与底环止漏环间隙安装值
Figure FDA0002959309200000032
以及机组安装时的轴系摆度数据为转轮中心处的轴系轨迹
Figure FDA0002959309200000033
得到转轮运行时的上止漏环实时间隙值
Figure FDA0002959309200000034
和转轮运行时的下止漏环实时间隙值
Figure FDA0002959309200000041
7.如权利要求1所述的水轮机转轮实时状态评估和生命周期预测系统,其特征在于:
所述转轮疲劳寿命的定期评估,包括转轮快速裂纹安全评估和转轮疲劳裂纹评估;
所述转轮快速裂纹安全评估是通过将激振力理论计算频率和激振力电站实测频率与转轮叶片在空气中的固有频率和转轮叶片在水环境中的固有频率进行比较判断是否共振;
所述转轮疲劳裂纹评估是通过对转轮疲劳裂纹预测、机组启停工况数据、转轮危险点应力时间频谱进行雨流法循环计算进行预测的。
CN202010135942.6A 2020-03-02 2020-03-02 水轮机转轮实时状态评估和生命周期管理预测系统 Active CN111365158B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010135942.6A CN111365158B (zh) 2020-03-02 2020-03-02 水轮机转轮实时状态评估和生命周期管理预测系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010135942.6A CN111365158B (zh) 2020-03-02 2020-03-02 水轮机转轮实时状态评估和生命周期管理预测系统

Publications (2)

Publication Number Publication Date
CN111365158A CN111365158A (zh) 2020-07-03
CN111365158B true CN111365158B (zh) 2021-06-18

Family

ID=71208457

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010135942.6A Active CN111365158B (zh) 2020-03-02 2020-03-02 水轮机转轮实时状态评估和生命周期管理预测系统

Country Status (1)

Country Link
CN (1) CN111365158B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111859566B (zh) * 2020-07-17 2023-11-17 重庆大学 基于数字孪生的表面粗糙度稳定方法
CN112196717A (zh) * 2020-08-18 2021-01-08 华电电力科学研究院有限公司 一种水轮机调速器液压系统状态监测与故障诊断方法
CN112308300A (zh) * 2020-10-19 2021-02-02 中国电子科技集团公司第三十八研究所 一种大型电子产品状态预测管理系统
CN112417619B (zh) * 2020-11-23 2021-10-08 江苏大学 一种基于数字孪生的泵机组优化运行调节系统及方法
CN112796920B (zh) * 2020-12-14 2022-08-05 华能澜沧江水电股份有限公司 立式混流式水轮发电机转轮贯穿性裂纹早期预警的方法
CN114687907A (zh) * 2020-12-29 2022-07-01 中国电建集团华东勘测设计研究院有限公司 水电站输水发电系统水力安全预警方法
CN113027658B (zh) * 2021-03-24 2022-08-02 华中科技大学 一种水轮机转轮实时状态评估方法及其应用
CN113090440B (zh) * 2021-04-13 2023-03-17 哈动国家水力发电设备工程技术研究中心有限公司 一种检验水泵水轮机转轮破损的方法
CN113124939A (zh) * 2021-05-07 2021-07-16 陆雨轩 一种轮机叶片在线评估系统及轮机叶片损伤预警方法
CN113833593A (zh) * 2021-08-26 2021-12-24 长江勘测规划设计研究有限责任公司 高水头多泥沙条件下立式同轴双级混流式水轮机检修周期量化方法
CN113898516A (zh) * 2021-11-09 2022-01-07 国网信通亿力科技有限责任公司 一种小水电运营样本生产数据现场采样收集方法
CN114280955B (zh) * 2021-11-18 2023-12-19 南方电网调峰调频发电有限公司 基于领域信息模型及数字孪生的智能水电厂体系系统及架构方法
CN114136573B (zh) * 2021-11-29 2024-01-23 湖南五凌电力科技有限公司 水轮发电机组关键部件的预警振幅的计算方法及相关设备
CN114704418B (zh) * 2021-12-22 2024-05-10 黄河水利水电开发集团有限公司 水轮发电机状态监测系统及其调速器协联关系优化方法
CN115471059B (zh) * 2022-09-03 2023-12-01 四川华电泸定水电有限公司 一种基于空蚀振动区规划的水轮机在线经济运行系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3339637A1 (en) * 2016-12-20 2018-06-27 General Electric Company Wind turbine tower placement and orientation
CN109268201A (zh) * 2018-09-06 2019-01-25 湖南江河机电自动化设备股份有限公司 一种水轮机转轮磨损大数据分析系统
CN110206676A (zh) * 2019-05-27 2019-09-06 浙江大学 水轮机远程监控、诊断预测与智能维护调控系统
CN110552832A (zh) * 2019-08-13 2019-12-10 大唐水电科学技术研究院有限公司 水轮机组故障诊断及健康状态评估方法
KR20200019107A (ko) * 2019-10-10 2020-02-21 (주) 플랜디 글로벌 대형 구조물 해체 공법
CN110826203A (zh) * 2019-10-25 2020-02-21 成都大汇智联科技有限公司 水轮机顶盖连接螺栓疲劳寿命的评估方法及监测系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20000019107U (ko) * 1999-04-02 2000-11-06 김종열 안경다리의 탄지장치

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3339637A1 (en) * 2016-12-20 2018-06-27 General Electric Company Wind turbine tower placement and orientation
CN109268201A (zh) * 2018-09-06 2019-01-25 湖南江河机电自动化设备股份有限公司 一种水轮机转轮磨损大数据分析系统
CN110206676A (zh) * 2019-05-27 2019-09-06 浙江大学 水轮机远程监控、诊断预测与智能维护调控系统
CN110552832A (zh) * 2019-08-13 2019-12-10 大唐水电科学技术研究院有限公司 水轮机组故障诊断及健康状态评估方法
KR20200019107A (ko) * 2019-10-10 2020-02-21 (주) 플랜디 글로벌 대형 구조물 해체 공법
CN110826203A (zh) * 2019-10-25 2020-02-21 成都大汇智联科技有限公司 水轮机顶盖连接螺栓疲劳寿命的评估方法及监测系统

Also Published As

Publication number Publication date
CN111365158A (zh) 2020-07-03

Similar Documents

Publication Publication Date Title
CN111365158B (zh) 水轮机转轮实时状态评估和生命周期管理预测系统
US8116990B2 (en) Method and system for real-time prognosis analysis and usage based residual life assessment of turbine engine components and display
CN106404403B (zh) 用于涡轮机的分析的方法和系统
JP5897824B2 (ja) ターボ機械のリスク解析を実施するシステム及び、コンピュータにターボ機械に対する検査勧告の作成を実施させる機械可読命令を格納する機械可読媒体
CN111539553A (zh) 基于svr算法和偏峰度的风电机组发电机故障预警方法
CN112417733B (zh) 超期服役水轮发电机组剩余寿命的综合评估方法
CN105184059A (zh) 一种基于海量数据的水电机组状态分析评估方法
CN110702394A (zh) 一种基于振动变化特征的汽轮发电机组振动故障诊断方法
CN113374652A (zh) 一种风力发电机组寿命评估方法
CN112329357B (zh) 一种清水离心泵振动故障的简易诊断方法及诊断系统
CN113158352B (zh) 基于关键部件的水轮发电机组剩余使用寿命的评估方法
CN110332080B (zh) 一种基于共振响应的风机叶片健康实时监测方法
CN116976191A (zh) 一种水轮发电机组轴系稳定性劣化趋势预测方法
Georgievskaia Analytical system for predicting cracks in hydraulic turbines
CN113446146A (zh) 一种在线水轮机效率试验方法
CN115165326A (zh) 通过机械传动链润滑油(脂)杂质分析的风机故障诊断方法
CN110993132B (zh) 一种核电厂支持疲劳监测功能的瞬态监测方法
CN110441063B (zh) 一种监测诊断大型高速转子轴裂纹的方法
Rademakers et al. Estimating costs of operation & maintenance for offshore wind farms
CN114295367A (zh) 一种风电机组齿轮箱工况在线监测方法
Nilsson Maintenance management of wind power systems: Cost effect analysis of condition monitoring systems
CN115471059B (zh) 一种基于空蚀振动区规划的水轮机在线经济运行系统
CN114611424B (zh) 一种融合声纹数据与cae算法的大型风机叶片寿命预测方法
CN114000989B (zh) 一种风力发电机组叶片气动性能衰减检测方法及测试系统
CN113065670A (zh) 一种发电厂离心式通风机状态检修工作方法

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