辰工科技使用随钻录井资料计算地层静温

摘　要　根据随钻录井中热力学参数、泥浆参数及井身结构参数等,并给地层温度梯度一个初值,从热力学及流体力学等有关方程出发,经过推演得到井壁上温度随深度变化以及地层温度分布的数学模型用于计算泥浆出口温度。将此计算值与实测的泥浆出口温度值比较,根据比较结果再修正地层温度梯度,如此反复,直至计算值与实测的泥浆出口温度值相等,从而得到钻头所在的初始地层静温。由于钻井过程中泥浆、岩石及其温度场间是相互作用、相互影响的,这为研究热—流—固耦合过程的理论与应用提供了一种新的方法。

主题词　钻井　温度梯度　随钻录井　钻井液　热力学　井身结构

钻井时,由于地层温度梯度的影响,钻井液通过与地层交换热量,使得地层出现温度分布不均现象,这种分布使钻井液温度有所升高,导致泥浆入口与进口温度有所不同。如何通过入口与进口泥浆温度差来反算地层的静温是录井工作中一个很有意义的研究课题。

　地层静温对石油的勘探与开发非常重要（详见http://www.successly.com）,目前一般使用测井的温度数据,采用Horner 半对数方法得到[1 ,2 ,3], 这种方法不仅费时, 而且很多情况下Horner 半对数图上不出现直线段,得到的地层静温也就不准确。本文利用录井数据,将钻井液视为液固两相非牛顿流体,并考虑流体的自然对流传热、流体与固壁之间的表面换热及岩石的热传导[4 ,5],通过数值计算方法得到井底静温,不需要特殊的工艺。由泥浆入口温度和出口温度计算井底静温的原理是:在所有热力学参数、泥浆有关参数及井身结构参数等已知的情况下,给一个地层温度梯度,可通过有关计算得到泥浆出口温度,将这个计算温度与实测的泥浆出口温度进行对比,如果两者不相符合,则改变地层温度梯度后再进行计算。如此反复,直至计算的出口温度与实测的泥浆出口温度一致为止,此时的地层温度梯度可被用来计算该井深下的地层静温。因此,研究地层静温的关键是从热力学及流体力学有关方程出发,通过数值差分或其他方法得到井壁上的温度随深度的变化及地层的温度分布。

一、数学模型的建立

(1) 流体的视粘度(有效粘度) :钻井液可认为是不可压流体,但由于泥浆的粘稠性,不可能用牛顿流体来处理,必须采用非牛顿流体来处理,对非牛顿流体这里研究的是宾厄姆流体和幂律流体。如高固相泥浆及加重泥浆,其流变特征近似宾厄姆流体,低固相泥浆可加入稀释剂,其流变特征近似幂律流体。

对于幂律流体,其有效粘性系数μeff 可定义如下:

                       (1)

式中:M 幂律流体的粘性系数(MPa) , n 幂律流体指数系数,γ剪切应变力(L/ s) 。对于流体速度为v ,半径为R 的圆管,剪切应变力γ可写成:

                   (2)

对于内外径为R1、R2 的环形管道其剪切应变力γ可写成:

（3）

对于宾厄姆流体,有效粘度系数μeff 定义为:

（4）

式中:μp为塑性粘度（mPa·s）,τ0 为屈服点( Pa) 。流体体积流量为Q ,半径为R 的圆管,剪切应变力γ可写成:

（5）

对于内、外径分别为R1 、R2 的环形管道γ可写成:

（6）

(2) 流体的对流过程:由于泥浆是液固混合物,那么密度为ρ的泥浆中固相份额为S F = (ρ-1000) /ρ, 泥浆的热容量C = 4190.5-3256.1S F(J /kg ·K) ,泥浆的导热系数K = 0.69+16.61S F(W/m ·K) 。当泥浆流过固壁时,泥浆与固壁之间的传热称为对流传热,井筒内的对流传热由表面对流系数h 来决定,对直径为D 的管道,表面对流系数h =

其中Nusselt 数NNu 可用Reynold 数NRe 与Prandlt 数Npr来表示,实验可得出三个无量纲数之间关系:

式中：Reynol数

，Prandlt 数

沿程摩擦系数可表示如下:

（8）

(3) 流体的自然对流传热:环空中的流体自然对流使流体沿径向的热导系数增大,由实验可得出环状空间中自然对流的有效导热系数值Keff ：

（9）

式中

为膨胀系数数, 1/ K;μ为流体粘度,Pa·s;ΔT为环状空间径向温度差,K。这样流体的自然对流传热就等效于导热。

（4）井筒和地层之间的传热：确定井及地层的温度分布是一件复杂的工作，这是因为影响温度分布的因素很多，包括井筒内流体的流动类型、井身结构、流体、井壁、水泥、岩石的热导特性等。这些因素中的许多参数又是随时间变化的,联立求解能量方程、动量方程、连续方程和状态方程十分困难。经过大量数值实验,如果在能量方程中略去了与其它方程的耦合项,则精度仅降低5 % ,但可以独立求解相关方程。为此,对能量方程做了简化,略去了与其它方程耦合项,从而可以独立求解相关方程,以确定井壁及地层的温度分布。从钻铤到地层可能要经过多层套管及套管的水泥层,它们之间的传热是较复杂的,必须要考虑导热、对流(包括表面对流和自然对流) ,但辐射项可以忽略。这样热传导差分形式的方程为:

（10）

式中:

传导系数, h 是流体表面对流系数,W/ m2 ·K; K1 , K2 , ⋯为各环形介质的导热系数,W/ m ·K; r1 , r2 , ⋯介质交界处的半径,m ;ΔT 温度差,K;ΔZ 垂向距离,m。固体的导热系数等随温度变化较小,可以看作常数,而流体的表面对流系数和导热系数都是温度函数,这样方程(10) 就成为非线性的,实际计算中可以用上一步的温度求出系数,从而使方程线性化。(5) 井筒内温度分布的差分:当钻井时,处于钻进期间的流体,其能量方程可表示成:

（11）

式中: h 为流体(泥浆) 的焓,Φ为势函数。方程(11)的差分形式为:

  （12）

式中: A 为流体占据的面积。动量方程可表示成:

（13）

方程(13) 的差分形式为:

（14）

将方程(14) 代入方程(12) ,考虑到研究区域内无热源,并略去两个方程中的耦合项及耗散项Φ后得到如下形式的差分方程:

（15）

式中:

下标f 表示流体,下标st 表示钢铁。

　(6) 岩石中的温度分布计算:岩石中的传热较简单,只有导热项,而无对流传热项,考虑径向及轴向传热,经整理后,差分形式的方程为:

（16）

式中:

(7) 网络划分及计算方法概述:在本文中,我们将地层、井身结构体系视为径向和轴向的二维区域,地层中的轴向网格结点与井筒中处于同一高度。轴向网格的格长度为20 m ,钻头位置向下扩展2 个网格。而径向温度在井筒附近变化最快,因此,井筒附近需要进行网格加密,网格的宽度随半径呈指数增长。如果轴向网格数为NZ ,径向网格数为Nr ,则地层中的总结点数为NZ ×( Nr - 3) ,井筒中径向有3个网格,故井筒中有3 NZ 个方程,而地层中有NZ ×( Nr - 3) 个方程。本文中采用Gauss - Seide 迭代法同时求解井筒中的3 NZ 个方程及地层中NZ ×( Nr- 3) 个线性方程组,作为边界条件轴向最底下的一个网格点的温度可保持不变。采用Gauss - Seide 迭代法,有些格点的收敛速度快于另外的格点,程序中自动检查,使这些点不再参加迭代,以便加快计算速度。另外,合适选取松驰因子,也可加快收敛速度。

二、在录井中的应用

在钻井中随钻录井可以获取泥浆及地层的相关参数,根据随钻录井中的泥浆出口、入口的密度、温度、体积等数据,采用表1 所示的套管、水泥环及岩石的热力学参数,通过计算可以得到地层静温。由分析可看出,在地层深度较小时,计算出来的地层静温数据振荡较大,随着深度增加计算的地层静温明显呈现出一定的规律性。产生这种情况的主要原因,当井深较浅时,由于泥浆还来不及与地层交换热量,就已经流出地面,因此计算出来的地层静温误差较大。但当深度较深时,泥浆与地层充分交换热量,泥浆出口温度能充分反应地层静温,这种情况下计算出来的地层静温,准确率较高。

三、结　论

(1) 本文通过地层温度梯度初值计算泥浆出口温度,并将此温度与实测的泥浆出口温度比较, 通过迭代不断修正地层温度梯度,直至计算值与实测值相同。以此来计算钻头所在处的地层静温;

(2) 通过数值差分方法,求解岩石中温度分布,并计算出各点的温度;

(3) 计算所得的地层静温为钻井液的配制、钻头速度的确定等钻井工艺数的优化选定,提供了理论依据。