豫东通许凸起回灌条件下新近系地热资源优化开采方案研究.pdf

地热井的优化开采对地热资源的集约高效利用意义重大。本文以豫东通许凸起中部地热地质单元为研究对象,采用F E F L OW软件建立了新近系热储三维水热耦合模型,利用回灌试验成果预测了全区在“仅开采模式”、“一采二灌(采灌比3 8.4 1%)”、“一采三灌(采灌比1 0 0%)”下自2 0 2 0-2 0 3 0年开采年限内的地热流体压力场和温度场变化情况,并计算了新近系热储不同开采模式下可采热能。预测结果表明:回灌率1 0 0%时,研究区地热流体压力水头回升更明显,回灌引起研究区新近系热储温度下降更微弱,获得可采地热能最大为0.8 5 21 01 6J。在此基础上提出了最优化的采灌布井模式,为通许凸起深部地热资源可持续开发提供理论依据。关键词:回灌;模型;新近系热储;温度场;压力场;可采热能中图分类号:P 3 1 4.1 文献标识码:文章编号:1 0 0 1-5 6 6 3(2 0 2 3)0 3-0 6 4 6-1 10 引言目前,国际国内普遍认为回灌式开采是一种新的地热开发利用方法。只有回灌,才能实现热储压力采灌平衡,才能使地热资源成为真正的可再生能源1。河南省地下水管理办法规定建设需要取水的地热能开发利用项目应实行同一含水层等量取水和回灌2。豫东通许凸起中部地热地质单元内地热水回灌式开发利用的主要途径为地热供暖3。对于地热供暖项目,其工作原理主要是提取地热水中的可采热能,开采回灌过程中水量基本不产生消耗。因此研究区内可回灌条件下的地热资源优化开采配置对于保证深部地热资源可持续发展意义重大4,5。2 0世纪7 0年代末,德国WA S Y水资源规划和系统研究所开发了基于有限单元法的F E F L OW软件,它是迄今为止功能最为齐全的地下水模拟软件包之一。在F E F L OW系统中,用户可以很方便迅速地产生空间有限单元网格,设置模型的参数和定义边界条件,运行数值模拟以及实时图形显示结果与成图6。F E F L OW发展过程中经过了大量的测试和检验,它成功地解决了一系列与地下水有关的实质性问题,如判断污染物迁移途径、追溯污染物的来源,开展地热的模拟等6。翟美静7利用F E F L OW软件开展了西安市三桥地区地热采灌系统数值模拟研究,根据研究区地热井实测水位和温度资料对模型进行了拟合验证,确定了模型的可靠性。程万庆等8利用F E F l o w求解地热采灌对井系统水热耦合模型,进行温度场模拟研究。本文以研究区施工的D Z K 1地热钻孔新近系热储层回灌试验数据为基础,开展区内新近系热储温度场和压力场模拟预测,为研究区合理开发利用深部地热资源提供依据。1 研究区地热地质概况研究区位于通许凸起构造单元中部(图1),区内地热收稿日期:2 0 2 3-0 1-3 1;修回日期:2 0 2 3-0 3-0 3基金项目:河南省自然资源厅2 0 1 9年度省财政地质勘查项目“河南省通许县岩溶热储地热资源调查”(4 9号)资助。第一作者:任静(1 9 8 3),女,硕士,高级工程师,主要从事水工环勘查,地热资源勘查工作。E-m a i l:2 5 0 5 3 4 0 3 7q q.c o m引文格式:任静,邓晓颖,郭彦.豫东通许凸起回灌条件下新近系地热资源优化开采方案研究J.矿产与地质,2 0 2 3,3 7(3):6 4 6-6 5 6.资源的形成是地球深部的热源以传导的形式对新生界及古生界地层影响的结果,地热资源类型属于沉降盆地传导型,平面上呈层状展布,根据热储成因模式、储层时代、分布及储水介质特征,划分为新近系明化镇组(Nm)、馆陶组(Ng)、寒武-奥陶系(O)三个热储层。新近系孔隙热储层是本区地热流体的主要开采层,具有分布面积大,产状倾角较缓,沉积厚度较大的特点。其储 水 介 质 岩 性 主 要 为 细 砂(岩)、粉 细 砂(岩),局部为中粒砂(岩)。隔水层岩性为厚层黏土及粉质黏土、泥岩。地热流体储层有砂(砂岩)层2 5-2 9层,单层厚度一般为3.55 5 m,累积总厚平均为5 0 08 0 0 m。根据地层岩性特征新近系热储分为新近系明化镇组(Nm)和馆陶组(Ng)。目前此层开采主要用于洗浴和城区小区地热供暖。图1 研究区地热地质简图F i g.1 G e o t h e r m a l a n d g e o l o g i c a l m a p o f t h e s t u d y a r e a1侏罗-白垩系未分 2三叠系潭庄组 3三叠系刘家沟组及二马营组 4二叠系上石盒子组 5二叠系上统未分 6二叠系山西组及下石盒子组 7石炭系 8奥陶系 9寒武系 1 0古元古界五指岭组 1 1断层及隐伏断层 1 2地层分界线 1 3岩相界线 1 4前新生界基岩深度线 1 5研究区边界 1 6三级构造单元分区线 1 7钻孔位置及底板埋深2 新近系热储层地热地质特征和回灌能力2.1 新近系热储地热地质特征根据D Z K 1井揭露地层情况,新近系明化镇组由浅黄色、黄褐色粉砂岩、中细砂岩、细砂岩与棕红、紫红色、浅棕红色泥岩互层组成;新近系馆陶组岩性为紫红色泥岩与浅黄色、黄褐色中细砂岩、细砂岩互层,与下伏地层呈不整合接触。两组岩层无显著隔水边界存在。全区新近系热储层有砂(砂岩)层2 52 9层,单层厚度一般为3.5 5 5 m,累积总厚平均为5 0 08 0 0 m;顶板埋深2 4 03 8 0 m,底板埋深9 0 01 2 0 0 m;孔隙度为2 7.1 0%。新近系热储层顶、底部均具有厚层黏746 第3 7卷 第3期 任静等:豫东通许凸起回灌条件下新近系地热资源优化开采方案研究土(泥岩)可作为相对隔水层,故垂向视为无越流作用;地热水径流微弱,开采主要消耗含水砂层弹性储量。热源供给主要为大地热流传导方式。2.2 新近系热储回灌能力确定本次利用区内最新施工的D Z K 1钻孔开展新近系热储回灌试验。D Z K 1地热勘探孔成井深度2 0 0 0 m,新近系取水段8 7 0 1 0 8 5 m(新近系馆陶组),成井工艺为筛管填砾成井。该井抽水试验结束后利用该井同层位进行无压重力自流回灌方式,采用浅层地下水作为回灌水源。本次试验采用流量表、水位测量仪计量流量、水位的变化,回灌试验自2 0 2 0年3月2 0日1 0:3 0开始至3月2 5日1 3:2 0结束,累计回灌时间持续7 3 1 0 m i n,持续回灌时间1 2 1.8 3 3 h,回灌流体温度1 7,稳定回灌量3 5.1 2 m3/h,回灌前静水位埋深6 0.2 6 m。回灌后稳定动水位9.2 3 m,稳定时间1 0 h。本次回灌试验基本数据见表1,历时曲线见图2。表1 D Z K 1井回灌试验基本数据T a b l e 1 L i s t o f b a s i c d a t a o f D Z K 1 w e l l r e i n j e c t i o n t e s t观测时间稳定灌量/(m3h-1)回灌水温/稳定动水位埋深/m稳定时间/h2 0 2 0年3月2 0日1 0:3 03月2 5日1 3:2 03 5.1 21 79.2 31 0图2 D Z K 1井回灌试验历时曲线图F i g.2 D Z K 1 W e l l r e i n j e c t i o n t e s t c u r v e3 新近系热储层模拟利用F E F L OW软件开展“直接开采”、“一采2灌”、“一采3灌”不同开采方案的新近系热储水热耦合模型预测。模型模拟流程见图3。3.1 边界条件概化(1)模拟范围根据研究区地质、地热地质条件,以研究区范围作为模拟范围,面积为1 7 5 3.9 8 k m2,见图4。(2)边界概化根据区内新近系热储层2 0 1 9年5月和2 0 2 0年7月的压(力)头流场(自南向北径流),将模拟区北部边界概化为零通量边界,西侧、南侧和东侧边界概化为补给边界,如图3所示。垂向上,新近系上部由于第四系岩层的覆盖,不考虑大气降水入渗补给量,加之图3 新近系热储水热耦合模型模拟流程tF i g.3 S i m u l a t i o n p r o g r a m o f c o u p l e d m o d e l o f w a t e r a n dh e a t i n N e o g e n e t h e r m a l r e s e r v o i r846矿 产 与 地 质2 0 2 3年 图4 边界概化示意图F i g.4 S c h e m a t i c d i a g r a m o f b o u n d a r y g e n e r a l i z a t i o n1零通量边界 2补给边界 3新近系热储压力等值线 4新近系地热水流向规模性黏土层的作用,概化为隔水边界面;下部与下伏二叠系热储层基本无水力联系,根据地层岩性分布同样概化为隔水边界面。(3)含水层概化根据D Z K 1井揭露地层情况,新近系明化镇组由浅黄色、黄褐色粉砂岩、中细砂岩、细砂岩与棕红、紫红色、浅棕红色泥岩互层组成。新近系馆陶组岩性为紫红色泥岩与浅黄色、黄褐色中细砂岩、细砂岩互层,与下伏地层呈不整合接触。两组岩层无显著隔水边界存在,因此将新近系地层统一概化为一个含水层,见图5。图5 热储层概化后示意图F i g.5 S c h e m a t i c d i a g r a m o f t h e t h e r m a l r e s e r v o i r a f t e r g e n e r a l i z a t i o n946 第3 7卷 第3期 任静等:豫东通许凸起回灌条件下新近系地热资源优化开采方案研究 (4)数学模型地下水流动模型在上述水文地质概念模型的基础上,根据达西渗流定律和渗流连续性方程,将研究区地下水流用以下方程和定解条件描述:xKx xHx +yKy yHy +zKz zHz +=s Ht,x,y,z ,t0;Kx xHx 2+Ky yHy 2-Kz zHz +W=dHt,x,y,z 0,t0;Knn|1=q x,y,z,t ,x,y,z 0,t0;H x,y,z,t|t=0=H0 x,y,t ,x,y,z 。式中:为模拟范围;H为热储层水头;Kx x、Ky y、Kz z为x,y,z方向上的渗透系数(m/d);Kn为边界法线方 向上的渗透 系数(m/d);s为 单 位 储 水系数(l/m);d为重力给水度;为源汇项(l/d);0为上边界;为第二类边界;n为研究区边界外法线方向;q(x,y,z,t)为第二类边界单宽流量(m/d)。热量运移方程tt-1Ce q(i jT)+CLCe qqfiT+QCQQCe q(T-Tq)=QTCe q;T(x,y,z,0)=T0,(x,y,z);T(x,y,z,t)|1=T1,(x,y,z)1;x xTx+x yTy+x zTz c o s(nx)+y xTx+y yTy+y zTz c o s(ny)+z xTx+z yTy+z zTz c o s(nz)|2=qr,(x,y,z)2。式中:Ce q为等效体积热容;Ce q=(C)/。其中表示多孔介质中水所占的体积分数,即孔隙度,C表示多孔介质中(包括水)总的体积热容,单位:J/(m3K);CL为wCw,表示水的体积热容,单位J/(m3K);w为水的密度(k g/m3);Cw为比热容,单位J/(k gK);为等效热传导系数,W/(mKd);QT为一般热源(W/m2);qT为边界2处流入的热通量(W/m2);其他符号意义同上。3.2 模型识别与验证(1)模型赋值在水文地质概念模型中,采用给定水头边界和流量边界来进行模型区含水层侧向补给与排泄的赋值。其中给定水头边界由软件根据周边压(力)头流场分布来自动计算边界地下水的补排量;流量边界所赋数值由人工计算而来,具体见表2。表2 水文地质概念模型边界参数值T a b l e 2 B o u n d a r y p a r a m e t e r v a l u e s o f h y d r o g e o l o g i c a lc o n c e p t u a l m o d e l边界 边界渗透系数 K(md-1)南边界-2.6 4 0 21 0-4北边界1.4 1 91 0-4 注:正值表示排泄流量,负值表示补给流量。(Kx、Ky、Kz)、弹性释水率(Ss)、孔隙度(n)。综合抽水试验资料及相关经验值,对各参数赋值初始值见表3。表3 水文地质参数初值T a b l e 3 I n i t i a l v a l u e o f h y d r o g e o l o g i c a l p a r a m e t e r s水文地质参数参数值参数单位Kx0.4m/dKy0.4m/dKz0.0 4m/dSs11 0-51/mn0.2无(2)水文地质参数选取建模过程中选用的水文地质参数有:渗透系数(3)模型识别与验证根据本次模拟研究的需要,以2 0 1 9年5月2 0 2 0年7月为模型综合识别验证期,以2 0 2 0年7月流场作为模拟 流 场 拟 合 验 证 的 目 标。经 过“对照-调参”的 参 数 循 环 优 化 过 程,模 拟 结 果 的 平 均相对误差小于1 0%,满足 地下水资源管理模型工作要求 中对水位拟合的要求,所建地下水模型可用于后续模拟研究。最终获得优化后的水文地质参数见表4。056矿 产 与 地 质2 0 2 3年 表4 水文地质参数优化值T a b l e 4 L i s t o f o p t i m u m v a l u e o f h y d r o g e o l o g i c a l p a r a m e t e r s水文地质参数参数值参数单位Kx0.4 3m/dKy0.4 3m/dKz0.0 4 3m/dSs11 0-61/mn0.2无3.3 热传导数值模型建立在已建立地下水流模型的基础上,耦合热传导参数,建立全区新近系热储热传导模型。(1)热量运移边界顶板边界研究区恒温带深度取2 3 m,恒温带温度确定为1 5.9(研究区年平均气温1 4.2)。而新近系顶层埋深在2 4 03 8 0 m,已处于增温带,其热量来源为下部热源的对流与弥散传导,因此以给定的温度场作为初始状态,由模型根据条件计算其温度变化。底板边界通过数据分析确定通许凸起中部地区的大地热流值6 2.2 mW/m29。(2)模拟校正检验研究区新近系热储的水温动态在2 0 1 92 0 2 0年处于稳定状态。在此基础上,当进行温度场模拟时,只要计算温度跟实测温度变化不大,则表示模型符合实际情况。在此对比模拟期前后温度场变化,如图6所示,可以看出两者差别极弱,说明本次所建地热模型稳定可靠,可用于后续模拟研究。图6 研究区新近系热储温度场模拟F i g.6 T e m p e r a t u r e f i e l d s i m u l a t i o n o f N e o g e n e t h e r m a l r e s e r v o i r(a)初始温度场 (b)期末温度场4 地热资源开采方案预测在已建立地热模型的基础上,结合研究区开采现状及规划,并在考虑供暖期、停暖期的情况下,进行水热耦合模拟,研究地热资源不同开发利用情景对整个区域压力场和温度场的影响1 0。方案一:直接开采方案。区内调查了4 7眼地热井开采新近系热储的资料,本方案设定为:按照现状分布的4 7眼地热井使其保持常年开采,其它条件不变1 1。以2 0 2 0年7月地热流体压力场为预测模型初始流场,以现状资料温度场为初始温度场,将模型模拟预测期定为2 0 2 0年7月2 0 3 0年7月,共1 0年,3 6 5 0天。方案二:一采二灌方案。基于调查的4 7眼开采井分布位置,采用“一采二灌”的地热井采灌组合形式,单井开采量设为1 6 8 0 m3/d(7 0 m3/h),单井回灌流量设为8 4 0 m3/d(3 5 m3/h),回灌率:3 8.4 1%。在研究区现有的各地热开采井附近设置两个回灌虚拟井点,回灌点至开采点的距离在1 0 0 01 5 0 0 m之间,模拟预测的开采和回灌时间段定为当年1 1月至次年3月,停止开采及回灌时间段为4月至1 0月,根据相关规范及资料选择2 5为回灌水温度。方案三:一采三灌方案。基于调查的4 7眼开采井分布位置,采用“一采三灌”的地热井采灌组合形式,单井开采量设为2 4 0 0 m3/d(1 0 0 m3/h),单井回灌流量设为8 4 0 m3/d(3 5 m3/h)。回灌率:1 0 0%。在研究区现有的各地热开采井附近设置三个回灌虚拟井点,回灌点至开采点的距离在1 0 0 01 5 0 0 m之156 第3 7卷 第3期 任静等:豫东通许凸起回灌条件下新近系地热资源优化开采方案研究间,模拟预测的开采和回灌时间段定为当年1 1月至次年3月,停止开采及回灌时间段为4月至1 0月,选择2 5为回灌水温度。4.1 方案一现状直接开采模式下研究区地热流体压力场模拟预测结果如图7所示,可以看出在地热井的持续开采下模拟区地热流体水头整体发生不同程度下降,并在北侧开采井集中区形成了小规模的降落漏斗,水头下降至约6 m。分析流场变化原因,主要为模拟区地热水流系统补给项较为单一,为侧向径流补给,人工地热井群共4 7眼,单井抽水量在2 4 08 4 0 m3/d,而热储层渗透系数较小,使得区域压力头持续下降,并在通许县城井群集中区发育水头降落漏斗。温度场预测结果如图8所示,可以看出在地热井图7 新近系热储层模拟压力场(2 0 3 0年)F i g.7 S i m u l a t e d p r e s s u r e f i e l d o f N e o g e n e t h e r m a lr e s e r v o i r i n 2 0 3 0图8 研究区2 0 3 0年温度模拟结果F i g.8 S i m u l a t e d t e m p e r a t u r e f i e l d o f t h e s t u d y a r e a 2 0 3 0(a)顶层模拟温度场(b)底层模拟温度场的持续开采下模拟期末温度场变化甚微。观测孔顶层点温度维持在2 62 6.0 3 5 5,底层点温度维持在4 5.2 6 8 44 5.2 4 9 3。在热能方面,按照现 状 开 采 条 件,考 虑 开 采 水 温 及 底 线 温 度(2 5),计算得到2 0 2 02 0 3 0年区内开采可获得热能为0.2 8 11 01 6J。综上所述,在该方案下维持现状地热井直接开采模式,会使得研究区区域地热流体压力持续下降,并在井群集中区发育水头降落漏斗。而温度场方面由于通许凸起为高温异常区,在现状地热开采条件下,可维持温度场的稳定。4.2 方案二在现状直接开采条件的基础上对区域性地热流体回灌做进一步的模拟研究,以探究在供暖期进行区域性回灌条件下,研究区新近系热储流场和温度场的变化情况。具体设置如下:基于现有开采井,采用“一采两灌”的地热井采灌组合形式,在研究区现有的各地热开采井附近设置两个回灌虚拟井点,根据目前模型三角网格节点剖分情况,选取回灌点至开采点的距离在1 0 0 01 5 0 0 m之间,以使回灌点与三角网格节点重合。具体点位分布如图9所示:256矿 产 与 地 质2 0 2 3年 图9 地热井采及回灌井分布简图F i g.9 D i s t r i b u t i o n d i a g r a m o f g e o t h e r m a l w e l lp r o d u c t i o n a n d r e i n j e c t i o n w e l l s模拟区压力场预测结果如图1 0所示。可以看出,该方案下研究区内地热开采井附近地热流体水头降落漏斗的规模有所增大,表现为平面上扩展有限,而垂向上水头值下降明显。为定量化研究该方案回灌对研究区流场的影响,并验证上述分析,使用“数值模拟+G I S地理空间分析”的方法,对模拟结果进行处理,得到方案二回灌条件下研究区新近系热储地热流体回灌成效分布图1 0。结合图1 1可以看出,回灌措施下以回灌井为中心产生规模不等的“水丘”,其水头回升值大多在0.1 50.7 8 m,在研究区北部的群井回灌中甚至产生了0.7 83.4 5 m的回升值,中部部分开采井水头则呈现出进一步下降的现象,下降值大多在0.1 7 1.2 9 m,局部达5.5 m,发现其“两灌一采”组合之间距离较远,且水头成效具有明显的分带性。从图1 2(a)可以看出,在方案二条件下进行区域性供暖期回灌时,2 0 3 0年研究区底层温度场整体保持稳定,但是在各回灌点附近小范围的温度下降更加明显。为定量该方案回灌对研究区温度场的影响,利用G I S地理空间分析方法对地热模拟结果进一步处理,得到由回灌引起的底层温度下降带分布图1 2(b)。可以看出,由回灌引起的温度下降大多在“0.0 5 90.1 1 4”区间,局部群井回灌则会引起热储温度进一步下降,最大可达1.0 8。受此影响,区域上温度下降主要在00.0 5 9。“一采两灌”开采时温度下降值分布特征明显,呈分别以两回灌点为中心向开采井递减,且开采井附近地温受此影响的下降值在“00.0 2 4”区间,即局部出现热失衡现象,该现象与采灌井间距布局有关。在地热能方面,2 0 2 02 0 3 0年考虑开采水温及底线温度(2 5),计算得到1 0年内一采二灌模式下可获得的热能为0.5 9 61 01 6J,比直接消耗式开采可获得的热能高0.3 1 51 01 6J,说明回灌有助于增加深部大地热流传导补给被消耗的地热资源。图1 0 新近系热储层模拟压力场(2 0 3 0年)F i g.1 0 S i m u l a t e d p r e s s u r e f i e l d o fN e o g e n e t h e r m a l r e s e r v o i r i n 2 0 3 0图1 1 地热流体开采、回灌成效分布图F i g.1 1 D i s t r i b u t i o n m a p o f g e o t h e r m a l f l u i de x p l o i t a t i o n a n d r e i n j e c t i o n e f f e c t356 第3 7卷 第3期 任静等:豫东通许凸起回灌条件下新近系地热资源优化开采方案研究图1 2 研究区2 0 3 0年温度模拟结果F i g.1 2 t e m p e r a t u r e s i m u l a t i o n r e s u l t o f t h e s t u d y a r e a i n 2 0 3 0(a)顶层温度模拟 (b)由回灌引起的底层温度下降分布4.3 方案三基于现有开采井,采用“一采三灌”的地热井采灌组合形式,在研究区现有各地热开采井附近设置三个回灌虚拟井点,根据目前模型三角网格节点剖分情况,选取回灌点至开采点的距离在1 0 0 01 5 0 0 m之间,以使回灌点与三角网格节点重合。具体点位分布如图1 3所示:图1 3 地热采井及回灌井分布简图F i g.1 3 D i s t r i b u t i o n d i a g r a m o f g e o t h e r m a l w e l lp r o d u c t i o n a n d r e i n j e c t i o n w e l l s模拟区流场预测结果如图1 4所示。可以看出,该方案下模拟区开采井附近地下水头降落漏斗的规模相较方案一进一步增大,表现为水头值下降更加明显。为定量化研究该方案回灌对研究区流场的影响,并验证上述分析,使用“数值模拟+G I S地理空间分析”的方法,对模拟结果进行处理,得到方案三回灌条件下研究区地下水回灌成效分布图1 5。结合图1 5可以看出,回灌措施下以回灌井为中心产生的“水丘”更高,其成效水头值大多在0.1 5 0.8 3 m,且在组合中三口回灌井的综合作用下,成效水头值大多已连成包围性面状,在研究区北部集中的群井回灌中甚至产生了0.8 33.5 m的回升值,而在开采井中心,受大流量开采作用产生的负面水头成效最高达到了8.8 8 m,但是受外围三口回灌井的综合影响,其扩展面积有限。从图1 6(a)可以看出,在方案三条件下进行区域性供暖期回灌时,2 0 3 0年研究区底层温度场整体保持稳定,但是在各回灌点附近小范围的温度下降更加明显。为定量该方案回灌注水对研究区温度场的影响,利用G I S地理空间分析方法对地热模拟结果进行进一步处理,得到由回灌引起的底层温度下降带分布图1 6(b)。可以看出,由 回灌引起的 温 度 下 降 大 多 在“0.0 32.0 9 9”区间,且 普 遍 呈 由 回 灌 井 为 中心四周自“0.6 0 1 2.0 9 9”区 间 向“0.0 30.0 6 7”区间递减,受此影响区域上温度下降主要在00.0 6 7。“三灌一采”组合之间温度下降值分456矿 产 与 地 质2 0 2 3年 图1 4 2 0 3 0年热储层模拟压力场F i g.1 4 S i m u l a t e d p r e s s u r e f i e l d o f N e o g e n e t h e r m a lr e s e r v o i r i n 2 0 3 0图1 5 地热流体开采、回灌成效分布图F i g.1 5 D i s t r i b u t i o n m a p o f g e o t h e r m a l f l u i de x p l o i t a t i o n a n d r e i n j e c t i o n e f f e c t图1 6 研究区2 0 3 0年温度模拟结果F i g.1 6 T e m p e r a t u r e s i m u l a t i o n r e s u l t i n 2 0 3 0(a)顶层温度模拟 (b)由回灌引起的底层温度下降分布布特征明显,呈分别以三回灌点为中心向开采井递减,且开采井附近地温受此影响的下降值在“00.0 3”区间。在地热能方面,考虑开采水温及底线温度(2 5),计算得到2 0 2 02 0 3 0年按照一采三灌模式区内可增加净热能0.8 5 2 1 01 6J,比直接消耗式开采获得的热能高0.5 7 1 01 6J。即该回灌方案相较方案二更有助于增加深部大地热流传导补给被消耗的地热资源。5 结论通过对通许凸起中部研究区地热数模型的建立和不同地热开采方案的预测分析,获得以下几点认识:(1)单井开采模式下,研究区深层地热流体压力水头持续下降,且在通许县城井群集中区出现水头降落漏斗,区内温度场变化不明显,持续性开采可获得热能0.2 8 11 01 6J。556 第3 7卷 第3期 任静等:豫东通许凸起回灌条件下新近系地热资源优化开采方案研究(2)采灌模式下,回灌率3 8.4 1%时,研究区北部水头出现回升,地下水回灌引起的温度下降主要在0.0 5 9 0.1 1 4区间。回灌式开采可获得的热能0.5 9 61 01 6J,比 直 接 消 耗 式 开 采 获 得 的 热 能 高0.3 1 51 01 6J;回灌率1 0 0%时,全区北部水头回升更明显,地下水回灌引起的温度下降主要在0.0 32.0 9 9区间,回灌式开采可获得地热能0.8 5 21 01 6J,比直接消耗式开采获得的热能高0.5 71 01 6J。(3)研究结果表明一采三灌模式相对最优化。回灌能力越大,深部大地热流传导补给作用越明显,区内可利用热能越大。回灌动态补给可以消减研究区内地热流体压力水头下降趋势,虽然会引起回灌点附近温度下降越明显,但以回灌点为中心向开采井呈递减趋势,说明回灌点周边一定区域内出现热失衡现象,但全区温度场整体保持稳定。研究结果可为通许地区规模性开发利用地热资源提供理论依据。参考文献:1 王连成.天津市新近系馆陶组地热流体回灌研究D.北京:中国地质大学(北京),2 0 1 4.2 河南省人民政府.河南省地下水管理办法P.河南日报.2 0 2 2-1 2-3 1.3 河南省地质矿产勘查开发局第二地质环境调查院.河南省通许县岩溶热储地热资源调查调查报告R.郑州,2 0 2 2.4 韩文艳,代春艳,郑金栋.我国地热资源的可持续开发利用J.地球,2 0 1 3,6(1 1):1.5 朱红丽,刘小满,杨芳,等.开封市深层地热水回灌试验分析与研究J.河南理工大学学报(自然科学版),2 0 1 1,3 0(2):2 1 5-2 1 9.6 贺国平,邵景力,崔亚莉,等.F E F L OW在地下水流模拟方面的应用J.成都理工大学学报(自然科学版),2 0 0 3,3 0(4):3 5 6-3 6 1.7 翟美静.西安市三桥地区地热采灌系统示踪试验及数值模拟研究D.西安:长安大学.2 0 1 9.8 程万庆,刘九龙,陈海波.地热采灌对井回灌温度场的模拟研究J.世界地质(自然科学版),2 0 1 1,3 0(3):4 8 6-4 9 2.9 聂新良,王心义,尚照顺.河南东部平原地热系统的基本特征J.焦作工学院学报,2 0 0 2,2 1(2):1 0 1-1 0 4.1 0 阮传霞.天津地区雾迷山组热储地热回灌研究D.北京:中国地质大学(北京),2 0 1 8.1 1 雷海燕,朱家玲.滨海新区东营组地热热储回灌开采优化方案的研究J.太阳能学报,2 0 0 9,3 0(1 2):1 6 7 8-1 6 8 2.1 2 魏永霞,唐仲华,左霖.数值模拟在地热资源评价及优化开采方案的应用研究J.安徽地质,2 0 1 9,2 9(1):7 4-8 0.1 3 张进平,王煜曦,刘桂宏,等.通州区地热资源优化开采模式动态研究J.工程力学,2 0 2 2,3 9(6):2 4 7-2 5 6.1 4 宋前进,王刚,许一川,等.深层地热能开发及其对地热水流场的影响:以兰考县深层地热能开发为例J.煤田地质与勘探,2 0 2 1,4 9(2):1 8 4-1 9 3.1 5 邹鹏飞,王彩会,杜建国,等.地热水系统采灌方案模拟优化研究:以苏北农村清洁能源供示范区为例J.水文地质工程地质,2 0 2 3,3 9(6):1-1 8.S t u d y o n t h e o p t i m a l e x p l o i t a t i o n s c h e m e o f N e o g e n e g e o t h e r m a l r e s o u r c e su n d e r t h e r e i n j e c t i o n c o n d i t i o n i n T o n g x u U p l i f t,E a s t H e n a nR E N N i n g1,D E NG X i a o y i n g1,GUO Y a n2,3(1.E c o l o g i c a l E n v i r o n m e n t G e o l o g i c a l S e r v i c e C e n t e r o f H e n a n P r o v i n c i a l G e o l o g i c a l B u r e a u,Z h e n g z h o u 4 5 0 0 5 3,H e n a n,C h i n a;2.Y e l l o w R i v e r I n s t i t u t e o f H y d r a u l i c R e s e a r c h,Y e l l o w R i v e r C o n s e r v a n c y C o mm i s s i o n,Z h e n g z h o u 4 5 0 0 0 3,H e n a n,C h i n a;3.K e y L a b o r a t o r y o f L o w e r Y e l l o w R i v e r C h a n n e l a n d E s t u a r y R e g u l a t i o n,M i n i s t r y o f W a t e r R e s o u r c e s,Z h e n g z h o u 4 5 0 0 0 3,H e n a n,C h i n a)A b s t r a c t:T h e r e i s a g r e a t s i g n i f i c a n c e o f t h e o p t i m a l e x p l o i t a t i o n o f g e o t h e r m a l w e l l s f o r t h e i n t e n s i v e a n d e f f i c i e n t u t i l i z a t i o n o f g e o t h e r m a l r e s o u r c e s.T h i s p a p e r t a k e s t h e g e o t h e r m a l g e o l o g i c a l u n i t a t t h e m i d d l e o f T o n g x u U p l i f t i n t h e e a s t o f H e n a n P r o v i n c e a s t h e r e s e a r c h o b j e c t,F E F L OW s o f t w a r e i s u s e d t o e s t a b l i s h a t h r e e-d i m e n s i o n a l u n-s t e a d y f l o w h e a t c o n d u c t i o n m o d e l o f N e o g e n e t h e r m a l r e s e r v o i r,a n d t h e r e s u l t o f t h e r e i n j e c t i o n t e s t i s a p p l i e d t o p r e d i c t t h e c h a n g e o f t h e g e o t h e r m a l f l u i d p r e s s u r e f i e l d a n d t e m p e r a t u r e f i e l d i n t h e w h o l e r e g i o n f r o m 2 0 2 0 t o 2 0 3 0 u n d e r t h e c o n d i t i o n s o f o n l y m i n i n g m o d e,o n e w e l l m i n i n g a n d t w o w e l l i r r i g a t i o n(m i n i n g a n d i r r i