豫东通许凸起回灌条件下新近系地热资源优化开采方案研究.pdf
《豫东通许凸起回灌条件下新近系地热资源优化开采方案研究.pdf》由会员分享,可在线阅读,更多相关《豫东通许凸起回灌条件下新近系地热资源优化开采方案研究.pdf(11页珍藏版)》请在咨信网上搜索。
1、地热井的优化开采对地热资源的集约高效利用意义重大。本文以豫东通许凸起中部地热地质单元为研究对象,采用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。在此基础上提出了最优化的采灌布井模式,为通
2、许凸起深部地热资源可持续开发提供理论依据。关键词:回灌;模型;新近系热储;温度场;压力场;可采热能中图分类号: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。对于地热供暖项目,其工作原理主要是提取地热水中的可采热能,开采回
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
4、 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年度省财政地质勘查项目
5、“河南省通许县岩溶热储地热资源调查”(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)、寒武-奥陶系(
6、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
7、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井揭露地层情况,新近系明化镇组由浅黄色、黄褐色粉砂岩、中细砂岩、细砂岩与棕红、紫红色、浅棕红色泥岩互层组成;新近系
8、馆陶组岩性为紫红色泥岩与浅黄色、黄褐色中细砂岩、细砂岩互层,与下伏地层呈不整合接触。两组岩层无显著隔水边界存在。全区新近系热储层有砂(砂岩)层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 新近系热储回灌
9、能力确定本次利用区内最新施工的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.
10、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
11、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所示。垂向上,新近系上部由于第四系岩层的覆盖,不考虑
12、大气降水入渗补给量,加之图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零通量边
13、界 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
14、 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
15、,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 y
16、Ty+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
17、)模型赋值在水文地质概念模型中,采用给定水头边界和流量边界来进行模型区含水层侧向补给与排泄的赋值。其中给定水头边界由软件根据周边压(力)头流场分布来自动计算边界地下水的补排量;流量边界所赋数值由人工计算而来,具体见表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 注:正
18、值表示排泄流量,负值表示补给流量。(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年
19、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水文地
20、质参数参数值参数单位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)模拟校正检验研究区新近系热
21、储的水温动态在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 地热资源开采方案预测在已建立地热模型的基础上,结合研
22、究区开采现状及规划,并在考虑供暖期、停暖期的情况下,进行水热耦合模拟,研究地热资源不同开发利用情景对整个区域压力场和温度场的影响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
23、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%。在研究区现有的各
24、地热开采井附近设置三个回灌虚拟井点,回灌点至开采点的距离在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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 豫东 通许 凸起 条件下 新近 地热资源 优化 开采 方案 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。