由于地下水不合理开采而导致诸如地面沉降等环境地质灾害,已经成为全球的普遍性问题;地面沉降不仅因为地面变形对沉降区域的建筑物产生严重影响,而且还因为地面失高导致沉降区防洪标准降低,在沿海地区还可能因为地面失高而导致大陆面积的减少;因此,关于地面沉降问题的研究,现已经受到科学界的高度关注。阜阳市地面沉降,是典型的因为开采含水层中地下水而导致相邻粘性土层产生固结的问题,在进行地面沉降模拟时,必须要对含水层中地下水流场变化进行模拟;而要模拟含水层中地下水流场的变化,又不可忽略相邻粘性土层固结而产生的“压密释水”。目前,关于地面沉降过程模拟与预测的有效途径之一是流固耦合理论,其中Biot流固耦合理论应用最为广泛;但由于Biot流固耦合模型涉及的参数过多,且有些参数难以确定,所以目前在实际中广泛应用还有一定困难;另外,当产生固结的粘性土层中的地下水流是非Darcy流时,Biot流固耦合理论还难以直接应用;除Biot理论外的另外一些流固耦合方法,大多都相对偏重固结的机理与变化过程,而将渗流模拟仅仅作为计算固结沉降的辅助过程,实际上,要实现这种条件下的地下水流系统模拟是比较困难的。本文在研究阜阳市地面沉降问题时,通过建立固结量与压密释水量的关系方程,并利用分步计算方法,将粘性土层的太沙基一维固结方程与含水层中的拟三维渗流方程实现了间接的耦合。将1997年1月至1999年12月的水位测验数据和地面沉降中心的实测沉降量用于模型的识别、检验,识别与检验过程中将1997年1月水位分布状态作为模型的初始流场,计算步长为4个月。利用识别检验后的耦合模型,预测了未来8年的沉降量。