一种多炉一注树枝状注汽管网分解计算方法与流程

文档序号:11155789阅读:398来源:国知局
一种多炉一注树枝状注汽管网分解计算方法与制造工艺
本发明属于石油测井
技术领域
,具体涉及一种多炉一注树枝状注汽管网分解计算方法。
背景技术
:稠油热采区块普遍采用蒸汽吞吐、蒸汽驱、SAGD等开发方式,根据原油的黏度随温度的升高而降低的特性,利用注入蒸汽所携带的热量加热稠油和地层,进而降低原油黏度增加其流动性,达到稠油开采的目的,目前注汽方式主要有一炉一注、多炉一注和一炉多注等,对应的注汽管网有单一注汽管线、星状管网、树枝状管网等,在计算地面管线热损失的过程中,目前只有一炉一注单一管线的计算方法,即从注汽锅炉到井口只有一条管线,管线不存在分支及汇入情况,但大多数注汽是采用星状或树枝状的管网结构,即多个注汽锅炉向主管线注汽,然后在分支注向单井,因此单一注汽管线热损失计算方法已不再适用。技术实现要素:为解决上述问题,本发明的目的是提供一种多炉一注树枝状注汽管网分解计算方法。为达到上述目的,本发明提供了一种多炉一注树枝状注汽管网分解计算方法,该方法包括以下步骤:S1:将树枝状注汽管网分解为若干单元沿着树枝状注汽管网的注气方向,以分支管线与主管线的汇合点作为分解点,将树枝状注汽管网分解成若干子单元;所述若干单元依次为第一子单元、第二子单元、……以及末端子单元;在所述若干单元中,末端子单元为一段主管线,其余的子单元均包括一段主管线以及与该段主管线末端相连的分支管线;S2:分别计算各子单元的干度和热损失分别对所述若干子单元的干度和热损失进行计算,其中,前一个子单元汇合点处的温度、压力、干度和流量为后一个子单元主管线段的起始参数;S3:根据步骤S2的结果获得整个管网的干度、热损失分布以及管网终点处的参数。本发明提供一种多炉一注树枝状注汽管网分解计算方法能将复杂的树枝状注汽管网结构分解成简单的分段计算模型,然后逐段计算,最终计算出整个管网的热损失及干度分布数据,以满足生产过程中井口干度的初判、管线保温效果评价的需求,及时采取措施,保证注汽效果。在上述一种多炉一注树枝状注汽管网分解计算方法中,末端子单元为一段主管线,属于一炉一注的简单计算模型,目前有已知的计算方法;包括一段主管线以及与该段主管线末端相连的分支管线的子单元属于爪状模型,目前尚未有针对该模型的计算方法。在上述一种多炉一注树枝状注汽管网分解计算方法中,优选地,在所述步骤S2中,对所述包括一段主管线以及与该段主管线末端相连的分支管线的子单元的干度和热损失进行计算的步骤包括:分别计算该段主管线和分支管线的干度和热损失;计算该段主管线和分支管线汇合点处的混合干度和压力值;输出汇合点处的温度、压力、干度和注汽量。在上述一种多炉一注树枝状注汽管网分解计算方法中,优选地,在所述计算该段主管线和分支管线汇合点处的混合干度和压力值的步骤中,计算的基本假设条件为:根据质量守恒原理,汇合点处的总注汽量等于汇合点处相连接各管线注汽量的总和,各个锅炉注入的蒸汽经过各段管线的压降后,在汇合点处的蒸汽的气相质量之和应等于主管线在此位置的气相质量。在上述一种多炉一注树枝状注汽管网分解计算方法中,优选地,在所述计算该段主管线和分支管线汇合点处的混合干度和压力值的步骤中,还包括计算汇合点处的局部阻力损失的步骤;在综合汇合点处的局部阻力损失后,再输出汇合点处的温度、压力、干度和注汽量。在上述一种多炉一注树枝状注汽管网分解计算方法中,优选地,对于由多段管道连接形成的主管线段或分支管线,采用逐步计算法计算主管线段或分支管线的干度和热损失,所述逐步计算法包括以下步骤:以实际注汽管网中的管道连接处为作为节点对管线进行分段;获取管线的入口参数和管线的基础参数;以管线的入口处为起点,采用迭代计算法逐段计算各段管道出口处的温度、压力、蒸汽干度以及热损失,直到最后一根管道,最终获得管线出口处的温度、压力、干度以及热损失。在上述一种多炉一注树枝状注汽管网分解计算方法中,优选地,当管线的入口处为锅炉时,管线的入口参数为锅炉出口参数;所述锅炉出口参数包括:每个锅炉的出口温度Tb、每个锅炉的出口压力Pb、每个锅炉的出口蒸汽干度Fb和每个锅炉的注汽量G;当管线的入口处为上一个子单元的汇合点处时,管线的入口参数为上一个子单元输出的温度、压力、干度和注汽量;所述管线的基础参数包括:每根管道类型、每根管道是否有保温层、每根管道处的空气温度Ta(i)、每根管道处的风速νa(i)、每根管道/阀门的导热系数λp(i)、每根管道/阀门内径ri(i)、每根管道/阀门外径ro(i)、每根管道/阀门长度z(i)、每根管壁/阀门外黑度ε(i)、每根管道内壁表面粗糙度Ra(i)、每根管道倾角θ(i)、每根管道修正系数、每根管道热阻修正系数、每根管道/阀门保温层导热系数λil(i)和每根管道保温层厚度Pipe(i,20)。在上述一种多炉一注树枝状注汽管网分解计算方法中,优选地,所述以管线的入口处为起点,采用迭代计算法逐段计算各段管道出口处的温度、压力、蒸汽干度以及热损失,直到最后一根管道,最终获得管线出口处的温度、压力、干度以及热损失的步骤包括:步骤一:以管线的入口为计算起点,计算第i段管道的相关参数;其中,入口处的管道编号为1,依次类推,管道编号i=1-n,管线的入口参数即为第1段管道的入口参数,按照管道的连接顺序依次计算,i=1时开始计算;步骤二:迭代计算第i段管道压力损失,同时计算第i段管道局部阻力损失;迭代计算第i段管道外表面温度和热损失;迭代计算第i段管道末端的蒸汽干度;输出第i段末端的相关参数,所述相关参数包括温度、压力、蒸汽干度、热损失;步骤三:按照计算第i段管道的蒸汽干度的方法计算后续管道的相关参数,直到最后一根管道,最后一根管道末端的温度、压力、蒸汽干度、热损失即为管线出口处的相关参数。本发明提供了多炉一注树枝状管网的分解方法及爪状计算模型的计算方法,利用该方法可以将复杂的管网结构简化成若干简单计算模型,然后采用节点计算的方法分别对各计算模型进行计算,最终得到整个管网的干度、热损失分布情况,根据计算结果可以直观的了解整个注汽管线的热损失状况,进而采取相应的保温措施、调整注汽参数、调整注汽方案,有效降低注汽运行成本。附图说明图1为实施例1中多炉一注树枝状管网结构示意图。具体实施方式为了对本发明的技术特征、目的和有益效果有更加清楚的理解,现对本发明的技术方案进行以下详细说明,但不能理解为对本发明的可实施范围的限定。实施例1本实施例提供了一种多炉一注树枝状注汽管网分解计算方法,该方法的具体方案是:S1:将树枝状注汽管网分解为四个单元多炉一注树枝状管网如图1所示,沿着树枝状注汽管网的注气方向,整个管线分成主管线和分支管线两部分组成,本着模型统一、简单、可计算的原则,可以将整个管网划分成第一子单元、第二子单元、第三子单元和末端子单元,其中,第一子单元、第二子单元、第三子单元为爪状模型,末端子单元为一炉一注的简单计算模型。S2:分别计算各子单元的干度和热损失分别对四个子单元的干度和热损失进行计算,其中,前一个子单元汇合点处的温度、压力、干度和流量为后一个子单元主管线段的起始参数;具体地,将第一子单元作为计算起点,输入两个起点的温度、压力、干度和注汽量,首先计算出第一子单元的干度及热损失分布,并输出交汇点1处的温度、压力、干度和流量;输入起点2的温度、压力、干度和流量,计算得到第二字单元的热损失和干度分布,并输出交汇点2处的温度、压力、干度和流量,然后输入起点3的温度、压力、干度和流量,计算得到第三子单元的热损失、干度分布,并输出交汇点3处的温度、压力、干度和流量;最后采用一炉一注单一管线的计算方法计算模型四的热损失和干度分布,最后输出终点的温度、压力、干度和注汽量。S3:根据步骤S2的结果获得整个管网的干度、热损失分布以及管网终点处的参数。在上述步骤S2中,在对属于爪状模型的第一子单元、第二子单元、第三子单元的干度、热损失进行计算时,汇合点处满足如下假设:根据质量守恒原理,汇合点的总注汽量等于各支路注汽量的总和,各个锅炉注入的蒸汽经过各段管线的压降后,在交汇点处的蒸汽的气相质量之和应等于主管线在此位置的气相质量,即:在此基础上,计算爪状模型的干度、热损失的步骤可概括为:(1)输入爪状模型中各管线的入口参数和管线的基础参数,当管线的入口处为锅炉时,管线的入口参数为锅炉出口参数;各参数见表3:(2)当Num_L=1时,计算对应管线第i段管道的相关参数;(3)迭代计算Num_L=1对应管线第i段管道压力损失;(4)计算Num_L=1对应管线第i段管道局部阻力损失;(5)迭代计算Num_L=1对应管线第i段管道外表面温度和热损失;(6)迭代计算Num_L=1对应管线第i段管道末端的干度;(7)计算Num_L=1对应管线第i+1段管道的干度;(8)判断Num_L=1对应管线是否计算到节点处,如果计算完毕,保存计算得到的温度、压力、干度、注汽量;如果没有则继续计算下一管道的参数;(9)计算下一管线Num_L=1+Num_L对应管线的参数,直至所有管线计算完毕,每计算完一根管线保存终点处的温度、压力、干度和注汽量;(10)计算管线汇合点处的混合干度和压力值;表3输入参数列表序号名称符号单位备注1每个锅炉的出口温度Tb℃输入2每个锅炉的出口压力PbMPa输入3每个锅炉的出口干度Fb%输入4每个锅炉的注汽量Gt/h输入5每根管道连接上一段管道编号输入6每根管道的编号输入7每根管道的类型输入8每根管道是否有保温层输入9每段管道处的空气温度Ta(i)℃输入10每段管道处的风速νa(i)m/s输入11每段管道/阀门的导热系数λp(i)W/(m.K)输入12每段管道/阀门内径ri(i)m输入13每段管道/阀门外径ro(i)m输入14每段管道/阀门长度z(i)m输入13每段管壁/阀门外黑度ε(i)输入14每段管道内壁表面粗糙度Ra(i)输入15每段管道倾角θ(i)°输入17每段管道修正系数输入18每段管道热阻修正系数输入20每段管道/阀门保温层导热系数λil(i)W/(m.K)输入22每段管道保温层厚度m输入23锅炉个数Num_B输入24管道数Num_P输入25井数Num_W输入26管线数Num_L输入(11)计算汇合点处的局部阻力损失;(12)输出汇合点处的温度、压力、干度和注汽量。在一种具体的实施方式中,上述计算爪状模型的干度、热损失的步骤具体为:步骤1、分别对子单元中的主管线和分支管线所涉及的单井注汽量G、管道的导热系数λp和绝热层导热系数λins进行单位换算;并分别计算第一子单元中的主管线和分支管线各段管道的空气导热系数λa和空气运动黏度υa;Gt/h=注汽量×1000/3600kg/sλp=0.859845×管线导热系数kcal/(h·m·℃)λins=0.859845×绝热层导热系数kcal/(h·m·℃)λa=(9×10-18×Ta6-3×10-14×Ta5+4×10-11×Ta4-2×10-8×Ta3+2×10-6×Ta2+0.0077×Ta+2.4313)×10-2υa=(3×10-16×Ta6-9×10-13×Ta5+9×10-10×Ta4-4×10-7×Ta3+0.0002×Ta2+0.0862×Ta+13.232)×10-6λa=0.859845×空气导热系数kcal/(h·m·℃)步骤2、分别计算子单元中的主管线和分支管线各段管道的R3;在上述公式中,R3为管线的管壁热阻;λp为管线的导热系数,kcal/(h·m·℃);ro为管线外半径,m;ri为管线内半径,m。步骤3、分别计算子单元中的主管线和分支管线各段管道对空气的对流换热系数hfc';Re=νaDs/υa上述公式中,λa为空气的导热系数,kcal/(h·m·℃);Re为雷诺数;νa为风速,m/s;υa为空气的运动粘度,m2/s;Ds为保温层外径,m;C和n根据Re按照表2的标准进行选值;表2C和n的选取标准Re5-8080-5×1035×103-5×104>5×104C0.810.6250.1970.023n0.400.460.60.8步骤4、分别计算子单元中的主管线和分支管线各段管道的管子内截面积AA=πri2在上述公式中,A为管线的内截面积,m2。步骤5、根据有无保温层情况分别计算子单元中的主管线和分支管线各段管道的外半径或者保温层外半径rins有保温层时:采用以下公式计算保温层外半径:rins=ro+Pipe(i,20)Ds=2rins在上述公式中,rins为管线的外半径,m;ro为管线外半径;Ds为保温层外径,m;无保温层时,采用以下公式计算管线外半径:rins=ro。步骤6、子单元中的主管线和分支管线各段管道分别取一段长度△z,干度降△xi。步骤7、子单元中的主管线和分支管线各段管道分别取压力降△pi。步骤8、计算子单元中的主管线和分支管线各段管道出口处压力pi、出口处温度Ti以及出口处蒸汽干度xi;pi=pi-1-△piTi=195.94pi0.225-17.8xi=xi-1-△xi。步骤9、计算子单元中的主管线和分支管线各段管道的平均压力pavi、平均温度Tavi以及平均蒸汽干度xavipavi=(pi-1+pi)/2Tavi=(Ti-1+Ti)/2xavi=(xi-1+xi)/2。步骤10、计算子单元中的主管线和分支管线各段管道的蒸汽汽相密度ρg和蒸汽液相密度ρl;采用以下经验关系式计算:ρl=(0.9967-4.615×10-5Tavi-3.063×10-6Tavi2)×103Zg=1.012-4.461×10-4Tavi+2.98×10-6Tavi2-1.663×10-8Tavi3上述公式中,ρl为第i段管道的蒸汽液相密度,kg/m3;ρg为第i段管道的蒸汽汽相密度,kg/m3。步骤11、计算子单元中的主管线和分支管线各段管道的蒸汽液相密度μl和蒸汽汽相密度μg;μg=(0.36Tavi+88.37)×10-4上述公式中,μl为第i段管道的蒸汽液相密度,mPa.s;μg为第i段管道的蒸汽汽相密度,mPa.s。步骤12、计算子单元中的主管线和分支管线各段管道的体积含气率Hg;步骤13、计算子单元中的主管线和分支管线各段管道的平均密度ρm;ρm=Hgρg+(1-Hg)ρl。步骤14、计算子单元中的主管线和分支管线各段管道的平均密度ρm;μm=Hgμg+(1-Hg)μl。步骤15、计算子单元中的主管线和分支管线各段管道的平均流速νm;上述公式中,νm为第i段管道的平均流速,m/s。步骤16、计算子单元中的主管线和分支管线各段管道的雷诺数Re步骤17、计算子单元中的主管线和分支管线各段管道的摩擦系数fm;fm的具体取值根据表1的标准选取;表1第i段管道的摩擦系数选取标准表1中,Ra为管壁粗糙度。步骤18、计算子单元中的主管线和分支管线各段管道pi-1和Ti-1下的蒸汽汽相密度ρgi-1和蒸汽液相密度ρli-1,以及体积含气率Hgi-1、平均密度ρmi-1和流速νi-1;ρli-1=(0.9967-4.615×10-5Ti-1-3.063×10-6Ti-12)×103Zgi-1=1.012-4.461×10-4Ti-1+2.98×10-6Ti-12-1.663×10-8Ti-13ρmi-1=Hgi-1ρgi-1+(1-Hgi-1)ρli-1步骤19、计算子单元中的主管线和分支管线各段管道pi和Ti下的蒸汽汽相密度ρgi和蒸汽液相密度ρli,以及体积含气率Hgi、平均密度ρmi以及流速νi;pi=pi-1-△piTi=195.94pi0.225-17.8ρli=(0.9967-4.615×10-5Ti-3.063×10-6Ti2)×103Zgi=1.012-4.461×10-4Ti+2.98×10-6Ti2-1.663×10-8Ti3ρmi=Hgiρgi+(1-Hgi)ρli步骤20、计算子单元中的主管线和分支管线各段管道的局部阻力△pj上述公式中,△pj为第i段管道局部压力降,MPa;△pw为液体单向流的局部压强损失,Pa;X为马蒂内利参数;ζ为局部阻力系数;当管径扩大时,B1和ζ的计算公式为:B1=1.0,当管径缩小时,B1和ζ的计算公式为:B1=1.0,当存在90°弯头时,B1和ζ的计算公式为:ζ=0.12;当存在阀门时,B1和ζ的计算公式为:闸阀:B1=1.5,ζ=0.2,球阀:B1=2.3,ζ=10.0,控制阀:B1=1.0,ζ=5.0s,在上述涉及B1和ζ的计算公式中,l为管子弯头部分的长度;A2为下游小管道的截面积;A1为上游大管道的截面积。步骤21、计算子单元中的主管线和分支管线各段管道的压降△pi';△pi'=pi-1-pi。步骤22、判断计算得到的△pi'与△pi,如果在误差范围内则进行下步计算,否则取△pi=△pi'返回步骤6;△pi'与△pi的误差范围可具体根据使用者对精度要求设定,在本实施例中设置为万分之一。步骤23、假定子单元中的主管线和分支管线各段管道的外表面温度Tw。步骤24、计算子单元中的主管线和分支管线各段管道对空气的强迫对流换热的热阻R5;绝热层外表面由于通过强迫对流方式与大气进行热交换,形成低速气膜层,其热阻可表示为:式中hfc为绝热层外表面上强迫对流热系数,kcal/(m2·h·℃)。其中hfc=hfc'+hfc",即hfc由对流换热系数hfc'和辐射换热系数hfc"组成;管外壁至大气的辐射换热系数hfc"由下式计算:在上述公式中,ε为管壁外黑度;Ta为空气平均温度,℃;Tw为绝热层外壁温度,℃。步骤25、根据有无保温层情况计算子单元中的主管线和分支管线各段管道的保温层热阻R4;并计算子单元中的主管线和分支管线各段管道△z上的单位长度、单位时间的热损失q;(1)有保温层时计算绝热层热阻在上述公式中,λins为绝热层的导热系数,kcal/(h·m·℃);rins为绝热层外半径,m。(2)无保温层时步骤26、计算子单元中的主管线和分支管线各段管道的绝热层外表面温度Tw'(1)有保温层时(2)无保温层步骤27、判断Tw和Tw',如果误差较大则Tw=Tw'返回步骤21,如果满足则进行下一步计算;Tw和Tw'的误差可具体根据使用者对精度要求设定,在本实施例中,设置为万分之五。步骤28、计算子单元中的主管线和分支管线各段管道热损失qi;qi=q×z×4.186/3600/G单位kJ/kg。步骤29、计算子单元中的主管线和分支管线各段管道热流密度Q和累计长度Z管线有保温层时,管线无保温层时,累计长度Z为从锅炉出口出至第i段管道的管线长度。步骤30、计算子单元中的主管线和分支管线各段管道的饱和蒸汽的焓hg以及饱和水的焓hl;并计算子单元中的主管线和分支管线各段管道的蒸汽干度xi;hg=(12500+1.88Tavi-3.7×10-6Tavi3.2)/4.186在上述公式中,hl为饱和水的热焓,kcal/kg;hg为饱和蒸汽的热焓,kcal/kg。C1=G(hg-hl)步骤31、判断子单元中的主管线和分支管线各段管道的(xi-xi-1)/△xi是否小于误差允许范围,如果不满足则△xi=xi-xi-1返回步骤5重新计算,满足则进行下一△z管道的计算,直到每个管道长度;误差可具体根据使用者对精度要求设定,在本实施例中,设置为万分之一。步骤32、计算每个节点的累计热损失量qi:qi=qi-1+q×△z×4.186/3600;步骤33、计算汇合点处主管线和分支管线计算得到压力的平均值:(1)算术平均值(2)几何平均值步骤34、计算汇合点处的干度值:根据质量守恒原理,各个锅炉注入的蒸汽经过各段管线的压降后,在交汇点处的蒸汽的气相质量之和应等于主管线在此位置的气相质量,即:步骤35、计算汇合点处各锅炉管线主管线和分支管线计算结果压力的最大值pmax和最小值pmin;步骤36、判断主管线和分支管线计算结果压力的最大值pmax和最小值pmin与平均压力pavi,如果满足要求则进行下一步计算,如果不满足,若pi大于pavi,加大对应管线的管线修正系数值,若pi小于pavi,减小对应管线的管线修正系数值,然后返回步骤1重新计算,直到pmax-pavi和pmin-pavi绝对值的最大值满足误差允许范围;步骤37、计算汇合点处的局部阻力△pj;步骤38、输出汇合点处的压力、温度、干度、总注汽量。在设置完成上述步骤1至步骤38的基础上,对某多炉一注树枝状注汽管网进行实际计算,其中,输入的锅炉出口参数见表4;本次计算取三个锅炉管线的参数相同,锅炉管线的基础参数见表5;锅炉1对应管线的计算参数见表6,锅炉2对应管线的计算参数见表7,锅炉3对应管线的计算参数见表8;井口计算结果见表9。表4输入的注汽参数名称锅炉1锅炉2锅炉3单位锅炉出口温度324331315℃锅炉出口压力11.412.910.5MPa锅炉出口干度0.7450.7510.755%注汽量44.56t/d表6锅炉1对应管线的计算参数长度压力温度干度热损失管线外壁温度热流密度mMPa℃%kJ/kg℃W/m^2011.43240.7450001.411.236319.8840.7454.748321.13211107.192.511.104318.990.7448.426318.63810953.824.710.835317.1350.7449.19959.514501.2269.710.189312.5350.74310.93759.117495.73714.79.489307.2860.74312.66455.839492.54616.19.282305.6790.74213.1645.642505.66817.59.073040.74213.65345.481502.79518.98.852302.2430.74114.14445.312499.79220.38.628300.40.74114.63145.136496.64421.78.396298.4610.7415.11644.95493.337表7锅炉2对应管线的计算参数长度压力温度干度热损失管线外壁温度热流密度mMPa℃%kJ/kg℃W/m^2012.93310.7510001.412.669329.1270.7514.413329.21611614.582.512.483327.9760.757.851327.71111518.934.712.1325.5620.758.55960.579516.0399.711.164319.3980.74910.14460.055508.74114.710.117312.0070.74811.70956.56502.37716.19.801309.660.74812.15746.068513.27717.59.472307.1570.74712.60145.83509.03618.99.13304.4760.74613.04145.577504.50420.38.772301.5890.74513.47745.304499.63621.78.396298.460.74413.90945.008494.377表8锅炉3对应管线的计算参数长度压力温度干度热损失管线外壁温度热流密度mMPa℃%kJ/kg℃W/m^2010.53150.7550001.410.379313.9120.7553.036313.67910653.192.510.283313.220.7545.409312.79310600.074.710.088311.7950.7545.91558.831491.789.79.626308.3370.7537.05458.53487.62914.79.136304.5250.7538.19155.405486.64316.18.994303.3870.7528.51945.394501.2517.58.848302.2130.7518.84645.281499.22918.98.701301.0020.759.17145.164497.14520.38.55299.7520.7499.49545.043494.99421.78.396298.460.7489.81844.919492.771表9井口计算结果压力温度干度注汽量MPa℃%t/h8.396298.460.74514.5表9列出了该爪状模型最终汇合处井口的计算结果。当前第1页1 2 3 
当前第1页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1