D0I:0.16076j.cki.cjhd22.B.05 A辑第17卷第3期 水动力学研究与进展 Ser.A.Vol 17.No.3 2002年6月 JOURNAL OF HYDRODYNAMICS June 2002 文章编号：10004874(2002)03038209 多层地层中的井筒及地层 温度解析解 卢德唐，曾亿山，郭永存 (中国科学技术大学力学系，安徽合肥230026) 摘要：本文采用分层地层的假设米近似实际情况下的地层垂向的非均质性.并在此基础 上建立了多孔介质中的热传导问题的数学模型通过对方程的无量纲化及Laplace变换给出了 Laplace空间上的油气井在注入及生产情况下方程的解.最后将Laplace空间的解解析反演到实空 间上. 关键词：无量纲温度，多孔介质：热传导，laplace变换：解析反演 中图分类号：0357.5 文献标识码：A 1前言 当井筒中有流体注入或有流体采出时，周围地层与流体之间就存在温度差，使得流体与地 层的温度重新分布。在石油的勘探与开发中，了解井筒中流体温度随井深，时间、产量的变化 关系是很重要的。例如如果给出井筒温度与流量的关系，那么，就可以利用井筒温度来反求流 体产量：给出井筒温度与时间的关系，就可以利用井筒温度反求地层的热力学参数等。 50年代以来，就有一些学者研究多孔介质热传导问题，其中R在并筒传热方 面的研究最为经典，他引起了综合传热系数并给出了综合传热系数的表达式，但Ram心y的研 究采用了过多的假设，这使得Ramey的井筒瞬时传热导问题的解仅适合时间较大时的情况. 由于多孔介质中的热传导问题非常复杂，对井筒或地层传热问题研究最多的是数值解4习，因 为数值解可考虑许多复杂的问题（如地层的热力学参数的非均质性等，.但数值模拟以往往过于 复杂，也需要高性能的计算机一般人也很难掌握。在实际的应用中，解析解更利于人们对问 题本质的了解。本文正是从这一目的出发，根据热传导问题的性质，采用较符合实际的假设 (将地层分成多层且地层热力学参数在每个小层中为常数)，给出地层热传导方程及井筒中的 流体流动方程。对方程无量纲化后，给出无量纲方程的解。 代鑫公背管#础研觉9方项日90205，国家自然科学基金资助项目(401020m0,安 自然科 基金资围 日(200 作者充C德唐CI男教授博悬fronic Publishing House..All rights reserved http://www
A 辑第17 卷第 3 期 水 动 力 学 研 究 与 进 展 Ser.A, Vol.17 , No.3 2002 年 6 月 JOURNAL OF HYDRODYNAMICS June , 2002 文章编号:1000-4874(2002)03-0382-09 多层地层中的井筒及地层 温度解析解 卢德唐 , 曾亿山 , 郭永存 (中国科学技术大学力学系 ,安徽合肥 230026) 摘 要: 本文采用分层地层的假设来近似实际情况下的地层垂向的非均质性, 并在此基础 上建立了多孔介质中的热传导问题的数学模型, 通过对方程的无量纲化及 Laplace 变换, 给出了 Laplace 空间上的油气井在注入及生产情况下方程的解。 最后将 Laplace 空间的解解析反演到实空 间上。 关 键 词: 无量纲温度;多孔介质;热传导;Laplace 变换;解析反演 中图分类号: O357.5 文献标识码:A 1 前言 当井筒中有流体注入或有流体采出时,周围地层与流体之间就存在温度差 ,使得流体与地 层的温度重新分布。在石油的勘探与开发中 ,了解井筒中流体温度随井深 、时间 、产量的变化 关系是很重要的 。例如如果给出井筒温度与流量的关系,那么, 就可以利用井筒温度来反求流 体产量;给出井筒温度与时间的关系,就可以利用井筒温度反求地层的热力学参数等 。 50 年代以来 , ,就有一些学者研究多孔介质热传导问题 [ 1, 2] ,其中 Ramey [ 3] 在井筒传热方 面的研究最为经典, 他引起了综合传热系数, 并给出了综合传热系数的表达式 ,但 Ramey 的研 究采用了过多的假设 ,这使得 Ramey 的井筒瞬时传热导问题的解仅适合时间较大时的情况 。 由于多孔介质中的热传导问题非常复杂 ,对井筒或地层传热问题研究最多的是数值解[ 4, 5] ,因 为数值解可考虑许多复杂的问题 (如地层的热力学参数的非均质性等), 但数值模拟往往过于 复杂, 也需要高性能的计算机, 一般人也很难掌握。在实际的应用中 ,解析解更利于人们对问 题本质的了解 。本文正是从这一目的出发 ,根据热传导问题的性质 ,采用较符合实际的假设 (将地层分成多层,且地层热力学参数在每个小层中为常数),给出地层热传导方程及井筒中的 流体流动方程。对方程无量纲化后 ,给出无量纲方程的解。 收稿日期: 1998-11-21 基金项目: 国家重点基础研究“ 973” 项目(G1999032805);国家自然科学基金资助项目(10102020);安徽 省教育自然科学研究基金资助项目(2001KJ228) 作者简介: 卢德唐(1966～ ), 男, 教授, 博导。 DOI :10.16076/j .cnki .cjhd.2002.03.015
卢德唐等：多层地层中的井筒及地层温度解析解 383 2数学模型及其解 考虑多层的井筒瞬时热传导问题的温度分布如图1、图2所示，井筒中的流体通过对流传 递热量，然后通过热传导进入地层地层是由个不同的热力学及物理性质的多孔介质层组 成。整个系数由井筒区、热表皮区（包括套管、环空、水泥环等）及地层三部分组成。根据对问 题的研究，本文采用如下的近似及假设： 地层 ”0的0和 6) 10m0 120 160 图】井筒和地层温度分布曲线 图2给定时间下生产井温度剖面图（实例 (1)井筒中的流体为一维、垂向流动，且流量为常数： (2)同流动中的流体热对流相比，井筒中的流体垂直方向的热传导可忽略不计： (3)每个小层中的热力学参数、物理性质参数及初始温度梯度为常数： (4)与水平方向的热流量相比，地层中垂直方向的热传导可忽略不计： (5)用热表皮处理地层与井筒之间的热流量，同时引进热量储存常数： 根据以上假设，地层第j层的热传导方程为： (1) 在井筒中，流体的控制方程 对液体有 29C交+9g空=4=2要 (2a) 对气体有 (2h) http://www
2 数学模型及其解 考虑多层的井筒瞬时热传导问题的温度分布如图 1 、图 2 所示, 井筒中的流体通过对流传 递热量 ,然后通过热传导进入地层, 地层是由 n 个不同的热力学及物理性质的多孔介质层组 成。整个系数由井筒区、热表皮区 (包括套管、环空、水泥环等)及地层三部分组成 。根据对问 题的研究 ,本文采用如下的近似及假设 : 图 1 井筒和地层温度分布曲线 图 2 给定时间下生产井温度剖面图(实例) (1) 井筒中的流体为一维 、垂向流动 ,且流量为常数 ; (2) 同流动中的流体热对流相比,井筒中的流体垂直方向的热传导可忽略不计; (3) 每个小层中的热力学参数 、物理性质参数及初始温度梯度为常数 ; (4) 与水平方向的热流量相比 ,地层中垂直方向的热传导可忽略不计 ; (5) 用热表皮处理地层与井筒之间的热流量 ,同时引进热量储存常数 ; 根据以上假设, 地层第 j 层的热传导方程为: 1 r r rλj Tj r =ρjCj Tj t (1) 在井筒中 ,流体的控制方程 对液体有 πr 2 wρfCf Twj t +QρfCf Twj z =uj =2πrw Tj r |r =r w (2a) 对气体有 πr 2 wρfCf Twj t +QρfCf Twj z ± g Cf =uj =2 πrw Tj r |r =r w (2b) 卢德唐等:多层地层中的井筒及地层温度解析解 383
384 水动力学研究与进展 2002年第3期 由Ramey定义的综合热传导系数U,可表示成： 4=2要=-UT-T)- 3) 初始及边界条件可写成 Tg(z,t=0)=T(n2,t=0)=Tm(z) 4) T1(z=0,t)=T6 5) (2=0 1)=Tc 注入井 (6a) Teli(2=01)=To 生产井 (6) Tj(r,z=zt)=Tti(r,z=z 1) 7) limTj(r,z,t)=Ti (z) (8) 式 T一地层温度，℃)： T一井筒温度.C): p一地层中岩石及流体的综合密度，(kgm3): C一地层中岩石及流体的综合定压比热容，/kg.C): 地层中岩石及流体的综合热传导系数(W/(m.°C): 一井筒单位长度的热通量，(W/m): U 地层与井筒之间的综合热传导系数，(W/(m2.C): 一流体注入时的注入温度，（℃）： T0一对于注入井表示地面温度.对生产井为井底温度，(C): 2一流体的体积流量，(m/s): -重力加速度，(m/s2): w一井筒半径，(m). 下标 j=1,2.3…n(1<2<z 一第j层的物理量及参数： 一井筒中的物理量及参数： ∫一与流体有关的物理量及参数： 一初始状态的物理量及参数。 定义如下的无量纲量： (1)无量纲地层及井筒温度T,TD定义为 ?1994-2016 China Academic Journal Electronic Publishing House.All rights reserved.http://www
由 Ramey 定义的综合热传导系数 Uj ,可表示成: uj =2πrw Tj r |r =r w =-Uj(Twj -Tj)|r =r w (3) 初始及边界条件可写成: Twj(z , t =0)=Tj(r, z , t =0)=Tji(z) (4) T1i(z =0 , t)=T0 (5) Tw1i(z =0 , t)=TC 注入井 (6a) Tw1i(z =0 , t)=T0 生产井 (6b) Tj(r , z =zj , t)=Tj+1(r ,z =zj+1 , t) (7) limr※∞ Tj(r , z , t)=Tji(z) (8) 式中 T ———地层温度 ,(°C); Tw —井筒温度,(°C); ρ———地层中岩石及流体的综合密度 ,(kg/m 3); C ———地层中岩石及流体的综合定压比热容 ,(J/kg .°C)); λ———地层中岩石及流体的综合热传导系数,(W/(m.°C)); u ———井筒单位长度的热通量 ,(W/m); U ———地层与井筒之间的综合热传导系数 ,(W/(m 2.°C)); TC ———流体注入时的注入温度 ,(°C); T0 ———对于注入井表示地面温度, 对生产井为井底温度 ,(°C); Q ———流体的体积流量,(m 3 /s); g ———重力加速度 ,(m/s 2); rw ———井筒半径 ,(m)。 下标 j =1 , 2 , 3 …n (zj-1 <z <zj ——— 第 j 层的物理量及参数 ; w ———井筒中的物理量及参数; f ———与流体有关的物理量及参数 ; i ———初始状态的物理量及参数。 定义如下的无量纲量 : (1)无量纲地层及井筒温度 TDj , TWDj 定义为 384 水 动 力 学 研 究 与 进 展 2002 年第 3 期
卢德唐等：多层地层中的井筒及地层温度解析解 385 (2)无量纲时间及无量纲距离D,rD,D定义为 o=Cm=片n=, 2πg (3)热表皮、无量纲热储存常数及热力学参数比S,?,m定义为 9=铲房=防=C (C)λ (4)Laplace空间上的无量纲地层及井筒温度Tg,T四r定义为 To=Ty(ro.t)exp(-uD) Twy=Tge,o)pdo 根据上述定义的无量纲量，可以给出Laplace空间上温度所满足的方程及定解条件 9) (10) 阳--s器 (11) (12) Tg(rp→o∞u)=0 (13) Twi(zo=0u)=To 注入井 (14a） u TD1(zD=0.u)=0 生产井 (14h) Ta (液体) 式中D= 一对应于液体与气体的常数，(Cm): 2I994-2i人Z及击8Sa气体)Electron Publishing House,.ll rihts reserved.hup/Awnw
TDj = 2πλj(Tj -Tji) QρfCfDj ;TWDj = 2πλj(Twj -Tji) QρfCfDj (2)无量纲时间及无量纲距离 tD , rD , zD 定义为 tD = λ1 t (ρC)1r 2 w ;rD = r rw ;zD = 2 πλjz Q(ρC)f (3)热表皮 、无量纲热储存常数及热力学参数比 Sj , βj , mj 定义为 Sj = 2πλj Uj ;βj = (ρC)fλ1 2(ρC)1λj ;mj = (ρC)jλ1 (ρC)1λj (4)Laplace 空间上的无量纲地层及井筒温度 Dj , WDj 定义为 Dj =∫ ∞ 0 TDj(rD , zD , tD)exp(-utD)dtD WDj =∫ ∞ 0 TWDj(zD , tD)exp(-utD)dtD 根据上述定义的无量纲量 ,可以给出 Laplace 空间上温度所满足的方程及定解条件 1 rD rD rD Dj rD =mju Dj (9) uβj WDj -1 u + WDj zD = Dj rD |r D =1 (10) WDj = Dj -Sj Dj rD r D =1 (11) Dj λj WDj(zDj , u)= Dj-1 λj-1 WD(j-1)(zD(j-1), u) (12) Dj(rD ※ ∞, u)=0 (13) WD1(zD =0 , u)= TCD u 注入井 (14a) WD1(zD =0 , u)=0 生产井 (14b) 式中 Dj = TGj (液体) TGj ±g/ Ct (气体) ——— 对应于液体与气体的常数, (°C/m); 卢德唐等:多层地层中的井筒及地层温度解析解 385
386 水动力学研究与进展 2002年第3期 TG一第j层的地层静温梯度，(Cm): 求解上述方程可以得到Laplace空间上的无量纲地层及井筒温度Tor:T四分别为 Tman,）=G(uexH-E,(w)'z0+g)·u (15) Ko(rD umi )TwDi(zD.u) T(m,2,》）=Ko(Jm)+Sm冰1(m网） (16) 式中 umK(anj) 6m)=略+Ko(Jm)+Smk(两】 Ci(a)-1-E 注入井 C(u)=-E(u)·u 生产井 C(u)=9C1(u)espl[E(u)-E-(u月·r}+ 2 E1u)·u-写()·exu)·z0 7m=2(- QGD Ko(x)一零阶虚宗量Bessel函数： K1(x)一一阶虚宗量Besd函数. 使用围道积分可以得到方程(15)、(16)的L即ace解析反演解，即无量纲井筒温度、地层 温度分布实空间的解Tm(zD,o)及T(D,D,tD)(Laplace解析反演过程非常复杂，本文 仅给出反演后的结果)。 对第一层，无量纲井筒温度TwD(zD,D)可表示成 (tD≤D) TiDI (zD,to)= (17) f1+f2+3 (D>D) 式中 w2 = dw G()BI- 1w2 -Fw)]'+4 ?1994-2016 China Academic Journal Electro 。pm ng r se.All fights reserved.http://www
TGj ——— 第 j 层的地层静温梯度 ,(°C/m); 求解上述方程, 可以得到 Laplace 空间上的无量纲地层及井筒温度 Dj , WDj 分别为 WDj(zD , u)=Cj(u)exp[ -Ej(u)· zD] + 1 Ej(u)· u (15) Dj(rD , zD , u)= K0(rD umj) WDj(zD , u) K0( umj)+Sj umjK1( umj) (16) 式中 Ej(u)= uβj + umjK1( umj) K0( umj)+Sj umjK1( umj) C1(u)= TCD u - 1 E1(u)· u 注入井 C1(u)=- 1 E1(u)· u 生产井 Cj(u)= ΨjCj-1(u)exp{[ Ej(u)-Ej-1(u)] · zD(j-1)}+ Ψj Ej-1(u)· u - 1 Ej(u)· u exp[ Ej(u)· zD(j-1)] Ψj = Dj-1λj Djλj-1 ; TCD = 2 πλ1(TC -T1i) QρfCfD1 K0(x)———零阶虚宗量 Bessel 函数; K1(x)———一阶虚宗量 Bessel 函数 。 使用围道积分可以得到方程 (15)、(16)的 Laplace 解析反演解 ,即无量纲井筒温度、地层 温度分布实空间的解 TWDj(zD , tD)及 TDj(rD , zD , tD )(Laplace 解析反演过程非常复杂 ,本文 仅给出反演后的结果)。 对第一层, 无量纲井筒温度 TWDj(zD , tD)可表示成 TWD1(zD , tD)= f1 (tD ≤zD) f 1 +f 2 +f 3 (tD >zD) (17) 式中 f 1 = 4 π 2∫ ∞ 0 G1(w) w exp - w 2 m1 tD -1 G1(w) β1 - w 2 m1 -F1(w) 2 + 4 π 2 dw 386 水 动 力 学 研 究 与 进 展 2002 年第 3 期
卢德唐等：多层地层中的井筒及地层温度解析解 387 6=2me21-em-片o-o】， w Guf1-es-(2) -[2+0[oma-周-r时}ar 对于j=2,3,4n层，无量纲井筒温度可表示成 T晒(D,D)=(o)+JiT-1ar,t)-(c月B,(D,to-tMr (18) 式中 e) [9B-周-fm可+3r 0t≤(-0r1) B(zD,0)=2 「2(z-z0】· 春-@2_o=o+gdwo>o-g) exp G(w) G(w)=[Bw/o(w)+(1-Bw2S)Ji(w)]2+[BwYo(w)+(I-Bw2S)Y1(w)2 F(w)=Jo(w)BMJo(w)+(1-Bws;)J1(w)+ Yo(w)BwYo(w)+(1-BA2S)Y1(w) Jo(x,Yo(x)一第一类及第二类零阶Bsel函数： J1(x),Y1(x)一第一类及第二类一阶Bsel函数. 地层无量纲温度分布可表示成 店mab)=m-tir (19) -01China eademie Joumal Electronic Publishing House.ll rights reserved.hup:/
f 2 = 2exp(-β1zD) π ∫ ∞ 0 1 -exp - w 2 m1 (tD -zD) w exp zDF1(w) G1(w) sin 2zD πG1(w) dw f3 = 2exp(-β1zD) π ∫ ∞ 0 G1{1 -exp[ - w 2 m1 (tD -zD)] } w{[ G1(w)(β1 - w 2 m1 )-F1(w)] 2 + 4 π 2} exp zDF 1(w) G1(w) · 2 π cos 2zD πG1(w) +sin 2zD πG1(w) G1(w) β1 -w 2 m1 -F1(w) dw 对于 j =2 , 3 , 4 …n 层 ,无量纲井筒温度可表示成 TWDj(zD , tD)=Aj(tD)+∫ t D 0 [ TWDj-1(zDj-1 , τ)-Aj(τ)] Bj(zD , tD -τ)d τ (18) 式中 Aj(tD)= 4 π 2∫ ∞ 0 Gj(w) w exp - w 2 mj tD -1 Gj(w) βj -w 2 mj -Fj(w) 2 + 4 π 2 dw Bj(zD , tD)= 0 tD ≤(zD -zDj-1) 2 πmj exp[ -βj(zD -zDj-1)]∫ ∞ 0 usin 2(zD -zDj-1) πGj(w) · exp (zD -zDj-1)Fj(w) Gj(w) - w 2 mj (tD -zD +zDj-1)dw (tD >zD -zDj-1) Gj(w)=[ βjwJ 0(w)+(1 -βjw 2 Sj)J 1(w)] 2 +[ βjwY0(w)+(1 -βjw 2 Sj)Y1(w)] 2 Fj(w)=J 0(w)[ βjwJ 0(w)+(1 -βjw 2 Sj)J 1(w)] + Y0(w)[ βjwY0(w)+(1 -βjw 2 Sj)Y 1(w)] J 0(x), Y0(x)———第一类及第二类零阶 Bessel 函数; J 1(x), Y1(x)———第一类及第二类一阶 Bessel 函数。 地层无量纲温度分布可表示成 TDj(rD , zD , tD)=∫ t D 0 TWDj(zD , τ)gj(rD , tD -τ)dτ (19) 式中 卢德唐等:多层地层中的井筒及地层温度解析解 387
388 水动力学研究与进展 2002年第3期 {Yo(WID)[BWJo(w)+(1-Bwsj)J1(w)] Jo(wrp)[BYo(w)+(1-Bw2S)Y1(w)]}dw 3计算结果 方程(15)、(16)给出了无量纲井筒温度及地层温度分布的Laplace空间的解，可以使用 L即le数值反演给出它们的解：方程(17)、(18)、(19)给出了无量纲井筒温度及地层温度分布 实空间的解，但表达式非常复杂：为了进行比较本文也给出该问题的差分结果（计算结果由 WS井筒热模拟软件给出)。本文实例是注入情况的井筒及地层温度分布，并且，假定地层分 为两层，第一层为砂岩，第二层为粘土，计算所使用的参数如表1所示。图3是采用三种解得 到的不同时间下，井筒温度随井深的关系曲线，从比较结果来看三种计算结果在时间较大时趋 于一致。但当D≤zD时，Laplace数值反演结果不准确。所以，时间较小时的井筒温度解析解 是非常重要的。 表】计算所需的参数列表 名称 数值 单位 名称 数值 单位 地热梯度 0.03 m 油井半径 m 注入温度 c 第一层密度 2200 kg/m3 注入流量 100 第一层比热容 740 J/(kg.'c) 流体密度 986 kg/m 第一层热导率 28 W/(m.'C) 综合热传导系数 978 W/(m2C) 第二层密度 150 kg/m 地表温度 20 第二层比热容 800 /g.c) 产层深度 1000 m 第二层热导率 1.4 W/(m.'C) 为了研究垂直方向的热传导对地层温度分布的贡献，本文给出了不同时间下400处，垂 直方向的温度梯度与水平方向温度梯度之比对径向距离图（如图4所示），从图中可以看出： 时间在很大的范围内，垂直方向的温度梯度与水平方向温度梯度比值总是小于1%。这证明 了许多作者在研究多孔介质热传导时忽略垂直方向的热传导是正确的，但对于多层，如果相邻 两层的地层热力学参数差别较大在两层交界处附近，垂直方向的温度梯度可能较大，垂直方 向的热传导是否可以忽略要视具体情况而定。对于一般工程问题，忽略垂直方向的热传导是 非常合理的. ?1994-2016 China Academie Journal Electronic Publishing House.All rights reserved. http://www
gj(rD , tD)= 2 πmjβj∫ ∞ 0 w Gj(w) exp - w 2 mj tD {Y0(wrD)[ βjwJ 0(w)+(1 -βjw 2 Sj)J 1(w)] · J 0(wrD)[ βjwY0(w)+(1 -βjw 2 Sj)Y1(w)] }dw 3 计算结果 方程 (15)、(16)给出了无量纲井筒温度及地层温度分布的 Laplace 空间的解 , 可以使用 Laplace 数值反演给出它们的解 ;方程 (17)、(18)、(19)给出了无量纲井筒温度及地层温度分布 实空间的解 ,但表达式非常复杂 ;为了进行比较, 本文也给出该问题的差分结果 (计算结果由 WTS 井筒热模拟软件给出)。本文实例是注入情况的井筒及地层温度分布, 并且,假定地层分 为两层,第一层为砂岩,第二层为粘土, 计算所使用的参数如表 1 所示 。图 3 是采用三种解得 到的不同时间下 ,井筒温度随井深的关系曲线 ,从比较结果来看三种计算结果在时间较大时趋 于一致。但当 tD ≤zD 时 , Laplace 数值反演结果不准确 。所以, 时间较小时的井筒温度解析解 是非常重要的。 表 1 计算所需的参数列表 名称 数值 单位 名称 数值 单位 地热梯度 0 .03 °C/m 油井半径 0.1 m 注入温度 20 °C 第一层密度 2200 kg/m 3 注入流量 100 m 3 /d 第一层比热容 740 J/(kg .°C) 流体密度 986 kg/m 3 第一层热导率 2.8 W/(m .°C) 综合热传导系数 978 W/(m 2·°C) 第二层密度 1500 kg/m 3 地表温度 20 °C 第二层比热容 800 J/(kg .°C) 产层深度 1000 m 第二层热导率 1.4 W/(m .°C) 为了研究垂直方向的热传导对地层温度分布的贡献 ,本文给出了不同时间下 400m 处 ,垂 直方向的温度梯度与水平方向温度梯度之比对径向距离图 (如图 4 所示),从图中可以看出 : 时间在很大的范围内 ,垂直方向的温度梯度与水平方向温度梯度比值总是小于 1 %。这证明 了许多作者在研究多孔介质热传导时忽略垂直方向的热传导是正确的,但对于多层 ,如果相邻 两层的地层热力学参数差别较大,在两层交界处附近 , 垂直方向的温度梯度可能较大, 垂直方 向的热传导是否可以忽略要视具体情况而定 。对于一般工程问题, 忽略垂直方向的热传导是 非常合理的。 388 水 动 力 学 研 究 与 进 展 2002 年第 3 期
卢德唐等：多层地层中的井筒及地层温度解析解 389 40 15 过 05 t=120A 数值解 15o 0 M2 01.0203.04.050 La(r.) 图3井筒温度随井深的关系曲线比较 图4垂向与径向温度梯度之比 图5是注入5d时地层三维温度分布图，其中径向最大的计算距离为R=15m从图中可 以看出：用20C的水.以流量为m3/的注入量注入5d.温度影响半径较小，即在离井筒几米 以外，地层温度几乎没有受到扰动，地层温度分布仍然可以近似为地热静温曲线。但在井筒附 近温度变化较大，尤其在井筒中，温度随深度的变化曲线远偏离地热静温曲线。研究表明：注入 (或产出)流量是影响井筒中温度随深度的变化曲线偏离地热静温曲线最主要的参量。所以， 我们可以通过分析井筒中温度随深度的变化曲线或产层处的温度随时间变化的曲线来反求注 入（可产出）流量以及地层热力学参数 图5注入5对时地层三维温度分布图 ?1994-2016 China Academic Journal Electronic Publishing House.All rights reserved.http://www
图 3 井筒温度随井深的关系曲线比较 图 4 垂向与径向温度梯度之比 图 5 是注入 5d 时地层三维温度分布图,其中径向最大的计算距离为 R =15m, 从图中可 以看出 :用 20°C 的水,以流量为 m 3 /d 的注入量注入 5d ,温度影响半径较小, 即在离井筒几米 以外 ,地层温度几乎没有受到扰动 ,地层温度分布仍然可以近似为地热静温曲线 。但在井筒附 近温度变化较大 ,尤其在井筒中,温度随深度的变化曲线远偏离地热静温曲线, 研究表明 :注入 (或产出)流量是影响井筒中温度随深度的变化曲线偏离地热静温曲线最主要的参量 。所以 , 我们可以通过分析井筒中温度随深度的变化曲线或产层处的温度随时间变化的曲线来反求注 入 (可产出)流量以及地层热力学参数。 图 5 注入 5d 时地层三维温度分布图 卢德唐等:多层地层中的井筒及地层温度解析解 389
390 水动力学研究与进展 2002年第3期 4结论 通过本文的研究我们可以得到以下的结论 (I)可以用解析或半解析(Laplace变换及Laplace数值反变换)方法给出多孔介质中的 热传导问题的解，但时间较小(to≤D)时Laplace数值反变换得到的结果不准确； (2)研究多孔介质地层温度分布可以忽略垂直方向的热传导： (3)流量是影响多孔介质地层温度分布及井筒温度最主要的参数，因此，可以通过分析 井筒中温度随深度的变化曲线或产层处的温度随时间变化的曲线来反求注入（或产出）流量 以及地层热力学参数 参考文献 [I IESSEM L.B et al.A method of cakulating the distribution of temperature in flowing gasw.Trans..AMME 19572(10):169172. MOS J T.et a.How to temperature proiles in a water injection wel.Oil &GaJ..1959.57(11) 7年178 4.m.196144:427-435 A and PA MA 1981 5102 11 10013 Analytical solution of temperature in wellbore and formation in multi-laver LU De-tang,ZENG Yi-shan,GUO Yong-cun (University of Science and Technokgy of China.Hefei 230026.China) The mathematical model for wellbore and fom tod in muti-lacer The ed as lay s.The Laplace tained in both the Lplace and real spaces Key words: 1994-2016 China Academic Joural Electronic Publishing House.All rights reserved.http://www
4 结论 通过本文的研究我们可以得到以下的结论 (1) 可以用解析或半解析 (Laplace 变换及 Laplace 数值反变换)方法给出多孔介质中的 热传导问题的解 ,但时间较小 (tD ≤zD )时 Laplace 数值反变换得到的结果不准确 ; (2) 研究多孔介质地层温度分布可以忽略垂直方向的热传导 ; (3) 流量是影响多孔介质地层温度分布及井筒温度最主要的参数 , 因此 ,可以通过分析 井筒中温度随深度的变化曲线或产层处的温度随时间变化的曲线来反求注入 (或产出)流量 以及地层热力学参数 。 参 考 文 献 : [ 1] LESSEM L B, et al.A method of calculating the distribution of temperature in flowing gas wells[ J] .Trans., AIME, 1957, 2(10):169-172 . [ 2] MOSS J T, et al .How to calculate temperature profiles in a water injection well[ J] .Oil &Gas J ., 1959 , 57(11): 174-178. [ 3] RAMEY H J JR.Wellbore heat transmission[ J] .JPT , 1962, 14(4):427-435 . [ 4] BARELLI A and PALAMA A.A new method for evaluation formation equilibrium temperature in holes during drilling [ J] .Geothermics, 1981, 14(4):95-102 . [ 5] WOOLEY G R.Computing downhole temperatures in circulation , injection, and production wells[ J] .JPT, 1980 , 32 (12):1509-1522. [ 6] STEHFEST H .Numerical inversion of Laplace transforms[ J] .ACM, Communications, 1970, 13(1):47-49 . Analytical solution of temperature in wellbore and formation in multi-layer LU De-tang , ZENG Yi-shan , GUO Yong-cun (University of Science and Technology of China , Hefei 230026 ,China) Abstract: The mathematical model for wellbore and formations heat transmission were presented in multi-layer.The vertical non-homogeneous formations were approximated as layed formations with different physical properties.The Laplace transformation was used to the dimensionless equations of heat transmission.The closed-form analytical solutions were obtained in both the Laplace and real spaces. Key words : dimensionless temperature ;multi-porous ;heat transmission ;Laplace transformation ; analytical inversion of Laplace transforms 390 水 动 力 学 研 究 与 进 展 2002 年第 3 期