《工程勘察》承压含水层非稳定流有限层分析-修改稿.docVIP

承压含水层非稳定流有限层分析 诸宏博1,2 王旭东1 宰金珉1 (1. 南京工业大学，江苏 南京 210009. 浙江省建筑科学设计研究院，浙江 杭州 310012) 摘 要：推导了以伽辽金法结合傅立叶函数为基础的各向异性承压含水层非稳定流有限层方程，编制了相应的计算程序。通过算例的有限层解与解析解的对比分析，在验证有限层法正确性的基础上，深入探讨了计算区域、计算项数等因素对计算结果的影响，得到了计算区域尺寸影响有限层解的收敛速度等结论，最后给出了一定精度条件下计算项数选用的近似估算公式。 关键词：有限层法；承压含水层；非稳定流；降深 0 引 言 Booker等运用了有限层法对层状弹性地基和土的固结问题做了较为系统的研究［2-4］。A.M. oskoorouchi和 B. Novrouzian将有限条方法推广应用于岩土工程分析中，对基坑开挖、隧道和土坝等问题进行了研究分析［5］。Stanley, Myron等将有限层理论应用于地下水分析［6］。张问清、赵锡宏、宰金珉等采用有限层方法探讨了任意力系作用下的层状弹性半空间问题［7］。宰金珉、王旭东等将有限层法应用于土与结构物的共同作用分析［8-9］。 地下水计算方法主要有解析法和数值法。解析法可以将描述地下水运动的各种物理量反映在一个表达式中，使用方便，但自然界地下水运动状态复杂，往往与解析解的假设条件相差甚远，其应用范围也受到一定的限制。对于一些复杂的地下水问题，就必须借助于数值模拟方法求其近似解。目前，地下水数值分析中主要有有限差分法、有限元法、边界元法等数值方法［10-12］。数值法能较为方便地处理解析解难以解决的复杂边界等问题，因此，得到广泛的应用。但是，数值法计算工作量大，尤其是在考虑三维模型的计算分析时，占用内存多。国内外学者不断尝试寻找有效的数值方法用于地下水模拟计算，薛禹群，叶淑君，谢春红等将多尺度有限元法应用于地下水模拟，节省了计算量提高了计算精度［13］。 本文吸取有限层法的优点，以伽辽金法结合傅立叶函数为基础，推导了承压含水层非稳定流有限层方程，编制了相应的计算程序，在验证方法的正确性基础上，探讨了计算区域、计算项数等因素对计算结果的影响。 1 承压含水层非稳定流模型的数学描述 (1) 式中，s(x,y,z,t)为降深，Kx、Ky、Kz分别为x，y，z方向的渗透系数，Ss为单位贮水系数，R(x,y,z,t)为源汇项，补给源取正值，排泄汇取负值。 (1)式可用微分算子表示为： (2) 式中，，L为微分算子。 图1 承压含水层物理模型与有限层层元 Fig.1 Configuration of confined aquifer physical model and finite layer models 2 有限层格式推导 有限层计算简图如图1所示。x和y方向的计算尺寸分别取为2a和2b，沿z轴方向将承压含水层离散成L个层元。 2.1 试探函数 设降深在x和y方向分别用解析函数表示，而z方向采用有限层元离散，层元内降深线性分布。则试探函数可表示为下列分离变量的函数项级数： (3) 式中，为降深试探函数；Φmnj (t)为待求系数；为满足定水头边界的已知函数，；M，N为级数项数；L为有限层层元数。 Ni(z)为线性插值基函数，其形式如下： (4) 2.2 伽辽金方程 根据加权余量法原理，可得有限层法伽辽金方程： (5) 式中,为权函数,。 将含水层离散为L个层元，则(5)式可改写为： (6) 将(3)式代入(6)式得： (7) 2.3 层元有限层方程 取一典型有限层层元l，其节面分别为l，l+1，层厚cl，如图2所示。 图2 典型有限层层元 Fig.2 Configuration of finite layer element models 由Amn(x, y)和Am′n′(x, y) 的正交性，可得： (8) (9) 故，当和时，由(7)式得： (10) 由基函数的性质可得： (11) 由(11)式可分别得到m、n阶的层元渗透矩阵系数、贮水矩阵系数和水量矩阵系数： (12a) (12b) (12c) 式中，i=l,l+1，j=l,l+1。 由(12a)式、(12b)式积分可得： (13a) (13b) 对于位于区域中心处的完整抽水井，由(12c)式积分可得： (13c) 式中，为抽水率。 层元l的有限层方程为： (14) 用矩阵可简化表示为： (15) 2.4 整体有限层方程 由层元有限层