技术领域
[0001] 本
发明涉及一种井口
温度计算方法，特别涉及一种地热井井口温度计算方法。
背景技术
[0002] 沉积盆地型
地热田多为中低温地热田，井口温度一般较低，其井口温度对可利用
能源的数量有着直接的影响，因而其值的确定在中低温地热系统评价中占有举足轻重的地位，完善井口温度的计算方法对沉积盆地型地热田的
地热资源开发利用有重要的实际意义。井口温度一般受地温场的
地温梯度、
地热水在井筒中的流动速度以及地下
岩石的导热系数的影响。其中地温梯度的高低决定是否有可利用的
地热能，地热水流速较快、岩石热导率较低时，热量散失少，因而井口温度相应较高。
[0003] 1.传统井口温度计算方法
[0004] 依据管道工程中的舒霍夫公式：
[0005] TK＝T0+(TH-T0)·ea·L (1-1)
[0006] 式中：
[0007] 其中TK表示水的终点温度，℃；TH表示水的起点温度，℃；T0表示管道周围的介质温度，℃；L表示管道长度，m；a为常数；G为水的
质量流量，kg/s；C为水的
比热，J/kg·℃；D为管道的外直径，m；K为管道的总传导系数，一般来讲。
[0008]
[0009] 其中α1为
流体与管道内壁的内部放热系数；α2为管道外壁至周围介质的外部放热系数；δ为管道的壁厚；λ为管道壁的热导率。对于不保温的输
热管道的总
传热系数K近似等于2α。由于沉积盆地中的地热资源往往存在于油气田中，进行地热研究的井多数为油井，且多为
水泥固井。水泥的热导率为1.43，与岩石接近，所以可以把它看做岩石的一部分，这样可以认为输水管道是不保温的，其总传热系数K近似等于2α，即K值近似等于管道外壁至周围介质的外部放热系数。
[0010] 具有一定温度的水在井筒内向上流动过程中，井筒周围介质温度是(地温场温度)逐渐降低，即公式1-1中的T0是变化的，因而必须采用分段计算的办法，逐步向井口
迭代。然而，这样建立的传热方程是在稳态条件下的虽然解决了从井底到井口热量损耗量的计算，但这种计算方法在反映井筒周围的瞬间温度场及井口温度随时间的变化过程还存在明显不足。
发明内容
[0011] 本发明的目的在于提供一种地热井井口温度计算方法，考虑到需要在多层热储中对地热资源的混采，通过模型计算井口温度，可以很明显的反映井筒周围的瞬间温度场及井口温度随时间的变化过程，以解决上述背景技术中提出的问题。
[0012] 为实现上述目的，本发明提供如下技术方案：
[0013] 一种地热井井口温度计算方法，通过模型计算井口温度，计算步骤包括：
[0014] 步骤1：根据已知
地层温度计算地温梯度；
[0015] 步骤2：计算总传热系数U及A值；
[0016] 步骤3：由于水的比热随深度而变，故由浅至深逐步迭代计算至目标深度的流体温度；
[0017] 步骤4：对于单产层而言，可直接迭代计算到井口，而对于多产层而言，首先把各单产层水的产量相加得到新的总产水量，然后根据水的比热、温度和产量参数值计算得到井口温度。
[0018] 进一步地，计算井口温度为地下目的层位的地下温度，减去热水开采过程中在井筒中的
热损失。
[0019] 进一步地，根据热动态平衡原理，流体的热损失等于向水泥环传递的热量：
[0020] dq＝-W·CdT1 (1-3)
[0021] 式中：dq＝2πr3U(T1-T2) (-dZ)。
[0022] 进一步地，从水泥环向地层传导的热量可表示为：
[0023]
[0024] 其中Te表示未受干扰的地温场温度：Te＝a·(z-z*)+b。
[0025] 与
现有技术相比，本发明的有益效果是：本地热井井口温度计算方法，考虑到需要在多层热储中对地热资源的混采，通过模型计算井口温度，根据已知地层温度计算地温梯度；计算总传热系数U及A值；由于水的比热随深度而变，故由浅至深逐步迭代计算至目标深度的流体温度；对于单产层而言，可直接迭代计算到井口，而对于多产层而言，首先把各单产层水的产量相加得到新的总产水量，然后根据水的比热、温度和产量参数值计算得到井口温度，可以很明显的反映井筒周围的瞬间温度场及井口温度随时间的变化过程。
附图说明
[0026] 图1为本发明的地热井井身结构示意图；
[0027] 图2为本发明的无限长径向系统的瞬时热传导图；
[0028] 图3为本发明的井口温度计算流程。
具体实施方式
[0029] 下面将结合本发明
实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0030] 一种地热井井口温度计算方法，通过模型计算井口温度，计算井口温度，其实质问题是将地下目的层位的地下温度，减去热水开采过程中在井筒中的热损失。因此，计算的重点即为井筒的热损失(根据如图1所示的井筒结构，其中Z表示深度，①为
套管，②为地面，尖头为流体流动方向)，计算步骤包括：
[0031] 步骤1：根据已知地层温度计算地温梯度；
[0032] 步骤2：计算总传热系数U及A值；
[0033] 步骤3：由于水的比热随深度而变，故由浅至深逐步迭代计算至目标深度的流体温度；
[0034] 步骤4：对于单产层而言，可直接迭代计算到井口，而对于多产层而言，首先把各单产层水的产量相加得到新的总产水量，然后根据水的比热、温度和产量参数值计算得到井口温度，可以很明显的反映井筒周围的瞬间温度场及井口温度随时间的变化过程，计算流程见图3。
[0035] 根据热动态平衡原理，流体的热损失等于向水泥环传递的热量：
[0036] dq＝-W·CdT1 (1-3)
[0037] 式中：dq＝2πr3U(T1-T2) (-dZ)。
[0038] 从水泥环向地层传导的热量可表示为：
[0039]
[0040] 其中Te表示未受干扰的地温场温度：Te＝a·(z-z*)+b。
[0041] 由公式(1-3、1-4)可建立如下方程：
[0042]
[0043] 式中：
[0044] 对方程(1-5)积分可得：
[0045]
[0046] 其中，T0为在深度z0处水的温度，℃；T为深度z处水的温度，℃；z表示深度，m；z*表示常温层深度，m；a为地温梯度，℃/m；b为常温层温度，℃；W为水的产量，kg/s；c为水的比热，J/kg·℃，λ为岩层热导率，W/(m·℃)；U为总传热系数，W/(m·℃)；f(t)为时间函数，见图2。
[0047] 当t大于7天时，可用公式(1-7)计算：
[0048]
[0049] 其中α表示岩石的扩散系数，m2/s。
[0050] 当日产量较大时，总传导系数： 其中λ0表示水泥环的热导率，W/(m·℃)；r2、r3的选取见图1。当日产量较低时，需结合井筒中流体的
密度、线速度及
粘度进行计算，而由于产量降低，需综合分析该井是否仍有地热资源的开发利用价值。
[0051] 根据图3所示，图中T为温度；Z为深度；λ为岩石热导率；Q为产水量；a、b、c分别为各产水层编号。对于a、b、c三层，已知Za、Ta、Qa，计算到b层时的温度Ta→b，得到Ta→b层及Qa，结合已知的Zb、Tb、Qb，计算a、b层混合温度Tab、水量Qa+Qb，然后计算a、b两层混合后到c层的温度Tab→c，由Tab→c层、Qa+Qb结合已知的Zc、Tc、Qc，迭代至井口的温度T0，总水量Q。
[0052] 综上所述：本地热井井口温度计算方法，考虑到需要在多层热储中对地热资源的混采，通过模型计算井口温度，根据已知地层温度计算地温梯度；计算总传热系数U及A值；由于水的比热随深度而变，故由浅至深逐步迭代计算至目标深度的流体温度；对于单产层而言，可直接迭代计算到井口，而对于多产层而言，首先把各单产层水的产量相加得到新的总产水量，然后根据水的比热、温度和产量参数值计算得到井口温度，可以很明显的反映井筒周围的瞬间温度场及井口温度随时间的变化过程。
[0053] 以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉
本技术领域的技术人员在本发明披露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。
