1管道总传热系数管道总传热系数是热油管道设计和运行管理中的重要参数。
在热油管道稳态运行方案的工艺计算中,温降和压降的计算至关重要,而管道总传热系数是影响温降计算的关键因素,同时它也通过温降影响压降的计算结果。
1.1 利用管道周围埋设介质热物性计算K 值管道总传热系数K 指油流与周围介质温差为1℃时,单位时间内通过管道单位传热表面所传递的热量,它表示油流至周围介质散热的强弱。
当考虑结蜡层的热阻对管道散热的影响时,根据热量平衡方程可得如下计算表达式:1112ln 111ln22i in enwiL LD D D KD D D D (1-1)式中:K ——总传热系数,W /(m 2·℃);e D ——计算直径,m ;(对于保温管路取保温层内外径的平均值,对于无保温埋地管路可取沥青层外径);n D ——管道内直径,m ; w D ——管道最外层直径,m ;1α——油流与管内壁放热系数,W/(m 2·℃); 2α——管外壁与周围介质的放热系数,W/(m 2·℃); i λ——第i 层相应的导热系数,W/(m·℃);i D ,1i D ——管道第i 层的内外直径,m ,其中1,2,3...i n ;L D ——结蜡后的管内径,m 。
为计算总传热系数K ,需分别计算内部放热系数1α、自管壁至管道最外径的导热热阻、管道外壁或最大外围至周围环境的放热系数2α。
(1)内部放热系数1α的确定放热强度决定于原油的物理性质及流动状态,可用1α与放热准数u N 、自然对流准数r G 和流体物理性质准数r P 间的数学关系式来表示[47]。
在层流状态(Re<2000),当Pr 500Gr <时:13.65ydNu (1-2)在层流状态(Re<2000),当Pr 500Gr >时:0.250.330.430.11Pr 0.15Re Pr Pr y yyy y bdNu Gr (1-3)在激烈的紊流状态(Re>104),Pr<2500时:0.250.80.441Pr 0.021Re Pr Pr y y ybd(1-4)在过渡区(2000<Re<104)(1-5)式中:u N ——放热准数,无因次;——流体物理性质准数,无因次;——自然对流准数,无因次; ——雷诺数;0(Re )f K f =——系数;d ——管道内径,m ;g ——重力加速度,g =9.81m/s 2;υ——定性温度下的流体运动粘度,m 2/s ; C ——定性温度下的流体比热容,J/(kg·K); v q ——流体体积流量,m 3/s ;ρ——定性温度下的流体密度,kg/m 3;β——定性温度下的流体体积膨胀系数,可查得,亦可按下式计算:(1-6)f λ——定性温度下的流体导热系数,原油的导热系数f λ约在0.1~0.16W/(m ·K)间,随温度变化的关系可用下式表示:(1-7)15f ρ——l5℃时的原油密度,kg/m 3;f t ——油(液)的平均温度,℃;b t ——管内壁平均温度,℃;204d ——20℃时原油的相对密度。
注:上面各式中,参数角标f 表示以管内油(液)的平均温度f t 为定性温度;角标b 表示以管壁温度为定性温度。
(2)各处管壁导热的热阻这部分热阻包括钢管、防腐层和保温层的热阻。
钢管的导热系数g λ约为45 W/(m·℃),其热阻可忽略不计;煤焦油瓷漆防腐层导热系数f λ约为 1.1 W/(m·℃) ,黄夹克保温材料的导热系数b λ约为0.04 W/(m·℃)。
对于壁厚g δ、外包f δ厚煤焦油瓷漆防腐层的非保温热油管道,钢管及防腐层对总传热系数的影响很小。
如忽略内外径的差值,则总传热系数可近似按下式计算:12111i iK(1-8)其中:对于保温管道,保温层的热阻起决定影响。
故对于壁厚g δ、外包b δ厚保温材料的保温热油管道:1ln (22)/(2)ln(/)22n bg n gi i ibD D D D (1-9)(3) 外部放热系数2α的确定在原油长输管道内,液体的流动状态绝大部分是紊流状态,出现层流状态极少。
因此,在热力计算中,确定K 值将主要使用公式(1-1)。
在公式(1-1)中关键的参数是与管道周围许多因素有关的2α,对于埋地敷设管道:当管道的埋设深度(管中心至地表面)小于2m 时,采用下面的公式计算:(1-10)(1-11)(1-12)(1-13)式中:s λ——土壤的导热系数,W/(m·℃);e D ——与土壤接触的管道外直径,m ;ta α——土壤至地表空气间的放热系数,W/(m 2·℃);0h ——管道埋深(管中心至地表面),m 。
该放热系数包括对流放热系数tac α和辐射放热系数taR α两部分。
tac α和taR α分别用下式确定:(1-14)(1-15)式中:a v ——地表面的平均风速,m/s ;ε——土壤表面折算黑度;R C ——辐射系数,可取5.7 W/(m 2·h 4);s t ——土壤表面温度,取当地一年中月平均的最低地面温度,℃; a t ——空气温度,取当地一年中月平均的最低空气温度,℃。
当管道理设深度大于2m 时,可采用下面的公式计算2α:(1-16)式中符号的意义同前。
从上述的公式中可以看出,确定出土壤导热系数是计算埋地管道2α的关键。
土壤的导热系数与组成土壤固体物质的导热系数、土壤中固体物质颗粒大小的分布、土壤含水率、土壤状态等许多因素有关。
用理论计算很难得到准确值,因此推荐采用理论计算与参考类似管道实测值相结合的方法。
(4)结蜡层厚度计算在计入原油蜡结晶析出的潜热后,长为d x 的微元管道上,热油管道的热量平衡关系式(1-1)可简化为[48]:101lnneLD KD b b D (1-19)其中: 112ln()112i i iwD d b dD (1-20)112Lb (1-21) 如取温降为1℃时,从单位质量的原油中析出并沉积到管表面的凝油质量为biTε∂∂,则在d τ时间内在轴向温降为dT 的dx 段上沉积的量为: biLdG GdTd T(1-22)因而使内径缩小了()L d D ,则:2LLL D dG D dT dx (1-23)将式(1-23)代入式(1-22)得:2biLL D G dT d d D TT dx(1-24) 0n yK T D T T dT dxG C TkT(1-25)将式(1-25)代入式(1-24)得:012lnbinL LLyT T D D b b d D d D TTC TkT(1-26) 积分后可得:220110101ln ln22222L n L n L LbiyD D D D b b b b b b D D T T TTC TkT(1-27)对于距泵站出口x 米处的管路而言,其清管后的运行时间τ可由下式计算:L L xV其中,22243ln 2nn LL nLLDD D D D D D ,则12Ln L D D 。
联解式(1-24)与式(1-27),可求出线路上热泵站出口x m 处,经清管后运行τ小时的结蜡层内径L D τ,从而求解出结蜡厚度。
式中:y,bi——分别为油温的函数,其规律可通过试验求得;L D ——结蜡后的管道内径,m ; L λ——结蜡层的导热系数,W/( m·℃);0L τ——从下一站收到清管器开始计算的时间,s ; L D τ——运行τ小时后的结蜡层内径,m ; L τδ——运行τ小时后的结蜡厚度,m ;L ——为站间距离,m ;V ——管内流速,m/s ; k ——蜡的结晶潜热,kJ/kg 。
1.2 总传热系数的反算法热油管道稳态运行时,根据各已知的运行参数,利用苏霍夫公式反算出埋地管道总传热系数,根据《油气集输设计规范》的规定,当管道长度30L ≥km 且管径D ≥300mm 时,输油管道的热力计算应考虑管道水力摩擦生热的影响,即按列宾宗公式进行热力计算:1020lnw T T bK D LT T bGc(1-28) (1-29)式中:0T —— 管外环境温度,℃(取管道中心埋深处地温);1T ——管道起点油温,℃;2T ——终点油品温度,℃;L —— 管道长度,m ; w D —— 管道外径,m ;K ——管道至周围介质的总传热系数,W/(m 2·℃); E —— 热功当量,E =102(kg·m )/kJ ; G —— 原油质量流量,kg/h ;c —— 原油比热容,kJ/(kg·℃);b —— 由于油流在管道内摩擦功转化为热量;i —— 管道水力坡降值,m/m 。
为了更好地反映热油管道在一段时间内的实际传热状况,可采用最小二乘法来推算总传热系数K 。
设在某一运行期内,某站间管段的n 组运行参数记录值为i G ,1i G ,2iG ,1~i n ,根据最小二乘法原理可构造一个关于变量K 的无约束优化问题。
因为最小二乘法能充分利用管道的实际运行参数,在一定程度上可以消除各种随机因素的影响,反算出的K 值比较真实可靠。
采用最小二乘法拟合K 值的基本原理是:求得一个合适的K 值,使得按此K 值计算的进站温度与相应的实际记录值的偏差平方和最小。
1.2.1总传热系数因素分析及结蜡对管道温降和摩阻的影响总传热系数影响因数分析:从管道总传热系数的定义来看,影响总传热系数的根本因素是管道结构、管道埋深和管道周围土壤的性质;从总传热系数的确定方法--运行参数反算法来看,管道运行的稳定性和运行参数的测定精度是影响总传热系数测试结果的主要因素。
1管道埋深的影响:根据埋地总传热系数计算公式,管道埋深越深,管道的散热热阻越大,总传热系数越小。
2土壤性质的影响:土壤物性参数中对总传热系数影响最大的是土壤导热系数,土壤导热系数系数越大,管道总传热系数就越大。
而土壤的导热系数主要受土壤含水量的影响,土壤含水量越高,土壤的导热系数越大,管线的总传热系数也就越大。
3 管道运行工况稳定性的影响:由于轴向温降公式的前提是稳定运行工况,因而管道运行工况的稳定性对总传热系数测试结果有重大的影响,运行工况不稳定可能会导致极不合理的总传热系数测试结果。
因此,在反算总传热系数时,应当选取管线稳定运行期间的运行参数。
4 站间温降的影响:站间温降越小,抵抗运行参数波动和测量误差的能力越差,总传热系数计算结果的误差就越大。
运行参数测量精度的影响:由计算公式可知,影响总传热系数计算结果的运行参数包括输量、进出站压力和进出站温度,其中影响最大、测量精度最难保证的是进出站温度。
目前大多数油田的输油管道仍然采用套管中插玻璃温度计的方法测量进出站温度,由于套管热阻、温度计本身误差和读数误差等原因,测量结果很难反映管道中的实际油温,误差常在1℃以上,当站间温降较小时,会给总传热系数测试结果带来巨大误差。