CN106126822A - 一种直流锅炉的汽水分离器建模方法及系统 - Google Patents

一种直流锅炉的汽水分离器建模方法及系统 Download PDF

Info

Publication number
CN106126822A
CN106126822A CN201610479586.3A CN201610479586A CN106126822A CN 106126822 A CN106126822 A CN 106126822A CN 201610479586 A CN201610479586 A CN 201610479586A CN 106126822 A CN106126822 A CN 106126822A
Authority
CN
China
Prior art keywords
steam
parameter
working medium
water
water separator
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
Application number
CN201610479586.3A
Other languages
English (en)
Other versions
CN106126822B (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.)
Electric Power Research Institute of Guangdong Power Grid Co Ltd
Original Assignee
Unigroup Beijing Intelligent Control Science & Technology Co ltd
Electric Power Research Institute of Guangdong Power Grid Co Ltd
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 Unigroup Beijing Intelligent Control Science & Technology Co ltd, Electric Power Research Institute of Guangdong Power Grid Co Ltd filed Critical Unigroup Beijing Intelligent Control Science & Technology Co ltd
Priority to CN201610479586.3A priority Critical patent/CN106126822B/zh
Publication of CN106126822A publication Critical patent/CN106126822A/zh
Application granted granted Critical
Publication of CN106126822B publication Critical patent/CN106126822B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本申请公开了一种直流锅炉的汽水分离器建模方法及系统,该方法包括:基于MUSE仿真环境,对汽水分离器进行仿真模型搭建,得到相应的汽水分离器仿真模型;依据质量守恒定律和能量守恒定律,并结合工质物理特性,对汽水分离器仿真模型进行数学建模,得到汽水分离器的机理模型;当汽水分离器仿真模型处于运行仿真过程,则实时分析汽水分离器仿真模型当前所处的仿真工况;在当前参数计算周期内,利用机理模型,对汽水分离器仿真模型当前所处的仿真工况进行参数计算,得到当前参数计算周期下的工质参数;上述仿真工况包括上水启动状态、湿态运行状态、干态运行状态、干湿态转换状态和超临界状态。本申请实现了对直流锅炉中汽水分离器的仿真建模。

Description

一种直流锅炉的汽水分离器建模方法及系统
技术领域
本发明涉及直流锅炉仿真技术领域,特别涉及一种直流锅炉的汽水分离器建模方法及系统。
背景技术
为了提高火力发电机组的效率和减少污染物排放,近年来许多电力企业大力发展了超临界及超超临界发电机组,产生了明显的经济和社会效益。
直流锅炉作为超临界及超超临界发电机组中的重要组成部件,对整个发电机组的性能产生了至关重要的影响。通过对直流锅炉进行仿真建模,可以有利于人们对直流锅炉的设计、制造、调试和运行过程进行后续的优化,并有利于提高对相关运维人员进行业务培训的效果。
而汽水分离器作为直流锅炉中的一个必不可少的组成部件,汽水分离器的仿真建模是直流锅炉仿真建模过程中的一个重要子过程。如何实现对直流锅炉中汽水分离器的仿真建模是目前有待解决的问题。
发明内容
有鉴于此,本发明的目的在于提供一种直流锅炉的汽水分离器建模方法及系统,实现了对直流锅炉中汽水分离器的仿真建模。其具体方案如下:
一种直流锅炉的汽水分离器建模方法,包括:
基于MUSE仿真环境,对汽水分离器进行仿真模型搭建,得到相应的汽水分离器仿真模型;
依据质量守恒定律和能量守恒定律,并结合工质物理特性,对所述汽水分离器仿真模型进行数学建模,得到所述汽水分离器的机理模型;
当所述汽水分离器仿真模型处于运行仿真过程,则实时分析所述汽水分离器仿真模型当前所处的仿真工况;
在当前参数计算周期内,利用所述机理模型,对所述汽水分离器仿真模型当前所处的仿真工况进行参数计算,得到当前参数计算周期下与该仿真工况对应的工质参数;
其中,所述汽水分离器仿真模型的仿真工况包括上水启动状态、湿态运行状态、干态运行状态、干湿态转换状态和超临界状态。
优选的,所述基于MUSE仿真环境,对汽水分离器进行仿真模型搭建的过程,包括:
利用所述MUSE仿真环境,分别对所述汽水分离器的容器主体、入口、出汽口、循环水出口、排污口和安全阀门进行仿真建模,得到所述汽水分离器仿真模型。
优选的,所述依据质量守恒定律和能量守恒定律,并结合工质物理特性,对所述汽水分离器仿真模型进行数学建模的过程,包括:
利用工质质量守恒方程、工质能量守恒方程、水位下工质质量方程和体积守恒方程,对所述汽水分离器仿真模型进行数学建模,得到所述机理模型。
优选的,所述工质质量守恒方程为:
d d τ [ ρ s V s + ρ w V w ] = q r - q s - q a f - q c - q n ;
所述工质能量守恒方程为:
d d τ [ ρ s h s V s + ρ w h w V w - pV t + m t C p t m ] = q r h r - q s h s - q a f h a f - q c h c - q n h n ;
所述水位下工质质量方程为:
d d τ [ ρ s V s d + ρ w V w d ] = x · q r · ( 1 - η ) + ( 1 - x ) · q r - q s d - q c ;
所述体积守恒方程为:
Vt=Vs+Vw
式中,下标s表示蒸汽,下标w表示水,下标r表示入口,下标c表示出水口,下标n表示排污口,下标af表示安全阀门;符号ρ表示工质密度,符号V表示工质体积,符号q表示工质流量,符号h表示工质焓;参数p表示压力参数,参数mt表示分离器金属质量,参数Cp表示金属比热,参数tm表示金属温度,参数Vsd表示水位下蒸汽体积,参数Vwd表示水位下水体积,参数qsd表示由温度过高引起的从水中溢出的饱和蒸汽流量,参数x表示入口工质干度,参数η表示与汽水分离器分离效率有关的预设系数,参数Vt表示汽水分离器的整体内部容积。
优选的,所述在当前参数计算周期内,利用所述机理模型,对所述汽水分离器仿真模型当前所处的仿真工况进行参数计算的过程,包括:
对所述工质质量守恒方程、所述工质能量守恒方程、所述水位下工质质量方程和所述体积守恒方程进行化简代换处理,得到压力变化率方程、水体积变化率方程和水位下蒸汽体积变化率方程;
利用所述压力变化率方程、所述水体积变化率方程和所述水位下蒸汽体积变化率方程,分别计算得到当前参数计算周期内的压力变化率、水体积变化率和水位下蒸汽体积变化率;
利用当前参数计算周期内的压力变化率、水体积变化率和水位下蒸汽体积变化率,并结合上一参数计算周期内得到的压力参数、水体积参数和水位下蒸汽体积参数,计算得到当前参数计算周期内的压力参数、水体积参数和水位下蒸汽体积参数。
优选的,所述汽水分离器建模方法,还包括:
利用当前参数计算周期内的水体积参数和水位下蒸汽体积参数,计算当前参数计算周期内的水位参数。
优选的,所述汽水分离器建模方法,还包括:
若当前参数计算周期内得到的压力参数为饱和压力,则直接利用该压力参数对工质特性参数表进行查表处理,得到相应的工质焓参数和工质密度参数;
若当前参数计算周期内得到的压力参数为过热压力,则利用该压力参数以及当前参数计算周期内采集到的温度参数对所述工质特性参数表进行查表处理,得到相应的工质焓参数和工质密度参数;
其中,所述工质特性参数表包括工程热力学中水蒸气特性参数表和水特性参数表。
优选的,在对所述工质质量守恒方程、所述工质能量守恒方程、所述水位下工质质量方程和所述体积守恒方程进行化简代换处理的过程之前,还包括:
从所述工质质量守恒方程、所述工质能量守恒方程、所述水位下工质质量方程和所述体积守恒方程中确定出在所述汽水分离器仿真模型当前所处的仿真工况下适合被简化的优化项;
利用所述优化项,对所述工质质量守恒方程、所述工质能量守恒方程、所述水位下工质质量方程和所述体积守恒方程中相应的方程进行简化处理。
本发明还公开了一种直流锅炉的汽水分离器建模系统,包括:
仿真模型搭建模块,用于基于MUSE仿真环境,对汽水分离器进行仿真模型搭建,得到相应的汽水分离器仿真模型;
机理模型构建模块,用于依据质量守恒定律和能量守恒定律,并结合工质物理特性,对所述汽水分离器仿真模型进行数学建模,得到所述汽水分离器的机理模型;
仿真工况分析模块,用于当所述汽水分离器仿真模型处于运行仿真过程,则实时分析所述汽水分离器仿真模型当前所处的仿真工况;
工质参数计算模块,用于在当前参数计算周期内,利用所述机理模型,对所述汽水分离器仿真模型当前所处的仿真工况进行参数计算,得到当前参数计算周期下与该仿真工况对应的工质参数;
其中,所述汽水分离器仿真模型的仿真工况包括上水启动状态、湿态运行状态、干态运行状态、干湿态转换状态和超临界状态。
本发明中,汽水分离器建模方法包括:基于MUSE仿真环境,对汽水分离器进行仿真模型搭建,得到相应的汽水分离器仿真模型;依据质量守恒定律和能量守恒定律,并结合工质物理特性,对汽水分离器仿真模型进行数学建模,得到汽水分离器的机理模型;当汽水分离器仿真模型处于运行仿真过程,则实时分析汽水分离器仿真模型当前所处的仿真工况;在当前参数计算周期内,利用机理模型,对汽水分离器仿真模型当前所处的仿真工况进行参数计算,得到当前参数计算周期下与该仿真工况对应的工质参数;其中,汽水分离器仿真模型的仿真工况包括上水启动状态、湿态运行状态、干态运行状态、干湿态转换状态和超临界状态。可见,本发明首先基于MUSE仿真环境,构建了汽水分离器仿真模型,然后对该仿真模型进行数学建模,得到汽水分离器的机理模型;利用该机理模型,可计算得到与上述汽水分离器仿真模型当前所处的仿真工况所对应的工质参数,其中,汽水分离器仿真模型的仿真工况包括上水启动状态、湿态运行状态、干态运行状态、干湿态转换状态和超临界状态,这几种工况能够全面地涵盖了汽水分离器实际工况的种类。由上可见,本发明实现了对直流锅炉中汽水分离器的仿真建模。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本发明实施例公开的一种直流锅炉的汽水分离器建模方法流程图;
图2为本发明实施例公开的一种直流锅炉的汽水分离器建模方法子流程图;
图3为本发明实施例公开的一种直流锅炉的汽水分离器建模系统结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例公开了一种直流锅炉的汽水分离器建模方法,参见图1所示,该方法包括:
步骤S11:基于MUSE仿真环境,对汽水分离器进行仿真模型搭建,得到相应的汽水分离器仿真模型。
需要说明的是,上述MUSE仿真环境具体为Multi-User Simulation Environment,也被称为一体化仿真支撑系统。
上述步骤S11的过程具体可以包括,利用MUSE仿真环境,分别对汽水分离器的容器主体、入口、出汽口、循环水出口、排污口和安全阀门进行仿真建模,得到汽水分离器仿真模型。
步骤S12:依据质量守恒定律和能量守恒定律,并结合工质物理特性,对汽水分离器仿真模型进行数学建模,得到汽水分离器的机理模型。
也即,本实施例在对汽水分离器仿真模型进行数学建模时,不仅结合了质量和能量守恒定律,还需结合工质本身的物理特性来展开建模,由此使得最终得到的机理模型能够更加贴合汽水分离器的实际情况。
步骤S13:当汽水分离器仿真模型处于运行仿真过程,则实时分析汽水分离器仿真模型当前所处的仿真工况。
步骤S14:在当前参数计算周期内,利用机理模型,对汽水分离器仿真模型当前所处的仿真工况进行参数计算,得到当前参数计算周期下与该仿真工况对应的工质参数。
可以理解的是,本实施例是基于上一参数计算周期内得到的工质参数来分析当前汽水分离器仿真模型所处的仿真工况的。
其中,汽水分离器仿真模型的仿真工况包括上水启动状态、湿态运行状态、干态运行状态、干湿态转换状态和超临界状态。其中,上水启动状态对应于蒸汽量为零的状态;湿态运行状态对应于汽水分离器仿真模型的入口工质和内部工质均为汽水混合物的状态,或者对应于入口工质仅为微过热蒸汽并且内部工质还包含水的状态;干态运行状态对应于入口工质和内部工质均仅为过热蒸汽的状态;干湿态转换状态对应于入口工质仅为过热蒸汽但内部工质还包含水的状态,或者对应于入口工质为汽水混合物并且内部工质不含水的状态;超临界状态对应于内部压力大于水的临界压力并且内部温度大于水的临界温度的状态。
需要说明的是,上述过热蒸汽、微过热蒸汽以及饱和蒸汽的温度大小关系为:过热蒸汽>微过热蒸汽>饱和蒸汽。
本发明实施例中,汽水分离器建模方法包括:基于MUSE仿真环境,对汽水分离器进行仿真模型搭建,得到相应的汽水分离器仿真模型;依据质量守恒定律和能量守恒定律,并结合工质物理特性,对汽水分离器仿真模型进行数学建模,得到汽水分离器的机理模型;当汽水分离器仿真模型处于运行仿真过程,则实时分析汽水分离器仿真模型当前所处的仿真工况;在当前参数计算周期内,利用机理模型,对汽水分离器仿真模型当前所处的仿真工况进行参数计算,得到当前参数计算周期下与该仿真工况对应的工质参数;其中,汽水分离器仿真模型的仿真工况包括上水启动状态、湿态运行状态、干态运行状态、干湿态转换状态和超临界状态。可见,本发明实施例首先基于MUSE仿真环境,构建了汽水分离器仿真模型,然后对该仿真模型进行数学建模,得到汽水分离器的机理模型;利用该机理模型,可计算得到与上述汽水分离器仿真模型当前所处的仿真工况所对应的工质参数,其中,汽水分离器仿真模型的仿真工况包括上水启动状态、湿态运行状态、干态运行状态、干湿态转换状态和超临界状态,这几种工况能够全面地涵盖了汽水分离器实际工况的种类。由上可见,本发明实施例实现了对直流锅炉中汽水分离器的仿真建模。
本发明实施例公开了一种具体的直流锅炉的汽水分离器建模方法,相对于上一实施例,本实施例对技术方案作了进一步的说明和优化。具体的:
上一实施例步骤S12中,依据质量守恒定律和能量守恒定律,并结合工质物理特性,对汽水分离器仿真模型进行数学建模的过程,具体可以包括:利用工质质量守恒方程、工质能量守恒方程、水位下工质质量方程和体积守恒方程,对汽水分离器仿真模型进行数学建模,得到机理模型。
其中,上述工质质量守恒方程为:
d d τ [ ρ s V s + ρ w V w ] = q r - q s - q a f - q c - q n ;
上述工质能量守恒方程为:
d d τ [ ρ s h s V s + ρ w h w V w - pV t + m t C p t m ] = q r h r - q s h s - q a f h a f - q c h c - q n h n ;
上述水位下工质质量方程为:
d d τ [ ρ s V s d + ρ w V w d ] = x · q r · ( 1 - η ) + ( 1 - x ) · q r - q s d - q c ;
上述体积守恒方程为:
Vt=Vs+Vw
式中,下标s表示蒸汽,下标w表示水,下标r表示入口,下标c表示出水口,下标n表示排污口,下标af表示安全阀门;符号ρ表示工质密度,符号V表示工质体积,符号q表示工质流量,符号h表示工质焓;也即,ρs表示蒸汽密度,Vs表示蒸汽体积,ρw表示水密度,Vw表示水体积,qr表示入口工质流量,qs表示蒸汽流量,qaf表示安全阀门工质流量,qc表示出水口工质流量,qn表示排污口工质流量,hs表示蒸汽焓,hw表示水焓,hr表示入口工质焓,haf表示安全阀门工质焓,hc表示出水口工质焓,hn表示排污口工质焓。另外,参数p表示压力参数,参数mt表示分离器金属质量,参数Cp表示金属比热,参数tm表示金属温度,参数Vsd表示水位下蒸汽体积,参数Vwd表示水位下水体积,参数qsd表示由温度过高引起的从水中溢出的饱和蒸汽流量,参数x表示入口工质干度,参数η表示与汽水分离器分离效率有关的预设系数,参数Vt表示汽水分离器的整体内部容积。
参见图2所示,上一实施例步骤S13中,在当前参数计算周期内,利用机理模型,对汽水分离器仿真模型当前所处的仿真工况进行参数计算的过程,包括下面步骤S131至133;具体的:
步骤S131:对上述工质质量守恒方程、工质能量守恒方程、水位下工质质量方程和体积守恒方程进行化简代换处理,得到压力变化率方程、水体积变化率方程和水位下蒸汽体积变化率方程。
其中,上述压力变化率方程具体为:
d p d τ = e 21 · A 2 - e 11 · A 1 - ( e 21 · e 13 - e 11 · e 23 ) dh s d τ ( e 21 · e 12 - e 11 · e 22 )
上述水体积变化率方程具体为:
dV w d τ = e 22 · A 2 - e 12 · A 1 - ( e 22 · e 13 - e 12 · e 23 ) dh s d τ ( e 22 · e 11 - e 12 · e 21 )
上述水位下蒸汽体积变化率方程具体为:
dV s d d τ = A 3 - ρ w dV w d τ - e 32 d p d τ - e 33 dh s d τ e 31
式中,A1=qrhr-qshs-qafhaf-qchc-qnhn;A2=qr-qs-qaf-qc-qn
A3=x·qr·(1-η)+(1-x)·qr-qsd-qc;e11=ρws
e21=ρwhwshs
e 23 = ρ s V s + h s V s ∂ ρ s ∂ h s + m t C p ∂ t m ∂ h s , e 32 = V s d ∂ ρ s ∂ p + V w d ∂ ρ w ∂ p , e 33 = V s d ∂ ρ s ∂ h s ;
Vt=Vs+Vw,Vwd=Vw
步骤S132:利用上述压力变化率方程、水体积变化率方程和水位下蒸汽体积变化率方程,分别计算得到当前参数计算周期内的压力变化率、水体积变化率和水位下蒸汽体积变化率。
步骤S133:利用当前参数计算周期内的压力变化率、水体积变化率和水位下蒸汽体积变化率,并结合上一参数计算周期内得到的压力参数、水体积参数和水位下蒸汽体积参数,计算得到当前参数计算周期内的压力参数、水体积参数和水位下蒸汽体积参数。
其中,上述步骤S133中涉及到的具体计算公式为:
V w , τ + 1 = V w , τ + dV w d τ · Δ τ ;
p τ + 1 = p τ + d p d τ · Δ τ ;
V s d , τ + 1 = V s d , τ + dV s d d τ · Δ τ .
式中,Vw,τ+1表示当前参数计算周期内的水体积参数,Vw,τ表示上一参数计算周期内的水体积参数,pτ+1表示当前参数计算周期内的压力参数,pτ表示上一参数计算周期内的压力参数,Vsd,τ+1表示当前参数计算周期内的水位下蒸汽体积参数,Vsd,τ表示上一参数计算周期内的水位下蒸汽体积参数。
进一步的,本实施例中的汽水分离器建模方法,还可以包括:利用当前参数计算周期内的水体积参数和水位下蒸汽体积参数,计算当前参数计算周期内的水位参数。具体的,利用当前参数计算周期内的水体积参数与水位下蒸汽体积参数之和除以汽水分离器仿真模型的横截面积,便可得到当前参数计算周期内的水位参数。
另外,本实施例中的汽水分离器建模方法,还可以进一步包括:
若当前参数计算周期内得到的压力参数为饱和压力,则直接利用该压力参数对工质特性参数表进行查表处理,得到相应的工质焓参数和工质密度参数;若当前参数计算周期内得到的压力参数为过热压力,则利用该压力参数以及当前参数计算周期内采集到的温度参数对工质特性参数表进行查表处理,得到相应的工质焓参数和工质密度参数;
其中,上述工质特性参数表是本领域技术人员常用的数据表,该表具体包括工程热力学中水蒸气特性参数表和水特性参数表。
进一步的,上述步骤S131中,在对工质质量守恒方程、工质能量守恒方程、水位下工质质量方程和体积守恒方程进行化简代换处理的过程之前,还包括:
从工质质量守恒方程、工质能量守恒方程、水位下工质质量方程和体积守恒方程中确定出在汽水分离器仿真模型当前所处的仿真工况下适合被简化的优化项;利用优化项,对工质质量守恒方程、工质能量守恒方程、水位下工质质量方程和体积守恒方程中相应的方程进行简化处理。
例如,若当前所处的仿真工况为上水启动状态,由于这种状态下蒸汽量为零,所以在进行参数计算时,可以将涉及到蒸汽部分的项从上述工质质量守恒方程、工质能量守恒方程、水位下工质质量方程和体积守恒方程中去除,具体为:令Vs=0,Vsd=0;
若当前所处的仿真工况为对应于汽水分离器仿真模型的入口工质和内部工质均为汽水混合物的状态,这种状态下汽水分离器仿真模型处于标准的湿态运行,此时工质参数均为饱和参数,为了简化参数获取过程,可以直接利用上一参数计算周期内得到的压力参数来查阅水蒸气特性参数表和水特性参数表,从而方便地得到与上一参数计算周期内得到的压力参数相对应的蒸汽密度、蒸汽焓、水密度和水焓参数;
若当前所处的仿真工况为对应于入口工质仅为微过热蒸汽并且内部工质还包含水的状态,此时可以在方程中考虑微过热蒸汽的参数变化情况,以此来简化参数获取过程,具体通过利用压力参数、工质焓、温度来查阅水蒸气特性参数表和水特性参数表,从而方便地得到与相应的蒸汽密度、蒸汽焓、水密度和水焓参数;
若当前所处的仿真工况为干态运行状态,由于这种状态下入口处和内部均没有水,所以在进行参数计算时,可以将涉及到水部分的项从上述工质质量守恒方程、工质能量守恒方程、水位下工质质量方程和体积守恒方程中去除,具体为:Vw=0,Vwd=0;
若当前所处的仿真工况为对应于入口工质仅为过热蒸汽但内部工质还包含水的状态,这种情况下,汽水分离器内部的蒸汽将从饱和蒸汽变化至过热蒸汽,为了保证焓的连续性,需令:
若当前所处的仿真工况为超临界状态,在这种情况下,当汽水分离器仿真模型内的工质焓低于饱和水焓,此时内部工质干度为0,当汽水分离器仿真模型内的工质焓等于或高于饱和汽焓,此时内部工质干度为1,其中,上述饱和水焓和饱和汽焓可以通过查找上述工质特性参数表得到。
相应的,本发明实施例还公开了一种直流锅炉的汽水分离器建模系统,参见图3所示,该系统包括:
仿真模型搭建模块31,用于基于MUSE仿真环境,对汽水分离器进行仿真模型搭建,得到相应的汽水分离器仿真模型;
机理模型构建模块32,用于依据质量守恒定律和能量守恒定律,并结合工质物理特性,对汽水分离器仿真模型进行数学建模,得到汽水分离器的机理模型;
仿真工况分析模块33,用于当汽水分离器仿真模型处于运行仿真过程,则实时分析汽水分离器仿真模型当前所处的仿真工况;
工质参数计算模块34,用于在当前参数计算周期内,利用机理模型,对汽水分离器仿真模型当前所处的仿真工况进行参数计算,得到当前参数计算周期下与该仿真工况对应的工质参数;
其中,汽水分离器仿真模型的仿真工况包括上水启动状态、湿态运行状态、干态运行状态、干湿态转换状态和超临界状态。
可见,本发明实施例首先基于MUSE仿真环境,构建了汽水分离器仿真模型,然后对该仿真模型进行数学建模,得到汽水分离器的机理模型;利用该机理模型,可计算得到与上述汽水分离器仿真模型当前所处的仿真工况所对应的工质参数,其中,汽水分离器仿真模型的仿真工况包括上水启动状态、湿态运行状态、干态运行状态、干湿态转换状态和超临界状态,这几种工况能够全面地涵盖了汽水分离器实际工况的种类。由上可见,本发明实施例实现了对直流锅炉中汽水分离器的仿真建模。
本发明实施例公开了一种具体的直流锅炉的汽水分离器建模系统,相对于上一实施例,本实施例对技术方案作了进一步的说明和优化。具体的
上一实施例中的仿真模型搭建模块,具体用于利用MUSE仿真环境,分别对汽水分离器的容器主体、入口、出汽口、循环水出口、排污口和安全阀门进行仿真建模,得到汽水分离器仿真模型。
上一实施例中的机理模型构建模块,具体用于利用工质质量守恒方程、工质能量守恒方程、水位下工质质量方程和体积守恒方程,对汽水分离器仿真模型进行数学建模,得到机理模型。
其中,上述工质质量守恒方程为:
d d τ [ ρ s V s + ρ w V w ] = q r - q s - q a f - q c - q n ;
上述工质能量守恒方程为:
d d τ [ ρ s h s V s + ρ w h w V w - pV t + m t C p t m ] = q r h r - q s h s - q a f h a f - q c h c - q n h n ;
上述水位下工质质量方程为:
d d τ [ ρ s V s d + ρ w V w d ] = x · q r · ( 1 - η ) + ( 1 - x ) · q r - q s d - q c ;
上述体积守恒方程为:
Vt=Vs+Vw
式中,下标s表示蒸汽,下标w表示水,下标r表示入口,下标c表示出水口,下标n表示排污口,下标af表示安全阀门;符号ρ表示工质密度,符号V表示工质体积,符号q表示工质流量,符号h表示工质焓;也即,ρs表示蒸汽密度,Vs表示蒸汽体积,ρw表示水密度,Vw表示水体积,qr表示入口工质流量,qs表示蒸汽流量,qaf表示安全阀门工质流量,qc表示出水口工质流量,qn表示排污口工质流量,hs表示蒸汽焓,hw表示水焓,hr表示入口工质焓,haf表示安全阀门工质焓,hc表示出水口工质焓,hn表示排污口工质焓。另外,参数p表示压力参数,参数mt表示分离器金属质量,参数Cp表示金属比热,参数tm表示金属温度,参数Vsd表示水位下蒸汽体积,参数Vwd表示水位下水体积,参数qsd表示由温度过高引起的从水中溢出的饱和蒸汽流量,参数x表示入口工质干度,参数η表示与汽水分离器分离效率有关的预设系数,参数Vt表示汽水分离器的整体内部容积。
另外,上一实施例的工质参数计算模块具体包括方程化简代换单元、变化率计算单元和参数计算单元;其中,
方程化简代换单元,用于对上述工质质量守恒方程、工质能量守恒方程、水位下工质质量方程和体积守恒方程进行化简代换处理,得到压力变化率方程、水体积变化率方程和水位下蒸汽体积变化率方程;
变化率计算单元,用于利用上述压力变化率方程、水体积变化率方程和水位下蒸汽体积变化率方程,分别计算得到当前参数计算周期内的压力变化率、水体积变化率和水位下蒸汽体积变化率;
参数计算单元,用于利用当前参数计算周期内的压力变化率、水体积变化率和水位下蒸汽体积变化率,并结合上一参数计算周期内得到的压力参数、水体积参数和水位下蒸汽体积参数,计算得到当前参数计算周期内的压力参数、水体积参数和水位下蒸汽体积参数。
进一步的,上述工质参数计算模块还可以包括水位计算单元,用于利用当前参数计算周期内的水体积参数和水位下蒸汽体积参数,计算当前参数计算周期内的水位参数。
另外,本实施例中的汽水分离器建模系统还可以包括:
第一查表模块,用于在当前参数计算周期内得到的压力参数为饱和压力的情况下,直接利用该压力参数对工质特性参数表进行查表处理,得到相应的工质焓参数和工质密度参数。
第二查表模块,用于在当前参数计算周期内得到的压力参数为过热压力的情况下,则利用该压力参数以及当前参数计算周期内采集到的温度参数对工质特性参数表进行查表处理,得到相应的工质焓参数和工质密度参数。
其中,上述工质特性参数表是本领域技术人员常用的数据表,该表具体包括工程热力学中水蒸气特性参数表和水特性参数表。
进一步的,本实施例中的汽水分离器建模系统还可以包括:方程优化模块,用于在上述方程化简代换单元对工质质量守恒方程、工质能量守恒方程、水位下工质质量方程和体积守恒方程进行化简代换处理的过程之前,从工质质量守恒方程、工质能量守恒方程、水位下工质质量方程和体积守恒方程中确定出在汽水分离器仿真模型当前所处的仿真工况下适合被简化的优化项,然后利用优化项,对工质质量守恒方程、工质能量守恒方程、水位下工质质量方程和体积守恒方程中相应的方程进行简化处理。
最后,还需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上对本发明所提供的一种直流锅炉的汽水分离器建模方法及系统进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。

Claims (9)

1.一种直流锅炉的汽水分离器建模方法,其特征在于,包括:
基于MUSE仿真环境,对汽水分离器进行仿真模型搭建,得到相应的汽水分离器仿真模型;
依据质量守恒定律和能量守恒定律,并结合工质物理特性,对所述汽水分离器仿真模型进行数学建模,得到所述汽水分离器的机理模型;
当所述汽水分离器仿真模型处于运行仿真过程,则实时分析所述汽水分离器仿真模型当前所处的仿真工况;
在当前参数计算周期内,利用所述机理模型,对所述汽水分离器仿真模型当前所处的仿真工况进行参数计算,得到当前参数计算周期下与该仿真工况对应的工质参数;
其中,所述汽水分离器仿真模型的仿真工况包括上水启动状态、湿态运行状态、干态运行状态、干湿态转换状态和超临界状态。
2.根据权利要求1所述的直流锅炉的汽水分离器建模方法,其特征在于,所述基于MUSE仿真环境,对汽水分离器进行仿真模型搭建的过程,包括:
利用所述MUSE仿真环境,分别对所述汽水分离器的容器主体、入口、出汽口、循环水出口、排污口和安全阀门进行仿真建模,得到所述汽水分离器仿真模型。
3.根据权利要求1所述的直流锅炉的汽水分离器建模方法,其特征在于,所述依据质量守恒定律和能量守恒定律,并结合工质物理特性,对所述汽水分离器仿真模型进行数学建模的过程,包括:
利用工质质量守恒方程、工质能量守恒方程、水位下工质质量方程和体积守恒方程,对所述汽水分离器仿真模型进行数学建模,得到所述机理模型。
4.根据权利要求3所述的直流锅炉的汽水分离器建模方法,其特征在于,
所述工质质量守恒方程为:
d d τ [ ρ s V s + ρ w V w ] = q r - q s - q a f - q c - q n ;
所述工质能量守恒方程为:
d d τ [ ρ s h s V s + ρ w h w V w - pV t + m t C p t m ] = q r h r - q s h s - q a f h a f - q c h c - q n h n ;
所述水位下工质质量方程为:
d d τ [ ρ s V s d + ρ w V w d ] = x · q r · ( 1 - η ) + ( 1 - x ) · q r - q s d - q c ;
所述体积守恒方程为:
Vt=Vs+Vw
式中,下标s表示蒸汽,下标w表示水,下标r表示入口,下标c表示出水口,下标n表示排污口,下标af表示安全阀门;符号ρ表示工质密度,符号V表示工质体积,符号q表示工质流量,符号h表示工质焓;参数p表示压力参数,参数mt表示分离器金属质量,参数Cp表示金属比热,参数tm表示金属温度,参数Vsd表示水位下蒸汽体积,参数Vwd表示水位下水体积,参数qsd表示由温度过高引起的从水中溢出的饱和蒸汽流量,参数x表示入口工质干度,参数η表示与汽水分离器分离效率有关的预设系数,参数Vt表示汽水分离器的整体内部容积。
5.根据权利要求4所述的直流锅炉的汽水分离器建模方法,其特征在于,所述在当前参数计算周期内,利用所述机理模型,对所述汽水分离器仿真模型当前所处的仿真工况进行参数计算的过程,包括:
对所述工质质量守恒方程、所述工质能量守恒方程、所述水位下工质质量方程和所述体积守恒方程进行化简代换处理,得到压力变化率方程、水体积变化率方程和水位下蒸汽体积变化率方程;
利用所述压力变化率方程、所述水体积变化率方程和所述水位下蒸汽体积变化率方程,分别计算得到当前参数计算周期内的压力变化率、水体积变化率和水位下蒸汽体积变化率;
利用当前参数计算周期内的压力变化率、水体积变化率和水位下蒸汽体积变化率,并结合上一参数计算周期内得到的压力参数、水体积参数和水位下蒸汽体积参数,计算得到当前参数计算周期内的压力参数、水体积参数和水位下蒸汽体积参数。
6.根据权利要求5所述的直流锅炉的汽水分离器建模方法,其特征在于,还包括:
利用当前参数计算周期内的水体积参数和水位下蒸汽体积参数,计算当前参数计算周期内的水位参数。
7.根据权利要求6所述的直流锅炉的汽水分离器建模方法,其特征在于,还包括:
若当前参数计算周期内得到的压力参数为饱和压力,则直接利用该压力参数对工质特性参数表进行查表处理,得到相应的工质焓参数和工质密度参数;
若当前参数计算周期内得到的压力参数为过热压力,则利用该压力参数以及当前参数计算周期内采集到的温度参数对所述工质特性参数表进行查表处理,得到相应的工质焓参数和工质密度参数;
其中,所述工质特性参数表包括工程热力学中水蒸气特性参数表和水特性参数表。
8.根据权利要求5所述的直流锅炉的汽水分离器建模方法,其特征在于,在对所述工质质量守恒方程、所述工质能量守恒方程、所述水位下工质质量方程和所述体积守恒方程进行化简代换处理的过程之前,还包括:
从所述工质质量守恒方程、所述工质能量守恒方程、所述水位下工质质量方程和所述体积守恒方程中确定出在所述汽水分离器仿真模型当前所处的仿真工况下适合被简化的优化项;
利用所述优化项,对所述工质质量守恒方程、所述工质能量守恒方程、所述水位下工质质量方程和所述体积守恒方程中相应的方程进行简化处理。
9.一种直流锅炉的汽水分离器建模系统,其特征在于,包括:
仿真模型搭建模块,用于基于MUSE仿真环境,对汽水分离器进行仿真模型搭建,得到相应的汽水分离器仿真模型;
机理模型构建模块,用于依据质量守恒定律和能量守恒定律,并结合工质物理特性,对所述汽水分离器仿真模型进行数学建模,得到所述汽水分离器的机理模型;
仿真工况分析模块,用于当所述汽水分离器仿真模型处于运行仿真过程,则实时分析所述汽水分离器仿真模型当前所处的仿真工况;
工质参数计算模块,用于在当前参数计算周期内,利用所述机理模型,对所述汽水分离器仿真模型当前所处的仿真工况进行参数计算,得到当前参数计算周期下与该仿真工况对应的工质参数;
其中,所述汽水分离器仿真模型的仿真工况包括上水启动状态、湿态运行状态、干态运行状态、干湿态转换状态和超临界状态。
CN201610479586.3A 2016-06-23 2016-06-23 一种直流锅炉的汽水分离器建模方法及系统 Active CN106126822B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610479586.3A CN106126822B (zh) 2016-06-23 2016-06-23 一种直流锅炉的汽水分离器建模方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610479586.3A CN106126822B (zh) 2016-06-23 2016-06-23 一种直流锅炉的汽水分离器建模方法及系统

Publications (2)

Publication Number Publication Date
CN106126822A true CN106126822A (zh) 2016-11-16
CN106126822B CN106126822B (zh) 2018-08-17

Family

ID=57266710

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610479586.3A Active CN106126822B (zh) 2016-06-23 2016-06-23 一种直流锅炉的汽水分离器建模方法及系统

Country Status (1)

Country Link
CN (1) CN106126822B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108469744A (zh) * 2018-02-11 2018-08-31 东南大学 一种建立核电机组蒸汽发生器机理模型的方法及其系统
CN113161583A (zh) * 2021-03-30 2021-07-23 中通客车股份有限公司 一种燃料电池汽水分离器性能评价方法及系统
CN114282389A (zh) * 2021-12-30 2022-04-05 西安交通大学 一种基于matlab的离心风机系统动态仿真方法和系统
CN114646051A (zh) * 2022-03-17 2022-06-21 国网湖南省电力有限公司 超临界火电机组湿态运行锅炉给水自动控制方法及系统
CN118274656A (zh) * 2024-06-03 2024-07-02 山东科技大学 基于参数辨识的一类汽水换热器报警监测方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102436177A (zh) * 2011-10-24 2012-05-02 长春工业大学 一种多效蒸馏水机的自动控制方法
JP2014077759A (ja) * 2012-10-12 2014-05-01 Bridgestone Corp シミュレーション方法及びシミュレーション装置
CN104833698A (zh) * 2015-04-27 2015-08-12 广东电网有限责任公司电力科学研究院 一种燃煤电站锅炉烟气酸露点检测系统及检测方法
CN104979025A (zh) * 2015-06-30 2015-10-14 浙江工业大学 一种反应堆压力容器ivr条件下结构完整性试验平台

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102436177A (zh) * 2011-10-24 2012-05-02 长春工业大学 一种多效蒸馏水机的自动控制方法
JP2014077759A (ja) * 2012-10-12 2014-05-01 Bridgestone Corp シミュレーション方法及びシミュレーション装置
CN104833698A (zh) * 2015-04-27 2015-08-12 广东电网有限责任公司电力科学研究院 一种燃煤电站锅炉烟气酸露点检测系统及检测方法
CN104979025A (zh) * 2015-06-30 2015-10-14 浙江工业大学 一种反应堆压力容器ivr条件下结构完整性试验平台

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
唐贞鹏 等: "压水堆核电机组一次调频动态仿真", 《电网技术》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108469744A (zh) * 2018-02-11 2018-08-31 东南大学 一种建立核电机组蒸汽发生器机理模型的方法及其系统
CN108469744B (zh) * 2018-02-11 2021-07-27 东南大学 一种建立核电机组蒸汽发生器机理模型的方法及其系统
CN113161583A (zh) * 2021-03-30 2021-07-23 中通客车股份有限公司 一种燃料电池汽水分离器性能评价方法及系统
CN114282389A (zh) * 2021-12-30 2022-04-05 西安交通大学 一种基于matlab的离心风机系统动态仿真方法和系统
CN114282389B (zh) * 2021-12-30 2024-01-16 西安交通大学 一种基于matlab的离心风机系统动态仿真方法和系统
CN114646051A (zh) * 2022-03-17 2022-06-21 国网湖南省电力有限公司 超临界火电机组湿态运行锅炉给水自动控制方法及系统
CN114646051B (zh) * 2022-03-17 2023-06-23 国网湖南省电力有限公司 超临界火电机组湿态运行锅炉给水自动控制方法及系统
CN118274656A (zh) * 2024-06-03 2024-07-02 山东科技大学 基于参数辨识的一类汽水换热器报警监测方法及系统

Also Published As

Publication number Publication date
CN106126822B (zh) 2018-08-17

Similar Documents

Publication Publication Date Title
CN106126822B (zh) 一种直流锅炉的汽水分离器建模方法及系统
Hoseinzadeh et al. Advanced energy, exergy, and environmental (3E) analyses and optimization of a coal-fired 400 MW thermal power plant
CN105303032B (zh) 影响发电机组能效的客观因素分析方法
Sindareh-Esfahani et al. Cold start-up condition model for heat recovery steam generators
CN109376858B (zh) 一种基于部分负荷率对冷凝式换热器性能预测的方法
CN105184395A (zh) 含余热利用系统的火电机组的初参数确定方法
CN102749156A (zh) 一种汽轮机排汽焓的检测方法
CA2479238A1 (en) Method for detecting heat exchanger tube failures and their location when using input/loss performance monitoring of a power plant
CN203376244U (zh) 抗高温水蒸气氧化性能模拟和加速测试装置
Jamnani et al. Energy, exergy, environmental (3E) and parametric assessment of a triple-pressure reheat combined-cycle power plant
CN106289754B (zh) 一种基于统计控制技术的火电机组部件性能变化判定方法
RU2621422C2 (ru) Система и способ тестирования показателя работы паровой турбины
Plis et al. Modelling and simulation of the effect of ambient parameters on the performance of a combined cycle gas turbine power plant
Duarte-Forero et al. Comparative analysis of intelligence optimization algorithms in the thermo-economic performance of an energy recovery system based on organic rankine cycle
CN104459542A (zh) 再热回热复合循环机组的热耗率测定方法
Plis et al. Adaptive simulation model of a double-pressure heat recovery steam generator for current optimization in control systems
Hercus et al. Robust shape optimization of uncertain dense gas flows through a plane turbine cascade
Pierobon et al. Application of unscented Kalman filter for condition monitoring of an organic Rankine cycle turbogenerator
CN111079070B (zh) 一种热力参数分析方法及装置
Santos et al. Capital formation, energy consumption and economic growth in Brazil: An ARDL bound testing approach and Granger causality
CN114424128A (zh) 非线性模型的建模方法、装置及存储介质
Azam et al. Evaluation of hybrid ocean thermal energy conversion system plantwide performance
CN114424196A (zh) 非线性模型的线性化处理方法、装置及存储介质
Kishore et al. Temperature Prediction Algorithms Using Machine Learning for Electric Vehicles
CN109885802B (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200306

Address after: 510080 Dongfeng East Road, Dongfeng, Guangdong, Guangzhou, Zhejiang Province, No. 8

Patentee after: ELECTRIC POWER RESEARCH INSTITUTE OF GUANGDONG POWER GRID Co.,Ltd.

Address before: 510080 Dongfeng East Road, Dongfeng, Guangdong, Guangzhou, Zhejiang Province, No. 8

Co-patentee before: UNIGROUP(BEIJING) INTELLIGENT CONTROL SCIENCE&TECHNOLOGY Co.,Ltd.

Patentee before: Electric Power Research Institute of Guangdong Power Grid Co.,Ltd.