热油管道温降计算
- 格式:ppt
- 大小:517.50 KB
- 文档页数:51
王霁:油田站间输油管道热力优化计算第11卷第6期(2021-06)温度是热油管道输送过程中最重要的参数之一,温度高低直接影响油品的黏度、比热、密度等物性。
出站温度过高会导致输送过程中的散热量增加,能量损失变大;而出站温度过低,会导致安全停输时间减少,增加安全隐患。
因此选择合理的出站温度对于热油管道安全经济输送具有重要意义。
在热油管道生产运行参数优化过程中,管道温降计算的精度对于方案的合理性影响巨大。
1现状为提高输油管道温降计算的精度,许多学者进行了大量研究,徐睿妤为解决现有热力模型上存在的不足,根据热阻原理编制了具备常见模型模块和自定义模块的OHTC计算软件,提高了管道运输上热力计算的精准度[1];庞海涛采用PIPESIM流体计算软件进行对K值影响因素的模拟,分析得出了各因素对其影响程度[2];魏立新,周刚等提出一种更适用于现场实际工况的基于相关向量机算法(RVM)的埋地热油管道温降预测的新方法[3];杨新明、蒋洪等在利用探针法测量土壤导热系数的基础上,结合管道设计数据,通过理论公式计算求解管道的总传热系数[4-5];杨加栋,张晓灵等以聚氨酯弹性体保温管道作为研究对象,设计可模拟深水环境下压力和温度的管道散热测试装置,有助于分析测验材料保温性能,为工程现场提供了有力参考[6]。
针对某油田站间原油输送管道,采用最优化拟合方法建立了管道热力计算模型,并对不同输量不同月份管道的合理出站温度进行了计算。
2热力计算模型某油田站间管道全长30km,采用加热输送方式。
该管道环境温度-8~15℃,油品比热容1.94~2.69J/(kg·℃),黏度23.5~43.8mPa·s,密度0.87kg/m3,进站温度39~49℃。
热油管道温降通常采用苏霍夫公式:TZ=T0+()T R-T0e-aLa=KπDGc(1)式中:G为油品的质量流量,kg/s;c为输油温度下油品的比热容,J/(kg·℃);K为管道总传热系数,W/(m2∙℃);D为管道外径,m;L为管道长度,m;TR为管道起点温度,℃;TZ为管道终点温度,℃;T0为管道周围土壤温度,℃;a为参数。
油气集输管线温降计算方法摘要:在油管散热以及沿程压降等因素的影响下,油气集输管线会出现温降等问题,科学计算温降有利于促进石油工业的发展,因此本文利用文献资料法等方法对油气集输管线温降计算方法进行了研究与探讨。
在探究过程中先分析了相关的数学模型,例如热力计算模型、传热系数计算模型等,之后分析了计算方法以及具体案例。
在分析后发现通过能量平衡微分方程等模型以及计算程序可以增强计算结果的准确性,减少计算误差,所以需要提高对这一计算方法的重视程度。
关键词:油气集输管线;温降;计算方法前言:温降是油气集输管线经常出现的问题,在进行油气开采管理时做好温降计算工作有利于增强管线运行的安全性与经济性。
但利用传统的苏霍夫温降公式进行估算会降低计算结果的准确性,所以需要探讨新的、精准的计算方法,从而增强温降计算结果的精准性,为管理工作提供依据。
1.油气集输管线温降相关数学模型1.1热力计算模型由于液相中含有不相溶的油和水这两种液体,所以石油工业中的油、气、水混合物流动属于多液多相流动。
但从实际情况来看,其流动的力学关系类似于气液两相流动相,所以将其归纳到气液两相流动的范畴中。
在这种情况下,若假设两相之间没有温度滑梯,且不计油品径向温度阶梯,就可以将气液两相混合物沿管线的能量微分方程设置为公式(1)【1】。
(1)公式(1)当中的q指的是垂直于管壁方向的热流量、H指的是混合物焓、V指的是混合物平均速度、g为重力加速度、θ指的是管轴线与水平面夹角。
但混合物焓会受到其压力和温度的影响,所以可以利用公式(2)表示。
(2)公式(2)中的为焦耳-汤普森系数、为混合物定压比热、为混合物焓自身的温度。
从公式(1)与公式(2)可以获取公式(3)。
(3)公式(3)中的P为混合物焓自身的压力,且=-2π(4)公式(4)中的负号表示散热,k为传热系数,为环境土壤温度。
从公式(3)与公式(4)可以得到公式(5)。
(5)公式(5)中的, - - 。
1管道总传热系数管道总传热系数是热油管道设计和运行管理中的重要参数。
在热油管道稳态运行方案的工艺计算中,温降和压降的计算至关重要,而管道总传热系数是影响温降计算的关键因素,同时它也通过温降影响压降的计算结果。
1.1 利用管道周围埋设介质热物性计算K 值管道总传热系数K 指油流与周围介质温差为1℃时,单位时间内通过管道单位传热表面所传递的热量,它表示油流至周围介质散热的强弱。
当考虑结蜡层的热阻对管道散热的影响时,根据热量平衡方程可得如下计算表达式:1112ln 111ln 22i i ne n w i L L D D D KD D D D a a l l -+轾骣犏琪桫犏=+++犏犏犏臌å (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 <时:1 3.65y d Nu a l== (1-2) 在层流状态(Re<2000),当Pr 500Gr >时: 0.250.330.430.11Pr 0.15Re Pr Pr y y y y y b d Nu Gr a l 骣琪==鬃琪桫(1-3) 在激烈的紊流状态(Re>104),Pr<2500时: 0.250.80.441Pr 0.021Re Pr Pr y y y b d l a 骣琪=鬃琪桫 (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℃时原油的相对密度。
管道温降计算1管道总传热系数管道总传热系数K 指油流与周围介质温差为1℃时,单位时间内通过管道单位传热表面所传递的热量,它表示油流至周围介质散热的强弱。
当考虑结蜡层的热阻对管道散热的影响时,根据热量平衡方程可得如下计算表达式:1112ln 111ln 22i i n e n wi L L D 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 ;L λ——所结蜡导热系数。
为计算总传热系数K ,需分别计算内部放热系数1α、自管壁至管道最外径的导热热阻、管道外壁或最大外围至周围环境的放热系数2α。
(1)内部放热系数1α的确定放热强度决定于原油的物理性质及流动状态,可用1α与放热准数u N 、自然对流准数r G 和流体物理性质准数r P 间的数学关系式来表示。
在层流状态(Re<2000),当500Pr1 3.65y dNu αλ== (1-2)在层流状态(Re<2000),当500Pr >?Gr 时:0.250.330.430.11Pr 0.15Re Pr Pr y y y y y b d Nu Gr αλ??==??(1-3)在激烈的紊流状态(Re>104),Pr<2500时:0.250.80.441Pr 0.021Re Pr Pr y y yb d λα??=?? ???(1-4)在过渡区(2000<re<104)< p="">25.043.001)Pr Pr (Pr bf f fdK ?λα=(1-5)式中:u N ——放热准数,无因次;λρυC =Pr ——流体物理性质准数,无因次; ()υβw f t t g d Gr -=3——自然对流准数,无因次;υπρd q vdv4Re ==——雷诺数; )(Re 0f f K =——系数;d ——管道内径,m ;g ——重力加速度,g =9.81m/s 2;υ——定性温度下的流体运动粘度,m 2/s ;C ——定性温度下的流体比热容,J/(kg·K);v q ——流体体积流量,m 3/s ;ρ——定性温度下的流体密度,kg/m 3;β——定性温度下的流体体积膨胀系数,可查得,亦可按下式计算:td d -+-=2042045965634023101β (1-6)f λ——定性温度下的流体导热系数,原油的导热系数f λ约在0.1~0.16W/(m·K)间,随温度变化的关系可用下式表示:153/)1054.01(137.0f t f t ρλ-?-= (1-7)15f ρ——l5℃时的原油密度,kg/m 3; f t ——油(液)的平均温度,℃; b t ——管内壁平均温度,℃; 204d ——20℃时原油的相对密度。
1.1.1 压降、温降计算公式根据《动力管道手册》压降计算公式:)(10)(10215.11232H H Ld L d w p -++⨯=∆ρλρ式中:1.15——安全系数;p ∆——介质沿管道内流动总阻力,Pa ;L ——为管道直线长度m ;Ld ——为管道局部阻力当量长度m ;W ——蒸汽管道平均流速m/s ;d ——管道内径mm ;ρ——蒸汽介质平均密度kg/m 3;λ——管道摩擦阻力系数,根据管道绝对粗糙度K 值选择,对过热蒸汽管道,按管道绝对粗糙度K=0.1mm 取用;H2-H1——管道终端与始端的高差,m 。
根据《设备及管道绝热设计导则》GB/T8175-2008 单层保温的管道单位热损失计算公式:Do Di Do In Ta T R R Ta T q ∙+-=+-=αλαπ21)(21 W/m.h 式中:T ——设备和管道的外表面温度(℃),T 应取管道蒸汽介质的平均温度即221t t T +=;t1——管道始端蒸汽温度℃;t2——管道终端蒸汽温度℃;Ta ——环境温度,根据工程情况定℃;R1——保温层热阻 对管道(m.K )/W ;对平面:(m 2.K)/W ; R2——保温层表面热阻 对管道(m.K )/W ;λ——保温材料制品在平均温度下导热系数W/(m.K ); Do ——保温层外径 m ;Di ——保温层内径 m ;α——保温层外表面与大气的换热系数 W/(m 2.K ),w 36α+= GB/T8175-2008规范推荐.K W/m .α26311= 此时风速w 为3.5m/s 。
管径计算是按照正常负荷计算管径,同时以最大负荷及最小负荷校核计算后综合选取的。
1管道总传热系数管道总传热系数K 指油流与周围介质温差为1℃时,单位时间内通过管道单位传热表面所传递的热量,它表示油流至周围介质散热的强弱。
当考虑结蜡层的热阻对管道散热的影响时,根据热量平衡方程可得如下计算表达式:1112ln 111ln 22i i n e n wi L L D 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 ;L λ——所结蜡导热系数。
为计算总传热系数K ,需分别计算内部放热系数1α、自管壁至管道最外径的导热热阻、管道外壁或最大外围至周围环境的放热系数2α。
(1)内部放热系数1α的确定放热强度决定于原油的物理性质及流动状态,可用1α与放热准数u N 、自然对流准数r G 和流体物理性质准数r P 间的数学关系式来表示。
在层流状态(Re<2000),当500Pr <⋅Gr 时:1 3.65y dNu αλ== (1-2) 在层流状态(Re<2000),当500Pr >⋅Gr 时:0.250.330.430.11Pr 0.15Re Pr Pr y y y y y b d Nu Gr αλ⎛⎫==⋅⋅⎪⎝⎭(1-3)在激烈的紊流状态(Re>104),Pr<2500时:0.250.80.441Pr 0.021Re Pr Pr y y y b d λα⎛⎫=⋅⋅ ⎪⎝⎭(1-4)在过渡区(2000<Re<104)25.043.001)Pr Pr (Pr bf f fdK ⋅λα= (1-5)式中:u N ——放热准数,无因次;λρυC =Pr ——流体物理性质准数,无因次; ()υβw f t t g d Gr -=3——自然对流准数,无因次;υπρd q vdv4Re ==——雷诺数; )(Re 0f f K =——系数;d ——管道内径,m ;g ——重力加速度,g =9.81m/s 2; υ——定性温度下的流体运动粘度,m 2/s ; C ——定性温度下的流体比热容,J/(kg·K); v q ——流体体积流量,m 3/s ;ρ——定性温度下的流体密度,kg/m 3;β——定性温度下的流体体积膨胀系数,可查得,亦可按下式计算:td d -+-=2042045965634023101β (1-6)f λ——定性温度下的流体导热系数,原油的导热系数f λ约在0.1~0.16W/(m·K)间,随温度变化的关系可用下式表示:153/)1054.01(137.0f t f t ρλ-⨯-= (1-7)15f ρ——l5℃时的原油密度,kg/m 3;f t ——油(液)的平均温度,℃;b t ——管内壁平均温度,℃;204d ——20℃时原油的相对密度。
1管道总传热系数管道总传热系数K 指油流与周围介质温差为1℃时,单位时间内通过管道单位传热表面所传递的热量,它表示油流至周围介质散热的强弱。
当考虑结蜡层的热阻对管道散热的影响时,根据热量平衡方程可得如下计算表达式:1112ln 111ln 22i i n e n w i L L D 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 ;L λ——所结蜡导热系数。
为计算总传热系数K ,需分别计算内部放热系数1α、自管壁至管道最外径的导热热阻、管道外壁或最大外围至周围环境的放热系数2α。
(1)内部放热系数1α的确定放热强度决定于原油的物理性质及流动状态,可用1α与放热准数u N 、自然对流准数r G 和流体物理性质准数r P 间的数学关系式来表示。
在层流状态(Re<2000),当500Pr <⋅Gr 时:1 3.65y dNu αλ== (1-2) 在层流状态(Re<2000),当500Pr >⋅Gr 时:0.250.330.430.11Pr 0.15Re Pr Pr y y y y y b d Nu Gr αλ⎛⎫==⋅⋅⎪⎝⎭(1-3)在激烈的紊流状态(Re>104),Pr<2500时:0.250.80.441Pr 0.021Re Pr Pr y y yb d λα⎛⎫=⋅⋅ ⎪⎝⎭(1-4)在过渡区(2000<Re<104)25.043.001)Pr Pr (Pr bf f fdK ⋅λα= (1-5)式中:u N ——放热准数,无因次;λρυC =Pr ——流体物理性质准数,无因次; ()υβw f t t g d Gr -=3——自然对流准数,无因次;υπρd q vdv4Re ==——雷诺数; )(Re 0f f K =——系数;d ——管道内径,m ;——重力加速度,g =9.81m/s 2;υ——定性温度下的流体运动粘度,m 2/s ;C ——定性温度下的流体比热容,J/(kg·K);v q ——流体体积流量,m 3/s ;ρ——定性温度下的流体密度,kg/m 3;——定性温度下的流体体积膨胀系数,可查得,亦可按下式计算:tdd-+-=2042045965634023101β (1-6)f λ——定性温度下的流体导热系数,原油的导热系数f λ约在0.1~0.16W/(m·K)间,随温度变化的关系可用下式表示:153/)1054.01(137.0f t f t ρλ-⨯-= (1-7)15f ρ——l5℃时的原油密度,kg/m 3;f t ——油(液)的平均温度,℃; b t ——管内壁平均温度,℃;204d ——20℃时原油的相对密度。
热油输送管路的温降计算热油在埋地管路输送过程中因无法做到完全绝热,它会沿管线向四周传热,下面仅以纵向温降进行研究计算。
1·设热油输送管道,管外径为D ,周围介质温度为T 0,总传热系数为K ,输量为G ,油品的比热为C ,出站油温为T Q ,油流流到距加热站出口X 米处时,温度降为T ℃。
注:(1)在稳定工况下:温度不随时间而变化,输量不随时间而变化;(2)油流至周围介质的总传热系数K 沿线为常数; (3)沿线地温和油品的比热C 为常数; (4)油品沿管轴线温度不变。
2·在距输油站为X 处取一微元段dx ,设X 处断面油温为T ,油流经过dx 段的温度变化为dt ,由能量方程推导温降公式,稳定流动的能量方程; dx dQ g dx dv v dxdPP h dx dT T h T p -=++⎪⎭⎫ ⎝⎛∂∂+⎪⎭⎫⎝⎛∂∂θsin 忽略高差和速度变化的影响,则: dx dQ dxdPP h dx dT T h T p -=⎪⎭⎫ ⎝⎛∂∂+⎪⎭⎫⎝⎛∂∂另外由热力学知识可知: hp T P T T h P h ⎪⎭⎫⎝⎛∂∂⎪⎭⎫ ⎝⎛∂∂-=⎪⎭⎫ ⎝⎛∂∂因此: dx dQ dx dPP T T h dx dT T h h p p -=⎪⎭⎫ ⎝⎛∂∂⎪⎭⎫ ⎝⎛∂∂-⎪⎭⎫⎝⎛∂∂ 由于: P p C T h =⎪⎭⎫ ⎝⎛∂∂ i h D P T =⎪⎭⎫⎝⎛∂∂ 则:dQ dp D C dT C i P P -=-故在L+dL 断面上油温为T+dT ,稳定传热时,dL 上的热平衡方程为:单位时间内管线向周围介质的散热量 = 油流温降放出的热量dQ 表示单位质量液体在单位管长上的热量损失,由传热学关系可知: ()dx MT T D K dQ 0-=π因此: ()d T C d p D C d x MT T D K P i P -=-0π令 PMCD K a π=, 则:()()dxdp D T T a dxT T d i=-+-00非齐次线性微分方程的通解为: dxedxdp D e Ce T T axi ax ax ---⎰+=-0由于: 0=x 时,Q T T =,所以:()dxedxdp eD eT T T T axaxi axQ ---⎰+-+=00在热油液流中不考虑节流效应,则得到苏霍夫公式: ()ax Q e T T T T --+=00单位质量下取:PGCD K a π=适用于流速低、温降大、摩阻热影响较小的情况下。
管道总传热系数算————————————————————————————————作者:————————————————————————————————日期:1管道总传热系数管道总传热系数是热油管道设计和运行管理中的重要参数。
在热油管道稳态运行方案的工艺计算中,温降和压降的计算至关重要,而管道总传热系数是影响温降计算的关键因素,同时它也通过温降影响压降的计算结果。
1.1 利用管道周围埋设介质热物性计算K 值管道总传热系数K 指油流与周围介质温差为1℃时,单位时间内通过管道单位传热表面所传递的热量,它表示油流至周围介质散热的强弱。
当考虑结蜡层的热阻对管道散热的影响时,根据热量平衡方程可得如下计算表达式:1112ln 111ln 22i i n e n w i L L D 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),当500Pr <⋅Gr 时:1 3.65y d Nu αλ== (1-2) 在层流状态(Re<2000),当500Pr >⋅Gr 时: 0.250.330.430.11Pr 0.15Re Pr Pr y y y y y b d Nu Gr αλ⎛⎫==⋅⋅ ⎪⎝⎭ (1-3)在激烈的紊流状态(Re>104),Pr<2500时:0.250.80.441Pr 0.021Re Pr Pr y y y b d λα⎛⎫=⋅⋅ ⎪⎝⎭ (1-4)在过渡区(2000<Re<104)25.043.001)Pr Pr (Pr b ff f d K ⋅λα= (1-5)式中:u N ——放热准数,无因次;λρυC =Pr ——流体物理性质准数,无因次; ()υβw f t t g d Gr -=3——自然对流准数,无因次; υπρd q vd v 4Re ==——雷诺数; )(Re 0f f K =——系数;d ——管道内径,m ;g ——重力加速度,g =9.81m/s 2;υ——定性温度下的流体运动粘度,m 2/s ;C ——定性温度下的流体比热容,J/(kg·K); v q ——流体体积流量,m 3/s ;ρ——定性温度下的流体密度,kg/m 3;β——定性温度下的流体体积膨胀系数,可查得,亦可按下式计算: t d d -+-=2042045965634023101β (1-6)f λ——定性温度下的流体导热系数,原油的导热系数f λ约在0.1~0.16 W/(m ·K)间,随温度变化的关系可用下式表示:153/)1054.01(137.0f t f t ρλ-⨯-= (1-7)15f ρ——l5℃时的原油密度,kg/m 3;f t ——油(液)的平均温度,℃;b t ——管内壁平均温度,℃;204d ——20℃时原油的相对密度。
第三章热油输送管道的工艺计算(Hot-oil Pipelines)随着世界能源需求的增长,易凝和高粘原油的产量不断地增加。
目前我国所产原油大多为这两种原油。
生产含蜡原油(waxy crude)(即易凝原油)的油田主要有:大庆油田、胜利油田、中原油田、华北油田、河南油田、长庆油田、克拉玛依油田。
生产稠油(thick oil ,heavey oil)的油田有:辽河油田、胜利的单家寺油田和孤岛油田等。
含蜡原油的特点是含蜡量高、凝固点高、低温下粘度高、高温下粘度低。
如大庆原油,凝固点为28~32℃,6,胜利原油凝固点为23~32℃,50℃运动粘度约为50℃运动粘度约为20~25×s102m/6。
稠油的特点是凝固点很低,通常低于0℃,但粘度很大,如孤岛原油凝80~90×sm/1026。
固点为-2.3~4.9℃,50℃运动粘度约为490×s102m/凝固点(Freezing point):是指在规定条件下(热力和剪切条件)所测得的油样不流动的最高温度。
我国常把它作为评价原油流动性的指标之一。
西方国家常用的是倾点(Pour point),它与凝固点有所不同。
倾点是指在规定条件下测得的油样刚开始流动的最低温度。
由于测量方法的不同,因而两者在数值上亦有差别。
对于同一种原油,倾点一般比凝固点低2~3℃。
原油的高含蜡、高凝固点和高粘度给储运工作带来以下几个方面的问题:1.由于原油的凝固点比较高,一般在环境温度下就失去流动性或流动性很差,因而不能直接常温输送。
2.在环境温度下,含蜡原油既使能够流动其表观粘度(Apparent Viscosity)也很高。
对于稠油,虽然在环境温度下并不凝固,但其粘度很大。
因此无论是高含蜡原油还是稠油,常温输送时摩阻损失都很大,是很不经济的。
3.高凝高粘原油给储运系统的运行管理也带来了某些特殊问题,主要有:①储罐和管道系统的结蜡问题②管道停输后的再启动问题。
对于易凝高粘问题,不能直接采用前面讲到的等温输送方法,必须在输入管道前采用降凝降粘措施。
1管道总传热系数管道总传热系数K 指油流与周围介质温差为1℃时,单位时间内通过管道单位传热表面所传递的热量,它表示油流至周围介质散热的强弱。
当考虑结蜡层的热阻对管道散热的影响时,根据热量平衡方程可得如下计算表达式:(1-1)1112ln 111ln 22i i n e n wi L L D D D KD D D D ααλλ-+⎡⎤⎛⎫ ⎪⎢⎥⎝⎭=+++⎢⎥⎢⎥⎢⎥⎣⎦∑式中:——总传热系数,W/(m 2·℃);K ——计算直径,m ;(对于保温管路取保温层内外径的平均值,对于e D 无保温埋地管路可取沥青层外径);——管道内直径,m ;n D ——管道最外层直径,m ;w D ——油流与管内壁放热系数,W/(m 2·℃);1α ——管外壁与周围介质的放热系数,W/(m 2·℃);2α ——第层相应的导热系数,W/(m·℃);i λi ,——管道第层的内外直径,m ,其中;i D 1i D +i 1,2,3...i n =——结蜡后的管内径,m ;L D ——所结蜡导热系数。
L λ为计算总传热系数,需分别计算内部放热系数、自管壁至管道最外径K 1α的导热热阻、管道外壁或最大外围至周围环境的放热系数。
2α(1)内部放热系数的确定1α放热强度决定于原油的物理性质及流动状态,可用与放热准数、自然1αu N 对流准数和流体物理性质准数间的数学关系式来表示。
r G r P 在层流状态(Re<2000),当时:500Pr <⋅Gr(1-2)1 3.65y dNu αλ==在层流状态(Re<2000),当时:500Pr >⋅Gr(1-3)0.250.330.430.11Pr 0.15Re Pr Pr y y y y y b d Nu Gr αλ⎛⎫==⋅⋅⎪⎝⎭在激烈的紊流状态(Re>104),Pr<2500时:(1-4)0.250.80.441Pr 0.021Re Pr Pr y y yb d λα⎛⎫=⋅⋅ ⎪⎝⎭在过渡区(2000<Re<104)(1-5)25.043.001Pr Pr (Prbf ffd K ⋅λα=式中:——放热准数,无因次;u N ——流体物理性质准数,无因次;λρυC =Pr ——自然对流准数,无因次;()υβw f t t g d Gr -=3——雷诺数;υπρd q vdv4Re ==——系数;)(Re 0f f K =——管道内径,m ;d ——重力加速度,=9.81m/s 2;g g ——定性温度下的流体运动粘度,m 2/s ;υ——定性温度下的流体比热容,J/(kg·K);C ——流体体积流量,m 3/s ;v q ——定性温度下的流体密度,kg/m 3;ρ——定性温度下的流体体积膨胀系数,可查得,亦可按下式计算:β(1-6)tdd-+-=2042045965634023101β——定性温度下的流体导热系数,原油的导热系数约在0.1~0.16f λf λW/(m·K)间,随温度变化的关系可用下式表示:(1-7)153/)1054.01(137.0f t f t ρλ-⨯-=——l5℃时的原油密度,kg/m 3;15f ρ——油(液)的平均温度,℃;f t ——管内壁平均温度,℃;b t ——20℃时原油的相对密度。
1管道总传热系数管道总传热系数是热油管道设计和运行管理中的重要参数。
在热油管道稳态运行方案的工艺计算中,温降和压降的计算至关重要,而管道总传热系数是影响温降计算的关键因素,同时它也通过温降影响压降的计算结果。
1.1 利用管道周围埋设介质热物性计算K 值管道总传热系数K 指油流与周围介质温差为1℃时,单位时间内通过管道单位传热表面所传递的热量,它表示油流至周围介质散热的强弱。
当考虑结蜡层的热阻对管道散热的影响时,根据热量平衡方程可得如下计算表达式:1112ln 111ln 22i i n e n w i L L D 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),当500Pr <⋅Gr 时:1 3.65y d Nu αλ== (1-2) 在层流状态(Re<2000),当500Pr >⋅Gr 时: 0.250.330.430.11Pr 0.15Re Pr Pr y y y y y b d Nu Gr αλ⎛⎫==⋅⋅ ⎪⎝⎭ (1-3)在激烈的紊流状态(Re>104),Pr<2500时:0.250.80.441Pr 0.021Re Pr Pr y y y b d λα⎛⎫=⋅⋅ ⎪⎝⎭ (1-4)在过渡区(2000<Re<104)25.043.001)Pr Pr (Pr b ff f d K ⋅λα= (1-5)式中:u N ——放热准数,无因次;λρυC =Pr ——流体物理性质准数,无因次; ()υβw f t t g d Gr -=3——自然对流准数,无因次; υπρd q vd v 4Re ==——雷诺数; )(Re 0f f K =——系数;d ——管道内径,m ;g ——重力加速度,g =9.81m/s 2;υ——定性温度下的流体运动粘度,m 2/s ;C ——定性温度下的流体比热容,J/(kg·K); v q ——流体体积流量,m 3/s ;ρ——定性温度下的流体密度,kg/m 3;β——定性温度下的流体体积膨胀系数,可查得,亦可按下式计算:t d d -+-=2042045965634023101β (1-6)f λ——定性温度下的流体导热系数,原油的导热系数f λ约在0.1~0.16W/(m ·K)间,随温度变化的关系可用下式表示:153/)1054.01(137.0f t f t ρλ-⨯-= (1-7)15f ρ——l5℃时的原油密度,kg/m 3;f t ——油(液)的平均温度,℃;b t ——管内壁平均温度,℃;204d ——20℃时原油的相对密度。