大型LNG储罐保冷性能的研究

1 概述

LNG储罐的保冷性能直接关系到储罐的安全运营，静态蒸发率是衡量储罐保冷性能的重要指标

。文献[2]建立了预测LNG储罐静态蒸发率和罐内压力变化的数学模型，分析了静态蒸发率的影响因素。文献[3]给出了多种LNG储罐静态蒸发率的测量和计算方法，并且详细论述了各种测量方法的适用情况。以往文献对大型LNG储罐保冷性能的研究主要集中在热流量和静态蒸发率计算方法的分析上，对储罐热流量影响因素的研究较少。本文以某3×10

m

双金属LNG储罐为例，给出了详细的热流量和静态蒸发率计算方法，并计算分析了环境温度、环境风速、充满率对储罐热流量的影响。

2 双金属LNG储罐基本结构

双金属LNG储罐的主体结构通常由不锈钢外罐、不锈钢内罐、热角保护层和保冷层等组成。本文以济南某双金属LNG储罐为例进行计算分析，图1和图2为该LNG储罐的保冷结构，内外罐体材质均为S30408不锈钢，内外罐之间填充保冷材料。本文约定，外罐罐壁板、内罐罐壁板以及它们之间的夹层空间为罐壁，罐壁的保冷主要依靠膨胀珍珠岩和弹性玻璃纤维毡的双层组合保冷材料；在外罐罐壁板内侧底部设有一定厚度的泡沫玻璃砖作为热角保护层；内罐罐底板以下的储罐部分为罐底，罐底的保冷材料主要为泡沫玻璃砖；穹顶、吊顶板以及它们之间的玻璃棉和BOG储存空间为罐顶，吊顶板和玻璃棉合称为吊顶，罐顶的保冷主要依靠吊顶的玻璃棉保冷层。

3 储罐热流量和静态蒸发率的计算

根据LNG储罐的结构特点，分析建立储罐热流量和静态蒸发率计算的数学模型。为了简化计算，对模型做出如下5个假设：储罐各层材料之间不存在空隙，忽略邻层之间的接触热阻。罐内处于气液平衡状态，罐内温度分布均匀且相等。将储罐与外界环境的传热过程简化为一维稳态传热过程。 LNG的组成为纯甲烷。大气压按101.325 kPa考虑。

3.1 无太阳辐射的罐壁热流量计算

① 储罐罐壁热流量的计算

在无太阳辐射时，罐壁的热流量主要是空气传递给储罐的热流量。由于LNG与BOG对内罐罐壁板内表面的对流传热强度不同，需要分别考虑BOG和LNG对应罐壁的热流量。吊顶以上罐壁内外温差很小，忽略热流量。根据圆筒壁的热流量计算公式

可推导得到无太阳辐射时储罐罐壁的热流量计算公式，见式(1)、(2)。

使用直接请求策略，主要有两种状态，其一：直接表达会话核心，只用中心行为语。其二：使用称呼语引起注意，然后表达核心内容，即起始行为语+中心行为语。如（1）“帮我拿那件衣服看看。”“我要这个。”（2）“爸爸，给我打点钱。”“爷爷，您能给我讲讲你们时代的事吗？”“xxx，清明一起去旅游吧？”采访发现，留学生直接请求策略使用频率相对不高，表意较委婉。

(1)

(2)

式中

——无太阳辐射时内罐LNG对应壁面的热流量，W

桔皮书是实施专利链接制度的基础，只有列入桔皮书中的专利才能适用这一制度。列入桔皮书的专利情况能够反映原研药企业对这一制度的利用情况，桔皮书中登记的专利的量和类型变化也反映了药品专利的发展情况和趋势。

没有选沈从文的作品，主要原因恐怕是鲁迅不认同沈从文的作品。但是，两个人的矛盾影响了选本编选的可能是存在的。1935年11月，沈从文写了《读〈中国新文学大系〉》的书评，没有谈及自己的作品是否应该入选，但认为鲁迅的选本“有抑彼扬此处”，“取舍之间不尽合理”，特别点明未选入王统照等几人的作品，对沉钟社、莽原社的评价过高，“皆与印行这套书籍的本意稍稍不合”。

——LNG液位高度，m

——环境温度，K

——LNG储罐运行温度，K

——空气对储罐外表面的表面传热系数，W/(m

·K)

——外罐外直径，m

——罐壁保冷材料的总层数

——罐壁保冷材料层数的编号

——罐壁第

层保冷材料的热导率，W/(m·K)

——罐壁第

层保冷材料的外直径，m

——罐壁第

层保冷材料的内直径，m

——LNG对罐壁内表面的表面传热系数，W/(m

·K)

——内罐内直径，m

二是成立国务院直属的跨部门、跨行政区的管理协调机构，推进一体化管理。长江大保护要一体化推进已经形成共识，建议成立国务院直属的跨部门、跨行政区的管理协调机构。早在20世纪，部分发达国家开始建立强有力的流域管理机构，它们以流域为对象，建立流域管理机构，并通过法律赋予有效的行政管理权。如美国田纳西流域管理局、法国塞纳河流域管理局、保护莱茵河国际委员会等。此外，日本“琵琶湖综合发展工程”也非常成功，均可借鉴。

——无太阳辐射时内罐BOG对应壁面的热流量，W

——内罐BOG所占空间高度，m

——BOG对罐壁内表面的表面传热系数，W/(m

·K)

空气对储罐外表面的表面传热系数随环境风速的增大而增大，可按式(3)

计算。

=6.02+3.46

（3）估计该家庭使用节水龙头后，一年能节省多少水？（一年按365天计算，同一组中的数据以这组数据所在区间中点的值作代表．）

(3)

式中

——环境风速，m/s

② BOG对罐壁内表面的表面传热系数计算

内罐BOG与罐壁内表面的对流传热可以近似看作竖直热面的大空间自然对流传热，设BOG对应罐壁内表面的温度为

，由于

和

温差较小，计算时定性温度取

，特征长度取

，其实验关联式见式(4)

。

随着社会的不断发展，以及对不动产测绘要求的提升，给不动产测绘技术也提出了新的要求。GPS技术作为新兴技术测绘手段，在不动产测绘工作中应用十分广泛，它可以实现对目标区域的精准定位，一方面提高了不动产测绘结果的准确性，另一方面也有效提高了测绘的整体效率，缩短了测绘工作的周期。同时，GPS技术在进行定位测量中，还可以根据目标区域的实际情况开展测量分析，确保测得的数据结果具有更高的可靠性。另外，不同于其他测绘技术，GPS技术在实际的应用过程中不会受到气候等外部因素的影响，可以实现全天候的测量，给不动产测绘工作的开展提供了巨大的便利。

(4)

式中

——内罐BOG与罐壁内表面对流传热过程的努塞尔数

——内罐BOG与罐壁内表面对流传热过程的格拉晓夫数

——内罐BOG的普朗特数

、

、

根据式(5)、(6)、(7)

进行计算。

2.3 4组患者治疗前后才藤氏吞咽障碍评定结果比较 治疗2周、4周后，观察3组后一时间节点治愈率较前一时间节点均有显著增加(均P<0.05)，同时，观察3组治愈率较对照组、观察1组、观察2组相同时间节点有显著增加(均P<0.05)；观察1组与对照组比较，治愈率差异无统计学意义；观察2组治疗1周、治疗2周后与对照组比较，治愈率差异无统计学意义，治疗4周后治愈率较对照组明显提高(P<0.05)。见表4。

(5)

(6)

(7)

(8)

(9)

式中

——内罐BOG的热导率，W/(m·K)

在初中信息技术教学中，教师要做的不仅是教给学生办公软件的应用、网页制造等方面的计算机知识，还要将更多的时间与精力放在学生网络学习能力的培养上来。我们当前正处在一个信息爆炸的时代，而对于中学生来说，网络无疑是他们学习知识、认识世界最有效、最便利的途径。因此，教师要在日常教学之余教给学生一些正确的快速搜索信息、处理信息的方法，使学生能够利用网络这一途径解决问题。这样不仅能够逐步改善信息技术教育的现状，使其朝着正确的方向发展，还能逐步增强学生的网络学习能力，促进学生综合素质的提高。

——重力加速度，m/s

，取9.8 m/s

——内罐BOG的体胀系数，K

——BOG对应罐壁内表面的温度，K

式中

——内罐LNG对罐底上表面对流传热过程的努塞尔数

——内罐BOG的密度，kg/m

——内罐BOG的黏度，Pa·s

g

——内罐BOG的比定压热容，J/(kg·K)

——BOG在温度为273.15 K和压力为0时的密度，kg/m

5.提升绩效人员素质，推动绩效审计不断发展。油田企业内部审计部门应结合企业特点，一方面，积极吸纳各专业技术人才加入审计队伍，并注重其财经审计知识、绩效审计方法和技术的培训；另一方面，尽快建立一套切实可行的培训制度，进一步督促审计人员对新业务的学习，快速掌握绩效审计的方法和技术，着力培养复合型人才，优化审计队伍的知识结构，提高审计人员素质，加快建设适应绩效审计要求的审计队伍，使之适应绩效审计发展的需要，组建一支具有多学科、复合型的审计队伍。

——与空气接触的承台下表面面积，m

使用云天化镁立硼复合肥+四全水溶复合肥的示范田与对照田比较，前期长势明显，中期病虫害少、落果少，后期肥劲长、不脱肥、着色好。云天化示范田与其它品牌对照田相比，同等的气候条件下，示范田樱桃亩增产28.5kg，增产率达8%。这表明，使用云天化系列化肥可以显著提高该地区樱桃的产量。

外部环境通过BOG对应罐壁传入储罐内部的热流量等于BOG从罐壁接收的热流量，即：

=π

(

-

)

(10)

联立式(4)和(10)，可求得

，然后代入式(2)可得到

。

③ LNG对罐壁内表面的表面传热系数计算

LNG与罐壁内表面的对流传热也可以近似看作竖直热面的大空间自然对流传热，计算时定性温度取

，特征长度取

，其实验关联式见式(11)

。

(11)

式中

——内罐LNG与罐壁内表面对流传热过程的努塞尔数

——内罐LNG与罐壁内表面对流传热过程的格拉晓夫数

联立式(16)和(17)，可求得

，然后代入式(14)和(15)可得到

和

。

——LNG在温度为

时的普朗特数

——LNG对应罐壁内表面的温度，K

由于

与

差别较小，在计算中校正因子(

/

)

近似取1；

根据LNG热导率计算的经验公式

进行计算；

按式(12)计算

，在

计算中LNG近似按纯甲烷处理，甲烷在不同温度下的密度可通过甲烷状态图

查得；LNG在温度为

时的黏度

和比定压热容

L

根据文献[12]中LNG黏度和比定压热容的计算方法进行计算。

(12)

式中

——LNG在

～

温度范围内的平均体胀系数，K

——LNG在温度为

时的密度，kg/m

——LNG在温度为

时的密度，kg/m

——地面温度，K

——LNG的温度范围下限，K

外部环境通过LNG对应罐壁传入储罐内部的热流量等于LNG从罐壁接收的热流量，即：

=π

(

-

)

(13)

联立式(11)和(13)，可求得

，然后代入式(1)可得到

。

3.2 罐底热流量计算

储罐罐底不受太阳辐射的影响。罐底的热流量主要包括2部分：空气以对流传热的形式传递给储罐的热流量，土壤以导热的形式传递给储罐的热流量。由于承重柱与周围环境温差较小，忽略其对流换热的热流量。地面处的承重柱温度可按地面温度进行取值，在已知地面温度的情况下可只对地面以上罐底部分进行分析。由平壁的热流量计算公式

可推导得到储罐罐底热流量的计算公式。

(14)

(15)

式中

——与空气接触的承台下表面对应罐底的热流量，W

——LNG对罐底上表面的表面传热系数，W/(m

·K)

LNG储罐内压力通常处于微正压状态，罐内LNG和BOG与压力有关的参数均按压力为0进行取值。

根据文献[6]

中天然气的热导率计算方法进行计算；

按式(8)计算

；BOG在273.15 K和压力为0时的密度

在文献[7]中查得；

根据文献[8]中天然气黏度的计算方法计算；BOG在温度为

时的比定压热容

g

在文献[9]中查得。

——罐底保冷材料的总层数

作为前明尼苏达州参议员(1964－1976)、卡特时期副总统(1977－1981)，蒙代尔仅赢得明尼苏达州和哥伦比亚特区的13张选举人票，创下选举人票最大差距纪录——13：525。败选后的蒙代尔离开政坛，重操旧业，加入明尼苏达州一家律所，做起了律师的老本行。克林顿政府时期重返政坛，担任美国驻日本大使(1993-1996)。之后从日本回到明尼苏达老家，重新在律所做事，当普通公民。

——罐底保冷材料层数的编号

——罐底第

层保冷材料的厚度，m

——罐底第

层保冷材料的热导率，W/(m·K)

10月10日，五师八十九团“棉思眠”农业合作社来了一位成都老板，他与“棉思眠”农业合作社社长陈建辉签订了2000套千层无网被合同，这是今年新棉下来后合作社签订的第一份订单合同。

——承重柱对应罐底的热流量，W

——LNG的温度范围上限，K

——承重柱总截面积，m

LNG与罐底上表面的对流传热可以看作水平热面向上的大空间自然对流传热，计算时定性温度取

，特征长度取0.25

，其实验关联式见式(16)

。

(16)

上庄乡正在向全省乃至省外大力推广本乡的舍饲养羊技术，希望各地都能通过该模式解决羊规模化养殖带来的生态平衡问题，同时加速羊只的生长速度，提高养殖户的整体经济效益，为新野县乃至河南省的畜牧产业长期可持续稳定发展创造有利条件。

——内罐LNG对罐底上表面对流传热过程的格拉晓夫数

外部环境通过罐底传入储罐内部的热流量等于LNG通过对流传热从罐底接收的热流量，即：

+

=(

+

)

(

-

)

(17)

式中

——罐底上表面的平均温度，K

——LNG在温度为

时的普朗特数

3.3 无太阳辐射的罐顶热流量计算

对于双金属LNG储罐，穹顶不具备保冷性能，在计算分析中可以近似认为穹顶内外表面温度相同且内外表面积相同。

① 传热热流量的关系

根据对流传热热流量基本计算公式

可得无太阳辐射时空气向穹顶传递的热流量为：

=

(

-

)

(18)

式中

——无太阳辐射时空气向穹顶传递的热流量，W

——穹顶外表面面积，m

——无太阳辐射时穹顶温度，K

同理可得穹顶向罐顶BOG传递的热流量为：

=

(

-

)

(19)

式中

——无太阳辐射时穹顶向罐顶BOG传递的热流量，W

——罐顶BOG对穹顶的表面传热系数，W/(m

·K)

——无太阳辐射时罐顶BOG温度，K

同理可得罐顶BOG向吊顶传递的热流量为：

=

(

-

)

(20)

式中

——无太阳辐射时罐顶BOG向吊顶传递的热流量，W

——吊顶上表面面积，m

——罐顶BOG对吊顶上表面的表面传热系数，W/(m

·K)

——无太阳辐射时吊顶上表面温度，K

根据两表面封闭系统的辐射传热热流量计算公式

，可得无太阳辐射时穹顶向吊顶辐射传热的热流量为：

(21)

式中

——无太阳辐射时穹顶向吊顶辐射传热的热流量，W

——斯忒藩-玻尔兹曼常量，W/(m

·K

)，取5.67×10

W/(m

·K

)

——穹顶内表面的发射率

——穹顶内表面对吊顶上表面的角系数

——吊顶上表面的发射率

根据平壁的热流量计算公式得到无太阳辐射时吊顶上表面向储罐内部传递的热流量为：

(22)

式中

——无太阳辐射时吊顶上表面向储罐内部传递的热流量，W

——吊顶保冷材料的总层数

——吊顶保冷材料层数的编号

——吊顶第

层保冷材料的厚度，m

——吊顶第

层保冷材料的热导率，W/(m·K)

——内罐BOG对吊顶下表面的表面传热系数，W/(m

·K)

根据对流传热热流量基本计算公式，可得无太阳辐射时吊顶下表面向内罐BOG传递的热流量为：

=

(

-

)

(23)

式中

——无太阳辐射时吊顶下表面向内罐BOG传递的热流量，W

——吊顶下表面面积，m

——吊顶下表面温度，K

② 对流传热系数的计算

穹顶近似按水平面处理，罐顶BOG与穹顶的对流传热可以看作水平热面向下的自然对流传热。由于穹顶不具有保冷性能，

和

均接近

，计算时定性温度取

，特征长度取

/(π

)，其实验关联式见式(24)

。

(24)

式中

——罐顶BOG对穹顶内表面对流传热过程的努塞尔数

——罐顶BOG对穹顶内表面对流传热过程的格拉晓夫数

——罐顶BOG在温度为

时的普朗特数

罐顶BOG与吊顶上表面的对流传热可以看作为水平冷面向上的自然对流传热，吊顶上表面和罐顶BOG温度虽然与环境温度有一定温差，但是温差较小，计算时定性温度近似取

，特征长度取

/4，其实验关联式见式(25)

。

(25)

式中

——罐顶BOG对吊顶上表面对流传热过程的努赛尔数

——罐顶BOG对吊顶上表面对流传热过程的格拉晓夫数

内罐BOG与吊顶下表面的对流传热可以看作为水平热面向下的自然对流传热，计算时定性温度近似取

，特征长度取

/4，其实验关联式见式(26)

。

(26)

式中

——罐顶BOG对吊顶上表面对流传热过程的努塞尔数

——罐顶BOG对吊顶上表面对流传热过程的格拉晓夫数

③ 吊顶向储罐内部传递的热流量求解

传入与流出穹顶的热流量相等，即：

=

+

(27)

穹顶向罐顶BOG传递的热流量等于罐顶BOG向吊顶传递的热流量，即：

=

(28)

从穹顶和罐顶BOG流入到吊顶上表面的热流量等于从吊顶上表面经吊顶各层流出到内罐BOG的热流量，见式(29)。

+

=

(29)

从吊顶上表面经吊顶各层流入到吊顶下表面的热流量等于从吊顶下表面流出到内罐BOG的热流量，见式(30)。

=

(30)

联立式(24)～(30)，可求得

和

，然后代入式(22)可求得

。

3.4 有太阳辐射的罐壁和罐顶热流量计算

在有太阳辐射时，储罐罐壁外表面和穹顶(穹顶内、外表面温度相同)在接收太阳辐射热后温度升高，增加储罐罐壁和罐顶的热流量。罐壁和穹顶外表面附近空气温度相比无太阳辐射时的温升见公式(31)

。以

表示不同朝向的储罐外表面，将储罐罐壁外表面按8个朝向各占罐壁外表面积1/8进行处理。以东朝向罐壁面为例，正东左右各偏22.5°(共45°)的罐壁圆柱面为东朝向罐壁面(

=1)，其他方向以此类推，西、南、北、东南、东北、西南、西北朝向罐壁面分别为

=2～8，穹顶外表面近似按水平表面(

=9)进行处理。储罐各朝向的日平均太阳辐射照度可在文献[13]

中查得。罐外表面的太阳辐射热吸收系数

参照GB 50176—2016《民用建筑热工设计规范》附录第B.5条选取。

(31)

式中 Δ

——太阳辐射导致的

朝向储罐外表面附近空气温度相比无太阳辐射时的温升，K

——储罐不同朝向外表面的编号

——

朝向日平均太阳辐射照度，W/m

——储罐外表面的太阳辐射热吸收系数

储罐罐壁外表面附近空气的平均温升为：

(32)

式中 Δ

——太阳辐射导致的罐壁外表面附近空气温度相比无太阳辐射时的平均温升，K

有太阳辐射的罐壁热流量计算方法与无辐射时的相同，只是计算过程中用到的环境温度不同。太阳辐射下罐壁热流量计算用到的环境温度

等于

加上温升Δ

。

同理，太阳辐射下罐顶热流量计算用到的环境温度

等于

加上温升Δ

。

3.5 LNG蒸发量的计算

本文定义

为储罐单日接收太阳辐射的时间占单日总时间(24 h)的比例，对储罐有太阳辐射和无太阳辐射时的热流量进行加权平均可得储罐总热流量的计算公式，见式(33)。

=

+

+(1-

)(

+

+

)+

(

+

+

)

(33)

式中

——储罐的总热流量，W

——储罐单日接收太阳辐射的时间占单日总时间(24 h)的比例

——有太阳辐射时LNG对应罐壁的热流量，W

——有太阳辐射时BOG对应罐壁的热流量，W

——有太阳辐射时吊顶向储罐内部传递的热流量，W

储罐24 h内LNG蒸发的质量可通过式(34)、(35)计算，其中LNG的气化潜热

可通过文献[6]

查得。

=

(34)

式中

——24 h内环境流入储罐的总热量，J

——时间，s，取86 400 s

储罐24 h内LNG蒸发的质量为：

(35)

式中

——24 h内LNG蒸发的质量，kg

——LNG的气化潜热，J/kg

3.6 静态蒸发率的计算

LNG储罐的静态蒸发率是指LNG储罐在额定充满率(额定充液高度时储存LNG的体积与LNG储罐内罐总容积之比)下，静置达到热平衡后，折算成在储罐外部环境为标准状态(1.013 25×10

Pa，293.15 K)下储罐24 h内自然蒸发产生的LNG蒸发气(BOG)的质量与储罐有效容积(额定充满率时储存LNG的体积)下LNG总质量的比例，其计算公式见式(36)。

(36)

式中

——LNG储罐的静态蒸发率

——

取293.15 K时储罐24 h内LNG蒸发的质量，kg

——LNG的密度，kg/m

——LNG储罐的有效容积，m

4 实例计算

4.1 热流量和静态蒸发率计算

应用上述计算方法对济南市某3×10

m

双金属LNG储罐的热流量和静态蒸发率进行计算分析。案例储罐的基本参数和保冷材料的物性参数见表1、2，表2中的参数参照GB 50264—2013《工业设备及管道绝热工程设计规范》附录A和GB/T 13350—2017《绝热用玻璃棉及其制品》第5.6.3条。

储罐热流量和静态蒸发率计算用到的其他参数如下：

=293.15 K；

=111.15 K；

=23.6 m，

=1.9 m；济南市全年平均风速为2.6 m/s，取

=2.6 m/s；

=0.717 4 kg/m

；

=3.24×10

Pa·s；

g

=2.202 kJ/(kg·K)；

=0.01 W/(m·K)；

=140 K时

=400 kg/m

，

=100 K时

=460 kg/m

；

=116.43×10

Pa·s；

L

=3.22 kJ/(kg·K)；

=0.48 W/(m·K)；

=293.15 K；罐底承重柱数量为144根，直径1 m，高出地面1.5 m，

=113.09 m

；

=0.8，

=0.7；

=0.03 W/(m·K)；

=10.3×10

Pa·s；

vg

=2.202 kJ/(kg·K)；

=0.48；济南不同朝向表面接收的日平均太阳辐射照度见表3

；假设储罐单日接收太阳辐射的时间为12 h，即

=0.5；在温度为111.15 K和压力为0时，

=509.74 kJ/kg，

=450.74 kg/m

。

基于以上已知参数，计算可得储罐各部分热流量，见表4。利用式(33)～(36)，计算可得储罐的静态日蒸发率

=0.074 8%。可见，储罐的静态蒸发率较小，符合国际上对大型LNG储罐的静态蒸发率控制要求：罐容为(1～5)×10

m

的大型LNG储罐的静态蒸发率≤0.08%

。

4.2 储罐热流量影响因素分析

由于储罐热流量在有无太阳辐射工况下的差别较小，并且受外部环境条件和运行工况影响的变化趋势相同，下面以无太阳辐射工况为例，对储罐各部分热流量的影响因素进行计算分析。

.

.

环境温度的影响

储罐在额定充满率和环境风速为2.6 m/s情况下各部分热流量随环境温度的变化见图3。罐壁、罐顶和罐底的热流量随环境温度的变化基本呈线性变化规律，其变化率分别为182.5、94.3、51.8 W/℃。储罐各部分热流量受环境温度影响较大。

.

.

环境风速的影响

LNG储罐外部的环境风速主要是通过改变储罐外表面的表面传热系数来影响储罐各部分的热流量，储罐在不同环境风速下的外表面表面传热系数可按式(3)计算。储罐在额定充满率和环境温度取济南全年平均气温15.15 ℃情况下各部分的热流量随环境风速的变化见图4。可见，随着环境风速的增大，储罐的各部分热流量都会增加，但是增加的幅度较小，环境风速对储罐各部分热流量的影响较小。

.

.

储罐充满率的影响

LNG储罐充满率主要影响罐壁的热流量，图5为LNG储罐在环境风速为2.6 m/s和环境温度为15.15 ℃时罐壁热流量随储罐充满率的变化曲线。可见，随着储罐充满率的增加，储罐罐壁热流量呈线性变化规律，平均每增加10%的充满率，BOG对应罐壁的热流量减少3 278.7 W，LNG对应罐壁的热流量增加3 325.9 W，储罐罐壁总热流量增加47.2 W，储罐充满率对储罐罐壁热流量的影响较小。

5 结论

根据LNG储罐的结构，给出了储罐热流量和静态蒸发率的计算方法，利用该方法对济南某3×10

m

双金属LNG储罐热流量和静态蒸发率进行计算，分析储罐热流量的影响因素，结果表明：

① 案例储罐的静态蒸发率为0.074 8%，静态蒸发率较小，储罐的保冷性能良好。

② 环境温度对储罐热流量的影响较大，环境风速和充满率对储罐热流量的影响较小。

[1] 彭明,丁乙. 全容式LNG储罐绝热性能及保冷系统研究[J]. 天然气工业,2012 (3):94-97.

[2] 王武昌,李玉星,孙法峰,等. 大型LNG储罐内压力及蒸发率的影响因素分析[J]. 天然气工业,2010 (7):87-92.

[3] 李晓明,许燕,姚淑婷,等. LNG储罐日蒸发率测量方法及计算[J]. 低温与超导,2013 (9):29-32.

[4] 杨世铭,陶文铨. 传热学[M]. 4版. 北京:高等教育出版社,2015.

[5] 李兆慈,郭保玲,严俊伟. LNG储罐温度场计算及影响因素分析[J]. 油气储运,2015 (3):244-247.

[6] 顾安忠，鲁雪生，金国强，等. 液化天然气技术手册[M]. 北京:机械工业出版社，2010.

[7] 崔永章，史永征，陈彬剑. 燃气气源[M]. 北京: 机械工业出版社，2013: 285-286.

[8] 高华伟，张树龙，李秀娟，等. 天然气黏度和热导率的计算[J]. 煤气与热力，2008 (5):B10-B13.

[ 9] 黄永华，陈国邦. 低温流体热物理性质[M]. 2版. 北京: 国防工业出版社，2014：428-432.

[10] 詹淑慧. 燃气供应[M]. 2版. 北京: 中国建筑工业出版社，2011:31-32.

[11] 段常贵. 燃气输配[M]. 5版. 北京: 中国建筑工业出版社，2015:26-27.

[12] 刘璐，刘宝玉，李少华，等. LNG热物性计算[J]. 当代化工，2010(6): 696-698.

[13] 电子工业部第十设计研究院. 空气调节设计手册[M]. 2版. 北京:中国建筑工业出版社,1995.

[14] 金明皇,许克军,程松民,等. 大型LNG储罐静态日蒸发率的计算方法[J]. 油气储运,2016 (4): 386-390.