Search results

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 simulate the tensile and shear types of fractures using the mixed fracture criteria considering the energy evolution based on the dual bilinear cohesive zone model and investigate the dynamic propagation of tensile and shear fractures induced by an impact load in rock. The propagation of tension and shear at different scales induced by the impact load is also an important aspect of this study.

In this study, based on the well-developed dual bilinear cohesive zone model and combined finite element-discrete element method, the dynamic propagation of tensile and shear fractures induced by the impact load in rock is investigated. Some key technologies, such as the governing partial differential equations, fracture criteria, numerical discretisation and detection and separation, are introduced to form the global algorithm and procedure. By comparing with the tensile and shear fractures induced by the impact load in rock disc in typical experiments, the effectiveness and reliability of the proposed method are well verified.

The dynamic propagation of tensile and shear fractures in the laboratory- and engineering-scale rock disc and rock strata are derived. The influence of mesh sensitivity, impact load velocities and load positions are investigated. The larger load velocities may induce larger fracture width and entire failure. When the impact load is applied near the left support constraint boundary, concentrated shear fractures appear around the loading region, as well as induced shear fracture band, which may induce local instability. The proposed method shows good applicability in studying the propagation of tensile and shear fractures under impact loads.

The proposed method can identify fracture propagation via the stress and energy evolution of rock masses under the impact load, which has potential to be extended into the investigation of the mixed fractures and disturbance of in-situ stresses during dynamic strata mining in deep energy development.

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.

With the development of fracturing technology, the research of multi-well hydrofracturing becomes the key issue. Frac-hits in multi-well hydrofracturing has an important effect on fracture propagation and final production of fractured well; in the process of hydrofracturing, there are many implement parameters that can affect frac-hits, and previous studies in this area have not systematically targeted the influence of a single parameter on multi-well hydrofracturing. Therefore, it is of great significance to study the occurrence rule and influence of frac-hits for optimizing the design of fracturing wells.

Based on the proposed numerical models, the effects of different fracturing implement parameters (perforation cluster spacing, well spacing and injection rate) on frac-hits are compared in numerical cases. Through the analysis of fracture network, stress field and microseismic, the effects of different fracturing implement parameters on frac-hits and connections are compared.

The simulation results show that the effect of perforation cluster spacing and well spacing on frac-hits is greater than that of injection rate. Smaller well spacing makes it easier for fractures between adjacent wells to interact with each other, which increases the risk of frac-hits and reduces the risk of fracture connections. Smaller perforation cluster spacing results in larger individual fracture lengths and greater deflection angles, which makes the possibility of frac-hits and connections greater. The lower the injection rate, the lower the probability of frac-hits.

In this study, the influence of different fracturing implement parameters on frac-hits and connections in multi-well hydrofracturing is studied, and the mechanism of frac-hits and connections is analyzed through fracture network, stress field and microseismic analysis. Different simulation results are compared to optimize fracturing well parameter design and provide reference for engineering application.

Thermal fractures initiated under cooling at the surfaces of a 2-D or 3-D structure propagate, arrest and coalesce, leading to its structural failure and material-property changes, while the same processes can happen in the rock mass between parallel hydraulic fractures filled with cold fluid, leading to enhanced fracture connectivity and permeability.

This study used a 2-D plane strain fracture model for mixed-mode thermal fractures from two parallel cooling surfaces. Fracture propagation was governed by the theory of linear elastic fracture mechanics, while the displacement and temperature fields were discretized using the adaptive finite element method. This model was validated using two numerical benchmarks with strong fracture curvature and then used to simulate the propagation and coalescence of thermal fractures in a long rock mass.

Modeling results show two regimes: (1) thermal fractures from a cooling surface propagate and arrest by following the theoretical solutions of half-plane fractures before the unfractured portion decreases to 20% rock-mass width and (2) some pairs of fractures from the opposite cooling surfaces tend to eventually coalesce. The fracture coalescence time is in a power law with rock-mass width.

These findings are relevant to both subsurface engineering and material engineering: structure failure is a key concern in the latter, while fracture coalescence can enhance the connectivity of thermal and hydraulic fractures and thus reservoir permeability in the former.

This paper aims to employ a multivariate nonlinear regression analysis to establish a predictive model for the final fracture area, while accounting for the impact of individual parameters.

This analysis is based on the numerical simulation data obtained, using the hybrid finite element–discrete element (FE–DE) method. The forecasting model was compared with the numerical results and the accuracy of the model was evaluated by the root mean square (RMS) and the RMS error, the mean absolute error and the mean absolute percentage error.

The multivariate nonlinear regression model can accurately predict the nonlinear relationships between injection rate, leakoff coefficient, elastic modulus, permeability, Poisson’s ratio, pore pressure and final fracture area. The regression equations obtained from the Newton iteration of the least squares method are strong in terms of the fit to the six sensitive parameters, and the model follow essentially the same trend with the numerical simulation data, with no systematic divergence detected. Least absolutely deviation has a significantly weaker performance than the least squares method. The percentage contribution of sensitive parameters to the final fracture area is available from the simulation results and forecast model. Injection rate, leakoff coefficient, permeability, elastic modulus, pore pressure and Poisson’s ratio contribute 43.4%, −19.4%, 24.8%, −19.2%, −21.3% and 10.1% to the final fracture area, respectively, as they increased gradually. In summary, (1) the fluid injection rate has the greatest influence on the final fracture area. (2)The multivariate nonlinear regression equation was optimally obtained after 59 iterations of the least squares-based Newton method and 27 derivative evaluations, with a decidability coefficient R2 = 0.711 representing the model reliability and the regression equations fit the four parameters of leakoff coefficient, permeability, elastic modulus and pore pressure very satisfactorily. The models follow essentially the identical trend with the numerical simulation data and there is no systematic divergence. The least absolute deviation has a significantly weaker fit than the least squares method. (3)The nonlinear forecasting model of physical parameters of hydraulic fracturing established in this paper can be applied as a standard for optimizing the fracturing strategy and predicting the fracturing efficiency in situ field and numerical simulation. Its effectiveness can be trained and optimized by experimental and simulation data, and taking into account more basic data and establishing regression equations, containing more fracturing parameters will be the further research interests.

The nonlinear forecasting model of physical parameters of hydraulic fracturing established in this paper can be applied as a standard for optimizing the fracturing strategy and predicting the fracturing efficiency in situ field and numerical simulation. Its effectiveness can be trained and optimized by experimental and simulation data, and taking into account more basic data and establishing regression equations, containing more fracturing parameters will be the further research interests.

The interaction between rock mass structural planes and dynamic stress levels is important to determine the stability of rock mass structures in underground geotechnical engineering. In this work, the authors aim to focus on the degradation effects of fracture geometric parameters and unloading stress paths on rock mechanical properties.

A three-dimensional Particle Flow Code (PFC3D) was used for a systematic numerical simulation of the strength failure and cracking behavior of granite specimens containing prefabricated cracks under conventional triaxial compression and triaxial unilateral unloading. The authors demonstrated the unique mechanical response of prefabricated fractured rock under two conditions. The crack initiation, propagation, and coalescence process of pre-fissured specimens were analyzed in detail.

The authors show that the prefabricated cracks and unilateral unloading conditions not only deteriorate the mechanical strength but also have significant differences in failure modes. The degrading effect of cracks on model strength increases linearly with the decrease of the dip angle. Under the condition of true triaxial unilateral unloading, the deterioration effect of peak strength of rock is very significant, and unloading plays a role in promoting the instability failure of rock after peak, making the rock earlier instability failure. Associating with the particle vector diagram and crack coalescence process, the authors find that model failure mode under unilateral loading conditions is obviously distinct from that in triaxial loading. The peak strain in the unloading direction increases sharply, resulting in a new shear slip.

This study is expected to improve the understanding of the strength failure and cracking behavior of fractured rock under unilateral unloading.

The determination of mode-I fracture toughness of brittle structural materials by means of the notched Brazilian disc configuration is studied. Advantage is taken of a recently introduced analytical solution and, also, of data provided by an experimental protocol with notched marble specimens under diametral compression using the loading device suggested by International Society for Rock Mechanics (ISRM) and also the three-dimensional digital image correlation (3D-DIC) technique.

The analytical solution highlighted the role of geometrical factors, like, for example, the width of the notch, which are usually disregarded. The data of the experimental protocol were comparatively considered with those concerning the response of the specific material under uniaxial tensile load.

This combined study provided interesting data concerning some open issues, as it is the exact crack initiation point and the level of the critical load causing crack initiation. It was definitely indicated that the crack initiation point is not a priori known (even for notched specimens) and, also, that the maximum recorded load does not correspond by default to the critical load responsible for the onset of catastrophic macroscopic fracture.

It was suggested that the load considered critical one for the determination of mode-I fracture toughness K_IC is erroneous. At a load equal to about 70% of the maximum one, a process zone is formed (zone of non-reversible phenomena) around the notch's crown, designating termination of the validity of any linear elastic solution used to determine the normalized stress intensity factors (SIFs). Moreover, at a load level equal to about 95% of the macroscopically observed fracture load, crack propagation has already begun. Therefore, the experimental procedure must be monitored with additional equipment, providing an overview of the displacement field developed during loading.

The main aim of the current paper is to find a numerical plan for hydraulic fracturing problem with application in extracting natural gases and oil.

First, time discretization is accomplished via Crank-Nicolson and semi-implicit techniques. At the second step, a high-order finite element method using quadratic triangular elements is proposed to derive the spatial discretization. The efficiency and time consuming of both obtained schemes will be investigated. In addition to the popular uniform mesh refinement strategy, an adaptive mesh refinement strategy will be employed to reduce computational costs.

Numerical results show a good agreement between the two schemes as well as the efficiency of the employed techniques to capture acceptable patterns of the model. In central single-crack mode, the experimental results demonstrate that maximal values of displacements in x- and y- directions are 0.1 and 0.08, respectively. They occur around both ends of the line and sides directly next to the line where pressure takes impact. Moreover, the pressure of injected fluid almost gained its initial value, i.e. 3,000 inside and close to the notch. Further, the results for non-central single-crack mode and bifurcated crack mode are depicted. In central single-crack mode and square computational area with a uniform mesh, computational times corresponding to the numerical schemes based on the high order finite element method for spatial discretization and Crank-Nicolson as well as semi-implicit techniques for temporal discretizations are 207.19s and 97.47s, respectively, with 2,048 elements, final time T = 0.2 and time step size τ = 0.01. Also, the simulations effectively illustrate a further decrease in computational time when the method is equipped with an adaptive mesh refinement strategy. The computational cost is reduced to 4.23s when the governed model is solved with the numerical scheme based on the adaptive high order finite element method and semi-implicit technique for spatial and temporal discretizations, respectively. Similarly, in other samples, the reduction of computational cost has been shown.

This is the first time that the high-order finite element method is employed to solve the model investigated in the current paper.

This paper aims to provide SIF calculation method for engineering application.

In this paper, the stress intensity factors (SIFs) calculation method is applied to the anisotropic Ni-based single crystal film cooling holes (FCHs) structure.

Based on contour integral, the anisotropic SIFs analysis finite element method (FEM) in Ni-based single crystal is proposed. The applicability and mesh independence of the method is assessed by comparing the calculated SIFs using mode of plate with an edge crack. Anisotropic SIFs can be calculated with excellent accuracy using the finite element contour integral approach. Then, the effect of crystal orientation and FCHs interference on the anisotropic SIFs is clarified. The SIFs of FCH edge crack in the [011] orientated Ni-based single crystal increases faster than the other two orientations. And the SIF of horizontal interference FCHs edge crack is also larger than that of the inclined interference one.

The SIFs of the FCH edge crack in the turbine air-cooled blade are innovatively computed using the sub-model method. Both the Mode I and II SIFs of FCHs edge crack in blade increase with crack growing.