周涛 李兵 齐实 马栋梁 黄彦平



模糊贴近度方法在超临界水堆换热关系式选择中应用

周涛1,2,3李兵1,2,3齐实1,2,3马栋梁1,2,3黄彦平4

1（华北电力大学核科学与工程学院 北京 102206）2（华北电力大学核热工安全与标准化研究所 北京 102206）3（华北电力大学非能动核能安全技术北京市重点实验室 北京 102206） 4（中国核工业集团公司核反应堆热工水力技术重点实验室 成都 610041）

以超临界水堆换热关系式为基础，利用模糊数学中贴近度的概念，并引入变异系数法确定各影响因素的权重。比较了确定工况与各经验公式适用范围的贴近程度，获得了该工况下最为适用的经验关系式。结果表明，采用Hamming贴近度、Euclid贴近度及最大最小贴近度等三种方法的计算结果基本一致，模糊数学中贴近度的应用，为超临界水堆的流动换热问题提供了一种可以在已知关系式的基础上，做出针对特定工况的最优关联式选择方法，使计算结果更加精确。

超临界水堆，关联式，模糊贴近度，变异系数法

超临界水冷堆(Supercritical Water-cooled Reactor, SCWR)作为第四代反应堆中唯一的水冷堆，堆芯工作压力25MPa，出口冷却剂温度500ºC左右，具有经济性好、系统简单和可持续性强等诸多优势[1‒2]。SCWR的热工水力过程是研究超临界水堆的关键环节，直接关系到反应堆安全性和经济性。随着超临界水传热和流动研究的发展，国内外诸多学者在超临界水传热关系式方面开展了大量工作[3‒8]。Bishop等[3]研究了近临界温度和超临界压力下的强迫对流换热关系式；Watts等[4]研究了竖直通道内超临界水核热耦合瞬态工况传热关系式；陈玮玮等[8]通过对大量实验数据的分类对比对超临界水竖直管内传热关系式计算的精度进行了研究。这些关系式在各自的实验条件下与实验符合良好，但是，随着工况的变化，出现了不同程度的误差。因此，针对实际工况选择合适的关系式，可以提高实际计算的准确度，减小计算误差。在影响换热的因素中，每一种因素会因其所处的条件差异表现出不同的特征，其中的关系较为复杂，不易确定，具有一定的模糊性。换热关系式的选择是依据各影响因素的标准，判断实际工况与哪个换热关系式标准最接近，则认为该工况最符合这个关系式。模糊贴近度是模糊数学中的一种比较实用的方法[9]，用来描述两个模糊集的贴近程度。它已成功地在水环境质量评价、土壤环境评价中得到广泛应用[10‒11]。以现有存在的传热关系式为基础，借助模糊数学中贴近度的概念，通过比较与各个经验公式适用范围的贴近程度，找到了实际工况下最为适用的经验关系式，可以减少计算误差。

1 计算方法

1.1 贴近度公式

贴近度[9]是用来度量两个模糊集接近程度的数量指标，式(1)‒(6)给出其定义。

Hamming贴近度：

Euclid贴近度：

最大最小贴近度：

1.2 归一化公式

为了更好地比较不同取值范围间向量的差异，计算前通过式(7)‒(8)对数据进行归一化，使其位于[0, 1]之间。

1.3 权重系数公式

由于每个影响因素对结果影响大小不同，需要考虑各影响因素的权重问题。权重系数的确定方法一般有层次分析法和变异系数法等，研究采用变异系数法来确定影响因素的权重指标[12‒13]。变异系数法认为某个影响因素数值之间的差异性越大，则该指标在综合评价中的作用越大，其权重也越大。所以，根据各个影响因素数值之间的差异性，利用变异系数这个参数来确定各影响因素的权重。变异系数法计算权重如式(9)、(10)所示。

1.4 选择原则

2 选择对象

表1 各关系式适用范围

3 选择流程

基于§1.2中的选择原则，进行换热关系式选择，计算流程如下：

1) 建立评价标准及整理样本：根据选择的不同换热关系式的适用范围，确定各关系式的特征指标。

2) 样本数据归一化：使计算结果方便、准确，对样本进行归一化处理。

3) 权重系数矩阵确定：压力(, MPa)、质量流密度(, kg·m−2·s−1)等影响因素对换热系数(, W·m−2·k−1)的影响程度不同，采用变异系数法进行权重计算。

4) 计算模糊贴近度：分别采用Hamming贴近度、Euclid贴近度、最大最小贴近度计算样本与各个关系式的贴近程度，对结果进行对比。

5) 确定待分类样本及验证：根据模糊择近原则，选择与样本贴近度最大的关系式进行计算。

4 选择计算及验证

4.1 评价标准及样本

根据现有的研究，超临界水换热系数的计算主要与压力、质量流密度、热流密度()有关，选取这三个参量作为考核标准，根据贴近度定义[9]描述某一点与某范围的贴近程度，则认为该点与区间中点越靠近，则越贴近区间；集本团队部分研究[14‒15]成果，选取表1中参数中点作为换热关联式的参数指标如表2所示。另外方便检验，选取Swenson等[7]的实验数据作为待选定数据。

表2 超临界水换热公式评价指标及待选定参数

4.2 数据归一化计算

为了更好地比较不同取值范围间向量的差异，计算前采用式(8)对表2中数据进行归一化，使其位于[0, 1]之间。计算结果如表3所示。表2与表3意义相同。

表3 归一化后超临界水换热公式评价指标及待选定参数

4.3 权重系数计算

由于每个因素影响换热系数的程度不相同，采用式(9)、(10)计算表3中的数据得到各影响因素权重：压力为0.1，质量流密度为0.4，热流密度为0.5。

4.4 贴近度结果计算

根据表3中的数据及各影响因素权重，分别采用式(4)‒(6)计算Hamming贴近度、Euclid贴近度、最大最小贴近度结果如表4所示。

表4 不同贴近度计算值

从表4可知，采用Hamming贴近度、Euclid贴近度计算结果一致，由好到差的顺序为：Swenson、Yamagata、Bishop、Shitsman和Watts公式。最大最小贴近度计算结果由好到差的顺序为：Swenson、Bishop、Yamagata、Shitsman和Watts公式；但Bishop和Yamagata公式计算结果大小接近。

引入平均绝对误差(Mean Absolute Deviation, MAD)MA、平均相对误差(Mean Relative Deviation, MRD)MR、标准偏差(Standard Deviation, SD)S来评价Swenson换热公式。

式中：MA表示公式计算值与实验值之间差距的大小；S表示公式计算值偏离算数平均值程度的大小，S越小，表示这些数值偏离平均值的程度越小，反之则偏离平均值的程度越大。Swenson换热公式的平均绝对相对误差MA12.8%、标准偏差S3.01%。

4.5 计算结果验证

通过表1中的换热关系式计算，得出换热系数随主流温度的变化情况，并与Swenson实验数据[7]进行比较，结果如图1所示。

图1 各关系式计算的换热系数随主流温度T的变化情况与实验值的比较

从图1中可知，以上5个换热关系式都可以预测出超临界水换热系数随着温度升高，出现先增大后减小的趋势。其中，Shitsman公式预测的换热系数比实验值高将近两倍。Watts公式预测结果偏低，且最大值出现位置提前。Bishop、Yamagata、Swenson公式预测效果较好。因此，在单一选用一种换热关联式进行计算时，在部分工况下与实验数据相比会出现较大偏差。贴近度计算就是针对不同工况，合理选择最符合的超临界换热关系式。

5 结语

1) 借助模糊贴近度的概念，并引入变异系数法综合考虑各影响因素权重系数，以具体的工况为例采用Hamming贴近度、Euclid贴近度及最大最小贴近度三种贴近度对该工况计算，比较分析了超临界水换热公式选择的步骤，获得了较好的效果。

2) 模糊贴近度综合考虑了各因素对超临界水换热系数的影响，根据确定工况与各公式的贴近程度进行选择。该方法评价原理直观、计算方法简单，是一种实用而准确的评价方法，为超临界水堆的流动换热问题提供了一种可以在已知关联式的基础上，做出针对实际工况的最优关联式选择方法，使实际计算结果更加准确。

3) 模糊数学中贴近度在超临界水堆中的应用是一种新的途径，随着对流动换热问题的深入研究，实验关联式增多，这种方法越精确，所做的选择与实际越相符。模糊评价指标选定是后续工作研究的重点。下一步打算采用权值的办法进行修正即选择中点、左四分之一点及右四分之三点，针对不同的情况，进行不同的权值修正，以得到更为准确的模糊评价指标。对于研究适用范围广的换热关系式，该方法还存在一些不足，改善模糊评价指标、提高结果精确度是后续工作研究的主要问题。

1 程旭, 刘晓晶. 超临界水冷堆国内外研发现状与趋 势[J]. 原子能科学技术, 2008, 42(2): 167‒172. CHENG Xu, LIU Xiaojing. Development status and trend of supercritical water cooled reactor at home and abroad[J]. Atomic Energy Science and Technology, 2008, 42(2): 167‒172.

2 Oka Y, Koshizuka S, Ishiwatari Y. Super light water reactors and super fast reactors: supercritical-pressure light water cooled reactors[M]. Springer, 2010.

3 Bishop A A, Sandberg R O, Tong L S. Forced convection heat transfer to water at near-critical temperatures and supercritical pressures[R]. Report WCAP-2056, Part IV, Pittsburgh, USA: Westinghouse Electric Corp, 1964.

4 Watts M J, Chou C T. Mixed convection heat transfer to supercritical pressure water[C]. Proceedings of the 7th IHTC, Munchen, Germany, 1982: 495‒500.

5 Yamagata K, Nishikawa K, Hasegawa S,. Forced convective heat transfer to supercritical water flowing in tubes[J]. International Journal of Heat and Mass Transfer, 1972, 15(12): 2575‒2593.

6 Shitsman M E. Impairment of the heat transmission at supercritical pressures[J]. High Temperature, 1963, 1: 237‒244.

7 Swenson H S, Carver J R, Kakarala C R. Heat transfer to supercritical water in smooth-bore tubes[J]. Journal of Heat Transfer, 1965, 87(4): 477‒483.

8 陈玮玮, 方贤德, 商辉, 等. 超临界压力下竖直管内水的传热关系式研究[J]. 工程热物理学报, 2016, 37(1): 104‒110. CHEN Weiwei, FANG Xiande, SHANG Hui,. Study on heat transfer formula of vertical pipe water under supercritical pressure[J]. Journal of Engineering Thermophysics, 2016, 37(1): 104‒110.

9 李安贵, 张志宏, 孟艳, 等. 模糊数学及其应用[M]. 2版. 北京: 冶金工业出版社, 2006: 31‒44. LI Angui, ZHANG Zhihong, MENG Yan,. Fuzzy mathematics and its application[M]. 2nd ed. Beijing: Metallurgical Industry Press, 2006: 31‒44.

10 Guan Q C, Zhao S X, Liu J F,. Application of fuzzy nearness method in groundwater quality evaluation[J]. Water Saving Irrigation, 2012, 7: 43‒46.

11 Wang J G, Yang L Z. Application of fuzzy mathematics to soil quality evaluation[J]. Acta Pedologica Sinica, 2001, 2: 176‒183.

12 Shi G X, Jiang Y X, Chen X S. Evaluation model of water environment quality assessment and its application based on coefficient of variability[J]. Ground Water, 2002, 34(3): 147‒148.

13 邵良杉, 赵琳琳, 温廷新, 等. 基于区间直觉模糊数的双向投影决策模型[J]. 计算机工程与应用, 2017, 53(1): 83‒86. SHAO Liangshan, ZHAO Linlin, WEN Tingxin,. Bi directional projection decision model based on interval valued intuitionistic fuzzy numbers[J]. Computer Engineering and Applications, 2017, 53(1): 83‒86.

14 胡雨, 周涛. 贴近度在超临界水堆实验关联式选择中的应用[J]. 核动力工程, 2009, 30(6): 46‒48. HU Yu, ZHOU Tao. Application of proximity degree in the choice of correlation experiment in supercritical water reactor[J]. Nuclear Power Engineering, 2009, 30(6): 46‒48.

15 Wang H D, Zhou T, Luo F. Supercritical heat transfer correlation select study based on fuzzy nearness principle[J]. Advanced Materials Research, 2011, 347‒353: 2692‒2696.

Application of fuzzy nearness method in the selection of heat transfer relationship of supercritical water reactor

ZHOU Tao1,2,3LI Bing1,2,3QI Shi1,2,3MA Dongliang1,2,3HUANG Yanping4

1(School of Nuclear Science and Engineering, North China Electric Power University, Beijing 102206, China)2(Institute of Nuclear Thermal-hydraulic Safety and Standardization, North China Electric Power University, Beijing 102206, China) 3(Beijing Key Laboratory of Passive Safety Technology for Nuclear Energy, North China Electric Power University, Beijing 102206, China) 4(Key Laboratory on Nuclear Reactor Thermal Hydraulics Technology, China National Nuclear Corporation, Chengdu610041, China)

Supercritical water-cooled reactor (SCWR) is the only water-cooled reactor among the fourth generation reactors.And its thermal hydraulic process is the key to the study of SCWR, which is directly related to the safety and economy of the reactor.The degree of closeness between the determined working condition and the applicable range of empirical formula is compared to get the most suitable empirical relation.Based on the concept of closeness degree in fuzzy mathematics, the coefficient of variation method is used to determine the weight of each factor.The calculation results by Hamming closeness degree, Euclid closeness degree and maximum minimum closeness degree are basically the same.For the problem of flow and heat transfer of supercritical water reactor, a kind of optimal correlation selection method is proposed based on the known correlation.

Supercritical water-cooled reactor, Relevance, Fuzzy nearness degree,Variation coefficient method

ZHOU Tao, male, born in 1965, graduated from Xi'an Jiaotong University with a doctoral degree in 2002, focusing on nuclear engineering,thermal hydraulics and safety

2017-02-24,

2017-05-25

TL339

10.11889/j.0253-3219.2017.hjs.40.120604

周涛，男，1965年出生，2002年于西安交通大学获博士学位，研究领域为核热工水力与安全

2017-02-24，

2017-05-25

Supported by the Key Laboratory of Thermal Nuclear Hydraulic Technology for Nuclear Reactors (No.2013B40), the National Marine Major Project (No.2016YFC1402501)

中核核反应堆热工水力技术重点实验室项目(No.2013B40)、国家海洋重大专项(No.2016YFC1402501)资助