CN111508564B - System and method for quantifying the mineralization process of hydrothermal fluid in quartz H2O-CO2-NaCl system - Google Patents

System and method for quantifying the mineralization process of hydrothermal fluid in quartz H2O-CO2-NaCl system Download PDF

Info

Publication number
CN111508564B
CN111508564B CN202010560360.2A CN202010560360A CN111508564B CN 111508564 B CN111508564 B CN 111508564B CN 202010560360 A CN202010560360 A CN 202010560360A CN 111508564 B CN111508564 B CN 111508564B
Authority
CN
China
Prior art keywords
quartz
solubility
fluid
temperature
nacl
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.)
Active
Application number
CN202010560360.2A
Other languages
Chinese (zh)
Other versions
CN111508564A (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.)
National Institute of Natural Hazards
Original Assignee
National Institute of Natural Hazards
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 National Institute of Natural Hazards filed Critical National Institute of Natural Hazards
Priority to CN202010560360.2A priority Critical patent/CN111508564B/en
Publication of CN111508564A publication Critical patent/CN111508564A/en
Priority to JP2020140828A priority patent/JP2021196336A/en
Application granted granted Critical
Publication of CN111508564B publication Critical patent/CN111508564B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10—Analysis or design of chemical reactions, syntheses or processes
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/70—Machine learning, data mining or chemometrics
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P10/00—Technologies related to metal processing
    • Y02P10/20—Recycling

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Data Mining & Analysis (AREA)
  • Analytical Chemistry (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Geometry (AREA)
  • Databases & Information Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Manufacture And Refinement Of Metals (AREA)
  • Investigating Or Analyzing Non-Biological Materials By The Use Of Chemical Means (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The invention discloses a quantitative quartz H 2 O‑CO 2 A system and a method for the hydrothermal fluid mineralization process of a NaCl system, wherein the system is used for researching phase boundary curves, quartz and other solubility curves and quartz and other solubility curve domain boundary curves of different components of mineralization fluids in a mining area fluid temperature and pressure range, and comprises the following steps: h 2 O‑CO 2 -a NaCl system phase boundary calculation module for calculating a liquid phase and gas-liquid mixture phase boundary curve; the solubility calculation module of quartz and the like is used for calculating a solubility curve of quartz and the like; the quartz-like solubility curve domain boundary calculating module is used for calculating the boundary curves of a quartz-like solubility curve domain temperature sensitive area, a quartz-like solubility curve domain pressure sensitive area and a quartz-like solubility curve domain; the curves obtained by the calculation modules are used for combining and constructing quartz and other solubility phase diagrams of different mineral-forming fluid components, and quartz H can be quantified based on the quartz and other solubility phase diagrams 2 O‑CO 2 The invention can provide a new view for understanding the fluid ore-forming process based on the quartz hydrothermal process.

Description

Quantification of quartz H 2 O-CO 2 System and method for hydrothermal fluid ore-forming process of NaCl system
Technical Field
The invention relates to the technical field of mineral deposit geochemistry and calculation geochemistry intersection, in particular to a quantitative quartz H 2 O-CO 2 -a system and method for hydrothermal fluid mineralization of NaCl system.
Background
The current method patents related to the ore forming field are only limited to the aspects of ore forming target area delineation, ore forming lithofacies area positioning, ore forming age prediction, ore forming depth estimation, ore forming space-time positioning, ore forming unit division, ore forming space modeling, ore forming mode construction, ore forming potential calculation and the like (CN 104865613B, CN102243628A, CN110187387A, CN105785466A, CN110060173A, CN107211585B, CN108573206A, CN108181669A, CN107765323A, CN109540929A, CN109270589B, CN 107782878A). There are only 6 patents relating to ore-forming fluids or ore deposit causes, which disclose the discrimination of the redox properties of ore-forming fluids of mainly porphyry deposits (CN 107655915B), ore-forming mode conditions or causes of uranium deposits, etc. (CN 107576996A, CN109752443A, CN111044599A, CN109752443A, CN 106324700B).
Therefore, the research on the ore-forming fluid in the prior art mainly aims at obtaining the physicochemical properties of the fluid in different ore-forming stages, and a system and a method related to quantitative dynamic evolution of the fluid and even ore-forming are not available, so how to create a new system and a method for quantitative dynamic evolution of the fluid and even ore-forming is one of the important research and development problems at present.
Disclosure of Invention
The invention aims to provide a quantitative quartz H 2 O-CO 2 The system and the method for the hydrothermal fluid mineralization process of the NaCl system can quantitatively evaluate the dynamic evolution of the mineralization fluid and even the mineralization process.
In order to solve the technical problems, the invention firstly provides a quantitative quartz H 2 O-CO 2 -a system for a hydrothermal fluid mineralization process of the NaCl system for studying phase boundary curves, quartz-like solubility curve domain boundary curves of different composition mineralization fluids in a range of mining area fluid temperature pressure, comprising:
H 2 O-CO 2 the NaCl system phase boundary calculation module is used for calculating and obtaining the temperature and pressure coordinates (T, P) of a fixed fluid component and a certain point on a phase boundary at a specified temperature, and repeating the calculation module at a selected temperature interval to obtain a series of temperature and pressure coordinates in the temperature and pressure range of the ore-forming fluid in the mining area, namely a liquid phase and gas-liquid mixed phase boundary curve;
H 2 O-CO 2 -a calculation module of the solubility of quartz and the like of NaCl system, which is used for calculating and obtaining the temperature and pressure coordinates (T, P) of a certain temperature when the solubility of the quartz is appointed and obtaining a series of temperature and pressure coordinates in the temperature and pressure range of the ore-forming fluid in the mining area by repeating the calculation module at a selected temperature interval, namely the solubility curve of quartz and the like;
H 2 O-CO 2 The calculating module of the boundary of the solubility curve domain of quartz and the like of the NaCl system is used for obtaining a series of temperature and pressure coordinates of the boundary of the solubility curve domain of quartz and the like under the condition of fixed fluid components, namely the boundary curves of the temperature sensitive area, the pressure sensitive area and the retrograde area of the solubility curve domain of quartz and the like;
the phase boundary curve, the solubility curve of quartz and the like and the solubility curve domain boundary curve of quartz and the like obtained by the calculation modules are used for combining and constructing the solubility phase diagrams of quartz and the like of different components of ore-forming fluid, and quartz H can be quantified based on the solubility phase diagrams of quartz and the like 2 O-CO 2 -a hydrothermal fluid mineralization process of the NaCl system.
As a further improvement of the invention, the sources of the input data of the system are:
(1) Observing the interpenetration relation between an opponent specimen and a quartz vein under a microscope, analyzing a quartz SEM-CL microstructure, and determining the ore forming stage of a deposit and Dan Yingqi times of each ore forming stage by photographing;
(2) Carrying out laser Raman spectrum research on a quartz fluid inclusion to obtain the chemical composition of the ore-forming fluid;
(3) Carrying out microscopic temperature measurement on fluid inclusion in quartz of different periods to obtain uniform temperature, uniform pressure and CO of ore-forming fluid in quartz of each period 2 And NaCl content;
(4) The temperature and pressure range of the ore deposit ore-forming fluid evolution of the research area obtained according to the reference literature is used as the upper limit range and the lower limit range of the temperature and the pressure of a solubility phase diagram such as quartz and the like, and is used as the input value of the temperature and the pressure range of the ore deposit ore-forming fluid of the system;
(5) Mining area mineralizing fluid CO obtained according to the microscopic temperature measurement of the fluid inclusion in (3) 2 And NaCl content range, and the mineral fluid components are freely combined, and a plurality of different fluid components are used as input values of the mineral fluid components of the system.
Further, the H 2 O-CO 2 The calculation flow of the NaCl system phase boundary calculation module is as follows:
101. input temperature T (DEG C), the mineral-forming fluid chemical composition comprises carbon dioxide CO 2 And sodium chloride NaCl, wherein CO 2 The content is expressed as mole percent (mol%) and the content of NaCl is expressed as mass percent (wt%).
102. Chemical composition of ore-forming fluid CO 2 ，NaCl，H 2 O equivalent is converted into mole fraction form and is xCO 2 xNaCl and xH 2 O represents the component of the fluid, and is converted to CO according to the component of the fluid 2 At H 2 Solubility in O-NaCl System in mCO 2 (mol/kg);
103. for pressure P, given its initial interval [ P ] 1 ,P 2 ]By CO 2 At H 2 O-NaCl system solubility model, calculated at pressure P 1 Time CO 2 Solubility mCO 2 P1 And at a pressure of P 2 Time CO 2 Solubility mCO 2 P2 And the interval boundary satisfies the condition: mCO 2 P1 <mCO 2 <mCO 2 P2 Pressure initiation range [ P 1 ,P 2 ]Should be within the applicable range of the model pressure;
104. let p= (P 1 +P 2 )/2；
105. Using CO as described in step 103 2 At H 2 Obtaining the CO by using an O-NaCl system solubility model 2 Solubility calculationValue mCO 2 cal ；
106. When |mCO 2 cal -mCO 2 |>10 -6 And mCO 2 cal <mCO 2 Let P 1 Let P, conversely, let P 2 =p, and jumps to step 104; when |mCO 2 cal -mCO 2 |≤10 -6 When it is, jump to step 107;
107. coordinates (T, P) of points on the boundary line of the liquid phase and the gas-liquid mixed phase are obtained for specifying the mineral-forming fluid component, and the procedure is ended.
Further, the H 2 O-CO 2 The calculation flow of the solubility calculation module of the NaCl system quartz and the like is as follows:
201. input temperature T (DEG C), the mineral-forming fluid chemical composition comprises carbon dioxide CO 2 And sodium chloride NaCl, wherein CO 2 Content is expressed in mole percent (mol%), content of NaCl is expressed in mass percent (wt%), and SiO 2 At H 2 O-CO 2 Given solubility mSiO of NaCl system 2 Expressed in terms of molar mass concentration (mol/kg);
202. to be mineralized fluid component CO 2 ，NaCl，H 2 O equivalent is converted into mole fraction form and is xCO 2 xNaCl and xH 2 O represents the composition of the fluid;
203. for pressure P, given its initial interval [ P ] 1 ,P 2 ]By SiO 2 At H 2 O-CO 2 -NaCl system solubility model calculated at pressure P 1 Quartz solubility mSiO 2 P1 And at a pressure of P 2 Quartz solubility mSiO 2 P2 And the interval boundary should satisfy the condition: mSiO 2 P1 <mSiO 2 <mSiO 2 P2 Pressure initiation range [ P 1 ,P 2 ]Should be included in the applicable range of model pressure;
204. let p= (P 1 +P 2 )/2；
205. Obtaining a calculated value mSiO of the quartz solubility by using the solubility model in the step 203 2 cal ；
206. When |mSiO 2 cal -mSiO 2 |>10 -9 And mSiO 2 cal <mSiO 2 At the time, let P 1 Let P, conversely, let P 2 =p, and jumps to step 204; when |mSiO 2 cal -mSiO 2 |≤10 -9 When it is, jump to step 207;
207. coordinates (T, P) of points satisfying the specified mineral-forming fluid components and quartz solubility are obtained, and the procedure ends.
Further, the H 2 O-CO 2 The calculation flow of the temperature and pressure sensitive area boundary calculation part in the solubility curve domain boundary calculation module of the NaCl system quartz and the like is as follows:
301. input study area SiO 2 At H 2 O-CO 2 Maximum NaCl system solubility mSiO 2 max The mineralizer chemical composition, expressed in terms of molar mass (mol/kg), comprises carbon dioxide CO 2 And sodium chloride NaCl, wherein CO 2 The content is expressed as mole percent (mol%) and the content of NaCl is expressed as mass percent (wt%).
302. To be mineralized fluid component CO 2 ，NaCl，H 2 O equivalent is converted into mole fraction form and is xCO 2 xNaCl and xH 2 O represents the composition of the fluid;
303. setting an initial value mSiO of quartz solubility according to a solubility diagram of quartz and the like 2 ；
304. For the temperature T, given its initial interval [ T ] 1 ,T 2 ]Calculating the temperature as T by using the solubility calculation module of quartz and the like 1 Pressure P at a constant solubility 1 At a temperature T 2 Pressure P at a constant solubility 2 And the interval boundary should satisfy the condition: p (P) 1 >P>P 2 Initial temperature range [ T ] 1 ,T 2 ]Should be within the applicable range of the temperature of the calculation module;
305. let t= (T 1 +T 2 ) 2, calculating to obtain pressure P by utilizing a solubility calculation module such as quartz;
306. according to SiO 2 At H 2 O-CO 2 The NaCl system solubility model shows that the quartz solubility is a function of temperature and pressure under the specified fluid composition, so that mSiO 2 =m (T, P), the partial derivative of the quartz solubility with respect to temperature and pressure at the point (T, P) is calculated using a numerical differentiation method under the current quartz solubility conditionsAnd->
307. When (when)And->At the time, let T 1 Let T, vice versa 2 =t, and jumps to step 305; when->When it is, jump to step 308;
308. obtaining coordinates (T, P) of points on the boundary satisfying the condition;
309. let mSiO 2 ＝mSiO 2 + calculating a step value;
310. when mSiO 2 <mSiO 2 max Jump to step 304; otherwise, go to step 311;
311. And obtaining the temperature and pressure sensitive area boundary of the equal solubility curve domain under the condition of the specified mineral fluid component, and ending the procedure.
Further, the H 2 O-CO 2 The calculation flow of the degenerative region boundary calculation part in the solubility curve region boundary calculation module of the NaCl system quartz and the like is as follows:
401. input study area SiO 2 At H 2 O-CO 2 Maximum NaCl system solubility mSiO 2 max And expressed in terms of molar mass concentration (mol/kg), the mineral-forming fluid component comprises carbon dioxide CO 2 And sodium chloride NaCl, wherein CO 2 The content is expressed as mole percent (mol%) and the content of NaCl is expressed as mass percent (wt%).
402. To be mineralized fluid component CO 2 ，NaCl，H 2 O equivalent is converted into mole fraction form and is xCO 2 xNaCl and xH 2 O represents the composition of the fluid;
403. let mSiO 2 ＝mSiO 2 max ；
404. Setting a temperature T calculation interval [ T ] according to a solubility diagram of quartz and the like 1 ,T 2 ]；
405. Let t= (T 1 +T 2 ) 2, calculating to obtain pressure P by utilizing a solubility calculation module such as quartz;
406. according to SiO 2 At H 2 O-CO 2 -NaCl system solubility model, quartz solubility as a function of temperature and pressure at a given fluid composition, let mSiO 2 =m (T, P), the partial derivative of the quartz solubility with respect to temperature at point (T, P) under the current quartz solubility condition was calculated using a numerical differentiation method
407. If |T 2 -T 1 |<10 -5 ，mSiO 2 ＝mSiO 2 -calculating a step value, otherwise jumping to step 409;
408. when mSiO 2 >When 0, jumping to step 404, otherwise, under the current mineral formation fluid composition condition, no equal solubility curve domain retrogression area exists, and ending the procedure;
409. when (when)And->At the time, let T 1 Let T, vice versa 2 =t, and jumps to step 405; when->At this time, the peak of the degenerative region boundary is obtained, at which the solubility mSiO 2 summit ＝mSiO 2 Vertex coordinates are (T, P);
410. initializing mSiO according to solubility diagram of quartz and the like 2 ；
411. According to the solubility diagram of quartz and the like, pre-calculating the left end point of a retrograde region under the current solubility condition, and setting a temperature T calculation interval [ T ] 1 ,T 2 ]The temperature T is required to contain the temperature corresponding to the left end point;
412. let t= (T 1 +T 2 ) 2, calculating to obtain pressure P by utilizing a solubility calculation module such as quartz;
413. calculating the partial derivative of the solubility of quartz with respect to temperature at point (T, P) under the current solubility of quartz using the method of step 406
414. When (when)And->At the time T 1 =t, vice versa 2 =t, and jumps to step 412; when->When the current solubility condition is obtained, the left end point of the degenerative region is obtained;
415. according to the solubility diagram of quartz and the like, pre-calculating the right end point of the degenerative region under the current solubility condition, and setting a temperature T calculation interval [ T ] 1 ,T 2 ]The temperature T is required to contain the temperature corresponding to the right endpoint;
416. let t= (T 1 +T 2 ) 2, calculating to obtain pressure P by utilizing a solubility calculation module such as quartz;
417、calculating the partial derivative of the solubility of quartz with respect to temperature at point (T, P) under the current solubility of quartz using the method of step 406
418. When (when)And->At the time T 2 =t, vice versa 1 =t, and jumps to step 416; when->Obtaining the right end point of the degenerative region under the current solubility condition;
419、mSiO 2 ＝mSiO 2 + calculating a step value;
420. if mSiO 2 <mSiO 2 summit When it is, jump to step 411; otherwise jump to step 421;
421. and obtaining the degenerative region boundary of the isosolubility curve domain under the condition of fixing the mineralized fluidized components, and ending the procedure.
Based on quantitative quartz H 2 O-CO 2 The invention also provides a system for quantifying quartz H in the hydrothermal fluid ore-forming process of the NaCl system 2 O-CO 2 -a method of hydrothermal fluid mineralization process of NaCl system, comprising:
experimental stage (a):
step (1): observing the interpenetration relation of quartz pulse systems in ore samples in different ore forming stages in a research area under a hand specimen and an optical microscope, and simultaneously researching a scanning electron microscope-cathode luminescence SEM-CL microstructure of quartz, and dividing corresponding fineness Dan Yingqi times in different ore forming stages;
step (2): carrying out laser Raman spectrum research on fluid inclusion in the secondary quartz in different periods of the research area to obtain the ore-forming fluid with the chemical composition of H 2 O-CO 2 -a hydrothermal fluid of NaCl;
step (3): microscopic temperature measurement is carried out on fluid inclusion in the secondary quartz in different periods of the research area to obtain gas phase component CO in the ore-forming fluid 2 Obtaining the content of NaCl as a salt substance and obtaining the uniform temperature and uniform pressure range of ore-forming fluid in quartz of different periods;
step (4): obtaining a temperature and pressure range of ore deposit ore-forming fluid evolution of a research area according to a reference document as upper and lower limits of a solubility phase diagram such as quartz and the like;
step (5): mining fluid CO from investigation region obtained from fluid inclusion temperature measurement 2 And NaCl content range, and taking various different fluid compositions as input values of the mineral fluid components in the calculation stage;
and (II) a calculation stage:
step (6): by H 2 O-CO 2 -NaCl system phase boundary calculation module for calculating H when mineral fluid component is 1 … n 2 O-CO 2 -a point on the boundary of the liquid phase and the gas-liquid mixed phase of the NaCl fluid system;
step (7): repeating the calculation in the step (6) with any temperature step length in the temperature range of the ore deposit ore forming fluid evolution in the research area to obtain a series of points (T, P), namely H 2 O-CO 2 -a phase boundary of a liquid phase and a gas-liquid mixed phase of the NaCl fluid system;
Step (8): by H 2 O-CO 2 -calculating module of solubility of quartz and the like of NaCl system, and calculating H when mineral fluid component is 1 … n respectively 2 O-CO 2 -points on the solubility curve of the NaCl fluid system quartz etc.;
step (9): repeating the calculation in the step (8) with any temperature step length in the temperature range of the ore deposit ore forming fluid evolution in the research area to obtain a series of points (T, P), namely H 2 O-CO 2 -an isosolubility curve for a NaCl fluid system at a given quartz solubility;
step (10): repeating the calculation of the steps (8) and (9) by using any quartz solubility step length to obtain H 2 O-CO 2 -solubility of quartz etc. under different solubility conditions of NaCl fluid systemA curve;
step (11): calculating the boundary of the temperature and pressure sensitive area of the solubility curve domain boundary calculation module of quartz and the like by using the boundary calculation part of the temperature and pressure sensitive area of the solubility curve domain boundary calculation module of quartz and the like when the mineral fluid component is 1 … n;
step (12): calculating the boundary of the solubility curve retrograde region of quartz and the like when the mineral fluid composition is 1 … n by utilizing a retrograde region boundary calculating part in a solubility curve region boundary calculating module of quartz and the like;
and (III) analysis stage:
step (13): by H 2 O-CO 2 -NaCl system phase boundary calculation module, H 2 O-CO 2 -a solubility calculation module of NaCl system quartz and the like, a phase boundary curve obtained by the solubility curve domain boundary calculation module of quartz and the like, a solubility curve of quartz and the like and a solubility curve domain boundary of quartz and the like, and constructing a solubility phase diagram of quartz and the like under the condition that the mineral deposit ore-forming fluid component is 1 … n;
Step (14): the uniform temperature and uniform pressure ranges of the ore-forming fluid in the quartz of different periods obtained in the experimental stage are put in the solubility phase diagrams of quartz and the like under the components of various fluids;
step (15): according to the solubility phase diagram of quartz and the like, the hydrothermal fluid mineralization process is analyzed.
By adopting the technical scheme, the invention has at least the following advantages:
according to the invention, by constructing a system which accords with a solubility phase diagram of quartz and the like of the physical and chemical properties of the ore-forming fluid, a path which is undergone by fluid temperature, pressure, composition and phase change which cause the solubility change of the quartz is utilized, so that the dynamic evolution of the temperature, pressure, composition and phase of the ore-forming fluid is quantitatively clarified, a hydrothermal fluid ore-forming process method is provided, the description of quantitative change of physical and chemical conditions of the fluid ore-forming process in research of traditional ore-deposit science is made up, and the method can provide a new visual angle for understanding the fluid ore-forming process based on the quartz hydrothermal process.
Drawings
The foregoing is merely an overview of the present invention, and the present invention is further described in detail below with reference to the accompanying drawings and detailed description.
FIG. 1 is a quantification of quartz H 2 O-CO 2 -a process diagram of a hydrothermal fluid mineralization process of the NaCl system;
FIG. 2 is a quartz phase sublevel diagram of the ore formation phase of the Dongfeng gold ore deposit;
FIG. 3 is H 2 O-CO 2 -a calculation flow chart of a NaCl system phase boundary calculation module;
FIG. 4 is H 2 O-CO 2 -a phase boundary diagram of a liquid phase and a gas-liquid mixed phase of the NaCl system; wherein, the solid line is the phase boundary of the liquid phase and the gas-liquid mixed phase when the component 1 is adopted, and the dotted line is the phase boundary of the liquid phase and the gas-liquid mixed phase when the component 2 is adopted;
FIG. 5 is H 2 O-CO 2 -a calculation flow chart of a solubility calculation module of quartz or the like of a NaCl system;
FIG. 6 is H 2 O-CO 2 -a solubility profile of NaCl fluid system quartz etc.; wherein the solid line is the solubility curve of quartz and the like in component 1, and the dotted line is the solubility curve of quartz and the like in component 2;
FIG. 7 is H 2 O-CO 2 -a calculation flow chart of a temperature and pressure sensitive area boundary calculation part of a solubility curve domain boundary calculation module of NaCl system quartz and the like;
FIG. 8 is a boundary diagram of a temperature and pressure sensitive region of a solubility curve such as quartz; wherein, the solid line is the boundary of the temperature and pressure sensitive area when the component 1 is used, and the dotted line is the boundary of the temperature and pressure sensitive area when the component 2 is used;
FIG. 9 is a calculation flow chart of the degenerative change region boundary calculation section of the solubility curve region boundary calculation module of quartz or the like;
FIG. 10 is a boundary diagram of a solubility curve domain retrograde region of quartz, etc.; wherein the solid line is the boundary of the degenerative region at component 1 and the dotted line is the boundary of the degenerative region at component 2;
FIG. 11 is a graph of the solubility phase diagram of Dongfeng gold ore bed quartz and the temperature and pressure ranges of quartz of different periods; wherein the solid line is the isosolubility phase diagram of quartz at component 1 and the solid line is the isosolubility phase diagram of quartz at component 2.
Detailed Description
The invention provides a research method of hydrothermal fluid mineralization process related to gangue quartz based on solubility phase diagrams of quartz and the like. To facilitate an understanding of the invention, the following terms are defined:
a. solubility diagram of quartz: in the representation H 2 O-CO 2 -lines connecting states of equal quartz solubility on a state diagram of the NaCl hydrothermal system temperature-pressure state;
b. solubility curve domain of quartz, etc.: the three areas covered by the solubility diagram of quartz and the like comprise a temperature sensitive area, a pressure sensitive area and a retrograde area;
c. pressure sensitive area: in a solubility diagram of quartz and the like, the change of quartz solubility caused by pressure change is larger than the region of quartz solubility change caused by temperature change;
d. temperature sensitive area: in a solubility diagram of quartz and the like, a change in solubility of quartz due to temperature change is larger than a region in which solubility of quartz is changed due to pressure change;
e. degenerative region: in the solubility diagram of quartz, etc., the solubility change of quartz is inversely related to the temperature change.
f. Quartz, etc. solubility phase diagram: in the representation H 2 O-CO 2 On the state diagram of the temperature-pressure state of the NaCl hydrothermal system, lines connecting the states with the same solubility of quartz are overlapped with the boundary of the liquid phase and the gas-liquid mixed phase, and meanwhile, the temperature sensitive area, the pressure sensitive area and the degenerative area of the equal solubility curve area are distinguished.
The invention relies on hydrothermal fluid deposit Dan Yingmai sample slice lithology observation, fluid inclusion laser Raman spectrum test and microscopic temperature measurement and other experimental parts combined with quartz and other solubility phase diagram calculation parts to establish an ore-forming fluid in different COs 2 And under the composition condition of NaCl components, the quartz of different periods is in the temperature and pressure range of the solubility phase diagram of quartz and the like, and the quartz of each period is in the temperature sensitive area, the pressure sensitive area, the retrograde area, the fluid different phase state area and the different fluid related to the solubility phase diagram of quartz and the likeThe evolution of the temperature and the pressure corresponding to the above region caused by the body composition, and the like, quantitatively discusses the evolution paths of the temperature, the pressure and the components possibly experienced by the ore-forming fluid in different periods of quartz and even different ore-forming stages, thereby achieving quantitative evaluation of quartz H 2 O-CO 2 -an mineralization process of a hydrothermal fluid of the NaCl system.
As shown in fig. 1, the experimental part of the mineral flakes contained: (1) Observing the interpenetration relation between an opponent specimen and a quartz vein under a microscope, analyzing a quartz SEM-CL microstructure, and determining the ore forming stage of a deposit and Dan Yingqi times of each ore forming stage by photographing; (2) Carrying out laser Raman spectrum research on a quartz fluid inclusion to obtain the chemical composition of the ore-forming fluid; (3) Carrying out microscopic temperature measurement on fluid inclusion in quartz of different periods to obtain uniform temperature, uniform pressure and CO of ore-forming fluid in quartz of each period 2 And NaCl content.
The calculation part specifically comprises: (1) By H 2 O-CO 2 The NaCl system phase boundary calculation module obtains the temperature and pressure coordinates (T, P) of a certain point on the phase boundary at a specified temperature, and repeats the calculation module at a selected temperature interval to obtain a series of temperature and pressure coordinates, namely the liquid phase and gas-liquid mixed phase boundary; (2) By H 2 O-CO 2 The solubility calculation module of quartz and the like of NaCl system obtains the temperature and pressure coordinates (T, P) of a certain temperature when the solubility of quartz is appointed, and repeats the calculation module at a selected temperature interval to obtain a series of temperature and pressure coordinates, namely the solubility curve of quartz and the like; (3) By H 2 O-CO 2 The calculation module of the solubility curve domain boundary of quartz and the like of the NaCl system obtains the temperature and pressure coordinates (T, P) of a certain point on the solubility curve domain boundary of quartz and the like at a specified temperature, and repeats the calculation module at a selected temperature interval to obtain a series of temperature and pressure coordinates, namely the boundaries of the temperature sensitive area, the pressure sensitive area and the degenerative area of the solubility curve domain of quartz and the like. Here, the calculating part should use quartz laser Raman spectrum test and microscopic temperature measurement to obtain the chemical composition (CO 2 NaCl and H 2 The content of O is expressed as mole fraction x i Expressed) and the chemical composition of the ore-forming fluid can be ground according to the experimental resultsThe method is characterized in that the method is arbitrarily set, and the fluid composition and the temperature and pressure range related to the calculation part of the method are set according to the temperature and the pressure upper limit of the ore formation evolution of the hydrothermal fluid deposit.
The ore forming process analysis part specifically comprises: (1) Obtaining a plurality of groups of solubility phase diagrams of quartz and the like composed of different fluids according to the calculation part; (2) And quantitatively analyzing the effect of the temperature, pressure and component evolution of the quartz hydrothermal fluid on mineral precipitation. It should be noted here that fluid mineralisation process analysis requires the experimental section to provide for each installment of Dan Yingcheng mine fluid chemistry and temperature pressure ranges.
The technical scheme of the invention is described in the following by a specific application example.
Example 1
The Jiaodong ore collection area is the most important primary gold production area in China, and the gold ore production of the Jiaodong ore collection area accounts for about 1/4 of the gold production of China. Many gold deposits in the zone are typically hydrothermal gold ores, quartz being the most important gangue mineral of these hydrothermal deposits. The research on the inclusion of the quartz fluid of the hydrothermal deposit and the quartz precipitation-dissolution behavior can provide a new view and knowledge for the hydrothermal fluid ore-forming process and even the hydrothermal fluid ore-forming mechanism of the hydrothermal deposit. In the embodiment, the hydrothermal solution type gold mine Dongfeng ore deposit in the Jiaodong gold mine collection area is taken as a research object, and the microscopic temperature measurement and the isosolubility phase diagram experimental simulation research of quartz inclusion with ore deposits with different fluid compositions are carried out to quantitatively evaluate quartz H 2 O-CO 2 -an mineralization process of a hydrothermal fluid of the NaCl system.
The embodiment provides quartz H with different mineral forming fluid compositions 2 O-CO 2 -a method of hydrothermal fluid mineralization of NaCl system, the method comprising:
step (1): observing the penetration relation of quartz vein systems in ore samples of different ore forming stages of the Dongfeng gold ore under a hand specimen and an optical microscope, and simultaneously researching a scanning electron microscope-cathode luminescence (SEM-CL) microstructure of quartz, and dividing the corresponding fine Dan Yingqi times of different ore forming stages;
in this embodiment, the sample of the Dongfeng gold ore is divided into three stages of V1, V2 and V3 before, during and after the ore formation, dan Yingqi times of refinement are Qz1, qz2 and Qz3, and the details are shown in fig. 2.
In the ore sample of this example, quartz Qz1 has a strong CL luminescence, qz2 has a poor CL luminescence, and Qz3 has a broken "spider-web" structure under CL, and Qz1 is interspersed with Qz2 and Qz 3.
Step (2): carrying out laser Raman spectrum research on fluid inclusion in different phases of secondary quartz in Dongfeng gold mine to obtain the ore-forming fluid with the chemical composition of H 2 O-CO 2 -a hydrothermal fluid of NaCl;
step (3): microscopic temperature measurement is carried out on fluid inclusion in the secondary quartz of different periods of the Dongfengjinning to obtain gas phase component CO in the ore-forming fluid 2 The content (expressed in terms of mole fraction) of salt substance NaC is obtained, and the uniform temperature and uniform pressure range of the ore-forming fluid in the quartz of different periods are obtained;
when microscopic temperature measurement is carried out on the fluid inclusion, the inclusion with the largest filling degree and the smallest filling degree of the gas-liquid two-phase inclusion captured by the heterogeneous system is selected as much as possible, and the obtained uniform temperature is the formation temperature of minerals, so that pressure correction is not needed; pressure correction is required for the fluid inclusion temperature measurement data of the intermediate filling degree.
Ore forming stage covered by quartz of different stages of Dongfeng gold ore and CO in fluid 2 And NaCl contents are shown in Table 1.
TABLE 1 sub-division of quartz phase of Dongfeng gold deposit in Jiaodong and corresponding fluid composition
Step (4): according to the reference document, the approximate temperature and pressure range of the evolution of the mining fluid of the Jiaodong hydrothermal type gold mine is obtained as the upper limit and the lower limit of the solubility phase diagram of quartz and the like drawn by the method, wherein the temperature range is 100-600 ℃ and the pressure range is 0-4000bar;
step (5): ore-forming fluid CO of Dongfengjingore obtained according to microscopic temperature measurement of fluid inclusion 2 And NaCl content range (10-15 mol% CO) 2 0.4-11.1% wt NaCl) to freely combine into the components of the ore stream. Here we choose CO 2 The mol percent is 10, the NaCl mass percent is 5 (component 1) and CO 2 The two fluid compositions with the mole percentage of 15 and the NaCl mass percentage of 5 (component 2) are used as input values for calculating partial mineral fluid components;
step (6): by H in combination with the illustration of FIG. 3 2 O-CO 2 -NaCl system phase boundary calculation module, calculating H at component 1 2 O-CO 2 -points on the boundary of the NaCl fluid system liquid phase and gas-liquid mixed phase:
(1) input temperature 100 ℃, CO 2 10 mole percent and 5 mass percent of NaCl;
(2) converting the chemical composition of the ore-forming fluid into a representation in mole fraction: xCO 2 ＝0.1,xNaCl＝0.0183,xH 2 O= 0.8817, converted to CO 2 Solubility mCO of (C) 2 ＝6.30mol/kg；
(3) According to Mao et al (Mao s.d., zhang d.h., li y.q.eta., animproved model for calculating CO) 2 solubility in aqueous NaCl solutions and the application to CO 2 –H 2 O-NaCl fluid inclusions, chemical biology, 2013) said CO 2 At H 2 O-NaCl system solubility model, pressure range setting initial interval end point P 1 ＝0，P 2 =4000, the condition can be satisfied: mCO 2 P1 <mCO 2 <mCO 2 P2 ；
(4) Let p= (P 1 +P 2 )/2；
(5) By CO 2 At H 2 O-NaCl system solubility model, and CO when calculating the current pressure P 2 Solubility mCO 2 cal ；
(6) When |mCO 2 cal -mCO 2 |>10 -6 And mCO 2 cal <mCO 2 At the time, let P 1 Let P, conversely, let P 2 =p, and jump to step (4); when the value |mCO 2 cal -mCO 2 |≤10 -6 When P= 2812.07bar is the requirementPressure.
Calculation of component 2H 2 O-CO 2 -points on the boundary of the NaCl fluid system liquid phase and gas-liquid mixed phase:
(1) input temperature 100 ℃, CO 2 15 mol percent and 5 mass percent of NaCl;
(2) converting the chemical composition of the ore-forming fluid into a representation in mole fraction: xCO 2 ＝0.15,xNaCl＝0.0194,xH 2 O= 0.8306, converted to CO 2 Solubility mCO of (C) 2 ＝10.03mol/kg；
(3) According to Mao et al (Mao s.d., zhang d.h., li y.q.et a., an improved model for calculating CO) 2 solubility in aqueous NaCl solutions and the application to CO 2 –H 2 O-NaCl fluid inclusions, chemical biology, 2013) said CO 2 At H 2 O-NaCl system solubility model, pressure range setting initial interval end point P 1 ＝0，P 2 =4000, the condition can be satisfied: mCO 2 P1 <mCO 2 <mCO 2 P2 ；
(4) Let p= (P 1 +P 2 )/2；
(5) By CO 2 At H 2 O-NaCl system solubility model formula, when calculating current pressure P, CO 2 Solubility mCO 2 cal ；
(6) When |mCO 2 cal -mCO 2 |>10 -6 And mCO 2 cal <mCO 2 At the time, let P 1 Let P, conversely, let P 2 =p, and jump to step (4); when the value |mCO 2 cal -mCO 2 |≤10 -6 P= 3040.01bar is the pressure required.
Step (7): repeating the calculation of step (6) at any temperature step (such as 1 ℃ for interval) within the temperature range of 100-600 ℃ to obtain a series of points (T, P), namely H 2 O-CO 2 -phase boundaries of the NaCl fluid system liquid phase and gas-liquid mixed phase (fig. 4);
Step (8): by H in combination with the illustration of FIG. 5 2 O-CO 2 -a module for calculating the solubility of quartz and the like of NaCl system, and calculating H in component 1 2 O-CO 2 Points on the solubility curve of quartz etc. of NaCl fluid system:
(1) input temperature 150 ℃, CO 2 10 mol percent, 5 mass percent of NaCl and SiO 2 At H 2 O-CO 2 Solubility of NaCl System mSiO 2 ＝0.002mol/kg；
(2) Converting the chemical composition of the ore-forming fluid into a representation in mole fraction: xCO 2 ＝0.1,xNaCl＝0.0183,xH 2 O＝0.8817；
(3) According to the use of Wei Qing et al (Wei Qing, duan Zhenhao, mao Shide. H 2 O-CO 2 -quartz solubility model of NaCl system (suitable for high temperature and high pressure environments up to 1000 ℃, 1.5 GPa), rock journal, 2012) of said SiO in 2 At H 2 O-CO 2 -NaCl system solubility (quartz solubility for short) model, setting initial interval endpoint P 1 ＝0，P 2 =4000, the condition can be satisfied: mSiO 2 P1 <mCO 2 <mSiO 2 P2 ；
(4) Let p= (P 1 +P 2 )/2；
(5) Calculating the current pressure P by using a quartz solubility model, and calculating SiO 2 Solubility mSiO 2 cal ；
(6) When |mSiO 2 cal -mSiO 2 |＞10 -9 And mSiO 2 cal <mSiO 2 At the time, let P 1 Let P, conversely, let P 2 =p, and jump to step (4); when |mSiO 2 cal -mSiO 2 |≤10 -9 P= 3408.33bar is the pressure required.
Calculation of component 2H 2 O-CO 2 Points on the solubility curve of quartz etc. of NaCl fluid system:
(1) input temperature 160 ℃, CO 2 15 mol percent, 5 mass percent of NaCl and SiO 2 At H 2 O-CO 2 Solubility of NaCl System mSiO 2 ＝0.002mol/kg；
(2) Converting the chemical composition of the ore-forming fluid into a representation in mole fraction: xCO 2 ＝0.15,xNaCl＝0.0194,xH 2 O＝0.8306；
(3) According to the use of Wei Qing et al (Wei Qing, duan Zhenhao, mao Shide. H 2 O-CO 2 -quartz solubility model of NaCl system (suitable for high temperature and high pressure environments up to 1000 ℃, 1.5 GPa), rock journal, 2012) of said SiO in 2 At H 2 O-CO 2 -NaCl system solubility (quartz solubility for short) model, setting initial interval endpoint P 1 ＝0，P 2 =4000, the condition can be satisfied: mSiO 2 P1 <mCO 2 <mSiO 2 P2 ；
(4) Let p= (P 1 +P 2 )/2；
(5) Calculating the current pressure P by using a quartz solubility model, and calculating SiO 2 Solubility mSiO 2 cal ；
(6) When |mSiO 2 cal -mSiO 2 |＞10 -9 And mSiO 2 cal <mSiO 2 At the time, let P 1 Let P, conversely, let P 2 =p, and jump to step (4); when |mSiO 2 cal -mSiO 2 |≤10 -9 P= 3831.83bar is the pressure required.
Step (9): repeating the calculation of step (8) at any temperature step (such as 1 ℃ for interval) within the temperature range of 100-600 ℃ to obtain a series of points (T, P), namely H 2 O-CO 2 -an isosolubility curve for a NaCl fluid system with a quartz solubility of 0.002 mol/kg;
step (10): repeating the calculation of (8) and (9) with arbitrary quartz solubility step length (such as solubility of 0.001mol/kg as interval) to obtain H 2 O-CO 2 -solubility curves of quartz etc. under different solubility conditions of NaCl fluid system (fig. 6);
Step (11): with reference to fig. 7, the boundary between the temperature and pressure sensitive area and the solubility curve temperature of quartz and the pressure sensitive area at the time of component 1 is calculated by using a solubility curve domain boundary calculating module of quartz and the like, namely a temperature and pressure sensitive area boundary calculating part:
(1) inputting SiO in the temperature and pressure range of the isosolubility phase diagram 2 At H 2 O-CO 2 Maximum solubility (quartz solubility) mSiO of NaCl system 2 max ＝0.087mol/kg，CO 2 10 mole percent and 5 mass percent of NaCl;
(2) converting the chemical composition of the ore-forming fluid into a representation in mole fraction: xCO 2 ＝0.1,xNaCl＝0.0183,xH 2 O＝0.8817；
(3) Let quartz solubility mSiO 2 ＝0.002molg/kg；
(4) Setting an initial interval endpoint T according to a solubility diagram of quartz and the like 1 ＝175，T 2 ＝600；
(5) Let t= (T 1 +T 2 ) 2, calculating to obtain pressure P by utilizing a solubility calculation module such as quartz;
(6) according to Wei Qing et al (Wei Qing, duan Zhenhao, mao Shide. H 2 O-CO 2 -quartz solubility model of NaCl system (suitable for high temperature and high pressure environments up to 1000 ℃, 1.5 GPa), rock journal, 2012) of said SiO in 2 At H 2 O-CO 2 The NaCl system solubility model shows that the quartz solubility is a function of temperature and pressure under the specified fluid composition, so that mSiO 2 =m (T, P). Calculating partial derivative of quartz solubility with respect to temperature and pressure at point (T, P) under current quartz solubility condition by using numerical differentiation method And->
(7) When (when)And->At the time, let T 1 Let T, vice versa 2 =t, and jump to step (5)The method comprises the steps of carrying out a first treatment on the surface of the When->When t= 183.00 ℃, p=14.12 bar is the coordinates (T, P) of the point on the boundary curve that satisfies the condition;
(8) let mSiO 2 ＝mSiO 2 +10 -4 ；
(9) If mSiO 2 <mSiO 2 max And (4) jumping to the step, otherwise ending the program, wherein the obtained series of points are boundaries of the temperature and pressure sensitive areas of the solubility curve of quartz and the like.
Calculating the boundary between the temperature of the solubility curve of quartz and the like and the pressure sensitive area in the component 2:
(1) inputting SiO in the temperature and pressure range of the isosolubility phase diagram 2 At H 2 O-CO 2 Maximum solubility (quartz solubility) mSiO of NaCl system 2 max ＝0.068mol/kg，CO 2 15 mol percent and 5 mass percent of NaCl;
(2) converting the chemical composition of the ore-forming fluid into a representation in mole fraction: xCO 2 ＝0.15,xNaCl＝0.0194,xH 2 O＝0.8306；
(3) Let quartz solubility mSiO 2 ＝0.002molg/kg；
(4) Setting an initial interval endpoint T according to a solubility diagram of quartz and the like 1 ＝190，T 2 ＝600；
(5) Let t= (T 1 +T 2 ) 2, calculating to obtain pressure P by utilizing a solubility calculation module such as quartz;
(6) according to Wei Qing et al (Wei Qing, duan Zhenhao, mao Shide. H 2 O-CO 2 -quartz solubility model of NaCl system (suitable for high temperature and high pressure environments up to 1000 ℃, 1.5 GPa), rock journal, 2012) of said SiO in 2 At H 2 O-CO 2 The NaCl system solubility model shows that the quartz solubility is a function of temperature and pressure under the specified fluid composition, so that mSiO 2 =m (T, P). Calculating by numerical differentiation to obtain the temperature and pressure of the quartz solubility at the point (T, P) under the current quartz solubility conditionPartial derivativeAnd->
(7) When (when)And->At the time, let T 1 Let T, vice versa 2 =t, and jump to step (5); when->When t= 222.50 ℃, p=26.76 bar is the coordinates (T, P) of the point on the boundary curve that satisfies the condition;
(8) let mSiO 2 ＝mSiO 2 +10 -4 ；
(9) If mSiO 2 <mSiO 2 max And (4) jumping to the step, otherwise ending the procedure, wherein the obtained series of points are the boundaries of the temperature and pressure sensitive areas of the solubility curve fields such as quartz and the like (fig. 8).
Step (12): with reference to fig. 9, the boundary of the solubility curve region of quartz and the like at the time of component 1 is calculated by using the solubility curve region boundary calculation module of quartz and the like—the region boundary calculation section of the regression:
1. inputting SiO in the temperature and pressure range of the isosolubility phase diagram 2 At H 2 O-CO 2 Maximum solubility (quartz solubility) mSiO of NaCl system 2 max ＝0.02mol/kg，CO 2 10 mole percent and 5 mass percent of NaCl;
2. converting the chemical composition of the ore-forming fluid into a representation in mole fraction: xCO 2 ＝0.1,xNaCl＝0.0183,xH 2 O＝0.8817；
3. Let quartz solubility mSiO 2 ＝mSiO 2 max
4. Setting an initial interval endpoint T according to a solubility diagram of quartz and the like 1 ＝450，T 2 ＝600；
5. Let t= (T 1 +T 2 ) 2, calculating to obtain pressure P by utilizing a solubility calculation module such as quartz;
6. According to Wei Qing et al (Wei Qing, duan Zhenhao, mao Shide. H 2 O-CO 2 -quartz solubility model of NaCl system (suitable for high temperature and high pressure environments up to 1000 ℃, 1.5 GPa), rock journal, 2012) of said SiO in 2 At H 2 O-CO 2 The NaCl system solubility model shows that the quartz solubility is a function of temperature and pressure under the specified fluid composition, so that mSiO 2 =m (T, P). Calculating partial derivative of quartz solubility with respect to temperature and pressure at point (T, P) under current quartz solubility condition by using numerical differentiation method
7. If |T 2 -T 1 |<10 -5 ，mSiO 2 ＝mSiO 2 -10 -4 (here 10 -4 In order to calculate the step length value, the step length value can be adjusted according to the actual calculation condition, otherwise, the step length value is jumped to 9;
8. when mSiO 2 >When 0, jumping to 4, otherwise, under the current mineral fluid composition condition, no equal solubility curve retrogression area exists, and ending the procedure;
9. when (when)And->At the time, let T 1 Let T, vice versa 2 =t, and jump to step (5); when->When mSiO 2 summit =0.0121 mol/kg, and t= 568.75 ℃, p= 649.79bar is the boundary satisfying the conditionCoordinates of vertices on the curve;
10. let mSiO 2 ＝0.002；
11. Calculating the left end point of the boundary of the degenerative region under the current solubility condition, and setting an initial region end point T according to a solubility diagram of quartz and the like 1 ＝180，T 2 ＝200；
12. Let t= (T 1 +T 2 ) 2, calculating to obtain pressure P by utilizing a solubility calculation module such as quartz;
13. obtaining the partial derivative of the quartz solubility with respect to the temperature and the pressure at the point (T, P) under the current quartz solubility condition by using the method of (6)
14. When (when)And->At the time T 1 =t, vice versa 2 =t, jump to 12; when->When t= 184.00 ℃, p=10.98 bar is the coordinate of the left end point of the boundary of the downward-retrogressive zone at the current mineral fluid composition.
15. Calculating the right end point of the boundary of the degenerative region under the current solubility condition, and setting an initial region end point T according to a solubility diagram of quartz and the like 1 ＝450，T 2 ＝600；
16. Let t= (T 1 +T 2 ) 2, calculating to obtain pressure P by utilizing a solubility calculation module such as quartz;
17. calculating the partial derivative of the quartz solubility with respect to the temperature and the pressure at the point (T, P) under the current quartz solubility condition by using the method of (6)
18. When (when)And->At the time T 2 =t, vice versa 1 =t, jump to 16; when (when)At t= 523.05 ℃, p= 350.65bar is the coordinates of the right end point of the zone boundary of the zone of retrogradation at a quartz solubility of 0.002mol/kg at the current mineralisation fluid composition;
19、mSiO 2 ＝mSiO 2 +10 -4 ；
20. if mSiO 2 <mSiO 2 max Jump to 11 if mSiO 2 ≥mSiO 2 max When this procedure is completed, the obtained series of points are boundaries of the solubility curve retrograde region of quartz or the like (fig. 10).
Calculating the boundary of the quartz-like solubility curve retrograde region at component 2:
1. inputting SiO in the temperature and pressure range of the isosolubility phase diagram 2 At H 2 O-CO 2 Maximum solubility (quartz solubility) mSiO of NaCl system 2 max ＝0.01mol/kg，CO 2 15 mol percent and 5 mass percent of NaCl;
2. converting the chemical composition of the ore-forming fluid into a representation in mole fraction: xCO 2 ＝0.15,xNaCl＝0.0194,xH 2 O＝0.8306；
3. Let quartz solubility mSiO 2 ＝mSiO 2 max ；
4. Setting an initial interval endpoint T according to a solubility diagram of quartz and the like 1 ＝450，T 2 ＝600；
5. Let t= (T 1 +T 2 ) 2, calculating to obtain pressure P by utilizing a solubility calculation module such as quartz;
6. according to Wei Qing et al (Wei Qing, duan Zhenhao, mao Shide. H 2 O-CO 2 -quartz solubility model of NaCl system (suitable for high temperature and high pressure environments up to 1000 ℃, 1.5 GPa), rock journal, 2012) of said SiO in 2 At H 2 O-CO 2 The NaCl system solubility model shows that the quartz solubility is a function of temperature and pressure under the specified fluid composition, so that mSiO 2 =m (T, P). Calculating partial derivative of quartz solubility with respect to temperature and pressure at point (T, P) under current quartz solubility condition by using numerical differentiation method
7. If |T 2 -T 1 |<10 -5 ，mSiO 2 ＝mSiO 2 -10 -4 (here 10 -4 In order to calculate the step length value, the step length value can be adjusted according to the actual calculation condition, otherwise, the step length value is jumped to 9;
8. When mSiO 2 >When 0, jumping to 4, otherwise, under the current mineral fluid composition condition, no equal solubility curve retrogression area exists, and ending the procedure;
9. when (when)And->At the time, let T 1 Let T, vice versa 2 =t, and jump to 5; when->When mSiO 2 =0.0057 mol/kg, and t=525.00 ℃, p= 548.22bar is the coordinates of the vertices on the boundary curve that satisfy the condition;
10. let mSiO 2 ＝0.002；
11. Calculating the left end point of the boundary of the degenerative region under the current solubility condition, and setting an initial region end point T according to a solubility diagram of quartz and the like 1 ＝180，T 2 ＝200；
12. Let t= (T 1 +T 2 ) 2, utilizing solubility of quartz or the likeThe calculation module is used for calculating to obtain pressure P;
13. obtaining the partial derivative of the quartz solubility with respect to the temperature and the pressure at the point (T, P) under the current quartz solubility condition by using the method of (6)
14. When (when)And->At the time T 1 =t, vice versa 2 =t, jump to 12; when (when)When t= 223.05 ℃, p=24.58 bar is the coordinate of the left end point of the boundary of the downward-retrogressive zone at the current mineral fluid composition.
15. Calculating the right end point of the boundary of the degenerative region under the current solubility condition, and setting an initial region end point T according to a solubility diagram of quartz and the like 1 ＝450，T 2 ＝550；
16. Let t= (T 1 +T 2 ) 2, calculating to obtain pressure P by utilizing a solubility calculation module such as quartz;
17. Calculating the partial derivative of the quartz solubility with respect to the temperature and the pressure at the point (T, P) under the current quartz solubility condition by using the method of (6)
18. When (when)And->At the time T 2 =t, vice versa 1 =t, jump to 16; when (when)When t= 528.83 ℃, p= 385.25bar is the right end point of the zone boundary where the current ore-forming fluid composition is down-graded;
19、mSiO 2 ＝mSiO 2 +10 -4 ；
20. if mSiO 2 <mSiO 2 max Jump to 11 if mSiO 2 ≥mSiO 2 max When this procedure is completed, the obtained series of points are boundaries of the solubility curve domain degenerated region of quartz or the like (fig. 10).
Step (13): by H 2 O-CO 2 -a NaCl system phase boundary calculation module, a quartz and other solubility curve domain boundary calculation module, a quartz and other solubility curve and quartz and other solubility curve domain boundary, and constructing a quartz and other solubility phase diagram under the condition that the fluid composition of the Dongfeng gold deposit is respectively component 1 and component 2;
step (14): the uniform temperature and uniform pressure range of the ore-forming fluid in quartz of different periods of the Dongfeng gold ore obtained in the experimental part are put in a solubility phase diagram of quartz and the like under the composition of two fluids, as shown in fig. 11;
step (15): the hydrothermal fluid mineralization process is analyzed according to the solubility phase diagram (figure 11) of the eastern galium ore bed quartz and the like.
In the embodiment, when the components of the ore-forming fluid of the downcast gold ore are two components, namely a component 1 and a component 2, qz1 falls in a temperature sensitive area of a liquid phase area of a solubility phase diagram of quartz and the like; when the component of the ore forming fluid is the component 1, qz2 falls in a temperature sensitive area of a liquid phase region of a solubility phase diagram of quartz and the like, but when the component of the ore forming fluid is the component 2, qz2 spans the temperature sensitive area of the liquid phase region of the solubility phase diagram of quartz and the like and a gas-liquid mixed phase; when the component of the ore-forming fluid is the component 1, qz3 spans three areas of a temperature sensitive area, a pressure sensitive area and a retrograde area of a gas-liquid two-phase area of a solubility phase diagram of quartz and the like, but when the component of the ore-forming fluid is the component 2, qz3 only spans the pressure sensitive area and the retrograde area of a gas-liquid mixed phase of the solubility phase diagram of quartz and the like.
For the embodiment of east windGold ore fluid evolves, and for a certain hypothesized temperature and pressure range (550-600 ℃,3200-3600 bar) corresponding to the initial fluid, qz1 is precipitated from the initial fluid and can be formed by cooling the ore-forming fluid from 600 ℃ to 320 ℃ in a liquid phase region; the quartz is changed from Qz1 to Qz2, and the temperature of the ore-forming fluid in the liquid phase region is reduced from 370 ℃ to 300 ℃ or the CO in the fluid 2 The mole fraction is formed by gas-liquid phase separation caused by the increase of ten percent to fifteen percent or the superposition of the gas and the liquid phase separation; quartz evolves from Qz2 to Qz3, and the pressure of the ore-forming fluid, which may be caused by gas-liquid phase separation of the ore-forming fluid, is reduced from 1600bar to 200bar, or the degenerative dissolution behavior of quartz, or a superposition of both.
Combining the formation processes of Qz1, qz2 and Qz3 according to Dongfengjingkou quartz H 2 O-CO 2 The solubility phase diagram of NaCl system, etc. shows that the liquid phase area of the Dongfeng gold ore is first cooled to 320 deg.c and small amount of metal matter is unloaded, and this stage is the pre-ore forming stage of Dongfeng gold ore (stage V1); further, with further decreases in the temperature of the ore-forming fluid (from about 370 ℃ to 300 ℃) and the nonpolar gas CO in the fluid 2 A fluid phase separation process caused by the increase of the content (the mole percentage of the system is increased from ten to fifteen), a large amount of metal substances are unloaded from the ore-forming fluid, and the stage is the main ore-forming stage (V2 stage) of the eastern anemone gold ore; with further reduction of the temperature and pressure of the ore-forming fluid (temperature reduction from 350 ℃ to 270 ℃ and pressure reduction from 1600bar to about 200 bar), the ore-forming fluid undergoes a significant phase separation stage (V3 stage), and the rapid release of pressure causes precipitation of the metal species, which also precipitates slightly with a small increase in temperature in the narrower temperature interval of the quartz retrograde zone (280 ℃ -330 ℃).
Comprehensive above, it is Dongfengjingore quantitative quartz H 2 O-CO 2 -an mineralization process of a hydrothermal fluid of the NaCl system.
The above description is only of the preferred embodiments of the present invention, and is not intended to limit the invention in any way, and some simple modifications, equivalent variations or modifications can be made by those skilled in the art using the teachings disclosed herein, which fall within the scope of the present invention.

Claims (3)

1. Quantification of quartz H 2 O-CO 2 -a system for a hydrothermal fluid mineralization process of the NaCl system, characterized in that it is used for studying phase boundary curves, solubility curves such as quartz, solubility curve domain boundary curves such as quartz of different components of the mineralized fluid in the temperature and pressure range of the fluid in the mining area, comprising:
H 2 O-CO 2 the NaCl system phase boundary calculation module is used for calculating and obtaining the temperature and pressure coordinates (T, P) of a certain point on a phase boundary at a fixed fluid component and a specified temperature, wherein T is temperature, P is pressure, and the calculation module is repeated at a selected temperature interval to obtain a series of temperature and pressure coordinates in the temperature and pressure range of the ore-forming fluid in a mining area, namely a liquid phase and gas-liquid mixed phase boundary curve;
H 2 O-CO 2 -a calculation module of the solubility of quartz and the like of a NaCl system, which is used for calculating and obtaining a temperature-pressure coordinate (T, P) of a certain temperature when the solubility of the quartz is appointed and a fixed fluid component, wherein T is the temperature, P is the pressure, and the calculation module is repeated at a selected temperature interval to obtain a series of temperature-pressure coordinates under the condition of the solubility of the quartz and the like in the temperature-pressure range of the ore-forming fluid in a mining area, namely a solubility curve of the quartz and the like;
H 2 O-CO 2 The calculating module of the boundary of the solubility curve domain of quartz and the like of the NaCl system is used for obtaining a series of temperature and pressure coordinates of the boundary of the solubility curve domain of quartz and the like under the condition of fixed fluid components, namely the boundary curves of the temperature sensitive area, the pressure sensitive area and the retrograde area of the solubility curve domain of quartz and the like;
the phase boundary curve, the solubility curve of quartz and the like and the solubility curve domain boundary curve of quartz and the like obtained by the calculation modules are used for combining and constructing the solubility phase diagrams of quartz and the like of different components of ore-forming fluid, and quartz H can be quantified based on the solubility phase diagrams of quartz and the like 2 O-CO 2 -a hydrothermal fluid mineralization process of the NaCl system;
the H is 2 O-CO 2 The calculation flow of the NaCl system phase boundary calculation module is as follows:
101. input temperature T (DEG C), the mineral-forming fluid chemical composition comprises carbon dioxide CO 2 And sodium chloride NaCl, wherein CO 2 The content is expressed as mole percent (mol%) and the content of NaCl is expressed as mass percent (wt%).
102. Chemical composition of ore-forming fluid CO 2 ，NaCl，H 2 O equivalent is converted into mole fraction form and is xCO 2 xNaCl and xH 2 O represents the component of the fluid, and is converted to CO according to the component of the fluid 2 At H 2 Solubility in O-NaCl System in mCO 2 (mol/kg); wherein x is a mole fraction and m is a molar mass concentration;
103. For pressure P, given its initial interval [ P ] 1 ,P 2 ]The method comprises the steps of carrying out a first treatment on the surface of the By CO 2 At H 2 O-NaCl system solubility model, calculated at pressure P 1 Time CO 2 Solubility mCO 2 P1 And at a pressure of P 2 Time CO 2 Solubility mCO 2 P2 And the interval boundary satisfies the condition: mCO 2 P1 <mCO 2 <mCO 2 P2 Pressure initiation range [ P 1 ,P 2 ]Should be within the applicable range of the model pressure;
104. let p= (P 1 +P 2 )/2；
105. Using CO as described in step 103 2 At H 2 Obtaining the CO by using an O-NaCl system solubility model 2 Solubility calculation mCO 2 cal ；
106. When |mCO 2 cal -mCO 2 |>10 -6 And mCO 2 cal <mCO 2 Let P 1 Let P, conversely, let P 2 =p, and jumps to step 104; when |mCO 2 cal -mCO 2 |≤10 -6 When it is, jump to step 107;
107. obtaining coordinates (T, P) of points on the boundary line of the liquid phase and the gas-liquid mixed phase, which specify the components of the ore-forming fluid, and ending the process;
the H is 2 O-CO 2 The calculation flow of the solubility calculation module of the NaCl system quartz and the like is as follows:
201. input temperature T (DEG C), the mineral-forming fluid chemical composition comprises carbon dioxide CO 2 And sodium chloride NaCl, wherein CO 2 Content is expressed in mole percent (mol%), content of NaCl is expressed in mass percent (wt%), and SiO 2 At H 2 O-CO 2 Given solubility mSiO of NaCl system 2 Expressed as molar mass (mol/kg), m is the molar mass;
202. to be mineralized fluid component CO 2 ，NaCl，H 2 O equivalent is converted into mole fraction form and is xCO 2 xNaCl and xH 2 O represents the composition of the fluid;
203. for pressure P, given its initial interval [ P ] 1 ,P 2 ]By SiO 2 At H 2 O-CO 2 -NaCl system solubility model calculated at pressure P 1 Quartz solubility mSiO 2 P1 And at a pressure of P 2 Quartz solubility mSiO 2 P2 And the interval boundary should satisfy the condition: mSiO 2 P1 <mSiO 2 <mSiO 2 P2 Pressure initiation range [ P 1 ,P 2 ]Should be included in the applicable range of model pressure;
204. let p= (P 1 +P 2 )/2；
205. Obtaining a calculated value mSiO of the quartz solubility by using the solubility model in the step 203 2 cal ；
206. When |mSiO 2 cal -mSiO 2 |>10 -9 And mSiO 2 cal <mSiO 2 At the time, let P 1 Let P, conversely, let P 2 =p, and jumps to step 204; when |mSiO 2 cal -mSiO 2 |≤10 -9 When it is, jump to step 207;
207. obtaining coordinates (T, P) of points satisfying the specified mineral-forming fluid composition and quartz solubility, and ending the procedure;
the H is 2 O-CO 2 NaCl-stoneThe calculation flow of the temperature and pressure sensitive area boundary calculation part in the British equivalent solubility curve domain boundary calculation module is as follows:
301. input study area SiO 2 At H 2 O-CO 2 Maximum NaCl system solubility mSiO 2 max The mineralizer chemical composition, expressed in terms of molar mass (mol/kg), comprises carbon dioxide CO 2 And sodium chloride NaCl, wherein CO 2 The content is expressed as mole percent (mol%) and the content of NaCl is expressed as mass percent (wt%).
302. To be mineralized fluid component CO 2 ，NaCl，H 2 O equivalent is converted into mole fraction form and is xCO 2 xNaCl and xH 2 O represents the composition of the fluid;
303. setting a quartz solubility initial value mSiO2 according to a quartz solubility diagram;
304. for the temperature T, given its initial interval [ T ] 1 ,T 2 ]The method comprises the steps of carrying out a first treatment on the surface of the Calculating the temperature to be T by using the solubility calculation module of quartz and the like 1 Pressure P at a constant solubility 1 At a temperature T 2 Pressure P at a constant solubility 2 And the interval boundary should satisfy the condition: p (P) 1 >P>P 2 Initial temperature range [ T ] 1 ,T 2 ]Should be within the applicable range of the temperature of the calculation module;
305. let t= (T 1 +T 2 ) 2, calculating to obtain pressure P by utilizing a solubility calculation module such as quartz;
306. according to SiO 2 At H 2 O-CO 2 The NaCl system solubility model shows that the quartz solubility is a function of temperature and pressure under the specified fluid composition, so that mSiO 2 =m (T, P); calculating partial derivative of quartz solubility with respect to temperature and pressure at point (T, P) under current quartz solubility condition by using numerical differentiation methodAnd->Said->The partial derivative of the quartz solubility function with respect to the temperature T under the condition of the fixed pressure P;The partial derivative of the quartz solubility function with respect to the pressure P at a fixed temperature T;
307. When (when)And->At the time, let T 1 Let T, vice versa 2 =t, and jumps to step 305; when->When it is, jump to step 308;
308. obtaining coordinates (T, P) of points on the boundary satisfying the condition;
309. let mSiO 2 ＝mSiO 2 + calculating a step value;
310. when mSiO 2 <mSiO 2 max Jump to step 304; otherwise, go to step 311;
311. obtaining the temperature and pressure sensitive area boundary of the equal solubility curve domain under the condition of appointed ore-forming fluid components, and ending the procedure;
the H is 2 O-CO 2 The calculation flow of the degenerative region boundary calculation part in the solubility curve region boundary calculation module of the NaCl system quartz and the like is as follows:
401. input study area SiO 2 At H 2 O-CO 2 Maximum NaCl system solubility mSiO 2 max And expressed in terms of molar mass concentration (mol/kg), the mineral-forming fluid component comprises carbon dioxide CO 2 And sodium chloride NaCl, whichMedium CO 2 The content is expressed as mole percent (mol%) and the content of NaCl is expressed as mass percent (wt%).
402. To be mineralized fluid component CO 2 ，NaCl，H 2 O equivalent is converted into mole fraction form and is xCO 2 xNaCl and xH 2 O represents the composition of the fluid;
403. let mSiO 2 ＝mSiO 2 max ；
404. Setting a temperature T calculation interval [ T ] according to a solubility diagram of quartz and the like 1 ,T 2 ]；
405. Let t= (T 1 +T 2 ) 2, calculating to obtain pressure P by utilizing a solubility calculation module such as quartz;
406. According to SiO 2 At H 2 O-CO 2 -NaCl system solubility model, quartz solubility as a function of temperature and pressure at a given fluid composition, let mSiO 2 =m (T, P), the partial derivative of the quartz solubility with respect to temperature at point (T, P) under the current quartz solubility condition was calculated using a numerical differentiation method407. If |T 2 -T 1 |<10 -5 ，mSiO 2 ＝mSiO 2 -calculating a step value, otherwise jumping to step 409;
408. when mSiO 2 >When 0, jumping to step 404, otherwise, under the current mineral formation fluid composition condition, no equal solubility curve domain retrogression area exists, and ending the procedure;
409. when (when)And->At the time, let T 1 Let T, vice versa 2 =t, and jumps to step 405; when->At this time, the peak of the degenerative region boundary is obtained, at which the solubility mSiO 2 summit ＝mSiO 2， Vertex coordinates are (T, P);
410. initializing mSiO according to solubility diagram of quartz and the like 2 ；
411. According to the solubility diagram of quartz and the like, pre-calculating the left end point of a retrograde region under the current solubility condition, and setting a temperature T calculation interval [ T ] 1 ,T 2 ]The temperature T is required to contain the temperature corresponding to the left end point;
412. let t= (T 1 +T 2 ) 2, calculating to obtain pressure P by utilizing a solubility calculation module such as quartz;
413. calculating the partial derivative of the solubility of quartz with respect to temperature at point (T, P) under the current solubility of quartz using the method of step 406
414. When (when)And->At the time T 1 =t, vice versa 2 =t, and jumps to step 412, when +.>When the current solubility condition is obtained, the left end point of the degenerative region is obtained;
415. according to the solubility diagram of quartz and the like, pre-calculating the right end point of the degenerative region under the current solubility condition, and setting a temperature T calculation interval [ T ] 1 ,T 2 ]The temperature T is required to contain the temperature corresponding to the right endpoint;
416. let t= (T 1 +T 2 ) 2, calculating to obtain pressure P by utilizing a solubility calculation module such as quartz;
417. calculating the dissolution of quartz at point (T, P) under the current quartz solubility conditions using the method of step 406Partial derivative of degree with respect to temperature
418. When (when)And->At the time T 2 =t, vice versa 1 =t, and jumps to step 416; when->Obtaining the right end point of the degenerative region under the current solubility condition;
419、mSiO 2 ＝mSiO 2 + calculating a step value;
420. if mSiO 2 <mSiO 2 summit When it is, jump to step 411; otherwise jump to step 421;
421. and obtaining the degenerative region boundary of the isosolubility curve domain under the condition of fixing the mineralized fluidized components, and ending the procedure.
2. Quantitative quartz H according to claim 1 2 O-CO 2 -a system of hydrothermal fluid mineralization process of NaCl system, characterized in that the source of input data of the system is:
(1) Observing the interpenetration relation between an opponent specimen and a quartz vein under a microscope, analyzing a quartz SEM-CL microstructure, and determining the ore forming stage of a deposit and Dan Yingqi times of each ore forming stage by photographing;
(2) Carrying out laser Raman spectrum research on a quartz fluid inclusion to obtain the chemical composition of the ore-forming fluid;
(3) Carrying out microscopic temperature measurement on fluid inclusion in quartz of different periods to obtain uniform temperature, uniform pressure and CO of ore-forming fluid in quartz of each period 2 And NaCl content;
(4) The temperature and pressure range of the ore deposit ore-forming fluid evolution of the research area obtained according to the reference literature is used as the upper limit range and the lower limit range of the temperature and the pressure of a solubility phase diagram such as quartz and the like, and is used as the input value of the temperature and the pressure range of the ore deposit ore-forming fluid of the system;
(5) Mining area mineralizing fluid CO obtained according to the microscopic temperature measurement of the fluid inclusion in (3) 2 And NaCl content range, and the mineral fluid components are freely combined, and a plurality of different fluid components are used as input values of the mineral fluid components of the system.
3. Quantitative quartz H 2 O-CO 2 Method for the hydrothermal fluid mineralization of the NaCl system, characterized in that the quantitative quartz H according to any of claims 1-2 is used 2 O-CO 2 -a system of hydrothermal fluid mineralization process of NaCl system, performing a computational analysis, comprising:
experimental stage (a):
step (1): observing the interpenetration relation of quartz pulse systems in ore samples in different ore forming stages in a research area under a hand specimen and an optical microscope, and simultaneously researching a scanning electron microscope-cathode luminescence SEM-CL microstructure of quartz, and dividing corresponding fineness Dan Yingqi times in different ore forming stages;
Step (2): carrying out laser Raman spectrum research on fluid inclusion in the secondary quartz in different periods of the research area to obtain the ore-forming fluid with the chemical composition of H 2 O-CO 2 -a hydrothermal fluid of NaCl;
step (3): microscopic temperature measurement is carried out on fluid inclusion in the secondary quartz in different periods of the research area to obtain gas phase component CO in the ore-forming fluid 2 Obtaining the content of NaCl as a salt substance and obtaining the uniform temperature and uniform pressure range of ore-forming fluid in quartz of different periods;
step (4): obtaining a temperature and pressure range of ore deposit ore-forming fluid evolution of a research area according to a reference document as upper and lower limits of a solubility phase diagram such as quartz and the like;
step (5): mining fluid CO from investigation region obtained from fluid inclusion temperature measurement 2 And NaCl content range, and using various fluid compositions as calculation stage ore-forming fluidInput values of the body composition;
and (II) a calculation stage:
step (6): by H 2 O-CO 2 -NaCl system phase boundary calculation module for calculating H when mineral fluid component is 1 … n 2 O-CO 2 -a point on the boundary of the liquid phase and the gas-liquid mixed phase of the NaCl fluid system;
step (7): repeating the calculation in the step (6) with any temperature step length in the temperature range of the ore deposit ore forming fluid evolution in the research area to obtain a series of points (T, P), namely H 2 O-CO 2 -a phase boundary of a liquid phase and a gas-liquid mixed phase of the NaCl fluid system;
step (8): by H 2 O-CO 2 -calculating module of solubility of quartz and the like of NaCl system, and calculating H when mineral fluid component is 1 … n respectively 2 O-CO 2 -points on the solubility curve of the NaCl fluid system quartz etc.;
step (9): repeating the calculation in the step (8) with any temperature step length in the temperature range of the ore deposit ore forming fluid evolution in the research area to obtain a series of points (T, P), namely H 2 O-CO 2 -an isosolubility curve for a NaCl fluid system at a given quartz solubility;
step (10): repeating the calculation of the steps (8) and (9) by using any quartz solubility step length to obtain H 2 O-CO 2 -solubility curves of quartz etc. under different solubility conditions of NaCl fluid system;
step (11): calculating the boundary of the temperature and pressure sensitive area of the solubility curve domain boundary calculation module of quartz and the like by using the boundary calculation part of the temperature and pressure sensitive area of the solubility curve domain boundary calculation module of quartz and the like when the mineral fluid component is 1 … n;
step (12): calculating the boundary of the solubility curve retrograde region of quartz and the like when the mineral fluid composition is 1 … n by utilizing a retrograde region boundary calculating part in a solubility curve region boundary calculating module of quartz and the like;
and (III) analysis stage:
step (13): by H 2 O-CO 2 -NaCl system phase boundary calculation module, H 2 O-CO 2 -NaCl system quartz etc. solubility calculation mouldA phase boundary curve, a solubility curve of quartz and the like and a solubility curve domain boundary of quartz and the like obtained by the solubility curve domain boundary calculation module of the block, the quartz and the like are used for constructing a solubility phase diagram of quartz and the like under the condition that the mineral fluid component of the ore deposit is 1 … n;
step (14): the uniform temperature and uniform pressure ranges of the ore-forming fluid in the quartz of different periods obtained in the experimental stage are put in the solubility phase diagrams of quartz and the like under the components of various fluids;
step (15): according to the solubility phase diagram of quartz and the like, the hydrothermal fluid mineralization process is analyzed.
CN202010560360.2A 2020-06-18 2020-06-18 System and method for quantifying the mineralization process of hydrothermal fluid in quartz H2O-CO2-NaCl system Active CN111508564B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010560360.2A CN111508564B (en) 2020-06-18 2020-06-18 System and method for quantifying the mineralization process of hydrothermal fluid in quartz H2O-CO2-NaCl system
JP2020140828A JP2021196336A (en) 2020-06-18 2020-08-24 THEORETICAL MODEL AND METHOD FOR QUANTIFYING HYDROTHERMAL FLUID MINERALIZATION PROCESS OF QUARTZ H2O-CO2-NaCl SYSTEM

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010560360.2A CN111508564B (en) 2020-06-18 2020-06-18 System and method for quantifying the mineralization process of hydrothermal fluid in quartz H2O-CO2-NaCl system

Publications (2)

Publication Number Publication Date
CN111508564A CN111508564A (en) 2020-08-07
CN111508564B true CN111508564B (en) 2023-10-03

Family

ID=71878825

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010560360.2A Active CN111508564B (en) 2020-06-18 2020-06-18 System and method for quantifying the mineralization process of hydrothermal fluid in quartz H2O-CO2-NaCl system

Country Status (2)

Country Link
JP (1) JP2021196336A (en)
CN (1) CN111508564B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113539387A (en) * 2021-07-09 2021-10-22 西南石油大学 A Method for Predicting CO2 Solubility in NaCl Aqueous Solution Based on CPA Equation of State
CN115406880B (en) * 2022-10-19 2023-03-24 核工业北京地质研究院 A method for judging the reducibility of hydrothermal uranium ore-forming fluid
CN116359263A (en) * 2023-03-10 2023-06-30 昆明理工大学 A Discrimination Method of Keatite Vein Generation Based on Quartz Cathodoluminescence Characteristics
CN116665801B (en) * 2023-04-25 2026-02-06 东南大学 MgO-MgSO (magnesium oxide-zinc oxide) determination method4-MgCl2-H2Method for O quaternary system cement phase diagram
CN119982470B (en) * 2024-12-27 2025-10-31 中国石油化工股份有限公司 Supercritical CO2Automatic control method for compressor

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1045762A (en) * 1988-12-09 1990-10-03 巴特尔纪念研究院 Adopt control stoichiometry and granularity to produce the method for the submicron ceramic powder of uhligite mixture
US6393906B1 (en) * 2001-01-31 2002-05-28 Exxonmobil Upstream Research Company Method to evaluate the hydrocarbon potential of sedimentary basins from fluid inclusions
CN104632158A (en) * 2013-11-13 2015-05-20 中国石油化工股份有限公司 Oil well response stage dividing and distinguishing method under carbon dioxide miscible displacement condition
CN108344855A (en) * 2018-02-06 2018-07-31 长江大学 The determination method of Carbonate Reservoir Diagenetic Facies boundary point
CN109488261A (en) * 2018-12-18 2019-03-19 青岛海洋地质研究所 A kind of virtual construction method of ocean gas hydrate three phase dynamic
CN109752443A (en) * 2018-12-24 2019-05-14 核工业北京地质研究院 A method for determining the metallogenic causes of sandstone-type uranium deposits

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1045762A (en) * 1988-12-09 1990-10-03 巴特尔纪念研究院 Adopt control stoichiometry and granularity to produce the method for the submicron ceramic powder of uhligite mixture
US6393906B1 (en) * 2001-01-31 2002-05-28 Exxonmobil Upstream Research Company Method to evaluate the hydrocarbon potential of sedimentary basins from fluid inclusions
CN104632158A (en) * 2013-11-13 2015-05-20 中国石油化工股份有限公司 Oil well response stage dividing and distinguishing method under carbon dioxide miscible displacement condition
CN108344855A (en) * 2018-02-06 2018-07-31 长江大学 The determination method of Carbonate Reservoir Diagenetic Facies boundary point
CN109488261A (en) * 2018-12-18 2019-03-19 青岛海洋地质研究所 A kind of virtual construction method of ocean gas hydrate three phase dynamic
CN109752443A (en) * 2018-12-24 2019-05-14 核工业北京地质研究院 A method for determining the metallogenic causes of sandstone-type uranium deposits

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
H2O-CO2-NaCl体系石英溶解度模型(适用于高达1000℃、1.5GPa的高温高压环境);卫清等;岩石学报;第28卷(第08期);全文 *
胶东寺庄金矿床成因: 流体包裹体与石英溶解度证据;卫清等;岩石学报;第031卷(第04期);全文 *

Also Published As

Publication number Publication date
JP2021196336A (en) 2021-12-27
CN111508564A (en) 2020-08-07

Similar Documents

Publication Publication Date Title
CN111508564B (en) System and method for quantifying the mineralization process of hydrothermal fluid in quartz H2O-CO2-NaCl system
von Eynatten et al. Modelling compositional change: the example of chemical weathering of granitoid rocks
Mackay et al. Dynamic froth stability of copper flotation tailings
Shi et al. A novel method to improve depressants actions on calcite flotation
Romashev et al. Adaptive approach formation using machine vision technology to determine the parameters of enrichment products deposition
Taboada et al. Evaluation of the quality of a granite quarry
Birken et al. Quantification of impurities, including carbonates speciation for phosphates beneficiation by flotation
Michaux et al. Water-saving strategies in the mining industry–the potential of mineral processing simulators as a tool for their implementation
Lewandowski et al. Application of Raman spectroscopy to the flotation process of fluorite
Atashev Decomposition of Magnesite-Sparing Waste in Sulfuric Acid with a High Concentration: Empirical Modeling and Determination of Optimal Conditions
de Oliveira et al. Understanding regional groundwater flow using water-rock interaction processes and anthropic signatures in a fractured volcanic aquifer
Févotte et al. Batch cooling solution crystallization of ammonium oxalate in the presence of impurities: study of solubility, supersaturation, and steady-state inhibition
Tian et al. Investigation on the flotation separation of fluorite from celestite using a novel depressant: Sodium polynaphthalene formaldehyde sulfonate
Spier et al. Mineralogy and genesis of rare Al-phosphate minerals in weathered itabirite and iron ore from the Quadrilátero Ferrífero, Minas Gerais, Brazil
Kleiv A flotation-based approach to the quantification of pyrrhotite in concrete aggregates
Ren et al. A novel fracture-cavity reservoir outcrop geological knowledge base construction method considering parameter collection and processing, mutual transformation of data-knowledge, application and update
Salehi et al. Iron leaching from bauxite ore in hydrochloric acid using response surface methodology
Ross Key aspects of bench flotation as a geometallurgical characterization tool
Vyalov et al. Rare-earth elements in dictyonema shales of the Baltic sedimentary paleobasin
Nzerem et al. Study of the potential of sodium carbonate extracted from trona as a drilling fluid additive
Senol-Arslan et al. Diffuse-layer surface potentials of colemanites mined in Turkey
Žůrková et al. “Hungarian Mine Green”, a Semi‐Natural Copper Pigment from Banská Bystrica Region (Slovakia)‐Analytical Evidence and Laboratory Replication
Deniz et al. Interpretation using sink-float data by simple recovery equations for two different pumices
Mâamri et al. Characterization of the Oum El Khacheb phosphorites (South Tunisia) and enrichment of big rejections by grinding
CN119028468B (en) Quantitative simulation method for ore-forming element enrichment during metamorphism-anatexis

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
CB02 Change of applicant information

Address after: 100085, Anning Road, Beijing, Haidian District, No. 1

Applicant after: National natural disaster prevention and Control Research Institute Ministry of emergency management

Address before: 100000 Xisanqi, Haidian District, Beijing

Applicant before: THE INSTITUTE OF CRUSTAL DYNAMICS, CHINA EARTHQUAKE ADMINISTRATION

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant