A new method for generating a pipe network to handle channelled flow in fractured rocks

The pipe network generation proposed here assumes that flow in a 3D fracture field is highly channelled, and that there exist invariant crux nodes corresponding to crossroads connecting preferential flow paths whatever the mean head gradient's direction. The method allows preferentially for these crux nodes by conditioning the simulated network first on a ‘seeding' of nodes and second, in generating 1D pipes between these nodes. Note that usually, nodes are the consequence (the intersection) of the pipe generation. Pipes can be also conditioned on their length and other features, but this conditioning is not strict, since all pipes must pass by the nodes. In the end, the network is simple without dead ends; each cluster is flowing and connected to the boundaries of the system. This type of network is expected to facilitate further calculations, for instance solving inverse problems in deforming the network by simply moving the crux node locations.