利用地下水流模拟代码模拟地下气流的标准指南（D5719-13）

📋 概述与适用范围

本项标准由美国材料与试验协会土壤与岩石分委员会（D18.21）发布，编号为D5719-13，是对地下气流模拟方法的重要补充。其核心思路是利用地下水流模拟代码来近似计算地下空气流动，因为控制地下水流动的微分方程在数学形式上与描述气流运动的方程高度相似。标准最早源于地下环境修复工程中对简易气流模拟工具的需求，经过多次修订形成当前版本。适用对象包括不饱和带中的气相迁移问题，如土壤蒸气抽提系统设计、垃圾填埋气收集评估及地下空间通风模拟等。标准并不推荐特定代码，而是提供一套通用的变量替换与结果转换技术，用户可基于已有的地下水模拟软件（如基于有限差分或有限元方法的代码）直接开展气流分析。该指南需与地下水模型应用指南（D5447）等配套使用，同时强调了模型建立过程中的专业判断与实测数据验证的重要性。

在实际工程中，许多场地涉及气液两相流动，但专门的气流模拟软件往往价格高昂且掌握难度大。本指南通过数学变换降低了技术门槛，使工程人员能够借助成熟的地下水模拟工具完成初步评估。标准明确指出，替换方法适用于达西定律成立的层流范围，且假设气体为理想气体、温度恒定。对于非达西流或高度非等温情形，则需要更复杂的模型。指南结构清晰，从符号定义、参数转换到结果输出均给出说明，为后续的模拟验证与应用提供了基础框架。

⚙️ 模拟原理与方法

地下水流动的基本方程（如布西涅斯克方程）基于达西定律和质量守恒，而空气在孔隙介质中的流动同样遵循达西形式，差异仅在于流体物性。本指南的核心原理是将空气流动方程改写为与地下水流动方程一致的形式，通过定义“等效水头”将气压变量替换为水头变量。具体步骤包括：首先选择合适的地下水模拟代码，建立空间离散网格并确定边界条件；然后将空气渗透率转换为等效水力传导系数，转换中需考虑空气密度与黏度随压力的变化；接着在代码中将气压边界条件转化为等效水头边界，并运行模拟；最后将计算得到的水头场反向推导为气压场或流速场。

关键近似在于假设空气为可压缩理想气体，且流动过程等温。标准中给出了明确的变量对应关系：例如，将气流方程中的气压平方项与地下水方程中的水头项对应，并引入中间变量以消除压缩性影响。对于各向异性介质，渗透率张量的处理同样适用。在实际操作中，用户需根据目标问题的压力水平预先评估压缩效应：当压力变化超过百分之十时，必须采用平方压力法；低压差时可直接使用线性近似。指南还强调网格独立性检验与稳态/瞬态选择的差异，瞬态模拟需合理估计储气系数，这一系数可通过孔隙度与流体压缩系数计算得到。

📊 技术参数与指标

标准中列出了一组关键符号及其单位，下表汇总了这些基本参数，便于用户在实际建模时统一量纲。

下表给出了指南中建议的地下水变量与气流等效变量之间的对应关系，这是实施模拟时进行参数替换的具体依据。

这些转换关系是保证模拟结果一致性的基础，用户在使用不同代码时应严格遵循上述量纲替换原则，避免因单位误算导致结果偏离。

🔬 工程应用与注意事项

本指南最典型的应用场景是土壤气相抽提（SVE）系统的设计优化。在污染场地修复中，通过模拟地下气流压力分布与流速场，可以合理确定抽气井的间距、深度及抽气速率。此外，垃圾填埋场的气体收集系统设计、地下停车场通风评估以及氡气迁移风险评估均可借助该指南进行初步分析。实际工程中，场地表征数据（如空气渗透率、孔隙度、饱和度）的准确性直接影响模拟质量。建议在建模前开展气压或示踪试验以获取现场参数，并用实测压力数据对模型进行校正。

使用中的常见问题包括：边界条件设定不当（如将地表设为定压边界而忽略大气压变化）、忽视近井区域的非达西流动、以及低估了含水饱和度对气体相对渗透率的影响。指南也提醒用户注意，地下水模拟代码通常不包含气体压缩项，因此对高真空或高压力工况需格外谨慎。当模拟范围包含部分饱和带时，气-水两相流动的耦合效应可能显著，此时单一气相假设的误差会增大，可参考更复杂的变饱和地下水流代码。质量控制方面，建议进行网格灵敏度分析、稳态与瞬态结果对比，并与解析解或现场实测数据比对。

另外，标准强调该指南不取代专业工程师的判断。模拟结果必须结合地质统计资料和工程经验进行解释，避免单纯依赖数值输出。对于涉及人身安全或重大环境影响的项目，应使用经充分验证的气流专用软件或现场试验数据进行复核。

❓ 常见问题解答

📥 标准文件下载