技术领域
[0001] 本
发明属于石油火驱开采技术领域，尤其涉及一种重力火驱中
水平生产井混合液的
温度计算、控制方法及装置。
背景技术
[0002] 在石油开采领域，火驱开采法是一种常用的且非常重要的石油开采方法。该方法通常是指将空气或
氧气注入到地下油层里面，使其与油层中的有机
燃料燃烧，用燃烧产生的热量加热待开采
原油，加热后原油在重力的驱动下卸到水平生产井的水平段，然后通过水平井的提升作用将石油开采出来的方法。
[0003] 火驱开采法是一种在油层内部产生热量的热力采油技术，在石油开采领域具有明显技术优势和潜力，它具有驱油效率高(一般达80％～90％)、单位成本低(与
蒸汽相当)、油藏适应范围广(从薄油层到厚油层、从浅油层到深油层、从稀油到稠油，及已开发油藏)等特点，是石油开采尤其是稠油开采的重要手段。常用的火驱开采方式是直井+水平井组合。直井为注气井，用于将燃烧产生的热量向前推进加热待开采的原油。加热后并达到一定流动性的原油在重力的驱动作用下卸到水平生产井的水平段，原油通过水平段的筛管空隙进入油管和筛管的环空之间，最后通过生产管柱将原油举升到地面。由于目前火驱开采主要针对
粘度较大的稠油、超稠油，因此在水平生产井内原油的温度一旦降低到一定
阀值，原油流动性就会变得非常差，对石油举升造成很大困难。水平井内原油的温度是后续传输过程中最重要的参数，直接影响到整个火驱采油效率及能耗
[0004] 在现场石油开采中常常需要对水平生产管的原油进行伴热，使其达到一定的温度，以保证原油在高温状态下的良好流动性。
现有技术中传统的油井测温装置采用的是
套管式温度测量，这种测试装置安装和使用都十分费时费力，且使用和维护成本较高。同时，在测定管道中
流体温度时，为了保护温度计，通常会在管道中装有测温套管，导致该测试装置由于存在
传热热阻，导致测量温度的准确性不高。
发明内容
[0005] 本
申请目的在于提供一种重力火驱生产井混合液的温度计算、控制方法及装置，能根据特定采集的参数及
算法计算得出生产井内原油的温度，并可以根据计算得出的温度进行相应的调控，降低垂直生产井内原油举升难度，同时降低使用和维护成本，省时省力。
[0006] 本申请提供的一种重力火驱生产井混合液的温度计算、控制方法及装置是这样实现的：
[0007] 一种重力火驱生产井混合液的温度计算、控制方法，所述方法包括：
[0008] S1：计算生产井井筒
热损失的上限值；
[0009] S2：根据所述井筒热损失的上限值设置脚跟处井筒分段单元径向热损失初始值；
[0010] S3：根据所述设置的热损失初始值
迭代计算得出所述生产井井筒从脚跟到脚尖的第一温度场分布；
[0011] S4：根据所述第一温度场分布计算脚尖油管出口处热水温度T2，并利用下述公式计算所述脚尖处油水混合液的初始温度T：
[0012] LvΔm汽+C水Δm水T2+C油Δm油T1＝CΔmT
[0013] 上述公式中：Lv为
汽化潜热，J/kg；C＝C汽fw+C油(1-fw)，C汽为蒸汽
比热容，J/(kg·℃)；C油为原油比
热容，J/(kg·℃)；Δm油为蒸汽与原油混合瞬间原油
质量，kg；Δm汽为蒸汽与原油混合瞬间蒸汽质量，kg；m为水、蒸汽与原油后混合液质量；T1为重力火驱卸油后流到脚尖处原有的初始温度，℃；T2为油管脚尖出口处热水温度，℃；fw为产出液含水率；
[0014] S5：根据所述第一温度场分布和所述脚尖处油水混合液的初始温度利用下述公式从所述井筒脚尖开始计算所述脚尖至脚跟的第二温度场分布：
[0015]
[0016] 上述公式中，Tti为油管内壁温度，Tto为油管内壁温度，Q液为油水混合液热损失，ρ为油水混合液
密度，ρ＝ρ水fw+ρ油(1-fw)，rci为套管内壁半径，rco为套管外壁半径，R1为蒸汽与管内壁之间的热
对流热阻，R2为油管内外壁之间的热传导热阻，R3为环空热对流热阻，Kcos为常数，T′l为前一段混合液温度，在初始时，所述T′l为所述计算得出的脚尖处油水混合液的初始温度T；
[0017] S6：根据计算得到的生产井混合液的第二温度场分布调节井口的注入参数，控制所述混合液的温度。
[0018] 一种生产井混合液的温度计算装置，所述装置包括：
[0019] 参数获取单元，用于获取计算混合液温度所需的预置的特定参数；
[0020] 热损失上限计算模
块，用于计算生产井井筒热损失的上限值；
[0021] 初始值设置模块，用于根据所述井筒热损失的上限值设置脚跟处井筒分段单元径向热损失初始值；
[0022] 第一温度场分布计算模块，用于根据所述设置的热损失初始值迭代计算得出所述生产井井筒从脚跟到脚尖的第一温度场分布；
[0023] 脚尖初始温度计算模块，用于根据所述第一温度场分布计算脚尖油管出口处热水温度，并计算所述脚尖处油水混合液的初始温度：
[0024] 第二温度场分布计算模块，用于根据所述第一温度场分布和所述脚尖处油水混合液的初始温度利用下述公式从所述井筒脚尖开始计算所述脚尖至脚跟的第二温度场分布：
[0025]
[0026] 上述公式中，Tti为油管内壁温度，Tto为油管内壁温度，Q液为油水混合液热损失，ρ为油水混合液密度，ρ＝ρ水fw+ρ油(1-fw)，rci为套管内壁半径，rco为套管外壁半径，R1为蒸汽与管内壁之间的热对流热阻，R2为油管内外壁之间的热传导热阻，R3为环空热对流热阻，Kcas为常数，T′l为前一段混合液温度，在初始时，所述T′l为所述计算得出的脚尖处油水混合液的初始温度T。
[0027] 本申请所述火驱生产井内油水混合液的计算、控制方法，计算热损失上限值和出脚跟处
指定位置的热阻及温度，采用分段迭代算法计算出从脚跟到脚尖油管的第一温度场分布，然后根据所述第一温度场分布以及由此计算得到的脚尖混合液初始温度反向计算得出所述混合液的第二温度场分布。根据计算得到的生产井混合液的第二温度场分布调节井口的注入参数，控制所述混合液的温度。本申请中以井筒单元径向热流量热损失作为迭代变量，选取略小于井筒单元径向热流量热损失上限值的数值作为迭代初始值，使得迭代过程收敛并且经过多次迭代使得计算结果更加接近真实数据。并且不需要在设置单独的混合液温度测试装置，过滤掉了上述温度测试装置的热阻，并且降低生产运营和维护成本。
附图说明
[0028] 为了更清楚地说明本申请
实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0029] 图1是本申请重力火驱采油水平生产井井筒的剖面结构示意图；
[0030] 图2是本申请重力火驱采油中水平生产井的位置示意图；
[0031] 图3是本申请提供的一种重力火驱生产井混合液的温度计算、控制方法一种实施例的流程示意图；
[0032] 图4是本申请提供的一种重力火驱生产井混合液的温度计算、控制方法另一种实施例的流程示意图；
[0033] 图5是本申请提供的一种重力火驱生产井混合液的温度计算、控制方法另一种实施例的流程示意图；
[0034] 图6是本申请提供的一种生产井混合液的温度计算装置的模块结构示意图；
[0035] 图7是本申请提供的一种生产井混合液的温度计算装置另一种实施例的模块结构示意图；
[0036] 图8本申请提供的一种生产井混合液的温度计算装置另一种实施例的模块结构示意图。
具体实施方式
[0037] 为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
[0038] 图1为本申请重力火驱采油水平生产井井筒的剖面结构示意图。如图1所示，水平段采用普通油管，油管外层是套管(筛管)，所述水平生产井空间位置及作业环境按
地层位置从上至下依次可以为蒸汽、油管、环空层、套管、地层。通常上层加热的蒸汽温度最高，由于井筒中各层之间热传递时存在热损失，在井筒的半径方向从上至下温度逐次降低。如图1所示，蒸汽流经直井段的隔
热管，然后流经水平段的油管，从尾管处流出。原油通过水平段的筛管空隙进入油管和筛管的环空之间。在油管尾部流出的蒸汽与热水与重力驱卸下来的原油混合，混合后的原油即为产出液，最后通过生产管柱将原油举升到地面。图2是本申请重力火驱采油中水平生产井的位置示意图，如图2所示，水平生产井在垂直下井水平拐弯处A端通常称为脚跟，另一端B端通常称为脚尖。
[0039] 下面是本申请一种重力火驱生产井混合液温度的计算方法的一个实施例。图3是所述一种重力火驱生产井混合液温度的计算、控制方法一种实施例的流程示意图，如图3所示所述方法包括：
[0040] S1：计算生产井井筒热损失的上限值。
[0041] 通常在采油过程中，与井筒热损失相关的参数主要包括三个方面：
[0042] 第一种是井筒的结构及其相关的物理性质参数，例如：油管内壁半径rti、油管外壁半径rto、套管内径rci、套管外径rto、套管导热系数、套管长度l等。
[0043] 第二种是地层的热物性参数，例如：地表温度Te、地层导热系数Ke等。
[0044] 第三种是井口注入参数，例如：注汽压力Ps、注汽干度Xi、注汽速率Ms、注汽时间t等。
[0045] 本实施例中可以采用下述方法计算所述井筒热损失的上限值：
[0046] S101：计算井筒径向指
定位置的热阻值，具体的可以包括R1、R2、R3、R4、R5。
[0047] 首先计算井筒各部分的热传导热阻或热对流热阻。通常生产井井筒的热损失是由热对流和热传导引起的，本申请中井筒热损失计算数学模型的详细描述可以参见文献《稠油注蒸汽热采工程》，这里仅列出为获取本申请的技术方案所用到的主要计算公式。在推导模型中可以将水平段管柱图划分成若干个单元，每段单元长度dl，分别计算所述井筒直径方向上指定位置的热阻值。井筒中径向指定的各个部分的热传导或热对流及相应的热阻值计算方法可以依次如下：
[0048] 蒸汽与油管内壁之间的热对流：
[0049]
[0050] 上式(1)中：Q为单位长度dl上的热损
失速度，W；Ts为注入蒸汽温度，℃；Tti为油管内壁温度，℃；R1为蒸汽与管内壁之间的热对流热阻，m·K/W；hf为水膜
传热系数，为已知量，W/(m2·K)；γti为油管内壁半径，m。
[0051] 油管内外壁之间的热传导：
[0052]
[0053] 上式(2)中：Tto为油管外壁温度，℃；R2为油管内外壁之间的热传导热阻，m·K/W；Ktub为油管导热系数，W/(m·K)；rto为油管外壁半径，m。
[0054] 环空中充满油水混合液体时，混合液体的热对流：
[0055]
[0056] 上式(3)中：Tci为套管内壁温度，℃；R3为环空热对流热阻，m·K/W；hc为环空内混合液的强迫对流传热系数，W/(m2·K)。
[0057] 套管壁的热传导：
[0058]
[0059] 上式(4)中：Tco为套管外壁温度，℃；R4为套管壁的热传导热阻，m·K/W；Kcas为套管导热系数，W/(m·K)；rci为套管内壁半径，m；rco为套管外壁半径，m。
[0060] 地层的热传导：
[0061]
[0062] 上式(6)中：Te为地层温度，℃；R5为筛管与地层交界处的热传导热阻，m·K/W；Ke为地层导热系数，W/(m·K)。
[0063] S102：计算地层温度Te。
[0064] 在本实施例中可以采用下述公式计算水平生产井下方地层的地层温度Te：
[0065] Te＝(b0+a1*h/100) (6)
[0066] 上式(7)中：b0为地表温度；℃；a1为
地温梯度，℃/100m；h为直井段井深，m。所述的地温梯度在地球物理勘探领域通常可以指从垂直方向向下，地层垂深每增加100米时相应的地层温度所增加的量。该地温梯度可以是根据作业现场实际的地理环境选取的经验值，是一个已知值。这里所述的直井段井深h可以表示所述水平生产井下方地层的深度。
[0067] S103：计算井筒热损失上限值。
[0068] 上述公式(1)-(7)中，只有公式(3)计算所需的参量环空内混合液的强迫对流传热系数hc为未知量，不能直接计算得到R3的值，其余均可以通过直接计算得到热阻的阻值。
[0069] 在本实施例中，可以将井筒从脚跟到脚尖划分为若干分段单元，可以在井筒的半径方向上将每分段dl长度上井筒分段单元的热损失设为Q，单位为千焦/小时，根据注蒸汽采油的原理，井筒半径方向上各部分Q值满足下式：
[0070]
[0071] 公式(7)中，Ts为注入蒸汽温度；Tti为油管内壁温度；Tto为油管外壁温度；Th为
水泥环外温度；Tco为套管外壁温度；Tci为套管内壁温度；Te为地层温度；dl表示当前井筒分段，米。
[0072] 根据公式(7)可以得到下面两式：
[0073] Tto＝Ts-(R1+R2)Q/dl (8)
[0074] Tci＝Te+(R4+R5)Q/dl (9)
[0075] 将上述(8)、(9)代入公式(3)得到所述井筒热损失上限值Qm：
[0076]
[0077] 上式中，R1为蒸汽与管内壁之间的热对流热阻，R2为油管内外壁之间的热传导热阻，R3为环空热对流热阻，R4为套管壁的热传导热阻，R5为筛管与地层交界处的热传导热阻，Ts为注入蒸汽温度可以直接测量得到，Te为地层温度。
[0078] 一般的情况下，在水平井井口处的热损失为最大，若环空层不存在热损失，即R3的值为0，则井口处的井筒分段单元径向热损失最大。本申请中可以将R3的值为0时的井筒热损失作为井筒分段单元径向热损失的上限值。
[0079] 获取井筒相应参数，计算井筒中半径方向环空部分以外各处的热阻值及地层温度，并根据所述计算得出的热阻值及地层温度计算井筒分段单元径向热损失上限值。
[0080] S2：根据所述井筒热损失的上限值设置脚跟处井筒分段单元径向热损失初始值。
[0081] 井筒分段单元径向热损失的初始值可以用Q0表示。在重力火驱开采中，由于脚跟处燃烧产生的水和蒸汽向脚尖方向流动，可以在油管尾部与重力驱动卸下的原油混合形成油水混合液。所述油水混合液通过油井与套管(筛管)之间的环空从所述脚尖回流到所述脚跟处，然后可以通过垂直生产井将油水混合液
抽取到地面。基于本实施例所述的生产方式，在实施例中可以设置一个脚跟处井筒的热损失初始值。由于井筒分段单元径向热损失的值一定不大于上限值Qm，因此井筒分段单元径向热损失的初始值Q0的取值为0～Qm。例如本实施例中可以根据实验的经验可以取Q0的值为0.9Qm。
[0082] S3：根据所述设置的热损失初始值迭代计算得出所述生产井井筒从脚跟到脚尖的第一温度场分布。
[0083] 该实施步骤中可以将所述水平生产井井筒总长度L划分为若干段，每一分段可以用dl表示。本实施例中从油管脚跟分段处开始计算每一分段的温度场分布数据，所述的dl可以理解为表示为计算完当前分段的温度场分布数据后需要当前分段增加dl计算的下一分段的温度场分布的预定步长。所述预定步长dl的取值范围可以为：0≤dl≤L。本实施例中可以分别用k和Qk表示第k次循环和第k次迭代后的井筒热损失，然后采用下述方法迭代计算所述油管从脚跟到脚尖的温度场分布数据：
[0084] S301：分别用k和Qk表示第k次循环和第k次迭代后的井筒热损失，所述生产井井筒脚跟处开始，令所述l＝0，k＝1，Qk＝NQm，0≤N≤1，如上所述，N可以取值为0.9，即井筒分段单元径向热损失的初始值用Q0可以为0.9倍的热损失上限值Qm。获取所述井筒脚跟处指定位置的热阻值和地层温度Te。
[0085] 上述所述指定位置的热阻值包括R1、R2、R3、R4、R5，所述R1为蒸汽与管内壁之间的热对流热阻，R2为油管内外壁之间的热传导热阻，R3为环空热对流热阻，取值为0，R4为套管壁的热传导热阻，R5为筛管与地层交界处的热传导热阻，其中环空热对流热阻R3初始取值为0。
[0086] S302：利用下述公式计算所述当前分段的蒸汽的平均温度Ts：
[0087] Ts＝195.94P0.225-17.8 (11)
[0088] 上式公式(12)中P为该段的平均压力。在本实施例中可以采用Beggs-Bill算法(具体计算详见张琪主编《采油工程原理与设计》P50)计算获取所述平均压力P：
[0089]
[0090] 上述公式(13)中：p──当前分段平均压力(绝对)，Pa；z──轴向流动的距离，m；ρl──液相密度，kg/m3；ρg──气相密度，kg/m3；Hl──持液率，m3/m3；g──重力
加速度，m/s2；θ──管道与水平方向的夹
角，°；λ──两相流动的沿程阻力系数，无因次；G──混合物的质量流量，kg/s；v──混合物的流速，m/s；vsg──气相的折算速度，m/s；D──管道直径，m；A──管道截面积，m2。相应的，通过公司(13)计算得出的当前分段的平均压力可以表示为Pk。
[0091] S303：计算当前所述分段下述指定位置的各处温度：
[0092] 油管内壁温度为：Tti＝Ts-R1Qk/dl (13)
[0093] 油管外壁温度为：Tto＝Tti-R2Qk/dl (14)
[0094] 筛管与地层交界处温度为：Tco＝Te+R5Qk/dl (15)
[0095] 筛管内壁温度为：Tci＝Tco+R4Qk/dl (16)
[0096] 利用上述公式(13)～(16)可以计算获取出当前分段井筒预定指定位置的温度数据。
[0097] S304：重新计算环空热对流热阻R3。
[0098] 上述中初始环空热对流热阻R3为预先假设的初始值，在此可以基于选取的当前分段热损失Qk和上述公式(15)、(16)计算得到的筛管外壁温度Tco、内壁温度Tci，根据公式(3)重新计算环空热对流热阻R3。
[0099] 本实施例中计算环空热对流热阻R3的公式(3)中的参数环空内混合液的强迫对流传热系数hc的计算公式可以为：
[0100]
[0101] 上述公式(17)中：
[0102]
[0103]
[0104] 其中，
[0105]
[0106] 上述公式(17)～(21)中：Gr为Grashof数；Pr为Prandtl数；Kha为环空流体的导热系数，W/(m·K)；g为
重力加速度，m/s2；ρan为环空流体在平均温度Tan下的密度，kg/m3；ρan为环空流体在平均温度Tan下的粘度，mPas；Can为环空流体在平均温度Tan下的热容，J(m3K)。
[0107] S305：井筒热损失校正。
[0108] 为保证迭代过程的
稳定性和数据的可靠性，本实施例中需要对当前分段井筒热损失进行校正。具体的校正方法可以包括：
[0109] S3051：基于重新计算获取的环空热对流热阻R3利用下式计算总热阻值R：
[0110] R＝R1+R2+R3+R4+R5 (22)
[0111] S3052：基于所述总热阻值R重新计算当前分段的热损失Qk。
[0112] 所述当前分段的热损失Qk的计算公式可以为：
[0113]
[0114] S3053：判断所述重新计算的当前分段的热损失Qk是否满足：Qk≥Qm，如果满足，则令Qk＝Qm。
[0115] S306：利用所述井筒校正后的井筒热损失Qk重新计算S303中所述的各处的温度及相应的热阻值R1、R2、R3、R4、R5。具体的可以将所述校正后的井筒热损失Qk代入公式(13)～(16)中重新计算S303中所述的各处的温度及相应的热阻值。
[0116] S307：令k＝k+1,l＝l+dl，Qk+1＝Qk，Pk+1＝Pk；以所述S306重新计算的所述指定位置的热阻值的作为所述井筒下一分段的初始热阻值计算所述下一分段的总热阻值，并以所述下一分段的总热阻值R作为所述井筒下一分段的初始热阻值利用公式： 计算所述下一分段的热损失Qk；以所述下一分段作为所述当前分段，重复步骤S303～S307进行迭代计算，直到l>L，得到所述水平生产井井筒脚跟至脚尖的第一温度场分布数据；L为所述井筒脚跟至脚尖长度。
[0117] 利用上述过程可以计算得出因蒸汽及由蒸汽形成的热水形成的第一温度场分布数据，在计算过程中迭代变量的修正受控于热损失速度上限值Qm,从而保证了迭代过程的稳定性。同时需要说明的是，在迭代过程中迭代次数k越多，Uto、Qk、Gr、Pr、Tti、Tto等值趋于稳定数值，计算结果更加精确。另外，设置的分段单元或步长dl越小，计算
精度也越高。
[0118] S4：根据所述第一温度场分布计算脚尖油管出口处热水温度T2，并利用下述公式计算所述脚尖处油水混合液的初始温度：
[0119] LvΔm汽+C水Δm水T2+C油Δm油T1＝CΔmT (24)
[0120] Δm＝q*dt (25)
[0121] q＝Q液ρ/24000 (26)
[0122] ρ＝ρ水fw+ρ油(1-fw) (27)
[0123] 上述公式(24)中：Lv为
汽化潜热，J/kg；C汽为蒸汽比热容，J/(kg·℃)；C油为原油比热容，J/(kg·℃)；Δm油为蒸汽与原油混合瞬间原油质量，kg；Δm汽为蒸汽与原油混合瞬间蒸汽质量，kg；m为水、蒸汽与原油后混合液质量；T1为重力火驱卸油后流到脚尖处原有的初始温度，℃；T2为油管脚尖出口处热水温度，℃；fw为产出液含水率；C＝C汽fw+C油(1-fw)；
[0124] 上式公式(25)～(27)中，Δm水＝q汽*dtΔm＝q*dt，Δm油＝Δm-Δm水，q汽为注汽速率，t/h；q为原油流动速率，t/h；Q液为生产井产液量，m3/d；ρ油为脱水原油密度，kg/m3；ρ水为水的密度，kg/m3；ρ为混合液密度，kg/m3。
[0125] 下表为本实施例给出的
饱和蒸汽密度与温度的关系表：
[0126] 表1 饱和蒸汽密度与温度的关系表
[0127]温度范围(℃) 密度对于温度的回归方程(kg/m3)
10～20 0.001(0.02075T2+0.166T+5.674)
20～30 0.01(0.00315T2-0.037T+1.0091)
3 2
30～50 0.001(0.000632T-0.02048T+1.1736T-3.473)
50～70 0.00001(0.0984T3-7.364T2+386.4T-4911)
70～100 0.0001(0.0152T3-1.8996T2+123.202T-2548.87)
100～150 0.0001(0.02487T3-4.9338T2+442.021T-13761.2)
150～200 0.0001(0.04362T3-13.4749T2+1745.043T-80320.5)
200～250 0.0001(0.008611T3-3.95083T2+707.2292T-44435.969)
250～300 0.0001(0.024638T3-16.2298T2+3846.47T-312234.219)
[0128] 相应的可以利用下述公式计算所述脚尖水的密度：
[0129] ρ水＝-0.000003T2-0.000108T+1.000937
[0130] 相应的可以利用下述公式计算所述脚尖原油的密度：
[0131] ρ油＝ρ0+(13.561-0.191x1)×10-3-(63.9-0.87x1)×10-5T1.02
[0132] x1＝INT[100(ρ0-0.8001)]
[0133] 上式中：ρ0为地面20℃时原油密度，g/cm3；INT为取整函数。
[0134] 采用本申请所述方法给出根据对辽河油田进行的火烧现场试验的杜66区块和其它稠油区块的大量现场原油进行取样分析，拟合出的辽河油田稠油比热容与温度关系式：
[0135] C油＝1000(A+BT+CT2+DT3+ET4+FT5) (28)
[0136]样品 A B C D E F
辽河普通稠油 2.5834 -0.0441 0.0007 -4×10-6 1×10-8 -1×10-11
杜66稠油 2.0297 -0.0054 0.0001 -5×10-7 9×10-10 -2×10-13
[0137] 本申请中可以根据上述给出的稠油比热容与温度关系式计算所述为原油比热容C油。
[0138] 另一种实施例中，由于水平生产井
散热较大，通常在所述脚尖油管末端不存在有蒸汽，基本都是热水。在此情况下，LvΔm汽与C汽取值为0，可以采用下述公式计算所述脚尖处油水混合液的初始温度T：
[0139] C水Δm水T2+C油Δm油T1＝CΔmT。
[0140] S5：根据所述第一温度场分布和所述脚尖处油水混合液的初始温度从所述井筒脚尖开始计算所述脚尖至脚跟的第二温度场分布。所述第二温度场分布即为本申请所要计算得到的生产井内油水混合液的温度。
[0141] 根据温度场分布1每段的步长dl从脚尖B端到脚跟A端分步计算每段混合液温度。
[0142]
[0143] 上述公式(29)中Tl为混合液温度，℃；C为混合液比热容，J/(kg·℃)；
[0144] 上式公式(29)
整理后的近似公式为：
[0145]
[0146] 上述公式中，Tti为油管内壁温度，Tto为油管内壁温度，Q液为油水混合液热损失，ρ为油水混合液密度，ρ＝ρ水fw+ρ油(1-fw)，rci为套管内壁半径，rco为套管外壁半径，R1为蒸汽与管内壁之间的热对流热阻，R2为油管内外壁之间的热传导热阻，R3为环空热对流热阻，Kcas为常数，T′l为前一段混合液温度。在初始计算所述第二温度场分布时，所述T′l可以取值为所述计算得出的脚尖处油水混合液的初始温度。
[0147] S6：根据计算得到的生产井混合液的第二温度场分布调节井口的注入参数，控制所述混合液的温度。
[0148] 根据所述生产井混合液的第二温度场分布可以计算得出生产井中各个位置混合液的温度值，可以根据生产井内各个位置的温度值调节相应的井口注入参数，例如可以调节注入
泵的注汽速率或注汽压力，控制所述混合液的温度，以适应作业生产要求。
[0149] 在另一种实施情况中，计算完井筒内蒸汽或水的第一温度场分布后可以得到B点处油管外壁温度Tto的值，此时可以再计算完油管脚尖末端热水与原油混合后混合液温度值Tl。通常所述油管外壁温度Tto的值大于所述油管末端混合液温度值Tl的值，本申请还给出了出现Tto＜Tl的情况时所述第二温度场的计算方法。图4是本申请提供的一种重力火驱生产井混合液温度的计算、控制方法另一种实施例的流程示意图，如图4所示，所述的一种重力火驱生产井混合液的温度计算、控制方法还可以包括：
[0150] S501：比较所述生产井中脚尖处油管的外部温度Tto是否小于计算得出的脚尖处混合液温度Tl，在出现所述出现Tto＜Tl的情况时采用下述公式(32)计算所述第二温度场分布：
[0151]
[0152] 整理后得到：
[0153]
[0154] 本申请提供的一种火驱生产井内油水混合液的计算、控制方法，计算热损失上限值和出脚跟处指定位置的热阻及温度，采用分段迭代算法计算出从脚跟到脚尖油管的第一温度场分布，然后根据所述第一温度场分布以及由此计算得到的脚尖混合液初始温度反向计算得出所述混合液的第二温度场分布。本申请中以井筒单元径向热流量热损失作为迭代变量，选取略小于井筒单元径向热流量热损失上限值的数值作为迭代初始值，使得迭代过程收敛并且经过多次迭代使得计算结果更加接近真实数据。
[0155] 在获取整个所述水平生产井内油水混合液的第二温度场分布后，可以根据需求获取不同段的混合液的温度，可以使作业人员准确掌握生产井内混合液的温度情况。不仅如此，利用本申请的另一个重要的工业应用为可以利用所述第二温度场分布获取脚跟处混合液的温度。在生产作业中可以通过垂直生产井在所述脚跟处将混合液抽取到地面，此时在所述脚跟处泵入口的混合液通常可以称为产出液。为保障将所述产出液顺利的抽取到地面，需要保证所述脚跟泵入口处的产出液具有一定的流动性，达到一定的温度阀值，该温度阀值通常大于100℃。因此可以利用本申请方法控制所述脚跟泵入口处产出液的温度。具体的，图5是本申请提供的一种重力火驱生产井混合液温度的计算方法另一种实施例的流程示意图，如图5所示，本申请所述一种重力火驱生产井混合液的温度计算、控制方法还可以包括：
[0156] S7：利用所述第二温度场分布计算所述脚跟泵入口处产出液的温度；如果所述脚跟泵入口出产出液的温度小于预定温度阀值，则相应的调节井口注入参数(例如注汽速率或注汽压力)以保障所述泵入口产出液的温度不小于所述预定温度阀值。
[0157] 利用本申请所述的生产井混合液的温度计算方法，还可以实时控制泵入口产出液的温度，使产出液具有良好的流动性，使垂直生产井内原油的举升难度降低。并且本申请的监测和控制方法获取常规
测井参数并通过相应的计算即可获取产出液温度，使用简单，维护成本低，省时省力。
[0158] 基于本申请所述的一种重力火驱生产井混合液的温度计算、控制方法，本申请还提供一种生产井混合液的温度计算装置，图6是本申请所述一种生产井混合液的温度计算装置的模块结构示意图，如图6所示，所述装置可以包括：
[0159] 参数获取单元101，可以用于获取计算混合液温度所需的预置的特定参数；
[0160] 热损失上限计算模块102，可以用于根据获取的参数计算生产井井筒热损失的上限值；
[0161] 初始值设置模块103，可以用于根据所述井筒热损失的上限值设置脚跟处井筒分段单元径向热损失初始值；
[0162] 第一温度场分布计算模块104，可以用于根据所述设置的热损失初始值迭代计算得出所述生产井井筒从脚跟到脚尖的第一温度场分布；
[0163] 脚尖初始温度计算模块105，可以用于根据所述第一温度场分布计算脚尖油管出口处热水温度T2，并利用下述公式计算所述脚尖处油水混合液的初始温度T：
[0164] LvΔm汽+C水Δm水T2+C油Δm油T1＝CΔmT
[0165] 上述公式中：Lv为汽化潜热，J/kg；C＝C汽fw+C油(1-fw)，C汽为蒸汽比热容，J/(kg·℃)；C油为原油比热容，J/(kg·℃)；Δm油为蒸汽与原油混合瞬间原油质量，kg；Δm汽为蒸汽与原油混合瞬间蒸汽质量，kg；m为水、蒸汽与原油后混合液质量；T1为重力火驱卸油后流到脚尖处原有的初始温度，℃；T2为油管脚尖出口处热水温度，℃；fw为产出液含水率；
[0166] 第二温度场分布计算模块106，可以用于根据所述第一温度场分布和所述脚尖处油水混合液的初始温度利用下述公式从所述井筒脚尖开始计算所述脚尖至脚跟的第二温度场分布：
[0167]
[0168] 上述公式中，Tti为油管内壁温度，Tto为油管内壁温度，Q液为油水混合液热损失，ρ为油水混合液密度，ρ＝ρ水fw+ρ油(1-fw)，rci为套管内壁半径，rco为套管外壁半径，R1为蒸汽与管内壁之间的热对流热阻，R2为油管内外壁之间的热传导热阻，R3为环空热对流热阻，Kcas为常数，T′l为前一段混合液温度，在初始时，所述T′l为所述计算得出的脚尖处油水混合液的初始温度。
[0169] 上述所述第一温度场分布计算模块104，可以包括：
[0170] 初始化模块，可以用于设置计算所述第一温度场分布的预置参数值；
[0171] 指定温度计算模块，可以用于计算当前所述分段下述指定位置的各处温度：
[0172] 环空热阻重计算模块，可以用于重新计算环空热对流热阻；
[0173] 可以利用下述公式重新计算环空热对流热阻R3：
[0174]
[0175] 热损失校正模块，可以用于基于预定规则校正当前分段井筒热损失；
[0176] 重计算模块，可以用于利用所述井筒校正后的井筒热损失重新计算指定温度计算模块和环空热阻重计算模块所述的指定位置各处的温度及相应的热阻值；
[0177] 迭代模块：可以用于循环迭代计算所述生产井从脚跟至脚尖处的第一温度场分布数据。
[0178] 图7是本申请所述一种生产井混合液的温度计算装置另一种实施例的模块结构示意图，如图7所示，本申请所述的生产井混合液的温度计算装置还可以包括：
[0179] 产出液温度计算模块107，可以用于利用所述第二温度场分布计算所述脚跟泵入口处产出液的温度。
[0180] 进一步的，图8是本申请所述一种生产井混合液的温度计算装置另一种实施例的模块结构示意图，如图8所示，本申请所述的生产井混合液的温度计算装置还可以包括：
[0181] 产出液温度调控模块108，可以用于基于所述计算得到的脚跟泵入口处产出液的温度调节相应的井口注入参数，控制所述产出液的温度。
[0182] 上述各个所述装置中各个模块中所涉及到的具体的计算过程可以参照本申请所述相应的方法，在此不做赘述。
[0183] 利用本申请所述的生产井混合液的温度计算装置可以有效、方便、准确的获取生产井内混合液的温度分布情况，并可以根据生产井内的混合液的温度调控注入参数，对重力火驱石油开采领域带来了很大的方便。
[0184] 本申请所述涉及的方法部分可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
[0185] 上述实施例阐明的装置或模块涉及到的方法部分，具体可以由计算
机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个
软件和/或
硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现。
[0186] 本领域技术人员也知道，除了以纯计算机可读程序代码方式实现
控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以
逻辑门、
开关、专用集成
电路、可编程逻辑控制器和嵌入
微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的
软件模块又可以是硬件部件内的结构。
[0187] 本
说明书中的各个实施例采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。本申请可用于众多通用或专用的
计算机系统环境或配置中。例如：个人计算机、
服务器计算机、多处理器系统、可编程的
电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
[0188] 虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多
变形和变化而不脱离本申请的精神，希望所附的
权利要求包括这些变形和变化而不脱离本申请的精神。
