Complex Geological Structures Research Articles

Correlation between Soil Moisture Change and Geological Disasters in E’bian Area (Sichuan, China)

E’bian Yi Autonomous County is a mineral-rich area located in a complex geological structure zone. The region experiences frequent geological disasters due to concentrated rainfall, steep terrain, and uneven vegetation cover. In particular, during the rainy season, large amounts of rainwater rapidly accumulate, increasing soil moisture and slope pressure, making landslides and debris flows more likely. Additionally, human activities such as mining, road construction, and building can alter the original geological structure, exacerbating the risk of geological disasters. According to publicly available data from the Leshan government, various types of geological disasters occurred in 2019, 2020, 2022, and 2023, resulting in economic losses and casualties. Although some studies have focused on geological disaster issues in E’bian, these studies are often limited to specific areas or types of disasters and lack comprehensive spatial and temporal analysis. Furthermore, due to constraints in technology, funding, and manpower, geophysical exploration, field geological exploration, and environmental ecological investigations have been challenging to carry out comprehensively, leading to insufficient and unsystematic data collection. To provide data support and monitoring for regional territorial spatial planning and geological disaster prevention and control, this paper proposes a new method to study the correlation between soil moisture changes and geological disasters. Six high-resolution Landsat remote sensing images were used as the main data sources to process the image band data, and terrain factors were extracted and classified using a digital elevation model (DEM). Meanwhile, a Normalized Difference Vegetation Index–Land Surface Temperature (NDVI-LST) feature space was constructed. The Temperature Vegetation Drought Index (TVDI) was calculated to analyze the variation trend and influencing factors of soil moisture in the study area. The research results showed that the variation in soil moisture in the study area was relatively stable, and the overall soil moisture content was high (0.18 < TVDI < 0.33). However, due to the large variation in topographic relief, it could provide power and be a source basis for geological disasters such as landslide and collapse, so the inversion value of TVDI was small. The minimum and maximum values of the correlation coefficient (R2) were 0.60 and 0.72, respectively, indicating that the surface water content was relatively large, which was in good agreement with the calculated results of vegetation coverage and conducive to the restoration of ecological stability. In general, based on the characteristics of remote sensing technology and the division of soil moisture critical values, the promoting and hindering effects of soil moisture on geological hazards can be accurately described, and the research results can provide effective guidance for the prevention and control of geological hazards in this region.

Seismic imaging of the complex geological structures in the southwestern edge of the Western limb, Bushveld Complex through focusing pre‐stack depth migration of legacy 2D seismic data

Seismic imaging of the complex geological structures in the southwestern edge of the Western limb, Bushveld Complex through focusing pre‐stack depth migration of legacy 2D seismic data

Impact fracturing of rock-like material using carbon dioxide under different temperatures and pressures

Unconventional resources (oil, gas, and geothermal) are often buried deep underground within dense rock strata and complex geological structures, making it increasingly difficult to create volumetric fractures through conventional hydraulic fracturing. This paper introduces a novel method of supercritical energetic fluid thermal shock fracturing. It pioneers a CO2 deflagration impact triaxial pneumatic fracturing experimental system, using high-strength similar materials to simulate deep, hard rock masses. The study investigates the rock-breaking process and crack propagation patterns under supercritical CO2 thermal shock, revealing and discussing the types of thermal shock-induced fractures, their formation conditions, and discrimination criteria. The research indicates that higher supercritical CO2 thermal shock pressures and faster pressure release rates facilitate the formation of radial branching fractures, circumferential cracks, and branch cracks. Typically, CO2 thermal shock generates 3–5 radial main cracks, which is significantly more than the single main crack formed by hydraulic fracturing. The formation of branched cracks is often caused by compression-shear failure and occurs under relatively harsh conditions, determined by the confining pressure, rock properties, peak thermal shock pressure, and the pressure sustained post-decompression. The findings are expected to offer a safe, efficient, and controllable shockwave method of supercritical fluid thermal shock fracturing for the exploitation of deep unconventional oil and gas resources.

UPDATE OF PETROPHYSICAL AND GEOLOGICAL MODEL OF MIDDLE JURASSIC DEPOSITS DURING WELL MONITORING

Tyumen deposits of the studied subsoil area of Western Siberia were put into operation from the middle of the last century, but active development of the YuS2 formation began not so long ago. Currently, the share of wells planned for the operation of this target in the coming years is growing. There are cases when the well performance does not meet the forecasts, which calls into question the correctness of the model used in the development of the target and the structure of the section itself.For high-quality planning of the production stock, it is necessary to prepare a reliable geological basis. The problem of updating the model and improving its reliability lies in the complex geological structure of the formation, which is characterized by high lateral and vertical heterogeneity, as well as in the absence of the basis for all geological models - seismic studies of the investigated subsoil area.The cross-section of the YuS2 reservoirs is complex for petrophysical and geological modeling and is characterized by low reservoir properties. The reservoirs are represented by thin interbeds of permeable interlayers and impermeable clay, carbonaceous and dense rocks.The article examines the specifics of building petrophysical models of YuS2 formations during development drilling using all available core and borehole data to detail the deposits of the Tyumen formation. The purpose of the work is to create own dependencies for interpretation of well geophysical survey data taking into account sedimentation conditions.To detail the Jurassic section, a comprehensive analysis of lithological and petrographic studies of core material from the point of view of geology and petrophysics was used.Facies analysis of the Tyumen deposits in the absence of seismic studies 3D acts as the main factor improving the quality of the predictive ability of geological models in support of development drilling. Adaptation of the petrophysical model to modern geological and technical conditions of well drilling is shown.An integrated approach to the study of the Tyumen formation deposits based on facies and petrophysical studies made it possible to improve the quality of the forecast of the distribution of the reservoir of formation YuS2 and obtain a detailed model. This will help optimize the placement of the production stock, as well as identify the prospects for the development of deposits in areas not covered by drilling.

Geotectonic architecture beneath Northern Vietnam revealed by local earthquake tomography combining seismic data from multiple networks

Extrusion tectonics driven by the collision of Indian-Eurasian continents and the relative motion among microcontinents (e.g. South China, Simao-Indochina blocks) has led to multi-phase and complex geological activity and structures in Northern Vietnam. While extensively studied, its detailed crustal architecture in association with multiple large-scale shear zones and plate boundaries remains unclear and is key to resolving long-lasting tectonic debates. This study constructs the first 3-D regional S-wave velocity (Vs) model with consistent P-wave velocity (Vp) and Vp/Vs ratio models by jointly using the crustal and head waves (e.g. Pg/Sg and Pn/Sn phases) from an augmented earthquake catalog combining the Vietnam National Earthquake Catalog (VNCAT) and the data from multiple adjacent seismic networks. The quality of phase picking and earthquake location incorporating seismic data from various sources are carefully examined and collected for the tomographic inversion. We use a stepwise joint inversion scheme to subsequently constrain the earthquake location, 1-D and 3-D velocity models. The output models suggest that the Song Ma suture represents a better candidate for the Indochina-South China boundary rather than the Red River Shear Zone. The latter seems to act as a crustal structure while the former and the Dien Bien Phu Fault are likely lithospheric ones. Several other geologic features such as the Tu Le volcanic basin, the Song Hong basin, and the Day Nui Con Voi metamorphic complex are also clearly depicted and discussed in the model. This new set of the Vp, Vs, and Vp/Vs models could further contribute to better assessment to seismic hazards in the northern Vietnam.

Discussion on Deep Geothermal Characteristics and Exploration Prospects in the Northern Jiangsu Basin

The Northern Jiangsu Basin (NJB), located at the northeast edge of the Yangtze block, is not only rich in oil and gas resources but also contains abundant geothermal resources. Nevertheless, the distribution of geothermal resources at medium depth in the NJB is still unclear due to its complex geological structure and tectonic–thermal evolution process, which restricts its exploitation and utilization. The characteristics of the geothermal field and distribution of geothermal reservoirs within the NJB are preliminarily analyzed based on available temperature measurements and geothermal exploration data. The prospective areas for the exploration of deep geothermal resources are discussed. The analysis results show that (1) Mesozoic–Paleozoic marine carbonate rocks are appropriate for use as principal geothermal reservoirs for the deep geothermal exploration and development within the NJB; (2) the geothermal field is evidently affected by the base fluctuation, and the high-temperature area is mainly concentrated at the junction of the Jianhu uplift and Dongtai depression; (3) the southeast margin of Jinhu sag, Lianbei sag, the east and west slope zone of Gaoyou sag, the low subuplifts within the depression such as Lingtangqiao–Liubao–Zheduo subuplifts, Xiaohai–Yuhua subuplifts and the west of Wubao low subuplift, have good prospects for deep geothermal exploration.

Изучение дизъюнктивной сети о. Кунашир (Курильские острова) с целью реконструкции особенностей его тектонического развития

The paper is devoted to the study of the faults on Kunashir Island, focusing on their classification, morphology, kinematics, timing, and causes of formation, which are key to understanding the tectonic evolution of the area. The magma-supporting role of these faults is noted, contributing to the understanding of the distribution of volcanic and plutonic formations in the Great Kuril Ridge. The ore-controlling role of the faults was evaluated, which may contribute to the discovery of new mineral deposits. Additionally, studying these faults is important for improving the prediction of earthquakes and tsunamis. The research methodology includes an integrated approach that combines traditional geological observations with modern remote sensing methods, allowing us to detail the map of the Kunashir Island rupture faults. The results of the work confirm the presence of three systems of tectonic structures that differ in depth of emplacement, length, timing, kinematics, and causes of formation. The study contributes to our understanding of the complex geological structure of Kunashir Island and the entire island arc and provides a basis for further research in this area."

Subsurface Imaging by a Post-Stimulation Walkaway Vertical Seismic Profile Using Distributed Acoustic Sensing at the Utah FORGE Enhanced Geothermal System Site

A 2D walkway vertical seismic profile (VSP) survey was conducted using a distributed acoustic sensing (DAS) system in southwest Utah, which is part of an enhanced geothermal system (EGS) project. The VSP was undertaken to obtain detailed structural information for a better understanding of the area’s subsurface geology and associated fracture development. By combining a 3D composite velocity model from previous studies and considering the complex geological structure beneath this region, we processed the data to create P-P depth image. We also modified the interval Q calculation using a moving window over the gauge-length corrected DAS record to generate the velocity profile and the comparable interval attenuation curve. The correlated P-P images from two DAS records successfully indicate not only the main contact between shallow unconsolidated sediments and the metamorphic basement rocks at 2650 ft (807.72 m) but also several distinct reflections related to the geological contacts. The refined velocity profiles and the depth images can provide baseline results for further seismic modeling and time-lapse imaging.

IORSim: A Mathematical Workflow for Field-Scale Geochemistry Simulations in Porous Media

Reservoir modeling consists of two key components: the reproduction of the historical performance and the prediction of the future reservoir performance. Industry-standard reservoir simulators must run fast on enormous and possibly unstructured grids while yet guaranteeing a reasonable representation of physical and chemical processes. However, computational demands limit simulators in capturing involved physical and geochemical mechanisms, especially when chemical reactions interfere with reservoir flow. This paper presents a mathematical workflow, implemented in IORSim, that makes it possible to add geochemical calculations to porous media flow simulators without access to the source code of the original host simulator. An industry-standard reservoir simulator calculates velocity fields of the fluid phases (e.g., water, oil, and gas), while IORSim calculates the transport and reaction of geochemical components. Depending on the simulation mode, the geochemical solver estimates updated relative and/or capillary pressure curves to modify the global fluid flow. As one of the key innovations of the coupling mechanism, IORSim uses a sorting algorithm to permute the grid cells along flow directions. Instead of solving an over-dimensionalized global matrix calling a Newton–Raphson solver, the geochemical software tool treats the species balance as a set of local nonlinear problems. Moreover, IORSim applies basis swapping and splay tree techniques to accelerate geochemical computations in complex full-field reservoir models. The presented work introduces the mathematical IORSim concept, verifies the chemical species advection, and demonstrates the IORSim computation efficiency. After validating the geochemical solver against reference software, IORSim is used to investigate the impact of seawater injection on the NCS Ekofisk reservoir chemistry.Article HighlightsThe IORSim sorting algorithm decouples the nonlinear geochemical reaction calculations into recurring one-dimensional problems to assure numerical stability and computation efficiency. To the best of our knowledge, this work presents the mathematical concept, implementation, and application of topological sorting for the first time on (industry) field-scale problems.IORSim combines topological sorting with basis swapping and splay trees to significantly reduce computation times. Moreover, a high-speed forward simulation mode was developed to allow the post-advection of chemical components to visualize species distribution, water chemistry, and mineral interactions. If the geochemical reactions interfere with the fluid flow, the IORSim backward mode uses relative permeability curves to update the global fluid flow at each time step.We validate the implemented topological scheme on a reservoir grid, show the computation efficiency, and compare the impact of explicit, implicit, and grid refinement on numerical dispersion.The decoupled flow simulator and geochemical reaction calculations allow seamless integration of full-field reservoir models that contain complex geological structures, a large number of wells, and long production histories.The computation capabilities of IORSim are demonstrated by simulating and reproducing the impact of seawater injection in the southern segment of the giant Ekofisk field (more than 50 years of injection and production history). IORSim shows that seawater injection changed the Ekofisk mineralogy and impacted the produced water chemistry. In the investigated Ekofisk case, seawater promoted calcite dissolution and led to the precipitation of magnesite and anhydrite. Moreover, surface complexation modeling revealed that sulfate is adsorbed on the calcite surface.

A Scientific Investigation of the Shangfang Mountain Yunshui Cave in Beijing Based on LiDAR Technology

The Yunshui Cave in Shangfang Mountain, Beijing, is a famous high-altitude karst cave in northern China. As the third scientific survey of Yunshui Cave in history, this is the first time to use the latest LiDAR technology to carry out a related detection survey. Traditional cave measurement methods are limited by natural conditions and make it difficult to reach the destination. Traditional methods mainly rely on experience and obtain data with strong subjectivity, making it difficult to conduct quantitative research and obtain reproducible results in the current information era. Applying LiDAR technology to cave measurement can obtain comprehensive and accurate digital measurement results within the same survey time and reveal many richer and more accurate features of Yunshui Cave. The obtained digital measurement results can be used for 3D modeling as well as provide a large amount of accurate basic data and preliminary materials for subsequent geological, environmental, and archaeological investigation and analysis, as well as cultural and tourism resource development. The rapid geological survey of Shangfang Mountain Yunshui Cave using LiDAR technology shows that LiDAR cave geological survey technology can achieve real-time collection of centimeter-level accuracy and generate billions of points of cloud data, greatly improving survey efficiency and accuracy. At the same time, digital survey results can be obtained. Through modeling and GIS technology, all on-site survey details can be easily moved back to the laboratory for real-scene reproduction, network sharing, and dissemination. This study provides a foundation for future explorations of the Yunshui cave and highlights the potential for LiDAR techniques to enhance our understanding of complex geological structures such as caves.

Assessing the application of ultradeep azimuthal electromagnetic resistivity tools in detecting igneous formations ahead of the transmitter

The Faroe–Shetland basin (FSB) is considered challenging for oil and gas exploration due to its complex geological structure and limited exploration drilling activity. The costs of drilling in the basin are the most expensive within the United Kingdom Continental Shelf (UKCS) due to varying geological and engineering challenges. The prevalence of igneous sill complexes within the basin present drilling hazards that contribute to drilling nonproductive time (NPT) and are documented in multiple studies. Recent advancements in ultradeep azimuthal resistivity (UDAR) technology which utilizes the look-ahead portion of the electromagnetic (EM) signal presents an opportunity to assess its suitability as a derisking tool to ‘look-ahead’ and detect igneous intrusions in volcanic basins. This paper utilizes resistivity log data from Well 214/28-1, which encountered multiple igneous intrusions and is used to validate the ability of these tools to detect the intrusions. The EM look-ahead technology detected with good confidence the top of an igneous intrusion 72 ft (22 m) true vertical depth (TVD) ahead of the transmitter and once drilled, detect the base of the intrusion 54 ft (16 m) TVD ahead of transmitter. The detection of the intrusions prior to drilling and information on intrusion thickness can enable proactive drilling decisions to reduce nonproductive time.

Hydrogeomorphology and recharge mechanism of a coastal aquifer in tectonic-controlled watersheds: A multi-proxy approach based on remote sensing and environmental isotopes

Coastal aquifers worldwide can be considered an essential water source required for preservation of the coastal ecosystems. However, these aquifers are vulnerable to seawater intrusions and over-pumping due to their proximity to the sea and human activities, respectively. Therefore, the investigation of recharge mechanism has special importance in the regions of tectonic rift activities where they exhibit complex geological structures and hydrogeological characteristics. In the current research, an integrated approach of remote sensing and environmental isotopes together with field investigation and subsurface datasets (aeromagnetic, gravity and well logs) are employed to identify the nature and factors affecting groundwater recharge of the Miocene coastal aquifer along the Red Sea-Gulf of Suez western margin (continental rift basin). The research findings reveal that: (1) The coastal Miocene aquifer consists of both clastic (sandstone, sand and gravels) and carbonate rocks (limestone and dolomite) with subsurface thickness ranges between 100 and 200 m. (2) Its groundwater shows contrast salinity values (expressed by total dissolved solids, TDS) between 2755 and 10,996 mg/l, due to the variation in recharge rates and the lithologic dissimilarity of the water bearing formations. (3) The environmental isotopes indicate that the Miocene groundwater has a mixed isotopic signature between modern meteoric (rainfall) and fossil waters. (4) The isotopic data of the Miocene aquifer and the enriched Gulf of Suez verifies that no seawater intrusions are affecting this aquifer with existing hydraulic barriers (clogged/sealed faults or impermeable massive blocks). (5) The hydrogeomorphological investigations, aeromagnetic and gravity data reveal an existence of two morphotectonic depressions (water collectors) which have paleo and recent recharge opportunities for the subsurface sedimentary layers (2–4 km thick). From the applied viewpoint, these two depressions have potential prospects for future groundwater exploration, which has significant impact on food security and land reclamation.

Stability Challenges and Remedial Practices in Himalayan Hydro-power Tunnels – A Review

The instability in tunnels is mainly affected by geological anomalies, rock mass quality, complex geological structures, active tectonics, and stress anisotropy. This review article presents challenges associated with stability and applied remedial measures prevailing in hydropower tunnels in the Himalayas. The review covers nine hydropower tunnels located in different parts of the Himalayas. The review found that rock bursting/spalling frequently occurs when the tunnel passes through a high overburden with good rock mass quality. On the other hand, plastic deformation (squeezing) occurs when a tunnel passes through the weak and schistose rock mass. It has been found through the review that the tunnel crew was able to successfully solve instability challenges. Effective planning, design, and selection of appropriate construction techniques help to complete tunneling projects in the Himalayas.

A Grand Challenge: Update on Improved Recovery From Tight/Shale Reservoirs

_ This is the third of a series of six articles on SPE’s Grand Challenges in Energy, formulated as the output of a 2023 workshop held by the SPE Research and Development Technical Section in Austin, Texas. Described in a JPT article last year, each of the challenges will be discussed separately in this series: geothermal energy; net-zero operations; improving recovery from tight/shale resources; carbon capture, storage, and utilization; digital transformation; and education and advocacy. _ The exploration and extraction of hydrocarbons from tight or shale formations have revolutionized the global energy landscape, unlocking vast oil and gas reserves previously considered inaccessible. According to estimates, substantial oil production has been reported for various tight plays in the US in 2024 (Fig. 1). However, the production process in tight or shale formations is not without its own set of challenges ranging from technological hurdles to environmental concerns. Therefore, understanding the current outlook is essential for stakeholders to navigate the complex terrain. To address the significant issues outlined by SPE in July 2023 (Halsey et al. 2023), this article intends to explore the complexities of hydrocarbon extraction from tight and shale reservoirs. Description of Tight/Shale Formations First, we present a succinct description of tight or shale formations. Tight or shale oil and gas are hydrocarbons found in oil-bearing mudstone, and shale gas is produced from gas shale or associated with tight oil (Boak and Kleinberg 2020). Tight unconventional formations are characterized by the presence of fine-grained sedimentary rocks that are high in organic matter, such as shale, where hydrocarbons are firmly embedded within the rock matrix. Challenges in Tight or Shale Oil Production The Existence of Complex Geometry. Tight and shale formations are known for their complex geological structures, which can make it difficult to predict the distribution and behavior of hydrocarbons within the reservoir. Additionally, shale formations can exhibit significant variability in their thickness, composition, and the occurrence of natural fractures. Hence, it is crucial to understand and navigate the intricate geology of unconventional shale formations to achieve successful shale oil production (Jiang et al. 2016). Ultralow Reservoir Permeability. The low permeability of tight formations restricts the flow of hydrocarbons, necessitating techniques such as hydraulic fracturing to create artificial pathways for extraction. However, the effective stimulation of fractures requires the use of cutting-edge drilling and extraction methods, such as multistage hydraulic fracturing (fracking) and horizontal drilling, which in turn require precise engineering and a thorough understanding of reservoir characteristics (Pokalai et al. 2015).

Space performance assessment of traditional houses in Akre

ABSTRACT The complex natural and geological structure of cities affects buildings’ spatial structure in terms of the forms of spatial configuration adopted in the design. Spatial structure affects buildings’ space performance in several aspects, such as privacy and accessibility, which are the most significant indicators of residential houses’ functional performance. Akre, a treasured, historical city in Iraqi Kurdistan, is characterised by a complex and distinct geological nature, and its residential part contains numerous traditional houses that deserve analysis. This study employed both a qualitative approach – via case study plan preparations and photographic documentation – and a quantitative approach, utilising space syntax tools. The results were analysed with the statistical program SPSS and show three types of spatial configuration in the residential part (linear, central and cluster); each type has a different and effective impact on the achievement of space performance in terms of privacy and accessibility. The study concludes that the linear type represents the most private typology, followed by the cluster and then the central type. The centralised typology is the most connected and accessible to different spaces within the residential unit, followed by the cluster type and the linear type, with a significant difference.

Elastic reverse time migration based on nested triangular mesh in 2D isotropic media

Abstract Reverse time migration (RTM) based on two-way wave equation is an important imaging tool, which has been widely used in exploration geophysics. It offers distinct advantages in handling complex geological structures. Moreover, elastic reverse time migration (ERTM) can produce more physically meaningful images for subsurface and has received increasing attention in recent years. However, rugged topography presents challenges in the reconstruction of subsurface structures for ERTM. To overcome this problem, we introduce the grid method for topography ERTM, which is an unstructured mesh based algorithm. However, the conventional gird method for RTM is implementing by discretizing model into small meshes, leading to high computational cost. To address this problem, we introduce the grid method based on nested unstructured meshes to implement ERTM with complex surface topography. Low-order small meshes, due to their fine scale, can accurately depict complex surface topography. However, their generation process is relatively time-consuming. In contrast, high-order large meshes contain many similar small meshes internally, which can be processed in parallel, resulting in higher computational efficiency. Therefore, we propose the following strategy: for the receiver wavefield, several layers on surface discreted by the small meshes are used to place receivers and remaining layers discreted by the large meshes to improve computational efficiency. For the source wavefield, only the large meshes are used because the spacing of shot is relatively large. In this way, we ensure the accurate description of complex surface topography for the implementation of ERTM with high computational efficiency. Finally, the effectiveness of our method is validated through experiments on three models with complex topography.

DC3DPAFEM: An efficient and accurate 3-D direct current resistivity anisotropic forward modeling software for complex geological settings

Nowadays, there is a growing trend that direct current (DC) field surveys are shifting towards challenging areas characterized by mountainous topography and electrical anisotropy. Given these complex geological settings, there is an urgent need for 3-D DC forward modeling software capable of effectively addressing large-scale problems and delivering accurate modeling results to interpret field data. However, most open-source software packages face certain limitations, such as the high numerical cost to handle complex surface topography, the lack of consideration for anisotropic conductivity, the absence of mesh refinement techniques to guarantee accuracy in forward modeling, and the lack of parallel computing techniques to solve large-scale problems. In this study, we develop an efficient and highly accurate 3-D DC anisotropic forward modeling software, namely DC3DPAFEM, using the adaptive finite element algorithm based on the unstructured tetrahedral mesh. Firstly, we construct a strong compatible boundary value problem (BVP) for 3-D anisotropic DC problems by adopting a specialized secondary potential approach to handle the surface topography efficiently. Then, we develop a goal-oriented adaptive mesh refinement (AMR) technique to ensure accurate forward modeling results, even with a coarse initial mesh. To ensure time and memory efficiency, we employ a robust conjugate gradient (CG) algorithm preconditioned by the algebraic multigrid (AMG) solver to solve the large-scale linear system of equations resulting from complex geological structures. We aim to investigate the performance of the AMG scheme in anisotropic DC cases. Furthermore, we incorporate the domain decomposition technique into the iterative solution scheme for further efficiency gains. This technique significantly improves computing efficiency for large-scale problems in parallel clusters. Finally, we conduct comprehensive performance tests for DC3DPAFEM using a two-layer anisotropic model and a 3-D complex model with undulating terrain. The results of both examples validate the accuracy of DC3DPAFEM, as they closely align with the analytical solutions and the solutions obtained from the existing 3-D DC forward modeling code. Compared to traditional direct solver MUMPS and ILU-preconditioned iterative solvers, DC3DPAFEM exhibits highly scalable performance for large-scale problems, offering significant advantages in terms of memory and time consumption. Overall, DC3DPAFEM demonstrates substantial advances in efficiency, accuracy, and practicality through a series of numerical examples. This open-source code provides an efficient and available tool for developing a 3-D DC inversion method that can deal with large-scale problems involving intricate topography and anisotropic media.

Dynamic analysis of field-scale rockslides based on three-dimensional discontinuous smoothed particle hydrodynamics: A case study of Tangjiashan rockslide

To improve the understanding of the dynamic disastrous process of field-scale rockslides, a novel numerical approach, Three-dimensional Discontinuous Smoothed Particle Hydrodynamics (3DSPH), was originally developed. This method comprehensively captures sequential stages of crack initiation and propagation, formation of contacts, frictional slip, catastrophic slides of rock masses, and final deposition, which was verified by three benchmark tests including the bouncing test of rigid balls, block sliding test and unconfined compression test of layered rock specimens. The approach was subsequently employed in the case of the Tangjiashan rockslide blocking valley event with particular attention paid to the influence of layered rock structure on the sliding and deposition processes. Field survey, geomorphological analysis and laboratory test of rock specimens were conducted to determine fundamental geological conditions and parameters required by the numerical simulation. Finally, 3D rockslide simulations with different rock layer thickness and strength subject to seismicity were conducted. The duration of the actual Tangjiashan rockslide's valley-blocking event (approximately 60 s) and the deposition area derived from the numerical simulation closely align with the field investigation. The rockslide mode is characterized by an ‘en masse’ motion with a peak sliding velocity of approximately 35–37 m/s. This single numerical code systematically elucidated the authentic attributes of the sliding process of large-scale rockslides, and realistically captured the characteristics of preservation of layered features within the sliding mass and ‘high-speed and short-distance’ movements with fluidization. These insights offer a fresh perspective for understanding the dynamics of large-scale rockslides with complex geological structures and subsequent accumulation processes.

Study of hydraulic cutting as a method to prevent coal burst disasters

This study explores the causes of coal bursts in the Xinzhou Kiln Mine, identifying key factors such as residual pillars, hard coal seams and/or roofs, stress concentration due to complex geological structures, and the stress distribution characteristics of the primary rock. A significant finding is that hydraulic cutting not only diminishes and redistributes the stress concentration region inside the coal seam but also mitigates the burst potential of the coal-rock mass, fundamentally reducing the likelihood of coal bursts. By taking Face No. 8937 in Xinzhou Kiln Mine as the test object, a coal burst prevention test was performed using hydraulic cutting. In combination with theoretical analysis and numerical simulation, the mechanism of hydraulic cutting for preventing coal burst was discussed, and reasonable cutting parameters were established. Onsite monitoring revealed that hydraulic cutting disrupts the integrity of the coal-rock mass, releases internal stress, and increases its water content, thereby weakening its burst tendency. Additionally, the deformation and fracturing of the cutting slots and the closure of boreholes shifted the stress concentration from the coal seam to deeper areas and to the two ribs. Post-cutting observations showed a significant reduction in both the frequency and impact energy of coal bursts; there was also a noticeable increase in the convergence of the roadway in the cutting area compared to non-cutting areas. Furthermore, displacement of the roof and floor increased by 78.9 % and that of the two ribs increased by 47.4 % after cutting, preventing the coal-rock mass from accumulating high stress. In conclusion, hydraulic cutting is a promising method for effectively preventing coal bursts and enhancing the safety of mining operations.

Stresses Induced by Magma Chamber Pressurization Altered by Mechanical Layering and Layer Dip

AbstractUnderstanding the stress distribution around shallow magma chambers is vital for forecasting eruption sites and magma propagation directions. To achieve accurate forecasts, comprehensive insight into the stress field surrounding magma chambers and near the surface is essential. Existing stress models for pressurized magma chambers often assume a homogenous elastic half‐space or a heterogeneous crust with varying mechanical properties in horizontal layers. However, as many volcanoes have complex, non‐horizontal, and heterogeneous layers, we enhance these assumptions by considering mechanically stratified layers with varying dips. We employed the Finite Element Method (FEM) to create numerical models simulating three chamber geometries: circular, sill‐like and prolate. The primary condition was a 10 MPa excess pressure within the magma chamber, generating the stress field. Layers dips by 20‐degree increments, with differing elastic moduli, represented by stiffness ratios of the successive layers (EU/EL) ranging from 0.01 to 100. Our findings validate prior research on heterogeneous crustal modeling, showing that high stiffness ratios disrupt stress within layers and induce local stress rotations at mismatched interfaces. Layer dip further influences stress fields, shifting the location of maximum stress concentration over varying distances. This study underscores the significance of accurately understanding mechanical properties, layer dip in volcanoes, and magma chamber geometry. Improving forecasting of future eruption vents in active volcanoes, particularly in the Andes with its deformed, folded, and non‐horizontal stratified crust, hinges on this knowledge. By expanding stress models to incorporate complex geological structures, we enhance our ability to forecast eruption sites and magma propagation paths.