3. LNG物性数据获取与处理

板翅式换热器的温差与压降校核，本质上是在求解一组耦合的传热与流动方程，而方程中几乎每一个系数都依赖流体物性。对LNG而言，这一点尤为敏感：它是以甲烷为主、含少量乙烷、丙烷、氮等的多组分低温混合物，温度通常跨越-160℃至-100℃的深冷区间，压力从常压储罐到8MPa以上的高压泵出口都有可能。在这样的范围内，密度可变化数倍，比热在临界区附近出现陡峰，粘度与导热系数也随温度和相态显著漂移。若物性取值偏差10%，换热面积或压降的校核结论就可能出现方向性错误。因此，本节要解决的核心问题是：物性从哪里来、在什么条件下取、如何组织成可被校核公式直接调用的形式。

一、先明确物性输入的三要素

任何一组LNG物性都不是孤立数字，而是由三个条件共同锁定的：组成、温度、压力。组成决定了临界参数和相平衡行为，温度压力决定了所处的相态与对比状态。工程上常见的错误是只给温度不给压力，或用一个平均组成代替实际组分。正确的做法是先建立组分表，再明确每个校核点的T、P，然后逐点求物性。

以一个典型示例组成为例（以下数据为教学假设，非实测）：甲烷92%、乙烷5%、丙烷2%、氮1%（摩尔分数）。该混合物在3MPa、-150℃下处于压缩液相，在0.3MPa、-160℃下则接近饱和液。两种状态下的密度可能相差20%以上，绝不能混用。

二、物性获取的三条路径

路径一：状态方程法。对LNG这类轻烃混合物，PR方程或SRK方程配合经典混合规则是工程主流。给定组成、T、P，可解出压缩因子Z、密度、焓与熵，再通过剩余性质或焓差法得到比热。优点是组成变化可追溯，缺点是近临界区精度下降，需要引入体积平移或二元交互系数修正。

路径二：物性数据库与图表。NIST REFPROP、GERG-2008等模型对天然气体系精度较高，适合作为基准。工程现场也常用物性手册中的LNG密度-温度曲线、比热曲线。使用图表时必须注意其基准组成和压力条件，外推到其他工况要谨慎。

路径三：实测与回归。若有在线色谱和密度计数据，可用实测值回归局部物性关联式，用于在线校核。这种方法最贴近实际，但需要保证取样代表性，避免分层或闪蒸导致组成失真。

三、把物性组织成可调用的曲线

换热校核不是单点计算，而是沿换热器长度逐段积分。因此物性应以“温度-物性”曲线或分段表格的形式提供。典型做法是：在LNG侧温度区间内取若干节点（如每5℃一个点），在每个节点上按当地压力求密度ρ、定压比热cp、导热系数λ、动力粘度μ，并计算普朗特数Pr=μcp/λ。板翅式换热器的翅片通道水力直径小，雷诺数常在数百到数千之间，Pr的变化直接影响努塞尔数关联式的选取。

下图展示了从组成与工况输入，到物性曲线生成，再到换热校核调用的完整数据链。图中同时标出了四个关键物性在深冷区间的典型变化趋势，帮助理解为何必须逐段取值而非取平均值。

四、物性曲线在换热校核中的具体用法

拿到物性曲线后，校核流程通常按以下方式展开。以热端到冷端分段为例，每段取进出口平均温度对应的物性，计算该段的雷诺数Re=ρ·u·d/μ、普朗特数Pr，再代入翅片通道的传热关联式（如锯齿翅片的j因子关联式）求努塞尔数和对流换热系数h。同时用摩擦因子关联式求压降梯度dP/dx。由于LNG侧温度变化大，段数不宜过少，工程上常取8~15段，使每段内物性变化控制在5%以内。

一个常见的判断依据是：若相邻两段算出的h相差超过15%，说明分段过粗或物性曲线节点不足，应加密节点重新计算。压降方面，若某段液相密度因接近饱和而骤降，该段压降会明显放大，这正是需要重点关注的闪蒸风险区。

五、一份可直接使用的物性数据表结构

为便于校核调用，建议将物性整理为如下字段的表格（示例数值为教学假设）：

表中Pr由μ·cp/λ计算得到，注意单位统一：μ用Pa·s、cp用J/(kg·K)、λ用W/(m·K)。这张表可直接作为校核程序的输入，也可用于手工估算。实际工程中还应补充饱和温度、临界压力等辅助字段，以便判断相态边界。

最后需要强调，物性数据的可靠性直接决定校核结论的可信度。在提交校核报告时，应同时记录物性来源、模型版本、组成基准和取值节点，使结果可复现、可追溯。这不仅是技术习惯，也是低温换热器设计评审中的基本要求。