随着国民经济的迅速发展，我国对水资源的需求量越来越大。由于地下水具有分布广、水质好且不容易被污染等优点，许多城市把地下水作为主要的供水水源，大量开采来满足城市生活和工农业用水的需要。但是过量开采地下水引发的地面沉降问题日益严重。因此，开采地下水的同时应对地面沉降进行合理预测，以便达到科学防控地面沉降的目的。 本文以比奥固结理论为基础，将地下水渗流场与土体应力场相耦合，建立了地下水渗流与地面沉降三维耦合模型，应用伽辽金有限单元法推导了其空间八节点六面体等参元的有限元离散方程。模型考虑了流量边界、井处理、土体非线性特性、土的渗透性随应力状态的变化。在有限元模型的基础上，本文采用Fortran语言开发了地下水渗流与地面沉降三维耦合程序。 最后本文将编制的程序用于预测江阴市浅层地下水开采出现的地面沉降，结果表明： (1)本文建立的地下水开采与地面沉降三维耦合模型，可一次性的求解出地下水位与土体变形位移量，能应用于实际工程。 (2)模型采取土体的弹塑性本构模型，考虑土体的非线性特征及渗透性与应力之间的相互影响，使得模型对土体变形过程的刻画更切合实际。 (3)优化网格节点编号时采用了AD算法，减小了程序运算时占用的带宽，加快了模型的运算速度。