南秦岭夏家店卡林型金矿床构造特征与成矿关系
Tectonic and mineralization of Xiajiadian Carlin-type gold deposit in South Qinling,Shaanxi
-
摘要: 从地质背景、矿床地质、控矿构造分析了夏家店金矿床特征及其形成过程及成矿机理;矿床产于南秦岭陡岭古隆起边缘寒武系-震旦系碳-泥-硅质板岩及白云岩中;控矿构造具有二次韧性剪切变形、二次构造角砾岩化及二次脉体矿化的特征,显示出式样多种、性质多样、形式多变、活动期次与阶段多次叠加的特点;具体阐述了各构造作用阶段与矿化蚀变之间的关系,建立了地层建造-剪切与张裂构造-热液渗流的成矿系统.属于沉积岩系、构造、热液活动3种地质作用叠加改造的卡林型金矿床.Abstract: Based on study of geological setting,ore-controlling tectonic,the paper probes into the geological characteristics and metallogenic mechanism of Xiajiadian gold deposit.The deposit occurs in the carbargillite-siliceous rock and ferroan carbonate rock of the Sinian-Cambrian system on the border of the Douling doming paleostructure in South Qinling.The tectonic system is become ore-control leading factor,this is manifold type and shape,varied nature,pulsution of many time of the tectonic activity.By study of relationship of structural analysis and mineralization-alteration has been set up remould met allogenic system of strat a control-shear and tension struct ure-hydr othermal convection.Xiajiadian gold deposit of like Carlin type is subjected to control of manifold control and reform of geological condition.
-
近年来,随着计算机技术的不断发展,磁法勘探资料反演研究已由原始的手工计算发展到以计算机为主的算法研究,由单一反演方法发展到多方法的综合反演技术,由二维平面向三维立体发展。磁法三维反演解释技术是近年来国内外重磁研究的热点之一。磁性反演主要有两类,一类是将场源的边界形态作为反问题求解的形态反演;一类是将场源的物性参数作为反问题求解的物性反演。形态反演的模拟是通过任意多边体或多面体来模拟复杂地质体,它可以利用人的经验对反演进程加以引导,通过人机交互模式将反演结果反馈,根据反馈结果进行约束,进而提高反演效率和形态。磁法建模基于磁法形态反演的发展,经历了的人机交互建模技术、三维可视化技术等,对于真实地质体的模拟能力有了较大提升。
通过各种地球物理观测数据建立地表以下地质模型,结合地球物理异常数据及地质剖面进行建模是了解地下深处构造单元的重要方法,Goleby等(2002)在澳大利亚Yilgarn东Norseman-Wiluna地区采用2.5D重力模拟,研究了该区绿岩的深度及与花岗岩之间的位置关系。Roy等(2000)利用反射地震、重磁2.5D和3D建模,对加拿大哥伦比亚省富含斑岩型铜钼矿的Guichon Cheek岩基进行了深部结构研究,建立了三维模型,岩基的边界、内部结构、岩浆通道清晰可见。Malehmir等(2006,2007)在瑞典北部Skellefte成矿带进行了的3D地质-地球物理建模,瑞典Upplasa大学在该成矿带西部的Kristineberg矿集区开展了高分辨率反射地震剖面和重、磁、震联合反演解释研究,建立了该地区3D地质-地球物理模型。向中林(2009)将三维地质建模及可视化用于矿床的成矿地质条件分析,应用于危机矿山找矿。吕庆田等(2010)在铜陵矿集区狮子山矿田利用反射地震成像及钻孔资料约束,进行重力3 D反演,建立了3D地质模型,初步实现了狮子山-铜官山矿田地壳结构“透明化”。祁光等(2012)以安徽泥河铁矿为例开展先验地质信息约束的三维地质重磁建模研究;严加永等(2014)在安徽省沙溪铜矿进行三维重磁反演,通过反演磁化率体和密度体,识别了4种主要的岩体类型,为三维地质建模提供借鉴。陈炳锦(2021)在陕西省龙王沟地区以磁测数据为基础建立了三维地质体模型,直观预测磁铁矿的空间分布和形态。王昊(2023)基于航磁数据开展位场分离获取重磁反演的异常数据,完成人机交互反演,通过剖面建立了朱溪矿矿区三维地质建模。随着信息化的发展,三维地质建模及可视化已呈现出较为明显的优势,也是对于开展深部找矿研究的重要手段。
1984年,法国地质物理学家Morlet首次提出“小波分析”的概念,引入了被命名为“Morlet小波”的一种基函数,用于时间和频率的局域变化,进而在信号中提取信息,通过伸缩和平移等功能实现对信号的多尺度细化分析,解决了许多难题。1987年,信号分析专家Mallat将计算机的多尺度分析思想引入到小波分析中,提出了多分辨分析的概念,研究了小波变换的离散化情形,提出了著名的分解与重构快速算法,即“Mallat”算法。Wickerhuaes(1990)、Beylkin(2010)等提出了“小波包”的概念,将Mallat进一步优化,提出了小波包算法,经过30多年的发展,小波理论不断完善,成熟,小波分析在信号处理、图像处理、语音识别等众多非线性学科领域取得了重大成就。小波多尺度分析法可以将磁异常分解到不同的尺度空间中,并且尺度的大小决定了异常所反映的地质体的规模和埋深情况。侯遵泽(1995,1997)、杨文采(2001,2004)等将小波变化应用于重力异常分析,李宗杰(1997a,1997b)在位场数据处理中引入小波变换,将异常信号分解、滤波、重建,达到了提取区域异常和局部异常的目的。张恒磊(2009)研究了小波分析的磁测数据处理流程,
宁津生(2010)等利用连续小波变换进行了干扰源分离和场源深度确定的模拟,发现小波变换边缘分析可有效确定场源深度。刘芳(2013)等通过二维小波多尺度分析,实现了不同场源深度密度体与重力异常分离的目的。刘天佑(2007)、张恒磊(2009)、王建复(2013)、尚世贵(2014)、宋小超(2016)、蒋勇平(2018)等将小波多尺度分析用于重磁法数据处理,精细处理携带了多尺度地质信息的位场信号,使得小波多尺度分析方法在资料的位场分离领域得到了广泛的应用。
根据频率与深度的对应关系可知:小尺度代表高频信号,大尺度代表低频信号,低频部分更多反映的是深部信息,高频部分反应浅部信息,通过小波多尺度因子变化将异常信号分解到不同的细节。借助功率谱计算,来反映不同细节对应的异常源深度。功率谱可以借助重磁异常的径向对数功率谱分析,来定量的确定重磁异常的场源深度。Bhattacharyya(1966)给出了频率域内单个均匀磁化棱柱体的总磁场场强表达式,提出了一种矩谱法来直接计算磁性体界面深度;Spector等(1970)提出了一种基于统计理论模型的“等效论”磁异常数据处理方法。侯重初(1985,1985a)等利用磁性体导数异常的对数功率谱直接计算磁性体下底深度推导出利用磁异常垂向一阶导数的对数功率谱计算二度板状体下底埋深的近似公式。申宁华(1985)提出先计算磁性体的顶深和质心深度,在“磁性体顶部距离质心与质心距离底部相等”的假设条件下,间接求出磁性体的底部埋深。张先(2007)等研究了不同尺度磁性体的场源深度。
党月辉(2008)研究了功率谱估算磁源深度的方法,给出了使用条件和范围。段瑞锋(2016)将功率谱计算用于银额盆地居延海坳陷磁测异常解释中; 笔者利用小波多尺度分解将陕西省略阳县金子山地区磁测数据分解为几个不同阶的细节,通过引用直立矩形棱柱体功率谱计算各阶细节异常对应的场源深度。结合区域地质背景,地形地貌及钻孔信息,利用Voxler平台,在虚拟环境下实现二维小波多尺度分解数据的三维离散可视化,建立三维地质体模型,预测岩浆热液通道部位、岩体空间形态,并结合时间域激电测深,圈定矿体的赋存空间,提高地球物理勘探解释的准确性、可靠性,实现地质、物探和地表及深部信息的集成展示,服务于分析研究和制定决策,为下一步钻孔设计提供依据。
1. 研究区区域地质背景
陕西勉(县)–略(阳)–阳(平关)三角区是重要的多(贵)金属矿集区,常见有中酸性沉积岩带和超基性–中基性–中酸性侵入岩及古生代碎屑岩和碳酸盐岩。整体上呈“东收西开”的扇骨形区域构造。南侧为汉江大断裂控制,北面为二里坝—铜厂–七里沟断裂,控制着区域内构造活动和岩浆活动。金子山金铜多金属矿即分布在二里坝—铜厂–七里沟断裂中段(图1)。
图 1 勉略阳地区区域地质构造略图(据西北有色地质勘查局七一一总队,2015)1.勉略康构造混杂岩带;2.震旦纪碎屑碳酸盐岩;3.太古界绿岩;4.中部火山岩浆岩带;5.中晚元古界中酸性火山岩;6.基性岩体;7.中下元古界基性火山岩;8.基底拼合主构造线;9.超基性岩;10.闪长岩;11.太古界鱼洞子岩群;12.中下元古界东沟坝组;13.中下元古界何家岩岩群;14.金矿床;15.铜矿床;16.多金属矿床;17.镍矿床;18.铁矿床;19.古基底缝合带;20.地名;21.研究区Figure 1. The Regional geological structure of Mianlveyang area2. 研究区地质地球物理特征
工作区出露地层主要为郭家沟组第一亚层、震旦系断头崖组及第四系(图2)。郭家沟组第一亚层主要岩性为细碧岩。震旦系断头崖组主体分布在测区西北部,与下伏火山岩呈不整合接触。是一套以灰岩、白云质灰岩、含碳绢云母板岩、碳质板岩、板岩、凝灰质板岩构成的碎屑化学沉积变质岩。区内断裂构造发育,分为近EW向、NW向、近SN向和NE向4组。控矿断裂为NE向断裂,其余为破矿构造。区内侵入岩发育,超基性岩大面积分布,基性、中酸性岩体多呈岩墙、岩脉或岩株产出。超基性岩第一期主要为菱镁岩、滑镁岩,磁性较弱,据研究成果表明MgO含量明显偏高,CaO和Al2O3偏低,属深源浅成岩浆类型,据纯橄岩年龄测试Rb-Sr等时线年龄(927±49)Ma,属晋宁期产物;第二期主要为蛇纹岩、碎裂状蛇纹岩、纯橄榄蛇纹岩(原岩:橄榄石)、纤胶蛇纹岩、叶蛇纹岩、绢云母化蛇纹岩、斜辉绿橄岩、含磁铁矿蛇纹岩,磁性较强,与晋宁期超基性岩带同位产出,构成复合超基性岩带,同位素年龄328~540 Ma(K-Ar法),为加里东期–华力西期产物。根据定向标本的测定结果显示,二里坝–铜厂背斜发生于该期次岩体产出之后,最后在金子山地区出露剩磁方向较为凌乱的含磁铁矿角砾状蛇纹岩,为负异常,金子山的金铜矿出露于含磁铁矿角砾状蛇纹岩中。对前期施工的探槽钻探工程,结合民坑调研工作重新进行综合研究后,认为金子山金铜矿体是加里东期–华力西期产物,后被含磁铁矿角砾状蛇纹岩推覆至地表。
根据测定的物性标本显示(表1),磁性最强的是磁铁矿和蛇纹岩,郭家沟组的细碧岩、震旦系的碎屑化学沉积变质岩、晋宁期的菱镁岩和滑镁岩、云英岩、斜长花岗岩磁性较弱。在本工作区蛇纹岩呈低阻低极化,滑镁岩呈高阻低极化特点,但是含磁铁角砾状蛇纹岩中含有大量的硫铁矿,呈低阻高极化特点。综上所述利用蛇纹岩与围岩的磁性、电性差异,预测浆浆热液通道部位、岩体空间形态是可行的。
表 1 工区岩(矿)石标本磁参数测定统计表Table 1. Area of rock (ore) specimen magnetic parameters measurement calculation table岩(矿)石名称 κ/(4π×10−6 SI) Mr /(10−3 A/m) 变化范围 常见值 变化范围 常见值 磁铁矿 34 18000 ~63047 33480 2 000~ 20000 12850 含磁铁角砾状蛇纹岩 40 14141 ~82764 40728 1863 ~27505 7546 蛇纹岩 54 1000 ~3600 1900 400~ 2600 1100 滑镁岩 39 100~800 433 150~600 226 细碧岩 30 80~300 175 78~425 179 斜长花岗岩 31 60~400 141 60~300 133 云英岩 30 100~600 244 100~500 268 金子山金铜矿前期勘查工作对成矿控制因素、矿化类型研究较少,致使勘查工程见矿不佳,对整个金子山地区的主攻矿化类型、主攻矿种缺少全面认识,进而在勘查靶区、靶位选取及工程布设方面存在一定的盲目性,造成进一步勘查工作难以深入,找矿无法取得突破。笔者主要对前期施工的槽探钻探工程结合民坑调研工作进行重新综合研究后,认为金子山金铜矿为与加里东–华力西期超基性岩有关的岩浆熔离型金铜镍多金属块状硫化物型,利用蛇纹岩与围岩的磁性、电性差异,确定蛇纹岩产生的异常与其对应三维空间的定量关系,把蛇纹岩产生的异常范围划分出来,从而约束反演结果,得到最优解。通过预测岩浆通道底部或者转折端位置,以期发现深部的盲矿体。
3. Voxler建模
根据小波多尺度分解塔式算法的低阶细节不变性,对略阳县金子山地区垂直磁异常ΔZ(单位:nT)进行小波多尺度分解处理,根据小波多尺度因子的变化,将磁异常信号分解到不同阶的细节,不同阶的细节代表了不同深度的局部场。
利用Voxler三维可视化科学制图软件,对离散的垂直磁异常数据进行小波多尺度分解,分解为1~6阶细节和逼近(图3)。
据党月辉(2008)通过对Bhattacharrya提出的功率谱表达式进行化极处理和求取平均径向功率谱后,推导出的有限延深直立矩形棱柱体平均径向功率谱表达式
$$ \mathrm{I}\mathrm{n}\mathrm{ }\mathrm{E}\left(\gamma \right)=A-2{h}_{t}\gamma +2In(1-{e}^{-tr}) $$ (1) 对公式中的r求极值,并取自然对数,得到:
$$ r=\frac{In\left(1+\dfrac{h}{t}\right)}{t} $$ (2) 其中,A为常数;r为径向频率;t为棱柱体的延伸;h为埋深。
由上式可知r存在极大值拐点。当已知棱柱体埋深时,可以根据径向频率的峰值位置估计其延深t,计算出$ 2In(1-{e}^{-tr}) $的结果,反过来可以对估算埋深进行延伸影响的改正。
根据功率谱分析得出:1阶细节异常场源似深度为22 m;2阶细节异常场源似深度为61 m;3阶细节异常场源似深度为124 m;4阶细节异常场源似深度为238 m;5阶细节异常场源似深度为474 m;6阶细节异常场源似深度为860 m;6阶逼近场源似深度为1 812 m。异常主要为大面积加里东期–华力西期磁性较强的蛇纹岩引起,磁性体较为简单,通过地形改正,用似深度代替真深度,将传统的二维小波多尺度分解数据整合为三维数据(图4)。
借助VolRender(形体渲染)模块编制形体渲染图(图5),根据形体渲染图件色彩变化特征、区内地质及钻孔等先验信息,提取大于50 nT的磁异常区间,在Isosurface(等值面)模块下编制空间等值面图(图6)。
该等值面基本反映了金子山测区蛇纹岩的分布特征,代表了岩浆在金子山地区的侵入特征。从图中可以看出,金子山地区的超基性岩体仅为岩浆主通道北侧的一部分,根部位于金子山地区的西南方向,岩浆主通道的北缘向北陡倾,按照同位产出的规律,金子山地区的含磁铁矿角砾状蛇纹岩为岩浆主通道的分支,推断向南陡倾,结果和时间域激电测深推断的产状一致(图7)。前期施工4个钻孔均未见矿(图8),ZK2501和ZK2502孔为交叉孔,钻孔岩性均为滑镁岩,未见矿原因分析认为矿体控制深度不够,均未穿透滑镁岩,达到含矿岩层。ZKA01和ZK2701孔未见矿,主要原因是未考虑控矿岩层的产状问题,从含矿岩层的背后和侧面穿过。
通过对金子山金铜矿资料重新进行综合研究,结合岩浆通道的三维模型,推断在岩浆通道底部或者转折端位置存在与超基性岩有关的岩浆熔离型金铜矿体,施工了ZK2601孔(图8),在岩浆分支通道的底端,标高约
1270 m处见金铜矿体,Au品位5.2×10−6,Cu品位0.11×10−2,厚1.0 m。4. 结论
(1)笔者通过对陕西省略阳县金子山地区磁异常数据进行小波多尺度分解,通过功率谱与磁异常深度对应关系计算出磁源似深度,依靠钻孔数据进行约束反演,借助Voxler平台建立三维磁性体模型,直观反映了不同异常源对应的地质体和构造信息,提高了磁异常的垂向分辨率,有助于提高对深部异常及矿体赋存位置的认识,提高深部找矿的工作效率。
(2)Voxler平台为磁法三维反演提供了良好的平台,其多个数据模块在实现交互式可视化三维地质建模重构中起到了较好的作用,通过平台的可视化重构,以直观预测了岩浆的通道位置,推断了岩浆通道底部或者转折端位置,大致圈定矿体的赋存空间,提高了地球物理勘探解释的准确性、可靠性。
(3)在功率谱计算时引用了有限延深直立矩形棱柱体功率谱公式,仅考虑直立矩形条件下,且将少量钻孔资料信息作为地质建模约束条件,条件过于单一,针对复杂模型还需进一步研究。
-
张复新,季军良,龙灵利,等.南秦岭卡林型-似卡林型金矿床综合地质地球化学特征[J].地质论评,2001,47(5):492-499.
陈柏林.金矿床和金成矿作用研究进展[J].地质论评,2001,47(1):111-112.
张复新,宗静婷,马建秦.秦岭卡林型金矿床及相关问题探讨[J].矿床地质,1998,17(2):172-184.
张复新,魏宽义,马建秦,等.南秦岭微细粒浸染型金矿床地质与找矿[M].西安:西北大学出版社,1997.
张复新,马建秦,魏宽义.镇安米粮地区微细粒浸染型金矿床与构造演化的关系[J].矿床地质,1995,14(1):26-34.
任涛.陕西山阳夏家店金矿床地质特征及前景分析[J].西北金属矿产地质,合刊(26-27).2001.
张复新.秦岭微细粒浸染型层控锑金矿床金赋存状态研究[J].地质论评,1996,42(6):541-549.
郑明华,等.喷流型与浊积型金矿床[M].成都:四川科学技术出版社,1994.
卢纪英,李作华,张复新,等.秦岭板块金矿床[M].陕西科学技术出版社,2001.
张本仁,等.秦巴岩石圈构造及成矿规律地球化学研究[M].武汉:中国地质大学出版社.1994.
张复新,马建秦,陈衍景.秦岭卡林型金矿床金、砷地球化学探讨[J].地球化学,1999,28(5):453-463.
计量
- 文章访问数: 2461
- HTML全文浏览量: 0
- PDF下载量: 3375