Fracture network characterization using 1D and 2D data of the Moragy Granite body, southern Hungary

A disposal system for low- and medium-level nuclear waste in Hungary is being constructed inside the fractured rock body of the Lower Carboniferous Moragy Granite. Previous studies proved that the granitoid massif is rather heterogeneous in terms of lithological composition, brittle structure and hydrodynamic behaviour. A significant part of the body consists of monzogranite, while other portions are more mafic in composition and are monzonites. As a result of at least three significant brittle deformation events, the area is at present crosscut by wide shear zones that separate intensively fractured zones and poorly deformed domains among them. Due to late mineralization processes, some of these fractured zones are totally sealed and cannot conduct fluids, while others are excellent migration pathways. The spatial distribution of these two types nevertheless does not show any systematics. Hydrodynamic behaviour clearly reflects this heterogeneous picture; in some places, hydraulic jumps as great as 25 m at compartment boundaries can be detected. In this study, the fracture network of the Moragy Granite body is evaluated from a geometric aspect using datasets measured at a wide range of scales. 2D digitized images of a hand specimen, one large (20 x 60 m) and 12 smaller subvertical wall rocks (outcrops) and 120 images from tunnel faces representing the ground level of the underground repository site were analysed. Moreover, 1D data from 13 wells that all penetrate the granitoid massif were studied. Based on measured geometric data (spatial position, length, orientation, and aperture) fracture networks are simulated to study connectivity relations and for computing the fractured porosity and permeability at different scales. The results prove the scale-invariant geometry of the fracture system. Geostatistical calculations indicate that measurable fracture geometry parameters behave as regionalized variables and so can be extended spatially. Estimated localities of connected subsystems fit very well with fault zones mapped previously. Moreover, the spatial position of regimes of different hydrodynamic behaviours can be explained by connectivity relations both regionally and within wells.