基金项目:国家自然科学基金(41702298,41572287); 陕西省自然科学基础研究计划(2017JQ4020); 陕西省教育厅专项科研计划项目(17JK0515)
通信作者:段 钊(1985-),男,陕西西安人,副教授,E-mail:landsliders@163.com
(1.西安科技大学 地质与环境学院,陕西 西安 710054; 2.国土资源部煤炭资源勘查与综合利用重点实验室,陕西 西安 710021)
(1.College of Geology and Environment,Xi'an University of Science and Technology,Xi'an 710054,China; 2.Key Experiment of Coal Resource Exploration and Comprehensive Utilization,Ministry of Land and Resources,Xi'an 710021,China)
U-shaped geothermal wells; formation temperature recovery; numerical simulation; best velocity
为研究U型地热井水平井段在采暖期间换热能力随注水流量的变化规律,并在非采暖期间不同流速且地层恢复至初始温度的前提下确定最优流速及最大取热量,以西安市深度为2 050 m的垂直埋深U型地热井水平井段取暖系统为研究对象,利用ANSYS FLUENT软件进行数值计算。根据恒定入口水温和流速的现场实际情况,建立了数值模拟的物理模型和数学模型,设置合理的边界条件,得到了U型地热井水平井段换热能力随流速的变化规律。研究表明在同一采暖期间,注水流量越大,出口水温越低,取热量越大。在流速为0.3 m/s即流量为33.85 t/h时地层在非取暖期间可以恢复至初始温度且水平井段最大取热量为1.81×1013J.
The purpose of this paper is to study the change rule of the heat transfer capacity of horizontal well section of u-shaped geothermal well with water injection flow during heating period and determine the optimal flow rate and maximum heat withdrawal under the premise of different flow rates during non-heating period and the formation recovery to the initial temperature.The u-shaped geothermal well heating system with a vertical buried depth of 2 050 m in Xi'an was taken as the research object,and the software ANSYS FLUENT was used for numerical calculation.According to the actual situation of constant inlet water temperature and flow velocity,the physical and mathematical models of numerical simulation were established,and reasonable boundary conditions were set up.The results show that during the same heating period,the greater the water injection flow is,the lower the outlet water temperature is,and the greater the heat extraction is.When the flow rate is 0.3 m/s,that is,the flow rate is 33.85 t/h,the formation can recover to the initial temperature during the non-heating period,and the maximum heat recovery of the horizontal well section is 1.81×1013J.
地热能作为一种可再生新型能源,已广泛应用到各大领域[1-2]。利用U型地热井对建筑物进行供暖是一种高效的地热能开发利用方式,在满足取暖要求的前提下,合理设计井身结构与水平段钻井长度,就可以达到最高的能量利用率[3-5]。因此掌握U型地热井水平井段换热能力随注水流速的变化规律对高效开发地热能起到关键作用。
目前,对地热井换热能力和井筒温度分布规律的研究有很多,郤保平等在试验研究、理论分析和数值模拟的基础上,提出高温岩体地热钻井施工中三大关键技术问题[6]; 于进洋利用反演法推导出原始地层温度分布以及井筒温度的解析[7]。何世明等根据传热学的基本原理,建立了井内液体与井筒之间热交换二维瞬态循环温度的数学模型对影响井下循环的参数进行了敏感性分析[8]; 刘杰等建立了开采井稳态传热模型,推导出单项流以及气液两相流情况下井筒温度的解析解[9]; 杨雪山等综合考虑了摩擦生热对井筒温度的影响,建立竖直、造斜和水平井段的井筒温度模型,并应用于实际工程,进行了验证[10]; 杨艳林等研究五点布井法在不同间距条件下观测生产井量和热提取率的变化特征[11]; 唐志伟等对生产井探测数据进行了核实并提出干热岩生产井液温度场影响因素的分析[12]; 冯绍航等在合理假设的基础上建立了采热井管道内外传热模型并对温度梯度、岩性变化进行了分析[13]; 袁益龙等对储层初始温度、裂隙间隔等影响因素进行分析得出渗透率对储层开采影响较小的观点[14]; 杨谋等分析了钻井过程中钻井液轴向导热和径向温度梯度对井筒的影响,得出钻井液径向温度梯度对井筒温度的影响规律[15-16]。秦立科等通过模拟实验揭示了岩石内部温度分布与演化的发展规律[17]。叶万军等对比雨水溶蚀前后试样粒度组成、矿物成分、离子含量及物理力学方面的差异[18]。刘之的等查明影响储层物的主控因素,提高油气勘探的成功率。结果表明井筒温度的分布随时间和空间具有一定的变化规律,井筒内的流体均向周围地层传热,其中环空部位为对流传热,而其他部位均为导热传热[19]。结果表明井筒温度的分布随时间和空间具有一定的变化规律,井筒内的流体均向周围地层传热,其中环空部位为对流传热,而其他部位均为导热传热。
ANSYS FLUENT目前是功能较为全面,对流体、热传递等适用最为广泛的流体力学软件。利用ANSYS FLUENT软件对关中地区U型地热井水平井段换热过程进行数值模拟,研究不同注水流量下水平井段换热能力的规律,并在非采暖期间不同流速和地层恢复至初始温度的前提下确定最优流速及最大取热量。
以西安市凹陷地区U型地热井水平井段为例,模拟研究水平井段在采暖期间换热能力随注水流量的变化规律并在非采暖期间不同流速且地层恢复至初始温度的前提下确定最优流速及最大取热量。西安凹陷是关中断陷盆地中的沉积中心之一,周边为4条深大断裂带所切围,凹陷内新生代地层厚为7 000 m,其中第四系地层厚达500~1 000 m.区内构造形迹主要表现为隐伏断裂构造,按其走向可分为EW向、NE向和NW向3组。
针对U型井的水平井段进行研究,根据研究区现场的实际工程情况,设置水平井段井长为120 m,井筒直径为200 mm,取水平井段周围50 m半径范围内的地层作为取热对象,水平井段埋深为2 050 m(图1)。
水平井段深度为2 050 m,为适当简化问题的研究,可忽略不同地层岩性对地温梯度的影响,根据关中地区平均的地温梯度设置水平井段周围地层温度为82 ℃,大位移水平井作为供水系统,水通过其垂直井段向下流动过程中吸收井筒外地层热量使自身温度不断升高,进入水平井段,但该段由于固井换热量并不高,直到水平井段的入口端,井筒内水和外界地层依然进行热传递使水的温度逐渐上升,设置模型水平段入口水温恒定为30 ℃.
对于U型地热井水平井段在取热过程中,井筒内水在流动的过程中与外界地层发生热传递并使井筒内水的温度逐渐上升,这种现象属于传热现象,水平井井筒内循环水与周围地层之间的传热模式包括水与管壁之间的对流传导和井筒周围地层之间的热传递。
井筒内水在流动的过程中,流体传热方程
(ρ·C)(∂T)/(∂t)+ρ·C·v·ΔT=·(KT)(1)
井筒内水在流动的过程中其温度场控制方程主要由动量方程、连续性方程和能量方程推到得出
k2T=ρc+Q(x,y,z)(2)
其初始条件为
T(x,y,z,0)=T0(3)
边界条件为
T|Γ1=Tw(x,y,z,t)Γ∈Γ1(4)
k(∂T)/(∂n)|Γ2=q(x,y,z,t)Γ∈Γ2(5)
式中 ρ为井筒内水的密度,kg·m-3; k为水的热传导系数,w·(m·K)-1; T为井筒内水的温度,℃; t为时间,s; c为水的比热容,J·(kg·℃)-1; Q为源热项,J; v为水流在井筒中的流速,m/s.
井筒内水的温度场方程、初始条件和边界条件构成了水在流动过程中与外界地层发生热传递的非稳态数学模型。
水平井段离地面深度为2 050 m,设置该深度地层温度为82 ℃,水平井段上覆地层和下伏地层大多是泥岩。水平井段中的水通过半径为200 mm,长度为120 m的水平井段与井筒外高温地层进行换热。利用Gambit建立模型(图1),且红色区域的边界EF定义为AXIS轴边界,边界CE定义为VELOCITY INLET、边界DF定义为PRESSURE OUTLET,边界 AC,AB,BD定义为WALL设置为地层边界,利用Quad网格型式及Map网格种类对模型进行网格划分且Interval size设置0.5如图2所示。因三维水平井段是轴对称旋转模型,所以通过对二维模型的研究表示三维水平井段,通过ANSYS FLUENT对模型数值模拟得出地层恢复到最佳状态下的最优流速及最大取热量。数值模拟中涉及到的地层、水的物性参数见表1.
通过对不同流速分别为0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.50,1.00,2.00 m/s情况下,分析时间与水平井段出口温度变化的关系,如图3所示。从图3可以看出,同一流速,随传热时间越长出口温度呈降低趋势且在同一时间段内,低流速的出口温度要高于高流速的出口温度。主要原因是:在同一流速下,随传热时间越长导致井筒中心处温度与管壁处的温度差较小,根据传热理论,流体温度与固体壁面温度的温度差越小,单位时间内的传热量越小,导致出口温度程逐渐降低的趋势。在水平井段内由于流速越大流量越大,在单位时间内携带传递的热量越大,导致同一时间段高流速的出口温度低于低流速的出口温度。
通过ANSYS FLUENT软件模拟在采暖期间不同初始流速条件下不同传热时间后水平井段出入口温度差,并得出以下数据见表2.
根据表2所得到的相关数据,利用ORIGIN拟合出在不同流速下,传热时间与水平井段出入口温度差的关系式如下所示。
v=0.10 m/s时,ΔT=51.78-0.002 4t+8.46×10-6t2(6)
式中 ΔT为出入口始末温度差,℃; v为水的初始流速,m/s; t为传热时间,s.
表2 不同传热时间后出入口的温度差
Table 2 Temperature difference table of entrance and exit after different heat transfer time℃
流速为0.15,0.2,0.25,0.3,0.35,0.4,0.5,1.0和2.0 m/s时的出入口温度差随传热时间的关系式均可由上述方法均拟合得出。
当水平井段内水的初始流速不同时,垂直埋深U型井工作一个采暖期(4个月即120 d)后,水平井段内出入口平均温度差及水携带的总热量关系式为
Δ(-overT)(∫1200f(t)dt)/(120)(7)
G=v·S·t·ρ(8)
Q=c·G·Δ(-overT)(9)
式中 Δ(-overT)为一个采暖期后水平井段出入口始末平均温度差,℃; t为时间,s; Q为一个采暖期后水平井段内水携带的总热量,J; G为水平井井筒内的循环水量,kg; c为水的比热容,J·(kg·℃)-1; v为水平井井筒入口端初始水流速,m/s; S为水平井井筒入口截面面积,m2; ρ为水的密度,kg·m-3.
利用公式(6)-(9)计算得出采暖期4个月后不同流速下出入口始末平均温度差和水平井段内水携带的总能量及总体取热量见表3,其对应曲线图如图4所示。
从图4可以看出,同一时间段,流速越快,水平井段的总取热量越大,主要原因是流速越快,相同截水断面下流量越大,根据湍流传热理论,水平井段里的水在单位时间内携带总体能量越大及取热量越大。
U型地热井在采暖期结束后,要进行地层温度的恢复阶段,以便在不破坏地层初始温度场的前提下,下个采暖期内U型地热井能够正常工作,能够获取最大取热量。通过ANSYS FLUENT软件模拟在不同流速下在非采暖期间地层恢复后的温度如图5所示。从图5能够看出,在非采暖期后地层恢复温度随采暖期注水流速的增长呈先水平后下降的趋势,流速为0.3 m/s时为拐点且地层恢复温度至81.845 ℃,最接近地层原始温度82 ℃,若高于0.3 m/s,非采暖期地层温度将无法恢复至初始温度。综上可以看出,为了不影响下一年度采暖期U型地热井的正常工作且单个取暖季可以获得最高取热量的前提下,应取最优流速为0.3 m/s,即33.85 t/h,对应最大取热量为1.81×1013J.
1)分析了流速对取热效果的影响,流速越慢,水平井段出口水温越高。且在同一流速下,随取热时间的增长水平井段出口水温呈下降趋势。
2)通过拟合公式计算采暖期(4个月)水平井段出入口平均温度差进而计算不同流速下的取热量,同一时间段内,U型井内水的初始流速越快,水平井段的取热量越大。
3)利用ANSYS FLUENT软件模拟计算了停止采热后非供暖期(240 d)内地层恢复效果。计算结果表明,经过240 d温度恢复后,当采暖期注水流速小于0.3 m/s时,地层均可恢复至初始温度,但为获得最大取热量应尽量选择高流速,因而最佳流速为0.3 m/s(即注水流量为33.85 t/h),对应水平井段最大取热量为1.81×1013J.