《南京工业大学学报》--承压含水层地下水稳定流的有限层分析-诸宏博.docVIP

承压含水层地下水稳定流的有限层分析 诸宏博,王旭东,宰金珉 ( 南京工业大学土木工程学院，江苏 南京 210009 )中图分类号： 文献标识码：A （1981－）(College of Civil Engineering , Nanjing University of Technology, Nanjing 210009, China) Abstract: The finite layer method (FLM) is a quasi-numerical and quasi-analytic method, which is to discretize one dimension of the spatial domain using finite elements, approximating variations in the other two dimensions using continuous function. The method therefore reduces three-dimensional problems to one-dimensional problems effectively. The FLM can resolve the disadvantages radically, such as the large amount of works, the much memory and more time, et al, by using common numerical methods to simulate three-dimensional groundwater flows. Taking advantages of the FLM, finite layer formulations are derived for calculating the groundwater drawdown in the heterogeneous and anisotropic confined aquifer, which based on Galerkin method and Bessel function, and a computer program using FORTRAN language is developed. Two numerical examples are presented and compared with analytical solutions to demonstrate the validity of finite layer method. Key words: finite layer method; anisotropic; confined aquifer; steady groundwater flow; drawdown 地下水运动由于受含水层的非均质、各向异性和边界条件等因素的影响，其运动状态十分复杂，对于非均质、各向异性问题很难用解析方法进行求解，必须借助于数值模拟方法求其近似解［1-2］。国内外学者不断尝试寻找有效的数值方法用于地下水模拟计算，薛禹群等将多尺度有限元法应用于地下水模拟，节省了计算量，提高了计算精度［3］。朱学愚等采用边界元法对二维稳定流进行了分析计算［4］。但是，数值法在进行三维问题的计算分析时，常常会遇到计算分析占用内存多、计算工作量大等问题。 Cheung于1968年首先提出了有限条法，并成功地将有限条法应用于结构工程问题的分析［5］。该方法能有效地将三维问题简化为一维问题求解，从而大大减少单元数、节省计算工作量和内存的需要量。有限层法正是利用了有限条法的特点，进一步发展和完善起来的。它是一种对空间某一方向进行数值离散，而在其余两方向采用连续函数的半数值半解析方法。Small和Booker等运用了有限层法对土的固结与流变耦合问题做了较为系统的研究［6］；Stanley、Myron等将有限层理论应用于地下水分析［7］；Southcott也采用该方法对竖向荷载下的桩土相互作用问题进行了探讨［8］；支喜兰、王秉纲利用有限层方法对弹性地基上受对称荷载作用的圆形厚板进行了力学分析［9］；宰金珉、王旭东等将有限层法应用于土与结构物的共同作用分析［10-11］。有限层法已在岩土工程数值计算中得到了广泛的应用。 本文利用有限层法的优点，以伽辽金法结合贝塞耳函数为基础，推导了承压含水层稳定流有限层方程，编制了相应的计算程序，通过算例验证了该方法的正确性。 1 承压含水层稳定流模型的数学描述 假设层状非均质、各向异性、等厚承压含水层，初始水位水平，定水头边界条件，如图1所示。根据水量均衡原理和达西定律，柱坐标下承