CN115078520A - Mineral geochemistry-based porphyry system mineralization evaluation method - Google Patents

Mineral geochemistry-based porphyry system mineralization evaluation method Download PDF

Info

Publication number
CN115078520A
CN115078520A CN202210664417.2A CN202210664417A CN115078520A CN 115078520 A CN115078520 A CN 115078520A CN 202210664417 A CN202210664417 A CN 202210664417A CN 115078520 A CN115078520 A CN 115078520A
Authority
CN
China
Prior art keywords
rock mass
rock
data
lean
sio
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210664417.2A
Other languages
Chinese (zh)
Other versions
CN115078520B (en
Inventor
吴松
吴健辉
林毅斌
林德才
易建洲
郑有业
刘鹏
沙仙武
豆孝芳
陈才贤
曾庆淋
张亚明
李健昌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Geosciences Beijing
Tibet Julong Copper Co Ltd
Original Assignee
China University of Geosciences Beijing
Tibet Julong Copper Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Geosciences Beijing, Tibet Julong Copper Co Ltd filed Critical China University of Geosciences Beijing
Priority to CN202210664417.2A priority Critical patent/CN115078520B/en
Publication of CN115078520A publication Critical patent/CN115078520A/en
Priority to NL2033632A priority patent/NL2033632A/en
Application granted granted Critical
Publication of CN115078520B publication Critical patent/CN115078520B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24Earth materials
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/62Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode
    • G01N27/64Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode using wave or particle radiation to ionise a gas, e.g. in an ionisation chamber

Abstract

The invention discloses a porphyry system mineralizing evaluation method based on mineral geochemistry, which indicates the mineralization potential of a rock mass by using the total rock geochemistry and the mineral chemical parameters. The method specifically comprises the following steps: the method is applied to the field of porphyry system exploration, has the advantages of short testing time, low cost, convenience and quickness, can effectively distinguish the mineralized rock mass from the lean mineralized rock mass in the porphyry system and shorten the mineral exploration period, is a novel indispensable exploration means and method, and has important popularization value.

Description

Mineral geochemistry-based porphyry system mineralization evaluation method
Technical Field
The invention belongs to the field of mineral exploration and evaluation, and particularly relates to a method for evaluating the mineralocorticity of a porphyry system based on mineral geochemistry.
Background
Porphyry deposits are characterized by typical large tonnage, low grade and large scale hydrothermal alteration and enrichment of metal sulfides. The formation of large porphyry deposits requires that the magma have high oxygen fugacity, water content and volatile elements such as S, Cl. The magma water has important control function on the migration of metal elements and the precipitation of ore-containing hot liquid, and the deep metal is released from a sulfide phase by high oxygen fugacity and is brought to a shallow ore-forming system. Therefore, finding out the magma water and oxygen fugacity characteristics of rock mass is key to the determination of the mineralised rock mass. In the process of forming the porphyry deposit, high water content and oxygen fugacity of the rock pulp can be recorded in the processes of rock pulp formation and mineral separation and crystallization, so that the rock mass and lean rock mass can be effectively distinguished by using the geochemical and mineralogical chemical components (zircon and amphibole) of the whole rock.
The traditional method for searching the region with porphyry mineralization potential needs to carry out a large proportion of geochemical survey, geophysical survey, regional geological survey and systematic test and comprehensive research, and has the following defects: the early period of the exploration evaluation is long, the cost is high, the accuracy is low, and the urgent need of rapid exploration evaluation cannot be met.
Aiming at the problems, the invention provides a quantitative index for rapidly identifying the ore-forming rock mass in the porphyry ore-forming system through a large amount of experimental research and ore-finding practice, solves the technical problem of searching porphyry ore deposit under large scale, and realizes the organic combination of mineral geochemistry and mineral-containing evaluation.
Disclosure of Invention
The invention aims to provide a method for evaluating the mineralization of a porphyry system based on mineral geochemistry, which has the advantages of short test time, low cost, convenience and rapidness, can quickly evaluate the mineralization potential of the porphyry under a large scale, effectively distinguishes the mineralization rock mass and lean rock mass in the porphyry system, and shortens the mineral exploration period.
In order to solve the technical problem, the invention adopts the following technical scheme:
the method for evaluating the mineralization of the porphyry system based on the mineral geochemistry comprises the following steps:
1) the target area of the ore is determined according to geological, geophysical prospecting, chemical prospecting and remote sensing data in the collecting area of the system;
2) the system collects all the acidic rock masses in the target area and describes the lithology, alteration and mineralization characteristics of each type of sample;
3) selecting a sample without alteration or weak alteration, and carrying out chemical analysis to obtain the total rock main trace element SiO 2 The contents of Sr, Y, V and Sc are represented as c (SiO) 2 ) C (Sr), c (Y), c (V), c (Sc); the contents of trace elements Ce, Nd, Y, Eu, Sm and Gd of zircon are marked as c (Ce), c (Nd), c (Y'), c (Eu), c (Sm) and c (Gd); the temperature T of the amphibole and the oxygen fugacity Delta FMQ;
4) taking F1-F4 as discrimination factors to discriminate the mineralized rock mass from the lean rock mass, specifically:
discrimination factor F1: f1 ═ 0.925 × c (SiO 2) 2 ) +113.75 (equation 1);
discrimination factor F2: f2 ═ 0.355 × c (SiO) 2 ) +34.15 (equation 2);
discrimination factor F3: f3 ═ 41.7 ═ c (ce)/c (nd)/c (Y') +4.707 (formula 3);
discrimination factor F4: f4 ═ 0.0025 ═ c (t) +3.8 (formula 4);
the total rock main trace element SiO obtained in the step 3) 2 Substituting the content into the above formula 1, calculating the discrimination factor F1 when c (Sr)/c (Y)>F1, judging the rock mass to be an ore-forming rock mass, otherwise, judging the rock mass to be a lean rock mass;
the total rock main trace element SiO obtained in the step 3) 2 Substituting the content into the above formula 2, calculating the discrimination factor F2 when c (V)/c (Sc)>F2, judging the rock mass to be an ore-forming rock mass, otherwise, judging the rock mass to be a lean rock mass;
substituting the ratio of zircon (c (Ce)/c (Nd))/c (Y ') obtained in step 3 into the above formula 3, calculating the discrimination factor F3 when 10000 (c (Eu)/c (Eu))/c (Y'))>F3, judging to be an ore-forming rock mass, otherwise, judging to be a lean rock mass, wherein
Figure BDA0003691065860000021
Substituting the amphibole temperature T obtained in the step 3) into the formula 4, and calculating a discrimination factor F4, wherein when the delta FMQ is greater than F4, the amphibole is discriminated as an ore-forming rock mass, and otherwise, the amphibole is discriminated as a lean rock mass;
when the four discrimination factors are all discriminated to obtain an ore-forming rock mass, the ore-forming rock mass is determined; and when at least one discrimination factor is used for discriminating the lean rock mass, the lean rock mass is determined.
When the individual sample is equal to a certain discrimination factor, comprehensive judgment and error correction can be carried out by referring to the results of other three discrimination factors, and under the condition of ensuring that the sample amount is sufficient, the discrimination of the ore-forming rock mass and the lean rock mass obtained on the basis of the statistical rule is credible and accurate.
According to the scheme, in the step 3), a sample without change or weak change is selected, and the specific steps of chemical analysis are as follows: and grinding the collected sample into powder, and grinding the powder into a probe sheet and a laser in-situ target, wherein the powder sample is used for testing main and trace components of the whole rock, and the probe sheet and the laser in-situ target are used for testing main and trace components of single mineral.
Preferably, the monominerals are zircon and amphibole.
More preferably, the sample is ground into powder, ground into a probe sheet and a laser in-situ target; wherein: the rock powder is used for analyzing the main and trace components of the whole rock to obtain the main and trace elements SiO of the whole rock 2 Sr, Y, V, Sc content; selecting hornblende from the sample, analyzing the components of the hornblende, determining the chemical components and types of the hornblende, and obtaining the temperature and oxygen loss of the hornblende, wherein the temperature of the hornblende is T (° C) (-151.487 × Si) +2041, and Si ═ Si + (Al) IV /15)-(2×Ti IV )-(Al VI /2)-(Ti VI /1.8)+(Fe 3+ /9)+(Fe 2+ /3.3)+(Mg/26)+( B Ca/5)+( B Na/1.3)-( A Na/15)+((1- A Na- A K)/2.3),Al IV And Ti IV Is the atomic number of Al and Ti in tetrahedron, Al VI And Ti VI Is the atomic number of Al and Ti in an octahedron, A na is the Na content of the A site, A k is the content of K at the A site, B na is the Na content of the B site, B ca is the Ca content of the B site, and the oxygen fugacity is represented by log fO 2 Obtained at-24441.9/t (k) +8.290(± 0.167); and (3) selecting a representative mineral (zircon) to perform laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) in-situ micro-area element analysis on the sample ground into the laser in-situ target, and obtaining the contents of trace elements Ce, Nd, Y, Eu, Sm and Gd of the zircon.
According to the above scheme, in the step 3), for a single mineral, in the analysis process, a small-particle mineral inclusion or inherited mineral is generally encountered, and in order to obtain more accurate data, the data needs to be preprocessed and interpreted, which is divided into the following three steps: firstly, importing the obtained data in the original csv format into ICPMSDataCal software, and eliminating the data which are hit to an inclusion or are hit by minerals according to the abnormality of an element integration curve of each analysis test point; secondly, regarding the eliminated data, the mixed dyeing data are further eliminated by using the following standards: la > 1.5ppm is regarded as mixed-dyeing apatite, Fe>5500ppm reacted the mixed dyeing of Fe oxide, Ti > 60ppm regarded as mixed dyeing of Ti oxide, Ba>10ppm was taken as a mixed dyeing of the fluid inclusions; (iii) use 206 Pb/ 238 U isotope age values exclude interference from inherited minerals, i.e. when data points are obtained 206 Pb/ 238 Weighted average of U isotope ages less than or greater than sample 206 Pb/ 238 At the age of the U isotope, the data point is considered invalid data; after the data are processed by the three steps, the rest data can be used for the next mining character evaluation.
According to the scheme, in the step 4), the steps of obtaining the discrimination factors F1-F4 are as follows:
1) respectively selecting more than 60 ore-forming rock mass samples without or with weak alteration and more than 60 lean rock mass samples;
2) carrying out chemical analysis on the sample selected in the step 1) to obtain the total rock main trace element SiO 2 The contents of Sr, Y, V and Sc are represented as c (SiO) 2 ) C (Sr), c (Y), c (V), c (Sc); the contents of trace elements Ce, Nd, Y, Eu, Sm and Gd of zircon are marked as c (Ce), c (Nd), c (Y'), c (Eu), c (Sm) and c (Gd); the temperature T of the amphibole and the oxygen fugacity Delta FMQ;
3) drawing a scatter diagram based on the data of a large number of mining areas developed ore-forming rock masses and lean rock masses obtained in the step 2), fitting a straight line according to a plurality of data points near the boundary point of the ore-forming rock masses and the lean rock masses, and solving the equation of the straight line as a discrimination factor, wherein the equation specifically comprises the following steps:
obtaining discrimination factors F1-F4 through fitting and calculation, specifically:
computing a discriminant factor F1
C (SiO) is used as the main trace element of the whole rock obtained in the step 2) 2 ) Projecting the diagram by using the abscissa c (Sr)/c (Y) as the ordinate, obtaining and fitting a boundary line between the mineralized rock mass and the lean rock mass according to the projection range, and calculating by using the following formula to obtain a discrimination factor F1: f1 ═ 0.925 × c (SiO 2) 2 ) +113.75 (equation 1);
② calculating a discrimination factor F2
To the product obtained in step 2)Total major trace elements of rock, with c (SiO) 2 ) C (V)/c (Sc) is used as an abscissa for projection, a boundary line of the mineralized rock mass and the lean rock mass is obtained according to the projection range and is fitted, and a discrimination factor F2 is obtained by utilizing the following formula: f2 ═ 0.355 × c (SiO) 2 ) +34.15 (equation 2);
computing discrimination factor F3
Projecting the zircon microelements obtained in the step 2) by taking (c (Ce)/c (Nd)/c (Y ') as abscissa and 10000 (c) (Eu)/c (Eu))/c (Y') as ordinate, wherein
Figure BDA0003691065860000041
Obtaining and fitting a boundary line of the mineralized rock mass and the lean rock mass according to the projection range, and calculating by using the following formula to obtain a discrimination factor F3: f3 ═ 41.7 ═ c (ce)/c (nd)/c (Y') +4.707 (formula 3);
fourthly, calculating a discrimination factor F4
Calculating the formation temperature and oxygen fugacity of the hornblende component obtained in the step 2), carrying out projection by taking T (temperature) as a horizontal ordinate and Delta FMQ (oxygen fugacity) as a vertical coordinate, obtaining a boundary line between an ore-forming rock mass and a lean rock mass according to the projection range, fitting, and calculating by using the following formula to obtain a discrimination factor F4: f4 ═ 0.0025 × T +3.8 (formula 4).
Preferably, in the step 2), for a single mineral, during the analysis process, a small-particle mineral inclusion or inherited mineral is generally encountered, and in order to obtain more accurate data, the data needs to be preprocessed and interpreted, which is divided into the following three steps: firstly, importing the obtained data in the original csv format into ICPMSDataCal software, and eliminating the data which are hit to an inclusion or are hit by minerals according to the abnormality of an element integration curve of each analysis test point; secondly, regarding the eliminated data, further eliminating mixed dyeing data by using the following standards: la > 1.5ppm is regarded as mixed-dyeing apatite, Fe>5500ppm reacted the mixed dyeing of Fe oxide, Ti > 60ppm regarded as mixed dyeing of Ti oxide, Ba>10ppm was taken as a mixed dyeing of the fluid inclusions; (iii) use 206 Pb/ 238 U isotope age values exclude interference from inherited minerals, i.e. when data points are obtained 206 Pb/ 238 Weighted average of U isotope ages less than or greater than sample 206 Pb/ 238 At the age of the U isotope, the data point is considered invalid data; after the data are processed by the three steps, the rest data can be used for the next mining character evaluation.
The invention provides a novel method for rapidly evaluating the mineralization potential of regional porphyry under a large scale, which utilizes the geochemistry and mineralogical chemical parameters of the whole rock to indicate the mineralization potential of a rock mass. The principle of obtaining the discrimination factors F1-F4 is as follows: for the porphyry deposit, the key to the mineralization is the high oxygen fugacity and water bearing property of the magma, for example in the case of H-rich ore deposits 2 Magma system of O (>6%H 2 O), first crystallizing amphibole, then plagioclase and finally magnetite, resulting in increased Sr and V elements in the residual melt and near uniform Eu/Eu values in the evolving aqueous melt due to the inhibition of early plagioclase and magnetite crystallization; the crystallization of amphibole leads to the reduction of Sc element and Y element, and finally leads to the rock mass with mineralization potential to have the characteristics of high whole rock c (Sr)/c (Y) and c (V)/c (Sc) ratio, high zircon 10000 (c) (Eu)/c (Eu))/c (Y '), (c (Ce)/c (Nd))/c (Y')) ratio and high amphibole oxygen escape degree. According to the characteristics, a scatter diagram is drawn based on data of a large number of mining areas developing diagenetic rocks and lean rocks, a straight line is fitted according to a plurality of data points near a boundary point of the diagenetic rocks and the lean rocks, an equation of the straight line is obtained as a discrimination factor, and 4 discrimination factors F1-F4 are obtained.
Although the relevant parameters of the discrimination factor may have certain errors, under the premise of comprehensively considering the mineralization characteristics of the porphyry deposit (the chemical properties of the mineralized rock mass and the lean rock mass of different porphyry deposits have no significant change), and under the condition of simultaneously predicting large sample data points (data points >10), based on the statistical law, the discrimination factor provided by the invention can limit the range of the lean rock mass and the mineralized rock mass, and can meet the rapid evaluation of the mineralization of the development rock mass in the regional scale.
The invention has the following beneficial effects:
the invention provides a method for evaluating the mineralization of a porphyry system based on mineral geochemistry, which is used for quantitatively identifying a discrimination factor of an mineralized rock mass in the porphyry system by utilizing data of whole-rock geochemistry, mineralogy and the like according to the mineralization characteristics of porphyry deposits.
Drawings
FIG. 1 is a scatter diagram and a fitting straight line drawn according to data of an ore-forming rock body and an ore-poor rock body in the embodiment of the invention, wherein (a) and (b) are main trace elements of a whole rock; FIG. (c) shows zircon trace element; FIG. d shows the scintillation thermometer and oxygen fugacity, and the equations in the graphs of FIGS. a-d correspond to discriminant factors F1-F4, respectively.
Fig. 2 is a technical flow chart for identifying an ore-forming rock mass in a porphyry ore-forming system in the embodiment of the invention.
Detailed Description
The invention is further illustrated by the following specific examples.
Example 1
Obtaining discrimination factors F1-F4, which comprises the following steps:
1) respectively selecting 61 pieces of mineralized rock mass samples and 66 pieces of lean rock mass samples which are not altered or weakly altered;
2) grinding the sample selected in the step 1) into powder, grinding into a probe sheet and laser in-situ target. The rock powder is used for analyzing the main and trace components of the whole rock. And (3) selecting the magma amphibole under a microscope for the sample ground into the probe piece, carrying out electronic probe component analysis, and further determining the chemical components and types of the amphibole. And selecting a representative mineral (zircon) for carrying out laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) in-situ micro-area element analysis on the sample ground into the laser in-situ target.
3) For single mineral, in the in-situ LA-ICP-MS analysis process, small-particle mineral inclusion or inherited mineral is generally encountered, and in order to obtain more accurate data, the data needs to be processedThe pretreatment and the interpretation are divided into the following three steps: firstly, importing the obtained data in the original csv format into ICPMSDataCal software, and eliminating the data which are hit to an inclusion or are hit by minerals according to the abnormality of an element integration curve of each analysis test point; secondly, regarding the eliminated data, further eliminating mixed dyeing data by using the following standards: la > 1.5ppm is regarded as mixed-dyeing apatite, Fe>5500ppm reacted the mixed dyeing of Fe oxide, Ti > 60ppm regarded as mixed dyeing of Ti oxide, Ba>10ppm was taken as a mixed dyeing of the fluid inclusions; (iii) use 206 Pb/ 238 U isotope age values exclude interference from inherited minerals, i.e. when data points are obtained 206 Pb/ 238 Weighted average of U isotope ages less than or greater than sample 206 Pb/ 238 At the age of the U isotope, the data point is considered to be invalid data. After the data are processed by the three steps, the rest data can be used for the next mining character evaluation.
4) Processing the data obtained in the step 3) by using Excel to obtain discrimination factors F1-F4. Wherein the amphibole temperature is determined by T (° C) (-151.487 × Si) +2041, and the oxygen fugacity is determined by log fO 2 Obtained at-24441.9/t (k) +8.290(± 0.167). Defining to obtain total rock main trace element SiO 2 The contents of Sr, Y, V and Sc are c (SiO) 2 ) C (Sr), c (Y), c (V), c (Sc), and the trace elements Ce, Nd, Y, Eu, Sm and Gd contained in the zircon in the amounts of c (Ce), c (Nd), c (Y'), c (Eu), c (Sm) and c (Gd). The method specifically comprises the following steps:
on the premise of collecting data of a large number of developed mineralised and lean rock masses in a mining area, drawing a scatter diagram shown in figure 1, fitting a plurality of data points near a demarcation point of the mineralised rock mass and the lean rock mass into a straight line, and solving an equation of the straight line as a discrimination factor, wherein:
calculating discrimination factor F1
C (SiO) is used as the main trace element of the whole rock obtained in the step 2) 2 ) Projecting the diagram by using the abscissa c (Sr)/c (Y) as the ordinate, obtaining and fitting a boundary line between the mineralized rock mass and the lean rock mass according to the projection range, and calculating by using the following formula to obtain a discrimination factor F1: f1 ═ 0.925 × c (SiO 2) 2 ) +113.75 (equation 1);
② calculating a discrimination factor F2
C (SiO) is used as the main trace element of the whole rock obtained in the step 2) 2 ) C (V)/c (Sc) is used as an abscissa for projection, a boundary line of the mineralized rock mass and the lean rock mass is obtained according to the projection range and is fitted, and a discrimination factor F2 is obtained by utilizing the following formula: f2 ═ 0.355 × c (SiO) 2 ) +34.15 (equation 2);
computing discrimination factor F3
Projection is carried out on the zircon trace elements obtained in the step 2 by taking (c (Ce)/c (Nd)/c (Y ') as abscissa and 10000 (c) (Eu)/c (Eu))/c (Y') as ordinate, a boundary line between the formed rock mass and the lean rock mass is obtained according to the projection range, fitting is carried out, and a discrimination factor F3 is obtained by utilizing the following formula: f3 ═ 41.7 ═ c (ce)/c (nd)/c (Y') +4.707 (formula 3);
fourthly, calculating a discriminant factor F4
Calculating the formation temperature and oxygen fugacity of the hornblende component obtained in the step 2), carrying out projection by taking T (temperature) as a horizontal ordinate and Delta FMQ (oxygen fugacity) as a vertical coordinate, obtaining a boundary line between an ore-forming rock mass and a lean rock mass according to the projection range, fitting, and calculating by using the following formula to obtain a discrimination factor F4: f4 ═ 0.0025 × T +3.8 (formula 4).
5) Method for judging mineralized rock mass and lean rock mass
The obtained total rock major element SiO 2 Substituting the content into the above formula 1, calculating the discrimination factor F1 when c (Sr)/c (Y)>And F1, judging the rock mass to be an ore-forming rock mass, and otherwise, judging the rock mass to be a lean rock mass.
The obtained total rock major element SiO 2 Substituting the content into the above formula 2, calculating the discrimination factor F2 when c (V)/c (Sc)>And F2, judging the rock mass to be an ore-forming rock mass, and otherwise, judging the rock mass to be a lean rock mass.
Substituting the obtained ratio of zircon (c (Ce)/c (Nd))/c (Y ') into the above formula 3, calculating the discrimination factor F3 when 10000 (c (Eu)/c (Eu))/c (Y'))>F3, judging the rock mass to be an ore-forming rock mass, otherwise, judging the rock mass to be a lean rock mass,
Figure BDA0003691065860000071
substituting the obtained amphibole temperature T into the formula 4, calculating a discrimination factor F4, and when c (delta FMQ) > F4, discriminating the amphibole rock mass as an ore-forming rock mass, otherwise, discriminating the amphibole rock mass as a lean rock mass.
When the four discrimination factors are all discriminated to obtain an ore-forming rock mass, the ore-forming rock mass is determined; and when at least one discrimination factor is used for discriminating the lean rock mass, the lean rock mass is determined.
When the individual sample is equal to a certain discrimination factor, comprehensive judgment and error correction can be carried out by referring to the results of other three discrimination factors, and under the condition of ensuring that the sample amount is sufficient, the discrimination of the ore-forming rock mass and the lean rock mass obtained on the basis of the statistical rule is credible and accurate.
Example 2
The method provides the mineralogical evaluation of the whole rock geochemistry and mineralogical chemistry porphyry system: taking a Cu deposit of cinnabar rock as an example, as shown in fig. 2, the method specifically comprises the following steps:
a. through geological mapping and drilling record, the acidic rock mass in the Zhunuo mining area is identified to comprise quartz porphyry, second-Chang granite and granite porphyry.
b. Collecting field samples: collecting weak-alteration or non-alteration medium-acidity rock mass in the Zhunuo earth surface and a borehole. During sampling, the following information was recorded in real detail, as shown in table 1:
TABLE 1 cinnunor mining area sampling record table
Sample number- X Y Lithology Alteration of hand specimen Mineralization of minerals Location of a site
ZN1501 527394 3268370 Quartz porphyry Weak serite petrochemistry Is free of Zhu Nuo
ZN1502 527359 3267072 Two long granite No alteration Is free of Zhu Nuo
ZN1503 526129 3266961 Granite porphyry Weak potassium long petrified stone Small amount of Brass mineralization Zhu Nuo
c. And (3) sample testing: and grinding the collected sample into powder or grinding the collected sample into a probe sheet and a laser in-situ target, wherein the powder sample is used for testing main and trace components of the whole rock, the probe sheet and the laser in-situ target are used for testing main and trace components of a single mineral, and the specific steps refer to step 2) in example 1.
d. Data processing: data processing is carried out by using ICPMSDataCal software, including data importing; (same as step 3 in example 1)); and thirdly, screening data. Finally obtaining the total rock main trace element SiO 2 The contents of Sr, Y, V and Sc are c (SiO) 2 ) C (Sr), c (Y), c (V), c (Sc), trace zircon elements Ce, Nd, Y, Eu, Sm, Gd in c (Ce), c (Nd), c (Y'), c (Eu), c (Sm), c (Gd), and goniamphibole temperature and oxygen loss, wherein the goniamphibole temperature is obtained from T (DEG C) (-151.487 × Si) +2041, and the oxygen loss is obtained from log fO 2 Obtained at-24441.9/t (k) +8.290(± 0.167).
e. And (3) mineral content evaluation: by using the final data (tables 2, 3 and 4) after Excel processing, judging that the collected Honunou di granite (porphyry) is an ore-forming rock mass and is matched with the actual ore-forming rock mass of the Honunou ore deposit according to the calculation results of the judgment factors F1, F2, F3 and F4, and further proving the effectiveness of the new method for evaluating the mineralization of the porphyry system based on mineral geochemistry (figure 1).
Table 2 partial whole rock elemental data results for cinnabar mine
Lithology Sample number SiO2 Sr Y V Sc
Two long granite 14-D01 70.12 604.01 8.75995 43.4805 4.78
Two long granite 14-D06 69.22 695.495 6.84855 67.4415 6.78
Two long granite ZX2-8 67.29 598.31 6.59775 72.3975 6.84
Two long granite 802-105.8 68.48 648.2 9.065 68.38 6.44
Two long granite 802-106.3 68.55 644.5 10.46 78.1 7.19
Two long granite 802-249.4 67.42 631.9 10.09 70.8 6.48
Granite of Erchang granite 702-276.2 71.05 375.3 4.863 44.6 3.64
Granite of Erchang granite 702-278.5 72.17 351.4 4.733 44.39 3.75
Granite of Erchang granite 702-336.5 71.21 420.9 6.314 46.41 3.73
Granite of Erchang granite 005-330.1 69.83 424.27 8.3304 46.08 4.75
Granite of Erchang granite 005-444.7 70.00 446.16 7.7058 45.25 3.77
Granite of Erchang granite 806-202.9 67.92 547.47 10.341 57.62 5.36
Two long granite 11-13 69.11 729.6 6.3483 48.95 3.92
Two long granite 11-14 69.62 748.03 6.19185 53.775 4.19
Two long granite 11-15 69.44 818.33 6.96465 58.375 4.96
Two long granite 11-16 69.29 743.755 6.05325 51.3875 4.24
Table 3 results of partial zircon data for cinnamyl mine area
Figure BDA0003691065860000081
Figure BDA0003691065860000091
Table 4 results of the angulometer composition data for the cinnunor mine area
Figure BDA0003691065860000092

Claims (6)

1. A porphyry system mineralization evaluation method based on mineral geochemistry is characterized by comprising the following steps:
1) delineating an ore-finding target area;
2) collecting all the medium-acid rock masses in the target area, selecting a sample without alteration or weak alteration, and carrying out chemical analysis to obtain the total-rock major trace element SiO 2 The contents of Sr, Y, V and Sc are denoted as c (SiO) 2 ) C (Sr), c (Y), c (V), c (Sc); the contents of the trace elements Ce, Nd, Y, Eu, Sm and Gd in zircon are marked as c (Ce), c (Nd), c (Y'), c (Eu), c (Sm) and c (Gd); the temperature T of the amphibole and the oxygen fugacity Delta FMQ;
3) taking F1-F4 as discrimination factors to discriminate the mineralized rock mass from the lean rock mass, specifically:
discrimination factor F1: f1 ═ 0.925 × c (SiO 2) 2 ) +113.75 (equation 1);
discrimination factor F2: f2 ═ 0.355 × c (SiO) 2 ) +34.15 (equation 2);
discrimination factor F3: f3 ═ 41.7 ═ c (ce)/c (nd)/c (Y') +4.707 (formula 3);
discrimination factor F4: f4 ═ 0.0025 ═ c (t) +3.8 (formula 4);
the total rock major trace element SiO obtained in the step 2) 2 Substituting the content into the above formula 1, calculating the discrimination factor F1 when c (Sr)/c (Y)>F1, judging the rock mass to be an ore-forming rock mass, otherwise, judging the rock mass to be a lean rock mass;
using the total rock major trace element SiO obtained in the step 2) 2 Substituting the content into the above formula 2, calculating the discrimination factor F2 when c (V)/c (Sc)>F2, judging the rock mass to be an ore-forming rock mass, otherwise, judging the rock mass to be a lean rock mass;
the zircon (c), (Ce)/c (c) (of)Substituting the Nd))/c (Y ') ratio into the above formula 3, calculating the discrimination factor F3 when 10000 (c (Eu)/c (Eu))/c (Y'))>F3, judging to be an ore-forming rock mass, otherwise, judging to be a lean rock mass, wherein
Figure FDA0003691065850000011
Substituting the amphibole temperature T obtained in the step 2) into the formula 4, and calculating a discrimination factor F4, wherein when the delta FMQ is greater than F4, the amphibole is discriminated as an ore-forming rock mass, and otherwise, the amphibole is discriminated as a lean rock mass;
when the four discrimination factors are all discriminated to obtain an ore-forming rock mass, the ore-forming rock mass is determined; and when at least one discrimination factor is used for discriminating the lean rock mass, the lean rock mass is determined.
2. The method according to claim 1, wherein in step 2), the sample without alteration or with weak alteration is selected, and the specific steps of performing chemical analysis are as follows: grinding the collected sample into powder, and grinding the powder into a probe sheet and a laser in-situ target, wherein the powder sample is used for testing main and trace components of the whole rock, and the probe sheet and the laser in-situ target are used for testing main and trace components of single minerals, wherein the single minerals are zircon and amphibole.
3. The method of claim 2, wherein the sample is ground into powder, ground into a probe tile and a laser in-situ target; wherein: the rock powder is used for analyzing the main and trace components of the whole rock to obtain the main and trace elements SiO of the whole rock 2 Sr, Y, V, Sc content; selecting magma amphibole from the sample ground into the probe piece, carrying out electronic probe component analysis, further determining the chemical components and types of the amphibole, and obtaining the temperature and the oxygen escape degree of the amphibole, wherein the temperature of the amphibole is obtained from T (° C) (-151.487 xSi) +2041, and the oxygen escape degree is obtained from log fO 2 Obtained at-24441.9/t (k) +8.290(± 0.167); and (3) selecting representative minerals to carry out laser ablation inductively coupled plasma mass spectrum in-situ micro-area element analysis on the sample ground into the laser in-situ target, thereby obtaining the contents of zircon trace elements Ce, Nd, Y, Eu, Sm and Gd.
4. The method according to claim 1, wherein in the step 2), for the single mineral, the preprocessing and interpretation of the data are divided into the following three steps: firstly, importing the obtained data in the original csv format into ICPMSDataCal software, and eliminating the data which are hit to an inclusion or are hit by minerals according to the abnormality of an element integration curve of each analysis test point; secondly, regarding the eliminated data, further eliminating mixed dyeing data by using the following standards: la > 1.5ppm is regarded as mixed-dyeing apatite, Fe>5500ppm reacted the mixed dyeing of Fe oxide, Ti > 60ppm regarded as mixed dyeing of Ti oxide, Ba>10ppm was taken as a mixed dyeing of the fluid inclusions; (iii) use 206 Pb/ 238 U isotope age values exclude interference from inherited minerals, i.e. when data points are obtained 206 Pb/ 238 Weighted average of U isotope ages less than or greater than sample 206 Pb/ 238 At the age of the U isotope, the data point is considered to be invalid data.
5. The method as claimed in claim 1, wherein in the step 3), the steps of obtaining the discriminant factors F1-F4 are as follows:
(1) respectively selecting more than 60 ore-forming rock mass samples without or with weak alteration and more than 60 lean rock mass samples;
(2) carrying out chemical analysis on the sample selected in the step (1) to obtain the total rock main trace element SiO 2 The contents of Sr, Y, V and Sc are represented as c (SiO) 2 ) C (Sr), c (Y), c (V), c (Sc); the contents of the trace elements Ce, Nd, Y, Eu, Sm and Gd in zircon are marked as c (Ce), c (Nd), c (Y'), c (Eu), c (Sm) and c (Gd); the temperature T of the amphibole and the oxygen fugacity Delta FMQ;
(3) drawing a scatter diagram based on the data of a large number of mining areas developed ore-forming rock masses and lean rock masses obtained in the step (2), fitting a plurality of data points near the demarcation point of the ore-forming rock masses and the lean rock masses into a straight line, and solving the equation of the straight line as a discrimination factor, wherein the equation specifically comprises the following steps:
obtaining discrimination factors F1-F4 through fitting and calculation, wherein the discrimination factors are as follows:
calculating discrimination factor F1
C (SiO) is used as the main trace element of the whole rock obtained in the step (2) 2 ) Projecting the diagram by using the abscissa c (Sr)/c (Y) as the ordinate, obtaining and fitting a boundary line between the mineralized rock mass and the lean rock mass according to the projection range, and calculating by using the following formula to obtain a discrimination factor F1: f1 ═ 0.925 × c (SiO 2) 2 ) +113.75 (equation 1);
② calculating a discrimination factor F2
C (SiO) is used as the main trace element of the whole rock obtained in the step (2) 2 ) C (V)/c (Sc) is used as an abscissa for projection, a boundary line of the mineralized rock mass and the lean rock mass is obtained according to the projection range and is fitted, and a discrimination factor F2 is obtained by utilizing the following formula: f2 ═ 0.355 × c (SiO) 2 ) +34.15 (equation 2);
computing discrimination factor F3
Projecting the zircon trace elements obtained in the step (2) by taking (c (Ce)/c (Nd)/c (Y ') as abscissa and 10000 (c) (Eu)/c (Eu))/c (Y') as ordinate, wherein
Figure FDA0003691065850000031
Obtaining and fitting a boundary line of the mineralized rock mass and the lean rock mass according to the projection range, and calculating by using the following formula to obtain a discrimination factor F3: f3 ═ 41.7 ═ c (ce)/c (nd)/c (Y') +4.707 (formula 3);
fourthly, calculating a discriminant factor F4
Calculating the formation temperature and oxygen fugacity of the hornblende component obtained in the step (2), carrying out projection by taking T (temperature) as an abscissa and delta FMQ (oxygen fugacity) as an ordinate, obtaining a boundary line between an ore-forming rock mass and a lean rock mass according to the projection range, fitting, and calculating by using the following formula to obtain a discrimination factor F4: f4 ═ 0.0025 × T +3.8 (formula 4).
6. The method according to claim 5, wherein in the step (2), for the single mineral, the data is preprocessed and interpreted, and the preprocessing and the interpretation are divided into the following three steps: firstly, importing the obtained data in the original csv format into ICPMSDataCal software, and eliminating and packaging the data according to the abnormity of the element integration curve of each analysis test pointData of casing or mineral being perforated; secondly, regarding the eliminated data, further eliminating mixed dyeing data by using the following standards: la > 1.5ppm is regarded as mixed-dyeing apatite, Fe>5500ppm reacted the mixed dyeing of Fe oxide, Ti > 60ppm regarded as mixed dyeing of Ti oxide, Ba>10ppm was taken as a mixed dyeing of the fluid inclusions; (iii) use 206 Pb/ 238 U isotope age values exclude interference from inherited minerals, i.e. when data points are obtained 206 Pb/ 238 Weighted average of U isotope ages less than or greater than sample 206 Pb/ 238 At the age of the U isotope, the data point is considered to be invalid data.
CN202210664417.2A 2022-06-13 2022-06-13 Mineral geochemistry-based method for evaluating mineral content of porphyry system Active CN115078520B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202210664417.2A CN115078520B (en) 2022-06-13 2022-06-13 Mineral geochemistry-based method for evaluating mineral content of porphyry system
NL2033632A NL2033632A (en) 2022-06-13 2022-11-28 Method for evaluating ore-bearing potential of porphyry system based on mineral geochemistry

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210664417.2A CN115078520B (en) 2022-06-13 2022-06-13 Mineral geochemistry-based method for evaluating mineral content of porphyry system

Publications (2)

Publication Number Publication Date
CN115078520A true CN115078520A (en) 2022-09-20
CN115078520B CN115078520B (en) 2023-07-04

Family

ID=83252027

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210664417.2A Active CN115078520B (en) 2022-06-13 2022-06-13 Mineral geochemistry-based method for evaluating mineral content of porphyry system

Country Status (2)

Country Link
CN (1) CN115078520B (en)
NL (1) NL2033632A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116297798A (en) * 2023-02-13 2023-06-23 西藏巨龙铜业有限公司 Method for efficiently judging ore-bearing porphyry ore prospect by using paramineral tourmaline
CN116773774A (en) * 2023-06-20 2023-09-19 西藏巨龙铜业有限公司 Method and system for rapidly distinguishing ore forming background of porphyry ore deposit based on tourmaline component

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1497602A1 (en) * 1987-06-26 1989-07-30 Институт Геологии Рудных Месторождений, Петрографии, Минералогии И Геохимии Ан Ссср Method of geochemical prospecting for deposits of rare earths of exogenous ore-forming systems
WO2018212680A1 (en) * 2017-05-17 2018-11-22 Mineral Exploration Network (Finland) Ltd. Geochemical method for searching mineral resource deposits
RU2691418C1 (en) * 2018-04-09 2019-06-13 Акционерная Компания "АЛРОСА" (публичное акционерное общество) (АК "АЛРОСА" (ПАО)) Method for diamond-containing ores preconcentration
CN111122822A (en) * 2019-12-12 2020-05-08 中色华宇(北京)生态科技有限公司 Comprehensive quantification method for formation condition of rock biotinization
CN114384100A (en) * 2021-12-10 2022-04-22 核工业北京地质研究院 Element geochemistry discrimination method beneficial to uranium mineralization granite

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1497602A1 (en) * 1987-06-26 1989-07-30 Институт Геологии Рудных Месторождений, Петрографии, Минералогии И Геохимии Ан Ссср Method of geochemical prospecting for deposits of rare earths of exogenous ore-forming systems
WO2018212680A1 (en) * 2017-05-17 2018-11-22 Mineral Exploration Network (Finland) Ltd. Geochemical method for searching mineral resource deposits
RU2691418C1 (en) * 2018-04-09 2019-06-13 Акционерная Компания "АЛРОСА" (публичное акционерное общество) (АК "АЛРОСА" (ПАО)) Method for diamond-containing ores preconcentration
CN111122822A (en) * 2019-12-12 2020-05-08 中色华宇(北京)生态科技有限公司 Comprehensive quantification method for formation condition of rock biotinization
CN114384100A (en) * 2021-12-10 2022-04-22 核工业北京地质研究院 Element geochemistry discrimination method beneficial to uranium mineralization granite

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DA WANG ETC.: "Two pulse of mineralization and genesis of the Zhaxikang Sb-Pb-Zn-Ag deposit in southern Tibet: Constraints from Fe-Zn isotopes", ORE GEOLOGY REVIEWS, vol. 84, pages 347 - 363, XP029940313, DOI: 10.1016/j.oregeorev.2016.12.030 *
吴松: "西藏冈底斯朱诺超大型斑岩铜矿床:岩浆与成矿", 中国博士学位论文工程科技Ⅰ辑 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116297798A (en) * 2023-02-13 2023-06-23 西藏巨龙铜业有限公司 Method for efficiently judging ore-bearing porphyry ore prospect by using paramineral tourmaline
CN116297798B (en) * 2023-02-13 2023-10-10 西藏巨龙铜业有限公司 Method for efficiently judging ore-bearing porphyry ore prospect by using paramineral tourmaline
CN116773774A (en) * 2023-06-20 2023-09-19 西藏巨龙铜业有限公司 Method and system for rapidly distinguishing ore forming background of porphyry ore deposit based on tourmaline component
CN116773774B (en) * 2023-06-20 2023-12-29 西藏巨龙铜业有限公司 Method and system for rapidly distinguishing ore forming background of porphyry ore deposit based on tourmaline component

Also Published As

Publication number Publication date
CN115078520B (en) 2023-07-04
NL2033632A (en) 2023-12-20

Similar Documents

Publication Publication Date Title
CN115078520A (en) Mineral geochemistry-based porphyry system mineralization evaluation method
Lapakko Metal mine rock and waste characterization tools: an overview
Whitehouse et al. What can Hadean detrital zircon really tell us? A critical evaluation of their geochronology with implications for the interpretation of oxygen and hafnium isotopes
CN114646682B (en) Mineral prospecting method based on trace elements of green-curtain stone
Parbhakar-Fox et al. Cost-effective means for identifying acid rock drainage risks—integration of the geochemistry-mineralogy-texture approach and geometallurgical techniques
Danišík Integration of fission-track thermochronology with other geochronologic methods on single crystals
Krestianinov et al. U-Pb systematics of uranium-rich apatite from Adirondacks: Inferences about regional geological and geochemical evolution, and evaluation of apatite reference materials for in situ dating
Han et al. Multi-scale exploration of giant Qulong porphyry deposit in a collisional setting
Arehart et al. Thermal and isotopic profiling of the Pipeline hydrothermal system: Application to exploration for Carlin-type gold deposits
Whitbread et al. Two lithogeochemical approaches to the identification of alteration patterns at the Elura Zn-Pb-Ag deposit, Cobar, New South Wales, Australia: Use of Pearce element ratio analysis and isocon analysis
RU2360108C1 (en) Method of evaluating porosity and penetrability of oil and gas deposits
Silva et al. A neural network approach for spatial variation assessment–A nepheline syenite case study
Kopp et al. Synchrotron X-ray fluorescence analysis of zoned carbonate gangue in Mississippi Valley-type deposits (USA)
Parbhakar-Fox et al. Predictive waste classification using field-based and environmental geometallurgy indicators, Mount Lyell, Tasmania
Thacker et al. The magnetic susceptibility of cherts: archaeological and geochemical implications of source variation
CN112525939A (en) Open-air PXRF core testing method capable of keeping data accuracy
Liangquan Geochemical prospecting
Downing et al. Trace Element Geochemistry and Acid Rock Drainage
CN115586155B (en) Rapid delineating vein-like gold ore deposit ore body and method for obtaining indication mark index thereof
CN116593407B (en) Rare earth metal mineral rapid investigation device and method
Doherty et al. Stream sediment geochemistry in mineral exploration: a review of fine-fraction, clay-fraction, bulk leach gold, heavy mineral concentrate and indicator mineral chemistry
Kim et al. Determination of Alteration Indexes using a Portable X-ray Fluorescence Analyzer: A Case Study on Hydrothermal Alteration of Limestones from Gangwon Province, South Korea
Urbano et al. Ore-waste and ore type classification using portable XRF: a case study of an iron mine from the Quadrilátero Ferrífero, Brazil.
Hokka et al. Lithogeochemical pXRF study on the Virtasalmi Cu deposit, eastern Finland
Manakhov et al. Speciation of 226 Ra and 232 Th in Soddy-Podzolic-Gley Concretion Soil

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant