Search results

The microseismic monitoring technique has great advantages on identifying the location, extent and the mechanism of damage process occurring in rock mass. This study aims to analyze distribution characteristics and the evolution law of excavation damage zone of surrounding rock based on microseismic monitoring data.

In situ test using microseismic monitoring technique is carried out in the large-span transition tunnel of Badaling Great Wall Station of Beijing-Zhangjiakou high-speed railway. An intelligent microseismic monitoring system is built with symmetry monitoring point layout both on the mountain surface and inside the tunnel to achieve three-dimensional and all-round monitoring results.

Microseismic events can be divided into high density area, medium density area and low density area according to the density distribution of microseismic events. The positions where the cumulative distribution frequencies of microseismic events are 60 and 80% are identified as the boundaries between high and medium density areas and between medium and low density areas, respectively. The high density area of microseismic events is regarded as the high excavation damage zone of surrounding rock, which is affected by the grade of surrounding rock and the span of tunnel. The prediction formulas for the depth of high excavation damage zone of surrounding rock at different tunnel positions are given considering these two parameters. The scale of the average moment magnitude parameters of microseismic events is adopted to describe the damage degree of surrounding rock. The strong positive correlation and multistage characteristics between the depth of excavation damage zone and deformation of surrounding rock are revealed. Based on the depth of high excavation damage zone of surrounding rock, the prestressed anchor cable (rod) is designed, and the safety of anchor cable (rod) design parameters is verified by the deformation results of surrounding rock.

The research provides a new method to predict the surrounding rock damage zone of large-span tunnel and also provides a reference basis for design parameters of prestressed anchor cable (rod).

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.

Optimized three-dimensional (3D) fracture networks are crucial for multistage hydrofracturing. To better understand the mechanisms controlling potential disasters as well as to predict them in 3D multistage hydrofracturing, some governing factors, such as fluid injection-induced stratal movement, compression between multiple hydraulic fractures, fracturing fluid flow, fracturing-induced microseismic damaged and contact slip events, must be properly simulated via numerical models. This study aims to analyze the stratal movement and microseismic behaviours induced by multistage propagation of 3D multiple hydraulic fractures.

Adaptive finite element–discrete element method was used to overcome the limitations of conventional finite element methods in simulating 3D fracture propagation. This new approach uses a local remeshing and coarsening strategy to ensure the accuracy of solutions, reliability of fracture propagation path and computational efficiency. Engineering-scale numerical models were proposed that account for the hydro-mechanical coupling and fracturing fluid leak-off, to simulate multistage propagation of 3D multiple hydraulic fractures, by which the evolution of the displacement, porosity and fracture fields, as well as the fracturing-induced microseismic events were computed.

Stratal movement and compression between 3D multiple hydraulic fractures intensify with increasing proximity to the propagating fractures. When the perforation cluster spaces are very narrow, alternate fracturing can improve fracturing effects over those achieved via sequential or simultaneous fracturing. Furthermore, the number and magnitude of microseismic events are directly proportional to the stratal movement and compression induced by multistage propagation of fracturing fracture networks.

Microseismic events induced by multistage propagation of 3D multiple hydraulic fractures and perforation cluster spaces and fracturing scenarios that impact the deformation and compression among fractures in porous rock matrices are well predicted and analyzed.

Simultaneous hydrofracturing of multiple perforation clusters in vertical wells has been applied in the stimulation of hydrocarbon resources reservoirs. This technology is significantly impeded due to the challenges in its application to the multilayered reservoirs that comprise multiple interlayers. One of the challenges is the accurate understanding and characterization of propagation and deflection of the multiple hydraulic fractures between reservoirs and embedded interlayers.

Numerical models of the tight multilayered reservoirs containing multiple interlayers were established to study hydrofracturing of multiple perforation clusters and its influencing factors on unstable propagation and deflection of hydraulic fractures. Brittle and plastic multilayered reservoirs fully considering the influences of different in situ stress ratio and physical attributes for reservoir and interlayer strata on propagations of hydraulic fractures were investigated. The combined finite element–discrete element method and mesh refinement strategy were adopted to guarantee the accuracy of stress solutions and reliability of fracture path in computation.

Results show that the shear stress fields between adjacent multiple hydraulic fractures are superposed to cause fractures deflection. Stress shadows induce the shielding effects of hydraulic fractures and inhibit fractures growth to emerge unstable propagation behaviors, and a main single fracture and several minor fractures develop. As the in situ stress ratio increases, hydraulic fractures more easily deflect toward the direction of maximum in situ stress, and stress shadow and mutual interaction effects between them are intensified. Compared to brittle reservoir, plastic-enhanced reservoir may limit fracture growth and cannot form long fracture length; nevertheless, plastic properties of reservoir are prone to induce more microseismic events with larger magnitude.

The obtained fracturing behaviors and mechanisms based on engineering-scale multilayered reservoir may provide effective schemes for controlling and estimating the unstable propagation of multiple hydraulic fractures.

This paper aims to propose a deployment optimization and efficient synchronous acquisition method for compressive stress sensors used by stress distribution law research based on the genetic algorithm and numerical simulations. The authors established a new method of collecting the mining compressive stress-strain distribution data to address the problem of the number of sensors and to optimize the sensor locations in physical similarity simulations to improve the efficiency and accuracy of data collection.

First, numerical simulations were used to obtain the compressive stress distribution curve under specific mining conditions. Second, by comparing the mean square error between a fitted curve and simulation data for different numbers of sensors, a genetic algorithm was used to optimize the three-dimensional (3D) spatial deployment of sensors. Third, the authors designed an efficient synchronous acquisition module to allow distributed sensors to achieve synchronous and efficient acquisition of hundreds of data points through a built-in on-board database and a synchronous sampling communication structure.

The sensor deployment scheme was established through the genetic algorithm, A synchronous and selective data acquisition method was established for reduced the amount of sensor data required under synchronous acquisition and improved the system acquisition efficiency. The authors obtained a 3D compressive stress distribution when the advancement was 200 m on a large-scale 3D physical similarity simulation platform.

The proposed method provides a new optimization method for sensor deployment in physical similarity simulations, which improves the efficiency and accuracy of system data acquisition, providing accurate acquisition data for experimental data analysis.

In this study, a physical similarity simulation plays a significant role in the study of crack evolution and the gas migration mechanism. A sensor is deployed inside a comparable artificial rock formation to assure the accuracy of the experiment results. During the building of the simulated rock formation, a huge volume of acidic gas is released, causing numerous sensor measurement mistakes. Additionally, the gas concentration estimation approach is subject to uncertainty because of the complex rock formation environment. As a result, the purpose of this study is to introduce an adaptive Kalman filter approach to reduce observation noise, increase the accuracy of the gas concentration estimation model and, finally, determine the gas migration law.

First, based on the process of gas floatation-diffusion and seepage, the gas migration model is established according to Fick’s second law, and a simplified modeling method using diffusion flux instead of gas concentration is presented. Second, an adaptive Kalman filter algorithm is introduced to establish a gas concentration estimation model, taking into account the model uncertainty and the unknown measurement noise. Finally, according to a large-scale physical similarity simulation platform, a thorough experiment about gas migration is carried out to extract gas concentration variation data with certain ventilation techniques and to create a gas chart of the time-changing trend.

This approach is used to determine the changing process of gas distribution for a certain ventilation mode. The results match the rock fissure distribution condition derived from the microseismic monitoring data, proving the effectiveness of the approach.

For the first time in large-scale three-dimensional physical similarity simulations, the adaptive Kalman filter data processing method based on the inverse Wishart probability density function is used to solve the problem of an inaccurate process and measurement noise, laying the groundwork for studying the gas migration law and determining the gas migration mechanism.

The unstable dynamic propagation of multistage hydrofracturing fractures leads to uneven development of the fracture network and research on the mechanism controlling this phenomenon indicates that the stress shadow effects around the fractures are the main mechanism causing this behaviour. Further studies and simulations of the stress shadow effects are necessary to understand the controlling mechanism and evaluate the fracturing effect.

In the process of stress-dependent unstable dynamic propagation of fractures, there are both continuous stress fields and discontinuous fractures; therefore, in order to study the stress-dependent unstable dynamic propagation of multistage fracture networks, a series of continuum-discontinuum numerical methods and models are reviewed, including the well-developed extended finite element method, displacement discontinuity method, boundary element method and finite element-discrete element method.

The superposition of the surrounding stress field during fracture propagation causes different degrees of stress shadow effects between fractures and the main controlling factors of stress shadow effects are fracture initiation sequence, perforation cluster spacing and well spacing. The perforation cluster spacing varies with the initiation sequence, resulting in different stress shadow effects between fractures; for example, the smaller the perforation cluster spacing and well spacing are, the stronger the stress shadow effects are and the more seriously the fracture propagation inhibition arises. Moreover, as the spacing of perforation clusters and well spacing increases, the stress shadow effects decrease and the fracture propagation follows an almost straight pattern. In addition, the computed results of the dynamic distribution of stress-dependent unstable dynamic propagation of fractures under different stress fields are summarised.

A state-of-art review of stress shadow effects and continuum-discontinuum methods for stress-dependent unstable dynamic propagation of multiple hydraulic fractures are well summarized and analysed. This paper can provide a reference for those engaged in the research of unstable dynamic propagation of multiple hydraulic structures and have a comprehensive grasp of the research in this field.

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.

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.

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.