Geomechanics and Geophysics for Geo-Energy and Geo-Resources最新文献

Studying crystal anisotropy is of great importance for understanding the thermal–mechanical coupling behavior of crystalline rocks in deep underground engineering. In this study, a microscopic parameter calibration method incorporating the size effect is proposed. Subsequently, a thermal–mechanical coupling model accounting for the quartz crystal anisotropy is established to investigate the thermal–mechanical coupling behavior of monomineral quartzite. The results show that thermal-induced microcracks are exclusively distributed along crystal boundaries, and initiate preferentially from crystal boundaries with a larger average linear thermal expansion coefficient, eventually leading to the formation of a crack network. With the increase in temperature, the peak strength of monomineral quartzite increases slightly at first and then decreases rapidly, and the transition threshold temperature is 200 °C. Both elastic modulus and Poisson’s ratio show a monotonic pattern, with abrupt changes occurring at 200 and 300 °C, respectively. The monomineral quartzite exhibits a significant compaction stage under uniaxial compression, and the ductile strengthening critical temperature for monomineral quartzite are between 400 and 500 °C. The quartz crystal anisotropy leads to an anisotropic distribution of inclination angles for tensile microcracks under high temperatures while having no obvious effect on the shear microcracks. In addition, the average size of fragments generated under uniaxial compression is influenced by thermal cracking, demonstrating an initial decrease followed by an increase, and the distribution of fragment sizes is solely correlated with the temperature, which is more concentrated with the increase in temperature.

This study presents a comprehensive procedure for determining, calibrating, and validating the Riedel–Hiermaier–Thoma (RHT) material model parameters for granite. The process involves collecting a comprehensive dataset of conventional mechanical tests conducted on various types of granite worldwide. Based on this dataset, one set of RHT material model parameters is determined. Additionally, a specific granite sample from Iran, known as Nehbandan granite, is characterized through physical and mechanical testing to obtain another set of parameters. The challenging task of determining the third set of parameters, which are difficult to obtain analytically and experimentally, is accomplished through a calibration process that iteratively adjusts the parameters based on comparisons between numerical simulation results and experimental data. To validate the determined parameters, a series of tests, including uniaxial compressive strength (UCS), triaxial compressive strength (TCS), Brazilian tensile strength, and dynamic Brazilian using split Hopkinson pressure bar (SHPB) tests, are conducted on the Nehbandan granite. These tests are also simulated using LS-Dyna software, and the numerical simulation results are compared with the corresponding experimental test results. The comparison between the numerical and experimental data serves as a means of validating and verifying the accuracy and reliability of the determined RHT material model parameters for granite. The results demonstrate the successful determination and calibration of the RHT material model parameters for granite. The model exhibits effectiveness in predicting the behavior of granite under various loading conditions. The validation process confirms the accuracy and reliability of the determined parameters through a close agreement between numerical simulations and experimental data. The findings contribute to a better understanding of granite's mechanical response and provide a reliable tool for simulating and predicting its behavior in engineering applications. The validated RHT material model parameters offer a robust framework for accurate numerical simulations, enabling engineers to make informed decisions in rock engineering projects.

Forming a fracture network through fracturing stimulation is significant to efficiently developing shale oil resources. However, the complex lithological characteristics and dense laminas of continental shale oil strongly shield fracture propagation. The concept of "cyclic fluid injection induces rock fatigue" was introduced into shale oil fracturing technology, and the cyclic pressure shock fracturing method was designed. The horizontal well fracturing simulation experiments used prepared shale rock samples from the Permian Lucaogou Formation shale oil reservoir outcrop in Jimusar Sag, Junggar Basin. Two pressurization states were obtained through constant injection and rapid release of accumulated high pressure, corresponding to conventional and pressure shock conditions. The characteristics of fracture propagation under different fracturing methods were analyzed by combining acoustic emission monitoring and injection pressure curve response. Research has found that the dense laminas with a certain original width near the wellbore significantly inhibit the vertical propagation of hydraulic fractures (HFs), and conventional constant-rate fracturing methods make it difficult to stimulate the reservoir effectively. Fatigue fracturing can increase the complexity of near-wellbore fractures, but the HFs still tend to be arrested by the laminas. The bottom hole pressure (BHP) is artificially increased to a value far exceeding the rock breakdown pressure near the wellbore by applying the cyclic pressure shock fracturing method. It can avoid the communication between micro-cracks and horizontal laminas during the BHP constant rate increase process and overcome the inhibition of weak layers on vertical propagation. Besides, the fracture height and number of activated laminas positively correlate with the number of cycles. When the shock pressure is about 30 MPa, the fracture height of the three cycles increases by 50% compared to a single shock. In addition, the shock pressure has a more significant effect on the fracture height, and the fracture height increases significantly with the increase of shock pressure. The shock pressure increased by about 57%, and the fracture height increased by 60%. High-pressure shock has a certain effect on the naturally weak surface of the far well, which may cause the slip of the weak surface to produce AE signals. This provides a new approach for improving fracture complexity and control volume in layered shale oil reservoirs with dense laminas.

Owing to the rapid decline of oil production from tight oil reservoirs after primary hydraulic fracturing treatment on the horizontal well, a refracturing stimulation is proposed for tight oil recovery. In this study, a fully coupled dynamic stress computational method with a finite-element method is presented, and depletion-induced dynamic stress is simulated by coupled numerical modeling. In addition, the extended finite-element method (XFEM) approach is used to investigate the effect of different parameters on fracture dynamic propagation in recompletion from refracturing. The results highlight the effects of production-induced stress changes following hydraulic fracture propagation. When multiple fractures are stimulated simultaneously from recompletion during refracturing, curved fractures are commonly observed, and the deflected fractures generally divert toward the primary stimulated area with low pore pressure. The results indicate that comprehensive factors can affect the hydraulic fractures propagation from recompletion. The optimal refracturing time window can be determined using the dynamic stress condition and stimulated area. The initial completion spacing, initial fracture length, and recompletion perforation cluster spacing can also affect the fracture geometry from the recompletion. A larger initial fracture length can induce a larger stress change area, whereas a larger distance between the new perforations in recompletion and the old perforations can decrease the depletion-induced stress effect. A high horizontal stress contrast can increase the depletion-induced stress effect because a long fracture extends the area. Owing to the nonuniform pressure and stress distributions, more nonuniform fractures are commonly generated in the refracturing treatment. Thus, temporary plugging injection and proppant inertia must be designed while reducing the number of perforations near the initial perforation positions. This can help decrease the possibility of strong curved fractures and screen out problems.

Roughness ubiquitously prevails in rock joints and controls the shear behaviours, permeability and damage characteristics of rock joints. A plethora of investigations have focused on the description of joint roughness; however, a detailed method for quantifying joint roughness and evaluating the shear strength has not yet been established. In this study, within the framework of fractal theory, an optical measurement scale was defined to depict the fractal characteristics of joint roughness, and a boundary measurement scale was used to identify first- and second-order asperities. A composite indicator (eta), including the fractal roughness coefficient ((R_{D})), the coefficient describing the roughness order ((mu)) and the anisotropy parameter ((lambda)), was proposed to quantify the surface morphology, which takes the asperity distribution and roughness anisotropy into account. The relationship between (eta) and JRC was established, and the JRC–JCS criterion was further updated. Moreover, representative examples were given to show the steps required to quantify the morphology of rough joint surfaces with the new quantitative parameter. Direct shear tests were conducted to validate the effectiveness of the proposed method in describing joint roughness and estimating joint shear strength; the results indicate that it is appropriate to use (eta) to estimate the joint roughness and that the proposed shear strength criterion can appropriately predict the shear strength within an acceptable error.

In this study, a discriminative model of fracture expansion through a composite reservoir interface was established. During fracture transformation, the model, based on the theory of energy release rate at the fracture tip, helps clarify the critical conditions and assess the law of hydraulic fracture expansion through coal composite reservoir interfaces. This model considers shear and tensile processes during fracturing at the coal–rock interface and the effects of interface material property differences on fracture extension. In addition, fracture extension pattern differentiation is achieved by comparing the hydraulic fracture tip penetration and bending energy release rates (G_p/G_d) with the ratio of reservoir fracture and interface fracture energies (Γ_R/Γ_F). The accuracy of this discriminant criterion was verified by comparing the predicted results of the discriminant model with fracturing test results under identical conditions. The effects of various factors on G_p/G_d and hydraulic fracture penetration extension were investigated by using the sand-mudstone composite reservoir of the Shibox Formation in the Linxing Block, China, as the geological background. The results indicate that for certain values of Γ_R/Γ_F, G_p/G_d increases logarithmically with an increase in water pressure in the wellbore. Moreover, the difference in ground stress and elastic modulus between the layers exponentially increases with increasing Poisson’s ratio difference and the angle between fracture and interface and decreases logarithmically with increasing fracture height. These results indicate that hydraulic fractures are likely to penetrate from reservoirs with a high elastic modulus and Poisson’s ratio into reservoirs with a low elastic modulus and Poisson’s ratio. Furthermore, hydraulic fractures are likely to penetrate formations with large differences in pinch angle and ground stress. Under specific geological conditions, high water pressure in the wellbore and small seam height are favourable for fracture penetration.

High-intensity coal mining has induced a series of ecological and environmental problems issues, including surface subsidence, the development of ground cracks, and the deterioration of vegetation. The disruption of water circulation systems induced by mining, such as perched groundwater, groundwater of aeration zone, and phreatic water, is the root cause of vegetation withering. The aeration zone serves as a crucial nexus in the process of water cycling and exerts a significant influence on soil fertility. To explore the characteristics of soil moisture transport in subsidence areas under the mining disturbance, on-site monitoring of the size and morphology characteristics of subsidence areas and ground cracks was conducted in typical mining areas in Inner Mongolia, China. Subsequently, a typical soil moisture transport model was constructed in subsidence areas, the soil moisture transport patterns under the influence of different types of subsidence and cracks were analyzed, and the influence law of soil damage on soil moisture transport in the aerated zone was clarified. The results indicate that (1) Based on the occurrence and distribution characteristics of subsidence cracks, the subsidence area can be divided into tension zone, compression zone, and neutral zone; the ground cracks are divided into permanent tension cracks and dynamic cracks. (2) The drought stress effect of soil in the subsidence area is significant. Under the influence of soil structure variation, the water-holding capacity of the soil in the subsidence area decreases, and the soil moisture dissipation is strong. The soil moisture transport rate in the aeration zone of the subsidence area is ranked as follows: tension zone > neutral zone > compression zone. (3) Ground cracks can exacerbate the soil moisture transport rate in the aeration zone. After 15 d of crack appearance, the soil moisture transport reaches a relatively stable state, and the soil moisture transport rate in the surface layer of the crack is the fastest, and the loss of soil moisture is the most significant. The crack effect is not significant beyond 100 cm from the crack. This study provides a theoretical and data support for soil and vegetation remediation in mining subsidence areas.

Lately, Saudi Arabia has been developing its mineral exploration. However, comprehensive studies of the collected data are not accessible. Thus, the purpose of this research is to identify and map the hydrothermal alteration zones and structural lineaments that regulate the mineral occurrences in the Nuqrah region of the Kingdom of Saudi Arabia using remote sensing and aeromagnetic data. To achieve the desired goal of the study, ASTER remote sensing data were employed, and they were processed in several ways, including principal component analysis, band ratio, and false color composites to reveal the zones of alteration and structure lineaments. In addition, aeromagnetic data was employed to map the lineaments controlling the mineralization. These datasets were integrated using GIS tools to generate a new mineralization potential map of the Nuqrah area, which was classified into three classes: low, moderate, and high probability mineralization. The results showed thirteen intriguing anomalies (high potential mineralization) dispersed over the research area to be prospected. Additionally, techniques such as residual, regional, first vertical derivative, and tilt derivative were utilized to detect the potential mineral-related geologic structures. The results were validated by plotting known mineralization sites on our maps. Six significant faulting trends have been found, according to the lineament map and rose diagrams from remote sensing: NE–SW, WNW–ESE, N–S, ENE–WSW, NNE–SSW, NW–SE, and E–W. The research region is most affected by the NW–SE, ENE–WSW, E–W, and N–S trends, which are organized in decreasing order of magnitude, according to the rose diagram of the aeromagnetic maps. The applied approach can be employed to map potential mineral deposits in Saudi Arabia and similar zones around the globe.

Using surfactants to extract oil, the anionic surfactant Karamay petroleum sulfonate (KPS), the zwitterionic surfactant octadecyl betaine (BS-18) and the nonionic surfactant coconut oil fatty acid diethanolamide (6501) were selected for adsorption experiments with minerals contained in the conglomerate reservoir with different mineral compositions to study the adsorption law of different types of surfactants. Zeolite and montmorillonite, which have the highest specific surface area and zeta potential among the minerals in the conglomerate reservoir, have the most obvious adsorption effect on surfactants, resulting in a large amount of adsorption of KPS and BS-18. The three types of surfactants were then used to conduct physical simulation oil recovery experiments with conglomerate core samples, and the results showed that 6501 had better overall performance, the best adsorption resistance, and a higher degree of recovery in oil recovery experiments, which provided a basis for the selection of surfactants in the process of chemical drive in conglomerate reservoirs.

Abandoned coal mine goaf is affected by air leakages and prone to spontaneous combustion, resulting in environmental pollution and geological disasters. Haizhou Open-pit Mine adopts both underground and open-pit mining methods. During the long-term mining process, the original stable stratum structure is constantly destroyed, and the slope slides, increasing cracks and severe air leakage around the goaf and roadway. The spontaneous combustion of coal is particularly prominent after the mine shut down. At present, there is no suitable indirect monitoring method to effectively explore the spontaneous combustion area in goaf. The study developed an all-weather monitoring plan and conducted multi-point continuous long-term measurements of the spontaneous combustion state in one abandoned coal mine goaf located in the eastern part of the Haizhou Open-pit Mine. We evaluated the dynamic correlation between surface CO₂ flux (SCF) and changes in the underground fire areas, determined the scope and evolution trend of the fire areas, and identified the distribution and change laws of SCF. The results show a significant positive correlation between SCF and soil temperature; moreover, the SCF value was found to reflect the CO₂ emission intensity of the goaf. The high SCF in the test area showed month-wise expansion and increase, while the CO₂ emission gradually increased monthly, and the calculated annual total emission was approximately 7017 t. Hence, the study can further provide guidance for the monitoring of spontaneous combustion in shallow coal seams, goaf and the assessment of CO₂ emissions from underground coal fires through the on-site monitoring and analysis results.