CO2水力模型构建及管径确定方法.docxVIP

超临界二氧化碳管道输运工艺计算方法 （1）开展管径Di计算的参数要开展管道内径D计算，需要给出的参数包括：令二氧化碳质量流量Qm，以每年的最大质量流量确定；令管道长度L；令管道入口压力P、出口压力P2令管道入口温度 （1）开展管径Di计算的参数 要开展管道内径D计算，需要给出的参数包括： 令二氧化碳质量流量Qm，以每年的最大质量流量确定； 令管道长度L； 令管道入口压力P、出口压力P2 令管道入口温度T1、出口温度72。 个 管道入口海拔高度”1、出口海拔高度”2。 （2）平均压力Pave、平均温度Tave计算方法 在对能量方程式中各项积分时，采用了平均压力Pave、平均温度Tave、平均 压缩因子Zave的概念。因此，必须确定Pave、Tave、Zave的数值。 由于压力、温度沿管道长度L的分布并非线性，因此，可采用下式计算管 道沿线平均压力Pave和平均温度Tave： %3=|(02+01_筌 3v2十产1 1 ave -十 .十, (19) (20) （3）摩擦系数2计算 采用经验公式计算摩擦系数2： 当 Re2000 时： c 64 A =— Re (33) 当Re在2000-10000之间时: (Z + B)3/2 (34) A = [2.4571n(-o；+注]- I (35) 8 =( 37530 16Re) (36) 式（35）中，2为管道粗糙度，对钢管，一般取0.046。Re为雷诺数，可由下式计算得到: Re =生 Re =生2 式中，〃为超临界二氧化碳粘度。 (37) 2.5管径计算流程框图 在采用公式(15)计算管道内径D时，由于摩擦系数2与雷诺数上有关， 而Re与管径D有关，故需要迭代计算。首先预估一个初始管径。o,继而计算 雷诺数Re、摩擦系数九 结合其他参数计算管径内径O],假设。与Qo的差 值在10-6之内，那么满足要求，所求的管道内径假设。与为偏差超过10- 6,那么以计算得到的Di作为初始值重新计算，直到Oj+i与Dj的偏差在10-6之 内，所求直径D=Dj+i。由于管道直径均为标准系列，因此求得D后，还需要 依据管道系列确定实际管道规格，并计算最终使用的管道内径Dr。 计算流程框图见图20 管径A=Oj+i —选择大于Di的标准化系列管径 图2管径计算流程图 2超临界C02状态方程及热、动力学特性参数计算方法 明确超临界C02状态方程及其热、动力学物性参数如压缩因子、密度、粘 度、比容热、导热系数等的计算方法，是后续构建输运过程水力模型、传热模型 等的基础。本章解决这些问题。 1超临界C02状态方程 已有大量研究说明，Peng-Robinson方程(简称PR方程，下同)最适合作 为CO2的状态方程。PR方程关系式如下： RmT aVm-b RmT a Vm-b *+2 匕%匕 2 _ 0.45724R备索 Pc (2.1) (2.2) (2.2) (2.5)(2.6)f = [1 + (0.37464 + 1.54226a - 0.26992户)(1 一 冏)[2(2.4) (2.5) (2.6) 70.0778RmTc b = Pc 式(2.1)?(2.6)中： P：介质绝对压力，Pa； T：介质温度，K； Rm：通用气体常数，IX 8.314 J/(mol-K)； Vm：介质摩尔容积，m3/mol； Tr：比照温度； Tc：介质临界温度,K,对CO2, ￡=304.2K。 Pc：介质临界压力，Pa,对 CO2, pc=7.375xl06Pao ①：偏心因子。对CO2, co=0.225o a、小、人均为计算过程中的中间参数。4为比照温度，和偏心因子切的函 数。通常在文献中，该参数用a表示。本文为了将。与。区分，故标记为短 通常在使用状态方程时，介质压力p和温度r为参数，需要求其他参 数如摩尔容积Vm或者压缩因子Z。由于式（2.1）为非线性方程，无法获得解析 解，故需要采用牛顿法迭代编程计算。 压力因子计算方法 利用气体状态方程（2.1）迭代计算得到CO2的摩尔容积心后，即可依据 压缩因子定义： (2.7) 得到计算压缩因子Z的公式: (2.8)密度计算方法 CO2的密度/与压力p、温度T有关。虽然通过气体状态方程（2.1）计算 摩尔热容片可以获得密度数值，但在本文选择一种更为便捷的计算方法。文献 通过大量数据验证，提出了压力p在7.6MPa?24.8MPa之间，温度7在-1.1? 82.2 ℃之间的CO2密度拟合公式。该拟合公式精度高且便于编程计算。 该文献指出，不同等级温度下的密度p与压力（公式中以％表示）关系的 拟合公式可采用式（2.9）的多项式形式，各温度等级对应的拟合公式系数见表 2.U当实际介质温度与表中各温度等级不同时，可选择与实际温度上、下临近 的温度等级对应的拟合公式