Physical Address
304 North Cardinal St.
Dorchester Center, MA 02124
Physical Address
304 North Cardinal St.
Dorchester Center, MA 02124
标准 D6033‑16 由美国材料与试验协会(ASTM)于 2016 年正式发布,其前身可追溯至 1996 年的最初版本,本次修订主要完善了功能组分类体系并补充了新兴过程模块。该标准的适用对象是所有用于模拟地下水流、溶质运移、生物及化学过程以及多相流等问题的计算代码,无论其维数、饱和状态或耦合程度如何。作为地下水建模系列指南之一,它与 D5447(地下水模型应用指南)、D5490(模型校准指南)、D5609、D5610、D5611 以及 D5718 等标准形成完整的技术支撑框架。这些标准并不孤立,而是共同构成从模型概念设计、代码选择、参数校核到结果评价的全流程规范。本标准的核心价值在于:它不替代专业判断,而是提供一套统一的术语和分类方法,使工程人员能够用同一种语言描述不同代码的能力边界,从而显著降低因功能误解导致选型错误的风险。
地下水建模代码本质上是非参数化的数学描述,通过一组控制方程表达地下系统中水流、化学物质及生物活动的因果关联。本标准提出的分类协议首先将代码功能拆解为若干主功能组,再在每个组内细化出具体的子功能和特征选项。例如“水流模拟”主功能组下包含饱和流、非饱和流、密度依赖流及多相流等子功能,每一子功能又可进一步描述其时间依赖(稳态/瞬态)、空间离散方式(有限差分/有限元)、求解算法(直接法/迭代法)等属性。描述过程需遵循系统化步骤:第一步,根据工程问题的物理过程列出所需功能组;第二步,对照标准提供的表格逐一比对各代码在这些组中的实现方式;第三步,记录代码特有或附加的功能,如参数自动率定、敏感性分析等。该方法论的核心优势在于将模糊的“代码能力”转化为可分解、可对比的技术条目,使不同背景的专家能够就功能是否满足需求达成一致。标准还强调,对于研究发展出的新型代码,若标准尚未涵盖其功能,仍应在文档中单独编写功能描述章节,以保持完整性。
标准通过大量表格对各功能组进行了层次化分解,并辅以典型代码的实现示例。下表根据标准附录整理出最核心的主功能组及其子功能对照关系,供实际工程描述时直接套用。
| 🟦 主功能组 | 📏 子功能类别 | 📐 功能描述要点 |
|---|---|---|
| 水流模拟 | 饱和水流、非饱和水流 | 基于达西定律或理查兹方程,描述地下水在多孔介质中的流动 |
| 溶质与生物运移 | 保守溶质、反应性溶质、微生物 | 模拟对流、弥散、吸附、降解及微生物生长衰减过程 |
| 多相流 | 非水相液体(NAPL)、气‑水‑油三相 | 考虑不同流体间的相界面交换及残余饱和度 |
| 热量传输 | 热传导、热对流、热弥散 | 耦合温度场与流场,模拟地热系统或热污染迁移 |
| 耦合过程 | 密度依赖、化学‑力学耦合、生物‑化学耦合 | 处理流体密度随浓度变化、土体变形等交互作用 |
| 🎯 功能特征 | ⚡ 选项/参数 | 📌 典型范围与说明 |
|---|---|---|
| 维数 | 一维、二维、三维 | 一维用于柱试验,三维用于复杂真实场地 |
| 时间模式 | 稳态、瞬态 | 瞬态需定义时间步长与存储格式 |
| 饱和状态 | 饱和、非饱和、变饱和 | 变饱和需提供水分特征曲线参数 |
| 密度依赖 | 是、否 | 若为是,则还需指定浓度‑密度关系式 |
| 耦合算法 | 全耦合、顺序耦合、迭代耦合 | 全耦合精度高但计算量大,顺序耦合更灵活 |
在实际工程中,本标准主要用于以下场景:在项目初期编写模型需求说明书时,利用标准分类框架逐项列出模拟必须包含的功能;在招标或采购阶段,比对各供应商所提供代码的功能描述是否完全覆盖需求;在成果评审时,检查模拟代码是否确实实现了其声称的功能。典型操作流程为——首先根据水文地质条件确定主要过程(如饱和流+反应性溶质运移),然后从标准表格中抽取对应子功能,并附加特定要求(例如必须支持三维变饱和流和全耦合传输),由此形成一份“功能需求清单”。再将候选代码的官方文档按照同一标准整理出“功能供给清单”,两者对比即可快速发现不匹配项。注意事项包括:不要仅关注主功能组而忽略子功能中的实现限制(如某些代码虽宣称能模拟非饱和流,但仅适用于特定边界条件);对于多过程耦合问题,应特别注意耦合方式的描述(顺序耦合可能导致质量守恒误差);标准框架无法涵盖代码的数值精度和稳定性,这部分需另通过验证标准(如 D6012)评估。质量控制方面,建议由两位独立专家分别完成功能描述,然后进行交叉核对,以减少主观遗漏。