第
34
卷第
4
期
2006
年
4
月
西北农林科技大学学报(自然科学版)
Vo
l.
34
No.4
April
2006
1our.
of
Northwest Sci-Tech
Univ.
of
Agn.
and
Fo
r.
(Na
t.
Sci.
Ed.)
黄土潜水层深井降水三维有限元模拟研究
康婷白应生
1
,郭增玉
2
o
空军工程大学工程学院,陕西西安
710038;2
西安理工大学水利水建学院,陕西西安
710048)
[摘
要]
视抽水井为定流量边界,并假定流量沿井壁的分布是均匀的,在此条件下建立了黄土厚潜水层深
井抽水时的三维渗流模型,并用
c++
语言编制了三维有限元程序,计算了单井抽水、双井抽水时水位随时间的变化
规律,计算结果与观测结果误差较小。另计算了少陵塘隧洞第一期降水方案的降水效果,表明其不能满足施工进度
的要求,并对此方案作了修改,在后期施工中取得了良好的降水效果。
[关键词]
黄土
i
潜水层;三维有限元;深井降水
[中图分类号]
TU46+3
[文献标识码]
A
随着高层建筑、地铁、隧道等工程的发展,降低
地下水位的问题日益突出。如何有放评价降水方案
的合理性及预测动水位,较好地指导施工开挖进度,
是工程实践中普遍存在的一大难题。目前,阵水问题
大多采用经验公式或按二维流问题来计算,但存在
与实际情况不相符合的问题。三维渗流模型在地下
水资源的优化管理
[1]
中比较多见,但并不能直接应
用于阵水工程,因为二者所要达到的目标不同,所考
虑的渗流场范围也不同。对地下水资源进行合理开
采时,影响范围内地下水位的下降受到严格控制,阵
深一般较小,抽水井常当作定水头边界。在降水工程
中,受施工场地、施工工期的限制,在短期内须将水
位降到预定高度,降深一般较大,井壁水位随抽水时
间的变化而变化。
用有限元方法模拟降水问题时,抽水井的处理
是一个难点
[2]
。在以往的二维有限元分析中,一般不
将排水井作为流量边界,而是将其作为网格的一个
结点来近似地反映排水作用。文献
[3J
指出,排水孔
直径大小对渗流场有很大影响,把排水井当作一个
结点不能完全反映排水井的尺寸效应,多数情况下
将导致较大的计算误差。且根据本人的计算经验,这
样做极易引起水位异常。
本文结合黑河续建工程少睦塘隧洞降水工程,
将抽水井当作定流量边界处理,井假定流量沿井壁
的分布是均匀的,参考文献
[4J
的方法建立了黄土区
三维潜水渗流数学模型,并参考文献
[5-7J
用
c++
语
言编制了三维有限元程序,分别计算了单井抽水和
[收稿日期]
2005-07-01
[文章编号]
1671-9387(2006)04-0139-04
双井抽水时地下水位随抽水时间的变化规律。同时
对少睦塘隧洞第一期降水方案作了修改,并在后期
施工中取得了良好的降水效果。说明本文所建立的
三维潜水渗流模型和所采用的计算方法具有一定的
可靠性,可以为黄土地区类似工程提供参考。
1
少陵塘区潜水渗流数学模型
1.
1
工程概况
黑河少陵塘隧洞位于陕西省西安市长安县少睦
塘区,根据水文地质与工程地质勘察报告,地下水位
埋深为地面以下
25
m
左右,而隧洞位于塘面以下
35~52
m
,即地下水位最大高出洞底以上
18
mo
为
了保证施工安全和施工质量,隧洞开挖前必须进行
人工降水。
1.
2
概念模型及渗流数学模型
根据少睦塘隧洞的水文地质和工程地质勘察报
告,少陵晾区黄土层地下水为孔隙潜水,地层岩性主
要为黄土、黄土状壤土夹古土壤,呈成层分布山。群
井抽水时,降水区域内的地下水流具有明显的三维
流特征,因此将其概化为三维潜水流动,地层概化为
均质各向异性,渗流基本偏微分方程为:
主
(hzit)
十三
(hvi
业)十主
(hzi
业)
= 0
dX dX
dy"dy
dZ dZ
(X
,
y
,
Z)
ε
U
(1)
式中,工
,
y
,
z
为迪卡尔坐标
;k
x
,ky
,
k
z
为
x
,
y
,
z
轴方
向上的渗透系数
(m/h)
;ψ
为水头
;U
为所选定的渗
流区域。
[作者简介]
康 婷
0977-)
,女,陕西佳县人,讲师,硕士,主要从事岩土工程研究。