CN108021742A - 一种考虑疏水模型的蒸汽供热管道稳态运行状态估计方法 - Google Patents
一种考虑疏水模型的蒸汽供热管道稳态运行状态估计方法 Download PDFInfo
- Publication number
- CN108021742A CN108021742A CN201711181127.8A CN201711181127A CN108021742A CN 108021742 A CN108021742 A CN 108021742A CN 201711181127 A CN201711181127 A CN 201711181127A CN 108021742 A CN108021742 A CN 108021742A
- Authority
- CN
- China
- Prior art keywords
- steam
- heat supply
- supply pipeline
- steam heat
- pipeline
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pipeline Systems (AREA)
Abstract
本发明涉及一种考虑疏水模型的蒸汽供热管道稳态运行状态估计方法，属于综合能源系统的运行和控制技术领域。首先将蒸汽供热管道和疏水阀等效成两条管道，建立考虑疏水模型的蒸汽供热管道的Γ型等效模型，再建立考虑疏水模型的蒸汽供热管道稳态运行状态估计的目标函数，根据目标函数对蒸汽供热管道稳态运行状态进行估计。本发明方法考虑了蒸汽供热管道中由于部分蒸汽凝结产生的疏水现象，可以有效的监测蒸汽供热管道的运行状况，对疏水量给出定量估计，在非全量测配置下补全量测，辨识坏数据，为能量管理系统和调度管理系统提供详实的数据支撑。
Description
技术领域
本发明涉及一种考虑疏水模型的蒸汽供热管道稳态运行状态估计方法，属于综合能源系统的运行和控制技术领域。
背景技术
热网是综合能源系统中非常重要的一部分，目前已有很多研究通过利用能源网络中热网的灵活性提高新能源的渗透率和能源利用率，在这些研究中均考虑热水为热网供热介质，但是在很多工业园区热网会选择高温高压蒸汽作为供热介质。相比较于热水管网，蒸汽管网的传输过程更加复杂，成为利用蒸汽管网的灵活性对综合能源系统进行组合分析和优化的一大障碍。
目前已经有一部分的研究涉及蒸汽管网的建模和状态估计，但是在目前的研究中还存在一定的问题，在蒸汽传输过程中有一部分蒸汽会发生凝结产生疏水现象。通常情况下，当蒸汽流速很快时疏水量可忽略不计，然而当管道足够长时凝结水可能会引起一些严重事故。目前关于蒸汽供热管道的建模和状态估计技术，大多未考虑疏水阀门的模型，或者无法对疏水量给出定量的估计，因而无法对蒸汽供热管道进行准确、有效的监控和评估，本发明针对目前蒸汽供热管道建模和状态估计技术对存在的疏水现象考虑不足的漏洞提出一种考虑疏水模型的蒸汽供热管道稳态运行状态估计方法，对疏水量给出定量估计，有效监测蒸汽供热管道的运行状况，在非全量测配置下补全量测，辨识坏数据。
发明内容
本发明的目的是提出一种考虑疏水模型的蒸汽供热管道稳态运行状态估计方法，在进行进行状态估计时，充分考虑蒸汽传输过程中由于部分蒸汽发生凝结而产生定疏水现象，对疏水量给出定量估计，以达到有效监测蒸汽供热管道的运行状况，在非全量测配置下补全量测，辨识坏数据的目的。
本发明提出的考虑疏水模型的蒸汽供热管道稳态运行状态估计方法，包括以下步骤：
(1)将蒸汽供热管道和疏水阀等效成两条管道，建立一个考虑疏水模型的蒸汽供热管道的Γ型等效模型，具体过程如下：
(1-1)根据蒸汽供热管道内蒸汽的流量变化，得到一个蒸汽供热管道流量守恒方程：
m1-m2＝md
其中，m1为蒸汽供热管道首端流量，m2为蒸汽供热管道末端流量，md为疏水阀等效成的管道内蒸汽流量；
(1-2)根据蒸汽供热管道内蒸汽的压力变化，得到一个蒸汽供热管道动量守恒方程：
其中，p1为蒸汽供热管道首端压力，p2为蒸汽供热管道末端压力，d为蒸汽供热管道的直径，f为蒸汽供热管道的摩擦阻力，从蒸汽管网设计说明书中获取，为供热管道内蒸汽的平均密度，从不同工况下蒸汽密度表中获取，为供热管道内蒸汽的平均流速， 为供热管道内蒸汽的平均流量，S为蒸汽供热管道的横截面积；
(1-3)根据蒸汽供热管道内蒸汽的焓变化，得到一个蒸汽供热管道能量守恒方程：
其中，H1为蒸汽供热管道首端蒸汽的焓，H2为蒸汽供热管道末端蒸汽的焓，l为蒸汽供热管道的长度，S为上述步骤(1-2)中蒸汽供热管道的横截面积，d为上述步骤(1-2)中蒸汽供热管道的直径，λ为蒸汽供热管道的散热系数，从相应的数据手册中获取，为蒸汽供热管道内平均温度，通过蒸汽供热管道首末端温度求出，T1为管道首端温度，T2为管道末端温度，TA为环境温度，Hd为进入疏水阀等效成的管道内的蒸汽的焓，Hd＝md·hd，md为上述步骤(1-1)中疏水阀等效成的管道内蒸汽流量，hd为进入疏水阀等效成的管道内的蒸汽的比焓，hd＝hd(pd,Td)，pd为疏水阀等效成的管道内蒸汽的压力，Td为疏水阀等效成的管道内蒸汽的温度，Qd为产生疏水现象过程中的热损失，Qd＝cp·md(T2-Td)+md·γd，cp为蒸汽的比热容，由流体的物性参数表获取，γd为蒸汽的潜热系数，由蒸汽汽化潜热值表获取；
令Td≈T2，pd＝p2，p2为蒸汽供热管道末端压力，得到蒸汽供热管道能量守恒方程的最终表达形式如下：
其中，
(2)建立一个考虑疏水模型的蒸汽供热管道稳态运行状态估计的目标函数，如下：
minJ(xh)＝min{[zh-f(xh)]TW-1[zh-f(xh)]}
其中，W为测量值的协方差矩阵，上标T表示矩阵转置，J(xh)表示目标函数表达式，zh为蒸汽供热管道量测值，包括蒸汽管道首端温度m1、首端压力p1、首端温度T1、末端温度m2、末端压力p2、末端温度T2，zh＝[m1,p1,T1,m2,p2,T2]T，上述量测值从热网的数据采集与监视控制系统中实时测量获取，xh为蒸汽供热管道状态量，xh＝[m1,p1,T1,md]T，m1、p1、T1、md分别为蒸汽供热管道首端流量、蒸汽供热管道首端压力、蒸汽供热管道首端温度和疏水阀内流量，f(x)为描述蒸汽管道状态量与测量值之间关系的量测函数，f(x)＝f(xh)，量测方程如下：
其中，h1为供热管道首端蒸汽的比焓，h1＝h1(p1,T1)，h2为供热管道末端蒸汽的比焓，h2＝h2(p2,T2)；
(3)对上述步骤(2)的目标函数进行状态估计，设定蒸汽供热管道状态估计的收敛精度δ和最大循环次数d，初始化时设循环次数k为0，选用无约束的牛拉法进行估计，该方法取最优解的条件为：
进行第k次迭代式为：
其中，J(xh)为上述步骤(2)中的目标函数，f(xh)为上述步骤(2)中的量测函数，xh为上述步骤(2)中的状态量，上标T表示矩阵转置，W为上述步骤(2)中测量值的协方差矩阵，zh为上述步骤(2)中的量测值；
(4)对上述步骤(3)的状态估计结果进行收敛性判断：
若循环次数k达到预设循环次数d，即k≥d，则将本次状态估计结果作为t时刻考虑疏水模型的蒸汽供热管道稳态运行状态估计结果；
若循环次数k未达到预设循环次数d，即k<d，则进一步根据蒸汽供热管道状态估计的精度δ对状态估计结果收敛性进行判断：若最近相邻两次状态估计结果中的状态量估计值x(k)和x(k-1)的差值小于状态估计精度δ，即则将本次状态估计结果作为t时刻考虑疏水模型的蒸汽供热管道稳态运行状态估计结果，若最近两次状态估计结果中的状态变量估计值x(k)和x(k-1)的差值大于状态估计精度δ，即max|x(k)-x(k-1)|≥δ，则更新状态变量，并根据本次状态估计所得的温度值更新蒸汽供热管道首端流量、压力、温度和疏水阀内流量，同时使k＝k+1，并返回步骤(3)，继续本次状态估计过程。
本发明提出的考虑疏水模型的蒸汽供热管道稳态运行状态估计方法，弥补了对于蒸汽供热管道状态估计方面的不足，形成了考虑蒸汽供热过程中由于部分蒸汽凝结产生的疏水现象的状态估计方法，可以很好的解决蒸汽供热管道状态估计问题，当蒸汽供热管道受调度控制而改变运行方式时，本发明能够较为精确地追踪系统温度等状态变量的变化，且有较好的收敛性。考虑疏水模型的蒸汽供热管道稳态运行状态估计，可以有效的监测蒸汽供热管道的运行状况，对疏水量给出定量估计，在非全量测配置下补全量测，辨识坏数据，为能量管理系统和调度管理系统提供详实的数据支撑。
附图说明
图1为本发明方法涉及到考虑疏水模型的蒸汽供热管道示意图。
图2为图1所示的考虑疏水模型的蒸汽供热管道的Γ型等效模型示意图。
具体实施方式
本发明提出的考虑疏水模型的蒸汽供热管道稳态运行状态估计方法，其中蒸汽供热管道示意图的结构示意图如图1所示，该方法包括以下步骤：
(1)将蒸汽供热管道和疏水阀等效成两条管道，建立一个考虑疏水模型的蒸汽供热管道的Γ型等效模型，如图2所示，具体过程如下：
(1-1)根据蒸汽供热管道内蒸汽的流量变化，得到一个蒸汽供热管道流量守恒方程：
m1-m2＝md
其中，m1为蒸汽供热管道首端流量，m2为蒸汽供热管道末端流量，md为疏水阀等效成的管道内蒸汽流量；
(1-2)根据蒸汽供热管道内蒸汽的压力变化，得到一个蒸汽供热管道动量守恒方程：
其中，p1为蒸汽供热管道首端压力，p2为蒸汽供热管道末端压力，d为蒸汽供热管道的直径，f为蒸汽供热管道的摩擦阻力，从蒸汽管网设计说明书中获取，为供热管道内蒸汽的平均密度，从不同工况下蒸汽密度表中获取，为供热管道内蒸汽的平均流速，
为供热管道内蒸汽的平均流量，S为蒸汽供热管道的横截面积；
(1-3)根据蒸汽供热管道内蒸汽的焓变化，得到一个蒸汽供热管道能量守恒方程：
其中，H1为蒸汽供热管道首端蒸汽的焓，H2为蒸汽供热管道末端蒸汽的焓，l为蒸汽供热管道的长度，S为上述步骤(1-2)中蒸汽供热管道的横截面积，d为上述步骤(1-2)中蒸汽供热管道的直径，λ为蒸汽供热管道的散热系数，从相应的数据手册中获取，为蒸汽供热管道内平均温度，通过蒸汽供热管道首末端温度求出，T1为管道首端温度，T2为管道末端温度，TA为环境温度，Hd为进入疏水阀等效成的管道内的蒸汽的焓，Hd＝md·hd，md为上述步骤(1-1)中疏水阀等效成的管道内蒸汽流量，hd为进入疏水阀等效成的管道内的蒸汽的比焓，hd＝hd(pd,Td)，pd为疏水阀等效成的管道内蒸汽的压力，Td为疏水阀等效成的管道内蒸汽的温度，Qd为产生疏水现象过程中的热损失，Qd＝cp·md(T2-Td)+md·γd，cp为蒸汽的比热容，由流体的物性参数表获取，γd为蒸汽的潜热系数，由蒸汽汽化潜热值表获取。
一般情况下，疏水阀很短且疏水阀内温度和压力变化很小，因此可以近似认为疏水阀等效成的管道内温度和压力与蒸汽供热管道末端的温度和压力相等，即令Td≈T2，pd＝p2，p2为蒸汽供热管道末端压力，得到蒸汽供热管道能量守恒方程的最终表达形式如下：
其中，
(2)建立一个考虑疏水模型的蒸汽供热管道稳态运行状态估计的目标函数，如下：
minJ(xh)＝min{[zh-f(xh)]TW-1[zh-f(xh)]}
其中，W为测量值的协方差矩阵，上标T表示矩阵转置，J(xh)表示目标函数表达式，zh为蒸汽供热管道量测值，包括蒸汽管道首端温度m1、首端压力p1、首端温度T1、末端温度m2、末端压力p2、末端温度T2，zh＝[m1,p1,T1,m2,p2,T2]T，上述量测值从热网的数据采集与监视控制系统中实时测量获取，xh为蒸汽供热管道状态量，xh＝[m1,p1,T1,md]T，m1、p1、T1、md分别为蒸汽供热管道首端流量、蒸汽供热管道首端压力、蒸汽供热管道首端温度和疏水阀内流量，f(x)为描述蒸汽管道状态量与测量值之间关系的量测函数，f(x)＝f(xh)，量测方程如下：
其中，h1为供热管道首端蒸汽的比焓，h1＝h1(p1,T1)，h2为供热管道末端蒸汽的比焓，h2＝h2(p2,T2)，其余参数与上述步骤(1)中含义相同；
(3)对上述步骤(2)的目标函数进行状态估计，设定蒸汽供热管道状态估计的收敛精度δ和最大循环次数d，初始化时设循环次数k为0，选用无约束的牛拉法进行估计，该方法取最优解的条件为：
进行第k次迭代式为：
其中，J(xh)为上述步骤(2)中的目标函数，f(xh)为上述步骤(2)中的量测函数，xh为上述步骤(2)中的状态量，上标T表示矩阵转置，W为上述步骤(2)中测量值的协方差矩阵，zh为上述步骤(2)中的量测值；
(4)对上述步骤(3)的状态估计结果进行收敛性判断：
若循环次数k达到预设循环次数d，即k≥d，则将本次状态估计结果作为t时刻考虑疏水模型的蒸汽供热管道稳态运行状态估计结果；
若循环次数k未达到预设循环次数d，即k<d，则进一步根据蒸汽供热管道状态估计的精度δ对状态估计结果收敛性进行判断：若最近相邻两次状态估计结果中的状态量估计值x(k)和x(k-1)的差值小于状态估计精度δ，即则将本次状态估计结果作为t时刻考虑疏水模型的蒸汽供热管道稳态运行状态估计结果，若最近两次状态估计结果中的状态变量估计值x(k)和x(k-1)的差值大于状态估计精度δ，即max|x(k)-x(k-1)|≥δ，则更新状态变量，并根据本次状态估计所得的温度值更新蒸汽供热管道首端流量、压力、温度和疏水阀内流量，同时使k＝k+1，并返回步骤(3)，继续本次状态估计过程。
Claims (1)
1.一种考虑疏水模型的蒸汽供热管道稳态运行状态估计方法，其特征在于该方法包括以下步骤：
(1)将蒸汽供热管道和疏水阀等效成两条管道，建立一个考虑疏水模型的蒸汽供热管道的Γ型等效模型，具体过程如下：
(1-1)根据蒸汽供热管道内蒸汽的流量变化，得到一个蒸汽供热管道流量守恒方程：
m1-m2＝md
其中，m1为蒸汽供热管道首端流量，m2为蒸汽供热管道末端流量，md为疏水阀等效成的管道内蒸汽流量；
(1-2)根据蒸汽供热管道内蒸汽的压力变化，得到一个蒸汽供热管道动量守恒方程：
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>f</mi>
<mo>&CenterDot;</mo>
<mover>
<mi>&rho;</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<mn>2</mn>
<mi>d</mi>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<msup>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>d</mi>
<mo>&CenterDot;</mo>
<msup>
<mi>S</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mi>f</mi>
<mo>&CenterDot;</mo>
<mfrac>
<msup>
<mover>
<mi>m</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msup>
<mover>
<mi>&rho;</mi>
<mo>&OverBar;</mo>
</mover>
</mfrac>
</mrow>
其中，p1为蒸汽供热管道首端压力，p2为蒸汽供热管道末端压力，d为蒸汽供热管道的直径，f为蒸汽供热管道的摩擦阻力，从蒸汽管网设计说明书中获取，为供热管道内蒸汽的平均密度，从不同工况下蒸汽密度表中获取，为供热管道内蒸汽的平均流速， 为供热管道内蒸汽的平均流量，S为蒸汽供热管道的横截面积；
(1-3)根据蒸汽供热管道内蒸汽的焓变化，得到一个蒸汽供热管道能量守恒方程：
<mrow>
<msub>
<mi>H</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>4</mn>
<mi>l</mi>
<mo>&CenterDot;</mo>
<mi>S</mi>
</mrow>
<mi>d</mi>
</mfrac>
<mo>&CenterDot;</mo>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>T</mi>
<mo>&OverBar;</mo>
</mover>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>A</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>H</mi>
<mi>d</mi>
</msub>
<mo>-</mo>
<msub>
<mi>Q</mi>
<mi>d</mi>
</msub>
</mrow>
其中，H1为蒸汽供热管道首端蒸汽的焓，H2为蒸汽供热管道末端蒸汽的焓，l为蒸汽供热管道的长度，S为上述步骤(1-2)中蒸汽供热管道的横截面积，d为上述步骤(1-2)中蒸汽供热管道的直径，λ为蒸汽供热管道的散热系数，从相应的数据手册中获取，为蒸汽供热管道内平均温度，通过蒸汽供热管道首末端温度求出，T1为管道首端温度，T2为管道末端温度，TA为环境温度，Hd为进入疏水阀等效成的管道内的蒸汽的焓，Hd＝md·hd，md为上述步骤(1-1)中疏水阀等效成的管道内蒸汽流量，hd为进入疏水阀等效成的管道内的蒸汽的比焓，hd＝hd(pd，Td)，pd为疏水阀等效成的管道内蒸汽的压力，Td为疏水阀等效成的管道内蒸汽的温度，Qd为产生疏水现象过程中的热损失，Qd＝cp·md(T2-Td)+md·γd，cp为蒸汽的比热容，由流体的物性参数表获取，γd为蒸汽的潜热系数，由蒸汽汽化潜热值表获取；
令Td≈T2，pd＝p2，p2为蒸汽供热管道末端压力，得到蒸汽供热管道能量守恒方程的最终表达形式如下：
<mrow>
<msub>
<mi>H</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>4</mn>
<mi>l</mi>
<mo>&CenterDot;</mo>
<mi>S</mi>
</mrow>
<mi>d</mi>
</mfrac>
<mo>&CenterDot;</mo>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>T</mi>
<mo>&OverBar;</mo>
</mover>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>A</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>H</mi>
<mo>~</mo>
</mover>
<mi>d</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>Q</mi>
<mo>~</mo>
</mover>
<mi>d</mi>
</msub>
</mrow>
其中，
(2)建立一个考虑疏水模型的蒸汽供热管道稳态运行状态估计的目标函数，如下：
min J(xh)＝min{[zh-f(xh)]TW-1[zh-f(xh)]}
其中，W为测量值的协方差矩阵，上标T表示矩阵转置，J(xh)表示目标函数表达式，zh为蒸汽供热管道量测值，包括蒸汽管道首端温度m1、首端压力p1、首端温度T1、末端温度m2、末端压力p2、末端温度T2，zh＝[m1，p1，T1，m2，p2，T2]T，上述量测值从热网的数据采集与监视控制系统中实时测量获取，xh为蒸汽供热管道状态量，xh＝[m1，p1，T1，md]T，m1、p1、T1、md分别为蒸汽供热管道首端流量、蒸汽供热管道首端压力、蒸汽供热管道首端温度和疏水阀内流量，f(x)为描述蒸汽管道状态量与测量值之间关系的量测函数，f(x)＝f(xh)，量测方程如下：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>m</mi>
<mi>d</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>d</mi>
<mo>&CenterDot;</mo>
<msup>
<mi>S</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mi>f</mi>
<mo>&CenterDot;</mo>
<mfrac>
<msubsup>
<mi>m</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mover>
<mi>&rho;</mi>
<mo>&OverBar;</mo>
</mover>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mn>4</mn>
<mi>l</mi>
<mo>&CenterDot;</mo>
<mi>S</mi>
</mrow>
<mi>d</mi>
</mfrac>
<mo>&CenterDot;</mo>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>T</mi>
<mo>&OverBar;</mo>
</mover>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>A</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>m</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&gamma;</mi>
<mi>d</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中，h1为供热管道首端蒸汽的比焓，h1＝h1(p1，T1)，h2为供热管道末端蒸汽的比焓，h2＝h2(p2，T2)；
(3)对上述步骤(2)的目标函数进行状态估计，设定蒸汽供热管道状态估计的收敛精度δ和最大循环次数d，初始化时设循环次数k为0，选用无约束的牛拉法进行估计，该方法取最优解的条件为：
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>J</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>h</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>h</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
</mrow>
进行第k次迭代式为：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&Delta;x</mi>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>G</mi>
<mi>T</mi>
</msup>
<msup>
<mi>W</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>G</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>G</mi>
<mi>T</mi>
</msup>
<msup>
<mi>W</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>&Delta;z</mi>
<mi>h</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>x</mi>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>x</mi>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&Delta;x</mi>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中，J(xh)为上述步骤(2)中的目标函数，f(xh)为上述步骤(2)中的量测函数，xh为上述步骤(2)中的状态量，上标T表示矩阵转置，W为上述步骤(2)中测量值的协方差矩阵，zh为上述步骤(2)中的量测值；
(4)对上述步骤(3)的状态估计结果进行收敛性判断：
若循环次数k达到预设循环次数d，即k≥d，则将本次状态估计结果作为t时刻考虑疏水模型的蒸汽供热管道稳态运行状态估计结果；
若循环次数k未达到预设循环次数d，即k＜d，则进一步根据蒸汽供热管道状态估计的精度δ对状态估计结果收敛性进行判断：若最近相邻两次状态估计结果中的状态量估计值x(k)和x(k-1)的差值小于状态估计精度δ，即max|x(k)-x(k-1)|＜δ，则将本次状态估计结果作为t时刻考虑疏水模型的蒸汽供热管道稳态运行状态估计结果，若最近两次状态估计结果中的状态变量估计值x(k)和x(k-1)的差值大于状态估计精度δ，即max|x(k)-x(k-1)|≥δ，则更新状态变量，并根据本次状态估计所得的温度值更新蒸汽供热管道首端流量、压力、温度和疏水阀内流量，同时使k＝k+1，并返回步骤(3)，继续本次状态估计过程。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711181127.8A CN108021742B (zh) | 2017-11-23 | 2017-11-23 | 一种考虑疏水模型的蒸汽供热管道稳态运行状态估计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711181127.8A CN108021742B (zh) | 2017-11-23 | 2017-11-23 | 一种考虑疏水模型的蒸汽供热管道稳态运行状态估计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108021742A true CN108021742A (zh) | 2018-05-11 |
| CN108021742B CN108021742B (zh) | 2019-11-08 |
Family
ID=62080758
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711181127.8A Active CN108021742B (zh) | 2017-11-23 | 2017-11-23 | 一种考虑疏水模型的蒸汽供热管道稳态运行状态估计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108021742B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109611814A (zh) * | 2018-12-29 | 2019-04-12 | 湖南鸿远高压阀门有限公司 | 管阀布置疏水系统热能损失的测量方法及测量系统 |
| CN110362893A (zh) * | 2019-06-28 | 2019-10-22 | 浙江大学 | 一种基于物联网的供热管道保温特性及疏水监测诊断系统 |
| CN112664817A (zh) * | 2020-12-17 | 2021-04-16 | 杭州英集动力科技有限公司 | 基于蒸汽管网模型的主动式疏水系统 |
| CN112883662A (zh) * | 2021-02-01 | 2021-06-01 | 清华大学 | 一种蒸汽供热网络动态运行水力状态估计方法及系统 |
| CN115079592A (zh) * | 2022-07-12 | 2022-09-20 | 哈尔滨工程大学 | 一种船舶核动力装置热力系统管网仿真方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101825502A (zh) * | 2010-04-16 | 2010-09-08 | 东南大学 | 汽机带疏水冷却器的加热器出水及疏水温度测算方法 |
| CN104500979A (zh) * | 2014-12-02 | 2015-04-08 | 常州英集动力科技有限公司 | 供热管网内蒸汽流动滞留段在线诊断消除系统及工作方法 |
-
2017
- 2017-11-23 CN CN201711181127.8A patent/CN108021742B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101825502A (zh) * | 2010-04-16 | 2010-09-08 | 东南大学 | 汽机带疏水冷却器的加热器出水及疏水温度测算方法 |
| CN104500979A (zh) * | 2014-12-02 | 2015-04-08 | 常州英集动力科技有限公司 | 供热管网内蒸汽流动滞留段在线诊断消除系统及工作方法 |
Non-Patent Citations (3)
| Title |
|---|
| TOMASZ PORZUCEK ET AL.: "Modeling and Evaluation Framework for FMC-QE", 《2009 16TH ANNUAL IEEE INTERNATIONAL CONFERENCE AND WORKSHOP ON THE ENGINEERING OF COMPUTER BASED SYSTEMS》 * |
| 刘运雷: "蒸汽管网散热损失计算分析与负荷预测研究", 《中国优秀硕士学位论文全文数据库（电子期刊）工程科技II辑》 * |
| 王旭光: "大型工业供热蒸汽管网运行状态分析及操作优化", 《中国优秀硕士学位论文全文数据库（电子期刊） 信息科技辑》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109611814A (zh) * | 2018-12-29 | 2019-04-12 | 湖南鸿远高压阀门有限公司 | 管阀布置疏水系统热能损失的测量方法及测量系统 |
| CN109611814B (zh) * | 2018-12-29 | 2019-12-20 | 湖南鸿远高压阀门有限公司 | 管阀布置疏水系统热能损失的测量方法及测量系统 |
| CN110362893A (zh) * | 2019-06-28 | 2019-10-22 | 浙江大学 | 一种基于物联网的供热管道保温特性及疏水监测诊断系统 |
| CN112664817A (zh) * | 2020-12-17 | 2021-04-16 | 杭州英集动力科技有限公司 | 基于蒸汽管网模型的主动式疏水系统 |
| CN112883662A (zh) * | 2021-02-01 | 2021-06-01 | 清华大学 | 一种蒸汽供热网络动态运行水力状态估计方法及系统 |
| CN112883662B (zh) * | 2021-02-01 | 2024-05-10 | 清华大学 | 一种蒸汽供热网络动态运行水力状态估计方法及系统 |
| CN115079592A (zh) * | 2022-07-12 | 2022-09-20 | 哈尔滨工程大学 | 一种船舶核动力装置热力系统管网仿真方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108021742B (zh) | 2019-11-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108021742A (zh) | 一种考虑疏水模型的蒸汽供热管道稳态运行状态估计方法 | |
| CN107817681B (zh) | 一种基于双侧等效模型的热网稳态运行状态估计方法 | |
| CN102338568B (zh) | 基于清洁系数指标的电厂凝汽器性能在线监测系统及方法 | |
| WO2020015237A1 (zh) | 基于滚动时域估计理论的热网动态调节运行参数估计方法 | |
| CN106706215A (zh) | 一种热力系统阀门内漏量监测方法 | |
| CN104749999B (zh) | 装配湿式冷却塔的汽轮发电机组冷端系统优化运行精确指导系统 | |
| CN105067661A (zh) | 气-液换热器传热系数测定装置 | |
| CN104500979A (zh) | 供热管网内蒸汽流动滞留段在线诊断消除系统及工作方法 | |
| WO2015171196A1 (en) | Virtual flow measurement system | |
| CN104697704A (zh) | 冷水机组排气压力故障检测方法 | |
| CN103822098A (zh) | 供水管网泄露监控方法 | |
| CN103063354A (zh) | 火电机组能耗评估和煤耗查定试验中汽轮机基准背压的确定 | |
| CN101477071B (zh) | 浅层岩土体热物理性质原位测试仪 | |
| CN203561596U (zh) | 地埋管岩土热响应测试仪 | |
| CN109726464A (zh) | 一种土石坝流热耦合模型参数敏感性分析构建方法 | |
| WO2022160681A1 (zh) | 一种蒸汽供热网络动态运行水力状态估计方法及系统 | |
| CN114117822B (zh) | 一种基于管壁温度条件检测瞬态增压供水管道结冰厚度方法 | |
| CN109340583A (zh) | 供热管网泄漏监测系统及方法 | |
| CN106768118B (zh) | 一种湿蒸汽流量计量装置及其计算方法 | |
| CN106840733B (zh) | 一种空调机组性能在线检测方法及其检测设备 | |
| CN107133468B (zh) | 一种间接空冷塔冷却扇段进风量的在线软测量方法 | |
| CN106894853B (zh) | 凝汽式汽轮机冷端节能诊断方法 | |
| CN112413948B (zh) | 一种空气源热泵除霜控制点实验室测定系统及方法 | |
| CN211476791U (zh) | 一种卧式高压加热器运行效率实时控制系统 | |
| CN103383365B (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 |