CN111508564A - 定量化石英H2O-CO2-NaCl体系水热流体成矿过程的理论模型及方法 - Google Patents

定量化石英H2O-CO2-NaCl体系水热流体成矿过程的理论模型及方法 Download PDF

Info

Publication number
CN111508564A
CN111508564A CN202010560360.2A CN202010560360A CN111508564A CN 111508564 A CN111508564 A CN 111508564A CN 202010560360 A CN202010560360 A CN 202010560360A CN 111508564 A CN111508564 A CN 111508564A
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.)
Granted
Application number
CN202010560360.2A
Other languages
English (en)
Other versions
CN111508564B (zh
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.)
Institute of Crustal Dynamics of China Earthquake Administration
Original Assignee
Institute of Crustal Dynamics of China Earthquake Administration
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 Institute of Crustal Dynamics of China Earthquake Administration filed Critical Institute of Crustal Dynamics of China Earthquake Administration
Priority to CN202010560360.2A priority Critical patent/CN111508564B/zh
Publication of CN111508564A publication Critical patent/CN111508564A/zh
Priority to JP2020140828A priority patent/JP2021196336A/ja
Application granted granted Critical
Publication of CN111508564B publication Critical patent/CN111508564B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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

本发明公开了一种定量化石英H2O‑CO2‑NaCl体系水热流体成矿过程的理论模型及方法，理论模型用于研究不同成分成矿流体在矿区流体温度压力范围内的相边界曲线、石英等溶解度曲线、石英等溶解度曲线域边界曲线，包括：H2O‑CO2‑NaCl体系相边界计算模块，用于计算液相与气液混合相边界曲线；石英等溶解度计算模块，用于计算石英等溶解度曲线；石英等溶解度曲线域边界计算模块，用于计算石英等溶解度曲线域温度敏感区、压力敏感区和退行性区域边界曲线；上述各计算模块获得的曲线用于组合构建不同成矿流体成分的石英等溶解度相图，基于所述石英等溶解度相图可定量化石英H2O‑CO2‑NaCl体系水热流体动态演化乃至成矿过程，本发明可为理解基于石英水热过程的流体成矿过程提供新的视角。

Description

定量化石英H2O-CO2-NaCl体系水热流体成矿过程的理论模型 及方法
技术领域
本发明涉及矿床地球化学与计算地球化学交叉技术领域，特别是涉及一种定量化石英H2O-CO2-NaCl体系水热流体成矿过程的理论模型及方法。
背景技术
经查证，目前关于成矿领域的方法类专利仅限于成矿靶区圈定、成矿岩相带定位、成矿年代预测、成矿深度估算、成矿时空定位、成矿单元划分、成矿空间建模、成矿模式构建、成矿潜力计算等方面(CN104865613B、CN102243628A、CN110187387A、CN105785466A、CN110060173A、 CN107211585B、CN108573206A、CN108181669A、CN107765323A、CN109540929A、CN109270589B、CN107782878A)。涉及成矿流体或矿床成因方面的的方法专利也仅有6例，其揭示的主要为斑岩型矿床成矿流体氧化还原性质的判别(CN107655915B)及铀矿床成矿模式条件或成因等 (CN107576996A、CN109752443A、CN111044599A、CN109752443A、 CN106324700B)。
由此可见，上述现有技术对于成矿流体的研究主要为获得不同成矿阶段流体的物理化学性质，尚缺乏涉及流体定量动态演化乃至成矿的理论模型及方法，因此，如何能创设一种新的流体定量动态演化乃至成矿的理论模型及方法，实属当前重要研发课题之一。
发明内容
本发明要解决的技术问题是提供一种定量化石英H2O-CO2-NaCl体系水热流体成矿过程的理论模型及方法，可定量评价成矿流体动态演化乃至成矿的过程。
为解决上述技术问题，本发明首先提供了一种定量化石英H2O-CO2-NaCl 体系水热流体成矿过程的理论模型，所述理论模型用于研究不同成分成矿流体在矿区流体温度压力范围内的相边界曲线、石英等溶解度曲线、石英等溶解度曲线域边界曲线，包括：
H2O-CO2-NaCl体系相边界计算模块，用于计算获得固定流体成分及指定温度下相边界上某一点的温度压力坐标(T,P)，并以选定温度间隔重复本计算模块，得到一系列的在矿区成矿流体温度压力范围内的温度压力坐标，此即为液相与气液混合相边界曲线；
H2O-CO2-NaCl体系石英等溶解度计算模块，用于计算获得固定流体成分及指定石英溶解度时某一温度的温度压力坐标(T,P)，并以选定温度间隔重复本计算模块，得到一系列的在矿区成矿流体温度压力范围内的温度压力坐标，此即为石英等溶解度曲线；
H2O-CO2-NaCl体系石英等溶解度曲线域边界计算模块，用于获得固定流体成分条件下石英等溶解度曲线域边界的一系列温度压力坐标，此即为石英等溶解度曲线域温度敏感区、压力敏感区和退行性区域边界曲线；
上述各计算模块获得的相边界曲线、石英等溶解度曲线和石英等溶解度曲线域边界曲线用于组合构建不同成分成矿流体的石英等溶解度相图，基于所述石英等溶解度相图可定量化石英H2O-CO2-NaCl体系水热流体成矿过程。
作为本发明进一步地改进，所述理论模型的输入数据的来源为：
(1)对手标本和显微镜下石英脉穿插关系观察及对石英SEM-CL显微结构分析并照相确定矿床的成矿阶段及各成矿阶段的石英期次；
(2)进行石英流体包裹体激光拉曼光谱研究获得成矿流体的化学组成；
(3)开展不同期次石英中流体包裹体显微测温得到各期次石英中成矿流体的均一温度、均一压力及CO2和NaCl的含量；
(4)根据查阅文献获得的研究区矿床成矿流体演化的温度压力范围作为绘制的石英等溶解度相图温度和压力的上、下限范围，作为理论模型的矿区成矿流体温度压力范围的输入值；
(5)根据(3)中流体包裹体显微测温获得的矿区成矿流体CO2和NaCl 含量范围自由组合成矿流体的成分，以多种不同流体组成作为理论模型的成矿流体成分的输入值。
进一步地，所述H2O-CO2-NaCl体系相边界计算模块的计算流程为：
101、输入温度T(℃)，成矿流体化学组成包括二氧化碳CO2和氯化钠 NaCl，其中CO2含量以摩尔百分比(mol％)表示，NaCl的含量以质量百分比 (wt％)表示；
102、将成矿流体化学组分CO2，NaCl，H2O等值换算成摩尔分数形式并以xCO2,xNaCl及xH2O表示流体的成分，根据流体成分换算得到CO2在 H2O-NaCl体系中的溶解度以mCO2(mol/kg)表示；
103、对于压力P，给定其初始区间[P1,P2],利用CO2在H2O-NaCl体系溶解度模型，计算在压力为P1时CO2溶解度mCO2 P1，及在压力为P2时CO2溶解度mCO2 P2，且该区间边界满足条件：mCO2 P1<mCO2<mCO2 P2，压力初始范围[P1,P2] 应在其模型压力适用范围内；
104、令P＝(P1+P2)/2；
105、利用步骤103所述的CO2在H2O-NaCl体系溶解度模型，得到CO2溶解度计算值mCO2 cal；
106、当|mCO2 cal-mCO2|>10-6，且mCO2 cal<mCO2，则令P1＝P，反之，则令 P2＝P，并跳转至步骤104；当|mCO2 cal-mCO2|≤10-6时，跳转至步骤107；
107、得到指定成矿流体成分条件液相与气液混合相边界线上的点的坐标(T,P)，程序结束。
4、根据权利要求1所述的定量化石英H2O-CO2-NaCl体系水热流体成矿过程的理论模型，其特征在于，所述H2O-CO2-NaCl体系石英等溶解度计算模块的计算流程为：
201、输入温度T(℃)，成矿流体化学组成包括二氧化碳CO2和氯化钠 NaCl，其中CO2含量以摩尔百分比(mol％)表示，NaCl的含量以质量百分比 (wt％)表示，以及SiO2在H2O-CO2-NaCl体系的给定溶解度mSiO2以质量摩尔浓度(mol/kg)表示；
202、将成矿流体成分CO2，NaCl，H2O等值换算成摩尔分数形式，并以 xCO2，xNaCl及xH2O表示流体的成分；
203、对于压力P，给定其初始区间[P1,P2],利用SiO2在H2O-CO2-NaCl 体系溶解度模型，计算在压力为P1时石英溶解度mSiO2 P1，及在压力为P2时石英溶解度mSiO2 P2，且区间边界应满足条件：mSiO2 P1<mSiO2<mSiO2 P2，压力初始范围[P1,P2]应包含于模型压力适用范围；
204、令P＝(P1+P2)/2；
205、利用步骤203所述的溶解度模型，得到石英溶解度的计算值 mSiO2 cal；
206、当|mSiO2 cal-mSiO2|>10-9，且mSiO2 cal<mSiO2时，令P1＝P，反之，则令P2＝P，并跳转至步骤204；当|mSiO2 cal-mSiO2|≤10-9时，跳转至步骤207；
207、得到满足指定成矿流体成分及石英溶解度的点的坐标(T,P)，程序结束。
进一步地，所述H2O-CO2-NaCl体系石英等溶解度曲线域边界计算模块中温度与压力敏感区边界计算部分的计算流程为：
301、输入研究区SiO2在H2O-CO2-NaCl体系溶解度最大值mSiO2 max，以质量摩尔浓度(mol/kg)表示，成矿流体化学组成包括二氧化碳CO2和氯化钠NaCl，其中CO2含量以摩尔百分比(mol％)表示，NaCl的含量以质量百分比 (wt％)表示；
302、将成矿流体成分CO2，NaCl，H2O等值换算成摩尔分数形式，并以 xCO2，xNaCl及xH2O表示流体的成分；
303、根据石英等溶解度图，设置石英溶解度初始值mSiO2；
304、对于温度T，给定其初始区间[T1,T2],利用所述石英等溶解度计算模块计算温度为T1时固定溶解度条件下的压力P1，及温度为T2时固定溶解度条件下的压力P2，且区间边界应满足条件：P1>P>P2，温度初始范围[T1,T2] 应在其计算模块温度的适用范围内；
305、令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
306、根据SiO2在H2O-CO2-NaCl体系溶解度模型可知，在指定流体成分下石英溶解度为温度与压力的函数，令mSiO2＝m(T,P)，利用数值微分方法计算得到在当前石英溶解度条件下，点(T,P)处石英溶解度对温度与压力的偏导数
Figure BDA0002546082990000051
与
Figure BDA0002546082990000052
307、当
Figure BDA0002546082990000053
且
Figure BDA0002546082990000054
时，令T1＝T，反之，则令T2＝T，并跳转至步骤305；当
Figure BDA0002546082990000055
时，跳转至步骤308；
308、得到满足条件的边界上的点的坐标(T,P)；
309、令mSiO2＝mSiO2+计算步长值；
310、当mSiO2<mSiO2 max，跳转至步骤304；反之，则跳转至步骤311；
311、得到指定成矿流体成分条件下等溶解度曲线域温度与压力敏感区边界，程序结束。
进一步地，所述H2O-CO2-NaCl体系石英等溶解度曲线域边界计算模块中退行性区边界计算部分的计算流程为：
401、输入研究区SiO2在H2O-CO2-NaCl体系溶解度最大值mSiO2 max，并以质量摩尔浓度(mol/kg)表示，成矿流体成分包括二氧化碳CO2和氯化钠 NaCl，其中CO2含量以摩尔百分比(mol％)表示，NaCl的含量以质量百分比 (wt％)表示；
402、将成矿流体成分CO2，NaCl，H2O等值换算成摩尔分数形式，并以 xCO2,xNaCl及xH2O表示流体的成分；
403、令mSiO2＝mSiO2 max；
404、根据石英等溶解度图，设置温度T计算区间[T1,T2]；
405、令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
406、根据SiO2在H2O-CO2-NaCl体系溶解度模型，在指定流体成分下，石英溶解度为温度与压力的函数，令mSiO2＝m(T,P)，利用数值微分方法计算当前石英溶解度条件下点(T,P)处石英溶解度对温度的偏导数
Figure DEST_PATH_FDA0002546082980000055
407、若|T2-T1|<10-5，mSiO2＝mSiO2-计算步长值，反之则跳转至步骤 409；
408、当mSiO2>0时，跳转至步骤404，反之，当前成矿流体成分条件下，无等溶解度曲线域退行性区，程序结束；
409、当
Figure BDA0002546082990000062
且
Figure BDA0002546082990000063
时，令T1＝T，反之，则令T2＝T，并跳转至步骤405；当
Figure BDA0002546082990000064
时，得到退行性区边界的顶点，此时溶解度mSiO2 summit＝mSiO2，顶点坐标为(T,P)；
410、根据石英等溶解度图，初始化mSiO2；
411、根据石英等溶解度图，预计算当前溶解度条件下退行性区左端点，设置温度T计算区间[T1,T2]，所求温度T应包含左端点所对应温度；
412、令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
413、利用步骤406的方法计算当前石英溶解度条件下点(T,P)处石英溶解度对温度的偏导数
Figure BDA0002546082990000065
414、当
Figure BDA0002546082990000066
且
Figure BDA0002546082990000067
时，T1＝T，反之T2＝T，并跳转至步骤412；当
Figure BDA0002546082990000068
时，得到当前溶解度条件下，退行性区左端点；
415、根据石英等溶解度图，预计算当前溶解度条件下退行性区右端点，设置温度T计算区间[T1,T2]，所求温度T应包含右端点所对应温度；
416、令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
417、利用步骤406的方法计算当前石英溶解度条件下点(T,P)处石英溶解度对温度的偏导数
Figure BDA0002546082990000071
418、当
Figure BDA0002546082990000072
且
Figure BDA0002546082990000073
时，T2＝T，反之T1＝T，并跳转至步骤416；当
Figure BDA0002546082990000074
时，得到当前溶解度条件下，退行性区右端点；
419、mSiO2＝mSiO2+计算步长值；
420、若mSiO2<mSiO2 summit时，跳转至步骤411；反之跳转至步骤421；
421、得到固定成矿流体化成分条件下等溶解度曲线域的退行性区边界，程序结束。
基于定量化石英H2O-CO2-NaCl体系水热流体成矿过程的理论模型，本发明还提供了一种定量化石英H2O-CO2-NaCl体系水热流体成矿过程的方法，包括：
(一)实验阶段：
步骤(1)：在手标本及光学显微镜下观察研究区不同成矿阶段矿石样品中石英脉系的穿插关系，同时研究石英的扫描电镜-阴极发光SEM-CL显微结构，划分不同成矿阶段对应的精细石英期次；
步骤(2)：对研究区不同期次石英中的流体包裹体进行激光拉曼光谱研究，获得成矿流体化学成分为H2O-CO2-NaCl的水热流体；
步骤(3)：对研究区不同期次石英中的流体包裹体进行显微测温，获得成矿流体中气相组分CO2的含量，获得盐类物质NaCl的含量，获得不同期次石英中成矿流体的均一温度和均一压力范围；
步骤(4)：根据查阅文献获得研究区矿床成矿流体演化的温度压力范围作为绘制的石英等溶解度相图温度和压力上、下限；
步骤(5)：根据流体包裹体测温获得的研究区成矿流体CO2和NaCl含量范围自由组合成矿流体的成分，以多种不同流体组成作为计算阶段成矿流体成分的输入值；
(二)计算阶段：
步骤(6)：利用H2O-CO2-NaCl体系相边界计算模块，分别计算成矿流体成分为1…n时H2O-CO2-NaCl流体体系液相与气液混合相的边界上的点；
步骤(7)：在研究区矿床成矿流体演化的温度范围内，以任意温度步长重复步骤(6)计算，得到一系列的点(T,P)，此即为H2O-CO2-NaCl流体体系液相与气液混合相的相边界；
步骤(8)：利用H2O-CO2-NaCl体系石英等溶解度计算模块，分别计算成矿流体成分为1…n时H2O-CO2-NaCl流体体系石英等溶解度曲线上的点；
步骤(9)：在研究区矿床成矿流体演化的温度范围内，以任意温度步长重复步骤(8)计算，得到一系列的点(T,P)，此即为H2O-CO2-NaCl流体体系在指定石英溶解度时的等溶解度曲线；
步骤(10)：以任意石英溶解度步长重复步骤(8)及(9)计算，得到 H2O-CO2-NaCl流体体系不同溶解度条件下的石英等溶解度曲线；
步骤(11)：利用石英等溶解度曲线域边界计算模块中温度与压力敏感区边界计算部分，分别计算成矿流体成分为1…n时石英等溶解度曲线温度与压力敏感区的边界；
步骤(12)：利用石英等溶解度曲线域边界计算模块中退行性区边界计算部分，分别计算成矿流体成分为1…n时石英等溶解度曲线退行性区的边界；
(三)分析阶段：
步骤(13)：利用H2O-CO2-NaCl体系相边界计算模块、H2O-CO2-NaCl体系石英等溶解度计算模块及石英等溶解度曲线域边界计算模块获得的相边界曲线、石英等溶解度曲线和石英等溶解度曲线域边界，构建矿床成矿流体成分为1…n情况下的石英等溶解度相图；
步骤(14)：将实验阶段获得的不同期次石英中成矿流体的均一温度和均一压力范围投置于多种流体成分下的石英等溶解度相图中；
步骤(15)：根据石英等溶解度相图，分析其水热流体成矿过程。
通过采用上述技术方案，本发明至少具有以下优点：
本发明通过构建符合成矿流体物理化学性质的石英等溶解度相图的理论模型，利用引起石英溶解度变化的流体温度、压力、成分及相态改变所经历的路径，进而定量阐明成矿流体温度、压力、成分及相的动态演化，提供一种水热流体成矿过程的方法，弥补传统矿床学在研究流体成矿过程物理化学条件定量变化的描述，这种方法可为理解基于石英水热过程的流体成矿过程提供新的视角。
附图说明
上述仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，以下结合附图与具体实施方式对本发明作进一步的详细说明。
图1是定量化石英H2O-CO2-NaCl体系水热流体成矿过程的方法图；
图2是东风金矿床成矿期石英期次图；
图3是H2O-CO2-NaCl体系相边界计算模块的计算流程图；
图4是H2O-CO2-NaCl体系液相与气液混合相的相边界图；其中，实线为组分1时液相与气液混合相的相边界，虚线为组分2时液相与气液混合相的相边界；
图5是H2O-CO2-NaCl体系石英等溶解度计算模块的计算流程图；
图6是H2O-CO2-NaCl流体体系石英等溶解度曲线图；其中，实线为组分1时石英等溶解度曲线，虚线为组分2时石英等溶解度曲线；
图7是H2O-CO2-NaCl体系石英等溶解度曲线域边界计算模块的温度与压力敏感区边界计算部分的计算流程图；
图8是石英等溶解度曲线域温度与压力敏感区的边界图；其中，实线为组分1时温度与压力敏感区的边界，虚线为组分2时温度与压力敏感区的边界；
图9是石英等溶解度曲线域边界计算模块的退行性变化区边界计算部分的计算流程图；
图10是石英等溶解度曲线域退行性区的边界图；其中，实线为组分1 时退行性区的边界，虚线为组分2时退行性区的边界；
图11为东风金矿床石英等溶解度相图及不同期次石英的温度压力范围图；其中，实线为组分1时石英的等溶解度相图，实线为组分2时石英的等溶解度相图。
具体实施方式
本发明提供了一种基于石英等溶解度相图的涉及脉石石英的水热流体成矿过程的研究方法。为了便于对本发明的理解，对于如下术语进行定义：
a、石英等溶解度图：在表示H2O-CO2-NaCl水热系统温度-压力状态的状态图上，连接石英溶解度相同的状态的线；
b、石英等溶解度曲线域：石英等溶解度图涵盖的三个区域，具体包括温度敏感区、压力敏感区及退行性区；
c、压力敏感区：在石英等溶解度图中，压力变化引起的石英溶解度变化大于温度变化引起的石英溶解度变化的区域；
d、温度敏感区：在石英等溶解度图中，温度变化引起的石英溶解度变化大于压力变化引起的石英溶解度变化的区域；
e、退行性区：石英等溶解度图中石英溶解度变化与温度变化反相关的区域。
f、石英等溶解度相图：在表示H2O-CO2-NaCl水热系统温度-压力状态的状态图上，连接石英溶解度相同的状态的线，叠加液相与气液混合相的边界，同时区分等溶解度曲线域的温度敏感区、压力敏感区及退行性区。
本发明依托水热流体矿床石英脉样品薄片岩相学观察、流体包裹体激光拉曼光谱测试和显微测温等实验部分联合石英等溶解度相图计算部分，建立一种成矿流体在不同CO2和NaCl成分组成条件下，不同期次石英在石英等溶解度相图所处的温度和压力区间，通过各期次石英在石英等溶解度相图涉及的温度敏感区、压力敏感区、退行性区、流体不同相态区及不同流体组成引起的以上区域对应的温度和压力的演化等，定量化讨论不同期次石英乃至不同成矿阶段成矿流体可能经历的温度、压力和成分演化路径，从而达到定量评价石英H2O-CO2-NaCl体系水热流体的成矿过程。
如图1所示，矿物薄片实验部分共包含：(1)对手标本和显微镜下石英脉穿插关系观察及对石英SEM-CL显微结构分析并照相确定矿床的成矿阶段及各成矿阶段的石英期次；(2)进行石英流体包裹体激光拉曼光谱研究获得成矿流体的化学组成；(3)开展不同期次石英中流体包裹体显微测温得到各期次石英中成矿流体的均一温度、均一压力及CO2和NaCl的含量等三个步骤。
计算部分具体包括：(1)利用H2O-CO2-NaCl体系相边界计算模块获得指定温度下相边界上某一点的温度压力坐标(T,P)，并以选定温度间隔重复本计算模块，得到一系列的温度压力坐标，此即为液相与气液混合相边界；(2)利用H2O-CO2-NaCl体系石英等溶解度计算模块获得指定石英溶解度时某一温度的温度压力坐标(T,P)，并以选定温度间隔重复本计算模块，得到一系列的温度压力坐标，此即为石英等溶解度曲线；(3)采用 H2O-CO2-NaCl体系石英等溶解度曲线域边界计算模块获得指定温度下石英等溶解度曲线域边界上某一点的温度压力坐标(T,P)，并以选定温度间隔重复本计算模块，得到一系列的温度压力坐标，此即为石英等溶解度曲线域温度敏感区、压力敏感区和退行性区边界。这里需要说明是，计算部分应采用石英激光拉曼光谱测试和显微测温获得成矿流体的化学组成(CO2、 NaCl及H2O的含量均以摩尔分数xi表示)的实验结果，且成矿流体化学组成可根据研究需要任意设定，并根据水热流体矿床成矿演化的温度和压力上限设置本方法计算部分涉及的流体成分及温度压力范围。
成矿过程分析部分具体包括：(1)根据计算部分获得多组不同流体组成的石英等溶解度相图；(2)定量化分析石英水热流体温度、压力及成分演化对矿质沉淀的作用。这里需要说明的是流体成矿过程分析需要实验部分提供各期次石英成矿流体化学组成和温度压力范围。
下面通过具体应用例对本发明的技术方案进行展开描述。
实施例1
胶东矿集区是我国最为重要的原生金产区，其金矿石产出约占我国金产量的1/4。区内的众多金矿床均为典型的热液型金矿，石英是这些热液矿床最为重要的脉石矿物。对热液矿床石英流体包裹体及石英沉淀-溶解行为的研究，可为水热流体成矿过程乃至热液型矿床流体成矿机理提供新的视角和认识。本实施例以胶东金矿集区热液型金矿东风矿床为研究对象，开展具有不同流体组成矿床的石英包裹体显微测温及等溶解度相图实验模拟研究，定量评价石英H2O-CO2-NaCl体系水热流体的成矿过程。
本实施例提供了一种不同成矿流体组成的石英H2O-CO2-NaCl体系水热流体成矿过程的方法，该方法包括：
步骤(1)：在手标本及光学显微镜下观察东风金矿不同成矿阶段矿石样品中石英脉系的穿插关系，同时研究石英的扫描电镜-阴极发光(SEM-CL) 显微结构，划分不同成矿阶段对应的精细石英期次；
本实施例中，将东风金矿矿石样品分为V1、V2及V3成矿前、成矿期和成矿后三个阶段，石英期次细化为Qz1、Qz2及Qz3三期，详见图2。
其中，本实施例矿石样品石英Qz1的CL发光性较强，Qz2的CL发光性差，Qz3在CL下呈现出破碎的“蛛网状”结构，Qz1被Qz2和Qz3穿插。
步骤(2)：对东风金矿不同期次石英中的流体包裹体进行激光拉曼光谱研究，获得成矿流体的化学成分为H2O-CO2-NaCl的水热流体；
步骤(3)：对东风金矿不同期次石英中的流体包裹体进行显微测温，获得成矿流体中气相组分CO2的含量(以摩尔分数形式表示)，获得盐类物质NaCl的含量(以质量百分数形式表示)，获得不同期次石英中成矿流体的均一温度和均一压力范围；
其中，对于流体包裹体进行显微测温时，应尽量选择非均匀体系捕获的气液两相包裹体充填度最大和最小的包裹体进行，所获得的均一温度即为矿物的形成温度，不必进行压力校正；中间充填度的流体包裹体测温数据需进行压力校正。
东风金矿不同期次石英涵盖的成矿阶段及流体中CO2及NaCl的含量见表1。
表1.胶东东风金矿床石英期次划分及其对应的流体成分
Figure BDA0002546082990000131
步骤(4)：根据查阅文献获得胶东热液型金矿床成矿流体演化的大致温度压力范围作为本方法绘制的石英等溶解度相图温度和压力上，下限，这里选择温度范围为100-600℃，压力范围为0-4000bar；
步骤(5)：根据流体包裹体显微测温获得的东风金矿成矿流体CO2和 NaCl含量范围(10-15％molCO2，0.4-11.1％wtNaCl)自由组合成矿流体的成分。这里我们选择CO2摩尔百分比为10，NaCl质量百分比为5(组分 1)及CO2摩尔百分比为15，NaCl质量百分比为5(组分2)的两种流体组成作为计算部分成矿流体成分的输入值；
步骤(6)：配合图3所示，利用H2O-CO2-NaCl体系相边界计算模块，计算组分1时H2O-CO2-NaCl流体体系液相与气液混合相的边界上的点：
①输入温度100℃，CO2摩尔百分数10，NaCl质量百分数5；
②将成矿流体的化学组成换算成以摩尔分数形式表示：xCO2＝0.1, xNaCl＝0.0183,xH2O＝0.8817，换算得到CO2的溶解度mCO2＝6.30mol/kg；
③根据Mao等人文献(Mao S.D.,Zhang D.H.,Li Y.Q.et a.,An improved modelfor calculating CO2 solubility in aqueous NaCl solutions and the applicationto CO2–H2O–NaCl fluid inclusions， Chemical Geology,2013)所述CO2在H2O-NaCl体系溶解度模型，压力范围设置初始区间端点P1＝0，P2＝4000，可满足条件：mCO2 P1<mCO2<mCO2 P2；
④令P＝(P1+P2)/2；
⑤利用CO2在H2O-NaCl体系溶解度模型，计算当前压力P时，CO2溶解度mCO2 cal；
⑥当|mCO2 cal-mCO2|>10-6，且mCO2 cal<mCO2时，令P1＝P，反之，则令P2＝P，并跳转至步骤④；当该值|mCO2 cal-mCO2|≤10-6时，P＝2812.07bar即为所求压力。
计算组分2时H2O-CO2-NaCl流体体系液相与气液混合相的边界上的点：
①输入温度100℃，CO2摩尔百分数15，NaCl质量百分数5；
②将成矿流体的化学组成换算成以摩尔分数形式表示：xCO2＝0.15, xNaCl＝0.0194,xH2O＝0.8306，换算得到CO2的溶解度mCO2＝10.03mol/kg；
③根据Mao等人文献(Mao S.D.,Zhang D.H.,Li Y.Q.et a.,An improved modelfor calculating CO2 solubility in aqueous NaCl solutions and the applicationto CO2–H2O–NaCl fluid inclusions， Chemical Geology,2013)所述CO2在H2O-NaCl体系溶解度模型，压力范围设置初始区间端点P1＝0，P2＝4000，可满足条件：mCO2 P1<mCO2<mCO2 P2；
④令P＝(P1+P2)/2；
⑤利用CO2在H2O-NaCl体系溶解度模型公式，计算当前压力P时，CO2溶解度mCO2 cal；
⑥当|mCO2 cal-mCO2|>10-6，且mCO2 cal<mCO2时，令P1＝P，反之，则令P2＝P，并跳转至步骤④；当该值|mCO2 cal-mCO2|≤10-6时，P＝3040.01bar即为所求压力。
步骤(7)：在100-600℃温度范围内，以任意温度步长(如以温度1℃为间隔)重复步骤(6)计算，得到一系列的点(T,P)，此即为H2O-CO2-NaCl 流体体系液相与气液混合相的相边界(图4)；
步骤(8)：配合图5所示，利用H2O-CO2-NaCl体系石英等溶解度计算模块，计算组分1时H2O-CO2-NaCl流体体系石英等溶解度曲线上的点：
①输入温度150℃，CO2摩尔百分数10，NaCl质量百分数5，SiO2在 H2O-CO2-NaCl体系的溶解度mSiO2＝0.002mol/kg；
②将成矿流体的化学组成换算成以摩尔分数形式表示：xCO2＝0.1, xNaCl＝0.0183,xH2O＝0.8817；
③根据利用卫清等人文献(卫清，段振豪，毛世德.H2O-CO2-NaCl体系石英溶解度模型(适用于高达1000℃、1.5GPa的高温高压环境)，岩石学报，2012)中所述SiO2在H2O-CO2-NaCl体系溶解度(简称：石英溶解度) 模型，设置初始区间端点P1＝0，P2＝4000，可满足条件：mSiO2 P1<mCO2<mSiO2 P2；
④令P＝(P1+P2)/2；
⑤利用石英溶解度模型，计算当前压力P时，SiO2溶解度mSiO2 cal；
⑥当|mSiO2 cal-mSiO2|＞10-9，且mSiO2 cal<mSiO2时，令P1＝P，反之，则令 P2＝P，并跳转至步骤④；当|mSiO2 cal-mSiO2|≤10-9时，P＝3408.33bar即为所求压力。
计算组分2时H2O-CO2-NaCl流体体系石英等溶解度曲线上的点：
①输入温度160℃，CO2摩尔百分数15，NaCl质量百分数5，SiO2在 H2O-CO2-NaCl体系的溶解度mSiO2＝0.002mol/kg；
②将成矿流体的化学组成换算成以摩尔分数形式表示：xCO2＝0.15, xNaCl＝0.0194,xH2O＝0.8306；
③根据利用卫清等人文献(卫清，段振豪，毛世德.H2O-CO2-NaCl体系石英溶解度模型(适用于高达1000℃、1.5GPa的高温高压环境)，岩石学报，2012)中所述SiO2在H2O-CO2-NaCl体系溶解度(简称：石英溶解度) 模型，设置初始区间端点P1＝0，P2＝4000，可满足条件：mSiO2 P1<mCO2<mSiO2 P2；
④令P＝(P1+P2)/2；
⑤利用石英溶解度模型，计算当前压力P时，SiO2溶解度mSiO2 cal；
⑥当|mSiO2 cal-mSiO2|＞10-9，且mSiO2 cal<mSiO2时，令P1＝P，反之，则令P2＝P，并跳转至步骤④；当|mSiO2 cal-mSiO2|≤10-9时，P＝3831.83bar即为所求压力。
步骤(9)：在100-600℃温度范围内，以任意温度步长(如以温度1℃为间隔)重复步骤(8)计算，得到一系列的点(T,P)，此即为H2O-CO2-NaCl 流体体系石英溶解度为0.002mol/kg时的等溶解度曲线；
步骤(10)：以任意石英溶解度步长(如以溶解度0.001mol/kg为间隔) 重复步骤(8)及(9)计算，得到H2O-CO2-NaCl流体体系不同溶解度条件下的石英等溶解度曲线(图6)；
步骤(11)：配合图7所示，利用石英等溶解度曲线域边界计算模块- 温度与压力敏感区边界计算部分，计算组分1时石英等溶解度曲线温度与压力敏感区的边界：
①输入等溶解度相图温度与压力范围内的SiO2在H2O-CO2-NaCl体系的溶解度(石英溶解度)最大值mSiO2 max＝0.087mol/kg，CO2摩尔百分数10，NaCl 质量百分数5；
②将成矿流体的化学组成换算成以摩尔分数形式表示：xCO2＝0.1, xNaCl＝0.0183,xH2O＝0.8817；
③令石英溶解度mSiO2＝0.002molg/kg；
④根据石英等溶解度图，设置初始区间端点T1＝175，T2＝600；
⑤令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
⑥根据卫清等人文献(卫清，段振豪，毛世德.H2O-CO2-NaCl体系石英溶解度模型(适用于高达1000℃、1.5GPa的高温高压环境)，岩石学报， 2012)中所述SiO2在H2O-CO2-NaCl体系溶解度模型可知，在指定流体成分下石英溶解度为温度与压力的函数，令mSiO2＝m(T,P)。利用数值微分方法计算得到在当前石英溶解度条件下，点(T,P)处石英溶解度对温度与压力的偏导数
Figure BDA0002546082990000161
与
Figure BDA0002546082990000162
⑦当
Figure BDA0002546082990000163
且
Figure BDA0002546082990000164
时，令T1＝T，反之，则令T2＝T，并跳转至步骤⑤；当
Figure BDA0002546082990000171
时， T＝183.00℃，P＝14.12bar即为满足条件的边界曲线上的点的坐标(T,P)；
⑧令mSiO2＝mSiO2+10-4；
⑨若mSiO2<mSiO2 max,则跳转至步骤④，否则结束程序，所求得的一系列点，即为石英等溶解度曲线温度与压力敏感区的边界。
计算组分2时石英等溶解度曲线温度与压力敏感区的边界：
①输入等溶解度相图温度与压力范围内的SiO2在H2O-CO2-NaCl体系的溶解度(石英溶解度)最大值mSiO2 max＝0.068mol/kg，CO2摩尔百分数15，NaCl 质量百分数5；
②将成矿流体的化学组成换算成以摩尔分数形式表示：xCO2＝0.15, xNaCl＝0.0194,xH2O＝0.8306；
③令石英溶解度mSiO2＝0.002molg/kg；
④根据石英等溶解度图，设置初始区间端点T1＝190，T2＝600；
⑤令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
⑥根据卫清等人文献(卫清，段振豪，毛世德.H2O-CO2-NaCl体系石英溶解度模型(适用于高达1000℃、1.5GPa的高温高压环境)，岩石学报， 2012)中所述SiO2在H2O-CO2-NaCl体系溶解度模型可知，在指定流体成分下石英溶解度为温度与压力的函数，令mSiO2＝m(T,P)。利用数值微分方法计算得到在当前石英溶解度条件下，点(T,P)处石英溶解度对温度与压力的偏导数
Figure BDA0002546082990000172
与
Figure BDA0002546082990000173
⑦当
Figure BDA0002546082990000174
且
Figure BDA0002546082990000175
时，令T1＝T，反之，则令T2＝T，并跳转至步骤⑤；当
Figure BDA0002546082990000176
时， T＝222.50℃，P＝26.76bar即为满足条件的边界曲线上的点的坐标(T,P)；
⑧令mSiO2＝mSiO2+10-4；
⑨若mSiO2<mSiO2 max,则跳转至步骤④，否则结束程序，所求得的一系列点，即为石英等溶解度曲线域温度与压力敏感区的边界(图8)。
步骤(12)：配合图9所示，利用石英等溶解度曲线域边界计算模块- 退行性区边界计算部分，计算组分1时石英等溶解度曲线退行性区的边界：
1、输入等溶解度相图温度与压力范围内的SiO2在H2O-CO2-NaCl体系的溶解度(石英溶解度)最大值mSiO2 max＝0.02mol/kg，CO2摩尔百分数10，NaCl 质量百分数5；
2、将成矿流体的化学组成换算成以摩尔分数形式表示：xCO2＝0.1, xNaCl＝0.0183,xH2O＝0.8817；
3、令石英溶解度mSiO2＝mSiO2 max
4、根据石英等溶解度图，设置初始区间端点T1＝450，T2＝600；
5、令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
6、根据卫清等人文献(卫清，段振豪，毛世德.H2O-CO2-NaCl体系石英溶解度模型(适用于高达1000℃、1.5GPa的高温高压环境)，岩石学报， 2012)中所述SiO2在H2O-CO2-NaCl体系溶解度模型可知，在指定流体成分下石英溶解度为温度与压力的函数，令mSiO2＝m(T,P)。利用数值微分方法计算得到在当前石英溶解度条件下，点(T,P)处石英溶解度对温度与压力的偏导数
Figure BDA0002546082990000181
7、若|T2-T1|<10-5，mSiO2＝mSiO2-10-4(这里10-4为计算步长值，步长值可根据实际计算条件调整)，反之则跳转至9；
8、当mSiO2>0时，跳转至4，反之，当前成矿流体成分条件下，无等溶解度曲线退行性区，程序结束；
9、当
Figure BDA0002546082990000182
且
Figure BDA0002546082990000183
时，令T1＝T，反之，则令T2＝T，并跳转至步骤⑤；当
Figure BDA0002546082990000184
时，mSiO2 summit＝0.0121mol/kg，且 T＝568.75℃，P＝649.79bar即为满足条件的边界曲线上的顶点的坐标；
10、令mSiO2＝0.002；
11、计算当前溶解度条件下，退行性区边界左端点,根据石英等溶解度图，设置初始区间端点T1＝180，T2＝200；
12、令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
13、利用⑥的方法得到在当前石英溶解度条件下，点(T,P)处石英溶解度对温度与压力的偏导数
Figure BDA0002546082990000191
14、当
Figure BDA0002546082990000192
且
Figure BDA0002546082990000193
时，T1＝T，反之T2＝T，跳转至12；当
Figure BDA0002546082990000194
时，T＝184.00℃，P＝10.98bar即为在当前成矿流体组成下退行性区边界左端点的坐标。
15、计算当前溶解度条件下，退行性区边界右端点，根据石英等溶解度图，设置初始区间端点T1＝450，T2＝600；
16、令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
17、利用⑥的方法计算得到在当前石英溶解度条件下，点(T,P)处石英溶解度对温度与压力的偏导数
Figure BDA0002546082990000195
18、当
Figure BDA0002546082990000196
且
Figure BDA0002546082990000197
时，T2＝T，反之T1＝T，跳转至16；当
Figure BDA0002546082990000198
时，T＝523.05℃，P＝350.65bar为在当前成矿流体组成下石英溶解度为0.002mol/kg时退行性区边界右端点的坐标；
19、mSiO2＝mSiO2+10-4；
20、若mSiO2<mSiO2 max,则跳转至11，若mSiO2≥mSiO2 max时，结束该程序，所求得的一系列点，即为石英等溶解度曲线退行性区的边界(图10)。
计算组分2时石英等溶解度曲线退行性区的边界：
1、输入等溶解度相图温度与压力范围内的SiO2在H2O-CO2-NaCl体系的溶解度(石英溶解度)最大值mSiO2 max＝0.01mol/kg，CO2摩尔百分数15，NaCl 质量百分数5；
2、将成矿流体的化学组成换算成以摩尔分数形式表示：xCO2＝0.15, xNaCl＝0.0194,xH2O＝0.8306；
3、令石英溶解度mSiO2＝mSiO2 max；
4、根据石英等溶解度图，设置初始区间端点T1＝450，T2＝600；
5、令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
6、根据卫清等人文献(卫清，段振豪，毛世德.H2O-CO2-NaCl体系石英溶解度模型(适用于高达1000℃、1.5GPa的高温高压环境)，岩石学报， 2012)中所述SiO2在H2O-CO2-NaCl体系溶解度模型可知，在指定流体成分下石英溶解度为温度与压力的函数，令mSiO2＝m(T,P)。利用数值微分方法计算得到在当前石英溶解度条件下，点(T,P)处石英溶解度对温度与压力的偏导数
Figure 621624DEST_PATH_FDA0002546082980000055
7、若|T2-T1|<10-5，mSiO2＝mSiO2-10-4(这里10-4为计算步长值，步长值可根据实际计算条件调整)，反之则跳转至9；
8、当mSiO2>0时，跳转至4，反之，当前成矿流体成分条件下，无等溶解度曲线退行性区，程序结束；
9、当
Figure BDA0002546082990000202
且
Figure BDA0002546082990000203
时，令T1＝T，反之，则令T2＝T，并跳转至5；当
Figure BDA0002546082990000204
时，mSiO2＝0.0057mol/kg，且T＝525.00 ℃，P＝548.22bar即为满足条件的边界曲线上的顶点的坐标；
10、令mSiO2＝0.002；
11、算当前溶解度条件下，退行性区边界左端点，根据石英等溶解度图，设置初始区间端点T1＝180，T2＝200；
12、令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
13、利用⑥的方法得到在当前石英溶解度条件下，点(T,P)处石英溶解度对温度与压力的偏导数
Figure BDA0002546082990000205
14、当
Figure BDA0002546082990000206
且
Figure BDA0002546082990000207
时，T1＝T，反之T2＝T，跳转至12；当
Figure BDA0002546082990000208
时，T＝223.05℃，P＝24.58bar即为在当前成矿流体组成下退行性区边界左端点的坐标。
15、计算当前溶解度条件下，退行性区边界右端点，根据石英等溶解度图，设置初始区间端点T1＝450，T2＝550；
16、令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
17、利用⑥的方法计算得到在当前石英溶解度条件下，点(T,P)处石英溶解度对温度与压力的偏导数
Figure BDA0002546082990000211
18、当
Figure BDA0002546082990000212
且
Figure BDA0002546082990000213
时，T2＝T，反之T1＝T，跳转至16；当
Figure BDA0002546082990000214
时，T＝528.83℃，P＝385.25bar为在当前成矿流体组成下退行性区边界右端点；
19、mSiO2＝mSiO2+10-4；
20、若mSiO2<mSiO2 max,则跳转至11，若mSiO2≥mSiO2 max时，结束该程序，所求得的一系列点，即为石英等溶解度曲线域退行性区的边界(图10)。
步骤(13)：利用H2O-CO2-NaCl体系相边界计算模块、石英等溶解度计算模块及石英等溶解度曲线域边界计算模块获得的相边界曲线、石英等溶解度曲线和石英等溶解度曲线域边界，构建东风金矿床流体组成分别为组分1和组分2两种情况下的石英等溶解度相图；
步骤(14)：将实验部分获得的东风金矿不同期次石英中成矿流体的均一温度和均一压力范围投置于两种流体组成下的石英等溶解度相图中，如图11所示；
步骤(15)：根据东风金矿床石英等溶解度相图(图11)，分析其水热流体成矿过程。
本实施例中，当东风金矿成矿流体组分为组分1和组分2两种情况时， Qz1均落在石英等溶解度相图液相区的温度敏感区；当成矿流体组分为组分 1时，Qz2落在石英等溶解度相图液相区的温度敏感区，但当成矿流体组分为组分2时，Qz2则横跨石英等溶解度相图液相和气液混合相的温度敏感区；当成矿流体组分为组分1时，Qz3横跨石英等溶解度相图气液两相区的温度敏感区，压力敏感区和退行性区三个区域，但当成矿流体组分为组分2时，Qz3则只横跨石英等溶解度相图气液混合相的压力敏感区和退行性区。
就本实施例东风金矿流体演化，对某假设初始流体对应的温度压力 (550-600℃，3200-3600bar)范围而言，从初始流体沉淀出Qz1，可由成矿流体在液相区由600℃降温至320℃而形成；石英由Qz1演变为Qz2，可由成矿流体在液相区由370℃降温至300℃，或是由于流体中CO2摩尔分数由百分之十增加至百分之十五造成气液相分离，亦或是二者叠加作用而形成；石英由Qz2演化至Qz3，可由成矿流体的气液相分离造成的成矿流体压力由1600bar降低至200bar，或是石英的退行性溶解行为，亦或是二者叠加作用而形成。
综合以上Qz1，Qz2和Qz3的形成过程，根据东风金矿石英H2O-CO2-NaCl 体系等溶解度相图可知，东风金矿成矿流体首先在液相区经历了降温(约从600℃降低至320℃)，卸载出少量的金属物质，此阶段为东风金矿成矿前阶段(V1阶段)；进而，随着成矿流体温度的进一步下降(从约370℃降至300℃)以及流体中非极性气体CO2含量增高(体系占比摩尔百分数由十增加至十五)造成的流体相分离过程，成矿流体中卸载出大量的金属物质，此阶段为东风金矿最为主要的成矿阶段(V2阶段)；随着成矿流体温度和压力的进一步降低(温度由350℃降至270℃，压力由1600bar下降至约 200bar)，成矿流体经历了显著的相分离阶段(V3阶段)，压力的急剧释放使得金属物质得以沉淀析出，另外，在石英退行性区较窄的温度区间内 (280℃-330℃)，金属物质也会随着温度的小幅增加而少量析出。
综合以上，即为东风金矿定量化石英H2O-CO2-NaCl体系水热流体的成矿过程。
以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，本领域技术人员利用上述揭示的技术内容做出些许简单修改、等同变化或修饰，均落在本发明的保护范围内。

Claims (7)

1.定量化石英H2O-CO2-NaCl体系水热流体成矿过程的理论模型，其特征在于，所述理论模型用于研究不同成分成矿流体在矿区流体温度压力范围内的相边界曲线、石英等溶解度曲线、石英等溶解度曲线域边界曲线，包括：
H2O-CO2-NaCl体系相边界计算模块，用于计算获得固定流体成分及指定温度下相边界上某一点的温度压力坐标(T,P)，并以选定温度间隔重复本计算模块，得到一系列的在矿区成矿流体温度压力范围内的温度压力坐标，此即为液相与气液混合相边界曲线；
H2O-CO2-NaCl体系石英等溶解度计算模块，用于计算获得固定流体成分及指定石英溶解度时某一温度的温度压力坐标(T,P)，并以选定温度间隔重复本计算模块，得到一系列的在矿区成矿流体温度压力范围内指定石英溶解度条件下的温度压力坐标，此即为石英等溶解度曲线；
H2O-CO2-NaCl体系石英等溶解度曲线域边界计算模块，用于获得固定流体成分条件下石英等溶解度曲线域边界的一系列温度压力坐标，此即为石英等溶解度曲线域温度敏感区、压力敏感区和退行性区域边界曲线；
上述各计算模块获得的相边界曲线、石英等溶解度曲线和石英等溶解度曲线域边界曲线用于组合构建不同成分成矿流体的石英等溶解度相图，基于所述石英等溶解度相图可定量化石英H2O-CO2-NaCl体系水热流体成矿过程。
2.根据权利要求1所述的定量化石英H2O-CO2-NaCl体系水热流体成矿过程的理论模型，其特征在于，所述理论模型的输入数据的来源为：
(1)对手标本和显微镜下石英脉穿插关系观察及对石英SEM-CL显微结构分析并照相确定矿床的成矿阶段及各成矿阶段的石英期次；
(2)进行石英流体包裹体激光拉曼光谱研究获得成矿流体的化学组成；
(3)开展不同期次石英中流体包裹体显微测温得到各期次石英中成矿流体的均一温度、均一压力及CO2和NaCl的含量；
(4)根据查阅文献获得的研究区矿床成矿流体演化的温度压力范围作为绘制的石英等溶解度相图温度和压力的上、下限范围，作为理论模型的矿区成矿流体温度压力范围的输入值；
(5)根据(3)中流体包裹体显微测温获得的矿区成矿流体CO2和NaCl含量范围自由组合成矿流体的成分，以多种不同流体组成作为理论模型的成矿流体成分的输入值。
3.根据权利要求1所述的定量化石英H2O-CO2-NaCl体系水热流体成矿过程的理论模型，其特征在于，所述H2O-CO2-NaCl体系相边界计算模块的计算流程为：
101、输入温度T(℃)，成矿流体化学组成包括二氧化碳CO2和氯化钠NaCl，其中CO2含量以摩尔百分比(mol％)表示，NaCl的含量以质量百分比(wt％)表示；
102、将成矿流体化学组分CO2，NaCl，H2O等值换算成摩尔分数形式并以xCO2,xNaCl及xH2O表示流体的成分，根据流体成分换算得到CO2在H2O-NaCl体系中的溶解度以mCO2(mol/kg)表示；
103、对于压力P，给定其初始区间[P1,P2],利用CO2在H2O-NaCl体系溶解度模型，计算在压力为P1时CO2溶解度mCO2 P1，及在压力为P2时CO2溶解度mCO2 P2，且该区间边界满足条件：mCO2 P1<mCO2<mCO2 P2，压力初始范围[P1,P2]应在其模型压力适用范围内；
104、令P＝(P1+P2)/2；
105、利用步骤103所述的CO2在H2O-NaCl体系溶解度模型，得到CO2溶解度计算值mCO2 cal；
106、当|mCO2 cal-mCO2|>10-6，且mCO2 cal<mCO2，则令P1＝P，反之，则令P2＝P，并跳转至步骤104；当|mCO2 cal-mCO2|≤10-6时，跳转至步骤107；
107、得到指定成矿流体成分条件液相与气液混合相边界线上的点的坐标(T,P)，程序结束。
4.根据权利要求1所述的定量化石英H2O-CO2-NaCl体系水热流体成矿过程的理论模型，其特征在于，所述H2O-CO2-NaCl体系石英等溶解度计算模块的计算流程为：
201、输入温度T(℃)，成矿流体化学组成包括二氧化碳CO2和氯化钠NaCl，其中CO2含量以摩尔百分比(mol％)表示，NaCl的含量以质量百分比(wt％)表示，以及SiO2在H2O-CO2-NaCl体系的给定溶解度mSiO2以质量摩尔浓度(mol/kg)表示；
202、将成矿流体成分CO2，NaCl，H2O等值换算成摩尔分数形式，并以xCO2，xNaCl及xH2O表示流体的成分；
203、对于压力P，给定其初始区间[P1,P2],利用SiO2在H2O-CO2-NaCl体系溶解度模型，计算在压力为P1时石英溶解度mSiO2 P1，及在压力为P2时石英溶解度mSiO2 P2，且区间边界应满足条件：mSiO2 P1<mSiO2<mSiO2 P2，压力初始范围[P1,P2]应包含于模型压力适用范围；
204、令P＝(P1+P2)/2；
205、利用步骤203所述的溶解度模型，得到石英溶解度的计算值mSiO2 cal；
206、当|mSiO2 cal-mSiO2|>10-9，且mSiO2 cal<mSiO2时，令P1＝P，反之，则令P2＝P，并跳转至步骤204；当|mSiO2 cal-mSiO2|≤10-9时，跳转至步骤207；
207、得到满足指定成矿流体成分及石英溶解度的点的坐标(T,P)，程序结束。
5.根据权利要求1-4任一项所述的定量化石英H2O-CO2-NaCl体系水热流体成矿过程的理论模型，其特征在于，所述H2O-CO2-NaCl体系石英等溶解度曲线域边界计算模块中温度与压力敏感区边界计算部分的计算流程为：
301、输入研究区SiO2在H2O-CO2-NaCl体系溶解度最大值mSiO2 max，以质量摩尔浓度(mol/kg)表示，成矿流体化学组成包括二氧化碳CO2和氯化钠NaCl，其中CO2含量以摩尔百分比(mol％)表示，NaCl的含量以质量百分比(wt％)表示；
302、将成矿流体成分CO2，NaCl，H2O等值换算成摩尔分数形式，并以xCO2，xNaCl及xH2O表示流体的成分；
303、根据石英等溶解度图，设置石英溶解度初始值mSiO2；
304、对于温度T，给定其初始区间[T1,T2],利用所述石英等溶解度计算模块计算温度为T1时固定溶解度条件下的压力P1，及温度为T2时固定溶解度条件下的压力P2，且区间边界应满足条件：P1>P>P2，温度初始范围[T1,T2]应在其计算模块温度的适用范围内；
305、令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
306、根据SiO2在H2O-CO2-NaCl体系溶解度模型可知，在指定流体成分下石英溶解度为温度与压力的函数，令mSiO2＝m(T,P)，利用数值微分方法计算得到在当前石英溶解度条件下，点(T,P)处石英溶解度对温度与压力的偏导数
Figure FDA0002546082980000041
与
Figure FDA0002546082980000042
307、当
Figure FDA0002546082980000043
且
Figure FDA0002546082980000044
时，令T1＝T，反之，则令T2＝T，并跳转至步骤305；当
Figure FDA0002546082980000045
时，跳转至步骤308；
308、得到满足条件的边界上的点的坐标(T,P)；
309、令mSiO2＝mSiO2+计算步长值；
310、当mSiO2<mSiO2 max，跳转至步骤304；反之，则跳转至步骤311；
311、得到指定成矿流体成分条件下等溶解度曲线域温度与压力敏感区边界，程序结束。
6.根据权利要求1-4任一项所述的定量化石英H2O-CO2-NaCl体系水热流体成矿过程的理论模型，其特征在于，所述H2O-CO2-NaCl体系石英等溶解度曲线域边界计算模块中退行性区边界计算部分的计算流程为：
401、输入研究区SiO2在H2O-CO2-NaCl体系溶解度最大值mSiO2 max，并以质量摩尔浓度(mol/kg)表示，成矿流体成分包括二氧化碳CO2和氯化钠NaCl，其中CO2含量以摩尔百分比(mol％)表示，NaCl的含量以质量百分比(wt％)表示；
402、将成矿流体成分CO2，NaCl，H2O等值换算成摩尔分数形式，并以xCO2,xNaCl及xH2O表示流体的成分；
403、令mSiO2＝mSiO2 max；
404、根据石英等溶解度图，设置温度T计算区间[T1,T2]；
405、令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
406、根据SiO2在H2O-CO2-NaCl体系溶解度模型，在指定流体成分下，石英溶解度为温度与压力的函数，令mSiO2＝m(T,P)，利用数值微分方法计算当前石英溶解度条件下点(T,P)处石英溶解度对温度的偏导数
Figure FDA0002546082980000055
407、若|T2-T1|<10-5，mSiO2＝mSiO2-计算步长值，反之则跳转至步骤409；
408、当mSiO2>0时，跳转至步骤404，反之，当前成矿流体成分条件下，无等溶解度曲线域退行性区，程序结束；
409、当
Figure FDA0002546082980000052
且
Figure FDA0002546082980000053
时，令T1＝T，反之，则令T2＝T，并跳转至步骤405；当
Figure FDA0002546082980000054
时，得到退行性区边界的顶点，此时溶解度mSiO2 summit＝mSiO2，顶点坐标为(T,P)；
410、根据石英等溶解度图，初始化mSiO2；
411、根据石英等溶解度图，预计算当前溶解度条件下退行性区左端点，设置温度T计算区间[T1,T2]，所求温度T应包含左端点所对应温度；
412、令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
413、利用步骤406的方法计算当前石英溶解度条件下点(T,P)处石英溶解度对温度的偏导数
Figure FDA0002546082980000055
414、当
Figure FDA0002546082980000056
且
Figure FDA0002546082980000057
时，T1＝T，反之T2＝T，并跳转至步骤412，当
Figure FDA0002546082980000058
时，得到当前溶解度条件下，退行性区左端点；
415、根据石英等溶解度图，预计算当前溶解度条件下退行性区右端点，设置温度T计算区间[T1,T2]，所求温度T应包含右端点所对应温度；
416、令T＝(T1+T2)/2，利用石英等溶解度计算模块，计算得到压力P；
417、利用步骤406的方法计算当前石英溶解度条件下点(T,P)处石英溶解度对温度的偏导数
Figure FDA0002546082980000061
418、当
Figure FDA0002546082980000062
且
Figure FDA0002546082980000063
时，T2＝T，反之T1＝T，并跳转至步骤416；当
Figure FDA0002546082980000064
时，得到当前溶解度条件下，退行性区右端点；
419、mSiO2＝mSiO2+计算步长值；
420、若mSiO2<mSiO2 summit时，跳转至步骤411；反之跳转至步骤421；
421、得到固定成矿流体化成分条件下等溶解度曲线域的退行性区边界，程序结束。
7.一种定量化石英H2O-CO2-NaCl体系水热流体成矿过程的方法，其特征在于，利用权利要求1-6任一项所述的定量化石英H2O-CO2-NaCl体系水热流体成矿过程的理论模型，进行计算分析，包括：
(一)实验阶段：
步骤(1)：在手标本及光学显微镜下观察研究区不同成矿阶段矿石样品中石英脉系的穿插关系，同时研究石英的扫描电镜-阴极发光SEM-CL显微结构，划分不同成矿阶段对应的精细石英期次；
步骤(2)：对研究区不同期次石英中的流体包裹体进行激光拉曼光谱研究，获得成矿流体化学成分为H2O-CO2-NaCl的水热流体；
步骤(3)：对研究区不同期次石英中的流体包裹体进行显微测温，获得成矿流体中气相组分CO2的含量，获得盐类物质NaCl的含量，获得不同期次石英中成矿流体的均一温度和均一压力范围；
步骤(4)：根据查阅文献获得研究区矿床成矿流体演化的温度压力范围作为绘制的石英等溶解度相图温度和压力上、下限；
步骤(5)：根据流体包裹体测温获得的研究区成矿流体CO2和NaCl含量范围自由组合成矿流体的成分，以多种不同流体组成作为计算阶段成矿流体成分的输入值；
(二)计算阶段：
步骤(6)：利用H2O-CO2-NaCl体系相边界计算模块，分别计算成矿流体成分为1…n时H2O-CO2-NaCl流体体系液相与气液混合相的边界上的点；
步骤(7)：在研究区矿床成矿流体演化的温度范围内，以任意温度步长重复步骤(6)计算，得到一系列的点(T,P)，此即为H2O-CO2-NaCl流体体系液相与气液混合相的相边界；
步骤(8)：利用H2O-CO2-NaCl体系石英等溶解度计算模块，分别计算成矿流体成分为1…n时H2O-CO2-NaCl流体体系石英等溶解度曲线上的点；
步骤(9)：在研究区矿床成矿流体演化的温度范围内，以任意温度步长重复步骤(8)计算，得到一系列的点(T,P)，此即为H2O-CO2-NaCl流体体系在指定石英溶解度时的等溶解度曲线；
步骤(10)：以任意石英溶解度步长重复步骤(8)及(9)计算，得到H2O-CO2-NaCl流体体系不同溶解度条件下的石英等溶解度曲线；
步骤(11)：利用石英等溶解度曲线域边界计算模块中温度与压力敏感区边界计算部分，分别计算成矿流体成分为1…n时石英等溶解度曲线温度与压力敏感区的边界；
步骤(12)：利用石英等溶解度曲线域边界计算模块中退行性区边界计算部分，分别计算成矿流体成分为1…n时石英等溶解度曲线退行性区的边界；
(三)分析阶段：
步骤(13)：利用H2O-CO2-NaCl体系相边界计算模块、H2O-CO2-NaCl体系石英等溶解度计算模块及石英等溶解度曲线域边界计算模块获得的相边界曲线、石英等溶解度曲线和石英等溶解度曲线域边界，构建矿床成矿流体成分为1…n情况下的石英等溶解度相图；
步骤(14)：将实验阶段获得的不同期次石英中成矿流体的均一温度和均一压力范围投置于多种流体成分下的石英等溶解度相图中；
步骤(15)：根据石英等溶解度相图，分析其水热流体成矿过程。
CN202010560360.2A 2020-06-18 2020-06-18 定量化石英H2O-CO2-NaCl体系水热流体成矿过程的系统及方法 Active CN111508564B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010560360.2A CN111508564B (zh) 2020-06-18 2020-06-18 定量化石英H2O-CO2-NaCl体系水热流体成矿过程的系统及方法
JP2020140828A JP2021196336A (ja) 2020-06-18 2020-08-24 石英Ｈ２Ｏ−ＣＯ２−ＮａＣｌ系熱水流体鉱化プロセスを定量化する理論モデル及び方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010560360.2A CN111508564B (zh) 2020-06-18 2020-06-18 定量化石英H2O-CO2-NaCl体系水热流体成矿过程的系统及方法

Publications (2)

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

Family

ID=71878825

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010560360.2A Active CN111508564B (zh) 2020-06-18 2020-06-18 定量化石英H2O-CO2-NaCl体系水热流体成矿过程的系统及方法

Country Status (2)

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

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113539387A (zh) * 2021-07-09 2021-10-22 西南石油大学 一种基于CPA状态方程预测NaCl水溶液中CO2溶解度的方法
CN116359263A (zh) * 2023-03-10 2023-06-30 昆明理工大学 一种基于石英阴极发光特征的热液石英脉世代的判别方法
CN121188574A (zh) * 2025-11-07 2025-12-23 中国地质科学院矿产资源研究所 一种基于流体参数和流体场建模的成矿预测方法、设备、介质及产品
CN121188574B (zh) * 2025-11-07 2026-05-01 中国地质科学院矿产资源研究所 一种基于流体参数和流体场建模的成矿预测方法、设备、介质及产品

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115406880B (zh) * 2022-10-19 2023-03-24 核工业北京地质研究院 一种判断热液型铀矿成矿流体的还原性的方法
CN116665801B (zh) * 2023-04-25 2026-02-06 东南大学 一种确定MgO-MgSO4-MgCl2-H2O四元体系水泥相图的方法
CN119982470B (zh) * 2024-12-27 2025-10-31 中国石油化工股份有限公司 一种超临界co2压缩机自动控制方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1045762A (zh) * 1988-12-09 1990-10-03 巴特尔纪念研究院 采用控制化学计量和颗粒度生产钙钛矿复合物的亚微型陶瓷粉末的方法
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 (zh) * 2013-11-13 2015-05-20 中国石油化工股份有限公司 二氧化碳混相驱条件下油井受效阶段划分及判别方法
CN108344855A (zh) * 2018-02-06 2018-07-31 长江大学 碳酸盐储层成岩相边界点的确定方法
CN109488261A (zh) * 2018-12-18 2019-03-19 青岛海洋地质研究所 一种海洋天然气水合物三相动态虚拟构建方法
CN109752443A (zh) * 2018-12-24 2019-05-14 核工业北京地质研究院 一种确定砂岩型铀矿床成矿原因的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1045762A (zh) * 1988-12-09 1990-10-03 巴特尔纪念研究院 采用控制化学计量和颗粒度生产钙钛矿复合物的亚微型陶瓷粉末的方法
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 (zh) * 2013-11-13 2015-05-20 中国石油化工股份有限公司 二氧化碳混相驱条件下油井受效阶段划分及判别方法
CN108344855A (zh) * 2018-02-06 2018-07-31 长江大学 碳酸盐储层成岩相边界点的确定方法
CN109488261A (zh) * 2018-12-18 2019-03-19 青岛海洋地质研究所 一种海洋天然气水合物三相动态虚拟构建方法
CN109752443A (zh) * 2018-12-24 2019-05-14 核工业北京地质研究院 一种确定砂岩型铀矿床成矿原因的方法

Non-Patent Citations (2)

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

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113539387A (zh) * 2021-07-09 2021-10-22 西南石油大学 一种基于CPA状态方程预测NaCl水溶液中CO2溶解度的方法
CN116359263A (zh) * 2023-03-10 2023-06-30 昆明理工大学 一种基于石英阴极发光特征的热液石英脉世代的判别方法
CN121188574A (zh) * 2025-11-07 2025-12-23 中国地质科学院矿产资源研究所 一种基于流体参数和流体场建模的成矿预测方法、设备、介质及产品
CN121188574B (zh) * 2025-11-07 2026-05-01 中国地质科学院矿产资源研究所 一种基于流体参数和流体场建模的成矿预测方法、设备、介质及产品

Also Published As

Publication number Publication date
JP2021196336A (ja) 2021-12-27
CN111508564B (zh) 2023-10-03

Similar Documents

Publication Publication Date Title
CN111508564A (zh) 定量化石英H2O-CO2-NaCl体系水热流体成矿过程的理论模型及方法
CN106324700B (zh) 沉积盆地红层中砂岩型铀矿成矿模式构建方法
von Eynatten et al. Modelling compositional change: the example of chemical weathering of granitoid rocks
Okewale Applicability of chemical indices to characterize weathering degrees in decomposed volcanic rocks
Liu et al. Chemo-sedimentary facies analysis of fine-grained sediment formations: An example from the Lucaogou Fm in the Jimusaer sag, Junggar Basin, NW China
Yang et al. Correlation between mineralization intensity and fluid–rock reaction in the Xinli gold deposit, Jiaodong Peninsula, China: constraints from petrographic and statistical approaches
Jiu et al. Effect of hydrothermal fluids on the ultra-deep Ordovician carbonate rocks in Tarim Basin, China
Zhu et al. Hydrogeochemical and isotopic analyses of deep geothermal fluids in the Wumishan formation in Xiong’an new area, China
Zeng et al. Analysis of the mineral compositions of swell-shrink clays from Guangxi province, China
CN113534285B (zh) 一种针对热液型钍矿床成矿模式构建的方法
Lewandowski et al. Application of Raman spectroscopy to the flotation process of fluorite
Li et al. In-situ monazite Nd and pyrite S isotopes as fingerprints for the source of ore-forming fluids in the Jiaodong gold province
Nie et al. Determination of the thickness and fracture toughness of shale interfacial transition zone near various micromineral aggregates using nano-scratch tests
de Oliveira et al. Understanding regional groundwater flow using water-rock interaction processes and anthropic signatures in a fractured volcanic aquifer
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
Bomeni et al. Application of geostatistical methods to estimate the mineral contents in the alluvial clay deposit, Monoum plain, West Cameroon
Hu et al. Quantitative characterization of vertical zonation of Mesozoic granite weathering reservoirs in the coastal area of eastern Fujian Province, China
Chen et al. Fluid inclusion and in situ sulfur isotope studies of the carbonaceous strata-hosted Eryiling gold deposit in the Beishan orogen, NW China
CN116482212A (zh) 一种利用断层型古油藏约束卡林型金矿床成矿时代的方法
WO2022099397A1 (pt) Método para caracterização da matriz do solo e sua pedogênese com indicador geofísico indireto
Liu et al. Local sulfur sources control the development of Mississippi Valley-type Pb-Zn deposits worldwide
CN120577887B (zh) 一种构造单元产出砂岩型铀矿潜力的确定方法
CN119028468B (zh) 变质-深熔过程中成矿元素富集的定量模拟方法
Chen et al. Fluid evolution and genesis of the Shipenggou gold deposit in Central Jilin Province, China: Constraints from C–H–O and in–situ sulphur isotopes of pyrite and fluid inclusions
CN119914279B (zh) 基于元素录井分析的碳酸盐岩判别方法及系统

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