地质—地球化学建模在湖南沅麻盆地的应用
杨保国
宁夏地球物理地球化学勘查院,宁夏 银川 750004

作者简介: 杨保国(1985-),男,硕士研究生学历,毕业于中国地质大学(武汉)地质学专业,主要从事地球化学研究工作。

摘要

以沅麻盆地1:20万区域地球化学水系沉积物数据为基础,对沅麻盆地九曲湾铜矿床的地质、地球化学特征及区域异常信息的判别与提取进行分析,建立了九曲湾典型铜矿床地质—地球化学找矿模型,并成功圈定了预测区。由此可知该方法可行,应用效果明显。

关键词: 沅麻盆地; 地质—地球化学模型; 异常信息; 找矿标志; 预测区圈定
中图分类号:P632 文献标志码:A 文章编号:1000-8918(2016)02-0272-07 doi: 10.11720/wtyht.2016.2.07
The application of geological-geochemical modeling to Yuanma Basin, Hunan Province
YANG Bao-Guo
Ningxia Institute of Geophysical and Geochemical Exploration, Yinchuan 750004, China
Abstract

Based on the data of 1:200 000 regional geochemical stream sediments in Yuanma Basin, the authors analyzed geological and geochemical characteristics as well as the discrimination and extraction of regional anomaly information of the Jiuquwan copper ore deposit in Yuanma basin. Then the authors established a typical geological-geochemical prospecting model for the Jiuquwan copper ore deposit, and successfully delineated the prediction area. It can thus be seen that the method is feasible and effective.

Keyword: Yuanma basin; geological-geochemical model; anomaly information; prospecting criteria; delineation of prediction area

就目前国内外所建立的地球化学模型而言, 大多是基于岩石地球化学和矿区地质资料建立的针对单个矿床的地球化学模型, 缺乏利用区域水系沉积物资料建立的地球化学模型, 因此也难从“ 面” 上把握整个区域矿产资源的分布情况[1, 2, 3, 4, 5]。此次研究主要以成矿成晕地质、地球化学理论为指导, 结合沅麻盆地区域地质、地球化学资料, 应用元素地球化学特征参数、衬值异常、综合异常、相似度等地球化学指标, 对区域成矿条件进行分析, 并在前人对九曲湾典型铜矿床研究的基础上, 归纳出一套显示找矿标志和矿床特征的地质— 地球化学参数, 进而建立典型矿床的铜矿地质— 地球化学找矿模型[6, 7]

1 地质概况
1.1 区域地质概况

沅麻盆地出露地层有侏罗系、白垩系及古近系等, 区域内分布范围最广的是白垩系, 其次为侏罗系, 古近系较少[8]。侏罗系主要在区域内南东侧出露, 岩性以暗紫红色、黄绿色的砂岩、粉砂岩及泥质粉砂岩为主。

沅麻盆地是受两次地质构造运动作用形成的山间断陷盆地, 形成于印支、燕山期。整个盆地被看成一个向西歪斜的不对称向斜[9, 10]。区内构造主要有断裂和褶皱。断裂主要出现在盆地的北段与南段, 多为大型断层。褶皱主要有燕山期后形成的晓坪复式背斜和莱草坡背斜等, 盆地被这些背斜分割成多个沉积中心[4]

1.2 典型矿区地质概况

九曲湾铜矿含铜岩系属白垩系下统锦江组中上部的浅色层和新近系中新统罗旧组底部的砂砾岩层。矿区白垩系中神皇山组含铜岩系主要由浅色层和深色层即红层交互组成[11, 12], 矿体大多常赋存在浅色层中, 岩石是一套粒级不同的碎屑岩, 岩石类型较丰富, 从细砾岩到泥岩均有发现。

九曲湾矿区处于茶溪— 吕家坪不对称向斜之南翼[13, 14]。矿区断层比较发育, 常见南北向纵断层和东西向横断层。它们常以群产出并相互切割。

含铜岩石(矿石)中, 造岩矿物主要有钾长石、斜长石、石英、黑云母。矿体形态大多为扁豆状和透镜状, 少数以似层状或层状赋存在灰绿色、灰白的含矿层之中[7]。矿石中的主要有用矿物为自然铜, 其次为辉铜矿、孔雀石、赤铜矿、蓝铜矿等, 黄铜矿、斑铜矿偶尔可见, 但含量较少。各种铜矿物产出形式各异, 有的单独产出, 而有的则密切共生。九曲湾铜矿床最大的特点就是以自然铜为主的含铜砂岩型铜矿[11, 15]

2 区域地球化学特征

根据已有的地球化学调查资料, 对区域水系沉积物数据中一些不符合正态分布的数据用统计软件SPASS剔除掉, 最终得到服从正态分布且能反映初始本底的元素丰度值(背景值)。在化探中以找矿视角来分析, 将背景区定义为未受成矿作用影响的地区, 将地球化学异常区理解为与背景区相比有显著差异的元素含量富集区或者贫化区, 其中异常下限就是背景上限值, 据此绘制元素地球化学异常图。

Cu地球化学异常主要出现在沅麻盆地的西南部、中部, 在东北部几乎无异常(图1)。在区内已知铜矿床点上, Cu异常与矿点所在位置对应较好, 异常较强且异常面积较大, 多呈面状和带状分布, 走向多呈北东向, 尤其在九曲湾Cu异常显示尤为强烈, 内带面积很大。这些Cu地球化学异常显示较强的区域对应的地层主要为白垩系和古近系, 亦有一小部分侏罗系和第四系地层。

图1 沅麻盆地地质概况及Cu地球化学异常分布K— 白垩系; J— 侏罗系; P— 二叠系; D— 泥盆系; O— 奥陶系; Э — 寒武系; Z— 震旦系; Nh— 南华系; Qb— 青白口系

Ag地球化学异常在整个沅麻盆地分布较为广泛。区内九曲湾铜矿床、麻阳县黄双铜矿、泸溪砂岩铜矿、泸溪上堡铜矿Ag异常显示较强, 内、中、外带清晰, 异常面积较大。区内其他矿种如古丈县曹家坪黄铁矿、沅陵县王公坪方铅矿亦有Ag的强异常显示。另外Ag在区内没有矿点的地方也有强异常。结合地质图可以发现, Ag地球化学高值区对应的地层以白垩系和古近系为主, 在寒武系和震旦系同样可发现Ag异常, 但异常面积较小, 异常强度较低, 仅有外、中带。

Cd地球化学异常在整个沅麻盆地均有显示, 相比其他主成矿及伴生元素, 其异常面积较大, 异常强度较高。Cd在九曲湾铜矿床上的异常较弱, 仅有中带出现, 且异常面积不大, 但在区内其他铜矿床如泸溪上堡铜矿、泸溪砂岩铜矿, 异常面积大且内带较强, 呈带状分布, 走向呈北东向。Hg地球化学异常在整个盆地均有显示, 但强度均较小, 在区内已知矿床点的异常显示较区内其他地方的异常强, 如典型铜矿床九曲湾, Hg异常面积较大, 走向呈近东西向; 在泸溪砂岩铜矿异常较强, 走向北东。Pb地球化学异常在区域内均有显示, 但较强异常主要出现在整个区域的东北部、南部。在麻阳九曲湾铜矿床上, Pb异常面积较大, 但强度较小, 走向与Hg地球化学异常走向一致, 呈近东西向。

纵观主成矿及伴生元素在沅麻盆地的地球化学异常显示可以发现, 这些元素异常显示较好的区域, 其对应的地层主要为白垩系和古近系, 而这两种地层恰为研究区内砂岩型铜矿的赋矿层位。

3 典型矿床地球化学特征

通常来讲, 衬度异常量可被用来判别一个矿床或者矿田的主成矿元素及主要伴生元素, 也可用来预测研究区的矿种类型。

在九曲湾铜矿, Cu的异常面积最大, Ag次之(表1)。Cu衬度异常量最大, 达到1 279.68。总体来看, 该矿区主成矿元素中Cu、Ag的衬度异常量明显高于其他元素, 能较好地反映典型矿床的矿种或矿化类型。因此Cu、Ag为矿区主成矿元素, Pb、Cd、Hg为伴生元素。

表1 九曲湾铜矿主成矿元素和主要伴生元素异常参数统计

从图2中可以看出, Cu、Ag、Pb、Cd、Hg在九曲湾铜矿具有不同程度的异常显示, 其中主成矿元素Cu面积大, 异常内、中、外分带清晰; Ag、Cd相对较弱, 浓集中心与矿区位置却对应较好; Pb、Hg异常面积较大, 异常呈近东西向条带状, 浓集中心倾向矿区外围。

图2 九曲湾铜矿主要成矿元素和伴生元素地球化学异常分布示意

从地球化学图中可以看出, 主成矿元素异常显示好, 浓集中心与矿区对应好, 伴生元素相对较弱。由此可知, 在1:20万水系沉积物地球化学数据中, 可以找到很多与成矿关系密切的相关地球化学信息。

4 地球化学异常信息提取
4.1 弱小异常识别

为了能够发现研究区低缓、弱小的异常, 制作了元素衬值图。主要采用平均圆滑滤波法, 把地球化学背景当作一个连续变化的曲面, 以窗口来替代子区, 用大窗口作局部背景、小窗口作局部噪声来模拟地球化学背景的变化轨迹。以小窗口和大窗口之比值为衬值, 以变化着的衬值来替代不变的异常下限, 采用累频分级来圈定异常, 分级频数为25%、50%、80%、90%、95%、98%(表2)。

表2 九曲湾铜矿主成矿元素和主要伴生元素衬值参数统计

衬值异常图相比于地球化学异常图不仅突出了弱异常, 而且强化了强异常。图3是传统“ C+2S” 制作的Cd地球化学图件, 图4是利用衬值计算制作的Cd衬值地球化学图。通过对比可以看出:①Cd地球化学异常总体较弱, 仅在沅陵县王公坪方铅矿的北东向上出现较强的异常, 其他区域异常较弱; 在Cd地球化学衬值图上, 已知铜矿床上的Cd衬值异常明显加强; ②Cd地球化学图中泸溪上堡铜矿Cd无异常显示, 在Cd地球化学衬值图中, 泸溪上堡铜矿出现了Cd异常, 异常强度虽相对其他矿区较弱, 但与成矿密切相关; ③古丈县曹家坪黄铁矿一带的Cd地球异常显示较弱, 在Cd地球化学衬值图中, 异常的浓集中心更加凸显。

图3 沅麻盆地Cd地球化学异常分布

图4 沅麻盆地Cd地球化学衬值分布

4.2 综合异常

为能够更好地体现研究区的一些成矿信息, 选取与Cu成矿关系密切的Ag作为伴生指示元素, 制作研究区Cu-Ag的综合异常图(图5)。

图5 沅麻盆地 Cu-Ag综合异常分布

从沅麻盆地Cu-Ag综合异常图中可以看出, 区域内已知铜矿床如泸溪上堡铜矿、泸溪砂岩铜矿及典型矿床九曲湾铜矿上综合异常较强, 综合异常与这些已知铜矿床对应较好。另外在整个沅麻盆地其他未知区域亦有较强的综合异常显示。这就为进一步研究未知区域是否成矿提供了地球化学依据。

4.3 相似度判别

在找矿过程中常用已知来推测未知, 对已发现的矿床进行深入研究, 找出其成矿规律, 运用已知矿床已有的必要成矿因素来寻找相类似的矿床。为此需引入相似系数这一参数, 以此判断未知矿区的样本与典型矿区的标准样本之间的相似程度, 通过样品间的相似度判别来把握典型矿床与预测区的相似程度, 具体作法如下。

找出最能代表典型矿床矿化信息的元素组合, 即标准样本(虚拟样本)。通常选取典型矿床1:20万化探数据中异常内带或中带较为清晰的元素组合作为标准样本的元素组合, 并且这些元素与预测矿种要有相关性(一般通过相关性分析及经验判断)。在标准样本元素组合确定后, 就要对选取的元素组合进行取值, 通常每个典型矿床采用各元素异常内带(中带)内的平均含量进行赋值, 并将该值作为标准样本的含量特征值 (表3)。

表3 标准样本元素组合及标型元素含量值

标准样本与样品的相似度一般通过距离函数进行表示:

D(si)为未知区实际样本与标准样本的相似距离值, s为已知典型矿床标准样本, i为被判别的未知区(实际样本), p为元素组合的元素个数, Xsk为标准样本的第k个元素(如Ag)的取值, Xik为第i个实际样本在第k个元素上的取值。Ck为第k个元素与Cu的相关系数。DmaxD(si)的最大值, R为未知区实际样本与标准样本的地球化学相似系数或相近系数, 则

R=1-D(si)/Dmax(2)

为了克服极高值和元素量纲的影响, 需对原始数据进行均一化处理, 即除以该元素背景值。通过式(1)可知, D(si)值越大, 相似程度就越低, 反之亦然。为了更能说明问题, 把距离转换成相似系数R, 由式(2)可知, 所有样点与标准样本之相似度值的范围都在(0, 1)区间, 其中典型矿床自身与标准样本最相似, 即相似度值是1, 除此之外都小于1。由公式可知R越大, 则实际样品与标准样本的相似程度就越高。将每个样点的相似度值求得后, 采用7级累频分级的方法进行成图, 分级频数为50%、80%、90%、95%、98%、99%(表4)。

表4 相似度系数分级

异常显示的地区主要是与典型矿床成矿较为相似的地区, 这是圈定预测区的主要依据, 相似度异常显示越高, 表明与典型矿床的相似程度就越高。从图6中可以看出, 在整体上, 相似度异常在沅麻盆地中有零星的小块状异常分布。从矿床(矿点)角度看, 盆地内已知的绝大多数砂岩型铜矿床(矿点)都位于相似度异常的高值区内或边缘; 除铜矿以外的其他矿种所在的矿床位于其相似度异常的低值区。从研究区内已知不同矿床的相似度异常显示强弱上就很好地印证了相似度判别的准确性。这就说明相似度判别的效果很明显。

图6 沅麻盆地地球化学相似度分布

综上所述, 在有相似度异常但未见已知矿床或矿点的区域, 很有可能成矿, 应作为预测区圈定的主要区域, 在预测区圈定时应考虑相似度的影响。相似度地球化学图不但对已知矿床(矿点)有很好的示踪作用, 而且可以根据相似程度之高低进一步优选区域异常, 进而判定其矿化类型及矿化强弱。

5 九曲湾典型铜矿床地质— 地球化学找矿模型

地质— 地球化学找矿模型的建立需要对同类型典型矿床的地质特征和地球化学特征进行全面的研究与总结。九曲湾典型铜矿床的找矿标志概括如下。

5.1 地质找矿标志

地层标志:九曲湾铜矿床严格受层位控制, 矿体主要赋存于下白垩统神皇山组浅色含砾中— 粗粒砂岩中, 下白垩统— 古近系浅色砂岩是直接找矿标志。

构造标志:向斜核部、背斜倾伏端、断层与含矿层交切部位是含矿热液迁移富集的有利部位。

围岩蚀变标志:褪色化, 常见紫色层褪色为浅色层。

地表矿体标志:地表氧化铜矿体是直接找矿标志。

5.2 地球化学找矿标志

在判别矿床的主成矿元素及伴生元素时, 衬度异常量具有明显效果, 故常用衬度异常量来判定预测区之矿种类型。

成矿元素铜的1:20万水系沉积物地球化学异常区在铜矿种预测中是最直接的找矿预测区。

以九曲湾铜矿Cu、Ag、Pb、Hg、Cd为标型元素绘制其相似度图, 使与典型矿床相类似的一些铜矿化信息得以充分地反映。故可根据相似程度高低进行区域异常的优选, 从而判定其矿化类型及矿化强弱。

元素的衬值异常在识别与矿化相关的低缓、弱小异常中效果显著, Cu、Ag、Pb、Cd、Hg元素组合衬值异常在优选与评价预测区中起到较好的约束作用, 选择元素组合多且套合程度高的未知区作为预测区, 可提高预测区圈定的可信度。

选取沅麻盆地内的主成矿及主要伴生元素Cu、Ag制作Cu-Ag综合异常图, 能更好地凸显与砂岩型铜矿化相关的矿化信息。

6 找矿预测区圈定

在典型矿床地质、地球化学特征研究及地球化学找矿模型建立的基础上, 对预测靶区圈定的地质地球化学指标进行制定, 概括为:

条件1:九曲湾矿床的相似度图中相似度异常明显(累频分级≥ 98%);

条件2:成矿的地质条件较好(下白垩统第三岩组(K3)是主要的赋矿层位; 向斜核部、背斜倾伏端、断层与含矿层交切部位是含矿热液迁移富集的有利部位);

条件3:已发现铜矿点或铜矿化点;

条件4:Cu-Ag综合异常显示较好(累频分级≥ 95%);

条件5:Cu衬值不低于1.15(累频分级≥ 80%);

条件6:Cu、Ag、Pb、Cd、Hg中至少有3个元素衬值不低于其衬值异常下限(累频分级≥ 90%);

条件7:Cu、Ag、Pb、Cd、Hg中至少有4个元素衬值不低于其衬值异常下限(累频分级≥ 90%)。

根据有无已知矿点或者矿化点可将预测区初步分为A级、B级两级; 在B级中如果没有相似度, 则被划分成C级。分级条件如下:

A级:符合条件1、2、3、4、5、6; B级:符合条件1、2、4、5、6(缺少3); C 级:符合条件2、4、5、7(缺少1和3, 但必须符合7)。其可信度级别由高到低为:A级、B级、C级。

按照以上分级原则进行圈区, 具体方法为:先检索出符合A、B级和C级的点, 该步主要使用Excel的条件检索功能来完成, 将检索出的样点利用GIS技术转化为点文件, 再将这些点投在地质矿产图及其地球化学图(综合异常图、元素组合异常图、相似度图等)上作优选评价(面异常), 此过程是由点到面的综合评价过程。从而优选出不同级别的预测区(图7)。

图7 沅麻盆地预测区分布

7 结论

地质— 地球化学建模是在对典型矿床成矿模式研究的基础上, 对九曲湾典型矿床的关键地质因素和地球化学信息加以综合, 形成一个较为完整的成矿系统, 并由此来确定其找矿方向。

这种找矿模型就是利用已知来预测未知, 在未知区域若能满足典型矿床找矿模型所必须的条件, 就可说明该区域和典型矿床有相似的地质、地球化学特征, 成矿潜力大。

通过九曲湾典型铜矿床成矿模型, 在沅麻盆地成功圈定预测区11个, 其中A级预测区2个, B级预测区5个, C级预测区4个。圈定的两个A级找矿预测区恰好是两个已知的砂岩型铜矿床, 说明找矿模型的推测与实际存在能够很好地吻合, 同时也说明该地质— 地球化学建模方法可行, 应用效果明显。

The authors have declared that no competing interests exist.

参考文献
[1] 曾键年, 范永香, 潭铁龙. 现代成矿预测若干理论述评[J]. 矿产与地质, 1996, 10(6): 361-367. [本文引用:1]
[2] 赵鹏大, 孟宪国. 地质异常与矿产预测[J]. 地球科学: 中国地质大学学报, 1993, 18(1): 390-47 [本文引用:1]
[3] 刘大文. 地球化学块体理论与方法技术应用于矿产资源评价的研究[D]. 北京: 中国地质科学院, 2002. [本文引用:1]
[4] 阳正熙. 西方国家的“成矿规律和成矿预测”的发展和现状[J]. 成都理工学院学报, 2000, 27(增刊): 259-264. [本文引用:2]
[5] 李随民, 姚书振, 周宗桂. 矿产资源定量预测的研究现状[J]. 地质找矿论丛, 2007, 22(1): 9-18. [本文引用:1]
[6] 马振东, 曾键年, 龚敏, . 矿产资源地球化学模型建立与定量预测方法技术二[R]. 中国地质调查局发展研究中心, 2010. [本文引用:1]
[7] 杨向荣, 张晓帆, 孙宝生. 基于勘查地球化学数据分析的大型矿集区预测模型讨论[J]. 中国地质, 2006, 33(3): 682-688. [本文引用:2]
[8] 曾勇. 沅麻盆地砂岩型铜矿控矿条件[J]. 国土资源导刊, 2007, 4(3): 24-26. [本文引用:1]
[9] 郑贵洲, 罗新民, 侯光久, . 沅麻盆地白垩系划分及岩石地层研究[J]. JOURNAL OF STRATIGRAPHY, 1998, 22(3): 201-205. [本文引用:1]
[10] 贾宝华. 湖南雪峰隆起区构造变形研究[J]. 中国区域地质, 1994, 52(1): 65-71. [本文引用:1]
[11] 许首权. 九曲湾砂岩铜矿床成因及富集规律研究[J]. 矿产与地质, 1988, 2(4): 25-29. [本文引用:2]
[12] 余左清. 麻阳九曲湾砂岩铜矿床地质特征及富集规律与找矿方向[J]. 湖南地质, 2001, 20(1): 22-24. [本文引用:1]
[13] 尹承忠. 九曲湾砂岩铜矿床矿体的空间分布规律及找矿方向[J]. 湖南冶金, 2001, 1: 21-23. [本文引用:1]
[14] 黄满湘. 湖南麻阳铜矿成矿机制探讨[J]. 大地构造与成矿, 1999, 23(1): 42-49. [本文引用:1]
[15] 王大鹏, 张乾, 祝朝辉, . 湖南麻阳铜矿自然铜形态特征及其成因意义[J]. 矿物学报, 2007, 27(3): 384-387. [本文引用:1]
[16] 蒋金明, 沙元成, 文亮先, . 江西瑞昌市仙姑台铜多金属矿综合勘查模型[J]. 物探与化探, 2015, 39(6): 1104-1112. [本文引用:1]
[17] 闫红圃, 刘文斌, 程兴国. 津巴布韦Mutawatawa地区Chelsea East金矿床多元信息综合找矿方法及找矿模型[J]. 物探与化探, 2015, 39(3): 516-524. [本文引用:1]
[18] 任爱琴, 张宏伟, 吴宏伟. 河南千鹅冲钼矿地球化学异常特征及找矿模型[J]. 物探与化探, 2014, 38(5): 865-871. [本文引用:1]
[19] 温德娟, 滕寿仁, 刘长纯, . 辽宁西部肖家营子钼矿地球化学找矿模型[J]. 物探与化探, 2013, 37(6): 998-1002, 1007. [本文引用:1]