CN107145705A - 一种获取循环温度场的方法及装置 - Google Patents
一种获取循环温度场的方法及装置 Download PDFInfo
- Publication number
- CN107145705A CN107145705A CN201710197241.3A CN201710197241A CN107145705A CN 107145705 A CN107145705 A CN 107145705A CN 201710197241 A CN201710197241 A CN 201710197241A CN 107145705 A CN107145705 A CN 107145705A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- msubsup
- msup
- 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
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Earth Drilling (AREA)
- Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
Abstract
一种获取循环温度场的方法及装置,包括:将井筒沿管柱径向划分为两个或两个以上组成部分;根据钻井状态参数,分别获得计算在钻进和不钻进过程中各组成部分的瞬态传热信息的传热微分方程;对传热微分方程进行离散和数值迭代处理,获得井筒的瞬态温度分布。本发明实施例根据钻井状态参数进行循环温度场的计算,实现了钻进和不钻进过程中循环温度场的获取,提升了获取循环温度场的工作效率。
Description
技术领域
本文涉及但不限于石油钻探技术,尤指一种获取循环温度场的方法及装置。
背景技术
井筒循环温度对钻井和固井工程的影响很大,它不仅关系到注水泥作业的成败和注水泥质量的高低,而且与井内压力平衡、井壁稳定、井内工作液体系选择、套管和钻柱强度设计等方面有关。因此,准确地确定井筒循环温度分布及其变化规律,对水泥浆体系设计、井控和安全快速钻进有重要的意义。
从20世纪60年代起,国外有很多学者针对井筒循环温度进行了研究,建立了不同的理论模型和算法,国内近二十多年来也对井下温度预测做了大量的研究,比较有代表性有针对地面钻井及注水泥过程建立的井内循环温度计算模型,井内循环温度计算模型包括以下几个部分:
管柱内流体:
管柱壁:
环空内液体:
地层:
计算公式(I)到(IV)中,q为排量,单位为立方每小时(m3/h);z为井深,单位为米(m);t为时间,单位为秒(s);rci为管柱内半径,单位为毫米(mm);rco为管柱外半径,单位为mm;rb为井眼半径,单位为mm;ρL为液体密度,单位为克每立方厘米(g/cm3);ρw为管柱材料的密度,单位为g/cm3;ρf为地层岩石的密度,单位为g/cm3;cL为液体比热,焦耳每千克(J/gK);cw为管柱材料的比热,单位为J/gK;cf为地层比热,单位为J/gK;kw为管柱材料的热导率,单位为瓦特/米开(W/mK);kf为地层岩石的热导率,单位为W/mK;Tc为管柱内流体的温度,单位为摄氏度(℃);Tw为管柱壁的温度,单位为℃;Ta为环空内液体的温度,单位为℃;Tf为地层温度,单位为℃;Tin为管柱入口的液体温度,单位为℃;Tout为环空出口液体温度,单位为℃;Ta为地表温度,单位为℃;G为地温梯度,单位为℃/m。hci、hco、hb分别为管柱内壁、管柱外壁和井壁的对流换热系数,单位为瓦特/平方米开(W/m2K);Qc、Qa分别为管柱内、环空内液体的热源,通常指液体流动摩擦生热。
采用上述计算公式,基于数值迭代方法,即可获得井筒循环温度;上述方法虽然简练,上述方法仅适用于地面钻井,不适用于深水钻井;且数值计算时钻进循环与不钻进循环考虑因素不一样,需要采用不同的数学方法分别进行计算;深水区域内通常为了增加泥浆携带量,都要在海底处向环空加注增压流量,因此整个环空区域内流体的流量是不均匀的,在上述计算公式中未考虑这一点,造成环空温度预测偏差较大。另外,深水区域内隔水管或管柱外侧为海水对流,而不是导热,上述计算公式也没有考虑海水对流的影响,计算获得深水钻井时的瞬态温度分布存在误差。
上述公开的通用计算方法,需要对模型离散后进行数值迭代,但钻进和不钻进过程所涉及的因素不同,计算方法也不同,目前没有公开的通用算法,无法完整地计算从钻进到不钻进过程中井筒的瞬态温度分布。
发明内容
以下是对本文详细描述的主题的概述。本概述并非是为了限制权利要求的保护范围。
本发明实施例提供一种获取循环温度场的方法及装置,能够实现完整的计算从钻进到不钻进过程中井筒的瞬态温度分布。
本发明实施例提供了一种获取循环温度场的方法,包括:
将井筒沿管柱径向划分为两个或两个以上组成部分;
根据钻井状态参数,分别获得计算在钻进和不钻进过程中各组成部分的瞬态传热信息的传热微分方程;
对传热微分方程进行离散和数值迭代处理,获得井筒的瞬态温度分布。
可选的,所述将井筒沿管柱径向划分为两个或两个以上组成部分包括:
将井筒沿管柱径向划分为管柱内流体、筒壁、环空流体、隔水管;
其中,划分的各所述组成部分分别包括对应的预设数个节点。
可选的,获得所述传热微分方程之前,所述方法还包括:
获取所述钻井状态参数;
其中,所述钻井状态参数包括:时间步长Δt、空间步长Δz、钻井液入口流量G、入口温度Tin、各时间节点对应的各空间节点位置、各时间节点的钻速、井身结构参数、物性参数、钻头机械磨损参数和/或初始时刻各空间节点温度值。
可选的,所述获取钻井状态参数包括:
获取井身结构参数包括:rco、rci、rgi、rgo、rwe;
物性参数包括:ρL、kL、cp、kw、ρw、cw、kg、ρg、cg、kf、ρf、μf;
钻头机械磨损参数包括:钻压Fdb、钻头平均直径Ddb;
水深Ls、循环时间t_end、地层井深Ld;
计算初始总井深LT为:LT=Ls+Ld;
空间步长:Δz=LT/(n-1);
时间步长:Δt=t_end/(m-1);
各时间节点对应的各空间节点位置:Lz(j)=(j-1)Δz,j=1~n;
各时间节点钻速:ud(i),i=1~m;
初始时刻井筒各节点温度:
Lz(j)≤Ls时,
Lz(j)>Ls时,(j=1~n;i=1);
井身结构参数包括:rco、rci、rgi、rgo、rwe;
物性参数包括:ρL、kL、cp、kw、ρw、cw、kg、ρg、cg、kf、ρf、μf;
钻头机械磨损参数包括:钻压Fdb、钻头平均直径Ddb;
其中,ud为钻进速度,钻进时大于零,不钻进时等于0,单位米每秒m/s;循环时间t_end单位为秒s;时间节点数为m;井筒节点数为n;海水沿井深的温度分布为Tf,单位为℃;井壁沿深度方向的初始温度分布为单位为℃;为第i时刻,管柱内流体在轴向方向上第j节点的温度,单位为℃;为第i时刻,筒壁在轴向方向上第j节点的温度,单位为℃;为第i时刻,环空流体在轴向方向上第j节点的温度,单位为℃;为第i时刻,隔水管在轴向方向上第j节点的温度,单位为℃;为第i时刻,井壁在轴向方向上第j节点的温度,单位为℃;rco为筒壁外壁面半径,单位为米m;rci为筒壁内壁面半径,单位为m;ρL为钻井液密度,单位为千克每立方米kg/m3;kL为钻井液导热系数,单位为焦耳每千克J/kgK;cp为钻井液比热,单位为焦耳每千克J/kgK;kw为筒壁的导热系数,单位为瓦每米开W/mk;ρw为筒壁密度,单位为kg/m3;cw为筒壁比热,单位为J/kgK;rwe为井壁内径,单位为m;ρg为隔水管密度,单位为kg/m3;kg为隔水管导热系数,单位为W/mK;cg为隔水管比热,单位为J/kgK;kf为海水导热系数,单位为W/mk;ρf为海水密度,单位为kg/m3;μf为海水粘度系数,单位为千克每米每秒kg/ms;rgi为隔水管内径,单位为米m;rgo为隔水管外径,单位为米m。
可选的,所述方法还包括:
进行所述数值迭代处理时,将各节点在前一个时间步长的温度值作为计算当前时间步长内温度的初始值;
对总井深、空间节点数量根据下述公式进行更新:
更新的总井深:LT=LT+ud Δt
更新的空间节点数量:
若LT-Lz(n)≥Δz,则n=n+1。
可选的,所述管柱内流体的传热微分方程为:
或
式(1)或式(2)中,Qhp为管柱内流体的热源,单位为瓦每米W/m;G为钻井液体积流量,单位为立方米每秒m3/s;hci为筒壁内壁面的对流换热系数,单位为瓦每平方米开W/m2k;Tp为筒内液体的温度,单位为摄氏度℃;Tw为筒壁的温度,单位为℃;z为轴向长度,单位为米m;t为时间,单位为秒s。
可选的,所述对传热微分方程进行离散处理包括:
对管柱内流体的传热微分方程进行离散处理,获得第i+1时刻,包含井身轴向方向上第j节点的管柱内流体温度的计算方程:
其中:
可选的,钻进过程中,所述管柱内流体的热源Qhp为:
Qhp=GPp+Qdb/△z+G△pb/△z (4)
式(4)中,GPp为各节点在钻进过程中均产生的发热量;Qdb/△z、G△pb/△z为钻头所对应的节点的发热量;Pp为单位长度的流动磨损,单位为帕每米Pa/m;Qdb为钻进过程中,由钻头机械摩擦做功转化的热损耗,Qdb=fdbFdburl,单位为瓦W,其中,url=πDdbω钻头旋转线速度,单位为米每秒m/s;fdb为摩擦系数,取值为0~1;Fdb为钻压,单位为牛N;Ddb为钻头平均直径,单位为m;ω为转速,单位为弧度每秒rad/s;△pb为泥浆经过钻头喷嘴的节流作用产生的局部压力损耗,△pb单位为帕Pa;C代表喷嘴流量系数,无因次,取值范围为0.914~0.98;Ab代表钻头水口总面积,单位为平方米m2;△z为钻头节点与轴向相邻节点间距离,单位为米m;
不钻进过程中,管柱内流体的热源Qhp为:
Qhp=GPp+G△pb/△z (5)。
可选的,钻进过程中,所述管柱内流体在旋转的管柱内螺旋流动,对应的对流换热系数hci为:
式(6)中,Nup为管柱内流体的努谢尔数,无量纲;等效流速单位为米每秒m/s,Reeff为等效流速所对应的等效雷诺数;up为管柱内流体的轴向流动速度,单位为m/s;Pr为流体的普朗特数,无量纲;α为旋转流动对换热影响的权重系数,取值范围约为0.25~1;系数Ah取值范围为0.01~0.03;系数γ取值范围为0~0.5;
不钻进过程中,对流换热系数hci通过非牛顿流体的计算公式获得。
可选的,所述筒壁的传热微分方程为:
或,
式(7)或式(8)中,hco为筒壁外壁面的对流换热系数,单位为W/m2k;Ta为环空流体的温度,单位为℃;Tw为筒壁的温度,单位为℃。
可选的,地层区域内,不钻进过程中,筒壁外侧的所述对流换热系数hco通过非牛顿流体的计算公式获得;
有隔水管的深水区域,不钻进过程中,筒壁外侧的所述对流换热系数hco通过非牛顿流体的计算公式获得;
地层区域内,钻进过程中,筒壁外侧对流换热系数hco为:
式(9)中,Nua为环空流体的努谢尔数,无量纲;ueff为等效流速,ueff单位为m/s,Reeff为根据等效流速ueff确定的等效雷诺数;ua为环空流体轴向流动速度,单位为m/s;α为旋转流动对换热影响的权重系数,取值范围约为0.25~1;系数Ah取值范围为0.01~0.03;系数γ取值范围为0~0.5;
有隔水管的深水区域,钻进过程中,筒壁外侧对流换热系数hco为:
无隔水管的深水区域,环空流体的温度Ta为海水的温度Tf;对应管柱壁外侧的对流换热系数hco为海水的流动换热系数,hco为:其中,海水努谢尔数Ref为海水粘度系数uf对应的雷诺数;Prf为海水普朗特数,无量纲;系数c的取值范围是0.024~0.88;n的取值范围是0.33~0.805。
可选的,所述对传热微分方程进行离散处理包括:
对筒壁的传热微分方程进行离散处理,获得第i+1时刻,包含井身轴向方向上第j节点的筒壁温度的计算方程:
其中:
可选的,所述环空流体为地层区域环空流体时,所述地层区域的环空流体的传热微分方程为:
或,
式(12)或(13)中,Qha为液体流动摩擦生热产生的环空流体的热源,单位为瓦每米W/m,Qha=GPa,Pa为单位长度流动磨损,单位为Pa/m,Pa通过非牛顿流体的流动摩阻计算公式获得;G为钻井液体积流量,单位为m3/s;hwe为井壁的对流换热系数,单位为W/mk;Twe为井壁温度,单位为℃。
可选的,所述对传热微分方程进行离散处理包括:
对地层区域的环空流体的传热微分方程进行离散处理,获得第i+1时刻,包含井身轴向方向上第j节点的地层区域的环空流体温度的计算方程:
其中:
可选的,地层区域内,不钻进过程中,筒壁外侧及井壁内侧的对流换热系数hco和hwe通过非牛顿流体的计算公式获得;
地层区域内,钻进过程中,筒壁外侧的对流换热系数hco:
井壁内侧的对流换热系数hwe与筒壁外侧对流换热系数hco相同;
井壁温度Twe由一维稳态传热模型或二维地层导热模型获得。
可选的,有隔水管的深水区域,所述环空流体的传热微分方程为:
或
式中,hgi为隔水管内壁面对流换热系数,单位为W/m2K;Tg为隔水管壁面温度,单位为℃;Gs为深水区域环空流体的体积流量:Gs=G+Gsa,其中,Gsa为海底增压流量,单位为m3/s;对有隔水管的深水区域,钻进过程中,hco=hgi;
对有隔水管的深水区域,钻进过程中,hco=hgi;
不钻进过程中,筒壁外侧及隔水管内侧的对流换热系数hco和hgi通过非牛顿流体的计算公式获得。
可选的,所述隔水层的传热微分方程为:
或,
式中,Tg为隔水管的温度,单位为℃;Tf为海水温度,单位为℃;hgo为隔水管外壁面的对流换热系数,单位为W/m2k:其中,海水努谢尔数Ref为海水粘度系数uf对应的雷诺数
可选的,所述对传热微分方程进行离散处理包括:
对隔水管的传热微分方程进行离散处理,获得第i+1时刻,包含井身轴向方向上第j节点的隔水管温度的计算方程:
其中:
另一方面,本发明实施例还提供一种获取循环温度场的装置,包括:划分单元、传热微分单元和离散迭代单元;其中,
划分单元用于,将井筒沿管柱径向划分为两个或两个以上组成部分;
传热微分单元用于,根据钻井状态参数,分别获得计算在钻进和不钻进过程中各组成部分的瞬态传热信息的传热微分方程;
离散迭代单元用于,对传热微分方程进行离散和数值迭代处理,获得井筒的瞬态温度分布。
可选的,所述划分单元具体用于:
将井筒沿管柱径向划分为管柱内流体、筒壁、环空流体、隔水管;
其中,划分的各所述组成部分分别包括对应的预设数个节点。
可选的,所述装置还包括获取单元,用于获得所述传热微分方程之前,获取所述钻井状态参数;
其中,所述钻井状态参数包括:时间步长Δt、空间步长Δz、钻井液入口流量G、入口温度Tin、各时间节点对应的各空间节点位置、各时间节点的钻速、井身结构参数、物性参数、钻头机械磨损参数和/或初始时刻各空间节点温度值。
与相关技术相比,本申请技术方案包括:将井筒沿管柱径向划分为两个或两个以上组成部分;根据钻井状态参数,分别获得计算在钻进和不钻进过程中各组成部分的瞬态传热信息的传热微分方程;对传热微分方程进行离散和数值迭代处理,获得井筒的瞬态温度分布。本发明实施例根据钻井状态参数进行循环温度场的计算,实现了钻进和不钻进过程中循环温度场的获取,提升了获取循环温度场的工作效率。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本发明技术方案的进一步理解,并且构成说明书的一部分,与本申请的实施例一起用于解释本发明的技术方案,并不构成对本发明技术方案的限制。
图1为本发明实施例获取循环温度场的方法的流程图;
图2为本发明实施例划分的组成部分的结构示意图;
图3为本发明实施例获取循环温度场的装置的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下文中将结合附图对本发明的实施例进行详细说明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互任意组合。
在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行。并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
图1为本发明实施例获取循环温度场的方法的流程图,如图1所示,包括:
步骤100、将井筒沿管柱径向划分为两个或两个以上组成部分;
可选的,本发明实施例将井筒沿管柱径向划分为两个或两个以上组成部分包括:
将井筒沿管柱径向划分为管柱内流体、筒壁、环空流体、隔水管;
其中,划分的各所述组成部分分别包括对应的预设数个节点。
需要说明的是,本发明实施例划分的每一个组成部分可以分别通过对应的两个或两个以上节点组成。节点个数根据包括井深、Δz在内的参数确定,各个组成部分包含的节点个数可以相同,划分节点的方法可以参照相关技术,在此不做赘述。图2为本发明实施例划分的组成部分的结构示意图,如图2所示,本发明实施例将井筒沿管柱径向划分为管柱内流体、筒壁、环空流体、隔水管四个组成部分。
步骤101、根据钻井状态参数,分别获得计算在钻进和不钻进过程中各组成部分的瞬态传热信息的传热微分方程;
可选的,步骤101获得传热微分方程之前,本发明实施例方法还包括:
获取所述钻井状态参数;
其中,钻井状态参数包括:时间步长Δt、空间步长Δz、钻井液入口流量G、入口温度Tin、各时间节点对应的各空间节点位置、各时间节点的钻速、井身结构参数、物性参数、钻头机械磨损参数和/或初始时刻各空间节点温度值。
获取钻井状态参数包括:
获取井身结构参数包括:rco、rci、rgi、rgo、rwe;
物性参数包括:ρL、kL、cp、kw、ρw、cw、kg、ρg、cg、kf、ρf、μf;
钻头机械磨损参数包括:钻压Fdb、钻头平均直径Ddb;
水深Ls、循环时间t_end、地层井深Ld;
计算初始总井深LT为:LT=Ls+Ld;
空间步长:Δz=LT/(n-1);
时间步长:Δt=t_end/(m-1);
各时间节点对应的各空间节点位置:Lz(j)=(j-1)Δz,j=1~n;
各时间节点钻速:ud(i),i=1~m;
初始时刻井筒各节点温度:
Lz(j)≤Ls时,
Lz(j)>Ls时,(j=1~n;i=1);
井身结构参数包括:rco、rci、rgi、rgo、rwe;
物性参数包括:ρL、kL、cp、kw、ρw、cw、kg、ρg、cg、kf、ρf、μf;
钻头机械磨损参数包括:钻压Fdb、钻头平均直径Ddb;
其中,ud为钻进速度,钻进时大于零,不钻进时等于0,单位米每秒m/s;循环时间t_end单位为秒s;时间节点数为m;井筒节点数为n;海水沿井深的温度分布为Tf,单位为℃;井壁沿深度方向的初始温度分布为单位为℃;为第i时刻,管柱内流体在轴向方向上第j节点的温度,单位为℃;为第i时刻,筒壁在轴向方向上第j节点的温度,单位为℃;为第i时刻,环空流体在轴向方向上第j节点的温度,单位为℃;为第i时刻,隔水管在轴向方向上第j节点的温度,单位为℃;为第i时刻,井壁在轴向方向上第j节点的温度,单位为℃;rco为筒壁外壁面半径,单位为米m;rci为筒壁内壁面半径,单位为m;ρL为钻井液密度,单位为千克每立方米kg/m3;kL为钻井液导热系数,单位为焦耳每千克J/kgK;cp为钻井液比热,单位为焦耳每千克J/kgK;kw为筒壁的导热系数,单位为瓦每米开W/mk;ρw为筒壁密度,单位为kg/m3;cw为筒壁比热,单位为J/kgK;rwe为井壁内径,单位为m;ρg为隔水管密度,单位为kg/m3;kg为隔水管导热系数,单位为W/mK;cg为隔水管比热,单位为J/kgK;kf为海水导热系数,单位为W/mk;ρf为海水密度,单位为kg/m3;μf为海水粘度系数,单位为千克每米每秒kg/ms;rgi为隔水管内径,单位为米m;rgo为隔水管外径,单位为米m。
需要说明的是,本发明实施例根据以下公式计算:
更新的总井深:LT=LT+udΔt;
更新的空间节点数量:若LT-Lz(n)≥Δz,则n=n+1;
更新的空间节点位置为:Lz(j)=(j-1)Δz,j=1~n。
可选的,本发明实施例方法还包括:
进行所述数值迭代处理时,将各节点在前一个时间步长的温度值作为计算当前时间步长内温度的初始值。
需要说明的是,本发明实施例无论地面钻井(Ls≤0)或深水钻井(Ls>0),计算方法都适用。井筒各节点临时温度赋值(j=1~n;i=2~m)可以为:
Lz(j)≤Ls时
Lz(j)>Ls时
可选的,本发明实施例管柱内流体的传热微分方程为:
或
式(1)或式(2)中,Qhp为管柱内流体的热源,单位为瓦每米W/m;G为钻井液体积流量,单位为立方米每秒m3/s;hci为筒壁内壁面的对流换热系数,单位为瓦每平方米开W/m2k;Tp为筒内液体的温度,单位为摄氏度℃;Tw为筒壁的温度,单位为℃;z为轴向长度,单位为米m;t为时间,单位为秒s。
可选的,本发明实施例:
钻进过程中,所述管柱内流体的热源Qhp为:
Qhp=GPp+Qdb/△z+G△pb/△z(4)
式(4)中,GPp为各节点在钻进过程中均产生的发热量;Qdb/△z、G△pb/△z为钻头所对应的节点的发热量;Pp为单位长度的流动磨损,单位为帕每米Pa/m;Qdb为钻进过程中,由钻头机械摩擦做功转化的热损耗,Qdb=fdbFdburl,单位为瓦W,其中,url=πDdbω钻头旋转线速度,单位为米每秒m/s;fdb为摩擦系数,取值为0~1;Fdb为钻压,单位为牛N;Ddb为钻头平均直径,单位为m;ω为转速,单位为弧度每秒rad/s;△pb为泥浆经过钻头喷嘴的节流作用产生的局部压力损耗,△pb单位为帕Pa;C代表喷嘴流量系数,无因次,取值范围为0.914~0.98;Ab代表钻头水口总面积,单位为平方米m2;△z为钻头节点与轴向相邻节点间距离,单位为米m;
不钻进过程中,管柱内流体的热源Qhp为:
Qhp=GPp+G△pb/△z (5)。
可选的,本发明实施例钻进过程中,所述管柱内流体在旋转的管柱内螺旋流动,对应的对流换热系数hci为:
式(6)中,Nup为管柱内流体的努谢尔数,无量纲;等效流速单位为米每秒m/s,Reeff为等效流速所对应的等效雷诺数;up为管柱内流体的轴向流动速度,单位为m/s;Pr为流体的普朗特数,无量纲;α为旋转流动对换热影响的权重系数,取值范围约为0.25~1;系数Ah取值范围为0.01~0.03;系数γ取值范围为0~0.5;
不钻进过程中,对流换热系数hci通过非牛顿流体的计算公式获得。
需要说明的是,本发明实施例非牛顿流体的计算公式包括由豪森(Hausen)等推荐的非牛顿流体的经验公式。
可选的,本发明实施例筒壁的传热微分方程为:
或,
式(7)或式(8)中,hco为筒壁外壁面的对流换热系数,单位为W/m2k;Ta为环空流体的温度,单位为℃;Tw为筒壁的温度,单位为℃。
可选的,本发明实施例,地层区域内,不钻进过程中,筒壁外侧的所述对流换热系数hco通过非牛顿流体的计算公式获得;
有隔水管的深水区域,不钻进过程中,筒壁外侧的所述对流换热系数hco通过非牛顿流体的计算公式获得;
地层区域内,钻进过程中,筒壁外侧对流换热系数hco为:
式(9)中,Nua为环空流体的努谢尔数,无量纲;ueff为等效流速,ueff单位为m/s,Reeff为根据等效流速ueff确定的等效雷诺数;ua为环空流体轴向流动速度,单位为m/s;α为旋转流动对换热影响的权重系数,取值范围约为0.25~1;系数Ah取值范围为0.01~0.03;系数γ取值范围为0~0.5;
有隔水管的深水区域,钻进过程中,筒壁外侧对流换热系数hco为:
无隔水管的深水区域,环空流体的温度Ta为海水的温度Tf;对应管柱壁外侧的对流换热系数hco为海水的流动换热系数,hco为:其中,海水努谢尔数Ref为海水粘度系数uf对应的雷诺数;Prf为海水普朗特数,无量纲;系数c的取值范围是0.024~0.88;n的取值范围是0.33~0.805。
可选的,所述环空流体为地层区域环空流体时,所述地层区域的环空流体的传热微分方程为:
或,
式(12)或(13)中,Qha为液体流动摩擦生热产生的环空流体的热源,单位为瓦每米W/m,Qha=GPa,Pa为单位长度流动磨损,单位为Pa/m,Pa通过非牛顿流体的流动摩阻计算公式获得;hwe为井壁的对流换热系数,单位为W/mk;Twe为井壁温度,单位为℃。
可选的,本发明实施例中:
地层区域内,不钻进过程中,筒壁外侧及井壁内侧的对流换热系数hco和hwe通过非牛顿流体的计算公式获得;
地层区域内,钻进过程中,筒壁外侧的对流换热系数hco:
井壁内侧的对流换热系数hwe与筒壁外侧对流换热系数hco相同;
井壁温度Twe由一维稳态传热模型或二维地层导热模型获得。
需要说明的时,本发明实施例井壁温度可由Hansan-Kabir(Hansan-Kabir为现有算法,不做赘述)一维稳态传热模型计算,也可通过二维地层导热模型计算得到。本发明实施例公式中未做说明的参数为本领域技术人员公知,在此不做赘述。
可选的,有隔水管的深水区域,所述环空流体的传热微分方程为:
或
式中hgi为隔水管内壁面对流换热系数,单位为W/m2K;Tg为隔水管壁面温度,单位为℃;Gs为深水区域环空流体的体积流量:Gs=G+Gsa,其中,Gsa为海底增压流量,单位为m3/s;
对有隔水管的深水区域,钻进过程中,hco=hgi;
不钻进过程中,筒壁外侧及隔水管内侧的对流换热系数hco和hgi通过非牛顿流体的计算公式获得。
可选的,本发明实施例所述隔水管的传热微分方程为:
或,
式中,Tg为隔水管的温度,单位为℃;Tf为海水温度,单位为℃;hgo为隔水管外壁面的对流换热系数,单位为W/m2k:其中,海水努谢尔数Ref为海水粘度系数uf对应的雷诺数
步骤102、对传热微分方程进行离散和数值迭代处理,获得井筒的瞬态温度分布。
可选的,对传热微分方程进行离散处理包括:
对管柱内流体的传热微分方程(1)进行离散处理,获得第i+1时刻,包含井身轴向方向上第j节点的管柱内流体温度的计算方程:
其中:
可选的,本发明实施例对传热微分方程进行离散处理包括:
对筒壁的传热微分方程(7)进行离散处理,获得第i+1时刻,包含井身轴向方向上第j节点的筒壁温度的计算方程:
其中:
可选的,本发明实施例对传热微分方程进行离散处理包括:
对地层区域的环空流体的传热微分方程进行离散处理,获得第i+1时刻,包含井身轴向方向上第j节点的地层区域的环空流体温度的计算方程:
其中:
可选的,本发明实施例对传热微分方程进行离散处理包括:
对隔水管的传热微分方程(17)进行离散处理,获得第i+1时刻,包含井身轴向方向上第j节点的隔水管温度的计算方程:
其中:
需要说明的是,本发明实施例迭代处理过程和相关技术中迭代处理过程原理相同,在此不做赘述。
本发明实施例根据钻井状态参数进行循环温度场的计算,实现了钻进和不钻进过程中循环温度场的获取,提升了获取循环温度场的工作效率。
图3为本发明实施例获取循环温度场的装置的结构框图,如图3所示,包括:划分单元、传热微分单元和离散迭代单元;其中,
划分单元用于,将井筒沿管柱径向划分为两个或两个以上组成部分;
可选的,本发明实施例划分单元具体用于:
将井筒沿管柱径向划分为管柱内流体、筒壁、环空流体、隔水管;
其中,划分的各所述组成部分分别包括对应的预设数值个微元体。
传热微分单元用于,根据钻井状态参数,分别获得计算在钻进和不钻进过程中各组成部分的瞬态传热信息的传热微分方程;
可选的,本发明实施例装置还包括获取单元,用于获得所述传热微分方程之前,获取所述钻井状态参数;
其中,所述钻井状态参数包括:时间步长Δt、空间步长Δz、钻井液入口流量G、入口温度Tin、各时间节点对应的各空间节点位置、各时间节点的钻速、井身结构参数、物性参数、钻头机械磨损参数和/或初始时刻各空间节点温度值。
离散迭代单元用于,对传热微分方程进行离散和数值迭代处理,获得井筒的瞬态温度分布。
本发明实施例传热微分单元获得的计算在钻进和不钻进过程中各组成部分的瞬态传热信息的传热微分方程所需的钻井状态参数包括:时间步长Δt、空间步长Δz、钻井液入口流量G、入口温度Tin、各时间节点对应的各空间节点位置、各时间节点的钻速、井身结构参数、物性参数、钻头机械磨损参数和/或初始时刻各空间节点温度值。
钻井状态参数可以通过以下方式获取:
从系统属性及参数信息中获取井身结构参数包括:rco、rci、rgi、rgo、rwe;
物性参数包括:ρL、kL、cp、kw、ρw、cw、kg、ρg、cg、kf、ρf、μf;
钻头机械磨损参数包括:钻压Fdb、钻头平均直径Ddb;
水深Ls、循环时间t_end、地层井深Ld;
计算初始总井深LT为:LT=Ls+Ld;
空间步长:Δz=LT/(n-1);
时间步长:Δt=t_end/(m-1);
各时间节点对应的各空间节点位置:Lz(j)=(j-1)Δz,j=1~n;
各时间节点钻速:ud(i),i=1~m;
初始时刻井筒各节点温度:
Lz(j)≤Ls时,
Lz(j)>Ls时,(j=1~n;i=1);
井身结构参数包括:rco、rci、rgi、rgo、rwe;
物性参数包括:ρL、kL、cp、kw、ρw、cw、kg、ρg、cg、kf、ρf、μf;
钻头机械磨损参数包括:钻压Fdb、钻头平均直径Ddb;
对总井深、空间节点数量、空间节点位置根据下述公式进行更新:
更新的总井深:LT=LT+udΔt;
更新的空间节点数量:
若LT-Lz(n)≥Δz,则n=n+1;
各空间节点位置:Lz(j)=(j-1)Δz,j=1~n;。
本发明实施例,进行所述数值迭代处理时,将各节点在前一个时间步长的温度值作为计算当前时间步长内温度的初始值。
本发明实施例管内流体入口温度边界为:
筒壁及环空井底节点温度边界条件:
本领域普通技术人员可以理解上述方法中的全部或部分步骤可通过程序来指令相关硬件(例如处理器)完成,所述程序可以存储于计算机可读存储介质中,如只读存储器、磁盘或光盘等。可选地,上述实施例的全部或部分步骤也可以使用一个或多个集成电路来实现。相应地,上述实施例中的每个模块/单元可以采用硬件的形式实现,例如通过集成电路来实现其相应功能,也可以采用软件功能模块的形式实现,例如通过处理器执行存储于存储器中的程序/指令来实现其相应功能。本发明不限制于任何特定形式的硬件和软件的结合。
虽然本发明所揭露的实施方式如上,但所述的内容仅为便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属领域内的技术人员,在不脱离本发明所揭露的精神和范围的前提下,可以在实施的形式及细节上进行任何的修改与变化,但本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (21)
1.一种获取循环温度场的方法,其特征在于,包括:
将井筒沿管柱径向划分为两个或两个以上组成部分;
根据钻井状态参数,分别获得计算在钻进和不钻进过程中各组成部分的瞬态传热信息的传热微分方程;
对传热微分方程进行离散和数值迭代处理,获得井筒的瞬态温度分布。
2.根据权利要求1所述的方法,其特征在于,所述将井筒沿管柱径向划分为两个或两个以上组成部分包括:
将井筒沿管柱径向划分为管柱内流体、筒壁、环空流体、隔水管;
其中,划分的各所述组成部分分别包括对应的预设数个节点。
3.根据权利要求1或2所述的方法,其特征在于,获得所述传热微分方程之前,所述方法还包括:
获取所述钻井状态参数;
其中,所述钻井状态参数包括:时间步长△t、空间步长△z、钻井液入口流量G、入口温度Tin、各时间节点对应的各空间节点位置、各时间节点的钻速、井身结构参数、物性参数、钻头机械磨损参数和/或初始时刻各空间节点温度值。
4.根据权利要求3所述的方法,其特征在于,所述获取钻井状态参数包括:
获取井身结构参数包括:rco、rci、rgi、rgo、rwe;
物性参数包括:ρL、kL、cp、kw、ρw、cw、kg、ρg、cg、kf、ρf、μf;
钻头机械磨损参数包括:钻压Fdb、钻头平均直径Ddb;
水深Ls、循环时间t_end、地层井深Ld;
计算初始总井深LT为:LT=Ls+Ld;
空间步长:△z=LT/(n-1);
时间步长:△t=t_end/(m-1);
各时间节点对应的各空间节点位置:Lz(j)=(j-1)△z,j=1~n;
各时间节点钻速:ud(i),i=1~m;
初始时刻井筒各节点温度:
Lz(j)≤Ls时,Tp,j i=Tw,j i=Ta,j i=Tg,j i=Tf;(j=1~n;i=1);
Lz(j)>Ls时,Tp,j i=Tw,j i=Ta,j i=Twe,,j i=Twe 0;(j=1~n;i=1);
井身结构参数包括:rco、rci、rgi、rgo、rwe;
物性参数包括:ρL、kL、cp、kw、ρw、cw、kg、ρg、cg、kf、ρf、μf;
钻头机械磨损参数包括:钻压Fdb、钻头平均直径Ddb;
其中,ud为钻进速度,钻进时大于零,不钻进时等于0,单位米每秒m/s;循环时间t_end单位为秒s;时间节点数为m;井筒节点数为n;海水沿井深的温度分布为Tf,单位为℃;井壁沿深度方向的初始温度分布为Twe 0,单位为℃;为第i时刻,管柱内流体在轴向方向上第j节点的温度,单位为℃;Tw,j i为第i时刻,筒壁在轴向方向上第j节点的温度,单位为℃;Ta,j i为第i时刻,环空流体在轴向方向上第j节点的温度,单位为℃;Tg,j i为第i时刻,隔水管在轴向方向上第j节点的温度,单位为℃;Twe,,j i为第i时刻,井壁在轴向方向上第j节点的温度,单位为℃;rco为筒壁外壁面半径,单位为米m;rci为筒壁内壁面半径,单位为m;ρL为钻井液密度,单位为千克每立方米kg/m3;kL为钻井液导热系数,单位为焦耳每千克J/kgK;cp为钻井液比热,单位为焦耳每千克J/kgK;kw为筒壁的导热系数,单位为瓦每米开W/mk;ρw为筒壁密度,单位为kg/m3;cw为筒壁比热,单位为J/kgK;rwe为井壁内径,单位为m;ρg为隔水管密度,单位为kg/m3;kg为隔水管导热系数,单位为W/mK;cg为隔水管比热,单位为J/kgK;kf为海水导热系数,单位为W/mk;ρf为海水密度,单位为kg/m3;μf为海水粘度系数,单位为千克每米每秒kg/ms;rgi为隔水管内径,单位为米m;rgo为隔水管外径,单位为米m。
5.根据权利要求4所述的方法,其特征在于,所述方法还包括:
进行所述数值迭代处理时,将各节点在前一个时间步长的温度值作为计算当前时间步长内温度的初始值;
对总井深、空间节点数量根据下述公式进行更新:
更新的总井深:LT=LT+ud △t;
更新的空间节点数量:
若LT-Lz(n)≥△z,则n=n+1。
6.根据权利要求2所述的方法,其特征在于,所述管柱内流体的传热微分方程为:
或
<mrow>
<msub>
<mi>Q</mi>
<mi>hp</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>Gc</mi>
<mi>p</mi>
</msub>
<mfrac>
<msub>
<mrow>
<mo>&PartialD;</mo>
<mi>T</mi>
</mrow>
<mi>p</mi>
</msub>
<mrow>
<mo>&PartialD;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mn>2</mn>
<mi>&pi;</mi>
<msub>
<mi>r</mi>
<mi>ci</mi>
</msub>
<msub>
<mi>h</mi>
<mi>ci</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
<msubsup>
<mi>&pi;r</mi>
<mi>ci</mi>
<mn>2</mn>
</msubsup>
<mfrac>
<mrow>
<msup>
<mo>&PartialD;</mo>
<mn>2</mn>
</msup>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
</mrow>
<msup>
<mrow>
<mo>&PartialD;</mo>
<mi>z</mi>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mo>=</mo>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<msup>
<msub>
<mi>&pi;r</mi>
<mi>ci</mi>
</msub>
<mn>2</mn>
</msup>
<mfrac>
<msub>
<mrow>
<mo>&PartialD;</mo>
<mi>T</mi>
</mrow>
<mi>p</mi>
</msub>
<mrow>
<mo>&PartialD;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
式(1)或式(2)中,Qhp为管柱内流体的热源,单位为瓦每米W/m;G为钻井液体积流量,单位为立方米每秒m3/s;hci为筒壁内壁面的对流换热系数,单位为瓦每平方米开W/m2k;Tp为筒内液体的温度,单位为摄氏度℃;Tw为筒壁的温度,单位为℃;z为轴向长度,单位为米m;t为时间,单位为秒s。
7.根据权利要求6所述的方法,其特征在于,所述对传热微分方程进行离散处理包括:
对管柱内流体的传热微分方程进行离散处理,获得第i+1时刻,包含井身轴向方向上第j节点的管柱内流体温度的计算方程:
<mrow>
<msub>
<mi>A</mi>
<mi>p</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>h</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>B</mi>
<mi>p</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>E</mi>
<mi>p</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>p</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:
<mrow>
<msub>
<mi>A</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>Gc</mi>
<mi>p</mi>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>z</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>l</mi>
<mi>n</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;k</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
<msubsup>
<mi>&pi;r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<msup>
<mi>&Delta;z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<msup>
<msub>
<mi>&pi;r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>B</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>Gc</mi>
<mi>p</mi>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>z</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
<msubsup>
<mi>&pi;r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<msup>
<mi>&Delta;z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>l</mi>
<mi>n</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;k</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
2
<mrow>
<msub>
<mi>D</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<msup>
<msub>
<mi>&pi;r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
<msup>
<msub>
<mi>&pi;r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<mi>&Delta;z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>.</mo>
</mrow>
8.根据权利要求6或7所述的方法,其特征在于,
钻进过程中,所述管柱内流体的热源Qhp为:
Qhp=GPp+Qdb/△z+G△pb/△z (4)
式(4)中,GPp为各节点在钻进过程中均产生的发热量;Qdb/△z、G△pb/△z为钻头所对应的节点的发热量;Pp为单位长度的流动磨损,单位为帕每米Pa/m;Qdb为钻进过程中,由钻头机械摩擦做功转化的热损耗,Qdb=fdbFdburl,单位为瓦W,其中,url=πDdbω钻头旋转线速度,单位为米每秒m/s;fdb为摩擦系数,取值为0~1;Fdb为钻压,单位为牛N;Ddb为钻头平均直径,单位为m;ω为转速,单位为弧度每秒rad/s;△pb为泥浆经过钻头喷嘴的节流作用产生的局部压力损耗,单位为帕Pa;C代表喷嘴流量系数,无因次,取值范围为0.914~0.98;Ab代表钻头水口总面积,单位为平方米m2;△z为钻头节点与轴向相邻节点间距离,单位为米m;
不钻进过程中,管柱内流体的热源Qhp为:
Qhp=GPp+G△pb/△z (5)。
9.根据权利要求6或7所述的方法,其特征在于,
钻进过程中,所述管柱内流体在旋转的管柱内螺旋流动,对应的对流换热系数hci为:
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>Nu</mi>
<mi>p</mi>
</msub>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mi>h</mi>
</msub>
<msubsup>
<mi>Re</mi>
<mrow>
<mi>e</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<msup>
<mi>Pr</mi>
<mi>&gamma;</mi>
</msup>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式(6)中,Nup为管柱内流体的努谢尔数,无量纲;等效流速单位为米每秒m/s,Reeff为等效流速所对应的等效雷诺数;up为管柱内流体的轴向流动速度,单位为m/s;Pr为流体的普朗特数,无量纲;α为旋转流动对换热影响的权重系数,取值范围约为0.25~1;系数Ah取值范围为0.01~0.03;系数γ取值范围为0~0.5;
不钻进过程中,对流换热系数hci通过非牛顿流体的计算公式获得。
10.根据权利要求2所述的方法,其特征在于,所述筒壁的传热微分方程为:
或,
<mrow>
<msub>
<mi>k</mi>
<mi>w</mi>
</msub>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>&rho;</mi>
<mi>w</mi>
</msub>
<msub>
<mi>c</mi>
<mi>w</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式(7)或式(8)中,hco为筒壁外壁面的对流换热系数,单位为W/m2k;Ta为环空流体的温度,单位为℃;Tw为筒壁的温度,单位为℃。
11.根据权利要求10所述的方法,其特征在于,
地层区域内,不钻进过程中,筒壁外侧的所述对流换热系数hco通过非牛顿流体的计算公式获得;
有隔水管的深水区域,不钻进过程中,筒壁外侧的所述对流换热系数hco通过非牛顿流体的计算公式获得;
地层区域内,钻进过程中,筒壁外侧对流换热系数hco为:
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>Nu</mi>
<mi>a</mi>
</msub>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mi>h</mi>
</msub>
<msubsup>
<mi>Re</mi>
<mrow>
<mi>e</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<msup>
<mi>Pr</mi>
<mi>&gamma;</mi>
</msup>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
式(9)中,Nua为环空流体的努谢尔数,无量纲;ueff为等效流速,ueff单位为m/s,Reeff为根据等效流速ueff确定的等效雷诺数;ua为环空流体轴向流动速度,单位为m/s;α为旋转流动对换热影响的权重系数,取值范围约为0.25~1;系数Ah取值范围为0.01~0.03;系数γ取值范围为0~0.5;
有隔水管的深水区域,钻进过程中,筒壁外侧对流换热系数hco为:
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>Nu</mi>
<mi>a</mi>
</msub>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mi>h</mi>
</msub>
<msubsup>
<mi>Re</mi>
<mrow>
<mi>e</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<msup>
<mi>Pr</mi>
<mi>&gamma;</mi>
</msup>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
式(10)中,rgi为隔水管内径,单位为米m;
无隔水管的深水区域,环空流体的温度Ta为海水的温度Tf;对应管柱壁外侧的对流换热系数hco为海水的流动换热系数,hco为:其中,海水努谢尔数Ref为海水粘度系数uf对应的雷诺数;Prf为海水普朗特数,无量纲;系数c的取值范围是0.024~0.88;n的取值范围是0.33~0.805。
12.根据权利要求10所述的方法,其特征在于,所述对传热微分方程进行离散处理包括:
对筒壁的传热微分方程进行离散处理,获得第i+1时刻,包含井身轴向方向上第j节点的筒壁温度的计算方程:
<mrow>
<msub>
<mi>A</mi>
<mi>w</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>B</mi>
<mi>w</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>E</mi>
<mi>w</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>w</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>F</mi>
<mi>w</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>w</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:
<mrow>
<msub>
<mi>A</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>k</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<msup>
<mi>&Delta;z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
<mo>{</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>ln</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>k</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
<mo>}</mo>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
<mo>{</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>ln</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>k</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
<mo>}</mo>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>w</mi>
</msub>
<msub>
<mi>c</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>B</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>k</mi>
<mi>w</mi>
</msub>
<mrow>
<msup>
<mi>&Delta;z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>C</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
<mo>{</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>ln</mi>
<mo>&lsqb;</mo>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>/</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>k</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
<mo>}</mo>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>D</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>w</mi>
</msub>
<msub>
<mi>c</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>k</mi>
<mi>w</mi>
</msub>
<mrow>
<msup>
<mi>&Delta;z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>F</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
<mo>{</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>l</mi>
<mi>n</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>k</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
<mo>}</mo>
</mrow>
</mfrac>
<mo>.</mo>
</mrow>
4
13.根据权利要求2所述的方法,其特征在于,所述环空流体为地层区域环空流体时,所述地层区域的环空流体的传热微分方程为:
或,
<mrow>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>&alpha;c</mi>
<mi>p</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&pi;r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&pi;r</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>h</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>r</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
式(12)或(13)中,Qha为液体流动摩擦生热产生的环空流体的热源,单位为瓦每米W/m,Qha=GPa,Pa为单位长度流动磨损,单位为Pa/m,Pa通过非牛顿流体的流动摩阻计算公式获得;hwe为井壁的对流换热系数,单位为W/mk;Twe为井壁温度,单位为℃。
14.根据权利要求13所述的方法,其特征在于,所述对传热微分方程进行离散处理包括:
对地层区域的环空流体的传热微分方程进行离散处理,获得第i+1时刻,包含井身轴向方向上第j节点的地层区域的环空流体温度的计算方程:
<mrow>
<msub>
<mi>A</mi>
<mi>a</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>B</mi>
<mi>a</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>E</mi>
<mi>a</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>a</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>F</mi>
<mi>a</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>a</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>-</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>h</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:
<mrow>
<msub>
<mi>A</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>Gc</mi>
<mi>p</mi>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>z</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>l</mi>
<mi>n</mi>
<mo>&lsqb;</mo>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>/</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;k</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mfrac>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&pi;r</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>&Delta;z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>r</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>B</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>&Delta;z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>C</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<mn>2</mn>
<msub>
<mi>&pi;r</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>;</mo>
<msub>
<mi>D</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>r</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>r</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>Gc</mi>
<mi>p</mi>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>z</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>&Delta;z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
<msub>
<mi>F</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>ln</mi>
<mo>&lsqb;</mo>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>/</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;k</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mfrac>
<mo>.</mo>
</mrow>
15.根据权利要求13所述的方法,其特征在于,
地层区域内,不钻进过程中,筒壁外侧及井壁内侧的对流换热系数hco和hwe通过非牛顿流体的计算公式获得;
地层区域内,钻进过程中,筒壁外侧的对流换热系数hco:
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>Nu</mi>
<mi>a</mi>
</msub>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mi>h</mi>
</msub>
<msubsup>
<mi>Re</mi>
<mrow>
<mi>e</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<msup>
<mi>Pr</mi>
<mi>&gamma;</mi>
</msup>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
井壁内侧的对流换热系数hwe与筒壁外侧对流换热系数hco相同;
井壁温度Twe由一维稳态传热模型或二维地层导热模型获得。
16.根据权利要求2所述的方法,其特征在于,有隔水管的深水区域,所述环空流体的传热微分方程为:
或
<mrow>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>G</mi>
<mi>s</mi>
</msub>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&pi;r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&pi;r</mi>
<mrow>
<mi>f</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>g</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>G</mi>
<mi>s</mi>
</msub>
<msub>
<mi>P</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,hgi为隔水管内壁面对流换热系数,单位为W/m2K;Tg为隔水管壁面温度,单位为℃;Gs为深水区域环空流体的体积流量:Gs=G+Gsa,其中,Gsa为海底增压流量,单位为m3/s;
对有隔水管的深水区域,钻进过程中,hco=hgi;
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>Nu</mi>
<mi>a</mi>
</msub>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mi>h</mi>
</msub>
<msubsup>
<mi>Re</mi>
<mrow>
<mi>e</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<msup>
<mi>Pr</mi>
<mi>&gamma;</mi>
</msup>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
不钻进过程中,筒壁外侧及隔水管内侧的对流换热系数hco和hgi通过非牛顿流体的计算公式获得。
17.根据权利要求2所述的方法,其特征在于,所述隔水管的传热微分方程为:
或,
<mrow>
<msub>
<mi>k</mi>
<mi>g</mi>
</msub>
<mfrac>
<mrow>
<msup>
<mo>&PartialD;</mo>
<mn>2</mn>
</msup>
<msub>
<mi>T</mi>
<mi>g</mi>
</msub>
</mrow>
<msup>
<mrow>
<mo>&PartialD;</mo>
<mi>z</mi>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mrow>
<mn>2</mn>
<mi>r</mi>
</mrow>
<mi>gi</mi>
</msub>
<msub>
<mi>h</mi>
<mi>gi</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>g</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>r</mi>
<mi>go</mi>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>r</mi>
<mi>gi</mi>
</msub>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mrow>
<mn>2</mn>
<mi>r</mi>
</mrow>
<mi>go</mi>
</msub>
<msub>
<mi>h</mi>
<mi>go</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>g</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>r</mi>
<mi>go</mi>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>r</mi>
<mi>gi</mi>
</msub>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>&rho;</mi>
<mi>g</mi>
</msub>
<msub>
<mi>c</mi>
<mi>g</mi>
</msub>
<mfrac>
<msub>
<mrow>
<mo>&PartialD;</mo>
<mi>T</mi>
</mrow>
<mi>g</mi>
</msub>
<mrow>
<mo>&PartialD;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
式中,Tg为隔水管的温度,单位为℃;Tf为海水温度,单位为℃;hgo为隔水管外壁面的对流换热系数,单位为W/m2k:其中,海水努谢尔数海水雷诺数Prf为海水普朗特数,无量纲;系数c的取值范围是0.024~0.88;n的取值范围是0.33~0.805。
18.根据权利要求17所述的方法,其特征在于,所述对传热微分方程进行离散处理包括:
对隔水管的传热微分方程进行离散处理,获得第i+1时刻,包含井身轴向方向上第j节点的隔水管温度的计算方程:
<mrow>
<msub>
<mi>A</mi>
<mi>g</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>g</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>B</mi>
<mi>g</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>g</mi>
<mo>,</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>E</mi>
<mi>g</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>g</mi>
<mo>,</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>g</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>f</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>F</mi>
<mi>g</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>g</mi>
</msub>
<msubsup>
<mi>T</mi>
<mrow>
<mi>g</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>19</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:
<mrow>
<msub>
<mi>A</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>k</mi>
<mi>g</mi>
</msub>
</mrow>
<mrow>
<msup>
<mi>&Delta;z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
<mo>{</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>ln</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>k</mi>
<mi>g</mi>
</msub>
</mrow>
</mfrac>
<mo>}</mo>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
<mo>{</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>ln</mi>
<mo>&lsqb;</mo>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>/</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>k</mi>
<mi>g</mi>
</msub>
</mrow>
</mfrac>
<mo>}</mo>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>g</mi>
</msub>
<msub>
<mi>c</mi>
<mi>g</mi>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>B</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>k</mi>
<mi>g</mi>
</msub>
<mrow>
<msup>
<mi>&Delta;z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>C</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
<mo>{</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>ln</mi>
<mo>&lsqb;</mo>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>/</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>k</mi>
<mi>g</mi>
</msub>
</mrow>
</mfrac>
<mo>}</mo>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>D</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>g</mi>
</msub>
<msub>
<mi>c</mi>
<mi>g</mi>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>k</mi>
<mi>g</mi>
</msub>
<mrow>
<msup>
<mi>&Delta;z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>F</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
<mo>{</mo>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>ln</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>r</mi>
<mrow>
<mi>g</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>k</mi>
<mi>g</mi>
</msub>
</mrow>
</mfrac>
</mrow>
<mo>}</mo>
</mrow>
</mfrac>
<mo>.</mo>
</mrow>
19.一种获取循环温度场的装置,其特征在于,包括:划分单元、传热微分单元和离散迭代单元;其中,
划分单元用于,将井筒沿管柱径向划分为两个或两个以上组成部分;
传热微分单元用于,根据钻井状态参数,分别获得计算在钻进和不钻进过程中各组成部分的瞬态传热信息的传热微分方程;
离散迭代单元用于,对传热微分方程进行离散和数值迭代处理,获得井筒的瞬态温度分布。
20.根据权利要求19所述的装置,其特征在于,所述划分单元具体用于:
将井筒沿管柱径向划分为管柱内流体、筒壁、环空流体、隔水管;
其中,划分的各所述组成部分分别包括对应的预设数个节点。
21.根据权利要求19或20所述的装置,其特征在于,所述装置还包括获取单元,用于获得所述传热微分方程之前,获取所述钻井状态参数;
其中,所述钻井状态参数包括:时间步长△t、空间步长△z、钻井液入口流量G、入口温度Tin、各时间节点对应的各空间节点位置、各时间节点的钻速、井身结构参数、物性参数、钻头机械磨损参数和/或初始时刻各空间节点温度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710197241.3A CN107145705B (zh) | 2017-03-29 | 2017-03-29 | 一种获取循环温度场的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710197241.3A CN107145705B (zh) | 2017-03-29 | 2017-03-29 | 一种获取循环温度场的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107145705A true CN107145705A (zh) | 2017-09-08 |
CN107145705B CN107145705B (zh) | 2020-08-21 |
Family
ID=59783642
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710197241.3A Active CN107145705B (zh) | 2017-03-29 | 2017-03-29 | 一种获取循环温度场的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107145705B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108319771A (zh) * | 2018-01-25 | 2018-07-24 | 西安石油大学 | 一种低渗透油藏井筒温度计算方法 |
CN109815636A (zh) * | 2019-03-01 | 2019-05-28 | 同济大学 | 一种可用于回转粘度仪的水泥净浆流变性迭代算法 |
CN109869289A (zh) * | 2017-12-05 | 2019-06-11 | 中车株洲电力机车研究所有限公司 | 用于高塔架塔筒的通风冷却性能配置方法、装置及塔筒 |
CN110836766A (zh) * | 2019-10-28 | 2020-02-25 | 江苏科技大学 | 一种双流体喷嘴的检测装置及检测方法 |
CN111814100A (zh) * | 2020-08-06 | 2020-10-23 | 西南石油大学 | 无隔水管海底泵举升钻井系统循环温度变化动态模拟方法 |
CN114048657A (zh) * | 2021-11-16 | 2022-02-15 | 西南石油大学 | 一种水平井钻井井筒温度场分布模型的构建方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008061129A2 (en) * | 2006-11-14 | 2008-05-22 | University Of Utah Research Foundation | Methods and compositions related to continuous flow thermal gradient pcr |
CN102855414A (zh) * | 2012-09-24 | 2013-01-02 | 湖南大学 | 一种竖直地埋管换热器参数化设计的高效计算方法 |
EP2673332A1 (de) * | 2011-02-08 | 2013-12-18 | Wintershall Holding GmbH | Mehrstufiges verfahren zur förderung von erdöl unter verwendung von mikroorganismen |
CN104594834A (zh) * | 2014-12-01 | 2015-05-06 | 中国石油大学(华东) | 一种深水油基钻井液钻井溢流情况的监测方法 |
-
2017
- 2017-03-29 CN CN201710197241.3A patent/CN107145705B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008061129A2 (en) * | 2006-11-14 | 2008-05-22 | University Of Utah Research Foundation | Methods and compositions related to continuous flow thermal gradient pcr |
EP2673332A1 (de) * | 2011-02-08 | 2013-12-18 | Wintershall Holding GmbH | Mehrstufiges verfahren zur förderung von erdöl unter verwendung von mikroorganismen |
CN102855414A (zh) * | 2012-09-24 | 2013-01-02 | 湖南大学 | 一种竖直地埋管换热器参数化设计的高效计算方法 |
CN104594834A (zh) * | 2014-12-01 | 2015-05-06 | 中国石油大学(华东) | 一种深水油基钻井液钻井溢流情况的监测方法 |
Non-Patent Citations (2)
Title |
---|
CHANGHONG HE等: "Method to Suppress the Temperature Dependency of Scale Factor in Measurement-While-Drilling Surveying Utilizing Fiber Optic Gyroscope", 《IEEE》 * |
张建国: "钻井循环过程中井内瞬态温度预测模型", 《钻采工艺》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109869289A (zh) * | 2017-12-05 | 2019-06-11 | 中车株洲电力机车研究所有限公司 | 用于高塔架塔筒的通风冷却性能配置方法、装置及塔筒 |
CN108319771A (zh) * | 2018-01-25 | 2018-07-24 | 西安石油大学 | 一种低渗透油藏井筒温度计算方法 |
CN109815636A (zh) * | 2019-03-01 | 2019-05-28 | 同济大学 | 一种可用于回转粘度仪的水泥净浆流变性迭代算法 |
CN109815636B (zh) * | 2019-03-01 | 2023-07-21 | 同济大学 | 一种可用于回转粘度仪的水泥净浆流变性迭代算法 |
CN110836766A (zh) * | 2019-10-28 | 2020-02-25 | 江苏科技大学 | 一种双流体喷嘴的检测装置及检测方法 |
CN110836766B (zh) * | 2019-10-28 | 2021-11-12 | 江苏科技大学 | 一种双流体喷嘴的检测装置及检测方法 |
CN111814100A (zh) * | 2020-08-06 | 2020-10-23 | 西南石油大学 | 无隔水管海底泵举升钻井系统循环温度变化动态模拟方法 |
CN111814100B (zh) * | 2020-08-06 | 2021-03-19 | 西南石油大学 | 无隔水管海底泵举升钻井系统循环温度变化动态模拟方法 |
CN114048657A (zh) * | 2021-11-16 | 2022-02-15 | 西南石油大学 | 一种水平井钻井井筒温度场分布模型的构建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107145705B (zh) | 2020-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107145705A (zh) | 一种获取循环温度场的方法及装置 | |
Alimonti et al. | Study of geothermal power generation from a very deep oil well with a wellbore heat exchanger | |
Keller et al. | Temperature distribution in circulating mud columns | |
Song et al. | Heat extraction performance of a downhole coaxial heat exchanger geothermal system by considering fluid flow in the reservoir | |
CN104806230B (zh) | 超临界注汽井的井筒温度场及压力场分布的计算方法 | |
CN103226641B (zh) | 深水气液两相流循环温度压力耦合计算方法 | |
CN106968667A (zh) | 一种温度场预测方法及装置 | |
Zhang et al. | Comparative analysis of heat transfer performance of coaxial pipe and U-type deep borehole heat exchangers | |
Xiao et al. | Research on wellbore temperature control and heat extraction methods while drilling in high-temperature wells | |
CN102682195B (zh) | 半潜式平台瞬态钻井井筒温度计算方法 | |
CN107575214A (zh) | 用于注采过程的井筒内温度与压力的预测方法 | |
Akhmadullin et al. | Numerical analysis of downhole heat exchanger designed for geothermal energy production | |
Khaled et al. | Downhole heat management for drilling shallow and ultra-deep high enthalpy geothermal wells | |
CN102877835B (zh) | 一种火驱过程中水平生产井井温预测方法 | |
CN106934106A (zh) | 一种获取固井循环温度的方法及装置 | |
Francke et al. | Thermal–hydraulic measurements and modelling of the brine circuit in a geothermal well | |
CN104866681B (zh) | 高温高压油气斜井关井过程中温度压力数值模拟方法 | |
Wang et al. | Development and application of wellbore heat transfer model considering variable mass flow | |
Lin et al. | Study on two-phase flow and heat transfer in offshore wells | |
Abdelhafiz et al. | Application of a thermal transient subsurface model to a coaxial borehole heat exchanger system | |
Xi et al. | Coupling Effect of Transient Temperature-Pressure on Casing String During Volume Fracturing in Shale Gas Wells | |
Liu et al. | Development of a transient method on predicting multi-annuli temperature of subsea wells | |
Wang et al. | Modeling the flow of carbon dioxide during the drilling of oil, gas, and geothermal energy | |
Fallah | An advanced hydraulic-thermal model for drilling & control of geothermal wells | |
CN105952443A (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 |