A pore-scale numerical model for flow through porous media

A pore-scale numerical model based on Smoothed Particle Hydrodynamics (SPH) is described for modelling fluid flow phenomena in porous media. Originally developed for astrophysics applications, SPH is extended to model incompressible flows of low Reynolds number as encountered in groundwater flow systems. In this paper, an overview of SPH is provided and the required modifications for modelling flow through porous media are described, including treatment of viscosity, equation of state, and no-slip boundary conditions. The performance of the model is demonstrated for two-dimensional flow through idealized porous media composed of spatially periodic square and hexagonal arrays of cylinders. The results are in close agreement with solutions obtained using the finite element method and published solutions in the literature. Copyright © 1999 John Wiley & Sons, Ltd.     

Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Theoretical and numerical study of the steady‐state flow through finite fractured porous media

This paper investigates the two‐dimensional flow problem through an anisotropic porous medium containing several intersecting curved fractures. First, the governing equations of steady‐state fluid flow in a fractured porous body are summarized. The flow follows Darcy's law in matrix and Poiseuille's law in fractures. An infinite transversal permeability is considered for the fractures. A multi‐region boundary element method is used to derive a general pressure solution as a function of discharge through the fractures and the pressure and the normal flux on the domain boundary. The obtained solution fully accounts for the interaction and the intersection between fractures. A numerical procedure based on collocation method is presented to compute the unknowns on the boundaries and on the fractures. The numerical solution is validated by comparing with finite element solution or the results obtained for an infinite matrix. Pressure fields in the matrix are illustrated for domains containing several interconnected fractures, and mass balance at the intersection points is also checked. Copyright © 2013 John Wiley & Sons, Ltd.     

Multiphase SPH modeling of free surface flow in porous media with variable porosity

This paper presents a numerical model for simulating free surface flow in porous media with spatially varying porosity. The governing equations are based on the mixture theory. The resistance forces between solid and fluid is assumed to be nonlinear. A multiphase SPH approach is presented to solve the governing equations. In the multiphase SPH, water is modeled as a weakly compressible fluid, and solid phase is discretized by fixed solid particles carrying information of porosity. The model is validated by several numerical examples including seepage through specimen, fast flow through rockfill dam and wave interaction with porous structure. Good agreements between numerical results and experimental data are obtained in terms of flow rate and evolution of free surface. Parameter study shows that (1) the nonlinear resistance law provides more accurate results; (2) particle size and porosity have significant influence on the porous flow.     

Construction of three-phase data to model multiphase flow in porous media: Comparing an optimization approach to the finite element approach

Multiphase flow modelling is a major issue in the assessment of groundwater pollution. Three-phase flows are commonly governed by mathematical models that associate a pressure equation with two saturation equations. These equations involve a number of secondary variables that reflect the fluid behaviour in a porous medium. To improve the computational efficiency of multiphase flow simulators, several simplified reformulations of three-phase flow equations have been proposed. However, they require the construction of new secondary variables adapted to the reformulated flow equations. In this article, two different approaches are compared to quantify these variables. A numerical example is given for a typical fine sand.     

Time‐scale analysis in unsaturated porous media under external wave loads

Extreme waves caused by tsunamis and storm surges can lead to soil failures in the near‐shore region, which may have severe impact on coastal environments and communities. Multiphase flows in deformable porous media involve several coupled processes and multiple time scales, which are challenging for numerical simulations. The objective of this study is to investigate the roles of the various processes and their interactions in multiphase flows in unsaturated soils under external wave loading, via theoretical time‐scale analysis and numerical simulations. A coupled geomechanics–multiphase flow model based on conservation laws is used. Theoretical analysis based on coupled and decoupled models demonstrates that transient and steady‐state responses are governed by pore pressure diffusion and saturation front propagation, respectively, and that the two processes are essentially decoupled. Numerical simulations suggest that the compressibility of the pore fluids and the deformation of the soil skeleton are important when the transient responses of the media are of concern, while the steady‐state responses are not sensitive to these factors. The responses obtained from the fully coupled numerical simulations are explained by a simplified time‐scale analysis based on coupled and decoupled models. Copyright © 2010 John Wiley & Sons, Ltd.     

LBM–DEM modeling of fluid–solid interaction in porous media

Three porous media flow problems, in which the fluid mechanical interactions are critical, are studied in a mesoscopic–microscopic coupling system. In this system, fluid flow in the pore space is explicitly modeled at mesoscopic level by the lattice Boltzmann method, the geometrical representation and the mechanical behavior of the solid skeleton are modeled at microscopic level by the particulate distinct element method (DEM), and the interfacial interaction between the fluid and the solids is resolved by an immersed boundary scheme. In the first benchmark problem, the well‐known and frequently utilized Ergun equation is validated in periodic particle and periodic pore models. In the second problem, the upward seepage problem is simulated over three stages: The settlement of the column of sphere under gravity loading is measured to illustrate the accuracy of the DEM scheme; the system is solved to hydrostatic state with pore space filled with fluid, showing that the buoyancy effect is captured correctly in the mesoscopic–microscopic coupling system; then, the flow with constant rate is supplied at the bottom of the column; the swelling of the ground surface and pore pressure development from the numerical simulation are compared with the predictions of the macroscopic consolidation theory. In the third problem, the fluid‐flow‐induced collapse of a sand arch inside a perforation cavity is tested to illustrate a more practical application of the developed system. Through comparing simulation results with analytical solutions, empirical law and physical laboratory observations, it is demonstrated that the developed lattice Boltzmann–distinct element coupling system is a powerful fundamental research tool for investigating hydromechanical physics in porous media flow. Copyright © 2012 John Wiley & Sons, Ltd.     

Effective strength of saturated double porous media with a Drucker–Prager solid phase

This study aims at determining the macroscopic strength of porous materials having a Drucker–Prager solid phase at microscale and two populations of saturated pores with different pressures at both micro and meso scales. To this end, and taking account of the available results by Maghous et al. (2009), we first derive a closed‐form expression of approximate criterion for a dry porous medium whose matrix obeys to a general elliptic criterion. The methodology to formulate this criterion is based on limit analysis of a hollow sphere subjected to a uniform strain rate boundary conditions. The obtained results are then implemented in a two‐step homogenization procedure, which interestingly delivers analytical expression of the macroscopic criterion for dry double porous media whose solid phase at microscale obeys to a Drucker–Prager criterion. After a brief discussion of the results, we propose an extension to double porous saturated media, allowing therefore to quantify the simultaneous effects of the different pore pressures applied on each voids population. The results are discussed in terms of the existence or not of effective stresses. Finally, they are assessed by comparing them to recently available results. Copyright © 2013 John Wiley & Sons, Ltd.     

Development of a discontinuous approach for modeling fluid flow in heterogeneous media using the numerical manifold method

In the numerical modeling of fluid flow in heterogeneous geological media, large material contrasts associated with complexly intersected material interfaces are challenging, not only related to mesh discretization but also for the accurate realization of the corresponding boundary constraints. To address these challenges, we developed a discontinuous approach for modeling fluid flow in heterogeneous media using the numerical manifold method (NMM) and the Lagrange multiplier method (LMM) for modeling boundary constraints. The advantages of NMM include meshing efficiency with fixed mathematical grids (covers), the convenience of increasing the approximation precision, and the high integration precision provided by simplex integration. In this discontinuous approach, the elements intersected by material interfaces are divided into different elements and linked together using the LMM. We derive and compare different forms of LMMs and arrive at a new LMM that is efficient in terms of not requiring additional Lagrange multiplier topology, yet stringently derived by physical principles, and accurate in numerical performance. To demonstrate the accuracy and efficiency of the NMM with the developed LMM for boundary constraints, we simulate a number of verification and demonstration examples, involving a Dirichlet boundary condition and dense and intersected material interfaces. Last, we applied the developed model for modeling fluid flow in heterogeneous media with several material zones containing a fault and an opening. We show that the developed discontinuous approach is very suitable for modeling fluid flow in strongly heterogeneous media with good accuracy for large material contrasts, complex Dirichlet boundary conditions, or complexly intersected material interfaces. Copyright © 2015 John Wiley & Sons, Ltd.     

Genetic significance of fluid inclusions in the CSA Cu–Pb–Zn deposit, Cobar, Australia

CSA mine exploits a ‘Cobar-type’ Cu–Pb–Zn±Au±Ag deposit within a cleaved and metamorphosed portion of the Cobar Supergroup, central New South Wales. The deposit comprises systems of ‘lenses’ that encompass veins, disseminations and semi-massive to massive Cu–Pb–Zn ores. The systems and contained lenses truncate bedding, are approximately coplanar with regional cleavage and similarly oriented shear zones and plunge parallel to the elongation lineation. Systems have extreme vertical continuity (>1000 m), short strike length (400 m) and narrow width (100 m), exhibit vertical and lateral ore-type variation and have alteration haloes. Models of ore formation include classical hydrothermalism, structurally controlled remobilisation and polymodal concepts; syntectonic emplacement now holds sway.Fluid inclusions were examined from quartz±sulphide veins adjacent to now-extracted ore, from coexisting quartz–sulphide within ore, and from vughs in barren quartz veins. Lack of early primary inclusions precluded direct determination of fluids associated with D₂–D₃ ore and vein emplacement. Similarly, decrepitation (by near-isobaric heating) of the two oldest secondary populations precluded direct determination of fluid phases immediately following D₂–D₃ ore and vein emplacement. Post-decrepitation outflow (late D₃ to early post-D₃) is recorded by monophase CH₄ inclusions. Entrained outflow of deeply circulated meteoric fluid modified the CH₄ system; modification is recorded by H₂O+CH₄ and H₂O+(trace CH₄) secondary populations and by an H₂O+(trace CH₄) primary population. The contractional tectonics (D₂–D₃) of ore emplacement was superseded by relaxational tectonics (D_4P) that facilitated meteoric water penetration and return flow.Under D₂ prograde metamorphism, entrapment temperatures (Tt) and pressures (Pt) for pre-decrepitation secondary inclusions are estimated as Tt300–330 °C and Pt1.5–2 kbar≈Plith (the lithostatic pressure). Decrepitation accompanied peak metamorphism (T350–380 °C) in mid- to late-D₃, while in late-D₃ to early post-D₃, essentially monophase CH₄ secondary inclusions were entrapped at Tt350 °C and Pt=1.5–2 kbar≈Plith. Subsequently, abundant CH₄ and entrained meteoric water were entrapped as H₂O+CH₄ secondaries under slowly decreasing temperature (Tt330–350 °C) and constant pressure (Pt1.5–2 kbar). Finally, with increasingly dominant meteoric outflow, H₂O+(trace CH₄) populations record decreasing temperatures (Tt>300 to <350 down to 275–300 °C) at pressures of Phydrostatic<Pt (1 kbar) <Plith (1.5 kbar).The populations of inclusions provide insight into fluid types, flow regimes and P–T conditions during parts of the deposit's evolution. They indirectly support the role of basin-derived CH₄ fluids in ore formation, but provide no insight into a basement-sourced ore-forming fluid. They fully support post-ore involvement of meteoric water. The poorly constrained entrapment history is believed to span 10 Ma from 395 to 385 Ma.     

Fault–valve action and vein development during strike–slip faulting: An example from the Ribeira Shear Zone, Southeastern Brazil

Fluid inclusion microthermometry and structural data are presented for quartz vein systems of a major dextral transcurrent shear zone of Neoproterozoic–Cambrian age in the Ribeira River Valley area, southeastern Brazil. Geometric and microstructural constraints indicate that foliation–parallel and extensional veins were formed during dextral strike–slip faulting. Both vein systems are formed essentially by quartz and lesser contents of sulfides and carbonates, and were crystallized in the presence of CO₂–CH₄ and H₂O–CO₂–CH₄–NaCl immiscible fluids following unmixing from a homogeneous parental fluid. Contrasting fluid entrapment conditions indicate that the two vein systems were formed in different structural levels. Foliation–parallel veins were precipitated beneath the seismogenic zone under pressure fluctuating from moderately sublithostatic to moderately subhydrostatic values (319–397 °C and 47–215 MPa), which is compatible with predicted fluid pressure cycle curves derived from fault–valve action. Growth of extensional veins occurred in shallower structural levels, under pressure fluctuating from near hydrostatic to moderately subhydrostatic values (207–218 °C and 18–74 MPa), which indicate that precipitation occurred within the near surface hydrostatically pressured seismogenic zone. Fluid immiscibility and precipitation of quartz in foliation–parallel veins resulted from fluid pressure drop immediately after earthquake rupture. Fluid immiscibility following a local pressure drop during extensional veining occurred in pre-seismic stages in response to the development of fracture porosity in the dilatant zone. Late stages of fluid circulation within the fault zone are represented dominantly by low to high salinity (0.2 to 44 wt.% equivalent NaCl) H₂O–NaCl–CaCl₂ fluid inclusions trapped in healed fractures mainly in foliation–parallel veins, which also exhibit subordinate H₂O–NaCl–CaCl₂, CO₂–(CH₄) and H₂O–CO₂–(CH₄)–NaCl fluid inclusions trapped under subsolvus conditions in single healed microcracks. Recurrent circulation of aqueous–carbonic fluids and aqueous fluids of highly contrasting salinities during veining and post-veining stages suggests that fluids of different reservoirs were pumped to the ruptured fault zone during faulting episodes. A fluid evolution trending toward CH₄ depletion for CO₂–CH₄–bearing fluids and salinity depletion and dilution (approximation of the system H₂O–NaCl) for aqueous–saline fluids occurred concomitantly with decrease in temperature and pressure related to fluid entrapment in progressively shallower structural levels reflecting the shear zone exhumation history.     

Behavior of sand–rubber particle mixtures: experimental observations and numerical simulations

A series of laboratory experiments and numerical simulations are conducted to explore the characteristics of mixtures composed of sand and rubber particles of the same median diameter. The mixtures are prepared with different volumetric sand fractions (sf = V_sand/V_total). The experiment focuses on assessing the strain level on the characteristics of the mixture with the volume fraction of each component. Numerical simulations using the discrete element method are performed to obtain insight into the microscale behavior and internal mechanism of the mixtures. The experimental results show that the behavior of the mixtures is dependent on the relative sand and rubber particles composition with variation in the strain levels. The numerical simulation reveals the effect of the soft rubber particle inclusion in the mixture on the micromechanical parameters. In low sand fraction mixtures, a high shear stress along the contact is mobilized, and the stress state is driven to a more anisotropic condition because of the relatively high particle friction angle of the rubber. The rubber particles play different roles with the strain level in the mixture, including increasing the coordination number and controlling plasticity of the mixture in a small strain, preventing buckling of the force chain in an intermediate strain, and leading to contractive behavior in a large strain. Copyright © 2014 John Wiley & Sons, Ltd.     

Evolution of bulk composition, mineralogy, strain style and fluid flow during an HT–LP metamorphic event: sillimanite zone of the Catalan Coastal Ranges Variscan basement, NE Iberia

Metamorphic rocks in the Osor complex (Guilleries massif, NE Iberian Peninsula) show the following structural and compositional features: strong differentiation into quartz-rich gneissic semipelitic and quartz-absent, mica-rich schistose bands, higher density of igneous (both basic and leucogranitic) and quartz veins in the schistose domains and strong strain partitioning in the pelitic bands. Garnet is present in both kinds of lithologies, showing also differential textural and chemical features interpreted to be dependent on bulk composition, deformation and fluid interaction histories. Textures, mineral composition and thermobarometry suggest the operation of concurrent mechanical, mass transfer and thermal phenomena such as: (1) variations in strain style, (2) fluid infiltration, (3) magmatic injection and (4) HT–LP metamorphic and metasomatic episodes. The following sequence of events is suggested: initially the cooling of syntectonic high-T basic quartz diorite sheets promoted high strain rates, low dP/dT thermobaric evolution, incipient anatexis in the pelitic bands and devolatilization through a pervasive to vein-channelized prograde fluid flow. The prograde flow enhanced an ongoing compositional tectono-metamorphic differentiation and produced metasomatism through depletion of the Osor rocks in alkalis and calcium. Later injection and cooling of peraluminous leucogranitoid sheets, preferentially along pelitic bands, increased the ratio of magmatic/metamorphic components in the fluids and strongly enriched them in alkalis producing a second metasomatic episode. During crystallization of quartz and leucogranitoid veins, the pelitic bands were strongly enriched back again in alkalis, promoting the blastesis of big crystals of post-peak muscovite and albite as well as the retrogression of garnet. The metasomatic mica-rich levels must have been the preferred locus for development of a new deformation style dominated by shear band fabrics in metapelites and related to a release of the gravitational instability originated previously due to crustal thickening. The increasing decompressional component of the retrograde P–T path also suggests that this style of deformation was prevalent during, if not responsible for, a phase of exhumation of the metamorphic complex. It is suggested that similar patterns of thermomechanical and mass transfer phenomena could well be a fundamental characteristic common to all HT–LP metamorphic belts.     

Fluid permeability of sedimentary rocks in a complete stress–strain process

The permeability of sedimentary rocks during triaxial compression tests was investigated to relate it to the complete strain–stress process. It was found that the permeability was not constant, but varied with the stress and strain states in the rocks. Prior to the peak strength, the permeability decreases with increasing load. A dramatic increase in permeability occurs during the strain softening period. In the present study, in situ measurements of fluid flow and pressure in floor strata was carried out in a double longwall mining face in the Yangzhuang colliery. These measurements show that both the strata pressure and the position with respect to the mining face influence the hydrogeologic properties. The permeability increased in the floor strata behind the mining face because those mining induced fractures opened as the strata pressure decreased. To better understand this change in hydraulic behavior around the mining faces, 3-D numerical modeling was carried out. The model provides the general picture of the stress distribution and failure zone both in the floor and roof strata. The field and model results demonstrate the importance of changes in the stress and strain states on the hydrogeology of a site.     

Genesis of the Weiquan Ag–Polymetallic Deposit in East Tianshan, China: Evidence from Zircon U–Pb Geochronology and C–H–O–S–Pb Isotope Systematics

The Weiquan Ag-polymetallic deposit is located on the southern margin of the Central Asian Orogenic Belt and in the western segment of the Aqishan-Yamansu arc belt in East Tianshan,northwestern China. Its orebodies, controlled by faults, occur in the lower Carboniferous volcanosedimentary rocks of the Yamansu Formation as irregular veins and lenses. Four stages of mineralization have been recognized on the basis of mineral assemblages, ore fabrics, and crosscutting relationships among the ore veins. Stage I is the skarn stage(garnet + pyroxene), Stage Ⅱ is the retrograde alteration stage(epidote + chlorite + magnetite ± hematite 士 actinolite ± quartz),Stage Ⅲ is the sulfide stage(Ag and Bi minerals + pyrite + chalcopyrite + galena + sphalerite + quartz ± calcite ± tetrahedrite),and Stage IV is the carbonate stage(quartz + calcite ± pyrite). Skarnization,silicification, carbonatization,epidotization,chloritization, sericitization, and actinolitization are the principal types of hydrothermal alteration. LAICP-MS U-Pb dating yielded ages of 326.5±4.5 and 298.5±1.5 Ma for zircons from the tuff and diorite porphyry, respectively. Given that the tuff is wall rock and that the orebodies are cut by a late diorite porphyry dike, the ages of the tuff and the diorite porphyry provide lower and upper time limits on the age of ore formation. The δ~(13)C values of the calcite samples range from-2.5‰ to 2.3‰, the δ~(18)O_(H2 O) and δD_(VSMOW) values of the sulfide stage(Stage Ⅲ) vary from 1.1‰ to 5.2‰ and-111.7‰ to-66.1‰, respectively,and the δ~(13)C, δ~(18)O_(H2 O) and δD_(V-SMOW) values of calcite in one Stage IV sample are 1.5‰,-0.3‰, and-115.6‰, respectively. Carbon, hydrogen, and oxygen isotopic compositions indicate that the ore-forming fluids evolved gradually from magmatic to meteoric sources. The δ~(34)S_(V-CDT) values of the sulfides have a large range from-6.9‰ to 1.4‰, with an average of-2.2‰, indicating a magmatic source, possibly with sedimentary contributions. The ~(206)Pb/~(204)Pb, ~(207)Pb/~(204)Pb, and ~(208)Pb/~(204)Pb ratios of the sulfides are 17.9848-18.2785,15.5188-15.6536, and 37.8125-38.4650, respectively, and one whole-rock sample at Weiquan yields~(206)Pb/~(204)Pb,~(207)Pb/~(204)Pb, and ~(208)Pb/~(204)Pb ratios of 18.2060, 15.5674, and 38.0511,respectively. Lead isotopic systems suggest that the ore-forming materials of the Weiquan deposit were derived from a mixed source involving mantle and crustal components. Based on geological features, zircon U-Pb dating, and C-H-OS-Pb isotopic data, it can be concluded that the Weiquan polymetallic deposit is a skarn type that formed in a tectonic setting spanning a period from subduction to post-collision. The ore materials were sourced from magmatic ore-forming fluids that mixed with components derived from host rocks during their ascent, and a gradual mixing with meteoric water took place in the later stages.     

Geochronological framework and Pb, Sr isotope geochemistry of the Qingchengzi Pb–Zn–Ag–Au orefield, Northeastern China

The Qingchengzi orefield in northeastern China, is a concentration of several Pb–Zn, Ag, and Au ore deposits. A combination of geochronological and Pb, Sr isotopic investigations was conducted. Zircon SHRIMP U–Pb ages of 225.3 ± 1.8 Ma and 184.5 ± 1.6 Ma were obtained for the Xinling and Yaojiagou granites, respectively. By step-dissolution Rb–Sr dating, ages of 221 ± 12 Ma and 138.7 ± 4.1 Ma were obtained for the sphalerite of the Zhenzigou Zn–Pb deposit and pyrargyrite of the Ag ore in the Gaojiabaozi Ag deposit, respectively. Pb isotopic ratios of the Ag ore at Gaojiabaozi (²⁰⁶Pb/²⁰⁴Pb = 18.38 to 18.53) are higher than those of the Pb–Zn ores (²⁰⁶Pb/²⁰⁴Pb = 17.66 to 17.96; Chen et al. [Chen, J.F., Yu, G., Xue, C.J., Qian, H., He, J.F., Xing, Z., Zhang, X., 2005. Pb isotope geochemistry of lead, zinc, gold and silver deposit clustered region, Liaodong rift zone, northeastern China. Science in China Series D 48, 467–476.]). Triassic granites show low Pb isotopic ratios (²⁰⁶Pb/²⁰⁴Pb = 17.12 to 17.41, ²⁰⁷Pb/²⁰⁴Pb = 15.47 to 15.54, ²⁰⁸Pb/²⁰⁴Pb = 37.51 to 37.89) and metamorphic rocks of the Liaohe Group have high ratios (²⁰⁶Pb/²⁰⁴Pb = 18.20 to 24.28 and 18.32 to 20.06, ²⁰⁷Pb/²⁰⁴Pb = 15.69 to 16.44 and 15.66 to 15.98, ²⁰⁸Pb/²⁰⁴Pb = 37.29 to 38.61 and 38.69 to 40.00 for the marble of the Dashiqiao Formation and schist of the Gaixian Formation, respectively).Magmatic activities at Qingchengzi and in adjacent regions took place in three stages, and each contained several magmatic pulses: ca. 220 to 225 Ma and 211 to 216 Ma in the Triassic; 179 to 185 Ma, 163 to 168 Ma, 155 Ma and 149 Ma in the Jurassic, as well as ca. 140 to 130 Ma in the Early Cretaceous. The Triassic magmatism was part of the Triassic magmatic belt along the northern margin of the North China Craton produced in a post-collisional extensional setting, and granites in it formed by crustal melting induced by mantle magma. The Jurassic and Early Cretaceous magmatism was related to the lithospheric delamination in eastern China. The Triassic is the most important metallogenic stage at Qingchengzi. The Pb–Zn deposits, the Pb–Zn–Ag ore at Gaojiabaozi, and the gold deposits were all formed in this stage. They are temporally and spatially associated with the Triassic magmatic activity. Mineralization is very weak in the Jurassic. Ag ore at Gaojiabaozi was formed in the Early Cretaceous, which is suggested by the young Rb–Sr isochron age, field relations, and significantly different Pb isotopic ratios between the Pb–Zn–Ag and Ag ores. Pb isotopic compositions of the Pb–Zn ores suggest binary mixing for the source of the deposits. The magmatic end-member is the Triassic granites and the other metamorphic rocks of the Liaohe Group. Slightly different proportions of the two end-members, or an involvement of materials from hidden Cretaceous granites with slightly different Pb isotopic ratios, is postulated to interpret the difference of Pb isotopic compositions between the Pb–Zn–(Ag) and Ag ores. Sr isotopic ratios support this conclusion. At the western part of the Qingchengzi orefield, hydrothermal fluid driven by the heat provided by the now exposed Triassic granites deposited ore-forming materials in the low and middle horizons of the marbles of the Dashiqiao Formation near the intrusions to form mesothermal Zn–Pb deposits. In the eastern part, hydrothermal fluids associated with deep, hidden Triassic intrusions moved upward along a regional fault over a long distance and then deposited the ore-forming materials to form epithermal Au and Pb–Zn–Ag ores. Young magmatic activities are all represented by dykes across the entire orefield, suggesting that the corresponding main intrusion bodies are situated in the deep part of the crust. Among these, only intrusions with age of ca. 140 Ma might have released sufficient amounts of fluid to be responsible for the formation of the Ag ore at Gaojiabaozi.Our age results support previous conclusions that sphalerite can provide a reliable Rb–Sr age as long as the fluid inclusion phase is effectively separated from the “sulfide” phase. Our work suggests that the separation can be achieved by a step-resolution technique. Moreover, we suggest that pyrargyrite is a promising mineral for Rb–Sr isochron dating.     

Genesis of Pb–Zn Mineralization Beneath the Xiangshan Uranium Orefield,South China: Constraints from H–O–S–Pb Isotopes and Rb–Sr Dating

The volcanic‐hosted Xiangshan uranium orefield is the largest uranium deposit in South China. Recent exploration has discovered extensive Pb–Zn mineralization beneath the uranium orebodies. Detailed geological investigation reveals that the major metallic minerals include pyrite, sphalerite, galena, and chalcopyrite, whilst the major non‐metallic minerals include quartz, sericite, and calcite. New δ¹⁸O_fluid and δD_fluid data indicate that the ore‐forming fluids were mainly derived from magmatic, and the sulfide δ³⁴S values (2.2–6.9‰) suggest a dominantly magmatic sulfur source. The Pb isotope compositions are homogeneous (²⁰⁶Pb/²⁰⁴Pb = 18.120–18.233, ²⁰⁷Pb/²⁰⁴Pb = 15.575–15.698, and ²⁰⁸Pb/²⁰⁴Pb = 37.047–38.446). The ⁸⁷Sr/⁸⁶Sr ratios of sulfide minerals range from 0.7197 to 0.7204, which is much higher than volcanic rocks and fall into the range of metamorphic basement. Lead and strontium isotopic compositions indicate that the metallogenic materials probably were derived from metamorphic basement. Pyrite Rb–Sr dating of the ores yielded 131.3 ± 4.0 Ma, indicating that the Pb–Zn mineralization occurred in the Early Cretaceous.     

Evolution of mineralizing fluids in the Zn–Pb–Cu(–Ag ± Au) skarn and epithermal deposits of the world-class San Martín district, Zacatecas, Mexico

Mineralizing fluids at the San Martín skarn show an evolution characterized by prograde and retrograde associations. The prograde mineral associations consist of (1) a massive garnet zone, (2) a tremolite ± garnet zone, and (3) a late association of quartz, sphalerite, calcite and fluorite lining the vugs in the garnet zone. The fluids of the prograde associations exhibit decreasing temperatures of homogenization (Th) and variable salinities. The fluids of the massive garnet zone have salinities of 36 wt.% NaCl equiv. and Th of 645 to 570 °C, corresponding to pressures of 1055 bar. At the tremolite ± garnet zone, Th range from 438 to 354 °C. In the late association at the endoskarn, the following evolution can be drawn: (a) salinities of 50 to 42 wt.% NaCl equiv., and Th of 455 to 346 °C in quartz, (b) salinities of 46 wt.% NaCl equiv., and Th of 415 to 410 °C in sphalerite, (c) salinities of 50 to 37 wt.% NaCl equiv., and Th of 479 to 310 °C in calcite, (d) salinities of 33 to 28 wt.% NaCl equiv. and of 24 to 22 wt.% KCl in fluorite, and (e) two types of fluids with salinities of 2 and 39 wt.% NaCl equiv. and Th 344 and 300 °C, respectively, in later saccharoidal quartz segregations. The retrograde mineral associations comprise pervasive propylitic alteration to carbonization, and mantos with sulfides. Fluids in epidote have salinities of 7.6 wt.% NaCl equiv. and Th of 287 to 252 °C, and in calcite have salinities of 9.2 to 1 wt.% NaCl equiv. and Th of 188 to 112 °C. Fluids in the sulfide assemblages in the mantos have salinities of 8 to 3 wt.% NaCl equiv. and Th 300 °C, with corresponding pressures of 94 bar. Fluids in late epithermal veins close to the intrusive body have salinities of 10 to 5 wt.% NaCl equiv. and Th of 275 to 200 °C, and distal veins show salinities of 2 to 1 wt.% NaCl equiv. and Th of 160 °C.