一种超超临界锅炉通用水动力计算方法与流程

本发明属于热能工程领域，涉及超超临界锅炉水冷壁系统的设计和安全性校核。特别涉及一种通用的可以使用到各种形式的超超临界锅炉水冷壁系统的设计和安全性校核方法。

背景技术：

超超临界锅炉水冷壁水动力特性研究是超超临界锅炉的关键技术之一，对于水冷壁优化设计和保证锅炉的安全可靠运行具有十分重要的意义。目前水动力计算方法主要包括弦截法计算模型、多项式拟合模型。弦截法采用弦截法求解非线性的动量守恒方程，为二阶收敛。从数学方法来看，因为需要求解方程组的偏导数，计算量较大，只能用于解决简单的回路问题。对工作在相同入口集箱和出口集箱之间的并联回路来说，多项式拟合方法是通过逐步计算，获得同一热负荷下各回路(包括节流元件、局部阻力元件、管子)在不同流量下的压降，得到这些流量下的压降－流量关系。然后，根据已经得到的压降－流量关系，采用多项式拟合的方法，得到压降随流量变化的计算多项式，也即得到了这些并联回路的水动力曲线。根据并联回路“压降相等，流量相加”的原则，可以得到并联管屏的总水动力曲线，并将并联管屏用一根当量管来代替。按照上述方法，水冷壁总是可以简化为由串联回路组成的管路系统，按照串联回路“流量相等，压降相加”的原则，可以得到水冷壁总的水动力曲线，即水冷壁系统总压降与进口流量的关系。最后，采用与上述计算相反的顺序，根据水冷壁入口流量已知这一条件，就可先后求出水冷壁总压降、各串联回路压降、并联管屏中各回路压降与流量等参数。在对并联管屏各回路的水动力曲线进行拟合时，一般采用具有较高精度的切比雪夫方法。但是这种多项式拟合计算模型，只能对简单的串并联回路组合进行计算。在复杂的水冷壁结构时，需要对其进行简化，大大的限制了模型的适用范围，也影响了模型的精度。

为了克服复杂超超临界锅炉的水动力计算的困难，不少学者使用计算机编写过水动力计算程序但是大多针对某一种特定的超超临界锅炉。也有学者提出了“一种锅炉水动力通用设计和校核方法”，这种方法采用部件压力法作为水动力的基本方法，该方法忽略了锅炉水冷壁吸热过程中工质热物性变化对阻力的影响且水冷壁壁温计算。随着现在超超临界锅炉参数的运行参数的提高，水冷壁出口汽温偏差及壁温的安全裕度变的变的越来越小，需要一种可应用于计算机且全方位考虑锅炉水冷壁系统内工质状态变化、系统质量守恒、系统能量守恒和系统动量守恒的通用超超临界锅炉水动力计算方法。

技术实现要素：

本发明的目是克服上述现有方法的不足，提出了在计算机上应用的非线性模型组成的流动网络系统法。本计算方法适用于各种锅炉整体布置形式(П型、塔式、T型)、各种管圈型式(螺旋管圈水冷壁、垂直水冷壁)、各种燃烧技术(切向燃烧方式、墙式燃烧方式、W火焰方式、循环流化床燃烧方式)、各种质量流速组合方式下的超超临界直流锅炉水动力计算和壁温安全分析。

为了实现上述目的，本发明的技术方案是：

一种超超临界锅炉通用水动力计算方法，该方法将水冷壁简化为由连接管、受热回路和节点等元件组成的流动网络系统，根据质量守恒方程、动量守恒方程和能量守恒方程建立流量分配的非线性计算模型。通过对非线性方程组的求解，确定各受热回路和连接管中的流量分配及节点压力分布，在此基础上，进行锅炉水冷壁流动和壁温安全校验，指导锅炉设计。包括下列步骤：

1)划分回路与管段：

回路划分：将相邻的几何结构和吸热偏差基本相同的水冷壁管子划分为同一个流动回路。在炉膛水平方向吸热偏差变化明显的或几何结构变化较大的部位，回路划分应密集，对应每回路分配管子根数较小，而在吸热偏差或几何结构变化平缓的区域，回路划分可以较为稀疏，对应每个回路分配的管子根数较多。

管段划分：将每一个回路沿工质流动方向划分为若干个管段。在炉膛热负荷变化剧烈或物性变化较大的部位管段划分应该密集，对应管段长度较短，而在远离煤粉燃烧器或循环流化床锅炉布风板等区域管段可以划分较为稀疏，对应管段的较长。

2)确定热负荷分布：

根据典型炉型的热负荷分布曲线对应拟合新锅炉热负荷曲线，可以分别得到水平环带热负荷q_avg(z)和尖峰热负荷q_pe仅关于炉膛高度的函数关系式。对任意受热回路的某一管段来说，取管段中心高度位置处的高度值z_m代入热负荷函数关系式便可计算得到该管段热负荷。其中水平环带热负荷用于吸热量和流量分配的计算，尖峰热负荷用于壁温的计算和膜态沸腾的判别。

3)简化水冷壁系统：

根据步骤1)划分回路的方法，将整个水冷壁简化为由连接管回路、受热回路和节点元件组成的流动网络系统，并对应各个回路和节点进行编号处理。简化的流动网络系统中，连接管回路代表水冷壁系统中的不受热的连接管，受热回路代表水冷壁系统中受热的水冷壁管子或者受热的悬吊屏，节点代表水冷壁系统中各种集箱和汽水分离器。

4)建立非线性模型：

根据质量守恒方程、动量守恒方程和能量守恒方程建立以回路单管流量和集箱压力为未知量的非线性方程组计算模型。质量守恒方程即为流入节点的流量等于流出节点的流量；动量守恒方程即回路各管段的压降总和等于该回路进出口节点压力差，管段压降为本管段内流动阻力、重位压降和局部阻力的代数和，各管段压降的代数和即管子压降；能量守恒方程即各管段吸收热量等于工质焓值增加量，能量守恒用于计算工质状态。

计算回路各管段压降时，首先要判断工质压力状态，超临界压力下，根据进出口焓值相对大小确定超临界水和蒸汽段长度分别计算其压降。临界压力以下，根据管段进口和出口的工质状态，分为五种不同的情况进行压降计算，五种情况分别是：第一种：进口过冷，出口过冷，第二种：进口过冷，出口两相。第三种：进口两相，出口两相，第四种：进口两相，出口过热；第五种：进口过热、出口过热；进出口含有相变点的管段，根据进出口焓值相对于饱和水焓值和饱和蒸汽焓值的大小可以得到两相工质长度和单相工质长度，然后分别计算管段在相变点前部分和相变点后部分的压降。

计算ΔP_i过程中局部压降、重位压降及加速压降参照国家水动力计算标准压降计算的内容；在超超临界压力下的内螺纹管单相区摩擦阻力系数按照下式计算：

5)非线性方程组求解：

通过弦割法求解非线性方程组，获得各受热回路和连接管中的流量分配及节点压力分布值。

6)壁温计算：

首先计算出各回路管段中心位置处的工质压力，然后按超临界压力和亚临界压力区的传热规律分别计算内壁与外壁温度、中间壁温以及鳍片温度；计算壁温需要分管型计算，管型包括光管与内螺纹管、垂直管圈与螺旋管圈；在亚临界压力区，管段中心位置处工质为汽水两相时，需要判断该管段是否发生膜态沸腾和干涸，并根据不同区域的传热规律进行壁温计算；超临界压力区的计算首先判断管段中心流体为超临界水或超临界蒸汽，然后判断是否有发生拟膜态沸腾的可能，并根据工质状态、管型和传热规律的不同组合计算壁温。

锅炉在高负荷运行时，工质水处于超临界区。在低负荷工作时，随着吸热的进行会经历单相水，汽水两相，及单相汽等几个阶段。

在压力大于22.115MPa的超临界区，换热公式按如下公式计算。

低焓值区：

高焓值区：

式中：Nu表示努塞尔数；Re为雷诺数；H为单位流体焓值，kJ/kg；Pr为普朗克数；v为工质比容，kJ/(kg℃)；小标w和f表示壁面及工质特性。

在临界压力附近，流动传热过程比较复杂。流动传热过程中可能发生传热恶化，传热恶化后的传热计算如下。

在压力为19～22.1MPa时，计算临界质量流速如下:

G_cr＝800.45+223.87ln(22.115-p) (4)

当G<G_cr时，计算临界热负荷如下：

q_cr＝3343.93(22.115-p)^0.4091×G^-0.3835(1-x)^0.6792 (5)

当G>G_cr时，计算临界热负荷如下：

q_cr＝2.2669(22.115-p)^0.1007×G^0.7385(1-x)^0.1888 (6)

式中：q为热流密度，kW/m2；p为压力，MPa；G为质量流速，kg/(m2s)；cr表示临界特性。

当管子内壁热负荷大于临界热负荷时，DNB发生，此时的传热关联式为：

式中：λ_cr为热力学临界点导热系数，值为0.914；p_cr为水的临界压力，值为22.115MPa；L和G分别表示液相和汽相。

换热计算模型对水冷壁金属壁温影响非常大，其他压力及情况下的传热按照《电站锅炉水动力计算方法》(JB/Z 201-83)计算。

管子正面内壁温度计算公式如下：

式中：t_f为壁温计算点处管内工质温度，℃；J_n为管子正面内壁热量均流系数；β为管子外径与内径的比值，β＝d_w/d_i；q_w为壁温计算点处正面外壁辐射热负荷，kW/m2；α₂为壁温计算点管子内壁与介质间的换热系数，kW/(m2℃)。

管子正面外壁温度计算公式如下：

式中：为管子正面沿厚度方向上的平均热量均流系数；δ为管壁厚度，m；λ为管子部分的金属导热系数，W/(m℃)。

鳍根温度计算公式如下：

式中：t_f为壁温计算点处管内工质温度，℃；q_wmax为壁温计算点处局部最大内壁辐射热负荷，kW/m2；μ_qg为鳍根处热流密度分流系数；R为外径，m；r为内径，m。

鳍端温度计算公式如下：

式中：μ_qd为鳍端处热流密度系数；s为相邻管之间的节距，m。

本发明的技术效果：

利用本发明提出的超超临界锅炉通用水动力计算方法，对各种锅炉整体布置形式(П型、塔式、T型)、各种管圈型式(螺旋管圈水冷壁、垂直水冷壁)、各种燃烧技术(切向燃烧方式、墙式燃烧方式、W火焰方式、循环流化床燃烧方式)、各种质量流速组合方式下的的超超临界直流锅炉进行水动力计算，可以获得锅炉在各种负荷下的水冷壁压降、水冷壁流量分配情况、炉膛出口汽温分布情况及水冷壁金属壁温。目前塘寨电厂600MW超超临界W型火焰锅炉现场能够采集到锅炉炉膛出口汽温分布的运行DCS(Distributed Control System)界面数据，将计算值与实际DCS数据值进行比较，对比结果如图13和图14所示，可以发现计算与实测数据能够很好地吻合。图13为600MW超超临界W型火焰锅炉BMCR负荷时下炉膛出口汽温实际值与计算值对比，图14为600MW超超临界W型火焰锅炉BMCR负荷时上炉膛出口汽实际值与计算值对比。

附图说明

图1为热负荷计算示意图；

图2为倾斜炉膛热负荷修正系数选取图；

图3为四角切圆煤粉锅炉水平环带热负荷曲线；

图4为切圆为顺时针方向时的水平方向吸热不均匀系数；

图5为W火焰锅炉水平环带热负荷曲线；

图6为W火焰锅炉水平方向吸热不均匀系数；

图7为循环流化床锅炉水平环带热负荷曲线；

图8为循环流化床锅炉水平方向吸热不均匀系数；

图9为简化的流动网络系统示意图；

图10为流动网络法应用到计算机上的计算流程图；

图11为单管压降计算流程图；

图12为壁温计算流程图；

图13为600MW超超临界W型火焰锅炉BMCR负荷时下炉膛出口汽温实际值与计算值对比图；

图14为600MW超超临界W型火焰锅炉BMCR负荷时上炉膛出口汽实温际值与计算值对比图。

为了更加清楚地理解本发明，以下结合附图和发明人给出的具体实例对本发明进行进一步详细描述。

具体实施方式

1、划分回路与管段：

回路划分：组成锅炉水冷壁的各根管子几何结构和吸热情况都不相同，存在着结构偏差和热力偏差。将相邻的几何结构和吸热偏差基本相同的管子划分在同一个回路，认为其中各单管的热力参数、几何结构完全相同。同时在炉膛水平方向吸热偏差变化明显的或几何结构变化较大的部位，回路划分应密集，对应每回路分配管子根数较小，而在吸热偏差或几何结构变化平缓的区域，回路划分可以较为稀疏，对应每个回路分配管子根数较多。

管段划分：为准确计算压力、温度及焓值沿管长的变化特性将每一个回路沿工质流动方向划分为若干个管段，认为在每一个管段中，认为各热力参数保持恒定不变。同时在在炉膛热负荷变化剧烈或物性变化较大的部位管段划分应该密集，对应管段长度较短，而在远离煤粉燃烧器或循环流化床锅炉布风板的区域管段可以划分的较长。

2、确定水平环带热负荷和尖峰热负荷：

根据典型炉型的热负荷分布曲线对应拟合新锅炉热负荷曲线，可以分别得到水平环带热负荷q_avg(z)和尖峰热负荷q_pe仅关于炉膛高度的函数关系式。对任意受热回路的某一管段来说，取管段中心高度位置z_m处的高度值代入热负荷函数关系式便可计算得到该管段热负荷。其中水平环带热负荷用于吸热量和流量分配的计算，尖峰热负荷用于壁温计算和膜态沸腾判别。

局部热负荷q_h(x,z)与尖峰热负荷q_p(z)的计算示意图见图1。

局部热负荷q_h(x,z)的计算公式为：

q_h(x,z)＝f_sp×q_avg(z)×M(x) (12)

式中：x—炉膛水平方向坐标，m；

z—炉膛高度方向坐标，m；

f_sp—倾斜炉膛相对于垂直炉膛的热负荷修正系数；

q_avg(z)—炉膛水平环带热负荷，kW/m²；

M(x)—炉膛水平方向吸热不均匀系数。

尖峰热负荷q_p(z)的计算公式为：

q_p(z)＝f_sp×q_pe(z) (13)

式中：q_pe—垂直炉膛尖峰热负荷，kW/m²。

q_pe也只和高度有关，取管段出口位置对应高度z₂处的值。

M(x)和水平位置有关，取管段中心水平位置x_m处对应的值。

对于垂直炉膛，f_sp＝1，П型锅炉倾斜炉膛的热负荷修正系数如图2所示。

其中四角切圆煤粉锅炉水平环带热负荷曲线如图3所示；切圆为顺时针方向时的水平方向吸热不均匀系数如图4所示；图5为W火焰锅炉水平环带热负荷曲线；图6为W火焰锅炉水平方向吸热不均匀系数；图7为循环流化床锅炉水平环带热负荷曲线；图8为循环流化床锅炉水平方向吸热不均匀系数。利用这些不同燃烧方式锅炉热负荷分布规律，该水动力计算方法可以准确计算多种锅炉水动力特性。

3、简化水冷壁系统：

根据步骤1)划分回路的方法，将整个水冷壁简化为由连接管回路、受热回路和节点元件组成的流动网络系统，并对应各个回路和节点进行编号处理。简化的流动网络系统中，连接管回路代表水冷壁系统中的不受热的连接管，受热回路代表水冷壁系统中受热的水冷壁管子或者受热的悬吊屏，节点代表水冷壁系统中各种集箱和汽水分离器。图9示出了一个简单的流动网络系统示意图。图中编号1-12为受热回路，13-15为连接管回路，16-20指的是集箱压力节点。图中箭头方向代表工质流动方向。

4、建立非线性模型：根据质量守恒方程、动量守恒方程和能量守恒方程建立流量分配和压力节点的非线性方程组计算模型。质量守恒方程即为流入节点(集箱)的流量等于流出节点的流量。动量守恒方程即，回路各管段的压降总和等于该回路进出口节点压力差。

能量守恒方程即各管段吸收热量等于工质焓值增加量，能量守恒用于计算工质状态。受热回路的能量守恒方程形式为:

式中：h_out—回路管段出口焓值，kJ/kg；h_in—回路管段进口焓值，kJ/kg；s—管子节距，m；l—管子长度，m；w—管子质量流速，kg/s。

以图9的流动网络系统为例，建立流动网络系统的守恒方程、动量守恒方程。

对于受热回路1-12，忽略集箱中的静压变化，可以列出12个动量守恒方程：

[P_in-x(17)]-ΔP_i(x(i))＝0 i＝1,2 (15)

[P_in-x(18)]-ΔP_i(x(i))＝0 i＝3,4 (16)

[x_in-x_out]-ΔP_i(x(i))＝0 i＝5,12 (17)

对于连接管13-15，可以列出3个动量守恒方程：

[P_in-x(18)]-ΔP_i(x(i))＝0 i＝13 (18)

[x(18)-x(19)]-ΔP_i(x(i))＝0 i＝14,15 (19)

式中：i—回路或节点编号，一个回路中包含有多个受热和结构相同的水冷壁单管。

x(i)—i回路中的单管流量或i节点压力，kg/s或Pa；

x_in—回路5-12进口处集箱压力，Pa；

x_out—回路5-12出口处集箱压力，Pa；

ΔP_i—i回路压降，包括摩擦压降、局部压降、重位压降及加速压降，Pa。

计算ΔP_i过程中局部压降、重位压降及加速压降参照国家水动力计算标准《电站锅炉水动力计算方法》(JB/Z 201-83)的内容；在超超临界压力下的内螺纹管单相区摩擦阻力系数按照下式计算：

节点16-20遵守的质量守恒方程为：

[x(1)+x(2)]-[x(5)+x(6)+x(9)+x(10)]＝0 i＝16 (21)

[x(3)+x(4)]-[x(7)+x(8)+x(11)+x(12)]＝0 i＝17 (22)

[x(13)+x(9)+x(10)+x(11)+x(12)]-[x(14)+x(15)]＝0 i＝18 (23)

G-[x(5)+x(6)+x(7)+x(8)+x(14)+x(15)]＝0 i＝19 (24)

以上12个受热回路和3个连接管遵守的15个动量守恒方程，以及4个节点遵守的4个质量守恒方程共19个方程组成了一个的封闭的非线性方程组。该非线性方程组包括15个流量未知数x(1),x(2),...,x(15)和4个压力未知数x(16),x(17),x(18),x(19)，共19个未知数。在进口集箱压力P_in和进口流量G已知的条件下，对非线性方程组进行数值迭代求解即可获得各回路、连接管的流量分配，以及各节点的压力。

在进行锅炉水动力计算时，流量G、进口集箱焓值h_in以及出口集箱焓值h₁₉通过热力计算已知。如果进口集箱压力P_in未知，而出口集箱压力P₁₉已知，则按上述方法同样可以列出包含P_in的封闭非线性方程组，从而求出各未知参数。

5)通过弦割发求解非线性方程，获得各受热回路和连接管中的流量分配及节点压力分布。如图10所示流动网络法应用到计算机上的计算流程图，图中HTEF—热平衡系数；k—迭代次数；i—回路编号；j—管段编号；I_cir—回路总数；J_cir—某一回路的管段总数；n—回路与节点总数；f_i—第i个非线性方程的残差；ε—动量守恒方程收敛准则；δ—质量守恒方程收敛准则；—第k+1次迭代时第i个非线性方程的迭代新解。整个计算对所建立的由未知的单管流量和集箱节点压力组成的非线性方程组进行迭代求解。在得到流量分配和集箱压力分布后进行壁温的求解计算。

6)壁温计算：如图12，壁温计算需要计算出管段中心位置处的工质压力，然后按超临界压力和亚临界压力区的传热规律分别计算内壁与外壁温度、中间壁温以及鳍片温度。计算壁温需要分管型计算，管型包括光管与内螺纹管、垂直管圈与螺旋管圈。在亚临界压力区，管段中心位置处工质为汽水两相时，需要判断该管段是否发生膜态沸腾和干涸，并根据不同区域的传热规律进行壁温计算。超临界压力区的计算首先判断管段中心流体为超临界水或超临界蒸汽，然后判断是否有发生拟膜态沸腾的可能，并根据工质状态、管型和传热规律的不同组合计算壁温。

锅炉在高负荷运行时，工质水处于超临界区。在低负荷工作时，随着吸热的进行会经历单相水，汽水两相，及单相汽等几个阶段。

在压力大于22.115MPa的超临界区，换热公式按如下公式计算。

低焓值区：

高焓值区：

式中：Nu表示努塞尔数；Re为雷诺数；H为单位流体焓值，kJ/kg；Pr为普朗克数；v为工质比容，kJ/(kg℃)；小标w和f表示壁面及工质特性。

在临界压力附近，流动传热过程比较复杂。流动传热过程中可能发生传热恶化，传热恶化后的传热计算如下。

在压力为19～22.1MPa时，计算临界质量流速如下:

G_cr＝800.45+223.87ln(22.115-p) (27)

当G<G_cr时，计算临界热负荷如下：

q_cr＝3343.93(22.115-p)^0.4091×G^-0.3835(1-x)^0.6792 (28)

当G>G_cr时，计算临界热负荷如下：

q_cr＝2.2669(22.115-p)^0.1007×G^0.7385(1-x)^0.1888 (29)

式中：q为热流密度，kW/m²；p为压力，MPa；G为质量流速，kg/(m²s)；cr表示临界特性。

当管子内壁热负荷大于临界热负荷时，DNB发生，此时的传热关联式为：

式中：λ_cr为热力学临界点导热系数，值为0.914；p_cr为水的临界压力，值为22.115MPa；L和G分别表示液相和汽相。

换热计算模型对水冷壁金属壁温影响非常大，其他压力及情况下的传热按照《电站锅炉水动力计算方法》(JB/Z 201-83)计算。

管子正面内壁温度计算公式如下：

式中：t_f为壁温计算点处管内工质温度，℃；J_n为管子正面内壁热量均流系数；β为管子外径与内径的比值，β＝d_w/d_i；q_w为壁温计算点处正面外壁辐射热负荷，kW/m²；α₂为壁温计算点管子内壁与介质间的换热系数，kW/(m²℃)。

管子正面外壁温度计算公式如下：

式中：为管子正面沿厚度方向上的平均热量均流系数；δ为管壁厚度，m；λ为管子部分的金属导热系数，W/(m℃)。

鳍根温度计算公式如下：

式中：t_f为壁温计算点处管内工质温度，℃；q_wmax为壁温计算点处局部最大内壁辐射热负荷，kW/m²；μ_qg为鳍根处热流密度分流系数；R为外径，m；r为内径，m。

鳍端温度计算公式如下：

式中：μ_qd为鳍端处热流密度系数；s为相邻管之间的节距，m。

超超临界锅炉通用水动力计算方法对适合对各种容量等级、各种蒸汽参数、各种炉型和燃烧技术的超超临界锅炉和汽包锅炉的水动力特性和壁温特性进行计算分析。能够计算出各种负荷工况下的水冷壁回路流量分配、压力沿流动方向的变化、出口汽温分布、工质温度与管壁金属温度沿高度方向的变化趋势。对水冷壁在各种负荷工况条件下的壁温安全特性进行校核计算，为锅炉设计和运行提供依据。图10示出了超超临界锅炉通用水动力计算方法应用的计算机上的计算流程。

单管压降的流程如图11。图中：m—管段编号；m＝m+1—表示m的新值等于m的当前值加1；P—管段压力，Pa；P_crit—临界压力，Pa；M_crit—回路最大管段数。单管压降计算过程与工质的状态有关，单管压降的计算程序编制遵循下列原则：

(1)按管子结构和受热状态不同划分管段，管段越多计算精度越高。

(2)根据管段进口和出口的工质状态，分为五种不同的情况进行压降计算。含有相变点的管段需首先确定相变点的位置，然后分别计算管段在相变点前部分和相变点后部分的压降。

(3)管段压降为本管段内流动阻力、重位压降和局部阻力的代数和，各管段压降的代数和即管子压降。

(4)管段的出口参数(压力、流量、温度、焓值)也是后一管段的入口参数，最后一个管段的出口参数即是管子出口参数。

壁温计算流程如图12。图中x—管段中心位置处的工质干度；—光管发生干涸时的临界干度；—内螺纹管发生干涸时的临界干度；—光管发生膜态沸腾时的临界热负荷；—内螺纹管发生膜态沸腾时的临界热负荷，壁温计算程序编制遵循下列原则：

(1)壁温计算模块需要计算出管段中心位置处的工质压力，然后按超临界压力和亚临界压力区的传热规律分别计算内壁与外壁温度、中间壁温以及鳍片温度。

(2)管型包括光管与内螺纹管、垂直管圈与螺旋管圈。

(3)在亚临界压力区，管段中心位置处工质为汽水两相时，需要判断该管段是否发生膜态沸腾和干涸，并根据不同区域的传热规律进行壁温计算。

(4)超临界压力区的计算首先判断管段中心流体为超临界水或超临界蒸汽，然后判断是否有发生拟膜态沸腾的可能，并根据工质状态、管型和传热规律的不同组合计算壁温。