Assessment of Method to Model Slip of Isolated, Non-Planar Fractures Using 3DEC

"The Swedish Nuclear Fuel and Waste Management Company are planning to construct a permanent storage facility for high-level nuclear waste at about 500 m depth in fractured crystalline rock. One potential concern for the long-term performance of the repository relates to large fracture shear displacements, since these may damage intersected waste containers. Although details in fracture surface geometries are known to influence slip magnitudes and distributions, fractures are commonly approximated as planar features with lab-scale properties due to lack of relevant field-scale data and to limitations in available modelling tools. The purpose of the work described in this paper is to explore the possibilities to use the distinct element code 3DEC to generate and analyze models with isolated, non-planar fractures and, based on the results of such models, get perspectives on the influence of basic disturbances such as ridges, steps and waves. Given the model set-up, fracture segments with orientations deviating from the critical one will contribute to stabilize the fracture. As a result, all analyzed non-planar fractures slip less than the critically oriented planar fracture. We demonstrate that the angle of deviation is the most influential parameter determining the slip magnitude for comparable surface geometries. Similarly inclined steps and ridges have almost equal slip reductions despite differences in surface area restricting the movement. Although it may be possible to translate the influence of such irregularities to an effective coefficient of friction that could be used as input to models with idealized planar fractures, additional analyses with randomly shaped and distributed asperities, different fracture sizes and with alternative stress loads, need to be carried out to devise a defensible strategy for this type of parameter variation. INTRODUCTIONThe Swedish Nuclear Fuel and Waste Management Company’s (SKB) concept for deep geological storage of high-level nuclear waste is the KBS-3 method according to which the waste will be encapsulated in copper canisters and deposited, surrounded by compacted bentonite clay, in vertical deposition holes at about 500 m depth in fractured crystalline rock (SKB 2011). A site has been selected in the Forsmark area, which is located some 120 km north of the Swedish capital Stockholm. One potential concern for the long-term performance of the repository that has been identified by SKB (SKB 2011) is related to large shear displacements along repository host rock fractures induced by a nearby earthquake, since these may cause mechanical damage to intersected waste containers. Due to lack of relevant data at the appropriate field scale and to limitations in available modelling tools, fractures are commonly approximated as planar features with lab-scale properties. The standard numerical tool, 3DEC (Itasca 2013), used in SKB’s rock mechanics analyses, for example, generates planar fractures by default. Results from dynamic 3DEC earthquake-models (Fälth et al. 2010) currently provide input to layout decisions, i.e., the diameter of the smallest planar fractures that must be avoided when deciding on waste container positions at different distances from potential earthquake faults. It is, however, well known that fracture surfaces are rough or undulating on a wide range of scales (e.g., Candela et al. 2012) and that such deviations from planarity influence the mechanical behaviour of the fractures (e.g., Dieterich and Smith 2009, Marshall and Morris 2012, Ritz and Pollard 2012). It will probably never be a realistic option to attempt to represent the geometrical complexity of a population of fractures of different size, shape and surface structure in any large-scale rock mechanical analysis relevant for the performance assessments of the nuclear waste repository. It may however be possible to find approximate ways to represent different styles and degrees of non-planarity by assigning mod"