地下水开采与地面沉降控制三维全耦合模型研究

地下水开采与地面沉降控制三维全耦合模型研究

骆祖江1,曾　峰1,李　颖2

1.河海大学水文地质研究所,南京　210098

2.长春工程学院勘查与测绘学院,长春　130021

摘要:针对第四纪松散沉积层中地下水开采所引起的地面沉降问题,以比奥固结理论为基础,考虑到土体的非线性特征及土体的渗透性随应力状态的动态变化,引入邓肯2张非线性模型和渗透率动态模型,将地下水渗流场和土体应力场耦合起来,并以吴江市浅层地下水开采为例,建立了浅层地下水开采与地面沉降三维全耦合数值模型。在对模型进行校正、识别的基础上,以微承压含水层未来10a 内地下水位埋深不低于含水层顶板、地面沉降量不超过50mm 为约束条件,预测了吴江市各镇的地下水可采资源量合计为1

203.59×104m 3/a 。

关键词:比奥固结理论;地面沉降;三维全耦合模型;非线性;地下水;开采

中图分类号:P641 文献标识码:A 文章编号:167125888(2009)0621080207

收稿日期:2009202219

基金项目:江苏省与国土资源部合作重大项目(200312300009)

作者简介:骆祖江(1964—

),男,江苏吴江人,教授,博士生导师,主要从事水文地质与工程地质研究,E 2mail :luozujiang @ 。

Study on Three 2dimensional Full Coupling Model of G roundw ater

Exploitation and Land 2Subsidence Control

L UO Zu 2jiang 1,ZEN G Feng 1,L I Y ing 2

1.Research I nstit ute of H y d rogeolog y ,Hohai Universit y ,N anj i ng 210098,Chi na

2.S chool of Pros pecting &S urveying Engineering ,Changchun I nstit ute of Technolog y ,Changchun 130021,China

Abstract :Aimed at t he land 2subsidence cased by groundwater exploitation in Quaternary loose se 2diment s ,based on Biot ’s consolidation t heory ,took t he shallow groundwater exploitation in Wujiang as t he case st udy ,a t hree dimensional f ull coupling numerical model between t he hydraulic action of t he seepage field and t he stress field was established ,in which Duncan 2Chang no n 2linear model and t he dy 2namic models of permeability were int roduced wit h considering t he non 2linear characteristics and t he dy 2namic variatio n of hydraulic conductivity of soils wit h st ress.Based on identification and verification of t he model ,wit h t he const raint condition t hat t he groundwater level would not be lower t han half of t he roof dept h of t he first confined aquifer and t he land 2subsidence wo und not exceed 50mm ,t he allowable exploitable quantity of groundwater in every town was p redicted ,t he exploitable quantity of groundwa 2ter is 1203.59×104m 3/a in t he whole area.

K ey w ords :Biot consolidation t heory ;land 2subsidence ;t hree dimensional f ull coupling model ;non 2linear ;groundwater ;exploitation

引言

地面沉降是一种由多种因素引起的地面高程缓慢降低的地质现象,其中过量开采地下水是造成地面沉降的最主要原因。江苏省吴江市地处长江三角

洲,经济发达,城市供水对地下水需求量大且开采历

　第39卷　第6期2009年11月 吉林大学学报(地球科学版)Journal of Jilin University (Earth Science Edition ) Vol.39　No.6Nov.2009　

史较早,已形成了区域性的地下水位降落漏斗,进而导致了区域性的地面沉降地质灾害问题;因此建立区域性的地下水开采与地面沉降模拟模型,准确预测基于地面沉降控制的地下水可采资源量已成为当务之急。目前,国内外关于地下水渗流与地面沉降问题的计算模型主要有水土分算模型[1]、部分耦合模型[2]和全耦合模型[3]

。但大多数模型并未考虑到土体的变形参数以及渗透性随着土体中应力场改变的动态变化。笔者根据比奥固结理论,并考虑土体的非线性特征及土的渗透性随应力状态的动态变化,结合吴江市地下水系统的水文地质机制,建立了吴江市浅层地下水开采与地面沉降的三维全耦合模型;同时计算出地下水开采过程中地下水位和地面沉降的动态变化,为地面沉降的科学防控提供决策依据。1

土水耦合模型

1.1

比奥(Biot)理论

假定饱和土体中,土体骨架是线弹性、变形微

小、地下水渗流符合达西定律且不可压缩,取z 坐标

向上为正,应力以压为正,则三维比奥固结方程如

下[4]:

-G 2

w x G 1-2v ・55x (5w x 5x +5w y 5y +5w z 5z )+5u 5x

=0,

-G 2w y G 1-2v ・55y (5w x 5x +5w y 5y +5w z 5z )+5u

5y =0,

-G 2w z

G 1-2v ・55z (5w x 5x +5w y 5y +5w z 5z )+5u

5z

=-γ。

(1)

5t (5w x x +5w y

y +w z z )1γw [

5x (k x 5u x )+55y (k y 5u 5y )+55z (k z (5u 5z

+γw ))

]=0。(2)式中:G 为剪切模量(Pa );v 为泊松比;w x 、w y 、w z 分别为x ,y ,z 方向上的位移分量(m );u 为孔隙水压力(Pa );k x ,k y ,k z 分别为x 、y 、z 方向上的渗透系数(m/d );γ为土的重度(N/m 3);γw 为水重度(N/

m 3)。

利用Galerkin 加权余量法离散方程,考虑到土体的非线性特性,取Δt 时间内的位移增量来代替位移,将式(1)、

(2)离散成增量形式[526]:K C C T

Δt H +

B Δr Δu =Δf Δt Q -Δt H u 。(3)式中:K,H 分别是固体刚度矩阵和流体“刚度”矩

阵;C 是耦合矩阵;Q 是流量矩阵;B 是关于自由面的积分矩阵;Δr 为结点位移增量(m );Δf 为等效节点力增量(Pa )。上述方程结合一定的定解条件即可运用伽辽金有限单元法进行求解。1.2方程的定解条件1.2.1初始条件地应力初始条件:采用土体的自重应力估算土

体的初始应力,σz 0=γ

z ,σx 0=σy 0=K 0γz 。(4)

式中:σx 0,σy 0为土体的初始水平向应力(Pa );σz 0为土体的初始垂向应力;z 为计算点深度(m );K 0为静止侧压力系数,K 0=1-sin φ′砂土

0.95-sin φ′粘土

,φ′为有效内摩擦角。位移初始条件:w (x ,y ,z ,t )|t =0=0。(5)

孔隙水压力初始条件:

u (x ,y ,z ,t )|

t =0

=u 0(x ,y ,z )。(6)

式中:u 0(x ,y ,z )为区域内已知初始孔隙水压力。

1.2.2边界条件

孔隙水压力边界条件Γ1:

u (x ,y ,z )|Γ1=u s 。

(7)式中:u s 为水头边界Γ1上的已知孔隙水压力。

流量边界条件Γ2:K

5H

5n |Γ2

=q L 。(8)

式中:q L 为边界Γ2上的已知单位面积流量。

自由面边界条件Γ4:

u =0;q =-μ

5u 5t

cos θ

。(9)

式中:μ为土体给水度;θ为自由面外法线方向与垂线的交角;q 为通过自由面边界Γ4的单位面积流量。1.2.3

自由面边界条件的处理

在地下水开采的计算过程中,自由面的位置是未知的,需要通过迭代才能求出[728]。应用改进的截止负压法[9210]将式(2)拓展到整个区域,使待定的自

由面边界固化,迭代求解。并引入罚函数M (P ),得到式(2)的弱表达式:

∫

Ω

( sk (p γ)+M (p ) sk z +

s 5εv 5t )d Ω-∫

Γ4

sq d Γ=0。(10)式中:

1

801　第6期 骆祖江,曾　峰,李　颖:地下水开采与地面沉降控制三维全耦合模型研究