ANALYSIS OF RADIONUCLIDE TRANSPORT THROUGH FRACTURE NETWORKS BY PERCOLATION THEORY

Presented are results of numerical simulations for radionuclide diffusion through fracture networks in geologic layers. Actual fracture networks are expressed as two-dimensional honeycomb percolation lattices. Random-walk simulations of diffusion on percolation lattices are made by the exact-enumeration method, and compared with those from Fickian diffusion with constant and decreasing diffusion coefficients. Mean-square displacement of a random-walker on percolation lattices increases more slowly with time than that for Fickian diffusion with the constant diffusion coefficient. Though the same relation of mean-square displacement vs. time as for the percolation lattices can be obtained for a continuum with decreasing diffusion coefficients, spatial distribution of probability densities of finding the random-walker on the percolation lattice differs from that on a continuum with the decreasing diffusion coefficient. The percolation model results in slow spreading near the origin and fast spreading in the outer region, whereas the decreasing-diffusion coefficient model shows the reverse because of smaller diffusion coefficient in the outer region. We could derive a general formula that can include both Fickian and anomalos diffusion in terms of fractal and fracton dimensionalities and the anomalous diffusion exponent.