CN105184035A - 一种确定地下热水可开采量的开采井优化布局方法 - Google Patents
一种确定地下热水可开采量的开采井优化布局方法 Download PDFInfo
- Publication number
- CN105184035A CN105184035A CN201510186652.3A CN201510186652A CN105184035A CN 105184035 A CN105184035 A CN 105184035A CN 201510186652 A CN201510186652 A CN 201510186652A CN 105184035 A CN105184035 A CN 105184035A
- Authority
- CN
- China
- Prior art keywords
- unit
- well
- exploitation
- max
- water
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Pipeline Systems (AREA)
Abstract
本发明公开了一种确定地下热水可开采量的地热井优化布局方法,包括有如下步骤:给定计算区面积F、开采井井管半径rω和地下热水开采年限t,根据计算区水文地质和经济技术条件给定设计水位降深Smax;依据单井非稳定流抽水试验资料确定导水系数T及弹性释水系数μ*;根据导水系数T、弹性释水系数μ*,选用绘制设计水位降深时单位可开采量Q采/Smax和地下热水开采井数n的关系曲线;选取(Q采/Smax)~n关系曲线变化陡峭和平缓的分界点,确定优化布局井数n采和单位可开采量Q采/Smax;根据优化布局井数n采、单位可开采量Q采/Smax、设计水位降深Smax,确定地下热水可开采量并应用于工程实践。
Description
技术领域
本发明涉及一种确定地下热水可开采量的开采井优化布局方法,是一种自行设计的保证同样设计降深和相同开采井数条件下,计算地下热水可开采量的优化布井法。根据地下热储含水层的导水系数、弹性释水系数、单井开采量、开采井影响半径、设计水位降深等参数,选择优化布井方案,实现自动确定地下热水可开采量的要求。本方法适用于大面积分布、顶底板水平、层间水力联系不密切、均质各向同性的孔隙型地下热水承压含水层。
背景技术
孔隙型地下热水作为一种洁净且具有“温、矿、水”综合特征的资源,其设计降深条件下可开采量的确定对于地下热水的开发利用及保护具有重要意义。传统的确定地下热水可开采量的均匀布井方法,在达到设计水位降深时由于单井开采量系人为给定,故确定的开采井数也表现出一定的随意性,所确定的地下热水可开采量受人为影响严重。
一种确定地下热水可开采量的开采井优化布局方法可以较好地解决上述问题,以均匀布井法为基础,依据计算区面积、开采井管半径、开采年限,综合考虑地下热储含水层的导水系数、弹性释水系数、设计水位降深等因素,引入井数与单位开采量的概念,达到优化布局开采井并自动确定地下热水可开采量的目的。
发明内容
本发明针对传统的均匀布井方法实际应用过程的随意性,提出了自行优化确定地下热水可开采量的实践方法。该方法应用地下热水含水层的导水系数、弹性释水系数、计算区面积、开采井管半径、开采年限等参数,在满足设计水位降深的基础上,实现了优化布局开采井条件下地下热水可开采量的确定。
本发明提供了一种确定地下热水可开采量的地热井优化布局方法。包括设计水位降深、导水系数、弹性释水系数、计算点水位降深、单位可开采量Q采/Smax、优化布局井数n采、地下热水可开采量确定。具体步骤如下:
步骤一,给定计算区面积F(单位为km2)、开采井管半径rω(单位为m)和地热水开采年限t(单位为d);根据计算区水文地质和经济技术条件给定设计水位降深Smax(单位为m);
步骤二,依据单井非稳定流抽水试验资料确定地下热水含水层的导水系数T(单位为m2/d)及弹性释水系数μ*;
步骤三,根据导水系数T、弹性释水系数μ*,选用承压完整井定流量非稳定流公式确定计算点的水位降深(单位为m),式中S为抽水影响范围内,任一点任一时刻的水位降深,m;Q为单个抽水井开采量,m3/d;r为计算点到开采井的距离(单位为m),当计算点位于开采井时取开采井管半径rω,t为地热水开采年限;
步骤四,设计计算区有n眼地下热水开采井,确定某计算点水位降深的公式为
步骤五,计算区所有水井的开采量相同时,某计算点的单位开采量计算公式为计算区地下热水可开采量为Q采=n×Q,其中n为地下热水开采井数;
步骤六,地下热水开采期末计算区中心井达到设计水位降深Smax时,单位可开采量计算公式为绘制地下热水含水层的(Q采/Smax)~n关系曲线,其中单位可开采量Q采/Smax的单位为(m3/d)/m;
步骤七,根据地下热水含水层的(Q采/Smax)~n关系曲线,选取随地下热水开采井数n增加单位可开采量Q采/Smax不再显著增加的点(即(Q采/Smax)~n关系曲线变化陡峭和平缓的分界点),相应的地下热水开采井数n和单位可开采量Q采/Smax为开采井优选布局条件下的优化布局井数n采和相应的单位可开采量Q采/Smax;
步骤八,根据优化布局井数n采、单位可开采量Q采/Smax、设计水位降深Smax,确定地下热水可开采量Q采(单位为m3/d)并应用于工程实践。
附图说明
图1是地下热水开采井正方形布局方案示意图,图中座标原点位于计算区中心的井上,横排沿X轴方向,以i=-M,…-2,-1,0,1,2,…,M编号;纵排沿y轴方向,以j=-K,…,-2,-1,0,1,2,…,K编号;任一井的编号为(i,j),任一井到中心井的距离式中r0为井排间距(单位为m),井排间距中心井的编号为(0,0),中心井r取开采井管半径rω;,F为计算区面积,n为地下热水开采井数(单位为眼);
图2是确定800~1000m地下热水含水层设计水位降深时单位可开采量Q采/Smax和地下热水开采井数n的关系曲线,图中纵坐标为单位可开采量Q采/Smax,横坐标为地下热水开采井数n;
图3是确定1000~1200m地下热水含水层设计水位降深时单位可开采量Q采/Smax和地下热水开采井数n的关系曲线,图中纵坐标为单位可开采量Q采/Smax,横坐标为地下热水开采井数n。
具体实施方式
本发明提供了一种确定地下热水可开采量的地热井优化布局方法。
下面以开封市主城区为计算区的800~1000m、1000~1200m地下热水含水层为例进行进一步说明。具体步骤如下:
(1)计算区面积165.23km2,地下热水开采井管半径0.08m;根据当地水文地质和经济技术条件,计算区的800~1000m、1000~1200m地下热水含水层设计水位降深50m,开采年限30年;
(2)根据单井非稳定流抽水试验资料,确定导水系数T和弹性释水系数μ*,见表1;地下热水开采井布局方案为正方形,见图1;
表1例区计算参数
(3)将已知的导水系数T、弹性释水系数μ*、开采期限t、某一开采井到中心井的距离r代入公式绘制(Q采/Smax)~n的关系曲线,见图2与图3;
(4)根据图2与图3选择曲线变化陡峭和平缓的分界点,得到优化布局井数n采、相应的井排间距r0和单位可开采量Q采/Smax,见表2;
表2开采井优化布局方案参数
(5)根据优化布局井数n采(井排间距r0)、单位可开采量Q采/Smax和设计水位降深Smax,确定地下热水可开采量并应用于工程实践。
Claims (5)
1.一种确定地下热水可开采量的开采井优化布局方法,其特征在于,确定地下热水可开采量的步骤如下:
步骤一,给定计算区面积F(单位为km2)、开采井井管半径rω(单位为m)和地热水开采年限t(单位为d);根据计算区水文地质和经济技术条件给定设计水位降深Smax(单位为m);
步骤二,依据单井非稳定流抽水试验资料确定地下热水含水层的导水系数T(单位为m2/d)及弹性释水系数μ*;
步骤三,根据导水系数T、弹性释水系数μ*,选用承压完整井定流量非稳定流公式确定计算点的水位降深(单位为m),式中的S为抽水影响范围内,任一点任一时刻的水位降深,m;Q为单个抽水井开采量,m3/d; r为计算点到开采井的距离(单位为m),当计算点位于开采井时取开采井管半径rω,t为地热水开采年限;
步骤四,设计计算区有n眼地下热水开采井,确定某计算点水位降深的公式为
步骤五,计算区所有水井的开采量相同时,某计算点单位开采量计算公式为计算区地下热水可开采量为Q采=n×Q;
步骤六,地下热水开采期末计算区中心井达到设计水位降深Smax时,单位可开采量计算公式为绘制地下热水含水层的(Q采/Smax)~n关系曲线,其中单位可开采量Q采/Smax的单位为(m3/d)/m;
步骤七,根据(Q采/Smax)~n关系曲线,选取随地下热水开采井数n增加Q采/Smax不再显著增加的点(即(Q采/Smax)~n关系曲线变化陡峭和平缓的分界点),相应的n和Q采/Smax为开采井优选布局条件下的优化布局井数n采和相应的单位可开采量Q采/Smax;
步骤八,根据优化布局井数n、单位可开采量Q采/Smax、设计水位降深Smax,确定地下热水可开采量(单位为m3/d)并应用于工程实践。
2.根据权利要求1所述的一种确定地下热水可开采量的开采井优化布局方法,其特征在于,在步骤一中,设计水位降深Smax需根据计算区水文地质和经济技术条件确定。
3.根据权利要求1所述的一种确定地下热水可开采量的开采井优化布局方法,其特征在于,在步骤二中,地下热水含水层的导水系数T(m2/d)及弹性释水系数μ*应依据单井非稳定流抽水试验资料确定。
4.根据权利要求1所述的一种确定地下热水可开采量的开采井优化布局方法,其特征在于,在步骤三中,式中的 其中:Q为单个抽水井涌水量(m3),S为水位降深(m),r为计算点到开采井的距离(m),当计算点位于开采井时取开采井管半径rω(m),t为地热水开采年限(d)。
5.根据权利要求1所述的一种确定地下热水可开采量的开采井优化布局方法,其特征在于,在步骤七中,(Q采/Smax)~n关系曲线的变化陡峭和平缓的分界点对应的n和Q采/Smax为优化布局井数和相应的单位可开采量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510186652.3A CN105184035A (zh) | 2015-04-20 | 2015-04-20 | 一种确定地下热水可开采量的开采井优化布局方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510186652.3A CN105184035A (zh) | 2015-04-20 | 2015-04-20 | 一种确定地下热水可开采量的开采井优化布局方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105184035A true CN105184035A (zh) | 2015-12-23 |
Family
ID=54906114
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510186652.3A Pending CN105184035A (zh) | 2015-04-20 | 2015-04-20 | 一种确定地下热水可开采量的开采井优化布局方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105184035A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106874560A (zh) * | 2017-01-17 | 2017-06-20 | 北京师范大学 | 一种傍河水源地优化布井方法 |
CN107180304A (zh) * | 2017-05-15 | 2017-09-19 | 山东省水利勘测设计院 | 一种基于袭夺量的沿河水源地评价方法 |
CN112146845A (zh) * | 2020-09-27 | 2020-12-29 | 武汉大学 | 地下水瞬时弹性释水系数和完全弹性释水系数的测试方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1273351A (zh) * | 1999-05-05 | 2000-11-15 | 李凤仪 | 地热区深层和非地热区超深层地热能源的开采利用方法 |
CN201965744U (zh) * | 2010-12-01 | 2011-09-07 | 天津地热勘查开发设计院 | 地热井监测数据远程采集装置 |
CN202562529U (zh) * | 2012-05-24 | 2012-11-28 | 天津地热勘查开发设计院 | 地热井水位-温度数字测量*** |
-
2015
- 2015-04-20 CN CN201510186652.3A patent/CN105184035A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1273351A (zh) * | 1999-05-05 | 2000-11-15 | 李凤仪 | 地热区深层和非地热区超深层地热能源的开采利用方法 |
CN201965744U (zh) * | 2010-12-01 | 2011-09-07 | 天津地热勘查开发设计院 | 地热井监测数据远程采集装置 |
CN202562529U (zh) * | 2012-05-24 | 2012-11-28 | 天津地热勘查开发设计院 | 地热井水位-温度数字测量*** |
Non-Patent Citations (3)
Title |
---|
刘小满, 王心义: "开封市区地热水可开采量评价", 《西部探矿工程》 * |
施凤英,王宇辉: "地下水开采量的合理设计探讨", 《山东冶金》 * |
王兴: "地下热水生产井与回灌井间距及其影响因素", 《地下水》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106874560A (zh) * | 2017-01-17 | 2017-06-20 | 北京师范大学 | 一种傍河水源地优化布井方法 |
CN106874560B (zh) * | 2017-01-17 | 2019-11-12 | 北京师范大学 | 一种傍河水源地优化布井方法 |
CN107180304A (zh) * | 2017-05-15 | 2017-09-19 | 山东省水利勘测设计院 | 一种基于袭夺量的沿河水源地评价方法 |
CN112146845A (zh) * | 2020-09-27 | 2020-12-29 | 武汉大学 | 地下水瞬时弹性释水系数和完全弹性释水系数的测试方法 |
CN112146845B (zh) * | 2020-09-27 | 2021-11-02 | 武汉大学 | 地下水瞬时弹性释水系数和完全弹性释水系数的测试方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sommer et al. | Thermal performance and heat transport in aquifer thermal energy storage | |
CN105808793B (zh) | 一种基于非结构网格的水平井分段压裂数值模拟方法 | |
Willems et al. | The impact of reduction of doublet well spacing on the Net Present Value and the life time of fluvial Hot Sedimentary Aquifer doublets | |
Jiang et al. | Heat extraction of novel underground well pattern systems for geothermal energy exploitation | |
CN104265270B (zh) | 一种碳酸盐岩缝洞钻井轨迹设计及控制方法 | |
LaMoreaux et al. | New development in theory and practice in mine water control in China | |
CN105184035A (zh) | 一种确定地下热水可开采量的开采井优化布局方法 | |
Surinaidu et al. | Assessment of groundwater inflows into Kuteshwar Limestone Mines through flow modeling study, Madhya Pradesh, India | |
González-Nicolás et al. | Pressure management via brine extraction in geological CO2 storage: Adaptive optimization strategies under poorly characterized reservoir conditions | |
CN102682195B (zh) | 半潜式平台瞬态钻井井筒温度计算方法 | |
Al-Zyoud et al. | Dynamic numerical modeling of the usage of groundwater for cooling in north east Jordan–A geothermal case study | |
CN102587887A (zh) | 一种气井井况预测方法 | |
Choe et al. | Method of simulation and estimation of SCW system considering hydrogeological conditions of aquifer | |
Somogyi et al. | Thermal impact assessment with hydrodynamics and transport modeling | |
Al-Turki et al. | On SAGD in oil sands reservoirs with no caprock and top water zone | |
Meng et al. | Coupled simulation–optimization model for draining confined aquifer via underground boreholes to prevent water inrush of coal mines | |
CN104564069A (zh) | 一种基于方格网法的地面动态沉陷预测与复垦方法 | |
CN104265185B (zh) | 一种三级油管柱设计方法及装置 | |
CN104533519A (zh) | 立井井筒通过强含水厚岩层时涌水水害的治理方法 | |
WANG et al. | Study on influence of freezing tubes layout on development law of temperature field of freezing wall | |
Zhang et al. | On Tunable Facture Conductivity to Improve Heat Extraction from Fractured Geothermal Systems | |
CN111119992A (zh) | 一种煤层顶板疏放水钻孔参数确定方法 | |
CN104563971B (zh) | 采矿***及资源重叠矿区的采矿方法 | |
CN214499020U (zh) | 水热型地热井工厂 | |
CN113982589B (zh) | 一种富油煤原位开采的温度控制方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151223 |