Fast method for simulation of radionuclide chain migration in dual porosity fracture rocks

In fractured rocks with a porous rock matrix such as granites, radionuclides will flow with the water in the fracture network. The nuclides will diffuse in and out the rock matrix where they can sorb and be considerably retarded compared to the water velocity. A water parcel entering the network will mix and split at the fracture intersections and parts of the original parcel will traverse a multitude of different fractures. The flowrates, velocities, sizes and apertures of the fractures can vary widely. Normally one must solve the transport equations for every fracture and use the effluent concentration as inlet condition to the next fracture and so on. It is shown that under some weakly simplified conditions it suffices to determine one single parameter group containing information on the flow wetted surface that a water parcel contacts along the entire path. It is also shown how this can be obtained. Then, solving the transport equations only once for time and location along the path gives the concentration and nuclide flux of every nuclide in the chain everywhere along a path. The same solution actually is valid for every path in the network. This dramatically reduces the computation effort. The same approach can be used for models based on streamtubes.