Characteristics and Mechanism of the Ore-Forming Fluids in the Shimensi Tungsten Polymetallic Deposit in Southeastern China

: The Shimensi super-large tungsten polymetallic deposit is located in the Late Jurassic– Early Cretaceous Porphyry–Skarn tungsten ore belt in the south Yangtze metallogenic belt. There are three types of mineralization: veinlet-disseminated type, thick quartz vein type and hydrothermal cryptoexplosive breccia type. Based on geological studies, this paper presents new petro-graphic, microthermometric, laser Raman spectroscopic and hydrogen and oxygen isotope research on the fluid inclusions from the deposit. The results show that there are five different types of fluid inclusions: liquid-rich inclusions, vapor-rich inclusions, pure liquid inclusions, pure vapor inclusions, and fluid inclusions containing a solid crystal. The homogenization temperatures of the fluid inclusion range from 140 °C to 270 °C, the salinities are 3 wt.%–5 wt.% NaCl eq and the densities of ore-forming fluid range from 0.64 g/cm 3 to 0.99 g/cm 3 . For the analyses of laser Raman spectroscopy, the ore-forming fluids can be approximated by a Ca 2+ -Na + -SO 42− -Cl − fluid system with small amounts of CO 2 , CH 4 and N 2 . Otherwise, the data of the pressure, pH and Eh show a fluid metallogenic environment of low pressure, weak acid and weak reduction. The values of the homogenization temperature in the three types of orebodies suggest that the mineralization is characterized by a decrease in temperature under the conditions of fluid immiscibility. The H-O isotope values are interpreted to indicate that the ore-forming fluids are mainly composed of magmatic water, and meteoric water is added with the process of magma rising.


Introduction
Jiangxi province is an important W producer in China and hosts several giant tungsten deposits in its southern region, such as Xihuashan and Dajishan [1][2][3][4].With the development of geological prospecting, a cluster of world class W deposits have been discovered in the northern Jiangxi province, including the Dahutang W-Cu-Mo, Xianglushan W, Yangchuling W-Mo and Zhuxi W-Cu-Mo deposits [5].The discovery of these tungsten polymetallic deposits has gradually changed the distribution patterns in Jiangxi province, and even in China.The Dahutang super-large tungsten polymetallic ore field is located in the Late Jurassic-Early Cretaceous porphyry-skarn tungsten ore zone of the Jiangnan Orogen between the Yangtze and Cathaysia Craton [6,7].The Dahutang ore field has estimated reserves of 2 million tons of WO3 with 0.152% W, accompanied by 500 Kt Cu with 0.12% Cu and 80.2 Kt Mo with 0.098% Mo [8].The discovery of the Dahutang polymetallic ore district, especially the veinlet-disseminated scheelite mineralization in the biotite granodiorite of the Shimensi deposit, has changed the exploration strategy of focusing on the quartz vein type wolframite mineralization [9].
The Shimensi tungsten deposit itself is situated at the junction of Wuning, Xiushui and Jing'an Counties in the north of the Dahutang ore field.By May 2012, the 916 Corps of Jiangxi Bureau of Geological Exploration for Mining Resources had invested 77 drills, completed a drilling footage of 190 km and submitted 743,000 tons of WO3 resources with an average grade of 0.195%, accompanied by 4000 Kt Cu and 28 Kt Mo [10].Although the Shimensi tungsten deposit has been studied extensively in recent years, most of the studies have focused on the geological characteristics, metallogenic chronology, mineralization background, ore-forming material sources, etc. [11][12][13][14][15][16][17][18][19][20][21][22][23][24][25], and little is known about the characteristics of ore-forming fluids and their genesis in the region.
Most mineralization in the Earth's system involves geological fluids, which exert an important influence on the genesis, transport and precipitation of mineralized materials [26,27].Therefore, the fluid effects associated with mineralization have attracted many monographs and in-depth studies by scholars.Based on previous results, this paper examines fluid inclusions in three main types of ore bodies in the area in detail and discusses the characteristics and sources of ore-forming fluids.This provides a new geochemical basis for studying the widespread enrichment of tungsten in the Shimensi tungsten deposit.

Regional Geological Background
Since the Mesoproterozoic era, the Jiangnan orogenic belt has long been under the control of the Yangtze and Cathaysia plates and their junction zones [28,29].After undergoing magmatic-sedimentary-metamorphic-tectonic-metallogenic events as well the Indosinian tectonic movement, the belt was transformed to a continent.During the Yanshanian tectonic movement, due to intraplate contraction and the interaction of the Kula-Pacific Ocean Plates with the neighboring plates, strong continental orogeny occurred in the region, and brought about an NE-NEE-dominated strike slip thrust extensional structure [30,31], which was superimposed and compounded with the ancient structure, bringing up rich nonferrous, rare and precious metal minerals (Figure 1).
The exposed strata of Dahutang ore field are mainly the Mesoproterozoic shallow metamorphic rock series of Anlelin Formation of Shuangqiaoshan Group, with intermediate basic ancient volcanic rock layers inside, which is one of the important metallogenic wall rocks in the area.The fracturing structures in the area are mainly near EW and NE-NNE trending, followed by NW and NS trending.Near EW trending fractures are mainly distributed on the north and south of the ore field, represented by the northern Luositang-Xin'anli fracture, which extends for more than several kilometers and is manifested as strong silicification zone and compression fracture zone.The NE-NNE trending fractures are the most developed within the region and run through the whole area, represented by the NNE trending strike slip thrust belt extensional structure of Wuning-Yifeng-Lianhuashan, which is dominated by compression and torsion.NW trending fractures are not widely distributed in the area, with two relatively larger faults distributed in Guanyintang-Maogongdong-Shimensi areas [10].A large area of intermediate-acid/acid rock is exposed in the area, mainly composed of Jinning stage medium coarse biotite granodiorite, which occurs as a batholith.There are two tectonic magmatic rock belts in Yanshanian rock mass, of which the Huangshaqiao-Meimaoshan magmatic rock belt in the west is the early Yanshanian monzonite, occurring as a batholith, with a total exposed area of more than 100 km 2 [10].Dozens of rock bodies (veins) are distributed in the eastern Jiulongjian Dahutang tectonic magmatic rock belt.From the north to the south, there are five rock bodies: Xin'anli, Shimensi, Dahutang, Shiweidong and Jiulongjian, with an area of 0.3-3 km 2 .Only a small part of the three rock bodies (Shimensi, Dahutang and Shiweidong) are exposed to the surface while most of them are hidden below.The rock body in the shape of a rock stem was formed in the early Yanshanian period.Its zircon U-Pb isotopic age is 136-138 Ma [14].Its lithology is mainly porphyritic biotite granite, fine biotite granite and granite porphyry.It is closely related to the mineralization of tungsten, tin and molybdenum.From north to south, it is manifested as W-Sn-Mo-Cu-Ag (Dahutang ore concentration area) → Mo-W-Cu (Yangshidian) → Cu-Mo-W (Jiulongjian) mineralization zoning (Figure 2).The main deposit types are quartz vein and fine vein disseminated altered granitecomposite deposit type (Shiweidong), quartz vein-fine vein disseminated cryptoexplosive breccia altered granite-composite deposit (Shimensi, Dalingshang), fine vein disseminated altered granite greisen-composite deposit (Dongdouya), quartz vein belt type (Yangshidian), fine vein disseminated altered granite blasting breccia composite deposit (Shiziyan) and tectonic altered rock type (Jiulongjian).All are post magmatic hydrothermal deposits related to Yanshanian rock mass.

Deposit Geology
The Shimensi tungsten deposit has experienced multiple periods and stages of magmatic activity.The biotite granodiorite formed in the late Jinning period is a part of Jiuling batholith.The porphyritic biotite granite (with porphyritic texture, coarse grained phenocrysts, medium fine grained matrix), fine grained biotite granite, and biotite granite porphyry (with porphyritic texture, medium fine phenocrysts and aphanitic matrix), which successively intruded in the middle Yanshan period, have a regularly refined particle size.This shows that their formation depth becomes shallower and their invasion time later, in turn [30].The Yanshanian porphyritic biotite granite is a W-Cu ore-forming rock body, while Jinningian biotite granodiorite is the main ore hosting body (Figure 3).The structure of the ore field is mainly manifested in three ways, namely, ductile shear zone, fracture and joint.According to the strike, it can be divided into four groups: NNE, NEE, NE and NW, with NW trending fractures most developed [32].
The formation of Shimensi deposit takes Yanshanian plutonic to hypabyssal granite body as the metallogenic parent rock, and its genetic type is post magmatic medium-low temperature deposit.The industrial types of ores can be divided into veinlet disseminated type, cryptoexplosive breccia type and quartz vein type, which are characterized by multiple mineralization and "one field, three types" [32].
The vein-disseminated ore bodies are distributed in the main sections of Shimensi deposit and occur in the inner and outer contact zones of porphyritic biotite granite in Mid-Yanshanian period and biotite granodiorite batholith in the Late-Jinning Period.The mode of occurrence is generally gentle, basically consistent with the mode of occurrence of the contact surface, and tends to change with the contact surface (Figure 4).Among them, thick industrial tungsten (copper) ore bodies can be seen in the outer contact zone with a good mineralization continuity.The ore bodies in the inner contact zone are generally thin and poor, with many inclusions and poor mineralization continuity.The alteration type is closely related to the lithology of wall rock.The K-feldsparization in porphyritic biotite granite in the inner contact zone is more obvious, and the biotite granodiorite in the outer contact zone is more obvious.The common alteration processes are greiseni-zation and chloritization.The cryptoexplosive breccia type ore bodies are mainly distributed in the middle of Shimensi ore field.The cryptoexplosive breccia is developed in the Middle Yanshanian fine-grained biotite granite and granite porphyry penetrated along the F20 fracture, mainly in breccia and net vein mineralization.The ore body minerals are various and the ore mineral combination is complex, mainly including wolframite, scheelite, chalcopyrite, molybdenite, etc.The quartz vein type ore bodies are mainly distributed in the range of about 0.6 km 2 in the middle of the ore district, cutting through all rock units and other types of ore bodies in the ore district.

Sampling and Analytical Methods
Three types of ore bodies are basically covered by the samples in this study.Samples are collected from different elevations, drills and tunnels (Table 1).The primary fluid inclusions in these samples are used to represent the ore-forming fluids in different mineralization stages in the Shimensi deposit.A total of 26 inclusion slices were prepared, and the inclusions of quartz in various types of ore bodies were studied in detail.Single minerals were selected for fluid inclusion composition and H-O isotope test.The micro temperature measurement and composition test of inclusions were completed in Beijing Geological Research Institute of China Nuclear Industry.The micro temperature measurement was carried out on the LINKAM THMS600 cold and hot stage, and the working temperature of the cold and hot stage range from 196 °C to 600 °C .The test accuracy is ±0.2 °C below 30 °C and ±2 °C above 30 °C .According to Campbel and Bodnar et al., the salinity is calculated from the melting temperature and freezing point temperature of the cage [33,34].
The vapor composition of fluid inclusion group is extracted by temperature interval explosion to obtain the fluid vapor composition of different mineralization stages.The Prisma TM QMS200 quadrupole mass spectrometer produced by RG202 produced by SHIMADZU in Kyoto, Japan and Anviz in Langenthal, Switzerland was used for the test.The precision of repeated measurement of the instrument is less than 5% [35].
The liquid composition and vapor composition of the fluid inclusion group uses similar means, and the ionic composition of the fluid inclusion group in different temperature ranges is extracted by the same temperature range explosion.The ion chromatograph produced by SHIMADZU in Kyoto, Japan was used to analyze the relative concentration of ions, and the repeated test accuracy of the instrument was less than 5%.The pretreatment of the extracted ion component sample is the same as that of the vapor component sample.
Laser Raman analysis of single fluid inclusion uses a LABHR-VIS LabRAM HR800 produced by HORIBA Jobin Yvon in France research grade micro laser Raman spectrometer, and the excitation wavelength is as follows: λ = 532 nm, Yag crystal frequency doubling, solid laser, laser beam spot ≥1 μm.The scanning time is 10 s to scan twice.
The samples used for H-O isotope analysis are mainly taken from hydrothermal cryptoexplosive breccia and quartz in quartz vein.Pure quartz shall be selected under the microscope with a purity of more than 99%.The method of oxygen isotope analysis is BrF5.First, the pure 12 mg quartz sample is reacted with BrF5 for 15 h, and the oxygen is extracted.The separated oxygen enters the CO2 conversion system with a temperature of 700 °C and a time of 12 min.Finally, CO2 is collected.The hydrogen isotope analysis adopts the explosion method to take water from the quartz inclusion.The test procedure is as follows: heat the quartz inclusion sample, make it explode to release volatile matter, extract water vapor, then make the water react with zinc at 400 °C to produce hydrogen, and then freeze it with liquid nitrogen, and collect it into the sample bottle with active carbon.The stable isotope test was completed in the Test Center of Beijing Geological Research Institute of Nuclear Industry with Finningan MAT251EM mass spectrometer.The international standard for hydrogen and oxygen isotopes is V-SMOW [36], the analysis accuracy of hydrogen isotopes is ±2, and the analysis accuracy of oxygen isotopes is ±0.2.

Distribution and Types of Fluid Inclusions
Microscopic observation shows that there are a large number of fluid inclusions in quartz in different ore bodies in the Shimensi tungsten polymetallic deposit.The types of inclusions are diverse, generally distributed in groups or isolated, and a few secondary inclusions are distributed along fractures.The size of inclusions varies greatly, ranging from 2 μm × 2 μm to 20 μm × 15 μm, but mainly in the range of 6-8 μm.The inclusions are of various shapes, mainly round, oval, strip and irregular.
The objects of this study are all primary inclusions.According to the phase behavior classification criteria of fluid inclusions at room temperature and the phase behavior changes during freezing and rewarming proposed by references [26,37], fluid inclusions can be divided into five types: liquid-rich inclusions (type Ⅰ), vapor-rich inclusions (type Ⅱ), pure liquid inclusions (type Ⅲ), pure vapor inclusions (type Ⅳ), and fluid inclusions containing a solid crystal (type Ⅴ).
Type I: liquid-rich inclusions.This type of inclusion is colorless gray, accounting for more than 70% of the total inclusion, and is the main type of fluid inclusion in the ore field.It is mainly round, oval, strip or irregular, composed of vapor and liquids at room temperature, but dominated by liquid.The vapor-liquid ratio is between 10% and 30%, most of which are between 15% and 20%, with a size of 3 to 20 μm.Most of these inclusions are between 5 and 10 μm (Figure 5A,B).Type II: vapor-rich inclusions.This type of inclusion is few in number, and is composed of vapor and liquids at room temperature, dominated by vapor.The vapor part accounts for more than 50% of the volume of the entire inclusion, and its filling degree is less than 50%.When the homogenization method is used for temperature measurement on the cold and hot stage, the bubble volume gradually increases with the rising temperature, and finally it turns out to be vapor by adapting the homogenization method.The inclusions are usually round, oval or irregular in shape, occurring in isolation, with a size of 2-10 μm.Such inclusions account for about 2% of the total fluid inclusions (Figure 5C).
Type III: pure liquid inclusions.This type of inclusion is less developed, colorless and transparent, mainly round or oval.Only single liquid can be seen at room temperature, and the individual inclusion is generally small (2-6 μm), generally distributed in isolation (Figure 5D).
Type IV: pure vapor inclusions.This type of inclusion is small in number, with a vapor volume content of more than 80%.It is dark gray, isolated, usually irregular and long in shape, with a size of 2-15 μm (Figure 5E).Type V: fluid inclusions containing a solid crystal.This type of inclusion is less developed and consists of vapor, liquid and solid mineral phases.The size of the inclusions is 6-15 μm.Generally, they occur in isolation, in elliptical or irregular shape, and the minerals are divided into opaque blocks and transparent blocks (Figure 5F).

Homogenization Temperature, Salinity and Fluid Density
After detailed microscopic observation of primary fluid inclusions in different mineralization stages, the combination type and abundance of fluid inclusions in different stages are determined, and a certain number of representative fluid-rich fluid inclusions are measured.The results of microthermometry measurement show that the rich fluid inclusions of veinlet-disseminated ore bodies are developed, and the uniform temperature distribution range is wide, from 119 °C to 378 °C , mainly concentrated between 140 °C and 270 °C , with an average of 218 °C .The salinities are from 0.88 wt.% to 9.47 wt.% NaCl equivalent, concentrated between 4.03 wt.% and 8.81 wt.% NaCl equivalent.The densities of fluid are between 0.64 g/cm 3 and 0.99 g/cm 3 .The homogenization temperatures of hydrothermal cryptoexplosive breccia ore body are 139 °C -354 °C , mainly concentrated between 170 °C and 270 °C , with an average of 229 °C .The salinities are distributed between 2.57 wt.% and 8.00 wt.% NaCl equivalent, concentrated between 3.06 wt.% and 6.45 wt.% NaCl equivalent, and the fluid densities are between 0.68 g/cm 3 and 0.96 g/cm 3 .The homogenization temperatures of quartz vein type ore body are between 153 °C and 335 °C , mainly between 190 °C and 270 °C , with an average of 238 °C .The salinities are distributed between 2.90 wt.% and 8.95 wt.% NaCl equivalent, concentrated between 4.0 wt.% and 7.0 wt.% NaCl equivalent, and the fluid densities are between 0.70 g/cm 3 and 0.96 g/cm 3 (Table 2).The salinity and homogenization temperature histograms of the fluid inclusions from different types of ore body are shown in Figure 6, and the homogenization temperature versus salinity diagram is shown in Figure 7.Among them, the densities of the inclusions are calculated by the density calculation formula [38].

Vapor Composition in Fluid Inclusions
Convert the unit of vapor composition analysis data of mineral inclusions μg/g (ω/10 −6 ) into its corresponding molar concentration percentage (mol%).It can be seen from the table that the vapor composition in fluid inclusions in the Shimensi tungsten polymetallic ore is dominated by H2O, which has an absolute advantage (the content is more than 99.9 mol%) and belongs to high water-bearing fluid.In addition, it contains a very small amount of CO2 (0.00005 mol%-0.00057mol%) and trace amounts of N2, H2, CO and CH4.

Laser Raman Spectroscopy Test of Fluid Inclusions
On the basis of micro temperature measurement, the fluid inclusions of quartz in three different types of ore bodies were also analyzed by laser Raman spectroscopy test.The results show that the composition of inclusions in three different types of ore bodies in the ore district is basically the same.The vapor-liquid inclusions show 2916 cm −1 , 2329 cm −1 and 1389 cm −1 spectral peaks, and the 1285 cm −1 peak is relatively strong (Figure 8), indicating that there are vapors of CH4, N2 and CO2.

H-O Isotope Characteristics
The results of hydrogen and oxygen isotope analysis of quartz samples are shown in Table 4.The hydrogen isotope variation of quartz samples is small, ranging from −77.8‰ to −60.6‰ (average value −69.25‰).
The δDV-SMOW values of veinlet-disseminated type ore bodies range from −69.3‰ to −64.2‰, with an average value of −66.93‰.The δ 18 OV-SMOW of values of quartz range from 12.5‰ to 13.5‰, with an average value of 12.95%, and the calculated δ 18 OH 2 O values range from 1.92‰ to 2.91‰, with an average value of 2.37%.
The δDV-SMOW values of hydrothermal cryptoexplosive breccia-type ore bodies range from −68.6‰ to −60.6‰, with an average value of −64.45‰.The δ 18 OV-SMOW of values of quartz range from 12.4‰ to 13.5‰, with an average value of 12.95%, and the calculated δ 18 OH 2 O values range from 1.82‰ to 2.91‰, with an average value of 2.37%.
The δDV-SMOW values of thick quartz vein-type ore bodies range from −77.8‰ to −66.7‰, with an average value of −71.85‰.The δ 18 OV-SMOW of values of quartz range from 12.2‰ to 15.3‰, with an average value of 13.02%, and the calculated δ 18 OH 2 O values range from 2.09‰ to 5.12‰, with an average value of 2.90%.

Ore-Forming Fluid System
The study of fluid inclusions in the Shimensi mine area shows that the inclusions are mainly liquid-rich inclusions, with a small number of vapor-rich inclusions, pure liquid inclusions, pure vapor inclusions and fluid inclusions containing a solid crystal, with formation temperatures concentrated in the range of 140 °C to 270 °C and salinities varying mainly in the range of 3 wt.%-5wt.% NaCl equivalent, which are medium to low temperature, low salinity magmatic fluids.The fluid temperature and salinity ranges of the finevein dipping, quartz vein type and hidden explosion breccia type ore bodies in the area are basically the same (Figure 7), indicating that the fluids of different mineralization types in the area are the same mineralizing fluid.
The anionic composition of the liquid of the inclusion body indicates that the mineralizing fluid is a hot aqueous solution rich in Na + , Ca 2+ , SO4 2− , Cl − and small amounts of F − and NO3 − (Figure 9), which has a strong ability to dissolve mineralized materials [20].Therefore, the ore-forming fluid of Shimensi tungsten ore is mainly a Ca 2+ -Na + -SO4 2− -Cl − type fluid.Mg 2+ and Ca 2+ are the most important alkaline metal ions and prevalent elements in hydrothermal fluids, and the Na + /K + of magmatic hydrothermal fluids is generally less than one, while the fluids from deeper genesis contain higher F − , Cl − and larger F − /Cl − [39,40].The analysis results show that the fluid inclusions have higher Ca 2+ content, indicating that Ca 2+ dominates in the mineralizing solution, which is conducive to the formation of scheelite (CaWO4).Recent studies have shown that the granitic magma may be the source of the ore-forming material [13−14].According to the ion concentration characteristics of the fluid components, the Na + /K + in the fluid inclusions within the deposit are all greater than one, while the F − /Cl − are all less than one.This is mainly due to the characteristics of the mineralizing fluid as a magmatic fluid, and with the mineralization proceeds, atmospheric precipitation enters a larger proportion of the fluid system [20].The low F − content in fluid inclusions within the Shimensi mine indicates that F does not play a dominant role in the mineralization process, or it is reduced due to the decrease of atmospheric precipitation.
The vapor in the inclusions is dominated by CO2 and N2, and contains a small amount of hydrocarbon vapor components.For the transport mechanism of tungsten in mineralizing fluids, some scholars believe that it may migrate as carbonate and bicarbonate under high temperature conditions [41][42][43].However, some others believe that tungsten mainly migrates as multiple complexes [44], and the effect of CO2 on tungsten in hydrothermal transport is minimal.Although the presence of CO2 was detected in the fluid composition analysis, its distribution is relatively small, and the mineralization temperature of the deposit is 200 °C to 300 °C , which is mainly low to medium temperature, and it does not meet the migration conditions of carbonate and bicarbonate under high temperature conditions, indicating that tungsten still mainly migrates as complexes.

Ore-Forming Pressure
Inferring the capture pressure based on fluid inclusion thermometry data is one of the important means to restore the original mineralizing environment [45].The compositional analysis of the mineralizing fluids indicates that the Shimensi tungsten ore-forming fluid system is the NaCl-H2O system.For the low and medium salinity systems, the inclusions equivalence equation was used to calculate the mean pressure [37], which ranged from 83.67 × 10 5 Pa to 365.09 × 10 5 Pa for the Shimensi deposit.The mean pressure of veindisseminated ore body ranges from 150.66 × 10 5 Pa to 365.09 × 10 5 Pa, hydrothermal cryptoexplosive breccia type ore body from 83.67 × 10 5 Pa to 365.09 × 10 5 Pa and quartz vein type ore body from 107.77 × 10 5 Pa to 349.59 × 10 5 Pa.The above three ore bodies show consistent uniform pressure characteristics with a wide range of variation, but the overall pressure values are small and belong to a low-pressure mineralization environment.
This may be caused by two factors.First, as the mineralization works, the mineralization system gradually transforms into an open system, and each mineralization stage contains multiple phases of mineralization fluids, such as the addition of atmospheric precipitation, which leads to a sudden decrease in mineralization pressure.Second, the fluid inclusions in each mineralization stage may be affected by late geological events and the inclusions are damaged, which cannot accurately reflect the mineralization pressure, and the mineralization pressure of some inclusions is lower than 100 × 10 5 Pa, which is obviously inconsistent with the geological phenomenon in the area.

pH Value of Ore-Forming Fluid
pH is a physical quantity that measures the acidity of the mineralizing solution.Generally, it cannot be measured directly from the inclusions, so calculations are made with the help of full analytical data of inclusion composition and chemical equilibrium reactions of mineral coeval combinations.After analysis, the inclusions in the Shimensi mine turn out to belong to the low-salinity NaCl-H2O fluid system, so the concentration approximation can be used instead of activity in the calculation.In addition to the H2O and NaCl components, five main ions are present in the NaCl-H2O hydrothermal fluid: H + , OH − , HCl, Cl − , and Na + , and the total mass balance is reached.Based on the calculating formula by Liu et al. [38], it can be concluded that the pH value of mineralizing fluids in the Shimensi ore district ranges from 5.14 to 5.34 in the corresponding temperature-pressure range.The pH values of the mineralizing fluids do not vary much in each stage and all show a weakly acidic mineralizing condition, which is basically consistent with the mineralizing fluids in regional porphyry tungsten ores [46].

Eh Value of Ore-Forming Fluid
Eh value represents the oxidation-reduction potential of mineralizing solution, which can quantitatively measure the degree of oxidation-reduction of mineralizing solution.The Eh value of the mineralizing fluid is closely related to the pH value of the solution, and its calculation is based on the vapor-liquid equilibrium relationship of fluid inclusions, where the equilibrium vapor compositions include CH4(g), CO2(g) and H2O(g), and the equilibrium liquid is mainly H2O(l).Using the Coulomb reaction formula, it is concluded that the Eh value of the ore-forming fluids in the Shimensi tungsten deposit is between −0.14 and −0.08 within the corresponding temperature-pressure range.It can be seen that the Eh values of its fluids do not vary much, showing a weakly reducing environment.

Genesis of Ore-Forming Fluid
The hydrogen isotope variation of the quartz samples is small, with an average value of −71.24‰, all in the range of normal magmatic water (−80‰ to −50‰) [47], indicating that the original fluid of the mineralizing fluid came from the magma itself.The δ 18 O value of quartz varies between 12.2‰ and 15.3‰, with an average value of 12.99‰, showing the granitic δ 18 O (7‰-13‰) characteristics [48].With reference to the measured fluid inclusion homogeneous temperatures of hydrothermal cryptogenic breccia-type ore bodies and quartz vein-type ore bodies in the Shimensi tungsten polymetallic deposit, the δ 18 O water values of hydrothermal water in equilibrium with quartz in both ore body types were calculated using the quartz-water fractionation equation at 229 °C and 238 °C .By virtue of δ 18 O values of quartz and the homogeneous temperatures of minerals [49], the calculated δ 18 O values of the mineralized fluids vary between 1.82‰ and 5.16‰, with an average value of 2.68‰.The δ 18 O values of mineralizing fluids in hydrothermal cryptoblastic breccia ore bodies vary from 1.82‰ to 2.91‰, with an average value of 2.37‰.The δ 18 O values of mineralizing fluids in quartz vein ore bodies vary from 2.19‰ to 5.16‰, with an average value of 2.90‰, all deviating from the range of magmatic water variation.In the δD-δ 18 OH 2 O water diagram (Figure 10), the sample points fall below and to the left of the magma water, i.e., the transition area between magma water and atmospheric precipitation.It has been shown that many factors can influence the hydrogen and oxygen isotopic compositions of mineralized fluids, such as mineralization temperature, type of water, W/R ratio during water-rock exchange, etc.The lack of hydrogen and oxygen isotope values for the surrounding rocks and the Yanshan-age atmospheric precipitation make it difficult to explore the specific water-rock reaction patterns in the mine area.However, quartz is prone to undergo isotope equilibrium re-exchange reactions with the water it contains, causing the measured oxygen isotopic compositions of the inclusions to not fully reflect the δ 18 OH 2 O water values of the original ore-bearing solution.But quartz contains almost no hydrogen atoms, so the exchange effect on the hydrogen isotopic compositions of the fluid inclusions is minimal [50].The surrounding rocks of the mine are biotite granodiorite, which also contains hydrogen-bearing minerals, but if the water-rock exchange reaction occurs, the hydrogen isotopic composition of the fluid will not change greatly after the exchange, i.e., this change is negligible [51].Therefore, the hydrogen isotopic composition represents the composition of the original solution.δD-δ 18 OH 2 O water values deviate from the normal magma water values, giving rise to an obvious "oxygen drift".The reason for this phenomenon may be the addition of atmospheric precipitation, which caused the drift of oxygen isotopes to atmospheric precipitation.
The H and O isotopic compositions of the Shimensi tungsten polymetallic deposit are mainly derived from magmatic fluids that were mixed with atmospheric precipitation during evolution of fluids, which is consistent with most current studies on the genesis of mineralizing fluids for tungsten deposits [7,43,51].

Mineral Precipitation Mechanism
Previous studies have shown that miscibility of fluids, cooling of the fluid system, boiling of the fluid system, and reduction of fluid pressure are important mechanisms for tungsten precipitation in fluids [41,[52][53][54][55][56].
The mineralizing fluids of tungsten deposits are considered by previous authors to be mainly magmatic water with the mixing of atmospheric precipitation.In this paper, the results of H and O isotope analysis indicate that the mineralizing fluid of Shimensi tungsten deposit is mainly derived from magmatic hydrothermal fluid with the mixing of atmospheric precipitation during the evolution.Therefore, the mixing of different fluids is the main factor of the precipitation of tungsten ore.
In general, boiling of fluids will result in a series of inclusions capturing different liquid and vapor ratios.The mineralizing fluids are apparently multicomponent fluids, in addition to the main equisolutes, the fluid also contains trace equicomponents, often in close symbiosis, showing that they are captured in an inhomogeneous fluid state [57].This phenomenon may be caused by single fluid immiscibility.The petrographic study of fluid inclusions shows that the proportion of "fluid boiling" phenomenon, where liquid-rich inclusions coexist with vapor-rich inclusions and pure vapor inclusions, is low at Shimensi tungsten mine.This indicates that overall, this phenomenon contributes little to mineralization.Therefore boiling is not a major factor in the precipitation of tungsten ore.
Sudden changes in pressure usually lead to boiling, accompanied by changes in fluid inclusion volume and density [58].Fluid inclusion tests show a wide range of pressure variation in the Shimensi mine, but the whole may have been damaged at a later stage and does not reflect the mineralization process.Petrographic observations of fluid inclusions in the area did not reveal significant changes in volume of inclusions, and density variations range from 0.64 g/cm 3 to 0.99 g/cm 3 , indicating that pressure reduction is not the main mechanism leading to mineral precipitation.
The cooling effect of the fluid system is also one of the main mechanisms of tungsten precipitation in fluids [52,[59][60][61].It was shown that the solubility of wolframite, scheelite and other tungsten minerals in the fluids increased significantly when the temperature increased [44], which shows that the effect of temperature change on tungsten precipitation is obvious.The present study shows that fluid inclusion thermometry data in the mine area show a wide range of mineralization temperature variation (119 °C to 378 °C ).Salinity decreases with falling homogeneous temperature.The three different types of ore bodies have common variation characteristics, so the decrease in temperature may also play a role in promoting the precipitation enrichment of the mineral matter in Shimenji ore district.

Conclusions
(1) The fluid inclusions of quartz in the Shimensi tungsten polymetallic deposit are classified into five types, including liquid-rich inclusions (type I), pure liquid inclusions (type II), vapor-rich inclusions (type III), pure vapor inclusions (type IV), and fluid inclusions containing a solid crystal (type V).Microthermometry results show that the homogeneous temperature distribution of fluid inclusions in each mineralization stage is wide, ranging from 119 °C to 378 °C.The salinities vary from 0.88 wt.% to 9.47 wt.% NaCl equivalent, and the densities of the fluid are between 0.64 g/cm 3 and 0.99 g/cm 3 , which belong to low-density fluid.(2) The liquid anion and cation components of fluid inclusions in the region are rich in Na + , Ca 2+ , SO4 2− , Cl − and a small amount of F − and NO3 − .The vapor components are mainly H2O, which belong to NaCl-H2O fluid system containing a small amount of CO2, CH4 and N2.The high Ca 2+ content in the fluid inclusions is favorable to the formation of scheelite (CaWO4).Pa to 349.59 × 10 5 Pa.The variation ranges widely, but the overall pressure value is small, which refers to a low-pressure metallogenic environment.The pH values in the corresponding temperature-pressure range are 5.14 to 5.34, indicating a weakly acidic metallogenic condition, which is basically consistent with the metallogenic fluid of regional porphyry tungsten ore.The Eh values are between −0.14 and −0.08, showing a weakly reducing environment.(4) The H-O isotope test results show that the δDV-SMOW values of the quartz samples range from −77.8‰ to −60.6‰, and the average δ 18 O value of the mineralizing fluid is 2.90‰, indicating that the mineralizing fluid of the Shimensi tungsten polymetallic deposit mainly originates from the magma and has the mixing effect of atmospheric precipitation during the ascending process.The lower temperature may play a certain role in the precipitation enrichment of tungsten elements.The lower temperature may have contributed to the precipitation enrichment of tungsten.

Figure 1 .
Figure 1.Map presenting the distribution of granites and volcanic rocks in South China, and also the location of the Jiangnan Orogen and the Shimensi ore deposit (modified after[25]).

Figure 3 .
Figure 3. Geological map of the Shimensi deposit.

Figure 6 .
Figure 6.Histograms of homogenization temperature and salinity in different ore types.(a) The homogenization temperature of the veinlet-disseminated ore body; (b) the salinity of the veinletdisseminated ore body; (c) the homogenization temperature of the cryptoexplosive breccia ore body; (d) the salinity of the cryptoexplosive breccia ore body; (e) the homogenization temperature of the quartz vein type ore body; (f) the salinity of the quartz vein type ore body.

Figure 7 .
Figure 7. Diagram of the homogenization temperature versus salinity of fluid inclusions in the Shimensi tungsten deposit.

Figure 8 .
Figure 8. Laser Raman spectra of fluid inclusions in the Shimensi tungsten deposit.(a,b) Vapor of type I fluid inclusion of veinlet-disseminated ore body; (c) vapor of type I fluid inclusion of hydrothermal cryptoexplosive breccia ore body; (d) vapor of type I fluid inclusion of thick quartz vein ore body.

Figure 9 .
Figure 9.The trigonometry of the liquid phase composition of the fluid inclusions in the Shimensi tungsten deposit (the base figure was modified after [24]).(a) Cation triangle diagram of fluid inclusions; (b) anion triangle diagram of fluid inclusions.

Table 1 .
Characteristics and locations of samples for microthermometric studies of the Shimensi tungsten deposit.

Table 2 .
Homogenization temperatures, salinities and densities of the primary fluid inclusions in the Shimensi tungsten deposit.

Table 3 .
Compositions of the fluid inclusions in the Shimensi tungsten deposit (units: μg/g).

Table 4 .
Hydrogen and oxygen isotope compositions of the Shimensi tungsten deposit.