CN104806230B - 超临界注汽井的井筒温度场及压力场分布的计算方法 - Google Patents

超临界注汽井的井筒温度场及压力场分布的计算方法 Download PDF

Info

Publication number
CN104806230B
CN104806230B CN201510084615.1A CN201510084615A CN104806230B CN 104806230 B CN104806230 B CN 104806230B CN 201510084615 A CN201510084615 A CN 201510084615A CN 104806230 B CN104806230 B CN 104806230B
Authority
CN
China
Prior art keywords
mrow
msub
mfrac
calculate
pit shaft
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510084615.1A
Other languages
English (en)
Other versions
CN104806230A (zh
Inventor
张成博
马振
孙振宇
曲绍刚
刘锦
于晓聪
崔加利
吕孝明
姜佳悦
李�杰
王文刚
卢丽丝
李辉
唐丽
刘佩衡
哈长鸣
康宸博
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201510084615.1A priority Critical patent/CN104806230B/zh
Publication of CN104806230A publication Critical patent/CN104806230A/zh
Application granted granted Critical
Publication of CN104806230B publication Critical patent/CN104806230B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
    • E21B43/24—Enhanced recovery methods for obtaining hydrocarbons using heat, e.g. steam injection

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Geology (AREA)
  • Theoretical Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Geometry (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Measuring Fluid Pressure (AREA)
  • Engine Equipment That Uses Special Cycles (AREA)

Abstract

本发明提供了一种超临界注汽井的井筒温度场及压力场分布的计算方法，涉及油藏开采技术领域。方法包括：获取井筒相关参数以及井口注入参数；获取预先设置的节点步长，将井筒从井口开始划分为多个节点；根据所述井筒相关参数以及井口注入参数，以井口为初始计算点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力；根据所述各节点下端超临界蒸汽温度和超临界蒸汽压力分别确定所述井筒的温度场分布及压力场分布。本能够精确计算超临界注汽井的井筒温度场及压力场。

Description

超临界注汽井的井筒温度场及压力场分布的计算方法
技术领域
本发明涉及油藏开采技术领域，尤其涉及一种超临界注汽井的井筒温度场及压力场分布的计算方法。
背景技术
当前，随着油气田的深入开发，轻、中质原油和浅层稠油区块的产量呈逐年下降的趋势，因此今后深层、超深层稠油油藏将逐渐成为主力开发油藏，但由于其埋深超过2000m，油层压力大于30MPa，通常使用的21MPa亚临界压力蒸汽发生器无法开采这类稠油。根据稠油热采的注汽参数，热采注汽的蒸汽热力学状态主要分为饱和蒸汽状态(包括亚临界注汽锅炉和普通注汽锅炉注汽)以及超临界状态(超临界压力锅炉注汽)。注入油层的热焓影响油层的加热效果，注入热焓越多，油层加热效果越好，油井产量越高。在锅炉发热量相同的情况下，超临界状态比饱和蒸汽状态注入油层的热焓多。如果使用超临界蒸汽发生器产生的高压介质，将注入压力提高到22.14MPa以上，温度提高到374℃以上时，保证蒸汽到达油层以后仍会保持汽相状态，由于超临界蒸汽具有更高的热焓值和穿透性，对油层的加热作用和渗透作用更好，因此能有效解决深层稠油热采的难题。
水通常有三相(固、液、汽)、五态(未饱和水、饱和水、湿饱和蒸汽、干饱和蒸汽和过热蒸汽)，在一般情况下水由液相变为气相都是要经过一个汽化过程，即水经过吸热首先变为饱合水，再经过吸热部分水变为蒸汽，继续吸热后水全部变为蒸汽形成饱合蒸汽，整个汽化要经过一段时间的两相共存过程，并且在湿饱和蒸汽和干饱合蒸汽状态时，增大压力可使水蒸汽重新变为液态。但是当压力大于22.14MPa时，水由液相向汽相的转化没有液、汽两相共存过程，而是在温度升到374℃时，水由液相全部转变为汽相，并且超过此温度后不管再加多大的压力也不能将它变为液相。此时的压力22.14MPa称为临界压力，此时的温度374℃称为临界温度。凡超过此压力、温度的状态，称为超临界状态。这种状态的水即称为超临界水(也称为超临界蒸汽)。临界流体的物性兼具液体与气体性质，其密度比一般气体要大很多，粘度比液体小，但扩散速度比液体快，因而有较好的流动性和传递性能。
当前，在超临界注汽条件下，蒸汽的热力学性质与湿饱和蒸汽差别很大，原有的注汽参数计算模型无法对超临界下井筒的温度场、压力场进行正确计算。
发明内容
本发明的实施例提供一种超临界注汽井的井筒温度场及压力场分布的计算方法，以精确计算超临界注汽井的井筒温度场及压力场。
为达到上述目的，本发明采用如下技术方案：
一种超临界注汽井的井筒温度场及压力场分布的计算方法，包括：
获取井筒相关参数以及井口注入参数；
获取预先设置的节点步长，将井筒从井口开始划分为多个节点；
根据所述井筒相关参数以及井口注入参数，以井口为初始计算点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力；
根据所述各节点下端超临界蒸汽温度和超临界蒸汽压力分别确定所述井筒的温度场分布及压力场分布。
具体的，所述井筒相关参数包括水泥导热系数Kcem；井眼半径rh；隔热管深度h；隔热管导热系数Kins；油管导热系数Ktub；套管导热系数Kcas；环空流体的导热系数Kha；隔热管内管内半径rti；隔热管内管外半径rto；隔热管外管内半径ri；隔热管外管外半径ro；套管内半径rci；套管外半径rco；地温梯度a1；地表温度b0；地层导热系数Ke。
具体的，所述井口注入参数包括井口注汽压力P0；井口蒸汽温度T0；注汽速率qsc；注汽时间tsc。
具体的，所述根据所述井筒相关参数以及井口注入参数，以井口为初始极端点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，包括：
计算蒸汽与隔热管内管内壁之间的热对流热阻R1：
计算隔热管内管内外壁之间的热传导热阻R2：
计算隔热层的热传导热阻R3：
计算隔热管外管壁的热传导热阻R4：
计算环空热对流热阻R5：
计算套管壁的热传导热阻R6：
计算水泥环的热传导热阻R7：
计算地层的热传导热阻R8：
其中，hf为水膜传热系数；hc为环空自然然对流传热系数；hr为环空辐射传热系数；f(t)是随时间变化的导热传热函数。
进一步的，所述根据所述井筒相关参数以及井口注入参数，以井口为初始极端点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
根据所述节点步长、井口蒸汽温度T0、地表温度b0、所述R1、R2、R3、R4、R6、R7、R8计算井口处热损失速度Qm：
其中，dl为所述节点步长。
进一步的，所述根据所述井筒相关参数以及井口注入参数，以井口为初始极端点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
根据所述P0和所述dl计算各节点下端的初始超临界蒸汽压力
其中，Pi为第i个节点上端的超临界蒸汽压力，其中，i≥0。
进一步的，所述根据所述井筒相关参数以及井口注入参数，以井口为初始极端点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
根据所述b0、a1、h确定地层温度Te：
Te＝(b0+a1*h/100)。
进一步的，所述根据所述井筒相关参数以及井口注入参数，以井口为初始极端点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
获取各节点上端的热损失速度Qi，并根据所述Qi、Te、R1、R2、R3、R4、R6、R7、R8、dl依次计算各节点下端的超临界蒸汽温度Ti+1：
Ti+1＝Te+(R1+R2+R3+R4+R6+R7+R8)*Qi/dl。
进一步的，所述根据所述井筒相关参数以及井口注入参数，以井口为初始极端点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
根据所述R1、R2、R3、R4、R5、R6、R7、R8计算一总热阻R：
R＝R1+R2+R3+R4+R5+R6+R7+R8
根据所述Ti+1校正各节点下端的热损失速度Qi+1：
进一步的，所述根据所述井筒相关参数以及井口注入参数，以井口为初始极端点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
计算各节点段的平均温度
计算各节点段的平均压力
进一步的，所述根据所述井筒相关参数以及井口注入参数，以井口为初始极端点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
获取井筒的超临界蒸汽相对密度γg；
计算井筒的拟临界温度Tpc：
Tpc＝92.2+176.6γg
计算井筒的拟对比温度Tpr：
计算井筒的超临界蒸汽粘度μg：
其中，
Ppc＝4.881-0.386γg
计算井筒的超临界蒸汽的体积系数Bg：
计算井筒内的超临界蒸汽流速v：
计算井筒内的超临界蒸汽ρ：
其中，为井筒内的超临界蒸汽在平均温度平均压力下的偏差系数。
进一步的，所述根据所述井筒相关参数以及井口注入参数，以井口为初始极端点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
计算井筒内的超临界蒸汽在平均温度平均压力下的摩阻系数f：
其中，e为绝对粗糙度。
进一步的，所述根据所述井筒相关参数以及井口注入参数，以井口为初始极端点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
计算各节点下端超临界蒸汽压力最终值Pi+1：
其中，
具体的，所述计算各节点下端超临界蒸汽压力最终值Pi+1，包括：
判断的值是否小于0.1；
若的值小于0.1，则确定所述Pi+1准确；
若的值大于等于0.1，则将设置为Pi+1，并重新确定所述
本发明实施例提供的超临界注汽井的井筒温度场及压力场分布的计算方法，通过获取井筒相关参数以及井口注入参数，并获取预先设置的节点步长，将井筒从井口开始划分为多个节点；之后根据所述井筒相关参数以及井口注入参数，以井口为初始计算点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力；从而根据所述各节点下端超临界蒸汽温度和超临界蒸汽压力分别确定所述井筒的温度场分布及压力场分布。可见本发明能够精确计算超临界注汽井的井筒温度场及压力场。解决了在超临界注汽条件下，蒸汽的热力学性质与湿饱和蒸汽差别很大，现有的注汽参数计算模型无法对超临界下井筒的温度场、压力场进行正确计算。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明提供一种超临界注汽井的井筒温度场及压力场分布的计算方法的流程图；
图2为本发明实施例中的超临界注汽井周围结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本发明提供一种超临界注汽井的井筒温度场及压力场分布的计算方法，包括：
步骤101、获取井筒相关参数以及井口注入参数。
步骤102、获取预先设置的节点步长，将井筒从井口开始划分为多个节点。
步骤103、根据所述井筒相关参数以及井口注入参数，以井口为初始计算点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力。
步骤104、根据所述各节点下端超临界蒸汽温度和超临界蒸汽压力分别确定所述井筒的温度场分布及压力场分布。
本发明实施例提供的超临界注汽井的井筒温度场及压力场分布的计算方法，通过获取井筒相关参数以及井口注入参数，并获取预先设置的节点步长，将井筒从井口开始划分为多个节点；之后根据所述井筒相关参数以及井口注入参数，以井口为初始计算点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力；从而根据所述各节点下端超临界蒸汽温度和超临界蒸汽压力分别确定所述井筒的温度场分布及压力场分布。可见本发明能够精确计算超临界注汽井的井筒温度场及压力场。解决了在超临界注汽条件下，蒸汽的热力学性质与湿饱和蒸汽差别很大，现有的注汽参数计算模型无法对超临界下井筒的温度场、压力场进行正确计算。
值得说明的是，该井筒相关参数包括水泥导热系数Kcem；井眼半径rh；隔热管深度h；隔热管导热系数Kins；油管导热系数Ktub；套管导热系数Kcas；环空流体的导热系数Kha；隔热管内管内半径rti；隔热管内管外半径rto；隔热管外管内半径ri；隔热管外管外半径ro；套管内半径rci；套管外半径rco；地温梯度a1；地表温度b0；地层导热系数Ke。如图2所示的超临界注汽井周围的结构，图中分别标注出了隔热管内管内半径rti；隔热管内管外半径rto；隔热管外管内半径ri；隔热管外管外半径ro；套管内半径rci；套管外半径rco等参数。
另外，该井口注入参数包括井口注汽压力P0；井口蒸汽温度T0；注汽速率qsc；注汽时间tsc。
上述的步骤103中的根据所述井筒相关参数以及井口注入参数，以井口为初始极端点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，可以通过如下方式来实现：
计算蒸汽与隔热管内管内壁之间的热对流热阻R1：
计算隔热管内管内外壁之间的热传导热阻R2：
计算隔热层的热传导热阻R3：
计算隔热管外管壁的热传导热阻R4：
计算环空热对流热阻R5：
计算套管壁的热传导热阻R6：
计算水泥环的热传导热阻R7：
计算地层的热传导热阻R8：
其中，hf为水膜传热系数；hc为环空自然然对流传热系数；hr为环空辐射传热系数；f(t)是随时间变化的导热传热函数。
根据所述节点步长、井口蒸汽温度T0、地表温度b0、所述R1、R2、R3、R4、R6、R7、R8计算井口处热损失速度Qm：
其中，dl为所述节点步长。
根据所述P0和所述dl计算各节点下端的初始超临界蒸汽压力
其中，Pi为第i个节点上端的超临界蒸汽压力，其中，i≥0。
根据所述b0、a1、h确定地层温度Te：
Te＝(b0+a1*h/100)。
获取各节点上端的热损失速度Qi，并根据所述Qi、Te、R1、R2、R3、R4、R6、R7、R8、dl依次计算各节点下端的超临界蒸汽温度Ti+1：
Ti+1＝Te+(R1+R2+R3+R4+R6+R7+R8)*Qi/dl。
根据所述R1、R2、R3、R4、R5、R6、R7、R8计算一总热阻R：
R＝R1+R2+R3+R4+R5+R6+R7+R8
根据所述Ti+1校正各节点下端的热损失速度Qi+1：
计算各节点段的平均温度
计算各节点段的平均压力
获取井筒的超临界蒸汽相对密度γg；
计算井筒的拟临界温度Tpc：
Tpc＝92.2+176.6γg
计算井筒的拟对比温度Tpr：
计算井筒的超临界蒸汽粘度μg：
其中，
Ppc＝4.881-0.386γg
计算井筒的超临界蒸汽的体积系数Bg：
计算井筒内的超临界蒸汽流速v：
计算井筒内的超临界蒸汽ρ：
其中，为井筒内的超临界蒸汽在平均温度平均压力下的偏差系数。
当平均压力时，使用Granmer法计算偏差系数
其中，ρpr为无因次对比密度。已知计算采用迭代法，首先计算拟临界压力Ppc＝4.881-0.386rg，和拟临界温度Tpc＝92.2+176.6rg；根据计算拟对比压力与拟对比温度赋值Z(0)＝1，计算ρpr，进而计算Z(1)，控制循环次数，当Z(1)-Z(0)小于等于0.0001时结束运算。
当平均压力＞35MPa时，使用Hall-Yarborough法计算偏差系数：
T为对比温度的倒数，t＝1/Tpr
y是对比密度；
通过(A)、(B)两式联立得(两端同乘y，并相减)：
使用牛顿迭代法求解y
y(k+1)＝yk-F(yk)/F’(yk)
计算步骤已知PprTpry (0)赋初值＝0.001或
y(0)＝0.06125Ppr×t×exp[-1.2(1-t)2]带入牛顿迭代公式当F(y)≈0或
|y(k+1)-yk|≤0.00001时输出结果
其中：e为绝对粗糙度，对于已经使用的旧油管，e＝0.00001524m(0.0006in)。
之后，计算井筒内的超临界蒸汽在平均温度平均压力下的摩阻系数f：
其中，e为绝对粗糙度。
计算各节点下端超临界蒸汽压力最终值Pi+1：
其中，
此外，该计算各节点下端超临界蒸汽压力最终值Pi+1，包括：
判断的值是否小于0.1；
若的值小于0.1，则确定所述Pi+1准确；
若的值大于等于0.1，则将设置为Pi+1，并重新确定所述
若要实现上述步骤101至步骤104：
则需要满足如下条件：
(1)流体流动状态为稳定流动；
(2)流体状态为气态单相流
(3)井筒内传热为一维稳态传热；
(4)地层传热为一维不稳态传热，且满足Ramey的无因次时间函数；
(5)油套管同心。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (9)

1.一种超临界注汽井的井筒温度场及压力场分布的计算方法，其特征在于，包括：
获取井筒相关参数以及井口注入参数；所述井筒相关参数包括水泥导热系数Kcem；井眼半径rh；隔热管深度h；隔热管导热系数Kins；油管导热系数Ktub；套管导热系数Kcas；环空流体的导热系数Kha；隔热管内管内半径rti；隔热管内管外半径rto；隔热管外管内半径ri；隔热管外管外半径ro；套管内半径rci；套管外半径rco；地温梯度a1；地表温度b0；地层导热系数Ke；所述井口注入参数包括井口注汽压力P0；井口蒸汽温度T0；注汽速率qsc；注汽时间tsc；
获取预先设置的节点步长，将井筒从井口开始划分为多个节点；
根据所述井筒相关参数以及井口注入参数，以井口为初始计算点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，包括：
计算蒸汽与隔热管内管内壁之间的热对流热阻R1：
<mrow> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;h</mi> <mi>f</mi> </msub> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mfrac> </mrow>
计算隔热管内管内外壁之间的热传导热阻R2：
<mrow> <msub> <mi>R</mi> <mn>2</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;K</mi> <mrow> <mi>t</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </mrow> </mfrac> <mi>ln</mi> <mfrac> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> </msub> <msub> <mi>t</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mfrac> </mrow>
计算隔热层的热传导热阻R3：
<mrow> <msub> <mi>R</mi> <mn>3</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;K</mi> <mrow> <mi>i</mi> <mi>n</mi> <mi>s</mi> </mrow> </msub> </mrow> </mfrac> <mi>ln</mi> <mfrac> <msub> <mi>r</mi> <mi>i</mi> </msub> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> </msub> </mfrac> </mrow>
计算隔热管外管壁的热传导热阻R4：
<mrow> <msub> <mi>R</mi> <mn>4</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;K</mi> <mrow> <mi>t</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </mrow> </mfrac> <mi>ln</mi> <mfrac> <msub> <mi>r</mi> <mi>o</mi> </msub> <msub> <mi>r</mi> <mi>i</mi> </msub> </mfrac> </mrow>
计算环空热对流热阻R5：
<mrow> <msub> <mi>R</mi> <mn>5</mn> </msub> <mo>=</mo> <mfrac> <msub> <mi>r</mi> <mi>o</mi> </msub> <mrow> <mn>2</mn> <mi>&amp;pi;</mi> <mrow> <mo>(</mo> <msub> <mi>h</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>h</mi> <mi>r</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>r</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> </mrow> </mfrac> </mrow>
计算套管壁的热传导热阻R6：
<mrow> <msub> <mi>R</mi> <mn>6</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;K</mi> <mrow> <mi>c</mi> <mi>a</mi> <mi>s</mi> </mrow> </msub> </mrow> </mfrac> <mi>ln</mi> <mfrac> <msub> <mi>r</mi> <mrow> <mi>c</mi> <mi>o</mi> </mrow> </msub> <msub> <mi>r</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> </mfrac> </mrow>
计算水泥环的热传导热阻R7：
<mrow> <msub> <mi>R</mi> <mn>7</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;K</mi> <mrow> <mi>c</mi> <mi>e</mi> <mi>m</mi> </mrow> </msub> </mrow> </mfrac> <mi>ln</mi> <mfrac> <msub> <mi>r</mi> <mi>h</mi> </msub> <msub> <mi>r</mi> <mrow> <mi>c</mi> <mi>o</mi> </mrow> </msub> </mfrac> </mrow>
计算地层的热传导热阻R8：
<mrow> <msub> <mi>R</mi> <mn>8</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;K</mi> <mi>e</mi> </msub> </mrow> </mfrac> </mrow>
其中，hf为水膜传热系数；hc为环空自然然对流传热系数；hr为环空辐射传热系数；f(t)是随时间变化的导热传热函数；
根据所述节点步长、井口蒸汽温度T0、地表温度b0、所述R1、R2、R3、R4、R6、R7、R8计算井口处热损失速度Qm：
<mrow> <msub> <mi>Q</mi> <mi>m</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>T</mi> <mn>0</mn> </msub> <mo>-</mo> <msub> <mi>b</mi> <mn>0</mn> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>R</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>R</mi> <mn>3</mn> </msub> <mo>+</mo> <msub> <mi>R</mi> <mn>4</mn> </msub> <mo>+</mo> <msub> <mi>R</mi> <mn>6</mn> </msub> <mo>+</mo> <msub> <mi>R</mi> <mn>7</mn> </msub> <mo>+</mo> <msub> <mi>R</mi> <mn>8</mn> </msub> </mrow> </mfrac> <mi>d</mi> <mi>l</mi> </mrow>
其中，dl为所述节点步长；根据所述P0和所述dl计算各节点下端的初始超临界蒸汽压力
<mrow> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>c</mi> </msubsup> <mo>=</mo> <msub> <mi>P</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>0.00008</mn> <mi>d</mi> <mi>l</mi> <mo>)</mo> </mrow> </mrow>
其中，Pi为第i个节点上端的超临界蒸汽压力，其中，i≥0；
根据所述各节点下端超临界蒸汽温度和超临界蒸汽压力分别确定所述井筒的温度场分布及压力场分布。
2.根据权利要求1所述的超临界注汽井的井筒温度场及压力场分布的计算方法，其特征在于，所述根据所述井筒相关参数以及井口注入参数，以井口为初始计算点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
根据所述b0、a1、h确定地层温度Te：
Te＝(b0+a1*h/100)。
3.根据权利要求2所述的超临界注汽井的井筒温度场及压力场分布的计算方法，其特征在于，所述根据所述井筒相关参数以及井口注入参数，以井口为初始计算点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
获取各节点上端的热损失速度Qi，并根据所述Qi、Te、R1、R2、R3、R4、R6、R7、R8、dl依次计算各节点下端的超临界蒸汽温度Ti+1：
Ti+1＝Te+(R1+R2+R3+R4+R6+R7+R8)*Qi/dl。
4.根据权利要求3所述的超临界注汽井的井筒温度场及压力场分布的计算方法，其特征在于，所述根据所述井筒相关参数以及井口注入参数，以井口为初始计算点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
根据所述R1、R2、R3、R4、R5、R6、R7、R8计算一总热阻R：
R＝R1+R2+R3+R4+R5+R6+R7+R8
根据所述Ti+1校正各节点下端的热损失速度Qi+1：
<mrow> <msub> <mi>Q</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>T</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>T</mi> <mi>e</mi> </msub> </mrow> <mi>R</mi> </mfrac> <mi>d</mi> <mi>l</mi> <mo>.</mo> </mrow>
5.根据权利要求4所述的超临界注汽井的井筒温度场及压力场分布的计算方法，其特征在于，所述根据所述井筒相关参数以及井口注入参数，以井口为初始计算点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
计算各节点段的平均温度
<mrow> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> <mo>+</mo> <msub> <mi>T</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mn>2</mn> </mrow>
计算各节点段的平均压力
<mrow> <mover> <mi>P</mi> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mi>i</mi> </msub> <mo>+</mo> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>0</mn> </msubsup> <mo>)</mo> </mrow> <mo>/</mo> <mn>2.</mn> </mrow>
6.根据权利要求5所述的超临界注汽井的井筒温度场及压力场分布的计算方法，其特征在于，所述根据所述井筒相关参数以及井口注入参数，以井口为初始计算点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
获取井筒的超临界蒸汽相对密度γg；
计算井筒的拟临界温度Tpc：
Tpc＝92.2+176.6γg
计算井筒的拟对比温度Tpr：
<mrow> <msub> <mi>T</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <msub> <mi>T</mi> <mrow> <mi>p</mi> <mi>c</mi> </mrow> </msub> </mfrac> </mrow>
计算井筒的超临界蒸汽粘度μg：
<mrow> <msub> <mi>&amp;mu;</mi> <mi>g</mi> </msub> <mo>=</mo> <msub> <mi>&amp;mu;</mi> <mrow> <mi>g</mi> <mn>1</mn> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mi>A</mi> <mo>+</mo> <msub> <mi>BT</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> </msub> <mo>+</mo> <msubsup> <mi>CT</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>DT</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> <mn>3</mn> </msubsup> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>/</mo> <msub> <mi>T</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> </msub> </mrow>
其中，
<mrow> <mi>A</mi> <mo>=</mo> <mo>-</mo> <mn>2.461182</mn> <mo>+</mo> <mn>2.97054714</mn> <msub> <mi>P</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> </msub> <mo>-</mo> <mn>0.286264054</mn> <msubsup> <mi>P</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> <mn>2</mn> </msubsup> <mo>-</mo> <mn>0.00805420522</mn> <msubsup> <mi>P</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> <mn>3</mn> </msubsup> </mrow>
<mrow> <mi>B</mi> <mo>=</mo> <mn>2.80860949</mn> <mo>-</mo> <mn>3.49803305</mn> <msub> <mi>P</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> </msub> <mo>+</mo> <mn>0.36037302</mn> <msubsup> <mi>P</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> <mn>2</mn> </msubsup> <mo>-</mo> <mn>0.01044324</mn> <msubsup> <mi>P</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> <mn>3</mn> </msubsup> </mrow>
<mrow> <mi>C</mi> <mo>=</mo> <mo>-</mo> <mn>0.793385684</mn> <mo>+</mo> <mn>1.39643306</mn> <msub> <mi>P</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> </msub> <mo>-</mo> <mn>0.149144925</mn> <msubsup> <mi>P</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> <mn>2</mn> </msubsup> <mo>+</mo> <mn>0.00441015512</mn> <msubsup> <mi>P</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> <mn>3</mn> </msubsup> </mrow>
<mrow> <mi>D</mi> <mo>=</mo> <mn>0.0839387178</mn> <mo>-</mo> <mn>0.186408848</mn> <msub> <mi>P</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> </msub> <mo>+</mo> <mn>0.023367881</mn> <msubsup> <mi>P</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> <mn>2</mn> </msubsup> <mo>-</mo> <mn>0.000609579263</mn> <msubsup> <mi>P</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> <mn>3</mn> </msubsup> </mrow>
<mrow> <msub> <mi>P</mi> <mrow> <mi>p</mi> <mi>r</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mover> <mi>P</mi> <mo>&amp;OverBar;</mo> </mover> <msub> <mi>P</mi> <mrow> <mi>p</mi> <mi>c</mi> </mrow> </msub> </mfrac> </mrow>
Ppc＝4.881-0.386γg
<mrow> <msub> <mi>&amp;mu;</mi> <mrow> <mi>g</mi> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mn>3.0762</mn> <mo>&amp;times;</mo> <msup> <mn>10</mn> <mrow> <mo>-</mo> <mn>5</mn> </mrow> </msup> <mo>-</mo> <mn>3.7116</mn> <mo>&amp;times;</mo> <msup> <mn>10</mn> <mn>6</mn> </msup> <msub> <mi>&amp;gamma;</mi> <mi>g</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mo>-</mo> <mn>510.9</mn> <mo>)</mo> </mrow> <mo>+</mo> <mn>8.188</mn> <mo>&amp;times;</mo> <msup> <mn>10</mn> <mrow> <mo>-</mo> <mn>3</mn> </mrow> </msup> <mo>-</mo> <mn>6.15</mn> <mo>&amp;times;</mo> <msup> <mn>10</mn> <mrow> <mo>-</mo> <mn>3</mn> </mrow> </msup> <mi>lg</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;gamma;</mi> <mi>g</mi> </msub> <mo>)</mo> </mrow> </mrow>
计算井筒的超临界蒸汽的体积系数Bg：
<mrow> <msub> <mi>B</mi> <mi>g</mi> </msub> <mo>=</mo> <mn>3.447</mn> <mo>&amp;times;</mo> <msup> <mn>10</mn> <mrow> <mo>-</mo> <mn>4</mn> </mrow> </msup> <mfrac> <mover> <mrow> <mi>Z</mi> <mi>T</mi> </mrow> <mo>&amp;OverBar;</mo> </mover> <mover> <mi>P</mi> <mo>&amp;OverBar;</mo> </mover> </mfrac> </mrow>
计算井筒内的超临界蒸汽流速v：
<mrow> <mi>v</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>q</mi> <mrow> <mi>s</mi> <mi>c</mi> </mrow> </msub> <msub> <mi>B</mi> <mi>g</mi> </msub> </mrow> <mrow> <mn>86400</mn> <msup> <msub> <mi>&amp;pi;r</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> <mn>2</mn> </msup> <mo>/</mo> <mn>4</mn> </mrow> </mfrac> </mrow>
计算井筒内的超临界蒸汽ρ：
<mrow> <mi>&amp;rho;</mi> <mo>=</mo> <mfrac> <mrow> <mn>3484.4</mn> <msub> <mi>&amp;gamma;</mi> <mi>g</mi> </msub> <mover> <mi>P</mi> <mo>&amp;OverBar;</mo> </mover> </mrow> <mrow> <mover> <mi>Z</mi> <mo>&amp;OverBar;</mo> </mover> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> </mrow> </mfrac> </mrow>
其中，为井筒内的超临界蒸汽在平均温度平均压力下的偏差系数。
7.根据权利要求6所述的超临界注汽井的井筒温度场及压力场分布的计算方法，其特征在于，所述根据所述井筒相关参数以及井口注入参数，以井口为初始计算点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
计算井筒内的超临界蒸汽在平均温度平均压力下的摩阻系数f：
<mrow> <mfrac> <mn>1</mn> <msqrt> <mi>f</mi> </msqrt> </mfrac> <mo>=</mo> <mn>1.14</mn> <mo>-</mo> <mn>2</mn> <mi>l</mi> <mi>g</mi> <mrow> <mo>(</mo> <mfrac> <mi>e</mi> <mi>d</mi> </mfrac> <mo>+</mo> <mfrac> <mn>21.25</mn> <msup> <mi>Re</mi> <mn>0.9</mn> </msup> </mfrac> <mo>)</mo> </mrow> </mrow>
其中，e为绝对粗糙度，d为管径，Re为雷诺数。
8.根据权利要求7所述的超临界注汽井的井筒温度场及压力场分布的计算方法，其特征在于，所述根据所述井筒相关参数以及井口注入参数，以井口为初始计算点，依次计算各节点下端超临界蒸汽温度和超临界蒸汽压力，还包括：
计算各节点下端超临界蒸汽压力最终值Pi+1：
<mrow> <msub> <mi>p</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msqrt> <mrow> <msubsup> <mi>p</mi> <mi>i</mi> <mn>2</mn> </msubsup> <msup> <mi>e</mi> <mrow> <mn>2</mn> <mi>S</mi> </mrow> </msup> <mo>-</mo> <mi>&amp;beta;</mi> <mrow> <mo>(</mo> <msup> <mi>e</mi> <mrow> <mn>2</mn> <mi>S</mi> </mrow> </msup> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msqrt> </mrow>
其中，
9.根据权利要求7所述的超临界注汽井的井筒温度场及压力场分布的计算方法，其特征在于，所述计算各节点下端超临界蒸汽压力最终值Pi+1，包括：
判断的值是否小于0.1；
若的值小于0.1，则确定所述Pi+1准确；
若的值大于等于0.1，则将设置为Pi+1，并重新确定所述
CN201510084615.1A 2015-02-16 2015-02-16 超临界注汽井的井筒温度场及压力场分布的计算方法 Active CN104806230B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510084615.1A CN104806230B (zh) 2015-02-16 2015-02-16 超临界注汽井的井筒温度场及压力场分布的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510084615.1A CN104806230B (zh) 2015-02-16 2015-02-16 超临界注汽井的井筒温度场及压力场分布的计算方法

Publications (2)

Publication Number Publication Date
CN104806230A CN104806230A (zh) 2015-07-29
CN104806230B true CN104806230B (zh) 2018-01-05

Family

ID=53691382

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510084615.1A Active CN104806230B (zh) 2015-02-16 2015-02-16 超临界注汽井的井筒温度场及压力场分布的计算方法

Country Status (1)

Country Link
CN (1) CN104806230B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110929447A (zh) * 2019-12-20 2020-03-27 西南石油大学 一种稠化酸酸化过程中井筒温度场数值计算方法

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105178930B (zh) * 2015-08-17 2018-02-02 中国石油天然气股份有限公司 地面注汽管线内蒸汽热力参数计算方法及装置
CN105134179B (zh) * 2015-08-21 2018-02-02 中国石油天然气股份有限公司 一种天然气井井筒压力和温度的分布数据的计算方法
CN105114060B (zh) * 2015-08-26 2018-02-02 中国石油天然气股份有限公司 一种计算注气井的井筒温度分布的方法及装置
CN105403590B (zh) * 2015-11-12 2018-06-01 中国石油天然气股份有限公司 隔热管导热系数测试方法及装置
CN105422084B (zh) * 2015-11-26 2018-07-13 中国石油天然气股份有限公司 一种热水循环加热降粘工艺的井筒温度场获取方法及装置
CN105370266B (zh) * 2015-12-01 2018-06-01 中国石油天然气股份有限公司 火烧油层分层电点火注气井井筒温度分布确定方法及装置
CN105370255B (zh) * 2015-12-01 2018-01-05 中国石油天然气股份有限公司 火烧油层笼统注气电点火井筒温度分布的确定方法及装置
CN106021958B (zh) * 2016-06-14 2019-04-12 中国石油天然气股份有限公司 一种注气井封隔器上下两端温度的确定方法及装置
CN106018472B (zh) * 2016-07-08 2019-03-12 中国石油天然气股份有限公司 注汽井隔热管导热系数测试方法及装置
CN106285600B (zh) * 2016-08-11 2018-08-14 中国石油天然气集团公司 获取油气压裂增产中井底温度场的方法及装置
CN106503299A (zh) * 2016-09-27 2017-03-15 中国石油天然气股份有限公司 注蒸汽井井筒总传热系数计算方法
CN106951599A (zh) * 2017-02-27 2017-07-14 南京航空航天大学 一种飞机凝结尾生成对全球地表温度影响的评估方法
CN110344796B (zh) * 2018-04-04 2021-08-06 中国石油化工股份有限公司 基于三参数原理的注汽参数预测方法
CN109826602B (zh) * 2019-01-24 2019-10-22 西南石油大学 注汽压力高的主控因素判断及治理手段优选方法
CN112796704B (zh) * 2019-10-28 2023-07-28 中国石油化工股份有限公司 油井热洗清蜡方式优选与参数优化方法
CN113051790B (zh) * 2021-03-04 2023-11-24 西安建筑科技大学 用于有限元模拟的蒸汽荷载加载方法、系统、设备及介质
CN113821998B (zh) * 2021-08-31 2024-03-29 中国船舶重工集团公司第七0三研究所 利用牛顿迭代法进行凝汽器实时动态仿真模型壳侧压力求解的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003078791A2 (en) * 2002-03-13 2003-09-25 Weatherford/Lamb, Inc. Method and apparatus for injecting steam into a geological formation
CN103590818A (zh) * 2013-10-21 2014-02-19 中国石油天然气股份有限公司 一种井筒内温度分布半解析确定方法及装置
CN103899295A (zh) * 2014-03-10 2014-07-02 中国石油天然气股份有限公司 稠油蒸汽采注中直井井筒热损失检测方法及检测装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003078791A2 (en) * 2002-03-13 2003-09-25 Weatherford/Lamb, Inc. Method and apparatus for injecting steam into a geological formation
CN103590818A (zh) * 2013-10-21 2014-02-19 中国石油天然气股份有限公司 一种井筒内温度分布半解析确定方法及装置
CN103899295A (zh) * 2014-03-10 2014-07-02 中国石油天然气股份有限公司 稠油蒸汽采注中直井井筒热损失检测方法及检测装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
注蒸汽水平井井筒内参数计算新模型;王一平等;《西南石油大学学报（自然科学版）》;20100815;第32卷(第4期);第127-132页 *
注超临界气体井筒温度压力场计算方法;刘洪波等;《石油大学学报（自然科学版）》;20040220;第28卷(第1期);第52-54页 *
注过热蒸气井筒物性参数计算综合数学模型;许安著等;《大庆石油学院学报》;20090229;第33卷(第1期);第29-35页 *
蒸汽驱注蒸汽井筒热损失分析;张军涛等;《中国石油和化工石油工程技术》;20120105;第55-57页 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110929447A (zh) * 2019-12-20 2020-03-27 西南石油大学 一种稠化酸酸化过程中井筒温度场数值计算方法
CN110929447B (zh) * 2019-12-20 2020-11-03 西南石油大学 一种稠化酸酸化过程中井筒温度场数值计算方法

Also Published As

Publication number Publication date
CN104806230A (zh) 2015-07-29

Similar Documents

Publication Publication Date Title
CN104806230B (zh) 超临界注汽井的井筒温度场及压力场分布的计算方法
Alimonti et al. The wellbore heat exchangers: A technical review
Song et al. Numerical analysis of the heat production performance of a closed loop geothermal system
Gupta et al. Developing empirical heat-transfer correlations for supercritical CO2 flowing in vertical bare tubes
Alimonti et al. Study of geothermal power generation from a very deep oil well with a wellbore heat exchanger
Vijayan et al. Steady state and stability characteristics of single-phase natural circulation in a rectangular loop with different heater and cooler orientations
CN104453861B (zh) 一种高压气井井筒温度分布的确定方法以及系统
CN103775058B (zh) 一种井筒热损失的确定方法
CN104481482B (zh) 水平井同心双管注气隔热分析方法及装置
CN109800527B (zh) 针对中深层地埋套管换热器内自循环换热量的计算方法
Han et al. Fast thermal simulation of a heated crude oil pipeline with a BFC-Based POD reduced-order model
Luo et al. Proposing stratified segmented finite line source (SS-FLS) method for dynamic simulation of medium-deep coaxial borehole heat exchanger in multiple ground layers
Nian et al. Study on the effect of wellbore heat capacity on steam injection well heat loss
CN107575214A (zh) 用于注采过程的井筒内温度与压力的预测方法
Sun et al. Geothermal exploitation by circulating supercritical CO2 in a closed horizontal wellbore
CN105134179B (zh) 一种天然气井井筒压力和温度的分布数据的计算方法
CN107145705A (zh) 一种获取循环温度场的方法及装置
Sivaramkrishnan et al. Predicting wellbore dynamics in a steam-assisted gravity drainage system: numeric and semi-analytic model, and validation
Jiashu et al. The influence of thermal boundary conditions of wellbore on the heat extraction performance of deep borehole heat exchangers
CN102777159B (zh) 一种注co2气井井筒流态确定及参数优化方法
Yao et al. Characteristics of phase-change flow and heat transfer in loop thermosyphon: Three-dimension CFD modeling and experimentation
CN109033012B (zh) 一种空心抽油杆注热水循环井筒温度场的确定方法及装置
Ma et al. Analysis and applications of a two-phase closed thermosyphon for improving the fluid temperature distribution in wellbores
Samir et al. Simulation of natural convection boiling heat transfer for refrigerant R-134a flow in a metal foam filled vertical tube
CN105403590B (zh) 隔热管导热系数测试方法及装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant