Model assessment of reductive dechlorination as a remediation technology for contaminant sources in fractured clay: Modeling tool

Appendix J Sensitivity analysis

J.1 Sensitivity to fracture aperture and spacing

The model response to change in fracture aperture and spacing is assessed for the three different degradation scenarios by analyzing the output parameters shown on Figure J.1: time to remove 90% of the initial total contaminant mass, average time for the fracture outlet concentration to be lower than a limit concentration and finally average peak concentrations for the two daughter products.

The different geometric configurations of Table G.1 are used with the transport and degradation parameters from the base case scenario.

Figure J.1 - Parameters for sensitivity analysis of fracture aperture and spacing – total mass on right axis (blue line)

Figure J.1 - Parameters for sensitivity analysis of fracture aperture and spacing – total mass on right axis (blue line)

Clim is defined for each of the chlorinated compounds as a function of the groundwater quality standard (Cstandard) recommended in Miljøstyrelsen [2005] and an assumed dilution factor of 10 between the concentration at the fracture outlet and the resulting concentration in the underlying aquifer. Hence Clim = Cstandard*10. The dilution factor of 10 is assumed here, but does not change the main conclusions from this part.

Table J.1 - Groundwater quality standard and output parameters from model

  Cstandard Clim
  µg/L µg/L µmol/L
TCE 1 10 0.08
cis-DCE 1 10 0.10
VC 0.2 2 0.03

Click here to see Figure J.2

Figure J.2 – Average time to reach Ci < Clim for no degradation (a), degradation in fracture (b), degradation in fracture and reaction zone in matrix (c) and degradation in fracture and the whole matrix (d), note the log vertical scale

Click here to see Figure J.3

Figure J.3 – Average time to remove 90% of the initial contaminant mass, for no degradation (a), degradation in fracture (b), degradation in fracture and reaction zone in matrix (c) and degradation in fracture and the whole matrix (d), note the log vertical scale

Click here to see Figure J.4

Figure J.4 – Average maximum concentration for the daughter products (DCE and VC), for degradation in fracture (a), degradation in fracture and reaction zone in matrix (b) and degradation in fracture and the whole matrix (c), note the log vertical and horizontal scale

The single fracture/matrix model is not sensitive to the fracture aperture, except for the maximum daughter products concentration in case of degradation in fracture (Figure J.4 (a)) and a limited sensitivity for the average time to reach an output concentration below Clim in case of degradation in the fracture (Figure J.2 (b)). In these two cases, aperture reduction results in a higher water velocity in the fracture and therefore the daughter products do not have the time to be produced. Otherwise the non-sensitivity of this parameter is explained by the definition of flow in the fracture, which depends on the fracture spacing only (see Section 5.1.2 in the main report). On the contrary the model results are very sensitive to the fracture spacing: the mass removal time increases with fracture spacing, as well as the time to reach an output concentration below Clim. Furthermore, assuming degradation in the whole matrix leads to a decrease of the clean-up times. In this case, the model is less sensitive to fracture spacing.

J.2 Global sensitivity analysis

The sensitivity analysis is performed on all independent parameters (Table J.2), of which the value is changed by +/- 20%. The resulting change in the three output parameters (time to remove 90% of the initial contaminant mass, average time to reach Ci < Clim and average peak concentration of daughter products) is normalized to calculate the sensitivity index (see Appendix D). The sensitivity analysis is performed on the third degradation scenario (degradation in the fracture and in a reaction zone in the matrix) with a homogenous TCE concentration as initial condition. In order to be able to compare the different simulations, this initial concentration is corrected in order to maintain the same initial total mass.

Table J.2 - Independent parameters for sensitivity analysis (in orange transport parameters, in green degradation parameters)

Parameter Symbol Value – base caase Unit
Net recharge I 0.1 m/year
Fracture spacing 2B 0.3 m
Fracture aperture 2b 7*10-4 m
Sorption coefficient TCE Kd_TCE 1 L/kg
Sorption coefficient DCE Kd_DCE 0.7 L/kg
Sorption coefficient VC Kd_VC 0.3 L/kg
Matrix porosity φ 0.33 -
Exponent p p 1 -
Longitudinal dispersivity in fracture aL 0.1 m
Max growth rate TCE µTCE 730 year-1
Max growth rate DCE µDCE 138.7 year-1
Max growth rate VC µVC 51.1 year-1
Specific yield Y 5.2*108 cell.µmol-1
Initial biomass XO 108 cell.L-1
Half velocity coefficient DCE KDCE 9.9 µmol.L-1

The four most sensitive independent parameters are the same for the different output considered, matrix porosity, fracture spacing, net recharge and TCE sorption coefficient. The most sensitive parameters are then the ones controlling transport, especially diffusion/sorption processes, and not dechlorination. The limiting process in this system is the counter diffusion out of the matrix, which is controlled by sorption, flushing of the fracture and fracture spacing.

The two least sensitive parameters are the fracture aperture and longitudinal dispersivity in the fracture. The low sensitivity of this last parameter is due to the fact that transport in the fracture is mainly advective and not dispersive.

Table J.3 - Sensitivity index for the three output parameters

Parameter M< 10%Mini Parameter Ci < Clim Parameter max Ci
Matrix porosity 55.0 Matrix porosity 115.0 Fracture spacing 31.0
Net recharge 52.5 Fracture spacing 95.0 Matrix porosity 24.8
Fracture spacing 42.5 Sorption coefficient TCE 61.7 Net recharge 21.9
Sorption coefficient TCE 35.0 Net recharge 58.3 Sorption coefficient TCE 18.7
Sorption coefficient DCE 20.0 Specific yield 47.5 Max growth rate TCE 17.6
Specific yield 20.0 Initial biomass 47.5 Initial biomass 16.3
Initial biomass 20.0 Sorption coefficient DCE 34.2 Specific yield 16.1
Exponent p 10.0 Max growth rate DCE 22.5 Sorption coefficient DCE 14.8
Max growth rate DCE 10.0 Exponent p 20.0 Exponent p 6.3
Half velocity coefficient DCE 7.5 Max growth rate TCE 16.7 Max growth rate VC 5.2
Sorption coefficient VC 5.0 Half velocity coefficient DCE 16.7 Max growth rate DCE 3.9
Max growth rate TCE 5.0 Max growth rate VC 7.5 Half velocity coefficient DCE 2.1
Max growth rate VC 5.0 Sorption coefficient VC 3.3 Sorption coefficient VC 2.0
Fracture aperture 0.0 Fracture aperture 1.7 Longitudinal dispersivity in fracture 1.6
Longitudinal dispersivity in fracture 0.0 Longitudinal dispersivity in fracture 0.8 Fracture aperture 0.3

Sensitivity analysis has also been performed by varying the parameters in the typical ranges found in the literature and the same conclusions can be done, concerning the most and least sensitive parameters.

 



Version 1.0 July 2009, © Danish Environmental Protection Agency