The development of the adjoint state method for the system of partial differential equations (PDEs) describing radionuclide transport in heterogeneous discrete fractured porous media in steady-state regime is presented in this study for both continuous and discrete adjoint approaches. The continuous adjoint approach has its own advantage in the sense that the adjoint system can be derived in terms of PDEs which can be solved analytically for specific performance measures. These PDEs can then be used to obtain a physical interpretation of the adjoint state functions. The discrete adjoint approach permits a simpler mathematical treatment based on the use of the same matrix (or its transpose) of the linear system of the discretized transport equations referring to from hereon as the forward problem. Analytical solutions for both the forward problem and its associated adjoint problem have been derived for radionuclide transport in layered formations including a fracture using an analytical inversion of tridiagonal matrices. A set of 14 formulas of sensitivity coefficients with respect to pertinent system parameters has been derived and presented in integral forms. Analytical results have been compared with numerical results and have shown excellent agreements. The discrete adjoint approach has been coupled with the differential sensitivity analysis method and applied to a real field case involving the performance assessment of a site currently being considered for underground nuclear storage in northern Switzerland, Zürich Nordost. The results show that when a fracture is present in the vicinity of waste repository, the state parameter controlling flow through the fracture, referred to herein as fracture transmissivity, is among the important parameters which should be considered in sensitivity analysis studies.
Read full abstract