CN102682195A - 半潜式平台瞬态钻井井筒温度计算方法 - Google Patents

Abstract

本发明公开了一种半潜式平台瞬态钻井井筒温度计算方法，包括以下步骤：1)应用初始条件海水区所有节点的初始温度为对应深度处海水温度；地层区所有节点的初始温度为对应深度处地层原始温度；2)根据钻井液总模拟循环时间和时间步长划分时间段，从初始时刻开始，累加时间，每一时间步的温度计算需要以上一时间步的温度场数据作为初始条件，迭代计算，直至收敛。同时提供了实现该方法的模拟器，本发明所述的模拟器是基于离散格式表示的二维全瞬态数学模型，并采用隐式的有限差分法进行求解，提高了计算的精度。

Description

半潜式平台瞬态钻井井筒温度计算方法

技术领域

本发明涉及一种温度计算方法，具体涉及一种半潜式平台瞬态钻井井筒温度计算方法。 

背景技术

半潜式钻井平台钻进作业时井筒的温度，特别是处于低温高压构造时，准确预测钻井液、水泥浆、以及井筒周围地层的温度对合理的钻井液设计至关重要。 

目前，国内外学者对井筒传热问题的研究主要集中在陆地井筒，即井筒内传热和地层传热问题。井筒换热研究成果主要有稳态模型和瞬态模型。全瞬态法原始模型由Raymond(1969)提出，其后，Keller(1973)在换热模型中加入了钻井液摩阻和机械能损失引起的内热源，Gary R.Wooley(1980)将全瞬态法应用于热流体注入和采油作业，David W.Marshall and Ramon G.Bentsen(1982)提出了采用直接接技术求解离散方程组进行全瞬态井筒与地层换热研究，B.Corre等(1984)进行了考虑井深随换热时间不断变化的全瞬态井筒温度分布研究，Beirute(1991)开发了用于估算循环和关井情况下温度计算的模拟器，Garcia(1998a)开发了一种热采数值模拟程序。 

半潜式平台钻井作业期间，钻井液不断循环，井筒不断加深，海水段井筒通过隔水管与海水换热，海水对钻井液的持续冷却和入口温度的持续改变，导致半潜式平台独特的瞬态井筒温度场。而现有温度场模拟器通常假设入口温度不变，尤其是没有描述海水段井筒换热，也没有考虑温度对钻井液热物性的影响，不能用于模拟半潜式平台钻井循环温度场。 

发明内容

为了解决上述技术问题，本发明提供一种半潜式平台瞬态钻井井筒温度计算方法。基于海洋钻井液循环瞬态换热特征，采用全隐式有限体积法离散格式，耦合多换热区域，考虑钻井热源，开发了WHTSubmersible模拟器，并结合现场数据验证了模拟器的有效性。 

其特征在于，包括以下步骤： 

1)应用初始条件 

海水区所有节点的初始温度为对应深度处海水温度；地层区所有节点的初始温度为对应深度处地层原始温度； 

2)根据钻井液总模拟循环时间和时间步长划分时间段，从初始时刻开始，累加时间， 每一时间步的温度计算需要以上一时间步的温度场数据作为初始条件，迭代计算，直至收敛。 

进一步优选，所述步骤2)的温度场计算步骤如下： 

A用二维数组Node记录上一时间步结束时的温度数据； 

B将节点温度数据复制一份保存至二维数组NodeNew中； 

C收敛状态记为“假”； 

D以Node中的数据作为初始值按照钻柱内钻井液、钻柱管体、钻柱外环空、井筒界面、隔水管、地层的顺序按列依次计算各传热对象内所有节点的新时间步的温度数据，每一列节点的每个节点计算结束时都要比较新计算出的节点温度与NodeNew中保存的温度数据的差值，并记录在变量TDiff二维数组中，然后将新计算出的温度数据存入二维数组NodeNew中； 

E找出数组TDiff中的最大值，若该值很小，则认为该时间步温度计算达到收敛，结束该时间步的迭代，将NodeNew中的温度数据复制到Node数组中，进入下一时间步；否则，重复步骤D直至收敛。 

一种本发明所述方法的WHTSubmersible模拟器，包括以下19个模块： 

DataInput函数是WHTSubmersible程序的数据输入函数，完成模拟井所有数据的输入； 

TPField函数是WHTSubmersible程序的总功能模块，完成深水井筒温度场计算及数据存储功能； 

GridGeneration函数根据模拟井的井身结构、钻具结构、海水深度对温度场求解域进行轴向和径向网格划分，存储网格节点的轴向和径向几何信息及介质信息； 

TOriginGeneration函数根据海水和地层的垂直温度分布数据插值产生节点轴向深度处的原始温度；模拟器有两种方式建立原始温度剖面：根据海水和地层的温度分布关系式；根据实测海水和地层温度数据建立温度剖面。海洋油井海水段存在温跃层，地温也存在多个梯度。 

Ini函数对求解域内网格节点变量应用初始条件，赋初值； 

TInDrillStem函数、TDrillStem函数、TAnnu函数、TInterface函数、TRiser函数和TFormation函数分别用于计算钻柱内钻井液、钻柱管体、环空钻井液、环空外壁处、隔水管绝热层内和地层内所有节点的温度； 

ThermalPhysics函数。该模块根据节点处钻井液的温度和压力计算节点处钻井液密度、钻井液表观粘度、钻井液比热、钻井液热导率和高温高压密度； 

HeatGeneration函数。该模块考虑钻井液摩阻压降生热、钻头破岩摩擦生热、钻柱与井 壁摩擦生热计算钻柱内和环空这两个换热区域的内热源； 

HPipe函数、HAnnu函数、HAcross函数的功能分别是计算管流、环空流与内外壁两个的、液体横掠圆管的强迫对流换热系数； 

ViscTM函数用于计算节点温度条件下钻井液的表观粘度； 

Kbound函数可用于计算轴向或径向相邻节点控制体界面处的介质热导； 

TriDiag函数用来求解具有三对角矩阵形式的代数方程组。 

本发明的有益效果： 

(1)本发明所述的模拟器是基于离散格式表示的二维全瞬态数学模型，并采用隐式的有限差分法进行求解，提高了计算的精度； 

(2)本发明的技术方案对每个传热对象按照有限体积法沿轴向和径向进行网格划分，确定了程序的求解过程，并介绍了模拟器相关模块的功能； 

(3)本发明中的模拟器WHTSubmersible不仅可以用于海洋半潜式平台钻井循环温度计算，也可以用来计算陆地井筒温度场，并校核温度场；而且模拟器的动态链接库技术可以将其与其他计算软件相结合，这样在一定程度上提高了模拟器的实用性、准确性和高效性； 

(4)利用本发明技术所述模拟器的计算数据与现场实例数据进行了对比，验证了该模拟器的准确性，计算误差不超过5％，可以投入现场使用。 

附图说明

图1为半潜式钻井平台作业时的物理模型图； 

图2为钻柱内和钻柱管体轴向网格划分图； 

图3为地层内径向网格划分图； 

图4为本发明方法计算流程图； 

图5为本发明模拟器的模块结构图。 

具体实施方式

下面结合附图具体实施方式对本发明的方法作进一步详细地说明。 

半潜式钻井平台作业时的物理模型如图1所示，钻柱内流体区。钻井流体在地面流入钻柱内，沿钻柱向下一直到井底；钻柱管体区；钻柱与井壁之间的环空区。钻井流体从井底进入环空，向上流动，直至地面；井壁区；地层区；隔水管区； 

物理模型的分析说明循环流程可以看作一个热交换系统。在这个过程中，流体沿钻柱流下，再由环空返出。 

将井筒传热区域看作轴对称二维求解域，任意节点由两个序号标记(i，j)，j为轴向序 号，自井口向井底节点轴向序号递增，井口节点轴向序号为0；i为节点径向序号，自井筒轴线向地层或隔水管方向递增，钻柱内钻井液径向节点序号为0。 

轴向网格划分时，依据井身结构和钻具结构进行轴向分段，每一段内涉及到的所有传热对象的几何尺寸只有一种。然后自上至下对每一段再根据段长进行轴向网格划分，网格节点轴向坐标位于轴向控制区域中心。轴向划分结束后，从上至下对每一段进行径向网格划分。钻柱内钻井液、钻柱管体、钻柱外环空三个传热对象内各分配一个径向节点，径向序号分别为0、1、2，节点径向坐标位于传热对象中心，环空外径处布置一个径向节点，径向序号为4，对环空外传热区域则根据每一段的套管层次情况增加相应的径向节点，直至进入地层，进入地层后，径向节点间距逐步增大，直至布置到径向坐标超过50米。 

海水段径向节点布置。半潜式平台作业时，环空钻井液通过隔水管与海水进行对流换热，环空节点与海水之间只有隔水管中心一个节点的控制体。 

依据半潜式平台钻井井筒换热机理，将换热对象分为：钻柱内钻井液、钻柱、环空钻井液、隔水管、地层(含套管及之间的水泥环)、界面(包括井壁与地层、井筒和海水)。 

1钻柱内换热模型： 

钻柱内(#1区)轴向节点和网格划分如图2所示，节点位于网格中心，温度控制方程如下： 

$A_{1,j}{T}_{1,j}^{n+1}=A_{1,j-1}{T}_{1,j-1}^{n+1}+A_{2,j}{T}_{2,j}^{n+1}{A}_{1,j}^0{T}_{1,j}^n+S_{1,j},(0\≤r\≤r_1)---\left(\text{1}\right)$

式中： $A_{1,j}=A_{2,j}+m_1{c}_{1,j+1/2}\mathrm{\Δt}+A_{1,j}^0;$ A_2，j＝2πr_1，jw₁Δz_jΔt；A_1，j-1＝m₁c_1，j-1/2Δt； 

$A_{1,j}^0={\mathrm{\ρ}}_{1,j}{\mathrm{\πr}}_{1,j}^2{c}_{1,j}{\mathrm{\Δz}}_j=0;$ $S_{1,j}={\stackrel{\·}{q}}_1{\mathrm{\Δz}}_j\mathrm{\Δt}=0$

上述方程中，T_j-1代表U界面处温度，用T_j代表D界面处温度。初始条件和边界条件为： 

初始条件： $T_{i,j}^1=T_w^{\∞}\left(z_{i,j}\right),z_{i,j}\≤Z_b;$ $T_{i,j}^1=T_e^{\∞}\left(z_{i,j}\right),z_{i,j}>Z_b$

边界条件： 

(T_E为钻井液入口温度，等于上一时刻环空钻井液出口温度 

)； 

$T_{1,j\max }^n=T_{2,j\max }^n=T_{3,j\max }^n;$ $\frac{\∂T}{\∂z}{|}_{z=Z_b}=0;$ $T_{i\max ,j}^n=T_f^{\∞}\left(z_{i\max ,j}\right)$

2钻柱管体换热模型： 

$A_{2,j-1}{T}_{2,j-1}^{n+1}+B_{2,j}{T}_{2,j}^{n+1}+C_{2,j+1}{T}_{2,j+1}^{n+1}=D_{1,j}{T}_{1,j}^{n+1}+D_{3,j}{T}_{3,j}^{n+1}+D_{2,j}^0{T}_{2,j}^n,(r_1\≤r\≤r_2)---\left(2\right)$

式中： 

$B_{2,j}=-A_{2,j-1}-C_{2,j+1}+D_{1,j}+D_{3,j}+D_{2,j}^0;$ $A_{2,j-1}=\frac{{-k}_{2,j-1/2}\mathrm{\Δt}}{\mathrm{\Δ}{Z}_{j-1/2}};$ $C_{2,j+1}=\frac{{-k}_{2,j+1/2}\mathrm{\Δt}}{\mathrm{\Δ}{Z}_{j+1/2}};$

$D_{1,j}=\frac{{2w}_{1,j}{r}_{1,j}{\mathrm{\ΔZ}}_j\mathrm{\Δt}}{r_{2,j}^2-r_{1,j}^2};$ $D_{3,j}=\frac{{2w}_{2,j}{r}_{2,j}{\mathrm{\ΔZ}}_j\mathrm{\Δt}}{r_{2,j}^2-r_{1,j}^2};$ $D_{2,j}^0={\mathrm{\ρ}}_2{c}_{p2}{\mathrm{\ΔZ}}_j$

上述方程中，U界面： 

D界面： 

(1)式的边界条件依然适用。 

3环空换热模型： 

$C_{3,j}{T}_{3,j}^{n+1}=C_{3,j+1}{T}_{3,j+1}^{n+1}+C_{2,j}{T}_{2,j}^{n+1}{C}_{4,j}{T}_{4,j}^{n+1}+{C_{3,j}^0{T}_{3,j}^n+S}_{3,j},(r_2\≤r\≤r_3)---\left(3\right)$

式中： $C_{3,j}=C_{2,j}+C_{4,j}-m_3{c}_{3,j-1/2}\mathrm{\Δt}+C_{3,j}^0;$ C_3，j+1＝m₃c_3，j+1/2Δt；C_2，j＝2πw_2，jr_2，jΔZ_jΔt； 

C_4，j＝2πw_3，jr_3，jΔZ_jΔt； $C_{3,j}^0=\mathrm{\π}(r_{3,j}^2-r_{2,j}^2){\mathrm{\ρ}}_3{c}_{p3}{\mathrm{\ΔZ}}_j;$ $S_{3,j}={\stackrel{\·}{q}}_{3,j}{\mathrm{\ΔZ}}_j\mathrm{\Δt}$

上述方程中，用T_j代表U界面处温度，用T_j+1代表D界面处温度。 

4隔水管换热模型： 

$A_{r,j-1}{T}_{r,j-1}^{n+1}+B_{r,j}{T}_{r,j}^{n+1}{+C}_{r,j+1}{T}_{r,j+1}^{n+1}=D_{3,j}{T}_{3,j}^{n+1}+D_{w,j}{T}_{w,j}^{n+1}+D_{r,j}^0{T}_{r,j}^n,(r_3<r\&le;r_{\mathrm{ins}})---\left(4\right)$

式中： $B_{r,j}=-A_{r,j-1}-C_{r,j+1}+D_{3,j}+D_{w,j}+D_{r,j}^0;$ $A_{r,j-1}=\frac{-k_{r,j-1/2}\mathrm{\&Delta;t}}{{\mathrm{\&Delta;Z}}_{j-1/2}};$ $C_{r,j+1}=\frac{{-k}_{r,j+1/2}\mathrm{\&Delta;t}}{{\mathrm{\&Delta;Z}}_{j+1/2}};$

$D_{3,j}=\frac{{2w}_{3,j}{r}_{\mathrm{ri},j}{\mathrm{\&Delta;Z}}_j\mathrm{\&Delta;t}}{r_{\mathrm{ro},j}^2-r_{\mathrm{ri},j}^2};$ $D_{w,j}=\frac{{2w}_{w,j}{r}_{\mathrm{ro},j}{\mathrm{\&Delta;Z}}_j\mathrm{\&Delta;t}}{r_{\mathrm{ro},j}^2-r_{\mathrm{ri},j}^2};$ $D_{r,j}^0={\mathrm{\&rho;}}_r{c}_{\mathrm{pr}}{\mathrm{\&Delta;Z}}_j$

海水段隔水管绝热层与海水界面 $\left({r=r}_{\mathrm{ins}}\right):-k_{\mathrm{ins}}\frac{(T_{i\max ,j}^{n+1}-T_{i\max -1,j}^{n+1})}{r_{i\max ,j}\ln (r_{i\max ,j}/r_{i\max -1,j})}=w_{w,j}(T_{i\max ,j}^{n+1}-T_{w,j}^{\&infin;})$

5地层换热模型： 

地层内节点(i，j)及控制体如图3所示。其与上下节点的控制体界面为U和D，与内外节点的控制体界面为W和E。采用交替方向法计算地层内节点温度，离散格式温度控制方程如下： 

(1)Z方向隐式，r方向显式 

$A_{i,j-1}{T}_{i,j-1}^{n+1/2}+B_{i,j}{T}_{i,j}^{n+1/2}+C_{i,j+1}{T}_{i,j+1}^{n+1/2}=D_{i+1,j}{T}_{i+1,j}^n+D_{i-1,j}{T}_{i-1,j}^n+D_{i,j}^0{T}_{i,j}^n---\left(5\right)$

式中： 

$B_{i,j}=D_{i-1,j}+D_{i+1,j}+A_{i,j-1}-C_{i,j+1}+D_{i,j}^0$

$D_{i-1,j}=\frac{k_W{\mathrm{\&Delta;Z}}_j\mathrm{\&Delta;t}}{\ln (r_{i,j}/r_{i-1,j})}$

$D_{i+1,j}=\frac{k_E{\mathrm{\&Delta;Z}}_j\mathrm{\&Delta;t}}{\ln (r_{i+1,j}/r_{i,j})}$

$A_{i,j-1}=\frac{{-k}_{i,j-1/2}(r_E+r_W){\mathrm{\&Delta;r}}_i\mathrm{\&Delta;t}}{2{\mathrm{\&Delta;Z}}_{j-1/2}}$

$C_{i,j+1}=\frac{{-k}_{i,j+1/2}(r_E+r_W){\mathrm{\&Delta;r}}_i\mathrm{\&Delta;t}}{2\mathrm{\&Delta;}{Z}_{j+1/2}}$

$D_{i,j}^0={\mathrm{\&rho;}}_{i,j}{c}_{i,j}{\mathrm{\&Delta;r}}_i\frac{(r_E+r_W)}{2}{\mathrm{\&Delta;Z}}_j$

(2)r方向隐式，Z方向显式 

$A_{i-1,j}{T}_{i-1,j}^{n+1}{B}_{i,j}{T}_{i,j}^{n+1}{C}_{i+1,j}{T}_{i+1,j}^{n+1}=D_{i,j+1}{T}_{i,j+1}^{n+1/2}+D_{i,j-1}{T}_{i,j-1}^{n+1/2}+D_{i,j}^0{T}_{i,j}^n---\left(6\right)$

式中： 

$B_{i,j}=-A_{i-1,j}+C_{i+1,j}+D_{i,j-1}+D_{i,j+1}+D_{i,j}^0$

$A_{i-1,j}=\frac{{-k}_W{\mathrm{\&Delta;Z}}_j\mathrm{\&Delta;t}}{\ln (r_{i,j}/r_{i-1,j})}$

$C_{i+1,j}=\frac{{-k}_E{\mathrm{\&Delta;Z}}_j\mathrm{\&Delta;t}}{\ln (r_{i+1,j}/r_{i,j})}$

$D_{i,j-1}=\frac{k_{i,j-1/2}(r_E+r_W){\mathrm{\&Delta;r}}_i\mathrm{\&Delta;t}}{2\mathrm{\&Delta;}{Z}_{j-1/2}}$

$D_{i,j+1}=\frac{k_{i,j+1/2}(r_E+r_W){\mathrm{\&Delta;r}}_i\mathrm{\&Delta;t}}{{2\mathrm{\&Delta;Z}}_{j+1/2}}$

$D_{i,j}^0={\mathrm{\&rho;}}_{i,j}{c}_{i,j}{\mathrm{\&Delta;r}}_i\frac{(r_E+r_W)}{2}{\mathrm{\&Delta;Z}}_j$

上述离散方程中，W界面处： $\frac{\&PartialD;T}{\&PartialD;r}=\frac{(T_{i,j}^{n+1}-T_{i-1,j}^{n+1})}{r_W\ln (r_{i,j}/r_{i-1,j})};$ E界面处： $\frac{\&PartialD;T}{\&PartialD;r}=\frac{(T_{i+1,j}^{n+1}-T_{i,j}^{n+1})}{r_E\ln (r_{i+1,j}/r_{i,j})};$

地层段井壁界面处节点(i＝4)： $-k_e\frac{(T_{5,j}^{n+1}-T_{4,j}^{n+1})}{r_{4,j}\ln (r_{5,j}/r_{4,j})}=w_{3,j}(T_{3,j}^{n+1}-T_{4,j}^{n+1}).$

6界面传热模型 

界面主要包括海水段井筒与海水交界面(即隔水管绝热层与海水交界面)和井筒与地层交界面。主要的模块即：井壁与地层界面换热模块TInterface和隔水管绝热层内节点温度模块Tins。 

海水段井壁界面处节点(i＝4)： $-k_{\mathrm{ms}}\frac{(T_{5,j}^{n+1}-T_{4,j}^{n+1})}{r_{4,j}\ln (r_{5,j}/r_{4,j})}=w_{3,j}(T_{3,j}^{n+1}-T_{4,j}^{n+1})$

海水段环空入口段井壁界面处节点(i＝4)： 

隔水管外海底节点(i，J_b)： $w_w^{\&prime;}(T_{w,J_b}^{\&infin;}-T_{i,J_b}^{n+1})=k_{i,J_b}\frac{T_{l,J_b}^{n+1}-T_{i,J_{b+1}}^{n+1}}{Z_{i,J_{b+1}}-Z_{i,J_b}}$

隔水管内海底节点(i，J_b)： $w_b(T_{3,J_b-1}^{n+1}-T_{i,J_b}^{n+1})=k_{i,J_b}\frac{T_{i,J_b}^{n+1}-T_{i,J_b+1}^{n+1}}{Z_{i,J_b+1}-Z_{i,J_b}}$

7钻井液热物性计算 

根据节点处钻井液的温度和压力计算节点处钻井液密度、钻井液表观粘度、钻井液比热、钻井液热导率。采用Karstad模型计算钻井液的高温高压密度： 

ρ_m＝ρ_m0e^Γ(P，T)

Γ(P，T)＝γ_P(P-P₀)+γ_PP(P-P₀)²+            (3) 

γ_T(T-T₀)+γ_TT(T-T₀)²+γ_PT(P-P₀)(T-T₀) 

公式中γ_p、γ_pp、γ_T、γ_TT和γ_pT为拟合系数。 

钻井液比热和热导率预测模型如下： 

c_p(o，T)＝c₀+c_co+c_T ^T

k(o，T)＝k₀-k_co+k_T ^T

式中c₀、c_c、c_T为钻井液定压比热拟合参数，k₀、k_c、k_T为热导率拟合参数。 

8热源计算 

模拟器考虑钻井液摩阻压降生热、钻头破岩摩擦生热、钻柱与井壁摩擦生热计算钻柱内和环空这两个换热区域的内热源。 

摩阻压降生热： ${\stackrel{\&CenterDot;}{q}}_f=Q\&times;{\mathrm{Ip}}_f$

钻头破岩生热： ${\stackrel{\&CenterDot;}{q}}_b=\mathrm{\&alpha;}\&times;2\mathrm{\&pi;M}\&times;\frac{\mathrm{RPM}}{60}$

钻柱与井壁摩擦生热： ${\stackrel{\&CenterDot;}{q}}_{i,j}=\mathrm{\&beta;}\&times;{\mathrm{IZ}}_j\&times;2\mathrm{\&pi;M}\&times;\frac{\frac{\mathrm{RPM}}{60}}{D}$

网格划分结束后，按照下述步骤求解海洋半潜式平台钻井井筒循环温度场。 

1)应用初始条件 

海水区(泥线以上)所有节点的初始温度为对应深度处海水温度；地层区(泥线以下)所有节点的初始温度为对应深度处地层原始温度。 

2)根据钻井液总模拟循环时间和时间步长划分时间段，从初始时刻开始，累加时间，每一时间步的温度计算需要以上一时间步的温度场数据作为初始条件，迭代计算，直至收敛。每一时间步的温度场计算步骤如下： 

A用二维数组Node记录上一时间步结束时的温度数据； 

B将节点温度数据复制一份保存至二维数组NodeNew中； 

C收敛状态记为“假”； 

D以Node中的数据作为初始值按照钻柱内钻井液、钻柱管体、钻柱外环空、井筒界面、隔水管、地层的顺序按列依次计算各传热对象内所有节点的新时间步的温度数据，每一列节点的每个节点计算结束时都要比较新计算出的节点温度与NodeNew中保存的温度数据的差值，并记录在变量TDiff二维数组中，然后将新计算出的温度数据存入二维数组NodeNew中。 

E找出数组TDiff中的最大值，若该值很小，则认为该时间步温度计算达到收敛，结束该时间步的迭代，将NodeNew中的温度数据复制到Node数组中，进入下一时间步；否则，重复步骤D直至收敛。计算流程如图4所示。 

模拟器动态链接库的入口函数为WHTSumbersible(char FileInput，char FileCir，char FileMonitor。其中FileInput为温度场分析输入的数据文件名，FileCir为温度场分析结果输出文件名，模拟器将温度场分析结果存储到FileCir指定的文本文件中，FileMonitor为钻井液循环期间监视位置处不同时刻的温度数据文件名。 

数据输入：模拟器从FileInput指定的文件读入基础数据，数据文件组织见表1。 

表1输入文件数据结构 

模拟器根据输入数据计算分析结束后，将井筒内与近井地层区域的节点深度和温度场数据输出到FileCir指定的文本文件中，将监视位置的温度变化数据输出到FileMonitor指定的文件中。温度场输出文件数据结构见表2，监视位置温度输出文件数据结构见表3。 

表2温度场文件数据结构 

表3监视位置数据文件结构 

用户可根据WHTSubmersible模拟器输出的这两个数据文件，利用连接的接口程序或者商业数据处理程序进一步处理数据和生成报表。 

WHTSubmersible模拟器由19个主要功能模块组成，模块结构如图5所示，各模块功能如下。 

DataInput函数是WHT程序的数据输入函数，完成模拟井所有数据的输入； 

TPField函数是WHT程序的总功能模块，完成深水井筒温度场计算及数据存储功能； 

GridGeneration函数根据模拟井的井身结构、钻具结构、海水深度对温度场求解域进行轴向和径向网格划分，存储网格节点的轴向和径向几何信息及介质信息； 

TOriginGeneration函数根据海水和地层的垂直温度分布数据插值产生节点轴向深度处的原始温度；模拟器有两种方式建立原始温度剖面：根据海水和地层的温度分布关系式；根据实测海水和地层温度数据建立温度剖面。海洋油井海水段存在温跃层，地温也存在多个梯度。 

Ini函数对求解域内网格节点变量应用初始条件，赋初值； 

TInDrillStem函数、TDrillStem函数、TAnnu函数、TInterface函数、TRiser函数和TFormation函数分别用于计算钻柱内钻井液、钻柱管体、环空钻井液、环空外壁处、隔水管绝热层内和地层内所有节点的温度； 

ThermalPhysics函数。该模块根据节点处钻井液的温度和压力计算节点处钻井液密度、钻井液表观粘度、钻井液比热、钻井液热导率和高温高压密度； 

HeatGeneration函数。该模块考虑钻井液摩阻压降生热、钻头破岩摩擦生热、钻柱与井壁摩擦生热计算钻柱内和环空这两个换热区域的内热源； 

HPipe函数、HAnnu函数、HAcross函数的功能分别是计算管流、环空流与内外壁两个的、液体横掠圆管的强迫对流换热系数； 

ViscTM函数用于计算节点温度条件下钻井液的表观粘度； 

Kbound函数可用于计算轴向或径向相邻节点控制体界面处的介质热导； 

TriDiag函数用来求解具有三对角矩阵形式的代数方程组。 

采用Visual C++动态链接库技术开发半潜式平台循环温度场计算模拟器。 

有效性分析 

现场数据如下：水深63m，海面温度为28℃，排量为30l/s，钻井液密度为1.80g/cm3，定压比热为2.00kJ/kg/K，热导率为，钻具组合为：8-3/8″PDC钻头+12根6-1/2″钻铤+15根5″加重钻杆+5″钻杆，井身结构见表4，地温分三段描述：泥线至1500m、1500m至2100m和2100m至3000m，地温梯度依次为4.82℃/100m、2.41℃/100m和4.43℃/100m。 

表4井身结构 

利用该模拟器计算了南海某半潜式平台井筒温度，计算结果与实际结果基本吻合(理论 值与实测值对比见表5)，误差不超过5％。 

表5理论与实测出口温度对比 

以上所述，仅为本发明较佳的具体实施方式，本发明的保护范围不限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可显而易见地得到的技术方案的简单变化或等效替换均落入本发明的保护范围内。 

Claims (3)

1.一种半潜式平台瞬态钻井井筒温度计算方法，其特征在于，包括以下步骤：

1)应用初始条件

海水区所有节点的初始温度为对应深度处海水温度；地层区所有节点的初始温度为对应深度处地层原始温度；

2)根据钻井液总模拟循环时间和时间步长划分时间段，从初始时刻开始，累加时间，每一时间步的温度计算需要以上一时间步的温度场数据作为初始条件，迭代计算，直至收敛。

2.根据权利要求1所述的半潜式平台瞬态钻井井筒温度计算方法，其特征在于，所述步骤2)的温度场计算步骤如下：

A用二维数组Node记录上一时间步结束时的温度数据；

B将节点温度数据复制一份保存至二维数组NodeNew中；

C收敛状态记为“假”；

D以Node中的数据作为初始值按照钻柱内钻井液、钻柱管体、钻柱外环空、井筒界面、隔水管、地层的顺序按列依次计算各传热对象内所有节点的新时间步的温度数据，每一列节点的每个节点计算结束时都要比较新计算出的节点温度与NodeNew中保存的温度数据的差值，并记录在变量TDiff二维数组中，然后将新计算出的温度数据存入二维数组NodeNew中；

E找出数组TDiff中的最大值，若该值很小，则认为该时间步温度计算达到收敛，结束该时间步的迭代，将NodeNew中的温度数据复制到Node数组中，进入下一时间步；否则，重复步骤D直至收敛。

3.一种实现权利要求1所述方法的WHTSubmersible模拟器，其特征在于，包括以下19个模块：

DataInput函数是WHTSubmersible程序的数据输入函数，完成模拟井所有数据的输入；

TPField函数是WHTSubmersible程序的总功能模块，完成深水井筒温度场计算及数据存储功能；

GridGeneration函数根据模拟井的井身结构、钻具结构、海水深度对温度场求解域进行轴向和径向网格划分，存储网格节点的轴向和径向几何信息及介质信息；

TOriginGeneration函数根据海水和地层的垂直温度分布数据插值产生节点轴向深度处的原始温度；模拟器有两种方式建立原始温度剖面：根据海水和地层的温度分布关系式；根据实测海水和地层温度数据建立温度剖面，海洋油井海水段存在温跃层，地温也存在多个梯度；

Ini函数对求解域内网格节点变量应用初始条件，赋初值；

TInDrillStem函数、TDrillStem函数、TAnnu函数、TInterface函数、TRiser函数和TFormation函数分别用于计算钻柱内钻井液、钻柱管体、环空钻井液、环空外壁处、隔水管和地层内所有节点的温度；

HPipe函数、HAnnu函数、HAcross函数的功能分别是计算管流、环空流与内外壁两个的、液体横掠圆管的强迫对流换热系数；

ViscTM函数用于计算节点温度条件下钻井液的表观粘度；

Kbound函数可用于计算轴向或径向相邻节点控制体界面处的介质热导；

TriDiag函数用来求解具有三对角矩阵形式的代数方程组。

Images