化探第13章++蒸发盐类沉积环境和历史的研究.docx

1、化探第13章+蒸发盐类沉积环境和历史的研究第十三章 蒸发盐类沉积环境和历史的研究海水中总是含有少量的Br。由于它在海水中的浓度低，不能形成Br的独立矿物，一般都是呈类质同相混入物形式分散到氯化物类矿物中去。布赖奇和赫尔曼（Braitsch and Hermann，1963）以实验方法测定了Br在海水与NaC1晶体之间的分配系数（T25C岩盐开始从海水中沉淀讨）：KBr BrNaClBr 海水0.15这意味着，一但岩盐开始沉淀，海水（卤水）Br含量的增长要比海水盐度增长快的多。同时还证明直到Br的浓度很高时，KBr值仍然恒定不变，但MgC12在溶液中的含量增高则影响KBr。应用Br的这种分配规律

2、，如果在蒸发岩沉积序列中发现了一个NaC1层，就可以根据NaC1中的Br含量来判断NaC1自其中沉淀的古海水盐度。蒸发盆地的历史往往是复杂的，图18表示出一个蒸发盆地发展的两种可能类型。图中下面的曲线说明一个封闭盆地进行简单蒸发的情况，而上面的曲线则是反映海水的补充刚好抵消由于蒸发而失去的水量的情况。而在后一种情况下，盆地中水的体积是恒定的，但其盐度则在逐渐增长。对于固定体积，“半开放” 的盆地来说，必然会形成具厚的盐层序列。德国斯特斯福特（Stassfurt）附近发现的盐类地层，称为泽施田系（ZeChst1einSeries）。盐类层序研究证明，该区的许多蒸发盆地过去都经过几个阶段，图19反

3、映了该岩系底部的情况，说明由于岩盐的形成而引起的卤水中盐度的最初变化，根据图中的曲线，考虑Br的分配规律，可以认为有很厚的岩盐曾是在海水不断补给盆地的情况下沉淀的。由于Br在NaC1中含量的增长来得慢些（对比图18与19），可以推断当时海水流入量必定超过由于蒸发而失去的水量。剖面的顶部盐类中，Br含量曲线反映出封闭盆地快速变化的情况，其结果是盐度迅速增高，而最终导致钾盐的沉淀。 图18 a岩盐中Br的含量与两种水体平衡情况盆地中沉积盐类厚度的关系（据Holserr,1966） b泽施田岩系盐类Br中的沉积盐类厚度的关系（据Holserr,1966）地质温度计：微迹元素服从亨利稀溶液定律，K基本

4、受T影响，而独立于组分的浓度，做地质温度计比较理想。KD与T的关系一般表达如下：1nKDAT1b通过实验确定了两相间某一微量元素的KD与T的关系式，则可作为地质温度计。例如，通过实验确定了Se在方铅矿和闪锌矿（两者平衡共生）之间的分配系数（KSeGLSph）与绝对温皮（T）的关系式为：1n KSeGLSph （1T0.00044） 0.00035根据共生方铅矿和闪锌矿中测定的Se含量，计量出KSeGLSph。（CSeGLCSeSph）的值，代入上式，即可得出这两种矿物形成的温度。这就是方铅矿闪锌矿温度计（适用范围：8721163K）。其他还有利用Rb在金云母和透长石之间分配的地质温度计：1n

5、KRb金透 0046.103T0.091利用Ni在橄榄石和辉石之间分配的地质温度计：1n KNi橄辉 16.8RT7.65还有利用Na在钾长石和斜长石之间分配的长石温度计等等。其他应用：分配定律的应用范围正在逐年扩大，目前，它还常被引用来解决下列间题：（1）检验变质作用是否达到平衡；（2）解释和确定变质区域中的等温线和等压线；（3）研究岩浆形成的部分熔融过程；（4）花岗岩的成因；（5）成矿溶液的成分；（6）岩浆与气液间金属元素的分配；（7）挥发元素的地球化学行为；（8）古海水含盐度；（9）利用固气相间元素分配规律研究太阳星云凝聚过程；（10）作为寻找盐类矿床的指示等。从1967年开始，分配定律

6、已被引入化探领域，例如前面介绍的巴尔苏科和奥洛索关所设计的利用F含量指示盲矿规模和埋深的方法中，就引用了F在多相平衡中的分配系数。前苏联学者又开创了应用分配系数随温度变化的规律（温度随深度增高，指示元素在脉和晕间的分配系数也随深度而变化），探索指示矿床侵蚀截面深度的新途径。第十四章 有关成矿成晕的某些基础理论、关于原生晕形成的补充材料在有关地球化学找矿的某些基础理论中，己指出在热液矿床原生晕的诸特征中最具有理论与实际意义的应为原生晕的分带结构，并且重点讨论了阐明原生晕分带的某些理论探索情况，例如：（1）元素搬运形式（络合物）的稳定性；（2）元素沉淀的条件；（3）类质同相；（4）吸咐作用；（5）

7、扩散；（6）渗滤和过滤效应：（7）其他因素。从以上7各方面等角度探讨了原生晕的分带间题。现在主要补充有关热力学研究的某些进展情况及原生晕研究必须注意的其他因素。1、关于热液中络合物稳定性的热力学探讨巴恩斯（Barnes，H，L，1962）曾运用热力学计算方法近似地标定成矿金属元素络合物（络离子，多具共价健）的相对稳定性。一个简单的络合物生成反应可写作：M2 （“.） nAz- （“） MAn2nx并证明可以用：DG0C（未知常数）来度量具有相同价态、不同元素的同类络合物的稳定性，其中G0是络合物形成反应的标准自由能。对于各种金属形成硫化络合物，巴恩斯赋于未知常数C以200千卡的值，然后计算这些

8、络合物的相对稳定性。例如： Hg2（aq）2S2HgS22DG0CG0 =G0 HgS22G0 Hg22G0 S211.639382196271.7千卡 DG0C71.7200271.7千卡克分子（本例与Barnes 的计算结果不一致，可能是由于采用的热力学参数值不同或络离子的形式不同）。通过这样的计算，巴恩斯按同类络合物稳定性减小顺序将不同金属做如下排列：元素 Hg Cd Pb Cu Zn Sn Ni Fe Co MnD值 227 158 154 134 132 126 83 82 81 78这一结果与热液矿床原生晕中观察到的轴向分带序列十分一致。研究证明，当电解质浓度增大时，络合物将变得越

9、来越重要；随着温度的上升，络合物的作用将更加明显。然而，由于缺乏高温下的数据，这方面的工作往往是从低温关系外推的。尽管如此，现有的工作仍足以证明，以热力学估计的络合物稳定性来阐明原生晕分带是可取的。2、元素沉淀条件的热力学研究及原生晕分带奥夫饮科夫和格里戈良在研究大量热液矿床原生晕的基础上，综合得出统一的原生晕分带序列。这有力的说明，不同类型、不同成因、不同地质条件下的成矿溶液应具有相似的阳离子阴离子成分和性质，反映出成矿过程中体系的物理化学条件在时间和空间中变化的范围是极窄的和单方面的。巴夫洛夫（A.JI.B，1976）通过热力学分析的方法证明，原始热液在其刚分出时应具有明显的碱性和还原特征

10、，这决定着成矿热液阳离子阴离子成分的特殊性。因此矿石沉淀的基本原因的最大可能是：体系物理化学条件，首先是PH和Eh的自发演化，即在氧化电位不断增长的基础上（由于热液从地壳深部向浅部富氧环境上升），随着温度和压力的降低，成矿溶液的酸度逐渐向增高方向自发演化。正是成矿溶液的酸碱度和氧化还原电位等条件，决定着矿物共生组合的成分及它们在时间和空间中的沉淀顺序。在上述思想的指导下，巴夫洛夫尝试把硫化物矿床原生晕轴向分带序列同金属元素硫化物和硫氧络离子由于原始碱性还原热液的酸性（和氧化）发生分解，而引起的硫化物和自然金属沉淀的顺序进行对比。为了进行这种对比，巴夫洛夫提出以下两点假定：（1）硫化物是同金属元

11、素一道在一种统一的成矿溶液中被搬运到矿石沉淀带中，并且溶液中硫化物硫（S2、HS、H2S）的活度足够高（101克分子升）；（2）硫化络离子（MSnx）和硫氢络离子（MHSnx）活度的下限（低于它时，溶液就被认为不能成矿或不育的）和决定着矿物稳定场与络合物溶液之间平衡的这两种络离子的活度都假定为相等，即都为106克分子升。为了确定反映PH值与溶液中硫化物和溶解络合物活度之间的函数关系，采取了如下方法。例如，对锡的硫化络离子（SnS32）在溶液酸化过程中的分解，就可以采取下列途径找出所要求的函数关系或方程。在含有络离子SnS32的碱性溶液酸化时（SnS32、G0=24.21千卡克分子），SnS2（

12、G0=47.55千卡克分子）的沉淀按下列反映式进行：G0=SnS32H SnS2HS在这种情况下，PH值和溶解组分活度之间的函数关系式可以推导如下：依据，G0反应=RTlnK=1.3641g （1） 得出：PH （2）由于G0反应47.55十3.01十24.210=20.33千卡克分子代入（2），得：PH=（20.331.364） （3）此即所求方程式。考虑到上述提出的假定：101克分子升， =106克分子升，在形成原生晕过程中难溶化合物SnS2的沉淀将在溶液PH值降低到10.0就已发生（PH15.01g1011g10615.51610.0）。利用己有的热力学常数，计算得出原始碱性热液酸化过程

13、中，主要通过硫化和硫氢络合物的分解顺序沉淀出硫化物、自然元素矿物或氧化物时的PH值（表1）。 矿物自热液中沉淀的PH值 表1序列迁移形式沉淀形式矿物沉淀的PH值备注BaAsSbHg1Cu1AgSn1Hg2AuCu2CdPbZnCuBiUMoCoNiSn2BaS、Ba2AsS2、AsS33、HAsS2SbS2-、HSb2S4HgS22-、Hg(HS)3-CuS、Cu(HS)3-AgS、AgHSSnS32-、SnS44-HgS22-、Hg(HS)3-AuS-CuS22-、Cu(HS)3-Cd(HS)3-PbS22-、Pb(HS)3-ZnS22-、Zn(HS)3-CuS22-、Cu(HS)3-BiS

14、33-、BiS2-UO(CO3)34-MoS42-Co(HS)2、Co(HS)Ni(HS)2、Co(HS)3SnS32-、SnO44-硫酸盐硫化物硫化物自然元素硫化物硫化物硫化物硫化物自然元素硫化物（不含铁）硫化物硫化物硫化物硫化物（含铁）硫化物氧化物硫化物硫化物硫化物硫化物、氧化物低PH值、高Eh值4.04.54.35.25.25.86.55.05.57.37.78.28.88.59.19.2（到5.0）9.510.07.79.88.59.910.510.0CuS型CuFeS2型分析表中资料，可以看出存在物理化学梯度的情况下，PH值相当可靠地控制着元素的带状分配。显然，自然条件下情况将更为复

15、杂，因为围岩结构构造特征的不同、存在温度和压力梯度等因素都在发生影响。但上述PH值与溶解络合物活度之间存在的函数关系，可以使我们有根据地认为，原生晕成分和结构方面存在的规律首先取决于热液PH值和Eh自变演化过程中硫化物、自然元素、氧化物和硫酸盐矿物的沉淀分带。正是由于存在这种情况，硫化物矿床原生晕的分带规律实际上才很少取决于围岩的成分。3、奈恩斯特（Nest）分配定律在研究原生晕中的应用梁布契科关（1965）从理论上证明，奈恩斯特分配定律也适用于元素在热液脉（与围岩有清楚的界限），与脉旁蚀变晕之间的分配，亦即适用于多项平衡的元素分配。这时分配系数KCPCoP。其中，CP脉中给定元素的平均含量或

16、线金属量：CoP晕中该元素的平均含量或线金属量。根据热力学原理，分配系数应为温度、其次是压力的函数。在温度和压力固定的条件下K值为固定常数，在温、压发生变化的条件下，K值也应发生变化。按照卢戈夫等人（C.JIy O B，1972）的资料，热液矿床矿石沉淀过程是在梯度达到40C100m的温度场中进行的，因此可以设想，随着矿脉倾向方向上成矿条件的总变化（主要是温度，其次是压力）元素在脉和晕之间的分配系数值（K值）也应发生规律的变化，从而为研究原生晕在空间的分带特征开创了新的途径。严格导出的KACoPCpBMOP，M和Mo分别为脉和晕的厚度。1967年巴尔苏科夫和奥洛夫（B.，1968）曾应用F在脉

17、和晕之间的分配系数，或其倒数（KP1K）预测苏联远东脉状锡矿化带深部锡的线金属量和矿体深部的上下界。以后的工作证明，根据Q=68.559.3KP关系式，预测下伏锡的线金属量（Q），方法是有效的。为了进一步探索应用奈恩斯特(Nest）分配定律研究原生晕和预测深部盲矿的可能，克拉夫钦柯、鲍洛琴柯夫和马克叶夫（.C.Kp，A.，H.Maee）又在前苏联两个锡矿床莫帕和少先队矿床开展了试验研究。莫帕矿床产于中生代酸性和和中酸性火山岩中，是一系列互相接近的含少量硫化物及锡石的石英脉。少先队矿床产于白垩纪粉砂岩中，为具含石英硫化物成分的脉状含锡石和黝锡矿的矿床。为了研究地球化学分带进行了800个刻槽和岩心

18、的取样，并在原生晕和矿体中用检块法采集了3000个岩石样品。所有样品都做了近似定量光谱分析。在莫帕矿床中，采用元素在脉和晕中平均含量的比值计算分配系数（KCpCop）。在几条勘查线剖面上，按34个不同标高的截面计算了分配系数值。兹举一条剖面为例来说明分配系数值的变化：中段（m） Pb Sn W Ag Li Mo Bi900 1.75 23.8 250 7.0 137 2.17 0.5800 0.86 11.5 1.13 108 1.01 2.70 0.61760 082 3.06 1.05 1.12 0.66 2.29 0.62680 074 2.06 0.9 0.84 0.28 54.13

19、3.43由上列数据可以看出，自矿体下部向上，分配系数值的总体变化是有规律的，但对不同元素来说变化的方向则是不同的：Pb、Sn、W、Li的分配系数值是向上增长，而Mo和Bi的分配系数值则为向上逐渐降低。可以认为，分带标志的这种变化差异应是元素地球化学个性（行为）在相同或近似相同条件下的反映。所以观察到的情况对于其他剖面来说也是稳定的，亦即在矿床不同部分分布特征基本不发生重大变化。对于少先矿床的矿体和矿体，采用元素在脉和晕中的线金属含量（QLP和QLOP）的比值为分配系擞（KQLPQLOP），做为分带标志。由下列数据可见，从矿体下部中段向上分配系数值同样是在方向性上具有变化的：矿体I中段（m） A

20、s Pb Cu Sn Bi Ni Ag900 9.0 2045 20.0 164.0 2.0 0.27 31.8 750 3.8 2.6 19.5 3.07 130 0.30 1.0510 606 2.7 0.95 0.80 0.09 0.80 矿体 中段（m） As Pb Cu Sn Bi Ni Ag920 0.23 0.39 0.18 1.1 1.4 9.18 1.1740 74 0.27 5.1 1.3 5.0 0.59 3.5460 16.6 12.6 10.6 3.2 12.6 0.72 10.3但是，同时应该看到在这里分配系数值变化的规律也不是简单划一的。因为，在两过矿体中对所有元

21、素（无例外）的分配系数值都是在相反方向中变化的。此外，在同一个矿体范围内的不同剖面中也发现存在同种元素分配系数值的变化是有规律的，但却又是相当复杂的，如果简单化地应用它作为地球化学分带标志，预测深部盲矿则可能导致错误结论。为了进一步查明元素分配系数值在矿体周围空间中的变化规律，克拉夫钦柯等人曾将Sn和Ag在脉和晕间的分配系数值以投影法表示于矿床和矿体的垂直剖面上（图1）。 图1 锡（a）和银（b）在矿脉和晕之间分配系数等值线在垂直平面上 的投影（少先队矿床之矿体）由图可以看出，两种元素分配系数等值线的轮廓总体上是相似的。这种共生应该说不是偶然的，分配系数值线在矿体中心部分圈出Sn和Ag优先在晕

22、中集中的地段（KQLPQLOP值低，QLOPQLP）由这个地段向上和向外，沉淀在矿脉（矿体）中的元素总量所占的比例逐渐增大。在这里围岩的性质是均一的，这就排除了分带规律受围岩成分和结构构造特征影响的可能性。由前述原理出发，分配系数值（分带标志）最大的变化方向应该就是成矿条件，首先是温度和压力变化最大的方向。因此，可以推测任何分带标志在脉状矿体整个延伸部分都将是不稳定的，但这并不排除应用这种分带规律来定量地预测深部盲矿的可能。为此就需制定选择标准剖面的程序和定量预测深部盲矿的方法。克拉柯钦夫等认为，制定程序应该同时研究分带在矿体平面中的一般情况，而选择剖面应该反映垂直分带的最典型情况。利用这些原

23、则，他们在莫帕和少先队矿床各选择一条勘查剖面作为标准剖面。所研究截面标高的确定是借助Q相关矩阵分析型因子分析的矩心法来实现。这种类型的分析预先研究截面对之间在一系列标志的测量值方面的相关关系，预测任务被归结为：确定矿化的可能剥蚀深度和该截面以下（与标准截面相比）的远景资源量。为了说明预测的正确性，克拉柯钦夫等人列出了他们通过计算得出的预测数据与勘探过程实际获得的数据（表2），对比两种数据后认为，所提出的方法可以解决定量评价成矿参数的任务，其精度是可以被接受的。因为，真值和计算值在顺序方面基本上是一致的。而且计算值落在标准剖面相同的区间内，根据勘探资料已知的真值亦如此。莫帕和少先队矿床中计算数据

24、与勘探获得的数据对比 表2矿床中段标高（m）计算出的中段标高（m）截面以下的真实资源量（m2%）截面以下计算的资源量（m2%）莫帕矿床800840770700820308251579056955290100038.502931529840980少先队矿床720640950800900750710257001090556901910307251042.311.41471624910326822421630201255504510综上所述，可以有根据地说，在脉状锡矿床中存在着由奈恩斯特分配定律所决定的地球化学分带。带状分配系数值的空间图象可能是复杂的，这种复杂的分带大概反映着在热液沉淀矿质的地段中

25、物理化学条件的变化。虽然利用元素分配系数的分带规律来预测矿产还存在许多困难和间题（因为目前这方面的研究还很肤浅），但初步研究所取得的成就己经表明这是一个很值得继续深入探讨的间题。4、成矿和成晕能量前苏联探矿方法和探矿技术研究所萨费隆诺夫等人运用热力学的基础理论，研究成矿成晕能量，取得一定进展。1978年出版了成矿能量与找矿，论述了成矿能量的概念计算方法及其在化探中的应用。各种金属成矿元素从原始分散状态到局部富集成矿成晕，需要消耗相应的能量，犹如选矿场把低品位的矿石，按照一定流程，通过各种手段，将其变为达到冶炼品级的情况，要消耗电能、机械能、化学能等各种形式的能量。尽管目前要把成矿成晕过程中能量

26、转换的全部形式表达出来是困难的，但可以设想成矿和成晕能量与矿石和原生晕的克拉克浓度之间有一定的关系。萨费隆诺夫等人从这一设想出发，以热力学为基础，推导出计算矿石和原生晕形成能量的公式：En式中：En为形成单位体积矿石或原生晕所消耗的能量（相对单位）：K为组成矿石或原生晕的任一元素的克拉克值。在计算成矿成晕能量时，必须首先求得与矿床有成因关系的某一岩浆体（即母岩）中每种元素的单位体积的克原子量，即克原子密度，其步骤如下：（1）根据地质及地球化学特征，确定形成矿床的母岩；（2）对母岩做硅酸盐全分析及有关元素的简相分析，获得母岩的平均化学成分：（3）将氧化物重量百分比换算为元素的重量百分比；（4）测

27、定母岩的密度（单位体积的重量）； （5）计算各元素在母岩中的密度，即母岩密度乘以各元素的重量百分比：（6）算出各元素的克原子密度，即用各元素的原子量除各元素的密度； （7）用同样的方法，计算出矿石或原生晕中各元素的克原子密度；（8）然后用矿石或原生晕中各元素的克原子密度与母岩中对应元素的克原子密度之比值，求出各元素的克拉克浓度Ki，代入公式： En即可算出单位体积矿石或原生晕的形成能量（相对单位）。现以某锡石石英脉矿床的数据为例，说明矿石形成能量的计算程序。成矿母岩和矿石的化学成分和计算步骤都列于表3中。在地球化学找矿中，可用原生晕形成的相对能量为指标，圈定异常，往往比用元素含量圈定异常更有效

28、、更直观。原生晕的形成能量用上述相同的方法进行计算，其中元素克拉克浓度值可与矿区内有成因关系的岩浆岩中的元素克拉克浓度值相比较，也可与地壳中相应岩浆岩相比较，还可与工作区域的元素区域背景值相比较。例如，某地区金属矿床，在矿区2.4km2的范围内，用5020m间距取样（岩石样）。分析原生晕的成分，求出每个取样点的综合能量，并按E10、E1050、50100、E100等间隔圈出等能量线（图2a）。在地球化学异常图上可以明显地圈出两个东西向延伸的形成能量较高的异常带。北边的异常带中包括了己知矿体，而南边的异常带可作为寻找盲矿体的远景区。对一些成矿作用复杂，多元素综合性矿化的地区，在原生晕找矿工作中，

29、用元素的含量圈定异常区，往住比较困难。但如用单位体积原生晕的形成能量作指标整理化探资料，可获得较为满意的结果（图2b）。在一个覆盖较厚的变质岩地区，用浅钻和地表工程取样，以成晕能量为指标，计算出各剖面上每一个钻孔和地表工程的相对成晕能量。从综合能量分布特点看，在变质岩和闪长岩中有两出成矿能量较高的地段，可能找到具有远景的金矿化。除了在原生晕研究中的应用外，萨费隆诺夫等人还提出成矿能量用于计算不同品级矿石的相对资源量，阐明矿石品级的结构变化，以预测各类品级矿石的相对远景资源量的方法和实例，说明这种研究对找矿勘探工作具有重要意义。 图2 a多金属矿区地球化学异常图（原生晕能量100；矿体）； b某金矿区一个剖面的综合能量分布图 （据.，1978）某锡石石英脉矿床矿石形成能量的计算 表3元素（g/m3）OSiAl Fe Ca Na K Mg Ti H CSn其他金属2680花岗岩（母岩）48.732.37.72.71.582.773.340.560.230.090.03310-40.492681锡石石英脉（矿石）51.4338.754.451.611.910.190.120.330.06