Rock Structure Research Articles

Sustainability-oriented construction materials for traditional residential buildings: From material characteristics to environmental suitability

The increasing replacement of traditional construction materials with modern alternatives often overlooks their superior environmental adaptability and sustainability, resulting in missed opportunities for energy efficiency and carbon emission reductions. This study aimed to address this issue by scientifically evaluating the material properties of volcanic rock from Hainan Island, China, to assess its potential for sustainable building practices in tropical climates. A series of material characterization techniques were employed to analyze the composition and structure of the volcanic rock, including optical profilometry (OP), polarized light microscopy (PLM), X-ray diffraction (XRD), X-ray fluorescence (XRF), scanning electron microscopy (SEM) with energy-dispersive spectroscopy (EDS), transient plane source technique (TPS), and mercury intrusion porosimetry (MIP). Our findings identified the rock as vesicular cryptocrystalline basalt, with its unique pore structure and mineral composition jointly contributing to a low thermal conductivity of 0.22 ± 0.04 W/(m·K). The pore structure was characterized by closed spherical cavities, needle-like lava fillings, nanoporous surfaces, complex pore-throat networks, and an extensive pore size distribution. These features significantly reduced gas convection and enhanced thermal resistance. Additionally, the rock composition, dominated by the amorphous volcanic glass and scattered hollow skeletal crystalline plagioclase, further diminished heat transfer. These characteristics indicated that Hainan volcanic rock possessed excellent thermal insulation properties, making it highly suitable for energy-efficient construction in hot climates. This study highlighted the potential for traditional materials like volcanic rock to contribute to sustainable building practices and the preservation of cultural heritage.

Synergistic application of digital outcrop characterization techniques and deep learning algorithms in geological exploration

In order to meet the needs of geologists for the analysis of data characterizing field outcrops (rock sections or formations exposed on the ground surface), this study developed a field digital outcrop visualization platform based on Cesium (a 3D geospatial visualization technology) digital outcrop characterization technology. The platform was developed based on WebGL (a protocol for rendering interactions on web pages), which overcame the shortcomings of traditional software in terms of visualization, cross-device, cross-platform, and ease of use. Firstly, UAV inclined photography is used for data collection, which transforms a large amount of geological data into an intuitive 3D geological model, while the visualization platform provides rich measurement and mapping tools for the identified features, which more intuitively displays the outcrop information, helps geological explorers to understand the geological conditions in the field more quickly and comprehensively, and improves the analysis efficiency and ease-of-use of outcrop characterization data. Combined with the improved VGG19 (a deep convolutional neural network architecture) algorithm model, it has excellent performance in dealing with the fine texture and complex structure of rocks, which significantly improves the accuracy of lithology identification. The synergistic application of this technology provides geologists with a faster and more comprehensive means to understand the geological conditions in the field. The reliability of combining the Cesium digital outcrop characterization technology with the VGG19 lithology identification algorithm in geological exploration is verified through case studies. The synergistic application of this technology will greatly enhance the efficiency and ease of analysis of outcrop characterization in the field, and provide new perspectives for future research in geosciences.

Localization of AE sources in rocks improved by enhanced arrival time localization

Accurate localization of Acoustic Emission (AE) events in rock structures is crucial for investigating rock failure mechanisms. This paper presents a comprehensive analysis of 3D AE localization in basalt-type rock samples, utilizing a novel combination of three powerful techniques: the Simulated Annealing Algorithm (SAA), an innovative AIC-EEMD picking method, and a velocity-free approach to wave propagation. We employed a localized Hsu–Nielsen source to generate elastic waves, which were captured by six AE sensors strategically placed on the rock surface. To address challenges in identifying P-wave arrivals, particularly in noisy or spurious data, we developed an innovative technique that integrates the Akaike Information Criterion (AIC) with Empirical Ensemble Mode Decomposition (EEMD) to enhance signal clarity. The performance of our proposal was validated achieving an average absolute distance error of 3.74 mm. Compared to traditional methods, our technique demonstrates superior accuracy and offers enhanced robustness in the presence of noise.

Study on influencing factors of controllable mechanical behavior of rock-like porous materials.

Due to the discrete and non-homogeneous of similar materials and the inability to realize large-size and original scale modeling, it is difficult to restore the structure and stress state of underground coal and rock mass in similar simulation tests. To solve this problem, a lightweight and suitable for large-scale modeling similar material, rock-like porous material has been developed. The quasi-static uniaxial compression experiment was carried out by using the large tonnage multi-module electronic control test system. And the influencing factors of controllable mechanical behavior of rock-like porous materials were studied. The results showed that, under uniaxial compression conditions, the material stress-strain curve exhibits three phases: elastic stage, failure stage, and platform stage. The uniaxial compressive strength, elasticity modulus, stress drop, and softening modulus of rock-like porous materials basically increase with the increase of density. The stress after peak strength changes from a slow decrease to a "stepped" or even "cliff like" downward trend. Polypropylene fibers have the effect of enhancing the uniaxial compressive strength, elasticity modulus, stress drop, softening modulus, shear deformation, and residual strength stability of rock-like porous materials. The rock-like porous material has a critical loading velocity, and it increases with density. At the critical loading velocity, the material shows obvious shear failure, and the shear inclination angle is the largest, and so is the uniaxial compressive strength. Through the experimental research, the influence laws of density, polypropylene fiber, and loading velocity on the failure mode, mechanical parameters, and mechanical behavior of the material are clarified, and the quantitative relationship between density and each mechanical parameter is obtained. The research is helpful to realize the accurate control of mechanical behavior of rock-like porous materials and further inverts the deformation and failure mechanism of underground coal and rock structures through indoor similar simulation tests.

Fault zone architecture and hydraulic properties characterization using earth tide analysis

Fault zone architecture and hydraulic properties, which determine subsurface fluid flow and storage, have received long-term research attention from hydrogeologists. Current field experimental techniques and numerical simulation methods are limited by the scale and complexity of the investigations and are inadequate for deep subsurface fault zones. Estimating the internal rock properties and structures of fracture zones requires a combination of methods from many different fields. In recent years, groundwater tidal response methods have been widely used to characterize subsurface hydrogeomechanical properties because of their advantages of being in situ and low cost. In this study, tidal analysis was employed for a well within the Blue Ridge Fault Zone, and we integrated the fracture zone properties revealed by a single well. Well W-03 revealed an aquifer with an overall permeability of ∼ 10–14 m2 and hydraulic diffusivity of 0.24 m2⋅s−1. The fault damage zone had a permeability of 10–14 ∼ 10–13 m2 while the host rock had a permeability of 10–16 ∼ 10–14 m2. In addition, we found that the heterogeneity of specific storage affects the assessment of poroelasticity parameters, we offer a method to estimate these parameters without relying on specific storage, yielding porosity estimates ranging from 0.05 to 0.06, aligning with field test results. Furthermore, the estimated dip and strike of the fault zone were also consistent with the outcrop and borehole observations.

Improving fluid identification in well logging using Continuous Wavelet Transform and Vision Transformers: An innovative approach

Well logging fluid prediction is one of the key steps in assessing oil and gas reserves. By analyzing downhole logging data, different types of fluids contained in underground rocks, such as crude oil, natural gas, and water, can be determined. This information is crucial for assessing the abundance and recoverable reserves of oil and gas resources and helps guide oil and gas exploration and development work. We have introduced a novel model called CWT (Continuous Wavelet Transform)-ViT (Vision Transformer). CWT can simultaneously provide frequency information at different scales, enabling the model to analyze downhole logging data more comprehensively and accurately at different scales. Underground rock structures often exhibit features at multiple scales, and CWT can effectively capture these features, aiding in better differentiation of different types of fluids. The ViT model utilizes the Transformer architecture, allowing for global attention over input sequences without being limited by sequence length. This enables the model to comprehensively understand the overall information of downhole logging data and extract richer features. For complex geological structures and fluid distributions in geological exploration, the global attention mechanism helps the model better grasp the overall situation, thereby improving the accuracy of fluid prediction. When we used the CWT-ViT method for well logging fluid prediction, we achieved a high accuracy rate of 97.50% in the first dataset, which further improved to 97.77% in the second dataset. These results demonstrate the significant robustness and efficiency of the CWT-ViT method in lithology prediction using well logging data. We also conducted blind well experiments, and our CWT-ViT model outperformed other models, achieving a blind well prediction accuracy of 97.36%. Therefore, the experiments indicate that the key to improving accuracy in well logging fluid prediction with CWT lies in its multiscale analysis capability, effectively capturing different fluid characteristic frequencies. Additionally, CWT enhances signal features and removes noise, increasing the precision of fluid identification. Finally, the integration with ViT further optimizes fluid prediction performance, making it outstanding in complex geological environments. The advantages of ViT in fluid prediction include its excellent sequence modeling capability, effective handling of long-distance dependencies, and enhanced ability to capture fluid characteristics in complex well logging data.

IDENTIFICATION OF AQUIFERS BASED ON THE VERTICAL ELECTRICAL SOUNDING (VES) METHOD SCHLUMBERGER CONFIGURATION CASE STUDY: PULAU BAAI KAMPUNG MELAYU SUB-DISTRICT, BENGKULU CITY, INDONESIA

Investigation of the groundwater potential in the Pulau Baai area, Kampung Melayu Sub-district, Bengkulu City, must be carried out in such a way that the activities and needs of the people in the area can be fulfilled and the needs of the population in the area can be met. This study aims to determine the status of groundwater using the Schlumberger configuration geoelectric method. Measurements were made using a resistivity meter, and the results for each configuration depended on changes in resistivity. Measurements for each configuration depend on changes in resistivity at depth, the vertical direction (sounding), and the lateral direction (mapping), so hydrogeological analysis in this activity aims to get the maximum use of groundwater / underground water in aquifers for raw water needs. The dominant rock structures in the study area are clay, alluvium, siltstone, and sandstone, as well as some rocks with suitable porosity and permeability as water carriers, such as sand and gravel. However, what appears to have considerable potential is that groundwater is found at depths of 4-53 meters in VES 1, VES 2, VES 3, VES 5, and VES 10. The results of the analysis show that the location of the Pulau Baai, Kampung Melayu Sub-district, Bengkulu City Priority Utilization Area is within the groundwater storage area, so it can be used to meet the raw water needs of the study area.

A pore-scale investigation of viscous coupling effect on immiscible two-phase flow in porous media

Viscous coupling effect plays a significant role in immiscible two-phase flow within porous media, while its influence on relative permeability remains uncertain. In this paper, an improved MRT-based viscosity-modified multicomponent multiphase (MCMP) pseudopotential lattice Boltzmann model, capable of handling high viscosity ratio, is employed to simulate two-phase flow with different viscosities in a cross-array circular structure and a real rock structure, respectively. The applicability of this model for two-phase flow with various viscosity ratios has been verified by some typical tests. Systematically, the effects of viscosity ratio, structural configuration, and wetting condition on the relative permeability curves are investigated in conjunction with their component distributions and velocity fields at different two-phase saturations. These results indicate that due to the two-phase flow competition under different structural conditions, the viscous coupling effect has varying degrees of impacts on the mobility of thin phase and viscous phase. Further, the mechanism of two-phase lubricating effect is also discussed under different wetting conditions at Darcy flow regime.

A review of multiscale numerical modeling of rock mechanics and rock engineering

AbstractRock is geometrically and mechanically multiscale in nature, and the traditional phenomenological laws at the macroscale cannot render a quantitative relationship between microscopic damage of rocks and overall rock structural degradation. This may lead to problems in the evaluation of rock structure stability and safe life. Multiscale numerical modeling is regarded as an effective way to gain insight into factors affecting rock properties from a cross‐scale view. This study compiles the history of theoretical developments and numerical techniques related to rock multiscale issues according to different modeling architectures, that is, the homogenization theory, the hierarchical approach, and the concurrent approach. For these approaches, their benefits, drawbacks, and application scope are underlined. Despite the considerable attempts that have been made, some key issues still result in multiple challenges. Therefore, this study points out the perspectives of rock multiscale issues so as to provide a research direction for the future. The review results show that, in addition to numerical techniques, for example, high‐performance computing, more attention should be paid to the development of an advanced constitutive model with consideration of fine geometrical descriptions of rock to facilitate solutions to multiscale problems in rock mechanics and rock engineering.

Characteristics of Host Rocks Manganese of The Anabanua Village Barru District South Sulawesi Province, Indonesia

One of the prospects for a wealth of geological resources that can be utilized for the benefit of mankind is manganese metal mineral resources. In Anabanua Village, especially at the research location, there is an indication of manganese mineralization that the type of rock carrying the mineralization is not yet known. The purpose of this study was to determine the chemical elements of metallic minerals, mineralogy, and types of mineralization-carrying rocks, as many as three samples in the form of chunks, which were analyzed using megascopic analysis methods identifying the texture and structure of carrier rocks, XRD (X-Ray Diffraction) analysis minerals, and XRF (X-Ray Fluorescence). The characteristics of manganese mineralization carrier rocks consist of metallic elements Mn (0.046%-20.455%), Fe (1.555%-3.673%), and nonmetallic elements SiO2 (11.403%-48.165%), K2O (0.398%-4.177%). Mineralogy of manganese mineralization carrier rocks are roeblingite (Pb2Ca6Mn2+(Si3O9)2(SO4)2(OH)2·4H2O), rhodonite (CaMn3Mn[Si5O15]), diopside (CaMgSi2O6), calcite (CaCO3), kieserite (Al, Ga)2(GeO4)(OH)2), zeolite (Mn2O.Al2O3.xSiO .yH2O), and palygorskite ((Mg, Al)2 Si4O 10(OH)·4H2O). Manganese mineralization host rocks are determined by referring to the SiO2 (%weight) and K2O(%weight) diagrams to produce basalt igneous rock types.

INVESTIGATION OF RARE AND RARE-EARTH ELEMENTS LEACHING PROCESSES FROM THE WEATHERING CRUST ORES OF THE KUNDYBAY DEPOSIT

The weathering crust ore of the Kundybay deposit is a hard-enriched ore that contains such minerals as cherchite, bastnesite, kaolinite, etc., which are embedded in the matrix structure of other silicate rocks. These minerals contain rare earth (RE) and rare earth elements (REE). The multi-component structure and dense packing of minerals in the aluminosilicate matrix of the ore make it difficult to be blocking out. Determination of the optimal ore blocking method in order to extract RE and REE is an urgent task for today. The purpose of this work is to develop an effective method of leaching RE, REE and "white soot" from the weathering crust ores of the Kundybay deposit. Methods. Processing with soluble ammonium, magnesium and aluminium salts will allow estimation of the state of RE and REE in the ore. The indicator of RE and REE blocking out during leaching with sodium hydroxide is the transfer of silica into solution, due to which RE and REE form their own hydroxides, which can be subsequently transferred into solution by acid treatment. Acid leaching with sodium pyrosulfite, hydrofluoric acid and sulfuric acid will allow the RE and REE to be immediately transferred into solution. Results and Discussion. The main characteristics of the Kundybay deposit weathering crust ore were determined, such as: silica content, moisture content, base metals, RE and REE. Conclusion. The degree of extraction of RE and REE by soluble salts of ammonium, magnesium and aluminium did not exceed 2%, which indicates the existence of RE and REE in the form of their own minerals. Leaching with sodium hydroxide showed low recovery of silica, with RE and REE remaining in a difficult-torecover form. Leaching with a mixture of sodium pyrosulfite, hydrofluoric and sulfuric acids is a promising method for RE and REE recovery due to the complex effect of the components of this mixture.

Sedimentary exhalative Pb-Zn deposit model of the Ban Lin – Phia Dam ore range, Vietnam

Purpose. The aim of this study is to construct a typical ore model for Pb-Zn ore range in Ban Lin – Phia Dam and to compare it with the existing ones. The model is constructed based on three main elements of the mineral generation process, such as formation environment, generated-ore fluid source, and deposited-ore mechanism. The obtained results are aimed at determining and predicting the resource of sedimentary exhalative (SEDEX) Pb-Zn ore deposit in Vietnam. Methods. To comprehend the characteristics of the SEDEX Pb-Zn ore deposits in the Ban Lin – Phia Dam area, we employ a combination of approaches that include mineralography to determine the composition and tectonic structure of rocks and ore, inclusion analysis to determine the temperature and inclusion composition, Pb isotopic analysis and stable S isotopic analysis to explore the environment and sources of ore-forming materials. Findings. The sedimentary exhalative (SEDEX) Pb-Zn ore deposit arises from the dissolution of metal-rich salt deposits due to heating at shallow depths (several kilometers). Subsequently, this fluid migrates through faults and fractures to preferred locations where galena and sphalerite precipitate alongside sedimentary basin deposition. Originality. SEDEX ore deposits have been extensively studied in various regions of the world. Although there are multiple theories regarding their origin, including explosion, biological debris accumulation, and surface replacement, each model has its own distinct advantages and limitations in explaining the genesis and development of SEDEX ore deposits. Practical implications. The acquired findings are intended to identify and predict the resources of sedimentary exhalative Pb-Zn ore deposits in Vietnam.

A rigid‐flexible coupled rocking structure with nonlinear stress distribution along its base: Analytical modeling, validation and parametric investigation

AbstractStructures allowing their bases to uplift and rock during an earthquake event usually experience less damage; thus, they are more resilient to seismic hazards. Accurate simulation is critical to design and seismic performance evaluation of these flexible rocking structures. Motivated from this, this article presents a new rigid‐flexible coupled rocking model for numerical evaluation of flexible rocking structures under both cyclic and dynamic loadings. The key idea in the model formulation is that the total motion of flexible rocking structure can be decomposed into a rigid‐body motion and an associated flexible deformation. The flexible deformation of the rocking body is described using displacement field of classical Timoshenko beam. Using this approach and appropriate degrees‐of‐freedom (DOFs), the governing equations‐of‐motion for flexible rocking structures are formulated using the variational principle. Multiple springs with appropriate constitutive model are distributed along rocking base to represent the inelastic behavior of rocking body. Potential post‐tensioned tendons are also considered in this model. In addition, nonlinear stress distribution along the rocking base, which cannot be considered in existing flexible rocking model, is also incorporated in the developed model. The proposed model is subsequently validated against published experimental data through quasi‐static and shake table tests, showing good accuracy when the rocking structure is relatively stockier. Finally, a parametric investigation on the effect of flexibility, nonlinear stress distribution and yield stress of rocking body on the rocking response is performed. This preliminary investigation shows that influence of flexibility and yield stress is significant while impact of nonlinear stress distribution is minor.

The mechanism and prevention of rockburst induced by the instability of the composite hard-roof coal structure and roof fractures

Despite the implementation of prevention and control measures, rockburst still occur. The characteristics of composite coal were analyzed by combining monitoring and early warning means. Two types of composite coal in the hard roof were distinguished based on their distinct initiation and destruction processes despite experiencing the same rockburst event. Rockburst caused by structural instability was studied in different mining stages, and the following conclusions were considered. 1) The mechanical model of the instability of composite coal in hard roof was established for different mining stages to reveal rockburst mechanisms. The high-speed advance of working faces changed the fractured structure of upper hard rocks in the pressure-relief protection range. The sudden fracture of long cantilever structures caused dynamic load disturbance to high-stress coal, which resulted in rockburst in the insufficient mining stage. The failure and deformation of coal under high-stress static loads created conditions for the instability and fractures of roof. The disturbance of low hard rock strata aggravated the deformations of damaged coal. Therefore, rockburst appeared in the full mining stage. 2) Structural instability and roof fracture caused by different factors were the main causes of rockburst for composite coal in hard roof in different mining stages. Relevant prevention and control measures were formulated to ensure the safety of working faces, which provided references for the prevention and control of rockburst in working faces under similar conditions.

A method for velocity modeling in a volcanic rock zone that integrates the amplitude envelope attribute and grid tomography: application to data from the North Jiangsu Basin

Abstract The complexity of volcanic rocks presents challenges in seismic data quality and fine characterization of structural morphology, making velocity modeling difficult. In this study, we propose a velocity modeling method for volcanic rock areas that combines the use of amplitude envelope attribute and grid tomography inversion. By exploiting the strong reflection characteristics of volcanic rocks, the amplitude envelope attribute is employed to identify volcanic rock structures, eliminating the need for labor-intensive manual interpretation. The grid tomography inversion technique is then used to derive background velocity fields for all strata. These are further refined by incorporating the residual velocity fields of volcanic rocks through fusion and superposition. This integrated approach significantly improves the efficiency and accuracy of velocity modeling in complex volcanic rock areas, overcoming the limitations associated with conventional volcanic rock structural models and resulting in more accurate velocity inversion. The effectiveness and feasibility of the proposed method for velocity modeling in volcanic rock zone are demonstrated through the application of seismic data from an oilfield in eastern China. The results highlight the method's strong industrial application value.

Mechanical behaviour and instability mechanism of sandstones with impact tendency under different loading paths

ABSTRACT To reveal the instability and failure mechanism of pillars with impact tendency, the Brazilian test, uniaxial compressive test, cyclic loading and unloading (CLU) tests were performed on sandstones which had impact tendency. The stress drop phenomenon was obvious in the fracture stage during the Brazilian test. The main tensile crack with penetration type appeared. Under the CLU condition, the CLU had a significant effect in enhancing the mechanical character of rocks. Under the combined loading, the character stress of rocks tested with the CLU was greater than the uniaxial loading. Moreover, the greater the first cyclic threshold, the greater the character stress of rocks. Under the uniaxial loading, the failed samples showed tensile fractures. The damage location of the samples was closely related to the Poisson effect and the ending face effect. With the electron microscope which had the same magnification, the crack deterioration degree of samples under the uniaxial loading was greater than the CLU condition. Moreover, with the first cycle threshold increasing, the surface roughness of rocks gradually decreased. The quantity of secondary micro cracks gradually decreased. The propagation and penetration phenomenon gradually disappeared. The micro structure and loading properties of rocks were affected by the bonding force and structural characters between the internal structure and matrix. Based on the micro-seismic signal of the damaging threshold, the seismic source damage location can be located. Moreover, the pillar-typed burst area can be forewarned. This research benefits predicting the pillar-typed burst disaster.

Crack propagation process in double-flawed granite under compression using digital image correlation method and numerical simulation

The existence of flaws seriously weakens the rock strength, directly affects the crack expansion morphology, and indirectly affects the slope stability. In this study, uniaxial compression tests were carried out on double-flawed granite to investigate the effects of flaw angle on its compressive strength and crack coalescence based on DIC technology and PFC2D numerical simulation. The result indicates that the UCS values of flawed specimens at angles of 15°, 30°, 45°, and 60° are approximately 25.5–51.7% lower than that of the intact specimen. The failure strength of the double-flawed sample increases with the increase of flaw inclination angle, and the fracture morphology shifts from no expansion to split expansion. When the flaw tip strain of each flaw sample exceeded 0.6%, the stress concentration was generated at the flaw tip, and the sample began to appear macroscopic large cracks. DIC technique can well observe the crack initiation and propagation process, recording inclined shear strain bands at flaw tips, tensile strain bands in the middle of flaws, and tensile shear O-shaped strain ring in rock bridge areas. Numerical simulations using PFC2D software were carried out and showed a good agreement with the physical results. In addition, the effect of structural inclination on specimen failure strength is clearly explained by the theory of compressive damage based on the projected size of flaws. The deformation of rock mass is not a simple material deformation but is composed of material deformation and structural deformation. These experimental and numerical results enhance our understanding of crack initiation and coalescence characteristics, aiding in the analysis of rock structure stability in scenarios such as excavated underground openings, slopes, and tunneling construction, where pre-existing cracks or step-path fractures are pivotal to structural integrity.

Effects of rock pore and micromorphology on electromagnetic radiation characteristics

Electromagnetic radiation (EMR) is a crucial tool for monitoring and early warning of underground engineering disasters. Investigating the inherent pore characteristics of rocks is essential for a comprehensive understanding of the EMR phenomenon. The EMR was monitored during various types of rock splitting failures. The pore structure and micromorphology of rocks are studied using quantitative methods such as mercury intrusion porosimetry, fractal analysis, and the Gray Level Co-occurrence Matrix (GLCM). Results indicate that the fractal dimension of red sandstone is significantly lower than the other three rocks. The fractal complexity increases sequentially from red sandstone to marble, granite, and limestone. As the fractal dimension decreases, the signal waveform characteristics of the four rocks become more complex before the main fracture, with a significant increase in signals during the compaction and elasticity stages. Higher fractal dimensions lead to a shift in energy and count from elasticity stage to the post-peak stage. The main fracture amplitudes of the four rocks generally exhibited a consistent pattern, following the sequence of granite > marble > limestone > red sandstone. The main fracture amplitude decreases with increasing complexity of the rock's pore micromorphology. Rock pore characteristics affect frequency domain characteristics by influencing rock strength and crack expansion. An increase in the average pore diameter tends to decrease both the main and center frequencies.

Evolution and Quantitative Characterization of Stress and Displacement of Surrounding Rock Structure due to the Multiple Layers Backfill Mining under Loose Aquifers

Evolution and Quantitative Characterization of Stress and Displacement of Surrounding Rock Structure due to the Multiple Layers Backfill Mining under Loose Aquifers

Fluid flow in the Katanga Supergroup: From Lufilian brittle tectonic stages to the post-Lufilian period (Democratic Republic of Congo)

The metasedimentary rock succession of the Neoproterozoic-Cambrian Katanga Supergroup in the Central Africa Copperbelt shows evidence of several complex tectonic events. The deformation of this supergroup started from the tectonic inversion at about 570 Ma and lasted up to today, but reached paroxysm at ∼550 Ma. This long period was characterized by folding and faulting throughout multiple compressive and extensional events, which controlled the regional fluid flow on the one hand, and played an important role during formation of the stratiform to stratabound Cu-Co (Ni, U) deposits and the polymetallic Cu-Zn-Pb (Ag, Ge, Mo, Cd) vein type deposits on the other hand. Based on the structural analysis and paleostress reconstruction, coupled with fluid inclusion characterization from mineralized structures in rocks from the Nguba, Kundelungu and Biano Groups, this study demonstrates that the composition of hydrothermal fluids changed during brittle tectonic deformation during the Lufilian orogeny and subsequent uplift and post-Lufilian faulting.During early brittle tectonic deformation along strike slip faults with sinistral and dextral movement related to a NE-SW transpression, the Cu-mineralizing fluid was hypersaline (27.9–31.1 eq. wt% NaCl) with moderate temperatures (Th = 128–216 °C). The subsequent Cu or Cu (Zn, Pb) mineralization formed within an E-W extensional stress regime, related to the late Lufilian orogenic collapse. The fluid inclusions present in the gangue minerals associated with this latter mineralization show a large range in Th (50–264 °C) and salinity (26.7–36.0 eq. wt% NaCl). The decrease in temperature is interpreted to be due to migration of the fluids at shallower depth in the subsurface after uplift and erosion of the orogen. The increased salinity of the fluid is related to the dissolution of evaporites, mainly NaCl. A second H2O-NaCl-CaCl2 fluid with a homogenization temperature below 55 °C has also been found associated with this brittle stage and mineralization phase, but only in rocks belonging to the Kundelungu Group. A third mineralization phase, also characterized by Cu or Cu (Zn, Pb), formed during the post-Lufilian period within a NW-SE transpressional inversion regime. The fluid inclusion in the gangue minerals of this mineralization phase have a smaller range in homogenization temperature (Th = 37–172 °C) and the largest range in salinity (0.71–30 eq. wt% NaCl), compared to the earlier fluid inclusions generations. This large range in salinity may be explained by the mixing of a high salinity fluid, already present during the earlier tectonic stages in the sedimentary basin, with meteoric water. During the more recent rift-related extension, a fluid with again a large and higher range in homogenization temperatures (Th = 47–257 °C) and with a typical low salinity (<10 eq. wt% NaCl) has been recognized in minerals filling NNE-SSW to NE-SW oriented faults and fractures. The upward migration of a relatively low-salinity fluid from deeper parts in the subsurface explains the variation in the temperatures observed with this tectonic event.