CN114705245B - 一种核电厂低加疏水管路u型水封消失及维持的诊断方法 - Google Patents

Abstract

本发明公开了一种核电厂低加疏水管路U型水封消失及维持的诊断方法，根据低加汽侧压力、凝汽器汽侧压力、低加疏水流量计算疏水管路中工质的热力学状态(压力、干度)，对疏水管路U型水封消失的原因进行诊断并提出解决方法。本发明诊断方法可以根据低加汽侧压力、凝汽器汽侧压力、低加疏水流量判断U型水封是否存在。在U型水封消失时，通过计算疏水管路不同位置处工质的热力学状态分析水封消失的原因，并计算维持U型水封的临界注水量。

Description

一种核电厂低加疏水管路U型水封消失及维持的诊断方法

技术领域

本发明涉及核电厂低加疏水管路技术领域，具体涉及一种核电厂低加疏水管路U型水封消失及维持的诊断方法，主要为核电厂低加疏水管路U型水封高度及水封消失的原因进行诊断，并对维持U型水封的临界注水量进行计算。

背景技术

在核电机组中，低加(低压加热器)的凝结疏水通过疏水管路流入凝汽器中，疏水管路设计有U型水封以防止“串汽”的发生。但是由于疏水管路布置不合理、低加传热恶化等原因，导致水封不能维持，部分抽气未被凝结直接进入凝汽器，造成热损失，使机组经济性降低。对疏水管路U型水封高度及水封消失的原因进行诊断是很多核电厂亟待解决的问题。

现有研究主要针对低加汽侧压力、凝汽器汽侧压力、疏水管路汇合点高度等参数进行热力计算，无法对疏水管路U型水封消失的机理进行研究。

现有研究对U型水封消失的原因分析，都是在疏水动能、U型水封两端压差、U型水封两端重力势能差求解水封高度，但有些工况即便满足以上条件，仍不能维持U型水封。

发明内容

为了克服上述现有技术存在的问题，本发明的目的在于提供一种核电厂低加疏水管路U型水封消失及维持的诊断方法，对核电厂低加疏水管路U型水封高度及沿疏水管道不同位置处工质的热力学状态进行计算，对U型水封消失的原因进行诊断，可有效分析U型水封消失的根本原因并提出解决方法。

为达到上述目的，本发明采用如下技术方案：

一种核电厂低加疏水管路U型水封消失及维持的诊断方法，包括如下步骤：

步骤1：根据核电厂低加疏水管路低加汽侧压力p_a、低加疏水流量Q_m、凝汽器汽侧压力p_b、水封测量高度h_测计算低加疏水管路阻力修正系数m；

低加整个疏水管路的压降：

式中：

p_a—低加汽侧压力，Pa；

p_b—凝汽器汽侧压力，Pa；

ρ—疏水密度，是关于当地压力的函数，kg/m³；

g—当地重力加速度，m/s²

R—沿程阻力函数，Pa/m；

a—低加疏水管路起点，低压加热器；

b—低加疏水管路终点，低加疏水立管；

将疏水管路按空间离散，采用一阶迎风格式，则

p_i-1-p_i＝ρ_i-1g(z_i-z_i-1)+R_i-1(l_i-l_i-1)

整理得：

式中：

l_i—疏水管路第i个节点长度，m；

l_i-1—疏水管路第i-1个节点长度，m；

z_i—疏水管路第i个节点高度，m；

z_i-1—疏水管路第i-1个节点高度，m；

p_i—疏水管路第i个节点压力，pa；

p_i-1—疏水管路第i-1个节点压力，pa；

R_i-1—疏水管路第i-1个节点沿程阻力函数，Pa/m；

ρ_i-1—疏水管路第i-1个节点密度，kg/m³；

对低加的整个疏水管道进行建模，已知低加疏水管路的入口压力为低加汽侧压力p_a、出口压力为凝汽器汽侧压力p_b以及低加疏水流量Q_m；整个疏水管路的压降p_a-p_b受阻力修正系数m、低加疏水流量Q_m和低加疏水管路水封高度h₁影响，故

低加疏水管路阻力修正系数：

m＝f(Q_m,p_a,p_b,h₁)

式中：

Q_m—低加疏水流量，kg/s；

h₁—低加疏水管路水封高度，m；

低加疏水管路阻力修正系数m求解步骤为：

1.输入低加汽侧压力p_a、凝汽器汽侧压力p_b、低加疏水流量Q_m、测量水封高度h_测；

2.将疏水管路按n个节点划分成n-1个微元段；

3.假设低加疏水管路阻力修正系数m；

4.假设低加疏水焓H_q，p₁＝p_a，i＝2，H₁＝H_q，Q₁＝Q_m；

5.根据节点焓值H_i-1、节点压力p_i-1、节点流量Q_i-1计算得到节点密度ρ_i-1、节点沿程阻力函数R_i-1；

6.计算节点压力

H_i＝H_q，Q_i＝Q_m；

7.根据H_i、p_i计算节点疏水干度x_i，若x_i＞0，x_i-1＜0，则该节点为汽液分界面，低加疏水管路水封高度h₁＝z_i；

8.计算下一个节点，即i＝i+1；

9.如果i＞n，则进行下一步，反之返回第5步；

10.若|p_n-p_b|≤0.001，则进行下一步，反之返回第4步重新假设低加疏水焓H_q；

11.若|h₁-h_测|≤0.001，则进行下一步，反之返回第3步重新假设低加疏水管路阻力修正系数m；

12.结束循环输出低加疏水管路阻力修正系数m、低加疏水焓H_q以及低加疏水管路不同位置处工质热力学状态；

步骤2：根据低加疏水管路阻力修正系数m、低加汽侧压力p_a、低加疏水流量Q_m、凝汽器汽侧压力p_b计算低加疏水管路水封高度h₁；

通过上述计算求解获得低加疏水管路阻力修正系数m之后，根据低压加热器汽侧压力p_a、凝汽器汽侧压力p_b、低加疏水流量Q_m，获得整个疏水管路压力、干度分布，通过干度获得低加疏水管路水封高度h₁；

低加疏水管路水封高度：

h₁＝g(Q_m,m,p_a,p_b)

低加疏水管路水封高度h₁求解步骤为：

1.输入低加汽侧压力p_a、凝汽器汽侧压力p_b、低加疏水流量Q_m，阻力修正系数m；

2.将疏水管路按n个节点划分成n-1个微元段；

3.假设低加疏水焓H_q，p₁＝p_a，i＝2，H₁＝H_q，Q₁＝Q_m；

4.根据节点焓值H_i-1、节点压力p_i-1、节点流量Q_i-1计算得到节点密度ρ_i-1、节点沿程阻力函数R_i-1；

5.计算节点压力

H_i＝H_q，Q_i＝Q_m；

6.根据H_i、p_i计算节点疏水干度x_i，若x_i＞0，x_i-1＜0，则该节点为汽液分界面，低加疏水管路水封高度h₁＝z_i；

7.计算下一个节点，即i＝i+1；

8.如果i＞n，则进行下一步，反之返回第4步；

9.若|p_n-p_b|≤0.001，则进行下一步，反之返回第3步重新假设低加疏水焓H_q；

10.结束循环输出低加疏水焓H_q、低加疏水管路水封高度h₁以及低加疏水管路不同位置处工质热力学状态；

通过计算不同工况下的水封高度，当计算水封高度大于疏水管路最低点位置时，则水封存在，当计算水封高度小于疏水管路最低点位置时，水封消失；通过分析水封消失时疏水管路不同位置处工质的热力学状态，对比低加疏水管路最低点压力与低加疏水焓值对应下的饱和压力，对水封消失的原因进行诊断；

步骤3：根据低加疏水管路阻力修正系数m、低加汽侧压力p_a、低加疏水焓H_q、凝汽器汽侧压力p_b、注水焓H_b以及疏水管路最低点高度h_min，计算维持U型水封的低加疏水管路临界注水量Q_b进行计算；

在实际运行过程中，低压加热器换热不足时，低加疏水焓值增加，整个疏水管路的压力、干度发生变化，水封消失，低压加热器抽气直接进入凝汽器，潜热未被利用，造成热损失，使机组经济性降低；在核电厂中，通过注入冷却水，降低低加疏水焓值，是维持水封存在最有效、经济的方法。在U型水封底部注入冷却水，可以有效降低疏水焓值，降低汽化压力。根据低加汽侧压力p_a、凝汽器汽侧压力p_b、低加疏水焓H_q，通过迭代计算求解维持U型水封的低加疏水管路临界注水量；

维持U型水封的低加疏水管路临界注水量Q_b求解步骤为：

1.将疏水管路按n个节点划分成n-1个微元段；

2.输入低加汽侧压力p_a、凝汽器汽侧压力p_b、低加疏水焓H_q、低加疏水管路阻力修正系数m、注水焓H_b，疏水管路最低点高度h_min；

3.假设低加临界疏水量Q_m；

4.假设注水量Q_b，h₁＝0，p₁＝p_a，i＝2，H₁＝H_q，Q₁＝Q_m；

5.根据节点焓值H_i-1、节点压力p_i-1、节点流量Q_i-1计算得到节点密度ρ_i-1、节点沿程阻力函数R_i-1；

6.计算节点压力

7.疏水管路补水前，疏水焓H_i＝H_q，疏水流量Q_i＝Q_m；疏水管路补水后，疏水焓

疏水流量Q_i＝Q_m+Q_b；

8.根据H_i、p_i计算节点疏水干度x_i，若x_i＞0，x_i-1＜0，则该节点为汽液分界面，低加疏水管路水封高度h₁＝z_i；

9.计算下一个节点，即i＝i+1；

10.如果i＞n，则进行下一步，反之返回第5步；

11.若|p_n-p_b|≤0.001，则进行下一步，反之返回第4步重新假设注水量Q_b；

12.若|h₁-h_min|≤0.001，则进行下一步，反之返回第3步重新假设低加疏水流量Q_m；

13.结束循环输出维持U型水封的低加疏水管路临界注水量Q_b、低加疏水流量Q_m以及低加疏水管路不同位置处工质热力学状态。

综上所述，本发明的目的在于提供一种以核电厂机组运行过程中，低压加热器疏水管路U型水封消失原因进行诊断，并对水封消失时临界注水量进行计算。首先通过计算低加疏水管路阻力系数，然后计算不同工况下的水封高度，当计算水封高度大于疏水管路最低点位置时，则水封存在；当计算水封高度小于疏水管路最低点位置时，水封消失。通过分析水封消失时疏水管路不同位置处工质的热力学状态，对比低加疏水管路最低点压力与低加疏水焓值对应下的饱和压力，对水封消失的原因进行诊断。当水封消失时，通过向疏水管路注入冷却水维持水封存在，计算维持U型水封的临界注水量。

附图说明

图1是核电厂低加疏水水封图。

图2是低加疏水管路阻力修正系数求解流程图。

图3是低加疏水管路水封高度求解流程图。

图4是低加疏水管路注水示意图。

图5是低加疏水管路临界注水量求解流程图。

具体实施方式

下面结合附图和具体实施方式对本发明做进一步详细说明，以使本发明的优点和特征能更易于被本领域技术人员理解，从而对本发明的保护范围作出更为清楚明确的界定。

沿程阻力函数R推导流程：

式中：

ε——低加疏水管路内壁粗糙度，m；

d——低加疏水管路内径，m；

A——低加疏水管路内截面积，m²；

m——低加疏水管路阻力修正系数；

α——两相混合物的体积含气率；

x——两相混合物的质量含气率；

Q——疏水管路中流量，kg·s^-1；

ω_l——液相流速，m·s^-1；

Re_l——液相雷诺数；

λ_l——液相沿程阻力系数；

ρ_l——液相密度，kg·m^-3；

ρ_g——蒸汽密度，kg·m^-3；

μ_l——液相动力粘度，Pa·s；

μ_g——蒸汽动力粘度，Pa·s；

φ_l ²——摩擦倍增因子；

如图1所示，低压加热器的疏水通过疏水管路进入低加疏水立管。在设计工况下低加疏水在疏水管路中形成U型水封，防止低加抽气未被凝结直接流入凝汽器，造成热损失。在本发明中，如图2所示，输入低加汽侧压力p_a、凝汽器汽侧压力p_b、低加疏水流量Q_m、水封测量高度h_测，按照图2所示流程图，求解得到阻力修正系数m。

低加疏水管路阻力修正系数m求解流程：

1.输入低加汽侧压力p_a、凝汽器汽侧压力p_b、低加疏水流量Q_m、测量水封高度h_测；

2.将疏水管路按n个节点划分成n-1个微元段；

3.假设低加疏水管路阻力修正系数m；

4.假设低加疏水焓H_q，p₁＝p_a，i＝2，H₁＝H_q，Q₁＝Q_m；

5.根据节点焓值H_i-1、节点压力p_i-1、节点流量Q_i-1计算得到节点密度ρ_i-1、节点沿程阻力函数R_i-1；

6.计算节点压力

H_i＝H_q，Q_i＝Q_m；

7.根据H_i、p_i计算节点疏水干度x_i，若x_i＞0，x_i-1＜0，则该节点为汽液分界面，低加疏水管路水封高度h₁＝z_i；

8.计算下一个节点，即i＝i+1；

9.如果i＞n，则进行下一步，反之返回第5步；

10.若|p_n-p_b|≤0.001，则进行下一步，反之返回第4步重新假设低加疏水焓H_q；

11.若|h₁-h_测|≤0.001，则进行下一步，反之返回第3步重新假设低加疏水管路阻力修正系数m；

12.结束循环输出低加疏水管路阻力修正系数m、低加疏水焓H_q以及低加疏水管路不同位置处工质热力学状态。

通过上述计算求解获得低加疏水管路阻力修正系数之后，按照图3所示流程图，根据低加汽侧压力p_a、凝汽器汽侧压力p_b、低加疏水流量Q_m，可以获得整个疏水管路压力、干度分布，通过干度判断获得低加疏水管路水封高度。

图3是低加疏水管路水封高度求解流程图，计算水封高度步骤为：

1.输入低加汽侧压力p_a、凝汽器汽侧压力p_b、低加疏水流量Q_m，阻力修正系数m；

2.将疏水管路按n个节点划分成n-1个微元段；

3.假设低加疏水焓H_q，p₁＝p_a，i＝2，H₁＝H_q，Q₁＝Q_m；

4.根据节点焓值H_i-1、节点压力p_i-1、节点流量Q_i-1计算得到节点密度ρ_i-1、节点沿程阻力函数R_i-1；

5.计算节点压力

H_i＝H_q，Q_i＝Q_m；

6.根据H_i、p_i计算节点疏水干度x_i，若x_i＞0，x_i-1＜0，则该节点为汽液分界面，低加疏水管路水封高度h₁＝z_i；

7.计算下一个节点，即i＝i+1；

8.如果i＞n，则进行下一步，反之返回第4步；

9.若|p_n-p_b|≤0.001，则进行下一步，反之返回第3步重新假设低加疏水焓H_q；

10.结束循环输出低加疏水焓H_q、低加疏水管路水封高度h₁以及低加疏水管路不同位置处工质热力学状态。

通过计算不同工况下的水封高度，当计算水封高度大于疏水管路最低点位置时，则水封存在，当计算水封高度小于疏水管路最低点位置时，水封消失。通过分析水封消失时疏水管路不同位置处工质的热力学状态，对比低加疏水管路最低点压力与低加疏水焓值对应下的饱和压力，对水封消失的原因进行诊断。

根据上述计算，若水封消失，在U型水封底部注入冷却水，可以有效降低疏水焓值，降低汽化压力，低加疏水管路注水示意图如图4所示。根据低加汽侧压力p_a、凝汽器汽侧压力p_b以及低加疏水焓H_q，通过迭代计算求解维持U型水封的低加疏水管路临界注水量。图5是低加疏水管路补水流量求解流程图。

维持U型水封的低加疏水管路临界注水量的计算流程为：

1.将疏水管路按n个节点划分成n-1个微元段；

2.输入低加汽侧压力p_a、凝汽器汽侧压力p_b、低加疏水焓H_q、低加疏水管路阻力修正系数m、注水焓H_b，疏水管路最低点高度h_min；

3.假设低加疏水流量Q_m；

4.假设注水量Q_b，h₁＝0，p₁＝p_a，i＝2，H₁＝H_q，Q₁＝Q_m；

5.根据节点焓值H_i-1、节点压力p_i-1、节点流量Q_i-1计算得到节点密度ρ_i-1、节点沿程阻力函数R_i-1；

6.计算节点压力

7.疏水管路补水前，疏水焓H_i＝H_q，疏水流量Q_i＝Q_m；疏水管路补水后，疏水焓

疏水流量Q_i＝Q_m+Q_b；

8.根据H_i、p_i计算节点疏水干度x_i，若x_i＞0，x_i-1＜0，则该节点为汽液分界面，低加疏水管路水封高度h₁＝z_i；

9.计算下一个节点，即i＝i+1；

10.如果i＞n，则进行下一步，反之返回第5步；

11.若|p_n-p_b|≤0.001，则进行下一步，反之返回第4步重新假设注水量Q_b；

12.若|h₁-h_min|≤0.001，则进行下一步，反之返回第3步重新假设低加疏水流量Q_m；

13.结束循环输出维持U型水封的低加疏水管路临界注水量Q_b、低加疏水流量Q_m以及低加疏水管路不同位置处工质热力学状态。

本发明计算不同工况下的水封高度，当计算水封高度大于疏水管路最低点位置时，则水封存在；当计算水封高度小于疏水管路最低点位置时，水封消失。通过分析水封消失时疏水管路不同位置处工质的热力学状态，对比低加疏水管路最低点压力与低加疏水焓值对应下的饱和压力，对水封消失的原因进行诊断。当水封消失时，通过向疏水管路注入冷却水维持水封存在，计算维持U型水封的临界注水量。

Claims (1)

1.一种核电厂低加疏水管路U型水封消失及维持的诊断方法，其特征在于：包括如下步骤：

步骤1：根据核电厂低加疏水管路低加汽侧压力p_a、低加疏水流量Q_m、凝汽器汽侧压力p_b、水封测量高度h_测计算低加疏水管路阻力修正系数m；

整个低加 疏水管路的压降：

式中：

p_a—低加汽侧压力，Pa；

p_b—凝汽器汽侧压力，Pa；

ρ—疏水密度，是关于当地压力的函数，kg/m³；

g—当地重力加速度，m/s²

R—沿程阻力函数，Pa/m；

a—低加疏水管路起点，低压加热器；

b—低加疏水管路终点，低加疏水立管；

将疏水管路按空间离散，采用一阶迎风格式，则

p_i-1-p_i＝ρ_i-1g(z_i-z_i-1)+R_i-1(l_i-l_i-1)

整理得：

式中：

l_i—疏水管路第i个节点长度，m；

l_i-1—疏水管路第i-1个节点长度，m；

z_i—疏水管路第i个节点高度，m；

z_i-1—疏水管路第i-1个节点高度，m；

p_i—疏水管路第i个节点压力，pa；

p_i-1—疏水管路第i-1个节点压力，pa；

R_i-1—疏水管路第i-1个节点沿程阻力函数，Pa/m；

ρ_i-1—疏水管路第i-1个节点密度，kg/m³；

对低加的整个疏水管道进行建模，已知低加疏水管路的入口压力为低加汽侧压力p_a、出口压力为凝汽器汽侧压力p_b以及低加疏水流量Q_m；整个疏水管路的压降p_a-p_b受阻力修正系数m、低加疏水流量Q_m和低加疏水管路水封高度h₁影响，故

低加疏水管路阻力修正系数：

m＝f(Q_m,p_a,p_b,h₁)

式中：

Q_m—低加疏水流量，kg/s；

h₁—低加疏水管路水封高度，m；

低加疏水管路阻力修正系数m求解步骤为：

1.输入低加汽侧压力p_a、凝汽器汽侧压力p_b、低加疏水流量Q_m、测量水封高度h_测；

2.将疏水管路按n个节点划分成n-1个微元段；

3.假设低加疏水管路阻力修正系数m；

4.假设低加疏水焓H_q，p₁＝p_a，i＝2，H₁＝H_q，Q₁＝Q_m；

5.根据节点焓值H_i-1、节点压力p_i-1、节点流量Q_i-1计算得到节点密度ρ_i-1、节点沿程阻力函数R_i-1；

6.计算节点压力

H_i＝H_q，Q_i＝Q_m；

7.根据H_i、p_i计算节点疏水干度x_i，若x_i＞0，x_i-1＜0，则该节点为汽液分界面，低加疏水管路水封高度h₁＝z_i；

8.计算下一个节点，即i＝i+1；

9.如果i＞n，则进行下一步，反之返回第5步；

10.若|p_n-p_b|≤0.001，则进行下一步，反之返回第4步重新假设低加疏水焓H_q；

11.若|h₁-h_测|≤0.001，则进行下一步，反之返回第3步重新假设低加疏水管路阻力修正系数m；

12.结束循环输出低加疏水管路阻力修正系数m、低加疏水焓H_q以及低加疏水管路不同位置处工质热力学状态；

步骤2：根据低加疏水管路阻力修正系数m、低加汽侧压力p_a、低加疏水流量Q_m、凝汽器汽侧压力p_b计算低加疏水管路水封高度h₁；

通过上述计算求解获得低加疏水管路阻力修正系数m之后，根据低压加热器汽侧压力p_a、凝汽器汽侧压力p_b、低加疏水流量Q_m，获得整个疏水管路压力、干度分布，通过干度获得低加疏水管路水封高度h₁；

低加疏水管路水封高度：

h₁＝g(Q_m,m,p_a,p_b)

低加疏水管路水封高度h₁求解步骤为：

1.输入低加汽侧压力p_a、凝汽器汽侧压力p_b、低加疏水流量Q_m，阻力修正系数m；

2.将疏水管路按n个节点划分成n-1个微元段；

3.假设低加疏水焓H_q，p₁＝p_a，i＝2，H₁＝H_q，Q₁＝Q_m；

4.根据节点焓值H_i-1、节点压力p_i-1、节点流量Q_i-1计算得到节点密度ρ_i-1、节点沿程阻力函数R_i-1；

5.计算节点压力

H_i＝H_q，Q_i＝Q_m；

6.根据H_i、p_i计算节点疏水干度x_i，若x_i＞0，x_i-1＜0，则该节点为汽液分界面，低加疏水管路水封高度h₁＝z_i；

7.计算下一个节点，即i＝i+1；

8.如果i＞n，则进行下一步，反之返回第4步；

9.若|p_n-p_b|≤0.001，则进行下一步，反之返回第3步重新假设低加疏水焓H_q；

10.结束循环输出低加疏水焓H_q、低加疏水管路水封高度h₁以及低加疏水管路不同位置处工质热力学状态；

通过计算不同工况下的水封高度，当计算水封高度大于疏水管路最低点位置时，则水封存在，当计算水封高度小于疏水管路最低点位置时，水封消失；通过分析水封消失时疏水管路不同位置处工质的热力学状态，对比低加疏水管路最低点压力与低加疏水焓值对应下的饱和压力，对水封消失的原因进行诊断；

步骤3：根据低加疏水管路阻力修正系数m、低加汽侧压力p_a、低加疏水焓H_q、凝汽器汽侧压力p_b、注水焓H_b以及疏水管路最低点高度h_min，计算维持U型水封的低加疏水管路临界注水量Q_b进行计算；

在实际运行过程中，低压加热器换热不足时，低加疏水焓值增加，整个疏水管路的压力、干度发生变化，水封消失，低压加热器抽气直接进入凝汽器，潜热未被利用，造成热损失，使机组经济性降低；在核电厂中，通过在U型水封底部注入冷却水，能够有效降低疏水焓值，降低汽化压力； 根据低加汽侧压力p_a、凝汽器汽侧压力p_b、低加疏水焓H_q，通过迭代计算求解维持U型水封的低加疏水管路临界注水量；

维持U型水封的低加疏水管路临界注水量Q_b求解步骤为：

1.将疏水管路按n个节点划分成n-1个微元段；

2.输入低加汽侧压力p_a、凝汽器汽侧压力p_b、低加疏水焓H_q、低加疏水管路阻力修正系数m、注水焓H_b，疏水管路最低点高度h_min；

3.假设低加临界疏水量Q_m；

4.假设注水量Q_b，h₁＝0，p₁＝p_a，i＝2，H₁＝H_q，Q₁＝Q_m；

5.根据节点焓值H_i-1、节点压力p_i-1、节点流量Q_i-1计算得到节点密度ρ_i-1、节点沿程阻力函数R_i-1；

6.计算节点压力

7.疏水管路补水前，疏水焓H_i＝H_q，疏水流量Q_i＝Q_m；疏水管路补水后，疏水焓

疏水流量Q_i＝Q_m+Q_b；

8.根据H_i、p_i计算节点疏水干度x_i，若x_i＞0，x_i-1＜0，则该节点为汽液分界面，低加疏水管路水封高度h₁＝z_i；

9.计算下一个节点，即i＝i+1；

10.如果i＞n，则进行下一步，反之返回第5步；

11.若|p_n-p_b|≤0.001，则进行下一步，反之返回第4步重新假设注水量Q_b；

12.若|h₁-h_min|≤0.001，则进行下一步，反之返回第3步重新假设低加疏水流量Q_m；

13.结束循环输出维持U型水封的低加疏水管路临界注水量Q_b、低加疏水流量Q_m以及低加疏水管路不同位置处工质热力学状态。

Images