CN112444470B - 基于稳定流群孔抽水试验确定地下含水层渗透系数的方法 - Google Patents

基于稳定流群孔抽水试验确定地下含水层渗透系数的方法 Download PDF

Info

Publication number
CN112444470B
CN112444470B CN201910809716.9A CN201910809716A CN112444470B CN 112444470 B CN112444470 B CN 112444470B CN 201910809716 A CN201910809716 A CN 201910809716A CN 112444470 B CN112444470 B CN 112444470B
Authority
CN
China
Prior art keywords
aquifer
water
well
pumping
formula
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.)
Active
Application number
CN201910809716.9A
Other languages
English (en)
Other versions
CN112444470A (zh
Inventor
王俊智
李清波
万伟锋
曾峰
张海丰
蔡金龙
杜亚军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yellow River Engineering Consulting Co Ltd
Original Assignee
Yellow River Engineering Consulting Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yellow River Engineering Consulting Co Ltd filed Critical Yellow River Engineering Consulting Co Ltd
Priority to CN201910809716.9A priority Critical patent/CN112444470B/zh
Publication of CN112444470A publication Critical patent/CN112444470A/zh
Application granted granted Critical
Publication of CN112444470B publication Critical patent/CN112444470B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials

Landscapes

  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明公开了一种基于稳定流群孔抽水试验确定地下含水层渗透系数的方法,1,水文地质模型概化:使用现有水文地质勘察方法查明研究区的水文地质条件,并对研究区进行模型概化;查明地下水含水系统边界的位置、水文地质性质,查明地下水含水岩组的岩性及空间分布;2,实施稳定流群孔抽水试验:根据研究区水文地质勘察结果,布置实施稳定流群孔抽水试验;3,模型建立:将概化后的水文地质模型转换为数学模型;4,数学模型计算:根据含水层类型、含水系统边界条件的组合方式,计算含水层的渗透系数。本发明有助于准确获取研究区地下含水层的渗透系数,为地下水资源评价、水利工程渗漏计算、采掘工程隧道涌水预测提供科学依据和技术支持。

Description

基于稳定流群孔抽水试验确定地下含水层渗透系数的方法
技术领域
本发明涉及水文地质勘察领域,尤其是涉及基于稳定流群孔抽水试验确定地下含水层渗透系数的方法。
背景技术
抽水试验是水文地质勘察的重要技术手段,它不仅可以用来探明地下水系统的变化规律,还可较为准确地确定研究区的水文地质参数。大型群孔抽水试验是水文地质勘察中查明研究区水文地质条件最直接、保证程度最高的技术手段之一。利用大型群孔抽水试验可以确定诸如地下水补给来源、抽水影响范围、含水层的富水性及供水保证能力、地下水之间及与地表水体的水力联系、断裂构造的水力性质等问题。
《地下水动力学》教材和各水文地质勘察设计规范中所推荐的以泰斯和裘布依为代表的经典单孔井流公式是水文地质生产和科研中刻画地下水流场、确定水文地质参数的常用技术手段。但是,由于群孔在抽水过程中会相互干扰,上述经典单孔井流公式理论上无法使用。虽然单孔井流公式经过保角变换或虚拟井叠加修正后可以用于群孔抽水,但是在抽水井数量增多、流量不一,抽水井、观测井位置布置分散等情况下,修正后的井流公式不仅形式上复杂、臃肿,使用起来也容易出错。
随着计算机技术的发展,数值模拟方法服务于水文地质的生产和研究。通过建立地下水模型仿真群孔抽水试验,可以反演得到含水层的水文地质参数。但是,数值模拟方法属于黑箱模型,输入所需参数后直接输出计算结果,无法像井流公式一样直观揭示各水文地质参数之间的内在联系。此外,数值模拟法对相关人员的综合素质要求较高,数值、软件的价格也相对不菲。
发明内容
本发明目的是提供一种基于稳定流群孔抽水试验确定地下含水层渗透系数的方法。
为实现上述目的,本发明采取下述技术方案:
本发明所述基于稳定流群孔抽水试验确定地下含水层渗透系数的方法,包括下述步骤:
步骤1,水文地质模型概化:使用现有水文地质勘察方法查明研究区的水文地质条件,并对所述研究区进行模型概化;确定研究区的范围,查明地下水含水系统边界的位置、水文地质性质,查明地下水含水岩组的岩性及空间分布;
步骤2,实施稳定流群孔抽水试验:根据研究区水文地质勘察结果,布置实施稳定流群孔抽水试验;抽水井及观测井需要揭露水文地质单元,抽水井、观测井的数量、布置方式及抽水井的抽水量根据设计要求设定;试验过程中按照相关规范记录每个观测井的降深和抽水井的抽水量;
步骤3,模型建立:将概化后的水文地质模型转换为数学模型;数学模型的长度为研究区长度,数学模型的宽度为研究区宽度,含水层厚度按已有抽水井及观测井揭露含水层厚度的平均值控制;抽水井、观测井在数学模型中应为其局部坐标系坐标;
步骤4,数学模型计算:根据含水层类型、含水系统边界条件的组合方式,选择相对应的井流公式描述稳定流群孔抽水试验;将抽水井的抽水量及观测井的降深输入渗透系数的求参公式,计算含水层的渗透系数;所述含水层类型包括承压含水层、潜水含水层。
本发明建立稳定流群孔抽水试验的所述井流公式和渗透系数的求参公式步骤为:
步骤4.1,引入如下无量纲变量:
,公式 (1) ;
步骤4.2,对于所述承压含水层,基于“双重傅立叶变换”和“分离变量法”,群孔抽水试验稳定流井流公式的通解形式为:
,公式 (2) ;
经变形后,承压含水层渗透系数的通解形式为:
,公式 (3) ;
式中,s表示降深,T=H·K表示承压含水层的导水系数,H表示承压含水层的厚度,K表示承压含水层的渗透系数,w表示抽(注)水井的数量,Q i 表示第i口抽(注)水井的流量,当Q i <0时表示注水,当Q i >0时表示抽水;x i y i 表示第i口抽(注)水井的位置,xy表示观测井的位置,L x L y 表示含水层的长和宽,表示井函数;
步骤4.3,对于所述潜水含水层,在裘布依假设的前提条件下,使用“双重傅立叶变换”和“分离变量法”,群孔抽水试验稳定流井流公式的通解形式为:
,公式(4) ;
式中,h表示潜水含水层的厚度,即潜水面相对于隔水底板的高程差,h 0表示其初始值;
当抽水引起的降深远小于潜水含水层的厚度时,即时,
, 公式(5) ;
经变形后,潜水含水层渗透系数的表达式为:
, 公式(6)。
本发明所述潜水含水层与所述承压含水层的井函数表达式一致。对于所述潜水含水层或所述承压含水层的边界条件组合为四个定水头边界的情形,井函数为:
,公式(7);
对于潜水含水层或承压含水层的边界条件组合为一个隔水边界、三个定水头边界的情形,井函数为:
,公式(8);
对于潜水含水层或承压含水层的边界条件组合为两个隔水边界、两个定水头边界(边界相对)的情形,井函数为:
, 公式(9);
其中,
对于潜水含水层或承压含水层的边界条件组合为两个隔水边界、两个定水头边界(边界相邻)的情形,井函数为:
,公式(10);
对于潜水含水层或承压含水层的边界条件组合为一个定水头边界、三个隔水边界的情形,井函数为:
,
公式 (11);
其中,
本发明优点体现在以下方面:
1、现有单孔井流公式在抽水井数量增多、流量不一、抽水井位置布置分散等情况下,修正后的井流公式不仅形式上复杂、臃肿,使用起来也容易出错。相比较而言,本发明方法更为实用、简洁、准确。
2、相较于地下水数值模型,本发明方法所确定的数学模型直观揭示了各水文地质参数之间的内在联系,对相关人员的综合素质要求低。
3、本发明方法有助于准确获取研究区地下水含水层的渗透系数,为地下水资源评价、地下水环境影响评价、水利工程渗漏计算、市政工程基坑降水方案设计、采掘工程隧道涌水预测等提供科学依据和技术支持。
附图说明
图1是本发明水文地质概念模型的示意图。
图2是本发明所述潜水含水层或承压含水层的边界条件组合为四个定水头边界的示意图。
图3是本发明所述潜水含水层或承压含水层的边界条件组合为三个定水头边界、一个隔水边界的示意图。
图4是本发明所述潜水含水层或承压含水层的边界条件组合为两个定水头边界、两个隔水边界的示意图(定水头边界位置相对)。
图5是本发明所述潜水含水层或承压含水层的边界条件组合为两个隔水边界、两个定水头边界的示意图(定水头边界位置相邻)。
图6是本发明所述潜水含水层或承压含水层的边界条件组个为一个定水头边界、三个隔水边界的示意图。
具体实施方式
下面结合附图对本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述实施例。
本发明所述基于稳定流群孔抽水试验确定地下含水层渗透系数的方法,包括下述步骤:
步骤1,水文地质模型概化:使用现有水文地质勘察方法查明研究区的水文地质条件,并对研究区进行模型概化;确定研究区的范围,查明地下水含水系统边界的位置,查明水文地质性质及地下水含水岩组的岩性及空间分布;
步骤2,实施稳定流群孔抽水试验:根据研究区水文地质勘察结果,布置实施稳定流群孔抽水试验;抽水井1及观测井2用于揭露水文地质单元,抽水井1、观测井2的数量、布置方式及抽水井1的抽水量Q i根据设计要求设定;试验过程中按照相关规范记录观测井2的降深s和抽水井1的抽水量Q i
步骤3,按图1所示将概化后的水文地质模型转化为数学模型:数学模型的长度为研究区的长度L x ,数学模型的宽度为研究区的宽度L y ,含水层厚度(承压含水层厚度为H,潜水含水层厚度为h 0)按已有抽水井1揭露含水层厚度的平均值控制;抽水井1、观测井2在数学模型中应为其局部坐标系坐标;抽水井1的位置用(x i y i )坐标表示,观测井2的位置用坐标(xy)表示;
步骤4,数学模型计算:根据含水层类型(承压含水层或潜水含水层)、含水系统边界条件的组合方式,选择相对应的井流公式描述稳定流群孔抽水试验;将抽水井1的抽水量Q i及观测井2的降深s及其空间坐标(抽水井:x i y i ;观测井:xy)输入渗透系数的求参公式,计算含水层的渗透系数K
本发明建立群孔抽水试验的所述井流公式和渗透系数的求参公式步骤为:
步骤4.1,引入如下无量纲变量:
公式 (1) ;
步骤4.2,对于矩形承压含水层,基于“双重傅立叶变换”和“分离变量法”,群孔抽水试验稳定流井流公式的通解形式为:
, 公式 (2) ;
经变形后,矩形承压含水层渗透系数的通解形式为
, 公式 (3) ;
式中,s表示降深,T=H·K表示承压含水层的导水系数,H表示承压含水层的厚度,K表示承压含水层的渗透系数,w表示抽水井1的数量,Q i 表示第i口抽水井1的流量,当Q i <0时表示注水,当Q i >0时表示抽水;x i y i 表示第i口抽水井1的位置,xy表示观测井的位置,L x L y 表示含水层的长和宽,表示井函数;
步骤4.3,对于矩形潜水含水层,在裘布依假设的前提条件下,基于“双重傅立叶变换”和“分离变量法”,群孔抽水试验稳定流井流公式的通解形式为:
, 公式(4) ;
式中,h表示潜水含水层的厚度,即潜水面相对于隔水底板的高程差,h 0表示其初始值;
当抽水引起的降深远小于潜水含水层的厚度时,即时,
, 公式(5) ;
经变形后,潜水含水层渗透系数的表达式为:
, 公式(6)。
本发明矩形潜水含水层与矩形承压含水层的井函数表达式一致,对于矩形潜水含水层或矩形承压含水层的边界条件组合为四个定水头边界3的情形,如图2所示,井函数为:
,公式(7);
对于矩形潜水含水层或矩形承压含水层的边界条件组合为一个隔水边界4、三个定水头边界3的情形,如图3所示,井函数为:
,公式(8);
对于矩形潜水含水层或矩形承压含水层的边界条件组合为两个隔水边界4、两个定水头边界3(两个定水头边界3位置相对)的情形,如图4所示,井函数为:
, 公式(9);
其中,
对于矩形潜水含水层或矩形承压含水层的边界条件组合为两个隔水边界4、两个定水头边界3(两个定水头边界3位置相邻)的情形,如图5所示,井函数为:
,公式(10);
对于矩形潜水含水层或矩形承压含水层的边界条件组合为一个定水头边界3、三个隔水边界4的情形,如图6所示,井函数为:
,
公式(11);
其中,

Claims (1)

1.一种基于稳定流群孔抽水试验确定地下含水层渗透系数的方法,其特征在于:包括下述步骤:
步骤1,水文地质模型概化:使用现有水文地质勘察方法查明研究区的水文地质条件,并对所述研究区进行模型概化;确定研究区的范围,查明地下水含水系统边界的位置、水文地质性质,查明地下水含水岩组的岩性及空间分布;
步骤2,实施稳定流群孔抽水试验:根据研究区水文地质勘察结果,布置实施稳定流群孔抽水试验;抽水井及观测井需要揭露水文地质单元,抽水井、观测井的数量、布置方式及抽水井的抽水量根据设计要求设定;试验过程中按照相关规范记录每个观测井的降深和抽水井的抽水量;
步骤3,模型建立:将概化后的水文地质模型转换为数学模型;所述数学模型的长度为研究区长度,数学模型的宽度为研究区宽度,含水层厚度按已有抽水井揭露含水层厚度的平均值控制;抽水井、观测井在数学模型中应为其局部坐标系坐标;
步骤4,数学模型计算:根据含水层类型、含水系统边界条件的组合方式,选择相对应的井流公式描述稳定流群孔抽水试验;将抽水井的抽水量及观测井的降深输入渗透系数的求参公式,计算含水层的渗透系数;所述含水层类型包括承压含水层、潜水含水层;
其中,建立稳定流群孔抽水试验的所述井流公式和渗透系数的求参公式步骤为:
步骤4.1,引入如下无量纲变量:
,公式 (1);
步骤4.2,对于所述承压含水层,基于“双重傅立叶变换”和“分离变量法”,群孔抽水试验稳定流井流公式的通解形式为:
,公式 (2);
经变形后,承压含水层渗透系数的通解形式为:
,公式 (3) ;
式中,s表示降深,T=H·K表示承压含水层的导水系数,H表示承压含水层的厚度,K表示承压含水层的渗透系数,w表示抽水井或注水井的数量,Q i 表示第i口抽水井或注水井的流量,当Q i <0时表示注水,当Q i >0时表示抽水;x i y i 表示第i口抽水井或注水井的位置,xy表示观测井的位置,L x L y 表示含水层的长度和宽度,表示井函数;
步骤4.3,对于所述潜水含水层,在裘布依假设的前提条件下,基于“双重傅立叶变换”和“分离变量法”,群孔抽水试验稳定流井流公式的通解形式为:
,公式(4) ;
式中,h表示潜水含水层的厚度,即潜水面相对于隔水底板的高程差,h 0表示其初始值;
当抽水引起的降深远小于潜水含水层的厚度时,即s≪ h0时,
, 公式(5) ;
经变形后,潜水含水层渗透系数的表达式为:
, 公式(6);
所述潜水含水层与所述承压含水层的井函数表达式一致;对于所述潜水含水层或所述承压含水层的边界条件组合为四个定水头边界的情形,井函数为:
,公式(7);
对于潜水含水层或承压含水层的边界条件组合为一个隔水边界、三个定水头边界的情形,井函数为:
,公式(8);
对于潜水含水层或承压含水层的边界条件组合为两个隔水边界、两个定水头边界且定水头边界位置相对的情形,井函数为:
,公式(9);
其中,
对于潜水含水层或承压含水层的边界条件组合为两个隔水边界、两个定水头边界且定水头边界位置相邻的情形,井函数为:
,公式(10);
对于潜水含水层或承压含水层的边界条件组合为一个定水头边界、三个隔水边界的情形,井函数为:
,公式(11);
其中,
CN201910809716.9A 2019-08-29 2019-08-29 基于稳定流群孔抽水试验确定地下含水层渗透系数的方法 Active CN112444470B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910809716.9A CN112444470B (zh) 2019-08-29 2019-08-29 基于稳定流群孔抽水试验确定地下含水层渗透系数的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910809716.9A CN112444470B (zh) 2019-08-29 2019-08-29 基于稳定流群孔抽水试验确定地下含水层渗透系数的方法

Publications (2)

Publication Number Publication Date
CN112444470A CN112444470A (zh) 2021-03-05
CN112444470B true CN112444470B (zh) 2023-09-26

Family

ID=74741553

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910809716.9A Active CN112444470B (zh) 2019-08-29 2019-08-29 基于稳定流群孔抽水试验确定地下含水层渗透系数的方法

Country Status (1)

Country Link
CN (1) CN112444470B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113202459B (zh) * 2021-06-11 2022-02-15 中国矿业大学 一种利用钻孔自然伽马测井预测断层渗透性的方法
CN114487347B (zh) * 2022-01-24 2022-10-25 河海大学 一种识别钻孔正薄壁效应并确定含水层水文地质参数的微水试验法
CN114660667A (zh) * 2022-04-08 2022-06-24 青海省环境地质勘查局 一种山前冲洪积扇大型水源地勘探方法
CN116777229B (zh) * 2023-05-10 2024-02-27 南方科技大学 一种地下水位预测方法、系统、终端及存储介质
CN116341280B (zh) * 2023-05-11 2023-08-18 中国地质大学(北京) 确定隔水层不稳定的多含水层分层渗透系数的方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108257488A (zh) * 2018-03-22 2018-07-06 广西大学 不稳定承压-无压井流模型及承压含水层参数反演方法
CN109242312A (zh) * 2018-09-11 2019-01-18 中国矿业大学(北京) 一种单井循环系统承压含水层渗流解析解的求解方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108257488A (zh) * 2018-03-22 2018-07-06 广西大学 不稳定承压-无压井流模型及承压含水层参数反演方法
CN109242312A (zh) * 2018-09-11 2019-01-18 中国矿业大学(北京) 一种单井循环系统承压含水层渗流解析解的求解方法

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
Chunhui Lu等.Steady state analytical solutions for pumping in a fully bounded rectangular aquifer .《Water Resources Research》.2015,全文. *
In Wook Yeo等.Analytical solution for arbitrarily located multiwells in an anisotropic homogeneous confined aquifer.《WATER RESOURCES RESEARCH》.2003,第第39卷卷(第第5期期),全文. *
Jun-Zhi Wang等.An analytical study on groundwater flow in drainage basins with horizontal wells.《Hydrogeology Journal》.2014,全文. *
史骥 ; 张建华 ; 夏岸雄 ; 周勇 ; .越流补给群孔抽水试验的含水层参数计算及应用.金属矿山.2016,(第04期),全文. *
周子俣 ; 秦志泉 ; .单孔稳定流抽水试验计算渗透系数对比研究.矿产勘查.2018,(第04期),全文. *
周志芳,汤瑞凉,汪斌.基于抽水试验资料确定含水层水文地质参数.河海大学学报(自然科学版).1999,(第03期),全文. *
李霞等.基于解析法和数值法的非稳定流抽水试验参数反演.《地球科学》.2017,第第 42 卷卷(第第 5 期期),全文. *
秦建甫等.群孔干扰抽水试验在哈头才当水源地的应用.《资源环境与工程》.2012,第第 26 卷卷(第第 5 期期),摘要和第1~2节. *
董一鹏 ; .煤矿水文地质勘探中群孔抽水试验的应用研究.能源与节能.2016,(第02期),全文. *

Also Published As

Publication number Publication date
CN112444470A (zh) 2021-03-05

Similar Documents

Publication Publication Date Title
CN112444470B (zh) 基于稳定流群孔抽水试验确定地下含水层渗透系数的方法
Hu et al. Simulation study on water quality based on sediment release flume experiment in Lake Taihu, China
Wang et al. Modeling of sediment and heavy metal transport in Taihu Lake, China
CN108241777A (zh) 基于非结构网格有限元法计算水合物沉积物中渗流速度场的方法
Shi et al. A statistical model of deformation during the construction of a concrete face rockfill dam
CN108301823A (zh) 一种识别储层油气甜点的方法
CN111898257B (zh) 区域暗管布局及排水排盐数值模拟方法和装置
CN111723505A (zh) 一种流域水质水量监测系统
CN112007943B (zh) 原位注入高压旋喷注射修复地下水污染注入药剂扩散半径的确定方法
CN106446438B (zh) 一种电厂温排水深水排放的近远区耦合数值模拟方法
CN115203945A (zh) 一种水文地质参数估计方法、装置、设备和存储介质
CN105910587A (zh) 一种基于潮汐参数反演的潮汐预测方法
CN106930749A (zh) 基于降压的海域天然气水合物层钻井等效渗透率计算方法
CN102720487B (zh) 一种获取气藏地层压力系统及其方法
CN116467958A (zh) 盐湖卤水数值模型构建与补水溶矿效率计算方法
Wang et al. Variable diffusion boundary layer and diffusion flux at sediment-water interface in response to dynamic forcing over an intertidal mudflat
CN111749664B (zh) 一种油田全程化学驱组合接替时机优化方法和系统
CN116432551B (zh) 基于卫星图像和cfd的河流交汇区水质分布仿真方法
CN103615239B (zh) 一种测井地层成像系统及其方法
Zhang et al. Analytical solutions for seawater intrusion and the maximum pumping rate of a well in unconfined coastal aquifers bounded by L-shaped coastlines
You et al. On improvement of water quality of a reservoir by optimizing water exchange
Wang et al. Scaling relationships for diffusive boundary layer thickness and diffusive flux based on in situ measurements in coastal seas
Miao et al. Research on seepage control of Jurong pumped storage hydroelectric power station
CN115563746A (zh) 一种非均质三维水文地质模型的动态构建方法
CN115017701A (zh) 一种南海北部三维物理-生态耦合模型的构建方法及其应用

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant