地面沉降是全球许多地区面临的严峻环境问题，尤其是在地下水资源过度开采的区域。一维固结（1D Consolidation）作为土力学中的核心概念，直接关系到含水层系统的压缩与沉降预测。那么，究竟有哪些程序或模型可以模拟一维固结，从而助力地下水与包气带含水层的沉降估算？本文将为您梳理主流工具与应用场景。
一、一维固结与沉降估算的关系
一维固结是指饱和土体在侧向受限条件下，因孔隙水压力消散而缓慢压缩变形的过程。在地下水开采中，含水层及弱透水层的水头下降导致有效应力增加，引发土体压缩，最终表现为地表沉降。包气带（非饱和区）的土层同样可能因水含量变化而变形。准确模拟一维固结，是评估长期抽水对地面稳定性影响的前提。
二、主流模拟程序与模型
目前，可用于一维固结模拟的软件和模型种类繁多，从专业有限元分析到水文地质专用模块，各具特色。
1. MODFLOW 的 SUB 系列模块
美国地质调查局（USGS）开发的 MODFLOW 是地下水模拟的行业标准。其 SUB（Subsidence）包和 SUB-WT（Water Table）包专门用于模拟含水层系统的固结沉降。SUB 基于太沙基一维固结理论，能够模拟弹性与不可逆压缩，适用于承压含水层。SUB-WT 则拓展至潜水含水层及包气带，可处理非饱和区的有效应力变化。这些模块与 MODFLOW 的流场耦合，是区域沉降估算的首选工具。
2. PLAXIS
PLAXIS 是专业岩土工程有限元软件，支持弹塑性本构模型和流固耦合分析。其“固结分析”（Consolidation Analysis）功能可模拟一维与二维条件下的饱和土固结。通过输入土层分层、渗透系数和压缩指数，PLAXIS 能精确预测沉降随时间的变化。适用于局部精细研究，如地铁隧道、基坑工程附近的含水层沉降。
3. ABAQUS
ABAQUS 作为通用有限元软件，凭借丰富的用户子程序可自定义固结本构模型。利用其孔隙流体扩散/应力分析耦合功能，可模拟一维固结问题。对于包气带非饱和固结，ABAQUS 支持采用 van Genuchten 模型描述土水特征曲线，实现水分运动与变形的耦合。该软件学术应用广泛，但操作门槛较高。
4. FLAC / FLAC3D
FLAC 是二维/三维显式有限差分软件，擅长模拟岩土大变形和流固耦合。其内置的“一维固结”功能可快速计算单层或多层土在恒定荷载下的沉降。通过 FISH 语言可自定义固结参数与边界条件。FLAC 在矿山、油田采水沉降研究中应用较多。
5. 其他专用模型与代码
- IBS (Interbed Storage) 模型：针对薄层弱透水层固结编写的解析解程序，常用于 MODFLOW 预前估算。
- Settle3：Rocscience 开发的一维固结沉降计算软件，界面友好，支持多种固结理论。
- OpenGeoSys：开源多物理场耦合模拟平台，可结合 THM（热-水-力）框架模拟非饱和固结。
- Python 自定义脚本：基于 SciPy 和 FEniCS 等库，可自行编写一维固结有限元求解器，灵活性强。
三、应用场景与实际意义
在超采地下水的农业灌区（如华北平原、美国加州中央谷地），采用 MODFLOW 的 SUB 包模拟多年抽水导致的沉降，已取得较好效果。在包气带厚度大的地区，SUB-WT 模块则能考虑非饱和区的滞后效应。而对于深层含水层因油气开采引起的沉降，PLAXIS 或 ABAQUS 则被用于预测单井附近的地面形变。
这些模型的核心价值在于：将水文学与土力学结合起来，帮助决策者制定合理的抽水方案，避免诱发地质灾害。例如，通过调整开采层位或回灌策略，利用模型模拟不同情景下的沉降量，从而达到“保水护土”的双重目标。
四、挑战与展望
尽管现有程序众多，但一维固结模拟仍面临参数获取困难（如压缩系数、渗透率的空间变异性）、非饱和固结理论不完善等挑战。未来，随着遥感数据（InSAR 地表位移）与数值模型的同化技术发展，一维固结论将更精准。同时，开源框架的普及也将降低建模门槛。
对于研究者或工程师而言，选择工具的关键在于问题尺度：区域宏观分析首选 MODFLOW 系列，局部精细研究推荐 PLAXIS 或有限元软件，快速估算可借助专用小程序。唯有因地制宜、综合运用，方能使一维固结模拟真正服务于地下水资源可持续利用与地面沉降防治。