地球化学在沉积学研究中发挥了重要作用。本研究以灵璧石国家地质公园新元古代张渠组灰岩为研究对象,在系统采样的基础上测试了其岩石地球化学组成,并进行了数理统计分析。结果表明,张渠组上下部灰岩中元素组成存在差异,主体表现出高的Ca含量,但其他元素含量相对较低,这种差异可能与其沉积环境变化有关。基于数理统计分析(相关性、聚类及因子分析),10种元素被划分为3类:代表陆相的Al-K-Ti-Fe、代表海相的Ca-P和代表其他的Mg-Mn-Zr-Ni。在此基础上,构建出了用于探讨陆源、海相特征及陆海相对贡献的Al/Fe、Ca/P和Al/P比值,并被应用于张渠组沉积环境分析,发现张渠组上下部的海相环境无明显差异,但陆源贡献上变化显著。 Geochemistry has played important roles for the research of sedimentary geology. In this study, the limestone from the Neoproterozoic Zhangqu Formation in the National Geological Park of the Lingbi Stone has been taken as a research object and the geochemistry of the rocks has been measured based on the systematic collection of samples, and then analyzed by statistical analyses. The results indicate that the geochemical compositions of the limestone rocks from the upper and lower formation are different with each other, and both of them are characterized by high Ca contents with low other elements, which is probably related to the change of the sedimentary environment. Based on the statistical analyses (including correlation, cluster and factor analyses), the 10 kinds of elements have been classified to be 3 categories: the Al-K-Ti-Fe (representing the terrigenous detrital), the Ca-P (representing the marine facies) and the Mg-Mn-Zr-Ni (other). And three ratios (Al/Fe, Ca/P and Al/P) have been constructed for the discussion of the terrigenous, marine characteristics and the contribution of terrigenous versus marine, which was then applied for the analyses of the sedimentary environment of the Zhangqu Formation, and then found that there was no significant change of the marine environment during the formation of the lower and upper Zhangqu Formation, but the terrigenous contribution changed significantly.
王梁冰,孙林华,曹慧丽,何玉鹏,王雨航
宿州学院资源与土木工程学院,安徽 宿州
收稿日期:2018年3月21日;录用日期:2018年4月9日;发布日期:2018年4月16日
地球化学在沉积学研究中发挥了重要作用。本研究以灵璧石国家地质公园新元古代张渠组灰岩为研究对象,在系统采样的基础上测试了其岩石地球化学组成,并进行了数理统计分析。结果表明,张渠组上下部灰岩中元素组成存在差异,主体表现出高的Ca含量,但其他元素含量相对较低,这种差异可能与其沉积环境变化有关。基于数理统计分析(相关性、聚类及因子分析),10种元素被划分为3类:代表陆相的Al-K-Ti-Fe、代表海相的Ca-P和代表其他的Mg-Mn-Zr-Ni。在此基础上,构建出了用于探讨陆源、海相特征及陆海相对贡献的Al/Fe、Ca/P和Al/P比值,并被应用于张渠组沉积环境分析,发现张渠组上下部的海相环境无明显差异,但陆源贡献上变化显著。
关键词 :地球化学,灰岩,沉积环境,数理统计分析,新元古代
Copyright © 2018 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
岩石地球化学是研究岩石的化学组成,包含其成分的来源、含量、分布、种类及化学变化的地球科学,是近代岩石学和地球化学的交叉学科。通过研究各类岩石中的主、微量元素和同位素组成,进而分析岩石的源区、成因、演化和构造环境等基础理论问题。在沉积学领域,地球化学在古地理环境分析、海陆相地层划分、物源成分分析以及古气候条件研究等多个方面发挥了重要作用 [
皖北地区新元古代地层因为其独特性(产出灵璧石等多种独特岩石),而吸引了很多学者的关注,一系列研究得以展开。如郭峰等 [
在这套地层中,张渠组因为其丰富的灵璧石种类产出,尤其是磬石和臼齿构造的存在且明显不同于其他层位而备受关注,但目前缺乏针对性的研究,尤其是地球化学研究,这无疑大大制约了对张渠组沉积环境乃至多种灵璧石,尤其是磬石和臼齿构造成因的理解。
本文以灵璧石国家地质公园核心园区磬云山(位于灵璧县渔沟镇东部)的张渠组地层为研究对象,在野外调查的基础上,对其上下部进行了系统采样,并利用X荧光光谱仪进行了元素地球化学组成测试,进而对数据进行了数理统计分析,目的是查明元素之间的相互关系,进而构建具有不同意义的元素地球化学比值,从而分析其沉积环境变化。之所以进行数理统计分析,主要是因为几个方面的原因:一是现有的可用于碳酸盐岩的地球化学元素比值或者图解极其有限,二是从元素成因的角度来看,一个元素可能存在多种成因(来源或者存在形式),以往构建的元素比值可能并不具有普适性。以沉积岩中的Ca元素为例,其可以受制于碳酸盐矿物(如方解石,可能与海相有关),也可以受制于硅酸盐矿物(如斜长石,可能与陆相碎屑有关),如果在分析沉积环境过程中简单地将其与海相关联,可能会导致分析结果的错误。
研究区域位于华北板块南缘的淮北凹陷带 [
研究样品采自渔沟镇灵璧石国家地质公园内张渠组下部和上部,各采集15个样品,共计30个样品。张渠组下部为薄层灰黑色微晶灰岩,上部为灰色-浅灰色微晶灰岩。全岩样品首先去除岩石表面的粉尘,然后采用切割机切割岩石样品并选取岩石样品平整的新鲜面清洗烘干后进行测试。样品测试分析采用X荧光光谱仪(天瑞EDX5500H)完成。
图1. 研究区地质图
样品测试结果见表1。从表1可以看出,张渠组下部灰岩中Ca含量相对较高,介于54.1%~56.5%之
样品 | Al | Ca | Fe | K | Mg | Mn | Ni | P | Ti | Zr |
---|---|---|---|---|---|---|---|---|---|---|
H1 | 0.046 | 54.1 | 0.142 | 0.412 | 0.470 | 0.005 | 0.027 | 0.334 | 0.000 | 0.053 |
H2 | 0.110 | 55.1 | 0.197 | 0.534 | 0.459 | 0.012 | 0.031 | 0.349 | 0.012 | 0.075 |
H3 | 0.135 | 55.6 | 0.582 | 0.499 | 0.520 | 0.006 | 0.030 | 0.351 | 0.008 | 0.007 |
H4 | 0.074 | 55.1 | 0.877 | 0.431 | 0.476 | 0.005 | 0.029 | 0.349 | 0.011 | 0.092 |
H5 | 0.000 | 56.1 | 0.128 | 0.423 | 0.504 | 0.013 | 0.036 | 0.356 | 0.000 | 0.087 |
H6 | 0.232 | 55.3 | 0.358 | 0.542 | 0.514 | 0.009 | 0.030 | 0.352 | 0.023 | 0.090 |
H7 | 0.116 | 55.7 | 0.384 | 0.483 | 0.500 | 0.008 | 0.028 | 0.350 | 0.008 | 0.028 |
H8 | 0.009 | 56.5 | 0.093 | 0.353 | 0.490 | 0.010 | 0.029 | 0.356 | 0.002 | 0.061 |
H9 | 0.000 | 55.6 | 0.135 | 0.430 | 0.526 | 0.013 | 0.033 | 0.361 | 0.000 | 0.118 |
H10 | 0.151 | 54.9 | 0.383 | 0.462 | 0.560 | 0.009 | 0.029 | 0.348 | 0.005 | 0.104 |
H11 | 0.067 | 55.4 | 0.187 | 0.504 | 0.505 | 0.009 | 0.029 | 0.359 | 0.000 | 0.043 |
H12 | 0.000 | 55.9 | 0.119 | 0.381 | 0.550 | 0.012 | 0.032 | 0.356 | 0.006 | 0.111 |
H13 | 0.177 | 54.8 | 0.489 | 0.454 | 0.503 | 0.012 | 0.026 | 0.343 | 0.019 | 0.088 |
H14 | 0.290 | 54.4 | 0.467 | 0.600 | 0.547 | 0.014 | 0.033 | 0.339 | 0.014 | 0.098 |
H15 | 0.000 | 56.4 | 0.115 | 0.534 | 0.517 | 0.013 | 0.030 | 0.356 | 0.000 | 0.084 |
平均值 | 0.094 | 55.4 | 0.310 | 0.469 | 0.509 | 0.010 | 0.030 | 0.351 | 0.007 | 0.076 |
BH1 | 1.32 | 48.8 | 2.29 | 1.36 | 0.604 | 0.037 | 0.029 | 0.309 | 0.057 | 0.045 |
BH2 | 0.356 | 53.9 | 0.686 | 0.633 | 0.538 | 0.025 | 0.028 | 0.338 | 0.024 | 0.074 |
BH3 | 0.090 | 55.6 | 0.326 | 0.484 | 0.522 | 0.009 | 0.026 | 0.351 | 0.018 | 0.080 |
BH4 | 1.04 | 49.6 | 1.34 | 0.773 | 0.671 | 0.049 | 0.029 | 0.308 | 0.050 | 0.106 |
BH5 | 1.17 | 48.7 | 1.37 | 0.883 | 0.659 | 0.049 | 0.023 | 0.308 | 0.056 | 0.101 |
BH6 | 0.282 | 54.5 | 0.449 | 0.656 | 0.535 | 0.011 | 0.033 | 0.350 | 0.016 | 0.084 |
BH7 | 1.18 | 49.1 | 1.41 | 0.834 | 0.681 | 0.048 | 0.024 | 0.308 | 0.056 | 0.106 |
BH8 | 1.42 | 47.3 | 1.51 | 1.06 | 0.719 | 0.047 | 0.023 | 0.293 | 0.072 | 0.096 |
BH9 | 1.19 | 49.0 | 1.36 | 1.19 | 0.687 | 0.046 | 0.029 | 0.311 | 0.057 | 0.104 |
BH10 | 0.331 | 54.5 | 0.369 | 0.693 | 0.560 | 0.014 | 0.025 | 0.349 | 0.024 | 0.083 |
BH11 | 0.364 | 54.3 | 0.299 | 0.629 | 0.546 | 0.012 | 0.023 | 0.340 | 0.028 | 0.077 |
BH12 | 0.847 | 51.8 | 0.628 | 0.685 | 0.573 | 0.010 | 0.028 | 0.328 | 0.036 | 0.065 |
BH13 | 0.296 | 54.9 | 0.395 | 0.586 | 0.525 | 0.012 | 0.024 | 0.357 | 0.028 | 0.081 |
BH14 | 0.317 | 54.5 | 0.466 | 0.600 | 0.556 | 0.010 | 0.029 | 0.340 | 0.008 | 0.084 |
BH15 | 0.160 | 55.1 | 0.334 | 0.483 | 0.477 | 0.010 | 0.028 | 0.345 | 0.023 | 0.078 |
平均值 | 0.691 | 52.1 | 0.881 | 0.770 | 0.590 | 0.026 | 0.027 | 0.329 | 0.037 | 0.084 |
表1. 微量元素组成(%),H和BH分别为下部和上部样品
间,平均为55.4%,其他元素的含量相对Ca而言较低,Al、Fe、K、Mg、Mn、Ni、P、Ti和Zr含量的平均值分别为0.094%、0.310%、0.469%、0.509%、0.010%、0.030%、0.351%、0.007%和0.076%。所有元素含量从高到低依次为Ca > Mg > K > P > Fe > Al > Zr > Ni > Mn > Ti。
相比较而言,张渠组上部薄层灰岩具有相对下部偏低的Ca含量,其平均值为52.1%,而Al、Fe、K、Mg、Mn、Ti和Zr则相对上部灰岩偏高,其含量分别为0.691%、0.881%、0.770%、0.590%、0.026%、0.037%和0.084%,Ni和P含量相对上部灰岩偏低,其平均值分别为0.027%和0.329%。
相关性分析是指对两个或多个具备相关性的变量元素进行分析,从而衡量两个变量因素的相关密切程度。在地球化学研究中,通常将相关性好的元素考虑成具有共同的来源 [
本文研究中对上下部岩石进行分开研究,其相关性分析结果见表2和表3。因为各有15个样品,对应显著性因子α = 0.05的相关系数临界值为:rα = 0.514。从表2可以看出,下部灰岩Al-K-Ti-Fe、Ca-P之间具有显著的正相关,而Al和Ca之间则表现出显著的负相关。对于上部灰岩而言(表3),所有的元素大概分为两组,且两组之间表现为明显的负相关关系,其中第一组Al-Mg-K-Ti-Mn-Fe-Zr之间具有明显的正相关,第二组Ca-P之间具有明显的正相关。综合来看,无论上部还是下部灰岩,Al-K-Ti-Fe之间都具有显著的正相关,而Ca-P之间也具有明显的正相关,且这种关系表明这些元素之间可能存在成因上的联系 [
聚类分析指基于某些变量将一组研究对象分为相对同质的群组的统计分析技术。在地球科学研究中,聚类分析应用广泛,如在污染物来源解析中,利用R型聚类分析将污染物划分为几个不同的类,每个类
Al | Ca | Fe | K | Mg | Mn | Ni | P | Ti | Zr | |
---|---|---|---|---|---|---|---|---|---|---|
Al | 1.000 | |||||||||
Ca | −0.604 | 1.000 | ||||||||
Fe | 0.544 | −0.388 | 1.000 | |||||||
K | 0.693 | −0.304 | 0.241 | 1.000 | ||||||
Mg | 0.230 | 0.079 | −0.024 | 0.148 | 1.000 | |||||
Mn | −0.023 | 0.278 | −0.474 | 0.216 | 0.377 | 1.000 | ||||
Ni | −0.180 | 0.321 | −0.286 | 0.090 | 0.317 | 0.581 | 1.000 | |||
P | −0.552 | 0.832 | −0.353 | −0.213 | 0.180 | 0.295 | 0.411 | 1.000 | ||
Ti | 0.800 | −0.382 | 0.563 | 0.445 | 0.002 | −0.003 | −0.220 | −0.366 | 1.000 | |
Zr | −0.028 | −0.073 | −0.113 | −0.091 | 0.383 | 0.570 | 0.396 | 0.100 | 0.126 | 1.000 |
表2. 相关性分析结果(下部样品)
Al | Ca | Fe | K | Mg | Mn | Ni | P | Ti | Zr | |
---|---|---|---|---|---|---|---|---|---|---|
Al | 1.000 | |||||||||
Ca | −0.994 | 1.000 | ||||||||
Fe | 0.917 | −0.911 | 1.000 | |||||||
K | 0.865 | −0.841 | 0.912 | 1.000 | ||||||
Mg | 0.915 | −0.934 | 0.770 | 0.725 | 1.000 | |||||
Mn | 0.887 | −0.923 | 0.853 | 0.729 | 0.918 | 1.000 | ||||
Ni | −0.171 | 0.178 | −0.016 | 0.023 | −0.245 | −0.204 | 1.000 | |||
P | −0.964 | 0.980 | −0.882 | −0.781 | −0.921 | −0.913 | 0.179 | 1.000 | ||
Ti | 0.954 | −0.960 | 0.864 | 0.811 | 0.887 | 0.896 | −0.335 | −0.927 | 1.000 | |
Zr | 0.272 | −0.345 | 0.074 | −0.004 | 0.589 | 0.545 | −0.239 | −0.357 | 0.327 | 1.000 |
表3. 相关性分析结果(上部样品)
代表不同的来源 [
在本次研究中,对所有元素进行了R型聚类分析(聚类方法为Ward,距离为Pearson),结果如图2所示。从图中可以看出,无论是下部灰岩还是上部灰岩,所有的元素都分为两个大类,但上下部灰岩之间存在一定的区别。其中下部灰岩Fe-Ti-Al-K为一类,而Mg-Zr-Mn-Ni-P-Ca为一类。上部灰岩Ca-P-Ni为一类,而K-Fe-Ti-Al-Mg-Mn-Zr为一类。这一结果与相关性分析结果基本上是一致的(除了相关性分析结果中缺少Ni)。
因子分析是基于降维的思想,把多个具有复杂关系的变量合并为几个综合因子的统计分析方法 [
本次研究对上下部灰岩分别进行了因子分析,结果如表4所示。从表中可以看出,在选择特征值大于1且进行正交旋转后,下部灰岩共提取出了3个因子,其总的方差解释率仅为56.1%,而上部灰岩仅提取出了2个因子,其总的方差解释率仅为64.0%。对于下部灰岩而言,因子1的方差解释率为19.1%,在Al、K和Ti上具有较高的正载荷,而在Fe上具有中等的正载荷;因子2的方差解释率为24.5%且在Mn和Zr上具有较高的正载荷,而在Mg和Ni上具有中等的正载荷;因子3的方差解释率为12.5%且在Ca和P上具有较高的正载荷。对于上部灰岩而言,因子1的方差解释率为46.0%,在Al、Fe、K、Mg、Mn和Ti上具有较高的正载荷,而在Ca、P上具有较高的负载荷;因子2的方差解释率为18.0%且在Zr上具有较高的正载荷,而在Ni上具有中等的负载荷。
从上述数理统计分析中可以发现两个方面的问题:1) 不同沉积环境下形成的灰岩在元素组合上存在明显的差异;2) 同一元素在不同沉积环境下形成的灰岩中其来源(或者影响因素)存在不同。因此,如果
图2. 聚类分析结果(左下部右上部样品)
F1 | F2 | F3 | F1* | F2* | |
---|---|---|---|---|---|
Al | 0.889 | 0.053 | −0.385 | 0.974 | 0.160 |
Ca | −0.303 | 0.076 | 0.877 | −0.969 | −0.220 |
Fe | 0.598 | −0.345 | −0.279 | 0.971 | −0.085 |
K | 0.828 | 0.151 | 0.050 | 0.931 | −0.181 |
Mg | 0.222 | 0.618 | 0.110 | 0.860 | 0.451 |
Mn | 0.013 | 0.856 | 0.222 | 0.879 | 0.388 |
Ni | −0.059 | 0.683 | 0.382 | −0.034 | −0.687 |
P | −0.251 | 0.210 | 0.850 | −0.941 | −0.244 |
Ti | 0.764 | 0.024 | −0.318 | 0.922 | 0.290 |
Zr | −0.153 | 0.823 | −0.294 | 0.170 | 0.848 |
特征值 | 1.91 | 2.45 | 1.25 | 4.60 | 1.80 |
方差解释率 | 19.1% | 24.5% | 12.5% | 46.0% | 18.0% |
表4. 因子分析结果(F1-F3下部样品,F1*和F2*上部样品)
利用元素(包括其比值等)来分析灰岩的沉积环境,固定地采用一些元素及其比值可能存在不适用的问题,反而自己根据研究对象的元素组合关系来选择合适的元素并构建新的元素比值可能更为合适。例如在本次研究对象中,下部灰岩中Mg、Zr、Mn和Ni可能主要与海相沉积环境有关,因为无论是相关性、聚类还是因子分析中都可以发现,这几个元素与Ca、P具有相对密切的关系,但对于上部灰岩而言,这几个元素则明显与陆相有着更为密切的关系(图2)。
从以上分析中可以发现,对于本研究中的10种元素而言,无论是在上部灰岩还是下部灰岩中,它们大体上可以分为3类:一是以Al-K-Ti-Fe为代表的陆相相关元素,二是以Ca-P为代表的海相相关元素,三是包括Mg-Mn-Zr-Ni在内的变化元素(其来源不稳定,可能受制于海相也可能受制于陆相)。因此,Al-K-Ti-Fe之间的比值应该可以用于分析陆相的特征,反之Ca-P之间的比值则可以用于分析海相的特征,而Al-K-Ti-Fe与Ca-P之间的比值则可能更多的反映了陆相和海相之间的贡献差异。基于这种思想,并结合前人的研究成果,我们构建了如下比值,包括Al/Fe、Ca/P和Al/P,并对其地质意义进行了分析讨论。分别如下:
1) Al/Fe:因为Al和Fe均与陆源碎屑有关,因此这一比值主要反映了陆源物质的特征。计算结果表明,下部灰岩中Al/Fe为0~0.648 (平均0.303),上部灰岩为0.276~1.35 (平均0.784)。这种情况说明,张渠组上下部灰岩中的陆相贡献部分在来源上可能发生了改变。
2) Ca/P:因为Ca和P均与海相有关,因此这一比值主要反映了海相物质的特征。计算结果表明,下部灰岩中Ca/P为154~162 (平均158);上部灰岩为154~161 (平均158),上下部是一致的。这种情况说明,张渠组上下部灰岩中的海相贡献部分并没有明显的不同,海水的环境(组成特征)相对是稳定的。
3) Al/P:该比值主要反映海相和陆相贡献的相对变化。下部灰岩中Al/P为0~0.855 (平均0.268),而上部灰岩的Al/P为0.256~4.85 (平均2.10),两者存在明显的差异,说明上部灰岩中陆相贡献比例明显增加。这与实际情况是一致的,下部薄层灰岩相对纯净,而张渠组上部中逐渐出现了大量的泥质(如泥灰岩)。
在以往的研究中,有诸多元素及其比值被应用于沉积环境分析,如Fe/Mn被用于分析水的深浅 [
通过对皖北地区新元古代张渠组上下部灰岩的岩石地球化学数理统计分析,取得了如下认识:
1) 张渠组上下部灰岩中Ca含量相对较高,其他元素含量相对较低且上下部灰岩在元素组成上存在差异,可能与沉积环境及物源变化有关;
2) 相关性、聚类及因子分析将元素分为3类,Al-K-Ti-Fe、Ca-P和Mg-Mn-Zr-Ni,分别代表陆相、海相和其他;
3) 构建了Al/Fe、Ca/P和Al/P,分别用于探讨陆源、海相特征及陆海相对贡献,并发现张渠组上下部的海相环境无明显差异,但陆源贡献上变化显著。
国家大学生创新创业项目(201710379033)。
王梁冰,孙林华,曹慧丽,何玉鹏,王雨航. 皖北新元古代张渠组灰岩地球化学数理统计分析 Statistical Analyses of Geochemistry of the Limestone from the Neoproterozoic Zhangqu Formation in Northern Anhui Province[J]. 地球科学前沿, 2018, 08(02): 266-274. https://doi.org/10.12677/AG.2018.82028