本文基于金川铜镍矿床II-1号矿体的三维地质模型,结合地质资料以及勘探数据,利用三维空间分析方法对II-1号矿体的控矿因素进行定量提取,提取了岩浆通道骨架中心线距离场,断层距离场以及顶底板相对距离三个控矿因素。分析发现II-1号矿体的铜镍品位与提取的三个控矿因素存在较强相关性,因此可对控矿因素进行线性变换或非线性变换转变成找矿信息指标。此提取方法可应用于整个金川铜镍矿床,对金川铜镍矿床隐伏矿体三维可视化预测有重要意义。 Based on the three-dimension geological model of Jinchuan Cu-Ni sulfide deposit II-1 orebody, com-bined with geology materials and exploration data, this study uses three-dimensional spatial analysis method quantitatively extracted three ore-controlling factors of Jinchuan Cu-Ni sulfide deposit II-1 orebody. Relative distance between rock’s roof and floor and Euclidean distance field of magma channel central line and fault are generated. Analysis shows that there is a certain correlation between Cu and Ni grade of orebodies and ore-controlling factors extracted, thus linear transformation or non-linear transformation can be applied to these ore-controlling factors to construct ore-pros- pecting factors. It has a significant value to the further three-dimensional prediction of hidden orebodies in Jinchuan Cu-Ni sulfide deposit that the method is appropriate for the whole deposit.
朱佳玮1,2,毛先成1,2*,赵莹1,2,邓浩1,2,刘占坤1,2,张苗苗1,2
1中南大学有色金属成矿预测教育部重点实验室,湖南 长沙
2中南大学地球科学与信息物理学院,湖南 长沙
收稿日期:2015年12月10日;录用日期:2015年12月28日;发布日期:2015年12月31日
本文基于金川铜镍矿床II-1号矿体的三维地质模型,结合地质资料以及勘探数据,利用三维空间分析方法对II-1号矿体的控矿因素进行定量提取,提取了岩浆通道骨架中心线距离场,断层距离场以及顶底板相对距离三个控矿因素。分析发现II-1号矿体的铜镍品位与提取的三个控矿因素存在较强相关性,因此可对控矿因素进行线性变换或非线性变换转变成找矿信息指标。此提取方法可应用于整个金川铜镍矿床,对金川铜镍矿床隐伏矿体三维可视化预测有重要意义。
关键词 :金川铜镍矿床,三维地质模型,断层距离场,岩浆通道距离场,顶底板相对距离
金川铜镍硫化物矿床是世界上目前正在开采的第三大镍矿,仅次于俄罗斯的Noril’sk-Talnakh矿床以及加拿大的Sudbury矿床。自1958年发现至今,金川铜镍矿床出矿颇多,并且仍有找矿前景,而且近年来,也未曾发现新的大型硫化镍矿床,因此对金川铜镍矿床找矿预测意义重大。汤中立等针对金川铜镍硫化物矿床成因特点和岩体富集规律提出了“小岩体成大矿”以及“深部熔离–多期贯入”的成矿模式[
金川铜镍硫化物矿床位于华北板块西南边部的阿拉善地块西南缘的龙首山隆起带内。金川镁铁-超镁铁质侵入体大致以10˚交角不整合侵位于前长城系白家嘴子组中,其呈NWW向岩墙状直接侵位于片麻岩、大理岩和条带状混合岩中 [
金川II矿区位于断层F16-1与F23之间,是金川铜镍矿床资源量最大的矿区。II矿区岩体长约3000 m,除东部隐伏于混合岩下,其余地段均出露地表。由于张扭性断裂F17的作用,使岩体发生了约200 m的水平位错。岩体自西向东逐渐变宽,最宽可达530 m。岩体整体走向NW50˚,倾向SW,倾角约50˚~80˚,西陡东缓。从横剖面上看,西段岩体呈“凹”字形漏斗状,东段岩体呈“V”字形漏斗状;纵剖面上看,西段岩体向下存在两个分支,主矿体位于下盘岩枝,东段岩体呈歪漏斗状,矿体靠近底板。岩体由多期岩浆前后侵入形成,存在着明显的岩相分带。含矿部位岩体,以矿体为中心往外二辉橄榄岩(Σ3)、橄榄辉石岩(Σ4)。
II-1号矿体为金川矿床最大的矿体,位于II矿区西北侧岩体的下部,全长约1900 m,平均厚度约70 m。矿体呈透镜状、板状赋存于超基性岩体中,倾向SW,倾角25˚~75˚。围岩主要为二辉橄榄岩(占60%),次为大理岩(占15%),再次为橄榄辉石岩及斜长角闪岩。矿体以超基性岩型富矿(Ni品位1%~3%)为主,占矿体的87.49%,富矿体位于超基性岩体中下部,贫矿(Ni品位0.3%~1%)分布于富矿周围,局部也见有富矿体直接与围岩接触。
目前,岩浆通道系统成矿已经被越来越多的学者认可,并用来解释金川铜镍硫化物矿床的成矿模式[
综上,本文将对岩浆通道距离场、断层距离场以及岩墙顶底板相对距离这三个控矿因素进行定量提取并分析与矿化指标之间的规律。
岩浆通道可分为水平通道和岩墙型通道[
图1. 金川铜镍硫化物矿床地质简图(据汤中立等[
铁质岩体与围岩多具穿、切层现象。因此,形成II-1号矿体的岩浆通道属岩墙型通道。岩浆通道成矿理论认为,成矿的岩浆房是一个开放系统,既接受新岩浆的注入也接受围岩、地层的混染,当岩浆达到硫过饱和时,发生硫化物的熔离作用,硫化物乳珠就地或沿通道随岩浆移动到适当地方沉降形成矿体 [
金川岩体的形成至少经历了4期岩浆的侵位,II-1号矿体与第三期岩浆侵入关系最为密切 [
在综合地球化学、地质因素以及地质体模型后使用主成分分析法,并结合地质认识构建了岩浆通道的骨架模型。岩浆通道骨架模型是一中空的实体模型,为开展对岩浆通道距离场的提取,首先需要提取岩浆通道骨架模型的中心线并对其体素化(图2),然后计算矿体单元距离中心线的最小欧式距离值。
本文将点到地质体的距离约定为点到地质体的最小距离,因此矿体距岩浆通道的距离即转换为矿体单元中心点距其最近的岩浆通道中心线单元的距离。在基于体素模型的空间分析中需要大量的距离量算,若使用传统方法遍历空间中所有矿体单元,计算其与所有岩浆通道中心线单元距离并求取距离最小的单元将带来巨大的计算量。因此需要寻求一种能够快速在栅格空间生成距离场的量算方法,相对传统方法可高达O(n6)的时间复杂度,蔺宏伟提出了三维带符号的欧氏距离变换方法(3 Dimensional Signed Euclidean Distance Transform, 3-SEDT)将全局搜索问题转化成一个局部搜索问题从而有效减少计算量 [
本文应用3-SEDT算法过程如下:
1) 计算岩浆通道中心线以及矿体实体的包围盒B,本研究将包围盒扩大至整个地质空间。
2) 将B离散为体素,将岩浆通道中心线所在体素设为特征体素,记特征体素组成的集合为T1,并令其距离向量初值
3) 用3-SEDT中8个模板对B进行扫描,计算出岩浆通道中心线的距离场并对距离向量值更新,对矿体体素与B进行求交,提取出B中矿体所在体素,对距离向量求模,最终提取出矿体距岩浆通道距离场。
断层既是容矿构造同时也是导矿构造,原始为矿浆提供通道空间的断层虽然受到了后期运动的改造,
图2. 岩浆通道中心线体素化过程
但是依然对岩浆通道位置具有一定指示作用。金川矿床形成后主要处于挤压环境中,在最大剪应力的方向上,发育了以F8、F16-1和F23断层为代表的北东向左行平移断层 [
控矿地质因素场与空间中某点到相关联的断层面的距离有关,即控矿地质因素场是到断层面距离的空间分布函数 [
记断层所在体素为特征体素,其集合为T1,距离向量初值
计算矿体距离向量的模即可得出矿体离断层的距离值。值得注意的是,断层一般由上盘和下盘之分,断层的上下盘的物理化学性质一般不同;这里为了体现这种不同,约定位于断层面上盘的单元为正距离,位于断层面下盘的单元为负距离(图3)。
为了判断正负号,可以令a为从断层面中心点出发到矿体单元的向量,b为断层面的法向量,则
其中
II-1号矿体所处岩体为一岩墙,岩墙形成初期比现在产状平缓 [
图3. 断层距离场正负示意图
的距离并不能完全体现矿体在岩体中的分布位置,分布在厚大岩体底部的矿体单元距底板的距离可能大于分布在厚度较小岩体顶部的矿体单元距底板的距离,所以需要一个不受岩体厚度影响且能体现矿体单元在岩体中分布位置的因子,因此对每一矿体体元分别求取与岩体顶板和底板在厚度方向上的距离,并求取距离比值。因此,本文选取矿体在岩墙顶底板中的相对距离作为一个因素,从大量的地质数据中探寻II-1号矿体的富集规律。
求矿化单元到顶(底)板的距离的实质即寻找该单元质心点到顶(底)板面某一三角形的距离。为求取厚度方向上的距离,首先定义一个三维平面为走向面,其走向、倾向、倾角均与岩体一致,岩体厚度方向即该走向面的法线方向(图4)。
矿体单元与顶底板相对距离计算方法如下:
1) 首先对每个属于矿体的体元
2) 以走向面法线方向
将(2)代入规范化后的(3)式,可以得到:
图4. 岩体顶底板与走向面示意图
其中
若满足
3) 然后以走向面法线反方向
4) 求矿体单元与顶底板距离(图6)
5) 获取距离比值
其中R值越大,矿体单元距离底板越近;R值为0时,矿体单元与顶板接触;R值为1时,矿体单元与底板接触。
对II-1号矿体求取岩浆通道骨架距离场、断层距离场以及顶底板相对距离结果后对矿体数据库进行更新,新建字段dis_MC、dis_Fault、dis_Ratio分别存储三个控矿因素值。将II-1号矿体数据导入Gocad后分别以dis_MC、dis_Fault、dis_Ratio值进行颜色渲染,如图7~图9。
绘制矿体Cu、Ni品位与三个控矿因素的直方图如图10~图15。
由图10、图11可知dMC值范围为(0,700),在0附近时Cu、Ni品位最高,Cu、Ni品位随着dMC值增大而减小,衰减规律非常明显,可以很有力的支撑岩浆通道成矿。在硫化物上侵过程中,早期岩浆形成的矿体处于岩浆通道边缘部位,围岩混染最严重、温度变化最快。而且早期矿浆硫化物浓度相对较
图5. 点在三角形内
图6. 顶底板相对距离示意图
图7. II-1号矿体岩浆通道中心线距离场
图8. F16-1断层距离场(II-1号矿体)
图9. II-1号矿体顶底板相对距离
图10. Cu平均品位直方图(dMC)
图11. Ni平均品位直方图(dMC)
图12. Cu平均品位直方图(dFault)
图13. Ni平均品位直方图(dFault)
图14. Cu平均品位直方图(dRatio)
图15. Ni平均品位直方图(dRatio)
低,因而形成的矿体品位也相对较低(Cu < 0.5%,Ni < 0.8%)。随后,矿浆不断沿通道中心上侵。由于通道中心,尤其是靠近岩浆通道口附近,温度下降较慢,拥有较好的成矿环境。因此,靠近通道中心的位置,矿体品位最高。但是在dMC等于200和380值附近时,出现了两次峰值。根据已有认识,图中出现多峰最主要的原因为,矿体为多期岩浆贯入形成的复式矿体。在后期岩浆上侵过程,早期形成的岩体(可能尚未冷凝结晶)被不断往外挤压使其偏离岩浆通道中心,因此在距岩浆通道一定距离时出现峰值。另外,矿体富集也可能受多因素共同控制,并非仅与岩浆通道有关。
dFault值范围为0-900,由于II-1号矿体在F16-1上盘,故其断层距离场皆为正值。Cu、Ni品位与断层距离场指标非线性相关,Cu、Ni品位峰值偏左,在dFault值为300左右时Cu、Ni品位最高,其位置与岩浆通道中心相重叠,因此此处矿体富集可能为岩浆通道影响,但在(300,900)范围时,矿体单元距断层越远,Cu、Ni品位越低,在900时达到最低。
dRatio值范围为0~1,其值越接近1表明矿体单元离岩体底板越近。由图可知,近底板矿体单元Cu、Ni品位总体大于近顶板矿体单元,在dRatio值为0.6左右时,Cu、Ni品位最高。由于岩浆通道中心成矿环境最好,硫化物最为富集,结合图14、图15可知,岩浆通道中心在岩墙中部偏下位置(dRatio = 0.6)。矿浆在通道中侵位和冷凝过程中,一直发生着硫化物固溶体分离作用。硫化物乳珠受重力影响将向下沉淀,因此,矿体富集部位应偏于底板。但由于主成矿期中存在多期岩浆上侵,而富矿岩浆为晚期沿通道中心侵位,所以,矿体富集部位更偏于岩浆通道中心。两种因素的叠加,导致了矿体在岩墙中下部(dRatio = 0.6)最为富集。
1) 本文在三维地质模型的基础上,使用三维空间分析方法求取了金川铜镍矿床II-1号矿体的三个控矿因素,分别为顶底板相对距离、岩浆通道中心线欧式距离场以及其距F16-1的断层欧式距离场,利用Datamine、Gocad等三维建模软件对数据进行了处理及显示。对矿化指标以及控矿因素进行分析后,发现岩浆通道距离场与矿化指标相关性显著,离岩浆通道越近,Cu、Ni品位越高;顶底板相对距离与Cu、Ni品位之间的规律也能指示矿体富集于岩体中下部;断层距离场与矿化指标之间规律受岩浆通道影响但仍具较强相关性。可对三个控矿因素和矿化指标进行回归分析找出数学规律,并用于预测,对未来II-1号矿体深边部找矿有指示性作用。
2) 鉴于金川铜镍矿床的成矿模式,可以将同样的方法应用于III-1号矿体,I-24号矿体以及II-2号矿体定量提取三个控矿因素:构建III-58号矿体,I-24号矿体和II-2号矿体的岩浆通道骨架模型并求取岩浆通道距离场;分别求取对应的F8、F16-1以及F17的断层距离场;根据III、I、II矿区超基性岩体形态求取顶底板相对距离。找出控矿因素与矿体富集之间存在规律并验证其有效性后对控矿因素进行转换变成找矿信息指标,对整个金川铜镍矿床深边部找矿提供支持。
朱佳玮,毛先成,赵莹,邓浩,刘占坤,张苗苗. 基于三维地质模型的控矿因素提取—以金川铜镍矿床II-1号矿体为例 Extraction of Ore-Controlling Factors Based on Three-Dimensional Geological Model—A Case Study of Jinchuan Cu-Ni Deposit II-1 Orebody[J]. 地球科学前沿, 2015, 05(06): 435-448. http://dx.doi.org/10.12677/AG.2015.56044