Search Results (127)

Basalt is an important geotechnical material for engineering construction in Southwest China. However, it has complicated structural features due to its special origin, particularly the widespread occurrence of hidden cracks. Such discontinuities significantly affect the mechanical properties and engineering stability of basalt, and related research is lacking and unsystematic. In this work, taking the underground caverns in the Baihetan Hydropower Station as the engineering background, the size-dependent mechanical behaviors and excavation responses of basalt with hidden cracks were systematically explored based on a synthetic rock mass (SRM) model combining the finite-discrete element method (FDEM) and discrete fracture network (DFN) method. The results showed that: (1) The DFN–FDEM model generated based on the statistical characteristics of the geometric parameters of hidden cracks can consider the real structural characteristics of basalt, whereby the mechanical behaviors found in laboratory tests and at the engineering site could be exactly reproduced. (2) The representative elementary volume (REV) size of basalt blocks containing hidden cracks was 0.5 m, and the mechanical properties obtained at this size were considered equivalent continuum properties. With an increase in the sample dimensions, the mechanical properties reflected in the stress–strain curves changed from elastic–brittle to elastic–plastic or ductile, the strength failure criterion changed from linear to nonlinear, and the failure modes changed from fragmentation failure to local structure-controlled failure and then to splitting failure. (3) The surrounding rock mass near the excavation face of underground caverns typically showed a spalling failure mode, mainly affected by the complex structural characteristics and high in situ stresses, i.e., a tensile fracture mechanism characterized by stress–structure coupling. The research findings not only shed new light on the failure mechanisms and size-dependent mechanical behaviors of hard brittle rocks represented by basalt but also further enrich the basic theory and technical methods for multi-scale analyses in geotechnical engineering, which could provide a reference for the design optimization, construction scheme formulation, and disaster prevention of deep engineering projects. Full article

The focus of this paper is the efficient numerical solution of the fluid flow in the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) reservoir. In this study, the public data available for Discrete Fracture Networks (DFN) around well 58-32 is used to represent the DFN. In this research, a novel computationally efficient method called Hele-Shaw (HS) approximation is used for modeling fluid flow in FORGE well. An analysis of the influence of fracture intensity in a network is carried out using the HS method. The HS method was validated by solving the full Navier–Stokes equations (NSE) for a network of eight fractures. A good agreement was observed between the evaluated results (average deviation of 0.76%). Full article

In order to fully explore the development degree and distribution law of the structural plane of a tunnel surrounding rock in three-dimensional space, this paper studies the geometric characteristic parameters of a structural plane in the study area through field investigation, data acquisition and statistical analysis. The structural plane is divided into three dominant groups by using DIPS. v5. 103 software. The probability distribution model of occurrence, trace length, diameter and spacing of the structural plane is established. This paper focuses on the error correction of structural plane occurrence and the estimation of average trace length based on the rectangular window method. The discrete fracture network model is generated by using MATLAB R2021b software, and the discrete fracture network model is verified from three aspects: structural plane occurrence, average trace length and area density. The verification results are compared with the measured data, and the simulation results are in line with the actual situation on site. Based on the discrete fracture network model, the volume joint number of rock mass is calculated. Based on the JSR index, BQ classification method and RQD classification, the development degree of fractures and surrounding rock classification in this area are evaluated. A method of surrounding rock classification based on three evaluation indexes is discussed to comprehensively and accurately classify the quality of rock mass in this area. Full article

This study performed a representative elementary volume (REV) and 3D equivalent continuum study of rock fractures based on fluid simulations of 3D discrete fracture networks (DFNs). A series of 3D DFNs with heterogeneous aperture distributions (the DFN-H model) and uniform apertures (the DFN-I model) were established, in which the fractures were oriented according to the geological field mapping of a high-level radioactive waste candidate site in China. The 3D DFNs of the different model sizes were extracted and rotated in a number of directions to check whether there was a tensor quality of the permeability at a certain scale. The results show that aperture heterogeneity increases the REV size and results in a necessarily larger model size to reach an equivalent continuum behavior, and this effect is more obvious when the fracture density is smaller. The shape of the 2D permeability contour is irregular, with some breaks when the model size is small. As the model size increases, its shape gradually tends to become smooth and approaches an ellipse. The shape of the permeability contours of the DFN-H model is slender compared to the DFN-I model, indicating a larger difference between the minimum and maximum values of the permeability. For the DFN-H model, there is no appropriate approximation for the equivalent permeability tensor over the studied model size range, whereas a good fit of the permeability ellipsoid is obtained for the DFN-I model, and the 3D directional permeability is calculated at this model scale. The corresponding magnitude and direction of the principal permeability are obtained, which can be viewed as the equivalent permeability tensor for the approximated continuum medium. Full article

In underground construction projects, traversing hard rock layers demands concentrated CO₂ fracturing energy and precise directional crack expansion. Due to the discontinuity of the rock mass at the tip of prefabricated directional fractures in CO₂ fracturing, traditional simulations assuming continuous media are limited. It is challenging to set boundary conditions for high strain rate and large deformation processes. The dynamic expansion mechanism of the 3D fracture network in CO₂ directional fracturing is not yet fully understood. By treating CO₂ fracturing stress waves as hemispherical resonance waves and using a particle expansion loading method along with dynamic boundary condition processing, a 3D numerical model of CO₂ fracturing is constructed. This model analyzes the dynamic propagation mechanism of 3D spatial fractures network in CO₂ directional fracturing rock materials. The results show that in undirected fracturing, the fracture network relies on the weak structures near the rock borehole, whereas in directional fracturing, the crack propagation is guided, extending the fracture’s range. Additionally, the tip of the directional crack is vital for the re-expansion of the rock mass by high-pressure CO₂ gas, leading to the formation of a symmetrical, umbrella-shaped structure with evenly developed fractures. The findings also demonstrate that the discrete element method (DEM) effectively reproduces the dynamic fracture network expansion at each stage of fracturing, providing a basis for studying the CO₂ directional rock cracking mechanism. Full article

Rock quality designation (RQD) is a parameter that describes rock mass quality in terms of percentage recovery of core pieces greater than 10 cm. The RQD represents a basic element of several classification systems. This paper studies scale effects for RQD measurements using synthetic rock masses generated using discrete fracture network (DFN) models. RQD measurements are performed for rock masses with varying fracture intensities and by changing the orientation of the simulated boreholes to account for orientation bias. The objective is to demonstrate the existence of a representative elementary length (REL, 1D analogue of a 3D representative elementary volume, or REV) above which RQD measurements would represent an average indicator of rock mass quality. For the synthetic rock masses, RQD measurements were calculated using the relationship proposed by Priest and Hudson and compared to the simulated RQD measurements along the boreholes. DFN models generated for a room-and-pillar mine using mapped field data were then used as an initial validation, and the conclusion of the study was further validated using the RQD calculation results directly obtained from the depth data collected at an iron cap deposit. The relationship between rock mass scale and assumed threshold length used to calculate RQD is also studied. Full article

To investigate the macroscopic mechanical properties and failure evolution mechanism of desulfurization gypsum–fly ash fluid lightweight soil, a microscale numerical model using PFC2D (Particle Flow Code) was constructed. Uniaxial compression tests were conducted to determine the microscopic parameters of the model, extracting information on the discrete fracture network type, quantity, age, and particle displacement trend. The crack morphology and propagation evolution of desulfurization gypsum–fly ash fluid lightweight soil were explored, and the destructive properties of desulfurization gypsum–fly ash fluid lightweight soil material were evaluated through energy indicators. The research findings suggest that the discrete element numerical model effectively simulates the stress–strain curve and failure characteristics of materials. Under uniaxial compression conditions, microcracks dominated by shear failure occur in the initial loading stage of desulfurization gypsum–fly ash fluid lightweight soil, with a through crack dominated by tensile failure appearing once the load exceeds the peak stress. The dissipated energy evolution in the flow state of desulfurization gypsum–fly ash fluid lightweight soil is relatively gentle, leading to delayed cracking after surpassing the peak stress point. Full article

Multi-cluster fracturing has emerged as an effective technique for enhancing the productivity of deep shale reservoirs. The presence of natural bedding planes in these reservoirs plays a significant role in shaping the evolution and development of multi-cluster hydraulic fractures. Therefore, conducting detailed research on the propagation mechanisms of multi-cluster hydraulic fractures in deep shale formations is crucial for optimizing reservoir transformation efficiency and achieving effective development outcomes. This study employs the finite discrete element method (FDEM) to construct a comprehensive three-dimensional simulation model of multi-cluster fracturing, considering the number of natural fractures present and the geo-mechanical characteristics of a target block. The propagation of hydraulic fractures is investigated in response to the number of natural fractures and the design of the multi-cluster fracturing operations. The simulation results show that, consistent with previous research on fracturing in shale oil and gas reservoirs, an increase in the number of fracturing clusters and natural fractures leads to a larger total area covered by artificial fractures and the development of more intricate fracture patterns. Furthermore, the present study highlights that an escalation in the number of fracturing clusters results in a notable reduction in the balanced expansion of the double wings of the main fracture within the reservoir. Instead, the effects of natural fractures, geo-stress, and other factors contribute to enhanced phenomena such as single-wing expansion, bifurcation, and the bending of different main fractures, facilitating the creation of complex artificial fracture networks. It is important to note that the presence of natural fractures can also significantly alter the failure mode of artificial fractures, potentially resulting in the formation of small opening shear fractures that necessitate careful evaluation of the overall renovation impact. Moreover, this study demonstrates that even in comparison to single-cluster fracturing, the presence of 40 natural main fractures in the region can lead to the development of multiple branching main fractures. This finding underscores the importance of considering natural fractures in deep reservoir fracturing operations. In conclusion, the findings of this study offer valuable insights for optimizing deep reservoir fracturing processes in scenarios where natural fractures play a vital role in shaping fracture development. Full article

Geothermal energy plays a crucial role in the large-scale deep decarbonisation process and the transition of energy structure in our country. Due to the complex reservoir environment of geothermal energy, characterised by low porosity and permeability, conventional fracturing methods struggle to create a complex network of fractures. Temporary plugging and diverting fracturing technology (TPDF) is a key technology to improve the efficiency of geothermal reservoir extraction. However, there is still a lack of knowledge about the migration and sealing law of temporary plugging agents in complex fractures. Therefore, in this study, two multiphase flow models of temporary plugging particle transport at the fracture slit and inside the complex fracture were established by using a Computational Fluid Dynamics (CFD)-Discrete Element Method (DEM) algorithm. The influence of fracturing fluid concentration, temperature, the concentration of temporary plugging particles, and particle size combinations on migration blocking in fractures was investigated. The simulation results indicate the following: High-viscosity fracturing fluid may cause plugging particles to adhere to each other to form clusters of plugging particles, reducing dispersion during transport and slowing down the velocity of the plugging particles. A particle concentration that is too high does not have a better temporary plugging effect. The use of different combinations of particle sizes is significantly better than using a single particle size, which is a key factor for the success of fracture plugging. The research findings are of great theoretical and practical significance for scaled-up, vibration-controlled fracturing technology in geothermal reservoirs. Full article

Fractures are widely distributed in karst areas, and when flow rates are high, they exhibit complex nonlinear behavior that cannot be accurately described by Darcy’s law. In this work, a hydro-mechanical coupling model based on a discrete fracture network is proposed to predict tunnel water inflow, accounting for the impact of non-Darcy flow. The model’s feasibility has been validated by comparing it with experimental results and the field measurements of flow rates at the Bodaoling Tunnel in Guizhou, China. The results show that Darcy flow tends to overestimate water inflow by approximately 25% compared to non-Darcy flow. The non-Darcy effect grows with the increase in initial fracture width and empirical constant q. When q exceeds 8.77 × 10⁻⁶, the growth rate of the Forchheimer number along the fracture width slowed down, and the inhibitory effect of non-Darcy flow on flow became gentle. Additionally, in a complex fracture network, the inflow rate limited by non-Darcy flow at one point drives the water flow through a connect fracture to another point, which increases the difficulty in water inflow prediction. This work highlights the importance of non-Darcy flow and fracture networks when accurately predicting water inflow in tunnels. Full article

Deformation monitoring data provide a direct representation of the structural behavior of reservoir bank rock slopes, and accurate deformation prediction is pivotal for slope safety monitoring and disaster warning. Among various deformation prediction models, hybrid models that integrate field monitoring data and numerical simulations stand out due to their well-defined physical and mechanical concepts, and their ability to make effective predictions with limited monitoring data. The predictive accuracy of hybrid models is closely tied to the precise determination of rock mass mechanical parameters in structural numerical simulations. However, rock masses in rock slopes are characterized by intersecting geological structural planes, resulting in reduced strength and the creation of multiple fracture flow channels. These factors contribute to the heterogeneous, anisotropic, and size-dependent properties of the macroscopic deformation parameters of the rock mass, influenced by the coupling of seepage and stress. To improve the predictive accuracy of the hybrid model, this study introduces the theory of equivalent continuous media. It proposes a method for determining the equivalent deformation parameters of fractured rock mass considering the coupling of seepage and stress. This method, based on a discrete fracture network (DFN) model, is integrated into the hybrid prediction model for rock slope deformation. Engineering case studies demonstrate that this approach achieves a high level of prediction accuracy and holds significant practical value. Full article

The potential of unconventional hydrocarbon resources has been unlocked since the hydraulic fracturing technique in combination with long horizontal wells was applied to develop this type of reservoir economically. The design and optimization of the fracturing treatment and the stimulated reservoir volume and the forecasting of production performance are crucial for the development and management of such resources. However, the production performance of tight gas reservoirs is a complicated nonlinear problem, described by many parameters loaded with uncertainty. The complexity of the problem influences and inspires the sophistication of the solution to be used. This paper proposed an artificial network model that allows for fast, extended, and accurate analyses of the production performance of multiple fractured unconventional tight gas wells. In the comprehensive approach developed, the reservoir rock parameters, the drainage area, and the hydraulic fracture parameters are treated as a variable input to the model. The analysis is no longer constrained by fixed “shoes box” geometry, and the values of the parameters defining the reservoir and stimulated volume are not limited to a few discrete values. The numerical experiment used to construct a database for model development was designed using a genetically optimized Latin hypercube sampling technique. A special approach was used in the preparation of “blind data”, which are crucial for truly reliable model verification. In the result, a developed tool offers an extended rock-fluid description, flexible model, and stimulated reservoir volume dimensioning and parameterization, as well as a high degree of applicability in sensitivity analysis and/or optimization. Full article

Understanding subsurface natural fracture systems is crucial to characterize well production dynamics and long-term productivity potential. In addition, the placement of future wells can benefit from in-depth fracture network connectivity investigations, vastly improving new wells’ profitability and life cycles if they are placed in dense, well-connected natural fracture zones. In this study, a novel natural fracture calibration workflow is proposed. This workflow starts with the extraction of sector geology and a natural fracture model from the pre-built full-field model. Then, a cross wellbore discrete fracture network (CW-DFN) is created using a novel CW-DFN generation tool, based on image log data. An innovative fracture network identification tool is developed to detect the interconnected regional fracture network (IcFN) with CW-DFN. The non-intrusive embedded discrete fracture model (EDFM) is utilized to numerically incorporate the complex IcFN and CW-DFN in a reservoir simulation, and it is history-matched by tuning their conductivities. This workflow is applied to a single vertical well within a natural fracture carbonate reservoir in Northwest China. The study results show that the number of CW-DFNs is 11, and the number of IcFNs is 72. The non-intersected natural fractures only account for 5.5% of the production, and thus can be removed to improve simulation efficiency. The history-matching absolute average relative deviation (AARD) is 15.16%. The calibrated effective fracture permeability is 280 millidarcy, with an aperture of 0.001 m, equating to a conductivity of 0.28 millidarcy-meter. The 30-year gas production forecast is estimated to be 1.66 billion cubic meters based on a history-matched model. Finally, if the well is drilled to the east of the sector, 30-year production declines to 1.33 billion cubic meters (a reduction of 20%). However, if the well is drilled to the west of the sector, 30-year production increases to 2 billion cubic meters (an improvement of 20.5%). Full article

As a geological barrier for high-level radioactive waste (HLW) disposal in China, granite is crucial for blocking nuclide migration into the biosphere. However, the high uncertainty associated with the 3D geological system, such as the stochastic discrete fracture networks in granite, significantly impedes practical safety assessments of HLW disposal. This study proposes a Monte Carlo simulation (MCS)-based simulation framework for evaluating the long-term barrier performance of nuclide migration in fractured rocks. Statistical data on fracture geometric parameters, on-site hydrogeological conditions, and relevant migration parameters are obtained from a research site in Northwestern China. The simulation models consider the migration of three key nuclides, Cs-135, Se-79, and Zr-93, in fractured granite, with mechanisms including adsorption, advection, diffusion, dispersion, and decay considered as factors. Subsequently, sixty MCS realizations are performed to conduct a sensitivity analysis using the open-source software OpenGeoSys-5 (OGS-5). The results reveal the maximum and minimum values of the nuclide breakthrough time T_t (12,000 and 3600 years, respectively) and the maximum and minimum values of the nuclide breakthrough concentration C_max (4.26 × 10⁻⁴ mSv/a and 2.64 × 10⁻⁵ mSv/a, respectively). These significant differences underscore the significant effect of the uncertainty in the discrete fracture network model on long-term barrier performance. After the failure of the waste tank (1000 years), nuclides are estimated to reach the outlet boundary 6480 years later. The individual effective dose in the biosphere initially increases and then decreases, reaching a peak value of C_max = 4.26 × 10⁻⁴ mSv/a around 350,000 years, which is below the critical dose of 0.01 mSv/a. These sensitivity analysis results concerning nuclide migration in discrete fractured granite can enhance the simulation and prediction accuracy for risk evaluation of HLW disposal. Full article

In the process of the large-scale hydraulic fracturing of a shale gas field in the Weiyuan area of Sichuan province, the quantitative description and evaluation of hydraulic fracture expansion morphology and the three-dimensional distribution law are the key points of evaluation of block fracturing transformation effect. Many scholars have used the finite element method, discrete element method, grid-free method and other numerical simulation methods to quantitatively characterize hydraulic fractures, but there are often the problems that the indoor physical simulation results are much different from the actual results and the accuracy of most quantitative studies is poor. Considering rock mechanics parameters and based on the displacement discontinuity method (DDM), a single-stage multi-cluster fracture propagation model of horizontal well was established. The effects of Young’s modulus, Poisson’s ratio, the in situ stress difference, the approximation angle, the perforation cluster number and the perforation spacing on the formation of complex fracture networks and on the geometrical parameters of hydraulic fractures were simulated. The research results can provide theoretical reference and practical guidance for the optimization of large-scale fracturing parameters and the quantitative post-fracturing evaluation of horizontal wells in unconventional reservoirs such as shale gas reservoirs. Full article