The contamination of groundwater aquifers by chlorinated solvents is a water-resource issue of great importance worldwide. Such contaminants are difficult to treat because they are often released as dense non-aqueous phase liquids (DNAPLs). Research shows that the application of remediation technologies to both the NAPL source and dissolved plume can lead to more efficient remediation, rather than to either alone. In these remediation efforts, analytical models that evaluate behavior and fate of contaminants do provide a better understanding of the performance of these remedial technologies. To the best of our knowledge, there exist no analytical model of simulating the plume migration of multiple contaminants with capabilities of accounting for both NAPL source and plume remediation simultaneously and different retardation for original chlorinated solvent contaminant and its degradation byproducts. In this study, we present a new analytical model for remediating both NAPL source and downgradient contaminant plume in groundwater at sites contaminated with chlorinated solvents and their degradation products with different retardation factors as well as considering both NAPL source and plume remediation simultaneously. A source model that accounts for the depletion of mass by the processes of dissolution or first-order decay reactions, corresponding with the removal or destruction of the source mass, is coupled to a plume reactive transport model. The source model is accounted for by relating source mass to the flux-averaged source discharge concentration through a power function. The developed analytical model considers 1-D advection, 3-D dispersion, first-order decay reactions and ingrowth as well as linear isothermal equilibrium sorption. The proposed analytical solution was derived through successive application of the Laplace transform in time and the double finite Fourier cosine transform regarding y and z. The correctness of the analytical model and its auxiliary FORTRAN computer program code are proved by showing excellent agreements between the simulated plume concentrations of all contaminants obtained from the derived analytical model and from a semi-analytical model available in the literature. Application of the proposed analytical solutions illustrates that the use of identical retardation factors for all contaminants may lead to underestimation or overestimation of the mobility of the contaminants, in cases when the retardation factors of the individual contaminants are greatly different from the identical retardation factor value adopted in all contaminants. From the experiments on six scenarios corresponding six remedial treatments, we found out that both the enhanced source decay and partial removal of source mass are main controlling factors at reducing the concentrations of all the contaminants, whereas plume decay leads to effective reduction in the concentrations of PCE, however, rather it causes unfavorable increases of the concentrations of the degradation byproducts.

Ultimately, the developed model is used to better understand the impacts of various possible combinations of remedial efforts and management decisions on remediation of the subsurface contamination and quantify the benefit of a certain remediation decision.

氯化溶剂对地下水含水层的污染是世界范围内非常重要的水资源问题。此类污染物难以处理，因为它们通常以致密非水相液体 (DNAPL) 的形式释放。研究表明，将修复技术应用于 NAPL 源和溶解羽流可以导致更有效的修复，而不是单独使用。在这些修复工作中，评估污染物行为和归宿的分析模型确实提供了对这些修复技术性能的更好理解。据我们所知，不存在模拟多种污染物羽流迁移的分析模型，能够同时考虑 NAPL 源和羽流修复以及原始氯化溶剂污染物及其降解副产物的不同延迟。在这项研究中，我们提出了一种新的分析模型，用于修复受氯化溶剂及其降解产物污染的场地地下水中的 NAPL 源和降级污染物羽流，并同时考虑 NAPL 源和羽流修复。通过溶解或一级衰变反应过程解释质量消耗的源模型，对应于源质量的去除或破坏，与羽流反应传输模型耦合。通过幂函数将源质量与通量平均源放电浓度相关联来说明源模型。开发的分析模型考虑了 1-D 平流、3-D 色散、一级衰减反应和向内生长以及线性等温平衡吸附。所提出的解析解是通过在时间上连续应用拉普拉斯变换和双有限傅里叶余弦变换得到的y和z. 分析模型及其辅助 FORTRAN 计算机程序代码的正确性通过显示从派生分析模型和文献中可用的半分析模型获得的所有污染物的模拟羽流浓度之间的良好一致性来证明。所提出的分析解决方案的应用表明，当单个污染物的延迟因子与采用的相同延迟因子值有很大差异时，对所有污染物使用相同的延迟因子可能会导致低估或高估污染物的迁移率在所有污染物中。从对应六种补救措施的六种情景的实验中，

最终，开发的模型用于更好地了解各种可能的补救措施和管理决策组合对地下污染修复的影响，并量化某个修复决策的收益。