Natural source zone depletion (NSZD) of light non-aqueous phase liquids (LNAPLs) may be a valid long-term management option at petroleum impacted sites. However, its future long-term reliability needs to be established. NSZD includes partitioning, biotic and abiotic degradation of LNAPL components plus multiphase fluid dynamics in the subsurface. Over time, LNAPL components are depleted and those partitioning to various phases change, as do those available for biodegradation. To accommodate these processes and predict trends and NSZD over decades to centuries, for the first time, we incorporated a multi-phase multi-component multi-microbe non-isothermal approach to representatively simulate NSZD at field scale. To validate the approach we successfully mimic data from the LNAPL release at the Bemidji site. We simulate the entire depth of saturated and unsaturated zones over the 27 years of post-release measurements. The study progresses the idea of creating a generic digital twin of NSZD processes and future trends. Outcomes show the feasibility and affordability of such detailed computational approaches to improve decision-making for site management and restoration strategies. The study provided a basis to progress a computational digital twin for complex subsurface systems.

轻质非水相液体 (LNAPL) 的自然源区耗竭 (NSZD) 可能是石油受影响地点的有效长期管理选项。然而，其未来的长期可靠性需要建立。NSZD 包括 LNAPL 组分的分配、生物和非生物降解以及地下的多相流体动力学。随着时间的推移，LNAPL 成分会耗尽，分配到不同相的成分会发生变化，可用于生物降解的成分也会发生变化。为了适应这些过程并预测数十年到数百年的趋势和 NSZD，我们首次采用了多相多组分多微生物非等温方法来代表性地模拟现场规模的 NSZD。为了验证该方法，我们成功地模拟了来自 Bemidji 站点的 LNAPL 版本的数据。我们模拟了 27 年的释放后测量中饱和和非饱和带的整个深度。该研究推进了创建 NSZD 流程和未来趋势的通用数字孪生的想法。结果表明，这种详细的计算方法在改善场地管理和恢复策略的决策方面具有可行性和可负担性。该研究为推进复杂地下系统的计算数字孪生提供了基础。