A Stochastic Study of Flow Anisotropy and Channelling in Open Rough Fractures

Journal Article · 16 July 2019 · Rock Mechanics and Rock Engineering

DOI: https://doi.org/10.1007/s00603-019-01907-4 · OSTI ID:1580874

Marchand, Sophie  ^[1]; Mersch, Olivier ^[1]; Selzer, Michael  ^[1]; Nitschke, Fabian ^[1]; Schoenball, Martin  ^[2]; Schmittbuhl, Jean  ^[3]; Nestler, Britta  ^[1]; Kohl, Thomas  ^[1]

---

1. Karlsruhe Inst. of Technology (Germany)
2. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
3. School of Observatory of Earth Sciences, Strasbourg (France)

The quantification of fluid flow in rough fractures is of high interest for reservoir engineering, especially for deep geothermal applications. Herein, rough self-affine fractures are stochastically generated with incremental shear displacement and geometrically described by two aperture definitions, the vertical aperture a_vert and the effective aperture a_eff. In order to compare their effect on fracture flow, such as anisotropy and channelling, Local Cubic Law (LCL) model-based 2D fluid flow is simulated. The particularity of this approach is the combination of a stochastic generation of self-affine fractures with a statistical analysis (560 individual realizations) of the impact of the LCL’s aperture constraint on fracture flow. The results show that aperture definition affects the quantitative interpretation of flow anisotropy and channeling as well as the aperture distribution of the fractures with shearing. Higher values of mean aperture for a given fracture are found using a_vert, whereas the aperture standard deviation is larger with a_eff. In addition, flow anisotropy is significantly sensitive to aperture definition for small shear displacements and shows a relative higher dispersion with a_eff. Thus, LCL prediction models based on a_vert are expected to lead to higher dispersion of anisotropy results with a higher uncertainty (factor ~ 2). Realizations based on avert lead to an enhanced clustering of high flow rates for higher shearing displacements. This channeling development results in higher total flow rates for these simulations. Finally, these findings support the direct calibration of pre-existing LCL anisotropy simulations based on avert towards more representative results using a_eff.

• Research Organization: Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
• Sponsoring Organization: USDOE Office of Science (SC)
• Grant/Contract Number: AC02-05CH11231
• OSTI ID: 1580874
• Journal Information: Rock Mechanics and Rock Engineering, Journal Name: Rock Mechanics and Rock Engineering Vol. 53; ISSN 0723-2632
• Publisher: SpringerCopyright Statement
• Country of Publication: United States
• Language: English