陈志华，黄金超，何永禹，王哲，刘红波

(1.天津大学 水利工程仿真与安全国家重点实验室，天津 300072； 2.天津大学 建筑工程学院，天津 300072)

深海管道稳态流动非均匀温度场研究

陈志华1,2，黄金超2，何永禹2，王哲2，刘红波1,2

(1.天津大学 水利工程仿真与安全国家重点实验室，天津 300072； 2.天津大学 建筑工程学院，天津 300072)

为了研究深海管道稳态流动非均匀温度场，利用GAMBIT软件建立了单层管道几何模型并划分网格，利用FLUENT软件对管道在稳态流动时的非均匀温度场进行了数值模拟，对网格独立性进行了分析，确定了合理的网格划分方式，得到了管道沿程温度曲线。通过对比数值模拟与理论公式计算得到的原油温度结果，对FLUENT模型进行了验证，最后进行了管道温度场参数化分析。研究结果表明，稳态流动时管道沿程温度逐渐减小，温度梯度也逐渐减小，原油入口温度、入口流速、保温层厚度及保温层外壁与海水的换热系数均影响管道温度场的分布。研究为管道的非均匀温度效应研究提供参考。

深海管道；稳态流动；数值模拟；非均匀温度场；参数化分析

随着经济和社会的持续快速发展，中国已成为全球最大的能源消费国。我国陆地及近海油气资源开发增长已较困难，石油开采从陆地走向海洋，从浅海走向深海已经成为必然趋势。

作为海洋油气开发系统的组成部分之一，水下油气输送系统包括深海管道、立管、管汇中心、管道终端等，长期承受极端载荷作用，包括来自系统外部的环境载荷，如高水压、内波、海流及海水腐蚀的联合作用，与来自系统内部的油气输送载荷，如高油温、高油压及油气腐蚀的作用，极易发生结构失效，甚至引发泄漏事故，造成巨额的经济损失和灾难性的环境污染[1]。

为了保障流动，避免原油中的石蜡在管壁上沉积及水合物的形成，必须提高原油输送温度和输送压力，而高温和高压产生的轴向压力是导致管道整体屈曲[2]的最主要原因。

目前的研究多针对于原油温度场，Jian Su[3-6]对某多层管道在流动开启和关闭过程中的瞬时热传导进行了分析，给出了开启和关闭过程中流体沿程温度变化曲线，M. Carr[7]和H. Rong[8]在相关研究中直接给出了管道在流动开启和关闭过程中沿程温度时变曲线；盛磊祥[9]、张园园[10]和孙伟栋[11]等对管道停输后原油径向温降进行了研究，忽略了停输过程中管道轴向温度梯度，将三维不稳定传热问题简化为横截面上的二维不稳定传热问题。此外，国内的研究主要针对埋地管道[10-11]，而对于深海管道来说，埋置或者挖槽并不经济，也不现实，故多是直接铺设于海床之上，其边界条件不同于埋地管道。

综上，目前缺乏对管道温度场的研究，因此对管道温度效应的研究也是考虑恒温作用，而实际上在原油输送过程中，原油不断向周围介质散热，原油沿程温度会有所降低，由于原油与管道之间的对流换热，管道的沿程温度也随之降低。因此，有必要对裸铺管道在稳态流动时的非均匀温度场进行研究。

本文利用GAMBIT建立了单层管道几何模型并划分网格，利用FLUENT软件对管道在稳态流动时的非均匀温度场进行了数值模拟，对网格独立性进行了分析，得到了管道沿程温度曲线，并对FLUENT模型进行了验证，最后进行了管道温度场参数化分析。

1 单层管道模型建立

参考某单层管道工程，管道尺寸取为219 mm×7 mm，长度取1 000 m，保温层厚度20 mm，原油入口流速0.5 m/s，入口温度80℃。深海水温保持在4℃。由于管道直接铺设于海床之上，与海床接触面积相比与海水接触面积很小，为简化计算并节约计算成本，假设与管道接触的海床温度与周围海水温度相同，因此可根据对称性，利用轴边界条件将实际三维模型简化为二维模型，仅取旋转面计算。采用GAMBIT建立几何模型并划分网格，管道中间任意一段的GAMBIT模型如图1所示，模型中包括原油、管道和保温层三部分，其物性参数选取见表1。鉴于管道和原油之间存在热交换，为了获得管道准确的温度场，需要划分边界层网格。

图1 管道GAMBIT模型Fig.1 The GAMBIT model of the pipeline

材料密度/(kg·m-3)比热容/(J·(K·kg)-1)导热系数/(W·(m·K)-1)动力粘度系数/(kg·(m·s)-1)原油钢管保温层85578506021005007000．14480．040．04——

本模型共选用四种不同的边界条件，其中水平对称轴选用轴边界条件，原油进口选用速度进口边界条件，原油出口选用出流边界条件，其余边界条件均为壁面边界条件。由于钢管内壁(钢管与原油的接触面)和外壁(钢管与保温层的接触面)属于双侧壁面(壁面两侧均为计算域)，可将其设为耦合壁面，而无需定义其他热力学参数，求解器可以直接利用邻近网格的流场变量计算壁面上的热交换，本文即是将钢管内壁和外壁均设置成了耦合壁面。

考虑到深海管道直接裸铺在海床上，深海水温保持在4℃，而管道最外层与海水之间的换热系数较大，使得两者的温度基本相等，因此，将管道保温层的外壁定义为温度边界条件。

2 单层管道网格独立性分析

2.1 网格划分

原油和钢管之间发生对流换热，在钢管和保温层中仅存在热传导，经计算发现，钢管和保温层的网格划分方式对管道温度场无影响，故本节仅讨论原油的不同网格划分方式对管道温度场的影响。

网格质量对计算精度和稳定性有很大的影响，好的网格质量能够提高收敛性和计算精度，减少计算时间。考虑到几何模型简单且规则，而且将网格布置得与流动方向正交，可以最小化数值扩散的影响，故使用结构化网格—四边形网格。

一般情况下，在流动核心区(远离壁面的区域)网格纵横比不超过5∶1，对于边界层内的四边形网格，纵横比不超过10∶1。

综合考虑计算精度与计算成本，钢管和保温层的径向网格数分别取1和3，轴向网格尺寸分别取25 mm和40 mm。对于边界层网格，需要定义首层网格高度、增长系数以及边界层层数，不同网格划分方式如表2所示。

2.2 求解设置

故计算模型选用粘性模型中的层流模型。由于需要求解温度，故激活能量方程。

考虑到流动的不可压缩性，故采用基于压力的求解器，压力—速度耦合算法选用SIMPLE，压力插值格式选用Standard，进行稳态求解。

2.3 计算结果

采用FLUENT对管道的非均匀温度场进行了数值模拟，提取管道各网格节点温度数据，绘制管道沿程温度曲线，六种不同网格计算模型的管道沿程温度计算结果如图2所示。

表2 不同网格划分方式

图2 不同网格划分方式的计算结果Fig.2 The calculation results of different mesh generation methods

从图2中可以看出，不同网格划分方式，管道沿程温度的变化趋势基本一致，管道各处的温度梯度不同，管道入口附近，温度梯度较大，温降较快，曲线较陡，随着流动进行，温度梯度减小，温降变慢，温度曲线变缓。

分别计算了各网格划分方式计算结果相比网格划分方式6的误差，最大相对误差如表3所示。

表3 各网格划分方式计算结果的最大相对误差

从表3可以看出，最大相对误差仅1.13%，说明各网格划分方式对计算结果的影响较小，综合考虑计算成本与精度，故选择网格划分方式6。

3 模型验证

本节通过对比理论公式计算与FLUENT数值模拟得到的原油沿程温度，来验证FLUENT模型的正确性，并对比了原油沿程温度与管道沿程温度，其中FLUENT数值模拟的原油沿程温度取截面平均温度。

3.1 原油与管道之间的换热系数

不考虑原油粘度随温度的变化，管内层流换热Sieder-Tate准则关系式为[12]

(1)

经计算，可知：

故满足该准则关系式的使用范围。

将计算结果代入，有

3.2 原油与海水之间的总传热系数

原油与海水之间的总传热系数公式为[3,12]

(2)

式中:h为原油与管道之间的换热系数，W/(m2·K)；λ1为钢管的导热系数，W/(m·K)；λ2为保温层的导热系数，W/(m·K)；ri(i=1,2,3)为管道半径，m，如图3所示。将各参数取值代入，得U=1.81W/(m2·K)。

图3 管道横截面Fig.3 The cross section of the pipeline

3.3 原油沿程温降理论计算

管道入口油温为TR，℃；距离管道入口L的油温为TL，℃；管道周围介质温度为T0，℃；原油与管道周围介质之间的总传热系数为U，W/(m2·K)；管道直径为D，m；原油质量流量为G，kg/s；原油比热容为c，J/(K·kg)；L为距管道入口长度，m。原油沿程温降计算如下所述。

设在距离管道入口l(l

UπDdl(T-T0)=-GcdT

(3)

将式(3)采用分离变量法积分，可得管道内原油沿程温降关系式为

(4)

式(4)即为管道内原油沿程温降计算公式，又称苏霍夫公式。该公式没有考虑管内油流摩擦生热对温降的影响。

3.4 FLUENT数值模拟与苏霍夫公式计算结果对比

通过FLUENT数值模拟得到的原油沿程温度与通过苏霍夫公式计算得到的原油沿程温度结果如表4所示，最大相对误差0.07%，满足精度要求，从而验证了FLUENT模型的正确性。

3.5 原油沿程温度与管道沿程温度对比

通过FLUENT数值模拟得到的原油沿程温度与管道沿程温度如图4所示，从图中可以看出，管道沿程温度低于原油沿程温度，二者差值随着沿管道长度增大而增大，最大相差27%。

提取距入口250、500、750和1 000 m的原油径向温度数值并绘制成曲线，如图5所示。

表4 FLUENT模拟油温与苏霍夫公式计算油温对比

图4 原油沿程温度与管道沿程温度Fig.4 The temperature of thecrude oil and pipe along the pipeline

图5 原油径向温度Fig.5 The radial temperature of the crude oil

从图中可见，距离管道壁面一定厚度范围内，原油温度显著降低，该温度发生显著变化的薄层称为温度边界层，随着沿管道长度的增大，温度边界层的厚度也逐渐增大。温度边界层的存在导致了管道沿程温度与原油沿程温度之间的较大差异。

4 稳态流动管道温度场参数化分析

为了进一步研究稳态流动管道温度场，本节选取了原油入口温度、原油入口流速和保温层厚度三个参数，对其进行参数化分析。

4.1 原油入口温度

原油入口温度分别取80℃、70℃和60℃，原油入口流速为0.5 m/s，保温层厚度为20 mm。不同油温下的管道沿程温度曲线如图6所示，管道温降如表5所示。

图6 不同油温下的管道沿程温度曲线Fig.6 The curve of the pipe temperature along the pipeline with different crude oil inlet temperature

表5 不同原油入口温度管道温降

从图6和表5可以看出，对于不同原油入口温度，管道温度梯度均逐渐减小。随着原油入口温度的增大，管道沿程温度升高，而管道温降也有所增加，这是因为在其他条件相同的情况下，原油入口温度越高，原油散热越多，原油沿程温降越大，因此管道温降也越大。

4.2 原油入口流速

原油入口流速分别取0.1、0.3和0.5 m/s，原油入口温度为80℃，保温层厚度为20 mm。不同原油入口流速下的管道沿程温度曲线如图7所示，管道温降如表6所示。

表6 不同原油入口流速管道温降

从图7和表6可以看出，对于不同原油入口流速，管道温度梯度均逐渐减小。随着原油入口流速的增大，管道沿程温度升高，管道温降也有所减小，这是因为原油入口流速的增大导致原油与管道的换热系数减小，原油对管道的散热减小，因此管道沿程温度升高。

图7 不同原油入口流速管道沿程温度曲线Fig.7 The curve ofthe pipe temperaturealong the pipeline with different crude oil inlet velocity

4.3 保温层厚度

保温层厚度分别取20、25和30 mm，原油入口温度为80℃，入口流速为0.5 m/s。不同保温层厚度下的管道沿程温度曲线如图8所示，管道温降如表7所示。

图8 不同保温层厚度管道沿程温度曲线Fig.8 The curves of the temperaturealong the pipeline with different thermal insulation layerthickness

表7 不同保温层厚度管道温降

从图8和表7可以看出，对于不同保温层厚度，管道温度梯度均逐渐减小。随着保温层厚度的增大，管道沿程温度升高，管道温降减小，这是因为保温层厚度增大使得保温层的热阻增大，总传热系数减小，原油散失的热量减小，因此管道沿程温度升高。

4.4 保温层外壁与海水换热系数

上述分析均是将管道保温层的外壁定义为恒温的温度边界条件，实际上海水与管道之间存在对流换热，本节探索表面传热系数(换热系数)对管道温度场的影响。

换热系数分别取500、100、50和10 W/(m2·K)，海水温度取4℃，原油入口温度为80℃，保温层厚度为20 mm，入口流速取0.5 m/s。不同换热系数(恒温边界条件相当于换热系数无穷大)下的管道沿程温度曲线如图9所示，管道温降如表8所示。

图9 不同换热系数层流管道沿程温度曲线Fig.9 The curve of the temperaturealong the pipeline with different heat transfer coefficient

表8 不同换热系数管道温降

注：恒温相当于换热系数无穷大

从图9可以看出，对于不同换热系数，层流时管道温度梯度均逐渐减小。随着换热系数减小，管道温度升高，管道温度梯度减小，管道温降减小，这是因为换热系数减小使得海水与管道保温层之间的热交换减弱，海水与原油之间的总传热系数减小，原油散失的热量减小，因此管道沿程温度升高。

经计算可知，换热系数为500 W/(m2·K)时的管道沿程温度与恒温时管道沿程温度的最大相对误差仅为0.18%，换热系数为100 W/(m2·K)时的管道沿程温度的最大相对误差也仅为0.75%，而水自然对流的换热系数一般大于100 W/(m2·K)，因此将管道保温层的外壁定义为恒温的温度边界条件满足计算精度要求。

5 结论

本文利用FLUENT软件对管道在稳态流动时的非均匀温度场进行了数值模拟，得到了管道沿程温度曲线，并对管道温度场进行了参数化分析，可得出以下结论：

1)稳态流动时管道沿程温度逐渐减小，温度梯度也逐渐减小；原油径向温度分布不均匀，管道沿程温度低于原油沿程温度。

2)管道沿程温度随着原油入口温度的升高而升高，随着原油入口流速的增大而升高，随着保温层厚度的增大而升高。

3)管道温降随着原油入口温度的升高而增大，随着原油入口流速的增大而减小，随着保温层厚度的增大而减小。

4)保温层外壁与海水之间换热系数变化时，管道沿程温度与恒温边界条件管道沿程温度十分接近，满足计算精度要求。因此对于管道非均匀温度场研究，可以将管道保温层的外壁简化为恒温的温度边界条件。

5)对管道温度效应的研究应考虑管道实际温度场，本文所得到的管道沿程温度曲线可为管道非均匀温度效应研究提供参考。

[1]WANG Zhe, YU Jianxing, LIU Hongbo, et al. Study of failure mechanisms of deep-sea oil and gas transport systems under extreme conditions[J]. Science, 2014, 346: 58-62.

[2]王哲, 马克俭, 陈志华, 等. 深海管道整体屈曲研究综述[J]. 天津大学学报: 自然科学与工程技术版, 2014, 47(S): 17-23. WANG Zhe, MA Kejian, CHEN Zhihua, et al. Overview of global buckling of deep-sea pipeline[J]. Journal of Tianjin university: science and technology, 2014, 47(S): 17-23.

[3]SU Jian, CERQUEIRA D R. Simulation of transient heat transfer in multilayered composite pipeline[C]//Proceeding of OMAE01, 20th International Conference on Offshore Mechanics and Arctic Engineering. Rio de Janeiro, Brazil: American Society of Mechanical Engineers, 2001: 299-304.

[4]SU Jian, CERQUEIRA D R, ESTEFEN S F. Simulation of transient heat transfer of sandwich pipes with active electrical heating[J]. Journal of offshore mechanics and arctic engineering, 2005, 127(4): 366-370.

[5]SU Jian, REGIS C R, BOTTO A. Thermal analysis of combined active heating and passive insulation of deepwater pipelines[C]//ASME 2002 21st International Conference on Offshore Mechanics and Arctic Engineering. Oslo, Norway: American Society of Mechanical Engineers, 2002: 185-190.

[6]SU Jian, CERQUEIRA D R, ESTEFEN S F. Thermal analysis of sandwich pipes with active electrical heating[C]//ASME 2003 22nd International Conference on Offshore Mechanics and Arctic Engineering. Cancun, Mexico: American Society of Mechanical Engineers, 2003: 809-814.

[7]CARR M, SINCLAIR F, BRUTON D. Pipeline Walking-understanding the field layout challenges and analytical solutions developed for the SAFEBUCK JIP[C]//Proceedings of the 2006 Offshore Technology Conference. Houston, Texas, USA: Offshore Technology Conference, 2006.

[8]RONG Haicheng, INGLIS R, BELL G, et al. Evaluation and Mitigation of axial walking with a focus on deep water flowlines[C]// Proceedings of the 2009 Offshore Technology Conference. Houston, Texas: Offshore Technology Conference, 2009.

[9]盛磊祥, 许亮斌, 蒋世全, 等. 基于Fluent的海管内停输管内原油温度变化过程分析[J]. 中国造船, 2009, 50(S): 333-337. SHENG Leixiang, XU Liangbin, JIANG Shiquan, et al. Analysis of oil temperature decreasing trend in shut down submarine delivery pipe based on fluent[J]. Shipbuilding of China, 2009, 50(S): 333-337.

[10]张园园. 热油管道停输降温过程的数值计算[D]. 青岛: 中国石油大学(华东), 2007. ZHANG Yuanyuan. Numerical calculation of hot oil pipeline temperature drop process after shutdown[D]. Qingdao: China University of Petroleum (East China), 2007.

[11]孙伟栋. 海底输油管道传热模拟计算[D]. 大庆: 大庆石油学院, 2007. SUN Weidong. The heat transfer calculation and analysis of submarine oil pipeline[D]. Daqing: Daqing Petroleum Institute, 2007.

[12]许国良, 王晓墨, 邬田华, 等. 工程传热学[M]. 北京: 中国电力出版社, 2005.

Non-uniform temperature fields of a deep-sea pipeline in steady flow

CHEN Zhihua1,2, HUANG Jinchao2, HE Yongyu2, WANG Zhe2, LIU Hongbo1,2

(1.State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300072，China; 2.School of Civil Engineering, Tianjin University, Tianjin 300072，China)

To study the non-uniform temperature field of a deep-sea pipeline in steady flow, this paper established a geometric model of a single-layer pipeline and meshed it with GAMBIT. Numerical simulation of the pipeline in steady flow was then conducted under the condition of a non-uniform temperature using FLUENT software. The mesh independence of the system was analyzed, a reasonable method of mesh division was adopted, and the curve of the temperature along the pipeline was obtained. The FLUENT model was verified by comparing the oil temperature results of the numerical simulation with those of the theoretical formula. Finally, parametric analysis of the temperature field of the pipeline was carried out. Study results showed that both the temperature along the pipeline in steady flow and the temperature gradient decreased gradually. The oil inlet temperature, oil inlet velocity, thermal insulation layer thickness, and heat transfer coefficient between the thermal insulator and the seawater influenced the distribution of the pipeline temperature field. This paper provides a reliable reference for the study of non-uniform temperature effects on the pipeline.

deep-sea pipeline; steady flow; numerical simulation; non-uniform temperature field; parametric analysis

2016-01-18.

日期：2016-11-16.

国家重点基础研究发展计划(2014CB046801).

陈志华(1966-), 男,教授,博士生导师； 刘红波(1983-), 男,副教授.

刘红波，E-mail:hb_liu2008@163.com.

10.11990/jheu.201601065

http://www.cnki.net/kcms/detail/23.1390.u.20161116.1613.018.html

TE973.1

A

1006-7043(2017)02-0189-07

陈志华，黄金超，何永禹，等. 深海管道稳态流动非均匀温度场研究[J]. 哈尔滨工程大学学报， 2017, 38(2): 189-194. CHEN Zhihua, HUANG Jinchao, HE Yongyu， et al. Non-uniform temperature fields of a deep-sea pipeline in steady flow[J]. Journal of Harbin Engineering University, 2017, 38(2): 189-194.