Natural fractured media are highly unpredictable because of existing complex structures at the fracture and at the network levels. Fractures are by themselves heterogeneous objects of broadly distributed sizes, shapes, and orientations that are interconnected in large correlated networks. With little field data and evidence, numerical modeling can provide important information on the underground hydraulic phenomena. However, it must overcome several barriers. First, the complex network structure produces a structure difficult to mesh. Second, the absence of a priori homogenization scale, along with the double fracture and network heterogeneity levels, requires the calculation of large but finely resolved fracture networks resulting in very large simulation domains. To tackle these two related issues, we reduce the highly complex geometry of the fractures by applying a local transformation that suppresses the cumbersome meshing configurations while keeping the networks fundamental, geological, and geometrical characteristics. We show that the flow properties are marginally affected while the problem complexity (i.e., memory capacity and resolution time) can be divided by orders of magnitude. The goal of this article is to propose a method of resolution which takes into account the geometrical complexity met in the networks and which makes it possible to treat a few thousand fractures. The principal aim of this article is to present a tool to slowly modify the structures of the fracture networks to have a good quality mesh with a marginal loss in precision.
Introduction.The underground waste repository projects and exploitation of hot, dry, rock geothermal energy have spurred studies in fracture networks transport properties. Many site-specific studies have flourished around various projects [7], which are generally based on a careful characterization of the structure of the fractured rock mass, with the classical difficulty of deducing three-dimensional information from one-or two-dimensional field data [3]. Then the hydraulic properties are computed by using reconstructed model networks, based on the experimental geometrical characteristics, and various flow models. The flow fluid modeling in underground media requires taking into account the very high geological heterogeneity. The complexity of the geological mediums comes from metamorphic and sedimentary processes, and mechanics which create structures having hydraulic properties varying on several orders of magnitude and correlated on a broad range of scales. A general overview of stochastic generation of fractured media problems is given in [11]. Numerically, the challenge is to integrate this range of heterogeneity in models of great extension and good resolution outcome, with linear systems from 10 6 to 10 8 unknown factors. In the fractured medium with a matrix of very low permeability, the fluid flow is focused in the highly heterogeneous fractures [21]. The objective of numerical modeling consists of simulating hydraulic phenomena in a large number of ne...