CN112697313A - 一种注蒸汽井筒传热过程中的热损失量评估方法及装置 - Google Patents

一种注蒸汽井筒传热过程中的热损失量评估方法及装置 Download PDF

Info

Publication number
CN112697313A
CN112697313A CN202011422405.6A CN202011422405A CN112697313A CN 112697313 A CN112697313 A CN 112697313A CN 202011422405 A CN202011422405 A CN 202011422405A CN 112697313 A CN112697313 A CN 112697313A
Authority
CN
China
Prior art keywords
heat transfer
transfer coefficient
temperature
coefficient
oil pipe
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.)
Granted
Application number
CN202011422405.6A
Other languages
English (en)
Other versions
CN112697313B (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.)
Yangtze University
Original Assignee
Yangtze University
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 Yangtze University filed Critical Yangtze University
Priority to CN202011422405.6A priority Critical patent/CN112697313B/zh
Publication of CN112697313A publication Critical patent/CN112697313A/zh
Application granted granted Critical
Publication of CN112697313B publication Critical patent/CN112697313B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K17/00—Measuring quantity of heat
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N25/00—Investigating or analyzing materials by the use of thermal means
    • G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/10—Geothermal energy

Landscapes

  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明涉及一种注蒸汽井筒传热过程中的热损失量评估方法、装置及计算机可读存储介质，所述方法包括以下步骤：获取井眼半径及地层热扩散系数，根据所述井眼半径及地层热扩散系数，获取无因次地层时间函数；获取油管的外壁温度、套管的内壁温度及环空液体的等效导热系数，根据油管的外壁温度、套管的内壁温度及环空液体的等效导热系数获取总传热系数；获取油管外半径、蒸汽温度及水泥环外缘温度，根据所述油管外半径、蒸汽温度、水泥环外缘温度、无因次地层时间函数及总传热系数，得到井筒传热过程中的热损失量。本发明所述的注蒸汽井筒传热过程中的热损失量评估方法，实现了对注蒸汽井筒传热过程中的热损失量的评估。

Description

一种注蒸汽井筒传热过程中的热损失量评估方法及装置
技术领域
本发明涉及注蒸汽采油技术领域，尤其涉及一种注蒸汽井筒传热过程中的热损失量评估方法、装置及计算机可读存储介质。
背景技术
世界上的稠油资源非常丰富，其地质储量远远超过常规原油储量。到目前为止，已探明储量达三千亿吨，而稀油资源只有1700亿吨可供开采，稠油是二十一世纪的重要资源。中国的稠油资源也十分丰富，分布广泛，目前中国已经成为四大稠油开采国之一，稠油资源是我国一项巨大的潜在资源，在今后的能源结构中起着至关重要的作用。
目前注蒸汽热采技术在国内外稠油油田被广泛地应用，由锅炉产生的高温高压的水蒸汽经地面管线和井筒进入油藏的过程中，必然伴随着热量损失，直接影响热采的效果。蒸汽从地面管道流出后，就会进入井筒直到最后进入各射孔层，在进入射孔层之前，蒸汽在井筒中还是会沿着井筒损失热量，这样就会引起蒸汽干度的变化，进入各射孔层的热量也会变化，因此，需要评估注汽过程中，蒸汽在井筒中损失热量的情况，从而判断吸汽效果的好坏，指导和调整注汽参数，为油田适时调整开发方案提供参考。现有技术中缺少评估注蒸汽井筒传热过程中的热损失量的方案。
发明内容
有鉴于此，有必要提供一种注蒸汽井筒传热过程中的热损失量评估方法、装置及计算机可读存储介质，用以解决现有技术中缺少评估注蒸汽井筒传热过程中的热损失量方案的问题。
本发明提供一种注蒸汽井筒传热过程中的热损失量评估方法，包括以下步骤：
获取井眼半径及地层热扩散系数，根据所述井眼半径及地层热扩散系数，获取无因次地层时间函数；
获取油管的外壁温度、套管的内壁温度及环空液体的等效导热系数，根据油管的外壁温度、套管的内壁温度及环空液体的等效导热系数获取总传热系数；
获取油管外半径、蒸汽温度及水泥环外缘温度，根据所述油管外半径、蒸汽温度、水泥环外缘温度、无因次地层时间函数及总传热系数，得到井筒传热过程中的热损失量。
进一步地，根据所述井眼半径及地层热扩散系数，获取无因次地层时间函数，具体包括，根据所述井眼半径、地层热扩散系数及无因次地层时间函数计算公式获取无因次地层时间函数，所述无因次地层时间函数计算公式为
Figure BDA0002823026710000021
其中，
Figure BDA0002823026710000022
α为地层热扩散系数，τ为计算时间，rh为井眼半径。
进一步地，根据油管的外壁温度、套管的内壁温度及环空液体的等效导热系数获取总传热系数，具体包括，根据所述油管的外壁温度及套管的内壁温度，确定环空辐射换热系数，根据环空液体的等效导热系数，确定自然对流传热系数，根据所述环空辐射换热系数及自然对流传热系数，获取总传热系数。
进一步地，根据所述油管的外壁温度及套管的内壁温度，确定环空辐射换热系数，具体包括，根据所述油管的外壁温度、套管的内壁温度及辐射换热系数公式，确定环空辐射换热系数，所述辐射换热系数公式为hr＝σFtci(Tto *2+Tci *2)(Tto *+Tci *)，其中，*为绝对温度，σ为常数，Ftci为油管外表面向套管内表面辐射散热有效系数，Tto为油管的外壁温度，Tci为套管的内壁温度。
进一步地，根据环空液体的等效导热系数，确定自然对流传热系数，具体包括，根据所述环空液体的等效导热系数及自然对流传热系数公式，确定自然对流传热系数，所述自然对流传热系数公式为
Figure BDA0002823026710000031
其中，Khc为环空液体的等效导热系数，rto为油管的外半径，rci为套管的内半径。
进一步地，根据所述环空辐射换热系数及自然对流传热系数，获取总传热系数，具体包括，根据所述环空辐射换热系数、自然对流传热系数及总传热系数计算公式，获取总传热系数，所述总传热系数计算公式为
Figure BDA0002823026710000032
其中hc、hr分别为环空辐射换热系数、自然对流传热系数,rco为套管的外半径，Kcem水泥环的导热系数,rh井眼半径。
进一步地，根据所述油管外半径、蒸汽温度、水泥环外缘温度、无因次地层时间函数及总传热系数，得到井筒传热过程中热损失量，具体包括，根据所述油管外半径、蒸汽温度、水泥环外缘温度、无因次地层时间函数及总传热系数，得到井筒内微元段dz的径向传热量，由井筒内微元段的径向传热量获取井筒传热过程中的热损失量。
进一步地，根据所述油管外半径、蒸汽温度、水泥环外缘温度、无因次地层时间函数、总传热系数及传热公式，得到井筒内微元段dz的径向传热量，所述传热公式为
Figure BDA0002823026710000041
其中，Ke地层导热系数，Te为地层温度，rto为油管外半径，Ts为蒸汽温度，Th为水泥环外缘温度，dQ为井筒内微元段dz的径向传热量。
本发明还提供了一种注蒸汽井筒传热过程中的热损失量评估装置，包括处理器以及存储器，所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行时，实现如上述任一技术方案所述的注蒸汽井筒传热过程中的热损失量评估方法。
本发明还提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机该程序被处理器执行时，实现如上述任一技术方案所述的注蒸汽井筒传热过程中的热损失量评估方法。
与现有技术相比，本发明的有益效果包括：通过获取井眼半径及地层热扩散系数，根据所述井眼半径及地层热扩散系数，获取无因次地层时间函数；获取油管的外壁温度、套管的内壁温度及环空液体的等效导热系数，根据油管的外壁温度、套管的内壁温度及环空液体的等效导热系数获取总传热系数；获取油管外半径、蒸汽温度及水泥环外缘温度，根据所述油管外半径、蒸汽温度、水泥环外缘温度、无因次地层时间函数及总传热系数，得到井筒传热过程中热损失量；实现了对注蒸汽井筒传热过程中的热损失量的评估。
附图说明
图1为本发明提供的注蒸汽井筒传热过程中的热损失量评估方法的流程示意图；
图2为本发明提供的井筒的径向结构示意图。
附图标记：1-油管；2-环空；3-套管；4-水泥环；5-地层。
具体实施方式
下面结合附图来具体描述本发明的优选实施例，其中，附图构成本申请一部分，并与本发明的实施例一起用于阐释本发明的原理，并非用于限定本发明的范围。
实施例1
本发明实施例提供了一种注蒸汽井筒传热过程中的热损失量评估方法，其流程示意图，如图1所示，其包括以下步骤：
S1、获取井眼半径及地层热扩散系数，根据所述井眼半径及地层热扩散系数，获取无因次地层时间函数；
S2、获取油管的外壁温度、套管的内壁温度及环空液体的等效导热系数，根据油管的外壁温度、套管的内壁温度及环空液体的等效导热系数获取总传热系数；
S3、获取油管外半径、蒸汽温度及水泥环外缘温度，根据所述油管外半径、蒸汽温度、水泥环外缘温度、无因次地层时间函数及总传热系数，得到井筒传热过程中的热损失量。
一个具体实施例中，井筒的径向结构示意图，如图2所示，注蒸汽井筒主要包括由油管、套管和水泥环；
需要说明的是，蒸汽的注入一般有两种方式，一种是从套管注入，另一种是从油管注入；当从套管注入时，由于蒸汽会直接给套管壁加热，热量损失会更多；而从油管注入时，由于油管和套管之间的环空区域充满隔热介质，热量损失会减少，一般采用第二种方法，在环空中一般会充填空气、氮气等气体，可以在环空中间加入隔热管，以最好的隔热效果；
具体实施时，把井筒中蒸汽的径向传热看成是两部分组成，即油管中心到水泥环外缘的一维稳定传热和水泥环外缘到地层之间的一维不稳定传热，同时径向热损失是沿井深方向，并随时间变化的，因此，井筒热损失的计算必须在井筒某一深度和时间上分段进行，假设在单位时间内，井筒某段长度dz上的热损失为dQ；
由于是稳定传热，可知井筒在dz段内蒸汽从油管中心到水泥环外缘在传热量为dQ＝2πrtoUto(Ts-Th)dz，其中，Uto为总传热系数，W/(m2·℃)，rto为油管外半径，m，Ts为蒸汽温度，℃，Th为水泥环外缘温度，℃；这一传热量等于由水泥环外缘至未受影响的地层的不稳定导热量；由于是非稳定的热传导，热量随时间发生变化，则有
Figure BDA0002823026710000061
其中，Ke为地层导热系数，W/(m·℃)，Te为地层温度，℃；f(τD)为无因次地层时间函数；
优选的，根据所述井眼半径及地层热扩散系数，获取无因次地层时间函数，具体包括，根据所述井眼半径、地层热扩散系数及无因次地层时间函数计算公式获取无因次地层时间函数，所述无因次地层时间函数计算公式为
Figure BDA0002823026710000062
其中，
Figure BDA0002823026710000063
α为地层热扩散系数(m2/s)，τ为计算时间(s)，rh为井眼半径(m)；
优选的，根据油管的外壁温度、套管的内壁温度及环空液体的等效导热系数获取总传热系数，具体包括，根据所述油管的外壁温度及套管的内壁温度，确定环空辐射换热系数，根据环空液体的等效导热系数，确定自然对流传热系数，根据所述环空辐射换热系数及自然对流传热系数，获取总传热系数。
优选的，根据所述油管的外壁温度及套管的内壁温度，确定环空辐射换热系数，具体包括，根据所述油管的外壁温度、套管的内壁温度及辐射换热系数公式，确定环空辐射换热系数，所述辐射换热系数公式为hr＝σFtci(Tto *2+Tci *2)(Tto *+Tci *)，其中，*为绝对温度，σ为常数，Ftci为油管外表面向套管内表面辐射散热有效系数，Tto为油管的外壁温度，Tci为套管的内壁温度；
一个具体实施例中，当油套环空中有气体时，辐射热量Qr由油管外壁的温度Tto与套管内壁温度Tci决定；则有
Figure BDA0002823026710000071
其中，*为绝对温度(T+273)；
σ为Stefan-Boltzman常数[1.173*10-19Btu/(ft2·h·℃)]；
Ftci为油管外表面Ato向套管内表面Aci辐射散热有效系数，它表示辐射吸收能力；Tto为油管的外壁温度；Tci为套管的内壁温度；对于井筒条件，
Figure BDA0002823026710000072
其中，εto为油管外壁辐射系数，无因次，εci为套管内壁辐射系数，无因次，
Figure BDA0002823026710000073
为两个表面间的总交换系数；对于井简传热条件，
Figure BDA0002823026710000074
取值为1.0；因此上式可简化为
Figure BDA0002823026710000075
根据式
Figure BDA0002823026710000076
因子分解，可得出辐射换热系数hr的计算公式
Figure BDA0002823026710000077
由于在上面已经计算出Tci，而Tto在计算时近似把它当做某一深度点的蒸汽温度，所以hr可以用hr的计算公式计算；
优选的，根据环空液体的等效导热系数，确定自然对流传热系数，具体包括，根据所述环空液体的等效导热系数及自然对流传热系数公式，确定自然对流传热系数，所述自然对流传热系数公式为
Figure BDA0002823026710000081
其中，Khc为环空液体的等效导热系数，rto为油管的外半径，rci为套管的内半径；
一个具体实施例中，油套环空的热传导及自然对流引起的径向传热速度为
Figure BDA0002823026710000082
其中，Qc为热传导及自然对流引起的径向传热速度，w；Khc为环空液体的等效导热系数，即在环空的平均温度和压力下，自然对流影响的环空液体的综合导热系数，W/(m2·℃)；当自然对流很小时，Khc＝Kha，Kha是环空液体或者气体的导热系数；因为
Qc＝2πrtohc(Tci-Tto)ΔL
所以hc可以由下面的公式确定：
Figure BDA0002823026710000083
在井筒条件下有
Figure BDA0002823026710000084
上式中的格拉绍夫数Gr和普朗特常数Pr分别由下面的计算式决定：
Figure BDA0002823026710000085
Figure BDA0002823026710000091
其中，g为重力加速度，1.2718×108m/h；ρan为流体在平均温度和压力下的密度，kg/m3；μan为在平均温度和压力下环空流体的粘度，kg/(m2·s)；β为环空流体的体积膨胀系数；Can为在平均温度和压力下环空流体的比热容，kJ/(kg·℃)；
公式Qc＝2πrtohc(Tci-Tto)ΔL的有效范围是5×104<GrPr<7.17×108，环空中为高压气体时，GrPr乘积在105～109之间；对于上面环空流体的参数，可以根据实际流体的物性表进行插值或者拟合；
9、优选的，根据所述环空辐射换热系数及自然对流传热系数，获取总传热系数，具体包括，根据所述环空辐射换热系数、自然对流传热系数及总传热系数计算公式，获取总传热系数，所述总传热系数计算公式为
Figure BDA0002823026710000092
其中，hc、hr分别为环空辐射换热系数、自然对流传热系数,rco为套管的外半径，Kcem水泥环的导热系数,rh井眼半径；
需要说明的是，当井筒中只是单层油管，下部带有封隔器，油管和套管之间的环形空间是液体或者气体时，井筒总的传热系数可由下式计算
Figure BDA0002823026710000093
上式中括号内各项依次表示油管内壁强迫对流传热热阻，油管壁热阻，环空液体或气体的热阻，套管壁热阻及水泥环的热阻；
由于井筒流体和油管内表面之间的传热系数hf是非常大的，其热阻基本可以忽略不计，油管和套管的导热系数Ktub和Kcas要比水泥环和地层的导热系数大得多，且油管壁和套管壁一般都很薄，因此他们的热阻很小，也可以忽略不计；这样可以把上式式简化为
Figure BDA0002823026710000101
其中，hc为环空气体或液体热传导和自然对流的传热系数，W/(m2·℃)；hr为环空气体或液体辐射传热系数，W/(m2·℃)；Kcem为水泥环的导热系数，W/(m2·℃)；rco为套管的外半径，m；rto为油管的外半径，m；当井筒中油管柱是双层隔热管，下端有封隔器，环空是液体或者气体时，井筒的总的传热系数由下面的公式计算
Figure BDA0002823026710000102
这个公式也可以简化成
Figure BDA0002823026710000103
其中，ri为隔热油管的外半径，m；Kins为隔热油管的导热系数，W/(m2·℃)；hc为环空气体或液体热传导和自然对流传热系数，W/(m2·℃)；hr＇为环空气体或液体辐射传热系数，W/(m2·℃)；这几个参数均与油管的外表面性质、液体的物理性质、油管外壁与套管内壁之间的温度与距离、套管内壁表面性质等有关；
具体实施例时，欲求某一时刻t下，Δz管长上的热损失ΔQ，则必须求得Th值，而Th值的大小又与Uto的取值有关，设定已知Uto，则可计算出Th及ΔQ值，从而得到新的计算值Ut'o，这样即可以将作为Uto迭代参数进行选代计；
Uto具体迭代计算的步骤如下：根据设定一个Uto值，计算f(τD)，计算出Th值，计算出Tci值，估算hr及hc，求出新的Uto值；一般情况下，计算出的Uto值与最初设定的Uto值会有差别，利用计算出的Uto，重复计算得到最新Uto值，这样进行迭代若干就可求准确的Uto值，从而确定最终的Uto值以后；
优选的，根据所述油管外半径、蒸汽温度、水泥环外缘温度、无因次地层时间函数及总传热系数，得到井筒传热过程中热损失量，具体包括，根据所述油管外半径、蒸汽温度、水泥环外缘温度、无因次地层时间函数及总传热系数，得到井筒内微元段dz的径向传热量，由井筒内微元段的径向传热量获取井筒传热过程中的热损失量。
优选的，根据所述油管外半径、蒸汽温度、水泥环外缘温度、无因次地层时间函数、总传热系数及传热公式，得到井筒内微元段dz的径向传热量，所述传热公式为
Figure BDA0002823026710000111
其中，Ke地层导热系数，Te为地层温度，rto为油管外半径，Ts为蒸汽温度，Th为水泥环外缘温度，dQ为井筒内微元段dz的径向传热量；
一个具体实施例中，由于井筒中热传递具有连续性，那么蒸汽从油管中心至水泥环外缘传递的热量等于水泥环外缘到地层传热量，因此水泥环外缘温度可表示为
Figure BDA0002823026710000112
同时可得到套管内表面温度
Figure BDA0002823026710000121
因此，可得井筒内一微元段dz的径向传热量为
Figure BDA0002823026710000122
实施例2
本发明实施例提供了一种注蒸汽井筒传热过程中的热损失量评估装置，包括处理器以及存储器，所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行时，实现如实施例1所述的注蒸汽井筒传热过程中的热损失量评估方法。
实施例3
本发明实施例提供了一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机该程序被处理器执行时，实现如实施例1所述的注蒸汽井筒传热过程中的热损失量评估方法。
本发明公开了一种注蒸汽井筒传热过程中的热损失量评估方法、装置及计算机可读存储介质；通过获取井眼半径及地层热扩散系数，根据所述井眼半径及地层热扩散系数，获取无因次地层时间函数；获取油管的外壁温度、套管的内壁温度及环空液体的等效导热系数，根据油管的外壁温度、套管的内壁温度及环空液体的等效导热系数获取总传热系数；获取油管外半径、蒸汽温度及水泥环外缘温度，根据所述油管外半径、蒸汽温度、水泥环外缘温度、无因次地层时间函数及总传热系数，得到井筒传热过程中热损失量；实现了对注蒸汽井筒传热过程中的热损失量的评估。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。

Claims (10)

1.一种注蒸汽井筒传热过程中的热损失量评估方法，其特征在于，包括以下步骤：
获取井眼半径及地层热扩散系数，根据所述井眼半径及地层热扩散系数，获取无因次地层时间函数；
获取油管的外壁温度、套管的内壁温度及环空液体的等效导热系数，根据油管的外壁温度、套管的内壁温度及环空液体的等效导热系数获取总传热系数；
获取油管外半径、蒸汽温度及水泥环外缘温度，根据所述油管外半径、蒸汽温度、水泥环外缘温度、无因次地层时间函数及总传热系数，得到井筒传热过程中的热损失量。
2.根据权利要求1所述的注蒸汽井筒传热过程中的热损失量评估方法，其特征在于，根据所述井眼半径及地层热扩散系数，获取无因次地层时间函数，具体包括，根据所述井眼半径、地层热扩散系数及无因次地层时间函数计算公式获取无因次地层时间函数，所述无因次地层时间函数计算公式为
Figure FDA0002823026700000011
其中，
Figure FDA0002823026700000012
α为地层热扩散系数，τ为计算时间，rh为井眼半径。
3.根据权利要求2所述的注蒸汽井筒传热过程中的热损失量评估方法，其特征在于，根据油管的外壁温度、套管的内壁温度及环空液体的等效导热系数获取总传热系数，具体包括，根据所述油管的外壁温度及套管的内壁温度，确定环空辐射换热系数，根据环空液体的等效导热系数，确定自然对流传热系数，根据所述环空辐射换热系数及自然对流传热系数，获取总传热系数。
4.根据权利要求3所述的注蒸汽井筒传热过程中的热损失量评估方法，其特征在于，根据所述油管的外壁温度及套管的内壁温度，确定环空辐射换热系数，具体包括，根据所述油管的外壁温度、套管的内壁温度及辐射换热系数公式，确定环空辐射换热系数，所述辐射换热系数公式为hr＝σFtci(Tto *2+Tci *2)(Tto *+Tci *)，其中，*为绝对温度，σ为常数，Ftci为油管外表面向套管内表面辐射散热有效系数，Tto为油管的外壁温度，Tci为套管的内壁温度。
5.根据权利要求4所述的注蒸汽井筒传热过程中的热损失量评估方法，其特征在于，根据环空液体的等效导热系数，确定自然对流传热系数，具体包括，根据所述环空液体的等效导热系数及自然对流传热系数公式，确定自然对流传热系数，所述自然对流传热系数公式为
Figure FDA0002823026700000021
其中，Khc为环空液体的等效导热系数，rto为油管的外半径，rci为套管的内半径。
6.根据权利要求5所述的注蒸汽井筒传热过程中的热损失量评估方法，其特征在于，根据所述环空辐射换热系数及自然对流传热系数，获取总传热系数，具体包括，根据所述环空辐射换热系数、自然对流传热系数及总传热系数计算公式，获取总传热系数，所述总传热系数计算公式为
Figure FDA0002823026700000022
其中，hc、hr分别为环空辐射换热系数、自然对流传热系数,rco为套管的外半径，Kcem水泥环的导热系数,rh井眼半径。
7.根据权利要求6所述的注蒸汽井筒传热过程中的热损失量评估方法，其特征在于，根据所述油管外半径、蒸汽温度、水泥环外缘温度、无因次地层时间函数及总传热系数，得到井筒传热过程中热损失量，具体包括，根据所述油管外半径、蒸汽温度、水泥环外缘温度、无因次地层时间函数及总传热系数，得到井筒内微元段dz的径向传热量，由井筒内微元段的径向传热量获取井筒传热过程中的热损失量。
8.根据权利要求7所述的注蒸汽井筒传热过程中的热损失量评估方法，其特征在于，根据所述油管外半径、蒸汽温度、水泥环外缘温度、无因次地层时间函数、总传热系数及传热公式，得到井筒内微元段dz的径向传热量，所述传热公式为
Figure FDA0002823026700000031
其中，Ke地层导热系数，Te为地层温度，rto为油管外半径，Ts为蒸汽温度，Th为水泥环外缘温度，dQ为井筒内微元段dz的径向传热量。
9.一种注蒸汽井筒传热过程中的热损失量评估装置，其特征在于，包括处理器以及存储器，所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行时，实现如权利要求1-8任一所述的注蒸汽井筒传热过程中的热损失量评估方法。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机该程序被处理器执行时，实现如权利要求1-8任一所述的注蒸汽井筒传热过程中的热损失量评估方法。
CN202011422405.6A 2020-12-08 2020-12-08 一种注蒸汽井筒传热过程中的热损失量评估方法及装置 Active CN112697313B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011422405.6A CN112697313B (zh) 2020-12-08 2020-12-08 一种注蒸汽井筒传热过程中的热损失量评估方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011422405.6A CN112697313B (zh) 2020-12-08 2020-12-08 一种注蒸汽井筒传热过程中的热损失量评估方法及装置

Publications (2)

Publication Number Publication Date
CN112697313A true CN112697313A (zh) 2021-04-23
CN112697313B CN112697313B (zh) 2023-07-18

Family

ID=75506890

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011422405.6A Active CN112697313B (zh) 2020-12-08 2020-12-08 一种注蒸汽井筒传热过程中的热损失量评估方法及装置

Country Status (1)

Country Link
CN (1) CN112697313B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2424420C1 (ru) * 2010-02-01 2011-07-20 Открытое акционерное общество "Газпром" Способ определения коэффициента теплопроводности теплоизоляции теплоизолированной лифтовой трубы в скважине
CN202956353U (zh) * 2012-11-26 2013-05-29 中国石油天然气股份有限公司 热水驱采油井筒环空介质隔热效果实验装置
CN103775058A (zh) * 2013-12-31 2014-05-07 中国石油天然气股份有限公司 一种井筒热损失的确定方法
CN103899295A (zh) * 2014-03-10 2014-07-02 中国石油天然气股份有限公司 稠油蒸汽采注中直井井筒热损失检测方法及检测装置
CN103900737A (zh) * 2014-03-10 2014-07-02 中国石油天然气股份有限公司 稠油蒸汽采注中地面管线热损失检测方法及检测装置
CN104850681A (zh) * 2015-04-16 2015-08-19 中国石油天然气股份有限公司 注入的蒸汽、氮气混合物的温度、压力、干度的计算方法
CN105160161A (zh) * 2015-08-17 2015-12-16 中国石油天然气股份有限公司 井筒内热力参数计算方法及装置
CN105370255A (zh) * 2015-12-01 2016-03-02 中国石油天然气股份有限公司 火烧油层笼统注气电点火井筒温度分布的确定方法及装置
CN105403590A (zh) * 2015-11-12 2016-03-16 中国石油天然气股份有限公司 隔热管导热系数测试方法及装置
CN106018472A (zh) * 2016-07-08 2016-10-12 中国石油天然气股份有限公司 注汽井隔热管导热系数测试方法及装置
CN106441637A (zh) * 2016-08-31 2017-02-22 中国石油天然气股份有限公司 确定注汽管热损失的方法及装置

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2424420C1 (ru) * 2010-02-01 2011-07-20 Открытое акционерное общество "Газпром" Способ определения коэффициента теплопроводности теплоизоляции теплоизолированной лифтовой трубы в скважине
CN202956353U (zh) * 2012-11-26 2013-05-29 中国石油天然气股份有限公司 热水驱采油井筒环空介质隔热效果实验装置
CN103775058A (zh) * 2013-12-31 2014-05-07 中国石油天然气股份有限公司 一种井筒热损失的确定方法
CN103899295A (zh) * 2014-03-10 2014-07-02 中国石油天然气股份有限公司 稠油蒸汽采注中直井井筒热损失检测方法及检测装置
CN103900737A (zh) * 2014-03-10 2014-07-02 中国石油天然气股份有限公司 稠油蒸汽采注中地面管线热损失检测方法及检测装置
CN104850681A (zh) * 2015-04-16 2015-08-19 中国石油天然气股份有限公司 注入的蒸汽、氮气混合物的温度、压力、干度的计算方法
CN105160161A (zh) * 2015-08-17 2015-12-16 中国石油天然气股份有限公司 井筒内热力参数计算方法及装置
CN105403590A (zh) * 2015-11-12 2016-03-16 中国石油天然气股份有限公司 隔热管导热系数测试方法及装置
CN105370255A (zh) * 2015-12-01 2016-03-02 中国石油天然气股份有限公司 火烧油层笼统注气电点火井筒温度分布的确定方法及装置
CN106018472A (zh) * 2016-07-08 2016-10-12 中国石油天然气股份有限公司 注汽井隔热管导热系数测试方法及装置
CN106441637A (zh) * 2016-08-31 2017-02-22 中国石油天然气股份有限公司 确定注汽管热损失的方法及装置

Also Published As

Publication number Publication date
CN112697313B (zh) 2023-07-18

Similar Documents

Publication Publication Date Title
CN109598099B (zh) 一种考虑油藏与井筒耦合的双管sagd长水平井均匀注汽数值模拟方法
Izgec et al. Transient fluid and heat flow modeling in coupled wellbore/reservoir systems
CN105160161B (zh) 井筒内热力参数确定方法及装置
CN102272418B (zh) 用于估计sagd过程特性的方法
CN104453861B (zh) 一种高压气井井筒温度分布的确定方法以及系统
CN104806230A (zh) 超临界注汽井的井筒温度场及压力场分布的计算方法
CN106968667A (zh) 一种温度场预测方法及装置
CN103775058A (zh) 一种井筒热损失的确定方法
CN103899295A (zh) 稠油蒸汽采注中直井井筒热损失检测方法及检测装置
CN105422084B (zh) 一种热水循环加热降粘工艺的井筒温度场获取方法及装置
CN102682195B (zh) 半潜式平台瞬态钻井井筒温度计算方法
CN106368633A (zh) 深水油气生产井环空密闭压力调控管柱系统
Yasunami et al. CO2 temperature prediction in injection tubing considering supercritical condition at Yubari ECBM Pilot-Test
CN102777159B (zh) 一种注co2气井井筒流态确定及参数优化方法
Benson et al. Nonisothermal effects during injection and falloff tests
CN109033012B (zh) 一种空心抽油杆注热水循环井筒温度场的确定方法及装置
CN105178930A (zh) 地面注汽管线内蒸汽热力参数计算方法及装置
CN112697313A (zh) 一种注蒸汽井筒传热过程中的热损失量评估方法及装置
CN105403590B (zh) 隔热管导热系数测试方法及装置
CN107130949A (zh) 深层稠油二氧化碳吞吐注入过程相态变化测试方法
CN106951666A (zh) 一种海洋天然气水合物层钻井井筒温度场计算方法
Kabir et al. Does gauge placement matter in downhole transient-data acquisition?
Aeschliman The Effect of Annulus Water on the Wellbore Heat Loss From a Steam Injection Well With Insulated Tubing
Abdelhafiz et al. Application of a thermal transient subsurface model to a coaxial borehole heat exchanger system
CN111520132B (zh) 一种确定地层中洞距离的方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant