Search

Grouting plays a crucial role in dam safety. Due to the concealment of grouting activities, complexity of fracture distribution in rock masses and rheological properties of cement grout, it is difficult to analyze the effects of grouting. In this paper, a computational fluid dynamics (CFD) simulation approach of dam foundation grouting based on a 3D fracture network model is proposed. In this approach, the 3D fracture network model, which is based on an improved bootstrap sampling method and established by VisualGeo software, can provide a reliable and accurate geometric model for CFD simulation of dam foundation grouting. Based on the model, a CFD simulation is performed, in which the Papanastasiou regularized model is used to express the grout rheological properties, and the volume of fluid technique is utilized to capture the grout fronts. Two sets of tests are performed to verify the effectiveness of the Papanastasiou regularized model. When applying the CFD simulation approach for dam foundation grouting, three technical issues can be solved: (1) collapsing potential of the fracture samples, (2) inconsistencies in the geometric model in actual fractures under complex geological conditions, and (3) inappropriate method of characterizing the rheological properties of cement grout. The applicability of the proposed approach is demonstrated by an illustrative case study—a hydropower station dam foundation in southwestern China. [ABSTRACT FROM AUTHOR]

Fracture network modelling plays an important role in many application areas in which the behaviour of a rock mass is of interest. These areas include mining, civil, petroleum, water and environmental engineering and geothermal systems modelling. The aim is to model the fractured rock to assess fluid flow or the stability of rock blocks. One important step in fracture network modelling is to estimate the number of fractures and the properties of individual fractures such as their size and orientation. Due to the lack of data and the complexity of the problem, there are significant uncertainties associated with fracture network modelling in practice. Our primary interest is the modelling of fracture networks in geothermal systems and, in this paper, we propose a general stochastic approach to fracture network modelling for this application. We focus on using the seismic point cloud detected during the fracture stimulation of a hot dry rock reservoir to create an enhanced geothermal system; these seismic points are the conditioning data in the modelling process. The seismic points can be used to estimate the geographical extent of the reservoir, the amount of fracturing and the detailed geometries of fractures within the reservoir. The objective is to determine a fracture model from the conditioning data by minimizing the sum of the distances of the points from the fitted fracture model. Fractures are represented as line segments connecting two points in two-dimensional applications or as ellipses in three-dimensional (3D) cases. The novelty of our model is twofold: (1) it comprises a comprehensive fracture modification scheme based on simulated annealing and (2) it introduces new spatial approaches, a goodness-of-fit measure for the fitted fracture model, a measure for fracture similarity and a clustering technique for proposing a locally optimal solution for fracture parameters. We use a simulated dataset to demonstrate the application of the proposed approach followed by a real 3D case study of the Habanero reservoir in the Cooper Basin, Australia. [ABSTRACT FROM AUTHOR]

Different failure modes during fracture shearing have been introduced including normal dilation or sliding, asperity cut-off and degradation. Attempts have been made to study these mechanisms using analytical, experimental and numerical methods. However, the majority of the existing models simplify the problem, which leads to unrealistic results. With this in mind, the aim of this paper is to simulate the mechanical behaviour of synthetic and rock fracture profiles during direct shear tests by using the two-dimensional particle flow computer code PFC2D. Correlations between the simulated peak shear strength and the fracture roughness parameter D recently proposed by Rasouli and Harrison () are developed. Shear test simulations are carried out with PFC2D and the effects of the geometrical features as well as the model micro-properties on the fracture shear behaviour are studied. The shear strength and asperity degradation processes of synthetic profiles including triangular, sinusoidal and randomly generated profiles are analysed. Different failure modes including asperity sliding, cut-off, and asperity degradation are explicitly observed and compared with the available models. The D parameter is applied to the analysis of synthetic and rock fracture profiles. Accordingly, correlations are developed between D and the peak shear strength obtained from simulations and by using analytical solutions. The results are shown to be in good agreement with the basic understanding of rock fracture shear behaviour and asperity contact degradation. [ABSTRACT FROM AUTHOR]

Rock is a heterogeneous geological material. When rock is subjected to internal hydraulic pressure and external mechanical loading, the fluid flow properties will be altered by closing, opening, or other interaction of pre-existing weaknesses or by induced new fractures. Meanwhile, the pore pressure can influence the fracture behavior on both a local and global scale. A finite element model that can consider the coupled effects of seepage, damage and stress field in heterogeneous rock is described. First, two series of numerical tests in relatively homogeneous and heterogeneous rocks were performed to investigate the influence of pore pressure magnitude and gradient on initiation and propagation of tensile fractures. Second, to examine the initiation of hydraulic fractures and their subsequent propagation, a series of numerical simulations of the behavior of two injection holes inside a saturated rock mass are carried out. The rock is subjected to different initial in situ stress ratios and to an internal injection (pore) pressure at the two injection holes. Numerically, simulated results indicate that tensile fracture is strongly influenced by both pore pressure magnitude and pore pressure gradient. In addition, the heterogeneity of rock, the initial in situ stress ratio ( K), the distance between two injection holes, and the difference of the pore pressure in the two injection holes all play important roles in the initiation and propagation of hydraulic fractures. At relatively close spacing and when the two principal stresses are of similar magnitude, the proximity of adjacent injection holes can cause fracturing to occur in a direction perpendicular to the maximum principal stress. [ABSTRACT FROM AUTHOR]

Fractured rock has often been conceptualized as a dual-continuum system for many practical applications. This study proposes a systematic approach to deal with multiphase flow in a dual-continuum system. Considering that fluid flow occurs in pore volumes (including fracture apertures), we first develop a so-called pore-space conservation equation for deformed fractured rock and then combine this equation with fluid mass balance to derive governing equations for multiphase flow associated with rock deformation. Constitutive relationships are also presented for describing stress dependence of hydraulic properties and effective mechanical parameters for bulk rock body (as a function of the corresponding parameters for fracture and matrix continua). Finally, we applied the developed approach to a CO geological sequestration problem to demonstrate the usefulness of the approach. [ABSTRACT FROM AUTHOR]

The failure mechanism of hydraulic fractures in heterogeneous geological materials is an important topic in mining and petroleum engineering. A three-dimensional (3D) finite element model that considers the coupled effects of seepage, damage, and the stress field is introduced. This model is based on a previously developed two-dimensional (2D) version of the model (RFPA2D-Rock Failure Process Analysis). The RFPA3D-Parallel model is developed using a parallel finite element method with a message-passing interface library. The constitutive law of this model considers strength and stiffness degradation, stress-dependent permeability for the pre-peak stage, and deformation-dependent permeability for the post-peak stage. Using this model, 3D modelling of progressive failure and associated fluid flow in rock are conducted and used to investigate the hydro-mechanical response of rock samples at laboratory scale. The responses investigated are the axial stress-axial strain together with permeability evolution and fracture patterns at various stages of loading. Then, the hydraulic fracturing process inside a rock specimen is numerically simulated. Three coupled processes are considered: (1) mechanical deformation of the solid medium induced by the fluid pressure acting on the fracture surfaces and the rock skeleton, (2) fluid flow within the fracture, and (3) propagation of the fracture. The numerically simulated results show that the fractures from a vertical wellbore propagate in the maximum principal stress direction without branching, turning, and twisting in the case of a large difference in the magnitude of the far-field stresses. Otherwise, the fracture initiates in a non-preferred direction and plane then turns and twists during propagation to become aligned with the preferred direction and plane. This pattern of fracturing is common when the rock formation contains multiple layers with different material properties. In addition, local heterogeneity of the rock matrix and macro-scale stress fluctuations due to the variability of material properties can cause the branching, turning, and twisting of fractures. [ABSTRACT FROM AUTHOR]