CN106050637A - 一种大型变转速给水泵运行状态的在线监测方法 - Google Patents
一种大型变转速给水泵运行状态的在线监测方法 Download PDFInfo
- Publication number
- CN106050637A CN106050637A CN201610386058.3A CN201610386058A CN106050637A CN 106050637 A CN106050637 A CN 106050637A CN 201610386058 A CN201610386058 A CN 201610386058A CN 106050637 A CN106050637 A CN 106050637A
- Authority
- CN
- China
- Prior art keywords
- feed pump
- lift
- speed
- curve
- flow
- 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.)
- Granted
Links
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04B—POSITIVE-DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS
- F04B51/00—Testing machines, pumps, or pumping installations
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04B—POSITIVE-DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS
- F04B49/00—Control, e.g. of pump delivery, or pump pressure of, or safety measures for, machines, pumps, or pumping installations, not otherwise provided for, or of interest apart from, groups F04B1/00 - F04B47/00
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04B—POSITIVE-DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS
- F04B49/00—Control, e.g. of pump delivery, or pump pressure of, or safety measures for, machines, pumps, or pumping installations, not otherwise provided for, or of interest apart from, groups F04B1/00 - F04B47/00
- F04B49/10—Other safety measures
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Control Of Non-Positive-Displacement Pumps (AREA)
- Control Of Positive-Displacement Pumps (AREA)
Abstract
本发明涉及一种大型变转速给水泵运行状态的在线监测方法,包括记录给水泵状态监测的运行参数;根据给水泵设计流量‑扬程曲线,确定给水泵正常运行区域的各边界的数值函数;收集给水泵新投运正常状态的运行数据,经修正后拟合出给水泵某一规定转速下的流量‑扬程性能曲线;根据监测状态下的给水泵运行参数,判断给水泵是否在正常运行区域内工作。如否,给出运行异常状态的报警提示;根据监测状态下的给水泵运行参数,修正至给水泵某一规定转速与正常状态曲线相比,如与正常状态曲线偏差超出一定程度,给出给水泵运行异常状态的报警提示,通过该方法可以对给水泵的日常运行状态进行监测,并给出给水泵运行状态偏离正常运行区域时的判断准则。
Description
技术领域
本发明涉及一种大型变转速给水泵运行状态的在线监测方法。
背景技术
大型给水泵作为电厂中重要的机械转动设备,承担着向锅炉提供高压力、大流量给水的功能。随着电力机组向大容量高参数发展,电厂中给水泵运行的可靠性对机组运行的影响越来越重要。以某一亚临界600MW汽轮机组配置两台给水泵为例,如一台给水泵异常停运,则相应的机组出力需要降低一半。目前大型变转速给水泵运行状态的监测,不仅需要判断给水泵是否运行在正常区域,还需要对给水泵长期的运行状态变化趋势进行监测,目前现场运行的DCS测量系统无法满足上述要求。
因此电厂迫切需要一种大型变转速给水泵运行状态的在线监测方法,来能够实现上述两个方面的要求。
发明内容
本发明的目的在于克服现有技术的缺陷,提供一种大型变转速给水泵运行状态的在线监测方法,通过该方法可以对给水泵的日常运行状态进行监测,并给出给水泵运行状态偏离正常运行区域时的判断准则。
为了实现上述目的,本发明采取的技术方案如下:
一种大型变转速给水泵运行状态的在线监测方法,其关键技术在于,其包括如下步骤:
(1)、记录给水泵状态监测的运行参数,具体包括:给水泵转速、给水泵进口压力、给水泵进口温度、给水泵出口压力、给水泵出口温度和给水泵进口流量;
(2)、根据给水泵设计流量-扬程曲线,确定给水泵正常运行区域的各边界的数值函数。
(3)、收集给水泵新投运正常状态的运行数据,经修正后拟合出给水泵某一规定转速下的流量-扬程性能曲线;
(4)、根据监测状态下的给水泵运行参数,判断给水泵是否在正常运行区域内工作。如否,给出运行异常状态的报警提示;
(5)、根据监测状态下的给水泵运行参数,修正至给水泵某一规定转速与正常状态曲线相比,如与正常状态曲线偏差超出一定程度,给出给水泵运行异常状态的报警提示。
进一步地,所述步骤(2)的具体步骤如下:
a). 额定转速下的设计流量-扬程曲线拟合,首先根据厂家提供的给水泵额定转速下的设计流量-扬程曲线,拟合出数值函数,H=f(Q),f为拟合多项式函数;
b). 最大转速曲线的确定
根据厂家提供的最大转速nmax,列出最大转速下的流量-扬程曲线如下:
H=(nmax/nd)2×f[Q×(nd/nmax)],式中nd为给水泵额定转速,nmax为给水泵最大转速,f为拟合多项式函数,与额定转速下的拟合函数相同;
c). 最小转速曲线的确定
根据厂家提供的最小转速nmin,列出最小转速下的流量-扬程曲线如下:
H=(nmin/nd)2×f[Q×(nd/nmin)],式中nd为给水泵额定转速,nmin为给水泵最小转速,f为拟合多项式函数,与额定转速下的拟合函数相同;
d). 最大流量曲线的确定
按照拟合得到的额定转速下的设计流量-扬程曲线,按Q1=1.2×Qd,计算得到相应的额定转速下的对应扬程H1:
H1=f(Q1)=f(1.2×Qd)
式中Qd为给水泵额定转速下的设计流量,Q1为额定转速下的给水泵最大流量,H1为额定转速下给水泵最大流量对应的扬程,f为拟合多项式函数,与额定转速下的拟合函数相同;
计算出Q1及H1后,最大流量曲线的函数如下:
H=f1(Q)= (H1/Q1 2)×Q2
式中,f1为获得的最大流量曲线的二次函数。
e). 最小流量曲线的确定
按照拟合得到的额定转速下的设计流量-扬程曲线,按Q2=0.25×Qd,计算得到相应的额定转速下的对应扬程H2:
H2=f(Q2)=f(0.25×Qd)
式中Qd为给水泵额定转速下的设计流量,Q2为额定转速下的给水泵最小流量,H2为额定转速下给水泵最小流量对应的扬程,f为拟合多项式函数,与额定转速下的拟合函数相同。
计算出Q2及H2后,最小流量曲线的函数如下:
H=f2(Q)=(H2/Q2 2)×Q2
式中,f2为获得的最小流量曲线的二次函数。
进一步地,所述步骤(3)的具体步骤如下:
a). 数据处理
收集给水泵新投运正常状态的运行数据,将收集的各工况运行参数见表1,分别按照公式(1)和公式(2)得到各工况的流量Q及扬程H;
表1
给水泵状态监测中的流量为容积流量Q,其计算公式为:
(1)
式中,
F为给水泵进口流量,单位t/h;
ρ为经过给水泵的给水密度(kg/m3),计算如下式;
经过给水泵的给水密度由公式ρ=(f1(pi,ti)+f1(po,to))/2计算求得,所述函数f1从国际水和水蒸汽特性协会(IAPWS)的相关计算公式中查得;
给水泵状态监测中的扬程H,其计算公式为:
H=(po-pi)×106/(ρg) (2)
式中:
po,pi为给水泵出口和进口处液体的压力(MPa),为实测量;
ρ为经过给水泵的给水密度(kg/m3),计算如上。
g为重力加速度常数,取9.8,单位N/kg;
b). 各工况数据修正至额定转速
按照下列公式,将各数据n、Q、H修正至额定转速nd:
Qc = Q×(nd/n)
Hc = H×(nd/n)2
c). 修正后流量、扬程数据的拟合
将修正后的各个工况点(Qc,Hc),进行多项式拟合,得到给水泵新投运后正常的实际运行特性曲线:H=f3(Q),f3为获得的给水泵新投运后额定转速下的流量扬程多项式拟合函数。
进一步地,所述步骤(4)的具体步骤如下:
a). 运行区域交点的计算。根据确定的给水泵运行范围四条拟合曲线,计算出交点的流量及扬程数据A(QA,HA),B(QB,HB),C(QC,HC),D(QD,HD);
b). 运行状态点处理。对给水泵每一组运行参数,包括n、pi、ti、po、to、F,按照计算原理部分所述的方法,计算得到给水泵转速n、流量Q及扬程H;
c). 运行异常的判断逻辑。如下所示判断逻辑:
IF Q<QA,THEN 报警提示“运行流量小于给水泵最小要求流量!”
ELSEIF Q<QB,THEN
IF H>=f2(Q) THEN 报警提示“运行流量小于给水泵最小要求流量!”(注:f2(Q)为最小流量曲线的二次函数计算值)
ELSEIF H>=(nmin/nd)2×f[Q×(nd/nmin)] THEN 提示“给水泵运行参数正常” (注:不等式右侧为最小转速流量扬程曲线拟合公式计算值)
ELSE报警提示“运行扬程小于给水泵最小要求扬程!”
END IF
ELSEIF Q<QD,THEN
IF H>=(nmax/nd)2×f[Q×(nd/nmax)] THEN提示“给水泵运行参数异常,扬程超限!”(注:不等式右侧为最大转速流量扬程曲线拟合公式计算值)
ELSEIF H>=(nmin/nd)2×f[Q×(nd/nmin)] THEN 提示“给水泵运行参数正常”(注:不等式右侧为最小转速流量扬程曲线拟合公式计算值)
ELSE报警提示“运行扬程小于给水泵最小要求扬程!”
END IF
ELSEIF Q<QC,THEN
IF H>=(nmax/nd)2×f[Q×(nd/nmax)] THEN提示“给水泵运行参数异常,扬程超限!”(注:不等式右侧为最大转速流量扬程曲线拟合公式计算值)
ELSEIF H>=f1(Q) THEN 提示“给水泵运行参数正常” (注:f1(Q)为最大流量曲线的二次函数计算值)
ELSE报警提示“运行流量超出给水泵最大流量!”
END IF
ELSE 报警提示“运行流量超出给水泵最大流量!”
END IF
进一步地,所述步骤(5)的具体步骤如下:
a). 运行状态点处理
对给水泵每一组运行参数,包括n、pi、ti、po、to、F,按照计算原理部分所述的方法,计算得到给水泵运行状态点参数包括给水泵转速n、流量Q及扬程H;
b). 运行状态点参数修正至额定转速nd:
修正至额定转速nd的状态参数:
Qc = Q×(nd/n)
Hc = H×(nd/n)2
c). 运行异常的判断逻辑。如下所示判断逻辑:
IF Hc<f3(Qc) AND ABS[(f3(Qc)-Hc)/f3(Qc)]>=δTHEN(注:f3(Qc)为根据给水泵新投运后额定转速下的流量扬程多项式拟合函数计算得到的扬程值)
报警提示“给水泵扬程下降超出规定范围值!”
ELSE
提示“给水泵运行参数正常”
END IF
上式中,δ为扬程降低的允许规定值,可设定为0.1。
与现有技术相比,本发明所取得的有益效果如下:
(1)本发明基于常用的给水泵运行参数,采取本发明所提供的判断准则及方法即可对正常运行变转速给水泵的状态进行判定,可以为日常运行中电厂监测给水泵状态、保证给水泵案例运行提供准确的基础数据;采用本专利方法还可以对给水泵的工作状态进行连续的状态监控,对给水泵工作状态的变化及时进行预警,满足了电厂状态检修的需要;
(2)本发明数据采用的测量参数现场容易获取,监测方法中的各种基准曲线(四条运行区域边界曲线,及正常运行状态额定转速下的流量-扬程曲线)容易计算拟合,降低了本专利在现场的实施难度。
本监测方法采用的数据容易测量,结构简单,监测判断方法可靠易行,便于在各个电厂的现场实施。
附图说明
附图1为典型给水泵运行区域示意图;
附图2为实施例1流量-扬程曲线。
具体实施方式
以下结合附图和具体实施例对本发明进行进一步详细的叙述。
本发明是基于以下原理实现的:
1、反映给水泵运行状态的主要运行参数是流量及扬程。正常运行状态下,给水泵的流量与扬程变化应符合一定的变化曲线,即流量-扬程曲线。给水泵运行状态异常时,首先会反映在给水泵流量与扬程的对应关系上。例如,当给水泵流量一定,扬程降低时,表示给水泵的密封泄漏量增加。
2、给水泵设备制造商提供的流量扬程-曲线中,相关的数据现场无法直接测得,必须经过一定的转换。
给水泵现场运行参数包括:给水泵转速、给水泵进口压力、给水泵进口温度、给水泵出口压力、给水泵出口温度、给水泵进口流量F。各参数的相应符号及单位如表1:
表1
参数 | 符号 | 单位 |
给水泵转速 | n | rpm |
给水泵进口压力 | pi | MPa |
给水泵进口温度 | ti | ℃ |
给水泵出口压力 | po | MPa |
给水泵出口温度 | to | ℃ |
给水泵进口流量 | F | t/h |
给水泵容积流量 | Q | m3/h |
给水泵扬程 | H | m |
给水泵状态监测中的流量为容积流量Q,其计算公式为:
(1)
式中,
F为给水泵进口流量,单位t/h;
ρ为经过给水泵的给水密度(kg/m3),计算如下式;
经过给水泵的给水密度由公式ρ=(f1(pi,ti)+f1(po,to))/2计算求得,所述函数f1从国际水和水蒸汽特性协会(IAPWS)的相关计算公式中查得;
给水泵状态监测中的扬程H,其计算公式为:
H=(po-pi)×106/(ρg) (2)
式中:
po,pi为给水泵出口和进口处液体的压力(MPa),为实测量;
ρ为经过给水泵的给水密度(kg/m3),计算如上。
g为重力加速度常数,取9.8,单位N/kg。
3、给水泵的正常运行区域如附图1所示,从图中可以看到,正常运行中给水泵的运行工况点处在四个曲线围的区域内。这四个曲线分别是最大转速曲线、最小转速曲线、最大流量曲线及最小流量曲线。
本发明具体实施的步骤如下:
(1)、记录给水泵状态监测的运行参数见表1,具体包括:给水泵转速、给水泵进口压力、给水泵进口温度、给水泵出口压力、给水泵出口温度和给水泵进口流量。
(2)、根据给水泵设计流量-扬程曲线,确定给水泵正常运行区域的各边界的数值函数。
a). 额定转速下的设计流量-扬程曲线拟合,首先根据厂家提供的给水泵额定转速下的设计流量-扬程曲线,拟合出数值函数,H=f(Q),f为拟合多项式函数;
b). 最大转速曲线的确定
根据厂家提供的最大转速nmax,列出最大转速下的流量-扬程曲线如下:
H=(nmax/nd)2×f[Q×(nd/nmax)],式中nd为给水泵额定转速,nmax为给水泵最大转速,f为拟合多项式函数,与额定转速下的拟合函数相同;
c). 最小转速曲线的确定
根据厂家提供的最小转速nmin,列出最小转速下的流量-扬程曲线如下:
H=(nmin/nd)2×f[Q×(nd/nmin)],式中nd为给水泵额定转速,nmin为给水泵最小转速,f为拟合多项式函数,与额定转速下的拟合函数相同;
d). 最大流量曲线的确定
按照拟合得到的额定转速下的设计流量-扬程曲线,按Q1=1.2×Qd,计算得到相应的额定转速下的对应扬程H1:
H1=f(Q1)=f(1.2×Qd)
式中Qd为给水泵额定转速下的设计流量,Q1为额定转速下的给水泵最大流量,H1为额定转速下给水泵最大流量对应的扬程,f为拟合多项式函数,与额定转速下的拟合函数相同;
计算出Q1及H1后,最大流量曲线的函数如下:
H=f1(Q)= (H1/Q1 2)×Q2
式中,f1为获得的最大流量曲线的二次函数。
e). 最小流量曲线的确定
按照拟合得到的额定转速下的设计流量-扬程曲线,按Q2=0.25×Qd,计算得到相应的额定转速下的对应扬程H2:
H2=f(Q2)=f(0.25×Qd)
式中Qd为给水泵额定转速下的设计流量,Q2为额定转速下的给水泵最小流量,H2为额定转速下给水泵最小流量对应的扬程,f为拟合多项式函数,与额定转速下的拟合函数相同。
计算出Q2及H2后,最小流量曲线的函数如下:
H=f2(Q)=(H2/Q2 2)×Q2
式中,f2为获得的最小流量曲线的二次函数。
(3)、收集给水泵新投运正常状态的运行数据,经修正后拟合出给水泵某一规定转速下的流量-扬程性能曲线;
a). 数据处理
收集给水泵新投运正常状态的运行数据,将收集的各工况运行参数见表1,分别按照公式(1)和公式(2)得到各工况的流量Q及扬程H;
b). 各工况数据修正至额定转速
按照下列公式,将各数据n、Q、H修正至额定转速nd:
Qc = Q×(nd/n)
Hc = H×(nd/n)2
c). 修正后流量、扬程数据的拟合
将修正后的各个工况点(Qc,Hc),进行多项式拟合,得到给水泵新投运后正常的实际运行特性曲线:H=f3(Q),f3为获得的给水泵新投运后额定转速下的流量扬程多项式拟合函数。
(4)、根据监测状态下的给水泵运行参数,判断给水泵是否在正常运行区域内工作。如否,给出运行异常状态的报警提示;
a). 运行区域交点的计算。根据确定的给水泵运行范围四条拟合曲线,计算出交点的流量及扬程数据A(QA,HA),B(QB,HB),C(QC,HC),D(QD,HD);
b). 运行状态点处理。对给水泵每一组运行参数,包括n、pi、ti、po、to、F,按照计算原理部分所述的方法,计算得到给水泵转速n、流量Q及扬程H;
c). 运行异常的判断逻辑。如下所示判断逻辑:
IF Q<QA,THEN 报警提示“运行流量小于给水泵最小要求流量!”
ELSEIF Q<QB,THEN
IF H>=f2(Q) THEN 报警提示“运行流量小于给水泵最小要求流量!”(注:f2(Q)为最小流量曲线的二次函数计算值)
ELSEIF H>=(nmin/nd)2×f[Q×(nd/nmin)] THEN 提示“给水泵运行参数正常” (注:不等式右侧为最小转速流量扬程曲线拟合公式计算值)
ELSE报警提示“运行扬程小于给水泵最小要求扬程!”
END IF
ELSEIF Q<QD,THEN
IF H>=(nmax/nd)2×f[Q×(nd/nmax)] THEN提示“给水泵运行参数异常,扬程超限!”(注:不等式右侧为最大转速流量扬程曲线拟合公式计算值)
ELSEIF H>=(nmin/nd)2×f[Q×(nd/nmin)] THEN 提示“给水泵运行参数正常”(注:不等式右侧为最小转速流量扬程曲线拟合公式计算值)
ELSE报警提示“运行扬程小于给水泵最小要求扬程!”
END IF
ELSEIF Q<QC,THEN
IF H>=(nmax/nd)2×f[Q×(nd/nmax)] THEN提示“给水泵运行参数异常,扬程超限!”(注:不等式右侧为最大转速流量扬程曲线拟合公式计算值)
ELSEIF H>=f1(Q) THEN 提示“给水泵运行参数正常” (注:f1(Q)为最大流量曲线的二次函数计算值)
ELSE报警提示“运行流量超出给水泵最大流量!”
END IF
ELSE 报警提示“运行流量超出给水泵最大流量!”
END IF
(5)、根据监测状态下的给水泵运行参数,修正至给水泵某一规定转速与正常状态曲线相比,如与正常状态曲线偏差超出一定程度,给出给水泵运行异常状态的报警提示。
a). 运行状态点处理
对给水泵每一组运行参数,包括n、pi、ti、po、to、F,按照计算原理部分所述的方法,计算得到给水泵运行状态点参数包括给水泵转速n、流量Q及扬程H;
b). 运行状态点参数修正至额定转速nd:
修正至额定转速nd的状态参数:
Qc = Q×(nd/n)
Hc = H×(nd/n)2
c). 运行异常的判断逻辑。如下所示判断逻辑:
IF Hc<f3(Qc) AND ABS[(f3(Qc)-Hc)/f3(Qc)]>=δTHEN(注:f3(Qc)为根据给水泵新投运后额定转速下的流量扬程多项式拟合函数计算得到的扬程值)
报警提示“给水泵扬程下降超出规定范围值!”
ELSE
提示“给水泵运行参数正常”
END IF
上式中,δ为扬程降低的允许规定值,可设定为0.1。
实施例1
步骤(1)、记录给水泵状态监测的运行参数见表1,具体包括:给水泵转速、给水泵进口压力、给水泵进口温度、给水泵出口压力、给水泵出口温度和给水泵进口流量。计算相应的流量及扬程参数。
给水泵在某稳定状态下的参数如下表。
参数 | 符号 | 单位 | 数据 |
给水泵转速 | n | rpm | 4453 |
给水泵进口压力 | pi | MPa | 1.942 |
给水泵进口温度 | ti | ℃ | 164.52 |
给水泵出口压力 | po | MPa | 18.386 |
给水泵出口温度 | to | ℃ | 167.85 |
给水泵进口流量 | F | t/h | 975.34 |
根据上述运行参数,可以计算得到此转速下的给水泵流量及扬程为:
Q=1075.373 m3/h
H=1848.168 m
步骤(2)、根据给水泵设计流量-扬程曲线,确定给水泵正常运行区域的各边界的数值函数。
某机组给水泵的额定转速为4956 rpm,额定流量1318m3/h,额定扬程2145m。其设计流量-扬程曲线如附图2所示。
根据图中曲线,可以拟合得到额定转速4956 rpm下的如下多项式函数。
H=a+bQ+cQ2+dQ3
其中,多项式的系数分别为:
a = 3095.4545
b = -0.21454286
c = -0.00021393421
d = -6.1584506e-008
此给水泵的最大转速为5054 rpm,此时,给水泵的流量扬程曲线函数为:
(nd/nmax)=4956/5054=0.9806;(nmax/nd)=5054/4956=1.0198
H=[a+b(Q×0.9806)+c(Q×0.9806)2+d(Q×0.9806)3] ×1.04
其中,多项式的系数分别为:
a = 3095.4545
b = -0.21454286
c = -0.00021393421
d = -6.1584506e-008
此给水泵的最小转速为1517 rpm,此时,给水泵的流量扬程曲线函数为:
(nd/nmin)=4956/1517=3.267;(nmin/nd)=1517/4956=0.306
H=[a+b(Q×3.267)+c(Q×3.267)2+d(Q×3.267)3] ×0.0937
其中,多项式的系数分别为:
a = 3095.4545
b = -0.21454286
c = -0.00021393421
d = -6.1584506e-008
给水泵运行最大流量为1.2倍额定流量,即1581.6m3/h。此时按照额定转速下的流量扬程拟合公式,可以计算得到对应的扬程为1978.15m。因此最大流量曲线的函数为:
H=(1978.15/1581.62)Q2=7.914×10-4Q2
给水泵运行最小流量为0.25倍额定流量,即329.5m3/h。此时按照额定转速下的流量扬程拟合公式,可以计算得到对应的扬程为2999.33m。因此最小流量曲线的函数为:
H=(2999.33/329.52)Q2=2.763×10-2Q2
步骤(3)、收集给水泵新投运正常状态的运行数据,经修正后拟合出给水泵某一规定转速下的流量-扬程性能曲线;
将收集数据计算并修正至额定转速下后得到的流量、扬程数据如下表所示:
流量 | 扬程 |
450 | 2814 |
600 | 2751 |
750 | 2658 |
900 | 2575 |
1050 | 2442 |
1200 | 2339 |
1350 | 2161 |
1500 | 2003 |
拟合得到的多项式公式为:
H=a+bQ+cQ2+dQ3
其中,多项式的系数分别为:
a = 2938.45
b = -0.16788
c = -0.00021393
d = -6.15845e-008
步骤(4)、根据监测状态下的给水泵运行参数,判断给水泵是否在正常运行区域内工作。如否,给出运行异常状态的报警提示;
根据步骤(2)中确定的给水泵四条运行边界曲线,可以确定出相应的四个交点(交点横坐标为流量,m3/h;纵坐标为扬程,m)为:
交点A参数,(100.85,281.03)
交点B参数,(336,3119.3)
交点C参数,(1612.3,2057.25)
交点D参数,(484,185.39)
运行异常判断:
步骤(1)获得的给水泵正常运行流量为1075.373 m3/h,且有QD<1075.373<QC,在此状态下:
对应1075.373m3/h流量时的最大转速流量扬程曲线拟合公式计算得到的扬程为2661.473m,对应1075.373m3/h流量时的最大流量二次曲线函数计算得到的扬程为915.196m,因为915.196<1848.168<2661.473,因此状态判断结果为“给水泵运行参数正常”。
步骤(5)、根据监测状态下的给水泵运行参数,修正至给水泵某一规定转速与正常状态曲线相比,如与正常状态曲线偏差超出一定程度,给出给水泵运行异常状态的报警提示。
根据步骤(1)计算的流量及扬程,修正至额定转速4956rpm时的给水泵流量及扬程为:
Q=1196.844 m3/h
H=2289.278 m
根据步骤(3)获得的给水泵新投运后额定转速下的流量扬程多项式拟合函数,计算得到对应1196.844 m3/h流量下的扬程值为2325.496m。
由于2289.278<2325.496,且(2325.496-2289.278)/2325.496=0.0156<0.1,因此状态判断结果为“给水泵运行参数正常”。
以上所述实施方式仅为本发明的优选实施例,而并非本发明可行实施的穷举。对于本领域一般技术人员而言,在不背离本发明原理和精神的前提下对其所作出的任何显而易见的改动,都应当被认为包含在本发明的权利要求保护范围之内。
Claims (5)
1.一种大型变转速给水泵运行状态的在线监测方法,其特征在于,其包括如下步骤:
(1)、记录给水泵状态监测的运行参数,具体包括:给水泵转速、给水泵进口压力、给水泵进口温度、给水泵出口压力、给水泵出口温度和给水泵进口流量;
(2)、根据给水泵设计流量-扬程曲线,确定给水泵正常运行区域的各边界的数值函数;
(3)、收集给水泵新投运正常状态的运行数据,经修正后拟合出给水泵某一规定转速下的流量-扬程性能曲线;
(4)、根据监测状态下的给水泵运行参数,判断给水泵是否在正常运行区域内工作;
如否,给出运行异常状态的报警提示;
(5)、根据监测状态下的给水泵运行参数,修正至给水泵某一规定转速与正常状态曲线相比,如与正常状态曲线偏差超出一定程度,给出给水泵运行异常状态的报警提示。
2.根据权利要求1所述的一种大型变转速给水泵运行状态的在线监测方法,其特征在于,所述步骤(2)的具体步骤如下:
a). 额定转速下的设计流量-扬程曲线拟合,首先根据厂家提供的给水泵额定转速下的设计流量-扬程曲线,拟合出数值函数,H=f(Q),f为拟合多项式函数;
b). 最大转速曲线的确定
根据厂家提供的最大转速nmax,列出最大转速下的流量-扬程曲线如下:
H=(nmax/nd)2×f[Q×(nd/nmax)],式中nd为给水泵额定转速,nmax为给水泵最大转速,f为拟合多项式函数,与额定转速下的拟合函数相同;
c). 最小转速曲线的确定
根据厂家提供的最小转速nmin,列出最小转速下的流量-扬程曲线如下:
H=(nmin/nd)2×f[Q×(nd/nmin)],式中nd为给水泵额定转速,nmin为给水泵最小转速,f为拟合多项式函数,与额定转速下的拟合函数相同;
d). 最大流量曲线的确定
按照拟合得到的额定转速下的设计流量-扬程曲线,按Q1=1.2×Qd,计算得到相应的额定转速下的对应扬程H1:
H1=f(Q1)=f(1.2×Qd)
式中Qd为给水泵额定转速下的设计流量,Q1为额定转速下的给水泵最大流量,H1为额定转速下给水泵最大流量对应的扬程,f为拟合多项式函数,与额定转速下的拟合函数相同;
计算出Q1及H1后,最大流量曲线的函数如下:
H=f1(Q)= (H1/Q1 2)×Q2
式中,f1为获得的最大流量曲线的二次函数;
e). 最小流量曲线的确定
按照拟合得到的额定转速下的设计流量-扬程曲线,按Q2=0.25×Qd,计算得到相应的额定转速下的对应扬程H2:
H2=f(Q2)=f(0.25×Qd)
式中Qd为给水泵额定转速下的设计流量,Q2为额定转速下的给水泵最小流量,H2为额定转速下给水泵最小流量对应的扬程,f为拟合多项式函数,与额定转速下的拟合函数相同;
计算出Q2及H2后,最小流量曲线的函数如下:
H=f2(Q)=(H2/Q2 2)×Q2
式中,f2为获得的最小流量曲线的二次函数。
3.根据权利要求1所述的一种大型变转速给水泵运行状态的在线监测方法,其特征在于,所述步骤(3)的具体步骤如下:
a). 数据处理
收集给水泵新投运正常状态的运行数据,将收集的各工况运行参数见表1,分别按照公式(1)和公式(2)得到各工况的流量Q及扬程H;
表1
给水泵状态监测中的流量为容积流量Q,其计算公式为:
(1)
式中,
F为给水泵进口流量,单位t/h;
ρ为经过给水泵的给水密度(kg/m3),计算如下式;
经过给水泵的给水密度由公式ρ=(f1(pi,ti)+f1(po,to))/2计算求得,所述函数f1从国际水和水蒸汽特性协会(IAPWS)的相关计算公式中查得;
给水泵状态监测中的扬程H,其计算公式为:
H=(po-pi)×106/(ρg) (2)
式中:
po,pi为给水泵出口和进口处液体的压力(MPa),为实测量;
ρ为经过给水泵的给水密度(kg/m3),计算如上;
g为重力加速度常数,取9.8,单位N/kg;
b). 各工况数据修正至额定转速
按照下列公式,将各数据n、Q、H修正至额定转速nd:
Qc = Q×(nd/n)
Hc = H×(nd/n)2
c). 修正后流量、扬程数据的拟合
将修正后的各个工况点(Qc,Hc),进行多项式拟合,得到给水泵新投运后正常的实际运行特性曲线:H=f3(Q),f3为获得的给水泵新投运后额定转速下的流量扬程多项式拟合函数。
4.根据权利要求1所述的一种大型变转速给水泵运行状态的在线监测方法,其特征在于,所述步骤(4)的具体步骤如下:
a). 运行区域交点的计算;
根据确定的给水泵运行范围四条拟合曲线,计算出交点的流量及扬程数据A(QA,HA),B(QB,HB),C(QC,HC),D(QD,HD);
b). 运行状态点处理;
对给水泵每一组运行参数,包括n、pi、ti、po、to、F,按照计算原理部分所述的方法,计算得到给水泵转速n、流量Q及扬程H;
c). 运行异常的判断逻辑;
如下所示判断逻辑:
IF Q<QA,THEN 报警提示“运行流量小于给水泵最小要求流量!”
ELSEIF Q<QB,THEN
IF H>=f2(Q) THEN 报警提示“运行流量小于给水泵最小要求流量!”
ELSEIF H>=(nmin/nd)2×f[Q×(nd/nmin)] THEN 提示“给水泵运行参数正常”
ELSE报警提示“运行扬程小于给水泵最小要求扬程!”
END IF
ELSEIF Q<QD,THEN
IF H>=(nmax/nd)2×f[Q×(nd/nmax)] THEN提示“给水泵运行参数异常,扬程超限!”
ELSEIF H>=(nmin/nd)2×f[Q×(nd/nmin)] THEN 提示“给水泵运行参数正常”
ELSE报警提示“运行扬程小于给水泵最小要求扬程!”
END IF
ELSEIF Q<QC,THEN
IF H>=(nmax/nd)2×f[Q×(nd/nmax)] THEN提示“给水泵运行参数异常,扬程超限!”
ELSEIF H>=f1(Q) THEN 提示“给水泵运行参数正常”
ELSE报警提示“运行流量超出给水泵最大流量!”
END IF
ELSE 报警提示“运行流量超出给水泵最大流量!”
END IF。
5.根据权利要求1所述的一种大型变转速给水泵运行状态的在线监测方法,其特征在于,所述步骤(5)的具体步骤如下:
a). 运行状态点处理
对给水泵每一组运行参数,包括n、pi、ti、po、to、F,按照计算原理部分所述的方法,计算得到给水泵运行状态点参数包括给水泵转速n、流量Q及扬程H;
b). 运行状态点参数修正至额定转速nd:
修正至额定转速nd的状态参数:
Qc = Q×(nd/n)
Hc = H×(nd/n)2
c). 运行异常的判断逻辑;
如下所示判断逻辑:
IF Hc<f3(Qc) AND ABS[(f3(Qc)-Hc)/f3(Qc)]>=δTHEN
报警提示“给水泵扬程下降超出规定范围值!”
ELSE
提示“给水泵运行参数正常”
END IF
上式中,δ为扬程降低的允许规定值,可设定为0.1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610386058.3A CN106050637B (zh) | 2016-06-03 | 2016-06-03 | 一种大型变转速给水泵运行状态的在线监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610386058.3A CN106050637B (zh) | 2016-06-03 | 2016-06-03 | 一种大型变转速给水泵运行状态的在线监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106050637A true CN106050637A (zh) | 2016-10-26 |
CN106050637B CN106050637B (zh) | 2018-09-21 |
Family
ID=57172155
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610386058.3A Active CN106050637B (zh) | 2016-06-03 | 2016-06-03 | 一种大型变转速给水泵运行状态的在线监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106050637B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108678936A (zh) * | 2018-04-10 | 2018-10-19 | 新界泵业(杭州)有限公司 | 泵运行管控方法 |
CN109578262A (zh) * | 2018-12-13 | 2019-04-05 | 保定申辰泵业有限公司 | 一种蠕动泵传输粘性液体的控制方法、装置及蠕动泵 |
CN111368374A (zh) * | 2020-03-03 | 2020-07-03 | 上海辛格林纳新时达电机有限公司 | 水泵参数的计算方法、系统、服务器、存储介质 |
CN113719446A (zh) * | 2021-08-31 | 2021-11-30 | 华能国际电力股份有限公司上海石洞口第一电厂 | 一种基于数据挖掘的汽动给水泵状态监测系统 |
CN113803283A (zh) * | 2020-06-16 | 2021-12-17 | 中国石油天然气股份有限公司 | 压缩机性能曲线的校正方法与系统 |
CN114251852A (zh) * | 2021-08-24 | 2022-03-29 | 佛山市顺德区美的饮水机制造有限公司 | 即热装置及其控制方法和控制装置、水处理装置和介质 |
CN115496628A (zh) * | 2022-10-21 | 2022-12-20 | 普杰软件(南京)有限公司 | 一种智慧酒店管理系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS57101911A (en) * | 1980-12-18 | 1982-06-24 | Toshiba Corp | Feed water control device |
JPS57163194A (en) * | 1981-03-31 | 1982-10-07 | Toshiba Corp | Speed increasing apparatus |
CN102536774A (zh) * | 2011-12-27 | 2012-07-04 | 酒泉钢铁(集团)有限责任公司 | 工业循环水泵组运行特性及运行能效分析装置 |
CN102562561A (zh) * | 2011-12-27 | 2012-07-11 | 酒泉钢铁(集团)有限责任公司 | 工业循环水系统的泵组运行能效分析方法 |
CN102734148A (zh) * | 2012-07-04 | 2012-10-17 | 杭州哲达科技股份有限公司 | 水泵型式检验方法 |
CN103195698A (zh) * | 2013-04-19 | 2013-07-10 | 哈尔滨工业大学 | 变水位取水泵站中水泵同步变速调节的节能控制方法 |
CN103575339A (zh) * | 2013-10-31 | 2014-02-12 | 无锡溥汇机械科技有限公司 | 流量测量方法及流量控制方法 |
CN105571844A (zh) * | 2015-12-15 | 2016-05-11 | 北京建筑大学 | 一种测量暖通空调系统水泵或阀门流量的方法 |
-
2016
- 2016-06-03 CN CN201610386058.3A patent/CN106050637B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS57101911A (en) * | 1980-12-18 | 1982-06-24 | Toshiba Corp | Feed water control device |
JPS57163194A (en) * | 1981-03-31 | 1982-10-07 | Toshiba Corp | Speed increasing apparatus |
CN102536774A (zh) * | 2011-12-27 | 2012-07-04 | 酒泉钢铁(集团)有限责任公司 | 工业循环水泵组运行特性及运行能效分析装置 |
CN102562561A (zh) * | 2011-12-27 | 2012-07-11 | 酒泉钢铁(集团)有限责任公司 | 工业循环水系统的泵组运行能效分析方法 |
CN102734148A (zh) * | 2012-07-04 | 2012-10-17 | 杭州哲达科技股份有限公司 | 水泵型式检验方法 |
CN103195698A (zh) * | 2013-04-19 | 2013-07-10 | 哈尔滨工业大学 | 变水位取水泵站中水泵同步变速调节的节能控制方法 |
CN103575339A (zh) * | 2013-10-31 | 2014-02-12 | 无锡溥汇机械科技有限公司 | 流量测量方法及流量控制方法 |
CN105571844A (zh) * | 2015-12-15 | 2016-05-11 | 北京建筑大学 | 一种测量暖通空调系统水泵或阀门流量的方法 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108678936A (zh) * | 2018-04-10 | 2018-10-19 | 新界泵业(杭州)有限公司 | 泵运行管控方法 |
CN109578262A (zh) * | 2018-12-13 | 2019-04-05 | 保定申辰泵业有限公司 | 一种蠕动泵传输粘性液体的控制方法、装置及蠕动泵 |
CN111368374A (zh) * | 2020-03-03 | 2020-07-03 | 上海辛格林纳新时达电机有限公司 | 水泵参数的计算方法、系统、服务器、存储介质 |
CN111368374B (zh) * | 2020-03-03 | 2023-04-14 | 上海辛格林纳新时达电机有限公司 | 水泵参数的计算方法、系统、服务器、存储介质 |
CN113803283A (zh) * | 2020-06-16 | 2021-12-17 | 中国石油天然气股份有限公司 | 压缩机性能曲线的校正方法与系统 |
CN113803283B (zh) * | 2020-06-16 | 2024-03-26 | 中国石油天然气股份有限公司 | 压缩机性能曲线的校正方法与系统 |
CN114251852A (zh) * | 2021-08-24 | 2022-03-29 | 佛山市顺德区美的饮水机制造有限公司 | 即热装置及其控制方法和控制装置、水处理装置和介质 |
CN113719446A (zh) * | 2021-08-31 | 2021-11-30 | 华能国际电力股份有限公司上海石洞口第一电厂 | 一种基于数据挖掘的汽动给水泵状态监测系统 |
CN115496628A (zh) * | 2022-10-21 | 2022-12-20 | 普杰软件(南京)有限公司 | 一种智慧酒店管理系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106050637B (zh) | 2018-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106050637A (zh) | 一种大型变转速给水泵运行状态的在线监测方法 | |
KR100975845B1 (ko) | 가스 기구 판별 장치와 그 방법 | |
CN106524613A (zh) | 一种变频风冷热泵机组及其控制方法、控制装置 | |
Schiavello et al. | Pump Cavitation: various NPSHR criteria, NPSHA margins, impeller life expectancy | |
CN104633458B (zh) | 一种火电站水泵实时在线汽蚀监测预警系统及方法 | |
CA2184882A1 (en) | Turbomachinery with variable-angle flow guiding vanes | |
CN108360608B (zh) | 一种供水系统输水管的爆管识别与定位方法 | |
CN105699062A (zh) | 一种阀门流量流阻测试系统及其进行小微流量测试的方法 | |
CN104597934B (zh) | 一种变压器用智能型冷却器控制系统及其控制方法 | |
CN106050722A (zh) | 基于相似原理的通用特性曲线喘振控制方法和系统 | |
CN112329357B (zh) | 一种清水离心泵振动故障的简易诊断方法及诊断系统 | |
CN103742425A (zh) | 水循环系统节能修正方法 | |
CN110056530A (zh) | 防喘振控制方法及系统 | |
CN109723609A (zh) | 一种风电机组变桨系统的故障预警方法及系统 | |
CN103850872A (zh) | 转桨式水轮机协联关系综合优化试验方法 | |
CN105201807A (zh) | 一种基于压力差和流量控制的压缩机运行控制方法和装置 | |
CN105260580A (zh) | 一种确定混流式水轮机尾水管压力脉动幅值的方法 | |
CN108870642A (zh) | 新风净化器滤网状态自校验方法 | |
CN104633457B (zh) | 一种火电站调节阀实时在线汽蚀监测预警系统及方法 | |
CN105865220B (zh) | 一种双背压凝汽器优化运行装置的运行方法 | |
CN104965496B (zh) | 基于防喘振控制系统的防喘振控制方法 | |
CN117076959A (zh) | 大数据可视化交互系统 | |
CN104678992B (zh) | 矿用主通风机稳定性与经济性实时检测方法 | |
CN116928112A (zh) | 一种潜水离心泵的故障监测系统和方法 | |
CN105604953A (zh) | 一种基于振动优化的多级泵叶轮错位布置方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |