一种海洋天然气水合物层钻井井筒温度场计算方法与流程

本发明属于海洋钻井技术领域，具体地，涉及海洋天然气水合物层钻井技术方法，特别是一种海洋天然气水合物层钻井井筒温度场计算方法。

背景技术：

天然气水合物是一种稳定存在于低温、高压环境条件下的晶体化合物，其资源量极为丰富，主要分布在陆地永久冻土带和陆缘外围的海底沉积物中，其中海洋天然气水合物资源量约为陆地冻土带的100倍以上。海洋天然气水合物层处于海底低温、高压环境下，保持较为稳定的赋存状态。随着海洋油气钻采技术的发展，海洋钻井日趋增多，海洋钻井钻进过程中，依次钻至地层、天然气水合物层、地层，其中，钻井液从钻柱中自井口即海平面处注入井筒，流经泥线以上即海水所处井段，至泥线位置即海底处后，通过泥线以下即地层所处井段，从钻头流出至环空中，再通过环空上返至井口，从环空井口返出井筒。

然而，海洋钻井钻遇天然气水合物层时，由于钻头处于海洋天然气水合物层，井筒中流动的钻井液温度较高，会改变海洋天然气水合物层的原始温度场，进而影响水合物的稳定状态，造成水合物的分解气化，对钻井施工安全造成重大威胁。另一方面，海洋天然气水合物层钻井中，进入井筒环空中的水合物钻屑颗粒不同于常规海洋钻井中的岩屑颗粒，其随环空钻井液上返过程中，会由于井筒温度升高、压力降低而发生分解，水合物钻屑颗粒的分解吸热作用会进一步影响井筒温度变化，整个井筒温度场会也会因此改变，进而影响钻井液流变性、井筒流动参数变化。因此，海洋天然气水合物层钻井井筒温度场的准确计算，对判断海洋天然气水合物层稳定状态、钻井液流变性、计算井筒流动参数，进而保障海洋钻井施工安全具有重要的指导意义。

目前，国内外对海洋天然气水合物层钻井井筒温度场研究较少，现有研究主要针对海洋钻井过程中的井筒温度场，并未研究钻遇海洋天然气水合物层时，环空中水合物钻屑颗粒随钻井液上返并发生分解下的井筒温度场。专利cn103226641a公开了一种深水气液两相流循环温度压力耦合计算方法，通过先钻柱内钻井液后环空钻井液的顺序，迭代计算钻柱和环空内钻井液节点温度和压力数据，最终得到深海气液两相流井筒温度和压力模拟结果，此方法可以应用在海洋钻井过程中井筒温度场的计算，但是不能准确反映海洋天然气水合物层钻井时的井筒温度场；专利cn102943620a公开了一种基于钻井环空井筒多相流动计算的控压钻井方法，井筒环空多相流动控制方程组求解中，所需的井筒温度有所体现，但具体求解方法并未提出，不能直接应用在海洋天然气水合物层钻井井筒温度场计算中。因此，迫切需要一种针对海洋天然气水合物层钻井的井筒温度场计算方法，当钻遇海洋天然气水合物层时考虑环空中水合物钻屑颗粒分解的影响，为判断海洋天然气水合物层稳定状态、钻井液流变性、计算井筒流动参数提供理论依据。

技术实现要素：

本发明的目的在于提供一种海洋天然气水合物层钻井井筒温度场计算方法，该方法原理可靠，便于操作，能够为判断海洋天然气水合物层稳定状态、钻井液流变性、计算井筒流动参数，进而保障海洋钻井施工安全提供理论依据，具有广阔的市场前景。

为了实现上述技术目的，本发明采用以下技术方案。

获取海洋天然气水合物层钻井参数和初始条件，在此基础上进行空间节点划分；考虑天然气水合物钻屑颗粒分解吸热的影响，建立海洋天然气水合物层钻井井筒温度场计算模型，得到计算网格内的温度变化；然后，根据已知节点处的温度，计算下一节点处的井筒温度，按照钻柱内和环空内的节点同时自井口至井底的顺序，迭代计算井筒温度，直至井底温度满足计算误差，则迭代结束；迭代计算所得所有节点处的井筒温度即为海洋天然气水合物层钻井井筒温度场。

一种海洋天然气水合物层钻井井筒温度场计算方法，依次包括以下步骤：

(1)根据钻井设计和储层参数，获取海洋天然气水合物层钻井参数和初始条件，钻井参数包括：井身结构、钻具组合、泵入参数、海水深度、海水温度、地层温度、储层中水合物丰度、机械钻速、井深，初始条件包括：计算时刻的钻井液注入温度tp(0)⁽ⁿ⁾、环空井口压力pa(0)⁽ⁿ⁾。

(2)进行空间节点划分，根据步骤(1)中海洋天然气水合物层钻井参数，空间域为整个井筒，自井口至井底，节点轴向序号从0开始依次递增；井口处的节点序号为0，井筒中任一计算节点序号用i表示，下一计算节点序号用i+1表示，井底处节点序号为k。

(3)根据传热学理论及能量守恒定律，建立海洋天然气水合物层钻井井筒温度场计算模型，计算网格内的温度变化表达式如下(高永海，孙宝江，王志远等.深水钻探井筒温度场的计算与分析[j].中国石油大学学报(自然科学版)，32(2003)：58-62)：

钻柱中：

环空中：

泥线(海底处)以上井段，即i·δh＜hsea时：

泥线(海底处)以下井段，即i·δh≥hsea时：

式中：i为计算节点；

n为计算时刻；

δh为计算网格长度，计算中设定为1m；

hsea为海水深度，m；

δtp(i)⁽ⁿ⁾、δta(i)⁽ⁿ⁾分别为钻柱中、环空中计算网格内的温度变化，k；

ρp(i)、ρa(i)分别为钻柱中、环空中计算节点i处的混合流体密度，kg/m³；

vp(i)、va(i)分别为钻柱中、环空中计算节点i处的混合流体流速，m/s；

cp(i)、ca(i)分别为钻柱中、环空中井筒中计算节点i处的混合流体比热容，j/(kg·k)；

dpi、dpo、dri、dci分别为钻柱内径、钻柱外径、隔水管内径、套管内径，m；

qap(i)为环空与钻柱内计算节点i处流体间热交换，w；

qwa(i)、qsa(i)分别为海水与环空内、地层与环空内计算节点i处流体间热交换，w；

qfp(i)、qfa(i)分别为钻柱中和环空中计算节点i处流动摩擦产生的热量，w；

qh(i)为环空中计算节点i处水合物分解吸收的热量，w。

(4)井口处即i＝0时，钻柱中井口处的流体温度即为计算时刻的钻井液注入温度tp(0)⁽ⁿ⁾，根据步骤(1)中的初始条件，为已知参数；而环空中井口处的流体温度ta(0)⁽ⁿ⁾即为从环空中返出流体的温度是未知参数，假设ta(0)⁽ⁿ⁾的值，设定假设范围：273k≤ta(0)⁽ⁿ⁾≤tp(0)⁽ⁿ⁾。

(5)按照钻柱内和环空内的节点同时自井口至井底的顺序，根据节点i处钻柱中流体温度tp(i)⁽ⁿ⁾和环空中流体温度ta(i)⁽ⁿ⁾，计算下一节点i+1处的井筒温度：

式中：tp(i)⁽ⁿ⁾、ta(i)⁽ⁿ⁾分别为计算节点i处钻柱中、环空中流体温度，k；

tp(i+1)⁽ⁿ⁾、ta(i+1)⁽ⁿ⁾分别为计算节点i+1处钻柱中、环空中流体温度，k。

根据公式(4)、(5)，自井口迭代计算直至井底节点k处，得到钻柱中井底处的流体温度tp(k)⁽ⁿ⁾和环空中井底处的流体温度ta(k)⁽ⁿ⁾。

(6)根据步骤(5)迭代计算得到的钻柱中、环空中井底处的流体温度，比较是否满足计算误差：

式中：tp(k)⁽ⁿ⁾为钻柱中井底处的流体温度，k；

ta(k)⁽ⁿ⁾为环空中井底处的流体温度，k；

γ为井底处流体温度计算误差，取1k。

如果公式(6)成立，则满足计算误差，通过步骤(5)迭代计算所得所有节点处钻柱中、环空中流体温度即为海洋天然气水合物层钻井井筒温度场。如果公式(6)不成立，则不满足计算误差，需对步骤(4)中的环空中井口处流体温度ta(0)⁽ⁿ⁾重新假设，并再次重复步骤(5)迭代计算，直至公式(6)成立。

所述步骤(3)中，环空与钻柱内计算节点i处流体间热交换qap(i)、海水与环空内计算节点i处流体间热交换qwa(i)、地层与环空内计算节点i处流体间热交换qsa(i)、钻柱中计算节点i处流动摩擦产生的热量qfp(i)、环空中计算节点i处流动摩擦产生的热量qfa(i)(z.m.wang,x.n.hao,x.q.wangetal.numericalsimulationondeepwaterdrillingwellboretemperatureandpressuredistribution[j].petroleumscienceandtechnology,28(2010):911–919)以及环空中计算节点i处水合物分解吸收的热量qh(i)(e.d.sloana,f.fleyfelb.hydratedissociationenthalpyandguestsize[j].fluidphaseequilibria,76(1992):123-140)的计算方法如下：

环空与钻柱内计算节点i处流体间热交换qap(i)计算如下

海水与环空内计算节点i处流体间的热交换qwa(i)计算如下

地层与环空内计算节点i处流体间的热交换qsa(i)计算如下

式(9)中的环空内流体与地层在计算节点i处的综合换热系数usa(i)计算如下

钻柱中计算节点i处流动摩擦产生的热量qfp(i)计算如下

环空中计算节点i处流动摩擦产生的热量qfa(i)计算如下

环空中计算节点i处水合物分解吸收的热量qh(i)计算如下

式(10)～(13)中：

tw(i)⁽ⁿ⁾、ts(i)⁽ⁿ⁾分别为计算节点i处的海水温度、地层温度，k；

dro、dco、dcso、dcsi分别为隔水管外径、套管外径、水泥环外径、水泥环内径，m；

αf1(i)、αf2(i)、αf3(i)分别为钻柱内表面、隔水管内表面、套管内表面上计算节点i处的受迫对流换热系数，w/(m²·k)；

αm1(i)、αm2(i)分别为钻柱外表面、隔水管外表面上计算节点i处的外绕面对流换热系数，w/(m²·k)；

λp(i)、λr(i)、λc(i)、λcs(i)、λs(i)分别为钻柱、隔水管、套管、水泥环、地层在计算节点i处的导热系数，w/(m·k)；

usa(i)为环空中流体与地层在计算节点i处的综合换热系数，w/(m²·k)；

td为瞬态传热函数，无因次；

vp(i)、va(i)分别为钻柱中、环空中混合流体在计算节点i处的流速，m/s；

mp(i)、ma(i)分别为钻柱中、环空中混合流体在计算节点i处的质量流量，kg/s；

dao为环空外径，计算节点位于泥线以上井段即为隔水管内径dri，位于泥线以下井段即为套管内径dci，m；

zgeq为节点处环空中温度、压力条件下的天然气压缩因子，无因次；

r为通用气体常数，j/(mol·k)；

peq(i)⁽ⁿ⁾为计算节点i处天然气水合物相平衡压力，pa。

与现有技术相比，本发明具有以下显著优点：

(1)本发明通过按照钻柱内和环空内的节点同时自井口至井底的顺序，迭代计算得到井筒温度场，计算便捷，误差较小，能够准确快速计算出井筒温度场。

(2)本发明能够针对海洋天然气水合物层钻井实现其井筒温度场计算，为判断海洋天然气水合物层稳定状态、钻井液流变性、计算井筒流动参数，进而保障海洋钻井施工安全提供理论依据。

附图说明

图1为海洋天然气水合物层钻井井筒换热示意图。

图2为海洋天然气水合物层钻井井筒温度场实例计算结果图。

具体实施方式

下面以某地区实际海洋天然气水合物层钻井为例，结合附图，对本发明作进一步详细地说明，但本发明并不局限于以下实施例。

海洋天然气水合物层钻井井筒换热示意图如图1所示，泥线以上井筒外部为海水，泥线以下井筒外部为地层。钻井液在井口从钻柱中注入，通过钻柱与环空中流体发生热交换。到达井底之后，钻井液携带天然气水合物钻屑颗粒从环空中上返，一方面，环空中流体与钻柱中流体发生热交换；另一方面，泥线以下井段环空中流体通过套管、水泥环与地层发生热交换，泥线以上井段环空中流体通过隔水管与海水发生热交换；上返过程中，天然气水合物钻屑颗粒随着井筒温度升高、压力降低而发生分解吸热，进而影响井筒温度场。

一种海洋天然气水合物层钻井井筒温度场计算方法，依次包括以下步骤：

(1)根据钻井设计和储层参数，获取海洋天然气水合物层钻井参数和初始条件：隔水管外径dro、套管外径dco均为0.508m；钻柱外径dpo为0.127m，钻头直径为0.445m；井深h为1600m，泵入钻井液排量30l/s，泵入钻井液密度ρp(0)为1030kg/m³；海水深度hsea为1500m；海面温度为298k；地层温度梯度为3℃/100m；储层中水合物丰度为70％；机械钻速为10m/h；计算时刻的钻井液注入温度tp(0)⁽ⁿ⁾为298k；环空井口压力pa(0)⁽ⁿ⁾为101300pa。

(2)进行空间节点划分，根据步骤(1)中海洋天然气水合物层钻井参数，空间域为整个井筒，自井口至井底，节点轴向序号从0开始依次递增；井口处的节点序号为0，井底处节点序号为1600。

(3)根据传热学理论及能量守恒定律，基于海洋天然气水合物层钻井井筒温度场计算模型，计算网格内的温度变化：钻柱中δtp(i)⁽ⁿ⁾和环空中δta(i)⁽ⁿ⁾。

(4)井口处即节点0处，进行计算时，钻柱中井口处的流体温度即为计算时刻的钻井液注入温度tp(0)⁽ⁿ⁾为298k；而环空中井口处的流体温度ta(0)⁽ⁿ⁾即为从环空中返出流体的温度是未知参数，进行假设ta(0)⁽ⁿ⁾的值，根据假设范围：273k≤ta(0)⁽ⁿ⁾≤tp(0)⁽ⁿ⁾假设为292k，并进行计算。

(5)按照钻柱内和环空内的节点同时自井口至井底的顺序，根据节点i处钻柱中流体温度tp(i)⁽ⁿ⁾和环空中流体温度ta(i)⁽ⁿ⁾，计算下一节点i+1处的井筒温度：钻柱中流体温度tp(i+1)⁽ⁿ⁾和环空中流体温度ta(i+1)⁽ⁿ⁾；自井口迭代计算至井底处，得到钻柱中井底处的流体温度tp(k)⁽ⁿ⁾＝281.06k，环空中井底处的流体温度ta(k)⁽ⁿ⁾＝281.41k。

(6)根据步骤(5)中井筒温度计算方法，得到的钻柱中井底处的流体温度tp(k)⁽ⁿ⁾、环空中井底处的流体温度ta(k)⁽ⁿ⁾，并根据计算误差γ＝1k，比较是否满足计算误差：则满足计算误差。

因此，环空中井口处的流体温度ta(0)⁽ⁿ⁾为292k时，通过计算，井底处流体温度满足计算误差，通过步骤(5)迭代计算所得所有节点处钻柱中、环空中流体温度即为海洋天然气水合物层钻井井筒温度场，如图2所示(图2为海洋天然气水合物层钻井井筒温度场实例计算结果图)。