Search results

The purpose of this paper is to extend the bridge scale method (BSM) developed for granular materials with only the solid phase to that taking into account the effects of wetting process in porous continuum. The granular material is modeled as partially saturated porous Cosserat continuum and discrete particle assembly in the coarse and fine scales, respectively.

Based on the mass and momentum conservation laws for the three phases, i.e. the solid skeleton, the pore water and the pore air, the governing equations for the unsaturated porous Biot-Cosserat continuum model in the coarse scale are derived. In light of the passive air pressure assumption, a reduced finite element model for the model is proposed. According to the decoupling of the fine and coarse scale calculations in the BSM, the unsaturated porous Cosserat continuum model using the finite element method and the discrete element model using the discrete element method for granular media are combined.

The numerical results for a 2D example problem of slope stability subjected to increasing rainfall along with mechanical loading demonstrate the applicability and performance of the present BSM. The microscopic mechanisms of macroscopic shear band developed in the slope are demonstrated.

Do not account for yet the effects of unsaturated pore water in the fine scale.

The novel BSM that couples the Biot-Cosserat porous continuum modeling and the discrete particle assembly modeling in both coarse and fine scales, respectively, is proposed to provide a micro-macro discrete-continuum two-scale modeling approach for numerical simulations of the hydro-mechanical coupling problems in unsaturated granular materials.

The purpose of this paper is to perform the thermo-hydro-mechanical (THM) numerical analysis in order to study the thawing process for frozen soil and to predict the thawing settlement.

A new one-dimensional multi-field physical coupled model was proposed here to describe the thawing process of saturated frozen soil, whereby the void ratio varied linearly with effective stress (Eq. 10) and hydraulic conductivity (Eq. 27b). The thawing process was simulated with different initial and boundary conditions in an open system with temperature variations. The mechanical behavior and water migration of the representative cases were also investigated.

The comparisons of representative cases with experimental data demonstrated that the model predicts thawing settlement well. It was found that the larger temperature gradient, higher overburden pressure and higher water content could lead to larger thawing settlement. The temperature was observed that to distribute height linearly in both frozen zone and unfrozen zone of the sample. Water migration forced to a decrease in the water content of the unfrozen zone and an increase in water content at the thawing front.

In this study, only the one-directional thawing processes along the frozen soil samples were investigated numerically and compared with test results, which can be extended to two-dimensional analysis of thawing process in frozen soil.

This study helps to understand the thawing process of frozen soil by coupled thermo-hydro-mechanical numerical simulation.

Large deformation problems are frequently encountered in various fields of geotechnical engineering. The particle finite element method (PFEM) has been proven to be a promising method to solve large deformation problems. This study aims to develop a computational framework for modelling the hydro-mechanical coupled porous media at large deformation based on the PFEM.

The PFEM is extended by adopting the linear and quadratic triangular elements for pore water pressure and displacements. A six-node triangular element is used for modelling two-dimensional problems instead of the low-order three-node triangular element. Thus, the numerical instability induced by volumetric locking is avoided. The Modified Cam Clay (MCC) model is used to describe the elasto-plastic soil behaviour.

The proposed approach is used for analysing several consolidation problems. The numerical results have demonstrated that large deformation consolidation problems with the proposed approach can be accomplished without numerical difficulties and loss of accuracy. The coupled PFEM provides a stable and robust numerical tool in solving large deformation consolidation problems. It is demonstrated that the proposed approach is intrinsically stable.

The PFEM is extended to consider large deformation-coupled hydro-mechanical problem. PFEM is enhanced by using a six-node quadratic triangular element for displacement and this is coupled with a four-node quadrilateral element for modelling excess pore pressure.

Natural gas hydrate (NGH) has been regarded as one of the most important resources due to NGH's large amounts of reserve. However, NGH development still faces many technical challenges, such as low production rate and reservoir instability resulting from NGH decomposition. Therefore, developing a fully coupled THMC model for simulating the hydrate decomposition and studying its mechanical behavior is very important and necessary. The purpose of this article is to develop and solve a multi-phase, strong nonlinearity and large-scale fully coupled thermal-hydro-mechanical–chemical (THMC) model for simulating the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and rock deformation during NGH decomposition.

In this paper, a multi-phase, strong nonlinearity and large-scale fully coupled THMC model is developed for simulating the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and rock deformation during NGH dissociation. The fully coupled THMC model is solved by using a fully implicit finite element method, in which the gas pressure, water pressure, temperature and displacement are taken as basic unknown variables. The proposed model is validated against with the experimental data, showing high accuracy and reliability.

A multi-phase, strong nonlinearity and large-scale fully coupled THMC model is developed for simulating the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and rock deformation during NGH decomposition. The proposed model is validated against with the experimental data, showing high accuracy and reliability.

Some assumptions are made to make the model tractable, including (1) the composition gas of hydrate is pure methane; (2) the gas-liquid multi-phase flow in the pore obeys Darcy's law; (3) hydrate occurs on the surface of soil particles, both of them form the composite consolidation material; (4) the small-strain assumption is applied to composite solid materials, which are treated as skeletons and cannot be moved; (5) momentum change caused by phase change is not considered.

NGH has been regarded as one of the most important resources due to its large amounts of reserve. However, NGH development still faces many technical challenges, such as low production rate and reservoir instability resulting from NGH decomposition. Most of the existing studies decouple the process with solid deformation and seepage behavior, but the accuracy of the numerical results will be sacrificed to certain extent. Therefore, it is very important and necessary to develop a fully coupled THMC model for simulating the hydrate decomposition and studying its mechanical behavior.

NGH, widely distributed in shallow seabed or permanent frozen region, has the characteristics of high energy density and high combustion efficiency (Yan et al., 2020). A total of around 7.5 × 1,018 m³ has been proved to exist around the world and 1 m³ of NGH can release about 160–180 m³ of natural gas (Kvenvolden and Lorenson) under normal conditions. Safely and sustainably extracting NGH commercially can effectively relieve global energy pressure and contribute to achieving carbon reduction goals.

The novelty of the present work lies in mainly two aspects. First, a fully coupled THMC model is developed for studying the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and solid deformation during NGH dissociation. Second, the numerical solution is obtained by using a fully implicit finite element method (FEM) and is validated against experimental data.

The objective of this paper is to provide a better understanding of the effect of irregular columnar jointed structure on the permeability and flow characteristics of rock masses.

An efficient numerical procedure is proposed to investigate the permeability and fluid flow in columnar jointed rock masses (CJRMs), of which the columnar jointed networks are generated by a modified constrained centroid Voronoi algorithm according to the field statistical results. The fractures are represented explicitly by using the lower-dimensional zero thickness elements. And the modeling scheme is validated by a benchmark test for flow in fractured porous media. The effective permeability and representative elementary volume (REV) size of CJRMs are estimated using finite element method (FEM). The influences of joint density and variation coefficient of columnar joint structure on the permeability of the rock mass are discussed.

The simulation results indicate that the permeability is scale-dependent and tends to be stable with increase of model size. The hydraulic REV size is determined as 3.5 m for CJRMs in the present study. Moreover, the joint density is a dominant factor affecting the permeability of CJRMs. The average permeability of columnar jointed structures increases linearly with the joint density under the same REV size, while the influence from the coefficient of variation can be neglected.

The present paper investigates the REV size of the CJRMs and the effect of joint parameters on the permeability. The proposed method and the results obtained are useful on understanding the hydraulic characteristic of the irregular CJRMs in engineering projects.

Accurately evaluating fluid flow behaviors and determining permeability for deforming porous media is time-consuming and remains challenging. This paper aims to propose a novel machine-learning method for the rapid estimation of permeability of porous media at different deformation stages constrained by hydro-mechanical coupling analysis.

A convolutional neural network (CNN) is proposed in this paper, which is guided by the results of finite element coupling analysis of equilibrium equation for mechanical deformation and Boltzmann equation for fluid dynamics during the hydro-mechanical coupling process [denoted as Finite element lattice Boltzmann model (FELBM) in this paper]. The FELBM ensures the Lattice Boltzmann analysis of coupled fluid flow with an unstructured mesh, which varies with the corresponding nodal displacement resulting from mechanical deformation. It provides reliable label data for permeability estimation at different stages using CNN.

The proposed CNN can rapidly and accurately estimate the permeability of deformable porous media, significantly reducing processing time. The application studies demonstrate high accuracy in predicting the permeability of deformable porous media for both the test and validation sets. The corresponding correlation coefficients (R²) is 0.93 for the validation set, and the R² for the test set A and test set B are 0.93 and 0.94, respectively.

This study proposes an innovative approach with the CNN to rapidly estimate permeability in porous media under dynamic deformations, guided by FELBM coupling analysis. The fast and accurate performance of CNN underscores its promising potential for future applications.

Multi-well hydrofracturing is an important technology to create new fractures and expand existing fractures to increase reservoir permeability. The propagation morphology of the fracture network is affected by the disturbance between the fractures initiation sequences and spacings between adjacent wells. However, it remains unclear how well spacing and initiation sequences lead to fracture propagation, deflection and connection.

In this study, the thermal-hydro-mechanical coupling effect in the hydrofracturing process was considered, to establish a finite element-discrete element model of multistage hydrofracturing in a horizontal well. Using typical cases, the unstable propagation of hydraulic fractures in multiple horizontal wells was investigated under varying well spacing and initiation sequences. Combined with the shear stress shadow caused by in situ stress disturbed by fracture tip propagation, the quantitative indexes of fracture propagation such as length, volume, displacement vector, deflection and unstable propagation behavior of the hydrofracturing fracture network were analyzed.

The results show that the shear stress disturbance caused by multiple hydraulic fractures is a significant factor in multi-well hydrofracturing. Reducing the spacing between multiple wells increases the stress shadow area and aggravates the mutual disturbance and deflection between the fractures. The quantitative analysis results show that a decrease of well spacing reduces the total length of hydraulic fractures but increases the total volume of the fracture; compared with sequential and simultaneous fracturing, alternate fracturing can effectively reduce stress shadow area, alleviate fracture disturbance and generate larger fracture propagation length and volume.

The numerical models and results of the unstable propagation and stress evolution of the hydraulic fracture network under thermal-hydro-mechanical coupling obtained in this study can provide useful guidance for the evaluation and design of rock mass fracture networks in deep unconventional oil and gas reservoirs.

The purpose of this study is to investigate the unstable propagation of parallel hydraulic fractures induced by interferences of adjacent perforation clusters and thermal diffusion. Fracture propagation in the process of multistage fracturing of a rock mass is deflected owing to various factors. Hydrofracturing of rock masses in deep tight reservoirs involves thermal diffusion, fluid flow and deformation of rock between the rock matrix and fluid in pores and fractures.

To study the unstable propagation behaviours of three-dimensional (3D) parallel hydraulic fractures induced by the interferences of adjacent perforation clusters and thermal diffusion, a 3D engineering-scale numerical model is established under different fracturing scenarios (sequential, simultaneous and alternate fracturing) and different perforation cluster spacings while considering the thermal-hydro-mechanical coupling effect. Stress disturbance region caused by fracture propagation in a deep tight rock mass is superimposed and overlaid with multiple fractures, resulting in a stress shadow effect and fracture deflection.

The results show that the size of the stress shadow areas and the interaction between fractures increase with decreasing multiple perforation cluster spacing in horizontal wells. Alternate fracturing can produce more fracture areas and improve the fracturing effect compared with those of sequential and simultaneous fracturing. The larger the temperature gradient between the fracturing fluid and rock matrix, the stronger the thermal diffusion effect, and the effect of thermal diffusion on the fracture propagation is significant.

This study focuses on the behaviours of the unstable dynamic propagation of 3D parallel hydraulic fractures induced by the interferences of adjacent perforation clusters and thermal diffusion. Further, the temperature field affects the fracture deflection requires could be investigated from the mechanisms; this paper is to study the unstable propagation of fractures in single horizontal well, which can provide a basis for fracture propagation and stress field disturbance in multiple horizontal wells.

Hydrofracturing technology has been widely used in tight oil and gas reservoir exploitation, and the fracture network formed by fracturing is crucial to determining the resources recovery rate. Due to the complexity of fracture network induced by the random morphology and type of fluid-driven fractures, controlling and optimising its mechanisms is challenging. This paper aims to study the types of multiscale mode I/II fractures, the fluid-driven propagation of multiscale tensile and shear fractures need to be studied.

A dual bilinear cohesive zone model (CZM) based on energy evolution was introduced to detect the initiation and propagation of fluid-driven tensile and shear fractures. The model overcomes the limitations of classical linear fracture mechanics, such as the stress singularity at the fracture tip, and considers the important role of fracture surface behaviour in the shear activation. The bilinear cohesive criterion based on the energy evolution criterion can reflect the formation mechanism of complex fracture networks objectively and accurately. Considering the hydro-mechanical (HM) coupling and leak-off effects, the combined finite element-discrete element-finite volume approach was introduced and implemented successfully, and the results showed that the models considering HM coupling and leak-off effects could form a more complex fracture network. The multiscale (laboratory- and engineering-scale) Mode I/II fractures can be simulated in hydrofracturing process.

Based on the proposed method, the accuracy and applicability of the algorithm were verified by comparing the analytical solution of KGD and PKN models. The effects of different in situ stresses and flow rates on the dynamic propagation of hydraulic fractures at laboratory and engineering scales were investigated. when the ratio of in situ stress is small, the fracture propagation direction is not affected, and the fracture morphology is a cross-type fracture. When the ratio of in situ stress is relatively large, the propagation direction of the fracture is affected by the maximum in situ stress, and it is more inclined to propagate along the direction of the maximum in situ stress, forming double wing-type fractures. Hydrofracturing tensile and shear fractures were identified, and the distribution and number of each type were obtained. There are fewer hydraulic shear fractures than tensile fractures, and shear fractures appear in the initial stage of fracture propagation and then propagate and distribute around the perforation.

The proposed dual bilinear CZM is effective for simulating the types of Mode I/II fractures and seizing the fluid-driven propagation of multiscale tensile and shear fractures. Practical fracturing process involves the multi-type and multiscale fluid-driven fracture propagation. This study introduces general fluid-driven fracture propagation, which can be extended to the fracture propagation analysis of potential fluid fracturing, such as other liquids or supercritical gases.

Multi-well hydrofracturing is a key technology in engineering, and the evaluation, control and optimization of the fracturing network determine the recovery rate of unconventional oil and gas production. In engineering terms, altering well spacing and perforation initiation sequences changes fracture propagation behavior. Fracture propagation can result in fracture-to-fracture and well-to-well interactions. This may be attributed to the interference between fractures caused by squeezing of the reservoir strata. Meanwhile, the stratal movement caused by the propagation of the fractures may lead to either the secondary fracturing of wells with primary fractures or perforation to begin fracturing. Besides, the stratal compression and squeeze of multi-well hydrofracturing will cause earthquakes; the fracture size is different owing to the different fracturing scenarios, and the occurrence of induced microseismic events is still unknown; microseismic events also affect fracture orientation and deflection. If the mechanism of the above mechanical behavior cannot be clarified, optimizing the fracture network and reduce the induced microseismic disaster becomes difficult.

In this study, combined finite element-discrete element models were used to simulate the multi-well hydrofracturing. Numerical cases compared the fracture network, dynamic stratal movement and microseismic events at 50, 75 and 100 m well spacings, respectively, and varying initiation sequence of multiple horizontal wells.

From the results, fracture propagation in multi-well hydrofracturing may simulate the propagation and deflection of adjacent fractures and induce fracture-to-fracture and well-to-well interactions. As the well spacing increases, the effect of fracturing-induced stratal movement and squeezing deformation decrease. In alternate fracturing, starting from a well located in the middle can effectively reduce the influence of stratal movement on fracturing, and the fracturing scenario with cross-perforation can minimize the influence of stratal movement. The stratal movement between multiple wells is positively correlated to microseismic events, which behaviors can be effectively weakened by reducing the strata movement.

The fracture network, thermal-hydro-mechanical coupling, fracturing-induced stratal movement and microseismic events were analyzed. This study analyzed the intersection and propagation behavior of fractures in multi-well hydrofracturing, which can be used to evaluate and study the mechanism of hydrofracturing fracture network propagation in multiple horizontal wells and conduct fracture optimization research to form an optimized hydrofracturing scheme by reasonably arranging the spacing between wells and initiation sequences of perforation clusters.