3-D Fracture Network Modelling in Hydropower Engineering Based on Optimal Monte Carlo Simulation

Abstract

Effectively and accurately modelling the spatial relation of fracture surfaces is crucial in the design and construction of large hydropower dams having a complex underlying geology. However, fracture surfaces are randomly formed and vary greatly with respect to their spatial distribution, which makes the construction of accurate 3-D models challenging. In this study, we use an optimal Monte Carlo simulation and dynamic conditioning to construct a fracture network model. We found the optimal Monte Carlo simulation to effectively reduce the error associated with the Monte Carlo method and use dynamic conditioning to ensure the consistency of the model with the actual distribution of fractures on the excavation faces and outcrops. We applied this novel approach to a hydropower station on the Jinshajiang River, China. The simulation results matched the real sampled values well, confirming that the model is capable of effectively and accurately simulating the spatial relations in a fracture network.

Keywords

Hydropower engineering / 3-D fracture model / Dynamic conditioning / Monte Carlo simulation

References

[1] Jafari A, Babadagli T. Estimation of equivalent fracture network permeability using fractal and statistical network properties. J Petrol Sci Eng, 2012, 92: 110-123.

[2] Khamforoush M, Shams K. Percolation thresholds of a group of anisotropic three-dimensional fracture networks. Phys A, 2007, 385: 407-420.

[3] Li LC, Tang CA, Li G, et al. Numerical simulation of 3D hydraulic fracturing based on an improved flow-stress-damage model and a parallel FEM technique. Rock Mech Rock Eng, 2012, 45(5): 801-818.

[4] Jia HB, Tang HM, Liu YR, et al. Theory and engineering application of 3-D network modeling of discontinues in rock mass, 2008, Beijing: Science Press (in Chinese)

[5] Starzec P, Andersson J. Application of two-level factorial design to sensitivity analysis of keyblock statistics from fracture geometry. Int J Rock Mech Min Sci, 2002, 39(2): 243-255.

[6] Liesenfeld R, Richard JF. Monte Carlo methods and Bayesian computation: importance sampling. Int Encyclopedia Soc Behav Sci, 2001, 8(12): 10000-10004.

[7] Wathugala DN, Kulatilake PHSW. A general procedure to correct sampling bias on joint orientation using a vector approach. Comput Geotech, 1990, 10(1): 1-31.

[8] Kulatilake P, Chen J, Teng J, et al. Discontinuity geometry characterization in a tunnel close to the proposed permanent shiplock area of the Three Gorges Dam site in China. Int J Rock Mech Min Sci Geomech Abstr, 1996, 33(3): 255-277.

[9] Mauldon M. Estimating mean fracture trace length and density from observations in convex windows[J]. Rock Mech Rock Eng, 1998, 31(4): 201-216.

[10] Mauldon M, Dunne WM, Rohrbaugh MB. Circular scanlines and circular windows: new tools for characterizing the geometry of fracture traces. J Struct Geol, 2001, 23(2): 247-258.

[11] Kulatilake PHSW, Um JG, Wang M, et al. Stochastic fracture geometry modeling in 3-D including validations for a part of Arrowhead East Tunnel, California, USA. Eng Geol, 2003, 70(1): 131-155.

[12] Jimenez-Rodriguez R, Sitar N. Inference of discontinuity trace length distributions using statistical graphical models. Int J Rock Mech Min Sci, 2006, 43(6): 877-893.

[13] Yang ZJ, Su XT, Chen JF, et al. Monte Carlo simulation of complex cohesive fracture in random heterogeneous quasi-brittle materials[J]. Int J Solids Struct, 2009, 46(17): 3222-3234.

[14] Most T. Stochastic crack growth simulation in reinforced concrete structures by means of coupled finite element and meshless methods, 2005, Weimar: Bauhaus University.

[15] Vorechovsky M. Interplay of size effects in concrete specimens under tension studied via computational stochastic fracture mechanics. Int J Solids Struct, 2007, 44(9): 2715-2731.

[16] Yang Z, Xu XF. A heterogeneous cohesive model for quasi-brittle materials considering spatially varying random fracture properties. Comput Methods Appl Mech Eng, 2008, 197(45): 4027-4039.

[17] Papadrakakis M, Papadopoulos V, Lagaros ND, et al. Vulnerability analysis of large concrete dams using the continuum strong discontinuity approach and neural networks. Struct Saf, 2008, 30(3): 217-235.

[18] Jing L, Stephansson O. The basics of fracture system characterization—field mapping and stochastic simulations. Dev Geotech Eng, 2007, 85(7): 147-177.

[19] Tivadar MT. Determination of geometric parameters of fracture networks using 1D data[J]. J Struct Geol, 2010, 32(7): 878-885.

[20] Pahl PJ. Estimating the mean length of discontinuity traces. Int J Rock Mech Min Sci Geomech Abstr, 1981, 18(3): 221-228.

[21] Kulatilake PHSW, Wu TH. Estimation of mean trace length of discontinuities. Rock Mech Rock Eng, 1984, 17(4): 215-232.