Surrounding Rock Influence Correction for Array Laterolog Responses with Borehole Eccentricities in Horizontal Wells
-
摘要: 为消除围岩对阵列侧向测井仪器响应的影响,使测量结果更为准确,应对其进行校正。为此,建立了水平井三维层状介质模型,应用三维有限元算法研究了不同井眼偏心距下阵列侧向测井的正演响应特性,得到了数值模拟围岩校正图版;对数值模拟结果进行了非线性分析,得到了围岩快速校正公式和校正图版。通过数值模拟发现,阵列侧向测井仪器的探测深度越深,受围岩的影响越大;当井眼偏心程度系数小于20%时,井眼偏心距变化对阵列侧向测井响应的影响可以忽略;当井眼偏心程度系数大于40%后,阵列侧向测井响应受井眼偏心距的影响明显,视电阻率变化幅度最大可达44.8%。对比分析发现,快速校正图版与数值模拟图版具有较好的一致性,采样点最大相对误差小于5.3%。综合研究表明,快速校正图版可用于水平井中阵列侧向测井响应的实时校正。Abstract: The response of the array lateral logging tool is influenced by the surrounding rock.Therefore,a surrounding rock correction needs to be performed.In this paper,a 3D layered medium model of horizontal well was established as the first step.The forward response characteristics of array laterolog under various borehole eccentricities were studied by using the 3D finite element method,and the correction plate of surrounding rock was obtained.The nonlinear analysis of the numerical simulation plate was carried out,and the fast correction formula and correction plate of the surrounding rock/layer thickness were obtained.Research showed that the deeper the instrument detection depth,the more they are affected by surrounding rock.When the borehole eccentricity was lower than 20%,borehole eccentricity changes had little effect on the response of array lateral logging.However,when the borehole eccentricity was greater than 40%,the response of array laterolog logging was significantly affected by borehole eccentricity,and the maximum change range of apparent resistivity could reach 44.8%.It was found that the fast correction plate was in good agreement with the numerical simulation plate and the relative error of sampling points between the two plates was lower than 5.3%.As shown by the comprehensive study,the fast correction plate could be used for fast and real-time correction of array lateral well response in horizontal wells.
-
Keywords:
- horizontal well /
- borehole /
- eccentricity of borehole /
- array laterolog logging /
- surrounding rock /
- correction
-
超深地层温度、压力高,导致超深井钻井过程中井筒温度高,这对钻井液、钻井工具和测量仪器等性能的影响很大。因此,需要采取相应的技术方法控制井筒温度,在保证钻井安全的同时,提高钻井效率[1–8]。准确预测超深井井筒温度演变规律是实现井筒温度高效调控的基础。目前国内外学者主要基于稳态法和瞬态法预测和分析钻井过程中的井筒温度分布特征。稳态方法假设井筒内流体温度为稳态分布,并忽略热源项对井筒温度分布的影响:R. Spindler和A. R. Hasan等人[9–10]基于井筒与地层间热力学平衡原理,构建了井筒传热解析模型;A. Q. Al Saedi等人[11]考虑钻井液流动摩擦生热和钻柱旋转做功生热,建立了井筒温度场解析模型;Wu Xingru、Xu Boyue、Peng Yu和赵金洲等人[12–15]针对关井、生产和压裂过程中的井筒传热特点,分别提出了适用于计算关井、生产和压裂时井筒温度剖面的半解析模型。瞬态方法认为井筒流体温度随循环时间变化,Yang Mou等人[16]考虑径向温度梯度和轴向热传导的影响,构建了地层–井筒瞬态传热模型;Mao Liangjie等人[17]考虑钻头破岩生热和钻井液热物性变化的影响构建了井筒温度瞬态预测模型;李军和杨谋等人[18–19]针对深水钻井和高温定向井特点建立了地层–井筒瞬态传热模型,并分析了溢流、漏失等对井筒温度场的影响规律。与稳态法相比,瞬态方法预测井筒温度场精度更高[20]。由于钻井液性能受温度压力影响显著,导致井筒温度和压力也会相互影响[21]。近年来,国内外学者建立了井筒温度–压力耦合模型来精确计算井筒温度场和压力场。Chen Xin等人[22]针对欠平衡钻井过程,建立了考虑循环流体热物性的井筒温度场和压力场全瞬态耦合预测模型。Xu Zhengming和孙宝江等人[23–25]考虑气侵、压裂和固井等情况,建立了井筒温压耦合模型,准确刻画井筒温度压力的演变规律。
井筒降温方法主要可以分为2种:一种是地面降温,另一种是井下降温[26]。地面降温是采用较多的井筒降温方法,在地面将钻井液温度降低后再进入井筒循环,从而达到降低井筒温度的目的[27–31]。该方法在现场应用中取得了较好的效果,但是需要进行大量的试验,并专门研发相应的降温设备,降温设备会占用大量地面空间,且降温幅度小,无法精确控制井筒温度。井下降温方法可分循环降温、井下节流降温和绝缘降温。循环降温的工艺相对简单,其原理是在钻井的同时控制井筒温度,每钻进一段后停钻,观察井筒温度,如果温度合适,则继续钻进;如果温度过高,则钻井液循环一段时间,待温度合适后再继续钻进[32–33]。该方法需要多次测量井筒温度,且深井、超深井对循环降温的限制条件较多,操作繁琐。井下节流降温,本质上是一种生产气井防止生成天然气水合物的方法,其通过建立气液两相嘴流预测模型、节流温降预测模型和气井井筒压力温度分布预测模型等,增强高含水气井临界携液能力,避免生成天然气水合物[34–36]。该方法在中间环节上具有井筒降温作用,具体降温效果还有待研究。绝缘降温主要通过在钻杆上喷涂隔热涂层并调整涂层厚度、长度、导热系数等参数,阻止热量由地层传入井筒造成井筒温度升高。国外采用该方法防止在井下生成天然气水合物[37–41];国内则采用该方法保护随钻测控仪器[42–44],保护随钻测控仪器在合适的温度下正常工作。
钻杆喷涂隔热涂层是一种适合井下控温的方法,在大幅度精确降温的同时,操作方法相比于其他方法更加简便,施工成本更低。考虑经济性与喷涂工艺等因素,隔热涂层覆盖钻杆的长度不宜过长,其厚度也不宜过厚。因此,应用绝缘降温技术时,需研究隔热涂层喷涂位置、喷涂长度、喷涂厚度等因素对井筒降温效果的影响。笔者建立了井筒温度场瞬态模型,依据模型分析钻杆隔热涂层参数对井筒温度的影响,定性得出井筒温度与钻杆隔热涂层的关系,以期对万米超深井钻井井筒温度场控制起到指导作用。
1. 超深井井筒温度场瞬态模型的建立
超深井井筒温度场控制单元主要包括钻杆内钻井液、钻杆、环空内钻井液、套管、水泥环和地层,基于钻井液循环过程以及井筒与地层间的传热特点,作如下假设:
1)忽略钻井液循环过程中热传导传递的热量,仅考虑其轴向上钻井液流动携带的热量以及径向上与钻杆/井壁之间的对流换热;
2)对于钻杆,忽略其径向上热传导产生的热量交换,只考虑其在轴向上与钻柱内和环空内钻井液的对流换热;
3)钻井液为单相不可压缩流体,忽略温度、压力对其热物性和流变性的影响;
4)地层具有均质性,其热物性始终保持不变;
5)井眼和井壁形状规则,钻柱居中,全井段不存在钻柱屈曲和偏心现象;
6)距离井筒径向一定距离地层的温度始终为地层原始温度。
基于以上假设,采用热力学第一定律分别建立了钻杆内钻井液、环空内钻井液、套管、水泥环和地层的温度场控制方程。
1.1 基本公式
物体内能变化为:
Q=ρcVΔθ (1) 式中:Q为热量,J;ρ为物体密度,kg/m3;c为物体比热容,J/(kg·℃);V为物体体积,m3;Δθ为温度变化量,℃。
计算对流换热的牛顿冷却公式为:
Q=hAΔθ (2) 式中:h为对流换热系数,W/(m2·℃);A为换热表面面积,m2。
计算热传导的傅里叶定律公式为:
Q=−λAΔθδ (3) 式中:λ为导热系数,W/(m·℃);δ为导热距离,m。
热力学第一定律为:
dU=dQ+dW (4) 式中:dU为热力学能改变量,J;dQ为传热量,J;dW为外界做功,J。
1.2 地层–井筒传热控制方程
使用微元法,将井筒分为多个微元体。基于以上传热学基本公式,针对钻杆内钻井液、环空内钻井液、套管、水泥环和地层,分析其微元体的传热情况,建立微元体的传热控制方程。为降低模型的复杂度,提高计算效率,笔者采用综合传热系数描述钻杆内钻井液与钻杆壁或隔热涂层间的对流换热、隔热涂层与钻杆壁的导热以及环空内钻井液与钻杆壁间对流换热的复杂传热过程。
1)钻杆内钻井液传热模型。对于钻杆内部钻井液,其存在钻井液向下流动携带热量以及钻井液通过钻杆与环空钻井液间的换热。这2部分传热组成了钻杆内钻井液微元体内能的变化,其传热控制方程为:
ρdcp,dπr2p,i∂θd∂t=−ρdvdπr2p,icp,d∂θd∂z+2πrp,iK(θan−θd) (5) 式中:ρd为钻井液密度,kg/m3;cp,d为钻杆内钻井液的比热容,J/(kg·℃);rp,i为钻杆内半径,m;θd为钻杆内钻井液温度,℃;t为时间,s;vd为钻杆内钻井液流速,m/s;K为钻杆的综合传热系数,W/(m2·℃);θan为环空内钻井液的温度,℃。
2)环空内钻井液传热模型。对于环空内钻井液,其存在钻井液向上流动携带热量、通过钻杆壁与钻杆内钻井液换热以及与套管的对流换热。这3部分传热组成了环空内钻井液微元体内能的变化,其传热控制方程为:
ρancp,anπ(r2ca1,i−r2p,o)∂θan∂t=ρanvanπ(r2ca1,i−r2p,o)cp,an⋅∂θan∂z+2πrp,oK(θd−θan)+2πrca1,ihan-ca1(θca1−θan) (6) 式中:ρan为环空内钻井液密度,kg/m3;cp,an为环空内钻井液的比热容,J/(kg·℃);rca1,i为由内向外第一层套管内半径或井眼半径,m;rp,o为钻杆外半径,m;van为钻杆内钻井液流速,m/s;han-ca1为环空内钻井液与套管或井壁的对流换热系数,W/(m2·℃);θca1为由内向外第一层套管或井眼的温度,℃。
3)套管传热模型。对于套管,其存在轴向热传导、与环空内钻井液对流换热以及与水泥环及后方地层的热传导。这3部分传热组成了套管微元体内能的变化,其传热控制方程为:
ρca1cp,ca1π(r2ca1,o−r2ca1,i)∂θca1∂t=λca1π(r2ca1,o−r2ca1,i)⋅∂2θca1∂z2+2πrca1,ihan - ca1(θan−θca1)+2πrca1,o2λca1-ce1rca2,i−rca1,i(θce1−θca1) (7) 式中:ρca1为套管的密度,kg/m3;cp,ca1为套管的比热容,J/(kg·℃);rca1,o为由内向外第一层套管的外半径,m;λca1为套管的导热系数,W/(m·℃);rca2,i为由内向外第二层套管内半径,m;λca1-ce1为套管与水泥环的综合导热系数,W/(m·℃);θce1为由内向外第一层水泥环的温度,℃。
4)水泥环及外侧地层传热模型。对于水泥环及外侧地层,采用柱坐标建模,其存在径向和轴向热传导。这2部分传热组成了水泥环及外侧地层微元体内能的变化,其传热控制方程为:
ρfcp,f∂θf∂t=λfr∂∂r(r∂θf∂r)+λf∂2θf∂z2 (8) 式中:ρf为水泥环或地层的密度,kg/m3;cp,f为水泥环或地层的比热容,J/(kg·℃);θf为水泥环或地层温度,℃;r为微元体距离井筒中心轴线的距离,m。
1.3 综合传热系数
对于钻杆壁部分,隔热涂层的导热系数远小于钻杆的导热系数,采用传热热阻形式计算钻杆的综合传热系数,如图1所示。
钻井液与隔热钻杆的传热过程包含4个环节:1)环空钻井液到钻杆外壁面的热量传递;2)钻杆外壁面到钻杆内壁面(隔热涂层外壁面)的热量传递;3)隔热层外壁面到隔热涂层内壁面的热量传递;4)隔热涂层内壁面到钻杆内钻井液的热量传递。串联的每个环节的热流量Φ是相同的[45]。Φ分别表示为:
Φ=Ahan(θan−θs,1) (9) Φ=Aλ2δ2(θs,1−θs,2) (10) Φ=Aλ1δ1(θs,2−θs,3) (11) Φ=Ahd(θs,3−θd) (12) 式中:Φ为热流量,W;han为环空内钻井液与钻杆外壁的对流换热系数,W/(m2·℃);θs,1为钻杆外壁面的温度,℃;λ2为钻杆的导热系数,W/(m·℃);δ2为钻杆的壁厚,m;θs,2为钻杆内壁面(隔热涂层外壁面)的温度,℃;λ1为隔热涂层的导热系数,W/m/℃;δ1为隔热涂层的厚度,m;θs,3为隔热涂层内壁面的温度,℃;hd为钻杆内钻井液与隔热涂层内壁面的对流换热系数,W/(m2·℃)。
将式(9)—式(12)改写为温压形式,并消去壁面温度θs,1,θs,2和θs,3,可得:
Φ=A(θan−θd)1han+δ2λ2+δ1λ1+1hd=AK(θan−θd) (13) 则综合传热系数K可表示为:
K=11han+δ2λ2+δ1λ1+1hd (14) 1.4 初始条件与边界条件
1)初始条件。假设钻井液开始循环前,井筒已静置一段时间,井筒温度与地层原始温度相等,以此作为模型的初始条件。
{θd(z,t=0)=θf(z)=θs+Gfzθp(z,t=0)=θf(z)=θs+Gfzθan(z,t=0)=θf(z)=θs+Gfzθca(z,t=0)=θf(z)=θs+Gfz……θf(z,t=0)=θf(z)=θs+Gfz (15) 式中:θs为地表的温度,℃;Gf为地层的温度梯度,℃/m。
2)边界条件。开始循环后,井口实时测量的入口温度保持不变,则井口边界条件为:
θd(z=0,t)=θin (16) 式中:θin为钻柱内钻井液的入口温度,℃。
在井底处,钻柱内钻井液与环空内钻井液的温度相等,则井底边界条件可表示为:
θd(z=Zb,t)=θan(z=Zb,t) (17) 式中:Zb为井底的深度,m。
径向上与井筒一定距离处地层的温度始终为地层原始温度。
2. 模型求解及验证
2.1 模型求解
基于有限差分法分别对钻杆内钻井液传热控制方程、环空内钻井液传热控制方程、套管传热控制方程和水泥环及后续地层传热控制方程进行离散。其中,一阶偏导项采用向后差分格式,二阶偏导项采用中心差分格式。
钻杆内钻井液传热控制方程的离散形式为:
(ρdcp,dπr2p,iΔt+ρdvdπr2p,icp,dΔz+2πrp,iK)θntd,j−ρdvdπr2p,icp,dΔzθntd,j−1−2πrp,iKθntan,j=ρdcp,dπr2p,iΔtθnt−1d,j (18) 式中:nt为时间步;j为空间步。
环空内钻井液传热控制方程的离散形式为:
[ρancp,anπ(r2ca1,i−r2p,o)Δt+ρanvanπ(r2ca1,i−r2p,o)cp,anΔz+2πrp,oK+2πrca1,ihan - ca1]θntan,j−ρanvanπ(r2ca1,i−r2p,o)cp,anΔzθntan,j+1−2πrp,oKθntd,j−2πrca1,ihan - ca1θntca1,j=−ρancp,anπ(r2ca1,i−r2p,o)Δtθnt−1an,j (19) 套管传热控制方程的离散形式为:
[ρca1cp,ca1π(r2ca1,o−r2ca1,i)Δt+2λca1π(r2ca1,o−r2ca1,i)Δz2+2πrca1,ihan - ca1+2πrca1,oλca1 - ce1rca2,i−rca1,i2]θntca1,j−λca1π(r2ca1,o−r2ca1,i)Δz2θntca1,j−1−λca1π(r2ca1,o−r2ca1,i)Δz2θntca1,j+1−2πrca1,ihan - ca1θntan,j−2πrca1,oλca1 - ce1rca2,i−rca1,i2θntce1,j=ρca1cp,ca1π(r2ca1,o−r2ca1,i)Δtθnt−1ca1,j (20) 水泥环及外侧地层传热控制方程的离散形式为:
[ρicp,iΔt+λ(i−1)iriΔrilnriri−1+λi(i+1)riΔrilnri+1r+2λiΔz2]θnti,j−λiΔz2θnti,j−1−λiΔz2θnti,j+1−λ(i−1)iriΔrilnriri−1θnti−1,j−λi(i+1)riΔrilnri+1riθnti+1,j=ρiCp,iΔtθnt−1i,j (21) 为了方便计算,离散后的控制方程合并同类项,可得:
Akθnti,j+Bkθnti,j−1+Ckθnti,j+1+Dkθnti−1,j+Ekθnti+1,j=Fkθnt−1i,j (22) 式中:Ak,Bk,Ck,Dk,Ek和Fk为各项系数。
为便于高效求解,联立所有传热控制方程,利用各项系数构建稀疏矩阵:
(A1+B1E1⋯⋱⋱⋱⋯B1A1E1⋯⋱⋱⋱⋯B1A1E1⋯D2A2+B2C2E2⋯⋱⋱⋱⋱⋱⋯D2B2A2C2E2⋯⋱⋱⋱⋱⋱⋯D2B2A2+C2E2⋯D3A3C3E3⋯⋱⋱⋱⋱⋯D3A3C3E3⋯⋱⋱⋱⋱⋯D3A3+C3E3⋯Dca3,1Aca3,1+Bca3,1Cca3,1⋯⋱⋱⋱⋱⋯Dca3,1Bca3,1Aca3,1Cca3,1⋯⋱⋱⋱⋱⋯Dca3,1Bca3,1Aca3,1+Cca3,1⋯⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋱)(θnd,1⋮θnd,j⋮θnd,Nθndp,1⋮θndp,j⋮θndp,Nθna,1⋮θna,j⋮θna,Nθnca3,1,1⋮θnca3,1,j⋮θnca3,1,N⋮)=(F1θn−1d,1⋮F1θn−1d,j⋮F1θn−1d,NF2θn−1dp,1⋮F2θn−1dp,j⋮F2θn−1dp,NF3θn−1a,1⋮F3θn−1a,j⋮F3θn−1a,NFca3,1θnca3,1,1⋮Fca3,1θnca3,1,j⋮Fca3,1θnca3,1,N⋮) (23) 该行列式是五对角行列式,其主对角线占优,说明迭代收敛,可以使用高斯–塞德尔迭代对各传热控制方程进行求解。
2.2 模型验证
为验证本文所建模型的正确性,将其与Holmes & Swift模型[11]进行了对比。Holmes & Swift模型是稳态解析模型,未考虑钻杆壁厚及井筒温度对周边地层的影响。为了能够对比这2个模型,需要调整本文模型,使循环时间足够大以达到稳态温度场,并忽略隔热涂层、钻杆壁厚和地层的影响。将以下参数输入Holmes & Swift模型和本文所建模型,井深4 572 m,井眼直径212.7 mm;地层比热容873.36 J/(kg·℃),地层密度2644.2 kg/m3,地温梯度2.31 ℃/100m,地层导热系数2.249 W/(m·℃);地表温度15.278 ℃;钻柱外径168.3 mm,钻柱内径151.5 mm,钻柱材料密度7 800 kg/m3,钻柱壁导热系数43.75 W/(m·℃),钻柱壁比热容400 J/(kg·℃);钻井液密度1 198.3 kg/m3,钻井液比热容1 674.7 J/(kg·℃),钻井液入口温度23.889 ℃;循环排量0.013 2 m3/s;钻井液与钻柱壁对流换热系数340.695 8 W/(m2·℃),钻井液与套管对流换热系数5.678 3 W/(m2·℃),计算得到钻柱内和环空内钻井液温度分布(见图2)。
由图2可知,本文模型计算出的钻柱内和环空内钻井液的温度分布曲线与Holmes & Swift模型计算出的温度分布曲线总体吻合程度高,验证了本文所建模型的正确性。
2.3 实例验证
为了进一步验证本文所建立井筒–地层瞬态传热模型在超深定向井实际井身结构条件下的精度和可靠性,利用塔里木盆地1口超深定向井及其井底实测循环温度对模型进行了验证。
X井为塔里木盆地的一口超深探井,主要钻探目的是探索奥陶系中统及中—下统储层发育情况与含油气性。该井完钻井深8 543 m,垂深8 296 m,补心高13.70 m,钻进8 425~8 516 m井段时,主要钻遇灰岩地层,正常循环钻进22 h后,随钻测量仪器测得井底循环温度为160 ℃。该井的实际井身结构和井眼扩大率情况见表1。
表 1 X井的井身结构与井眼扩大率Table 1. Casing program and borehole expansion rate of Well X开次 钻头直径/mm 井深/m 套管外径/mm 套管下深/m 水泥返深/m 井眼扩大率,% 导管 660.4 105.00 508.0 105.00 0 一开 444.5 1 199.50 339.7 1199.15 0 二开 311.1 4 363.00 250.8 4 362.00 0 3.94 三开 215.9 7 728.00 177.8 7 726.77 4 162 4.14 四开 149.2 8 543.00 3.43 该井段钻进所采用的钻具组合为ϕ149.2 mm钻头+ϕ120.0 mm螺杆+ϕ127.0 mm浮阀+ϕ120.6 mm无磁钻铤+ϕ121.0 mm无磁悬挂短节+ϕ88.9 mm钻杆+ϕ88.9 mm加重钻杆+ϕ88.9 mm钻杆+ϕ127.0 mm钻杆+ϕ139.7 mm短钻杆。钻井液、钻柱、套管、水泥环和地层等传热介质的热物性参数见表2。
表 2 X井的传热介质热物性参数Table 2. Thermal property parameters of heat transfer medium of Well X介质 密度/
(kg·m−3)比热容/
(J·kg−1·℃−1)导热系数/
(W·m−1·℃−1)钻井液 1 600 1 600 1.200 钻柱 7 800 500 48.000 套管 7 800 500 48.000 水泥环 2 140 2 000 0.700 地层岩石 2 655 985 2.021 8 425~8 516 m井段时采用聚磺钻井液钻进,黏度22 mPa·s,排量15 L/s,钻井液入口温度34 ℃,出口温度35 ℃,考虑实际井身结构、钻具组合和平均井径扩大率的影响,以完钻静置一段时间后的电测温度为地层原始温度,构建X井井筒循环温度场预测模型,预测该井井筒循环温度场,结果如图3所示。
由图3可知:模型预测X井井底循环温度为163.82 ℃,出口温度为33.68 ℃;与实测结果相比,预测井底循环温度的绝对误差为−3.82 ℃,相对误差为2.39%,预测出口温度的绝对误差为1.32 ℃,相对误差为3.77%;预测结果与实测结果较为吻合,二者相对误差在5%以内,满足工程上相对误差不大于10%的要求,进一步验证了上文所建立井筒–地层瞬态温度场预测模型的有效性和可靠性。
3. 钻杆内隔热涂层对井筒温度场的影响
采用1口井深为8 000 m的模拟直井,分析钻杆内隔热涂层导热系数、厚度和长度对井筒温度场的影响规律。模拟井的参数:井深8 000 m;井眼直径0.213 m;地层比热容837 J/(kg·℃),地层密度2 645 kg/m3,地温梯度0.023 1℃/m,地表温度15.3 ℃,地层导热系数2.25 W/(m·℃);钻柱外径0.168 m,钻柱内径0.151 m,钻柱材料密度7 800 kg/m3,钻柱导热系数43.75 W/(m·℃),钻柱比热容400 J/(kg·℃);钻井液密度1 200 kg/m3,钻井液黏度45.4 mPa·s,钻井液比热容1 675 J/(kg·℃),钻井液入口温度24.0 ℃,钻井液导热系数1.73 W/(m·℃);循环排量13.2 L/s。
3.1 隔热涂层导热系数
钻井参数与钻井液性能保持不变,钻井液循环时间为100 h,钻杆内隔热涂层的厚度为2.0 mm,8 000 m钻柱均喷涂内隔热涂层,模拟计算喷涂不同导热系数隔热涂层下的井筒环空温度剖面,结果如图4所示。
由图4可知:随着钻杆内隔热涂层导热系数增大,井筒环空温度逐渐升高,且导热系数越大,对井筒环空温度的影响越小;此外,环空温度最高点的位置逐渐变深,而出口温度则随着钻杆内隔热涂层导热系数增大而降低。在相同钻井参数和循环时间条件下,钻杆内隔热涂层的导热系数越小,钻杆对其内钻井液的保温效果越好,则井底处的钻井液温度越低;钻井液上返过程中,由于在井底流入环空时的温度较低,所以在吸收相同热量条件下,采用隔热钻杆时的环空钻井液温度低于无隔热钻杆;在井口附近,由于隔热钻杆的保温作用,钻杆内钻井液从环空钻井液中吸收的热量减少,热交换作用减弱,且钻杆隔热涂层的导热系数越小,环空出口温度越高。当全部钻柱喷涂厚度为2.0 mm、导热系数为0.02 W/(m·℃)的隔热涂层时,相比于地层静态温度,井底循环温度降低了148 ℃,降低率高达74%;相比于无隔热涂层,井底循环温度降低了105 ℃,降低率为67%。
3.2 隔热涂层厚度
钻井参数和钻井液性能保持不变,钻井液循环时间为100 h,钻杆内隔热涂层的导热系数为0.06 W/m/℃,8 000 m钻柱内均喷涂隔热涂层,模拟计算钻柱喷涂不同厚度隔热涂层时的井筒环空温度剖面,结果如图5所示。
由图5可知:随着钻杆内隔热涂层厚度增大,井筒环空温度逐渐降低,且厚度增加到一定程度后,对井筒环空温度的影响程度降低;此外,出口温度随着钻杆内隔热涂层厚度增大而升高,而环空温度最高点的位置则随着隔热涂层厚度增大而不断变浅。在相同钻井参数、循环时间及隔热涂层参数条件下,钻杆内隔热涂层的厚度越大,钻杆的综合传热系数越小,隔热钻杆对钻杆内钻井液的保温效果越好,井底处的钻井液温度越低;在钻井液上返过程中,由于井底钻井液流入环空时的温度较低,所以在吸收相同热量的条件下,钻杆内隔热涂层厚度大时环空钻井液的温度低于无隔热钻杆;在井口附近,由于隔热钻杆的保温作用,隔热涂层厚度越大,钻杆内钻井液从环空钻井液中吸收的热量越少,热交换作用越弱,致使出口温度随着隔热涂层厚度增厚而升高。井下钻柱内喷涂厚度为4.0 mm、导热系数为0.06 W/(m·℃)的隔热涂层时,相比于地层静态温度,井底循环温度降低了138 ℃,降低率达69%;相比于钻杆内无隔热涂层,井底循环温度降低了95 ℃,降低率达60%。
3.3 隔热涂层长度
钻井参数与钻井液性能保持不变,钻杆内隔热涂层的导热系数为0.06 W/(m·℃)、隔热涂层喷涂厚度为2.0 mm,钻井液循环时间为100 h,模拟计算钻柱内喷涂不同长度隔热涂层时的井筒环空温度剖面,结果如图6所示。
由图6可知:随着钻杆内隔热涂层长度增加,井筒环空温度逐渐降低,隔热涂层越长,井底循环温度越低;此外,环空温度最高点的位置则随着隔热涂层长度增大而不断变浅,而出口温度则随着隔热涂层长度增加而升高,当井口钻杆有隔热涂层时,出口温度升高最显著。在相同钻井参数、循环时间、钻杆内隔热涂层厚度和导热系数条件下,隔热涂层长度越大,钻杆对其内钻井液的保温效果越好,井底循环温度越低;钻井液上返过程中,由于流入环空钻井液的温度低,所以在吸收相同热量的条件下,隔热涂层长度大的环空钻井液温度低于无隔热钻杆。在井口附近,当井口钻杆没有隔热涂层时,热交换作用显著,钻杆内钻井液的温度迅速升高;当井口钻杆具有隔热涂层时,由于隔热涂层的保温作用,钻杆内钻井液从环空钻井液中吸收的热量少,热交换作用弱,致使钻杆内温度升高幅度很小。当井内全部管柱喷涂厚度2.0 mm、导热系数0.06 W/(m·℃)的隔热涂层时,相比于地层静态温度,井底循环温度降低了116 ℃,降低率高达58%,相比于无隔热涂层的井筒循环温度,井底循环温度降低了74 ℃,降低率47%;而当仅井底800 m喷涂厚度2.0 mm、导热系数0.06 W/(m·℃)的隔热涂层时,相比于地层静态温度,井底循环温度仅降低了45 ℃,降低率23%,相比于无隔热涂层的井底循环温度,降低了3 ℃左右,降低率不到2%。
4. 结 论
1)针对超深井钻井开泵循环钻井液的过程,基于能量守恒原理,构建了考虑钻杆内喷涂隔热涂层时的超深井井筒–地层瞬态传热模型;与Holmes & Swift模型和超深井钻井现场实际数据进行了对比,验证了模型的可靠性和准确性。
2)钻杆内隔热涂层的导热系数对井底循环温度影响显著,随着钻杆内隔热涂层导热系数减小,井筒环空温度迅速降低,出口温度升高,环空最高温度点的位置上移。
3)钻杆内隔热涂层的厚度和长度对井筒温度有重要影响,井筒环空温度随着钻杆内隔热涂层厚度增大而降低,出口温度升高,环空最高温度点的位置上移。
4)建立的考虑隔热钻杆的超深井井筒–地层瞬态传热模型进一步拓展了传统的井筒传热模型,为超深井、特深井钻井过程中井筒温度的调控与隔热钻杆参数的优选提供了理论基础。
-
[1] 刘振华,胡启.阵列侧向测井响应的计算及其特征[J].西安石油学院学报(自然科学版),2002,17(1):53-57. LIU Zhenhua,HU Qi.Calculation and characteristics of array laterolog responses[J].Journal of Xi’an Petroleum Institute(Natural Science Edition),2002,17(1):53-57. [2] 仵杰,谢尉尉,解茜草,等.阵列侧向测井仪器的正演响应分析[J].西安石油大学学报(自然科学版),2008,23(1):73-76,80. WU Jie,XIE Weiwei,XIE Xicao,et al.Forward response analysis of array lateral logging tool[J].Journal of Xi’an Shiyou University(Natural Science Edition),2008,23(1):73-76, 80. [3] 潘克家,王文娟,汤井田,等.高分辨率阵列侧向测井的数学模型及有限元快速正演[J].地球物理学报,2013,56(9):3197-3211. PAN Kejia,WANG Wenjuan,TANG Jingtian,et al.Mathematical model and fast finite element modeling of high resolution array lateral logging[J].Chinese Journal of Geophysics,2013,56(9):3197-3211. [4] 范宜仁,蒋建亮,邓少贵,等.高分辨率阵列侧向测井响应数值模拟[J].测井技术,2009,33(4):333-336. FAN Yiren,JIANG Jianliang,DENG Shaogui,et al.Numerical simulation of high resolution array lateral logging responses[J].Well Logging Technology,2009,33(4):333-336. [5] 冯琳伟,汪德刚,贺飞,等.薄互层和倾斜地层阵列侧向测井响应计算分析[J].测井技术,2013,37(1):80-84. FENG Linwei,WANG Degang,HE Fei,et al.On response characteristics analysis of HAL tool in thin inter-beds and deviated hole formation[J].Well Logging Technology,2013,37(1):80-84. [6] 秦绪英.电阻率测井泥浆侵入影响校正方法研究:以鄂尔多斯盆地天然气储层为例[J].石油物探,2006,45(5):537-541. QIN Xuying.Slurry invasion correction in resistivity logging:case study of natural gas reservoir in Ordos Basin[J].Geophysical Prospecting for Petroleum,2006,45(5):537-541. [7] 邓少贵,李智强,范宜仁,等.斜井泥浆侵入仿真及其阵列侧向测井响应数值模拟[J].地球物理学报,2010,53(4):994-1000. DENG Shaogui,LI Zhiqiang,FAN Yiren,et al.Numerical simulation of mud invasion and its array laterolog response in deviated wells[J].Chinese Journal of Geophysics,2010,53(4):994-1000. [8] 徐建华,翟丽华,朱德怀,等.水平井中双侧向测井的偏心井眼校正[J].石油地球物理勘探,1994,29(4):520-525. XU Jianhua,ZHAI Lihua,ZHU Dehuai,et al.Eccentric borehole correction for dual laterolog in horizontal wells[J].Petroleum Geophysical Exploration,1994,29(4):520-525. [9] 邓少贵,李竹强,李智强.水平井双侧向测井响应及层厚/围岩影响快速校正[J].石油勘探与开发,2009,36(6):725-729. DENG Shaogui,LI Zhuqiang,LI Zhiqiang.Response of dual laterolog and fast correction for layer thickness and shoulder bed in horizontal wells[J].Petroleum Exploration and Development,2009,36(6):725-729. [10] 谭茂金,高杰,邹友龙,等.盐水泥浆条件下定向井双侧向测井环境校正方法研究[J].地球物理学报,2012,55(4):1422-1432. TAN Maojin,GAO Jie,ZOU Youlong,et al.Study on correction method of dual laterolog environment for directional wells under saltwater mud[J].Chinese Journal of Geophysics,2012,55(4):1422-1432. [11] 吴宝玉,夏宏泉,张智勇.随钻电阻率测井的围岩影响及校正方法研究[J].西南石油学院学报,2006,28(6):20-23. WU Baoyu,XIA Hongquan,ZHANG Zhiyong.Effect of shoulder beds on resistivity logging while drilling and its correction method[J].Journal of Southwest Petroleum Institute,2006,28(6):20-23. [12] 付恩玲,杨林,张淑芝,等.斜井中双感应测井曲线围岩校正方法[J].测井技术,2004,28(2):124-127,179. FU Enling,YANG Lin,ZHANG Shuzhi,et al.Correction method of double induction log curves in inclined shafts[J].Well Logging Technology,2004,28(2):124-127,179. [13] 于新娟,刘迪仁,蔡明,等.围岩对井间电磁测井的影响[J].石油地球物理勘探,2011,46(4):654-656. YU Xinjuan,LIU Diren,CAI Ming,et al.Influence of surrounding rock of crosshole electromagnetic logging[J].Oil Geophysical Prospecting,2011,46(4):654-656. [14] 舒心,柯式镇,许巍,等.高分辨率阵列侧向测井井眼影响自动校正研究[J].断块油气田,2016,23(4):470-475. SHU Xin,KE Shizhen,XU Wei,et al.Auto-correction of borehole effect for high-resolution array laterolog[J].Fault-Block Oil Gas Field,2016,23(4):470-475. [15] 倪小威,徐观佑,敖旋峰,等.阵列侧向测井曲线极化角影响因素研究[J].石油钻探技术,2018,46(2):120-126. NI Xiaowei,XU Guanyou,AO Xuanfeng,et al.The influencing factors on the polarizing angle of array laterolog curves[J].Petroleum Drilling Techniques,2018,46(2):120-126. [16] 刘迪仁,万文春,夏培,等.不同井眼条件下的水平井双侧向测井响应[J].物探与化探,2012,36(3):422-425. LIU Diren,WAN Wenchun,XIA Pei,et al.The response of dual laterolog to different borehole conditions in the horizontal well[J].Geophysical and Geochemical Exploration,2012,36(3):422-425. [17] 张庚骥.电法测井:上册[M].北京:石油工业出版社,1984:72-77. ZHANG Gengji.Electrical logging:Ⅰ[M].Beijing:Petroleum Industry Press,1984:72-77. [18] 夏培,刘迪仁,万文春,等.水平井各向异性地层双侧向测井响应数值模拟[J].石油天然气学报,2011,33(8):104-106. XIA Pei,LIU Diren,WAN Wenchun,et al.Numerical simulation of dual laterolog response in anisotropic formation of horizontal wells[J].Journal of Oil and Gas Technology,2011,33(8):104-106. [19] 倪小威,徐观佑,冯加明,等.斜井中各向异性储层阵列侧向测井正演响应特性[J].断块油气田,2017,24(5):637-641. NI Xiaowei,XU Guanyou,FENG Jiaming,et al.Forward response of array lateral logging in anisotropic reservoir of inclined shaft[J].Fault-Block Oil Gas Field,2017,24(5):637-641. [20] 陈光,任志良,孙海柱.最小二乘曲线拟合及Matlab实现[J].兵工自动化,2005,24(3):107-108. CHEN Guang,REN Zhiliang,SUN Haizhu.Curve fitting in least-square method and its realization with Matlab[J].Ordnance Industry Automation,2005,24(3):107-108. -
期刊类型引用(5)
1. 孔德政,罗姜民,肖冰清,姜翠苹,陈建国,李艳君. 基于反演技术的碎屑岩高精度三维观测系统研究. 非常规油气. 2024(06): 17-24 . 百度学术
2. 谢关宝,李永杰,吴海燕,黎有炎. 近井眼洞穴型地层双侧向测井敏感因素分析. 石油钻探技术. 2020(01): 120-126 . 本站查看
3. 张永健,汤子余,卢美月,吴承美. 基于曲线重构的复杂砂砾岩体岩性识别研究——以北16井区梧桐沟组为例. 西安石油大学学报(自然科学版). 2020(02): 20-25+34 . 百度学术
4. 林香亮,袁瑞,孙玉秋,王超,陈长胜. 支持向量机的基本理论和研究进展. 长江大学学报(自科版). 2018(17): 48-53+6 . 百度学术
5. 孔强夫,杨松. 塔河碎屑岩储层测井解释评价研究进展与挑战. 石油地质与工程. 2016(03): 81-83 . 百度学术
其他类型引用(6)
计量
- 文章访问数: 3999
- HTML全文浏览量: 93
- PDF下载量: 5403
- 被引次数: 11