CN104481482B - 水平井同心双管注气隔热分析方法及装置 - Google Patents

Abstract

本发明公开一种水平井同心双管注气隔热分析方法及装置，包括：采集井身结构及热物理性质参数、井口注入参数以及地层热物性参数，确定井身结构热阻，进一步确定井口处的最大热损失、单位长度的地层平均温度、隔热管向井壁的单位长度的热损失速度、内管的单位长度热损失速度、隔热管与内管环空的蒸汽干度以及内管的蒸汽干度，以确定环空对流热阻，根据环空对流热阻对单位长度的隔热管向井壁的热损失速度进行校正，根据校正后的隔热管向井壁的单位长度的热损失速度进行迭代分析生成水平井同心双管注气隔热分析结果。本方案具有良好的计算稳定性和较高的计算精度，可以对水平井同心双管注汽工艺参数进行预测和调整，对提高油田的稠油水平井产量具有实意义。

Description

水平井同心双管注气隔热分析方法及装置

技术领域

本发明专利涉及稠油注蒸汽热采领域，具体的讲是一种水平井同心双管注气隔热分析方法及装置。

背景技术

从世界范围看，水平井的开发应用已经成为油藏开发的主流，现场试验证明，稠油油藏采取水平井同心双管注汽方式，能够提高水平段吸汽的均匀程度，从而提高蒸汽的波及面积，是进一步提高稠油油藏采收率的有效开采方式。

由于以往同心双管注汽隔热分析理论计算研究工作不完整，同心双管注汽工艺参数优化设计研究基本仍属空白，目前急需配套集方案设计、系统参数调配、注采参数优化的软件平台，以指导现场进行施工方案设计和最佳开采工艺实施。

发明内容

为建立集工艺方案设计、系统参数调配、注采参数优化的水平井双管注汽工艺参数优化设计模型，以指导现场进行施工方案设计和最佳开采工艺实施，本发明实施例提供了一种水平井同心双管注气隔热分析方法，包括：

步骤1，采集井身结构及热物理性质参数、井口注入参数以及地层热物性参数；

步骤2，根据所述井身结构及热物理性质参数、井口注入参数以及地层热物性参数确定井身结构热阻；

步骤3，根据井口注入参数、所述井身结构热阻以及地层热物性参数确定井口处的最大热损失、单位长度的地层平均温度、隔热管向井壁的单位长度的热损失速度、内管的单位长度热损失速度、隔热管与内管环空的蒸汽干度以及内管的蒸汽干度；

步骤4，根据贝格斯-比尔算法确定单位长度的内管与隔热管的环空平均压力、单位长度的内管平均压力；

步骤5，根据所述的单位长度的内管与隔热管的环空平均压力、单位长度的内管平均压力确定单位长度的内管和隔热管环空的蒸汽平均温度、内管的蒸汽平均温度；

步骤6、根据所述单位长度的内管和隔热管环空的蒸汽平均温度、单位长度的地层平均温度以及内管的蒸汽平均温度确定井身结构温度；

步骤7，根据确定的井身结构温度确定环空辐射传热系数和环空对流传热系数；

步骤8，根据确定的环空辐射传热系数和环空对流传热系数确定环空对流热阻；

步骤9，根据确定的环空对流热阻对所述单位长度的隔热管向井壁的热损失速度进行校正生成校正后的隔热管向井壁的单位长度的热损失速度；

步骤10，根据所述校正后的隔热管向井壁的单位长度的热损失速度进行迭代分析并重复步骤3～步骤10生成水平井同心双管注气隔热分析结果。

此外，本发明还公开一种水平井同心双管注气隔热分析装置，包括：

参数采集模块，用于采集井身结构及热物理性质参数、井口注入参数以及地层热物性参数；

井身结构热阻确定模块，用于根据所述井身结构及热物理性质参数、井口注入参数以及地层热物性参数确定井身结构热阻；

温度干度确定模块，用于根据井口注入参数、所述井身结构热阻以及地层热物性参数确定井口处的最大热损失、单位长度的地层平均温度、隔热管向井壁的单位长度的热损失速度、内管的单位长度热损失速度、隔热管与内管环空的蒸汽干度以及内管的蒸汽干度；

压力确定模块，根据贝格斯-比尔算法确定单位长度的内管与隔热管的环空平均压力、单位长度的内管平均压力；

平均温度确定模块，用于根据所述的单位长度的内管与隔热管的环空平均压力、单位长度的内管平均压力确定单位长度的内管和隔热管环空的蒸汽平均温度、内管的蒸汽平均温度；

井身结构温度确定模块，用于根据所述单位长度的内管和隔热管环空的蒸汽平均温度、单位长度的地层平均温度以及内管的蒸汽平均温度确定井身结构温度；

传热系数确定模块，用于根据确定的井身结构温度确定环空辐射传热系数和环空对流传热系数；

环空对流热阻确定模块，用于根据确定的环空辐射传热系数和环空对流传热系数确定环空对流热阻；

校正热损失确定模块，根据确定的环空对流热阻对所述单位长度的隔热管向井壁的热损失速度进行校正生成校正后的隔热管向井壁的单位长度的热损失速度；

隔热分结果生成模块，根据所述校正后的隔热管向井壁的单位长度的热损失速度进行迭代分析生成水平井同心双管注气隔热分析结果。

本发明公开的技术方案具有良好的计算稳定性和较高的计算精度，通过本发明的方案可以很好的对水平井同心双管注汽工艺参数进行预测和调整。本发明对提高油田的稠油水平井产量具有极大的现实意义。

为让本发明的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明公开的一种水平井同心双管注气隔热分析方法的流程图；

图2为本发明实施例的示意图；

图3为本发明公开的一种水平井同心双管注气隔热分析装置的框图；

图4为本发明实施例的示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

如图1所示，为本发明公开的一种水平井同心双管注气隔热分析方法，包括：

步骤S101，采集井身结构及热物理性质参数、井口注入参数以及地层热物性参数；

步骤S102，根据所述井身结构及热物理性质参数、井口注入参数以及地层热物性参数确定井身结构热阻；

步骤S103，根据井口注入参数、所述井身结构热阻以及地层热物性参数确定井口处的最大热损失、单位长度的地层平均温度、隔热管向井壁的单位长度的热损失速度、内管的单位长度热损失速度、隔热管与内管环空的蒸汽干度以及内管的蒸汽干度；

步骤S104，根据贝格斯-比尔算法确定单位长度的内管与隔热管的环空平均压力、单位长度的内管平均压力；

步骤S105，根据所述的单位长度的内管与隔热管的环空平均压力、单位长度的内管平均压力确定单位长度的内管和隔热管环空的蒸汽平均温度、内管的蒸汽平均温度；

步骤S106、根据所述单位长度的内管和隔热管环空的蒸汽平均温度、单位长度的地层平均温度以及内管的蒸汽平均温度确定井身结构温度；

步骤S107，根据确定的井身结构温度确定环空辐射传热系数和环空对流传热系数；

步骤S108，根据确定的环空辐射传热系数和环空对流传热系数确定环空对流热阻；

步骤S109，根据确定的环空对流热阻对所述单位长度的隔热管向井壁的热损失速度进行校正生成校正后的隔热管向井壁的单位长度的热损失速度；

步骤S110，根据所述校正后的隔热管向井壁的单位长度的热损失速度进行迭代分析并重复步骤S103～步骤S110生成水平井同心双管注气隔热分析结果。

水平井同心双管注汽，如图2所示，注汽内管201一般采用油管，注汽外管202般采用隔热管。井口两套注汽调节装置，分别控制注汽内管和注汽外管的注汽压力、注汽速率、注汽时间。内管长度大于外管长度，外管下深到脚跟处，内管下深到脚尖处。水平段是筛管方式完井。

由于隔热管保温效果好，内管又相当于一个加热体，因此外管段热量损失较小，外管出口处有较高的干度和温度。由于普通油管导热系数一般都在150～160kJ/(m·h·℃)，热损失较大，但是外管对内管起到了很好的保温作用，因此内管在直井段热损失较小。内管在水平段热损失较大，但内管的散热对水平段油层起到了很好的预热效果。

此外，如图3所示，本发明还提供了一种水平井同心双管注气隔热分析装置，包括：

参数采集模块301，用于采集井身结构及热物理性质参数、井口注入参数以及地层热物性参数；

井身结构热阻确定模块302，用于根据所述井身结构及热物理性质参数、井口注入参数以及地层热物性参数确定井身结构热阻；

温度干度确定模块303，用于根据井口注入参数、所述井身结构热阻以及地层热物性参数确定井口处的最大热损失、单位长度的地层平均温度、隔热管向井壁的单位长度的热损失速度、内管的单位长度热损失速度、隔热管与内管环空的蒸汽干度以及内管的蒸汽干度；

压力确定模块304，根据贝格斯-比尔算法确定单位长度的内管与隔热管的环空平均压力、单位长度的内管平均压力；

平均温度确定模块305，用于根据所述的单位长度的内管与隔热管的环空平均压力、单位长度的内管平均压力确定单位长度的内管和隔热管环空的蒸汽平均温度、内管的蒸汽平均温度；

井身结构温度确定模块306，用于根据所述单位长度的内管和隔热管环空的蒸汽平均温度、单位长度的地层平均温度以及内管的蒸汽平均温度确定井身结构温度；

传热系数确定模块307，用于根据确定的井身结构温度确定环空辐射传热系数和环空对流传热系数；

环空对流热阻确定模块308，用于根据确定的环空辐射传热系数和环空对流传热系数确定环空对流热阻；

校正热损失确定模块309，根据确定的环空对流热阻对所述单位长度的隔热管向井壁的热损失速度进行校正生成校正后的隔热管向井壁的单位长度的热损失速度；

隔热分结果生成模块310，根据所述校正后的隔热管向井壁的单位长度的热损失速度进行迭代分析生成水平井同心双管注气隔热分析结果。

下面结合具体实施例对本发明的技术方案作进一步详细描述，本实施例中仅列出为推导新算法用到的主要计算公式。注汽管柱划分成若干个单元，每段单元长度dl，一段单元从直井段注汽内管中心线截开，侧面图如图4所示。

蒸汽与内管内壁之间的热对流：

式中：Q为井筒单元径向热损失，KJ/h；T_{s_i}为内管注入蒸汽温度，℃；T_{ti_i}为内管内壁温度，℃；R₁为蒸汽与内管内壁之间的热对流热阻，m·K/W；h_f为水膜传热系数，W/(m²·K)；r_{ti_i}为内管内壁半径，m。

内管内外壁之间的热传导：

式中：T_{to_i}为内管外壁温度，℃；R₂为内管内外壁之间的热传导热阻，m·K/W；K_tub为油管导热系数，W/(m·K)；r_{to_i}为内管外壁半径，m。

蒸汽与内管外壁之间的热对流：

式中：T_s为隔热管注入蒸汽温度，℃；R₃为蒸汽与内管外壁之间的热对流热阻，m·K/W。

蒸汽与隔热管内管内壁之间的热对流：

式中：T_s为注入蒸汽温度，℃；T_ti为隔热管内管内壁温度，℃；R₄为蒸汽与隔热管内管内壁之间的热对流热阻，m·K/W；r_ti为隔热管内管内壁半径，m。

隔热管内管内外壁之间的热传导：

式中：T_to为隔热管内管外壁温度，℃；R₅为隔热管内管内外壁之间的热传导热阻，m·K/W；r_to为隔热管内管外壁半径，m。

隔热层的热传导：

式中：T_i为隔热管外管内壁温度，℃；R₆为绝热层的热传导热阻，m·K/W；K_ins为隔热管导热系数，W/(m·K)；r_i为隔热管外管内壁半径，m。

隔热管外壁的热传导：

式中：T_o为隔热管外壁温度，℃；R₇为隔热管外管内外壁之间的热传导热阻，m·K/W；r_o为隔热管外管外壁半径，m。

环空热对流：

式中：T_ci为套管内壁温度，℃；R₈为竖井段的环空热对流热阻，m·K/W；hr为环空辐射传热系数，W/(m²·K)；hc为环空自然然对流传热系数，W/(m²·K)。

套管壁的热传导：

式中：T_co为套管外壁温度，℃；R₉为套管壁的热传导热阻，m·K/W；K_cas为套管导热系数，W/(m·K)；r_ci为套管内壁半径，m；r_co为套管外壁半径，m。

水泥环的热传导：

式中：T_h为水泥环与地层交界处温度，℃；R₁₀为水泥环的热传导热阻，m·K/W；K_cem为水泥环导热系数，W/(m·K)；r_h为井眼半径，m。

地层的热传导：

式中：T_e为地层平均温度，℃；R₁₁为地层的热传导热阻，m·K/W；K_e为地层导热系数，W/(m·K)。

地层温度b_k的计算：

b_k＝(b_k-1+a1*dl/100) (12)

式中：地层温度初始值b₀为地表温度，℃；a1为地温梯度，℃/100m。

压力计算公式使用Beggs-Bill算法：

式中

p──混合物的压力(绝对)，Pa；

z──轴向流动的距离，m；

ρ_l──液相密度，kg/m³；

ρ_g──气相密度，kg/m³；

H_l──持液率，m³/m³；

g──重力加速度，m/s²；

θ──管道与水平方向的夹角，°；

λ──两相流动的沿程阻力系数，无因次；

G──混合物的质量流量，kg/s；

v──混合物的流速，m/s；

v_sg──气相的折算速度，m/s；

D──管道直径，m；

A──管道截面积，m²。(具体计算详见张琪主编《采油工程原理与设计》P50)干度计算公式：

其中

c₁＝G(h_s-h_w)

式中

G──注入蒸汽的质量流量，kg/s；

dq──井段长度dz上的热损失，W；

ρ_m──混合物的密度，kg/m³；

h_w──饱和水的焓，J/kg；(由水蒸汽的热力学性质，用插值法求得)

h_s──干饱和蒸汽的焓，J/kg；(由水蒸汽的热力学性质，用插值法求得)

x──蒸汽干度。

经整理，可简化为

其中

c₁＝G(h_s-h_w)

下面是新的算法：

由式(4)～(7)得：

T_o＝T_s-(R₄+R₅+R₆+R₇)Q/dl (15)

由式(9)～(11)得：

T_ci＝T_e+(R₉+R₁₀+R₁₁)Q/dl (16)

将(14)、(15)代入(8)得隔热管向井壁的单位长度热损失：

该式给出了隔热管向井壁的单位长度热损失速度的上限值。

1、直井段注汽隔热分析计算的迭代算法：

(1)计算热阻R₁，R₂，R₃，R₄，R₅，R₆，R₇，R₈，R₉，R₁₀，R₁₁

(2)从井口开始计算，井口处热损失最大，l＝0,k＝1,令Q_k＝0.9Q_m。

(3)通过公式(12)计算地层温度平均温度T_e＝(b_k+b_k-1)/2

(4)通过公式(13)Beggs-Bill算法计算该段内管和隔热管环空的平均压力P和内管的平均压力P′

(5)通过公式T_s＝195.94P^0.225-17.8计算该段内管和隔热管环空的蒸汽平均温度T_s和内管的蒸汽平均温度T_{s_i}

(6)计算以下温度：

隔热管内管内壁温度为：T_ti＝T_s-R₄Q_k/dl (18)

隔热管内管外壁温度为：T_to＝T_ti-R₅Q_k/dl (19)

隔热管外管内壁温度为：T_i＝T_to-R₆Q_k/dl (20)

隔热管外管外壁温度为：T_o＝T_i-R₇Q_k/dl (21)

水泥环外温度为：T_h＝T_e+R₁₁Q_k/dl (22)

套管外壁温度为：T_co＝T_h+R₁₀Q_k/dl (23)

套管内壁温度为：T_ci＝T_co+R₉Q_k/dl (24)

(7)计算环空辐射传热系数：

式中：δ为Stefan-Boltzmann常数，2.189×10^-8W/(m²·K)；F_tci为油管或绝热管外壁表面向套管内壁表面辐射有效系数；ε_o为绝热管外壁黑度；ε_ci为套管内壁黑度。

(8)计算环空自然对流传热系数：

式中：Gr为Grashof数；Pr为Prandtl数；K_ha为环空流体的导热系数，W/(m·K)；g为重力加速度，m/s²；ρ_an为环空流体在平均温度T_an下的密度，kg/m³；ρ_an为环空流体在平均温度T_an下的粘度，mPa·s；C_an为环空流体在平均温度T_an下的热容，J(m³·K)。

(9)计算隔热管向井壁的单位长度总热阻R：

R＝R₄+R₅+R₆+R₇+R₈+R₉+R₁₀+R₁₁ (31)

(10)隔热管向井壁的单位长度热损失速度校正：

若Q_k≥Q_m,则令Q_k＝Q_m。

(11)由式(1)～(3)得,内管的竖井段的单位长度热损失速度为：

(12)通过公式(14)计算该段隔热管与内管环空的蒸汽干度x和内管的蒸汽干度x′

(13)k＝k+1,令Q_k-1＝Q_k,b_k-1＝b_k,P_k-1＝P_k,x_k-1＝x_k，x′_k-1＝x′_k，l＝l+dl,返回第(3)步继续迭代计算。若l>＝h,则迭代结束。(h为隔热管下入深度)。

2、一般内管长度大于外管长度，外管出口在脚跟处，水平井段注汽隔热分析接着上面直井段最后计算的结果继续进行迭代计算：

(1)通过公式(12)计算地层温度T_e＝b₀+a1*h/100

(2)通过公式(13)Beggs-Bill算法计算该段内管的平均压力P′

(3)通过公式(34)计算该段内管的蒸汽平均温度T_{s_i}

x₀＝-0.02447P+0.81666 (35)

T_s＝41.41(P×1014)^0.225-17.78 (36)

式中：x为按照静态计算得到蒸汽干度，T_s为按照饱和蒸汽温度－压力关系计算得到的温度，T_w为按照等质量单相液态水计算得到温度。

(4)计算以下温度：

内管内壁温度为：T_{ti_i}＝T_{s_i}-R₁Q_k/dl (37)

内管外壁温度为：T_{to_i}＝T_{ti_i}-R₂Q_k/dl (38)

筛管外管温度为：T_co＝T_e+R₁₁Q_k/dl (39)

筛管内管温度为：T_ci＝T_co+R′₉Q_k/dl (40)

(5)计算环空辐射传热系数：

式中：F_tci为筛管外壁表面向筛管内壁表面辐射有效系数；ε_o为绝热管外壁黑度；ε_ci为筛管内壁黑度。

(6)计算环空自然对流传热系数：

(7)通过公式(46)计算该段内管的热损失速度为：

环空热对流热阻：

筛管壁的热传导热阻：

式中：K′_cas为筛管导热系数，W/(m·K)；r′_ci为筛管内壁半径，m；r′_co为筛管外壁半径，m。

(8)通过公式(14)计算该段内管的蒸汽干度x′

(9)k＝k+1,令Q′_k-1＝Q′_k,P′_k-1＝P′_k,x′_k-1＝x′_k，l＝l+dl,返回第(2)步继续迭代计算。若l>＝H,则迭代结束。(H为内管下入深度)。

综上，本发明具有良好的计算稳定性和较高的计算精度，通过本发明的算法和解释工具可以很好的对水平井同心双管注汽工艺参数进行预测和调整。对完善国内水平井双管注汽工艺参数优化设计研究领域，对水平井同心双管注汽的推广和应用提供了坚实的理论基础，对提高油田的稠油水平井产量具有极大的现实意义。

本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (24)

1.一种水平井同心双管注气隔热分析方法，其特征在于，所述的方法包括：

步骤1，采集井身结构及热物理性质参数、井口注入参数以及地层热物性参数；

步骤2，根据所述井身结构及热物理性质参数、井口注入参数以及地层热物性参数确定井身结构热阻；

步骤3，根据井口注入参数、所述井身结构热阻以及地层热物性参数确定井口处的最大热损失、单位长度的地层平均温度、隔热管向井壁的单位长度的热损失速度、内管的单位长度热损失速度、隔热管与内管环空的蒸汽干度以及内管的蒸汽干度；

步骤4，根据贝格斯-比尔算法确定单位长度的内管与隔热管的环空平均压力、单位长度的内管平均压力；

步骤5，根据所述的单位长度的内管与隔热管的环空平均压力、单位长度的内管平均压力确定单位长度的内管和隔热管环空的蒸汽平均温度、内管的蒸汽平均温度；

步骤6、根据所述单位长度的内管和隔热管环空的蒸汽平均温度、单位长度的地层平均温度以及内管的蒸汽平均温度确定井身结构温度；

步骤7，根据确定的井身结构温度确定环空辐射传热系数和环空对流传热系数；

步骤8，根据确定的环空辐射传热系数和环空对流传热系数确定环空对流热阻；

步骤9，根据确定的环空对流热阻对所述单位长度的隔热管向井壁的热损失速度进行校正生成校正后的隔热管向井壁的单位长度的热损失速度；

步骤10，根据所述校正后的隔热管向井壁的单位长度的热损失速度进行迭代分析并重复步骤3～步骤10生成水平井同心双管注气隔热分析结果。

2.如权利要求1所述的水平井同心双管注气隔热分析方法，其特征在于，所述的井身结构热阻包括：蒸汽与内管内壁之间的热对流热阻、内管内外壁热传导热阻、蒸汽与内管外壁之间的热对流热阻、蒸汽与隔热管内管内壁之间的热对流热阻、隔热管的内管内外壁之间的热传导热阻、隔热层的热传导热阻、隔热管外管内外壁之间的热传导热阻、套管壁的热传导热阻、水泥环的热传导热阻以及地层的热传导热阻。

3.如权利要求1所述的水平井同心双管注气隔热分析方法，其特征在于，所述的步骤10中生成的水平井同心双管注气隔热分析结果包括：

生成直井段注汽隔热分析结果和水平井段注汽隔热分析结果。

4.如权利要求3所述的水平井同心双管注气隔热分析方法，其特征在于，所述的井身结构温度包括：隔热管内管内壁温度、隔热管内管外壁温度、隔热管外管内壁温度、隔热管外管外壁温度、水泥环外温度、套管外壁温度、套管内壁温度、内管内壁温度、筛管外管温度、筛管内管温度。

5.如权利要求1所述的水平井同心双管注气隔热分析方法，其特征在于，

所述的井身结构及热物理性质参数包括：水膜传热系数、水泥导热系数、井眼半径、内管内壁半径、内管外壁半径、内管下入深度、内管导热系数、隔热管下入深度、隔热管导热系数、隔热管内管内径、隔热管内管外径、隔热管外管外径、隔热管外管内径、套管内径、套管外径、筛管内径、筛管外径、油管导热系数、套管导热系数、筛管导热系数；

所述的井口注入参数包括：内管注汽压力；内管注汽干度；内管注汽速率；内管注汽时间；外管注汽压力；外管注汽干度；外管注汽速率；外管注汽时间；

所述的地层热物性参数包括：地温梯度；地表温度；地层导热系数。

6.如权利要求5所述的水平井同心双管注气隔热分析方法，其特征在于，所述的根据所述井身结构及热物理性质参数、井口注入参数以及地层热物性参数确定井身结构热阻包括：

根据所述的水膜传热系数、内管内壁半径以及式(1)确定蒸汽与内管内壁之间的热对流热阻；

$R_1=\frac{1}{2{\mathrm{\πh}}_f{r}_{ti\_i}}---\left(1\right)$

其中，R₁为蒸汽与内管内壁之间的热对流热阻，m·K/W；h_f为水膜传热系数，2800～11136W/(m²·K)；r_{ti_i}为内管内壁半径，m；

根据所述油管导热系数、内管内壁半径、内管外壁半径及式(2)确定内管内外壁之间的热传导热阻；

$R_2=\frac{1}{2{\mathrm{\πK}}_{tub}}\ln \frac{r_{to\_i}}{r_{ti\_i}}---\left(2\right)$

式中：R₂为内管内外壁之间的热传导热阻，m·K/W；K_tub为油管导热系数，W/(m·K)；r_{to_i}为内管外壁半径，m；

所述的水膜传热系数、内管内壁半径以及式(3)确定蒸汽与内管外壁之间的热对流热阻；

$R_3=\frac{1}{2{\mathrm{\πh}}_f{r}_{to\_i}}---\left(3\right)$

式中：R₃为蒸汽与内管外壁之间的热对流热阻，m·K/W；

根据水膜传热系数、隔热管内管内壁半径以及式(4)确定蒸汽与隔热管内管内壁之间的热对流热阻；

$R_4=\frac{1}{2{\mathrm{\πh}}_f{r}_{ti}}---\left(4\right)$

式中：R₄为蒸汽与隔热管内管内壁之间的热对流热阻，m·K/W；r_ti为隔热管内管内壁半径，m；

根据油管导热系数、隔热管内管内壁半径、隔热管内管外壁半径和式(5)确定隔热管内管内外壁之间的热传导热阻：

$R_5=\frac{1}{2{\mathrm{\πK}}_{tub}}ln\frac{r_{to}}{r_{ti}}---\left(5\right)$

式中：R₅为隔热管内管内外壁之间的热传导热阻，m·K/W，K_tub为油管导热系数，W/(m·K)；r_to为隔热管内管外壁半径，m；

根据所述的隔热管导热系数、隔热管内管外壁半径、隔热管外管内壁半径以及式(6)确定隔热层的热传导热阻；

$R_6=\frac{1}{2{\mathrm{\πK}}_{ins}}ln\frac{r_i}{r_{to}}---\left(6\right)$

式中：R₆为绝热层的热传导热阻，m·K/W；K_ins为隔热管导热系数，W/(m·K)；r_i为隔热管外管内壁半径，m；

隔热管外管内外壁之间的热传导热阻：

$R_7=\frac{1}{2{\mathrm{\πK}}_{tub}}ln\frac{r_o}{r_i}---\left(7\right)$

式中：R₇为隔热管外管内外壁之间的热传导热阻，m·K/W，r_o为隔热管外管外壁半径，m，

根据所述套管导热系数、套管内壁半径、套管外壁半径以及式(9)确定套管壁的热传导热阻；

$R_9=\frac{1}{2{\mathrm{\πK}}_{cas}}ln\frac{r_{co}}{r_{ci}}---\left(9\right)$

式中：R₉为套管壁的热传导热阻，m·K/W，K_cas为套管导热系数，W/(m·K)，r_ci为套管内壁半径，m，r_co为套管外壁半径，m；

根据水泥环导热系数、井眼半径、套管外壁半径以及式(10)确定水泥环的热传导热阻：

$R_{10}=\frac{1}{2{\mathrm{\πK}}_{cem}}ln\frac{r_h}{r_{co}}---\left(10\right)$

式中：R₁₀为水泥环的热传导热阻，m·K/W；K_cem为水泥环导热系数，W/(m·K)；r_h为井眼半径，m；

根据瑞曼时间函数、地层导入系数以及式(11)确定地层的热传导热阻：

$R_{11}=\frac{f\left(t\right)}{2{\mathrm{\πK}}_e}---\left(11\right)$

式中：R₁₁为地层的热传导热阻，m·K/W；K_e为地层导热系数，W/(m·K)，f(t)为瑞曼时间函数。

7.如权利要求6所述的水平井同心双管注气隔热分析方法，其特征在于，所述的根据井口注入参数、所述井身结构热阻以及地层热物性参数确定井口处的最大热损失、单位长度的地层平均温度、隔热管向井壁的单位长度的热损失速度、内管的单位长度热损失速度包括：

根据地层热物性参数和式(12)确定单位长度对应的地层温度，

b_k＝(b_k-1+a1*dl/100) (12)

式中：b_k为地层温度，地层温度初始值b₀为地表温度，℃；a1为地层热物性参数中的地温梯度，℃/100m，dl为单位长度；

根据式(12-1)确定地层平均温度；

T_e＝(b_k+b_k-1)/2 (12-1)

根据井口注入蒸汽温度、所述地层平均温度、蒸汽与内管内壁之间的热对流热阻、内管内外壁热传导热阻、蒸汽与内管外壁之间的热对流热阻、蒸汽与隔热管内管内壁之间的热对流热阻、隔热管的内管内外壁之间的热传导热阻、隔热层的热传导热阻、隔热管外管内外壁之间的热传导热阻、套管壁的热传导热阻、水泥环的热传导热阻以及地层的热传导热阻以及式(17)确定井口处的最大热损失，其中，

$Q_m=\frac{T_s-T_e}{R_4+R_5+R_6+R_7+R_9+R_{10}+R_{11}}dl---\left(17\right)$

其中，Q_m为井口处的最大热损失，T_s为内管和隔热管环空的蒸汽平均温度，T_e为地层平均温度，式(17)中确定井口处的最大热损失时，T_e取在井口注入的蒸汽温度；

根据所述井口最大热损失确定竖井段的隔热管向井壁的单位长度的热损失速度为Q_k＝0.9Q_m；

根据内管的蒸汽平均温度、内管和隔热管环空的蒸汽平均温度以及式(33)确定竖井段的内管单位长度热损失速度；其中，

${Q_s}^{\′}=\frac{T_{s\_i}-T_s}{R_1+R_2+R_3}dl---\left(33\right)$

Q_s′为竖井段的内管单位长度热损失速度；T_{s_i}为内管的蒸汽平均温度。

8.如权利要求4所述的水平井同心双管注气隔热分析方法，其特征在于，所述的根据井身结构热阻、确定的单位长度的内管和隔热管环空的蒸汽平均温度、内管的蒸汽平均温度以及单位长度隔热管向井壁的热损失速度确定单位长度的井身结构温度包括：

根据直井段的井身结构热阻、单位长度的内管和隔热管环空的蒸汽平均温度、地层平均温度确定竖井段的隔热管内管内壁温度、隔热管内管外壁温度、隔热管外管内壁温度、隔热管外管外壁温度、水泥环外温度、套管外壁温度、套管内壁温度；

根据所述竖井段的单位长度的内管和隔热管环空的蒸汽平均温度、所述的内管的蒸汽平均温度、内管的单位长度热损失速度以及井身结构热阻确定水平井段的内管内壁温度、内管外壁温度、筛管外管温度以及筛管内管温度。

9.如权利要求8所述的水平井同心双管注气隔热分析方法，其特征在于，所述的根据确定的井身结构温度确定环空辐射传热系数和环空对流传热系数包括：

根据所述竖井段的隔热管外管外壁温度和套管内壁温度确定竖井段的环空辐射传热系数以及竖井段环空自然对流传热系数；

根据所述水平井段的内管外壁温度和筛管内管温度确定水平井段的环空辐射传热系数以及水平井段的自然对流传热系数。

10.如权利要求9所述的水平井同心双管注气隔热分析方法，其特征在于，所述的根据确定的环空辐射传热系数和环空对流传热系数确定环空对流热阻包括：

根据确定的竖井段的环空辐射传热系数、竖井段的环空自然对流传热系数以及式(8)确定竖井段的环空对流热阻，

$R_8=\frac{1}{2\mathrm{\π}(h_c+h_r)r_o}---\left(8\right)$

式中：R₈为竖井段的环空热对流热阻，m·K/W；hr为环空辐射传热系数，W/(m²·K)；hc为环空自然对流传热系数，W/(m²·K)，r_o为隔热管外管外壁半径，m；

根据所述水平井段的环空辐射传热系数、水平井段的自然对流传热系数以及式(47)、式(48)确定水平井段的环空对流热阻和筛管壁的热传导热阻，

水平段的环空热对流热阻R′₈：

$R_8^{\′}=\frac{1}{2\mathrm{\π}(h_c+h_r)r_{to\_i}}---\left(47\right)$

其中，hr为水平段的环空辐射传热系数，W/(m²·K)；hc为水平段环空自然对流传热系数，W/(m²·K)；r_{to_i}为内管外壁半径；

筛管壁的热传导热阻R′₉：

$R_9^{\′}=\frac{1}{2{\mathrm{\πK}}_{cas}^{\′}}ln\frac{r_{co}^{\′}}{r_{ci}^{\′}}---\left(48\right)$

式中：K′_cas为筛管导热系数，W/(m·K)；r′_ci为筛管内壁半径，m；r′_co为筛管外壁半径，m。

11.如权利要求10所述的水平井同心双管注气隔热分析方法，其特征在于，所述的根据确定的环空对流热阻对所述单位长度的隔热管向井壁的热损失速度进行校正生成校正后的隔热管向井壁的单位长度的热损失速度包括：

根据确定的竖井段的环空对流热阻和式(32)校正隔热管向井壁的单位长度热损失速度，其中，

$Q_k=\frac{T_s-T_e}{R}dl---\left(32\right)$

其中，R＝R₄+R₅+R₆+R₇+R₈+R₉+R₁₀+R₁₁，R₄为蒸汽与隔热管内管内壁之间的热对流热阻；R₅为隔热管内管内外壁之间的热传导热阻；R₆为绝热层的热传导热阻；R₇为隔热管外管内外壁之间的热传导热阻；R₉为套管壁的热传导热阻；R₁₀为水泥环的热传导热阻；R₁₁为地层的热传导热阻；dl为单位长度；T_s为内管和隔热管环空的蒸汽平均温度；T_e为地层平均温度；Q_k为竖井段的隔热管向井壁的单位长度的热损失速度。

12.如权利要求7或10所述的水平井同心双管注气隔热分析方法，其特征在于，根据井口注入参数、所述井身结构热阻以及地层热物性参数确定井口处的最大热损失、单位长度的地层平均温度、隔热管向井壁的单位长度的热损失速度、内管的单位长度热损失速度包括：

根据水平井段的内管的蒸汽平均温度、地层平均温度以及式(46)确定水平井段内管的热损失速度，其中，

$Q^{\′}=\frac{T_{s\_i}-T_e}{R_1+R_2+R_8^{\′}+R_9^{\′}+R_{11}}dl---\left(46\right)$

Q′为水平井段内管的热损失速度；T_{s_i}为内管的蒸汽平均温度；其中，

水平井段的环空热对流热阻R′₈：

$R_8^{\′}=\frac{1}{2\mathrm{\π}(h_c+h_r)r_{to\_i}}---\left(47\right)$

水平井段筛管壁的热传导热阻R′₉：

$R_9^{\′}=\frac{1}{2{\mathrm{\πK}}_{cas}^{\′}}ln\frac{r_{co}^{\′}}{r_{ci}^{\′}}---\left(48\right)$

式中：K′_cas为筛管导热系数，W/(m·K)；r′_ci为筛管内壁半径，m；r′_co为筛管外壁半径，m；hc为环空自然对流传热系数，W/(m²·K)；hr为环空辐射传热系数，W/(m²·K)；R₁为蒸汽与内管内壁之间的热对流热阻，m·K/W；R₂为内管内外壁之间的热传导热阻，m·K/W；R₁₁为地层的热传导热阻，m·K/W；dl为单位长度；T_e为地层平均温度。

13.一种水平井同心双管注气隔热分析装置，其特征在于，所述的装置包括：

参数采集模块，用于采集井身结构及热物理性质参数、井口注入参数以及地层热物性参数；

井身结构热阻确定模块，用于根据所述井身结构及热物理性质参数、井口注入参数以及地层热物性参数确定井身结构热阻；

温度干度确定模块，用于根据井口注入参数、所述井身结构热阻以及地层热物性参数确定井口处的最大热损失、单位长度的地层平均温度、隔热管向井壁的单位长度的热损失速度、内管的单位长度热损失速度、隔热管与内管环空的蒸汽干度以及内管的蒸汽干度；

压力确定模块，根据贝格斯-比尔算法确定单位长度的内管与隔热管的环空平均压力、单位长度的内管平均压力；

平均温度确定模块，用于根据所述的单位长度的内管与隔热管的环空平均压力、单位长度的内管平均压力确定单位长度的内管和隔热管环空的蒸汽平均温度、内管的蒸汽平均温度；

井身结构温度确定模块，用于根据所述单位长度的内管和隔热管环空的蒸汽平均温度、单位长度的地层平均温度以及内管的蒸汽平均温度确定井身结构温度；

传热系数确定模块，用于根据确定的井身结构温度确定环空辐射传热系数和环空对流传热系数；

环空对流热阻确定模块，用于根据确定的环空辐射传热系数和环空对流传热系数确定环空对流热阻；

校正热损失确定模块，根据确定的环空对流热阻对所述单位长度的隔热管向井壁的热损失速度进行校正生成校正后的隔热管向井壁的单位长度的热损失速度；

隔热分结果生成模块，根据所述校正后的隔热管向井壁的单位长度的热损失速度进行迭代分析生成水平井同心双管注气隔热分析结果。

14.如权利要求13所述的水平井同心双管注气隔热分析装置，其特征在于，所述的井身结构热阻包括：蒸汽与内管内壁之间的热对流热阻、内管内外壁热传导热阻、蒸汽与内管外壁之间的热对流热阻、蒸汽与隔热管内管内壁之间的热对流热阻、隔热管的内管内外壁之间的热传导热阻、隔热层的热传导热阻、隔热管外管内外壁之间的热传导热阻、套管壁的热传导热阻、水泥环的热传导热阻以及地层的热传导热阻。

15.如权利要求13所述的水平井同心双管注气隔热分析装置，其特征在于，所述隔热分结果生成模块中生成的水平井同心双管注气隔热分析结果包括：

生成直井段注汽隔热分析结果和水平井段注汽隔热分析结果。

16.如权利要求15所述的水平井同心双管注气隔热分析装置，其特征在于，所述的井身结构温度包括：隔热管内管内壁温度、隔热管内管外壁温度、隔热管外管内壁温度、隔热管外管外壁温度、水泥环外温度、套管外壁温度、套管内壁温度、内管内壁温度、筛管外管温度、筛管内管温度。

17.如权利要求13所述的水平井同心双管注气隔热分析装置，其特征在于，

所述的井身结构及热物理性质参数包括：水膜传热系数、水泥导热系数、井眼半径、内管内壁半径、内管外壁半径、内管下入深度、内管导热系数、隔热管下入深度、隔热管导热系数、隔热管内管内径、隔热管内管外径、隔热管外管外径、隔热管外管内径、套管内径、套管外径、筛管内径、筛管外径、油管导热系数、套管导热系数、筛管导热系数；

所述的井口注入参数包括：内管注汽压力；内管注汽干度；内管注汽速率；内管注汽时间；外管注汽压力；外管注汽干度；外管注汽速率；外管注汽时间；

所述的地层热物性参数包括：地温梯度；地表温度；地层导热系数。

18.如权利要求17所述的水平井同心双管注气隔热分析装置，其特征在于，所述的根据所述井身结构及热物理性质参数、井口注入参数以及地层热物性参数确定井身结构热阻包括：

根据所述的水膜传热系数、内管内壁半径以及式(1)确定蒸汽与内管内壁之间的热对流热阻；

$R_1=\frac{1}{2{\mathrm{\πh}}_f{r}_{ti\_i}}---\left(1\right)$

其中，R₁为蒸汽与内管内壁之间的热对流热阻，m·K/W；h_f为水膜传热系数，W/(m²·K)；r_{ti_i}为内管内壁半径，m；

根据所述油管导热系数、内管内壁半径、内管外壁半径及式(2)确定内管内外壁之间的热传导热阻；

$R_2=\frac{1}{2{\mathrm{\πK}}_{tub}}\ln \frac{r_{to\_i}}{r_{ti\_i}}---\left(2\right)$

式中：R₂为内管内外壁之间的热传导热阻，m·K/W；K_tub为油管导热系数，W/(m·K)；r_{to_i}为内管外壁半径，m；

所述的水膜传热系数、内管内壁半径以及式(3)确定蒸汽与内管外壁之间的热对流热阻；

$R_3=\frac{1}{2{\mathrm{\πh}}_f{r}_{to\_i}}---\left(3\right)$

式中：R₃为蒸汽与内管外壁之间的热对流热阻，m·K/W；

根据水膜传热系数、隔热管内管内壁半径以及式(4)确定蒸汽与隔热管内管内壁之间的热对流热阻；

$R_4=\frac{1}{2{\mathrm{\πh}}_f{r}_{ti}}---\left(4\right)$

式中：R₄为蒸汽与隔热管内管内壁之间的热对流热阻，m·K/W；r_ti为隔热管内管内壁半径，m；

根据油管导热系数、隔热管内管内壁半径、隔热管内管外壁半径和式(5)确定隔热管内管内外壁之间的热传导热阻：

$R_5=\frac{1}{2{\mathrm{\πK}}_{tub}}ln\frac{r_{to}}{r_{ti}}---\left(5\right)$

式中：R₅为隔热管内管内外壁之间的热传导热阻，m·K/W，K_tub为油管导热系数，W/(m·K)；r_to为隔热管内管外壁半径，m；

根据所述的隔热管导热系数、隔热管内管外壁半径、隔热管外管内壁半径以及式(6)确定隔热层的热传导热阻；

$R_6=\frac{1}{2{\mathrm{\πK}}_{ins}}ln\frac{r_i}{r_{to}}---\left(6\right)$

式中：R₆为绝热层的热传导热阻，m·K/W；K_ins为隔热管导热系数，W/(m·K)；r_i为隔热管外管内壁半径，m；

隔热管外管内外壁之间的热传导热阻：

$R_7=\frac{1}{2{\mathrm{\πK}}_{tub}}ln\frac{r_o}{r_i}---\left(7\right)$

式中：R₇为隔热管外管内外壁之间的热传导热阻，m·K/W，r_o为隔热管外管外壁半径，m，

根据所述套管导热系数、套管内壁半径、套管外壁半径以及式(9)确定套管壁的热传导热阻；

$R_9=\frac{1}{2{\mathrm{\πK}}_{cas}}ln\frac{r_{co}}{r_{ci}}---\left(9\right)$

式中：R₉为套管壁的热传导热阻，m·K/W，K_cas为套管导热系数，W/(m·K)，r_ci为套管内壁半径，m，r_co为套管外壁半径，m；

根据水泥环导热系数、井眼半径、套管外壁半径以及式(10)确定水泥环的热传导热阻：

$R_{10}=\frac{1}{2{\mathrm{\πK}}_{cem}}ln\frac{r_h}{r_{co}}---\left(10\right)$

式中：R₁₀为水泥环的热传导热阻，m·K/W；K_cem为水泥环导热系数，W/(m·K)；r_h为井眼半径，m；

根据瑞曼时间函数、地层导入系数以及式(11)确定地层的热传导热阻：

$R_{11}=\frac{f\left(t\right)}{2{\mathrm{\πK}}_e}---\left(11\right)$

式中：R₁₁为地层的热传导热阻，m·K/W；K_e为地层导热系数，W/(m·K)，f(t)为瑞曼时间函数。

19.如权利要求18所述的水平井同心双管注气隔热分析装置，其特征在于，所述的根据井口注入参数、所述井身结构热阻以及地层热物性参数确定井口处的最大热损失、单位长度的地层平均温度、隔热管向井壁的单位长度的热损失速度、内管的单位长度热损失速度包括：

根据地层热物性参数和式(12)确定单位长度对应的地层温度，

b_k＝(b_k-1+a1*dl/100) (12)

式中：b_k为地层温度，地层温度初始值b₀为地表温度，℃；a1为地层热物性参数中的地温梯度，℃/100m，dl为单位长度；

根据式(12-1)确定地层平均温度；

T_e＝(b_k+b_k-1)/2 (12-1)

根据井口注入蒸汽温度、所述地层平均温度、蒸汽与内管内壁之间的热对流热阻、内管内外壁热传导热阻、蒸汽与内管外壁之间的热对流热阻、蒸汽与隔热管内管内壁之间的热对流热阻、隔热管的内管内外壁之间的热传导热阻、隔热层的热传导热阻、隔热管外管内外壁之间的热传导热阻、套管壁的热传导热阻、水泥环的热传导热阻以及地层的热传导热阻以及式(17)确定井口处的最大热损失，其中，

$Q_m=\frac{T_s-T_e}{R_4+R_5+R_6+R_7+R_9+R_{10}+R_{11}}dl---\left(17\right)$

其中，Q_m为井口处的最大热损失，T_s为内管和隔热管环空的蒸汽平均温度，T_e为地层平均温度，式(17)中确定井口处的最大热损失时，T_e取在井口注入的蒸汽温度；

根据所述井口最大热损失确定竖井段的隔热管向井壁的单位长度的热损失速度为Q_k＝0.9Q_m；

根据内管的蒸汽平均温度、内管和隔热管环空的蒸汽平均温度以及式(33)确定竖井段的内管单位长度热损失速度；其中，

${Q_s}^{\′}=\frac{T_{s\_i}-T_s}{R_1+R_2+R_3}dl---\left(33\right)$

Q_s′为竖井段的内管单位长度热损失速度；T_{s_i}为内管的蒸汽平均温度。

20.如权利要求16所述的水平井同心双管注气隔热分析装置，其特征在于，所述的根据井身结构热阻、确定的单位长度的内管和隔热管环空的蒸汽平均温度、内管的蒸汽平均温度以及单位长度隔热管向井壁的热损失速度确定单位长度的井身结构温度包括：

根据直井段的井身结构热阻、单位长度的内管和隔热管环空的蒸汽平均温度、地层平均温度确定竖井段的隔热管内管内壁温度、隔热管内管外壁温度、隔热管外管内壁温度、隔热管外管外壁温度、水泥环外温度、套管外壁温度、套管内壁温度；

根据所述竖井段的单位长度的内管和隔热管环空的蒸汽平均温度、所述的内管的蒸汽平均温度、内管的单位长度热损失速度以及井身结构热阻确定水平井段的内管内壁温度、内管外壁温度、筛管外管温度以及筛管内管温度。

21.如权利要求20所述的水平井同心双管注气隔热分析装置，其特征在于，所述的根据确定的井身结构温度确定环空辐射传热系数和环空对流传热系数包括：

根据所述竖井段的隔热管外管外壁温度和套管内壁温度确定竖井段的环空辐射传热系数以及竖井段环空自然对流传热系数；

根据所述水平井段的内管外壁温度和筛管内管温度确定水平井段的环空辐射传热系数以及水平井段的自然对流传热系数。

22.如权利要求21所述的水平井同心双管注气隔热分析装置，其特征在于，所述的根据确定的环空辐射传热系数和环空对流传热系数确定环空对流热阻包括：

根据确定的竖井段的环空辐射传热系数、竖井段的环空自然对流传热系数以及式(8)确定竖井段的环空对流热阻，

$R_8=\frac{1}{2\mathrm{\π}(h_c+h_r)r_o}---\left(8\right)$

式中：R₈为竖井段的环空热对流热阻，m·K/W；hr为环空辐射传热系数，W/(m²·K)；hc为环空自然对流传热系数，W/(m²·K)，r_o为隔热管外管外壁半径，m；

根据所述水平井段的环空辐射传热系数、水平井段的自然对流传热系数以及式(47)、式(48)确定水平井段的环空对流热阻和筛管壁的热传导热阻，

水平段的环空热对流热阻R′₈：

$R_8^{\′}=\frac{1}{2\mathrm{\π}(h_c+h_r)r_{to\_i}}---\left(47\right)$

其中，hr为水平段的环空辐射传热系数，W/(m²·K)；hc为水平段环空自然然对流传热系数，W/(m²·K)；r_{to_i}为内管外壁半径；

筛管壁的热传导热阻R′₉：

$R_9^{\′}=\frac{1}{2{\mathrm{\πK}}_{cas}^{\′}}ln\frac{r_{co}^{\′}}{r_{ci}^{\′}}---\left(48\right)$

式中：K′_cas为筛管导热系数，W/(m·K)；r′_ci为筛管内壁半径，m；r′_co为筛管外壁半径，m。

23.如权利要求22所述的水平井同心双管注气隔热分析装置，其特征在于，所述的根据确定的环空对流热阻对所述单位长度的隔热管向井壁的热损失速度进行校正生成校正后的隔热管向井壁的单位长度的热损失速度包括：

根据确定的竖井段的环空对流热阻和式(32)校正隔热管向井壁的单位长度热损失速度，其中，

$Q_k=\frac{T_s-T_e}{R}dl---\left(32\right)$

其中，R＝R₄+R₅+R₆+R₇+R₈+R₉+R₁₀+R₁₁，R₄为蒸汽与隔热管内管内壁之间的热对流热阻；R₅为隔热管内管内外壁之间的热传导热阻；R₆为绝热层的热传导热阻；R₇为隔热管外管内外壁之间的热传导热阻；R₉为套管壁的热传导热阻；R₁₀为水泥环的热传导热阻；R₁₁为地层的热传导热阻；dl为单位长度；T_s为内管和隔热管环空的蒸汽平均温度；T_e为地层平均温度；Q_k为竖井段的隔热管向井壁的单位长度的热损失速度。

24.如权利要求19或22所述的水平井同心双管注气隔热分析装置，其特征在于，根据井口注入参数、所述井身结构热阻以及地层热物性参数确定井口处的最大热损失、单位长度的地层平均温度、隔热管向井壁的单位长度的热损失速度、内管的单位长度热损失速度包括：

根据水平井段的内管的蒸汽平均温度、地层平均温度以及式(46)确定水平井段内管的热损失速度，其中，

$Q^{\′}=\frac{T_{s\_i}-T_e}{R_1+R_2+R_8^{\′}+R_9^{\′}+R_{11}}dl---\left(46\right)$

Q′为水平井段内管的热损失速度；T_{s_i}为内管的蒸汽平均温度；其中，

水平井段的环空热对流热阻R′₈：

$R_8^{\′}=\frac{1}{2\mathrm{\π}(h_c+h_r)r_{to\_i}}---\left(47\right)$

水平井段筛管壁的热传导热阻R′₉：

$R_9^{\′}=\frac{1}{2{\mathrm{\πK}}_{cas}^{\′}}ln\frac{r_{co}^{\′}}{r_{ci}^{\′}}---\left(48\right)$

式中：K′_cas为筛管导热系数，W/(m·K)；r′_ci为筛管内壁半径，m；r′_co为筛管外壁半径，m；hc为环空自然对流传热系数，W/(m²·K)；hr为环空辐射传热系数，W/(m²·K)；R₁为蒸汽与内管内壁之间的热对流热阻，m·K/W；R₂为内管内外壁之间的热传导热阻，m·K/W；R₁₁为地层的热传导热阻，m·K/W；dl为单位长度；T_e为地层平均温度。