试验方法确定水文地质参数
水文地质试验

水文地质试验为测定水文地质参数和了解地下水的运动规律而进行的试验工作,内容包括抽水、注水、压水、渗水、连通、流速和弥散系数测定等。
其中抽水试验是最主要的手段。
抽水试验利用井(孔)抽取地下水,以了解井的涌水量及其与水位下降的历时变化关系。
抽水试验按地下水流态可分为稳定流和非稳定流抽水。
按抽水井与观测孔的关系可分为单孔抽水和多孔抽水;按井孔贯穿含水层的程度可分为完整井抽水和非完整井抽水;按抽水井与含水层的关系可分为分层抽水和混合抽水等。
①稳定流抽水。
抽水时流量和水位降同时保持不变,适用于抽水量小于补给量的地区,这种抽水一般需进行三次水位降。
其最大降深值,潜水应介于其含水层厚度的1/3~1/2之间;承压水不得大于其承压水头。
稳定时间一般为8~24小时当水质和水量发生突然变化时则要延长稳定时间。
②非稳定流抽水。
保持抽水量为常量,观测水位随时间的变化,在抽水量大于补给量或抽水过程中水位一直持续下降的地区更为适用。
抽水时间视其目的、水文地质特征、水位降与时间关系曲线类型和选用计算参数的公式而定。
一般为12~24小时。
稳定流与非稳定流抽水可结合进行,观测孔兼顾两者的计算要求布设,既满足后者对水量、水位的观测精度,又达到前者的延续时间,互相校正,以获得较理想的成果。
抽水试验的设备通常为空气压缩机或深井泵。
当地下水最大动水位深度小于7.5米时,可采用卧式离心泵。
若是非稳定流抽水,则宜采用电动离心泵或深井泵。
抽水试验过程中,为便于发现和及时处理异常现象,确定抽水试验延续时间,应根据试验要求并作为成果绘制和提交下列资料:当进行稳定流抽水时,绘制涌水量、水位降-历时(、-)曲线、涌水量-水位降关系[=()]曲线(图1[地下水水位及流量历时曲线])及单位涌水量-水位降关系[=()]曲线。
当进行非稳定流抽水时,应绘制抽水井水位降与时间,观测孔水位降与抽水井距离()、水位恢复与时间的对数关系曲线,即-lg(图2[水位下降-时间对数关系曲线])、-lg、-lg(1+/)(图[kg2]3[水位恢复-时间对数关系曲线])曲线。
科技成果——现场快速测定水文地质参数的微水试验方法

科技成果——现场快速测定水文地质参数的微
水试验方法
技术简介
该成果通过瞬时抽注水、气压泵、振荡棒等激发手段,测量井孔内水位发生瞬时微量变化,根据测量到的水位随时间变化数据推导岩土体渗透性参数。
与传统的抽水试验、压水试验等方法相比,微水试验方法不需要用水、用电,也不从含水体中抽水或注水,具有设备简便、操作方法简单、试验周期短的特点,对地下水环境扰动极小,几乎不造成任何污染。
适用于水利水电工程环境、水文地质勘察、地下水资源评价中含水层水文地质参数获取。
技术特点
1、与传统的抽水试验、压水试验等方法相比,微水试验方法不需要用水、用电,也不从含水体中抽水或注水,具有设备简便、操作方法简单、试验周期短的特点对地下水环境扰动小;
2、强渗透地层只需几分钟就可完成一组微水试验,对于弱透水地层也能缩短一半时间以上,可用于各类勘察中水文地质参数的确定工作,丰富水文地质参数原位测试的技术方法
3、微水试验基本消除了人为误差以及越流的影响,精度较高,可提高工作效率,节约生产成本。
知识产权情况
实用新型专利3项,软件著作权1项
应用情况
从2017年至今,该成果在山东、河北、河南等省份水文地质勘察、地下水环境影响评价、生态治理与科研项目工作中应用,通过建立大型物理试验平台,分析和讨论不同激发方式与激发强度、不同含水体特征的微水试验技术手段,总结和凝练试验过程与参数求解方法,与其他原位试验成果进行比较,具有设备简便易携带、场地适应性好和数据精度可靠的特点,对周边环境扰动小。
水文地质参数-渗透系数和导水系数的确定

s = a1Q + a 2 Q 2 + KK + a n Q n
的
Q 以 1/ a1代换,分别进行计算。 H − h2
2
B B
3. 当 s / Q (或 Δh / Q )~Q关系曲线呈直线时,可采用作图截距法求出a1后,按上述
2
B B
方法计算。 二、单孔稳定流抽水试验观测孔水位下降资料求渗透系数 当利用观测孔中的水位下降资料计算渗透系数时,若观测孔中的值 s(或 Δh )在 s(或
2
(1)承压水完整孔: K = (2)承压水非完整孔:
Q R ln 2πsM r Q R M − l 1.12M (ln + ln ) 2πsM r πr l
(6-9)
当 M﹥150r, l/M﹥0.1 时, K =
(6-10)
当过滤器位于含水层的顶部或底部时, K =
Q R M −l M [ln + ln(1 + 0.2 )] (6-11) 2πsM r r l
(6-12)
(3)潜水完整孔: K = (4)潜水非完整孔:
Q R ln 2 2 π (H − h ) r
l/ h ﹥0.1 时,K = 当 h ﹥150r,
1.12h Q R h −l (ln + ⋅ ln ) 2 r l πr π (滤器位于含水层的顶部或底部时, K =
K=
t Q ln(1 + k ) 2 2 tT 2π ( H − h )
(6-24)
; tk—抽水孔从开始到停止的时间(min)
B B
; tT—抽水停止时算起的恢复时间(min)
B B
s—水位恢复时的剩余下降值(m) ; h—水位恢复时的潜水含水层的厚度(m) 。
第三章水文地质参数计算

内,s~t/r2曲线和W(u)~1/u标准曲线在形状上是相同的,只是
纵坐标平移了 Q 4T,横坐标平移了
* 。 4T
将二曲线重合,任选一匹配点,记下对应的坐标值,代入
Theis公式可求。
②具体步骤
在双对数坐标纸上绘制W(u)—1/u标准曲线。
在另一张模数相同的透明双对数纸上绘制实测的s—t/r2曲 线。
流计算。其降深s的计算公式为:
s
Q
2KM
r K0( B)
因为:r r B
B
对二式两边取对数,得:
lg
s
lg
k0
r B
lg
Q
2T
lg
r
lg
r B
lg
B
式中,lg Q 和 lg B是常数。曲线 lg s ~ lg r与曲线
2T
相似l,gk只0 是Br 坐~标lg Br平移了
lg
Q
2T
在另一张模数相同的透明双对数纸上绘制实测的 s—t曲线。
将实际曲线置于标准曲线上,在保持对应坐标轴彼 此平行的条件下相对平移,直至两曲线重合为止。
任取一配点(在曲线上或曲线外均可),记下匹配 点的对应坐标:W(u)、1/u、s、t,代入下式求参数:
T
Q
4 s
W
* 4T t
r
2
1 u
u
配线法的优点:
)
,
K M
T B2
B
§3-3 利用地下水动态资料确定 水文地质参数
利用地下水动态长期观测资料来确定水文地质参数 是一种比较经济的方法,并且确定参数的范围比前 者更为广泛,可以求出一些用抽水试验不能求得的 一些参数。
本节主要介绍给水度、降水入渗系数和潜水蒸发系 数的确定方法。
07第七章水文地质参数的计算

07第七章⽔⽂地质参数的计算第七章⽔⽂地质参数的计算⽔⽂地质参数是表征含⽔介质⽔⽂地质性能的数量指标,是地下⽔资源评价的重要基础资料,主要包括含⽔介质的渗透系数和导⽔系数、承压含⽔层的储⽔系数、潜⽔含⽔层的重⼒给⽔度、弱透⽔层的越流系数及⽔动⼒弥散系数等,还有表征与岩⼟性质、⽔⽂⽓象等因素的有关参数,如降⽔⼊渗系数、潜⽔蒸发强度、灌溉⼊渗补给系数等。
⽔⽂地质参数常通过野外试验、实验室测试及根据地下⽔动态观测资料采⽤有关理论公式计算求取,或采取数值法反演求参等。
第⼀节给⽔度⼀、影响给⽔度的主要因素给⽔度(µ)是表征潜⽔含⽔层给⽔能⼒或储⽔能⼒的⼀个指标,给⽔度和饱⽔带的岩性有关,随排⽔时间、潜⽔埋深、⽔位变化幅度及⽔质的变化⽽变化。
不同岩性给⽔度经验值见表7.l。
⼆、给⽔度的确定⽅法确定给⽔度的⽅法除⾮稳定流抽⽔试验法(参考《地下⽔动⼒学》等⽂献)外,还常⽤下列⽅法:1.根据抽⽔前后包⽓带上层天然温度的变化来确定p 值根据包⽓带中⾮饱和流的运移和分带规律知,抽⽔前包⽓带内⼟层的天然湿度分布应如图 7.1中的 Oacd 线所⽰。
抽⽔后,潜⽔⾯由 A 下降到 B (下降⽔头⾼度为功),故⽑细⽔带将下移,由aa '段下移到bb '段,此时的⼟层天然湿度分布线则变为图中的Oacd 。
对⽐抽⽔前后的两条湿度分布线可知,由于抽⽔使⽔位下降,⽔位变动带将给出⼀定量的⽔。
根据⽔均衡原理,抽⽔前后包⽓带内湿度之差,应等于潜⽔位下降Δh 时包⽓带(主要是⽑细⽔带)所给出之⽔量(µΔh )即h W W Z i i n i i=-∑=µ)(121故给⽔度为h W W Z i i n i i-=∑=)(121µ (7.1)式中:△Z i ——包⽓带天然湿度测定分段长度(m );△h ——抽⽔产⽣的潜⽔⾯下移深度(m );W 1i ,W 2i ;——抽⽔前后△Z i 段内的⼟层天然湿度(%);n ——取样数。
抽水试验确定水文地质参数的方法

在进行水 资源 分析工作时 ,常常使 用到地质水文参数 ,
目前 ,对 于 地 质 水 文 参 数 的推 求 方 法 有 很 多 种 ,其 中利 用抽
1 、观 6 ) ,其 中主 孔及 观 测 孔各 项 指 标 详 见 表 1 。
表 1 测 井 基 本 情 况 表
水实验来 确定地下水文参数是其 中之一 ,本文主要探 求通过
三 、各 项 参 数 计 算
域 内主要为粉砂一 粘土结构 ,该类结构岩性单一 ,层 次划 分较
为 明显 。 以主 孔 岩 性 可 知 ,0 - 5 . 7 m 为粉 砂 、 5 . 7  ̄ 7 . 6 m 为
依据实验区域的实际水文地质特点 , 其地下水为浅层地下 水 ,抽水底板是一层较厚的粘土层 ,因此是完全井观测。其中 对各项参数 的计算一般使 用以下 4中方法 ,即稳定流计算法 、 降深一 距离配线法、漏斗疏干法及直线解析法等方法 。
数据 的记 录 ,及做好安全保 障工作 ,抽 水结束要观测水位恢 复情况并做好记录工作 。
抽 水 过 程 中使 用 1 9 5柴 油 机 两 台 ,3台 自吸 泵 ,在 抽 水
参数 。该抽水试验实例首先进行非稳定流单 孔抽水 ,然后详细 的记录实验时各阶段 的信息 , 最后以这些资料信息来计算该 区 域的水文地质参数,本文详细的罗列了实验中应该收集的数据 及资料 ,并且详细的阐述 了如何计算水文地质参数 的方法 。
对 静 水 水 位 的观 测 ,准 备 需 要 使 用 的测 具 ,调 试 机 械 设 备 , 预 抽 及 对 观 测 人 员 的 培 训 等 ,抽 水 过 程 中 ,主 要 包 括 对 各 项
握某河流某段地质水文特点 ,因此选取了某河流流域某段抽水
水文地质参数计算与评价实验报告

水文地质参数计算与评价实验报告实验报告:水文地质参数的计算与评价一、引言水文地质参数是指描述水文地质条件的物理参数,对于水文地质调查和水文地质工程设计具有重要的意义。
本实验通过实地勘察和实验室测试的方法,对水文地质参数进行计算和评价。
本报告将详细介绍实验的目的、原理、方法、结果和讨论。
二、实验目的1.理解水文地质参数的概念和重要性;2.学会使用实地勘察和实验室测试的方法计算水文地质参数;3.掌握水文地质参数的评价方法。
三、实验原理1.水负荷试验:通过向井或孔隙中注入一定量的水,观察水位上升的情况,根据注入的水量和孔隙容积计算孔隙度和渗透系数。
2.介质颗粒分析:采用筛分和沉降法,将不同粒度的颗粒分离出来,计算颗粒组成和含水率。
3.渗透试验:在实验室中制备模型,通过施加一定的压力差和时间,测量渗透流量,并计算渗透系数。
四、实验方法1.实地勘察:选择一片地块,选取观测点,在井内注入一定量的水,观察水位上升的情况,并记录注水量和孔隙容积。
2.实验室测试:收集地块中的土样,进行介质颗粒分析和渗透试验,得到颗粒组成、含水率和渗透系数。
五、实验结果和讨论1.水负荷试验:根据实地勘察得到的数据,计算出孔隙度和渗透系数,用于评价地块的水文地质条件。
2.介质颗粒分析:通过实验室测试得到的颗粒组成和含水率,分析土壤的结构和水分状况,对水文地质条件进行评价。
3.渗透试验:根据实验室测试得到的渗透系数,评价土壤的渗透性能,为地下水运动和水文地质工程设计提供参考。
六、结论通过实地勘察和实验室测试的方法,成功计算和评价了水文地质参数。
根据计算和评价结果,可以得到地块的孔隙度、渗透系数、颗粒组成等参数,为水文地质调查和水文地质工程设计提供了重要的依据。
此外,本实验还掌握了水文地质参数的计算和评价方法,对于进一步研究水文地质领域具有一定的参考价值。
1.水文地质参数计算与评价实验指导书2.XXX等.水文地质学.北京:科学出版社,2024.。
通过抽水试验计算水文地质参数

一
、
工作 区概况
水 文地质钻孔 w1 为傍河钻孔 , 本次试验为潜水非完整井稳定 流
根据现场 收集 的原始数据 , 包括 流量 、 水位 、 降深等 , 绘 乌拉 陶勒盖矿 区地貌成 因类 型属剥蚀堆积类 型 ,矿 区植被较 发 单孔抽水试验 , 育, 属草原低丘陵地形地貌 , 地形起伏较小 , 切割深度 中等。矿区内地 制得相关抽水试验 Q、 s — t 曲线。 下水类型主要为基岩裂隙水 , 对矿床开采有一定影响的含水岩组。基 岩裂隙水又分为基岩裂隙风化层水和深部基岩裂 隙水 。
《 论文天地 l
的抽水试验 。同时说 明该 目的层 的补给条件较好 , 涌水 量与渗透系数 的数据关系符合客观实际 。
t( h )
3 0 4 0 5 0
0 2 0
孔 含 水 层 有 试段深 静 止 水位 试验时 试 号 效 厚 度 ㈣ 度 )
7 0 _ 3 8
的层 ) 的水文地质参数 。 三、 试验数据及处理计算
静水位 。 Q、 s - t 曲线形态符合地下水 涌水量与降深之间的关系 , 即涌水 量随着水位降深的增加而增大 , 其 增加量越来越小 。水位恢复 曲线呈 基本平滑 曲线 , 并需要很长时间 , 说 明该 区域含水层的补给条 件较 差 ,
渗透系数较小 。 该抽水试验求解水文地质参数为渗透系数 K, 计算公式采用 以下 公式进行 :
其中, K为渗透系数 ; Q为涌水量 ; s为抽水试验稳定流降深值 ; H 为含水层厚度 , 即降深值 s与含水层余 值 l 之和 ; b为抽水 钻: f L N河 流 的垂直距离 ; r 为钻孔半径 。 经计算 ,得到 本次 该工 作 区抽水 试验 的渗 透 系数 , K为 O . 0 0 3
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
3. 现场降水试验
现场成井施工,共完成了 2 眼抽水井和 6 眼水位观测井。现场试验分为 2 类水文地质试验,即 Slug Test(也称微水试验)与承压水抽水试验。各抽水井与观测井结构参数详见表 1。各抽水井与观 测井的平面剖面布置如图 2 所示
试验 类别
Slug Test
承压水 抽水试 验
孔号
K2 G5、G6 K1 G1、G2、G3 G4
孔径 (mm) 650
550 650 550 650
表 1 井结构参数表
井径 孔 深 滤管位置
mm m
m
273
34
30-33
219 34 273 24 219 21 273 21
30-33 17-23 17-20 17-20
井管长 m
30 30 17 17 17
填砾位 置m
28-34 28-34 14-24 14-21 14-21
作者简介:韩传梅(1982-),女(回族),新疆塔城人,硕士研究生。
-5-
本文以上海轨道交通七号线新村路站深基坑降水为例,详细论述了在第四纪巨厚沉积层地区含 水层结构,由于基坑维护连续墙未达到含水层底板,深基坑非完整井降水存在三维地下水流场。本 文在建立上海地铁七号线新村路站深基坑水文地质概念模型的基础上,为地下施工提供了科学决策 依据 [1] 。
2. 研究区概况
拟建上海轨道交通 7 号线一期工程,北起外环路站,南至东安路零陵路站,线路总长约 20 ㎞。 新村路车站位于新村路与岚皋路相交处,车站主体位于岚皋路。车站长度为 168.90 米,标准段宽度 19.70m,开挖深度约 15.30m,端头井开挖深度约 16.90m。车站主体采用二层二跨现浇钢筋混凝土结 构,采用地下连续墙作围护结构,设计初定标准段地墙深度为 29.0m,端头井地墙为 31.0m,明挖顺 作法施工,基坑等级为一级。
a (㎡/d)
0.301
9.04
8.53E-03
1.06E+03
0.607
18.22
2.77E-03
6.58E+03
0.435
3.260
3.80E-03
8.59E+02
B (m)
16 65 16.2
K ' (m/d)
0.23 0.016 0.019
已知越流因数 B ,可以转化为弱透水层的渗透系数 K ' 。由地质剖面图可知 G1,G2 井的 补给层为第④层和第⑧1 层。而 G3 井的补给层为第④层和第⑤3 层。由于第④层为淤泥质粘土
-3-
4.2 承压含水层抽水试验
本次抽水试验只针对第⑤2 层进行抽水试验,因基坑的南北地质情况不同,在两处布置两 组试验井,所以抽水试验共进行了两次。按下列公式计算参数:
⎧ ⎪T ⎪
=
Q
4 π [s ]W
(u )
⎪
4 ⋅ T ⋅ [t ] ⋅ u
⎨S = ⎪
Confirm the hydrogeological parameters whit test method
Han Chuanmei
State Key Lab. of Hydrology-Water Resources and Hydraulic Engineering,Hohai University, Nanjing (210098)
层位 第 ⑦1-1 层
第⑤2 层
Slug Test 试验在第⑦1-1 层中的 3 眼试验井(K2、G5、G6)中进行。承压含水层抽水试验 井针对第⑤2 层进行,共布置了 5 眼非稳定流抽水试验井(K1、G1、G2、G3、G4)。
5
6
7
8
9 10
18 19 20 21 22
图 2 抽水试验井群平面布置图
1. 引言
地下空间的开发已日益得到人们的重视,且不断向着大和深的方向发展。在地下空间建设过程 中会出现流沙以及承压水对基坑底部的突涌等一系列显现。所以基坑开挖的深度超过该区的地下水 位,必须采取降水措施。由于地质条件复杂,确定研究区的水文地质参数在深基坑降水方案中显得尤 为重要,因此,用不同方法确定水文地质参数是科学有效地对深基坑降水引起的地下水位变化进行 模拟预测的前提。
试验方法确定水文地质参数
韩传梅
河海大学水文水资源与水利工程科学国家重点实验室,南京(210098)
E-mail:hanmei424@
摘 要:针对长江三角洲第四纪松散沉积层厚度大、地下水位高、地下水量丰富,地下水流动状态 复杂等特点,本文以上海地铁七号线新村路站深基坑降水为例,利用深基坑降水试验推求水文地质 参数。 关键词:水文地质参数,深基坑降水,试验方法
r 2 ⋅1440
⎪ ⎪a ⎩
=
S
(2)
式中,T -导水系数 ( m2 / d );S -释水系数;s -水位降深( m );u -给水度; a -压力
传导系数( m2 / d )。采用曲线拟合方法得出计算结果(表 3)。
井号
G1 G2 G3
表 3 ⑤2 层水文地质参数汇总表
K(m/d) T (㎡/d) S (无量纲)
Fig 2 Pumping well design of groundwater withdrawal
-2-
4. 水文地质参数的推求
4.1 潜水含水层Slug Test试验
Slug Test参数计算公式如下:
⎧ ⎪K
=
(Tt
/
rc2 )
rc2 /(tl)
本工程场地地质情况相当复杂。主要有浅部土层中的潜水和深部粉性土层中的(微)承压水(图 1)。据区域资料,微承压水和承压水位,一般均低于潜水位,浅部土层中的潜水位埋深,一般离地 表面 0.3m~1.5m,年平均地下水位离地表面 0.5m~0.7m,低水位埋深为 1.5m;第⑤2 微承压水位埋 深在 3~6m,第⑦层承压水位埋深为 3~11m。潜水位和承压水位随季节、气候等因素有所变化。江 河边一定距离范围内,特别是有浅层粉性土分布区,其潜水位受潮汐影响较明显。受古河道切割影 响,4-22 轴分布有第⑤2 层粉性土,厚度为 6.0m~30.0m;在场区北部 1-4 轴,上覆⑤2 层缺失, 正常分布有第⑥层和第⑦层,第⑥层厚度约为 5.0m,第⑦1-1 层顶面埋深约 29.80m,层厚约 1.80m~ 3.50m。
Abstract The hydrogeological characteristics of the Yangtze delta region are large thickness in quaternary sediment layer, high groundwater table, abundance in groundwater resources, and complicated groundwater flow etc. The research is based on the Xincun Station along 7th line of Shanghai subway. A groundwater depression test at the deep pit was used to determine the hydrogeological parameters. Keywords: Hydrogeological parameter, deep foundation pit dewatering, test method
伏含水层组的水文地质参数建议值,如表 4 所示。
表 4 水文地质参数建议值
序号 地层 渗透系数 K ( m / d ) 贮水系数 S
1
⑤2 0.447
5.03E-03
2
⑦1-1 6.20E-03
1.06E-02
越流系数 B( m ) 32 /
所有试验数据均由水位传感器和数据采集仪自动采集并输入微机处理后绘制 lgs-lgt 或
⎨
⎪⎩S = a rc2 / rs2
(1)
式中, K -渗透系数( m / d ); rc -套管半径( m ); rs -滤管处钻孔半径( m );l -进水滤 管长度( m ); a -计算系数; t -瞬间取水开始后起算的时间( d ); S -贮水系数;T -导
水系数( m2 / d )。
采用 Cooper-Jacob 现行回归法 [2] ,在 K2、G5 和 G6 孔内分别进行。试验数据采用全自动 数据采集仪和振弦式孔隙水压力传感器采集,并与计算机连接进行数据处理。利用试验资料,
-1-
12 34 5
6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 图 1 地质剖面图
Fig1 Geological Section
依据本工程地质勘察资料,拟建场地地面标高为 3.28m~3.56m。⑤2 层承压含水层顶板距基坑 底面仅 1.0~2.0m。基坑开挖至设计标高后,坑内承压含水层顶板以上的覆土的自重压力远小于承压 水压力,必须采用深井减压降水措施。
(H / H0 ) − lg t 曲线,在微机上与标准曲线拟合后计算参数,最大限度地降低了人为的误差, 试验前传感器均进行了率定,所以试验得出的参数可以作为深基坑降水设计模型的依据。
5. 结论
在与地下水有关的工程中,确定地下含水层(尤其是承压含水层)水文地质参数中的渗 透系数K 对工程设计至关重要。通过实地抽水试验及稳定流计算方法检验,可以得到承压含 水层的渗透系数和蓄水系数等。
-4-
参考文献
[1] 骆祖江,武永霞. 第四纪松散沉积层深基坑降水三维非稳定流数值模拟.沈阳建筑大学学报(自然科学版) , 2006, 2: 181-185
[2] Cooper, H.H. and C.E. Jacob, “A generalized graphical method for evaluating formation constants and summarized well field history,” Am. Geophys. Union Trans., Vol.27,p.526-534(1946). [3] 薛禹群. 地下水动力学(第二版)[ M].北京:地质出版社,1997,115-116